----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 00:14:48]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) Number of symmetry operations in supercell: 384 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.725673860443566 0.000000000000000 0.000000000000000 b -1.862836930221783 3.226528209359767 0.000000000000000 c 0.000000000000000 0.000000000000000 14.704334960000001 Atomic positions (fractional): *1 S 0.33333333333333 0.66666666666667 0.09646753023511 32.065 2 S 0.66666666666667 0.33333333333333 0.59646753023511 32.065 3 S 0.33333333333333 0.66666666666667 0.40353246976489 32.065 4 S 0.66666666666667 0.33333333333333 0.90353246976489 32.065 *5 Sc 0.00000000000000 0.00000000000000 0.50000000000000 44.956 6 Sc 0.00000000000000 0.00000000000000 0.00000000000000 44.956 *7 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 8 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.725673860443566 0.000000000000000 0.000000000000000 b -1.862836930221783 3.226528209359767 0.000000000000000 c 0.000000000000000 0.000000000000000 14.704334960000001 Atomic positions (fractional): *1 S 0.33333333333333 0.66666666666667 0.09646753023511 32.065 > 1 2 S 0.66666666666667 0.33333333333333 0.59646753023511 32.065 > 2 3 S 0.33333333333333 0.66666666666667 0.40353246976489 32.065 > 3 4 S 0.66666666666667 0.33333333333333 0.90353246976489 32.065 > 4 *5 Sc 0.00000000000000 0.00000000000000 0.50000000000000 44.956 > 5 6 Sc 0.00000000000000 0.00000000000000 0.00000000000000 44.956 > 6 *7 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 > 7 8 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 14.902695441774265 0.000000000000000 0.000000000000000 b -7.451347720887132 12.906112837439068 0.000000000000000 c 0.000000000000000 0.000000000000000 14.704334960000001 Atomic positions (fractional): *1 S 0.08333333333333 0.16666666666667 0.09646753023511 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.09646753023511 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.09646753023511 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.09646753023511 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.09646753023511 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.09646753023511 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.09646753023511 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.09646753023511 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.09646753023511 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.09646753023511 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.09646753023511 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.09646753023511 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.09646753023511 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.09646753023511 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.09646753023511 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.09646753023511 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.59646753023511 32.065 > 2 18 S 0.41666666666667 0.08333333333333 0.59646753023511 32.065 > 2 19 S 0.66666666666667 0.08333333333333 0.59646753023511 32.065 > 2 20 S 0.91666666666667 0.08333333333333 0.59646753023511 32.065 > 2 21 S 0.16666666666667 0.33333333333333 0.59646753023511 32.065 > 2 22 S 0.41666666666667 0.33333333333333 0.59646753023511 32.065 > 2 23 S 0.66666666666667 0.33333333333333 0.59646753023511 32.065 > 2 24 S 0.91666666666667 0.33333333333333 0.59646753023511 32.065 > 2 25 S 0.16666666666667 0.58333333333333 0.59646753023511 32.065 > 2 26 S 0.41666666666667 0.58333333333333 0.59646753023511 32.065 > 2 27 S 0.66666666666667 0.58333333333333 0.59646753023511 32.065 > 2 28 S 0.91666666666667 0.58333333333333 0.59646753023511 32.065 > 2 29 S 0.16666666666667 0.83333333333333 0.59646753023511 32.065 > 2 30 S 0.41666666666667 0.83333333333333 0.59646753023511 32.065 > 2 31 S 0.66666666666667 0.83333333333333 0.59646753023511 32.065 > 2 32 S 0.91666666666667 0.83333333333333 0.59646753023511 32.065 > 2 33 S 0.08333333333333 0.16666666666667 0.40353246976489 32.065 > 3 34 S 0.33333333333333 0.16666666666667 0.40353246976489 32.065 > 3 35 S 0.58333333333333 0.16666666666667 0.40353246976489 32.065 > 3 36 S 0.83333333333333 0.16666666666667 0.40353246976489 32.065 > 3 37 S 0.08333333333333 0.41666666666667 0.40353246976489 32.065 > 3 38 S 0.33333333333333 0.41666666666667 0.40353246976489 32.065 > 3 39 S 0.58333333333333 0.41666666666667 0.40353246976489 32.065 > 3 40 S 0.83333333333333 0.41666666666667 0.40353246976489 32.065 > 3 41 S 0.08333333333333 0.66666666666667 0.40353246976489 32.065 > 3 42 S 0.33333333333333 0.66666666666667 0.40353246976489 32.065 > 3 43 S 0.58333333333333 0.66666666666667 0.40353246976489 32.065 > 3 44 S 0.83333333333333 0.66666666666667 0.40353246976489 32.065 > 3 45 S 0.08333333333333 0.91666666666667 0.40353246976489 32.065 > 3 46 S 0.33333333333333 0.91666666666667 0.40353246976489 32.065 > 3 47 S 0.58333333333333 0.91666666666667 0.40353246976489 32.065 > 3 48 S 0.83333333333333 0.91666666666667 0.40353246976489 32.065 > 3 49 S 0.16666666666667 0.08333333333333 0.90353246976489 32.065 > 4 50 S 0.41666666666667 0.08333333333333 0.90353246976489 32.065 > 4 51 S 0.66666666666667 0.08333333333333 0.90353246976489 32.065 > 4 52 S 0.91666666666667 0.08333333333333 0.90353246976489 32.065 > 4 53 S 0.16666666666667 0.33333333333333 0.90353246976489 32.065 > 4 54 S 0.41666666666667 0.33333333333333 0.90353246976489 32.065 > 4 55 S 0.66666666666667 0.33333333333333 0.90353246976489 32.065 > 4 56 S 0.91666666666667 0.33333333333333 0.90353246976489 32.065 > 4 57 S 0.16666666666667 0.58333333333333 0.90353246976489 32.065 > 4 58 S 0.41666666666667 0.58333333333333 0.90353246976489 32.065 > 4 59 S 0.66666666666667 0.58333333333333 0.90353246976489 32.065 > 4 60 S 0.91666666666667 0.58333333333333 0.90353246976489 32.065 > 4 61 S 0.16666666666667 0.83333333333333 0.90353246976489 32.065 > 4 62 S 0.41666666666667 0.83333333333333 0.90353246976489 32.065 > 4 63 S 0.66666666666667 0.83333333333333 0.90353246976489 32.065 > 4 64 S 0.91666666666667 0.83333333333333 0.90353246976489 32.065 > 4 *65 Sc 0.00000000000000 0.00000000000000 0.50000000000000 44.956 > 5 66 Sc 0.25000000000000 0.00000000000000 0.50000000000000 44.956 > 5 67 Sc 0.50000000000000 0.00000000000000 0.50000000000000 44.956 > 5 68 Sc 0.75000000000000 0.00000000000000 0.50000000000000 44.956 > 5 69 Sc 0.00000000000000 0.25000000000000 0.50000000000000 44.956 > 5 70 Sc 0.25000000000000 0.25000000000000 0.50000000000000 44.956 > 5 71 Sc 0.50000000000000 0.25000000000000 0.50000000000000 44.956 > 5 72 Sc 0.75000000000000 0.25000000000000 0.50000000000000 44.956 > 5 73 Sc 0.00000000000000 0.50000000000000 0.50000000000000 44.956 > 5 74 Sc 0.25000000000000 0.50000000000000 0.50000000000000 44.956 > 5 75 Sc 0.50000000000000 0.50000000000000 0.50000000000000 44.956 > 5 76 Sc 0.75000000000000 0.50000000000000 0.50000000000000 44.956 > 5 77 Sc 0.00000000000000 0.75000000000000 0.50000000000000 44.956 > 5 78 Sc 0.25000000000000 0.75000000000000 0.50000000000000 44.956 > 5 79 Sc 0.50000000000000 0.75000000000000 0.50000000000000 44.956 > 5 80 Sc 0.75000000000000 0.75000000000000 0.50000000000000 44.956 > 5 81 Sc 0.00000000000000 0.00000000000000 0.00000000000000 44.956 > 6 82 Sc 0.25000000000000 0.00000000000000 0.00000000000000 44.956 > 6 83 Sc 0.50000000000000 0.00000000000000 0.00000000000000 44.956 > 6 84 Sc 0.75000000000000 0.00000000000000 0.00000000000000 44.956 > 6 85 Sc 0.00000000000000 0.25000000000000 0.00000000000000 44.956 > 6 86 Sc 0.25000000000000 0.25000000000000 0.00000000000000 44.956 > 6 87 Sc 0.50000000000000 0.25000000000000 0.00000000000000 44.956 > 6 88 Sc 0.75000000000000 0.25000000000000 0.00000000000000 44.956 > 6 89 Sc 0.00000000000000 0.50000000000000 0.00000000000000 44.956 > 6 90 Sc 0.25000000000000 0.50000000000000 0.00000000000000 44.956 > 6 91 Sc 0.50000000000000 0.50000000000000 0.00000000000000 44.956 > 6 92 Sc 0.75000000000000 0.50000000000000 0.00000000000000 44.956 > 6 93 Sc 0.00000000000000 0.75000000000000 0.00000000000000 44.956 > 6 94 Sc 0.25000000000000 0.75000000000000 0.00000000000000 44.956 > 6 95 Sc 0.50000000000000 0.75000000000000 0.00000000000000 44.956 > 6 96 Sc 0.75000000000000 0.75000000000000 0.00000000000000 44.956 > 6 *97 Tl 0.16666666666667 0.08333333333333 0.25000000000000 204.383 > 7 98 Tl 0.41666666666667 0.08333333333333 0.25000000000000 204.383 > 7 99 Tl 0.66666666666667 0.08333333333333 0.25000000000000 204.383 > 7 100 Tl 0.91666666666667 0.08333333333333 0.25000000000000 204.383 > 7 101 Tl 0.16666666666667 0.33333333333333 0.25000000000000 204.383 > 7 102 Tl 0.41666666666667 0.33333333333333 0.25000000000000 204.383 > 7 103 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 > 7 104 Tl 0.91666666666667 0.33333333333333 0.25000000000000 204.383 > 7 105 Tl 0.16666666666667 0.58333333333333 0.25000000000000 204.383 > 7 106 Tl 0.41666666666667 0.58333333333333 0.25000000000000 204.383 > 7 107 Tl 0.66666666666667 0.58333333333333 0.25000000000000 204.383 > 7 108 Tl 0.91666666666667 0.58333333333333 0.25000000000000 204.383 > 7 109 Tl 0.16666666666667 0.83333333333333 0.25000000000000 204.383 > 7 110 Tl 0.41666666666667 0.83333333333333 0.25000000000000 204.383 > 7 111 Tl 0.66666666666667 0.83333333333333 0.25000000000000 204.383 > 7 112 Tl 0.91666666666667 0.83333333333333 0.25000000000000 204.383 > 7 113 Tl 0.08333333333333 0.16666666666667 0.75000000000000 204.383 > 8 114 Tl 0.33333333333333 0.16666666666667 0.75000000000000 204.383 > 8 115 Tl 0.58333333333333 0.16666666666667 0.75000000000000 204.383 > 8 116 Tl 0.83333333333333 0.16666666666667 0.75000000000000 204.383 > 8 117 Tl 0.08333333333333 0.41666666666667 0.75000000000000 204.383 > 8 118 Tl 0.33333333333333 0.41666666666667 0.75000000000000 204.383 > 8 119 Tl 0.58333333333333 0.41666666666667 0.75000000000000 204.383 > 8 120 Tl 0.83333333333333 0.41666666666667 0.75000000000000 204.383 > 8 121 Tl 0.08333333333333 0.66666666666667 0.75000000000000 204.383 > 8 122 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 > 8 123 Tl 0.58333333333333 0.66666666666667 0.75000000000000 204.383 > 8 124 Tl 0.83333333333333 0.66666666666667 0.75000000000000 204.383 > 8 125 Tl 0.08333333333333 0.91666666666667 0.75000000000000 204.383 > 8 126 Tl 0.33333333333333 0.91666666666667 0.75000000000000 204.383 > 8 127 Tl 0.58333333333333 0.91666666666667 0.75000000000000 204.383 > 8 128 Tl 0.83333333333333 0.91666666666667 0.75000000000000 204.383 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.4378485 0.0000000 0.0000000 0.0000000 9.4378485 0.0000000 0.0000000 0.0000000 10.5854917 -------------------------- Born effective charges -------------------------- 1 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 2 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 3 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 4 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 5 Sc 4.5835394 -0.0000000 0.0000000 0.0000000 4.5835394 0.0000000 0.0000000 0.0000000 2.7214515 6 Sc 4.5835394 -0.0000000 0.0000000 0.0000000 4.5835394 0.0000000 0.0000000 0.0000000 2.7214515 7 Tl 1.6756646 0.0000000 0.0000000 0.0000000 1.6756646 0.0000000 0.0000000 0.0000000 3.7313297 8 Tl 1.6756646 0.0000000 0.0000000 0.0000000 1.6756646 0.0000000 0.0000000 0.0000000 3.7313297 ---------------------------------------------------------------------------- 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: 384/384 Permutation basis: 5136/5136 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 259 Number of blocks in projector: 186 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 123 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 76 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 60 Use standard eigh solver. Tree of FC basis block matrices: - (259, 253), data: False |-- (60, 57), data: True |-- (76, 76), data: True |-- (123, 120), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.008 Solver_block: 80 / 80 - Time: 0.196 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.206 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 384/384 Permutation basis: 5136/5136 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 259 Number of blocks in projector: 186 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 123 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 76 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 60 Use standard eigh solver. Tree of FC basis block matrices: - (259, 253), data: False |-- (60, 57), data: True |-- (76, 76), data: True |-- (123, 120), data: True Max drift after symmetrization by symfc projector: 0.00000000 (yy) 0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 00:14:53]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 00:14:53]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.725673860443566 0.000000000000000 0.000000000000000 b -1.862836930221783 3.226528209359767 0.000000000000000 c 0.000000000000000 0.000000000000000 14.704334960000001 Atomic positions (fractional): 1 S 0.33333333333333 0.66666666666667 0.09646753023511 32.065 2 S 0.66666666666667 0.33333333333333 0.59646753023511 32.065 3 S 0.33333333333333 0.66666666666667 0.40353246976489 32.065 4 S 0.66666666666667 0.33333333333333 0.90353246976489 32.065 5 Sc 0.00000000000000 0.00000000000000 0.50000000000000 44.956 6 Sc 0.00000000000000 0.00000000000000 0.00000000000000 44.956 7 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 8 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.902695441774265 0.000000000000000 0.000000000000000 b -7.451347720887132 12.906112837439068 0.000000000000000 c 0.000000000000000 0.000000000000000 14.704334960000001 Atomic positions (fractional): 1 S 0.08333333333333 0.16666666666667 0.09646753023511 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.09646753023511 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.09646753023511 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.09646753023511 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.09646753023511 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.09646753023511 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.09646753023511 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.09646753023511 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.09646753023511 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.09646753023511 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.09646753023511 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.09646753023511 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.09646753023511 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.09646753023511 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.09646753023511 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.09646753023511 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.59646753023511 32.065 > 17 18 S 0.41666666666667 0.08333333333333 0.59646753023511 32.065 > 17 19 S 0.66666666666667 0.08333333333333 0.59646753023511 32.065 > 17 20 S 0.91666666666667 0.08333333333333 0.59646753023511 32.065 > 17 21 S 0.16666666666667 0.33333333333333 0.59646753023511 32.065 > 17 22 S 0.41666666666667 0.33333333333333 0.59646753023511 32.065 > 17 23 S 0.66666666666667 0.33333333333333 0.59646753023511 32.065 > 17 24 S 0.91666666666667 0.33333333333333 0.59646753023511 32.065 > 17 25 S 0.16666666666667 0.58333333333333 0.59646753023511 32.065 > 17 26 S 0.41666666666667 0.58333333333333 0.59646753023511 32.065 > 17 27 S 0.66666666666667 0.58333333333333 0.59646753023511 32.065 > 17 28 S 0.91666666666667 0.58333333333333 0.59646753023511 32.065 > 17 29 S 0.16666666666667 0.83333333333333 0.59646753023511 32.065 > 17 30 S 0.41666666666667 0.83333333333333 0.59646753023511 32.065 > 17 31 S 0.66666666666667 0.83333333333333 0.59646753023511 32.065 > 17 32 S 0.91666666666667 0.83333333333333 0.59646753023511 32.065 > 17 33 S 0.08333333333333 0.16666666666667 0.40353246976489 32.065 > 33 34 S 0.33333333333333 0.16666666666667 0.40353246976489 32.065 > 33 35 S 0.58333333333333 0.16666666666667 0.40353246976489 32.065 > 33 36 S 0.83333333333333 0.16666666666667 0.40353246976489 32.065 > 33 37 S 0.08333333333333 0.41666666666667 0.40353246976489 32.065 > 33 38 S 0.33333333333333 0.41666666666667 0.40353246976489 32.065 > 33 39 S 0.58333333333333 0.41666666666667 0.40353246976489 32.065 > 33 40 S 0.83333333333333 0.41666666666667 0.40353246976489 32.065 > 33 41 S 0.08333333333333 0.66666666666667 0.40353246976489 32.065 > 33 42 S 0.33333333333333 0.66666666666667 0.40353246976489 32.065 > 33 43 S 0.58333333333333 0.66666666666667 0.40353246976489 32.065 > 33 44 S 0.83333333333333 0.66666666666667 0.40353246976489 32.065 > 33 45 S 0.08333333333333 0.91666666666667 0.40353246976489 32.065 > 33 46 S 0.33333333333333 0.91666666666667 0.40353246976489 32.065 > 33 47 S 0.58333333333333 0.91666666666667 0.40353246976489 32.065 > 33 48 S 0.83333333333333 0.91666666666667 0.40353246976489 32.065 > 33 49 S 0.16666666666667 0.08333333333333 0.90353246976489 32.065 > 49 50 S 0.41666666666667 0.08333333333333 0.90353246976489 32.065 > 49 51 S 0.66666666666667 0.08333333333333 0.90353246976489 32.065 > 49 52 S 0.91666666666667 0.08333333333333 0.90353246976489 32.065 > 49 53 S 0.16666666666667 0.33333333333333 0.90353246976489 32.065 > 49 54 S 0.41666666666667 0.33333333333333 0.90353246976489 32.065 > 49 55 S 0.66666666666667 0.33333333333333 0.90353246976489 32.065 > 49 56 S 0.91666666666667 0.33333333333333 0.90353246976489 32.065 > 49 57 S 0.16666666666667 0.58333333333333 0.90353246976489 32.065 > 49 58 S 0.41666666666667 0.58333333333333 0.90353246976489 32.065 > 49 59 S 0.66666666666667 0.58333333333333 0.90353246976489 32.065 > 49 60 S 0.91666666666667 0.58333333333333 0.90353246976489 32.065 > 49 61 S 0.16666666666667 0.83333333333333 0.90353246976489 32.065 > 49 62 S 0.41666666666667 0.83333333333333 0.90353246976489 32.065 > 49 63 S 0.66666666666667 0.83333333333333 0.90353246976489 32.065 > 49 64 S 0.91666666666667 0.83333333333333 0.90353246976489 32.065 > 49 65 Sc 0.00000000000000 0.00000000000000 0.50000000000000 44.956 > 65 66 Sc 0.25000000000000 0.00000000000000 0.50000000000000 44.956 > 65 67 Sc 0.50000000000000 0.00000000000000 0.50000000000000 44.956 > 65 68 Sc 0.75000000000000 0.00000000000000 0.50000000000000 44.956 > 65 69 Sc 0.00000000000000 0.25000000000000 0.50000000000000 44.956 > 65 70 Sc 0.25000000000000 0.25000000000000 0.50000000000000 44.956 > 65 71 Sc 0.50000000000000 0.25000000000000 0.50000000000000 44.956 > 65 72 Sc 0.75000000000000 0.25000000000000 0.50000000000000 44.956 > 65 73 Sc 0.00000000000000 0.50000000000000 0.50000000000000 44.956 > 65 74 Sc 0.25000000000000 0.50000000000000 0.50000000000000 44.956 > 65 75 Sc 0.50000000000000 0.50000000000000 0.50000000000000 44.956 > 65 76 Sc 0.75000000000000 0.50000000000000 0.50000000000000 44.956 > 65 77 Sc 0.00000000000000 0.75000000000000 0.50000000000000 44.956 > 65 78 Sc 0.25000000000000 0.75000000000000 0.50000000000000 44.956 > 65 79 Sc 0.50000000000000 0.75000000000000 0.50000000000000 44.956 > 65 80 Sc 0.75000000000000 0.75000000000000 0.50000000000000 44.956 > 65 81 Sc 0.00000000000000 0.00000000000000 0.00000000000000 44.956 > 81 82 Sc 0.25000000000000 0.00000000000000 0.00000000000000 44.956 > 81 83 Sc 0.50000000000000 0.00000000000000 0.00000000000000 44.956 > 81 84 Sc 0.75000000000000 0.00000000000000 0.00000000000000 44.956 > 81 85 Sc 0.00000000000000 0.25000000000000 0.00000000000000 44.956 > 81 86 Sc 0.25000000000000 0.25000000000000 0.00000000000000 44.956 > 81 87 Sc 0.50000000000000 0.25000000000000 0.00000000000000 44.956 > 81 88 Sc 0.75000000000000 0.25000000000000 0.00000000000000 44.956 > 81 89 Sc 0.00000000000000 0.50000000000000 0.00000000000000 44.956 > 81 90 Sc 0.25000000000000 0.50000000000000 0.00000000000000 44.956 > 81 91 Sc 0.50000000000000 0.50000000000000 0.00000000000000 44.956 > 81 92 Sc 0.75000000000000 0.50000000000000 0.00000000000000 44.956 > 81 93 Sc 0.00000000000000 0.75000000000000 0.00000000000000 44.956 > 81 94 Sc 0.25000000000000 0.75000000000000 0.00000000000000 44.956 > 81 95 Sc 0.50000000000000 0.75000000000000 0.00000000000000 44.956 > 81 96 Sc 0.75000000000000 0.75000000000000 0.00000000000000 44.956 > 81 97 Tl 0.16666666666667 0.08333333333333 0.25000000000000 204.383 > 97 98 Tl 0.41666666666667 0.08333333333333 0.25000000000000 204.383 > 97 99 Tl 0.66666666666667 0.08333333333333 0.25000000000000 204.383 > 97 100 Tl 0.91666666666667 0.08333333333333 0.25000000000000 204.383 > 97 101 Tl 0.16666666666667 0.33333333333333 0.25000000000000 204.383 > 97 102 Tl 0.41666666666667 0.33333333333333 0.25000000000000 204.383 > 97 103 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 > 97 104 Tl 0.91666666666667 0.33333333333333 0.25000000000000 204.383 > 97 105 Tl 0.16666666666667 0.58333333333333 0.25000000000000 204.383 > 97 106 Tl 0.41666666666667 0.58333333333333 0.25000000000000 204.383 > 97 107 Tl 0.66666666666667 0.58333333333333 0.25000000000000 204.383 > 97 108 Tl 0.91666666666667 0.58333333333333 0.25000000000000 204.383 > 97 109 Tl 0.16666666666667 0.83333333333333 0.25000000000000 204.383 > 97 110 Tl 0.41666666666667 0.83333333333333 0.25000000000000 204.383 > 97 111 Tl 0.66666666666667 0.83333333333333 0.25000000000000 204.383 > 97 112 Tl 0.91666666666667 0.83333333333333 0.25000000000000 204.383 > 97 113 Tl 0.08333333333333 0.16666666666667 0.75000000000000 204.383 > 113 114 Tl 0.33333333333333 0.16666666666667 0.75000000000000 204.383 > 113 115 Tl 0.58333333333333 0.16666666666667 0.75000000000000 204.383 > 113 116 Tl 0.83333333333333 0.16666666666667 0.75000000000000 204.383 > 113 117 Tl 0.08333333333333 0.41666666666667 0.75000000000000 204.383 > 113 118 Tl 0.33333333333333 0.41666666666667 0.75000000000000 204.383 > 113 119 Tl 0.58333333333333 0.41666666666667 0.75000000000000 204.383 > 113 120 Tl 0.83333333333333 0.41666666666667 0.75000000000000 204.383 > 113 121 Tl 0.08333333333333 0.66666666666667 0.75000000000000 204.383 > 113 122 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 > 113 123 Tl 0.58333333333333 0.66666666666667 0.75000000000000 204.383 > 113 124 Tl 0.83333333333333 0.66666666666667 0.75000000000000 204.383 > 113 125 Tl 0.08333333333333 0.91666666666667 0.75000000000000 204.383 > 113 126 Tl 0.33333333333333 0.91666666666667 0.75000000000000 204.383 > 113 127 Tl 0.58333333333333 0.91666666666667 0.75000000000000 204.383 > 113 128 Tl 0.83333333333333 0.91666666666667 0.75000000000000 204.383 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.4378485 0.0000000 0.0000000 0.0000000 9.4378485 0.0000000 0.0000000 0.0000000 10.5854917 -------------------------- Born effective charges -------------------------- 1 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 2 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 3 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 4 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 5 Sc 4.5835394 -0.0000000 0.0000000 0.0000000 4.5835394 0.0000000 0.0000000 0.0000000 2.7214515 6 Sc 4.5835394 -0.0000000 0.0000000 0.0000000 4.5835394 0.0000000 0.0000000 0.0000000 2.7214515 7 Tl 1.6756646 0.0000000 0.0000000 0.0000000 1.6756646 0.0000000 0.0000000 0.0000000 3.7313297 8 Tl 1.6756646 0.0000000 0.0000000 0.0000000 1.6756646 0.0000000 0.0000000 0.0000000 3.7313297 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 97, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000009 (yzy) 0.00000009 (yzy) 0.00000009 (yyz) 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-08 00:15:00]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 00:15:00]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.725673860443566 0.000000000000000 0.000000000000000 b -1.862836930221783 3.226528209359767 0.000000000000000 c 0.000000000000000 0.000000000000000 14.704334960000001 Atomic positions (fractional): 1 S 0.33333333333333 0.66666666666667 0.09646753023511 32.065 2 S 0.66666666666667 0.33333333333333 0.59646753023511 32.065 3 S 0.33333333333333 0.66666666666667 0.40353246976489 32.065 4 S 0.66666666666667 0.33333333333333 0.90353246976489 32.065 5 Sc 0.00000000000000 0.00000000000000 0.50000000000000 44.956 6 Sc 0.00000000000000 0.00000000000000 0.00000000000000 44.956 7 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 8 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.902695441774265 0.000000000000000 0.000000000000000 b -7.451347720887132 12.906112837439068 0.000000000000000 c 0.000000000000000 0.000000000000000 14.704334960000001 Atomic positions (fractional): 1 S 0.08333333333333 0.16666666666667 0.09646753023511 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.09646753023511 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.09646753023511 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.09646753023511 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.09646753023511 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.09646753023511 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.09646753023511 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.09646753023511 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.09646753023511 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.09646753023511 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.09646753023511 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.09646753023511 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.09646753023511 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.09646753023511 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.09646753023511 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.09646753023511 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.59646753023511 32.065 > 17 18 S 0.41666666666667 0.08333333333333 0.59646753023511 32.065 > 17 19 S 0.66666666666667 0.08333333333333 0.59646753023511 32.065 > 17 20 S 0.91666666666667 0.08333333333333 0.59646753023511 32.065 > 17 21 S 0.16666666666667 0.33333333333333 0.59646753023511 32.065 > 17 22 S 0.41666666666667 0.33333333333333 0.59646753023511 32.065 > 17 23 S 0.66666666666667 0.33333333333333 0.59646753023511 32.065 > 17 24 S 0.91666666666667 0.33333333333333 0.59646753023511 32.065 > 17 25 S 0.16666666666667 0.58333333333333 0.59646753023511 32.065 > 17 26 S 0.41666666666667 0.58333333333333 0.59646753023511 32.065 > 17 27 S 0.66666666666667 0.58333333333333 0.59646753023511 32.065 > 17 28 S 0.91666666666667 0.58333333333333 0.59646753023511 32.065 > 17 29 S 0.16666666666667 0.83333333333333 0.59646753023511 32.065 > 17 30 S 0.41666666666667 0.83333333333333 0.59646753023511 32.065 > 17 31 S 0.66666666666667 0.83333333333333 0.59646753023511 32.065 > 17 32 S 0.91666666666667 0.83333333333333 0.59646753023511 32.065 > 17 33 S 0.08333333333333 0.16666666666667 0.40353246976489 32.065 > 33 34 S 0.33333333333333 0.16666666666667 0.40353246976489 32.065 > 33 35 S 0.58333333333333 0.16666666666667 0.40353246976489 32.065 > 33 36 S 0.83333333333333 0.16666666666667 0.40353246976489 32.065 > 33 37 S 0.08333333333333 0.41666666666667 0.40353246976489 32.065 > 33 38 S 0.33333333333333 0.41666666666667 0.40353246976489 32.065 > 33 39 S 0.58333333333333 0.41666666666667 0.40353246976489 32.065 > 33 40 S 0.83333333333333 0.41666666666667 0.40353246976489 32.065 > 33 41 S 0.08333333333333 0.66666666666667 0.40353246976489 32.065 > 33 42 S 0.33333333333333 0.66666666666667 0.40353246976489 32.065 > 33 43 S 0.58333333333333 0.66666666666667 0.40353246976489 32.065 > 33 44 S 0.83333333333333 0.66666666666667 0.40353246976489 32.065 > 33 45 S 0.08333333333333 0.91666666666667 0.40353246976489 32.065 > 33 46 S 0.33333333333333 0.91666666666667 0.40353246976489 32.065 > 33 47 S 0.58333333333333 0.91666666666667 0.40353246976489 32.065 > 33 48 S 0.83333333333333 0.91666666666667 0.40353246976489 32.065 > 33 49 S 0.16666666666667 0.08333333333333 0.90353246976489 32.065 > 49 50 S 0.41666666666667 0.08333333333333 0.90353246976489 32.065 > 49 51 S 0.66666666666667 0.08333333333333 0.90353246976489 32.065 > 49 52 S 0.91666666666667 0.08333333333333 0.90353246976489 32.065 > 49 53 S 0.16666666666667 0.33333333333333 0.90353246976489 32.065 > 49 54 S 0.41666666666667 0.33333333333333 0.90353246976489 32.065 > 49 55 S 0.66666666666667 0.33333333333333 0.90353246976489 32.065 > 49 56 S 0.91666666666667 0.33333333333333 0.90353246976489 32.065 > 49 57 S 0.16666666666667 0.58333333333333 0.90353246976489 32.065 > 49 58 S 0.41666666666667 0.58333333333333 0.90353246976489 32.065 > 49 59 S 0.66666666666667 0.58333333333333 0.90353246976489 32.065 > 49 60 S 0.91666666666667 0.58333333333333 0.90353246976489 32.065 > 49 61 S 0.16666666666667 0.83333333333333 0.90353246976489 32.065 > 49 62 S 0.41666666666667 0.83333333333333 0.90353246976489 32.065 > 49 63 S 0.66666666666667 0.83333333333333 0.90353246976489 32.065 > 49 64 S 0.91666666666667 0.83333333333333 0.90353246976489 32.065 > 49 65 Sc 0.00000000000000 0.00000000000000 0.50000000000000 44.956 > 65 66 Sc 0.25000000000000 0.00000000000000 0.50000000000000 44.956 > 65 67 Sc 0.50000000000000 0.00000000000000 0.50000000000000 44.956 > 65 68 Sc 0.75000000000000 0.00000000000000 0.50000000000000 44.956 > 65 69 Sc 0.00000000000000 0.25000000000000 0.50000000000000 44.956 > 65 70 Sc 0.25000000000000 0.25000000000000 0.50000000000000 44.956 > 65 71 Sc 0.50000000000000 0.25000000000000 0.50000000000000 44.956 > 65 72 Sc 0.75000000000000 0.25000000000000 0.50000000000000 44.956 > 65 73 Sc 0.00000000000000 0.50000000000000 0.50000000000000 44.956 > 65 74 Sc 0.25000000000000 0.50000000000000 0.50000000000000 44.956 > 65 75 Sc 0.50000000000000 0.50000000000000 0.50000000000000 44.956 > 65 76 Sc 0.75000000000000 0.50000000000000 0.50000000000000 44.956 > 65 77 Sc 0.00000000000000 0.75000000000000 0.50000000000000 44.956 > 65 78 Sc 0.25000000000000 0.75000000000000 0.50000000000000 44.956 > 65 79 Sc 0.50000000000000 0.75000000000000 0.50000000000000 44.956 > 65 80 Sc 0.75000000000000 0.75000000000000 0.50000000000000 44.956 > 65 81 Sc 0.00000000000000 0.00000000000000 0.00000000000000 44.956 > 81 82 Sc 0.25000000000000 0.00000000000000 0.00000000000000 44.956 > 81 83 Sc 0.50000000000000 0.00000000000000 0.00000000000000 44.956 > 81 84 Sc 0.75000000000000 0.00000000000000 0.00000000000000 44.956 > 81 85 Sc 0.00000000000000 0.25000000000000 0.00000000000000 44.956 > 81 86 Sc 0.25000000000000 0.25000000000000 0.00000000000000 44.956 > 81 87 Sc 0.50000000000000 0.25000000000000 0.00000000000000 44.956 > 81 88 Sc 0.75000000000000 0.25000000000000 0.00000000000000 44.956 > 81 89 Sc 0.00000000000000 0.50000000000000 0.00000000000000 44.956 > 81 90 Sc 0.25000000000000 0.50000000000000 0.00000000000000 44.956 > 81 91 Sc 0.50000000000000 0.50000000000000 0.00000000000000 44.956 > 81 92 Sc 0.75000000000000 0.50000000000000 0.00000000000000 44.956 > 81 93 Sc 0.00000000000000 0.75000000000000 0.00000000000000 44.956 > 81 94 Sc 0.25000000000000 0.75000000000000 0.00000000000000 44.956 > 81 95 Sc 0.50000000000000 0.75000000000000 0.00000000000000 44.956 > 81 96 Sc 0.75000000000000 0.75000000000000 0.00000000000000 44.956 > 81 97 Tl 0.16666666666667 0.08333333333333 0.25000000000000 204.383 > 97 98 Tl 0.41666666666667 0.08333333333333 0.25000000000000 204.383 > 97 99 Tl 0.66666666666667 0.08333333333333 0.25000000000000 204.383 > 97 100 Tl 0.91666666666667 0.08333333333333 0.25000000000000 204.383 > 97 101 Tl 0.16666666666667 0.33333333333333 0.25000000000000 204.383 > 97 102 Tl 0.41666666666667 0.33333333333333 0.25000000000000 204.383 > 97 103 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 > 97 104 Tl 0.91666666666667 0.33333333333333 0.25000000000000 204.383 > 97 105 Tl 0.16666666666667 0.58333333333333 0.25000000000000 204.383 > 97 106 Tl 0.41666666666667 0.58333333333333 0.25000000000000 204.383 > 97 107 Tl 0.66666666666667 0.58333333333333 0.25000000000000 204.383 > 97 108 Tl 0.91666666666667 0.58333333333333 0.25000000000000 204.383 > 97 109 Tl 0.16666666666667 0.83333333333333 0.25000000000000 204.383 > 97 110 Tl 0.41666666666667 0.83333333333333 0.25000000000000 204.383 > 97 111 Tl 0.66666666666667 0.83333333333333 0.25000000000000 204.383 > 97 112 Tl 0.91666666666667 0.83333333333333 0.25000000000000 204.383 > 97 113 Tl 0.08333333333333 0.16666666666667 0.75000000000000 204.383 > 113 114 Tl 0.33333333333333 0.16666666666667 0.75000000000000 204.383 > 113 115 Tl 0.58333333333333 0.16666666666667 0.75000000000000 204.383 > 113 116 Tl 0.83333333333333 0.16666666666667 0.75000000000000 204.383 > 113 117 Tl 0.08333333333333 0.41666666666667 0.75000000000000 204.383 > 113 118 Tl 0.33333333333333 0.41666666666667 0.75000000000000 204.383 > 113 119 Tl 0.58333333333333 0.41666666666667 0.75000000000000 204.383 > 113 120 Tl 0.83333333333333 0.41666666666667 0.75000000000000 204.383 > 113 121 Tl 0.08333333333333 0.66666666666667 0.75000000000000 204.383 > 113 122 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 > 113 123 Tl 0.58333333333333 0.66666666666667 0.75000000000000 204.383 > 113 124 Tl 0.83333333333333 0.66666666666667 0.75000000000000 204.383 > 113 125 Tl 0.08333333333333 0.91666666666667 0.75000000000000 204.383 > 113 126 Tl 0.33333333333333 0.91666666666667 0.75000000000000 204.383 > 113 127 Tl 0.58333333333333 0.91666666666667 0.75000000000000 204.383 > 113 128 Tl 0.83333333333333 0.91666666666667 0.75000000000000 204.383 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 9.4378485 0.0000000 0.0000000 0.0000000 9.4378485 0.0000000 0.0000000 0.0000000 10.5854917 -------------------------- Born effective charges -------------------------- 1 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 2 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 3 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 4 S -3.1296020 0.0000000 0.0000000 0.0000000 -3.1296020 0.0000000 0.0000000 0.0000000 -3.2263906 5 Sc 4.5835394 -0.0000000 0.0000000 0.0000000 4.5835394 0.0000000 0.0000000 0.0000000 2.7214515 6 Sc 4.5835394 -0.0000000 0.0000000 0.0000000 4.5835394 0.0000000 0.0000000 0.0000000 2.7214515 7 Tl 1.6756646 0.0000000 0.0000000 0.0000000 1.6756646 0.0000000 0.0000000 0.0000000 3.7313297 8 Tl 1.6756646 0.0000000 0.0000000 0.0000000 1.6756646 0.0000000 0.0000000 0.0000000 3.7313297 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000009 (yzy) 0.00000009 (yzy) 0.00000009 (yyz) 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 ... [ 15 15 3 ] fc3-r2q-transformation over three atoms: True --------------------------- Phonon calculations ---------------------------- Use NAC by Gonze et al. (no real space sum in current implementation) PRB 50, 13035(R) (1994), PRB 55, 10355 (1997) G-cutoff distance: 0.74, Number of G-points: 291, Lambda: 0.24 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/54) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 54 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 0.000 ( 0.000 0.000 0.000) 0.000 1.070 ( 0.000 0.000 0.000) 0.000 1.070 ( 0.000 0.000 0.000) 0.000 1.769 ( 0.000 0.000 0.000) 0.000 1.769 ( 0.000 0.000 0.000) 0.000 1.857 ( 0.000 0.000 0.000) 0.000 1.857 ( 0.000 0.000 0.000) 0.000 2.338 ( 0.000 0.000 0.000) 0.000 2.694 ( 0.000 0.000 0.000) 0.000 3.525 ( 0.000 0.000 0.000) 0.000 6.324 ( 0.000 0.000 0.000) 0.000 6.324 ( 0.000 0.000 0.000) 0.000 6.398 ( 0.000 0.000 0.000) 0.000 6.398 ( 0.000 0.000 0.000) 0.000 6.651 ( 0.000 0.000 0.000) 0.000 6.651 ( 0.000 0.000 0.000) 0.000 6.822 ( 0.000 0.000 0.000) 0.000 6.822 ( 0.000 0.000 0.000) 0.000 8.582 ( 0.000 0.000 0.000) 0.000 8.638 ( 0.000 0.000 0.000) 0.000 9.401 ( 0.000 0.000 0.000) 0.000 10.825 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/54) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.472 ( 18.220 10.519 0.000) 21.039 0.488 ( 18.106 10.454 0.000) 20.907 0.794 ( 26.562 15.336 0.000) 30.671 1.080 ( 0.832 0.480 0.000) 0.960 1.118 ( 3.471 2.004 0.000) 4.007 1.902 ( 1.150 0.664 0.000) 1.328 1.933 ( 12.769 7.372 0.000) 14.744 1.952 ( 8.018 4.629 0.000) 9.258 2.233 ( 21.004 12.127 0.000) 24.254 2.423 ( 15.570 8.990 0.000) 17.979 2.760 ( 5.571 3.216 0.000) 6.433 3.386 ( -9.704 -5.603 0.000) 11.205 6.318 ( -0.479 -0.276 0.000) 0.553 6.396 ( -0.137 -0.079 0.000) 0.158 6.444 ( 8.936 5.159 0.000) 10.318 6.449 ( 4.186 2.417 0.000) 4.833 6.660 ( 0.701 0.405 0.000) 0.809 6.830 ( 0.597 0.345 0.000) 0.690 7.067 ( 17.166 9.911 0.000) 19.822 8.532 ( -4.134 -2.387 0.000) 4.773 8.673 ( 3.143 1.815 0.000) 3.630 9.292 ( -8.897 -5.137 0.000) 10.273 10.051 ( -2.975 -1.718 0.000) 3.435 10.590 ( -18.502 -10.682 0.000) 21.364 ======================= Grid point 2 (3/54) ======================= q-point: ( 0.13 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.833 ( 10.715 6.186 0.000) 12.373 0.841 ( 12.362 7.137 0.000) 14.275 1.109 ( 1.570 0.906 0.000) 1.813 1.193 ( 2.234 1.290 0.000) 2.580 1.240 ( 12.273 7.086 0.000) 14.172 1.739 ( -10.320 -5.958 0.000) 11.916 2.231 ( 14.811 8.551 0.000) 17.102 2.317 ( 18.185 10.499 0.000) 20.998 2.911 ( 28.917 16.695 0.000) 33.390 2.954 ( 10.433 6.024 0.000) 12.047 3.081 ( 33.537 19.363 0.000) 38.725 3.208 ( 1.023 0.590 0.000) 1.181 6.304 ( -0.624 -0.361 0.000) 0.721 6.396 ( 0.241 0.139 0.000) 0.278 6.585 ( 6.797 3.924 0.000) 7.848 6.680 ( 0.866 0.500 0.000) 1.000 6.689 ( 10.601 6.121 0.000) 12.241 6.844 ( 0.433 0.250 0.000) 0.500 7.457 ( 13.059 7.540 0.000) 15.079 8.391 ( -7.449 -4.300 0.000) 8.601 8.780 ( 2.921 1.686 0.000) 3.372 9.017 ( -10.784 -6.226 0.000) 12.453 9.798 ( -24.953 -14.406 0.000) 28.813 10.144 ( -10.759 -6.212 0.000) 12.423 ======================= Grid point 3 (4/54) ======================= q-point: ( 0.20 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.020 ( 5.550 3.205 0.000) 6.409 1.055 ( 5.576 3.219 0.000) 6.439 1.153 ( 2.053 1.185 0.000) 2.370 1.210 ( 0.019 0.011 0.000) 0.022 1.465 ( 7.267 4.196 0.000) 8.392 1.501 ( -8.959 -5.172 0.000) 10.345 2.623 ( 17.346 10.015 0.000) 20.029 2.759 ( 18.346 10.592 0.000) 21.184 3.127 ( 5.453 3.148 0.000) 6.296 3.232 ( 12.231 7.062 0.000) 14.124 3.859 ( 35.783 20.659 0.000) 41.319 3.886 ( 32.696 18.877 0.000) 37.754 6.288 ( -0.765 -0.441 0.000) 0.883 6.406 ( 0.450 0.260 0.000) 0.520 6.696 ( 0.349 0.201 0.000) 0.403 6.754 ( 6.949 4.012 0.000) 8.024 6.847 ( -0.284 -0.164 0.000) 0.328 6.929 ( 9.500 5.485 0.000) 10.969 7.665 ( 5.012 2.894 0.000) 5.787 8.187 ( -9.448 -5.455 0.000) 10.910 8.514 ( -17.806 -10.280 0.000) 20.560 9.002 ( -37.306 -21.539 0.000) 43.077 9.071 ( 7.843 4.528 0.000) 9.057 9.994 ( -5.213 -3.010 0.000) 6.020 ======================= Grid point 4 (5/54) ======================= q-point: ( 0.27 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.110 ( -0.797 -0.460 0.000) 0.920 1.121 ( 3.313 1.913 0.000) 3.826 1.203 ( 2.014 1.163 0.000) 2.326 1.242 ( 3.428 1.979 0.000) 3.958 1.331 ( -5.074 -2.929 0.000) 5.858 1.603 ( 4.410 2.546 0.000) 5.092 3.038 ( 17.189 9.924 0.000) 19.848 3.179 ( 16.801 9.700 0.000) 19.400 3.337 ( 11.423 6.595 0.000) 13.190 3.505 ( 10.004 5.776 0.000) 11.551 4.612 ( 27.934 16.128 0.000) 32.256 4.674 ( 31.890 18.412 0.000) 36.824 6.266 ( -1.173 -0.677 0.000) 1.354 6.412 ( -0.005 -0.003 0.000) 0.006 6.690 ( -0.945 -0.545 0.000) 1.091 6.828 ( -1.352 -0.781 0.000) 1.561 6.900 ( 4.967 2.867 0.000) 5.735 7.141 ( 8.214 4.742 0.000) 9.484 7.697 ( -3.744 -2.161 0.000) 4.323 7.913 ( -18.793 -10.850 0.000) 21.700 8.151 ( -27.116 -15.655 0.000) 31.311 8.153 ( -9.149 -5.282 0.000) 10.564 9.220 ( 3.984 2.300 0.000) 4.601 9.847 ( -7.478 -4.317 0.000) 8.635 ======================= Grid point 5 (6/54) ======================= q-point: ( 0.33 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.053 ( -2.881 -1.664 0.000) 3.327 1.188 ( 2.355 1.360 0.000) 2.720 1.237 ( -3.640 -2.101 0.000) 4.203 1.246 ( 1.602 0.925 0.000) 1.850 1.377 ( 7.365 4.252 0.000) 8.505 1.679 ( 2.086 1.204 0.000) 2.409 3.436 ( 15.966 9.218 0.000) 18.436 3.557 ( 14.806 8.548 0.000) 17.096 3.632 ( 12.478 7.204 0.000) 14.409 3.692 ( 5.523 3.189 0.000) 6.377 5.216 ( 22.603 13.050 0.000) 26.100 5.353 ( 24.274 14.014 0.000) 28.029 6.231 ( -1.785 -1.030 0.000) 2.061 6.401 ( -1.000 -0.577 0.000) 1.154 6.648 ( -2.536 -1.464 0.000) 2.928 6.782 ( -2.352 -1.358 0.000) 2.715 6.967 ( -0.080 -0.046 0.000) 0.092 7.178 ( -32.375 -18.692 0.000) 37.384 7.314 ( 5.740 3.314 0.000) 6.628 7.514 ( -7.665 -4.426 0.000) 8.851 7.905 ( -1.656 -0.956 0.000) 1.912 8.121 ( 3.694 2.133 0.000) 4.266 9.253 ( -1.049 -0.605 0.000) 1.211 9.622 ( -11.592 -6.692 0.000) 13.385 ======================= Grid point 6 (7/54) ======================= q-point: ( 0.40 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.997 ( -1.741 -1.005 0.000) 2.010 1.144 ( -3.863 -2.230 0.000) 4.460 1.233 ( 1.415 0.817 0.000) 1.633 1.279 ( 1.172 0.677 0.000) 1.354 1.556 ( 6.634 3.830 0.000) 7.660 1.712 ( 1.086 0.627 0.000) 1.254 3.774 ( 1.745 1.008 0.000) 2.015 3.783 ( 12.618 7.285 0.000) 14.570 3.874 ( 11.310 6.530 0.000) 13.060 3.899 ( 9.308 5.374 0.000) 10.748 5.683 ( 16.279 9.398 0.000) 18.797 5.799 ( 12.745 7.359 0.000) 14.717 6.181 ( -2.304 -1.330 0.000) 2.660 6.363 ( -2.136 -1.233 0.000) 2.466 6.503 ( -23.001 -13.280 0.000) 26.559 6.579 ( -2.901 -1.675 0.000) 3.350 6.726 ( -2.074 -1.198 0.000) 2.395 6.888 ( -5.395 -3.115 0.000) 6.230 7.427 ( 3.914 2.260 0.000) 4.520 7.432 ( 0.232 0.134 0.000) 0.268 8.007 ( 10.007 5.778 0.000) 11.555 8.237 ( 5.099 2.944 0.000) 5.887 9.185 ( -4.154 -2.398 0.000) 4.796 9.288 ( -16.329 -9.428 0.000) 18.855 ======================= Grid point 7 (8/54) ======================= q-point: ( 0.47 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.967 ( -0.793 -0.458 0.000) 0.916 1.074 ( -1.596 -0.922 0.000) 1.843 1.255 ( 0.431 0.249 0.000) 0.498 1.300 ( 0.482 0.278 0.000) 0.557 1.649 ( 1.249 0.721 0.000) 1.442 1.751 ( 1.734 1.001 0.000) 2.002 3.793 ( 0.207 0.119 0.000) 0.239 4.004 ( 5.114 2.952 0.000) 5.905 4.053 ( 3.318 1.915 0.000) 3.831 4.070 ( 4.499 2.598 0.000) 5.195 5.965 ( 6.571 3.794 0.000) 7.588 5.985 ( 4.259 2.459 0.000) 4.917 6.114 ( -9.486 -5.477 0.000) 10.953 6.133 ( -1.329 -0.767 0.000) 1.535 6.315 ( -1.377 -0.795 0.000) 1.589 6.527 ( -1.224 -0.707 0.000) 1.414 6.692 ( -0.702 -0.405 0.000) 0.810 6.775 ( -3.033 -1.751 0.000) 3.503 7.461 ( 1.149 0.664 0.000) 1.327 7.500 ( 1.823 1.053 0.000) 2.105 8.333 ( 2.432 1.404 0.000) 2.808 8.346 ( 16.718 9.652 0.000) 19.304 8.862 ( -18.173 -10.492 0.000) 20.985 9.093 ( -2.626 -1.516 0.000) 3.032 ======================= Grid point 16 (9/54) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.755 ( 8.122 14.068 0.000) 16.245 0.773 ( 7.832 13.566 0.000) 15.665 1.095 ( 2.201 3.812 0.000) 4.402 1.159 ( 7.423 12.856 0.000) 14.845 1.167 ( 1.285 2.226 0.000) 2.570 1.800 ( -5.600 -9.699 0.000) 11.199 2.155 ( 8.414 14.574 0.000) 16.828 2.220 ( 10.887 18.857 0.000) 21.774 2.712 ( 18.782 32.531 0.000) 37.563 2.860 ( 17.855 30.926 0.000) 35.711 2.891 ( 5.361 9.286 0.000) 10.722 3.217 ( -4.697 -8.135 0.000) 9.393 6.321 ( 0.072 0.124 0.000) 0.143 6.403 ( 0.329 0.569 0.000) 0.657 6.525 ( 2.781 4.817 0.000) 5.562 6.585 ( 3.900 6.755 0.000) 7.801 6.700 ( 2.509 4.345 0.000) 5.017 6.860 ( 1.332 2.308 0.000) 2.665 7.362 ( 8.846 15.321 0.000) 17.691 8.436 ( -3.869 -6.701 0.000) 7.738 8.755 ( 3.600 6.235 0.000) 7.199 9.088 ( -8.118 -14.061 0.000) 16.236 9.995 ( -2.557 -4.429 0.000) 5.114 10.167 ( -16.232 -28.115 0.000) 32.465 ======================= Grid point 17 (10/54) ======================= q-point: ( 0.13 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.995 ( 4.441 8.089 0.000) 9.228 0.997 ( 6.112 6.952 0.000) 9.257 1.143 ( 0.763 2.229 0.000) 2.356 1.196 ( 1.437 -0.515 0.000) 1.526 1.392 ( 6.471 6.881 0.000) 9.445 1.575 ( -7.403 -9.094 0.000) 11.726 2.519 ( 12.117 18.359 0.000) 21.997 2.649 ( 13.173 20.395 0.000) 24.279 3.115 ( -0.475 2.829 0.000) 2.868 3.129 ( 9.134 10.519 0.000) 13.931 3.496 ( 29.931 27.244 0.000) 40.474 3.576 ( 26.281 27.091 0.000) 37.744 6.322 ( -1.065 1.925 0.000) 2.200 6.417 ( -0.064 1.687 0.000) 1.688 6.636 ( 5.458 1.315 0.000) 5.614 6.662 ( 2.195 -0.650 0.000) 2.289 6.861 ( 6.276 8.550 0.000) 10.606 6.917 ( 0.092 6.132 0.000) 6.133 7.614 ( 5.957 6.934 0.000) 9.141 8.263 ( -6.718 -7.886 0.000) 10.359 8.708 ( -15.597 -11.411 0.000) 19.326 8.943 ( 8.589 1.627 0.000) 8.741 9.428 ( -32.493 -22.187 0.000) 39.345 9.929 ( 1.691 -12.454 0.000) 12.568 ======================= Grid point 18 (11/54) ======================= q-point: ( 0.20 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.110 ( 1.405 1.848 0.000) 2.321 1.113 ( 2.269 4.208 0.000) 4.781 1.169 ( 0.897 1.588 0.000) 1.824 1.241 ( 2.033 3.357 0.000) 3.925 1.374 ( -5.392 -6.744 0.000) 8.635 1.546 ( 4.934 3.395 0.000) 5.989 2.935 ( 12.267 18.934 0.000) 22.561 3.085 ( 11.852 19.982 0.000) 23.233 3.233 ( 7.329 7.118 0.000) 10.217 3.406 ( 9.499 9.429 0.000) 13.384 4.307 ( 26.659 21.394 0.000) 34.183 4.311 ( 30.926 22.824 0.000) 38.437 6.322 ( -2.315 3.868 0.000) 4.508 6.432 ( -0.604 2.219 0.000) 2.299 6.673 ( 0.742 -2.319 0.000) 2.435 6.724 ( 4.578 -2.496 0.000) 5.214 6.991 ( -0.769 9.419 0.000) 9.450 7.062 ( 7.054 6.711 0.000) 9.737 7.730 ( 0.838 3.563 0.000) 3.661 8.044 ( -8.374 -8.102 0.000) 11.652 8.291 ( -15.090 -11.086 0.000) 18.724 8.582 ( -36.695 -18.637 0.000) 41.156 9.106 ( 7.997 -1.254 0.000) 8.095 9.812 ( 0.087 -13.624 0.000) 13.624 ======================= Grid point 19 (12/54) ======================= q-point: ( 0.27 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.091 ( -2.938 -1.294 0.000) 3.210 1.190 ( 2.104 3.476 0.000) 4.063 1.219 ( 1.859 2.587 0.000) 3.186 1.252 ( -2.225 -4.647 0.000) 5.152 1.334 ( 3.438 5.783 0.000) 6.728 1.642 ( 3.236 1.204 0.000) 3.453 3.338 ( 11.790 17.073 0.000) 20.748 3.477 ( 9.928 17.489 0.000) 20.110 3.490 ( 11.047 9.310 0.000) 14.447 3.643 ( 6.565 8.444 0.000) 10.696 4.956 ( 23.396 17.106 0.000) 28.982 5.044 ( 27.087 17.203 0.000) 32.088 6.315 ( -4.049 5.428 0.000) 6.772 6.436 ( -1.685 2.126 0.000) 2.713 6.633 ( -1.202 -5.259 0.000) 5.394 6.745 ( 1.488 -3.643 0.000) 3.935 7.073 ( -1.950 10.061 0.000) 10.248 7.242 ( 6.348 4.077 0.000) 7.545 7.645 ( -11.393 -2.547 0.000) 11.674 7.670 ( -32.283 -8.109 0.000) 33.286 7.935 ( -9.161 -7.948 0.000) 12.129 8.078 ( 0.729 -2.119 0.000) 2.241 9.185 ( 4.617 -5.530 0.000) 7.204 9.640 ( -3.137 -15.160 0.000) 15.481 ======================= Grid point 20 (13/54) ======================= q-point: ( 0.33 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.023 ( -2.468 -1.312 0.000) 2.795 1.166 ( -2.492 -4.563 0.000) 5.199 1.251 ( 0.194 4.329 0.000) 4.333 1.274 ( 1.245 2.219 0.000) 2.544 1.486 ( 6.098 5.757 0.000) 8.387 1.690 ( 1.724 -0.042 0.000) 1.725 3.695 ( 11.570 10.516 0.000) 15.636 3.737 ( 4.516 4.755 0.000) 6.558 3.793 ( 9.067 13.257 0.000) 16.061 3.870 ( 5.232 17.805 0.000) 18.557 5.490 ( 18.315 13.417 0.000) 22.703 5.601 ( 18.409 10.056 0.000) 20.976 6.284 ( -5.720 5.324 0.000) 7.814 6.412 ( -3.069 1.487 0.000) 3.410 6.491 ( -6.057 -12.135 0.000) 13.562 6.711 ( -0.792 -3.596 0.000) 3.682 7.035 ( -28.170 3.835 0.000) 28.430 7.093 ( -8.879 10.489 0.000) 13.743 7.385 ( 3.079 3.872 0.000) 4.947 7.438 ( -4.415 -3.303 0.000) 5.514 7.880 ( 6.319 -1.383 0.000) 6.469 8.158 ( 6.003 1.386 0.000) 6.161 9.161 ( 0.821 -7.983 0.000) 8.025 9.376 ( -8.211 -16.516 0.000) 18.444 ======================= Grid point 21 (14/54) ======================= q-point: ( 0.40 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.980 ( -1.298 -0.728 0.000) 1.488 1.079 ( -1.722 -3.820 0.000) 4.190 1.286 ( -1.094 4.055 0.000) 4.200 1.309 ( 0.113 2.185 0.000) 2.188 1.627 ( 3.275 2.729 0.000) 4.263 1.719 ( 2.160 0.104 0.000) 2.163 3.807 ( 0.219 2.624 0.000) 2.633 3.946 ( 8.755 1.955 0.000) 8.971 4.057 ( 2.214 16.343 0.000) 16.492 4.121 ( 1.347 16.610 0.000) 16.664 5.870 ( 10.838 8.371 0.000) 13.695 5.896 ( 7.825 2.771 0.000) 8.301 6.190 ( -8.873 -11.376 0.000) 14.427 6.229 ( -6.611 4.523 0.000) 8.010 6.362 ( -3.533 1.152 0.000) 3.716 6.654 ( -1.006 -3.688 0.000) 3.823 6.761 ( -15.256 14.006 0.000) 20.710 7.001 ( -11.890 11.680 0.000) 16.668 7.419 ( 2.742 -1.662 0.000) 3.207 7.470 ( 2.889 1.751 0.000) 3.378 8.114 ( 16.971 4.568 0.000) 17.575 8.277 ( 5.159 1.030 0.000) 5.261 9.005 ( -13.395 -17.423 0.000) 21.977 9.071 ( -0.546 -8.223 0.000) 8.241 ======================= Grid point 22 (15/54) ======================= q-point: ( 0.47 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.962 ( -0.001 0.001 0.000) 0.002 1.043 ( 1.138 -1.971 0.000) 2.276 1.297 ( -2.066 3.578 0.000) 4.132 1.321 ( -0.997 1.727 0.000) 1.994 1.659 ( -0.187 0.323 0.000) 0.373 1.752 ( 0.555 -0.961 0.000) 1.110 3.820 ( -1.383 2.395 0.000) 2.765 4.036 ( 1.873 -3.244 0.000) 3.745 4.168 ( -7.288 12.623 0.000) 14.576 4.225 ( -6.894 11.941 0.000) 13.788 6.002 ( 0.926 -1.604 0.000) 1.853 6.004 ( 1.262 -2.186 0.000) 2.524 6.017 ( -0.652 1.129 0.000) 1.303 6.194 ( -3.029 5.247 0.000) 6.059 6.333 ( -1.414 2.450 0.000) 2.829 6.619 ( 1.520 -2.634 0.000) 3.041 6.704 ( -9.418 16.313 0.000) 18.836 6.952 ( -8.314 14.400 0.000) 16.628 7.442 ( 1.496 -2.591 0.000) 2.992 7.508 ( 0.224 -0.388 0.000) 0.448 8.335 ( 0.722 -1.251 0.000) 1.444 8.474 ( 1.082 -1.873 0.000) 2.163 8.616 ( 4.916 -8.514 0.000) 9.831 9.015 ( 3.346 -5.796 0.000) 6.692 ======================= Grid point 33 (16/54) ======================= q-point: ( 0.13 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.105 ( 1.912 3.311 0.000) 3.823 1.108 ( 1.285 2.226 0.000) 2.570 1.208 ( 1.874 3.245 0.000) 3.748 1.214 ( 2.453 4.249 0.000) 4.907 1.387 ( -5.007 -8.672 0.000) 10.014 1.516 ( 2.935 5.084 0.000) 5.870 2.931 ( 11.536 19.981 0.000) 23.072 3.099 ( 12.667 21.940 0.000) 25.334 3.217 ( 4.757 8.239 0.000) 9.513 3.369 ( 7.056 12.222 0.000) 14.113 4.129 ( 18.701 32.391 0.000) 37.401 4.164 ( 16.887 29.248 0.000) 33.773 6.377 ( 1.861 3.224 0.000) 3.722 6.456 ( 1.094 1.895 0.000) 2.188 6.644 ( -0.805 -1.394 0.000) 1.609 6.663 ( 0.841 1.456 0.000) 1.682 7.033 ( 4.614 7.991 0.000) 9.227 7.053 ( 3.687 6.386 0.000) 7.373 7.741 ( 3.191 5.527 0.000) 6.382 8.075 ( -5.790 -10.028 0.000) 11.579 8.395 ( -9.855 -17.070 0.000) 19.711 8.879 ( -17.336 -30.026 0.000) 34.671 8.999 ( 1.397 2.420 0.000) 2.794 9.686 ( -6.243 -10.813 0.000) 12.486 ======================= Grid point 34 (17/54) ======================= q-point: ( 0.20 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.125 ( -1.581 -0.402 0.000) 1.631 1.157 ( 2.767 2.294 0.000) 3.594 1.234 ( -2.776 -6.543 0.000) 7.108 1.265 ( 0.142 5.322 0.000) 5.324 1.335 ( 3.756 5.611 0.000) 6.752 1.603 ( 2.611 2.155 0.000) 3.385 3.287 ( 9.588 11.740 0.000) 15.157 3.484 ( 8.743 19.576 0.000) 21.440 3.505 ( 8.253 16.431 0.000) 18.387 3.639 ( 8.330 15.915 0.000) 17.963 4.765 ( 19.657 22.824 0.000) 30.122 4.796 ( 22.373 23.885 0.000) 32.727 6.442 ( -1.332 7.350 0.000) 7.470 6.487 ( -0.510 2.772 0.000) 2.819 6.595 ( -0.859 -5.358 0.000) 5.427 6.690 ( 1.515 -0.899 0.000) 1.762 7.168 ( 1.435 6.859 0.000) 7.008 7.183 ( 4.593 4.333 0.000) 6.314 7.822 ( -3.574 5.398 0.000) 6.474 7.866 ( -6.112 -8.781 0.000) 10.698 8.070 ( -6.907 -9.882 0.000) 12.057 8.246 ( -26.279 -13.609 0.000) 29.593 9.027 ( 3.707 -5.972 0.000) 7.029 9.474 ( -1.176 -18.132 0.000) 18.170 ======================= Grid point 35 (18/54) ======================= q-point: ( 0.27 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.058 ( -3.264 -1.822 0.000) 3.738 1.146 ( -1.285 -5.443 0.000) 5.593 1.249 ( 3.423 3.392 0.000) 4.819 1.320 ( -0.355 5.520 0.000) 5.531 1.464 ( 3.844 6.433 0.000) 7.494 1.655 ( 1.974 0.185 0.000) 1.982 3.555 ( 10.730 5.963 0.000) 12.275 3.722 ( 4.328 7.021 0.000) 8.247 3.873 ( 7.643 25.093 0.000) 26.231 3.994 ( 6.435 23.833 0.000) 24.687 5.319 ( 18.193 17.837 0.000) 25.478 5.385 ( 19.447 15.747 0.000) 25.023 6.419 ( -5.521 -3.418 0.000) 6.494 6.490 ( -2.549 2.914 0.000) 3.872 6.529 ( -4.407 2.959 0.000) 5.309 6.690 ( 0.005 -1.880 0.000) 1.880 7.249 ( 1.078 4.695 0.000) 4.817 7.258 ( 1.426 -2.313 0.000) 2.717 7.633 ( -16.288 3.210 0.000) 16.601 7.665 ( -17.820 5.422 0.000) 18.627 7.791 ( -6.136 -4.702 0.000) 7.731 8.061 ( 5.700 0.384 0.000) 5.713 9.008 ( 3.992 -10.965 0.000) 11.669 9.252 ( -1.190 -21.468 0.000) 21.501 ======================= Grid point 36 (19/54) ======================= q-point: ( 0.33 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.997 ( -1.726 -1.287 0.000) 2.153 1.062 ( -1.694 -5.028 0.000) 5.306 1.328 ( 1.191 3.607 0.000) 3.799 1.363 ( -1.218 5.452 0.000) 5.586 1.597 ( 3.917 4.537 0.000) 5.994 1.682 ( 1.596 -0.617 0.000) 1.711 3.789 ( 10.053 2.228 0.000) 10.297 3.833 ( 1.396 5.063 0.000) 5.252 4.246 ( 3.235 26.643 0.000) 26.839 4.324 ( 1.589 24.875 0.000) 24.926 5.760 ( 13.043 12.342 0.000) 17.957 5.779 ( 10.743 7.238 0.000) 12.954 6.218 ( -7.083 -10.251 0.000) 12.460 6.443 ( -10.494 5.885 0.000) 12.032 6.455 ( -4.323 2.410 0.000) 4.949 6.651 ( -1.740 -2.379 0.000) 2.947 7.191 ( -9.730 8.144 0.000) 12.689 7.308 ( 1.907 1.754 0.000) 2.591 7.420 ( -13.951 9.179 0.000) 16.699 7.495 ( -7.418 7.579 0.000) 10.605 7.858 ( 12.958 -0.317 0.000) 12.962 8.181 ( 6.270 1.018 0.000) 6.352 8.944 ( 2.536 -12.427 0.000) 12.683 8.972 ( -5.539 -21.651 0.000) 22.348 ======================= Grid point 37 (20/54) ======================= q-point: ( 0.40 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.966 ( -0.696 -0.654 0.000) 0.955 1.001 ( 0.209 -3.425 0.000) 3.431 1.366 ( -0.873 3.281 0.000) 3.395 1.387 ( -2.331 5.207 0.000) 5.705 1.663 ( 0.496 0.850 0.000) 0.985 1.716 ( 2.328 -0.545 0.000) 2.391 3.883 ( -1.286 4.540 0.000) 4.719 3.937 ( 5.583 -1.768 0.000) 5.856 4.496 ( -6.348 24.809 0.000) 25.608 4.534 ( -6.504 22.567 0.000) 23.486 5.938 ( 3.587 1.341 0.000) 3.829 6.014 ( 4.477 5.238 0.000) 6.891 6.022 ( -4.439 -5.092 0.000) 6.756 6.328 ( -6.599 4.385 0.000) 7.923 6.406 ( -3.370 2.938 0.000) 4.471 6.595 ( -1.155 -2.074 0.000) 2.374 7.131 ( -13.635 20.867 0.000) 24.927 7.318 ( -12.603 18.033 0.000) 22.001 7.354 ( 3.868 -4.784 0.000) 6.152 7.482 ( 0.915 -1.255 0.000) 1.553 8.189 ( 21.217 2.088 0.000) 21.320 8.284 ( 3.882 -0.264 0.000) 3.891 8.606 ( -11.166 -19.686 0.000) 22.632 8.864 ( 3.309 -11.280 0.000) 11.756 ======================= Grid point 49 (21/54) ======================= q-point: ( 0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.092 ( -1.761 -3.049 0.000) 3.521 1.114 ( -2.658 -4.603 0.000) 5.316 1.232 ( 2.727 4.723 0.000) 5.454 1.358 ( 2.207 3.823 0.000) 4.414 1.464 ( 3.707 6.420 0.000) 7.413 1.636 ( 0.661 1.145 0.000) 1.322 3.485 ( 5.147 8.914 0.000) 10.294 3.720 ( 4.376 7.580 0.000) 8.753 3.949 ( 13.571 23.505 0.000) 27.141 4.071 ( 12.915 22.369 0.000) 25.829 5.241 ( 13.133 22.747 0.000) 26.266 5.282 ( 12.956 22.440 0.000) 25.912 6.439 ( -5.164 -8.944 0.000) 10.328 6.536 ( 0.960 1.662 0.000) 1.920 6.611 ( 4.225 7.318 0.000) 8.450 6.677 ( -0.342 -0.593 0.000) 0.685 7.216 ( -0.770 -1.334 0.000) 1.541 7.251 ( 0.605 1.047 0.000) 1.209 7.704 ( -3.383 -5.859 0.000) 6.765 7.885 ( -2.984 -5.168 0.000) 5.968 7.964 ( 3.074 5.324 0.000) 6.148 8.030 ( -4.385 -7.595 0.000) 8.770 8.887 ( -3.960 -6.858 0.000) 7.919 9.103 ( -9.784 -16.947 0.000) 19.568 ======================= Grid point 50 (22/54) ======================= q-point: ( 0.27 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.014 ( -2.395 -2.696 0.000) 3.606 1.041 ( -1.409 -4.161 0.000) 4.393 1.334 ( 2.729 4.705 0.000) 5.439 1.426 ( 0.514 4.607 0.000) 4.635 1.585 ( 2.991 4.901 0.000) 5.742 1.654 ( 0.872 -0.268 0.000) 0.912 3.669 ( 6.939 5.442 0.000) 8.819 3.855 ( 2.839 6.112 0.000) 6.739 4.428 ( 9.771 27.875 0.000) 29.538 4.512 ( 8.218 25.443 0.000) 26.738 5.675 ( 12.160 16.034 0.000) 20.123 5.682 ( 10.526 12.414 0.000) 16.276 6.233 ( -6.138 -9.883 0.000) 11.634 6.545 ( -2.288 2.061 0.000) 3.079 6.655 ( -0.970 -1.755 0.000) 2.005 6.685 ( -6.700 7.198 0.000) 9.834 7.173 ( 3.241 -6.284 0.000) 7.071 7.234 ( 0.661 -4.424 0.000) 4.474 7.683 ( 4.980 0.752 0.000) 5.037 7.788 ( -9.815 9.655 0.000) 13.767 7.890 ( -13.418 9.246 0.000) 16.295 8.090 ( 5.475 2.410 0.000) 5.982 8.763 ( 1.130 -11.902 0.000) 11.956 8.791 ( -2.585 -22.092 0.000) 22.242 ======================= Grid point 51 (23/54) ======================= q-point: ( 0.33 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.962 ( -0.440 -2.499 0.000) 2.537 0.979 ( -1.398 -2.344 0.000) 2.729 1.412 ( 0.419 4.421 0.000) 4.441 1.471 ( -0.925 4.827 0.000) 4.915 1.659 ( 0.963 0.894 0.000) 1.314 1.669 ( 2.046 -0.213 0.000) 2.057 3.828 ( 5.933 1.791 0.000) 6.197 3.946 ( 0.449 5.580 0.000) 5.598 4.844 ( 1.547 30.647 0.000) 30.686 4.866 ( -0.064 26.948 0.000) 26.948 5.893 ( 4.206 3.572 0.000) 5.518 5.976 ( 6.938 7.834 0.000) 10.464 6.049 ( -4.563 -5.660 0.000) 7.271 6.501 ( -4.046 1.813 0.000) 4.434 6.534 ( -12.607 2.664 0.000) 12.886 6.605 ( -2.034 -2.178 0.000) 2.979 7.216 ( 3.921 -6.651 0.000) 7.720 7.245 ( 10.104 -4.794 0.000) 11.184 7.714 ( -16.996 12.543 0.000) 21.123 7.727 ( -11.665 15.746 0.000) 19.596 7.894 ( 12.783 4.138 0.000) 13.437 8.205 ( 4.970 1.410 0.000) 5.167 8.533 ( -3.096 -19.283 0.000) 19.530 8.675 ( 3.236 -12.910 0.000) 13.309 ======================= Grid point 52 (24/54) ======================= q-point: ( 0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.945 ( 1.088 -1.885 0.000) 2.176 0.953 ( 0.328 -0.568 0.000) 0.656 1.438 ( -2.073 3.591 0.000) 4.146 1.489 ( -2.555 4.425 0.000) 5.110 1.667 ( 0.241 -0.417 0.000) 0.481 1.692 ( 0.992 -1.718 0.000) 1.984 3.894 ( 1.186 -2.054 0.000) 2.372 3.980 ( -2.547 4.411 0.000) 5.093 5.016 ( -13.241 22.935 0.000) 26.483 5.037 ( -15.001 25.983 0.000) 30.002 5.956 ( -0.296 0.513 0.000) 0.593 5.961 ( 0.662 -1.146 0.000) 1.324 6.076 ( -0.352 0.609 0.000) 0.703 6.412 ( -2.455 4.252 0.000) 4.910 6.472 ( -1.790 3.100 0.000) 3.580 6.566 ( 0.431 -0.747 0.000) 0.862 7.221 ( 4.425 -7.665 0.000) 8.851 7.352 ( 6.676 -11.563 0.000) 13.352 7.586 ( -11.856 20.535 0.000) 23.712 7.703 ( -10.422 18.051 0.000) 20.843 8.124 ( 5.235 -9.068 0.000) 10.470 8.267 ( 0.755 -1.308 0.000) 1.511 8.325 ( 2.950 -5.110 0.000) 5.901 8.629 ( 6.299 -10.910 0.000) 12.598 ======================= Grid point 66 (25/54) ======================= q-point: ( 0.27 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.951 ( -1.472 -2.550 0.000) 2.944 0.982 ( -1.270 -2.200 0.000) 2.540 1.433 ( 2.623 4.543 0.000) 5.246 1.507 ( 1.833 3.175 0.000) 3.666 1.647 ( -0.157 -0.272 0.000) 0.314 1.654 ( 0.975 1.689 0.000) 1.950 3.781 ( 2.963 5.132 0.000) 5.926 3.972 ( 2.890 5.006 0.000) 5.781 5.001 ( 15.542 26.920 0.000) 31.085 5.019 ( 13.254 22.956 0.000) 26.508 5.871 ( 3.125 5.412 0.000) 6.249 5.949 ( 5.558 9.627 0.000) 11.117 6.062 ( -3.455 -5.984 0.000) 6.910 6.558 ( -0.675 -1.169 0.000) 1.350 6.612 ( -1.386 -2.400 0.000) 2.772 6.759 ( -1.533 -2.656 0.000) 3.067 7.010 ( -3.669 -6.355 0.000) 7.338 7.132 ( -2.521 -4.366 0.000) 5.042 7.794 ( 5.177 8.967 0.000) 10.354 7.956 ( 3.267 5.659 0.000) 6.535 8.026 ( 1.234 2.138 0.000) 2.469 8.158 ( 2.293 3.971 0.000) 4.585 8.395 ( -8.323 -14.416 0.000) 16.647 8.549 ( -4.751 -8.229 0.000) 9.502 ======================= Grid point 67 (26/54) ======================= q-point: ( 0.33 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.921 ( 0.051 -0.977 0.000) 0.978 0.944 ( -0.566 -1.454 0.000) 1.560 1.505 ( -0.276 4.480 0.000) 4.488 1.557 ( -0.085 3.446 0.000) 3.447 1.645 ( 0.800 -1.206 0.000) 1.447 1.660 ( 0.184 -1.630 0.000) 1.640 3.862 ( 1.790 1.384 0.000) 2.263 4.048 ( 0.103 3.855 0.000) 3.857 5.410 ( -0.040 24.917 0.000) 24.917 5.471 ( 1.303 28.010 0.000) 28.040 5.920 ( 1.340 -1.268 0.000) 1.845 6.005 ( -1.350 2.996 0.000) 3.286 6.048 ( 1.569 -2.073 0.000) 2.600 6.520 ( -1.400 0.136 0.000) 1.406 6.542 ( -12.199 -3.180 0.000) 12.606 6.561 ( -1.255 -1.974 0.000) 2.339 7.001 ( 11.642 -16.248 0.000) 19.988 7.077 ( 2.146 -5.835 0.000) 6.217 7.945 ( -6.274 8.505 0.000) 10.569 8.001 ( 1.526 4.849 0.000) 5.083 8.020 ( -5.784 12.205 0.000) 13.507 8.233 ( 1.372 -9.036 0.000) 9.139 8.236 ( 2.482 1.371 0.000) 2.836 8.432 ( 2.396 -10.419 0.000) 10.691 ======================= Grid point 82 (27/54) ======================= q-point: ( 0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 54 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.920 ( -0.368 -0.638 0.000) 0.737 0.920 ( 0.368 0.638 0.000) 0.737 1.577 ( -0.000 -0.000 0.000) 0.000 1.614 ( -0.000 -0.000 0.000) 0.000 1.615 ( -0.281 -0.486 0.000) 0.561 1.615 ( 0.281 0.486 0.000) 0.562 3.878 ( -0.000 -0.000 0.000) 0.000 4.090 ( -0.000 -0.000 0.000) 0.000 5.839 ( -3.689 -6.390 0.000) 7.379 5.839 ( 3.689 6.390 0.000) 7.379 5.853 ( -8.289 -14.357 0.000) 16.578 5.853 ( 8.289 14.357 0.000) 16.578 6.280 ( -10.387 -17.990 0.000) 20.774 6.280 ( 10.387 17.990 0.000) 20.773 6.528 ( -0.514 -0.891 0.000) 1.029 6.528 ( 0.514 0.891 0.000) 1.029 6.801 ( 0.000 0.000 0.000) 0.000 7.012 ( -0.000 -0.000 0.000) 0.000 8.043 ( -0.596 -1.033 0.000) 1.192 8.043 ( 0.596 1.033 0.000) 1.192 8.151 ( 0.000 0.000 0.000) 0.000 8.235 ( -4.537 -7.857 0.000) 9.073 8.235 ( 4.537 7.857 0.000) 9.073 8.253 ( -0.000 -0.000 0.000) 0.000 ======================= Grid point 241 (28/54) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 54 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.516 ( -0.000 -0.000 20.971) 20.971 0.516 ( 0.000 0.000 20.971) 20.971 0.912 ( 0.000 -0.000 -13.124) 13.124 0.912 ( 0.000 0.000 -13.124) 13.124 0.989 ( 0.000 0.000 42.268) 42.268 1.800 ( 0.000 0.000 2.124) 2.124 1.800 ( 0.000 -0.000 2.124) 2.124 1.842 ( 0.000 -0.000 -1.331) 1.331 1.842 ( 0.000 0.000 -1.331) 1.331 1.871 ( 0.000 -0.000 -33.911) 33.911 3.816 ( 0.000 0.000 18.472) 18.472 4.176 ( 0.000 -0.000 -11.729) 11.729 6.334 ( 0.000 0.000 0.925) 0.925 6.334 ( 0.000 -0.000 0.925) 0.925 6.368 ( 0.000 -0.000 -1.928) 1.928 6.368 ( 0.000 0.000 -1.928) 1.928 6.705 ( -0.000 -0.000 3.846) 3.846 6.705 ( 0.000 -0.000 3.846) 3.846 6.788 ( -0.000 0.000 -2.887) 2.887 6.788 ( -0.000 0.000 -2.887) 2.887 8.596 ( 0.000 0.000 1.130) 1.130 8.624 ( 0.000 -0.000 -1.134) 1.134 10.843 ( 0.000 0.000 1.432) 1.432 10.879 ( 0.000 0.000 -1.405) 1.405 ======================= Grid point 242 (29/54) ======================= q-point: ( 0.07 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.685 ( 10.760 6.213 13.468) 18.324 0.724 ( 10.174 5.874 14.262) 18.477 0.966 ( 3.629 2.095 -9.667) 10.536 1.004 ( 5.467 3.156 -9.482) 11.391 1.157 ( 11.391 6.576 27.024) 30.055 1.732 ( -3.982 -2.299 -18.116) 18.691 1.939 ( 11.234 6.486 0.473) 12.981 1.949 ( 8.926 5.153 -0.302) 10.311 2.169 ( 25.802 14.897 -0.337) 29.795 2.240 ( 27.206 15.708 -8.113) 32.446 3.360 ( -27.801 -16.051 21.100) 38.415 3.500 ( -18.640 -10.762 4.488) 21.987 6.329 ( -0.372 -0.215 1.018) 1.105 6.365 ( -0.188 -0.109 -1.989) 2.001 6.464 ( 7.630 4.405 1.173) 8.888 6.468 ( 4.879 2.817 0.928) 5.710 6.714 ( 0.643 0.371 3.819) 3.891 6.796 ( 0.604 0.348 -2.890) 2.972 7.182 ( 25.712 14.845 12.118) 32.068 7.795 ( 40.568 23.422 -47.577) 66.767 8.630 ( 4.534 2.618 -0.567) 5.266 8.654 ( 2.888 1.667 -1.429) 3.628 10.537 ( -22.556 -13.022 8.196) 27.304 10.597 ( -19.041 -10.993 0.087) 21.986 ======================= Grid point 243 (30/54) ======================= q-point: ( 0.13 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.916 ( 7.801 4.504 6.185) 10.927 0.928 ( 7.415 4.281 6.949) 11.026 1.052 ( 3.361 1.940 -4.827) 6.194 1.104 ( 2.601 1.502 -7.087) 7.697 1.381 ( 6.515 3.761 10.650) 13.039 1.626 ( -4.716 -2.723 -9.205) 10.695 2.251 ( 15.666 9.044 1.633) 18.163 2.294 ( 17.360 10.023 -1.786) 20.125 2.924 ( 30.751 17.754 2.942) 35.630 3.017 ( 32.150 18.562 -4.222) 37.363 3.093 ( 4.782 2.761 7.187) 9.063 3.200 ( -1.388 -0.801 -1.844) 2.443 6.319 ( -0.381 -0.220 1.342) 1.413 6.363 ( 0.061 0.035 -2.247) 2.249 6.615 ( 7.302 4.216 2.281) 8.735 6.666 ( 9.085 5.245 -1.851) 10.652 6.731 ( 0.728 0.420 3.673) 3.768 6.811 ( 0.521 0.301 -2.800) 2.864 7.662 ( 11.901 6.871 17.492) 22.245 8.163 ( -0.231 -0.133 -21.624) 21.626 8.750 ( 1.916 1.106 -1.662) 2.767 8.788 ( 4.094 2.364 -6.854) 8.327 9.959 ( -25.891 -14.948 8.462) 31.071 10.102 ( -17.316 -9.997 -3.997) 20.390 ======================= Grid point 244 (31/54) ======================= q-point: ( 0.20 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.058 ( 4.477 2.585 2.924) 5.939 1.075 ( 4.505 2.601 1.874) 5.529 1.125 ( 2.751 1.588 -2.404) 3.983 1.145 ( 1.500 0.866 -4.088) 4.440 1.488 ( 3.123 1.803 1.639) 3.961 1.512 ( -4.577 -2.643 -0.013) 5.285 2.656 ( 17.621 10.174 2.681) 20.523 2.724 ( 18.119 10.461 -2.747) 21.102 3.161 ( 6.265 3.617 2.557) 7.672 3.214 ( 9.635 5.563 -1.682) 11.252 3.857 ( 35.314 20.389 0.166) 40.777 3.872 ( 34.163 19.724 -0.973) 39.460 6.310 ( -0.438 -0.253 1.889) 1.956 6.366 ( 0.169 0.097 -2.727) 2.734 6.743 ( 0.181 0.105 3.385) 3.392 6.801 ( 7.944 4.587 3.656) 9.875 6.816 ( -0.135 -0.078 -2.555) 2.560 6.888 ( 9.302 5.370 -3.357) 11.253 7.790 ( 0.552 0.319 10.283) 10.303 8.059 ( -6.847 -3.953 -10.897) 13.463 8.544 ( -16.767 -9.680 3.463) 19.668 8.696 ( -17.061 -9.850 -10.455) 22.303 9.474 ( -9.412 -5.434 18.637) 21.575 9.840 ( -6.243 -3.605 -12.615) 14.530 ======================= Grid point 245 (32/54) ======================= q-point: ( 0.27 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.127 ( -0.038 -0.022 1.875) 1.876 1.144 ( 2.914 1.682 1.743) 3.789 1.184 ( 2.271 1.311 -1.521) 3.031 1.189 ( 1.475 0.852 -3.089) 3.527 1.431 ( -1.268 -0.732 6.189) 6.360 1.555 ( 2.747 1.586 -4.104) 5.186 3.073 ( 17.088 9.866 2.804) 19.929 3.143 ( 16.892 9.753 -2.853) 19.713 3.380 ( 11.152 6.438 3.461) 13.334 3.465 ( 10.474 6.047 -3.274) 12.530 4.624 ( 28.964 16.723 1.126) 33.464 4.657 ( 30.834 17.802 -1.406) 35.632 6.295 ( -0.880 -0.508 2.486) 2.685 6.366 ( -0.314 -0.181 -3.289) 3.309 6.735 ( -1.016 -0.586 3.117) 3.330 6.801 ( -1.236 -0.714 -2.285) 2.694 6.978 ( 6.360 3.672 5.937) 9.443 7.102 ( 8.364 4.829 -3.797) 10.378 7.700 ( -9.519 -5.496 1.442) 11.085 7.804 ( -18.634 -10.758 -7.150) 22.673 8.141 ( -14.249 -8.227 1.194) 16.497 8.169 ( -14.943 -8.628 -0.648) 17.267 9.409 ( -0.649 -0.375 13.778) 13.798 9.714 ( -5.631 -3.251 -11.037) 12.809 ======================= Grid point 246 (33/54) ======================= q-point: ( 0.33 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.092 ( -2.399 -1.385 3.205) 4.236 1.179 ( -2.256 -1.302 -3.878) 4.672 1.203 ( 2.134 1.232 1.216) 2.748 1.233 ( 1.760 1.016 -1.118) 2.319 1.477 ( 4.550 2.627 6.798) 8.592 1.621 ( 2.650 1.530 -4.969) 5.836 3.465 ( 15.672 9.048 2.404) 18.256 3.526 ( 15.091 8.713 -2.470) 17.600 3.652 ( 10.901 6.294 1.583) 12.687 3.684 ( 7.534 4.350 -0.872) 8.743 5.245 ( 22.959 13.256 2.481) 26.627 5.313 ( 23.764 13.720 -3.000) 27.604 6.266 ( -1.635 -0.944 2.921) 3.478 6.347 ( -1.288 -0.744 -3.782) 4.064 6.693 ( -2.392 -1.381 3.079) 4.136 6.757 ( -2.345 -1.354 -2.162) 3.466 7.016 ( -6.767 -3.907 4.027) 8.791 7.126 ( -21.500 -12.413 -4.719) 25.271 7.337 ( 0.137 0.079 3.598) 3.601 7.453 ( -4.358 -2.516 -5.103) 7.167 7.969 ( -0.940 -0.542 4.844) 4.964 8.076 ( 2.107 1.216 -3.775) 4.491 9.355 ( -4.268 -2.464 7.833) 9.255 9.538 ( -9.416 -5.436 -6.884) 12.869 ======================= Grid point 247 (34/54) ======================= q-point: ( 0.40 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.033 ( -2.235 -1.290 2.881) 3.868 1.106 ( -3.257 -1.880 -2.984) 4.801 1.245 ( 1.342 0.775 0.945) 1.815 1.268 ( 1.222 0.705 -0.899) 1.673 1.601 ( 4.903 2.831 3.406) 6.607 1.678 ( 2.194 1.267 -2.846) 3.811 3.805 ( 12.297 7.100 1.788) 14.312 3.814 ( 3.594 2.075 2.978) 5.108 3.851 ( 11.644 6.722 -1.850) 13.571 3.877 ( 7.358 4.248 -2.019) 8.733 5.709 ( 15.588 9.000 2.167) 18.130 5.767 ( 13.871 8.008 -2.456) 16.204 6.217 ( -2.334 -1.348 3.039) 4.062 6.303 ( -2.304 -1.330 -4.061) 4.856 6.569 ( -19.498 -11.257 5.806) 23.251 6.630 ( -2.573 -1.485 3.387) 4.505 6.699 ( -2.214 -1.278 -2.309) 3.445 6.751 ( -11.323 -6.538 -9.247) 16.014 7.435 ( 4.002 2.310 2.510) 5.258 7.482 ( 2.904 1.677 -0.080) 3.355 8.066 ( 8.961 5.174 5.286) 11.620 8.191 ( 6.879 3.972 -4.248) 9.008 9.206 ( -7.854 -4.535 1.885) 9.263 9.259 ( -13.855 -7.999 -2.290) 16.161 ======================= Grid point 248 (35/54) ======================= q-point: ( 0.47 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.993 ( -1.003 -0.579 2.091) 2.390 1.046 ( -1.403 -0.810 -2.199) 2.732 1.267 ( 0.442 0.255 0.916) 1.049 1.289 ( 0.468 0.270 -0.885) 1.037 1.677 ( 1.324 0.765 2.128) 2.620 1.727 ( 1.568 0.906 -1.938) 2.653 3.864 ( 0.874 0.504 5.508) 5.600 3.994 ( 2.411 1.392 -4.870) 5.609 4.020 ( 4.966 2.867 1.304) 5.880 4.053 ( 4.659 2.690 -1.346) 5.545 5.954 ( 4.326 2.498 0.562) 5.027 5.982 ( 4.677 2.700 -1.260) 5.546 6.168 ( -1.376 -0.794 2.950) 3.350 6.253 ( -1.417 -0.818 -4.120) 4.433 6.254 ( -6.354 -3.668 11.400) 13.557 6.568 ( -3.742 -2.160 -14.492) 15.122 6.585 ( -1.046 -0.604 3.807) 3.994 6.662 ( -0.803 -0.464 -2.597) 2.757 7.496 ( 1.209 0.698 2.603) 2.954 7.534 ( 1.211 0.699 0.371) 1.447 8.336 ( 8.531 4.925 2.018) 10.055 8.375 ( 9.594 5.539 -0.110) 11.078 8.907 ( -13.777 -7.954 3.996) 16.403 9.024 ( -5.731 -3.309 -5.273) 8.462 ======================= Grid point 257 (36/54) ======================= q-point: ( 0.07 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.874 ( 5.597 9.694 7.442) 13.441 0.875 ( 4.472 7.746 8.816) 12.559 1.038 ( 2.374 4.112 -5.645) 7.376 1.078 ( 1.765 3.056 -7.384) 8.184 1.331 ( 4.521 7.830 13.832) 16.525 1.652 ( -2.811 -4.869 -12.106) 13.347 2.170 ( 9.032 15.644 1.276) 18.109 2.203 ( 10.288 17.819 -1.367) 20.621 2.705 ( 19.698 34.118 1.672) 39.431 2.785 ( 19.638 34.014 -4.699) 39.556 3.097 ( -2.615 -4.529 10.859) 12.052 3.232 ( -6.141 -10.636 -1.008) 12.322 6.330 ( 0.059 0.101 0.948) 0.955 6.367 ( 0.174 0.302 -2.200) 2.228 6.552 ( 3.027 5.242 1.934) 6.355 6.586 ( 3.809 6.598 -0.666) 7.648 6.744 ( 1.570 2.719 3.493) 4.696 6.824 ( 1.305 2.261 -2.950) 3.939 7.571 ( 9.707 16.814 18.456) 26.787 8.151 ( 2.962 5.131 -27.719) 28.345 8.731 ( 3.252 5.633 -1.554) 6.687 8.754 ( 3.788 6.561 -7.575) 10.713 10.153 ( -11.876 -20.571 5.109) 24.296 10.184 ( -15.232 -26.383 0.868) 30.477 ======================= Grid point 258 (37/54) ======================= q-point: ( 0.13 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.021 ( 4.109 4.855 2.408) 6.801 1.040 ( 3.866 5.519 2.921) 7.344 1.110 ( 0.985 3.023 -2.794) 4.233 1.134 ( 2.435 2.214 -3.940) 5.134 1.450 ( 3.117 2.757 4.396) 6.053 1.544 ( -3.481 -4.891 -2.992) 6.707 2.552 ( 12.366 18.985 2.634) 22.810 2.618 ( 12.854 20.025 -2.601) 23.937 3.130 ( 2.423 4.362 1.039) 5.097 3.140 ( 4.435 5.071 0.388) 6.748 3.510 ( 29.033 28.085 1.145) 40.410 3.548 ( 27.516 28.257 -1.974) 39.490 6.332 ( -0.627 1.361 1.096) 1.856 6.376 ( -0.069 1.139 -2.566) 2.808 6.668 ( 5.251 1.422 2.314) 5.912 6.695 ( 3.577 0.054 0.990) 3.712 6.847 ( 4.072 7.919 0.061) 8.905 6.887 ( 1.050 6.294 -2.370) 6.807 7.783 ( 2.740 4.118 13.477) 14.356 8.117 ( -4.119 -4.158 -12.874) 14.142 8.705 ( -10.220 -7.118 1.656) 12.564 8.813 ( -2.505 -1.373 -7.214) 7.759 9.599 ( -17.741 -18.542 10.994) 27.918 9.828 ( -4.151 -14.339 -8.177) 17.021 ======================= Grid point 259 (38/54) ======================= q-point: ( 0.20 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.109 ( 1.826 1.656 0.584) 2.533 1.130 ( 2.789 1.309 -0.271) 3.093 1.159 ( 0.626 4.042 0.680) 4.146 1.202 ( 1.274 3.396 -2.264) 4.276 1.444 ( -2.111 -3.245 4.150) 5.676 1.520 ( 2.715 1.402 -2.299) 3.824 2.973 ( 12.089 19.286 3.032) 22.963 3.048 ( 11.846 19.782 -2.990) 23.251 3.274 ( 7.986 7.340 3.376) 11.360 3.361 ( 9.146 8.462 -3.566) 12.960 4.305 ( 28.327 22.409 0.228) 36.119 4.313 ( 29.126 22.515 -0.222) 36.814 6.334 ( -1.526 2.578 1.295) 3.263 6.385 ( -0.603 1.635 -2.967) 3.441 6.719 ( 0.912 -1.609 2.755) 3.318 6.747 ( 3.119 -2.169 0.771) 3.877 6.991 ( 2.371 9.160 0.795) 9.495 7.036 ( 6.414 7.712 -2.226) 10.274 7.799 ( -3.657 0.741 5.954) 7.027 7.964 ( -7.934 -4.950 -6.897) 11.620 8.320 ( -17.401 -11.662 3.122) 21.179 8.441 ( -25.322 -13.731 -7.489) 29.763 9.349 ( 1.431 -7.347 15.731) 17.421 9.675 ( 0.192 -11.855 -11.412) 16.457 ======================= Grid point 260 (39/54) ======================= q-point: ( 0.27 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.112 ( -1.513 -1.269 2.045) 2.843 1.172 ( 0.850 -1.314 -2.591) 3.028 1.211 ( 0.512 4.117 0.922) 4.250 1.237 ( 0.129 2.238 -1.034) 2.468 1.439 ( 2.554 2.128 7.021) 7.768 1.585 ( 2.916 1.364 -4.915) 5.875 3.373 ( 11.244 17.243 2.730) 20.765 3.440 ( 10.325 16.976 -2.691) 20.051 3.532 ( 10.317 9.473 3.045) 14.334 3.607 ( 8.089 8.539 -2.950) 12.126 4.977 ( 24.109 17.350 1.763) 29.755 5.022 ( 25.913 17.398 -1.795) 31.264 6.327 ( -2.908 3.467 1.373) 4.729 6.383 ( -1.729 1.635 -3.278) 4.051 6.687 ( -1.533 -3.895 3.353) 5.363 6.742 ( -0.409 -3.442 -1.014) 3.611 7.128 ( 0.872 9.046 4.583) 10.178 7.223 ( 5.910 5.956 -2.422) 8.733 7.593 ( -20.999 -5.059 -2.499) 21.744 7.606 ( -22.953 -5.189 -3.542) 23.798 7.990 ( -6.255 -6.150 3.529) 9.455 8.056 ( -1.717 -3.591 -1.979) 4.445 9.315 ( 1.921 -8.527 9.813) 13.142 9.539 ( -1.730 -13.180 -8.270) 15.656 ======================= Grid point 261 (40/54) ======================= q-point: ( 0.33 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.056 ( -2.235 -2.027 2.668) 4.028 1.126 ( -2.102 -3.546 -2.984) 5.089 1.260 ( 0.121 3.942 0.656) 3.998 1.273 ( 0.510 2.826 -0.308) 2.888 1.546 ( 4.765 3.799 4.515) 7.584 1.647 ( 2.552 1.058 -3.637) 4.567 3.710 ( 10.247 8.955 1.017) 13.647 3.728 ( 7.417 5.873 -0.562) 9.477 3.818 ( 7.329 14.262 1.566) 16.112 3.852 ( 5.830 17.053 -1.367) 18.074 5.517 ( 18.280 12.942 2.224) 22.508 5.574 ( 18.322 11.329 -2.265) 21.661 6.299 ( -4.419 3.791 1.479) 6.007 6.353 ( -3.177 0.996 -3.147) 4.582 6.562 ( -5.658 -9.452 4.336) 11.838 6.660 ( -3.027 -6.298 -3.880) 7.993 7.041 ( -23.324 5.649 1.210) 24.029 7.086 ( -19.892 8.300 -1.794) 21.629 7.367 ( 5.302 1.115 2.087) 5.806 7.442 ( 1.984 0.932 -2.472) 3.304 7.960 ( 5.460 -0.166 6.323) 8.356 8.102 ( 6.000 1.304 -4.896) 7.853 9.217 ( -1.966 -10.571 4.405) 11.620 9.324 ( -6.390 -14.756 -4.203) 16.620 ======================= Grid point 262 (41/54) ======================= q-point: ( 0.40 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.003 ( -1.386 -1.527 1.904) 2.807 1.053 ( -1.588 -3.062 -2.049) 4.012 1.293 ( -0.862 3.632 0.517) 3.768 1.304 ( -0.279 2.705 -0.386) 2.746 1.652 ( 2.890 1.943 1.975) 4.003 1.698 ( 2.329 0.655 -1.727) 2.972 3.845 ( 2.466 1.925 2.976) 4.318 3.915 ( 6.690 1.363 -2.623) 7.314 4.075 ( 1.761 17.091 1.457) 17.243 4.109 ( 1.544 16.858 -1.152) 16.968 5.873 ( 9.218 6.082 0.677) 11.064 5.894 ( 8.746 4.590 -0.675) 9.900 6.251 ( -5.050 3.759 2.122) 6.643 6.279 ( -4.688 -4.348 4.549) 7.847 6.337 ( -6.833 -2.577 0.971) 7.367 6.533 ( -2.365 -6.267 -9.699) 11.787 6.804 ( -14.126 13.171 3.306) 19.594 6.903 ( -13.101 12.496 -5.453) 18.908 7.453 ( 3.537 -0.631 2.964) 4.658 7.505 ( 2.591 0.480 -0.178) 2.641 8.169 ( 14.041 4.765 4.652) 15.540 8.262 ( 8.523 2.779 -2.230) 9.238 9.009 ( -10.729 -15.639 0.635) 18.976 9.042 ( -4.156 -10.967 -2.014) 11.900 ======================= Grid point 263 (42/54) ======================= q-point: ( 0.47 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.981 ( 0.290 -0.503 1.543) 1.648 1.021 ( 0.857 -1.485 -1.683) 2.403 1.304 ( -1.824 3.159 0.536) 3.686 1.316 ( -1.295 2.243 -0.439) 2.627 1.684 ( 0.033 -0.057 1.901) 1.902 1.730 ( 0.395 -0.684 -1.802) 1.967 3.875 ( -0.441 0.764 4.333) 4.422 3.982 ( 1.147 -1.986 -4.281) 4.857 4.186 ( -7.296 12.636 1.356) 14.654 4.215 ( -7.061 12.230 -0.928) 14.152 5.972 ( 0.301 -0.522 0.319) 0.682 5.992 ( 0.644 -1.115 -0.902) 1.572 6.193 ( 0.763 -1.322 11.398) 11.500 6.217 ( -2.671 4.627 2.136) 5.754 6.285 ( -2.045 3.541 -3.339) 5.279 6.479 ( 1.555 -2.694 -11.609) 12.019 6.752 ( -8.766 15.183 3.664) 17.911 6.854 ( -8.555 14.818 -5.460) 17.961 7.481 ( 1.228 -2.127 2.933) 3.826 7.531 ( 0.603 -1.045 -0.318) 1.248 8.370 ( 0.675 -1.169 3.146) 3.423 8.444 ( 0.927 -1.606 -2.605) 3.197 8.723 ( 4.424 -7.663 8.197) 12.062 8.920 ( 3.635 -6.296 -7.666) 10.565 ======================= Grid point 274 (43/54) ======================= q-point: ( 0.13 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.105 ( 1.690 2.927 -0.085) 3.381 1.106 ( 1.048 1.816 -0.220) 2.108 1.193 ( 2.377 4.117 -0.528) 4.783 1.199 ( 2.020 3.499 -0.636) 4.090 1.444 ( -2.332 -4.039 3.326) 5.729 1.500 ( 1.333 2.308 -1.513) 3.065 2.975 ( 11.873 20.565 3.542) 24.009 3.060 ( 12.449 21.563 -3.235) 25.108 3.248 ( 5.120 8.869 2.622) 10.571 3.323 ( 6.292 10.898 -3.403) 13.035 4.142 ( 18.367 31.812 0.886) 36.744 4.159 ( 17.479 30.274 -0.498) 34.961 6.369 ( 1.178 2.040 0.047) 2.356 6.402 ( 0.713 1.234 -2.843) 3.181 6.698 ( 0.038 0.067 2.496) 2.497 6.700 ( 0.948 1.643 2.026) 2.776 7.015 ( 4.687 8.118 -0.808) 9.408 7.027 ( 3.978 6.890 -1.642) 8.124 7.836 ( 0.802 1.389 7.327) 7.500 8.005 ( -3.618 -6.267 -6.134) 9.487 8.445 ( -9.866 -17.089 5.040) 20.366 8.631 ( -10.267 -17.783 -11.059) 23.323 9.286 ( -5.565 -9.638 14.412) 18.209 9.570 ( -6.003 -10.397 -9.667) 15.413 ======================= Grid point 275 (44/54) ======================= q-point: ( 0.20 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.124 ( -0.873 -0.297 0.358) 0.989 1.143 ( 1.107 -0.437 -0.993) 1.550 1.242 ( 0.289 1.966 0.816) 2.148 1.262 ( 0.729 3.667 -0.688) 3.801 1.423 ( 1.790 2.171 6.315) 6.913 1.553 ( 2.189 1.942 -4.307) 5.207 3.342 ( 9.181 12.482 4.420) 16.113 3.451 ( 8.221 14.994 -4.167) 17.600 3.519 ( 9.119 18.929 2.649) 21.177 3.594 ( 8.805 17.253 -3.357) 19.659 4.777 ( 20.206 23.270 0.906) 30.832 4.793 ( 21.594 23.728 -0.355) 32.085 6.406 ( -0.922 3.667 -1.408) 4.034 6.422 ( -0.117 1.702 -2.786) 3.266 6.681 ( -1.171 -2.225 4.237) 4.927 6.723 ( 0.229 -0.324 0.871) 0.957 7.160 ( 3.328 6.525 0.063) 7.325 7.176 ( 4.926 5.223 -0.805) 7.225 7.831 ( -6.277 2.410 1.308) 6.850 7.867 ( -6.059 -4.144 -0.856) 7.390 8.081 ( -10.158 -10.680 1.713) 14.838 8.164 ( -19.773 -12.273 -5.112) 23.828 9.165 ( 1.204 -10.165 9.919) 14.254 9.381 ( -0.672 -15.788 -7.719) 17.586 ======================= Grid point 276 (45/54) ======================= q-point: ( 0.27 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.074 ( -2.399 -2.352 1.390) 3.636 1.115 ( -1.219 -3.859 -1.976) 4.503 1.286 ( 1.395 3.428 1.949) 4.182 1.314 ( 0.031 4.822 -0.618) 4.862 1.515 ( 3.525 4.567 4.193) 7.132 1.613 ( 2.420 1.398 -3.563) 4.529 3.597 ( 9.256 6.231 3.340) 11.647 3.680 ( 5.977 6.911 -3.298) 9.714 3.904 ( 7.552 24.579 2.392) 25.824 3.963 ( 6.911 24.106 -2.396) 25.191 5.339 ( 18.383 17.508 1.587) 25.436 5.374 ( 19.040 16.447 -1.097) 25.184 6.392 ( -2.686 -0.166 -0.164) 2.696 6.437 ( -2.378 4.742 -3.638) 6.432 6.598 ( -4.921 -3.404 4.537) 7.509 6.679 ( -4.016 -2.565 -1.988) 5.163 7.268 ( 1.744 1.189 1.396) 2.531 7.279 ( 4.465 -0.415 0.927) 4.579 7.585 ( -20.300 5.086 -1.276) 20.967 7.618 ( -18.904 6.571 -3.183) 20.265 7.891 ( -0.911 -2.767 6.150) 6.805 8.013 ( 4.202 -0.363 -4.086) 5.872 9.075 ( 2.191 -14.085 5.163) 15.161 9.196 ( -0.294 -19.229 -4.581) 19.769 ======================= Grid point 277 (46/54) ======================= q-point: ( 0.33 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.012 ( -1.682 -2.179 1.209) 3.006 1.044 ( -1.650 -4.022 -1.398) 4.566 1.342 ( 0.265 3.999 0.911) 4.110 1.358 ( -0.851 4.963 -0.459) 5.056 1.621 ( 3.287 3.120 1.878) 4.906 1.664 ( 2.088 0.530 -1.565) 2.663 3.801 ( 7.877 2.864 0.942) 8.434 3.823 ( 3.553 4.316 -0.832) 5.652 4.268 ( 2.975 26.095 1.682) 26.318 4.306 ( 2.174 25.250 -1.424) 25.384 5.765 ( 11.866 10.189 0.568) 15.651 5.779 ( 11.370 8.584 -0.334) 14.251 6.298 ( -4.322 -5.405 6.546) 9.526 6.425 ( -5.740 6.118 -2.196) 8.672 6.459 ( -7.785 0.922 0.708) 7.871 6.540 ( -5.179 -5.026 -7.327) 10.284 7.202 ( -9.579 6.472 1.932) 11.720 7.292 ( -8.881 4.613 -5.703) 11.518 7.396 ( -4.178 7.426 5.389) 10.082 7.503 ( -4.265 6.021 -2.851) 7.910 7.976 ( 10.177 2.037 8.348) 13.320 8.135 ( 7.770 2.025 -4.550) 9.229 8.942 ( -0.073 -15.320 0.143) 15.321 8.958 ( -4.084 -19.863 -1.016) 20.304 ======================= Grid point 278 (47/54) ======================= q-point: ( 0.40 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.973 ( -0.471 -1.315 0.636) 1.535 0.991 ( -0.020 -2.694 -0.761) 2.800 1.373 ( -1.285 3.781 0.515) 4.026 1.383 ( -2.011 4.748 -0.308) 5.166 1.676 ( 0.904 0.433 1.084) 1.476 1.703 ( 1.818 -0.264 -1.044) 2.113 3.896 ( 0.375 2.914 1.038) 3.116 3.923 ( 3.776 -0.228 -1.123) 3.946 4.510 ( -6.294 24.186 0.976) 25.010 4.528 ( -6.337 23.045 -0.522) 23.906 5.941 ( 3.237 1.425 0.592) 3.586 5.963 ( 2.631 1.120 -1.273) 3.130 6.190 ( -1.779 -1.596 10.277) 10.551 6.334 ( -6.048 2.485 0.278) 6.545 6.385 ( -4.204 4.277 -1.428) 6.164 6.466 ( -1.231 0.276 -9.813) 9.893 7.134 ( -11.200 16.456 1.140) 19.938 7.208 ( -10.630 14.596 -5.404) 18.848 7.420 ( 1.575 -1.014 4.687) 5.048 7.504 ( -0.035 0.453 -1.446) 1.516 8.259 ( 16.193 3.702 4.256) 17.147 8.300 ( 7.429 0.631 0.520) 7.474 8.643 ( -8.034 -18.771 3.952) 20.797 8.780 ( 0.448 -13.446 -6.532) 14.955 ======================= Grid point 290 (48/54) ======================= q-point: ( 0.20 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.091 ( -1.719 -2.978 -0.007) 3.439 1.098 ( -1.833 -3.175 -0.689) 3.730 1.290 ( 1.953 3.382 3.287) 5.105 1.343 ( 2.121 3.674 -1.406) 4.469 1.508 ( 2.878 4.985 3.703) 6.845 1.598 ( 1.276 2.210 -3.293) 4.166 3.542 ( 5.013 8.683 4.678) 11.063 3.661 ( 4.619 8.000 -4.770) 10.396 3.978 ( 13.418 23.240 2.349) 26.938 4.039 ( 13.096 22.682 -2.540) 26.314 5.257 ( 13.094 22.680 1.180) 26.215 5.278 ( 12.997 22.511 -0.472) 25.997 6.412 ( -1.427 -2.471 0.323) 2.872 6.486 ( 1.625 2.814 -5.839) 6.682 6.631 ( -1.346 -2.331 5.471) 6.097 6.717 ( -0.510 -0.884 -1.005) 1.433 7.226 ( -0.265 -0.459 0.702) 0.879 7.240 ( 0.558 0.967 -0.564) 1.251 7.771 ( -2.576 -4.462 5.105) 7.253 7.897 ( -0.458 -0.793 -5.165) 5.246 7.913 ( -1.822 -3.155 2.602) 4.477 7.988 ( -1.337 -2.316 -3.260) 4.217 8.950 ( -5.794 -10.036 4.746) 12.522 9.057 ( -8.599 -14.894 -3.859) 17.626 ======================= Grid point 291 (49/54) ======================= q-point: ( 0.27 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.019 ( -2.092 -2.966 0.441) 3.656 1.032 ( -1.582 -3.671 -0.637) 4.047 1.367 ( 1.783 4.302 2.262) 5.178 1.411 ( 0.785 4.393 -1.348) 4.662 1.604 ( 2.497 3.624 1.501) 4.650 1.639 ( 1.394 0.964 -1.258) 2.111 3.716 ( 5.951 5.634 3.761) 9.017 3.809 ( 3.873 5.975 -3.676) 8.013 4.448 ( 9.482 27.173 1.622) 28.825 4.490 ( 8.698 25.975 -1.764) 27.450 5.679 ( 11.140 14.299 0.312) 18.129 5.684 ( 10.905 13.075 -0.031) 17.025 6.314 ( -3.784 -5.709 6.901) 9.723 6.511 ( -0.973 0.590 -8.987) 9.058 6.585 ( -3.679 1.241 2.693) 4.725 6.661 ( -6.810 1.862 -3.035) 7.684 7.191 ( 2.033 -6.186 1.548) 6.693 7.227 ( 1.608 -4.897 -1.039) 5.258 7.708 ( -4.341 4.241 2.933) 6.741 7.798 ( -10.404 8.942 -3.614) 14.186 7.925 ( -0.709 6.230 6.355) 8.927 8.057 ( 3.641 4.920 -3.704) 7.154 8.766 ( -0.339 -14.988 0.384) 14.996 8.781 ( -2.127 -19.952 -0.778) 20.080 ======================= Grid point 292 (50/54) ======================= q-point: ( 0.33 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.965 ( -0.674 -2.420 0.311) 2.532 0.974 ( -1.151 -2.336 -0.396) 2.634 1.430 ( -0.031 4.447 1.346) 4.647 1.459 ( -0.684 4.667 -1.020) 4.826 1.663 ( 1.154 0.574 0.278) 1.318 1.668 ( 1.696 0.015 -0.155) 1.703 3.858 ( 4.507 2.742 2.390) 5.791 3.917 ( 1.754 4.647 -2.341) 5.491 4.850 ( 1.325 29.516 0.418) 29.549 4.860 ( 0.524 27.705 -0.434) 27.713 5.897 ( 4.241 3.435 0.464) 5.478 5.921 ( 4.418 3.436 -1.838) 5.891 6.211 ( -2.238 -1.816 10.383) 10.776 6.446 ( -5.344 -3.460 -8.496) 10.617 6.516 ( -5.653 3.266 1.119) 6.624 6.555 ( -6.620 5.498 -2.278) 8.902 7.198 ( 4.912 -6.065 0.305) 7.810 7.236 ( 6.320 -6.384 -1.964) 9.195 7.674 ( -10.815 13.860 0.014) 17.580 7.706 ( -10.983 11.759 -1.605) 16.171 8.054 ( 8.092 5.514 8.658) 13.071 8.184 ( 5.984 2.469 -2.576) 6.967 8.549 ( -2.443 -18.204 1.769) 18.452 8.621 ( 1.136 -14.669 -3.943) 15.232 ======================= Grid point 293 (51/54) ======================= q-point: ( 0.40 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.947 ( 0.883 -1.530 0.131) 1.772 0.951 ( 0.505 -0.874 -0.184) 1.026 1.453 ( -2.190 3.794 1.131) 4.524 1.478 ( -2.433 4.215 -0.900) 4.949 1.673 ( 0.415 -0.718 0.500) 0.969 1.686 ( 0.793 -1.374 -0.517) 1.669 3.915 ( 0.212 -0.368 1.704) 1.756 3.958 ( -1.655 2.867 -1.727) 3.734 5.022 ( -13.530 23.435 0.449) 27.064 5.032 ( -14.397 24.937 -0.382) 28.797 5.947 ( 0.309 -0.536 0.056) 0.621 5.956 ( 0.645 -1.117 -0.412) 1.354 6.186 ( -0.106 0.184 7.703) 7.706 6.351 ( -0.181 0.314 -5.539) 5.551 6.485 ( -2.952 5.112 1.213) 6.026 6.529 ( -2.624 4.546 -2.440) 5.788 7.214 ( 4.007 -6.940 0.427) 8.025 7.277 ( 4.550 -7.880 -4.772) 10.275 7.627 ( -10.029 17.371 2.700) 20.239 7.675 ( -8.254 14.297 -1.530) 16.579 8.247 ( 4.045 -7.006 4.516) 9.265 8.302 ( 3.095 -5.361 -1.779) 6.441 8.322 ( 2.228 -3.858 7.829) 9.008 8.535 ( 5.321 -9.216 -7.909) 13.259 ======================= Grid point 307 (52/54) ======================= q-point: ( 0.27 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.958 ( -1.403 -2.429 0.578) 2.864 0.973 ( -1.300 -2.252 -0.645) 2.679 1.455 ( 2.326 4.029 1.674) 4.944 1.492 ( 1.948 3.374 -1.264) 4.096 1.651 ( 0.097 0.168 0.225) 0.297 1.654 ( 0.666 1.153 -0.052) 1.333 3.830 ( 2.949 5.107 3.882) 7.060 3.925 ( 2.908 5.038 -3.766) 6.930 5.003 ( 14.900 25.807 0.179) 29.800 5.012 ( 13.768 23.847 -0.540) 27.542 5.873 ( 3.014 5.220 0.304) 6.035 5.895 ( 3.116 5.397 -1.927) 6.523 6.224 ( -1.299 -2.250 10.593) 10.907 6.472 ( -1.946 -3.371 -10.163) 10.883 6.604 ( -0.178 -0.308 3.869) 3.885 6.709 ( 0.186 0.322 -4.434) 4.449 7.048 ( -3.368 -5.833 3.030) 7.386 7.116 ( -2.820 -4.884 -2.028) 5.994 7.794 ( 2.650 4.589 2.736) 5.964 7.918 ( 1.501 2.600 -6.804) 7.437 8.084 ( 4.243 7.350 7.336) 11.218 8.182 ( 2.986 5.172 -0.603) 6.002 8.421 ( -7.754 -13.430 2.242) 15.668 8.496 ( -5.808 -10.060 -3.892) 12.252 ======================= Grid point 308 (53/54) ======================= q-point: ( 0.33 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.927 ( -0.102 -1.086 0.454) 1.182 0.938 ( -0.410 -1.325 -0.481) 1.468 1.521 ( -0.265 4.188 1.143) 4.350 1.546 ( -0.170 3.681 -0.928) 3.800 1.649 ( 0.630 -1.329 0.345) 1.511 1.657 ( 0.318 -1.536 -0.279) 1.594 3.910 ( 1.348 2.046 3.777) 4.502 4.003 ( 0.499 3.282 -3.656) 4.939 5.421 ( 0.413 25.624 0.989) 25.647 5.451 ( 1.086 27.241 -1.468) 27.302 5.920 ( 1.221 -1.069 0.114) 1.626 5.937 ( 0.700 -0.248 -1.864) 2.007 6.184 ( -1.304 -2.738 7.473) 8.065 6.366 ( -4.148 -5.364 -8.543) 10.906 6.582 ( -2.353 2.884 4.436) 5.791 6.644 ( -6.121 2.800 0.474) 6.748 7.014 ( 8.264 -11.892 0.913) 14.510 7.051 ( 3.092 -8.539 -1.905) 9.279 7.876 ( -0.401 5.641 -1.237) 5.789 7.919 ( -5.564 6.711 -3.169) 9.276 8.184 ( -1.086 6.039 4.421) 7.563 8.207 ( -0.514 -1.772 -0.149) 1.851 8.281 ( 1.405 -5.371 3.421) 6.521 8.386 ( 2.850 -7.058 -3.946) 8.574 ======================= Grid point 323 (54/54) ======================= q-point: ( 0.33 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 2.00e-05 2.00e-05 Number of triplets: 73 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.920 ( -0.184 -0.319 -0.010) 0.368 0.920 ( 0.184 0.319 -0.010) 0.368 1.587 ( -0.000 -0.000 0.808) 0.808 1.606 ( -0.000 -0.000 -0.675) 0.675 1.616 ( -0.140 -0.243 0.052) 0.286 1.616 ( 0.140 0.243 0.052) 0.286 3.932 ( -0.000 -0.000 4.340) 4.340 4.039 ( -0.000 -0.000 -4.159) 4.159 5.845 ( -4.016 -6.956 0.571) 8.053 5.845 ( 4.016 6.956 0.571) 8.053 5.855 ( -6.714 -11.629 -0.040) 13.428 5.855 ( 6.714 11.629 -0.040) 13.428 6.187 ( -4.578 -7.929 -3.877) 9.943 6.187 ( 4.578 7.929 -3.877) 9.943 6.643 ( -1.436 -2.486 5.087) 5.841 6.643 ( 1.436 2.486 5.087) 5.841 6.849 ( 0.000 0.000 3.970) 3.970 6.954 ( 0.000 0.000 -4.486) 4.486 7.947 ( -1.066 -1.846 -4.331) 4.827 7.947 ( 1.066 1.846 -4.331) 4.827 8.181 ( -0.000 -0.000 2.247) 2.247 8.231 ( -0.000 -0.000 -1.805) 1.805 8.300 ( -1.666 -2.885 2.582) 4.215 8.300 ( 1.666 2.885 2.582) 4.215 =================== 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/16200 10.0 529.533 529.533 86.199 -0.000 0.000 0.000 3/16200 20.0 157.890 157.890 55.223 -0.000 0.000 0.000 3/16200 30.0 72.691 72.691 34.946 -0.000 0.000 0.000 3/16200 40.0 43.796 43.796 24.668 -0.000 0.000 0.000 3/16200 50.0 30.407 30.407 18.721 -0.000 0.000 0.000 3/16200 60.0 23.039 23.039 14.879 -0.000 0.000 0.000 3/16200 70.0 18.538 18.538 12.238 -0.000 0.000 0.000 3/16200 80.0 15.561 15.561 10.347 -0.000 0.000 0.000 3/16200 90.0 13.460 13.460 8.946 -0.000 0.000 0.000 3/16200 100.0 11.900 11.900 7.876 -0.000 0.000 0.000 3/16200 110.0 10.693 10.693 7.037 -0.000 0.000 0.000 3/16200 120.0 9.727 9.727 6.363 -0.000 0.000 0.000 3/16200 130.0 8.935 8.935 5.812 -0.000 0.000 0.000 3/16200 140.0 8.271 8.271 5.352 -0.000 0.000 0.000 3/16200 150.0 7.704 7.704 4.962 -0.000 0.000 0.000 3/16200 160.0 7.214 7.214 4.628 -0.000 0.000 0.000 3/16200 170.0 6.786 6.786 4.338 -0.000 0.000 0.000 3/16200 180.0 6.407 6.407 4.084 -0.000 0.000 0.000 3/16200 190.0 6.070 6.070 3.859 -0.000 0.000 0.000 3/16200 200.0 5.767 5.767 3.659 -0.000 0.000 0.000 3/16200 210.0 5.494 5.494 3.479 -0.000 0.000 0.000 3/16200 220.0 5.246 5.246 3.317 -0.000 0.000 0.000 3/16200 230.0 5.020 5.020 3.170 -0.000 0.000 0.000 3/16200 240.0 4.813 4.813 3.036 -0.000 0.000 0.000 3/16200 250.0 4.622 4.622 2.913 -0.000 0.000 0.000 3/16200 260.0 4.446 4.446 2.800 -0.000 0.000 0.000 3/16200 270.0 4.283 4.283 2.695 -0.000 0.000 0.000 3/16200 280.0 4.132 4.132 2.599 -0.000 0.000 0.000 3/16200 290.0 3.991 3.991 2.509 -0.000 0.000 0.000 3/16200 300.0 3.859 3.859 2.426 -0.000 0.000 0.000 3/16200 310.0 3.736 3.736 2.348 -0.000 0.000 0.000 3/16200 320.0 3.621 3.621 2.275 -0.000 0.000 0.000 3/16200 330.0 3.512 3.512 2.206 -0.000 0.000 0.000 3/16200 340.0 3.410 3.410 2.142 -0.000 0.000 0.000 3/16200 350.0 3.313 3.313 2.081 -0.000 0.000 0.000 3/16200 360.0 3.222 3.222 2.024 -0.000 0.000 0.000 3/16200 370.0 3.136 3.136 1.970 -0.000 0.000 0.000 3/16200 380.0 3.054 3.054 1.918 -0.000 0.000 0.000 3/16200 390.0 2.976 2.976 1.870 -0.000 0.000 0.000 3/16200 400.0 2.903 2.903 1.824 -0.000 0.000 0.000 3/16200 410.0 2.832 2.832 1.780 -0.000 0.000 0.000 3/16200 420.0 2.766 2.766 1.738 -0.000 0.000 0.000 3/16200 430.0 2.702 2.702 1.698 -0.000 0.000 0.000 3/16200 440.0 2.641 2.641 1.660 -0.000 0.000 0.000 3/16200 450.0 2.583 2.583 1.624 -0.000 0.000 0.000 3/16200 460.0 2.527 2.527 1.589 -0.000 0.000 0.000 3/16200 470.0 2.473 2.473 1.555 -0.000 0.000 0.000 3/16200 480.0 2.422 2.422 1.524 -0.000 0.000 0.000 3/16200 490.0 2.373 2.373 1.493 -0.000 0.000 0.000 3/16200 500.0 2.326 2.326 1.464 -0.000 0.000 0.000 3/16200 510.0 2.281 2.281 1.435 -0.000 0.000 0.000 3/16200 520.0 2.237 2.237 1.408 -0.000 0.000 0.000 3/16200 530.0 2.195 2.195 1.382 -0.000 0.000 0.000 3/16200 540.0 2.155 2.155 1.357 -0.000 0.000 0.000 3/16200 550.0 2.116 2.116 1.332 -0.000 0.000 0.000 3/16200 560.0 2.078 2.078 1.309 -0.000 0.000 0.000 3/16200 570.0 2.042 2.042 1.286 -0.000 0.000 0.000 3/16200 580.0 2.007 2.007 1.265 -0.000 0.000 0.000 3/16200 590.0 1.973 1.973 1.244 -0.000 0.000 0.000 3/16200 600.0 1.940 1.940 1.223 -0.000 0.000 0.000 3/16200 610.0 1.909 1.909 1.203 -0.000 0.000 0.000 3/16200 620.0 1.878 1.878 1.184 -0.000 0.000 0.000 3/16200 630.0 1.849 1.849 1.166 -0.000 0.000 0.000 3/16200 640.0 1.820 1.820 1.148 -0.000 0.000 0.000 3/16200 650.0 1.792 1.792 1.131 -0.000 0.000 0.000 3/16200 660.0 1.765 1.765 1.114 -0.000 0.000 0.000 3/16200 670.0 1.739 1.739 1.097 -0.000 0.000 0.000 3/16200 680.0 1.713 1.713 1.081 -0.000 0.000 0.000 3/16200 690.0 1.688 1.688 1.066 -0.000 0.000 0.000 3/16200 700.0 1.664 1.664 1.051 -0.000 0.000 0.000 3/16200 710.0 1.641 1.641 1.036 -0.000 0.000 0.000 3/16200 720.0 1.618 1.618 1.022 -0.000 0.000 0.000 3/16200 730.0 1.596 1.596 1.009 -0.000 0.000 0.000 3/16200 740.0 1.575 1.575 0.995 -0.000 0.000 0.000 3/16200 750.0 1.554 1.554 0.982 -0.000 0.000 0.000 3/16200 760.0 1.534 1.534 0.969 -0.000 0.000 0.000 3/16200 770.0 1.514 1.514 0.957 -0.000 0.000 0.000 3/16200 780.0 1.494 1.494 0.945 -0.000 0.000 0.000 3/16200 790.0 1.476 1.476 0.933 -0.000 0.000 0.000 3/16200 800.0 1.457 1.457 0.922 -0.000 0.000 0.000 3/16200 810.0 1.439 1.439 0.910 -0.000 0.000 0.000 3/16200 820.0 1.422 1.422 0.899 -0.000 0.000 0.000 3/16200 830.0 1.405 1.405 0.889 -0.000 0.000 0.000 3/16200 840.0 1.388 1.388 0.878 -0.000 0.000 0.000 3/16200 850.0 1.372 1.372 0.868 -0.000 0.000 0.000 3/16200 860.0 1.356 1.356 0.858 -0.000 0.000 0.000 3/16200 870.0 1.340 1.340 0.849 -0.000 0.000 0.000 3/16200 880.0 1.325 1.325 0.839 -0.000 0.000 0.000 3/16200 890.0 1.310 1.310 0.830 -0.000 0.000 0.000 3/16200 900.0 1.296 1.296 0.821 -0.000 0.000 0.000 3/16200 910.0 1.282 1.282 0.812 -0.000 0.000 0.000 3/16200 920.0 1.268 1.268 0.803 -0.000 0.000 0.000 3/16200 930.0 1.254 1.254 0.795 -0.000 0.000 0.000 3/16200 940.0 1.241 1.241 0.786 -0.000 0.000 0.000 3/16200 950.0 1.228 1.228 0.778 -0.000 0.000 0.000 3/16200 960.0 1.215 1.215 0.770 -0.000 0.000 0.000 3/16200 970.0 1.203 1.203 0.762 -0.000 0.000 0.000 3/16200 980.0 1.190 1.190 0.755 -0.000 0.000 0.000 3/16200 990.0 1.178 1.178 0.747 -0.000 0.000 0.000 3/16200 1000.0 1.167 1.167 0.740 -0.000 0.000 0.000 3/16200 Thermal conductivity related properties were written into "kappa-m15153.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 00:15:15]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|