----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 00:18:04]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [1 1 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I-43d (220) Number of symmetry operations in supercell: 48 ------------------------------ primitive cell ------------------------------ Lattice vectors: a -4.895274240000000 4.895274240000000 4.895274240000000 b 4.895274240000000 -4.895274240000000 4.895274240000000 c 4.895274240000000 4.895274240000000 -4.895274240000000 Atomic positions (fractional): *1 Y 0.37500000000000 0.62500000000000 0.75000000000000 88.906 2 Y 0.87500000000000 0.25000000000000 0.12500000000000 88.906 3 Y 0.62500000000000 0.75000000000000 0.37500000000000 88.906 4 Y 0.25000000000000 0.12500000000000 0.87500000000000 88.906 5 Y 0.75000000000000 0.37500000000000 0.62500000000000 88.906 6 Y 0.12500000000000 0.87500000000000 0.25000000000000 88.906 *7 Sb 0.17733449843421 0.50000000000000 0.00000000000000 121.760 8 Sb 0.50000000000000 0.00000000000000 0.17733449843421 121.760 9 Sb 0.32266550156579 0.32266550156579 0.32266550156579 121.760 10 Sb 0.50000000000000 0.00000000000000 0.67733449843421 121.760 11 Sb 0.00000000000000 0.67733449843421 0.50000000000000 121.760 12 Sb 0.82266550156579 0.82266550156579 0.82266550156579 121.760 13 Sb 0.00000000000000 0.17733449843421 0.50000000000000 121.760 14 Sb 0.67733449843421 0.50000000000000 0.00000000000000 121.760 *15 Au 0.62500000000000 0.37500000000000 0.25000000000000 196.967 16 Au 0.25000000000000 0.62500000000000 0.37500000000000 196.967 17 Au 0.12500000000000 0.75000000000000 0.87500000000000 196.967 18 Au 0.37500000000000 0.25000000000000 0.62500000000000 196.967 19 Au 0.87500000000000 0.12500000000000 0.75000000000000 196.967 20 Au 0.75000000000000 0.87500000000000 0.12500000000000 196.967 -------------------------------- unit cell --------------------------------- Lattice vectors: a 9.790548480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 9.790548480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 9.790548480000000 Atomic positions (fractional): *1 Y 0.50000000000000 0.25000000000000 0.12500000000000 88.906 > 1 2 Y 0.75000000000000 0.37500000000000 0.50000000000000 88.906 > 2 3 Y 0.25000000000000 0.12500000000000 0.50000000000000 88.906 > 3 4 Y 0.37500000000000 0.50000000000000 0.75000000000000 88.906 > 4 5 Y 0.62500000000000 0.00000000000000 0.75000000000000 88.906 > 5 6 Y 0.00000000000000 0.25000000000000 0.87500000000000 88.906 > 6 7 Y 0.00000000000000 0.75000000000000 0.62500000000000 88.906 > 1 8 Y 0.25000000000000 0.87500000000000 0.00000000000000 88.906 > 2 9 Y 0.75000000000000 0.62500000000000 0.00000000000000 88.906 > 3 10 Y 0.87500000000000 0.00000000000000 0.25000000000000 88.906 > 4 11 Y 0.12500000000000 0.50000000000000 0.25000000000000 88.906 > 5 12 Y 0.50000000000000 0.75000000000000 0.37500000000000 88.906 > 6 *13 Sb 0.16133275078289 0.83866724921711 0.33866724921711 121.760 > 7 14 Sb 0.33866724921711 0.83866724921711 0.66133275078289 121.760 > 8 15 Sb 0.16133275078289 0.16133275078289 0.16133275078289 121.760 > 9 16 Sb 0.08866724921711 0.58866724921711 0.91133275078289 121.760 > 10 17 Sb 0.08866724921711 0.41133275078289 0.58866724921711 121.760 > 11 18 Sb 0.41133275078289 0.41133275078289 0.41133275078289 121.760 > 12 19 Sb 0.33866724921711 0.16133275078289 0.83866724921711 121.760 > 13 20 Sb 0.91133275078289 0.08866724921711 0.58866724921711 121.760 > 14 21 Sb 0.66133275078289 0.33866724921711 0.83866724921711 121.760 > 7 22 Sb 0.83866724921711 0.33866724921711 0.16133275078289 121.760 > 8 23 Sb 0.66133275078289 0.66133275078289 0.66133275078289 121.760 > 9 24 Sb 0.58866724921711 0.08866724921711 0.41133275078289 121.760 > 10 25 Sb 0.58866724921711 0.91133275078289 0.08866724921711 121.760 > 11 26 Sb 0.91133275078289 0.91133275078289 0.91133275078289 121.760 > 12 27 Sb 0.83866724921711 0.66133275078289 0.33866724921711 121.760 > 13 28 Sb 0.41133275078289 0.58866724921711 0.08866724921711 121.760 > 14 *29 Au 0.00000000000000 0.25000000000000 0.37500000000000 196.967 > 15 30 Au 0.87500000000000 0.50000000000000 0.75000000000000 196.967 > 16 31 Au 0.75000000000000 0.12500000000000 0.00000000000000 196.967 > 17 32 Au 0.25000000000000 0.37500000000000 0.00000000000000 196.967 > 18 33 Au 0.50000000000000 0.25000000000000 0.62500000000000 196.967 > 19 34 Au 0.12500000000000 0.00000000000000 0.75000000000000 196.967 > 20 35 Au 0.50000000000000 0.75000000000000 0.87500000000000 196.967 > 15 36 Au 0.37500000000000 0.00000000000000 0.25000000000000 196.967 > 16 37 Au 0.25000000000000 0.62500000000000 0.50000000000000 196.967 > 17 38 Au 0.75000000000000 0.87500000000000 0.50000000000000 196.967 > 18 39 Au 0.00000000000000 0.75000000000000 0.12500000000000 196.967 > 19 40 Au 0.62500000000000 0.50000000000000 0.25000000000000 196.967 > 20 -------------------------------- super cell -------------------------------- Lattice vectors: a 9.790548480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 9.790548480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 9.790548480000000 Atomic positions (fractional): *1 Y 0.50000000000000 0.25000000000000 0.12500000000000 88.906 > 1 2 Y 0.75000000000000 0.37500000000000 0.50000000000000 88.906 > 2 3 Y 0.25000000000000 0.12500000000000 0.50000000000000 88.906 > 3 4 Y 0.37500000000000 0.50000000000000 0.75000000000000 88.906 > 4 5 Y 0.62500000000000 0.00000000000000 0.75000000000000 88.906 > 5 6 Y 0.00000000000000 0.25000000000000 0.87500000000000 88.906 > 6 7 Y 0.00000000000000 0.75000000000000 0.62500000000000 88.906 > 1 8 Y 0.25000000000000 0.87500000000000 0.00000000000000 88.906 > 2 9 Y 0.75000000000000 0.62500000000000 0.00000000000000 88.906 > 3 10 Y 0.87500000000000 0.00000000000000 0.25000000000000 88.906 > 4 11 Y 0.12500000000000 0.50000000000000 0.25000000000000 88.906 > 5 12 Y 0.50000000000000 0.75000000000000 0.37500000000000 88.906 > 6 *13 Sb 0.16133275078289 0.83866724921711 0.33866724921711 121.760 > 7 14 Sb 0.33866724921711 0.83866724921711 0.66133275078289 121.760 > 8 15 Sb 0.16133275078289 0.16133275078289 0.16133275078289 121.760 > 9 16 Sb 0.08866724921711 0.58866724921711 0.91133275078289 121.760 > 10 17 Sb 0.08866724921711 0.41133275078289 0.58866724921711 121.760 > 11 18 Sb 0.41133275078289 0.41133275078289 0.41133275078289 121.760 > 12 19 Sb 0.33866724921711 0.16133275078289 0.83866724921711 121.760 > 13 20 Sb 0.91133275078289 0.08866724921711 0.58866724921711 121.760 > 14 21 Sb 0.66133275078289 0.33866724921711 0.83866724921711 121.760 > 7 22 Sb 0.83866724921711 0.33866724921711 0.16133275078289 121.760 > 8 23 Sb 0.66133275078289 0.66133275078289 0.66133275078289 121.760 > 9 24 Sb 0.58866724921711 0.08866724921711 0.41133275078289 121.760 > 10 25 Sb 0.58866724921711 0.91133275078289 0.08866724921711 121.760 > 11 26 Sb 0.91133275078289 0.91133275078289 0.91133275078289 121.760 > 12 27 Sb 0.83866724921711 0.66133275078289 0.33866724921711 121.760 > 13 28 Sb 0.41133275078289 0.58866724921711 0.08866724921711 121.760 > 14 *29 Au 0.00000000000000 0.25000000000000 0.37500000000000 196.967 > 15 30 Au 0.87500000000000 0.50000000000000 0.75000000000000 196.967 > 16 31 Au 0.75000000000000 0.12500000000000 0.00000000000000 196.967 > 17 32 Au 0.25000000000000 0.37500000000000 0.00000000000000 196.967 > 18 33 Au 0.50000000000000 0.25000000000000 0.62500000000000 196.967 > 19 34 Au 0.12500000000000 0.00000000000000 0.75000000000000 196.967 > 20 35 Au 0.50000000000000 0.75000000000000 0.87500000000000 196.967 > 15 36 Au 0.37500000000000 0.00000000000000 0.25000000000000 196.967 > 16 37 Au 0.25000000000000 0.62500000000000 0.50000000000000 196.967 > 17 38 Au 0.75000000000000 0.87500000000000 0.50000000000000 196.967 > 18 39 Au 0.00000000000000 0.75000000000000 0.12500000000000 196.967 > 19 40 Au 0.62500000000000 0.50000000000000 0.25000000000000 196.967 > 20 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 26.0865967 0.0000000 0.0000000 0.0000000 26.0865967 0.0000000 0.0000000 0.0000000 26.0865967 -------------------------- Born effective charges -------------------------- 1 Y 2.4344523 -0.0045948 0.0000000 0.0045948 2.4344523 0.0000000 0.0000000 0.0000000 4.2560064 2 Y 2.4344523 0.0000000 -0.0045948 0.0000000 4.2560064 0.0000000 0.0045948 0.0000000 2.4344523 3 Y 2.4344523 0.0000000 0.0045948 0.0000000 4.2560064 0.0000000 -0.0045948 0.0000000 2.4344523 4 Y 4.2560064 0.0000000 0.0000000 0.0000000 2.4344523 0.0045948 0.0000000 -0.0045948 2.4344523 5 Y 4.2560064 0.0000000 0.0000000 0.0000000 2.4344523 -0.0045948 0.0000000 0.0045948 2.4344523 6 Y 2.4344523 0.0045948 0.0000000 -0.0045948 2.4344523 0.0000000 0.0000000 0.0000000 4.2560064 7 Sb -1.4350617 0.0377863 0.5083919 0.5083919 -1.4350617 -0.0377863 0.0377863 -0.5083919 -1.4350617 8 Sb -1.4350617 -0.0377863 0.5083919 -0.5083919 -1.4350617 0.0377863 0.0377863 0.5083919 -1.4350617 9 Sb -1.4350617 -0.0377863 -0.5083919 -0.5083919 -1.4350617 -0.0377863 -0.0377863 -0.5083919 -1.4350617 10 Sb -1.4350617 -0.5083919 0.0377863 -0.0377863 -1.4350617 0.5083919 0.5083919 0.0377863 -1.4350617 11 Sb -1.4350617 0.5083919 -0.0377863 0.0377863 -1.4350617 0.5083919 -0.5083919 0.0377863 -1.4350617 12 Sb -1.4350617 -0.5083919 -0.0377863 -0.0377863 -1.4350617 -0.5083919 -0.5083919 -0.0377863 -1.4350617 13 Sb -1.4350617 0.0377863 -0.5083919 0.5083919 -1.4350617 0.0377863 -0.0377863 0.5083919 -1.4350617 14 Sb -1.4350617 0.5083919 0.0377863 0.0377863 -1.4350617 -0.5083919 0.5083919 -0.0377863 -1.4350617 15 Au -1.1699786 0.2249693 0.0000000 -0.2249693 -1.1699786 0.0000000 0.0000000 0.0000000 -1.0447069 16 Au -1.0447069 0.0000000 0.0000000 0.0000000 -1.1699786 0.2249693 0.0000000 -0.2249693 -1.1699786 17 Au -1.1699786 0.0000000 0.2249693 0.0000000 -1.0447069 0.0000000 -0.2249693 0.0000000 -1.1699786 18 Au -1.1699786 0.0000000 -0.2249693 0.0000000 -1.0447069 0.0000000 0.2249693 0.0000000 -1.1699786 19 Au -1.1699786 -0.2249693 0.0000000 0.2249693 -1.1699786 0.0000000 0.0000000 0.0000000 -1.0447069 20 Au -1.0447069 0.0000000 0.0000000 0.0000000 -1.1699786 -0.2249693 0.0000000 0.2249693 -1.1699786 ---------------------------------------------------------------------------- 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: 120/120 Permutation basis: 4182/4182 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 163 Number of blocks in projector: 163 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: 94 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 69 Use standard eigh solver. Tree of FC basis block matrices: - (163, 154), data: False |-- (69, 64), data: True |-- (94, 90), data: True ----- Solver_atoms: 1 -- 40 / 40 Time (Solver_compr_matrix_reshape): 0.001 Solver_block: 100 / 160 - Time: 0.041 Solver_block: 160 / 160 - Time: 0.016 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.058 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) Permutation basis: 120/120 Permutation basis: 4182/4182 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 163 Number of blocks in projector: 163 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: 94 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 69 Use standard eigh solver. Tree of FC basis block matrices: - (163, 154), data: False |-- (69, 64), data: True |-- (94, 90), data: True Max drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 00:18:07]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 00:18:07]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [1 1 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I-43d (220) ------------------------------ primitive cell ------------------------------ Lattice vectors: a -4.895274240000000 4.895274240000000 4.895274240000000 b 4.895274240000000 -4.895274240000000 4.895274240000000 c 4.895274240000000 4.895274240000000 -4.895274240000000 Atomic positions (fractional): 1 Y 0.37500000000000 0.62500000000000 0.75000000000000 88.906 2 Y 0.87500000000000 0.25000000000000 0.12500000000000 88.906 3 Y 0.62500000000000 0.75000000000000 0.37500000000000 88.906 4 Y 0.25000000000000 0.12500000000000 0.87500000000000 88.906 5 Y 0.75000000000000 0.37500000000000 0.62500000000000 88.906 6 Y 0.12500000000000 0.87500000000000 0.25000000000000 88.906 7 Sb 0.17733449843421 0.50000000000000 0.00000000000000 121.760 8 Sb 0.50000000000000 0.00000000000000 0.17733449843421 121.760 9 Sb 0.32266550156579 0.32266550156579 0.32266550156579 121.760 10 Sb 0.50000000000000 0.00000000000000 0.67733449843421 121.760 11 Sb 0.00000000000000 0.67733449843421 0.50000000000000 121.760 12 Sb 0.82266550156579 0.82266550156579 0.82266550156579 121.760 13 Sb 0.00000000000000 0.17733449843421 0.50000000000000 121.760 14 Sb 0.67733449843421 0.50000000000000 0.00000000000000 121.760 15 Au 0.62500000000000 0.37500000000000 0.25000000000000 196.967 16 Au 0.25000000000000 0.62500000000000 0.37500000000000 196.967 17 Au 0.12500000000000 0.75000000000000 0.87500000000000 196.967 18 Au 0.37500000000000 0.25000000000000 0.62500000000000 196.967 19 Au 0.87500000000000 0.12500000000000 0.75000000000000 196.967 20 Au 0.75000000000000 0.87500000000000 0.12500000000000 196.967 -------------------------------- supercell --------------------------------- Lattice vectors: a 9.790548480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 9.790548480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 9.790548480000000 Atomic positions (fractional): 1 Y 0.50000000000000 0.25000000000000 0.12500000000000 88.906 > 1 2 Y 0.75000000000000 0.37500000000000 0.50000000000000 88.906 > 2 3 Y 0.25000000000000 0.12500000000000 0.50000000000000 88.906 > 3 4 Y 0.37500000000000 0.50000000000000 0.75000000000000 88.906 > 4 5 Y 0.62500000000000 0.00000000000000 0.75000000000000 88.906 > 5 6 Y 0.00000000000000 0.25000000000000 0.87500000000000 88.906 > 6 7 Y 0.00000000000000 0.75000000000000 0.62500000000000 88.906 > 1 8 Y 0.25000000000000 0.87500000000000 0.00000000000000 88.906 > 2 9 Y 0.75000000000000 0.62500000000000 0.00000000000000 88.906 > 3 10 Y 0.87500000000000 0.00000000000000 0.25000000000000 88.906 > 4 11 Y 0.12500000000000 0.50000000000000 0.25000000000000 88.906 > 5 12 Y 0.50000000000000 0.75000000000000 0.37500000000000 88.906 > 6 13 Sb 0.16133275078289 0.83866724921711 0.33866724921711 121.760 > 13 14 Sb 0.33866724921711 0.83866724921711 0.66133275078289 121.760 > 14 15 Sb 0.16133275078289 0.16133275078289 0.16133275078289 121.760 > 15 16 Sb 0.08866724921711 0.58866724921711 0.91133275078289 121.760 > 16 17 Sb 0.08866724921711 0.41133275078289 0.58866724921711 121.760 > 17 18 Sb 0.41133275078289 0.41133275078289 0.41133275078289 121.760 > 18 19 Sb 0.33866724921711 0.16133275078289 0.83866724921711 121.760 > 19 20 Sb 0.91133275078289 0.08866724921711 0.58866724921711 121.760 > 20 21 Sb 0.66133275078289 0.33866724921711 0.83866724921711 121.760 > 13 22 Sb 0.83866724921711 0.33866724921711 0.16133275078289 121.760 > 14 23 Sb 0.66133275078289 0.66133275078289 0.66133275078289 121.760 > 15 24 Sb 0.58866724921711 0.08866724921711 0.41133275078289 121.760 > 16 25 Sb 0.58866724921711 0.91133275078289 0.08866724921711 121.760 > 17 26 Sb 0.91133275078289 0.91133275078289 0.91133275078289 121.760 > 18 27 Sb 0.83866724921711 0.66133275078289 0.33866724921711 121.760 > 19 28 Sb 0.41133275078289 0.58866724921711 0.08866724921711 121.760 > 20 29 Au 0.00000000000000 0.25000000000000 0.37500000000000 196.967 > 29 30 Au 0.87500000000000 0.50000000000000 0.75000000000000 196.967 > 30 31 Au 0.75000000000000 0.12500000000000 0.00000000000000 196.967 > 31 32 Au 0.25000000000000 0.37500000000000 0.00000000000000 196.967 > 32 33 Au 0.50000000000000 0.25000000000000 0.62500000000000 196.967 > 33 34 Au 0.12500000000000 0.00000000000000 0.75000000000000 196.967 > 34 35 Au 0.50000000000000 0.75000000000000 0.87500000000000 196.967 > 29 36 Au 0.37500000000000 0.00000000000000 0.25000000000000 196.967 > 30 37 Au 0.25000000000000 0.62500000000000 0.50000000000000 196.967 > 31 38 Au 0.75000000000000 0.87500000000000 0.50000000000000 196.967 > 32 39 Au 0.00000000000000 0.75000000000000 0.12500000000000 196.967 > 33 40 Au 0.62500000000000 0.50000000000000 0.25000000000000 196.967 > 34 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 26.0865967 0.0000000 0.0000000 0.0000000 26.0865967 0.0000000 0.0000000 0.0000000 26.0865967 -------------------------- Born effective charges -------------------------- 1 Y 2.4344523 -0.0045948 0.0000000 0.0045948 2.4344523 0.0000000 0.0000000 0.0000000 4.2560064 2 Y 2.4344523 0.0000000 -0.0045948 0.0000000 4.2560064 0.0000000 0.0045948 0.0000000 2.4344523 3 Y 2.4344523 0.0000000 0.0045948 0.0000000 4.2560064 0.0000000 -0.0045948 0.0000000 2.4344523 4 Y 4.2560064 0.0000000 0.0000000 0.0000000 2.4344523 0.0045948 0.0000000 -0.0045948 2.4344523 5 Y 4.2560064 0.0000000 0.0000000 0.0000000 2.4344523 -0.0045948 0.0000000 0.0045948 2.4344523 6 Y 2.4344523 0.0045948 0.0000000 -0.0045948 2.4344523 0.0000000 0.0000000 0.0000000 4.2560064 7 Sb -1.4350617 0.0377863 0.5083919 0.5083919 -1.4350617 -0.0377863 0.0377863 -0.5083919 -1.4350617 8 Sb -1.4350617 -0.0377863 0.5083919 -0.5083919 -1.4350617 0.0377863 0.0377863 0.5083919 -1.4350617 9 Sb -1.4350617 -0.0377863 -0.5083919 -0.5083919 -1.4350617 -0.0377863 -0.0377863 -0.5083919 -1.4350617 10 Sb -1.4350617 -0.5083919 0.0377863 -0.0377863 -1.4350617 0.5083919 0.5083919 0.0377863 -1.4350617 11 Sb -1.4350617 0.5083919 -0.0377863 0.0377863 -1.4350617 0.5083919 -0.5083919 0.0377863 -1.4350617 12 Sb -1.4350617 -0.5083919 -0.0377863 -0.0377863 -1.4350617 -0.5083919 -0.5083919 -0.0377863 -1.4350617 13 Sb -1.4350617 0.0377863 -0.5083919 0.5083919 -1.4350617 0.0377863 -0.0377863 0.5083919 -1.4350617 14 Sb -1.4350617 0.5083919 0.0377863 0.0377863 -1.4350617 -0.5083919 0.5083919 -0.0377863 -1.4350617 15 Au -1.1699786 0.2249693 0.0000000 -0.2249693 -1.1699786 0.0000000 0.0000000 0.0000000 -1.0447069 16 Au -1.0447069 0.0000000 0.0000000 0.0000000 -1.1699786 0.2249693 0.0000000 -0.2249693 -1.1699786 17 Au -1.1699786 0.0000000 0.2249693 0.0000000 -1.0447069 0.0000000 -0.2249693 0.0000000 -1.1699786 18 Au -1.1699786 0.0000000 -0.2249693 0.0000000 -1.0447069 0.0000000 0.2249693 0.0000000 -1.1699786 19 Au -1.1699786 -0.2249693 0.0000000 0.2249693 -1.1699786 0.0000000 0.0000000 0.0000000 -1.0447069 20 Au -1.0447069 0.0000000 0.0000000 0.0000000 -1.1699786 -0.2249693 0.0000000 0.2249693 -1.1699786 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 13, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 29, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000000 (xyz) 0.00000000 (xyz) 0.00000000 (xzy) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (zx) 0.00000000 (zx) 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-08 00:18:10]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 00:18:10]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [1 1 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I-43d (220) ------------------------------ primitive cell ------------------------------ Lattice vectors: a -4.895274240000000 4.895274240000000 4.895274240000000 b 4.895274240000000 -4.895274240000000 4.895274240000000 c 4.895274240000000 4.895274240000000 -4.895274240000000 Atomic positions (fractional): 1 Y 0.37500000000000 0.62500000000000 0.75000000000000 88.906 2 Y 0.87500000000000 0.25000000000000 0.12500000000000 88.906 3 Y 0.62500000000000 0.75000000000000 0.37500000000000 88.906 4 Y 0.25000000000000 0.12500000000000 0.87500000000000 88.906 5 Y 0.75000000000000 0.37500000000000 0.62500000000000 88.906 6 Y 0.12500000000000 0.87500000000000 0.25000000000000 88.906 7 Sb 0.17733449843421 0.50000000000000 0.00000000000000 121.760 8 Sb 0.50000000000000 0.00000000000000 0.17733449843421 121.760 9 Sb 0.32266550156579 0.32266550156579 0.32266550156579 121.760 10 Sb 0.50000000000000 0.00000000000000 0.67733449843421 121.760 11 Sb 0.00000000000000 0.67733449843421 0.50000000000000 121.760 12 Sb 0.82266550156579 0.82266550156579 0.82266550156579 121.760 13 Sb 0.00000000000000 0.17733449843421 0.50000000000000 121.760 14 Sb 0.67733449843421 0.50000000000000 0.00000000000000 121.760 15 Au 0.62500000000000 0.37500000000000 0.25000000000000 196.967 16 Au 0.25000000000000 0.62500000000000 0.37500000000000 196.967 17 Au 0.12500000000000 0.75000000000000 0.87500000000000 196.967 18 Au 0.37500000000000 0.25000000000000 0.62500000000000 196.967 19 Au 0.87500000000000 0.12500000000000 0.75000000000000 196.967 20 Au 0.75000000000000 0.87500000000000 0.12500000000000 196.967 -------------------------------- supercell --------------------------------- Lattice vectors: a 9.790548480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 9.790548480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 9.790548480000000 Atomic positions (fractional): 1 Y 0.50000000000000 0.25000000000000 0.12500000000000 88.906 > 1 2 Y 0.75000000000000 0.37500000000000 0.50000000000000 88.906 > 2 3 Y 0.25000000000000 0.12500000000000 0.50000000000000 88.906 > 3 4 Y 0.37500000000000 0.50000000000000 0.75000000000000 88.906 > 4 5 Y 0.62500000000000 0.00000000000000 0.75000000000000 88.906 > 5 6 Y 0.00000000000000 0.25000000000000 0.87500000000000 88.906 > 6 7 Y 0.00000000000000 0.75000000000000 0.62500000000000 88.906 > 1 8 Y 0.25000000000000 0.87500000000000 0.00000000000000 88.906 > 2 9 Y 0.75000000000000 0.62500000000000 0.00000000000000 88.906 > 3 10 Y 0.87500000000000 0.00000000000000 0.25000000000000 88.906 > 4 11 Y 0.12500000000000 0.50000000000000 0.25000000000000 88.906 > 5 12 Y 0.50000000000000 0.75000000000000 0.37500000000000 88.906 > 6 13 Sb 0.16133275078289 0.83866724921711 0.33866724921711 121.760 > 13 14 Sb 0.33866724921711 0.83866724921711 0.66133275078289 121.760 > 14 15 Sb 0.16133275078289 0.16133275078289 0.16133275078289 121.760 > 15 16 Sb 0.08866724921711 0.58866724921711 0.91133275078289 121.760 > 16 17 Sb 0.08866724921711 0.41133275078289 0.58866724921711 121.760 > 17 18 Sb 0.41133275078289 0.41133275078289 0.41133275078289 121.760 > 18 19 Sb 0.33866724921711 0.16133275078289 0.83866724921711 121.760 > 19 20 Sb 0.91133275078289 0.08866724921711 0.58866724921711 121.760 > 20 21 Sb 0.66133275078289 0.33866724921711 0.83866724921711 121.760 > 13 22 Sb 0.83866724921711 0.33866724921711 0.16133275078289 121.760 > 14 23 Sb 0.66133275078289 0.66133275078289 0.66133275078289 121.760 > 15 24 Sb 0.58866724921711 0.08866724921711 0.41133275078289 121.760 > 16 25 Sb 0.58866724921711 0.91133275078289 0.08866724921711 121.760 > 17 26 Sb 0.91133275078289 0.91133275078289 0.91133275078289 121.760 > 18 27 Sb 0.83866724921711 0.66133275078289 0.33866724921711 121.760 > 19 28 Sb 0.41133275078289 0.58866724921711 0.08866724921711 121.760 > 20 29 Au 0.00000000000000 0.25000000000000 0.37500000000000 196.967 > 29 30 Au 0.87500000000000 0.50000000000000 0.75000000000000 196.967 > 30 31 Au 0.75000000000000 0.12500000000000 0.00000000000000 196.967 > 31 32 Au 0.25000000000000 0.37500000000000 0.00000000000000 196.967 > 32 33 Au 0.50000000000000 0.25000000000000 0.62500000000000 196.967 > 33 34 Au 0.12500000000000 0.00000000000000 0.75000000000000 196.967 > 34 35 Au 0.50000000000000 0.75000000000000 0.87500000000000 196.967 > 29 36 Au 0.37500000000000 0.00000000000000 0.25000000000000 196.967 > 30 37 Au 0.25000000000000 0.62500000000000 0.50000000000000 196.967 > 31 38 Au 0.75000000000000 0.87500000000000 0.50000000000000 196.967 > 32 39 Au 0.00000000000000 0.75000000000000 0.12500000000000 196.967 > 33 40 Au 0.62500000000000 0.50000000000000 0.25000000000000 196.967 > 34 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 26.0865967 0.0000000 0.0000000 0.0000000 26.0865967 0.0000000 0.0000000 0.0000000 26.0865967 -------------------------- Born effective charges -------------------------- 1 Y 2.4344523 -0.0045948 0.0000000 0.0045948 2.4344523 0.0000000 0.0000000 0.0000000 4.2560064 2 Y 2.4344523 0.0000000 -0.0045948 0.0000000 4.2560064 0.0000000 0.0045948 0.0000000 2.4344523 3 Y 2.4344523 0.0000000 0.0045948 0.0000000 4.2560064 0.0000000 -0.0045948 0.0000000 2.4344523 4 Y 4.2560064 0.0000000 0.0000000 0.0000000 2.4344523 0.0045948 0.0000000 -0.0045948 2.4344523 5 Y 4.2560064 0.0000000 0.0000000 0.0000000 2.4344523 -0.0045948 0.0000000 0.0045948 2.4344523 6 Y 2.4344523 0.0045948 0.0000000 -0.0045948 2.4344523 0.0000000 0.0000000 0.0000000 4.2560064 7 Sb -1.4350617 0.0377863 0.5083919 0.5083919 -1.4350617 -0.0377863 0.0377863 -0.5083919 -1.4350617 8 Sb -1.4350617 -0.0377863 0.5083919 -0.5083919 -1.4350617 0.0377863 0.0377863 0.5083919 -1.4350617 9 Sb -1.4350617 -0.0377863 -0.5083919 -0.5083919 -1.4350617 -0.0377863 -0.0377863 -0.5083919 -1.4350617 10 Sb -1.4350617 -0.5083919 0.0377863 -0.0377863 -1.4350617 0.5083919 0.5083919 0.0377863 -1.4350617 11 Sb -1.4350617 0.5083919 -0.0377863 0.0377863 -1.4350617 0.5083919 -0.5083919 0.0377863 -1.4350617 12 Sb -1.4350617 -0.5083919 -0.0377863 -0.0377863 -1.4350617 -0.5083919 -0.5083919 -0.0377863 -1.4350617 13 Sb -1.4350617 0.0377863 -0.5083919 0.5083919 -1.4350617 0.0377863 -0.0377863 0.5083919 -1.4350617 14 Sb -1.4350617 0.5083919 0.0377863 0.0377863 -1.4350617 -0.5083919 0.5083919 -0.0377863 -1.4350617 15 Au -1.1699786 0.2249693 0.0000000 -0.2249693 -1.1699786 0.0000000 0.0000000 0.0000000 -1.0447069 16 Au -1.0447069 0.0000000 0.0000000 0.0000000 -1.1699786 0.2249693 0.0000000 -0.2249693 -1.1699786 17 Au -1.1699786 0.0000000 0.2249693 0.0000000 -1.0447069 0.0000000 -0.2249693 0.0000000 -1.1699786 18 Au -1.1699786 0.0000000 -0.2249693 0.0000000 -1.0447069 0.0000000 0.2249693 0.0000000 -1.1699786 19 Au -1.1699786 -0.2249693 0.0000000 0.2249693 -1.1699786 0.0000000 0.0000000 0.0000000 -1.0447069 20 Au -1.0447069 0.0000000 0.0000000 0.0000000 -1.1699786 -0.2249693 0.0000000 0.2249693 -1.1699786 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000000 (xyz) 0.00000000 (xyz) 0.00000000 (xzy) Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 7 7 7 ] 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.53, Number of G-points: 321, Lambda: 0.28 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/20) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 20 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 1.621 ( -0.000 -0.000 0.000) 0.000 1.621 ( 0.000 0.000 0.000) 0.000 1.621 ( -0.000 0.000 0.000) 0.000 1.649 ( 0.000 0.000 0.000) 0.000 1.649 ( 0.000 -0.000 0.000) 0.000 1.649 ( -0.000 -0.000 0.000) 0.000 1.873 ( -0.000 -0.000 0.000) 0.000 1.873 ( 0.000 0.000 0.000) 0.000 1.873 ( 0.000 0.000 0.000) 0.000 1.991 ( 0.000 0.000 0.000) 0.000 2.148 ( 0.000 0.000 0.000) 0.000 2.148 ( -0.000 -0.000 0.000) 0.000 2.187 ( 0.000 -0.000 0.000) 0.000 2.187 ( 0.000 0.000 0.000) 0.000 2.187 ( 0.000 -0.000 0.000) 0.000 2.929 ( -0.000 -0.000 0.000) 0.000 3.026 ( -0.000 0.000 0.000) 0.000 3.026 ( 0.000 0.000 0.000) 0.000 3.054 ( 0.000 0.000 0.000) 0.000 3.054 ( -0.000 -0.000 0.000) 0.000 3.054 ( -0.000 0.000 0.000) 0.000 3.363 ( 0.000 0.000 0.000) 0.000 3.363 ( 0.000 0.000 0.000) 0.000 3.363 ( 0.000 0.000 0.000) 0.000 3.648 ( 0.000 -0.000 0.000) 0.000 3.648 ( 0.000 -0.000 0.000) 0.000 3.648 ( 0.000 0.000 0.000) 0.000 3.726 ( 0.000 -0.000 0.000) 0.000 3.726 ( 0.000 -0.000 0.000) 0.000 3.726 ( 0.000 0.000 0.000) 0.000 4.096 ( 0.000 -0.000 0.000) 0.000 4.096 ( -0.000 0.000 0.000) 0.000 4.096 ( 0.000 0.000 0.000) 0.000 4.137 ( 0.000 -0.000 0.000) 0.000 4.153 ( -0.000 0.000 0.000) 0.000 4.153 ( 0.000 0.000 0.000) 0.000 4.153 ( 0.000 -0.000 0.000) 0.000 4.189 ( -0.000 -0.000 0.000) 0.000 4.227 ( 0.000 0.000 0.000) 0.000 4.227 ( 0.000 0.000 0.000) 0.000 4.227 ( 0.000 -0.000 0.000) 0.000 4.577 ( -0.000 0.000 0.000) 0.000 4.577 ( -0.000 -0.000 0.000) 0.000 4.592 ( 0.000 0.000 0.000) 0.000 4.592 ( 0.000 0.000 0.000) 0.000 4.592 ( -0.000 -0.000 0.000) 0.000 4.715 ( 0.000 0.000 0.000) 0.000 4.715 ( -0.000 0.000 0.000) 0.000 4.715 ( -0.000 -0.000 0.000) 0.000 4.858 ( -0.000 -0.000 0.000) 0.000 4.858 ( -0.000 0.000 0.000) 0.000 4.858 ( -0.000 -0.000 0.000) 0.000 5.219 ( 0.000 -0.000 0.000) 0.000 5.219 ( 0.000 0.000 0.000) 0.000 5.219 ( 0.000 0.000 0.000) 0.000 5.348 ( -0.000 -0.000 0.000) 0.000 5.348 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/20) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 64 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.416 ( -0.000 13.937 13.937) 19.710 0.422 ( -0.000 14.101 14.101) 19.942 0.784 ( -0.000 25.977 25.977) 36.738 1.616 ( 0.000 -0.427 -0.427) 0.604 1.617 ( 0.000 -1.236 -1.236) 1.748 1.639 ( 0.000 -0.549 -0.549) 0.776 1.650 ( -0.000 0.201 0.201) 0.285 1.660 ( -0.000 1.372 1.372) 1.940 1.676 ( -0.000 1.639 1.639) 2.318 1.895 ( -0.000 0.436 0.436) 0.617 1.904 ( -0.000 1.819 1.819) 2.572 1.904 ( -0.000 1.770 1.770) 2.503 2.002 ( -0.000 0.698 0.698) 0.987 2.119 ( 0.000 -1.626 -1.626) 2.300 2.142 ( 0.000 -0.235 -0.235) 0.333 2.169 ( 0.000 -1.158 -1.158) 1.638 2.188 ( 0.000 -0.068 -0.068) 0.096 2.192 ( -0.000 0.239 0.239) 0.338 2.914 ( 0.000 -1.105 -1.105) 1.562 2.945 ( 0.000 -2.602 -2.602) 3.680 2.980 ( 0.000 -1.687 -1.687) 2.386 3.050 ( 0.000 -0.344 -0.344) 0.487 3.104 ( -0.000 1.956 1.956) 2.766 3.151 ( -0.000 3.584 3.584) 5.068 3.340 ( 0.000 -1.552 -1.552) 2.195 3.357 ( 0.000 -0.485 -0.485) 0.686 3.362 ( 0.000 -0.204 -0.204) 0.288 3.553 ( 0.000 -3.998 -3.998) 5.654 3.558 ( 0.000 -3.403 -3.403) 4.813 3.654 ( 0.000 -4.204 -4.204) 5.945 3.689 ( -0.000 0.437 0.437) 0.618 3.790 ( -0.000 1.949 1.949) 2.757 3.793 ( -0.000 2.295 2.295) 3.246 4.022 ( 0.000 -3.348 -3.348) 4.735 4.026 ( 0.000 -2.727 -2.727) 3.857 4.066 ( 0.000 -1.685 -1.685) 2.382 4.080 ( 0.000 -2.203 -2.203) 3.116 4.162 ( 0.000 -1.566 -1.566) 2.215 4.197 ( 0.000 -0.304 -0.304) 0.431 4.209 ( -0.000 1.404 1.404) 1.985 4.245 ( -0.000 3.335 3.335) 4.716 4.250 ( -0.000 3.275 3.275) 4.631 4.252 ( -0.000 1.780 1.780) 2.517 4.445 ( 0.000 -3.010 -3.010) 4.256 4.527 ( 0.000 -1.925 -1.925) 2.722 4.548 ( 0.000 -1.298 -1.298) 1.836 4.588 ( 0.000 -1.772 -1.772) 2.506 4.602 ( 0.000 -0.151 -0.151) 0.213 4.647 ( -0.000 1.846 1.846) 2.611 4.695 ( 0.000 -0.708 -0.708) 1.001 4.704 ( 0.000 -0.332 -0.332) 0.469 4.717 ( -0.000 0.089 0.089) 0.126 4.847 ( 0.000 -0.531 -0.531) 0.751 4.866 ( -0.000 0.467 0.467) 0.660 4.879 ( -0.000 1.488 1.488) 2.105 5.228 ( 0.000 -0.318 -0.318) 0.450 5.234 ( -0.000 0.749 0.749) 1.059 5.242 ( -0.000 1.316 1.316) 1.862 5.341 ( 0.000 -0.266 -0.266) 0.376 5.360 ( -0.000 0.723 0.723) 1.023 ======================= Grid point 2 (3/20) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 64 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.804 ( -0.000 12.545 12.545) 17.741 0.813 ( -0.000 12.529 12.529) 17.718 1.462 ( -0.000 15.906 15.906) 22.494 1.546 ( 0.000 -3.574 -3.574) 5.054 1.584 ( -0.000 0.802 0.802) 1.134 1.635 ( -0.000 0.913 0.913) 1.292 1.654 ( 0.000 -0.088 -0.088) 0.125 1.717 ( -0.000 1.995 1.995) 2.821 1.727 ( -0.000 1.764 1.764) 2.495 1.908 ( -0.000 0.649 0.649) 0.918 1.939 ( -0.000 0.088 0.088) 0.124 1.979 ( -0.000 2.809 2.809) 3.972 2.030 ( -0.000 0.692 0.692) 0.978 2.080 ( 0.000 -0.394 -0.394) 0.558 2.121 ( 0.000 -0.806 -0.806) 1.139 2.157 ( -0.000 0.495 0.495) 0.701 2.183 ( 0.000 -0.215 -0.215) 0.304 2.206 ( -0.000 0.972 0.972) 1.375 2.862 ( 0.000 -2.058 -2.058) 2.911 2.897 ( 0.000 -0.762 -0.762) 1.078 2.938 ( 0.000 -1.367 -1.367) 1.933 3.027 ( 0.000 -1.286 -1.286) 1.819 3.149 ( -0.000 1.036 1.036) 1.465 3.220 ( -0.000 0.412 0.412) 0.583 3.281 ( 0.000 -2.279 -2.279) 3.224 3.331 ( 0.000 -2.458 -2.458) 3.476 3.338 ( 0.000 -0.479 -0.479) 0.677 3.474 ( 0.000 -1.088 -1.088) 1.539 3.518 ( -0.000 1.288 1.288) 1.821 3.527 ( 0.000 -3.474 -3.474) 4.913 3.701 ( -0.000 0.226 0.226) 0.319 3.802 ( 0.000 -1.170 -1.170) 1.655 3.826 ( 0.000 -0.135 -0.135) 0.191 3.909 ( 0.000 -4.453 -4.453) 6.298 3.959 ( 0.000 -1.809 -1.809) 2.558 4.012 ( 0.000 -2.511 -2.511) 3.551 4.026 ( 0.000 -0.598 -0.598) 0.846 4.107 ( 0.000 -1.527 -1.527) 2.159 4.162 ( 0.000 -1.619 -1.619) 2.290 4.216 ( 0.000 -0.837 -0.837) 1.184 4.296 ( -0.000 0.672 0.672) 0.950 4.324 ( -0.000 1.011 1.011) 1.430 4.346 ( -0.000 3.168 3.168) 4.481 4.398 ( -0.000 0.346 0.346) 0.489 4.499 ( -0.000 0.349 0.349) 0.493 4.515 ( 0.000 -0.797 -0.797) 1.126 4.517 ( 0.000 -2.147 -2.147) 3.037 4.565 ( 0.000 -2.117 -2.117) 2.994 4.680 ( -0.000 0.805 0.805) 1.138 4.681 ( 0.000 -0.401 -0.401) 0.568 4.712 ( -0.000 0.357 0.357) 0.505 4.717 ( 0.000 -0.128 -0.128) 0.181 4.836 ( 0.000 -0.151 -0.151) 0.214 4.875 ( -0.000 0.012 0.012) 0.017 4.942 ( -0.000 2.732 2.732) 3.864 5.195 ( 0.000 -1.948 -1.948) 2.755 5.263 ( -0.000 1.022 1.022) 1.445 5.279 ( -0.000 1.033 1.033) 1.460 5.345 ( -0.000 0.601 0.601) 0.850 5.380 ( -0.000 0.524 0.524) 0.741 ======================= Grid point 3 (4/20) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 64 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.136 ( -0.000 9.813 9.813) 13.877 1.143 ( -0.000 9.718 9.718) 13.743 1.402 ( 0.000 -6.326 -6.326) 8.947 1.410 ( 0.000 -6.463 -6.463) 9.140 1.631 ( -0.000 0.250 0.250) 0.354 1.642 ( 0.000 -0.561 -0.561) 0.794 1.749 ( -0.000 0.863 0.863) 1.221 1.756 ( -0.000 0.511 0.511) 0.723 1.899 ( -0.000 2.321 2.321) 3.283 1.926 ( 0.000 -0.597 -0.597) 0.844 1.964 ( -0.000 3.276 3.276) 4.633 2.012 ( 0.000 -1.223 -1.223) 1.730 2.074 ( -0.000 4.447 4.447) 6.288 2.124 ( -0.000 0.170 0.170) 0.240 2.126 ( -0.000 3.148 3.148) 4.452 2.163 ( 0.000 -0.288 -0.288) 0.407 2.187 ( -0.000 1.237 1.237) 1.749 2.240 ( -0.000 0.307 0.307) 0.434 2.825 ( -0.000 0.229 0.229) 0.324 2.880 ( 0.000 -2.623 -2.623) 3.710 2.909 ( -0.000 1.791 1.791) 2.533 2.980 ( 0.000 -2.125 -2.125) 3.006 3.168 ( -0.000 0.500 0.500) 0.708 3.189 ( 0.000 -1.347 -1.347) 1.905 3.221 ( 0.000 -1.529 -1.529) 2.163 3.232 ( 0.000 -3.016 -3.016) 4.265 3.348 ( -0.000 1.134 1.134) 1.604 3.404 ( 0.000 -2.705 -2.705) 3.826 3.548 ( -0.000 3.433 3.433) 4.855 3.613 ( -0.000 4.009 4.009) 5.670 3.667 ( 0.000 -2.724 -2.724) 3.853 3.763 ( 0.000 -4.874 -4.874) 6.892 3.775 ( -0.000 0.514 0.514) 0.728 3.807 ( 0.000 -1.248 -1.248) 1.765 3.920 ( 0.000 -0.799 -0.799) 1.129 3.935 ( 0.000 -2.079 -2.079) 2.940 4.040 ( -0.000 1.372 1.372) 1.940 4.086 ( 0.000 -1.122 -1.122) 1.587 4.128 ( -0.000 0.327 0.327) 0.462 4.170 ( 0.000 -1.914 -1.914) 2.707 4.308 ( -0.000 0.422 0.422) 0.597 4.322 ( 0.000 -0.493 -0.493) 0.697 4.419 ( -0.000 1.699 1.699) 2.402 4.432 ( -0.000 1.028 1.028) 1.454 4.480 ( 0.000 -0.417 -0.417) 0.589 4.491 ( 0.000 -1.215 -1.215) 1.718 4.517 ( 0.000 -0.415 -0.415) 0.587 4.517 ( -0.000 0.153 0.153) 0.216 4.680 ( -0.000 0.526 0.526) 0.744 4.697 ( -0.000 0.060 0.060) 0.084 4.709 ( 0.000 -0.230 -0.230) 0.325 4.713 ( 0.000 -0.147 -0.147) 0.208 4.842 ( -0.000 0.580 0.580) 0.820 4.864 ( 0.000 -0.725 -0.725) 1.025 5.030 ( -0.000 3.153 3.153) 4.458 5.121 ( 0.000 -2.963 -2.963) 4.190 5.289 ( -0.000 0.655 0.655) 0.927 5.297 ( -0.000 0.150 0.150) 0.212 5.369 ( -0.000 0.818 0.818) 1.156 5.385 ( 0.000 -0.230 -0.230) 0.326 ======================= Grid point 8 (5/20) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.705 ( 7.734 7.734 14.854) 18.446 0.710 ( 7.659 7.659 14.870) 18.397 1.308 ( 12.398 12.398 26.063) 31.412 1.587 ( -1.142 -1.142 -2.187) 2.719 1.599 ( -1.232 -1.232 -1.206) 2.119 1.628 ( -0.567 -0.567 -0.289) 0.853 1.629 ( -1.264 -1.264 0.297) 1.811 1.708 ( 1.612 1.612 1.176) 2.566 1.716 ( 1.466 1.466 2.459) 3.216 1.910 ( 0.634 0.634 1.121) 1.436 1.939 ( 0.935 0.935 1.151) 1.753 1.948 ( 0.904 0.904 2.372) 2.694 2.011 ( 0.351 0.351 -0.065) 0.501 2.098 ( -0.206 -0.206 -0.467) 0.550 2.131 ( -0.973 -0.973 0.231) 1.396 2.158 ( -0.247 -0.247 -0.471) 0.586 2.167 ( 0.275 0.275 -0.712) 0.811 2.208 ( -0.781 -0.781 1.944) 2.236 2.886 ( -0.902 -0.902 -1.917) 2.303 2.903 ( -0.166 -0.166 -2.138) 2.151 2.945 ( -2.007 -2.007 0.009) 2.838 3.047 ( 0.346 0.346 -1.586) 1.660 3.123 ( 2.759 2.759 -2.013) 4.390 3.216 ( -1.153 -1.153 4.402) 4.694 3.313 ( -0.436 -0.436 -2.711) 2.780 3.315 ( -0.336 -0.336 -2.997) 3.034 3.360 ( -1.981 -1.981 0.552) 2.856 3.498 ( -0.343 -0.343 -1.662) 1.731 3.498 ( -0.532 -0.532 -1.207) 1.423 3.556 ( -3.031 -3.031 -4.653) 6.326 3.706 ( 1.354 1.354 -0.041) 1.915 3.812 ( 0.889 0.889 -0.865) 1.527 3.824 ( 1.085 1.085 0.069) 1.536 3.950 ( -1.209 -1.209 -3.904) 4.262 3.958 ( -1.606 -1.606 -3.825) 4.449 4.022 ( -1.838 -1.838 -3.324) 4.220 4.044 ( -0.930 -0.930 -1.933) 2.338 4.120 ( -1.225 -1.225 -2.634) 3.152 4.188 ( -0.337 -0.337 0.254) 0.540 4.194 ( -1.665 -1.665 -0.064) 2.355 4.297 ( 1.269 1.269 0.934) 2.023 4.313 ( 2.488 2.488 1.375) 3.778 4.338 ( 2.063 2.063 4.265) 5.167 4.393 ( -1.894 -1.894 0.185) 2.685 4.481 ( -0.532 -0.532 -2.341) 2.460 4.537 ( -2.873 -2.873 2.458) 4.749 4.539 ( -1.887 -1.887 -2.718) 3.809 4.563 ( -2.036 -2.036 -0.833) 2.997 4.680 ( -0.103 -0.103 -1.462) 1.469 4.686 ( 1.786 1.786 -0.575) 2.591 4.698 ( 0.825 0.825 -0.475) 1.259 4.717 ( -1.130 -1.130 0.839) 1.805 4.855 ( 1.737 1.737 -1.302) 2.780 4.862 ( 1.958 1.958 -0.900) 2.911 4.928 ( 0.291 0.291 4.022) 4.043 5.224 ( -0.679 -0.679 0.124) 0.968 5.241 ( 0.026 0.026 -0.727) 0.728 5.267 ( 0.208 0.208 1.836) 1.859 5.338 ( 1.570 1.570 -0.942) 2.412 5.380 ( -0.437 -0.437 1.960) 2.055 ======================= Grid point 9 (6/20) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.028 ( 4.295 8.351 12.513) 15.645 1.033 ( 3.981 8.440 12.352) 15.480 1.499 ( 0.983 -4.961 -3.791) 6.321 1.506 ( 0.826 -5.010 -3.353) 6.084 1.603 ( -2.742 0.640 -0.558) 2.871 1.624 ( -3.020 0.018 -0.303) 3.035 1.736 ( 0.868 0.949 0.750) 1.489 1.751 ( 0.373 0.788 1.663) 1.878 1.819 ( 0.935 5.753 11.037) 12.481 1.938 ( 0.975 -0.007 -1.164) 1.518 1.951 ( 2.072 1.377 1.965) 3.171 2.001 ( -0.993 0.596 -2.190) 2.478 2.017 ( 1.297 4.479 4.487) 6.471 2.109 ( -0.447 -1.011 1.888) 2.188 2.117 ( -0.973 -0.697 2.196) 2.502 2.176 ( 0.866 0.318 1.732) 1.962 2.187 ( 0.027 -0.001 2.350) 2.351 2.232 ( -0.297 0.130 2.230) 2.254 2.843 ( 0.808 -0.866 -1.014) 1.559 2.893 ( 0.688 1.079 -0.213) 1.298 2.922 ( 0.398 -1.320 -0.471) 1.457 3.009 ( 1.157 -0.692 -2.651) 2.974 3.126 ( -0.925 2.394 -2.269) 3.425 3.211 ( -2.543 -3.498 0.024) 4.325 3.245 ( 0.226 -1.548 -2.955) 3.343 3.263 ( -0.856 -0.922 -4.332) 4.511 3.333 ( -1.951 -1.521 0.917) 2.638 3.436 ( -1.625 -3.249 -3.210) 4.848 3.528 ( 2.103 2.104 2.865) 4.130 3.558 ( 0.971 3.782 3.810) 5.455 3.715 ( 2.372 0.170 -1.332) 2.726 3.782 ( 0.955 -0.391 -2.393) 2.606 3.811 ( 0.347 -0.994 -2.604) 2.808 3.849 ( 1.174 -3.074 -4.049) 5.217 3.903 ( -2.022 0.651 -3.294) 3.920 3.941 ( -1.931 -1.110 -4.391) 4.924 4.017 ( -1.583 -0.117 0.206) 1.601 4.071 ( -1.609 -0.485 -1.644) 2.351 4.147 ( -1.296 -2.269 -1.341) 2.937 4.190 ( -0.152 -1.147 0.150) 1.167 4.297 ( 0.943 1.806 -1.732) 2.674 4.337 ( 2.173 0.968 -1.693) 2.920 4.403 ( 0.826 0.226 2.739) 2.870 4.419 ( 0.281 0.711 2.039) 2.177 4.463 ( -2.288 1.110 -1.303) 2.857 4.500 ( -0.165 -1.129 0.691) 1.334 4.510 ( -1.438 -1.357 -0.317) 2.002 4.555 ( -0.803 -2.873 3.259) 4.418 4.656 ( -2.302 -0.767 -0.245) 2.439 4.683 ( -0.770 0.397 -0.795) 1.176 4.701 ( 0.304 0.271 -0.802) 0.900 4.711 ( -1.802 -1.250 0.432) 2.235 4.861 ( 3.692 1.880 -1.099) 4.286 4.887 ( 2.774 2.006 -0.711) 3.496 5.000 ( 0.938 1.706 3.870) 4.332 5.175 ( 1.414 -3.146 -2.188) 4.085 5.247 ( -2.832 1.833 -0.772) 3.460 5.292 ( -0.558 -0.530 1.756) 1.917 5.360 ( 1.301 1.705 0.249) 2.159 5.395 ( 0.011 -0.613 0.894) 1.084 ======================= Grid point 10 (7/20) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.259 ( 1.442 -1.143 1.143) 2.166 1.261 ( 1.453 -1.228 1.228) 2.264 1.360 ( 2.973 -1.114 1.114) 3.365 1.368 ( 3.435 -1.575 1.575) 4.094 1.601 ( -3.436 1.271 -1.271) 3.877 1.621 ( -2.703 0.264 -0.264) 2.729 1.746 ( -0.575 -0.053 0.053) 0.580 1.761 ( -0.378 -0.632 0.632) 0.970 1.915 ( -0.021 0.111 -0.111) 0.159 1.916 ( 0.189 0.025 -0.025) 0.192 1.976 ( -1.552 1.154 -1.154) 2.253 1.978 ( -1.210 0.662 -0.662) 1.530 2.096 ( -1.356 -0.721 0.721) 1.696 2.120 ( -1.582 -1.682 1.682) 2.856 2.183 ( 2.075 -1.425 1.425) 2.892 2.194 ( 1.388 -1.774 1.774) 2.868 2.235 ( 0.384 -0.989 0.989) 1.450 2.240 ( 0.583 -1.209 1.209) 1.807 2.862 ( 2.384 -2.063 -2.043) 3.757 2.862 ( 2.384 2.043 2.063) 3.757 2.943 ( 0.744 0.958 -0.958) 1.546 2.958 ( 1.746 0.102 -0.102) 1.752 3.147 ( -1.748 0.693 -0.693) 2.004 3.148 ( -2.619 0.533 -0.533) 2.726 3.207 ( -0.891 0.429 -0.429) 1.077 3.209 ( 0.499 0.528 -0.528) 0.899 3.348 ( -2.367 -0.061 0.061) 2.368 3.368 ( -1.558 -1.118 1.118) 2.220 3.617 ( 2.348 0.464 -0.464) 2.438 3.634 ( 2.994 -0.648 0.648) 3.131 3.692 ( 0.780 0.710 -0.710) 1.271 3.707 ( 1.618 -0.338 0.338) 1.687 3.785 ( 0.968 1.208 -1.208) 1.963 3.798 ( 1.067 0.476 -0.476) 1.262 3.884 ( -1.900 2.123 -2.123) 3.554 3.887 ( -1.716 2.103 -2.103) 3.434 4.036 ( -2.628 0.224 -0.224) 2.647 4.050 ( -1.948 -0.622 0.622) 2.137 4.113 ( -1.695 2.048 -2.048) 3.356 4.177 ( 1.738 -1.574 1.574) 2.824 4.309 ( 1.190 1.674 -1.674) 2.650 4.318 ( 1.576 1.121 -1.121) 2.235 4.430 ( -1.296 0.018 -0.018) 1.297 4.435 ( -0.726 -1.028 1.028) 1.625 4.464 ( -0.620 0.232 -0.232) 0.702 4.485 ( 0.325 0.189 -0.189) 0.421 4.528 ( 0.064 -1.706 1.706) 2.413 4.542 ( 0.627 -1.799 1.799) 2.620 4.662 ( -2.536 0.125 -0.125) 2.542 4.679 ( -2.047 -0.713 0.713) 2.282 4.690 ( -1.386 0.425 -0.425) 1.511 4.701 ( -1.913 -0.796 0.796) 2.220 4.870 ( 3.875 1.747 -1.747) 4.595 4.895 ( 3.699 0.391 -0.391) 3.740 5.087 ( 1.260 -0.192 0.192) 1.289 5.088 ( 1.442 -0.535 0.535) 1.629 5.256 ( -3.098 1.521 -1.521) 3.772 5.301 ( -0.740 -1.362 1.362) 2.064 5.382 ( 0.593 0.216 -0.216) 0.668 5.390 ( 0.633 -0.278 0.278) 0.745 ======================= Grid point 17 (8/20) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.254 ( 3.846 3.846 8.608) 10.182 1.257 ( 3.882 3.882 8.626) 10.225 1.461 ( -0.553 -0.553 -2.372) 2.498 1.464 ( -1.221 -1.221 -2.430) 2.981 1.528 ( -1.933 -1.933 -3.444) 4.397 1.555 ( -2.335 -2.335 -2.361) 4.060 1.754 ( -0.502 -0.502 1.595) 1.746 1.758 ( -1.189 -1.189 1.140) 2.032 1.904 ( 0.465 0.465 2.251) 2.345 1.930 ( 0.540 0.540 -1.411) 1.604 1.946 ( -0.613 -0.613 -2.452) 2.601 1.956 ( -0.639 -0.639 -0.248) 0.937 2.096 ( -0.760 -0.760 2.420) 2.648 2.119 ( -2.315 -2.315 1.696) 3.687 2.211 ( 2.022 2.022 5.344) 6.061 2.220 ( -0.014 -0.014 2.438) 2.438 2.242 ( 0.128 0.128 2.936) 2.942 2.253 ( -0.618 -0.618 1.418) 1.665 2.873 ( 2.487 2.487 1.763) 3.935 2.911 ( 1.545 1.545 -0.187) 2.194 2.926 ( 1.167 1.167 -1.503) 2.232 2.982 ( 0.874 0.874 -2.712) 2.981 3.101 ( -0.183 -0.183 1.342) 1.367 3.154 ( -3.245 -3.245 -0.888) 4.674 3.185 ( -0.539 -0.539 -3.899) 3.973 3.211 ( 0.823 0.823 -2.843) 3.072 3.325 ( -1.616 -1.616 2.296) 3.239 3.369 ( -2.940 -2.940 -1.011) 4.279 3.623 ( 3.257 3.257 3.581) 5.834 3.650 ( 2.933 2.933 3.512) 5.435 3.714 ( 0.447 0.447 -2.979) 3.045 3.757 ( 0.729 0.729 -1.375) 1.718 3.782 ( 0.296 0.296 -2.598) 2.632 3.808 ( 0.667 0.667 -2.550) 2.718 3.827 ( -0.254 -0.254 -4.130) 4.145 3.845 ( -0.989 -0.989 -4.333) 4.553 4.013 ( -0.059 -0.059 1.573) 1.576 4.036 ( -0.673 -0.673 -0.352) 1.015 4.078 ( -2.013 -2.013 -2.902) 4.065 4.208 ( -0.985 -0.985 1.872) 2.334 4.300 ( 2.729 2.729 -1.198) 4.041 4.331 ( 1.802 1.802 -3.244) 4.126 4.407 ( -0.965 -0.965 -0.922) 1.647 4.443 ( 0.071 0.071 2.008) 2.010 4.459 ( -0.486 -0.486 1.649) 1.787 4.483 ( -0.460 -0.460 -0.332) 0.730 4.541 ( 0.003 0.003 2.683) 2.683 4.570 ( -1.667 -1.667 0.588) 2.430 4.645 ( -1.175 -1.175 2.855) 3.304 4.656 ( -0.528 -0.528 -0.688) 1.015 4.681 ( -0.784 -0.784 -1.218) 1.647 4.686 ( -2.754 -2.754 1.059) 4.036 4.895 ( 3.440 3.440 -1.596) 5.121 4.913 ( 2.533 2.533 -1.847) 4.030 5.079 ( 2.073 2.073 3.744) 4.755 5.162 ( -0.386 -0.386 -2.611) 2.668 5.184 ( -1.331 -1.331 -1.390) 2.340 5.307 ( -1.516 -1.516 2.090) 2.994 5.383 ( 1.056 1.056 0.315) 1.527 5.399 ( -0.111 -0.111 -0.253) 0.298 ======================= Grid point 64 (9/20) ======================= q-point: ( 0.14 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 44 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.966 ( 9.355 9.355 9.355) 16.202 0.966 ( 9.355 9.355 9.355) 16.202 1.554 ( -1.797 -1.797 -1.797) 3.113 1.554 ( -1.797 -1.797 -1.797) 3.113 1.569 ( 0.395 0.395 0.395) 0.685 1.583 ( -1.107 -1.107 -1.107) 1.917 1.724 ( 8.983 8.983 8.983) 15.558 1.742 ( 0.180 0.180 0.180) 0.312 1.742 ( 0.180 0.180 0.180) 0.312 1.963 ( 0.210 0.210 0.210) 0.363 1.963 ( 0.210 0.210 0.210) 0.363 1.982 ( 5.251 5.251 5.251) 9.096 1.996 ( -0.979 -0.979 -0.979) 1.696 2.100 ( -0.904 -0.904 -0.904) 1.566 2.100 ( -0.904 -0.904 -0.904) 1.566 2.164 ( 1.005 1.005 1.005) 1.742 2.164 ( 1.005 1.005 1.005) 1.742 2.215 ( 0.928 0.928 0.928) 1.608 2.860 ( -0.573 -0.573 -0.573) 0.992 2.913 ( 0.407 0.407 0.407) 0.706 2.913 ( 0.407 0.407 0.407) 0.706 3.044 ( -0.714 -0.714 -0.714) 1.237 3.160 ( -1.389 -1.389 -1.389) 2.405 3.160 ( -1.389 -1.389 -1.389) 2.405 3.290 ( -1.741 -1.741 -1.741) 3.015 3.290 ( -1.741 -1.741 -1.741) 3.015 3.292 ( -1.292 -1.292 -1.292) 2.238 3.450 ( -3.673 -3.673 -3.673) 6.362 3.529 ( 2.611 2.611 2.611) 4.523 3.529 ( 2.611 2.611 2.611) 4.523 3.747 ( 1.080 1.080 1.080) 1.871 3.829 ( -0.583 -0.583 -0.583) 1.010 3.829 ( -0.583 -0.583 -0.583) 1.010 3.901 ( -1.854 -1.854 -1.854) 3.211 3.901 ( -1.854 -1.854 -1.854) 3.211 3.953 ( -2.685 -2.685 -2.685) 4.651 4.007 ( -1.768 -1.768 -1.768) 3.062 4.077 ( -1.605 -1.605 -1.605) 2.780 4.151 ( -1.238 -1.238 -1.238) 2.145 4.151 ( -1.238 -1.238 -1.238) 2.145 4.355 ( 0.178 0.178 0.178) 0.309 4.368 ( 1.218 1.218 1.218) 2.110 4.368 ( 1.218 1.218 1.218) 2.110 4.432 ( 3.141 3.141 3.141) 5.440 4.456 ( -1.024 -1.024 -1.024) 1.773 4.456 ( -1.024 -1.024 -1.024) 1.773 4.510 ( -0.710 -0.710 -0.710) 1.230 4.510 ( -0.710 -0.710 -0.710) 1.230 4.650 ( -1.900 -1.900 -1.900) 3.291 4.684 ( -0.911 -0.911 -0.911) 1.578 4.684 ( -0.911 -0.911 -0.911) 1.578 4.715 ( -0.320 -0.320 -0.320) 0.553 4.913 ( 1.461 1.461 1.461) 2.530 4.913 ( 1.461 1.461 1.461) 2.530 4.971 ( 2.769 2.769 2.769) 4.795 5.212 ( -0.314 -0.314 -0.314) 0.543 5.246 ( -0.589 -0.589 -0.589) 1.020 5.246 ( -0.589 -0.589 -0.589) 1.020 5.378 ( 0.689 0.689 0.689) 1.194 5.378 ( 0.689 0.689 0.689) 1.194 ======================= Grid point 65 (10/20) ======================= q-point: ( 0.29 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.218 ( 4.649 7.094 7.094) 11.057 1.219 ( 4.814 7.142 7.142) 11.189 1.461 ( 3.218 -4.052 -4.052) 6.572 1.462 ( 3.033 -4.166 -4.166) 6.626 1.568 ( -3.456 0.123 0.123) 3.460 1.574 ( -3.512 -0.142 -0.142) 3.518 1.728 ( -2.058 -0.484 -0.484) 2.169 1.729 ( -2.278 -0.592 -0.592) 2.427 1.879 ( 0.175 2.755 2.755) 3.900 1.938 ( 0.252 -1.530 -1.530) 2.178 1.958 ( -1.404 -0.292 -0.292) 1.464 1.968 ( -2.355 -0.966 -0.966) 2.722 2.080 ( -0.891 -0.125 -0.125) 0.908 2.080 ( -0.722 -0.102 -0.102) 0.736 2.182 ( 4.015 4.155 4.155) 7.116 2.207 ( 1.466 2.323 2.323) 3.597 2.209 ( 1.568 0.751 0.751) 1.894 2.242 ( 0.056 0.681 0.681) 0.964 2.879 ( 3.429 2.053 2.053) 4.493 2.930 ( 2.858 -0.523 -0.523) 2.952 2.943 ( 1.907 2.103 2.103) 3.532 3.008 ( 1.006 -1.568 -1.568) 2.435 3.105 ( -3.034 -1.021 -1.021) 3.360 3.121 ( -3.669 -1.665 -1.665) 4.360 3.227 ( 0.572 -1.453 -1.453) 2.132 3.233 ( 0.111 -1.874 -1.874) 2.653 3.292 ( -3.183 1.027 1.027) 3.499 3.351 ( -3.058 -2.901 -2.901) 5.117 3.621 ( 4.009 3.684 3.684) 6.574 3.643 ( 2.333 4.124 4.124) 6.282 3.739 ( 3.014 -2.410 -2.410) 4.550 3.776 ( 0.755 -3.349 -3.349) 4.796 3.799 ( 0.479 0.559 0.559) 0.924 3.841 ( 1.129 -1.729 -1.729) 2.694 3.870 ( -2.299 -0.670 -0.670) 2.487 3.881 ( -3.149 -1.746 -1.746) 4.002 3.988 ( -2.073 0.784 0.784) 2.351 4.030 ( -0.888 -1.674 -1.674) 2.528 4.113 ( -2.952 -0.165 -0.165) 2.962 4.132 ( -3.520 -0.942 -0.942) 3.764 4.352 ( 2.891 0.152 0.152) 2.899 4.362 ( 1.773 -0.525 -0.525) 1.922 4.412 ( -1.534 0.672 0.672) 1.804 4.426 ( 1.762 -0.557 -0.557) 1.930 4.465 ( 0.216 1.440 1.440) 2.048 4.497 ( 1.192 0.110 0.110) 1.203 4.517 ( -0.167 0.561 0.561) 0.811 4.525 ( 0.134 0.832 0.832) 1.185 4.606 ( -2.779 0.082 0.082) 2.782 4.648 ( -1.993 -1.589 -1.589) 3.003 4.666 ( -2.543 -0.145 -0.145) 2.551 4.686 ( -1.439 -1.315 -1.315) 2.352 4.939 ( 3.265 0.658 0.658) 3.395 4.945 ( 2.717 0.381 0.381) 2.770 5.065 ( 2.440 3.388 3.388) 5.377 5.160 ( 2.470 -2.974 -2.974) 4.877 5.236 ( -2.576 0.476 0.476) 2.662 5.248 ( -2.773 -0.190 -0.190) 2.786 5.391 ( 0.775 0.401 0.401) 0.960 5.396 ( 0.342 0.164 0.164) 0.413 ======================= Grid point 69 (11/20) ======================= q-point: (-0.14 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 40 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.577 ( 18.982 -0.000 -0.000) 18.982 0.577 ( 18.982 -0.000 -0.000) 18.982 1.101 ( 35.193 -0.000 -0.000) 35.193 1.611 ( -0.829 0.000 0.000) 0.829 1.611 ( -0.829 0.000 0.000) 0.829 1.655 ( 0.462 -0.000 -0.000) 0.462 1.655 ( 0.462 -0.000 -0.000) 0.462 1.669 ( 0.794 -0.000 -0.000) 0.794 1.683 ( 2.038 -0.000 -0.000) 2.038 1.909 ( 1.050 -0.000 -0.000) 1.050 1.926 ( 3.086 -0.000 -0.000) 3.086 1.926 ( 3.086 -0.000 -0.000) 3.086 1.997 ( -1.478 0.000 0.000) 1.478 2.110 ( 0.498 -0.000 -0.000) 0.498 2.155 ( 0.707 -0.000 -0.000) 0.707 2.156 ( -1.520 0.000 0.000) 1.520 2.156 ( -1.520 0.000 0.000) 1.520 2.224 ( 2.099 -0.000 -0.000) 2.099 2.901 ( -1.807 0.000 0.000) 1.807 2.910 ( -2.614 0.000 0.000) 2.614 3.014 ( -0.307 0.000 0.000) 0.307 3.038 ( -1.746 0.000 0.000) 1.746 3.038 ( -1.746 0.000 0.000) 1.746 3.231 ( 6.172 -0.000 -0.000) 6.172 3.318 ( -2.921 0.000 0.000) 2.921 3.318 ( -2.921 0.000 0.000) 2.921 3.392 ( 1.616 -0.000 -0.000) 1.616 3.512 ( -4.353 0.000 0.000) 4.353 3.512 ( -4.353 0.000 0.000) 4.353 3.605 ( -6.126 0.000 0.000) 6.126 3.688 ( 0.284 -0.000 -0.000) 0.284 3.796 ( 0.331 -0.000 -0.000) 0.331 3.796 ( 0.331 -0.000 -0.000) 0.331 3.970 ( -4.911 0.000 0.000) 4.911 3.970 ( -4.911 0.000 0.000) 4.911 4.049 ( -4.045 0.000 0.000) 4.045 4.079 ( -1.248 0.000 0.000) 1.248 4.141 ( -2.864 0.000 0.000) 2.864 4.185 ( -0.662 0.000 0.000) 0.662 4.185 ( -0.662 0.000 0.000) 0.662 4.266 ( 2.522 -0.000 -0.000) 2.522 4.314 ( 5.303 -0.000 -0.000) 5.303 4.314 ( 5.303 -0.000 -0.000) 5.303 4.453 ( -4.964 0.000 0.000) 4.964 4.486 ( -2.661 0.000 0.000) 2.661 4.593 ( -0.651 0.000 0.000) 0.651 4.593 ( -0.651 0.000 0.000) 0.651 4.593 ( 3.284 -0.000 -0.000) 3.284 4.613 ( -0.211 0.000 0.000) 0.211 4.674 ( -1.402 0.000 0.000) 1.402 4.674 ( -1.402 0.000 0.000) 1.402 4.736 ( 0.725 -0.000 -0.000) 0.725 4.829 ( -1.937 0.000 0.000) 1.937 4.829 ( -1.937 0.000 0.000) 1.937 4.928 ( 3.907 -0.000 -0.000) 3.907 5.244 ( 1.340 -0.000 -0.000) 1.340 5.244 ( 1.340 -0.000 -0.000) 1.340 5.252 ( 1.264 -0.000 -0.000) 1.264 5.307 ( -2.543 0.000 0.000) 2.543 5.390 ( 2.045 -0.000 -0.000) 2.045 ======================= Grid point 72 (12/20) ======================= q-point: ( 0.29 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 44 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.318 ( 0.457 0.457 -0.517) 0.828 1.318 ( 0.497 0.497 -0.437) 0.828 1.482 ( 1.031 1.031 -1.031) 1.786 1.502 ( 0.007 0.007 -0.007) 0.012 1.505 ( 1.213 1.213 -1.354) 2.185 1.505 ( 1.307 1.307 -1.166) 2.185 1.701 ( -2.321 -2.321 2.106) 3.900 1.701 ( -2.178 -2.178 2.393) 3.900 1.920 ( -0.706 -0.706 -0.696) 1.217 1.920 ( 0.229 0.229 1.173) 1.217 1.925 ( -0.649 -0.649 0.649) 1.124 1.935 ( -0.287 -0.287 0.287) 0.496 2.077 ( -0.460 -0.460 0.264) 0.701 2.077 ( -0.329 -0.329 0.525) 0.701 2.227 ( 0.515 0.515 -0.515) 0.893 2.240 ( 0.019 0.019 -0.019) 0.033 2.245 ( -0.036 -0.036 0.034) 0.061 2.245 ( -0.034 -0.034 0.037) 0.061 2.959 ( 1.435 1.435 -2.985) 3.610 2.959 ( 2.469 2.469 -0.918) 3.610 2.970 ( 1.484 1.484 -1.484) 2.571 3.006 ( 0.413 0.413 -0.413) 0.715 3.064 ( -2.328 -2.328 1.259) 3.524 3.064 ( -1.615 -1.615 2.684) 3.524 3.204 ( 1.709 1.709 -1.709) 2.961 3.243 ( 0.742 0.742 -0.742) 1.285 3.282 ( -2.536 -2.536 0.558) 3.630 3.282 ( -1.217 -1.217 3.196) 3.630 3.708 ( 0.779 0.779 -0.929) 1.441 3.708 ( 0.879 0.879 -0.729) 1.441 3.730 ( 2.205 2.205 -2.205) 3.819 3.787 ( -0.162 -0.162 0.162) 0.280 3.792 ( -0.225 -0.225 -0.926) 0.979 3.792 ( 0.543 0.543 0.608) 0.979 3.842 ( -0.062 -0.062 -0.877) 0.881 3.842 ( 0.564 0.564 0.375) 0.881 4.009 ( -0.586 -0.586 -0.803) 1.154 4.009 ( 0.340 0.340 1.049) 1.154 4.035 ( -0.360 -0.360 0.360) 0.623 4.117 ( -4.348 -4.348 4.348) 7.531 4.370 ( 0.865 0.865 -1.378) 1.842 4.370 ( 1.207 1.207 -0.694) 1.842 4.418 ( 1.668 1.668 -1.668) 2.888 4.458 ( -0.456 -0.456 -2.061) 2.160 4.458 ( 1.222 1.222 1.296) 2.160 4.491 ( 0.101 0.101 -0.101) 0.175 4.536 ( -0.446 -0.446 0.241) 0.675 4.536 ( -0.309 -0.309 0.514) 0.675 4.601 ( -1.114 -1.114 1.114) 1.930 4.633 ( -1.881 -1.881 -0.255) 2.673 4.633 ( -0.457 -0.457 2.593) 2.673 4.639 ( -0.835 -0.835 0.835) 1.446 4.968 ( 1.109 1.109 -1.367) 2.080 4.968 ( 1.281 1.281 -1.022) 2.080 5.155 ( 1.184 1.184 -1.184) 2.051 5.158 ( -0.372 -0.372 -2.716) 2.767 5.158 ( 1.687 1.687 1.401) 2.767 5.242 ( -2.707 -2.707 2.707) 4.689 5.400 ( 0.085 0.085 -0.218) 0.249 5.400 ( 0.174 0.174 -0.041) 0.249 ======================= Grid point 77 (13/20) ======================= q-point: (-0.29 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.886 ( 16.193 -5.605 -0.000) 17.136 0.890 ( 16.066 -5.939 -0.000) 17.128 1.555 ( 2.213 2.733 -0.000) 3.517 1.567 ( -1.072 3.372 0.000) 3.538 1.633 ( 1.577 1.596 -0.000) 2.244 1.638 ( 8.740 -3.834 -0.000) 9.544 1.666 ( 7.399 -2.433 -0.000) 7.789 1.704 ( 1.145 -2.196 -0.000) 2.477 1.746 ( 4.742 -2.488 -0.000) 5.355 1.924 ( 1.467 0.473 -0.000) 1.541 1.935 ( -1.548 0.456 0.000) 1.614 1.977 ( 1.399 -1.748 -0.000) 2.239 1.986 ( -1.476 -1.717 0.000) 2.264 2.129 ( 3.081 1.314 -0.000) 3.349 2.138 ( 1.492 0.954 -0.000) 1.771 2.165 ( 1.806 -0.359 -0.000) 1.842 2.180 ( 2.690 -0.897 -0.000) 2.835 2.251 ( 2.829 1.266 -0.000) 3.099 2.861 ( -1.860 1.615 0.000) 2.464 2.884 ( -0.700 -0.014 0.000) 0.700 2.955 ( -1.446 2.742 0.000) 3.100 2.999 ( -3.454 -0.202 0.000) 3.460 3.073 ( -1.853 -4.582 -0.000) 4.942 3.256 ( -3.926 1.095 0.000) 4.076 3.269 ( -3.652 -0.206 0.000) 3.657 3.284 ( 0.826 2.121 0.000) 2.276 3.387 ( 0.787 2.899 0.000) 3.004 3.468 ( -1.167 0.579 0.000) 1.302 3.501 ( 0.762 -2.863 -0.000) 2.963 3.517 ( -1.906 0.910 0.000) 2.112 3.692 ( -0.214 -0.251 0.000) 0.329 3.779 ( -2.586 0.163 0.000) 2.591 3.796 ( -2.169 -1.596 0.000) 2.693 3.884 ( -4.905 2.052 0.000) 5.317 3.906 ( -4.644 -0.872 0.000) 4.725 3.970 ( -5.343 1.104 0.000) 5.456 4.041 ( -1.147 1.217 0.000) 1.672 4.089 ( -2.713 0.578 0.000) 2.773 4.169 ( -1.645 0.130 0.000) 1.650 4.189 ( -1.684 -1.557 0.000) 2.294 4.293 ( -0.910 -1.308 0.000) 1.594 4.312 ( -2.539 3.413 0.000) 4.254 4.392 ( 5.066 0.013 -0.000) 5.066 4.411 ( 1.696 -0.732 -0.000) 1.847 4.469 ( 0.069 -0.757 -0.000) 0.760 4.527 ( -0.821 3.509 0.000) 3.603 4.566 ( -1.100 1.724 0.000) 2.045 4.616 ( 3.006 3.820 -0.000) 4.861 4.654 ( -0.306 -2.984 -0.000) 3.000 4.679 ( -0.620 -1.632 -0.000) 1.746 4.680 ( 1.035 -0.253 -0.000) 1.066 4.736 ( 0.487 0.601 -0.000) 0.774 4.808 ( -1.878 -1.438 0.000) 2.365 4.823 ( -2.117 -3.007 0.000) 3.678 4.985 ( 3.985 0.049 -0.000) 3.985 5.214 ( -1.222 3.532 0.000) 3.737 5.266 ( 1.002 -0.506 -0.000) 1.123 5.277 ( 0.422 -0.743 -0.000) 0.855 5.313 ( 0.328 -2.913 -0.000) 2.932 5.409 ( 1.202 0.803 -0.000) 1.445 ======================= Grid point 85 (14/20) ======================= q-point: (-0.43 0.43 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.183 ( 11.911 -5.814 -0.000) 13.255 1.187 ( 11.740 -5.930 -0.000) 13.153 1.434 ( -3.684 7.580 0.000) 8.428 1.442 ( -4.003 7.933 0.000) 8.886 1.629 ( 0.296 -0.121 -0.000) 0.320 1.642 ( -0.506 0.615 0.000) 0.796 1.738 ( 0.619 -1.552 -0.000) 1.671 1.759 ( -0.129 0.009 0.000) 0.129 1.898 ( 0.391 -1.366 -0.000) 1.421 1.911 ( -1.430 -0.860 0.000) 1.669 1.969 ( -0.219 -0.880 -0.000) 0.906 1.974 ( -1.925 -1.174 0.000) 2.254 2.098 ( 7.677 -2.774 -0.000) 8.163 2.150 ( 2.550 2.369 -0.000) 3.481 2.188 ( 5.355 0.145 -0.000) 5.357 2.216 ( 5.956 0.173 -0.000) 5.958 2.221 ( 2.415 1.953 -0.000) 3.106 2.274 ( 2.587 1.791 -0.000) 3.146 2.829 ( 0.644 0.082 -0.000) 0.649 2.876 ( -1.793 1.749 0.000) 2.505 2.909 ( 0.332 -1.002 -0.000) 1.056 2.957 ( -3.094 -0.101 0.000) 3.095 3.112 ( -0.631 -3.612 -0.000) 3.666 3.177 ( -5.022 0.769 0.000) 5.081 3.211 ( -2.830 -0.397 0.000) 2.858 3.230 ( -1.479 2.704 0.000) 3.082 3.372 ( 2.063 1.222 -0.000) 2.397 3.411 ( -1.185 2.953 0.000) 3.182 3.557 ( 3.685 -2.646 -0.000) 4.536 3.607 ( 3.393 -4.161 -0.000) 5.369 3.664 ( -1.995 1.864 0.000) 2.730 3.746 ( -3.202 0.948 0.000) 3.340 3.764 ( -2.238 0.749 0.000) 2.360 3.787 ( -2.682 0.149 0.000) 2.686 3.862 ( -3.967 -2.173 0.000) 4.523 3.881 ( -5.124 -0.485 0.000) 5.147 4.035 ( -0.161 -1.960 -0.000) 1.967 4.064 ( -2.090 -1.136 0.000) 2.379 4.136 ( -1.491 0.441 0.000) 1.555 4.150 ( -2.752 0.947 0.000) 2.910 4.269 ( -1.841 -2.434 0.000) 3.052 4.288 ( -2.767 -1.428 0.000) 3.114 4.439 ( 1.603 0.158 -0.000) 1.611 4.447 ( 2.153 0.829 -0.000) 2.308 4.477 ( 1.606 0.314 -0.000) 1.636 4.518 ( 0.157 2.844 0.000) 2.848 4.533 ( 3.247 0.323 -0.000) 3.263 4.583 ( 0.489 3.648 0.000) 3.680 4.678 ( -0.586 -0.791 0.000) 0.984 4.680 ( -1.030 -0.923 0.000) 1.383 4.700 ( 1.126 -0.772 -0.000) 1.366 4.728 ( -0.469 1.011 0.000) 1.115 4.823 ( -0.207 -2.054 -0.000) 2.064 4.837 ( -1.745 -2.074 0.000) 2.710 5.048 ( 3.862 -1.110 -0.000) 4.018 5.131 ( -2.093 3.968 0.000) 4.486 5.285 ( 0.310 -0.606 -0.000) 0.681 5.289 ( -0.154 -0.481 -0.000) 0.505 5.363 ( 1.321 -1.869 -0.000) 2.288 5.400 ( -0.330 1.311 0.000) 1.352 ======================= Grid point 86 (15/20) ======================= q-point: (-0.29 0.43 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.138 ( 13.558 -3.503 3.503) 14.434 1.138 ( 13.435 -3.443 3.443) 14.290 1.514 ( -2.870 2.438 -2.438) 4.486 1.524 ( -2.995 2.181 -2.181) 4.299 1.605 ( -0.159 3.148 -3.148) 4.454 1.626 ( -0.536 2.253 -2.253) 3.232 1.734 ( 1.021 -1.521 1.521) 2.381 1.761 ( 1.150 -0.960 0.960) 1.779 1.888 ( 1.716 -0.884 0.884) 2.122 1.898 ( -2.768 -1.159 1.159) 3.217 1.960 ( 0.086 0.002 -0.002) 0.086 1.962 ( -2.554 0.231 -0.231) 2.575 2.071 ( 11.320 -2.447 2.447) 11.837 2.176 ( 3.817 1.208 -1.208) 4.182 2.179 ( 5.797 0.428 -0.428) 5.829 2.223 ( 4.165 0.113 -0.113) 4.168 2.238 ( 5.056 -0.237 0.237) 5.067 2.286 ( 3.318 1.438 -1.438) 3.892 2.841 ( -0.109 0.200 -0.200) 0.304 2.883 ( 0.406 -0.186 0.186) 0.483 2.915 ( -2.531 0.354 -0.354) 2.580 2.955 ( -3.911 -1.106 1.106) 4.212 3.074 ( 0.108 -1.737 1.737) 2.459 3.188 ( -4.855 0.038 -0.038) 4.855 3.199 ( -4.722 -0.101 0.101) 4.724 3.253 ( -2.576 2.406 -2.406) 4.268 3.380 ( 1.793 2.305 -2.305) 3.720 3.444 ( -0.510 2.107 -2.107) 3.023 3.543 ( 3.336 -2.427 2.427) 4.787 3.563 ( 3.834 -2.642 2.642) 5.354 3.690 ( -1.051 -0.440 0.440) 1.222 3.741 ( -2.778 -0.309 0.309) 2.812 3.767 ( -2.669 -1.045 1.045) 3.051 3.817 ( -4.961 0.799 -0.799) 5.088 3.828 ( -4.367 0.066 -0.066) 4.368 3.879 ( -5.860 0.897 -0.897) 5.996 4.013 ( -1.082 0.301 -0.301) 1.163 4.045 ( -2.119 0.251 -0.251) 2.148 4.127 ( -4.404 0.824 -0.824) 4.556 4.181 ( -1.621 -2.166 2.166) 3.465 4.248 ( -3.445 -0.510 0.510) 3.520 4.278 ( -2.800 -1.433 1.433) 3.456 4.442 ( 2.225 0.717 -0.717) 2.445 4.465 ( 4.065 0.052 -0.052) 4.066 4.465 ( 1.257 1.674 -1.674) 2.680 4.538 ( 2.780 0.069 -0.069) 2.782 4.539 ( 0.512 2.226 -2.226) 3.189 4.623 ( 0.172 1.862 -1.862) 2.639 4.656 ( -1.039 -0.928 0.928) 1.674 4.677 ( -0.788 -0.805 0.805) 1.385 4.683 ( 0.661 1.453 -1.453) 2.159 4.744 ( 1.316 0.886 -0.886) 1.817 4.795 ( -2.316 -2.767 2.767) 4.547 4.829 ( -1.321 -3.026 3.026) 4.478 5.045 ( 3.969 -0.341 0.341) 3.999 5.188 ( -1.617 1.561 -1.561) 2.737 5.231 ( -0.037 1.815 -1.815) 2.567 5.320 ( 1.579 -0.564 0.564) 1.769 5.335 ( 0.545 -2.169 2.169) 3.115 5.411 ( -0.289 0.628 -0.628) 0.934 ======================= Grid point 92 (16/20) ======================= q-point: (-0.43 0.57 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.323 ( 3.049 3.049 1.353) 4.519 1.325 ( 2.967 2.967 1.273) 4.384 1.419 ( 2.103 2.103 1.434) 3.302 1.427 ( 1.681 1.681 1.644) 2.891 1.577 ( 0.370 0.370 -6.034) 6.057 1.602 ( -0.147 -0.147 -4.861) 4.865 1.763 ( 0.604 0.604 1.557) 1.776 1.767 ( -0.220 -0.220 1.252) 1.290 1.894 ( -1.402 -1.402 0.993) 2.217 1.897 ( -1.425 -1.425 1.067) 2.280 1.930 ( -1.516 -1.516 -1.464) 2.596 1.960 ( -0.576 -0.576 0.124) 0.823 2.163 ( 3.363 3.363 -1.477) 4.980 2.184 ( 2.712 2.712 -1.709) 4.199 2.253 ( 2.696 2.696 0.648) 3.867 2.260 ( 2.548 2.548 -0.038) 3.604 2.294 ( 2.707 2.707 -1.247) 4.027 2.298 ( 2.536 2.536 -1.171) 3.773 2.857 ( -0.650 -0.650 1.601) 1.846 2.872 ( 0.704 0.704 2.125) 2.347 2.910 ( -0.753 -0.753 -0.465) 1.162 2.924 ( -2.115 -2.115 0.904) 3.125 3.108 ( -2.033 -2.033 0.802) 2.985 3.130 ( -0.773 -0.773 1.788) 2.096 3.174 ( -1.462 -1.462 -1.449) 2.525 3.187 ( -0.957 -0.957 -1.665) 2.146 3.386 ( 2.587 2.587 -1.882) 4.115 3.410 ( 2.000 2.000 -2.028) 3.481 3.615 ( 0.458 0.458 2.318) 2.407 3.652 ( 0.377 0.377 3.320) 3.362 3.667 ( -0.924 -0.924 -0.156) 1.316 3.707 ( -0.498 -0.498 1.805) 1.938 3.746 ( -1.114 -1.114 1.065) 1.901 3.760 ( -2.104 -2.104 -0.177) 2.981 3.803 ( -2.915 -2.915 -1.471) 4.377 3.804 ( -2.690 -2.690 -1.234) 3.999 4.013 ( -0.943 -0.943 -0.028) 1.334 4.037 ( -1.412 -1.412 0.261) 2.014 4.062 ( -2.123 -2.123 -2.345) 3.810 4.180 ( -1.455 -1.455 4.011) 4.508 4.245 ( -2.626 -2.626 1.375) 3.960 4.259 ( -2.878 -2.878 1.568) 4.362 4.442 ( 1.734 1.734 -2.347) 3.395 4.475 ( 1.997 1.997 0.034) 2.824 4.499 ( 1.717 1.717 -1.868) 3.063 4.522 ( 2.280 2.280 -1.347) 3.494 4.578 ( 1.399 1.399 -0.079) 1.980 4.590 ( 1.287 1.287 -0.151) 1.826 4.652 ( -0.665 -0.665 -0.583) 1.107 4.674 ( -0.581 -0.581 0.111) 0.830 4.680 ( -0.323 -0.323 -1.424) 1.495 4.744 ( 1.110 1.110 -0.419) 1.624 4.808 ( -2.297 -2.297 3.075) 4.473 4.840 ( -3.039 -3.039 2.086) 4.778 5.104 ( 1.267 1.267 0.048) 1.793 5.132 ( 2.397 2.397 1.778) 3.828 5.222 ( -0.727 -0.727 -4.622) 4.735 5.332 ( 0.830 0.830 1.869) 2.207 5.376 ( -0.361 -0.361 -0.122) 0.525 5.388 ( -0.572 -0.572 0.130) 0.820 ======================= Grid point 137 (17/20) ======================= q-point: (-0.29 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 40 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.083 ( 15.162 -0.000 -0.000) 15.162 1.083 ( 15.162 -0.000 -0.000) 15.162 1.561 ( -2.632 0.000 0.000) 2.632 1.561 ( -2.632 0.000 0.000) 2.632 1.681 ( 1.345 -0.000 -0.000) 1.345 1.681 ( 1.345 -0.000 -0.000) 1.345 1.682 ( 0.422 -0.000 -0.000) 0.422 1.720 ( 1.789 -0.000 -0.000) 1.789 1.865 ( 3.859 -0.000 -0.000) 3.859 1.880 ( -4.573 0.000 0.000) 4.573 1.962 ( -1.529 0.000 0.000) 1.529 1.962 ( -1.529 0.000 0.000) 1.529 2.027 ( 15.002 -0.000 -0.000) 15.002 2.192 ( 4.822 -0.000 -0.000) 4.822 2.192 ( 4.822 -0.000 -0.000) 4.822 2.224 ( 4.964 -0.000 -0.000) 4.964 2.231 ( 5.688 -0.000 -0.000) 5.688 2.311 ( 3.789 -0.000 -0.000) 3.789 2.854 ( -0.318 0.000 0.000) 0.318 2.890 ( 1.253 -0.000 -0.000) 1.253 2.929 ( -5.209 0.000 0.000) 5.209 2.929 ( -5.209 0.000 0.000) 5.209 3.034 ( 1.867 -0.000 -0.000) 1.867 3.191 ( -5.479 0.000 0.000) 5.479 3.191 ( -5.479 0.000 0.000) 5.479 3.291 ( -3.933 0.000 0.000) 3.933 3.442 ( 2.109 -0.000 -0.000) 2.109 3.481 ( 2.117 -0.000 -0.000) 2.117 3.481 ( 2.117 -0.000 -0.000) 2.117 3.516 ( 2.083 -0.000 -0.000) 2.083 3.696 ( 0.457 -0.000 -0.000) 0.457 3.738 ( -3.353 0.000 0.000) 3.353 3.738 ( -3.353 0.000 0.000) 3.353 3.830 ( -4.598 0.000 0.000) 4.598 3.830 ( -4.598 0.000 0.000) 4.598 3.891 ( -6.032 0.000 0.000) 6.032 4.024 ( -2.300 0.000 0.000) 2.300 4.058 ( -2.351 0.000 0.000) 2.351 4.131 ( -3.247 0.000 0.000) 3.247 4.131 ( -3.247 0.000 0.000) 3.247 4.251 ( -8.075 0.000 0.000) 8.075 4.257 ( -2.663 0.000 0.000) 2.663 4.466 ( 4.729 -0.000 -0.000) 4.729 4.466 ( 4.729 -0.000 -0.000) 4.729 4.473 ( 1.825 -0.000 -0.000) 1.825 4.575 ( 0.594 -0.000 -0.000) 0.594 4.575 ( 0.594 -0.000 -0.000) 0.594 4.605 ( 0.026 -0.000 -0.000) 0.026 4.661 ( -0.373 0.000 0.000) 0.373 4.661 ( -0.373 0.000 0.000) 0.373 4.717 ( 0.518 -0.000 -0.000) 0.518 4.757 ( -2.701 0.000 0.000) 2.701 4.757 ( -2.701 0.000 0.000) 2.701 4.757 ( 2.343 -0.000 -0.000) 2.343 5.043 ( 3.998 -0.000 -0.000) 3.998 5.239 ( -0.977 0.000 0.000) 0.977 5.277 ( 0.742 -0.000 -0.000) 0.742 5.277 ( 0.742 -0.000 -0.000) 0.742 5.292 ( 1.116 -0.000 -0.000) 1.116 5.422 ( -0.307 0.000 0.000) 0.307 ======================= Grid point 144 (18/20) ======================= q-point: (-0.43 0.43 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.297 ( 10.883 -2.034 -0.000) 11.072 1.298 ( 10.866 -2.146 -0.000) 11.076 1.485 ( -2.763 4.656 0.000) 5.414 1.489 ( -3.105 4.395 0.000) 5.381 1.659 ( 1.611 2.633 0.000) 3.087 1.663 ( 1.208 2.573 0.000) 2.842 1.727 ( 1.002 -1.462 -0.000) 1.772 1.752 ( -0.479 -0.228 0.000) 0.530 1.848 ( -2.141 -2.963 0.000) 3.655 1.868 ( -2.384 -1.289 0.000) 2.710 1.936 ( -2.038 -0.666 0.000) 2.144 1.936 ( -2.136 -0.717 0.000) 2.253 2.204 ( 7.313 0.907 -0.000) 7.369 2.244 ( 3.972 3.276 -0.000) 5.149 2.272 ( 5.346 0.321 -0.000) 5.356 2.287 ( 4.009 0.813 -0.000) 4.090 2.326 ( 5.978 1.093 -0.000) 6.077 2.350 ( 3.103 1.966 -0.000) 3.673 2.843 ( -1.101 0.422 0.000) 1.180 2.857 ( -3.974 -0.395 0.000) 3.993 2.891 ( 0.948 -1.507 -0.000) 1.781 2.918 ( 0.981 0.370 -0.000) 1.049 3.069 ( -0.206 -1.410 -0.000) 1.425 3.110 ( -4.878 0.256 0.000) 4.884 3.142 ( -2.279 -3.282 0.000) 3.996 3.207 ( -5.282 0.849 0.000) 5.349 3.448 ( 4.348 3.094 -0.000) 5.336 3.479 ( 2.680 4.622 0.000) 5.343 3.571 ( 3.827 -3.164 -0.000) 4.965 3.597 ( 3.406 -3.528 -0.000) 4.903 3.664 ( -1.885 1.714 0.000) 2.547 3.703 ( -2.569 -1.397 0.000) 2.924 3.723 ( -1.927 -1.203 0.000) 2.271 3.753 ( -3.383 1.420 0.000) 3.669 3.775 ( -3.700 -1.355 0.000) 3.941 3.804 ( -4.913 0.369 0.000) 4.927 3.995 ( -1.335 -1.002 0.000) 1.669 4.019 ( -2.336 0.085 0.000) 2.338 4.086 ( -3.381 -1.151 0.000) 3.572 4.093 ( -4.294 -1.218 0.000) 4.463 4.178 ( -4.158 -4.502 0.000) 6.129 4.214 ( -4.716 -1.539 0.000) 4.961 4.494 ( 3.488 1.930 -0.000) 3.986 4.510 ( 2.701 2.830 -0.000) 3.912 4.529 ( 3.979 0.190 -0.000) 3.983 4.569 ( 1.703 1.885 -0.000) 2.540 4.583 ( 2.258 0.619 -0.000) 2.341 4.622 ( -0.457 1.600 0.000) 1.663 4.650 ( -0.443 -1.182 -0.000) 1.262 4.654 ( -0.975 -0.541 0.000) 1.115 4.696 ( -0.902 -0.368 0.000) 0.974 4.717 ( -2.189 0.190 0.000) 2.198 4.761 ( -1.999 -3.843 -0.000) 4.332 4.782 ( -1.924 -1.659 0.000) 2.541 5.110 ( 5.397 0.179 -0.000) 5.400 5.187 ( 0.328 5.258 0.000) 5.268 5.284 ( 0.306 -0.117 -0.000) 0.327 5.285 ( 0.184 -0.108 -0.000) 0.214 5.344 ( 1.003 -3.333 -0.000) 3.481 5.399 ( -2.087 0.768 0.000) 2.224 ======================= Grid point 152 (19/20) ======================= q-point: (-0.43 -0.43 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 44 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.414 ( 1.970 1.970 -2.075) 3.473 1.414 ( 2.040 2.040 -1.935) 3.473 1.442 ( 0.450 0.450 -0.699) 0.945 1.442 ( 0.616 0.616 -0.366) 0.945 1.648 ( 2.964 2.964 -2.964) 5.133 1.653 ( 4.035 4.035 -4.035) 6.988 1.756 ( -1.123 -1.123 0.044) 1.589 1.756 ( -0.403 -0.403 1.483) 1.589 1.844 ( -2.180 -2.180 1.795) 3.567 1.844 ( -1.923 -1.923 2.308) 3.567 1.895 ( -0.583 -0.583 0.583) 1.009 1.930 ( -1.548 -1.548 1.548) 2.681 2.267 ( 2.643 2.643 -3.794) 5.326 2.267 ( 3.410 3.410 -2.259) 5.326 2.318 ( 1.628 1.628 -1.628) 2.820 2.330 ( 1.900 1.900 -1.900) 3.291 2.366 ( 1.689 1.689 -1.914) 3.061 2.366 ( 1.839 1.839 -1.615) 3.061 2.814 ( -1.612 -1.612 1.612) 2.793 2.835 ( -2.668 -2.668 2.668) 4.620 2.925 ( 1.244 1.244 -2.390) 2.967 2.925 ( 2.008 2.008 -0.862) 2.967 3.045 ( -1.623 -1.623 1.623) 2.810 3.091 ( -1.727 -1.727 1.727) 2.990 3.139 ( -2.394 -2.394 0.340) 3.403 3.139 ( -1.025 -1.025 3.079) 3.403 3.492 ( 3.362 3.362 -4.753) 6.723 3.492 ( 4.290 4.290 -2.899) 6.723 3.636 ( -1.246 -1.246 0.422) 1.812 3.636 ( -0.697 -0.697 1.521) 1.812 3.646 ( 0.697 0.697 -0.697) 1.208 3.692 ( -1.715 -1.715 0.273) 2.441 3.692 ( -0.753 -0.753 2.196) 2.441 3.714 ( -1.091 -1.091 1.091) 1.890 3.750 ( -1.311 -1.311 -0.021) 1.854 3.750 ( -0.423 -0.423 1.755) 1.854 3.997 ( -0.652 -0.652 0.652) 1.130 3.998 ( -1.000 -1.000 0.575) 1.526 3.998 ( -0.717 -0.717 1.141) 1.526 4.115 ( -2.797 -2.797 2.797) 4.844 4.154 ( -4.142 -4.142 2.533) 6.382 4.154 ( -3.069 -3.069 4.678) 6.382 4.516 ( 2.884 2.884 -2.884) 4.995 4.550 ( 1.698 1.698 -3.122) 3.939 4.550 ( 2.647 2.647 -1.224) 3.939 4.581 ( 1.452 1.452 -1.452) 2.514 4.611 ( 0.149 0.149 -1.512) 1.527 4.611 ( 1.058 1.058 0.305) 1.527 4.635 ( -0.347 -0.347 0.347) 0.601 4.661 ( -0.912 -0.912 -0.281) 1.321 4.661 ( -0.117 -0.117 1.310) 1.321 4.742 ( -3.270 -3.270 1.122) 4.758 4.742 ( -1.838 -1.838 3.986) 4.758 4.743 ( -1.362 -1.362 1.362) 2.359 5.169 ( 2.732 2.732 -2.732) 4.733 5.214 ( 0.300 0.300 -4.246) 4.267 5.214 ( 2.931 2.931 1.015) 4.267 5.351 ( -0.499 -0.499 0.499) 0.865 5.352 ( -1.916 -1.916 -0.384) 2.736 5.352 ( -0.382 -0.382 2.682) 2.736 ======================= Grid point 208 (20/20) ======================= q-point: (-0.43 0.43 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 6.61e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 40 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.417 ( 6.554 -0.000 -0.000) 6.554 1.417 ( 6.554 -0.000 -0.000) 6.554 1.481 ( -2.041 0.000 0.000) 2.041 1.481 ( -2.041 0.000 0.000) 2.041 1.718 ( 2.273 -0.000 -0.000) 2.273 1.735 ( 2.347 -0.000 -0.000) 2.347 1.735 ( 2.347 -0.000 -0.000) 2.347 1.747 ( 0.635 -0.000 -0.000) 0.635 1.781 ( -2.074 0.000 0.000) 2.074 1.810 ( -3.573 0.000 0.000) 3.573 1.891 ( -2.711 0.000 0.000) 2.711 1.891 ( -2.711 0.000 0.000) 2.711 2.311 ( 5.591 -0.000 -0.000) 5.591 2.338 ( 2.988 -0.000 -0.000) 2.988 2.345 ( 3.910 -0.000 -0.000) 3.910 2.345 ( 3.910 -0.000 -0.000) 3.910 2.387 ( 0.886 -0.000 -0.000) 0.886 2.395 ( 0.440 -0.000 -0.000) 0.440 2.789 ( -3.350 0.000 0.000) 3.350 2.789 ( -3.350 0.000 0.000) 3.350 2.943 ( 5.614 -0.000 -0.000) 5.614 2.972 ( 3.938 -0.000 -0.000) 3.938 3.052 ( -3.025 0.000 0.000) 3.025 3.052 ( -3.025 0.000 0.000) 3.025 3.068 ( -1.294 0.000 0.000) 1.294 3.121 ( -6.343 0.000 0.000) 6.343 3.583 ( 3.949 -0.000 -0.000) 3.949 3.583 ( 3.949 -0.000 -0.000) 3.949 3.585 ( 7.475 -0.000 -0.000) 7.475 3.624 ( 4.203 -0.000 -0.000) 4.203 3.642 ( -2.379 0.000 0.000) 2.379 3.642 ( -2.379 0.000 0.000) 2.379 3.707 ( -0.505 0.000 0.000) 0.505 3.715 ( -3.357 0.000 0.000) 3.357 3.715 ( -3.357 0.000 0.000) 3.357 3.743 ( -4.051 0.000 0.000) 4.051 3.975 ( -0.014 0.000 0.000) 0.014 4.012 ( -0.663 0.000 0.000) 0.663 4.031 ( -2.473 0.000 0.000) 2.473 4.031 ( -2.473 0.000 0.000) 2.473 4.047 ( -4.241 0.000 0.000) 4.241 4.111 ( -7.090 0.000 0.000) 7.090 4.575 ( 4.695 -0.000 -0.000) 4.695 4.585 ( 3.434 -0.000 -0.000) 3.434 4.585 ( 3.434 -0.000 -0.000) 3.434 4.613 ( 1.218 -0.000 -0.000) 1.218 4.613 ( 1.218 -0.000 -0.000) 1.218 4.625 ( 1.315 -0.000 -0.000) 1.315 4.638 ( -0.760 0.000 0.000) 0.760 4.638 ( -0.760 0.000 0.000) 0.760 4.676 ( -2.029 0.000 0.000) 2.029 4.681 ( -2.608 0.000 0.000) 2.608 4.681 ( -2.608 0.000 0.000) 2.608 4.720 ( -4.548 0.000 0.000) 4.548 5.201 ( 6.322 -0.000 -0.000) 6.322 5.265 ( 1.795 -0.000 -0.000) 1.795 5.288 ( 0.132 -0.000 -0.000) 0.132 5.288 ( 0.132 -0.000 -0.000) 0.132 5.302 ( -0.495 0.000 0.000) 0.495 5.359 ( -4.012 0.000 0.000) 4.012 =================== 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/20580 10.0 749.918 749.918 749.918 0.000 -0.000 -0.000 3/20580 20.0 122.898 122.898 122.898 0.000 -0.000 -0.000 3/20580 30.0 43.363 43.363 43.363 0.000 -0.000 -0.000 3/20580 40.0 22.149 22.149 22.149 0.000 -0.000 -0.000 3/20580 50.0 14.024 14.024 14.024 0.000 -0.000 -0.000 3/20580 60.0 10.094 10.094 10.094 0.000 -0.000 -0.000 3/20580 70.0 7.866 7.866 7.866 0.000 -0.000 -0.000 3/20580 80.0 6.453 6.453 6.453 0.000 -0.000 -0.000 3/20580 90.0 5.482 5.482 5.482 0.000 -0.000 -0.000 3/20580 100.0 4.775 4.775 4.775 0.000 -0.000 -0.000 3/20580 110.0 4.236 4.236 4.236 0.000 -0.000 -0.000 3/20580 120.0 3.812 3.812 3.812 0.000 -0.000 -0.000 3/20580 130.0 3.468 3.468 3.468 0.000 -0.000 -0.000 3/20580 140.0 3.183 3.183 3.183 0.000 -0.000 -0.000 3/20580 150.0 2.944 2.944 2.944 0.000 -0.000 -0.000 3/20580 160.0 2.739 2.739 2.739 0.000 -0.000 -0.000 3/20580 170.0 2.562 2.562 2.562 0.000 -0.000 -0.000 3/20580 180.0 2.407 2.407 2.407 0.000 -0.000 -0.000 3/20580 190.0 2.270 2.270 2.270 0.000 -0.000 -0.000 3/20580 200.0 2.148 2.148 2.148 0.000 -0.000 -0.000 3/20580 210.0 2.039 2.039 2.039 0.000 -0.000 -0.000 3/20580 220.0 1.941 1.941 1.941 0.000 -0.000 -0.000 3/20580 230.0 1.852 1.852 1.852 0.000 -0.000 -0.000 3/20580 240.0 1.771 1.771 1.771 0.000 -0.000 -0.000 3/20580 250.0 1.697 1.697 1.697 0.000 -0.000 -0.000 3/20580 260.0 1.629 1.629 1.629 0.000 -0.000 -0.000 3/20580 270.0 1.567 1.567 1.567 0.000 -0.000 -0.000 3/20580 280.0 1.509 1.509 1.509 0.000 -0.000 -0.000 3/20580 290.0 1.455 1.455 1.455 0.000 -0.000 -0.000 3/20580 300.0 1.405 1.405 1.405 0.000 -0.000 -0.000 3/20580 310.0 1.358 1.358 1.358 0.000 -0.000 -0.000 3/20580 320.0 1.315 1.315 1.315 0.000 -0.000 -0.000 3/20580 330.0 1.274 1.274 1.274 0.000 -0.000 -0.000 3/20580 340.0 1.236 1.236 1.236 0.000 -0.000 -0.000 3/20580 350.0 1.200 1.200 1.200 0.000 -0.000 -0.000 3/20580 360.0 1.166 1.166 1.166 0.000 -0.000 -0.000 3/20580 370.0 1.133 1.133 1.133 0.000 -0.000 -0.000 3/20580 380.0 1.103 1.103 1.103 0.000 -0.000 -0.000 3/20580 390.0 1.074 1.074 1.074 0.000 -0.000 -0.000 3/20580 400.0 1.047 1.047 1.047 0.000 -0.000 -0.000 3/20580 410.0 1.021 1.021 1.021 0.000 -0.000 -0.000 3/20580 420.0 0.996 0.996 0.996 0.000 -0.000 -0.000 3/20580 430.0 0.973 0.973 0.973 0.000 -0.000 -0.000 3/20580 440.0 0.951 0.951 0.951 0.000 -0.000 -0.000 3/20580 450.0 0.929 0.929 0.929 0.000 -0.000 -0.000 3/20580 460.0 0.909 0.909 0.909 0.000 -0.000 -0.000 3/20580 470.0 0.889 0.889 0.889 0.000 -0.000 -0.000 3/20580 480.0 0.870 0.870 0.870 0.000 -0.000 -0.000 3/20580 490.0 0.852 0.852 0.852 0.000 -0.000 -0.000 3/20580 500.0 0.835 0.835 0.835 0.000 -0.000 -0.000 3/20580 510.0 0.819 0.819 0.819 0.000 -0.000 -0.000 3/20580 520.0 0.803 0.803 0.803 0.000 -0.000 -0.000 3/20580 530.0 0.788 0.788 0.788 0.000 -0.000 -0.000 3/20580 540.0 0.773 0.773 0.773 0.000 -0.000 -0.000 3/20580 550.0 0.759 0.759 0.759 0.000 -0.000 -0.000 3/20580 560.0 0.745 0.745 0.745 0.000 -0.000 -0.000 3/20580 570.0 0.732 0.732 0.732 0.000 -0.000 -0.000 3/20580 580.0 0.719 0.719 0.719 0.000 -0.000 -0.000 3/20580 590.0 0.707 0.707 0.707 0.000 -0.000 -0.000 3/20580 600.0 0.695 0.695 0.695 0.000 -0.000 -0.000 3/20580 610.0 0.684 0.684 0.684 0.000 -0.000 -0.000 3/20580 620.0 0.673 0.673 0.673 0.000 -0.000 -0.000 3/20580 630.0 0.662 0.662 0.662 0.000 -0.000 -0.000 3/20580 640.0 0.651 0.651 0.651 0.000 -0.000 -0.000 3/20580 650.0 0.641 0.641 0.641 0.000 -0.000 -0.000 3/20580 660.0 0.632 0.632 0.632 0.000 -0.000 -0.000 3/20580 670.0 0.622 0.622 0.622 0.000 -0.000 -0.000 3/20580 680.0 0.613 0.613 0.613 0.000 -0.000 -0.000 3/20580 690.0 0.604 0.604 0.604 0.000 -0.000 -0.000 3/20580 700.0 0.595 0.595 0.595 0.000 -0.000 -0.000 3/20580 710.0 0.587 0.587 0.587 0.000 -0.000 -0.000 3/20580 720.0 0.579 0.579 0.579 0.000 -0.000 -0.000 3/20580 730.0 0.571 0.571 0.571 0.000 -0.000 -0.000 3/20580 740.0 0.563 0.563 0.563 0.000 -0.000 -0.000 3/20580 750.0 0.556 0.556 0.556 0.000 -0.000 -0.000 3/20580 760.0 0.548 0.548 0.548 0.000 -0.000 -0.000 3/20580 770.0 0.541 0.541 0.541 0.000 -0.000 -0.000 3/20580 780.0 0.534 0.534 0.534 0.000 -0.000 -0.000 3/20580 790.0 0.527 0.527 0.527 0.000 -0.000 -0.000 3/20580 800.0 0.521 0.521 0.521 0.000 -0.000 -0.000 3/20580 810.0 0.514 0.514 0.514 0.000 -0.000 -0.000 3/20580 820.0 0.508 0.508 0.508 0.000 -0.000 -0.000 3/20580 830.0 0.502 0.502 0.502 0.000 -0.000 -0.000 3/20580 840.0 0.496 0.496 0.496 0.000 -0.000 -0.000 3/20580 850.0 0.490 0.490 0.490 0.000 -0.000 -0.000 3/20580 860.0 0.484 0.484 0.484 0.000 -0.000 -0.000 3/20580 870.0 0.479 0.479 0.479 0.000 -0.000 -0.000 3/20580 880.0 0.473 0.473 0.473 0.000 -0.000 -0.000 3/20580 890.0 0.468 0.468 0.468 0.000 -0.000 -0.000 3/20580 900.0 0.463 0.463 0.463 0.000 -0.000 -0.000 3/20580 910.0 0.458 0.458 0.458 0.000 -0.000 -0.000 3/20580 920.0 0.453 0.453 0.453 0.000 -0.000 -0.000 3/20580 930.0 0.448 0.448 0.448 0.000 -0.000 -0.000 3/20580 940.0 0.443 0.443 0.443 0.000 -0.000 -0.000 3/20580 950.0 0.438 0.438 0.438 0.000 -0.000 -0.000 3/20580 960.0 0.434 0.434 0.434 0.000 -0.000 -0.000 3/20580 970.0 0.429 0.429 0.429 0.000 -0.000 -0.000 3/20580 980.0 0.425 0.425 0.425 0.000 -0.000 -0.000 3/20580 990.0 0.421 0.421 0.421 0.000 -0.000 -0.000 3/20580 1000.0 0.416 0.416 0.416 0.000 -0.000 -0.000 3/20580 Thermal conductivity related properties were written into "kappa-m777.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 00:18:43]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|