----------------------------- ------- 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: [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) Number of symmetry operations in supercell: 1296 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.765565529999999 0.000000000000000 0.000000000000000 b 0.000000000000000 3.765565529999999 0.000000000000000 c 0.000000000000000 0.000000000000000 3.765565529999999 Atomic positions (fractional): *1 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 2 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 3 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 *4 Re 0.00000000000000 0.00000000000000 0.00000000000000 186.207 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.765565529999999 0.000000000000000 0.000000000000000 b 0.000000000000000 3.765565529999999 0.000000000000000 c 0.000000000000000 0.000000000000000 3.765565529999999 Atomic positions (fractional): *1 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 > 1 2 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 2 3 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 > 3 *4 Re 0.00000000000000 0.00000000000000 0.00000000000000 186.207 > 4 -------------------------------- super cell -------------------------------- Lattice vectors: a 11.296696589999996 0.000000000000000 0.000000000000000 b 0.000000000000000 11.296696589999996 0.000000000000000 c 0.000000000000000 0.000000000000000 11.296696589999996 Atomic positions (fractional): *1 O 0.16666666666667 0.00000000000000 0.00000000000000 15.999 > 1 2 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 > 1 3 O 0.83333333333333 0.00000000000000 0.00000000000000 15.999 > 1 4 O 0.16666666666667 0.33333333333333 0.00000000000000 15.999 > 1 5 O 0.50000000000000 0.33333333333333 0.00000000000000 15.999 > 1 6 O 0.83333333333333 0.33333333333333 0.00000000000000 15.999 > 1 7 O 0.16666666666667 0.66666666666667 0.00000000000000 15.999 > 1 8 O 0.50000000000000 0.66666666666667 0.00000000000000 15.999 > 1 9 O 0.83333333333333 0.66666666666667 0.00000000000000 15.999 > 1 10 O 0.16666666666667 0.00000000000000 0.33333333333333 15.999 > 1 11 O 0.50000000000000 0.00000000000000 0.33333333333333 15.999 > 1 12 O 0.83333333333333 0.00000000000000 0.33333333333333 15.999 > 1 13 O 0.16666666666667 0.33333333333333 0.33333333333333 15.999 > 1 14 O 0.50000000000000 0.33333333333333 0.33333333333333 15.999 > 1 15 O 0.83333333333333 0.33333333333333 0.33333333333333 15.999 > 1 16 O 0.16666666666667 0.66666666666667 0.33333333333333 15.999 > 1 17 O 0.50000000000000 0.66666666666667 0.33333333333333 15.999 > 1 18 O 0.83333333333333 0.66666666666667 0.33333333333333 15.999 > 1 19 O 0.16666666666667 0.00000000000000 0.66666666666667 15.999 > 1 20 O 0.50000000000000 0.00000000000000 0.66666666666667 15.999 > 1 21 O 0.83333333333333 0.00000000000000 0.66666666666667 15.999 > 1 22 O 0.16666666666667 0.33333333333333 0.66666666666667 15.999 > 1 23 O 0.50000000000000 0.33333333333333 0.66666666666667 15.999 > 1 24 O 0.83333333333333 0.33333333333333 0.66666666666667 15.999 > 1 25 O 0.16666666666667 0.66666666666667 0.66666666666667 15.999 > 1 26 O 0.50000000000000 0.66666666666667 0.66666666666667 15.999 > 1 27 O 0.83333333333333 0.66666666666667 0.66666666666667 15.999 > 1 28 O 0.00000000000000 0.00000000000000 0.16666666666667 15.999 > 2 29 O 0.33333333333333 0.00000000000000 0.16666666666667 15.999 > 2 30 O 0.66666666666667 0.00000000000000 0.16666666666667 15.999 > 2 31 O 0.00000000000000 0.33333333333333 0.16666666666667 15.999 > 2 32 O 0.33333333333333 0.33333333333333 0.16666666666667 15.999 > 2 33 O 0.66666666666667 0.33333333333333 0.16666666666667 15.999 > 2 34 O 0.00000000000000 0.66666666666667 0.16666666666667 15.999 > 2 35 O 0.33333333333333 0.66666666666667 0.16666666666667 15.999 > 2 36 O 0.66666666666667 0.66666666666667 0.16666666666667 15.999 > 2 37 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 2 38 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 2 39 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 2 40 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 2 41 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 2 42 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 2 43 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 2 44 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 2 45 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 2 46 O 0.00000000000000 0.00000000000000 0.83333333333333 15.999 > 2 47 O 0.33333333333333 0.00000000000000 0.83333333333333 15.999 > 2 48 O 0.66666666666667 0.00000000000000 0.83333333333333 15.999 > 2 49 O 0.00000000000000 0.33333333333333 0.83333333333333 15.999 > 2 50 O 0.33333333333333 0.33333333333333 0.83333333333333 15.999 > 2 51 O 0.66666666666667 0.33333333333333 0.83333333333333 15.999 > 2 52 O 0.00000000000000 0.66666666666667 0.83333333333333 15.999 > 2 53 O 0.33333333333333 0.66666666666667 0.83333333333333 15.999 > 2 54 O 0.66666666666667 0.66666666666667 0.83333333333333 15.999 > 2 55 O 0.00000000000000 0.16666666666667 0.00000000000000 15.999 > 3 56 O 0.33333333333333 0.16666666666667 0.00000000000000 15.999 > 3 57 O 0.66666666666667 0.16666666666667 0.00000000000000 15.999 > 3 58 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 > 3 59 O 0.33333333333333 0.50000000000000 0.00000000000000 15.999 > 3 60 O 0.66666666666667 0.50000000000000 0.00000000000000 15.999 > 3 61 O 0.00000000000000 0.83333333333333 0.00000000000000 15.999 > 3 62 O 0.33333333333333 0.83333333333333 0.00000000000000 15.999 > 3 63 O 0.66666666666667 0.83333333333333 0.00000000000000 15.999 > 3 64 O 0.00000000000000 0.16666666666667 0.33333333333333 15.999 > 3 65 O 0.33333333333333 0.16666666666667 0.33333333333333 15.999 > 3 66 O 0.66666666666667 0.16666666666667 0.33333333333333 15.999 > 3 67 O 0.00000000000000 0.50000000000000 0.33333333333333 15.999 > 3 68 O 0.33333333333333 0.50000000000000 0.33333333333333 15.999 > 3 69 O 0.66666666666667 0.50000000000000 0.33333333333333 15.999 > 3 70 O 0.00000000000000 0.83333333333333 0.33333333333333 15.999 > 3 71 O 0.33333333333333 0.83333333333333 0.33333333333333 15.999 > 3 72 O 0.66666666666667 0.83333333333333 0.33333333333333 15.999 > 3 73 O 0.00000000000000 0.16666666666667 0.66666666666667 15.999 > 3 74 O 0.33333333333333 0.16666666666667 0.66666666666667 15.999 > 3 75 O 0.66666666666667 0.16666666666667 0.66666666666667 15.999 > 3 76 O 0.00000000000000 0.50000000000000 0.66666666666667 15.999 > 3 77 O 0.33333333333333 0.50000000000000 0.66666666666667 15.999 > 3 78 O 0.66666666666667 0.50000000000000 0.66666666666667 15.999 > 3 79 O 0.00000000000000 0.83333333333333 0.66666666666667 15.999 > 3 80 O 0.33333333333333 0.83333333333333 0.66666666666667 15.999 > 3 81 O 0.66666666666667 0.83333333333333 0.66666666666667 15.999 > 3 *82 Re 0.00000000000000 0.00000000000000 0.00000000000000 186.207 > 4 83 Re 0.33333333333333 0.00000000000000 0.00000000000000 186.207 > 4 84 Re 0.66666666666667 0.00000000000000 0.00000000000000 186.207 > 4 85 Re 0.00000000000000 0.33333333333333 0.00000000000000 186.207 > 4 86 Re 0.33333333333333 0.33333333333333 0.00000000000000 186.207 > 4 87 Re 0.66666666666667 0.33333333333333 0.00000000000000 186.207 > 4 88 Re 0.00000000000000 0.66666666666667 0.00000000000000 186.207 > 4 89 Re 0.33333333333333 0.66666666666667 0.00000000000000 186.207 > 4 90 Re 0.66666666666667 0.66666666666667 0.00000000000000 186.207 > 4 91 Re 0.00000000000000 0.00000000000000 0.33333333333333 186.207 > 4 92 Re 0.33333333333333 0.00000000000000 0.33333333333333 186.207 > 4 93 Re 0.66666666666667 0.00000000000000 0.33333333333333 186.207 > 4 94 Re 0.00000000000000 0.33333333333333 0.33333333333333 186.207 > 4 95 Re 0.33333333333333 0.33333333333333 0.33333333333333 186.207 > 4 96 Re 0.66666666666667 0.33333333333333 0.33333333333333 186.207 > 4 97 Re 0.00000000000000 0.66666666666667 0.33333333333333 186.207 > 4 98 Re 0.33333333333333 0.66666666666667 0.33333333333333 186.207 > 4 99 Re 0.66666666666667 0.66666666666667 0.33333333333333 186.207 > 4 100 Re 0.00000000000000 0.00000000000000 0.66666666666667 186.207 > 4 101 Re 0.33333333333333 0.00000000000000 0.66666666666667 186.207 > 4 102 Re 0.66666666666667 0.00000000000000 0.66666666666667 186.207 > 4 103 Re 0.00000000000000 0.33333333333333 0.66666666666667 186.207 > 4 104 Re 0.33333333333333 0.33333333333333 0.66666666666667 186.207 > 4 105 Re 0.66666666666667 0.33333333333333 0.66666666666667 186.207 > 4 106 Re 0.00000000000000 0.66666666666667 0.66666666666667 186.207 > 4 107 Re 0.33333333333333 0.66666666666667 0.66666666666667 186.207 > 4 108 Re 0.66666666666667 0.66666666666667 0.66666666666667 186.207 > 4 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 1443.6115334 0.0000000 0.0000000 0.0000000 1443.6115334 0.0000000 0.0000000 0.0000000 1443.6115334 -------------------------- Born effective charges -------------------------- 1 O -0.4222211 0.0000000 0.0000000 0.0000000 -3.6068373 0.0000000 0.0000000 0.0000000 -3.6068373 2 O -3.6068373 0.0000000 0.0000000 0.0000000 -3.6068373 0.0000000 0.0000000 0.0000000 -0.4222211 3 O -3.6068373 0.0000000 0.0000000 0.0000000 -0.4222211 0.0000000 0.0000000 0.0000000 -3.6068373 4 Re 7.6358957 0.0000000 0.0000000 0.0000000 7.6358957 0.0000000 0.0000000 0.0000000 7.6358957 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 324/324 Permutation basis: 2406/2406 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 68 Number of blocks in projector: 68 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: 48 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 20 Use standard eigh solver. Tree of FC basis block matrices: - (68, 65), data: False |-- (20, 20), data: True |-- (48, 45), data: True ----- Solver_atoms: 1 -- 108 / 108 Time (Solver_compr_matrix_reshape): 0.003 Solver_block: 80 / 80 - Time: 0.020 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.023 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 324/324 Permutation basis: 2406/2406 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 68 Number of blocks in projector: 68 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: 48 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 20 Use standard eigh solver. Tree of FC basis block matrices: - (68, 65), data: False |-- (20, 20), data: True |-- (48, 45), data: True Max drift after symmetrization by symfc projector: 0.00000000 (xx) 0.00000000 (xx) 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:23]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ 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): [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.765565529999999 0.000000000000000 0.000000000000000 b 0.000000000000000 3.765565529999999 0.000000000000000 c 0.000000000000000 0.000000000000000 3.765565529999999 Atomic positions (fractional): 1 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 2 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 3 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 4 Re 0.00000000000000 0.00000000000000 0.00000000000000 186.207 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.296696589999996 0.000000000000000 0.000000000000000 b 0.000000000000000 11.296696589999996 0.000000000000000 c 0.000000000000000 0.000000000000000 11.296696589999996 Atomic positions (fractional): 1 O 0.16666666666667 0.00000000000000 0.00000000000000 15.999 > 1 2 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 > 1 3 O 0.83333333333333 0.00000000000000 0.00000000000000 15.999 > 1 4 O 0.16666666666667 0.33333333333333 0.00000000000000 15.999 > 1 5 O 0.50000000000000 0.33333333333333 0.00000000000000 15.999 > 1 6 O 0.83333333333333 0.33333333333333 0.00000000000000 15.999 > 1 7 O 0.16666666666667 0.66666666666667 0.00000000000000 15.999 > 1 8 O 0.50000000000000 0.66666666666667 0.00000000000000 15.999 > 1 9 O 0.83333333333333 0.66666666666667 0.00000000000000 15.999 > 1 10 O 0.16666666666667 0.00000000000000 0.33333333333333 15.999 > 1 11 O 0.50000000000000 0.00000000000000 0.33333333333333 15.999 > 1 12 O 0.83333333333333 0.00000000000000 0.33333333333333 15.999 > 1 13 O 0.16666666666667 0.33333333333333 0.33333333333333 15.999 > 1 14 O 0.50000000000000 0.33333333333333 0.33333333333333 15.999 > 1 15 O 0.83333333333333 0.33333333333333 0.33333333333333 15.999 > 1 16 O 0.16666666666667 0.66666666666667 0.33333333333333 15.999 > 1 17 O 0.50000000000000 0.66666666666667 0.33333333333333 15.999 > 1 18 O 0.83333333333333 0.66666666666667 0.33333333333333 15.999 > 1 19 O 0.16666666666667 0.00000000000000 0.66666666666667 15.999 > 1 20 O 0.50000000000000 0.00000000000000 0.66666666666667 15.999 > 1 21 O 0.83333333333333 0.00000000000000 0.66666666666667 15.999 > 1 22 O 0.16666666666667 0.33333333333333 0.66666666666667 15.999 > 1 23 O 0.50000000000000 0.33333333333333 0.66666666666667 15.999 > 1 24 O 0.83333333333333 0.33333333333333 0.66666666666667 15.999 > 1 25 O 0.16666666666667 0.66666666666667 0.66666666666667 15.999 > 1 26 O 0.50000000000000 0.66666666666667 0.66666666666667 15.999 > 1 27 O 0.83333333333333 0.66666666666667 0.66666666666667 15.999 > 1 28 O 0.00000000000000 0.00000000000000 0.16666666666667 15.999 > 28 29 O 0.33333333333333 0.00000000000000 0.16666666666667 15.999 > 28 30 O 0.66666666666667 0.00000000000000 0.16666666666667 15.999 > 28 31 O 0.00000000000000 0.33333333333333 0.16666666666667 15.999 > 28 32 O 0.33333333333333 0.33333333333333 0.16666666666667 15.999 > 28 33 O 0.66666666666667 0.33333333333333 0.16666666666667 15.999 > 28 34 O 0.00000000000000 0.66666666666667 0.16666666666667 15.999 > 28 35 O 0.33333333333333 0.66666666666667 0.16666666666667 15.999 > 28 36 O 0.66666666666667 0.66666666666667 0.16666666666667 15.999 > 28 37 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 28 38 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 28 39 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 28 40 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 28 41 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 28 42 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 28 43 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 28 44 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 28 45 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 28 46 O 0.00000000000000 0.00000000000000 0.83333333333333 15.999 > 28 47 O 0.33333333333333 0.00000000000000 0.83333333333333 15.999 > 28 48 O 0.66666666666667 0.00000000000000 0.83333333333333 15.999 > 28 49 O 0.00000000000000 0.33333333333333 0.83333333333333 15.999 > 28 50 O 0.33333333333333 0.33333333333333 0.83333333333333 15.999 > 28 51 O 0.66666666666667 0.33333333333333 0.83333333333333 15.999 > 28 52 O 0.00000000000000 0.66666666666667 0.83333333333333 15.999 > 28 53 O 0.33333333333333 0.66666666666667 0.83333333333333 15.999 > 28 54 O 0.66666666666667 0.66666666666667 0.83333333333333 15.999 > 28 55 O 0.00000000000000 0.16666666666667 0.00000000000000 15.999 > 55 56 O 0.33333333333333 0.16666666666667 0.00000000000000 15.999 > 55 57 O 0.66666666666667 0.16666666666667 0.00000000000000 15.999 > 55 58 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 > 55 59 O 0.33333333333333 0.50000000000000 0.00000000000000 15.999 > 55 60 O 0.66666666666667 0.50000000000000 0.00000000000000 15.999 > 55 61 O 0.00000000000000 0.83333333333333 0.00000000000000 15.999 > 55 62 O 0.33333333333333 0.83333333333333 0.00000000000000 15.999 > 55 63 O 0.66666666666667 0.83333333333333 0.00000000000000 15.999 > 55 64 O 0.00000000000000 0.16666666666667 0.33333333333333 15.999 > 55 65 O 0.33333333333333 0.16666666666667 0.33333333333333 15.999 > 55 66 O 0.66666666666667 0.16666666666667 0.33333333333333 15.999 > 55 67 O 0.00000000000000 0.50000000000000 0.33333333333333 15.999 > 55 68 O 0.33333333333333 0.50000000000000 0.33333333333333 15.999 > 55 69 O 0.66666666666667 0.50000000000000 0.33333333333333 15.999 > 55 70 O 0.00000000000000 0.83333333333333 0.33333333333333 15.999 > 55 71 O 0.33333333333333 0.83333333333333 0.33333333333333 15.999 > 55 72 O 0.66666666666667 0.83333333333333 0.33333333333333 15.999 > 55 73 O 0.00000000000000 0.16666666666667 0.66666666666667 15.999 > 55 74 O 0.33333333333333 0.16666666666667 0.66666666666667 15.999 > 55 75 O 0.66666666666667 0.16666666666667 0.66666666666667 15.999 > 55 76 O 0.00000000000000 0.50000000000000 0.66666666666667 15.999 > 55 77 O 0.33333333333333 0.50000000000000 0.66666666666667 15.999 > 55 78 O 0.66666666666667 0.50000000000000 0.66666666666667 15.999 > 55 79 O 0.00000000000000 0.83333333333333 0.66666666666667 15.999 > 55 80 O 0.33333333333333 0.83333333333333 0.66666666666667 15.999 > 55 81 O 0.66666666666667 0.83333333333333 0.66666666666667 15.999 > 55 82 Re 0.00000000000000 0.00000000000000 0.00000000000000 186.207 > 82 83 Re 0.33333333333333 0.00000000000000 0.00000000000000 186.207 > 82 84 Re 0.66666666666667 0.00000000000000 0.00000000000000 186.207 > 82 85 Re 0.00000000000000 0.33333333333333 0.00000000000000 186.207 > 82 86 Re 0.33333333333333 0.33333333333333 0.00000000000000 186.207 > 82 87 Re 0.66666666666667 0.33333333333333 0.00000000000000 186.207 > 82 88 Re 0.00000000000000 0.66666666666667 0.00000000000000 186.207 > 82 89 Re 0.33333333333333 0.66666666666667 0.00000000000000 186.207 > 82 90 Re 0.66666666666667 0.66666666666667 0.00000000000000 186.207 > 82 91 Re 0.00000000000000 0.00000000000000 0.33333333333333 186.207 > 82 92 Re 0.33333333333333 0.00000000000000 0.33333333333333 186.207 > 82 93 Re 0.66666666666667 0.00000000000000 0.33333333333333 186.207 > 82 94 Re 0.00000000000000 0.33333333333333 0.33333333333333 186.207 > 82 95 Re 0.33333333333333 0.33333333333333 0.33333333333333 186.207 > 82 96 Re 0.66666666666667 0.33333333333333 0.33333333333333 186.207 > 82 97 Re 0.00000000000000 0.66666666666667 0.33333333333333 186.207 > 82 98 Re 0.33333333333333 0.66666666666667 0.33333333333333 186.207 > 82 99 Re 0.66666666666667 0.66666666666667 0.33333333333333 186.207 > 82 100 Re 0.00000000000000 0.00000000000000 0.66666666666667 186.207 > 82 101 Re 0.33333333333333 0.00000000000000 0.66666666666667 186.207 > 82 102 Re 0.66666666666667 0.00000000000000 0.66666666666667 186.207 > 82 103 Re 0.00000000000000 0.33333333333333 0.66666666666667 186.207 > 82 104 Re 0.33333333333333 0.33333333333333 0.66666666666667 186.207 > 82 105 Re 0.66666666666667 0.33333333333333 0.66666666666667 186.207 > 82 106 Re 0.00000000000000 0.66666666666667 0.66666666666667 186.207 > 82 107 Re 0.33333333333333 0.66666666666667 0.66666666666667 186.207 > 82 108 Re 0.66666666666667 0.66666666666667 0.66666666666667 186.207 > 82 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 1443.6115334 0.0000000 0.0000000 0.0000000 1443.6115334 0.0000000 0.0000000 0.0000000 1443.6115334 -------------------------- Born effective charges -------------------------- 1 O -0.4222211 0.0000000 0.0000000 0.0000000 -3.6068373 0.0000000 0.0000000 0.0000000 -3.6068373 2 O -3.6068373 0.0000000 0.0000000 0.0000000 -3.6068373 0.0000000 0.0000000 0.0000000 -0.4222211 3 O -3.6068373 0.0000000 0.0000000 0.0000000 -0.4222211 0.0000000 0.0000000 0.0000000 -3.6068373 4 Re 7.6358957 0.0000000 0.0000000 0.0000000 7.6358957 0.0000000 0.0000000 0.0000000 7.6358957 ---------------------------------------------------------------------------- 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.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 82, 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.00000000 (xxx) 0.00000000 (xxx) 0.00000000 (xxx) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): 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:25]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:30:26]------------------------- 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): [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.765565529999999 0.000000000000000 0.000000000000000 b 0.000000000000000 3.765565529999999 0.000000000000000 c 0.000000000000000 0.000000000000000 3.765565529999999 Atomic positions (fractional): 1 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 2 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 3 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 4 Re 0.00000000000000 0.00000000000000 0.00000000000000 186.207 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.296696589999996 0.000000000000000 0.000000000000000 b 0.000000000000000 11.296696589999996 0.000000000000000 c 0.000000000000000 0.000000000000000 11.296696589999996 Atomic positions (fractional): 1 O 0.16666666666667 0.00000000000000 0.00000000000000 15.999 > 1 2 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 > 1 3 O 0.83333333333333 0.00000000000000 0.00000000000000 15.999 > 1 4 O 0.16666666666667 0.33333333333333 0.00000000000000 15.999 > 1 5 O 0.50000000000000 0.33333333333333 0.00000000000000 15.999 > 1 6 O 0.83333333333333 0.33333333333333 0.00000000000000 15.999 > 1 7 O 0.16666666666667 0.66666666666667 0.00000000000000 15.999 > 1 8 O 0.50000000000000 0.66666666666667 0.00000000000000 15.999 > 1 9 O 0.83333333333333 0.66666666666667 0.00000000000000 15.999 > 1 10 O 0.16666666666667 0.00000000000000 0.33333333333333 15.999 > 1 11 O 0.50000000000000 0.00000000000000 0.33333333333333 15.999 > 1 12 O 0.83333333333333 0.00000000000000 0.33333333333333 15.999 > 1 13 O 0.16666666666667 0.33333333333333 0.33333333333333 15.999 > 1 14 O 0.50000000000000 0.33333333333333 0.33333333333333 15.999 > 1 15 O 0.83333333333333 0.33333333333333 0.33333333333333 15.999 > 1 16 O 0.16666666666667 0.66666666666667 0.33333333333333 15.999 > 1 17 O 0.50000000000000 0.66666666666667 0.33333333333333 15.999 > 1 18 O 0.83333333333333 0.66666666666667 0.33333333333333 15.999 > 1 19 O 0.16666666666667 0.00000000000000 0.66666666666667 15.999 > 1 20 O 0.50000000000000 0.00000000000000 0.66666666666667 15.999 > 1 21 O 0.83333333333333 0.00000000000000 0.66666666666667 15.999 > 1 22 O 0.16666666666667 0.33333333333333 0.66666666666667 15.999 > 1 23 O 0.50000000000000 0.33333333333333 0.66666666666667 15.999 > 1 24 O 0.83333333333333 0.33333333333333 0.66666666666667 15.999 > 1 25 O 0.16666666666667 0.66666666666667 0.66666666666667 15.999 > 1 26 O 0.50000000000000 0.66666666666667 0.66666666666667 15.999 > 1 27 O 0.83333333333333 0.66666666666667 0.66666666666667 15.999 > 1 28 O 0.00000000000000 0.00000000000000 0.16666666666667 15.999 > 28 29 O 0.33333333333333 0.00000000000000 0.16666666666667 15.999 > 28 30 O 0.66666666666667 0.00000000000000 0.16666666666667 15.999 > 28 31 O 0.00000000000000 0.33333333333333 0.16666666666667 15.999 > 28 32 O 0.33333333333333 0.33333333333333 0.16666666666667 15.999 > 28 33 O 0.66666666666667 0.33333333333333 0.16666666666667 15.999 > 28 34 O 0.00000000000000 0.66666666666667 0.16666666666667 15.999 > 28 35 O 0.33333333333333 0.66666666666667 0.16666666666667 15.999 > 28 36 O 0.66666666666667 0.66666666666667 0.16666666666667 15.999 > 28 37 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 28 38 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 28 39 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 28 40 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 28 41 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 28 42 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 28 43 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 28 44 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 28 45 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 28 46 O 0.00000000000000 0.00000000000000 0.83333333333333 15.999 > 28 47 O 0.33333333333333 0.00000000000000 0.83333333333333 15.999 > 28 48 O 0.66666666666667 0.00000000000000 0.83333333333333 15.999 > 28 49 O 0.00000000000000 0.33333333333333 0.83333333333333 15.999 > 28 50 O 0.33333333333333 0.33333333333333 0.83333333333333 15.999 > 28 51 O 0.66666666666667 0.33333333333333 0.83333333333333 15.999 > 28 52 O 0.00000000000000 0.66666666666667 0.83333333333333 15.999 > 28 53 O 0.33333333333333 0.66666666666667 0.83333333333333 15.999 > 28 54 O 0.66666666666667 0.66666666666667 0.83333333333333 15.999 > 28 55 O 0.00000000000000 0.16666666666667 0.00000000000000 15.999 > 55 56 O 0.33333333333333 0.16666666666667 0.00000000000000 15.999 > 55 57 O 0.66666666666667 0.16666666666667 0.00000000000000 15.999 > 55 58 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 > 55 59 O 0.33333333333333 0.50000000000000 0.00000000000000 15.999 > 55 60 O 0.66666666666667 0.50000000000000 0.00000000000000 15.999 > 55 61 O 0.00000000000000 0.83333333333333 0.00000000000000 15.999 > 55 62 O 0.33333333333333 0.83333333333333 0.00000000000000 15.999 > 55 63 O 0.66666666666667 0.83333333333333 0.00000000000000 15.999 > 55 64 O 0.00000000000000 0.16666666666667 0.33333333333333 15.999 > 55 65 O 0.33333333333333 0.16666666666667 0.33333333333333 15.999 > 55 66 O 0.66666666666667 0.16666666666667 0.33333333333333 15.999 > 55 67 O 0.00000000000000 0.50000000000000 0.33333333333333 15.999 > 55 68 O 0.33333333333333 0.50000000000000 0.33333333333333 15.999 > 55 69 O 0.66666666666667 0.50000000000000 0.33333333333333 15.999 > 55 70 O 0.00000000000000 0.83333333333333 0.33333333333333 15.999 > 55 71 O 0.33333333333333 0.83333333333333 0.33333333333333 15.999 > 55 72 O 0.66666666666667 0.83333333333333 0.33333333333333 15.999 > 55 73 O 0.00000000000000 0.16666666666667 0.66666666666667 15.999 > 55 74 O 0.33333333333333 0.16666666666667 0.66666666666667 15.999 > 55 75 O 0.66666666666667 0.16666666666667 0.66666666666667 15.999 > 55 76 O 0.00000000000000 0.50000000000000 0.66666666666667 15.999 > 55 77 O 0.33333333333333 0.50000000000000 0.66666666666667 15.999 > 55 78 O 0.66666666666667 0.50000000000000 0.66666666666667 15.999 > 55 79 O 0.00000000000000 0.83333333333333 0.66666666666667 15.999 > 55 80 O 0.33333333333333 0.83333333333333 0.66666666666667 15.999 > 55 81 O 0.66666666666667 0.83333333333333 0.66666666666667 15.999 > 55 82 Re 0.00000000000000 0.00000000000000 0.00000000000000 186.207 > 82 83 Re 0.33333333333333 0.00000000000000 0.00000000000000 186.207 > 82 84 Re 0.66666666666667 0.00000000000000 0.00000000000000 186.207 > 82 85 Re 0.00000000000000 0.33333333333333 0.00000000000000 186.207 > 82 86 Re 0.33333333333333 0.33333333333333 0.00000000000000 186.207 > 82 87 Re 0.66666666666667 0.33333333333333 0.00000000000000 186.207 > 82 88 Re 0.00000000000000 0.66666666666667 0.00000000000000 186.207 > 82 89 Re 0.33333333333333 0.66666666666667 0.00000000000000 186.207 > 82 90 Re 0.66666666666667 0.66666666666667 0.00000000000000 186.207 > 82 91 Re 0.00000000000000 0.00000000000000 0.33333333333333 186.207 > 82 92 Re 0.33333333333333 0.00000000000000 0.33333333333333 186.207 > 82 93 Re 0.66666666666667 0.00000000000000 0.33333333333333 186.207 > 82 94 Re 0.00000000000000 0.33333333333333 0.33333333333333 186.207 > 82 95 Re 0.33333333333333 0.33333333333333 0.33333333333333 186.207 > 82 96 Re 0.66666666666667 0.33333333333333 0.33333333333333 186.207 > 82 97 Re 0.00000000000000 0.66666666666667 0.33333333333333 186.207 > 82 98 Re 0.33333333333333 0.66666666666667 0.33333333333333 186.207 > 82 99 Re 0.66666666666667 0.66666666666667 0.33333333333333 186.207 > 82 100 Re 0.00000000000000 0.00000000000000 0.66666666666667 186.207 > 82 101 Re 0.33333333333333 0.00000000000000 0.66666666666667 186.207 > 82 102 Re 0.66666666666667 0.00000000000000 0.66666666666667 186.207 > 82 103 Re 0.00000000000000 0.33333333333333 0.66666666666667 186.207 > 82 104 Re 0.33333333333333 0.33333333333333 0.66666666666667 186.207 > 82 105 Re 0.66666666666667 0.33333333333333 0.66666666666667 186.207 > 82 106 Re 0.00000000000000 0.66666666666667 0.66666666666667 186.207 > 82 107 Re 0.33333333333333 0.66666666666667 0.66666666666667 186.207 > 82 108 Re 0.66666666666667 0.66666666666667 0.66666666666667 186.207 > 82 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 1443.6115334 0.0000000 0.0000000 0.0000000 1443.6115334 0.0000000 0.0000000 0.0000000 1443.6115334 -------------------------- Born effective charges -------------------------- 1 O -0.4222211 0.0000000 0.0000000 0.0000000 -3.6068373 0.0000000 0.0000000 0.0000000 -3.6068373 2 O -3.6068373 0.0000000 0.0000000 0.0000000 -3.6068373 0.0000000 0.0000000 0.0000000 -0.4222211 3 O -3.6068373 0.0000000 0.0000000 0.0000000 -0.4222211 0.0000000 0.0000000 0.0000000 -3.6068373 4 Re 7.6358957 0.0000000 0.0000000 0.0000000 7.6358957 0.0000000 0.0000000 0.0000000 7.6358957 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000000 (xxx) 0.00000000 (xxx) 0.00000000 (xxx) Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) --------------------------- 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 ... [ 13 13 13 ] fc3-r2q-transformation over three atoms: True --------------------------- Phonon calculations ---------------------------- Use NAC by Gonze et al. (no real space sum in current implementation) PRB 50, 13035(R) (1994), PRB 55, 10355 (1997) G-cutoff distance: 1.10, Number of G-points: 305, Lambda: 4.37 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/84) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 84 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 10.756 ( 0.000 0.000 0.000) 0.000 10.756 ( 0.000 0.000 0.000) 0.000 10.756 ( 0.000 0.000 0.000) 0.000 11.344 ( 0.000 0.000 0.000) 0.000 11.344 ( 0.000 0.000 0.000) 0.000 11.344 ( 0.000 0.000 0.000) 0.000 21.228 ( 0.000 0.000 0.000) 0.000 21.228 ( 0.000 0.000 0.000) 0.000 21.228 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/84) ======================= q-point: ( 0.08 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.661 ( 30.765 0.000 0.000) 30.765 0.661 ( 30.765 0.000 0.000) 30.765 2.179 ( 100.762 0.000 0.000) 100.762 10.717 ( -3.815 0.000 0.000) 3.815 10.717 ( -3.815 0.000 0.000) 3.815 10.838 ( 3.895 0.000 0.000) 3.895 11.313 ( -2.695 0.000 0.000) 2.695 11.313 ( -2.695 0.000 0.000) 2.695 11.336 ( -0.783 0.000 0.000) 0.783 21.191 ( -3.574 0.000 0.000) 3.574 21.191 ( -3.574 0.000 0.000) 3.574 21.290 ( 5.848 0.000 0.000) 5.848 ======================= Grid point 2 (3/84) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.227 ( 23.955 0.000 0.000) 23.955 1.227 ( 23.955 0.000 0.000) 23.955 4.017 ( 77.437 0.000 0.000) 77.437 10.601 ( -7.285 0.000 0.000) 7.285 10.601 ( -7.285 0.000 0.000) 7.285 11.005 ( 13.791 0.000 0.000) 13.791 11.251 ( -3.078 0.000 0.000) 3.078 11.251 ( -3.078 0.000 0.000) 3.078 11.313 ( -1.390 0.000 0.000) 1.390 21.084 ( -6.736 0.000 0.000) 6.736 21.084 ( -6.736 0.000 0.000) 6.736 21.460 ( 10.307 0.000 0.000) 10.307 ======================= Grid point 3 (4/84) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.620 ( 14.241 0.000 0.000) 14.241 1.620 ( 14.241 0.000 0.000) 14.241 5.307 ( 49.015 0.000 0.000) 49.015 10.442 ( -7.663 0.000 0.000) 7.663 10.442 ( -7.663 0.000 0.000) 7.663 11.194 ( -2.503 0.000 0.000) 2.503 11.194 ( -2.503 0.000 0.000) 2.503 11.281 ( -1.681 0.000 0.000) 1.681 11.458 ( 31.483 0.000 0.000) 31.483 20.923 ( -8.883 0.000 0.000) 8.883 20.923 ( -8.883 0.000 0.000) 8.883 21.674 ( 9.423 0.000 0.000) 9.423 ======================= Grid point 4 (5/84) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.809 ( 4.646 0.000 0.000) 4.646 1.809 ( 4.646 0.000 0.000) 4.646 6.062 ( 26.444 0.000 0.000) 26.444 10.314 ( -4.435 0.000 0.000) 4.435 10.314 ( -4.435 0.000 0.000) 4.435 11.145 ( -2.391 0.000 0.000) 2.391 11.145 ( -2.391 0.000 0.000) 2.391 11.247 ( -1.588 0.000 0.000) 1.588 12.279 ( 47.282 0.000 0.000) 47.282 20.734 ( -9.216 0.000 0.000) 9.216 20.734 ( -9.216 0.000 0.000) 9.216 21.796 ( 1.664 0.000 0.000) 1.664 ======================= Grid point 5 (6/84) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.836 ( -1.087 0.000 0.000) 1.087 1.836 ( -1.087 0.000 0.000) 1.087 6.450 ( 12.763 0.000 0.000) 12.763 10.263 ( -0.758 0.000 0.000) 0.758 10.263 ( -0.758 0.000 0.000) 0.758 11.096 ( -2.270 0.000 0.000) 2.270 11.096 ( -2.270 0.000 0.000) 2.270 11.219 ( -1.129 0.000 0.000) 1.129 13.274 ( 46.235 0.000 0.000) 46.235 20.563 ( -7.097 0.000 0.000) 7.097 20.563 ( -7.097 0.000 0.000) 7.097 21.737 ( -6.664 0.000 0.000) 6.664 ======================= Grid point 6 (7/84) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.804 ( -1.273 0.000 0.000) 1.273 1.804 ( -1.273 0.000 0.000) 1.273 6.615 ( 3.851 0.000 0.000) 3.851 10.265 ( 0.352 0.000 0.000) 0.352 10.265 ( 0.352 0.000 0.000) 0.352 11.060 ( -1.053 0.000 0.000) 1.053 11.060 ( -1.053 0.000 0.000) 1.053 11.203 ( -0.408 0.000 0.000) 0.408 13.993 ( 20.104 0.000 0.000) 20.104 20.460 ( -2.686 0.000 0.000) 2.686 20.460 ( -2.686 0.000 0.000) 2.686 21.595 ( -4.996 0.000 0.000) 4.996 ======================= Grid point 14 (8/84) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.935 ( 21.741 21.741 0.000) 30.747 2.232 ( 51.842 51.842 0.000) 73.316 2.294 ( 51.551 51.551 0.000) 72.904 10.486 ( -13.159 -13.159 0.000) 18.610 10.683 ( -3.394 -3.394 0.000) 4.800 11.047 ( 11.373 11.373 0.000) 16.084 11.181 ( -7.077 -7.077 0.000) 10.008 11.279 ( -3.052 -3.052 0.000) 4.317 11.451 ( 5.666 5.666 0.000) 8.012 21.154 ( -3.611 -3.611 0.000) 5.107 21.160 ( -3.332 -3.332 0.000) 4.712 21.347 ( 5.884 5.884 0.000) 8.321 ======================= Grid point 15 (9/84) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.393 ( 21.060 14.541 0.000) 25.592 2.459 ( 10.428 85.192 0.000) 85.828 4.042 ( 75.618 2.291 0.000) 75.652 10.228 ( -11.667 -24.178 0.000) 26.846 10.578 ( -6.704 -2.340 0.000) 7.101 11.046 ( -5.672 -11.067 0.000) 12.436 11.207 ( -3.542 -3.927 0.000) 5.288 11.288 ( 12.308 14.374 0.000) 18.923 11.569 ( 6.771 17.582 0.000) 18.841 21.046 ( -6.803 -3.715 0.000) 7.751 21.055 ( -6.728 -2.965 0.000) 7.353 21.526 ( 11.150 6.645 0.000) 12.980 ======================= Grid point 16 (10/84) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.748 ( 13.133 11.486 0.000) 17.447 2.651 ( 7.737 76.544 0.000) 76.934 5.305 ( 48.144 -0.129 0.000) 48.144 10.020 ( -8.682 -30.754 0.000) 31.956 10.429 ( -7.200 -1.261 0.000) 7.310 10.964 ( -2.171 -13.315 0.000) 13.491 11.143 ( -2.814 -4.734 0.000) 5.507 11.491 ( 5.732 9.041 0.000) 10.705 11.840 ( 22.997 29.643 0.000) 37.517 20.883 ( -8.968 -3.868 0.000) 9.766 20.892 ( -9.109 -2.936 0.000) 9.571 21.768 ( 11.381 9.200 0.000) 14.635 ======================= Grid point 17 (11/84) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.923 ( 4.290 10.282 0.000) 11.141 2.760 ( 2.846 72.564 0.000) 72.620 6.047 ( 25.964 -1.415 0.000) 26.002 9.870 ( -6.079 -34.583 0.000) 35.113 10.310 ( -4.066 -0.427 0.000) 4.088 10.953 ( 0.711 -11.147 0.000) 11.169 11.088 ( -2.606 -5.254 0.000) 5.865 11.554 ( 1.382 15.689 0.000) 15.749 12.498 ( 39.739 19.677 0.000) 44.344 20.692 ( -9.709 -3.800 0.000) 10.426 20.692 ( -9.717 -3.795 0.000) 10.432 21.945 ( 5.121 14.094 0.000) 14.996 ======================= Grid point 18 (12/84) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.947 ( -1.096 9.984 0.000) 10.044 2.778 ( -0.569 72.302 0.000) 72.304 6.426 ( 12.383 -2.382 0.000) 12.610 9.771 ( -3.556 -38.602 0.000) 38.766 10.267 ( -0.436 0.278 0.000) 0.517 10.978 ( 1.470 -6.989 0.000) 7.142 11.035 ( -2.475 -5.678 0.000) 6.194 11.566 ( 0.085 21.017 0.000) 21.017 13.358 ( 40.905 7.865 0.000) 41.655 20.492 ( -9.084 -5.468 0.000) 10.603 20.520 ( -7.161 -4.193 0.000) 8.298 21.970 ( -2.098 21.138 0.000) 21.241 ======================= Grid point 19 (13/84) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.916 ( -1.227 10.036 0.000) 10.111 2.759 ( -0.736 73.207 0.000) 73.210 6.586 ( 3.686 -3.066 0.000) 4.794 9.724 ( -1.133 -41.712 0.000) 41.727 10.273 ( 0.509 0.806 0.000) 0.953 10.995 ( -1.166 -6.038 0.000) 6.150 11.002 ( 0.694 -4.464 0.000) 4.518 11.564 ( -0.129 23.968 0.000) 23.969 14.002 ( 18.154 0.856 0.000) 18.175 20.348 ( -4.106 -7.856 0.000) 8.864 20.416 ( -2.710 -4.285 0.000) 5.070 21.909 ( -2.412 27.331 0.000) 27.437 ======================= Grid point 28 (14/84) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.732 ( 16.801 16.801 0.000) 23.761 4.079 ( 35.211 35.211 0.000) 49.797 4.130 ( 39.909 39.909 0.000) 56.440 9.758 ( -20.929 -20.929 0.000) 29.598 10.503 ( -5.018 -5.018 0.000) 7.096 10.832 ( -9.281 -9.281 0.000) 13.126 11.112 ( -4.888 -4.888 0.000) 6.913 11.503 ( 8.088 8.088 0.000) 11.438 11.987 ( 21.955 21.955 0.000) 31.049 20.936 ( -6.995 -6.995 0.000) 9.892 20.961 ( -6.169 -6.169 0.000) 8.724 21.734 ( 13.463 13.463 0.000) 19.039 ======================= Grid point 29 (15/84) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.025 ( 11.102 14.093 0.000) 17.941 4.176 ( 3.469 67.129 0.000) 67.219 5.308 ( 45.786 0.701 0.000) 45.791 9.387 ( -15.568 -29.720 0.000) 33.550 10.387 ( -5.741 -2.935 0.000) 6.447 10.701 ( -2.650 -11.257 0.000) 11.565 11.022 ( -3.821 -6.443 0.000) 7.491 11.653 ( 6.991 6.975 0.000) 9.875 12.449 ( 23.522 28.368 0.000) 36.852 20.768 ( -9.214 -7.282 0.000) 11.744 20.808 ( -8.714 -5.126 0.000) 10.110 22.047 ( 16.107 17.830 0.000) 24.028 ======================= Grid point 30 (16/84) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.173 ( 3.504 12.722 0.000) 13.196 4.235 ( 2.114 66.043 0.000) 66.077 6.008 ( 23.977 -2.351 0.000) 24.092 9.106 ( -12.233 -38.359 0.000) 40.262 10.294 ( -2.921 -1.228 0.000) 3.169 10.730 ( 4.767 -8.982 0.000) 10.168 10.951 ( -3.279 -7.523 0.000) 8.207 11.797 ( 7.008 8.701 0.000) 11.173 12.969 ( 27.777 23.793 0.000) 36.575 20.572 ( -9.553 -7.614 0.000) 12.216 20.610 ( -10.486 -4.179 0.000) 11.288 22.347 ( 12.301 24.516 0.000) 27.429 ======================= Grid point 31 (17/84) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.189 ( -1.223 12.202 0.000) 12.263 4.261 ( 0.573 67.223 0.000) 67.225 6.345 ( 10.337 -5.980 0.000) 11.942 8.891 ( -8.462 -45.444 0.000) 46.225 10.272 ( 0.578 0.198 0.000) 0.611 10.846 ( 5.361 -4.348 0.000) 6.903 10.885 ( -3.128 -8.377 0.000) 8.942 11.928 ( 5.446 13.827 0.000) 14.861 13.564 ( 28.749 11.312 0.000) 30.894 20.392 ( -10.390 -3.842 0.000) 11.078 20.395 ( -7.354 -7.914 0.000) 10.803 22.530 ( 5.664 32.308 0.000) 32.801 ======================= Grid point 32 (18/84) ======================= q-point: ( 0.46 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.157 ( -1.187 12.110 0.000) 12.168 4.265 ( -0.016 68.682 0.000) 68.682 6.471 ( 2.701 -9.160 0.000) 9.550 8.772 ( -2.985 -49.034 0.000) 49.125 10.296 ( 1.047 1.366 0.000) 1.721 10.833 ( -1.569 -9.204 0.000) 9.337 10.923 ( 1.947 -1.660 0.000) 2.559 12.006 ( 1.985 18.307 0.000) 18.415 14.027 ( 13.374 1.548 0.000) 13.463 20.223 ( -4.920 -3.880 0.000) 6.266 20.288 ( -2.783 -8.096 0.000) 8.561 22.595 ( 1.227 37.960 0.000) 37.980 ======================= Grid point 42 (19/84) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.274 ( 9.407 9.407 0.000) 13.304 5.222 ( 22.197 22.197 0.000) 31.392 5.412 ( 21.979 21.979 0.000) 31.083 8.823 ( -25.619 -25.619 0.000) 36.231 10.316 ( -3.587 -3.587 0.000) 5.073 10.561 ( -1.259 -1.259 0.000) 1.780 10.900 ( -5.251 -5.251 0.000) 7.426 11.787 ( 6.573 6.573 0.000) 9.295 12.976 ( 23.247 23.247 0.000) 32.877 20.593 ( -9.594 -9.594 0.000) 13.567 20.684 ( -6.988 -6.988 0.000) 9.882 22.467 ( 22.264 22.264 0.000) 31.486 ======================= Grid point 43 (20/84) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.395 ( 2.597 8.237 0.000) 8.637 5.349 ( 3.276 42.129 0.000) 42.257 5.952 ( 18.585 -3.138 0.000) 18.848 8.340 ( -21.154 -35.856 0.000) 41.631 10.262 ( -1.233 -1.645 0.000) 2.056 10.659 ( 9.722 3.404 0.000) 10.301 10.806 ( -4.198 -6.314 0.000) 7.583 11.940 ( 8.715 5.970 0.000) 10.564 13.409 ( 19.591 18.755 0.000) 27.121 20.389 ( -9.951 -10.045 0.000) 14.140 20.521 ( -8.948 -4.754 0.000) 10.132 22.907 ( 19.619 29.101 0.000) 35.097 ======================= Grid point 44 (21/84) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.398 ( -1.588 7.536 0.000) 7.701 5.415 ( 3.153 44.273 0.000) 44.385 6.152 ( 2.351 -13.721 0.000) 13.921 7.992 ( -12.247 -40.122 0.000) 41.949 10.273 ( 2.115 -0.098 0.000) 2.117 10.722 ( -4.073 -7.051 0.000) 8.143 10.876 ( 9.537 8.773 0.000) 12.958 12.141 ( 10.212 7.577 0.000) 12.716 13.787 ( 17.146 9.970 0.000) 19.834 20.204 ( -7.666 -10.467 0.000) 12.974 20.328 ( -9.455 -2.949 0.000) 9.904 23.236 ( 12.150 35.439 0.000) 37.463 ======================= Grid point 45 (22/84) ======================= q-point: ( 0.46 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.362 ( -1.263 7.229 0.000) 7.338 5.468 ( 1.605 48.261 0.000) 48.287 6.132 ( -1.978 -26.791 0.000) 26.864 7.839 ( -3.312 -38.438 0.000) 38.581 10.326 ( 2.144 1.487 0.000) 2.609 10.649 ( -2.443 -8.149 0.000) 8.507 11.008 ( 3.184 11.437 0.000) 11.872 12.305 ( 4.642 11.165 0.000) 12.092 14.063 ( 8.077 1.925 0.000) 8.303 20.093 ( -2.903 -10.729 0.000) 11.114 20.173 ( -4.552 -1.909 0.000) 4.936 23.399 ( 3.908 39.414 0.000) 39.607 ======================= Grid point 56 (23/84) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.496 ( 1.689 1.689 0.000) 2.389 5.814 ( -5.635 -5.635 0.000) 7.970 6.003 ( 16.924 16.924 0.000) 23.934 7.694 ( -25.387 -25.387 0.000) 35.903 10.242 ( 0.123 0.123 0.000) 0.173 10.696 ( -4.610 -4.610 0.000) 6.519 10.875 ( 16.334 16.334 0.000) 23.100 12.068 ( 7.292 7.292 0.000) 10.313 13.733 ( 13.178 13.178 0.000) 18.637 20.175 ( -10.441 -10.441 0.000) 14.766 20.405 ( -6.824 -6.824 0.000) 9.651 23.481 ( 25.819 25.819 0.000) 36.514 ======================= Grid point 57 (24/84) ======================= q-point: ( 0.38 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.483 ( -2.162 0.927 0.000) 2.352 5.564 ( -16.021 -32.452 0.000) 36.191 6.230 ( 5.468 24.693 0.000) 25.291 7.364 ( -7.492 -19.918 0.000) 21.281 10.276 ( 3.116 0.696 0.000) 3.193 10.607 ( -4.343 -4.309 0.000) 6.118 11.219 ( 14.469 23.768 0.000) 27.826 12.261 ( 11.551 4.844 0.000) 12.526 13.960 ( 9.261 6.943 0.000) 11.574 19.981 ( -8.061 -10.914 0.000) 13.568 20.247 ( -8.137 -5.386 0.000) 9.757 23.922 ( 16.721 30.469 0.000) 34.756 ======================= Grid point 58 (25/84) ======================= q-point: ( 0.46 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.440 ( -1.457 0.486 0.000) 1.536 5.281 ( -8.696 -52.252 0.000) 52.971 6.287 ( 1.056 30.768 0.000) 30.786 7.303 ( -0.517 -15.172 0.000) 15.181 10.354 ( 3.630 1.240 0.000) 3.836 10.520 ( -3.655 -4.326 0.000) 5.663 11.410 ( 4.312 27.163 0.000) 27.503 12.474 ( 6.779 5.755 0.000) 8.893 14.102 ( 4.045 1.839 0.000) 4.443 19.864 ( -3.057 -11.213 0.000) 11.622 20.111 ( -4.090 -4.803 0.000) 6.308 24.151 ( 5.630 33.045 0.000) 33.521 ======================= Grid point 70 (26/84) ======================= q-point: ( 0.38 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.457 ( -2.769 -2.769 0.000) 3.917 4.773 ( -39.256 -39.256 0.000) 55.516 6.622 ( 12.964 12.964 0.000) 18.333 7.139 ( -4.184 -4.184 0.000) 5.917 10.310 ( 2.707 2.707 0.000) 3.828 10.533 ( -3.256 -3.256 0.000) 4.604 11.740 ( 23.193 23.193 0.000) 32.800 12.372 ( 7.171 7.171 0.000) 10.141 14.074 ( 4.274 4.274 0.000) 6.044 19.778 ( -8.453 -8.453 0.000) 11.954 20.108 ( -7.680 -7.680 0.000) 10.861 24.441 ( 19.690 19.690 0.000) 27.846 ======================= Grid point 71 (27/84) ======================= q-point: ( 0.46 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.404 ( -1.692 -3.169 0.000) 3.592 4.153 ( -17.498 -53.408 0.000) 56.202 6.775 ( 3.008 17.424 0.000) 17.682 7.118 ( 0.527 -4.800 0.000) 4.829 10.379 ( 3.673 1.377 0.000) 3.923 10.465 ( -3.442 -1.516 0.000) 3.761 12.034 ( 5.728 30.784 0.000) 31.312 12.555 ( 7.624 2.545 0.000) 8.037 14.135 ( 1.639 1.307 0.000) 2.097 19.655 ( -3.212 -8.705 0.000) 9.279 19.975 ( -4.090 -7.884 0.000) 8.881 24.711 ( 6.674 21.207 0.000) 22.232 ======================= Grid point 84 (28/84) ======================= q-point: ( 0.46 0.46 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.345 ( -1.856 -1.856 0.000) 2.625 3.301 ( -24.625 -24.625 0.000) 34.825 7.007 ( 5.139 5.139 0.000) 7.268 7.066 ( -0.606 -0.606 0.000) 0.856 10.411 ( 1.587 1.587 0.000) 2.244 10.440 ( -1.200 -1.200 0.000) 1.698 12.512 ( 11.465 11.465 0.000) 16.214 12.599 ( 3.215 3.215 0.000) 4.546 14.154 ( 0.481 0.481 0.000) 0.681 19.529 ( -3.313 -3.313 0.000) 4.686 19.835 ( -4.365 -4.365 0.000) 6.173 25.003 ( 7.190 7.190 0.000) 10.168 ======================= Grid point 183 (29/84) ======================= q-point: ( 0.08 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.318 ( 35.556 35.556 35.556) 61.584 2.318 ( 35.556 35.556 35.556) 61.584 2.399 ( 35.045 35.045 35.045) 60.700 10.457 ( -9.362 -9.362 -9.362) 16.215 10.457 ( -9.362 -9.362 -9.362) 16.215 10.991 ( -11.227 -11.227 -11.227) 19.446 11.252 ( 13.925 13.925 13.925) 24.119 11.434 ( 3.441 3.441 3.441) 5.960 11.434 ( 3.441 3.441 3.441) 5.960 21.123 ( -3.403 -3.403 -3.403) 5.895 21.123 ( -3.403 -3.403 -3.403) 5.895 21.405 ( 5.980 5.980 5.980) 10.358 ======================= Grid point 184 (30/84) ======================= q-point: ( 0.15 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.516 ( 10.429 44.354 44.354) 63.588 2.556 ( 9.891 44.333 44.333) 63.472 4.067 ( 73.823 2.390 2.390) 73.900 10.156 ( -12.794 -12.038 -12.038) 21.296 10.331 ( -9.671 -11.676 -11.676) 19.136 10.796 ( -6.844 -15.973 -15.973) 23.604 11.382 ( -2.274 7.108 7.108) 10.307 11.576 ( 17.046 17.742 17.742) 30.334 11.583 ( 6.241 6.416 6.416) 11.013 21.017 ( -6.990 -3.254 -3.254) 8.369 21.019 ( -6.421 -3.178 -3.178) 7.838 21.591 ( 11.795 6.598 6.598) 15.039 ======================= Grid point 185 (31/84) ======================= q-point: ( 0.23 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.707 ( 7.670 40.086 40.086) 57.207 2.739 ( 7.348 40.192 40.192) 57.312 5.304 ( 47.224 -0.073 -0.073) 47.224 9.967 ( -5.265 -15.423 -15.423) 22.437 10.080 ( -13.749 -14.536 -14.536) 24.731 10.713 ( -1.909 -16.915 -16.915) 23.997 11.329 ( -2.615 7.028 7.028) 10.277 11.691 ( 4.414 11.329 11.329) 16.618 12.005 ( 25.756 19.056 19.056) 37.278 20.842 ( -9.992 -3.511 -3.511) 11.157 20.869 ( -7.999 -2.573 -2.573) 8.788 21.855 ( 12.961 8.682 8.682) 17.853 ======================= Grid point 186 (32/84) ======================= q-point: ( 0.31 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.815 ( 2.876 38.057 38.057) 53.898 2.840 ( 2.498 38.155 38.155) 54.017 6.032 ( 25.479 -1.420 -1.420) 25.558 9.809 ( -12.050 -19.317 -19.317) 29.858 9.934 ( 1.394 -14.789 -14.789) 20.961 10.697 ( -0.036 -16.350 -16.350) 23.122 11.278 ( -2.397 7.171 7.171) 10.420 11.765 ( 2.849 15.400 15.400) 21.964 12.646 ( 36.521 14.290 14.290) 41.739 20.618 ( -11.699 -4.465 -4.465) 13.294 20.704 ( -7.820 -1.392 -1.392) 8.064 22.075 ( 7.755 12.663 12.663) 19.515 ======================= Grid point 187 (33/84) ======================= q-point: ( 0.38 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.836 ( -0.390 37.949 37.949) 53.669 2.852 ( -0.905 37.773 37.773) 53.427 6.403 ( 12.044 -2.342 -2.342) 12.491 9.605 ( -7.669 -24.043 -24.043) 34.856 9.992 ( 3.535 -12.013 -12.013) 17.353 10.703 ( 0.481 -15.879 -15.879) 22.461 11.232 ( -1.996 7.525 7.525) 10.827 11.809 ( 1.530 18.368 18.368) 26.021 13.427 ( 37.017 6.468 6.468) 38.130 20.381 ( -10.910 -6.097 -6.097) 13.906 20.563 ( -5.732 0.048 0.048) 5.732 22.160 ( 0.924 18.031 18.031) 25.517 ======================= Grid point 188 (34/84) ======================= q-point: ( 0.46 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.821 ( -0.630 38.470 38.470) 54.409 2.828 ( -0.881 37.991 37.991) 53.735 6.557 ( 3.546 -2.946 -2.946) 5.471 9.500 ( -2.574 -27.053 -27.053) 38.346 10.051 ( 1.715 -10.150 -10.150) 14.456 10.711 ( 0.256 -15.690 -15.690) 22.191 11.202 ( -0.850 7.940 7.940) 11.261 11.829 ( 0.469 19.961 19.961) 28.233 14.011 ( 16.496 0.813 0.813) 16.536 20.208 ( -4.972 -7.545 -7.545) 11.772 20.481 ( -2.109 1.060 1.060) 2.588 22.149 ( -0.950 22.359 22.359) 31.635 ======================= Grid point 197 (35/84) ======================= q-point: ( 0.15 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.715 ( 9.525 9.525 74.024) 75.240 4.107 ( 34.523 34.523 2.617) 48.893 4.158 ( 38.948 38.948 2.524) 55.139 9.762 ( -20.325 -20.325 0.274) 28.745 10.215 ( -5.631 -5.631 -22.506) 23.873 10.482 ( -12.747 -12.747 -22.442) 28.785 11.431 ( -0.854 -0.854 19.501) 19.538 11.699 ( 7.977 7.977 11.900) 16.398 12.040 ( 22.754 22.754 7.243) 32.984 20.914 ( -6.577 -6.577 -2.272) 9.575 20.926 ( -6.084 -6.084 -3.356) 9.236 21.801 ( 13.707 13.707 6.737) 20.522 ======================= Grid point 198 (36/84) ======================= q-point: ( 0.23 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.890 ( 6.933 8.816 67.511) 68.436 4.205 ( 3.458 65.512 2.654) 65.656 5.308 ( 44.768 0.729 0.016) 44.774 9.396 ( -15.603 -29.311 0.841) 33.216 10.025 ( -10.471 -0.448 -28.100) 29.991 10.347 ( -1.794 -17.713 -22.832) 28.953 11.400 ( -1.634 1.141 22.504) 22.592 11.867 ( 8.193 7.523 15.593) 19.154 12.510 ( 23.757 26.247 6.910) 36.070 20.746 ( -9.803 -5.796 -2.210) 11.601 20.785 ( -7.268 -5.469 -2.011) 9.315 22.125 ( 17.013 17.515 7.844) 25.647 ======================= Grid point 199 (37/84) ======================= q-point: ( 0.31 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.984 ( 2.255 8.198 64.422) 64.981 4.262 ( 2.000 64.465 2.514) 64.545 5.992 ( 23.496 -2.361 -1.437) 23.658 9.099 ( -13.816 -38.039 -1.087) 40.485 9.890 ( -2.036 0.252 -31.477) 31.544 10.345 ( 0.905 -16.839 -22.506) 28.123 11.373 ( -0.983 3.389 22.965) 23.234 12.027 ( 7.320 9.805 17.948) 21.722 13.029 ( 27.291 21.028 6.207) 35.007 20.521 ( -11.951 -4.926 -3.154) 13.305 20.639 ( -6.811 -4.690 1.105) 8.344 22.450 ( 13.828 23.443 10.259) 29.086 ======================= Grid point 200 (38/84) ======================= q-point: ( 0.38 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.993 ( -0.892 7.754 64.198) 64.671 4.286 ( 0.465 65.656 2.285) 65.697 6.323 ( 10.145 -5.853 -2.180) 11.913 8.844 ( -10.512 -43.664 -5.740) 45.277 9.913 ( 3.057 -0.174 -28.805) 28.967 10.368 ( 1.116 -15.815 -23.549) 28.388 11.358 ( -0.480 5.956 22.556) 23.334 12.155 ( 4.964 14.115 17.999) 23.406 13.602 ( 27.130 9.940 3.695) 29.129 20.276 ( -11.397 -4.102 -4.386) 12.882 20.517 ( -4.950 -4.125 5.026) 8.172 22.666 ( 7.211 30.300 13.364) 33.893 ======================= Grid point 201 (39/84) ======================= q-point: ( 0.46 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.970 ( -0.828 7.538 64.900) 65.342 4.288 ( -0.052 67.124 2.150) 67.159 6.447 ( 2.660 -8.809 -2.477) 9.530 8.693 ( -3.851 -46.344 -9.021) 47.371 9.970 ( 1.781 -0.525 -26.310) 26.376 10.384 ( 0.371 -15.378 -23.929) 28.447 11.352 ( -0.158 7.705 22.219) 23.517 12.224 ( 1.704 17.257 17.440) 24.594 14.034 ( 12.355 1.451 0.677) 12.459 20.095 ( -5.186 -3.453 -5.118) 8.063 20.446 ( -1.832 -3.926 7.887) 8.999 22.755 ( 1.890 35.216 15.630) 38.575 ======================= Grid point 211 (40/84) ======================= q-point: ( 0.23 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.050 ( 6.224 6.224 62.154) 62.774 5.228 ( 21.724 21.724 0.639) 30.729 5.407 ( 21.398 21.398 -0.426) 30.264 8.831 ( -26.007 -26.007 0.683) 36.786 10.013 ( -3.900 -3.900 -24.802) 25.408 10.054 ( -6.567 -6.567 -30.622) 31.999 11.417 ( 1.150 1.150 31.195) 31.237 12.019 ( 7.617 7.617 17.778) 20.787 13.012 ( 22.505 22.505 3.781) 32.051 20.603 ( -8.401 -8.401 0.853) 11.912 20.660 ( -6.446 -6.446 -2.174) 9.371 22.545 ( 22.574 22.574 7.827) 32.869 ======================= Grid point 212 (41/84) ======================= q-point: ( 0.31 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.130 ( 1.656 5.546 59.661) 59.941 5.349 ( 3.135 41.093 0.113) 41.213 5.937 ( 18.222 -3.126 -1.454) 18.545 8.335 ( -22.131 -36.036 -0.784) 42.297 9.906 ( -1.682 1.322 -30.806) 30.880 10.061 ( 1.990 -9.915 -29.000) 30.713 11.458 ( 2.812 5.731 34.548) 35.133 12.188 ( 9.038 6.609 19.776) 22.726 13.434 ( 19.205 17.737 2.558) 26.267 20.408 ( -10.346 -6.521 1.657) 12.341 20.530 ( -6.117 -5.500 1.233) 8.318 22.997 ( 20.487 28.892 9.014) 36.548 ======================= Grid point 213 (42/84) ======================= q-point: ( 0.38 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.129 ( -1.285 5.005 59.580) 59.804 5.414 ( 3.073 43.220 -0.108) 43.329 6.134 ( 2.491 -13.386 -1.718) 13.724 7.961 ( -13.634 -39.965 -3.411) 42.364 9.927 ( 2.665 1.088 -30.490) 30.625 10.103 ( 1.755 -9.199 -30.973) 32.358 11.518 ( 2.668 10.589 34.111) 35.816 12.378 ( 8.875 8.103 19.152) 22.611 13.803 ( 16.601 9.194 1.600) 19.044 20.197 ( -9.760 -4.336 0.785) 10.709 20.416 ( -4.823 -5.195 7.132) 10.056 23.345 ( 13.085 34.783 10.856) 38.716 ======================= Grid point 214 (43/84) ======================= q-point: ( 0.46 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.101 ( -0.959 4.706 60.251) 60.442 5.465 ( 1.594 47.262 -0.180) 47.289 6.119 ( -1.833 -26.163 -1.368) 26.263 7.785 ( -4.007 -38.187 -5.860) 38.842 9.979 ( 1.733 0.937 -30.025) 30.089 10.124 ( 0.359 -9.248 -30.798) 32.159 11.556 ( 0.969 13.213 33.108) 35.661 12.514 ( 3.683 11.139 17.034) 20.683 14.068 ( 7.624 1.752 0.415) 7.834 20.043 ( -4.392 -2.604 0.127) 5.108 20.345 ( -1.882 -5.388 11.607) 12.934 23.523 ( 4.309 38.480 12.228) 40.606 ======================= Grid point 225 (44/84) ======================= q-point: ( 0.31 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.198 ( 1.049 1.049 57.677) 57.696 5.798 ( -5.514 -5.514 -1.479) 7.937 5.990 ( 16.421 16.421 -1.286) 23.258 7.681 ( -26.014 -26.014 -1.480) 36.820 9.922 ( -0.437 -0.437 -27.033) 27.040 9.972 ( 1.154 1.154 -38.216) 38.250 11.620 ( 9.709 9.709 39.641) 41.951 12.324 ( 7.380 7.380 20.728) 23.207 13.744 ( 12.742 12.742 1.134) 18.056 20.252 ( -8.414 -8.414 7.384) 14.004 20.422 ( -5.282 -5.282 1.916) 7.712 23.574 ( 26.299 26.299 9.299) 38.337 ======================= Grid point 226 (45/84) ======================= q-point: ( 0.38 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.185 ( -1.720 0.549 57.941) 57.969 5.555 ( -15.693 -32.143 -0.951) 35.782 6.209 ( 5.379 23.885 -1.934) 24.559 7.331 ( -8.656 -20.118 -3.471) 22.174 9.940 ( 2.192 0.192 -28.780) 28.864 10.012 ( 1.673 -0.283 -38.325) 38.363 11.811 ( 7.759 17.871 37.644) 42.386 12.507 ( 10.260 5.076 20.035) 23.074 13.965 ( 9.007 6.540 0.464) 11.140 20.083 ( -7.700 -6.678 9.635) 14.025 20.314 ( -4.960 -4.949 6.867) 9.811 24.028 ( 17.385 30.702 10.468) 36.803 ======================= Grid point 227 (46/84) ======================= q-point: ( 0.46 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.151 ( -1.118 0.254 58.764) 58.775 5.276 ( -8.568 -52.093 -0.474) 52.795 6.267 ( 1.088 29.956 -1.931) 30.038 7.250 ( -1.094 -15.065 -5.544) 16.090 9.991 ( 2.379 -0.206 -31.129) 31.220 10.024 ( -0.579 -0.567 -35.248) 35.257 11.913 ( 2.318 21.572 35.009) 41.187 12.687 ( 5.506 6.090 17.481) 19.313 14.103 ( 3.888 1.615 0.020) 4.210 19.962 ( -3.395 -5.276 9.758) 11.601 20.237 ( -2.152 -5.314 12.011) 13.310 24.268 ( 5.929 33.150 11.455) 35.571 ======================= Grid point 239 (47/84) ======================= q-point: ( 0.38 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.164 ( -2.082 -2.082 58.636) 58.710 4.767 ( -39.145 -39.145 -0.545) 55.362 6.588 ( 12.489 12.489 -3.468) 18.000 7.100 ( -4.620 -4.620 -4.026) 7.674 9.962 ( 1.948 1.948 -30.168) 30.294 10.029 ( 0.976 0.976 -38.133) 38.158 12.191 ( 16.538 16.538 32.704) 40.207 12.616 ( 6.590 6.590 19.955) 22.024 14.072 ( 4.067 4.067 -0.126) 5.752 19.944 ( -6.321 -6.321 15.912) 18.251 20.208 ( -5.293 -5.293 9.851) 12.372 24.556 ( 20.215 20.215 11.247) 30.721 ======================= Grid point 240 (48/84) ======================= q-point: ( 0.46 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.125 ( -1.257 -2.294 59.735) 59.792 4.148 ( -17.474 -53.465 -0.423) 56.250 6.738 ( 3.096 16.604 -3.785) 17.309 7.068 ( 0.031 -4.712 -5.276) 7.074 10.004 ( 1.712 1.549 -32.810) 32.891 10.037 ( -0.243 0.839 -34.617) 34.628 12.394 ( 3.842 23.525 27.846) 36.655 12.775 ( 6.366 2.863 18.123) 19.421 14.131 ( 1.564 1.116 -0.403) 1.963 19.846 ( -2.734 -5.463 18.267) 19.262 20.121 ( -2.539 -5.894 14.169) 15.554 24.835 ( 6.928 21.669 12.043) 25.741 ======================= Grid point 253 (49/84) ======================= q-point: ( 0.46 0.46 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.082 ( -1.337 -1.337 61.046) 61.075 3.295 ( -24.680 -24.680 -0.556) 34.907 6.957 ( 4.929 4.929 -5.159) 8.672 7.016 ( -0.752 -0.752 -5.255) 5.361 10.036 ( 1.185 1.185 -32.800) 32.843 10.047 ( 0.114 0.114 -33.901) 33.901 12.761 ( 8.653 8.653 20.193) 23.611 12.821 ( 2.896 2.896 18.324) 18.776 14.147 ( 0.407 0.407 -0.676) 0.888 19.762 ( -2.365 -2.365 22.233) 22.483 20.022 ( -2.947 -2.947 17.892) 18.371 25.135 ( 7.424 7.424 12.753) 16.519 ======================= Grid point 366 (50/84) ======================= q-point: ( 0.15 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.156 ( 22.272 22.272 22.272) 38.576 4.223 ( 25.628 25.628 25.628) 44.389 4.223 ( 25.628 25.628 25.628) 44.389 9.747 ( -12.444 -12.444 -12.444) 21.554 9.747 ( -12.444 -12.444 -12.444) 21.554 10.026 ( -19.710 -19.710 -19.710) 34.140 11.813 ( 9.180 9.180 9.180) 15.900 11.813 ( 9.180 9.180 9.180) 15.900 12.338 ( 19.506 19.506 19.506) 33.786 20.835 ( -5.434 -5.434 -5.434) 9.412 20.835 ( -5.434 -5.434 -5.434) 9.412 22.019 ( 14.623 14.623 14.623) 25.328 ======================= Grid point 367 (51/84) ======================= q-point: ( 0.23 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.258 ( 3.759 31.438 31.438) 44.618 4.292 ( 2.976 34.408 34.408) 48.752 5.313 ( 42.161 0.663 0.663) 42.171 9.416 ( -10.575 -11.380 -11.380) 19.258 9.424 ( -23.412 -16.686 -16.686) 33.241 9.864 ( -0.282 -23.327 -23.327) 32.991 11.806 ( -0.808 16.261 16.261) 23.010 12.134 ( 13.771 10.264 10.264) 20.008 12.764 ( 21.826 17.283 17.283) 32.768 20.682 ( -9.188 -3.641 -3.641) 10.532 20.726 ( -5.050 -3.989 -3.989) 7.571 22.378 ( 19.557 16.968 16.968) 30.956 ======================= Grid point 368 (52/84) ======================= q-point: ( 0.31 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.315 ( 1.676 31.307 31.307) 44.306 4.340 ( 1.730 33.346 33.346) 47.189 5.952 ( 21.597 -2.406 -2.406) 21.863 8.943 ( -22.565 -22.346 -22.346) 38.831 9.321 ( 0.632 -14.447 -14.447) 20.440 9.878 ( 1.057 -22.161 -22.161) 31.358 11.792 ( -0.438 18.170 18.170) 25.700 12.371 ( 9.409 13.906 13.906) 21.801 13.230 ( 23.712 13.061 13.061) 30.057 20.467 ( -11.566 -1.856 -1.856) 11.861 20.638 ( -3.495 -1.351 -1.351) 3.984 22.771 ( 17.806 21.053 21.053) 34.692 ======================= Grid point 369 (53/84) ======================= q-point: ( 0.38 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.326 ( -0.347 31.858 31.858) 45.056 4.364 ( 0.717 33.668 33.668) 47.619 6.250 ( 8.799 -5.369 -5.369) 11.622 8.537 ( -16.415 -27.589 -27.589) 42.329 9.395 ( 5.207 -16.254 -16.254) 23.570 9.899 ( 0.845 -21.844 -21.844) 30.903 11.789 ( 0.074 19.841 19.841) 28.059 12.519 ( 5.238 15.895 15.895) 23.081 13.709 ( 21.985 6.544 6.544) 23.854 20.229 ( -10.993 -0.049 -0.049) 10.993 20.585 ( -1.746 1.566 1.566) 2.820 23.071 ( 11.124 25.845 25.845) 38.205 ======================= Grid point 370 (54/84) ======================= q-point: ( 0.46 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.315 ( -0.475 32.372 32.372) 45.783 4.373 ( 0.173 34.429 34.429) 48.690 6.354 ( 2.172 -7.439 -7.439) 10.743 8.301 ( -6.067 -31.094 -31.094) 44.390 9.483 ( 2.578 -16.763 -16.763) 23.847 9.911 ( 0.310 -21.789 -21.789) 30.816 11.792 ( 0.117 20.908 20.908) 29.569 12.589 ( 1.648 16.693 16.693) 23.665 14.053 ( 9.687 1.170 1.170) 9.827 20.057 ( -4.886 1.372 1.372) 5.256 20.563 ( -0.484 3.482 3.482) 4.948 23.220 ( 3.561 29.203 29.203) 41.453 ======================= Grid point 380 (55/84) ======================= q-point: ( 0.23 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.349 ( 3.276 3.276 58.747) 58.930 5.253 ( 20.508 20.508 2.031) 29.073 5.396 ( 19.923 19.923 -0.674) 28.184 8.838 ( -27.794 -27.794 -0.841) 39.316 9.430 ( -8.962 -8.962 -29.719) 32.308 9.467 ( -5.640 -5.640 -26.211) 27.398 12.034 ( 6.453 6.453 28.294) 29.729 12.354 ( 11.959 11.959 13.085) 21.383 13.148 ( 19.005 19.005 9.768) 28.597 20.614 ( -4.545 -4.545 -1.775) 6.668 20.618 ( -4.838 -4.838 0.351) 6.851 22.800 ( 23.461 23.461 17.351) 37.442 ======================= Grid point 381 (56/84) ======================= q-point: ( 0.31 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.398 ( 1.411 3.080 58.253) 58.351 5.356 ( 2.737 38.520 0.678) 38.623 5.896 ( 16.769 -3.069 -2.444) 17.222 8.283 ( -26.020 -35.307 -5.617) 44.217 9.235 ( -4.227 -0.289 -31.463) 31.747 9.485 ( 1.381 -15.328 -27.336) 31.371 12.123 ( 3.198 14.367 29.745) 33.188 12.594 ( 10.880 8.807 18.013) 22.812 13.520 ( 17.372 14.158 5.934) 23.184 20.453 ( -8.906 0.463 3.016) 9.414 20.584 ( -0.596 -3.852 3.833) 5.467 23.286 ( 22.745 28.255 19.371) 41.121 ======================= Grid point 382 (57/84) ======================= q-point: ( 0.38 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.409 ( -0.145 2.846 59.703) 59.771 5.415 ( 2.895 40.548 0.344) 40.653 6.077 ( 2.431 -12.165 -4.105) 13.067 7.814 ( -18.791 -37.457 -12.676) 43.782 9.233 ( 2.680 -2.809 -32.754) 32.983 9.509 ( 0.949 -15.533 -27.634) 31.714 12.178 ( 2.166 18.351 29.768) 35.037 12.781 ( 7.121 9.827 18.357) 22.006 13.852 ( 14.660 6.953 3.128) 16.524 20.269 ( -8.497 3.568 6.554) 11.309 20.585 ( 0.522 -1.307 8.978) 9.088 23.684 ( 15.493 32.875 22.278) 42.627 ======================= Grid point 383 (58/84) ======================= q-point: ( 0.46 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.402 ( -0.289 2.693 61.262) 61.321 5.465 ( 1.597 44.710 0.302) 44.740 6.068 ( -1.471 -23.600 -4.041) 23.989 7.550 ( -6.646 -35.884 -19.027) 41.156 9.290 ( 1.909 -4.321 -32.764) 33.103 9.522 ( 0.334 -15.585 -27.916) 31.973 12.209 ( 0.805 19.953 29.835) 35.901 12.879 ( 2.414 11.565 17.035) 20.731 14.079 ( 6.367 1.291 0.706) 6.535 20.136 ( -3.729 5.779 9.018) 11.341 20.596 ( 0.357 0.077 12.464) 12.470 23.899 ( 5.342 35.795 24.410) 43.654 ======================= Grid point 394 (59/84) ======================= q-point: ( 0.31 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.444 ( 1.302 1.302 58.211) 58.240 5.755 ( -5.280 -5.280 -2.685) 7.935 5.956 ( 14.688 14.688 -1.801) 20.850 7.617 ( -27.957 -27.957 -5.515) 39.920 9.221 ( -2.051 -2.051 -34.744) 34.865 9.292 ( -2.823 -2.823 -32.239) 32.485 12.379 ( 10.459 10.459 33.672) 36.777 12.763 ( 8.262 8.262 20.176) 23.315 13.781 ( 11.227 11.227 2.531) 16.078 20.464 ( -2.595 -2.595 12.781) 13.298 20.503 ( -0.828 -0.828 6.625) 6.728 23.869 ( 27.382 27.382 19.637) 43.418 ======================= Grid point 395 (60/84) ======================= q-point: ( 0.38 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.453 ( -0.261 1.316 60.133) 60.148 5.526 ( -14.809 -31.254 -1.857) 34.634 6.152 ( 4.972 21.025 -3.669) 21.914 7.197 ( -13.166 -20.427 -10.651) 26.534 9.207 ( 1.306 -0.486 -37.743) 37.769 9.285 ( 0.426 -6.219 -34.115) 34.680 12.540 ( 5.348 16.956 32.297) 36.867 12.937 ( 8.255 5.913 20.050) 22.475 13.979 ( 8.176 5.343 0.963) 9.815 20.355 ( -5.676 4.215 16.465) 17.919 20.542 ( 2.330 -2.569 15.485) 15.869 24.352 ( 18.988 31.014 21.257) 42.122 ======================= Grid point 396 (61/84) ======================= q-point: ( 0.46 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.444 ( -0.332 1.305 61.989) 62.003 5.261 ( -8.219 -51.455 -1.037) 52.118 6.210 ( 1.286 26.751 -3.568) 27.019 7.037 ( -3.540 -14.303 -17.020) 22.512 9.245 ( 1.436 -0.963 -38.535) 38.574 9.292 ( 0.157 -6.832 -34.862) 35.525 12.608 ( 1.500 18.960 31.141) 36.489 13.065 ( 3.542 6.886 17.816) 19.426 14.103 ( 3.437 1.012 0.048) 3.583 20.265 ( -2.549 6.194 19.640) 20.750 20.578 ( 0.968 -1.629 20.585) 20.672 24.618 ( 6.665 33.061 22.679) 40.642 ======================= Grid point 408 (62/84) ======================= q-point: ( 0.38 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.463 ( -0.161 -0.161 62.676) 62.677 4.751 ( -38.814 -38.814 -1.025) 54.900 6.471 ( 10.261 10.261 -8.299) 16.717 6.949 ( -6.362 -6.362 -11.939) 14.950 9.206 ( 0.698 0.698 -39.983) 39.995 9.233 ( -0.303 -0.303 -38.420) 38.423 12.842 ( 11.164 11.164 28.968) 32.991 13.047 ( 5.755 5.755 20.317) 21.887 14.069 ( 3.467 3.467 -0.173) 4.906 20.401 ( -0.667 -0.667 27.747) 27.763 20.514 ( 0.762 0.762 19.986) 20.015 24.898 ( 21.373 21.373 22.078) 37.431 ======================= Grid point 409 (63/84) ======================= q-point: ( 0.46 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.136 ( -17.394 -53.588 -0.779) 56.346 4.456 ( -0.324 -0.074 64.911) 64.912 6.605 ( 3.362 12.652 -9.885) 16.404 6.868 ( -2.132 -4.285 -15.851) 16.558 9.224 ( 0.431 -0.198 -40.114) 40.117 9.240 ( 0.705 -0.111 -40.462) 40.469 12.967 ( 2.191 15.475 25.982) 30.321 13.169 ( 4.333 3.540 18.667) 19.488 14.120 ( 1.356 0.578 -0.683) 1.625 20.366 ( -1.437 3.171 31.218) 31.411 20.543 ( 1.089 -1.273 26.693) 26.746 25.197 ( 7.514 22.619 23.104) 33.195 ======================= Grid point 422 (64/84) ======================= q-point: ( 0.46 0.46 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.279 ( -24.825 -24.825 -0.996) 35.122 4.450 ( -0.268 -0.268 67.445) 67.446 6.764 ( 3.736 3.736 -15.206) 16.098 6.816 ( -1.369 -1.369 -15.865) 15.982 9.239 ( 0.578 0.578 -41.537) 41.545 9.241 ( 0.336 0.336 -41.362) 41.364 13.194 ( 4.986 4.986 20.369) 21.555 13.221 ( 2.415 2.415 19.021) 19.325 14.127 ( 0.189 0.189 -1.187) 1.217 20.391 ( -0.006 -0.006 37.800) 37.800 20.536 ( 0.264 0.264 31.424) 31.426 25.514 ( 7.944 7.944 23.946) 26.450 ======================= Grid point 549 (65/84) ======================= q-point: ( 0.23 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.254 ( 15.151 15.151 15.151) 26.242 5.384 ( 12.015 12.015 12.015) 20.810 5.384 ( 12.015 12.015 12.015) 20.810 8.626 ( -25.733 -25.733 -25.733) 44.570 9.004 ( -12.108 -12.108 -12.108) 20.971 9.004 ( -12.108 -12.108 -12.108) 20.971 12.545 ( 13.836 13.836 13.836) 23.965 12.545 ( 13.836 13.836 13.836) 23.965 13.394 ( 13.468 13.468 13.468) 23.327 20.606 ( -0.179 -0.179 -0.179) 0.310 20.606 ( -0.179 -0.179 -0.179) 0.310 23.242 ( 25.109 25.109 25.109) 43.491 ======================= Grid point 550 (66/84) ======================= q-point: ( 0.31 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.366 ( 3.825 20.641 20.641) 29.440 5.390 ( 1.102 16.894 16.894) 23.917 5.842 ( 13.033 -2.629 -2.629) 13.553 7.998 ( -32.083 -24.171 -24.171) 46.880 8.755 ( -6.166 -13.257 -13.257) 19.736 8.978 ( 0.415 -21.286 -21.286) 30.106 12.646 ( 3.980 20.956 20.956) 29.902 12.896 ( 13.194 11.812 11.812) 21.287 13.671 ( 13.364 8.495 8.495) 17.970 20.549 ( -3.760 6.646 6.646) 10.123 20.667 ( 4.406 3.901 3.901) 7.060 23.775 ( 25.626 27.721 27.721) 46.836 ======================= Grid point 551 (67/84) ======================= q-point: ( 0.38 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.432 ( 2.501 22.432 22.432) 31.822 5.433 ( 3.008 17.092 17.092) 24.358 5.947 ( -0.661 -8.991 -8.991) 12.732 7.385 ( -26.544 -28.300 -28.300) 48.024 8.727 ( 2.098 -15.608 -15.608) 22.172 8.989 ( 0.534 -22.601 -22.601) 31.967 12.714 ( 2.646 22.499 22.499) 31.929 13.095 ( 6.632 12.242 12.242) 18.539 13.926 ( 11.058 3.911 3.911) 12.365 20.463 ( -4.240 12.173 12.173) 17.730 20.767 ( 4.727 8.412 8.412) 12.801 24.235 ( 18.321 30.662 30.662) 47.074 ======================= Grid point 552 (68/84) ======================= q-point: ( 0.46 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.467 ( 0.883 24.874 24.874) 35.188 5.496 ( 2.280 19.229 19.229) 27.289 5.901 ( -2.106 -15.196 -15.196) 21.594 6.981 ( -11.161 -33.240 -33.240) 48.316 8.781 ( 1.953 -16.201 -16.201) 22.995 8.997 ( 0.197 -23.024 -23.024) 32.562 12.751 ( 0.948 23.004 23.004) 32.546 13.180 ( 1.961 12.258 12.258) 17.447 14.094 ( 4.608 0.668 0.668) 4.704 20.395 ( -1.945 15.687 15.687) 22.269 20.839 ( 1.987 11.068 11.068) 15.778 24.493 ( 6.525 32.645 32.645) 46.626 ======================= Grid point 563 (69/84) ======================= q-point: ( 0.31 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.421 ( 2.408 2.408 33.695) 33.866 5.692 ( -5.380 -5.380 -3.372) 8.322 5.935 ( 12.127 12.127 1.672) 17.232 7.401 ( -30.479 -30.479 -17.328) 46.457 8.600 ( -5.764 -5.764 -24.269) 25.601 8.673 ( -6.739 -6.739 -26.824) 28.467 12.964 ( 10.141 10.141 23.207) 27.281 13.111 ( 9.973 9.973 13.728) 19.682 13.847 ( 8.369 8.369 3.923) 12.468 20.704 ( 5.653 5.653 12.857) 15.140 20.745 ( 5.841 5.841 13.935) 16.200 24.365 ( 28.521 28.521 28.235) 49.235 ======================= Grid point 564 (70/84) ======================= q-point: ( 0.38 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.465 ( -5.912 -13.593 16.241) 21.988 5.488 ( -4.926 -13.163 15.244) 20.734 6.079 ( 3.413 16.396 -0.409) 16.752 6.843 ( -23.236 -20.337 -24.734) 39.563 8.540 ( 0.632 -4.739 -25.496) 25.940 8.622 ( -0.726 -12.366 -29.480) 31.977 13.098 ( 4.029 15.269 22.054) 27.125 13.289 ( 6.852 7.173 14.220) 17.339 14.003 ( 6.645 3.502 1.347) 7.631 20.742 ( 0.338 14.253 21.117) 25.479 20.917 ( 8.190 6.147 19.964) 22.437 24.881 ( 20.705 30.975 29.686) 47.638 ======================= Grid point 565 (71/84) ======================= q-point: ( 0.46 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.228 ( -7.655 -49.107 -2.380) 49.757 5.534 ( 1.894 3.318 40.118) 40.300 6.135 ( 2.379 19.631 -1.972) 19.873 6.482 ( -10.866 -12.465 -37.346) 40.843 8.577 ( 1.674 -5.476 -24.427) 25.089 8.615 ( -0.127 -13.309 -30.397) 33.183 13.149 ( 1.163 16.122 21.644) 27.014 13.382 ( 2.249 7.717 12.985) 15.271 14.104 ( 2.751 0.246 -0.039) 2.762 20.741 ( -0.175 17.115 26.015) 31.141 21.039 ( 3.227 8.305 23.485) 25.118 25.174 ( 7.460 32.431 30.915) 45.423 ======================= Grid point 577 (72/84) ======================= q-point: ( 0.38 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.726 ( -38.239 -38.239 -1.390) 54.097 5.545 ( 3.759 3.759 35.278) 35.676 6.267 ( 3.694 3.694 -6.252) 8.147 6.549 ( -10.935 -10.935 -28.644) 32.552 8.486 ( -0.549 -0.549 -28.390) 28.401 8.492 ( -1.849 -1.849 -31.690) 31.798 13.337 ( 7.640 7.640 19.326) 22.141 13.408 ( 5.255 5.255 14.771) 16.535 14.066 ( 2.557 2.557 -0.155) 3.619 21.004 ( 7.852 7.852 27.169) 29.351 21.020 ( 6.701 6.701 31.282) 32.685 25.439 ( 22.407 22.407 30.138) 43.732 ======================= Grid point 578 (73/84) ======================= q-point: ( 0.46 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.117 ( -17.286 -53.700 -1.031) 56.423 5.617 ( 2.495 4.642 41.736) 42.068 6.295 ( -0.354 -0.558 -25.893) 25.901 6.373 ( -4.891 -2.582 -24.451) 25.069 8.472 ( -0.351 -1.363 -31.726) 31.757 8.511 ( 1.415 -1.728 -27.757) 27.847 13.417 ( 1.423 10.117 17.844) 20.562 13.503 ( 2.829 4.284 13.714) 14.644 14.104 ( 1.060 -0.165 -0.810) 1.345 21.055 ( 0.984 12.828 34.833) 37.133 21.174 ( 4.204 4.922 33.581) 34.199 25.757 ( 8.084 23.327 30.886) 39.540 ======================= Grid point 591 (74/84) ======================= q-point: ( 0.46 0.46 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.256 ( -25.016 -25.016 -1.222) 35.399 5.706 ( 3.229 3.229 50.738) 50.944 6.273 ( -1.766 -1.766 -30.943) 31.044 6.283 ( -3.119 -3.119 -38.001) 38.256 8.489 ( 0.703 0.703 -29.374) 29.391 8.493 ( 0.349 0.349 -28.647) 28.651 13.556 ( 2.832 2.832 14.789) 15.322 13.563 ( 2.162 2.162 14.123) 14.451 14.100 ( -0.138 -0.138 -1.415) 1.429 21.227 ( 2.828 2.828 41.922) 42.112 21.252 ( 3.459 3.459 37.215) 37.535 26.088 ( 8.409 8.409 31.404) 33.580 ======================= Grid point 732 (75/84) ======================= q-point: ( 0.31 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.622 ( -5.747 -5.747 -5.747) 9.955 5.622 ( -5.747 -5.747 -5.747) 9.955 6.175 ( 15.515 15.515 15.515) 26.873 6.892 ( -30.690 -30.690 -30.690) 53.156 8.266 ( -10.534 -10.534 -10.534) 18.246 8.266 ( -10.534 -10.534 -10.534) 18.246 13.330 ( 10.533 10.533 10.533) 18.243 13.330 ( 10.533 10.533 10.533) 18.243 13.937 ( 4.724 4.724 4.724) 8.182 21.006 ( 13.518 13.518 13.518) 23.414 21.006 ( 13.518 13.518 13.518) 23.414 24.973 ( 29.702 29.702 29.702) 51.446 ======================= Grid point 733 (76/84) ======================= q-point: ( 0.38 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.342 ( -16.476 -16.479 -16.479) 28.540 5.556 ( -2.964 -13.164 -13.164) 18.851 6.267 ( -27.565 -27.596 -27.596) 47.779 6.382 ( 5.511 20.430 20.430) 29.413 8.151 ( -3.007 -14.799 -14.799) 21.143 8.186 ( 1.832 -10.315 -10.315) 14.702 13.450 ( 3.947 12.841 12.841) 18.584 13.520 ( 6.110 8.590 8.590) 13.599 14.033 ( 4.342 1.490 1.490) 4.827 21.191 ( 6.760 21.787 21.787) 31.544 21.309 ( 12.764 17.498 17.498) 27.844 25.516 ( 22.065 30.962 30.962) 49.032 ======================= Grid point 734 (77/84) ======================= q-point: ( 0.46 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.064 ( -8.246 -24.544 -24.544) 35.676 5.503 ( -1.963 -30.210 -30.210) 42.768 5.864 ( -10.241 -13.696 -13.696) 21.910 6.445 ( 1.258 22.545 22.545) 31.908 8.116 ( -0.725 -16.270 -16.270) 23.021 8.247 ( 2.340 -9.579 -9.579) 13.748 13.500 ( 1.136 13.134 13.134) 18.610 13.597 ( 1.755 8.287 8.287) 11.850 14.100 ( 1.862 -0.406 -0.406) 1.949 21.282 ( 2.204 25.726 25.726) 36.449 21.494 ( 4.825 20.200 20.200) 28.972 25.831 ( 8.061 31.841 31.841) 45.746 ======================= Grid point 746 (78/84) ======================= q-point: ( 0.38 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.694 ( -37.186 -37.186 -1.850) 52.622 5.408 ( -5.197 -5.197 -33.806) 34.596 5.786 ( -18.914 -18.914 -43.765) 51.292 6.657 ( 6.902 6.902 29.881) 31.435 8.017 ( -1.500 -1.500 -13.973) 14.133 8.074 ( -0.450 -0.450 -12.891) 12.906 13.640 ( 5.793 5.793 10.808) 13.563 13.653 ( 5.050 5.050 9.370) 11.781 14.063 ( 1.443 1.443 -0.081) 2.042 21.575 ( 13.718 13.718 27.302) 33.492 21.621 ( 13.776 13.776 26.399) 32.809 26.082 ( 23.069 23.069 31.370) 45.259 ======================= Grid point 747 (79/84) ======================= q-point: ( 0.46 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.094 ( -17.146 -53.371 -1.262) 56.072 5.280 ( -5.551 -2.658 -54.591) 54.937 5.550 ( -4.379 -16.464 -38.850) 42.421 6.728 ( 1.217 5.597 34.611) 35.082 7.991 ( -0.741 2.362 -14.933) 15.137 8.126 ( 2.594 -3.273 -11.454) 12.192 13.703 ( 1.338 7.096 10.583) 12.812 13.730 ( 1.890 4.840 8.628) 10.072 14.087 ( 0.708 -0.809 -0.889) 1.395 21.743 ( 3.889 18.421 31.502) 36.700 21.839 ( 5.889 13.184 29.832) 33.143 26.412 ( 8.462 23.653 31.840) 40.557 ======================= Grid point 760 (80/84) ======================= q-point: ( 0.46 0.46 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.231 ( -25.196 -25.196 -1.175) 35.652 5.245 ( -2.283 -2.283 -55.491) 55.584 5.285 ( -6.240 -6.240 -56.961) 57.641 6.768 ( 0.218 0.218 41.242) 41.243 8.074 ( 2.122 2.122 -12.424) 12.781 8.092 ( 0.528 0.528 -11.994) 12.018 13.797 ( 2.058 2.058 8.886) 9.350 13.799 ( 1.954 1.954 9.095) 9.506 14.072 ( -0.496 -0.496 -1.326) 1.500 21.996 ( 5.760 5.760 34.205) 35.162 22.029 ( 5.404 5.404 35.214) 36.034 26.750 ( 8.679 8.679 32.110) 34.376 ======================= Grid point 915 (81/84) ======================= q-point: ( 0.38 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.616 ( -24.481 -24.481 -24.481) 42.402 4.616 ( -24.481 -24.481 -24.481) 42.402 4.926 ( -32.350 -32.350 -32.350) 56.032 7.123 ( 14.467 14.467 14.467) 25.058 7.904 ( -1.741 -1.741 -1.741) 3.015 7.904 ( -1.741 -1.741 -1.741) 3.015 13.799 ( 5.112 5.112 5.112) 8.853 13.799 ( 5.112 5.112 5.112) 8.853 14.063 ( 0.130 0.130 0.130) 0.225 22.067 ( 18.320 18.320 18.320) 31.731 22.067 ( 18.320 18.320 18.320) 31.731 26.658 ( 23.562 23.562 23.562) 40.810 ======================= Grid point 916 (82/84) ======================= q-point: ( 0.46 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.003 ( -17.394 -25.368 -25.368) 39.870 4.180 ( -16.607 -30.539 -30.539) 46.272 4.708 ( -1.463 -39.735 -39.735) 56.213 7.282 ( 2.889 17.589 17.589) 25.042 7.844 ( -1.299 0.022 0.022) 1.299 7.979 ( 3.290 -4.010 -4.010) 6.556 13.864 ( 1.516 5.484 5.484) 7.903 13.865 ( 1.516 4.865 4.865) 7.045 14.067 ( 0.175 -1.019 -1.019) 1.452 22.299 ( 5.714 21.817 21.817) 31.378 22.348 ( 7.172 19.467 19.467) 28.449 26.996 ( 8.699 23.895 23.895) 34.894 ======================= Grid point 929 (83/84) ======================= q-point: ( 0.46 0.46 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.210 ( -25.304 -25.304 -0.862) 35.796 4.071 ( -7.668 -7.668 -55.174) 56.229 4.098 ( -10.292 -10.292 -55.362) 57.243 7.458 ( 2.445 2.445 25.765) 25.996 7.918 ( 2.430 2.430 -4.101) 5.351 7.938 ( 0.783 0.783 -4.268) 4.409 13.933 ( 1.864 1.864 4.711) 5.398 13.937 ( 1.692 1.692 4.661) 5.240 14.048 ( -0.720 -0.720 -0.997) 1.424 22.599 ( 7.159 7.159 23.723) 25.794 22.624 ( 7.277 7.277 22.448) 24.695 27.339 ( 8.835 8.835 24.101) 27.147 ======================= Grid point 1098 (84/84) ======================= q-point: ( 0.46 0.46 0.46) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 2.71e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.188 ( -16.454 -16.454 -16.454) 28.499 3.188 ( -16.454 -16.454 -16.454) 28.499 3.224 ( -18.692 -18.692 -18.692) 32.376 7.796 ( 6.239 6.239 6.239) 10.806 7.888 ( 0.276 0.276 0.276) 0.478 7.888 ( 0.276 0.276 0.276) 0.478 13.997 ( 1.551 1.551 1.551) 2.686 13.997 ( 1.551 1.551 1.551) 2.686 14.031 ( -0.619 -0.619 -0.619) 1.072 22.933 ( 8.032 8.032 8.032) 13.912 22.933 ( 8.032 8.032 8.032) 13.912 27.686 ( 8.928 8.928 8.928) 15.464 =================== 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/26364 10.0 753.858 753.858 753.858 -0.000 -0.000 0.000 3/26364 20.0 273.067 273.067 273.067 -0.000 -0.000 0.000 3/26364 30.0 366.879 366.879 366.879 -0.000 -0.000 0.000 3/26364 40.0 322.650 322.650 322.650 -0.000 -0.000 0.000 3/26364 50.0 256.411 256.411 256.411 -0.000 -0.000 0.000 3/26364 60.0 201.674 201.674 201.674 -0.000 -0.000 0.000 3/26364 70.0 161.108 161.108 161.108 0.000 -0.000 0.000 3/26364 80.0 131.688 131.688 131.688 0.000 -0.000 0.000 3/26364 90.0 110.169 110.169 110.169 0.000 -0.000 0.000 3/26364 100.0 94.114 94.114 94.114 0.000 -0.000 0.000 3/26364 110.0 81.850 81.850 81.850 0.000 -0.000 0.000 3/26364 120.0 72.259 72.259 72.259 0.000 -0.000 0.000 3/26364 130.0 64.589 64.589 64.589 0.000 -0.000 0.000 3/26364 140.0 58.335 58.335 58.335 0.000 -0.000 0.000 3/26364 150.0 53.147 53.147 53.147 0.000 -0.000 0.000 3/26364 160.0 48.777 48.777 48.777 0.000 -0.000 0.000 3/26364 170.0 45.049 45.049 45.049 0.000 -0.000 0.000 3/26364 180.0 41.833 41.833 41.833 0.000 -0.000 0.000 3/26364 190.0 39.032 39.032 39.032 0.000 -0.000 0.000 3/26364 200.0 36.572 36.572 36.572 0.000 -0.000 0.000 3/26364 210.0 34.395 34.395 34.395 0.000 -0.000 0.000 3/26364 220.0 32.456 32.456 32.456 0.000 -0.000 0.000 3/26364 230.0 30.719 30.719 30.719 0.000 -0.000 0.000 3/26364 240.0 29.156 29.156 29.156 0.000 -0.000 0.000 3/26364 250.0 27.740 27.740 27.740 0.000 -0.000 0.000 3/26364 260.0 26.454 26.454 26.454 0.000 -0.000 0.000 3/26364 270.0 25.281 25.281 25.281 0.000 -0.000 0.000 3/26364 280.0 24.206 24.206 24.206 0.000 -0.000 0.000 3/26364 290.0 23.219 23.219 23.219 0.000 -0.000 0.000 3/26364 300.0 22.309 22.309 22.309 0.000 -0.000 0.000 3/26364 310.0 21.467 21.467 21.467 0.000 -0.000 0.000 3/26364 320.0 20.687 20.687 20.687 0.000 -0.000 0.000 3/26364 330.0 19.961 19.961 19.961 0.000 -0.000 0.000 3/26364 340.0 19.286 19.286 19.286 0.000 -0.000 0.000 3/26364 350.0 18.654 18.654 18.654 0.000 -0.000 0.000 3/26364 360.0 18.063 18.063 18.063 0.000 -0.000 0.000 3/26364 370.0 17.509 17.509 17.509 0.000 -0.000 0.000 3/26364 380.0 16.988 16.988 16.988 0.000 -0.000 0.000 3/26364 390.0 16.498 16.498 16.498 0.000 -0.000 0.000 3/26364 400.0 16.036 16.036 16.036 0.000 -0.000 0.000 3/26364 410.0 15.599 15.599 15.599 0.000 -0.000 0.000 3/26364 420.0 15.185 15.185 15.185 0.000 -0.000 0.000 3/26364 430.0 14.794 14.794 14.794 0.000 -0.000 0.000 3/26364 440.0 14.422 14.422 14.422 0.000 -0.000 0.000 3/26364 450.0 14.069 14.069 14.069 0.000 -0.000 0.000 3/26364 460.0 13.733 13.733 13.733 0.000 -0.000 0.000 3/26364 470.0 13.413 13.413 13.413 0.000 -0.000 0.000 3/26364 480.0 13.108 13.108 13.108 0.000 -0.000 0.000 3/26364 490.0 12.817 12.817 12.817 0.000 -0.000 0.000 3/26364 500.0 12.539 12.539 12.539 0.000 -0.000 0.000 3/26364 510.0 12.272 12.272 12.272 0.000 -0.000 0.000 3/26364 520.0 12.017 12.017 12.017 0.000 -0.000 0.000 3/26364 530.0 11.773 11.773 11.773 0.000 -0.000 0.000 3/26364 540.0 11.539 11.539 11.539 0.000 -0.000 0.000 3/26364 550.0 11.313 11.313 11.313 0.000 -0.000 0.000 3/26364 560.0 11.097 11.097 11.097 0.000 -0.000 0.000 3/26364 570.0 10.889 10.889 10.889 0.000 -0.000 0.000 3/26364 580.0 10.689 10.689 10.689 0.000 -0.000 0.000 3/26364 590.0 10.496 10.496 10.496 0.000 -0.000 0.000 3/26364 600.0 10.310 10.310 10.310 0.000 -0.000 0.000 3/26364 610.0 10.130 10.130 10.130 0.000 -0.000 0.000 3/26364 620.0 9.957 9.957 9.957 0.000 -0.000 0.000 3/26364 630.0 9.790 9.790 9.790 0.000 -0.000 0.000 3/26364 640.0 9.628 9.628 9.628 0.000 -0.000 0.000 3/26364 650.0 9.472 9.472 9.472 0.000 -0.000 0.000 3/26364 660.0 9.321 9.321 9.321 0.000 -0.000 0.000 3/26364 670.0 9.175 9.175 9.175 0.000 -0.000 0.000 3/26364 680.0 9.033 9.033 9.033 0.000 -0.000 0.000 3/26364 690.0 8.896 8.896 8.896 0.000 -0.000 0.000 3/26364 700.0 8.762 8.762 8.762 0.000 -0.000 0.000 3/26364 710.0 8.633 8.633 8.633 0.000 -0.000 0.000 3/26364 720.0 8.508 8.508 8.508 0.000 -0.000 0.000 3/26364 730.0 8.386 8.386 8.386 0.000 -0.000 0.000 3/26364 740.0 8.268 8.268 8.268 0.000 -0.000 0.000 3/26364 750.0 8.153 8.153 8.153 0.000 -0.000 0.000 3/26364 760.0 8.041 8.041 8.041 0.000 -0.000 0.000 3/26364 770.0 7.933 7.933 7.933 0.000 -0.000 0.000 3/26364 780.0 7.827 7.827 7.827 0.000 -0.000 0.000 3/26364 790.0 7.724 7.724 7.724 0.000 -0.000 0.000 3/26364 800.0 7.624 7.624 7.624 0.000 -0.000 0.000 3/26364 810.0 7.526 7.526 7.526 0.000 -0.000 0.000 3/26364 820.0 7.431 7.431 7.431 0.000 -0.000 0.000 3/26364 830.0 7.339 7.339 7.339 0.000 -0.000 0.000 3/26364 840.0 7.248 7.248 7.248 0.000 -0.000 0.000 3/26364 850.0 7.160 7.160 7.160 0.000 -0.000 0.000 3/26364 860.0 7.074 7.074 7.074 0.000 -0.000 0.000 3/26364 870.0 6.990 6.990 6.990 0.000 -0.000 0.000 3/26364 880.0 6.908 6.908 6.908 0.000 -0.000 0.000 3/26364 890.0 6.828 6.828 6.828 0.000 -0.000 0.000 3/26364 900.0 6.750 6.750 6.750 0.000 -0.000 0.000 3/26364 910.0 6.674 6.674 6.674 0.000 -0.000 0.000 3/26364 920.0 6.599 6.599 6.599 0.000 -0.000 0.000 3/26364 930.0 6.526 6.526 6.526 0.000 -0.000 0.000 3/26364 940.0 6.455 6.455 6.455 0.000 -0.000 0.000 3/26364 950.0 6.385 6.385 6.385 0.000 -0.000 0.000 3/26364 960.0 6.317 6.317 6.317 0.000 -0.000 0.000 3/26364 970.0 6.250 6.250 6.250 0.000 -0.000 0.000 3/26364 980.0 6.185 6.185 6.185 0.000 -0.000 0.000 3/26364 990.0 6.121 6.121 6.121 0.000 -0.000 0.000 3/26364 1000.0 6.058 6.058 6.058 0.000 -0.000 0.000 3/26364 Thermal conductivity related properties were written into "kappa-m131313.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 18:30:32]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|