----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 17:37:03]------------------------- Compiled with OpenMP support (max 32 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [4 4 1] Primitive matrix: [ 1. -0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) Number of symmetry operations in supercell: 192 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.222590011931395 0.000000000000000 0.000000000000000 b -2.111295005965697 3.656870220099024 0.000000000000000 c 0.000000000000000 0.000000000000000 20.593476540000001 Atomic positions (fractional): *1 Cd 0.66666666666667 0.33333333333333 0.33272442880812 112.411 2 Cd 0.33333333333333 0.66666666666667 0.66727557119188 112.411 *3 Cd 0.00000000000000 0.00000000000000 0.00000000000000 112.411 *4 I 0.33333333333333 0.66666666666667 0.08276777243055 126.904 *5 I 0.33333333333333 0.66666666666667 0.41547370484215 126.904 6 I 0.66666666666667 0.33333333333333 0.91723222756945 126.904 7 I 0.66666666666667 0.33333333333333 0.58452629515785 126.904 *8 I 0.00000000000000 0.00000000000000 0.24996113425933 126.904 9 I 0.00000000000000 0.00000000000000 0.75003886574067 126.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.222590011931395 0.000000000000000 0.000000000000000 b -2.111295005965697 3.656870220099024 0.000000000000000 c 0.000000000000000 0.000000000000000 20.593476540000001 Atomic positions (fractional): *1 Cd 0.66666666666667 0.33333333333333 0.33272442880812 112.411 > 1 2 Cd 0.33333333333333 0.66666666666667 0.66727557119188 112.411 > 2 *3 Cd 0.00000000000000 0.00000000000000 0.00000000000000 112.411 > 3 *4 I 0.33333333333333 0.66666666666667 0.08276777243055 126.904 > 4 *5 I 0.33333333333333 0.66666666666667 0.41547370484215 126.904 > 5 6 I 0.66666666666667 0.33333333333333 0.91723222756945 126.904 > 6 7 I 0.66666666666667 0.33333333333333 0.58452629515785 126.904 > 7 *8 I 0.00000000000000 0.00000000000000 0.24996113425933 126.904 > 8 9 I 0.00000000000000 0.00000000000000 0.75003886574067 126.904 > 9 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.890360047725579 0.000000000000000 0.000000000000000 b -8.445180023862788 14.627480880396096 0.000000000000000 c 0.000000000000000 0.000000000000000 20.593476540000001 Atomic positions (fractional): *1 Cd 0.16666666666667 0.08333333333333 0.33272442880812 112.411 > 1 2 Cd 0.41666666666667 0.08333333333333 0.33272442880812 112.411 > 1 3 Cd 0.66666666666667 0.08333333333333 0.33272442880812 112.411 > 1 4 Cd 0.91666666666667 0.08333333333333 0.33272442880812 112.411 > 1 5 Cd 0.16666666666667 0.33333333333333 0.33272442880812 112.411 > 1 6 Cd 0.41666666666667 0.33333333333333 0.33272442880812 112.411 > 1 7 Cd 0.66666666666667 0.33333333333333 0.33272442880812 112.411 > 1 8 Cd 0.91666666666667 0.33333333333333 0.33272442880812 112.411 > 1 9 Cd 0.16666666666667 0.58333333333333 0.33272442880812 112.411 > 1 10 Cd 0.41666666666667 0.58333333333333 0.33272442880812 112.411 > 1 11 Cd 0.66666666666667 0.58333333333333 0.33272442880812 112.411 > 1 12 Cd 0.91666666666667 0.58333333333333 0.33272442880812 112.411 > 1 13 Cd 0.16666666666667 0.83333333333333 0.33272442880812 112.411 > 1 14 Cd 0.41666666666667 0.83333333333333 0.33272442880812 112.411 > 1 15 Cd 0.66666666666667 0.83333333333333 0.33272442880812 112.411 > 1 16 Cd 0.91666666666667 0.83333333333333 0.33272442880812 112.411 > 1 17 Cd 0.08333333333333 0.16666666666667 0.66727557119188 112.411 > 2 18 Cd 0.33333333333333 0.16666666666667 0.66727557119188 112.411 > 2 19 Cd 0.58333333333333 0.16666666666667 0.66727557119188 112.411 > 2 20 Cd 0.83333333333333 0.16666666666667 0.66727557119188 112.411 > 2 21 Cd 0.08333333333333 0.41666666666667 0.66727557119188 112.411 > 2 22 Cd 0.33333333333333 0.41666666666667 0.66727557119188 112.411 > 2 23 Cd 0.58333333333333 0.41666666666667 0.66727557119188 112.411 > 2 24 Cd 0.83333333333333 0.41666666666667 0.66727557119188 112.411 > 2 25 Cd 0.08333333333333 0.66666666666667 0.66727557119188 112.411 > 2 26 Cd 0.33333333333333 0.66666666666667 0.66727557119188 112.411 > 2 27 Cd 0.58333333333333 0.66666666666667 0.66727557119188 112.411 > 2 28 Cd 0.83333333333333 0.66666666666667 0.66727557119188 112.411 > 2 29 Cd 0.08333333333333 0.91666666666667 0.66727557119188 112.411 > 2 30 Cd 0.33333333333333 0.91666666666667 0.66727557119188 112.411 > 2 31 Cd 0.58333333333333 0.91666666666667 0.66727557119188 112.411 > 2 32 Cd 0.83333333333333 0.91666666666667 0.66727557119188 112.411 > 2 *33 Cd 0.00000000000000 0.00000000000000 0.00000000000000 112.411 > 3 34 Cd 0.25000000000000 0.00000000000000 0.00000000000000 112.411 > 3 35 Cd 0.50000000000000 0.00000000000000 0.00000000000000 112.411 > 3 36 Cd 0.75000000000000 0.00000000000000 0.00000000000000 112.411 > 3 37 Cd 0.00000000000000 0.25000000000000 0.00000000000000 112.411 > 3 38 Cd 0.25000000000000 0.25000000000000 0.00000000000000 112.411 > 3 39 Cd 0.50000000000000 0.25000000000000 0.00000000000000 112.411 > 3 40 Cd 0.75000000000000 0.25000000000000 0.00000000000000 112.411 > 3 41 Cd 0.00000000000000 0.50000000000000 0.00000000000000 112.411 > 3 42 Cd 0.25000000000000 0.50000000000000 0.00000000000000 112.411 > 3 43 Cd 0.50000000000000 0.50000000000000 0.00000000000000 112.411 > 3 44 Cd 0.75000000000000 0.50000000000000 0.00000000000000 112.411 > 3 45 Cd 0.00000000000000 0.75000000000000 0.00000000000000 112.411 > 3 46 Cd 0.25000000000000 0.75000000000000 0.00000000000000 112.411 > 3 47 Cd 0.50000000000000 0.75000000000000 0.00000000000000 112.411 > 3 48 Cd 0.75000000000000 0.75000000000000 0.00000000000000 112.411 > 3 *49 I 0.08333333333333 0.16666666666667 0.08276777243055 126.904 > 4 50 I 0.33333333333333 0.16666666666667 0.08276777243055 126.904 > 4 51 I 0.58333333333333 0.16666666666667 0.08276777243055 126.904 > 4 52 I 0.83333333333333 0.16666666666667 0.08276777243055 126.904 > 4 53 I 0.08333333333333 0.41666666666667 0.08276777243055 126.904 > 4 54 I 0.33333333333333 0.41666666666667 0.08276777243055 126.904 > 4 55 I 0.58333333333333 0.41666666666667 0.08276777243055 126.904 > 4 56 I 0.83333333333333 0.41666666666667 0.08276777243055 126.904 > 4 57 I 0.08333333333333 0.66666666666667 0.08276777243055 126.904 > 4 58 I 0.33333333333333 0.66666666666667 0.08276777243055 126.904 > 4 59 I 0.58333333333333 0.66666666666667 0.08276777243055 126.904 > 4 60 I 0.83333333333333 0.66666666666667 0.08276777243055 126.904 > 4 61 I 0.08333333333333 0.91666666666667 0.08276777243055 126.904 > 4 62 I 0.33333333333333 0.91666666666667 0.08276777243055 126.904 > 4 63 I 0.58333333333333 0.91666666666667 0.08276777243055 126.904 > 4 64 I 0.83333333333333 0.91666666666667 0.08276777243055 126.904 > 4 *65 I 0.08333333333333 0.16666666666667 0.41547370484215 126.904 > 5 66 I 0.33333333333333 0.16666666666667 0.41547370484215 126.904 > 5 67 I 0.58333333333333 0.16666666666667 0.41547370484215 126.904 > 5 68 I 0.83333333333333 0.16666666666667 0.41547370484215 126.904 > 5 69 I 0.08333333333333 0.41666666666667 0.41547370484215 126.904 > 5 70 I 0.33333333333333 0.41666666666667 0.41547370484215 126.904 > 5 71 I 0.58333333333333 0.41666666666667 0.41547370484215 126.904 > 5 72 I 0.83333333333333 0.41666666666667 0.41547370484215 126.904 > 5 73 I 0.08333333333333 0.66666666666667 0.41547370484215 126.904 > 5 74 I 0.33333333333333 0.66666666666667 0.41547370484215 126.904 > 5 75 I 0.58333333333333 0.66666666666667 0.41547370484215 126.904 > 5 76 I 0.83333333333333 0.66666666666667 0.41547370484215 126.904 > 5 77 I 0.08333333333333 0.91666666666667 0.41547370484215 126.904 > 5 78 I 0.33333333333333 0.91666666666667 0.41547370484215 126.904 > 5 79 I 0.58333333333333 0.91666666666667 0.41547370484215 126.904 > 5 80 I 0.83333333333333 0.91666666666667 0.41547370484215 126.904 > 5 81 I 0.16666666666667 0.08333333333333 0.91723222756945 126.904 > 6 82 I 0.41666666666667 0.08333333333333 0.91723222756945 126.904 > 6 83 I 0.66666666666667 0.08333333333333 0.91723222756945 126.904 > 6 84 I 0.91666666666667 0.08333333333333 0.91723222756945 126.904 > 6 85 I 0.16666666666667 0.33333333333333 0.91723222756945 126.904 > 6 86 I 0.41666666666667 0.33333333333333 0.91723222756945 126.904 > 6 87 I 0.66666666666667 0.33333333333333 0.91723222756945 126.904 > 6 88 I 0.91666666666667 0.33333333333333 0.91723222756945 126.904 > 6 89 I 0.16666666666667 0.58333333333333 0.91723222756945 126.904 > 6 90 I 0.41666666666667 0.58333333333333 0.91723222756945 126.904 > 6 91 I 0.66666666666667 0.58333333333333 0.91723222756945 126.904 > 6 92 I 0.91666666666667 0.58333333333333 0.91723222756945 126.904 > 6 93 I 0.16666666666667 0.83333333333333 0.91723222756945 126.904 > 6 94 I 0.41666666666667 0.83333333333333 0.91723222756945 126.904 > 6 95 I 0.66666666666667 0.83333333333333 0.91723222756945 126.904 > 6 96 I 0.91666666666667 0.83333333333333 0.91723222756945 126.904 > 6 97 I 0.16666666666667 0.08333333333333 0.58452629515785 126.904 > 7 98 I 0.41666666666667 0.08333333333333 0.58452629515785 126.904 > 7 99 I 0.66666666666667 0.08333333333333 0.58452629515785 126.904 > 7 100 I 0.91666666666667 0.08333333333333 0.58452629515785 126.904 > 7 101 I 0.16666666666667 0.33333333333333 0.58452629515785 126.904 > 7 102 I 0.41666666666667 0.33333333333333 0.58452629515785 126.904 > 7 103 I 0.66666666666667 0.33333333333333 0.58452629515785 126.904 > 7 104 I 0.91666666666667 0.33333333333333 0.58452629515785 126.904 > 7 105 I 0.16666666666667 0.58333333333333 0.58452629515785 126.904 > 7 106 I 0.41666666666667 0.58333333333333 0.58452629515785 126.904 > 7 107 I 0.66666666666667 0.58333333333333 0.58452629515785 126.904 > 7 108 I 0.91666666666667 0.58333333333333 0.58452629515785 126.904 > 7 109 I 0.16666666666667 0.83333333333333 0.58452629515785 126.904 > 7 110 I 0.41666666666667 0.83333333333333 0.58452629515785 126.904 > 7 111 I 0.66666666666667 0.83333333333333 0.58452629515785 126.904 > 7 112 I 0.91666666666667 0.83333333333333 0.58452629515785 126.904 > 7 *113 I 0.00000000000000 0.00000000000000 0.24996113425933 126.904 > 8 114 I 0.25000000000000 0.00000000000000 0.24996113425933 126.904 > 8 115 I 0.50000000000000 0.00000000000000 0.24996113425933 126.904 > 8 116 I 0.75000000000000 0.00000000000000 0.24996113425933 126.904 > 8 117 I 0.00000000000000 0.25000000000000 0.24996113425933 126.904 > 8 118 I 0.25000000000000 0.25000000000000 0.24996113425933 126.904 > 8 119 I 0.50000000000000 0.25000000000000 0.24996113425933 126.904 > 8 120 I 0.75000000000000 0.25000000000000 0.24996113425933 126.904 > 8 121 I 0.00000000000000 0.50000000000000 0.24996113425933 126.904 > 8 122 I 0.25000000000000 0.50000000000000 0.24996113425933 126.904 > 8 123 I 0.50000000000000 0.50000000000000 0.24996113425933 126.904 > 8 124 I 0.75000000000000 0.50000000000000 0.24996113425933 126.904 > 8 125 I 0.00000000000000 0.75000000000000 0.24996113425933 126.904 > 8 126 I 0.25000000000000 0.75000000000000 0.24996113425933 126.904 > 8 127 I 0.50000000000000 0.75000000000000 0.24996113425933 126.904 > 8 128 I 0.75000000000000 0.75000000000000 0.24996113425933 126.904 > 8 129 I 0.00000000000000 0.00000000000000 0.75003886574067 126.904 > 9 130 I 0.25000000000000 0.00000000000000 0.75003886574067 126.904 > 9 131 I 0.50000000000000 0.00000000000000 0.75003886574067 126.904 > 9 132 I 0.75000000000000 0.00000000000000 0.75003886574067 126.904 > 9 133 I 0.00000000000000 0.25000000000000 0.75003886574067 126.904 > 9 134 I 0.25000000000000 0.25000000000000 0.75003886574067 126.904 > 9 135 I 0.50000000000000 0.25000000000000 0.75003886574067 126.904 > 9 136 I 0.75000000000000 0.25000000000000 0.75003886574067 126.904 > 9 137 I 0.00000000000000 0.50000000000000 0.75003886574067 126.904 > 9 138 I 0.25000000000000 0.50000000000000 0.75003886574067 126.904 > 9 139 I 0.50000000000000 0.50000000000000 0.75003886574067 126.904 > 9 140 I 0.75000000000000 0.50000000000000 0.75003886574067 126.904 > 9 141 I 0.00000000000000 0.75000000000000 0.75003886574067 126.904 > 9 142 I 0.25000000000000 0.75000000000000 0.75003886574067 126.904 > 9 143 I 0.50000000000000 0.75000000000000 0.75003886574067 126.904 > 9 144 I 0.75000000000000 0.75000000000000 0.75003886574067 126.904 > 9 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 6.0028127 -0.0000000 0.0000000 0.0000000 6.0028127 0.0000000 0.0000000 0.0000000 4.8454435 -------------------------- Born effective charges -------------------------- 1 Cd 3.0166357 -0.0000000 0.0000000 0.0000000 3.0166357 0.0000000 0.0000000 0.0000000 1.6823478 2 Cd 3.0166357 -0.0000000 0.0000000 0.0000000 3.0166357 0.0000000 0.0000000 0.0000000 1.6823478 3 Cd 3.0061572 -0.0000000 0.0000000 0.0000000 3.0061572 0.0000000 0.0000000 0.0000000 1.6589443 4 I -1.5130415 0.0000000 0.0000000 0.0000000 -1.5130415 0.0000000 0.0000000 0.0000000 -0.8380966 5 I -1.4892058 0.0000000 0.0000000 0.0000000 -1.4892058 0.0000000 0.0000000 0.0000000 -0.8488195 6 I -1.5130415 0.0000000 0.0000000 0.0000000 -1.5130415 0.0000000 0.0000000 0.0000000 -0.8380966 7 I -1.4892058 0.0000000 0.0000000 0.0000000 -1.4892058 0.0000000 0.0000000 0.0000000 -0.8488195 8 I -1.5174671 0.0000000 0.0000000 0.0000000 -1.5174671 0.0000000 0.0000000 0.0000000 -0.8249038 9 I -1.5174671 0.0000000 0.0000000 0.0000000 -1.5174671 0.0000000 0.0000000 0.0000000 -0.8249038 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 432/432 Permutation basis: 6426/6426 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 578 Number of blocks in projector: 444 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 255 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 203 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 120 Use standard eigh solver. Tree of FC basis block matrices: - (578, 568), data: False |-- (120, 115), data: True |-- (203, 203), data: True |-- (255, 250), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.017 Solver_block: 100 / 160 - Time: 0.255 Solver_block: 160 / 160 - Time: 0.183 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.459 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) Permutation basis: 432/432 Permutation basis: 6426/6426 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 578 Number of blocks in projector: 444 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 255 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 203 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 120 Use standard eigh solver. Tree of FC basis block matrices: - (578, 568), data: False |-- (120, 115), data: True |-- (203, 203), data: True |-- (255, 250), data: True Max drift after symmetrization by symfc projector: 0.00000000 (zz) 0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 17:37:14]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 17:37:15]------------------------- Compiled with OpenMP support (max 32 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 1] Primitive matrix: [ 1. -0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.222590011931395 0.000000000000000 0.000000000000000 b -2.111295005965697 3.656870220099024 0.000000000000000 c 0.000000000000000 0.000000000000000 20.593476540000001 Atomic positions (fractional): 1 Cd 0.66666666666667 0.33333333333333 0.33272442880812 112.411 2 Cd 0.33333333333333 0.66666666666667 0.66727557119188 112.411 3 Cd 0.00000000000000 0.00000000000000 0.00000000000000 112.411 4 I 0.33333333333333 0.66666666666667 0.08276777243055 126.904 5 I 0.33333333333333 0.66666666666667 0.41547370484215 126.904 6 I 0.66666666666667 0.33333333333333 0.91723222756945 126.904 7 I 0.66666666666667 0.33333333333333 0.58452629515785 126.904 8 I 0.00000000000000 0.00000000000000 0.24996113425933 126.904 9 I 0.00000000000000 0.00000000000000 0.75003886574067 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.890360047725579 0.000000000000000 0.000000000000000 b -8.445180023862788 14.627480880396096 0.000000000000000 c 0.000000000000000 0.000000000000000 20.593476540000001 Atomic positions (fractional): 1 Cd 0.16666666666667 0.08333333333333 0.33272442880812 112.411 > 1 2 Cd 0.41666666666667 0.08333333333333 0.33272442880812 112.411 > 1 3 Cd 0.66666666666667 0.08333333333333 0.33272442880812 112.411 > 1 4 Cd 0.91666666666667 0.08333333333333 0.33272442880812 112.411 > 1 5 Cd 0.16666666666667 0.33333333333333 0.33272442880812 112.411 > 1 6 Cd 0.41666666666667 0.33333333333333 0.33272442880812 112.411 > 1 7 Cd 0.66666666666667 0.33333333333333 0.33272442880812 112.411 > 1 8 Cd 0.91666666666667 0.33333333333333 0.33272442880812 112.411 > 1 9 Cd 0.16666666666667 0.58333333333333 0.33272442880812 112.411 > 1 10 Cd 0.41666666666667 0.58333333333333 0.33272442880812 112.411 > 1 11 Cd 0.66666666666667 0.58333333333333 0.33272442880812 112.411 > 1 12 Cd 0.91666666666667 0.58333333333333 0.33272442880812 112.411 > 1 13 Cd 0.16666666666667 0.83333333333333 0.33272442880812 112.411 > 1 14 Cd 0.41666666666667 0.83333333333333 0.33272442880812 112.411 > 1 15 Cd 0.66666666666667 0.83333333333333 0.33272442880812 112.411 > 1 16 Cd 0.91666666666667 0.83333333333333 0.33272442880812 112.411 > 1 17 Cd 0.08333333333333 0.16666666666667 0.66727557119188 112.411 > 17 18 Cd 0.33333333333333 0.16666666666667 0.66727557119188 112.411 > 17 19 Cd 0.58333333333333 0.16666666666667 0.66727557119188 112.411 > 17 20 Cd 0.83333333333333 0.16666666666667 0.66727557119188 112.411 > 17 21 Cd 0.08333333333333 0.41666666666667 0.66727557119188 112.411 > 17 22 Cd 0.33333333333333 0.41666666666667 0.66727557119188 112.411 > 17 23 Cd 0.58333333333333 0.41666666666667 0.66727557119188 112.411 > 17 24 Cd 0.83333333333333 0.41666666666667 0.66727557119188 112.411 > 17 25 Cd 0.08333333333333 0.66666666666667 0.66727557119188 112.411 > 17 26 Cd 0.33333333333333 0.66666666666667 0.66727557119188 112.411 > 17 27 Cd 0.58333333333333 0.66666666666667 0.66727557119188 112.411 > 17 28 Cd 0.83333333333333 0.66666666666667 0.66727557119188 112.411 > 17 29 Cd 0.08333333333333 0.91666666666667 0.66727557119188 112.411 > 17 30 Cd 0.33333333333333 0.91666666666667 0.66727557119188 112.411 > 17 31 Cd 0.58333333333333 0.91666666666667 0.66727557119188 112.411 > 17 32 Cd 0.83333333333333 0.91666666666667 0.66727557119188 112.411 > 17 33 Cd 0.00000000000000 0.00000000000000 0.00000000000000 112.411 > 33 34 Cd 0.25000000000000 0.00000000000000 0.00000000000000 112.411 > 33 35 Cd 0.50000000000000 0.00000000000000 0.00000000000000 112.411 > 33 36 Cd 0.75000000000000 0.00000000000000 0.00000000000000 112.411 > 33 37 Cd 0.00000000000000 0.25000000000000 0.00000000000000 112.411 > 33 38 Cd 0.25000000000000 0.25000000000000 0.00000000000000 112.411 > 33 39 Cd 0.50000000000000 0.25000000000000 0.00000000000000 112.411 > 33 40 Cd 0.75000000000000 0.25000000000000 0.00000000000000 112.411 > 33 41 Cd 0.00000000000000 0.50000000000000 0.00000000000000 112.411 > 33 42 Cd 0.25000000000000 0.50000000000000 0.00000000000000 112.411 > 33 43 Cd 0.50000000000000 0.50000000000000 0.00000000000000 112.411 > 33 44 Cd 0.75000000000000 0.50000000000000 0.00000000000000 112.411 > 33 45 Cd 0.00000000000000 0.75000000000000 0.00000000000000 112.411 > 33 46 Cd 0.25000000000000 0.75000000000000 0.00000000000000 112.411 > 33 47 Cd 0.50000000000000 0.75000000000000 0.00000000000000 112.411 > 33 48 Cd 0.75000000000000 0.75000000000000 0.00000000000000 112.411 > 33 49 I 0.08333333333333 0.16666666666667 0.08276777243055 126.904 > 49 50 I 0.33333333333333 0.16666666666667 0.08276777243055 126.904 > 49 51 I 0.58333333333333 0.16666666666667 0.08276777243055 126.904 > 49 52 I 0.83333333333333 0.16666666666667 0.08276777243055 126.904 > 49 53 I 0.08333333333333 0.41666666666667 0.08276777243055 126.904 > 49 54 I 0.33333333333333 0.41666666666667 0.08276777243055 126.904 > 49 55 I 0.58333333333333 0.41666666666667 0.08276777243055 126.904 > 49 56 I 0.83333333333333 0.41666666666667 0.08276777243055 126.904 > 49 57 I 0.08333333333333 0.66666666666667 0.08276777243055 126.904 > 49 58 I 0.33333333333333 0.66666666666667 0.08276777243055 126.904 > 49 59 I 0.58333333333333 0.66666666666667 0.08276777243055 126.904 > 49 60 I 0.83333333333333 0.66666666666667 0.08276777243055 126.904 > 49 61 I 0.08333333333333 0.91666666666667 0.08276777243055 126.904 > 49 62 I 0.33333333333333 0.91666666666667 0.08276777243055 126.904 > 49 63 I 0.58333333333333 0.91666666666667 0.08276777243055 126.904 > 49 64 I 0.83333333333333 0.91666666666667 0.08276777243055 126.904 > 49 65 I 0.08333333333333 0.16666666666667 0.41547370484215 126.904 > 65 66 I 0.33333333333333 0.16666666666667 0.41547370484215 126.904 > 65 67 I 0.58333333333333 0.16666666666667 0.41547370484215 126.904 > 65 68 I 0.83333333333333 0.16666666666667 0.41547370484215 126.904 > 65 69 I 0.08333333333333 0.41666666666667 0.41547370484215 126.904 > 65 70 I 0.33333333333333 0.41666666666667 0.41547370484215 126.904 > 65 71 I 0.58333333333333 0.41666666666667 0.41547370484215 126.904 > 65 72 I 0.83333333333333 0.41666666666667 0.41547370484215 126.904 > 65 73 I 0.08333333333333 0.66666666666667 0.41547370484215 126.904 > 65 74 I 0.33333333333333 0.66666666666667 0.41547370484215 126.904 > 65 75 I 0.58333333333333 0.66666666666667 0.41547370484215 126.904 > 65 76 I 0.83333333333333 0.66666666666667 0.41547370484215 126.904 > 65 77 I 0.08333333333333 0.91666666666667 0.41547370484215 126.904 > 65 78 I 0.33333333333333 0.91666666666667 0.41547370484215 126.904 > 65 79 I 0.58333333333333 0.91666666666667 0.41547370484215 126.904 > 65 80 I 0.83333333333333 0.91666666666667 0.41547370484215 126.904 > 65 81 I 0.16666666666667 0.08333333333333 0.91723222756945 126.904 > 81 82 I 0.41666666666667 0.08333333333333 0.91723222756945 126.904 > 81 83 I 0.66666666666667 0.08333333333333 0.91723222756945 126.904 > 81 84 I 0.91666666666667 0.08333333333333 0.91723222756945 126.904 > 81 85 I 0.16666666666667 0.33333333333333 0.91723222756945 126.904 > 81 86 I 0.41666666666667 0.33333333333333 0.91723222756945 126.904 > 81 87 I 0.66666666666667 0.33333333333333 0.91723222756945 126.904 > 81 88 I 0.91666666666667 0.33333333333333 0.91723222756945 126.904 > 81 89 I 0.16666666666667 0.58333333333333 0.91723222756945 126.904 > 81 90 I 0.41666666666667 0.58333333333333 0.91723222756945 126.904 > 81 91 I 0.66666666666667 0.58333333333333 0.91723222756945 126.904 > 81 92 I 0.91666666666667 0.58333333333333 0.91723222756945 126.904 > 81 93 I 0.16666666666667 0.83333333333333 0.91723222756945 126.904 > 81 94 I 0.41666666666667 0.83333333333333 0.91723222756945 126.904 > 81 95 I 0.66666666666667 0.83333333333333 0.91723222756945 126.904 > 81 96 I 0.91666666666667 0.83333333333333 0.91723222756945 126.904 > 81 97 I 0.16666666666667 0.08333333333333 0.58452629515785 126.904 > 97 98 I 0.41666666666667 0.08333333333333 0.58452629515785 126.904 > 97 99 I 0.66666666666667 0.08333333333333 0.58452629515785 126.904 > 97 100 I 0.91666666666667 0.08333333333333 0.58452629515785 126.904 > 97 101 I 0.16666666666667 0.33333333333333 0.58452629515785 126.904 > 97 102 I 0.41666666666667 0.33333333333333 0.58452629515785 126.904 > 97 103 I 0.66666666666667 0.33333333333333 0.58452629515785 126.904 > 97 104 I 0.91666666666667 0.33333333333333 0.58452629515785 126.904 > 97 105 I 0.16666666666667 0.58333333333333 0.58452629515785 126.904 > 97 106 I 0.41666666666667 0.58333333333333 0.58452629515785 126.904 > 97 107 I 0.66666666666667 0.58333333333333 0.58452629515785 126.904 > 97 108 I 0.91666666666667 0.58333333333333 0.58452629515785 126.904 > 97 109 I 0.16666666666667 0.83333333333333 0.58452629515785 126.904 > 97 110 I 0.41666666666667 0.83333333333333 0.58452629515785 126.904 > 97 111 I 0.66666666666667 0.83333333333333 0.58452629515785 126.904 > 97 112 I 0.91666666666667 0.83333333333333 0.58452629515785 126.904 > 97 113 I 0.00000000000000 0.00000000000000 0.24996113425933 126.904 > 113 114 I 0.25000000000000 0.00000000000000 0.24996113425933 126.904 > 113 115 I 0.50000000000000 0.00000000000000 0.24996113425933 126.904 > 113 116 I 0.75000000000000 0.00000000000000 0.24996113425933 126.904 > 113 117 I 0.00000000000000 0.25000000000000 0.24996113425933 126.904 > 113 118 I 0.25000000000000 0.25000000000000 0.24996113425933 126.904 > 113 119 I 0.50000000000000 0.25000000000000 0.24996113425933 126.904 > 113 120 I 0.75000000000000 0.25000000000000 0.24996113425933 126.904 > 113 121 I 0.00000000000000 0.50000000000000 0.24996113425933 126.904 > 113 122 I 0.25000000000000 0.50000000000000 0.24996113425933 126.904 > 113 123 I 0.50000000000000 0.50000000000000 0.24996113425933 126.904 > 113 124 I 0.75000000000000 0.50000000000000 0.24996113425933 126.904 > 113 125 I 0.00000000000000 0.75000000000000 0.24996113425933 126.904 > 113 126 I 0.25000000000000 0.75000000000000 0.24996113425933 126.904 > 113 127 I 0.50000000000000 0.75000000000000 0.24996113425933 126.904 > 113 128 I 0.75000000000000 0.75000000000000 0.24996113425933 126.904 > 113 129 I 0.00000000000000 0.00000000000000 0.75003886574067 126.904 > 129 130 I 0.25000000000000 0.00000000000000 0.75003886574067 126.904 > 129 131 I 0.50000000000000 0.00000000000000 0.75003886574067 126.904 > 129 132 I 0.75000000000000 0.00000000000000 0.75003886574067 126.904 > 129 133 I 0.00000000000000 0.25000000000000 0.75003886574067 126.904 > 129 134 I 0.25000000000000 0.25000000000000 0.75003886574067 126.904 > 129 135 I 0.50000000000000 0.25000000000000 0.75003886574067 126.904 > 129 136 I 0.75000000000000 0.25000000000000 0.75003886574067 126.904 > 129 137 I 0.00000000000000 0.50000000000000 0.75003886574067 126.904 > 129 138 I 0.25000000000000 0.50000000000000 0.75003886574067 126.904 > 129 139 I 0.50000000000000 0.50000000000000 0.75003886574067 126.904 > 129 140 I 0.75000000000000 0.50000000000000 0.75003886574067 126.904 > 129 141 I 0.00000000000000 0.75000000000000 0.75003886574067 126.904 > 129 142 I 0.25000000000000 0.75000000000000 0.75003886574067 126.904 > 129 143 I 0.50000000000000 0.75000000000000 0.75003886574067 126.904 > 129 144 I 0.75000000000000 0.75000000000000 0.75003886574067 126.904 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 6.0028127 -0.0000000 0.0000000 0.0000000 6.0028127 0.0000000 0.0000000 0.0000000 4.8454435 -------------------------- Born effective charges -------------------------- 1 Cd 3.0166357 -0.0000000 0.0000000 0.0000000 3.0166357 0.0000000 0.0000000 0.0000000 1.6823478 2 Cd 3.0166357 -0.0000000 0.0000000 0.0000000 3.0166357 0.0000000 0.0000000 0.0000000 1.6823478 3 Cd 3.0061572 -0.0000000 0.0000000 0.0000000 3.0061572 0.0000000 0.0000000 0.0000000 1.6589443 4 I -1.5130415 0.0000000 0.0000000 0.0000000 -1.5130415 0.0000000 0.0000000 0.0000000 -0.8380966 5 I -1.4892058 0.0000000 0.0000000 0.0000000 -1.4892058 0.0000000 0.0000000 0.0000000 -0.8488195 6 I -1.5130415 0.0000000 0.0000000 0.0000000 -1.5130415 0.0000000 0.0000000 0.0000000 -0.8380966 7 I -1.4892058 0.0000000 0.0000000 0.0000000 -1.4892058 0.0000000 0.0000000 0.0000000 -0.8488195 8 I -1.5174671 0.0000000 0.0000000 0.0000000 -1.5174671 0.0000000 0.0000000 0.0000000 -0.8249038 9 I -1.5174671 0.0000000 0.0000000 0.0000000 -1.5174671 0.0000000 0.0000000 0.0000000 -0.8249038 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 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[ 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] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000003 (yzy) 0.00000003 (yzy) 0.00000003 (yyz) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 17:37:30]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 17:37:30]------------------------- Compiled with OpenMP support (max 32 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 1] Primitive matrix: [ 1. -0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.222590011931395 0.000000000000000 0.000000000000000 b -2.111295005965697 3.656870220099024 0.000000000000000 c 0.000000000000000 0.000000000000000 20.593476540000001 Atomic positions (fractional): 1 Cd 0.66666666666667 0.33333333333333 0.33272442880812 112.411 2 Cd 0.33333333333333 0.66666666666667 0.66727557119188 112.411 3 Cd 0.00000000000000 0.00000000000000 0.00000000000000 112.411 4 I 0.33333333333333 0.66666666666667 0.08276777243055 126.904 5 I 0.33333333333333 0.66666666666667 0.41547370484215 126.904 6 I 0.66666666666667 0.33333333333333 0.91723222756945 126.904 7 I 0.66666666666667 0.33333333333333 0.58452629515785 126.904 8 I 0.00000000000000 0.00000000000000 0.24996113425933 126.904 9 I 0.00000000000000 0.00000000000000 0.75003886574067 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.890360047725579 0.000000000000000 0.000000000000000 b -8.445180023862788 14.627480880396096 0.000000000000000 c 0.000000000000000 0.000000000000000 20.593476540000001 Atomic positions (fractional): 1 Cd 0.16666666666667 0.08333333333333 0.33272442880812 112.411 > 1 2 Cd 0.41666666666667 0.08333333333333 0.33272442880812 112.411 > 1 3 Cd 0.66666666666667 0.08333333333333 0.33272442880812 112.411 > 1 4 Cd 0.91666666666667 0.08333333333333 0.33272442880812 112.411 > 1 5 Cd 0.16666666666667 0.33333333333333 0.33272442880812 112.411 > 1 6 Cd 0.41666666666667 0.33333333333333 0.33272442880812 112.411 > 1 7 Cd 0.66666666666667 0.33333333333333 0.33272442880812 112.411 > 1 8 Cd 0.91666666666667 0.33333333333333 0.33272442880812 112.411 > 1 9 Cd 0.16666666666667 0.58333333333333 0.33272442880812 112.411 > 1 10 Cd 0.41666666666667 0.58333333333333 0.33272442880812 112.411 > 1 11 Cd 0.66666666666667 0.58333333333333 0.33272442880812 112.411 > 1 12 Cd 0.91666666666667 0.58333333333333 0.33272442880812 112.411 > 1 13 Cd 0.16666666666667 0.83333333333333 0.33272442880812 112.411 > 1 14 Cd 0.41666666666667 0.83333333333333 0.33272442880812 112.411 > 1 15 Cd 0.66666666666667 0.83333333333333 0.33272442880812 112.411 > 1 16 Cd 0.91666666666667 0.83333333333333 0.33272442880812 112.411 > 1 17 Cd 0.08333333333333 0.16666666666667 0.66727557119188 112.411 > 17 18 Cd 0.33333333333333 0.16666666666667 0.66727557119188 112.411 > 17 19 Cd 0.58333333333333 0.16666666666667 0.66727557119188 112.411 > 17 20 Cd 0.83333333333333 0.16666666666667 0.66727557119188 112.411 > 17 21 Cd 0.08333333333333 0.41666666666667 0.66727557119188 112.411 > 17 22 Cd 0.33333333333333 0.41666666666667 0.66727557119188 112.411 > 17 23 Cd 0.58333333333333 0.41666666666667 0.66727557119188 112.411 > 17 24 Cd 0.83333333333333 0.41666666666667 0.66727557119188 112.411 > 17 25 Cd 0.08333333333333 0.66666666666667 0.66727557119188 112.411 > 17 26 Cd 0.33333333333333 0.66666666666667 0.66727557119188 112.411 > 17 27 Cd 0.58333333333333 0.66666666666667 0.66727557119188 112.411 > 17 28 Cd 0.83333333333333 0.66666666666667 0.66727557119188 112.411 > 17 29 Cd 0.08333333333333 0.91666666666667 0.66727557119188 112.411 > 17 30 Cd 0.33333333333333 0.91666666666667 0.66727557119188 112.411 > 17 31 Cd 0.58333333333333 0.91666666666667 0.66727557119188 112.411 > 17 32 Cd 0.83333333333333 0.91666666666667 0.66727557119188 112.411 > 17 33 Cd 0.00000000000000 0.00000000000000 0.00000000000000 112.411 > 33 34 Cd 0.25000000000000 0.00000000000000 0.00000000000000 112.411 > 33 35 Cd 0.50000000000000 0.00000000000000 0.00000000000000 112.411 > 33 36 Cd 0.75000000000000 0.00000000000000 0.00000000000000 112.411 > 33 37 Cd 0.00000000000000 0.25000000000000 0.00000000000000 112.411 > 33 38 Cd 0.25000000000000 0.25000000000000 0.00000000000000 112.411 > 33 39 Cd 0.50000000000000 0.25000000000000 0.00000000000000 112.411 > 33 40 Cd 0.75000000000000 0.25000000000000 0.00000000000000 112.411 > 33 41 Cd 0.00000000000000 0.50000000000000 0.00000000000000 112.411 > 33 42 Cd 0.25000000000000 0.50000000000000 0.00000000000000 112.411 > 33 43 Cd 0.50000000000000 0.50000000000000 0.00000000000000 112.411 > 33 44 Cd 0.75000000000000 0.50000000000000 0.00000000000000 112.411 > 33 45 Cd 0.00000000000000 0.75000000000000 0.00000000000000 112.411 > 33 46 Cd 0.25000000000000 0.75000000000000 0.00000000000000 112.411 > 33 47 Cd 0.50000000000000 0.75000000000000 0.00000000000000 112.411 > 33 48 Cd 0.75000000000000 0.75000000000000 0.00000000000000 112.411 > 33 49 I 0.08333333333333 0.16666666666667 0.08276777243055 126.904 > 49 50 I 0.33333333333333 0.16666666666667 0.08276777243055 126.904 > 49 51 I 0.58333333333333 0.16666666666667 0.08276777243055 126.904 > 49 52 I 0.83333333333333 0.16666666666667 0.08276777243055 126.904 > 49 53 I 0.08333333333333 0.41666666666667 0.08276777243055 126.904 > 49 54 I 0.33333333333333 0.41666666666667 0.08276777243055 126.904 > 49 55 I 0.58333333333333 0.41666666666667 0.08276777243055 126.904 > 49 56 I 0.83333333333333 0.41666666666667 0.08276777243055 126.904 > 49 57 I 0.08333333333333 0.66666666666667 0.08276777243055 126.904 > 49 58 I 0.33333333333333 0.66666666666667 0.08276777243055 126.904 > 49 59 I 0.58333333333333 0.66666666666667 0.08276777243055 126.904 > 49 60 I 0.83333333333333 0.66666666666667 0.08276777243055 126.904 > 49 61 I 0.08333333333333 0.91666666666667 0.08276777243055 126.904 > 49 62 I 0.33333333333333 0.91666666666667 0.08276777243055 126.904 > 49 63 I 0.58333333333333 0.91666666666667 0.08276777243055 126.904 > 49 64 I 0.83333333333333 0.91666666666667 0.08276777243055 126.904 > 49 65 I 0.08333333333333 0.16666666666667 0.41547370484215 126.904 > 65 66 I 0.33333333333333 0.16666666666667 0.41547370484215 126.904 > 65 67 I 0.58333333333333 0.16666666666667 0.41547370484215 126.904 > 65 68 I 0.83333333333333 0.16666666666667 0.41547370484215 126.904 > 65 69 I 0.08333333333333 0.41666666666667 0.41547370484215 126.904 > 65 70 I 0.33333333333333 0.41666666666667 0.41547370484215 126.904 > 65 71 I 0.58333333333333 0.41666666666667 0.41547370484215 126.904 > 65 72 I 0.83333333333333 0.41666666666667 0.41547370484215 126.904 > 65 73 I 0.08333333333333 0.66666666666667 0.41547370484215 126.904 > 65 74 I 0.33333333333333 0.66666666666667 0.41547370484215 126.904 > 65 75 I 0.58333333333333 0.66666666666667 0.41547370484215 126.904 > 65 76 I 0.83333333333333 0.66666666666667 0.41547370484215 126.904 > 65 77 I 0.08333333333333 0.91666666666667 0.41547370484215 126.904 > 65 78 I 0.33333333333333 0.91666666666667 0.41547370484215 126.904 > 65 79 I 0.58333333333333 0.91666666666667 0.41547370484215 126.904 > 65 80 I 0.83333333333333 0.91666666666667 0.41547370484215 126.904 > 65 81 I 0.16666666666667 0.08333333333333 0.91723222756945 126.904 > 81 82 I 0.41666666666667 0.08333333333333 0.91723222756945 126.904 > 81 83 I 0.66666666666667 0.08333333333333 0.91723222756945 126.904 > 81 84 I 0.91666666666667 0.08333333333333 0.91723222756945 126.904 > 81 85 I 0.16666666666667 0.33333333333333 0.91723222756945 126.904 > 81 86 I 0.41666666666667 0.33333333333333 0.91723222756945 126.904 > 81 87 I 0.66666666666667 0.33333333333333 0.91723222756945 126.904 > 81 88 I 0.91666666666667 0.33333333333333 0.91723222756945 126.904 > 81 89 I 0.16666666666667 0.58333333333333 0.91723222756945 126.904 > 81 90 I 0.41666666666667 0.58333333333333 0.91723222756945 126.904 > 81 91 I 0.66666666666667 0.58333333333333 0.91723222756945 126.904 > 81 92 I 0.91666666666667 0.58333333333333 0.91723222756945 126.904 > 81 93 I 0.16666666666667 0.83333333333333 0.91723222756945 126.904 > 81 94 I 0.41666666666667 0.83333333333333 0.91723222756945 126.904 > 81 95 I 0.66666666666667 0.83333333333333 0.91723222756945 126.904 > 81 96 I 0.91666666666667 0.83333333333333 0.91723222756945 126.904 > 81 97 I 0.16666666666667 0.08333333333333 0.58452629515785 126.904 > 97 98 I 0.41666666666667 0.08333333333333 0.58452629515785 126.904 > 97 99 I 0.66666666666667 0.08333333333333 0.58452629515785 126.904 > 97 100 I 0.91666666666667 0.08333333333333 0.58452629515785 126.904 > 97 101 I 0.16666666666667 0.33333333333333 0.58452629515785 126.904 > 97 102 I 0.41666666666667 0.33333333333333 0.58452629515785 126.904 > 97 103 I 0.66666666666667 0.33333333333333 0.58452629515785 126.904 > 97 104 I 0.91666666666667 0.33333333333333 0.58452629515785 126.904 > 97 105 I 0.16666666666667 0.58333333333333 0.58452629515785 126.904 > 97 106 I 0.41666666666667 0.58333333333333 0.58452629515785 126.904 > 97 107 I 0.66666666666667 0.58333333333333 0.58452629515785 126.904 > 97 108 I 0.91666666666667 0.58333333333333 0.58452629515785 126.904 > 97 109 I 0.16666666666667 0.83333333333333 0.58452629515785 126.904 > 97 110 I 0.41666666666667 0.83333333333333 0.58452629515785 126.904 > 97 111 I 0.66666666666667 0.83333333333333 0.58452629515785 126.904 > 97 112 I 0.91666666666667 0.83333333333333 0.58452629515785 126.904 > 97 113 I 0.00000000000000 0.00000000000000 0.24996113425933 126.904 > 113 114 I 0.25000000000000 0.00000000000000 0.24996113425933 126.904 > 113 115 I 0.50000000000000 0.00000000000000 0.24996113425933 126.904 > 113 116 I 0.75000000000000 0.00000000000000 0.24996113425933 126.904 > 113 117 I 0.00000000000000 0.25000000000000 0.24996113425933 126.904 > 113 118 I 0.25000000000000 0.25000000000000 0.24996113425933 126.904 > 113 119 I 0.50000000000000 0.25000000000000 0.24996113425933 126.904 > 113 120 I 0.75000000000000 0.25000000000000 0.24996113425933 126.904 > 113 121 I 0.00000000000000 0.50000000000000 0.24996113425933 126.904 > 113 122 I 0.25000000000000 0.50000000000000 0.24996113425933 126.904 > 113 123 I 0.50000000000000 0.50000000000000 0.24996113425933 126.904 > 113 124 I 0.75000000000000 0.50000000000000 0.24996113425933 126.904 > 113 125 I 0.00000000000000 0.75000000000000 0.24996113425933 126.904 > 113 126 I 0.25000000000000 0.75000000000000 0.24996113425933 126.904 > 113 127 I 0.50000000000000 0.75000000000000 0.24996113425933 126.904 > 113 128 I 0.75000000000000 0.75000000000000 0.24996113425933 126.904 > 113 129 I 0.00000000000000 0.00000000000000 0.75003886574067 126.904 > 129 130 I 0.25000000000000 0.00000000000000 0.75003886574067 126.904 > 129 131 I 0.50000000000000 0.00000000000000 0.75003886574067 126.904 > 129 132 I 0.75000000000000 0.00000000000000 0.75003886574067 126.904 > 129 133 I 0.00000000000000 0.25000000000000 0.75003886574067 126.904 > 129 134 I 0.25000000000000 0.25000000000000 0.75003886574067 126.904 > 129 135 I 0.50000000000000 0.25000000000000 0.75003886574067 126.904 > 129 136 I 0.75000000000000 0.25000000000000 0.75003886574067 126.904 > 129 137 I 0.00000000000000 0.50000000000000 0.75003886574067 126.904 > 129 138 I 0.25000000000000 0.50000000000000 0.75003886574067 126.904 > 129 139 I 0.50000000000000 0.50000000000000 0.75003886574067 126.904 > 129 140 I 0.75000000000000 0.50000000000000 0.75003886574067 126.904 > 129 141 I 0.00000000000000 0.75000000000000 0.75003886574067 126.904 > 129 142 I 0.25000000000000 0.75000000000000 0.75003886574067 126.904 > 129 143 I 0.50000000000000 0.75000000000000 0.75003886574067 126.904 > 129 144 I 0.75000000000000 0.75000000000000 0.75003886574067 126.904 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 6.0028127 -0.0000000 0.0000000 0.0000000 6.0028127 0.0000000 0.0000000 0.0000000 4.8454435 -------------------------- Born effective charges -------------------------- 1 Cd 3.0166357 -0.0000000 0.0000000 0.0000000 3.0166357 0.0000000 0.0000000 0.0000000 1.6823478 2 Cd 3.0166357 -0.0000000 0.0000000 0.0000000 3.0166357 0.0000000 0.0000000 0.0000000 1.6823478 3 Cd 3.0061572 -0.0000000 0.0000000 0.0000000 3.0061572 0.0000000 0.0000000 0.0000000 1.6589443 4 I -1.5130415 0.0000000 0.0000000 0.0000000 -1.5130415 0.0000000 0.0000000 0.0000000 -0.8380966 5 I -1.4892058 0.0000000 0.0000000 0.0000000 -1.4892058 0.0000000 0.0000000 0.0000000 -0.8488195 6 I -1.5130415 0.0000000 0.0000000 0.0000000 -1.5130415 0.0000000 0.0000000 0.0000000 -0.8380966 7 I -1.4892058 0.0000000 0.0000000 0.0000000 -1.4892058 0.0000000 0.0000000 0.0000000 -0.8488195 8 I -1.5174671 0.0000000 0.0000000 0.0000000 -1.5174671 0.0000000 0.0000000 0.0000000 -0.8249038 9 I -1.5174671 0.0000000 0.0000000 0.0000000 -1.5174671 0.0000000 0.0000000 0.0000000 -0.8249038 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000003 (yzy) 0.00000003 (yzy) 0.00000003 (yyz) Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 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.336 ( 0.000 0.000 0.000) 0.000 0.336 ( 0.000 0.000 0.000) 0.000 0.351 ( 0.000 0.000 0.000) 0.000 0.351 ( 0.000 0.000 0.000) 0.000 0.640 ( 0.000 0.000 0.000) 0.000 0.650 ( 0.000 0.000 0.000) 0.000 1.533 ( 0.000 0.000 0.000) 0.000 1.533 ( 0.000 0.000 0.000) 0.000 1.574 ( 0.000 0.000 0.000) 0.000 1.574 ( 0.000 0.000 0.000) 0.000 1.623 ( 0.000 0.000 0.000) 0.000 1.623 ( 0.000 0.000 0.000) 0.000 2.583 ( 0.000 0.000 0.000) 0.000 2.583 ( 0.000 0.000 0.000) 0.000 2.603 ( 0.000 0.000 0.000) 0.000 2.603 ( 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 3.147 ( 0.000 0.000 0.000) 0.000 3.151 ( 0.000 0.000 0.000) 0.000 3.213 ( 0.000 0.000 0.000) 0.000 4.029 ( 0.000 0.000 0.000) 0.000 4.405 ( 0.000 0.000 0.000) 0.000 4.407 ( 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.642 4.412 -0.224) 8.827 0.284 ( 12.284 7.092 0.161) 14.185 0.442 ( 7.995 4.616 0.380) 9.239 0.451 ( 7.699 4.445 -0.486) 8.904 0.527 ( 22.298 12.874 -0.063) 25.748 0.559 ( 6.068 3.503 -0.359) 7.016 0.563 ( 7.771 4.486 0.672) 8.998 0.670 ( 9.722 5.613 2.861) 11.585 0.674 ( 10.245 5.915 -3.012) 12.208 1.545 ( 0.952 0.550 0.065) 1.101 1.583 ( 0.750 0.433 -0.051) 0.867 1.634 ( 0.860 0.497 -0.010) 0.993 1.697 ( 12.679 7.320 0.184) 14.642 1.718 ( 11.179 6.454 -0.660) 12.925 1.746 ( 10.540 6.086 0.597) 12.186 2.594 ( 1.041 0.601 0.141) 1.210 2.620 ( 1.392 0.804 -0.677) 1.744 2.625 ( 1.188 0.686 0.537) 1.473 2.820 ( 14.188 8.192 -4.347) 16.950 2.820 ( 13.929 8.042 4.368) 16.666 3.157 ( 1.953 1.128 -0.071) 2.257 3.161 ( 1.947 1.124 0.057) 2.249 3.188 ( -2.021 -1.167 0.007) 2.333 3.855 ( 0.535 0.309 4.187) 4.232 3.957 ( -6.121 -3.534 -4.319) 8.283 4.290 ( -9.612 -5.549 0.822) 11.129 4.294 ( -9.176 -5.298 -0.661) 10.616 ======================= 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.337 ( 9.514 5.493 -0.292) 10.990 0.548 ( 11.019 6.362 -0.142) 12.724 0.586 ( 1.139 0.657 4.087) 4.293 0.600 ( 1.684 0.973 -3.694) 4.175 0.642 ( 9.102 5.255 2.139) 10.726 0.652 ( 9.420 5.439 -2.001) 11.060 1.013 ( 19.770 11.414 -0.242) 22.829 1.047 ( 17.287 9.981 -0.273) 19.963 1.051 ( 17.885 10.326 0.385) 20.655 1.568 ( 0.863 0.498 0.434) 1.087 1.603 ( 0.880 0.508 -0.653) 1.208 1.654 ( 0.769 0.444 0.227) 0.916 2.010 ( 14.143 8.165 -0.201) 16.332 2.029 ( 15.416 8.901 0.149) 17.802 2.052 ( 15.078 8.705 0.045) 17.411 2.633 ( 2.401 1.386 -0.013) 2.772 2.664 ( 2.351 1.357 0.400) 2.743 2.665 ( 2.444 1.411 -0.386) 2.848 2.992 ( 1.297 0.749 2.471) 2.890 2.999 ( 1.487 0.859 -2.497) 3.031 3.128 ( -3.061 -1.767 0.047) 3.535 3.315 ( 9.682 5.590 3.131) 11.610 3.325 ( 10.757 6.211 -3.795) 12.988 3.757 ( -9.729 -5.617 7.287) 13.391 3.845 ( -4.866 -2.809 -7.627) 9.474 4.020 ( -12.711 -7.339 2.278) 14.854 4.066 ( -7.875 -4.547 -1.221) 9.175 ======================= 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.553 ( 9.372 5.411 -0.208) 10.824 0.641 ( 3.777 2.181 1.714) 4.686 0.669 ( 4.367 2.522 -1.388) 5.231 0.776 ( 9.122 5.267 -0.370) 10.540 0.838 ( 8.143 4.701 1.806) 9.575 0.855 ( 8.364 4.829 -1.447) 9.766 1.392 ( 13.778 7.955 -1.252) 15.959 1.401 ( 13.872 8.009 1.156) 16.059 1.437 ( 15.778 9.110 -0.021) 18.219 1.578 ( -0.052 -0.030 0.750) 0.752 1.617 ( 0.312 0.180 -1.334) 1.382 1.663 ( -0.020 -0.012 0.595) 0.595 2.316 ( 12.080 6.975 0.527) 13.959 2.354 ( 12.666 7.313 -1.533) 14.706 2.381 ( 13.337 7.700 1.014) 15.433 2.699 ( 3.290 1.900 -0.201) 3.805 2.723 ( 2.938 1.696 0.670) 3.459 2.733 ( 3.352 1.935 -0.471) 3.899 2.990 ( -0.381 -0.220 0.981) 1.075 3.004 ( 0.004 0.002 -1.090) 1.090 3.059 ( -2.919 -1.685 0.152) 3.374 3.379 ( -6.122 -3.534 1.038) 7.145 3.425 ( -5.117 -2.954 -2.988) 6.621 3.522 ( -12.865 -7.428 4.178) 15.431 3.720 ( 1.757 1.015 -6.946) 7.236 3.832 ( -1.367 -0.789 6.298) 6.493 4.012 ( 0.495 0.286 -1.578) 1.678 ======================= 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.734 ( 6.449 3.723 -0.128) 7.447 0.768 ( 7.262 4.193 0.902) 8.434 0.801 ( 7.611 4.394 -0.701) 8.817 0.957 ( 6.971 4.025 -0.386) 8.059 1.004 ( 6.443 3.720 1.409) 7.572 1.023 ( 6.456 3.727 -1.032) 7.525 1.566 ( -0.948 -0.547 0.797) 1.354 1.616 ( -0.349 -0.201 -1.742) 1.788 1.649 ( 9.139 5.276 -2.245) 10.789 1.655 ( -0.670 -0.387 0.956) 1.230 1.663 ( 9.374 5.412 2.399) 11.087 1.738 ( 10.721 6.190 -0.243) 12.382 2.505 ( 3.308 1.910 0.838) 3.910 2.562 ( 4.438 2.562 -2.204) 5.578 2.610 ( 5.430 3.135 1.432) 6.432 2.774 ( 3.263 1.884 -0.331) 3.782 2.792 ( 2.980 1.720 0.623) 3.496 2.809 ( 3.302 1.907 -0.294) 3.825 2.992 ( 0.732 0.423 -0.530) 0.998 3.005 ( -1.449 -0.836 0.446) 1.731 3.017 ( 1.330 0.768 0.116) 1.540 3.214 ( -5.104 -2.947 -0.210) 5.898 3.239 ( -8.717 -5.033 1.094) 10.125 3.245 ( -7.322 -4.227 -0.810) 8.493 3.821 ( 4.442 2.564 -3.973) 6.488 3.882 ( 3.220 1.859 5.077) 6.293 4.028 ( 0.632 0.365 -1.234) 1.434 ======================= 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.887 ( 7.200 4.157 -1.166) 8.395 0.933 ( 7.352 4.245 2.251) 8.783 0.986 ( 7.953 4.592 -1.063) 9.245 1.091 ( 4.921 2.841 -0.335) 5.693 1.125 ( 4.289 2.476 1.092) 5.071 1.145 ( 4.322 2.496 -0.780) 5.052 1.538 ( -1.413 -0.816 0.748) 1.795 1.606 ( -0.423 -0.244 -2.053) 2.110 1.637 ( -0.847 -0.489 1.330) 1.651 1.813 ( 5.665 3.271 -1.972) 6.833 1.831 ( 5.729 3.308 2.236) 6.983 1.919 ( 5.499 3.175 -0.296) 6.357 2.449 ( -6.782 -3.916 0.769) 7.869 2.518 ( -7.009 -4.047 -2.608) 8.503 2.567 ( -7.748 -4.473 1.894) 9.144 2.840 ( 2.447 1.413 -0.311) 2.843 2.852 ( 2.325 1.343 0.483) 2.728 2.876 ( 2.474 1.428 -0.171) 2.862 3.030 ( 3.127 1.805 0.462) 3.640 3.039 ( 3.280 1.894 -0.793) 3.870 3.069 ( 3.001 1.733 0.488) 3.499 3.153 ( 0.931 0.538 -0.230) 1.099 3.188 ( 1.570 0.907 -0.347) 1.846 3.194 ( 1.281 0.740 0.396) 1.532 3.897 ( 2.401 1.386 -2.506) 3.737 3.938 ( 1.626 0.939 3.388) 3.873 4.035 ( -0.055 -0.032 -0.904) 0.907 ======================= 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.043 ( 6.079 3.510 -1.015) 7.093 1.086 ( 5.729 3.308 1.851) 6.870 1.133 ( 4.540 2.621 -0.838) 5.309 1.179 ( 2.791 1.612 -0.346) 3.242 1.196 ( 1.997 1.153 0.790) 2.438 1.218 ( 2.187 1.262 -0.496) 2.573 1.507 ( -1.181 -0.682 0.515) 1.458 1.601 ( -0.073 -0.042 -2.075) 2.077 1.619 ( -0.644 -0.372 1.612) 1.775 1.916 ( 3.734 2.156 -1.136) 4.458 1.933 ( 3.532 2.039 1.419) 4.318 2.000 ( 1.946 1.123 -0.236) 2.259 2.278 ( -7.192 -4.152 0.294) 8.310 2.343 ( -6.861 -3.961 -2.024) 8.177 2.374 ( -7.672 -4.429 1.688) 9.018 2.882 ( 1.251 0.722 -0.148) 1.452 2.893 ( 1.253 0.723 0.224) 1.464 2.918 ( 1.277 0.737 -0.075) 1.477 3.096 ( 2.382 1.375 0.349) 2.773 3.109 ( 2.380 1.374 -0.570) 2.806 3.130 ( 2.021 1.167 0.248) 2.347 3.209 ( 2.664 1.538 -0.033) 3.076 3.238 ( 2.103 1.214 -0.042) 2.428 3.242 ( 2.053 1.185 0.050) 2.371 3.936 ( 1.160 0.670 -1.784) 2.231 3.958 ( 0.310 0.179 2.265) 2.293 4.028 ( -0.411 -0.237 -0.493) 0.684 ======================= 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.122 ( 0.000 0.000 0.000) 0.000 1.169 ( -0.000 -0.000 -0.000) 0.000 1.172 ( 0.000 0.000 0.000) 0.000 1.213 ( 0.000 0.000 -0.000) 0.000 1.217 ( -0.000 -0.000 0.000) 0.000 1.243 ( 0.000 0.000 -0.000) 0.000 1.492 ( 0.000 0.000 0.000) 0.000 1.601 ( 0.000 0.000 0.000) 0.000 1.611 ( 0.000 0.000 -0.000) 0.000 1.974 ( -0.000 -0.000 0.000) 0.000 1.976 ( 0.000 0.000 -0.000) 0.000 2.020 ( 0.000 0.000 0.000) 0.000 2.178 ( 0.000 0.000 0.000) 0.000 2.257 ( 0.000 0.000 -0.000) 0.000 2.277 ( 0.000 0.000 0.000) 0.000 2.896 ( 0.000 0.000 -0.000) 0.000 2.908 ( 0.000 0.000 0.000) 0.000 2.933 ( 0.000 0.000 -0.000) 0.000 3.127 ( 0.000 0.000 0.000) 0.000 3.135 ( 0.000 0.000 -0.000) 0.000 3.153 ( 0.000 0.000 0.000) 0.000 3.242 ( 0.000 0.000 -0.000) 0.000 3.263 ( 0.000 0.000 -0.000) 0.000 3.266 ( 0.000 0.000 -0.000) 0.000 3.951 ( 0.000 0.000 0.000) 0.000 3.958 ( 0.000 0.000 -0.000) 0.000 4.023 ( 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.276 ( 5.211 9.026 0.000) 10.422 0.505 ( 7.555 13.086 0.000) 15.111 0.576 ( 0.824 1.427 0.000) 1.648 0.578 ( 1.006 1.743 0.000) 2.013 0.607 ( 5.629 9.750 0.000) 11.259 0.615 ( 5.281 9.147 0.000) 10.562 0.883 ( 11.702 20.268 0.000) 23.404 0.928 ( 9.887 17.125 0.000) 19.775 0.931 ( 9.769 16.921 0.000) 19.539 1.612 ( 2.365 4.096 0.000) 4.730 1.634 ( 2.091 3.622 0.000) 4.182 1.685 ( 2.095 3.629 0.000) 4.190 1.885 ( 6.856 11.875 0.000) 13.712 1.902 ( 6.120 10.600 0.000) 12.239 1.904 ( 6.530 11.310 0.000) 13.059 2.630 ( 1.637 2.835 0.000) 3.274 2.657 ( 1.873 3.245 0.000) 3.746 2.660 ( 1.766 3.060 0.000) 3.533 2.977 ( 2.285 3.957 0.000) 4.569 2.982 ( 2.334 4.042 0.000) 4.668 3.147 ( -1.673 -2.898 0.000) 3.346 3.257 ( 5.918 10.251 0.000) 11.837 3.261 ( 6.007 10.404 0.000) 12.014 3.823 ( -5.519 -9.559 0.000) 11.037 3.874 ( 0.570 0.986 0.000) 1.139 4.100 ( -7.557 -13.089 0.000) 15.114 4.101 ( -7.608 -13.177 0.000) 15.216 ======================= 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.470 ( 7.472 7.862 -0.142) 10.847 0.608 ( 2.150 1.915 2.438) 3.772 0.620 ( 3.142 1.621 -2.228) 4.179 0.753 ( 6.112 13.233 -0.178) 14.577 0.817 ( 5.413 11.666 1.571) 12.956 0.825 ( 6.008 11.345 -1.395) 12.913 1.256 ( 12.056 12.008 -0.541) 17.024 1.261 ( 12.150 11.217 0.321) 16.540 1.277 ( 14.109 12.396 0.154) 18.782 1.668 ( -3.145 8.522 0.599) 9.103 1.691 ( -2.151 7.601 -0.879) 7.948 1.738 ( -2.520 7.397 0.274) 7.819 2.118 ( 12.311 2.389 0.127) 12.541 2.142 ( 13.185 2.412 -0.436) 13.411 2.155 ( 13.957 1.703 0.320) 14.064 2.709 ( 2.575 6.392 -0.032) 6.891 2.742 ( 2.217 6.662 0.692) 7.056 2.744 ( 2.474 6.604 -0.673) 7.084 2.995 ( -0.390 -0.235 1.389) 1.462 3.002 ( -0.033 -0.406 -1.411) 1.469 3.084 ( -2.319 -2.634 0.055) 3.510 3.441 ( 0.583 5.954 2.230) 6.385 3.463 ( 2.971 6.375 -3.193) 7.724 3.617 ( -5.613 -8.485 1.914) 10.352 3.773 ( -12.012 -5.162 -2.913) 13.395 3.838 ( -5.874 -10.157 3.095) 12.134 3.945 ( 3.904 -6.108 -1.129) 7.337 ======================= 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.666 ( 6.744 5.818 0.322) 8.913 0.693 ( 5.395 3.437 0.210) 6.400 0.719 ( 6.023 3.105 -0.464) 6.792 0.962 ( 3.475 12.900 -0.335) 13.364 1.007 ( 3.276 11.760 1.370) 12.284 1.020 ( 3.657 11.417 -1.028) 12.032 1.533 ( 10.197 6.138 -1.910) 12.054 1.543 ( 10.637 6.133 1.976) 12.437 1.602 ( 12.185 7.365 -0.140) 14.239 1.677 ( -5.715 9.303 0.671) 10.939 1.709 ( -4.725 8.537 -1.357) 9.851 1.746 ( -5.068 7.827 0.680) 9.349 2.321 ( 12.116 -5.054 0.424) 13.135 2.353 ( 12.725 -5.514 -1.084) 13.911 2.373 ( 13.650 -6.334 0.682) 15.063 2.837 ( 2.581 10.021 -0.037) 10.348 2.855 ( 1.250 9.660 0.528) 9.755 2.861 ( 1.661 9.517 -0.493) 9.673 2.989 ( 0.002 0.205 0.198) 0.285 3.004 ( 0.783 0.225 -0.368) 0.894 3.024 ( -2.024 -1.716 0.196) 2.660 3.309 ( -9.775 -3.548 0.401) 10.407 3.361 ( -10.920 -3.221 -2.350) 11.625 3.395 ( -13.884 -5.317 2.449) 15.068 3.738 ( 6.236 0.292 -5.043) 8.025 3.804 ( 5.559 -2.075 5.676) 8.211 3.969 ( 3.859 -4.066 -1.141) 5.721 ======================= 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.809 ( 6.326 4.104 -0.838) 7.587 0.844 ( 7.188 3.741 1.793) 8.299 0.889 ( 8.323 4.644 -0.925) 9.576 1.121 ( 1.006 11.985 -0.334) 12.032 1.154 ( 0.918 10.961 1.080) 11.052 1.168 ( 1.114 10.584 -0.739) 10.668 1.645 ( -5.792 7.289 0.451) 9.321 1.680 ( -2.719 5.097 -1.346) 5.932 1.705 ( -1.588 4.068 0.812) 4.442 1.753 ( 3.384 6.307 -1.319) 7.278 1.763 ( 4.714 5.748 1.535) 7.590 1.834 ( 8.139 3.669 -0.180) 8.929 2.396 ( 4.369 -9.051 0.611) 10.069 2.436 ( 5.487 -10.211 -1.440) 11.681 2.467 ( 6.566 -10.764 0.824) 12.635 2.932 ( -2.792 7.226 0.158) 7.748 2.945 ( -3.712 7.061 -0.416) 7.989 2.951 ( -2.607 7.034 0.313) 7.509 3.009 ( 2.540 3.584 0.467) 4.418 3.018 ( 2.326 3.776 -0.691) 4.489 3.046 ( 2.486 3.334 0.254) 4.166 3.172 ( -5.249 -1.657 0.109) 5.505 3.181 ( -2.692 -1.256 -0.136) 2.973 3.198 ( -3.652 -1.224 0.013) 3.851 3.832 ( 5.055 -1.276 -2.913) 5.972 3.877 ( 4.721 -2.241 3.818) 6.472 3.987 ( 3.022 -4.280 -0.950) 5.325 ======================= 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.964 ( 6.941 3.764 -1.145) 7.978 1.003 ( 7.136 3.182 1.972) 8.058 1.060 ( 6.825 3.062 -0.825) 7.526 1.231 ( -1.110 10.614 -0.304) 10.676 1.250 ( -1.452 9.535 0.771) 9.675 1.265 ( -1.107 9.177 -0.478) 9.256 1.611 ( -6.120 7.227 0.416) 9.479 1.667 ( -3.820 5.166 -1.706) 6.647 1.685 ( -3.976 4.315 1.301) 6.010 1.875 ( 4.196 3.655 -1.034) 5.660 1.892 ( 4.367 3.557 1.305) 5.782 1.960 ( 4.180 1.251 -0.268) 4.372 2.311 ( -3.234 -7.773 0.665) 8.445 2.350 ( -2.315 -9.927 -1.639) 10.324 2.388 ( -2.416 -10.138 0.952) 10.465 2.931 ( -2.475 4.266 -0.028) 4.932 2.940 ( -1.930 4.287 0.248) 4.708 2.949 ( -2.185 3.997 -0.188) 4.559 3.096 ( 0.555 5.350 0.155) 5.381 3.100 ( 1.206 4.640 -0.234) 4.800 3.121 ( 1.346 4.148 0.220) 4.367 3.160 ( 3.611 -0.102 -0.229) 3.620 3.195 ( 3.011 -0.369 0.090) 3.035 3.200 ( 2.855 -0.332 -0.005) 2.875 3.886 ( 3.592 -2.390 -1.892) 4.711 3.916 ( 3.046 -3.019 2.533) 4.981 3.989 ( 2.288 -4.486 -0.657) 5.078 ======================= 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.090 ( 3.912 1.341 -0.381) 4.153 1.127 ( 4.323 1.282 0.813) 4.582 1.156 ( 2.142 0.002 -0.437) 2.186 1.288 ( -3.287 8.653 -0.339) 9.262 1.292 ( -3.544 8.045 0.485) 8.804 1.311 ( -3.062 7.541 -0.160) 8.141 1.587 ( -5.210 7.710 0.180) 9.307 1.659 ( -2.821 4.907 -1.820) 5.946 1.666 ( -3.200 4.308 1.653) 5.616 1.966 ( 2.256 3.361 -0.608) 4.093 1.979 ( 1.605 3.139 0.910) 3.641 2.013 ( 1.288 0.773 -0.219) 1.518 2.181 ( -3.393 -4.921 0.033) 5.977 2.229 ( -1.013 -6.379 -0.903) 6.522 2.253 ( -1.763 -6.642 0.788) 6.917 2.939 ( -1.654 3.742 -0.109) 4.092 2.948 ( -1.720 3.456 0.247) 3.868 2.959 ( -0.953 2.511 -0.133) 2.689 3.139 ( -0.235 3.105 0.159) 3.118 3.146 ( -0.173 2.570 -0.183) 2.583 3.166 ( -0.299 2.933 -0.007) 2.949 3.213 ( 2.490 -1.094 0.050) 2.720 3.238 ( 2.297 -1.247 0.176) 2.620 3.240 ( 2.129 -1.327 -0.191) 2.516 3.912 ( 2.683 -2.958 -1.245) 4.183 3.925 ( 1.953 -3.384 1.476) 4.177 3.982 ( 2.184 -4.369 -0.238) 4.890 ======================= 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.634 ( 4.988 8.639 0.000) 9.975 0.666 ( 2.317 4.014 0.000) 4.635 0.677 ( 2.581 4.471 0.000) 5.162 1.017 ( 7.500 12.991 0.000) 15.000 1.057 ( 6.810 11.795 0.000) 13.619 1.058 ( 6.950 12.038 0.000) 13.900 1.464 ( 5.531 9.580 0.000) 11.062 1.467 ( 5.566 9.641 0.000) 11.133 1.518 ( 6.729 11.656 0.000) 13.459 1.840 ( 4.123 7.142 0.000) 8.246 1.845 ( 3.849 6.667 0.000) 7.699 1.892 ( 3.797 6.577 0.000) 7.594 2.131 ( 0.303 0.524 0.000) 0.606 2.157 ( 0.225 0.390 0.000) 0.450 2.163 ( 0.643 1.114 0.000) 1.287 2.860 ( 4.749 8.225 0.000) 9.498 2.891 ( 4.332 7.503 0.000) 8.664 2.894 ( 4.612 7.988 0.000) 9.224 2.995 ( 0.088 0.153 0.000) 0.177 2.995 ( 0.164 0.284 0.000) 0.328 3.033 ( -1.348 -2.335 0.000) 2.697 3.416 ( -5.298 -9.177 0.000) 10.596 3.483 ( -3.722 -6.446 0.000) 7.443 3.489 ( -3.628 -6.284 0.000) 7.257 3.693 ( -0.129 -0.224 0.000) 0.258 3.701 ( -0.122 -0.211 0.000) 0.243 3.881 ( -0.562 -0.974 0.000) 1.124 ======================= 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.118 4.958 -0.399) 6.457 0.781 ( 5.009 5.009 0.853) 7.135 0.814 ( 6.372 6.483 -0.436) 9.101 1.239 ( 2.595 14.509 -0.190) 14.740 1.260 ( 2.477 13.319 0.587) 13.560 1.270 ( 2.785 13.352 -0.386) 13.645 1.629 ( 6.060 3.773 -2.056) 7.429 1.636 ( 6.666 3.531 2.151) 7.844 1.718 ( 7.210 4.458 -0.135) 8.478 1.896 ( -7.370 12.058 0.258) 14.134 1.909 ( -6.116 10.600 -0.681) 12.257 1.933 ( -7.755 10.196 0.443) 12.818 2.190 ( 9.278 -7.039 0.315) 11.651 2.211 ( 9.378 -7.980 -0.668) 12.331 2.223 ( 9.633 -7.685 0.349) 12.328 2.982 ( 0.571 1.590 -0.192) 1.700 2.990 ( 1.223 2.018 0.135) 2.363 2.997 ( -0.981 -0.926 0.179) 1.361 3.016 ( 2.937 5.244 0.018) 6.010 3.035 ( 2.887 6.187 -0.202) 6.830 3.045 ( 2.637 4.614 0.075) 5.315 3.253 ( -6.410 -1.446 0.007) 6.571 3.310 ( -8.242 -1.485 -1.374) 8.487 3.321 ( -8.685 -1.836 1.446) 8.994 3.731 ( 3.807 -0.790 -3.281) 5.088 3.757 ( 5.140 -2.502 3.702) 6.811 3.876 ( 2.925 -4.933 -0.509) 5.758 ======================= 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.897 ( 5.660 4.909 -0.891) 7.546 0.925 ( 6.229 4.664 1.445) 7.914 0.981 ( 6.500 4.755 -0.544) 8.072 1.386 ( -1.633 14.078 -0.243) 14.174 1.395 ( -1.622 12.716 0.380) 12.825 1.406 ( -1.776 12.853 -0.128) 12.976 1.727 ( 2.063 1.097 -1.659) 2.865 1.741 ( 2.153 0.824 1.787) 2.917 1.804 ( -2.973 3.281 -0.178) 4.431 1.907 ( -2.221 7.830 0.196) 8.141 1.908 ( -0.374 6.350 -0.034) 6.361 1.920 ( -2.012 8.234 -0.134) 8.478 2.242 ( 3.244 -5.443 0.524) 6.358 2.257 ( 3.386 -6.201 -0.510) 7.084 2.284 ( 4.424 -6.637 -0.021) 7.977 2.977 ( -0.894 0.139 0.235) 0.935 2.986 ( -0.865 0.623 -0.466) 1.163 3.002 ( -0.471 0.486 0.323) 0.750 3.095 ( 2.725 -0.937 -0.238) 2.892 3.114 ( 2.959 0.320 0.290) 2.990 3.117 ( 3.105 0.008 -0.139) 3.108 3.199 ( -2.878 5.897 -0.358) 6.572 3.224 ( -5.181 7.588 0.948) 9.237 3.232 ( -4.423 6.594 -0.561) 7.960 3.781 ( 4.525 -3.635 -1.855) 6.093 3.807 ( 4.780 -4.489 2.399) 6.983 3.877 ( 3.715 -6.581 -0.569) 7.578 ======================= 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.035 ( 5.026 3.445 -0.535) 6.117 1.066 ( 5.204 3.258 0.944) 6.212 1.111 ( 4.159 2.209 -0.409) 4.727 1.457 ( -4.116 10.355 0.118) 11.144 1.462 ( -4.245 11.581 -0.229) 12.336 1.466 ( -4.257 10.329 0.113) 11.172 1.736 ( -2.052 2.596 -0.615) 3.366 1.743 ( -1.718 1.644 0.329) 2.400 1.762 ( -4.301 3.618 0.231) 5.626 1.979 ( 2.228 6.130 0.984) 6.596 1.991 ( 2.417 4.917 -0.641) 5.516 2.002 ( 1.842 6.091 -0.327) 6.372 2.195 ( -2.319 -4.686 -0.253) 5.234 2.209 ( -3.725 -1.482 0.853) 4.099 2.237 ( -3.039 -3.931 -0.575) 5.002 2.975 ( 0.106 0.798 0.319) 0.866 2.983 ( -0.125 0.637 -0.409) 0.767 2.996 ( -0.677 0.601 0.101) 0.911 3.129 ( 3.634 -2.451 -0.105) 4.384 3.153 ( 3.016 -2.319 0.106) 3.806 3.162 ( 3.752 -2.813 0.031) 4.689 3.224 ( -2.262 6.469 -0.389) 6.864 3.238 ( -3.024 6.982 0.642) 7.635 3.249 ( -2.790 7.150 -0.273) 7.680 3.813 ( 4.159 -4.760 -1.025) 6.404 3.832 ( 3.803 -5.196 1.384) 6.586 3.873 ( 3.693 -6.971 -0.371) 7.897 ======================= 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.102 ( 0.008 -0.014 0.000) 0.016 1.143 ( 0.116 -0.201 0.000) 0.232 1.146 ( 0.167 -0.289 -0.000) 0.334 1.471 ( -5.370 9.302 -0.000) 10.741 1.475 ( -4.968 8.604 0.000) 9.935 1.480 ( -5.401 9.354 0.000) 10.802 1.728 ( -1.779 3.082 0.000) 3.559 1.739 ( -2.455 4.253 -0.000) 4.911 1.742 ( -1.298 2.248 0.000) 2.596 2.059 ( -2.637 4.568 0.000) 5.275 2.060 ( -3.127 5.415 -0.000) 6.253 2.072 ( -2.618 4.534 -0.000) 5.235 2.114 ( 0.258 -0.447 0.000) 0.516 2.146 ( 0.104 -0.180 -0.000) 0.208 2.162 ( 0.555 -0.962 -0.000) 1.111 2.986 ( -0.299 0.517 0.000) 0.597 2.986 ( -0.225 0.389 -0.000) 0.449 2.990 ( -0.348 0.602 -0.000) 0.695 3.154 ( 2.240 -3.880 0.000) 4.481 3.166 ( 1.722 -2.982 0.000) 3.443 3.181 ( 2.501 -4.331 0.000) 5.001 3.239 ( -3.476 6.020 -0.000) 6.951 3.244 ( -2.770 4.798 0.000) 5.540 3.261 ( -3.344 5.792 -0.000) 6.688 3.827 ( 3.188 -5.521 -0.000) 6.375 3.837 ( 3.103 -5.374 0.000) 6.206 3.869 ( 3.924 -6.797 0.000) 7.849 ======================= 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.873 ( 3.567 6.178 0.000) 7.133 0.893 ( 3.748 6.492 0.000) 7.497 0.947 ( 3.944 6.831 0.000) 7.888 1.491 ( 5.783 10.017 0.000) 11.566 1.499 ( 6.041 10.464 0.000) 12.082 1.509 ( 5.818 10.078 0.000) 11.637 1.687 ( 1.275 2.209 0.000) 2.551 1.690 ( 1.509 2.613 0.000) 3.018 1.781 ( 1.280 2.218 0.000) 2.561 2.010 ( -1.737 -3.008 0.000) 3.474 2.034 ( -1.901 -3.293 0.000) 3.803 2.067 ( -0.828 -1.433 0.000) 1.655 2.140 ( 3.266 5.657 0.000) 6.533 2.142 ( 3.219 5.576 0.000) 6.438 2.166 ( 3.773 6.535 0.000) 7.546 2.986 ( -0.174 -0.302 0.000) 0.349 2.998 ( 0.191 0.331 0.000) 0.382 3.003 ( 0.159 0.275 0.000) 0.318 3.065 ( -0.150 -0.260 0.000) 0.300 3.077 ( -0.606 -1.050 0.000) 1.213 3.082 ( -0.296 -0.513 0.000) 0.592 3.278 ( 1.827 3.165 0.000) 3.655 3.327 ( 1.457 2.523 0.000) 2.914 3.332 ( 1.334 2.311 0.000) 2.669 3.715 ( -0.669 -1.159 0.000) 1.339 3.718 ( -0.420 -0.728 0.000) 0.840 3.793 ( -1.935 -3.352 0.000) 3.870 ======================= 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| 0.998 ( 4.167 5.203 -0.299) 6.673 1.024 ( 4.415 5.199 0.489) 6.838 1.071 ( 3.904 4.325 -0.189) 5.829 1.603 ( -2.203 7.337 0.060) 7.661 1.618 ( -2.511 7.654 -0.644) 8.081 1.628 ( -2.215 8.175 0.684) 8.497 1.730 ( 2.274 -0.230 -1.457) 2.710 1.741 ( 1.925 -0.412 1.484) 2.466 1.803 ( 0.438 -1.780 -0.189) 1.843 1.983 ( 1.316 2.685 0.170) 2.995 2.008 ( 1.205 3.425 -0.014) 3.631 2.048 ( 0.600 3.839 -0.142) 3.888 2.229 ( -0.606 5.167 0.035) 5.202 2.257 ( 1.216 4.125 0.809) 4.376 2.274 ( 0.453 5.034 -0.802) 5.118 2.977 ( -0.534 -0.142 0.165) 0.577 2.990 ( -0.954 -0.201 -0.391) 1.050 3.000 ( -0.449 -0.596 0.273) 0.795 3.049 ( 1.538 -2.904 -0.316) 3.302 3.067 ( 2.438 -3.608 0.631) 4.400 3.075 ( 2.065 -3.074 -0.359) 3.721 3.330 ( -1.386 6.045 -0.127) 6.203 3.366 ( -2.369 6.266 0.553) 6.722 3.374 ( -2.043 6.537 -0.409) 6.861 3.703 ( 2.446 -4.040 -0.651) 4.767 3.718 ( 2.892 -4.221 0.862) 5.189 3.748 ( 2.008 -6.216 -0.226) 6.536 ======================= 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.096 ( 2.052 2.725 -0.044) 3.411 1.125 ( 2.503 2.590 0.105) 3.603 1.145 ( 1.065 1.271 -0.061) 1.659 1.608 ( -2.949 4.889 -0.144) 5.712 1.621 ( -3.396 4.802 -0.179) 5.884 1.627 ( -3.573 4.421 0.364) 5.696 1.740 ( 0.074 -1.226 0.423) 1.299 1.749 ( 0.922 -0.491 -0.833) 1.336 1.776 ( -0.607 -1.699 0.375) 1.843 2.078 ( 2.475 6.220 -0.303) 6.701 2.105 ( 1.996 6.392 0.516) 6.716 2.138 ( 1.541 6.377 -0.375) 6.571 2.221 ( -5.944 6.849 0.286) 9.073 2.256 ( -5.395 4.887 0.712) 7.314 2.271 ( -5.972 5.357 -0.848) 8.067 2.974 ( 0.686 -0.680 0.199) 0.987 2.979 ( 0.301 -0.851 -0.209) 0.927 2.987 ( 0.106 -1.203 -0.000) 1.207 3.060 ( 3.274 -4.289 -0.097) 5.396 3.077 ( 3.319 -4.858 0.475) 5.903 3.081 ( 2.982 -4.968 -0.364) 5.805 3.354 ( -2.855 6.443 -0.289) 7.053 3.371 ( -4.064 6.393 0.538) 7.594 3.389 ( -3.620 6.803 -0.241) 7.710 3.711 ( 3.659 -5.413 -0.121) 6.535 3.725 ( 3.428 -5.736 0.165) 6.684 3.733 ( 3.643 -6.791 -0.051) 7.707 ======================= 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.091 ( 2.331 4.037 0.000) 4.661 1.116 ( 2.195 3.803 0.000) 4.391 1.141 ( 1.442 2.497 0.000) 2.883 1.664 ( 0.153 0.265 0.000) 0.305 1.675 ( -0.192 -0.333 0.000) 0.385 1.686 ( 0.266 0.460 0.000) 0.531 1.725 ( -0.548 -0.949 0.000) 1.095 1.725 ( -0.787 -1.363 0.000) 1.574 1.764 ( -1.191 -2.062 0.000) 2.381 2.095 ( 4.751 8.230 0.000) 9.503 2.128 ( 4.679 8.104 0.000) 9.358 2.161 ( 4.435 7.681 0.000) 8.870 2.327 ( 1.142 1.978 0.000) 2.284 2.329 ( 1.346 2.331 0.000) 2.691 2.362 ( 1.299 2.250 0.000) 2.599 2.966 ( -0.554 -0.959 0.000) 1.107 2.976 ( -0.724 -1.254 0.000) 1.449 2.979 ( -0.855 -1.481 0.000) 1.710 3.001 ( -0.992 -1.718 0.000) 1.984 3.005 ( -1.032 -1.787 0.000) 2.063 3.020 ( -1.264 -2.189 0.000) 2.528 3.425 ( 2.113 3.660 0.000) 4.226 3.460 ( 1.845 3.195 0.000) 3.690 3.472 ( 1.998 3.461 0.000) 3.996 3.634 ( -1.546 -2.677 0.000) 3.091 3.643 ( -2.131 -3.691 0.000) 4.262 3.655 ( -1.529 -2.648 0.000) 3.057 ======================= 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.137 ( -0.806 1.396 0.000) 1.612 1.155 ( -0.114 0.198 -0.000) 0.229 1.163 ( -0.441 0.763 0.000) 0.881 1.667 ( -0.760 1.316 -0.000) 1.520 1.672 ( -0.835 1.446 0.000) 1.670 1.675 ( -1.322 2.290 -0.000) 2.644 1.705 ( 1.087 -1.882 0.000) 2.173 1.727 ( 1.436 -2.487 -0.000) 2.872 1.729 ( 0.912 -1.580 0.000) 1.825 2.222 ( -3.616 6.264 0.000) 7.233 2.228 ( -3.042 5.268 -0.000) 6.083 2.276 ( -3.441 5.960 -0.000) 6.882 2.316 ( -0.354 0.613 0.000) 0.708 2.339 ( -1.541 2.670 -0.000) 3.083 2.352 ( -1.065 1.844 0.000) 2.129 2.955 ( 0.510 -0.884 -0.000) 1.021 2.959 ( 0.486 -0.842 0.000) 0.972 2.961 ( 0.581 -1.006 -0.000) 1.161 2.986 ( 1.477 -2.559 -0.000) 2.955 2.991 ( 1.621 -2.807 0.000) 3.241 2.996 ( 1.651 -2.860 -0.000) 3.303 3.470 ( -2.983 5.166 0.000) 5.965 3.483 ( -2.894 5.012 -0.000) 5.788 3.508 ( -2.991 5.181 0.000) 5.983 3.602 ( 3.314 -5.741 0.000) 6.629 3.611 ( 2.637 -4.568 -0.000) 5.274 3.627 ( 2.680 -4.643 0.000) 5.361 ======================= 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.184 ( 0.000 0.000 0.000) 0.000 0.184 ( 0.000 -0.000 0.000) 0.000 0.209 ( 0.000 0.000 -0.000) 0.000 0.209 ( 0.000 -0.000 -0.000) 0.000 0.365 ( -0.000 0.000 0.000) 0.000 0.367 ( 0.000 0.000 -0.000) 0.000 0.398 ( 0.000 0.000 0.000) 0.000 0.398 ( 0.000 -0.000 0.000) 0.000 0.751 ( 0.000 0.000 0.000) 0.000 1.525 ( 0.000 0.000 0.000) 0.000 1.525 ( 0.000 -0.000 0.000) 0.000 1.593 ( 0.000 0.000 -0.000) 0.000 1.593 ( 0.000 -0.000 -0.000) 0.000 1.612 ( 0.000 0.000 0.000) 0.000 1.612 ( -0.000 -0.000 0.000) 0.000 2.584 ( 0.000 0.000 0.000) 0.000 2.584 ( 0.000 -0.000 0.000) 0.000 2.601 ( 0.000 0.000 -0.000) 0.000 2.601 ( 0.000 -0.000 -0.000) 0.000 2.613 ( 0.000 0.000 0.000) 0.000 2.613 ( 0.000 -0.000 0.000) 0.000 3.126 ( 0.000 -0.000 0.000) 0.000 3.191 ( -0.000 -0.000 0.000) 0.000 3.194 ( -0.000 -0.000 0.000) 0.000 4.399 ( 0.000 -0.000 0.000) 0.000 4.401 ( 0.000 -0.000 -0.000) 0.000 4.409 ( 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.467 6.043 -1.708) 12.206 0.347 ( 1.944 1.122 4.918) 5.406 0.350 ( 9.714 5.608 1.585) 11.328 0.357 ( 2.911 1.680 -4.945) 5.979 0.490 ( 7.178 4.144 0.071) 8.289 0.571 ( 20.143 11.630 2.011) 23.346 0.577 ( 19.846 11.458 -1.981) 23.001 0.648 ( 17.115 9.881 -0.226) 19.763 0.707 ( -2.664 -1.538 0.288) 3.089 1.536 ( 0.878 0.507 -0.042) 1.014 1.606 ( 0.994 0.574 0.033) 1.149 1.620 ( 0.688 0.397 0.006) 0.794 1.685 ( 12.580 7.263 -0.028) 14.526 1.733 ( 10.653 6.151 0.689) 12.320 1.749 ( 11.048 6.379 -0.776) 12.781 2.598 ( 1.254 0.724 -0.296) 1.477 2.610 ( 0.962 0.556 0.431) 1.192 2.631 ( 1.407 0.812 -0.137) 1.630 2.770 ( 11.218 6.477 0.004) 12.953 3.066 ( 14.318 8.266 14.169) 21.774 3.068 ( 14.050 8.112 -14.262) 21.600 3.140 ( 1.983 1.145 0.007) 2.289 3.221 ( 9.288 5.362 1.271) 10.799 3.228 ( 9.760 5.635 -1.568) 11.378 4.241 ( -11.777 -6.799 4.167) 14.223 4.250 ( -10.709 -6.183 -3.890) 12.964 4.301 ( -9.011 -5.202 0.072) 10.405 ======================= 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.427 ( 5.436 3.138 2.766) 6.859 0.455 ( 5.836 3.370 -3.064) 7.403 0.579 ( 10.095 5.828 2.373) 11.895 0.585 ( 10.571 6.103 -2.618) 12.484 0.662 ( -0.595 -0.344 0.192) 0.713 0.678 ( 8.836 5.102 0.249) 10.206 1.021 ( 18.348 10.593 0.629) 21.196 1.031 ( 19.211 11.091 -0.380) 22.186 1.050 ( 17.756 10.251 -0.114) 20.503 1.558 ( 0.920 0.531 -0.243) 1.090 1.627 ( 0.645 0.372 1.311) 1.508 1.639 ( 0.947 0.547 -1.076) 1.534 2.014 ( 15.298 8.832 0.450) 17.670 2.022 ( 14.119 8.152 -0.409) 16.308 2.059 ( 15.118 8.728 -0.034) 17.457 2.639 ( 2.332 1.347 0.044) 2.693 2.649 ( 2.482 1.433 -0.063) 2.867 2.674 ( 2.380 1.374 0.019) 2.749 2.950 ( 3.226 1.863 -0.168) 3.729 3.087 ( -2.409 -1.391 1.776) 3.301 3.092 ( -2.237 -1.291 -1.630) 3.054 3.277 ( 9.715 5.609 -0.091) 11.218 3.486 ( 5.203 3.004 10.633) 12.213 3.522 ( 7.649 4.416 -12.113) 14.991 3.943 ( -13.952 -8.055 4.426) 16.707 4.002 ( -9.946 -5.742 -4.112) 12.199 4.063 ( -9.455 -5.459 1.205) 10.984 ======================= 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.577 ( 7.479 4.318 0.970) 8.690 0.604 ( 6.943 4.008 -1.601) 8.175 0.685 ( 2.811 1.623 0.513) 3.286 0.792 ( 8.676 5.009 1.501) 10.130 0.807 ( 8.955 5.170 -1.887) 10.511 0.870 ( 7.996 4.617 0.396) 9.242 1.374 ( 13.190 7.615 0.304) 15.234 1.422 ( 15.174 8.761 -0.336) 17.525 1.431 ( 15.215 8.784 0.150) 17.569 1.571 ( 0.147 0.085 -0.508) 0.536 1.631 ( -0.268 -0.155 1.658) 1.687 1.655 ( 0.361 0.208 -1.160) 1.232 2.322 ( 11.887 6.863 -0.902) 13.755 2.345 ( 12.989 7.499 1.535) 15.077 2.387 ( 13.191 7.616 -0.639) 15.245 2.701 ( 3.067 1.771 0.301) 3.554 2.717 ( 3.433 1.982 -0.518) 3.998 2.736 ( 3.078 1.777 0.219) 3.561 2.982 ( 0.801 0.462 -0.339) 0.985 3.031 ( -2.250 -1.299 1.375) 2.940 3.041 ( -1.990 -1.149 -1.080) 2.539 3.378 ( -4.959 -2.863 -0.907) 5.797 3.428 ( -8.920 -5.150 3.008) 10.730 3.523 ( -8.883 -5.128 -4.631) 11.254 3.687 ( -1.261 -0.728 4.969) 5.178 3.900 ( 0.919 0.531 -6.875) 6.957 3.970 ( 0.202 0.117 4.434) 4.440 ======================= 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.742 ( 5.633 3.252 0.603) 6.532 0.754 ( 8.076 4.663 -1.134) 9.394 0.808 ( 7.578 4.375 0.458) 8.762 0.967 ( 6.795 3.923 0.957) 7.905 0.986 ( 6.826 3.941 -1.384) 8.002 1.032 ( 6.255 3.611 0.435) 7.236 1.564 ( -0.695 -0.401 -0.755) 1.102 1.617 ( -0.996 -0.575 1.661) 2.020 1.623 ( 8.991 5.191 0.430) 10.391 1.655 ( -0.280 -0.162 -0.917) 0.972 1.708 ( 10.095 5.828 -2.040) 11.834 1.719 ( 10.221 5.901 1.698) 11.924 2.510 ( 3.535 2.041 -1.113) 4.231 2.554 ( 4.053 2.340 2.240) 5.188 2.614 ( 5.495 3.173 -1.193) 6.457 2.772 ( 3.090 1.784 0.249) 3.577 2.796 ( 3.357 1.938 -0.699) 3.939 2.807 ( 3.096 1.788 0.451) 3.604 2.990 ( -0.483 -0.279 0.395) 0.683 3.010 ( 0.862 0.497 -0.330) 1.048 3.013 ( 0.274 0.158 -0.098) 0.332 3.217 ( -6.857 -3.959 0.461) 7.931 3.225 ( -5.085 -2.936 -0.592) 5.901 3.258 ( -9.324 -5.383 0.057) 10.766 3.789 ( 5.159 2.979 1.677) 6.189 3.944 ( 1.990 1.149 -5.120) 5.612 3.997 ( 1.289 0.744 3.574) 3.872 ======================= 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.885 ( 7.045 4.067 1.073) 8.205 0.936 ( 7.842 4.528 -2.236) 9.328 0.985 ( 7.624 4.402 1.142) 8.877 1.097 ( 4.740 2.737 0.652) 5.512 1.115 ( 4.687 2.706 -1.089) 5.520 1.150 ( 4.105 2.370 0.461) 4.763 1.544 ( -0.977 -0.564 -1.064) 1.551 1.589 ( -1.328 -0.767 1.695) 2.286 1.647 ( -0.377 -0.218 -0.656) 0.787 1.788 ( 5.938 3.428 0.417) 6.869 1.881 ( 5.446 3.144 -2.531) 6.779 1.894 ( 5.503 3.177 2.145) 6.707 2.461 ( -6.558 -3.786 -1.541) 7.727 2.499 ( -7.322 -4.227 2.799) 8.906 2.574 ( -7.662 -4.423 -1.313) 8.944 2.835 ( 2.368 1.367 0.145) 2.738 2.863 ( 2.533 1.463 -0.748) 3.019 2.870 ( 2.346 1.354 0.602) 2.775 3.026 ( 3.263 1.884 -0.168) 3.772 3.048 ( 2.998 1.731 0.804) 3.554 3.064 ( 3.146 1.816 -0.788) 3.717 3.160 ( 1.088 0.628 0.330) 1.299 3.174 ( 1.144 0.661 -0.051) 1.323 3.201 ( 1.558 0.899 -0.105) 1.802 3.878 ( 2.748 1.586 1.062) 3.346 3.977 ( 0.951 0.549 -3.431) 3.602 4.015 ( 0.281 0.162 2.390) 2.412 ======================= 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.043 ( 6.282 3.627 0.927) 7.313 1.092 ( 5.464 3.155 -2.031) 6.629 1.129 ( 4.593 2.652 1.106) 5.418 1.179 ( 2.506 1.447 0.344) 2.914 1.197 ( 2.593 1.497 -0.869) 3.118 1.217 ( 1.826 1.054 0.577) 2.186 1.525 ( -0.615 -0.355 -1.194) 1.389 1.561 ( -1.058 -0.611 1.482) 1.921 1.641 ( -0.176 -0.102 -0.340) 0.396 1.901 ( 4.178 2.412 0.296) 4.834 1.968 ( 2.562 1.479 -1.692) 3.408 1.981 ( 2.480 1.432 1.349) 3.165 2.295 ( -6.847 -3.953 -1.327) 8.017 2.316 ( -7.394 -4.269 2.256) 8.831 2.384 ( -7.484 -4.321 -0.887) 8.687 2.876 ( 1.241 0.717 0.058) 1.434 2.907 ( 1.345 0.777 -0.736) 1.719 2.910 ( 1.195 0.690 0.677) 1.537 3.096 ( 2.551 1.473 -0.300) 2.961 3.112 ( 2.245 1.296 0.645) 2.671 3.127 ( 1.972 1.139 -0.371) 2.307 3.213 ( 2.428 1.402 0.042) 2.804 3.227 ( 2.532 1.462 -0.008) 2.924 3.248 ( 1.876 1.083 -0.011) 2.166 3.920 ( 1.121 0.647 0.597) 1.426 3.991 ( 0.385 0.222 -2.309) 2.352 4.012 ( -0.449 -0.259 1.724) 1.800 ======================= 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.128 ( -0.000 -0.000 -0.000) 0.000 1.156 ( 0.000 0.000 0.000) 0.000 1.178 ( 0.000 0.000 0.000) 0.000 1.208 ( 0.000 0.000 0.000) 0.000 1.230 ( 0.000 0.000 -0.000) 0.000 1.233 ( -0.000 -0.000 0.000) 0.000 1.518 ( -0.000 -0.000 0.000) 0.000 1.548 ( 0.000 0.000 0.000) 0.000 1.639 ( 0.000 0.000 0.000) 0.000 1.962 ( 0.000 0.000 0.000) 0.000 2.001 ( 0.000 0.000 0.000) 0.000 2.007 ( 0.000 0.000 -0.000) 0.000 2.202 ( 0.000 0.000 0.000) 0.000 2.222 ( 0.000 0.000 0.000) 0.000 2.290 ( 0.000 0.000 0.000) 0.000 2.890 ( 0.000 0.000 0.000) 0.000 2.923 ( 0.000 0.000 -0.000) 0.000 2.923 ( -0.000 -0.000 0.000) 0.000 3.128 ( 0.000 0.000 0.000) 0.000 3.140 ( 0.000 0.000 0.000) 0.000 3.147 ( -0.000 -0.000 0.000) 0.000 3.242 ( 0.000 0.000 -0.000) 0.000 3.261 ( 0.000 0.000 0.000) 0.000 3.268 ( 0.000 0.000 0.000) 0.000 3.932 ( 0.000 0.000 0.000) 0.000 3.999 ( 0.000 0.000 -0.000) 0.000 4.002 ( -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.395 ( 2.496 4.324 0.000) 4.993 0.416 ( 2.718 4.707 -0.000) 5.435 0.539 ( 6.991 12.108 0.000) 13.981 0.541 ( 6.918 11.983 -0.000) 13.836 0.635 ( 5.196 9.000 0.000) 10.392 0.667 ( -0.394 -0.682 0.000) 0.788 0.899 ( 11.101 19.227 0.000) 22.202 0.905 ( 11.022 19.091 -0.000) 22.045 0.927 ( 9.629 16.678 0.000) 19.258 1.599 ( 2.380 4.122 0.000) 4.759 1.664 ( 1.910 3.308 0.000) 3.819 1.666 ( 2.292 3.970 0.000) 4.584 1.881 ( 7.057 12.222 0.000) 14.113 1.906 ( 6.316 10.939 0.000) 12.631 1.909 ( 6.013 10.415 0.000) 12.026 2.634 ( 1.721 2.980 0.000) 3.441 2.647 ( 1.668 2.890 -0.000) 3.337 2.666 ( 1.887 3.268 0.000) 3.774 2.923 ( 3.307 5.729 0.000) 6.615 3.101 ( -1.118 -1.937 -0.000) 2.236 3.104 ( -1.049 -1.817 0.000) 2.098 3.222 ( 4.978 8.623 0.000) 9.957 3.457 ( 6.399 11.084 -0.000) 12.798 3.458 ( 6.723 11.645 0.000) 13.446 4.048 ( -6.897 -11.946 0.000) 13.794 4.051 ( -6.804 -11.785 -0.000) 13.609 4.110 ( -7.711 -13.356 0.000) 15.423 ======================= 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.513 ( 5.386 5.632 1.029) 7.861 0.538 ( 5.540 5.217 -1.337) 7.726 0.654 ( 1.417 0.148 0.239) 1.445 0.773 ( 5.616 12.716 1.627) 13.996 0.780 ( 6.213 12.530 -1.802) 14.101 0.843 ( 5.675 10.931 0.177) 12.317 1.243 ( 11.432 11.057 0.176) 15.905 1.270 ( 13.410 12.509 0.159) 18.339 1.276 ( 13.626 12.210 -0.269) 18.298 1.658 ( -2.821 8.536 -0.257) 8.993 1.713 ( -2.759 7.354 1.348) 7.969 1.725 ( -2.208 7.627 -1.084) 8.014 2.125 ( 12.299 2.187 -0.478) 12.501 2.133 ( 12.992 2.808 0.658) 13.308 2.158 ( 14.120 1.449 -0.191) 14.196 2.716 ( 2.520 6.452 0.088) 6.928 2.727 ( 2.447 6.507 -0.135) 6.954 2.752 ( 2.328 6.707 0.059) 7.100 2.974 ( 0.724 0.667 -0.159) 0.997 3.051 ( -1.908 -2.081 0.980) 2.989 3.058 ( -1.721 -2.025 -0.854) 2.791 3.418 ( 3.207 7.784 -0.426) 8.430 3.526 ( -4.783 -0.127 4.344) 6.463 3.569 ( -0.524 -1.393 -4.698) 4.928 3.772 ( -13.524 -8.906 2.154) 16.336 3.868 ( -1.478 -6.271 -4.223) 7.704 3.915 ( 1.590 -8.012 2.845) 8.650 ======================= 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.673 ( 6.771 5.155 -0.607) 8.531 0.682 ( 5.230 4.112 0.199) 6.656 0.724 ( 6.055 3.023 0.339) 6.776 0.973 ( 3.254 12.618 1.084) 13.076 0.986 ( 3.661 12.345 -1.409) 12.953 1.030 ( 3.515 11.103 0.317) 11.651 1.509 ( 10.044 5.935 0.308) 11.671 1.579 ( 11.273 6.881 -1.191) 13.261 1.588 ( 11.673 6.878 0.956) 13.582 1.673 ( -5.299 9.497 -0.497) 10.887 1.716 ( -5.497 8.143 1.450) 9.931 1.742 ( -4.586 8.060 -0.945) 9.321 2.325 ( 11.985 -5.020 -0.628) 13.009 2.347 ( 12.909 -5.586 1.146) 14.112 2.375 ( 13.571 -6.359 -0.540) 14.997 2.841 ( 2.258 10.058 0.134) 10.309 2.848 ( 2.036 9.647 -0.321) 9.865 2.866 ( 1.242 9.589 0.189) 9.671 2.992 ( 0.843 0.653 -0.304) 1.109 3.006 ( -1.211 -1.034 0.714) 1.745 3.018 ( -0.964 -1.022 -0.436) 1.471 3.317 ( -9.438 -3.081 -1.119) 9.991 3.337 ( -11.392 -4.194 1.887) 12.285 3.417 ( -14.627 -4.985 -1.321) 15.510 3.694 ( 7.734 0.449 2.150) 8.039 3.881 ( 4.299 -2.308 -5.595) 7.424 3.930 ( 4.552 -3.743 4.006) 7.126 ======================= 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.806 ( 6.287 3.524 0.652) 7.236 0.848 ( 7.515 4.940 -1.670) 9.147 0.888 ( 8.029 4.019 0.987) 9.033 1.127 ( 0.930 11.760 0.714) 11.818 1.143 ( 1.085 11.427 -1.061) 11.527 1.174 ( 1.047 10.338 0.341) 10.396 1.648 ( -6.029 7.799 -0.652) 9.879 1.669 ( -0.675 3.721 1.113) 3.942 1.704 ( -2.169 5.461 -0.183) 5.879 1.743 ( 0.444 6.986 -0.029) 7.000 1.804 ( 6.861 4.398 -2.015) 8.395 1.812 ( 7.708 3.845 1.813) 8.803 2.400 ( 4.556 -9.264 -0.857) 10.360 2.430 ( 5.250 -9.805 1.586) 11.235 2.468 ( 6.587 -10.981 -0.724) 12.826 2.933 ( -3.081 7.561 -0.244) 8.168 2.940 ( -3.046 6.220 0.311) 6.932 2.955 ( -3.022 7.381 -0.123) 7.977 3.007 ( 2.605 4.342 -0.208) 5.068 3.028 ( 2.210 3.505 0.931) 4.247 3.038 ( 2.519 2.959 -0.752) 3.958 3.167 ( -3.708 -1.455 -0.045) 3.983 3.186 ( -3.922 -1.366 0.081) 4.154 3.198 ( -3.940 -1.270 -0.021) 4.140 3.808 ( 5.654 -0.875 1.167) 5.840 3.924 ( 3.679 -3.055 -3.836) 6.131 3.964 ( 3.518 -3.855 2.714) 5.883 ======================= 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.959 ( 7.156 3.588 0.793) 8.044 1.016 ( 6.909 3.684 -2.114) 8.110 1.052 ( 6.841 2.734 1.319) 7.484 1.233 ( -1.269 10.351 0.409) 10.436 1.247 ( -1.112 10.011 -0.781) 10.102 1.266 ( -1.299 8.949 0.384) 9.051 1.621 ( -5.598 7.292 -0.921) 9.239 1.645 ( -5.001 5.339 1.364) 7.442 1.696 ( -3.225 4.043 -0.453) 5.192 1.860 ( 4.210 4.532 0.275) 6.192 1.928 ( 4.098 2.122 -1.931) 5.003 1.940 ( 4.346 1.903 1.652) 5.024 2.316 ( -2.704 -8.503 -0.974) 8.976 2.347 ( -3.157 -8.522 1.975) 9.301 2.386 ( -2.102 -10.852 -0.980) 11.097 2.933 ( -2.756 4.822 0.145) 5.556 2.937 ( -1.432 3.448 -0.298) 3.745 2.951 ( -2.402 4.280 0.121) 4.909 3.093 ( 0.731 5.339 -0.040) 5.389 3.104 ( 1.186 4.257 0.143) 4.422 3.120 ( 1.288 4.603 -0.235) 4.786 3.166 ( 3.238 -0.323 0.353) 3.273 3.181 ( 3.317 -0.164 -0.203) 3.327 3.208 ( 2.829 -0.376 -0.015) 2.853 3.871 ( 3.688 -2.107 0.781) 4.319 3.946 ( 2.879 -3.597 -2.563) 5.272 3.973 ( 2.363 -4.195 1.797) 5.139 ======================= 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.091 ( 4.458 1.566 0.383) 4.740 1.132 ( 3.262 1.080 -1.170) 3.630 1.150 ( 2.665 -0.033 0.792) 2.780 1.285 ( -3.462 8.598 0.078) 9.269 1.300 ( -3.087 8.180 -0.564) 8.761 1.306 ( -3.399 7.481 0.500) 8.232 1.607 ( -4.586 7.508 -0.941) 8.849 1.620 ( -4.174 5.618 1.085) 7.083 1.686 ( -2.411 3.754 -0.157) 4.464 1.961 ( 2.224 4.235 0.258) 4.791 1.994 ( 1.619 1.645 -1.184) 2.594 2.004 ( 1.181 1.487 0.844) 2.078 2.194 ( -2.412 -5.290 -0.353) 5.825 2.214 ( -2.657 -5.167 1.032) 5.901 2.253 ( -0.983 -7.549 -0.597) 7.636 2.939 ( -1.992 4.033 0.087) 4.499 2.951 ( -0.841 2.789 -0.378) 2.937 2.956 ( -1.499 2.886 0.286) 3.264 3.141 ( -0.261 3.304 -0.214) 3.321 3.148 ( 0.020 2.460 0.265) 2.474 3.162 ( -0.473 2.781 -0.022) 2.821 3.214 ( 2.428 -1.352 -0.048) 2.780 3.234 ( 2.558 -0.727 -0.024) 2.659 3.244 ( 1.936 -1.517 0.039) 2.460 3.898 ( 2.367 -2.772 0.302) 3.658 3.956 ( 2.622 -3.720 -1.632) 4.835 3.966 ( 1.832 -4.222 1.337) 4.792 ======================= 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.639 ( 3.907 6.768 0.000) 7.815 0.657 ( 3.963 6.864 -0.000) 7.926 0.683 ( 1.926 3.336 0.000) 3.853 1.030 ( 7.286 12.619 0.000) 14.571 1.032 ( 7.303 12.649 -0.000) 14.606 1.070 ( 6.679 11.568 0.000) 13.357 1.441 ( 5.296 9.173 0.000) 10.592 1.502 ( 6.290 10.894 -0.000) 12.580 1.505 ( 6.305 10.921 0.000) 12.611 1.830 ( 4.226 7.319 0.000) 8.451 1.863 ( 3.657 6.334 0.000) 7.314 1.883 ( 3.865 6.694 0.000) 7.729 2.140 ( 0.480 0.831 -0.000) 0.960 2.148 ( 0.168 0.292 0.000) 0.337 2.166 ( 0.528 0.914 0.000) 1.056 2.866 ( 4.706 8.151 0.000) 9.412 2.877 ( 4.595 7.959 0.000) 9.190 2.903 ( 4.499 7.792 0.000) 8.998 2.987 ( 0.497 0.860 0.000) 0.993 3.014 ( -0.867 -1.501 -0.000) 1.734 3.021 ( -0.869 -1.504 0.000) 1.737 3.427 ( -4.783 -8.284 0.000) 9.565 3.439 ( -4.798 -8.310 -0.000) 9.596 3.566 ( -6.944 -12.027 0.000) 13.887 3.585 ( 3.558 6.162 -0.000) 7.116 3.821 ( -0.156 -0.270 -0.000) 0.312 3.823 ( -0.278 -0.481 0.000) 0.556 ======================= 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.754 ( 4.218 4.466 0.105) 6.144 0.794 ( 5.039 6.484 -0.909) 8.262 0.810 ( 6.223 5.474 0.786) 8.325 1.244 ( 2.379 14.238 0.532) 14.445 1.252 ( 2.806 14.027 -0.728) 14.324 1.272 ( 2.676 12.935 0.186) 13.210 1.603 ( 6.336 3.665 0.209) 7.323 1.686 ( 6.508 4.108 -1.651) 7.871 1.694 ( 7.054 3.989 1.481) 8.238 1.896 ( -7.068 12.109 -0.214) 14.023 1.911 ( -7.206 10.633 0.699) 12.864 1.932 ( -6.801 10.152 -0.505) 12.230 2.194 ( 8.670 -7.327 -0.573) 11.366 2.206 ( 9.792 -7.340 0.822) 12.265 2.223 ( 9.721 -8.067 -0.244) 12.635 2.984 ( 0.028 0.872 0.404) 0.961 2.987 ( 0.202 0.741 -0.393) 0.863 3.003 ( 0.059 0.521 -0.171) 0.552 3.012 ( 3.545 5.979 0.017) 6.951 3.036 ( 2.774 5.040 0.263) 5.759 3.043 ( 2.668 5.562 -0.132) 6.170 3.267 ( -6.742 -1.407 -0.294) 6.894 3.276 ( -7.273 -1.463 0.522) 7.437 3.344 ( -9.647 -2.009 -0.311) 9.859 3.694 ( 5.494 -0.343 0.728) 5.553 3.823 ( 2.697 -3.207 -3.524) 5.475 3.845 ( 4.001 -4.550 2.887) 6.712 ======================= 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.887 ( 5.873 4.825 0.407) 7.612 0.947 ( 5.958 5.145 -1.722) 8.058 0.970 ( 6.559 4.360 1.306) 7.984 1.386 ( -1.659 13.624 0.218) 13.726 1.396 ( -1.681 13.576 -0.409) 13.685 1.405 ( -1.690 12.463 0.182) 12.578 1.710 ( 2.977 0.865 0.408) 3.127 1.774 ( -0.716 2.367 -1.239) 2.766 1.788 ( -0.482 1.451 0.883) 1.766 1.896 ( -2.601 7.772 -0.153) 8.197 1.909 ( -0.397 7.333 0.114) 7.345 1.931 ( -2.162 7.952 0.009) 8.241 2.236 ( 3.183 -5.895 -0.295) 6.706 2.264 ( 3.883 -6.397 0.284) 7.488 2.281 ( 4.001 -6.130 0.018) 7.320 2.975 ( -0.884 0.391 -0.158) 0.979 2.986 ( -0.585 0.367 0.334) 0.767 3.004 ( -0.866 0.473 -0.269) 1.023 3.096 ( 3.206 -0.594 0.303) 3.275 3.108 ( 2.520 -0.656 -0.214) 2.613 3.121 ( 3.239 0.552 -0.002) 3.285 3.201 ( -3.992 7.088 0.576) 8.155 3.213 ( -3.118 5.821 -0.704) 6.640 3.241 ( -5.465 7.285 0.099) 9.108 3.765 ( 5.058 -3.210 0.692) 6.030 3.838 ( 3.867 -5.357 -2.373) 7.020 3.862 ( 4.127 -6.140 1.705) 7.592 ======================= 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.030 ( 5.409 3.684 0.350) 6.554 1.082 ( 4.643 3.032 -1.436) 5.728 1.100 ( 4.333 2.198 1.085) 4.979 1.456 ( -4.106 11.067 -0.118) 11.804 1.461 ( -4.116 10.020 0.183) 10.834 1.467 ( -4.416 11.200 -0.067) 12.039 1.727 ( -1.334 1.820 0.237) 2.269 1.753 ( -3.884 3.476 0.274) 5.220 1.760 ( -2.891 2.535 -0.454) 3.872 1.972 ( 2.502 5.752 -0.537) 6.296 1.997 ( 2.658 4.667 0.131) 5.373 2.006 ( 1.441 6.941 0.392) 7.099 2.191 ( -2.638 -3.022 0.132) 4.013 2.213 ( -2.528 -5.090 -0.882) 5.751 2.234 ( -3.976 -2.202 0.723) 4.602 2.974 ( 0.186 0.738 -0.235) 0.797 2.985 ( -0.234 0.834 0.373) 0.943 2.994 ( -0.629 0.453 -0.149) 0.789 3.136 ( 3.659 -2.575 0.368) 4.489 3.142 ( 3.297 -2.348 -0.357) 4.063 3.167 ( 3.505 -2.763 -0.044) 4.463 3.222 ( -2.783 6.955 0.365) 7.500 3.236 ( -2.322 6.518 -0.488) 6.936 3.252 ( -3.047 7.243 0.145) 7.859 3.805 ( 4.102 -4.323 0.472) 5.978 3.848 ( 3.973 -6.057 -1.433) 7.384 3.864 ( 3.586 -6.550 0.974) 7.531 ======================= 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.107 ( 0.017 -0.030 0.000) 0.034 1.140 ( 0.069 -0.119 0.000) 0.138 1.144 ( 0.201 -0.348 0.000) 0.402 1.471 ( -5.002 8.663 0.000) 10.004 1.477 ( -5.629 9.749 -0.000) 11.258 1.478 ( -5.126 8.879 0.000) 10.253 1.723 ( -2.148 3.720 0.000) 4.296 1.736 ( -1.118 1.937 0.000) 2.236 1.750 ( -2.259 3.913 -0.000) 4.519 2.052 ( -2.754 4.771 0.000) 5.509 2.064 ( -2.305 3.993 0.000) 4.611 2.082 ( -4.197 7.269 -0.000) 8.394 2.128 ( 0.056 -0.097 -0.000) 0.112 2.132 ( 2.037 -3.528 0.000) 4.074 2.155 ( -0.289 0.501 -0.000) 0.579 2.985 ( -0.249 0.432 0.000) 0.499 2.988 ( -0.254 0.439 0.000) 0.507 2.989 ( -0.365 0.632 0.000) 0.730 3.159 ( 2.355 -4.079 -0.000) 4.710 3.160 ( 1.867 -3.234 0.000) 3.734 3.183 ( 2.351 -4.072 -0.000) 4.702 3.232 ( -3.350 5.803 0.000) 6.700 3.249 ( -3.239 5.610 0.000) 6.478 3.261 ( -3.112 5.390 0.000) 6.224 3.819 ( 2.892 -5.009 0.000) 5.784 3.854 ( 3.695 -6.400 -0.000) 7.391 3.859 ( 3.629 -6.286 0.000) 7.258 ======================= 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.860 ( 3.638 6.302 0.000) 7.276 0.924 ( 3.767 6.525 -0.000) 7.535 0.930 ( 3.856 6.679 0.000) 7.713 1.493 ( 5.900 10.219 0.000) 11.799 1.498 ( 5.725 9.915 0.000) 11.449 1.508 ( 6.057 10.491 0.000) 12.114 1.660 ( 1.477 2.558 0.000) 2.953 1.745 ( 1.276 2.209 -0.000) 2.551 1.750 ( 1.235 2.140 0.000) 2.471 2.016 ( -1.803 -3.123 0.000) 3.607 2.031 ( -1.494 -2.588 0.000) 2.988 2.073 ( -1.881 -3.259 -0.000) 3.763 2.128 ( 4.103 7.106 -0.000) 8.205 2.139 ( 3.407 5.901 0.000) 6.814 2.174 ( 3.500 6.063 0.000) 7.001 2.986 ( -0.004 -0.007 0.000) 0.008 2.991 ( 0.028 0.049 -0.000) 0.057 3.013 ( 0.142 0.246 0.000) 0.285 3.059 ( -0.283 -0.489 0.000) 0.565 3.078 ( -0.329 -0.571 0.000) 0.659 3.084 ( -0.429 -0.743 0.000) 0.858 3.291 ( 1.753 3.036 0.000) 3.505 3.296 ( 1.670 2.893 0.000) 3.340 3.350 ( 1.168 2.022 0.000) 2.335 3.688 ( -0.011 -0.019 0.000) 0.022 3.767 ( -1.555 -2.694 -0.000) 3.110 3.770 ( -1.429 -2.475 0.000) 2.858 ======================= 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| 0.990 ( 4.442 5.504 0.150) 7.074 1.046 ( 3.856 4.845 -1.062) 6.283 1.057 ( 4.189 4.388 0.911) 6.135 1.605 ( -2.429 7.397 -0.174) 7.787 1.611 ( -2.162 7.385 0.391) 7.705 1.636 ( -2.532 8.885 -0.334) 9.244 1.711 ( 2.663 -0.614 0.352) 2.755 1.776 ( 1.163 -0.834 -1.661) 2.192 1.783 ( 1.092 -1.379 1.486) 2.303 1.986 ( 1.170 3.064 -0.200) 3.286 2.007 ( 1.044 2.809 0.047) 2.997 2.048 ( 0.865 3.853 0.141) 3.951 2.236 ( 0.255 4.709 0.585) 4.752 2.240 ( -0.170 5.232 -1.011) 5.331 2.282 ( 0.930 4.503 0.383) 4.614 2.978 ( -0.752 -0.043 -0.253) 0.795 2.986 ( -0.341 -0.287 0.386) 0.590 3.002 ( -0.978 -0.613 -0.181) 1.168 3.050 ( 2.189 -3.173 0.420) 3.878 3.063 ( 1.579 -2.912 -0.541) 3.356 3.078 ( 2.405 -3.491 0.167) 4.243 3.339 ( -1.809 6.151 0.580) 6.438 3.346 ( -1.477 6.189 -0.653) 6.396 3.385 ( -2.530 6.502 0.055) 6.977 3.697 ( 2.940 -3.460 0.305) 4.550 3.730 ( 1.999 -5.339 -0.883) 5.769 3.742 ( 2.429 -5.675 0.593) 6.202 ======================= 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.097 ( 2.462 2.986 0.044) 3.870 1.132 ( 1.743 1.936 -0.294) 2.621 1.138 ( 1.408 1.692 0.250) 2.215 1.607 ( -3.178 4.635 0.162) 5.622 1.615 ( -2.962 4.822 -0.042) 5.659 1.636 ( -3.987 4.852 -0.165) 6.282 1.730 ( 1.154 -0.972 -0.049) 1.509 1.757 ( -0.843 -1.327 0.542) 1.663 1.775 ( 0.286 -1.324 -0.453) 1.428 2.085 ( 2.672 6.400 0.666) 6.967 2.099 ( 2.252 6.315 -0.871) 6.761 2.137 ( 0.932 6.055 0.361) 6.137 2.234 ( -5.971 6.810 -1.178) 9.133 2.236 ( -5.029 5.674 0.631) 7.608 2.278 ( -6.148 4.814 0.403) 7.818 2.974 ( 0.653 -0.593 -0.135) 0.893 2.982 ( 0.478 -0.927 0.232) 1.069 2.984 ( -0.021 -1.217 -0.088) 1.220 3.064 ( 3.378 -4.431 0.295) 5.580 3.069 ( 2.990 -4.574 -0.400) 5.479 3.086 ( 3.186 -5.097 0.093) 6.012 3.354 ( -3.422 6.367 0.315) 7.235 3.369 ( -3.065 6.605 -0.508) 7.299 3.391 ( -4.050 6.648 0.185) 7.787 3.712 ( 3.556 -5.129 0.148) 6.243 3.724 ( 3.738 -6.234 -0.186) 7.272 3.733 ( 3.439 -6.568 0.045) 7.414 ======================= 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.090 ( 2.504 4.337 0.000) 5.008 1.125 ( 1.717 2.974 -0.000) 3.435 1.133 ( 1.764 3.055 0.000) 3.528 1.666 ( 0.057 0.099 -0.000) 0.114 1.672 ( 0.187 0.324 0.000) 0.374 1.692 ( -0.170 -0.295 0.000) 0.340 1.706 ( -0.417 -0.723 0.000) 0.834 1.748 ( -1.133 -1.963 -0.000) 2.266 1.754 ( -0.836 -1.448 0.000) 1.672 2.104 ( 4.821 8.350 0.000) 9.641 2.114 ( 4.514 7.819 -0.000) 9.029 2.166 ( 4.552 7.885 0.000) 9.105 2.318 ( 1.118 1.936 0.000) 2.236 2.338 ( 1.402 2.428 -0.000) 2.804 2.362 ( 1.239 2.146 0.000) 2.478 2.968 ( -0.608 -1.053 0.000) 1.215 2.973 ( -0.612 -1.059 0.000) 1.223 2.979 ( -0.926 -1.604 0.000) 1.852 2.998 ( -0.860 -1.490 0.000) 1.721 3.014 ( -1.292 -2.238 0.000) 2.584 3.014 ( -1.116 -1.934 -0.000) 2.233 3.434 ( 2.009 3.479 -0.000) 4.018 3.441 ( 2.064 3.575 0.000) 4.128 3.482 ( 1.876 3.249 0.000) 3.752 3.633 ( -1.837 -3.182 0.000) 3.675 3.647 ( -1.449 -2.510 0.000) 2.899 3.653 ( -1.915 -3.316 -0.000) 3.829 ======================= 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.140 ( -0.731 1.266 0.000) 1.462 1.153 ( -0.213 0.370 0.000) 0.427 1.164 ( -0.434 0.752 0.000) 0.869 1.664 ( -1.141 1.976 0.000) 2.282 1.674 ( -1.008 1.746 0.000) 2.016 1.676 ( -0.444 0.769 -0.000) 0.888 1.705 ( 0.695 -1.204 -0.000) 1.390 1.719 ( 1.167 -2.021 0.000) 2.334 1.735 ( 1.266 -2.193 -0.000) 2.532 2.217 ( -3.429 5.939 0.000) 6.857 2.240 ( -3.027 5.242 0.000) 6.053 2.263 ( -3.404 5.897 0.000) 6.809 2.320 ( -0.868 1.504 -0.000) 1.736 2.343 ( -0.902 1.562 0.000) 1.804 2.348 ( -1.423 2.465 0.000) 2.846 2.956 ( 0.536 -0.928 0.000) 1.071 2.959 ( 0.548 -0.949 -0.000) 1.096 2.960 ( 0.475 -0.823 0.000) 0.950 2.987 ( 1.610 -2.788 -0.000) 3.220 2.990 ( 1.395 -2.416 0.000) 2.789 2.996 ( 1.754 -3.038 0.000) 3.508 3.468 ( -3.006 5.206 0.000) 6.011 3.486 ( -2.920 5.058 0.000) 5.840 3.506 ( -2.845 4.928 0.000) 5.690 3.602 ( 3.033 -5.253 0.000) 6.065 3.615 ( 2.989 -5.178 0.000) 5.979 3.624 ( 2.517 -4.359 -0.000) 5.034 =================== 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 70.574 70.574 1.959 -0.000 -0.000 -0.000 3/10584 20.0 30.207 30.207 0.911 -0.000 -0.000 -0.000 3/10584 30.0 15.940 15.940 0.518 -0.000 -0.000 -0.000 3/10584 40.0 10.271 10.271 0.357 -0.000 -0.000 -0.000 3/10584 50.0 7.514 7.514 0.274 -0.000 -0.000 -0.000 3/10584 60.0 5.930 5.930 0.225 -0.000 -0.000 -0.000 3/10584 70.0 4.909 4.909 0.191 -0.000 -0.000 -0.000 3/10584 80.0 4.197 4.197 0.166 -0.000 -0.000 -0.000 3/10584 90.0 3.670 3.670 0.148 -0.000 -0.000 -0.000 3/10584 100.0 3.265 3.265 0.133 -0.000 -0.000 -0.000 3/10584 110.0 2.942 2.942 0.121 -0.000 -0.000 -0.000 3/10584 120.0 2.679 2.679 0.111 -0.000 -0.000 -0.000 3/10584 130.0 2.460 2.460 0.102 -0.000 -0.000 -0.000 3/10584 140.0 2.275 2.275 0.095 -0.000 -0.000 -0.000 3/10584 150.0 2.116 2.116 0.089 -0.000 -0.000 -0.000 3/10584 160.0 1.979 1.979 0.083 -0.000 -0.000 -0.000 3/10584 170.0 1.858 1.858 0.078 -0.000 -0.000 -0.000 3/10584 180.0 1.751 1.751 0.074 -0.000 -0.000 -0.000 3/10584 190.0 1.657 1.657 0.070 -0.000 -0.000 -0.000 3/10584 200.0 1.572 1.572 0.067 -0.000 -0.000 -0.000 3/10584 210.0 1.495 1.495 0.064 -0.000 -0.000 -0.000 3/10584 220.0 1.426 1.426 0.061 -0.000 -0.000 -0.000 3/10584 230.0 1.363 1.363 0.058 -0.000 -0.000 -0.000 3/10584 240.0 1.305 1.305 0.056 -0.000 -0.000 -0.000 3/10584 250.0 1.252 1.252 0.053 -0.000 -0.000 -0.000 3/10584 260.0 1.203 1.203 0.051 -0.000 -0.000 -0.000 3/10584 270.0 1.158 1.158 0.049 -0.000 -0.000 -0.000 3/10584 280.0 1.116 1.116 0.048 -0.000 -0.000 -0.000 3/10584 290.0 1.077 1.077 0.046 -0.000 -0.000 -0.000 3/10584 300.0 1.041 1.041 0.045 -0.000 -0.000 -0.000 3/10584 310.0 1.007 1.007 0.043 -0.000 -0.000 -0.000 3/10584 320.0 0.975 0.975 0.042 -0.000 -0.000 -0.000 3/10584 330.0 0.945 0.945 0.040 -0.000 -0.000 -0.000 3/10584 340.0 0.917 0.917 0.039 -0.000 -0.000 -0.000 3/10584 350.0 0.891 0.891 0.038 -0.000 -0.000 -0.000 3/10584 360.0 0.866 0.866 0.037 -0.000 -0.000 -0.000 3/10584 370.0 0.842 0.842 0.036 -0.000 -0.000 -0.000 3/10584 380.0 0.820 0.820 0.035 -0.000 -0.000 -0.000 3/10584 390.0 0.799 0.799 0.034 -0.000 -0.000 -0.000 3/10584 400.0 0.779 0.779 0.033 -0.000 -0.000 -0.000 3/10584 410.0 0.760 0.760 0.033 -0.000 -0.000 -0.000 3/10584 420.0 0.741 0.741 0.032 -0.000 -0.000 -0.000 3/10584 430.0 0.724 0.724 0.031 -0.000 -0.000 -0.000 3/10584 440.0 0.708 0.708 0.030 -0.000 -0.000 -0.000 3/10584 450.0 0.692 0.692 0.030 -0.000 -0.000 -0.000 3/10584 460.0 0.677 0.677 0.029 -0.000 -0.000 -0.000 3/10584 470.0 0.662 0.662 0.028 -0.000 -0.000 -0.000 3/10584 480.0 0.648 0.648 0.028 -0.000 -0.000 -0.000 3/10584 490.0 0.635 0.635 0.027 -0.000 -0.000 -0.000 3/10584 500.0 0.622 0.622 0.027 -0.000 -0.000 -0.000 3/10584 510.0 0.610 0.610 0.026 -0.000 -0.000 -0.000 3/10584 520.0 0.598 0.598 0.026 -0.000 -0.000 -0.000 3/10584 530.0 0.587 0.587 0.025 -0.000 -0.000 -0.000 3/10584 540.0 0.576 0.576 0.025 -0.000 -0.000 -0.000 3/10584 550.0 0.566 0.566 0.024 -0.000 -0.000 -0.000 3/10584 560.0 0.555 0.555 0.024 -0.000 -0.000 -0.000 3/10584 570.0 0.546 0.546 0.023 -0.000 -0.000 -0.000 3/10584 580.0 0.536 0.536 0.023 -0.000 -0.000 -0.000 3/10584 590.0 0.527 0.527 0.023 -0.000 -0.000 -0.000 3/10584 600.0 0.518 0.518 0.022 -0.000 -0.000 -0.000 3/10584 610.0 0.510 0.510 0.022 -0.000 -0.000 -0.000 3/10584 620.0 0.502 0.502 0.022 -0.000 -0.000 -0.000 3/10584 630.0 0.494 0.494 0.021 -0.000 -0.000 -0.000 3/10584 640.0 0.486 0.486 0.021 -0.000 -0.000 -0.000 3/10584 650.0 0.478 0.478 0.021 -0.000 -0.000 -0.000 3/10584 660.0 0.471 0.471 0.020 -0.000 -0.000 -0.000 3/10584 670.0 0.464 0.464 0.020 -0.000 -0.000 -0.000 3/10584 680.0 0.457 0.457 0.020 -0.000 -0.000 -0.000 3/10584 690.0 0.451 0.451 0.019 -0.000 -0.000 -0.000 3/10584 700.0 0.444 0.444 0.019 -0.000 -0.000 -0.000 3/10584 710.0 0.438 0.438 0.019 -0.000 -0.000 -0.000 3/10584 720.0 0.432 0.432 0.019 -0.000 -0.000 -0.000 3/10584 730.0 0.426 0.426 0.018 -0.000 -0.000 -0.000 3/10584 740.0 0.420 0.420 0.018 -0.000 -0.000 -0.000 3/10584 750.0 0.414 0.414 0.018 -0.000 -0.000 -0.000 3/10584 760.0 0.409 0.409 0.018 -0.000 -0.000 -0.000 3/10584 770.0 0.404 0.404 0.017 -0.000 -0.000 -0.000 3/10584 780.0 0.398 0.398 0.017 -0.000 -0.000 -0.000 3/10584 790.0 0.393 0.393 0.017 -0.000 -0.000 -0.000 3/10584 800.0 0.389 0.389 0.017 -0.000 -0.000 -0.000 3/10584 810.0 0.384 0.384 0.017 -0.000 -0.000 -0.000 3/10584 820.0 0.379 0.379 0.016 -0.000 -0.000 -0.000 3/10584 830.0 0.374 0.374 0.016 -0.000 -0.000 -0.000 3/10584 840.0 0.370 0.370 0.016 -0.000 -0.000 -0.000 3/10584 850.0 0.366 0.366 0.016 -0.000 -0.000 -0.000 3/10584 860.0 0.361 0.361 0.016 -0.000 -0.000 -0.000 3/10584 870.0 0.357 0.357 0.015 -0.000 -0.000 -0.000 3/10584 880.0 0.353 0.353 0.015 -0.000 -0.000 -0.000 3/10584 890.0 0.349 0.349 0.015 -0.000 -0.000 -0.000 3/10584 900.0 0.345 0.345 0.015 -0.000 -0.000 -0.000 3/10584 910.0 0.342 0.342 0.015 -0.000 -0.000 -0.000 3/10584 920.0 0.338 0.338 0.015 -0.000 -0.000 -0.000 3/10584 930.0 0.334 0.334 0.014 -0.000 -0.000 -0.000 3/10584 940.0 0.331 0.331 0.014 -0.000 -0.000 -0.000 3/10584 950.0 0.327 0.327 0.014 -0.000 -0.000 -0.000 3/10584 960.0 0.324 0.324 0.014 -0.000 -0.000 -0.000 3/10584 970.0 0.320 0.320 0.014 -0.000 -0.000 -0.000 3/10584 980.0 0.317 0.317 0.014 -0.000 -0.000 -0.000 3/10584 990.0 0.314 0.314 0.014 -0.000 -0.000 -0.000 3/10584 1000.0 0.311 0.311 0.013 -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 17:37:49]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|