----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 19:10:33]------------------------- 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.563962938187290 0.000000000000000 0.000000000000000 b -1.781981469093645 3.086482442616422 0.000000000000000 c 0.000000000000000 0.000000000000000 15.916651770000000 Atomic positions (fractional): *1 S 0.33333333333333 0.66666666666667 0.10531804128313 32.065 2 S 0.66666666666667 0.33333333333333 0.60531804128313 32.065 3 S 0.66666666666667 0.33333333333333 0.89468195871687 32.065 4 S 0.33333333333333 0.66666666666667 0.39468195871687 32.065 *5 Ga 0.33333333333333 0.66666666666667 0.67359408206119 69.723 6 Ga 0.66666666666667 0.33333333333333 0.17359408206119 69.723 7 Ga 0.66666666666667 0.33333333333333 0.32640591793881 69.723 8 Ga 0.33333333333333 0.66666666666667 0.82640591793881 69.723 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.563962938187290 0.000000000000000 0.000000000000000 b -1.781981469093645 3.086482442616422 0.000000000000000 c 0.000000000000000 0.000000000000000 15.916651770000000 Atomic positions (fractional): *1 S 0.33333333333333 0.66666666666667 0.10531804128313 32.065 > 1 2 S 0.66666666666667 0.33333333333333 0.60531804128313 32.065 > 2 3 S 0.66666666666667 0.33333333333333 0.89468195871687 32.065 > 3 4 S 0.33333333333333 0.66666666666667 0.39468195871687 32.065 > 4 *5 Ga 0.33333333333333 0.66666666666667 0.67359408206119 69.723 > 5 6 Ga 0.66666666666667 0.33333333333333 0.17359408206119 69.723 > 6 7 Ga 0.66666666666667 0.33333333333333 0.32640591793881 69.723 > 7 8 Ga 0.33333333333333 0.66666666666667 0.82640591793881 69.723 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 14.255851752749161 0.000000000000000 0.000000000000000 b -7.127925876374580 12.345929770465688 0.000000000000000 c 0.000000000000000 0.000000000000000 15.916651770000000 Atomic positions (fractional): *1 S 0.08333333333333 0.16666666666667 0.10531804128313 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.10531804128313 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.10531804128313 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.10531804128313 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.10531804128313 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.10531804128313 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.10531804128313 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.10531804128313 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.10531804128313 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.10531804128313 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.10531804128313 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.10531804128313 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.10531804128313 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.10531804128313 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.10531804128313 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.10531804128313 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.60531804128313 32.065 > 2 18 S 0.41666666666667 0.08333333333333 0.60531804128313 32.065 > 2 19 S 0.66666666666667 0.08333333333333 0.60531804128313 32.065 > 2 20 S 0.91666666666667 0.08333333333333 0.60531804128313 32.065 > 2 21 S 0.16666666666667 0.33333333333333 0.60531804128313 32.065 > 2 22 S 0.41666666666667 0.33333333333333 0.60531804128313 32.065 > 2 23 S 0.66666666666667 0.33333333333333 0.60531804128313 32.065 > 2 24 S 0.91666666666667 0.33333333333333 0.60531804128313 32.065 > 2 25 S 0.16666666666667 0.58333333333333 0.60531804128313 32.065 > 2 26 S 0.41666666666667 0.58333333333333 0.60531804128313 32.065 > 2 27 S 0.66666666666667 0.58333333333333 0.60531804128313 32.065 > 2 28 S 0.91666666666667 0.58333333333333 0.60531804128313 32.065 > 2 29 S 0.16666666666667 0.83333333333333 0.60531804128313 32.065 > 2 30 S 0.41666666666667 0.83333333333333 0.60531804128313 32.065 > 2 31 S 0.66666666666667 0.83333333333333 0.60531804128313 32.065 > 2 32 S 0.91666666666667 0.83333333333333 0.60531804128313 32.065 > 2 33 S 0.16666666666667 0.08333333333333 0.89468195871687 32.065 > 3 34 S 0.41666666666667 0.08333333333333 0.89468195871687 32.065 > 3 35 S 0.66666666666667 0.08333333333333 0.89468195871687 32.065 > 3 36 S 0.91666666666667 0.08333333333333 0.89468195871687 32.065 > 3 37 S 0.16666666666667 0.33333333333333 0.89468195871687 32.065 > 3 38 S 0.41666666666667 0.33333333333333 0.89468195871687 32.065 > 3 39 S 0.66666666666667 0.33333333333333 0.89468195871687 32.065 > 3 40 S 0.91666666666667 0.33333333333333 0.89468195871687 32.065 > 3 41 S 0.16666666666667 0.58333333333333 0.89468195871687 32.065 > 3 42 S 0.41666666666667 0.58333333333333 0.89468195871687 32.065 > 3 43 S 0.66666666666667 0.58333333333333 0.89468195871687 32.065 > 3 44 S 0.91666666666667 0.58333333333333 0.89468195871687 32.065 > 3 45 S 0.16666666666667 0.83333333333333 0.89468195871687 32.065 > 3 46 S 0.41666666666667 0.83333333333333 0.89468195871687 32.065 > 3 47 S 0.66666666666667 0.83333333333333 0.89468195871687 32.065 > 3 48 S 0.91666666666667 0.83333333333333 0.89468195871687 32.065 > 3 49 S 0.08333333333333 0.16666666666667 0.39468195871687 32.065 > 4 50 S 0.33333333333333 0.16666666666667 0.39468195871687 32.065 > 4 51 S 0.58333333333333 0.16666666666667 0.39468195871687 32.065 > 4 52 S 0.83333333333333 0.16666666666667 0.39468195871687 32.065 > 4 53 S 0.08333333333333 0.41666666666667 0.39468195871687 32.065 > 4 54 S 0.33333333333333 0.41666666666667 0.39468195871687 32.065 > 4 55 S 0.58333333333333 0.41666666666667 0.39468195871687 32.065 > 4 56 S 0.83333333333333 0.41666666666667 0.39468195871687 32.065 > 4 57 S 0.08333333333333 0.66666666666667 0.39468195871687 32.065 > 4 58 S 0.33333333333333 0.66666666666667 0.39468195871687 32.065 > 4 59 S 0.58333333333333 0.66666666666667 0.39468195871687 32.065 > 4 60 S 0.83333333333333 0.66666666666667 0.39468195871687 32.065 > 4 61 S 0.08333333333333 0.91666666666667 0.39468195871687 32.065 > 4 62 S 0.33333333333333 0.91666666666667 0.39468195871687 32.065 > 4 63 S 0.58333333333333 0.91666666666667 0.39468195871687 32.065 > 4 64 S 0.83333333333333 0.91666666666667 0.39468195871687 32.065 > 4 *65 Ga 0.08333333333333 0.16666666666667 0.67359408206119 69.723 > 5 66 Ga 0.33333333333333 0.16666666666667 0.67359408206119 69.723 > 5 67 Ga 0.58333333333333 0.16666666666667 0.67359408206119 69.723 > 5 68 Ga 0.83333333333333 0.16666666666667 0.67359408206119 69.723 > 5 69 Ga 0.08333333333333 0.41666666666667 0.67359408206119 69.723 > 5 70 Ga 0.33333333333333 0.41666666666667 0.67359408206119 69.723 > 5 71 Ga 0.58333333333333 0.41666666666667 0.67359408206119 69.723 > 5 72 Ga 0.83333333333333 0.41666666666667 0.67359408206119 69.723 > 5 73 Ga 0.08333333333333 0.66666666666667 0.67359408206119 69.723 > 5 74 Ga 0.33333333333333 0.66666666666667 0.67359408206119 69.723 > 5 75 Ga 0.58333333333333 0.66666666666667 0.67359408206119 69.723 > 5 76 Ga 0.83333333333333 0.66666666666667 0.67359408206119 69.723 > 5 77 Ga 0.08333333333333 0.91666666666667 0.67359408206119 69.723 > 5 78 Ga 0.33333333333333 0.91666666666667 0.67359408206119 69.723 > 5 79 Ga 0.58333333333333 0.91666666666667 0.67359408206119 69.723 > 5 80 Ga 0.83333333333333 0.91666666666667 0.67359408206119 69.723 > 5 81 Ga 0.16666666666667 0.08333333333333 0.17359408206119 69.723 > 6 82 Ga 0.41666666666667 0.08333333333333 0.17359408206119 69.723 > 6 83 Ga 0.66666666666667 0.08333333333333 0.17359408206119 69.723 > 6 84 Ga 0.91666666666667 0.08333333333333 0.17359408206119 69.723 > 6 85 Ga 0.16666666666667 0.33333333333333 0.17359408206119 69.723 > 6 86 Ga 0.41666666666667 0.33333333333333 0.17359408206119 69.723 > 6 87 Ga 0.66666666666667 0.33333333333333 0.17359408206119 69.723 > 6 88 Ga 0.91666666666667 0.33333333333333 0.17359408206119 69.723 > 6 89 Ga 0.16666666666667 0.58333333333333 0.17359408206119 69.723 > 6 90 Ga 0.41666666666667 0.58333333333333 0.17359408206119 69.723 > 6 91 Ga 0.66666666666667 0.58333333333333 0.17359408206119 69.723 > 6 92 Ga 0.91666666666667 0.58333333333333 0.17359408206119 69.723 > 6 93 Ga 0.16666666666667 0.83333333333333 0.17359408206119 69.723 > 6 94 Ga 0.41666666666667 0.83333333333333 0.17359408206119 69.723 > 6 95 Ga 0.66666666666667 0.83333333333333 0.17359408206119 69.723 > 6 96 Ga 0.91666666666667 0.83333333333333 0.17359408206119 69.723 > 6 97 Ga 0.16666666666667 0.08333333333333 0.32640591793881 69.723 > 7 98 Ga 0.41666666666667 0.08333333333333 0.32640591793881 69.723 > 7 99 Ga 0.66666666666667 0.08333333333333 0.32640591793881 69.723 > 7 100 Ga 0.91666666666667 0.08333333333333 0.32640591793881 69.723 > 7 101 Ga 0.16666666666667 0.33333333333333 0.32640591793881 69.723 > 7 102 Ga 0.41666666666667 0.33333333333333 0.32640591793881 69.723 > 7 103 Ga 0.66666666666667 0.33333333333333 0.32640591793881 69.723 > 7 104 Ga 0.91666666666667 0.33333333333333 0.32640591793881 69.723 > 7 105 Ga 0.16666666666667 0.58333333333333 0.32640591793881 69.723 > 7 106 Ga 0.41666666666667 0.58333333333333 0.32640591793881 69.723 > 7 107 Ga 0.66666666666667 0.58333333333333 0.32640591793881 69.723 > 7 108 Ga 0.91666666666667 0.58333333333333 0.32640591793881 69.723 > 7 109 Ga 0.16666666666667 0.83333333333333 0.32640591793881 69.723 > 7 110 Ga 0.41666666666667 0.83333333333333 0.32640591793881 69.723 > 7 111 Ga 0.66666666666667 0.83333333333333 0.32640591793881 69.723 > 7 112 Ga 0.91666666666667 0.83333333333333 0.32640591793881 69.723 > 7 113 Ga 0.08333333333333 0.16666666666667 0.82640591793881 69.723 > 8 114 Ga 0.33333333333333 0.16666666666667 0.82640591793881 69.723 > 8 115 Ga 0.58333333333333 0.16666666666667 0.82640591793881 69.723 > 8 116 Ga 0.83333333333333 0.16666666666667 0.82640591793881 69.723 > 8 117 Ga 0.08333333333333 0.41666666666667 0.82640591793881 69.723 > 8 118 Ga 0.33333333333333 0.41666666666667 0.82640591793881 69.723 > 8 119 Ga 0.58333333333333 0.41666666666667 0.82640591793881 69.723 > 8 120 Ga 0.83333333333333 0.41666666666667 0.82640591793881 69.723 > 8 121 Ga 0.08333333333333 0.66666666666667 0.82640591793881 69.723 > 8 122 Ga 0.33333333333333 0.66666666666667 0.82640591793881 69.723 > 8 123 Ga 0.58333333333333 0.66666666666667 0.82640591793881 69.723 > 8 124 Ga 0.83333333333333 0.66666666666667 0.82640591793881 69.723 > 8 125 Ga 0.08333333333333 0.91666666666667 0.82640591793881 69.723 > 8 126 Ga 0.33333333333333 0.91666666666667 0.82640591793881 69.723 > 8 127 Ga 0.58333333333333 0.91666666666667 0.82640591793881 69.723 > 8 128 Ga 0.83333333333333 0.91666666666667 0.82640591793881 69.723 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.3160896 0.0000000 0.0000000 -0.0000000 7.3160896 0.0000000 0.0000000 0.0000000 5.3074479 -------------------------- Born effective charges -------------------------- 1 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 2 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 3 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 4 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 5 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 6 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 7 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 8 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 ---------------------------------------------------------------------------- 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: 252 Number of blocks in projector: 192 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: 116 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 56 Use standard eigh solver. Tree of FC basis block matrices: - (252, 248), data: False |-- (56, 54), data: True |-- (80, 80), data: True |-- (116, 114), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.019 Solver_block: 80 / 80 - Time: 0.275 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.318 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) 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: 252 Number of blocks in projector: 192 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: 116 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 56 Use standard eigh solver. Tree of FC basis block matrices: - (252, 248), data: False |-- (56, 54), data: True |-- (80, 80), data: True |-- (116, 114), data: True Max drift after symmetrization by symfc projector: 0.00000000 (zz) 0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 19:10:39]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 19:10:40]------------------------- 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.563962938187290 0.000000000000000 0.000000000000000 b -1.781981469093645 3.086482442616422 0.000000000000000 c 0.000000000000000 0.000000000000000 15.916651770000000 Atomic positions (fractional): 1 S 0.33333333333333 0.66666666666667 0.10531804128313 32.065 2 S 0.66666666666667 0.33333333333333 0.60531804128313 32.065 3 S 0.66666666666667 0.33333333333333 0.89468195871687 32.065 4 S 0.33333333333333 0.66666666666667 0.39468195871687 32.065 5 Ga 0.33333333333333 0.66666666666667 0.67359408206119 69.723 6 Ga 0.66666666666667 0.33333333333333 0.17359408206119 69.723 7 Ga 0.66666666666667 0.33333333333333 0.32640591793881 69.723 8 Ga 0.33333333333333 0.66666666666667 0.82640591793881 69.723 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.255851752749161 0.000000000000000 0.000000000000000 b -7.127925876374580 12.345929770465688 0.000000000000000 c 0.000000000000000 0.000000000000000 15.916651770000000 Atomic positions (fractional): 1 S 0.08333333333333 0.16666666666667 0.10531804128313 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.10531804128313 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.10531804128313 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.10531804128313 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.10531804128313 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.10531804128313 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.10531804128313 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.10531804128313 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.10531804128313 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.10531804128313 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.10531804128313 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.10531804128313 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.10531804128313 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.10531804128313 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.10531804128313 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.10531804128313 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.60531804128313 32.065 > 17 18 S 0.41666666666667 0.08333333333333 0.60531804128313 32.065 > 17 19 S 0.66666666666667 0.08333333333333 0.60531804128313 32.065 > 17 20 S 0.91666666666667 0.08333333333333 0.60531804128313 32.065 > 17 21 S 0.16666666666667 0.33333333333333 0.60531804128313 32.065 > 17 22 S 0.41666666666667 0.33333333333333 0.60531804128313 32.065 > 17 23 S 0.66666666666667 0.33333333333333 0.60531804128313 32.065 > 17 24 S 0.91666666666667 0.33333333333333 0.60531804128313 32.065 > 17 25 S 0.16666666666667 0.58333333333333 0.60531804128313 32.065 > 17 26 S 0.41666666666667 0.58333333333333 0.60531804128313 32.065 > 17 27 S 0.66666666666667 0.58333333333333 0.60531804128313 32.065 > 17 28 S 0.91666666666667 0.58333333333333 0.60531804128313 32.065 > 17 29 S 0.16666666666667 0.83333333333333 0.60531804128313 32.065 > 17 30 S 0.41666666666667 0.83333333333333 0.60531804128313 32.065 > 17 31 S 0.66666666666667 0.83333333333333 0.60531804128313 32.065 > 17 32 S 0.91666666666667 0.83333333333333 0.60531804128313 32.065 > 17 33 S 0.16666666666667 0.08333333333333 0.89468195871687 32.065 > 33 34 S 0.41666666666667 0.08333333333333 0.89468195871687 32.065 > 33 35 S 0.66666666666667 0.08333333333333 0.89468195871687 32.065 > 33 36 S 0.91666666666667 0.08333333333333 0.89468195871687 32.065 > 33 37 S 0.16666666666667 0.33333333333333 0.89468195871687 32.065 > 33 38 S 0.41666666666667 0.33333333333333 0.89468195871687 32.065 > 33 39 S 0.66666666666667 0.33333333333333 0.89468195871687 32.065 > 33 40 S 0.91666666666667 0.33333333333333 0.89468195871687 32.065 > 33 41 S 0.16666666666667 0.58333333333333 0.89468195871687 32.065 > 33 42 S 0.41666666666667 0.58333333333333 0.89468195871687 32.065 > 33 43 S 0.66666666666667 0.58333333333333 0.89468195871687 32.065 > 33 44 S 0.91666666666667 0.58333333333333 0.89468195871687 32.065 > 33 45 S 0.16666666666667 0.83333333333333 0.89468195871687 32.065 > 33 46 S 0.41666666666667 0.83333333333333 0.89468195871687 32.065 > 33 47 S 0.66666666666667 0.83333333333333 0.89468195871687 32.065 > 33 48 S 0.91666666666667 0.83333333333333 0.89468195871687 32.065 > 33 49 S 0.08333333333333 0.16666666666667 0.39468195871687 32.065 > 49 50 S 0.33333333333333 0.16666666666667 0.39468195871687 32.065 > 49 51 S 0.58333333333333 0.16666666666667 0.39468195871687 32.065 > 49 52 S 0.83333333333333 0.16666666666667 0.39468195871687 32.065 > 49 53 S 0.08333333333333 0.41666666666667 0.39468195871687 32.065 > 49 54 S 0.33333333333333 0.41666666666667 0.39468195871687 32.065 > 49 55 S 0.58333333333333 0.41666666666667 0.39468195871687 32.065 > 49 56 S 0.83333333333333 0.41666666666667 0.39468195871687 32.065 > 49 57 S 0.08333333333333 0.66666666666667 0.39468195871687 32.065 > 49 58 S 0.33333333333333 0.66666666666667 0.39468195871687 32.065 > 49 59 S 0.58333333333333 0.66666666666667 0.39468195871687 32.065 > 49 60 S 0.83333333333333 0.66666666666667 0.39468195871687 32.065 > 49 61 S 0.08333333333333 0.91666666666667 0.39468195871687 32.065 > 49 62 S 0.33333333333333 0.91666666666667 0.39468195871687 32.065 > 49 63 S 0.58333333333333 0.91666666666667 0.39468195871687 32.065 > 49 64 S 0.83333333333333 0.91666666666667 0.39468195871687 32.065 > 49 65 Ga 0.08333333333333 0.16666666666667 0.67359408206119 69.723 > 65 66 Ga 0.33333333333333 0.16666666666667 0.67359408206119 69.723 > 65 67 Ga 0.58333333333333 0.16666666666667 0.67359408206119 69.723 > 65 68 Ga 0.83333333333333 0.16666666666667 0.67359408206119 69.723 > 65 69 Ga 0.08333333333333 0.41666666666667 0.67359408206119 69.723 > 65 70 Ga 0.33333333333333 0.41666666666667 0.67359408206119 69.723 > 65 71 Ga 0.58333333333333 0.41666666666667 0.67359408206119 69.723 > 65 72 Ga 0.83333333333333 0.41666666666667 0.67359408206119 69.723 > 65 73 Ga 0.08333333333333 0.66666666666667 0.67359408206119 69.723 > 65 74 Ga 0.33333333333333 0.66666666666667 0.67359408206119 69.723 > 65 75 Ga 0.58333333333333 0.66666666666667 0.67359408206119 69.723 > 65 76 Ga 0.83333333333333 0.66666666666667 0.67359408206119 69.723 > 65 77 Ga 0.08333333333333 0.91666666666667 0.67359408206119 69.723 > 65 78 Ga 0.33333333333333 0.91666666666667 0.67359408206119 69.723 > 65 79 Ga 0.58333333333333 0.91666666666667 0.67359408206119 69.723 > 65 80 Ga 0.83333333333333 0.91666666666667 0.67359408206119 69.723 > 65 81 Ga 0.16666666666667 0.08333333333333 0.17359408206119 69.723 > 81 82 Ga 0.41666666666667 0.08333333333333 0.17359408206119 69.723 > 81 83 Ga 0.66666666666667 0.08333333333333 0.17359408206119 69.723 > 81 84 Ga 0.91666666666667 0.08333333333333 0.17359408206119 69.723 > 81 85 Ga 0.16666666666667 0.33333333333333 0.17359408206119 69.723 > 81 86 Ga 0.41666666666667 0.33333333333333 0.17359408206119 69.723 > 81 87 Ga 0.66666666666667 0.33333333333333 0.17359408206119 69.723 > 81 88 Ga 0.91666666666667 0.33333333333333 0.17359408206119 69.723 > 81 89 Ga 0.16666666666667 0.58333333333333 0.17359408206119 69.723 > 81 90 Ga 0.41666666666667 0.58333333333333 0.17359408206119 69.723 > 81 91 Ga 0.66666666666667 0.58333333333333 0.17359408206119 69.723 > 81 92 Ga 0.91666666666667 0.58333333333333 0.17359408206119 69.723 > 81 93 Ga 0.16666666666667 0.83333333333333 0.17359408206119 69.723 > 81 94 Ga 0.41666666666667 0.83333333333333 0.17359408206119 69.723 > 81 95 Ga 0.66666666666667 0.83333333333333 0.17359408206119 69.723 > 81 96 Ga 0.91666666666667 0.83333333333333 0.17359408206119 69.723 > 81 97 Ga 0.16666666666667 0.08333333333333 0.32640591793881 69.723 > 97 98 Ga 0.41666666666667 0.08333333333333 0.32640591793881 69.723 > 97 99 Ga 0.66666666666667 0.08333333333333 0.32640591793881 69.723 > 97 100 Ga 0.91666666666667 0.08333333333333 0.32640591793881 69.723 > 97 101 Ga 0.16666666666667 0.33333333333333 0.32640591793881 69.723 > 97 102 Ga 0.41666666666667 0.33333333333333 0.32640591793881 69.723 > 97 103 Ga 0.66666666666667 0.33333333333333 0.32640591793881 69.723 > 97 104 Ga 0.91666666666667 0.33333333333333 0.32640591793881 69.723 > 97 105 Ga 0.16666666666667 0.58333333333333 0.32640591793881 69.723 > 97 106 Ga 0.41666666666667 0.58333333333333 0.32640591793881 69.723 > 97 107 Ga 0.66666666666667 0.58333333333333 0.32640591793881 69.723 > 97 108 Ga 0.91666666666667 0.58333333333333 0.32640591793881 69.723 > 97 109 Ga 0.16666666666667 0.83333333333333 0.32640591793881 69.723 > 97 110 Ga 0.41666666666667 0.83333333333333 0.32640591793881 69.723 > 97 111 Ga 0.66666666666667 0.83333333333333 0.32640591793881 69.723 > 97 112 Ga 0.91666666666667 0.83333333333333 0.32640591793881 69.723 > 97 113 Ga 0.08333333333333 0.16666666666667 0.82640591793881 69.723 > 113 114 Ga 0.33333333333333 0.16666666666667 0.82640591793881 69.723 > 113 115 Ga 0.58333333333333 0.16666666666667 0.82640591793881 69.723 > 113 116 Ga 0.83333333333333 0.16666666666667 0.82640591793881 69.723 > 113 117 Ga 0.08333333333333 0.41666666666667 0.82640591793881 69.723 > 113 118 Ga 0.33333333333333 0.41666666666667 0.82640591793881 69.723 > 113 119 Ga 0.58333333333333 0.41666666666667 0.82640591793881 69.723 > 113 120 Ga 0.83333333333333 0.41666666666667 0.82640591793881 69.723 > 113 121 Ga 0.08333333333333 0.66666666666667 0.82640591793881 69.723 > 113 122 Ga 0.33333333333333 0.66666666666667 0.82640591793881 69.723 > 113 123 Ga 0.58333333333333 0.66666666666667 0.82640591793881 69.723 > 113 124 Ga 0.83333333333333 0.66666666666667 0.82640591793881 69.723 > 113 125 Ga 0.08333333333333 0.91666666666667 0.82640591793881 69.723 > 113 126 Ga 0.33333333333333 0.91666666666667 0.82640591793881 69.723 > 113 127 Ga 0.58333333333333 0.91666666666667 0.82640591793881 69.723 > 113 128 Ga 0.83333333333333 0.91666666666667 0.82640591793881 69.723 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.3160896 0.0000000 0.0000000 -0.0000000 7.3160896 0.0000000 0.0000000 0.0000000 5.3074479 -------------------------- Born effective charges -------------------------- 1 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 2 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 3 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 4 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 5 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 6 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 7 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 8 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 ---------------------------------------------------------------------------- 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] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000002 (zzz) 0.00000002 (zzz) 0.00000002 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 19:10:46]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 19:10:47]------------------------- 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.563962938187290 0.000000000000000 0.000000000000000 b -1.781981469093645 3.086482442616422 0.000000000000000 c 0.000000000000000 0.000000000000000 15.916651770000000 Atomic positions (fractional): 1 S 0.33333333333333 0.66666666666667 0.10531804128313 32.065 2 S 0.66666666666667 0.33333333333333 0.60531804128313 32.065 3 S 0.66666666666667 0.33333333333333 0.89468195871687 32.065 4 S 0.33333333333333 0.66666666666667 0.39468195871687 32.065 5 Ga 0.33333333333333 0.66666666666667 0.67359408206119 69.723 6 Ga 0.66666666666667 0.33333333333333 0.17359408206119 69.723 7 Ga 0.66666666666667 0.33333333333333 0.32640591793881 69.723 8 Ga 0.33333333333333 0.66666666666667 0.82640591793881 69.723 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.255851752749161 0.000000000000000 0.000000000000000 b -7.127925876374580 12.345929770465688 0.000000000000000 c 0.000000000000000 0.000000000000000 15.916651770000000 Atomic positions (fractional): 1 S 0.08333333333333 0.16666666666667 0.10531804128313 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.10531804128313 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.10531804128313 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.10531804128313 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.10531804128313 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.10531804128313 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.10531804128313 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.10531804128313 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.10531804128313 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.10531804128313 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.10531804128313 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.10531804128313 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.10531804128313 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.10531804128313 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.10531804128313 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.10531804128313 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.60531804128313 32.065 > 17 18 S 0.41666666666667 0.08333333333333 0.60531804128313 32.065 > 17 19 S 0.66666666666667 0.08333333333333 0.60531804128313 32.065 > 17 20 S 0.91666666666667 0.08333333333333 0.60531804128313 32.065 > 17 21 S 0.16666666666667 0.33333333333333 0.60531804128313 32.065 > 17 22 S 0.41666666666667 0.33333333333333 0.60531804128313 32.065 > 17 23 S 0.66666666666667 0.33333333333333 0.60531804128313 32.065 > 17 24 S 0.91666666666667 0.33333333333333 0.60531804128313 32.065 > 17 25 S 0.16666666666667 0.58333333333333 0.60531804128313 32.065 > 17 26 S 0.41666666666667 0.58333333333333 0.60531804128313 32.065 > 17 27 S 0.66666666666667 0.58333333333333 0.60531804128313 32.065 > 17 28 S 0.91666666666667 0.58333333333333 0.60531804128313 32.065 > 17 29 S 0.16666666666667 0.83333333333333 0.60531804128313 32.065 > 17 30 S 0.41666666666667 0.83333333333333 0.60531804128313 32.065 > 17 31 S 0.66666666666667 0.83333333333333 0.60531804128313 32.065 > 17 32 S 0.91666666666667 0.83333333333333 0.60531804128313 32.065 > 17 33 S 0.16666666666667 0.08333333333333 0.89468195871687 32.065 > 33 34 S 0.41666666666667 0.08333333333333 0.89468195871687 32.065 > 33 35 S 0.66666666666667 0.08333333333333 0.89468195871687 32.065 > 33 36 S 0.91666666666667 0.08333333333333 0.89468195871687 32.065 > 33 37 S 0.16666666666667 0.33333333333333 0.89468195871687 32.065 > 33 38 S 0.41666666666667 0.33333333333333 0.89468195871687 32.065 > 33 39 S 0.66666666666667 0.33333333333333 0.89468195871687 32.065 > 33 40 S 0.91666666666667 0.33333333333333 0.89468195871687 32.065 > 33 41 S 0.16666666666667 0.58333333333333 0.89468195871687 32.065 > 33 42 S 0.41666666666667 0.58333333333333 0.89468195871687 32.065 > 33 43 S 0.66666666666667 0.58333333333333 0.89468195871687 32.065 > 33 44 S 0.91666666666667 0.58333333333333 0.89468195871687 32.065 > 33 45 S 0.16666666666667 0.83333333333333 0.89468195871687 32.065 > 33 46 S 0.41666666666667 0.83333333333333 0.89468195871687 32.065 > 33 47 S 0.66666666666667 0.83333333333333 0.89468195871687 32.065 > 33 48 S 0.91666666666667 0.83333333333333 0.89468195871687 32.065 > 33 49 S 0.08333333333333 0.16666666666667 0.39468195871687 32.065 > 49 50 S 0.33333333333333 0.16666666666667 0.39468195871687 32.065 > 49 51 S 0.58333333333333 0.16666666666667 0.39468195871687 32.065 > 49 52 S 0.83333333333333 0.16666666666667 0.39468195871687 32.065 > 49 53 S 0.08333333333333 0.41666666666667 0.39468195871687 32.065 > 49 54 S 0.33333333333333 0.41666666666667 0.39468195871687 32.065 > 49 55 S 0.58333333333333 0.41666666666667 0.39468195871687 32.065 > 49 56 S 0.83333333333333 0.41666666666667 0.39468195871687 32.065 > 49 57 S 0.08333333333333 0.66666666666667 0.39468195871687 32.065 > 49 58 S 0.33333333333333 0.66666666666667 0.39468195871687 32.065 > 49 59 S 0.58333333333333 0.66666666666667 0.39468195871687 32.065 > 49 60 S 0.83333333333333 0.66666666666667 0.39468195871687 32.065 > 49 61 S 0.08333333333333 0.91666666666667 0.39468195871687 32.065 > 49 62 S 0.33333333333333 0.91666666666667 0.39468195871687 32.065 > 49 63 S 0.58333333333333 0.91666666666667 0.39468195871687 32.065 > 49 64 S 0.83333333333333 0.91666666666667 0.39468195871687 32.065 > 49 65 Ga 0.08333333333333 0.16666666666667 0.67359408206119 69.723 > 65 66 Ga 0.33333333333333 0.16666666666667 0.67359408206119 69.723 > 65 67 Ga 0.58333333333333 0.16666666666667 0.67359408206119 69.723 > 65 68 Ga 0.83333333333333 0.16666666666667 0.67359408206119 69.723 > 65 69 Ga 0.08333333333333 0.41666666666667 0.67359408206119 69.723 > 65 70 Ga 0.33333333333333 0.41666666666667 0.67359408206119 69.723 > 65 71 Ga 0.58333333333333 0.41666666666667 0.67359408206119 69.723 > 65 72 Ga 0.83333333333333 0.41666666666667 0.67359408206119 69.723 > 65 73 Ga 0.08333333333333 0.66666666666667 0.67359408206119 69.723 > 65 74 Ga 0.33333333333333 0.66666666666667 0.67359408206119 69.723 > 65 75 Ga 0.58333333333333 0.66666666666667 0.67359408206119 69.723 > 65 76 Ga 0.83333333333333 0.66666666666667 0.67359408206119 69.723 > 65 77 Ga 0.08333333333333 0.91666666666667 0.67359408206119 69.723 > 65 78 Ga 0.33333333333333 0.91666666666667 0.67359408206119 69.723 > 65 79 Ga 0.58333333333333 0.91666666666667 0.67359408206119 69.723 > 65 80 Ga 0.83333333333333 0.91666666666667 0.67359408206119 69.723 > 65 81 Ga 0.16666666666667 0.08333333333333 0.17359408206119 69.723 > 81 82 Ga 0.41666666666667 0.08333333333333 0.17359408206119 69.723 > 81 83 Ga 0.66666666666667 0.08333333333333 0.17359408206119 69.723 > 81 84 Ga 0.91666666666667 0.08333333333333 0.17359408206119 69.723 > 81 85 Ga 0.16666666666667 0.33333333333333 0.17359408206119 69.723 > 81 86 Ga 0.41666666666667 0.33333333333333 0.17359408206119 69.723 > 81 87 Ga 0.66666666666667 0.33333333333333 0.17359408206119 69.723 > 81 88 Ga 0.91666666666667 0.33333333333333 0.17359408206119 69.723 > 81 89 Ga 0.16666666666667 0.58333333333333 0.17359408206119 69.723 > 81 90 Ga 0.41666666666667 0.58333333333333 0.17359408206119 69.723 > 81 91 Ga 0.66666666666667 0.58333333333333 0.17359408206119 69.723 > 81 92 Ga 0.91666666666667 0.58333333333333 0.17359408206119 69.723 > 81 93 Ga 0.16666666666667 0.83333333333333 0.17359408206119 69.723 > 81 94 Ga 0.41666666666667 0.83333333333333 0.17359408206119 69.723 > 81 95 Ga 0.66666666666667 0.83333333333333 0.17359408206119 69.723 > 81 96 Ga 0.91666666666667 0.83333333333333 0.17359408206119 69.723 > 81 97 Ga 0.16666666666667 0.08333333333333 0.32640591793881 69.723 > 97 98 Ga 0.41666666666667 0.08333333333333 0.32640591793881 69.723 > 97 99 Ga 0.66666666666667 0.08333333333333 0.32640591793881 69.723 > 97 100 Ga 0.91666666666667 0.08333333333333 0.32640591793881 69.723 > 97 101 Ga 0.16666666666667 0.33333333333333 0.32640591793881 69.723 > 97 102 Ga 0.41666666666667 0.33333333333333 0.32640591793881 69.723 > 97 103 Ga 0.66666666666667 0.33333333333333 0.32640591793881 69.723 > 97 104 Ga 0.91666666666667 0.33333333333333 0.32640591793881 69.723 > 97 105 Ga 0.16666666666667 0.58333333333333 0.32640591793881 69.723 > 97 106 Ga 0.41666666666667 0.58333333333333 0.32640591793881 69.723 > 97 107 Ga 0.66666666666667 0.58333333333333 0.32640591793881 69.723 > 97 108 Ga 0.91666666666667 0.58333333333333 0.32640591793881 69.723 > 97 109 Ga 0.16666666666667 0.83333333333333 0.32640591793881 69.723 > 97 110 Ga 0.41666666666667 0.83333333333333 0.32640591793881 69.723 > 97 111 Ga 0.66666666666667 0.83333333333333 0.32640591793881 69.723 > 97 112 Ga 0.91666666666667 0.83333333333333 0.32640591793881 69.723 > 97 113 Ga 0.08333333333333 0.16666666666667 0.82640591793881 69.723 > 113 114 Ga 0.33333333333333 0.16666666666667 0.82640591793881 69.723 > 113 115 Ga 0.58333333333333 0.16666666666667 0.82640591793881 69.723 > 113 116 Ga 0.83333333333333 0.16666666666667 0.82640591793881 69.723 > 113 117 Ga 0.08333333333333 0.41666666666667 0.82640591793881 69.723 > 113 118 Ga 0.33333333333333 0.41666666666667 0.82640591793881 69.723 > 113 119 Ga 0.58333333333333 0.41666666666667 0.82640591793881 69.723 > 113 120 Ga 0.83333333333333 0.41666666666667 0.82640591793881 69.723 > 113 121 Ga 0.08333333333333 0.66666666666667 0.82640591793881 69.723 > 113 122 Ga 0.33333333333333 0.66666666666667 0.82640591793881 69.723 > 113 123 Ga 0.58333333333333 0.66666666666667 0.82640591793881 69.723 > 113 124 Ga 0.83333333333333 0.66666666666667 0.82640591793881 69.723 > 113 125 Ga 0.08333333333333 0.91666666666667 0.82640591793881 69.723 > 113 126 Ga 0.33333333333333 0.91666666666667 0.82640591793881 69.723 > 113 127 Ga 0.58333333333333 0.91666666666667 0.82640591793881 69.723 > 113 128 Ga 0.83333333333333 0.91666666666667 0.82640591793881 69.723 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 7.3160896 0.0000000 0.0000000 -0.0000000 7.3160896 0.0000000 0.0000000 0.0000000 5.3074479 -------------------------- Born effective charges -------------------------- 1 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 2 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 3 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 4 S -2.3137276 -0.0000000 0.0000000 0.0000000 -2.3137276 0.0000000 0.0000000 0.0000000 -0.8873719 5 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 6 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 7 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 8 Ga 2.3137276 0.0000000 0.0000000 -0.0000000 2.3137276 0.0000000 0.0000000 0.0000000 0.8873719 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000002 (zzz) 0.00000002 (zzz) 0.00000002 (zzz) Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 16 16 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: 305, Lambda: 0.20 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/60) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 60 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 0.465 ( 0.000 0.000 0.000) 0.000 0.465 ( 0.000 0.000 0.000) 0.000 0.731 ( 0.000 0.000 0.000) 0.000 2.050 ( 0.000 0.000 0.000) 0.000 2.050 ( 0.000 0.000 0.000) 0.000 2.078 ( 0.000 0.000 0.000) 0.000 2.078 ( 0.000 0.000 0.000) 0.000 5.331 ( 0.000 0.000 0.000) 0.000 5.421 ( 0.000 0.000 0.000) 0.000 8.517 ( 0.000 0.000 0.000) 0.000 8.517 ( 0.000 0.000 0.000) 0.000 8.518 ( 0.000 0.000 0.000) 0.000 8.518 ( 0.000 0.000 0.000) 0.000 8.655 ( 0.000 0.000 0.000) 0.000 8.655 ( 0.000 0.000 0.000) 0.000 8.663 ( 0.000 0.000 0.000) 0.000 8.663 ( 0.000 0.000 0.000) 0.000 9.612 ( 0.000 0.000 0.000) 0.000 10.018 ( 0.000 0.000 0.000) 0.000 10.781 ( 0.000 0.000 0.000) 0.000 10.797 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/60) ======================= q-point: ( 0.06 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 144 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.302 ( 16.883 9.748 0.000) 19.495 0.592 ( 24.883 14.366 0.000) 28.732 0.757 ( 19.676 11.360 0.000) 22.719 0.770 ( 5.326 3.075 0.000) 6.149 1.089 ( 44.882 25.913 0.000) 51.826 1.221 ( 42.352 24.452 0.000) 48.904 2.143 ( 7.646 4.415 0.000) 8.829 2.173 ( 7.693 4.442 0.000) 8.883 2.369 ( 22.528 13.007 0.000) 26.014 2.391 ( 24.901 14.376 0.000) 28.753 5.305 ( -1.915 -1.106 0.000) 2.211 5.400 ( -1.624 -0.937 0.000) 1.875 8.520 ( 0.432 0.249 0.000) 0.499 8.530 ( 0.830 0.479 0.000) 0.958 8.545 ( 1.443 0.833 0.000) 1.666 8.575 ( 2.996 1.730 0.000) 3.460 8.661 ( 0.490 0.283 0.000) 0.566 8.672 ( 0.727 0.420 0.000) 0.839 8.952 ( 21.706 12.532 0.000) 25.064 9.600 ( -0.268 -0.155 0.000) 0.310 10.011 ( -0.209 -0.121 0.000) 0.242 10.428 ( 0.462 0.267 0.000) 0.533 10.720 ( -4.783 -2.761 0.000) 5.523 10.738 ( -4.840 -2.794 0.000) 5.589 ======================= Grid point 2 (3/60) ======================= q-point: ( 0.12 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.752 ( 20.581 11.883 0.000) 23.765 0.997 ( 13.048 7.533 0.000) 15.067 1.151 ( 22.673 13.091 0.000) 26.181 1.243 ( 20.721 11.963 0.000) 23.927 2.056 ( 36.797 21.245 0.000) 42.489 2.142 ( 34.624 19.990 0.000) 39.980 2.380 ( 11.916 6.879 0.000) 13.759 2.408 ( 11.687 6.748 0.000) 13.495 2.989 ( 28.345 16.365 0.000) 32.731 3.052 ( 29.307 16.920 0.000) 33.841 5.269 ( -0.399 -0.230 0.000) 0.461 5.361 ( -1.193 -0.689 0.000) 1.378 8.528 ( -3.578 -2.066 0.000) 4.131 8.537 ( 0.996 0.575 0.000) 1.150 8.555 ( 1.259 0.727 0.000) 1.454 8.574 ( -3.634 -2.098 0.000) 4.196 8.678 ( 0.877 0.506 0.000) 1.013 8.695 ( 1.133 0.654 0.000) 1.308 9.507 ( 22.650 13.077 0.000) 26.154 9.648 ( 4.995 2.884 0.000) 5.768 10.026 ( 1.527 0.881 0.000) 1.763 10.403 ( -4.389 -2.534 0.000) 5.068 10.588 ( -6.657 -3.844 0.000) 7.687 10.616 ( -2.442 -1.410 0.000) 2.820 ======================= Grid point 3 (4/60) ======================= q-point: ( 0.19 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 144 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.235 ( 20.268 11.702 0.000) 23.404 1.336 ( 15.374 8.876 0.000) 17.753 1.639 ( 18.740 10.819 0.000) 21.639 1.696 ( 17.678 10.207 0.000) 20.413 2.669 ( 12.231 7.061 0.000) 14.123 2.689 ( 11.898 6.869 0.000) 13.739 2.754 ( 21.197 12.238 0.000) 24.476 2.795 ( 21.457 12.388 0.000) 24.776 3.620 ( 24.607 14.207 0.000) 28.414 3.695 ( 24.855 14.350 0.000) 28.700 5.325 ( 5.884 3.397 0.000) 6.794 5.378 ( 3.416 1.972 0.000) 3.945 8.404 ( -5.605 -3.236 0.000) 6.472 8.445 ( -5.773 -3.333 0.000) 6.666 8.565 ( 1.362 0.786 0.000) 1.573 8.587 ( 1.443 0.833 0.000) 1.666 8.701 ( 1.107 0.639 0.000) 1.278 8.723 ( 1.223 0.706 0.000) 1.412 9.812 ( 7.917 4.571 0.000) 9.141 9.885 ( 6.866 3.964 0.000) 7.928 10.068 ( 1.692 0.977 0.000) 1.954 10.208 ( -10.874 -6.278 0.000) 12.556 10.495 ( 1.446 0.835 0.000) 1.670 10.642 ( 3.319 1.916 0.000) 3.832 ======================= Grid point 4 (5/60) ======================= q-point: ( 0.25 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.673 ( 16.186 9.345 0.000) 18.691 1.704 ( 16.369 9.451 0.000) 18.901 2.017 ( 13.483 7.785 0.000) 15.569 2.057 ( 13.055 7.538 0.000) 15.075 2.932 ( 9.916 5.725 0.000) 11.450 2.946 ( 9.807 5.662 0.000) 11.325 3.036 ( 4.532 2.617 0.000) 5.233 3.140 ( 7.606 4.391 0.000) 8.782 4.113 ( 17.452 10.076 0.000) 20.152 4.194 ( 17.708 10.223 0.000) 20.447 5.537 ( 10.787 6.228 0.000) 12.456 5.557 ( 12.669 7.314 0.000) 14.629 8.327 ( -0.112 -0.065 0.000) 0.129 8.370 ( 0.255 0.147 0.000) 0.294 8.599 ( 1.451 0.838 0.000) 1.676 8.621 ( 1.439 0.831 0.000) 1.662 8.728 ( 1.190 0.687 0.000) 1.374 8.751 ( 1.257 0.726 0.000) 1.451 9.832 ( -8.086 -4.669 0.000) 9.337 9.937 ( -11.790 -6.807 0.000) 13.614 9.972 ( 5.070 2.927 0.000) 5.855 10.095 ( 0.520 0.300 0.000) 0.601 10.656 ( 9.225 5.326 0.000) 10.652 10.737 ( 4.440 2.564 0.000) 5.127 ======================= Grid point 5 (6/60) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 144 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.999 ( 12.256 7.076 0.000) 14.152 2.066 ( 13.679 7.898 0.000) 15.795 2.266 ( 7.873 4.546 0.000) 9.092 2.303 ( 7.929 4.578 0.000) 9.155 3.036 ( -3.233 -1.867 0.000) 3.733 3.125 ( 6.587 3.803 0.000) 7.606 3.141 ( 6.801 3.927 0.000) 7.853 3.187 ( -2.325 -1.343 0.000) 2.685 4.451 ( 11.951 6.900 0.000) 13.800 4.536 ( 12.041 6.952 0.000) 13.904 5.831 ( 13.629 7.868 0.000) 15.737 5.894 ( 14.987 8.653 0.000) 17.306 8.412 ( 7.075 4.085 0.000) 8.169 8.468 ( 7.769 4.486 0.000) 8.971 8.630 ( 1.224 0.707 0.000) 1.413 8.654 ( 1.373 0.793 0.000) 1.585 8.756 ( 1.137 0.657 0.000) 1.313 8.782 ( 1.407 0.812 0.000) 1.624 9.609 ( -10.046 -5.800 0.000) 11.600 9.674 ( -10.426 -6.020 0.000) 12.039 10.029 ( -0.165 -0.095 0.000) 0.190 10.094 ( -0.683 -0.394 0.000) 0.788 10.830 ( 4.426 2.556 0.000) 5.111 10.848 ( 5.585 3.225 0.000) 6.449 ======================= Grid point 6 (7/60) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.247 ( 8.860 5.115 0.000) 10.230 2.330 ( 8.861 5.116 0.000) 10.232 2.394 ( 3.384 1.954 0.000) 3.907 2.433 ( 3.396 1.961 0.000) 3.921 2.936 ( -4.530 -2.615 0.000) 5.231 3.084 ( -5.513 -3.183 0.000) 6.365 3.243 ( 3.682 2.126 0.000) 4.252 3.264 ( 3.741 2.160 0.000) 4.319 4.691 ( 8.778 5.068 0.000) 10.136 4.775 ( 8.541 4.931 0.000) 9.863 6.137 ( 11.938 6.893 0.000) 13.785 6.219 ( 12.217 7.054 0.000) 14.107 8.627 ( 10.469 6.044 0.000) 12.089 8.655 ( 0.880 0.508 0.000) 1.017 8.684 ( 1.109 0.640 0.000) 1.280 8.699 ( 11.058 6.384 0.000) 12.769 8.781 ( 1.070 0.618 0.000) 1.236 8.815 ( 1.318 0.761 0.000) 1.522 9.389 ( -8.359 -4.826 0.000) 9.652 9.459 ( -7.810 -4.509 0.000) 9.018 9.984 ( -3.212 -1.855 0.000) 3.709 10.050 ( -3.168 -1.829 0.000) 3.658 10.908 ( 2.535 1.464 0.000) 2.928 10.945 ( 2.567 1.482 0.000) 2.964 ======================= Grid point 7 (8/60) ======================= q-point: ( 0.44 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 144 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.410 ( 4.920 2.841 0.000) 5.682 2.442 ( 1.016 0.587 0.000) 1.173 2.476 ( 0.626 0.361 0.000) 0.722 2.482 ( 4.172 2.409 0.000) 4.818 2.849 ( -2.545 -1.470 0.000) 2.939 2.963 ( -4.261 -2.460 0.000) 4.920 3.305 ( 1.714 0.990 0.000) 1.979 3.321 ( 1.320 0.762 0.000) 1.524 4.857 ( 5.112 2.952 0.000) 5.903 4.933 ( 4.812 2.778 0.000) 5.557 6.366 ( 7.191 4.152 0.000) 8.303 6.442 ( 6.557 3.786 0.000) 7.571 8.671 ( 0.495 0.286 0.000) 0.572 8.704 ( 0.598 0.345 0.000) 0.691 8.804 ( 0.744 0.430 0.000) 0.859 8.841 ( 0.792 0.457 0.000) 0.915 8.858 ( 8.261 4.769 0.000) 9.539 8.940 ( 8.547 4.935 0.000) 9.869 9.233 ( -4.679 -2.701 0.000) 5.403 9.316 ( -4.311 -2.489 0.000) 4.978 9.898 ( -3.645 -2.104 0.000) 4.209 9.960 ( -3.861 -2.229 0.000) 4.458 10.952 ( 1.219 0.704 0.000) 1.407 10.978 ( 0.526 0.304 0.000) 0.607 ======================= Grid point 8 (9/60) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 82 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.453 ( 0.000 0.000 0.000) 0.000 2.469 ( 0.000 0.000 0.000) 0.000 2.479 ( 0.000 0.000 0.000) 0.000 2.530 ( 0.000 0.000 0.000) 0.000 2.820 ( 0.000 0.000 0.000) 0.000 2.909 ( 0.000 0.000 0.000) 0.000 3.326 ( 0.000 0.000 0.000) 0.000 3.334 ( 0.000 0.000 0.000) 0.000 4.918 ( 0.000 0.000 0.000) 0.000 4.991 ( 0.000 0.000 0.000) 0.000 6.454 ( 0.000 0.000 0.000) 0.000 6.519 ( 0.000 0.000 0.000) 0.000 8.677 ( 0.000 0.000 0.000) 0.000 8.711 ( 0.000 0.000 0.000) 0.000 8.813 ( 0.000 0.000 0.000) 0.000 8.850 ( 0.000 0.000 0.000) 0.000 8.964 ( -0.000 -0.000 0.000) 0.000 9.050 ( 0.000 0.000 0.000) 0.000 9.178 ( 0.000 0.000 0.000) 0.000 9.264 ( 0.000 0.000 0.000) 0.000 9.848 ( 0.000 0.000 0.000) 0.000 9.907 ( 0.000 0.000 0.000) 0.000 10.967 ( 0.000 0.000 0.000) 0.000 10.981 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 18 (10/60) ======================= q-point: ( 0.06 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 144 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.619 ( 11.367 19.688 0.000) 22.734 0.915 ( 6.459 11.188 0.000) 12.918 1.058 ( 15.351 26.589 0.000) 30.703 1.194 ( 14.047 24.331 0.000) 28.095 1.793 ( 21.847 37.840 0.000) 43.693 1.856 ( 20.454 35.427 0.000) 40.908 2.327 ( 7.412 12.838 0.000) 14.824 2.357 ( 7.423 12.857 0.000) 14.846 2.801 ( 15.519 26.880 0.000) 31.039 2.849 ( 16.235 28.120 0.000) 32.471 5.273 ( -0.800 -1.385 0.000) 1.599 5.369 ( -0.987 -1.710 0.000) 1.974 8.517 ( -0.810 -1.403 0.000) 1.620 8.544 ( 0.450 0.780 0.000) 0.901 8.565 ( 0.286 0.496 0.000) 0.572 8.593 ( -0.559 -0.969 0.000) 1.119 8.678 ( 0.727 1.258 0.000) 1.453 8.693 ( 0.908 1.573 0.000) 1.816 9.357 ( 13.848 23.985 0.000) 27.696 9.620 ( 1.775 3.075 0.000) 3.550 10.017 ( 0.516 0.894 0.000) 1.033 10.437 ( -0.069 -0.119 0.000) 0.137 10.613 ( -3.992 -6.914 0.000) 7.984 10.636 ( -3.755 -6.504 0.000) 7.510 ======================= Grid point 19 (11/60) ======================= q-point: ( 0.12 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.050 ( 15.860 16.669 0.000) 23.009 1.196 ( 11.469 11.658 0.000) 16.354 1.577 ( 12.521 25.842 0.000) 28.715 1.660 ( 10.148 25.152 0.000) 27.122 2.498 ( 22.514 20.202 0.000) 30.249 2.526 ( 23.129 17.353 0.000) 28.915 2.622 ( 8.447 16.040 0.000) 18.128 2.647 ( 8.040 15.869 0.000) 17.789 3.374 ( 20.909 20.264 0.000) 29.117 3.441 ( 21.541 20.244 0.000) 29.560 5.281 ( 2.335 1.980 0.000) 3.061 5.354 ( 0.856 0.535 0.000) 1.009 8.437 ( -3.637 -5.621 0.000) 6.695 8.484 ( -4.050 -5.729 0.000) 7.016 8.581 ( 0.050 2.726 0.000) 2.726 8.597 ( 0.486 2.757 0.000) 2.799 8.709 ( 0.606 2.588 0.000) 2.658 8.731 ( 0.726 2.864 0.000) 2.955 9.736 ( 5.838 5.579 0.000) 8.075 9.799 ( 11.494 13.680 0.000) 17.868 10.045 ( 1.454 0.789 0.000) 1.654 10.334 ( -9.257 -5.090 0.000) 10.564 10.484 ( -2.310 -5.446 0.000) 5.916 10.574 ( 3.192 -1.894 0.000) 3.712 ======================= Grid point 20 (12/60) ======================= q-point: ( 0.19 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.489 ( 16.426 12.898 0.000) 20.885 1.531 ( 13.815 10.075 0.000) 17.099 2.002 ( 7.484 22.979 0.000) 24.167 2.048 ( 6.451 22.226 0.000) 23.143 2.918 ( 5.908 16.650 0.000) 17.667 2.919 ( 11.823 4.023 0.000) 12.488 2.934 ( 5.595 16.315 0.000) 17.247 2.988 ( 14.657 6.103 0.000) 15.876 3.899 ( 19.027 12.772 0.000) 22.916 3.973 ( 19.476 12.664 0.000) 23.231 5.421 ( 8.999 6.159 0.000) 10.905 5.442 ( 6.879 4.544 0.000) 8.244 8.330 ( -1.912 -3.190 0.000) 3.719 8.369 ( -1.997 -3.468 0.000) 4.002 8.617 ( 0.038 4.014 0.000) 4.014 8.639 ( 0.123 4.092 0.000) 4.093 8.749 ( 0.023 4.027 0.000) 4.027 8.773 ( -0.025 4.219 0.000) 4.219 9.898 ( 6.272 3.540 0.000) 7.202 9.919 ( -4.962 -0.966 0.000) 5.055 10.072 ( 1.367 -0.739 0.000) 1.554 10.079 ( -11.887 -6.631 0.000) 13.612 10.526 ( 8.954 2.179 0.000) 9.216 10.649 ( 5.991 -0.808 0.000) 6.045 ======================= Grid point 21 (13/60) ======================= q-point: ( 0.25 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.848 ( 13.076 8.008 0.000) 15.333 1.893 ( 15.082 8.976 0.000) 17.551 2.296 ( 2.416 17.870 0.000) 18.033 2.328 ( 2.424 17.403 0.000) 17.571 3.034 ( 1.687 -1.855 0.000) 2.507 3.153 ( 2.245 15.237 0.000) 15.401 3.162 ( 3.726 -1.581 0.000) 4.047 3.165 ( 2.427 15.164 0.000) 15.357 4.293 ( 14.322 7.748 0.000) 16.283 4.374 ( 14.693 7.677 0.000) 16.578 5.671 ( 12.377 6.869 0.000) 14.155 5.711 ( 14.382 7.615 0.000) 16.274 8.338 ( 4.172 1.410 0.000) 4.404 8.380 ( 4.870 0.993 0.000) 4.971 8.661 ( -0.434 5.096 0.000) 5.115 8.682 ( -0.481 5.136 0.000) 5.159 8.790 ( -0.749 5.285 0.000) 5.338 8.816 ( -0.693 5.491 0.000) 5.534 9.731 ( -10.354 -5.221 0.000) 11.596 9.805 ( -11.576 -6.260 0.000) 13.161 10.000 ( 2.959 -0.262 0.000) 2.970 10.077 ( 0.862 -2.069 0.000) 2.242 10.724 ( 9.350 1.294 0.000) 9.439 10.752 ( 6.529 -0.757 0.000) 6.572 ======================= Grid point 22 (14/60) ======================= q-point: ( 0.31 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.127 ( 10.738 5.769 0.000) 12.190 2.203 ( 11.531 5.805 0.000) 12.909 2.453 ( -1.288 11.632 0.000) 11.703 2.485 ( -1.100 11.575 0.000) 11.627 2.995 ( -3.547 -0.558 0.000) 3.591 3.135 ( -3.078 -2.526 0.000) 3.982 3.306 ( -1.107 13.362 0.000) 13.408 3.323 ( -0.843 13.373 0.000) 13.400 4.575 ( 10.627 5.474 0.000) 11.954 4.658 ( 10.662 5.279 0.000) 11.898 5.973 ( 13.584 6.015 0.000) 14.856 6.043 ( 14.633 6.026 0.000) 15.826 8.491 ( 9.377 3.531 0.000) 10.020 8.538 ( 9.780 2.555 0.000) 10.108 8.706 ( -0.716 6.269 0.000) 6.309 8.728 ( -0.633 6.418 0.000) 6.449 8.830 ( -1.466 6.445 0.000) 6.610 8.861 ( -1.222 6.690 0.000) 6.801 9.494 ( -9.915 -5.483 0.000) 11.330 9.558 ( -9.594 -5.435 0.000) 11.026 9.997 ( -0.932 -2.875 0.000) 3.023 10.056 ( -0.547 -3.199 0.000) 3.245 10.840 ( 5.288 -1.273 0.000) 5.439 10.869 ( 6.216 -1.112 0.000) 6.315 ======================= Grid point 23 (15/60) ======================= q-point: ( 0.38 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.336 ( 7.254 3.920 0.000) 8.245 2.411 ( 6.962 3.155 0.000) 7.644 2.505 ( -2.515 6.447 0.000) 6.920 2.535 ( -2.940 6.328 0.000) 6.978 2.929 ( -5.320 3.850 0.000) 6.566 3.040 ( -5.697 0.331 0.000) 5.707 3.391 ( -3.518 11.789 0.000) 12.302 3.409 ( -3.598 11.554 0.000) 12.101 4.781 ( 7.398 3.795 0.000) 8.315 4.860 ( 7.161 3.495 0.000) 7.968 6.240 ( 10.878 3.354 0.000) 11.383 6.317 ( 10.906 2.734 0.000) 11.243 8.671 ( 4.226 2.158 0.000) 4.745 8.694 ( 4.370 2.158 0.000) 4.874 8.781 ( 3.224 6.504 0.000) 7.259 8.840 ( 6.755 6.352 0.000) 9.273 8.870 ( -2.252 7.791 0.000) 8.109 8.907 ( -2.135 8.046 0.000) 8.324 9.290 ( -6.928 -4.862 0.000) 8.464 9.365 ( -6.283 -4.772 0.000) 7.890 9.929 ( -2.780 -3.367 0.000) 4.366 9.990 ( -2.644 -3.608 0.000) 4.473 10.902 ( 3.679 -1.922 0.000) 4.150 10.933 ( 3.435 -2.451 0.000) 4.220 ======================= Grid point 24 (16/60) ======================= q-point: ( 0.44 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.453 ( 2.627 1.584 0.000) 3.068 2.512 ( 2.330 0.813 0.000) 2.467 2.517 ( -1.987 4.348 0.000) 4.780 2.530 ( -2.957 3.631 0.000) 4.682 2.894 ( -5.044 7.194 0.000) 8.786 2.965 ( -4.809 3.862 0.000) 6.168 3.429 ( -4.990 10.666 0.000) 11.776 3.438 ( -5.381 10.270 0.000) 11.595 4.900 ( 2.867 1.244 0.000) 3.125 4.972 ( 2.751 1.020 0.000) 2.934 6.401 ( 5.256 -0.461 0.000) 5.277 6.466 ( 4.942 -1.233 0.000) 5.093 8.711 ( -1.212 3.430 0.000) 3.638 8.741 ( -0.955 3.124 0.000) 3.266 8.903 ( -4.015 9.231 0.000) 10.066 8.919 ( 4.837 1.060 0.000) 4.952 8.943 ( -4.134 9.453 0.000) 10.318 9.001 ( 5.271 0.812 0.000) 5.333 9.166 ( -1.610 -4.018 0.000) 4.328 9.251 ( -1.297 -4.049 0.000) 4.252 9.854 ( -1.730 -1.900 0.000) 2.569 9.916 ( -1.811 -1.778 0.000) 2.538 10.932 ( 2.355 -2.605 0.000) 3.512 10.951 ( 1.935 -2.814 0.000) 3.415 ======================= Grid point 35 (17/60) ======================= q-point: ( 0.12 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.403 ( 10.240 17.736 0.000) 20.480 1.455 ( 7.866 13.624 0.000) 15.731 2.093 ( 13.478 23.344 0.000) 26.955 2.149 ( 12.543 21.724 0.000) 25.085 2.800 ( 6.049 10.477 0.000) 12.098 2.853 ( 7.658 13.264 0.000) 15.316 2.974 ( 10.062 17.427 0.000) 20.123 2.991 ( 9.711 16.820 0.000) 19.422 3.777 ( 11.008 19.066 0.000) 22.016 3.844 ( 11.091 19.211 0.000) 22.183 5.371 ( 4.211 7.294 0.000) 8.422 5.404 ( 2.829 4.901 0.000) 5.659 8.330 ( -2.361 -4.090 0.000) 4.723 8.369 ( -2.687 -4.655 0.000) 5.375 8.644 ( 1.968 3.410 0.000) 3.937 8.665 ( 2.175 3.767 0.000) 4.350 8.775 ( 2.133 3.694 0.000) 4.265 8.802 ( 2.235 3.871 0.000) 4.470 9.855 ( 3.246 5.622 0.000) 6.492 9.965 ( 0.331 0.573 0.000) 0.662 10.052 ( -0.085 -0.146 0.000) 0.169 10.172 ( -6.101 -10.567 0.000) 12.202 10.442 ( 1.810 3.135 0.000) 3.620 10.570 ( 0.958 1.660 0.000) 1.917 ======================= Grid point 36 (18/60) ======================= q-point: ( 0.19 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.735 ( 9.871 10.938 0.000) 14.734 1.755 ( 11.925 12.537 0.000) 17.302 2.471 ( 3.219 20.601 0.000) 20.851 2.506 ( 2.896 20.973 0.000) 21.172 2.950 ( 5.548 0.617 0.000) 5.582 3.044 ( 6.911 0.351 0.000) 6.920 3.288 ( 4.407 18.543 0.000) 19.060 3.297 ( 4.579 18.233 0.000) 18.799 4.137 ( 12.514 11.036 0.000) 16.685 4.209 ( 12.922 10.911 0.000) 16.912 5.562 ( 7.805 7.622 0.000) 10.909 5.578 ( 9.501 8.935 0.000) 13.042 8.299 ( 1.374 0.416 0.000) 1.435 8.328 ( 1.384 -0.251 0.000) 1.407 8.717 ( 1.076 5.619 0.000) 5.721 8.743 ( 0.972 5.908 0.000) 5.987 8.851 ( 0.999 5.643 0.000) 5.731 8.880 ( 0.954 5.832 0.000) 5.910 9.850 ( -8.042 -5.585 0.000) 9.791 9.934 ( -9.795 -7.690 0.000) 12.453 9.945 ( 2.804 1.050 0.000) 2.994 10.041 ( 0.183 -2.196 0.000) 2.204 10.577 ( 7.697 2.628 0.000) 8.133 10.630 ( 4.857 -0.564 0.000) 4.890 ======================= Grid point 37 (19/60) ======================= q-point: ( 0.25 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 260 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.009 ( 10.125 8.185 0.000) 13.020 2.068 ( 11.258 8.550 0.000) 14.137 2.612 ( -3.185 8.848 0.000) 9.404 2.666 ( -3.733 12.497 0.000) 13.042 3.041 ( 1.396 5.866 0.000) 6.030 3.132 ( 2.718 1.353 0.000) 3.036 3.501 ( -1.072 18.133 0.000) 18.164 3.515 ( -0.771 18.252 0.000) 18.268 4.437 ( 11.463 6.747 0.000) 13.301 4.515 ( 11.839 6.465 0.000) 13.489 5.812 ( 11.417 7.109 0.000) 13.449 5.862 ( 12.683 7.400 0.000) 14.684 8.386 ( 6.108 3.304 0.000) 6.945 8.415 ( 6.793 2.499 0.000) 7.238 8.791 ( -0.241 7.351 0.000) 7.355 8.816 ( -0.439 7.610 0.000) 7.623 8.921 ( -0.286 6.711 0.000) 6.717 8.952 ( -0.175 7.090 0.000) 7.092 9.627 ( -10.509 -4.670 0.000) 11.500 9.689 ( -10.633 -4.858 0.000) 11.690 9.967 ( 0.643 -2.700 0.000) 2.775 10.017 ( 0.462 -3.522 0.000) 3.552 10.712 ( 6.422 -1.974 0.000) 6.719 10.724 ( 7.446 -1.813 0.000) 7.664 ======================= Grid point 38 (20/60) ======================= q-point: ( 0.31 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.247 ( 8.353 6.324 0.000) 10.476 2.315 ( 8.499 5.560 0.000) 10.156 2.593 ( -2.963 0.565 0.000) 3.016 2.662 ( -4.410 3.008 0.000) 5.339 3.108 ( -4.555 12.360 0.000) 13.173 3.170 ( -3.853 8.020 0.000) 8.897 3.625 ( -4.974 17.051 0.000) 17.761 3.642 ( -4.984 17.030 0.000) 17.744 4.676 ( 9.199 4.573 0.000) 10.273 4.754 ( 9.207 4.150 0.000) 10.100 6.078 ( 11.415 4.638 0.000) 12.321 6.143 ( 11.965 4.150 0.000) 12.665 8.556 ( 8.722 3.045 0.000) 9.238 8.587 ( 8.976 2.427 0.000) 9.299 8.855 ( -1.252 7.777 0.000) 7.877 8.883 ( -0.872 7.990 0.000) 8.037 8.978 ( -0.600 6.288 0.000) 6.316 9.019 ( -0.400 7.235 0.000) 7.247 9.408 ( -10.427 -1.690 0.000) 10.563 9.472 ( -10.052 -1.894 0.000) 10.229 9.926 ( -1.170 -3.799 0.000) 3.975 9.979 ( -0.621 -3.918 0.000) 3.967 10.790 ( 5.581 -3.441 0.000) 6.556 10.813 ( 5.972 -4.009 0.000) 7.193 ======================= Grid point 39 (21/60) ======================= q-point: ( 0.38 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 260 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.417 ( 4.562 4.302 0.000) 6.271 2.472 ( 4.246 3.165 0.000) 5.296 2.551 ( -0.985 -1.271 0.000) 1.608 2.607 ( -2.652 -0.345 0.000) 2.674 3.134 ( -7.817 14.769 0.000) 16.710 3.165 ( -7.837 11.863 0.000) 14.217 3.682 ( -7.258 15.886 0.000) 17.465 3.694 ( -7.599 15.637 0.000) 17.385 4.843 ( 5.539 2.312 0.000) 6.003 4.914 ( 5.416 1.813 0.000) 5.711 6.283 ( 7.964 1.211 0.000) 8.055 6.343 ( 7.758 0.235 0.000) 7.761 8.718 ( 5.877 1.886 0.000) 6.172 8.740 ( 5.141 1.771 0.000) 5.438 8.903 ( -0.942 5.521 0.000) 5.601 8.951 ( 1.120 4.837 0.000) 4.965 9.030 ( 0.948 3.482 0.000) 3.608 9.088 ( 0.801 5.836 0.000) 5.891 9.238 ( -9.812 3.577 0.000) 10.444 9.302 ( -9.329 2.049 0.000) 9.552 9.863 ( -1.734 -2.709 0.000) 3.217 9.925 ( -1.454 -2.437 0.000) 2.837 10.836 ( 4.254 -4.500 0.000) 6.193 10.857 ( 4.238 -4.794 0.000) 6.399 ======================= Grid point 40 (22/60) ======================= q-point: ( 0.44 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 144 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.481 ( -0.793 1.374 0.000) 1.586 2.525 ( -0.403 0.699 0.000) 0.807 2.537 ( 0.596 -1.032 0.000) 1.192 2.579 ( -0.092 0.159 0.000) 0.184 3.142 ( -8.975 15.545 0.000) 17.950 3.151 ( -7.739 13.405 0.000) 15.478 3.699 ( -8.615 14.921 0.000) 17.230 3.704 ( -8.668 15.013 0.000) 17.336 4.906 ( 0.391 -0.677 0.000) 0.781 4.972 ( 0.592 -1.026 0.000) 1.185 6.363 ( 1.668 -2.888 0.000) 3.335 6.412 ( 2.020 -3.499 0.000) 4.040 8.779 ( -1.013 1.754 0.000) 2.025 8.794 ( -0.146 0.252 0.000) 0.291 8.927 ( -0.913 1.582 0.000) 1.827 8.993 ( 0.049 -0.084 0.000) 0.097 9.070 ( 2.714 -4.700 0.000) 5.427 9.147 ( -1.167 2.022 0.000) 2.335 9.161 ( -4.230 7.327 0.000) 8.461 9.200 ( -8.461 14.655 0.000) 16.922 9.830 ( 0.263 -0.455 0.000) 0.525 9.896 ( 0.128 -0.221 0.000) 0.256 10.850 ( 3.062 -5.304 0.000) 6.125 10.869 ( 2.988 -5.175 0.000) 5.976 ======================= Grid point 53 (23/60) ======================= q-point: ( 0.19 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 144 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.958 ( 6.487 11.236 0.000) 12.974 2.009 ( 7.065 12.236 0.000) 14.129 2.693 ( 0.174 0.301 0.000) 0.348 2.783 ( 1.929 3.342 0.000) 3.859 3.064 ( 6.602 11.436 0.000) 13.205 3.096 ( 4.608 7.981 0.000) 9.215 3.628 ( 8.175 14.160 0.000) 16.350 3.630 ( 8.201 14.205 0.000) 16.402 4.357 ( 6.230 10.791 0.000) 12.461 4.427 ( 6.209 10.754 0.000) 12.418 5.744 ( 5.879 10.182 0.000) 11.758 5.783 ( 6.421 11.122 0.000) 12.843 8.350 ( 2.626 4.548 0.000) 5.252 8.365 ( 2.237 3.874 0.000) 4.473 8.835 ( 3.325 5.759 0.000) 6.650 8.865 ( 3.357 5.815 0.000) 6.715 8.958 ( 2.558 4.431 0.000) 5.117 8.991 ( 2.731 4.731 0.000) 5.463 9.712 ( -4.326 -7.492 0.000) 8.651 9.774 ( -4.530 -7.847 0.000) 9.061 9.943 ( -0.689 -1.194 0.000) 1.378 9.992 ( -1.392 -2.411 0.000) 2.784 10.629 ( 1.459 2.527 0.000) 2.918 10.636 ( 0.671 1.161 0.000) 1.341 ======================= Grid point 54 (24/60) ======================= q-point: ( 0.25 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.186 ( 7.290 9.264 0.000) 11.788 2.244 ( 7.390 8.842 0.000) 11.524 2.626 ( -3.025 -3.905 0.000) 4.939 2.729 ( -3.935 -3.806 0.000) 5.474 3.278 ( 1.816 13.700 0.000) 13.819 3.295 ( 2.785 12.101 0.000) 12.417 3.849 ( 0.270 15.295 0.000) 15.298 3.861 ( 0.909 15.074 0.000) 15.102 4.571 ( 8.001 6.580 0.000) 10.359 4.642 ( 8.215 6.207 0.000) 10.296 5.963 ( 8.228 7.827 0.000) 11.356 6.015 ( 8.667 7.790 0.000) 11.653 8.469 ( 5.145 4.704 0.000) 6.972 8.481 ( 6.114 4.053 0.000) 7.335 8.938 ( 0.961 6.683 0.000) 6.752 8.968 ( 0.853 6.799 0.000) 6.852 9.016 ( 0.257 1.653 0.000) 1.673 9.060 ( 0.743 2.551 0.000) 2.657 9.571 ( -6.665 -0.176 0.000) 6.667 9.631 ( -6.708 -0.265 0.000) 6.713 9.904 ( -0.867 -3.160 0.000) 3.277 9.949 ( -0.429 -2.860 0.000) 2.892 10.670 ( 3.449 -1.916 0.000) 3.945 10.679 ( 4.036 -2.389 0.000) 4.690 ======================= Grid point 55 (25/60) ======================= q-point: ( 0.31 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.386 ( 5.525 7.326 0.000) 9.176 2.435 ( 5.244 6.248 0.000) 8.157 2.545 ( -1.499 -3.559 0.000) 3.862 2.625 ( -2.647 -4.356 0.000) 5.098 3.394 ( -3.926 13.182 0.000) 13.754 3.412 ( -3.795 13.467 0.000) 13.991 3.958 ( -4.897 14.924 0.000) 15.707 3.973 ( -4.982 14.891 0.000) 15.702 4.759 ( 7.043 3.554 0.000) 7.889 4.825 ( 6.970 2.892 0.000) 7.547 6.170 ( 7.645 4.632 0.000) 8.939 6.223 ( 7.628 3.894 0.000) 8.564 8.614 ( 6.676 2.579 0.000) 7.156 8.636 ( 7.115 2.288 0.000) 7.474 8.996 ( -1.598 5.711 0.000) 5.930 9.005 ( 0.037 -3.333 0.000) 3.334 9.027 ( -1.382 5.854 0.000) 6.015 9.069 ( 0.623 -2.706 0.000) 2.777 9.487 ( -8.019 8.971 0.000) 12.032 9.544 ( -8.180 8.830 0.000) 12.036 9.859 ( -0.868 -2.374 0.000) 2.528 9.918 ( -0.405 -1.823 0.000) 1.868 10.707 ( 4.335 -4.601 0.000) 6.322 10.721 ( 4.825 -4.840 0.000) 6.834 ======================= Grid point 56 (26/60) ======================= q-point: ( 0.38 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.506 ( 0.432 -2.399 0.000) 2.437 2.509 ( 1.130 4.615 0.000) 4.751 2.542 ( 1.046 3.644 0.000) 3.791 2.560 ( -0.205 -3.172 0.000) 3.178 3.426 ( -6.775 12.761 0.000) 14.448 3.451 ( -7.338 14.279 0.000) 16.054 3.994 ( -7.291 14.097 0.000) 15.871 4.006 ( -7.676 14.447 0.000) 16.360 4.873 ( 3.153 0.622 0.000) 3.214 4.931 ( 3.254 -0.166 0.000) 3.258 6.303 ( 3.628 0.964 0.000) 3.754 6.343 ( 3.466 0.033 0.000) 3.466 8.732 ( 4.835 -0.500 0.000) 4.861 8.747 ( 3.987 -1.120 0.000) 4.142 8.974 ( 1.898 -5.220 0.000) 5.554 9.004 ( -2.850 4.473 0.000) 5.304 9.042 ( -2.313 4.105 0.000) 4.712 9.045 ( 2.293 -5.808 0.000) 6.244 9.460 ( -8.614 13.982 0.000) 16.422 9.514 ( -8.979 14.371 0.000) 16.945 9.831 ( -0.630 -0.356 0.000) 0.724 9.899 ( -0.485 -0.147 0.000) 0.507 10.727 ( 4.030 -6.024 0.000) 7.247 10.746 ( 4.262 -6.074 0.000) 7.420 ======================= Grid point 70 (27/60) ======================= q-point: ( 0.25 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.375 ( 5.206 9.018 0.000) 10.413 2.420 ( 4.714 8.165 0.000) 9.428 2.539 ( -2.236 -3.873 0.000) 4.472 2.625 ( -3.053 -5.288 0.000) 6.106 3.500 ( 4.517 7.825 0.000) 9.035 3.520 ( 5.182 8.976 0.000) 10.365 4.092 ( 4.875 8.444 0.000) 9.750 4.107 ( 5.144 8.910 0.000) 10.289 4.704 ( 3.726 6.453 0.000) 7.452 4.767 ( 3.465 6.002 0.000) 6.930 6.125 ( 4.547 7.876 0.000) 9.094 6.173 ( 4.317 7.478 0.000) 8.634 8.566 ( 2.779 4.814 0.000) 5.558 8.579 ( 3.078 5.331 0.000) 6.156 8.986 ( -2.097 -3.632 0.000) 4.194 9.043 ( 1.969 3.411 0.000) 3.939 9.046 ( -1.813 -3.140 0.000) 3.626 9.074 ( 2.011 3.483 0.000) 4.022 9.618 ( 2.238 3.877 0.000) 4.477 9.677 ( 2.290 3.966 0.000) 4.580 9.853 ( -0.950 -1.646 0.000) 1.901 9.910 ( -0.539 -0.934 0.000) 1.079 10.635 ( -0.937 -1.623 0.000) 1.874 10.640 ( -0.677 -1.172 0.000) 1.354 ======================= Grid point 71 (28/60) ======================= q-point: ( 0.31 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.483 ( -0.584 -2.253 0.000) 2.327 2.530 ( 3.012 6.422 0.000) 7.093 2.537 ( -1.444 -3.817 0.000) 4.081 2.558 ( 2.616 5.444 0.000) 6.040 3.604 ( -1.063 7.487 0.000) 7.562 3.647 ( -0.816 8.917 0.000) 8.954 4.201 ( -1.795 8.786 0.000) 8.967 4.223 ( -1.901 9.399 0.000) 9.589 4.822 ( 3.841 2.772 0.000) 4.737 4.872 ( 3.530 1.700 0.000) 3.918 6.264 ( 3.714 4.432 0.000) 5.783 6.300 ( 3.301 3.490 0.000) 4.804 8.659 ( 3.554 1.838 0.000) 4.001 8.676 ( 3.378 1.580 0.000) 3.729 8.912 ( -0.733 -4.550 0.000) 4.608 8.971 ( -0.810 -5.534 0.000) 5.593 9.079 ( -0.945 2.586 0.000) 2.753 9.110 ( -0.924 2.306 0.000) 2.484 9.691 ( -2.330 9.673 0.000) 9.949 9.755 ( -2.406 10.445 0.000) 10.719 9.838 ( -0.145 0.184 0.000) 0.235 9.905 ( 0.005 0.320 0.000) 0.320 10.612 ( 1.801 -4.448 0.000) 4.799 10.624 ( 2.088 -4.614 0.000) 5.064 ======================= Grid point 72 (29/60) ======================= q-point: ( 0.38 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.469 ( 0.676 -1.170 0.000) 1.351 2.503 ( 1.257 -2.178 0.000) 2.515 2.592 ( -1.871 3.240 0.000) 3.742 2.609 ( -1.513 2.621 0.000) 3.026 3.630 ( -4.001 6.930 0.000) 8.002 3.683 ( -4.638 8.033 0.000) 9.276 4.226 ( -4.818 8.346 0.000) 9.637 4.251 ( -5.205 9.016 0.000) 10.411 4.873 ( 0.194 -0.335 0.000) 0.387 4.912 ( 0.828 -1.434 0.000) 1.656 6.324 ( -0.587 1.017 0.000) 1.174 6.347 ( -0.173 0.300 0.000) 0.346 8.710 ( 1.025 -1.775 0.000) 2.050 8.711 ( 0.832 -1.442 0.000) 1.665 8.882 ( 1.890 -3.274 0.000) 3.781 8.932 ( 2.721 -4.714 0.000) 5.443 9.081 ( -1.679 2.908 0.000) 3.358 9.108 ( -1.273 2.204 0.000) 2.545 9.721 ( -6.091 10.549 0.000) 12.181 9.789 ( -6.638 11.497 0.000) 13.276 9.840 ( -0.476 0.825 0.000) 0.953 9.908 ( -0.393 0.682 0.000) 0.787 10.606 ( 3.106 -5.380 0.000) 6.212 10.619 ( 3.439 -5.957 0.000) 6.878 ======================= Grid point 88 (30/60) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 144 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.458 ( -0.124 -0.214 0.000) 0.247 2.479 ( -1.046 -1.812 0.000) 2.093 2.628 ( 1.606 2.781 0.000) 3.211 2.637 ( 1.169 2.024 0.000) 2.337 3.700 ( 1.296 2.245 0.000) 2.593 3.765 ( 1.645 2.849 0.000) 3.290 4.314 ( 1.520 2.633 0.000) 3.040 4.347 ( 1.672 2.897 0.000) 3.345 4.870 ( 1.001 1.734 0.000) 2.002 4.892 ( 0.112 0.193 0.000) 0.223 6.336 ( 1.310 2.269 0.000) 2.620 6.351 ( 0.676 1.170 0.000) 1.351 8.687 ( 0.548 0.949 0.000) 1.096 8.696 ( 0.161 0.279 0.000) 0.322 8.851 ( -0.680 -1.178 0.000) 1.361 8.878 ( -1.782 -3.086 0.000) 3.563 9.113 ( 0.601 1.041 0.000) 1.202 9.129 ( -0.026 -0.044 0.000) 0.051 9.838 ( 2.377 4.118 0.000) 4.755 9.850 ( 0.326 0.565 0.000) 0.653 9.917 ( 0.272 0.471 0.000) 0.544 9.920 ( 2.778 4.812 0.000) 5.557 10.542 ( -1.178 -2.040 0.000) 2.355 10.546 ( -1.436 -2.487 0.000) 2.871 ======================= Grid point 271 (31/60) ======================= 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 60 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.230 ( -0.000 0.000 10.018) 10.018 0.230 ( 0.000 0.000 10.018) 10.018 0.362 ( 0.000 0.000 15.778) 15.778 0.401 ( 0.000 0.000 -5.916) 5.916 0.401 ( 0.000 0.000 -5.916) 5.916 0.631 ( 0.000 0.000 -9.294) 9.294 2.057 ( 0.000 0.000 0.627) 0.627 2.057 ( 0.000 0.000 0.627) 0.627 2.071 ( 0.000 0.000 -0.589) 0.589 2.071 ( 0.000 0.000 -0.589) 0.589 5.354 ( -0.000 0.000 1.953) 1.953 5.399 ( -0.000 0.000 -1.921) 1.921 8.518 ( 0.000 0.000 -0.001) 0.001 8.518 ( -0.000 0.000 -0.001) 0.001 8.518 ( -0.000 0.000 -0.001) 0.001 8.518 ( -0.000 0.000 -0.001) 0.001 8.657 ( 0.000 0.000 0.168) 0.168 8.657 ( 0.000 0.000 0.168) 0.168 8.661 ( -0.000 0.000 -0.165) 0.165 8.661 ( -0.000 -0.000 -0.165) 0.165 9.975 ( 0.000 0.000 1.245) 1.245 10.004 ( 0.000 -0.000 -1.262) 1.262 10.785 ( 0.000 0.000 0.368) 0.368 10.793 ( 0.000 -0.000 -0.342) 0.342 ======================= Grid point 272 (32/60) ======================= q-point: ( 0.06 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 208 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.455 ( 10.822 6.248 11.180) 16.767 0.636 ( 23.214 13.402 3.721) 27.062 0.675 ( 6.652 3.841 -8.484) 11.445 0.718 ( 20.674 11.936 -3.400) 24.113 1.123 ( 43.943 25.370 2.895) 50.824 1.188 ( 42.723 24.666 -2.738) 49.408 2.151 ( 7.658 4.421 0.658) 8.867 2.166 ( 7.681 4.435 -0.619) 8.891 2.385 ( 23.168 13.376 0.828) 26.765 2.392 ( 24.310 14.035 0.008) 28.070 5.329 ( -1.848 -1.067 2.052) 2.960 5.376 ( -1.690 -0.976 -2.060) 2.838 8.522 ( 0.524 0.302 0.212) 0.641 8.527 ( 0.722 0.417 -0.226) 0.864 8.553 ( 1.762 1.017 0.679) 2.145 8.568 ( 2.533 1.463 -0.641) 2.995 8.664 ( 0.558 0.322 0.243) 0.688 8.670 ( 0.675 0.390 -0.230) 0.813 9.040 ( 24.694 14.257 8.979) 29.895 9.356 ( 21.071 12.165 -19.841) 31.395 10.013 ( 1.248 0.721 0.618) 1.568 10.099 ( 9.241 5.335 -11.164) 15.443 10.724 ( -4.912 -2.836 0.411) 5.687 10.733 ( -4.857 -2.804 -0.409) 5.623 ======================= Grid point 273 (33/60) ======================= q-point: ( 0.12 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 209 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.814 ( 18.202 10.509 5.325) 21.682 0.936 ( 14.538 8.394 -5.263) 17.593 1.174 ( 22.158 12.793 1.989) 25.663 1.220 ( 21.184 12.230 -1.973) 24.540 2.077 ( 36.301 20.959 1.828) 41.957 2.120 ( 35.236 20.343 -1.875) 40.730 2.387 ( 11.854 6.844 0.612) 13.702 2.401 ( 11.740 6.778 -0.578) 13.569 3.010 ( 28.364 16.376 1.643) 32.793 3.041 ( 28.906 16.689 -1.070) 33.395 5.292 ( -0.611 -0.353 1.980) 2.102 5.338 ( -1.005 -0.580 -1.996) 2.309 8.539 ( -3.565 -2.058 1.019) 4.241 8.541 ( 1.050 0.606 0.368) 1.267 8.550 ( 1.180 0.681 -0.410) 1.422 8.563 ( -3.592 -2.074 -0.980) 4.262 8.682 ( 0.955 0.551 0.389) 1.169 8.691 ( 1.080 0.624 -0.347) 1.295 9.558 ( 18.091 10.445 4.098) 21.288 9.630 ( 8.469 4.890 -2.042) 9.990 10.076 ( 3.717 2.146 5.402) 6.899 10.264 ( 3.000 1.732 -11.050) 11.580 10.587 ( -6.099 -3.521 0.059) 7.042 10.596 ( -5.066 -2.925 -0.843) 5.910 ======================= Grid point 274 (34/60) ======================= q-point: ( 0.19 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 208 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.257 ( 19.035 10.990 2.007) 22.072 1.308 ( 16.591 9.579 -2.395) 19.307 1.653 ( 18.473 10.665 1.239) 21.367 1.682 ( 17.943 10.359 -1.249) 20.756 2.674 ( 12.143 7.011 0.450) 14.029 2.684 ( 11.977 6.915 -0.428) 13.836 2.765 ( 21.306 12.301 0.947) 24.620 2.786 ( 21.429 12.372 -0.832) 24.758 3.641 ( 24.592 14.198 1.735) 28.449 3.679 ( 24.725 14.275 -1.509) 28.590 5.338 ( 5.253 3.033 1.105) 6.165 5.364 ( 4.018 2.320 -1.171) 4.785 8.415 ( -5.653 -3.264 0.930) 6.594 8.436 ( -5.737 -3.313 -0.868) 6.682 8.570 ( 1.375 0.794 0.440) 1.648 8.581 ( 1.414 0.816 -0.511) 1.710 8.707 ( 1.145 0.661 0.501) 1.414 8.718 ( 1.201 0.694 -0.430) 1.452 9.825 ( 7.661 4.423 1.267) 8.936 9.860 ( 7.186 4.149 -1.832) 8.497 10.117 ( -1.936 -1.118 3.716) 4.336 10.185 ( -8.431 -4.868 -2.271) 9.997 10.520 ( 2.779 1.604 2.415) 4.016 10.592 ( 4.019 2.320 -3.912) 6.070 ======================= Grid point 275 (35/60) ======================= q-point: ( 0.25 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 209 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.678 ( 16.288 9.404 0.522) 18.815 1.694 ( 16.387 9.461 -0.802) 18.939 2.027 ( 13.377 7.723 0.865) 15.471 2.047 ( 13.163 7.600 -0.875) 15.225 2.935 ( 9.887 5.708 0.321) 11.421 2.943 ( 9.833 5.677 -0.308) 11.358 3.063 ( 5.256 3.035 2.290) 6.487 3.115 ( 6.794 3.923 -2.194) 8.146 4.135 ( 17.518 10.114 1.817) 20.310 4.175 ( 17.645 10.188 -1.668) 20.443 5.541 ( 11.278 6.512 0.409) 13.030 5.551 ( 12.220 7.055 -0.461) 14.118 8.338 ( -0.035 -0.020 0.953) 0.954 8.360 ( 0.149 0.086 -0.921) 0.937 8.603 ( 1.445 0.834 0.442) 1.726 8.615 ( 1.438 0.830 -0.525) 1.742 8.735 ( 1.211 0.699 0.541) 1.499 8.746 ( 1.243 0.718 -0.457) 1.506 9.847 ( -7.664 -4.425 1.414) 8.962 9.889 ( -7.144 -4.124 -2.541) 8.631 10.025 ( 0.199 0.115 2.963) 2.972 10.076 ( 0.293 0.169 -1.767) 1.799 10.675 ( 8.110 4.682 1.662) 9.510 10.715 ( 5.716 3.300 -1.816) 6.846 ======================= Grid point 276 (36/60) ======================= q-point: ( 0.31 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 208 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.015 ( 12.648 7.302 1.375) 14.669 2.048 ( 13.357 7.712 -1.521) 15.498 2.275 ( 7.886 4.553 0.784) 9.140 2.294 ( 7.914 4.569 -0.794) 9.173 3.073 ( -3.028 -1.748 3.272) 4.788 3.129 ( 6.642 3.835 0.358) 7.678 3.137 ( 6.749 3.897 -0.345) 7.801 3.149 ( -2.573 -1.485 -3.280) 4.425 4.474 ( 11.970 6.911 1.925) 13.955 4.516 ( 12.015 6.937 -1.766) 13.986 5.846 ( 13.981 8.072 1.364) 16.201 5.878 ( 14.660 8.464 -1.364) 16.982 8.426 ( 7.228 4.173 1.221) 8.435 8.454 ( 7.574 4.373 -1.235) 8.832 8.635 ( 1.247 0.720 0.456) 1.511 8.647 ( 1.317 0.760 -0.562) 1.621 8.763 ( 1.223 0.706 0.634) 1.548 8.777 ( 1.353 0.781 -0.528) 1.649 9.624 ( -10.027 -5.789 1.369) 11.658 9.657 ( -10.198 -5.888 -1.463) 11.867 10.047 ( -0.444 -0.256 1.458) 1.545 10.079 ( -0.686 -0.396 -1.338) 1.555 10.835 ( 4.761 2.749 0.376) 5.511 10.843 ( 5.338 3.082 -0.387) 6.176 ======================= Grid point 277 (37/60) ======================= q-point: ( 0.38 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 209 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.267 ( 8.872 5.122 1.757) 10.394 2.309 ( 8.870 5.121 -1.847) 10.408 2.404 ( 3.387 1.955 0.827) 3.997 2.423 ( 3.393 1.959 -0.840) 4.006 2.973 ( -4.768 -2.753 3.215) 6.376 3.047 ( -5.258 -3.036 -3.232) 6.878 3.248 ( 3.697 2.134 0.451) 4.293 3.259 ( 3.726 2.151 -0.435) 4.325 4.713 ( 8.699 5.023 1.868) 10.217 4.755 ( 8.583 4.955 -1.743) 10.063 6.158 ( 12.007 6.932 1.785) 13.979 6.199 ( 12.146 7.013 -1.772) 14.137 8.644 ( 10.569 6.102 1.514) 12.298 8.661 ( 0.921 0.532 0.536) 1.191 8.675 ( 1.029 0.594 -0.691) 1.374 8.680 ( 10.856 6.268 -1.616) 12.639 8.791 ( 1.155 0.667 0.802) 1.556 8.808 ( 1.272 0.735 -0.647) 1.605 9.407 ( -8.157 -4.709 1.571) 9.549 9.442 ( -7.889 -4.555 -1.471) 9.228 10.000 ( -3.212 -1.854 1.395) 3.962 10.033 ( -3.189 -1.841 -1.435) 3.952 10.918 ( 2.545 1.470 0.820) 3.051 10.936 ( 2.561 1.479 -0.770) 3.056 ======================= Grid point 278 (38/60) ======================= q-point: ( 0.44 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 208 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.427 ( 4.742 2.738 1.530) 5.686 2.450 ( 0.921 0.532 0.729) 1.290 2.463 ( 4.369 2.522 -1.596) 5.291 2.467 ( 0.726 0.419 -0.738) 1.117 2.878 ( -2.974 -1.717 2.471) 4.231 2.935 ( -3.829 -2.211 -2.471) 5.066 3.309 ( 1.612 0.931 0.363) 1.897 3.317 ( 1.415 0.817 -0.352) 1.671 4.876 ( 5.024 2.901 1.693) 6.043 4.915 ( 4.875 2.814 -1.614) 5.856 6.385 ( 7.031 4.059 1.666) 8.287 6.423 ( 6.714 3.876 -1.659) 7.928 8.678 ( 0.517 0.299 0.607) 0.851 8.694 ( 0.567 0.328 -0.792) 1.028 8.815 ( 0.761 0.440 0.888) 1.250 8.833 ( 0.784 0.453 -0.704) 1.147 8.876 ( 8.210 4.740 1.612) 9.616 8.916 ( 8.297 4.790 -1.949) 9.776 9.257 ( -4.411 -2.547 1.942) 5.451 9.298 ( -4.283 -2.473 -1.602) 5.199 9.913 ( -3.690 -2.131 1.335) 4.466 9.945 ( -3.797 -2.192 -1.363) 4.592 10.958 ( 1.037 0.598 0.590) 1.334 10.972 ( 0.691 0.399 -0.557) 0.973 ======================= Grid point 279 (39/60) ======================= q-point: (-0.50 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 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.459 ( 0.000 0.000 0.579) 0.579 2.473 ( 0.000 0.000 -0.581) 0.581 2.484 ( 0.000 0.000 1.287) 1.287 2.514 ( 0.000 0.000 -1.335) 1.335 2.843 ( 0.000 0.000 1.917) 1.917 2.887 ( 0.000 0.000 -1.920) 1.920 3.328 ( 0.000 0.000 0.170) 0.170 3.332 ( 0.000 0.000 -0.165) 0.165 4.937 ( 0.000 0.000 1.598) 1.598 4.973 ( 0.000 0.000 -1.536) 1.536 6.470 ( 0.000 0.000 1.416) 1.416 6.503 ( 0.000 0.000 -1.410) 1.410 8.684 ( 0.000 0.000 0.628) 0.628 8.701 ( 0.000 0.000 -0.816) 0.816 8.824 ( 0.000 0.000 0.894) 0.894 8.842 ( 0.000 0.000 -0.707) 0.707 8.980 ( 0.000 0.000 1.492) 1.492 9.021 ( 0.000 0.000 -2.131) 2.131 9.207 ( 0.000 0.000 2.140) 2.140 9.248 ( 0.000 0.000 -1.502) 1.502 9.862 ( 0.000 0.000 1.273) 1.273 9.892 ( 0.000 0.000 -1.288) 1.288 10.970 ( 0.000 0.000 0.334) 0.334 10.978 ( 0.000 0.000 -0.314) 0.314 ======================= Grid point 289 (40/60) ======================= q-point: ( 0.06 0.06 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 208 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.697 ( 9.657 16.726 6.565) 20.398 0.843 ( 7.350 12.731 -6.201) 15.954 1.091 ( 15.024 26.022 2.794) 30.177 1.158 ( 14.419 24.974 -3.054) 28.999 1.810 ( 21.485 37.213 1.488) 42.995 1.842 ( 20.775 35.983 -1.288) 41.569 2.336 ( 7.423 12.857 0.697) 14.862 2.350 ( 7.429 12.867 -0.580) 14.869 2.819 ( 15.542 26.919 1.357) 31.114 2.842 ( 15.956 27.636 -0.723) 31.920 5.297 ( -0.849 -1.470 2.072) 2.678 5.345 ( -0.940 -1.628 -2.079) 2.803 8.525 ( -0.507 -0.878 0.682) 1.222 8.539 ( 0.155 0.268 -0.501) 0.589 8.570 ( 0.096 0.167 0.517) 0.552 8.584 ( -0.356 -0.617 -0.713) 1.008 8.682 ( 0.778 1.347 0.342) 1.593 8.689 ( 0.869 1.504 -0.317) 1.766 9.435 ( 11.801 20.440 6.756) 24.550 9.574 ( 4.969 8.606 -4.794) 11.033 10.052 ( 1.900 3.292 4.221) 5.680 10.240 ( 3.622 6.274 -13.360) 15.198 10.619 ( -3.980 -6.894 0.488) 7.975 10.630 ( -3.918 -6.786 -0.484) 7.850 ======================= Grid point 290 (41/60) ======================= q-point: ( 0.12 0.06 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.083 ( 14.716 15.335 2.972) 21.461 1.156 ( 12.529 12.840 -3.354) 18.251 1.597 ( 11.996 25.635 1.756) 28.357 1.638 ( 10.820 25.293 -1.852) 27.572 2.506 ( 22.639 19.571 0.688) 29.934 2.520 ( 22.882 18.126 -0.552) 29.196 2.629 ( 8.345 16.022 0.602) 18.075 2.642 ( 8.125 15.932 -0.498) 17.891 3.393 ( 21.011 20.134 1.569) 29.143 3.427 ( 21.334 20.140 -1.308) 29.367 5.299 ( 1.948 1.606 1.571) 2.973 5.336 ( 1.209 0.883 -1.611) 2.199 8.449 ( -3.722 -5.637 1.047) 6.836 8.472 ( -3.927 -5.690 -1.000) 6.985 8.585 ( 0.135 2.709 0.310) 2.730 8.593 ( 0.348 2.718 -0.361) 2.764 8.715 ( 0.649 2.678 0.495) 2.800 8.726 ( 0.706 2.813 -0.429) 2.931 9.750 ( 6.455 6.517 1.214) 9.253 9.780 ( 8.855 9.956 -1.483) 13.406 10.124 ( 1.267 1.948 6.755) 7.143 10.274 ( -4.891 -1.893 -5.899) 7.893 10.492 ( -1.520 -4.715 0.927) 5.040 10.530 ( 1.642 -2.588 -2.766) 4.128 ======================= Grid point 291 (42/60) ======================= q-point: ( 0.19 0.06 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.496 ( 15.818 12.235 0.732) 20.011 1.517 ( 14.519 10.825 -1.079) 18.142 2.013 ( 7.251 22.787 0.982) 23.933 2.036 ( 6.734 22.411 -1.006) 23.423 2.921 ( 6.297 15.718 0.351) 16.936 2.929 ( 6.401 15.127 -0.229) 16.427 2.939 ( 11.397 6.481 1.422) 13.188 2.972 ( 13.777 5.761 -1.360) 14.995 3.918 ( 19.139 12.713 1.678) 23.038 3.956 ( 19.364 12.662 -1.545) 23.188 5.426 ( 8.473 5.767 0.408) 10.258 5.436 ( 7.412 4.959 -0.470) 8.930 8.340 ( -1.941 -3.264 0.872) 3.896 8.360 ( -1.984 -3.403 -0.820) 4.023 8.622 ( 0.060 4.013 0.404) 4.033 8.632 ( 0.104 4.047 -0.505) 4.080 8.756 ( 0.006 4.095 0.570) 4.135 8.768 ( -0.017 4.188 -0.475) 4.215 9.901 ( 3.577 2.402 0.356) 4.323 9.913 ( -1.315 0.389 -0.584) 1.491 10.076 ( -4.360 -3.160 0.423) 5.402 10.082 ( -7.565 -4.724 -0.004) 8.919 10.553 ( 8.567 1.690 2.463) 9.073 10.615 ( 7.064 0.188 -2.845) 7.618 ======================= Grid point 292 (43/60) ======================= q-point: ( 0.25 0.06 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.857 ( 13.620 8.288 0.880) 15.968 1.880 ( 14.621 8.770 -1.076) 17.084 2.304 ( 2.420 17.779 0.684) 17.956 2.320 ( 2.424 17.544 -0.671) 17.723 3.066 ( 2.169 -1.795 2.754) 3.938 3.129 ( 3.190 -1.142 -2.701) 4.333 3.157 ( 2.266 14.777 0.222) 14.952 3.162 ( 2.368 15.048 -0.236) 15.235 4.314 ( 14.431 7.718 1.826) 16.467 4.355 ( 14.615 7.683 -1.688) 16.598 5.681 ( 12.894 7.070 0.860) 14.730 5.701 ( 13.896 7.443 -0.877) 15.788 8.349 ( 4.335 1.308 0.914) 4.619 8.370 ( 4.685 1.100 -0.885) 4.893 8.665 ( -0.448 5.081 0.401) 5.117 8.676 ( -0.471 5.097 -0.521) 5.146 8.797 ( -0.737 5.354 0.605) 5.438 8.810 ( -0.709 5.453 -0.503) 5.522 9.747 ( -10.313 -5.304 1.368) 11.677 9.782 ( -10.742 -5.732 -1.772) 12.304 10.024 ( 1.891 -0.996 1.871) 2.840 10.061 ( 1.022 -1.810 -1.418) 2.516 10.731 ( 8.688 0.812 0.562) 8.744 10.744 ( 7.278 -0.214 -0.634) 7.309 ======================= Grid point 293 (44/60) ======================= q-point: ( 0.31 0.06 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.145 ( 10.952 5.791 1.590) 12.491 2.183 ( 11.346 5.807 -1.701) 12.859 2.461 ( -1.255 11.687 0.719) 11.777 2.477 ( -1.160 11.656 -0.667) 11.733 3.029 ( -3.420 -1.115 3.008) 4.689 3.099 ( -3.184 -2.087 -3.084) 4.900 3.311 ( -1.047 13.365 0.375) 13.412 3.319 ( -0.911 13.371 -0.352) 13.407 4.597 ( 10.631 5.406 1.879) 12.074 4.638 ( 10.649 5.310 -1.746) 12.027 5.991 ( 13.850 6.019 1.538) 15.179 6.026 ( 14.374 6.025 -1.526) 15.660 8.503 ( 9.517 3.312 1.060) 10.133 8.527 ( 9.733 2.824 -1.012) 10.185 8.710 ( -0.791 6.255 0.360) 6.315 8.721 ( -0.751 6.317 -0.586) 6.388 8.839 ( -1.381 6.533 0.735) 6.718 8.854 ( -1.268 6.647 -0.599) 6.794 9.510 ( -9.763 -5.457 1.385) 11.270 9.542 ( -9.600 -5.432 -1.388) 11.117 10.012 ( -0.899 -2.968 1.296) 3.361 10.042 ( -0.703 -3.129 -1.276) 3.451 10.847 ( 5.544 -1.227 0.639) 5.714 10.862 ( 6.006 -1.147 -0.611) 6.145 ======================= Grid point 294 (45/60) ======================= q-point: ( 0.38 0.06 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.354 ( 7.187 3.730 1.585) 8.251 2.391 ( 7.036 3.347 -1.673) 7.969 2.514 ( -2.654 6.521 0.708) 7.076 2.529 ( -2.864 6.453 -0.600) 7.085 2.956 ( -5.373 2.877 2.349) 6.532 3.012 ( -5.564 1.128 -2.463) 6.188 3.396 ( -3.541 11.734 0.398) 12.263 3.405 ( -3.580 11.616 -0.375) 12.161 4.802 ( 7.327 3.703 1.755) 8.395 4.841 ( 7.209 3.554 -1.663) 8.208 6.260 ( 10.883 3.196 1.672) 11.465 6.298 ( 10.897 2.886 -1.664) 11.394 8.678 ( 3.985 2.128 0.561) 4.553 8.690 ( 4.056 2.106 -0.436) 4.591 8.791 ( 4.260 6.480 1.003) 7.820 8.820 ( 5.861 6.308 -1.497) 8.740 8.882 ( -2.037 7.921 0.877) 8.226 8.899 ( -2.128 8.018 -0.675) 8.323 9.310 ( -6.652 -4.823 1.703) 8.391 9.348 ( -6.350 -4.784 -1.543) 8.099 9.944 ( -2.752 -3.411 1.327) 4.579 9.975 ( -2.683 -3.531 -1.350) 4.636 10.910 ( 3.615 -2.062 0.683) 4.217 10.926 ( 3.493 -2.326 -0.645) 4.246 ======================= Grid point 295 (46/60) ======================= q-point: ( 0.44 0.06 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.467 ( 2.566 1.395 1.249) 3.177 2.496 ( 2.409 1.035 -1.308) 2.930 2.521 ( -2.275 4.220 0.337) 4.806 2.527 ( -2.751 3.850 -0.226) 4.738 2.911 ( -4.945 6.300 1.468) 8.142 2.946 ( -4.829 4.633 -1.586) 6.877 3.431 ( -5.094 10.573 0.215) 11.738 3.436 ( -5.289 10.376 -0.196) 11.648 4.918 ( 2.834 1.178 1.597) 3.460 4.954 ( 2.777 1.067 -1.537) 3.348 6.418 ( 5.179 -0.657 1.393) 5.403 6.450 ( 5.021 -1.043 -1.390) 5.314 8.718 ( -1.186 3.399 0.595) 3.649 8.732 ( -1.069 3.261 -0.701) 3.503 8.914 ( -3.802 9.112 0.852) 9.910 8.929 ( 0.441 5.499 0.041) 5.517 8.940 ( 0.274 4.788 0.641) 4.839 8.975 ( 4.707 0.850 -1.950) 5.165 9.193 ( -1.258 -3.904 2.052) 4.586 9.234 ( -1.215 -3.976 -1.538) 4.432 9.870 ( -1.752 -1.863 1.323) 2.879 9.900 ( -1.792 -1.802 -1.335) 2.870 10.937 ( 2.246 -2.659 0.423) 3.506 10.946 ( 2.037 -2.763 -0.399) 3.455 ======================= Grid point 306 (47/60) ======================= q-point: ( 0.12 0.12 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 209 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.412 ( 9.678 16.763 0.945) 19.380 1.438 ( 8.491 14.707 -1.315) 17.033 2.106 ( 13.260 22.968 1.180) 26.547 2.134 ( 12.794 22.159 -1.272) 25.619 2.815 ( 6.434 11.144 1.242) 12.928 2.841 ( 7.244 12.547 -1.042) 14.525 2.979 ( 9.963 17.257 0.417) 19.931 2.988 ( 9.791 16.958 -0.326) 19.584 3.795 ( 11.011 19.072 1.513) 22.075 3.828 ( 11.055 19.147 -1.412) 22.154 5.379 ( 3.868 6.699 0.693) 7.767 5.396 ( 3.177 5.502 -0.749) 6.397 8.340 ( -2.445 -4.234 0.881) 4.968 8.360 ( -2.608 -4.517 -0.824) 5.280 8.648 ( 2.006 3.474 0.403) 4.032 8.659 ( 2.106 3.648 -0.526) 4.245 8.783 ( 2.172 3.763 0.629) 4.390 8.796 ( 2.220 3.846 -0.510) 4.470 9.872 ( 2.858 4.950 1.578) 5.930 9.919 ( 1.772 3.070 -2.800) 4.517 10.105 ( -2.586 -4.480 3.240) 6.104 10.157 ( -5.231 -9.060 -1.532) 10.573 10.468 ( 1.983 3.434 2.406) 4.638 10.532 ( 1.569 2.718 -3.149) 4.446 ======================= Grid point 307 (48/60) ======================= q-point: ( 0.19 0.12 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.738 ( 10.422 11.393 0.335) 15.444 1.748 ( 11.448 12.193 -0.556) 16.734 2.480 ( 3.179 20.771 0.767) 21.027 2.497 ( 3.017 20.947 -0.742) 21.177 2.974 ( 5.819 0.453 2.064) 6.191 3.021 ( 6.504 0.334 -2.028) 6.821 3.291 ( 4.438 18.460 0.211) 18.987 3.295 ( 4.525 18.311 -0.149) 18.863 4.156 ( 12.632 10.986 1.603) 16.818 4.192 ( 12.835 10.925 -1.521) 16.924 5.566 ( 8.240 7.968 0.327) 11.467 5.573 ( 9.088 8.625 -0.352) 12.534 8.306 ( 1.371 0.244 0.653) 1.538 8.321 ( 1.376 -0.089 -0.610) 1.508 8.722 ( 1.054 5.665 0.475) 5.782 8.735 ( 1.004 5.802 -0.637) 5.922 8.860 ( 0.977 5.721 0.683) 5.844 8.874 ( 0.957 5.808 -0.537) 5.910 9.859 ( -7.339 -5.244 0.965) 9.072 9.890 ( -5.926 -4.609 -1.979) 7.763 9.992 ( -1.296 -2.331 2.233) 3.478 10.029 ( -0.318 -2.269 -1.159) 2.568 10.589 ( 7.028 1.874 1.106) 7.357 10.616 ( 5.605 0.276 -1.219) 5.742 ======================= Grid point 308 (49/60) ======================= q-point: ( 0.25 0.12 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 386 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.023 ( 10.425 8.303 1.234) 13.384 2.052 ( 10.990 8.484 -1.353) 13.949 2.627 ( -3.331 9.841 1.256) 10.466 2.654 ( -3.604 11.671 -1.075) 12.262 3.062 ( 1.721 4.657 1.892) 5.313 3.108 ( 2.383 2.397 -2.071) 3.964 3.505 ( -1.009 18.165 0.320) 18.196 3.512 ( -0.859 18.227 -0.282) 18.249 4.457 ( 11.574 6.651 1.728) 13.460 4.496 ( 11.760 6.512 -1.633) 13.542 5.825 ( 11.738 7.187 1.093) 13.807 5.850 ( 12.371 7.332 -1.085) 14.421 8.394 ( 6.278 3.102 0.631) 7.031 8.408 ( 6.621 2.699 -0.600) 7.175 8.796 ( -0.297 7.377 0.446) 7.396 8.808 ( -0.395 7.491 -0.641) 7.529 8.930 ( -0.266 6.859 0.754) 6.905 8.945 ( -0.209 7.034 -0.592) 7.061 9.641 ( -10.416 -4.689 1.240) 11.489 9.672 ( -10.447 -4.774 -1.419) 11.573 9.981 ( 0.432 -2.952 1.188) 3.212 10.006 ( 0.372 -3.355 -0.993) 3.519 10.715 ( 6.707 -1.920 0.239) 6.981 10.721 ( 7.218 -1.841 -0.256) 7.453 ======================= Grid point 309 (50/60) ======================= q-point: ( 0.31 0.12 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.263 ( 8.391 6.137 1.430) 10.493 2.297 ( 8.461 5.755 -1.519) 10.345 2.612 ( -3.311 1.118 1.573) 3.833 2.646 ( -4.038 2.345 -1.392) 4.873 3.122 ( -4.392 11.349 1.265) 12.235 3.153 ( -4.037 9.173 -1.444) 10.126 3.629 ( -4.978 17.048 0.389) 17.764 3.638 ( -4.982 17.037 -0.363) 17.755 4.696 ( 9.201 4.443 1.719) 10.361 4.735 ( 9.205 4.233 -1.637) 10.263 6.094 ( 11.550 4.511 1.420) 12.481 6.127 ( 11.825 4.267 -1.413) 12.651 8.565 ( 8.811 2.883 0.698) 9.297 8.580 ( 8.942 2.572 -0.651) 9.327 8.859 ( -1.229 7.771 0.428) 7.879 8.872 ( -1.052 7.835 -0.755) 7.941 8.991 ( -0.522 6.657 0.992) 6.751 9.011 ( -0.434 7.087 -0.747) 7.140 9.423 ( -10.272 -1.782 1.372) 10.515 9.455 ( -10.083 -1.886 -1.402) 10.353 9.939 ( -1.074 -3.799 1.173) 4.118 9.966 ( -0.796 -3.861 -1.127) 4.100 10.796 ( 5.689 -3.590 0.521) 6.747 10.808 ( 5.884 -3.873 -0.499) 7.062 ======================= Grid point 310 (51/60) ======================= q-point: ( 0.38 0.12 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 386 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.430 ( 4.488 4.003 1.144) 6.121 2.458 ( 4.326 3.437 -1.235) 5.662 2.566 ( -1.391 -1.081 1.285) 2.181 2.594 ( -2.228 -0.621 -1.145) 2.581 3.141 ( -7.834 14.111 0.624) 16.151 3.156 ( -7.844 12.658 -0.739) 14.909 3.685 ( -7.343 15.828 0.291) 17.451 3.691 ( -7.513 15.703 -0.260) 17.410 4.862 ( 5.505 2.170 1.563) 6.120 4.897 ( 5.443 1.921 -1.512) 5.967 6.298 ( 7.912 0.961 1.298) 8.075 6.328 ( 7.810 0.473 -1.301) 7.931 8.725 ( 5.705 1.850 0.518) 6.020 8.735 ( 5.339 1.790 -0.413) 5.646 8.910 ( -0.532 5.252 0.686) 5.324 8.932 ( 0.468 4.782 -1.325) 4.984 9.050 ( 0.950 4.424 1.441) 4.749 9.077 ( 0.841 5.464 -0.984) 5.615 9.254 ( -9.582 3.053 1.402) 10.154 9.286 ( -9.349 2.302 -1.367) 9.725 9.879 ( -1.687 -2.616 1.347) 3.391 9.909 ( -1.544 -2.481 -1.323) 3.208 10.841 ( 4.252 -4.576 0.464) 6.264 10.852 ( 4.244 -4.723 -0.438) 6.364 ======================= Grid point 311 (52/60) ======================= q-point: ( 0.44 0.12 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 208 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.491 ( -0.667 1.155 0.907) 1.613 2.514 ( -0.477 0.826 -1.002) 1.383 2.549 ( 0.417 -0.722 0.972) 1.280 2.570 ( 0.078 -0.135 -0.852) 0.866 3.143 ( -8.690 15.052 0.145) 17.381 3.148 ( -8.073 13.983 -0.229) 16.147 3.700 ( -8.631 14.950 0.120) 17.263 3.703 ( -8.658 14.996 -0.086) 17.316 4.923 ( 0.447 -0.774 1.463) 1.715 4.956 ( 0.547 -0.948 -1.427) 1.798 6.375 ( 1.759 -3.046 1.071) 3.677 6.400 ( 1.935 -3.351 -1.076) 4.017 8.783 ( -0.860 1.489 0.363) 1.757 8.791 ( -0.429 0.743 -0.292) 0.907 8.936 ( -0.569 0.985 0.974) 1.497 8.966 ( -0.004 0.008 -1.766) 1.766 9.100 ( 2.623 -4.543 2.052) 5.633 9.137 ( 2.663 -4.612 -1.260) 5.472 9.167 ( -8.174 14.158 0.865) 16.371 9.189 ( -8.394 14.539 -0.970) 16.816 9.847 ( 0.216 -0.374 1.437) 1.500 9.880 ( 0.150 -0.260 -1.423) 1.454 10.855 ( 3.043 -5.270 0.432) 6.101 10.865 ( 3.006 -5.206 -0.407) 6.025 ======================= Grid point 324 (53/60) ======================= q-point: ( 0.19 0.19 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 208 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.970 ( 6.649 11.516 1.031) 13.337 1.995 ( 6.937 12.015 -1.150) 13.921 2.717 ( 0.545 0.944 2.067) 2.337 2.762 ( 1.425 2.468 -1.871) 3.409 3.071 ( 6.165 10.677 0.587) 12.343 3.086 ( 5.167 8.949 -0.766) 10.361 3.629 ( 8.181 14.170 0.075) 16.362 3.630 ( 8.195 14.194 -0.007) 16.389 4.375 ( 6.219 10.772 1.540) 12.533 4.410 ( 6.209 10.754 -1.493) 12.507 5.753 ( 6.018 10.424 0.856) 12.067 5.773 ( 6.290 10.894 -0.849) 12.608 8.354 ( 2.527 4.377 0.330) 5.065 8.361 ( 2.332 4.040 -0.298) 4.674 8.841 ( 3.312 5.737 0.515) 6.645 8.855 ( 3.318 5.747 -0.749) 6.678 8.968 ( 2.631 4.558 0.803) 5.324 8.984 ( 2.707 4.689 -0.592) 5.447 9.724 ( -4.289 -7.429 1.157) 8.656 9.755 ( -4.346 -7.528 -1.528) 8.826 9.959 ( -1.006 -1.743 1.243) 2.365 9.983 ( -1.312 -2.272 -0.849) 2.758 10.630 ( 1.275 2.209 0.125) 2.554 10.634 ( 0.881 1.526 -0.166) 1.769 ======================= Grid point 325 (54/60) ======================= q-point: ( 0.25 0.19 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.199 ( 7.316 9.161 1.232) 11.789 2.229 ( 7.365 8.951 -1.316) 11.666 2.653 ( -3.235 -3.917 2.266) 5.563 2.704 ( -3.689 -3.871 -2.203) 5.784 3.282 ( 2.044 13.352 0.341) 13.512 3.290 ( 2.529 12.558 -0.385) 12.816 3.852 ( 0.414 15.240 0.287) 15.248 3.858 ( 0.734 15.129 -0.240) 15.149 4.589 ( 8.068 6.467 1.548) 10.455 4.625 ( 8.173 6.281 -1.508) 10.418 5.976 ( 8.335 7.812 1.146) 11.481 6.002 ( 8.554 7.794 -1.139) 11.628 8.472 ( 5.389 4.547 0.280) 7.056 8.478 ( 5.873 4.221 -0.247) 7.237 8.942 ( 0.896 6.530 0.403) 6.604 8.956 ( 0.809 6.385 -0.796) 6.485 9.031 ( 0.448 2.264 1.102) 2.558 9.052 ( 0.656 2.510 -0.737) 2.697 9.584 ( -6.630 -0.228 1.205) 6.743 9.614 ( -6.639 -0.281 -1.397) 6.790 9.917 ( -0.821 -3.050 1.057) 3.331 9.939 ( -0.589 -2.908 -0.865) 3.091 10.672 ( 3.606 -2.033 0.180) 4.143 10.677 ( 3.899 -2.270 -0.185) 4.516 ======================= Grid point 326 (55/60) ======================= q-point: ( 0.31 0.19 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.398 ( 5.439 7.014 0.995) 8.932 2.422 ( 5.298 6.483 -1.106) 8.445 2.566 ( -1.757 -3.726 1.778) 4.487 2.606 ( -2.337 -4.134 -1.683) 5.039 3.398 ( -3.907 13.277 0.375) 13.845 3.407 ( -3.841 13.420 -0.405) 13.965 3.962 ( -4.918 14.917 0.350) 15.711 3.969 ( -4.959 14.901 -0.319) 15.708 4.776 ( 7.026 3.372 1.458) 7.929 4.809 ( 6.989 3.042 -1.425) 7.755 6.184 ( 7.638 4.439 1.137) 8.907 6.210 ( 7.629 4.070 -1.146) 8.722 8.620 ( 6.791 2.525 0.505) 7.263 8.631 ( 7.012 2.378 -0.451) 7.418 8.995 ( -1.423 4.492 -0.050) 4.712 8.998 ( -0.888 1.422 -0.333) 1.709 9.042 ( -0.384 1.265 1.017) 1.668 9.061 ( 0.358 -1.678 -0.678) 1.845 9.500 ( -8.015 8.873 1.147) 12.011 9.528 ( -8.084 8.784 -1.308) 12.009 9.875 ( -0.813 -2.155 1.344) 2.667 9.904 ( -0.570 -1.898 -1.194) 2.314 10.711 ( 4.466 -4.662 0.296) 6.463 10.717 ( 4.710 -4.782 -0.285) 6.718 ======================= Grid point 327 (56/60) ======================= q-point: ( 0.38 0.19 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.509 ( 0.870 1.270 0.642) 1.668 2.527 ( 0.732 1.211 0.503) 1.502 2.530 ( 0.719 1.670 -0.313) 1.845 2.551 ( 0.083 -1.466 -0.846) 1.695 3.432 ( -6.920 13.148 0.527) 14.867 3.445 ( -7.201 13.905 -0.559) 15.669 3.998 ( -7.386 14.192 0.281) 16.001 4.004 ( -7.577 14.366 -0.237) 16.244 4.888 ( 3.179 0.414 1.267) 3.447 4.917 ( 3.230 0.020 -1.250) 3.463 6.313 ( 3.590 0.724 0.862) 3.762 6.333 ( 3.509 0.259 -0.882) 3.627 8.737 ( 4.638 -0.632 0.372) 4.696 8.744 ( 4.217 -0.942 -0.266) 4.330 8.979 ( 0.900 -3.321 0.473) 3.473 8.993 ( -1.176 0.896 -0.826) 1.694 9.045 ( -1.368 2.133 0.213) 2.544 9.048 ( 0.639 -2.190 0.053) 2.282 9.472 ( -8.661 14.001 1.072) 16.498 9.499 ( -8.830 14.172 -1.226) 16.743 9.850 ( -0.658 -0.205 1.537) 1.684 9.883 ( -0.573 -0.123 -1.409) 1.526 10.732 ( 4.092 -6.038 0.424) 7.306 10.741 ( 4.208 -6.063 -0.400) 7.391 ======================= Grid point 341 (57/60) ======================= q-point: ( 0.25 0.25 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 209 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.385 ( 5.052 8.751 0.920) 10.147 2.408 ( 4.812 8.335 -1.038) 9.681 2.562 ( -2.404 -4.164 1.891) 5.166 2.604 ( -2.818 -4.881 -1.803) 5.917 3.505 ( 4.682 8.110 0.448) 9.375 3.515 ( 5.014 8.685 -0.449) 10.039 4.096 ( 4.938 8.553 0.338) 9.882 4.103 ( 5.074 8.788 -0.304) 10.152 4.720 ( 3.660 6.339 1.372) 7.447 4.751 ( 3.529 6.113 -1.353) 7.187 6.137 ( 4.484 7.767 1.046) 9.029 6.161 ( 4.369 7.568 -1.058) 8.802 8.570 ( 2.869 4.968 0.316) 5.746 8.577 ( 3.017 5.226 -0.250) 6.040 8.994 ( -1.749 -3.030 0.811) 3.591 9.018 ( -0.840 -1.455 -1.312) 2.132 9.064 ( 0.922 1.598 0.641) 1.953 9.073 ( 1.710 2.961 -0.193) 3.425 9.630 ( 2.217 3.839 1.113) 4.571 9.659 ( 2.220 3.845 -1.417) 4.661 9.870 ( -0.790 -1.368 1.388) 2.103 9.899 ( -0.607 -1.051 -1.082) 1.626 10.637 ( -0.870 -1.506 0.107) 1.743 10.639 ( -0.740 -1.281 -0.107) 1.483 ======================= Grid point 342 (58/60) ======================= q-point: ( 0.31 0.25 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 384 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.494 ( -0.531 -2.076 0.981) 2.357 2.520 ( -0.723 -2.339 -1.307) 2.775 2.540 ( 2.644 5.601 0.776) 6.242 2.555 ( 2.215 4.632 -0.457) 5.154 3.614 ( -1.009 7.831 0.926) 7.950 3.636 ( -0.887 8.544 -0.954) 8.642 4.207 ( -1.812 8.949 0.504) 9.144 4.218 ( -1.866 9.254 -0.463) 9.452 4.834 ( 3.759 2.499 1.080) 4.641 4.859 ( 3.604 1.963 -1.070) 4.241 6.273 ( 3.609 4.190 0.762) 5.582 6.291 ( 3.402 3.719 -0.791) 5.102 8.665 ( 3.524 1.844 0.456) 4.004 8.673 ( 3.437 1.711 -0.302) 3.851 8.923 ( -0.765 -4.703 1.059) 4.881 8.952 ( -0.802 -5.084 -1.478) 5.355 9.090 ( -0.939 2.249 0.801) 2.565 9.104 ( -0.929 2.217 -0.515) 2.459 9.703 ( -2.291 9.592 1.099) 9.923 9.733 ( -2.263 9.668 -1.590) 10.056 9.861 ( -0.238 0.805 1.647) 1.848 9.892 ( -0.099 0.563 -1.161) 1.294 10.615 ( 1.877 -4.492 0.258) 4.875 10.621 ( 2.021 -4.574 -0.249) 5.007 ======================= Grid point 343 (59/60) ======================= q-point: ( 0.38 0.25 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 209 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.477 ( 0.803 -1.390 0.696) 1.749 2.494 ( 1.090 -1.888 -0.789) 2.318 2.597 ( -1.768 3.062 0.413) 3.560 2.606 ( -1.586 2.746 -0.318) 3.187 3.643 ( -4.153 7.193 1.134) 8.383 3.669 ( -4.470 7.743 -1.182) 9.019 4.233 ( -4.920 8.522 0.557) 9.856 4.245 ( -5.113 8.855 -0.495) 10.237 4.883 ( 0.354 -0.613 0.841) 1.099 4.903 ( 0.671 -1.162 -0.840) 1.583 6.329 ( -0.481 0.833 0.488) 1.078 6.341 ( -0.274 0.474 -0.522) 0.756 8.712 ( 0.935 -1.620 0.127) 1.875 8.712 ( 0.841 -1.456 0.102) 1.685 8.891 ( 2.067 -3.581 0.886) 4.228 8.916 ( 2.451 -4.245 -1.257) 5.060 9.089 ( -1.475 2.555 0.676) 3.027 9.103 ( -1.302 2.255 -0.508) 2.653 9.732 ( -5.991 10.377 1.093) 12.032 9.763 ( -5.905 10.228 -1.681) 11.929 9.865 ( -1.058 1.832 1.680) 2.701 9.896 ( -0.657 1.138 -1.105) 1.717 10.609 ( 3.194 -5.532 0.286) 6.394 10.616 ( 3.360 -5.820 -0.274) 6.726 ======================= Grid point 359 (60/60) ======================= q-point: ( 0.31 0.31 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.97e-04 1.97e-04 1.97e-04 1.97e-04 Number of triplets: 208 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.463 ( -0.343 -0.594 0.421) 0.805 2.473 ( -0.803 -1.390 -0.481) 1.676 2.631 ( 1.488 2.578 0.243) 2.986 2.635 ( 1.268 2.196 -0.170) 2.541 3.716 ( 1.376 2.384 1.379) 3.079 3.748 ( 1.550 2.685 -1.446) 3.421 4.323 ( 1.566 2.712 0.734) 3.216 4.339 ( 1.641 2.842 -0.662) 3.348 4.875 ( 0.776 1.344 0.474) 1.623 4.886 ( 0.331 0.574 -0.476) 0.816 6.339 ( 1.150 1.992 0.294) 2.319 6.347 ( 0.833 1.443 -0.334) 1.699 8.692 ( 0.487 0.844 0.344) 1.034 8.696 ( 0.295 0.510 -0.005) 0.589 8.854 ( -0.959 -1.661 0.421) 1.964 8.868 ( -1.505 -2.607 -0.769) 3.106 9.118 ( 0.406 0.704 0.367) 0.892 9.126 ( 0.099 0.171 -0.333) 0.388 9.841 ( 1.902 3.295 0.263) 3.813 9.847 ( 0.880 1.524 -0.264) 1.780 9.917 ( 0.861 1.492 0.074) 1.724 9.919 ( 2.111 3.657 -0.072) 4.223 10.543 ( -1.242 -2.151 0.090) 2.486 10.545 ( -1.371 -2.375 -0.086) 2.743 =================== 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/18432 10.0 413.709 413.709 18.255 0.000 0.000 0.000 3/18432 20.0 469.101 469.101 11.769 -0.000 0.000 0.000 3/18432 30.0 392.956 392.956 8.024 -0.000 0.000 0.000 3/18432 40.0 325.343 325.343 6.012 -0.000 0.000 0.000 3/18432 50.0 273.083 273.083 4.803 -0.000 0.000 0.000 3/18432 60.0 231.351 231.351 3.980 0.000 0.000 0.000 3/18432 70.0 197.610 197.610 3.379 0.000 0.000 0.000 3/18432 80.0 170.527 170.527 2.926 0.000 0.000 0.000 3/18432 90.0 148.917 148.917 2.576 0.000 0.000 0.000 3/18432 100.0 131.640 131.640 2.301 0.000 0.000 0.000 3/18432 110.0 117.715 117.715 2.080 0.000 0.000 0.000 3/18432 120.0 106.360 106.360 1.900 0.000 0.000 0.000 3/18432 130.0 96.982 96.982 1.750 0.000 0.000 0.000 3/18432 140.0 89.135 89.135 1.623 0.000 0.000 0.000 3/18432 150.0 82.489 82.489 1.515 0.000 0.000 0.000 3/18432 160.0 76.795 76.795 1.421 0.000 0.000 0.000 3/18432 170.0 71.866 71.866 1.339 0.000 0.000 0.000 3/18432 180.0 67.559 67.559 1.266 0.000 0.000 0.000 3/18432 190.0 63.763 63.763 1.201 0.000 0.000 0.000 3/18432 200.0 60.392 60.392 1.143 0.000 0.000 0.000 3/18432 210.0 57.377 57.377 1.091 0.000 0.000 0.000 3/18432 220.0 54.664 54.664 1.043 0.000 0.000 0.000 3/18432 230.0 52.209 52.209 0.999 0.000 0.000 0.000 3/18432 240.0 49.977 49.977 0.959 0.000 0.000 0.000 3/18432 250.0 47.936 47.936 0.922 0.000 0.000 0.000 3/18432 260.0 46.064 46.064 0.888 0.000 0.000 0.000 3/18432 270.0 44.340 44.340 0.857 0.000 0.000 0.000 3/18432 280.0 42.746 42.746 0.827 0.000 0.000 0.000 3/18432 290.0 41.267 41.267 0.800 0.000 0.000 0.000 3/18432 300.0 39.892 39.892 0.774 0.000 0.000 0.000 3/18432 310.0 38.609 38.609 0.750 0.000 0.000 0.000 3/18432 320.0 37.410 37.410 0.728 0.000 0.000 0.000 3/18432 330.0 36.285 36.285 0.707 0.000 0.000 0.000 3/18432 340.0 35.229 35.229 0.687 0.000 0.000 0.000 3/18432 350.0 34.235 34.235 0.668 0.000 0.000 0.000 3/18432 360.0 33.297 33.297 0.650 0.000 0.000 0.000 3/18432 370.0 32.411 32.411 0.633 0.000 0.000 0.000 3/18432 380.0 31.572 31.572 0.617 0.000 0.000 0.000 3/18432 390.0 30.777 30.777 0.602 0.000 0.000 0.000 3/18432 400.0 30.023 30.023 0.587 0.000 0.000 0.000 3/18432 410.0 29.305 29.305 0.573 0.000 0.000 0.000 3/18432 420.0 28.622 28.622 0.560 0.000 0.000 0.000 3/18432 430.0 27.971 27.971 0.548 0.000 0.000 0.000 3/18432 440.0 27.349 27.349 0.536 0.000 0.000 0.000 3/18432 450.0 26.756 26.756 0.524 0.000 0.000 0.000 3/18432 460.0 26.188 26.188 0.513 0.000 0.000 0.000 3/18432 470.0 25.644 25.644 0.502 0.000 0.000 0.000 3/18432 480.0 25.123 25.123 0.492 0.000 0.000 0.000 3/18432 490.0 24.623 24.623 0.483 0.000 0.000 0.000 3/18432 500.0 24.143 24.143 0.473 0.000 0.000 0.000 3/18432 510.0 23.682 23.682 0.464 0.000 0.000 0.000 3/18432 520.0 23.238 23.238 0.456 0.000 0.000 0.000 3/18432 530.0 22.811 22.811 0.447 0.000 0.000 0.000 3/18432 540.0 22.400 22.400 0.439 0.000 0.000 0.000 3/18432 550.0 22.004 22.004 0.431 0.000 0.000 0.000 3/18432 560.0 21.621 21.621 0.424 0.000 0.000 0.000 3/18432 570.0 21.252 21.252 0.417 0.000 0.000 0.000 3/18432 580.0 20.896 20.896 0.410 0.000 0.000 0.000 3/18432 590.0 20.551 20.551 0.403 0.000 0.000 0.000 3/18432 600.0 20.218 20.218 0.396 0.000 0.000 0.000 3/18432 610.0 19.895 19.895 0.390 0.000 0.000 0.000 3/18432 620.0 19.583 19.583 0.384 0.000 0.000 0.000 3/18432 630.0 19.281 19.281 0.378 0.000 0.000 0.000 3/18432 640.0 18.988 18.988 0.372 0.000 0.000 0.000 3/18432 650.0 18.704 18.704 0.367 0.000 0.000 0.000 3/18432 660.0 18.428 18.428 0.361 0.000 0.000 0.000 3/18432 670.0 18.160 18.160 0.356 0.000 0.000 0.000 3/18432 680.0 17.900 17.900 0.351 0.000 0.000 0.000 3/18432 690.0 17.648 17.648 0.346 0.000 0.000 0.000 3/18432 700.0 17.403 17.403 0.341 0.000 0.000 0.000 3/18432 710.0 17.164 17.164 0.336 0.000 0.000 0.000 3/18432 720.0 16.932 16.932 0.332 0.000 0.000 0.000 3/18432 730.0 16.707 16.707 0.327 0.000 0.000 0.000 3/18432 740.0 16.487 16.487 0.323 0.000 0.000 0.000 3/18432 750.0 16.273 16.273 0.319 0.000 0.000 0.000 3/18432 760.0 16.064 16.064 0.315 0.000 0.000 0.000 3/18432 770.0 15.861 15.861 0.311 0.000 0.000 0.000 3/18432 780.0 15.663 15.663 0.307 0.000 0.000 0.000 3/18432 790.0 15.470 15.470 0.303 0.000 0.000 0.000 3/18432 800.0 15.282 15.282 0.299 0.000 0.000 0.000 3/18432 810.0 15.098 15.098 0.296 0.000 0.000 0.000 3/18432 820.0 14.919 14.919 0.292 0.000 0.000 0.000 3/18432 830.0 14.744 14.744 0.289 0.000 0.000 0.000 3/18432 840.0 14.573 14.573 0.285 0.000 0.000 0.000 3/18432 850.0 14.406 14.406 0.282 0.000 0.000 0.000 3/18432 860.0 14.243 14.243 0.279 0.000 0.000 0.000 3/18432 870.0 14.083 14.083 0.276 0.000 0.000 0.000 3/18432 880.0 13.927 13.927 0.272 0.000 0.000 0.000 3/18432 890.0 13.775 13.775 0.269 0.000 0.000 0.000 3/18432 900.0 13.625 13.625 0.267 0.000 0.000 0.000 3/18432 910.0 13.479 13.479 0.264 0.000 0.000 0.000 3/18432 920.0 13.336 13.336 0.261 0.000 0.000 0.000 3/18432 930.0 13.197 13.197 0.258 0.000 0.000 0.000 3/18432 940.0 13.060 13.060 0.255 0.000 0.000 0.000 3/18432 950.0 12.926 12.926 0.253 0.000 0.000 0.000 3/18432 960.0 12.794 12.794 0.250 0.000 0.000 0.000 3/18432 970.0 12.666 12.666 0.248 0.000 0.000 0.000 3/18432 980.0 12.539 12.539 0.245 0.000 0.000 0.000 3/18432 990.0 12.416 12.416 0.243 0.000 0.000 0.000 3/18432 1000.0 12.295 12.295 0.240 0.000 0.000 0.000 3/18432 Thermal conductivity related properties were written into "kappa-m16163.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 19:11:26]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|