----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 18:41:18]------------------------- Compiled with OpenMP support (max 48 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.882200730000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.882200730000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.882200730000000 Atomic positions (fractional): *1 Zn 0.50000000000000 0.00000000000000 0.50000000000000 65.380 2 Zn 0.00000000000000 0.50000000000000 0.50000000000000 65.380 3 Zn 0.50000000000000 0.50000000000000 0.00000000000000 65.380 *4 Nb 0.00000000000000 0.00000000000000 0.00000000000000 92.906 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.882200730000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.882200730000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.882200730000000 Atomic positions (fractional): *1 Zn 0.50000000000000 0.00000000000000 0.50000000000000 65.380 > 1 2 Zn 0.00000000000000 0.50000000000000 0.50000000000000 65.380 > 2 3 Zn 0.50000000000000 0.50000000000000 0.00000000000000 65.380 > 3 *4 Nb 0.00000000000000 0.00000000000000 0.00000000000000 92.906 > 4 -------------------------------- super cell -------------------------------- Lattice vectors: a 11.646602189999999 0.000000000000000 0.000000000000000 b 0.000000000000000 11.646602189999999 0.000000000000000 c 0.000000000000000 0.000000000000000 11.646602189999999 Atomic positions (fractional): *1 Zn 0.16666666666667 0.00000000000000 0.16666666666667 65.380 > 1 2 Zn 0.50000000000000 0.00000000000000 0.16666666666667 65.380 > 1 3 Zn 0.83333333333333 0.00000000000000 0.16666666666667 65.380 > 1 4 Zn 0.16666666666667 0.33333333333333 0.16666666666667 65.380 > 1 5 Zn 0.50000000000000 0.33333333333333 0.16666666666667 65.380 > 1 6 Zn 0.83333333333333 0.33333333333333 0.16666666666667 65.380 > 1 7 Zn 0.16666666666667 0.66666666666667 0.16666666666667 65.380 > 1 8 Zn 0.50000000000000 0.66666666666667 0.16666666666667 65.380 > 1 9 Zn 0.83333333333333 0.66666666666667 0.16666666666667 65.380 > 1 10 Zn 0.16666666666667 0.00000000000000 0.50000000000000 65.380 > 1 11 Zn 0.50000000000000 0.00000000000000 0.50000000000000 65.380 > 1 12 Zn 0.83333333333333 0.00000000000000 0.50000000000000 65.380 > 1 13 Zn 0.16666666666667 0.33333333333333 0.50000000000000 65.380 > 1 14 Zn 0.50000000000000 0.33333333333333 0.50000000000000 65.380 > 1 15 Zn 0.83333333333333 0.33333333333333 0.50000000000000 65.380 > 1 16 Zn 0.16666666666667 0.66666666666667 0.50000000000000 65.380 > 1 17 Zn 0.50000000000000 0.66666666666667 0.50000000000000 65.380 > 1 18 Zn 0.83333333333333 0.66666666666667 0.50000000000000 65.380 > 1 19 Zn 0.16666666666667 0.00000000000000 0.83333333333333 65.380 > 1 20 Zn 0.50000000000000 0.00000000000000 0.83333333333333 65.380 > 1 21 Zn 0.83333333333333 0.00000000000000 0.83333333333333 65.380 > 1 22 Zn 0.16666666666667 0.33333333333333 0.83333333333333 65.380 > 1 23 Zn 0.50000000000000 0.33333333333333 0.83333333333333 65.380 > 1 24 Zn 0.83333333333333 0.33333333333333 0.83333333333333 65.380 > 1 25 Zn 0.16666666666667 0.66666666666667 0.83333333333333 65.380 > 1 26 Zn 0.50000000000000 0.66666666666667 0.83333333333333 65.380 > 1 27 Zn 0.83333333333333 0.66666666666667 0.83333333333333 65.380 > 1 28 Zn 0.00000000000000 0.16666666666667 0.16666666666667 65.380 > 2 29 Zn 0.33333333333333 0.16666666666667 0.16666666666667 65.380 > 2 30 Zn 0.66666666666667 0.16666666666667 0.16666666666667 65.380 > 2 31 Zn 0.00000000000000 0.50000000000000 0.16666666666667 65.380 > 2 32 Zn 0.33333333333333 0.50000000000000 0.16666666666667 65.380 > 2 33 Zn 0.66666666666667 0.50000000000000 0.16666666666667 65.380 > 2 34 Zn 0.00000000000000 0.83333333333333 0.16666666666667 65.380 > 2 35 Zn 0.33333333333333 0.83333333333333 0.16666666666667 65.380 > 2 36 Zn 0.66666666666667 0.83333333333333 0.16666666666667 65.380 > 2 37 Zn 0.00000000000000 0.16666666666667 0.50000000000000 65.380 > 2 38 Zn 0.33333333333333 0.16666666666667 0.50000000000000 65.380 > 2 39 Zn 0.66666666666667 0.16666666666667 0.50000000000000 65.380 > 2 40 Zn 0.00000000000000 0.50000000000000 0.50000000000000 65.380 > 2 41 Zn 0.33333333333333 0.50000000000000 0.50000000000000 65.380 > 2 42 Zn 0.66666666666667 0.50000000000000 0.50000000000000 65.380 > 2 43 Zn 0.00000000000000 0.83333333333333 0.50000000000000 65.380 > 2 44 Zn 0.33333333333333 0.83333333333333 0.50000000000000 65.380 > 2 45 Zn 0.66666666666667 0.83333333333333 0.50000000000000 65.380 > 2 46 Zn 0.00000000000000 0.16666666666667 0.83333333333333 65.380 > 2 47 Zn 0.33333333333333 0.16666666666667 0.83333333333333 65.380 > 2 48 Zn 0.66666666666667 0.16666666666667 0.83333333333333 65.380 > 2 49 Zn 0.00000000000000 0.50000000000000 0.83333333333333 65.380 > 2 50 Zn 0.33333333333333 0.50000000000000 0.83333333333333 65.380 > 2 51 Zn 0.66666666666667 0.50000000000000 0.83333333333333 65.380 > 2 52 Zn 0.00000000000000 0.83333333333333 0.83333333333333 65.380 > 2 53 Zn 0.33333333333333 0.83333333333333 0.83333333333333 65.380 > 2 54 Zn 0.66666666666667 0.83333333333333 0.83333333333333 65.380 > 2 55 Zn 0.16666666666667 0.16666666666667 0.00000000000000 65.380 > 3 56 Zn 0.50000000000000 0.16666666666667 0.00000000000000 65.380 > 3 57 Zn 0.83333333333333 0.16666666666667 0.00000000000000 65.380 > 3 58 Zn 0.16666666666667 0.50000000000000 0.00000000000000 65.380 > 3 59 Zn 0.50000000000000 0.50000000000000 0.00000000000000 65.380 > 3 60 Zn 0.83333333333333 0.50000000000000 0.00000000000000 65.380 > 3 61 Zn 0.16666666666667 0.83333333333333 0.00000000000000 65.380 > 3 62 Zn 0.50000000000000 0.83333333333333 0.00000000000000 65.380 > 3 63 Zn 0.83333333333333 0.83333333333333 0.00000000000000 65.380 > 3 64 Zn 0.16666666666667 0.16666666666667 0.33333333333333 65.380 > 3 65 Zn 0.50000000000000 0.16666666666667 0.33333333333333 65.380 > 3 66 Zn 0.83333333333333 0.16666666666667 0.33333333333333 65.380 > 3 67 Zn 0.16666666666667 0.50000000000000 0.33333333333333 65.380 > 3 68 Zn 0.50000000000000 0.50000000000000 0.33333333333333 65.380 > 3 69 Zn 0.83333333333333 0.50000000000000 0.33333333333333 65.380 > 3 70 Zn 0.16666666666667 0.83333333333333 0.33333333333333 65.380 > 3 71 Zn 0.50000000000000 0.83333333333333 0.33333333333333 65.380 > 3 72 Zn 0.83333333333333 0.83333333333333 0.33333333333333 65.380 > 3 73 Zn 0.16666666666667 0.16666666666667 0.66666666666667 65.380 > 3 74 Zn 0.50000000000000 0.16666666666667 0.66666666666667 65.380 > 3 75 Zn 0.83333333333333 0.16666666666667 0.66666666666667 65.380 > 3 76 Zn 0.16666666666667 0.50000000000000 0.66666666666667 65.380 > 3 77 Zn 0.50000000000000 0.50000000000000 0.66666666666667 65.380 > 3 78 Zn 0.83333333333333 0.50000000000000 0.66666666666667 65.380 > 3 79 Zn 0.16666666666667 0.83333333333333 0.66666666666667 65.380 > 3 80 Zn 0.50000000000000 0.83333333333333 0.66666666666667 65.380 > 3 81 Zn 0.83333333333333 0.83333333333333 0.66666666666667 65.380 > 3 *82 Nb 0.00000000000000 0.00000000000000 0.00000000000000 92.906 > 4 83 Nb 0.33333333333333 0.00000000000000 0.00000000000000 92.906 > 4 84 Nb 0.66666666666667 0.00000000000000 0.00000000000000 92.906 > 4 85 Nb 0.00000000000000 0.33333333333333 0.00000000000000 92.906 > 4 86 Nb 0.33333333333333 0.33333333333333 0.00000000000000 92.906 > 4 87 Nb 0.66666666666667 0.33333333333333 0.00000000000000 92.906 > 4 88 Nb 0.00000000000000 0.66666666666667 0.00000000000000 92.906 > 4 89 Nb 0.33333333333333 0.66666666666667 0.00000000000000 92.906 > 4 90 Nb 0.66666666666667 0.66666666666667 0.00000000000000 92.906 > 4 91 Nb 0.00000000000000 0.00000000000000 0.33333333333333 92.906 > 4 92 Nb 0.33333333333333 0.00000000000000 0.33333333333333 92.906 > 4 93 Nb 0.66666666666667 0.00000000000000 0.33333333333333 92.906 > 4 94 Nb 0.00000000000000 0.33333333333333 0.33333333333333 92.906 > 4 95 Nb 0.33333333333333 0.33333333333333 0.33333333333333 92.906 > 4 96 Nb 0.66666666666667 0.33333333333333 0.33333333333333 92.906 > 4 97 Nb 0.00000000000000 0.66666666666667 0.33333333333333 92.906 > 4 98 Nb 0.33333333333333 0.66666666666667 0.33333333333333 92.906 > 4 99 Nb 0.66666666666667 0.66666666666667 0.33333333333333 92.906 > 4 100 Nb 0.00000000000000 0.00000000000000 0.66666666666667 92.906 > 4 101 Nb 0.33333333333333 0.00000000000000 0.66666666666667 92.906 > 4 102 Nb 0.66666666666667 0.00000000000000 0.66666666666667 92.906 > 4 103 Nb 0.00000000000000 0.33333333333333 0.66666666666667 92.906 > 4 104 Nb 0.33333333333333 0.33333333333333 0.66666666666667 92.906 > 4 105 Nb 0.66666666666667 0.33333333333333 0.66666666666667 92.906 > 4 106 Nb 0.00000000000000 0.66666666666667 0.66666666666667 92.906 > 4 107 Nb 0.33333333333333 0.66666666666667 0.66666666666667 92.906 > 4 108 Nb 0.66666666666667 0.66666666666667 0.66666666666667 92.906 > 4 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 142.4427050 0.0000000 0.0000000 0.0000000 142.4427050 0.0000000 0.0000000 0.0000000 142.4427050 -------------------------- Born effective charges -------------------------- 1 Zn 0.2920830 0.0000000 0.0000000 0.0000000 3.3230961 0.0000000 0.0000000 0.0000000 0.2920830 2 Zn 3.3230961 0.0000000 0.0000000 0.0000000 0.2920830 0.0000000 0.0000000 0.0000000 0.2920830 3 Zn 0.2920830 0.0000000 0.0000000 0.0000000 0.2920830 0.0000000 0.0000000 0.0000000 3.3230961 4 Nb -3.9072622 0.0000000 0.0000000 0.0000000 -3.9072622 0.0000000 0.0000000 0.0000000 -3.9072622 ---------------------------------------------------------------------------- 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.005 Solver_block: 80 / 80 - Time: 0.029 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.035 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 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 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 18:41:23]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:41:24]------------------------- Compiled with OpenMP support (max 48 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.882200730000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.882200730000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.882200730000000 Atomic positions (fractional): 1 Zn 0.50000000000000 0.00000000000000 0.50000000000000 65.380 2 Zn 0.00000000000000 0.50000000000000 0.50000000000000 65.380 3 Zn 0.50000000000000 0.50000000000000 0.00000000000000 65.380 4 Nb 0.00000000000000 0.00000000000000 0.00000000000000 92.906 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.646602189999999 0.000000000000000 0.000000000000000 b 0.000000000000000 11.646602189999999 0.000000000000000 c 0.000000000000000 0.000000000000000 11.646602189999999 Atomic positions (fractional): 1 Zn 0.16666666666667 0.00000000000000 0.16666666666667 65.380 > 1 2 Zn 0.50000000000000 0.00000000000000 0.16666666666667 65.380 > 1 3 Zn 0.83333333333333 0.00000000000000 0.16666666666667 65.380 > 1 4 Zn 0.16666666666667 0.33333333333333 0.16666666666667 65.380 > 1 5 Zn 0.50000000000000 0.33333333333333 0.16666666666667 65.380 > 1 6 Zn 0.83333333333333 0.33333333333333 0.16666666666667 65.380 > 1 7 Zn 0.16666666666667 0.66666666666667 0.16666666666667 65.380 > 1 8 Zn 0.50000000000000 0.66666666666667 0.16666666666667 65.380 > 1 9 Zn 0.83333333333333 0.66666666666667 0.16666666666667 65.380 > 1 10 Zn 0.16666666666667 0.00000000000000 0.50000000000000 65.380 > 1 11 Zn 0.50000000000000 0.00000000000000 0.50000000000000 65.380 > 1 12 Zn 0.83333333333333 0.00000000000000 0.50000000000000 65.380 > 1 13 Zn 0.16666666666667 0.33333333333333 0.50000000000000 65.380 > 1 14 Zn 0.50000000000000 0.33333333333333 0.50000000000000 65.380 > 1 15 Zn 0.83333333333333 0.33333333333333 0.50000000000000 65.380 > 1 16 Zn 0.16666666666667 0.66666666666667 0.50000000000000 65.380 > 1 17 Zn 0.50000000000000 0.66666666666667 0.50000000000000 65.380 > 1 18 Zn 0.83333333333333 0.66666666666667 0.50000000000000 65.380 > 1 19 Zn 0.16666666666667 0.00000000000000 0.83333333333333 65.380 > 1 20 Zn 0.50000000000000 0.00000000000000 0.83333333333333 65.380 > 1 21 Zn 0.83333333333333 0.00000000000000 0.83333333333333 65.380 > 1 22 Zn 0.16666666666667 0.33333333333333 0.83333333333333 65.380 > 1 23 Zn 0.50000000000000 0.33333333333333 0.83333333333333 65.380 > 1 24 Zn 0.83333333333333 0.33333333333333 0.83333333333333 65.380 > 1 25 Zn 0.16666666666667 0.66666666666667 0.83333333333333 65.380 > 1 26 Zn 0.50000000000000 0.66666666666667 0.83333333333333 65.380 > 1 27 Zn 0.83333333333333 0.66666666666667 0.83333333333333 65.380 > 1 28 Zn 0.00000000000000 0.16666666666667 0.16666666666667 65.380 > 28 29 Zn 0.33333333333333 0.16666666666667 0.16666666666667 65.380 > 28 30 Zn 0.66666666666667 0.16666666666667 0.16666666666667 65.380 > 28 31 Zn 0.00000000000000 0.50000000000000 0.16666666666667 65.380 > 28 32 Zn 0.33333333333333 0.50000000000000 0.16666666666667 65.380 > 28 33 Zn 0.66666666666667 0.50000000000000 0.16666666666667 65.380 > 28 34 Zn 0.00000000000000 0.83333333333333 0.16666666666667 65.380 > 28 35 Zn 0.33333333333333 0.83333333333333 0.16666666666667 65.380 > 28 36 Zn 0.66666666666667 0.83333333333333 0.16666666666667 65.380 > 28 37 Zn 0.00000000000000 0.16666666666667 0.50000000000000 65.380 > 28 38 Zn 0.33333333333333 0.16666666666667 0.50000000000000 65.380 > 28 39 Zn 0.66666666666667 0.16666666666667 0.50000000000000 65.380 > 28 40 Zn 0.00000000000000 0.50000000000000 0.50000000000000 65.380 > 28 41 Zn 0.33333333333333 0.50000000000000 0.50000000000000 65.380 > 28 42 Zn 0.66666666666667 0.50000000000000 0.50000000000000 65.380 > 28 43 Zn 0.00000000000000 0.83333333333333 0.50000000000000 65.380 > 28 44 Zn 0.33333333333333 0.83333333333333 0.50000000000000 65.380 > 28 45 Zn 0.66666666666667 0.83333333333333 0.50000000000000 65.380 > 28 46 Zn 0.00000000000000 0.16666666666667 0.83333333333333 65.380 > 28 47 Zn 0.33333333333333 0.16666666666667 0.83333333333333 65.380 > 28 48 Zn 0.66666666666667 0.16666666666667 0.83333333333333 65.380 > 28 49 Zn 0.00000000000000 0.50000000000000 0.83333333333333 65.380 > 28 50 Zn 0.33333333333333 0.50000000000000 0.83333333333333 65.380 > 28 51 Zn 0.66666666666667 0.50000000000000 0.83333333333333 65.380 > 28 52 Zn 0.00000000000000 0.83333333333333 0.83333333333333 65.380 > 28 53 Zn 0.33333333333333 0.83333333333333 0.83333333333333 65.380 > 28 54 Zn 0.66666666666667 0.83333333333333 0.83333333333333 65.380 > 28 55 Zn 0.16666666666667 0.16666666666667 0.00000000000000 65.380 > 55 56 Zn 0.50000000000000 0.16666666666667 0.00000000000000 65.380 > 55 57 Zn 0.83333333333333 0.16666666666667 0.00000000000000 65.380 > 55 58 Zn 0.16666666666667 0.50000000000000 0.00000000000000 65.380 > 55 59 Zn 0.50000000000000 0.50000000000000 0.00000000000000 65.380 > 55 60 Zn 0.83333333333333 0.50000000000000 0.00000000000000 65.380 > 55 61 Zn 0.16666666666667 0.83333333333333 0.00000000000000 65.380 > 55 62 Zn 0.50000000000000 0.83333333333333 0.00000000000000 65.380 > 55 63 Zn 0.83333333333333 0.83333333333333 0.00000000000000 65.380 > 55 64 Zn 0.16666666666667 0.16666666666667 0.33333333333333 65.380 > 55 65 Zn 0.50000000000000 0.16666666666667 0.33333333333333 65.380 > 55 66 Zn 0.83333333333333 0.16666666666667 0.33333333333333 65.380 > 55 67 Zn 0.16666666666667 0.50000000000000 0.33333333333333 65.380 > 55 68 Zn 0.50000000000000 0.50000000000000 0.33333333333333 65.380 > 55 69 Zn 0.83333333333333 0.50000000000000 0.33333333333333 65.380 > 55 70 Zn 0.16666666666667 0.83333333333333 0.33333333333333 65.380 > 55 71 Zn 0.50000000000000 0.83333333333333 0.33333333333333 65.380 > 55 72 Zn 0.83333333333333 0.83333333333333 0.33333333333333 65.380 > 55 73 Zn 0.16666666666667 0.16666666666667 0.66666666666667 65.380 > 55 74 Zn 0.50000000000000 0.16666666666667 0.66666666666667 65.380 > 55 75 Zn 0.83333333333333 0.16666666666667 0.66666666666667 65.380 > 55 76 Zn 0.16666666666667 0.50000000000000 0.66666666666667 65.380 > 55 77 Zn 0.50000000000000 0.50000000000000 0.66666666666667 65.380 > 55 78 Zn 0.83333333333333 0.50000000000000 0.66666666666667 65.380 > 55 79 Zn 0.16666666666667 0.83333333333333 0.66666666666667 65.380 > 55 80 Zn 0.50000000000000 0.83333333333333 0.66666666666667 65.380 > 55 81 Zn 0.83333333333333 0.83333333333333 0.66666666666667 65.380 > 55 82 Nb 0.00000000000000 0.00000000000000 0.00000000000000 92.906 > 82 83 Nb 0.33333333333333 0.00000000000000 0.00000000000000 92.906 > 82 84 Nb 0.66666666666667 0.00000000000000 0.00000000000000 92.906 > 82 85 Nb 0.00000000000000 0.33333333333333 0.00000000000000 92.906 > 82 86 Nb 0.33333333333333 0.33333333333333 0.00000000000000 92.906 > 82 87 Nb 0.66666666666667 0.33333333333333 0.00000000000000 92.906 > 82 88 Nb 0.00000000000000 0.66666666666667 0.00000000000000 92.906 > 82 89 Nb 0.33333333333333 0.66666666666667 0.00000000000000 92.906 > 82 90 Nb 0.66666666666667 0.66666666666667 0.00000000000000 92.906 > 82 91 Nb 0.00000000000000 0.00000000000000 0.33333333333333 92.906 > 82 92 Nb 0.33333333333333 0.00000000000000 0.33333333333333 92.906 > 82 93 Nb 0.66666666666667 0.00000000000000 0.33333333333333 92.906 > 82 94 Nb 0.00000000000000 0.33333333333333 0.33333333333333 92.906 > 82 95 Nb 0.33333333333333 0.33333333333333 0.33333333333333 92.906 > 82 96 Nb 0.66666666666667 0.33333333333333 0.33333333333333 92.906 > 82 97 Nb 0.00000000000000 0.66666666666667 0.33333333333333 92.906 > 82 98 Nb 0.33333333333333 0.66666666666667 0.33333333333333 92.906 > 82 99 Nb 0.66666666666667 0.66666666666667 0.33333333333333 92.906 > 82 100 Nb 0.00000000000000 0.00000000000000 0.66666666666667 92.906 > 82 101 Nb 0.33333333333333 0.00000000000000 0.66666666666667 92.906 > 82 102 Nb 0.66666666666667 0.00000000000000 0.66666666666667 92.906 > 82 103 Nb 0.00000000000000 0.33333333333333 0.66666666666667 92.906 > 82 104 Nb 0.33333333333333 0.33333333333333 0.66666666666667 92.906 > 82 105 Nb 0.66666666666667 0.33333333333333 0.66666666666667 92.906 > 82 106 Nb 0.00000000000000 0.66666666666667 0.66666666666667 92.906 > 82 107 Nb 0.33333333333333 0.66666666666667 0.66666666666667 92.906 > 82 108 Nb 0.66666666666667 0.66666666666667 0.66666666666667 92.906 > 82 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 142.4427050 0.0000000 0.0000000 0.0000000 142.4427050 0.0000000 0.0000000 0.0000000 142.4427050 -------------------------- Born effective charges -------------------------- 1 Zn 0.2920830 0.0000000 0.0000000 0.0000000 3.3230961 0.0000000 0.0000000 0.0000000 0.2920830 2 Zn 3.3230961 0.0000000 0.0000000 0.0000000 0.2920830 0.0000000 0.0000000 0.0000000 0.2920830 3 Zn 0.2920830 0.0000000 0.0000000 0.0000000 0.2920830 0.0000000 0.0000000 0.0000000 3.3230961 4 Nb -3.9072622 0.0000000 0.0000000 0.0000000 -3.9072622 0.0000000 0.0000000 0.0000000 -3.9072622 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 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 (yyx) -0.00000000 (yyx) -0.00000000 (yxy) 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:41:26]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:41:26]------------------------- Compiled with OpenMP support (max 48 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.882200730000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.882200730000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.882200730000000 Atomic positions (fractional): 1 Zn 0.50000000000000 0.00000000000000 0.50000000000000 65.380 2 Zn 0.00000000000000 0.50000000000000 0.50000000000000 65.380 3 Zn 0.50000000000000 0.50000000000000 0.00000000000000 65.380 4 Nb 0.00000000000000 0.00000000000000 0.00000000000000 92.906 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.646602189999999 0.000000000000000 0.000000000000000 b 0.000000000000000 11.646602189999999 0.000000000000000 c 0.000000000000000 0.000000000000000 11.646602189999999 Atomic positions (fractional): 1 Zn 0.16666666666667 0.00000000000000 0.16666666666667 65.380 > 1 2 Zn 0.50000000000000 0.00000000000000 0.16666666666667 65.380 > 1 3 Zn 0.83333333333333 0.00000000000000 0.16666666666667 65.380 > 1 4 Zn 0.16666666666667 0.33333333333333 0.16666666666667 65.380 > 1 5 Zn 0.50000000000000 0.33333333333333 0.16666666666667 65.380 > 1 6 Zn 0.83333333333333 0.33333333333333 0.16666666666667 65.380 > 1 7 Zn 0.16666666666667 0.66666666666667 0.16666666666667 65.380 > 1 8 Zn 0.50000000000000 0.66666666666667 0.16666666666667 65.380 > 1 9 Zn 0.83333333333333 0.66666666666667 0.16666666666667 65.380 > 1 10 Zn 0.16666666666667 0.00000000000000 0.50000000000000 65.380 > 1 11 Zn 0.50000000000000 0.00000000000000 0.50000000000000 65.380 > 1 12 Zn 0.83333333333333 0.00000000000000 0.50000000000000 65.380 > 1 13 Zn 0.16666666666667 0.33333333333333 0.50000000000000 65.380 > 1 14 Zn 0.50000000000000 0.33333333333333 0.50000000000000 65.380 > 1 15 Zn 0.83333333333333 0.33333333333333 0.50000000000000 65.380 > 1 16 Zn 0.16666666666667 0.66666666666667 0.50000000000000 65.380 > 1 17 Zn 0.50000000000000 0.66666666666667 0.50000000000000 65.380 > 1 18 Zn 0.83333333333333 0.66666666666667 0.50000000000000 65.380 > 1 19 Zn 0.16666666666667 0.00000000000000 0.83333333333333 65.380 > 1 20 Zn 0.50000000000000 0.00000000000000 0.83333333333333 65.380 > 1 21 Zn 0.83333333333333 0.00000000000000 0.83333333333333 65.380 > 1 22 Zn 0.16666666666667 0.33333333333333 0.83333333333333 65.380 > 1 23 Zn 0.50000000000000 0.33333333333333 0.83333333333333 65.380 > 1 24 Zn 0.83333333333333 0.33333333333333 0.83333333333333 65.380 > 1 25 Zn 0.16666666666667 0.66666666666667 0.83333333333333 65.380 > 1 26 Zn 0.50000000000000 0.66666666666667 0.83333333333333 65.380 > 1 27 Zn 0.83333333333333 0.66666666666667 0.83333333333333 65.380 > 1 28 Zn 0.00000000000000 0.16666666666667 0.16666666666667 65.380 > 28 29 Zn 0.33333333333333 0.16666666666667 0.16666666666667 65.380 > 28 30 Zn 0.66666666666667 0.16666666666667 0.16666666666667 65.380 > 28 31 Zn 0.00000000000000 0.50000000000000 0.16666666666667 65.380 > 28 32 Zn 0.33333333333333 0.50000000000000 0.16666666666667 65.380 > 28 33 Zn 0.66666666666667 0.50000000000000 0.16666666666667 65.380 > 28 34 Zn 0.00000000000000 0.83333333333333 0.16666666666667 65.380 > 28 35 Zn 0.33333333333333 0.83333333333333 0.16666666666667 65.380 > 28 36 Zn 0.66666666666667 0.83333333333333 0.16666666666667 65.380 > 28 37 Zn 0.00000000000000 0.16666666666667 0.50000000000000 65.380 > 28 38 Zn 0.33333333333333 0.16666666666667 0.50000000000000 65.380 > 28 39 Zn 0.66666666666667 0.16666666666667 0.50000000000000 65.380 > 28 40 Zn 0.00000000000000 0.50000000000000 0.50000000000000 65.380 > 28 41 Zn 0.33333333333333 0.50000000000000 0.50000000000000 65.380 > 28 42 Zn 0.66666666666667 0.50000000000000 0.50000000000000 65.380 > 28 43 Zn 0.00000000000000 0.83333333333333 0.50000000000000 65.380 > 28 44 Zn 0.33333333333333 0.83333333333333 0.50000000000000 65.380 > 28 45 Zn 0.66666666666667 0.83333333333333 0.50000000000000 65.380 > 28 46 Zn 0.00000000000000 0.16666666666667 0.83333333333333 65.380 > 28 47 Zn 0.33333333333333 0.16666666666667 0.83333333333333 65.380 > 28 48 Zn 0.66666666666667 0.16666666666667 0.83333333333333 65.380 > 28 49 Zn 0.00000000000000 0.50000000000000 0.83333333333333 65.380 > 28 50 Zn 0.33333333333333 0.50000000000000 0.83333333333333 65.380 > 28 51 Zn 0.66666666666667 0.50000000000000 0.83333333333333 65.380 > 28 52 Zn 0.00000000000000 0.83333333333333 0.83333333333333 65.380 > 28 53 Zn 0.33333333333333 0.83333333333333 0.83333333333333 65.380 > 28 54 Zn 0.66666666666667 0.83333333333333 0.83333333333333 65.380 > 28 55 Zn 0.16666666666667 0.16666666666667 0.00000000000000 65.380 > 55 56 Zn 0.50000000000000 0.16666666666667 0.00000000000000 65.380 > 55 57 Zn 0.83333333333333 0.16666666666667 0.00000000000000 65.380 > 55 58 Zn 0.16666666666667 0.50000000000000 0.00000000000000 65.380 > 55 59 Zn 0.50000000000000 0.50000000000000 0.00000000000000 65.380 > 55 60 Zn 0.83333333333333 0.50000000000000 0.00000000000000 65.380 > 55 61 Zn 0.16666666666667 0.83333333333333 0.00000000000000 65.380 > 55 62 Zn 0.50000000000000 0.83333333333333 0.00000000000000 65.380 > 55 63 Zn 0.83333333333333 0.83333333333333 0.00000000000000 65.380 > 55 64 Zn 0.16666666666667 0.16666666666667 0.33333333333333 65.380 > 55 65 Zn 0.50000000000000 0.16666666666667 0.33333333333333 65.380 > 55 66 Zn 0.83333333333333 0.16666666666667 0.33333333333333 65.380 > 55 67 Zn 0.16666666666667 0.50000000000000 0.33333333333333 65.380 > 55 68 Zn 0.50000000000000 0.50000000000000 0.33333333333333 65.380 > 55 69 Zn 0.83333333333333 0.50000000000000 0.33333333333333 65.380 > 55 70 Zn 0.16666666666667 0.83333333333333 0.33333333333333 65.380 > 55 71 Zn 0.50000000000000 0.83333333333333 0.33333333333333 65.380 > 55 72 Zn 0.83333333333333 0.83333333333333 0.33333333333333 65.380 > 55 73 Zn 0.16666666666667 0.16666666666667 0.66666666666667 65.380 > 55 74 Zn 0.50000000000000 0.16666666666667 0.66666666666667 65.380 > 55 75 Zn 0.83333333333333 0.16666666666667 0.66666666666667 65.380 > 55 76 Zn 0.16666666666667 0.50000000000000 0.66666666666667 65.380 > 55 77 Zn 0.50000000000000 0.50000000000000 0.66666666666667 65.380 > 55 78 Zn 0.83333333333333 0.50000000000000 0.66666666666667 65.380 > 55 79 Zn 0.16666666666667 0.83333333333333 0.66666666666667 65.380 > 55 80 Zn 0.50000000000000 0.83333333333333 0.66666666666667 65.380 > 55 81 Zn 0.83333333333333 0.83333333333333 0.66666666666667 65.380 > 55 82 Nb 0.00000000000000 0.00000000000000 0.00000000000000 92.906 > 82 83 Nb 0.33333333333333 0.00000000000000 0.00000000000000 92.906 > 82 84 Nb 0.66666666666667 0.00000000000000 0.00000000000000 92.906 > 82 85 Nb 0.00000000000000 0.33333333333333 0.00000000000000 92.906 > 82 86 Nb 0.33333333333333 0.33333333333333 0.00000000000000 92.906 > 82 87 Nb 0.66666666666667 0.33333333333333 0.00000000000000 92.906 > 82 88 Nb 0.00000000000000 0.66666666666667 0.00000000000000 92.906 > 82 89 Nb 0.33333333333333 0.66666666666667 0.00000000000000 92.906 > 82 90 Nb 0.66666666666667 0.66666666666667 0.00000000000000 92.906 > 82 91 Nb 0.00000000000000 0.00000000000000 0.33333333333333 92.906 > 82 92 Nb 0.33333333333333 0.00000000000000 0.33333333333333 92.906 > 82 93 Nb 0.66666666666667 0.00000000000000 0.33333333333333 92.906 > 82 94 Nb 0.00000000000000 0.33333333333333 0.33333333333333 92.906 > 82 95 Nb 0.33333333333333 0.33333333333333 0.33333333333333 92.906 > 82 96 Nb 0.66666666666667 0.33333333333333 0.33333333333333 92.906 > 82 97 Nb 0.00000000000000 0.66666666666667 0.33333333333333 92.906 > 82 98 Nb 0.33333333333333 0.66666666666667 0.33333333333333 92.906 > 82 99 Nb 0.66666666666667 0.66666666666667 0.33333333333333 92.906 > 82 100 Nb 0.00000000000000 0.00000000000000 0.66666666666667 92.906 > 82 101 Nb 0.33333333333333 0.00000000000000 0.66666666666667 92.906 > 82 102 Nb 0.66666666666667 0.00000000000000 0.66666666666667 92.906 > 82 103 Nb 0.00000000000000 0.33333333333333 0.66666666666667 92.906 > 82 104 Nb 0.33333333333333 0.33333333333333 0.66666666666667 92.906 > 82 105 Nb 0.66666666666667 0.33333333333333 0.66666666666667 92.906 > 82 106 Nb 0.00000000000000 0.66666666666667 0.66666666666667 92.906 > 82 107 Nb 0.33333333333333 0.66666666666667 0.66666666666667 92.906 > 82 108 Nb 0.66666666666667 0.66666666666667 0.66666666666667 92.906 > 82 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 142.4427050 0.0000000 0.0000000 0.0000000 142.4427050 0.0000000 0.0000000 0.0000000 142.4427050 -------------------------- Born effective charges -------------------------- 1 Zn 0.2920830 0.0000000 0.0000000 0.0000000 3.3230961 0.0000000 0.0000000 0.0000000 0.2920830 2 Zn 3.3230961 0.0000000 0.0000000 0.0000000 0.2920830 0.0000000 0.0000000 0.0000000 0.2920830 3 Zn 0.2920830 0.0000000 0.0000000 0.0000000 0.2920830 0.0000000 0.0000000 0.0000000 3.3230961 4 Nb -3.9072622 0.0000000 0.0000000 0.0000000 -3.9072622 0.0000000 0.0000000 0.0000000 -3.9072622 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (yyx) -0.00000000 (yyx) -0.00000000 (yxy) Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 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.07, Number of G-points: 305, Lambda: 1.33 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: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 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 4.430 ( 0.000 0.000 0.000) 0.000 4.430 ( 0.000 0.000 0.000) 0.000 4.430 ( 0.000 0.000 0.000) 0.000 4.673 ( 0.000 0.000 0.000) 0.000 4.673 ( 0.000 0.000 0.000) 0.000 4.673 ( 0.000 0.000 0.000) 0.000 7.531 ( 0.000 0.000 0.000) 0.000 7.531 ( 0.000 0.000 0.000) 0.000 7.531 ( 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: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.648 ( 32.535 0.000 0.000) 32.535 0.648 ( 32.535 0.000 0.000) 32.535 0.996 ( 49.552 0.000 0.000) 49.552 4.437 ( 0.476 0.000 0.000) 0.476 4.437 ( 0.476 0.000 0.000) 0.476 4.557 ( -1.957 0.000 0.000) 1.957 4.703 ( 3.161 0.000 0.000) 3.161 4.703 ( 3.161 0.000 0.000) 3.161 4.712 ( 3.871 0.000 0.000) 3.871 7.488 ( -7.380 0.000 0.000) 7.380 7.509 ( -2.187 0.000 0.000) 2.187 7.509 ( -2.187 0.000 0.000) 2.187 ======================= Grid point 2 (3/84) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.286 ( 31.731 0.000 0.000) 31.731 1.286 ( 31.731 0.000 0.000) 31.731 1.950 ( 46.454 0.000 0.000) 46.454 4.435 ( -0.938 0.000 0.000) 0.938 4.435 ( -0.938 0.000 0.000) 0.938 4.500 ( -3.742 0.000 0.000) 3.742 4.802 ( 6.726 0.000 0.000) 6.726 4.802 ( 6.726 0.000 0.000) 6.726 4.820 ( 6.702 0.000 0.000) 6.702 7.276 ( -13.746 0.000 0.000) 13.746 7.446 ( -4.111 0.000 0.000) 4.111 7.446 ( -4.111 0.000 0.000) 4.111 ======================= Grid point 3 (4/84) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.901 ( 30.253 0.000 0.000) 30.253 1.901 ( 30.253 0.000 0.000) 30.253 2.825 ( 41.526 0.000 0.000) 41.526 4.393 ( -3.497 0.000 0.000) 3.497 4.393 ( -3.497 0.000 0.000) 3.497 4.414 ( -4.780 0.000 0.000) 4.780 4.961 ( 8.995 0.000 0.000) 8.995 4.961 ( 8.995 0.000 0.000) 8.995 4.967 ( 7.845 0.000 0.000) 7.845 6.954 ( -18.389 0.000 0.000) 18.389 7.350 ( -5.516 0.000 0.000) 5.516 7.350 ( -5.516 0.000 0.000) 5.516 ======================= Grid point 4 (5/84) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.478 ( 27.796 0.000 0.000) 27.796 2.478 ( 27.796 0.000 0.000) 27.796 3.576 ( 33.145 0.000 0.000) 33.145 4.293 ( -6.674 0.000 0.000) 6.674 4.293 ( -6.674 0.000 0.000) 6.674 4.332 ( -2.247 0.000 0.000) 2.247 5.118 ( 7.170 0.000 0.000) 7.170 5.144 ( 9.133 0.000 0.000) 9.133 5.144 ( 9.133 0.000 0.000) 9.133 6.560 ( -21.133 0.000 0.000) 21.133 7.233 ( -6.037 0.000 0.000) 6.037 7.233 ( -6.037 0.000 0.000) 6.037 ======================= Grid point 5 (6/84) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.993 ( 23.842 0.000 0.000) 23.842 2.993 ( 23.842 0.000 0.000) 23.842 3.996 ( 5.925 0.000 0.000) 5.925 4.126 ( -10.131 0.000 0.000) 10.131 4.126 ( -10.131 0.000 0.000) 10.131 4.475 ( 20.156 0.000 0.000) 20.156 5.241 ( 4.966 0.000 0.000) 4.966 5.308 ( 7.021 0.000 0.000) 7.021 5.308 ( 7.021 0.000 0.000) 7.021 6.129 ( -22.032 0.000 0.000) 22.032 7.120 ( -5.048 0.000 0.000) 5.048 7.120 ( -5.048 0.000 0.000) 5.048 ======================= Grid point 6 (7/84) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.397 ( 15.076 0.000 0.000) 15.076 3.397 ( 15.076 0.000 0.000) 15.076 3.908 ( -10.300 0.000 0.000) 10.300 3.908 ( -10.300 0.000 0.000) 10.300 4.008 ( -0.854 0.000 0.000) 0.854 4.916 ( 18.284 0.000 0.000) 18.284 5.309 ( 1.769 0.000 0.000) 1.769 5.408 ( 2.703 0.000 0.000) 2.703 5.408 ( 2.703 0.000 0.000) 2.703 5.722 ( -16.447 0.000 0.000) 16.447 7.047 ( -2.047 0.000 0.000) 2.047 7.047 ( -2.047 0.000 0.000) 2.047 ======================= Grid point 14 (8/84) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.820 ( 20.319 20.319 0.000) 28.735 0.919 ( 23.176 23.176 0.000) 32.776 1.460 ( 36.119 36.119 0.000) 51.079 4.387 ( -2.237 -2.237 0.000) 3.164 4.449 ( 0.886 0.886 0.000) 1.253 4.584 ( -0.073 -0.073 0.000) 0.104 4.721 ( 2.227 2.227 0.000) 3.150 4.728 ( 2.595 2.595 0.000) 3.669 4.767 ( 4.825 4.825 0.000) 6.823 7.408 ( -7.580 -7.580 0.000) 10.720 7.488 ( -2.101 -2.101 0.000) 2.971 7.493 ( -1.929 -1.929 0.000) 2.728 ======================= Grid point 15 (9/84) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.330 ( 28.585 7.014 0.000) 29.433 1.445 ( 28.372 15.024 0.000) 32.105 2.238 ( 40.493 25.156 0.000) 47.671 4.347 ( -2.337 -5.388 0.000) 5.873 4.457 ( -0.460 1.920 0.000) 1.975 4.550 ( -3.020 1.792 0.000) 3.512 4.797 ( 5.334 0.316 0.000) 5.344 4.813 ( 6.056 1.214 0.000) 6.176 4.891 ( 7.508 5.924 0.000) 9.564 7.198 ( -13.407 -7.553 0.000) 15.388 7.427 ( -3.973 -1.877 0.000) 4.395 7.432 ( -4.215 -1.664 0.000) 4.531 ======================= Grid point 16 (10/84) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.904 ( 28.758 1.631 0.000) 28.804 2.012 ( 28.466 10.757 0.000) 30.431 3.021 ( 37.841 17.897 0.000) 41.860 4.284 ( -4.032 -6.691 0.000) 7.812 4.422 ( -3.175 2.820 0.000) 4.247 4.470 ( -4.985 1.696 0.000) 5.265 4.925 ( 7.270 -1.528 0.000) 7.428 4.961 ( 8.614 0.072 0.000) 8.615 5.055 ( 8.878 6.442 0.000) 10.969 6.888 ( -17.436 -6.439 0.000) 18.587 7.331 ( -5.835 -1.908 0.000) 6.139 7.334 ( -5.379 -1.604 0.000) 5.613 ======================= Grid point 17 (11/84) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.456 ( 26.740 -1.374 0.000) 26.775 2.560 ( 26.490 7.927 0.000) 27.651 3.703 ( 29.243 11.519 0.000) 31.430 4.198 ( -3.920 -6.729 0.000) 7.788 4.328 ( -6.328 3.518 0.000) 7.240 4.360 ( -5.192 0.534 0.000) 5.219 5.072 ( 7.292 -2.894 0.000) 7.846 5.140 ( 9.009 -0.451 0.000) 9.020 5.232 ( 8.746 6.806 0.000) 11.083 6.521 ( -19.361 -3.864 0.000) 19.743 7.207 ( -6.495 -2.545 0.000) 6.976 7.220 ( -5.941 -1.374 0.000) 6.098 ======================= Grid point 18 (12/84) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.952 ( 22.976 -3.542 0.000) 23.248 3.050 ( 22.602 5.506 0.000) 23.263 3.988 ( -0.573 -2.270 0.000) 2.342 4.109 ( -8.001 -1.508 0.000) 8.142 4.170 ( -9.493 4.480 0.000) 10.497 4.476 ( 19.818 1.123 0.000) 19.849 5.200 ( 5.269 -3.505 0.000) 6.329 5.303 ( 7.042 -0.544 0.000) 7.063 5.395 ( 7.512 7.553 0.000) 10.652 6.134 ( -19.330 0.510 0.000) 19.337 7.084 ( -5.522 -3.413 0.000) 6.491 7.108 ( -5.008 -1.243 0.000) 5.160 ======================= Grid point 19 (13/84) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.337 ( 13.900 -5.547 0.000) 14.966 3.425 ( 13.376 2.724 0.000) 13.650 3.883 ( -8.318 -2.458 0.000) 8.674 3.972 ( -8.812 6.477 0.000) 10.936 3.974 ( -2.631 -3.909 0.000) 4.712 4.845 ( 12.676 -4.593 0.000) 13.482 5.271 ( 1.798 -3.770 0.000) 4.176 5.403 ( 2.734 -0.473 0.000) 2.775 5.538 ( 7.812 10.107 0.000) 12.774 5.777 ( -15.663 5.860 0.000) 16.723 7.003 ( -2.270 -4.134 0.000) 4.716 7.035 ( -2.042 -1.206 0.000) 2.371 ======================= Grid point 28 (14/84) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.591 ( 18.290 18.290 0.000) 25.866 1.828 ( 22.490 22.490 0.000) 31.805 2.823 ( 31.932 31.932 0.000) 45.158 4.247 ( -4.843 -4.843 0.000) 6.849 4.496 ( 1.466 1.466 0.000) 2.073 4.549 ( -1.961 -1.961 0.000) 2.773 4.823 ( 2.453 2.453 0.000) 3.468 4.859 ( 3.747 3.747 0.000) 5.299 5.034 ( 8.287 8.287 0.000) 11.720 6.993 ( -12.645 -12.645 0.000) 17.883 7.373 ( -3.617 -3.617 0.000) 5.115 7.375 ( -4.047 -4.047 0.000) 5.723 ======================= Grid point 29 (15/84) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.006 ( 22.624 8.912 0.000) 24.316 2.299 ( 24.364 17.410 0.000) 29.945 3.450 ( 30.329 23.618 0.000) 38.440 4.147 ( -5.330 -6.994 0.000) 8.793 4.473 ( -5.545 -1.424 0.000) 5.725 4.496 ( -1.633 4.348 0.000) 4.645 4.895 ( 4.726 -1.375 0.000) 4.922 4.968 ( 7.095 0.759 0.000) 7.135 5.211 ( 9.393 8.993 0.000) 13.004 6.714 ( -15.040 -10.635 0.000) 18.420 7.274 ( -6.057 -3.818 0.000) 7.160 7.286 ( -5.020 -3.176 0.000) 5.940 ======================= Grid point 30 (16/84) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.458 ( 22.503 2.126 0.000) 22.603 2.773 ( 23.067 13.069 0.000) 26.512 3.939 ( 10.041 6.983 0.000) 12.230 4.072 ( 3.812 -2.925 0.000) 4.806 4.356 ( -3.568 0.157 0.000) 3.572 4.430 ( -4.999 6.629 0.000) 8.303 5.000 ( 5.500 -4.550 0.000) 7.138 5.127 ( 8.504 -0.837 0.000) 8.545 5.400 ( 9.469 9.802 0.000) 13.629 6.416 ( -14.642 -6.310 0.000) 15.944 7.141 ( -7.159 -3.859 0.000) 8.133 7.178 ( -5.688 -2.800 0.000) 6.340 ======================= Grid point 31 (17/84) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.881 ( 19.663 -3.060 0.000) 19.899 3.197 ( 19.137 8.865 0.000) 21.090 3.883 ( -5.317 -7.618 0.000) 9.290 4.049 ( -9.540 -5.042 0.000) 10.791 4.303 ( -7.489 8.902 0.000) 11.633 4.546 ( 16.152 5.546 0.000) 17.077 5.099 ( 4.184 -6.794 0.000) 7.979 5.286 ( 7.070 -1.188 0.000) 7.169 5.579 ( 8.400 10.429 0.000) 13.391 6.149 ( -11.882 1.174 0.000) 11.940 7.001 ( -6.539 -4.453 0.000) 7.911 7.070 ( -4.906 -2.590 0.000) 5.547 ======================= Grid point 32 (18/84) ======================= q-point: ( 0.46 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.203 ( 10.998 -7.244 0.000) 13.169 3.497 ( 9.574 4.408 0.000) 10.540 3.801 ( -2.695 -6.766 0.000) 7.282 3.854 ( -7.578 -6.986 0.000) 10.307 4.160 ( -5.543 12.189 0.000) 13.391 4.790 ( 7.001 -0.953 0.000) 7.066 5.154 ( 1.297 -8.138 0.000) 8.241 5.388 ( 2.822 -1.082 0.000) 3.022 5.718 ( 4.907 7.421 0.000) 8.897 5.963 ( -6.185 12.310 0.000) 13.776 6.902 ( -2.886 -5.428 0.000) 6.148 6.998 ( -2.031 -2.537 0.000) 3.250 ======================= Grid point 42 (19/84) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.242 ( 14.237 14.237 0.000) 20.134 2.672 ( 19.587 19.587 0.000) 27.701 3.896 ( 19.479 19.479 0.000) 27.547 4.007 ( -7.137 -7.137 0.000) 10.093 4.416 ( -3.621 -3.621 0.000) 5.120 4.572 ( 2.598 2.598 0.000) 3.675 4.879 ( -0.059 -0.059 0.000) 0.083 5.005 ( 3.439 3.439 0.000) 4.864 5.404 ( 10.259 10.259 0.000) 14.509 6.491 ( -11.251 -11.251 0.000) 15.911 7.181 ( -5.483 -5.483 0.000) 7.754 7.209 ( -4.569 -4.569 0.000) 6.462 ======================= Grid point 43 (20/84) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.545 ( 15.765 6.408 0.000) 17.017 3.054 ( 18.423 14.654 0.000) 23.540 3.859 ( -7.581 -7.819 0.000) 10.891 4.081 ( -3.691 -4.214 0.000) 5.602 4.449 ( 9.655 9.293 0.000) 13.400 4.581 ( -1.713 8.179 0.000) 8.357 4.896 ( 1.685 -5.595 0.000) 5.843 5.111 ( 6.941 -0.537 0.000) 6.962 5.609 ( 10.125 10.950 0.000) 14.914 6.293 ( -8.115 -5.371 0.000) 9.731 7.060 ( -6.673 -4.284 0.000) 7.930 7.109 ( -5.374 -4.185 0.000) 6.811 ======================= Grid point 44 (21/84) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.846 ( 14.131 -0.247 0.000) 14.133 3.384 ( 14.360 9.646 0.000) 17.299 3.714 ( -6.687 -8.955 0.000) 11.176 3.907 ( -10.760 -9.181 0.000) 14.145 4.517 ( -4.244 12.320 0.000) 13.030 4.658 ( 9.109 4.954 0.000) 10.369 4.933 ( 1.688 -9.608 0.000) 9.755 5.256 ( 6.957 -1.829 0.000) 7.193 5.779 ( 6.155 8.932 0.000) 10.847 6.194 ( -1.384 3.849 0.000) 4.090 6.923 ( -6.807 -3.227 0.000) 7.533 7.005 ( -4.794 -3.988 0.000) 6.236 ======================= Grid point 45 (22/84) ======================= q-point: ( 0.46 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.073 ( 7.322 -5.464 0.000) 9.136 3.595 ( 6.046 5.344 0.000) 8.069 3.615 ( -2.654 -10.414 0.000) 10.747 3.720 ( -6.534 -7.432 0.000) 9.896 4.440 ( -2.647 15.597 0.000) 15.820 4.786 ( 3.504 -0.129 0.000) 3.506 4.952 ( 0.247 -11.878 0.000) 11.880 5.359 ( 2.946 -1.800 0.000) 3.453 5.841 ( 0.870 5.239 0.000) 5.311 6.221 ( 2.323 12.968 0.000) 13.174 6.811 ( -3.625 -3.350 0.000) 4.936 6.934 ( -2.031 -3.967 0.000) 4.456 ======================= Grid point 56 (23/84) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.699 ( 8.676 8.676 0.000) 12.270 3.335 ( 13.342 13.342 0.000) 18.868 3.700 ( -8.036 -8.036 0.000) 11.365 3.915 ( -11.074 -11.074 0.000) 15.661 4.629 ( 7.464 7.464 0.000) 10.556 4.725 ( 5.365 5.365 0.000) 7.587 4.795 ( -4.186 -4.186 0.000) 5.919 5.124 ( 2.529 2.529 0.000) 3.576 5.809 ( 8.238 8.238 0.000) 11.650 6.234 ( 0.236 0.236 0.000) 0.333 6.971 ( -4.693 -4.693 0.000) 6.638 7.015 ( -5.127 -5.127 0.000) 7.251 ======================= Grid point 57 (24/84) ======================= q-point: ( 0.38 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.871 ( 8.186 2.529 0.000) 8.568 3.543 ( -7.695 -8.065 0.000) 11.147 3.567 ( 9.717 8.654 0.000) 13.012 3.693 ( -10.191 -12.139 0.000) 15.849 4.714 ( 0.068 -2.904 0.000) 2.905 4.754 ( 1.183 -4.336 0.000) 4.494 4.773 ( -0.057 12.861 0.000) 12.861 5.217 ( 6.076 -1.822 0.000) 6.343 5.908 ( 2.190 4.136 0.000) 4.680 6.316 ( 7.190 8.206 0.000) 10.910 6.872 ( -5.251 -2.107 0.000) 5.658 6.914 ( -4.755 -5.056 0.000) 6.941 ======================= Grid point 58 (25/84) ======================= q-point: ( 0.46 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.998 ( 3.893 -1.963 0.000) 4.360 3.420 ( -3.693 -8.930 0.000) 9.664 3.545 ( -4.196 -10.208 0.000) 11.037 3.704 ( 3.702 5.528 0.000) 6.654 4.702 ( -0.336 -11.670 0.000) 11.675 4.757 ( -0.818 15.834 0.000) 15.855 4.779 ( 0.569 -1.546 0.000) 1.647 5.318 ( 3.072 -2.327 0.000) 3.854 5.927 ( 0.274 3.488 0.000) 3.499 6.458 ( 5.358 10.595 0.000) 11.872 6.774 ( -3.881 -0.620 0.000) 3.930 6.843 ( -2.071 -5.129 0.000) 5.531 ======================= Grid point 70 (26/84) ======================= q-point: ( 0.38 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.935 ( 3.537 3.537 0.000) 5.002 3.397 ( -6.934 -6.934 0.000) 9.806 3.457 ( -10.281 -10.281 0.000) 14.539 3.717 ( 6.264 6.264 0.000) 8.859 4.579 ( -5.932 -5.932 0.000) 8.390 4.752 ( -0.094 -0.094 0.000) 0.133 4.986 ( 7.148 7.148 0.000) 10.109 5.205 ( 1.556 1.556 0.000) 2.201 5.961 ( 1.622 1.622 0.000) 2.293 6.489 ( 8.546 8.546 0.000) 12.086 6.813 ( -4.863 -4.863 0.000) 6.878 6.830 ( -2.374 -2.374 0.000) 3.358 ======================= Grid point 71 (27/84) ======================= q-point: ( 0.46 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.991 ( 1.660 0.954 0.000) 1.915 3.269 ( -4.104 -6.193 0.000) 7.429 3.333 ( -2.832 -10.283 0.000) 10.666 3.805 ( 2.355 4.417 0.000) 5.005 4.500 ( -1.955 -7.989 0.000) 8.225 4.745 ( -0.339 -1.862 0.000) 1.892 5.043 ( 0.360 12.431 0.000) 12.436 5.272 ( 3.065 -2.114 0.000) 3.724 5.981 ( 0.490 1.998 0.000) 2.057 6.634 ( 5.397 7.054 0.000) 8.882 6.739 ( -2.182 -5.058 0.000) 5.508 6.774 ( -3.000 0.334 0.000) 3.019 ======================= Grid point 84 (28/84) ======================= q-point: ( 0.46 0.46 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.018 ( 1.189 1.189 0.000) 1.681 3.177 ( -3.529 -3.529 0.000) 4.991 3.180 ( -3.514 -3.514 0.000) 4.970 3.868 ( 1.731 1.731 0.000) 2.448 4.395 ( -2.657 -2.657 0.000) 3.757 4.715 ( -0.928 -0.928 0.000) 1.312 5.216 ( 3.513 3.513 0.000) 4.968 5.247 ( 0.538 0.538 0.000) 0.760 6.008 ( 0.665 0.665 0.000) 0.940 6.660 ( -2.323 -2.323 0.000) 3.286 6.736 ( 3.295 3.295 0.000) 4.659 6.774 ( -0.609 -0.609 0.000) 0.861 ======================= Grid point 183 (29/84) ======================= q-point: ( 0.08 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.034 ( 16.785 16.785 16.785) 29.073 1.034 ( 16.785 16.785 16.785) 29.073 1.819 ( 30.184 30.184 30.184) 52.281 4.388 ( -1.637 -1.637 -1.637) 2.835 4.388 ( -1.637 -1.637 -1.637) 2.835 4.643 ( 2.413 2.413 2.413) 4.180 4.729 ( 1.645 1.645 1.645) 2.850 4.781 ( 3.333 3.333 3.333) 5.773 4.781 ( 3.333 3.333 3.333) 5.773 7.333 ( -7.342 -7.342 -7.342) 12.717 7.471 ( -2.070 -2.070 -2.070) 3.586 7.471 ( -2.070 -2.070 -2.070) 3.586 ======================= Grid point 184 (30/84) ======================= q-point: ( 0.15 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.452 ( 25.674 9.198 9.198) 28.782 1.493 ( 25.656 9.560 9.560) 29.000 2.499 ( 36.359 23.097 23.097) 48.876 4.330 ( -3.109 -2.981 -2.981) 5.239 4.380 ( -1.088 -2.802 -2.802) 4.109 4.641 ( -2.139 5.408 5.408) 7.942 4.788 ( 4.275 0.391 0.391) 4.310 4.874 ( 6.965 2.874 2.874) 8.064 4.897 ( 6.999 2.582 2.582) 7.894 7.127 ( -13.122 -6.964 -6.964) 16.407 7.405 ( -4.219 -2.038 -2.038) 5.109 7.413 ( -3.915 -2.016 -2.016) 4.843 ======================= Grid point 185 (31/84) ======================= q-point: ( 0.23 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.989 ( 27.611 4.631 4.631) 28.377 2.016 ( 26.590 5.519 5.519) 27.712 3.210 ( 34.500 17.056 17.056) 42.095 4.254 ( -4.466 -4.167 -4.167) 7.394 4.337 ( -3.452 -2.872 -2.872) 5.330 4.572 ( -4.721 5.903 5.903) 9.591 4.895 ( 6.323 -1.280 -1.280) 6.577 5.044 ( 9.950 3.652 3.652) 11.211 5.047 ( 7.807 1.906 1.906) 8.259 6.827 ( -16.740 -6.070 -6.070) 18.813 7.308 ( -5.497 -2.030 -2.030) 6.202 7.317 ( -5.670 -1.836 -1.836) 6.236 ======================= Grid point 186 (32/84) ======================= q-point: ( 0.31 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.524 ( 25.941 2.157 2.157) 26.119 2.532 ( 25.175 2.631 2.631) 25.448 3.821 ( 24.468 10.297 10.297) 28.474 4.174 ( -1.717 -4.008 -4.008) 5.923 4.239 ( -6.505 -2.751 -2.751) 7.580 4.459 ( -6.101 5.703 5.703) 10.113 5.028 ( 6.870 -2.762 -2.762) 7.903 5.195 ( 6.829 0.755 0.755) 6.912 5.256 ( 11.109 5.331 5.331) 13.426 6.480 ( -17.885 -4.015 -4.015) 18.765 7.193 ( -6.776 -2.005 -2.005) 7.345 7.194 ( -5.750 -1.800 -1.800) 6.288 ======================= Grid point 187 (33/84) ======================= q-point: ( 0.38 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.001 ( 21.868 0.435 0.435) 21.877 3.002 ( 21.855 0.224 0.224) 21.857 3.973 ( -2.890 -2.663 -2.663) 4.747 4.078 ( -9.688 -2.431 -2.431) 10.280 4.208 ( -6.209 3.658 3.658) 8.082 4.515 ( 17.018 4.000 4.000) 17.933 5.155 ( 5.520 -3.572 -3.572) 7.483 5.308 ( 4.458 -0.769 -0.769) 4.589 5.470 ( 10.004 7.539 7.539) 14.620 6.134 ( -16.639 -0.039 -0.039) 16.639 7.059 ( -6.301 -2.819 -2.819) 7.456 7.090 ( -4.554 -1.318 -1.318) 4.921 ======================= Grid point 188 (34/84) ======================= q-point: ( 0.46 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.353 ( 11.629 -2.094 -2.094) 12.000 3.371 ( 13.629 -1.235 -1.235) 13.741 3.873 ( -9.384 -1.745 -1.745) 9.703 3.922 ( -0.929 -3.827 -3.827) 5.492 4.031 ( -8.253 4.850 4.850) 10.731 4.823 ( 10.002 -1.405 -1.405) 10.198 5.234 ( 2.198 -3.651 -3.651) 5.612 5.367 ( 1.476 -2.136 -2.136) 3.362 5.641 ( 7.351 7.439 7.439) 12.835 5.845 ( -11.562 7.088 7.088) 15.302 6.963 ( -2.805 -3.945 -3.945) 6.244 7.025 ( -1.771 -0.859 -0.859) 2.147 ======================= Grid point 197 (35/84) ======================= q-point: ( 0.15 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.691 ( 16.545 16.545 9.741) 25.345 1.825 ( 20.645 20.645 1.666) 29.244 3.038 ( 29.439 29.439 19.282) 45.881 4.240 ( -5.000 -5.000 -0.724) 7.108 4.354 ( -1.355 -1.355 -7.473) 7.715 4.720 ( 1.612 1.612 10.431) 10.677 4.815 ( 2.284 2.284 -0.085) 3.230 4.905 ( 3.161 3.161 3.772) 5.850 5.023 ( 8.072 8.072 -1.139) 11.472 6.932 ( -12.318 -12.318 -6.164) 18.479 7.349 ( -3.767 -3.767 -2.271) 5.791 7.352 ( -3.966 -3.966 -2.251) 6.044 ======================= Grid point 198 (36/84) ======================= q-point: ( 0.23 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.078 ( 21.443 7.658 6.987) 23.817 2.265 ( 23.145 15.986 -1.829) 28.188 3.611 ( 27.254 21.667 14.384) 37.671 4.136 ( -5.543 -6.759 -1.244) 8.829 4.306 ( -3.604 -1.510 -8.146) 9.035 4.690 ( -4.037 5.943 10.715) 12.901 4.879 ( 4.088 -0.474 -1.036) 4.243 5.019 ( 7.465 -1.538 4.452) 8.827 5.205 ( 10.251 9.445 -0.468) 13.947 6.658 ( -14.789 -10.550 -5.663) 19.029 7.254 ( -5.574 -3.468 -1.950) 6.848 7.263 ( -5.094 -3.457 -2.159) 6.524 ======================= Grid point 199 (37/84) ======================= q-point: ( 0.31 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.510 ( 21.612 1.081 5.060) 22.223 2.721 ( 22.362 11.943 -4.188) 25.695 3.972 ( -1.821 -1.590 0.740) 2.528 4.108 ( 8.018 1.520 4.091) 9.128 4.231 ( 0.701 2.831 -5.756) 6.452 4.587 ( -5.319 7.336 9.622) 13.217 4.971 ( 4.987 -3.256 -2.282) 6.378 5.165 ( 6.718 -2.595 2.555) 7.642 5.429 ( 12.194 10.915 3.756) 16.791 6.368 ( -13.926 -6.942 -4.815) 16.288 7.131 ( -6.860 -3.562 -1.032) 7.798 7.156 ( -5.423 -2.452 -2.006) 6.281 ======================= Grid point 200 (38/84) ======================= q-point: ( 0.38 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.916 ( 18.853 -3.844 3.427) 19.543 3.133 ( 18.745 8.086 -5.602) 21.170 3.865 ( -5.529 -7.455 -2.121) 9.521 4.014 ( -9.419 -4.635 -3.291) 11.002 4.314 ( -4.809 7.138 1.761) 8.785 4.624 ( 10.835 6.285 6.408) 14.070 5.065 ( 4.238 -5.666 -3.082) 7.718 5.273 ( 4.106 -2.478 -1.478) 5.018 5.663 ( 10.183 11.131 8.117) 17.131 6.131 ( -9.148 -0.127 -1.966) 9.358 6.987 ( -7.296 -3.817 -1.401) 8.352 7.059 ( -4.093 -1.834 -0.812) 4.558 ======================= Grid point 201 (39/84) ======================= q-point: ( 0.46 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.222 ( 10.222 -7.536 1.772) 12.823 3.430 ( 9.705 4.046 -6.167) 12.190 3.787 ( -2.108 -7.646 -1.955) 8.169 3.824 ( -7.787 -4.940 -2.538) 9.565 4.179 ( -5.592 9.916 2.029) 11.563 4.813 ( 5.759 0.357 2.120) 6.147 5.125 ( 1.563 -7.481 -2.956) 8.194 5.329 ( 1.494 -1.848 -4.721) 5.285 5.773 ( 0.868 5.511 4.059) 6.899 6.047 ( 0.141 12.320 8.011) 14.696 6.867 ( -3.832 -5.020 -3.058) 7.017 7.002 ( -1.535 -1.439 0.336) 2.131 ======================= Grid point 211 (40/84) ======================= q-point: ( 0.23 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.291 ( 13.236 13.236 4.800) 19.324 2.614 ( 18.612 18.612 -4.462) 26.697 3.997 ( -7.116 -7.116 -0.995) 10.113 4.005 ( 15.626 15.626 9.340) 23.992 4.266 ( -1.358 -1.358 -8.633) 8.844 4.779 ( 1.457 1.457 13.969) 14.120 4.873 ( 0.112 0.112 -0.521) 0.545 5.015 ( 2.369 2.369 0.892) 3.467 5.418 ( 12.023 12.023 1.852) 17.103 6.432 ( -11.581 -11.581 -5.958) 17.428 7.169 ( -5.005 -5.005 -1.200) 7.179 7.187 ( -4.217 -4.217 -1.961) 6.278 ======================= Grid point 212 (41/84) ======================= q-point: ( 0.31 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.577 ( 15.022 5.480 3.112) 16.291 2.982 ( 17.988 13.884 -6.200) 23.554 3.849 ( -7.650 -7.666 -1.036) 10.880 4.061 ( -6.884 -6.413 -3.187) 9.934 4.396 ( 10.601 10.433 0.354) 14.878 4.737 ( -3.657 6.986 9.402) 12.271 4.894 ( 1.659 -4.441 -0.105) 4.742 5.109 ( 5.927 -2.504 -0.156) 6.436 5.672 ( 12.987 13.129 6.787) 19.675 6.227 ( -8.436 -6.659 -6.768) 12.701 7.059 ( -6.095 -3.615 -0.044) 7.087 7.098 ( -4.585 -3.351 -0.875) 5.746 ======================= Grid point 213 (42/84) ======================= q-point: ( 0.38 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.865 ( 13.598 -1.050 1.856) 13.764 3.307 ( 14.252 9.050 -7.004) 18.278 3.703 ( -6.620 -8.477 -1.150) 10.816 3.883 ( -9.926 -8.643 -2.499) 13.396 4.487 ( -2.344 9.991 -1.194) 10.332 4.736 ( 4.499 4.387 5.606) 8.421 4.927 ( 1.356 -7.976 -0.767) 8.127 5.219 ( 4.812 -2.825 -3.388) 6.528 5.864 ( 3.608 7.271 6.879) 10.640 6.156 ( 3.558 4.180 -2.845) 6.183 6.926 ( -7.246 -2.186 0.566) 7.590 7.016 ( -3.467 -2.536 1.032) 4.418 ======================= Grid point 214 (43/84) ======================= q-point: ( 0.46 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.084 ( 7.039 -5.978 0.983) 9.287 3.513 ( 5.338 3.339 -6.936) 9.368 3.614 ( -1.677 -8.250 -1.169) 8.499 3.708 ( -6.262 -7.216 -1.313) 9.644 4.411 ( -3.037 13.065 -2.431) 13.632 4.826 ( 3.057 0.559 3.787) 4.899 4.940 ( 0.044 -10.695 -1.259) 10.769 5.291 ( 2.112 -2.029 -6.055) 6.726 5.854 ( -1.532 2.467 -0.300) 2.919 6.299 ( 6.474 12.536 7.507) 15.982 6.791 ( -5.083 -2.184 -0.612) 5.566 6.968 ( -1.276 -1.992 2.678) 3.573 ======================= Grid point 225 (44/84) ======================= q-point: ( 0.31 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.714 ( 7.941 7.941 1.460) 11.325 3.252 ( 13.010 13.010 -7.501) 19.870 3.694 ( -7.906 -7.906 -0.666) 11.201 3.887 ( -10.488 -10.488 -2.918) 15.116 4.585 ( 8.213 8.213 -0.426) 11.623 4.793 ( -4.154 -4.154 -0.188) 5.878 4.836 ( 1.269 1.269 6.639) 6.877 5.095 ( 1.702 1.702 -2.834) 3.718 5.901 ( 7.129 7.129 8.214) 13.004 6.153 ( 1.779 1.779 -7.169) 7.597 6.987 ( -3.805 -3.805 1.553) 5.600 7.026 ( -3.761 -3.761 1.102) 5.432 ======================= Grid point 226 (45/84) ======================= q-point: ( 0.38 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.874 ( 7.757 1.765 0.282) 7.960 3.477 ( 8.170 7.311 -7.884) 13.504 3.543 ( -6.107 -7.082 -0.667) 9.375 3.678 ( -9.695 -11.663 -1.585) 15.249 4.686 ( -0.488 6.583 -3.831) 7.632 4.730 ( 0.442 -3.687 1.222) 3.909 4.838 ( -0.367 1.431 4.035) 4.297 5.167 ( 4.719 -2.159 -4.687) 6.992 5.915 ( -2.543 -0.866 -3.056) 4.069 6.321 ( 11.535 10.840 3.647) 16.244 6.898 ( -5.317 -0.824 2.968) 6.145 6.960 ( -2.621 -3.009 3.835) 5.535 ======================= Grid point 227 (46/84) ======================= q-point: ( 0.46 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.997 ( 3.816 -2.649 -0.232) 4.651 3.414 ( -3.569 -8.302 -0.681) 9.062 3.534 ( -4.151 -10.065 -1.146) 10.947 3.622 ( 3.664 4.119 -7.702) 9.472 4.647 ( -1.400 2.435 -5.344) 6.037 4.744 ( 0.016 0.533 -0.544) 0.762 4.842 ( 0.447 0.031 5.623) 5.640 5.248 ( 2.557 -2.298 -6.413) 7.276 5.868 ( -1.502 -0.875 -7.240) 7.446 6.527 ( 8.020 10.264 7.178) 14.873 6.783 ( -5.518 1.020 3.301) 6.510 6.925 ( -0.915 -2.239 5.892) 6.369 ======================= Grid point 239 (47/84) ======================= q-point: ( 0.38 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.925 ( 3.058 3.058 -0.981) 4.434 3.394 ( -6.905 -6.905 -0.277) 9.769 3.449 ( -9.992 -9.992 -0.841) 14.156 3.619 ( 5.512 5.512 -9.197) 12.056 4.580 ( -5.764 -5.764 0.137) 8.153 4.823 ( 1.981 1.981 6.171) 6.778 4.912 ( 5.135 5.135 -6.159) 9.521 5.150 ( 1.053 1.053 -5.256) 5.463 5.874 ( -2.645 -2.645 -10.285) 10.944 6.531 ( 9.799 9.799 6.608) 15.353 6.881 ( -1.570 -1.570 4.896) 5.376 6.902 ( -2.386 -2.386 7.214) 7.964 ======================= Grid point 240 (48/84) ======================= q-point: ( 0.46 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.977 ( 1.625 0.411 -1.471) 2.230 3.268 ( -4.068 -6.249 -0.166) 7.458 3.326 ( -2.955 -10.101 -0.800) 10.554 3.695 ( 1.997 3.006 -10.084) 10.710 4.503 ( -1.967 -7.384 0.213) 7.644 4.839 ( 0.228 -0.146 8.413) 8.418 4.964 ( 0.567 11.784 -7.346) 13.898 5.203 ( 2.659 -2.055 -6.394) 7.224 5.827 ( -1.601 -2.985 -14.627) 15.014 6.697 ( 6.752 6.725 7.386) 12.057 6.815 ( -4.307 1.889 6.927) 8.373 6.882 ( -0.343 -2.021 8.733) 8.971 ======================= Grid point 253 (49/84) ======================= q-point: ( 0.46 0.46 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.997 ( 1.061 1.061 -2.101) 2.582 3.173 ( -3.541 -3.541 -0.491) 5.032 3.176 ( -3.581 -3.581 -0.454) 5.084 3.736 ( 1.039 1.039 -11.744) 11.835 4.404 ( -2.514 -2.514 0.877) 3.662 4.836 ( -0.072 -0.072 10.558) 10.558 5.138 ( 4.115 4.115 -7.304) 9.338 5.178 ( 0.365 0.365 -6.431) 6.452 5.769 ( -2.149 -2.149 -20.322) 20.548 6.793 ( 3.049 3.049 9.364) 10.309 6.844 ( -0.402 -0.402 10.730) 10.745 6.848 ( -0.302 -0.302 6.985) 6.998 ======================= Grid point 366 (50/84) ======================= q-point: ( 0.15 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.953 ( 13.704 13.704 13.704) 23.737 1.953 ( 13.704 13.704 13.704) 23.737 3.486 ( 24.056 24.056 24.056) 41.666 4.217 ( -4.157 -4.157 -4.157) 7.201 4.217 ( -4.157 -4.157 -4.157) 7.201 4.826 ( 1.255 1.255 1.255) 2.173 4.938 ( 8.934 8.934 8.934) 15.474 4.994 ( 3.013 3.013 3.013) 5.219 4.994 ( 3.013 3.013 3.013) 5.219 6.753 ( -11.594 -11.594 -11.594) 20.081 7.290 ( -3.738 -3.738 -3.738) 6.474 7.290 ( -3.738 -3.738 -3.738) 6.474 ======================= Grid point 367 (51/84) ======================= q-point: ( 0.23 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.255 ( 18.134 8.326 8.326) 21.621 2.300 ( 18.542 8.150 8.150) 21.832 3.936 ( 19.495 16.522 16.522) 30.431 4.089 ( -5.494 -3.749 -3.749) 7.635 4.166 ( -3.484 -5.619 -5.619) 8.676 4.864 ( 2.509 -0.389 -0.389) 2.568 4.891 ( -5.617 8.810 8.810) 13.666 5.095 ( 4.829 0.203 0.203) 4.837 5.246 ( 16.252 8.075 8.075) 19.863 6.493 ( -13.958 -10.662 -10.662) 20.548 7.204 ( -4.582 -2.918 -2.918) 6.166 7.208 ( -4.739 -3.235 -3.235) 6.587 ======================= Grid point 368 (52/84) ======================= q-point: ( 0.31 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.635 ( 19.472 3.268 3.268) 20.013 2.671 ( 18.603 4.045 4.045) 19.462 3.930 ( -7.416 -4.477 -4.477) 9.751 4.074 ( -5.883 -5.506 -5.506) 9.760 4.298 ( 12.950 9.001 9.001) 18.159 4.765 ( -5.997 7.779 7.779) 12.529 4.923 ( 3.328 -2.432 -2.432) 4.786 5.177 ( 3.301 -1.875 -1.875) 4.234 5.587 ( 17.438 12.325 12.325) 24.656 6.228 ( -12.091 -9.060 -9.060) 17.616 7.099 ( -6.322 -2.071 -2.071) 6.967 7.116 ( -4.092 -1.725 -1.725) 4.764 ======================= Grid point 369 (53/84) ======================= q-point: ( 0.38 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.999 ( 16.428 -0.508 -0.508) 16.443 3.024 ( 16.639 0.537 0.537) 16.656 3.789 ( -6.203 -5.433 -5.433) 9.875 3.933 ( -8.257 -4.905 -4.905) 10.784 4.379 ( -5.088 4.602 4.602) 8.260 4.754 ( 5.561 6.031 6.031) 10.182 4.989 ( 3.172 -4.438 -4.438) 7.032 5.222 ( 1.301 -3.648 -3.648) 5.320 5.864 ( 5.357 9.667 9.667) 14.683 6.079 ( 1.894 -1.511 -1.511) 2.855 6.958 ( -7.766 -1.020 -1.020) 7.899 7.047 ( -2.714 -0.357 -0.357) 2.761 ======================= Grid point 370 (54/84) ======================= q-point: ( 0.46 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.246 ( 7.250 -3.290 -3.290) 8.615 3.303 ( 9.973 -2.183 -2.183) 10.440 3.704 ( -2.059 -5.977 -5.977) 8.700 3.767 ( -7.054 -3.738 -3.738) 8.815 4.244 ( -5.297 4.688 4.688) 8.486 4.866 ( 3.695 3.057 3.057) 5.687 5.038 ( 1.447 -5.818 -5.818) 8.354 5.234 ( 0.107 -4.534 -4.534) 6.413 5.825 ( -2.820 0.252 0.252) 2.843 6.252 ( 7.907 12.027 12.027) 18.757 6.813 ( -5.523 -1.472 -1.472) 5.902 7.011 ( -0.931 0.577 0.577) 1.238 ======================= Grid point 380 (55/84) ======================= q-point: ( 0.23 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.423 ( 10.538 10.538 8.297) 17.057 2.555 ( 14.878 14.878 -0.493) 21.046 3.966 ( -7.091 -7.091 -2.173) 10.261 4.065 ( -4.555 -4.555 -6.949) 9.475 4.258 ( 13.704 13.704 10.429) 22.008 4.861 ( 0.088 0.088 -0.595) 0.607 5.017 ( -0.937 -0.937 7.548) 7.664 5.033 ( -0.090 -0.090 0.807) 0.817 5.521 ( 17.587 17.587 10.179) 26.873 6.258 ( -12.369 -12.369 -11.251) 20.799 7.138 ( -3.609 -3.609 -1.686) 5.375 7.144 ( -3.142 -3.142 -2.106) 4.917 ======================= Grid point 381 (56/84) ======================= q-point: ( 0.31 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.662 ( 12.956 2.909 5.293) 14.294 2.859 ( 15.342 11.035 -5.135) 19.584 3.814 ( -7.975 -6.746 -2.573) 10.758 3.956 ( -6.506 -6.717 -6.594) 11.442 4.480 ( 7.572 8.881 6.834) 13.524 4.852 ( -2.519 -1.416 1.394) 3.209 4.905 ( -3.236 2.881 1.542) 4.599 5.098 ( 2.542 -5.546 -1.356) 6.250 5.888 ( 18.417 17.036 14.255) 28.855 6.030 ( -9.970 -10.198 -12.198) 18.767 7.061 ( -4.411 -1.633 0.357) 4.717 7.085 ( -2.623 -1.441 -0.234) 3.002 ======================= Grid point 382 (57/84) ======================= q-point: ( 0.38 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.915 ( 11.962 -3.332 2.897) 12.751 3.145 ( 12.950 6.993 -8.429) 16.961 3.664 ( -6.754 -6.945 -2.937) 10.123 3.810 ( -7.973 -7.719 -4.772) 12.080 4.490 ( -5.893 6.424 1.342) 8.820 4.839 ( 3.089 0.568 4.224) 5.264 4.902 ( 1.002 -2.371 -1.263) 2.868 5.135 ( 1.316 -4.865 -4.889) 7.022 5.860 ( -6.673 -6.972 -9.576) 13.596 6.237 ( 15.887 14.194 12.701) 24.802 6.958 ( -6.042 0.967 2.774) 6.718 7.044 ( -1.492 0.016 1.817) 2.351 ======================= Grid point 383 (58/84) ======================= q-point: ( 0.46 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.104 ( 5.963 -7.405 0.808) 9.541 3.343 ( 6.139 2.787 -9.446) 11.606 3.570 ( -2.294 -7.230 -3.191) 8.229 3.663 ( -5.615 -6.844 -3.400) 9.483 4.367 ( -4.303 7.529 -1.531) 8.806 4.891 ( 1.079 -7.581 -1.889) 7.887 4.931 ( 1.477 2.095 4.654) 5.313 5.154 ( 0.628 -3.431 -7.102) 7.912 5.767 ( -2.520 -5.797 -8.990) 10.990 6.499 ( 9.408 12.427 12.098) 19.731 6.830 ( -5.937 3.017 4.718) 8.162 7.025 ( -0.459 0.870 3.067) 3.221 ======================= Grid point 394 (59/84) ======================= q-point: ( 0.31 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.753 ( 5.863 5.863 2.352) 8.619 3.082 ( 11.068 11.068 -8.711) 17.912 3.671 ( -7.574 -7.574 -1.721) 10.849 3.802 ( -8.705 -8.705 -5.499) 13.484 4.632 ( 5.353 5.353 4.407) 8.760 4.789 ( -4.041 -4.041 -0.232) 5.720 4.935 ( -1.440 -1.440 3.159) 3.759 5.016 ( -0.288 -0.288 -4.966) 4.983 5.845 ( -8.224 -8.224 -15.142) 19.093 6.226 ( 16.373 16.373 15.425) 27.822 7.032 ( -1.585 -1.585 2.927) 3.687 7.059 ( -1.056 -1.056 2.232) 2.686 ======================= Grid point 395 (60/84) ======================= q-point: ( 0.38 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.879 ( 6.440 -0.333 0.135) 6.450 3.280 ( 8.400 6.415 -11.018) 15.268 3.522 ( -6.993 -7.184 -1.560) 10.147 3.627 ( -8.308 -10.454 -3.553) 13.818 4.618 ( -5.485 5.329 -2.978) 8.207 4.744 ( -0.683 -5.792 0.074) 5.833 4.937 ( 2.408 2.722 5.905) 6.934 5.049 ( 2.478 -3.312 -6.835) 7.989 5.707 ( -5.767 -7.789 -16.363) 19.018 6.519 ( 12.913 13.367 14.462) 23.549 6.987 ( -3.043 1.656 5.706) 6.676 7.047 ( -0.271 0.221 4.678) 4.691 ======================= Grid point 396 (61/84) ======================= q-point: ( 0.46 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.984 ( 3.396 -4.436 -1.134) 5.701 3.369 ( -0.479 -2.350 -6.294) 6.735 3.447 ( 0.821 -3.052 -6.614) 7.330 3.497 ( -4.099 -9.514 -2.937) 10.767 4.524 ( -2.848 6.802 -6.371) 9.745 4.729 ( -0.547 -6.292 -0.984) 6.392 4.980 ( 1.143 2.125 7.031) 7.434 5.099 ( 1.947 -2.088 -7.134) 7.684 5.622 ( -2.458 -8.262 -16.949) 19.015 6.728 ( 7.720 10.368 12.237) 17.800 6.911 ( -4.294 4.530 8.805) 10.793 7.045 ( -0.016 1.057 6.016) 6.108 ======================= Grid point 408 (62/84) ======================= q-point: ( 0.38 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.896 ( 1.765 1.765 -1.999) 3.198 3.382 ( -6.807 -6.807 -1.087) 9.688 3.390 ( 4.067 4.067 -12.617) 13.867 3.421 ( -9.012 -9.012 -2.462) 12.980 4.585 ( -5.298 -5.298 0.368) 7.501 4.751 ( 3.130 3.130 -8.226) 9.342 4.993 ( 2.901 2.901 8.787) 9.698 5.015 ( 0.093 0.093 -7.864) 7.865 5.562 ( -6.676 -6.676 -19.649) 21.799 6.750 ( 9.736 9.736 13.886) 19.555 7.009 ( 0.215 0.215 7.606) 7.612 7.054 ( 0.562 0.562 7.743) 7.783 ======================= Grid point 409 (63/84) ======================= q-point: ( 0.46 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.933 ( 1.407 -0.872 -2.954) 3.386 3.259 ( -3.877 -6.241 -0.894) 7.401 3.299 ( -3.231 -9.405 -1.974) 10.139 3.451 ( 1.599 0.881 -13.569) 13.691 4.508 ( -2.070 -5.331 0.224) 5.723 4.783 ( 0.464 8.743 -10.285) 13.507 5.010 ( 0.208 0.922 0.698) 1.175 5.073 ( 3.005 -0.565 1.249) 3.303 5.454 ( -3.536 -8.467 -21.658) 23.522 6.901 ( 5.393 6.890 12.041) 14.884 6.990 ( -2.072 3.116 10.076) 10.748 7.065 ( 0.382 0.966 9.263) 9.321 ======================= Grid point 422 (64/84) ======================= q-point: ( 0.46 0.46 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.937 ( 0.746 0.746 -3.842) 3.984 3.156 ( -3.502 -3.502 -1.211) 5.098 3.161 ( -3.698 -3.698 -1.143) 5.354 3.463 ( 0.291 0.291 -14.790) 14.795 4.428 ( -2.156 -2.156 1.579) 3.434 4.944 ( 5.646 5.646 -11.661) 14.132 5.018 ( 0.062 0.062 -9.029) 9.029 5.078 ( 1.079 1.079 12.083) 12.179 5.300 ( -6.056 -6.056 -25.120) 26.540 6.999 ( 3.029 3.029 11.282) 12.068 7.025 ( 0.343 0.343 10.261) 10.273 7.081 ( 0.478 0.478 11.874) 11.893 ======================= Grid point 549 (65/84) ======================= q-point: ( 0.23 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.609 ( 7.925 7.925 7.925) 13.726 2.609 ( 7.925 7.925 7.925) 13.726 3.909 ( -5.974 -5.974 -5.974) 10.347 3.909 ( -5.974 -5.974 -5.974) 10.347 4.473 ( 10.318 10.318 10.318) 17.871 4.850 ( -0.523 -0.523 -0.523) 0.906 5.043 ( -1.901 -1.901 -1.901) 3.292 5.043 ( -1.901 -1.901 -1.901) 3.292 5.800 ( 13.160 13.160 13.160) 22.794 6.033 ( -7.590 -7.590 -7.590) 13.146 7.108 ( -1.443 -1.443 -1.443) 2.500 7.108 ( -1.443 -1.443 -1.443) 2.500 ======================= Grid point 550 (66/84) ======================= q-point: ( 0.31 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.764 ( 9.305 2.986 2.986) 10.219 2.816 ( 10.647 2.795 2.795) 11.356 3.742 ( -8.411 -4.823 -4.823) 10.829 3.818 ( -5.478 -7.123 -7.123) 11.467 4.629 ( 3.129 7.532 7.532) 11.102 4.840 ( -0.562 -1.802 -1.802) 2.610 4.926 ( -4.966 0.444 0.444) 5.006 5.038 ( -0.446 -4.986 -4.986) 7.065 5.744 ( -8.717 -14.307 -14.307) 22.031 6.229 ( 16.910 17.774 17.774) 30.295 7.076 ( -1.813 1.098 1.098) 2.387 7.090 ( -0.476 0.709 0.709) 1.110 ======================= Grid point 551 (67/84) ======================= q-point: ( 0.38 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.950 ( 8.776 -1.992 -1.992) 9.217 3.023 ( 9.964 -0.901 -0.901) 10.045 3.586 ( -6.925 -4.848 -4.848) 9.745 3.694 ( -6.797 -6.907 -6.907) 11.900 4.527 ( -9.974 2.110 2.110) 10.411 4.828 ( -0.599 -3.507 -3.507) 4.996 4.954 ( 5.420 4.609 4.609) 8.477 5.029 ( -0.402 -5.758 -5.758) 8.153 5.602 ( -5.596 -14.355 -14.355) 21.059 6.536 ( 13.505 15.689 15.689) 25.974 7.030 ( -2.951 4.180 4.180) 6.607 7.087 ( 0.059 2.477 2.477) 3.503 ======================= Grid point 552 (68/84) ======================= q-point: ( 0.46 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.080 ( 3.661 -5.080 -5.080) 8.063 3.191 ( 5.895 -3.619 -3.619) 7.806 3.488 ( -2.588 -5.038 -5.038) 7.581 3.568 ( -4.851 -6.225 -6.225) 10.051 4.362 ( -5.022 1.143 1.143) 5.275 4.818 ( -0.281 -4.749 -4.749) 6.722 5.024 ( -0.128 -5.754 -5.754) 8.138 5.047 ( 2.927 5.946 5.946) 8.903 5.522 ( -2.240 -14.454 -14.454) 20.563 6.750 ( 7.419 12.542 12.542) 19.227 6.959 ( -3.736 7.527 7.527) 11.281 7.089 ( 0.096 3.415 3.415) 4.830 ======================= Grid point 563 (69/84) ======================= q-point: ( 0.31 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.804 ( 2.850 2.850 2.799) 4.908 2.936 ( 7.051 7.051 -5.617) 11.445 3.621 ( -7.222 -7.222 -3.440) 10.777 3.678 ( -7.055 -7.055 -6.888) 12.125 4.730 ( -0.636 -0.636 4.823) 4.907 4.785 ( -3.686 -3.686 -0.185) 5.216 4.909 ( -2.515 -2.515 -5.549) 6.591 4.959 ( 1.051 1.051 -1.419) 2.055 5.501 ( -9.596 -9.596 -17.130) 21.854 6.564 ( 15.409 15.409 17.233) 27.783 7.099 ( 1.010 1.010 3.714) 3.979 7.112 ( 1.431 1.431 2.962) 3.588 ======================= Grid point 564 (70/84) ======================= q-point: ( 0.38 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.878 ( 4.183 -2.941 -0.303) 5.122 3.072 ( 6.349 3.412 -9.365) 11.818 3.475 ( -6.966 -6.276 -3.373) 9.965 3.537 ( -6.950 -8.793 -5.465) 12.470 4.565 ( -10.504 1.546 -2.537) 10.916 4.737 ( -2.031 -6.071 -0.737) 6.444 4.913 ( 0.941 -5.071 -6.439) 8.250 5.035 ( 5.572 2.998 1.051) 6.414 5.356 ( -5.676 -9.117 -15.119) 18.545 6.824 ( 10.587 12.850 15.244) 22.574 7.111 ( 0.102 3.651 6.436) 7.400 7.142 ( 1.382 3.000 4.805) 5.831 ======================= Grid point 565 (71/84) ======================= q-point: ( 0.46 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.948 ( 2.272 -6.175 -2.648) 7.092 3.172 ( 3.232 -0.009 -10.411) 10.902 3.372 ( -2.784 -6.561 -3.438) 7.914 3.420 ( -4.091 -8.501 -4.959) 10.658 4.408 ( -4.414 3.262 -5.342) 7.659 4.707 ( -0.805 -6.346 -1.246) 6.517 4.926 ( 0.300 -4.007 -7.030) 8.098 5.146 ( 5.120 1.404 1.563) 5.534 5.257 ( -4.555 -9.012 -13.856) 17.145 6.977 ( 4.601 10.058 12.182) 16.454 7.101 ( -0.921 6.235 9.770) 11.626 7.162 ( 0.550 3.930 5.751) 6.987 ======================= Grid point 577 (72/84) ======================= q-point: ( 0.38 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.850 ( 0.012 0.012 -2.455) 2.455 3.135 ( 2.759 2.759 -12.248) 12.855 3.345 ( -6.577 -6.577 -2.758) 9.702 3.363 ( -8.041 -8.041 -3.721) 11.966 4.579 ( -2.610 -2.610 -9.785) 10.459 4.595 ( -4.652 -4.652 0.657) 6.611 4.858 ( -0.651 -0.651 -7.556) 7.612 5.048 ( -1.642 -1.642 -13.396) 13.596 5.265 ( -0.688 -0.688 1.181) 1.530 7.033 ( 8.038 8.038 13.832) 17.904 7.167 ( 2.064 2.064 7.993) 8.509 7.202 ( 2.798 2.798 7.078) 8.109 ======================= Grid point 578 (73/84) ======================= q-point: ( 0.46 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.864 ( 0.838 -2.215 -3.826) 4.500 3.172 ( 0.753 -0.059 -12.252) 12.275 3.236 ( -2.882 -6.869 -2.977) 8.022 3.249 ( -3.436 -7.949 -3.410) 9.307 4.475 ( -3.126 2.575 -7.954) 8.926 4.588 ( -0.059 -4.493 -6.066) 7.549 4.865 ( 0.730 -2.107 -8.023) 8.327 5.007 ( -1.618 -8.229 -19.547) 21.271 5.274 ( 0.641 3.495 9.889) 10.508 7.140 ( 2.791 6.231 11.499) 13.374 7.192 ( 0.486 3.154 9.827) 10.332 7.243 ( 1.200 3.860 8.698) 9.591 ======================= Grid point 591 (74/84) ======================= q-point: ( 0.46 0.46 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.851 ( 0.282 0.282 -4.730) 4.747 3.124 ( -3.272 -3.272 -2.079) 5.073 3.129 ( -3.603 -3.603 -2.093) 5.509 3.174 ( -0.096 -0.096 -14.017) 14.017 4.464 ( -1.736 -1.736 1.940) 3.129 4.582 ( 1.734 1.734 -26.594) 26.707 4.842 ( -0.183 -0.183 -8.365) 8.369 4.884 ( -3.915 -3.915 -13.304) 14.410 5.324 ( 1.332 1.332 11.514) 11.667 7.222 ( 2.137 2.137 10.851) 11.264 7.233 ( 0.993 0.993 10.323) 10.418 7.298 ( 1.443 1.443 9.902) 10.110 ======================= Grid point 732 (75/84) ======================= q-point: ( 0.31 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.864 ( 0.809 0.809 0.809) 1.401 2.864 ( 0.809 0.809 0.809) 1.401 3.533 ( -6.485 -6.485 -6.485) 11.232 3.533 ( -6.485 -6.485 -6.485) 11.232 4.772 ( -1.926 -1.926 -1.926) 3.335 4.796 ( -5.478 -5.478 -5.478) 9.489 4.816 ( -4.474 -4.474 -4.474) 7.750 4.816 ( -4.474 -4.474 -4.474) 7.750 5.287 ( -3.015 -3.015 -3.015) 5.222 6.877 ( 13.815 13.815 13.815) 23.928 7.173 ( 3.332 3.332 3.332) 5.771 7.173 ( 3.332 3.332 3.332) 5.771 ======================= Grid point 733 (76/84) ======================= q-point: ( 0.38 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.857 ( 0.641 -2.302 -2.302) 3.317 2.937 ( 3.971 -3.369 -3.369) 6.202 3.387 ( -6.766 -5.345 -5.345) 10.145 3.414 ( -6.308 -6.860 -6.860) 11.572 4.497 ( -14.224 -5.422 -5.422) 16.159 4.713 ( -2.818 -2.386 -2.386) 4.396 4.803 ( 0.094 -5.234 -5.234) 7.402 4.839 ( 0.880 -12.975 -12.975) 18.371 5.301 ( 2.295 4.795 4.795) 7.160 7.096 ( 8.071 11.835 11.835) 18.581 7.232 ( 2.853 5.426 5.426) 8.186 7.235 ( 2.553 4.680 4.680) 7.093 ======================= Grid point 734 (77/84) ======================= q-point: ( 0.46 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.873 ( 0.594 -4.809 -4.809) 6.827 3.011 ( 2.778 -5.376 -5.376) 8.094 3.286 ( -2.898 -5.242 -5.242) 7.960 3.305 ( -3.787 -6.451 -6.451) 9.878 4.297 ( -5.349 -6.577 -6.577) 10.730 4.675 ( -0.975 -2.635 -2.635) 3.852 4.805 ( 0.104 -4.981 -4.981) 7.045 4.848 ( 0.187 -14.796 -14.796) 20.926 5.339 ( 1.143 7.258 7.258) 10.328 7.193 ( 2.061 9.391 9.391) 13.441 7.270 ( 0.919 5.397 5.397) 7.688 7.281 ( 1.690 7.793 7.793) 11.149 ======================= Grid point 746 (78/84) ======================= q-point: ( 0.38 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.806 ( -1.867 -1.867 -1.762) 3.174 2.922 ( 0.723 0.723 -8.937) 8.995 3.272 ( -6.046 -6.046 -4.578) 9.699 3.276 ( -6.566 -6.566 -4.945) 10.520 4.310 ( -11.697 -11.697 -17.978) 24.430 4.610 ( -3.983 -3.983 0.777) 5.685 4.727 ( -1.108 -1.108 -5.470) 5.689 4.748 ( -1.786 -1.786 -10.745) 11.038 5.401 ( 4.431 4.431 7.698) 9.926 7.277 ( 5.937 5.937 10.452) 13.407 7.315 ( 3.548 3.548 6.721) 8.387 7.333 ( 4.705 4.705 6.112) 9.034 ======================= Grid point 747 (79/84) ======================= q-point: ( 0.46 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.788 ( -0.212 -3.019 -3.613) 4.713 2.944 ( 1.002 -1.991 -9.940) 10.187 3.172 ( -2.955 -6.016 -3.823) 7.716 3.175 ( -3.223 -6.108 -4.077) 8.019 4.134 ( -4.864 -9.030 -22.230) 24.482 4.568 ( -0.871 -5.640 0.725) 5.752 4.711 ( -0.480 -1.540 -8.919) 9.064 4.732 ( -0.219 -2.626 -6.381) 6.904 5.462 ( 1.623 4.948 8.495) 9.964 7.343 ( 1.471 5.644 8.869) 10.615 7.363 ( 1.214 4.045 7.357) 8.483 7.408 ( 2.196 4.790 7.660) 9.297 ======================= Grid point 760 (80/84) ======================= q-point: ( 0.46 0.46 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.758 ( -0.423 -0.423 -4.456) 4.495 2.921 ( -0.458 -0.458 -11.258) 11.276 3.076 ( -2.803 -2.803 -2.622) 4.753 3.079 ( -3.044 -3.044 -2.797) 5.133 3.995 ( -3.941 -3.941 -30.221) 30.731 4.502 ( -1.383 -1.383 1.851) 2.693 4.697 ( -0.379 -0.379 -6.108) 6.131 4.700 ( -0.614 -0.614 -6.758) 6.814 5.529 ( 1.736 1.736 8.996) 9.325 7.417 ( 1.934 1.934 8.596) 9.021 7.422 ( 1.521 1.521 8.454) 8.724 7.470 ( 1.507 1.507 7.436) 7.736 ======================= Grid point 915 (81/84) ======================= q-point: ( 0.38 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.786 ( -2.531 -2.531 -2.531) 4.385 2.786 ( -2.531 -2.531 -2.531) 4.385 3.173 ( -5.074 -5.074 -5.074) 8.789 3.173 ( -5.074 -5.074 -5.074) 8.789 3.941 ( -16.798 -16.798 -16.798) 29.095 4.619 ( -2.023 -2.023 -2.023) 3.504 4.619 ( -2.023 -2.023 -2.023) 3.504 4.648 ( -2.019 -2.019 -2.019) 3.497 5.531 ( 5.230 5.230 5.230) 9.058 7.427 ( 4.501 4.501 4.501) 7.796 7.427 ( 4.501 4.501 4.501) 7.796 7.453 ( 5.877 5.877 5.877) 10.179 ======================= Grid point 916 (82/84) ======================= q-point: ( 0.46 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.725 ( -1.426 -2.744 -2.744) 4.134 2.792 ( 0.470 -5.176 -5.176) 7.335 3.090 ( -2.587 -4.157 -4.157) 6.423 3.093 ( -2.478 -4.030 -4.030) 6.214 3.706 ( -6.168 -18.861 -18.861) 27.377 4.576 ( -0.972 -1.662 -1.662) 2.543 4.595 ( -0.890 -1.513 -1.513) 2.318 4.642 ( 0.173 -3.201 -3.201) 4.531 5.602 ( 1.824 5.479 5.479) 7.960 7.486 ( 1.526 5.051 5.051) 7.304 7.487 ( 1.566 5.516 5.516) 7.956 7.532 ( 2.004 4.621 4.621) 6.836 ======================= Grid point 929 (83/84) ======================= q-point: ( 0.46 0.46 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.681 ( -1.360 -1.360 -3.094) 3.643 2.736 ( -0.982 -0.982 -7.213) 7.346 3.025 ( -2.105 -2.105 -2.525) 3.903 3.025 ( -2.131 -2.131 -2.334) 3.812 3.445 ( -6.729 -6.729 -23.864) 25.691 4.534 ( -1.140 -1.140 1.317) 2.082 4.598 ( -0.155 -0.155 -3.632) 3.639 4.602 ( -0.541 -0.541 -3.505) 3.587 5.675 ( 1.885 1.885 5.670) 6.266 7.558 ( 2.049 2.049 5.543) 6.255 7.561 ( 1.891 1.891 5.468) 6.087 7.589 ( 1.254 1.254 4.524) 4.859 ======================= Grid point 1098 (84/84) ======================= q-point: ( 0.46 0.46 0.46) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.98e-04 5.98e-04 5.98e-04 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.639 ( -1.781 -1.781 -1.781) 3.085 2.639 ( -1.781 -1.781 -1.781) 3.085 2.988 ( -1.224 -1.224 -1.224) 2.120 2.988 ( -1.224 -1.224 -1.224) 2.120 3.108 ( -8.883 -8.883 -8.883) 15.385 4.552 ( -0.467 -0.467 -0.467) 0.809 4.552 ( -0.467 -0.467 -0.467) 0.809 4.558 ( -0.850 -0.850 -0.850) 1.472 5.751 ( 1.933 1.933 1.933) 3.349 7.635 ( 2.017 2.017 2.017) 3.493 7.635 ( 2.017 2.017 2.017) 3.493 7.647 ( 1.279 1.279 1.279) 2.216 =================== 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 1165.429 1165.429 1165.429 0.000 -0.000 0.000 3/26364 20.0 361.779 361.779 361.779 0.000 -0.000 0.000 3/26364 30.0 112.110 112.110 112.110 0.000 -0.000 0.000 3/26364 40.0 55.014 55.014 55.014 0.000 -0.000 0.000 3/26364 50.0 35.932 35.932 35.932 0.000 -0.000 0.000 3/26364 60.0 27.224 27.224 27.224 0.000 -0.000 0.000 3/26364 70.0 22.345 22.345 22.345 0.000 -0.000 0.000 3/26364 80.0 19.219 19.219 19.219 0.000 -0.000 0.000 3/26364 90.0 17.021 17.021 17.021 0.000 -0.000 0.000 3/26364 100.0 15.371 15.371 15.371 0.000 -0.000 0.000 3/26364 110.0 14.072 14.072 14.072 0.000 -0.000 0.000 3/26364 120.0 13.011 13.011 13.011 0.000 -0.000 0.000 3/26364 130.0 12.122 12.122 12.122 0.000 -0.000 0.000 3/26364 140.0 11.362 11.362 11.362 0.000 -0.000 0.000 3/26364 150.0 10.702 10.702 10.702 0.000 -0.000 0.000 3/26364 160.0 10.121 10.121 10.121 0.000 -0.000 0.000 3/26364 170.0 9.605 9.605 9.605 0.000 -0.000 0.000 3/26364 180.0 9.143 9.143 9.143 0.000 -0.000 0.000 3/26364 190.0 8.725 8.725 8.725 0.000 -0.000 0.000 3/26364 200.0 8.346 8.346 8.346 0.000 -0.000 0.000 3/26364 210.0 8.000 8.000 8.000 0.000 -0.000 0.000 3/26364 220.0 7.683 7.683 7.683 0.000 -0.000 0.000 3/26364 230.0 7.391 7.391 7.391 0.000 -0.000 0.000 3/26364 240.0 7.121 7.121 7.121 0.000 -0.000 0.000 3/26364 250.0 6.871 6.871 6.871 0.000 -0.000 0.000 3/26364 260.0 6.638 6.638 6.638 0.000 -0.000 0.000 3/26364 270.0 6.421 6.421 6.421 0.000 -0.000 0.000 3/26364 280.0 6.218 6.218 6.218 0.000 -0.000 0.000 3/26364 290.0 6.028 6.028 6.028 0.000 -0.000 0.000 3/26364 300.0 5.849 5.849 5.849 0.000 -0.000 0.000 3/26364 310.0 5.681 5.681 5.681 0.000 -0.000 0.000 3/26364 320.0 5.522 5.522 5.522 0.000 -0.000 0.000 3/26364 330.0 5.372 5.372 5.372 0.000 -0.000 0.000 3/26364 340.0 5.231 5.231 5.231 0.000 -0.000 0.000 3/26364 350.0 5.096 5.096 5.096 0.000 -0.000 0.000 3/26364 360.0 4.969 4.969 4.969 0.000 -0.000 0.000 3/26364 370.0 4.848 4.848 4.848 0.000 -0.000 0.000 3/26364 380.0 4.732 4.732 4.732 0.000 -0.000 0.000 3/26364 390.0 4.622 4.622 4.622 0.000 -0.000 0.000 3/26364 400.0 4.518 4.518 4.518 0.000 -0.000 0.000 3/26364 410.0 4.418 4.418 4.418 0.000 -0.000 0.000 3/26364 420.0 4.322 4.322 4.322 0.000 -0.000 0.000 3/26364 430.0 4.230 4.230 4.230 0.000 -0.000 0.000 3/26364 440.0 4.143 4.143 4.143 0.000 -0.000 0.000 3/26364 450.0 4.059 4.059 4.059 0.000 -0.000 0.000 3/26364 460.0 3.978 3.978 3.978 0.000 -0.000 0.000 3/26364 470.0 3.900 3.900 3.900 0.000 -0.000 0.000 3/26364 480.0 3.826 3.826 3.826 0.000 -0.000 0.000 3/26364 490.0 3.754 3.754 3.754 0.000 -0.000 0.000 3/26364 500.0 3.685 3.685 3.685 0.000 -0.000 0.000 3/26364 510.0 3.618 3.618 3.618 0.000 -0.000 0.000 3/26364 520.0 3.554 3.554 3.554 0.000 -0.000 0.000 3/26364 530.0 3.492 3.492 3.492 0.000 -0.000 0.000 3/26364 540.0 3.433 3.433 3.433 0.000 -0.000 0.000 3/26364 550.0 3.375 3.375 3.375 0.000 -0.000 0.000 3/26364 560.0 3.319 3.319 3.319 0.000 -0.000 0.000 3/26364 570.0 3.265 3.265 3.265 0.000 -0.000 0.000 3/26364 580.0 3.213 3.213 3.213 0.000 -0.000 0.000 3/26364 590.0 3.162 3.162 3.162 0.000 -0.000 0.000 3/26364 600.0 3.114 3.114 3.114 0.000 -0.000 0.000 3/26364 610.0 3.066 3.066 3.066 0.000 -0.000 0.000 3/26364 620.0 3.020 3.020 3.020 0.000 -0.000 0.000 3/26364 630.0 2.975 2.975 2.975 0.000 -0.000 0.000 3/26364 640.0 2.932 2.932 2.932 0.000 -0.000 0.000 3/26364 650.0 2.890 2.890 2.890 0.000 -0.000 0.000 3/26364 660.0 2.849 2.849 2.849 0.000 -0.000 0.000 3/26364 670.0 2.809 2.809 2.809 0.000 -0.000 0.000 3/26364 680.0 2.771 2.771 2.771 0.000 -0.000 0.000 3/26364 690.0 2.733 2.733 2.733 0.000 -0.000 0.000 3/26364 700.0 2.697 2.697 2.697 0.000 -0.000 0.000 3/26364 710.0 2.661 2.661 2.661 0.000 -0.000 0.000 3/26364 720.0 2.626 2.626 2.626 0.000 -0.000 0.000 3/26364 730.0 2.593 2.593 2.593 0.000 -0.000 0.000 3/26364 740.0 2.560 2.560 2.560 0.000 -0.000 0.000 3/26364 750.0 2.528 2.528 2.528 0.000 -0.000 0.000 3/26364 760.0 2.496 2.496 2.496 0.000 -0.000 0.000 3/26364 770.0 2.466 2.466 2.466 0.000 -0.000 0.000 3/26364 780.0 2.436 2.436 2.436 0.000 -0.000 0.000 3/26364 790.0 2.407 2.407 2.407 0.000 -0.000 0.000 3/26364 800.0 2.379 2.379 2.379 0.000 -0.000 0.000 3/26364 810.0 2.351 2.351 2.351 0.000 -0.000 0.000 3/26364 820.0 2.324 2.324 2.324 0.000 -0.000 0.000 3/26364 830.0 2.297 2.297 2.297 0.000 -0.000 0.000 3/26364 840.0 2.272 2.272 2.272 0.000 -0.000 0.000 3/26364 850.0 2.246 2.246 2.246 0.000 -0.000 0.000 3/26364 860.0 2.222 2.222 2.222 0.000 -0.000 0.000 3/26364 870.0 2.198 2.198 2.198 0.000 -0.000 0.000 3/26364 880.0 2.174 2.174 2.174 0.000 -0.000 0.000 3/26364 890.0 2.151 2.151 2.151 0.000 -0.000 0.000 3/26364 900.0 2.128 2.128 2.128 0.000 -0.000 0.000 3/26364 910.0 2.106 2.106 2.106 0.000 -0.000 0.000 3/26364 920.0 2.084 2.084 2.084 0.000 -0.000 0.000 3/26364 930.0 2.063 2.063 2.063 0.000 -0.000 0.000 3/26364 940.0 2.042 2.042 2.042 0.000 -0.000 0.000 3/26364 950.0 2.022 2.022 2.022 0.000 -0.000 0.000 3/26364 960.0 2.002 2.002 2.002 0.000 -0.000 0.000 3/26364 970.0 1.982 1.982 1.982 0.000 -0.000 0.000 3/26364 980.0 1.963 1.963 1.963 0.000 -0.000 0.000 3/26364 990.0 1.944 1.944 1.944 0.000 -0.000 0.000 3/26364 1000.0 1.926 1.926 1.926 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:41:36]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|