----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 08:17:04]------------------------- Compiled with OpenMP support (max 48 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) Number of symmetry operations in supercell: 192 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.871573573743976 0.000000000000000 0.000000000000000 b -1.935786786871988 3.352881067482788 0.000000000000000 c 0.000000000000000 0.000000000000000 16.139780439999999 Atomic positions (fractional): *1 Mg 0.33333333333333 0.66666666666667 0.40180910014792 24.305 2 Mg 0.66666666666667 0.33333333333333 0.59819089985208 24.305 *3 Al 0.33333333333333 0.66666666666667 0.82775598530114 26.982 4 Al 0.66666666666667 0.33333333333333 0.17224401469886 26.982 *5 Se 0.66666666666667 0.33333333333333 0.89425500396964 78.960 6 Se 0.33333333333333 0.66666666666667 0.10574499603036 78.960 *7 Se 0.66666666666667 0.33333333333333 0.31447675906161 78.960 8 Se 0.33333333333333 0.66666666666667 0.68552324093839 78.960 *9 Se 0.00000000000000 0.00000000000000 0.50000000000000 78.960 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.871573573743976 0.000000000000000 0.000000000000000 b -1.935786786871988 3.352881067482788 0.000000000000000 c 0.000000000000000 0.000000000000000 16.139780439999999 Atomic positions (fractional): *1 Mg 0.33333333333333 0.66666666666667 0.40180910014792 24.305 > 1 2 Mg 0.66666666666667 0.33333333333333 0.59819089985208 24.305 > 2 *3 Al 0.33333333333333 0.66666666666667 0.82775598530114 26.982 > 3 4 Al 0.66666666666667 0.33333333333333 0.17224401469886 26.982 > 4 *5 Se 0.66666666666667 0.33333333333333 0.89425500396964 78.960 > 5 6 Se 0.33333333333333 0.66666666666667 0.10574499603036 78.960 > 6 *7 Se 0.66666666666667 0.33333333333333 0.31447675906161 78.960 > 7 8 Se 0.33333333333333 0.66666666666667 0.68552324093839 78.960 > 8 *9 Se 0.00000000000000 0.00000000000000 0.50000000000000 78.960 > 9 -------------------------------- super cell -------------------------------- Lattice vectors: a 15.486294294975904 0.000000000000000 0.000000000000000 b -7.743147147487952 13.411524269931151 0.000000000000000 c 0.000000000000000 0.000000000000000 16.139780439999999 Atomic positions (fractional): *1 Mg 0.08333333333333 0.16666666666667 0.40180910014792 24.305 > 1 2 Mg 0.33333333333333 0.16666666666667 0.40180910014792 24.305 > 1 3 Mg 0.58333333333333 0.16666666666667 0.40180910014792 24.305 > 1 4 Mg 0.83333333333333 0.16666666666667 0.40180910014792 24.305 > 1 5 Mg 0.08333333333333 0.41666666666667 0.40180910014792 24.305 > 1 6 Mg 0.33333333333333 0.41666666666667 0.40180910014792 24.305 > 1 7 Mg 0.58333333333333 0.41666666666667 0.40180910014792 24.305 > 1 8 Mg 0.83333333333333 0.41666666666667 0.40180910014792 24.305 > 1 9 Mg 0.08333333333333 0.66666666666667 0.40180910014792 24.305 > 1 10 Mg 0.33333333333333 0.66666666666667 0.40180910014792 24.305 > 1 11 Mg 0.58333333333333 0.66666666666667 0.40180910014792 24.305 > 1 12 Mg 0.83333333333333 0.66666666666667 0.40180910014792 24.305 > 1 13 Mg 0.08333333333333 0.91666666666667 0.40180910014792 24.305 > 1 14 Mg 0.33333333333333 0.91666666666667 0.40180910014792 24.305 > 1 15 Mg 0.58333333333333 0.91666666666667 0.40180910014792 24.305 > 1 16 Mg 0.83333333333333 0.91666666666667 0.40180910014792 24.305 > 1 17 Mg 0.16666666666667 0.08333333333333 0.59819089985208 24.305 > 2 18 Mg 0.41666666666667 0.08333333333333 0.59819089985208 24.305 > 2 19 Mg 0.66666666666667 0.08333333333333 0.59819089985208 24.305 > 2 20 Mg 0.91666666666667 0.08333333333333 0.59819089985208 24.305 > 2 21 Mg 0.16666666666667 0.33333333333333 0.59819089985208 24.305 > 2 22 Mg 0.41666666666667 0.33333333333333 0.59819089985208 24.305 > 2 23 Mg 0.66666666666667 0.33333333333333 0.59819089985208 24.305 > 2 24 Mg 0.91666666666667 0.33333333333333 0.59819089985208 24.305 > 2 25 Mg 0.16666666666667 0.58333333333333 0.59819089985208 24.305 > 2 26 Mg 0.41666666666667 0.58333333333333 0.59819089985208 24.305 > 2 27 Mg 0.66666666666667 0.58333333333333 0.59819089985208 24.305 > 2 28 Mg 0.91666666666667 0.58333333333333 0.59819089985208 24.305 > 2 29 Mg 0.16666666666667 0.83333333333333 0.59819089985208 24.305 > 2 30 Mg 0.41666666666667 0.83333333333333 0.59819089985208 24.305 > 2 31 Mg 0.66666666666667 0.83333333333333 0.59819089985208 24.305 > 2 32 Mg 0.91666666666667 0.83333333333333 0.59819089985208 24.305 > 2 *33 Al 0.08333333333333 0.16666666666667 0.82775598530114 26.982 > 3 34 Al 0.33333333333333 0.16666666666667 0.82775598530114 26.982 > 3 35 Al 0.58333333333333 0.16666666666667 0.82775598530114 26.982 > 3 36 Al 0.83333333333333 0.16666666666667 0.82775598530114 26.982 > 3 37 Al 0.08333333333333 0.41666666666667 0.82775598530114 26.982 > 3 38 Al 0.33333333333333 0.41666666666667 0.82775598530114 26.982 > 3 39 Al 0.58333333333333 0.41666666666667 0.82775598530114 26.982 > 3 40 Al 0.83333333333333 0.41666666666667 0.82775598530114 26.982 > 3 41 Al 0.08333333333333 0.66666666666667 0.82775598530114 26.982 > 3 42 Al 0.33333333333333 0.66666666666667 0.82775598530114 26.982 > 3 43 Al 0.58333333333333 0.66666666666667 0.82775598530114 26.982 > 3 44 Al 0.83333333333333 0.66666666666667 0.82775598530114 26.982 > 3 45 Al 0.08333333333333 0.91666666666667 0.82775598530114 26.982 > 3 46 Al 0.33333333333333 0.91666666666667 0.82775598530114 26.982 > 3 47 Al 0.58333333333333 0.91666666666667 0.82775598530114 26.982 > 3 48 Al 0.83333333333333 0.91666666666667 0.82775598530114 26.982 > 3 49 Al 0.16666666666667 0.08333333333333 0.17224401469886 26.982 > 4 50 Al 0.41666666666667 0.08333333333333 0.17224401469886 26.982 > 4 51 Al 0.66666666666667 0.08333333333333 0.17224401469886 26.982 > 4 52 Al 0.91666666666667 0.08333333333333 0.17224401469886 26.982 > 4 53 Al 0.16666666666667 0.33333333333333 0.17224401469886 26.982 > 4 54 Al 0.41666666666667 0.33333333333333 0.17224401469886 26.982 > 4 55 Al 0.66666666666667 0.33333333333333 0.17224401469886 26.982 > 4 56 Al 0.91666666666667 0.33333333333333 0.17224401469886 26.982 > 4 57 Al 0.16666666666667 0.58333333333333 0.17224401469886 26.982 > 4 58 Al 0.41666666666667 0.58333333333333 0.17224401469886 26.982 > 4 59 Al 0.66666666666667 0.58333333333333 0.17224401469886 26.982 > 4 60 Al 0.91666666666667 0.58333333333333 0.17224401469886 26.982 > 4 61 Al 0.16666666666667 0.83333333333333 0.17224401469886 26.982 > 4 62 Al 0.41666666666667 0.83333333333333 0.17224401469886 26.982 > 4 63 Al 0.66666666666667 0.83333333333333 0.17224401469886 26.982 > 4 64 Al 0.91666666666667 0.83333333333333 0.17224401469886 26.982 > 4 *65 Se 0.16666666666667 0.08333333333333 0.89425500396964 78.960 > 5 66 Se 0.41666666666667 0.08333333333333 0.89425500396964 78.960 > 5 67 Se 0.66666666666667 0.08333333333333 0.89425500396964 78.960 > 5 68 Se 0.91666666666667 0.08333333333333 0.89425500396964 78.960 > 5 69 Se 0.16666666666667 0.33333333333333 0.89425500396964 78.960 > 5 70 Se 0.41666666666667 0.33333333333333 0.89425500396964 78.960 > 5 71 Se 0.66666666666667 0.33333333333333 0.89425500396964 78.960 > 5 72 Se 0.91666666666667 0.33333333333333 0.89425500396964 78.960 > 5 73 Se 0.16666666666667 0.58333333333333 0.89425500396964 78.960 > 5 74 Se 0.41666666666667 0.58333333333333 0.89425500396964 78.960 > 5 75 Se 0.66666666666667 0.58333333333333 0.89425500396964 78.960 > 5 76 Se 0.91666666666667 0.58333333333333 0.89425500396964 78.960 > 5 77 Se 0.16666666666667 0.83333333333333 0.89425500396964 78.960 > 5 78 Se 0.41666666666667 0.83333333333333 0.89425500396964 78.960 > 5 79 Se 0.66666666666667 0.83333333333333 0.89425500396964 78.960 > 5 80 Se 0.91666666666667 0.83333333333333 0.89425500396964 78.960 > 5 81 Se 0.08333333333333 0.16666666666667 0.10574499603036 78.960 > 6 82 Se 0.33333333333333 0.16666666666667 0.10574499603036 78.960 > 6 83 Se 0.58333333333333 0.16666666666667 0.10574499603036 78.960 > 6 84 Se 0.83333333333333 0.16666666666667 0.10574499603036 78.960 > 6 85 Se 0.08333333333333 0.41666666666667 0.10574499603036 78.960 > 6 86 Se 0.33333333333333 0.41666666666667 0.10574499603036 78.960 > 6 87 Se 0.58333333333333 0.41666666666667 0.10574499603036 78.960 > 6 88 Se 0.83333333333333 0.41666666666667 0.10574499603036 78.960 > 6 89 Se 0.08333333333333 0.66666666666667 0.10574499603036 78.960 > 6 90 Se 0.33333333333333 0.66666666666667 0.10574499603036 78.960 > 6 91 Se 0.58333333333333 0.66666666666667 0.10574499603036 78.960 > 6 92 Se 0.83333333333333 0.66666666666667 0.10574499603036 78.960 > 6 93 Se 0.08333333333333 0.91666666666667 0.10574499603036 78.960 > 6 94 Se 0.33333333333333 0.91666666666667 0.10574499603036 78.960 > 6 95 Se 0.58333333333333 0.91666666666667 0.10574499603036 78.960 > 6 96 Se 0.83333333333333 0.91666666666667 0.10574499603036 78.960 > 6 *97 Se 0.16666666666667 0.08333333333333 0.31447675906161 78.960 > 7 98 Se 0.41666666666667 0.08333333333333 0.31447675906161 78.960 > 7 99 Se 0.66666666666667 0.08333333333333 0.31447675906161 78.960 > 7 100 Se 0.91666666666667 0.08333333333333 0.31447675906161 78.960 > 7 101 Se 0.16666666666667 0.33333333333333 0.31447675906161 78.960 > 7 102 Se 0.41666666666667 0.33333333333333 0.31447675906161 78.960 > 7 103 Se 0.66666666666667 0.33333333333333 0.31447675906161 78.960 > 7 104 Se 0.91666666666667 0.33333333333333 0.31447675906161 78.960 > 7 105 Se 0.16666666666667 0.58333333333333 0.31447675906161 78.960 > 7 106 Se 0.41666666666667 0.58333333333333 0.31447675906161 78.960 > 7 107 Se 0.66666666666667 0.58333333333333 0.31447675906161 78.960 > 7 108 Se 0.91666666666667 0.58333333333333 0.31447675906161 78.960 > 7 109 Se 0.16666666666667 0.83333333333333 0.31447675906161 78.960 > 7 110 Se 0.41666666666667 0.83333333333333 0.31447675906161 78.960 > 7 111 Se 0.66666666666667 0.83333333333333 0.31447675906161 78.960 > 7 112 Se 0.91666666666667 0.83333333333333 0.31447675906161 78.960 > 7 113 Se 0.08333333333333 0.16666666666667 0.68552324093839 78.960 > 8 114 Se 0.33333333333333 0.16666666666667 0.68552324093839 78.960 > 8 115 Se 0.58333333333333 0.16666666666667 0.68552324093839 78.960 > 8 116 Se 0.83333333333333 0.16666666666667 0.68552324093839 78.960 > 8 117 Se 0.08333333333333 0.41666666666667 0.68552324093839 78.960 > 8 118 Se 0.33333333333333 0.41666666666667 0.68552324093839 78.960 > 8 119 Se 0.58333333333333 0.41666666666667 0.68552324093839 78.960 > 8 120 Se 0.83333333333333 0.41666666666667 0.68552324093839 78.960 > 8 121 Se 0.08333333333333 0.66666666666667 0.68552324093839 78.960 > 8 122 Se 0.33333333333333 0.66666666666667 0.68552324093839 78.960 > 8 123 Se 0.58333333333333 0.66666666666667 0.68552324093839 78.960 > 8 124 Se 0.83333333333333 0.66666666666667 0.68552324093839 78.960 > 8 125 Se 0.08333333333333 0.91666666666667 0.68552324093839 78.960 > 8 126 Se 0.33333333333333 0.91666666666667 0.68552324093839 78.960 > 8 127 Se 0.58333333333333 0.91666666666667 0.68552324093839 78.960 > 8 128 Se 0.83333333333333 0.91666666666667 0.68552324093839 78.960 > 8 *129 Se 0.00000000000000 0.00000000000000 0.50000000000000 78.960 > 9 130 Se 0.25000000000000 0.00000000000000 0.50000000000000 78.960 > 9 131 Se 0.50000000000000 0.00000000000000 0.50000000000000 78.960 > 9 132 Se 0.75000000000000 0.00000000000000 0.50000000000000 78.960 > 9 133 Se 0.00000000000000 0.25000000000000 0.50000000000000 78.960 > 9 134 Se 0.25000000000000 0.25000000000000 0.50000000000000 78.960 > 9 135 Se 0.50000000000000 0.25000000000000 0.50000000000000 78.960 > 9 136 Se 0.75000000000000 0.25000000000000 0.50000000000000 78.960 > 9 137 Se 0.00000000000000 0.50000000000000 0.50000000000000 78.960 > 9 138 Se 0.25000000000000 0.50000000000000 0.50000000000000 78.960 > 9 139 Se 0.50000000000000 0.50000000000000 0.50000000000000 78.960 > 9 140 Se 0.75000000000000 0.50000000000000 0.50000000000000 78.960 > 9 141 Se 0.00000000000000 0.75000000000000 0.50000000000000 78.960 > 9 142 Se 0.25000000000000 0.75000000000000 0.50000000000000 78.960 > 9 143 Se 0.50000000000000 0.75000000000000 0.50000000000000 78.960 > 9 144 Se 0.75000000000000 0.75000000000000 0.50000000000000 78.960 > 9 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 6.8093533 0.0000000 0.0000000 0.0000000 6.8093533 0.0000000 0.0000000 0.0000000 6.0776553 -------------------------- Born effective charges -------------------------- 1 Mg 2.4325488 -0.0000000 0.0000000 0.0000000 2.4325488 0.0000000 0.0000000 0.0000000 1.9913658 2 Mg 2.4325488 -0.0000000 0.0000000 0.0000000 2.4325488 0.0000000 0.0000000 0.0000000 1.9913658 3 Al 2.5785105 -0.0000000 0.0000000 0.0000000 2.5785105 0.0000000 0.0000000 0.0000000 2.0509402 4 Al 2.5785105 -0.0000000 0.0000000 0.0000000 2.5785105 0.0000000 0.0000000 0.0000000 2.0509402 5 Se -2.2486373 0.0000000 0.0000000 0.0000000 -2.2486373 0.0000000 0.0000000 0.0000000 -0.8865659 6 Se -2.2486373 0.0000000 0.0000000 0.0000000 -2.2486373 0.0000000 0.0000000 0.0000000 -0.8865659 7 Se -1.4550392 -0.0000000 0.0000000 0.0000000 -1.4550392 0.0000000 0.0000000 0.0000000 -2.1607832 8 Se -1.4550392 -0.0000000 0.0000000 0.0000000 -1.4550392 0.0000000 0.0000000 0.0000000 -2.1607832 9 Se -2.6147657 0.0000000 0.0000000 0.0000000 -2.6147657 0.0000000 0.0000000 0.0000000 -1.9899137 ---------------------------------------------------------------------------- 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: 432/432 Permutation basis: 6426/6426 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 578 Number of blocks in projector: 485 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: 255 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 203 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 120 Use standard eigh solver. Tree of FC basis block matrices: - (578, 568), data: False |-- (120, 115), data: True |-- (203, 203), data: True |-- (255, 250), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.012 Solver_block: 100 / 160 - Time: 0.451 Solver_block: 160 / 160 - Time: 0.307 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.774 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) Permutation basis: 432/432 Permutation basis: 6426/6426 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 578 Number of blocks in projector: 485 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: 255 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 203 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 120 Use standard eigh solver. Tree of FC basis block matrices: - (578, 568), data: False |-- (120, 115), data: True |-- (203, 203), data: True |-- (255, 250), 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-08 08:17:12]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 08:17:12]------------------------- Compiled with OpenMP support (max 48 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.871573573743976 0.000000000000000 0.000000000000000 b -1.935786786871988 3.352881067482788 0.000000000000000 c 0.000000000000000 0.000000000000000 16.139780439999999 Atomic positions (fractional): 1 Mg 0.33333333333333 0.66666666666667 0.40180910014792 24.305 2 Mg 0.66666666666667 0.33333333333333 0.59819089985208 24.305 3 Al 0.33333333333333 0.66666666666667 0.82775598530114 26.982 4 Al 0.66666666666667 0.33333333333333 0.17224401469886 26.982 5 Se 0.66666666666667 0.33333333333333 0.89425500396964 78.960 6 Se 0.33333333333333 0.66666666666667 0.10574499603036 78.960 7 Se 0.66666666666667 0.33333333333333 0.31447675906161 78.960 8 Se 0.33333333333333 0.66666666666667 0.68552324093839 78.960 9 Se 0.00000000000000 0.00000000000000 0.50000000000000 78.960 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.486294294975904 0.000000000000000 0.000000000000000 b -7.743147147487952 13.411524269931151 0.000000000000000 c 0.000000000000000 0.000000000000000 16.139780439999999 Atomic positions (fractional): 1 Mg 0.08333333333333 0.16666666666667 0.40180910014792 24.305 > 1 2 Mg 0.33333333333333 0.16666666666667 0.40180910014792 24.305 > 1 3 Mg 0.58333333333333 0.16666666666667 0.40180910014792 24.305 > 1 4 Mg 0.83333333333333 0.16666666666667 0.40180910014792 24.305 > 1 5 Mg 0.08333333333333 0.41666666666667 0.40180910014792 24.305 > 1 6 Mg 0.33333333333333 0.41666666666667 0.40180910014792 24.305 > 1 7 Mg 0.58333333333333 0.41666666666667 0.40180910014792 24.305 > 1 8 Mg 0.83333333333333 0.41666666666667 0.40180910014792 24.305 > 1 9 Mg 0.08333333333333 0.66666666666667 0.40180910014792 24.305 > 1 10 Mg 0.33333333333333 0.66666666666667 0.40180910014792 24.305 > 1 11 Mg 0.58333333333333 0.66666666666667 0.40180910014792 24.305 > 1 12 Mg 0.83333333333333 0.66666666666667 0.40180910014792 24.305 > 1 13 Mg 0.08333333333333 0.91666666666667 0.40180910014792 24.305 > 1 14 Mg 0.33333333333333 0.91666666666667 0.40180910014792 24.305 > 1 15 Mg 0.58333333333333 0.91666666666667 0.40180910014792 24.305 > 1 16 Mg 0.83333333333333 0.91666666666667 0.40180910014792 24.305 > 1 17 Mg 0.16666666666667 0.08333333333333 0.59819089985208 24.305 > 17 18 Mg 0.41666666666667 0.08333333333333 0.59819089985208 24.305 > 17 19 Mg 0.66666666666667 0.08333333333333 0.59819089985208 24.305 > 17 20 Mg 0.91666666666667 0.08333333333333 0.59819089985208 24.305 > 17 21 Mg 0.16666666666667 0.33333333333333 0.59819089985208 24.305 > 17 22 Mg 0.41666666666667 0.33333333333333 0.59819089985208 24.305 > 17 23 Mg 0.66666666666667 0.33333333333333 0.59819089985208 24.305 > 17 24 Mg 0.91666666666667 0.33333333333333 0.59819089985208 24.305 > 17 25 Mg 0.16666666666667 0.58333333333333 0.59819089985208 24.305 > 17 26 Mg 0.41666666666667 0.58333333333333 0.59819089985208 24.305 > 17 27 Mg 0.66666666666667 0.58333333333333 0.59819089985208 24.305 > 17 28 Mg 0.91666666666667 0.58333333333333 0.59819089985208 24.305 > 17 29 Mg 0.16666666666667 0.83333333333333 0.59819089985208 24.305 > 17 30 Mg 0.41666666666667 0.83333333333333 0.59819089985208 24.305 > 17 31 Mg 0.66666666666667 0.83333333333333 0.59819089985208 24.305 > 17 32 Mg 0.91666666666667 0.83333333333333 0.59819089985208 24.305 > 17 33 Al 0.08333333333333 0.16666666666667 0.82775598530114 26.982 > 33 34 Al 0.33333333333333 0.16666666666667 0.82775598530114 26.982 > 33 35 Al 0.58333333333333 0.16666666666667 0.82775598530114 26.982 > 33 36 Al 0.83333333333333 0.16666666666667 0.82775598530114 26.982 > 33 37 Al 0.08333333333333 0.41666666666667 0.82775598530114 26.982 > 33 38 Al 0.33333333333333 0.41666666666667 0.82775598530114 26.982 > 33 39 Al 0.58333333333333 0.41666666666667 0.82775598530114 26.982 > 33 40 Al 0.83333333333333 0.41666666666667 0.82775598530114 26.982 > 33 41 Al 0.08333333333333 0.66666666666667 0.82775598530114 26.982 > 33 42 Al 0.33333333333333 0.66666666666667 0.82775598530114 26.982 > 33 43 Al 0.58333333333333 0.66666666666667 0.82775598530114 26.982 > 33 44 Al 0.83333333333333 0.66666666666667 0.82775598530114 26.982 > 33 45 Al 0.08333333333333 0.91666666666667 0.82775598530114 26.982 > 33 46 Al 0.33333333333333 0.91666666666667 0.82775598530114 26.982 > 33 47 Al 0.58333333333333 0.91666666666667 0.82775598530114 26.982 > 33 48 Al 0.83333333333333 0.91666666666667 0.82775598530114 26.982 > 33 49 Al 0.16666666666667 0.08333333333333 0.17224401469886 26.982 > 49 50 Al 0.41666666666667 0.08333333333333 0.17224401469886 26.982 > 49 51 Al 0.66666666666667 0.08333333333333 0.17224401469886 26.982 > 49 52 Al 0.91666666666667 0.08333333333333 0.17224401469886 26.982 > 49 53 Al 0.16666666666667 0.33333333333333 0.17224401469886 26.982 > 49 54 Al 0.41666666666667 0.33333333333333 0.17224401469886 26.982 > 49 55 Al 0.66666666666667 0.33333333333333 0.17224401469886 26.982 > 49 56 Al 0.91666666666667 0.33333333333333 0.17224401469886 26.982 > 49 57 Al 0.16666666666667 0.58333333333333 0.17224401469886 26.982 > 49 58 Al 0.41666666666667 0.58333333333333 0.17224401469886 26.982 > 49 59 Al 0.66666666666667 0.58333333333333 0.17224401469886 26.982 > 49 60 Al 0.91666666666667 0.58333333333333 0.17224401469886 26.982 > 49 61 Al 0.16666666666667 0.83333333333333 0.17224401469886 26.982 > 49 62 Al 0.41666666666667 0.83333333333333 0.17224401469886 26.982 > 49 63 Al 0.66666666666667 0.83333333333333 0.17224401469886 26.982 > 49 64 Al 0.91666666666667 0.83333333333333 0.17224401469886 26.982 > 49 65 Se 0.16666666666667 0.08333333333333 0.89425500396964 78.960 > 65 66 Se 0.41666666666667 0.08333333333333 0.89425500396964 78.960 > 65 67 Se 0.66666666666667 0.08333333333333 0.89425500396964 78.960 > 65 68 Se 0.91666666666667 0.08333333333333 0.89425500396964 78.960 > 65 69 Se 0.16666666666667 0.33333333333333 0.89425500396964 78.960 > 65 70 Se 0.41666666666667 0.33333333333333 0.89425500396964 78.960 > 65 71 Se 0.66666666666667 0.33333333333333 0.89425500396964 78.960 > 65 72 Se 0.91666666666667 0.33333333333333 0.89425500396964 78.960 > 65 73 Se 0.16666666666667 0.58333333333333 0.89425500396964 78.960 > 65 74 Se 0.41666666666667 0.58333333333333 0.89425500396964 78.960 > 65 75 Se 0.66666666666667 0.58333333333333 0.89425500396964 78.960 > 65 76 Se 0.91666666666667 0.58333333333333 0.89425500396964 78.960 > 65 77 Se 0.16666666666667 0.83333333333333 0.89425500396964 78.960 > 65 78 Se 0.41666666666667 0.83333333333333 0.89425500396964 78.960 > 65 79 Se 0.66666666666667 0.83333333333333 0.89425500396964 78.960 > 65 80 Se 0.91666666666667 0.83333333333333 0.89425500396964 78.960 > 65 81 Se 0.08333333333333 0.16666666666667 0.10574499603036 78.960 > 81 82 Se 0.33333333333333 0.16666666666667 0.10574499603036 78.960 > 81 83 Se 0.58333333333333 0.16666666666667 0.10574499603036 78.960 > 81 84 Se 0.83333333333333 0.16666666666667 0.10574499603036 78.960 > 81 85 Se 0.08333333333333 0.41666666666667 0.10574499603036 78.960 > 81 86 Se 0.33333333333333 0.41666666666667 0.10574499603036 78.960 > 81 87 Se 0.58333333333333 0.41666666666667 0.10574499603036 78.960 > 81 88 Se 0.83333333333333 0.41666666666667 0.10574499603036 78.960 > 81 89 Se 0.08333333333333 0.66666666666667 0.10574499603036 78.960 > 81 90 Se 0.33333333333333 0.66666666666667 0.10574499603036 78.960 > 81 91 Se 0.58333333333333 0.66666666666667 0.10574499603036 78.960 > 81 92 Se 0.83333333333333 0.66666666666667 0.10574499603036 78.960 > 81 93 Se 0.08333333333333 0.91666666666667 0.10574499603036 78.960 > 81 94 Se 0.33333333333333 0.91666666666667 0.10574499603036 78.960 > 81 95 Se 0.58333333333333 0.91666666666667 0.10574499603036 78.960 > 81 96 Se 0.83333333333333 0.91666666666667 0.10574499603036 78.960 > 81 97 Se 0.16666666666667 0.08333333333333 0.31447675906161 78.960 > 97 98 Se 0.41666666666667 0.08333333333333 0.31447675906161 78.960 > 97 99 Se 0.66666666666667 0.08333333333333 0.31447675906161 78.960 > 97 100 Se 0.91666666666667 0.08333333333333 0.31447675906161 78.960 > 97 101 Se 0.16666666666667 0.33333333333333 0.31447675906161 78.960 > 97 102 Se 0.41666666666667 0.33333333333333 0.31447675906161 78.960 > 97 103 Se 0.66666666666667 0.33333333333333 0.31447675906161 78.960 > 97 104 Se 0.91666666666667 0.33333333333333 0.31447675906161 78.960 > 97 105 Se 0.16666666666667 0.58333333333333 0.31447675906161 78.960 > 97 106 Se 0.41666666666667 0.58333333333333 0.31447675906161 78.960 > 97 107 Se 0.66666666666667 0.58333333333333 0.31447675906161 78.960 > 97 108 Se 0.91666666666667 0.58333333333333 0.31447675906161 78.960 > 97 109 Se 0.16666666666667 0.83333333333333 0.31447675906161 78.960 > 97 110 Se 0.41666666666667 0.83333333333333 0.31447675906161 78.960 > 97 111 Se 0.66666666666667 0.83333333333333 0.31447675906161 78.960 > 97 112 Se 0.91666666666667 0.83333333333333 0.31447675906161 78.960 > 97 113 Se 0.08333333333333 0.16666666666667 0.68552324093839 78.960 > 113 114 Se 0.33333333333333 0.16666666666667 0.68552324093839 78.960 > 113 115 Se 0.58333333333333 0.16666666666667 0.68552324093839 78.960 > 113 116 Se 0.83333333333333 0.16666666666667 0.68552324093839 78.960 > 113 117 Se 0.08333333333333 0.41666666666667 0.68552324093839 78.960 > 113 118 Se 0.33333333333333 0.41666666666667 0.68552324093839 78.960 > 113 119 Se 0.58333333333333 0.41666666666667 0.68552324093839 78.960 > 113 120 Se 0.83333333333333 0.41666666666667 0.68552324093839 78.960 > 113 121 Se 0.08333333333333 0.66666666666667 0.68552324093839 78.960 > 113 122 Se 0.33333333333333 0.66666666666667 0.68552324093839 78.960 > 113 123 Se 0.58333333333333 0.66666666666667 0.68552324093839 78.960 > 113 124 Se 0.83333333333333 0.66666666666667 0.68552324093839 78.960 > 113 125 Se 0.08333333333333 0.91666666666667 0.68552324093839 78.960 > 113 126 Se 0.33333333333333 0.91666666666667 0.68552324093839 78.960 > 113 127 Se 0.58333333333333 0.91666666666667 0.68552324093839 78.960 > 113 128 Se 0.83333333333333 0.91666666666667 0.68552324093839 78.960 > 113 129 Se 0.00000000000000 0.00000000000000 0.50000000000000 78.960 > 129 130 Se 0.25000000000000 0.00000000000000 0.50000000000000 78.960 > 129 131 Se 0.50000000000000 0.00000000000000 0.50000000000000 78.960 > 129 132 Se 0.75000000000000 0.00000000000000 0.50000000000000 78.960 > 129 133 Se 0.00000000000000 0.25000000000000 0.50000000000000 78.960 > 129 134 Se 0.25000000000000 0.25000000000000 0.50000000000000 78.960 > 129 135 Se 0.50000000000000 0.25000000000000 0.50000000000000 78.960 > 129 136 Se 0.75000000000000 0.25000000000000 0.50000000000000 78.960 > 129 137 Se 0.00000000000000 0.50000000000000 0.50000000000000 78.960 > 129 138 Se 0.25000000000000 0.50000000000000 0.50000000000000 78.960 > 129 139 Se 0.50000000000000 0.50000000000000 0.50000000000000 78.960 > 129 140 Se 0.75000000000000 0.50000000000000 0.50000000000000 78.960 > 129 141 Se 0.00000000000000 0.75000000000000 0.50000000000000 78.960 > 129 142 Se 0.25000000000000 0.75000000000000 0.50000000000000 78.960 > 129 143 Se 0.50000000000000 0.75000000000000 0.50000000000000 78.960 > 129 144 Se 0.75000000000000 0.75000000000000 0.50000000000000 78.960 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 6.8093533 0.0000000 0.0000000 0.0000000 6.8093533 0.0000000 0.0000000 0.0000000 6.0776553 -------------------------- Born effective charges -------------------------- 1 Mg 2.4325488 -0.0000000 0.0000000 0.0000000 2.4325488 0.0000000 0.0000000 0.0000000 1.9913658 2 Mg 2.4325488 -0.0000000 0.0000000 0.0000000 2.4325488 0.0000000 0.0000000 0.0000000 1.9913658 3 Al 2.5785105 -0.0000000 0.0000000 0.0000000 2.5785105 0.0000000 0.0000000 0.0000000 2.0509402 4 Al 2.5785105 -0.0000000 0.0000000 0.0000000 2.5785105 0.0000000 0.0000000 0.0000000 2.0509402 5 Se -2.2486373 0.0000000 0.0000000 0.0000000 -2.2486373 0.0000000 0.0000000 0.0000000 -0.8865659 6 Se -2.2486373 0.0000000 0.0000000 0.0000000 -2.2486373 0.0000000 0.0000000 0.0000000 -0.8865659 7 Se -1.4550392 -0.0000000 0.0000000 0.0000000 -1.4550392 0.0000000 0.0000000 0.0000000 -2.1607832 8 Se -1.4550392 -0.0000000 0.0000000 0.0000000 -1.4550392 0.0000000 0.0000000 0.0000000 -2.1607832 9 Se -2.6147657 0.0000000 0.0000000 0.0000000 -2.6147657 0.0000000 0.0000000 0.0000000 -1.9899137 ---------------------------------------------------------------------------- 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[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 97, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 129, 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.00000003 (zzz) -0.00000003 (zzz) -0.00000003 (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-08 08:17:23]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 08:17:24]------------------------- Compiled with OpenMP support (max 48 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: 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: P-3m1 (164) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.871573573743976 0.000000000000000 0.000000000000000 b -1.935786786871988 3.352881067482788 0.000000000000000 c 0.000000000000000 0.000000000000000 16.139780439999999 Atomic positions (fractional): 1 Mg 0.33333333333333 0.66666666666667 0.40180910014792 24.305 2 Mg 0.66666666666667 0.33333333333333 0.59819089985208 24.305 3 Al 0.33333333333333 0.66666666666667 0.82775598530114 26.982 4 Al 0.66666666666667 0.33333333333333 0.17224401469886 26.982 5 Se 0.66666666666667 0.33333333333333 0.89425500396964 78.960 6 Se 0.33333333333333 0.66666666666667 0.10574499603036 78.960 7 Se 0.66666666666667 0.33333333333333 0.31447675906161 78.960 8 Se 0.33333333333333 0.66666666666667 0.68552324093839 78.960 9 Se 0.00000000000000 0.00000000000000 0.50000000000000 78.960 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.486294294975904 0.000000000000000 0.000000000000000 b -7.743147147487952 13.411524269931151 0.000000000000000 c 0.000000000000000 0.000000000000000 16.139780439999999 Atomic positions (fractional): 1 Mg 0.08333333333333 0.16666666666667 0.40180910014792 24.305 > 1 2 Mg 0.33333333333333 0.16666666666667 0.40180910014792 24.305 > 1 3 Mg 0.58333333333333 0.16666666666667 0.40180910014792 24.305 > 1 4 Mg 0.83333333333333 0.16666666666667 0.40180910014792 24.305 > 1 5 Mg 0.08333333333333 0.41666666666667 0.40180910014792 24.305 > 1 6 Mg 0.33333333333333 0.41666666666667 0.40180910014792 24.305 > 1 7 Mg 0.58333333333333 0.41666666666667 0.40180910014792 24.305 > 1 8 Mg 0.83333333333333 0.41666666666667 0.40180910014792 24.305 > 1 9 Mg 0.08333333333333 0.66666666666667 0.40180910014792 24.305 > 1 10 Mg 0.33333333333333 0.66666666666667 0.40180910014792 24.305 > 1 11 Mg 0.58333333333333 0.66666666666667 0.40180910014792 24.305 > 1 12 Mg 0.83333333333333 0.66666666666667 0.40180910014792 24.305 > 1 13 Mg 0.08333333333333 0.91666666666667 0.40180910014792 24.305 > 1 14 Mg 0.33333333333333 0.91666666666667 0.40180910014792 24.305 > 1 15 Mg 0.58333333333333 0.91666666666667 0.40180910014792 24.305 > 1 16 Mg 0.83333333333333 0.91666666666667 0.40180910014792 24.305 > 1 17 Mg 0.16666666666667 0.08333333333333 0.59819089985208 24.305 > 17 18 Mg 0.41666666666667 0.08333333333333 0.59819089985208 24.305 > 17 19 Mg 0.66666666666667 0.08333333333333 0.59819089985208 24.305 > 17 20 Mg 0.91666666666667 0.08333333333333 0.59819089985208 24.305 > 17 21 Mg 0.16666666666667 0.33333333333333 0.59819089985208 24.305 > 17 22 Mg 0.41666666666667 0.33333333333333 0.59819089985208 24.305 > 17 23 Mg 0.66666666666667 0.33333333333333 0.59819089985208 24.305 > 17 24 Mg 0.91666666666667 0.33333333333333 0.59819089985208 24.305 > 17 25 Mg 0.16666666666667 0.58333333333333 0.59819089985208 24.305 > 17 26 Mg 0.41666666666667 0.58333333333333 0.59819089985208 24.305 > 17 27 Mg 0.66666666666667 0.58333333333333 0.59819089985208 24.305 > 17 28 Mg 0.91666666666667 0.58333333333333 0.59819089985208 24.305 > 17 29 Mg 0.16666666666667 0.83333333333333 0.59819089985208 24.305 > 17 30 Mg 0.41666666666667 0.83333333333333 0.59819089985208 24.305 > 17 31 Mg 0.66666666666667 0.83333333333333 0.59819089985208 24.305 > 17 32 Mg 0.91666666666667 0.83333333333333 0.59819089985208 24.305 > 17 33 Al 0.08333333333333 0.16666666666667 0.82775598530114 26.982 > 33 34 Al 0.33333333333333 0.16666666666667 0.82775598530114 26.982 > 33 35 Al 0.58333333333333 0.16666666666667 0.82775598530114 26.982 > 33 36 Al 0.83333333333333 0.16666666666667 0.82775598530114 26.982 > 33 37 Al 0.08333333333333 0.41666666666667 0.82775598530114 26.982 > 33 38 Al 0.33333333333333 0.41666666666667 0.82775598530114 26.982 > 33 39 Al 0.58333333333333 0.41666666666667 0.82775598530114 26.982 > 33 40 Al 0.83333333333333 0.41666666666667 0.82775598530114 26.982 > 33 41 Al 0.08333333333333 0.66666666666667 0.82775598530114 26.982 > 33 42 Al 0.33333333333333 0.66666666666667 0.82775598530114 26.982 > 33 43 Al 0.58333333333333 0.66666666666667 0.82775598530114 26.982 > 33 44 Al 0.83333333333333 0.66666666666667 0.82775598530114 26.982 > 33 45 Al 0.08333333333333 0.91666666666667 0.82775598530114 26.982 > 33 46 Al 0.33333333333333 0.91666666666667 0.82775598530114 26.982 > 33 47 Al 0.58333333333333 0.91666666666667 0.82775598530114 26.982 > 33 48 Al 0.83333333333333 0.91666666666667 0.82775598530114 26.982 > 33 49 Al 0.16666666666667 0.08333333333333 0.17224401469886 26.982 > 49 50 Al 0.41666666666667 0.08333333333333 0.17224401469886 26.982 > 49 51 Al 0.66666666666667 0.08333333333333 0.17224401469886 26.982 > 49 52 Al 0.91666666666667 0.08333333333333 0.17224401469886 26.982 > 49 53 Al 0.16666666666667 0.33333333333333 0.17224401469886 26.982 > 49 54 Al 0.41666666666667 0.33333333333333 0.17224401469886 26.982 > 49 55 Al 0.66666666666667 0.33333333333333 0.17224401469886 26.982 > 49 56 Al 0.91666666666667 0.33333333333333 0.17224401469886 26.982 > 49 57 Al 0.16666666666667 0.58333333333333 0.17224401469886 26.982 > 49 58 Al 0.41666666666667 0.58333333333333 0.17224401469886 26.982 > 49 59 Al 0.66666666666667 0.58333333333333 0.17224401469886 26.982 > 49 60 Al 0.91666666666667 0.58333333333333 0.17224401469886 26.982 > 49 61 Al 0.16666666666667 0.83333333333333 0.17224401469886 26.982 > 49 62 Al 0.41666666666667 0.83333333333333 0.17224401469886 26.982 > 49 63 Al 0.66666666666667 0.83333333333333 0.17224401469886 26.982 > 49 64 Al 0.91666666666667 0.83333333333333 0.17224401469886 26.982 > 49 65 Se 0.16666666666667 0.08333333333333 0.89425500396964 78.960 > 65 66 Se 0.41666666666667 0.08333333333333 0.89425500396964 78.960 > 65 67 Se 0.66666666666667 0.08333333333333 0.89425500396964 78.960 > 65 68 Se 0.91666666666667 0.08333333333333 0.89425500396964 78.960 > 65 69 Se 0.16666666666667 0.33333333333333 0.89425500396964 78.960 > 65 70 Se 0.41666666666667 0.33333333333333 0.89425500396964 78.960 > 65 71 Se 0.66666666666667 0.33333333333333 0.89425500396964 78.960 > 65 72 Se 0.91666666666667 0.33333333333333 0.89425500396964 78.960 > 65 73 Se 0.16666666666667 0.58333333333333 0.89425500396964 78.960 > 65 74 Se 0.41666666666667 0.58333333333333 0.89425500396964 78.960 > 65 75 Se 0.66666666666667 0.58333333333333 0.89425500396964 78.960 > 65 76 Se 0.91666666666667 0.58333333333333 0.89425500396964 78.960 > 65 77 Se 0.16666666666667 0.83333333333333 0.89425500396964 78.960 > 65 78 Se 0.41666666666667 0.83333333333333 0.89425500396964 78.960 > 65 79 Se 0.66666666666667 0.83333333333333 0.89425500396964 78.960 > 65 80 Se 0.91666666666667 0.83333333333333 0.89425500396964 78.960 > 65 81 Se 0.08333333333333 0.16666666666667 0.10574499603036 78.960 > 81 82 Se 0.33333333333333 0.16666666666667 0.10574499603036 78.960 > 81 83 Se 0.58333333333333 0.16666666666667 0.10574499603036 78.960 > 81 84 Se 0.83333333333333 0.16666666666667 0.10574499603036 78.960 > 81 85 Se 0.08333333333333 0.41666666666667 0.10574499603036 78.960 > 81 86 Se 0.33333333333333 0.41666666666667 0.10574499603036 78.960 > 81 87 Se 0.58333333333333 0.41666666666667 0.10574499603036 78.960 > 81 88 Se 0.83333333333333 0.41666666666667 0.10574499603036 78.960 > 81 89 Se 0.08333333333333 0.66666666666667 0.10574499603036 78.960 > 81 90 Se 0.33333333333333 0.66666666666667 0.10574499603036 78.960 > 81 91 Se 0.58333333333333 0.66666666666667 0.10574499603036 78.960 > 81 92 Se 0.83333333333333 0.66666666666667 0.10574499603036 78.960 > 81 93 Se 0.08333333333333 0.91666666666667 0.10574499603036 78.960 > 81 94 Se 0.33333333333333 0.91666666666667 0.10574499603036 78.960 > 81 95 Se 0.58333333333333 0.91666666666667 0.10574499603036 78.960 > 81 96 Se 0.83333333333333 0.91666666666667 0.10574499603036 78.960 > 81 97 Se 0.16666666666667 0.08333333333333 0.31447675906161 78.960 > 97 98 Se 0.41666666666667 0.08333333333333 0.31447675906161 78.960 > 97 99 Se 0.66666666666667 0.08333333333333 0.31447675906161 78.960 > 97 100 Se 0.91666666666667 0.08333333333333 0.31447675906161 78.960 > 97 101 Se 0.16666666666667 0.33333333333333 0.31447675906161 78.960 > 97 102 Se 0.41666666666667 0.33333333333333 0.31447675906161 78.960 > 97 103 Se 0.66666666666667 0.33333333333333 0.31447675906161 78.960 > 97 104 Se 0.91666666666667 0.33333333333333 0.31447675906161 78.960 > 97 105 Se 0.16666666666667 0.58333333333333 0.31447675906161 78.960 > 97 106 Se 0.41666666666667 0.58333333333333 0.31447675906161 78.960 > 97 107 Se 0.66666666666667 0.58333333333333 0.31447675906161 78.960 > 97 108 Se 0.91666666666667 0.58333333333333 0.31447675906161 78.960 > 97 109 Se 0.16666666666667 0.83333333333333 0.31447675906161 78.960 > 97 110 Se 0.41666666666667 0.83333333333333 0.31447675906161 78.960 > 97 111 Se 0.66666666666667 0.83333333333333 0.31447675906161 78.960 > 97 112 Se 0.91666666666667 0.83333333333333 0.31447675906161 78.960 > 97 113 Se 0.08333333333333 0.16666666666667 0.68552324093839 78.960 > 113 114 Se 0.33333333333333 0.16666666666667 0.68552324093839 78.960 > 113 115 Se 0.58333333333333 0.16666666666667 0.68552324093839 78.960 > 113 116 Se 0.83333333333333 0.16666666666667 0.68552324093839 78.960 > 113 117 Se 0.08333333333333 0.41666666666667 0.68552324093839 78.960 > 113 118 Se 0.33333333333333 0.41666666666667 0.68552324093839 78.960 > 113 119 Se 0.58333333333333 0.41666666666667 0.68552324093839 78.960 > 113 120 Se 0.83333333333333 0.41666666666667 0.68552324093839 78.960 > 113 121 Se 0.08333333333333 0.66666666666667 0.68552324093839 78.960 > 113 122 Se 0.33333333333333 0.66666666666667 0.68552324093839 78.960 > 113 123 Se 0.58333333333333 0.66666666666667 0.68552324093839 78.960 > 113 124 Se 0.83333333333333 0.66666666666667 0.68552324093839 78.960 > 113 125 Se 0.08333333333333 0.91666666666667 0.68552324093839 78.960 > 113 126 Se 0.33333333333333 0.91666666666667 0.68552324093839 78.960 > 113 127 Se 0.58333333333333 0.91666666666667 0.68552324093839 78.960 > 113 128 Se 0.83333333333333 0.91666666666667 0.68552324093839 78.960 > 113 129 Se 0.00000000000000 0.00000000000000 0.50000000000000 78.960 > 129 130 Se 0.25000000000000 0.00000000000000 0.50000000000000 78.960 > 129 131 Se 0.50000000000000 0.00000000000000 0.50000000000000 78.960 > 129 132 Se 0.75000000000000 0.00000000000000 0.50000000000000 78.960 > 129 133 Se 0.00000000000000 0.25000000000000 0.50000000000000 78.960 > 129 134 Se 0.25000000000000 0.25000000000000 0.50000000000000 78.960 > 129 135 Se 0.50000000000000 0.25000000000000 0.50000000000000 78.960 > 129 136 Se 0.75000000000000 0.25000000000000 0.50000000000000 78.960 > 129 137 Se 0.00000000000000 0.50000000000000 0.50000000000000 78.960 > 129 138 Se 0.25000000000000 0.50000000000000 0.50000000000000 78.960 > 129 139 Se 0.50000000000000 0.50000000000000 0.50000000000000 78.960 > 129 140 Se 0.75000000000000 0.50000000000000 0.50000000000000 78.960 > 129 141 Se 0.00000000000000 0.75000000000000 0.50000000000000 78.960 > 129 142 Se 0.25000000000000 0.75000000000000 0.50000000000000 78.960 > 129 143 Se 0.50000000000000 0.75000000000000 0.50000000000000 78.960 > 129 144 Se 0.75000000000000 0.75000000000000 0.50000000000000 78.960 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 6.8093533 0.0000000 0.0000000 0.0000000 6.8093533 0.0000000 0.0000000 0.0000000 6.0776553 -------------------------- Born effective charges -------------------------- 1 Mg 2.4325488 -0.0000000 0.0000000 0.0000000 2.4325488 0.0000000 0.0000000 0.0000000 1.9913658 2 Mg 2.4325488 -0.0000000 0.0000000 0.0000000 2.4325488 0.0000000 0.0000000 0.0000000 1.9913658 3 Al 2.5785105 -0.0000000 0.0000000 0.0000000 2.5785105 0.0000000 0.0000000 0.0000000 2.0509402 4 Al 2.5785105 -0.0000000 0.0000000 0.0000000 2.5785105 0.0000000 0.0000000 0.0000000 2.0509402 5 Se -2.2486373 0.0000000 0.0000000 0.0000000 -2.2486373 0.0000000 0.0000000 0.0000000 -0.8865659 6 Se -2.2486373 0.0000000 0.0000000 0.0000000 -2.2486373 0.0000000 0.0000000 0.0000000 -0.8865659 7 Se -1.4550392 -0.0000000 0.0000000 0.0000000 -1.4550392 0.0000000 0.0000000 0.0000000 -2.1607832 8 Se -1.4550392 -0.0000000 0.0000000 0.0000000 -1.4550392 0.0000000 0.0000000 0.0000000 -2.1607832 9 Se -2.6147657 0.0000000 0.0000000 0.0000000 -2.6147657 0.0000000 0.0000000 0.0000000 -1.9899137 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000003 (zzz) -0.00000003 (zzz) -0.00000003 (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 ... [ 15 15 3 ] fc3-r2q-transformation over three atoms: True --------------------------- Phonon calculations ---------------------------- Use NAC by Gonze et al. (no real space sum in current implementation) PRB 50, 13035(R) (1994), PRB 55, 10355 (1997) G-cutoff distance: 0.70, Number of G-points: 305, Lambda: 0.19 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/73) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 73 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.000 ( 0.000 0.000 0.000) 0.000 0.000 ( 0.000 0.000 0.000) 0.000 0.000 ( 0.000 0.000 0.000) 0.000 1.008 ( 0.000 0.000 0.000) 0.000 1.008 ( 0.000 0.000 0.000) 0.000 1.262 ( 0.000 0.000 0.000) 0.000 1.262 ( 0.000 0.000 0.000) 0.000 1.959 ( 0.000 0.000 0.000) 0.000 1.959 ( 0.000 0.000 0.000) 0.000 1.980 ( 0.000 0.000 0.000) 0.000 2.878 ( 0.000 0.000 0.000) 0.000 2.878 ( 0.000 0.000 0.000) 0.000 3.738 ( 0.000 0.000 0.000) 0.000 5.453 ( 0.000 0.000 0.000) 0.000 5.682 ( 0.000 0.000 0.000) 0.000 6.804 ( 0.000 0.000 0.000) 0.000 6.804 ( 0.000 0.000 0.000) 0.000 6.908 ( 0.000 0.000 0.000) 0.000 6.908 ( 0.000 0.000 0.000) 0.000 8.310 ( 0.000 0.000 0.000) 0.000 8.310 ( 0.000 0.000 0.000) 0.000 8.312 ( 0.000 0.000 0.000) 0.000 8.312 ( 0.000 0.000 0.000) 0.000 8.394 ( 0.000 0.000 0.000) 0.000 9.590 ( 0.000 0.000 0.000) 0.000 13.465 ( 0.000 0.000 0.000) 0.000 13.981 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/73) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.334 ( 18.740 10.820 0.310) 21.642 0.484 ( 20.527 11.851 4.345) 24.097 0.914 ( 36.286 20.950 -2.281) 41.961 1.152 ( 11.567 6.678 -2.404) 13.571 1.384 ( 9.959 5.750 -0.309) 11.504 1.437 ( 29.290 16.910 0.987) 33.835 1.576 ( 21.362 12.333 0.313) 24.668 1.986 ( 2.419 1.397 0.586) 2.854 2.049 ( 8.487 4.900 0.660) 9.822 2.107 ( 12.551 7.246 -1.012) 14.528 2.905 ( 2.281 1.317 -0.018) 2.634 3.118 ( 13.344 7.704 0.255) 15.410 3.764 ( 2.453 1.416 -0.019) 2.832 5.393 ( -4.861 -2.806 0.107) 5.614 5.606 ( -6.210 -3.586 -0.108) 7.172 6.801 ( -0.053 -0.030 0.005) 0.061 6.914 ( 0.425 0.245 -0.002) 0.491 6.923 ( 0.623 0.360 -0.074) 0.723 7.447 ( 11.020 6.363 2.380) 12.946 8.034 ( -12.032 -6.947 -1.786) 14.008 8.252 ( -3.490 -2.015 -0.199) 4.034 8.313 ( -1.531 -0.884 0.193) 1.778 8.770 ( 19.888 11.482 2.338) 23.083 9.406 ( -15.132 -8.737 0.011) 17.473 10.047 ( 2.284 1.319 -3.250) 4.185 13.355 ( -9.130 -5.271 0.804) 10.573 13.828 ( -12.555 -7.249 -0.101) 14.498 ======================= Grid point 2 (3/73) ======================= q-point: ( 0.13 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.814 ( 22.139 12.782 -0.360) 25.566 0.928 ( 17.779 10.265 2.727) 20.710 1.476 ( 15.100 8.718 -2.292) 17.586 1.529 ( 14.922 8.615 -0.435) 17.236 1.677 ( 14.617 8.439 -0.745) 16.894 2.022 ( 18.447 10.651 -0.024) 21.301 2.072 ( 5.199 3.002 1.070) 6.098 2.097 ( 25.482 14.712 0.539) 29.429 2.479 ( 26.891 15.526 1.865) 31.107 2.563 ( 26.785 15.464 -2.150) 31.004 2.976 ( 3.662 2.115 -0.030) 4.229 3.474 ( 15.655 9.039 0.048) 18.077 3.867 ( 7.046 4.068 0.054) 8.136 5.253 ( -6.772 -3.910 0.259) 7.824 5.426 ( -8.787 -5.073 -0.241) 10.150 6.812 ( 1.150 0.664 0.003) 1.327 6.883 ( -5.272 -3.044 0.126) 6.089 6.923 ( 0.308 0.178 0.008) 0.356 7.666 ( 3.965 2.289 0.788) 4.645 7.826 ( -4.079 -2.355 -0.239) 4.716 8.159 ( -4.471 -2.581 -0.322) 5.173 8.238 ( -4.881 -2.818 0.306) 5.644 8.919 ( -19.731 -11.392 -0.821) 22.798 9.271 ( 17.693 10.215 4.750) 20.975 10.118 ( 3.787 2.186 -5.342) 6.903 13.081 ( -13.839 -7.990 1.567) 16.057 13.455 ( -18.733 -10.816 -0.584) 21.639 ======================= Grid point 3 (4/73) ======================= q-point: ( 0.20 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.286 ( 13.085 7.555 0.926) 15.138 1.320 ( 21.572 12.455 -0.886) 24.925 1.738 ( 6.724 3.882 1.268) 7.867 1.780 ( 10.349 5.975 -0.337) 11.955 2.026 ( 15.405 8.894 -1.743) 17.873 2.236 ( 9.153 5.284 1.340) 10.653 2.292 ( 7.449 4.301 0.785) 8.637 2.735 ( 26.908 15.535 -1.411) 31.102 3.065 ( 4.007 2.313 -0.013) 4.627 3.088 ( 22.641 13.072 -2.701) 26.282 3.192 ( 25.470 14.705 2.167) 29.490 3.775 ( 10.874 6.278 0.521) 12.567 4.089 ( 10.984 6.342 0.197) 12.685 5.107 ( -5.487 -3.168 -0.562) 6.361 5.248 ( -5.325 -3.075 0.398) 6.162 6.682 ( -10.892 -6.288 0.007) 12.576 6.855 ( 2.466 1.424 -0.015) 2.848 6.930 ( 0.323 0.186 0.015) 0.373 7.595 ( -12.388 -7.152 -0.302) 14.308 7.721 ( -5.474 -3.160 0.682) 6.357 8.053 ( -4.616 -2.665 -0.135) 5.332 8.098 ( -6.957 -4.017 0.132) 8.035 8.597 ( -3.918 -2.262 -0.952) 4.623 9.648 ( 13.622 7.865 5.926) 16.809 10.212 ( 4.165 2.405 -6.146) 7.804 12.760 ( -13.330 -7.696 2.208) 15.550 13.021 ( -18.080 -10.439 -1.302) 20.918 ======================= Grid point 4 (5/73) ======================= q-point: ( 0.27 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.520 ( 7.237 4.178 -0.134) 8.358 1.782 ( 17.527 10.119 -2.390) 20.379 1.906 ( 9.023 5.209 3.185) 10.895 1.952 ( 5.366 3.098 0.829) 6.252 2.360 ( 13.114 7.571 -2.617) 15.367 2.430 ( 5.448 3.146 0.359) 6.302 2.463 ( 9.553 5.515 1.927) 11.198 3.156 ( 3.815 2.203 0.010) 4.405 3.189 ( 8.666 5.003 -2.149) 10.234 3.471 ( 13.013 7.513 -3.504) 15.429 3.665 ( 15.219 8.786 3.250) 17.871 4.027 ( 10.969 6.333 1.585) 12.764 4.273 ( 3.368 1.944 0.078) 3.890 5.044 ( 1.681 0.970 -1.871) 2.696 5.233 ( 3.938 2.274 1.312) 4.733 6.456 ( -7.513 -4.338 -0.055) 8.676 6.916 ( 2.533 1.462 -0.026) 2.925 6.942 ( 0.960 0.554 0.026) 1.109 7.124 ( -21.333 -12.317 0.421) 24.637 7.695 ( 0.524 0.303 -0.254) 0.656 7.935 ( -6.471 -3.736 -0.141) 7.473 7.955 ( -4.156 -2.399 0.139) 4.801 8.649 ( 5.384 3.109 -1.159) 6.325 9.887 ( 7.465 4.310 7.024) 11.120 10.290 ( 2.065 1.192 -6.617) 7.034 12.496 ( -9.308 -5.374 2.584) 11.054 12.659 ( -12.941 -7.472 -1.899) 15.063 ======================= Grid point 5 (6/73) ======================= q-point: ( 0.33 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.625 ( 2.333 1.347 -0.271) 2.708 2.043 ( 8.043 4.644 -2.660) 9.661 2.050 ( 3.356 1.937 0.618) 3.924 2.219 ( 14.834 8.564 3.078) 17.403 2.571 ( 6.941 4.007 0.079) 8.015 2.603 ( 7.412 4.279 -1.838) 8.754 2.651 ( 7.072 4.083 1.439) 8.292 3.177 ( -5.235 -3.022 -2.823) 6.672 3.237 ( 3.167 1.828 0.055) 3.657 3.655 ( 1.187 0.685 1.862) 2.312 3.904 ( 4.547 2.625 -0.607) 5.285 4.170 ( -2.569 -1.483 0.558) 3.019 4.314 ( 6.481 3.742 0.658) 7.512 5.198 ( 10.115 5.840 -1.990) 11.848 5.388 ( 8.638 4.987 1.560) 10.096 6.365 ( -0.373 -0.215 -0.019) 0.431 6.762 ( -9.576 -5.528 0.134) 11.058 6.955 ( 1.313 0.758 0.000) 1.516 6.984 ( 2.083 1.202 0.009) 2.405 7.713 ( 0.944 0.545 -0.048) 1.091 7.808 ( -4.224 -2.439 -0.333) 4.889 7.873 ( -3.045 -1.758 0.317) 3.531 8.781 ( 5.612 3.240 -1.336) 6.616 10.011 ( 3.777 2.181 8.346) 9.417 10.290 ( -2.112 -1.219 -7.577) 7.960 12.339 ( -4.426 -2.556 2.513) 5.696 12.434 ( -6.778 -3.913 -2.056) 8.093 ======================= Grid point 6 (7/73) ======================= q-point: ( 0.40 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.648 ( 0.075 0.043 -0.123) 0.150 2.108 ( 1.707 0.986 0.150) 1.977 2.211 ( 6.474 3.738 -2.789) 7.978 2.504 ( 9.602 5.544 2.025) 11.271 2.700 ( 2.467 1.424 -0.406) 2.877 2.704 ( 2.200 1.270 1.186) 2.803 2.789 ( 4.498 2.597 -1.377) 5.373 3.108 ( 1.820 1.051 -0.676) 2.207 3.295 ( 1.760 1.016 0.172) 2.040 3.583 ( -5.684 -3.282 2.153) 6.907 3.846 ( -6.912 -3.991 -0.757) 8.017 4.034 ( -7.517 -4.340 0.181) 8.682 4.522 ( 8.943 5.163 0.033) 10.326 5.427 ( 8.653 4.996 -1.602) 10.119 5.591 ( 8.180 4.723 1.409) 9.550 6.416 ( 3.869 2.234 0.005) 4.467 6.650 ( -1.575 -0.910 -0.017) 1.819 6.988 ( 1.504 0.868 0.013) 1.737 7.021 ( 1.094 0.632 0.003) 1.263 7.735 ( 0.964 0.556 0.037) 1.113 7.745 ( -1.388 -0.801 -0.429) 1.659 7.814 ( -2.204 -1.273 0.405) 2.577 8.898 ( 4.456 2.573 -0.796) 5.207 10.072 ( 1.700 0.982 8.585) 8.807 10.208 ( -4.429 -2.557 -8.133) 9.607 12.280 ( -1.093 -0.631 1.881) 2.266 12.333 ( -2.445 -1.412 -1.642) 3.266 ======================= Grid point 7 (8/73) ======================= q-point: ( 0.47 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.644 ( -0.191 -0.111 -0.032) 0.223 2.132 ( 0.457 0.264 -0.011) 0.528 2.322 ( 2.737 1.580 -1.402) 3.457 2.619 ( 0.336 0.194 -0.358) 0.528 2.730 ( 0.398 0.230 0.708) 0.844 2.736 ( 2.479 1.431 0.659) 2.937 2.855 ( 1.375 0.794 -0.784) 1.771 3.223 ( 4.853 2.802 0.824) 5.664 3.319 ( 0.452 0.261 0.121) 0.536 3.460 ( -3.692 -2.131 0.327) 4.275 3.719 ( -2.522 -1.456 -0.089) 2.914 3.860 ( -6.194 -3.576 -0.044) 7.152 4.669 ( 3.395 1.960 -0.082) 3.921 5.567 ( 3.191 1.842 -0.659) 3.743 5.729 ( 3.293 1.901 0.601) 3.850 6.499 ( 2.494 1.440 0.007) 2.880 6.636 ( -0.245 -0.142 -0.029) 0.285 7.019 ( 1.059 0.611 0.012) 1.222 7.033 ( -0.055 -0.032 -0.001) 0.063 7.736 ( 0.369 0.213 -0.409) 0.590 7.753 ( 0.446 0.258 0.019) 0.516 7.772 ( -1.347 -0.777 0.395) 1.604 8.975 ( 1.901 1.097 -0.131) 2.199 10.091 ( -0.150 -0.087 6.054) 6.057 10.125 ( -1.764 -1.019 -6.027) 6.362 12.271 ( -0.029 -0.016 0.728) 0.729 12.303 ( -0.473 -0.273 -0.660) 0.857 ======================= Grid point 16 (9/73) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.664 ( 12.321 21.340 0.000) 24.641 0.896 ( 13.065 22.629 0.000) 26.129 1.351 ( 12.740 22.066 0.000) 25.480 1.446 ( 8.608 14.910 0.000) 17.217 1.599 ( 9.512 16.475 0.000) 19.024 1.884 ( 12.206 21.141 0.000) 24.412 1.919 ( 10.220 17.701 0.000) 20.439 2.090 ( 6.972 12.076 0.000) 13.945 2.337 ( 10.788 18.685 0.000) 21.576 2.351 ( 14.086 24.398 0.000) 28.173 2.973 ( 2.913 5.045 0.000) 5.826 3.364 ( 8.981 15.556 0.000) 17.962 3.823 ( 2.768 4.794 0.000) 5.535 5.295 ( -3.767 -6.524 0.000) 7.534 5.478 ( -5.050 -8.747 0.000) 10.101 6.822 ( 0.830 1.437 0.000) 1.659 6.901 ( 0.056 0.096 0.000) 0.111 6.909 ( -2.221 -3.847 0.000) 4.442 7.639 ( 5.927 10.267 0.000) 11.855 7.853 ( -4.998 -8.656 0.000) 9.995 8.231 ( -2.015 -3.490 0.000) 4.030 8.234 ( -1.956 -3.389 0.000) 3.913 9.089 ( -12.287 -21.282 0.000) 24.575 9.117 ( 11.461 19.851 0.000) 22.922 10.086 ( 1.479 2.562 0.000) 2.958 13.162 ( -7.672 -13.287 0.000) 15.343 13.567 ( -10.330 -17.892 0.000) 20.659 ======================= Grid point 17 (10/73) ======================= q-point: ( 0.13 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.124 ( 17.111 18.191 -0.374) 24.977 1.294 ( 6.970 22.190 0.992) 23.280 1.636 ( 7.117 4.393 0.368) 8.372 1.789 ( 8.240 19.920 -0.937) 21.577 1.973 ( 11.169 20.268 -0.517) 23.148 2.135 ( 6.427 6.973 -0.593) 9.502 2.203 ( 11.078 6.210 1.358) 12.773 2.530 ( 16.772 23.543 0.320) 28.908 2.822 ( 20.455 16.172 -2.029) 26.154 2.930 ( 22.724 20.911 1.373) 30.911 3.093 ( 3.219 9.834 -0.107) 10.348 3.672 ( 9.987 10.378 0.082) 14.403 3.966 ( 8.094 5.378 0.093) 9.719 5.159 ( -4.644 -4.852 -0.041) 6.716 5.289 ( -5.882 -8.132 0.009) 10.036 6.742 ( -6.780 -10.185 0.038) 12.236 6.864 ( 0.070 3.427 0.000) 3.428 6.918 ( 1.388 -0.015 -0.022) 1.389 7.707 ( -4.415 1.484 -0.667) 4.705 7.763 ( -2.598 -5.055 0.970) 5.765 8.126 ( -5.096 -1.933 -0.009) 5.450 8.160 ( -3.936 -4.089 -0.323) 5.685 8.674 ( -11.170 -11.106 -0.096) 15.751 9.497 ( 13.068 11.560 2.943) 17.693 10.146 ( 3.590 0.444 -3.272) 4.878 12.864 ( -10.614 -13.072 1.084) 16.874 13.168 ( -14.400 -17.124 -0.563) 22.381 ======================= Grid point 18 (11/73) ======================= q-point: ( 0.20 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.573 ( 4.529 19.055 0.207) 19.587 1.594 ( 15.169 14.716 -1.063) 21.161 1.795 ( 7.101 4.132 1.555) 8.362 2.005 ( -0.396 12.184 0.098) 12.191 2.331 ( 8.715 6.796 1.266) 11.124 2.353 ( 6.869 21.501 -2.455) 22.705 2.460 ( 7.296 13.429 1.405) 15.347 3.004 ( 15.904 9.533 -0.652) 18.554 3.140 ( 6.298 4.927 -1.400) 8.118 3.356 ( 10.649 13.146 -2.688) 17.130 3.470 ( 18.736 12.735 2.898) 22.839 3.909 ( 9.156 6.943 0.945) 11.529 4.170 ( 9.220 1.395 0.111) 9.326 5.055 ( -2.963 -1.573 -1.364) 3.621 5.179 ( 1.470 -3.379 1.062) 3.835 6.512 ( -6.864 -9.349 -0.042) 11.598 6.911 ( 1.051 2.844 0.012) 3.032 6.934 ( -0.000 0.557 -0.021) 0.558 7.418 ( -21.510 -9.319 0.656) 23.451 7.690 ( -0.444 0.034 -0.497) 0.667 8.031 ( -4.942 -0.393 0.165) 4.961 8.040 ( -6.886 -0.920 -0.543) 6.968 8.576 ( 3.476 -0.609 -0.132) 3.531 9.772 ( 9.566 4.294 4.603) 11.451 10.212 ( 4.886 -2.587 -4.729) 7.276 12.571 ( -8.869 -10.824 1.855) 14.116 12.773 ( -12.507 -13.855 -1.272) 18.708 ======================= Grid point 19 (12/73) ======================= q-point: ( 0.27 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.751 ( -3.169 17.551 -0.339) 17.838 1.932 ( 7.360 5.139 -1.961) 9.188 2.033 ( 14.092 10.316 1.717) 17.549 2.096 ( 1.894 7.401 1.199) 7.733 2.518 ( 5.734 7.454 0.694) 9.430 2.643 ( 2.007 18.641 -1.941) 18.849 2.692 ( 1.991 16.768 1.525) 16.954 3.190 ( 1.120 -2.510 -2.246) 3.550 3.251 ( 2.150 6.204 -0.464) 6.583 3.603 ( 6.612 4.717 -0.231) 8.125 3.807 ( 10.029 5.038 1.083) 11.276 4.124 ( 6.893 2.325 1.149) 7.364 4.255 ( 3.292 -1.707 0.101) 3.709 5.092 ( 6.732 3.517 -1.896) 7.828 5.264 ( 8.626 0.846 1.433) 8.785 6.364 ( -1.944 -4.134 -0.038) 4.569 6.835 ( -6.872 -8.951 0.126) 11.285 6.975 ( 0.597 4.208 -0.017) 4.250 7.087 ( -10.011 5.884 0.153) 11.613 7.693 ( 0.839 -0.461 -0.221) 0.983 7.910 ( -7.652 1.013 -0.357) 7.727 7.946 ( -5.283 1.032 0.304) 5.391 8.676 ( 7.400 -0.684 -0.783) 7.472 9.923 ( 6.475 -0.729 6.168) 8.972 10.243 ( 2.993 -5.763 -5.801) 8.708 12.361 ( -4.485 -7.854 2.189) 9.306 12.484 ( -7.152 -9.863 -1.727) 12.305 ======================= Grid point 20 (13/73) ======================= q-point: ( 0.33 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.822 ( -8.088 17.232 -0.242) 19.038 2.070 ( 4.260 3.230 -1.039) 5.446 2.197 ( 2.814 7.610 -1.296) 8.217 2.372 ( 12.315 6.924 2.368) 14.325 2.681 ( 4.056 7.017 0.101) 8.105 2.786 ( -2.485 12.494 0.522) 12.749 2.856 ( -1.551 14.328 -0.514) 14.421 3.119 ( -2.785 -0.204 -2.005) 3.438 3.349 ( -0.636 9.169 0.171) 9.193 3.638 ( -3.184 -1.718 2.395) 4.339 3.884 ( -1.716 -4.246 -1.226) 4.740 4.085 ( -4.077 -6.013 0.353) 7.274 4.414 ( 9.322 5.745 0.337) 10.956 5.298 ( 10.716 3.734 -1.704) 11.476 5.453 ( 10.908 1.505 1.432) 11.104 6.360 ( 3.413 -0.189 -0.011) 3.418 6.664 ( -3.179 -4.237 0.025) 5.297 7.018 ( -0.430 4.373 -0.004) 4.394 7.068 ( -3.107 7.237 0.029) 7.876 7.707 ( 1.646 -1.060 -0.002) 1.958 7.810 ( -5.457 2.010 -0.406) 5.830 7.872 ( -4.484 1.122 0.392) 4.639 8.794 ( 7.582 -2.046 -0.857) 7.899 9.998 ( 4.977 -3.557 7.214) 9.459 10.196 ( -0.274 -7.784 -6.736) 10.298 12.257 ( -0.226 -5.470 1.933) 5.806 12.327 ( -1.868 -6.661 -1.638) 7.109 ======================= Grid point 21 (14/73) ======================= q-point: ( 0.40 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.835 ( -9.955 17.437 -0.062) 20.078 2.142 ( 0.437 3.165 -0.454) 3.227 2.313 ( 2.412 5.498 -1.809) 6.271 2.593 ( 6.884 3.895 0.993) 7.972 2.740 ( -3.154 3.291 -0.341) 4.572 2.831 ( -3.257 9.233 0.919) 9.834 2.947 ( -3.126 12.541 -0.551) 12.937 3.168 ( 3.665 5.851 0.229) 6.908 3.403 ( -3.422 9.028 0.158) 9.656 3.522 ( -6.107 -2.460 1.083) 6.672 3.762 ( -4.691 -3.264 -0.247) 5.720 3.930 ( -6.620 -5.424 -0.056) 8.558 4.604 ( 6.574 2.844 -0.072) 7.163 5.483 ( 7.533 0.558 -1.125) 7.637 5.632 ( 8.353 -0.562 0.995) 8.430 6.442 ( 4.831 0.354 0.007) 4.844 6.618 ( 0.479 -2.484 -0.034) 2.530 7.049 ( -0.909 4.060 0.013) 4.161 7.086 ( -2.907 6.394 0.007) 7.023 7.726 ( 1.794 -1.418 0.053) 2.287 7.765 ( -2.385 2.192 -0.491) 3.277 7.815 ( -3.309 1.051 0.462) 3.503 8.888 ( 6.365 -3.725 -0.324) 7.382 10.032 ( 4.091 -5.007 6.817) 9.396 10.103 ( -0.789 -7.423 -6.689) 10.023 12.228 ( 2.020 -4.462 1.089) 5.017 12.269 ( 1.214 -4.843 -0.960) 5.085 ======================= Grid point 22 (15/73) ======================= q-point: ( 0.47 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.835 ( -10.248 17.750 0.000) 20.496 2.159 ( -1.513 2.621 0.000) 3.026 2.366 ( -1.384 2.398 -0.000) 2.769 2.674 ( -3.087 5.347 -0.000) 6.174 2.717 ( 0.541 -0.938 -0.000) 1.083 2.857 ( -5.079 8.797 0.000) 10.158 2.982 ( -6.482 11.227 0.000) 12.964 3.253 ( -0.851 1.473 0.000) 1.701 3.420 ( -4.813 8.337 -0.000) 9.627 3.439 ( -0.113 0.195 -0.000) 0.225 3.714 ( -0.671 1.162 0.000) 1.342 3.816 ( 0.068 -0.118 0.000) 0.136 4.677 ( 0.678 -1.174 -0.000) 1.356 5.553 ( 1.816 -3.146 0.000) 3.632 5.704 ( 2.437 -4.222 -0.000) 4.875 6.494 ( 1.094 -1.895 -0.000) 2.188 6.609 ( 1.397 -2.419 0.000) 2.793 7.063 ( -1.776 3.077 0.000) 3.553 7.093 ( -3.564 6.173 -0.000) 7.128 7.735 ( 1.085 -1.879 0.000) 2.169 7.761 ( -1.048 1.815 0.000) 2.096 7.787 ( -1.272 2.204 -0.000) 2.545 8.927 ( 3.400 -5.888 -0.000) 6.799 10.046 ( 2.385 -4.132 -0.000) 4.771 10.049 ( 3.756 -6.505 0.000) 7.511 12.225 ( 2.517 -4.359 -0.000) 5.034 12.258 ( 2.375 -4.114 -0.000) 4.751 ======================= Grid point 33 (16/73) ======================= q-point: ( 0.13 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.508 ( 11.455 19.840 0.000) 22.909 1.693 ( 9.350 16.194 0.000) 18.699 1.743 ( 3.936 6.818 0.000) 7.872 2.143 ( 7.552 13.081 0.000) 15.105 2.238 ( 3.153 5.461 0.000) 6.305 2.442 ( 11.444 19.822 0.000) 22.889 2.443 ( 13.147 22.772 0.000) 26.294 2.923 ( 8.258 14.303 0.000) 16.516 3.066 ( 5.171 8.957 0.000) 10.342 3.349 ( 11.608 20.106 0.000) 23.217 3.354 ( 8.427 14.596 0.000) 16.854 3.870 ( 5.383 9.323 0.000) 10.766 4.071 ( 2.916 5.051 0.000) 5.832 5.089 ( -1.041 -1.803 0.000) 2.081 5.134 ( -3.569 -6.181 0.000) 7.137 6.521 ( -6.144 -10.642 0.000) 12.288 6.919 ( -0.186 -0.323 0.000) 0.373 6.940 ( 2.112 3.658 0.000) 4.224 7.562 ( -8.594 -14.885 0.000) 17.187 7.718 ( -0.213 -0.370 0.000) 0.427 8.063 ( -1.397 -2.419 0.000) 2.793 8.109 ( -1.340 -2.322 0.000) 2.681 8.545 ( -1.193 -2.067 0.000) 2.387 9.699 ( 4.776 8.273 0.000) 9.552 10.140 ( -0.430 -0.746 0.000) 0.861 12.592 ( -7.806 -13.520 0.000) 15.611 12.812 ( -10.213 -17.690 0.000) 20.426 ======================= Grid point 34 (17/73) ======================= q-point: ( 0.20 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.878 ( 7.538 12.174 -2.774) 14.585 1.909 ( 9.416 11.649 2.412) 15.172 1.942 ( 4.595 11.641 0.584) 12.528 2.248 ( -6.163 15.239 -0.387) 16.442 2.467 ( 11.337 7.069 0.834) 13.387 2.830 ( 4.611 22.911 -0.016) 23.371 2.850 ( 5.081 24.375 -0.026) 24.899 3.094 ( 6.028 -0.591 -0.571) 6.084 3.223 ( 6.371 4.264 -1.066) 7.740 3.589 ( 3.784 9.334 -0.270) 10.076 3.687 ( 10.132 8.186 0.892) 13.056 4.054 ( 5.196 7.715 0.607) 9.322 4.149 ( 5.252 -2.514 -0.089) 5.823 5.060 ( 0.154 2.235 -1.493) 2.692 5.133 ( 5.820 -0.766 1.304) 6.014 6.356 ( -3.036 -5.780 -0.022) 6.529 6.850 ( -4.580 -7.372 0.038) 8.679 7.011 ( 1.148 4.782 -0.007) 4.918 7.301 ( -8.456 -1.442 0.221) 8.581 7.700 ( -1.663 0.465 -0.272) 1.748 8.023 ( -2.804 -0.072 -0.572) 2.863 8.048 ( -4.569 0.356 0.272) 4.591 8.555 ( 3.626 -1.553 0.284) 3.955 9.818 ( 4.125 0.410 1.926) 4.571 10.123 ( 2.025 -5.762 -2.062) 6.446 12.353 ( -5.570 -10.513 0.882) 11.930 12.495 ( -7.946 -13.502 -0.681) 15.681 ======================= Grid point 35 (18/73) ======================= q-point: ( 0.27 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.042 ( 3.203 7.738 -1.573) 8.521 2.116 ( 2.784 13.176 -0.472) 13.476 2.231 ( 8.473 9.013 2.731) 12.668 2.307 ( -3.389 17.461 -0.936) 17.811 2.689 ( 5.085 8.379 0.809) 9.834 2.952 ( -3.487 3.838 -0.950) 5.272 3.070 ( -3.295 18.268 0.912) 18.585 3.214 ( -0.177 14.748 -1.864) 14.867 3.414 ( 6.716 9.415 0.202) 11.566 3.668 ( -1.139 2.815 1.248) 3.283 3.832 ( 3.075 -2.444 -0.717) 3.993 4.060 ( -1.134 -6.628 0.328) 6.733 4.325 ( 8.587 7.447 0.325) 11.371 5.167 ( 7.552 3.707 -1.435) 8.534 5.273 ( 9.985 0.249 1.211) 10.061 6.307 ( 1.435 -1.646 -0.024) 2.184 6.675 ( -4.575 -7.069 0.037) 8.420 7.072 ( 0.053 5.265 0.004) 5.265 7.237 ( -5.111 8.137 0.048) 9.609 7.680 ( -0.329 -0.769 -0.118) 0.845 7.947 ( -6.690 1.732 -0.570) 6.934 7.970 ( -4.889 0.775 0.513) 4.976 8.621 ( 7.235 -4.640 -0.152) 8.596 9.864 ( 4.287 -4.851 3.335) 7.282 10.085 ( 2.003 -9.263 -3.294) 10.033 12.201 ( -1.532 -7.787 1.232) 8.032 12.283 ( -3.142 -9.868 -1.029) 10.408 ======================= Grid point 36 (19/73) ======================= q-point: ( 0.33 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.156 ( 0.488 8.552 -1.206) 8.650 2.248 ( -7.321 18.776 -0.167) 20.154 2.345 ( 0.821 7.818 -0.569) 7.882 2.511 ( 6.870 8.288 1.057) 10.817 2.824 ( 2.618 6.128 0.753) 6.707 2.910 ( -2.050 0.366 -0.906) 2.271 3.170 ( -4.523 17.111 1.187) 17.739 3.322 ( -6.210 17.204 -1.157) 18.327 3.560 ( -2.612 9.074 0.281) 9.447 3.625 ( -5.799 2.611 1.315) 6.494 3.789 ( -1.487 -3.126 -0.758) 3.544 3.966 ( -3.187 -4.669 -0.006) 5.653 4.520 ( 6.646 4.382 0.081) 7.961 5.340 ( 8.832 0.500 -1.151) 8.921 5.452 ( 10.537 -1.461 0.973) 10.683 6.355 ( 4.679 -0.436 -0.006) 4.699 6.581 ( -0.345 -4.369 -0.015) 4.383 7.115 ( -1.067 5.139 0.015) 5.248 7.242 ( -4.899 9.444 -0.017) 10.639 7.676 ( 1.404 -1.877 0.068) 2.345 7.855 ( -5.319 1.644 -0.419) 5.583 7.896 ( -4.463 0.806 0.399) 4.553 8.699 ( 8.344 -7.360 -0.227) 11.129 9.880 ( 4.948 -7.798 3.440) 9.855 10.011 ( 0.977 -9.968 -3.383) 10.571 12.135 ( 1.670 -6.363 0.871) 6.636 12.183 ( 1.068 -7.465 -0.747) 7.578 ======================= Grid point 37 (20/73) ======================= q-point: ( 0.40 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.225 ( -3.607 11.091 -0.287) 11.666 2.261 ( -10.349 16.075 -0.081) 19.119 2.415 ( -0.484 4.875 -0.728) 4.953 2.687 ( 4.177 6.943 0.222) 8.106 2.826 ( -6.525 4.624 -0.249) 8.001 2.959 ( 2.638 2.477 0.971) 3.746 3.243 ( -6.712 14.829 0.227) 16.279 3.368 ( -3.958 12.341 -0.073) 12.960 3.484 ( -6.576 0.222 -0.395) 6.591 3.586 ( -4.962 6.154 0.469) 7.920 3.786 ( -2.658 7.914 0.054) 8.349 3.866 ( -5.361 0.803 -0.165) 5.423 4.629 ( 3.523 -0.393 -0.032) 3.545 5.453 ( 5.484 -3.374 -0.478) 6.456 5.575 ( 6.751 -4.929 0.406) 8.369 6.429 ( 4.087 -1.602 0.006) 4.390 6.550 ( 1.328 -4.293 -0.020) 4.494 7.137 ( -2.024 4.628 0.009) 5.052 7.249 ( -5.067 9.194 -0.009) 10.498 7.684 ( 1.879 -2.602 0.050) 3.210 7.806 ( -2.179 1.420 -0.271) 2.615 7.844 ( -2.626 1.540 0.258) 3.055 8.751 ( 7.313 -9.703 -0.048) 12.151 9.883 ( 5.330 -9.126 1.525) 10.678 9.943 ( 2.327 -8.178 -1.546) 8.642 12.117 ( 3.156 -6.137 0.208) 6.904 12.154 ( 3.091 -6.299 -0.170) 7.019 ======================= Grid point 49 (21/73) ======================= q-point: ( 0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.083 ( 5.882 10.188 0.000) 11.764 2.084 ( 3.700 6.409 0.000) 7.401 2.198 ( 7.995 13.848 0.000) 15.991 2.581 ( 5.664 9.810 0.000) 11.328 2.612 ( 7.178 12.432 0.000) 14.355 2.977 ( -3.275 -5.673 0.000) 6.551 3.281 ( 5.135 8.893 0.000) 10.269 3.294 ( 11.464 19.856 0.000) 22.928 3.340 ( 8.772 15.194 0.000) 17.545 3.723 ( 2.243 3.885 0.000) 4.486 3.762 ( 0.121 0.209 0.000) 0.241 4.032 ( -3.321 -5.752 0.000) 6.641 4.281 ( 6.517 11.287 0.000) 13.034 5.137 ( 2.822 4.889 0.000) 5.645 5.160 ( 2.166 3.752 0.000) 4.332 6.285 ( -0.834 -1.445 0.000) 1.669 6.680 ( -4.744 -8.217 0.000) 9.488 7.102 ( 2.431 4.210 0.000) 4.862 7.328 ( 1.642 2.844 0.000) 3.284 7.696 ( -0.552 -0.956 0.000) 1.103 8.011 ( -0.863 -1.494 0.000) 1.726 8.030 ( -1.470 -2.546 0.000) 2.940 8.521 ( -0.903 -1.564 0.000) 1.806 9.797 ( -1.235 -2.139 0.000) 2.470 10.007 ( -3.187 -5.520 0.000) 6.374 12.173 ( -4.181 -7.241 0.000) 8.362 12.258 ( -5.730 -9.925 0.000) 11.461 ======================= Grid point 50 (22/73) ======================= q-point: ( 0.27 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.203 ( 2.563 8.328 -0.777) 8.748 2.306 ( 7.585 8.878 -0.151) 11.678 2.436 ( 4.227 10.389 0.803) 11.244 2.675 ( -3.692 12.218 -0.375) 12.769 2.802 ( 0.880 4.788 0.843) 4.941 2.940 ( 2.674 2.453 -0.821) 3.720 3.357 ( -0.809 8.128 0.639) 8.193 3.593 ( -3.592 15.768 -0.430) 16.178 3.652 ( 5.703 14.735 0.423) 15.806 3.729 ( -0.107 3.125 0.617) 3.187 3.754 ( 0.083 -3.725 -0.900) 3.833 3.955 ( 0.380 -1.130 0.029) 1.192 4.476 ( 4.834 6.579 0.093) 8.165 5.232 ( 4.166 2.683 -0.829) 5.024 5.278 ( 7.712 0.302 0.741) 7.753 6.290 ( 2.069 -0.162 -0.008) 2.075 6.552 ( -2.140 -5.353 -0.000) 5.765 7.172 ( 0.560 4.549 0.010) 4.583 7.384 ( -0.761 6.238 -0.034) 6.284 7.667 ( -1.666 -0.485 0.035) 1.736 7.949 ( -4.646 -2.028 -0.229) 5.074 7.961 ( -2.981 -1.695 0.192) 3.434 8.505 ( 3.868 -6.472 0.039) 7.540 9.751 ( 0.950 -5.856 0.748) 5.980 9.906 ( -0.670 -8.248 -0.802) 8.314 12.067 ( -0.990 -5.386 0.350) 5.487 12.108 ( -1.673 -7.393 -0.303) 7.586 ======================= Grid point 51 (23/73) ======================= q-point: ( 0.33 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.327 ( 0.819 9.285 -0.655) 9.344 2.469 ( 0.637 7.236 -0.838) 7.312 2.524 ( -2.162 6.419 0.700) 6.809 2.755 ( 2.285 8.604 -0.026) 8.903 2.891 ( 3.076 4.576 0.956) 5.596 2.984 ( -4.074 8.633 -0.349) 9.552 3.420 ( -0.719 6.840 0.511) 6.896 3.531 ( -8.123 0.940 -0.019) 8.177 3.620 ( -5.772 -2.314 -0.265) 6.225 3.780 ( -5.189 14.400 -0.242) 15.308 3.917 ( 1.663 13.939 0.764) 14.059 3.992 ( -2.400 8.690 -0.577) 9.033 4.570 ( 2.488 0.139 0.020) 2.492 5.323 ( 5.430 -1.941 -0.487) 5.787 5.399 ( 7.660 -3.760 0.397) 8.542 6.341 ( 3.975 -0.915 0.001) 4.079 6.489 ( 0.166 -4.768 -0.012) 4.771 7.211 ( -1.063 4.256 0.014) 4.387 7.420 ( -3.318 7.996 -0.087) 8.658 7.640 ( -0.102 -1.614 0.138) 1.623 7.855 ( -3.352 -1.607 0.007) 3.717 7.894 ( -2.623 -1.147 -0.042) 2.863 8.518 ( 6.792 -10.087 0.014) 12.161 9.720 ( 3.115 -7.276 0.449) 7.927 9.823 ( 1.109 -8.582 -0.505) 8.668 12.019 ( 1.165 -5.011 0.033) 5.145 12.047 ( 1.765 -5.864 -0.004) 6.123 ======================= Grid point 52 (24/73) ======================= q-point: ( 0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.383 ( -4.662 8.074 0.000) 9.323 2.513 ( -2.936 5.086 0.000) 5.873 2.524 ( -3.748 6.492 -0.000) 7.496 2.880 ( -1.180 2.043 -0.000) 2.359 2.882 ( -0.762 1.321 0.000) 1.525 3.023 ( -7.222 12.510 0.000) 14.445 3.432 ( -2.243 3.885 0.000) 4.486 3.473 ( -1.383 2.395 -0.000) 2.765 3.535 ( -0.676 1.171 -0.000) 1.352 3.806 ( -9.306 16.119 0.000) 18.612 4.001 ( -6.972 12.076 -0.000) 13.944 4.020 ( -6.816 11.805 -0.000) 13.631 4.587 ( 2.158 -3.738 0.000) 4.316 5.363 ( 3.011 -5.216 0.000) 6.023 5.449 ( 4.240 -7.344 -0.000) 8.480 6.382 ( 1.744 -3.020 -0.000) 3.487 6.460 ( 2.430 -4.209 0.000) 4.860 7.223 ( -2.153 3.728 -0.000) 4.305 7.430 ( -4.905 8.495 -0.000) 9.810 7.632 ( 1.364 -2.363 0.000) 2.728 7.817 ( 0.085 -0.148 -0.000) 0.171 7.863 ( -0.133 0.230 -0.000) 0.265 8.528 ( 6.795 -11.768 0.000) 13.589 9.712 ( 4.172 -7.226 0.000) 8.344 9.787 ( 4.129 -7.151 0.000) 8.258 12.005 ( 2.708 -4.690 -0.000) 5.416 12.036 ( 3.043 -5.271 0.000) 6.086 ======================= Grid point 66 (25/73) ======================= q-point: ( 0.27 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.368 ( 4.759 8.243 0.000) 9.518 2.490 ( 4.988 8.640 0.000) 9.977 2.620 ( 4.023 6.968 0.000) 8.045 2.681 ( -0.928 -1.608 0.000) 1.857 2.927 ( 0.526 0.911 0.000) 1.052 3.095 ( 6.314 10.936 0.000) 12.628 3.456 ( 1.915 3.317 0.000) 3.830 3.615 ( -2.746 -4.757 0.000) 5.493 3.654 ( -3.309 -5.731 0.000) 6.618 3.907 ( 5.764 9.983 0.000) 11.527 3.973 ( 8.030 13.908 0.000) 16.060 4.072 ( 5.607 9.711 0.000) 11.213 4.553 ( 0.352 0.609 0.000) 0.703 5.268 ( 0.487 0.843 0.000) 0.973 5.292 ( 0.718 1.244 0.000) 1.436 6.298 ( 0.567 0.983 0.000) 1.135 6.461 ( -2.194 -3.801 0.000) 4.389 7.244 ( 1.499 2.596 0.000) 2.998 7.482 ( 2.241 3.881 0.000) 4.481 7.654 ( -0.661 -1.145 0.000) 1.322 7.879 ( -2.469 -4.277 0.000) 4.939 7.908 ( -1.875 -3.247 0.000) 3.749 8.391 ( -2.562 -4.437 0.000) 5.123 9.655 ( -2.017 -3.493 0.000) 4.033 9.765 ( -3.355 -5.811 0.000) 6.710 11.988 ( -1.607 -2.784 0.000) 3.214 11.995 ( -2.246 -3.890 0.000) 4.492 ======================= Grid point 67 (26/73) ======================= q-point: ( 0.33 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.517 ( 1.294 9.530 -0.088) 9.618 2.605 ( -1.055 6.299 -0.369) 6.397 2.621 ( -3.830 3.718 -0.078) 5.339 2.778 ( 7.638 -1.630 0.826) 7.853 2.884 ( -1.748 -4.950 -0.295) 5.258 3.285 ( 0.558 14.902 0.077) 14.912 3.497 ( -3.145 -0.621 -0.404) 3.231 3.508 ( -2.359 1.218 0.026) 2.654 3.556 ( -2.389 -2.966 0.333) 3.823 4.069 ( -1.738 13.444 -0.149) 13.556 4.190 ( 1.856 11.541 0.620) 11.705 4.219 ( -1.720 10.946 -0.500) 11.092 4.508 ( -0.547 -6.538 -0.002) 6.561 5.276 ( 1.534 -2.322 -0.041) 2.783 5.314 ( 3.263 -4.283 0.027) 5.384 6.327 ( 1.970 -0.261 0.002) 1.987 6.400 ( -0.372 -3.943 -0.003) 3.961 7.276 ( -0.428 2.186 0.003) 2.227 7.548 ( 0.538 4.772 -0.091) 4.803 7.617 ( -2.283 -0.683 0.096) 2.385 7.809 ( -0.990 -2.367 0.142) 2.569 7.855 ( -0.482 -2.314 -0.151) 2.368 8.336 ( 2.730 -7.213 0.016) 7.713 9.615 ( 0.997 -3.003 -0.048) 3.164 9.679 ( 0.320 -5.704 0.034) 5.713 11.944 ( 0.080 -2.498 -0.100) 2.501 11.957 ( 0.951 -3.145 0.102) 3.287 ======================= Grid point 82 (27/73) ======================= q-point: ( 0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 73 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.655 ( -0.830 -1.437 0.000) 1.659 2.655 ( 0.830 1.437 0.000) 1.659 2.726 ( -2.605 -4.512 0.000) 5.210 2.726 ( 2.605 4.512 0.000) 5.210 2.771 ( 0.000 0.000 0.000) 0.000 3.479 ( -0.991 -1.716 0.000) 1.982 3.479 ( 0.991 1.716 0.000) 1.982 3.483 ( -0.000 -0.000 0.000) 0.000 3.539 ( -0.000 -0.000 0.000) 0.000 4.316 ( -0.000 -0.000 0.000) 0.000 4.320 ( -5.589 -9.681 0.000) 11.178 4.320 ( 5.589 9.681 0.000) 11.178 4.339 ( -0.000 -0.000 0.000) 0.000 5.256 ( -0.464 -0.804 0.000) 0.928 5.256 ( 0.464 0.804 0.000) 0.928 6.341 ( -0.912 -1.579 0.000) 1.823 6.341 ( 0.912 1.579 0.000) 1.823 7.298 ( -0.000 -0.000 0.000) 0.000 7.605 ( -0.537 -0.930 0.000) 1.074 7.605 ( 0.537 0.930 0.000) 1.074 7.783 ( -0.000 -0.000 0.000) 0.000 7.826 ( -0.000 -0.000 0.000) 0.000 8.260 ( -0.000 -0.000 0.000) 0.000 9.604 ( -0.865 -1.498 0.000) 1.730 9.604 ( 0.865 1.498 0.000) 1.730 11.922 ( -0.169 -0.292 0.000) 0.337 11.922 ( 0.169 0.292 0.000) 0.337 ======================= Grid point 241 (28/73) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 73 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.288 ( 0.000 -0.000 10.054) 10.054 0.288 ( -0.000 0.000 10.054) 10.054 0.376 ( 0.000 0.000 11.734) 11.734 0.892 ( 0.000 -0.000 -6.792) 6.792 0.892 ( 0.000 0.000 -6.792) 6.792 1.321 ( 0.000 -0.000 2.761) 2.761 1.321 ( 0.000 -0.000 2.761) 2.761 1.908 ( 0.000 -0.000 -4.374) 4.374 1.955 ( 0.000 -0.000 -0.251) 0.251 1.955 ( 0.000 0.000 -0.251) 0.251 2.878 ( 0.000 -0.000 0.003) 0.003 2.878 ( 0.000 0.000 0.003) 0.003 3.765 ( 0.000 0.000 1.517) 1.517 5.447 ( 0.000 0.000 -0.369) 0.369 5.728 ( 0.000 -0.000 0.044) 0.044 6.804 ( 0.000 0.000 0.006) 0.006 6.804 ( 0.000 0.000 0.006) 0.006 6.908 ( 0.000 -0.000 -0.006) 0.006 6.908 ( 0.000 0.000 -0.006) 0.006 8.309 ( 0.000 0.000 -0.046) 0.046 8.309 ( 0.000 0.000 -0.046) 0.046 8.313 ( 0.000 0.000 0.047) 0.047 8.313 ( 0.000 -0.000 0.047) 0.047 9.226 ( 0.000 0.000 -0.047) 0.047 9.591 ( 0.000 0.000 0.044) 0.044 13.959 ( 0.000 0.000 0.073) 0.073 13.979 ( 0.000 -0.000 -0.073) 0.073 ======================= Grid point 242 (29/73) ======================= q-point: ( 0.07 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.451 ( 10.967 6.332 5.808) 13.932 0.588 ( 16.844 9.725 1.742) 19.528 0.914 ( 33.688 19.450 2.198) 38.962 1.037 ( 12.826 7.405 -3.809) 15.292 1.382 ( 31.073 17.940 -3.499) 36.050 1.433 ( 9.323 5.383 2.691) 11.097 1.627 ( 20.174 11.648 1.624) 23.352 1.988 ( 2.697 1.557 -0.488) 3.152 2.010 ( 11.933 6.890 -2.727) 14.046 2.094 ( 11.526 6.655 0.141) 13.310 2.905 ( 2.263 1.307 0.007) 2.613 3.099 ( 14.974 8.645 -1.051) 17.322 3.787 ( 2.235 1.291 1.404) 2.938 5.365 ( -4.979 -2.875 -1.274) 5.889 5.637 ( -6.975 -4.027 1.315) 8.161 6.801 ( -0.055 -0.032 -0.001) 0.064 6.914 ( 0.428 0.247 -0.004) 0.494 6.924 ( 0.872 0.504 0.122) 1.014 7.358 ( 23.779 13.729 -6.195) 28.148 8.038 ( -14.559 -8.405 1.622) 16.889 8.253 ( -3.457 -1.996 0.276) 4.001 8.311 ( -1.563 -0.903 -0.268) 1.825 8.741 ( 29.229 16.875 -3.963) 33.983 9.404 ( -15.199 -8.775 -0.342) 17.554 9.572 ( 16.199 9.353 -13.927) 23.320 13.639 ( -19.124 -11.041 11.221) 24.770 13.800 ( -14.112 -8.148 -2.552) 16.494 ======================= Grid point 243 (30/73) ======================= q-point: ( 0.13 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.820 ( 19.131 11.045 0.672) 22.101 0.976 ( 16.234 9.373 -0.153) 18.746 1.396 ( 16.701 9.642 -1.843) 19.372 1.506 ( 15.761 9.100 -0.972) 18.226 1.710 ( 13.906 8.029 2.324) 16.225 2.012 ( 21.189 12.233 -0.224) 24.468 2.080 ( 5.402 3.119 -0.592) 6.266 2.129 ( 24.371 14.071 1.065) 28.162 2.433 ( 23.305 13.455 -3.329) 27.115 2.581 ( 29.698 17.146 1.997) 34.351 2.975 ( 3.633 2.098 -0.015) 4.195 3.469 ( 15.810 9.128 -0.296) 18.258 3.884 ( 6.682 3.858 0.961) 7.775 5.239 ( -5.817 -3.358 -0.999) 6.791 5.437 ( -9.651 -5.572 0.776) 11.171 6.812 ( 1.146 0.662 -0.004) 1.323 6.887 ( -5.259 -3.036 0.137) 6.074 6.923 ( 0.308 0.178 -0.011) 0.355 7.693 ( 2.779 1.605 0.904) 3.334 7.796 ( -2.414 -1.394 -1.667) 3.247 8.160 ( -4.448 -2.568 0.431) 5.154 8.237 ( -4.901 -2.829 -0.412) 5.674 8.926 ( -21.458 -12.389 1.154) 24.805 9.390 ( 22.794 13.160 2.352) 26.425 9.846 ( 8.947 5.165 -10.196) 14.515 13.218 ( -17.498 -10.103 7.262) 21.471 13.403 ( -19.191 -11.080 -3.581) 22.447 ======================= Grid point 244 (31/73) ======================= q-point: ( 0.20 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.278 ( 20.200 11.662 -1.632) 23.382 1.307 ( 12.199 7.043 0.313) 14.090 1.739 ( 12.259 7.078 -1.793) 14.269 1.757 ( 8.840 5.104 -0.106) 10.208 2.039 ( 14.231 8.216 2.195) 16.578 2.246 ( 9.204 5.314 -0.706) 10.652 2.308 ( 7.334 4.234 0.081) 8.469 2.744 ( 26.418 15.252 2.076) 30.575 3.012 ( 22.838 13.185 -2.282) 26.469 3.064 ( 3.997 2.308 -0.061) 4.616 3.242 ( 25.653 14.811 1.102) 29.642 3.775 ( 11.035 6.371 -0.587) 12.755 4.097 ( 10.603 6.121 0.276) 12.246 5.104 ( -5.915 -3.415 0.395) 6.841 5.249 ( -4.860 -2.806 -0.303) 5.621 6.684 ( -11.051 -6.380 0.102) 12.761 6.855 ( 2.469 1.425 0.013) 2.851 6.930 ( 0.322 0.186 -0.021) 0.373 7.580 ( -11.554 -6.671 -0.546) 13.353 7.741 ( -5.950 -3.435 0.508) 6.890 8.055 ( -4.618 -2.666 0.279) 5.340 8.096 ( -6.957 -4.017 -0.269) 8.038 8.580 ( -4.684 -2.704 -0.076) 5.409 9.811 ( 13.746 7.936 2.295) 16.037 10.023 ( 6.710 3.874 -3.632) 8.557 12.843 ( -14.813 -8.552 2.946) 17.357 12.967 ( -17.914 -10.343 -2.164) 20.798 ======================= Grid point 245 (32/73) ======================= q-point: ( 0.27 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.524 ( 6.652 3.840 0.388) 7.690 1.723 ( 17.790 10.271 -1.556) 20.600 1.949 ( 6.523 3.766 -0.999) 7.598 1.963 ( 9.668 5.582 0.662) 11.184 2.340 ( 11.522 6.652 0.991) 13.342 2.439 ( 5.171 2.986 0.177) 5.974 2.483 ( 10.542 6.086 -0.277) 12.176 3.154 ( 3.797 2.192 -0.094) 4.385 3.167 ( 6.655 3.842 0.820) 7.728 3.445 ( 17.135 9.893 1.547) 19.846 3.690 ( 12.430 7.177 -1.278) 14.410 4.033 ( 11.683 6.745 -1.207) 13.545 4.275 ( 3.231 1.866 0.008) 3.731 5.022 ( 1.109 0.640 0.466) 1.363 5.253 ( 4.525 2.613 0.018) 5.226 6.456 ( -7.558 -4.364 0.020) 8.728 6.916 ( 2.512 1.451 0.021) 2.901 6.942 ( 0.992 0.573 -0.016) 1.146 7.128 ( -21.194 -12.236 -0.224) 24.474 7.699 ( 0.293 0.169 0.527) 0.627 7.934 ( -6.387 -3.688 0.071) 7.376 7.956 ( -4.247 -2.452 -0.073) 4.904 8.624 ( 5.405 3.120 -0.344) 6.250 10.030 ( 5.466 3.156 -4.368) 7.676 10.154 ( 4.698 2.712 4.399) 6.984 12.553 ( -10.253 -5.919 -0.404) 11.846 12.612 ( -12.478 -7.204 0.316) 14.412 ======================= Grid point 246 (33/73) ======================= q-point: ( 0.33 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.621 ( 2.205 1.273 0.041) 2.547 2.039 ( 9.520 5.497 2.259) 11.223 2.061 ( 3.487 2.013 0.030) 4.027 2.235 ( 13.594 7.848 -1.975) 15.821 2.555 ( 7.139 4.122 -1.182) 8.328 2.576 ( 6.871 3.967 0.233) 7.937 2.694 ( 7.356 4.247 1.274) 8.589 3.138 ( -5.250 -3.031 0.403) 6.076 3.235 ( 3.135 1.810 -0.188) 3.624 3.686 ( 0.506 0.292 0.242) 0.632 3.868 ( 4.979 2.875 -1.652) 5.982 4.178 ( -3.241 -1.871 -0.070) 3.743 4.334 ( 7.381 4.261 0.463) 8.536 5.173 ( 10.195 5.886 0.288) 11.776 5.413 ( 8.582 4.955 0.108) 9.910 6.365 ( -0.340 -0.197 -0.013) 0.393 6.764 ( -9.707 -5.604 -0.003) 11.209 6.955 ( 1.321 0.762 -0.017) 1.525 6.984 ( 2.064 1.192 0.013) 2.383 7.714 ( 0.838 0.484 0.134) 0.977 7.811 ( -4.073 -2.352 0.512) 4.731 7.871 ( -3.184 -1.839 -0.497) 3.710 8.763 ( 6.217 3.589 0.287) 7.184 10.076 ( -0.796 -0.459 -8.965) 9.012 10.236 ( 2.212 1.277 8.818) 9.180 12.371 ( -5.707 -3.295 -3.168) 7.312 12.404 ( -5.734 -3.311 2.852) 7.209 ======================= Grid point 247 (34/73) ======================= q-point: ( 0.40 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.644 ( 0.138 0.080 -0.110) 0.194 2.117 ( 1.511 0.872 0.388) 1.788 2.212 ( 6.169 3.562 2.850) 7.672 2.511 ( 9.284 5.360 -1.629) 10.843 2.673 ( 3.343 1.930 -2.159) 4.423 2.702 ( 2.043 1.180 0.555) 2.424 2.816 ( 3.421 1.975 2.128) 4.487 3.080 ( 3.070 1.773 -1.074) 3.704 3.293 ( 1.834 1.059 -0.280) 2.136 3.580 ( -6.875 -3.969 -2.068) 8.204 3.856 ( -5.065 -2.924 1.288) 5.989 4.031 ( -8.029 -4.636 -0.378) 9.279 4.533 ( 8.291 4.787 0.590) 9.591 5.408 ( 9.135 5.274 0.313) 10.553 5.612 ( 7.908 4.565 -0.016) 9.131 6.416 ( 3.876 2.238 -0.002) 4.476 6.650 ( -1.628 -0.940 0.016) 1.880 6.988 ( 1.514 0.874 -0.019) 1.749 7.021 ( 1.086 0.627 -0.002) 1.254 7.734 ( 0.853 0.492 -0.101) 0.990 7.748 ( -1.524 -0.880 0.662) 1.881 7.811 ( -2.075 -1.198 -0.630) 2.477 8.896 ( 5.016 2.896 0.666) 5.830 10.029 ( -2.513 -1.451 -7.322) 7.877 10.253 ( -0.518 -0.299 6.988) 7.013 12.285 ( -2.060 -1.189 -2.050) 3.140 12.327 ( -1.611 -0.930 1.704) 2.523 ======================= Grid point 248 (35/73) ======================= q-point: ( 0.47 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.642 ( -0.141 -0.082 -0.113) 0.198 2.137 ( 0.383 0.221 0.340) 0.558 2.325 ( 3.366 1.943 1.610) 4.207 2.616 ( 0.427 0.247 0.128) 0.510 2.716 ( 0.555 0.320 -1.344) 1.489 2.728 ( 1.616 0.933 -1.208) 2.223 2.866 ( 1.248 0.721 1.246) 1.905 3.225 ( 5.529 3.192 -0.711) 6.424 3.319 ( 0.455 0.262 -0.142) 0.544 3.449 ( -3.444 -1.988 -0.937) 4.086 3.740 ( -3.104 -1.792 1.393) 3.845 3.848 ( -6.370 -3.677 -0.761) 7.394 4.671 ( 3.275 1.891 0.213) 3.787 5.561 ( 3.774 2.179 0.263) 4.366 5.742 ( 2.884 1.665 0.199) 3.336 6.499 ( 2.475 1.429 -0.018) 2.858 6.636 ( -0.232 -0.134 0.007) 0.268 7.019 ( 1.064 0.614 -0.005) 1.229 7.033 ( -0.054 -0.031 0.001) 0.062 7.733 ( 0.054 0.031 0.176) 0.187 7.750 ( 0.391 0.226 -0.202) 0.495 7.775 ( -1.041 -0.601 -0.170) 1.214 8.979 ( 1.935 1.117 0.372) 2.265 9.987 ( -0.965 -0.557 -4.382) 4.521 10.230 ( -0.855 -0.494 4.359) 4.469 12.260 ( -0.484 -0.280 -1.042) 1.182 12.310 ( -0.083 -0.048 0.739) 0.745 ======================= Grid point 249 (36/73) ======================= q-point: (-0.47 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.642 ( 0.191 0.110 -0.145) 0.264 2.137 ( -0.390 -0.225 0.329) 0.557 2.351 ( -1.345 -0.777 0.219) 1.568 2.621 ( -0.936 -0.540 -0.241) 1.107 2.707 ( -1.215 -0.702 -0.678) 1.559 2.715 ( -1.732 -1.000 -0.551) 2.075 2.876 ( -0.483 -0.279 0.505) 0.752 3.211 ( -6.069 -3.504 0.114) 7.009 3.316 ( -0.605 -0.349 -0.021) 0.699 3.443 ( 2.666 1.539 -0.613) 3.139 3.741 ( 3.039 1.755 1.308) 3.745 3.848 ( 6.479 3.741 -0.807) 7.525 4.673 ( -3.189 -1.841 0.132) 3.685 5.571 ( -2.905 -1.677 -0.398) 3.378 5.732 ( -3.665 -2.116 0.801) 4.308 6.499 ( -2.481 -1.432 -0.011) 2.865 6.636 ( 0.263 0.152 -0.022) 0.304 7.019 ( -1.072 -0.619 0.007) 1.237 7.033 ( 0.047 0.027 0.000) 0.055 7.740 ( 0.366 0.212 -0.259) 0.496 7.749 ( -0.418 -0.241 -0.183) 0.516 7.768 ( 0.638 0.368 0.250) 0.778 8.981 ( -1.643 -0.949 0.241) 1.913 9.973 ( -0.628 -0.362 -3.158) 3.241 10.244 ( 2.319 1.339 3.163) 4.144 12.252 ( -0.201 -0.116 -0.464) 0.519 12.317 ( 0.660 0.381 0.230) 0.796 ======================= Grid point 250 (37/73) ======================= q-point: (-0.40 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.646 ( -0.048 -0.028 -0.233) 0.239 2.114 ( -1.766 -1.020 0.538) 2.109 2.267 ( -5.582 -3.223 0.048) 6.445 2.468 ( -10.985 -6.342 0.410) 12.691 2.662 ( -2.739 -1.582 -1.289) 3.416 2.709 ( -2.380 -1.374 0.151) 2.753 2.829 ( -3.927 -2.268 1.067) 4.659 3.091 ( -0.709 -0.409 -1.761) 1.942 3.290 ( -1.770 -1.022 -0.109) 2.047 3.546 ( 5.594 3.230 0.070) 6.460 3.868 ( 6.039 3.487 0.555) 6.995 4.028 ( 7.657 4.421 -0.200) 8.844 4.532 ( -8.657 -4.998 0.623) 10.016 5.433 ( -8.752 -5.053 -1.308) 10.190 5.590 ( -8.164 -4.713 1.412) 9.532 6.416 ( -3.867 -2.233 0.003) 4.465 6.650 ( 1.574 0.909 -0.001) 1.817 6.988 ( -1.509 -0.871 -0.006) 1.742 7.021 ( -1.087 -0.628 0.000) 1.256 7.734 ( -0.841 -0.485 -0.065) 0.973 7.757 ( 1.447 0.835 0.270) 1.692 7.802 ( 2.150 1.241 -0.262) 2.496 8.910 ( -4.416 -2.550 -0.130) 5.101 9.957 ( -0.592 -0.342 -1.258) 1.432 10.317 ( 3.213 1.855 1.376) 3.957 12.258 ( 1.059 0.611 0.121) 1.228 12.349 ( 2.462 1.421 -0.228) 2.852 ======================= Grid point 251 (38/73) ======================= q-point: (-0.33 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.626 ( -2.111 -1.219 -0.230) 2.448 2.050 ( -3.875 -2.237 0.648) 4.521 2.103 ( -9.387 -5.419 -1.199) 10.905 2.164 ( -14.309 -8.261 1.902) 16.632 2.569 ( -5.704 -3.293 -2.279) 6.969 2.575 ( -7.072 -4.083 0.311) 8.172 2.686 ( -8.564 -4.945 1.972) 10.084 3.186 ( 5.806 3.352 -2.447) 7.136 3.234 ( -3.062 -1.768 -0.133) 3.538 3.652 ( 2.012 1.161 2.057) 3.103 3.881 ( -7.123 -4.112 -2.183) 8.510 4.170 ( 2.483 1.434 0.483) 2.908 4.323 ( -7.141 -4.123 1.124) 8.322 5.203 ( -10.040 -5.797 -1.710) 11.719 5.390 ( -8.516 -4.917 1.675) 9.975 6.365 ( 0.367 0.212 -0.032) 0.425 6.761 ( 9.531 5.503 0.131) 11.006 6.955 ( -1.306 -0.754 -0.016) 1.508 6.984 ( -2.074 -1.197 0.023) 2.395 7.715 ( -0.721 -0.416 0.086) 0.837 7.818 ( 3.948 2.279 0.196) 4.563 7.864 ( 3.302 1.906 -0.197) 3.817 8.786 ( -6.128 -3.538 -1.051) 7.153 9.939 ( -1.450 -0.837 2.212) 2.774 10.360 ( 0.082 0.047 -1.589) 1.592 12.321 ( 4.788 2.765 1.009) 5.620 12.447 ( 6.489 3.746 -0.868) 7.543 ======================= Grid point 252 (39/73) ======================= q-point: (-0.27 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.526 ( -7.054 -4.073 0.254) 8.150 1.744 ( -18.796 -10.852 -3.345) 21.960 1.929 ( -8.619 -4.976 3.247) 10.469 1.935 ( -6.215 -3.588 -0.169) 7.179 2.382 ( -11.102 -6.410 -2.030) 12.979 2.433 ( -5.441 -3.141 0.536) 6.305 2.452 ( -10.920 -6.305 2.052) 12.775 3.154 ( -3.784 -2.185 -0.085) 4.370 3.204 ( -8.022 -4.632 -1.334) 9.359 3.510 ( -16.232 -9.371 -2.572) 18.919 3.630 ( -11.575 -6.683 2.589) 13.614 4.005 ( -11.559 -6.673 0.381) 13.352 4.273 ( -3.395 -1.960 0.086) 3.921 5.052 ( -1.481 -0.855 -1.421) 2.224 5.233 ( -4.314 -2.491 1.345) 5.160 6.457 ( 7.571 4.371 -0.035) 8.742 6.916 ( -2.515 -1.452 -0.005) 2.904 6.942 ( -0.983 -0.567 0.010) 1.135 7.121 ( 20.993 12.120 0.197) 24.242 7.703 ( -0.162 -0.094 0.273) 0.332 7.936 ( 6.064 3.501 -0.073) 7.002 7.954 ( 4.562 2.634 0.069) 5.268 8.644 ( -5.687 -3.283 -1.503) 6.736 9.865 ( -5.604 -3.235 5.144) 8.267 10.312 ( -4.115 -2.376 -4.706) 6.688 12.493 ( 10.325 5.961 2.352) 12.152 12.661 ( 12.243 7.069 -1.755) 14.246 ======================= Grid point 253 (40/73) ======================= q-point: (-0.20 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.290 ( -13.450 -7.765 1.233) 15.579 1.291 ( -20.316 -11.730 -2.505) 23.592 1.738 ( -8.184 -4.725 1.148) 9.519 1.748 ( -10.730 -6.195 -2.154) 12.576 2.070 ( -15.093 -8.714 0.447) 17.434 2.221 ( -8.866 -5.119 0.668) 10.259 2.295 ( -7.560 -4.365 0.870) 8.773 2.768 ( -26.623 -15.371 0.677) 30.749 3.035 ( -25.262 -14.585 -4.217) 29.473 3.064 ( -3.977 -2.296 -0.074) 4.593 3.227 ( -23.879 -13.787 2.488) 27.685 3.765 ( -10.117 -5.841 -0.066) 11.682 4.093 ( -10.571 -6.103 0.473) 12.215 5.114 ( 4.689 2.707 -0.168) 5.417 5.242 ( 5.824 3.362 0.096) 6.725 6.683 ( 10.939 6.316 0.109) 12.631 6.855 ( -2.482 -1.433 -0.001) 2.866 6.929 ( -0.318 -0.183 -0.006) 0.367 7.584 ( 13.607 7.856 -0.829) 15.734 7.731 ( 3.654 2.109 1.170) 4.378 8.058 ( 4.884 2.820 0.149) 5.641 8.094 ( 6.703 3.870 -0.143) 7.741 8.595 ( 4.713 2.721 -1.028) 5.538 9.658 ( -12.689 -7.326 6.966) 16.224 10.179 ( -7.408 -4.277 -8.522) 12.074 12.791 ( 15.431 8.909 4.490) 18.375 13.003 ( 17.198 9.929 -2.803) 20.055 ======================= Grid point 254 (41/73) ======================= q-point: (-0.13 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.827 ( -19.604 -11.318 0.313) 22.639 0.924 ( -17.863 -10.313 2.556) 20.784 1.437 ( -15.655 -9.038 -4.044) 18.523 1.512 ( -14.113 -8.148 -1.403) 16.357 1.726 ( -14.288 -8.249 1.502) 16.567 2.015 ( -19.304 -11.145 -0.232) 22.291 2.061 ( -5.092 -2.940 0.482) 5.900 2.120 ( -27.259 -15.738 1.594) 31.516 2.413 ( -25.037 -14.455 -2.003) 28.980 2.604 ( -28.004 -16.168 0.375) 32.339 2.975 ( -3.633 -2.097 -0.045) 4.195 3.468 ( -15.863 -9.158 -0.246) 18.318 3.884 ( -6.569 -3.793 1.014) 7.653 5.235 ( 5.774 3.334 -0.750) 6.710 5.440 ( 9.663 5.579 0.545) 11.172 6.812 ( -1.155 -0.667 -0.000) 1.334 6.885 ( 5.293 3.056 0.263) 6.118 6.923 ( -0.300 -0.173 -0.003) 0.346 7.676 ( -7.280 -4.203 1.581) 8.554 7.803 ( 6.777 3.912 -1.796) 8.028 8.166 ( 4.457 2.573 0.114) 5.147 8.231 ( 4.893 2.825 -0.110) 5.651 8.942 ( 21.745 12.555 0.305) 25.111 9.275 ( -20.755 -11.983 6.100) 24.730 9.967 ( -11.436 -6.602 -14.490) 19.605 13.184 ( 18.450 10.652 8.193) 22.825 13.421 ( 18.407 10.627 -3.548) 21.549 ======================= Grid point 255 (42/73) ======================= q-point: (-0.07 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.447 ( -11.107 -6.413 6.106) 14.205 0.506 ( -17.305 -9.991 5.902) 20.836 0.956 ( -32.668 -18.861 -0.108) 37.722 1.085 ( -13.557 -7.827 -5.931) 16.740 1.359 ( -32.137 -18.554 -2.408) 37.186 1.442 ( -9.649 -5.571 2.278) 11.373 1.625 ( -19.497 -11.257 1.895) 22.593 1.978 ( -2.266 -1.308 0.098) 2.619 1.997 ( -10.392 -6.000 -2.165) 12.194 2.110 ( -13.128 -7.579 -0.788) 15.179 2.905 ( -2.277 -1.315 -0.011) 2.629 3.096 ( -14.943 -8.627 -0.806) 17.273 3.788 ( -2.243 -1.295 1.384) 2.937 5.364 ( 5.124 2.958 -1.182) 6.033 5.638 ( 6.863 3.962 1.219) 8.018 6.801 ( 0.060 0.034 0.004) 0.069 6.914 ( -0.424 -0.245 -0.006) 0.490 6.925 ( -0.719 -0.415 0.048) 0.832 7.328 ( -22.450 -12.962 -3.953) 26.223 8.062 ( 14.063 8.119 -0.061) 16.238 8.257 ( 3.290 1.900 0.079) 3.800 8.308 ( 1.720 0.993 -0.078) 1.988 8.703 ( -27.025 -15.603 -1.553) 31.244 9.405 ( 15.465 8.929 -0.242) 17.859 9.617 ( -19.987 -11.539 -17.158) 28.758 13.629 ( 20.063 11.583 11.666) 25.939 13.804 ( 13.673 7.894 -2.460) 15.978 ======================= Grid point 257 (43/73) ======================= q-point: ( 0.07 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.697 ( 10.206 17.951 1.826) 20.730 0.924 ( 15.141 19.704 1.990) 24.929 1.294 ( 10.519 23.182 -4.049) 25.777 1.413 ( 7.741 15.662 -1.075) 17.503 1.635 ( 9.444 17.185 1.488) 19.666 1.880 ( 12.059 22.944 0.804) 25.932 1.947 ( 11.509 17.810 0.664) 21.215 2.062 ( 7.839 10.100 -1.711) 12.899 2.319 ( 13.487 17.894 -0.168) 22.409 2.361 ( 11.261 26.257 -0.096) 28.570 2.972 ( 2.837 4.944 -0.011) 5.700 3.357 ( 9.241 15.849 -0.380) 18.350 3.843 ( 2.603 4.500 1.162) 5.327 5.272 ( -3.176 -5.810 -1.151) 6.720 5.496 ( -5.626 -9.468 0.913) 11.052 6.824 ( 0.907 1.539 0.079) 1.788 6.900 ( 0.051 0.054 -0.004) 0.074 6.913 ( -2.195 -3.863 0.263) 4.451 7.627 ( 8.217 11.255 -0.705) 13.953 7.843 ( -6.373 -8.988 -0.536) 11.031 8.229 ( -3.545 -2.554 0.108) 4.370 8.247 ( -0.204 -3.823 0.546) 3.867 9.084 ( -12.445 -18.360 -0.091) 22.181 9.176 ( 17.176 19.005 4.142) 25.949 9.827 ( 2.095 12.869 -14.856) 19.766 13.312 ( -9.682 -19.150 9.187) 23.342 13.525 ( -11.107 -18.004 -3.652) 21.467 ======================= Grid point 258 (44/73) ======================= q-point: ( 0.13 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.101 ( 15.473 16.846 -1.021) 22.897 1.314 ( 6.907 20.904 0.195) 22.017 1.635 ( 7.995 5.837 -0.442) 9.909 1.732 ( 10.420 21.436 -2.093) 23.926 2.004 ( 9.894 20.300 2.066) 22.677 2.124 ( 6.625 6.646 -0.117) 9.385 2.231 ( 10.843 5.647 0.247) 12.228 2.536 ( 18.181 22.705 0.134) 29.087 2.784 ( 17.495 19.969 -0.537) 26.554 2.957 ( 24.837 19.429 0.549) 31.538 3.085 ( 2.619 8.791 -0.433) 9.183 3.669 ( 10.157 10.470 -0.259) 14.590 3.978 ( 7.727 5.082 0.640) 9.271 5.157 ( -4.125 -4.184 -0.083) 5.876 5.290 ( -6.319 -8.664 0.025) 10.724 6.746 ( -6.981 -10.387 0.182) 12.516 6.864 ( 0.072 3.480 0.011) 3.481 6.917 ( 1.395 -0.062 -0.024) 1.396 7.695 ( -4.868 0.562 -0.149) 4.902 7.767 ( -1.033 -4.010 -0.598) 4.184 8.133 ( -5.163 -0.984 0.440) 5.275 8.157 ( -4.928 -5.623 0.152) 7.478 8.675 ( -11.477 -10.412 0.116) 15.497 9.623 ( 15.869 9.542 5.510) 19.319 9.963 ( 3.702 5.571 -8.496) 10.813 12.953 ( -11.428 -15.428 4.917) 19.819 13.118 ( -14.811 -16.605 -3.169) 22.475 ======================= Grid point 259 (45/73) ======================= q-point: ( 0.20 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.539 ( 16.669 14.241 -2.076) 22.022 1.586 ( 1.940 18.732 0.407) 18.837 1.828 ( 8.278 5.033 0.372) 9.695 1.999 ( 1.077 13.250 -0.461) 13.302 2.318 ( 7.209 12.863 -0.583) 14.757 2.362 ( 7.731 14.179 1.262) 16.199 2.480 ( 7.425 13.612 -0.023) 15.505 3.007 ( 15.015 9.038 0.872) 17.547 3.128 ( 6.638 5.079 0.693) 8.387 3.317 ( 10.639 15.520 -0.380) 18.820 3.502 ( 18.812 10.934 -0.312) 21.761 3.910 ( 9.365 7.186 -0.831) 11.834 4.174 ( 8.954 1.205 0.126) 9.035 5.043 ( -4.007 -2.192 0.528) 4.598 5.190 ( 2.444 -2.734 -0.268) 3.677 6.511 ( -6.935 -9.435 0.034) 11.709 6.911 ( 1.049 2.805 -0.002) 2.995 6.933 ( 0.046 0.610 -0.010) 0.612 7.423 ( -21.659 -8.627 -0.364) 23.317 7.690 ( -0.082 -0.827 0.504) 0.972 8.028 ( -5.657 -1.773 -0.271) 5.935 8.038 ( -5.623 0.049 0.358) 5.635 8.567 ( 2.069 0.373 -0.364) 2.133 9.904 ( 10.588 1.435 2.072) 10.884 10.072 ( 5.250 0.648 -2.426) 5.820 12.634 ( -9.323 -11.867 1.948) 15.217 12.725 ( -12.564 -13.371 -1.641) 18.421 ======================= Grid point 260 (46/73) ======================= q-point: ( 0.27 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.750 ( -3.445 17.274 0.267) 17.616 1.918 ( 10.901 8.486 0.720) 13.833 2.027 ( 8.627 7.127 -1.430) 11.281 2.125 ( 4.903 7.369 0.269) 8.855 2.528 ( 5.825 8.068 -0.247) 9.954 2.605 ( 1.093 17.941 -0.533) 17.982 2.727 ( 2.524 17.100 0.881) 17.308 3.158 ( 0.095 -3.033 0.437) 3.066 3.249 ( 2.392 6.912 0.199) 7.316 3.623 ( 8.741 6.143 1.941) 10.859 3.787 ( 7.475 3.131 -2.680) 8.536 4.136 ( 6.957 1.995 -0.391) 7.248 4.260 ( 4.188 -0.961 0.189) 4.300 5.067 ( 6.579 3.532 0.267) 7.472 5.287 ( 8.785 0.830 0.138) 8.825 6.363 ( -1.928 -4.125 -0.011) 4.554 6.837 ( -6.793 -8.934 -0.001) 11.223 6.975 ( 0.573 4.221 0.019) 4.260 7.089 ( -10.229 5.800 -0.047) 11.759 7.694 ( 0.987 -0.643 0.250) 1.204 7.908 ( -7.181 0.829 0.186) 7.231 7.947 ( -5.615 1.143 -0.186) 5.734 8.660 ( 7.064 0.061 -0.150) 7.066 10.021 ( 4.574 -3.719 -4.904) 7.668 10.151 ( 5.310 -3.076 4.939) 7.877 12.408 ( -5.229 -8.297 -1.080) 9.867 12.442 ( -6.678 -9.595 0.943) 11.728 ======================= Grid point 261 (47/73) ======================= q-point: ( 0.33 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.817 ( -8.037 17.094 -0.060) 18.889 2.078 ( 4.061 2.780 1.513) 5.149 2.199 ( 3.067 7.809 1.349) 8.498 2.380 ( 12.337 6.641 -1.922) 14.142 2.689 ( 3.458 7.353 0.379) 8.135 2.751 ( -1.973 12.866 -1.805) 13.140 2.888 ( -1.691 13.621 1.563) 13.814 3.087 ( -2.678 0.685 0.018) 2.765 3.348 ( -0.696 9.251 -0.236) 9.280 3.652 ( -4.756 -2.730 -1.529) 5.693 3.869 ( 0.547 -2.766 0.440) 2.853 4.088 ( -4.468 -6.279 -0.204) 7.709 4.430 ( 8.925 5.171 0.598) 10.332 5.276 ( 10.921 4.064 0.220) 11.655 5.474 ( 10.814 1.179 0.018) 10.878 6.360 ( 3.439 -0.170 0.002) 3.444 6.665 ( -3.253 -4.308 0.010) 5.398 7.018 ( -0.443 4.394 -0.007) 4.417 7.068 ( -3.113 7.223 -0.004) 7.866 7.707 ( 1.570 -1.136 -0.011) 1.938 7.812 ( -5.267 1.847 0.540) 5.607 7.870 ( -4.627 1.275 -0.526) 4.828 8.786 ( 7.979 -1.515 0.417) 8.132 10.010 ( 0.602 -5.821 -7.208) 9.284 10.188 ( 3.931 -5.866 6.983) 9.931 12.279 ( -1.688 -5.825 -2.503) 6.561 12.306 ( -0.598 -6.353 2.242) 6.764 ======================= Grid point 262 (48/73) ======================= q-point: ( 0.40 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.831 ( -9.823 17.371 -0.174) 19.956 2.142 ( 0.522 2.591 0.511) 2.692 2.313 ( 2.425 5.564 1.776) 6.324 2.596 ( 6.206 3.683 -0.823) 7.263 2.743 ( -3.732 3.366 0.512) 5.052 2.808 ( -2.274 9.756 -1.952) 10.206 2.967 ( -4.088 12.489 1.429) 13.218 3.157 ( 4.900 6.425 -0.965) 8.138 3.403 ( -3.314 9.052 -0.158) 9.641 3.511 ( -6.426 -2.435 -1.542) 7.043 3.784 ( -4.303 -3.317 1.439) 5.620 3.922 ( -7.017 -5.602 -0.470) 8.991 4.608 ( 6.285 2.574 0.314) 6.799 5.471 ( 8.018 0.955 0.319) 8.081 5.648 ( 8.145 -0.925 0.008) 8.197 6.442 ( 4.818 0.350 -0.005) 4.830 6.618 ( 0.472 -2.494 0.009) 2.539 7.049 ( -0.901 4.068 -0.015) 4.166 7.087 ( -2.913 6.419 0.016) 7.049 7.723 ( 1.688 -1.427 -0.203) 2.220 7.764 ( -2.556 1.829 0.429) 3.172 7.816 ( -3.134 1.369 -0.385) 3.441 8.890 ( 6.594 -3.615 0.467) 7.534 9.959 ( 0.697 -5.300 -5.153) 7.425 10.175 ( 2.576 -7.175 4.994) 9.113 12.228 ( 1.161 -4.444 -1.156) 4.736 12.267 ( 1.944 -4.844 0.892) 5.295 ======================= Grid point 263 (49/73) ======================= q-point: ( 0.47 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.832 ( -10.182 17.666 -0.176) 20.391 2.159 ( -1.103 2.317 0.004) 2.566 2.378 ( -0.338 2.801 0.756) 2.920 2.675 ( -2.574 5.638 0.065) 6.198 2.707 ( -0.276 -1.005 -0.653) 1.230 2.840 ( -5.965 8.996 -0.947) 10.836 2.995 ( -6.120 11.302 0.669) 12.870 3.252 ( -1.312 0.920 -0.061) 1.603 3.422 ( -4.365 8.148 0.176) 9.245 3.432 ( -1.051 0.647 -0.515) 1.337 3.729 ( -0.324 0.895 0.990) 1.374 3.804 ( -0.014 0.039 -0.820) 0.821 4.678 ( 0.770 -1.206 0.079) 1.433 5.553 ( 2.273 -2.900 -0.002) 3.685 5.711 ( 2.095 -4.555 0.407) 5.030 6.494 ( 1.085 -1.896 -0.019) 2.184 6.609 ( 1.422 -2.408 -0.010) 2.796 7.063 ( -1.787 3.073 -0.000) 3.554 7.093 ( -3.581 6.200 0.029) 7.160 7.732 ( 1.057 -1.919 -0.202) 2.200 7.756 ( -0.686 1.989 -0.241) 2.117 7.791 ( -1.616 2.063 0.225) 2.630 8.931 ( 3.483 -5.887 0.240) 6.844 9.936 ( 1.904 -4.444 -3.265) 5.834 10.160 ( 4.257 -6.210 3.302) 8.221 12.215 ( 1.939 -4.016 -0.498) 4.488 12.264 ( 2.863 -4.403 0.277) 5.259 ======================= Grid point 265 (50/73) ======================= q-point: (-0.47 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.832 ( -10.047 17.237 -0.237) 19.953 2.150 ( -2.479 1.523 0.054) 2.910 2.348 ( -4.801 0.366 -0.031) 4.815 2.576 ( -8.019 -4.500 0.169) 9.197 2.750 ( -0.830 5.294 0.178) 5.362 2.795 ( -6.938 7.199 -1.120) 10.060 2.972 ( -8.344 8.320 0.956) 11.822 3.153 ( -7.319 0.639 -0.734) 7.383 3.401 ( -6.391 7.452 0.001) 9.817 3.494 ( 4.095 3.711 -0.475) 5.547 3.787 ( 5.501 2.353 1.209) 6.104 3.923 ( 8.127 3.231 -0.527) 8.762 4.609 ( -5.427 -4.240 0.242) 6.891 5.489 ( -4.174 -6.261 -0.817) 7.569 5.632 ( -3.843 -7.652 1.014) 8.623 6.442 ( -2.708 -3.991 0.002) 4.823 6.619 ( 1.920 -1.665 -0.026) 2.542 7.049 ( -3.070 2.818 -0.002) 4.167 7.087 ( -4.111 5.727 0.023) 7.049 7.722 ( 0.394 -2.196 -0.151) 2.236 7.774 ( -0.180 3.126 -0.075) 3.132 7.807 ( 0.257 3.393 0.092) 3.404 8.896 ( 0.326 -7.419 0.143) 7.427 9.925 ( 1.760 -4.178 -2.181) 5.031 10.207 ( 7.091 -4.944 2.150) 8.908 12.213 ( 2.394 -3.558 -0.032) 4.288 12.280 ( 3.968 -3.807 -0.104) 5.500 ======================= Grid point 266 (51/73) ======================= q-point: (-0.40 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.821 ( -10.706 15.743 -0.302) 19.041 2.094 ( -4.310 -1.122 0.468) 4.478 2.235 ( -8.107 0.022 -0.007) 8.107 2.327 ( -13.358 -7.507 0.513) 15.331 2.688 ( -8.120 -0.166 0.489) 8.136 2.746 ( -10.179 8.669 -1.462) 13.450 2.892 ( -10.769 7.696 1.217) 13.292 3.121 ( 2.104 3.390 -2.001) 4.464 3.345 ( -7.480 5.181 -0.064) 9.100 3.611 ( 4.266 2.980 0.891) 5.280 3.890 ( 2.429 -1.654 -0.799) 3.045 4.082 ( 7.521 0.698 0.150) 7.555 4.424 ( -9.486 -5.345 0.935) 10.928 5.302 ( -8.655 -7.534 -1.498) 11.572 5.453 ( -6.601 -8.628 1.463) 10.962 6.360 ( -1.559 -3.044 -0.009) 3.420 6.664 ( 5.265 0.632 0.036) 5.303 7.018 ( -3.564 2.583 -0.011) 4.402 7.068 ( -4.726 6.290 0.025) 7.867 7.707 ( 0.287 -1.848 -0.013) 1.870 7.821 ( 0.922 5.413 0.155) 5.493 7.862 ( 1.316 4.701 -0.154) 4.885 8.801 ( -2.169 -8.002 -0.440) 8.302 9.917 ( 1.527 -3.442 0.467) 3.794 10.273 ( 6.002 -6.080 -0.213) 8.546 12.242 ( 4.707 -1.923 0.659) 5.127 12.338 ( 6.802 -2.175 -0.625) 7.169 ======================= Grid point 267 (52/73) ======================= q-point: (-0.33 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.756 ( -13.463 11.825 -0.071) 17.919 1.961 ( -15.243 -6.976 -3.105) 17.049 1.985 ( -9.327 -4.542 2.189) 10.603 2.106 ( -8.188 2.241 1.431) 8.609 2.516 ( -9.958 -1.465 0.476) 10.076 2.623 ( -14.128 7.419 -1.906) 16.071 2.717 ( -17.473 6.608 1.810) 18.768 3.195 ( 0.912 -2.464 -1.726) 3.144 3.257 ( -5.604 1.893 -0.361) 5.926 3.634 ( -4.333 -2.640 1.994) 5.452 3.762 ( -11.119 -6.456 -1.869) 12.992 4.117 ( -6.391 -5.071 0.757) 8.193 4.258 ( -0.573 -3.783 0.292) 3.837 5.096 ( -6.187 -4.200 -1.632) 7.654 5.266 ( -5.101 -7.038 1.573) 8.833 6.363 ( 4.558 -0.386 -0.049) 4.575 6.835 ( 11.126 1.501 0.125) 11.228 6.975 ( -3.930 1.603 0.002) 4.244 7.086 ( -0.152 11.584 0.106) 11.585 7.697 ( 0.219 -0.965 0.029) 0.990 7.914 ( 2.740 6.780 -0.186) 7.315 7.942 ( 2.005 5.445 0.132) 5.804 8.673 ( -3.491 -7.094 -0.933) 7.961 9.892 ( -1.167 -3.366 3.527) 5.013 10.274 ( 2.173 -7.999 -3.125) 8.858 12.355 ( 9.391 0.966 1.718) 9.595 12.487 ( 11.849 0.509 -1.392) 11.942 ======================= Grid point 268 (53/73) ======================= q-point: (-0.27 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.550 ( -20.850 -7.726 -2.877) 22.421 1.582 ( -17.800 7.459 0.485) 19.305 1.809 ( -8.235 -4.121 1.797) 9.382 1.997 ( -11.811 6.841 -0.366) 13.654 2.317 ( -9.935 -3.992 -0.178) 10.708 2.382 ( -20.945 4.211 -0.472) 21.369 2.457 ( -15.279 0.335 1.524) 15.359 3.019 ( -15.110 -9.740 0.225) 17.979 3.151 ( -7.828 -1.486 -0.711) 7.999 3.351 ( -20.886 -4.569 -3.144) 21.610 3.464 ( -16.965 -9.517 2.662) 19.633 3.894 ( -10.403 -3.640 0.114) 11.022 4.173 ( -5.658 -7.106 0.237) 9.086 5.066 ( 3.243 1.980 -0.879) 3.900 5.172 ( 1.609 -3.368 0.838) 3.825 6.512 ( 11.620 1.261 -0.008) 11.688 6.911 ( -2.871 0.429 0.010) 2.903 6.934 ( -0.606 0.344 -0.031) 0.697 7.412 ( 18.062 14.141 0.293) 22.941 7.698 ( 0.723 -0.056 0.006) 0.725 8.028 ( 3.109 4.148 -0.125) 5.185 8.044 ( 4.104 5.421 -0.166) 6.801 8.569 ( -1.178 -2.928 -0.496) 3.195 9.785 ( -7.661 -3.745 5.718) 10.267 10.193 ( -1.878 -8.529 -6.199) 10.709 12.589 ( 14.890 3.739 3.246) 15.692 12.761 ( 17.680 2.907 -2.356) 18.072 ======================= Grid point 269 (54/73) ======================= q-point: (-0.20 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.108 ( -22.358 -5.552 -1.403) 23.080 1.295 ( -22.793 5.526 1.183) 23.483 1.629 ( -8.490 -3.439 -0.063) 9.160 1.748 ( -22.682 1.113 -3.003) 22.907 2.015 ( -22.657 -0.143 1.532) 22.709 2.133 ( -8.495 -0.099 -0.724) 8.527 2.207 ( -10.851 -6.379 1.610) 12.690 2.532 ( -29.047 -5.813 0.454) 29.627 2.808 ( -28.317 -5.117 -2.382) 28.874 2.943 ( -26.623 -11.663 1.738) 29.117 3.087 ( -9.538 1.812 -0.542) 9.724 3.667 ( -13.927 -3.338 -0.176) 14.322 3.977 ( -8.210 -4.013 0.732) 9.168 5.158 ( 5.185 1.310 -0.124) 5.349 5.290 ( 11.077 1.242 0.033) 11.146 6.745 ( 12.373 0.888 0.220) 12.407 6.864 ( -3.018 1.655 0.011) 3.442 6.918 ( -0.661 -1.221 -0.046) 1.389 7.706 ( 2.900 4.043 -0.864) 5.050 7.752 ( 3.033 -0.269 0.419) 3.073 8.133 ( 5.055 3.502 0.424) 6.164 8.163 ( 5.111 1.793 -0.164) 5.419 8.676 ( 15.487 4.140 0.020) 16.031 9.546 ( -16.722 -2.103 7.152) 18.308 10.044 ( -6.097 -7.584 -10.465) 14.290 12.928 ( 19.148 4.451 5.557) 20.429 13.135 ( 21.556 2.980 -3.290) 22.008 ======================= Grid point 274 (55/73) ======================= q-point: ( 0.13 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.465 ( 10.679 19.117 -2.383) 22.027 1.697 ( 9.294 16.001 0.302) 18.507 1.764 ( 5.320 7.636 0.986) 9.359 2.129 ( 7.608 15.239 -0.807) 17.052 2.253 ( 4.450 4.566 0.785) 6.424 2.413 ( 10.105 19.729 -0.687) 22.177 2.474 ( 13.081 21.936 0.861) 25.555 2.911 ( 8.491 13.821 -0.460) 16.228 3.100 ( 3.222 10.026 1.879) 10.698 3.334 ( 12.176 18.737 -0.929) 22.365 3.342 ( 9.109 15.896 -0.708) 18.335 3.869 ( 5.830 9.340 -0.050) 11.010 4.076 ( 2.762 4.668 0.307) 5.433 5.101 ( -3.557 -0.973 0.976) 3.815 5.124 ( -1.121 -6.889 -0.907) 7.038 6.521 ( -6.227 -10.744 0.015) 12.418 6.918 ( -0.137 -0.318 -0.061) 0.351 6.940 ( 2.062 3.674 0.028) 4.213 7.573 ( -7.860 -15.172 0.693) 17.102 7.708 ( -0.623 0.113 -0.671) 0.922 8.056 ( -1.873 -2.660 -0.457) 3.286 8.104 ( -1.713 -2.121 -0.176) 2.732 8.558 ( -0.618 -2.194 0.705) 2.386 9.769 ( 6.850 5.355 5.334) 10.200 10.052 ( -2.122 3.358 -6.421) 7.550 12.640 ( -7.301 -15.274 3.483) 17.284 12.778 ( -10.852 -16.760 -2.673) 20.144 ======================= Grid point 275 (56/73) ======================= q-point: ( 0.20 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.828 ( 10.510 14.427 -1.266) 17.894 1.933 ( 4.776 12.903 0.432) 13.765 1.959 ( 6.629 8.401 0.127) 10.702 2.252 ( -5.742 15.212 0.560) 16.269 2.480 ( 11.268 7.128 -0.103) 13.334 2.798 ( 4.334 23.668 -1.225) 24.092 2.873 ( 5.731 23.185 0.840) 23.897 3.078 ( 5.300 -1.040 -0.236) 5.406 3.234 ( 5.163 5.494 1.563) 7.700 3.592 ( 4.997 10.536 0.710) 11.683 3.677 ( 9.705 6.632 -1.755) 11.885 4.061 ( 5.372 7.702 -0.199) 9.392 4.150 ( 5.497 -2.459 0.167) 6.024 5.037 ( 0.098 1.934 -0.104) 1.939 5.156 ( 5.768 -0.438 0.283) 5.792 6.355 ( -3.040 -5.769 -0.029) 6.521 6.851 ( -4.478 -7.274 0.028) 8.542 7.011 ( 1.146 4.781 0.023) 4.917 7.307 ( -8.860 -1.736 0.167) 9.030 7.693 ( -1.308 0.501 -0.119) 1.405 8.012 ( -3.012 -0.021 -0.143) 3.015 8.050 ( -3.816 0.370 -0.095) 3.835 8.562 ( 2.639 -1.228 0.117) 2.913 9.885 ( 5.249 -2.702 2.061) 6.253 10.050 ( 1.424 -2.581 -2.250) 3.709 12.395 ( -5.226 -11.621 1.984) 12.896 12.460 ( -8.408 -12.699 -1.790) 15.335 ======================= Grid point 276 (57/73) ======================= q-point: ( 0.27 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.042 ( 2.995 7.773 1.511) 8.466 2.114 ( 2.249 13.184 0.306) 13.377 2.240 ( 12.230 9.908 -2.204) 15.894 2.303 ( -6.156 16.355 0.826) 17.495 2.702 ( 5.024 8.218 -0.078) 9.632 2.919 ( -3.769 4.171 -0.929) 5.698 3.087 ( -2.951 16.957 -0.042) 17.212 3.201 ( -0.379 16.629 1.118) 16.671 3.418 ( 6.221 9.481 0.086) 11.340 3.689 ( -1.190 1.902 -0.045) 2.244 3.806 ( 3.351 -1.624 -0.782) 3.805 4.063 ( -1.089 -6.757 -0.111) 6.845 4.337 ( 8.814 7.085 0.372) 11.315 5.147 ( 7.642 4.121 0.055) 8.682 5.292 ( 9.859 -0.171 0.094) 9.861 6.306 ( 1.458 -1.613 0.004) 2.174 6.676 ( -4.638 -7.157 0.023) 8.528 7.072 ( 0.046 5.271 0.011) 5.271 7.237 ( -5.210 8.068 -0.015) 9.604 7.678 ( -0.145 -0.762 0.040) 0.776 7.939 ( -5.999 1.343 -0.026) 6.148 7.976 ( -5.370 1.137 -0.041) 5.489 8.619 ( 6.846 -4.162 0.026) 8.012 9.904 ( 3.047 -7.347 -2.205) 8.254 10.045 ( 3.469 -7.003 2.174) 8.112 12.237 ( -1.457 -8.507 0.379) 8.639 12.249 ( -3.340 -9.261 -0.420) 9.854 ======================= Grid point 277 (58/73) ======================= q-point: ( 0.33 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.146 ( -0.016 8.662 0.568) 8.681 2.242 ( -6.144 17.293 -0.100) 18.352 2.359 ( -1.345 9.046 1.437) 9.258 2.512 ( 7.937 7.885 -1.051) 11.237 2.834 ( 4.032 5.146 -0.564) 6.561 2.885 ( -3.130 2.366 -0.191) 3.928 3.187 ( -4.076 16.786 -0.353) 17.277 3.315 ( -6.732 17.867 0.840) 19.112 3.553 ( -3.098 7.289 -0.783) 7.959 3.626 ( -6.057 3.446 -1.069) 7.050 3.798 ( -0.161 -2.515 1.174) 2.780 3.964 ( -3.613 -4.800 -0.148) 6.010 4.526 ( 6.390 3.971 0.263) 7.528 5.326 ( 9.007 1.003 0.242) 9.066 5.465 ( 10.514 -1.940 -0.084) 10.692 6.356 ( 4.682 -0.425 0.014) 4.701 6.580 ( -0.364 -4.399 -0.009) 4.414 7.115 ( -1.078 5.147 -0.007) 5.259 7.242 ( -4.847 9.410 0.016) 10.584 7.676 ( 1.303 -1.865 -0.117) 2.278 7.851 ( -5.140 1.277 0.131) 5.298 7.899 ( -4.585 1.158 -0.133) 4.731 8.699 ( 8.470 -7.164 0.243) 11.096 9.864 ( 2.069 -8.174 -3.620) 9.176 10.025 ( 3.839 -9.687 3.530) 11.002 12.150 ( 0.744 -6.839 -0.627) 6.908 12.168 ( 1.878 -7.003 0.506) 7.268 ======================= Grid point 278 (59/73) ======================= q-point: ( 0.40 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.215 ( -3.039 10.682 -0.268) 11.110 2.256 ( -10.488 16.058 -0.214) 19.181 2.418 ( -0.376 5.057 0.861) 5.144 2.686 ( 3.824 6.833 -0.277) 7.836 2.827 ( -6.812 4.482 0.274) 8.159 2.961 ( 2.769 3.857 -0.720) 4.802 3.256 ( -7.317 14.055 0.539) 15.855 3.358 ( -3.858 12.765 -0.493) 13.344 3.471 ( -5.684 0.332 -0.383) 5.707 3.598 ( -5.070 6.258 0.219) 8.057 3.793 ( -2.988 6.904 0.389) 7.532 3.859 ( -5.458 1.275 -0.245) 5.611 4.629 ( 3.469 -0.529 0.045) 3.509 5.448 ( 5.830 -3.050 0.184) 6.582 5.583 ( 6.613 -5.305 0.069) 8.478 6.429 ( 4.065 -1.604 -0.008) 4.370 6.549 ( 1.346 -4.284 -0.005) 4.490 7.137 ( -2.029 4.628 -0.005) 5.053 7.250 ( -5.037 9.202 0.052) 10.490 7.682 ( 1.791 -2.604 -0.211) 3.168 7.799 ( -2.190 1.260 -0.136) 2.531 7.850 ( -2.625 1.700 0.138) 3.131 8.754 ( 7.379 -9.755 0.204) 12.233 9.831 ( 3.097 -7.172 -2.645) 8.248 9.995 ( 4.625 -10.087 2.653) 11.410 12.120 ( 2.753 -5.915 -0.047) 6.525 12.150 ( 3.403 -6.471 -0.080) 7.312 ======================= Grid point 279 (60/73) ======================= q-point: (-0.53 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.220 ( -7.934 9.499 -0.554) 12.389 2.258 ( -7.875 15.940 -0.298) 17.782 2.430 ( -3.030 2.883 0.135) 4.184 2.682 ( -8.501 0.188 -0.057) 8.504 2.831 ( -0.104 8.315 0.017) 8.315 2.946 ( -5.443 -1.198 0.258) 5.579 3.251 ( -9.406 12.872 0.754) 15.960 3.360 ( -7.975 9.715 -0.552) 12.581 3.477 ( 1.875 6.082 -0.770) 6.411 3.590 ( -3.277 7.152 0.679) 7.896 3.792 ( -4.406 5.695 0.437) 7.213 3.862 ( 1.717 5.679 -0.405) 5.947 4.630 ( -1.239 -3.285 0.012) 3.512 5.456 ( 0.446 -6.346 -0.295) 6.369 5.576 ( 0.672 -8.552 0.477) 8.592 6.429 ( -0.646 -4.322 -0.002) 4.370 6.550 ( 3.061 -3.304 -0.026) 4.504 7.137 ( -3.006 4.064 0.004) 5.055 7.250 ( -5.451 8.995 0.042) 10.517 7.681 ( 1.309 -2.917 -0.160) 3.201 7.803 ( 0.401 2.497 -0.394) 2.559 7.846 ( -0.541 3.163 0.384) 3.232 8.755 ( 4.890 -11.272 0.157) 12.288 9.822 ( 3.530 -6.494 -1.886) 7.628 10.004 ( 7.541 -8.793 1.873) 11.734 12.116 ( 3.212 -5.364 0.160) 6.254 12.153 ( 4.364 -6.172 -0.250) 7.563 ======================= Grid point 280 (61/73) ======================= q-point: (-0.47 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.165 ( -7.134 5.195 -0.641) 8.848 2.245 ( -11.788 14.361 -0.257) 18.581 2.374 ( -6.683 2.645 0.912) 7.244 2.490 ( -12.399 -1.522 -0.046) 12.493 2.818 ( -5.896 -1.426 0.271) 6.072 2.901 ( 0.002 4.571 -1.168) 4.718 3.168 ( -12.555 11.974 0.852) 17.370 3.335 ( -11.214 14.509 -0.345) 18.341 3.546 ( -4.889 5.791 -0.497) 7.595 3.607 ( -1.647 7.564 0.218) 7.745 3.809 ( 3.556 -0.712 0.437) 3.652 3.964 ( 5.422 0.694 -0.155) 5.468 4.524 ( -6.827 -3.674 0.343) 7.761 5.345 ( -4.934 -7.469 -0.929) 9.000 5.450 ( -3.940 -9.911 0.908) 10.704 6.356 ( -1.958 -4.255 0.009) 4.684 6.581 ( 3.969 -1.904 -0.024) 4.402 7.115 ( -3.911 3.504 0.009) 5.251 7.242 ( -5.766 8.946 -0.001) 10.643 7.674 ( 1.049 -2.067 -0.049) 2.319 7.858 ( 1.637 4.815 -0.305) 5.094 7.893 ( 1.106 4.837 0.283) 4.970 8.703 ( 2.270 -11.181 0.016) 11.409 9.825 ( 3.695 -5.180 -0.552) 6.386 10.063 ( 8.626 -8.709 0.520) 12.269 12.131 ( 4.328 -4.098 0.510) 5.982 12.185 ( 6.235 -5.052 -0.507) 8.040 ======================= Grid point 281 (62/73) ======================= q-point: (-0.40 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.074 ( -7.015 0.770 -0.073) 7.057 2.124 ( -13.651 6.224 -0.308) 15.006 2.190 ( -14.407 -5.523 0.766) 15.448 2.316 ( -11.891 12.487 -0.197) 17.244 2.688 ( -9.545 -0.517 0.733) 9.587 2.932 ( -3.496 6.131 -1.796) 7.283 3.077 ( -12.117 10.408 0.823) 15.995 3.231 ( -13.274 7.613 -0.793) 15.323 3.415 ( -10.484 -0.418 0.296) 10.497 3.664 ( -0.400 2.966 1.160) 3.210 3.821 ( -1.586 -4.310 -1.455) 4.817 4.058 ( 6.293 -2.290 0.216) 6.700 4.331 ( -10.704 -3.880 0.697) 11.406 5.168 ( -7.107 -4.986 -1.375) 8.789 5.275 ( -4.957 -8.331 1.300) 9.780 6.307 ( 0.679 -2.073 -0.020) 2.182 6.676 ( 8.454 0.452 0.060) 8.466 7.072 ( -4.572 2.593 0.015) 5.256 7.236 ( -4.494 8.522 0.032) 9.634 7.680 ( 0.936 -0.183 -0.078) 0.957 7.946 ( 2.028 5.911 -0.604) 6.279 7.970 ( 1.584 5.270 0.480) 5.524 8.622 ( 0.221 -8.674 -0.126) 8.678 9.842 ( 2.779 -3.318 1.595) 4.612 10.106 ( 6.350 -9.156 -1.585) 11.255 12.201 ( 7.649 -1.635 1.295) 7.929 12.282 ( 10.003 -2.950 -1.134) 10.490 ======================= Grid point 282 (63/73) ======================= q-point: (-0.33 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.843 ( -19.261 -3.545 -2.581) 19.754 1.929 ( -11.845 1.693 1.313) 12.037 1.945 ( -10.977 1.005 0.782) 11.050 2.258 ( -10.318 13.055 0.172) 16.641 2.465 ( -11.968 -5.930 0.735) 13.377 2.799 ( -21.929 7.338 -1.255) 23.158 2.873 ( -22.925 7.294 0.826) 24.072 3.089 ( -3.018 -6.618 -0.830) 7.321 3.251 ( -6.365 -2.411 0.515) 6.826 3.602 ( -9.523 -0.194 0.699) 9.551 3.655 ( -12.665 -3.333 -1.118) 13.143 4.051 ( -9.430 -0.069 0.408) 9.439 4.152 ( -0.743 -5.868 0.079) 5.915 5.056 ( -1.460 -0.121 -1.508) 2.102 5.140 ( -2.736 -4.526 1.497) 5.496 6.355 ( 6.524 -0.285 -0.051) 6.530 6.850 ( 8.562 0.307 0.066) 8.568 7.011 ( -4.701 1.389 0.017) 4.902 7.303 ( 5.635 7.023 0.388) 9.013 7.698 ( 0.400 1.081 -0.391) 1.216 8.022 ( 1.456 1.993 -0.765) 2.584 8.045 ( 1.935 3.952 0.227) 4.406 8.557 ( -0.458 -3.000 0.401) 3.061 9.838 ( -1.254 -1.019 3.517) 3.870 10.099 ( 2.576 -7.410 -3.842) 8.735 12.371 ( 12.516 0.913 2.359) 12.769 12.480 ( 15.254 -0.945 -1.964) 15.409 ======================= Grid point 290 (64/73) ======================= q-point: ( 0.20 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.076 ( 2.795 8.828 -0.016) 9.260 2.116 ( 6.455 9.264 2.040) 11.474 2.173 ( 8.647 12.750 -2.054) 15.542 2.577 ( 3.568 10.711 -0.149) 11.290 2.626 ( 9.462 11.025 0.721) 14.546 2.955 ( -3.829 -5.124 -1.217) 6.512 3.275 ( 13.925 17.851 -0.585) 22.648 3.307 ( 2.680 11.866 0.954) 12.202 3.351 ( 8.165 13.976 0.591) 16.197 3.727 ( 4.093 2.001 0.549) 4.589 3.747 ( -1.227 2.210 -1.167) 2.784 4.032 ( -3.149 -5.769 0.028) 6.573 4.287 ( 6.622 11.153 0.342) 12.975 5.118 ( 2.114 5.730 -0.786) 6.157 5.180 ( 2.749 2.918 0.795) 4.087 6.285 ( -0.827 -1.403 -0.012) 1.629 6.681 ( -4.769 -8.301 0.054) 9.573 7.102 ( 2.433 4.207 0.014) 4.860 7.329 ( 1.592 2.650 0.081) 3.093 7.693 ( -0.548 -0.707 -0.182) 0.913 8.002 ( -1.467 -1.723 -0.488) 2.315 8.035 ( -0.842 -2.051 0.218) 2.228 8.526 ( -0.885 -1.872 0.295) 2.092 9.817 ( -0.534 -3.826 1.297) 4.074 9.984 ( -3.902 -3.662 -1.461) 5.547 12.193 ( -3.628 -8.267 1.605) 9.169 12.240 ( -6.265 -9.057 -1.472) 11.111 ======================= Grid point 291 (65/73) ======================= q-point: ( 0.27 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.198 ( 2.293 8.285 0.428) 8.607 2.325 ( 6.227 9.319 1.503) 11.308 2.427 ( 5.862 9.888 -1.518) 11.595 2.667 ( -3.201 11.193 -0.122) 11.642 2.803 ( -0.776 5.243 -0.786) 5.358 2.930 ( 3.779 3.460 0.286) 5.132 3.376 ( -1.125 8.231 0.433) 8.319 3.590 ( -2.352 16.666 0.397) 16.835 3.642 ( 3.202 11.336 -1.145) 11.835 3.724 ( 2.598 2.975 -0.579) 3.992 3.759 ( -1.007 -1.449 0.900) 1.981 3.958 ( 0.332 -0.886 0.135) 0.955 4.480 ( 4.809 6.242 0.154) 7.881 5.222 ( 4.300 3.346 0.139) 5.450 5.286 ( 7.553 -0.321 -0.108) 7.560 6.291 ( 2.072 -0.141 0.022) 2.077 6.552 ( -2.163 -5.397 -0.024) 5.814 7.172 ( 0.562 4.545 0.006) 4.580 7.382 ( -0.799 6.126 -0.042) 6.178 7.667 ( -1.596 -0.353 -0.014) 1.634 7.933 ( -3.921 -2.211 -0.527) 4.532 7.975 ( -3.622 -1.418 0.476) 3.919 8.507 ( 3.774 -6.530 0.054) 7.542 9.759 ( 0.339 -6.554 -0.334) 6.571 9.896 ( -0.027 -7.495 0.312) 7.502 12.084 ( -0.952 -6.478 1.056) 6.632 12.093 ( -1.746 -6.368 -1.048) 6.686 ======================= Grid point 292 (66/73) ======================= q-point: ( 0.33 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.315 ( 0.549 9.294 -0.017) 9.310 2.466 ( -0.554 7.725 0.489) 7.761 2.531 ( -1.494 5.510 -0.188) 5.712 2.750 ( 2.610 8.597 -0.231) 8.988 2.902 ( 3.900 6.439 -0.389) 7.538 2.979 ( -4.208 7.585 0.186) 8.676 3.435 ( -0.813 6.608 0.324) 6.665 3.520 ( -8.194 0.134 -0.560) 8.214 3.626 ( -6.073 -0.557 0.628) 6.131 3.766 ( -4.909 13.190 -0.606) 14.087 3.922 ( 2.148 13.314 -0.450) 13.493 3.996 ( -3.034 9.654 0.767) 10.148 4.570 ( 2.402 -0.056 -0.020) 2.403 5.319 ( 5.436 -1.578 0.210) 5.664 5.403 ( 7.752 -4.117 -0.117) 8.778 6.341 ( 3.963 -0.913 0.015) 4.067 6.488 ( 0.177 -4.760 -0.023) 4.764 7.211 ( -1.069 4.250 -0.005) 4.382 7.418 ( -3.165 7.916 0.019) 8.526 7.640 ( -0.301 -1.554 -0.138) 1.589 7.845 ( -3.134 -1.858 -0.468) 3.673 7.903 ( -2.829 -0.799 0.515) 2.985 8.519 ( 6.849 -10.187 0.038) 12.275 9.711 ( 1.960 -6.651 -0.879) 6.989 9.832 ( 2.286 -9.148 0.909) 9.473 12.025 ( 1.211 -5.358 0.357) 5.505 12.042 ( 1.674 -5.519 -0.393) 5.780 ======================= Grid point 293 (67/73) ======================= q-point: ( 0.40 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.372 ( -4.348 8.173 -0.616) 9.278 2.514 ( -1.375 6.657 0.207) 6.800 2.521 ( -5.094 5.032 -0.309) 7.167 2.879 ( -2.848 -0.044 0.050) 2.848 2.890 ( 0.144 4.696 0.357) 4.712 3.027 ( -7.152 11.032 0.263) 13.150 3.422 ( -1.518 3.384 -0.578) 3.754 3.480 ( -2.006 2.833 0.515) 3.509 3.544 ( -0.842 1.512 0.479) 1.796 3.790 ( -8.858 15.590 -0.894) 17.953 4.007 ( -7.607 10.780 0.250) 13.196 4.022 ( -6.531 13.557 0.186) 15.049 4.585 ( 2.194 -3.776 -0.083) 4.368 5.363 ( 3.207 -5.097 -0.004) 6.022 5.450 ( 4.151 -7.554 0.096) 8.620 6.382 ( 1.733 -3.014 -0.000) 3.476 6.460 ( 2.439 -4.201 -0.013) 4.858 7.223 ( -2.157 3.722 -0.000) 4.302 7.431 ( -4.843 8.495 0.044) 9.778 7.629 ( 1.268 -2.366 -0.165) 2.690 7.807 ( 0.246 -0.435 -0.489) 0.699 7.874 ( -0.339 0.630 0.540) 0.896 8.529 ( 6.866 -11.911 0.044) 13.748 9.696 ( 3.276 -5.838 -0.806) 6.743 9.804 ( 5.050 -8.516 0.838) 9.936 12.011 ( 2.765 -4.623 0.529) 5.413 12.028 ( 2.950 -5.300 -0.565) 6.092 ======================= Grid point 295 (68/73) ======================= q-point: (-0.53 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.326 ( -7.835 4.005 -0.671) 8.825 2.480 ( -5.928 4.467 -0.451) 7.436 2.518 ( -4.898 3.716 0.612) 6.179 2.750 ( -7.932 1.160 -0.251) 8.020 2.886 ( -8.480 1.070 0.584) 8.567 2.986 ( -3.861 7.490 -0.185) 8.429 3.426 ( -5.993 4.109 0.831) 7.314 3.521 ( 3.407 7.363 -0.579) 8.134 3.631 ( 4.875 4.033 0.371) 6.337 3.769 ( -9.666 11.853 -0.849) 15.318 3.911 ( -12.228 4.511 0.302) 13.037 4.004 ( -7.111 7.356 0.199) 10.233 4.569 ( -1.213 -2.151 0.000) 2.469 5.327 ( -0.859 -5.668 -0.284) 5.739 5.396 ( -0.728 -8.637 0.287) 8.672 6.341 ( -1.184 -3.884 0.016) 4.060 6.488 ( 4.037 -2.547 -0.035) 4.773 7.211 ( -3.154 3.053 0.009) 4.390 7.420 ( -5.294 6.766 -0.068) 8.591 7.638 ( 1.499 -0.559 -0.000) 1.600 7.845 ( 3.342 1.541 -0.465) 3.709 7.904 ( 1.953 2.312 0.477) 3.064 8.519 ( 5.431 -11.093 0.052) 12.351 9.705 ( 3.990 -4.610 -0.452) 6.114 9.839 ( 7.587 -6.961 0.427) 10.305 12.024 ( 3.437 -3.464 0.379) 4.895 12.042 ( 4.499 -4.452 -0.386) 6.341 ======================= Grid point 296 (69/73) ======================= q-point: (-0.47 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.210 ( -7.930 1.167 -0.342) 8.023 2.333 ( -11.276 -2.075 1.524) 11.566 2.409 ( -12.162 1.964 -0.892) 12.352 2.672 ( -7.539 7.984 -0.497) 10.992 2.791 ( -4.905 3.700 0.068) 6.144 2.942 ( -4.195 -2.008 -0.543) 4.683 3.365 ( -6.656 5.914 1.069) 8.968 3.601 ( -9.707 10.429 0.109) 14.248 3.631 ( -14.894 2.437 -0.864) 15.117 3.719 ( -4.743 0.937 -0.122) 4.836 3.769 ( 3.030 -1.466 0.160) 3.370 3.958 ( 0.159 -0.711 0.164) 0.747 4.479 ( -7.892 -0.920 0.247) 7.949 5.235 ( -4.831 -2.413 -0.740) 5.450 5.275 ( -3.688 -6.478 0.684) 7.486 6.291 ( -0.905 -1.875 0.014) 2.082 6.552 ( 5.744 -0.829 -0.024) 5.803 7.172 ( -4.213 1.801 0.016) 4.582 7.383 ( -4.971 3.714 -0.075) 6.206 7.667 ( 1.226 1.266 0.022) 1.762 7.934 ( 4.119 2.107 -0.637) 4.670 7.974 ( 2.803 2.588 0.549) 3.854 8.506 ( 3.753 -6.581 0.093) 7.576 9.746 ( 4.653 -2.125 0.423) 5.133 9.910 ( 7.363 -5.205 -0.499) 9.031 12.072 ( 5.383 -1.376 0.763) 5.609 12.104 ( 7.045 -2.629 -0.707) 7.553 ======================= Grid point 307 (70/73) ======================= q-point: ( 0.27 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.364 ( 4.319 8.254 -0.173) 9.317 2.505 ( 4.441 8.246 0.927) 9.411 2.614 ( 4.068 7.751 -0.406) 8.763 2.666 ( -0.044 -1.886 -0.878) 2.081 2.940 ( 0.703 0.774 0.773) 1.300 3.086 ( 6.132 11.194 -0.533) 12.774 3.468 ( 2.113 2.932 0.738) 3.688 3.601 ( -2.708 -4.458 -0.715) 5.265 3.666 ( -3.347 -5.876 0.576) 6.787 3.910 ( 4.761 10.370 0.268) 11.414 3.946 ( 8.626 12.590 -1.479) 15.333 4.092 ( 5.918 10.769 0.990) 12.328 4.552 ( 0.292 0.442 -0.050) 0.532 5.271 ( -0.205 0.938 0.254) 0.993 5.288 ( 1.367 1.173 -0.265) 1.821 6.299 ( 0.564 0.979 0.024) 1.130 6.461 ( -2.190 -3.786 -0.038) 4.374 7.244 ( 1.502 2.584 0.006) 2.989 7.479 ( 2.146 3.884 -0.157) 4.440 7.655 ( -0.555 -1.176 0.105) 1.304 7.860 ( -2.405 -4.391 -0.755) 5.063 7.927 ( -1.890 -2.974 0.831) 3.620 8.390 ( -2.593 -4.578 -0.068) 5.262 9.657 ( -1.958 -3.564 0.162) 4.070 9.763 ( -3.429 -5.683 -0.152) 6.639 11.989 ( -2.028 -3.159 0.031) 3.754 11.995 ( -1.823 -3.538 -0.026) 3.980 ======================= Grid point 308 (71/73) ======================= q-point: ( 0.33 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.507 ( 1.027 9.528 -0.468) 9.595 2.606 ( -4.158 5.030 0.112) 6.527 2.616 ( -0.571 4.744 -0.003) 4.779 2.788 ( 8.264 -0.947 -0.214) 8.321 2.892 ( -2.390 -5.385 0.760) 5.940 3.281 ( 0.944 15.033 -0.341) 15.066 3.491 ( -6.003 -1.412 -0.103) 6.167 3.514 ( 0.328 2.175 0.390) 2.234 3.566 ( -2.430 -3.228 0.282) 4.050 4.059 ( -2.727 13.924 -0.527) 14.198 4.163 ( 2.768 10.037 -1.151) 10.475 4.252 ( -1.808 12.061 1.403) 12.276 4.506 ( -0.497 -6.569 -0.124) 6.589 5.275 ( 1.510 -2.303 0.010) 2.754 5.314 ( 3.317 -4.341 -0.017) 5.463 6.327 ( 1.959 -0.267 0.011) 1.977 6.399 ( -0.360 -3.923 -0.012) 3.940 7.276 ( -0.429 2.173 -0.006) 2.215 7.546 ( 0.701 4.783 -0.049) 4.834 7.618 ( -2.492 -0.681 -0.044) 2.584 7.792 ( -0.706 -2.948 -0.858) 3.151 7.876 ( -0.788 -1.560 1.067) 2.047 8.334 ( 2.802 -7.399 -0.154) 7.913 9.615 ( 0.701 -2.782 0.050) 2.869 9.679 ( 0.622 -5.910 -0.014) 5.943 11.945 ( 0.417 -2.569 0.231) 2.613 11.955 ( 0.605 -3.070 -0.238) 3.138 ======================= Grid point 309 (72/73) ======================= q-point: (-0.60 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.509 ( -8.432 3.694 -0.563) 9.223 2.610 ( -1.557 3.958 -0.410) 4.273 2.618 ( -4.354 5.459 0.085) 6.983 2.775 ( -3.379 -7.575 0.602) 8.316 2.896 ( 5.405 -0.963 0.469) 5.510 3.279 ( -13.360 6.801 -0.264) 14.993 3.497 ( 4.503 4.713 -0.551) 6.542 3.513 ( -2.495 0.534 0.481) 2.596 3.561 ( 3.976 0.504 0.596) 4.052 4.061 ( -10.416 9.314 -0.679) 13.989 4.158 ( -10.493 2.812 -0.824) 10.895 4.254 ( -9.322 7.383 1.200) 11.953 4.506 ( 5.928 -2.859 -0.126) 6.583 5.276 ( 1.393 -2.552 -0.031) 2.908 5.314 ( 1.979 -4.969 0.011) 5.349 6.327 ( -0.749 -1.827 0.013) 1.975 6.399 ( 3.581 -1.654 -0.015) 3.944 7.276 ( -1.674 1.462 -0.003) 2.222 7.548 ( -4.466 1.695 -0.139) 4.778 7.616 ( 1.787 1.925 0.051) 2.628 7.790 ( 2.812 -0.840 -0.756) 3.031 7.877 ( 1.863 -0.133 0.957) 2.099 8.333 ( 4.980 -6.117 -0.138) 7.889 9.616 ( 2.062 -1.998 0.003) 2.871 9.679 ( 4.832 -3.509 0.020) 5.972 11.949 ( 2.321 -1.595 0.312) 2.834 11.952 ( 2.045 -2.106 -0.317) 2.953 ======================= Grid point 323 (73/73) ======================= q-point: ( 0.33 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 116 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.639 ( -0.646 -1.293 -0.799) 1.652 2.639 ( 0.646 1.293 -0.799) 1.652 2.744 ( -2.461 -4.922 0.848) 5.568 2.744 ( 2.461 4.922 0.848) 5.568 2.773 ( -0.000 0.000 0.077) 0.077 3.478 ( -0.000 0.000 -0.244) 0.244 3.480 ( -1.122 -2.246 0.045) 2.511 3.480 ( 1.123 2.246 0.045) 2.511 3.548 ( -0.000 -0.000 0.459) 0.459 4.267 ( 0.000 -0.000 -1.779) 1.779 4.319 ( -4.981 -9.963 -0.061) 11.139 4.319 ( 4.981 9.963 -0.061) 11.139 4.383 ( -0.000 -0.000 1.503) 1.503 5.256 ( -0.439 -0.878 -0.026) 0.982 5.256 ( 0.439 0.878 -0.026) 0.982 6.341 ( -0.812 -1.624 0.004) 1.816 6.341 ( 0.812 1.624 0.004) 1.816 7.297 ( -0.000 -0.000 -0.011) 0.011 7.605 ( -0.523 -1.045 -0.014) 1.169 7.605 ( 0.523 1.045 -0.014) 1.169 7.756 ( -0.000 -0.000 -1.092) 1.092 7.858 ( -0.000 -0.000 1.374) 1.374 8.256 ( -0.000 -0.000 -0.271) 0.271 9.604 ( -0.719 -1.437 0.012) 1.607 9.604 ( 0.719 1.437 0.012) 1.607 11.922 ( -0.092 -0.184 -0.001) 0.206 11.922 ( 0.092 0.184 -0.001) 0.206 =================== 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/18225 10.0 441.370 441.370 74.753 0.000 -0.000 0.000 3/18225 20.0 198.733 198.733 16.967 0.000 -0.000 0.000 3/18225 30.0 115.350 115.350 7.317 0.000 -0.000 0.000 3/18225 40.0 80.325 80.325 4.370 0.000 -0.000 0.000 3/18225 50.0 60.742 60.742 3.051 0.000 -0.000 0.000 3/18225 60.0 48.006 48.006 2.321 0.000 -0.000 0.000 3/18225 70.0 39.176 39.176 1.864 0.000 -0.000 0.000 3/18225 80.0 32.825 32.825 1.554 0.000 -0.000 0.000 3/18225 90.0 28.122 28.122 1.333 0.000 -0.000 0.000 3/18225 100.0 24.543 24.543 1.168 0.000 -0.000 0.000 3/18225 110.0 21.753 21.753 1.041 0.000 -0.000 0.000 3/18225 120.0 19.528 19.528 0.940 0.000 -0.000 0.000 3/18225 130.0 17.720 17.720 0.858 0.000 -0.000 0.000 3/18225 140.0 16.223 16.223 0.790 0.000 -0.000 0.000 3/18225 150.0 14.967 14.967 0.733 0.000 -0.000 0.000 3/18225 160.0 13.896 13.896 0.684 0.000 -0.000 0.000 3/18225 170.0 12.974 12.974 0.641 0.000 -0.000 0.000 3/18225 180.0 12.171 12.171 0.604 0.000 -0.000 0.000 3/18225 190.0 11.466 11.466 0.571 0.000 -0.000 0.000 3/18225 200.0 10.841 10.841 0.541 0.000 -0.000 0.000 3/18225 210.0 10.283 10.283 0.515 0.000 -0.000 0.000 3/18225 220.0 9.782 9.782 0.491 0.000 -0.000 0.000 3/18225 230.0 9.329 9.329 0.470 0.000 -0.000 0.000 3/18225 240.0 8.918 8.918 0.450 0.000 -0.000 0.000 3/18225 250.0 8.542 8.542 0.432 0.000 -0.000 0.000 3/18225 260.0 8.199 8.199 0.415 0.000 -0.000 0.000 3/18225 270.0 7.882 7.882 0.400 0.000 -0.000 0.000 3/18225 280.0 7.590 7.590 0.385 0.000 -0.000 0.000 3/18225 290.0 7.319 7.319 0.372 0.000 -0.000 0.000 3/18225 300.0 7.068 7.068 0.360 0.000 -0.000 0.000 3/18225 310.0 6.834 6.834 0.348 0.000 -0.000 0.000 3/18225 320.0 6.615 6.615 0.337 0.000 -0.000 0.000 3/18225 330.0 6.410 6.410 0.327 0.000 -0.000 0.000 3/18225 340.0 6.218 6.218 0.318 0.000 -0.000 0.000 3/18225 350.0 6.037 6.037 0.309 0.000 -0.000 0.000 3/18225 360.0 5.866 5.866 0.300 0.000 -0.000 0.000 3/18225 370.0 5.705 5.705 0.292 0.000 -0.000 0.000 3/18225 380.0 5.553 5.553 0.284 0.000 -0.000 0.000 3/18225 390.0 5.409 5.409 0.277 0.000 -0.000 0.000 3/18225 400.0 5.273 5.273 0.270 0.000 -0.000 0.000 3/18225 410.0 5.143 5.143 0.264 0.000 -0.000 0.000 3/18225 420.0 5.019 5.019 0.257 0.000 -0.000 0.000 3/18225 430.0 4.902 4.902 0.252 0.000 -0.000 0.000 3/18225 440.0 4.790 4.790 0.246 0.000 -0.000 0.000 3/18225 450.0 4.682 4.682 0.240 0.000 -0.000 0.000 3/18225 460.0 4.580 4.580 0.235 0.000 -0.000 0.000 3/18225 470.0 4.482 4.482 0.230 0.000 -0.000 0.000 3/18225 480.0 4.389 4.389 0.226 0.000 -0.000 0.000 3/18225 490.0 4.299 4.299 0.221 0.000 -0.000 0.000 3/18225 500.0 4.213 4.213 0.217 0.000 -0.000 0.000 3/18225 510.0 4.130 4.130 0.212 0.000 -0.000 0.000 3/18225 520.0 4.050 4.050 0.208 0.000 -0.000 0.000 3/18225 530.0 3.974 3.974 0.204 0.000 -0.000 0.000 3/18225 540.0 3.900 3.900 0.201 0.000 -0.000 0.000 3/18225 550.0 3.829 3.829 0.197 0.000 -0.000 0.000 3/18225 560.0 3.761 3.761 0.193 0.000 -0.000 0.000 3/18225 570.0 3.695 3.695 0.190 0.000 -0.000 0.000 3/18225 580.0 3.631 3.631 0.187 0.000 -0.000 0.000 3/18225 590.0 3.570 3.570 0.184 0.000 -0.000 0.000 3/18225 600.0 3.510 3.510 0.181 0.000 -0.000 0.000 3/18225 610.0 3.453 3.453 0.178 0.000 -0.000 0.000 3/18225 620.0 3.397 3.397 0.175 0.000 -0.000 0.000 3/18225 630.0 3.344 3.344 0.172 0.000 -0.000 0.000 3/18225 640.0 3.291 3.291 0.169 0.000 -0.000 0.000 3/18225 650.0 3.241 3.241 0.167 0.000 -0.000 0.000 3/18225 660.0 3.192 3.192 0.164 0.000 -0.000 0.000 3/18225 670.0 3.144 3.144 0.162 0.000 -0.000 0.000 3/18225 680.0 3.098 3.098 0.160 0.000 -0.000 0.000 3/18225 690.0 3.054 3.054 0.157 0.000 -0.000 0.000 3/18225 700.0 3.010 3.010 0.155 0.000 -0.000 0.000 3/18225 710.0 2.968 2.968 0.153 0.000 -0.000 0.000 3/18225 720.0 2.927 2.927 0.151 0.000 -0.000 0.000 3/18225 730.0 2.887 2.887 0.149 0.000 -0.000 0.000 3/18225 740.0 2.848 2.848 0.147 0.000 -0.000 0.000 3/18225 750.0 2.810 2.810 0.145 0.000 -0.000 0.000 3/18225 760.0 2.773 2.773 0.143 0.000 -0.000 0.000 3/18225 770.0 2.737 2.737 0.141 0.000 -0.000 0.000 3/18225 780.0 2.702 2.702 0.139 0.000 -0.000 0.000 3/18225 790.0 2.668 2.668 0.137 0.000 -0.000 0.000 3/18225 800.0 2.635 2.635 0.136 0.000 -0.000 0.000 3/18225 810.0 2.603 2.603 0.134 0.000 -0.000 0.000 3/18225 820.0 2.571 2.571 0.132 0.000 -0.000 0.000 3/18225 830.0 2.540 2.540 0.131 0.000 -0.000 0.000 3/18225 840.0 2.510 2.510 0.129 0.000 -0.000 0.000 3/18225 850.0 2.481 2.481 0.128 0.000 -0.000 0.000 3/18225 860.0 2.452 2.452 0.126 0.000 -0.000 0.000 3/18225 870.0 2.424 2.424 0.125 0.000 -0.000 0.000 3/18225 880.0 2.397 2.397 0.123 0.000 -0.000 0.000 3/18225 890.0 2.370 2.370 0.122 0.000 -0.000 0.000 3/18225 900.0 2.344 2.344 0.121 0.000 -0.000 0.000 3/18225 910.0 2.318 2.318 0.119 0.000 -0.000 0.000 3/18225 920.0 2.293 2.293 0.118 0.000 -0.000 0.000 3/18225 930.0 2.268 2.268 0.117 0.000 -0.000 0.000 3/18225 940.0 2.244 2.244 0.116 0.000 -0.000 0.000 3/18225 950.0 2.221 2.221 0.114 0.000 -0.000 0.000 3/18225 960.0 2.198 2.198 0.113 0.000 -0.000 0.000 3/18225 970.0 2.175 2.175 0.112 0.000 -0.000 0.000 3/18225 980.0 2.153 2.153 0.111 0.000 -0.000 0.000 3/18225 990.0 2.132 2.132 0.110 0.000 -0.000 0.000 3/18225 1000.0 2.110 2.110 0.109 0.000 -0.000 0.000 3/18225 Thermal conductivity related properties were written into "kappa-m15153.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 08:17:46]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|