----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 21:02:19]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [1 1 1] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fd-3m (227) Number of symmetry operations in supercell: 192 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 4.175817520000000 4.175817520000000 b 4.175817520000000 0.000000000000000 4.175817520000000 c 4.175817520000000 4.175817520000000 0.000000000000000 Atomic positions (fractional): *1 O 0.36360202051417 0.36360202051417 0.90919393845749 15.999 2 O 0.34080606154251 0.88639797948583 0.88639797948583 15.999 3 O 0.88639797948583 0.34080606154251 0.88639797948583 15.999 4 O 0.88639797948583 0.88639797948583 0.88639797948583 15.999 5 O 0.36360202051417 0.90919393845749 0.36360202051417 15.999 6 O 0.36360202051417 0.36360202051417 0.36360202051417 15.999 7 O 0.88639797948583 0.88639797948583 0.34080606154251 15.999 8 O 0.90919393845749 0.36360202051417 0.36360202051417 15.999 *9 Zn 0.75000000000000 0.75000000000000 0.75000000000000 65.380 10 Zn 0.50000000000000 0.50000000000000 0.50000000000000 65.380 *11 Ga 0.62500000000000 0.12500000000000 0.12500000000000 69.723 12 Ga 0.12500000000000 0.12500000000000 0.62500000000000 69.723 13 Ga 0.12500000000000 0.12500000000000 0.12500000000000 69.723 14 Ga 0.12500000000000 0.62500000000000 0.12500000000000 69.723 -------------------------------- unit cell --------------------------------- Lattice vectors: a 8.351635040000000 0.000000000000000 0.000000000000000 b 0.000000000000000 8.351635040000000 0.000000000000000 c 0.000000000000000 0.000000000000000 8.351635040000000 Atomic positions (fractional): *1 O 0.63639797948583 0.13639797948583 0.86360202051417 15.999 > 1 2 O 0.38639797948583 0.11360202051417 0.61360202051417 15.999 > 2 3 O 0.11360202051417 0.38639797948583 0.61360202051417 15.999 > 3 4 O 0.38639797948583 0.38639797948583 0.88639797948583 15.999 > 4 5 O 0.13639797948583 0.86360202051417 0.63639797948583 15.999 > 5 6 O 0.36360202051417 0.86360202051417 0.86360202051417 15.999 > 6 7 O 0.11360202051417 0.11360202051417 0.88639797948583 15.999 > 7 8 O 0.86360202051417 0.13639797948583 0.63639797948583 15.999 > 8 9 O 0.63639797948583 0.63639797948583 0.36360202051417 15.999 > 1 10 O 0.38639797948583 0.61360202051417 0.11360202051417 15.999 > 2 11 O 0.11360202051417 0.88639797948583 0.11360202051417 15.999 > 3 12 O 0.38639797948583 0.88639797948583 0.38639797948583 15.999 > 4 13 O 0.13639797948583 0.36360202051417 0.13639797948583 15.999 > 5 14 O 0.36360202051417 0.36360202051417 0.36360202051417 15.999 > 6 15 O 0.11360202051417 0.61360202051417 0.38639797948583 15.999 > 7 16 O 0.86360202051417 0.63639797948583 0.13639797948583 15.999 > 8 17 O 0.13639797948583 0.13639797948583 0.36360202051417 15.999 > 1 18 O 0.88639797948583 0.11360202051417 0.11360202051417 15.999 > 2 19 O 0.61360202051417 0.38639797948583 0.11360202051417 15.999 > 3 20 O 0.88639797948583 0.38639797948583 0.38639797948583 15.999 > 4 21 O 0.63639797948583 0.86360202051417 0.13639797948583 15.999 > 5 22 O 0.86360202051417 0.86360202051417 0.36360202051417 15.999 > 6 23 O 0.61360202051417 0.11360202051417 0.38639797948583 15.999 > 7 24 O 0.36360202051417 0.13639797948583 0.13639797948583 15.999 > 8 25 O 0.13639797948583 0.63639797948583 0.86360202051417 15.999 > 1 26 O 0.88639797948583 0.61360202051417 0.61360202051417 15.999 > 2 27 O 0.61360202051417 0.88639797948583 0.61360202051417 15.999 > 3 28 O 0.88639797948583 0.88639797948583 0.88639797948583 15.999 > 4 29 O 0.63639797948583 0.36360202051417 0.63639797948583 15.999 > 5 30 O 0.86360202051417 0.36360202051417 0.86360202051417 15.999 > 6 31 O 0.61360202051417 0.61360202051417 0.88639797948583 15.999 > 7 32 O 0.36360202051417 0.63639797948583 0.63639797948583 15.999 > 8 *33 Zn 0.25000000000000 0.25000000000000 0.75000000000000 65.380 > 9 34 Zn 0.50000000000000 0.00000000000000 0.00000000000000 65.380 > 10 35 Zn 0.25000000000000 0.75000000000000 0.25000000000000 65.380 > 9 36 Zn 0.50000000000000 0.50000000000000 0.50000000000000 65.380 > 10 37 Zn 0.75000000000000 0.25000000000000 0.25000000000000 65.380 > 9 38 Zn 0.00000000000000 0.00000000000000 0.50000000000000 65.380 > 10 39 Zn 0.75000000000000 0.75000000000000 0.75000000000000 65.380 > 9 40 Zn 0.00000000000000 0.50000000000000 0.00000000000000 65.380 > 10 *41 Ga 0.12500000000000 0.87500000000000 0.87500000000000 69.723 > 11 42 Ga 0.37500000000000 0.87500000000000 0.62500000000000 69.723 > 12 43 Ga 0.12500000000000 0.62500000000000 0.62500000000000 69.723 > 13 44 Ga 0.37500000000000 0.62500000000000 0.87500000000000 69.723 > 14 45 Ga 0.12500000000000 0.37500000000000 0.37500000000000 69.723 > 11 46 Ga 0.37500000000000 0.37500000000000 0.12500000000000 69.723 > 12 47 Ga 0.12500000000000 0.12500000000000 0.12500000000000 69.723 > 13 48 Ga 0.37500000000000 0.12500000000000 0.37500000000000 69.723 > 14 49 Ga 0.62500000000000 0.87500000000000 0.37500000000000 69.723 > 11 50 Ga 0.87500000000000 0.87500000000000 0.12500000000000 69.723 > 12 51 Ga 0.62500000000000 0.62500000000000 0.12500000000000 69.723 > 13 52 Ga 0.87500000000000 0.62500000000000 0.37500000000000 69.723 > 14 53 Ga 0.62500000000000 0.37500000000000 0.87500000000000 69.723 > 11 54 Ga 0.87500000000000 0.37500000000000 0.62500000000000 69.723 > 12 55 Ga 0.62500000000000 0.12500000000000 0.62500000000000 69.723 > 13 56 Ga 0.87500000000000 0.12500000000000 0.87500000000000 69.723 > 14 -------------------------------- super cell -------------------------------- Lattice vectors: a 8.351635040000000 0.000000000000000 0.000000000000000 b 0.000000000000000 8.351635040000000 0.000000000000000 c 0.000000000000000 0.000000000000000 8.351635040000000 Atomic positions (fractional): *1 O 0.63639797948583 0.13639797948583 0.86360202051417 15.999 > 1 2 O 0.38639797948583 0.11360202051417 0.61360202051417 15.999 > 2 3 O 0.11360202051417 0.38639797948583 0.61360202051417 15.999 > 3 4 O 0.38639797948583 0.38639797948583 0.88639797948583 15.999 > 4 5 O 0.13639797948583 0.86360202051417 0.63639797948583 15.999 > 5 6 O 0.36360202051417 0.86360202051417 0.86360202051417 15.999 > 6 7 O 0.11360202051417 0.11360202051417 0.88639797948583 15.999 > 7 8 O 0.86360202051417 0.13639797948583 0.63639797948583 15.999 > 8 9 O 0.63639797948583 0.63639797948583 0.36360202051417 15.999 > 1 10 O 0.38639797948583 0.61360202051417 0.11360202051417 15.999 > 2 11 O 0.11360202051417 0.88639797948583 0.11360202051417 15.999 > 3 12 O 0.38639797948583 0.88639797948583 0.38639797948583 15.999 > 4 13 O 0.13639797948583 0.36360202051417 0.13639797948583 15.999 > 5 14 O 0.36360202051417 0.36360202051417 0.36360202051417 15.999 > 6 15 O 0.11360202051417 0.61360202051417 0.38639797948583 15.999 > 7 16 O 0.86360202051417 0.63639797948583 0.13639797948583 15.999 > 8 17 O 0.13639797948583 0.13639797948583 0.36360202051417 15.999 > 1 18 O 0.88639797948583 0.11360202051417 0.11360202051417 15.999 > 2 19 O 0.61360202051417 0.38639797948583 0.11360202051417 15.999 > 3 20 O 0.88639797948583 0.38639797948583 0.38639797948583 15.999 > 4 21 O 0.63639797948583 0.86360202051417 0.13639797948583 15.999 > 5 22 O 0.86360202051417 0.86360202051417 0.36360202051417 15.999 > 6 23 O 0.61360202051417 0.11360202051417 0.38639797948583 15.999 > 7 24 O 0.36360202051417 0.13639797948583 0.13639797948583 15.999 > 8 25 O 0.13639797948583 0.63639797948583 0.86360202051417 15.999 > 1 26 O 0.88639797948583 0.61360202051417 0.61360202051417 15.999 > 2 27 O 0.61360202051417 0.88639797948583 0.61360202051417 15.999 > 3 28 O 0.88639797948583 0.88639797948583 0.88639797948583 15.999 > 4 29 O 0.63639797948583 0.36360202051417 0.63639797948583 15.999 > 5 30 O 0.86360202051417 0.36360202051417 0.86360202051417 15.999 > 6 31 O 0.61360202051417 0.61360202051417 0.88639797948583 15.999 > 7 32 O 0.36360202051417 0.63639797948583 0.63639797948583 15.999 > 8 *33 Zn 0.25000000000000 0.25000000000000 0.75000000000000 65.380 > 9 34 Zn 0.50000000000000 0.00000000000000 0.00000000000000 65.380 > 10 35 Zn 0.25000000000000 0.75000000000000 0.25000000000000 65.380 > 9 36 Zn 0.50000000000000 0.50000000000000 0.50000000000000 65.380 > 10 37 Zn 0.75000000000000 0.25000000000000 0.25000000000000 65.380 > 9 38 Zn 0.00000000000000 0.00000000000000 0.50000000000000 65.380 > 10 39 Zn 0.75000000000000 0.75000000000000 0.75000000000000 65.380 > 9 40 Zn 0.00000000000000 0.50000000000000 0.00000000000000 65.380 > 10 *41 Ga 0.12500000000000 0.87500000000000 0.87500000000000 69.723 > 11 42 Ga 0.37500000000000 0.87500000000000 0.62500000000000 69.723 > 12 43 Ga 0.12500000000000 0.62500000000000 0.62500000000000 69.723 > 13 44 Ga 0.37500000000000 0.62500000000000 0.87500000000000 69.723 > 14 45 Ga 0.12500000000000 0.37500000000000 0.37500000000000 69.723 > 11 46 Ga 0.37500000000000 0.37500000000000 0.12500000000000 69.723 > 12 47 Ga 0.12500000000000 0.12500000000000 0.12500000000000 69.723 > 13 48 Ga 0.37500000000000 0.12500000000000 0.37500000000000 69.723 > 14 49 Ga 0.62500000000000 0.87500000000000 0.37500000000000 69.723 > 11 50 Ga 0.87500000000000 0.87500000000000 0.12500000000000 69.723 > 12 51 Ga 0.62500000000000 0.62500000000000 0.12500000000000 69.723 > 13 52 Ga 0.87500000000000 0.62500000000000 0.37500000000000 69.723 > 14 53 Ga 0.62500000000000 0.37500000000000 0.87500000000000 69.723 > 11 54 Ga 0.87500000000000 0.37500000000000 0.62500000000000 69.723 > 12 55 Ga 0.62500000000000 0.12500000000000 0.62500000000000 69.723 > 13 56 Ga 0.87500000000000 0.12500000000000 0.87500000000000 69.723 > 14 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 4.1888409 0.0000000 0.0000000 0.0000000 4.1888409 0.0000000 0.0000000 0.0000000 4.1888409 -------------------------- Born effective charges -------------------------- 1 O -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 2 O -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 3 O -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 4 O -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 5 O -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 6 O -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 7 O -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 8 O -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 9 Zn 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 10 Zn 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 11 Ga 3.0620933 -0.2846752 -0.2846752 -0.2846752 3.0620933 0.2846752 -0.2846752 0.2846752 3.0620933 12 Ga 3.0620933 0.2846752 -0.2846752 0.2846752 3.0620933 -0.2846752 -0.2846752 -0.2846752 3.0620933 13 Ga 3.0620933 0.2846752 0.2846752 0.2846752 3.0620933 0.2846752 0.2846752 0.2846752 3.0620933 14 Ga 3.0620933 -0.2846752 0.2846752 -0.2846752 3.0620933 -0.2846752 0.2846752 -0.2846752 3.0620933 ---------------------------------------------------------------------------- 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: 168/168 Permutation basis: 4641/4641 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 103 Number of blocks in projector: 103 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 56 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 47 Use standard eigh solver. Tree of FC basis block matrices: - (103, 98), data: False |-- (47, 44), data: True |-- (56, 54), data: True ----- Solver_atoms: 1 -- 56 / 56 Time (Solver_compr_matrix_reshape): 0.001 Solver_block: 80 / 80 - Time: 0.016 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.018 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (xx) -0.00000000 (xx) Permutation basis: 168/168 Permutation basis: 4641/4641 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 103 Number of blocks in projector: 103 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 56 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 47 Use standard eigh solver. Tree of FC basis block matrices: - (103, 98), data: False |-- (47, 44), data: True |-- (56, 54), data: True Max drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 21:02:21]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 21:02:22]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [1 1 1] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fd-3m (227) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 4.175817520000000 4.175817520000000 b 4.175817520000000 0.000000000000000 4.175817520000000 c 4.175817520000000 4.175817520000000 0.000000000000000 Atomic positions (fractional): 1 O 0.36360202051417 0.36360202051417 0.90919393845749 15.999 2 O 0.34080606154251 0.88639797948583 0.88639797948583 15.999 3 O 0.88639797948583 0.34080606154251 0.88639797948583 15.999 4 O 0.88639797948583 0.88639797948583 0.88639797948583 15.999 5 O 0.36360202051417 0.90919393845749 0.36360202051417 15.999 6 O 0.36360202051417 0.36360202051417 0.36360202051417 15.999 7 O 0.88639797948583 0.88639797948583 0.34080606154251 15.999 8 O 0.90919393845749 0.36360202051417 0.36360202051417 15.999 9 Zn 0.75000000000000 0.75000000000000 0.75000000000000 65.380 10 Zn 0.50000000000000 0.50000000000000 0.50000000000000 65.380 11 Ga 0.62500000000000 0.12500000000000 0.12500000000000 69.723 12 Ga 0.12500000000000 0.12500000000000 0.62500000000000 69.723 13 Ga 0.12500000000000 0.12500000000000 0.12500000000000 69.723 14 Ga 0.12500000000000 0.62500000000000 0.12500000000000 69.723 -------------------------------- supercell --------------------------------- Lattice vectors: a 8.351635040000000 0.000000000000000 0.000000000000000 b 0.000000000000000 8.351635040000000 0.000000000000000 c 0.000000000000000 0.000000000000000 8.351635040000000 Atomic positions (fractional): 1 O 0.63639797948583 0.13639797948583 0.86360202051417 15.999 > 1 2 O 0.38639797948583 0.11360202051417 0.61360202051417 15.999 > 2 3 O 0.11360202051417 0.38639797948583 0.61360202051417 15.999 > 3 4 O 0.38639797948583 0.38639797948583 0.88639797948583 15.999 > 4 5 O 0.13639797948583 0.86360202051417 0.63639797948583 15.999 > 5 6 O 0.36360202051417 0.86360202051417 0.86360202051417 15.999 > 6 7 O 0.11360202051417 0.11360202051417 0.88639797948583 15.999 > 7 8 O 0.86360202051417 0.13639797948583 0.63639797948583 15.999 > 8 9 O 0.63639797948583 0.63639797948583 0.36360202051417 15.999 > 1 10 O 0.38639797948583 0.61360202051417 0.11360202051417 15.999 > 2 11 O 0.11360202051417 0.88639797948583 0.11360202051417 15.999 > 3 12 O 0.38639797948583 0.88639797948583 0.38639797948583 15.999 > 4 13 O 0.13639797948583 0.36360202051417 0.13639797948583 15.999 > 5 14 O 0.36360202051417 0.36360202051417 0.36360202051417 15.999 > 6 15 O 0.11360202051417 0.61360202051417 0.38639797948583 15.999 > 7 16 O 0.86360202051417 0.63639797948583 0.13639797948583 15.999 > 8 17 O 0.13639797948583 0.13639797948583 0.36360202051417 15.999 > 1 18 O 0.88639797948583 0.11360202051417 0.11360202051417 15.999 > 2 19 O 0.61360202051417 0.38639797948583 0.11360202051417 15.999 > 3 20 O 0.88639797948583 0.38639797948583 0.38639797948583 15.999 > 4 21 O 0.63639797948583 0.86360202051417 0.13639797948583 15.999 > 5 22 O 0.86360202051417 0.86360202051417 0.36360202051417 15.999 > 6 23 O 0.61360202051417 0.11360202051417 0.38639797948583 15.999 > 7 24 O 0.36360202051417 0.13639797948583 0.13639797948583 15.999 > 8 25 O 0.13639797948583 0.63639797948583 0.86360202051417 15.999 > 1 26 O 0.88639797948583 0.61360202051417 0.61360202051417 15.999 > 2 27 O 0.61360202051417 0.88639797948583 0.61360202051417 15.999 > 3 28 O 0.88639797948583 0.88639797948583 0.88639797948583 15.999 > 4 29 O 0.63639797948583 0.36360202051417 0.63639797948583 15.999 > 5 30 O 0.86360202051417 0.36360202051417 0.86360202051417 15.999 > 6 31 O 0.61360202051417 0.61360202051417 0.88639797948583 15.999 > 7 32 O 0.36360202051417 0.63639797948583 0.63639797948583 15.999 > 8 33 Zn 0.25000000000000 0.25000000000000 0.75000000000000 65.380 > 33 34 Zn 0.50000000000000 0.00000000000000 0.00000000000000 65.380 > 34 35 Zn 0.25000000000000 0.75000000000000 0.25000000000000 65.380 > 33 36 Zn 0.50000000000000 0.50000000000000 0.50000000000000 65.380 > 34 37 Zn 0.75000000000000 0.25000000000000 0.25000000000000 65.380 > 33 38 Zn 0.00000000000000 0.00000000000000 0.50000000000000 65.380 > 34 39 Zn 0.75000000000000 0.75000000000000 0.75000000000000 65.380 > 33 40 Zn 0.00000000000000 0.50000000000000 0.00000000000000 65.380 > 34 41 Ga 0.12500000000000 0.87500000000000 0.87500000000000 69.723 > 41 42 Ga 0.37500000000000 0.87500000000000 0.62500000000000 69.723 > 42 43 Ga 0.12500000000000 0.62500000000000 0.62500000000000 69.723 > 43 44 Ga 0.37500000000000 0.62500000000000 0.87500000000000 69.723 > 44 45 Ga 0.12500000000000 0.37500000000000 0.37500000000000 69.723 > 41 46 Ga 0.37500000000000 0.37500000000000 0.12500000000000 69.723 > 42 47 Ga 0.12500000000000 0.12500000000000 0.12500000000000 69.723 > 43 48 Ga 0.37500000000000 0.12500000000000 0.37500000000000 69.723 > 44 49 Ga 0.62500000000000 0.87500000000000 0.37500000000000 69.723 > 41 50 Ga 0.87500000000000 0.87500000000000 0.12500000000000 69.723 > 42 51 Ga 0.62500000000000 0.62500000000000 0.12500000000000 69.723 > 43 52 Ga 0.87500000000000 0.62500000000000 0.37500000000000 69.723 > 44 53 Ga 0.62500000000000 0.37500000000000 0.87500000000000 69.723 > 41 54 Ga 0.87500000000000 0.37500000000000 0.62500000000000 69.723 > 42 55 Ga 0.62500000000000 0.12500000000000 0.62500000000000 69.723 > 43 56 Ga 0.87500000000000 0.12500000000000 0.87500000000000 69.723 > 44 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 4.1888409 0.0000000 0.0000000 0.0000000 4.1888409 0.0000000 0.0000000 0.0000000 4.1888409 -------------------------- Born effective charges -------------------------- 1 O -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 2 O -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 3 O -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 4 O -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 5 O -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 6 O -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 7 O -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 8 O -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 9 Zn 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 10 Zn 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 11 Ga 3.0620933 -0.2846752 -0.2846752 -0.2846752 3.0620933 0.2846752 -0.2846752 0.2846752 3.0620933 12 Ga 3.0620933 0.2846752 -0.2846752 0.2846752 3.0620933 -0.2846752 -0.2846752 -0.2846752 3.0620933 13 Ga 3.0620933 0.2846752 0.2846752 0.2846752 3.0620933 0.2846752 0.2846752 0.2846752 3.0620933 14 Ga 3.0620933 -0.2846752 0.2846752 -0.2846752 3.0620933 -0.2846752 0.2846752 -0.2846752 3.0620933 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 41, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000003 (xxx) -0.00000003 (xxx) -0.00000003 (xxx) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 21:02:23]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 21:02:24]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [1 1 1] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fd-3m (227) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 4.175817520000000 4.175817520000000 b 4.175817520000000 0.000000000000000 4.175817520000000 c 4.175817520000000 4.175817520000000 0.000000000000000 Atomic positions (fractional): 1 O 0.36360202051417 0.36360202051417 0.90919393845749 15.999 2 O 0.34080606154251 0.88639797948583 0.88639797948583 15.999 3 O 0.88639797948583 0.34080606154251 0.88639797948583 15.999 4 O 0.88639797948583 0.88639797948583 0.88639797948583 15.999 5 O 0.36360202051417 0.90919393845749 0.36360202051417 15.999 6 O 0.36360202051417 0.36360202051417 0.36360202051417 15.999 7 O 0.88639797948583 0.88639797948583 0.34080606154251 15.999 8 O 0.90919393845749 0.36360202051417 0.36360202051417 15.999 9 Zn 0.75000000000000 0.75000000000000 0.75000000000000 65.380 10 Zn 0.50000000000000 0.50000000000000 0.50000000000000 65.380 11 Ga 0.62500000000000 0.12500000000000 0.12500000000000 69.723 12 Ga 0.12500000000000 0.12500000000000 0.62500000000000 69.723 13 Ga 0.12500000000000 0.12500000000000 0.12500000000000 69.723 14 Ga 0.12500000000000 0.62500000000000 0.12500000000000 69.723 -------------------------------- supercell --------------------------------- Lattice vectors: a 8.351635040000000 0.000000000000000 0.000000000000000 b 0.000000000000000 8.351635040000000 0.000000000000000 c 0.000000000000000 0.000000000000000 8.351635040000000 Atomic positions (fractional): 1 O 0.63639797948583 0.13639797948583 0.86360202051417 15.999 > 1 2 O 0.38639797948583 0.11360202051417 0.61360202051417 15.999 > 2 3 O 0.11360202051417 0.38639797948583 0.61360202051417 15.999 > 3 4 O 0.38639797948583 0.38639797948583 0.88639797948583 15.999 > 4 5 O 0.13639797948583 0.86360202051417 0.63639797948583 15.999 > 5 6 O 0.36360202051417 0.86360202051417 0.86360202051417 15.999 > 6 7 O 0.11360202051417 0.11360202051417 0.88639797948583 15.999 > 7 8 O 0.86360202051417 0.13639797948583 0.63639797948583 15.999 > 8 9 O 0.63639797948583 0.63639797948583 0.36360202051417 15.999 > 1 10 O 0.38639797948583 0.61360202051417 0.11360202051417 15.999 > 2 11 O 0.11360202051417 0.88639797948583 0.11360202051417 15.999 > 3 12 O 0.38639797948583 0.88639797948583 0.38639797948583 15.999 > 4 13 O 0.13639797948583 0.36360202051417 0.13639797948583 15.999 > 5 14 O 0.36360202051417 0.36360202051417 0.36360202051417 15.999 > 6 15 O 0.11360202051417 0.61360202051417 0.38639797948583 15.999 > 7 16 O 0.86360202051417 0.63639797948583 0.13639797948583 15.999 > 8 17 O 0.13639797948583 0.13639797948583 0.36360202051417 15.999 > 1 18 O 0.88639797948583 0.11360202051417 0.11360202051417 15.999 > 2 19 O 0.61360202051417 0.38639797948583 0.11360202051417 15.999 > 3 20 O 0.88639797948583 0.38639797948583 0.38639797948583 15.999 > 4 21 O 0.63639797948583 0.86360202051417 0.13639797948583 15.999 > 5 22 O 0.86360202051417 0.86360202051417 0.36360202051417 15.999 > 6 23 O 0.61360202051417 0.11360202051417 0.38639797948583 15.999 > 7 24 O 0.36360202051417 0.13639797948583 0.13639797948583 15.999 > 8 25 O 0.13639797948583 0.63639797948583 0.86360202051417 15.999 > 1 26 O 0.88639797948583 0.61360202051417 0.61360202051417 15.999 > 2 27 O 0.61360202051417 0.88639797948583 0.61360202051417 15.999 > 3 28 O 0.88639797948583 0.88639797948583 0.88639797948583 15.999 > 4 29 O 0.63639797948583 0.36360202051417 0.63639797948583 15.999 > 5 30 O 0.86360202051417 0.36360202051417 0.86360202051417 15.999 > 6 31 O 0.61360202051417 0.61360202051417 0.88639797948583 15.999 > 7 32 O 0.36360202051417 0.63639797948583 0.63639797948583 15.999 > 8 33 Zn 0.25000000000000 0.25000000000000 0.75000000000000 65.380 > 33 34 Zn 0.50000000000000 0.00000000000000 0.00000000000000 65.380 > 34 35 Zn 0.25000000000000 0.75000000000000 0.25000000000000 65.380 > 33 36 Zn 0.50000000000000 0.50000000000000 0.50000000000000 65.380 > 34 37 Zn 0.75000000000000 0.25000000000000 0.25000000000000 65.380 > 33 38 Zn 0.00000000000000 0.00000000000000 0.50000000000000 65.380 > 34 39 Zn 0.75000000000000 0.75000000000000 0.75000000000000 65.380 > 33 40 Zn 0.00000000000000 0.50000000000000 0.00000000000000 65.380 > 34 41 Ga 0.12500000000000 0.87500000000000 0.87500000000000 69.723 > 41 42 Ga 0.37500000000000 0.87500000000000 0.62500000000000 69.723 > 42 43 Ga 0.12500000000000 0.62500000000000 0.62500000000000 69.723 > 43 44 Ga 0.37500000000000 0.62500000000000 0.87500000000000 69.723 > 44 45 Ga 0.12500000000000 0.37500000000000 0.37500000000000 69.723 > 41 46 Ga 0.37500000000000 0.37500000000000 0.12500000000000 69.723 > 42 47 Ga 0.12500000000000 0.12500000000000 0.12500000000000 69.723 > 43 48 Ga 0.37500000000000 0.12500000000000 0.37500000000000 69.723 > 44 49 Ga 0.62500000000000 0.87500000000000 0.37500000000000 69.723 > 41 50 Ga 0.87500000000000 0.87500000000000 0.12500000000000 69.723 > 42 51 Ga 0.62500000000000 0.62500000000000 0.12500000000000 69.723 > 43 52 Ga 0.87500000000000 0.62500000000000 0.37500000000000 69.723 > 44 53 Ga 0.62500000000000 0.37500000000000 0.87500000000000 69.723 > 41 54 Ga 0.87500000000000 0.37500000000000 0.62500000000000 69.723 > 42 55 Ga 0.62500000000000 0.12500000000000 0.62500000000000 69.723 > 43 56 Ga 0.87500000000000 0.12500000000000 0.87500000000000 69.723 > 44 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 4.1888409 0.0000000 0.0000000 0.0000000 4.1888409 0.0000000 0.0000000 0.0000000 4.1888409 -------------------------- Born effective charges -------------------------- 1 O -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 2 O -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 3 O -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 4 O -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 5 O -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 0.1251226 -0.1251226 0.1251226 -2.0914624 6 O -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 -0.1251226 -0.1251226 -0.1251226 -2.0914624 7 O -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 8 O -2.0914624 0.1251226 0.1251226 0.1251226 -2.0914624 -0.1251226 0.1251226 -0.1251226 -2.0914624 9 Zn 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 10 Zn 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 0.0000000 0.0000000 0.0000000 2.2416631 11 Ga 3.0620933 -0.2846752 -0.2846752 -0.2846752 3.0620933 0.2846752 -0.2846752 0.2846752 3.0620933 12 Ga 3.0620933 0.2846752 -0.2846752 0.2846752 3.0620933 -0.2846752 -0.2846752 -0.2846752 3.0620933 13 Ga 3.0620933 0.2846752 0.2846752 0.2846752 3.0620933 0.2846752 0.2846752 0.2846752 3.0620933 14 Ga 3.0620933 -0.2846752 0.2846752 -0.2846752 3.0620933 -0.2846752 0.2846752 -0.2846752 3.0620933 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000003 (xxx) -0.00000003 (xxx) -0.00000003 (xxx) Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 10 10 10 ] 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.79, Number of G-points: 307, Lambda: 0.17 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/47) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 47 Calculating Gamma of ph-isotope with tetrahedron method Frequency 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 3.774 ( -0.000 -0.000 0.000) 0.000 3.774 ( 0.000 0.000 0.000) 0.000 3.774 ( -0.000 0.000 0.000) 0.000 5.045 ( 0.000 0.000 -0.000) 0.000 5.045 ( -0.000 0.000 -0.000) 0.000 5.045 ( 0.000 0.000 -0.000) 0.000 5.353 ( 0.000 -0.000 0.000) 0.000 5.353 ( -0.000 0.000 -0.000) 0.000 5.353 ( 0.000 0.000 -0.000) 0.000 6.246 ( 0.000 0.000 0.000) 0.000 6.246 ( 0.000 0.000 -0.000) 0.000 9.571 ( 0.000 -0.000 -0.000) 0.000 9.571 ( 0.000 -0.000 -0.000) 0.000 9.571 ( 0.000 -0.000 0.000) 0.000 9.662 ( 0.000 0.000 -0.000) 0.000 9.662 ( 0.000 -0.000 -0.000) 0.000 9.662 ( 0.000 0.000 -0.000) 0.000 10.570 ( -0.000 0.000 0.000) 0.000 10.570 ( -0.000 -0.000 0.000) 0.000 11.414 ( -0.000 -0.000 0.000) 0.000 11.414 ( 0.000 -0.000 0.000) 0.000 11.414 ( 0.000 -0.000 0.000) 0.000 11.743 ( 0.000 0.000 -0.000) 0.000 12.144 ( 0.000 -0.000 0.000) 0.000 12.144 ( -0.000 0.000 0.000) 0.000 12.144 ( 0.000 0.000 0.000) 0.000 13.526 ( -0.000 0.000 0.000) 0.000 13.526 ( 0.000 -0.000 0.000) 0.000 13.526 ( -0.000 0.000 -0.000) 0.000 15.870 ( 0.000 0.000 0.000) 0.000 15.870 ( 0.000 -0.000 0.000) 0.000 16.279 ( -0.000 0.000 0.000) 0.000 16.279 ( 0.000 0.000 -0.000) 0.000 16.279 ( -0.000 0.000 -0.000) 0.000 17.438 ( 0.000 -0.000 0.000) 0.000 17.438 ( 0.000 0.000 0.000) 0.000 17.438 ( 0.000 -0.000 0.000) 0.000 20.124 ( -0.000 0.000 0.000) 0.000 20.461 ( -0.000 0.000 -0.000) 0.000 ======================= Grid point 1 (2/47) ======================= q-point: ( 0.10 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 110 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.710 ( -19.308 19.308 19.308) 33.442 0.710 ( -19.308 19.308 19.308) 33.442 1.510 ( -41.403 41.403 41.403) 71.712 3.756 ( 0.963 -0.963 -0.963) 1.668 3.835 ( -2.943 2.943 2.943) 5.097 3.835 ( -2.943 2.943 2.943) 5.097 5.074 ( -1.722 1.722 1.722) 2.982 5.074 ( -1.722 1.722 1.722) 2.982 5.174 ( -4.085 4.085 4.085) 7.075 5.227 ( 5.842 -5.842 -5.842) 10.118 5.227 ( 5.842 -5.842 -5.842) 10.118 5.344 ( 0.383 -0.383 -0.383) 0.663 6.371 ( -5.936 5.936 5.936) 10.281 6.371 ( -5.936 5.936 5.936) 10.281 9.567 ( 2.503 -2.503 -2.503) 4.335 9.567 ( 2.503 -2.503 -2.503) 4.335 9.658 ( -4.679 4.679 4.679) 8.105 9.662 ( -1.920 1.920 1.920) 3.326 9.662 ( -1.920 1.920 1.920) 3.326 9.820 ( 5.800 -5.800 -5.800) 10.046 10.622 ( -3.009 3.009 3.009) 5.212 10.622 ( -3.009 3.009 3.009) 5.212 11.431 ( -0.931 0.931 0.931) 1.613 11.431 ( -0.931 0.931 0.931) 1.613 11.636 ( 5.728 -5.728 -5.728) 9.921 12.117 ( 1.528 -1.528 -1.528) 2.647 12.214 ( -3.736 3.736 3.736) 6.470 12.214 ( -3.736 3.736 3.736) 6.470 13.592 ( -3.458 3.458 3.458) 5.990 13.594 ( -3.690 3.690 3.690) 6.390 13.594 ( -3.690 3.690 3.690) 6.390 15.172 ( -4.239 4.239 4.239) 7.343 15.833 ( 2.098 -2.098 -2.098) 3.633 15.833 ( 2.098 -2.098 -2.098) 3.633 16.273 ( 0.174 -0.174 -0.174) 0.302 16.273 ( 0.174 -0.174 -0.174) 0.302 17.411 ( 1.539 -1.539 -1.539) 2.665 17.411 ( 1.539 -1.539 -1.539) 2.665 17.501 ( -3.375 3.375 3.375) 5.845 19.303 ( -0.352 0.352 0.352) 0.610 20.107 ( 0.843 -0.843 -0.843) 1.460 20.457 ( 0.201 -0.201 -0.201) 0.348 ======================= Grid point 2 (3/47) ======================= q-point: ( 0.20 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 116 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.373 ( -17.403 17.403 17.403) 30.143 1.373 ( -17.403 17.403 17.403) 30.143 2.948 ( -38.185 38.185 38.185) 66.138 3.710 ( 1.487 -1.487 -1.487) 2.576 3.925 ( -1.022 1.022 1.022) 1.770 3.925 ( -1.022 1.022 1.022) 1.770 5.005 ( 5.409 -5.409 -5.409) 9.368 5.005 ( 5.409 -5.409 -5.409) 9.368 5.197 ( -5.068 5.068 5.068) 8.779 5.197 ( -5.068 5.068 5.068) 8.779 5.282 ( 0.286 -0.286 -0.286) 0.495 5.425 ( -7.014 7.014 7.014) 12.148 6.607 ( -6.695 6.695 6.695) 11.595 6.607 ( -6.695 6.695 6.695) 11.595 9.410 ( 5.603 -5.603 -5.603) 9.704 9.410 ( 5.603 -5.603 -5.603) 9.704 9.508 ( 11.549 -11.549 -11.549) 20.004 9.746 ( -1.714 1.714 1.714) 2.969 9.746 ( -1.714 1.714 1.714) 2.969 9.892 ( -8.045 8.045 8.045) 13.934 10.804 ( -7.251 7.251 7.251) 12.560 10.804 ( -7.251 7.251 7.251) 12.560 11.362 ( 8.999 -8.999 -8.999) 15.586 11.479 ( -1.722 1.722 1.722) 2.982 11.479 ( -1.722 1.722 1.722) 2.982 12.034 ( 3.054 -3.054 -3.054) 5.290 12.391 ( -5.702 5.702 5.702) 9.876 12.391 ( -5.702 5.702 5.702) 9.876 13.751 ( -5.011 5.011 5.011) 8.680 13.783 ( -6.694 6.694 6.694) 11.594 13.783 ( -6.694 6.694 6.694) 11.594 15.388 ( -7.520 7.520 7.520) 13.026 15.712 ( 4.716 -4.716 -4.716) 8.169 15.712 ( 4.716 -4.716 -4.716) 8.169 16.281 ( -0.842 0.842 0.842) 1.458 16.281 ( -0.842 0.842 0.842) 1.458 17.325 ( 3.242 -3.242 -3.242) 5.615 17.325 ( 3.242 -3.242 -3.242) 5.615 17.667 ( -5.625 5.625 5.625) 9.742 19.306 ( 0.510 -0.510 -0.510) 0.884 20.075 ( 0.674 -0.674 -0.674) 1.167 20.443 ( 0.635 -0.635 -0.635) 1.099 ======================= Grid point 3 (4/47) ======================= q-point: ( 0.30 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 110 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.950 ( -14.609 14.609 14.609) 25.304 1.950 ( -14.609 14.609 14.609) 25.304 3.656 ( 1.411 -1.411 -1.411) 2.444 3.863 ( 4.346 -4.346 -4.346) 7.528 3.863 ( 4.346 -4.346 -4.346) 7.528 4.194 ( -29.561 29.561 29.561) 51.201 4.860 ( 2.655 -2.655 -2.655) 4.599 4.860 ( 2.655 -2.655 -2.655) 4.599 5.266 ( -0.532 0.532 0.532) 0.922 5.444 ( -8.141 8.141 8.141) 14.101 5.444 ( -8.141 8.141 8.141) 14.101 5.754 ( -11.193 11.193 11.193) 19.387 6.831 ( -5.633 5.633 5.633) 9.756 6.831 ( -5.633 5.633 5.633) 9.756 8.993 ( 17.031 -17.031 -17.031) 29.499 9.180 ( 6.980 -6.980 -6.980) 12.089 9.180 ( 6.980 -6.980 -6.980) 12.089 9.767 ( 0.263 -0.263 -0.263) 0.456 9.767 ( 0.263 -0.263 -0.263) 0.456 10.209 ( -9.158 9.158 9.158) 15.862 11.031 ( 8.818 -8.818 -8.818) 15.273 11.117 ( -9.332 9.332 9.332) 16.164 11.117 ( -9.332 9.332 9.332) 16.164 11.551 ( -2.253 2.253 2.253) 3.903 11.551 ( -2.253 2.253 2.253) 3.903 11.900 ( 4.287 -4.287 -4.287) 7.426 12.595 ( -5.305 5.305 5.305) 9.189 12.595 ( -5.305 5.305 5.305) 9.189 13.928 ( -4.564 4.564 4.564) 7.905 14.056 ( -8.211 8.211 8.211) 14.223 14.056 ( -8.211 8.211 8.211) 14.223 15.498 ( 6.993 -6.993 -6.993) 12.112 15.498 ( 6.993 -6.993 -6.993) 12.112 15.681 ( -8.260 8.260 8.260) 14.307 16.345 ( -2.809 2.809 2.809) 4.866 16.345 ( -2.809 2.809 2.809) 4.866 17.178 ( 4.931 -4.931 -4.931) 8.540 17.178 ( 4.931 -4.931 -4.931) 8.540 17.887 ( -6.380 6.380 6.380) 11.050 19.248 ( 2.792 -2.792 -2.792) 4.835 20.076 ( -0.878 0.878 0.878) 1.521 20.409 ( 1.327 -1.327 -1.327) 2.298 ======================= Grid point 4 (5/47) ======================= q-point: ( 0.40 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 116 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.404 ( -10.190 10.190 10.190) 17.649 2.404 ( -10.190 10.190 10.190) 17.649 3.614 ( 0.839 -0.839 -0.839) 1.452 3.659 ( 6.043 -6.043 -6.043) 10.467 3.659 ( 6.043 -6.043 -6.043) 10.467 4.799 ( 0.945 -0.945 -0.945) 1.638 4.799 ( 0.945 -0.945 -0.945) 1.638 4.867 ( -6.868 6.868 6.868) 11.895 5.418 ( -7.507 7.507 7.507) 13.003 5.721 ( -6.314 6.314 6.314) 10.937 5.721 ( -6.314 6.314 6.314) 10.937 6.283 ( -18.628 18.628 18.628) 32.265 7.002 ( -3.704 3.704 3.704) 6.416 7.002 ( -3.704 3.704 3.704) 6.416 8.308 ( 20.271 -20.271 -20.271) 35.110 8.940 ( 5.738 -5.738 -5.738) 9.938 8.940 ( 5.738 -5.738 -5.738) 9.938 9.751 ( 0.456 -0.456 -0.456) 0.789 9.751 ( 0.456 -0.456 -0.456) 0.789 10.511 ( -6.991 6.991 6.991) 12.108 10.768 ( 5.395 -5.395 -5.395) 9.345 11.398 ( -5.343 5.343 5.343) 9.254 11.398 ( -5.343 5.343 5.343) 9.254 11.641 ( -2.667 2.667 2.667) 4.619 11.641 ( -2.667 2.667 2.667) 4.619 11.743 ( 4.036 -4.036 -4.036) 6.990 12.751 ( -3.187 3.187 3.187) 5.519 12.751 ( -3.187 3.187 3.187) 5.519 14.062 ( -2.670 2.670 2.670) 4.625 14.347 ( -7.439 7.439 7.439) 12.885 14.347 ( -7.439 7.439 7.439) 12.885 15.233 ( 7.188 -7.188 -7.188) 12.451 15.233 ( 7.188 -7.188 -7.188) 12.451 15.938 ( -5.461 5.461 5.461) 9.459 16.482 ( -4.585 4.585 4.585) 7.941 16.482 ( -4.585 4.585 4.585) 7.941 16.980 ( 5.797 -5.797 -5.797) 10.040 16.980 ( 5.797 -5.797 -5.797) 10.040 18.102 ( -5.135 5.135 5.135) 8.894 19.120 ( 3.749 -3.749 -3.749) 6.493 20.135 ( -2.146 2.146 2.146) 3.718 20.349 ( 1.940 -1.940 -1.940) 3.361 ======================= Grid point 5 (6/47) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 58 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.611 ( -0.000 0.000 0.000) 0.000 2.611 ( -0.000 0.000 0.000) 0.000 3.524 ( -0.000 0.000 0.000) 0.000 3.524 ( -0.000 0.000 0.000) 0.000 3.599 ( -0.000 0.000 0.000) 0.000 4.783 ( -0.000 0.000 0.000) 0.000 4.783 ( -0.000 0.000 0.000) 0.000 4.946 ( -0.000 0.000 0.000) 0.000 5.576 ( -0.000 0.000 0.000) 0.000 5.842 ( -0.000 0.000 0.000) 0.000 5.842 ( -0.000 0.000 0.000) 0.000 6.806 ( 0.000 -0.000 -0.000) 0.000 7.074 ( -0.000 0.000 0.000) 0.000 7.074 ( -0.000 0.000 0.000) 0.000 7.778 ( -0.000 0.000 0.000) 0.000 8.825 ( -0.000 0.000 0.000) 0.000 8.825 ( -0.000 0.000 0.000) 0.000 9.742 ( -0.000 0.000 0.000) 0.000 9.742 ( -0.000 0.000 0.000) 0.000 10.650 ( -0.000 0.000 0.000) 0.000 10.668 ( -0.000 0.000 0.000) 0.000 11.487 ( -0.000 0.000 0.000) 0.000 11.487 ( -0.000 0.000 0.000) 0.000 11.658 ( -0.000 0.000 0.000) 0.000 11.702 ( -0.000 0.000 0.000) 0.000 11.702 ( -0.000 0.000 0.000) 0.000 12.810 ( -0.000 0.000 0.000) 0.000 12.810 ( -0.000 0.000 0.000) 0.000 14.111 ( -0.000 0.000 0.000) 0.000 14.513 ( -0.000 0.000 0.000) 0.000 14.513 ( -0.000 0.000 0.000) 0.000 15.070 ( 0.000 -0.000 -0.000) 0.000 15.070 ( 0.000 -0.000 -0.000) 0.000 16.041 ( -0.000 0.000 0.000) 0.000 16.611 ( -0.000 0.000 0.000) 0.000 16.611 ( -0.000 0.000 0.000) 0.000 16.828 ( -0.000 0.000 0.000) 0.000 16.828 ( -0.000 0.000 0.000) 0.000 18.207 ( -0.000 0.000 0.000) 0.000 19.039 ( 0.000 -0.000 -0.000) 0.000 20.193 ( 0.000 -0.000 -0.000) 0.000 20.296 ( -0.000 0.000 0.000) 0.000 ======================= Grid point 12 (7/47) ======================= q-point: ( 0.10 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 99 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.977 ( 0.000 -0.000 39.732) 39.732 0.977 ( 0.000 -0.000 39.732) 39.732 1.570 ( 0.000 -0.000 64.474) 64.474 3.803 ( 0.000 -0.000 1.917) 1.917 3.803 ( 0.000 -0.000 1.917) 1.917 3.860 ( 0.000 -0.000 7.047) 7.047 5.055 ( -0.000 0.000 -17.835) 17.835 5.110 ( 0.000 -0.000 5.572) 5.572 5.110 ( 0.000 -0.000 5.572) 5.572 5.148 ( 0.000 -0.000 4.140) 4.140 5.352 ( 0.000 -0.000 0.114) 0.114 5.352 ( 0.000 -0.000 0.114) 0.114 6.203 ( -0.000 0.000 -3.533) 3.533 6.543 ( 0.000 -0.000 17.867) 17.867 9.571 ( -0.000 0.000 -6.085) 6.085 9.571 ( -0.000 0.000 -6.085) 6.085 9.589 ( 0.000 -0.000 1.473) 1.473 9.712 ( 0.000 -0.000 9.479) 9.479 9.712 ( 0.000 -0.000 9.479) 9.479 9.811 ( -0.000 0.000 -8.721) 8.721 10.491 ( -0.000 0.000 -5.881) 5.881 10.797 ( 0.000 -0.000 18.057) 18.057 11.440 ( 0.000 -0.000 2.415) 2.415 11.440 ( 0.000 -0.000 2.415) 2.415 11.597 ( -0.000 0.000 -11.901) 11.901 12.127 ( -0.000 0.000 -1.502) 1.502 12.127 ( -0.000 0.000 -1.502) 1.502 12.321 ( 0.000 -0.000 13.918) 13.918 13.638 ( 0.000 -0.000 9.356) 9.356 13.638 ( 0.000 -0.000 9.356) 9.356 13.671 ( 0.000 -0.000 11.656) 11.656 15.171 ( 0.000 -0.000 6.271) 6.271 15.810 ( -0.000 0.000 -4.969) 4.969 15.837 ( -0.000 0.000 -2.619) 2.619 16.259 ( -0.000 0.000 -1.508) 1.508 16.259 ( -0.000 0.000 -1.508) 1.508 17.415 ( -0.000 0.000 -2.132) 2.132 17.415 ( -0.000 0.000 -2.132) 2.132 17.496 ( 0.000 -0.000 4.745) 4.745 19.254 ( -0.000 0.000 -3.347) 3.347 20.100 ( -0.000 0.000 -1.955) 1.955 20.455 ( -0.000 0.000 -0.651) 0.651 ======================= Grid point 13 (8/47) ======================= q-point: ( 0.20 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.417 ( -1.334 1.334 37.227) 37.275 1.511 ( -5.576 5.576 34.961) 35.840 2.691 ( -25.993 25.993 52.308) 63.933 3.781 ( 2.928 -2.928 2.190) 4.684 3.871 ( -3.139 3.139 1.581) 4.712 3.972 ( -0.176 0.176 6.708) 6.712 4.844 ( 0.125 -0.125 -15.244) 15.245 5.137 ( 0.730 -0.730 4.839) 4.948 5.176 ( 3.470 -3.470 1.446) 5.116 5.270 ( -2.032 2.032 1.573) 3.277 5.291 ( 2.148 -2.148 9.320) 9.802 5.418 ( -4.042 4.042 8.251) 10.038 6.287 ( -9.851 9.851 -4.998) 14.800 6.768 ( -1.456 1.456 15.776) 15.910 9.420 ( 3.819 -3.819 -10.968) 12.225 9.428 ( 3.212 -3.212 -10.056) 11.034 9.607 ( 6.675 -6.675 -10.638) 14.222 9.669 ( -4.024 4.024 1.559) 5.901 9.842 ( 1.057 -1.057 11.179) 11.279 9.876 ( -2.377 2.377 11.929) 12.393 10.505 ( -7.474 7.474 -5.582) 11.954 11.038 ( 0.571 -0.571 24.213) 24.227 11.364 ( 4.211 -4.211 -15.127) 16.256 11.485 ( -0.583 0.583 3.992) 4.077 11.508 ( -1.013 1.013 4.216) 4.453 12.071 ( 2.485 -2.485 -2.695) 4.429 12.175 ( -5.844 5.844 -2.443) 8.618 12.538 ( -1.489 1.489 17.979) 18.101 13.712 ( 1.367 -1.367 10.932) 11.101 13.813 ( -2.872 2.872 13.744) 14.331 13.895 ( -2.083 2.083 16.554) 16.814 15.329 ( -5.135 5.135 8.786) 11.399 15.703 ( 2.650 -2.650 -7.935) 8.775 15.770 ( 2.530 -2.530 -3.695) 5.144 16.253 ( -1.049 1.049 -1.448) 2.073 16.253 ( -1.171 1.171 -1.159) 2.022 17.341 ( 2.541 -2.541 -3.964) 5.350 17.356 ( 1.938 -1.938 -4.035) 4.878 17.622 ( -4.047 4.047 6.439) 8.615 19.248 ( -2.600 2.600 -4.558) 5.856 20.064 ( 0.436 -0.436 -2.688) 2.757 20.440 ( 0.205 -0.205 -1.534) 1.561 ======================= Grid point 14 (9/47) ======================= q-point: ( 0.30 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 280 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.920 ( -5.570 5.570 31.351) 32.326 2.039 ( -6.196 6.196 28.700) 30.008 3.730 ( 3.312 -3.312 1.937) 5.069 3.811 ( -14.007 14.007 16.119) 25.539 3.966 ( 5.405 -5.405 2.129) 7.934 4.005 ( -11.614 11.614 22.522) 27.875 4.727 ( -3.010 3.010 -9.128) 10.072 5.010 ( 6.768 -6.768 -0.504) 9.585 5.203 ( -2.397 2.397 2.393) 4.149 5.308 ( -2.401 2.401 0.993) 3.537 5.512 ( -4.171 4.171 15.386) 16.478 5.711 ( -7.356 7.356 15.151) 18.378 6.501 ( -12.231 12.231 -4.232) 17.807 6.974 ( -1.619 1.619 11.986) 12.203 9.089 ( 9.858 -9.858 -19.345) 23.845 9.194 ( 4.380 -4.380 -12.993) 14.394 9.292 ( 8.648 -8.648 -8.281) 14.770 9.758 ( -1.107 1.107 1.563) 2.212 9.901 ( 4.453 -4.453 9.970) 11.793 10.114 ( -7.006 7.006 9.536) 13.751 10.718 ( -13.610 13.610 -2.624) 19.426 11.044 ( 5.244 -5.244 -16.723) 18.294 11.352 ( -0.856 0.856 24.691) 24.721 11.560 ( -0.969 0.969 5.194) 5.372 11.588 ( 0.159 -0.159 6.601) 6.605 11.948 ( 4.089 -4.089 -4.488) 7.320 12.338 ( -9.269 9.269 -0.936) 13.142 12.792 ( -0.336 0.336 19.485) 19.491 13.851 ( -2.493 2.493 9.531) 10.162 14.089 ( -4.444 4.444 16.689) 17.833 14.133 ( -0.706 0.706 16.742) 16.771 15.495 ( 5.011 -5.011 -10.969) 13.059 15.598 ( -6.320 6.320 8.404) 12.268 15.625 ( 4.576 -4.576 -2.947) 7.111 16.288 ( -2.838 2.838 0.246) 4.021 16.298 ( -3.086 3.086 0.356) 4.379 17.207 ( 3.639 -3.639 -5.978) 7.888 17.224 ( 3.758 -3.758 -6.507) 8.402 17.822 ( -5.560 5.560 7.091) 10.589 19.235 ( -0.627 0.627 -5.177) 5.252 20.034 ( -1.065 1.065 -2.667) 3.063 20.405 ( 0.551 -0.551 -2.822) 2.928 ======================= Grid point 15 (10/47) ======================= q-point: ( 0.40 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.403 ( -5.132 5.132 26.907) 27.869 2.473 ( -4.310 4.310 21.496) 22.343 3.664 ( 3.236 -3.236 -0.920) 4.667 3.784 ( 6.563 -6.563 -0.772) 9.313 3.801 ( 8.662 -8.662 2.150) 12.437 4.730 ( -7.236 7.236 5.252) 11.502 4.761 ( -7.479 7.479 5.293) 11.827 4.879 ( 4.058 -4.058 0.046) 5.740 5.346 ( -6.576 6.576 7.134) 11.721 5.477 ( -9.413 9.413 2.801) 13.603 5.760 ( -3.073 3.073 8.329) 9.395 6.146 ( -12.656 12.656 20.074) 26.895 6.771 ( -11.299 11.299 2.116) 16.119 7.115 ( -0.161 0.161 7.008) 7.011 8.503 ( 16.813 -16.813 -22.664) 32.848 8.932 ( 3.509 -3.509 -13.548) 14.428 9.012 ( 6.806 -6.806 -7.101) 11.961 9.772 ( 0.958 -0.958 1.633) 2.121 9.890 ( 5.763 -5.763 9.661) 12.639 10.416 ( -8.463 8.463 7.826) 14.301 10.749 ( 2.976 -2.976 -14.411) 15.013 11.037 ( -14.650 14.650 -2.788) 20.906 11.563 ( 2.776 -2.776 9.824) 10.579 11.648 ( -0.817 0.817 5.268) 5.393 11.685 ( -0.931 0.931 9.282) 9.375 11.779 ( 4.349 -4.349 -5.862) 8.497 12.589 ( -8.404 8.404 6.733) 13.660 12.987 ( 3.465 -3.465 17.951) 18.608 14.036 ( -3.445 3.445 8.198) 9.536 14.320 ( -0.966 0.966 11.413) 11.494 14.395 ( -3.739 3.739 16.325) 17.160 15.226 ( 5.322 -5.322 -11.341) 13.611 15.396 ( 8.116 -8.116 -4.859) 12.464 15.883 ( -5.998 5.998 8.380) 11.924 16.405 ( -5.072 5.072 3.345) 7.914 16.410 ( -5.031 5.031 1.725) 7.322 17.026 ( 4.682 -4.682 -7.228) 9.802 17.027 ( 4.847 -4.847 -8.651) 11.038 18.040 ( -5.445 5.445 6.181) 9.875 19.148 ( 2.290 -2.290 -5.286) 6.199 20.061 ( -3.444 3.444 -2.083) 5.298 20.346 ( 0.907 -0.907 -4.044) 4.243 ======================= Grid point 16 (11/47) ======================= q-point: ( 0.50 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 280 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.700 ( 2.536 -2.536 8.951) 9.643 2.768 ( 0.978 -0.978 22.369) 22.412 3.549 ( 2.436 -2.436 -0.808) 3.538 3.619 ( 6.112 -6.112 2.546) 9.011 3.699 ( 4.804 -4.804 6.764) 9.587 4.786 ( -0.839 0.839 1.037) 1.575 4.808 ( 2.083 -2.083 -0.011) 2.946 4.914 ( -1.899 1.899 -1.052) 2.884 5.565 ( -4.376 4.376 2.099) 6.535 5.723 ( -6.334 6.334 1.530) 9.087 5.815 ( 0.248 -0.248 -4.943) 4.955 6.741 ( -11.771 11.771 14.754) 22.244 7.057 ( -5.728 5.728 8.743) 11.919 7.119 ( 3.443 -3.443 -1.374) 5.059 7.858 ( 12.049 -12.049 -13.745) 21.893 8.752 ( -1.258 1.258 -9.600) 9.764 8.833 ( 2.336 -2.336 -2.540) 4.168 9.769 ( 1.004 -1.004 2.159) 2.584 9.870 ( 5.456 -5.456 9.662) 12.365 10.588 ( -1.866 1.866 -9.219) 9.590 10.641 ( -2.821 2.821 3.850) 5.544 11.317 ( -11.596 11.596 -3.456) 16.760 11.512 ( 3.754 -3.754 -3.990) 6.641 11.624 ( 1.456 -1.456 -4.977) 5.387 11.724 ( 0.101 -0.101 3.994) 3.996 11.764 ( 2.029 -2.029 5.131) 5.878 12.855 ( -2.694 2.694 12.813) 13.367 13.024 ( 8.925 -8.925 12.190) 17.547 14.167 ( 1.057 -1.057 6.929) 7.088 14.472 ( -3.993 3.993 4.521) 7.233 14.562 ( 2.506 -2.506 3.043) 4.671 15.074 ( -0.551 0.551 1.962) 2.111 15.149 ( 7.192 -7.192 -3.444) 10.738 16.048 ( -1.185 1.185 3.530) 3.908 16.551 ( -4.265 4.265 -0.781) 6.081 16.587 ( -4.365 4.365 4.861) 7.857 16.823 ( 2.447 -2.447 -6.503) 7.367 16.833 ( 3.186 -3.186 -6.045) 7.540 18.192 ( -2.022 2.022 2.162) 3.585 19.046 ( 1.195 -1.195 -1.999) 2.618 20.124 ( -3.496 3.496 -3.393) 5.997 20.280 ( 0.157 -0.157 -3.222) 3.229 ======================= Grid point 17 (12/47) ======================= q-point: (-0.40 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.530 ( 11.687 -11.687 -5.417) 17.393 2.796 ( 14.939 -14.939 13.921) 25.301 3.615 ( -3.410 3.410 6.504) 8.097 3.650 ( -2.406 2.406 9.888) 10.457 3.694 ( 2.467 -2.467 6.894) 7.727 4.772 ( 0.667 -0.667 -0.585) 1.111 4.811 ( -1.105 1.105 -0.267) 1.585 4.897 ( 1.866 -1.866 -3.905) 4.713 5.552 ( 6.585 -6.585 -0.245) 9.315 5.652 ( 3.298 -3.298 -13.512) 14.294 5.788 ( -0.016 0.016 -4.239) 4.239 6.563 ( 19.195 -19.195 -13.291) 30.225 6.915 ( 4.518 -4.518 -12.304) 13.864 7.193 ( 3.136 -3.136 7.820) 8.990 8.059 ( -18.559 18.559 17.192) 31.376 8.798 ( -8.153 8.153 -1.967) 11.697 8.860 ( -4.414 4.414 2.191) 6.616 9.780 ( 0.578 -0.578 2.864) 2.979 9.862 ( 4.840 -4.840 9.396) 11.625 10.595 ( 5.972 -5.972 -4.086) 9.382 10.626 ( -7.302 7.302 -3.087) 10.779 11.295 ( 4.597 -4.597 -14.896) 16.253 11.494 ( -3.067 3.067 0.597) 4.379 11.624 ( -6.262 6.262 -2.475) 9.195 11.725 ( 3.241 -3.241 1.643) 4.869 11.732 ( 4.433 -4.433 3.530) 7.195 12.856 ( 12.059 -12.059 3.234) 17.357 12.974 ( 4.910 -4.910 10.656) 12.719 14.135 ( 5.263 -5.263 0.847) 7.490 14.389 ( 6.745 -6.745 -7.343) 12.038 14.550 ( 4.107 -4.107 4.313) 7.234 15.068 ( -4.204 4.204 -0.068) 5.946 15.260 ( -4.386 4.386 12.824) 14.246 16.011 ( 4.594 -4.594 -2.534) 6.974 16.476 ( 2.459 -2.459 -9.309) 9.937 16.548 ( 4.815 -4.815 -4.580) 8.207 16.880 ( -6.288 6.288 2.641) 9.276 16.906 ( -6.442 6.442 5.445) 10.613 18.157 ( 3.840 -3.840 -3.587) 6.508 19.080 ( -2.972 2.972 3.170) 5.264 20.102 ( -0.432 0.432 -6.584) 6.612 20.302 ( -2.808 2.808 0.098) 3.972 ======================= Grid point 18 (13/47) ======================= q-point: (-0.30 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 280 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.097 ( 15.647 -15.647 -10.199) 24.366 2.473 ( 22.443 -22.443 9.806) 33.219 3.708 ( 1.363 -1.363 6.197) 6.489 3.810 ( -4.005 4.005 4.877) 7.474 3.858 ( -2.950 2.950 8.989) 9.910 4.550 ( 21.042 -21.042 -16.253) 33.908 4.761 ( -1.201 1.201 -1.872) 2.528 4.848 ( -3.526 3.526 -1.904) 5.338 5.345 ( 7.625 -7.625 2.757) 11.130 5.386 ( 4.266 -4.266 -13.376) 14.674 5.612 ( 7.273 -7.273 -7.361) 12.648 6.043 ( 15.710 -15.710 -0.011) 22.217 6.649 ( 1.837 -1.837 -17.369) 17.562 7.117 ( 8.590 -8.590 4.545) 12.971 8.727 ( -19.292 19.292 13.203) 30.310 9.038 ( -10.112 10.112 2.523) 14.522 9.064 ( -8.721 8.721 4.321) 13.068 9.809 ( 0.343 -0.343 3.609) 3.641 9.855 ( 4.574 -4.574 7.624) 9.999 10.323 ( 9.423 -9.423 -7.870) 15.476 10.844 ( -11.023 11.023 0.988) 15.620 10.971 ( 4.625 -4.625 -17.958) 19.113 11.459 ( 7.834 -7.834 1.879) 11.237 11.613 ( 4.928 -4.928 -0.233) 6.974 11.652 ( 5.145 -5.145 2.681) 7.754 11.808 ( -7.607 7.607 1.335) 10.840 12.579 ( 9.366 -9.366 -4.427) 13.966 12.886 ( 9.871 -9.871 5.528) 15.015 13.962 ( 6.250 -6.250 -5.172) 10.240 14.126 ( 7.810 -7.810 -6.709) 12.923 14.404 ( 12.769 -12.769 7.863) 19.695 15.287 ( -10.858 10.858 1.383) 15.417 15.512 ( -4.133 4.133 10.744) 12.231 15.793 ( 8.345 -8.345 -6.330) 13.392 16.328 ( 1.135 -1.135 -6.868) 7.053 16.389 ( 3.669 -3.669 -3.724) 6.387 17.072 ( -6.588 6.588 2.354) 9.609 17.125 ( -5.670 5.670 5.473) 9.708 17.972 ( 6.405 -6.405 -4.944) 10.319 19.205 ( -3.241 3.241 2.968) 5.461 20.040 ( -0.833 0.833 -4.912) 5.052 20.372 ( -2.497 2.497 0.320) 3.545 ======================= Grid point 19 (14/47) ======================= q-point: (-0.20 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.566 ( 18.440 -18.440 -9.439) 27.734 1.988 ( 26.321 -26.321 7.209) 37.915 3.512 ( 35.917 -35.917 -20.087) 54.622 3.751 ( 0.824 -0.824 5.369) 5.494 3.890 ( -0.534 0.534 -1.808) 1.960 4.023 ( 3.988 -3.988 4.974) 7.520 4.812 ( -5.527 5.527 -1.692) 7.998 4.952 ( -6.820 6.820 -1.191) 9.718 5.177 ( 2.957 -2.957 -6.105) 7.399 5.268 ( 3.719 -3.719 2.627) 5.879 5.368 ( 5.801 -5.801 -4.331) 9.277 5.715 ( 15.479 -15.479 5.731) 22.628 6.441 ( -0.081 0.081 -13.456) 13.456 6.918 ( 11.027 -11.027 2.067) 15.731 9.226 ( -13.201 13.201 3.288) 18.956 9.293 ( -8.248 8.248 2.379) 11.905 9.368 ( -11.049 11.049 6.608) 16.965 9.806 ( 4.966 -4.966 2.648) 7.505 9.832 ( 1.972 -1.972 3.597) 4.551 10.004 ( 8.981 -8.981 -6.966) 14.486 10.681 ( 2.914 -2.914 -11.634) 12.342 11.162 ( -12.314 12.314 2.268) 17.561 11.175 ( 15.069 -15.069 1.256) 21.348 11.529 ( 2.832 -2.832 0.044) 4.005 11.561 ( 4.356 -4.356 1.652) 6.378 11.985 ( -5.001 5.001 2.527) 7.510 12.335 ( 5.172 -5.172 -6.564) 9.828 12.667 ( 11.377 -11.377 1.730) 16.183 13.751 ( 5.383 -5.383 -5.232) 9.237 13.872 ( 7.720 -7.720 -3.966) 11.616 14.155 ( 14.916 -14.916 6.516) 22.078 15.506 ( 8.904 -8.904 -5.282) 13.655 15.553 ( -9.373 9.373 1.133) 13.304 15.698 ( -3.356 3.356 5.089) 6.959 16.265 ( -0.364 0.364 -2.510) 2.562 16.292 ( 1.694 -1.694 -1.773) 2.981 17.237 ( -5.307 5.307 1.092) 7.584 17.293 ( -3.921 3.921 3.129) 6.367 17.761 ( 6.706 -6.706 -3.296) 10.041 19.278 ( -1.339 1.339 0.097) 1.896 20.036 ( -1.951 1.951 -1.604) 3.192 20.423 ( -1.487 1.487 0.157) 2.108 ======================= Grid point 20 (15/47) ======================= q-point: (-0.10 0.10 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.024 ( 20.684 -20.684 0.000) 29.252 1.383 ( 28.151 -28.151 0.000) 39.812 2.357 ( 46.195 -46.195 0.000) 65.329 3.781 ( 0.384 -0.384 0.000) 0.542 3.833 ( 1.906 -1.906 0.000) 2.695 3.946 ( 4.952 -4.952 0.000) 7.003 4.991 ( -9.602 9.602 0.000) 13.579 5.089 ( 1.779 -1.779 0.000) 2.516 5.120 ( -5.694 5.694 0.000) 8.053 5.176 ( 3.002 -3.002 0.000) 4.246 5.249 ( 3.746 -3.746 0.000) 5.297 5.459 ( 8.163 -8.163 0.000) 11.544 6.336 ( 1.831 -1.831 0.000) 2.590 6.651 ( 11.922 -11.922 0.000) 16.861 9.480 ( -5.999 5.999 0.000) 8.484 9.481 ( -7.198 7.198 0.000) 10.180 9.673 ( -9.695 9.695 0.000) 13.711 9.698 ( 4.140 -4.140 0.000) 5.855 9.749 ( 7.130 -7.130 0.000) 10.084 9.785 ( 4.919 -4.919 0.000) 6.957 10.572 ( 0.692 -0.692 0.000) 0.978 10.852 ( 11.665 -11.665 0.000) 16.497 11.429 ( -5.179 5.179 0.000) 7.324 11.466 ( 2.294 -2.294 0.000) 3.245 11.513 ( -2.490 2.490 0.000) 3.521 12.097 ( -2.547 2.547 0.000) 3.602 12.194 ( 2.556 -2.556 0.000) 3.614 12.414 ( 9.987 -9.987 0.000) 14.123 13.607 ( 3.564 -3.564 0.000) 5.040 13.679 ( 6.044 -6.044 0.000) 8.547 13.859 ( 12.725 -12.725 0.000) 17.996 15.277 ( 7.099 -7.099 0.000) 10.040 15.744 ( -5.847 5.847 0.000) 8.269 15.798 ( -2.666 2.666 0.000) 3.770 16.262 ( 0.091 -0.091 0.000) 0.129 16.267 ( -0.553 0.553 0.000) 0.782 17.351 ( -3.623 3.623 0.000) 5.123 17.386 ( -2.314 2.314 0.000) 3.272 17.589 ( 5.662 -5.662 0.000) 8.007 19.291 ( -0.287 0.287 0.000) 0.405 20.078 ( -1.843 1.843 0.000) 2.606 20.449 ( -0.657 0.657 0.000) 0.929 ======================= Grid point 23 (16/47) ======================= q-point: ( 0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 102 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.879 ( 0.000 -0.000 35.089) 35.089 1.879 ( 0.000 -0.000 35.089) 35.089 3.061 ( 0.000 -0.000 59.422) 59.422 3.833 ( 0.000 -0.000 0.064) 0.064 3.833 ( 0.000 -0.000 0.064) 0.064 4.099 ( 0.000 -0.000 12.590) 12.590 4.693 ( -0.000 0.000 -10.637) 10.637 5.235 ( 0.000 -0.000 0.693) 0.693 5.235 ( 0.000 -0.000 0.693) 0.693 5.280 ( 0.000 -0.000 6.373) 6.373 5.457 ( 0.000 -0.000 13.049) 13.049 5.457 ( 0.000 -0.000 13.049) 13.049 6.077 ( -0.000 0.000 -7.033) 7.033 6.913 ( 0.000 -0.000 11.413) 11.413 9.311 ( -0.000 0.000 -15.458) 15.458 9.311 ( -0.000 0.000 -15.458) 15.458 9.543 ( -0.000 0.000 -12.798) 12.798 9.640 ( 0.000 -0.000 2.725) 2.725 10.036 ( 0.000 -0.000 16.438) 16.438 10.036 ( 0.000 -0.000 16.438) 16.438 10.329 ( -0.000 0.000 -7.014) 7.014 11.194 ( -0.000 0.000 -21.222) 21.222 11.390 ( 0.000 -0.000 30.452) 30.452 11.538 ( 0.000 -0.000 5.721) 5.721 11.538 ( 0.000 -0.000 5.721) 5.721 12.066 ( -0.000 0.000 -3.788) 3.788 12.066 ( -0.000 0.000 -3.788) 3.788 12.757 ( 0.000 -0.000 21.426) 21.426 13.975 ( 0.000 -0.000 18.680) 18.680 13.975 ( 0.000 -0.000 18.680) 18.680 14.041 ( 0.000 -0.000 17.657) 17.657 15.381 ( 0.000 -0.000 10.759) 10.759 15.632 ( -0.000 0.000 -9.845) 9.845 15.754 ( -0.000 0.000 -3.873) 3.873 16.215 ( -0.000 0.000 -1.961) 1.961 16.215 ( -0.000 0.000 -1.961) 1.961 17.317 ( -0.000 0.000 -6.549) 6.549 17.317 ( -0.000 0.000 -6.549) 6.549 17.655 ( 0.000 -0.000 8.261) 8.261 19.139 ( -0.000 0.000 -6.222) 6.222 20.034 ( -0.000 0.000 -3.387) 3.387 20.417 ( -0.000 0.000 -2.842) 2.842 ======================= Grid point 24 (17/47) ======================= q-point: ( 0.30 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.282 ( -0.385 0.385 33.701) 33.705 2.305 ( -2.067 2.067 30.580) 30.719 3.793 ( -4.436 4.436 17.214) 18.321 3.802 ( 1.863 -1.863 -0.957) 2.803 4.004 ( -12.608 12.608 30.069) 34.958 4.193 ( 5.541 -5.541 12.026) 14.353 4.594 ( 0.047 -0.047 -4.873) 4.873 5.146 ( 5.526 -5.526 -2.646) 8.250 5.245 ( -1.860 1.860 -1.327) 2.947 5.301 ( 3.085 -3.085 6.961) 8.215 5.638 ( -0.195 0.195 15.944) 15.947 5.745 ( -5.832 5.832 18.359) 20.127 6.160 ( -11.905 11.905 -3.598) 17.216 7.040 ( -1.764 1.764 6.203) 6.685 9.034 ( 5.416 -5.416 -20.524) 21.907 9.076 ( 1.889 -1.889 -19.066) 19.252 9.335 ( 4.528 -4.528 -11.349) 13.031 9.726 ( -3.399 3.399 2.996) 5.664 10.149 ( 1.140 -1.140 6.608) 6.802 10.184 ( 4.626 -4.626 16.695) 17.931 10.437 ( -9.618 9.618 3.802) 14.123 10.892 ( 2.190 -2.190 -23.188) 23.394 11.610 ( -0.177 0.177 11.998) 12.000 11.618 ( -0.371 0.371 6.578) 6.599 11.732 ( 2.820 -2.820 19.080) 19.493 11.976 ( 2.829 -2.829 -5.437) 6.750 12.137 ( -8.604 8.604 3.743) 12.730 13.034 ( -0.533 0.533 22.748) 22.761 14.082 ( 5.933 -5.933 18.684) 20.482 14.247 ( -1.855 1.855 22.207) 22.362 14.325 ( -1.603 1.603 17.879) 18.022 15.459 ( 3.263 -3.263 -12.291) 13.128 15.577 ( -4.674 4.674 11.321) 13.110 15.677 ( 2.925 -2.925 -3.164) 5.207 16.221 ( -2.266 2.266 -0.926) 3.335 16.222 ( -2.363 2.363 -1.376) 3.614 17.192 ( 2.142 -2.142 -8.969) 9.467 17.197 ( 1.924 -1.924 -9.641) 10.018 17.806 ( -3.401 3.401 8.843) 10.066 19.101 ( -3.087 3.087 -7.584) 8.750 19.987 ( 0.201 -0.201 -3.459) 3.470 20.370 ( 0.213 -0.213 -4.472) 4.482 ======================= Grid point 25 (18/47) ======================= q-point: ( 0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 282 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.683 ( -1.592 1.592 24.183) 24.287 2.708 ( -1.761 1.761 32.316) 32.412 3.722 ( 3.118 -3.118 -2.490) 5.064 3.893 ( 0.537 -0.537 1.830) 1.981 4.120 ( 10.923 -10.923 10.437) 18.643 4.602 ( -4.746 4.746 -0.600) 6.739 4.680 ( -7.860 7.860 17.602) 20.818 4.985 ( 5.102 -5.102 -1.226) 7.319 5.297 ( -2.075 2.075 0.141) 2.938 5.414 ( -4.897 4.897 15.615) 17.082 5.773 ( -1.172 1.172 1.997) 2.596 6.139 ( -8.150 8.150 20.687) 23.681 6.531 ( -14.365 14.365 10.521) 22.878 7.138 ( -2.231 2.231 1.309) 3.416 8.558 ( 11.467 -11.467 -24.623) 29.484 8.789 ( 1.733 -1.733 -20.456) 20.602 9.089 ( 6.219 -6.219 -6.826) 11.133 9.812 ( -0.360 0.360 2.587) 2.637 10.220 ( 8.003 -8.003 15.245) 18.986 10.285 ( -8.201 8.201 3.444) 12.098 10.578 ( 0.980 -0.980 -20.921) 20.966 10.744 ( -9.765 9.765 6.157) 15.120 11.707 ( -0.379 0.379 6.322) 6.345 11.733 ( 4.769 -4.769 1.285) 6.865 11.767 ( -1.603 1.603 8.325) 8.628 11.808 ( 4.431 -4.431 -7.235) 9.571 12.538 ( -9.057 9.057 21.550) 25.069 13.294 ( 1.909 -1.909 21.857) 22.023 14.175 ( 3.076 -3.076 15.013) 15.630 14.496 ( 2.636 -2.636 12.839) 13.370 14.563 ( -1.391 1.391 21.922) 22.010 15.205 ( 4.519 -4.519 -12.503) 14.041 15.515 ( 6.605 -6.605 -3.543) 9.991 15.854 ( -5.904 5.904 10.530) 13.438 16.296 ( -4.682 4.682 -0.638) 6.652 16.319 ( -5.156 5.156 2.778) 7.803 16.999 ( 3.280 -3.280 -12.367) 13.209 17.014 ( 2.802 -2.802 -10.504) 11.226 18.006 ( -4.290 4.290 8.361) 10.330 19.072 ( -1.997 1.997 -8.146) 8.622 19.961 ( -1.619 1.619 -3.009) 3.781 20.300 ( 0.359 -0.359 -5.907) 5.929 ======================= Grid point 26 (19/47) ======================= q-point: ( 0.50 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.914 ( 3.126 -3.126 13.948) 14.632 3.116 ( -0.269 0.269 30.402) 30.404 3.609 ( 3.131 -3.131 -1.079) 4.557 3.833 ( 6.856 -6.856 4.450) 10.668 3.979 ( 9.900 -9.900 9.759) 17.066 4.728 ( -4.364 4.364 -1.247) 6.297 4.846 ( -2.580 2.580 0.034) 3.649 4.882 ( 3.567 -3.567 1.122) 5.168 5.363 ( -6.197 6.197 -4.815) 10.000 5.738 ( 0.351 -0.351 -9.070) 9.084 5.764 ( -2.303 2.303 17.909) 18.203 6.650 ( -12.073 12.073 19.523) 25.936 7.016 ( -10.438 10.438 17.108) 22.596 7.150 ( 0.773 -0.773 -3.784) 3.939 7.951 ( 12.779 -12.779 -21.975) 28.451 8.544 ( -1.845 1.845 -18.415) 18.599 8.900 ( 5.255 -5.255 -1.101) 7.513 9.826 ( 1.099 -1.099 2.380) 2.843 10.179 ( 6.759 -6.759 8.595) 12.854 10.391 ( -1.232 1.232 -9.781) 9.935 10.567 ( -9.358 9.358 4.401) 13.946 11.017 ( -8.730 8.730 1.774) 12.472 11.555 ( 5.908 -5.908 -8.617) 12.002 11.610 ( 3.332 -3.332 -7.621) 8.960 11.785 ( 0.084 -0.084 5.145) 5.146 11.876 ( 0.505 -0.505 6.508) 6.547 12.986 ( -4.401 4.401 23.122) 23.945 13.421 ( 8.175 -8.175 17.537) 21.005 14.293 ( -0.624 0.624 9.532) 9.573 14.520 ( 4.564 -4.564 6.291) 9.013 14.688 ( 5.561 -5.561 0.069) 7.864 15.103 ( -1.802 1.802 8.886) 9.245 15.293 ( 7.916 -7.916 -2.855) 11.553 16.066 ( -2.405 2.405 6.385) 7.234 16.410 ( -5.324 5.324 -2.695) 7.998 16.521 ( -6.613 6.613 6.075) 11.152 16.774 ( 2.211 -2.211 -11.305) 11.729 16.816 ( 3.004 -3.004 -10.262) 11.106 18.181 ( -2.385 2.385 5.517) 6.466 19.014 ( -0.924 0.924 -5.355) 5.512 19.987 ( -3.276 3.276 -3.829) 6.011 20.225 ( -0.330 0.330 -5.448) 5.467 ======================= Grid point 27 (20/47) ======================= q-point: (-0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 282 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.802 ( 12.553 -12.553 0.576) 17.762 3.357 ( 11.085 -11.085 23.535) 28.278 3.626 ( -1.914 1.914 5.608) 6.227 3.763 ( 1.931 -1.931 8.848) 9.260 3.893 ( 6.245 -6.245 7.723) 11.732 4.722 ( -1.527 1.527 -7.851) 8.142 4.828 ( 3.148 -3.148 2.518) 5.114 4.901 ( 0.108 -0.108 0.309) 0.345 5.479 ( -8.730 8.730 -6.987) 14.186 5.558 ( 2.571 -2.571 -13.779) 14.251 5.901 ( 5.949 -5.949 13.530) 15.933 6.844 ( 16.220 -16.220 -11.186) 25.521 6.884 ( 2.728 -2.728 -9.505) 10.258 7.344 ( -0.526 0.526 12.602) 12.624 7.837 ( -13.658 13.658 8.723) 21.194 8.504 ( -10.198 10.198 -10.098) 17.606 8.824 ( 0.894 -0.894 1.195) 1.740 9.831 ( 0.754 -0.754 2.443) 2.666 10.146 ( 4.779 -4.779 9.158) 11.382 10.388 ( -4.988 4.988 -3.320) 7.796 10.703 ( 3.355 -3.355 1.513) 4.980 11.243 ( -12.187 12.187 -2.498) 17.414 11.257 ( 6.908 -6.908 -14.584) 17.554 11.520 ( -1.434 1.434 -3.326) 3.896 11.817 ( 1.784 -1.784 3.207) 4.081 11.893 ( 3.739 -3.739 4.894) 7.204 13.235 ( 3.696 -3.696 15.332) 16.199 13.292 ( 15.894 -15.894 9.181) 24.281 14.348 ( 8.579 -8.579 6.195) 13.623 14.471 ( 6.373 -6.373 -6.558) 11.146 14.529 ( -6.744 6.744 0.191) 9.540 15.112 ( 5.840 -5.840 1.405) 8.378 15.345 ( -2.735 2.735 16.280) 16.733 16.093 ( 3.385 -3.385 -0.113) 4.788 16.375 ( 0.129 -0.129 -12.512) 12.513 16.597 ( 3.924 -3.924 -5.340) 7.702 16.749 ( -5.729 5.729 -0.328) 8.108 16.797 ( -7.499 7.499 4.212) 11.412 18.209 ( 2.754 -2.754 -0.406) 3.916 19.039 ( -3.375 3.375 1.237) 4.931 19.983 ( -1.695 1.695 -6.968) 7.369 20.222 ( -3.050 3.050 -1.535) 4.578 ======================= Grid point 28 (21/47) ======================= q-point: (-0.30 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.419 ( 16.884 -16.884 -2.854) 24.048 3.140 ( 21.120 -21.120 10.732) 31.738 3.761 ( -3.287 3.287 4.049) 6.164 3.855 ( 3.322 -3.322 4.709) 6.652 3.927 ( -1.684 1.684 10.203) 10.477 4.676 ( -3.357 3.357 -7.872) 9.193 4.777 ( 3.032 -3.032 1.070) 4.419 4.785 ( 8.370 -8.370 -4.687) 12.731 5.333 ( 3.212 -3.212 -9.878) 10.872 5.545 ( -0.194 0.194 -2.252) 2.268 5.798 ( 15.041 -15.041 7.306) 22.491 6.310 ( 8.069 -8.069 -8.390) 14.164 6.616 ( 10.892 -10.892 -8.722) 17.701 7.351 ( 6.717 -6.717 4.014) 10.312 8.397 ( -19.517 19.517 8.567) 28.900 8.754 ( -14.131 14.131 -1.256) 20.023 8.898 ( -6.202 6.202 0.692) 8.798 9.850 ( -0.017 0.017 2.162) 2.162 10.125 ( 7.411 -7.411 9.719) 14.294 10.497 ( 10.608 -10.608 -2.935) 15.286 10.560 ( -11.294 11.294 -1.234) 16.020 10.915 ( 6.550 -6.550 -13.256) 16.172 11.517 ( -5.240 5.240 1.349) 7.533 11.562 ( -8.606 8.606 -2.351) 12.395 11.761 ( 5.180 -5.180 1.002) 7.394 11.823 ( 5.675 -5.675 2.645) 8.451 12.926 ( 18.591 -18.591 1.864) 26.357 13.188 ( 9.525 -9.525 5.408) 14.516 14.127 ( 10.352 -10.352 -0.789) 14.662 14.251 ( 7.232 -7.232 -3.290) 10.743 14.697 ( 0.414 -0.414 5.036) 5.070 15.081 ( -3.707 3.707 1.286) 5.398 15.565 ( -1.853 1.853 9.947) 10.286 15.927 ( 7.654 -7.654 -3.515) 11.381 16.239 ( -0.658 0.658 -8.352) 8.404 16.441 ( 4.192 -4.192 -3.027) 6.657 16.924 ( -7.365 7.365 0.917) 10.456 17.028 ( -7.049 7.049 3.395) 10.531 18.077 ( 5.929 -5.929 -2.210) 8.672 19.162 ( -3.676 3.676 2.646) 5.833 19.942 ( -1.690 1.690 -4.878) 5.432 20.300 ( -3.582 3.582 -0.143) 5.068 ======================= Grid point 29 (22/47) ======================= q-point: (-0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 171 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.967 ( 18.683 -18.683 0.000) 26.422 2.656 ( 24.519 -24.519 0.000) 34.675 3.826 ( 1.744 -1.744 0.000) 2.466 3.852 ( -1.948 1.948 0.000) 2.755 3.955 ( 6.500 -6.500 0.000) 9.192 4.352 ( 21.758 -21.758 0.000) 30.771 4.733 ( -0.002 0.002 0.000) 0.002 4.759 ( -6.834 6.834 0.000) 9.665 5.203 ( 2.754 -2.754 0.000) 3.895 5.461 ( 5.733 -5.733 0.000) 8.108 5.470 ( 10.399 -10.399 0.000) 14.707 6.176 ( 17.620 -17.620 0.000) 24.919 6.299 ( -3.077 3.077 0.000) 4.351 7.176 ( 9.346 -9.346 0.000) 13.217 8.890 ( -16.950 16.950 0.000) 23.971 9.068 ( -11.721 11.721 0.000) 16.576 9.122 ( -11.325 11.325 0.000) 16.015 9.873 ( -0.311 0.311 0.000) 0.440 9.997 ( 7.957 -7.957 0.000) 11.254 10.205 ( 11.206 -11.206 0.000) 15.847 10.691 ( 4.529 -4.529 0.000) 6.404 10.860 ( -13.315 13.315 0.000) 18.830 11.482 ( 7.460 -7.460 0.000) 10.550 11.614 ( 5.390 -5.390 0.000) 7.623 11.697 ( 5.791 -5.791 0.000) 8.190 11.820 ( -9.687 9.687 0.000) 13.699 12.529 ( 13.742 -13.742 0.000) 19.435 12.958 ( 11.544 -11.544 0.000) 16.325 13.889 ( 8.379 -8.379 0.000) 11.850 14.044 ( 8.272 -8.272 0.000) 11.698 14.546 ( 12.895 -12.895 0.000) 18.236 15.303 ( -11.404 11.404 0.000) 16.128 15.667 ( -2.176 2.176 0.000) 3.077 15.693 ( 9.245 -9.245 0.000) 13.075 16.221 ( -1.395 1.395 0.000) 1.973 16.335 ( 3.041 -3.041 0.000) 4.301 17.100 ( -6.659 6.659 0.000) 9.418 17.200 ( -5.533 5.533 0.000) 7.825 17.905 ( 6.863 -6.863 0.000) 9.706 19.245 ( -1.981 1.981 0.000) 2.801 19.966 ( -2.502 2.502 0.000) 3.539 20.375 ( -2.628 2.628 0.000) 3.717 ======================= Grid point 34 (23/47) ======================= q-point: ( 0.30 0.30 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 99 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.637 ( 0.000 -0.000 27.889) 27.889 2.637 ( 0.000 -0.000 27.889) 27.889 3.802 ( -0.000 0.000 -2.457) 2.457 3.802 ( -0.000 0.000 -2.457) 2.457 4.372 ( 0.000 -0.000 48.478) 48.478 4.444 ( 0.000 -0.000 15.823) 15.823 4.587 ( 0.000 -0.000 1.884) 1.884 5.179 ( -0.000 0.000 -3.970) 3.970 5.179 ( -0.000 0.000 -3.970) 3.970 5.434 ( 0.000 -0.000 6.662) 6.662 5.867 ( -0.000 0.000 -10.409) 10.409 5.930 ( 0.000 -0.000 24.500) 24.500 5.930 ( 0.000 -0.000 24.500) 24.500 7.051 ( -0.000 0.000 -0.173) 0.173 8.840 ( -0.000 0.000 -23.411) 23.411 8.840 ( -0.000 0.000 -23.411) 23.411 9.255 ( -0.000 0.000 -9.786) 9.786 9.717 ( 0.000 -0.000 3.628) 3.628 10.170 ( -0.000 0.000 -6.156) 6.156 10.461 ( 0.000 -0.000 18.392) 18.392 10.461 ( 0.000 -0.000 18.392) 18.392 10.611 ( -0.000 0.000 -26.746) 26.746 11.693 ( 0.000 -0.000 6.565) 6.565 11.693 ( 0.000 -0.000 6.565) 6.565 11.937 ( -0.000 0.000 -7.030) 7.030 11.937 ( -0.000 0.000 -7.030) 7.030 12.201 ( 0.000 -0.000 36.142) 36.142 13.311 ( 0.000 -0.000 24.249) 24.249 14.397 ( 0.000 -0.000 8.636) 8.636 14.520 ( 0.000 -0.000 26.361) 26.361 14.520 ( 0.000 -0.000 26.361) 26.361 15.340 ( -0.000 0.000 -14.492) 14.492 15.658 ( 0.000 -0.000 11.692) 11.692 15.702 ( 0.000 -0.000 1.155) 1.155 16.172 ( -0.000 0.000 -1.627) 1.627 16.172 ( -0.000 0.000 -1.627) 1.627 17.083 ( -0.000 0.000 -13.211) 13.211 17.083 ( -0.000 0.000 -13.211) 13.211 17.882 ( 0.000 -0.000 10.533) 10.533 18.960 ( -0.000 0.000 -8.650) 8.650 19.951 ( -0.000 0.000 -3.147) 3.147 20.309 ( -0.000 0.000 -6.116) 6.116 ======================= Grid point 35 (24/47) ======================= q-point: ( 0.40 0.30 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.948 ( -0.193 0.193 22.700) 22.702 2.995 ( -3.525 3.525 25.356) 25.842 3.746 ( 1.680 -1.680 -3.310) 4.074 3.826 ( -4.455 4.455 -2.184) 6.668 4.523 ( 9.048 -9.048 14.488) 19.330 4.617 ( 1.837 -1.837 6.788) 7.268 4.828 ( 2.173 -2.173 18.756) 19.006 5.074 ( 3.763 -3.763 -2.599) 5.923 5.216 ( -4.466 4.466 0.378) 6.327 5.555 ( -2.383 2.383 18.638) 18.940 5.728 ( 0.031 -0.031 -8.585) 8.586 6.275 ( -3.043 3.043 25.117) 25.483 6.399 ( -11.602 11.602 23.079) 28.317 7.048 ( -2.682 2.682 -5.633) 6.791 8.467 ( 6.266 -6.266 -26.162) 27.622 8.537 ( 0.456 -0.456 -25.294) 25.303 9.133 ( 2.764 -2.764 -3.757) 5.422 9.806 ( -2.726 2.726 3.578) 5.260 10.132 ( -3.518 3.518 -3.981) 6.372 10.278 ( 0.450 -0.450 -26.142) 26.150 10.616 ( 5.057 -5.057 17.525) 18.928 10.728 ( -3.804 3.804 15.938) 16.822 11.771 ( -0.328 0.328 5.643) 5.662 11.791 ( -0.347 0.347 1.910) 1.972 11.803 ( 1.598 -1.598 -3.782) 4.406 11.805 ( 3.032 -3.032 -8.836) 9.821 12.702 ( -5.016 5.016 33.690) 34.429 13.598 ( 0.638 -0.638 23.839) 23.856 14.392 ( 3.213 -3.213 0.701) 4.597 14.788 ( 4.654 -4.654 23.708) 24.604 14.843 ( 0.771 -0.771 25.099) 25.123 15.127 ( 1.897 -1.897 -13.733) 13.993 15.661 ( 6.283 -6.283 2.694) 9.285 15.853 ( -4.614 4.614 10.579) 12.430 16.188 ( -2.843 2.843 -1.607) 4.330 16.222 ( -4.971 4.971 1.209) 7.133 16.886 ( 1.510 -1.510 -16.142) 16.283 16.898 ( 0.554 -0.554 -15.448) 15.468 18.043 ( -2.289 2.289 10.907) 11.377 18.890 ( -3.040 3.040 -9.959) 10.848 19.917 ( 0.052 -0.052 -1.838) 1.839 20.226 ( 0.142 -0.142 -7.313) 7.316 ======================= Grid point 36 (25/47) ======================= q-point: ( 0.50 0.30 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 280 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.157 ( 2.558 -2.558 14.552) 14.995 3.395 ( -4.594 4.594 23.816) 24.687 3.656 ( 2.215 -2.215 -2.387) 3.938 3.930 ( -3.890 3.890 0.262) 5.507 4.377 ( 13.306 -13.306 8.763) 20.758 4.656 ( 1.273 -1.273 4.240) 4.606 4.840 ( 1.627 -1.627 1.451) 2.720 4.998 ( 3.282 -3.282 3.806) 6.003 5.225 ( -0.521 0.521 -6.475) 6.517 5.615 ( -0.163 0.163 -11.825) 11.828 5.975 ( -3.601 3.601 26.876) 27.355 6.669 ( -5.936 5.936 21.231) 22.831 6.983 ( -12.459 12.459 23.955) 29.737 7.039 ( -4.021 4.021 -8.481) 10.211 7.940 ( 10.494 -10.494 -25.863) 29.819 8.240 ( -1.405 1.405 -24.763) 24.842 9.046 ( 4.964 -4.964 3.973) 8.066 9.869 ( 0.287 -0.287 1.896) 1.939 10.006 ( -3.129 3.129 -22.202) 22.639 10.298 ( -12.316 12.316 -0.980) 17.445 10.659 ( 8.192 -8.192 16.902) 20.491 10.989 ( -3.050 3.050 12.710) 13.422 11.579 ( 5.743 -5.743 -10.911) 13.602 11.594 ( 4.299 -4.299 -10.416) 12.061 11.839 ( -0.330 0.330 4.094) 4.120 11.919 ( -2.377 2.377 4.508) 5.623 13.199 ( -3.938 3.938 29.228) 29.754 13.810 ( 5.066 -5.066 20.598) 21.808 14.299 ( 0.143 -0.143 -7.619) 7.622 14.763 ( 5.923 -5.923 -12.641) 15.164 14.861 ( 10.094 -10.094 18.236) 23.159 15.245 ( -1.606 1.606 22.368) 22.483 15.477 ( 9.163 -9.163 -0.014) 12.958 16.079 ( -3.734 3.734 7.770) 9.395 16.248 ( -4.208 4.208 -4.612) 7.529 16.416 ( -7.434 7.434 4.479) 11.427 16.656 ( 0.920 -0.920 -14.482) 14.540 16.710 ( 0.039 -0.039 -13.878) 13.878 18.223 ( -1.568 1.568 9.430) 9.687 18.858 ( -4.343 4.343 -9.303) 11.148 19.911 ( -0.635 0.635 -0.767) 1.181 20.140 ( -0.503 0.503 -6.804) 6.842 ======================= Grid point 37 (26/47) ======================= q-point: (-0.40 0.30 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.107 ( 10.426 -10.426 2.934) 15.033 3.641 ( -1.116 1.116 2.334) 2.818 3.748 ( -2.513 2.513 18.849) 19.182 3.943 ( 5.619 -5.619 3.265) 8.591 4.162 ( 10.260 -10.260 4.100) 15.079 4.595 ( 0.696 -0.696 -7.938) 7.999 4.914 ( -1.467 1.467 5.779) 6.140 4.977 ( 4.356 -4.356 6.375) 8.864 5.177 ( -6.088 6.088 -9.131) 12.550 5.452 ( 1.446 -1.446 -11.480) 11.661 6.259 ( 4.398 -4.398 22.805) 23.638 6.775 ( 1.752 -1.752 -15.628) 15.823 7.104 ( 8.023 -8.023 4.775) 12.310 7.480 ( -3.811 3.811 19.218) 19.959 7.570 ( -9.014 9.014 -7.777) 14.932 8.081 ( -10.301 10.301 -18.789) 23.775 8.954 ( 6.917 -6.917 3.713) 10.464 9.869 ( 0.725 -0.725 1.005) 1.436 9.956 ( -9.446 9.446 -10.393) 16.926 10.623 ( -12.847 12.847 0.748) 18.184 10.634 ( 9.251 -9.251 13.409) 18.734 11.119 ( 2.661 -2.661 5.283) 6.486 11.274 ( 7.734 -7.734 -11.672) 15.996 11.426 ( -1.521 1.521 -6.643) 6.982 11.877 ( 0.465 -0.465 2.098) 2.199 11.991 ( 0.306 -0.306 2.611) 2.647 13.509 ( 3.082 -3.082 16.558) 17.122 13.777 ( 14.202 -14.202 10.303) 22.573 14.281 ( -7.980 7.980 -6.882) 13.218 14.509 ( 3.668 -3.668 -7.915) 9.463 14.724 ( 13.902 -13.902 6.931) 20.847 15.290 ( 6.713 -6.713 1.934) 9.689 15.522 ( -0.181 0.181 18.124) 18.125 16.165 ( 1.731 -1.731 1.626) 2.938 16.199 ( -1.310 1.310 -13.373) 13.501 16.584 ( 0.270 -0.270 -6.513) 6.524 16.601 ( -5.046 5.046 -2.445) 7.544 16.652 ( -8.010 8.010 3.270) 11.791 18.288 ( 2.684 -2.684 2.830) 4.735 18.915 ( -6.568 6.568 -2.416) 9.597 19.901 ( 0.104 -0.104 -2.542) 2.547 20.124 ( -3.085 3.085 -2.376) 4.967 ======================= Grid point 39 (27/47) ======================= q-point: (-0.30 0.30 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.798 ( 15.377 -15.377 0.000) 21.746 3.672 ( 16.798 -16.798 0.000) 23.756 3.712 ( -2.778 2.778 0.000) 3.928 3.956 ( -2.615 2.615 0.000) 3.697 3.994 ( 6.067 -6.067 0.000) 8.580 4.541 ( -1.934 1.934 0.000) 2.735 4.901 ( 5.332 -5.332 0.000) 7.541 4.903 ( 3.686 -3.686 0.000) 5.213 5.329 ( 2.202 -2.202 0.000) 3.114 5.367 ( -10.823 10.823 0.000) 15.306 6.335 ( 6.508 -6.508 0.000) 9.204 6.345 ( 16.776 -16.776 0.000) 23.725 6.838 ( 10.279 -10.279 0.000) 14.536 7.528 ( 5.412 -5.412 0.000) 7.654 7.978 ( -19.619 19.619 0.000) 27.745 8.355 ( -17.362 17.362 0.000) 24.554 8.849 ( 2.098 -2.098 0.000) 2.967 9.864 ( 0.256 -0.256 0.000) 0.362 10.159 ( -10.230 10.230 0.000) 14.467 10.511 ( 7.611 -7.611 0.000) 10.763 10.720 ( 6.953 -6.953 0.000) 9.833 10.997 ( 7.821 -7.821 0.000) 11.060 11.211 ( -12.100 12.100 0.000) 17.112 11.478 ( -3.430 3.430 0.000) 4.850 11.858 ( 2.326 -2.326 0.000) 3.289 11.957 ( 3.779 -3.779 0.000) 5.344 13.413 ( 19.470 -19.470 0.000) 27.535 13.446 ( 8.365 -8.365 0.000) 11.829 14.376 ( 4.758 -4.758 0.000) 6.729 14.418 ( 13.593 -13.593 0.000) 19.224 14.514 ( -12.365 12.365 0.000) 17.487 15.168 ( 5.492 -5.492 0.000) 7.767 15.630 ( 0.652 -0.652 0.000) 0.923 16.070 ( 5.760 -5.760 0.000) 8.146 16.138 ( -1.917 1.917 0.000) 2.711 16.513 ( 2.981 -2.981 0.000) 4.216 16.750 ( -7.322 7.322 0.000) 10.355 16.866 ( -8.022 8.022 0.000) 11.344 18.199 ( 5.093 -5.093 0.000) 7.203 19.066 ( -5.815 5.815 0.000) 8.224 19.881 ( -0.595 0.595 0.000) 0.841 20.204 ( -4.106 4.106 0.000) 5.807 ======================= Grid point 46 (28/47) ======================= q-point: ( 0.40 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 102 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.198 ( 0.000 -0.000 18.629) 18.629 3.198 ( 0.000 -0.000 18.629) 18.629 3.720 ( -0.000 0.000 -4.445) 4.445 3.720 ( -0.000 0.000 -4.445) 4.445 4.756 ( 0.000 -0.000 11.073) 11.073 4.838 ( 0.000 -0.000 16.720) 16.720 5.077 ( -0.000 0.000 -3.896) 3.896 5.077 ( -0.000 0.000 -3.896) 3.896 5.137 ( 0.000 -0.000 8.415) 8.415 5.581 ( -0.000 0.000 -13.438) 13.438 5.771 ( 0.000 -0.000 27.184) 27.184 6.570 ( 0.000 -0.000 27.621) 27.621 6.570 ( 0.000 -0.000 27.621) 27.621 6.901 ( -0.000 0.000 -12.445) 12.445 8.225 ( -0.000 0.000 -26.759) 26.759 8.225 ( -0.000 0.000 -26.759) 26.759 9.151 ( 0.000 -0.000 2.511) 2.511 9.811 ( 0.000 -0.000 4.201) 4.201 9.958 ( -0.000 0.000 -26.662) 26.662 10.035 ( -0.000 0.000 -5.227) 5.227 10.886 ( 0.000 -0.000 16.355) 16.355 10.886 ( 0.000 -0.000 16.355) 16.355 11.731 ( -0.000 0.000 -9.840) 9.840 11.731 ( -0.000 0.000 -9.840) 9.840 11.825 ( 0.000 -0.000 4.028) 4.028 11.825 ( 0.000 -0.000 4.028) 4.028 13.065 ( 0.000 -0.000 34.953) 34.953 13.895 ( 0.000 -0.000 24.151) 24.151 14.337 ( -0.000 0.000 -13.446) 13.446 14.942 ( -0.000 0.000 -18.687) 18.687 15.194 ( 0.000 -0.000 28.530) 28.530 15.194 ( 0.000 -0.000 28.530) 28.530 15.838 ( 0.000 -0.000 8.320) 8.320 15.899 ( 0.000 -0.000 7.500) 7.500 16.131 ( -0.000 0.000 -2.331) 2.331 16.131 ( -0.000 0.000 -2.331) 2.331 16.693 ( -0.000 0.000 -18.373) 18.373 16.693 ( -0.000 0.000 -18.373) 18.373 18.153 ( 0.000 -0.000 12.011) 12.011 18.726 ( -0.000 0.000 -10.812) 10.812 19.908 ( 0.000 -0.000 0.095) 0.095 20.136 ( -0.000 0.000 -7.790) 7.790 ======================= Grid point 47 (29/47) ======================= q-point: ( 0.50 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.370 ( 1.730 -1.730 11.825) 12.076 3.480 ( -6.544 6.544 14.897) 17.538 3.649 ( 0.736 -0.736 -4.626) 4.741 3.730 ( -5.974 5.974 -6.248) 10.508 4.786 ( 13.149 -13.149 4.721) 19.185 4.861 ( 4.372 -4.372 11.071) 12.681 4.941 ( 5.012 -5.012 0.452) 7.102 5.113 ( 0.073 -0.073 8.175) 8.176 5.143 ( 0.497 -0.497 -8.089) 8.119 5.443 ( -1.762 1.762 -13.203) 13.436 6.199 ( -4.942 4.942 28.229) 29.081 6.769 ( -3.834 3.834 -17.393) 18.219 6.879 ( -0.644 0.644 21.907) 21.925 7.006 ( -9.480 9.480 23.490) 27.047 7.845 ( 6.131 -6.131 -22.524) 24.135 7.920 ( -0.088 0.088 -23.403) 23.404 9.215 ( 1.755 -1.755 10.818) 11.099 9.668 ( -1.667 1.667 -22.828) 22.949 9.889 ( -0.803 0.803 2.933) 3.145 10.037 ( -6.537 6.537 -3.221) 9.790 11.003 ( 4.736 -4.736 12.528) 14.206 11.090 ( -2.207 2.207 12.380) 12.767 11.569 ( 3.678 -3.678 -9.519) 10.848 11.570 ( 3.691 -3.691 -9.182) 10.562 11.867 ( -0.372 0.372 2.126) 2.191 11.894 ( -2.432 2.432 2.159) 4.060 13.492 ( -1.801 1.801 30.686) 30.792 14.125 ( -0.275 0.275 -20.351) 20.355 14.155 ( 2.334 -2.334 22.311) 22.554 14.670 ( 2.806 -2.806 -19.392) 19.793 15.359 ( 10.972 -10.972 19.504) 24.924 15.519 ( -0.072 0.072 23.763) 23.763 15.724 ( 10.865 -10.865 0.362) 15.370 16.057 ( -5.444 5.444 5.908) 9.705 16.108 ( -1.666 1.666 -7.120) 7.499 16.263 ( -8.799 8.799 1.127) 12.494 16.489 ( -0.387 0.387 -12.837) 12.849 16.506 ( -1.674 1.674 -13.993) 14.192 18.321 ( -1.178 1.178 11.844) 11.961 18.630 ( -3.395 3.395 -11.420) 12.388 19.922 ( 0.237 -0.237 2.515) 2.537 20.047 ( -0.227 0.227 -6.829) 6.836 ======================= Grid point 48 (30/47) ======================= q-point: (-0.40 -0.60 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 171 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.339 ( 6.695 -6.695 0.000) 9.469 3.625 ( -0.916 0.916 0.000) 1.295 3.758 ( -4.996 4.996 0.000) 7.066 3.880 ( -10.417 10.417 0.000) 14.732 4.477 ( 13.679 -13.679 0.000) 19.345 4.689 ( 9.411 -9.411 0.000) 13.309 4.993 ( -0.611 0.611 0.000) 0.864 5.006 ( 1.345 -1.345 0.000) 1.903 5.135 ( 3.790 -3.790 0.000) 5.360 5.378 ( -0.738 0.738 0.000) 1.044 6.531 ( -2.530 2.530 0.000) 3.578 6.643 ( 1.988 -1.988 0.000) 2.811 7.159 ( -11.389 11.389 0.000) 16.107 7.465 ( 11.817 -11.817 0.000) 16.711 7.511 ( -16.243 16.243 0.000) 22.972 7.707 ( 2.216 -2.216 0.000) 3.133 9.182 ( 8.189 -8.189 0.000) 11.581 9.642 ( -9.097 9.097 0.000) 12.865 9.893 ( 0.786 -0.786 0.000) 1.111 10.281 ( -13.698 13.698 0.000) 19.372 10.939 ( 7.804 -7.804 0.000) 11.036 11.267 ( -4.460 4.460 0.000) 6.307 11.349 ( 7.700 -7.700 0.000) 10.889 11.372 ( 6.838 -6.838 0.000) 9.670 11.890 ( -0.327 0.327 0.000) 0.463 11.974 ( -2.602 2.602 0.000) 3.680 13.739 ( 3.852 -3.852 0.000) 5.448 14.008 ( -7.415 7.415 0.000) 10.486 14.176 ( 11.341 -11.341 0.000) 16.039 14.475 ( -0.143 0.143 0.000) 0.202 15.142 ( 15.632 -15.632 0.000) 22.107 15.479 ( 8.584 -8.584 0.000) 12.140 15.685 ( 0.982 -0.982 0.000) 1.389 16.055 ( -1.358 1.358 0.000) 1.921 16.193 ( -1.369 1.369 0.000) 1.936 16.461 ( -5.934 5.934 0.000) 8.392 16.471 ( -3.719 3.719 0.000) 5.259 16.503 ( -5.796 5.796 0.000) 8.197 18.383 ( 2.625 -2.625 0.000) 3.712 18.700 ( -8.452 8.452 0.000) 11.953 19.918 ( 1.685 -1.685 0.000) 2.383 20.034 ( -2.582 2.582 0.000) 3.652 ======================= Grid point 60 (31/47) ======================= q-point: (-0.50 -0.50 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 63 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.486 ( -0.000 0.000 -0.000) 0.000 3.486 ( -0.000 0.000 -0.000) 0.000 3.607 ( 0.000 -0.000 0.000) 0.000 3.607 ( 0.000 -0.000 0.000) 0.000 5.024 ( 0.000 -0.000 0.000) 0.000 5.024 ( 0.000 -0.000 0.000) 0.000 5.040 ( -0.000 0.000 -10.752) 10.752 5.040 ( 0.000 -0.000 10.752) 10.752 5.230 ( -0.000 0.000 -15.726) 15.726 5.230 ( 0.000 -0.000 15.726) 15.726 6.454 ( -0.000 0.000 -24.831) 24.831 6.454 ( 0.000 -0.000 24.831) 24.831 7.049 ( -0.000 0.000 -0.000) 0.000 7.049 ( -0.000 0.000 -0.000) 0.000 7.757 ( 0.000 -0.000 0.000) 0.000 7.757 ( 0.000 -0.000 0.000) 0.000 9.403 ( -0.000 0.000 -18.083) 18.083 9.403 ( 0.000 -0.000 18.083) 18.083 9.917 ( -0.000 0.000 -4.636) 4.636 9.917 ( 0.000 -0.000 4.636) 4.636 11.143 ( -0.000 0.000 -0.000) 0.000 11.143 ( -0.000 0.000 -0.000) 0.000 11.553 ( 0.000 -0.000 0.000) 0.000 11.553 ( 0.000 -0.000 0.000) 0.000 11.875 ( 0.000 -0.000 0.000) 0.000 11.875 ( 0.000 -0.000 0.000) 0.000 13.827 ( -0.000 0.000 -27.661) 27.661 13.827 ( 0.000 -0.000 27.661) 27.661 14.452 ( -0.000 0.000 -22.076) 22.076 14.452 ( 0.000 -0.000 22.076) 22.076 15.703 ( -0.000 0.000 -0.000) 0.000 15.703 ( -0.000 0.000 -0.000) 0.000 15.976 ( -0.000 0.000 -1.624) 1.624 15.976 ( 0.000 -0.000 1.624) 1.624 16.021 ( 0.000 -0.000 0.000) 0.000 16.021 ( 0.000 -0.000 0.000) 0.000 16.392 ( -0.000 0.000 -0.000) 0.000 16.392 ( -0.000 0.000 -0.000) 0.000 18.448 ( -0.000 0.000 -12.255) 12.255 18.448 ( 0.000 -0.000 12.255) 12.255 19.972 ( -0.000 0.000 -5.210) 5.210 19.972 ( 0.000 -0.000 5.210) 5.210 ======================= Grid point 132 (32/47) ======================= q-point: ( 0.30 0.20 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.765 ( 0.000 -2.107 35.091) 35.154 2.118 ( 0.000 18.375 31.015) 36.050 3.475 ( 0.000 26.666 42.723) 50.362 3.846 ( 0.000 0.125 6.123) 6.124 3.861 ( 0.000 1.361 0.061) 1.362 4.077 ( 0.000 1.283 6.935) 7.053 4.682 ( -0.000 0.179 -10.520) 10.522 5.076 ( -0.000 -14.269 -1.603) 14.359 5.185 ( 0.000 -1.578 5.369) 5.596 5.309 ( 0.000 4.014 2.059) 4.511 5.373 ( 0.000 0.826 13.003) 13.029 5.744 ( 0.000 17.550 15.550) 23.448 6.261 ( 0.000 7.636 -5.738) 9.552 6.971 ( 0.000 5.737 11.219) 12.601 9.203 ( -0.000 -8.933 -16.077) 18.392 9.246 ( -0.000 -5.168 -14.886) 15.758 9.413 ( -0.000 -10.034 -11.036) 14.915 9.756 ( 0.000 9.068 2.578) 9.428 9.941 ( 0.000 -4.661 8.251) 9.477 10.065 ( 0.000 2.603 16.557) 16.760 10.502 ( 0.000 10.143 -3.552) 10.747 11.105 ( -0.000 -7.134 -18.666) 19.982 11.350 ( 0.000 1.189 24.965) 24.994 11.572 ( 0.000 0.316 6.897) 6.904 11.584 ( 0.000 3.525 5.573) 6.594 12.001 ( -0.000 -5.273 -5.496) 7.616 12.212 ( 0.000 11.384 0.709) 11.406 12.794 ( 0.000 3.387 19.999) 20.283 13.856 ( 0.000 -4.998 16.384) 17.129 13.929 ( 0.000 0.906 14.298) 14.327 14.240 ( 0.000 11.733 18.168) 21.627 15.489 ( 0.000 7.443 10.409) 12.797 15.544 ( -0.000 -7.607 -10.563) 13.017 15.716 ( -0.000 -2.046 -3.951) 4.449 16.233 ( -0.000 0.555 -1.871) 1.951 16.261 ( 0.000 3.054 -0.051) 3.054 17.240 ( -0.000 -5.151 -6.223) 8.079 17.289 ( -0.000 -2.261 -6.264) 6.660 17.746 ( 0.000 6.146 7.458) 9.663 19.214 ( -0.000 4.005 -5.651) 6.926 20.016 ( -0.000 -1.509 -3.165) 3.507 20.410 ( -0.000 -0.636 -2.850) 2.920 ======================= Grid point 133 (33/47) ======================= q-point: ( 0.40 0.20 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 500 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.195 ( -5.420 1.241 30.627) 31.127 2.639 ( 8.170 19.165 26.580) 33.772 3.785 ( 5.398 0.934 0.549) 5.506 3.867 ( 0.033 -1.437 0.609) 1.561 4.011 ( 5.504 -5.863 6.887) 10.587 4.438 ( -8.681 18.787 26.477) 33.606 4.635 ( -4.512 4.143 -5.227) 8.053 4.903 ( -1.095 -11.079 -1.211) 11.199 5.259 ( -4.023 1.495 3.187) 5.346 5.409 ( 5.109 9.740 9.425) 14.485 5.595 ( -7.751 0.492 11.316) 13.725 6.083 ( 6.536 14.387 13.648) 20.880 6.405 ( -16.309 7.215 1.935) 17.938 7.154 ( 1.587 7.886 6.667) 10.447 8.803 ( 7.315 -13.761 -21.256) 26.357 8.972 ( 0.007 -6.383 -17.814) 18.924 9.142 ( 3.333 -10.713 -8.388) 14.009 9.858 ( 5.760 5.040 2.672) 8.107 9.994 ( -1.391 -3.119 8.454) 9.118 10.287 ( 1.168 7.051 13.137) 14.955 10.685 ( -11.837 10.244 -0.296) 15.657 10.787 ( 0.170 -6.284 -18.056) 19.119 11.609 ( 1.033 -3.574 14.365) 14.839 11.644 ( 0.080 2.552 7.112) 7.557 11.700 ( 1.140 4.030 6.809) 7.994 11.825 ( 0.243 -9.298 -8.284) 12.455 12.471 ( -1.926 17.639 10.070) 20.403 13.065 ( 2.074 2.397 20.104) 20.353 14.005 ( 1.096 0.133 15.213) 15.253 14.159 ( -6.492 -0.800 14.900) 16.272 14.537 ( 7.141 9.961 15.958) 20.121 15.288 ( 0.877 -10.088 -12.118) 15.792 15.607 ( 6.854 -1.901 -3.860) 8.093 15.741 ( -4.748 7.565 10.331) 13.656 16.246 ( -5.167 0.346 -0.854) 5.248 16.333 ( -2.247 4.922 1.930) 5.745 17.079 ( 1.020 -5.832 -8.415) 10.289 17.147 ( 3.204 -3.195 -8.558) 9.681 17.931 ( -2.819 6.160 7.571) 10.160 19.181 ( 0.564 2.016 -6.407) 6.740 19.975 ( -2.624 -0.895 -2.925) 4.030 20.355 ( 0.090 -1.246 -4.345) 4.521 ======================= Grid point 134 (34/47) ======================= q-point: ( 0.50 0.20 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 500 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.620 ( -4.614 1.882 25.208) 25.696 3.009 ( 13.460 15.260 21.125) 29.331 3.690 ( 8.089 -0.709 -2.352) 8.454 3.796 ( 6.738 -4.834 1.816) 8.489 3.895 ( 10.177 -6.381 7.634) 14.233 4.731 ( -5.254 4.507 -0.171) 6.924 4.793 ( -3.520 -2.758 -0.942) 4.571 4.857 ( -2.798 1.620 4.042) 5.176 5.394 ( -8.292 6.727 -1.799) 10.828 5.638 ( -0.208 3.587 11.350) 11.905 5.753 ( -1.678 2.479 2.689) 4.024 6.445 ( -10.190 15.380 19.406) 26.777 6.778 ( -14.165 6.537 10.543) 18.829 7.256 ( 5.764 6.438 1.776) 8.822 8.216 ( 14.069 -16.200 -22.513) 31.100 8.697 ( -3.048 -5.796 -17.211) 18.415 8.914 ( 3.442 -7.222 -3.512) 8.737 9.849 ( 5.239 0.577 3.247) 6.191 10.044 ( 6.986 -1.532 9.789) 12.123 10.503 ( -4.681 -0.770 -10.724) 11.727 10.560 ( -1.707 6.231 4.106) 7.655 10.997 ( -14.112 11.671 -0.116) 18.313 11.544 ( 0.304 -11.119 -7.856) 13.618 11.644 ( -0.668 -7.150 -0.096) 7.182 11.756 ( 2.458 3.059 5.052) 6.396 11.818 ( 1.993 2.276 7.381) 7.976 12.855 ( 1.698 16.538 15.932) 23.026 13.247 ( 8.026 -1.507 17.944) 19.714 14.188 ( 0.981 4.034 12.243) 12.928 14.378 ( -4.903 -0.720 10.440) 11.557 14.659 ( 6.630 -0.585 5.669) 8.743 15.070 ( 0.347 -3.416 -3.298) 4.762 15.410 ( 12.671 -2.797 -4.113) 13.612 15.988 ( -3.769 5.124 7.587) 9.901 16.352 ( -10.225 1.408 0.506) 10.334 16.478 ( -4.598 6.062 4.157) 8.669 16.883 ( 1.784 -5.368 -9.643) 11.180 16.941 ( 4.879 -4.312 -9.507) 11.523 18.121 ( -3.145 4.596 6.064) 8.233 19.092 ( 2.195 -0.759 -5.857) 6.301 19.995 ( -6.000 0.728 -2.696) 6.618 20.279 ( -0.187 -1.560 -5.130) 5.365 ======================= Grid point 135 (35/47) ======================= q-point: (-0.40 0.20 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.855 ( 9.844 0.000 9.844) 13.921 3.142 ( 16.012 -0.000 16.013) 22.645 3.612 ( 5.031 0.000 5.031) 7.114 3.695 ( 5.311 0.000 5.311) 7.511 3.849 ( 8.979 0.000 8.979) 12.699 4.773 ( -1.818 0.000 -1.818) 2.570 4.806 ( -0.610 -0.000 -0.610) 0.863 4.896 ( -0.425 0.000 -0.425) 0.601 5.512 ( -7.579 0.000 -7.579) 10.718 5.664 ( -6.943 0.000 -6.943) 9.819 5.820 ( 8.231 0.000 8.231) 11.640 6.874 ( 3.050 -0.000 3.050) 4.313 7.080 ( 0.166 -0.000 0.166) 0.235 7.174 ( 2.224 0.000 2.224) 3.145 7.768 ( -1.233 0.000 -1.233) 1.743 8.561 ( -10.990 0.000 -10.990) 15.542 8.827 ( 1.534 0.000 1.534) 2.169 9.840 ( 4.122 0.000 4.122) 5.830 10.047 ( 8.965 -0.000 8.965) 12.679 10.415 ( -9.224 0.000 -9.224) 13.044 10.705 ( 3.469 -0.000 3.469) 4.905 11.269 ( -7.393 -0.000 -7.393) 10.456 11.296 ( -7.968 0.000 -7.968) 11.269 11.602 ( 0.183 0.000 0.183) 0.258 11.790 ( 3.375 0.000 3.375) 4.772 11.864 ( 4.925 0.000 4.925) 6.965 13.087 ( 9.726 -0.000 9.726) 13.755 13.269 ( 14.938 0.000 14.938) 21.125 14.310 ( 8.164 0.000 8.164) 11.546 14.520 ( -0.620 -0.000 -0.620) 0.878 14.532 ( -0.832 0.000 -0.832) 1.177 15.153 ( 4.800 0.000 4.800) 6.789 15.186 ( 4.907 -0.000 4.907) 6.939 16.093 ( 1.922 0.000 1.922) 2.719 16.464 ( -5.958 0.000 -5.958) 8.425 16.637 ( 0.544 0.000 0.544) 0.770 16.743 ( -3.474 0.000 -3.474) 4.913 16.758 ( -3.013 -0.000 -3.013) 4.261 18.215 ( 1.126 0.000 1.126) 1.593 19.025 ( -1.320 -0.000 -1.320) 1.867 20.033 ( -5.123 0.000 -5.123) 7.246 20.229 ( -2.681 0.000 -2.681) 3.791 ======================= Grid point 144 (36/47) ======================= q-point: ( 0.40 0.30 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 280 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.572 ( 0.000 -3.774 31.110) 31.338 2.806 ( 0.000 13.256 25.548) 28.782 3.810 ( 0.000 1.826 -0.627) 1.931 3.844 ( 0.000 2.136 -1.239) 2.469 4.280 ( 0.000 -10.664 13.377) 17.107 4.544 ( 0.000 10.718 30.653) 32.472 4.559 ( -0.000 -0.965 0.600) 1.136 5.004 ( -0.000 -11.720 -4.344) 12.499 5.263 ( 0.000 1.897 0.829) 2.070 5.455 ( 0.000 8.963 12.463) 15.351 5.698 ( 0.000 -4.815 7.447) 8.868 6.182 ( 0.000 8.435 5.138) 9.876 6.189 ( 0.000 16.795 21.298) 27.123 7.121 ( 0.000 6.258 0.867) 6.318 8.725 ( -0.000 -9.326 -23.242) 25.043 8.783 ( -0.000 -4.512 -23.074) 23.511 9.170 ( -0.000 -6.772 -7.889) 10.397 9.828 ( 0.000 8.952 3.326) 9.550 10.114 ( -0.000 -5.469 3.392) 6.436 10.476 ( 0.000 0.618 10.427) 10.445 10.508 ( 0.000 5.447 6.479) 8.465 10.600 ( -0.000 -0.632 -16.338) 16.350 11.709 ( 0.000 1.653 7.317) 7.502 11.736 ( 0.000 3.301 6.370) 7.175 11.821 ( -0.000 -6.712 3.959) 7.793 11.831 ( -0.000 -8.470 -8.675) 12.124 12.434 ( 0.000 15.695 23.417) 28.190 13.319 ( 0.000 0.726 23.120) 23.131 14.288 ( 0.000 -3.901 12.018) 12.635 14.358 ( 0.000 -9.461 24.706) 26.456 14.684 ( 0.000 8.099 18.783) 20.455 15.250 ( -0.000 -7.571 -13.533) 15.507 15.640 ( -0.000 -2.657 -1.421) 3.013 15.755 ( 0.000 6.380 11.217) 12.904 16.190 ( -0.000 0.778 -1.564) 1.747 16.270 ( 0.000 6.261 1.105) 6.357 17.025 ( -0.000 -4.260 -12.044) 12.775 17.067 ( -0.000 -1.313 -12.395) 12.464 17.949 ( 0.000 4.875 9.549) 10.721 19.037 ( -0.000 5.064 -8.897) 10.237 19.939 ( -0.000 -0.969 -2.837) 2.998 20.304 ( -0.000 -0.555 -6.071) 6.096 ======================= Grid point 145 (37/47) ======================= q-point: ( 0.50 0.30 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 500 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.913 ( -1.499 -1.407 27.402) 27.479 3.209 ( 5.021 12.672 20.258) 24.417 3.752 ( 7.485 2.033 -2.354) 8.105 3.901 ( -2.367 -0.049 1.680) 2.903 4.234 ( 8.239 -11.851 11.324) 18.346 4.605 ( -3.517 2.485 2.710) 5.088 4.798 ( 0.368 -3.409 0.584) 3.478 4.886 ( -4.291 -4.579 3.740) 7.306 5.277 ( -1.368 0.725 -2.499) 2.939 5.660 ( -4.199 -6.052 -6.631) 9.911 5.833 ( 6.954 16.963 23.360) 29.695 6.435 ( -5.971 8.027 20.457) 22.773 6.692 ( -10.897 10.566 18.351) 23.815 7.197 ( 1.533 9.440 -3.211) 10.089 8.239 ( 9.319 -12.313 -25.147) 29.510 8.467 ( -4.276 -4.947 -23.860) 24.740 9.017 ( 3.284 -6.186 -0.364) 7.013 9.927 ( 6.861 6.267 2.674) 9.669 10.179 ( -7.519 -3.108 2.066) 8.394 10.269 ( -4.688 0.187 -18.145) 18.742 10.573 ( 7.730 5.527 8.633) 12.838 10.832 ( -7.061 4.924 12.400) 15.096 11.570 ( -0.303 -13.426 -11.502) 17.682 11.725 ( 2.123 -4.805 -2.635) 5.877 11.814 ( 2.115 2.617 5.064) 6.080 11.869 ( -1.076 3.459 5.770) 6.812 12.925 ( -0.901 13.082 25.822) 28.961 13.571 ( 3.533 -1.333 21.314) 21.646 14.321 ( 2.166 -1.084 4.181) 4.831 14.563 ( -1.381 -7.959 20.646) 22.170 14.816 ( 4.130 -7.976 -1.845) 9.169 15.098 ( 1.515 2.295 6.994) 7.515 15.534 ( 9.303 -2.645 -1.905) 9.857 15.981 ( -3.873 5.167 9.167) 11.213 16.224 ( -7.366 0.460 -1.488) 7.529 16.400 ( -2.940 7.286 3.810) 8.731 16.815 ( 0.822 -3.804 -13.630) 14.174 16.882 ( 1.956 -1.276 -13.109) 13.315 18.131 ( -1.215 4.396 9.335) 10.389 18.987 ( -0.984 3.627 -9.569) 10.281 19.915 ( -1.971 -0.201 -1.640) 2.572 20.217 ( -0.237 -0.723 -6.908) 6.949 ======================= Grid point 146 (38/47) ======================= q-point: ( 0.60 0.30 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 500 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.132 ( 11.384 0.083 13.089) 17.348 3.452 ( 2.704 0.040 16.167) 16.392 3.712 ( 5.718 6.640 5.249) 10.215 3.885 ( 5.552 -6.410 5.116) 9.904 4.120 ( 11.073 -6.493 8.830) 15.580 4.690 ( -1.332 2.247 -5.829) 6.387 4.782 ( -7.154 -2.979 3.043) 8.325 4.912 ( 1.021 -1.340 2.257) 2.816 5.277 ( -5.819 4.284 -7.697) 10.557 5.522 ( -5.824 -6.443 -11.880) 14.716 6.106 ( 12.707 4.222 20.497) 24.483 6.895 ( -1.665 10.342 11.799) 15.778 7.072 ( -8.322 1.768 7.507) 11.347 7.222 ( 1.257 0.788 2.556) 2.955 7.725 ( 4.823 -4.560 -15.517) 16.877 8.240 ( -11.115 -0.717 -20.256) 23.117 8.935 ( 6.804 -3.102 4.988) 8.989 9.929 ( 7.044 4.186 2.656) 8.613 10.063 ( -8.438 0.837 -14.430) 16.737 10.354 ( -2.900 -6.564 8.722) 11.294 10.728 ( 0.274 13.150 5.566) 14.282 11.055 ( -2.749 2.638 5.734) 6.884 11.262 ( -2.664 -14.918 -13.005) 19.969 11.597 ( 3.756 0.324 -5.395) 6.582 11.865 ( 2.564 1.205 3.839) 4.771 11.958 ( 1.211 0.717 4.164) 4.396 13.312 ( 3.836 4.340 19.873) 20.700 13.666 ( 12.510 -5.336 15.783) 20.835 14.348 ( -4.435 4.777 -3.392) 7.348 14.562 ( 1.095 -9.592 -8.990) 13.192 14.643 ( 10.071 -5.507 13.313) 17.578 15.299 ( -0.265 1.957 17.200) 17.313 15.342 ( 12.838 -2.401 -0.347) 13.066 16.135 ( -0.236 1.023 4.097) 4.229 16.294 ( -7.863 0.364 -7.271) 10.716 16.587 ( -6.416 4.773 3.073) 8.567 16.643 ( -2.625 0.649 -7.266) 7.753 16.698 ( 2.331 -1.581 -9.388) 9.801 18.264 ( 2.101 1.285 5.676) 6.187 18.945 ( -3.957 2.157 -5.777) 7.327 19.928 ( -2.457 0.504 -2.567) 3.589 20.151 ( -2.490 0.298 -4.805) 5.420 ======================= Grid point 147 (39/47) ======================= q-point: (-0.30 0.30 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 500 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.934 ( 24.327 -3.837 0.028) 24.628 3.534 ( 7.716 -12.263 11.259) 18.349 3.780 ( 7.746 1.239 5.781) 9.744 3.861 ( -1.682 0.458 7.247) 7.453 4.047 ( 10.229 2.119 5.621) 11.862 4.631 ( 0.376 3.038 -8.376) 8.918 4.794 ( -2.936 -7.651 1.052) 8.262 4.903 ( 2.140 -2.228 -0.047) 3.090 5.326 ( -4.899 -2.237 -4.855) 7.251 5.452 ( -7.162 9.224 -8.059) 14.189 6.150 ( 16.050 -7.405 18.748) 25.767 6.663 ( 19.564 -10.491 -21.512) 30.912 6.876 ( 11.001 -6.143 -1.923) 12.746 7.418 ( -2.133 -7.116 7.735) 10.724 7.905 ( -20.256 14.959 6.761) 26.072 8.313 ( -21.154 8.076 -8.057) 24.034 8.894 ( 5.573 2.399 2.680) 6.632 9.933 ( 6.534 4.949 2.622) 8.606 10.124 ( -12.386 3.662 -3.442) 13.367 10.344 ( -1.072 -9.894 4.105) 10.765 10.791 ( 10.923 -5.581 4.804) 13.173 10.987 ( 1.883 -10.141 -11.191) 15.219 11.221 ( -12.788 13.170 -1.562) 18.423 11.575 ( 1.858 7.291 -1.301) 7.636 11.865 ( 3.730 -1.118 2.015) 4.384 11.948 ( 4.162 -2.939 1.904) 5.439 13.290 ( 13.787 -15.790 5.406) 21.648 13.562 ( 18.654 -1.746 7.670) 20.244 14.374 ( 1.788 -4.926 -2.106) 5.648 14.503 ( 17.800 -3.353 1.877) 18.210 14.539 ( -9.554 12.021 -0.750) 15.374 15.184 ( 10.357 -3.497 2.875) 11.303 15.517 ( -1.978 -0.370 13.031) 13.185 16.092 ( 5.753 -4.039 -1.623) 7.214 16.217 ( -1.777 0.676 -10.838) 11.003 16.534 ( 3.260 -3.769 -3.959) 6.365 16.717 ( -9.651 4.347 -0.027) 10.585 16.814 ( -10.105 4.817 2.844) 11.550 18.225 ( 6.825 -3.126 0.041) 7.507 19.017 ( -8.071 1.523 0.238) 8.217 19.907 ( -0.927 0.419 -3.781) 3.915 20.186 ( -5.611 1.961 -0.914) 6.013 ======================= Grid point 148 (40/47) ======================= q-point: (-0.20 0.30 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.558 ( 28.318 -4.225 -0.000) 28.631 3.213 ( 20.413 -19.384 -0.000) 28.150 3.823 ( 0.775 4.158 -0.000) 4.230 3.884 ( 1.028 -3.509 0.000) 3.657 4.057 ( 6.109 6.399 -0.000) 8.847 4.594 ( -2.129 -1.390 0.000) 2.542 4.771 ( 10.656 -7.489 -0.000) 13.024 4.784 ( -1.596 2.406 -0.000) 2.887 5.274 ( 1.776 -1.253 -0.000) 2.173 5.581 ( 0.375 6.374 -0.000) 6.385 5.861 ( 21.932 -21.195 -0.000) 30.500 6.274 ( 5.943 -0.845 -0.000) 6.003 6.574 ( 15.812 -13.635 -0.000) 20.879 7.333 ( 2.800 -9.949 0.000) 10.335 8.414 ( -22.134 16.167 0.000) 27.410 8.639 ( -22.083 7.588 0.000) 23.350 8.961 ( -2.604 9.356 -0.000) 9.712 9.948 ( 6.223 6.601 -0.000) 9.072 10.111 ( 2.717 -12.966 0.000) 13.248 10.439 ( -14.253 3.649 0.000) 14.713 10.618 ( 13.956 -2.495 -0.000) 14.177 10.730 ( -0.211 -11.941 0.000) 11.943 11.548 ( -12.727 13.995 -0.000) 18.917 11.629 ( 4.478 7.623 -0.000) 8.841 11.791 ( 6.625 -3.202 -0.000) 7.358 11.846 ( 6.023 -5.024 -0.000) 7.843 12.883 ( 17.211 -21.258 0.000) 27.352 13.351 ( 19.780 -2.221 -0.000) 19.904 14.237 ( 19.781 -1.006 -0.000) 19.807 14.270 ( 8.246 -2.506 -0.000) 8.618 14.765 ( 0.635 5.062 -0.000) 5.102 15.090 ( -1.153 1.198 -0.000) 1.662 15.615 ( -2.009 -0.177 0.000) 2.017 15.916 ( 8.973 -6.796 -0.000) 11.256 16.182 ( -1.544 1.429 0.000) 2.104 16.413 ( 3.573 -5.226 0.000) 6.331 16.904 ( -10.080 5.693 0.000) 11.577 17.010 ( -10.401 3.906 0.000) 11.111 18.078 ( 7.822 -5.333 -0.000) 9.467 19.124 ( -7.773 -1.031 0.000) 7.841 19.908 ( -1.925 1.445 0.000) 2.407 20.276 ( -5.600 1.937 0.000) 5.925 ======================= Grid point 155 (41/47) ======================= q-point: ( 0.50 0.40 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.218 ( 0.000 1.140 22.251) 22.280 3.313 ( 0.000 8.950 16.449) 18.726 3.748 ( -0.000 2.546 -4.023) 4.761 3.794 ( 0.000 5.242 -3.476) 6.290 4.598 ( -0.000 -14.317 11.759) 18.527 4.701 ( 0.000 -4.012 10.128) 10.894 4.887 ( -0.000 -11.305 -4.652) 12.224 4.971 ( 0.000 -3.898 8.876) 9.694 5.228 ( -0.000 0.072 -5.059) 5.059 5.573 ( -0.000 -3.261 -12.469) 12.888 5.975 ( 0.000 16.649 28.100) 32.662 6.627 ( 0.000 3.037 25.602) 25.781 6.737 ( 0.000 10.793 23.540) 25.897 7.006 ( -0.000 9.033 -10.350) 13.737 8.127 ( -0.000 -7.886 -25.576) 26.764 8.171 ( -0.000 -4.246 -27.085) 27.416 9.116 ( 0.000 -2.723 4.632) 5.374 9.910 ( 0.000 7.301 2.009) 7.572 9.962 ( -0.000 0.897 -23.801) 23.818 10.095 ( 0.000 4.084 -3.192) 5.183 10.793 ( 0.000 -6.274 14.699) 15.982 10.932 ( 0.000 3.587 16.144) 16.538 11.586 ( -0.000 -11.522 -10.853) 15.829 11.710 ( -0.000 -1.805 -9.005) 9.184 11.857 ( 0.000 2.438 4.200) 4.856 11.868 ( 0.000 3.344 3.899) 5.137 13.149 ( 0.000 6.689 31.685) 32.383 13.873 ( 0.000 -1.728 22.719) 22.785 14.309 ( -0.000 -1.794 -10.881) 11.028 14.835 ( -0.000 -8.532 -15.635) 17.811 14.984 ( 0.000 -12.425 25.867) 28.697 15.256 ( 0.000 3.757 22.628) 22.938 15.654 ( -0.000 -8.111 1.165) 8.194 16.001 ( 0.000 5.394 9.135) 10.609 16.143 ( -0.000 0.448 -3.565) 3.593 16.323 ( 0.000 9.716 3.121) 10.205 16.668 ( -0.000 -1.874 -17.053) 17.155 16.713 ( -0.000 1.483 -15.835) 15.904 18.204 ( 0.000 3.722 11.729) 12.306 18.793 ( -0.000 4.657 -11.383) 12.299 19.903 ( -0.000 -0.368 0.350) 0.508 20.133 ( -0.000 -0.323 -7.667) 7.674 ======================= Grid point 156 (42/47) ======================= q-point: ( 0.60 0.40 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 500 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.401 ( 12.255 3.524 10.875) 16.759 3.604 ( -2.588 1.971 11.681) 12.126 3.706 ( 0.992 6.121 0.427) 6.216 3.897 ( -8.775 5.813 -3.536) 11.104 4.457 ( 11.907 -11.793 5.105) 17.518 4.672 ( 2.380 -6.826 -0.805) 7.273 4.833 ( -2.566 -4.368 9.011) 10.337 5.018 ( 0.532 -6.313 6.301) 8.936 5.124 ( -0.582 -0.046 -8.539) 8.559 5.382 ( -5.637 -6.670 -13.036) 15.691 6.433 ( 9.061 11.911 24.147) 28.409 6.917 ( 9.813 8.112 -7.814) 14.939 6.998 ( -10.380 6.911 17.401) 21.408 7.234 ( -13.106 6.752 19.257) 24.253 7.651 ( 8.996 -8.629 -21.970) 25.260 7.866 ( -7.897 -3.873 -24.984) 26.487 9.156 ( 7.570 -2.543 10.485) 13.180 9.722 ( -8.073 2.167 -19.001) 20.759 9.979 ( 8.124 6.546 1.109) 10.492 10.259 ( -13.863 9.557 0.407) 16.843 10.839 ( 5.082 -6.835 12.251) 14.921 11.136 ( -2.403 0.376 5.438) 5.957 11.310 ( 2.146 -13.141 -3.663) 13.809 11.546 ( 3.994 -0.282 -10.647) 11.375 11.908 ( 2.107 2.482 2.205) 3.932 11.966 ( -2.045 2.917 2.097) 4.134 13.564 ( 1.586 3.716 24.887) 25.213 14.050 ( 6.894 -5.867 17.341) 19.562 14.146 ( -3.191 1.850 -15.487) 15.920 14.533 ( -0.733 -8.307 -14.420) 16.658 15.060 ( 10.648 -10.469 16.490) 22.246 15.443 ( 1.090 -5.520 1.159) 5.745 15.578 ( 7.601 -0.254 15.983) 17.700 16.098 ( -6.339 -0.477 -10.089) 11.925 16.174 ( -1.007 2.203 6.268) 6.719 16.478 ( -3.432 -0.517 -11.346) 11.865 16.493 ( -3.473 7.972 2.457) 9.036 16.594 ( 0.042 3.695 -8.141) 8.940 18.370 ( 3.546 2.420 9.499) 10.425 18.743 ( -6.079 4.723 -9.600) 12.305 19.914 ( 0.772 -0.323 1.596) 1.802 20.055 ( -2.095 0.317 -5.608) 5.995 ======================= Grid point 157 (43/47) ======================= q-point: (-0.30 0.40 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 275 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.251 ( 21.221 -1.441 -0.000) 21.270 3.671 ( 0.606 1.386 -0.000) 1.513 3.856 ( -5.119 2.356 0.000) 5.635 3.963 ( 2.763 -3.350 0.000) 4.343 4.231 ( 12.443 -5.018 -0.000) 13.417 4.569 ( 4.124 0.777 -0.000) 4.197 4.878 ( -5.269 -11.110 0.000) 12.296 4.955 ( 0.840 -2.527 0.000) 2.663 5.154 ( -5.823 8.065 -0.000) 9.948 5.298 ( -6.053 -4.891 0.000) 7.782 6.637 ( 18.044 1.318 -0.000) 18.092 6.643 ( 10.868 1.974 -0.000) 11.046 7.144 ( 14.038 -12.790 -0.000) 18.991 7.480 ( -17.226 8.929 0.000) 19.403 7.594 ( -3.918 -2.345 0.000) 4.567 7.846 ( -23.707 10.895 0.000) 26.091 9.043 ( 12.875 -2.946 -0.000) 13.207 9.826 ( -16.380 5.991 0.000) 17.441 9.964 ( 7.765 6.175 -0.000) 9.921 10.493 ( -8.903 -5.023 0.000) 10.222 10.843 ( 1.633 8.526 -0.000) 8.681 11.025 ( -3.792 -14.660 0.000) 15.142 11.178 ( 9.837 -5.903 -0.000) 11.472 11.483 ( 1.067 7.498 -0.000) 7.573 11.915 ( 2.718 1.137 -0.000) 2.946 12.009 ( 0.399 -0.181 -0.000) 0.438 13.625 ( 10.066 -7.407 -0.000) 12.497 13.903 ( 17.352 -8.801 -0.000) 19.456 14.213 ( -12.437 9.283 0.000) 15.519 14.415 ( -3.102 -3.137 0.000) 4.412 14.838 ( 19.791 -8.296 -0.000) 21.459 15.344 ( 7.164 -2.043 -0.000) 7.450 15.638 ( 1.554 -3.307 0.000) 3.654 16.082 ( -4.168 1.710 0.000) 4.505 16.179 ( 3.605 -2.264 -0.000) 4.257 16.512 ( -6.858 1.673 0.000) 7.059 16.588 ( -2.201 2.261 -0.000) 3.156 16.662 ( -7.346 5.575 0.000) 9.222 18.348 ( 7.942 -1.182 -0.000) 8.030 18.861 ( -10.740 4.288 0.000) 11.565 19.893 ( 1.785 -0.429 -0.000) 1.836 20.087 ( -5.318 1.808 0.000) 5.617 ======================= Grid point 166 (44/47) ======================= q-point: (-0.40 -0.50 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 155 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.578 ( 0.000 6.818 -0.263) 6.823 3.578 ( 0.000 6.818 0.263) 6.823 3.653 ( 0.000 3.459 -4.105) 5.368 3.653 ( 0.000 3.459 4.105) 5.368 4.770 ( 0.000 -13.725 -3.141) 14.079 4.770 ( 0.000 -13.725 3.141) 14.079 5.008 ( -0.000 -2.605 -12.056) 12.335 5.008 ( 0.000 -2.605 12.056) 12.335 5.238 ( -0.000 -4.098 -13.814) 14.409 5.238 ( 0.000 -4.098 13.814) 14.409 6.603 ( -0.000 11.924 -22.005) 25.028 6.603 ( 0.000 11.924 22.005) 25.028 7.144 ( 0.000 6.326 -1.855) 6.593 7.144 ( 0.000 6.326 1.855) 6.593 7.673 ( 0.000 -6.821 -0.153) 6.823 7.673 ( 0.000 -6.821 0.153) 6.823 9.407 ( -0.000 0.303 -18.667) 18.670 9.407 ( 0.000 0.303 18.667) 18.670 10.006 ( 0.000 7.121 -3.150) 7.787 10.006 ( 0.000 7.121 3.150) 7.787 11.106 ( 0.000 -4.408 -7.234) 8.471 11.106 ( 0.000 -4.408 7.234) 8.471 11.472 ( -0.000 -4.770 -7.130) 8.578 11.472 ( 0.000 -4.770 7.130) 8.578 11.912 ( 0.000 2.876 -0.095) 2.877 11.912 ( 0.000 2.876 0.095) 2.877 13.850 ( -0.000 1.803 -25.509) 25.572 13.850 ( 0.000 1.803 25.509) 25.572 14.392 ( -0.000 -4.961 -20.262) 20.860 14.392 ( 0.000 -4.961 20.262) 20.860 15.476 ( 0.000 -12.865 -11.525) 17.273 15.476 ( 0.000 -12.865 11.525) 17.273 15.873 ( 0.000 -1.183 -18.366) 18.404 15.873 ( 0.000 -1.183 18.366) 18.404 16.253 ( 0.000 5.754 -12.315) 13.593 16.253 ( 0.000 5.754 12.315) 13.593 16.444 ( 0.000 5.405 -4.988) 7.355 16.444 ( 0.000 5.405 4.988) 7.355 18.501 ( 0.000 3.833 -12.481) 13.057 18.501 ( 0.000 3.833 12.481) 13.057 19.970 ( -0.000 -0.117 -5.137) 5.138 19.970 ( 0.000 -0.117 5.137) 5.138 ======================= Grid point 262 (45/47) ======================= q-point: ( 0.60 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 282 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.235 ( 0.000 0.365 26.654) 26.656 3.566 ( -0.000 10.153 9.029) 13.587 3.834 ( -0.000 3.649 -4.028) 5.435 3.893 ( 0.000 1.975 2.628) 3.287 4.292 ( 0.000 -9.467 11.581) 14.958 4.660 ( 0.000 1.464 6.670) 6.829 4.672 ( 0.000 -4.274 -5.479) 6.949 4.917 ( -0.000 -0.076 3.707) 3.707 5.222 ( 0.000 0.991 -6.770) 6.842 5.399 ( 0.000 -10.332 -12.866) 16.501 6.426 ( 0.000 17.889 23.550) 29.574 6.662 ( 0.000 -0.117 23.286) 23.286 6.952 ( 0.000 6.086 15.958) 17.079 7.288 ( 0.000 12.699 -6.684) 14.351 7.869 ( 0.000 -13.148 -22.333) 25.916 8.042 ( -0.000 -6.307 -27.070) 27.795 9.044 ( 0.000 -2.448 9.958) 10.255 9.925 ( 0.000 -2.400 -21.253) 21.388 10.120 ( 0.000 -4.194 6.657) 7.868 10.198 ( 0.000 13.937 0.998) 13.973 10.746 ( 0.000 5.752 3.007) 6.490 11.038 ( 0.000 4.438 16.662) 17.243 11.229 ( 0.000 -16.715 -12.994) 21.172 11.635 ( 0.000 -4.750 -4.835) 6.778 11.913 ( 0.000 1.554 4.149) 4.430 11.959 ( 0.000 3.595 3.037) 4.706 13.371 ( 0.000 11.381 22.634) 25.334 13.829 ( 0.000 -1.041 19.092) 19.120 14.289 ( 0.000 0.597 -8.653) 8.673 14.559 ( 0.000 -12.901 -12.184) 17.745 14.776 ( 0.000 -3.352 23.488) 23.726 15.322 ( 0.000 1.248 13.714) 13.770 15.569 ( -0.000 0.940 -0.665) 1.151 16.107 ( 0.000 3.580 7.118) 7.967 16.128 ( 0.000 -1.637 -4.704) 4.980 16.497 ( 0.000 4.097 2.366) 4.731 16.620 ( -0.000 -1.689 -12.436) 12.550 16.734 ( 0.000 -0.535 -10.828) 10.841 18.297 ( 0.000 3.309 10.705) 11.205 18.887 ( 0.000 2.043 -11.676) 11.853 19.893 ( 0.000 -0.356 1.047) 1.106 20.119 ( 0.000 -0.880 -7.035) 7.090 ======================= Grid point 263 (46/47) ======================= q-point: (-0.30 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.464 ( 14.048 0.000 14.048) 19.866 3.619 ( 5.803 0.000 5.803) 8.207 3.821 ( -1.249 0.000 -1.249) 1.767 3.946 ( -0.147 -0.000 -0.147) 0.208 4.301 ( 8.588 0.000 8.588) 12.146 4.613 ( -4.332 -0.000 -4.332) 6.126 4.752 ( 3.397 0.000 3.397) 4.805 4.932 ( 1.080 -0.000 1.080) 1.527 5.156 ( -5.252 0.000 -5.252) 7.428 5.240 ( -9.696 -0.000 -9.696) 13.712 6.605 ( 21.170 -0.000 21.170) 29.938 7.062 ( 1.338 0.000 1.338) 1.892 7.153 ( 5.992 -0.000 5.992) 8.475 7.211 ( 0.321 0.000 0.321) 0.454 7.554 ( -7.477 0.000 -7.477) 10.574 7.799 ( -19.782 0.000 -19.782) 27.976 9.121 ( 10.754 0.000 10.754) 15.209 9.750 ( -14.158 0.000 -14.158) 20.022 10.114 ( 6.669 0.000 6.669) 9.432 10.345 ( -0.876 -0.000 -0.876) 1.238 10.859 ( -1.756 -0.000 -1.756) 2.483 10.951 ( -5.422 0.000 -5.422) 7.668 11.203 ( 9.024 0.000 9.024) 12.762 11.547 ( -3.706 0.000 -3.706) 5.241 11.946 ( 2.603 0.000 2.603) 3.681 12.004 ( 0.571 0.000 0.571) 0.808 13.628 ( 12.336 -0.000 12.336) 17.445 13.951 ( 11.957 0.000 11.957) 16.909 14.203 ( -10.384 0.000 -10.384) 14.685 14.375 ( -5.768 0.000 -5.768) 8.158 14.931 ( 14.041 0.000 14.041) 19.856 15.381 ( -1.098 0.000 -1.098) 1.552 15.581 ( 11.586 0.000 11.586) 16.385 16.088 ( -8.696 0.000 -8.696) 12.299 16.197 ( 2.645 0.000 2.645) 3.740 16.471 ( -6.883 0.000 -6.883) 9.734 16.602 ( -1.682 0.000 -1.682) 2.378 16.634 ( -1.488 -0.000 -1.488) 2.105 18.405 ( 7.081 0.000 7.081) 10.014 18.806 ( -7.933 0.000 -7.933) 11.219 19.910 ( 1.025 0.000 1.025) 1.449 20.059 ( -4.050 0.000 -4.050) 5.727 ======================= Grid point 274 (47/47) ======================= q-point: (-0.30 -0.50 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 3.36e-05 5.98e-04 5.98e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 155 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.715 ( 0.000 2.670 -2.282) 3.512 3.715 ( 0.000 2.670 2.281) 3.512 3.766 ( 0.000 4.926 -12.247) 13.201 3.766 ( 0.000 4.926 12.247) 13.201 4.520 ( 0.000 -5.762 -5.988) 8.310 4.520 ( 0.000 -5.762 5.988) 8.310 4.925 ( -0.000 -3.521 -10.727) 11.290 4.925 ( 0.000 -3.521 10.727) 11.290 5.106 ( -0.000 -3.910 -7.769) 8.697 5.106 ( 0.000 -3.910 7.769) 8.697 6.930 ( -0.000 12.268 -15.734) 19.951 6.930 ( 0.000 12.268 15.734) 19.951 7.269 ( 0.000 2.795 -2.729) 3.906 7.269 ( 0.000 2.795 2.729) 3.906 7.463 ( 0.000 -8.734 -3.404) 9.374 7.463 ( 0.000 -8.734 3.404) 9.374 9.414 ( 0.000 0.183 -18.979) 18.980 9.414 ( 0.000 0.183 18.979) 18.980 10.212 ( 0.000 7.764 -0.197) 7.766 10.212 ( 0.000 7.764 0.197) 7.766 10.915 ( -0.000 -8.465 -10.723) 13.662 10.915 ( 0.000 -8.465 10.723) 13.662 11.414 ( 0.000 -0.381 -12.615) 12.621 11.414 ( 0.000 -0.381 12.615) 12.621 11.982 ( 0.000 2.027 -0.969) 2.247 11.982 ( 0.000 2.027 0.969) 2.247 13.907 ( 0.000 2.526 -20.392) 20.548 13.907 ( 0.000 2.526 20.392) 20.548 14.241 ( 0.000 -5.851 -14.667) 15.791 14.241 ( 0.000 -5.851 14.667) 15.791 15.262 ( 0.000 -4.327 -14.368) 15.005 15.262 ( 0.000 -4.327 14.368) 15.005 15.857 ( 0.000 -0.356 -16.671) 16.674 15.857 ( 0.000 -0.356 16.671) 16.674 16.304 ( 0.000 0.448 -10.133) 10.143 16.304 ( 0.000 0.448 10.133) 10.143 16.577 ( 0.000 3.495 -2.863) 4.518 16.577 ( 0.000 3.495 2.863) 4.518 18.585 ( 0.000 2.281 -12.429) 12.636 18.585 ( 0.000 2.281 12.429) 12.636 19.967 ( 0.000 -0.078 -4.870) 4.870 19.967 ( 0.000 -0.078 4.870) 4.870 =================== 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/42000 10.0 2557.092 2557.092 2557.092 -0.000 -0.000 0.000 3/42000 20.0 2800.671 2800.671 2800.671 -0.000 -0.000 0.000 3/42000 30.0 1198.070 1198.070 1198.070 -0.000 -0.000 0.000 3/42000 40.0 633.971 633.971 633.971 -0.000 -0.000 0.000 3/42000 50.0 397.178 397.178 397.178 -0.000 -0.000 0.000 3/42000 60.0 275.393 275.393 275.393 -0.000 -0.000 0.000 3/42000 70.0 203.325 203.325 203.325 -0.000 -0.000 0.000 3/42000 80.0 156.706 156.706 156.706 -0.000 -0.000 0.000 3/42000 90.0 124.775 124.775 124.775 -0.000 -0.000 0.000 3/42000 100.0 102.006 102.006 102.006 -0.000 -0.000 0.000 3/42000 110.0 85.258 85.258 85.258 -0.000 -0.000 0.000 3/42000 120.0 72.616 72.616 72.616 -0.000 -0.000 0.000 3/42000 130.0 62.858 62.858 62.858 -0.000 -0.000 0.000 3/42000 140.0 55.176 55.176 55.176 -0.000 -0.000 0.000 3/42000 150.0 49.021 49.021 49.021 -0.000 -0.000 0.000 3/42000 160.0 44.011 44.011 44.011 -0.000 -0.000 0.000 3/42000 170.0 39.874 39.874 39.874 -0.000 -0.000 0.000 3/42000 180.0 36.415 36.415 36.415 -0.000 -0.000 0.000 3/42000 190.0 33.488 33.488 33.488 -0.000 -0.000 0.000 3/42000 200.0 30.985 30.985 30.985 -0.000 -0.000 0.000 3/42000 210.0 28.825 28.825 28.825 -0.000 -0.000 0.000 3/42000 220.0 26.945 26.945 26.945 -0.000 -0.000 0.000 3/42000 230.0 25.294 25.294 25.294 -0.000 -0.000 0.000 3/42000 240.0 23.836 23.836 23.836 -0.000 -0.000 0.000 3/42000 250.0 22.538 22.538 22.538 -0.000 -0.000 0.000 3/42000 260.0 21.377 21.377 21.377 -0.000 -0.000 0.000 3/42000 270.0 20.332 20.332 20.332 -0.000 -0.000 0.000 3/42000 280.0 19.387 19.387 19.387 -0.000 -0.000 0.000 3/42000 290.0 18.528 18.528 18.528 -0.000 -0.000 0.000 3/42000 300.0 17.744 17.744 17.744 -0.000 -0.000 0.000 3/42000 310.0 17.026 17.026 17.026 -0.000 -0.000 0.000 3/42000 320.0 16.366 16.366 16.366 -0.000 -0.000 0.000 3/42000 330.0 15.757 15.757 15.757 -0.000 -0.000 0.000 3/42000 340.0 15.193 15.193 15.193 -0.000 -0.000 0.000 3/42000 350.0 14.669 14.669 14.669 -0.000 -0.000 0.000 3/42000 360.0 14.182 14.182 14.182 -0.000 -0.000 0.000 3/42000 370.0 13.727 13.727 13.727 -0.000 -0.000 0.000 3/42000 380.0 13.302 13.302 13.302 -0.000 -0.000 0.000 3/42000 390.0 12.903 12.903 12.903 -0.000 -0.000 0.000 3/42000 400.0 12.529 12.529 12.529 -0.000 -0.000 0.000 3/42000 410.0 12.176 12.176 12.176 -0.000 -0.000 0.000 3/42000 420.0 11.843 11.843 11.843 -0.000 -0.000 0.000 3/42000 430.0 11.529 11.529 11.529 -0.000 -0.000 0.000 3/42000 440.0 11.232 11.232 11.232 -0.000 -0.000 0.000 3/42000 450.0 10.950 10.950 10.950 -0.000 -0.000 0.000 3/42000 460.0 10.682 10.682 10.682 -0.000 -0.000 0.000 3/42000 470.0 10.428 10.428 10.428 -0.000 -0.000 0.000 3/42000 480.0 10.186 10.186 10.186 -0.000 -0.000 0.000 3/42000 490.0 9.955 9.955 9.955 -0.000 -0.000 0.000 3/42000 500.0 9.735 9.735 9.735 -0.000 -0.000 0.000 3/42000 510.0 9.524 9.524 9.524 -0.000 -0.000 0.000 3/42000 520.0 9.323 9.323 9.323 -0.000 -0.000 0.000 3/42000 530.0 9.131 9.131 9.131 -0.000 -0.000 0.000 3/42000 540.0 8.946 8.946 8.946 -0.000 -0.000 0.000 3/42000 550.0 8.769 8.769 8.769 -0.000 -0.000 0.000 3/42000 560.0 8.600 8.600 8.600 -0.000 -0.000 0.000 3/42000 570.0 8.436 8.436 8.436 -0.000 -0.000 0.000 3/42000 580.0 8.279 8.279 8.279 -0.000 -0.000 0.000 3/42000 590.0 8.128 8.128 8.128 -0.000 -0.000 0.000 3/42000 600.0 7.983 7.983 7.983 -0.000 -0.000 0.000 3/42000 610.0 7.843 7.843 7.843 -0.000 -0.000 0.000 3/42000 620.0 7.707 7.707 7.707 -0.000 -0.000 0.000 3/42000 630.0 7.577 7.577 7.577 -0.000 -0.000 0.000 3/42000 640.0 7.451 7.451 7.451 -0.000 -0.000 0.000 3/42000 650.0 7.329 7.329 7.329 -0.000 -0.000 0.000 3/42000 660.0 7.211 7.211 7.211 -0.000 -0.000 0.000 3/42000 670.0 7.097 7.097 7.097 -0.000 -0.000 0.000 3/42000 680.0 6.987 6.987 6.987 -0.000 -0.000 0.000 3/42000 690.0 6.880 6.880 6.880 -0.000 -0.000 0.000 3/42000 700.0 6.776 6.776 6.776 -0.000 -0.000 0.000 3/42000 710.0 6.676 6.676 6.676 -0.000 -0.000 0.000 3/42000 720.0 6.578 6.578 6.578 -0.000 -0.000 0.000 3/42000 730.0 6.484 6.484 6.484 -0.000 -0.000 0.000 3/42000 740.0 6.392 6.392 6.392 -0.000 -0.000 0.000 3/42000 750.0 6.303 6.303 6.303 -0.000 -0.000 0.000 3/42000 760.0 6.216 6.216 6.216 -0.000 -0.000 0.000 3/42000 770.0 6.132 6.132 6.132 -0.000 -0.000 0.000 3/42000 780.0 6.050 6.050 6.050 -0.000 -0.000 0.000 3/42000 790.0 5.970 5.970 5.970 -0.000 -0.000 0.000 3/42000 800.0 5.892 5.892 5.892 -0.000 -0.000 0.000 3/42000 810.0 5.817 5.817 5.817 -0.000 -0.000 0.000 3/42000 820.0 5.743 5.743 5.743 -0.000 -0.000 0.000 3/42000 830.0 5.671 5.671 5.671 -0.000 -0.000 0.000 3/42000 840.0 5.601 5.601 5.601 -0.000 -0.000 0.000 3/42000 850.0 5.533 5.533 5.533 -0.000 -0.000 0.000 3/42000 860.0 5.466 5.466 5.466 -0.000 -0.000 0.000 3/42000 870.0 5.401 5.401 5.401 -0.000 -0.000 0.000 3/42000 880.0 5.338 5.338 5.338 -0.000 -0.000 0.000 3/42000 890.0 5.276 5.276 5.276 -0.000 -0.000 0.000 3/42000 900.0 5.215 5.215 5.215 -0.000 -0.000 0.000 3/42000 910.0 5.156 5.156 5.156 -0.000 -0.000 0.000 3/42000 920.0 5.098 5.098 5.098 -0.000 -0.000 0.000 3/42000 930.0 5.042 5.042 5.042 -0.000 -0.000 0.000 3/42000 940.0 4.987 4.987 4.987 -0.000 -0.000 0.000 3/42000 950.0 4.933 4.933 4.933 -0.000 -0.000 0.000 3/42000 960.0 4.880 4.880 4.880 -0.000 -0.000 0.000 3/42000 970.0 4.828 4.828 4.828 -0.000 -0.000 0.000 3/42000 980.0 4.778 4.778 4.778 -0.000 -0.000 0.000 3/42000 990.0 4.728 4.728 4.728 -0.000 -0.000 0.000 3/42000 1000.0 4.680 4.680 4.680 -0.000 -0.000 0.000 3/42000 Thermal conductivity related properties were written into "kappa-m101010.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 21:02:50]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|