----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 18:55: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 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) Number of symmetry operations in supercell: 576 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.382603147621698 0.000000000000000 0.000000000000000 b -2.191301573810849 3.795445660546032 0.000000000000000 c 0.000000000000000 0.000000000000000 7.059576650000000 Atomic positions (fractional): *1 Yb 0.00000000000000 0.00000000000000 0.00000000000000 173.054 *2 I 0.66666666666667 0.33333333333333 0.75210898676006 126.904 3 I 0.33333333333333 0.66666666666667 0.24789101323994 126.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.382603147621698 0.000000000000000 0.000000000000000 b -2.191301573810849 3.795445660546032 0.000000000000000 c 0.000000000000000 0.000000000000000 7.059576650000000 Atomic positions (fractional): *1 Yb 0.00000000000000 0.00000000000000 0.00000000000000 173.054 > 1 *2 I 0.66666666666667 0.33333333333333 0.75210898676006 126.904 > 2 3 I 0.33333333333333 0.66666666666667 0.24789101323994 126.904 > 3 -------------------------------- super cell -------------------------------- Lattice vectors: a 17.530412590486790 0.000000000000000 0.000000000000000 b -8.765206295243395 15.181782642184128 0.000000000000000 c 0.000000000000000 0.000000000000000 21.178729950000001 Atomic positions (fractional): *1 Yb 0.00000000000000 0.00000000000000 0.00000000000000 173.054 > 1 2 Yb 0.25000000000000 0.00000000000000 0.00000000000000 173.054 > 1 3 Yb 0.50000000000000 0.00000000000000 0.00000000000000 173.054 > 1 4 Yb 0.75000000000000 0.00000000000000 0.00000000000000 173.054 > 1 5 Yb 0.00000000000000 0.25000000000000 0.00000000000000 173.054 > 1 6 Yb 0.25000000000000 0.25000000000000 0.00000000000000 173.054 > 1 7 Yb 0.50000000000000 0.25000000000000 0.00000000000000 173.054 > 1 8 Yb 0.75000000000000 0.25000000000000 0.00000000000000 173.054 > 1 9 Yb 0.00000000000000 0.50000000000000 0.00000000000000 173.054 > 1 10 Yb 0.25000000000000 0.50000000000000 0.00000000000000 173.054 > 1 11 Yb 0.50000000000000 0.50000000000000 0.00000000000000 173.054 > 1 12 Yb 0.75000000000000 0.50000000000000 0.00000000000000 173.054 > 1 13 Yb 0.00000000000000 0.75000000000000 0.00000000000000 173.054 > 1 14 Yb 0.25000000000000 0.75000000000000 0.00000000000000 173.054 > 1 15 Yb 0.50000000000000 0.75000000000000 0.00000000000000 173.054 > 1 16 Yb 0.75000000000000 0.75000000000000 0.00000000000000 173.054 > 1 17 Yb 0.00000000000000 0.00000000000000 0.33333333333333 173.054 > 1 18 Yb 0.25000000000000 0.00000000000000 0.33333333333333 173.054 > 1 19 Yb 0.50000000000000 0.00000000000000 0.33333333333333 173.054 > 1 20 Yb 0.75000000000000 0.00000000000000 0.33333333333333 173.054 > 1 21 Yb 0.00000000000000 0.25000000000000 0.33333333333333 173.054 > 1 22 Yb 0.25000000000000 0.25000000000000 0.33333333333333 173.054 > 1 23 Yb 0.50000000000000 0.25000000000000 0.33333333333333 173.054 > 1 24 Yb 0.75000000000000 0.25000000000000 0.33333333333333 173.054 > 1 25 Yb 0.00000000000000 0.50000000000000 0.33333333333333 173.054 > 1 26 Yb 0.25000000000000 0.50000000000000 0.33333333333333 173.054 > 1 27 Yb 0.50000000000000 0.50000000000000 0.33333333333333 173.054 > 1 28 Yb 0.75000000000000 0.50000000000000 0.33333333333333 173.054 > 1 29 Yb 0.00000000000000 0.75000000000000 0.33333333333333 173.054 > 1 30 Yb 0.25000000000000 0.75000000000000 0.33333333333333 173.054 > 1 31 Yb 0.50000000000000 0.75000000000000 0.33333333333333 173.054 > 1 32 Yb 0.75000000000000 0.75000000000000 0.33333333333333 173.054 > 1 33 Yb 0.00000000000000 0.00000000000000 0.66666666666667 173.054 > 1 34 Yb 0.25000000000000 0.00000000000000 0.66666666666667 173.054 > 1 35 Yb 0.50000000000000 0.00000000000000 0.66666666666667 173.054 > 1 36 Yb 0.75000000000000 0.00000000000000 0.66666666666667 173.054 > 1 37 Yb 0.00000000000000 0.25000000000000 0.66666666666667 173.054 > 1 38 Yb 0.25000000000000 0.25000000000000 0.66666666666667 173.054 > 1 39 Yb 0.50000000000000 0.25000000000000 0.66666666666667 173.054 > 1 40 Yb 0.75000000000000 0.25000000000000 0.66666666666667 173.054 > 1 41 Yb 0.00000000000000 0.50000000000000 0.66666666666667 173.054 > 1 42 Yb 0.25000000000000 0.50000000000000 0.66666666666667 173.054 > 1 43 Yb 0.50000000000000 0.50000000000000 0.66666666666667 173.054 > 1 44 Yb 0.75000000000000 0.50000000000000 0.66666666666667 173.054 > 1 45 Yb 0.00000000000000 0.75000000000000 0.66666666666667 173.054 > 1 46 Yb 0.25000000000000 0.75000000000000 0.66666666666667 173.054 > 1 47 Yb 0.50000000000000 0.75000000000000 0.66666666666667 173.054 > 1 48 Yb 0.75000000000000 0.75000000000000 0.66666666666667 173.054 > 1 *49 I 0.16666666666667 0.08333333333333 0.25070299558669 126.904 > 2 50 I 0.41666666666667 0.08333333333333 0.25070299558669 126.904 > 2 51 I 0.66666666666667 0.08333333333333 0.25070299558669 126.904 > 2 52 I 0.91666666666667 0.08333333333333 0.25070299558669 126.904 > 2 53 I 0.16666666666667 0.33333333333333 0.25070299558669 126.904 > 2 54 I 0.41666666666667 0.33333333333333 0.25070299558669 126.904 > 2 55 I 0.66666666666667 0.33333333333333 0.25070299558669 126.904 > 2 56 I 0.91666666666667 0.33333333333333 0.25070299558669 126.904 > 2 57 I 0.16666666666667 0.58333333333333 0.25070299558669 126.904 > 2 58 I 0.41666666666667 0.58333333333333 0.25070299558669 126.904 > 2 59 I 0.66666666666667 0.58333333333333 0.25070299558669 126.904 > 2 60 I 0.91666666666667 0.58333333333333 0.25070299558669 126.904 > 2 61 I 0.16666666666667 0.83333333333333 0.25070299558669 126.904 > 2 62 I 0.41666666666667 0.83333333333333 0.25070299558669 126.904 > 2 63 I 0.66666666666667 0.83333333333333 0.25070299558669 126.904 > 2 64 I 0.91666666666667 0.83333333333333 0.25070299558669 126.904 > 2 65 I 0.16666666666667 0.08333333333333 0.58403632892002 126.904 > 2 66 I 0.41666666666667 0.08333333333333 0.58403632892002 126.904 > 2 67 I 0.66666666666667 0.08333333333333 0.58403632892002 126.904 > 2 68 I 0.91666666666667 0.08333333333333 0.58403632892002 126.904 > 2 69 I 0.16666666666667 0.33333333333333 0.58403632892002 126.904 > 2 70 I 0.41666666666667 0.33333333333333 0.58403632892002 126.904 > 2 71 I 0.66666666666667 0.33333333333333 0.58403632892002 126.904 > 2 72 I 0.91666666666667 0.33333333333333 0.58403632892002 126.904 > 2 73 I 0.16666666666667 0.58333333333333 0.58403632892002 126.904 > 2 74 I 0.41666666666667 0.58333333333333 0.58403632892002 126.904 > 2 75 I 0.66666666666667 0.58333333333333 0.58403632892002 126.904 > 2 76 I 0.91666666666667 0.58333333333333 0.58403632892002 126.904 > 2 77 I 0.16666666666667 0.83333333333333 0.58403632892002 126.904 > 2 78 I 0.41666666666667 0.83333333333333 0.58403632892002 126.904 > 2 79 I 0.66666666666667 0.83333333333333 0.58403632892002 126.904 > 2 80 I 0.91666666666667 0.83333333333333 0.58403632892002 126.904 > 2 81 I 0.16666666666667 0.08333333333333 0.91736966225335 126.904 > 2 82 I 0.41666666666667 0.08333333333333 0.91736966225335 126.904 > 2 83 I 0.66666666666667 0.08333333333333 0.91736966225335 126.904 > 2 84 I 0.91666666666667 0.08333333333333 0.91736966225335 126.904 > 2 85 I 0.16666666666667 0.33333333333333 0.91736966225335 126.904 > 2 86 I 0.41666666666667 0.33333333333333 0.91736966225335 126.904 > 2 87 I 0.66666666666667 0.33333333333333 0.91736966225335 126.904 > 2 88 I 0.91666666666667 0.33333333333333 0.91736966225335 126.904 > 2 89 I 0.16666666666667 0.58333333333333 0.91736966225335 126.904 > 2 90 I 0.41666666666667 0.58333333333333 0.91736966225335 126.904 > 2 91 I 0.66666666666667 0.58333333333333 0.91736966225335 126.904 > 2 92 I 0.91666666666667 0.58333333333333 0.91736966225335 126.904 > 2 93 I 0.16666666666667 0.83333333333333 0.91736966225335 126.904 > 2 94 I 0.41666666666667 0.83333333333333 0.91736966225335 126.904 > 2 95 I 0.66666666666667 0.83333333333333 0.91736966225335 126.904 > 2 96 I 0.91666666666667 0.83333333333333 0.91736966225335 126.904 > 2 97 I 0.08333333333333 0.16666666666667 0.08263033774665 126.904 > 3 98 I 0.33333333333333 0.16666666666667 0.08263033774665 126.904 > 3 99 I 0.58333333333333 0.16666666666667 0.08263033774665 126.904 > 3 100 I 0.83333333333333 0.16666666666667 0.08263033774665 126.904 > 3 101 I 0.08333333333333 0.41666666666667 0.08263033774665 126.904 > 3 102 I 0.33333333333333 0.41666666666667 0.08263033774665 126.904 > 3 103 I 0.58333333333333 0.41666666666667 0.08263033774665 126.904 > 3 104 I 0.83333333333333 0.41666666666667 0.08263033774665 126.904 > 3 105 I 0.08333333333333 0.66666666666667 0.08263033774665 126.904 > 3 106 I 0.33333333333333 0.66666666666667 0.08263033774665 126.904 > 3 107 I 0.58333333333333 0.66666666666667 0.08263033774665 126.904 > 3 108 I 0.83333333333333 0.66666666666667 0.08263033774665 126.904 > 3 109 I 0.08333333333333 0.91666666666667 0.08263033774665 126.904 > 3 110 I 0.33333333333333 0.91666666666667 0.08263033774665 126.904 > 3 111 I 0.58333333333333 0.91666666666667 0.08263033774665 126.904 > 3 112 I 0.83333333333333 0.91666666666667 0.08263033774665 126.904 > 3 113 I 0.08333333333333 0.16666666666667 0.41596367107998 126.904 > 3 114 I 0.33333333333333 0.16666666666667 0.41596367107998 126.904 > 3 115 I 0.58333333333333 0.16666666666667 0.41596367107998 126.904 > 3 116 I 0.83333333333333 0.16666666666667 0.41596367107998 126.904 > 3 117 I 0.08333333333333 0.41666666666667 0.41596367107998 126.904 > 3 118 I 0.33333333333333 0.41666666666667 0.41596367107998 126.904 > 3 119 I 0.58333333333333 0.41666666666667 0.41596367107998 126.904 > 3 120 I 0.83333333333333 0.41666666666667 0.41596367107998 126.904 > 3 121 I 0.08333333333333 0.66666666666667 0.41596367107998 126.904 > 3 122 I 0.33333333333333 0.66666666666667 0.41596367107998 126.904 > 3 123 I 0.58333333333333 0.66666666666667 0.41596367107998 126.904 > 3 124 I 0.83333333333333 0.66666666666667 0.41596367107998 126.904 > 3 125 I 0.08333333333333 0.91666666666667 0.41596367107998 126.904 > 3 126 I 0.33333333333333 0.91666666666667 0.41596367107998 126.904 > 3 127 I 0.58333333333333 0.91666666666667 0.41596367107998 126.904 > 3 128 I 0.83333333333333 0.91666666666667 0.41596367107998 126.904 > 3 129 I 0.08333333333333 0.16666666666667 0.74929700441331 126.904 > 3 130 I 0.33333333333333 0.16666666666667 0.74929700441331 126.904 > 3 131 I 0.58333333333333 0.16666666666667 0.74929700441331 126.904 > 3 132 I 0.83333333333333 0.16666666666667 0.74929700441331 126.904 > 3 133 I 0.08333333333333 0.41666666666667 0.74929700441331 126.904 > 3 134 I 0.33333333333333 0.41666666666667 0.74929700441331 126.904 > 3 135 I 0.58333333333333 0.41666666666667 0.74929700441331 126.904 > 3 136 I 0.83333333333333 0.41666666666667 0.74929700441331 126.904 > 3 137 I 0.08333333333333 0.66666666666667 0.74929700441331 126.904 > 3 138 I 0.33333333333333 0.66666666666667 0.74929700441331 126.904 > 3 139 I 0.58333333333333 0.66666666666667 0.74929700441331 126.904 > 3 140 I 0.83333333333333 0.66666666666667 0.74929700441331 126.904 > 3 141 I 0.08333333333333 0.91666666666667 0.74929700441331 126.904 > 3 142 I 0.33333333333333 0.91666666666667 0.74929700441331 126.904 > 3 143 I 0.58333333333333 0.91666666666667 0.74929700441331 126.904 > 3 144 I 0.83333333333333 0.91666666666667 0.74929700441331 126.904 > 3 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 4.1225517 0.0000000 0.0000000 0.0000000 4.1225517 0.0000000 0.0000000 0.0000000 3.5941550 -------------------------- Born effective charges -------------------------- 1 Yb 2.5320143 0.0000000 0.0000000 0.0000000 2.5320143 0.0000000 0.0000000 0.0000000 1.4835980 2 I -1.2660071 -0.0000000 0.0000000 0.0000000 -1.2660071 0.0000000 0.0000000 0.0000000 -0.7417990 3 I -1.2660071 -0.0000000 0.0000000 0.0000000 -1.2660071 0.0000000 0.0000000 0.0000000 -0.7417990 ---------------------------------------------------------------------------- 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: 2574/2574 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 215 Number of blocks in projector: 165 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: 99 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 68 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 48 Use standard eigh solver. Tree of FC basis block matrices: - (215, 211), data: False |-- (48, 46), data: True |-- (68, 68), data: True |-- (99, 97), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.015 Solver_block: 80 / 80 - Time: 0.171 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.188 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (yy) 0.00000000 (yy) Permutation basis: 432/432 Permutation basis: 2574/2574 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 215 Number of blocks in projector: 165 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: 99 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 68 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 48 Use standard eigh solver. Tree of FC basis block matrices: - (215, 211), data: False |-- (48, 46), data: True |-- (68, 68), data: True |-- (99, 97), data: True Max drift after symmetrization by symfc projector: 0.00000000 (yy) 0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 18:55:38]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 18:55:39]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.382603147621698 0.000000000000000 0.000000000000000 b -2.191301573810849 3.795445660546032 0.000000000000000 c 0.000000000000000 0.000000000000000 7.059576650000000 Atomic positions (fractional): 1 Yb 0.00000000000000 0.00000000000000 0.00000000000000 173.054 2 I 0.66666666666667 0.33333333333333 0.75210898676006 126.904 3 I 0.33333333333333 0.66666666666667 0.24789101323994 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.530412590486790 0.000000000000000 0.000000000000000 b -8.765206295243395 15.181782642184128 0.000000000000000 c 0.000000000000000 0.000000000000000 21.178729950000001 Atomic positions (fractional): 1 Yb 0.00000000000000 0.00000000000000 0.00000000000000 173.054 > 1 2 Yb 0.25000000000000 0.00000000000000 0.00000000000000 173.054 > 1 3 Yb 0.50000000000000 0.00000000000000 0.00000000000000 173.054 > 1 4 Yb 0.75000000000000 0.00000000000000 0.00000000000000 173.054 > 1 5 Yb 0.00000000000000 0.25000000000000 0.00000000000000 173.054 > 1 6 Yb 0.25000000000000 0.25000000000000 0.00000000000000 173.054 > 1 7 Yb 0.50000000000000 0.25000000000000 0.00000000000000 173.054 > 1 8 Yb 0.75000000000000 0.25000000000000 0.00000000000000 173.054 > 1 9 Yb 0.00000000000000 0.50000000000000 0.00000000000000 173.054 > 1 10 Yb 0.25000000000000 0.50000000000000 0.00000000000000 173.054 > 1 11 Yb 0.50000000000000 0.50000000000000 0.00000000000000 173.054 > 1 12 Yb 0.75000000000000 0.50000000000000 0.00000000000000 173.054 > 1 13 Yb 0.00000000000000 0.75000000000000 0.00000000000000 173.054 > 1 14 Yb 0.25000000000000 0.75000000000000 0.00000000000000 173.054 > 1 15 Yb 0.50000000000000 0.75000000000000 0.00000000000000 173.054 > 1 16 Yb 0.75000000000000 0.75000000000000 0.00000000000000 173.054 > 1 17 Yb 0.00000000000000 0.00000000000000 0.33333333333333 173.054 > 1 18 Yb 0.25000000000000 0.00000000000000 0.33333333333333 173.054 > 1 19 Yb 0.50000000000000 0.00000000000000 0.33333333333333 173.054 > 1 20 Yb 0.75000000000000 0.00000000000000 0.33333333333333 173.054 > 1 21 Yb 0.00000000000000 0.25000000000000 0.33333333333333 173.054 > 1 22 Yb 0.25000000000000 0.25000000000000 0.33333333333333 173.054 > 1 23 Yb 0.50000000000000 0.25000000000000 0.33333333333333 173.054 > 1 24 Yb 0.75000000000000 0.25000000000000 0.33333333333333 173.054 > 1 25 Yb 0.00000000000000 0.50000000000000 0.33333333333333 173.054 > 1 26 Yb 0.25000000000000 0.50000000000000 0.33333333333333 173.054 > 1 27 Yb 0.50000000000000 0.50000000000000 0.33333333333333 173.054 > 1 28 Yb 0.75000000000000 0.50000000000000 0.33333333333333 173.054 > 1 29 Yb 0.00000000000000 0.75000000000000 0.33333333333333 173.054 > 1 30 Yb 0.25000000000000 0.75000000000000 0.33333333333333 173.054 > 1 31 Yb 0.50000000000000 0.75000000000000 0.33333333333333 173.054 > 1 32 Yb 0.75000000000000 0.75000000000000 0.33333333333333 173.054 > 1 33 Yb 0.00000000000000 0.00000000000000 0.66666666666667 173.054 > 1 34 Yb 0.25000000000000 0.00000000000000 0.66666666666667 173.054 > 1 35 Yb 0.50000000000000 0.00000000000000 0.66666666666667 173.054 > 1 36 Yb 0.75000000000000 0.00000000000000 0.66666666666667 173.054 > 1 37 Yb 0.00000000000000 0.25000000000000 0.66666666666667 173.054 > 1 38 Yb 0.25000000000000 0.25000000000000 0.66666666666667 173.054 > 1 39 Yb 0.50000000000000 0.25000000000000 0.66666666666667 173.054 > 1 40 Yb 0.75000000000000 0.25000000000000 0.66666666666667 173.054 > 1 41 Yb 0.00000000000000 0.50000000000000 0.66666666666667 173.054 > 1 42 Yb 0.25000000000000 0.50000000000000 0.66666666666667 173.054 > 1 43 Yb 0.50000000000000 0.50000000000000 0.66666666666667 173.054 > 1 44 Yb 0.75000000000000 0.50000000000000 0.66666666666667 173.054 > 1 45 Yb 0.00000000000000 0.75000000000000 0.66666666666667 173.054 > 1 46 Yb 0.25000000000000 0.75000000000000 0.66666666666667 173.054 > 1 47 Yb 0.50000000000000 0.75000000000000 0.66666666666667 173.054 > 1 48 Yb 0.75000000000000 0.75000000000000 0.66666666666667 173.054 > 1 49 I 0.16666666666667 0.08333333333333 0.25070299558669 126.904 > 49 50 I 0.41666666666667 0.08333333333333 0.25070299558669 126.904 > 49 51 I 0.66666666666667 0.08333333333333 0.25070299558669 126.904 > 49 52 I 0.91666666666667 0.08333333333333 0.25070299558669 126.904 > 49 53 I 0.16666666666667 0.33333333333333 0.25070299558669 126.904 > 49 54 I 0.41666666666667 0.33333333333333 0.25070299558669 126.904 > 49 55 I 0.66666666666667 0.33333333333333 0.25070299558669 126.904 > 49 56 I 0.91666666666667 0.33333333333333 0.25070299558669 126.904 > 49 57 I 0.16666666666667 0.58333333333333 0.25070299558669 126.904 > 49 58 I 0.41666666666667 0.58333333333333 0.25070299558669 126.904 > 49 59 I 0.66666666666667 0.58333333333333 0.25070299558669 126.904 > 49 60 I 0.91666666666667 0.58333333333333 0.25070299558669 126.904 > 49 61 I 0.16666666666667 0.83333333333333 0.25070299558669 126.904 > 49 62 I 0.41666666666667 0.83333333333333 0.25070299558669 126.904 > 49 63 I 0.66666666666667 0.83333333333333 0.25070299558669 126.904 > 49 64 I 0.91666666666667 0.83333333333333 0.25070299558669 126.904 > 49 65 I 0.16666666666667 0.08333333333333 0.58403632892002 126.904 > 49 66 I 0.41666666666667 0.08333333333333 0.58403632892002 126.904 > 49 67 I 0.66666666666667 0.08333333333333 0.58403632892002 126.904 > 49 68 I 0.91666666666667 0.08333333333333 0.58403632892002 126.904 > 49 69 I 0.16666666666667 0.33333333333333 0.58403632892002 126.904 > 49 70 I 0.41666666666667 0.33333333333333 0.58403632892002 126.904 > 49 71 I 0.66666666666667 0.33333333333333 0.58403632892002 126.904 > 49 72 I 0.91666666666667 0.33333333333333 0.58403632892002 126.904 > 49 73 I 0.16666666666667 0.58333333333333 0.58403632892002 126.904 > 49 74 I 0.41666666666667 0.58333333333333 0.58403632892002 126.904 > 49 75 I 0.66666666666667 0.58333333333333 0.58403632892002 126.904 > 49 76 I 0.91666666666667 0.58333333333333 0.58403632892002 126.904 > 49 77 I 0.16666666666667 0.83333333333333 0.58403632892002 126.904 > 49 78 I 0.41666666666667 0.83333333333333 0.58403632892002 126.904 > 49 79 I 0.66666666666667 0.83333333333333 0.58403632892002 126.904 > 49 80 I 0.91666666666667 0.83333333333333 0.58403632892002 126.904 > 49 81 I 0.16666666666667 0.08333333333333 0.91736966225335 126.904 > 49 82 I 0.41666666666667 0.08333333333333 0.91736966225335 126.904 > 49 83 I 0.66666666666667 0.08333333333333 0.91736966225335 126.904 > 49 84 I 0.91666666666667 0.08333333333333 0.91736966225335 126.904 > 49 85 I 0.16666666666667 0.33333333333333 0.91736966225335 126.904 > 49 86 I 0.41666666666667 0.33333333333333 0.91736966225335 126.904 > 49 87 I 0.66666666666667 0.33333333333333 0.91736966225335 126.904 > 49 88 I 0.91666666666667 0.33333333333333 0.91736966225335 126.904 > 49 89 I 0.16666666666667 0.58333333333333 0.91736966225335 126.904 > 49 90 I 0.41666666666667 0.58333333333333 0.91736966225335 126.904 > 49 91 I 0.66666666666667 0.58333333333333 0.91736966225335 126.904 > 49 92 I 0.91666666666667 0.58333333333333 0.91736966225335 126.904 > 49 93 I 0.16666666666667 0.83333333333333 0.91736966225335 126.904 > 49 94 I 0.41666666666667 0.83333333333333 0.91736966225335 126.904 > 49 95 I 0.66666666666667 0.83333333333333 0.91736966225335 126.904 > 49 96 I 0.91666666666667 0.83333333333333 0.91736966225335 126.904 > 49 97 I 0.08333333333333 0.16666666666667 0.08263033774665 126.904 > 97 98 I 0.33333333333333 0.16666666666667 0.08263033774665 126.904 > 97 99 I 0.58333333333333 0.16666666666667 0.08263033774665 126.904 > 97 100 I 0.83333333333333 0.16666666666667 0.08263033774665 126.904 > 97 101 I 0.08333333333333 0.41666666666667 0.08263033774665 126.904 > 97 102 I 0.33333333333333 0.41666666666667 0.08263033774665 126.904 > 97 103 I 0.58333333333333 0.41666666666667 0.08263033774665 126.904 > 97 104 I 0.83333333333333 0.41666666666667 0.08263033774665 126.904 > 97 105 I 0.08333333333333 0.66666666666667 0.08263033774665 126.904 > 97 106 I 0.33333333333333 0.66666666666667 0.08263033774665 126.904 > 97 107 I 0.58333333333333 0.66666666666667 0.08263033774665 126.904 > 97 108 I 0.83333333333333 0.66666666666667 0.08263033774665 126.904 > 97 109 I 0.08333333333333 0.91666666666667 0.08263033774665 126.904 > 97 110 I 0.33333333333333 0.91666666666667 0.08263033774665 126.904 > 97 111 I 0.58333333333333 0.91666666666667 0.08263033774665 126.904 > 97 112 I 0.83333333333333 0.91666666666667 0.08263033774665 126.904 > 97 113 I 0.08333333333333 0.16666666666667 0.41596367107998 126.904 > 97 114 I 0.33333333333333 0.16666666666667 0.41596367107998 126.904 > 97 115 I 0.58333333333333 0.16666666666667 0.41596367107998 126.904 > 97 116 I 0.83333333333333 0.16666666666667 0.41596367107998 126.904 > 97 117 I 0.08333333333333 0.41666666666667 0.41596367107998 126.904 > 97 118 I 0.33333333333333 0.41666666666667 0.41596367107998 126.904 > 97 119 I 0.58333333333333 0.41666666666667 0.41596367107998 126.904 > 97 120 I 0.83333333333333 0.41666666666667 0.41596367107998 126.904 > 97 121 I 0.08333333333333 0.66666666666667 0.41596367107998 126.904 > 97 122 I 0.33333333333333 0.66666666666667 0.41596367107998 126.904 > 97 123 I 0.58333333333333 0.66666666666667 0.41596367107998 126.904 > 97 124 I 0.83333333333333 0.66666666666667 0.41596367107998 126.904 > 97 125 I 0.08333333333333 0.91666666666667 0.41596367107998 126.904 > 97 126 I 0.33333333333333 0.91666666666667 0.41596367107998 126.904 > 97 127 I 0.58333333333333 0.91666666666667 0.41596367107998 126.904 > 97 128 I 0.83333333333333 0.91666666666667 0.41596367107998 126.904 > 97 129 I 0.08333333333333 0.16666666666667 0.74929700441331 126.904 > 97 130 I 0.33333333333333 0.16666666666667 0.74929700441331 126.904 > 97 131 I 0.58333333333333 0.16666666666667 0.74929700441331 126.904 > 97 132 I 0.83333333333333 0.16666666666667 0.74929700441331 126.904 > 97 133 I 0.08333333333333 0.41666666666667 0.74929700441331 126.904 > 97 134 I 0.33333333333333 0.41666666666667 0.74929700441331 126.904 > 97 135 I 0.58333333333333 0.41666666666667 0.74929700441331 126.904 > 97 136 I 0.83333333333333 0.41666666666667 0.74929700441331 126.904 > 97 137 I 0.08333333333333 0.66666666666667 0.74929700441331 126.904 > 97 138 I 0.33333333333333 0.66666666666667 0.74929700441331 126.904 > 97 139 I 0.58333333333333 0.66666666666667 0.74929700441331 126.904 > 97 140 I 0.83333333333333 0.66666666666667 0.74929700441331 126.904 > 97 141 I 0.08333333333333 0.91666666666667 0.74929700441331 126.904 > 97 142 I 0.33333333333333 0.91666666666667 0.74929700441331 126.904 > 97 143 I 0.58333333333333 0.91666666666667 0.74929700441331 126.904 > 97 144 I 0.83333333333333 0.91666666666667 0.74929700441331 126.904 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 4.1225517 0.0000000 0.0000000 0.0000000 4.1225517 0.0000000 0.0000000 0.0000000 3.5941550 -------------------------- Born effective charges -------------------------- 1 Yb 2.5320143 0.0000000 0.0000000 0.0000000 2.5320143 0.0000000 0.0000000 0.0000000 1.4835980 2 I -1.2660071 -0.0000000 0.0000000 0.0000000 -1.2660071 0.0000000 0.0000000 0.0000000 -0.7417990 3 I -1.2660071 -0.0000000 0.0000000 0.0000000 -1.2660071 0.0000000 0.0000000 0.0000000 -0.7417990 ---------------------------------------------------------------------------- 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[ 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] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000002 (yyy) 0.00000002 (yyy) 0.00000002 (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 18:55:44]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 18:55: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: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-3m1 (164) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.382603147621698 0.000000000000000 0.000000000000000 b -2.191301573810849 3.795445660546032 0.000000000000000 c 0.000000000000000 0.000000000000000 7.059576650000000 Atomic positions (fractional): 1 Yb 0.00000000000000 0.00000000000000 0.00000000000000 173.054 2 I 0.66666666666667 0.33333333333333 0.75210898676006 126.904 3 I 0.33333333333333 0.66666666666667 0.24789101323994 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.530412590486790 0.000000000000000 0.000000000000000 b -8.765206295243395 15.181782642184128 0.000000000000000 c 0.000000000000000 0.000000000000000 21.178729950000001 Atomic positions (fractional): 1 Yb 0.00000000000000 0.00000000000000 0.00000000000000 173.054 > 1 2 Yb 0.25000000000000 0.00000000000000 0.00000000000000 173.054 > 1 3 Yb 0.50000000000000 0.00000000000000 0.00000000000000 173.054 > 1 4 Yb 0.75000000000000 0.00000000000000 0.00000000000000 173.054 > 1 5 Yb 0.00000000000000 0.25000000000000 0.00000000000000 173.054 > 1 6 Yb 0.25000000000000 0.25000000000000 0.00000000000000 173.054 > 1 7 Yb 0.50000000000000 0.25000000000000 0.00000000000000 173.054 > 1 8 Yb 0.75000000000000 0.25000000000000 0.00000000000000 173.054 > 1 9 Yb 0.00000000000000 0.50000000000000 0.00000000000000 173.054 > 1 10 Yb 0.25000000000000 0.50000000000000 0.00000000000000 173.054 > 1 11 Yb 0.50000000000000 0.50000000000000 0.00000000000000 173.054 > 1 12 Yb 0.75000000000000 0.50000000000000 0.00000000000000 173.054 > 1 13 Yb 0.00000000000000 0.75000000000000 0.00000000000000 173.054 > 1 14 Yb 0.25000000000000 0.75000000000000 0.00000000000000 173.054 > 1 15 Yb 0.50000000000000 0.75000000000000 0.00000000000000 173.054 > 1 16 Yb 0.75000000000000 0.75000000000000 0.00000000000000 173.054 > 1 17 Yb 0.00000000000000 0.00000000000000 0.33333333333333 173.054 > 1 18 Yb 0.25000000000000 0.00000000000000 0.33333333333333 173.054 > 1 19 Yb 0.50000000000000 0.00000000000000 0.33333333333333 173.054 > 1 20 Yb 0.75000000000000 0.00000000000000 0.33333333333333 173.054 > 1 21 Yb 0.00000000000000 0.25000000000000 0.33333333333333 173.054 > 1 22 Yb 0.25000000000000 0.25000000000000 0.33333333333333 173.054 > 1 23 Yb 0.50000000000000 0.25000000000000 0.33333333333333 173.054 > 1 24 Yb 0.75000000000000 0.25000000000000 0.33333333333333 173.054 > 1 25 Yb 0.00000000000000 0.50000000000000 0.33333333333333 173.054 > 1 26 Yb 0.25000000000000 0.50000000000000 0.33333333333333 173.054 > 1 27 Yb 0.50000000000000 0.50000000000000 0.33333333333333 173.054 > 1 28 Yb 0.75000000000000 0.50000000000000 0.33333333333333 173.054 > 1 29 Yb 0.00000000000000 0.75000000000000 0.33333333333333 173.054 > 1 30 Yb 0.25000000000000 0.75000000000000 0.33333333333333 173.054 > 1 31 Yb 0.50000000000000 0.75000000000000 0.33333333333333 173.054 > 1 32 Yb 0.75000000000000 0.75000000000000 0.33333333333333 173.054 > 1 33 Yb 0.00000000000000 0.00000000000000 0.66666666666667 173.054 > 1 34 Yb 0.25000000000000 0.00000000000000 0.66666666666667 173.054 > 1 35 Yb 0.50000000000000 0.00000000000000 0.66666666666667 173.054 > 1 36 Yb 0.75000000000000 0.00000000000000 0.66666666666667 173.054 > 1 37 Yb 0.00000000000000 0.25000000000000 0.66666666666667 173.054 > 1 38 Yb 0.25000000000000 0.25000000000000 0.66666666666667 173.054 > 1 39 Yb 0.50000000000000 0.25000000000000 0.66666666666667 173.054 > 1 40 Yb 0.75000000000000 0.25000000000000 0.66666666666667 173.054 > 1 41 Yb 0.00000000000000 0.50000000000000 0.66666666666667 173.054 > 1 42 Yb 0.25000000000000 0.50000000000000 0.66666666666667 173.054 > 1 43 Yb 0.50000000000000 0.50000000000000 0.66666666666667 173.054 > 1 44 Yb 0.75000000000000 0.50000000000000 0.66666666666667 173.054 > 1 45 Yb 0.00000000000000 0.75000000000000 0.66666666666667 173.054 > 1 46 Yb 0.25000000000000 0.75000000000000 0.66666666666667 173.054 > 1 47 Yb 0.50000000000000 0.75000000000000 0.66666666666667 173.054 > 1 48 Yb 0.75000000000000 0.75000000000000 0.66666666666667 173.054 > 1 49 I 0.16666666666667 0.08333333333333 0.25070299558669 126.904 > 49 50 I 0.41666666666667 0.08333333333333 0.25070299558669 126.904 > 49 51 I 0.66666666666667 0.08333333333333 0.25070299558669 126.904 > 49 52 I 0.91666666666667 0.08333333333333 0.25070299558669 126.904 > 49 53 I 0.16666666666667 0.33333333333333 0.25070299558669 126.904 > 49 54 I 0.41666666666667 0.33333333333333 0.25070299558669 126.904 > 49 55 I 0.66666666666667 0.33333333333333 0.25070299558669 126.904 > 49 56 I 0.91666666666667 0.33333333333333 0.25070299558669 126.904 > 49 57 I 0.16666666666667 0.58333333333333 0.25070299558669 126.904 > 49 58 I 0.41666666666667 0.58333333333333 0.25070299558669 126.904 > 49 59 I 0.66666666666667 0.58333333333333 0.25070299558669 126.904 > 49 60 I 0.91666666666667 0.58333333333333 0.25070299558669 126.904 > 49 61 I 0.16666666666667 0.83333333333333 0.25070299558669 126.904 > 49 62 I 0.41666666666667 0.83333333333333 0.25070299558669 126.904 > 49 63 I 0.66666666666667 0.83333333333333 0.25070299558669 126.904 > 49 64 I 0.91666666666667 0.83333333333333 0.25070299558669 126.904 > 49 65 I 0.16666666666667 0.08333333333333 0.58403632892002 126.904 > 49 66 I 0.41666666666667 0.08333333333333 0.58403632892002 126.904 > 49 67 I 0.66666666666667 0.08333333333333 0.58403632892002 126.904 > 49 68 I 0.91666666666667 0.08333333333333 0.58403632892002 126.904 > 49 69 I 0.16666666666667 0.33333333333333 0.58403632892002 126.904 > 49 70 I 0.41666666666667 0.33333333333333 0.58403632892002 126.904 > 49 71 I 0.66666666666667 0.33333333333333 0.58403632892002 126.904 > 49 72 I 0.91666666666667 0.33333333333333 0.58403632892002 126.904 > 49 73 I 0.16666666666667 0.58333333333333 0.58403632892002 126.904 > 49 74 I 0.41666666666667 0.58333333333333 0.58403632892002 126.904 > 49 75 I 0.66666666666667 0.58333333333333 0.58403632892002 126.904 > 49 76 I 0.91666666666667 0.58333333333333 0.58403632892002 126.904 > 49 77 I 0.16666666666667 0.83333333333333 0.58403632892002 126.904 > 49 78 I 0.41666666666667 0.83333333333333 0.58403632892002 126.904 > 49 79 I 0.66666666666667 0.83333333333333 0.58403632892002 126.904 > 49 80 I 0.91666666666667 0.83333333333333 0.58403632892002 126.904 > 49 81 I 0.16666666666667 0.08333333333333 0.91736966225335 126.904 > 49 82 I 0.41666666666667 0.08333333333333 0.91736966225335 126.904 > 49 83 I 0.66666666666667 0.08333333333333 0.91736966225335 126.904 > 49 84 I 0.91666666666667 0.08333333333333 0.91736966225335 126.904 > 49 85 I 0.16666666666667 0.33333333333333 0.91736966225335 126.904 > 49 86 I 0.41666666666667 0.33333333333333 0.91736966225335 126.904 > 49 87 I 0.66666666666667 0.33333333333333 0.91736966225335 126.904 > 49 88 I 0.91666666666667 0.33333333333333 0.91736966225335 126.904 > 49 89 I 0.16666666666667 0.58333333333333 0.91736966225335 126.904 > 49 90 I 0.41666666666667 0.58333333333333 0.91736966225335 126.904 > 49 91 I 0.66666666666667 0.58333333333333 0.91736966225335 126.904 > 49 92 I 0.91666666666667 0.58333333333333 0.91736966225335 126.904 > 49 93 I 0.16666666666667 0.83333333333333 0.91736966225335 126.904 > 49 94 I 0.41666666666667 0.83333333333333 0.91736966225335 126.904 > 49 95 I 0.66666666666667 0.83333333333333 0.91736966225335 126.904 > 49 96 I 0.91666666666667 0.83333333333333 0.91736966225335 126.904 > 49 97 I 0.08333333333333 0.16666666666667 0.08263033774665 126.904 > 97 98 I 0.33333333333333 0.16666666666667 0.08263033774665 126.904 > 97 99 I 0.58333333333333 0.16666666666667 0.08263033774665 126.904 > 97 100 I 0.83333333333333 0.16666666666667 0.08263033774665 126.904 > 97 101 I 0.08333333333333 0.41666666666667 0.08263033774665 126.904 > 97 102 I 0.33333333333333 0.41666666666667 0.08263033774665 126.904 > 97 103 I 0.58333333333333 0.41666666666667 0.08263033774665 126.904 > 97 104 I 0.83333333333333 0.41666666666667 0.08263033774665 126.904 > 97 105 I 0.08333333333333 0.66666666666667 0.08263033774665 126.904 > 97 106 I 0.33333333333333 0.66666666666667 0.08263033774665 126.904 > 97 107 I 0.58333333333333 0.66666666666667 0.08263033774665 126.904 > 97 108 I 0.83333333333333 0.66666666666667 0.08263033774665 126.904 > 97 109 I 0.08333333333333 0.91666666666667 0.08263033774665 126.904 > 97 110 I 0.33333333333333 0.91666666666667 0.08263033774665 126.904 > 97 111 I 0.58333333333333 0.91666666666667 0.08263033774665 126.904 > 97 112 I 0.83333333333333 0.91666666666667 0.08263033774665 126.904 > 97 113 I 0.08333333333333 0.16666666666667 0.41596367107998 126.904 > 97 114 I 0.33333333333333 0.16666666666667 0.41596367107998 126.904 > 97 115 I 0.58333333333333 0.16666666666667 0.41596367107998 126.904 > 97 116 I 0.83333333333333 0.16666666666667 0.41596367107998 126.904 > 97 117 I 0.08333333333333 0.41666666666667 0.41596367107998 126.904 > 97 118 I 0.33333333333333 0.41666666666667 0.41596367107998 126.904 > 97 119 I 0.58333333333333 0.41666666666667 0.41596367107998 126.904 > 97 120 I 0.83333333333333 0.41666666666667 0.41596367107998 126.904 > 97 121 I 0.08333333333333 0.66666666666667 0.41596367107998 126.904 > 97 122 I 0.33333333333333 0.66666666666667 0.41596367107998 126.904 > 97 123 I 0.58333333333333 0.66666666666667 0.41596367107998 126.904 > 97 124 I 0.83333333333333 0.66666666666667 0.41596367107998 126.904 > 97 125 I 0.08333333333333 0.91666666666667 0.41596367107998 126.904 > 97 126 I 0.33333333333333 0.91666666666667 0.41596367107998 126.904 > 97 127 I 0.58333333333333 0.91666666666667 0.41596367107998 126.904 > 97 128 I 0.83333333333333 0.91666666666667 0.41596367107998 126.904 > 97 129 I 0.08333333333333 0.16666666666667 0.74929700441331 126.904 > 97 130 I 0.33333333333333 0.16666666666667 0.74929700441331 126.904 > 97 131 I 0.58333333333333 0.16666666666667 0.74929700441331 126.904 > 97 132 I 0.83333333333333 0.16666666666667 0.74929700441331 126.904 > 97 133 I 0.08333333333333 0.41666666666667 0.74929700441331 126.904 > 97 134 I 0.33333333333333 0.41666666666667 0.74929700441331 126.904 > 97 135 I 0.58333333333333 0.41666666666667 0.74929700441331 126.904 > 97 136 I 0.83333333333333 0.41666666666667 0.74929700441331 126.904 > 97 137 I 0.08333333333333 0.66666666666667 0.74929700441331 126.904 > 97 138 I 0.33333333333333 0.66666666666667 0.74929700441331 126.904 > 97 139 I 0.58333333333333 0.66666666666667 0.74929700441331 126.904 > 97 140 I 0.83333333333333 0.66666666666667 0.74929700441331 126.904 > 97 141 I 0.08333333333333 0.91666666666667 0.74929700441331 126.904 > 97 142 I 0.33333333333333 0.91666666666667 0.74929700441331 126.904 > 97 143 I 0.58333333333333 0.91666666666667 0.74929700441331 126.904 > 97 144 I 0.83333333333333 0.91666666666667 0.74929700441331 126.904 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 4.1225517 0.0000000 0.0000000 0.0000000 4.1225517 0.0000000 0.0000000 0.0000000 3.5941550 -------------------------- Born effective charges -------------------------- 1 Yb 2.5320143 0.0000000 0.0000000 0.0000000 2.5320143 0.0000000 0.0000000 0.0000000 1.4835980 2 I -1.2660071 -0.0000000 0.0000000 0.0000000 -1.2660071 0.0000000 0.0000000 0.0000000 -0.7417990 3 I -1.2660071 -0.0000000 0.0000000 0.0000000 -1.2660071 0.0000000 0.0000000 0.0000000 -0.7417990 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000002 (yyy) 0.00000002 (yyy) 0.00000002 (yyy) Max drift of fc2: 0.00000000 (yy) 0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 13 13 7 ] 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.85, Number of G-points: 311, Lambda: 0.18 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/126) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 126 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 2.258 ( 0.000 0.000 0.000) 0.000 2.258 ( 0.000 0.000 0.000) 0.000 2.831 ( 0.000 0.000 0.000) 0.000 2.897 ( 0.000 0.000 0.000) 0.000 2.897 ( 0.000 0.000 0.000) 0.000 3.798 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/126) ======================= q-point: ( 0.08 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.143 ( 6.695 3.866 0.532) 7.749 0.278 ( 11.580 6.686 1.156) 13.421 0.453 ( 18.937 10.933 -0.848) 21.883 2.251 ( -0.515 -0.297 -0.278) 0.656 2.294 ( 2.961 1.710 -0.782) 3.508 2.817 ( -1.133 -0.654 0.003) 1.309 2.905 ( 0.598 0.345 0.031) 0.691 3.729 ( -4.148 -2.395 -20.180) 20.740 3.764 ( -0.708 -0.409 20.664) 20.680 ======================= Grid point 2 (3/126) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.317 ( 8.105 4.679 0.664) 9.382 0.536 ( 10.376 5.991 0.989) 12.022 0.876 ( 17.006 9.818 -0.510) 19.643 2.237 ( -0.641 -0.370 -0.529) 0.909 2.383 ( 4.242 2.449 -1.093) 5.019 2.781 ( -1.798 -1.038 -0.023) 2.077 2.922 ( 0.794 0.458 0.152) 0.929 3.574 ( -8.275 -4.778 -5.733) 11.143 3.761 ( 0.003 0.002 6.665) 6.665 ======================= Grid point 3 (4/126) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.515 ( 8.527 4.923 1.389) 9.943 0.759 ( 8.582 4.955 0.996) 9.959 1.240 ( 13.835 7.987 0.084) 15.975 2.223 ( -0.553 -0.319 -0.639) 0.903 2.475 ( 3.219 1.858 -0.956) 3.837 2.741 ( -1.415 -0.817 -0.280) 1.657 2.940 ( 0.694 0.401 0.240) 0.837 3.367 ( -8.814 -5.089 -2.364) 10.449 3.758 ( -0.301 -0.174 3.680) 3.697 ======================= Grid point 4 (5/126) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.706 ( 7.576 4.374 2.257) 9.034 0.935 ( 6.421 3.707 1.154) 7.504 1.515 ( 9.541 5.509 0.873) 11.052 2.210 ( -0.545 -0.315 -0.562) 0.843 2.518 ( 0.370 0.213 -0.967) 1.057 2.724 ( 0.059 0.034 -0.808) 0.811 2.954 ( 0.556 0.321 0.166) 0.663 3.188 ( -6.019 -3.475 -0.653) 6.981 3.745 ( -0.871 -0.503 2.321) 2.530 ======================= Grid point 5 (6/126) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.859 ( 5.299 3.059 2.328) 6.546 1.057 ( 4.020 2.321 1.085) 4.767 1.684 ( 4.952 2.859 1.331) 5.871 2.197 ( -0.539 -0.311 -0.365) 0.721 2.499 ( -1.664 -0.961 -1.323) 2.333 2.741 ( 1.233 0.712 -0.945) 1.708 2.966 ( 0.416 0.240 0.019) 0.481 3.094 ( -2.202 -1.271 0.212) 2.551 3.720 ( -1.185 -0.684 1.444) 1.990 ======================= Grid point 6 (7/126) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.945 ( 1.921 1.109 1.019) 2.441 1.121 ( 1.383 0.798 0.464) 1.663 1.755 ( 1.400 0.808 0.680) 1.754 2.187 ( -0.249 -0.144 -0.125) 0.313 2.462 ( -1.003 -0.579 -0.699) 1.352 2.770 ( 0.867 0.500 -0.446) 1.096 2.973 ( 0.166 0.096 -0.022) 0.192 3.068 ( -0.414 -0.239 0.202) 0.518 3.697 ( -0.596 -0.344 0.526) 0.866 ======================= Grid point 14 (8/126) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.256 ( 4.149 7.186 0.000) 8.298 0.489 ( 6.943 12.026 0.000) 13.886 0.758 ( 9.809 16.989 0.000) 19.617 2.265 ( 0.240 0.416 0.000) 0.481 2.328 ( 1.530 2.651 0.000) 3.061 2.792 ( -0.998 -1.729 0.000) 1.997 2.922 ( 0.709 1.227 0.000) 1.417 3.648 ( -3.837 -6.645 0.000) 7.673 3.738 ( -0.489 -0.847 0.000) 0.977 ======================= Grid point 15 (9/126) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.423 ( 6.979 5.930 0.599) 9.178 0.726 ( 6.008 11.526 0.544) 13.010 1.099 ( 12.018 11.448 -0.119) 16.599 2.259 ( -1.968 2.066 -0.312) 2.870 2.400 ( 4.647 -0.358 -0.547) 4.693 2.752 ( -1.335 -1.725 -0.064) 2.182 2.950 ( 0.432 2.111 0.088) 2.157 3.476 ( -8.132 -5.104 -2.343) 9.883 3.728 ( 1.136 -2.681 2.972) 4.161 ======================= Grid point 16 (10/126) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.609 ( 7.774 4.580 1.498) 9.146 0.925 ( 3.801 10.597 0.789) 11.285 1.395 ( 10.681 7.030 0.412) 12.793 2.244 ( -2.113 2.158 -0.425) 3.050 2.464 ( 3.815 -2.463 -0.657) 4.589 2.719 ( -0.483 -1.238 -0.356) 1.376 2.978 ( -0.162 2.992 0.059) 2.997 3.284 ( -8.293 -3.100 -1.225) 8.938 3.718 ( 1.210 -3.501 2.347) 4.385 ======================= Grid point 17 (11/126) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.778 ( 6.831 3.004 2.146) 7.765 1.071 ( 1.460 9.024 0.874) 9.183 1.607 ( 7.254 3.512 1.091) 8.133 2.229 ( -1.960 1.985 -0.290) 2.805 2.479 ( 1.349 -3.358 -1.066) 3.773 2.716 ( 1.148 -0.648 -0.697) 1.491 3.003 ( -0.854 3.678 -0.103) 3.778 3.141 ( -5.257 -1.216 -0.158) 5.398 3.699 ( 0.827 -3.852 1.646) 4.270 ======================= Grid point 18 (12/126) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.899 ( 4.305 1.135 1.689) 4.762 1.159 ( -0.823 7.118 0.586) 7.190 1.722 ( 3.435 0.983 1.078) 3.732 2.217 ( -1.714 2.248 -0.101) 2.828 2.454 ( -0.059 -2.831 -1.108) 3.041 2.740 ( 1.950 -0.608 -0.641) 2.140 3.020 ( -1.847 4.209 -0.225) 4.602 3.074 ( -1.560 -0.723 0.348) 1.754 3.674 ( 0.855 -3.679 0.921) 3.888 ======================= Grid point 19 (13/126) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.944 ( 0.615 -1.065 -0.000) 1.229 1.187 ( -3.144 5.445 -0.000) 6.287 1.756 ( 0.431 -0.747 0.000) 0.863 2.214 ( -1.527 2.645 -0.000) 3.054 2.437 ( 0.992 -1.719 0.000) 1.985 2.756 ( 0.861 -1.492 0.000) 1.723 3.025 ( -2.550 4.416 -0.000) 5.099 3.059 ( 0.362 -0.628 0.000) 0.725 3.661 ( 1.762 -3.052 -0.000) 3.524 ======================= Grid point 29 (14/126) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.558 ( 4.156 7.199 0.000) 8.312 0.961 ( 6.311 10.931 0.000) 12.622 1.324 ( 6.087 10.543 0.000) 12.174 2.299 ( 0.750 1.299 0.000) 1.500 2.389 ( -0.150 -0.259 0.000) 0.299 2.718 ( -0.872 -1.510 0.000) 1.744 2.996 ( 1.295 2.243 0.000) 2.590 3.361 ( -3.524 -6.104 0.000) 7.048 3.673 ( -1.416 -2.453 0.000) 2.833 ======================= Grid point 30 (15/126) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.705 ( 5.333 4.922 0.894) 7.312 1.152 ( 2.535 10.990 0.283) 11.282 1.519 ( 6.752 5.138 0.440) 8.496 2.299 ( -2.782 2.972 0.064) 4.072 2.400 ( 3.052 -3.255 -0.563) 4.498 2.698 ( 0.051 -0.715 -0.228) 0.752 3.034 ( 0.477 1.954 -0.034) 2.011 3.241 ( -5.538 -1.194 -0.424) 5.681 3.631 ( 0.686 -4.682 0.964) 4.830 ======================= Grid point 31 (16/126) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.834 ( 4.690 2.575 1.321) 5.511 1.271 ( -1.033 10.018 0.236) 10.074 1.657 ( 5.071 1.498 0.907) 5.365 2.281 ( -2.266 2.937 0.249) 3.718 2.412 ( 1.511 -2.715 -1.070) 3.286 2.705 ( 1.204 -0.438 -0.446) 1.357 3.049 ( 0.046 -0.074 0.029) 0.092 3.158 ( -4.816 3.185 -0.200) 5.777 3.600 ( 1.656 -5.522 0.898) 5.835 ======================= Grid point 32 (17/126) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.911 ( 2.242 0.201 0.654) 2.344 1.323 ( -3.553 8.582 0.066) 9.289 1.725 ( 2.231 -0.634 0.504) 2.374 2.281 ( -1.869 3.707 0.220) 4.157 2.402 ( 0.532 -2.022 -0.666) 2.195 2.722 ( 1.330 -1.153 -0.261) 1.780 3.042 ( 0.359 -1.343 0.105) 1.394 3.131 ( -3.517 5.259 -0.092) 6.327 3.580 ( 2.385 -5.260 0.363) 5.787 ======================= Grid point 43 (18/126) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.803 ( 2.640 4.572 0.000) 5.280 1.347 ( 4.488 7.774 0.000) 8.976 1.607 ( 2.149 3.723 0.000) 4.299 2.344 ( -0.468 -0.810 0.000) 0.935 2.357 ( 0.764 1.324 0.000) 1.529 2.692 ( 0.018 0.031 0.000) 0.036 3.046 ( -0.354 -0.613 0.000) 0.707 3.232 ( 0.073 0.127 0.000) 0.146 3.543 ( -2.131 -3.691 0.000) 4.262 ======================= Grid point 44 (19/126) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.882 ( 2.168 2.036 0.380) 2.998 1.457 ( -0.909 7.937 -0.109) 7.990 1.671 ( 2.876 -0.085 0.424) 2.908 2.335 ( -0.542 2.132 0.607) 2.282 2.382 ( 0.481 -0.146 -0.861) 0.997 2.695 ( 0.520 -0.611 -0.156) 0.817 3.024 ( -0.358 -1.746 0.086) 1.784 3.236 ( -2.139 3.895 -0.168) 4.447 3.489 ( 1.026 -5.195 0.304) 5.304 ======================= Grid point 45 (20/126) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.912 ( 0.016 -0.028 -0.000) 0.032 1.487 ( -4.090 7.084 -0.000) 8.180 1.697 ( 1.219 -2.112 0.000) 2.438 2.352 ( -1.609 2.786 -0.000) 3.217 2.377 ( 0.238 -0.413 0.000) 0.477 2.696 ( 0.692 -1.199 0.000) 1.384 3.011 ( 0.853 -1.478 0.000) 1.707 3.236 ( -2.836 4.911 -0.000) 5.671 3.472 ( 2.966 -5.137 -0.000) 5.932 ======================= Grid point 58 (21/126) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.912 ( 0.476 0.825 0.000) 0.953 1.580 ( 2.277 3.943 0.000) 4.553 1.656 ( -0.706 -1.222 0.000) 1.411 2.370 ( 0.824 1.427 0.000) 1.648 2.388 ( -0.023 -0.040 0.000) 0.046 2.682 ( -0.266 -0.461 0.000) 0.533 2.993 ( -0.564 -0.977 0.000) 1.128 3.307 ( 1.693 2.933 0.000) 3.387 3.399 ( -1.976 -3.423 0.000) 3.953 ======================= Grid point 181 (22/126) ======================= q-point: ( 0.00 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.129 ( 0.000 -0.000 5.974) 5.974 0.129 ( 0.000 0.000 5.974) 5.974 0.239 ( 0.000 0.000 11.153) 11.153 2.253 ( -0.000 -0.000 -0.411) 0.411 2.253 ( 0.000 0.000 -0.411) 0.411 2.819 ( 0.000 -0.000 -1.138) 1.138 2.899 ( 0.000 -0.000 0.164) 0.164 2.899 ( 0.000 0.000 0.164) 0.164 4.084 ( 0.000 -0.000 0.191) 0.191 ======================= Grid point 182 (23/126) ======================= q-point: ( 0.08 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.260 ( 3.764 2.173 7.997) 9.102 0.323 ( 10.472 6.046 2.931) 12.442 0.465 ( 17.541 10.127 1.922) 20.345 2.242 ( -0.698 -0.403 -0.557) 0.980 2.289 ( 2.579 1.489 0.099) 2.980 2.804 ( -1.136 -0.656 -1.242) 1.806 2.907 ( 0.639 0.369 0.161) 0.756 3.286 ( 14.773 8.529 -16.445) 23.694 3.969 ( -6.965 -4.021 3.361) 8.717 ======================= Grid point 183 (24/126) ======================= q-point: ( 0.15 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.379 ( 6.520 3.764 4.757) 8.905 0.565 ( 9.907 5.720 1.673) 11.561 0.874 ( 16.813 9.707 0.241) 19.415 2.224 ( -0.709 -0.409 -0.600) 1.015 2.371 ( 4.149 2.395 -0.003) 4.790 2.769 ( -1.746 -1.008 -1.150) 2.321 2.926 ( 0.858 0.495 0.175) 1.006 3.433 ( -0.615 -0.355 -7.191) 7.226 3.835 ( -3.761 -2.171 1.136) 4.488 ======================= Grid point 184 (25/126) ======================= q-point: ( 0.23 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.555 ( 8.144 4.702 2.295) 9.679 0.780 ( 8.389 4.843 0.985) 9.736 1.237 ( 13.958 8.058 -0.393) 16.122 2.211 ( -0.424 -0.245 -0.426) 0.649 2.464 ( 3.369 1.945 0.012) 3.890 2.729 ( -1.451 -0.838 -0.828) 1.869 2.944 ( 0.654 0.378 0.134) 0.767 3.330 ( -6.621 -3.823 -1.302) 7.756 3.788 ( -0.985 -0.569 -0.623) 1.297 ======================= Grid point 185 (26/126) ======================= q-point: ( 0.31 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.743 ( 7.552 4.360 1.198) 8.802 0.953 ( 6.307 3.641 0.550) 7.303 1.518 ( 9.954 5.747 -0.567) 11.507 2.203 ( -0.290 -0.168 -0.066) 0.342 2.510 ( 0.251 0.145 0.190) 0.346 2.710 ( -0.033 -0.019 -0.515) 0.516 2.956 ( 0.411 0.237 -0.009) 0.474 3.182 ( -5.212 -3.009 0.038) 6.019 3.766 ( -1.131 -0.653 -0.332) 1.348 ======================= Grid point 186 (27/126) ======================= q-point: ( 0.38 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.891 ( 4.733 2.733 0.633) 5.502 1.071 ( 3.664 2.116 0.164) 4.234 1.697 ( 5.248 3.030 -0.118) 6.061 2.196 ( -0.285 -0.165 0.304) 0.448 2.483 ( -1.940 -1.120 -0.104) 2.243 2.728 ( 1.453 0.839 -0.243) 1.696 2.964 ( 0.315 0.182 -0.159) 0.397 3.100 ( -1.980 -1.143 0.295) 2.305 3.733 ( -1.610 -0.929 -0.203) 1.870 ======================= Grid point 187 (28/126) ======================= q-point: ( 0.46 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.953 ( 0.504 0.291 -0.255) 0.636 1.122 ( 0.700 0.404 -0.364) 0.886 1.765 ( 0.733 0.423 0.239) 0.880 2.192 ( -0.040 -0.023 0.560) 0.562 2.451 ( -0.282 -0.163 -0.362) 0.487 2.766 ( 1.441 0.832 0.103) 1.667 2.970 ( 0.184 0.106 -0.202) 0.293 3.073 ( -0.630 -0.364 0.279) 0.780 3.697 ( -1.311 -0.757 -0.527) 1.603 ======================= Grid point 188 (29/126) ======================= q-point: (-0.46 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.917 ( -3.387 -1.956 -1.552) 4.208 1.106 ( -2.011 -1.161 -0.945) 2.507 1.741 ( -2.563 -1.480 -0.622) 3.024 2.196 ( 0.411 0.237 0.709) 0.854 2.476 ( 2.147 1.240 0.524) 2.534 2.782 ( -0.221 -0.128 0.661) 0.709 2.971 ( -0.144 -0.083 -0.175) 0.241 3.065 ( 0.120 0.069 0.009) 0.139 3.678 ( -0.231 -0.133 -1.189) 1.218 ======================= Grid point 189 (30/126) ======================= q-point: (-0.38 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.808 ( -5.681 -3.280 -2.402) 6.986 1.033 ( -4.162 -2.403 -1.199) 4.953 1.650 ( -5.273 -3.044 -1.744) 6.333 2.209 ( 0.629 0.363 0.751) 1.045 2.529 ( 1.825 1.053 1.470) 2.569 2.762 ( -1.185 -0.684 0.967) 1.675 2.964 ( -0.488 -0.282 -0.187) 0.593 3.092 ( 2.631 1.519 0.041) 3.039 3.682 ( 0.327 0.189 -2.051) 2.085 ======================= Grid point 190 (31/126) ======================= q-point: (-0.31 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.662 ( -6.565 -3.791 -1.761) 7.783 0.912 ( -6.119 -3.533 -0.905) 7.123 1.487 ( -8.810 -5.086 -1.678) 10.310 2.223 ( 0.524 0.302 0.638) 0.879 2.542 ( -0.820 -0.473 1.178) 1.511 2.739 ( -0.617 -0.356 0.518) 0.881 2.950 ( -0.624 -0.360 -0.220) 0.754 3.208 ( 7.326 4.230 1.197) 8.543 3.682 ( -0.499 -0.288 -3.451) 3.498 ======================= Grid point 191 (32/126) ======================= q-point: (-0.23 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.506 ( -6.665 -3.848 0.615) 7.720 0.745 ( -8.195 -4.731 -0.257) 9.466 1.234 ( -12.817 -7.400 -0.464) 14.807 2.234 ( 0.403 0.232 0.372) 0.595 2.494 ( -3.159 -1.824 0.704) 3.716 2.740 ( 0.768 0.443 -0.400) 0.973 2.935 ( -0.631 -0.365 -0.159) 0.746 3.422 ( 10.464 6.041 2.941) 12.435 3.650 ( -2.297 -1.326 -6.311) 6.845 ======================= Grid point 192 (33/126) ======================= q-point: (-0.15 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.358 ( -5.755 -3.323 4.211) 7.867 0.531 ( -10.029 -5.790 0.506) 11.591 0.893 ( -16.115 -9.304 0.995) 18.635 2.243 ( 0.406 0.235 0.052) 0.472 2.403 ( -4.328 -2.499 0.568) 5.030 2.770 ( 1.587 0.916 -1.085) 2.129 2.920 ( -0.652 -0.376 -0.001) 0.752 3.540 ( -5.595 -3.230 -9.848) 11.778 3.697 ( 10.598 6.119 2.843) 12.563 ======================= Grid point 193 (34/126) ======================= q-point: (-0.08 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.250 ( -3.250 -1.877 8.214) 9.031 0.286 ( -10.440 -6.028 1.885) 12.202 0.493 ( -17.544 -10.129 2.725) 20.441 2.252 ( 0.279 0.161 -0.216) 0.388 2.305 ( -3.716 -2.145 -0.081) 4.291 2.804 ( 1.140 0.658 -1.242) 1.810 2.906 ( -0.522 -0.301 0.132) 0.617 3.296 ( -16.654 -9.615 -17.494) 25.997 3.944 ( 9.668 5.582 4.161) 11.913 ======================= Grid point 195 (35/126) ======================= q-point: ( 0.08 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.325 ( 3.540 4.657 5.799) 8.237 0.504 ( 8.104 10.820 1.296) 13.580 0.768 ( 8.452 16.848 0.927) 18.872 2.257 ( -0.284 0.715 -0.675) 1.023 2.337 ( 0.403 3.444 0.639) 3.525 2.779 ( -0.967 -1.651 -1.192) 2.254 2.923 ( 0.827 1.128 0.116) 1.403 3.463 ( 0.414 6.775 -11.132) 13.038 3.814 ( -0.926 -9.981 2.880) 10.429 ======================= Grid point 196 (36/126) ======================= q-point: ( 0.15 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.465 ( 6.663 4.783 3.124) 8.777 0.743 ( 6.337 10.810 0.994) 12.570 1.095 ( 11.810 11.293 -0.268) 16.343 2.248 ( -2.057 2.238 -0.621) 3.103 2.400 ( 3.891 0.425 0.457) 3.941 2.742 ( -1.309 -1.544 -0.927) 2.226 2.952 ( 0.511 1.977 0.103) 2.045 3.424 ( -6.763 -1.088 -2.580) 7.320 3.747 ( 1.980 -5.655 -0.791) 6.043 ======================= Grid point 197 (37/126) ======================= q-point: ( 0.23 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.643 ( 8.026 4.102 1.610) 9.156 0.940 ( 3.955 10.119 0.565) 10.879 1.392 ( 11.086 6.980 -0.680) 13.118 2.236 ( -2.017 2.447 -0.283) 3.184 2.460 ( 3.431 -2.158 0.306) 4.065 2.710 ( -0.645 -0.984 -0.510) 1.282 2.979 ( -0.204 2.802 0.005) 2.809 3.270 ( -7.619 -2.137 -0.187) 7.916 3.734 ( 1.740 -4.337 -0.773) 4.737 ======================= Grid point 198 (38/126) ======================= q-point: ( 0.31 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.812 ( 6.886 2.660 1.014) 7.451 1.084 ( 1.513 8.571 0.253) 8.708 1.614 ( 7.832 3.664 -0.421) 8.657 2.227 ( -1.820 2.277 0.083) 2.916 2.468 ( 0.903 -3.586 0.051) 3.699 2.705 ( 0.955 -0.361 -0.248) 1.051 3.001 ( -0.946 3.530 -0.125) 3.656 3.142 ( -4.755 -0.967 0.260) 4.860 3.713 ( 0.834 -4.268 -0.290) 4.359 ======================= Grid point 199 (39/126) ======================= q-point: ( 0.38 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.922 ( 3.402 0.570 0.436) 3.477 1.165 ( -1.093 6.773 -0.002) 6.861 1.736 ( 3.379 0.886 0.208) 3.499 2.220 ( -1.592 2.345 0.309) 2.851 2.437 ( 0.111 -2.810 -0.450) 2.848 2.731 ( 2.134 -0.409 -0.147) 2.178 3.016 ( -1.733 4.108 -0.145) 4.461 3.082 ( -1.593 -0.610 0.372) 1.746 3.680 ( 0.408 -4.031 -0.290) 4.062 ======================= Grid point 200 (40/126) ======================= q-point: ( 0.46 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.937 ( -1.021 -1.651 -0.630) 2.041 1.184 ( -3.799 5.452 -0.358) 6.655 1.756 ( -0.624 -1.265 0.012) 1.410 2.218 ( -1.320 2.401 0.389) 2.767 2.435 ( 2.119 -1.190 -0.142) 2.435 2.758 ( 1.547 -1.311 0.131) 2.032 3.025 ( -2.316 4.640 0.030) 5.186 3.061 ( -0.131 -0.791 0.179) 0.821 3.653 ( 0.963 -3.379 -0.734) 3.590 ======================= Grid point 202 (41/126) ======================= q-point: (-0.46 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.863 ( -4.353 -3.060 -1.689) 5.582 1.144 ( -6.466 4.644 -0.765) 7.997 1.698 ( -3.514 -2.726 -1.135) 4.591 2.224 ( -0.716 2.241 0.472) 2.400 2.476 ( 3.298 -1.250 0.881) 3.635 2.754 ( 0.049 -2.082 0.633) 2.176 3.025 ( -2.942 4.009 0.176) 4.976 3.069 ( 1.491 1.082 -0.075) 1.843 3.648 ( 1.865 -2.790 -1.461) 3.660 ======================= Grid point 203 (42/126) ======================= q-point: (-0.38 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.736 ( -6.080 -3.331 -1.772) 7.156 1.053 ( -8.827 3.852 -0.839) 9.667 1.576 ( -6.587 -4.158 -1.783) 7.991 2.237 ( -0.448 2.318 0.497) 2.412 2.504 ( 2.133 -3.341 1.226) 4.149 2.730 ( -0.119 -1.796 0.619) 1.904 3.005 ( -2.996 2.638 0.002) 3.992 3.149 ( 4.460 4.709 0.626) 6.516 3.654 ( 1.787 -2.739 -2.461) 4.093 ======================= Grid point 204 (43/126) ======================= q-point: (-0.31 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.590 ( -6.859 -2.786 -0.284) 7.409 0.912 ( -10.981 2.713 -0.421) 11.319 1.378 ( -10.621 -5.259 -1.216) 11.914 2.252 ( -0.704 2.410 0.295) 2.528 2.481 ( 0.130 -4.885 0.870) 4.964 2.722 ( 0.866 -0.736 -0.040) 1.137 2.977 ( -2.665 1.783 -0.076) 3.207 3.318 ( 8.149 5.939 1.992) 10.279 3.647 ( 0.413 -2.513 -4.229) 4.936 ======================= Grid point 205 (44/126) ======================= q-point: (-0.23 0.08 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.445 ( -6.982 -1.508 2.493) 7.566 0.724 ( -12.813 1.482 0.295) 12.902 1.100 ( -15.007 -4.751 -0.062) 15.741 2.260 ( -0.821 1.857 -0.163) 2.037 2.416 ( -1.841 -4.544 0.766) 4.963 2.744 ( 1.820 0.074 -0.814) 1.995 2.949 ( -2.050 0.873 -0.007) 2.229 3.510 ( 6.667 -0.750 0.455) 6.725 3.640 ( 1.527 3.438 -4.507) 5.870 ======================= Grid point 210 (45/126) ======================= q-point: ( 0.15 0.15 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.573 ( 4.999 5.767 1.312) 7.744 0.964 ( 6.889 10.367 0.315) 12.452 1.314 ( 6.165 10.081 -0.954) 11.856 2.293 ( 0.421 1.552 -0.537) 1.696 2.400 ( -0.646 0.065 0.984) 1.178 2.713 ( -0.939 -1.089 -0.401) 1.493 2.996 ( 1.307 2.211 0.008) 2.569 3.370 ( -5.463 -4.499 0.861) 7.129 3.644 ( 1.874 -3.785 -2.668) 4.996 ======================= Grid point 211 (46/126) ======================= q-point: ( 0.23 0.15 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.725 ( 6.242 4.063 0.899) 7.502 1.157 ( 2.734 10.632 0.227) 10.981 1.514 ( 7.565 5.005 -0.870) 9.113 2.298 ( -2.480 3.304 -0.216) 4.136 2.398 ( 2.068 -3.377 0.404) 3.980 2.695 ( -0.188 -0.390 -0.054) 0.436 3.033 ( 0.401 1.995 -0.036) 2.035 3.241 ( -5.887 -0.670 0.463) 5.943 3.631 ( 1.975 -5.409 -1.013) 5.847 ======================= Grid point 212 (47/126) ======================= q-point: ( 0.31 0.15 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.858 ( 4.904 1.934 0.829) 5.337 1.275 ( -0.970 9.778 0.183) 9.827 1.667 ( 5.669 1.530 -0.058) 5.872 2.285 ( -2.379 3.240 0.087) 4.021 2.394 ( 1.211 -3.102 -0.506) 3.368 2.699 ( 0.948 -0.265 -0.080) 0.987 3.049 ( 0.066 0.192 -0.004) 0.203 3.159 ( -4.710 3.078 0.328) 5.636 3.606 ( 1.863 -5.875 -0.318) 6.171 ======================= Grid point 213 (48/126) ======================= q-point: ( 0.38 0.15 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.923 ( 1.343 -0.254 0.412) 1.427 1.326 ( -3.655 8.635 0.235) 9.379 1.735 ( 1.773 -0.909 0.397) 2.031 2.282 ( -2.140 3.432 -0.102) 4.046 2.388 ( 1.377 -1.780 -0.689) 2.353 2.716 ( 1.546 -1.009 -0.215) 1.859 3.044 ( 0.268 -1.330 0.093) 1.359 3.134 ( -3.396 5.296 0.325) 6.300 3.581 ( 1.927 -5.464 -0.283) 5.801 ======================= Grid point 214 (49/126) ======================= q-point: (-0.54 0.15 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.900 ( -2.606 -1.643 -0.381) 3.105 1.324 ( -6.111 7.749 0.132) 9.870 1.717 ( -1.609 -2.486 -0.260) 2.972 2.275 ( -1.996 2.799 -0.273) 3.449 2.411 ( 2.416 -1.103 0.070) 2.657 2.725 ( 0.905 -1.694 0.083) 1.923 3.040 ( 0.869 -1.084 -0.019) 1.390 3.137 ( -2.762 5.938 0.457) 6.565 3.568 ( 2.572 -4.755 -0.750) 5.457 ======================= Grid point 215 (50/126) ======================= q-point: (-0.46 0.15 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.811 ( -5.266 -1.625 -0.819) 5.571 1.267 ( -8.715 6.486 -0.111) 10.865 1.635 ( -4.392 -3.284 -1.213) 5.616 2.278 ( -0.903 2.927 -0.059) 3.064 2.431 ( 1.886 -3.018 0.645) 3.617 2.715 ( -0.039 -1.564 0.472) 1.634 3.048 ( 0.247 -0.132 -0.039) 0.283 3.167 ( -0.293 6.192 0.671) 6.236 3.574 ( 3.083 -3.950 -1.510) 5.233 ======================= Grid point 216 (51/126) ======================= q-point: (-0.38 0.15 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.694 ( -6.706 -0.421 -0.238) 6.724 1.147 ( -10.972 3.983 -0.104) 11.673 1.498 ( -7.701 -2.614 -1.450) 8.261 2.297 ( -0.887 3.258 -0.103) 3.378 2.416 ( 1.210 -4.605 0.873) 4.841 2.703 ( 0.338 -0.835 0.232) 0.930 3.035 ( -1.931 0.519 0.002) 1.999 3.258 ( 4.344 3.928 1.236) 5.985 3.595 ( 2.682 -1.746 -2.426) 4.016 ======================= Grid point 224 (52/126) ======================= q-point: ( 0.23 0.23 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.807 ( 3.824 3.984 0.332) 5.532 1.350 ( 4.579 7.916 0.297) 9.150 1.599 ( 3.314 3.564 -0.752) 4.924 2.338 ( 2.377 -1.712 -0.459) 2.965 2.361 ( -3.030 2.129 0.261) 3.712 2.694 ( -0.378 0.202 0.205) 0.476 3.047 ( -0.291 -0.646 0.003) 0.709 3.240 ( -0.490 0.331 0.714) 0.927 3.532 ( -0.869 -3.927 -1.041) 4.155 ======================= Grid point 225 (53/126) ======================= q-point: ( 0.31 0.23 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.894 ( 2.445 1.703 0.786) 3.081 1.460 ( -0.937 8.103 0.461) 8.170 1.679 ( 3.408 -0.197 0.330) 3.429 2.342 ( -0.160 1.365 -0.017) 1.374 2.359 ( -0.259 0.400 -1.251) 1.339 2.691 ( 0.263 -0.562 -0.159) 0.640 3.026 ( -0.259 -1.807 0.066) 1.827 3.239 ( -2.284 4.056 0.421) 4.674 3.490 ( 1.351 -5.301 -0.174) 5.474 ======================= Grid point 226 (54/126) ======================= q-point: ( 0.38 0.23 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.920 ( -0.572 0.025 0.688) 0.895 1.494 ( -4.013 7.334 0.645) 8.384 1.702 ( 0.612 -2.266 0.443) 2.389 2.339 ( -2.298 1.637 -1.045) 3.009 2.374 ( 1.495 0.462 -0.490) 1.640 2.693 ( 0.916 -1.171 -0.300) 1.516 3.011 ( 0.723 -1.535 0.051) 1.698 3.242 ( -2.644 5.039 0.486) 5.712 3.471 ( 2.514 -5.118 -0.090) 5.703 ======================= Grid point 228 (55/126) ======================= q-point: (-0.54 0.23 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.881 ( -3.469 -0.112 0.325) 3.486 1.465 ( -6.688 4.890 0.625) 8.308 1.664 ( -2.052 -2.082 -0.255) 2.934 2.320 ( -1.213 1.648 -0.773) 2.188 2.390 ( 0.355 -1.195 -0.157) 1.256 2.697 ( 0.621 -0.969 0.020) 1.152 3.023 ( 1.643 -0.446 -0.029) 1.703 3.245 ( -2.173 3.727 0.674) 4.366 3.479 ( 3.370 -3.115 -0.589) 4.627 ======================= Grid point 239 (56/126) ======================= q-point: ( 0.31 0.31 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.923 ( 0.665 0.890 0.977) 1.479 1.588 ( 2.059 4.106 0.791) 4.661 1.663 ( -0.264 -1.217 0.580) 1.374 2.357 ( 1.300 1.060 -1.128) 2.021 2.379 ( -0.746 0.102 -0.871) 1.151 2.678 ( -0.374 -0.533 -0.358) 0.743 2.994 ( -0.505 -0.968 0.072) 1.094 3.312 ( 1.489 2.950 0.479) 3.340 3.400 ( -1.659 -3.320 0.134) 3.714 ======================= Grid point 362 (57/126) ======================= q-point: ( 0.00 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.234 ( -0.000 0.000 4.219) 4.219 0.234 ( 0.000 0.000 4.219) 4.219 0.437 ( 0.000 0.000 8.033) 8.033 2.243 ( 0.000 -0.000 -0.516) 0.516 2.243 ( 0.000 0.000 -0.516) 0.516 2.790 ( 0.000 -0.000 -1.486) 1.486 2.903 ( 0.000 -0.000 0.196) 0.196 2.903 ( 0.000 0.000 0.196) 0.196 4.089 ( 0.000 -0.000 0.229) 0.229 ======================= Grid point 363 (58/126) ======================= q-point: ( 0.08 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.380 ( 8.865 5.118 2.316) 10.496 0.410 ( 4.543 2.623 6.447) 8.312 0.519 ( 12.648 7.303 2.960) 14.902 2.231 ( -0.687 -0.397 -0.420) 0.898 2.287 ( 3.040 1.755 -0.305) 3.524 2.773 ( -1.174 -0.678 -1.632) 2.122 2.911 ( 0.613 0.354 0.165) 0.727 3.090 ( 12.024 6.942 -4.780) 14.684 4.002 ( -6.493 -3.749 0.748) 7.535 ======================= Grid point 364 (59/126) ======================= q-point: ( 0.15 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.481 ( 3.614 2.087 4.709) 6.292 0.595 ( 9.066 5.234 1.093) 10.526 0.883 ( 15.746 9.091 0.552) 18.190 2.215 ( -0.561 -0.324 -0.221) 0.684 2.375 ( 4.257 2.458 0.223) 4.920 2.741 ( -1.437 -0.830 -1.441) 2.198 2.928 ( 0.802 0.463 0.069) 0.929 3.314 ( 5.461 3.153 -4.229) 7.592 3.837 ( -6.462 -3.731 -0.546) 7.482 ======================= Grid point 365 (60/126) ======================= q-point: ( 0.23 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.595 ( 6.018 3.475 1.412) 7.091 0.793 ( 7.774 4.488 0.230) 8.980 1.225 ( 13.210 7.627 -0.710) 15.270 2.208 ( -0.114 -0.066 0.110) 0.171 2.472 ( 3.573 2.063 0.550) 4.162 2.712 ( -0.841 -0.485 -0.714) 1.205 2.944 ( 0.549 0.317 -0.065) 0.637 3.318 ( -4.334 -2.502 0.240) 5.010 3.747 ( -1.261 -0.728 -3.068) 3.396 ======================= Grid point 366 (61/126) ======================= q-point: ( 0.31 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.747 ( 6.586 3.803 -0.816) 7.649 0.954 ( 5.864 3.385 -0.476) 6.788 1.495 ( 9.787 5.651 -1.593) 11.413 2.208 ( 0.055 0.032 0.491) 0.495 2.524 ( 0.505 0.291 1.043) 1.195 2.707 ( 0.433 0.250 0.202) 0.539 2.954 ( 0.294 0.170 -0.175) 0.382 3.190 ( -5.215 -3.011 0.718) 6.064 3.735 ( -0.493 -0.285 -2.526) 2.590 ======================= Grid point 367 (62/126) ======================= q-point: ( 0.38 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.881 ( 4.369 2.523 -1.551) 5.278 1.063 ( 3.361 1.940 -0.882) 3.980 1.680 ( 5.916 3.416 -1.442) 6.982 2.208 ( -0.062 -0.036 0.749) 0.753 2.494 ( -2.456 -1.418 1.136) 3.055 2.733 ( 1.670 0.964 0.657) 2.038 2.960 ( 0.256 0.148 -0.223) 0.370 3.105 ( -2.055 -1.187 0.184) 2.381 3.713 ( -1.360 -0.785 -1.671) 2.293 ======================= Grid point 368 (63/126) ======================= q-point: ( 0.46 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.936 ( 0.159 0.092 -1.339) 1.352 1.108 ( 0.475 0.274 -0.916) 1.068 1.764 ( 1.108 0.640 -0.354) 1.327 2.207 ( 0.046 0.027 0.816) 0.818 2.449 ( -0.579 -0.334 0.210) 0.701 2.774 ( 1.549 0.894 0.598) 1.886 2.966 ( 0.184 0.106 -0.235) 0.316 3.077 ( -0.660 -0.381 0.130) 0.773 3.678 ( -1.389 -0.802 -1.178) 1.990 ======================= Grid point 369 (64/126) ======================= q-point: (-0.46 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.890 ( -3.819 -2.205 -0.908) 4.502 1.088 ( -2.114 -1.220 -0.712) 2.543 1.734 ( -3.376 -1.949 -0.063) 3.899 2.212 ( 0.407 0.235 0.756) 0.890 2.481 ( 2.890 1.669 -0.088) 3.339 2.793 ( -0.007 -0.004 0.405) 0.405 2.967 ( -0.136 -0.079 -0.245) 0.291 3.068 ( -0.014 -0.008 0.201) 0.202 3.655 ( -0.548 -0.317 -0.948) 1.140 ======================= Grid point 370 (65/126) ======================= q-point: (-0.38 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.774 ( -5.654 -3.264 -0.619) 6.558 1.015 ( -3.989 -2.303 -0.407) 4.625 1.622 ( -5.955 -3.438 -0.833) 6.926 2.224 ( 0.499 0.288 0.582) 0.819 2.550 ( 2.213 1.278 0.495) 2.603 2.775 ( -1.349 -0.779 0.253) 1.578 2.959 ( -0.475 -0.274 -0.217) 0.590 3.096 ( 3.042 1.756 0.299) 3.526 3.648 ( -0.322 -0.186 -1.070) 1.133 ======================= Grid point 371 (66/126) ======================= q-point: (-0.31 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.644 ( -5.234 -3.022 0.083) 6.044 0.903 ( -5.627 -3.249 0.027) 6.497 1.456 ( -8.418 -4.860 -1.253) 9.801 2.233 ( 0.242 0.140 0.243) 0.371 2.560 ( -1.385 -0.799 0.492) 1.673 2.742 ( -1.224 -0.707 -0.189) 1.426 2.947 ( -0.556 -0.321 -0.105) 0.651 3.230 ( 8.363 4.829 0.772) 9.688 3.625 ( -1.917 -1.107 -1.723) 2.805 ======================= Grid point 372 (67/126) ======================= q-point: (-0.23 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.536 ( -3.920 -2.263 2.147) 5.010 0.750 ( -7.505 -4.333 0.674) 8.692 1.221 ( -11.715 -6.764 -0.747) 13.548 2.236 ( 0.078 0.045 -0.173) 0.195 2.499 ( -3.584 -2.069 -0.192) 4.143 2.727 ( 0.024 0.014 -0.807) 0.808 2.934 ( -0.509 -0.294 0.046) 0.590 3.464 ( 10.029 5.790 -0.716) 11.602 3.549 ( -3.706 -2.140 -1.352) 4.488 ======================= Grid point 373 (68/126) ======================= q-point: (-0.15 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.461 ( -2.494 -1.440 5.278) 6.012 0.553 ( -9.115 -5.263 1.578) 10.643 0.907 ( -14.951 -8.632 0.227) 17.266 2.239 ( 0.180 0.104 -0.462) 0.507 2.400 ( -4.664 -2.693 -0.733) 5.435 2.743 ( 1.191 0.687 -1.414) 1.972 2.922 ( -0.543 -0.314 0.147) 0.644 3.364 ( -9.758 -5.634 -6.521) 13.019 3.753 ( 11.598 6.696 2.370) 13.600 ======================= Grid point 374 (69/126) ======================= q-point: (-0.08 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.340 ( -8.415 -4.858 2.989) 10.166 0.413 ( -2.744 -1.584 7.212) 7.877 0.545 ( -14.265 -8.236 1.960) 16.588 2.243 ( 0.169 0.098 -0.548) 0.581 2.293 ( -4.033 -2.328 -0.865) 4.736 2.773 ( 1.165 0.673 -1.623) 2.108 2.909 ( -0.475 -0.274 0.186) 0.579 3.090 ( -12.490 -7.211 -4.894) 15.230 3.991 ( 7.931 4.579 1.236) 9.242 ======================= Grid point 376 (70/126) ======================= q-point: ( 0.08 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.451 ( 1.829 2.063 5.952) 6.559 0.535 ( 7.871 9.740 1.608) 12.625 0.787 ( 7.458 15.639 0.661) 17.339 2.242 ( -0.580 0.795 -0.684) 1.198 2.344 ( 1.148 3.720 -0.020) 3.893 2.750 ( -0.865 -1.418 -1.543) 2.267 2.926 ( 0.814 1.049 0.127) 1.334 3.285 ( 4.593 10.310 -5.995) 12.780 3.851 ( -3.312 -10.501 1.117) 11.067 ======================= Grid point 377 (71/126) ======================= q-point: ( 0.15 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.531 ( 4.511 2.959 3.002) 6.174 0.761 ( 5.912 10.237 0.663) 11.840 1.089 ( 11.181 10.388 -0.378) 15.266 2.237 ( -1.498 2.057 -0.403) 2.576 2.412 ( 3.578 0.998 0.581) 3.759 2.720 ( -0.911 -1.048 -1.051) 1.741 2.953 ( 0.457 1.885 0.010) 1.939 3.381 ( -2.676 2.826 -1.498) 4.170 3.716 ( 0.002 -7.484 -1.897) 7.720 ======================= Grid point 378 (72/126) ======================= q-point: ( 0.23 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.664 ( 6.615 3.405 0.347) 7.448 0.945 ( 3.506 9.817 -0.062) 10.424 1.371 ( 10.840 6.550 -1.309) 12.733 2.234 ( -1.475 2.370 0.055) 2.792 2.473 ( 3.342 -1.865 0.839) 3.918 2.702 ( -0.161 -0.450 -0.178) 0.510 2.978 ( -0.292 2.739 -0.098) 2.756 3.278 ( -7.373 -1.493 0.979) 7.586 3.694 ( 2.548 -4.153 -2.933) 5.687 ======================= Grid point 379 (73/126) ======================= q-point: ( 0.31 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.813 ( 6.156 2.732 -0.981) 6.806 1.081 ( 1.194 8.590 -0.510) 8.688 1.593 ( 8.153 3.963 -1.583) 9.202 2.232 ( -1.387 2.119 0.392) 2.562 2.480 ( 0.792 -3.959 1.035) 4.168 2.707 ( 1.116 -0.147 0.468) 1.219 2.999 ( -0.992 3.585 -0.095) 3.721 3.150 ( -5.005 -0.897 0.515) 5.111 3.689 ( 1.344 -3.958 -1.951) 4.613 ======================= Grid point 380 (74/126) ======================= q-point: ( 0.38 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.915 ( 2.820 0.923 -1.142) 3.180 1.159 ( -1.442 7.146 -0.586) 7.313 1.729 ( 4.022 1.464 -0.824) 4.359 2.229 ( -1.259 1.969 0.526) 2.396 2.438 ( -0.292 -3.457 0.513) 3.507 2.734 ( 2.183 -0.643 0.443) 2.318 3.015 ( -1.723 4.326 0.039) 4.657 3.087 ( -1.734 -0.506 0.147) 1.813 3.663 ( 0.444 -3.817 -1.274) 4.048 ======================= Grid point 381 (75/126) ======================= q-point: ( 0.46 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.922 ( -1.738 -1.227 -0.748) 2.255 1.175 ( -4.309 5.993 -0.483) 7.397 1.757 ( -0.971 -1.298 0.043) 1.622 2.228 ( -0.989 1.894 0.545) 2.206 2.431 ( 2.541 -1.282 -0.270) 2.859 2.761 ( 1.887 -1.628 0.153) 2.497 3.027 ( -2.401 5.066 0.138) 5.608 3.065 ( -0.289 -0.850 0.135) 0.908 3.635 ( 0.645 -3.265 -0.913) 3.451 ======================= Grid point 383 (76/126) ======================= q-point: (-0.46 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.840 ( -5.023 -2.448 -0.405) 5.602 1.132 ( -6.876 5.192 -0.350) 8.623 1.682 ( -4.540 -3.048 -0.341) 5.479 2.234 ( -0.453 1.791 0.505) 1.915 2.485 ( 4.202 -1.295 -0.018) 4.397 2.762 ( 0.456 -2.250 0.141) 2.300 3.026 ( -3.214 4.108 -0.035) 5.216 3.072 ( 1.529 1.449 0.301) 2.128 3.622 ( 1.274 -2.853 -0.883) 3.247 ======================= Grid point 384 (77/126) ======================= q-point: (-0.38 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.715 ( -5.954 -2.253 -0.108) 6.367 1.042 ( -8.827 4.390 -0.125) 9.859 1.545 ( -6.937 -4.152 -1.113) 8.161 2.246 ( -0.435 1.872 0.301) 1.946 2.522 ( 2.367 -3.698 0.452) 4.414 2.738 ( -0.388 -2.022 0.154) 2.065 3.003 ( -3.124 2.685 -0.149) 4.122 3.163 ( 5.101 5.443 0.612) 7.485 3.612 ( 0.718 -3.145 -1.336) 3.492 ======================= Grid point 385 (78/126) ======================= q-point: (-0.31 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.598 ( -5.224 -1.268 0.996) 5.467 0.911 ( -10.579 3.235 0.302) 11.067 1.352 ( -9.863 -5.026 -1.190) 11.134 2.254 ( -0.858 1.823 -0.107) 2.018 2.494 ( -0.287 -4.807 0.398) 4.832 2.718 ( 0.145 -1.080 -0.336) 1.140 2.975 ( -2.645 1.876 -0.084) 3.244 3.355 ( 9.217 6.902 1.262) 11.584 3.571 ( -1.813 -3.680 -2.610) 4.862 ======================= Grid point 386 (79/126) ======================= q-point: (-0.23 0.08 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.508 ( -4.045 -0.266 3.370) 5.271 0.737 ( -12.164 1.936 0.916) 12.351 1.095 ( -13.695 -4.692 -0.469) 14.484 2.253 ( -0.922 1.052 -0.530) 1.496 2.425 ( -2.343 -3.835 0.125) 4.495 2.724 ( 1.156 -0.163 -1.054) 1.573 2.949 ( -1.950 0.953 0.046) 2.170 3.465 ( -5.680 -4.180 -3.950) 8.083 3.608 ( 11.066 6.344 0.931) 12.790 ======================= Grid point 391 (80/126) ======================= q-point: ( 0.15 0.15 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.602 ( 4.392 4.046 1.326) 6.117 0.972 ( 6.865 9.953 0.437) 12.099 1.290 ( 5.886 9.390 -1.265) 11.154 2.280 ( 0.281 1.721 -0.609) 1.847 2.421 ( -0.458 0.232 0.976) 1.103 2.704 ( -0.616 -0.368 -0.386) 0.815 2.996 ( 1.293 2.154 -0.052) 2.513 3.397 ( -6.874 -2.461 1.795) 7.518 3.577 ( 4.220 -4.559 -3.640) 7.200 ======================= Grid point 392 (81/126) ======================= q-point: ( 0.23 0.15 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.735 ( 5.889 3.669 0.002) 6.938 1.160 ( 2.470 10.733 0.097) 11.014 1.488 ( 7.976 5.053 -1.512) 9.563 2.291 ( -1.737 2.920 -0.378) 3.418 2.414 ( 1.277 -3.386 1.002) 3.755 2.697 ( -0.011 -0.017 0.277) 0.277 3.032 ( 0.367 2.037 -0.067) 2.071 3.258 ( -6.605 -0.401 1.155) 6.717 3.595 ( 3.228 -5.124 -2.297) 6.477 ======================= Grid point 393 (82/126) ======================= q-point: ( 0.31 0.15 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.863 ( 4.593 2.305 -0.356) 5.152 1.278 ( -1.127 10.236 0.057) 10.298 1.655 ( 6.344 2.199 -1.004) 6.789 2.284 ( -2.087 2.834 -0.140) 3.523 2.394 ( 0.777 -3.896 0.451) 3.998 2.702 ( 0.763 -0.411 0.356) 0.937 3.049 ( 0.094 0.315 -0.024) 0.330 3.169 ( -4.973 3.170 0.616) 5.929 3.589 ( 2.311 -5.487 -1.298) 6.093 ======================= Grid point 394 (83/126) ======================= q-point: ( 0.38 0.15 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.926 ( 0.759 0.418 -0.112) 0.874 1.331 ( -3.939 9.266 0.209) 10.070 1.739 ( 1.947 -0.485 -0.049) 2.007 2.279 ( -1.846 2.667 -0.122) 3.246 2.377 ( 1.425 -2.250 -0.374) 2.690 2.714 ( 1.502 -1.308 -0.038) 1.992 3.045 ( 0.323 -1.312 0.060) 1.352 3.143 ( -3.485 5.500 0.483) 6.529 3.570 ( 1.786 -5.080 -0.722) 5.433 ======================= Grid point 395 (84/126) ======================= q-point: (-0.54 0.15 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.898 ( -3.402 -0.877 0.196) 3.519 1.328 ( -6.668 8.394 0.223) 10.723 1.716 ( -2.451 -2.584 0.154) 3.565 2.272 ( -1.424 2.143 0.017) 2.574 2.404 ( 3.118 -1.291 -0.690) 3.445 2.724 ( 1.357 -1.778 -0.173) 2.243 3.041 ( 0.806 -1.166 0.117) 1.422 3.147 ( -2.828 6.245 0.462) 6.871 3.554 ( 1.988 -4.555 -0.573) 5.003 ======================= Grid point 396 (85/126) ======================= q-point: (-0.46 0.15 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.805 ( -5.785 -0.751 0.235) 5.838 1.267 ( -9.196 7.115 0.163) 11.628 1.615 ( -5.341 -3.379 -0.622) 6.350 2.279 ( -0.560 2.472 0.108) 2.537 2.435 ( 2.696 -3.315 -0.182) 4.276 2.721 ( 0.260 -1.612 0.157) 1.641 3.047 ( 0.206 -0.127 -0.011) 0.242 3.181 ( -0.113 6.677 0.617) 6.707 3.547 ( 2.162 -4.081 -0.972) 4.720 ======================= Grid point 397 (86/126) ======================= q-point: (-0.38 0.15 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.695 ( -6.217 0.483 0.378) 6.247 1.149 ( -11.039 4.506 0.250) 11.926 1.468 ( -7.849 -2.558 -1.358) 8.366 2.294 ( -0.964 2.400 -0.248) 2.599 2.432 ( 1.631 -4.234 0.624) 4.580 2.707 ( -0.101 -1.001 0.176) 1.021 3.034 ( -1.956 0.550 -0.086) 2.034 3.286 ( 5.108 4.310 1.320) 6.813 3.545 ( 1.162 -1.932 -2.194) 3.146 ======================= Grid point 405 (87/126) ======================= q-point: ( 0.23 0.23 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.814 ( 4.244 4.067 0.293) 5.885 1.358 ( 4.714 8.284 0.512) 9.545 1.580 ( 4.070 4.170 -1.010) 5.914 2.328 ( 1.631 -0.479 -0.472) 1.764 2.366 ( -3.017 -0.018 0.233) 3.026 2.699 ( -0.570 0.085 0.274) 0.638 3.046 ( -0.259 -0.622 -0.018) 0.674 3.259 ( -0.937 0.243 1.024) 1.409 3.506 ( 0.053 -3.294 -1.423) 3.588 ======================= Grid point 406 (88/126) ======================= q-point: ( 0.31 0.23 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.909 ( 2.480 2.223 0.573) 3.380 1.474 ( -1.019 8.618 0.759) 8.712 1.682 ( 4.091 0.449 -0.121) 4.117 2.336 ( -1.522 1.770 -0.723) 2.444 2.337 ( 0.717 -1.301 -0.531) 1.578 2.689 ( -0.024 -0.915 -0.082) 0.919 3.027 ( -0.212 -1.789 0.012) 1.802 3.251 ( -2.515 4.145 0.729) 4.903 3.482 ( 1.779 -4.760 -0.560) 5.112 ======================= Grid point 407 (89/126) ======================= q-point: ( 0.38 0.23 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.936 ( -1.029 0.607 0.853) 1.468 1.510 ( -4.175 7.746 0.832) 8.839 1.713 ( 0.245 -2.040 0.498) 2.114 2.319 ( -1.776 0.966 -0.754) 2.158 2.356 ( 1.671 0.166 -1.136) 2.027 2.685 ( 1.031 -1.315 -0.375) 1.712 3.012 ( 0.702 -1.533 0.069) 1.687 3.254 ( -2.589 5.109 0.612) 5.760 3.468 ( 2.150 -4.715 -0.136) 5.184 ======================= Grid point 409 (90/126) ======================= q-point: (-0.54 0.23 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.893 ( -4.121 0.377 0.752) 4.206 1.479 ( -7.130 5.274 0.749) 8.900 1.661 ( -3.079 -2.107 0.006) 3.731 2.308 ( -0.478 1.463 -0.329) 1.574 2.377 ( 1.032 -1.556 -1.023) 2.129 2.695 ( 1.145 -0.983 -0.162) 1.518 3.023 ( 1.584 -0.440 0.058) 1.645 3.260 ( -2.026 3.888 0.704) 4.440 3.467 ( 2.500 -3.038 -0.468) 3.962 ======================= Grid point 420 (91/126) ======================= q-point: ( 0.31 0.31 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.946 ( 0.857 1.166 1.190) 1.874 1.608 ( 1.950 4.176 1.067) 4.730 1.676 ( 0.187 -0.774 0.651) 1.028 2.332 ( 1.046 0.709 -1.140) 1.702 2.355 ( -0.906 -0.192 -1.342) 1.631 2.670 ( -0.489 -0.718 -0.440) 0.974 2.996 ( -0.473 -0.919 0.076) 1.036 3.324 ( 1.305 2.782 0.610) 3.133 3.403 ( -1.296 -2.857 0.125) 3.140 ======================= Grid point 543 (92/126) ======================= q-point: ( 0.00 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.293 ( 0.000 0.000 1.530) 1.530 0.293 ( 0.000 0.000 1.530) 1.530 0.551 ( 0.000 0.000 2.966) 2.966 2.235 ( 0.000 -0.000 -0.231) 0.231 2.235 ( 0.000 0.000 -0.231) 0.231 2.766 ( 0.000 -0.000 -0.688) 0.688 2.906 ( 0.000 -0.000 0.084) 0.084 2.906 ( 0.000 0.000 0.084) 0.084 4.092 ( 0.000 -0.000 0.099) 0.099 ======================= Grid point 544 (93/126) ======================= q-point: ( 0.08 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.408 ( 7.754 4.477 0.325) 8.960 0.513 ( 10.258 5.922 3.519) 12.357 0.561 ( 4.830 2.788 0.685) 5.619 2.227 ( -0.492 -0.284 0.029) 0.569 2.280 ( 3.401 1.963 -0.314) 3.939 2.747 ( -1.223 -0.706 -0.741) 1.595 2.913 ( 0.546 0.315 0.056) 0.632 3.039 ( 9.557 5.518 -0.969) 11.078 4.010 ( -6.558 -3.786 0.094) 7.573 ======================= Grid point 545 (94/126) ======================= q-point: ( 0.15 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.552 ( 1.176 0.679 1.938) 2.366 0.604 ( 8.432 4.868 -0.250) 9.740 0.893 ( 14.888 8.596 0.390) 17.196 2.217 ( -0.305 -0.176 0.319) 0.475 2.378 ( 4.641 2.680 0.094) 5.360 2.719 ( -1.050 -0.606 -0.587) 1.347 2.928 ( 0.670 0.387 -0.077) 0.778 3.264 ( 8.040 4.642 -0.801) 9.318 3.820 ( -8.673 -5.008 -1.058) 10.071 ======================= Grid point 546 (95/126) ======================= q-point: ( 0.23 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.606 ( 3.540 2.044 -0.370) 4.104 0.788 ( 7.181 4.146 -0.696) 8.321 1.211 ( 12.192 7.039 -0.514) 14.088 2.215 ( 0.153 0.088 0.564) 0.591 2.482 ( 3.953 2.282 0.463) 4.588 2.704 ( -0.063 -0.036 -0.032) 0.079 2.941 ( 0.451 0.260 -0.213) 0.562 3.341 ( -2.618 -1.511 2.026) 3.639 3.675 ( -2.262 -1.306 -3.865) 4.665 ======================= Grid point 547 (96/126) ======================= q-point: ( 0.31 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.713 ( 5.342 3.084 -2.320) 6.590 0.936 ( 5.425 3.132 -1.154) 6.370 1.462 ( 9.258 5.345 -1.464) 10.790 2.221 ( 0.228 0.132 0.676) 0.725 2.546 ( 1.075 0.621 1.014) 1.603 2.717 ( 1.100 0.635 0.737) 1.469 2.950 ( 0.289 0.167 -0.215) 0.397 3.210 ( -6.199 -3.579 1.194) 7.257 3.674 ( 0.554 0.320 -3.128) 3.193 ======================= Grid point 548 (97/126) ======================= q-point: ( 0.38 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.835 ( 4.553 2.629 -2.683) 5.902 1.040 ( 3.345 1.931 -1.272) 4.067 1.646 ( 6.446 3.721 -1.669) 7.628 2.223 ( -0.047 -0.027 0.621) 0.624 2.523 ( -2.726 -1.574 1.475) 3.477 2.752 ( 1.709 0.986 1.078) 2.248 2.956 ( 0.286 0.165 -0.120) 0.351 3.106 ( -2.476 -1.429 -0.057) 2.860 3.673 ( -0.594 -0.343 -1.942) 2.059 ======================= Grid point 549 (98/126) ======================= q-point: ( 0.46 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.906 ( 1.179 0.681 -1.440) 1.982 1.090 ( 0.879 0.508 -0.782) 1.281 1.753 ( 2.249 1.298 -0.709) 2.692 2.221 ( -0.050 -0.029 0.456) 0.459 2.459 ( -1.680 -0.970 0.659) 2.049 2.787 ( 1.114 0.643 0.636) 1.435 2.962 ( 0.168 0.097 -0.089) 0.213 3.077 ( -0.480 -0.277 -0.131) 0.570 3.655 ( -0.762 -0.440 -0.949) 1.294 ======================= Grid point 550 (99/126) ======================= q-point: (-0.46 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.885 ( -2.832 -1.635 0.445) 3.300 1.081 ( -1.611 -0.930 0.059) 1.861 1.739 ( -3.253 -1.878 0.541) 3.795 2.223 ( 0.249 0.144 0.241) 0.375 2.473 ( 2.702 1.560 -0.628) 3.183 2.796 ( -0.421 -0.243 -0.176) 0.517 2.963 ( -0.147 -0.085 -0.128) 0.213 3.073 ( 0.167 0.096 0.269) 0.331 3.645 ( -0.105 -0.061 0.008) 0.122 ======================= Grid point 551 (100/126) ======================= q-point: (-0.38 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.787 ( -5.160 -2.979 1.749) 6.210 1.018 ( -3.628 -2.094 0.700) 4.247 1.621 ( -6.470 -3.735 0.679) 7.501 2.230 ( 0.244 0.141 -0.024) 0.283 2.547 ( 2.589 1.495 -0.779) 3.089 2.771 ( -1.577 -0.911 -0.673) 1.941 2.956 ( -0.384 -0.222 -0.078) 0.451 3.103 ( 2.986 1.724 0.293) 3.461 3.644 ( -0.188 -0.108 0.750) 0.781 ======================= Grid point 552 (101/126) ======================= q-point: (-0.31 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.666 ( -4.687 -2.706 1.997) 5.769 0.913 ( -5.319 -3.071 0.952) 6.215 1.444 ( -8.680 -5.012 0.181) 10.025 2.232 ( -0.094 -0.054 -0.342) 0.359 2.561 ( -1.510 -0.872 -0.389) 1.787 2.733 ( -1.455 -0.840 -0.713) 1.825 2.946 ( -0.404 -0.233 0.092) 0.476 3.231 ( 7.779 4.491 -0.672) 9.007 3.623 ( -1.758 -1.015 1.508) 2.529 ======================= Grid point 553 (102/126) ======================= q-point: (-0.23 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.581 ( -2.501 -1.444 1.964) 3.493 0.769 ( -7.051 -4.071 1.097) 8.216 1.209 ( -11.525 -6.654 -0.320) 13.312 2.227 ( -0.171 -0.099 -0.590) 0.622 2.491 ( -4.054 -2.340 -0.444) 4.702 2.711 ( -0.273 -0.158 -0.646) 0.719 2.937 ( -0.429 -0.247 0.209) 0.537 3.399 ( 2.386 1.378 -3.600) 4.534 3.596 ( 2.828 1.633 3.763) 4.982 ======================= Grid point 554 (103/126) ======================= q-point: (-0.15 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.546 ( -0.654 -0.378 2.578) 2.686 0.586 ( -8.423 -4.863 1.421) 9.830 0.902 ( -14.568 -8.411 -0.549) 16.831 2.227 ( 0.135 0.078 -0.622) 0.641 2.384 ( -4.839 -2.794 -0.616) 5.622 2.720 ( 0.927 0.535 -0.701) 1.280 2.925 ( -0.554 -0.320 0.173) 0.663 3.278 ( -9.401 -5.427 -2.232) 11.082 3.793 ( 10.418 6.015 1.630) 12.140 ======================= Grid point 555 (104/126) ======================= q-point: (-0.08 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 322 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.392 ( -7.491 -4.325 1.854) 8.846 0.524 ( -8.393 -4.846 3.304) 10.239 0.562 ( -6.728 -3.885 -0.262) 7.773 2.232 ( 0.263 0.152 -0.461) 0.553 2.279 ( -3.662 -2.114 -0.349) 4.243 2.747 ( 1.212 0.700 -0.746) 1.586 2.912 ( -0.487 -0.281 0.095) 0.570 3.039 ( -9.647 -5.570 -0.956) 11.180 4.007 ( 7.000 4.041 0.453) 8.095 ======================= Grid point 557 (105/126) ======================= q-point: ( 0.08 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.544 ( 0.446 0.465 2.668) 2.744 0.561 ( 6.433 9.506 0.768) 11.504 0.791 ( 7.819 14.656 -0.097) 16.612 2.232 ( -0.183 0.503 -0.262) 0.596 2.340 ( 1.994 3.822 -0.200) 4.315 2.725 ( -0.728 -1.202 -0.700) 1.570 2.928 ( 0.690 1.040 0.053) 1.249 3.213 ( 5.506 10.138 -1.574) 11.643 3.865 ( -4.883 -9.961 0.338) 11.099 ======================= Grid point 558 (106/126) ======================= q-point: ( 0.15 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.573 ( 2.150 1.542 0.948) 2.810 0.767 ( 5.024 10.260 -0.126) 11.424 1.082 ( 10.716 9.481 -0.206) 14.310 2.234 ( -0.790 1.688 0.129) 1.869 2.420 ( 3.846 1.330 0.234) 4.076 2.705 ( -0.390 -0.609 -0.371) 0.813 2.952 ( 0.309 1.880 -0.076) 1.907 3.369 ( 1.363 5.114 0.406) 5.308 3.677 ( -3.758 -8.617 -1.915) 9.594 ======================= Grid point 559 (107/126) ======================= q-point: ( 0.23 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.655 ( 4.529 2.890 -1.131) 5.490 0.938 ( 2.776 9.972 -0.631) 10.371 1.346 ( 10.162 6.113 -1.003) 11.901 2.239 ( -0.920 2.023 0.401) 2.258 2.489 ( 3.686 -1.628 0.626) 4.078 2.703 ( 0.584 -0.030 0.251) 0.636 2.976 ( -0.357 2.823 -0.118) 2.848 3.309 ( -7.938 -1.542 2.017) 8.334 3.626 ( 3.134 -3.404 -3.539) 5.825 ======================= Grid point 560 (108/126) ======================= q-point: ( 0.31 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.777 ( 5.240 3.231 -2.292) 6.569 1.065 ( 0.691 9.088 -0.914) 9.160 1.558 ( 8.041 4.196 -1.582) 9.207 2.241 ( -1.016 1.702 0.459) 2.035 2.503 ( 1.165 -4.187 1.137) 4.492 2.721 ( 1.473 -0.170 0.797) 1.684 2.998 ( -0.927 3.792 0.017) 3.903 3.162 ( -5.991 -1.125 0.546) 6.120 3.643 ( 2.055 -3.155 -2.294) 4.409 ======================= Grid point 561 (109/126) ======================= q-point: ( 0.38 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.882 ( 3.042 2.015 -1.874) 4.102 1.144 ( -1.627 7.925 -0.764) 8.126 1.707 ( 4.855 2.232 -1.229) 5.484 2.238 ( -1.024 1.488 0.362) 1.842 2.455 ( -0.851 -4.291 1.056) 4.500 2.747 ( 2.031 -1.174 0.684) 2.444 3.018 ( -1.749 4.682 0.229) 5.003 3.087 ( -1.989 -0.440 -0.205) 2.047 3.635 ( 0.960 -3.194 -1.326) 3.589 ======================= Grid point 562 (110/126) ======================= q-point: ( 0.46 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.911 ( -1.023 -0.108 -0.319) 1.077 1.167 ( -4.241 6.671 -0.228) 7.908 1.757 ( -0.308 -0.815 0.029) 0.871 2.237 ( -0.841 1.544 0.265) 1.778 2.426 ( 1.945 -1.965 -0.152) 2.769 2.763 ( 1.672 -2.189 0.061) 2.755 3.029 ( -2.873 5.472 0.105) 6.182 3.066 ( 0.109 -0.864 0.022) 0.872 3.621 ( 1.058 -2.793 -0.405) 3.014 ======================= Grid point 564 (111/126) ======================= q-point: (-0.46 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.848 ( -4.558 -1.821 1.202) 5.053 1.132 ( -6.672 5.528 0.350) 8.672 1.686 ( -4.905 -3.217 0.710) 5.909 2.241 ( -0.505 1.537 0.119) 1.622 2.476 ( 4.583 -1.373 -0.852) 4.860 2.759 ( 0.457 -2.361 -0.435) 2.444 3.023 ( -3.330 4.008 -0.243) 5.217 3.080 ( 1.480 1.668 0.434) 2.272 3.617 ( 1.467 -2.692 0.376) 3.088 ======================= Grid point 565 (112/126) ======================= q-point: (-0.38 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.733 ( -5.680 -2.131 1.767) 6.319 1.048 ( -8.560 4.462 0.668) 9.676 1.537 ( -7.425 -4.477 0.401) 8.679 2.248 ( -0.672 1.594 -0.154) 1.737 2.521 ( 2.751 -3.587 -0.521) 4.551 2.735 ( -0.585 -1.796 -0.452) 1.942 3.000 ( -3.052 2.705 -0.158) 4.081 3.169 ( 4.893 5.404 -0.074) 7.291 3.608 ( 0.723 -3.424 0.909) 3.615 ======================= Grid point 566 (113/126) ======================= q-point: (-0.31 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.626 ( -4.244 -1.216 1.611) 4.699 0.923 ( -10.173 3.201 0.779) 10.693 1.337 ( -9.837 -5.253 -0.174) 11.152 2.248 ( -1.106 1.536 -0.439) 1.943 2.497 ( -0.593 -4.424 -0.153) 4.466 2.710 ( -0.395 -0.950 -0.417) 1.110 2.974 ( -2.451 1.834 0.008) 3.061 3.350 ( 7.382 7.097 -1.667) 10.375 3.565 ( -0.620 -4.615 1.938) 5.044 ======================= Grid point 567 (114/126) ======================= q-point: (-0.23 0.08 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.564 ( -2.005 -0.208 1.809) 2.708 0.756 ( -11.544 1.632 0.865) 11.691 1.084 ( -13.176 -4.591 -0.442) 13.960 2.241 ( -1.021 0.961 -0.512) 1.493 2.424 ( -3.022 -3.053 -0.167) 4.299 2.707 ( 0.637 -0.177 -0.544) 0.856 2.951 ( -1.817 0.853 0.086) 2.009 3.397 ( -7.168 -1.172 -2.395) 7.647 3.639 ( 11.430 2.384 1.820) 11.817 ======================= Grid point 572 (115/126) ======================= q-point: ( 0.15 0.15 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.621 ( 2.802 3.277 0.482) 4.338 0.980 ( 6.267 10.015 0.234) 11.817 1.270 ( 5.481 9.085 -0.611) 10.628 2.271 ( 0.564 1.562 -0.257) 1.680 2.435 ( 0.002 0.355 0.360) 0.505 2.699 ( -0.110 0.125 -0.132) 0.212 2.995 ( 1.246 2.111 -0.043) 2.452 3.442 ( -7.731 -0.857 2.689) 8.230 3.505 ( 4.887 -5.090 -3.380) 7.824 ======================= Grid point 573 (116/126) ======================= q-point: ( 0.23 0.15 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.726 ( 4.567 4.090 -0.826) 6.186 1.160 ( 1.925 11.242 -0.110) 11.406 1.460 ( 7.687 5.255 -1.072) 9.374 2.286 ( -1.156 2.343 -0.073) 2.613 2.432 ( 1.304 -3.427 0.678) 3.729 2.704 ( 0.436 0.098 0.351) 0.568 3.031 ( 0.402 2.063 -0.029) 2.102 3.285 ( -7.512 -0.754 1.287) 7.658 3.547 ( 3.856 -3.978 -2.166) 5.948 ======================= Grid point 574 (117/126) ======================= q-point: ( 0.31 0.15 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.845 ( 4.064 3.473 -1.238) 5.487 1.277 ( -1.476 11.057 -0.141) 11.156 1.631 ( 6.602 2.933 -1.207) 7.325 2.282 ( -1.668 2.163 -0.073) 2.732 2.408 ( 0.644 -4.539 0.853) 4.663 2.711 ( 0.740 -0.816 0.516) 1.217 3.048 ( 0.145 0.238 -0.031) 0.281 3.182 ( -5.508 3.358 0.540) 6.474 3.560 ( 2.751 -4.646 -1.380) 5.573 ======================= Grid point 575 (118/126) ======================= q-point: ( 0.38 0.15 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.919 ( 0.935 1.694 -0.542) 2.009 1.334 ( -4.214 10.035 0.013) 10.884 1.733 ( 2.677 0.305 -0.467) 2.734 2.277 ( -1.341 2.146 -0.073) 2.532 2.375 ( 0.733 -3.236 0.254) 3.328 2.715 ( 1.252 -1.804 0.156) 2.201 3.046 ( 0.482 -1.336 -0.027) 1.421 3.151 ( -3.761 5.733 0.308) 6.864 3.555 ( 2.092 -4.400 -0.629) 4.912 ======================= Grid point 576 (119/126) ======================= q-point: (-0.54 0.15 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.907 ( -3.109 -0.149 0.597) 3.169 1.332 ( -6.870 8.828 0.167) 11.188 1.723 ( -2.450 -2.470 0.491) 3.514 2.275 ( -1.158 1.914 0.180) 2.245 2.387 ( 3.228 -1.803 -0.854) 3.795 2.720 ( 1.361 -1.909 -0.268) 2.360 3.044 ( 0.846 -1.186 0.132) 1.462 3.153 ( -2.984 6.352 0.105) 7.019 3.548 ( 2.079 -4.231 0.051) 4.714 ======================= Grid point 577 (120/126) ======================= q-point: (-0.46 0.15 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.819 ( -5.728 -0.858 1.141) 5.904 1.272 ( -9.230 7.293 0.282) 11.767 1.613 ( -5.955 -3.826 0.473) 7.093 2.281 ( -0.674 2.263 0.074) 2.362 2.425 ( 3.467 -3.200 -0.735) 4.775 2.720 ( 0.476 -1.375 -0.282) 1.482 3.047 ( -0.055 -0.094 0.020) 0.110 3.188 ( -0.040 6.837 -0.004) 6.837 3.541 ( 1.969 -4.454 0.383) 4.885 ======================= Grid point 578 (121/126) ======================= q-point: (-0.38 0.15 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.708 ( -5.876 -0.218 0.829) 5.938 1.156 ( -10.929 4.520 0.352) 11.832 1.451 ( -8.168 -3.203 -0.207) 8.776 2.288 ( -1.205 1.914 -0.252) 2.276 2.439 ( 2.077 -3.449 0.024) 4.026 2.709 ( -0.362 -0.773 -0.062) 0.856 3.032 ( -1.994 0.622 -0.082) 2.090 3.301 ( 5.371 5.658 -0.076) 7.802 3.521 ( 0.479 -3.896 0.109) 3.927 ======================= Grid point 586 (122/126) ======================= q-point: ( 0.23 0.23 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.818 ( 3.617 4.806 0.088) 6.016 1.367 ( 4.869 8.560 0.292) 9.852 1.564 ( 3.728 5.088 -0.487) 6.326 2.320 ( 0.803 0.574 -0.219) 1.012 2.370 ( -2.146 -2.062 0.116) 2.979 2.704 ( -0.437 -0.277 0.127) 0.533 3.046 ( -0.268 -0.541 -0.016) 0.604 3.276 ( -0.763 -0.269 0.537) 0.971 3.482 ( -0.212 -2.087 -0.712) 2.215 ======================= Grid point 587 (123/126) ======================= q-point: ( 0.31 0.23 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.915 ( 2.281 3.198 -0.046) 3.928 1.487 ( -1.190 9.140 0.496) 9.231 1.675 ( 4.455 1.333 -0.506) 4.678 2.320 ( -1.396 1.257 -0.731) 2.016 2.337 ( 0.492 -2.191 0.422) 2.284 2.688 ( -0.137 -1.394 0.061) 1.402 3.026 ( -0.234 -1.712 -0.047) 1.728 3.265 ( -2.715 4.070 0.531) 4.921 3.470 ( 2.046 -3.968 -0.559) 4.499 ======================= Grid point 588 (124/126) ======================= q-point: ( 0.38 0.23 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.950 ( -1.022 1.255 0.378) 1.662 1.523 ( -4.492 8.027 0.381) 9.206 1.720 ( 0.466 -1.573 0.199) 1.653 2.312 ( -1.255 0.783 0.042) 1.480 2.334 ( 1.418 -0.619 -0.873) 1.777 2.680 ( 0.977 -1.496 -0.167) 1.795 3.014 ( 0.788 -1.495 0.031) 1.690 3.263 ( -2.742 5.050 0.276) 5.753 3.466 ( 2.175 -4.226 -0.070) 4.753 ======================= Grid point 590 (125/126) ======================= q-point: (-0.54 0.23 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.908 ( -4.332 0.244 0.615) 4.382 1.490 ( -7.400 5.670 0.233) 9.325 1.665 ( -3.667 -2.667 0.376) 4.549 2.308 ( -0.082 1.412 0.319) 1.450 2.353 ( 1.583 -1.562 -1.121) 2.491 2.691 ( 1.429 -0.813 -0.194) 1.656 3.025 ( 1.565 -0.539 0.087) 1.658 3.269 ( -2.003 4.226 0.136) 4.678 3.463 ( 2.062 -3.369 0.069) 3.950 ======================= Grid point 601 (126/126) ======================= q-point: ( 0.31 0.31 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 8.55e-05 0.00e+00 0.00e+00 Number of triplets: 592 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.965 ( 0.912 1.440 0.521) 1.782 1.626 ( 2.066 4.012 0.517) 4.543 1.686 ( 0.274 -0.104 0.250) 0.385 2.316 ( 0.475 0.232 -0.334) 0.625 2.332 ( -0.608 -0.332 -0.769) 1.035 2.663 ( -0.525 -0.856 -0.193) 1.023 2.997 ( -0.487 -0.883 0.029) 1.009 3.334 ( 1.296 2.506 0.280) 2.835 3.405 ( -1.179 -2.343 0.037) 2.624 =================== 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/10647 10.0 67.240 67.240 8.597 0.000 0.000 0.000 3/10647 20.0 39.470 39.470 3.764 0.000 0.000 0.000 3/10647 30.0 24.772 24.772 2.273 0.000 0.000 0.000 3/10647 40.0 17.381 17.381 1.629 0.000 0.000 0.000 3/10647 50.0 13.296 13.296 1.283 0.000 0.000 0.000 3/10647 60.0 10.769 10.769 1.066 0.000 0.000 0.000 3/10647 70.0 9.062 9.062 0.915 0.000 0.000 0.000 3/10647 80.0 7.831 7.831 0.802 0.000 0.000 0.000 3/10647 90.0 6.901 6.901 0.715 0.000 0.000 0.000 3/10647 100.0 6.172 6.172 0.646 0.000 0.000 0.000 3/10647 110.0 5.585 5.585 0.588 0.000 0.000 0.000 3/10647 120.0 5.102 5.102 0.540 0.000 0.000 0.000 3/10647 130.0 4.696 4.696 0.500 0.000 0.000 0.000 3/10647 140.0 4.352 4.352 0.465 0.000 0.000 0.000 3/10647 150.0 4.054 4.054 0.434 0.000 0.000 0.000 3/10647 160.0 3.796 3.796 0.408 0.000 0.000 0.000 3/10647 170.0 3.568 3.568 0.384 0.000 0.000 0.000 3/10647 180.0 3.367 3.367 0.363 0.000 0.000 0.000 3/10647 190.0 3.187 3.187 0.344 0.000 0.000 0.000 3/10647 200.0 3.026 3.026 0.327 0.000 0.000 0.000 3/10647 210.0 2.880 2.880 0.312 0.000 0.000 0.000 3/10647 220.0 2.748 2.748 0.298 0.000 0.000 0.000 3/10647 230.0 2.627 2.627 0.285 0.000 0.000 0.000 3/10647 240.0 2.517 2.517 0.273 0.000 0.000 0.000 3/10647 250.0 2.415 2.415 0.262 0.000 0.000 0.000 3/10647 260.0 2.322 2.322 0.252 0.000 0.000 0.000 3/10647 270.0 2.235 2.235 0.243 0.000 0.000 0.000 3/10647 280.0 2.155 2.155 0.234 0.000 0.000 0.000 3/10647 290.0 2.080 2.080 0.226 0.000 0.000 0.000 3/10647 300.0 2.011 2.011 0.219 0.000 0.000 0.000 3/10647 310.0 1.945 1.945 0.212 0.000 0.000 0.000 3/10647 320.0 1.884 1.884 0.205 0.000 0.000 0.000 3/10647 330.0 1.827 1.827 0.199 0.000 0.000 0.000 3/10647 340.0 1.773 1.773 0.193 0.000 0.000 0.000 3/10647 350.0 1.722 1.722 0.188 0.000 0.000 0.000 3/10647 360.0 1.674 1.674 0.183 0.000 0.000 0.000 3/10647 370.0 1.629 1.629 0.178 0.000 0.000 0.000 3/10647 380.0 1.586 1.586 0.173 0.000 0.000 0.000 3/10647 390.0 1.545 1.545 0.169 0.000 0.000 0.000 3/10647 400.0 1.506 1.506 0.164 0.000 0.000 0.000 3/10647 410.0 1.470 1.470 0.160 0.000 0.000 0.000 3/10647 420.0 1.435 1.435 0.157 0.000 0.000 0.000 3/10647 430.0 1.401 1.401 0.153 0.000 0.000 0.000 3/10647 440.0 1.369 1.369 0.150 0.000 0.000 0.000 3/10647 450.0 1.339 1.339 0.146 0.000 0.000 0.000 3/10647 460.0 1.310 1.310 0.143 0.000 0.000 0.000 3/10647 470.0 1.282 1.282 0.140 0.000 0.000 0.000 3/10647 480.0 1.255 1.255 0.137 0.000 0.000 0.000 3/10647 490.0 1.229 1.229 0.134 0.000 0.000 0.000 3/10647 500.0 1.205 1.205 0.132 0.000 0.000 0.000 3/10647 510.0 1.181 1.181 0.129 0.000 0.000 0.000 3/10647 520.0 1.158 1.158 0.127 0.000 0.000 0.000 3/10647 530.0 1.136 1.136 0.124 0.000 0.000 0.000 3/10647 540.0 1.115 1.115 0.122 0.000 0.000 0.000 3/10647 550.0 1.095 1.095 0.120 0.000 0.000 0.000 3/10647 560.0 1.075 1.075 0.118 0.000 0.000 0.000 3/10647 570.0 1.057 1.057 0.116 0.000 0.000 0.000 3/10647 580.0 1.038 1.038 0.114 0.000 0.000 0.000 3/10647 590.0 1.021 1.021 0.112 0.000 0.000 0.000 3/10647 600.0 1.004 1.004 0.110 0.000 0.000 0.000 3/10647 610.0 0.987 0.987 0.108 0.000 0.000 0.000 3/10647 620.0 0.971 0.971 0.106 0.000 0.000 0.000 3/10647 630.0 0.956 0.956 0.105 0.000 0.000 0.000 3/10647 640.0 0.941 0.941 0.103 0.000 0.000 0.000 3/10647 650.0 0.926 0.926 0.101 0.000 0.000 0.000 3/10647 660.0 0.912 0.912 0.100 0.000 0.000 0.000 3/10647 670.0 0.899 0.899 0.098 0.000 0.000 0.000 3/10647 680.0 0.886 0.886 0.097 0.000 0.000 0.000 3/10647 690.0 0.873 0.873 0.095 0.000 0.000 0.000 3/10647 700.0 0.860 0.860 0.094 0.000 0.000 0.000 3/10647 710.0 0.848 0.848 0.093 0.000 0.000 0.000 3/10647 720.0 0.836 0.836 0.091 0.000 0.000 0.000 3/10647 730.0 0.825 0.825 0.090 0.000 0.000 0.000 3/10647 740.0 0.814 0.814 0.089 0.000 0.000 0.000 3/10647 750.0 0.803 0.803 0.088 0.000 0.000 0.000 3/10647 760.0 0.792 0.792 0.087 0.000 0.000 0.000 3/10647 770.0 0.782 0.782 0.086 0.000 0.000 0.000 3/10647 780.0 0.772 0.772 0.084 0.000 0.000 0.000 3/10647 790.0 0.762 0.762 0.083 0.000 0.000 0.000 3/10647 800.0 0.753 0.753 0.082 0.000 0.000 0.000 3/10647 810.0 0.743 0.743 0.081 0.000 0.000 0.000 3/10647 820.0 0.734 0.734 0.080 0.000 0.000 0.000 3/10647 830.0 0.725 0.725 0.079 0.000 0.000 0.000 3/10647 840.0 0.717 0.717 0.078 0.000 0.000 0.000 3/10647 850.0 0.708 0.708 0.078 0.000 0.000 0.000 3/10647 860.0 0.700 0.700 0.077 0.000 0.000 0.000 3/10647 870.0 0.692 0.692 0.076 0.000 0.000 0.000 3/10647 880.0 0.684 0.684 0.075 0.000 0.000 0.000 3/10647 890.0 0.677 0.677 0.074 0.000 0.000 0.000 3/10647 900.0 0.669 0.669 0.073 0.000 0.000 0.000 3/10647 910.0 0.662 0.662 0.072 0.000 0.000 0.000 3/10647 920.0 0.654 0.654 0.072 0.000 0.000 0.000 3/10647 930.0 0.647 0.647 0.071 0.000 0.000 0.000 3/10647 940.0 0.641 0.641 0.070 0.000 0.000 0.000 3/10647 950.0 0.634 0.634 0.069 0.000 0.000 0.000 3/10647 960.0 0.627 0.627 0.069 0.000 0.000 0.000 3/10647 970.0 0.621 0.621 0.068 0.000 0.000 0.000 3/10647 980.0 0.614 0.614 0.067 0.000 0.000 0.000 3/10647 990.0 0.608 0.608 0.067 0.000 0.000 0.000 3/10647 1000.0 0.602 0.602 0.066 0.000 0.000 0.000 3/10647 Thermal conductivity related properties were written into "kappa-m13137.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 18:55:51]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|