----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 18:56:04]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P3m1 (156) Number of symmetry operations in supercell: 96 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.221951047522676 0.000000000000000 0.000000000000000 b -2.110975523761338 3.656316860688959 0.000000000000000 c 0.000000000000000 0.000000000000000 20.549632639999999 Atomic positions (fractional): *1 Cd 0.00000000000000 0.00000000000000 0.08329303494589 112.411 *2 Cd 0.33333333333333 0.66666666666667 0.41671967309688 112.411 *3 Cd 0.33333333333333 0.66666666666667 0.74995362118498 112.411 *4 I 0.00000000000000 0.00000000000000 0.66702061783345 126.904 *5 I 0.66666666666667 0.33333333333333 0.83292910146147 126.904 *6 I 0.66666666666667 0.33333333333333 0.16628926553231 126.904 *7 I 0.00000000000000 0.00000000000000 0.33375366569986 126.904 *8 I 0.33333333333333 0.66666666666667 0.00035804939328 126.904 *9 I 0.66666666666667 0.33333333333333 0.49968299085190 126.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.221951047522676 0.000000000000000 0.000000000000000 b -2.110975523761338 3.656316860688959 0.000000000000000 c 0.000000000000000 0.000000000000000 20.549632639999999 Atomic positions (fractional): *1 Cd 0.00000000000000 0.00000000000000 0.08329303494589 112.411 > 1 *2 Cd 0.33333333333333 0.66666666666667 0.41671967309688 112.411 > 2 *3 Cd 0.33333333333333 0.66666666666667 0.74995362118498 112.411 > 3 *4 I 0.00000000000000 0.00000000000000 0.66702061783345 126.904 > 4 *5 I 0.66666666666667 0.33333333333333 0.83292910146147 126.904 > 5 *6 I 0.66666666666667 0.33333333333333 0.16628926553231 126.904 > 6 *7 I 0.00000000000000 0.00000000000000 0.33375366569986 126.904 > 7 *8 I 0.33333333333333 0.66666666666667 0.00035804939328 126.904 > 8 *9 I 0.66666666666667 0.33333333333333 0.49968299085190 126.904 > 9 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.887804190090705 0.000000000000000 0.000000000000000 b -8.443902095045352 14.625267442755836 0.000000000000000 c 0.000000000000000 0.000000000000000 20.549632639999999 Atomic positions (fractional): *1 Cd 0.00000000000000 0.00000000000000 0.08329303494589 112.411 > 1 2 Cd 0.25000000000000 0.00000000000000 0.08329303494589 112.411 > 1 3 Cd 0.50000000000000 0.00000000000000 0.08329303494589 112.411 > 1 4 Cd 0.75000000000000 0.00000000000000 0.08329303494589 112.411 > 1 5 Cd 0.00000000000000 0.25000000000000 0.08329303494589 112.411 > 1 6 Cd 0.25000000000000 0.25000000000000 0.08329303494589 112.411 > 1 7 Cd 0.50000000000000 0.25000000000000 0.08329303494589 112.411 > 1 8 Cd 0.75000000000000 0.25000000000000 0.08329303494589 112.411 > 1 9 Cd 0.00000000000000 0.50000000000000 0.08329303494589 112.411 > 1 10 Cd 0.25000000000000 0.50000000000000 0.08329303494589 112.411 > 1 11 Cd 0.50000000000000 0.50000000000000 0.08329303494589 112.411 > 1 12 Cd 0.75000000000000 0.50000000000000 0.08329303494589 112.411 > 1 13 Cd 0.00000000000000 0.75000000000000 0.08329303494589 112.411 > 1 14 Cd 0.25000000000000 0.75000000000000 0.08329303494589 112.411 > 1 15 Cd 0.50000000000000 0.75000000000000 0.08329303494589 112.411 > 1 16 Cd 0.75000000000000 0.75000000000000 0.08329303494589 112.411 > 1 *17 Cd 0.08333333333333 0.16666666666667 0.41671967309688 112.411 > 2 18 Cd 0.33333333333333 0.16666666666667 0.41671967309688 112.411 > 2 19 Cd 0.58333333333333 0.16666666666667 0.41671967309688 112.411 > 2 20 Cd 0.83333333333333 0.16666666666667 0.41671967309688 112.411 > 2 21 Cd 0.08333333333333 0.41666666666667 0.41671967309688 112.411 > 2 22 Cd 0.33333333333333 0.41666666666667 0.41671967309688 112.411 > 2 23 Cd 0.58333333333333 0.41666666666667 0.41671967309688 112.411 > 2 24 Cd 0.83333333333333 0.41666666666667 0.41671967309688 112.411 > 2 25 Cd 0.08333333333333 0.66666666666667 0.41671967309688 112.411 > 2 26 Cd 0.33333333333333 0.66666666666667 0.41671967309688 112.411 > 2 27 Cd 0.58333333333333 0.66666666666667 0.41671967309688 112.411 > 2 28 Cd 0.83333333333333 0.66666666666667 0.41671967309688 112.411 > 2 29 Cd 0.08333333333333 0.91666666666667 0.41671967309688 112.411 > 2 30 Cd 0.33333333333333 0.91666666666667 0.41671967309688 112.411 > 2 31 Cd 0.58333333333333 0.91666666666667 0.41671967309688 112.411 > 2 32 Cd 0.83333333333333 0.91666666666667 0.41671967309688 112.411 > 2 *33 Cd 0.08333333333333 0.16666666666667 0.74995362118498 112.411 > 3 34 Cd 0.33333333333333 0.16666666666667 0.74995362118498 112.411 > 3 35 Cd 0.58333333333333 0.16666666666667 0.74995362118498 112.411 > 3 36 Cd 0.83333333333333 0.16666666666667 0.74995362118498 112.411 > 3 37 Cd 0.08333333333333 0.41666666666667 0.74995362118498 112.411 > 3 38 Cd 0.33333333333333 0.41666666666667 0.74995362118498 112.411 > 3 39 Cd 0.58333333333333 0.41666666666667 0.74995362118498 112.411 > 3 40 Cd 0.83333333333333 0.41666666666667 0.74995362118498 112.411 > 3 41 Cd 0.08333333333333 0.66666666666667 0.74995362118498 112.411 > 3 42 Cd 0.33333333333333 0.66666666666667 0.74995362118498 112.411 > 3 43 Cd 0.58333333333333 0.66666666666667 0.74995362118498 112.411 > 3 44 Cd 0.83333333333333 0.66666666666667 0.74995362118498 112.411 > 3 45 Cd 0.08333333333333 0.91666666666667 0.74995362118498 112.411 > 3 46 Cd 0.33333333333333 0.91666666666667 0.74995362118498 112.411 > 3 47 Cd 0.58333333333333 0.91666666666667 0.74995362118498 112.411 > 3 48 Cd 0.83333333333333 0.91666666666667 0.74995362118498 112.411 > 3 *49 I 0.00000000000000 0.00000000000000 0.66702061783345 126.904 > 4 50 I 0.25000000000000 0.00000000000000 0.66702061783345 126.904 > 4 51 I 0.50000000000000 0.00000000000000 0.66702061783345 126.904 > 4 52 I 0.75000000000000 0.00000000000000 0.66702061783345 126.904 > 4 53 I 0.00000000000000 0.25000000000000 0.66702061783345 126.904 > 4 54 I 0.25000000000000 0.25000000000000 0.66702061783345 126.904 > 4 55 I 0.50000000000000 0.25000000000000 0.66702061783345 126.904 > 4 56 I 0.75000000000000 0.25000000000000 0.66702061783345 126.904 > 4 57 I 0.00000000000000 0.50000000000000 0.66702061783345 126.904 > 4 58 I 0.25000000000000 0.50000000000000 0.66702061783345 126.904 > 4 59 I 0.50000000000000 0.50000000000000 0.66702061783345 126.904 > 4 60 I 0.75000000000000 0.50000000000000 0.66702061783345 126.904 > 4 61 I 0.00000000000000 0.75000000000000 0.66702061783345 126.904 > 4 62 I 0.25000000000000 0.75000000000000 0.66702061783345 126.904 > 4 63 I 0.50000000000000 0.75000000000000 0.66702061783345 126.904 > 4 64 I 0.75000000000000 0.75000000000000 0.66702061783345 126.904 > 4 *65 I 0.16666666666667 0.08333333333333 0.83292910146147 126.904 > 5 66 I 0.41666666666667 0.08333333333333 0.83292910146147 126.904 > 5 67 I 0.66666666666667 0.08333333333333 0.83292910146147 126.904 > 5 68 I 0.91666666666667 0.08333333333333 0.83292910146147 126.904 > 5 69 I 0.16666666666667 0.33333333333333 0.83292910146147 126.904 > 5 70 I 0.41666666666667 0.33333333333333 0.83292910146147 126.904 > 5 71 I 0.66666666666667 0.33333333333333 0.83292910146147 126.904 > 5 72 I 0.91666666666667 0.33333333333333 0.83292910146147 126.904 > 5 73 I 0.16666666666667 0.58333333333333 0.83292910146147 126.904 > 5 74 I 0.41666666666667 0.58333333333333 0.83292910146147 126.904 > 5 75 I 0.66666666666667 0.58333333333333 0.83292910146147 126.904 > 5 76 I 0.91666666666667 0.58333333333333 0.83292910146147 126.904 > 5 77 I 0.16666666666667 0.83333333333333 0.83292910146147 126.904 > 5 78 I 0.41666666666667 0.83333333333333 0.83292910146147 126.904 > 5 79 I 0.66666666666667 0.83333333333333 0.83292910146147 126.904 > 5 80 I 0.91666666666667 0.83333333333333 0.83292910146147 126.904 > 5 *81 I 0.16666666666667 0.08333333333333 0.16628926553231 126.904 > 6 82 I 0.41666666666667 0.08333333333333 0.16628926553231 126.904 > 6 83 I 0.66666666666667 0.08333333333333 0.16628926553231 126.904 > 6 84 I 0.91666666666667 0.08333333333333 0.16628926553231 126.904 > 6 85 I 0.16666666666667 0.33333333333333 0.16628926553231 126.904 > 6 86 I 0.41666666666667 0.33333333333333 0.16628926553231 126.904 > 6 87 I 0.66666666666667 0.33333333333333 0.16628926553231 126.904 > 6 88 I 0.91666666666667 0.33333333333333 0.16628926553231 126.904 > 6 89 I 0.16666666666667 0.58333333333333 0.16628926553231 126.904 > 6 90 I 0.41666666666667 0.58333333333333 0.16628926553231 126.904 > 6 91 I 0.66666666666667 0.58333333333333 0.16628926553231 126.904 > 6 92 I 0.91666666666667 0.58333333333333 0.16628926553231 126.904 > 6 93 I 0.16666666666667 0.83333333333333 0.16628926553231 126.904 > 6 94 I 0.41666666666667 0.83333333333333 0.16628926553231 126.904 > 6 95 I 0.66666666666667 0.83333333333333 0.16628926553231 126.904 > 6 96 I 0.91666666666667 0.83333333333333 0.16628926553231 126.904 > 6 *97 I 0.00000000000000 0.00000000000000 0.33375366569986 126.904 > 7 98 I 0.25000000000000 0.00000000000000 0.33375366569986 126.904 > 7 99 I 0.50000000000000 0.00000000000000 0.33375366569986 126.904 > 7 100 I 0.75000000000000 0.00000000000000 0.33375366569986 126.904 > 7 101 I 0.00000000000000 0.25000000000000 0.33375366569986 126.904 > 7 102 I 0.25000000000000 0.25000000000000 0.33375366569986 126.904 > 7 103 I 0.50000000000000 0.25000000000000 0.33375366569986 126.904 > 7 104 I 0.75000000000000 0.25000000000000 0.33375366569986 126.904 > 7 105 I 0.00000000000000 0.50000000000000 0.33375366569986 126.904 > 7 106 I 0.25000000000000 0.50000000000000 0.33375366569986 126.904 > 7 107 I 0.50000000000000 0.50000000000000 0.33375366569986 126.904 > 7 108 I 0.75000000000000 0.50000000000000 0.33375366569986 126.904 > 7 109 I 0.00000000000000 0.75000000000000 0.33375366569986 126.904 > 7 110 I 0.25000000000000 0.75000000000000 0.33375366569986 126.904 > 7 111 I 0.50000000000000 0.75000000000000 0.33375366569986 126.904 > 7 112 I 0.75000000000000 0.75000000000000 0.33375366569986 126.904 > 7 *113 I 0.08333333333333 0.16666666666667 0.00035804939328 126.904 > 8 114 I 0.33333333333333 0.16666666666667 0.00035804939328 126.904 > 8 115 I 0.58333333333333 0.16666666666667 0.00035804939328 126.904 > 8 116 I 0.83333333333333 0.16666666666667 0.00035804939328 126.904 > 8 117 I 0.08333333333333 0.41666666666667 0.00035804939328 126.904 > 8 118 I 0.33333333333333 0.41666666666667 0.00035804939328 126.904 > 8 119 I 0.58333333333333 0.41666666666667 0.00035804939328 126.904 > 8 120 I 0.83333333333333 0.41666666666667 0.00035804939328 126.904 > 8 121 I 0.08333333333333 0.66666666666667 0.00035804939328 126.904 > 8 122 I 0.33333333333333 0.66666666666667 0.00035804939328 126.904 > 8 123 I 0.58333333333333 0.66666666666667 0.00035804939328 126.904 > 8 124 I 0.83333333333333 0.66666666666667 0.00035804939328 126.904 > 8 125 I 0.08333333333333 0.91666666666667 0.00035804939328 126.904 > 8 126 I 0.33333333333333 0.91666666666667 0.00035804939328 126.904 > 8 127 I 0.58333333333333 0.91666666666667 0.00035804939328 126.904 > 8 128 I 0.83333333333333 0.91666666666667 0.00035804939328 126.904 > 8 *129 I 0.16666666666667 0.08333333333333 0.49968299085190 126.904 > 9 130 I 0.41666666666667 0.08333333333333 0.49968299085190 126.904 > 9 131 I 0.66666666666667 0.08333333333333 0.49968299085190 126.904 > 9 132 I 0.91666666666667 0.08333333333333 0.49968299085190 126.904 > 9 133 I 0.16666666666667 0.33333333333333 0.49968299085190 126.904 > 9 134 I 0.41666666666667 0.33333333333333 0.49968299085190 126.904 > 9 135 I 0.66666666666667 0.33333333333333 0.49968299085190 126.904 > 9 136 I 0.91666666666667 0.33333333333333 0.49968299085190 126.904 > 9 137 I 0.16666666666667 0.58333333333333 0.49968299085190 126.904 > 9 138 I 0.41666666666667 0.58333333333333 0.49968299085190 126.904 > 9 139 I 0.66666666666667 0.58333333333333 0.49968299085190 126.904 > 9 140 I 0.91666666666667 0.58333333333333 0.49968299085190 126.904 > 9 141 I 0.16666666666667 0.83333333333333 0.49968299085190 126.904 > 9 142 I 0.41666666666667 0.83333333333333 0.49968299085190 126.904 > 9 143 I 0.66666666666667 0.83333333333333 0.49968299085190 126.904 > 9 144 I 0.91666666666667 0.83333333333333 0.49968299085190 126.904 > 9 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 6.0550211 -0.0000000 0.0000000 0.0000000 6.0550211 0.0000000 0.0000000 0.0000000 4.8840961 -------------------------- Born effective charges -------------------------- 1 Cd 3.0193209 -0.0000000 0.0000000 0.0000000 3.0193209 0.0000000 0.0000000 0.0000000 1.6871870 2 Cd 3.0367361 -0.0000000 0.0000000 0.0000000 3.0367361 0.0000000 0.0000000 0.0000000 1.6919068 3 Cd 3.0236693 -0.0000000 0.0000000 0.0000000 3.0236693 0.0000000 0.0000000 0.0000000 1.6833891 4 I -1.5328396 0.0000000 0.0000000 0.0000000 -1.5328396 0.0000000 0.0000000 0.0000000 -0.8455897 5 I -1.5012836 0.0000000 0.0000000 0.0000000 -1.5012836 0.0000000 0.0000000 0.0000000 -0.8456941 6 I -1.5027788 0.0000000 0.0000000 0.0000000 -1.5027788 0.0000000 0.0000000 0.0000000 -0.8497911 7 I -1.5031141 0.0000000 0.0000000 0.0000000 -1.5031141 0.0000000 0.0000000 0.0000000 -0.8375722 8 I -1.5161627 0.0000000 0.0000000 0.0000000 -1.5161627 0.0000000 0.0000000 0.0000000 -0.8360942 9 I -1.5235474 0.0000000 0.0000000 0.0000000 -1.5235474 0.0000000 0.0000000 0.0000000 -0.8477416 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 432/432 Permutation basis: 6426/6426 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 1089 Number of blocks in projector: 837 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 468 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 405 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 216 Use standard eigh solver. Tree of FC basis block matrices: - (1089, 1071), data: False |-- (216, 207), data: True |-- (405, 405), data: True |-- (468, 459), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.013 Solver_block: 100 / 240 - Time: 1.452 Solver_block: 200 / 240 - Time: 1.472 Solver_block: 240 / 240 - Time: 0.679 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 3.634 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 432/432 Permutation basis: 6426/6426 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 1089 Number of blocks in projector: 837 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 468 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 405 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 216 Use standard eigh solver. Tree of FC basis block matrices: - (1089, 1071), data: False |-- (216, 207), data: True |-- (405, 405), data: True |-- (468, 459), data: True Max drift after symmetrization by symfc projector: 0.00000000 (xx) 0.00000000 (xx) 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 18:56:20]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 18:56:21]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P3m1 (156) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.221951047522676 0.000000000000000 0.000000000000000 b -2.110975523761338 3.656316860688959 0.000000000000000 c 0.000000000000000 0.000000000000000 20.549632639999999 Atomic positions (fractional): 1 Cd 0.00000000000000 0.00000000000000 0.08329303494589 112.411 2 Cd 0.33333333333333 0.66666666666667 0.41671967309688 112.411 3 Cd 0.33333333333333 0.66666666666667 0.74995362118498 112.411 4 I 0.00000000000000 0.00000000000000 0.66702061783345 126.904 5 I 0.66666666666667 0.33333333333333 0.83292910146147 126.904 6 I 0.66666666666667 0.33333333333333 0.16628926553231 126.904 7 I 0.00000000000000 0.00000000000000 0.33375366569986 126.904 8 I 0.33333333333333 0.66666666666667 0.00035804939328 126.904 9 I 0.66666666666667 0.33333333333333 0.49968299085190 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.887804190090705 0.000000000000000 0.000000000000000 b -8.443902095045352 14.625267442755836 0.000000000000000 c 0.000000000000000 0.000000000000000 20.549632639999999 Atomic positions (fractional): 1 Cd 0.00000000000000 0.00000000000000 0.08329303494589 112.411 > 1 2 Cd 0.25000000000000 0.00000000000000 0.08329303494589 112.411 > 1 3 Cd 0.50000000000000 0.00000000000000 0.08329303494589 112.411 > 1 4 Cd 0.75000000000000 0.00000000000000 0.08329303494589 112.411 > 1 5 Cd 0.00000000000000 0.25000000000000 0.08329303494589 112.411 > 1 6 Cd 0.25000000000000 0.25000000000000 0.08329303494589 112.411 > 1 7 Cd 0.50000000000000 0.25000000000000 0.08329303494589 112.411 > 1 8 Cd 0.75000000000000 0.25000000000000 0.08329303494589 112.411 > 1 9 Cd 0.00000000000000 0.50000000000000 0.08329303494589 112.411 > 1 10 Cd 0.25000000000000 0.50000000000000 0.08329303494589 112.411 > 1 11 Cd 0.50000000000000 0.50000000000000 0.08329303494589 112.411 > 1 12 Cd 0.75000000000000 0.50000000000000 0.08329303494589 112.411 > 1 13 Cd 0.00000000000000 0.75000000000000 0.08329303494589 112.411 > 1 14 Cd 0.25000000000000 0.75000000000000 0.08329303494589 112.411 > 1 15 Cd 0.50000000000000 0.75000000000000 0.08329303494589 112.411 > 1 16 Cd 0.75000000000000 0.75000000000000 0.08329303494589 112.411 > 1 17 Cd 0.08333333333333 0.16666666666667 0.41671967309688 112.411 > 17 18 Cd 0.33333333333333 0.16666666666667 0.41671967309688 112.411 > 17 19 Cd 0.58333333333333 0.16666666666667 0.41671967309688 112.411 > 17 20 Cd 0.83333333333333 0.16666666666667 0.41671967309688 112.411 > 17 21 Cd 0.08333333333333 0.41666666666667 0.41671967309688 112.411 > 17 22 Cd 0.33333333333333 0.41666666666667 0.41671967309688 112.411 > 17 23 Cd 0.58333333333333 0.41666666666667 0.41671967309688 112.411 > 17 24 Cd 0.83333333333333 0.41666666666667 0.41671967309688 112.411 > 17 25 Cd 0.08333333333333 0.66666666666667 0.41671967309688 112.411 > 17 26 Cd 0.33333333333333 0.66666666666667 0.41671967309688 112.411 > 17 27 Cd 0.58333333333333 0.66666666666667 0.41671967309688 112.411 > 17 28 Cd 0.83333333333333 0.66666666666667 0.41671967309688 112.411 > 17 29 Cd 0.08333333333333 0.91666666666667 0.41671967309688 112.411 > 17 30 Cd 0.33333333333333 0.91666666666667 0.41671967309688 112.411 > 17 31 Cd 0.58333333333333 0.91666666666667 0.41671967309688 112.411 > 17 32 Cd 0.83333333333333 0.91666666666667 0.41671967309688 112.411 > 17 33 Cd 0.08333333333333 0.16666666666667 0.74995362118498 112.411 > 33 34 Cd 0.33333333333333 0.16666666666667 0.74995362118498 112.411 > 33 35 Cd 0.58333333333333 0.16666666666667 0.74995362118498 112.411 > 33 36 Cd 0.83333333333333 0.16666666666667 0.74995362118498 112.411 > 33 37 Cd 0.08333333333333 0.41666666666667 0.74995362118498 112.411 > 33 38 Cd 0.33333333333333 0.41666666666667 0.74995362118498 112.411 > 33 39 Cd 0.58333333333333 0.41666666666667 0.74995362118498 112.411 > 33 40 Cd 0.83333333333333 0.41666666666667 0.74995362118498 112.411 > 33 41 Cd 0.08333333333333 0.66666666666667 0.74995362118498 112.411 > 33 42 Cd 0.33333333333333 0.66666666666667 0.74995362118498 112.411 > 33 43 Cd 0.58333333333333 0.66666666666667 0.74995362118498 112.411 > 33 44 Cd 0.83333333333333 0.66666666666667 0.74995362118498 112.411 > 33 45 Cd 0.08333333333333 0.91666666666667 0.74995362118498 112.411 > 33 46 Cd 0.33333333333333 0.91666666666667 0.74995362118498 112.411 > 33 47 Cd 0.58333333333333 0.91666666666667 0.74995362118498 112.411 > 33 48 Cd 0.83333333333333 0.91666666666667 0.74995362118498 112.411 > 33 49 I 0.00000000000000 0.00000000000000 0.66702061783345 126.904 > 49 50 I 0.25000000000000 0.00000000000000 0.66702061783345 126.904 > 49 51 I 0.50000000000000 0.00000000000000 0.66702061783345 126.904 > 49 52 I 0.75000000000000 0.00000000000000 0.66702061783345 126.904 > 49 53 I 0.00000000000000 0.25000000000000 0.66702061783345 126.904 > 49 54 I 0.25000000000000 0.25000000000000 0.66702061783345 126.904 > 49 55 I 0.50000000000000 0.25000000000000 0.66702061783345 126.904 > 49 56 I 0.75000000000000 0.25000000000000 0.66702061783345 126.904 > 49 57 I 0.00000000000000 0.50000000000000 0.66702061783345 126.904 > 49 58 I 0.25000000000000 0.50000000000000 0.66702061783345 126.904 > 49 59 I 0.50000000000000 0.50000000000000 0.66702061783345 126.904 > 49 60 I 0.75000000000000 0.50000000000000 0.66702061783345 126.904 > 49 61 I 0.00000000000000 0.75000000000000 0.66702061783345 126.904 > 49 62 I 0.25000000000000 0.75000000000000 0.66702061783345 126.904 > 49 63 I 0.50000000000000 0.75000000000000 0.66702061783345 126.904 > 49 64 I 0.75000000000000 0.75000000000000 0.66702061783345 126.904 > 49 65 I 0.16666666666667 0.08333333333333 0.83292910146147 126.904 > 65 66 I 0.41666666666667 0.08333333333333 0.83292910146147 126.904 > 65 67 I 0.66666666666667 0.08333333333333 0.83292910146147 126.904 > 65 68 I 0.91666666666667 0.08333333333333 0.83292910146147 126.904 > 65 69 I 0.16666666666667 0.33333333333333 0.83292910146147 126.904 > 65 70 I 0.41666666666667 0.33333333333333 0.83292910146147 126.904 > 65 71 I 0.66666666666667 0.33333333333333 0.83292910146147 126.904 > 65 72 I 0.91666666666667 0.33333333333333 0.83292910146147 126.904 > 65 73 I 0.16666666666667 0.58333333333333 0.83292910146147 126.904 > 65 74 I 0.41666666666667 0.58333333333333 0.83292910146147 126.904 > 65 75 I 0.66666666666667 0.58333333333333 0.83292910146147 126.904 > 65 76 I 0.91666666666667 0.58333333333333 0.83292910146147 126.904 > 65 77 I 0.16666666666667 0.83333333333333 0.83292910146147 126.904 > 65 78 I 0.41666666666667 0.83333333333333 0.83292910146147 126.904 > 65 79 I 0.66666666666667 0.83333333333333 0.83292910146147 126.904 > 65 80 I 0.91666666666667 0.83333333333333 0.83292910146147 126.904 > 65 81 I 0.16666666666667 0.08333333333333 0.16628926553231 126.904 > 81 82 I 0.41666666666667 0.08333333333333 0.16628926553231 126.904 > 81 83 I 0.66666666666667 0.08333333333333 0.16628926553231 126.904 > 81 84 I 0.91666666666667 0.08333333333333 0.16628926553231 126.904 > 81 85 I 0.16666666666667 0.33333333333333 0.16628926553231 126.904 > 81 86 I 0.41666666666667 0.33333333333333 0.16628926553231 126.904 > 81 87 I 0.66666666666667 0.33333333333333 0.16628926553231 126.904 > 81 88 I 0.91666666666667 0.33333333333333 0.16628926553231 126.904 > 81 89 I 0.16666666666667 0.58333333333333 0.16628926553231 126.904 > 81 90 I 0.41666666666667 0.58333333333333 0.16628926553231 126.904 > 81 91 I 0.66666666666667 0.58333333333333 0.16628926553231 126.904 > 81 92 I 0.91666666666667 0.58333333333333 0.16628926553231 126.904 > 81 93 I 0.16666666666667 0.83333333333333 0.16628926553231 126.904 > 81 94 I 0.41666666666667 0.83333333333333 0.16628926553231 126.904 > 81 95 I 0.66666666666667 0.83333333333333 0.16628926553231 126.904 > 81 96 I 0.91666666666667 0.83333333333333 0.16628926553231 126.904 > 81 97 I 0.00000000000000 0.00000000000000 0.33375366569986 126.904 > 97 98 I 0.25000000000000 0.00000000000000 0.33375366569986 126.904 > 97 99 I 0.50000000000000 0.00000000000000 0.33375366569986 126.904 > 97 100 I 0.75000000000000 0.00000000000000 0.33375366569986 126.904 > 97 101 I 0.00000000000000 0.25000000000000 0.33375366569986 126.904 > 97 102 I 0.25000000000000 0.25000000000000 0.33375366569986 126.904 > 97 103 I 0.50000000000000 0.25000000000000 0.33375366569986 126.904 > 97 104 I 0.75000000000000 0.25000000000000 0.33375366569986 126.904 > 97 105 I 0.00000000000000 0.50000000000000 0.33375366569986 126.904 > 97 106 I 0.25000000000000 0.50000000000000 0.33375366569986 126.904 > 97 107 I 0.50000000000000 0.50000000000000 0.33375366569986 126.904 > 97 108 I 0.75000000000000 0.50000000000000 0.33375366569986 126.904 > 97 109 I 0.00000000000000 0.75000000000000 0.33375366569986 126.904 > 97 110 I 0.25000000000000 0.75000000000000 0.33375366569986 126.904 > 97 111 I 0.50000000000000 0.75000000000000 0.33375366569986 126.904 > 97 112 I 0.75000000000000 0.75000000000000 0.33375366569986 126.904 > 97 113 I 0.08333333333333 0.16666666666667 0.00035804939328 126.904 > 113 114 I 0.33333333333333 0.16666666666667 0.00035804939328 126.904 > 113 115 I 0.58333333333333 0.16666666666667 0.00035804939328 126.904 > 113 116 I 0.83333333333333 0.16666666666667 0.00035804939328 126.904 > 113 117 I 0.08333333333333 0.41666666666667 0.00035804939328 126.904 > 113 118 I 0.33333333333333 0.41666666666667 0.00035804939328 126.904 > 113 119 I 0.58333333333333 0.41666666666667 0.00035804939328 126.904 > 113 120 I 0.83333333333333 0.41666666666667 0.00035804939328 126.904 > 113 121 I 0.08333333333333 0.66666666666667 0.00035804939328 126.904 > 113 122 I 0.33333333333333 0.66666666666667 0.00035804939328 126.904 > 113 123 I 0.58333333333333 0.66666666666667 0.00035804939328 126.904 > 113 124 I 0.83333333333333 0.66666666666667 0.00035804939328 126.904 > 113 125 I 0.08333333333333 0.91666666666667 0.00035804939328 126.904 > 113 126 I 0.33333333333333 0.91666666666667 0.00035804939328 126.904 > 113 127 I 0.58333333333333 0.91666666666667 0.00035804939328 126.904 > 113 128 I 0.83333333333333 0.91666666666667 0.00035804939328 126.904 > 113 129 I 0.16666666666667 0.08333333333333 0.49968299085190 126.904 > 129 130 I 0.41666666666667 0.08333333333333 0.49968299085190 126.904 > 129 131 I 0.66666666666667 0.08333333333333 0.49968299085190 126.904 > 129 132 I 0.91666666666667 0.08333333333333 0.49968299085190 126.904 > 129 133 I 0.16666666666667 0.33333333333333 0.49968299085190 126.904 > 129 134 I 0.41666666666667 0.33333333333333 0.49968299085190 126.904 > 129 135 I 0.66666666666667 0.33333333333333 0.49968299085190 126.904 > 129 136 I 0.91666666666667 0.33333333333333 0.49968299085190 126.904 > 129 137 I 0.16666666666667 0.58333333333333 0.49968299085190 126.904 > 129 138 I 0.41666666666667 0.58333333333333 0.49968299085190 126.904 > 129 139 I 0.66666666666667 0.58333333333333 0.49968299085190 126.904 > 129 140 I 0.91666666666667 0.58333333333333 0.49968299085190 126.904 > 129 141 I 0.16666666666667 0.83333333333333 0.49968299085190 126.904 > 129 142 I 0.41666666666667 0.83333333333333 0.49968299085190 126.904 > 129 143 I 0.66666666666667 0.83333333333333 0.49968299085190 126.904 > 129 144 I 0.91666666666667 0.83333333333333 0.49968299085190 126.904 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 6.0550211 -0.0000000 0.0000000 0.0000000 6.0550211 0.0000000 0.0000000 0.0000000 4.8840961 -------------------------- Born effective charges -------------------------- 1 Cd 3.0193209 -0.0000000 0.0000000 0.0000000 3.0193209 0.0000000 0.0000000 0.0000000 1.6871870 2 Cd 3.0367361 -0.0000000 0.0000000 0.0000000 3.0367361 0.0000000 0.0000000 0.0000000 1.6919068 3 Cd 3.0236693 -0.0000000 0.0000000 0.0000000 3.0236693 0.0000000 0.0000000 0.0000000 1.6833891 4 I -1.5328396 0.0000000 0.0000000 0.0000000 -1.5328396 0.0000000 0.0000000 0.0000000 -0.8455897 5 I -1.5012836 0.0000000 0.0000000 0.0000000 -1.5012836 0.0000000 0.0000000 0.0000000 -0.8456941 6 I -1.5027788 0.0000000 0.0000000 0.0000000 -1.5027788 0.0000000 0.0000000 0.0000000 -0.8497911 7 I -1.5031141 0.0000000 0.0000000 0.0000000 -1.5031141 0.0000000 0.0000000 0.0000000 -0.8375722 8 I -1.5161627 0.0000000 0.0000000 0.0000000 -1.5161627 0.0000000 0.0000000 0.0000000 -0.8360942 9 I -1.5235474 0.0000000 0.0000000 0.0000000 -1.5235474 0.0000000 0.0000000 0.0000000 -0.8477416 ---------------------------------------------------------------------------- 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[ 17, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 49, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 81, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 97, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 113, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 129, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000009 (zzz) -0.00000009 (zzz) -0.00000009 (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 18:56:45]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 18:56:46]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P3m1 (156) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.221951047522676 0.000000000000000 0.000000000000000 b -2.110975523761338 3.656316860688959 0.000000000000000 c 0.000000000000000 0.000000000000000 20.549632639999999 Atomic positions (fractional): 1 Cd 0.00000000000000 0.00000000000000 0.08329303494589 112.411 2 Cd 0.33333333333333 0.66666666666667 0.41671967309688 112.411 3 Cd 0.33333333333333 0.66666666666667 0.74995362118498 112.411 4 I 0.00000000000000 0.00000000000000 0.66702061783345 126.904 5 I 0.66666666666667 0.33333333333333 0.83292910146147 126.904 6 I 0.66666666666667 0.33333333333333 0.16628926553231 126.904 7 I 0.00000000000000 0.00000000000000 0.33375366569986 126.904 8 I 0.33333333333333 0.66666666666667 0.00035804939328 126.904 9 I 0.66666666666667 0.33333333333333 0.49968299085190 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.887804190090705 0.000000000000000 0.000000000000000 b -8.443902095045352 14.625267442755836 0.000000000000000 c 0.000000000000000 0.000000000000000 20.549632639999999 Atomic positions (fractional): 1 Cd 0.00000000000000 0.00000000000000 0.08329303494589 112.411 > 1 2 Cd 0.25000000000000 0.00000000000000 0.08329303494589 112.411 > 1 3 Cd 0.50000000000000 0.00000000000000 0.08329303494589 112.411 > 1 4 Cd 0.75000000000000 0.00000000000000 0.08329303494589 112.411 > 1 5 Cd 0.00000000000000 0.25000000000000 0.08329303494589 112.411 > 1 6 Cd 0.25000000000000 0.25000000000000 0.08329303494589 112.411 > 1 7 Cd 0.50000000000000 0.25000000000000 0.08329303494589 112.411 > 1 8 Cd 0.75000000000000 0.25000000000000 0.08329303494589 112.411 > 1 9 Cd 0.00000000000000 0.50000000000000 0.08329303494589 112.411 > 1 10 Cd 0.25000000000000 0.50000000000000 0.08329303494589 112.411 > 1 11 Cd 0.50000000000000 0.50000000000000 0.08329303494589 112.411 > 1 12 Cd 0.75000000000000 0.50000000000000 0.08329303494589 112.411 > 1 13 Cd 0.00000000000000 0.75000000000000 0.08329303494589 112.411 > 1 14 Cd 0.25000000000000 0.75000000000000 0.08329303494589 112.411 > 1 15 Cd 0.50000000000000 0.75000000000000 0.08329303494589 112.411 > 1 16 Cd 0.75000000000000 0.75000000000000 0.08329303494589 112.411 > 1 17 Cd 0.08333333333333 0.16666666666667 0.41671967309688 112.411 > 17 18 Cd 0.33333333333333 0.16666666666667 0.41671967309688 112.411 > 17 19 Cd 0.58333333333333 0.16666666666667 0.41671967309688 112.411 > 17 20 Cd 0.83333333333333 0.16666666666667 0.41671967309688 112.411 > 17 21 Cd 0.08333333333333 0.41666666666667 0.41671967309688 112.411 > 17 22 Cd 0.33333333333333 0.41666666666667 0.41671967309688 112.411 > 17 23 Cd 0.58333333333333 0.41666666666667 0.41671967309688 112.411 > 17 24 Cd 0.83333333333333 0.41666666666667 0.41671967309688 112.411 > 17 25 Cd 0.08333333333333 0.66666666666667 0.41671967309688 112.411 > 17 26 Cd 0.33333333333333 0.66666666666667 0.41671967309688 112.411 > 17 27 Cd 0.58333333333333 0.66666666666667 0.41671967309688 112.411 > 17 28 Cd 0.83333333333333 0.66666666666667 0.41671967309688 112.411 > 17 29 Cd 0.08333333333333 0.91666666666667 0.41671967309688 112.411 > 17 30 Cd 0.33333333333333 0.91666666666667 0.41671967309688 112.411 > 17 31 Cd 0.58333333333333 0.91666666666667 0.41671967309688 112.411 > 17 32 Cd 0.83333333333333 0.91666666666667 0.41671967309688 112.411 > 17 33 Cd 0.08333333333333 0.16666666666667 0.74995362118498 112.411 > 33 34 Cd 0.33333333333333 0.16666666666667 0.74995362118498 112.411 > 33 35 Cd 0.58333333333333 0.16666666666667 0.74995362118498 112.411 > 33 36 Cd 0.83333333333333 0.16666666666667 0.74995362118498 112.411 > 33 37 Cd 0.08333333333333 0.41666666666667 0.74995362118498 112.411 > 33 38 Cd 0.33333333333333 0.41666666666667 0.74995362118498 112.411 > 33 39 Cd 0.58333333333333 0.41666666666667 0.74995362118498 112.411 > 33 40 Cd 0.83333333333333 0.41666666666667 0.74995362118498 112.411 > 33 41 Cd 0.08333333333333 0.66666666666667 0.74995362118498 112.411 > 33 42 Cd 0.33333333333333 0.66666666666667 0.74995362118498 112.411 > 33 43 Cd 0.58333333333333 0.66666666666667 0.74995362118498 112.411 > 33 44 Cd 0.83333333333333 0.66666666666667 0.74995362118498 112.411 > 33 45 Cd 0.08333333333333 0.91666666666667 0.74995362118498 112.411 > 33 46 Cd 0.33333333333333 0.91666666666667 0.74995362118498 112.411 > 33 47 Cd 0.58333333333333 0.91666666666667 0.74995362118498 112.411 > 33 48 Cd 0.83333333333333 0.91666666666667 0.74995362118498 112.411 > 33 49 I 0.00000000000000 0.00000000000000 0.66702061783345 126.904 > 49 50 I 0.25000000000000 0.00000000000000 0.66702061783345 126.904 > 49 51 I 0.50000000000000 0.00000000000000 0.66702061783345 126.904 > 49 52 I 0.75000000000000 0.00000000000000 0.66702061783345 126.904 > 49 53 I 0.00000000000000 0.25000000000000 0.66702061783345 126.904 > 49 54 I 0.25000000000000 0.25000000000000 0.66702061783345 126.904 > 49 55 I 0.50000000000000 0.25000000000000 0.66702061783345 126.904 > 49 56 I 0.75000000000000 0.25000000000000 0.66702061783345 126.904 > 49 57 I 0.00000000000000 0.50000000000000 0.66702061783345 126.904 > 49 58 I 0.25000000000000 0.50000000000000 0.66702061783345 126.904 > 49 59 I 0.50000000000000 0.50000000000000 0.66702061783345 126.904 > 49 60 I 0.75000000000000 0.50000000000000 0.66702061783345 126.904 > 49 61 I 0.00000000000000 0.75000000000000 0.66702061783345 126.904 > 49 62 I 0.25000000000000 0.75000000000000 0.66702061783345 126.904 > 49 63 I 0.50000000000000 0.75000000000000 0.66702061783345 126.904 > 49 64 I 0.75000000000000 0.75000000000000 0.66702061783345 126.904 > 49 65 I 0.16666666666667 0.08333333333333 0.83292910146147 126.904 > 65 66 I 0.41666666666667 0.08333333333333 0.83292910146147 126.904 > 65 67 I 0.66666666666667 0.08333333333333 0.83292910146147 126.904 > 65 68 I 0.91666666666667 0.08333333333333 0.83292910146147 126.904 > 65 69 I 0.16666666666667 0.33333333333333 0.83292910146147 126.904 > 65 70 I 0.41666666666667 0.33333333333333 0.83292910146147 126.904 > 65 71 I 0.66666666666667 0.33333333333333 0.83292910146147 126.904 > 65 72 I 0.91666666666667 0.33333333333333 0.83292910146147 126.904 > 65 73 I 0.16666666666667 0.58333333333333 0.83292910146147 126.904 > 65 74 I 0.41666666666667 0.58333333333333 0.83292910146147 126.904 > 65 75 I 0.66666666666667 0.58333333333333 0.83292910146147 126.904 > 65 76 I 0.91666666666667 0.58333333333333 0.83292910146147 126.904 > 65 77 I 0.16666666666667 0.83333333333333 0.83292910146147 126.904 > 65 78 I 0.41666666666667 0.83333333333333 0.83292910146147 126.904 > 65 79 I 0.66666666666667 0.83333333333333 0.83292910146147 126.904 > 65 80 I 0.91666666666667 0.83333333333333 0.83292910146147 126.904 > 65 81 I 0.16666666666667 0.08333333333333 0.16628926553231 126.904 > 81 82 I 0.41666666666667 0.08333333333333 0.16628926553231 126.904 > 81 83 I 0.66666666666667 0.08333333333333 0.16628926553231 126.904 > 81 84 I 0.91666666666667 0.08333333333333 0.16628926553231 126.904 > 81 85 I 0.16666666666667 0.33333333333333 0.16628926553231 126.904 > 81 86 I 0.41666666666667 0.33333333333333 0.16628926553231 126.904 > 81 87 I 0.66666666666667 0.33333333333333 0.16628926553231 126.904 > 81 88 I 0.91666666666667 0.33333333333333 0.16628926553231 126.904 > 81 89 I 0.16666666666667 0.58333333333333 0.16628926553231 126.904 > 81 90 I 0.41666666666667 0.58333333333333 0.16628926553231 126.904 > 81 91 I 0.66666666666667 0.58333333333333 0.16628926553231 126.904 > 81 92 I 0.91666666666667 0.58333333333333 0.16628926553231 126.904 > 81 93 I 0.16666666666667 0.83333333333333 0.16628926553231 126.904 > 81 94 I 0.41666666666667 0.83333333333333 0.16628926553231 126.904 > 81 95 I 0.66666666666667 0.83333333333333 0.16628926553231 126.904 > 81 96 I 0.91666666666667 0.83333333333333 0.16628926553231 126.904 > 81 97 I 0.00000000000000 0.00000000000000 0.33375366569986 126.904 > 97 98 I 0.25000000000000 0.00000000000000 0.33375366569986 126.904 > 97 99 I 0.50000000000000 0.00000000000000 0.33375366569986 126.904 > 97 100 I 0.75000000000000 0.00000000000000 0.33375366569986 126.904 > 97 101 I 0.00000000000000 0.25000000000000 0.33375366569986 126.904 > 97 102 I 0.25000000000000 0.25000000000000 0.33375366569986 126.904 > 97 103 I 0.50000000000000 0.25000000000000 0.33375366569986 126.904 > 97 104 I 0.75000000000000 0.25000000000000 0.33375366569986 126.904 > 97 105 I 0.00000000000000 0.50000000000000 0.33375366569986 126.904 > 97 106 I 0.25000000000000 0.50000000000000 0.33375366569986 126.904 > 97 107 I 0.50000000000000 0.50000000000000 0.33375366569986 126.904 > 97 108 I 0.75000000000000 0.50000000000000 0.33375366569986 126.904 > 97 109 I 0.00000000000000 0.75000000000000 0.33375366569986 126.904 > 97 110 I 0.25000000000000 0.75000000000000 0.33375366569986 126.904 > 97 111 I 0.50000000000000 0.75000000000000 0.33375366569986 126.904 > 97 112 I 0.75000000000000 0.75000000000000 0.33375366569986 126.904 > 97 113 I 0.08333333333333 0.16666666666667 0.00035804939328 126.904 > 113 114 I 0.33333333333333 0.16666666666667 0.00035804939328 126.904 > 113 115 I 0.58333333333333 0.16666666666667 0.00035804939328 126.904 > 113 116 I 0.83333333333333 0.16666666666667 0.00035804939328 126.904 > 113 117 I 0.08333333333333 0.41666666666667 0.00035804939328 126.904 > 113 118 I 0.33333333333333 0.41666666666667 0.00035804939328 126.904 > 113 119 I 0.58333333333333 0.41666666666667 0.00035804939328 126.904 > 113 120 I 0.83333333333333 0.41666666666667 0.00035804939328 126.904 > 113 121 I 0.08333333333333 0.66666666666667 0.00035804939328 126.904 > 113 122 I 0.33333333333333 0.66666666666667 0.00035804939328 126.904 > 113 123 I 0.58333333333333 0.66666666666667 0.00035804939328 126.904 > 113 124 I 0.83333333333333 0.66666666666667 0.00035804939328 126.904 > 113 125 I 0.08333333333333 0.91666666666667 0.00035804939328 126.904 > 113 126 I 0.33333333333333 0.91666666666667 0.00035804939328 126.904 > 113 127 I 0.58333333333333 0.91666666666667 0.00035804939328 126.904 > 113 128 I 0.83333333333333 0.91666666666667 0.00035804939328 126.904 > 113 129 I 0.16666666666667 0.08333333333333 0.49968299085190 126.904 > 129 130 I 0.41666666666667 0.08333333333333 0.49968299085190 126.904 > 129 131 I 0.66666666666667 0.08333333333333 0.49968299085190 126.904 > 129 132 I 0.91666666666667 0.08333333333333 0.49968299085190 126.904 > 129 133 I 0.16666666666667 0.33333333333333 0.49968299085190 126.904 > 129 134 I 0.41666666666667 0.33333333333333 0.49968299085190 126.904 > 129 135 I 0.66666666666667 0.33333333333333 0.49968299085190 126.904 > 129 136 I 0.91666666666667 0.33333333333333 0.49968299085190 126.904 > 129 137 I 0.16666666666667 0.58333333333333 0.49968299085190 126.904 > 129 138 I 0.41666666666667 0.58333333333333 0.49968299085190 126.904 > 129 139 I 0.66666666666667 0.58333333333333 0.49968299085190 126.904 > 129 140 I 0.91666666666667 0.58333333333333 0.49968299085190 126.904 > 129 141 I 0.16666666666667 0.83333333333333 0.49968299085190 126.904 > 129 142 I 0.41666666666667 0.83333333333333 0.49968299085190 126.904 > 129 143 I 0.66666666666667 0.83333333333333 0.49968299085190 126.904 > 129 144 I 0.91666666666667 0.83333333333333 0.49968299085190 126.904 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 6.0550211 -0.0000000 0.0000000 0.0000000 6.0550211 0.0000000 0.0000000 0.0000000 4.8840961 -------------------------- Born effective charges -------------------------- 1 Cd 3.0193209 -0.0000000 0.0000000 0.0000000 3.0193209 0.0000000 0.0000000 0.0000000 1.6871870 2 Cd 3.0367361 -0.0000000 0.0000000 0.0000000 3.0367361 0.0000000 0.0000000 0.0000000 1.6919068 3 Cd 3.0236693 -0.0000000 0.0000000 0.0000000 3.0236693 0.0000000 0.0000000 0.0000000 1.6833891 4 I -1.5328396 0.0000000 0.0000000 0.0000000 -1.5328396 0.0000000 0.0000000 0.0000000 -0.8455897 5 I -1.5012836 0.0000000 0.0000000 0.0000000 -1.5012836 0.0000000 0.0000000 0.0000000 -0.8456941 6 I -1.5027788 0.0000000 0.0000000 0.0000000 -1.5027788 0.0000000 0.0000000 0.0000000 -0.8497911 7 I -1.5031141 0.0000000 0.0000000 0.0000000 -1.5031141 0.0000000 0.0000000 0.0000000 -0.8375722 8 I -1.5161627 0.0000000 0.0000000 0.0000000 -1.5161627 0.0000000 0.0000000 0.0000000 -0.8360942 9 I -1.5235474 0.0000000 0.0000000 0.0000000 -1.5235474 0.0000000 0.0000000 0.0000000 -0.8477416 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000009 (zzz) -0.00000009 (zzz) -0.00000009 (zzz) Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) --------------------------- 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 ... [ 14 14 2 ] 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.61, Number of G-points: 319, Lambda: 0.15 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/48) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 0.353 ( 0.000 0.000 0.000) 0.000 0.353 ( 0.000 0.000 -0.000) 0.000 0.354 ( 0.000 0.000 -0.000) 0.000 0.354 ( 0.000 0.000 -0.000) 0.000 0.651 ( 0.000 0.000 0.000) 0.000 0.657 ( 0.000 0.000 -0.000) 0.000 1.557 ( 0.000 0.000 -0.000) 0.000 1.557 ( 0.000 0.000 -0.000) 0.000 1.573 ( 0.000 0.000 0.000) 0.000 1.573 ( 0.000 0.000 0.000) 0.000 1.662 ( 0.000 0.000 -0.000) 0.000 1.662 ( 0.000 0.000 -0.000) 0.000 2.582 ( 0.000 0.000 0.000) 0.000 2.582 ( 0.000 0.000 0.000) 0.000 2.612 ( 0.000 0.000 0.000) 0.000 2.612 ( 0.000 0.000 0.000) 0.000 2.625 ( 0.000 0.000 0.000) 0.000 2.625 ( 0.000 0.000 -0.000) 0.000 3.149 ( 0.000 0.000 -0.000) 0.000 3.151 ( 0.000 0.000 -0.000) 0.000 3.217 ( 0.000 0.000 0.000) 0.000 4.025 ( 0.000 0.000 -0.000) 0.000 4.405 ( 0.000 0.000 0.000) 0.000 4.408 ( 0.000 0.000 -0.000) 0.000 ======================= Grid point 1 (2/48) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.139 ( 7.519 4.341 0.244) 8.685 0.286 ( 12.368 7.141 0.292) 14.285 0.448 ( 7.567 4.369 -3.809) 9.531 0.459 ( 7.813 4.511 3.559) 9.699 0.526 ( 22.241 12.841 -0.205) 25.682 0.560 ( 7.292 4.210 6.329) 10.533 0.565 ( 6.400 3.695 -6.649) 9.941 0.683 ( 9.966 5.754 0.622) 11.525 0.685 ( 9.613 5.550 -0.194) 11.101 1.565 ( 0.649 0.375 -1.404) 1.591 1.587 ( 1.127 0.650 1.729) 2.164 1.671 ( 0.744 0.429 -0.336) 0.922 1.712 ( 12.223 7.057 0.513) 14.123 1.729 ( 11.982 6.918 -0.336) 13.840 1.778 ( 9.795 5.655 -0.368) 11.316 2.596 ( 1.211 0.699 0.118) 1.404 2.626 ( 1.203 0.694 -0.359) 1.434 2.638 ( 1.216 0.702 0.243) 1.425 2.823 ( 13.873 8.009 -2.034) 16.147 2.831 ( 13.655 7.884 2.117) 15.909 3.159 ( 1.922 1.110 -0.102) 2.221 3.161 ( 1.899 1.096 0.111) 2.195 3.193 ( -2.074 -1.198 -0.009) 2.395 3.856 ( 0.585 0.338 -4.343) 4.395 3.953 ( -6.042 -3.488 4.433) 8.266 4.291 ( -9.595 -5.539 -0.829) 11.110 4.295 ( -9.166 -5.292 0.670) 10.606 ======================= Grid point 2 (3/48) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.334 ( 9.362 5.405 0.354) 10.816 0.551 ( 10.993 6.347 0.259) 12.697 0.593 ( 1.229 0.709 -5.886) 6.055 0.601 ( 1.670 0.964 5.533) 5.860 0.645 ( 9.190 5.306 -2.588) 10.923 0.659 ( 9.260 5.346 2.348) 10.947 1.015 ( 20.140 11.628 -0.101) 23.256 1.051 ( 17.356 10.020 0.253) 20.043 1.055 ( 17.187 9.923 -0.085) 19.846 1.581 ( 0.630 0.364 -1.552) 1.715 1.615 ( 1.168 0.674 2.163) 2.549 1.688 ( 0.633 0.365 -0.631) 0.966 2.032 ( 14.530 8.389 -0.221) 16.779 2.038 ( 14.719 8.498 0.390) 17.001 2.067 ( 14.743 8.512 -0.200) 17.025 2.637 ( 2.459 1.420 0.297) 2.855 2.666 ( 2.269 1.310 -0.772) 2.731 2.680 ( 2.503 1.445 0.479) 2.930 2.994 ( 1.219 0.704 -3.462) 3.737 3.001 ( 1.320 0.762 3.571) 3.883 3.132 ( -3.095 -1.787 -0.070) 3.575 3.313 ( 9.534 5.504 -3.483) 11.547 3.325 ( 10.756 6.210 4.193) 13.108 3.756 ( -9.694 -5.597 -6.874) 13.136 3.847 ( -4.756 -2.746 7.086) 8.965 4.022 ( -12.639 -7.297 -2.369) 14.786 4.066 ( -7.940 -4.584 1.374) 9.271 ======================= Grid point 3 (4/48) ======================= q-point: ( 0.21 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.546 ( 9.194 5.308 0.491) 10.628 0.651 ( 4.006 2.313 -3.026) 5.528 0.671 ( 4.506 2.601 2.513) 5.778 0.778 ( 9.087 5.246 0.299) 10.497 0.844 ( 8.167 4.715 -2.004) 9.641 0.860 ( 8.373 4.834 1.728) 9.821 1.399 ( 13.548 7.822 0.819) 15.666 1.407 ( 13.790 7.962 -0.901) 15.949 1.431 ( 16.059 9.271 0.122) 18.543 1.587 ( -0.195 -0.112 -1.369) 1.388 1.634 ( 0.471 0.272 2.299) 2.363 1.695 ( -0.108 -0.062 -0.959) 0.967 2.331 ( 11.485 6.631 -0.635) 13.277 2.360 ( 12.824 7.404 1.111) 14.849 2.391 ( 13.135 7.583 -0.484) 15.174 2.703 ( 3.299 1.905 0.510) 3.843 2.726 ( 2.997 1.730 -0.940) 3.586 2.748 ( 3.334 1.925 0.435) 3.874 2.991 ( -0.485 -0.280 -1.241) 1.361 3.003 ( -0.145 -0.084 1.495) 1.505 3.063 ( -2.736 -1.580 -0.234) 3.168 3.386 ( -5.353 -3.091 -0.963) 6.256 3.423 ( -5.219 -3.013 2.818) 6.653 3.521 ( -12.855 -7.422 -4.110) 15.402 3.718 ( 1.496 0.864 6.980) 7.191 3.833 ( -1.599 -0.923 -6.449) 6.708 4.012 ( 0.553 0.319 1.687) 1.804 ======================= Grid point 4 (5/48) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.729 ( 6.987 4.034 0.736) 8.102 0.775 ( 6.642 3.835 -1.694) 7.854 0.810 ( 7.855 4.535 0.955) 9.120 0.960 ( 7.023 4.055 0.449) 8.122 1.007 ( 6.260 3.614 -1.678) 7.421 1.030 ( 6.608 3.815 1.254) 7.733 1.572 ( -1.079 -0.623 -1.083) 1.651 1.637 ( -0.186 -0.107 2.412) 2.422 1.651 ( 9.017 5.206 1.956) 10.595 1.669 ( 9.420 5.439 -2.277) 11.113 1.685 ( -0.684 -0.395 -1.361) 1.574 1.737 ( 10.925 6.307 0.335) 12.619 2.505 ( 2.752 1.589 -0.992) 3.329 2.572 ( 4.511 2.604 2.594) 5.819 2.618 ( 5.410 3.123 -1.604) 6.449 2.779 ( 3.280 1.894 0.648) 3.843 2.795 ( 3.024 1.746 -0.930) 3.613 2.824 ( 3.261 1.883 0.288) 3.776 2.987 ( 0.246 0.142 -0.183) 0.338 3.004 ( 0.229 0.132 0.497) 0.563 3.028 ( 0.364 0.210 -0.330) 0.535 3.218 ( -6.901 -3.984 0.282) 7.974 3.244 ( -7.971 -4.602 -0.917) 9.249 3.253 ( -5.963 -3.443 0.525) 6.906 3.817 ( 4.472 2.582 4.095) 6.591 3.880 ( 3.199 1.847 -5.254) 6.423 4.030 ( 0.648 0.374 1.284) 1.486 ======================= Grid point 5 (6/48) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.888 ( 7.266 4.195 1.475) 8.518 0.933 ( 7.183 4.147 -2.417) 8.639 0.995 ( 7.835 4.523 0.956) 9.097 1.095 ( 4.934 2.848 0.626) 5.731 1.125 ( 4.173 2.409 -1.419) 5.023 1.155 ( 4.390 2.534 0.821) 5.135 1.543 ( -1.372 -0.792 -0.804) 1.777 1.629 ( -0.375 -0.217 2.419) 2.457 1.668 ( -0.785 -0.453 -1.649) 1.882 1.815 ( 5.704 3.293 1.993) 6.881 1.838 ( 5.800 3.349 -2.519) 7.156 1.922 ( 5.616 3.243 0.533) 6.507 2.448 ( -6.270 -3.620 -1.055) 7.316 2.525 ( -7.210 -4.163 3.039) 8.863 2.574 ( -7.802 -4.505 -2.020) 9.233 2.846 ( 2.499 1.443 0.651) 2.958 2.857 ( 2.305 1.331 -0.819) 2.785 2.889 ( 2.425 1.400 0.172) 2.805 3.027 ( 3.195 1.845 -0.339) 3.705 3.047 ( 3.321 1.917 0.841) 3.926 3.071 ( 2.792 1.612 -0.670) 3.293 3.152 ( 1.326 0.766 0.756) 1.708 3.188 ( 1.305 0.754 -0.991) 1.804 3.217 ( 1.289 0.744 0.415) 1.545 3.895 ( 2.548 1.471 2.531) 3.881 3.937 ( 1.780 1.028 -3.448) 4.014 4.036 ( -0.133 -0.077 0.936) 0.948 ======================= Grid point 6 (7/48) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.047 ( 6.293 3.633 1.269) 7.377 1.081 ( 5.335 3.080 -1.916) 6.451 1.137 ( 4.233 2.444 0.654) 4.932 1.182 ( 2.738 1.581 0.649) 3.227 1.196 ( 2.073 1.197 -1.020) 2.602 1.227 ( 2.033 1.174 0.407) 2.383 1.515 ( -0.971 -0.560 -0.491) 1.224 1.623 ( -0.145 -0.084 1.931) 1.938 1.652 ( -0.555 -0.320 -1.478) 1.611 1.920 ( 3.838 2.216 1.142) 4.577 1.943 ( 3.663 2.115 -1.620) 4.529 2.003 ( 1.806 1.043 0.481) 2.140 2.297 ( -5.898 -3.405 -0.864) 6.865 2.343 ( -7.432 -4.291 2.051) 8.823 2.380 ( -7.721 -4.458 -1.188) 8.994 2.889 ( 1.329 0.767 0.517) 1.619 2.896 ( 1.183 0.683 -0.600) 1.492 2.931 ( 1.244 0.718 0.085) 1.439 3.095 ( 2.521 1.456 -0.239) 2.921 3.116 ( 2.365 1.365 0.670) 2.812 3.127 ( 1.810 1.045 -0.430) 2.134 3.215 ( 2.995 1.729 0.411) 3.483 3.237 ( 2.168 1.251 -0.931) 2.670 3.258 ( 1.625 0.938 0.506) 1.943 3.938 ( 1.339 0.773 1.620) 2.240 3.962 ( 0.524 0.302 -2.134) 2.218 4.026 ( -0.606 -0.350 0.523) 0.874 ======================= Grid point 7 (8/48) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 64 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.131 ( 0.000 0.000 0.000) 0.000 1.151 ( 0.000 0.000 -0.000) 0.000 1.174 ( 0.000 0.000 -0.000) 0.000 1.215 ( 0.000 0.000 0.000) 0.000 1.218 ( -0.000 -0.000 -0.000) 0.000 1.249 ( 0.000 0.000 0.000) 0.000 1.503 ( 0.000 0.000 -0.000) 0.000 1.622 ( 0.000 0.000 0.000) 0.000 1.645 ( 0.000 0.000 -0.000) 0.000 1.977 ( 0.000 0.000 -0.000) 0.000 1.993 ( 0.000 0.000 -0.000) 0.000 2.015 ( 0.000 0.000 -0.000) 0.000 2.224 ( 0.000 0.000 -0.000) 0.000 2.244 ( 0.000 0.000 -0.000) 0.000 2.282 ( 0.000 0.000 -0.000) 0.000 2.904 ( 0.000 0.000 -0.000) 0.000 2.910 ( 0.000 0.000 0.000) 0.000 2.945 ( 0.000 0.000 -0.000) 0.000 3.129 ( 0.000 0.000 0.000) 0.000 3.141 ( 0.000 0.000 0.000) 0.000 3.148 ( 0.000 0.000 -0.000) 0.000 3.254 ( 0.000 0.000 0.000) 0.000 3.265 ( 0.000 0.000 0.000) 0.000 3.275 ( 0.000 0.000 0.000) 0.000 3.956 ( 0.000 0.000 -0.000) 0.000 3.966 ( 0.000 0.000 0.000) 0.000 4.017 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 16 (9/48) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.274 ( 5.103 8.839 0.000) 10.206 0.506 ( 7.417 12.847 0.000) 14.834 0.581 ( 1.097 1.901 0.000) 2.195 0.584 ( 0.531 0.920 0.000) 1.062 0.613 ( 5.756 9.970 0.000) 11.513 0.619 ( 5.470 9.473 0.000) 10.939 0.880 ( 11.685 20.239 0.000) 23.369 0.936 ( 9.696 16.794 0.000) 19.392 0.937 ( 9.737 16.865 0.000) 19.474 1.625 ( 2.270 3.931 0.000) 4.540 1.645 ( 2.254 3.904 0.000) 4.508 1.719 ( 1.971 3.414 0.000) 3.942 1.904 ( 6.820 11.813 0.000) 13.640 1.909 ( 6.409 11.100 0.000) 12.817 1.927 ( 5.884 10.191 0.000) 11.768 2.634 ( 1.804 3.124 0.000) 3.607 2.664 ( 1.852 3.208 0.000) 3.705 2.670 ( 1.662 2.879 0.000) 3.324 2.982 ( 2.269 3.929 0.000) 4.537 2.984 ( 2.161 3.743 0.000) 4.322 3.150 ( -1.704 -2.951 0.000) 3.407 3.258 ( 5.912 10.240 0.000) 11.824 3.260 ( 5.908 10.233 0.000) 11.816 3.821 ( -5.446 -9.433 0.000) 10.892 3.876 ( 0.590 1.022 0.000) 1.180 4.101 ( -7.546 -13.069 0.000) 15.091 4.101 ( -7.602 -13.168 0.000) 15.205 ======================= Grid point 17 (10/48) ======================= q-point: ( 0.14 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.465 ( 7.383 7.730 0.204) 10.691 0.613 ( 2.506 1.787 -2.871) 4.209 0.623 ( 2.974 1.864 2.665) 4.407 0.756 ( 6.039 13.187 0.156) 14.505 0.822 ( 5.557 11.625 -1.474) 12.968 0.831 ( 5.836 11.421 1.304) 12.892 1.262 ( 12.312 11.827 0.170) 17.073 1.267 ( 12.191 11.048 -0.228) 16.454 1.271 ( 13.802 12.568 0.092) 18.667 1.680 ( -3.194 8.531 -0.987) 9.163 1.705 ( -2.122 7.680 1.402) 8.091 1.767 ( -2.597 7.028 -0.421) 7.505 2.131 ( 12.544 1.430 -0.084) 12.626 2.154 ( 13.226 2.559 -0.005) 13.471 2.165 ( 13.303 2.068 0.079) 13.463 2.719 ( 2.583 6.718 0.262) 7.203 2.743 ( 2.044 6.394 -0.692) 6.748 2.756 ( 2.636 6.608 0.437) 7.127 2.997 ( -0.512 -0.251 -2.014) 2.093 3.002 ( -0.045 -0.430 2.105) 2.149 3.086 ( -2.302 -2.669 -0.086) 3.525 3.439 ( 0.861 6.118 -2.265) 6.581 3.462 ( 3.116 6.327 3.247) 7.764 3.616 ( -5.859 -8.528 -1.856) 10.512 3.774 ( -12.067 -5.357 2.820) 13.500 3.841 ( -5.891 -10.153 -3.206) 12.168 3.945 ( 3.899 -5.961 1.242) 7.230 ======================= Grid point 18 (11/48) ======================= q-point: ( 0.21 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.658 ( 6.815 5.804 0.276) 8.956 0.704 ( 5.476 3.418 -1.247) 6.575 0.722 ( 6.127 3.092 0.961) 6.930 0.965 ( 3.469 12.924 0.276) 13.384 1.012 ( 3.186 11.771 -1.443) 12.280 1.025 ( 3.744 11.381 1.164) 12.038 1.537 ( 10.182 6.032 1.564) 11.938 1.548 ( 10.753 6.129 -1.749) 12.500 1.599 ( 12.245 7.505 0.236) 14.364 1.687 ( -6.012 9.447 -0.957) 11.239 1.725 ( -4.537 8.326 1.763) 9.645 1.772 ( -4.832 7.426 -0.844) 8.899 2.328 ( 11.877 -5.308 -0.283) 13.012 2.366 ( 12.426 -5.243 0.574) 13.499 2.379 ( 13.637 -6.316 -0.282) 15.031 2.845 ( 2.285 9.959 0.324) 10.223 2.855 ( 1.670 9.522 -0.516) 9.681 2.874 ( 1.460 9.554 0.174) 9.666 2.989 ( -0.214 0.205 -0.392) 0.491 3.002 ( 0.295 0.125 0.698) 0.768 3.031 ( -1.281 -1.515 -0.281) 2.004 3.317 ( -9.534 -3.828 -0.219) 10.276 3.363 ( -10.684 -2.948 2.006) 11.264 3.395 ( -13.772 -5.170 -2.326) 14.893 3.735 ( 6.153 0.347 5.170) 8.044 3.802 ( 5.440 -2.162 -5.827) 8.260 3.971 ( 3.855 -4.023 1.177) 5.695 ======================= Grid point 19 (12/48) ======================= q-point: ( 0.29 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.810 ( 6.469 4.264 0.968) 7.808 0.849 ( 6.942 3.903 -1.801) 8.165 0.892 ( 8.527 4.068 0.836) 9.485 1.125 ( 1.041 12.020 0.431) 12.073 1.156 ( 0.676 11.000 -1.191) 11.085 1.175 ( 1.347 10.434 0.763) 10.548 1.655 ( -6.530 7.993 -0.611) 10.340 1.698 ( -0.864 4.030 1.238) 4.304 1.709 ( 0.766 3.169 0.337) 3.278 1.765 ( 4.248 5.816 -1.185) 7.299 1.780 ( 0.801 7.213 -0.103) 7.258 1.835 ( 8.267 3.762 0.321) 9.088 2.397 ( 4.058 -8.873 -0.773) 9.787 2.445 ( 5.311 -10.269 1.687) 11.684 2.475 ( 6.652 -10.544 -0.876) 12.497 2.940 ( -2.593 7.097 0.141) 7.557 2.945 ( -2.978 6.840 -0.056) 7.461 2.957 ( -3.134 6.090 -0.154) 6.851 3.010 ( 2.302 4.786 -0.383) 5.325 3.025 ( 2.330 3.722 0.766) 4.458 3.051 ( 2.378 3.201 -0.400) 4.008 3.163 ( -4.594 -1.583 0.350) 4.871 3.192 ( -3.932 -1.220 -0.695) 4.175 3.213 ( -2.654 -1.268 0.355) 2.963 3.830 ( 5.089 -1.139 2.996) 6.015 3.875 ( 4.784 -2.118 -3.926) 6.541 3.989 ( 2.995 -4.302 0.973) 5.331 ======================= Grid point 20 (13/48) ======================= q-point: ( 0.36 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.966 ( 6.986 3.858 1.331) 8.091 1.005 ( 6.764 3.287 -2.148) 7.821 1.060 ( 6.911 2.475 0.826) 7.387 1.234 ( -1.164 10.564 0.500) 10.639 1.251 ( -1.512 9.769 -0.888) 9.925 1.272 ( -1.037 8.828 0.389) 8.897 1.622 ( -6.191 7.811 -0.355) 9.973 1.686 ( -3.722 4.795 1.299) 6.208 1.705 ( -3.038 3.234 -0.841) 4.516 1.881 ( 3.914 3.968 0.721) 5.620 1.903 ( 4.251 4.025 -1.336) 6.004 1.964 ( 4.221 1.315 0.515) 4.451 2.316 ( -2.499 -7.904 -0.930) 8.342 2.358 ( -2.575 -9.758 1.821) 10.255 2.398 ( -2.649 -9.765 -0.881) 10.156 2.938 ( -2.637 4.332 0.034) 5.072 2.941 ( -2.199 4.116 -0.088) 4.667 2.957 ( -1.262 3.491 0.031) 3.713 3.095 ( 0.596 5.360 -0.094) 5.393 3.109 ( 1.162 4.722 0.386) 4.878 3.122 ( 1.093 4.044 -0.433) 4.211 3.164 ( 3.855 0.394 0.648) 3.929 3.195 ( 2.849 -0.226 -0.979) 3.021 3.220 ( 2.688 -0.645 0.466) 2.803 3.887 ( 3.727 -2.227 1.822) 4.709 3.918 ( 3.198 -2.833 -2.493) 4.947 3.988 ( 2.148 -4.571 0.686) 5.097 ======================= Grid point 21 (14/48) ======================= q-point: ( 0.43 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.099 ( 4.401 1.840 0.581) 4.806 1.119 ( 3.601 1.101 -0.863) 3.863 1.150 ( 2.098 -0.542 0.283) 2.185 1.291 ( -3.280 8.707 0.410) 9.314 1.296 ( -3.614 8.318 -0.531) 9.085 1.315 ( -3.027 7.179 0.125) 7.792 1.604 ( -5.313 8.333 -0.115) 9.883 1.675 ( -2.829 4.373 0.572) 5.239 1.690 ( -2.452 3.383 -0.452) 4.203 1.972 ( 2.215 3.448 0.291) 4.108 1.994 ( 1.794 3.526 -0.720) 4.021 2.014 ( 0.619 0.699 0.398) 1.015 2.209 ( -1.245 -4.970 -0.741) 5.177 2.225 ( -2.202 -6.314 1.155) 6.786 2.262 ( -2.082 -6.335 -0.387) 6.680 2.944 ( -1.587 3.387 0.051) 3.741 2.949 ( -1.537 3.418 -0.068) 3.748 2.969 ( -0.869 2.383 0.011) 2.537 3.140 ( -0.116 3.134 -0.119) 3.138 3.153 ( -0.247 2.501 0.123) 2.516 3.164 ( -0.464 3.094 0.085) 3.130 3.225 ( 2.599 -0.671 0.142) 2.688 3.237 ( 2.316 -1.139 -0.704) 2.675 3.251 ( 1.832 -1.658 0.468) 2.515 3.917 ( 2.745 -2.853 0.964) 4.074 3.932 ( 2.112 -3.193 -1.226) 4.020 3.977 ( 2.007 -4.499 0.268) 4.934 ======================= Grid point 31 (15/48) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.626 ( 4.869 8.434 0.000) 9.739 0.671 ( 2.512 4.350 0.000) 5.023 0.684 ( 2.538 4.395 0.000) 5.075 1.020 ( 7.482 12.959 0.000) 14.963 1.059 ( 6.863 11.888 0.000) 13.727 1.068 ( 6.898 11.947 0.000) 13.796 1.468 ( 5.527 9.573 0.000) 11.054 1.470 ( 5.566 9.640 0.000) 11.132 1.516 ( 6.763 11.715 0.000) 13.527 1.848 ( 3.982 6.896 0.000) 7.963 1.864 ( 3.945 6.833 0.000) 7.890 1.915 ( 3.616 6.264 0.000) 7.232 2.128 ( 0.155 0.268 0.000) 0.309 2.164 ( 0.347 0.602 0.000) 0.695 2.189 ( 0.592 1.026 0.000) 1.185 2.871 ( 4.615 7.993 0.000) 9.230 2.886 ( 4.293 7.435 0.000) 8.586 2.910 ( 4.693 8.129 0.000) 9.386 2.995 ( 0.099 0.172 0.000) 0.199 2.997 ( 0.198 0.343 0.000) 0.396 3.035 ( -1.346 -2.331 0.000) 2.692 3.417 ( -5.329 -9.230 0.000) 10.658 3.486 ( -3.564 -6.174 0.000) 7.129 3.489 ( -3.486 -6.038 0.000) 6.973 3.694 ( -0.205 -0.355 0.000) 0.409 3.696 ( -0.187 -0.324 0.000) 0.375 3.884 ( -0.572 -0.991 0.000) 1.145 ======================= Grid point 32 (16/48) ======================= q-point: ( 0.21 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.764 ( 4.318 5.336 0.236) 6.868 0.790 ( 4.788 5.032 -0.639) 6.975 0.809 ( 6.475 5.965 0.402) 8.813 1.242 ( 2.660 14.464 0.098) 14.707 1.266 ( 2.706 13.240 -0.311) 13.518 1.273 ( 2.417 13.358 0.201) 13.576 1.631 ( 6.009 3.748 1.855) 7.321 1.642 ( 6.888 3.567 -2.060) 8.026 1.717 ( 7.292 4.568 0.244) 8.608 1.908 ( -7.542 12.004 -0.188) 14.178 1.921 ( -6.318 10.705 0.611) 12.446 1.949 ( -7.064 9.577 -0.467) 11.909 2.196 ( 9.398 -6.956 -0.218) 11.694 2.226 ( 10.188 -7.756 0.023) 12.804 2.232 ( 8.229 -7.662 0.220) 11.246 2.981 ( 0.791 1.710 0.028) 1.884 2.987 ( 0.867 1.692 0.010) 1.901 3.001 ( -0.492 -0.475 -0.096) 0.691 3.032 ( 2.558 4.954 0.059) 5.575 3.038 ( 2.508 5.029 0.033) 5.620 3.053 ( 3.003 5.977 -0.031) 6.689 3.253 ( -6.379 -1.979 0.016) 6.679 3.317 ( -8.094 -1.275 1.939) 8.420 3.323 ( -8.530 -1.702 -2.057) 8.938 3.730 ( 3.814 -0.622 3.395) 5.144 3.756 ( 5.195 -2.330 -3.795) 6.842 3.878 ( 2.887 -4.990 0.494) 5.786 ======================= Grid point 33 (17/48) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 200 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.901 ( 5.645 5.090 0.903) 7.654 0.932 ( 6.033 4.631 -1.529) 7.757 0.974 ( 6.651 4.343 0.631) 7.969 1.390 ( -1.586 13.906 0.094) 13.997 1.402 ( -1.758 12.999 -0.241) 13.119 1.406 ( -1.663 12.483 0.132) 12.594 1.731 ( 2.621 0.540 1.438) 3.037 1.749 ( 2.222 0.857 -1.823) 2.999 1.811 ( -2.291 3.400 0.518) 4.132 1.912 ( -1.273 7.636 -0.226) 7.745 1.924 ( -1.167 7.320 0.149) 7.414 1.942 ( -2.769 7.991 -0.045) 8.457 2.239 ( 2.867 -6.021 -0.889) 6.728 2.265 ( 3.448 -6.537 1.223) 7.491 2.302 ( 4.510 -5.716 -0.321) 7.288 2.979 ( -0.647 0.335 -0.222) 0.762 2.991 ( -0.713 0.432 0.521) 0.983 3.002 ( -0.790 0.495 -0.342) 0.993 3.097 ( 2.567 -0.833 0.413) 2.731 3.113 ( 2.983 -0.054 -0.577) 3.039 3.134 ( 3.024 0.357 0.223) 3.053 3.196 ( -3.057 6.393 -0.058) 7.087 3.231 ( -4.550 6.560 0.333) 7.991 3.243 ( -4.454 6.908 -0.326) 8.226 3.782 ( 4.533 -3.404 1.818) 5.954 3.809 ( 4.843 -4.241 -2.373) 6.861 3.877 ( 3.685 -6.683 0.580) 7.654 ======================= Grid point 34 (18/48) ======================= q-point: ( 0.36 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.041 ( 5.226 3.736 0.764) 6.470 1.067 ( 4.752 3.094 -1.316) 5.821 1.102 ( 4.181 1.896 0.561) 4.625 1.459 ( -3.925 10.168 -0.004) 10.899 1.466 ( -4.238 10.770 0.152) 11.574 1.468 ( -4.368 11.288 -0.158) 12.105 1.743 ( -0.929 1.132 0.524) 1.556 1.752 ( -1.887 1.717 -0.597) 2.620 1.780 ( -4.308 4.280 0.079) 6.073 1.986 ( 1.567 6.460 -0.059) 6.648 2.004 ( 2.051 5.035 -0.417) 5.453 2.017 ( 2.283 5.830 0.529) 6.284 2.190 ( -2.381 -3.557 -0.804) 4.355 2.216 ( -2.841 -3.235 1.010) 4.422 2.261 ( -3.497 -3.051 -0.258) 4.648 2.980 ( 0.163 0.687 -0.055) 0.708 2.983 ( -0.411 0.572 0.092) 0.711 2.997 ( -0.298 0.587 -0.059) 0.661 3.133 ( 3.586 -2.181 0.387) 4.215 3.157 ( 3.603 -2.439 -0.773) 4.419 3.174 ( 3.139 -2.754 0.388) 4.194 3.220 ( -2.247 6.283 -0.046) 6.672 3.246 ( -2.796 7.068 0.411) 7.612 3.259 ( -3.007 7.131 -0.360) 7.747 3.818 ( 4.248 -4.577 0.869) 6.305 3.838 ( 3.945 -4.997 -1.290) 6.495 3.870 ( 3.490 -7.077 0.432) 7.903 ======================= Grid point 35 (19/48) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.119 ( -0.013 0.023 -0.000) 0.026 1.132 ( 0.069 -0.119 0.000) 0.138 1.135 ( 0.183 -0.317 0.000) 0.366 1.473 ( -4.754 8.235 0.000) 9.508 1.480 ( -5.489 9.507 0.000) 10.978 1.484 ( -5.555 9.622 -0.000) 11.110 1.741 ( -0.897 1.554 0.000) 1.794 1.743 ( -1.256 2.176 -0.000) 2.512 1.763 ( -2.823 4.890 0.000) 5.647 2.055 ( -3.176 5.501 0.000) 6.352 2.067 ( -2.039 3.531 -0.000) 4.078 2.088 ( -2.182 3.779 0.000) 4.364 2.132 ( -0.325 0.562 0.000) 0.650 2.152 ( -0.068 0.117 -0.000) 0.136 2.181 ( 0.282 -0.489 -0.000) 0.564 2.983 ( -0.276 0.478 0.000) 0.552 2.989 ( -0.226 0.392 -0.000) 0.452 2.997 ( -0.194 0.337 0.000) 0.389 3.159 ( 1.999 -3.463 -0.000) 3.998 3.178 ( 2.163 -3.746 -0.000) 4.326 3.185 ( 2.018 -3.496 -0.000) 4.037 3.233 ( -2.707 4.688 0.000) 5.413 3.258 ( -3.317 5.745 -0.000) 6.634 3.268 ( -3.424 5.930 0.000) 6.848 3.833 ( 3.125 -5.413 -0.000) 6.251 3.847 ( 3.040 -5.265 -0.000) 6.080 3.862 ( 3.963 -6.864 -0.000) 7.926 ======================= Grid point 47 (20/48) ======================= q-point: ( 0.21 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.878 ( 3.671 6.359 0.000) 7.342 0.902 ( 3.655 6.331 0.000) 7.310 0.936 ( 3.849 6.667 0.000) 7.698 1.491 ( 5.709 9.888 0.000) 11.417 1.502 ( 5.910 10.236 0.000) 11.819 1.513 ( 5.776 10.005 0.000) 11.553 1.692 ( 1.491 2.582 0.000) 2.981 1.695 ( 1.411 2.444 0.000) 2.822 1.782 ( 1.328 2.301 0.000) 2.657 2.036 ( -1.230 -2.130 0.000) 2.459 2.046 ( -1.430 -2.477 0.000) 2.860 2.068 ( -1.549 -2.683 0.000) 3.098 2.142 ( 3.258 5.643 0.000) 6.516 2.161 ( 3.042 5.269 0.000) 6.085 2.177 ( 3.765 6.522 0.000) 7.531 2.989 ( -0.103 -0.179 0.000) 0.207 2.998 ( 0.182 0.315 0.000) 0.364 3.005 ( 0.263 0.456 0.000) 0.526 3.075 ( -0.209 -0.363 0.000) 0.419 3.078 ( -0.378 -0.655 0.000) 0.756 3.094 ( -0.683 -1.184 0.000) 1.367 3.275 ( 1.914 3.315 0.000) 3.828 3.330 ( 1.363 2.361 0.000) 2.726 3.341 ( 1.460 2.529 0.000) 2.920 3.719 ( -0.471 -0.816 0.000) 0.943 3.721 ( -0.327 -0.566 0.000) 0.654 3.792 ( -2.034 -3.523 0.000) 4.068 ======================= Grid point 48 (21/48) ======================= q-point: ( 0.29 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.005 ( 4.159 5.369 0.399) 6.803 1.028 ( 4.243 5.074 -0.720) 6.654 1.058 ( 3.993 4.176 0.326) 5.787 1.600 ( -2.177 7.148 -0.135) 7.474 1.621 ( -2.231 7.650 0.834) 8.013 1.632 ( -1.966 8.088 -0.860) 8.368 1.731 ( 1.558 0.079 1.652) 2.272 1.751 ( 2.142 -0.196 -1.855) 2.840 1.809 ( 1.121 -1.946 0.370) 2.276 2.016 ( 0.738 2.787 0.005) 2.883 2.024 ( 0.621 2.240 -0.154) 2.329 2.042 ( 1.377 4.420 0.167) 4.633 2.237 ( -0.597 5.752 -0.662) 5.820 2.255 ( 0.202 4.826 0.856) 4.906 2.294 ( 1.460 3.925 -0.223) 4.193 2.982 ( -0.535 -0.019 -0.118) 0.548 2.991 ( -0.877 -0.389 0.205) 0.981 3.002 ( -0.702 -0.459 -0.117) 0.847 3.054 ( 1.514 -2.896 0.296) 3.281 3.073 ( 2.210 -2.844 -0.594) 3.650 3.084 ( 2.478 -3.756 0.334) 4.512 3.327 ( -1.448 5.910 -0.103) 6.085 3.371 ( -2.130 6.460 0.563) 6.826 3.385 ( -2.009 6.405 -0.477) 6.729 3.711 ( 2.228 -3.626 0.601) 4.298 3.724 ( 2.869 -4.127 -0.895) 5.105 3.745 ( 2.206 -6.409 0.305) 6.785 ======================= Grid point 49 (22/48) ======================= q-point: ( 0.36 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.107 ( 2.460 2.784 0.109) 3.717 1.123 ( 1.871 2.478 -0.355) 3.125 1.134 ( 1.164 1.421 0.253) 1.854 1.607 ( -3.018 5.049 -0.080) 5.882 1.627 ( -3.075 4.726 0.749) 5.688 1.635 ( -3.663 4.579 -0.775) 5.915 1.743 ( 0.697 -0.707 1.117) 1.494 1.754 ( 0.020 -1.127 -1.085) 1.564 1.791 ( 0.032 -2.014 0.053) 2.015 2.088 ( 0.978 5.289 0.003) 5.378 2.104 ( 1.919 5.914 -0.198) 6.221 2.155 ( 2.547 7.780 0.441) 8.198 2.236 ( -5.871 7.026 -0.749) 9.186 2.259 ( -4.744 5.416 0.789) 7.243 2.287 ( -6.130 4.219 -0.276) 7.447 2.976 ( 0.201 -1.075 0.012) 1.094 2.979 ( 0.550 -0.708 -0.008) 0.897 2.990 ( 0.280 -1.059 -0.010) 1.095 3.066 ( 3.377 -4.423 0.109) 5.566 3.086 ( 3.193 -4.715 -0.316) 5.704 3.091 ( 2.996 -4.678 0.220) 5.559 3.345 ( -3.116 6.167 -0.057) 6.910 3.384 ( -3.645 6.634 0.239) 7.573 3.398 ( -3.702 6.797 -0.187) 7.742 3.717 ( 3.757 -5.893 0.053) 6.989 3.727 ( 3.528 -6.351 -0.096) 7.266 3.738 ( 3.382 -5.539 0.045) 6.490 ======================= Grid point 62 (23/48) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.099 ( 2.255 3.905 0.000) 4.509 1.119 ( 2.150 3.724 0.000) 4.300 1.129 ( 1.602 2.775 0.000) 3.205 1.671 ( 0.414 0.717 0.000) 0.828 1.678 ( 0.022 0.038 0.000) 0.044 1.680 ( 0.104 0.181 0.000) 0.209 1.736 ( -0.758 -1.313 0.000) 1.516 1.738 ( -0.721 -1.249 0.000) 1.442 1.765 ( -1.183 -2.049 0.000) 2.366 2.117 ( 4.309 7.463 0.000) 8.617 2.125 ( 4.508 7.807 0.000) 9.015 2.174 ( 4.779 8.278 0.000) 9.559 2.329 ( 1.156 2.003 0.000) 2.313 2.342 ( 1.138 1.971 0.000) 2.275 2.373 ( 1.462 2.531 0.000) 2.923 2.971 ( -0.918 -1.590 0.000) 1.836 2.974 ( -0.574 -0.995 0.000) 1.149 2.982 ( -0.842 -1.459 0.000) 1.685 3.000 ( -1.076 -1.863 0.000) 2.152 3.020 ( -1.268 -2.196 0.000) 2.536 3.032 ( -0.656 -1.135 0.000) 1.311 3.419 ( 2.007 3.476 0.000) 4.013 3.468 ( 1.941 3.362 0.000) 3.883 3.483 ( 2.045 3.543 0.000) 4.091 3.636 ( -2.363 -4.093 0.000) 4.726 3.650 ( -1.452 -2.515 0.000) 2.904 3.661 ( -1.367 -2.368 0.000) 2.734 ======================= Grid point 63 (24/48) ======================= q-point: ( 0.36 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.142 ( -0.443 0.768 0.000) 0.887 1.154 ( -0.454 0.786 0.000) 0.907 1.160 ( -0.538 0.931 -0.000) 1.075 1.673 ( -1.062 1.840 -0.000) 2.124 1.676 ( -1.053 1.824 0.000) 2.106 1.681 ( -0.953 1.650 0.000) 1.906 1.714 ( 1.085 -1.879 0.000) 2.170 1.720 ( 1.051 -1.821 -0.000) 2.102 1.740 ( 1.525 -2.641 0.000) 3.049 2.223 ( -3.771 6.532 -0.000) 7.542 2.237 ( -3.498 6.059 0.000) 6.997 2.290 ( -2.663 4.612 -0.000) 5.326 2.326 ( -0.291 0.505 0.000) 0.583 2.339 ( -0.632 1.095 -0.000) 1.264 2.369 ( -2.075 3.594 0.000) 4.150 2.951 ( 0.667 -1.155 -0.000) 1.334 2.960 ( 0.497 -0.861 -0.000) 0.995 2.967 ( 0.633 -1.096 0.000) 1.266 2.983 ( 1.642 -2.844 -0.000) 3.284 2.998 ( 1.618 -2.803 0.000) 3.236 3.023 ( 1.197 -2.073 -0.000) 2.393 3.458 ( -2.936 5.085 0.000) 5.871 3.499 ( -2.841 4.921 0.000) 5.682 3.518 ( -3.042 5.269 -0.000) 6.084 3.594 ( 3.226 -5.587 -0.000) 6.451 3.624 ( 2.598 -4.501 0.000) 5.197 3.637 ( 2.745 -4.754 0.000) 5.489 ======================= Grid point 211 (25/48) ======================= q-point: ( 0.00 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.196 ( 0.000 -0.000 -0.000) 0.000 0.196 ( 0.000 0.000 -0.000) 0.000 0.200 ( 0.000 -0.000 0.000) 0.000 0.200 ( 0.000 0.000 0.000) 0.000 0.371 ( 0.000 0.000 -0.000) 0.000 0.374 ( 0.000 0.000 0.000) 0.000 0.415 ( 0.000 -0.000 0.000) 0.000 0.415 ( 0.000 0.000 0.000) 0.000 0.761 ( -0.000 0.000 0.000) 0.000 1.532 ( 0.000 -0.000 0.000) 0.000 1.532 ( 0.000 0.000 0.000) 0.000 1.616 ( 0.000 0.000 0.000) 0.000 1.616 ( 0.000 0.000 0.000) 0.000 1.642 ( 0.000 0.000 -0.000) 0.000 1.642 ( 0.000 -0.000 -0.000) 0.000 2.591 ( 0.000 -0.000 0.000) 0.000 2.591 ( 0.000 0.000 0.000) 0.000 2.598 ( 0.000 -0.000 0.000) 0.000 2.598 ( 0.000 -0.000 0.000) 0.000 2.630 ( -0.000 -0.000 0.000) 0.000 2.630 ( -0.000 0.000 0.000) 0.000 3.127 ( 0.000 0.000 0.000) 0.000 3.193 ( -0.000 0.000 0.000) 0.000 3.198 ( 0.000 0.000 -0.000) 0.000 4.399 ( 0.000 0.000 0.000) 0.000 4.400 ( 0.000 0.000 -0.000) 0.000 4.410 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 213 (26/48) ======================= q-point: ( 0.07 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.342 ( 10.094 5.828 -4.076) 12.348 0.343 ( 2.329 1.344 -10.529) 10.867 0.346 ( 2.780 1.605 10.374) 10.860 0.353 ( 10.191 5.884 4.303) 12.529 0.503 ( 6.897 3.982 -0.267) 7.969 0.577 ( 19.535 11.278 3.739) 22.865 0.585 ( 19.620 11.328 -3.790) 22.971 0.659 ( 16.828 9.716 0.426) 19.436 0.718 ( -2.584 -1.492 -0.298) 2.998 1.544 ( 0.946 0.546 0.392) 1.160 1.624 ( 0.597 0.345 -1.263) 1.439 1.655 ( 0.979 0.565 0.883) 1.435 1.695 ( 12.923 7.461 -0.190) 14.923 1.752 ( 10.397 6.003 -0.069) 12.006 1.779 ( 10.512 6.069 0.442) 12.146 2.604 ( 1.158 0.668 -0.466) 1.416 2.612 ( 1.263 0.729 0.559) 1.562 2.644 ( 1.208 0.698 -0.095) 1.398 2.780 ( 10.748 6.205 0.103) 12.411 3.066 ( 14.560 8.406 -13.442) 21.526 3.071 ( 14.626 8.445 13.271) 21.479 3.140 ( 1.940 1.120 0.001) 2.241 3.223 ( 8.959 5.172 -2.321) 10.602 3.227 ( 9.153 5.284 2.695) 10.907 4.241 ( -11.699 -6.754 -4.137) 14.128 4.250 ( -10.693 -6.174 3.875) 12.941 4.302 ( -8.993 -5.192 -0.071) 10.384 ======================= Grid point 215 (27/48) ======================= q-point: ( 0.14 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.426 ( 5.400 3.117 -5.935) 8.608 0.439 ( 5.827 3.364 6.091) 9.076 0.579 ( 10.347 5.974 -2.261) 12.160 0.591 ( 10.351 5.976 2.601) 12.232 0.678 ( -0.447 -0.258 -0.152) 0.539 0.686 ( 8.716 5.032 -0.358) 10.071 1.029 ( 19.182 11.075 0.918) 22.169 1.034 ( 18.903 10.914 -0.988) 21.850 1.051 ( 16.937 9.779 0.002) 19.557 1.568 ( 1.007 0.581 0.754) 1.386 1.638 ( 0.499 0.288 -1.856) 1.943 1.678 ( 0.932 0.538 1.123) 1.555 2.027 ( 14.703 8.489 0.022) 16.978 2.044 ( 14.991 8.655 -0.190) 17.311 2.068 ( 14.170 8.181 0.198) 16.363 2.643 ( 2.268 1.309 -0.658) 2.700 2.656 ( 2.590 1.495 0.898) 3.123 2.685 ( 2.371 1.369 -0.244) 2.749 2.950 ( 2.902 1.676 0.189) 3.356 3.090 ( -2.383 -1.376 -2.333) 3.607 3.096 ( -2.267 -1.309 2.107) 3.360 3.277 ( 9.625 5.557 0.185) 11.115 3.484 ( 5.341 3.083 -10.665) 12.320 3.522 ( 7.720 4.457 11.990) 14.940 3.945 ( -14.040 -8.106 -4.315) 16.777 4.001 ( -10.069 -5.813 4.056) 12.314 4.066 ( -9.273 -5.354 -1.180) 10.772 ======================= Grid point 217 (28/48) ======================= q-point: ( 0.21 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.573 ( 7.229 4.174 -2.158) 8.622 0.596 ( 7.670 4.428 2.608) 9.233 0.701 ( 2.653 1.532 -0.428) 3.093 0.796 ( 8.725 5.038 -1.590) 10.200 0.809 ( 8.915 5.147 2.011) 10.489 0.877 ( 7.979 4.607 -0.444) 9.224 1.385 ( 12.946 7.474 -0.195) 14.950 1.421 ( 15.184 8.766 0.484) 17.539 1.428 ( 15.373 8.876 -0.330) 17.755 1.583 ( 0.239 0.138 1.116) 1.149 1.641 ( -0.335 -0.194 -2.197) 2.231 1.692 ( 0.273 0.157 1.110) 1.154 2.331 ( 11.776 6.799 0.603) 13.611 2.364 ( 12.186 7.035 -1.279) 14.129 2.388 ( 13.452 7.766 0.685) 15.548 2.704 ( 3.067 1.771 -0.548) 3.584 2.725 ( 3.411 1.969 0.945) 4.050 2.748 ( 3.149 1.818 -0.402) 3.658 2.977 ( 0.604 0.349 0.281) 0.752 3.034 ( -2.262 -1.306 -1.570) 3.048 3.046 ( -1.834 -1.059 1.268) 2.469 3.382 ( -4.345 -2.509 0.675) 5.062 3.432 ( -8.587 -4.958 -2.833) 10.313 3.522 ( -9.306 -5.373 4.777) 11.759 3.684 ( -1.409 -0.813 -5.100) 5.353 3.898 ( 0.932 0.538 6.923) 7.006 3.973 ( 0.077 0.044 -4.404) 4.405 ======================= Grid point 219 (29/48) ======================= q-point: ( 0.29 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.732 ( 6.541 3.776 -0.898) 7.606 0.770 ( 7.691 4.440 1.709) 9.044 0.812 ( 7.223 4.170 -0.807) 8.379 0.970 ( 6.663 3.847 -1.131) 7.776 0.989 ( 7.003 4.043 1.719) 8.268 1.039 ( 6.222 3.592 -0.613) 7.211 1.578 ( -0.589 -0.340 1.545) 1.688 1.623 ( -1.131 -0.653 -2.457) 2.782 1.629 ( 8.873 5.123 -0.457) 10.256 1.692 ( -0.224 -0.129 0.943) 0.978 1.706 ( 10.058 5.807 1.991) 11.783 1.721 ( 10.472 6.046 -1.549) 12.191 2.514 ( 3.150 1.819 1.603) 3.975 2.557 ( 3.634 2.098 -2.743) 5.013 2.624 ( 5.849 3.377 1.142) 6.850 2.775 ( 3.110 1.796 -0.381) 3.612 2.802 ( 3.307 1.909 0.811) 3.904 2.821 ( 3.145 1.816 -0.436) 3.658 2.986 ( 0.198 0.114 0.157) 0.277 3.006 ( 0.263 0.152 -0.531) 0.611 3.027 ( 0.417 0.241 0.390) 0.620 3.223 ( -7.176 -4.143 -0.627) 8.310 3.235 ( -6.298 -3.636 0.907) 7.329 3.259 ( -7.522 -4.343 -0.168) 8.687 3.784 ( 5.193 2.998 -1.704) 6.234 3.943 ( 2.022 1.167 5.237) 5.734 3.998 ( 1.264 0.730 -3.658) 3.939 ======================= Grid point 221 (30/48) ======================= q-point: ( 0.36 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.882 ( 6.979 4.030 -1.051) 8.128 0.946 ( 7.685 4.437 2.463) 9.209 0.988 ( 7.612 4.395 -1.426) 8.905 1.097 ( 4.605 2.659 -0.727) 5.367 1.123 ( 4.812 2.778 1.477) 5.750 1.155 ( 4.070 2.350 -0.778) 4.764 1.561 ( -0.844 -0.487 2.028) 2.250 1.594 ( -1.323 -0.764 -2.688) 3.092 1.685 ( -0.355 -0.205 0.694) 0.805 1.794 ( 5.961 3.442 -0.558) 6.906 1.879 ( 5.493 3.171 2.601) 6.855 1.901 ( 5.651 3.263 -2.050) 6.840 2.460 ( -6.467 -3.734 1.873) 7.699 2.504 ( -6.826 -3.941 -3.124) 8.478 2.584 ( -7.968 -4.600 1.287) 9.290 2.838 ( 2.390 1.380 -0.234) 2.769 2.869 ( 2.480 1.432 0.580) 2.922 2.884 ( 2.359 1.362 -0.350) 2.746 3.029 ( 3.281 1.894 0.422) 3.812 3.046 ( 2.956 1.707 -0.991) 3.555 3.069 ( 3.041 1.755 0.724) 3.585 3.156 ( 1.302 0.752 -0.899) 1.752 3.185 ( 1.351 0.780 1.126) 1.924 3.217 ( 1.287 0.743 -0.393) 1.538 3.876 ( 2.978 1.720 -1.073) 3.603 3.977 ( 0.956 0.552 3.479) 3.650 4.016 ( 0.270 0.156 -2.424) 2.444 ======================= Grid point 223 (31/48) ======================= q-point: ( 0.43 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.038 ( 6.313 3.645 -0.716) 7.325 1.098 ( 5.292 3.055 1.962) 6.418 1.128 ( 4.256 2.457 -1.253) 5.071 1.178 ( 2.549 1.472 -0.319) 2.961 1.205 ( 2.462 1.422 1.056) 3.033 1.221 ( 1.798 1.038 -0.772) 2.215 1.545 ( -0.408 -0.235 2.339) 2.386 1.567 ( -1.009 -0.582 -2.688) 2.929 1.678 ( -0.218 -0.126 0.387) 0.462 1.908 ( 4.306 2.486 -0.430) 4.991 1.968 ( 2.699 1.558 1.860) 3.630 1.989 ( 2.340 1.351 -1.433) 3.058 2.301 ( -6.385 -3.687 1.103) 7.455 2.337 ( -6.620 -3.822 -2.048) 7.914 2.384 ( -8.092 -4.672 0.945) 9.392 2.880 ( 1.254 0.724 -0.113) 1.453 2.912 ( 1.301 0.751 0.327) 1.538 2.924 ( 1.200 0.693 -0.216) 1.402 3.099 ( 2.618 1.511 0.454) 3.056 3.109 ( 2.236 1.291 -0.674) 2.668 3.130 ( 1.894 1.093 0.219) 2.198 3.216 ( 2.869 1.657 -0.476) 3.347 3.235 ( 2.206 1.274 0.953) 2.720 3.259 ( 1.672 0.965 -0.464) 1.986 3.925 ( 1.451 0.838 -0.617) 1.786 3.990 ( 0.317 0.183 2.194) 2.225 4.012 ( -0.513 -0.296 -1.587) 1.694 ======================= Grid point 225 (32/48) ======================= q-point: (-0.50 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 57 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.127 ( 0.000 0.000 0.000) 0.000 1.159 ( 0.000 0.000 0.000) 0.000 1.171 ( 0.000 0.000 0.000) 0.000 1.208 ( 0.000 0.000 0.000) 0.000 1.234 ( 0.000 0.000 0.000) 0.000 1.239 ( 0.000 0.000 -0.000) 0.000 1.542 ( 0.000 0.000 0.000) 0.000 1.554 ( 0.000 0.000 -0.000) 0.000 1.676 ( 0.000 0.000 0.000) 0.000 1.973 ( 0.000 0.000 0.000) 0.000 2.002 ( -0.000 -0.000 -0.000) 0.000 2.010 ( 0.000 0.000 0.000) 0.000 2.218 ( 0.000 0.000 0.000) 0.000 2.256 ( 0.000 0.000 -0.000) 0.000 2.277 ( 0.000 0.000 0.000) 0.000 2.894 ( 0.000 0.000 0.000) 0.000 2.926 ( 0.000 0.000 0.000) 0.000 2.938 ( 0.000 0.000 0.000) 0.000 3.133 ( 0.000 0.000 0.000) 0.000 3.136 ( 0.000 0.000 0.000) 0.000 3.149 ( 0.000 0.000 0.000) 0.000 3.253 ( 0.000 0.000 -0.000) 0.000 3.263 ( 0.000 0.000 0.000) 0.000 3.277 ( -0.000 -0.000 -0.000) 0.000 3.942 ( 0.000 0.000 0.000) 0.000 3.995 ( 0.000 0.000 0.000) 0.000 4.001 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 243 (33/48) ======================= q-point: ( 0.07 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.392 ( 2.426 4.202 0.000) 4.852 0.400 ( 2.661 4.609 -0.000) 5.322 0.541 ( 6.961 12.057 0.000) 13.923 0.545 ( 7.028 12.172 0.000) 14.055 0.647 ( 5.450 9.441 0.000) 10.901 0.680 ( -0.665 -1.151 0.000) 1.329 0.902 ( 10.961 18.985 -0.000) 21.922 0.905 ( 11.002 19.056 0.000) 22.004 0.936 ( 9.540 16.525 0.000) 19.081 1.609 ( 2.423 4.196 0.000) 4.845 1.674 ( 2.019 3.496 0.000) 4.037 1.706 ( 2.098 3.634 -0.000) 4.197 1.897 ( 7.235 12.531 0.000) 14.470 1.915 ( 5.970 10.341 0.000) 11.940 1.931 ( 5.759 9.974 0.000) 11.517 2.642 ( 1.733 3.001 0.000) 3.465 2.648 ( 1.869 3.237 0.000) 3.738 2.678 ( 1.721 2.981 0.000) 3.442 2.924 ( 3.049 5.281 0.000) 6.098 3.104 ( -1.080 -1.871 0.000) 2.161 3.109 ( -1.087 -1.883 -0.000) 2.174 3.222 ( 4.964 8.597 0.000) 9.927 3.454 ( 6.579 11.395 -0.000) 13.158 3.457 ( 6.538 11.324 0.000) 13.076 4.050 ( -6.838 -11.843 0.000) 13.675 4.052 ( -6.783 -11.748 -0.000) 13.566 4.111 ( -7.708 -13.351 0.000) 15.417 ======================= Grid point 245 (34/48) ======================= q-point: ( 0.14 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.508 ( 5.425 5.302 -1.810) 7.798 0.524 ( 5.787 5.512 1.983) 8.234 0.672 ( 1.365 0.230 -0.169) 1.394 0.776 ( 5.748 12.772 -1.481) 14.084 0.783 ( 6.122 12.512 1.663) 14.028 0.850 ( 5.564 10.954 -0.168) 12.288 1.253 ( 11.559 10.615 -0.049) 15.694 1.269 ( 13.399 12.537 0.053) 18.350 1.274 ( 13.471 12.403 -0.038) 18.311 1.670 ( -2.789 8.614 0.428) 9.064 1.724 ( -2.832 7.575 -1.155) 8.169 1.759 ( -2.260 7.058 0.732) 7.447 2.138 ( 12.281 2.240 0.159) 12.484 2.148 ( 13.383 1.019 -0.088) 13.422 2.166 ( 13.356 2.747 -0.062) 13.636 2.721 ( 2.318 6.435 -0.378) 6.850 2.738 ( 2.644 6.968 0.644) 7.480 2.759 ( 2.297 6.342 -0.273) 6.751 2.970 ( 0.619 0.549 0.144) 0.840 3.056 ( -1.972 -1.993 -1.605) 3.230 3.061 ( -1.637 -2.091 1.456) 3.028 3.417 ( 3.413 7.773 0.440) 8.500 3.526 ( -4.701 -0.077 -4.387) 6.430 3.569 ( -0.348 -1.418 4.681) 4.903 3.773 ( -13.767 -8.914 -2.120) 16.537 3.868 ( -1.524 -6.096 4.218) 7.568 3.919 ( 1.535 -8.157 -2.812) 8.763 ======================= Grid point 247 (35/48) ======================= q-point: ( 0.21 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.664 ( 6.411 5.120 -0.552) 8.223 0.687 ( 6.859 4.926 0.923) 8.495 0.734 ( 5.115 2.229 -0.359) 5.591 0.977 ( 3.223 12.641 -1.165) 13.098 0.989 ( 3.722 12.359 1.482) 12.992 1.036 ( 3.455 11.077 -0.315) 11.608 1.517 ( 10.040 5.789 -0.290) 11.593 1.578 ( 11.351 6.879 1.312) 13.337 1.588 ( 11.868 7.028 -1.071) 13.834 1.685 ( -5.536 9.434 0.832) 10.970 1.727 ( -5.598 8.341 -1.671) 10.184 1.772 ( -4.259 7.458 0.875) 8.633 2.336 ( 11.675 -4.962 0.599) 12.700 2.354 ( 12.910 -6.004 -0.710) 14.256 2.384 ( 13.331 -5.946 0.101) 14.598 2.844 ( 2.303 10.025 -0.201) 10.288 2.861 ( 1.545 9.799 0.555) 9.936 2.871 ( 1.583 9.304 -0.337) 9.444 2.984 ( 0.422 0.408 0.188) 0.616 3.011 ( -1.196 -0.878 -0.752) 1.664 3.026 ( -0.495 -0.819 0.539) 1.099 3.325 ( -9.191 -3.181 0.742) 9.754 3.341 ( -10.891 -4.257 -1.495) 11.788 3.415 ( -14.786 -4.727 1.350) 15.582 3.690 ( 7.626 0.466 -2.199) 7.950 3.880 ( 4.253 -2.219 5.688) 7.441 3.931 ( 4.528 -3.809 -4.067) 7.179 ======================= Grid point 249 (36/48) ======================= q-point: ( 0.29 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.809 ( 6.022 4.121 -0.887) 7.351 0.853 ( 7.736 4.211 1.889) 9.008 0.889 ( 8.167 3.898 -1.004) 9.106 1.130 ( 0.724 11.828 -0.781) 11.876 1.148 ( 1.310 11.387 1.255) 11.531 1.179 ( 1.028 10.239 -0.477) 10.302 1.663 ( -5.277 7.419 1.069) 9.167 1.686 ( -3.313 5.286 -1.289) 6.371 1.700 ( 2.091 3.581 -0.569) 4.186 1.770 ( -0.284 6.765 0.858) 6.825 1.806 ( 5.858 4.778 1.373) 7.684 1.818 ( 7.609 4.162 -1.440) 8.792 2.404 ( 4.132 -9.228 1.210) 10.183 2.437 ( 5.183 -9.498 -1.959) 10.996 2.477 ( 6.697 -10.946 0.711) 12.852 2.939 ( -2.929 7.480 -0.088) 8.034 2.949 ( -2.906 5.919 -0.138) 6.595 2.954 ( -2.867 6.503 0.295) 7.114 3.010 ( 2.816 5.047 0.312) 5.788 3.029 ( 2.045 3.831 -0.883) 4.431 3.047 ( 2.089 2.839 0.586) 3.573 3.169 ( -4.619 -1.370 -0.579) 4.853 3.187 ( -3.523 -1.431 0.860) 3.899 3.215 ( -3.015 -1.182 -0.289) 3.251 3.805 ( 5.740 -0.695 -1.180) 5.901 3.924 ( 3.660 -2.998 3.913) 6.140 3.964 ( 3.527 -3.852 -2.777) 5.915 ======================= Grid point 251 (37/48) ======================= q-point: ( 0.36 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.960 ( 6.816 3.961 -0.895) 7.935 1.018 ( 7.040 3.062 2.169) 7.978 1.053 ( 6.800 2.598 -1.283) 7.392 1.233 ( -1.424 10.561 -0.404) 10.664 1.254 ( -1.023 9.749 0.924) 9.846 1.270 ( -1.272 8.860 -0.520) 8.966 1.639 ( -5.440 7.285 1.438) 9.205 1.655 ( -5.207 5.944 -1.744) 8.092 1.718 ( -2.075 2.492 0.201) 3.249 1.873 ( 3.652 5.232 -0.248) 6.385 1.927 ( 4.155 2.111 1.807) 4.999 1.949 ( 4.356 2.098 -1.458) 5.050 2.318 ( -2.536 -8.471 1.084) 8.909 2.357 ( -2.664 -8.630 -2.034) 9.258 2.396 ( -2.511 -10.353 0.940) 10.694 2.937 ( -2.884 4.503 0.003) 5.347 2.946 ( -1.842 3.917 0.075) 4.330 2.954 ( -1.392 3.498 -0.056) 3.765 3.094 ( 0.831 5.051 0.109) 5.120 3.105 ( 0.874 4.696 -0.283) 4.785 3.125 ( 1.328 4.586 0.315) 4.785 3.167 ( 3.474 0.226 -0.794) 3.571 3.192 ( 3.061 -0.274 1.090) 3.260 3.220 ( 2.680 -0.608 -0.430) 2.782 3.873 ( 3.933 -1.857 -0.797) 4.422 3.945 ( 2.827 -3.591 2.523) 5.220 3.974 ( 2.320 -4.187 -1.739) 5.093 ======================= Grid point 253 (38/48) ======================= q-point: ( 0.43 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.094 ( 4.495 2.105 -0.322) 4.974 1.130 ( 3.412 0.542 0.908) 3.572 1.145 ( 2.196 -0.246 -0.587) 2.286 1.287 ( -3.533 8.901 -0.110) 9.577 1.304 ( -3.070 7.898 0.519) 8.490 1.310 ( -3.342 7.425 -0.413) 8.153 1.626 ( -4.381 7.041 1.022) 8.356 1.634 ( -4.431 6.591 -1.111) 8.020 1.709 ( -1.715 2.413 0.084) 2.962 1.974 ( 2.230 4.469 -0.347) 5.006 1.994 ( 1.595 1.507 0.833) 2.347 2.013 ( 0.764 1.782 -0.455) 1.992 2.204 ( -1.694 -5.230 0.360) 5.509 2.236 ( -1.290 -5.926 -1.178) 6.178 2.255 ( -2.543 -6.537 0.791) 7.059 2.941 ( -1.871 3.846 -0.018) 4.277 2.956 ( -1.131 2.733 0.046) 2.959 2.966 ( -0.997 2.608 -0.022) 2.792 3.143 ( -0.019 3.168 0.287) 3.181 3.149 ( -0.213 2.902 -0.295) 2.925 3.165 ( -0.599 2.669 -0.081) 2.737 3.225 ( 2.527 -0.638 -0.152) 2.610 3.235 ( 2.317 -1.313 0.595) 2.729 3.253 ( 1.910 -1.513 -0.349) 2.461 3.907 ( 2.584 -2.547 -0.332) 3.643 3.953 ( 2.500 -3.796 1.355) 4.742 3.966 ( 1.783 -4.205 -1.028) 4.681 ======================= Grid point 273 (39/48) ======================= q-point: ( 0.14 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.631 ( 4.084 7.074 0.000) 8.168 0.652 ( 4.209 7.290 0.000) 8.418 0.699 ( 1.595 2.763 0.000) 3.190 1.032 ( 7.266 12.586 0.000) 14.533 1.037 ( 7.313 12.666 0.000) 14.625 1.078 ( 6.666 11.546 0.000) 13.332 1.446 ( 5.268 9.124 0.000) 10.536 1.501 ( 6.320 10.946 0.000) 12.639 1.505 ( 6.308 10.926 -0.000) 12.616 1.841 ( 4.170 7.222 0.000) 8.339 1.877 ( 3.700 6.409 0.000) 7.401 1.908 ( 3.657 6.334 -0.000) 7.314 2.137 ( 0.166 0.288 -0.000) 0.332 2.154 ( 0.415 0.719 0.000) 0.830 2.192 ( 0.521 0.902 0.000) 1.042 2.870 ( 4.637 8.031 0.000) 9.273 2.894 ( 4.643 8.042 0.000) 9.286 2.904 ( 4.346 7.527 0.000) 8.691 2.981 ( 0.453 0.785 0.000) 0.907 3.021 ( -0.780 -1.352 0.000) 1.561 3.024 ( -0.768 -1.331 -0.000) 1.537 3.432 ( -4.787 -8.291 0.000) 9.573 3.438 ( -4.749 -8.226 -0.000) 9.498 3.568 ( -6.551 -11.346 0.000) 13.101 3.582 ( 3.264 5.653 -0.000) 6.528 3.821 ( -0.229 -0.397 0.000) 0.458 3.824 ( -0.236 -0.409 -0.000) 0.473 ======================= Grid point 275 (40/48) ======================= q-point: ( 0.21 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.769 ( 4.023 5.284 -0.415) 6.654 0.786 ( 5.373 5.443 0.757) 7.686 0.809 ( 6.176 5.598 -0.340) 8.342 1.248 ( 2.688 14.053 -0.370) 14.312 1.255 ( 2.606 14.088 0.435) 14.334 1.279 ( 2.480 12.932 -0.052) 13.168 1.608 ( 6.423 3.628 -0.270) 7.381 1.686 ( 6.422 4.208 1.993) 7.932 1.696 ( 7.295 4.012 -1.761) 8.510 1.907 ( -7.150 12.108 0.112) 14.061 1.924 ( -7.355 10.590 -0.562) 12.906 1.949 ( -6.292 9.640 0.495) 11.522 2.206 ( 8.535 -7.163 0.752) 11.168 2.213 ( 10.388 -7.110 -0.645) 12.605 2.235 ( 8.834 -8.112 -0.132) 11.994 2.986 ( -0.146 0.319 -0.081) 0.360 2.991 ( 0.164 0.784 -0.052) 0.803 2.996 ( 0.810 1.599 0.217) 1.806 3.017 ( 3.248 5.436 -0.074) 6.333 3.048 ( 2.554 5.458 -0.077) 6.026 3.054 ( 2.596 5.295 0.063) 5.898 3.272 ( -6.626 -1.475 0.834) 6.839 3.274 ( -7.112 -2.038 -0.989) 7.464 3.350 ( -9.583 -1.570 0.261) 9.714 3.692 ( 5.531 -0.106 -0.731) 5.580 3.824 ( 2.688 -3.176 3.560) 5.476 3.845 ( 4.005 -4.537 -2.927) 6.723 ======================= Grid point 277 (41/48) ======================= q-point: ( 0.29 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.898 ( 5.446 5.044 -0.696) 7.456 0.939 ( 6.285 4.636 1.572) 7.967 0.970 ( 6.590 4.382 -0.881) 7.963 1.393 ( -1.678 13.753 -0.264) 13.857 1.397 ( -1.562 13.169 0.312) 13.265 1.408 ( -1.774 12.480 -0.032) 12.605 1.717 ( 3.360 0.599 -0.478) 3.446 1.776 ( 0.096 1.539 1.825) 2.389 1.796 ( -0.474 2.190 -1.481) 2.686 1.910 ( -1.539 8.216 0.282) 8.364 1.922 ( -1.522 7.614 -0.208) 7.767 1.947 ( -2.601 7.620 0.048) 8.052 2.237 ( 2.849 -6.588 0.674) 7.210 2.273 ( 3.584 -6.000 -1.231) 7.096 2.297 ( 4.415 -5.727 0.545) 7.252 2.980 ( -0.557 0.388 0.282) 0.735 2.989 ( -0.606 0.406 -0.548) 0.912 3.002 ( -1.050 0.425 0.310) 1.175 3.096 ( 2.933 -0.558 -0.345) 3.006 3.116 ( 2.689 -0.327 0.563) 2.766 3.133 ( 2.905 0.495 -0.276) 2.960 3.204 ( -3.400 6.231 0.212) 7.101 3.213 ( -3.348 6.444 -0.289) 7.268 3.252 ( -5.231 7.090 0.127) 8.812 3.767 ( 5.125 -2.882 -0.705) 5.922 3.838 ( 3.820 -5.339 2.332) 6.966 3.863 ( 4.153 -6.113 -1.653) 7.573 ======================= Grid point 279 (42/48) ======================= q-point: ( 0.36 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.038 ( 5.212 3.943 -0.558) 6.560 1.073 ( 4.760 2.635 1.290) 5.592 1.099 ( 4.179 2.139 -0.742) 4.753 1.461 ( -3.962 10.161 0.056) 10.906 1.462 ( -4.070 10.584 -0.096) 11.341 1.470 ( -4.509 11.484 0.050) 12.337 1.739 ( -0.384 0.597 -0.263) 0.757 1.760 ( -3.109 2.910 0.445) 4.282 1.775 ( -3.583 3.650 -0.188) 5.118 1.991 ( 1.652 5.560 0.040) 5.800 1.997 ( 1.873 6.230 0.260) 6.510 2.021 ( 2.418 5.589 -0.353) 6.100 2.185 ( -2.289 -3.802 0.558) 4.473 2.225 ( -3.160 -3.169 -0.871) 4.560 2.257 ( -3.348 -2.942 0.365) 4.471 2.978 ( -0.095 0.618 0.014) 0.626 2.989 ( -0.038 0.656 -0.204) 0.688 2.993 ( -0.380 0.544 0.212) 0.697 3.137 ( 3.706 -2.156 -0.653) 4.337 3.153 ( 3.389 -2.368 1.027) 4.260 3.173 ( 3.040 -2.566 -0.374) 3.996 3.226 ( -2.367 6.385 0.166) 6.811 3.237 ( -2.259 6.298 -0.364) 6.701 3.263 ( -3.268 7.537 0.192) 8.217 3.813 ( 4.328 -4.054 -0.520) 5.953 3.847 ( 3.854 -6.100 1.298) 7.332 3.866 ( 3.508 -6.500 -0.790) 7.428 ======================= Grid point 281 (43/48) ======================= q-point: ( 0.43 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.118 ( -0.060 0.105 0.000) 0.121 1.132 ( 0.113 -0.195 0.000) 0.226 1.135 ( 0.192 -0.332 0.000) 0.383 1.474 ( -4.933 8.544 -0.000) 9.866 1.478 ( -5.048 8.744 0.000) 10.097 1.485 ( -5.827 10.093 0.000) 11.654 1.739 ( -0.365 0.633 0.000) 0.730 1.746 ( -2.208 3.824 0.000) 4.416 1.763 ( -2.407 4.169 0.000) 4.814 2.052 ( -2.796 4.843 0.000) 5.593 2.067 ( -2.810 4.867 0.000) 5.620 2.108 ( -3.943 6.830 -0.000) 7.886 2.115 ( 1.815 -3.143 0.000) 3.629 2.150 ( 0.136 -0.235 0.000) 0.271 2.184 ( 0.095 -0.164 0.000) 0.190 2.984 ( -0.221 0.382 0.000) 0.441 2.990 ( -0.309 0.535 0.000) 0.618 2.996 ( -0.158 0.273 0.000) 0.315 3.164 ( 2.075 -3.595 -0.000) 4.151 3.173 ( 2.133 -3.694 0.000) 4.266 3.183 ( 1.943 -3.365 0.000) 3.886 3.238 ( -3.105 5.379 0.000) 6.211 3.250 ( -2.739 4.744 0.000) 5.477 3.272 ( -3.579 6.200 0.000) 7.159 3.832 ( 2.830 -4.902 0.000) 5.661 3.850 ( 3.709 -6.425 0.000) 7.419 3.860 ( 3.590 -6.218 0.000) 7.180 ======================= Grid point 305 (44/48) ======================= q-point: ( 0.21 0.21 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.877 ( 3.537 6.126 0.000) 7.074 0.906 ( 3.794 6.571 0.000) 7.587 0.932 ( 3.842 6.655 0.000) 7.685 1.491 ( 5.755 9.968 0.000) 11.510 1.501 ( 5.666 9.814 0.000) 11.332 1.515 ( 6.025 10.435 0.000) 12.050 1.666 ( 1.503 2.603 0.000) 3.006 1.750 ( 1.385 2.400 0.000) 2.771 1.751 ( 1.285 2.225 -0.000) 2.570 2.035 ( -1.008 -1.747 0.000) 2.017 2.047 ( -1.611 -2.791 0.000) 3.222 2.077 ( -2.422 -4.194 0.000) 4.843 2.131 ( 3.906 6.766 0.000) 7.812 2.157 ( 3.534 6.122 0.000) 7.069 2.184 ( 3.465 6.002 0.000) 6.931 2.988 ( 0.027 0.047 0.000) 0.054 2.995 ( 0.095 0.165 0.000) 0.190 3.011 ( 0.169 0.293 0.000) 0.339 3.065 ( -0.075 -0.129 0.000) 0.149 3.088 ( -0.410 -0.710 0.000) 0.820 3.092 ( -0.738 -1.278 0.000) 1.476 3.288 ( 1.743 3.019 -0.000) 3.486 3.298 ( 1.762 3.053 0.000) 3.525 3.361 ( 1.202 2.082 0.000) 2.404 3.693 ( 0.195 0.339 0.000) 0.391 3.769 ( -1.549 -2.682 0.000) 3.097 3.769 ( -1.450 -2.511 -0.000) 2.900 ======================= Grid point 307 (45/48) ======================= q-point: ( 0.29 0.21 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.003 ( 4.209 5.553 -0.330) 6.975 1.031 ( 4.094 4.615 0.690) 6.207 1.058 ( 4.084 4.444 -0.364) 6.047 1.603 ( -2.177 7.207 0.343) 7.536 1.611 ( -2.072 7.371 -0.552) 7.677 1.642 ( -2.467 9.068 0.407) 9.406 1.714 ( 2.435 -0.767 -0.664) 2.638 1.778 ( 1.076 -0.619 1.799) 2.186 1.795 ( 1.668 -1.317 -1.339) 2.512 2.016 ( 0.791 2.058 0.006) 2.205 2.027 ( 0.911 3.899 0.211) 4.009 2.039 ( 0.986 3.289 -0.235) 3.441 2.230 ( -0.506 5.833 0.328) 5.864 2.265 ( 0.512 4.641 -0.604) 4.708 2.291 ( 1.093 4.122 0.305) 4.275 2.983 ( -0.599 -0.121 0.137) 0.626 2.992 ( -0.546 -0.156 -0.265) 0.627 3.000 ( -1.021 -0.621 0.155) 1.205 3.058 ( 1.830 -2.574 -0.502) 3.199 3.069 ( 1.909 -3.383 0.732) 3.953 3.085 ( 2.506 -3.518 -0.264) 4.327 3.338 ( -1.656 6.148 0.488) 6.386 3.348 ( -1.547 6.068 -0.609) 6.292 3.397 ( -2.404 6.562 0.139) 6.990 3.706 ( 2.782 -3.044 -0.330) 4.137 3.730 ( 1.923 -5.366 0.804) 5.756 3.743 ( 2.631 -5.752 -0.486) 6.344 ======================= Grid point 309 (46/48) ======================= q-point: ( 0.36 0.21 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.108 ( 2.393 2.953 -0.156) 3.804 1.119 ( 2.041 2.129 0.322) 2.967 1.136 ( 1.051 1.605 -0.172) 1.926 1.610 ( -2.937 4.951 0.295) 5.764 1.616 ( -3.262 4.881 -0.454) 5.888 1.645 ( -3.832 4.855 0.282) 6.191 1.727 ( 0.880 -1.356 -0.348) 1.654 1.772 ( 0.016 -1.198 0.417) 1.269 1.785 ( 0.156 -1.649 -0.172) 1.665 2.087 ( 1.430 4.865 0.010) 5.071 2.117 ( 1.485 6.807 0.486) 6.984 2.139 ( 1.912 7.284 -0.702) 7.563 2.239 ( -4.636 7.040 0.698) 8.458 2.261 ( -5.616 5.208 -0.826) 7.704 2.286 ( -5.897 4.464 0.325) 7.403 2.974 ( 0.376 -0.885 -0.006) 0.961 2.983 ( 0.210 -1.127 -0.007) 1.147 2.987 ( 0.474 -0.840 0.019) 0.965 3.072 ( 3.345 -4.399 -0.482) 5.547 3.077 ( 3.199 -4.755 0.537) 5.756 3.093 ( 2.982 -4.662 -0.068) 5.535 3.356 ( -3.280 6.356 0.228) 7.156 3.366 ( -3.276 6.395 -0.294) 7.191 3.407 ( -3.900 6.850 0.070) 7.883 3.718 ( 3.735 -6.224 -0.091) 7.259 3.724 ( 3.587 -5.698 0.125) 6.734 3.738 ( 3.351 -5.862 -0.036) 6.752 ======================= Grid point 335 (47/48) ======================= q-point: ( 0.29 0.29 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.102 ( 2.342 4.056 0.000) 4.684 1.111 ( 1.955 3.387 0.000) 3.910 1.133 ( 1.713 2.967 0.000) 3.426 1.671 ( 0.323 0.559 -0.000) 0.646 1.675 ( 0.405 0.701 0.000) 0.809 1.695 ( -0.268 -0.465 0.000) 0.536 1.708 ( -0.622 -1.077 0.000) 1.244 1.758 ( -1.026 -1.777 -0.000) 2.052 1.761 ( -0.947 -1.640 0.000) 1.894 2.111 ( 4.345 7.526 0.000) 8.691 2.144 ( 4.519 7.827 0.000) 9.038 2.162 ( 4.764 8.252 -0.000) 9.529 2.323 ( 1.154 1.999 0.000) 2.308 2.352 ( 1.345 2.330 0.000) 2.691 2.369 ( 1.223 2.118 0.000) 2.446 2.970 ( -0.779 -1.348 0.000) 1.557 2.976 ( -0.851 -1.474 0.000) 1.703 2.980 ( -0.674 -1.167 0.000) 1.348 3.005 ( -1.278 -2.213 0.000) 2.555 3.013 ( -0.958 -1.659 0.000) 1.916 3.034 ( -0.793 -1.373 0.000) 1.585 3.433 ( 2.026 3.509 -0.000) 4.052 3.442 ( 2.055 3.559 0.000) 4.110 3.495 ( 1.910 3.308 0.000) 3.820 3.639 ( -2.138 -3.703 0.000) 4.276 3.647 ( -1.899 -3.289 0.000) 3.798 3.662 ( -1.136 -1.967 0.000) 2.271 ======================= Grid point 337 (48/48) ======================= q-point: ( 0.36 0.29 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.72e-04 2.72e-04 2.72e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 105 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.141 ( -0.409 0.709 0.000) 0.819 1.154 ( -0.474 0.821 0.000) 0.949 1.161 ( -0.563 0.975 0.000) 1.126 1.673 ( -1.190 2.061 0.000) 2.379 1.678 ( -0.895 1.550 0.000) 1.790 1.681 ( -0.647 1.121 0.000) 1.294 1.701 ( 0.567 -0.982 -0.000) 1.134 1.732 ( 1.321 -2.287 0.000) 2.641 1.739 ( 1.454 -2.518 0.000) 2.908 2.216 ( -3.691 6.393 0.000) 7.382 2.256 ( -3.267 5.659 0.000) 6.535 2.282 ( -3.321 5.752 0.000) 6.642 2.320 ( -0.247 0.427 0.000) 0.493 2.349 ( -1.233 2.136 0.000) 2.467 2.360 ( -1.151 1.993 0.000) 2.301 2.951 ( 0.675 -1.170 0.000) 1.351 2.960 ( 0.528 -0.915 0.000) 1.057 2.965 ( 0.593 -1.027 0.000) 1.186 2.985 ( 1.653 -2.863 -0.000) 3.306 2.998 ( 1.559 -2.701 0.000) 3.118 3.022 ( 1.220 -2.112 0.000) 2.439 3.470 ( -3.029 5.246 -0.000) 6.057 3.481 ( -2.967 5.139 0.000) 5.934 3.524 ( -2.818 4.881 0.000) 5.636 3.599 ( 3.152 -5.460 0.000) 6.305 3.612 ( 2.971 -5.147 0.000) 5.943 3.643 ( 2.440 -4.226 0.000) 4.879 =================== 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/10584 10.0 71.703 71.703 4.587 -0.000 -0.000 0.000 3/10584 20.0 29.965 29.965 1.953 -0.000 -0.000 0.000 3/10584 30.0 15.636 15.636 1.097 -0.000 -0.000 0.000 3/10584 40.0 10.039 10.039 0.746 -0.000 -0.000 0.000 3/10584 50.0 7.334 7.334 0.567 -0.000 -0.000 0.000 3/10584 60.0 5.785 5.785 0.459 -0.000 -0.000 0.000 3/10584 70.0 4.788 4.788 0.386 -0.000 -0.000 0.000 3/10584 80.0 4.092 4.092 0.335 -0.000 -0.000 0.000 3/10584 90.0 3.578 3.578 0.295 -0.000 -0.000 0.000 3/10584 100.0 3.183 3.183 0.264 -0.000 -0.000 0.000 3/10584 110.0 2.868 2.868 0.240 -0.000 -0.000 0.000 3/10584 120.0 2.612 2.612 0.219 -0.000 -0.000 0.000 3/10584 130.0 2.398 2.398 0.202 -0.000 -0.000 0.000 3/10584 140.0 2.218 2.218 0.187 -0.000 -0.000 0.000 3/10584 150.0 2.063 2.063 0.174 -0.000 -0.000 0.000 3/10584 160.0 1.929 1.929 0.163 -0.000 -0.000 0.000 3/10584 170.0 1.811 1.811 0.154 -0.000 -0.000 0.000 3/10584 180.0 1.707 1.707 0.145 -0.000 -0.000 0.000 3/10584 190.0 1.615 1.615 0.137 -0.000 -0.000 0.000 3/10584 200.0 1.532 1.532 0.130 -0.000 -0.000 0.000 3/10584 210.0 1.457 1.457 0.124 -0.000 -0.000 0.000 3/10584 220.0 1.390 1.390 0.119 -0.000 -0.000 0.000 3/10584 230.0 1.328 1.328 0.113 -0.000 -0.000 0.000 3/10584 240.0 1.272 1.272 0.109 -0.000 -0.000 0.000 3/10584 250.0 1.220 1.220 0.104 -0.000 -0.000 0.000 3/10584 260.0 1.173 1.173 0.100 -0.000 -0.000 0.000 3/10584 270.0 1.129 1.129 0.096 -0.000 -0.000 0.000 3/10584 280.0 1.088 1.088 0.093 -0.000 -0.000 0.000 3/10584 290.0 1.050 1.050 0.090 -0.000 -0.000 0.000 3/10584 300.0 1.014 1.014 0.087 -0.000 -0.000 0.000 3/10584 310.0 0.981 0.981 0.084 -0.000 -0.000 0.000 3/10584 320.0 0.950 0.950 0.081 -0.000 -0.000 0.000 3/10584 330.0 0.921 0.921 0.079 -0.000 -0.000 0.000 3/10584 340.0 0.894 0.894 0.077 -0.000 -0.000 0.000 3/10584 350.0 0.868 0.868 0.074 -0.000 -0.000 0.000 3/10584 360.0 0.844 0.844 0.072 -0.000 -0.000 0.000 3/10584 370.0 0.821 0.821 0.070 -0.000 -0.000 0.000 3/10584 380.0 0.799 0.799 0.068 -0.000 -0.000 0.000 3/10584 390.0 0.779 0.779 0.067 -0.000 -0.000 0.000 3/10584 400.0 0.759 0.759 0.065 -0.000 -0.000 0.000 3/10584 410.0 0.740 0.740 0.063 -0.000 -0.000 0.000 3/10584 420.0 0.723 0.723 0.062 -0.000 -0.000 0.000 3/10584 430.0 0.706 0.706 0.061 -0.000 -0.000 0.000 3/10584 440.0 0.690 0.690 0.059 -0.000 -0.000 0.000 3/10584 450.0 0.674 0.674 0.058 -0.000 -0.000 0.000 3/10584 460.0 0.660 0.660 0.057 -0.000 -0.000 0.000 3/10584 470.0 0.645 0.645 0.055 -0.000 -0.000 0.000 3/10584 480.0 0.632 0.632 0.054 -0.000 -0.000 0.000 3/10584 490.0 0.619 0.619 0.053 -0.000 -0.000 0.000 3/10584 500.0 0.607 0.607 0.052 -0.000 -0.000 0.000 3/10584 510.0 0.595 0.595 0.051 -0.000 -0.000 0.000 3/10584 520.0 0.583 0.583 0.050 -0.000 -0.000 0.000 3/10584 530.0 0.572 0.572 0.049 -0.000 -0.000 0.000 3/10584 540.0 0.562 0.562 0.048 -0.000 -0.000 0.000 3/10584 550.0 0.551 0.551 0.047 -0.000 -0.000 0.000 3/10584 560.0 0.541 0.541 0.046 -0.000 -0.000 0.000 3/10584 570.0 0.532 0.532 0.046 -0.000 -0.000 0.000 3/10584 580.0 0.523 0.523 0.045 -0.000 -0.000 0.000 3/10584 590.0 0.514 0.514 0.044 -0.000 -0.000 0.000 3/10584 600.0 0.505 0.505 0.043 -0.000 -0.000 0.000 3/10584 610.0 0.497 0.497 0.043 -0.000 -0.000 0.000 3/10584 620.0 0.489 0.489 0.042 -0.000 -0.000 0.000 3/10584 630.0 0.481 0.481 0.041 -0.000 -0.000 0.000 3/10584 640.0 0.474 0.474 0.041 -0.000 -0.000 0.000 3/10584 650.0 0.466 0.466 0.040 -0.000 -0.000 0.000 3/10584 660.0 0.459 0.459 0.039 -0.000 -0.000 0.000 3/10584 670.0 0.452 0.452 0.039 -0.000 -0.000 0.000 3/10584 680.0 0.446 0.446 0.038 -0.000 -0.000 0.000 3/10584 690.0 0.439 0.439 0.038 -0.000 -0.000 0.000 3/10584 700.0 0.433 0.433 0.037 -0.000 -0.000 0.000 3/10584 710.0 0.427 0.427 0.037 -0.000 -0.000 0.000 3/10584 720.0 0.421 0.421 0.036 -0.000 -0.000 0.000 3/10584 730.0 0.415 0.415 0.036 -0.000 -0.000 0.000 3/10584 740.0 0.409 0.409 0.035 -0.000 -0.000 0.000 3/10584 750.0 0.404 0.404 0.035 -0.000 -0.000 0.000 3/10584 760.0 0.399 0.399 0.034 -0.000 -0.000 0.000 3/10584 770.0 0.394 0.394 0.034 -0.000 -0.000 0.000 3/10584 780.0 0.388 0.388 0.033 -0.000 -0.000 0.000 3/10584 790.0 0.384 0.384 0.033 -0.000 -0.000 0.000 3/10584 800.0 0.379 0.379 0.033 -0.000 -0.000 0.000 3/10584 810.0 0.374 0.374 0.032 -0.000 -0.000 0.000 3/10584 820.0 0.369 0.369 0.032 -0.000 -0.000 0.000 3/10584 830.0 0.365 0.365 0.031 -0.000 -0.000 0.000 3/10584 840.0 0.361 0.361 0.031 -0.000 -0.000 0.000 3/10584 850.0 0.356 0.356 0.031 -0.000 -0.000 0.000 3/10584 860.0 0.352 0.352 0.030 -0.000 -0.000 0.000 3/10584 870.0 0.348 0.348 0.030 -0.000 -0.000 0.000 3/10584 880.0 0.344 0.344 0.030 -0.000 -0.000 0.000 3/10584 890.0 0.340 0.340 0.029 -0.000 -0.000 0.000 3/10584 900.0 0.337 0.337 0.029 -0.000 -0.000 0.000 3/10584 910.0 0.333 0.333 0.029 -0.000 -0.000 0.000 3/10584 920.0 0.329 0.329 0.028 -0.000 -0.000 0.000 3/10584 930.0 0.326 0.326 0.028 -0.000 -0.000 0.000 3/10584 940.0 0.322 0.322 0.028 -0.000 -0.000 0.000 3/10584 950.0 0.319 0.319 0.027 -0.000 -0.000 0.000 3/10584 960.0 0.316 0.316 0.027 -0.000 -0.000 0.000 3/10584 970.0 0.312 0.312 0.027 -0.000 -0.000 0.000 3/10584 980.0 0.309 0.309 0.027 -0.000 -0.000 0.000 3/10584 990.0 0.306 0.306 0.026 -0.000 -0.000 0.000 3/10584 1000.0 0.303 0.303 0.026 -0.000 -0.000 0.000 3/10584 Thermal conductivity related properties were written into "kappa-m14142.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 18:57:02]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|