----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 18:30: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: [2 2 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pa-3 (205) Number of symmetry operations in supercell: 192 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 5.318512559999999 0.000000000000000 0.000000000000000 b 0.000000000000000 5.318512559999999 0.000000000000000 c 0.000000000000000 0.000000000000000 5.318512559999999 Atomic positions (fractional): *1 S 0.88176571948806 0.61823428051194 0.38176571948806 32.065 2 S 0.61823428051194 0.38176571948806 0.88176571948806 32.065 3 S 0.38176571948806 0.88176571948806 0.61823428051194 32.065 4 S 0.11823428051194 0.11823428051194 0.11823428051194 32.065 5 S 0.88176571948806 0.88176571948806 0.88176571948806 32.065 6 S 0.61823428051194 0.11823428051194 0.38176571948806 32.065 7 S 0.38176571948806 0.61823428051194 0.11823428051194 32.065 8 S 0.11823428051194 0.38176571948806 0.61823428051194 32.065 *9 Fe 0.50000000000000 0.00000000000000 0.00000000000000 55.845 10 Fe 0.00000000000000 0.00000000000000 0.50000000000000 55.845 11 Fe 0.00000000000000 0.50000000000000 0.00000000000000 55.845 12 Fe 0.50000000000000 0.50000000000000 0.50000000000000 55.845 -------------------------------- unit cell --------------------------------- Lattice vectors: a 5.318512559999999 0.000000000000000 0.000000000000000 b 0.000000000000000 5.318512559999999 0.000000000000000 c 0.000000000000000 0.000000000000000 5.318512559999999 Atomic positions (fractional): *1 S 0.88176571948806 0.61823428051194 0.38176571948806 32.065 > 1 2 S 0.61823428051194 0.38176571948806 0.88176571948806 32.065 > 2 3 S 0.38176571948806 0.88176571948806 0.61823428051194 32.065 > 3 4 S 0.11823428051194 0.11823428051194 0.11823428051194 32.065 > 4 5 S 0.88176571948806 0.88176571948806 0.88176571948806 32.065 > 5 6 S 0.61823428051194 0.11823428051194 0.38176571948806 32.065 > 6 7 S 0.38176571948806 0.61823428051194 0.11823428051194 32.065 > 7 8 S 0.11823428051194 0.38176571948806 0.61823428051194 32.065 > 8 *9 Fe 0.50000000000000 0.00000000000000 0.00000000000000 55.845 > 9 10 Fe 0.00000000000000 0.00000000000000 0.50000000000000 55.845 > 10 11 Fe 0.00000000000000 0.50000000000000 0.00000000000000 55.845 > 11 12 Fe 0.50000000000000 0.50000000000000 0.50000000000000 55.845 > 12 -------------------------------- super cell -------------------------------- Lattice vectors: a 10.637025119999999 0.000000000000000 0.000000000000000 b 0.000000000000000 10.637025119999999 0.000000000000000 c 0.000000000000000 0.000000000000000 10.637025119999999 Atomic positions (fractional): *1 S 0.44088285974403 0.30911714025597 0.19088285974403 32.065 > 1 2 S 0.94088285974403 0.30911714025597 0.19088285974403 32.065 > 1 3 S 0.44088285974403 0.80911714025597 0.19088285974403 32.065 > 1 4 S 0.94088285974403 0.80911714025597 0.19088285974403 32.065 > 1 5 S 0.44088285974403 0.30911714025597 0.69088285974403 32.065 > 1 6 S 0.94088285974403 0.30911714025597 0.69088285974403 32.065 > 1 7 S 0.44088285974403 0.80911714025597 0.69088285974403 32.065 > 1 8 S 0.94088285974403 0.80911714025597 0.69088285974403 32.065 > 1 9 S 0.30911714025597 0.19088285974403 0.44088285974403 32.065 > 2 10 S 0.80911714025597 0.19088285974403 0.44088285974403 32.065 > 2 11 S 0.30911714025597 0.69088285974403 0.44088285974403 32.065 > 2 12 S 0.80911714025597 0.69088285974403 0.44088285974403 32.065 > 2 13 S 0.30911714025597 0.19088285974403 0.94088285974403 32.065 > 2 14 S 0.80911714025597 0.19088285974403 0.94088285974403 32.065 > 2 15 S 0.30911714025597 0.69088285974403 0.94088285974403 32.065 > 2 16 S 0.80911714025597 0.69088285974403 0.94088285974403 32.065 > 2 17 S 0.19088285974403 0.44088285974403 0.30911714025597 32.065 > 3 18 S 0.69088285974403 0.44088285974403 0.30911714025597 32.065 > 3 19 S 0.19088285974403 0.94088285974403 0.30911714025597 32.065 > 3 20 S 0.69088285974403 0.94088285974403 0.30911714025597 32.065 > 3 21 S 0.19088285974403 0.44088285974403 0.80911714025597 32.065 > 3 22 S 0.69088285974403 0.44088285974403 0.80911714025597 32.065 > 3 23 S 0.19088285974403 0.94088285974403 0.80911714025597 32.065 > 3 24 S 0.69088285974403 0.94088285974403 0.80911714025597 32.065 > 3 25 S 0.05911714025597 0.05911714025597 0.05911714025597 32.065 > 4 26 S 0.55911714025597 0.05911714025597 0.05911714025597 32.065 > 4 27 S 0.05911714025597 0.55911714025597 0.05911714025597 32.065 > 4 28 S 0.55911714025597 0.55911714025597 0.05911714025597 32.065 > 4 29 S 0.05911714025597 0.05911714025597 0.55911714025597 32.065 > 4 30 S 0.55911714025597 0.05911714025597 0.55911714025597 32.065 > 4 31 S 0.05911714025597 0.55911714025597 0.55911714025597 32.065 > 4 32 S 0.55911714025597 0.55911714025597 0.55911714025597 32.065 > 4 33 S 0.44088285974403 0.44088285974403 0.44088285974403 32.065 > 5 34 S 0.94088285974403 0.44088285974403 0.44088285974403 32.065 > 5 35 S 0.44088285974403 0.94088285974403 0.44088285974403 32.065 > 5 36 S 0.94088285974403 0.94088285974403 0.44088285974403 32.065 > 5 37 S 0.44088285974403 0.44088285974403 0.94088285974403 32.065 > 5 38 S 0.94088285974403 0.44088285974403 0.94088285974403 32.065 > 5 39 S 0.44088285974403 0.94088285974403 0.94088285974403 32.065 > 5 40 S 0.94088285974403 0.94088285974403 0.94088285974403 32.065 > 5 41 S 0.30911714025597 0.05911714025597 0.19088285974403 32.065 > 6 42 S 0.80911714025597 0.05911714025597 0.19088285974403 32.065 > 6 43 S 0.30911714025597 0.55911714025597 0.19088285974403 32.065 > 6 44 S 0.80911714025597 0.55911714025597 0.19088285974403 32.065 > 6 45 S 0.30911714025597 0.05911714025597 0.69088285974403 32.065 > 6 46 S 0.80911714025597 0.05911714025597 0.69088285974403 32.065 > 6 47 S 0.30911714025597 0.55911714025597 0.69088285974403 32.065 > 6 48 S 0.80911714025597 0.55911714025597 0.69088285974403 32.065 > 6 49 S 0.19088285974403 0.30911714025597 0.05911714025597 32.065 > 7 50 S 0.69088285974403 0.30911714025597 0.05911714025597 32.065 > 7 51 S 0.19088285974403 0.80911714025597 0.05911714025597 32.065 > 7 52 S 0.69088285974403 0.80911714025597 0.05911714025597 32.065 > 7 53 S 0.19088285974403 0.30911714025597 0.55911714025597 32.065 > 7 54 S 0.69088285974403 0.30911714025597 0.55911714025597 32.065 > 7 55 S 0.19088285974403 0.80911714025597 0.55911714025597 32.065 > 7 56 S 0.69088285974403 0.80911714025597 0.55911714025597 32.065 > 7 57 S 0.05911714025597 0.19088285974403 0.30911714025597 32.065 > 8 58 S 0.55911714025597 0.19088285974403 0.30911714025597 32.065 > 8 59 S 0.05911714025597 0.69088285974403 0.30911714025597 32.065 > 8 60 S 0.55911714025597 0.69088285974403 0.30911714025597 32.065 > 8 61 S 0.05911714025597 0.19088285974403 0.80911714025597 32.065 > 8 62 S 0.55911714025597 0.19088285974403 0.80911714025597 32.065 > 8 63 S 0.05911714025597 0.69088285974403 0.80911714025597 32.065 > 8 64 S 0.55911714025597 0.69088285974403 0.80911714025597 32.065 > 8 *65 Fe 0.25000000000000 0.00000000000000 0.00000000000000 55.845 > 9 66 Fe 0.75000000000000 0.00000000000000 0.00000000000000 55.845 > 9 67 Fe 0.25000000000000 0.50000000000000 0.00000000000000 55.845 > 9 68 Fe 0.75000000000000 0.50000000000000 0.00000000000000 55.845 > 9 69 Fe 0.25000000000000 0.00000000000000 0.50000000000000 55.845 > 9 70 Fe 0.75000000000000 0.00000000000000 0.50000000000000 55.845 > 9 71 Fe 0.25000000000000 0.50000000000000 0.50000000000000 55.845 > 9 72 Fe 0.75000000000000 0.50000000000000 0.50000000000000 55.845 > 9 73 Fe 0.00000000000000 0.00000000000000 0.25000000000000 55.845 > 10 74 Fe 0.50000000000000 0.00000000000000 0.25000000000000 55.845 > 10 75 Fe 0.00000000000000 0.50000000000000 0.25000000000000 55.845 > 10 76 Fe 0.50000000000000 0.50000000000000 0.25000000000000 55.845 > 10 77 Fe 0.00000000000000 0.00000000000000 0.75000000000000 55.845 > 10 78 Fe 0.50000000000000 0.00000000000000 0.75000000000000 55.845 > 10 79 Fe 0.00000000000000 0.50000000000000 0.75000000000000 55.845 > 10 80 Fe 0.50000000000000 0.50000000000000 0.75000000000000 55.845 > 10 81 Fe 0.00000000000000 0.25000000000000 0.00000000000000 55.845 > 11 82 Fe 0.50000000000000 0.25000000000000 0.00000000000000 55.845 > 11 83 Fe 0.00000000000000 0.75000000000000 0.00000000000000 55.845 > 11 84 Fe 0.50000000000000 0.75000000000000 0.00000000000000 55.845 > 11 85 Fe 0.00000000000000 0.25000000000000 0.50000000000000 55.845 > 11 86 Fe 0.50000000000000 0.25000000000000 0.50000000000000 55.845 > 11 87 Fe 0.00000000000000 0.75000000000000 0.50000000000000 55.845 > 11 88 Fe 0.50000000000000 0.75000000000000 0.50000000000000 55.845 > 11 89 Fe 0.25000000000000 0.25000000000000 0.25000000000000 55.845 > 12 90 Fe 0.75000000000000 0.25000000000000 0.25000000000000 55.845 > 12 91 Fe 0.25000000000000 0.75000000000000 0.25000000000000 55.845 > 12 92 Fe 0.75000000000000 0.75000000000000 0.25000000000000 55.845 > 12 93 Fe 0.25000000000000 0.25000000000000 0.75000000000000 55.845 > 12 94 Fe 0.75000000000000 0.25000000000000 0.75000000000000 55.845 > 12 95 Fe 0.25000000000000 0.75000000000000 0.75000000000000 55.845 > 12 96 Fe 0.75000000000000 0.75000000000000 0.75000000000000 55.845 > 12 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 20.9275513 0.0000000 0.0000000 0.0000000 20.9275513 0.0000000 0.0000000 0.0000000 20.9275513 -------------------------- Born effective charges -------------------------- 1 S 2.6271555 -0.1002041 -0.0089271 0.0089271 2.6271555 -0.1002041 0.1002041 0.0089271 2.6271555 2 S 2.6271555 -0.1002041 0.0089271 0.0089271 2.6271555 0.1002041 -0.1002041 -0.0089271 2.6271555 3 S 2.6271555 0.1002041 0.0089271 -0.0089271 2.6271555 -0.1002041 -0.1002041 0.0089271 2.6271555 4 S 2.6271555 0.1002041 -0.0089271 -0.0089271 2.6271555 0.1002041 0.1002041 -0.0089271 2.6271555 5 S 2.6271555 0.1002041 -0.0089271 -0.0089271 2.6271555 0.1002041 0.1002041 -0.0089271 2.6271555 6 S 2.6271555 0.1002041 0.0089271 -0.0089271 2.6271555 -0.1002041 -0.1002041 0.0089271 2.6271555 7 S 2.6271555 -0.1002041 0.0089271 0.0089271 2.6271555 0.1002041 -0.1002041 -0.0089271 2.6271555 8 S 2.6271555 -0.1002041 -0.0089271 0.0089271 2.6271555 -0.1002041 0.1002041 0.0089271 2.6271555 9 Fe -5.2543110 0.4643544 -0.0071481 0.0071481 -5.2543110 0.4643544 -0.4643544 0.0071481 -5.2543110 10 Fe -5.2543110 0.4643544 0.0071481 0.0071481 -5.2543110 -0.4643544 0.4643544 -0.0071481 -5.2543110 11 Fe -5.2543110 -0.4643544 0.0071481 -0.0071481 -5.2543110 0.4643544 0.4643544 0.0071481 -5.2543110 12 Fe -5.2543110 -0.4643544 -0.0071481 -0.0071481 -5.2543110 -0.4643544 -0.4643544 -0.0071481 -5.2543110 ---------------------------------------------------------------------------- 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: 288/288 Permutation basis: 5544/5544 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 234 Number of blocks in projector: 234 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: 148 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 86 Use standard eigh solver. Tree of FC basis block matrices: - (234, 228), data: False |-- (86, 84), data: True |-- (148, 144), data: True ----- Solver_atoms: 1 -- 96 / 96 Time (Solver_compr_matrix_reshape): 0.002 Solver_block: 80 / 80 - Time: 0.151 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.155 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 288/288 Permutation basis: 5544/5544 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 234 Number of blocks in projector: 234 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: 148 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 86 Use standard eigh solver. Tree of FC basis block matrices: - (234, 228), data: False |-- (86, 84), data: True |-- (148, 144), data: True Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 18:30:22]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:30:23]------------------------- 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): [2 2 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pa-3 (205) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 5.318512559999999 0.000000000000000 0.000000000000000 b 0.000000000000000 5.318512559999999 0.000000000000000 c 0.000000000000000 0.000000000000000 5.318512559999999 Atomic positions (fractional): 1 S 0.88176571948806 0.61823428051194 0.38176571948806 32.065 2 S 0.61823428051194 0.38176571948806 0.88176571948806 32.065 3 S 0.38176571948806 0.88176571948806 0.61823428051194 32.065 4 S 0.11823428051194 0.11823428051194 0.11823428051194 32.065 5 S 0.88176571948806 0.88176571948806 0.88176571948806 32.065 6 S 0.61823428051194 0.11823428051194 0.38176571948806 32.065 7 S 0.38176571948806 0.61823428051194 0.11823428051194 32.065 8 S 0.11823428051194 0.38176571948806 0.61823428051194 32.065 9 Fe 0.50000000000000 0.00000000000000 0.00000000000000 55.845 10 Fe 0.00000000000000 0.00000000000000 0.50000000000000 55.845 11 Fe 0.00000000000000 0.50000000000000 0.00000000000000 55.845 12 Fe 0.50000000000000 0.50000000000000 0.50000000000000 55.845 -------------------------------- supercell --------------------------------- Lattice vectors: a 10.637025119999999 0.000000000000000 0.000000000000000 b 0.000000000000000 10.637025119999999 0.000000000000000 c 0.000000000000000 0.000000000000000 10.637025119999999 Atomic positions (fractional): 1 S 0.44088285974403 0.30911714025597 0.19088285974403 32.065 > 1 2 S 0.94088285974403 0.30911714025597 0.19088285974403 32.065 > 1 3 S 0.44088285974403 0.80911714025597 0.19088285974403 32.065 > 1 4 S 0.94088285974403 0.80911714025597 0.19088285974403 32.065 > 1 5 S 0.44088285974403 0.30911714025597 0.69088285974403 32.065 > 1 6 S 0.94088285974403 0.30911714025597 0.69088285974403 32.065 > 1 7 S 0.44088285974403 0.80911714025597 0.69088285974403 32.065 > 1 8 S 0.94088285974403 0.80911714025597 0.69088285974403 32.065 > 1 9 S 0.30911714025597 0.19088285974403 0.44088285974403 32.065 > 9 10 S 0.80911714025597 0.19088285974403 0.44088285974403 32.065 > 9 11 S 0.30911714025597 0.69088285974403 0.44088285974403 32.065 > 9 12 S 0.80911714025597 0.69088285974403 0.44088285974403 32.065 > 9 13 S 0.30911714025597 0.19088285974403 0.94088285974403 32.065 > 9 14 S 0.80911714025597 0.19088285974403 0.94088285974403 32.065 > 9 15 S 0.30911714025597 0.69088285974403 0.94088285974403 32.065 > 9 16 S 0.80911714025597 0.69088285974403 0.94088285974403 32.065 > 9 17 S 0.19088285974403 0.44088285974403 0.30911714025597 32.065 > 17 18 S 0.69088285974403 0.44088285974403 0.30911714025597 32.065 > 17 19 S 0.19088285974403 0.94088285974403 0.30911714025597 32.065 > 17 20 S 0.69088285974403 0.94088285974403 0.30911714025597 32.065 > 17 21 S 0.19088285974403 0.44088285974403 0.80911714025597 32.065 > 17 22 S 0.69088285974403 0.44088285974403 0.80911714025597 32.065 > 17 23 S 0.19088285974403 0.94088285974403 0.80911714025597 32.065 > 17 24 S 0.69088285974403 0.94088285974403 0.80911714025597 32.065 > 17 25 S 0.05911714025597 0.05911714025597 0.05911714025597 32.065 > 25 26 S 0.55911714025597 0.05911714025597 0.05911714025597 32.065 > 25 27 S 0.05911714025597 0.55911714025597 0.05911714025597 32.065 > 25 28 S 0.55911714025597 0.55911714025597 0.05911714025597 32.065 > 25 29 S 0.05911714025597 0.05911714025597 0.55911714025597 32.065 > 25 30 S 0.55911714025597 0.05911714025597 0.55911714025597 32.065 > 25 31 S 0.05911714025597 0.55911714025597 0.55911714025597 32.065 > 25 32 S 0.55911714025597 0.55911714025597 0.55911714025597 32.065 > 25 33 S 0.44088285974403 0.44088285974403 0.44088285974403 32.065 > 33 34 S 0.94088285974403 0.44088285974403 0.44088285974403 32.065 > 33 35 S 0.44088285974403 0.94088285974403 0.44088285974403 32.065 > 33 36 S 0.94088285974403 0.94088285974403 0.44088285974403 32.065 > 33 37 S 0.44088285974403 0.44088285974403 0.94088285974403 32.065 > 33 38 S 0.94088285974403 0.44088285974403 0.94088285974403 32.065 > 33 39 S 0.44088285974403 0.94088285974403 0.94088285974403 32.065 > 33 40 S 0.94088285974403 0.94088285974403 0.94088285974403 32.065 > 33 41 S 0.30911714025597 0.05911714025597 0.19088285974403 32.065 > 41 42 S 0.80911714025597 0.05911714025597 0.19088285974403 32.065 > 41 43 S 0.30911714025597 0.55911714025597 0.19088285974403 32.065 > 41 44 S 0.80911714025597 0.55911714025597 0.19088285974403 32.065 > 41 45 S 0.30911714025597 0.05911714025597 0.69088285974403 32.065 > 41 46 S 0.80911714025597 0.05911714025597 0.69088285974403 32.065 > 41 47 S 0.30911714025597 0.55911714025597 0.69088285974403 32.065 > 41 48 S 0.80911714025597 0.55911714025597 0.69088285974403 32.065 > 41 49 S 0.19088285974403 0.30911714025597 0.05911714025597 32.065 > 49 50 S 0.69088285974403 0.30911714025597 0.05911714025597 32.065 > 49 51 S 0.19088285974403 0.80911714025597 0.05911714025597 32.065 > 49 52 S 0.69088285974403 0.80911714025597 0.05911714025597 32.065 > 49 53 S 0.19088285974403 0.30911714025597 0.55911714025597 32.065 > 49 54 S 0.69088285974403 0.30911714025597 0.55911714025597 32.065 > 49 55 S 0.19088285974403 0.80911714025597 0.55911714025597 32.065 > 49 56 S 0.69088285974403 0.80911714025597 0.55911714025597 32.065 > 49 57 S 0.05911714025597 0.19088285974403 0.30911714025597 32.065 > 57 58 S 0.55911714025597 0.19088285974403 0.30911714025597 32.065 > 57 59 S 0.05911714025597 0.69088285974403 0.30911714025597 32.065 > 57 60 S 0.55911714025597 0.69088285974403 0.30911714025597 32.065 > 57 61 S 0.05911714025597 0.19088285974403 0.80911714025597 32.065 > 57 62 S 0.55911714025597 0.19088285974403 0.80911714025597 32.065 > 57 63 S 0.05911714025597 0.69088285974403 0.80911714025597 32.065 > 57 64 S 0.55911714025597 0.69088285974403 0.80911714025597 32.065 > 57 65 Fe 0.25000000000000 0.00000000000000 0.00000000000000 55.845 > 65 66 Fe 0.75000000000000 0.00000000000000 0.00000000000000 55.845 > 65 67 Fe 0.25000000000000 0.50000000000000 0.00000000000000 55.845 > 65 68 Fe 0.75000000000000 0.50000000000000 0.00000000000000 55.845 > 65 69 Fe 0.25000000000000 0.00000000000000 0.50000000000000 55.845 > 65 70 Fe 0.75000000000000 0.00000000000000 0.50000000000000 55.845 > 65 71 Fe 0.25000000000000 0.50000000000000 0.50000000000000 55.845 > 65 72 Fe 0.75000000000000 0.50000000000000 0.50000000000000 55.845 > 65 73 Fe 0.00000000000000 0.00000000000000 0.25000000000000 55.845 > 73 74 Fe 0.50000000000000 0.00000000000000 0.25000000000000 55.845 > 73 75 Fe 0.00000000000000 0.50000000000000 0.25000000000000 55.845 > 73 76 Fe 0.50000000000000 0.50000000000000 0.25000000000000 55.845 > 73 77 Fe 0.00000000000000 0.00000000000000 0.75000000000000 55.845 > 73 78 Fe 0.50000000000000 0.00000000000000 0.75000000000000 55.845 > 73 79 Fe 0.00000000000000 0.50000000000000 0.75000000000000 55.845 > 73 80 Fe 0.50000000000000 0.50000000000000 0.75000000000000 55.845 > 73 81 Fe 0.00000000000000 0.25000000000000 0.00000000000000 55.845 > 81 82 Fe 0.50000000000000 0.25000000000000 0.00000000000000 55.845 > 81 83 Fe 0.00000000000000 0.75000000000000 0.00000000000000 55.845 > 81 84 Fe 0.50000000000000 0.75000000000000 0.00000000000000 55.845 > 81 85 Fe 0.00000000000000 0.25000000000000 0.50000000000000 55.845 > 81 86 Fe 0.50000000000000 0.25000000000000 0.50000000000000 55.845 > 81 87 Fe 0.00000000000000 0.75000000000000 0.50000000000000 55.845 > 81 88 Fe 0.50000000000000 0.75000000000000 0.50000000000000 55.845 > 81 89 Fe 0.25000000000000 0.25000000000000 0.25000000000000 55.845 > 89 90 Fe 0.75000000000000 0.25000000000000 0.25000000000000 55.845 > 89 91 Fe 0.25000000000000 0.75000000000000 0.25000000000000 55.845 > 89 92 Fe 0.75000000000000 0.75000000000000 0.25000000000000 55.845 > 89 93 Fe 0.25000000000000 0.25000000000000 0.75000000000000 55.845 > 89 94 Fe 0.75000000000000 0.25000000000000 0.75000000000000 55.845 > 89 95 Fe 0.25000000000000 0.75000000000000 0.75000000000000 55.845 > 89 96 Fe 0.75000000000000 0.75000000000000 0.75000000000000 55.845 > 89 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 20.9275513 0.0000000 0.0000000 0.0000000 20.9275513 0.0000000 0.0000000 0.0000000 20.9275513 -------------------------- Born effective charges -------------------------- 1 S 2.6271555 -0.1002041 -0.0089271 0.0089271 2.6271555 -0.1002041 0.1002041 0.0089271 2.6271555 2 S 2.6271555 -0.1002041 0.0089271 0.0089271 2.6271555 0.1002041 -0.1002041 -0.0089271 2.6271555 3 S 2.6271555 0.1002041 0.0089271 -0.0089271 2.6271555 -0.1002041 -0.1002041 0.0089271 2.6271555 4 S 2.6271555 0.1002041 -0.0089271 -0.0089271 2.6271555 0.1002041 0.1002041 -0.0089271 2.6271555 5 S 2.6271555 0.1002041 -0.0089271 -0.0089271 2.6271555 0.1002041 0.1002041 -0.0089271 2.6271555 6 S 2.6271555 0.1002041 0.0089271 -0.0089271 2.6271555 -0.1002041 -0.1002041 0.0089271 2.6271555 7 S 2.6271555 -0.1002041 0.0089271 0.0089271 2.6271555 0.1002041 -0.1002041 -0.0089271 2.6271555 8 S 2.6271555 -0.1002041 -0.0089271 0.0089271 2.6271555 -0.1002041 0.1002041 0.0089271 2.6271555 9 Fe -5.2543110 0.4643544 -0.0071481 0.0071481 -5.2543110 0.4643544 -0.4643544 0.0071481 -5.2543110 10 Fe -5.2543110 0.4643544 0.0071481 0.0071481 -5.2543110 -0.4643544 0.4643544 -0.0071481 -5.2543110 11 Fe -5.2543110 -0.4643544 0.0071481 -0.0071481 -5.2543110 0.4643544 0.4643544 0.0071481 -5.2543110 12 Fe -5.2543110 -0.4643544 -0.0071481 -0.0071481 -5.2543110 -0.4643544 -0.4643544 -0.0071481 -5.2543110 ---------------------------------------------------------------------------- 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[ 65, 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.00000010 (yyy) -0.00000010 (yyy) -0.00000010 (yyy) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000038 (yx) 0.00000038 (yx) 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 18:30:27]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:30:27]------------------------- 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): [2 2 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pa-3 (205) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 5.318512559999999 0.000000000000000 0.000000000000000 b 0.000000000000000 5.318512559999999 0.000000000000000 c 0.000000000000000 0.000000000000000 5.318512559999999 Atomic positions (fractional): 1 S 0.88176571948806 0.61823428051194 0.38176571948806 32.065 2 S 0.61823428051194 0.38176571948806 0.88176571948806 32.065 3 S 0.38176571948806 0.88176571948806 0.61823428051194 32.065 4 S 0.11823428051194 0.11823428051194 0.11823428051194 32.065 5 S 0.88176571948806 0.88176571948806 0.88176571948806 32.065 6 S 0.61823428051194 0.11823428051194 0.38176571948806 32.065 7 S 0.38176571948806 0.61823428051194 0.11823428051194 32.065 8 S 0.11823428051194 0.38176571948806 0.61823428051194 32.065 9 Fe 0.50000000000000 0.00000000000000 0.00000000000000 55.845 10 Fe 0.00000000000000 0.00000000000000 0.50000000000000 55.845 11 Fe 0.00000000000000 0.50000000000000 0.00000000000000 55.845 12 Fe 0.50000000000000 0.50000000000000 0.50000000000000 55.845 -------------------------------- supercell --------------------------------- Lattice vectors: a 10.637025119999999 0.000000000000000 0.000000000000000 b 0.000000000000000 10.637025119999999 0.000000000000000 c 0.000000000000000 0.000000000000000 10.637025119999999 Atomic positions (fractional): 1 S 0.44088285974403 0.30911714025597 0.19088285974403 32.065 > 1 2 S 0.94088285974403 0.30911714025597 0.19088285974403 32.065 > 1 3 S 0.44088285974403 0.80911714025597 0.19088285974403 32.065 > 1 4 S 0.94088285974403 0.80911714025597 0.19088285974403 32.065 > 1 5 S 0.44088285974403 0.30911714025597 0.69088285974403 32.065 > 1 6 S 0.94088285974403 0.30911714025597 0.69088285974403 32.065 > 1 7 S 0.44088285974403 0.80911714025597 0.69088285974403 32.065 > 1 8 S 0.94088285974403 0.80911714025597 0.69088285974403 32.065 > 1 9 S 0.30911714025597 0.19088285974403 0.44088285974403 32.065 > 9 10 S 0.80911714025597 0.19088285974403 0.44088285974403 32.065 > 9 11 S 0.30911714025597 0.69088285974403 0.44088285974403 32.065 > 9 12 S 0.80911714025597 0.69088285974403 0.44088285974403 32.065 > 9 13 S 0.30911714025597 0.19088285974403 0.94088285974403 32.065 > 9 14 S 0.80911714025597 0.19088285974403 0.94088285974403 32.065 > 9 15 S 0.30911714025597 0.69088285974403 0.94088285974403 32.065 > 9 16 S 0.80911714025597 0.69088285974403 0.94088285974403 32.065 > 9 17 S 0.19088285974403 0.44088285974403 0.30911714025597 32.065 > 17 18 S 0.69088285974403 0.44088285974403 0.30911714025597 32.065 > 17 19 S 0.19088285974403 0.94088285974403 0.30911714025597 32.065 > 17 20 S 0.69088285974403 0.94088285974403 0.30911714025597 32.065 > 17 21 S 0.19088285974403 0.44088285974403 0.80911714025597 32.065 > 17 22 S 0.69088285974403 0.44088285974403 0.80911714025597 32.065 > 17 23 S 0.19088285974403 0.94088285974403 0.80911714025597 32.065 > 17 24 S 0.69088285974403 0.94088285974403 0.80911714025597 32.065 > 17 25 S 0.05911714025597 0.05911714025597 0.05911714025597 32.065 > 25 26 S 0.55911714025597 0.05911714025597 0.05911714025597 32.065 > 25 27 S 0.05911714025597 0.55911714025597 0.05911714025597 32.065 > 25 28 S 0.55911714025597 0.55911714025597 0.05911714025597 32.065 > 25 29 S 0.05911714025597 0.05911714025597 0.55911714025597 32.065 > 25 30 S 0.55911714025597 0.05911714025597 0.55911714025597 32.065 > 25 31 S 0.05911714025597 0.55911714025597 0.55911714025597 32.065 > 25 32 S 0.55911714025597 0.55911714025597 0.55911714025597 32.065 > 25 33 S 0.44088285974403 0.44088285974403 0.44088285974403 32.065 > 33 34 S 0.94088285974403 0.44088285974403 0.44088285974403 32.065 > 33 35 S 0.44088285974403 0.94088285974403 0.44088285974403 32.065 > 33 36 S 0.94088285974403 0.94088285974403 0.44088285974403 32.065 > 33 37 S 0.44088285974403 0.44088285974403 0.94088285974403 32.065 > 33 38 S 0.94088285974403 0.44088285974403 0.94088285974403 32.065 > 33 39 S 0.44088285974403 0.94088285974403 0.94088285974403 32.065 > 33 40 S 0.94088285974403 0.94088285974403 0.94088285974403 32.065 > 33 41 S 0.30911714025597 0.05911714025597 0.19088285974403 32.065 > 41 42 S 0.80911714025597 0.05911714025597 0.19088285974403 32.065 > 41 43 S 0.30911714025597 0.55911714025597 0.19088285974403 32.065 > 41 44 S 0.80911714025597 0.55911714025597 0.19088285974403 32.065 > 41 45 S 0.30911714025597 0.05911714025597 0.69088285974403 32.065 > 41 46 S 0.80911714025597 0.05911714025597 0.69088285974403 32.065 > 41 47 S 0.30911714025597 0.55911714025597 0.69088285974403 32.065 > 41 48 S 0.80911714025597 0.55911714025597 0.69088285974403 32.065 > 41 49 S 0.19088285974403 0.30911714025597 0.05911714025597 32.065 > 49 50 S 0.69088285974403 0.30911714025597 0.05911714025597 32.065 > 49 51 S 0.19088285974403 0.80911714025597 0.05911714025597 32.065 > 49 52 S 0.69088285974403 0.80911714025597 0.05911714025597 32.065 > 49 53 S 0.19088285974403 0.30911714025597 0.55911714025597 32.065 > 49 54 S 0.69088285974403 0.30911714025597 0.55911714025597 32.065 > 49 55 S 0.19088285974403 0.80911714025597 0.55911714025597 32.065 > 49 56 S 0.69088285974403 0.80911714025597 0.55911714025597 32.065 > 49 57 S 0.05911714025597 0.19088285974403 0.30911714025597 32.065 > 57 58 S 0.55911714025597 0.19088285974403 0.30911714025597 32.065 > 57 59 S 0.05911714025597 0.69088285974403 0.30911714025597 32.065 > 57 60 S 0.55911714025597 0.69088285974403 0.30911714025597 32.065 > 57 61 S 0.05911714025597 0.19088285974403 0.80911714025597 32.065 > 57 62 S 0.55911714025597 0.19088285974403 0.80911714025597 32.065 > 57 63 S 0.05911714025597 0.69088285974403 0.80911714025597 32.065 > 57 64 S 0.55911714025597 0.69088285974403 0.80911714025597 32.065 > 57 65 Fe 0.25000000000000 0.00000000000000 0.00000000000000 55.845 > 65 66 Fe 0.75000000000000 0.00000000000000 0.00000000000000 55.845 > 65 67 Fe 0.25000000000000 0.50000000000000 0.00000000000000 55.845 > 65 68 Fe 0.75000000000000 0.50000000000000 0.00000000000000 55.845 > 65 69 Fe 0.25000000000000 0.00000000000000 0.50000000000000 55.845 > 65 70 Fe 0.75000000000000 0.00000000000000 0.50000000000000 55.845 > 65 71 Fe 0.25000000000000 0.50000000000000 0.50000000000000 55.845 > 65 72 Fe 0.75000000000000 0.50000000000000 0.50000000000000 55.845 > 65 73 Fe 0.00000000000000 0.00000000000000 0.25000000000000 55.845 > 73 74 Fe 0.50000000000000 0.00000000000000 0.25000000000000 55.845 > 73 75 Fe 0.00000000000000 0.50000000000000 0.25000000000000 55.845 > 73 76 Fe 0.50000000000000 0.50000000000000 0.25000000000000 55.845 > 73 77 Fe 0.00000000000000 0.00000000000000 0.75000000000000 55.845 > 73 78 Fe 0.50000000000000 0.00000000000000 0.75000000000000 55.845 > 73 79 Fe 0.00000000000000 0.50000000000000 0.75000000000000 55.845 > 73 80 Fe 0.50000000000000 0.50000000000000 0.75000000000000 55.845 > 73 81 Fe 0.00000000000000 0.25000000000000 0.00000000000000 55.845 > 81 82 Fe 0.50000000000000 0.25000000000000 0.00000000000000 55.845 > 81 83 Fe 0.00000000000000 0.75000000000000 0.00000000000000 55.845 > 81 84 Fe 0.50000000000000 0.75000000000000 0.00000000000000 55.845 > 81 85 Fe 0.00000000000000 0.25000000000000 0.50000000000000 55.845 > 81 86 Fe 0.50000000000000 0.25000000000000 0.50000000000000 55.845 > 81 87 Fe 0.00000000000000 0.75000000000000 0.50000000000000 55.845 > 81 88 Fe 0.50000000000000 0.75000000000000 0.50000000000000 55.845 > 81 89 Fe 0.25000000000000 0.25000000000000 0.25000000000000 55.845 > 89 90 Fe 0.75000000000000 0.25000000000000 0.25000000000000 55.845 > 89 91 Fe 0.25000000000000 0.75000000000000 0.25000000000000 55.845 > 89 92 Fe 0.75000000000000 0.75000000000000 0.25000000000000 55.845 > 89 93 Fe 0.25000000000000 0.25000000000000 0.75000000000000 55.845 > 89 94 Fe 0.75000000000000 0.25000000000000 0.75000000000000 55.845 > 89 95 Fe 0.25000000000000 0.75000000000000 0.75000000000000 55.845 > 89 96 Fe 0.75000000000000 0.75000000000000 0.75000000000000 55.845 > 89 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 20.9275513 0.0000000 0.0000000 0.0000000 20.9275513 0.0000000 0.0000000 0.0000000 20.9275513 -------------------------- Born effective charges -------------------------- 1 S 2.6271555 -0.1002041 -0.0089271 0.0089271 2.6271555 -0.1002041 0.1002041 0.0089271 2.6271555 2 S 2.6271555 -0.1002041 0.0089271 0.0089271 2.6271555 0.1002041 -0.1002041 -0.0089271 2.6271555 3 S 2.6271555 0.1002041 0.0089271 -0.0089271 2.6271555 -0.1002041 -0.1002041 0.0089271 2.6271555 4 S 2.6271555 0.1002041 -0.0089271 -0.0089271 2.6271555 0.1002041 0.1002041 -0.0089271 2.6271555 5 S 2.6271555 0.1002041 -0.0089271 -0.0089271 2.6271555 0.1002041 0.1002041 -0.0089271 2.6271555 6 S 2.6271555 0.1002041 0.0089271 -0.0089271 2.6271555 -0.1002041 -0.1002041 0.0089271 2.6271555 7 S 2.6271555 -0.1002041 0.0089271 0.0089271 2.6271555 0.1002041 -0.1002041 -0.0089271 2.6271555 8 S 2.6271555 -0.1002041 -0.0089271 0.0089271 2.6271555 -0.1002041 0.1002041 0.0089271 2.6271555 9 Fe -5.2543110 0.4643544 -0.0071481 0.0071481 -5.2543110 0.4643544 -0.4643544 0.0071481 -5.2543110 10 Fe -5.2543110 0.4643544 0.0071481 0.0071481 -5.2543110 -0.4643544 0.4643544 -0.0071481 -5.2543110 11 Fe -5.2543110 -0.4643544 0.0071481 -0.0071481 -5.2543110 0.4643544 0.4643544 0.0071481 -5.2543110 12 Fe -5.2543110 -0.4643544 -0.0071481 -0.0071481 -5.2543110 -0.4643544 -0.4643544 -0.0071481 -5.2543110 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000010 (yyy) -0.00000010 (yyy) -0.00000010 (yyy) Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) --------------------------- 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 ... [ 9 9 9 ] 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.78, Number of G-points: 305, Lambda: 0.37 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/45) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 45 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 6.658 ( 0.000 0.000 0.000) 0.000 6.691 ( 0.000 0.000 0.000) 0.000 6.691 ( 0.000 0.000 0.000) 0.000 6.691 ( 0.000 0.000 0.000) 0.000 7.509 ( 0.000 0.000 0.000) 0.000 7.509 ( 0.000 0.000 0.000) 0.000 9.204 ( 0.000 0.000 0.000) 0.000 9.204 ( 0.000 0.000 0.000) 0.000 9.204 ( 0.000 0.000 0.000) 0.000 10.256 ( 0.000 0.000 0.000) 0.000 10.482 ( 0.000 0.000 0.000) 0.000 10.482 ( 0.000 0.000 0.000) 0.000 10.695 ( 0.000 0.000 0.000) 0.000 10.695 ( 0.000 0.000 0.000) 0.000 10.695 ( 0.000 0.000 0.000) 0.000 10.898 ( 0.000 0.000 0.000) 0.000 10.898 ( 0.000 0.000 0.000) 0.000 10.898 ( 0.000 0.000 0.000) 0.000 11.334 ( 0.000 0.000 0.000) 0.000 11.436 ( 0.000 0.000 0.000) 0.000 11.436 ( 0.000 0.000 0.000) 0.000 11.436 ( 0.000 0.000 0.000) 0.000 11.627 ( 0.000 0.000 0.000) 0.000 11.627 ( 0.000 0.000 0.000) 0.000 11.989 ( 0.000 0.000 0.000) 0.000 11.989 ( 0.000 0.000 0.000) 0.000 11.989 ( 0.000 0.000 0.000) 0.000 12.347 ( 0.000 0.000 0.000) 0.000 12.347 ( 0.000 0.000 0.000) 0.000 12.347 ( 0.000 0.000 0.000) 0.000 13.243 ( 0.000 0.000 0.000) 0.000 13.243 ( 0.000 0.000 0.000) 0.000 13.243 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/45) ======================= q-point: ( 0.11 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 125 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.004 ( 47.984 0.000 0.000) 47.984 1.021 ( 48.481 0.000 0.000) 48.481 1.938 ( 91.264 0.000 0.000) 91.264 6.546 ( -14.385 0.000 0.000) 14.385 6.652 ( -3.569 0.000 0.000) 3.569 6.781 ( 8.497 0.000 0.000) 8.497 6.863 ( 15.936 0.000 0.000) 15.936 7.438 ( -3.278 0.000 0.000) 3.278 7.583 ( 6.852 0.000 0.000) 6.852 9.162 ( -3.842 0.000 0.000) 3.842 9.174 ( -7.558 0.000 0.000) 7.558 9.181 ( -2.608 0.000 0.000) 2.608 10.205 ( -4.656 0.000 0.000) 4.656 10.388 ( -10.013 0.000 0.000) 10.013 10.538 ( 4.865 0.000 0.000) 4.865 10.614 ( -9.733 0.000 0.000) 9.733 10.653 ( -7.687 0.000 0.000) 7.687 10.778 ( -7.237 0.000 0.000) 7.237 10.800 ( 8.094 0.000 0.000) 8.094 10.824 ( -7.042 0.000 0.000) 7.042 11.054 ( 3.996 0.000 0.000) 3.996 11.109 ( -7.750 0.000 0.000) 7.750 11.236 ( -9.815 0.000 0.000) 9.815 11.309 ( -6.948 0.000 0.000) 6.948 11.465 ( -10.828 0.000 0.000) 10.828 11.578 ( -4.360 0.000 0.000) 4.360 11.771 ( 17.037 0.000 0.000) 17.037 12.084 ( 5.743 0.000 0.000) 5.743 12.202 ( 7.206 0.000 0.000) 7.206 12.220 ( 16.342 0.000 0.000) 16.342 12.326 ( -2.037 0.000 0.000) 2.037 12.366 ( 7.770 0.000 0.000) 7.770 13.181 ( -5.480 0.000 0.000) 5.480 13.218 ( -2.381 0.000 0.000) 2.381 13.280 ( 3.298 0.000 0.000) 3.298 14.102 ( -5.558 0.000 0.000) 5.558 ======================= Grid point 2 (3/45) ======================= q-point: ( 0.22 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 125 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.003 ( 47.610 0.000 0.000) 47.610 2.019 ( 46.917 0.000 0.000) 46.917 3.777 ( 83.564 0.000 0.000) 83.564 6.135 ( -23.713 0.000 0.000) 23.713 6.357 ( -24.873 0.000 0.000) 24.873 7.035 ( 15.466 0.000 0.000) 15.466 7.317 ( 26.181 0.000 0.000) 26.181 7.544 ( 13.824 0.000 0.000) 13.824 7.779 ( 11.446 0.000 0.000) 11.446 8.965 ( -12.041 0.000 0.000) 12.041 9.056 ( -5.829 0.000 0.000) 5.829 9.065 ( -9.489 0.000 0.000) 9.489 10.070 ( -8.124 0.000 0.000) 8.124 10.079 ( -17.691 0.000 0.000) 17.691 10.231 ( -24.510 0.000 0.000) 24.510 10.313 ( -20.136 0.000 0.000) 20.136 10.658 ( 5.874 0.000 0.000) 5.874 10.681 ( -6.122 0.000 0.000) 6.122 10.729 ( -0.333 0.000 0.000) 0.333 10.863 ( -5.984 0.000 0.000) 5.984 10.999 ( -1.836 0.000 0.000) 1.836 11.141 ( -1.090 0.000 0.000) 1.090 11.209 ( -2.981 0.000 0.000) 2.981 11.230 ( 10.291 0.000 0.000) 10.291 11.337 ( -1.678 0.000 0.000) 1.678 11.501 ( 1.067 0.000 0.000) 1.067 12.102 ( 14.183 0.000 0.000) 14.183 12.156 ( -5.710 0.000 0.000) 5.710 12.175 ( 3.152 0.000 0.000) 3.152 12.258 ( -4.524 0.000 0.000) 4.524 12.530 ( 11.035 0.000 0.000) 11.035 12.686 ( 16.246 0.000 0.000) 16.246 13.059 ( -4.091 0.000 0.000) 4.091 13.141 ( -4.920 0.000 0.000) 4.920 13.361 ( 3.898 0.000 0.000) 3.898 13.940 ( -9.630 0.000 0.000) 9.630 ======================= Grid point 3 (4/45) ======================= q-point: ( 0.33 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 125 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.975 ( 44.496 0.000 0.000) 44.496 2.990 ( 46.778 0.000 0.000) 46.778 5.171 ( 17.977 0.000 0.000) 17.977 5.697 ( -36.666 0.000 0.000) 36.666 5.812 ( 24.404 0.000 0.000) 24.404 7.409 ( 19.734 0.000 0.000) 19.734 7.883 ( 25.526 0.000 0.000) 25.526 7.926 ( 19.649 0.000 0.000) 19.649 8.044 ( 14.545 0.000 0.000) 14.545 8.678 ( -15.517 0.000 0.000) 15.517 8.741 ( -22.559 0.000 0.000) 22.559 8.950 ( -3.179 0.000 0.000) 3.179 9.688 ( -26.085 0.000 0.000) 26.085 9.729 ( -14.258 0.000 0.000) 14.258 9.865 ( -11.848 0.000 0.000) 11.848 9.900 ( -18.221 0.000 0.000) 18.221 10.558 ( -5.518 0.000 0.000) 5.518 10.642 ( -11.673 0.000 0.000) 11.673 10.731 ( 0.577 0.000 0.000) 0.577 10.753 ( 2.709 0.000 0.000) 2.709 11.003 ( 2.121 0.000 0.000) 2.121 11.134 ( -0.156 0.000 0.000) 0.156 11.168 ( -1.387 0.000 0.000) 1.387 11.241 ( -5.621 0.000 0.000) 5.621 11.464 ( 11.649 0.000 0.000) 11.649 11.685 ( 11.247 0.000 0.000) 11.247 12.032 ( -5.602 0.000 0.000) 5.602 12.129 ( -7.958 0.000 0.000) 7.958 12.252 ( 5.442 0.000 0.000) 5.442 12.351 ( 9.444 0.000 0.000) 9.444 12.614 ( -1.705 0.000 0.000) 1.705 12.895 ( -2.354 0.000 0.000) 2.354 13.093 ( 6.122 0.000 0.000) 6.122 13.101 ( 7.519 0.000 0.000) 7.519 13.419 ( 1.190 0.000 0.000) 1.190 13.718 ( -11.115 0.000 0.000) 11.115 ======================= Grid point 4 (5/45) ======================= q-point: ( 0.44 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 125 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.872 ( 41.191 0.000 0.000) 41.191 3.952 ( 45.143 0.000 0.000) 45.143 4.684 ( -35.959 0.000 0.000) 35.959 4.867 ( -42.157 0.000 0.000) 42.157 6.982 ( 58.321 0.000 0.000) 58.321 7.802 ( 9.455 0.000 0.000) 9.455 8.119 ( -29.592 0.000 0.000) 29.592 8.200 ( 7.680 0.000 0.000) 7.680 8.241 ( 6.216 0.000 0.000) 6.216 8.314 ( -19.095 0.000 0.000) 19.095 8.476 ( 24.174 0.000 0.000) 24.174 8.909 ( -9.686 0.000 0.000) 9.686 9.279 ( -2.448 0.000 0.000) 2.448 9.537 ( -3.288 0.000 0.000) 3.288 9.554 ( -17.721 0.000 0.000) 17.721 9.602 ( -9.304 0.000 0.000) 9.304 10.463 ( -3.017 0.000 0.000) 3.017 10.466 ( -4.107 0.000 0.000) 4.107 10.751 ( 1.165 0.000 0.000) 1.165 10.771 ( -0.490 0.000 0.000) 0.490 11.070 ( 3.653 0.000 0.000) 3.653 11.130 ( 0.046 0.000 0.000) 0.046 11.134 ( -2.229 0.000 0.000) 2.229 11.153 ( -2.507 0.000 0.000) 2.507 11.705 ( 11.255 0.000 0.000) 11.255 11.855 ( 4.950 0.000 0.000) 4.950 11.933 ( -10.446 0.000 0.000) 10.446 11.935 ( -3.729 0.000 0.000) 3.729 12.438 ( 12.037 0.000 0.000) 12.037 12.494 ( 4.520 0.000 0.000) 4.520 12.561 ( -2.541 0.000 0.000) 2.541 12.704 ( -12.263 0.000 0.000) 12.263 13.278 ( 8.078 0.000 0.000) 8.078 13.288 ( 9.446 0.000 0.000) 9.446 13.398 ( -3.344 0.000 0.000) 3.344 13.492 ( -10.284 0.000 0.000) 10.284 ======================= Grid point 10 (6/45) ======================= q-point: ( 0.11 0.11 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.436 ( 34.838 33.929 0.000) 48.629 1.907 ( 43.655 45.272 0.000) 62.891 2.415 ( 56.855 56.131 0.000) 79.895 6.632 ( -13.239 7.405 0.000) 15.169 6.671 ( -0.256 -0.589 0.000) 0.642 6.719 ( 15.308 -12.700 0.000) 19.890 6.955 ( 8.242 15.579 0.000) 17.625 7.397 ( -5.741 -4.099 0.000) 7.054 7.594 ( 5.745 4.473 0.000) 7.281 9.071 ( -5.903 -5.538 0.000) 8.094 9.147 ( -1.980 -3.161 0.000) 3.730 9.169 ( -4.581 -3.843 0.000) 5.979 10.152 ( -5.041 -4.884 0.000) 7.019 10.353 ( -5.076 -10.597 0.000) 11.750 10.403 ( -15.763 -8.326 0.000) 17.827 10.563 ( -3.700 -2.334 0.000) 4.375 10.598 ( -12.851 -12.522 0.000) 17.943 10.747 ( 4.708 -1.785 0.000) 5.035 10.814 ( 1.026 4.194 0.000) 4.318 10.850 ( 3.174 -4.835 0.000) 5.784 10.950 ( -1.999 6.498 0.000) 6.798 11.074 ( -0.437 -1.223 0.000) 1.298 11.136 ( -7.215 -7.675 0.000) 10.534 11.398 ( -3.313 0.088 0.000) 3.314 11.430 ( -3.445 -6.651 0.000) 7.490 11.548 ( -3.206 -1.980 0.000) 3.768 11.739 ( 8.827 8.265 0.000) 12.092 12.137 ( 2.374 1.279 0.000) 2.697 12.272 ( -3.602 -1.141 0.000) 3.779 12.291 ( 7.712 1.629 0.000) 7.882 12.385 ( 2.717 11.139 0.000) 11.465 12.465 ( 16.256 11.841 0.000) 20.111 13.141 ( -3.236 -5.980 0.000) 6.799 13.216 ( -5.456 2.968 0.000) 6.211 13.259 ( 3.410 -1.924 0.000) 3.915 14.097 ( -2.974 -3.078 0.000) 4.280 ======================= Grid point 11 (7/45) ======================= q-point: ( 0.22 0.11 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.264 ( 42.227 22.141 0.000) 47.679 2.647 ( 35.289 50.906 0.000) 61.941 3.962 ( 77.992 18.683 0.000) 80.198 6.231 ( -23.363 9.001 0.000) 25.037 6.465 ( -21.708 9.195 0.000) 23.575 7.172 ( 12.555 11.901 0.000) 17.299 7.174 ( 26.343 -13.045 0.000) 29.396 7.424 ( 11.212 -11.819 0.000) 16.291 7.807 ( 13.361 5.109 0.000) 14.304 8.926 ( -8.142 -4.014 0.000) 9.077 9.022 ( -8.596 -3.592 0.000) 9.317 9.046 ( -8.667 -1.764 0.000) 8.845 9.967 ( -22.940 -14.868 0.000) 27.337 9.987 ( -11.166 -7.562 0.000) 13.486 10.124 ( -15.888 -3.578 0.000) 16.285 10.246 ( -17.084 -6.923 0.000) 18.434 10.539 ( -0.559 -4.577 0.000) 4.611 10.755 ( -4.676 0.394 0.000) 4.693 10.840 ( -1.406 8.236 0.000) 8.356 10.901 ( -0.166 3.594 0.000) 3.598 11.007 ( -0.627 -11.075 0.000) 11.093 11.007 ( 3.354 -0.668 0.000) 3.420 11.115 ( 3.413 -0.609 0.000) 3.467 11.306 ( -0.787 3.946 0.000) 4.023 11.412 ( 0.224 0.563 0.000) 0.606 11.457 ( -1.292 -2.762 0.000) 3.049 12.013 ( 11.735 0.406 0.000) 11.743 12.129 ( -2.726 -3.334 0.000) 4.307 12.198 ( -1.656 -1.637 0.000) 2.329 12.365 ( 0.524 7.982 0.000) 8.000 12.556 ( 9.282 2.963 0.000) 9.743 12.800 ( 14.498 9.291 0.000) 17.220 13.049 ( -5.312 -7.287 0.000) 9.018 13.089 ( -4.536 2.671 0.000) 5.264 13.348 ( 4.387 -1.356 0.000) 4.591 13.968 ( -8.720 0.736 0.000) 8.751 ======================= Grid point 12 (8/45) ======================= q-point: ( 0.33 0.11 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.153 ( 42.285 16.473 0.000) 45.380 3.427 ( 38.959 37.841 0.000) 54.312 5.291 ( 17.347 11.433 0.000) 20.776 5.856 ( -34.268 14.544 0.000) 37.227 5.904 ( 23.038 8.964 0.000) 24.720 7.492 ( 17.534 6.388 0.000) 18.662 7.741 ( 25.386 -12.558 0.000) 28.322 7.752 ( 17.400 -14.875 0.000) 22.891 8.123 ( 16.469 7.771 0.000) 18.210 8.714 ( -12.762 2.123 0.000) 12.937 8.742 ( -21.393 0.029 0.000) 21.393 8.861 ( -5.268 -8.594 0.000) 10.080 9.500 ( -21.124 -12.615 0.000) 24.604 9.684 ( -16.614 -11.552 0.000) 20.235 9.769 ( -15.841 0.491 0.000) 15.848 9.934 ( -13.466 -2.037 0.000) 13.619 10.502 ( -2.867 -4.037 0.000) 4.951 10.580 ( -9.499 -3.348 0.000) 10.071 10.827 ( -0.032 5.320 0.000) 5.320 10.862 ( -1.142 8.331 0.000) 8.409 11.027 ( -0.069 -8.718 0.000) 8.719 11.075 ( 3.493 1.721 0.000) 3.894 11.120 ( -2.093 2.475 0.000) 3.242 11.316 ( -4.577 3.448 0.000) 5.730 11.422 ( 9.578 -3.521 0.000) 10.205 11.630 ( 10.982 -5.123 0.000) 12.118 12.040 ( -5.551 -5.442 0.000) 7.774 12.044 ( -3.785 1.443 0.000) 4.050 12.343 ( 10.260 -0.702 0.000) 10.284 12.367 ( 0.493 7.897 0.000) 7.913 12.635 ( -1.128 2.245 0.000) 2.512 12.858 ( -8.674 -4.118 0.000) 9.602 13.112 ( 10.965 3.303 0.000) 11.451 13.119 ( 7.133 1.435 0.000) 7.276 13.418 ( 1.836 -0.255 0.000) 1.854 13.757 ( -10.968 2.316 0.000) 11.210 ======================= Grid point 13 (9/45) ======================= q-point: ( 0.44 0.11 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.014 ( 39.815 13.317 0.000) 41.983 4.254 ( 39.513 27.046 0.000) 47.883 4.812 ( -34.731 12.084 0.000) 36.774 5.087 ( -38.574 20.048 0.000) 43.473 7.040 ( 56.742 5.677 0.000) 57.026 7.839 ( 7.933 3.218 0.000) 8.561 8.069 ( 5.073 -11.526 0.000) 12.593 8.083 ( 4.034 -10.128 0.000) 10.902 8.145 ( -21.650 -0.856 0.000) 21.666 8.397 ( -16.530 -0.035 0.000) 16.530 8.500 ( 18.420 5.443 0.000) 19.207 8.814 ( -5.031 -6.300 0.000) 8.062 9.156 ( -4.993 -8.060 0.000) 9.481 9.350 ( -14.940 -12.263 0.000) 19.328 9.556 ( -2.913 -0.348 0.000) 2.934 9.660 ( -11.486 -0.921 0.000) 11.522 10.441 ( -3.066 -1.711 0.000) 3.511 10.441 ( -2.261 -1.856 0.000) 2.925 10.831 ( 0.411 4.028 0.000) 4.049 10.844 ( -0.763 6.274 0.000) 6.320 11.026 ( 0.669 -5.955 0.000) 5.993 11.062 ( -2.562 -3.021 0.000) 3.961 11.160 ( 4.031 5.217 0.000) 6.593 11.238 ( -3.425 4.586 0.000) 5.724 11.634 ( 9.915 -5.007 0.000) 11.107 11.779 ( -0.204 -4.910 0.000) 4.914 11.929 ( -2.225 -3.013 0.000) 3.746 11.967 ( -3.282 0.979 0.000) 3.425 12.436 ( 6.852 0.772 0.000) 6.896 12.509 ( 5.197 0.889 0.000) 5.272 12.593 ( -2.256 3.235 0.000) 3.944 12.644 ( -10.611 -3.468 0.000) 11.164 13.308 ( 7.269 2.538 0.000) 7.699 13.308 ( 9.907 1.438 0.000) 10.011 13.412 ( -2.446 1.195 0.000) 2.722 13.527 ( -10.792 2.432 0.000) 11.063 ======================= Grid point 19 (10/45) ======================= q-point: ( 0.11 0.22 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.257 ( 22.769 42.621 0.000) 48.322 2.663 ( 50.844 35.286 0.000) 61.889 3.956 ( 18.304 77.985 0.000) 80.104 6.317 ( 16.708 -22.884 0.000) 28.334 6.436 ( 6.999 -23.324 0.000) 24.352 6.895 ( -12.962 15.785 0.000) 20.424 7.361 ( 2.823 20.461 0.000) 20.655 7.499 ( -8.652 17.249 0.000) 19.298 7.767 ( 5.522 11.249 0.000) 12.531 8.900 ( -5.342 -11.057 0.000) 12.280 9.032 ( -3.712 -9.724 0.000) 10.409 9.060 ( -0.438 -4.788 0.000) 4.808 9.953 ( -12.896 -20.484 0.000) 24.205 10.000 ( -7.159 -10.531 0.000) 12.734 10.192 ( -5.209 -22.100 0.000) 22.705 10.204 ( -7.613 -14.025 0.000) 15.958 10.590 ( -4.682 3.601 0.000) 5.906 10.670 ( 2.235 -4.893 0.000) 5.379 10.811 ( 6.014 0.224 0.000) 6.018 10.862 ( -2.317 -3.305 0.000) 4.036 11.021 ( 2.682 -0.899 0.000) 2.828 11.065 ( -4.965 -1.699 0.000) 5.248 11.168 ( -4.841 10.430 0.000) 11.498 11.298 ( 4.335 -5.174 0.000) 6.750 11.304 ( -2.121 -4.410 0.000) 4.893 11.523 ( 0.194 2.520 0.000) 2.527 12.036 ( -1.067 14.140 0.000) 14.181 12.095 ( -2.035 -4.319 0.000) 4.775 12.216 ( -3.412 -2.710 0.000) 4.357 12.318 ( 9.255 1.872 0.000) 9.442 12.689 ( 0.021 15.119 0.000) 15.119 12.693 ( 11.772 6.692 0.000) 13.541 13.029 ( -2.520 -2.089 0.000) 3.273 13.184 ( 3.860 -4.807 0.000) 6.165 13.297 ( -5.570 3.763 0.000) 6.722 13.964 ( 0.608 -9.036 0.000) 9.056 ======================= Grid point 20 (11/45) ======================= q-point: ( 0.22 0.22 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.870 ( 34.100 34.291 0.000) 48.360 3.674 ( 38.981 40.425 0.000) 56.158 4.658 ( 50.282 49.389 0.000) 70.481 6.492 ( -20.511 14.439 0.000) 25.084 6.569 ( 0.374 -8.660 0.000) 8.668 6.818 ( 25.808 -18.651 0.000) 31.842 7.163 ( -10.531 -0.091 0.000) 10.531 7.539 ( 8.661 18.288 0.000) 20.235 8.013 ( 15.684 14.422 0.000) 21.307 8.768 ( -7.046 -11.924 0.000) 13.850 8.902 ( -8.525 -7.892 0.000) 11.618 8.995 ( -6.773 -2.882 0.000) 7.361 9.597 ( -18.743 -17.452 0.000) 25.610 9.763 ( -14.687 -13.976 0.000) 20.274 10.017 ( -11.207 -8.853 0.000) 14.282 10.030 ( -9.335 -12.420 0.000) 15.537 10.492 ( -4.105 0.618 0.000) 4.151 10.720 ( -9.974 -3.651 0.000) 10.621 10.778 ( 5.820 -7.186 0.000) 9.247 10.937 ( 4.627 -0.827 0.000) 4.700 10.995 ( -1.995 -0.375 0.000) 2.030 11.053 ( -4.595 6.115 0.000) 7.649 11.116 ( 5.761 1.276 0.000) 5.901 11.315 ( 4.299 -7.081 0.000) 8.284 11.323 ( -0.721 -1.067 0.000) 1.288 11.481 ( -3.570 4.731 0.000) 5.927 12.045 ( -3.481 -4.522 0.000) 5.707 12.099 ( -0.508 -6.184 0.000) 6.205 12.180 ( 6.220 12.315 0.000) 13.796 12.508 ( 6.524 5.152 0.000) 8.313 12.694 ( 1.227 10.020 0.000) 10.095 12.878 ( 0.160 -6.636 0.000) 6.638 13.020 ( 6.620 9.394 0.000) 11.492 13.152 ( -5.621 2.666 0.000) 6.221 13.300 ( 5.554 -3.266 0.000) 6.443 13.918 ( -5.158 -5.692 0.000) 7.681 ======================= Grid point 21 (12/45) ======================= q-point: ( 0.33 0.22 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.627 ( 37.418 27.911 0.000) 46.681 4.353 ( 29.661 45.803 0.000) 54.568 5.647 ( 9.722 22.612 0.000) 24.613 6.195 ( 22.664 19.334 0.000) 29.790 6.247 ( -26.438 20.454 0.000) 33.426 7.266 ( 16.047 -23.386 0.000) 28.362 7.404 ( 26.951 -18.676 0.000) 32.789 7.737 ( 9.478 9.842 0.000) 13.664 8.367 ( 17.248 15.332 0.000) 23.078 8.564 ( -14.667 -18.988 0.000) 23.993 8.734 ( -4.525 -1.392 0.000) 4.735 8.739 ( -18.537 -0.372 0.000) 18.540 9.251 ( -13.200 -10.479 0.000) 16.854 9.432 ( -16.320 -11.752 0.000) 20.111 9.741 ( -14.164 -3.503 0.000) 14.591 9.804 ( -12.306 -9.789 0.000) 15.725 10.427 ( -2.270 -2.406 0.000) 3.307 10.516 ( -8.324 -3.013 0.000) 8.852 10.853 ( 1.496 -7.429 0.000) 7.578 10.942 ( -2.493 4.776 0.000) 5.388 10.964 ( -1.273 -0.395 0.000) 1.333 11.002 ( -1.606 -4.514 0.000) 4.791 11.242 ( 4.340 4.249 0.000) 6.074 11.363 ( 6.463 -0.635 0.000) 6.494 11.409 ( -2.751 6.401 0.000) 6.967 11.478 ( 10.198 -9.149 0.000) 13.700 11.938 ( -6.687 -4.125 0.000) 7.857 12.014 ( -5.181 -6.771 0.000) 8.526 12.391 ( 9.918 7.101 0.000) 12.198 12.539 ( -2.056 8.168 0.000) 8.423 12.720 ( -0.208 6.069 0.000) 6.072 12.776 ( -6.724 -2.199 0.000) 7.074 13.145 ( 5.716 0.630 0.000) 5.750 13.215 ( 8.978 5.458 0.000) 10.507 13.401 ( 3.570 -1.507 0.000) 3.875 13.763 ( -9.144 -2.332 0.000) 9.437 ======================= Grid point 22 (13/45) ======================= q-point: ( 0.44 0.22 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.406 ( 36.618 23.517 0.000) 43.520 4.952 ( 21.368 35.431 0.000) 41.376 5.186 ( -24.209 24.514 0.000) 34.453 5.636 ( -30.502 30.470 0.000) 43.113 7.227 ( 51.346 12.418 0.000) 52.826 7.657 ( 18.476 -24.724 0.000) 30.865 7.818 ( 9.803 -15.844 0.000) 18.632 7.905 ( 6.190 1.605 0.000) 6.395 8.180 ( -34.465 3.592 0.000) 34.651 8.182 ( -19.801 -18.383 0.000) 27.019 8.706 ( 13.360 12.634 0.000) 18.388 8.761 ( 6.294 -0.019 0.000) 6.294 9.044 ( -7.232 -2.666 0.000) 7.708 9.111 ( -13.368 -9.346 0.000) 16.311 9.499 ( -6.303 -4.230 0.000) 7.591 9.550 ( -9.959 -9.019 0.000) 13.436 10.393 ( -0.929 -2.613 0.000) 2.773 10.403 ( -2.463 -2.128 0.000) 3.255 10.874 ( 1.154 -5.292 0.000) 5.416 10.903 ( -1.354 2.896 0.000) 3.197 10.939 ( -0.859 -1.067 0.000) 1.370 10.950 ( -2.285 -6.883 0.000) 7.252 11.311 ( 2.932 7.291 0.000) 7.859 11.363 ( -2.043 7.140 0.000) 7.427 11.549 ( 9.678 -2.473 0.000) 9.989 11.675 ( 4.909 -6.302 0.000) 7.989 11.792 ( -4.135 -7.704 0.000) 8.744 11.892 ( -6.660 -8.987 0.000) 11.186 12.498 ( -0.312 5.030 0.000) 5.039 12.519 ( 0.666 1.913 0.000) 2.026 12.673 ( -0.833 5.804 0.000) 5.864 12.701 ( -1.427 6.962 0.000) 7.107 13.320 ( 9.163 -0.780 0.000) 9.196 13.361 ( 4.620 2.009 0.000) 5.038 13.439 ( 0.398 1.091 0.000) 1.162 13.557 ( -10.027 -0.052 0.000) 10.027 ======================= Grid point 28 (14/45) ======================= q-point: ( 0.11 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.172 ( 16.720 44.349 0.000) 47.396 3.427 ( 39.021 37.372 0.000) 54.031 5.330 ( 14.516 27.101 0.000) 30.744 5.795 ( 9.227 -35.894 0.000) 37.061 5.962 ( 14.507 14.362 0.000) 20.413 7.278 ( -12.098 20.388 0.000) 23.707 7.784 ( -9.332 20.303 0.000) 22.345 7.916 ( -6.898 15.443 0.000) 16.914 8.096 ( 11.866 22.013 0.000) 25.008 8.613 ( -6.106 -16.433 0.000) 17.531 8.720 ( -2.013 -21.228 0.000) 21.323 8.976 ( 1.338 -2.054 0.000) 2.451 9.589 ( -11.778 -13.135 0.000) 17.642 9.603 ( -7.584 -24.260 0.000) 25.418 9.857 ( -1.373 -12.422 0.000) 12.497 9.861 ( -3.912 -17.347 0.000) 17.783 10.560 ( 2.182 -5.114 0.000) 5.560 10.667 ( -4.320 3.465 0.000) 5.538 10.709 ( 2.990 -9.170 0.000) 9.645 10.821 ( 4.044 -0.131 0.000) 4.046 11.017 ( 2.043 1.361 0.000) 2.455 11.022 ( -8.148 -1.733 0.000) 8.331 11.192 ( -4.168 -6.048 0.000) 7.345 11.219 ( 3.224 -2.812 0.000) 4.278 11.402 ( -5.080 11.450 0.000) 12.526 11.679 ( -1.107 9.316 0.000) 9.382 11.996 ( -2.091 -4.603 0.000) 5.055 12.055 ( -5.988 -7.490 0.000) 9.589 12.360 ( 2.335 11.044 0.000) 11.288 12.389 ( 10.791 5.543 0.000) 12.132 12.692 ( 5.292 -3.607 0.000) 6.404 12.941 ( 2.291 4.190 0.000) 4.776 13.083 ( 0.650 -0.433 0.000) 0.781 13.124 ( 2.363 8.624 0.000) 8.942 13.354 ( -5.792 1.244 0.000) 5.924 13.746 ( 1.722 -11.207 0.000) 11.339 ======================= Grid point 29 (15/45) ======================= q-point: ( 0.22 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.646 ( 27.566 39.131 0.000) 47.866 4.372 ( 46.057 29.909 0.000) 54.916 5.753 ( 25.879 36.684 0.000) 44.894 6.063 ( 15.647 -33.114 0.000) 36.625 6.405 ( 26.509 -1.119 0.000) 26.533 6.936 ( -19.048 21.369 0.000) 28.627 7.429 ( -20.833 18.758 0.000) 28.033 7.884 ( 0.858 12.866 0.000) 12.894 8.390 ( 15.553 19.765 0.000) 25.151 8.420 ( -12.586 -20.743 0.000) 24.262 8.667 ( -2.517 -16.074 0.000) 16.269 8.951 ( -4.797 0.262 0.000) 4.804 9.308 ( -13.560 -9.058 0.000) 16.307 9.414 ( -9.428 -18.389 0.000) 20.665 9.748 ( -6.563 -15.851 0.000) 17.156 9.778 ( -6.676 -11.617 0.000) 13.398 10.555 ( -6.136 4.671 0.000) 7.711 10.605 ( -3.159 -5.480 0.000) 6.325 10.695 ( -0.407 -2.805 0.000) 2.834 10.885 ( 1.121 -3.877 0.000) 4.036 10.898 ( -3.072 -6.434 0.000) 7.130 11.097 ( 5.085 1.280 0.000) 5.244 11.145 ( 1.498 -8.881 0.000) 9.007 11.212 ( -4.241 2.552 0.000) 4.950 11.346 ( 1.353 5.336 0.000) 5.505 11.632 ( -2.994 8.084 0.000) 8.621 11.926 ( -5.997 -9.608 0.000) 11.326 11.952 ( -2.578 -4.182 0.000) 4.913 12.457 ( 6.423 12.297 0.000) 13.873 12.626 ( 10.040 6.886 0.000) 12.174 12.773 ( 1.913 -2.966 0.000) 3.530 12.913 ( -3.271 9.937 0.000) 10.461 13.152 ( 0.044 -4.301 0.000) 4.301 13.204 ( 5.098 7.946 0.000) 9.441 13.237 ( 0.135 -1.420 0.000) 1.427 13.749 ( -1.941 -9.727 0.000) 9.919 ======================= Grid point 30 (16/45) ======================= q-point: ( 0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.280 ( 32.303 33.874 0.000) 46.808 5.194 ( 31.554 33.601 0.000) 46.094 6.173 ( -7.402 23.193 0.000) 24.345 6.380 ( 10.517 -19.277 0.000) 21.960 6.727 ( 23.466 28.421 0.000) 36.856 7.029 ( -13.069 12.134 0.000) 17.833 7.036 ( 28.237 -10.189 0.000) 30.019 7.895 ( -0.488 3.069 0.000) 3.108 8.089 ( -18.523 -25.359 0.000) 31.403 8.646 ( 1.374 -8.909 0.000) 9.014 8.689 ( 4.936 7.255 0.000) 8.775 8.778 ( -4.601 12.022 0.000) 12.872 9.072 ( -8.411 -5.343 0.000) 9.965 9.219 ( -9.684 -8.617 0.000) 12.963 9.570 ( -12.841 -11.289 0.000) 17.098 9.594 ( -8.086 -10.979 0.000) 13.635 10.431 ( -5.143 3.174 0.000) 6.043 10.453 ( -7.144 -2.650 0.000) 7.620 10.716 ( 0.122 -6.411 0.000) 6.412 10.856 ( -3.299 -7.055 0.000) 7.788 10.880 ( 1.103 -4.949 0.000) 5.071 11.075 ( -5.138 5.638 0.000) 7.628 11.248 ( 7.210 -11.999 0.000) 13.999 11.265 ( 7.323 -0.399 0.000) 7.334 11.411 ( 4.584 3.458 0.000) 5.742 11.573 ( -2.287 7.902 0.000) 8.226 11.800 ( -6.336 -11.657 0.000) 13.267 11.869 ( -5.639 -2.708 0.000) 6.256 12.585 ( 4.475 9.627 0.000) 10.616 12.704 ( -2.798 7.696 0.000) 8.189 12.785 ( 0.318 2.533 0.000) 2.553 12.877 ( -0.721 8.164 0.000) 8.196 13.132 ( 4.171 -1.812 0.000) 4.548 13.316 ( 4.898 4.059 0.000) 6.361 13.350 ( 5.821 -3.395 0.000) 6.739 13.664 ( -5.931 -6.628 0.000) 8.894 ======================= Grid point 31 (17/45) ======================= q-point: ( 0.44 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.966 ( 32.731 29.428 0.000) 44.015 5.584 ( -16.822 27.022 0.000) 31.830 5.846 ( 10.438 33.548 0.000) 35.135 6.253 ( -18.676 24.299 0.000) 30.647 7.125 ( 17.690 -22.606 0.000) 28.705 7.418 ( 8.970 -20.387 0.000) 22.274 7.543 ( 26.182 10.809 0.000) 28.326 7.775 ( 2.498 -13.406 0.000) 13.637 7.834 ( -6.181 -9.262 0.000) 11.135 8.320 ( -27.185 9.970 0.000) 28.956 8.762 ( 10.127 -1.678 0.000) 10.265 8.944 ( 3.961 5.989 0.000) 7.180 8.967 ( -1.686 -4.209 0.000) 4.534 9.010 ( -10.434 1.549 0.000) 10.548 9.331 ( -5.016 -10.035 0.000) 11.219 9.415 ( -8.513 -4.436 0.000) 9.599 10.355 ( -2.130 -1.512 0.000) 2.612 10.355 ( -1.996 -0.230 0.000) 2.009 10.727 ( 1.482 -8.387 0.000) 8.517 10.780 ( -3.418 -8.575 0.000) 9.231 10.914 ( 2.030 -0.452 0.000) 2.080 10.976 ( -4.018 3.747 0.000) 5.494 11.394 ( 4.184 -1.041 0.000) 4.312 11.397 ( 5.168 -11.420 0.000) 12.535 11.543 ( 7.982 1.623 0.000) 8.146 11.551 ( 1.117 3.504 0.000) 3.677 11.657 ( -7.146 -10.177 0.000) 12.436 11.720 ( -8.276 -0.696 0.000) 8.305 12.615 ( -0.469 6.668 0.000) 6.684 12.629 ( -2.510 7.097 0.000) 7.528 12.820 ( 2.421 6.650 0.000) 7.077 12.857 ( -1.120 6.648 0.000) 6.742 13.281 ( 7.806 -2.105 0.000) 8.085 13.378 ( 0.024 -0.266 0.000) 0.267 13.442 ( 3.079 -1.144 0.000) 3.285 13.526 ( -6.224 -2.587 0.000) 6.740 ======================= Grid point 37 (18/45) ======================= q-point: ( 0.11 0.44 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.094 ( 13.240 43.677 0.000) 45.639 4.205 ( 29.828 36.701 0.000) 47.294 4.941 ( 23.364 -32.910 0.000) 40.360 4.984 ( 10.931 -41.154 0.000) 42.581 7.034 ( 5.112 56.041 0.000) 56.274 7.704 ( -9.456 16.037 0.000) 18.617 8.114 ( -0.398 -33.122 0.000) 33.125 8.119 ( -11.396 8.901 0.000) 14.460 8.149 ( -4.788 6.723 0.000) 8.254 8.226 ( -7.912 -20.508 0.000) 21.981 8.547 ( 6.867 19.960 0.000) 21.108 8.947 ( 3.032 -11.817 0.000) 12.199 9.235 ( -3.604 2.348 0.000) 4.302 9.439 ( -8.211 -0.920 0.000) 8.263 9.548 ( -0.710 -17.559 0.000) 17.573 9.549 ( -4.945 -11.024 0.000) 12.082 10.490 ( 3.515 -0.672 0.000) 3.579 10.545 ( 4.842 -5.687 0.000) 7.469 10.734 ( -1.724 3.005 0.000) 3.465 10.792 ( 1.073 -2.376 0.000) 2.607 11.013 ( -9.083 1.254 0.000) 9.169 11.077 ( -6.696 -4.615 0.000) 8.132 11.080 ( 1.490 4.004 0.000) 4.273 11.159 ( 2.007 -3.298 0.000) 3.861 11.638 ( -5.633 11.053 0.000) 12.406 11.831 ( -2.130 5.020 0.000) 5.453 11.863 ( -6.044 -10.329 0.000) 11.967 11.912 ( -1.770 -3.510 0.000) 3.931 12.538 ( 4.382 5.952 0.000) 7.392 12.564 ( 10.404 10.769 0.000) 14.973 12.624 ( 5.008 -2.987 0.000) 5.832 12.801 ( 8.293 -10.964 0.000) 13.747 13.225 ( -4.802 7.530 0.000) 8.931 13.313 ( 2.119 9.210 0.000) 9.450 13.337 ( -5.545 -3.073 0.000) 6.340 13.516 ( 1.788 -10.438 0.000) 10.590 ======================= Grid point 38 (19/45) ======================= q-point: ( 0.22 0.44 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.479 ( 22.818 40.124 0.000) 46.158 4.985 ( 41.341 29.051 0.000) 50.528 5.305 ( 19.180 -38.491 0.000) 43.005 5.575 ( 34.787 -26.422 0.000) 43.684 7.195 ( 8.900 46.613 0.000) 47.455 7.425 ( -14.853 24.690 0.000) 28.813 7.806 ( -17.137 15.635 0.000) 23.198 7.944 ( -15.944 -23.248 0.000) 28.190 8.011 ( -8.534 -2.176 0.000) 8.807 8.166 ( 2.966 -32.066 0.000) 32.203 8.720 ( 8.161 12.282 0.000) 14.746 9.006 ( 1.916 -12.523 0.000) 12.668 9.161 ( -2.892 12.450 0.000) 12.782 9.246 ( -9.204 3.084 0.000) 9.707 9.423 ( -6.197 -13.306 0.000) 14.678 9.508 ( -3.764 -15.334 0.000) 15.790 10.539 ( -1.473 -0.022 0.000) 1.473 10.605 ( -0.158 -5.356 0.000) 5.358 10.666 ( -4.818 5.684 0.000) 7.451 10.784 ( -2.253 -5.485 0.000) 5.930 10.845 ( -4.344 2.369 0.000) 4.948 10.963 ( -1.805 -7.950 0.000) 8.153 11.145 ( 4.458 2.959 0.000) 5.350 11.202 ( 1.589 -2.148 0.000) 2.672 11.514 ( -4.967 9.275 0.000) 10.521 11.717 ( -7.150 -10.034 0.000) 12.320 11.776 ( -2.871 5.477 0.000) 6.184 11.870 ( -2.299 -3.889 0.000) 4.518 12.657 ( 6.170 6.490 0.000) 8.955 12.735 ( 5.063 -1.682 0.000) 5.335 12.799 ( 10.185 9.361 0.000) 13.834 12.993 ( 8.048 -8.189 0.000) 11.482 13.101 ( -6.032 6.150 0.000) 8.614 13.201 ( -5.281 -3.044 0.000) 6.096 13.361 ( 2.219 7.635 0.000) 7.951 13.541 ( 0.093 -9.618 0.000) 9.619 ======================= Grid point 39 (20/45) ======================= q-point: ( 0.33 0.44 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.015 ( 27.720 35.866 0.000) 45.329 5.756 ( 23.260 -34.468 0.000) 41.582 5.796 ( 33.284 25.317 0.000) 41.818 6.256 ( 24.512 -16.399 0.000) 29.492 7.008 ( -13.788 20.950 0.000) 25.081 7.440 ( -17.317 20.770 0.000) 27.042 7.514 ( -6.007 -12.438 0.000) 13.813 7.611 ( 6.373 23.988 0.000) 24.820 7.800 ( -11.433 -12.167 0.000) 16.696 8.284 ( 8.676 -26.007 0.000) 27.416 8.801 ( -2.676 8.134 0.000) 8.563 8.997 ( -3.381 -7.162 0.000) 7.919 9.087 ( -5.481 6.643 0.000) 8.612 9.145 ( 1.992 11.370 0.000) 11.543 9.315 ( -3.948 -13.925 0.000) 14.474 9.367 ( -9.983 -8.804 0.000) 13.311 10.437 ( -5.459 2.102 0.000) 5.850 10.545 ( -6.269 7.028 0.000) 9.417 10.554 ( -2.888 -8.226 0.000) 8.718 10.702 ( -5.110 -7.617 0.000) 9.172 10.859 ( 3.733 3.481 0.000) 5.104 11.008 ( 4.750 -10.111 0.000) 11.171 11.181 ( -3.394 4.303 0.000) 5.481 11.246 ( 3.530 -1.760 0.000) 3.944 11.480 ( 2.112 3.536 0.000) 4.118 11.589 ( -4.336 -7.646 0.000) 8.790 11.716 ( -2.878 5.633 0.000) 6.325 11.809 ( -3.806 -3.384 0.000) 5.093 12.756 ( 2.503 6.300 0.000) 6.779 12.827 ( 3.724 -0.086 0.000) 3.725 12.860 ( -3.885 7.077 0.000) 8.073 12.974 ( -5.173 -2.007 0.000) 5.549 13.131 ( 7.862 4.611 0.000) 9.114 13.256 ( 7.202 -5.647 0.000) 9.152 13.398 ( 1.019 4.384 0.000) 4.501 13.515 ( -2.338 -6.844 0.000) 7.232 ======================= Grid point 40 (21/45) ======================= q-point: ( 0.44 0.44 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.608 ( 28.337 31.323 0.000) 42.238 6.132 ( -15.582 24.363 0.000) 28.920 6.228 ( 19.703 -25.172 0.000) 31.966 6.380 ( 5.515 10.657 0.000) 12.000 7.032 ( -0.184 -8.871 0.000) 8.873 7.111 ( -12.766 19.518 0.000) 23.322 7.248 ( 5.082 -20.428 0.000) 21.051 7.536 ( -14.043 -18.425 0.000) 23.166 8.090 ( 26.040 23.701 0.000) 35.211 8.542 ( 15.628 -18.025 0.000) 23.856 8.575 ( -18.350 13.731 0.000) 22.919 8.856 ( -11.319 -11.602 0.000) 16.209 9.029 ( 0.189 9.128 0.000) 9.130 9.118 ( -8.728 5.999 0.000) 10.591 9.218 ( 1.785 -2.776 0.000) 3.300 9.257 ( -1.887 -10.736 0.000) 10.901 10.375 ( -0.116 4.578 0.000) 4.579 10.426 ( -4.586 6.589 0.000) 8.028 10.534 ( 1.195 -9.462 0.000) 9.537 10.596 ( -4.437 -8.900 0.000) 9.945 10.956 ( 5.307 4.929 0.000) 7.243 11.075 ( -5.886 6.229 0.000) 8.570 11.129 ( 6.393 -11.469 0.000) 13.130 11.247 ( -3.843 -9.847 0.000) 10.570 11.572 ( 2.742 0.048 0.000) 2.743 11.574 ( 5.850 0.615 0.000) 5.883 11.648 ( -3.766 4.030 0.000) 5.515 11.705 ( -6.030 -1.413 0.000) 6.193 12.771 ( -0.219 7.469 0.000) 7.473 12.784 ( -2.150 7.373 0.000) 7.680 12.882 ( 1.456 -1.747 0.000) 2.274 12.905 ( -1.549 -2.867 0.000) 3.259 13.290 ( 6.447 3.618 0.000) 7.393 13.385 ( 4.693 -4.205 0.000) 6.301 13.387 ( -2.538 2.134 0.000) 3.315 13.458 ( -2.836 -3.787 0.000) 4.731 ======================= Grid point 91 (22/45) ======================= q-point: ( 0.11 0.11 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 125 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.151 ( 33.358 33.358 33.358) 57.778 2.151 ( 33.358 33.358 33.358) 57.778 2.813 ( 43.769 43.769 43.769) 75.811 6.655 ( -0.507 -0.507 -0.507) 0.878 6.817 ( 3.855 3.855 3.855) 6.676 6.817 ( 3.855 3.855 3.855) 6.676 6.840 ( 4.596 4.596 4.596) 7.961 7.470 ( -0.988 -0.988 -0.988) 1.712 7.470 ( -0.988 -0.988 -0.988) 1.712 9.049 ( -4.313 -4.313 -4.313) 7.470 9.049 ( -4.313 -4.313 -4.313) 7.470 9.179 ( -2.206 -2.206 -2.206) 3.820 10.094 ( -6.085 -6.085 -6.085) 10.539 10.261 ( -16.892 -16.892 -16.892) 29.259 10.345 ( -6.381 -6.381 -6.381) 11.052 10.345 ( -6.381 -6.381 -6.381) 11.052 10.686 ( -1.495 -1.495 -1.495) 2.589 10.686 ( -1.495 -1.495 -1.495) 2.589 10.855 ( -0.558 -0.558 -0.558) 0.966 10.855 ( -0.558 -0.558 -0.558) 0.966 10.868 ( 2.609 2.609 2.609) 4.519 11.072 ( 0.265 0.265 0.265) 0.460 11.212 ( -3.024 -3.024 -3.024) 5.238 11.212 ( -3.024 -3.024 -3.024) 5.238 11.517 ( -1.371 -1.371 -1.371) 2.375 11.517 ( -1.371 -1.371 -1.371) 2.375 11.771 ( 6.945 6.945 6.945) 12.029 12.203 ( -1.030 -1.030 -1.030) 1.784 12.203 ( -1.030 -1.030 -1.030) 1.784 12.392 ( 5.411 5.411 5.411) 9.372 12.392 ( 5.411 5.411 5.411) 9.372 12.590 ( 9.703 9.703 9.703) 16.805 13.169 ( -1.454 -1.454 -1.454) 2.518 13.198 ( -1.437 -1.437 -1.437) 2.489 13.198 ( -1.437 -1.437 -1.437) 2.489 14.092 ( -2.197 -2.197 -2.197) 3.806 ======================= Grid point 92 (23/45) ======================= q-point: ( 0.22 0.11 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.747 ( 33.447 30.344 29.131) 53.741 2.917 ( 33.904 31.856 33.726) 57.460 4.145 ( 72.443 19.004 18.581) 77.165 6.385 ( -22.052 6.667 13.081) 26.492 6.558 ( -20.491 10.346 9.362) 24.791 7.039 ( 12.761 12.983 -12.187) 21.907 7.237 ( 22.594 -11.710 5.591) 26.055 7.484 ( 13.861 -6.666 1.291) 15.435 7.667 ( 11.500 -2.132 -7.549) 13.921 8.891 ( -9.077 -1.503 -3.823) 9.964 8.959 ( -6.516 -5.064 -3.672) 9.033 9.074 ( -7.756 -0.710 0.070) 7.789 9.785 ( -20.617 -14.609 -14.334) 29.051 9.974 ( -11.132 -4.126 -4.543) 12.711 10.071 ( -17.826 -7.462 -4.962) 19.952 10.158 ( -13.344 -6.162 -7.542) 16.520 10.561 ( -2.505 -2.452 1.592) 3.850 10.724 ( -0.572 2.505 -6.954) 7.413 10.772 ( -5.180 -2.622 1.028) 5.896 10.885 ( 1.257 0.547 -2.373) 2.741 10.988 ( 4.018 -5.305 -1.926) 6.928 11.062 ( -1.560 1.849 3.180) 3.996 11.162 ( -1.015 -0.809 3.411) 3.650 11.209 ( 4.442 2.377 -3.872) 6.355 11.428 ( -5.122 2.622 0.384) 5.767 11.510 ( 2.389 -1.273 1.597) 3.143 12.017 ( 11.710 1.835 1.940) 12.010 12.113 ( -5.004 -0.239 -2.032) 5.406 12.184 ( -0.278 -2.824 -1.123) 3.051 12.413 ( 0.760 7.955 5.744) 9.841 12.646 ( 10.254 -0.646 6.608) 12.216 12.797 ( 9.191 6.138 -0.038) 11.052 13.078 ( -2.455 3.800 -1.166) 4.672 13.109 ( -4.723 -5.551 5.385) 9.062 13.284 ( 4.498 -1.314 -5.486) 7.215 13.981 ( -8.085 0.229 0.003) 8.088 ======================= Grid point 93 (24/45) ======================= q-point: ( 0.33 0.11 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.498 ( 37.564 21.510 22.990) 49.013 3.651 ( 36.292 27.336 27.341) 53.027 5.446 ( 25.005 11.187 14.341) 30.920 5.930 ( -25.398 11.841 7.841) 29.099 6.070 ( 6.111 10.851 15.012) 19.505 7.366 ( 17.933 7.405 -11.689) 22.650 7.694 ( 20.300 -10.767 -1.148) 23.008 7.897 ( 15.902 -14.215 7.251) 22.528 7.951 ( 20.030 3.656 -11.889) 23.578 8.662 ( -12.864 2.473 -4.223) 13.763 8.718 ( -18.751 -0.258 -2.476) 18.915 8.913 ( -4.913 -5.861 3.230) 8.302 9.381 ( -17.688 -12.607 -9.979) 23.903 9.652 ( -17.777 -1.205 -9.063) 19.990 9.696 ( -16.260 -9.406 -0.841) 18.804 9.890 ( -12.583 -2.305 -4.445) 13.543 10.510 ( -3.121 -3.261 1.450) 4.742 10.614 ( -7.062 -5.979 1.319) 9.347 10.748 ( 1.746 4.132 -5.723) 7.271 10.896 ( -0.541 2.306 -1.740) 2.939 10.991 ( -2.296 -2.926 -3.669) 5.225 11.027 ( -1.833 -1.604 0.451) 2.478 11.139 ( -0.255 5.481 1.780) 5.768 11.326 ( -3.718 4.298 1.227) 5.813 11.386 ( 9.798 -1.312 -2.708) 10.250 11.645 ( 8.446 -3.251 0.081) 9.051 11.995 ( -5.571 -0.533 -2.668) 6.200 12.043 ( -4.913 -0.592 -2.046) 5.355 12.360 ( 11.636 0.149 2.455) 11.893 12.436 ( 2.110 4.500 6.846) 8.460 12.698 ( -2.505 0.456 5.390) 5.961 12.915 ( -5.076 -3.259 3.312) 6.881 13.065 ( 7.984 0.226 -1.706) 8.167 13.162 ( 7.564 3.321 2.419) 8.608 13.356 ( 1.965 0.032 -5.330) 5.680 13.777 ( -10.861 1.798 1.122) 11.065 ======================= Grid point 94 (25/45) ======================= q-point: ( 0.44 0.11 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.287 ( 37.219 13.689 20.887) 44.821 4.422 ( 37.127 23.574 19.074) 47.937 5.039 ( -32.881 10.456 19.782) 39.772 5.212 ( -36.888 19.693 12.603) 43.674 7.095 ( 54.432 5.900 5.335) 55.011 7.740 ( 14.005 3.136 -9.488) 17.204 8.017 ( 6.087 -9.537 -4.041) 12.014 8.086 ( -6.288 -6.562 -1.605) 9.229 8.142 ( -19.038 1.509 -2.053) 19.207 8.246 ( -6.728 -0.500 -10.369) 12.370 8.535 ( 10.475 0.910 2.849) 10.893 8.850 ( -9.658 -6.836 2.629) 12.121 9.101 ( 0.238 -8.632 -4.377) 9.681 9.337 ( -14.757 -11.219 -1.558) 18.603 9.466 ( -1.159 -0.881 -7.371) 7.513 9.615 ( -12.792 -1.232 -4.101) 13.490 10.460 ( -0.542 -2.265 2.423) 3.361 10.499 ( -3.874 -3.695 3.247) 6.261 10.800 ( 2.944 3.066 -4.322) 6.062 10.861 ( -2.589 3.405 -1.565) 4.555 10.953 ( -0.506 -3.053 -5.125) 5.987 10.977 ( -2.462 -4.165 -3.412) 5.921 11.175 ( 3.233 6.884 1.344) 7.724 11.250 ( -3.608 5.772 1.138) 6.902 11.599 ( 9.936 -2.938 -3.089) 10.812 11.755 ( -1.557 -4.323 -2.256) 5.119 11.893 ( -1.545 0.070 -3.335) 3.676 11.942 ( -3.990 1.284 -2.575) 4.919 12.527 ( 6.854 -2.159 7.727) 10.552 12.560 ( 7.255 1.666 4.718) 8.813 12.640 ( -3.107 0.809 4.911) 5.868 12.733 ( -9.822 -4.143 7.144) 12.833 13.254 ( 7.438 2.452 -4.478) 9.022 13.336 ( 8.816 1.946 1.227) 9.111 13.355 ( -1.985 1.298 -4.445) 5.038 13.547 ( -10.798 2.096 1.466) 11.097 ======================= Grid point 101 (26/45) ======================= q-point: ( 0.22 0.22 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.297 ( 28.758 28.709 36.084) 54.344 3.794 ( 37.629 38.504 11.277) 55.006 4.830 ( 47.794 47.348 17.049) 69.403 6.539 ( -8.511 0.963 5.051) 9.944 6.743 ( -10.511 4.501 14.844) 18.737 6.914 ( 23.032 -17.690 8.976) 30.397 7.238 ( -9.568 -0.826 7.027) 11.900 7.413 ( 9.162 18.102 -11.365) 23.255 7.820 ( 14.437 13.061 -16.661) 25.624 8.788 ( -5.286 -9.051 0.645) 10.501 8.849 ( -7.462 -5.768 -4.047) 10.263 9.012 ( -6.898 -4.783 0.288) 8.399 9.438 ( -16.865 -16.417 -12.744) 26.764 9.771 ( -13.905 -14.584 -0.563) 20.158 9.927 ( -11.230 -7.804 -7.576) 15.634 9.955 ( -9.402 -12.188 -7.928) 17.315 10.515 ( -4.642 -1.009 1.573) 5.004 10.669 ( -4.359 -5.532 -5.870) 9.169 10.717 ( -2.427 -3.346 -3.710) 5.555 10.841 ( -0.302 -3.931 -7.842) 8.777 10.987 ( 0.412 2.782 0.489) 2.855 11.111 ( 0.277 1.402 4.399) 4.626 11.192 ( 3.689 1.993 5.032) 6.550 11.230 ( 3.489 1.014 -0.589) 3.680 11.439 ( 0.986 -2.774 3.336) 4.449 11.493 ( -1.982 2.992 0.855) 3.689 12.050 ( -3.650 -4.181 0.659) 5.589 12.101 ( -0.918 -5.146 -0.217) 5.232 12.185 ( 7.328 11.707 0.854) 13.838 12.556 ( 4.757 3.340 5.351) 7.900 12.710 ( -0.109 7.718 1.715) 7.907 12.884 ( 8.199 -0.445 -2.323) 8.533 13.022 ( -5.125 1.342 2.701) 5.946 13.182 ( -1.498 4.818 1.371) 5.229 13.241 ( 6.330 -2.678 -4.136) 8.022 13.930 ( -5.117 -5.333 0.279) 7.396 ======================= Grid point 102 (27/45) ======================= q-point: ( 0.33 0.22 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.950 ( 33.043 24.723 28.198) 49.982 4.469 ( 29.701 42.678 11.107) 53.168 5.797 ( 13.438 22.177 13.986) 29.462 6.243 ( -0.035 17.917 6.023) 18.903 6.415 ( -8.450 20.469 14.313) 26.367 7.292 ( 16.472 -19.609 2.183) 25.702 7.491 ( 23.328 -19.247 7.305) 31.113 7.630 ( 10.522 8.192 -9.109) 16.149 8.132 ( 13.234 9.266 -17.358) 23.712 8.612 ( -12.781 -8.452 -0.324) 15.326 8.717 ( -6.155 -0.198 -2.147) 6.522 8.767 ( -12.631 -5.023 2.964) 13.912 9.146 ( -10.560 -9.603 -8.134) 16.428 9.449 ( -16.145 -11.450 0.330) 19.796 9.621 ( -15.393 -3.470 -9.947) 18.653 9.756 ( -11.296 -9.921 -5.004) 15.845 10.442 ( -2.367 -3.092 0.964) 4.012 10.501 ( -7.231 -4.504 -1.389) 8.632 10.726 ( 0.407 -5.435 -9.697) 11.124 10.826 ( -1.308 -6.795 -10.703) 12.745 10.977 ( -1.147 2.396 2.740) 3.816 11.075 ( -3.123 2.368 5.448) 6.711 11.280 ( 3.206 5.694 2.631) 7.044 11.356 ( 8.058 -0.265 0.609) 8.085 11.428 ( -1.537 5.723 1.372) 6.083 11.546 ( 6.975 -6.094 3.414) 9.872 11.942 ( -6.382 -3.665 0.419) 7.371 12.020 ( -4.913 -4.214 0.087) 6.473 12.406 ( 10.434 5.560 1.678) 11.942 12.559 ( -2.458 6.606 2.278) 7.407 12.699 ( -1.476 -0.183 1.498) 2.111 12.872 ( -6.929 1.348 4.642) 8.448 13.094 ( 10.611 1.477 -3.236) 11.191 13.243 ( 5.515 3.818 1.304) 6.833 13.349 ( 3.675 -0.697 -3.854) 5.371 13.777 ( -9.048 -2.301 0.747) 9.366 ======================= Grid point 103 (28/45) ======================= q-point: ( 0.44 0.22 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.651 ( 33.445 21.372 21.662) 45.217 5.074 ( 24.983 34.522 11.486) 44.135 5.371 ( -26.054 21.338 16.727) 37.602 5.742 ( -29.847 29.222 10.270) 43.014 7.288 ( 48.866 12.668 5.886) 50.824 7.682 ( 19.072 -21.797 1.378) 28.995 7.813 ( 5.384 -0.001 -6.509) 8.448 7.839 ( 8.392 -12.573 -0.139) 15.117 8.113 ( -16.416 -14.513 -6.116) 22.749 8.200 ( -32.013 3.967 1.402) 32.289 8.609 ( 12.090 6.332 -6.770) 15.234 8.793 ( 3.025 1.855 1.255) 3.764 8.974 ( -4.931 -3.775 -5.098) 8.034 9.127 ( -13.737 -7.046 0.518) 15.448 9.387 ( -3.772 -5.979 -9.305) 11.686 9.498 ( -11.714 -9.360 -4.706) 15.715 10.406 ( -1.155 -2.961 0.656) 3.245 10.410 ( -1.603 -4.572 0.332) 4.857 10.741 ( 1.433 -6.847 -10.747) 12.823 10.786 ( -2.423 -7.495 -11.274) 13.754 10.963 ( 0.534 2.807 3.757) 4.720 11.005 ( -3.100 2.698 4.796) 6.316 11.338 ( 2.618 7.572 1.814) 8.214 11.386 ( -1.889 6.970 1.740) 7.428 11.548 ( 8.884 -1.357 -0.193) 8.989 11.673 ( 1.024 -3.799 -0.641) 3.987 11.825 ( -1.539 -6.221 2.088) 6.740 11.908 ( -5.698 -5.714 0.718) 8.101 12.514 ( 1.092 1.518 2.267) 2.939 12.557 ( -0.416 -1.055 3.699) 3.869 12.677 ( 2.248 5.194 1.450) 5.842 12.759 ( -4.170 5.071 3.674) 7.523 13.301 ( 7.662 1.160 -3.185) 8.378 13.362 ( 4.112 0.352 1.108) 4.273 13.392 ( 2.538 2.024 -3.743) 4.955 13.571 ( -10.228 -0.188 1.044) 10.283 ======================= Grid point 110 (29/45) ======================= q-point: ( 0.22 0.33 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.961 ( 24.725 34.210 27.525) 50.391 4.482 ( 43.011 29.888 10.583) 53.435 5.869 ( 25.905 32.348 11.159) 42.918 6.171 ( 14.434 -25.888 9.752) 31.203 6.510 ( 25.911 -4.314 10.097) 28.141 7.034 ( -18.623 18.753 8.764) 27.844 7.497 ( -19.781 18.457 6.126) 27.739 7.769 ( 0.908 13.933 -10.167) 17.272 8.116 ( 7.317 9.376 -17.209) 20.918 8.561 ( -3.924 -8.016 4.402) 9.952 8.655 ( -2.642 -14.331 -1.152) 14.618 8.926 ( -3.463 -1.891 -2.598) 4.724 9.148 ( -10.902 -11.306 -11.102) 19.234 9.423 ( -10.971 -16.324 -1.552) 19.729 9.651 ( -4.112 -16.019 -7.967) 18.358 9.749 ( -8.968 -9.915 -3.870) 13.918 10.524 ( -2.902 0.152 -3.009) 4.183 10.556 ( -6.202 -3.491 -4.631) 8.491 10.647 ( -3.228 -3.086 -4.570) 6.390 10.756 ( -6.076 -3.677 -7.183) 10.101 10.987 ( 0.409 -2.671 2.319) 3.561 11.093 ( 3.940 -3.118 3.697) 6.238 11.190 ( 4.250 0.235 5.670) 7.090 11.325 ( -1.506 3.451 5.123) 6.358 11.378 ( 2.631 0.682 3.164) 4.171 11.618 ( -2.304 6.992 -1.429) 7.499 11.930 ( -3.658 -7.330 0.659) 8.219 11.978 ( -3.785 -5.764 1.830) 7.134 12.451 ( 5.877 12.196 -0.274) 13.540 12.610 ( 7.767 3.355 -0.498) 8.475 12.793 ( 0.398 -0.852 1.727) 1.966 12.938 ( 2.746 6.035 -1.264) 6.750 13.069 ( -1.459 -0.001 -2.493) 2.888 13.204 ( 1.179 1.187 -1.509) 2.253 13.276 ( 1.480 3.991 3.339) 5.409 13.771 ( -2.199 -9.313 1.221) 9.647 ======================= Grid point 111 (30/45) ======================= q-point: ( 0.33 0.33 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.533 ( 29.375 30.490 22.635) 48.010 5.267 ( 31.094 32.667 6.945) 45.631 6.283 ( -6.163 16.968 9.919) 20.598 6.481 ( 10.517 -10.410 9.610) 17.645 6.829 ( 20.001 20.817 9.184) 30.294 7.086 ( 4.967 7.996 6.292) 11.323 7.121 ( 9.242 -3.529 7.229) 12.253 7.804 ( 1.128 4.872 -8.365) 9.746 8.041 ( -13.514 -18.255 -5.451) 23.358 8.587 ( 3.504 2.659 -6.611) 7.941 8.674 ( 3.341 -2.647 -0.273) 4.271 8.751 ( -9.141 4.658 -1.818) 10.419 8.969 ( -5.676 -5.516 -7.791) 11.105 9.223 ( -8.856 -10.207 -1.200) 13.566 9.508 ( -12.568 -7.378 -7.224) 16.265 9.521 ( -9.713 -11.829 -4.492) 15.951 10.394 ( -6.208 -1.626 -3.312) 7.222 10.429 ( -5.404 -1.689 -2.128) 6.048 10.600 ( -1.916 -6.390 -9.471) 11.585 10.671 ( -3.074 -6.930 -11.146) 13.480 10.995 ( 0.489 -1.182 5.324) 5.476 11.119 ( -1.634 2.533 5.672) 6.423 11.322 ( 7.420 -1.007 4.010) 8.494 11.362 ( 5.182 -3.268 6.539) 8.960 11.442 ( 3.698 0.489 3.219) 4.927 11.569 ( -1.905 6.859 -0.343) 7.127 11.838 ( -5.214 -9.596 2.836) 11.283 11.878 ( -5.996 -4.314 1.472) 7.531 12.568 ( 3.870 8.381 -1.348) 9.329 12.693 ( -2.667 6.296 -1.281) 6.957 12.735 ( 0.566 3.687 -1.892) 4.182 12.955 ( -5.186 4.502 2.686) 7.374 13.102 ( 9.386 -0.362 -1.451) 9.504 13.286 ( 5.263 -0.606 -2.254) 5.758 13.352 ( 3.637 1.932 0.750) 4.186 13.683 ( -5.941 -6.260 1.179) 8.711 ======================= Grid point 112 (31/45) ======================= q-point: ( 0.44 0.33 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.162 ( 30.268 27.009 17.737) 44.275 5.713 ( -8.159 26.511 11.337) 29.966 5.940 ( 4.023 30.391 9.184) 32.002 6.335 ( -19.263 23.710 7.757) 31.518 7.181 ( 17.017 -21.761 4.988) 28.071 7.457 ( 9.072 -19.608 3.545) 21.894 7.604 ( 19.207 6.064 5.169) 20.794 7.757 ( 3.420 -8.840 -1.660) 9.623 7.800 ( -1.217 -4.586 -3.097) 5.666 8.329 ( -25.156 8.689 0.793) 26.626 8.712 ( 9.381 2.980 -6.063) 11.560 8.860 ( 0.462 -1.663 -5.795) 6.046 8.941 ( 4.485 2.787 -0.964) 5.368 9.039 ( -7.599 -1.191 0.619) 7.717 9.235 ( -8.288 -10.648 -8.553) 15.976 9.310 ( -10.243 -4.711 -8.245) 13.967 10.319 ( -0.905 -3.473 -2.634) 4.452 10.343 ( -2.829 -2.579 -1.634) 4.162 10.580 ( 0.394 -8.338 -11.491) 14.202 10.613 ( -2.436 -8.380 -12.017) 14.852 11.013 ( 1.528 1.838 6.588) 7.009 11.063 ( -3.029 3.826 7.085) 8.603 11.444 ( 4.017 1.677 3.139) 5.366 11.487 ( 1.622 -4.950 5.207) 7.366 11.539 ( 6.940 1.032 0.488) 7.033 11.573 ( 4.343 0.111 1.560) 4.616 11.711 ( -6.723 -10.634 4.888) 13.497 11.729 ( -7.841 -2.013 1.435) 8.222 12.588 ( -0.334 5.626 -1.987) 5.976 12.613 ( -2.814 5.572 -1.698) 6.469 12.803 ( 3.925 5.372 -1.256) 6.771 12.874 ( -3.009 4.485 0.513) 5.425 13.289 ( 7.457 -1.198 -0.062) 7.553 13.367 ( 0.619 0.049 0.147) 0.638 13.416 ( 3.962 0.097 -1.973) 4.427 13.540 ( -7.187 -2.530 1.062) 7.693 ======================= Grid point 119 (32/45) ======================= q-point: ( 0.22 0.44 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.701 ( 21.455 36.046 19.582) 46.294 5.092 ( 38.310 28.054 10.497) 48.630 5.458 ( 19.608 -33.703 13.386) 41.226 5.698 ( 31.829 -27.217 12.121) 43.597 7.262 ( 9.252 44.902 6.494) 46.303 7.465 ( -14.535 21.769 3.517) 26.411 7.850 ( -14.083 9.687 1.665) 17.174 7.897 ( -13.686 -6.463 -1.530) 15.213 8.003 ( -8.572 -6.542 -3.338) 11.288 8.194 ( 2.995 -29.248 2.057) 29.473 8.639 ( 6.357 8.879 -6.401) 12.658 8.884 ( 0.113 -13.970 -9.631) 16.968 9.021 ( -2.663 10.521 -7.765) 13.344 9.242 ( -7.568 -3.168 -6.032) 10.183 9.337 ( -3.744 -7.847 -4.134) 9.627 9.509 ( -5.802 -13.021 -2.278) 14.436 10.484 ( -2.006 -0.055 -4.412) 4.847 10.530 ( -0.875 -3.198 -5.931) 6.794 10.630 ( -9.947 2.587 -4.638) 11.276 10.696 ( -8.886 -2.779 -6.526) 11.370 10.942 ( 1.503 -0.039 5.865) 6.055 10.999 ( 3.726 -4.739 5.204) 7.964 11.226 ( 3.748 2.638 6.562) 8.004 11.286 ( 2.035 -2.950 6.622) 7.530 11.524 ( -3.329 8.709 0.723) 9.351 11.692 ( -3.579 -4.526 -1.891) 6.072 11.805 ( -4.585 0.460 2.021) 5.032 11.877 ( -3.216 -4.202 0.778) 5.349 12.655 ( 4.952 6.931 -0.417) 8.528 12.709 ( 6.191 4.280 -2.314) 7.874 12.779 ( 7.467 0.655 -1.441) 7.633 12.897 ( 7.079 -7.092 -7.869) 12.741 13.146 ( -4.518 5.536 3.788) 8.088 13.222 ( -5.307 -1.261 2.350) 5.939 13.387 ( 2.152 7.093 1.703) 7.605 13.566 ( -0.033 -9.684 1.693) 9.831 ======================= Grid point 120 (33/45) ======================= q-point: ( 0.33 0.44 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.200 ( 25.830 32.855 16.812) 45.048 5.833 ( 28.936 3.250 6.866) 29.917 5.918 ( 25.049 -9.459 11.506) 29.143 6.342 ( 24.184 -17.950 8.144) 31.199 7.065 ( -14.304 19.592 5.199) 24.809 7.469 ( -13.477 16.849 2.758) 21.751 7.544 ( -15.242 -16.773 2.151) 22.765 7.676 ( 12.147 30.456 6.208) 33.371 7.800 ( -9.774 -8.088 -0.478) 12.695 8.314 ( 8.744 -23.614 2.626) 25.317 8.725 ( -0.588 5.778 -6.741) 8.898 8.876 ( 0.162 -3.486 -11.966) 12.465 8.969 ( -3.041 5.784 -5.454) 8.512 9.110 ( -3.721 1.360 -7.534) 8.512 9.267 ( -3.600 -8.328 -2.080) 9.308 9.355 ( -9.661 -7.555 -2.196) 12.459 10.373 ( -6.217 0.938 -5.525) 8.370 10.446 ( -6.158 -2.380 -7.524) 10.010 10.471 ( -4.935 -0.321 -6.800) 8.409 10.555 ( -5.462 -4.780 -9.872) 12.253 10.994 ( 3.391 2.381 8.643) 9.585 11.094 ( 5.045 -5.687 7.210) 10.478 11.290 ( -0.614 1.760 7.643) 7.867 11.331 ( 3.668 -1.111 7.006) 7.985 11.498 ( 1.602 3.614 1.601) 4.265 11.611 ( -3.752 -6.917 1.552) 8.021 11.719 ( -3.525 3.939 0.770) 5.342 11.805 ( -3.950 -3.621 -0.079) 5.360 12.727 ( 1.532 6.421 -2.452) 7.042 12.795 ( 1.923 0.522 -1.993) 2.818 12.826 ( -2.287 6.114 -4.265) 7.798 12.946 ( -2.845 -4.929 -5.121) 7.656 13.153 ( 6.518 4.712 3.456) 8.754 13.240 ( 5.993 -3.307 0.381) 6.855 13.419 ( 0.719 4.446 1.442) 4.729 13.540 ( -2.389 -6.826 1.658) 7.419 ======================= Grid point 121 (34/45) ======================= q-point: ( 0.44 0.44 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.754 ( 26.546 29.056 13.410) 41.579 6.235 ( -11.977 22.545 9.408) 27.207 6.323 ( 18.850 -19.810 8.476) 28.628 6.462 ( 3.874 8.485 7.922) 12.237 7.067 ( 0.402 -7.231 3.277) 7.949 7.153 ( -12.907 15.197 3.836) 20.304 7.267 ( 3.532 -19.488 1.734) 19.881 7.537 ( -13.492 -16.503 -0.014) 21.317 8.143 ( 23.674 22.693 4.650) 33.122 8.549 ( -2.118 0.917 -0.063) 2.309 8.578 ( 2.099 -4.360 -0.208) 4.844 8.836 ( -8.495 -11.635 -2.327) 14.593 8.931 ( 3.277 9.127 -8.432) 12.851 9.091 ( -9.754 0.884 -3.130) 10.282 9.116 ( 0.575 -0.294 -7.128) 7.158 9.194 ( -4.326 -6.642 -5.133) 9.443 10.296 ( -0.625 1.985 -6.415) 6.745 10.340 ( -4.269 2.804 -6.565) 8.318 10.406 ( -0.186 -7.499 -10.060) 12.549 10.452 ( -4.045 -6.881 -10.750) 13.389 11.080 ( 4.714 5.314 9.138) 11.574 11.165 ( -1.347 4.884 7.869) 9.359 11.270 ( 2.735 -8.790 9.903) 13.521 11.358 ( -2.696 -8.918 8.468) 12.590 11.578 ( 3.847 0.289 0.556) 3.897 11.582 ( 3.327 -0.118 0.905) 3.450 11.648 ( -3.503 3.566 0.163) 5.001 11.703 ( -5.725 -1.602 0.009) 5.945 12.731 ( -0.177 7.154 -3.412) 7.928 12.751 ( -2.195 7.082 -3.603) 8.243 12.843 ( 1.576 -2.137 -3.460) 4.362 12.883 ( -2.364 -3.940 -3.000) 5.488 13.306 ( 6.429 2.957 1.476) 7.229 13.367 ( 4.440 -1.295 0.499) 4.652 13.409 ( -1.260 1.059 0.698) 1.787 13.473 ( -3.850 -3.702 1.329) 5.504 ======================= Grid point 182 (35/45) ======================= q-point: ( 0.22 0.22 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 125 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.127 ( 29.317 29.317 29.317) 50.778 4.127 ( 29.317 29.317 29.317) 50.778 5.404 ( 37.979 37.979 37.979) 65.781 6.641 ( 0.498 0.498 0.498) 0.863 7.106 ( 4.253 4.253 4.253) 7.366 7.106 ( 4.253 4.253 4.253) 7.366 7.189 ( 5.939 5.939 5.939) 10.287 7.444 ( 0.764 0.764 0.764) 1.323 7.444 ( 0.764 0.764 0.764) 1.323 8.759 ( -4.121 -4.121 -4.121) 7.138 8.759 ( -4.121 -4.121 -4.121) 7.138 8.929 ( -9.665 -9.665 -9.665) 16.741 9.178 ( -9.874 -9.874 -9.874) 17.103 9.696 ( -12.266 -12.266 -12.266) 21.246 9.696 ( -12.266 -12.266 -12.266) 21.246 9.755 ( -6.627 -6.627 -6.627) 11.478 10.527 ( -3.423 -3.423 -3.423) 5.929 10.527 ( -3.423 -3.423 -3.423) 5.929 10.651 ( -6.201 -6.201 -6.201) 10.741 10.651 ( -6.201 -6.201 -6.201) 10.741 11.019 ( 1.859 1.859 1.859) 3.220 11.214 ( 3.847 3.847 3.847) 6.663 11.264 ( 2.836 2.836 2.836) 4.912 11.264 ( 2.836 2.836 2.836) 4.912 11.491 ( 0.914 0.914 0.914) 1.584 11.491 ( 0.914 0.914 0.914) 1.584 12.058 ( -2.508 -2.508 -2.508) 4.343 12.058 ( -2.508 -2.508 -2.508) 4.343 12.266 ( 7.731 7.731 7.731) 13.390 12.720 ( 2.804 2.804 2.804) 4.857 12.720 ( 2.804 2.804 2.804) 4.857 12.819 ( -0.905 -0.905 -0.905) 1.567 13.106 ( -0.137 -0.137 -0.137) 0.237 13.106 ( -0.137 -0.137 -0.137) 0.237 13.265 ( 3.562 3.562 3.562) 6.170 13.895 ( -3.965 -3.965 -3.965) 6.867 ======================= Grid point 183 (36/45) ======================= q-point: ( 0.33 0.22 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.609 ( 25.588 28.702 28.886) 48.093 4.843 ( 30.920 25.908 26.524) 48.278 6.209 ( 12.515 20.132 24.274) 33.929 6.444 ( -4.001 16.286 12.900) 21.157 6.815 ( -7.288 18.083 22.995) 30.147 7.240 ( 11.540 6.663 -12.518) 18.283 7.460 ( 14.977 -10.839 0.905) 18.510 7.699 ( 11.259 -11.047 0.878) 15.798 7.727 ( 19.322 -0.132 -11.884) 22.685 8.568 ( -12.588 -5.386 -4.483) 14.407 8.627 ( -10.597 -0.254 -4.297) 11.438 8.787 ( -0.245 -8.417 -2.255) 8.718 8.994 ( -8.385 -4.641 -5.806) 11.206 9.335 ( -17.153 -6.161 -11.752) 21.686 9.457 ( -13.072 -10.189 -4.514) 17.177 9.615 ( -8.429 -7.738 -7.390) 13.621 10.438 ( -3.031 -5.535 -2.667) 6.851 10.448 ( -5.504 -4.490 -3.577) 7.953 10.519 ( -2.783 -6.737 -8.982) 11.567 10.590 ( -2.798 -8.528 -10.419) 13.752 11.040 ( 0.213 3.071 3.058) 4.340 11.208 ( -2.960 4.806 6.758) 8.805 11.316 ( 2.695 6.126 0.979) 6.764 11.386 ( 6.956 1.388 1.864) 7.334 11.469 ( -1.198 4.781 2.660) 5.601 11.578 ( 5.297 -2.162 0.066) 5.721 11.941 ( -6.587 -0.363 -0.955) 6.665 11.999 ( -4.492 -0.580 -2.248) 5.056 12.466 ( 9.691 2.207 3.865) 10.664 12.629 ( -2.965 2.083 3.501) 5.038 12.786 ( -1.854 -6.609 5.540) 8.821 12.852 ( 0.224 0.524 -6.592) 6.616 13.099 ( 4.928 2.903 3.591) 6.753 13.208 ( 5.631 1.969 -2.525) 6.478 13.335 ( 2.800 2.387 2.217) 4.296 13.767 ( -7.897 -1.968 -2.062) 8.395 ======================= Grid point 184 (37/45) ======================= q-point: ( 0.44 0.22 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.174 ( 27.852 22.714 24.611) 43.559 5.451 ( 25.641 24.049 25.277) 43.299 5.798 ( -24.346 18.481 21.613) 37.435 6.079 ( -27.154 23.116 22.207) 42.010 7.469 ( 38.647 13.157 9.518) 41.920 7.544 ( 15.826 2.212 -12.700) 20.412 7.732 ( 9.952 -11.085 -0.994) 14.930 7.824 ( -0.404 -10.587 -3.463) 11.147 7.917 ( -2.885 -4.502 -12.263) 13.378 8.250 ( -26.013 3.852 3.531) 26.533 8.494 ( 9.201 -4.865 -2.464) 10.695 8.774 ( -11.599 -0.170 -3.553) 12.132 8.895 ( 5.461 -3.736 -1.923) 6.890 9.066 ( -8.565 -7.701 -8.466) 14.295 9.226 ( -2.870 -7.426 -2.684) 8.401 9.390 ( -12.333 -8.367 -4.640) 15.609 10.376 ( -0.494 -5.875 -4.082) 7.171 10.400 ( -1.856 -5.927 -2.329) 6.633 10.492 ( 0.660 -7.796 -11.778) 14.140 10.534 ( -2.632 -8.577 -11.266) 14.402 11.055 ( 1.891 4.910 4.647) 7.020 11.125 ( -4.349 5.431 6.222) 9.335 11.377 ( 3.038 7.097 2.118) 8.005 11.434 ( -2.262 6.612 2.807) 7.531 11.539 ( 7.404 0.803 -0.569) 7.469 11.653 ( -0.889 -1.664 -1.104) 2.186 11.831 ( -1.424 -0.208 -1.741) 2.258 11.893 ( -5.071 -0.148 -2.264) 5.556 12.590 ( 3.492 -4.684 3.941) 7.047 12.663 ( 1.041 -6.487 5.068) 8.297 12.725 ( 0.114 0.206 2.756) 2.766 12.807 ( -3.624 1.005 -0.272) 3.770 13.230 ( 5.607 3.709 -1.783) 6.955 13.294 ( 0.368 3.959 -3.470) 5.277 13.412 ( 6.378 0.644 2.243) 6.791 13.584 ( -9.105 -0.204 -0.045) 9.108 ======================= Grid point 192 (38/45) ======================= q-point: ( 0.33 0.33 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.129 ( 23.865 24.140 31.938) 46.608 5.476 ( 29.528 30.143 12.718) 44.071 6.537 ( -3.148 2.956 12.972) 13.672 6.784 ( 8.441 6.140 19.149) 21.809 7.057 ( 10.876 2.735 11.022) 15.724 7.309 ( -10.557 11.188 11.586) 19.257 7.316 ( 22.959 1.569 13.975) 26.924 7.559 ( 5.087 9.151 -14.239) 17.674 7.842 ( -0.052 -2.661 -13.375) 13.638 8.423 ( -4.123 -2.399 -7.777) 9.123 8.638 ( 0.340 -6.095 -4.245) 7.435 8.691 ( -5.024 3.138 -3.889) 7.086 8.850 ( -5.645 -7.014 -2.970) 9.481 9.132 ( -7.647 -13.896 -6.493) 17.139 9.291 ( -11.303 -4.985 -12.121) 17.307 9.469 ( -7.956 -6.671 -0.314) 10.387 10.311 ( -6.210 -5.805 -4.586) 9.659 10.356 ( -4.956 -4.923 -5.620) 8.966 10.389 ( -4.380 -5.139 -9.125) 11.352 10.436 ( -4.899 -5.358 -10.730) 12.955 11.097 ( 2.421 2.255 4.603) 5.669 11.278 ( -0.944 2.455 8.430) 8.831 11.388 ( 6.494 -0.446 2.259) 6.891 11.443 ( 3.833 5.002 1.688) 6.523 11.533 ( 3.459 -1.692 4.076) 5.608 11.568 ( -0.658 4.426 0.642) 4.520 11.888 ( -4.383 -5.133 2.114) 7.073 11.929 ( -4.327 -5.526 1.901) 7.272 12.542 ( 2.364 4.249 -0.972) 4.959 12.651 ( -3.969 0.343 -2.524) 4.716 12.718 ( -1.597 1.571 -0.964) 2.439 12.856 ( 0.299 -1.772 -8.910) 9.089 13.174 ( 5.000 4.372 6.300) 9.155 13.263 ( 6.225 2.714 0.548) 6.813 13.381 ( 1.576 2.188 1.903) 3.300 13.688 ( -5.263 -5.261 -0.922) 7.498 ======================= Grid point 193 (39/45) ======================= q-point: ( 0.44 0.33 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.638 ( 24.516 22.579 25.800) 42.148 5.998 ( 12.368 25.681 15.212) 32.309 6.226 ( -10.411 21.845 17.205) 29.692 6.571 ( -21.046 20.896 14.593) 33.053 7.305 ( 15.075 -18.392 5.943) 24.512 7.545 ( 7.765 -16.810 3.849) 18.913 7.642 ( 1.622 -0.573 -9.046) 9.208 7.659 ( -8.679 -14.209 -1.842) 16.751 7.837 ( 25.092 16.874 7.321) 31.112 8.334 ( -19.229 4.097 -1.441) 19.714 8.537 ( 12.173 7.451 -7.544) 16.143 8.742 ( -0.213 -4.843 -4.994) 6.960 8.882 ( 4.446 -0.057 -8.988) 10.028 8.966 ( -9.819 -1.901 -10.677) 14.629 9.051 ( -3.696 -7.235 0.199) 8.126 9.242 ( -13.018 -3.898 2.652) 13.845 10.244 ( -0.086 -6.252 -4.973) 7.989 10.285 ( -3.120 -5.288 -4.295) 7.493 10.321 ( -0.470 -8.308 -12.311) 14.859 10.355 ( -2.993 -7.919 -11.448) 14.238 11.153 ( 3.184 4.700 6.494) 8.625 11.227 ( -3.497 5.036 8.024) 10.098 11.492 ( 3.286 3.044 1.543) 4.738 11.536 ( 0.281 3.795 1.294) 4.019 11.567 ( 4.378 1.061 1.551) 4.764 11.615 ( 3.012 -2.157 1.700) 4.076 11.782 ( -3.875 -3.899 3.160) 6.341 11.823 ( -5.474 -6.883 4.216) 9.752 12.543 ( -0.860 0.999 -2.566) 2.885 12.562 ( -2.600 -0.609 -3.068) 4.067 12.761 ( 3.472 1.002 -3.451) 4.997 12.827 ( -2.506 -0.480 -5.403) 5.975 13.296 ( 5.437 2.976 2.108) 6.547 13.369 ( -0.480 3.138 -0.323) 3.191 13.419 ( 5.777 0.139 2.602) 6.337 13.560 ( -6.743 -1.936 0.726) 7.053 ======================= Grid point 201 (40/45) ======================= q-point: ( 0.33 0.44 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.653 ( 22.471 25.775 24.733) 42.202 6.036 ( 27.384 18.033 12.453) 35.074 6.232 ( 19.799 -16.616 17.190) 31.042 6.585 ( 21.640 -21.543 14.539) 33.819 7.207 ( -14.689 15.886 7.703) 22.967 7.539 ( -11.208 1.594 2.965) 11.702 7.583 ( -15.005 -2.659 1.109) 15.279 7.706 ( -6.700 -0.381 -6.764) 9.528 7.876 ( 15.433 27.623 8.624) 32.796 8.353 ( 6.493 -16.124 -0.937) 17.407 8.532 ( 5.066 5.722 -10.059) 12.633 8.674 ( -5.935 -4.929 -3.655) 8.537 8.866 ( 0.560 0.334 -4.597) 4.643 8.882 ( -7.517 -4.175 -11.939) 14.713 9.196 ( -2.486 -1.257 -5.028) 5.749 9.310 ( -4.143 -8.383 -0.694) 9.376 10.235 ( -7.435 -0.833 -7.215) 10.394 10.266 ( -7.773 -3.220 -8.960) 12.291 10.325 ( -5.839 -0.771 -7.294) 9.375 10.353 ( -6.958 -2.933 -8.959) 11.717 11.155 ( 5.035 3.877 6.644) 9.194 11.256 ( 6.232 -4.757 7.299) 10.712 11.425 ( 1.940 4.250 5.379) 7.124 11.494 ( 2.451 -1.729 7.079) 7.688 11.538 ( 1.196 3.954 2.003) 4.591 11.628 ( -1.369 -1.538 0.216) 2.070 11.768 ( -3.930 -2.023 3.300) 5.517 11.823 ( -4.295 -4.524 2.024) 6.558 12.636 ( 0.289 5.336 -7.579) 9.273 12.702 ( -2.349 2.665 -6.360) 7.285 12.744 ( -2.811 -0.235 -6.351) 6.949 12.775 ( 0.908 -4.837 -5.679) 7.515 13.262 ( 4.055 3.300 5.918) 7.896 13.297 ( 4.075 0.041 4.442) 6.028 13.446 ( 0.788 4.555 1.223) 4.782 13.564 ( -1.856 -6.219 0.587) 6.516 ======================= Grid point 202 (41/45) ======================= q-point: ( 0.44 0.44 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 6.125 ( 22.300 23.712 20.862) 38.662 6.488 ( 0.406 20.723 13.590) 24.785 6.550 ( 13.905 -6.722 12.889) 20.116 6.706 ( -1.268 -0.010 14.887) 14.941 7.160 ( 1.741 -4.725 5.334) 7.336 7.254 ( -11.825 6.298 5.238) 14.385 7.313 ( -1.134 -16.387 2.515) 16.617 7.524 ( -10.853 -10.943 -1.616) 15.497 8.220 ( 14.429 17.237 -0.736) 22.491 8.456 ( -6.984 6.971 -9.286) 13.550 8.541 ( 6.370 -10.432 -5.729) 13.499 8.696 ( -2.941 -11.908 -12.318) 17.384 8.779 ( 8.097 11.693 -2.600) 14.459 8.972 ( 8.045 -3.205 -0.604) 8.681 9.062 ( -10.839 5.589 0.137) 12.196 9.172 ( -8.754 -4.363 4.274) 10.674 10.146 ( -0.485 -2.187 -7.626) 7.948 10.172 ( -0.916 -4.416 -10.507) 11.434 10.203 ( -5.043 -2.388 -7.733) 9.536 10.223 ( -4.992 -4.214 -10.265) 12.168 11.267 ( 5.475 6.553 8.222) 11.854 11.348 ( -0.166 5.316 8.844) 10.320 11.452 ( 3.043 -5.580 6.450) 9.055 11.524 ( -0.875 -4.365 6.552) 7.921 11.599 ( 3.568 0.714 1.540) 3.952 11.606 ( 0.947 -0.116 1.220) 1.548 11.680 ( -3.809 -1.560 3.571) 5.449 11.718 ( -5.395 -2.696 1.814) 6.298 12.619 ( -1.054 5.391 -8.903) 10.462 12.630 ( -2.205 5.917 -6.572) 9.114 12.730 ( 0.704 -3.988 -8.556) 9.466 12.762 ( -1.736 -5.464 -7.589) 9.511 13.363 ( 4.735 3.093 4.067) 6.966 13.410 ( 5.537 -1.230 3.121) 6.474 13.435 ( -1.993 3.411 2.257) 4.550 13.507 ( -3.657 -3.112 1.903) 5.166 ======================= Grid point 273 (42/45) ======================= q-point: ( 0.33 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 125 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.787 ( 23.393 23.393 23.393) 40.518 5.787 ( 23.393 23.393 23.393) 40.518 6.753 ( 4.102 4.102 4.102) 7.105 7.193 ( 9.679 9.679 9.679) 16.764 7.261 ( 1.322 1.322 1.322) 2.290 7.261 ( 1.322 1.322 1.322) 2.290 7.507 ( -0.501 -0.501 -0.501) 0.869 7.507 ( -0.501 -0.501 -0.501) 0.869 7.716 ( 13.691 13.691 13.691) 23.713 8.294 ( -2.014 -2.014 -2.014) 3.488 8.581 ( -4.008 -4.008 -4.008) 6.943 8.581 ( -4.008 -4.008 -4.008) 6.943 8.724 ( -8.677 -8.677 -8.677) 15.028 9.060 ( -4.216 -4.216 -4.216) 7.303 9.060 ( -4.216 -4.216 -4.216) 7.303 9.467 ( -1.806 -1.806 -1.806) 3.128 10.185 ( -7.378 -7.378 -7.378) 12.779 10.185 ( -7.378 -7.378 -7.378) 12.779 10.262 ( -5.366 -5.366 -5.366) 9.295 10.262 ( -5.366 -5.366 -5.366) 9.295 11.190 ( 4.353 4.353 4.353) 7.540 11.433 ( 2.672 2.672 2.672) 4.628 11.433 ( 2.672 2.672 2.672) 4.628 11.465 ( 3.991 3.991 3.991) 6.913 11.594 ( 1.713 1.713 1.713) 2.966 11.594 ( 1.713 1.713 1.713) 2.966 11.913 ( -1.991 -1.991 -1.991) 3.448 11.913 ( -1.991 -1.991 -1.991) 3.448 12.528 ( -0.971 -0.971 -0.971) 1.682 12.596 ( -5.130 -5.130 -5.130) 8.885 12.596 ( -5.130 -5.130 -5.130) 8.885 12.760 ( -1.463 -1.463 -1.463) 2.534 13.308 ( 4.751 4.751 4.751) 8.229 13.308 ( 4.751 4.751 4.751) 8.229 13.422 ( 1.830 1.830 1.830) 3.170 13.645 ( -3.062 -3.062 -3.062) 5.303 ======================= Grid point 274 (43/45) ======================= q-point: ( 0.44 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 6.162 ( 18.615 22.627 22.949) 37.218 6.353 ( 20.956 18.036 19.106) 33.608 6.594 ( -11.765 16.792 16.854) 26.542 6.912 ( -22.364 15.207 16.258) 31.555 7.336 ( 10.961 0.804 -7.374) 13.235 7.476 ( 6.191 -6.805 -2.206) 9.460 7.545 ( -2.830 -6.952 -5.210) 9.137 7.575 ( -0.439 -4.953 -6.951) 8.546 7.944 ( 13.444 -0.637 2.323) 13.658 8.275 ( -14.562 -4.140 -2.773) 15.391 8.492 ( 8.253 0.717 -0.885) 8.332 8.578 ( -1.992 -9.923 -16.046) 18.972 8.658 ( 7.147 -0.419 2.144) 7.473 8.908 ( -9.112 8.237 5.498) 13.458 9.122 ( 6.810 -0.847 3.710) 7.801 9.330 ( -10.742 1.775 3.747) 11.515 10.079 ( -1.798 -7.042 -9.021) 11.584 10.089 ( -2.328 -7.999 -9.038) 12.292 10.165 ( -2.019 -7.531 -6.569) 10.196 10.185 ( -2.918 -7.441 -6.589) 10.358 11.287 ( 5.062 6.432 6.613) 10.523 11.389 ( -4.100 5.961 7.532) 10.444 11.520 ( 4.679 3.438 1.334) 5.957 11.571 ( 2.170 2.170 2.243) 3.801 11.605 ( 0.816 2.971 1.769) 3.552 11.638 ( 0.945 0.591 0.502) 1.222 11.818 ( -2.714 0.613 -0.275) 2.796 11.861 ( -4.205 0.028 -0.765) 4.274 12.467 ( -1.435 -6.050 -4.536) 7.696 12.496 ( -3.238 -9.649 -2.652) 10.518 12.621 ( 3.778 -5.562 -9.020) 11.251 12.708 ( -3.727 -4.168 -5.431) 7.795 13.377 ( 2.648 4.921 4.777) 7.352 13.399 ( 1.264 4.822 3.149) 5.896 13.473 ( 4.048 1.307 2.102) 4.745 13.567 ( -4.276 -0.241 -0.090) 4.284 ======================= Grid point 283 (44/45) ======================= q-point: ( 0.44 0.44 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 365 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 6.582 ( 17.639 18.029 22.061) 33.509 6.764 ( 16.301 18.894 12.843) 28.065 6.855 ( 3.844 1.769 15.435) 16.004 7.042 ( -4.990 -4.307 16.089) 17.387 7.273 ( 3.141 -3.288 4.788) 6.603 7.344 ( -5.641 -0.395 2.680) 6.257 7.370 ( -8.846 -11.232 3.021) 14.613 7.460 ( -4.143 -2.592 -4.439) 6.602 8.033 ( 6.931 8.679 -14.318) 18.121 8.215 ( -6.856 -3.100 -12.311) 14.428 8.306 ( -4.188 -7.968 -14.656) 17.200 8.374 ( -6.689 -8.862 -17.062) 20.356 8.922 ( 15.471 15.038 13.840) 25.633 9.082 ( 6.751 2.188 9.335) 11.727 9.155 ( 0.551 6.006 6.934) 9.190 9.313 ( -6.604 -5.096 7.962) 11.531 9.961 ( -1.553 -2.965 -8.993) 9.596 9.972 ( -1.491 -2.965 -9.061) 9.649 10.028 ( -5.191 -3.925 -7.483) 9.917 10.052 ( -5.588 -4.883 -7.842) 10.797 11.431 ( 6.725 7.104 8.085) 12.691 11.519 ( 0.423 4.608 7.266) 8.614 11.548 ( 4.172 0.464 3.730) 5.616 11.621 ( -1.032 -0.274 2.963) 3.150 11.637 ( 2.928 1.411 2.376) 4.026 11.646 ( 1.862 0.532 2.270) 2.984 11.762 ( -3.738 -4.117 3.817) 6.745 11.797 ( -4.743 -4.937 4.850) 8.390 12.396 ( -0.832 0.934 -11.368) 11.437 12.454 ( -6.028 4.269 -10.102) 12.514 12.489 ( 2.565 -6.483 -12.893) 14.657 12.586 ( -5.898 -6.992 -9.559) 13.231 13.459 ( 2.892 2.705 4.677) 6.128 13.482 ( 3.216 1.574 3.491) 5.001 13.497 ( 0.653 2.080 3.119) 3.806 13.549 ( -1.849 -1.705 1.886) 3.143 ======================= Grid point 364 (45/45) ======================= q-point: ( 0.44 0.44 0.44) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.65e-04 8.24e-05 8.24e-05 8.24e-05 8.24e-05 Number of triplets: 125 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 7.021 ( 15.192 15.192 15.192) 26.314 7.021 ( 15.192 15.192 15.192) 26.314 7.129 ( 7.166 7.166 7.166) 12.411 7.297 ( -0.905 -0.905 -0.905) 1.567 7.321 ( 0.627 0.627 0.627) 1.085 7.321 ( 0.627 0.627 0.627) 1.085 7.401 ( -1.375 -1.375 -1.375) 2.382 7.401 ( -1.375 -1.375 -1.375) 2.382 7.807 ( -1.059 -1.059 -1.059) 1.835 7.990 ( -10.022 -10.022 -10.022) 17.358 8.032 ( -10.588 -10.588 -10.588) 18.339 8.032 ( -10.588 -10.588 -10.588) 18.339 9.226 ( 11.982 11.982 11.982) 20.753 9.286 ( 6.534 6.534 6.534) 11.318 9.286 ( 6.534 6.534 6.534) 11.318 9.425 ( -0.229 -0.229 -0.229) 0.397 9.821 ( -3.399 -3.399 -3.399) 5.887 9.821 ( -3.399 -3.399 -3.399) 5.887 9.890 ( -5.972 -5.972 -5.972) 10.343 9.890 ( -5.972 -5.972 -5.972) 10.343 11.620 ( 8.956 8.956 8.956) 15.513 11.620 ( 2.832 2.832 2.832) 4.905 11.620 ( 2.832 2.832 2.832) 4.905 11.664 ( 0.487 0.487 0.487) 0.844 11.664 ( 0.487 0.487 0.487) 0.844 11.735 ( 4.683 4.683 4.683) 8.112 11.834 ( -0.270 -0.270 -0.270) 0.468 11.834 ( -0.270 -0.270 -0.270) 0.468 12.237 ( -5.085 -5.085 -5.085) 8.807 12.237 ( -5.085 -5.085 -5.085) 8.807 12.265 ( -5.734 -5.734 -5.734) 9.932 12.376 ( -9.738 -9.738 -9.738) 16.867 13.540 ( 2.097 2.097 2.097) 3.633 13.540 ( 2.097 2.097 2.097) 3.633 13.547 ( 1.585 1.585 1.585) 2.746 13.572 ( 0.002 0.002 0.002) 0.004 =================== 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/26244 10.0 6067.376 6067.376 6067.376 0.000 0.000 0.000 3/26244 20.0 9598.858 9598.858 9598.858 0.000 0.000 0.000 3/26244 30.0 5400.760 5400.760 5400.760 0.000 0.000 0.000 3/26244 40.0 2898.448 2898.448 2898.448 0.000 0.000 0.000 3/26244 50.0 1547.158 1547.158 1547.158 0.000 0.000 0.000 3/26244 60.0 883.321 883.321 883.321 0.000 0.000 0.000 3/26244 70.0 558.049 558.049 558.049 0.000 0.000 0.000 3/26244 80.0 387.548 387.548 387.548 0.000 0.000 0.000 3/26244 90.0 290.381 290.381 290.381 0.000 0.000 0.000 3/26244 100.0 230.517 230.517 230.517 0.000 0.000 0.000 3/26244 110.0 191.080 191.080 191.080 0.000 0.000 0.000 3/26244 120.0 163.602 163.602 163.602 0.000 0.000 0.000 3/26244 130.0 143.544 143.544 143.544 0.000 0.000 0.000 3/26244 140.0 128.323 128.323 128.323 0.000 0.000 0.000 3/26244 150.0 116.393 116.393 116.393 0.000 0.000 0.000 3/26244 160.0 106.787 106.787 106.787 0.000 0.000 0.000 3/26244 170.0 98.875 98.875 98.875 0.000 0.000 0.000 3/26244 180.0 92.231 92.231 92.231 0.000 0.000 0.000 3/26244 190.0 86.561 86.561 86.561 0.000 0.000 0.000 3/26244 200.0 81.655 81.655 81.655 0.000 0.000 0.000 3/26244 210.0 77.359 77.359 77.359 0.000 0.000 0.000 3/26244 220.0 73.558 73.558 73.558 0.000 0.000 0.000 3/26244 230.0 70.165 70.165 70.165 0.000 0.000 0.000 3/26244 240.0 67.114 67.114 67.114 0.000 0.000 0.000 3/26244 250.0 64.350 64.350 64.350 0.000 0.000 0.000 3/26244 260.0 61.833 61.833 61.833 0.000 0.000 0.000 3/26244 270.0 59.528 59.528 59.528 0.000 0.000 0.000 3/26244 280.0 57.407 57.407 57.407 0.000 0.000 0.000 3/26244 290.0 55.447 55.447 55.447 0.000 0.000 0.000 3/26244 300.0 53.629 53.629 53.629 0.000 0.000 0.000 3/26244 310.0 51.937 51.937 51.937 0.000 0.000 0.000 3/26244 320.0 50.358 50.358 50.358 0.000 0.000 0.000 3/26244 330.0 48.880 48.880 48.880 0.000 0.000 0.000 3/26244 340.0 47.492 47.492 47.492 0.000 0.000 0.000 3/26244 350.0 46.187 46.187 46.187 0.000 0.000 0.000 3/26244 360.0 44.956 44.956 44.956 0.000 0.000 0.000 3/26244 370.0 43.793 43.793 43.793 0.000 0.000 0.000 3/26244 380.0 42.692 42.692 42.692 0.000 0.000 0.000 3/26244 390.0 41.649 41.649 41.649 0.000 0.000 0.000 3/26244 400.0 40.658 40.658 40.658 0.000 0.000 0.000 3/26244 410.0 39.715 39.715 39.715 0.000 0.000 0.000 3/26244 420.0 38.817 38.817 38.817 0.000 0.000 0.000 3/26244 430.0 37.961 37.961 37.961 0.000 0.000 0.000 3/26244 440.0 37.143 37.143 37.143 0.000 0.000 0.000 3/26244 450.0 36.362 36.362 36.362 0.000 0.000 0.000 3/26244 460.0 35.613 35.613 35.613 0.000 0.000 0.000 3/26244 470.0 34.896 34.896 34.896 0.000 0.000 0.000 3/26244 480.0 34.209 34.209 34.209 0.000 0.000 0.000 3/26244 490.0 33.549 33.549 33.549 0.000 0.000 0.000 3/26244 500.0 32.914 32.914 32.914 0.000 0.000 0.000 3/26244 510.0 32.304 32.304 32.304 0.000 0.000 0.000 3/26244 520.0 31.717 31.717 31.717 0.000 0.000 0.000 3/26244 530.0 31.152 31.152 31.152 0.000 0.000 0.000 3/26244 540.0 30.607 30.607 30.607 0.000 0.000 0.000 3/26244 550.0 30.081 30.081 30.081 0.000 0.000 0.000 3/26244 560.0 29.573 29.573 29.573 0.000 0.000 0.000 3/26244 570.0 29.083 29.083 29.083 0.000 0.000 0.000 3/26244 580.0 28.609 28.609 28.609 0.000 0.000 0.000 3/26244 590.0 28.150 28.150 28.150 0.000 0.000 0.000 3/26244 600.0 27.706 27.706 27.706 0.000 0.000 0.000 3/26244 610.0 27.277 27.277 27.277 0.000 0.000 0.000 3/26244 620.0 26.860 26.860 26.860 0.000 0.000 0.000 3/26244 630.0 26.457 26.457 26.457 0.000 0.000 0.000 3/26244 640.0 26.066 26.066 26.066 0.000 0.000 0.000 3/26244 650.0 25.686 25.686 25.686 0.000 0.000 0.000 3/26244 660.0 25.317 25.317 25.317 0.000 0.000 0.000 3/26244 670.0 24.959 24.959 24.959 0.000 0.000 0.000 3/26244 680.0 24.612 24.612 24.612 0.000 0.000 0.000 3/26244 690.0 24.273 24.273 24.273 0.000 0.000 0.000 3/26244 700.0 23.945 23.945 23.945 0.000 0.000 0.000 3/26244 710.0 23.625 23.625 23.625 0.000 0.000 0.000 3/26244 720.0 23.313 23.313 23.313 0.000 0.000 0.000 3/26244 730.0 23.010 23.010 23.010 0.000 0.000 0.000 3/26244 740.0 22.715 22.715 22.715 0.000 0.000 0.000 3/26244 750.0 22.427 22.427 22.427 0.000 0.000 0.000 3/26244 760.0 22.147 22.147 22.147 0.000 0.000 0.000 3/26244 770.0 21.874 21.874 21.874 0.000 0.000 0.000 3/26244 780.0 21.607 21.607 21.607 0.000 0.000 0.000 3/26244 790.0 21.347 21.347 21.347 0.000 0.000 0.000 3/26244 800.0 21.093 21.093 21.093 0.000 0.000 0.000 3/26244 810.0 20.845 20.845 20.845 0.000 0.000 0.000 3/26244 820.0 20.603 20.603 20.603 0.000 0.000 0.000 3/26244 830.0 20.367 20.367 20.367 0.000 0.000 0.000 3/26244 840.0 20.136 20.136 20.136 0.000 0.000 0.000 3/26244 850.0 19.910 19.910 19.910 0.000 0.000 0.000 3/26244 860.0 19.690 19.690 19.690 0.000 0.000 0.000 3/26244 870.0 19.474 19.474 19.474 0.000 0.000 0.000 3/26244 880.0 19.263 19.263 19.263 0.000 0.000 0.000 3/26244 890.0 19.056 19.056 19.056 0.000 0.000 0.000 3/26244 900.0 18.854 18.854 18.854 0.000 0.000 0.000 3/26244 910.0 18.657 18.657 18.657 0.000 0.000 0.000 3/26244 920.0 18.463 18.463 18.463 0.000 0.000 0.000 3/26244 930.0 18.273 18.273 18.273 0.000 0.000 0.000 3/26244 940.0 18.088 18.088 18.088 0.000 0.000 0.000 3/26244 950.0 17.906 17.906 17.906 0.000 0.000 0.000 3/26244 960.0 17.727 17.727 17.727 0.000 0.000 0.000 3/26244 970.0 17.552 17.552 17.552 0.000 0.000 0.000 3/26244 980.0 17.381 17.381 17.381 0.000 0.000 0.000 3/26244 990.0 17.213 17.213 17.213 0.000 0.000 0.000 3/26244 1000.0 17.048 17.048 17.048 0.000 0.000 0.000 3/26244 Thermal conductivity related properties were written into "kappa-m999.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 18:31:07]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|