----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 10:50:34]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [3 3 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P4/mmm (123) Number of symmetry operations in supercell: 288 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.895413370000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.895413370000000 0.000000000000000 c 0.000000000000000 0.000000000000000 9.483492630000001 Atomic positions (fractional): *1 Bi 0.50000000000000 0.50000000000000 0.73986464999698 208.980 2 Bi 0.50000000000000 0.50000000000000 0.26013535000302 208.980 *3 O 0.00000000000000 0.50000000000000 0.85671206518205 15.999 4 O 0.00000000000000 0.50000000000000 0.14328793481795 15.999 5 O 0.50000000000000 0.00000000000000 0.85671206518205 15.999 6 O 0.50000000000000 0.00000000000000 0.14328793481795 15.999 *7 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 *8 I 0.00000000000000 0.00000000000000 0.50000000000000 126.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.895413370000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.895413370000000 0.000000000000000 c 0.000000000000000 0.000000000000000 9.483492630000001 Atomic positions (fractional): *1 Bi 0.50000000000000 0.50000000000000 0.73986464999698 208.980 > 1 2 Bi 0.50000000000000 0.50000000000000 0.26013535000302 208.980 > 2 *3 O 0.00000000000000 0.50000000000000 0.85671206518205 15.999 > 3 4 O 0.00000000000000 0.50000000000000 0.14328793481795 15.999 > 4 5 O 0.50000000000000 0.00000000000000 0.85671206518205 15.999 > 5 6 O 0.50000000000000 0.00000000000000 0.14328793481795 15.999 > 6 *7 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 > 7 *8 I 0.00000000000000 0.00000000000000 0.50000000000000 126.904 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 11.686240110000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.686240110000000 0.000000000000000 c 0.000000000000000 0.000000000000000 18.966985260000001 Atomic positions (fractional): *1 Bi 0.16666666666667 0.16666666666667 0.36993232499849 208.980 > 1 2 Bi 0.50000000000000 0.16666666666667 0.36993232499849 208.980 > 1 3 Bi 0.83333333333333 0.16666666666667 0.36993232499849 208.980 > 1 4 Bi 0.16666666666667 0.50000000000000 0.36993232499849 208.980 > 1 5 Bi 0.50000000000000 0.50000000000000 0.36993232499849 208.980 > 1 6 Bi 0.83333333333333 0.50000000000000 0.36993232499849 208.980 > 1 7 Bi 0.16666666666667 0.83333333333333 0.36993232499849 208.980 > 1 8 Bi 0.50000000000000 0.83333333333333 0.36993232499849 208.980 > 1 9 Bi 0.83333333333333 0.83333333333333 0.36993232499849 208.980 > 1 10 Bi 0.16666666666667 0.16666666666667 0.86993232499849 208.980 > 1 11 Bi 0.50000000000000 0.16666666666667 0.86993232499849 208.980 > 1 12 Bi 0.83333333333333 0.16666666666667 0.86993232499849 208.980 > 1 13 Bi 0.16666666666667 0.50000000000000 0.86993232499849 208.980 > 1 14 Bi 0.50000000000000 0.50000000000000 0.86993232499849 208.980 > 1 15 Bi 0.83333333333333 0.50000000000000 0.86993232499849 208.980 > 1 16 Bi 0.16666666666667 0.83333333333333 0.86993232499849 208.980 > 1 17 Bi 0.50000000000000 0.83333333333333 0.86993232499849 208.980 > 1 18 Bi 0.83333333333333 0.83333333333333 0.86993232499849 208.980 > 1 19 Bi 0.16666666666667 0.16666666666667 0.13006767500151 208.980 > 2 20 Bi 0.50000000000000 0.16666666666667 0.13006767500151 208.980 > 2 21 Bi 0.83333333333333 0.16666666666667 0.13006767500151 208.980 > 2 22 Bi 0.16666666666667 0.50000000000000 0.13006767500151 208.980 > 2 23 Bi 0.50000000000000 0.50000000000000 0.13006767500151 208.980 > 2 24 Bi 0.83333333333333 0.50000000000000 0.13006767500151 208.980 > 2 25 Bi 0.16666666666667 0.83333333333333 0.13006767500151 208.980 > 2 26 Bi 0.50000000000000 0.83333333333333 0.13006767500151 208.980 > 2 27 Bi 0.83333333333333 0.83333333333333 0.13006767500151 208.980 > 2 28 Bi 0.16666666666667 0.16666666666667 0.63006767500151 208.980 > 2 29 Bi 0.50000000000000 0.16666666666667 0.63006767500151 208.980 > 2 30 Bi 0.83333333333333 0.16666666666667 0.63006767500151 208.980 > 2 31 Bi 0.16666666666667 0.50000000000000 0.63006767500151 208.980 > 2 32 Bi 0.50000000000000 0.50000000000000 0.63006767500151 208.980 > 2 33 Bi 0.83333333333333 0.50000000000000 0.63006767500151 208.980 > 2 34 Bi 0.16666666666667 0.83333333333333 0.63006767500151 208.980 > 2 35 Bi 0.50000000000000 0.83333333333333 0.63006767500151 208.980 > 2 36 Bi 0.83333333333333 0.83333333333333 0.63006767500151 208.980 > 2 *37 O 0.00000000000000 0.16666666666667 0.42835603259102 15.999 > 3 38 O 0.33333333333333 0.16666666666667 0.42835603259102 15.999 > 3 39 O 0.66666666666667 0.16666666666667 0.42835603259102 15.999 > 3 40 O 0.00000000000000 0.50000000000000 0.42835603259102 15.999 > 3 41 O 0.33333333333333 0.50000000000000 0.42835603259102 15.999 > 3 42 O 0.66666666666667 0.50000000000000 0.42835603259102 15.999 > 3 43 O 0.00000000000000 0.83333333333333 0.42835603259102 15.999 > 3 44 O 0.33333333333333 0.83333333333333 0.42835603259102 15.999 > 3 45 O 0.66666666666667 0.83333333333333 0.42835603259102 15.999 > 3 46 O 0.00000000000000 0.16666666666667 0.92835603259102 15.999 > 3 47 O 0.33333333333333 0.16666666666667 0.92835603259102 15.999 > 3 48 O 0.66666666666667 0.16666666666667 0.92835603259102 15.999 > 3 49 O 0.00000000000000 0.50000000000000 0.92835603259102 15.999 > 3 50 O 0.33333333333333 0.50000000000000 0.92835603259102 15.999 > 3 51 O 0.66666666666667 0.50000000000000 0.92835603259102 15.999 > 3 52 O 0.00000000000000 0.83333333333333 0.92835603259102 15.999 > 3 53 O 0.33333333333333 0.83333333333333 0.92835603259102 15.999 > 3 54 O 0.66666666666667 0.83333333333333 0.92835603259102 15.999 > 3 55 O 0.00000000000000 0.16666666666667 0.07164396740898 15.999 > 4 56 O 0.33333333333333 0.16666666666667 0.07164396740898 15.999 > 4 57 O 0.66666666666667 0.16666666666667 0.07164396740898 15.999 > 4 58 O 0.00000000000000 0.50000000000000 0.07164396740898 15.999 > 4 59 O 0.33333333333333 0.50000000000000 0.07164396740898 15.999 > 4 60 O 0.66666666666667 0.50000000000000 0.07164396740898 15.999 > 4 61 O 0.00000000000000 0.83333333333333 0.07164396740898 15.999 > 4 62 O 0.33333333333333 0.83333333333333 0.07164396740898 15.999 > 4 63 O 0.66666666666667 0.83333333333333 0.07164396740898 15.999 > 4 64 O 0.00000000000000 0.16666666666667 0.57164396740898 15.999 > 4 65 O 0.33333333333333 0.16666666666667 0.57164396740898 15.999 > 4 66 O 0.66666666666667 0.16666666666667 0.57164396740898 15.999 > 4 67 O 0.00000000000000 0.50000000000000 0.57164396740898 15.999 > 4 68 O 0.33333333333333 0.50000000000000 0.57164396740898 15.999 > 4 69 O 0.66666666666667 0.50000000000000 0.57164396740898 15.999 > 4 70 O 0.00000000000000 0.83333333333333 0.57164396740898 15.999 > 4 71 O 0.33333333333333 0.83333333333333 0.57164396740898 15.999 > 4 72 O 0.66666666666667 0.83333333333333 0.57164396740898 15.999 > 4 73 O 0.16666666666667 0.00000000000000 0.42835603259102 15.999 > 5 74 O 0.50000000000000 0.00000000000000 0.42835603259102 15.999 > 5 75 O 0.83333333333333 0.00000000000000 0.42835603259102 15.999 > 5 76 O 0.16666666666667 0.33333333333333 0.42835603259102 15.999 > 5 77 O 0.50000000000000 0.33333333333333 0.42835603259102 15.999 > 5 78 O 0.83333333333333 0.33333333333333 0.42835603259102 15.999 > 5 79 O 0.16666666666667 0.66666666666667 0.42835603259102 15.999 > 5 80 O 0.50000000000000 0.66666666666667 0.42835603259102 15.999 > 5 81 O 0.83333333333333 0.66666666666667 0.42835603259102 15.999 > 5 82 O 0.16666666666667 0.00000000000000 0.92835603259102 15.999 > 5 83 O 0.50000000000000 0.00000000000000 0.92835603259102 15.999 > 5 84 O 0.83333333333333 0.00000000000000 0.92835603259102 15.999 > 5 85 O 0.16666666666667 0.33333333333333 0.92835603259102 15.999 > 5 86 O 0.50000000000000 0.33333333333333 0.92835603259102 15.999 > 5 87 O 0.83333333333333 0.33333333333333 0.92835603259102 15.999 > 5 88 O 0.16666666666667 0.66666666666667 0.92835603259102 15.999 > 5 89 O 0.50000000000000 0.66666666666667 0.92835603259102 15.999 > 5 90 O 0.83333333333333 0.66666666666667 0.92835603259102 15.999 > 5 91 O 0.16666666666667 0.00000000000000 0.07164396740898 15.999 > 6 92 O 0.50000000000000 0.00000000000000 0.07164396740898 15.999 > 6 93 O 0.83333333333333 0.00000000000000 0.07164396740898 15.999 > 6 94 O 0.16666666666667 0.33333333333333 0.07164396740898 15.999 > 6 95 O 0.50000000000000 0.33333333333333 0.07164396740898 15.999 > 6 96 O 0.83333333333333 0.33333333333333 0.07164396740898 15.999 > 6 97 O 0.16666666666667 0.66666666666667 0.07164396740898 15.999 > 6 98 O 0.50000000000000 0.66666666666667 0.07164396740898 15.999 > 6 99 O 0.83333333333333 0.66666666666667 0.07164396740898 15.999 > 6 100 O 0.16666666666667 0.00000000000000 0.57164396740898 15.999 > 6 101 O 0.50000000000000 0.00000000000000 0.57164396740898 15.999 > 6 102 O 0.83333333333333 0.00000000000000 0.57164396740898 15.999 > 6 103 O 0.16666666666667 0.33333333333333 0.57164396740898 15.999 > 6 104 O 0.50000000000000 0.33333333333333 0.57164396740898 15.999 > 6 105 O 0.83333333333333 0.33333333333333 0.57164396740898 15.999 > 6 106 O 0.16666666666667 0.66666666666667 0.57164396740898 15.999 > 6 107 O 0.50000000000000 0.66666666666667 0.57164396740898 15.999 > 6 108 O 0.83333333333333 0.66666666666667 0.57164396740898 15.999 > 6 *109 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 > 7 110 Y 0.33333333333333 0.00000000000000 0.00000000000000 88.906 > 7 111 Y 0.66666666666667 0.00000000000000 0.00000000000000 88.906 > 7 112 Y 0.00000000000000 0.33333333333333 0.00000000000000 88.906 > 7 113 Y 0.33333333333333 0.33333333333333 0.00000000000000 88.906 > 7 114 Y 0.66666666666667 0.33333333333333 0.00000000000000 88.906 > 7 115 Y 0.00000000000000 0.66666666666667 0.00000000000000 88.906 > 7 116 Y 0.33333333333333 0.66666666666667 0.00000000000000 88.906 > 7 117 Y 0.66666666666667 0.66666666666667 0.00000000000000 88.906 > 7 118 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 > 7 119 Y 0.33333333333333 0.00000000000000 0.50000000000000 88.906 > 7 120 Y 0.66666666666667 0.00000000000000 0.50000000000000 88.906 > 7 121 Y 0.00000000000000 0.33333333333333 0.50000000000000 88.906 > 7 122 Y 0.33333333333333 0.33333333333333 0.50000000000000 88.906 > 7 123 Y 0.66666666666667 0.33333333333333 0.50000000000000 88.906 > 7 124 Y 0.00000000000000 0.66666666666667 0.50000000000000 88.906 > 7 125 Y 0.33333333333333 0.66666666666667 0.50000000000000 88.906 > 7 126 Y 0.66666666666667 0.66666666666667 0.50000000000000 88.906 > 7 *127 I 0.00000000000000 0.00000000000000 0.25000000000000 126.904 > 8 128 I 0.33333333333333 0.00000000000000 0.25000000000000 126.904 > 8 129 I 0.66666666666667 0.00000000000000 0.25000000000000 126.904 > 8 130 I 0.00000000000000 0.33333333333333 0.25000000000000 126.904 > 8 131 I 0.33333333333333 0.33333333333333 0.25000000000000 126.904 > 8 132 I 0.66666666666667 0.33333333333333 0.25000000000000 126.904 > 8 133 I 0.00000000000000 0.66666666666667 0.25000000000000 126.904 > 8 134 I 0.33333333333333 0.66666666666667 0.25000000000000 126.904 > 8 135 I 0.66666666666667 0.66666666666667 0.25000000000000 126.904 > 8 136 I 0.00000000000000 0.00000000000000 0.75000000000000 126.904 > 8 137 I 0.33333333333333 0.00000000000000 0.75000000000000 126.904 > 8 138 I 0.66666666666667 0.00000000000000 0.75000000000000 126.904 > 8 139 I 0.00000000000000 0.33333333333333 0.75000000000000 126.904 > 8 140 I 0.33333333333333 0.33333333333333 0.75000000000000 126.904 > 8 141 I 0.66666666666667 0.33333333333333 0.75000000000000 126.904 > 8 142 I 0.00000000000000 0.66666666666667 0.75000000000000 126.904 > 8 143 I 0.33333333333333 0.66666666666667 0.75000000000000 126.904 > 8 144 I 0.66666666666667 0.66666666666667 0.75000000000000 126.904 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.7021339 0.0000000 0.0000000 0.0000000 7.7021339 0.0000000 0.0000000 0.0000000 6.1468052 -------------------------- Born effective charges -------------------------- 1 Bi 5.4072979 0.0000000 0.0000000 0.0000000 5.4072979 0.0000000 0.0000000 0.0000000 2.9131594 2 Bi 5.4072979 0.0000000 0.0000000 0.0000000 5.4072979 0.0000000 0.0000000 0.0000000 2.9131594 3 O -4.0362949 0.0000000 0.0000000 0.0000000 -2.1634310 0.0000000 0.0000000 0.0000000 -2.4666391 4 O -4.0362949 0.0000000 0.0000000 0.0000000 -2.1634310 0.0000000 0.0000000 0.0000000 -2.4666391 5 O -2.1634310 0.0000000 0.0000000 0.0000000 -4.0362949 0.0000000 0.0000000 0.0000000 -2.4666391 6 O -2.1634310 0.0000000 0.0000000 0.0000000 -4.0362949 0.0000000 0.0000000 0.0000000 -2.4666391 7 Y 3.7163555 0.0000000 0.0000000 0.0000000 3.7163555 0.0000000 0.0000000 0.0000000 5.8539902 8 I -2.1314996 0.0000000 0.0000000 0.0000000 -2.1314996 0.0000000 0.0000000 0.0000000 -1.8137529 ---------------------------------------------------------------------------- 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: 5784/5784 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 428 Number of blocks in projector: 428 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 8 --- Eigsh_solver_block: 1 / 8 --- Block_size: 164 Use standard eigh solver. --- Eigsh_solver_block: 2 / 8 --- Block_size: 124 Use standard eigh solver. --- Eigsh_solver_block: 3 / 8 --- Block_size: 99 Use standard eigh solver. --- Eigsh_solver_block: 4 / 8 --- Block_size: 37 Use standard eigh solver. --- Eigsh_solver_block: 5 / 8 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 6 / 8 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 7 / 8 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 8 / 8 --- Block_size: 1 Use standard eigh solver. Tree of FC basis block matrices: - (428, 419), data: False |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (37, 37), data: True |-- (99, 95), data: True |-- (124, 124), data: True |-- (164, 159), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.008 Solver_block: 80 / 80 - Time: 0.359 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.370 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 432/432 Permutation basis: 5784/5784 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 428 Number of blocks in projector: 428 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 8 --- Eigsh_solver_block: 1 / 8 --- Block_size: 164 Use standard eigh solver. --- Eigsh_solver_block: 2 / 8 --- Block_size: 124 Use standard eigh solver. --- Eigsh_solver_block: 3 / 8 --- Block_size: 99 Use standard eigh solver. --- Eigsh_solver_block: 4 / 8 --- Block_size: 37 Use standard eigh solver. --- Eigsh_solver_block: 5 / 8 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 6 / 8 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 7 / 8 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 8 / 8 --- Block_size: 1 Use standard eigh solver. Tree of FC basis block matrices: - (428, 419), data: False |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (37, 37), data: True |-- (99, 95), data: True |-- (124, 124), data: True |-- (164, 159), 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 10:50:39]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 10:50:39]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [3 3 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P4/mmm (123) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.895413370000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.895413370000000 0.000000000000000 c 0.000000000000000 0.000000000000000 9.483492630000001 Atomic positions (fractional): 1 Bi 0.50000000000000 0.50000000000000 0.73986464999698 208.980 2 Bi 0.50000000000000 0.50000000000000 0.26013535000302 208.980 3 O 0.00000000000000 0.50000000000000 0.85671206518205 15.999 4 O 0.00000000000000 0.50000000000000 0.14328793481795 15.999 5 O 0.50000000000000 0.00000000000000 0.85671206518205 15.999 6 O 0.50000000000000 0.00000000000000 0.14328793481795 15.999 7 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 8 I 0.00000000000000 0.00000000000000 0.50000000000000 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.686240110000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.686240110000000 0.000000000000000 c 0.000000000000000 0.000000000000000 18.966985260000001 Atomic positions (fractional): 1 Bi 0.16666666666667 0.16666666666667 0.36993232499849 208.980 > 1 2 Bi 0.50000000000000 0.16666666666667 0.36993232499849 208.980 > 1 3 Bi 0.83333333333333 0.16666666666667 0.36993232499849 208.980 > 1 4 Bi 0.16666666666667 0.50000000000000 0.36993232499849 208.980 > 1 5 Bi 0.50000000000000 0.50000000000000 0.36993232499849 208.980 > 1 6 Bi 0.83333333333333 0.50000000000000 0.36993232499849 208.980 > 1 7 Bi 0.16666666666667 0.83333333333333 0.36993232499849 208.980 > 1 8 Bi 0.50000000000000 0.83333333333333 0.36993232499849 208.980 > 1 9 Bi 0.83333333333333 0.83333333333333 0.36993232499849 208.980 > 1 10 Bi 0.16666666666667 0.16666666666667 0.86993232499849 208.980 > 1 11 Bi 0.50000000000000 0.16666666666667 0.86993232499849 208.980 > 1 12 Bi 0.83333333333333 0.16666666666667 0.86993232499849 208.980 > 1 13 Bi 0.16666666666667 0.50000000000000 0.86993232499849 208.980 > 1 14 Bi 0.50000000000000 0.50000000000000 0.86993232499849 208.980 > 1 15 Bi 0.83333333333333 0.50000000000000 0.86993232499849 208.980 > 1 16 Bi 0.16666666666667 0.83333333333333 0.86993232499849 208.980 > 1 17 Bi 0.50000000000000 0.83333333333333 0.86993232499849 208.980 > 1 18 Bi 0.83333333333333 0.83333333333333 0.86993232499849 208.980 > 1 19 Bi 0.16666666666667 0.16666666666667 0.13006767500151 208.980 > 19 20 Bi 0.50000000000000 0.16666666666667 0.13006767500151 208.980 > 19 21 Bi 0.83333333333333 0.16666666666667 0.13006767500151 208.980 > 19 22 Bi 0.16666666666667 0.50000000000000 0.13006767500151 208.980 > 19 23 Bi 0.50000000000000 0.50000000000000 0.13006767500151 208.980 > 19 24 Bi 0.83333333333333 0.50000000000000 0.13006767500151 208.980 > 19 25 Bi 0.16666666666667 0.83333333333333 0.13006767500151 208.980 > 19 26 Bi 0.50000000000000 0.83333333333333 0.13006767500151 208.980 > 19 27 Bi 0.83333333333333 0.83333333333333 0.13006767500151 208.980 > 19 28 Bi 0.16666666666667 0.16666666666667 0.63006767500151 208.980 > 19 29 Bi 0.50000000000000 0.16666666666667 0.63006767500151 208.980 > 19 30 Bi 0.83333333333333 0.16666666666667 0.63006767500151 208.980 > 19 31 Bi 0.16666666666667 0.50000000000000 0.63006767500151 208.980 > 19 32 Bi 0.50000000000000 0.50000000000000 0.63006767500151 208.980 > 19 33 Bi 0.83333333333333 0.50000000000000 0.63006767500151 208.980 > 19 34 Bi 0.16666666666667 0.83333333333333 0.63006767500151 208.980 > 19 35 Bi 0.50000000000000 0.83333333333333 0.63006767500151 208.980 > 19 36 Bi 0.83333333333333 0.83333333333333 0.63006767500151 208.980 > 19 37 O 0.00000000000000 0.16666666666667 0.42835603259102 15.999 > 37 38 O 0.33333333333333 0.16666666666667 0.42835603259102 15.999 > 37 39 O 0.66666666666667 0.16666666666667 0.42835603259102 15.999 > 37 40 O 0.00000000000000 0.50000000000000 0.42835603259102 15.999 > 37 41 O 0.33333333333333 0.50000000000000 0.42835603259102 15.999 > 37 42 O 0.66666666666667 0.50000000000000 0.42835603259102 15.999 > 37 43 O 0.00000000000000 0.83333333333333 0.42835603259102 15.999 > 37 44 O 0.33333333333333 0.83333333333333 0.42835603259102 15.999 > 37 45 O 0.66666666666667 0.83333333333333 0.42835603259102 15.999 > 37 46 O 0.00000000000000 0.16666666666667 0.92835603259102 15.999 > 37 47 O 0.33333333333333 0.16666666666667 0.92835603259102 15.999 > 37 48 O 0.66666666666667 0.16666666666667 0.92835603259102 15.999 > 37 49 O 0.00000000000000 0.50000000000000 0.92835603259102 15.999 > 37 50 O 0.33333333333333 0.50000000000000 0.92835603259102 15.999 > 37 51 O 0.66666666666667 0.50000000000000 0.92835603259102 15.999 > 37 52 O 0.00000000000000 0.83333333333333 0.92835603259102 15.999 > 37 53 O 0.33333333333333 0.83333333333333 0.92835603259102 15.999 > 37 54 O 0.66666666666667 0.83333333333333 0.92835603259102 15.999 > 37 55 O 0.00000000000000 0.16666666666667 0.07164396740898 15.999 > 55 56 O 0.33333333333333 0.16666666666667 0.07164396740898 15.999 > 55 57 O 0.66666666666667 0.16666666666667 0.07164396740898 15.999 > 55 58 O 0.00000000000000 0.50000000000000 0.07164396740898 15.999 > 55 59 O 0.33333333333333 0.50000000000000 0.07164396740898 15.999 > 55 60 O 0.66666666666667 0.50000000000000 0.07164396740898 15.999 > 55 61 O 0.00000000000000 0.83333333333333 0.07164396740898 15.999 > 55 62 O 0.33333333333333 0.83333333333333 0.07164396740898 15.999 > 55 63 O 0.66666666666667 0.83333333333333 0.07164396740898 15.999 > 55 64 O 0.00000000000000 0.16666666666667 0.57164396740898 15.999 > 55 65 O 0.33333333333333 0.16666666666667 0.57164396740898 15.999 > 55 66 O 0.66666666666667 0.16666666666667 0.57164396740898 15.999 > 55 67 O 0.00000000000000 0.50000000000000 0.57164396740898 15.999 > 55 68 O 0.33333333333333 0.50000000000000 0.57164396740898 15.999 > 55 69 O 0.66666666666667 0.50000000000000 0.57164396740898 15.999 > 55 70 O 0.00000000000000 0.83333333333333 0.57164396740898 15.999 > 55 71 O 0.33333333333333 0.83333333333333 0.57164396740898 15.999 > 55 72 O 0.66666666666667 0.83333333333333 0.57164396740898 15.999 > 55 73 O 0.16666666666667 0.00000000000000 0.42835603259102 15.999 > 73 74 O 0.50000000000000 0.00000000000000 0.42835603259102 15.999 > 73 75 O 0.83333333333333 0.00000000000000 0.42835603259102 15.999 > 73 76 O 0.16666666666667 0.33333333333333 0.42835603259102 15.999 > 73 77 O 0.50000000000000 0.33333333333333 0.42835603259102 15.999 > 73 78 O 0.83333333333333 0.33333333333333 0.42835603259102 15.999 > 73 79 O 0.16666666666667 0.66666666666667 0.42835603259102 15.999 > 73 80 O 0.50000000000000 0.66666666666667 0.42835603259102 15.999 > 73 81 O 0.83333333333333 0.66666666666667 0.42835603259102 15.999 > 73 82 O 0.16666666666667 0.00000000000000 0.92835603259102 15.999 > 73 83 O 0.50000000000000 0.00000000000000 0.92835603259102 15.999 > 73 84 O 0.83333333333333 0.00000000000000 0.92835603259102 15.999 > 73 85 O 0.16666666666667 0.33333333333333 0.92835603259102 15.999 > 73 86 O 0.50000000000000 0.33333333333333 0.92835603259102 15.999 > 73 87 O 0.83333333333333 0.33333333333333 0.92835603259102 15.999 > 73 88 O 0.16666666666667 0.66666666666667 0.92835603259102 15.999 > 73 89 O 0.50000000000000 0.66666666666667 0.92835603259102 15.999 > 73 90 O 0.83333333333333 0.66666666666667 0.92835603259102 15.999 > 73 91 O 0.16666666666667 0.00000000000000 0.07164396740898 15.999 > 91 92 O 0.50000000000000 0.00000000000000 0.07164396740898 15.999 > 91 93 O 0.83333333333333 0.00000000000000 0.07164396740898 15.999 > 91 94 O 0.16666666666667 0.33333333333333 0.07164396740898 15.999 > 91 95 O 0.50000000000000 0.33333333333333 0.07164396740898 15.999 > 91 96 O 0.83333333333333 0.33333333333333 0.07164396740898 15.999 > 91 97 O 0.16666666666667 0.66666666666667 0.07164396740898 15.999 > 91 98 O 0.50000000000000 0.66666666666667 0.07164396740898 15.999 > 91 99 O 0.83333333333333 0.66666666666667 0.07164396740898 15.999 > 91 100 O 0.16666666666667 0.00000000000000 0.57164396740898 15.999 > 91 101 O 0.50000000000000 0.00000000000000 0.57164396740898 15.999 > 91 102 O 0.83333333333333 0.00000000000000 0.57164396740898 15.999 > 91 103 O 0.16666666666667 0.33333333333333 0.57164396740898 15.999 > 91 104 O 0.50000000000000 0.33333333333333 0.57164396740898 15.999 > 91 105 O 0.83333333333333 0.33333333333333 0.57164396740898 15.999 > 91 106 O 0.16666666666667 0.66666666666667 0.57164396740898 15.999 > 91 107 O 0.50000000000000 0.66666666666667 0.57164396740898 15.999 > 91 108 O 0.83333333333333 0.66666666666667 0.57164396740898 15.999 > 91 109 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 > 109 110 Y 0.33333333333333 0.00000000000000 0.00000000000000 88.906 > 109 111 Y 0.66666666666667 0.00000000000000 0.00000000000000 88.906 > 109 112 Y 0.00000000000000 0.33333333333333 0.00000000000000 88.906 > 109 113 Y 0.33333333333333 0.33333333333333 0.00000000000000 88.906 > 109 114 Y 0.66666666666667 0.33333333333333 0.00000000000000 88.906 > 109 115 Y 0.00000000000000 0.66666666666667 0.00000000000000 88.906 > 109 116 Y 0.33333333333333 0.66666666666667 0.00000000000000 88.906 > 109 117 Y 0.66666666666667 0.66666666666667 0.00000000000000 88.906 > 109 118 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 > 109 119 Y 0.33333333333333 0.00000000000000 0.50000000000000 88.906 > 109 120 Y 0.66666666666667 0.00000000000000 0.50000000000000 88.906 > 109 121 Y 0.00000000000000 0.33333333333333 0.50000000000000 88.906 > 109 122 Y 0.33333333333333 0.33333333333333 0.50000000000000 88.906 > 109 123 Y 0.66666666666667 0.33333333333333 0.50000000000000 88.906 > 109 124 Y 0.00000000000000 0.66666666666667 0.50000000000000 88.906 > 109 125 Y 0.33333333333333 0.66666666666667 0.50000000000000 88.906 > 109 126 Y 0.66666666666667 0.66666666666667 0.50000000000000 88.906 > 109 127 I 0.00000000000000 0.00000000000000 0.25000000000000 126.904 > 127 128 I 0.33333333333333 0.00000000000000 0.25000000000000 126.904 > 127 129 I 0.66666666666667 0.00000000000000 0.25000000000000 126.904 > 127 130 I 0.00000000000000 0.33333333333333 0.25000000000000 126.904 > 127 131 I 0.33333333333333 0.33333333333333 0.25000000000000 126.904 > 127 132 I 0.66666666666667 0.33333333333333 0.25000000000000 126.904 > 127 133 I 0.00000000000000 0.66666666666667 0.25000000000000 126.904 > 127 134 I 0.33333333333333 0.66666666666667 0.25000000000000 126.904 > 127 135 I 0.66666666666667 0.66666666666667 0.25000000000000 126.904 > 127 136 I 0.00000000000000 0.00000000000000 0.75000000000000 126.904 > 127 137 I 0.33333333333333 0.00000000000000 0.75000000000000 126.904 > 127 138 I 0.66666666666667 0.00000000000000 0.75000000000000 126.904 > 127 139 I 0.00000000000000 0.33333333333333 0.75000000000000 126.904 > 127 140 I 0.33333333333333 0.33333333333333 0.75000000000000 126.904 > 127 141 I 0.66666666666667 0.33333333333333 0.75000000000000 126.904 > 127 142 I 0.00000000000000 0.66666666666667 0.75000000000000 126.904 > 127 143 I 0.33333333333333 0.66666666666667 0.75000000000000 126.904 > 127 144 I 0.66666666666667 0.66666666666667 0.75000000000000 126.904 > 127 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.7021339 0.0000000 0.0000000 0.0000000 7.7021339 0.0000000 0.0000000 0.0000000 6.1468052 -------------------------- Born effective charges -------------------------- 1 Bi 5.4072979 0.0000000 0.0000000 0.0000000 5.4072979 0.0000000 0.0000000 0.0000000 2.9131594 2 Bi 5.4072979 0.0000000 0.0000000 0.0000000 5.4072979 0.0000000 0.0000000 0.0000000 2.9131594 3 O -4.0362949 0.0000000 0.0000000 0.0000000 -2.1634310 0.0000000 0.0000000 0.0000000 -2.4666391 4 O -4.0362949 0.0000000 0.0000000 0.0000000 -2.1634310 0.0000000 0.0000000 0.0000000 -2.4666391 5 O -2.1634310 0.0000000 0.0000000 0.0000000 -4.0362949 0.0000000 0.0000000 0.0000000 -2.4666391 6 O -2.1634310 0.0000000 0.0000000 0.0000000 -4.0362949 0.0000000 0.0000000 0.0000000 -2.4666391 7 Y 3.7163555 0.0000000 0.0000000 0.0000000 3.7163555 0.0000000 0.0000000 0.0000000 5.8539902 8 I -2.1314996 0.0000000 0.0000000 0.0000000 -2.1314996 0.0000000 0.0000000 0.0000000 -1.8137529 ---------------------------------------------------------------------------- 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[ 37, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 109, 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[ 127, 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.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (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 10:50:49]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 10:50:49]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [3 3 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P4/mmm (123) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.895413370000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.895413370000000 0.000000000000000 c 0.000000000000000 0.000000000000000 9.483492630000001 Atomic positions (fractional): 1 Bi 0.50000000000000 0.50000000000000 0.73986464999698 208.980 2 Bi 0.50000000000000 0.50000000000000 0.26013535000302 208.980 3 O 0.00000000000000 0.50000000000000 0.85671206518205 15.999 4 O 0.00000000000000 0.50000000000000 0.14328793481795 15.999 5 O 0.50000000000000 0.00000000000000 0.85671206518205 15.999 6 O 0.50000000000000 0.00000000000000 0.14328793481795 15.999 7 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 8 I 0.00000000000000 0.00000000000000 0.50000000000000 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.686240110000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.686240110000000 0.000000000000000 c 0.000000000000000 0.000000000000000 18.966985260000001 Atomic positions (fractional): 1 Bi 0.16666666666667 0.16666666666667 0.36993232499849 208.980 > 1 2 Bi 0.50000000000000 0.16666666666667 0.36993232499849 208.980 > 1 3 Bi 0.83333333333333 0.16666666666667 0.36993232499849 208.980 > 1 4 Bi 0.16666666666667 0.50000000000000 0.36993232499849 208.980 > 1 5 Bi 0.50000000000000 0.50000000000000 0.36993232499849 208.980 > 1 6 Bi 0.83333333333333 0.50000000000000 0.36993232499849 208.980 > 1 7 Bi 0.16666666666667 0.83333333333333 0.36993232499849 208.980 > 1 8 Bi 0.50000000000000 0.83333333333333 0.36993232499849 208.980 > 1 9 Bi 0.83333333333333 0.83333333333333 0.36993232499849 208.980 > 1 10 Bi 0.16666666666667 0.16666666666667 0.86993232499849 208.980 > 1 11 Bi 0.50000000000000 0.16666666666667 0.86993232499849 208.980 > 1 12 Bi 0.83333333333333 0.16666666666667 0.86993232499849 208.980 > 1 13 Bi 0.16666666666667 0.50000000000000 0.86993232499849 208.980 > 1 14 Bi 0.50000000000000 0.50000000000000 0.86993232499849 208.980 > 1 15 Bi 0.83333333333333 0.50000000000000 0.86993232499849 208.980 > 1 16 Bi 0.16666666666667 0.83333333333333 0.86993232499849 208.980 > 1 17 Bi 0.50000000000000 0.83333333333333 0.86993232499849 208.980 > 1 18 Bi 0.83333333333333 0.83333333333333 0.86993232499849 208.980 > 1 19 Bi 0.16666666666667 0.16666666666667 0.13006767500151 208.980 > 19 20 Bi 0.50000000000000 0.16666666666667 0.13006767500151 208.980 > 19 21 Bi 0.83333333333333 0.16666666666667 0.13006767500151 208.980 > 19 22 Bi 0.16666666666667 0.50000000000000 0.13006767500151 208.980 > 19 23 Bi 0.50000000000000 0.50000000000000 0.13006767500151 208.980 > 19 24 Bi 0.83333333333333 0.50000000000000 0.13006767500151 208.980 > 19 25 Bi 0.16666666666667 0.83333333333333 0.13006767500151 208.980 > 19 26 Bi 0.50000000000000 0.83333333333333 0.13006767500151 208.980 > 19 27 Bi 0.83333333333333 0.83333333333333 0.13006767500151 208.980 > 19 28 Bi 0.16666666666667 0.16666666666667 0.63006767500151 208.980 > 19 29 Bi 0.50000000000000 0.16666666666667 0.63006767500151 208.980 > 19 30 Bi 0.83333333333333 0.16666666666667 0.63006767500151 208.980 > 19 31 Bi 0.16666666666667 0.50000000000000 0.63006767500151 208.980 > 19 32 Bi 0.50000000000000 0.50000000000000 0.63006767500151 208.980 > 19 33 Bi 0.83333333333333 0.50000000000000 0.63006767500151 208.980 > 19 34 Bi 0.16666666666667 0.83333333333333 0.63006767500151 208.980 > 19 35 Bi 0.50000000000000 0.83333333333333 0.63006767500151 208.980 > 19 36 Bi 0.83333333333333 0.83333333333333 0.63006767500151 208.980 > 19 37 O 0.00000000000000 0.16666666666667 0.42835603259102 15.999 > 37 38 O 0.33333333333333 0.16666666666667 0.42835603259102 15.999 > 37 39 O 0.66666666666667 0.16666666666667 0.42835603259102 15.999 > 37 40 O 0.00000000000000 0.50000000000000 0.42835603259102 15.999 > 37 41 O 0.33333333333333 0.50000000000000 0.42835603259102 15.999 > 37 42 O 0.66666666666667 0.50000000000000 0.42835603259102 15.999 > 37 43 O 0.00000000000000 0.83333333333333 0.42835603259102 15.999 > 37 44 O 0.33333333333333 0.83333333333333 0.42835603259102 15.999 > 37 45 O 0.66666666666667 0.83333333333333 0.42835603259102 15.999 > 37 46 O 0.00000000000000 0.16666666666667 0.92835603259102 15.999 > 37 47 O 0.33333333333333 0.16666666666667 0.92835603259102 15.999 > 37 48 O 0.66666666666667 0.16666666666667 0.92835603259102 15.999 > 37 49 O 0.00000000000000 0.50000000000000 0.92835603259102 15.999 > 37 50 O 0.33333333333333 0.50000000000000 0.92835603259102 15.999 > 37 51 O 0.66666666666667 0.50000000000000 0.92835603259102 15.999 > 37 52 O 0.00000000000000 0.83333333333333 0.92835603259102 15.999 > 37 53 O 0.33333333333333 0.83333333333333 0.92835603259102 15.999 > 37 54 O 0.66666666666667 0.83333333333333 0.92835603259102 15.999 > 37 55 O 0.00000000000000 0.16666666666667 0.07164396740898 15.999 > 55 56 O 0.33333333333333 0.16666666666667 0.07164396740898 15.999 > 55 57 O 0.66666666666667 0.16666666666667 0.07164396740898 15.999 > 55 58 O 0.00000000000000 0.50000000000000 0.07164396740898 15.999 > 55 59 O 0.33333333333333 0.50000000000000 0.07164396740898 15.999 > 55 60 O 0.66666666666667 0.50000000000000 0.07164396740898 15.999 > 55 61 O 0.00000000000000 0.83333333333333 0.07164396740898 15.999 > 55 62 O 0.33333333333333 0.83333333333333 0.07164396740898 15.999 > 55 63 O 0.66666666666667 0.83333333333333 0.07164396740898 15.999 > 55 64 O 0.00000000000000 0.16666666666667 0.57164396740898 15.999 > 55 65 O 0.33333333333333 0.16666666666667 0.57164396740898 15.999 > 55 66 O 0.66666666666667 0.16666666666667 0.57164396740898 15.999 > 55 67 O 0.00000000000000 0.50000000000000 0.57164396740898 15.999 > 55 68 O 0.33333333333333 0.50000000000000 0.57164396740898 15.999 > 55 69 O 0.66666666666667 0.50000000000000 0.57164396740898 15.999 > 55 70 O 0.00000000000000 0.83333333333333 0.57164396740898 15.999 > 55 71 O 0.33333333333333 0.83333333333333 0.57164396740898 15.999 > 55 72 O 0.66666666666667 0.83333333333333 0.57164396740898 15.999 > 55 73 O 0.16666666666667 0.00000000000000 0.42835603259102 15.999 > 73 74 O 0.50000000000000 0.00000000000000 0.42835603259102 15.999 > 73 75 O 0.83333333333333 0.00000000000000 0.42835603259102 15.999 > 73 76 O 0.16666666666667 0.33333333333333 0.42835603259102 15.999 > 73 77 O 0.50000000000000 0.33333333333333 0.42835603259102 15.999 > 73 78 O 0.83333333333333 0.33333333333333 0.42835603259102 15.999 > 73 79 O 0.16666666666667 0.66666666666667 0.42835603259102 15.999 > 73 80 O 0.50000000000000 0.66666666666667 0.42835603259102 15.999 > 73 81 O 0.83333333333333 0.66666666666667 0.42835603259102 15.999 > 73 82 O 0.16666666666667 0.00000000000000 0.92835603259102 15.999 > 73 83 O 0.50000000000000 0.00000000000000 0.92835603259102 15.999 > 73 84 O 0.83333333333333 0.00000000000000 0.92835603259102 15.999 > 73 85 O 0.16666666666667 0.33333333333333 0.92835603259102 15.999 > 73 86 O 0.50000000000000 0.33333333333333 0.92835603259102 15.999 > 73 87 O 0.83333333333333 0.33333333333333 0.92835603259102 15.999 > 73 88 O 0.16666666666667 0.66666666666667 0.92835603259102 15.999 > 73 89 O 0.50000000000000 0.66666666666667 0.92835603259102 15.999 > 73 90 O 0.83333333333333 0.66666666666667 0.92835603259102 15.999 > 73 91 O 0.16666666666667 0.00000000000000 0.07164396740898 15.999 > 91 92 O 0.50000000000000 0.00000000000000 0.07164396740898 15.999 > 91 93 O 0.83333333333333 0.00000000000000 0.07164396740898 15.999 > 91 94 O 0.16666666666667 0.33333333333333 0.07164396740898 15.999 > 91 95 O 0.50000000000000 0.33333333333333 0.07164396740898 15.999 > 91 96 O 0.83333333333333 0.33333333333333 0.07164396740898 15.999 > 91 97 O 0.16666666666667 0.66666666666667 0.07164396740898 15.999 > 91 98 O 0.50000000000000 0.66666666666667 0.07164396740898 15.999 > 91 99 O 0.83333333333333 0.66666666666667 0.07164396740898 15.999 > 91 100 O 0.16666666666667 0.00000000000000 0.57164396740898 15.999 > 91 101 O 0.50000000000000 0.00000000000000 0.57164396740898 15.999 > 91 102 O 0.83333333333333 0.00000000000000 0.57164396740898 15.999 > 91 103 O 0.16666666666667 0.33333333333333 0.57164396740898 15.999 > 91 104 O 0.50000000000000 0.33333333333333 0.57164396740898 15.999 > 91 105 O 0.83333333333333 0.33333333333333 0.57164396740898 15.999 > 91 106 O 0.16666666666667 0.66666666666667 0.57164396740898 15.999 > 91 107 O 0.50000000000000 0.66666666666667 0.57164396740898 15.999 > 91 108 O 0.83333333333333 0.66666666666667 0.57164396740898 15.999 > 91 109 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 > 109 110 Y 0.33333333333333 0.00000000000000 0.00000000000000 88.906 > 109 111 Y 0.66666666666667 0.00000000000000 0.00000000000000 88.906 > 109 112 Y 0.00000000000000 0.33333333333333 0.00000000000000 88.906 > 109 113 Y 0.33333333333333 0.33333333333333 0.00000000000000 88.906 > 109 114 Y 0.66666666666667 0.33333333333333 0.00000000000000 88.906 > 109 115 Y 0.00000000000000 0.66666666666667 0.00000000000000 88.906 > 109 116 Y 0.33333333333333 0.66666666666667 0.00000000000000 88.906 > 109 117 Y 0.66666666666667 0.66666666666667 0.00000000000000 88.906 > 109 118 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 > 109 119 Y 0.33333333333333 0.00000000000000 0.50000000000000 88.906 > 109 120 Y 0.66666666666667 0.00000000000000 0.50000000000000 88.906 > 109 121 Y 0.00000000000000 0.33333333333333 0.50000000000000 88.906 > 109 122 Y 0.33333333333333 0.33333333333333 0.50000000000000 88.906 > 109 123 Y 0.66666666666667 0.33333333333333 0.50000000000000 88.906 > 109 124 Y 0.00000000000000 0.66666666666667 0.50000000000000 88.906 > 109 125 Y 0.33333333333333 0.66666666666667 0.50000000000000 88.906 > 109 126 Y 0.66666666666667 0.66666666666667 0.50000000000000 88.906 > 109 127 I 0.00000000000000 0.00000000000000 0.25000000000000 126.904 > 127 128 I 0.33333333333333 0.00000000000000 0.25000000000000 126.904 > 127 129 I 0.66666666666667 0.00000000000000 0.25000000000000 126.904 > 127 130 I 0.00000000000000 0.33333333333333 0.25000000000000 126.904 > 127 131 I 0.33333333333333 0.33333333333333 0.25000000000000 126.904 > 127 132 I 0.66666666666667 0.33333333333333 0.25000000000000 126.904 > 127 133 I 0.00000000000000 0.66666666666667 0.25000000000000 126.904 > 127 134 I 0.33333333333333 0.66666666666667 0.25000000000000 126.904 > 127 135 I 0.66666666666667 0.66666666666667 0.25000000000000 126.904 > 127 136 I 0.00000000000000 0.00000000000000 0.75000000000000 126.904 > 127 137 I 0.33333333333333 0.00000000000000 0.75000000000000 126.904 > 127 138 I 0.66666666666667 0.00000000000000 0.75000000000000 126.904 > 127 139 I 0.00000000000000 0.33333333333333 0.75000000000000 126.904 > 127 140 I 0.33333333333333 0.33333333333333 0.75000000000000 126.904 > 127 141 I 0.66666666666667 0.33333333333333 0.75000000000000 126.904 > 127 142 I 0.00000000000000 0.66666666666667 0.75000000000000 126.904 > 127 143 I 0.33333333333333 0.66666666666667 0.75000000000000 126.904 > 127 144 I 0.66666666666667 0.66666666666667 0.75000000000000 126.904 > 127 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 7.7021339 0.0000000 0.0000000 0.0000000 7.7021339 0.0000000 0.0000000 0.0000000 6.1468052 -------------------------- Born effective charges -------------------------- 1 Bi 5.4072979 0.0000000 0.0000000 0.0000000 5.4072979 0.0000000 0.0000000 0.0000000 2.9131594 2 Bi 5.4072979 0.0000000 0.0000000 0.0000000 5.4072979 0.0000000 0.0000000 0.0000000 2.9131594 3 O -4.0362949 0.0000000 0.0000000 0.0000000 -2.1634310 0.0000000 0.0000000 0.0000000 -2.4666391 4 O -4.0362949 0.0000000 0.0000000 0.0000000 -2.1634310 0.0000000 0.0000000 0.0000000 -2.4666391 5 O -2.1634310 0.0000000 0.0000000 0.0000000 -4.0362949 0.0000000 0.0000000 0.0000000 -2.4666391 6 O -2.1634310 0.0000000 0.0000000 0.0000000 -4.0362949 0.0000000 0.0000000 0.0000000 -2.4666391 7 Y 3.7163555 0.0000000 0.0000000 0.0000000 3.7163555 0.0000000 0.0000000 0.0000000 5.8539902 8 I -2.1314996 0.0000000 0.0000000 0.0000000 -2.1314996 0.0000000 0.0000000 0.0000000 -1.8137529 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (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 ... [ 13 13 5 ] 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.79, Number of G-points: 299, Lambda: 0.22 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/84) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 84 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 0.000 ( 0.000 0.000 0.000) 0.000 1.760 ( 0.000 0.000 0.000) 0.000 1.760 ( 0.000 0.000 0.000) 0.000 1.874 ( 0.000 0.000 0.000) 0.000 1.874 ( 0.000 0.000 0.000) 0.000 2.181 ( 0.000 0.000 0.000) 0.000 3.614 ( 0.000 0.000 0.000) 0.000 3.895 ( 0.000 0.000 0.000) 0.000 3.895 ( 0.000 0.000 0.000) 0.000 6.873 ( 0.000 0.000 0.000) 0.000 8.289 ( 0.000 0.000 0.000) 0.000 8.289 ( 0.000 0.000 0.000) 0.000 8.449 ( 0.000 0.000 0.000) 0.000 8.449 ( 0.000 0.000 0.000) 0.000 9.158 ( 0.000 0.000 0.000) 0.000 11.876 ( 0.000 0.000 0.000) 0.000 12.656 ( 0.000 0.000 0.000) 0.000 13.430 ( 0.000 0.000 0.000) 0.000 13.430 ( 0.000 0.000 0.000) 0.000 13.604 ( 0.000 0.000 0.000) 0.000 13.604 ( 0.000 0.000 0.000) 0.000 15.768 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/84) ======================= q-point: ( 0.08 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.330 ( 18.166 0.000 0.000) 18.166 0.605 ( 28.983 0.000 0.000) 28.983 0.888 ( 43.602 0.000 0.000) 43.602 1.801 ( 3.895 0.000 0.000) 3.895 1.824 ( 5.058 0.000 0.000) 5.058 1.880 ( 0.680 0.000 0.000) 0.680 2.230 ( 5.575 0.000 0.000) 5.575 2.273 ( 6.520 0.000 0.000) 6.520 3.511 ( -6.954 0.000 0.000) 6.954 3.900 ( 0.538 0.000 0.000) 0.538 4.166 ( 18.152 0.000 0.000) 18.152 6.672 ( -17.268 0.000 0.000) 17.268 8.403 ( 11.369 0.000 0.000) 11.369 8.574 ( 12.465 0.000 0.000) 12.465 8.637 ( 14.432 0.000 0.000) 14.432 9.131 ( -1.113 0.000 0.000) 1.113 11.514 ( 3.769 0.000 0.000) 3.769 11.831 ( -4.226 0.000 0.000) 4.226 12.639 ( -1.560 0.000 0.000) 1.560 13.389 ( -4.234 0.000 0.000) 4.234 13.533 ( 9.669 0.000 0.000) 9.669 13.578 ( -2.598 0.000 0.000) 2.598 15.510 ( -8.238 0.000 0.000) 8.238 15.789 ( 0.220 0.000 0.000) 0.220 ======================= Grid point 2 (3/84) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.723 ( 21.409 0.000 0.000) 21.409 1.117 ( 22.410 0.000 0.000) 22.410 1.691 ( 36.795 0.000 0.000) 36.795 1.897 ( 1.766 0.000 0.000) 1.766 1.899 ( 5.572 0.000 0.000) 5.572 1.904 ( 1.835 0.000 0.000) 1.835 2.422 ( 13.381 0.000 0.000) 13.381 2.441 ( 9.769 0.000 0.000) 9.769 3.400 ( -3.518 0.000 0.000) 3.518 3.914 ( 0.812 0.000 0.000) 0.812 4.596 ( 24.486 0.000 0.000) 24.486 6.305 ( -18.038 0.000 0.000) 18.038 8.725 ( 20.762 0.000 0.000) 20.762 8.737 ( -7.175 0.000 0.000) 7.175 8.927 ( 22.695 0.000 0.000) 22.695 9.293 ( 18.878 0.000 0.000) 18.878 11.611 ( 5.674 0.000 0.000) 5.674 11.726 ( -6.063 0.000 0.000) 6.063 12.609 ( -0.927 0.000 0.000) 0.927 13.264 ( -8.431 0.000 0.000) 8.431 13.502 ( -5.113 0.000 0.000) 5.113 13.773 ( 13.647 0.000 0.000) 13.647 15.326 ( -9.963 0.000 0.000) 9.963 15.743 ( -4.951 0.000 0.000) 4.951 ======================= Grid point 3 (4/84) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.161 ( 22.658 0.000 0.000) 22.658 1.473 ( 13.352 0.000 0.000) 13.352 1.907 ( -0.147 0.000 0.000) 0.147 1.961 ( 4.206 0.000 0.000) 4.206 2.007 ( 5.106 0.000 0.000) 5.106 2.290 ( 22.255 0.000 0.000) 22.255 2.630 ( 8.630 0.000 0.000) 8.630 2.705 ( 13.839 0.000 0.000) 13.839 3.419 ( 7.048 0.000 0.000) 7.048 3.930 ( 0.744 0.000 0.000) 0.744 5.102 ( 25.989 0.000 0.000) 25.989 6.001 ( -12.341 0.000 0.000) 12.341 8.489 ( -15.557 0.000 0.000) 15.557 9.197 ( 26.234 0.000 0.000) 26.234 9.446 ( 29.226 0.000 0.000) 29.226 9.720 ( 21.967 0.000 0.000) 21.967 11.600 ( -6.604 0.000 0.000) 6.604 11.717 ( 4.552 0.000 0.000) 4.552 12.624 ( 3.001 0.000 0.000) 3.001 13.056 ( -12.562 0.000 0.000) 12.562 13.378 ( -7.361 0.000 0.000) 7.361 14.038 ( 12.512 0.000 0.000) 12.512 15.123 ( -10.679 0.000 0.000) 10.679 15.600 ( -9.216 0.000 0.000) 9.216 ======================= Grid point 4 (5/84) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.605 ( 21.894 0.000 0.000) 21.894 1.651 ( 5.349 0.000 0.000) 5.349 1.909 ( 0.717 0.000 0.000) 0.717 2.067 ( 6.050 0.000 0.000) 6.050 2.094 ( 3.679 0.000 0.000) 3.679 2.540 ( 4.081 0.000 0.000) 4.081 2.782 ( 7.400 0.000 0.000) 7.400 2.926 ( 7.752 0.000 0.000) 7.752 3.688 ( 18.147 0.000 0.000) 18.147 3.942 ( 0.479 0.000 0.000) 0.479 5.592 ( 22.861 0.000 0.000) 22.861 5.825 ( -5.624 0.000 0.000) 5.624 8.174 ( -15.381 0.000 0.000) 15.381 9.726 ( 26.425 0.000 0.000) 26.425 10.054 ( 31.699 0.000 0.000) 31.699 10.116 ( 17.518 0.000 0.000) 17.518 11.472 ( -6.147 0.000 0.000) 6.147 11.781 ( 2.031 0.000 0.000) 2.031 12.734 ( 7.776 0.000 0.000) 7.776 12.768 ( -16.536 0.000 0.000) 16.536 13.216 ( -8.818 0.000 0.000) 8.818 14.252 ( 9.007 0.000 0.000) 9.007 14.904 ( -11.350 0.000 0.000) 11.350 15.400 ( -10.454 0.000 0.000) 10.454 ======================= Grid point 5 (6/84) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.716 ( 1.818 0.000 0.000) 1.818 1.888 ( 2.580 0.000 0.000) 2.580 2.064 ( 18.569 0.000 0.000) 18.569 2.151 ( 2.116 0.000 0.000) 2.116 2.176 ( 4.487 0.000 0.000) 4.487 2.545 ( -1.601 0.000 0.000) 1.601 2.911 ( 5.154 0.000 0.000) 5.154 2.995 ( -0.798 0.000 0.000) 0.798 3.949 ( 0.219 0.000 0.000) 0.219 4.038 ( 15.547 0.000 0.000) 15.547 5.763 ( -1.189 0.000 0.000) 1.189 5.978 ( 15.623 0.000 0.000) 15.623 7.911 ( -10.713 0.000 0.000) 10.713 10.200 ( 20.530 0.000 0.000) 20.530 10.396 ( 10.690 0.000 0.000) 10.690 10.672 ( 30.220 0.000 0.000) 30.220 11.366 ( -4.360 0.000 0.000) 4.360 11.806 ( 0.683 0.000 0.000) 0.683 12.407 ( -19.769 0.000 0.000) 19.769 12.907 ( 8.968 0.000 0.000) 8.968 13.043 ( -8.202 0.000 0.000) 8.202 14.391 ( 5.153 0.000 0.000) 5.153 14.685 ( -10.384 0.000 0.000) 10.384 15.211 ( -8.135 0.000 0.000) 8.135 ======================= Grid point 6 (7/84) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.737 ( 0.454 0.000 0.000) 0.454 1.910 ( 0.475 0.000 0.000) 0.475 2.179 ( 0.681 0.000 0.000) 0.681 2.236 ( 1.556 0.000 0.000) 1.556 2.394 ( 11.725 0.000 0.000) 11.725 2.516 ( -0.939 0.000 0.000) 0.939 2.919 ( -4.971 0.000 0.000) 4.971 2.979 ( 1.758 0.000 0.000) 1.758 3.952 ( 0.056 0.000 0.000) 0.056 4.256 ( 5.832 0.000 0.000) 5.832 5.756 ( 0.040 0.000 0.000) 0.040 6.190 ( 5.543 0.000 0.000) 5.543 7.766 ( -3.743 0.000 0.000) 3.743 10.493 ( 8.050 0.000 0.000) 8.050 10.536 ( 3.484 0.000 0.000) 3.484 11.199 ( 20.143 0.000 0.000) 20.143 11.306 ( -1.552 0.000 0.000) 1.552 11.814 ( 0.191 0.000 0.000) 0.191 12.023 ( -16.431 0.000 0.000) 16.431 12.919 ( -3.638 0.000 0.000) 3.638 13.048 ( 4.267 0.000 0.000) 4.267 14.458 ( 1.652 0.000 0.000) 1.652 14.526 ( -4.684 0.000 0.000) 4.684 15.097 ( -3.069 0.000 0.000) 3.069 ======================= Grid point 14 (8/84) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.473 ( 13.231 13.231 0.000) 18.711 0.568 ( 14.192 14.192 0.000) 20.071 1.331 ( 28.049 28.049 0.000) 39.667 1.824 ( 2.506 2.506 0.000) 3.545 1.912 ( 7.179 7.179 0.000) 10.153 1.931 ( 2.989 2.989 0.000) 4.227 2.261 ( 3.895 3.895 0.000) 5.508 2.362 ( 8.456 8.456 0.000) 11.958 3.418 ( -6.631 -6.631 0.000) 9.377 4.078 ( 8.923 8.923 0.000) 12.620 4.189 ( 9.119 9.119 0.000) 12.896 6.552 ( -11.538 -11.538 0.000) 16.318 8.289 ( -0.103 -0.103 0.000) 0.146 8.360 ( -4.653 -4.653 0.000) 6.580 9.057 ( 23.789 23.789 0.000) 33.643 9.083 ( -3.793 -3.793 0.000) 5.364 11.657 ( 8.786 8.786 0.000) 12.425 11.805 ( -2.425 -2.425 0.000) 3.429 12.574 ( -4.516 -4.516 0.000) 6.387 13.282 ( -7.478 -7.478 0.000) 10.576 13.480 ( -5.998 -5.998 0.000) 8.483 13.689 ( 12.028 12.028 0.000) 17.010 15.495 ( -5.274 -5.274 0.000) 7.459 15.860 ( 4.142 4.142 0.000) 5.858 ======================= Grid point 15 (9/84) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.805 ( 19.502 8.341 0.000) 21.210 1.021 ( 25.665 -3.536 0.000) 25.907 1.813 ( 16.694 10.838 0.000) 19.903 1.880 ( 2.552 0.989 0.000) 2.737 2.019 ( 9.617 6.685 0.000) 11.712 2.040 ( 5.013 10.658 0.000) 11.778 2.421 ( 12.313 0.119 0.000) 12.314 2.571 ( 11.637 12.922 0.000) 17.390 3.307 ( -3.841 -6.683 0.000) 7.708 4.131 ( 1.637 18.910 0.000) 18.980 4.585 ( 24.783 -0.595 0.000) 24.791 6.267 ( -15.751 -4.161 0.000) 16.291 8.395 ( 5.867 -21.079 0.000) 21.881 8.472 ( 16.701 -17.791 0.000) 24.402 8.952 ( -6.578 -2.075 0.000) 6.898 9.557 ( 24.922 18.461 0.000) 31.015 11.745 ( -3.468 2.006 0.000) 4.006 11.803 ( 6.040 15.448 0.000) 16.586 12.498 ( -2.612 -10.408 0.000) 10.731 13.112 ( -9.883 -12.127 0.000) 15.644 13.365 ( -6.135 -10.315 0.000) 12.002 13.957 ( 14.415 15.138 0.000) 20.904 15.337 ( -9.601 -1.886 0.000) 9.785 15.879 ( -2.276 11.334 0.000) 11.560 ======================= Grid point 16 (10/84) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.217 ( 21.790 5.841 0.000) 22.559 1.467 ( 18.160 -0.548 0.000) 18.168 1.910 ( 0.553 0.799 0.000) 0.972 1.966 ( 4.594 3.573 0.000) 5.820 2.129 ( 4.520 11.079 0.000) 11.966 2.326 ( 15.554 1.736 0.000) 15.651 2.695 ( 13.861 -0.877 0.000) 13.888 2.771 ( 7.341 12.451 0.000) 14.454 3.306 ( 4.760 -7.450 0.000) 8.841 4.180 ( 3.481 20.998 0.000) 21.285 5.098 ( 26.302 -0.150 0.000) 26.302 5.981 ( -12.302 -2.011 0.000) 12.465 8.409 ( -7.988 -10.200 0.000) 12.956 8.889 ( 24.453 -25.958 0.000) 35.662 8.975 ( 11.971 -26.006 0.000) 28.629 10.026 ( 22.288 12.096 0.000) 25.358 11.674 ( -3.578 7.744 0.000) 8.531 11.883 ( 1.757 12.886 0.000) 13.005 12.492 ( 2.114 -12.329 0.000) 12.508 12.885 ( -13.208 -14.790 0.000) 19.829 13.233 ( -7.014 -11.244 0.000) 13.252 14.234 ( 13.086 16.979 0.000) 21.437 15.135 ( -10.594 -0.537 0.000) 10.608 15.783 ( -6.950 16.142 0.000) 17.574 ======================= Grid point 17 (11/84) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.649 ( 21.385 4.631 0.000) 21.881 1.715 ( 7.323 5.797 0.000) 9.340 1.918 ( 0.954 1.070 0.000) 1.434 2.077 ( 6.636 2.018 0.000) 6.936 2.217 ( 4.187 11.503 0.000) 12.242 2.544 ( 6.066 0.292 0.000) 6.073 2.865 ( 3.693 7.376 0.000) 8.249 2.921 ( 8.148 -0.376 0.000) 8.156 3.480 ( 10.780 -11.136 0.000) 15.499 4.276 ( 6.375 23.724 0.000) 24.566 5.592 ( 23.026 0.206 0.000) 23.026 5.798 ( -6.204 -2.647 0.000) 6.745 8.165 ( -14.020 -1.300 0.000) 14.080 9.328 ( 20.787 -42.791 0.000) 47.572 9.397 ( 25.990 -29.333 0.000) 39.190 10.424 ( 17.683 1.633 0.000) 17.758 11.612 ( -2.456 14.220 0.000) 14.430 11.885 ( -0.764 5.267 0.000) 5.322 12.544 ( 1.155 -15.772 0.000) 15.815 12.592 ( -16.265 -15.529 0.000) 22.488 13.114 ( -3.761 -7.479 0.000) 8.372 14.460 ( 9.604 18.380 0.000) 20.738 14.923 ( -10.823 0.842 0.000) 10.855 15.625 ( -8.480 19.729 0.000) 21.474 ======================= Grid point 18 (12/84) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.799 ( 2.196 7.977 0.000) 8.274 1.898 ( 1.932 1.042 0.000) 2.195 2.101 ( 18.806 4.409 0.000) 19.316 2.199 ( 5.036 2.690 0.000) 5.710 2.290 ( 3.213 12.645 0.000) 13.047 2.589 ( -0.187 3.706 0.000) 3.710 2.946 ( 4.023 3.625 0.000) 5.415 2.996 ( -0.598 0.338 0.000) 0.687 3.661 ( 6.774 -13.845 0.000) 15.413 4.422 ( 7.760 23.406 0.000) 24.659 5.726 ( -1.607 -3.719 0.000) 4.051 5.980 ( 15.690 0.361 0.000) 15.694 7.915 ( -10.478 0.187 0.000) 10.480 9.724 ( 17.982 -45.657 0.000) 49.071 9.875 ( 21.292 -29.425 0.000) 36.320 10.714 ( 11.533 -14.977 0.000) 18.903 11.592 ( 0.901 20.546 0.000) 20.566 11.880 ( 0.668 1.271 0.000) 1.436 12.255 ( -17.411 -12.881 0.000) 21.657 12.473 ( -8.047 -27.886 0.000) 29.024 13.123 ( 4.083 -0.222 0.000) 4.090 14.610 ( 5.611 19.465 0.000) 20.257 14.719 ( -9.326 3.075 0.000) 9.819 15.470 ( -6.726 22.343 0.000) 23.334 ======================= Grid point 19 (13/84) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.824 ( 0.524 8.387 0.000) 8.403 1.917 ( 0.416 0.679 0.000) 0.796 2.265 ( 1.676 3.234 0.000) 3.642 2.317 ( 1.535 13.235 0.000) 13.324 2.451 ( 11.813 5.710 0.000) 13.121 2.577 ( -0.542 5.463 0.000) 5.490 2.919 ( -5.488 0.083 0.000) 5.489 3.004 ( 1.600 2.619 0.000) 3.069 3.745 ( 1.944 -13.167 0.000) 13.310 4.543 ( 3.562 21.104 0.000) 21.403 5.714 ( -0.076 -4.260 0.000) 4.261 6.193 ( 5.560 0.406 0.000) 5.575 7.772 ( -3.742 0.287 0.000) 3.753 9.986 ( 7.373 -41.050 0.000) 41.707 10.185 ( 8.748 -27.723 0.000) 29.070 10.871 ( 4.136 -33.646 0.000) 33.899 11.663 ( 6.159 23.498 0.000) 24.292 11.915 ( 2.249 5.465 0.000) 5.910 11.937 ( -13.263 -6.651 0.000) 14.837 12.300 ( -6.861 -41.099 0.000) 41.667 13.209 ( 3.017 1.266 0.000) 3.272 14.580 ( -4.015 5.330 0.000) 6.673 14.683 ( 1.820 20.064 0.000) 20.146 15.376 ( -2.543 23.783 0.000) 23.918 ======================= Grid point 28 (14/84) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.045 ( 15.507 15.507 0.000) 21.931 1.130 ( 14.501 14.501 0.000) 20.507 1.912 ( 1.772 1.772 0.000) 2.507 2.056 ( 10.228 10.228 0.000) 14.464 2.112 ( 6.027 6.027 0.000) 8.523 2.273 ( 10.524 10.524 0.000) 14.884 2.459 ( 5.700 5.700 0.000) 8.061 2.856 ( 12.716 12.716 0.000) 17.984 3.200 ( -3.105 -3.105 0.000) 4.390 4.555 ( 14.574 14.574 0.000) 20.610 4.621 ( 13.048 13.048 0.000) 18.452 6.131 ( -9.546 -9.546 0.000) 13.500 8.081 ( -9.184 -9.184 0.000) 12.989 8.259 ( -1.844 -1.844 0.000) 2.608 8.848 ( -8.416 -8.416 0.000) 11.901 9.881 ( 13.852 13.852 0.000) 19.590 11.806 ( 4.025 4.025 0.000) 5.692 12.101 ( 12.470 12.470 0.000) 17.635 12.248 ( -12.687 -12.687 0.000) 17.942 12.846 ( -14.538 -14.538 0.000) 20.560 13.175 ( -8.687 -8.687 0.000) 12.285 14.287 ( 17.258 17.258 0.000) 24.407 15.241 ( -6.888 -6.888 0.000) 9.742 16.063 ( 5.743 5.743 0.000) 8.121 ======================= Grid point 29 (15/84) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.395 ( 19.497 12.111 0.000) 22.952 1.503 ( 20.459 5.841 0.000) 21.276 1.937 ( 0.631 1.636 0.000) 1.754 2.126 ( 2.072 11.316 0.000) 11.504 2.347 ( 12.422 0.770 0.000) 12.446 2.401 ( 4.643 15.473 0.000) 16.154 2.680 ( 13.809 -0.197 0.000) 13.810 2.988 ( 1.339 5.625 0.000) 5.782 3.230 ( 5.489 2.039 0.000) 5.856 4.660 ( 5.274 26.389 0.000) 26.911 5.113 ( 26.529 2.168 0.000) 26.618 5.915 ( -11.305 -4.948 0.000) 12.341 8.035 ( 3.290 -23.335 0.000) 23.566 8.386 ( 13.465 -21.360 0.000) 25.250 8.635 ( -11.986 -11.205 0.000) 16.408 10.108 ( 9.230 -2.234 0.000) 9.496 11.890 ( 3.963 12.983 0.000) 13.574 11.949 ( -13.913 -13.767 0.000) 19.573 12.349 ( 8.571 5.498 0.000) 10.183 12.543 ( -15.964 -19.229 0.000) 24.992 13.021 ( -6.425 -9.859 0.000) 11.768 14.612 ( 15.170 19.859 0.000) 24.990 15.090 ( -8.074 -3.529 0.000) 8.811 16.101 ( -1.408 13.614 0.000) 13.687 ======================= Grid point 30 (16/84) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.790 ( 18.787 9.250 0.000) 20.941 1.850 ( 13.006 7.638 0.000) 15.083 1.952 ( 2.514 2.520 0.000) 3.559 2.185 ( 4.468 8.980 0.000) 10.031 2.501 ( 5.389 15.979 0.000) 16.863 2.553 ( 8.047 0.810 0.000) 8.088 2.921 ( 9.383 0.780 0.000) 9.415 3.002 ( 1.267 4.415 0.000) 4.593 3.346 ( 5.621 -1.522 0.000) 5.823 4.781 ( 6.681 26.653 0.000) 27.478 5.610 ( 23.050 1.893 0.000) 23.127 5.721 ( -7.702 -5.111 0.000) 9.243 8.082 ( -4.005 -9.922 0.000) 10.699 8.502 ( 3.989 -36.754 0.000) 36.970 8.741 ( 21.344 -33.584 0.000) 39.793 10.236 ( 3.322 -16.793 0.000) 17.119 11.691 ( -11.948 -25.065 0.000) 27.767 11.927 ( -2.121 11.604 0.000) 11.796 12.251 ( -11.246 -14.754 0.000) 18.552 12.407 ( -3.733 3.124 0.000) 4.868 12.963 ( 1.372 -8.146 0.000) 8.261 14.874 ( 11.125 21.721 0.000) 24.404 14.934 ( -7.440 0.711 0.000) 7.474 16.033 ( -4.781 19.027 0.000) 19.619 ======================= Grid point 31 (17/84) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.925 ( 0.863 1.605 0.000) 1.822 2.003 ( 3.721 11.919 0.000) 12.487 2.238 ( 20.022 9.151 0.000) 22.014 2.295 ( 5.550 7.177 0.000) 9.073 2.601 ( 4.457 17.561 0.000) 18.117 2.661 ( 3.135 3.273 0.000) 4.532 3.021 ( 0.324 2.688 0.000) 2.708 3.045 ( 2.700 5.625 0.000) 6.239 3.448 ( 4.525 -7.417 0.000) 8.689 4.910 ( 5.911 25.440 0.000) 26.118 5.618 ( -2.973 -7.144 0.000) 7.738 5.998 ( 15.625 1.550 0.000) 15.702 7.900 ( -10.351 -2.495 0.000) 10.648 8.771 ( 18.976 -49.006 0.000) 52.551 9.180 ( 21.692 -38.137 0.000) 43.874 10.217 ( -5.618 -30.877 0.000) 31.384 11.479 ( -9.550 -36.831 0.000) 38.049 11.764 ( -11.154 -14.245 0.000) 18.092 12.169 ( 0.601 13.345 0.000) 13.358 12.244 ( -10.571 0.649 0.000) 10.591 13.043 ( 5.115 -7.465 0.000) 9.049 14.805 ( -5.394 5.728 0.000) 7.868 15.048 ( 6.555 22.888 0.000) 23.808 15.936 ( -4.522 22.306 0.000) 22.760 ======================= Grid point 32 (18/84) ======================= q-point: ( 0.46 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.936 ( 0.276 1.238 0.000) 1.269 2.042 ( 0.763 12.950 0.000) 12.973 2.371 ( 1.961 7.494 0.000) 7.746 2.597 ( 14.514 10.135 0.000) 17.702 2.665 ( 1.741 18.722 0.000) 18.802 2.694 ( 0.609 5.643 0.000) 5.676 2.940 ( -7.286 3.100 0.000) 7.918 3.089 ( 1.313 5.969 0.000) 6.112 3.513 ( 1.801 -9.970 0.000) 10.131 4.996 ( 2.443 24.045 0.000) 24.169 5.588 ( -0.497 -8.436 0.000) 8.451 6.209 ( 5.519 1.316 0.000) 5.674 7.749 ( -4.149 -4.072 0.000) 5.813 9.129 ( 13.531 -43.162 0.000) 45.233 9.522 ( 10.481 -37.332 0.000) 38.775 10.047 ( -8.606 -46.259 0.000) 47.053 11.328 ( -4.826 -47.608 0.000) 47.852 11.600 ( -4.429 -19.472 0.000) 19.970 12.071 ( -5.238 2.097 0.000) 5.642 12.191 ( 0.827 20.936 0.000) 20.952 13.121 ( 2.264 -9.153 0.000) 9.428 14.730 ( -2.008 9.553 0.000) 9.762 15.134 ( 2.138 23.440 0.000) 23.537 15.871 ( -1.794 23.856 0.000) 23.923 ======================= Grid point 42 (19/84) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.683 ( 16.461 16.461 0.000) 23.280 1.717 ( 14.981 14.981 0.000) 21.186 1.965 ( 1.224 1.224 0.000) 1.731 2.348 ( 5.975 5.975 0.000) 8.450 2.383 ( 7.390 7.390 0.000) 10.451 2.700 ( 6.593 6.593 0.000) 9.324 2.715 ( 11.689 11.689 0.000) 16.531 3.006 ( -0.467 -0.467 0.000) 0.660 3.336 ( 5.162 5.162 0.000) 7.300 5.172 ( 15.916 15.916 0.000) 22.509 5.204 ( 15.571 15.571 0.000) 22.021 5.780 ( -8.604 -8.604 0.000) 12.168 7.674 ( -10.881 -10.881 0.000) 15.388 8.118 ( -5.285 -5.285 0.000) 7.474 8.393 ( -14.680 -14.680 0.000) 20.760 10.014 ( -5.950 -5.950 0.000) 8.415 11.550 ( -22.321 -22.321 0.000) 31.566 12.134 ( -21.521 -21.521 0.000) 30.436 12.142 ( 11.362 11.362 0.000) 16.068 12.512 ( 6.418 6.418 0.000) 9.077 12.839 ( -8.184 -8.184 0.000) 11.574 14.983 ( 17.115 17.115 0.000) 24.204 15.017 ( -3.451 -3.451 0.000) 4.881 16.287 ( 5.336 5.336 0.000) 7.547 ======================= Grid point 43 (20/84) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.940 ( 3.048 3.441 0.000) 4.597 2.032 ( 15.071 11.180 0.000) 18.766 2.077 ( 16.518 12.252 0.000) 20.566 2.396 ( 2.554 11.179 0.000) 11.467 2.584 ( 9.150 2.579 0.000) 9.507 2.798 ( 6.084 12.658 0.000) 14.044 2.983 ( 11.672 6.815 0.000) 13.516 3.027 ( 3.001 -0.196 0.000) 3.008 3.392 ( 0.810 3.953 0.000) 4.035 5.302 ( 5.887 25.464 0.000) 26.136 5.597 ( -9.131 -7.554 0.000) 11.851 5.672 ( 22.218 4.326 0.000) 22.635 7.634 ( 6.411 -27.998 0.000) 28.723 8.060 ( -17.962 -14.232 0.000) 22.917 8.169 ( 9.762 -22.685 0.000) 24.696 9.879 ( -8.590 -17.501 0.000) 19.495 11.120 ( -20.873 -29.882 0.000) 36.450 11.689 ( -22.263 -24.963 0.000) 33.449 12.344 ( 8.514 13.051 0.000) 15.583 12.483 ( -9.699 2.624 0.000) 10.048 12.775 ( 2.040 -11.344 0.000) 11.526 14.963 ( -1.848 2.416 0.000) 3.042 15.275 ( 12.314 18.312 0.000) 22.068 16.331 ( -0.226 10.912 0.000) 10.914 ======================= Grid point 44 (21/84) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.958 ( 0.366 1.620 0.000) 1.661 2.247 ( 6.103 12.431 0.000) 13.848 2.449 ( 18.999 11.887 0.000) 22.411 2.472 ( 4.876 10.020 0.000) 11.143 2.726 ( 5.200 3.720 0.000) 6.394 2.917 ( 5.192 12.874 0.000) 13.882 3.119 ( 5.746 1.281 0.000) 5.887 3.146 ( 4.337 11.291 0.000) 12.095 3.386 ( -1.036 1.099 0.000) 1.510 5.406 ( 4.420 24.227 0.000) 24.627 5.449 ( -5.406 -9.812 0.000) 11.203 6.044 ( 14.947 3.149 0.000) 15.275 7.690 ( -13.095 -29.453 0.000) 32.233 7.944 ( 18.905 -23.371 0.000) 30.060 8.458 ( 17.835 -33.173 0.000) 37.664 9.625 ( -18.240 -26.598 0.000) 32.251 10.733 ( -18.022 -36.159 0.000) 40.401 11.303 ( -15.493 -28.935 0.000) 32.822 12.242 ( -12.590 -1.908 0.000) 12.734 12.478 ( 5.076 14.378 0.000) 15.248 12.830 ( 2.567 -14.129 0.000) 14.361 14.942 ( -0.446 7.887 0.000) 7.900 15.467 ( 7.157 18.895 0.000) 20.205 16.303 ( -1.980 14.553 0.000) 14.687 ======================= Grid point 45 (22/84) ======================= q-point: ( 0.46 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.964 ( 0.172 1.538 0.000) 1.547 2.308 ( 1.104 13.410 0.000) 13.455 2.550 ( 2.235 9.964 0.000) 10.212 2.784 ( 13.656 8.091 0.000) 15.873 2.793 ( 1.636 4.742 0.000) 5.016 2.966 ( -1.546 11.574 0.000) 11.677 3.160 ( -1.543 19.446 0.000) 19.507 3.220 ( 3.421 5.963 0.000) 6.874 3.362 ( -1.099 -4.522 0.000) 4.654 5.384 ( -1.422 -12.081 0.000) 12.164 5.468 ( 1.709 23.138 0.000) 23.201 6.247 ( 5.249 2.400 0.000) 5.772 7.485 ( -5.773 -26.353 0.000) 26.978 8.420 ( 24.822 -25.397 0.000) 35.513 8.775 ( 10.892 -36.925 0.000) 38.498 9.174 ( -24.117 -38.569 0.000) 45.489 10.444 ( -9.351 -40.024 0.000) 41.102 11.104 ( -4.780 -30.426 0.000) 30.799 12.049 ( -5.633 -4.220 0.000) 7.039 12.545 ( 1.697 14.804 0.000) 14.901 12.865 ( 0.894 -16.600 0.000) 16.624 14.939 ( 0.003 11.099 0.000) 11.099 15.560 ( 2.316 19.108 0.000) 19.248 16.271 ( -0.985 16.295 0.000) 16.325 ======================= Grid point 56 (23/84) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.974 ( 0.882 0.882 0.000) 1.247 2.278 ( 12.759 12.759 0.000) 18.044 2.373 ( 16.466 16.466 0.000) 23.287 2.568 ( 5.128 5.128 0.000) 7.252 2.681 ( 7.694 7.694 0.000) 10.881 2.988 ( 8.202 8.202 0.000) 11.600 3.045 ( 2.509 2.509 0.000) 3.549 3.195 ( 12.594 12.594 0.000) 17.811 3.446 ( 1.089 1.089 0.000) 1.540 5.424 ( -9.670 -9.670 0.000) 13.676 5.757 ( 13.043 13.043 0.000) 18.445 5.775 ( 12.659 12.659 0.000) 17.903 7.207 ( -14.800 -14.800 0.000) 20.930 7.771 ( -13.543 -13.543 0.000) 19.153 7.886 ( -5.588 -5.588 0.000) 7.902 9.568 ( -14.845 -14.845 0.000) 20.994 10.561 ( -26.177 -26.177 0.000) 37.020 11.153 ( -27.852 -27.852 0.000) 39.389 12.500 ( -7.031 -7.031 0.000) 9.943 12.507 ( -9.098 -9.098 0.000) 12.866 12.568 ( 9.296 9.296 0.000) 13.146 15.027 ( 3.806 3.806 0.000) 5.383 15.584 ( 12.868 12.868 0.000) 18.198 16.474 ( 4.093 4.093 0.000) 5.789 ======================= Grid point 57 (24/84) ======================= q-point: ( 0.38 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.987 ( 0.494 1.270 0.000) 1.363 2.483 ( 6.821 11.150 0.000) 13.071 2.648 ( 4.221 6.766 0.000) 7.975 2.695 ( 15.274 12.741 0.000) 19.890 2.830 ( 6.237 7.320 0.000) 9.617 3.130 ( 5.751 0.842 0.000) 5.812 3.142 ( 6.925 11.642 0.000) 13.546 3.424 ( 9.762 14.746 0.000) 17.684 3.439 ( -1.524 2.922 0.000) 3.295 5.237 ( -8.450 -11.443 0.000) 14.225 5.845 ( 3.157 19.749 0.000) 20.000 6.113 ( 13.760 3.453 0.000) 14.187 6.972 ( -11.321 -34.463 0.000) 36.275 7.669 ( 6.567 -10.868 0.000) 12.698 7.927 ( 8.167 -19.721 0.000) 21.345 9.220 ( -22.441 -14.764 0.000) 26.862 10.101 ( -19.994 -27.091 0.000) 33.670 10.661 ( -19.875 -35.180 0.000) 40.406 12.174 ( -16.253 -4.534 0.000) 16.874 12.490 ( -0.826 -19.821 0.000) 19.838 12.711 ( 5.274 9.301 0.000) 10.692 15.101 ( 3.337 7.776 0.000) 8.462 15.783 ( 7.343 13.027 0.000) 14.954 16.515 ( 0.618 7.263 0.000) 7.289 ======================= Grid point 58 (25/84) ======================= q-point: ( 0.46 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.994 ( 0.204 1.478 0.000) 1.492 2.553 ( 1.230 11.045 0.000) 11.114 2.723 ( 2.322 6.336 0.000) 6.748 2.913 ( 2.111 8.222 0.000) 8.489 2.928 ( 6.852 7.500 0.000) 10.159 3.241 ( 4.313 -1.640 0.000) 4.614 3.243 ( 2.914 16.361 0.000) 16.619 3.397 ( -2.167 5.474 0.000) 5.887 3.551 ( 2.940 18.685 0.000) 18.915 5.119 ( -3.141 -14.552 0.000) 14.887 5.888 ( 1.136 18.857 0.000) 18.891 6.298 ( 4.773 2.579 0.000) 5.425 6.783 ( -5.827 -39.484 0.000) 39.911 8.021 ( 26.228 -18.369 0.000) 32.020 8.110 ( 7.069 -29.580 0.000) 30.413 8.669 ( -31.203 -11.816 0.000) 33.365 9.802 ( -8.881 -24.018 0.000) 25.607 10.411 ( -5.820 -39.282 0.000) 39.711 11.926 ( -7.267 -7.639 0.000) 10.544 12.478 ( -0.365 -22.132 0.000) 22.135 12.780 ( 1.707 9.174 0.000) 9.332 15.148 ( 1.293 9.681 0.000) 9.767 15.878 ( 2.352 13.041 0.000) 13.251 16.516 ( -0.132 8.893 0.000) 8.894 ======================= Grid point 70 (26/84) ======================= q-point: ( 0.38 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.008 ( 0.852 0.852 0.000) 1.205 2.665 ( 6.565 6.565 0.000) 9.284 2.740 ( 3.408 3.408 0.000) 4.819 2.934 ( 10.750 10.750 0.000) 15.203 2.990 ( 7.538 7.538 0.000) 10.660 3.168 ( 2.942 2.942 0.000) 4.160 3.388 ( 12.571 12.571 0.000) 17.779 3.481 ( 1.272 1.272 0.000) 1.799 3.692 ( 12.028 12.028 0.000) 17.010 5.007 ( -11.475 -11.475 0.000) 16.228 6.165 ( 7.497 7.497 0.000) 10.602 6.174 ( 7.443 7.443 0.000) 10.525 6.393 ( -23.760 -23.760 0.000) 33.601 7.567 ( 2.810 2.810 0.000) 3.973 7.681 ( -6.077 -6.077 0.000) 8.594 8.963 ( -14.357 -14.357 0.000) 20.304 9.684 ( -14.872 -14.872 0.000) 21.032 9.983 ( -30.466 -30.466 0.000) 43.086 12.070 ( -12.606 -12.606 0.000) 17.828 12.098 ( -11.474 -11.474 0.000) 16.227 12.852 ( 5.051 5.051 0.000) 7.144 15.235 ( 5.474 5.474 0.000) 7.741 15.983 ( 7.336 7.336 0.000) 10.375 16.608 ( 2.699 2.699 0.000) 3.817 ======================= Grid point 71 (27/84) ======================= q-point: ( 0.46 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 255 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.020 ( 0.315 1.064 0.000) 1.110 2.733 ( 1.127 6.991 0.000) 7.081 2.801 ( 2.067 2.293 0.000) 3.087 3.088 ( 4.150 8.029 0.000) 9.038 3.092 ( 2.534 8.513 0.000) 8.882 3.226 ( 2.070 0.149 0.000) 2.076 3.494 ( 0.211 3.932 0.000) 3.938 3.592 ( 6.440 18.434 0.000) 19.527 3.873 ( 5.293 13.523 0.000) 14.522 4.819 ( -6.022 -15.485 0.000) 16.615 6.061 ( -8.983 -31.541 0.000) 32.795 6.198 ( 0.726 12.240 0.000) 12.262 6.341 ( 4.420 1.688 0.000) 4.731 7.631 ( -0.433 -18.575 0.000) 18.580 7.784 ( 15.648 -2.670 0.000) 15.874 8.575 ( -21.471 -1.017 0.000) 21.495 9.505 ( -3.945 -6.702 0.000) 7.777 9.588 ( -8.823 -42.620 0.000) 43.524 11.775 ( -9.068 -6.927 0.000) 11.411 12.026 ( -1.666 -22.225 0.000) 22.288 12.916 ( 1.590 4.840 0.000) 5.095 15.309 ( 1.946 6.373 0.000) 6.663 16.077 ( 2.334 7.300 0.000) 7.664 16.638 ( 0.585 3.837 0.000) 3.881 ======================= Grid point 84 (28/84) ======================= q-point: ( 0.46 0.46 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 147 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.034 ( 0.387 0.387 0.000) 0.547 2.824 ( 1.841 1.841 0.000) 2.604 2.830 ( 1.154 1.154 0.000) 1.632 3.208 ( 3.437 3.437 0.000) 4.861 3.211 ( 3.042 3.042 0.000) 4.302 3.236 ( 0.617 0.617 0.000) 0.873 3.547 ( 1.421 1.421 0.000) 2.009 3.933 ( 13.586 13.586 0.000) 19.214 4.068 ( 5.768 5.768 0.000) 8.157 4.528 ( -12.326 -12.326 0.000) 17.431 5.611 ( -12.359 -12.359 0.000) 17.478 6.362 ( 2.595 2.595 0.000) 3.669 6.365 ( 2.321 2.321 0.000) 3.283 7.383 ( -6.409 -6.409 0.000) 9.063 7.936 ( 15.173 15.173 0.000) 21.458 8.539 ( -6.026 -6.026 0.000) 8.522 8.843 ( -25.575 -25.575 0.000) 36.168 9.472 ( 0.922 0.922 0.000) 1.304 11.669 ( -5.993 -5.993 0.000) 8.476 11.685 ( -7.283 -7.283 0.000) 10.299 12.978 ( 1.494 1.494 0.000) 2.113 15.395 ( 2.201 2.201 0.000) 3.113 16.171 ( 2.317 2.317 0.000) 3.276 16.683 ( 0.989 0.989 0.000) 1.399 ======================= Grid point 169 (29/84) ======================= q-point: ( 0.00 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 84 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.403 ( 0.000 0.000 18.256) 18.256 0.403 ( -0.000 0.000 18.256) 18.256 0.664 ( -0.000 -0.000 29.805) 29.805 1.601 ( 0.000 0.000 -10.669) 10.669 1.601 ( 0.000 0.000 -10.669) 10.669 1.964 ( 0.000 0.000 4.397) 4.397 1.964 ( 0.000 0.000 4.397) 4.397 2.658 ( 0.000 0.000 -13.674) 13.674 3.675 ( 0.000 0.000 4.389) 4.389 3.901 ( 0.000 0.000 0.537) 0.537 3.901 ( 0.000 0.000 0.537) 0.537 6.982 ( 0.000 0.000 0.049) 0.049 8.300 ( 0.000 0.000 0.988) 0.988 8.300 ( 0.000 0.000 0.988) 0.988 8.433 ( 0.000 0.000 -1.464) 1.464 8.433 ( 0.000 0.000 -1.464) 1.464 9.167 ( -0.000 0.000 0.745) 0.745 12.652 ( -0.000 0.000 -0.328) 0.328 13.423 ( 0.000 0.000 -0.588) 0.588 13.423 ( 0.000 0.000 -0.588) 0.588 13.613 ( 0.000 0.000 0.758) 0.758 13.613 ( -0.000 -0.000 0.758) 0.758 15.424 ( 0.000 -0.000 2.377) 2.377 15.766 ( 0.000 0.000 -0.235) 0.235 ======================= Grid point 170 (30/84) ======================= q-point: ( 0.08 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.532 ( 7.548 0.000 15.504) 17.244 0.724 ( 23.908 0.000 9.992) 25.912 1.050 ( 33.757 0.000 12.003) 35.827 1.642 ( 3.917 0.000 -10.192) 10.918 1.747 ( 9.923 0.000 -8.396) 12.999 1.974 ( 1.079 0.000 4.279) 4.413 2.041 ( 6.873 0.000 -1.191) 6.975 2.452 ( -3.490 0.000 -0.190) 3.496 3.544 ( -8.089 0.000 2.807) 8.562 3.906 ( 0.467 0.000 0.476) 0.667 4.131 ( 18.730 0.000 -1.928) 18.829 6.723 ( -18.439 0.000 3.208) 18.716 8.414 ( 11.386 0.000 1.002) 11.430 8.558 ( 12.512 0.000 -1.424) 12.592 8.640 ( 16.740 0.000 0.205) 16.741 9.140 ( -1.418 0.000 0.763) 1.610 9.937 ( 81.726 0.000 -50.280) 95.954 12.612 ( -4.008 0.000 0.016) 4.008 13.070 ( -33.642 0.000 28.034) 43.792 13.383 ( -4.058 0.000 -0.453) 4.083 13.586 ( -2.774 0.000 0.622) 2.843 13.628 ( 9.890 0.000 10.335) 14.305 15.251 ( -8.141 0.000 -16.003) 17.955 15.729 ( -2.911 0.000 -1.553) 3.299 ======================= Grid point 171 (31/84) ======================= q-point: ( 0.15 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.771 ( 16.653 0.000 4.369) 17.217 1.179 ( 20.579 0.000 5.556) 21.316 1.728 ( 29.700 0.000 0.241) 29.701 1.743 ( 6.005 0.000 -8.564) 10.460 1.886 ( 6.520 0.000 1.400) 6.669 2.007 ( 2.240 0.000 3.707) 4.331 2.252 ( 13.874 0.000 -7.133) 15.600 2.555 ( 10.869 0.000 2.637) 11.184 3.416 ( -4.318 0.000 1.332) 4.519 3.918 ( 0.709 0.000 0.328) 0.781 4.569 ( 24.784 0.000 -2.093) 24.873 6.345 ( -18.420 0.000 3.010) 18.665 8.736 ( 20.731 0.000 0.972) 20.754 8.788 ( -6.021 0.000 4.194) 7.337 8.912 ( 22.831 0.000 -1.253) 22.865 9.294 ( 19.430 0.000 0.104) 19.430 10.993 ( 33.376 0.000 -24.920) 41.653 12.358 ( -24.926 0.000 23.698) 34.393 12.616 ( -6.018 0.000 2.549) 6.535 13.262 ( -8.159 0.000 -0.094) 8.159 13.505 ( -5.394 0.000 0.258) 5.400 13.872 ( 13.251 0.000 9.828) 16.498 15.152 ( -3.635 0.000 -17.297) 17.675 15.645 ( -5.802 0.000 -5.282) 7.846 ======================= Grid point 172 (32/84) ======================= q-point: ( 0.23 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.155 ( 21.416 0.000 -0.504) 21.422 1.502 ( 11.757 0.000 2.591) 12.039 1.877 ( 7.700 0.000 -5.211) 9.298 1.886 ( 0.982 0.000 -1.756) 2.011 2.060 ( 2.995 0.000 2.496) 3.899 2.314 ( 21.393 0.000 1.886) 21.476 2.508 ( 10.501 0.000 -6.062) 12.125 2.809 ( 13.564 0.000 4.687) 14.351 3.414 ( 5.539 0.000 -0.299) 5.547 3.932 ( 0.657 0.000 0.170) 0.679 5.081 ( 26.349 0.000 -1.663) 26.401 6.029 ( -13.134 0.000 2.221) 13.320 8.544 ( -16.225 0.000 4.545) 16.850 9.207 ( 26.142 0.000 0.836) 26.155 9.435 ( 29.444 0.000 -0.927) 29.459 9.730 ( 22.082 0.000 0.828) 22.098 11.446 ( 14.111 0.000 -6.941) 15.726 11.921 ( -17.196 0.000 11.357) 20.608 12.587 ( 2.003 0.000 -3.243) 3.812 13.061 ( -12.278 0.000 0.367) 12.283 13.375 ( -7.672 0.000 -0.232) 7.675 14.122 ( 11.686 0.000 7.734) 14.013 15.061 ( -6.181 0.000 -6.463) 8.942 15.501 ( -8.566 0.000 -7.474) 11.368 ======================= Grid point 173 (33/84) ======================= q-point: ( 0.31 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.593 ( 22.383 0.000 -1.122) 22.411 1.655 ( 4.638 0.000 0.355) 4.652 1.901 ( 0.922 0.000 -0.497) 1.048 2.035 ( 7.618 0.000 -1.883) 7.847 2.119 ( 2.807 0.000 1.257) 3.075 2.569 ( 4.145 0.000 2.260) 4.721 2.676 ( 8.082 0.000 -5.973) 10.050 3.035 ( 8.276 0.000 6.325) 10.416 3.647 ( 16.263 0.000 -3.393) 16.614 3.943 ( 0.431 0.000 0.059) 0.435 5.579 ( 23.261 0.000 -1.039) 23.284 5.837 ( -6.263 0.000 1.005) 6.343 8.198 ( -17.851 0.000 1.910) 17.953 9.734 ( 26.312 0.000 0.638) 26.320 10.048 ( 31.940 0.000 -0.520) 31.945 10.124 ( 17.341 0.000 0.604) 17.352 11.508 ( -7.303 0.000 3.427) 8.067 11.767 ( 1.601 0.000 -1.439) 2.153 12.693 ( 8.431 0.000 -3.349) 9.072 12.778 ( -16.292 0.000 0.813) 16.313 13.207 ( -9.109 0.000 -0.739) 9.139 14.326 ( 8.815 0.000 6.539) 10.976 14.900 ( -10.005 0.000 -0.484) 10.017 15.323 ( -8.954 0.000 -6.696) 11.181 ======================= Grid point 174 (34/84) ======================= q-point: ( 0.38 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.714 ( 1.756 0.000 -0.163) 1.763 1.895 ( 2.130 0.000 0.613) 2.217 2.053 ( 20.014 0.000 -1.055) 20.041 2.151 ( 3.186 0.000 -0.003) 3.186 2.175 ( 3.500 0.000 -0.091) 3.502 2.580 ( -1.055 0.000 2.863) 3.051 2.795 ( 3.230 0.000 -5.423) 6.312 3.123 ( 0.860 0.000 6.913) 6.967 3.949 ( 0.202 0.000 0.009) 0.202 3.963 ( 13.933 0.000 -6.459) 15.357 5.771 ( -0.994 0.000 0.591) 1.156 5.972 ( 15.943 0.000 -0.447) 15.949 7.877 ( -13.716 0.000 -2.771) 13.993 10.206 ( 20.458 0.000 0.472) 20.464 10.404 ( 10.825 0.000 0.528) 10.838 10.670 ( 30.364 0.000 -0.176) 30.365 11.353 ( -6.608 0.000 -0.904) 6.670 11.814 ( 1.845 0.000 0.610) 1.943 12.421 ( -19.560 0.000 1.184) 19.596 12.889 ( 10.263 0.000 -1.535) 10.377 13.029 ( -8.429 0.000 -1.177) 8.511 14.466 ( 5.292 0.000 6.535) 8.409 14.686 ( -10.784 0.000 0.119) 10.785 15.164 ( -6.761 0.000 -4.272) 7.998 ======================= Grid point 175 (35/84) ======================= q-point: ( 0.46 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.735 ( 0.472 0.000 -0.164) 0.499 1.913 ( 0.388 0.000 0.246) 0.459 2.182 ( 0.650 0.000 0.302) 0.717 2.232 ( 1.578 0.000 -0.369) 1.621 2.417 ( 14.230 0.000 0.730) 14.248 2.556 ( -1.331 0.000 4.504) 4.697 2.786 ( -3.719 0.000 -6.229) 7.255 3.107 ( -1.119 0.000 6.066) 6.169 3.952 ( 0.053 0.000 -0.005) 0.053 4.156 ( 5.128 0.000 -8.438) 9.874 5.771 ( 0.367 0.000 1.042) 1.105 6.189 ( 5.667 0.000 -0.087) 5.668 7.686 ( -5.134 0.000 -6.506) 8.288 10.498 ( 8.034 0.000 0.402) 8.044 10.547 ( 3.677 0.000 0.962) 3.801 11.197 ( 19.897 0.000 -0.208) 19.899 11.264 ( -2.293 0.000 -3.401) 4.102 11.836 ( 0.552 0.000 1.798) 1.881 12.042 ( -16.056 0.000 1.623) 16.138 12.901 ( -3.734 0.000 -1.456) 4.007 13.051 ( 4.957 0.000 0.223) 4.962 14.511 ( -4.020 0.000 -0.364) 4.037 14.541 ( 0.604 0.000 6.275) 6.304 15.069 ( -2.547 0.000 -2.650) 3.676 ======================= Grid point 183 (36/84) ======================= q-point: ( 0.08 0.08 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.622 ( 6.467 6.467 11.565) 14.745 0.698 ( 11.618 11.618 10.783) 19.653 1.417 ( 25.203 25.203 6.612) 36.251 1.729 ( 5.879 5.879 -10.458) 13.361 1.787 ( 4.836 4.836 -4.017) 7.931 2.043 ( 4.220 4.220 4.057) 7.216 2.099 ( 7.063 7.063 -5.295) 11.305 2.453 ( 2.949 2.949 0.199) 4.176 3.446 ( -7.181 -7.181 2.225) 10.396 4.085 ( 8.944 8.944 0.590) 12.663 4.160 ( 9.475 9.475 -1.952) 13.541 6.598 ( -11.857 -11.857 3.246) 17.080 8.299 ( -0.512 -0.512 0.845) 1.112 8.360 ( -3.581 -3.581 -0.119) 5.066 9.064 ( 24.824 24.824 0.564) 35.111 9.095 ( -3.699 -3.699 1.035) 5.332 10.604 ( 43.423 43.423 -38.579) 72.522 12.562 ( -4.875 -4.875 -0.983) 6.965 12.812 ( -22.609 -22.609 34.504) 47.041 13.280 ( -7.236 -7.236 -0.157) 10.234 13.487 ( -6.071 -6.071 0.582) 8.605 13.786 ( 12.213 12.213 10.332) 20.126 15.279 ( 0.543 0.543 -20.480) 20.494 15.725 ( -0.011 -0.011 -5.013) 5.013 ======================= Grid point 184 (37/84) ======================= q-point: ( 0.15 0.08 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.837 ( 15.866 7.042 2.853) 17.591 1.096 ( 23.126 -2.886 6.488) 24.191 1.794 ( 6.978 9.286 -5.059) 12.669 1.864 ( 4.496 3.524 -1.445) 5.892 1.942 ( 17.465 7.215 0.625) 18.907 2.126 ( 4.975 8.746 3.268) 10.579 2.316 ( 14.263 6.455 -6.848) 17.088 2.606 ( 11.288 6.394 1.084) 13.018 3.321 ( -4.695 -6.996 1.139) 8.502 4.133 ( 1.660 18.779 0.164) 18.853 4.566 ( 24.909 0.089 -1.450) 24.952 6.306 ( -16.199 -4.253 2.987) 17.012 8.406 ( 8.217 -21.620 0.758) 23.141 8.489 ( 15.928 -17.361 1.449) 23.605 8.966 ( -7.096 -2.062 1.251) 7.494 9.576 ( 25.177 19.907 1.587) 32.135 11.283 ( 25.272 22.817 -17.955) 38.492 12.313 ( -20.544 -3.900 19.388) 28.516 12.498 ( -5.398 -11.105 2.086) 12.522 13.118 ( -9.533 -11.981 0.424) 15.317 13.369 ( -6.303 -10.111 0.358) 11.920 14.050 ( 13.732 14.578 9.103) 21.999 15.237 ( -4.553 3.450 -11.261) 12.627 15.718 ( -1.303 8.000 -10.708) 13.430 ======================= Grid point 185 (38/84) ======================= q-point: ( 0.23 0.08 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.209 ( 20.928 5.687 -0.702) 21.699 1.504 ( 16.622 0.072 3.405) 16.967 1.890 ( 2.478 3.301 -4.079) 5.803 1.916 ( 4.282 2.906 -1.644) 5.430 2.156 ( 4.555 9.404 2.120) 10.662 2.356 ( 15.092 2.089 2.202) 15.395 2.594 ( 11.953 6.951 -5.899) 15.033 2.840 ( 11.232 4.129 3.033) 12.345 3.296 ( 3.235 -8.034 -0.683) 8.688 4.181 ( 3.311 20.944 0.106) 21.205 5.082 ( 26.528 0.339 -1.232) 26.558 6.008 ( -13.073 -2.137 2.127) 13.416 8.461 ( -7.762 -10.854 4.223) 13.996 8.895 ( 23.648 -26.173 0.627) 35.280 8.971 ( 12.446 -25.453 -0.426) 28.336 10.040 ( 21.490 12.820 1.118) 25.048 11.626 ( 9.335 14.648 -2.768) 17.589 11.979 ( -11.887 5.386 6.987) 14.803 12.451 ( 0.664 -12.536 -3.577) 13.053 12.896 ( -12.943 -14.478 0.850) 19.439 13.234 ( -7.211 -10.905 0.071) 13.073 14.310 ( 12.295 16.301 6.944) 21.566 15.110 ( -8.081 1.522 -2.650) 8.640 15.659 ( -4.669 15.284 -10.241) 18.981 ======================= Grid point 186 (39/84) ======================= q-point: ( 0.31 0.08 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.637 ( 21.643 4.821 -1.042) 22.198 1.727 ( 6.683 6.332 1.053) 9.266 1.912 ( 1.038 1.260 -0.572) 1.730 2.049 ( 7.803 1.900 -2.181) 8.322 2.232 ( 3.422 11.011 1.253) 11.598 2.576 ( 6.473 0.476 2.522) 6.963 2.756 ( 5.428 7.214 -5.508) 10.576 3.028 ( 7.364 -0.569 5.329) 9.108 3.445 ( 9.959 -11.386 -3.021) 15.425 4.269 ( 5.658 23.732 -0.531) 24.403 5.582 ( 23.332 0.495 -0.790) 23.350 5.810 ( -6.805 -2.686 0.971) 7.380 8.189 ( -16.381 -1.194 1.975) 16.543 9.327 ( 21.513 -43.003 -0.035) 48.084 9.407 ( 25.718 -29.315 0.788) 39.005 10.417 ( 16.560 1.382 -0.589) 16.628 11.654 ( -3.641 14.244 3.757) 15.175 11.872 ( -0.667 5.489 -1.220) 5.662 12.498 ( 1.659 -15.863 -3.572) 16.344 12.611 ( -15.451 -15.122 1.125) 21.649 13.109 ( -4.221 -7.299 -0.434) 8.443 14.526 ( 9.383 17.697 5.853) 20.868 14.925 ( -10.304 1.188 0.226) 10.375 15.545 ( -6.445 19.849 -6.989) 22.008 ======================= Grid point 187 (40/84) ======================= q-point: ( 0.38 0.08 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.798 ( 2.067 8.062 -0.099) 8.324 1.905 ( 1.580 0.981 0.542) 1.937 2.096 ( 20.286 4.498 -0.695) 20.790 2.184 ( 5.370 2.428 -1.184) 6.011 2.289 ( 2.198 11.945 0.168) 12.147 2.629 ( 0.285 4.119 3.381) 5.336 2.841 ( 2.690 4.806 -4.951) 7.406 3.110 ( 0.890 -1.143 5.856) 6.033 3.617 ( 6.607 -14.217 -3.943) 16.166 4.395 ( 6.536 24.217 -2.143) 25.175 5.733 ( -1.453 -3.756 0.561) 4.066 5.976 ( 15.957 0.516 -0.313) 15.969 7.885 ( -13.273 0.585 -2.435) 13.507 9.731 ( 18.111 -45.787 0.574) 49.242 9.883 ( 21.225 -29.326 0.620) 36.206 10.690 ( 10.967 -15.554 -2.031) 19.139 11.605 ( -0.111 21.493 1.154) 21.524 11.881 ( 1.441 0.543 0.110) 1.544 12.269 ( -17.282 -12.932 1.219) 21.619 12.460 ( -7.139 -27.084 -1.178) 28.034 13.112 ( 4.217 -0.037 -0.895) 4.311 14.668 ( 3.175 15.765 4.544) 16.712 14.729 ( -7.185 6.034 1.368) 9.482 15.423 ( -5.448 22.425 -4.192) 23.455 ======================= Grid point 188 (41/84) ======================= q-point: ( 0.46 0.08 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.821 ( 0.524 8.372 -0.176) 8.390 1.920 ( 0.336 0.686 0.244) 0.802 2.253 ( 1.705 2.876 -0.965) 3.480 2.317 ( 0.741 12.410 -0.244) 12.434 2.464 ( 14.847 4.744 0.602) 15.598 2.618 ( -1.696 5.681 4.219) 7.277 2.830 ( -3.630 4.439 -4.525) 7.305 3.095 ( -0.995 -0.954 4.908) 5.098 3.700 ( 1.935 -14.178 -4.102) 14.886 4.496 ( 2.934 22.856 -3.724) 23.342 5.727 ( 0.214 -4.386 0.958) 4.494 6.193 ( 5.666 0.496 -0.009) 5.688 7.699 ( -5.022 1.045 -5.896) 7.816 9.992 ( 7.248 -41.246 0.497) 41.881 10.192 ( 8.729 -27.655 0.554) 29.005 10.842 ( 4.100 -33.100 -2.466) 33.444 11.668 ( 6.181 24.119 0.437) 24.902 11.928 ( 1.860 3.440 1.069) 4.054 11.951 ( -12.611 -5.993 1.168) 14.012 12.296 ( -6.832 -40.253 -0.474) 40.832 13.205 ( 3.322 1.205 -0.347) 3.551 14.567 ( -4.266 5.388 -1.007) 6.946 14.757 ( 1.820 19.402 6.356) 20.498 15.346 ( -2.119 23.663 -2.795) 23.921 ======================= Grid point 197 (42/84) ======================= q-point: ( 0.15 0.15 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.051 ( 14.292 14.292 0.417) 20.216 1.202 ( 13.582 13.582 6.345) 20.229 1.909 ( 1.892 1.892 0.135) 2.680 2.005 ( 7.525 7.525 -8.357) 13.530 2.088 ( 9.741 9.741 2.453) 13.992 2.310 ( 9.161 9.161 2.046) 13.117 2.455 ( 8.474 8.474 -0.313) 11.988 2.813 ( 11.941 11.941 -3.613) 17.269 3.200 ( -4.405 -4.405 0.132) 6.230 4.563 ( 14.544 14.544 0.587) 20.577 4.607 ( 13.404 13.404 -1.134) 18.990 6.165 ( -9.952 -9.952 2.695) 14.330 8.092 ( -8.738 -8.738 0.771) 12.382 8.280 ( -1.997 -1.997 1.839) 3.370 8.861 ( -8.488 -8.488 1.072) 12.051 9.920 ( 14.091 14.091 3.195) 20.182 11.703 ( 17.349 17.349 -5.950) 25.247 12.223 ( -12.900 -12.900 -2.093) 18.363 12.257 ( -2.842 -2.842 10.200) 10.964 12.853 ( -14.393 -14.393 0.534) 20.361 13.182 ( -8.623 -8.623 0.601) 12.209 14.365 ( 16.367 16.367 7.308) 24.272 15.213 ( -4.725 -4.725 -3.002) 7.326 15.898 ( 7.781 7.781 -13.559) 17.463 ======================= Grid point 198 (43/84) ======================= q-point: ( 0.23 0.15 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.387 ( 19.246 12.247 -0.715) 22.823 1.550 ( 19.111 6.102 4.260) 20.509 1.935 ( 0.680 1.626 0.011) 1.763 2.079 ( 2.766 11.802 -4.807) 13.040 2.327 ( 11.159 3.084 0.059) 11.577 2.453 ( 6.881 12.026 3.020) 14.181 2.685 ( 12.958 2.577 0.580) 13.225 2.972 ( 3.265 5.741 -1.730) 6.827 3.190 ( 3.933 -0.762 -3.134) 5.087 4.661 ( 5.154 26.389 0.048) 26.887 5.108 ( 26.582 2.613 -0.397) 26.714 5.938 ( -11.968 -5.176 1.870) 13.173 8.061 ( 4.264 -24.527 1.946) 24.971 8.399 ( 12.688 -21.043 1.257) 24.604 8.642 ( -12.471 -10.871 0.557) 16.554 10.133 ( 7.727 -2.122 2.049) 8.271 11.893 ( -4.102 1.455 -0.894) 4.443 11.970 ( -4.397 -3.035 3.231) 6.243 12.314 ( 5.830 5.540 -3.148) 8.637 12.555 ( -15.408 -19.253 0.953) 24.678 13.028 ( -6.638 -9.674 0.524) 11.744 14.673 ( 14.432 19.047 5.484) 24.519 15.089 ( -7.276 -2.854 -0.085) 7.816 15.986 ( 1.195 15.018 -9.815) 17.981 ======================= Grid point 199 (44/84) ======================= q-point: ( 0.31 0.15 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.786 ( 19.846 9.892 -0.396) 22.179 1.870 ( 11.610 8.046 1.925) 14.257 1.948 ( 1.859 2.279 -0.170) 2.946 2.153 ( 5.344 9.058 -3.070) 10.956 2.484 ( 5.540 10.549 0.013) 11.915 2.610 ( 7.840 5.756 3.632) 10.382 2.895 ( 6.928 5.434 -1.874) 9.002 3.023 ( 3.877 0.006 1.570) 4.183 3.299 ( 5.799 -2.419 -4.145) 7.528 4.777 ( 6.294 26.812 -0.358) 27.543 5.607 ( 23.171 2.218 -0.259) 23.278 5.733 ( -8.184 -5.143 0.906) 9.708 8.112 ( -6.036 -9.574 2.438) 11.577 8.488 ( 5.065 -37.647 -1.095) 38.002 8.755 ( 20.873 -33.455 1.009) 39.446 10.229 ( 1.889 -16.630 -0.568) 16.747 11.672 ( -11.484 -24.948 -1.557) 27.509 11.963 ( -3.922 9.741 3.021) 10.927 12.265 ( -9.959 -13.028 1.176) 16.441 12.386 ( -2.620 3.464 -1.942) 4.758 12.959 ( 0.730 -8.150 -0.295) 8.188 14.906 ( 3.079 12.473 2.448) 13.078 14.959 ( 0.305 9.196 2.560) 9.550 15.963 ( -2.953 19.798 -6.060) 20.914 ======================= Grid point 200 (45/84) ======================= q-point: ( 0.38 0.15 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.931 ( 0.694 1.589 0.450) 1.791 2.005 ( 3.215 12.150 0.198) 12.570 2.224 ( 18.555 8.144 -1.120) 20.294 2.282 ( 7.636 7.930 -1.401) 11.097 2.569 ( 3.165 13.307 -0.817) 13.703 2.729 ( 4.154 7.385 4.277) 9.492 2.975 ( 2.081 7.938 -2.568) 8.599 3.091 ( 1.380 -0.187 2.692) 3.031 3.401 ( 4.370 -7.360 -4.109) 9.495 4.896 ( 5.402 25.980 -1.096) 26.558 5.625 ( -2.894 -7.132 0.541) 7.716 5.997 ( 15.775 1.735 -0.022) 15.871 7.880 ( -12.551 -2.020 -1.650) 12.819 8.775 ( 19.306 -49.051 0.396) 52.715 9.189 ( 21.544 -38.173 0.679) 43.839 10.193 ( -5.899 -30.549 -1.992) 31.177 11.469 ( -9.060 -36.528 -0.788) 37.643 11.778 ( -11.469 -14.742 1.140) 18.713 12.186 ( -0.591 12.134 1.499) 12.241 12.247 ( -9.067 2.341 0.224) 9.367 13.034 ( 5.201 -7.430 -0.721) 9.098 14.798 ( -5.727 5.650 -0.413) 8.056 15.105 ( 6.559 21.990 4.791) 23.443 15.893 ( -3.520 22.660 -3.707) 23.230 ======================= Grid point 201 (46/84) ======================= q-point: ( 0.46 0.15 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.939 ( 0.219 1.290 0.273) 1.337 2.039 ( 0.697 12.915 -0.200) 12.935 2.351 ( 2.024 7.297 -1.578) 7.735 2.563 ( 7.090 8.503 -4.179) 11.833 2.645 ( 9.133 14.063 2.485) 16.951 2.775 ( 0.053 11.408 4.234) 12.169 2.945 ( -5.381 6.303 0.359) 8.295 3.094 ( 0.101 2.131 0.459) 2.182 3.463 ( 1.658 -9.658 -4.416) 10.748 4.974 ( 2.176 24.953 -1.746) 25.108 5.599 ( -0.287 -8.522 0.822) 8.566 6.211 ( 5.585 1.418 0.159) 5.764 7.698 ( -5.027 -2.673 -4.177) 7.062 9.132 ( 13.307 -43.271 0.303) 45.272 9.528 ( 10.401 -37.464 0.495) 38.884 10.026 ( -8.278 -46.074 -1.790) 46.845 11.326 ( -4.522 -46.876 -0.038) 47.093 11.611 ( -4.514 -19.716 0.837) 20.243 12.076 ( -5.030 1.660 0.388) 5.311 12.207 ( 0.493 21.260 1.316) 21.306 13.115 ( 2.365 -9.278 -0.498) 9.588 14.717 ( -2.194 9.494 -1.024) 9.798 15.192 ( 2.186 22.544 4.980) 23.190 15.841 ( -1.480 24.013 -2.650) 24.204 ======================= Grid point 211 (47/84) ======================= q-point: ( 0.23 0.23 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.679 ( 16.814 16.814 -0.345) 23.781 1.761 ( 14.492 14.492 4.072) 20.895 1.961 ( 1.079 1.079 -0.219) 1.542 2.307 ( 7.752 7.752 -6.452) 12.721 2.370 ( 5.700 5.700 1.784) 8.256 2.732 ( 11.763 11.763 1.192) 16.678 2.747 ( 6.792 6.792 4.269) 10.510 2.995 ( -0.701 -0.701 -0.955) 1.377 3.268 ( 5.257 5.257 -6.035) 9.575 5.177 ( 15.866 15.866 0.444) 22.442 5.199 ( 15.678 15.678 -0.389) 22.176 5.798 ( -8.954 -8.954 1.425) 12.742 7.686 ( -11.289 -11.289 0.832) 15.986 8.131 ( -5.522 -5.522 1.221) 7.905 8.402 ( -14.757 -14.757 0.675) 20.881 10.030 ( -6.793 -6.793 1.321) 9.697 11.528 ( -21.951 -21.951 -1.812) 31.096 12.141 ( -21.675 -21.675 0.527) 30.657 12.185 ( 10.850 10.850 3.795) 15.806 12.480 ( 6.831 6.831 -2.961) 10.105 12.848 ( -8.218 -8.218 0.695) 11.642 15.001 ( 4.065 4.065 1.527) 5.949 15.051 ( 9.036 9.036 3.093) 13.148 16.205 ( 7.008 7.008 -7.026) 12.148 ======================= Grid point 212 (48/84) ======================= q-point: ( 0.31 0.23 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.948 ( 2.663 3.168 0.687) 4.196 2.061 ( 14.398 11.282 2.122) 18.414 2.066 ( 16.664 12.980 -0.294) 21.125 2.374 ( 3.166 12.160 -2.771) 12.867 2.560 ( 9.336 1.395 -1.127) 9.507 2.843 ( 5.908 12.075 3.769) 13.961 2.996 ( 9.841 6.547 0.438) 11.828 3.012 ( 4.289 0.121 -0.528) 4.323 3.336 ( 1.656 4.164 -5.000) 6.714 5.300 ( 5.675 25.591 -0.125) 26.213 5.607 ( -9.376 -7.542 0.872) 12.065 5.673 ( 22.242 4.387 0.078) 22.670 7.639 ( 6.016 -29.229 0.320) 29.844 8.058 ( -17.457 -13.633 0.089) 22.150 8.183 ( 8.722 -22.835 0.976) 24.463 9.878 ( -9.357 -17.031 -0.073) 19.432 11.107 ( -20.319 -29.398 -1.050) 35.752 11.694 ( -22.263 -25.428 0.433) 33.800 12.368 ( 6.152 11.443 1.948) 13.137 12.481 ( -6.705 4.346 -0.232) 7.994 12.770 ( 1.350 -11.313 -0.331) 11.398 14.960 ( -2.155 2.202 -0.089) 3.083 15.318 ( 12.014 17.627 3.598) 21.633 16.278 ( 1.010 11.822 -4.513) 12.694 ======================= Grid point 213 (49/84) ======================= q-point: ( 0.38 0.23 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.964 ( 0.336 1.652 0.445) 1.744 2.253 ( 5.069 12.561 0.577) 13.557 2.417 ( 13.346 11.361 -1.537) 17.594 2.486 ( 10.760 11.946 -0.345) 16.081 2.709 ( 5.856 3.299 -0.712) 6.759 2.954 ( 4.684 11.351 3.033) 12.649 3.097 ( 5.611 1.007 -1.946) 6.024 3.162 ( 4.414 11.564 1.261) 12.442 3.345 ( -0.464 1.642 -3.602) 3.986 5.399 ( 3.932 23.758 -0.484) 24.086 5.458 ( -5.110 -8.979 0.657) 10.352 6.047 ( 15.002 3.184 0.167) 15.337 7.677 ( -13.309 -28.567 -1.077) 31.533 7.948 ( 18.334 -24.184 0.327) 30.350 8.464 ( 17.615 -33.324 0.529) 37.698 9.612 ( -18.665 -25.672 -1.044) 31.757 10.732 ( -17.386 -35.644 -0.077) 39.658 11.311 ( -15.308 -29.228 0.615) 33.000 12.234 ( -12.087 -1.960 -0.570) 12.258 12.509 ( 4.383 14.312 2.476) 15.172 12.822 ( 2.552 -14.082 -0.721) 14.330 14.932 ( -0.803 7.610 -0.799) 7.694 15.507 ( 7.095 18.092 3.390) 19.727 16.269 ( -1.271 15.050 -2.890) 15.378 ======================= Grid point 214 (50/84) ======================= q-point: ( 0.46 0.23 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.969 ( 0.152 1.646 0.399) 1.701 2.304 ( 0.965 13.311 -0.299) 13.349 2.534 ( 2.367 10.698 -1.097) 11.012 2.723 ( 5.237 8.138 -3.614) 10.330 2.842 ( 10.155 5.743 2.788) 11.995 2.994 ( -1.934 9.071 2.414) 9.584 3.175 ( -1.128 18.035 0.982) 18.097 3.199 ( 3.176 7.234 -1.616) 8.064 3.326 ( -1.100 -3.371 -3.093) 4.705 5.391 ( -1.265 -10.688 0.653) 10.782 5.462 ( 1.564 22.162 -0.586) 22.224 6.249 ( 5.274 2.377 0.236) 5.789 7.466 ( -5.951 -24.658 -1.504) 25.410 8.418 ( 24.769 -25.941 -0.226) 35.868 8.778 ( 10.736 -37.161 0.200) 38.681 9.156 ( -24.398 -38.208 -1.512) 45.359 10.454 ( -8.942 -39.492 0.847) 40.500 11.114 ( -4.689 -30.344 0.867) 30.717 12.045 ( -5.586 -4.570 -0.308) 7.224 12.568 ( 1.509 15.092 1.882) 15.284 12.856 ( 0.891 -16.709 -0.707) 16.747 14.924 ( -0.135 10.863 -1.274) 10.938 15.600 ( 2.313 18.183 3.390) 18.641 16.246 ( -0.757 16.628 -2.131) 16.781 ======================= Grid point 225 (51/84) ======================= q-point: ( 0.31 0.31 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.980 ( 0.923 0.923 0.462) 1.385 2.305 ( 12.192 12.192 2.498) 17.422 2.371 ( 16.645 16.645 -0.167) 23.540 2.578 ( 4.978 4.978 0.918) 7.100 2.624 ( 8.574 8.574 -4.853) 13.061 3.017 ( 5.412 5.412 0.552) 7.674 3.025 ( 4.417 4.417 0.457) 6.263 3.216 ( 12.535 12.535 1.597) 17.798 3.402 ( 1.949 1.949 -3.861) 4.744 5.435 ( -9.678 -9.678 0.855) 13.714 5.760 ( 12.940 12.940 0.228) 18.301 5.773 ( 12.669 12.669 -0.206) 17.918 7.200 ( -14.936 -14.936 -0.628) 21.132 7.772 ( -13.934 -13.934 0.075) 19.706 7.894 ( -5.807 -5.807 0.624) 8.235 9.578 ( -14.368 -14.368 0.799) 20.335 10.558 ( -25.620 -25.620 -0.198) 36.232 11.152 ( -28.021 -28.021 -0.102) 39.628 12.465 ( -5.272 -5.272 -2.074) 7.739 12.511 ( -9.290 -9.290 0.315) 13.142 12.617 ( 7.582 7.582 3.104) 11.163 15.018 ( 3.421 3.421 -0.704) 4.890 15.616 ( 12.425 12.425 2.691) 17.776 16.438 ( 4.848 4.848 -3.055) 7.505 ======================= Grid point 226 (52/84) ======================= q-point: ( 0.38 0.31 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.994 ( 0.545 1.380 0.566) 1.588 2.489 ( 5.600 10.932 0.475) 12.292 2.638 ( 6.674 8.691 -0.388) 10.964 2.722 ( 13.796 11.348 1.840) 17.958 2.794 ( 6.793 7.859 -2.958) 10.801 3.108 ( 6.041 1.081 -1.842) 6.408 3.158 ( 6.118 10.859 1.464) 12.549 3.411 ( -0.782 3.633 -2.352) 4.398 3.441 ( 9.384 14.465 1.254) 17.288 5.248 ( -8.426 -11.288 0.854) 14.112 5.842 ( 3.003 19.779 -0.242) 20.007 6.114 ( 13.764 3.276 0.042) 14.149 6.969 ( -10.870 -34.132 -0.208) 35.822 7.657 ( 5.992 -11.618 -0.900) 13.103 7.930 ( 7.901 -19.881 0.262) 21.395 9.232 ( -22.753 -13.322 0.990) 26.384 10.110 ( -19.401 -26.608 0.734) 32.938 10.663 ( -19.388 -35.467 0.140) 40.421 12.164 ( -16.148 -4.655 -0.796) 16.825 12.479 ( -0.755 -19.671 -0.826) 19.703 12.745 ( 4.862 9.153 2.646) 10.697 15.084 ( 3.019 7.427 -1.354) 8.130 15.809 ( 7.097 12.368 2.142) 14.419 16.491 ( 1.097 7.739 -2.010) 8.071 ======================= Grid point 227 (53/84) ======================= q-point: ( 0.46 0.31 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.002 ( 0.211 1.619 0.610) 1.743 2.547 ( 1.097 10.916 -0.504) 10.983 2.727 ( 2.392 6.761 0.550) 7.192 2.881 ( 2.150 8.600 -2.581) 9.233 2.950 ( 6.916 7.832 1.694) 10.585 3.224 ( 4.277 1.276 -1.184) 4.618 3.250 ( 2.967 13.831 0.378) 14.151 3.380 ( -1.945 5.725 -1.471) 6.223 3.561 ( 2.684 18.176 0.742) 18.388 5.130 ( -3.103 -14.434 0.923) 14.792 5.883 ( 1.081 18.886 -0.414) 18.922 6.298 ( 4.767 2.358 0.034) 5.318 6.787 ( -5.699 -38.894 0.262) 39.310 8.004 ( 25.913 -19.265 -1.444) 32.322 8.107 ( 6.832 -29.802 -0.191) 30.576 8.672 ( -31.328 -10.133 0.324) 32.927 9.822 ( -8.496 -23.632 1.589) 25.163 10.423 ( -5.443 -39.203 1.019) 39.592 11.917 ( -7.259 -7.804 -0.733) 10.683 12.467 ( -0.361 -22.172 -0.863) 22.192 12.808 ( 1.572 9.353 2.305) 9.761 15.127 ( 1.174 9.361 -1.754) 9.596 15.900 ( 2.273 12.237 1.901) 12.591 16.499 ( 0.031 9.252 -1.480) 9.370 ======================= Grid point 239 (54/84) ======================= q-point: ( 0.38 0.38 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.017 ( 0.935 0.935 0.734) 1.512 2.663 ( 6.051 6.051 -0.212) 8.560 2.734 ( 3.369 3.369 -0.488) 4.790 2.954 ( 10.963 10.963 1.755) 15.604 2.984 ( 8.749 8.749 -0.397) 12.380 3.154 ( 3.686 3.686 -1.202) 5.350 3.391 ( 11.995 11.995 0.221) 16.965 3.463 ( 1.528 1.528 -1.491) 2.626 3.700 ( 11.492 11.492 0.612) 16.264 5.020 ( -11.392 -11.392 1.060) 16.146 6.161 ( 7.320 7.320 -0.263) 10.355 6.171 ( 7.365 7.365 -0.305) 10.420 6.398 ( -23.384 -23.384 0.381) 33.071 7.543 ( 2.350 2.350 -1.890) 3.823 7.680 ( -6.319 -6.319 -0.135) 8.937 8.998 ( -13.763 -13.763 2.892) 19.678 9.701 ( -14.568 -14.568 1.461) 20.654 9.982 ( -30.130 -30.130 -0.073) 42.611 12.053 ( -12.427 -12.427 -1.373) 17.627 12.093 ( -11.670 -11.670 -0.374) 16.509 12.883 ( 4.984 4.984 2.604) 7.514 15.211 ( 5.192 5.192 -1.912) 7.587 15.998 ( 6.905 6.905 1.230) 9.843 16.593 ( 3.032 3.032 -1.306) 4.483 ======================= Grid point 240 (55/84) ======================= q-point: ( 0.46 0.38 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.030 ( 0.339 1.175 0.820) 1.472 2.725 ( 1.092 6.955 -0.647) 7.070 2.797 ( 2.097 1.831 -0.360) 2.807 3.080 ( 2.385 10.014 0.156) 10.295 3.124 ( 4.319 8.289 2.012) 9.561 3.232 ( 2.887 1.085 0.488) 3.122 3.476 ( 0.081 3.742 -1.500) 4.032 3.588 ( 6.326 18.237 -0.330) 19.306 3.872 ( 4.984 12.934 -0.156) 13.862 4.833 ( -6.021 -15.376 1.117) 16.550 6.071 ( -8.814 -31.321 0.859) 32.549 6.192 ( 0.671 12.140 -0.487) 12.168 6.337 ( 4.388 1.456 -0.361) 4.637 7.623 ( -0.536 -18.806 -0.658) 18.825 7.756 ( 15.129 -2.965 -2.324) 15.591 8.608 ( -21.605 0.191 2.718) 21.776 9.523 ( -4.169 -10.551 1.756) 11.480 9.608 ( -7.607 -38.592 1.412) 39.360 11.764 ( -9.025 -6.959 -0.898) 11.432 12.015 ( -1.710 -22.204 -0.885) 22.288 12.947 ( 1.559 4.929 2.545) 5.762 15.282 ( 1.840 6.112 -2.255) 6.769 16.086 ( 2.184 6.755 0.757) 7.139 16.627 ( 0.704 4.101 -0.924) 4.262 ======================= Grid point 253 (56/84) ======================= q-point: ( 0.46 0.46 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.046 ( 0.427 0.427 0.941) 1.118 2.816 ( 1.874 1.874 -0.642) 2.726 2.823 ( 1.146 1.146 -0.595) 1.726 3.230 ( 3.411 3.411 1.925) 5.194 3.233 ( 3.330 3.330 1.851) 5.060 3.260 ( 1.280 1.280 1.978) 2.682 3.525 ( 1.268 1.268 -1.951) 2.649 3.927 ( 13.429 13.429 -0.568) 19.000 4.057 ( 5.461 5.461 -0.987) 7.786 4.544 ( -12.149 -12.149 1.321) 17.231 5.624 ( -12.299 -12.299 1.070) 17.426 6.354 ( 2.555 2.555 -0.619) 3.666 6.358 ( 2.249 2.249 -0.590) 3.235 7.371 ( -6.537 -6.537 -1.000) 9.299 7.902 ( 14.457 14.457 -2.926) 20.654 8.588 ( -5.881 -5.881 4.063) 9.257 8.856 ( -24.548 -24.548 1.361) 34.743 9.503 ( 0.930 0.930 2.402) 2.738 11.658 ( -5.903 -5.903 -0.949) 8.402 11.675 ( -7.348 -7.348 -0.818) 10.424 13.009 ( 1.518 1.518 2.645) 3.407 15.364 ( 2.103 2.103 -2.564) 3.926 16.173 ( 2.124 2.124 0.141) 3.007 16.675 ( 1.085 1.085 -0.616) 1.654 ======================= Grid point 338 (57/84) ======================= q-point: ( 0.00 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 84 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.733 ( -0.000 -0.000 11.053) 11.053 0.733 ( -0.000 -0.000 11.053) 11.053 1.190 ( 0.000 0.000 16.820) 16.820 1.377 ( 0.000 0.000 -8.532) 8.532 1.377 ( -0.000 -0.000 -8.532) 8.532 2.026 ( 0.000 -0.000 1.472) 1.472 2.026 ( 0.000 0.000 1.472) 1.472 2.349 ( 0.000 0.000 -11.868) 11.868 3.749 ( 0.000 0.000 2.034) 2.034 3.912 ( 0.000 0.000 0.330) 0.330 3.912 ( -0.000 -0.000 0.330) 0.330 6.983 ( 0.000 0.000 0.030) 0.030 8.326 ( 0.000 0.000 1.253) 1.253 8.326 ( 0.000 0.000 1.253) 1.253 8.397 ( 0.000 0.000 -1.548) 1.548 8.397 ( 0.000 0.000 -1.548) 1.548 9.182 ( 0.000 0.000 0.462) 0.462 12.645 ( 0.000 0.000 -0.204) 0.204 13.412 ( 0.000 0.000 -0.315) 0.315 13.412 ( 0.000 0.000 -0.315) 0.315 13.628 ( -0.000 0.000 0.419) 0.419 13.628 ( 0.000 0.000 0.419) 0.419 15.473 ( 0.000 0.000 1.635) 1.635 15.759 ( -0.000 -0.000 -0.314) 0.314 ======================= Grid point 339 (58/84) ======================= q-point: ( 0.08 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.879 ( -1.395 0.000 15.612) 15.674 0.952 ( 18.562 0.000 9.217) 20.724 1.173 ( 18.061 0.000 -4.783) 18.683 1.426 ( 4.620 0.000 -8.620) 9.780 1.634 ( 18.269 0.000 0.534) 18.276 2.034 ( 0.830 0.000 1.409) 1.635 2.101 ( 0.838 0.000 5.833) 5.893 2.307 ( 5.606 0.000 -10.392) 11.808 3.598 ( -9.816 0.000 1.592) 9.944 3.916 ( 0.353 0.000 0.293) 0.459 4.108 ( 17.100 0.000 -0.428) 17.106 6.770 ( -17.636 0.000 1.122) 17.671 8.440 ( 11.295 0.000 1.130) 11.351 8.524 ( 12.707 0.000 -1.391) 12.783 8.642 ( 20.344 0.000 0.031) 20.344 9.156 ( -1.714 0.000 0.537) 1.797 9.372 ( 77.055 0.000 -10.228) 77.731 12.603 ( -3.701 0.000 -0.373) 3.720 13.302 ( -13.559 0.000 2.544) 13.796 13.374 ( -3.822 0.000 -0.249) 3.831 13.598 ( -3.011 0.000 0.354) 3.031 13.846 ( 19.585 0.000 6.132) 20.522 15.038 ( -32.777 0.000 -4.377) 33.068 15.713 ( -4.449 0.000 -0.276) 4.457 ======================= Grid point 340 (59/84) ======================= q-point: ( 0.15 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.876 ( 6.721 0.000 4.016) 7.830 1.333 ( 17.538 0.000 8.372) 19.433 1.550 ( 8.196 0.000 -9.164) 12.294 1.642 ( 20.219 0.000 -4.701) 20.758 2.003 ( 17.302 0.000 6.474) 18.473 2.058 ( 1.599 0.000 1.199) 1.999 2.145 ( 9.690 0.000 -3.589) 10.333 2.558 ( 15.422 0.000 -1.148) 15.465 3.442 ( -5.447 0.000 0.793) 5.505 3.925 ( 0.543 0.000 0.202) 0.580 4.531 ( 24.757 0.000 -1.077) 24.780 6.398 ( -18.557 0.000 1.484) 18.616 8.759 ( 20.527 0.000 0.876) 20.545 8.871 ( -4.220 0.000 2.585) 4.948 8.884 ( 23.207 0.000 -1.049) 23.231 9.296 ( 20.389 0.000 0.049) 20.389 10.650 ( 48.713 0.000 -7.868) 49.344 12.523 ( -4.217 0.000 0.180) 4.221 12.805 ( -33.916 0.000 7.497) 34.735 13.261 ( -7.739 0.000 -0.055) 7.739 13.510 ( -5.827 0.000 0.156) 5.829 14.195 ( 10.345 0.000 20.073) 22.582 14.682 ( 1.934 0.000 -24.291) 24.368 15.584 ( -8.616 0.000 -1.092) 8.684 ======================= Grid point 341 (60/84) ======================= q-point: ( 0.23 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.145 ( 18.911 0.000 -0.321) 18.914 1.568 ( 6.854 0.000 2.927) 7.453 1.773 ( 13.171 0.000 -3.848) 13.721 1.850 ( 4.048 0.000 -1.149) 4.208 2.096 ( 2.087 0.000 0.845) 2.251 2.351 ( 16.783 0.000 1.297) 16.833 2.419 ( 13.103 0.000 -2.286) 13.300 2.870 ( 15.664 0.000 1.245) 15.713 3.411 ( 3.246 0.000 -0.020) 3.246 3.936 ( 0.517 0.000 0.105) 0.527 5.050 ( 26.888 0.000 -0.948) 26.905 6.070 ( -14.147 0.000 1.209) 14.199 8.634 ( -17.287 0.000 2.789) 17.510 9.225 ( 25.881 0.000 0.627) 25.889 9.415 ( 29.908 0.000 -0.683) 29.916 9.746 ( 22.128 0.000 0.482) 22.133 11.352 ( 24.620 0.000 -2.071) 24.707 12.115 ( -29.868 0.000 5.757) 30.418 12.512 ( 1.291 0.000 -2.884) 3.160 13.068 ( -11.789 0.000 0.235) 11.791 13.370 ( -8.204 0.000 -0.151) 8.205 14.330 ( 6.849 0.000 9.911) 12.047 14.851 ( 7.674 0.000 -10.858) 13.296 15.382 ( -11.362 0.000 -2.779) 11.697 ======================= Grid point 342 (61/84) ======================= q-point: ( 0.31 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.568 ( 22.994 0.000 -0.852) 23.010 1.663 ( 3.342 0.000 0.243) 3.351 1.895 ( 1.236 0.000 -0.077) 1.239 2.004 ( 9.674 0.000 -0.860) 9.712 2.137 ( 2.022 0.000 0.449) 2.072 2.572 ( 5.177 0.000 -2.026) 5.559 2.625 ( 8.479 0.000 1.051) 8.544 3.136 ( 10.015 0.000 2.717) 10.377 3.577 ( 12.286 0.000 -2.333) 12.506 3.944 ( 0.352 0.000 0.037) 0.354 5.559 ( 23.895 0.000 -0.576) 23.902 5.856 ( -7.218 0.000 0.596) 7.242 8.235 ( -21.946 0.000 1.116) 21.974 9.747 ( 26.081 0.000 0.428) 26.085 10.037 ( 32.379 0.000 -0.355) 32.381 10.133 ( 16.919 0.000 0.223) 16.920 11.603 ( -12.116 0.000 4.774) 13.023 11.717 ( 3.664 0.000 -3.174) 4.847 12.625 ( 9.806 0.000 -2.212) 10.052 12.795 ( -15.856 0.000 0.543) 15.865 13.192 ( -9.622 0.000 -0.497) 9.635 14.482 ( 7.907 0.000 6.206) 10.052 14.880 ( -5.331 0.000 -1.295) 5.486 15.173 ( -8.419 0.000 -5.352) 9.977 ======================= Grid point 343 (62/84) ======================= q-point: ( 0.38 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.711 ( 1.654 0.000 -0.101) 1.658 1.908 ( 0.905 0.000 0.477) 1.023 2.028 ( 22.123 0.000 -0.965) 22.144 2.151 ( 5.131 0.000 -0.001) 5.131 2.173 ( 1.688 0.000 -0.057) 1.688 2.621 ( 0.936 0.000 0.290) 0.980 2.741 ( 3.123 0.000 0.082) 3.124 3.239 ( 0.561 0.000 3.490) 3.534 3.823 ( 11.055 0.000 -4.943) 12.110 3.949 ( 0.175 0.000 0.005) 0.175 5.781 ( -0.980 0.000 0.265) 1.015 5.964 ( 16.425 0.000 -0.233) 16.427 7.822 ( -18.361 0.000 -1.746) 18.444 10.216 ( 20.329 0.000 0.301) 20.331 10.411 ( 11.135 0.000 0.136) 11.135 10.667 ( 30.604 0.000 -0.121) 30.605 11.340 ( -10.410 0.000 -0.281) 10.414 11.825 ( 3.697 0.000 0.351) 3.713 12.446 ( -19.193 0.000 0.803) 19.210 12.858 ( 12.372 0.000 -0.955) 12.409 13.005 ( -8.822 0.000 -0.796) 8.857 14.620 ( 5.329 0.000 6.156) 8.142 14.688 ( -11.356 0.000 -0.066) 11.356 15.055 ( -4.426 0.000 -4.658) 6.425 ======================= Grid point 344 (63/84) ======================= q-point: ( 0.46 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.731 ( 0.500 0.000 -0.101) 0.510 1.918 ( 0.248 0.000 0.165) 0.298 2.189 ( 0.526 0.000 0.281) 0.596 2.223 ( 1.690 0.000 -0.322) 1.720 2.410 ( 14.107 0.000 -0.650) 14.122 2.635 ( 1.211 0.000 0.343) 1.258 2.736 ( -3.828 0.000 2.037) 4.336 3.200 ( -2.624 0.000 2.455) 3.594 3.952 ( 0.049 0.000 -0.003) 0.049 3.977 ( 4.117 0.000 -6.152) 7.403 5.785 ( 0.599 0.000 0.259) 0.652 6.187 ( 5.846 0.000 -0.035) 5.846 7.560 ( -7.150 0.000 -3.825) 8.109 10.506 ( 8.008 0.000 0.254) 8.012 10.566 ( 4.230 0.000 0.567) 4.267 11.192 ( 19.509 0.000 -0.140) 19.510 11.199 ( -3.761 0.000 -2.002) 4.261 11.871 ( 1.105 0.000 1.037) 1.515 12.075 ( -15.451 0.000 1.089) 15.490 12.871 ( -3.893 0.000 -0.979) 4.015 13.055 ( 6.059 0.000 0.117) 6.060 14.493 ( -5.957 0.000 -0.623) 5.990 14.700 ( 1.971 0.000 6.507) 6.798 14.991 ( -1.782 0.000 -3.921) 4.307 ======================= Grid point 352 (64/84) ======================= q-point: ( 0.08 0.08 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.854 ( -3.075 -3.075 7.985) 9.092 0.936 ( 8.957 8.957 9.080) 15.586 1.481 ( 21.992 21.992 -1.964) 31.164 1.512 ( 6.303 6.303 -8.263) 12.155 1.736 ( 10.802 10.802 1.032) 15.311 2.091 ( -0.045 -0.045 2.778) 2.779 2.097 ( 3.783 3.783 1.242) 5.493 2.369 ( 8.241 8.241 -5.706) 12.976 3.488 ( -8.065 -8.065 1.277) 11.477 4.096 ( 8.979 8.979 0.363) 12.703 4.130 ( 9.020 9.020 -0.776) 12.780 6.651 ( -11.552 -11.552 1.358) 16.394 8.321 ( -2.052 -2.052 1.041) 3.083 8.351 ( -1.024 -1.024 -0.649) 1.587 9.075 ( 26.454 26.454 0.344) 37.413 9.116 ( -3.481 -3.481 0.705) 4.974 10.112 ( 53.250 53.250 -10.317) 76.010 12.543 ( -5.459 -5.459 -0.616) 7.745 13.195 ( -15.294 -15.294 6.192) 22.498 13.277 ( -6.895 -6.895 -0.100) 9.752 13.499 ( -6.148 -6.148 0.358) 8.702 14.097 ( 17.941 17.941 14.390) 29.170 14.844 ( -11.765 -11.765 -16.033) 23.106 15.678 ( -3.477 -3.477 -0.684) 4.965 ======================= Grid point 353 (65/84) ======================= q-point: ( 0.15 0.08 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.899 ( 9.242 3.632 2.174) 10.165 1.254 ( 19.640 -1.887 6.821) 20.876 1.628 ( 5.782 8.821 -7.510) 12.948 1.831 ( 9.004 9.425 -1.182) 13.088 2.019 ( 16.409 3.097 3.812) 17.128 2.169 ( 4.818 8.576 0.984) 9.886 2.199 ( 12.904 7.430 -3.605) 15.321 2.602 ( 13.610 4.601 -0.787) 14.389 3.344 ( -5.841 -7.406 0.706) 9.459 4.136 ( 1.763 18.488 0.058) 18.572 4.540 ( 24.826 1.060 -0.744) 24.860 6.359 ( -16.508 -4.272 1.504) 17.118 8.418 ( 10.622 -22.300 0.328) 24.703 8.517 ( 15.867 -16.539 0.865) 22.935 8.993 ( -7.754 -2.280 0.925) 8.135 9.607 ( 25.607 22.070 0.979) 33.819 11.035 ( 37.461 30.455 -5.746) 48.619 12.420 ( -6.355 -9.705 -0.467) 11.610 12.695 ( -30.475 -10.887 7.497) 33.218 13.126 ( -9.105 -11.830 0.290) 14.930 13.376 ( -6.574 -9.778 0.236) 11.785 14.345 ( 7.128 10.923 17.884) 22.135 14.851 ( 10.543 14.935 -21.715) 28.386 15.581 ( -6.092 0.277 -2.482) 6.584 ======================= Grid point 354 (66/84) ======================= q-point: ( 0.23 0.08 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.194 ( 19.260 5.460 -0.457) 20.024 1.606 ( 13.401 1.384 6.084) 14.782 1.753 ( 7.840 4.485 -7.316) 11.623 1.914 ( 2.443 2.774 0.312) 3.710 2.198 ( 3.731 9.295 1.395) 10.113 2.398 ( 14.581 1.339 1.479) 14.717 2.491 ( 13.195 9.739 -3.221) 16.714 2.876 ( 13.715 0.786 0.611) 13.751 3.286 ( 0.744 -8.881 -0.233) 8.915 4.183 ( 3.097 20.861 0.074) 21.090 5.059 ( 26.868 1.099 -0.693) 26.899 6.048 ( -14.081 -2.283 1.173) 14.313 8.544 ( -6.897 -12.194 2.587) 14.246 8.916 ( 22.359 -25.694 1.238) 34.083 8.954 ( 12.676 -25.061 -1.114) 28.107 10.062 ( 20.178 13.964 0.691) 24.548 11.588 ( 18.872 19.244 -0.854) 26.967 12.117 ( -22.924 0.622 4.915) 23.453 12.365 ( -0.937 -13.634 -3.596) 14.132 12.911 ( -12.645 -14.112 0.443) 18.954 13.235 ( -7.553 -10.342 0.050) 12.807 14.490 ( 8.355 13.939 7.883) 18.062 15.013 ( 2.418 9.826 -5.594) 11.563 15.461 ( -5.256 10.133 -5.860) 12.832 ======================= Grid point 355 (67/84) ======================= q-point: ( 0.31 0.08 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.615 ( 22.217 5.160 -0.743) 22.821 1.753 ( 4.233 8.163 1.049) 9.255 1.895 ( 2.944 0.692 -0.843) 3.140 2.011 ( 8.559 0.826 -1.002) 8.657 2.258 ( 2.480 11.323 0.830) 11.621 2.625 ( 7.871 0.754 1.637) 8.075 2.669 ( 5.687 8.318 -2.372) 10.352 3.113 ( 9.240 -2.131 2.418) 9.787 3.378 ( 7.499 -12.168 -2.396) 14.493 4.259 ( 4.587 23.716 -0.269) 24.157 5.567 ( 23.818 0.943 -0.428) 23.841 5.829 ( -7.714 -2.745 0.576) 8.208 8.227 ( -20.321 -1.021 1.158) 20.379 9.327 ( 22.358 -43.539 0.043) 48.944 9.421 ( 25.621 -29.021 0.396) 38.715 10.405 ( 14.734 1.070 -0.381) 14.778 11.747 ( -7.486 13.383 4.261) 15.916 11.835 ( 0.867 7.050 -2.265) 7.456 12.425 ( 4.648 -16.163 -2.402) 16.989 12.630 ( -15.720 -14.915 0.543) 21.677 13.100 ( -4.952 -6.941 -0.270) 8.530 14.664 ( 8.697 16.195 5.306) 19.133 14.929 ( -9.253 1.935 0.109) 9.454 15.385 ( -3.041 20.319 -5.982) 21.398 ======================= Grid point 356 (68/84) ======================= q-point: ( 0.38 0.08 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.796 ( 1.802 8.238 -0.046) 8.433 1.916 ( 0.753 0.829 0.391) 1.186 2.074 ( 21.557 4.692 -0.923) 22.081 2.162 ( 5.774 1.524 -0.676) 6.010 2.298 ( 1.753 11.987 0.443) 12.123 2.700 ( 1.380 6.426 2.334) 6.974 2.769 ( 3.510 3.691 -1.708) 5.372 3.212 ( 0.636 -2.631 3.274) 4.248 3.521 ( 5.990 -14.761 -3.822) 16.382 4.356 ( 4.759 25.237 -1.152) 25.707 5.743 ( -1.449 -3.780 0.262) 4.056 5.971 ( 16.358 0.755 -0.155) 16.376 7.837 ( -17.607 1.215 -1.545) 17.717 9.742 ( 18.305 -45.930 0.354) 49.444 9.895 ( 21.133 -29.172 0.355) 36.024 10.649 ( 10.042 -16.323 -1.265) 19.206 11.630 ( -1.936 22.848 0.843) 22.945 11.883 ( 2.720 -0.612 0.049) 2.789 12.296 ( -17.242 -12.808 0.986) 21.502 12.433 ( -5.371 -26.118 -0.984) 26.683 13.094 ( 4.398 0.276 -0.574) 4.444 14.714 ( -10.514 2.950 0.396) 10.927 14.818 ( 5.888 17.118 4.677) 18.697 15.322 ( -3.257 23.141 -4.019) 23.712 ======================= Grid point 357 (69/84) ======================= q-point: ( 0.46 0.08 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.818 ( 0.521 8.354 -0.107) 8.371 1.924 ( 0.207 0.681 0.158) 0.729 2.235 ( 1.778 1.956 -0.522) 2.694 2.305 ( 0.930 10.083 -0.566) 10.141 2.467 ( 14.006 6.438 -0.114) 15.416 2.717 ( 0.822 6.887 2.844) 7.497 2.774 ( -3.349 4.625 -0.165) 5.712 3.173 ( -2.652 -2.612 2.123) 4.286 3.600 ( 1.945 -15.397 -3.887) 15.999 4.428 ( 2.075 25.225 -1.968) 25.386 5.741 ( 0.442 -4.429 0.280) 4.460 6.193 ( 5.820 0.638 0.012) 5.855 7.585 ( -6.897 2.180 -3.493) 8.033 10.002 ( 7.074 -41.582 0.319) 42.181 10.203 ( 8.690 -27.553 0.314) 28.893 10.792 ( 3.946 -32.214 -1.585) 32.493 11.678 ( 6.114 25.221 0.336) 25.954 11.950 ( 3.143 2.594 0.758) 4.146 11.975 ( -13.290 -7.255 0.789) 15.162 12.284 ( -6.871 -39.075 -0.467) 39.677 13.197 ( 3.839 1.099 -0.230) 4.000 14.548 ( -4.981 5.398 -0.616) 7.370 14.901 ( 2.167 17.859 5.450) 18.797 15.271 ( -1.498 24.063 -3.286) 24.332 ======================= Grid point 366 (70/84) ======================= q-point: ( 0.15 0.15 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.056 ( 12.405 12.405 0.095) 17.544 1.358 ( 12.010 12.010 6.605) 18.224 1.820 ( 8.603 8.603 -7.248) 14.162 1.921 ( 1.726 1.726 0.667) 2.531 2.134 ( 8.738 8.738 1.399) 12.437 2.339 ( 8.480 8.480 0.653) 12.010 2.447 ( 14.663 14.663 -0.322) 20.739 2.739 ( 8.797 8.797 -2.420) 12.674 3.206 ( -5.917 -5.917 0.243) 8.371 4.574 ( 14.502 14.502 0.366) 20.512 4.585 ( 13.946 13.946 -0.686) 19.734 6.215 ( -10.391 -10.391 1.439) 14.765 8.104 ( -8.180 -8.180 0.288) 11.572 8.319 ( -2.006 -2.006 1.260) 3.104 8.883 ( -8.682 -8.682 0.746) 12.301 9.982 ( 14.539 14.539 1.893) 20.649 11.616 ( 25.881 25.881 -2.084) 36.660 12.182 ( -13.217 -13.217 -1.284) 18.735 12.426 ( -14.083 -14.083 4.550) 20.429 12.862 ( -14.176 -14.176 0.275) 20.050 13.194 ( -8.531 -8.531 0.384) 12.071 14.567 ( 11.696 11.696 9.548) 19.098 15.096 ( 6.032 6.032 -7.203) 11.164 15.636 ( 6.091 6.091 -7.675) 11.537 ======================= Grid point 367 (71/84) ======================= q-point: ( 0.23 0.15 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.372 ( 18.900 12.445 -0.471) 22.634 1.662 ( 16.711 6.449 5.235) 18.661 1.931 ( 1.950 7.106 -2.317) 7.724 1.956 ( 3.603 5.697 -3.342) 7.524 2.344 ( 10.710 2.658 0.846) 11.067 2.501 ( 8.332 11.587 1.218) 14.324 2.710 ( 11.118 10.185 1.544) 15.157 2.916 ( 8.275 3.554 -2.765) 9.421 3.136 ( -0.353 -5.503 -1.398) 5.689 4.662 ( 5.035 26.364 0.030) 26.840 5.101 ( 26.649 3.304 -0.206) 26.854 5.974 ( -12.886 -5.471 1.063) 14.040 8.094 ( 5.910 -26.320 0.894) 26.990 8.431 ( 12.297 -20.910 1.239) 24.290 8.651 ( -14.228 -10.069 0.191) 17.432 10.173 ( 5.347 -1.921 1.258) 5.819 11.865 ( -9.371 -7.082 -1.049) 11.793 12.068 ( 7.912 4.540 4.785) 10.301 12.230 ( -2.819 5.362 -4.178) 7.359 12.572 ( -14.956 -19.402 0.478) 24.502 13.038 ( -6.972 -9.379 0.331) 11.692 14.805 ( 12.119 16.799 5.163) 21.348 15.085 ( -5.192 -1.068 -0.219) 5.306 15.772 ( 5.606 17.477 -7.618) 19.872 ======================= Grid point 368 (72/84) ======================= q-point: ( 0.31 0.15 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.776 ( 21.202 10.983 -0.344) 23.880 1.917 ( 5.244 5.070 1.483) 7.443 1.960 ( 3.361 7.339 1.828) 8.277 2.068 ( 8.620 6.775 -4.178) 11.733 2.501 ( 5.180 8.939 0.876) 10.369 2.675 ( 8.652 7.423 1.977) 11.570 2.860 ( 4.261 9.222 -1.127) 10.222 3.064 ( 6.905 -2.476 2.204) 7.660 3.198 ( 5.412 -5.004 -4.299) 8.533 4.770 ( 5.697 27.036 -0.185) 27.630 5.602 ( 23.376 2.704 -0.124) 23.533 5.750 ( -8.944 -5.199 0.537) 10.360 8.159 ( -9.822 -8.764 1.480) 13.247 8.468 ( 6.729 -39.132 -0.550) 39.710 8.772 ( 20.626 -33.399 0.481) 39.258 10.218 ( -0.408 -16.374 -0.348) 16.383 11.642 ( -10.770 -24.930 -0.911) 27.172 12.028 ( -8.466 5.248 2.306) 10.224 12.292 ( -3.876 -4.388 1.308) 5.999 12.336 ( -3.100 -0.253 -2.321) 3.881 12.953 ( -0.341 -8.137 -0.150) 8.145 14.943 ( -7.632 0.939 0.865) 7.739 15.035 ( 10.080 18.958 3.268) 21.718 15.829 ( 0.447 21.602 -4.758) 22.124 ======================= Grid point 369 (73/84) ======================= q-point: ( 0.38 0.15 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.940 ( 0.367 1.514 0.301) 1.586 2.010 ( 2.197 12.638 0.172) 12.829 2.207 ( 18.411 8.365 -0.222) 20.224 2.242 ( 8.253 7.296 -1.789) 11.160 2.576 ( 3.219 13.105 0.736) 13.515 2.813 ( 4.994 7.293 3.219) 9.407 2.921 ( 2.494 9.629 -2.329) 10.216 3.152 ( 0.939 -2.639 2.828) 3.980 3.300 ( 4.312 -7.328 -4.436) 9.590 4.875 ( 4.653 26.793 -0.609) 27.200 5.635 ( -2.903 -7.085 0.279) 7.662 5.997 ( 16.002 2.007 0.010) 16.127 7.847 ( -15.924 -1.237 -1.065) 16.008 8.784 ( 19.679 -49.114 0.277) 52.911 9.202 ( 21.350 -38.222 0.385) 43.782 10.154 ( -6.376 -30.033 -1.216) 30.727 11.455 ( -8.197 -36.200 -0.410) 37.118 11.800 ( -12.066 -15.494 0.675) 19.649 12.221 ( -3.975 7.652 1.530) 8.758 12.246 ( -5.109 7.707 -0.496) 9.259 13.020 ( 5.320 -7.383 -0.465) 9.112 14.793 ( -6.720 5.247 -0.106) 8.526 15.208 ( 6.908 20.376 3.593) 21.813 15.810 ( -1.819 23.679 -3.060) 23.945 ======================= Grid point 370 (74/84) ======================= q-point: ( 0.46 0.15 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.945 ( 0.123 1.354 0.176) 1.371 2.035 ( 0.589 12.875 -0.129) 12.889 2.322 ( 2.027 7.196 -0.862) 7.525 2.491 ( 5.472 8.236 -2.016) 10.092 2.699 ( 10.906 14.661 1.671) 18.349 2.860 ( -1.248 9.676 3.378) 10.324 2.943 ( -0.828 10.698 -0.743) 10.755 3.109 ( -3.244 -3.169 0.814) 4.607 3.360 ( 1.647 -8.849 -3.904) 9.811 4.942 ( 1.827 26.313 -0.947) 26.393 5.613 ( -0.083 -8.522 0.332) 8.529 6.215 ( 5.681 1.571 0.111) 5.896 7.617 ( -6.355 -0.505 -2.513) 6.852 9.138 ( 12.966 -43.467 0.182) 45.360 9.537 ( 10.263 -37.657 0.282) 39.031 9.990 ( -7.827 -45.705 -1.142) 46.385 11.328 ( -3.949 -46.244 0.118) 46.413 11.626 ( -4.725 -19.652 0.427) 20.216 12.084 ( -4.850 0.842 0.284) 4.931 12.232 ( 0.132 21.945 0.763) 21.959 13.105 ( 2.543 -9.476 -0.320) 9.816 14.697 ( -2.607 9.380 -0.631) 9.756 15.301 ( 2.380 20.657 3.862) 21.149 15.778 ( -1.004 24.729 -2.438) 24.869 ======================= Grid point 380 (75/84) ======================= q-point: ( 0.23 0.23 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.671 ( 17.372 17.372 -0.250) 24.569 1.870 ( 13.716 13.716 5.164) 20.073 1.959 ( 0.719 0.719 -0.019) 1.018 2.155 ( 8.226 8.226 -6.382) 13.269 2.405 ( 5.472 5.472 1.077) 7.813 2.751 ( 11.965 11.965 0.497) 16.929 2.867 ( 6.602 6.602 6.285) 11.255 2.972 ( 0.694 0.694 -1.236) 1.578 3.117 ( 4.025 4.025 -6.750) 8.830 5.186 ( 15.780 15.780 0.280) 22.319 5.191 ( 15.855 15.855 -0.228) 22.423 5.825 ( -9.471 -9.471 0.834) 13.420 7.699 ( -11.812 -11.812 0.358) 16.709 8.159 ( -5.921 -5.921 1.028) 8.437 8.413 ( -15.014 -15.014 0.304) 21.236 10.056 ( -8.157 -8.157 0.819) 11.565 11.492 ( -21.372 -21.372 -1.098) 30.244 12.150 ( -21.883 -21.883 0.290) 30.948 12.286 ( 9.455 9.455 5.296) 14.382 12.395 ( 8.060 8.060 -4.819) 12.375 12.861 ( -8.276 -8.276 0.427) 11.712 15.029 ( -3.284 -3.284 0.827) 4.718 15.127 ( 14.911 14.911 3.018) 21.303 16.053 ( 10.270 10.270 -5.294) 15.459 ======================= Grid point 381 (76/84) ======================= q-point: ( 0.31 0.23 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.964 ( 1.021 1.939 0.571) 2.265 2.045 ( 18.879 14.870 -0.721) 24.043 2.148 ( 11.768 11.555 5.139) 17.274 2.275 ( 5.927 12.047 -5.891) 14.662 2.559 ( 8.650 0.547 0.479) 8.680 2.927 ( 5.847 12.610 3.122) 14.246 2.969 ( 3.339 2.162 -1.480) 4.244 3.051 ( 8.614 3.409 3.191) 9.798 3.208 ( 3.869 4.823 -6.071) 8.666 5.298 ( 5.375 25.787 -0.062) 26.341 5.624 ( -9.857 -7.564 0.520) 12.436 5.675 ( 22.345 4.490 0.051) 22.792 7.644 ( 5.391 -30.893 0.099) 31.360 8.065 ( -18.448 -13.021 0.356) 22.584 8.198 ( 8.814 -23.019 0.390) 24.652 9.877 ( -10.601 -16.270 -0.047) 19.419 11.086 ( -19.432 -28.652 -0.635) 34.626 11.703 ( -22.233 -26.164 0.264) 34.335 12.408 ( -1.875 5.958 1.392) 6.399 12.473 ( 2.638 10.126 -0.394) 10.472 12.764 ( -0.043 -11.349 -0.184) 11.351 14.963 ( -2.999 1.574 0.204) 3.394 15.393 ( 11.578 16.425 2.511) 20.252 16.181 ( 3.288 13.662 -3.305) 14.436 ======================= Grid point 382 (77/84) ======================= q-point: ( 0.38 0.23 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.972 ( 0.287 1.701 0.279) 1.747 2.267 ( 2.691 12.805 0.516) 13.095 2.406 ( 12.757 12.045 0.715) 17.560 2.447 ( 11.998 12.673 -2.680) 17.656 2.715 ( 7.379 2.293 0.768) 7.765 3.022 ( 3.991 9.487 2.644) 10.627 3.054 ( 5.031 2.266 -1.556) 5.733 3.191 ( 4.049 10.093 1.417) 10.966 3.258 ( 1.155 3.308 -3.850) 5.206 5.391 ( 3.824 24.725 -0.207) 25.020 5.469 ( -5.323 -9.379 0.308) 10.789 6.050 ( 15.089 3.223 0.108) 15.430 7.656 ( -13.681 -27.182 -0.653) 30.438 7.954 ( 17.404 -25.459 0.185) 30.840 8.475 ( 17.304 -33.524 0.317) 37.728 9.591 ( -19.345 -24.187 -0.646) 30.978 10.731 ( -16.368 -34.861 -0.038) 38.513 11.323 ( -15.011 -29.669 0.384) 33.252 12.225 ( -11.678 -2.298 -0.238) 11.904 12.557 ( 3.656 14.675 1.478) 15.195 12.806 ( 2.582 -14.188 -0.523) 14.431 14.917 ( -1.579 7.083 -0.408) 7.268 15.578 ( 7.094 16.554 2.388) 18.167 16.207 ( -0.064 16.155 -2.158) 16.299 ======================= Grid point 383 (78/84) ======================= q-point: ( 0.46 0.23 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.977 ( 0.116 1.825 0.248) 1.846 2.298 ( 0.727 13.148 -0.183) 13.169 2.519 ( 2.285 12.401 -0.250) 12.612 2.656 ( 4.685 8.383 -2.120) 9.835 2.897 ( 11.109 5.003 1.862) 12.325 3.050 ( -2.386 5.320 2.212) 6.236 3.151 ( 3.698 6.909 -1.702) 8.019 3.202 ( -1.674 18.761 0.729) 18.850 3.256 ( -0.901 -0.739 -2.730) 2.969 5.407 ( -1.086 -10.439 0.665) 10.516 5.447 ( 1.394 22.605 -0.639) 22.657 6.254 ( 5.313 2.328 0.150) 5.803 7.437 ( -6.241 -21.954 -0.925) 22.842 8.413 ( 24.729 -26.851 -0.147) 36.504 8.782 ( 10.456 -37.494 0.128) 38.924 9.126 ( -24.836 -37.566 -0.927) 45.043 10.471 ( -8.303 -38.721 0.517) 39.605 11.132 ( -4.543 -30.173 0.538) 30.518 12.040 ( -5.540 -5.144 -0.150) 7.561 12.605 ( 1.212 15.621 1.152) 15.710 12.842 ( 0.928 -16.915 -0.455) 16.947 14.898 ( -0.428 10.460 -0.797) 10.499 15.672 ( 2.366 16.384 2.452) 16.734 16.199 ( -0.396 17.476 -1.683) 17.562 ======================= Grid point 394 (79/84) ======================= q-point: ( 0.31 0.31 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.989 ( 0.986 0.986 0.281) 1.422 2.367 ( 16.871 16.871 -0.117) 23.859 2.379 ( 9.911 9.911 4.249) 14.646 2.506 ( 11.248 11.248 -5.522) 16.838 2.599 ( 4.804 4.804 0.744) 6.834 2.980 ( 2.222 2.222 -1.359) 3.423 3.098 ( 5.378 5.378 3.470) 8.360 3.244 ( 12.547 12.547 0.800) 17.762 3.307 ( 4.138 4.138 -4.019) 7.100 5.451 ( -9.702 -9.702 0.518) 13.731 5.764 ( 12.776 12.776 0.141) 18.068 5.768 ( 12.697 12.697 -0.131) 17.957 7.187 ( -15.152 -15.152 -0.377) 21.432 7.774 ( -14.622 -14.622 0.048) 20.678 7.906 ( -6.111 -6.111 0.372) 8.650 9.594 ( -13.593 -13.593 0.496) 19.229 10.554 ( -24.719 -24.719 -0.122) 34.958 11.150 ( -28.296 -28.296 -0.061) 40.017 12.432 ( -5.818 -5.818 -0.837) 8.271 12.517 ( -9.598 -9.598 0.193) 13.574 12.669 ( 8.230 8.230 1.454) 11.730 15.006 ( 2.642 2.642 -0.296) 3.748 15.671 ( 11.646 11.646 1.791) 16.567 16.374 ( 6.254 6.254 -2.140) 9.100 ======================= Grid point 395 (80/84) ======================= q-point: ( 0.38 0.31 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.005 ( 0.636 1.550 0.350) 1.711 2.498 ( 3.192 10.269 0.260) 10.757 2.645 ( 8.936 8.594 0.928) 12.432 2.719 ( 9.123 14.132 -2.458) 17.000 2.766 ( 11.949 7.100 0.905) 13.929 3.069 ( 6.627 1.512 -1.291) 6.919 3.195 ( 4.161 9.174 1.496) 10.185 3.358 ( 0.937 5.118 -2.048) 5.592 3.463 ( 8.950 14.164 0.666) 16.768 5.264 ( -8.384 -11.045 0.517) 13.877 5.837 ( 2.759 19.840 -0.151) 20.031 6.114 ( 13.770 3.001 0.030) 14.093 6.965 ( -10.141 -33.633 -0.138) 35.129 7.640 ( 5.047 -12.692 -0.530) 13.669 7.935 ( 7.478 -20.249 0.141) 21.586 9.251 ( -23.254 -10.979 0.611) 25.723 10.124 ( -18.462 -25.798 0.452) 31.727 10.666 ( -18.574 -35.965 0.096) 40.478 12.149 ( -15.985 -4.961 -0.464) 16.744 12.465 ( -0.981 -19.730 -0.413) 19.759 12.795 ( 4.571 9.345 1.515) 10.513 15.058 ( 2.382 6.789 -0.802) 7.239 15.852 ( 6.731 11.176 1.424) 13.124 16.449 ( 1.926 8.675 -1.391) 8.994 ======================= Grid point 396 (81/84) ======================= q-point: ( 0.46 0.31 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.014 ( 0.225 1.847 0.381) 1.899 2.537 ( 0.921 10.698 -0.319) 10.742 2.750 ( 2.498 5.302 1.523) 6.056 2.821 ( 2.242 12.348 -2.608) 12.818 2.983 ( 6.821 7.223 1.017) 9.987 3.200 ( 5.042 1.504 -0.733) 5.312 3.258 ( 2.107 14.107 0.281) 14.266 3.349 ( -1.354 6.336 -1.002) 6.556 3.574 ( 2.295 17.428 0.380) 17.583 5.148 ( -3.039 -14.252 0.561) 14.583 5.875 ( 0.994 18.948 -0.255) 18.976 6.299 ( 4.757 2.018 0.025) 5.168 6.792 ( -5.487 -37.949 0.159) 38.343 7.975 ( 25.263 -20.604 -0.919) 32.613 8.103 ( 6.484 -30.268 -0.117) 30.955 8.679 ( -31.374 -7.427 0.250) 32.242 9.852 ( -7.927 -22.982 0.945) 24.329 10.443 ( -4.832 -39.074 0.651) 39.377 11.902 ( -7.253 -8.111 -0.430) 10.889 12.450 ( -0.422 -22.272 -0.538) 22.282 12.854 ( 1.438 9.718 1.427) 9.927 15.092 ( 0.935 8.811 -1.114) 8.930 15.940 ( 2.177 10.771 1.299) 11.065 16.468 ( 0.298 10.004 -1.045) 10.063 ======================= Grid point 408 (82/84) ======================= q-point: ( 0.38 0.38 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.031 ( 1.071 1.071 0.452) 1.581 2.657 ( 5.279 5.279 -0.243) 7.470 2.724 ( 3.305 3.305 -0.299) 4.684 2.979 ( 10.622 10.622 -0.066) 15.021 2.995 ( 11.610 11.610 1.517) 16.489 3.126 ( 4.682 4.682 -1.098) 6.712 3.397 ( 10.660 10.660 0.337) 15.080 3.431 ( 2.240 2.240 -1.195) 3.385 3.711 ( 10.649 10.649 0.302) 15.063 5.040 ( -11.254 -11.254 0.643) 15.928 6.156 ( 7.059 7.059 -0.139) 9.984 6.165 ( 7.246 7.246 -0.185) 10.249 6.405 ( -22.762 -22.762 0.215) 32.191 7.508 ( 1.726 1.726 -1.068) 2.664 7.676 ( -6.890 -6.890 -0.188) 9.746 9.056 ( -12.792 -12.792 1.802) 18.180 9.732 ( -13.756 -13.756 1.017) 19.480 9.979 ( -29.874 -29.874 -0.152) 42.249 12.026 ( -12.227 -12.227 -0.827) 17.311 12.086 ( -11.999 -11.999 -0.235) 16.970 12.935 ( 4.980 4.980 1.607) 7.224 15.173 ( 4.670 4.670 -1.200) 6.713 16.022 ( 6.187 6.187 0.792) 8.785 16.566 ( 3.627 3.627 -0.861) 5.202 ======================= Grid point 409 (83/84) ======================= q-point: ( 0.46 0.38 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 423 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.046 ( 0.379 1.344 0.507) 1.486 2.712 ( 1.102 6.913 -0.400) 7.012 2.789 ( 2.121 1.173 -0.265) 2.438 3.090 ( 2.003 12.463 0.473) 12.632 3.162 ( 4.404 9.093 1.180) 10.172 3.241 ( 4.635 2.467 0.254) 5.257 3.445 ( -0.197 3.233 -0.985) 3.385 3.580 ( 6.070 17.759 -0.296) 18.770 3.867 ( 4.466 11.918 -0.185) 12.729 4.855 ( -6.018 -15.183 0.681) 16.347 6.088 ( -8.525 -30.969 0.531) 32.125 6.182 ( 0.575 12.013 -0.288) 12.031 6.330 ( 4.347 1.105 -0.199) 4.490 7.610 ( -0.887 -19.271 -0.379) 19.295 7.709 ( 14.424 -3.505 -1.528) 14.922 8.662 ( -21.769 2.112 1.702) 21.937 9.568 ( -3.630 -8.898 1.847) 9.786 9.627 ( -6.550 -39.784 0.161) 40.320 11.746 ( -8.948 -7.056 -0.566) 11.409 11.997 ( -1.821 -22.186 -0.566) 22.268 12.998 ( 1.541 5.122 1.619) 5.588 15.237 ( 1.640 5.658 -1.450) 6.067 16.102 ( 1.950 5.817 0.496) 6.155 16.608 ( 0.902 4.588 -0.603) 4.715 ======================= Grid point 422 (84/84) ======================= q-point: ( 0.46 0.46 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 3.36e-05 3.36e-05 3.36e-05 3.36e-05 0.00e+00 0.00e+00 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.064 ( 0.486 0.486 0.576) 0.897 2.803 ( 1.952 1.952 -0.389) 2.787 2.811 ( 1.134 1.134 -0.354) 1.643 3.272 ( 3.817 3.817 1.270) 5.545 3.276 ( 3.224 3.224 1.875) 4.930 3.299 ( 3.010 3.010 1.205) 4.424 3.480 ( 0.533 0.533 -1.743) 1.899 3.912 ( 13.072 13.072 -0.545) 18.495 4.035 ( 4.911 4.911 -0.763) 6.988 4.570 ( -11.800 -11.800 0.838) 16.709 5.645 ( -12.189 -12.189 0.673) 17.251 6.342 ( 2.485 2.485 -0.347) 3.531 6.346 ( 2.149 2.149 -0.340) 3.058 7.351 ( -6.760 -6.760 -0.631) 9.581 7.841 ( 13.039 13.039 -2.022) 18.551 8.668 ( -5.652 -5.652 2.530) 8.384 8.891 ( -22.718 -22.718 1.293) 32.154 9.546 ( 1.052 1.052 1.244) 1.939 11.638 ( -5.767 -5.767 -0.621) 8.180 11.658 ( -7.470 -7.470 -0.537) 10.578 13.063 ( 1.576 1.576 1.709) 2.808 15.312 ( 1.929 1.929 -1.660) 3.194 16.176 ( 1.805 1.805 0.094) 2.554 16.663 ( 1.249 1.249 -0.388) 1.809 =================== 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/20280 10.0 1123.044 1123.044 441.765 0.000 -0.000 0.000 3/20280 20.0 332.137 332.137 47.965 0.000 -0.000 0.000 3/20280 30.0 205.485 205.485 22.982 0.000 -0.000 0.000 3/20280 40.0 142.350 142.350 15.372 0.000 -0.000 0.000 3/20280 50.0 104.693 104.693 11.358 0.000 -0.000 0.000 3/20280 60.0 80.648 80.648 8.832 0.000 -0.000 0.000 3/20280 70.0 64.425 64.425 7.115 0.000 -0.000 0.000 3/20280 80.0 52.992 52.992 5.892 0.000 -0.000 0.000 3/20280 90.0 44.647 44.647 4.991 0.000 -0.000 0.000 3/20280 100.0 38.375 38.375 4.308 0.000 -0.000 0.000 3/20280 110.0 33.541 33.541 3.779 0.000 -0.000 0.000 3/20280 120.0 29.732 29.732 3.359 0.000 -0.000 0.000 3/20280 130.0 26.671 26.671 3.020 0.000 -0.000 0.000 3/20280 140.0 24.168 24.168 2.741 0.000 -0.000 0.000 3/20280 150.0 22.090 22.090 2.509 0.000 -0.000 0.000 3/20280 160.0 20.340 20.340 2.314 0.000 -0.000 0.000 3/20280 170.0 18.850 18.850 2.146 0.000 -0.000 0.000 3/20280 180.0 17.565 17.565 2.002 0.000 -0.000 0.000 3/20280 190.0 16.448 16.448 1.876 0.000 -0.000 0.000 3/20280 200.0 15.467 15.467 1.766 0.000 -0.000 0.000 3/20280 210.0 14.599 14.599 1.668 0.000 -0.000 0.000 3/20280 220.0 13.826 13.826 1.581 0.000 -0.000 0.000 3/20280 230.0 13.134 13.134 1.502 0.000 -0.000 0.000 3/20280 240.0 12.509 12.509 1.431 0.000 -0.000 0.000 3/20280 250.0 11.942 11.942 1.367 0.000 -0.000 0.000 3/20280 260.0 11.427 11.427 1.309 0.000 -0.000 0.000 3/20280 270.0 10.955 10.955 1.255 0.000 -0.000 0.000 3/20280 280.0 10.522 10.522 1.206 0.000 -0.000 0.000 3/20280 290.0 10.123 10.123 1.161 0.000 -0.000 0.000 3/20280 300.0 9.753 9.753 1.119 0.000 -0.000 0.000 3/20280 310.0 9.411 9.411 1.080 0.000 -0.000 0.000 3/20280 320.0 9.092 9.092 1.043 0.000 -0.000 0.000 3/20280 330.0 8.795 8.795 1.010 0.000 -0.000 0.000 3/20280 340.0 8.517 8.517 0.978 0.000 -0.000 0.000 3/20280 350.0 8.257 8.257 0.948 0.000 -0.000 0.000 3/20280 360.0 8.012 8.012 0.920 0.000 -0.000 0.000 3/20280 370.0 7.782 7.782 0.894 0.000 -0.000 0.000 3/20280 380.0 7.565 7.565 0.869 0.000 -0.000 0.000 3/20280 390.0 7.360 7.360 0.846 0.000 -0.000 0.000 3/20280 400.0 7.166 7.166 0.824 0.000 -0.000 0.000 3/20280 410.0 6.982 6.982 0.803 0.000 -0.000 0.000 3/20280 420.0 6.808 6.808 0.783 0.000 -0.000 0.000 3/20280 430.0 6.642 6.642 0.764 0.000 -0.000 0.000 3/20280 440.0 6.484 6.484 0.746 0.000 -0.000 0.000 3/20280 450.0 6.334 6.334 0.729 0.000 -0.000 0.000 3/20280 460.0 6.190 6.190 0.712 0.000 -0.000 0.000 3/20280 470.0 6.054 6.054 0.696 0.000 -0.000 0.000 3/20280 480.0 5.923 5.923 0.681 0.000 -0.000 0.000 3/20280 490.0 5.797 5.797 0.667 0.000 -0.000 0.000 3/20280 500.0 5.677 5.677 0.653 0.000 -0.000 0.000 3/20280 510.0 5.562 5.562 0.640 0.000 -0.000 0.000 3/20280 520.0 5.452 5.452 0.627 0.000 -0.000 0.000 3/20280 530.0 5.346 5.346 0.615 0.000 -0.000 0.000 3/20280 540.0 5.244 5.244 0.604 0.000 -0.000 0.000 3/20280 550.0 5.146 5.146 0.592 0.000 -0.000 0.000 3/20280 560.0 5.051 5.051 0.582 0.000 -0.000 0.000 3/20280 570.0 4.960 4.960 0.571 0.000 -0.000 0.000 3/20280 580.0 4.873 4.873 0.561 0.000 -0.000 0.000 3/20280 590.0 4.788 4.788 0.551 0.000 -0.000 0.000 3/20280 600.0 4.706 4.706 0.542 0.000 -0.000 0.000 3/20280 610.0 4.627 4.627 0.533 0.000 -0.000 0.000 3/20280 620.0 4.551 4.551 0.524 0.000 -0.000 0.000 3/20280 630.0 4.477 4.477 0.516 0.000 -0.000 0.000 3/20280 640.0 4.406 4.406 0.507 0.000 -0.000 0.000 3/20280 650.0 4.337 4.337 0.499 0.000 -0.000 0.000 3/20280 660.0 4.269 4.269 0.492 0.000 -0.000 0.000 3/20280 670.0 4.205 4.205 0.484 0.000 -0.000 0.000 3/20280 680.0 4.142 4.142 0.477 0.000 -0.000 0.000 3/20280 690.0 4.080 4.080 0.470 0.000 -0.000 0.000 3/20280 700.0 4.021 4.021 0.463 0.000 -0.000 0.000 3/20280 710.0 3.963 3.963 0.457 0.000 -0.000 0.000 3/20280 720.0 3.907 3.907 0.450 0.000 -0.000 0.000 3/20280 730.0 3.853 3.853 0.444 0.000 -0.000 0.000 3/20280 740.0 3.800 3.800 0.438 0.000 -0.000 0.000 3/20280 750.0 3.749 3.749 0.432 0.000 -0.000 0.000 3/20280 760.0 3.699 3.699 0.426 0.000 -0.000 0.000 3/20280 770.0 3.650 3.650 0.421 0.000 -0.000 0.000 3/20280 780.0 3.602 3.602 0.415 0.000 -0.000 0.000 3/20280 790.0 3.556 3.556 0.410 0.000 -0.000 0.000 3/20280 800.0 3.511 3.511 0.405 0.000 -0.000 0.000 3/20280 810.0 3.467 3.467 0.400 0.000 -0.000 0.000 3/20280 820.0 3.424 3.424 0.395 0.000 -0.000 0.000 3/20280 830.0 3.383 3.383 0.390 0.000 -0.000 0.000 3/20280 840.0 3.342 3.342 0.385 0.000 -0.000 0.000 3/20280 850.0 3.302 3.302 0.381 0.000 -0.000 0.000 3/20280 860.0 3.263 3.263 0.376 0.000 -0.000 0.000 3/20280 870.0 3.225 3.225 0.372 0.000 -0.000 0.000 3/20280 880.0 3.188 3.188 0.367 0.000 -0.000 0.000 3/20280 890.0 3.152 3.152 0.363 0.000 -0.000 0.000 3/20280 900.0 3.117 3.117 0.359 0.000 -0.000 0.000 3/20280 910.0 3.082 3.082 0.355 0.000 -0.000 0.000 3/20280 920.0 3.048 3.048 0.351 0.000 -0.000 0.000 3/20280 930.0 3.015 3.015 0.348 0.000 -0.000 0.000 3/20280 940.0 2.983 2.983 0.344 0.000 -0.000 0.000 3/20280 950.0 2.951 2.951 0.340 0.000 -0.000 0.000 3/20280 960.0 2.920 2.920 0.337 0.000 -0.000 0.000 3/20280 970.0 2.890 2.890 0.333 0.000 -0.000 0.000 3/20280 980.0 2.860 2.860 0.330 0.000 -0.000 0.000 3/20280 990.0 2.831 2.831 0.326 0.000 -0.000 0.000 3/20280 1000.0 2.802 2.802 0.323 0.000 -0.000 0.000 3/20280 Thermal conductivity related properties were written into "kappa-m13135.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 10:51:12]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|