----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 06:38:34]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [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.188586112463025 0.000000000000000 0.000000000000000 b -2.094293056231512 3.627421979331683 0.000000000000000 c 0.000000000000000 0.000000000000000 20.803385620000000 Atomic positions (fractional): *1 Ge 0.00000000000000 0.00000000000000 0.00000000000000 72.640 *2 Ge 0.33333333333333 0.66666666666667 0.33328191863494 72.640 3 Ge 0.66666666666667 0.33333333333333 0.66671808136506 72.640 *4 I 0.00000000000000 0.00000000000000 0.24981113415455 126.904 5 I 0.00000000000000 0.00000000000000 0.75018886584544 126.904 *6 I 0.33333333333333 0.66666666666667 0.08347332948300 126.904 7 I 0.66666666666667 0.33333333333333 0.91652667051700 126.904 *8 I 0.33333333333333 0.66666666666667 0.58324225270597 126.904 9 I 0.66666666666667 0.33333333333333 0.41675774729402 126.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.188586112463025 0.000000000000000 0.000000000000000 b -2.094293056231512 3.627421979331683 0.000000000000000 c 0.000000000000000 0.000000000000000 20.803385620000000 Atomic positions (fractional): *1 Ge 0.00000000000000 0.00000000000000 0.00000000000000 72.640 > 1 *2 Ge 0.33333333333333 0.66666666666667 0.33328191863494 72.640 > 2 3 Ge 0.66666666666667 0.33333333333333 0.66671808136506 72.640 > 3 *4 I 0.00000000000000 0.00000000000000 0.24981113415455 126.904 > 4 5 I 0.00000000000000 0.00000000000000 0.75018886584544 126.904 > 5 *6 I 0.33333333333333 0.66666666666667 0.08347332948300 126.904 > 6 7 I 0.66666666666667 0.33333333333333 0.91652667051700 126.904 > 7 *8 I 0.33333333333333 0.66666666666667 0.58324225270597 126.904 > 8 9 I 0.66666666666667 0.33333333333333 0.41675774729402 126.904 > 9 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.754344449852098 0.000000000000000 0.000000000000000 b -8.377172224926047 14.509687917326731 0.000000000000000 c 0.000000000000000 0.000000000000000 20.803385620000000 Atomic positions (fractional): *1 Ge 0.00000000000000 0.00000000000000 0.00000000000000 72.640 > 1 2 Ge 0.25000000000000 0.00000000000000 0.00000000000000 72.640 > 1 3 Ge 0.50000000000000 0.00000000000000 0.00000000000000 72.640 > 1 4 Ge 0.75000000000000 0.00000000000000 0.00000000000000 72.640 > 1 5 Ge 0.00000000000000 0.25000000000000 0.00000000000000 72.640 > 1 6 Ge 0.25000000000000 0.25000000000000 0.00000000000000 72.640 > 1 7 Ge 0.50000000000000 0.25000000000000 0.00000000000000 72.640 > 1 8 Ge 0.75000000000000 0.25000000000000 0.00000000000000 72.640 > 1 9 Ge 0.00000000000000 0.50000000000000 0.00000000000000 72.640 > 1 10 Ge 0.25000000000000 0.50000000000000 0.00000000000000 72.640 > 1 11 Ge 0.50000000000000 0.50000000000000 0.00000000000000 72.640 > 1 12 Ge 0.75000000000000 0.50000000000000 0.00000000000000 72.640 > 1 13 Ge 0.00000000000000 0.75000000000000 0.00000000000000 72.640 > 1 14 Ge 0.25000000000000 0.75000000000000 0.00000000000000 72.640 > 1 15 Ge 0.50000000000000 0.75000000000000 0.00000000000000 72.640 > 1 16 Ge 0.75000000000000 0.75000000000000 0.00000000000000 72.640 > 1 *17 Ge 0.08333333333333 0.16666666666667 0.33328191863494 72.640 > 2 18 Ge 0.33333333333333 0.16666666666667 0.33328191863494 72.640 > 2 19 Ge 0.58333333333333 0.16666666666667 0.33328191863494 72.640 > 2 20 Ge 0.83333333333333 0.16666666666667 0.33328191863494 72.640 > 2 21 Ge 0.08333333333333 0.41666666666667 0.33328191863494 72.640 > 2 22 Ge 0.33333333333333 0.41666666666667 0.33328191863494 72.640 > 2 23 Ge 0.58333333333333 0.41666666666667 0.33328191863494 72.640 > 2 24 Ge 0.83333333333333 0.41666666666667 0.33328191863494 72.640 > 2 25 Ge 0.08333333333333 0.66666666666667 0.33328191863494 72.640 > 2 26 Ge 0.33333333333333 0.66666666666667 0.33328191863494 72.640 > 2 27 Ge 0.58333333333333 0.66666666666667 0.33328191863494 72.640 > 2 28 Ge 0.83333333333333 0.66666666666667 0.33328191863494 72.640 > 2 29 Ge 0.08333333333333 0.91666666666667 0.33328191863494 72.640 > 2 30 Ge 0.33333333333333 0.91666666666667 0.33328191863494 72.640 > 2 31 Ge 0.58333333333333 0.91666666666667 0.33328191863494 72.640 > 2 32 Ge 0.83333333333333 0.91666666666667 0.33328191863494 72.640 > 2 33 Ge 0.16666666666667 0.08333333333333 0.66671808136506 72.640 > 3 34 Ge 0.41666666666667 0.08333333333333 0.66671808136506 72.640 > 3 35 Ge 0.66666666666667 0.08333333333333 0.66671808136506 72.640 > 3 36 Ge 0.91666666666667 0.08333333333333 0.66671808136506 72.640 > 3 37 Ge 0.16666666666667 0.33333333333333 0.66671808136506 72.640 > 3 38 Ge 0.41666666666667 0.33333333333333 0.66671808136506 72.640 > 3 39 Ge 0.66666666666667 0.33333333333333 0.66671808136506 72.640 > 3 40 Ge 0.91666666666667 0.33333333333333 0.66671808136506 72.640 > 3 41 Ge 0.16666666666667 0.58333333333333 0.66671808136506 72.640 > 3 42 Ge 0.41666666666667 0.58333333333333 0.66671808136506 72.640 > 3 43 Ge 0.66666666666667 0.58333333333333 0.66671808136506 72.640 > 3 44 Ge 0.91666666666667 0.58333333333333 0.66671808136506 72.640 > 3 45 Ge 0.16666666666667 0.83333333333333 0.66671808136506 72.640 > 3 46 Ge 0.41666666666667 0.83333333333333 0.66671808136506 72.640 > 3 47 Ge 0.66666666666667 0.83333333333333 0.66671808136506 72.640 > 3 48 Ge 0.91666666666667 0.83333333333333 0.66671808136506 72.640 > 3 *49 I 0.00000000000000 0.00000000000000 0.24981113415455 126.904 > 4 50 I 0.25000000000000 0.00000000000000 0.24981113415455 126.904 > 4 51 I 0.50000000000000 0.00000000000000 0.24981113415455 126.904 > 4 52 I 0.75000000000000 0.00000000000000 0.24981113415455 126.904 > 4 53 I 0.00000000000000 0.25000000000000 0.24981113415455 126.904 > 4 54 I 0.25000000000000 0.25000000000000 0.24981113415455 126.904 > 4 55 I 0.50000000000000 0.25000000000000 0.24981113415455 126.904 > 4 56 I 0.75000000000000 0.25000000000000 0.24981113415455 126.904 > 4 57 I 0.00000000000000 0.50000000000000 0.24981113415455 126.904 > 4 58 I 0.25000000000000 0.50000000000000 0.24981113415455 126.904 > 4 59 I 0.50000000000000 0.50000000000000 0.24981113415455 126.904 > 4 60 I 0.75000000000000 0.50000000000000 0.24981113415455 126.904 > 4 61 I 0.00000000000000 0.75000000000000 0.24981113415455 126.904 > 4 62 I 0.25000000000000 0.75000000000000 0.24981113415455 126.904 > 4 63 I 0.50000000000000 0.75000000000000 0.24981113415455 126.904 > 4 64 I 0.75000000000000 0.75000000000000 0.24981113415455 126.904 > 4 65 I 0.00000000000000 0.00000000000000 0.75018886584544 126.904 > 5 66 I 0.25000000000000 0.00000000000000 0.75018886584544 126.904 > 5 67 I 0.50000000000000 0.00000000000000 0.75018886584544 126.904 > 5 68 I 0.75000000000000 0.00000000000000 0.75018886584544 126.904 > 5 69 I 0.00000000000000 0.25000000000000 0.75018886584544 126.904 > 5 70 I 0.25000000000000 0.25000000000000 0.75018886584544 126.904 > 5 71 I 0.50000000000000 0.25000000000000 0.75018886584544 126.904 > 5 72 I 0.75000000000000 0.25000000000000 0.75018886584544 126.904 > 5 73 I 0.00000000000000 0.50000000000000 0.75018886584544 126.904 > 5 74 I 0.25000000000000 0.50000000000000 0.75018886584544 126.904 > 5 75 I 0.50000000000000 0.50000000000000 0.75018886584544 126.904 > 5 76 I 0.75000000000000 0.50000000000000 0.75018886584544 126.904 > 5 77 I 0.00000000000000 0.75000000000000 0.75018886584544 126.904 > 5 78 I 0.25000000000000 0.75000000000000 0.75018886584544 126.904 > 5 79 I 0.50000000000000 0.75000000000000 0.75018886584544 126.904 > 5 80 I 0.75000000000000 0.75000000000000 0.75018886584544 126.904 > 5 *81 I 0.08333333333333 0.16666666666667 0.08347332948300 126.904 > 6 82 I 0.33333333333333 0.16666666666667 0.08347332948300 126.904 > 6 83 I 0.58333333333333 0.16666666666667 0.08347332948300 126.904 > 6 84 I 0.83333333333333 0.16666666666667 0.08347332948300 126.904 > 6 85 I 0.08333333333333 0.41666666666667 0.08347332948300 126.904 > 6 86 I 0.33333333333333 0.41666666666667 0.08347332948300 126.904 > 6 87 I 0.58333333333333 0.41666666666667 0.08347332948300 126.904 > 6 88 I 0.83333333333333 0.41666666666667 0.08347332948300 126.904 > 6 89 I 0.08333333333333 0.66666666666667 0.08347332948300 126.904 > 6 90 I 0.33333333333333 0.66666666666667 0.08347332948300 126.904 > 6 91 I 0.58333333333333 0.66666666666667 0.08347332948300 126.904 > 6 92 I 0.83333333333333 0.66666666666667 0.08347332948300 126.904 > 6 93 I 0.08333333333333 0.91666666666667 0.08347332948300 126.904 > 6 94 I 0.33333333333333 0.91666666666667 0.08347332948300 126.904 > 6 95 I 0.58333333333333 0.91666666666667 0.08347332948300 126.904 > 6 96 I 0.83333333333333 0.91666666666667 0.08347332948300 126.904 > 6 97 I 0.16666666666667 0.08333333333333 0.91652667051700 126.904 > 7 98 I 0.41666666666667 0.08333333333333 0.91652667051700 126.904 > 7 99 I 0.66666666666667 0.08333333333333 0.91652667051700 126.904 > 7 100 I 0.91666666666667 0.08333333333333 0.91652667051700 126.904 > 7 101 I 0.16666666666667 0.33333333333333 0.91652667051700 126.904 > 7 102 I 0.41666666666667 0.33333333333333 0.91652667051700 126.904 > 7 103 I 0.66666666666667 0.33333333333333 0.91652667051700 126.904 > 7 104 I 0.91666666666667 0.33333333333333 0.91652667051700 126.904 > 7 105 I 0.16666666666667 0.58333333333333 0.91652667051700 126.904 > 7 106 I 0.41666666666667 0.58333333333333 0.91652667051700 126.904 > 7 107 I 0.66666666666667 0.58333333333333 0.91652667051700 126.904 > 7 108 I 0.91666666666667 0.58333333333333 0.91652667051700 126.904 > 7 109 I 0.16666666666667 0.83333333333333 0.91652667051700 126.904 > 7 110 I 0.41666666666667 0.83333333333333 0.91652667051700 126.904 > 7 111 I 0.66666666666667 0.83333333333333 0.91652667051700 126.904 > 7 112 I 0.91666666666667 0.83333333333333 0.91652667051700 126.904 > 7 *113 I 0.08333333333333 0.16666666666667 0.58324225270597 126.904 > 8 114 I 0.33333333333333 0.16666666666667 0.58324225270597 126.904 > 8 115 I 0.58333333333333 0.16666666666667 0.58324225270597 126.904 > 8 116 I 0.83333333333333 0.16666666666667 0.58324225270597 126.904 > 8 117 I 0.08333333333333 0.41666666666667 0.58324225270597 126.904 > 8 118 I 0.33333333333333 0.41666666666667 0.58324225270597 126.904 > 8 119 I 0.58333333333333 0.41666666666667 0.58324225270597 126.904 > 8 120 I 0.83333333333333 0.41666666666667 0.58324225270597 126.904 > 8 121 I 0.08333333333333 0.66666666666667 0.58324225270597 126.904 > 8 122 I 0.33333333333333 0.66666666666667 0.58324225270597 126.904 > 8 123 I 0.58333333333333 0.66666666666667 0.58324225270597 126.904 > 8 124 I 0.83333333333333 0.66666666666667 0.58324225270597 126.904 > 8 125 I 0.08333333333333 0.91666666666667 0.58324225270597 126.904 > 8 126 I 0.33333333333333 0.91666666666667 0.58324225270597 126.904 > 8 127 I 0.58333333333333 0.91666666666667 0.58324225270597 126.904 > 8 128 I 0.83333333333333 0.91666666666667 0.58324225270597 126.904 > 8 129 I 0.16666666666667 0.08333333333333 0.41675774729402 126.904 > 9 130 I 0.41666666666667 0.08333333333333 0.41675774729402 126.904 > 9 131 I 0.66666666666667 0.08333333333333 0.41675774729402 126.904 > 9 132 I 0.91666666666667 0.08333333333333 0.41675774729402 126.904 > 9 133 I 0.16666666666667 0.33333333333333 0.41675774729402 126.904 > 9 134 I 0.41666666666667 0.33333333333333 0.41675774729402 126.904 > 9 135 I 0.66666666666667 0.33333333333333 0.41675774729402 126.904 > 9 136 I 0.91666666666667 0.33333333333333 0.41675774729402 126.904 > 9 137 I 0.16666666666667 0.58333333333333 0.41675774729402 126.904 > 9 138 I 0.41666666666667 0.58333333333333 0.41675774729402 126.904 > 9 139 I 0.66666666666667 0.58333333333333 0.41675774729402 126.904 > 9 140 I 0.91666666666667 0.58333333333333 0.41675774729402 126.904 > 9 141 I 0.16666666666667 0.83333333333333 0.41675774729402 126.904 > 9 142 I 0.41666666666667 0.83333333333333 0.41675774729402 126.904 > 9 143 I 0.66666666666667 0.83333333333333 0.41675774729402 126.904 > 9 144 I 0.91666666666667 0.83333333333333 0.41675774729402 126.904 > 9 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.1809291 0.0000000 0.0000000 0.0000000 9.1809291 0.0000000 0.0000000 0.0000000 5.3931443 -------------------------- Born effective charges -------------------------- 1 Ge 4.5750768 0.0000000 0.0000000 0.0000000 4.5750768 0.0000000 0.0000000 0.0000000 1.4957523 2 Ge 4.5751309 0.0000000 0.0000000 0.0000000 4.5751309 0.0000000 0.0000000 0.0000000 1.4963455 3 Ge 4.5751309 0.0000000 0.0000000 0.0000000 4.5751309 0.0000000 0.0000000 0.0000000 1.4963455 4 I -2.2876063 -0.0000000 0.0000000 0.0000000 -2.2876063 0.0000000 0.0000000 0.0000000 -0.7477991 5 I -2.2876063 -0.0000000 0.0000000 0.0000000 -2.2876063 0.0000000 0.0000000 0.0000000 -0.7477991 6 I -2.2875520 -0.0000000 0.0000000 0.0000000 -2.2875520 0.0000000 0.0000000 0.0000000 -0.7478594 7 I -2.2875520 -0.0000000 0.0000000 0.0000000 -2.2875520 0.0000000 0.0000000 0.0000000 -0.7478594 8 I -2.2875110 -0.0000000 0.0000000 0.0000000 -2.2875110 0.0000000 0.0000000 0.0000000 -0.7485632 9 I -2.2875110 -0.0000000 0.0000000 0.0000000 -2.2875110 0.0000000 0.0000000 0.0000000 -0.7485632 ---------------------------------------------------------------------------- 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: 484 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.016 Solver_block: 100 / 160 - Time: 0.595 Solver_block: 160 / 160 - Time: 0.409 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 1.025 --------------------------------- 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: 484 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 06:38:44]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 06:38:44]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.188586112463025 0.000000000000000 0.000000000000000 b -2.094293056231512 3.627421979331683 0.000000000000000 c 0.000000000000000 0.000000000000000 20.803385620000000 Atomic positions (fractional): 1 Ge 0.00000000000000 0.00000000000000 0.00000000000000 72.640 2 Ge 0.33333333333333 0.66666666666667 0.33328191863494 72.640 3 Ge 0.66666666666667 0.33333333333333 0.66671808136506 72.640 4 I 0.00000000000000 0.00000000000000 0.24981113415455 126.904 5 I 0.00000000000000 0.00000000000000 0.75018886584544 126.904 6 I 0.33333333333333 0.66666666666667 0.08347332948300 126.904 7 I 0.66666666666667 0.33333333333333 0.91652667051700 126.904 8 I 0.33333333333333 0.66666666666667 0.58324225270597 126.904 9 I 0.66666666666667 0.33333333333333 0.41675774729402 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.754344449852098 0.000000000000000 0.000000000000000 b -8.377172224926047 14.509687917326731 0.000000000000000 c 0.000000000000000 0.000000000000000 20.803385620000000 Atomic positions (fractional): 1 Ge 0.00000000000000 0.00000000000000 0.00000000000000 72.640 > 1 2 Ge 0.25000000000000 0.00000000000000 0.00000000000000 72.640 > 1 3 Ge 0.50000000000000 0.00000000000000 0.00000000000000 72.640 > 1 4 Ge 0.75000000000000 0.00000000000000 0.00000000000000 72.640 > 1 5 Ge 0.00000000000000 0.25000000000000 0.00000000000000 72.640 > 1 6 Ge 0.25000000000000 0.25000000000000 0.00000000000000 72.640 > 1 7 Ge 0.50000000000000 0.25000000000000 0.00000000000000 72.640 > 1 8 Ge 0.75000000000000 0.25000000000000 0.00000000000000 72.640 > 1 9 Ge 0.00000000000000 0.50000000000000 0.00000000000000 72.640 > 1 10 Ge 0.25000000000000 0.50000000000000 0.00000000000000 72.640 > 1 11 Ge 0.50000000000000 0.50000000000000 0.00000000000000 72.640 > 1 12 Ge 0.75000000000000 0.50000000000000 0.00000000000000 72.640 > 1 13 Ge 0.00000000000000 0.75000000000000 0.00000000000000 72.640 > 1 14 Ge 0.25000000000000 0.75000000000000 0.00000000000000 72.640 > 1 15 Ge 0.50000000000000 0.75000000000000 0.00000000000000 72.640 > 1 16 Ge 0.75000000000000 0.75000000000000 0.00000000000000 72.640 > 1 17 Ge 0.08333333333333 0.16666666666667 0.33328191863494 72.640 > 17 18 Ge 0.33333333333333 0.16666666666667 0.33328191863494 72.640 > 17 19 Ge 0.58333333333333 0.16666666666667 0.33328191863494 72.640 > 17 20 Ge 0.83333333333333 0.16666666666667 0.33328191863494 72.640 > 17 21 Ge 0.08333333333333 0.41666666666667 0.33328191863494 72.640 > 17 22 Ge 0.33333333333333 0.41666666666667 0.33328191863494 72.640 > 17 23 Ge 0.58333333333333 0.41666666666667 0.33328191863494 72.640 > 17 24 Ge 0.83333333333333 0.41666666666667 0.33328191863494 72.640 > 17 25 Ge 0.08333333333333 0.66666666666667 0.33328191863494 72.640 > 17 26 Ge 0.33333333333333 0.66666666666667 0.33328191863494 72.640 > 17 27 Ge 0.58333333333333 0.66666666666667 0.33328191863494 72.640 > 17 28 Ge 0.83333333333333 0.66666666666667 0.33328191863494 72.640 > 17 29 Ge 0.08333333333333 0.91666666666667 0.33328191863494 72.640 > 17 30 Ge 0.33333333333333 0.91666666666667 0.33328191863494 72.640 > 17 31 Ge 0.58333333333333 0.91666666666667 0.33328191863494 72.640 > 17 32 Ge 0.83333333333333 0.91666666666667 0.33328191863494 72.640 > 17 33 Ge 0.16666666666667 0.08333333333333 0.66671808136506 72.640 > 33 34 Ge 0.41666666666667 0.08333333333333 0.66671808136506 72.640 > 33 35 Ge 0.66666666666667 0.08333333333333 0.66671808136506 72.640 > 33 36 Ge 0.91666666666667 0.08333333333333 0.66671808136506 72.640 > 33 37 Ge 0.16666666666667 0.33333333333333 0.66671808136506 72.640 > 33 38 Ge 0.41666666666667 0.33333333333333 0.66671808136506 72.640 > 33 39 Ge 0.66666666666667 0.33333333333333 0.66671808136506 72.640 > 33 40 Ge 0.91666666666667 0.33333333333333 0.66671808136506 72.640 > 33 41 Ge 0.16666666666667 0.58333333333333 0.66671808136506 72.640 > 33 42 Ge 0.41666666666667 0.58333333333333 0.66671808136506 72.640 > 33 43 Ge 0.66666666666667 0.58333333333333 0.66671808136506 72.640 > 33 44 Ge 0.91666666666667 0.58333333333333 0.66671808136506 72.640 > 33 45 Ge 0.16666666666667 0.83333333333333 0.66671808136506 72.640 > 33 46 Ge 0.41666666666667 0.83333333333333 0.66671808136506 72.640 > 33 47 Ge 0.66666666666667 0.83333333333333 0.66671808136506 72.640 > 33 48 Ge 0.91666666666667 0.83333333333333 0.66671808136506 72.640 > 33 49 I 0.00000000000000 0.00000000000000 0.24981113415455 126.904 > 49 50 I 0.25000000000000 0.00000000000000 0.24981113415455 126.904 > 49 51 I 0.50000000000000 0.00000000000000 0.24981113415455 126.904 > 49 52 I 0.75000000000000 0.00000000000000 0.24981113415455 126.904 > 49 53 I 0.00000000000000 0.25000000000000 0.24981113415455 126.904 > 49 54 I 0.25000000000000 0.25000000000000 0.24981113415455 126.904 > 49 55 I 0.50000000000000 0.25000000000000 0.24981113415455 126.904 > 49 56 I 0.75000000000000 0.25000000000000 0.24981113415455 126.904 > 49 57 I 0.00000000000000 0.50000000000000 0.24981113415455 126.904 > 49 58 I 0.25000000000000 0.50000000000000 0.24981113415455 126.904 > 49 59 I 0.50000000000000 0.50000000000000 0.24981113415455 126.904 > 49 60 I 0.75000000000000 0.50000000000000 0.24981113415455 126.904 > 49 61 I 0.00000000000000 0.75000000000000 0.24981113415455 126.904 > 49 62 I 0.25000000000000 0.75000000000000 0.24981113415455 126.904 > 49 63 I 0.50000000000000 0.75000000000000 0.24981113415455 126.904 > 49 64 I 0.75000000000000 0.75000000000000 0.24981113415455 126.904 > 49 65 I 0.00000000000000 0.00000000000000 0.75018886584544 126.904 > 65 66 I 0.25000000000000 0.00000000000000 0.75018886584544 126.904 > 65 67 I 0.50000000000000 0.00000000000000 0.75018886584544 126.904 > 65 68 I 0.75000000000000 0.00000000000000 0.75018886584544 126.904 > 65 69 I 0.00000000000000 0.25000000000000 0.75018886584544 126.904 > 65 70 I 0.25000000000000 0.25000000000000 0.75018886584544 126.904 > 65 71 I 0.50000000000000 0.25000000000000 0.75018886584544 126.904 > 65 72 I 0.75000000000000 0.25000000000000 0.75018886584544 126.904 > 65 73 I 0.00000000000000 0.50000000000000 0.75018886584544 126.904 > 65 74 I 0.25000000000000 0.50000000000000 0.75018886584544 126.904 > 65 75 I 0.50000000000000 0.50000000000000 0.75018886584544 126.904 > 65 76 I 0.75000000000000 0.50000000000000 0.75018886584544 126.904 > 65 77 I 0.00000000000000 0.75000000000000 0.75018886584544 126.904 > 65 78 I 0.25000000000000 0.75000000000000 0.75018886584544 126.904 > 65 79 I 0.50000000000000 0.75000000000000 0.75018886584544 126.904 > 65 80 I 0.75000000000000 0.75000000000000 0.75018886584544 126.904 > 65 81 I 0.08333333333333 0.16666666666667 0.08347332948300 126.904 > 81 82 I 0.33333333333333 0.16666666666667 0.08347332948300 126.904 > 81 83 I 0.58333333333333 0.16666666666667 0.08347332948300 126.904 > 81 84 I 0.83333333333333 0.16666666666667 0.08347332948300 126.904 > 81 85 I 0.08333333333333 0.41666666666667 0.08347332948300 126.904 > 81 86 I 0.33333333333333 0.41666666666667 0.08347332948300 126.904 > 81 87 I 0.58333333333333 0.41666666666667 0.08347332948300 126.904 > 81 88 I 0.83333333333333 0.41666666666667 0.08347332948300 126.904 > 81 89 I 0.08333333333333 0.66666666666667 0.08347332948300 126.904 > 81 90 I 0.33333333333333 0.66666666666667 0.08347332948300 126.904 > 81 91 I 0.58333333333333 0.66666666666667 0.08347332948300 126.904 > 81 92 I 0.83333333333333 0.66666666666667 0.08347332948300 126.904 > 81 93 I 0.08333333333333 0.91666666666667 0.08347332948300 126.904 > 81 94 I 0.33333333333333 0.91666666666667 0.08347332948300 126.904 > 81 95 I 0.58333333333333 0.91666666666667 0.08347332948300 126.904 > 81 96 I 0.83333333333333 0.91666666666667 0.08347332948300 126.904 > 81 97 I 0.16666666666667 0.08333333333333 0.91652667051700 126.904 > 97 98 I 0.41666666666667 0.08333333333333 0.91652667051700 126.904 > 97 99 I 0.66666666666667 0.08333333333333 0.91652667051700 126.904 > 97 100 I 0.91666666666667 0.08333333333333 0.91652667051700 126.904 > 97 101 I 0.16666666666667 0.33333333333333 0.91652667051700 126.904 > 97 102 I 0.41666666666667 0.33333333333333 0.91652667051700 126.904 > 97 103 I 0.66666666666667 0.33333333333333 0.91652667051700 126.904 > 97 104 I 0.91666666666667 0.33333333333333 0.91652667051700 126.904 > 97 105 I 0.16666666666667 0.58333333333333 0.91652667051700 126.904 > 97 106 I 0.41666666666667 0.58333333333333 0.91652667051700 126.904 > 97 107 I 0.66666666666667 0.58333333333333 0.91652667051700 126.904 > 97 108 I 0.91666666666667 0.58333333333333 0.91652667051700 126.904 > 97 109 I 0.16666666666667 0.83333333333333 0.91652667051700 126.904 > 97 110 I 0.41666666666667 0.83333333333333 0.91652667051700 126.904 > 97 111 I 0.66666666666667 0.83333333333333 0.91652667051700 126.904 > 97 112 I 0.91666666666667 0.83333333333333 0.91652667051700 126.904 > 97 113 I 0.08333333333333 0.16666666666667 0.58324225270597 126.904 > 113 114 I 0.33333333333333 0.16666666666667 0.58324225270597 126.904 > 113 115 I 0.58333333333333 0.16666666666667 0.58324225270597 126.904 > 113 116 I 0.83333333333333 0.16666666666667 0.58324225270597 126.904 > 113 117 I 0.08333333333333 0.41666666666667 0.58324225270597 126.904 > 113 118 I 0.33333333333333 0.41666666666667 0.58324225270597 126.904 > 113 119 I 0.58333333333333 0.41666666666667 0.58324225270597 126.904 > 113 120 I 0.83333333333333 0.41666666666667 0.58324225270597 126.904 > 113 121 I 0.08333333333333 0.66666666666667 0.58324225270597 126.904 > 113 122 I 0.33333333333333 0.66666666666667 0.58324225270597 126.904 > 113 123 I 0.58333333333333 0.66666666666667 0.58324225270597 126.904 > 113 124 I 0.83333333333333 0.66666666666667 0.58324225270597 126.904 > 113 125 I 0.08333333333333 0.91666666666667 0.58324225270597 126.904 > 113 126 I 0.33333333333333 0.91666666666667 0.58324225270597 126.904 > 113 127 I 0.58333333333333 0.91666666666667 0.58324225270597 126.904 > 113 128 I 0.83333333333333 0.91666666666667 0.58324225270597 126.904 > 113 129 I 0.16666666666667 0.08333333333333 0.41675774729402 126.904 > 129 130 I 0.41666666666667 0.08333333333333 0.41675774729402 126.904 > 129 131 I 0.66666666666667 0.08333333333333 0.41675774729402 126.904 > 129 132 I 0.91666666666667 0.08333333333333 0.41675774729402 126.904 > 129 133 I 0.16666666666667 0.33333333333333 0.41675774729402 126.904 > 129 134 I 0.41666666666667 0.33333333333333 0.41675774729402 126.904 > 129 135 I 0.66666666666667 0.33333333333333 0.41675774729402 126.904 > 129 136 I 0.91666666666667 0.33333333333333 0.41675774729402 126.904 > 129 137 I 0.16666666666667 0.58333333333333 0.41675774729402 126.904 > 129 138 I 0.41666666666667 0.58333333333333 0.41675774729402 126.904 > 129 139 I 0.66666666666667 0.58333333333333 0.41675774729402 126.904 > 129 140 I 0.91666666666667 0.58333333333333 0.41675774729402 126.904 > 129 141 I 0.16666666666667 0.83333333333333 0.41675774729402 126.904 > 129 142 I 0.41666666666667 0.83333333333333 0.41675774729402 126.904 > 129 143 I 0.66666666666667 0.83333333333333 0.41675774729402 126.904 > 129 144 I 0.91666666666667 0.83333333333333 0.41675774729402 126.904 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.1809291 0.0000000 0.0000000 0.0000000 9.1809291 0.0000000 0.0000000 0.0000000 5.3931443 -------------------------- Born effective charges -------------------------- 1 Ge 4.5750768 0.0000000 0.0000000 0.0000000 4.5750768 0.0000000 0.0000000 0.0000000 1.4957523 2 Ge 4.5751309 0.0000000 0.0000000 0.0000000 4.5751309 0.0000000 0.0000000 0.0000000 1.4963455 3 Ge 4.5751309 0.0000000 0.0000000 0.0000000 4.5751309 0.0000000 0.0000000 0.0000000 1.4963455 4 I -2.2876063 -0.0000000 0.0000000 0.0000000 -2.2876063 0.0000000 0.0000000 0.0000000 -0.7477991 5 I -2.2876063 -0.0000000 0.0000000 0.0000000 -2.2876063 0.0000000 0.0000000 0.0000000 -0.7477991 6 I -2.2875520 -0.0000000 0.0000000 0.0000000 -2.2875520 0.0000000 0.0000000 0.0000000 -0.7478594 7 I -2.2875520 -0.0000000 0.0000000 0.0000000 -2.2875520 0.0000000 0.0000000 0.0000000 -0.7478594 8 I -2.2875110 -0.0000000 0.0000000 0.0000000 -2.2875110 0.0000000 0.0000000 0.0000000 -0.7485632 9 I -2.2875110 -0.0000000 0.0000000 0.0000000 -2.2875110 0.0000000 0.0000000 0.0000000 -0.7485632 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 17, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 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[ 81, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 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.00000007 (yyy) 0.00000007 (yyy) 0.00000007 (yyy) 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 06:38:58]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 06:38:59]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.188586112463025 0.000000000000000 0.000000000000000 b -2.094293056231512 3.627421979331683 0.000000000000000 c 0.000000000000000 0.000000000000000 20.803385620000000 Atomic positions (fractional): 1 Ge 0.00000000000000 0.00000000000000 0.00000000000000 72.640 2 Ge 0.33333333333333 0.66666666666667 0.33328191863494 72.640 3 Ge 0.66666666666667 0.33333333333333 0.66671808136506 72.640 4 I 0.00000000000000 0.00000000000000 0.24981113415455 126.904 5 I 0.00000000000000 0.00000000000000 0.75018886584544 126.904 6 I 0.33333333333333 0.66666666666667 0.08347332948300 126.904 7 I 0.66666666666667 0.33333333333333 0.91652667051700 126.904 8 I 0.33333333333333 0.66666666666667 0.58324225270597 126.904 9 I 0.66666666666667 0.33333333333333 0.41675774729402 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.754344449852098 0.000000000000000 0.000000000000000 b -8.377172224926047 14.509687917326731 0.000000000000000 c 0.000000000000000 0.000000000000000 20.803385620000000 Atomic positions (fractional): 1 Ge 0.00000000000000 0.00000000000000 0.00000000000000 72.640 > 1 2 Ge 0.25000000000000 0.00000000000000 0.00000000000000 72.640 > 1 3 Ge 0.50000000000000 0.00000000000000 0.00000000000000 72.640 > 1 4 Ge 0.75000000000000 0.00000000000000 0.00000000000000 72.640 > 1 5 Ge 0.00000000000000 0.25000000000000 0.00000000000000 72.640 > 1 6 Ge 0.25000000000000 0.25000000000000 0.00000000000000 72.640 > 1 7 Ge 0.50000000000000 0.25000000000000 0.00000000000000 72.640 > 1 8 Ge 0.75000000000000 0.25000000000000 0.00000000000000 72.640 > 1 9 Ge 0.00000000000000 0.50000000000000 0.00000000000000 72.640 > 1 10 Ge 0.25000000000000 0.50000000000000 0.00000000000000 72.640 > 1 11 Ge 0.50000000000000 0.50000000000000 0.00000000000000 72.640 > 1 12 Ge 0.75000000000000 0.50000000000000 0.00000000000000 72.640 > 1 13 Ge 0.00000000000000 0.75000000000000 0.00000000000000 72.640 > 1 14 Ge 0.25000000000000 0.75000000000000 0.00000000000000 72.640 > 1 15 Ge 0.50000000000000 0.75000000000000 0.00000000000000 72.640 > 1 16 Ge 0.75000000000000 0.75000000000000 0.00000000000000 72.640 > 1 17 Ge 0.08333333333333 0.16666666666667 0.33328191863494 72.640 > 17 18 Ge 0.33333333333333 0.16666666666667 0.33328191863494 72.640 > 17 19 Ge 0.58333333333333 0.16666666666667 0.33328191863494 72.640 > 17 20 Ge 0.83333333333333 0.16666666666667 0.33328191863494 72.640 > 17 21 Ge 0.08333333333333 0.41666666666667 0.33328191863494 72.640 > 17 22 Ge 0.33333333333333 0.41666666666667 0.33328191863494 72.640 > 17 23 Ge 0.58333333333333 0.41666666666667 0.33328191863494 72.640 > 17 24 Ge 0.83333333333333 0.41666666666667 0.33328191863494 72.640 > 17 25 Ge 0.08333333333333 0.66666666666667 0.33328191863494 72.640 > 17 26 Ge 0.33333333333333 0.66666666666667 0.33328191863494 72.640 > 17 27 Ge 0.58333333333333 0.66666666666667 0.33328191863494 72.640 > 17 28 Ge 0.83333333333333 0.66666666666667 0.33328191863494 72.640 > 17 29 Ge 0.08333333333333 0.91666666666667 0.33328191863494 72.640 > 17 30 Ge 0.33333333333333 0.91666666666667 0.33328191863494 72.640 > 17 31 Ge 0.58333333333333 0.91666666666667 0.33328191863494 72.640 > 17 32 Ge 0.83333333333333 0.91666666666667 0.33328191863494 72.640 > 17 33 Ge 0.16666666666667 0.08333333333333 0.66671808136506 72.640 > 33 34 Ge 0.41666666666667 0.08333333333333 0.66671808136506 72.640 > 33 35 Ge 0.66666666666667 0.08333333333333 0.66671808136506 72.640 > 33 36 Ge 0.91666666666667 0.08333333333333 0.66671808136506 72.640 > 33 37 Ge 0.16666666666667 0.33333333333333 0.66671808136506 72.640 > 33 38 Ge 0.41666666666667 0.33333333333333 0.66671808136506 72.640 > 33 39 Ge 0.66666666666667 0.33333333333333 0.66671808136506 72.640 > 33 40 Ge 0.91666666666667 0.33333333333333 0.66671808136506 72.640 > 33 41 Ge 0.16666666666667 0.58333333333333 0.66671808136506 72.640 > 33 42 Ge 0.41666666666667 0.58333333333333 0.66671808136506 72.640 > 33 43 Ge 0.66666666666667 0.58333333333333 0.66671808136506 72.640 > 33 44 Ge 0.91666666666667 0.58333333333333 0.66671808136506 72.640 > 33 45 Ge 0.16666666666667 0.83333333333333 0.66671808136506 72.640 > 33 46 Ge 0.41666666666667 0.83333333333333 0.66671808136506 72.640 > 33 47 Ge 0.66666666666667 0.83333333333333 0.66671808136506 72.640 > 33 48 Ge 0.91666666666667 0.83333333333333 0.66671808136506 72.640 > 33 49 I 0.00000000000000 0.00000000000000 0.24981113415455 126.904 > 49 50 I 0.25000000000000 0.00000000000000 0.24981113415455 126.904 > 49 51 I 0.50000000000000 0.00000000000000 0.24981113415455 126.904 > 49 52 I 0.75000000000000 0.00000000000000 0.24981113415455 126.904 > 49 53 I 0.00000000000000 0.25000000000000 0.24981113415455 126.904 > 49 54 I 0.25000000000000 0.25000000000000 0.24981113415455 126.904 > 49 55 I 0.50000000000000 0.25000000000000 0.24981113415455 126.904 > 49 56 I 0.75000000000000 0.25000000000000 0.24981113415455 126.904 > 49 57 I 0.00000000000000 0.50000000000000 0.24981113415455 126.904 > 49 58 I 0.25000000000000 0.50000000000000 0.24981113415455 126.904 > 49 59 I 0.50000000000000 0.50000000000000 0.24981113415455 126.904 > 49 60 I 0.75000000000000 0.50000000000000 0.24981113415455 126.904 > 49 61 I 0.00000000000000 0.75000000000000 0.24981113415455 126.904 > 49 62 I 0.25000000000000 0.75000000000000 0.24981113415455 126.904 > 49 63 I 0.50000000000000 0.75000000000000 0.24981113415455 126.904 > 49 64 I 0.75000000000000 0.75000000000000 0.24981113415455 126.904 > 49 65 I 0.00000000000000 0.00000000000000 0.75018886584544 126.904 > 65 66 I 0.25000000000000 0.00000000000000 0.75018886584544 126.904 > 65 67 I 0.50000000000000 0.00000000000000 0.75018886584544 126.904 > 65 68 I 0.75000000000000 0.00000000000000 0.75018886584544 126.904 > 65 69 I 0.00000000000000 0.25000000000000 0.75018886584544 126.904 > 65 70 I 0.25000000000000 0.25000000000000 0.75018886584544 126.904 > 65 71 I 0.50000000000000 0.25000000000000 0.75018886584544 126.904 > 65 72 I 0.75000000000000 0.25000000000000 0.75018886584544 126.904 > 65 73 I 0.00000000000000 0.50000000000000 0.75018886584544 126.904 > 65 74 I 0.25000000000000 0.50000000000000 0.75018886584544 126.904 > 65 75 I 0.50000000000000 0.50000000000000 0.75018886584544 126.904 > 65 76 I 0.75000000000000 0.50000000000000 0.75018886584544 126.904 > 65 77 I 0.00000000000000 0.75000000000000 0.75018886584544 126.904 > 65 78 I 0.25000000000000 0.75000000000000 0.75018886584544 126.904 > 65 79 I 0.50000000000000 0.75000000000000 0.75018886584544 126.904 > 65 80 I 0.75000000000000 0.75000000000000 0.75018886584544 126.904 > 65 81 I 0.08333333333333 0.16666666666667 0.08347332948300 126.904 > 81 82 I 0.33333333333333 0.16666666666667 0.08347332948300 126.904 > 81 83 I 0.58333333333333 0.16666666666667 0.08347332948300 126.904 > 81 84 I 0.83333333333333 0.16666666666667 0.08347332948300 126.904 > 81 85 I 0.08333333333333 0.41666666666667 0.08347332948300 126.904 > 81 86 I 0.33333333333333 0.41666666666667 0.08347332948300 126.904 > 81 87 I 0.58333333333333 0.41666666666667 0.08347332948300 126.904 > 81 88 I 0.83333333333333 0.41666666666667 0.08347332948300 126.904 > 81 89 I 0.08333333333333 0.66666666666667 0.08347332948300 126.904 > 81 90 I 0.33333333333333 0.66666666666667 0.08347332948300 126.904 > 81 91 I 0.58333333333333 0.66666666666667 0.08347332948300 126.904 > 81 92 I 0.83333333333333 0.66666666666667 0.08347332948300 126.904 > 81 93 I 0.08333333333333 0.91666666666667 0.08347332948300 126.904 > 81 94 I 0.33333333333333 0.91666666666667 0.08347332948300 126.904 > 81 95 I 0.58333333333333 0.91666666666667 0.08347332948300 126.904 > 81 96 I 0.83333333333333 0.91666666666667 0.08347332948300 126.904 > 81 97 I 0.16666666666667 0.08333333333333 0.91652667051700 126.904 > 97 98 I 0.41666666666667 0.08333333333333 0.91652667051700 126.904 > 97 99 I 0.66666666666667 0.08333333333333 0.91652667051700 126.904 > 97 100 I 0.91666666666667 0.08333333333333 0.91652667051700 126.904 > 97 101 I 0.16666666666667 0.33333333333333 0.91652667051700 126.904 > 97 102 I 0.41666666666667 0.33333333333333 0.91652667051700 126.904 > 97 103 I 0.66666666666667 0.33333333333333 0.91652667051700 126.904 > 97 104 I 0.91666666666667 0.33333333333333 0.91652667051700 126.904 > 97 105 I 0.16666666666667 0.58333333333333 0.91652667051700 126.904 > 97 106 I 0.41666666666667 0.58333333333333 0.91652667051700 126.904 > 97 107 I 0.66666666666667 0.58333333333333 0.91652667051700 126.904 > 97 108 I 0.91666666666667 0.58333333333333 0.91652667051700 126.904 > 97 109 I 0.16666666666667 0.83333333333333 0.91652667051700 126.904 > 97 110 I 0.41666666666667 0.83333333333333 0.91652667051700 126.904 > 97 111 I 0.66666666666667 0.83333333333333 0.91652667051700 126.904 > 97 112 I 0.91666666666667 0.83333333333333 0.91652667051700 126.904 > 97 113 I 0.08333333333333 0.16666666666667 0.58324225270597 126.904 > 113 114 I 0.33333333333333 0.16666666666667 0.58324225270597 126.904 > 113 115 I 0.58333333333333 0.16666666666667 0.58324225270597 126.904 > 113 116 I 0.83333333333333 0.16666666666667 0.58324225270597 126.904 > 113 117 I 0.08333333333333 0.41666666666667 0.58324225270597 126.904 > 113 118 I 0.33333333333333 0.41666666666667 0.58324225270597 126.904 > 113 119 I 0.58333333333333 0.41666666666667 0.58324225270597 126.904 > 113 120 I 0.83333333333333 0.41666666666667 0.58324225270597 126.904 > 113 121 I 0.08333333333333 0.66666666666667 0.58324225270597 126.904 > 113 122 I 0.33333333333333 0.66666666666667 0.58324225270597 126.904 > 113 123 I 0.58333333333333 0.66666666666667 0.58324225270597 126.904 > 113 124 I 0.83333333333333 0.66666666666667 0.58324225270597 126.904 > 113 125 I 0.08333333333333 0.91666666666667 0.58324225270597 126.904 > 113 126 I 0.33333333333333 0.91666666666667 0.58324225270597 126.904 > 113 127 I 0.58333333333333 0.91666666666667 0.58324225270597 126.904 > 113 128 I 0.83333333333333 0.91666666666667 0.58324225270597 126.904 > 113 129 I 0.16666666666667 0.08333333333333 0.41675774729402 126.904 > 129 130 I 0.41666666666667 0.08333333333333 0.41675774729402 126.904 > 129 131 I 0.66666666666667 0.08333333333333 0.41675774729402 126.904 > 129 132 I 0.91666666666667 0.08333333333333 0.41675774729402 126.904 > 129 133 I 0.16666666666667 0.33333333333333 0.41675774729402 126.904 > 129 134 I 0.41666666666667 0.33333333333333 0.41675774729402 126.904 > 129 135 I 0.66666666666667 0.33333333333333 0.41675774729402 126.904 > 129 136 I 0.91666666666667 0.33333333333333 0.41675774729402 126.904 > 129 137 I 0.16666666666667 0.58333333333333 0.41675774729402 126.904 > 129 138 I 0.41666666666667 0.58333333333333 0.41675774729402 126.904 > 129 139 I 0.66666666666667 0.58333333333333 0.41675774729402 126.904 > 129 140 I 0.91666666666667 0.58333333333333 0.41675774729402 126.904 > 129 141 I 0.16666666666667 0.83333333333333 0.41675774729402 126.904 > 129 142 I 0.41666666666667 0.83333333333333 0.41675774729402 126.904 > 129 143 I 0.66666666666667 0.83333333333333 0.41675774729402 126.904 > 129 144 I 0.91666666666667 0.83333333333333 0.41675774729402 126.904 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 9.1809291 0.0000000 0.0000000 0.0000000 9.1809291 0.0000000 0.0000000 0.0000000 5.3931443 -------------------------- Born effective charges -------------------------- 1 Ge 4.5750768 0.0000000 0.0000000 0.0000000 4.5750768 0.0000000 0.0000000 0.0000000 1.4957523 2 Ge 4.5751309 0.0000000 0.0000000 0.0000000 4.5751309 0.0000000 0.0000000 0.0000000 1.4963455 3 Ge 4.5751309 0.0000000 0.0000000 0.0000000 4.5751309 0.0000000 0.0000000 0.0000000 1.4963455 4 I -2.2876063 -0.0000000 0.0000000 0.0000000 -2.2876063 0.0000000 0.0000000 0.0000000 -0.7477991 5 I -2.2876063 -0.0000000 0.0000000 0.0000000 -2.2876063 0.0000000 0.0000000 0.0000000 -0.7477991 6 I -2.2875520 -0.0000000 0.0000000 0.0000000 -2.2875520 0.0000000 0.0000000 0.0000000 -0.7478594 7 I -2.2875520 -0.0000000 0.0000000 0.0000000 -2.2875520 0.0000000 0.0000000 0.0000000 -0.7478594 8 I -2.2875110 -0.0000000 0.0000000 0.0000000 -2.2875110 0.0000000 0.0000000 0.0000000 -0.7485632 9 I -2.2875110 -0.0000000 0.0000000 0.0000000 -2.2875110 0.0000000 0.0000000 0.0000000 -0.7485632 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000007 (yyy) 0.00000007 (yyy) 0.00000007 (yyy) 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.18 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: 5.87e-04 5.87e-04 5.87e-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.325 ( 0.000 0.000 0.000) 0.000 0.325 ( 0.000 0.000 0.000) 0.000 0.327 ( 0.000 0.000 0.000) 0.000 0.327 ( 0.000 0.000 0.000) 0.000 0.594 ( 0.000 0.000 0.000) 0.000 0.606 ( 0.000 0.000 0.000) 0.000 2.395 ( 0.000 0.000 0.000) 0.000 2.395 ( 0.000 0.000 0.000) 0.000 2.395 ( 0.000 0.000 0.000) 0.000 2.395 ( 0.000 0.000 0.000) 0.000 2.405 ( 0.000 0.000 0.000) 0.000 2.405 ( 0.000 0.000 0.000) 0.000 2.555 ( 0.000 0.000 0.000) 0.000 2.555 ( 0.000 0.000 0.000) 0.000 2.555 ( 0.000 0.000 0.000) 0.000 2.555 ( 0.000 0.000 0.000) 0.000 2.580 ( 0.000 0.000 0.000) 0.000 2.580 ( 0.000 0.000 0.000) 0.000 2.978 ( 0.000 0.000 0.000) 0.000 2.978 ( 0.000 0.000 0.000) 0.000 3.037 ( 0.000 0.000 0.000) 0.000 5.279 ( 0.000 0.000 0.000) 0.000 5.562 ( 0.000 0.000 0.000) 0.000 5.562 ( 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: 5.87e-04 5.87e-04 5.87e-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.201 ( 9.225 5.326 0.038) 10.652 0.357 ( 15.062 8.696 5.292) 18.180 0.488 ( 12.019 6.939 -5.612) 14.970 0.519 ( 12.273 7.086 1.544) 14.256 0.574 ( 6.453 3.726 -4.848) 8.890 0.580 ( 24.305 14.033 -3.439) 28.275 0.586 ( 8.773 5.065 6.390) 11.977 0.685 ( 14.122 8.153 4.010) 16.793 0.715 ( 15.950 9.209 -2.277) 18.558 2.370 ( -1.213 -0.700 -0.189) 1.414 2.402 ( -1.166 -0.673 1.074) 1.722 2.412 ( 0.486 0.280 -0.895) 1.056 2.534 ( -1.601 -0.924 -0.175) 1.857 2.539 ( -0.583 -0.337 0.805) 1.049 2.552 ( -2.230 -1.288 -0.597) 2.644 2.555 ( -0.491 -0.283 -0.614) 0.836 2.556 ( 0.208 0.120 0.135) 0.275 2.557 ( -1.736 -1.002 0.361) 2.036 2.951 ( -0.995 -0.574 1.209) 1.668 2.951 ( -1.040 -0.600 -1.202) 1.699 3.011 ( -2.185 -1.261 0.025) 2.523 3.075 ( 41.153 23.759 14.761) 49.759 3.090 ( 41.038 23.693 -14.839) 49.655 4.769 ( -0.157 -0.090 -3.392) 3.397 5.184 ( -8.043 -4.643 3.793) 10.032 5.428 ( -11.280 -6.513 -0.871) 13.054 5.441 ( -9.690 -5.594 0.337) 11.194 ======================= Grid point 2 (3/48) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.426 ( 10.609 6.125 -0.601) 12.265 0.670 ( 4.138 2.389 -3.588) 5.975 0.673 ( 12.531 7.235 3.299) 14.840 0.686 ( 3.725 2.151 4.419) 6.167 0.766 ( 11.363 6.561 -4.889) 14.003 0.791 ( 10.790 6.230 1.837) 12.594 1.086 ( 19.847 11.459 -2.558) 23.059 1.110 ( 18.227 10.523 3.493) 21.335 1.133 ( 17.505 10.107 -1.232) 20.251 2.323 ( -3.344 -1.931 -0.994) 3.987 2.326 ( -4.920 -2.840 1.574) 5.895 2.389 ( -3.889 -2.245 -0.478) 4.516 2.497 ( -3.236 -1.868 0.222) 3.743 2.515 ( 0.591 0.341 0.371) 0.776 2.516 ( -2.862 -1.652 -0.478) 3.339 2.520 ( 0.846 0.488 -1.280) 1.610 2.539 ( -2.304 -1.330 0.158) 2.666 2.553 ( 1.603 0.926 0.843) 2.034 2.895 ( -2.855 -1.648 -1.086) 3.471 2.897 ( -2.585 -1.492 0.868) 3.108 2.944 ( -3.505 -2.024 0.152) 4.050 3.838 ( 21.559 12.447 -13.667) 28.400 3.861 ( 25.227 14.565 15.851) 33.164 4.708 ( -8.287 -4.784 -10.013) 13.850 4.973 ( -7.187 -4.149 11.348) 14.058 5.088 ( -17.716 -10.228 -2.807) 20.648 5.176 ( -12.164 -7.023 -0.659) 14.061 ======================= Grid point 3 (4/48) ======================= q-point: ( 0.21 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.682 ( 11.727 6.771 -0.792) 13.565 0.773 ( 5.063 2.923 0.583) 5.875 0.779 ( 4.788 2.765 0.363) 5.541 0.920 ( 9.123 5.267 0.956) 10.578 0.983 ( 7.640 4.411 -2.757) 9.243 0.994 ( 6.994 4.038 1.811) 8.276 1.463 ( 12.509 7.222 2.246) 14.618 1.466 ( 11.323 6.537 -1.041) 13.116 1.471 ( 13.899 8.024 -1.401) 16.110 2.215 ( -4.169 -2.407 0.813) 4.882 2.233 ( -3.837 -2.215 -0.946) 4.530 2.261 ( -6.184 -3.570 0.176) 7.143 2.414 ( -3.920 -2.263 0.281) 4.535 2.443 ( -3.151 -1.819 -0.542) 3.678 2.458 ( -4.125 -2.382 0.307) 4.773 2.557 ( 2.567 1.482 0.223) 2.973 2.572 ( 2.779 1.605 -0.890) 3.331 2.593 ( 1.545 0.892 0.641) 1.896 2.838 ( -1.851 -1.069 -0.168) 2.144 2.853 ( -1.043 -0.602 -0.062) 1.205 2.867 ( -2.934 -1.694 0.174) 3.393 4.031 ( -4.868 -2.811 -3.759) 6.763 4.197 ( 1.609 0.929 11.204) 11.358 4.342 ( -21.377 -12.342 -6.566) 25.542 4.701 ( -11.901 -6.871 -5.296) 14.728 4.921 ( 0.409 0.236 8.957) 8.969 4.975 ( -3.737 -2.157 -4.499) 6.234 ======================= Grid point 4 (5/48) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.911 ( 7.285 4.206 3.862) 9.256 0.918 ( 7.500 4.330 -3.763) 9.443 0.950 ( 11.616 6.707 -0.103) 13.414 1.089 ( 5.955 3.438 -0.492) 6.894 1.121 ( 4.953 2.860 -0.127) 5.721 1.122 ( 4.790 2.765 0.631) 5.567 1.645 ( 4.889 2.823 -0.497) 5.668 1.677 ( 6.546 3.779 0.094) 7.559 1.714 ( 7.415 4.281 0.379) 8.571 2.148 ( -3.067 -1.770 0.553) 3.584 2.158 ( -0.734 -0.424 -0.634) 1.058 2.173 ( -1.144 -0.661 0.044) 1.322 2.333 ( -2.706 -1.562 -0.028) 3.125 2.386 ( -1.801 -1.040 0.109) 2.083 2.389 ( -1.505 -0.869 -0.165) 1.745 2.605 ( 1.123 0.649 0.028) 1.297 2.611 ( -0.114 -0.066 0.072) 0.151 2.619 ( 0.920 0.531 -0.079) 1.066 2.820 ( 0.338 0.195 1.239) 1.299 2.822 ( -0.909 -0.525 -0.166) 1.063 2.849 ( 0.491 0.283 -0.975) 1.128 3.702 ( -20.897 -12.065 0.843) 24.145 3.806 ( -24.936 -14.397 -4.084) 29.081 3.869 ( -23.940 -13.822 3.699) 27.890 4.675 ( 4.243 2.450 -2.254) 5.393 4.939 ( 0.791 0.457 7.343) 7.400 4.982 ( 2.406 1.389 -5.546) 6.203 ======================= Grid point 5 (6/48) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.102 ( 9.018 5.206 3.810) 11.088 1.103 ( 8.173 4.718 -4.323) 10.380 1.194 ( 9.351 5.399 0.391) 10.805 1.206 ( 4.687 2.706 -0.384) 5.426 1.227 ( 4.506 2.601 0.985) 5.295 1.231 ( 5.018 2.897 -0.542) 5.820 1.721 ( 2.443 1.410 -0.133) 2.824 1.782 ( 3.257 1.881 -1.840) 4.187 1.816 ( 2.004 1.157 2.275) 3.245 2.131 ( 1.331 0.768 0.525) 1.624 2.173 ( 1.707 0.985 -1.436) 2.439 2.181 ( 1.637 0.945 0.815) 2.059 2.304 ( 0.105 0.061 -1.333) 1.338 2.355 ( -1.175 -0.678 1.392) 1.944 2.381 ( 0.329 0.190 -0.382) 0.538 2.585 ( -2.120 -1.224 -0.346) 2.472 2.596 ( -1.978 -1.142 0.161) 2.289 2.606 ( -2.007 -1.159 0.239) 2.330 2.821 ( 0.628 0.363 -0.518) 0.891 2.846 ( 1.633 0.943 2.032) 2.772 2.865 ( 0.713 0.411 -1.411) 1.633 3.188 ( -22.626 -13.063 1.032) 26.147 3.247 ( -23.188 -13.388 -2.756) 26.917 3.294 ( -24.809 -14.324 1.757) 28.701 4.767 ( 3.508 2.025 -2.550) 4.786 4.947 ( -0.209 -0.121 6.672) 6.676 5.035 ( 1.899 1.096 -4.186) 4.725 ======================= Grid point 6 (7/48) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.254 ( 4.577 2.642 -2.898) 6.027 1.288 ( 6.549 3.781 2.962) 8.122 1.305 ( 3.757 2.169 0.526) 4.370 1.320 ( 3.329 1.922 0.202) 3.850 1.342 ( 4.247 2.452 -0.702) 4.954 1.352 ( 4.172 2.409 -0.143) 4.820 1.771 ( 1.914 1.105 -0.404) 2.247 1.828 ( -0.557 -0.322 2.375) 2.460 1.840 ( 1.812 1.046 -1.890) 2.820 2.187 ( 2.973 1.717 0.670) 3.498 2.217 ( 1.802 1.041 -1.308) 2.458 2.236 ( 2.871 1.657 0.558) 3.362 2.326 ( 1.453 0.839 -1.582) 2.306 2.329 ( -0.866 -0.500 1.482) 1.787 2.383 ( -0.423 -0.244 0.064) 0.493 2.524 ( -3.857 -2.227 0.372) 4.470 2.526 ( -2.647 -1.528 -0.104) 3.058 2.540 ( -3.232 -1.866 -0.197) 3.737 2.745 ( -14.788 -8.538 0.560) 17.085 2.802 ( -13.950 -8.054 -2.106) 16.245 2.809 ( -14.805 -8.547 1.341) 17.147 2.844 ( 0.520 0.300 -0.055) 0.603 2.879 ( -0.910 -0.525 -1.144) 1.554 2.883 ( 1.288 0.744 1.376) 2.027 4.835 ( 2.600 1.501 -3.537) 4.639 4.929 ( -1.345 -0.776 5.836) 6.039 5.064 ( 0.764 0.441 -2.299) 2.462 ======================= Grid point 7 (8/48) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.305 ( 0.000 0.000 0.000) 0.000 1.356 ( -0.000 -0.000 -0.000) 0.000 1.358 ( 0.000 0.000 0.000) 0.000 1.381 ( -1.486 -0.858 -1.387) 2.207 1.381 ( 1.486 0.858 1.387) 2.207 1.396 ( 0.000 0.000 0.000) 0.000 1.805 ( -0.000 -0.000 0.000) 0.000 1.805 ( 0.000 0.000 -0.000) 0.000 1.860 ( 0.000 0.000 -0.000) 0.000 2.235 ( 0.001 0.000 0.001) 0.002 2.235 ( -0.001 -0.000 -0.001) 0.002 2.281 ( 0.000 0.000 0.000) 0.000 2.322 ( 0.000 0.000 -0.000) 0.000 2.360 ( 0.000 0.000 -0.000) 0.001 2.361 ( -0.000 -0.000 0.000) 0.001 2.464 ( -0.000 -0.000 0.000) 0.000 2.492 ( -0.324 -0.187 -0.528) 0.647 2.492 ( 0.324 0.187 0.528) 0.647 2.566 ( -0.000 -0.000 -0.000) 0.000 2.630 ( -0.000 -0.000 -0.000) 0.000 2.631 ( 0.000 0.000 0.000) 0.000 2.866 ( -0.826 -0.477 -0.507) 1.080 2.866 ( 0.826 0.477 0.507) 1.080 2.899 ( 0.000 0.000 0.000) 0.000 4.888 ( -0.000 -0.000 0.000) 0.000 4.889 ( 0.000 0.000 -0.000) 0.000 5.073 ( 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: 5.87e-04 5.87e-04 5.87e-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.355 ( 5.569 9.646 0.000) 11.139 0.622 ( 8.826 15.287 0.000) 17.652 0.641 ( 2.200 3.811 0.000) 4.400 0.647 ( 2.156 3.734 0.000) 4.312 0.742 ( 7.994 13.846 0.000) 15.988 0.743 ( 7.942 13.756 0.000) 15.884 0.945 ( 11.582 20.060 0.000) 23.163 0.990 ( 9.993 17.309 0.000) 19.987 0.990 ( 9.993 17.308 0.000) 19.985 2.382 ( -0.360 -0.623 0.000) 0.720 2.382 ( -0.358 -0.620 0.000) 0.715 2.423 ( 0.345 0.597 0.000) 0.689 2.480 ( -1.603 -2.776 0.000) 3.206 2.480 ( -1.599 -2.769 0.000) 3.197 2.509 ( -1.590 -2.753 0.000) 3.179 2.514 ( -1.737 -3.009 0.000) 3.474 2.560 ( -0.529 -0.916 0.000) 1.057 2.560 ( -0.524 -0.908 0.000) 1.049 2.911 ( -1.671 -2.893 0.000) 3.341 2.911 ( -1.669 -2.890 0.000) 3.337 2.964 ( -1.967 -3.408 0.000) 3.935 3.693 ( 17.263 29.900 0.000) 34.525 3.694 ( 17.266 29.906 0.000) 34.533 4.772 ( -0.088 -0.152 0.000) 0.176 4.998 ( -7.783 -13.480 0.000) 15.566 5.218 ( -8.643 -14.970 0.000) 17.286 5.219 ( -8.646 -14.975 0.000) 17.291 ======================= Grid point 17 (10/48) ======================= q-point: ( 0.14 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.569 ( 8.924 8.331 -0.483) 12.218 0.723 ( 3.699 3.042 -0.276) 4.797 0.727 ( 3.588 2.771 0.851) 4.612 0.904 ( 6.381 14.531 1.131) 15.911 0.983 ( 4.374 13.701 -3.567) 14.818 0.999 ( 4.787 12.953 2.389) 14.014 1.311 ( 12.248 9.025 2.399) 15.402 1.315 ( 13.220 10.854 -1.187) 17.147 1.322 ( 12.289 8.323 -1.262) 14.896 2.310 ( -6.523 1.234 1.023) 6.717 2.337 ( -6.283 2.547 -1.062) 6.863 2.373 ( -7.148 0.693 0.240) 7.186 2.442 ( 1.322 -3.615 -1.004) 3.978 2.448 ( -2.710 -2.636 0.287) 3.791 2.476 ( 0.516 -3.713 0.875) 3.849 2.507 ( 4.359 -1.969 -1.044) 4.896 2.527 ( -0.264 -0.592 0.537) 0.841 2.546 ( 1.781 -1.687 0.171) 2.459 2.849 ( -1.898 -2.797 -0.625) 3.437 2.854 ( -1.262 -2.715 0.474) 3.032 2.887 ( -2.616 -3.686 0.104) 4.521 4.090 ( 3.004 11.661 -8.239) 14.590 4.155 ( 11.154 13.109 11.678) 20.800 4.595 ( -19.265 -8.054 -8.218) 22.440 4.832 ( -15.901 -15.080 -2.718) 22.082 4.832 ( 4.240 -8.595 9.021) 13.162 4.948 ( -4.076 -15.880 -1.495) 16.463 ======================= Grid point 18 (11/48) ======================= q-point: ( 0.21 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.814 ( 10.707 6.915 -0.526) 12.757 0.832 ( 5.928 3.330 -1.356) 6.933 0.835 ( 5.678 3.509 1.957) 6.956 1.116 ( 2.455 13.333 0.182) 13.559 1.152 ( 0.945 11.697 -1.311) 11.808 1.161 ( 0.541 11.685 1.042) 11.744 1.556 ( 8.614 2.834 -0.609) 9.088 1.564 ( 10.098 3.164 0.985) 10.627 1.596 ( 10.618 4.605 -0.428) 11.581 2.228 ( -5.145 3.081 -0.054) 5.997 2.248 ( -7.627 1.295 0.578) 7.757 2.257 ( -5.881 3.610 -0.446) 6.915 2.376 ( -3.305 -0.668 0.270) 3.382 2.411 ( -1.677 -0.922 -0.432) 1.961 2.419 ( -2.595 -1.026 0.084) 2.792 2.560 ( 3.549 -0.868 0.309) 3.666 2.569 ( 4.170 -1.666 -0.567) 4.526 2.582 ( 2.793 -1.415 0.272) 3.143 2.806 ( -0.045 -2.051 0.535) 2.120 2.818 ( -1.030 -3.158 -0.050) 3.323 2.826 ( 0.815 -2.069 -0.462) 2.272 3.975 ( -19.294 -2.973 0.145) 19.522 4.114 ( -24.725 -10.944 -4.641) 27.435 4.182 ( -21.750 -2.937 6.307) 22.836 4.552 ( 6.349 -8.033 -3.098) 10.697 4.853 ( 5.303 -6.488 6.844) 10.819 4.857 ( 7.299 -9.656 -5.543) 13.314 ======================= Grid point 19 (12/48) ======================= q-point: ( 0.29 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.994 ( 8.331 4.383 3.626) 10.088 0.994 ( 8.449 3.940 -3.602) 9.994 1.058 ( 10.475 5.066 -0.188) 11.637 1.258 ( 0.044 11.880 -0.085) 11.880 1.263 ( -0.200 10.610 0.691) 10.635 1.269 ( -0.311 10.860 -0.520) 10.877 1.683 ( 4.322 1.649 -0.114) 4.628 1.722 ( 6.119 1.559 -0.936) 6.384 1.766 ( 5.781 1.689 1.201) 6.141 2.160 ( -3.627 1.182 0.558) 3.856 2.200 ( -1.942 3.540 -1.293) 4.239 2.214 ( -3.006 3.021 0.690) 4.318 2.340 ( -2.022 2.938 -0.435) 3.592 2.385 ( -1.621 1.504 0.507) 2.269 2.395 ( -0.736 1.677 -0.271) 1.851 2.594 ( 1.047 -1.212 0.363) 1.642 2.595 ( 0.042 -1.059 -0.288) 1.099 2.599 ( 1.206 -2.055 -0.040) 2.383 2.795 ( 1.240 -2.197 -0.429) 2.559 2.810 ( 2.198 -1.219 1.682) 3.024 2.835 ( 2.032 -1.661 -1.159) 2.869 3.509 ( -26.348 -7.328 1.454) 27.387 3.567 ( -26.651 -9.529 -3.213) 28.485 3.643 ( -29.252 -9.072 1.918) 30.686 4.650 ( 8.612 -4.591 -1.608) 9.891 4.874 ( 4.591 -6.699 6.396) 10.337 4.918 ( 8.213 -7.468 -4.948) 12.153 ======================= Grid point 20 (13/48) ======================= q-point: ( 0.36 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.172 ( 7.525 2.833 -3.276) 8.682 1.186 ( 8.698 3.623 3.125) 9.927 1.255 ( 7.439 3.363 -0.167) 8.165 1.360 ( -0.389 9.625 0.772) 9.663 1.361 ( -0.417 9.508 0.526) 9.531 1.371 ( 0.088 9.875 -1.038) 9.930 1.753 ( 2.090 2.019 -0.103) 2.908 1.814 ( 2.981 1.579 -1.977) 3.910 1.829 ( 0.716 0.310 2.305) 2.434 2.164 ( 2.439 1.840 -0.252) 3.066 2.215 ( 0.905 0.752 1.391) 1.822 2.226 ( 0.980 2.994 -1.318) 3.414 2.357 ( -1.264 5.222 -0.875) 5.444 2.381 ( -2.061 4.005 0.999) 4.613 2.409 ( -0.927 2.884 -0.207) 3.036 2.566 ( -2.408 -0.267 -0.449) 2.464 2.568 ( -2.528 -1.134 0.601) 2.835 2.572 ( -1.847 -1.679 -0.103) 2.499 2.809 ( 2.213 -1.474 -0.574) 2.720 2.844 ( 2.587 -1.123 1.897) 3.399 2.853 ( 1.727 -1.633 -1.299) 2.709 3.006 ( -23.600 -4.768 0.639) 24.086 3.056 ( -23.137 -5.702 -2.033) 23.916 3.082 ( -25.572 -6.622 1.390) 26.452 4.736 ( 7.042 -4.882 -2.647) 8.969 4.873 ( 3.247 -7.080 5.901) 9.772 4.971 ( 6.398 -7.231 -3.277) 10.196 ======================= Grid point 21 (14/48) ======================= q-point: ( 0.43 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.281 ( 3.027 0.459 -1.292) 3.323 1.334 ( 5.087 1.385 1.811) 5.575 1.365 ( 2.596 0.389 -0.573) 2.686 1.418 ( -2.500 7.088 0.780) 7.557 1.422 ( -2.576 7.281 -0.350) 7.731 1.450 ( -1.636 7.595 -0.406) 7.780 1.795 ( 0.987 1.273 -0.812) 1.804 1.820 ( -1.045 -0.404 1.780) 2.103 1.855 ( 0.750 0.405 -0.961) 1.285 2.258 ( 3.965 -1.411 0.703) 4.267 2.263 ( 5.634 2.857 -0.601) 6.345 2.290 ( 3.249 2.153 -0.125) 3.900 2.370 ( -4.572 5.731 -0.096) 7.332 2.391 ( -3.679 7.142 0.295) 8.039 2.393 ( -6.064 3.983 -0.233) 7.259 2.506 ( -3.319 0.575 0.504) 3.406 2.526 ( -2.200 1.617 0.164) 2.736 2.526 ( -1.965 0.825 -0.592) 2.212 2.656 ( -12.316 0.692 0.203) 12.337 2.709 ( -11.766 -0.804 1.283) 11.863 2.710 ( -11.215 -0.615 -1.486) 11.329 2.840 ( 1.682 -0.783 0.169) 1.863 2.856 ( 0.579 -1.758 -0.851) 2.038 2.879 ( 1.527 -0.969 0.663) 1.926 4.799 ( 6.125 -5.045 -3.815) 8.805 4.846 ( 2.032 -7.429 4.988) 9.176 4.995 ( 4.882 -7.324 -1.170) 8.879 ======================= Grid point 31 (15/48) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.752 ( 5.864 10.157 0.000) 11.728 0.798 ( 2.810 4.866 0.000) 5.619 0.800 ( 2.759 4.779 0.000) 5.518 1.186 ( 7.538 13.057 0.000) 15.076 1.243 ( 6.427 11.132 0.000) 12.854 1.244 ( 6.430 11.137 0.000) 12.860 1.462 ( 3.867 6.698 0.000) 7.734 1.462 ( 3.875 6.712 0.000) 7.750 1.505 ( 5.049 8.745 0.000) 10.098 2.341 ( -0.805 -1.394 0.000) 1.610 2.341 ( -0.804 -1.393 0.000) 1.609 2.344 ( -2.594 -4.493 0.000) 5.188 2.423 ( 0.177 0.307 0.000) 0.354 2.423 ( 0.172 0.297 0.000) 0.343 2.450 ( 0.253 0.439 0.000) 0.507 2.461 ( 0.313 0.543 0.000) 0.627 2.518 ( 0.130 0.225 0.000) 0.259 2.518 ( 0.129 0.224 0.000) 0.259 2.798 ( -1.333 -2.308 0.000) 2.666 2.799 ( -1.331 -2.305 0.000) 2.662 2.809 ( -2.226 -3.856 0.000) 4.452 4.274 ( -0.129 -0.224 0.000) 0.258 4.274 ( -0.145 -0.252 0.000) 0.291 4.327 ( -10.849 -18.790 0.000) 21.697 4.568 ( -7.762 -13.443 0.000) 15.523 4.569 ( -7.754 -13.430 0.000) 15.507 4.721 ( -1.557 -2.697 0.000) 3.114 ======================= Grid point 32 (16/48) ======================= q-point: ( 0.21 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.919 ( 5.411 5.382 -2.462) 8.020 0.921 ( 5.506 5.268 2.604) 8.053 0.959 ( 7.755 7.712 -0.180) 10.938 1.392 ( 1.049 13.566 0.034) 13.606 1.402 ( -0.882 12.913 -0.310) 12.947 1.408 ( -0.676 12.578 0.234) 12.599 1.600 ( 6.993 2.169 0.337) 7.330 1.609 ( 8.291 1.990 -0.511) 8.542 1.668 ( 7.145 3.155 0.231) 7.814 2.239 ( -4.778 -2.539 0.330) 5.421 2.282 ( -3.761 1.033 -1.239) 4.092 2.299 ( -3.499 -0.604 1.049) 3.703 2.423 ( -3.265 4.590 0.399) 5.647 2.432 ( -0.888 2.383 0.297) 2.560 2.433 ( -0.958 2.984 -0.875) 3.254 2.527 ( 4.539 -1.951 -0.544) 4.970 2.546 ( 2.712 0.087 1.062) 2.913 2.559 ( 2.851 -0.415 -0.510) 2.926 2.751 ( -0.027 -3.233 -0.197) 3.240 2.765 ( 0.475 -1.810 1.015) 2.129 2.778 ( 1.364 -2.544 -0.786) 2.991 3.902 ( -20.739 -10.946 -2.391) 23.572 3.906 ( -25.144 -4.443 2.854) 25.692 4.050 ( -23.845 -11.137 -0.996) 26.336 4.424 ( 7.660 -3.806 2.701) 8.970 4.621 ( 13.443 -13.099 -5.815) 19.650 4.697 ( 4.490 -8.429 3.650) 10.224 ======================= Grid point 33 (17/48) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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| 1.077 ( 6.987 4.342 -2.687) 8.654 1.084 ( 7.482 4.571 2.608) 9.147 1.157 ( 7.801 4.788 -0.039) 9.153 1.495 ( -2.125 11.965 0.637) 12.169 1.503 ( -2.280 12.052 0.234) 12.268 1.503 ( -1.868 11.977 -0.814) 12.149 1.725 ( 4.451 2.756 0.226) 5.241 1.757 ( 6.066 2.079 -1.279) 6.539 1.795 ( 4.402 1.223 1.222) 4.729 2.159 ( -1.584 -1.055 -0.309) 1.928 2.230 ( -2.891 -1.509 1.731) 3.692 2.241 ( -1.383 -0.192 -1.517) 2.062 2.435 ( -1.294 5.437 -0.006) 5.589 2.448 ( -1.062 3.709 -0.168) 3.862 2.452 ( -0.332 3.365 0.083) 3.383 2.561 ( 1.312 -1.305 -0.367) 1.887 2.577 ( 0.469 -0.309 1.200) 1.325 2.585 ( 0.234 0.103 -0.796) 0.837 2.746 ( 2.394 -2.243 -0.419) 3.307 2.779 ( 2.818 -1.660 1.447) 3.576 2.789 ( 2.654 -2.633 -1.026) 3.877 3.426 ( -26.580 -0.674 1.883) 26.655 3.438 ( -24.262 -3.390 -2.225) 24.599 3.522 ( -28.353 -2.735 0.380) 28.487 4.514 ( 9.819 -8.646 -0.528) 13.094 4.694 ( 6.073 -11.019 4.628) 13.406 4.710 ( 11.798 -12.907 -4.131) 17.968 ======================= Grid point 34 (18/48) ======================= q-point: ( 0.36 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.218 ( 5.230 1.838 -1.799) 5.828 1.246 ( 6.833 2.315 1.792) 7.433 1.303 ( 5.149 1.292 -0.045) 5.309 1.561 ( -3.467 10.113 0.660) 10.711 1.563 ( -3.011 10.137 -0.075) 10.575 1.580 ( -2.376 10.590 -0.554) 10.868 1.797 ( 1.026 2.279 -0.360) 2.525 1.831 ( -0.274 0.035 1.601) 1.625 1.842 ( 1.854 1.067 -1.177) 2.441 2.179 ( 4.077 -0.578 -1.177) 4.282 2.191 ( 1.042 -2.533 1.378) 3.066 2.234 ( 1.755 -2.144 -0.286) 2.785 2.467 ( -0.631 2.959 -0.052) 3.026 2.473 ( -0.602 4.836 -0.166) 4.876 2.479 ( -0.239 2.792 0.294) 2.818 2.555 ( -0.942 -0.139 -0.371) 1.022 2.556 ( -1.855 -0.549 0.720) 2.064 2.575 ( -1.430 0.450 -0.359) 1.542 2.775 ( 2.598 -1.536 -0.368) 3.041 2.804 ( 1.837 -2.689 -0.712) 3.333 2.811 ( 2.149 -1.771 1.032) 2.970 3.042 ( -21.555 9.044 0.851) 23.391 3.052 ( -21.143 6.305 -1.181) 22.094 3.083 ( -24.273 7.661 0.341) 25.456 4.588 ( 9.526 -9.863 -2.293) 13.903 4.680 ( 5.596 -12.181 4.488) 14.136 4.766 ( 9.762 -13.101 -2.207) 16.486 ======================= Grid point 35 (19/48) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.274 ( 0.649 -1.124 -0.000) 1.298 1.338 ( 1.027 -1.779 0.000) 2.054 1.338 ( 1.021 -1.768 -0.000) 2.042 1.585 ( -5.003 8.665 0.000) 10.006 1.585 ( -4.987 8.637 -0.000) 9.974 1.616 ( -4.972 8.611 -0.000) 9.943 1.817 ( -0.404 0.699 -0.000) 0.807 1.818 ( -0.406 0.703 0.000) 0.812 1.856 ( 0.117 -0.202 0.000) 0.233 2.197 ( 2.086 -3.612 0.000) 4.171 2.235 ( 3.362 -3.201 -1.029) 4.754 2.236 ( 1.091 -4.511 1.029) 4.754 2.481 ( -0.760 1.317 -0.000) 1.521 2.494 ( -1.339 2.319 0.000) 2.677 2.494 ( -1.338 2.317 -0.000) 2.676 2.521 ( -0.444 0.769 -0.000) 0.888 2.553 ( -0.495 0.858 0.000) 0.990 2.553 ( -0.495 0.858 -0.000) 0.991 2.792 ( -0.734 -0.305 0.165) 0.812 2.792 ( 0.630 0.485 -0.165) 0.812 2.806 ( -1.976 3.423 0.000) 3.952 2.903 ( -9.556 16.550 -0.000) 19.111 2.903 ( -9.574 16.584 0.000) 19.149 2.913 ( -8.147 14.110 -0.000) 16.293 4.644 ( 6.698 -11.601 -0.001) 13.396 4.645 ( 6.710 -11.623 0.001) 13.421 4.785 ( 7.910 -13.701 -0.000) 15.820 ======================= Grid point 47 (20/48) ======================= q-point: ( 0.21 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.038 ( 3.740 6.479 0.000) 7.481 1.038 ( 3.748 6.491 0.000) 7.496 1.110 ( 4.246 7.355 0.000) 8.493 1.615 ( 4.715 8.166 0.000) 9.429 1.631 ( 4.912 8.507 0.000) 9.824 1.631 ( 4.905 8.496 0.000) 9.810 1.656 ( 2.602 4.507 0.000) 5.204 1.657 ( 2.600 4.503 0.000) 5.200 1.733 ( 2.240 3.880 0.000) 4.480 2.175 ( -1.966 -3.406 0.000) 3.933 2.256 ( -1.836 -3.180 0.000) 3.672 2.256 ( -1.831 -3.171 0.000) 3.661 2.478 ( 0.846 1.466 0.000) 1.693 2.478 ( 0.851 1.474 0.000) 1.701 2.502 ( 0.896 1.553 0.000) 1.793 2.512 ( 1.011 1.752 0.000) 2.022 2.570 ( 0.755 1.309 0.000) 1.511 2.570 ( 0.757 1.311 0.000) 1.514 2.710 ( -0.397 -0.688 0.000) 0.794 2.743 ( -0.265 -0.459 0.000) 0.530 2.743 ( -0.264 -0.457 0.000) 0.528 3.672 ( -7.134 -12.356 0.000) 14.267 3.771 ( -7.388 -12.796 0.000) 14.775 3.771 ( -7.388 -12.796 0.000) 14.775 4.403 ( -1.903 -3.295 0.000) 3.805 4.403 ( -1.908 -3.305 0.000) 3.817 4.551 ( -3.471 -6.012 0.000) 6.942 ======================= Grid point 48 (21/48) ======================= q-point: ( 0.29 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.161 ( 4.228 4.135 -1.632) 6.135 1.169 ( 5.160 4.039 1.539) 6.731 1.240 ( 4.305 3.493 0.082) 5.545 1.700 ( -2.020 8.014 0.326) 8.271 1.705 ( -2.080 7.985 0.297) 8.257 1.711 ( -1.151 8.321 -0.621) 8.423 1.789 ( 3.799 3.689 -0.186) 5.299 1.801 ( 4.954 2.388 -0.366) 5.512 1.814 ( 2.988 0.894 0.620) 3.180 2.131 ( 0.340 -1.482 -0.572) 1.625 2.179 ( -2.453 -2.992 1.555) 4.170 2.198 ( -0.651 -3.527 -1.058) 3.739 2.495 ( -0.273 1.160 -0.627) 1.346 2.497 ( 0.051 1.109 0.562) 1.244 2.515 ( -0.711 2.272 0.219) 2.390 2.549 ( 1.182 -0.268 -0.557) 1.334 2.570 ( -1.049 -0.973 1.110) 1.811 2.582 ( -0.029 -1.096 -0.679) 1.290 2.723 ( 1.940 0.052 -0.145) 1.946 2.751 ( 1.530 -0.777 -0.498) 1.787 2.755 ( 2.015 -0.558 0.619) 2.181 3.437 ( -13.932 3.159 -1.083) 14.326 3.487 ( -17.541 6.324 2.483) 18.811 3.538 ( -16.226 3.799 -1.443) 16.727 4.335 ( 4.913 -8.956 0.752) 10.243 4.437 ( 10.538 -14.153 -3.029) 17.904 4.461 ( 3.129 -12.099 2.319) 12.711 ======================= Grid point 49 (22/48) ======================= q-point: ( 0.36 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.246 ( 2.175 1.085 -0.672) 2.522 1.277 ( 4.040 0.990 0.742) 4.226 1.308 ( 1.595 -0.564 -0.057) 1.693 1.732 ( -2.736 6.478 -0.526) 7.052 1.732 ( -2.965 6.610 0.763) 7.285 1.762 ( -2.198 7.254 -0.237) 7.583 1.837 ( -0.487 0.962 0.492) 1.185 1.837 ( -0.616 1.581 -0.499) 1.768 1.851 ( 0.271 -0.123 -0.005) 0.297 2.143 ( 0.803 -1.820 0.543) 2.062 2.154 ( 3.053 -1.425 -0.991) 3.512 2.176 ( 1.315 -2.965 0.424) 3.271 2.496 ( -0.545 0.598 -0.385) 0.896 2.502 ( 0.112 -0.026 0.431) 0.446 2.512 ( -0.733 -0.295 -0.096) 0.796 2.530 ( -1.021 -1.643 0.436) 1.983 2.554 ( 0.231 -0.381 -0.777) 0.895 2.561 ( -0.247 -1.407 0.434) 1.493 2.758 ( 1.739 -0.281 -0.123) 1.765 2.768 ( 1.010 -1.014 -0.091) 1.435 2.782 ( 1.663 -1.103 0.192) 2.004 3.303 ( -12.975 17.737 -0.551) 21.983 3.326 ( -15.069 17.962 1.024) 23.468 3.347 ( -17.270 17.020 -0.464) 24.251 4.359 ( 10.007 -13.085 -1.687) 16.559 4.406 ( 5.974 -15.379 2.750) 16.726 4.469 ( 10.784 -16.818 -1.076) 20.008 ======================= Grid point 62 (23/48) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.237 ( 1.761 3.050 0.000) 3.522 1.238 ( 1.766 3.060 0.000) 3.533 1.289 ( 0.788 1.364 0.000) 1.575 1.798 ( 0.951 1.647 0.000) 1.902 1.803 ( 1.367 2.368 0.000) 2.734 1.803 ( 1.360 2.356 0.000) 2.721 1.841 ( 1.058 1.833 0.000) 2.116 1.852 ( 1.072 1.856 0.000) 2.144 1.852 ( 1.072 1.857 0.000) 2.145 2.118 ( 0.221 0.384 0.000) 0.443 2.138 ( -0.749 -1.297 0.000) 1.498 2.138 ( -0.739 -1.279 0.000) 1.477 2.504 ( -0.090 -0.156 0.000) 0.180 2.505 ( -0.096 -0.167 0.000) 0.193 2.522 ( -0.734 -1.272 0.000) 1.469 2.540 ( -1.086 -1.881 0.000) 2.171 2.540 ( -1.088 -1.884 0.000) 2.176 2.544 ( -0.107 -0.186 0.000) 0.215 2.740 ( 0.731 1.267 0.000) 1.463 2.756 ( 0.349 0.605 0.000) 0.698 2.756 ( 0.346 0.598 0.000) 0.691 3.549 ( 4.354 7.542 0.000) 8.708 3.635 ( 3.496 6.055 0.000) 6.991 3.635 ( 3.499 6.060 0.000) 6.998 4.182 ( -4.855 -8.410 0.000) 9.711 4.182 ( -4.860 -8.418 0.000) 9.721 4.236 ( -5.975 -10.350 0.000) 11.951 ======================= Grid point 63 (24/48) ======================= q-point: ( 0.36 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.265 ( -0.576 0.998 -0.000) 1.153 1.288 ( 0.315 -0.545 0.000) 0.630 1.288 ( 0.311 -0.541 -0.000) 0.624 1.812 ( -1.153 1.996 -0.000) 2.305 1.813 ( -1.158 2.006 0.000) 2.317 1.839 ( 0.438 -0.758 -0.000) 0.876 1.857 ( -1.281 2.219 0.000) 2.562 1.863 ( -0.596 1.032 -0.000) 1.191 1.864 ( -0.587 1.017 0.000) 1.174 2.125 ( 0.063 -0.110 -0.000) 0.127 2.132 ( 0.587 -1.017 0.000) 1.174 2.132 ( 0.588 -1.017 -0.000) 1.175 2.498 ( -0.213 -0.492 0.001) 0.536 2.498 ( 0.532 -0.061 -0.001) 0.536 2.503 ( 0.530 -0.918 -0.000) 1.060 2.510 ( -0.625 1.082 -0.000) 1.249 2.537 ( 0.243 -1.312 -0.473) 1.416 2.537 ( 1.015 -0.866 0.473) 1.416 2.758 ( -0.369 -0.103 0.080) 0.392 2.759 ( 0.273 0.269 -0.080) 0.392 2.767 ( 0.324 -0.561 -0.000) 0.647 3.684 ( -9.997 17.292 -0.000) 19.974 3.684 ( -9.973 17.298 0.000) 19.967 3.691 ( -10.260 17.770 -0.000) 20.519 4.088 ( 6.111 -16.168 0.817) 17.304 4.088 ( 10.946 -13.376 -0.817) 17.304 4.117 ( 10.504 -18.194 -0.000) 21.008 ======================= Grid point 211 (25/48) ======================= q-point: ( 0.00 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.191 ( 0.000 0.000 0.000) 0.000 0.191 ( 0.000 -0.000 0.000) 0.000 0.193 ( -0.000 0.000 -0.000) 0.000 0.193 ( 0.000 0.000 -0.000) 0.000 0.346 ( -0.000 0.000 0.000) 0.000 0.354 ( 0.000 0.000 -0.000) 0.000 0.373 ( 0.000 -0.000 0.000) 0.000 0.373 ( 0.000 0.000 0.000) 0.000 0.688 ( 0.000 0.000 0.000) 0.000 2.392 ( 0.000 0.000 0.000) 0.000 2.392 ( 0.000 0.000 0.000) 0.000 2.402 ( 0.000 -0.000 -0.244) 0.244 2.402 ( -0.000 0.000 0.244) 0.244 2.402 ( 0.000 0.000 -0.244) 0.244 2.402 ( 0.000 0.000 0.244) 0.244 2.547 ( 0.000 0.000 0.000) 0.000 2.547 ( 0.000 -0.000 0.000) 0.000 2.571 ( -0.000 0.000 -0.000) 0.000 2.571 ( 0.000 0.000 -0.000) 0.000 2.571 ( 0.000 0.000 0.000) 0.000 2.571 ( -0.000 0.000 0.000) 0.000 2.960 ( 0.000 0.000 0.000) 0.000 3.016 ( 0.000 0.000 -0.000) 0.000 3.017 ( 0.000 0.000 0.000) 0.000 5.561 ( -0.000 0.000 0.000) 0.000 5.563 ( 0.000 0.000 -0.000) 0.000 5.564 ( 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: 5.87e-04 5.87e-04 5.87e-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.332 ( 13.825 7.982 -4.708) 16.644 0.383 ( 5.653 3.264 3.359) 7.341 0.389 ( 5.262 3.038 -3.718) 7.123 0.471 ( 13.324 7.692 3.169) 15.708 0.549 ( 12.226 7.059 0.218) 14.119 0.559 ( 21.410 12.361 3.345) 24.947 0.654 ( 21.770 12.569 -2.446) 25.256 0.689 ( 11.900 6.870 -3.230) 14.115 0.715 ( 9.941 5.740 2.829) 11.823 2.379 ( -1.487 -0.859 0.781) 1.886 2.387 ( -0.338 -0.195 -1.046) 1.117 2.422 ( 0.227 0.131 0.277) 0.381 2.528 ( -1.101 -0.635 -0.073) 1.273 2.538 ( -2.285 -1.319 -0.346) 2.661 2.550 ( 0.025 0.014 0.222) 0.224 2.554 ( -0.941 -0.544 -0.152) 1.098 2.555 ( -1.304 -0.753 0.385) 1.554 2.565 ( -1.127 -0.651 0.049) 1.303 2.902 ( 27.414 15.827 -0.365) 31.657 2.947 ( 6.647 3.838 -0.136) 7.676 2.991 ( -2.152 -1.243 1.332) 2.820 2.992 ( -2.059 -1.189 -1.333) 2.726 3.750 ( 46.329 26.748 -42.500) 68.324 3.754 ( 48.434 27.963 44.102) 71.223 5.380 ( -13.964 -8.062 -4.143) 16.648 5.407 ( -10.898 -6.292 3.480) 13.056 5.440 ( -10.061 -5.809 -0.267) 11.620 ======================= Grid point 215 (27/48) ======================= q-point: ( 0.14 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.532 ( 7.526 4.345 7.092) 11.217 0.542 ( 8.284 4.783 -6.496) 11.562 0.653 ( 13.433 7.755 -2.634) 15.733 0.735 ( 1.358 0.784 -0.840) 1.778 0.747 ( 10.643 6.145 2.260) 12.496 0.825 ( 10.888 6.286 0.121) 12.573 1.042 ( 19.657 11.349 0.449) 22.703 1.116 ( 18.249 10.536 -0.255) 21.073 1.152 ( 17.949 10.363 0.098) 20.726 2.310 ( -4.355 -2.515 -0.088) 5.030 2.358 ( -3.256 -1.880 -1.732) 4.140 2.372 ( -4.729 -2.730 1.717) 5.724 2.506 ( 0.652 0.376 0.130) 0.764 2.517 ( -2.936 -1.695 0.721) 3.466 2.518 ( -2.506 -1.447 -1.508) 3.262 2.520 ( -2.736 -1.580 0.884) 3.281 2.527 ( 1.211 0.699 0.891) 1.658 2.552 ( 1.351 0.780 -0.954) 1.828 2.884 ( -2.313 -1.336 0.124) 2.674 2.924 ( -3.382 -1.952 1.254) 4.101 2.929 ( -3.273 -1.889 -1.311) 4.000 3.670 ( 26.567 15.339 0.125) 30.677 4.280 ( 6.720 3.880 -21.118) 22.499 4.383 ( 14.350 8.285 24.912) 29.919 4.996 ( -19.169 -11.067 -5.200) 22.737 5.143 ( -14.988 -8.654 -1.906) 17.412 5.156 ( -9.675 -5.586 3.135) 11.604 ======================= Grid point 217 (28/48) ======================= q-point: ( 0.21 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.728 ( 9.652 5.572 3.405) 11.653 0.747 ( 9.348 5.397 -3.021) 11.209 0.766 ( 1.867 1.078 -0.538) 2.222 0.923 ( 10.120 5.843 -1.507) 11.782 0.952 ( 7.484 4.321 1.545) 8.779 1.021 ( 6.270 3.620 -0.048) 7.240 1.430 ( 14.247 8.226 -0.585) 16.461 1.466 ( 12.235 7.064 1.057) 14.167 1.489 ( 11.335 6.545 -0.280) 13.092 2.212 ( -3.506 -2.024 -0.565) 4.087 2.243 ( -5.759 -3.325 1.195) 6.756 2.253 ( -4.970 -2.870 -0.673) 5.779 2.440 ( -3.991 -2.304 -1.683) 4.907 2.441 ( -3.494 -2.017 1.023) 4.163 2.442 ( -3.465 -2.000 0.614) 4.047 2.558 ( 2.992 1.728 -0.276) 3.466 2.571 ( 2.200 1.270 0.899) 2.694 2.593 ( 1.730 0.999 -0.597) 2.086 2.847 ( -0.727 -0.420 0.582) 1.022 2.855 ( -2.407 -1.390 0.208) 2.787 2.856 ( -2.745 -1.585 -0.733) 3.254 4.019 ( 2.017 1.165 3.265) 4.011 4.171 ( -14.568 -8.411 -7.161) 18.282 4.439 ( -18.514 -10.689 2.210) 21.492 4.665 ( 0.980 0.566 7.234) 7.322 4.842 ( -8.443 -4.875 -6.094) 11.498 5.031 ( -1.263 -0.729 0.506) 1.544 ======================= Grid point 219 (29/48) ======================= q-point: ( 0.29 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.849 ( 6.008 3.469 -0.284) 6.943 0.962 ( 10.013 5.781 0.136) 11.562 0.963 ( 10.180 5.877 0.152) 11.756 1.097 ( 5.623 3.246 1.043) 6.576 1.110 ( 6.470 3.736 -0.864) 7.521 1.125 ( 3.580 2.067 -0.192) 4.138 1.668 ( 5.622 3.246 2.142) 6.836 1.668 ( 4.913 2.836 -0.719) 5.718 1.688 ( 8.548 4.935 -1.401) 9.969 2.145 ( -2.376 -1.372 -0.288) 2.759 2.165 ( -2.319 -1.339 0.643) 2.753 2.170 ( -0.220 -0.127 -0.318) 0.407 2.365 ( -2.964 -1.711 1.659) 3.804 2.366 ( -1.924 -1.111 -1.692) 2.792 2.387 ( -1.142 -0.659 0.116) 1.324 2.609 ( 0.774 0.447 0.167) 0.909 2.611 ( 1.154 0.666 -0.445) 1.405 2.615 ( -0.009 -0.005 0.257) 0.257 2.812 ( -0.903 -0.521 -0.436) 1.130 2.826 ( -0.284 -0.164 -0.501) 0.598 2.852 ( 1.032 0.596 0.839) 1.458 3.723 ( -22.939 -13.244 -2.413) 26.598 3.763 ( -20.664 -11.930 4.025) 24.198 3.890 ( -26.180 -15.115 -2.066) 30.300 4.742 ( 3.132 1.808 7.236) 8.090 4.811 ( 3.100 1.790 -7.647) 8.443 5.043 ( 1.215 0.701 0.862) 1.646 ======================= Grid point 221 (30/48) ======================= q-point: ( 0.36 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.032 ( 9.305 5.372 0.149) 10.746 1.176 ( 8.851 5.110 1.535) 10.335 1.188 ( 8.535 4.928 -1.561) 9.978 1.213 ( 4.446 2.567 -0.299) 5.142 1.221 ( 5.382 3.107 1.130) 6.316 1.230 ( 4.373 2.525 -0.890) 5.127 1.742 ( 2.378 1.373 -1.409) 3.086 1.748 ( 2.070 1.195 2.460) 3.430 1.825 ( 3.712 2.143 -1.350) 4.494 2.138 ( 1.468 0.848 -1.101) 2.022 2.156 ( 1.390 0.802 1.247) 2.032 2.192 ( 1.837 1.061 -0.049) 2.122 2.313 ( -1.417 -0.818 1.725) 2.378 2.354 ( 0.391 0.226 -1.987) 2.038 2.379 ( -0.102 -0.059 0.581) 0.593 2.590 ( -2.076 -1.199 0.745) 2.511 2.596 ( -1.897 -1.095 -0.311) 2.212 2.600 ( -2.143 -1.237 -0.488) 2.522 2.814 ( 0.912 0.527 0.308) 1.098 2.835 ( 0.789 0.456 -0.830) 1.233 2.881 ( 1.203 0.695 0.418) 1.451 3.193 ( -22.433 -12.952 -1.424) 25.943 3.235 ( -23.564 -13.605 2.623) 27.336 3.301 ( -24.572 -14.186 -1.233) 28.400 4.795 ( 1.593 0.920 4.653) 5.004 4.892 ( 3.356 1.938 -6.722) 7.759 5.062 ( 0.231 0.133 2.133) 2.150 ======================= Grid point 223 (31/48) ======================= q-point: ( 0.43 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.228 ( 7.097 4.097 1.153) 8.276 1.313 ( 2.798 1.616 -0.749) 3.317 1.320 ( 4.445 2.566 -0.572) 5.164 1.325 ( 3.235 1.868 -2.276) 4.374 1.334 ( 4.111 2.373 1.294) 4.920 1.340 ( 5.067 2.925 1.202) 5.973 1.780 ( 0.792 0.457 1.266) 1.562 1.795 ( 2.303 1.330 -1.567) 3.087 1.865 ( 0.008 0.004 0.219) 0.219 2.189 ( 2.400 1.386 -0.797) 2.884 2.213 ( 3.181 1.837 1.248) 3.880 2.239 ( 2.036 1.176 -0.370) 2.380 2.306 ( 0.751 0.434 0.162) 0.883 2.366 ( -1.031 -0.595 1.104) 1.624 2.366 ( 0.397 0.229 -1.229) 1.312 2.526 ( -3.710 -2.142 -0.631) 4.330 2.530 ( -3.464 -2.000 -0.026) 4.000 2.533 ( -2.548 -1.471 0.585) 3.000 2.754 ( -13.631 -7.870 -1.417) 15.803 2.772 ( -15.176 -8.762 1.541) 17.592 2.828 ( -13.805 -7.971 0.095) 15.941 2.852 ( 1.050 0.606 0.849) 1.480 2.856 ( -0.943 -0.545 -0.831) 1.370 2.900 ( -0.039 -0.023 -0.208) 0.212 4.819 ( 0.622 0.359 2.357) 2.464 4.958 ( 2.597 1.499 -5.739) 6.475 5.050 ( -1.189 -0.687 3.381) 3.650 ======================= Grid point 225 (32/48) ======================= q-point: (-0.50 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.334 ( -0.001 -0.000 -0.001) 0.001 1.336 ( 0.001 0.000 0.001) 0.001 1.347 ( 0.000 0.000 0.000) 0.000 1.381 ( -0.000 -0.000 0.001) 0.001 1.382 ( 0.000 0.000 -0.001) 0.001 1.397 ( 0.000 0.000 0.000) 0.000 1.787 ( 0.000 0.000 -0.000) 0.000 1.841 ( -0.001 -0.000 0.000) 0.001 1.842 ( 0.001 0.000 -0.000) 0.001 2.221 ( 0.000 0.000 -0.000) 0.000 2.265 ( -0.490 -0.283 -1.182) 1.311 2.265 ( 0.490 0.283 1.182) 1.311 2.335 ( 0.000 0.000 -0.000) 0.000 2.335 ( -0.000 -0.000 0.000) 0.000 2.373 ( 0.000 0.000 0.000) 0.000 2.475 ( 0.000 0.000 -0.000) 0.000 2.475 ( -0.000 -0.000 0.000) 0.000 2.498 ( 0.000 0.000 -0.000) 0.000 2.589 ( 0.000 0.000 -0.000) 0.000 2.589 ( -0.000 -0.000 0.000) 0.000 2.649 ( 0.000 0.000 -0.000) 0.000 2.861 ( -0.000 -0.000 0.000) 0.000 2.885 ( 0.000 0.000 0.000) 0.000 2.886 ( -0.000 -0.000 -0.000) 0.000 4.826 ( 0.000 0.000 0.000) 0.000 5.012 ( 2.093 1.209 -4.606) 5.202 5.012 ( -2.093 -1.209 4.606) 5.202 ======================= Grid point 243 (33/48) ======================= q-point: ( 0.07 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.484 ( 3.907 6.767 0.000) 7.813 0.486 ( 3.816 6.610 -0.000) 7.632 0.645 ( 8.762 15.176 0.000) 17.523 0.647 ( 8.765 15.182 -0.000) 17.531 0.709 ( 1.680 2.910 0.000) 3.360 0.799 ( 7.128 12.347 0.000) 14.257 0.949 ( 10.682 18.503 0.000) 21.365 0.950 ( 10.680 18.498 -0.000) 21.360 1.017 ( 10.236 17.729 0.000) 20.471 2.375 ( -0.340 -0.588 0.000) 0.679 2.410 ( 0.021 0.037 -0.000) 0.042 2.410 ( 0.030 0.052 0.000) 0.060 2.463 ( -1.887 -3.269 0.000) 3.775 2.501 ( -1.492 -2.585 -0.000) 2.984 2.501 ( -1.491 -2.583 0.000) 2.982 2.537 ( -1.121 -1.941 0.000) 2.241 2.538 ( -1.125 -1.949 -0.000) 2.251 2.558 ( -0.469 -0.812 0.000) 0.938 2.897 ( -1.511 -2.617 0.000) 3.022 2.945 ( -1.896 -3.283 -0.000) 3.791 2.946 ( -1.896 -3.284 0.000) 3.793 3.498 ( 18.310 31.714 0.000) 36.620 4.254 ( 10.115 17.520 -0.003) 20.230 4.254 ( 10.110 17.510 0.003) 20.219 5.159 ( -8.756 -15.169 0.001) 17.515 5.160 ( -8.765 -15.179 -0.001) 17.528 5.230 ( -8.554 -14.815 0.000) 17.107 ======================= Grid point 245 (34/48) ======================= q-point: ( 0.14 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.638 ( 6.551 6.486 4.721) 10.357 0.650 ( 7.529 6.145 -4.549) 10.730 0.745 ( 1.385 0.821 -0.264) 1.632 0.907 ( 5.702 15.599 -1.862) 16.713 0.942 ( 6.504 12.819 2.017) 14.515 1.032 ( 3.828 12.879 -0.107) 13.437 1.278 ( 12.127 11.296 -0.570) 16.583 1.313 ( 13.297 8.941 0.630) 16.036 1.343 ( 12.126 8.151 -0.013) 14.610 2.311 ( -6.608 2.266 -0.905) 7.044 2.348 ( -6.885 0.495 1.665) 7.100 2.362 ( -6.703 1.673 -0.952) 6.974 2.445 ( 1.713 -3.979 1.196) 4.494 2.466 ( -0.805 -2.824 -0.261) 2.948 2.470 ( -1.827 -2.972 -1.338) 3.736 2.501 ( 2.011 -1.145 1.198) 2.606 2.534 ( 3.891 -2.210 -1.040) 4.594 2.536 ( 0.435 -0.863 0.415) 1.052 2.846 ( -1.093 -2.382 0.215) 2.630 2.872 ( -2.364 -3.312 0.886) 4.164 2.873 ( -2.371 -3.559 -1.054) 4.405 3.998 ( 10.061 15.950 1.144) 18.893 4.347 ( -9.576 1.973 -11.772) 15.303 4.510 ( 8.248 2.845 16.093) 18.306 4.744 ( -20.222 -13.182 -4.578) 24.570 4.897 ( -10.593 -16.061 -2.485) 19.399 4.953 ( 1.469 -14.114 1.570) 14.277 ======================= Grid point 247 (35/48) ======================= q-point: ( 0.21 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.797 ( 3.556 2.602 -0.150) 4.409 0.834 ( 9.289 5.596 1.343) 10.927 0.849 ( 8.983 5.255 -1.267) 10.484 1.124 ( 2.929 13.199 -0.805) 13.544 1.133 ( 1.715 12.737 1.107) 12.900 1.173 ( -0.571 10.858 -0.216) 10.875 1.561 ( 11.173 5.003 -1.271) 12.307 1.569 ( 9.383 3.615 1.705) 10.199 1.574 ( 8.836 2.271 -0.385) 9.131 2.235 ( -6.889 1.766 0.513) 7.131 2.239 ( -4.941 3.961 -0.835) 6.388 2.259 ( -6.835 2.324 0.247) 7.224 2.401 ( -2.929 -0.939 -1.539) 3.440 2.406 ( -3.016 -0.917 1.216) 3.379 2.410 ( -1.698 -0.837 0.402) 1.935 2.558 ( 4.302 -1.417 -0.209) 4.534 2.571 ( 3.290 -0.810 0.527) 3.429 2.581 ( 3.108 -1.783 -0.334) 3.599 2.809 ( -1.026 -3.000 -0.586) 3.225 2.817 ( -0.388 -2.583 -0.090) 2.614 2.825 ( 1.116 -1.743 0.654) 2.171 4.016 ( -22.724 -7.429 -3.227) 24.125 4.033 ( -17.204 -0.090 4.696) 17.833 4.218 ( -26.317 -11.444 -3.170) 28.872 4.639 ( 7.270 -3.257 9.157) 12.137 4.694 ( 5.191 -10.082 -7.222) 13.444 4.933 ( 6.968 -8.729 -0.247) 11.172 ======================= Grid point 249 (36/48) ======================= q-point: ( 0.29 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.933 ( 7.735 5.084 0.003) 9.256 1.051 ( 10.014 4.131 -0.783) 10.861 1.055 ( 9.476 4.330 0.943) 10.461 1.255 ( -1.141 10.278 -0.264) 10.345 1.268 ( 0.325 11.846 0.465) 11.859 1.269 ( 0.334 11.225 -0.286) 11.234 1.700 ( 4.587 1.155 -0.980) 4.830 1.710 ( 4.548 1.759 2.410) 5.439 1.754 ( 7.284 2.274 -1.577) 7.792 2.168 ( -2.899 2.404 -1.083) 3.919 2.188 ( -3.616 1.672 1.386) 4.219 2.220 ( -2.034 3.796 -0.262) 4.315 2.360 ( -2.724 2.122 1.530) 3.777 2.373 ( -1.079 2.128 -1.598) 2.871 2.394 ( -0.677 1.538 0.266) 1.701 2.592 ( 1.289 -1.993 -0.266) 2.388 2.594 ( 0.350 -1.623 0.230) 1.676 2.601 ( 0.673 -0.742 0.001) 1.002 2.787 ( 1.268 -1.881 0.072) 2.269 2.807 ( 1.570 -1.778 -0.821) 2.511 2.844 ( 2.569 -1.444 0.655) 3.019 3.508 ( -26.096 -8.289 -1.369) 27.415 3.571 ( -27.668 -7.461 3.325) 28.849 3.639 ( -28.288 -10.192 -2.114) 30.143 4.709 ( 5.869 -4.837 5.916) 9.636 4.765 ( 9.133 -6.151 -6.757) 12.919 4.969 ( 6.214 -7.757 1.000) 9.989 ======================= Grid point 251 (37/48) ======================= q-point: ( 0.36 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.130 ( 8.691 4.598 0.569) 9.849 1.238 ( 7.474 2.436 -1.615) 8.026 1.244 ( 7.922 2.821 1.362) 8.519 1.349 ( 0.060 10.138 -0.308) 10.143 1.365 ( -0.890 9.196 -0.426) 9.249 1.379 ( -0.258 9.804 0.476) 9.819 1.772 ( 2.367 1.804 -1.119) 3.179 1.773 ( 1.351 1.401 1.547) 2.486 1.852 ( 2.318 0.756 -0.655) 2.525 2.179 ( 1.715 0.475 1.353) 2.235 2.187 ( 1.765 3.141 -1.538) 3.917 2.238 ( 0.728 1.911 0.367) 2.078 2.358 ( -2.080 5.456 0.798) 5.894 2.389 ( -0.970 3.607 -1.361) 3.975 2.403 ( -1.343 3.008 0.648) 3.358 2.563 ( -2.195 -1.835 -0.303) 2.877 2.565 ( -1.952 -0.823 0.386) 2.153 2.576 ( -2.613 -0.466 -0.133) 2.657 2.807 ( 2.275 -1.104 0.737) 2.634 2.825 ( 2.059 -1.365 -0.845) 2.611 2.871 ( 2.167 -1.666 0.084) 2.734 3.013 ( -22.647 -4.852 -1.193) 23.192 3.039 ( -24.967 -5.684 1.872) 25.675 3.092 ( -24.632 -6.603 -0.676) 25.511 4.748 ( 4.706 -5.523 3.510) 8.061 4.850 ( 7.687 -5.801 -5.886) 11.286 4.982 ( 4.280 -7.862 2.398) 9.267 ======================= Grid point 253 (38/48) ======================= q-point: ( 0.43 0.07 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.287 ( 5.168 1.958 1.601) 5.753 1.331 ( 1.723 -0.476 -2.194) 2.830 1.363 ( 3.716 0.630 0.650) 3.825 1.413 ( -2.812 6.766 -0.308) 7.334 1.433 ( -1.442 8.020 -0.660) 8.175 1.443 ( -2.322 7.370 0.941) 7.784 1.790 ( -0.057 0.515 0.501) 0.721 1.823 ( 1.660 1.363 -1.414) 2.572 1.857 ( -0.981 -0.741 0.906) 1.527 2.252 ( 5.496 0.060 -0.204) 5.500 2.279 ( 4.093 4.276 -0.614) 5.951 2.280 ( 3.133 -0.752 0.841) 3.330 2.379 ( -3.926 7.564 0.627) 8.545 2.382 ( -5.244 3.577 -0.630) 6.379 2.394 ( -5.124 5.865 0.038) 7.788 2.507 ( -2.298 0.601 -0.672) 2.469 2.520 ( -3.391 0.987 0.451) 3.560 2.531 ( -1.709 1.384 0.143) 2.203 2.671 ( -11.335 0.360 -1.374) 11.424 2.675 ( -12.265 -0.063 1.247) 12.329 2.728 ( -11.677 -1.035 0.128) 11.724 2.843 ( 1.269 -1.221 -0.316) 1.789 2.855 ( 2.067 -0.412 1.039) 2.350 2.878 ( 0.508 -1.863 -0.703) 2.055 4.764 ( 3.832 -5.862 1.189) 7.104 4.915 ( 6.614 -5.778 -4.899) 10.057 4.961 ( 2.599 -8.156 3.707) 9.329 ======================= Grid point 273 (39/48) ======================= q-point: ( 0.14 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.779 ( 1.848 3.201 0.000) 3.696 0.786 ( 4.685 8.115 0.000) 9.370 0.787 ( 4.617 7.996 -0.000) 9.233 1.199 ( 7.309 12.660 0.000) 14.618 1.200 ( 7.294 12.633 -0.000) 14.588 1.273 ( 5.922 10.257 0.000) 11.844 1.472 ( 3.312 5.737 0.000) 6.624 1.473 ( 4.787 8.291 -0.000) 9.574 1.474 ( 4.795 8.306 0.000) 9.590 2.346 ( -1.965 -3.404 -0.000) 3.930 2.346 ( -1.967 -3.407 0.000) 3.934 2.350 ( -0.507 -0.877 0.000) 1.013 2.397 ( 0.446 0.773 0.000) 0.892 2.443 ( 0.197 0.341 -0.000) 0.394 2.443 ( 0.194 0.335 0.000) 0.387 2.490 ( 0.130 0.225 0.000) 0.259 2.491 ( 0.138 0.239 -0.000) 0.275 2.523 ( 0.292 0.506 0.000) 0.584 2.801 ( -1.950 -3.331 -0.042) 3.860 2.802 ( -1.165 -2.065 0.042) 2.371 2.802 ( -1.802 -3.122 0.000) 3.605 4.252 ( 4.754 8.234 0.000) 9.508 4.313 ( -18.408 -7.982 -1.220) 20.101 4.313 ( 2.291 -19.932 1.220) 20.101 4.528 ( -11.353 -19.664 -0.000) 22.706 4.664 ( -2.736 -4.744 -0.001) 5.476 4.664 ( -2.737 -4.737 0.001) 5.471 ======================= Grid point 275 (40/48) ======================= q-point: ( 0.21 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.878 ( 4.839 5.594 0.035) 7.397 0.956 ( 6.773 6.666 -0.477) 9.515 0.961 ( 6.956 6.182 0.480) 9.318 1.396 ( 0.844 13.470 -0.264) 13.499 1.399 ( 0.854 13.253 0.469) 13.289 1.409 ( -2.134 12.301 -0.164) 12.485 1.602 ( 7.729 1.268 -0.271) 7.837 1.630 ( 6.629 3.044 2.107) 7.593 1.639 ( 8.128 3.285 -1.891) 8.969 2.252 ( -4.635 -0.923 -1.135) 4.860 2.266 ( -4.082 -2.221 1.520) 4.889 2.306 ( -3.656 1.138 -0.530) 3.866 2.422 ( 0.022 2.281 0.215) 2.292 2.429 ( -2.168 3.594 -0.556) 4.234 2.442 ( -2.971 4.078 0.533) 5.074 2.523 ( 4.068 -1.500 0.449) 4.359 2.543 ( 3.627 -1.605 -0.728) 4.033 2.563 ( 2.750 0.513 0.265) 2.810 2.750 ( -0.193 -2.594 0.195) 2.608 2.759 ( 0.637 -2.930 -0.557) 3.050 2.785 ( 1.330 -2.048 0.330) 2.464 3.872 ( -23.217 -7.543 0.004) 24.411 3.978 ( -19.106 -13.205 -3.709) 23.520 3.998 ( -27.074 -4.412 4.176) 27.747 4.472 ( 13.849 -10.868 -5.549) 18.458 4.561 ( 1.490 -3.863 6.921) 8.065 4.719 ( 9.921 -11.971 -1.846) 15.657 ======================= Grid point 277 (41/48) ======================= q-point: ( 0.29 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.042 ( 7.194 5.621 0.254) 9.133 1.136 ( 7.534 4.244 -1.738) 8.820 1.138 ( 7.582 4.044 1.603) 8.741 1.487 ( -2.295 12.234 -0.190) 12.449 1.504 ( -2.288 11.768 -0.059) 11.988 1.512 ( -1.726 11.989 0.191) 12.114 1.733 ( 5.028 2.349 -0.725) 5.597 1.750 ( 3.901 2.201 1.826) 4.837 1.791 ( 6.115 1.551 -1.266) 6.434 2.183 ( -2.630 -1.805 1.895) 3.711 2.193 ( -1.078 -0.293 -2.115) 2.392 2.256 ( -2.253 -0.739 0.311) 2.391 2.446 ( -1.760 5.056 0.484) 5.375 2.446 ( -0.770 4.488 -0.548) 4.587 2.448 ( -0.250 3.187 0.160) 3.200 2.556 ( 1.120 -1.314 0.216) 1.741 2.569 ( 0.807 -0.707 -0.394) 1.143 2.595 ( 0.218 0.194 0.135) 0.322 2.749 ( 2.335 -1.619 0.747) 2.939 2.763 ( 2.572 -2.285 -0.944) 3.568 2.801 ( 2.924 -2.555 0.193) 3.888 3.406 ( -24.880 -1.868 -0.291) 24.952 3.483 ( -28.475 -0.916 2.458) 28.596 3.497 ( -25.682 -4.014 -2.204) 26.087 4.571 ( 6.173 -8.682 4.670) 11.631 4.590 ( 12.332 -10.942 -5.023) 17.235 4.759 ( 9.024 -12.957 0.384) 15.794 ======================= Grid point 279 (42/48) ======================= q-point: ( 0.36 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.208 ( 6.514 3.145 0.931) 7.294 1.274 ( 4.821 1.097 -2.214) 5.417 1.284 ( 5.900 1.207 1.334) 6.168 1.558 ( -3.655 9.787 -0.253) 10.450 1.567 ( -2.321 11.004 -0.364) 11.252 1.581 ( -2.936 10.169 0.585) 10.601 1.799 ( 0.166 1.284 0.725) 1.484 1.813 ( 1.719 2.150 -0.995) 2.927 1.856 ( 0.759 -0.237 0.205) 0.821 2.168 ( 2.810 -1.417 0.342) 3.165 2.212 ( 3.428 -0.889 -1.350) 3.790 2.223 ( 0.676 -2.886 1.092) 3.159 2.471 ( -0.346 2.577 0.269) 2.614 2.472 ( -0.791 4.141 -0.102) 4.217 2.480 ( -0.477 3.997 -0.240) 4.033 2.542 ( -1.232 -0.402 -0.115) 1.301 2.566 ( -0.786 0.516 -0.347) 1.002 2.575 ( -2.121 -0.526 0.470) 2.235 2.782 ( 2.407 -1.100 0.988) 2.825 2.789 ( 2.385 -1.878 -0.623) 3.099 2.819 ( 1.793 -2.928 -0.318) 3.448 3.033 ( -20.542 7.749 -0.231) 21.956 3.066 ( -23.402 9.041 0.987) 25.107 3.078 ( -22.984 6.203 -0.768) 23.819 4.589 ( 6.759 -10.299 2.369) 12.545 4.678 ( 10.884 -11.274 -4.488) 16.300 4.767 ( 7.222 -13.566 2.131) 15.516 ======================= Grid point 281 (43/48) ======================= q-point: ( 0.43 0.14 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.302 ( 0.731 -1.266 -0.000) 1.462 1.303 ( 0.744 -1.288 0.000) 1.488 1.344 ( 1.246 -2.158 0.000) 2.492 1.575 ( -4.981 8.628 -0.000) 9.962 1.606 ( -5.024 8.701 0.001) 10.047 1.607 ( -5.020 8.696 -0.001) 10.041 1.802 ( -0.489 0.847 0.000) 0.979 1.843 ( -0.082 0.141 0.000) 0.164 1.844 ( -0.076 0.132 -0.000) 0.152 2.210 ( 2.159 -3.751 0.004) 4.328 2.210 ( 2.162 -3.733 -0.004) 4.314 2.250 ( 2.204 -3.817 0.000) 4.407 2.487 ( -0.948 1.645 0.002) 1.898 2.487 ( -0.951 1.643 -0.002) 1.899 2.498 ( -1.527 2.646 0.000) 3.055 2.532 ( -0.422 0.730 0.000) 0.843 2.533 ( -0.425 0.737 -0.000) 0.851 2.560 ( -0.567 0.982 0.000) 1.133 2.796 ( 0.369 -0.639 0.000) 0.738 2.797 ( -1.191 2.062 -0.000) 2.382 2.797 ( -1.203 2.084 0.000) 2.407 2.902 ( -9.866 17.089 0.000) 19.732 2.909 ( -8.749 15.151 -0.000) 17.495 2.909 ( -8.739 15.139 0.000) 17.480 4.596 ( 6.306 -10.922 0.000) 12.612 4.739 ( 9.519 -11.842 -3.520) 15.595 4.739 ( 5.496 -14.164 3.520) 15.595 ======================= Grid point 305 (44/48) ======================= q-point: ( 0.21 0.21 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.007 ( 4.107 7.113 0.000) 8.213 1.089 ( 3.854 6.675 -0.000) 7.708 1.089 ( 3.835 6.643 0.000) 7.671 1.622 ( 4.839 8.409 -0.477) 9.714 1.622 ( 4.863 8.395 0.477) 9.714 1.632 ( 4.798 8.309 0.000) 9.595 1.641 ( 2.574 4.458 0.000) 5.148 1.701 ( 2.453 4.248 -0.000) 4.905 1.702 ( 2.459 4.259 0.000) 4.918 2.204 ( -1.938 -3.357 0.000) 3.876 2.205 ( -1.935 -3.352 -0.000) 3.870 2.279 ( -1.840 -3.186 0.000) 3.679 2.468 ( 0.984 1.705 0.000) 1.969 2.498 ( 2.287 0.765 -0.708) 2.513 2.498 ( -0.481 2.363 0.708) 2.513 2.529 ( 0.809 1.401 0.000) 1.618 2.530 ( 0.797 1.380 -0.000) 1.593 2.589 ( 0.760 1.316 0.000) 1.519 2.720 ( -0.264 -0.457 0.000) 0.528 2.720 ( -0.259 -0.449 -0.000) 0.518 2.756 ( -0.381 -0.660 0.000) 0.762 3.700 ( -7.137 -12.357 0.001) 14.270 3.701 ( -7.133 -12.360 -0.001) 14.271 3.817 ( -7.818 -13.541 0.000) 15.636 4.342 ( -1.052 -1.822 0.000) 2.104 4.505 ( 3.050 -8.753 -3.543) 9.923 4.505 ( -9.105 -1.735 3.543) 9.923 ======================= Grid point 307 (45/48) ======================= q-point: ( 0.29 0.21 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.145 ( 4.924 4.824 0.304) 6.900 1.212 ( 4.319 3.395 -1.584) 5.717 1.212 ( 4.443 3.429 1.290) 5.758 1.699 ( -1.147 8.211 -0.246) 8.294 1.706 ( -2.770 8.244 -0.049) 8.697 1.717 ( -1.626 8.316 0.289) 8.478 1.792 ( 3.130 2.044 0.694) 3.803 1.794 ( 4.229 3.871 -0.220) 5.738 1.813 ( 4.728 0.611 -0.537) 4.797 2.141 ( -1.370 -1.877 1.290) 2.658 2.162 ( 0.617 -2.380 -1.697) 2.987 2.205 ( -2.044 -3.714 0.481) 4.266 2.489 ( 0.061 1.010 0.066) 1.014 2.511 ( -0.384 1.697 0.389) 1.783 2.517 ( -0.814 1.787 -0.772) 2.110 2.542 ( 0.147 -0.337 0.765) 0.849 2.562 ( 1.087 -0.473 -0.741) 1.398 2.587 ( -0.930 -1.539 0.264) 1.818 2.734 ( 1.854 0.189 0.919) 2.078 2.736 ( 1.924 -0.182 -0.725) 2.064 2.760 ( 1.702 -1.236 -0.171) 2.111 3.439 ( -15.546 4.971 1.208) 16.366 3.484 ( -13.864 2.509 -2.534) 14.315 3.538 ( -18.359 5.857 1.368) 19.320 4.361 ( 9.743 -11.597 -2.677) 15.381 4.389 ( 1.508 -9.322 3.212) 9.975 4.484 ( 7.397 -14.350 -0.578) 16.155 ======================= Grid point 309 (46/48) ======================= q-point: ( 0.36 0.21 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.254 ( 3.570 1.605 1.032) 4.048 1.282 ( 1.509 -0.140 -1.852) 2.393 1.294 ( 2.725 -0.031 0.807) 2.842 1.725 ( -3.347 6.684 -0.066) 7.475 1.749 ( -1.944 6.839 -0.748) 7.149 1.755 ( -2.706 6.913 0.815) 7.468 1.824 ( -0.586 1.355 -0.007) 1.476 1.848 ( 0.087 0.104 0.122) 0.183 1.849 ( -0.210 0.996 -0.102) 1.023 2.138 ( 2.113 -1.162 -0.236) 2.423 2.160 ( 0.389 -2.786 0.789) 2.921 2.175 ( 2.675 -2.249 -0.529) 3.534 2.494 ( 0.039 0.274 0.195) 0.339 2.508 ( -1.303 0.036 -0.323) 1.343 2.511 ( -0.271 -0.416 0.211) 0.539 2.532 ( 0.058 -0.677 -0.754) 1.015 2.546 ( -1.050 -1.829 0.751) 2.239 2.565 ( 0.299 -0.678 -0.122) 0.751 2.765 ( 1.429 -0.521 -0.257) 1.543 2.767 ( 1.720 -0.434 0.810) 1.950 2.776 ( 1.260 -1.392 -0.531) 1.951 3.304 ( -13.095 18.168 0.607) 22.404 3.326 ( -15.014 17.203 -1.121) 22.861 3.346 ( -17.166 17.354 0.504) 24.415 4.353 ( 6.929 -13.383 1.213) 15.119 4.419 ( 11.791 -14.830 -2.753) 19.145 4.463 ( 8.014 -17.072 1.553) 18.924 ======================= Grid point 335 (47/48) ======================= q-point: ( 0.29 0.29 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.229 ( 2.027 3.511 0.000) 4.055 1.267 ( 1.122 1.944 -0.000) 2.244 1.267 ( 1.126 1.951 0.000) 2.252 1.796 ( 1.484 2.570 0.000) 2.968 1.814 ( 0.774 1.340 0.000) 1.547 1.815 ( 0.776 1.346 -0.000) 1.553 1.829 ( 0.723 2.756 -0.430) 2.882 1.829 ( 2.026 2.004 0.430) 2.882 1.866 ( 1.198 2.075 0.000) 2.397 2.126 ( 0.935 -0.774 -0.539) 1.328 2.126 ( -1.137 0.423 0.539) 1.328 2.144 ( -1.072 -1.856 0.000) 2.143 2.497 ( -0.024 -0.041 0.000) 0.048 2.518 ( -0.655 -1.134 -0.000) 1.309 2.518 ( -0.653 -1.131 0.000) 1.306 2.539 ( -1.835 0.297 -0.210) 1.871 2.539 ( 1.174 -1.441 0.210) 1.871 2.542 ( -1.299 -2.250 0.000) 2.598 2.749 ( 0.577 1.046 -0.517) 1.302 2.749 ( 0.617 1.023 0.517) 1.302 2.755 ( 0.255 0.442 0.000) 0.511 3.577 ( 5.714 6.193 -2.127) 8.690 3.577 ( 2.506 8.044 2.127) 8.690 3.665 ( 3.117 5.399 0.000) 6.234 4.162 ( -4.417 -7.650 0.000) 8.833 4.219 ( -5.628 -9.752 -0.001) 11.260 4.219 ( -5.630 -9.748 0.001) 11.257 ======================= Grid point 337 (48/48) ======================= q-point: ( 0.36 0.29 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 5.87e-04 5.87e-04 5.87e-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.275 ( -0.118 0.205 0.000) 0.236 1.276 ( -0.121 0.210 -0.000) 0.242 1.288 ( 0.305 -0.528 0.000) 0.609 1.813 ( -1.363 2.360 0.000) 2.725 1.826 ( -0.513 0.889 -0.000) 1.027 1.826 ( -0.518 0.897 0.000) 1.036 1.854 ( -0.837 1.450 0.000) 1.675 1.866 ( -0.574 0.997 -0.000) 1.151 1.866 ( -0.575 0.993 0.000) 1.148 2.127 ( 0.201 -0.349 0.000) 0.403 2.127 ( 0.201 -0.347 -0.000) 0.400 2.137 ( 0.836 -1.449 0.000) 1.673 2.497 ( 0.229 -0.397 -0.000) 0.458 2.497 ( 0.164 -0.280 -0.000) 0.325 2.497 ( 0.155 -0.272 0.000) 0.313 2.523 ( 0.076 -0.131 -0.000) 0.152 2.523 ( 0.081 -0.140 0.000) 0.162 2.543 ( 0.826 -1.431 -0.000) 1.652 2.761 ( 0.024 -0.042 0.000) 0.048 2.762 ( 0.093 -0.161 -0.000) 0.186 2.762 ( 0.096 -0.166 0.000) 0.191 3.684 ( -9.884 17.120 0.000) 19.769 3.688 ( -10.180 17.610 -0.001) 20.341 3.688 ( -10.160 17.620 0.001) 20.340 4.077 ( 7.871 -13.634 0.000) 15.743 4.108 ( 7.655 -18.310 0.708) 19.858 4.108 ( 12.030 -15.784 -0.708) 19.858 =================== 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 73.940 73.940 6.714 -0.000 0.000 -0.000 3/10584 20.0 32.700 32.700 2.111 -0.000 0.000 -0.000 3/10584 30.0 19.483 19.483 1.162 -0.000 0.000 -0.000 3/10584 40.0 14.191 14.191 0.858 -0.000 0.000 -0.000 3/10584 50.0 11.412 11.412 0.721 -0.000 0.000 -0.000 3/10584 60.0 9.650 9.650 0.638 -0.000 0.000 -0.000 3/10584 70.0 8.399 8.399 0.577 -0.000 0.000 -0.000 3/10584 80.0 7.449 7.449 0.528 -0.000 0.000 -0.000 3/10584 90.0 6.696 6.696 0.486 -0.000 0.000 -0.000 3/10584 100.0 6.082 6.082 0.449 -0.000 0.000 -0.000 3/10584 110.0 5.572 5.572 0.417 -0.000 0.000 -0.000 3/10584 120.0 5.139 5.139 0.389 -0.000 0.000 -0.000 3/10584 130.0 4.768 4.768 0.365 -0.000 0.000 -0.000 3/10584 140.0 4.447 4.447 0.343 -0.000 0.000 -0.000 3/10584 150.0 4.166 4.166 0.323 -0.000 0.000 -0.000 3/10584 160.0 3.918 3.918 0.305 -0.000 0.000 -0.000 3/10584 170.0 3.697 3.697 0.289 -0.000 0.000 -0.000 3/10584 180.0 3.500 3.500 0.275 -0.000 0.000 -0.000 3/10584 190.0 3.322 3.322 0.262 -0.000 0.000 -0.000 3/10584 200.0 3.162 3.162 0.250 -0.000 0.000 -0.000 3/10584 210.0 3.016 3.016 0.239 -0.000 0.000 -0.000 3/10584 220.0 2.883 2.883 0.229 -0.000 0.000 -0.000 3/10584 230.0 2.761 2.761 0.220 -0.000 0.000 -0.000 3/10584 240.0 2.648 2.648 0.211 -0.000 0.000 -0.000 3/10584 250.0 2.545 2.545 0.203 -0.000 0.000 -0.000 3/10584 260.0 2.449 2.449 0.196 -0.000 0.000 -0.000 3/10584 270.0 2.361 2.361 0.189 -0.000 0.000 -0.000 3/10584 280.0 2.278 2.278 0.182 -0.000 0.000 -0.000 3/10584 290.0 2.201 2.201 0.176 -0.000 0.000 -0.000 3/10584 300.0 2.129 2.129 0.171 -0.000 0.000 -0.000 3/10584 310.0 2.061 2.061 0.165 -0.000 0.000 -0.000 3/10584 320.0 1.998 1.998 0.161 -0.000 0.000 -0.000 3/10584 330.0 1.938 1.938 0.156 -0.000 0.000 -0.000 3/10584 340.0 1.882 1.882 0.151 -0.000 0.000 -0.000 3/10584 350.0 1.829 1.829 0.147 -0.000 0.000 -0.000 3/10584 360.0 1.779 1.779 0.143 -0.000 0.000 -0.000 3/10584 370.0 1.732 1.732 0.139 -0.000 0.000 -0.000 3/10584 380.0 1.687 1.687 0.136 -0.000 0.000 -0.000 3/10584 390.0 1.644 1.644 0.133 -0.000 0.000 -0.000 3/10584 400.0 1.603 1.603 0.129 -0.000 0.000 -0.000 3/10584 410.0 1.564 1.564 0.126 -0.000 0.000 -0.000 3/10584 420.0 1.528 1.528 0.123 -0.000 0.000 -0.000 3/10584 430.0 1.492 1.492 0.120 -0.000 0.000 -0.000 3/10584 440.0 1.459 1.459 0.118 -0.000 0.000 -0.000 3/10584 450.0 1.427 1.427 0.115 -0.000 0.000 -0.000 3/10584 460.0 1.396 1.396 0.113 -0.000 0.000 -0.000 3/10584 470.0 1.367 1.367 0.110 -0.000 0.000 -0.000 3/10584 480.0 1.338 1.338 0.108 -0.000 0.000 -0.000 3/10584 490.0 1.311 1.311 0.106 -0.000 0.000 -0.000 3/10584 500.0 1.285 1.285 0.104 -0.000 0.000 -0.000 3/10584 510.0 1.260 1.260 0.102 -0.000 0.000 -0.000 3/10584 520.0 1.236 1.236 0.100 -0.000 0.000 -0.000 3/10584 530.0 1.213 1.213 0.098 -0.000 0.000 -0.000 3/10584 540.0 1.191 1.191 0.096 -0.000 0.000 -0.000 3/10584 550.0 1.169 1.169 0.095 -0.000 0.000 -0.000 3/10584 560.0 1.149 1.149 0.093 -0.000 0.000 -0.000 3/10584 570.0 1.129 1.129 0.091 -0.000 0.000 -0.000 3/10584 580.0 1.109 1.109 0.090 -0.000 0.000 -0.000 3/10584 590.0 1.091 1.091 0.088 -0.000 0.000 -0.000 3/10584 600.0 1.072 1.072 0.087 -0.000 0.000 -0.000 3/10584 610.0 1.055 1.055 0.085 -0.000 0.000 -0.000 3/10584 620.0 1.038 1.038 0.084 -0.000 0.000 -0.000 3/10584 630.0 1.022 1.022 0.083 -0.000 0.000 -0.000 3/10584 640.0 1.006 1.006 0.081 -0.000 0.000 -0.000 3/10584 650.0 0.990 0.990 0.080 -0.000 0.000 -0.000 3/10584 660.0 0.976 0.976 0.079 -0.000 0.000 -0.000 3/10584 670.0 0.961 0.961 0.078 -0.000 0.000 -0.000 3/10584 680.0 0.947 0.947 0.077 -0.000 0.000 -0.000 3/10584 690.0 0.933 0.933 0.076 -0.000 0.000 -0.000 3/10584 700.0 0.920 0.920 0.075 -0.000 0.000 -0.000 3/10584 710.0 0.907 0.907 0.074 -0.000 0.000 -0.000 3/10584 720.0 0.895 0.895 0.073 -0.000 0.000 -0.000 3/10584 730.0 0.882 0.882 0.072 -0.000 0.000 -0.000 3/10584 740.0 0.871 0.871 0.071 -0.000 0.000 -0.000 3/10584 750.0 0.859 0.859 0.070 -0.000 0.000 -0.000 3/10584 760.0 0.848 0.848 0.069 -0.000 0.000 -0.000 3/10584 770.0 0.837 0.837 0.068 -0.000 0.000 -0.000 3/10584 780.0 0.826 0.826 0.067 -0.000 0.000 -0.000 3/10584 790.0 0.816 0.816 0.066 -0.000 0.000 -0.000 3/10584 800.0 0.806 0.806 0.065 -0.000 0.000 -0.000 3/10584 810.0 0.796 0.796 0.065 -0.000 0.000 -0.000 3/10584 820.0 0.786 0.786 0.064 -0.000 0.000 -0.000 3/10584 830.0 0.777 0.777 0.063 -0.000 0.000 -0.000 3/10584 840.0 0.767 0.767 0.062 -0.000 0.000 -0.000 3/10584 850.0 0.758 0.758 0.062 -0.000 0.000 -0.000 3/10584 860.0 0.750 0.750 0.061 -0.000 0.000 -0.000 3/10584 870.0 0.741 0.741 0.060 -0.000 0.000 -0.000 3/10584 880.0 0.733 0.733 0.059 -0.000 0.000 -0.000 3/10584 890.0 0.724 0.724 0.059 -0.000 0.000 -0.000 3/10584 900.0 0.716 0.716 0.058 -0.000 0.000 -0.000 3/10584 910.0 0.709 0.709 0.057 -0.000 0.000 -0.000 3/10584 920.0 0.701 0.701 0.057 -0.000 0.000 -0.000 3/10584 930.0 0.693 0.693 0.056 -0.000 0.000 -0.000 3/10584 940.0 0.686 0.686 0.056 -0.000 0.000 -0.000 3/10584 950.0 0.679 0.679 0.055 -0.000 0.000 -0.000 3/10584 960.0 0.672 0.672 0.055 -0.000 0.000 -0.000 3/10584 970.0 0.665 0.665 0.054 -0.000 0.000 -0.000 3/10584 980.0 0.658 0.658 0.053 -0.000 0.000 -0.000 3/10584 990.0 0.651 0.651 0.053 -0.000 0.000 -0.000 3/10584 1000.0 0.645 0.645 0.052 -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 06:39:15]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|