----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-09 04:29:48]------------------------- 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 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) Number of symmetry operations in supercell: 384 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.142169519522996 0.000000000000000 0.000000000000000 b -2.071084759761498 3.587224030688497 0.000000000000000 c 0.000000000000000 0.000000000000000 6.798646230000000 Atomic positions (fractional): *1 In 0.33333333333333 0.66666666666667 0.00024968091662 114.818 2 In 0.66666666666667 0.33333333333333 0.50024968091662 114.818 *3 P 0.33333333333333 0.66666666666667 0.37475131908338 30.974 4 P 0.66666666666667 0.33333333333333 0.87475131908338 30.974 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.142169519522996 0.000000000000000 0.000000000000000 b -2.071084759761498 3.587224030688497 0.000000000000000 c 0.000000000000000 0.000000000000000 6.798646230000000 Atomic positions (fractional): *1 In 0.33333333333333 0.66666666666667 0.00024968091662 114.818 > 1 2 In 0.66666666666667 0.33333333333333 0.50024968091662 114.818 > 2 *3 P 0.33333333333333 0.66666666666667 0.37475131908338 30.974 > 3 4 P 0.66666666666667 0.33333333333333 0.87475131908338 30.974 > 4 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.568678078091985 0.000000000000000 0.000000000000000 b -8.284339039045992 14.348896122753988 0.000000000000000 c 0.000000000000000 0.000000000000000 13.597292460000000 Atomic positions (fractional): *1 In 0.08333333333333 0.16666666666667 0.00012484045831 114.818 > 1 2 In 0.33333333333333 0.16666666666667 0.00012484045831 114.818 > 1 3 In 0.58333333333333 0.16666666666667 0.00012484045831 114.818 > 1 4 In 0.83333333333333 0.16666666666667 0.00012484045831 114.818 > 1 5 In 0.08333333333333 0.41666666666667 0.00012484045831 114.818 > 1 6 In 0.33333333333333 0.41666666666667 0.00012484045831 114.818 > 1 7 In 0.58333333333333 0.41666666666667 0.00012484045831 114.818 > 1 8 In 0.83333333333333 0.41666666666667 0.00012484045831 114.818 > 1 9 In 0.08333333333333 0.66666666666667 0.00012484045831 114.818 > 1 10 In 0.33333333333333 0.66666666666667 0.00012484045831 114.818 > 1 11 In 0.58333333333333 0.66666666666667 0.00012484045831 114.818 > 1 12 In 0.83333333333333 0.66666666666667 0.00012484045831 114.818 > 1 13 In 0.08333333333333 0.91666666666667 0.00012484045831 114.818 > 1 14 In 0.33333333333333 0.91666666666667 0.00012484045831 114.818 > 1 15 In 0.58333333333333 0.91666666666667 0.00012484045831 114.818 > 1 16 In 0.83333333333333 0.91666666666667 0.00012484045831 114.818 > 1 17 In 0.08333333333333 0.16666666666667 0.50012484045831 114.818 > 1 18 In 0.33333333333333 0.16666666666667 0.50012484045831 114.818 > 1 19 In 0.58333333333333 0.16666666666667 0.50012484045831 114.818 > 1 20 In 0.83333333333333 0.16666666666667 0.50012484045831 114.818 > 1 21 In 0.08333333333333 0.41666666666667 0.50012484045831 114.818 > 1 22 In 0.33333333333333 0.41666666666667 0.50012484045831 114.818 > 1 23 In 0.58333333333333 0.41666666666667 0.50012484045831 114.818 > 1 24 In 0.83333333333333 0.41666666666667 0.50012484045831 114.818 > 1 25 In 0.08333333333333 0.66666666666667 0.50012484045831 114.818 > 1 26 In 0.33333333333333 0.66666666666667 0.50012484045831 114.818 > 1 27 In 0.58333333333333 0.66666666666667 0.50012484045831 114.818 > 1 28 In 0.83333333333333 0.66666666666667 0.50012484045831 114.818 > 1 29 In 0.08333333333333 0.91666666666667 0.50012484045831 114.818 > 1 30 In 0.33333333333333 0.91666666666667 0.50012484045831 114.818 > 1 31 In 0.58333333333333 0.91666666666667 0.50012484045831 114.818 > 1 32 In 0.83333333333333 0.91666666666667 0.50012484045831 114.818 > 1 33 In 0.16666666666667 0.08333333333333 0.25012484045831 114.818 > 2 34 In 0.41666666666667 0.08333333333333 0.25012484045831 114.818 > 2 35 In 0.66666666666667 0.08333333333333 0.25012484045831 114.818 > 2 36 In 0.91666666666667 0.08333333333333 0.25012484045831 114.818 > 2 37 In 0.16666666666667 0.33333333333333 0.25012484045831 114.818 > 2 38 In 0.41666666666667 0.33333333333333 0.25012484045831 114.818 > 2 39 In 0.66666666666667 0.33333333333333 0.25012484045831 114.818 > 2 40 In 0.91666666666667 0.33333333333333 0.25012484045831 114.818 > 2 41 In 0.16666666666667 0.58333333333333 0.25012484045831 114.818 > 2 42 In 0.41666666666667 0.58333333333333 0.25012484045831 114.818 > 2 43 In 0.66666666666667 0.58333333333333 0.25012484045831 114.818 > 2 44 In 0.91666666666667 0.58333333333333 0.25012484045831 114.818 > 2 45 In 0.16666666666667 0.83333333333333 0.25012484045831 114.818 > 2 46 In 0.41666666666667 0.83333333333333 0.25012484045831 114.818 > 2 47 In 0.66666666666667 0.83333333333333 0.25012484045831 114.818 > 2 48 In 0.91666666666667 0.83333333333333 0.25012484045831 114.818 > 2 49 In 0.16666666666667 0.08333333333333 0.75012484045831 114.818 > 2 50 In 0.41666666666667 0.08333333333333 0.75012484045831 114.818 > 2 51 In 0.66666666666667 0.08333333333333 0.75012484045831 114.818 > 2 52 In 0.91666666666667 0.08333333333333 0.75012484045831 114.818 > 2 53 In 0.16666666666667 0.33333333333333 0.75012484045831 114.818 > 2 54 In 0.41666666666667 0.33333333333333 0.75012484045831 114.818 > 2 55 In 0.66666666666667 0.33333333333333 0.75012484045831 114.818 > 2 56 In 0.91666666666667 0.33333333333333 0.75012484045831 114.818 > 2 57 In 0.16666666666667 0.58333333333333 0.75012484045831 114.818 > 2 58 In 0.41666666666667 0.58333333333333 0.75012484045831 114.818 > 2 59 In 0.66666666666667 0.58333333333333 0.75012484045831 114.818 > 2 60 In 0.91666666666667 0.58333333333333 0.75012484045831 114.818 > 2 61 In 0.16666666666667 0.83333333333333 0.75012484045831 114.818 > 2 62 In 0.41666666666667 0.83333333333333 0.75012484045831 114.818 > 2 63 In 0.66666666666667 0.83333333333333 0.75012484045831 114.818 > 2 64 In 0.91666666666667 0.83333333333333 0.75012484045831 114.818 > 2 *65 P 0.08333333333333 0.16666666666667 0.18737565954169 30.974 > 3 66 P 0.33333333333333 0.16666666666667 0.18737565954169 30.974 > 3 67 P 0.58333333333333 0.16666666666667 0.18737565954169 30.974 > 3 68 P 0.83333333333333 0.16666666666667 0.18737565954169 30.974 > 3 69 P 0.08333333333333 0.41666666666667 0.18737565954169 30.974 > 3 70 P 0.33333333333333 0.41666666666667 0.18737565954169 30.974 > 3 71 P 0.58333333333333 0.41666666666667 0.18737565954169 30.974 > 3 72 P 0.83333333333333 0.41666666666667 0.18737565954169 30.974 > 3 73 P 0.08333333333333 0.66666666666667 0.18737565954169 30.974 > 3 74 P 0.33333333333333 0.66666666666667 0.18737565954169 30.974 > 3 75 P 0.58333333333333 0.66666666666667 0.18737565954169 30.974 > 3 76 P 0.83333333333333 0.66666666666667 0.18737565954169 30.974 > 3 77 P 0.08333333333333 0.91666666666667 0.18737565954169 30.974 > 3 78 P 0.33333333333333 0.91666666666667 0.18737565954169 30.974 > 3 79 P 0.58333333333333 0.91666666666667 0.18737565954169 30.974 > 3 80 P 0.83333333333333 0.91666666666667 0.18737565954169 30.974 > 3 81 P 0.08333333333333 0.16666666666667 0.68737565954169 30.974 > 3 82 P 0.33333333333333 0.16666666666667 0.68737565954169 30.974 > 3 83 P 0.58333333333333 0.16666666666667 0.68737565954169 30.974 > 3 84 P 0.83333333333333 0.16666666666667 0.68737565954169 30.974 > 3 85 P 0.08333333333333 0.41666666666667 0.68737565954169 30.974 > 3 86 P 0.33333333333333 0.41666666666667 0.68737565954169 30.974 > 3 87 P 0.58333333333333 0.41666666666667 0.68737565954169 30.974 > 3 88 P 0.83333333333333 0.41666666666667 0.68737565954169 30.974 > 3 89 P 0.08333333333333 0.66666666666667 0.68737565954169 30.974 > 3 90 P 0.33333333333333 0.66666666666667 0.68737565954169 30.974 > 3 91 P 0.58333333333333 0.66666666666667 0.68737565954169 30.974 > 3 92 P 0.83333333333333 0.66666666666667 0.68737565954169 30.974 > 3 93 P 0.08333333333333 0.91666666666667 0.68737565954169 30.974 > 3 94 P 0.33333333333333 0.91666666666667 0.68737565954169 30.974 > 3 95 P 0.58333333333333 0.91666666666667 0.68737565954169 30.974 > 3 96 P 0.83333333333333 0.91666666666667 0.68737565954169 30.974 > 3 97 P 0.16666666666667 0.08333333333333 0.43737565954169 30.974 > 4 98 P 0.41666666666667 0.08333333333333 0.43737565954169 30.974 > 4 99 P 0.66666666666667 0.08333333333333 0.43737565954169 30.974 > 4 100 P 0.91666666666667 0.08333333333333 0.43737565954169 30.974 > 4 101 P 0.16666666666667 0.33333333333333 0.43737565954169 30.974 > 4 102 P 0.41666666666667 0.33333333333333 0.43737565954169 30.974 > 4 103 P 0.66666666666667 0.33333333333333 0.43737565954169 30.974 > 4 104 P 0.91666666666667 0.33333333333333 0.43737565954169 30.974 > 4 105 P 0.16666666666667 0.58333333333333 0.43737565954169 30.974 > 4 106 P 0.41666666666667 0.58333333333333 0.43737565954169 30.974 > 4 107 P 0.66666666666667 0.58333333333333 0.43737565954169 30.974 > 4 108 P 0.91666666666667 0.58333333333333 0.43737565954169 30.974 > 4 109 P 0.16666666666667 0.83333333333333 0.43737565954169 30.974 > 4 110 P 0.41666666666667 0.83333333333333 0.43737565954169 30.974 > 4 111 P 0.66666666666667 0.83333333333333 0.43737565954169 30.974 > 4 112 P 0.91666666666667 0.83333333333333 0.43737565954169 30.974 > 4 113 P 0.16666666666667 0.08333333333333 0.93737565954169 30.974 > 4 114 P 0.41666666666667 0.08333333333333 0.93737565954169 30.974 > 4 115 P 0.66666666666667 0.08333333333333 0.93737565954169 30.974 > 4 116 P 0.91666666666667 0.08333333333333 0.93737565954169 30.974 > 4 117 P 0.16666666666667 0.33333333333333 0.93737565954169 30.974 > 4 118 P 0.41666666666667 0.33333333333333 0.93737565954169 30.974 > 4 119 P 0.66666666666667 0.33333333333333 0.93737565954169 30.974 > 4 120 P 0.91666666666667 0.33333333333333 0.93737565954169 30.974 > 4 121 P 0.16666666666667 0.58333333333333 0.93737565954169 30.974 > 4 122 P 0.41666666666667 0.58333333333333 0.93737565954169 30.974 > 4 123 P 0.66666666666667 0.58333333333333 0.93737565954169 30.974 > 4 124 P 0.91666666666667 0.58333333333333 0.93737565954169 30.974 > 4 125 P 0.16666666666667 0.83333333333333 0.93737565954169 30.974 > 4 126 P 0.41666666666667 0.83333333333333 0.93737565954169 30.974 > 4 127 P 0.66666666666667 0.83333333333333 0.93737565954169 30.974 > 4 128 P 0.91666666666667 0.83333333333333 0.93737565954169 30.974 > 4 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 21.4794297 -0.0000000 0.0000000 0.0000000 21.4794296 0.0000000 0.0000000 0.0000000 18.8542021 -------------------------- Born effective charges -------------------------- 1 In 2.5737552 -0.0000000 0.0000000 0.0000000 2.5737552 0.0000000 0.0000000 0.0000000 2.8011319 2 In 2.5737552 -0.0000000 0.0000000 0.0000000 2.5737552 0.0000000 0.0000000 0.0000000 2.8011319 3 P -2.5737552 0.0000000 0.0000000 0.0000000 -2.5737552 0.0000000 0.0000000 0.0000000 -2.8011319 4 P -2.5737552 0.0000000 0.0000000 0.0000000 -2.5737552 0.0000000 0.0000000 0.0000000 -2.8011319 ---------------------------------------------------------------------------- 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: 384/384 Permutation basis: 2856/2856 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 224 Number of blocks in projector: 156 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: 98 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 46 Use standard eigh solver. Tree of FC basis block matrices: - (224, 220), data: False |-- (46, 44), data: True |-- (80, 80), data: True |-- (98, 96), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.009 Solver_block: 80 / 80 - Time: 0.167 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.177 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 384/384 Permutation basis: 2856/2856 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 224 Number of blocks in projector: 156 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: 98 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 46 Use standard eigh solver. Tree of FC basis block matrices: - (224, 220), data: False |-- (46, 44), data: True |-- (80, 80), data: True |-- (98, 96), data: True Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-09 04:29:53]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:29:53]------------------------- 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 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.142169519522996 0.000000000000000 0.000000000000000 b -2.071084759761498 3.587224030688497 0.000000000000000 c 0.000000000000000 0.000000000000000 6.798646230000000 Atomic positions (fractional): 1 In 0.33333333333333 0.66666666666667 0.00024968091662 114.818 2 In 0.66666666666667 0.33333333333333 0.50024968091662 114.818 3 P 0.33333333333333 0.66666666666667 0.37475131908338 30.974 4 P 0.66666666666667 0.33333333333333 0.87475131908338 30.974 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.568678078091985 0.000000000000000 0.000000000000000 b -8.284339039045992 14.348896122753988 0.000000000000000 c 0.000000000000000 0.000000000000000 13.597292460000000 Atomic positions (fractional): 1 In 0.08333333333333 0.16666666666667 0.00012484045831 114.818 > 1 2 In 0.33333333333333 0.16666666666667 0.00012484045831 114.818 > 1 3 In 0.58333333333333 0.16666666666667 0.00012484045831 114.818 > 1 4 In 0.83333333333333 0.16666666666667 0.00012484045831 114.818 > 1 5 In 0.08333333333333 0.41666666666667 0.00012484045831 114.818 > 1 6 In 0.33333333333333 0.41666666666667 0.00012484045831 114.818 > 1 7 In 0.58333333333333 0.41666666666667 0.00012484045831 114.818 > 1 8 In 0.83333333333333 0.41666666666667 0.00012484045831 114.818 > 1 9 In 0.08333333333333 0.66666666666667 0.00012484045831 114.818 > 1 10 In 0.33333333333333 0.66666666666667 0.00012484045831 114.818 > 1 11 In 0.58333333333333 0.66666666666667 0.00012484045831 114.818 > 1 12 In 0.83333333333333 0.66666666666667 0.00012484045831 114.818 > 1 13 In 0.08333333333333 0.91666666666667 0.00012484045831 114.818 > 1 14 In 0.33333333333333 0.91666666666667 0.00012484045831 114.818 > 1 15 In 0.58333333333333 0.91666666666667 0.00012484045831 114.818 > 1 16 In 0.83333333333333 0.91666666666667 0.00012484045831 114.818 > 1 17 In 0.08333333333333 0.16666666666667 0.50012484045831 114.818 > 1 18 In 0.33333333333333 0.16666666666667 0.50012484045831 114.818 > 1 19 In 0.58333333333333 0.16666666666667 0.50012484045831 114.818 > 1 20 In 0.83333333333333 0.16666666666667 0.50012484045831 114.818 > 1 21 In 0.08333333333333 0.41666666666667 0.50012484045831 114.818 > 1 22 In 0.33333333333333 0.41666666666667 0.50012484045831 114.818 > 1 23 In 0.58333333333333 0.41666666666667 0.50012484045831 114.818 > 1 24 In 0.83333333333333 0.41666666666667 0.50012484045831 114.818 > 1 25 In 0.08333333333333 0.66666666666667 0.50012484045831 114.818 > 1 26 In 0.33333333333333 0.66666666666667 0.50012484045831 114.818 > 1 27 In 0.58333333333333 0.66666666666667 0.50012484045831 114.818 > 1 28 In 0.83333333333333 0.66666666666667 0.50012484045831 114.818 > 1 29 In 0.08333333333333 0.91666666666667 0.50012484045831 114.818 > 1 30 In 0.33333333333333 0.91666666666667 0.50012484045831 114.818 > 1 31 In 0.58333333333333 0.91666666666667 0.50012484045831 114.818 > 1 32 In 0.83333333333333 0.91666666666667 0.50012484045831 114.818 > 1 33 In 0.16666666666667 0.08333333333333 0.25012484045831 114.818 > 33 34 In 0.41666666666667 0.08333333333333 0.25012484045831 114.818 > 33 35 In 0.66666666666667 0.08333333333333 0.25012484045831 114.818 > 33 36 In 0.91666666666667 0.08333333333333 0.25012484045831 114.818 > 33 37 In 0.16666666666667 0.33333333333333 0.25012484045831 114.818 > 33 38 In 0.41666666666667 0.33333333333333 0.25012484045831 114.818 > 33 39 In 0.66666666666667 0.33333333333333 0.25012484045831 114.818 > 33 40 In 0.91666666666667 0.33333333333333 0.25012484045831 114.818 > 33 41 In 0.16666666666667 0.58333333333333 0.25012484045831 114.818 > 33 42 In 0.41666666666667 0.58333333333333 0.25012484045831 114.818 > 33 43 In 0.66666666666667 0.58333333333333 0.25012484045831 114.818 > 33 44 In 0.91666666666667 0.58333333333333 0.25012484045831 114.818 > 33 45 In 0.16666666666667 0.83333333333333 0.25012484045831 114.818 > 33 46 In 0.41666666666667 0.83333333333333 0.25012484045831 114.818 > 33 47 In 0.66666666666667 0.83333333333333 0.25012484045831 114.818 > 33 48 In 0.91666666666667 0.83333333333333 0.25012484045831 114.818 > 33 49 In 0.16666666666667 0.08333333333333 0.75012484045831 114.818 > 33 50 In 0.41666666666667 0.08333333333333 0.75012484045831 114.818 > 33 51 In 0.66666666666667 0.08333333333333 0.75012484045831 114.818 > 33 52 In 0.91666666666667 0.08333333333333 0.75012484045831 114.818 > 33 53 In 0.16666666666667 0.33333333333333 0.75012484045831 114.818 > 33 54 In 0.41666666666667 0.33333333333333 0.75012484045831 114.818 > 33 55 In 0.66666666666667 0.33333333333333 0.75012484045831 114.818 > 33 56 In 0.91666666666667 0.33333333333333 0.75012484045831 114.818 > 33 57 In 0.16666666666667 0.58333333333333 0.75012484045831 114.818 > 33 58 In 0.41666666666667 0.58333333333333 0.75012484045831 114.818 > 33 59 In 0.66666666666667 0.58333333333333 0.75012484045831 114.818 > 33 60 In 0.91666666666667 0.58333333333333 0.75012484045831 114.818 > 33 61 In 0.16666666666667 0.83333333333333 0.75012484045831 114.818 > 33 62 In 0.41666666666667 0.83333333333333 0.75012484045831 114.818 > 33 63 In 0.66666666666667 0.83333333333333 0.75012484045831 114.818 > 33 64 In 0.91666666666667 0.83333333333333 0.75012484045831 114.818 > 33 65 P 0.08333333333333 0.16666666666667 0.18737565954169 30.974 > 65 66 P 0.33333333333333 0.16666666666667 0.18737565954169 30.974 > 65 67 P 0.58333333333333 0.16666666666667 0.18737565954169 30.974 > 65 68 P 0.83333333333333 0.16666666666667 0.18737565954169 30.974 > 65 69 P 0.08333333333333 0.41666666666667 0.18737565954169 30.974 > 65 70 P 0.33333333333333 0.41666666666667 0.18737565954169 30.974 > 65 71 P 0.58333333333333 0.41666666666667 0.18737565954169 30.974 > 65 72 P 0.83333333333333 0.41666666666667 0.18737565954169 30.974 > 65 73 P 0.08333333333333 0.66666666666667 0.18737565954169 30.974 > 65 74 P 0.33333333333333 0.66666666666667 0.18737565954169 30.974 > 65 75 P 0.58333333333333 0.66666666666667 0.18737565954169 30.974 > 65 76 P 0.83333333333333 0.66666666666667 0.18737565954169 30.974 > 65 77 P 0.08333333333333 0.91666666666667 0.18737565954169 30.974 > 65 78 P 0.33333333333333 0.91666666666667 0.18737565954169 30.974 > 65 79 P 0.58333333333333 0.91666666666667 0.18737565954169 30.974 > 65 80 P 0.83333333333333 0.91666666666667 0.18737565954169 30.974 > 65 81 P 0.08333333333333 0.16666666666667 0.68737565954169 30.974 > 65 82 P 0.33333333333333 0.16666666666667 0.68737565954169 30.974 > 65 83 P 0.58333333333333 0.16666666666667 0.68737565954169 30.974 > 65 84 P 0.83333333333333 0.16666666666667 0.68737565954169 30.974 > 65 85 P 0.08333333333333 0.41666666666667 0.68737565954169 30.974 > 65 86 P 0.33333333333333 0.41666666666667 0.68737565954169 30.974 > 65 87 P 0.58333333333333 0.41666666666667 0.68737565954169 30.974 > 65 88 P 0.83333333333333 0.41666666666667 0.68737565954169 30.974 > 65 89 P 0.08333333333333 0.66666666666667 0.68737565954169 30.974 > 65 90 P 0.33333333333333 0.66666666666667 0.68737565954169 30.974 > 65 91 P 0.58333333333333 0.66666666666667 0.68737565954169 30.974 > 65 92 P 0.83333333333333 0.66666666666667 0.68737565954169 30.974 > 65 93 P 0.08333333333333 0.91666666666667 0.68737565954169 30.974 > 65 94 P 0.33333333333333 0.91666666666667 0.68737565954169 30.974 > 65 95 P 0.58333333333333 0.91666666666667 0.68737565954169 30.974 > 65 96 P 0.83333333333333 0.91666666666667 0.68737565954169 30.974 > 65 97 P 0.16666666666667 0.08333333333333 0.43737565954169 30.974 > 97 98 P 0.41666666666667 0.08333333333333 0.43737565954169 30.974 > 97 99 P 0.66666666666667 0.08333333333333 0.43737565954169 30.974 > 97 100 P 0.91666666666667 0.08333333333333 0.43737565954169 30.974 > 97 101 P 0.16666666666667 0.33333333333333 0.43737565954169 30.974 > 97 102 P 0.41666666666667 0.33333333333333 0.43737565954169 30.974 > 97 103 P 0.66666666666667 0.33333333333333 0.43737565954169 30.974 > 97 104 P 0.91666666666667 0.33333333333333 0.43737565954169 30.974 > 97 105 P 0.16666666666667 0.58333333333333 0.43737565954169 30.974 > 97 106 P 0.41666666666667 0.58333333333333 0.43737565954169 30.974 > 97 107 P 0.66666666666667 0.58333333333333 0.43737565954169 30.974 > 97 108 P 0.91666666666667 0.58333333333333 0.43737565954169 30.974 > 97 109 P 0.16666666666667 0.83333333333333 0.43737565954169 30.974 > 97 110 P 0.41666666666667 0.83333333333333 0.43737565954169 30.974 > 97 111 P 0.66666666666667 0.83333333333333 0.43737565954169 30.974 > 97 112 P 0.91666666666667 0.83333333333333 0.43737565954169 30.974 > 97 113 P 0.16666666666667 0.08333333333333 0.93737565954169 30.974 > 97 114 P 0.41666666666667 0.08333333333333 0.93737565954169 30.974 > 97 115 P 0.66666666666667 0.08333333333333 0.93737565954169 30.974 > 97 116 P 0.91666666666667 0.08333333333333 0.93737565954169 30.974 > 97 117 P 0.16666666666667 0.33333333333333 0.93737565954169 30.974 > 97 118 P 0.41666666666667 0.33333333333333 0.93737565954169 30.974 > 97 119 P 0.66666666666667 0.33333333333333 0.93737565954169 30.974 > 97 120 P 0.91666666666667 0.33333333333333 0.93737565954169 30.974 > 97 121 P 0.16666666666667 0.58333333333333 0.93737565954169 30.974 > 97 122 P 0.41666666666667 0.58333333333333 0.93737565954169 30.974 > 97 123 P 0.66666666666667 0.58333333333333 0.93737565954169 30.974 > 97 124 P 0.91666666666667 0.58333333333333 0.93737565954169 30.974 > 97 125 P 0.16666666666667 0.83333333333333 0.93737565954169 30.974 > 97 126 P 0.41666666666667 0.83333333333333 0.93737565954169 30.974 > 97 127 P 0.66666666666667 0.83333333333333 0.93737565954169 30.974 > 97 128 P 0.91666666666667 0.83333333333333 0.93737565954169 30.974 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 21.4794297 -0.0000000 0.0000000 0.0000000 21.4794296 0.0000000 0.0000000 0.0000000 18.8542021 -------------------------- Born effective charges -------------------------- 1 In 2.5737552 -0.0000000 0.0000000 0.0000000 2.5737552 0.0000000 0.0000000 0.0000000 2.8011319 2 In 2.5737552 -0.0000000 0.0000000 0.0000000 2.5737552 0.0000000 0.0000000 0.0000000 2.8011319 3 P -2.5737552 0.0000000 0.0000000 0.0000000 -2.5737552 0.0000000 0.0000000 0.0000000 -2.8011319 4 P -2.5737552 0.0000000 0.0000000 0.0000000 -2.5737552 0.0000000 0.0000000 0.0000000 -2.8011319 ---------------------------------------------------------------------------- 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[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000118 (zzz) 0.00000118 (zzz) 0.00000118 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 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-09 04:29:58]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:29:59]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.142169519522996 0.000000000000000 0.000000000000000 b -2.071084759761498 3.587224030688497 0.000000000000000 c 0.000000000000000 0.000000000000000 6.798646230000000 Atomic positions (fractional): 1 In 0.33333333333333 0.66666666666667 0.00024968091662 114.818 2 In 0.66666666666667 0.33333333333333 0.50024968091662 114.818 3 P 0.33333333333333 0.66666666666667 0.37475131908338 30.974 4 P 0.66666666666667 0.33333333333333 0.87475131908338 30.974 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.568678078091985 0.000000000000000 0.000000000000000 b -8.284339039045992 14.348896122753988 0.000000000000000 c 0.000000000000000 0.000000000000000 13.597292460000000 Atomic positions (fractional): 1 In 0.08333333333333 0.16666666666667 0.00012484045831 114.818 > 1 2 In 0.33333333333333 0.16666666666667 0.00012484045831 114.818 > 1 3 In 0.58333333333333 0.16666666666667 0.00012484045831 114.818 > 1 4 In 0.83333333333333 0.16666666666667 0.00012484045831 114.818 > 1 5 In 0.08333333333333 0.41666666666667 0.00012484045831 114.818 > 1 6 In 0.33333333333333 0.41666666666667 0.00012484045831 114.818 > 1 7 In 0.58333333333333 0.41666666666667 0.00012484045831 114.818 > 1 8 In 0.83333333333333 0.41666666666667 0.00012484045831 114.818 > 1 9 In 0.08333333333333 0.66666666666667 0.00012484045831 114.818 > 1 10 In 0.33333333333333 0.66666666666667 0.00012484045831 114.818 > 1 11 In 0.58333333333333 0.66666666666667 0.00012484045831 114.818 > 1 12 In 0.83333333333333 0.66666666666667 0.00012484045831 114.818 > 1 13 In 0.08333333333333 0.91666666666667 0.00012484045831 114.818 > 1 14 In 0.33333333333333 0.91666666666667 0.00012484045831 114.818 > 1 15 In 0.58333333333333 0.91666666666667 0.00012484045831 114.818 > 1 16 In 0.83333333333333 0.91666666666667 0.00012484045831 114.818 > 1 17 In 0.08333333333333 0.16666666666667 0.50012484045831 114.818 > 1 18 In 0.33333333333333 0.16666666666667 0.50012484045831 114.818 > 1 19 In 0.58333333333333 0.16666666666667 0.50012484045831 114.818 > 1 20 In 0.83333333333333 0.16666666666667 0.50012484045831 114.818 > 1 21 In 0.08333333333333 0.41666666666667 0.50012484045831 114.818 > 1 22 In 0.33333333333333 0.41666666666667 0.50012484045831 114.818 > 1 23 In 0.58333333333333 0.41666666666667 0.50012484045831 114.818 > 1 24 In 0.83333333333333 0.41666666666667 0.50012484045831 114.818 > 1 25 In 0.08333333333333 0.66666666666667 0.50012484045831 114.818 > 1 26 In 0.33333333333333 0.66666666666667 0.50012484045831 114.818 > 1 27 In 0.58333333333333 0.66666666666667 0.50012484045831 114.818 > 1 28 In 0.83333333333333 0.66666666666667 0.50012484045831 114.818 > 1 29 In 0.08333333333333 0.91666666666667 0.50012484045831 114.818 > 1 30 In 0.33333333333333 0.91666666666667 0.50012484045831 114.818 > 1 31 In 0.58333333333333 0.91666666666667 0.50012484045831 114.818 > 1 32 In 0.83333333333333 0.91666666666667 0.50012484045831 114.818 > 1 33 In 0.16666666666667 0.08333333333333 0.25012484045831 114.818 > 33 34 In 0.41666666666667 0.08333333333333 0.25012484045831 114.818 > 33 35 In 0.66666666666667 0.08333333333333 0.25012484045831 114.818 > 33 36 In 0.91666666666667 0.08333333333333 0.25012484045831 114.818 > 33 37 In 0.16666666666667 0.33333333333333 0.25012484045831 114.818 > 33 38 In 0.41666666666667 0.33333333333333 0.25012484045831 114.818 > 33 39 In 0.66666666666667 0.33333333333333 0.25012484045831 114.818 > 33 40 In 0.91666666666667 0.33333333333333 0.25012484045831 114.818 > 33 41 In 0.16666666666667 0.58333333333333 0.25012484045831 114.818 > 33 42 In 0.41666666666667 0.58333333333333 0.25012484045831 114.818 > 33 43 In 0.66666666666667 0.58333333333333 0.25012484045831 114.818 > 33 44 In 0.91666666666667 0.58333333333333 0.25012484045831 114.818 > 33 45 In 0.16666666666667 0.83333333333333 0.25012484045831 114.818 > 33 46 In 0.41666666666667 0.83333333333333 0.25012484045831 114.818 > 33 47 In 0.66666666666667 0.83333333333333 0.25012484045831 114.818 > 33 48 In 0.91666666666667 0.83333333333333 0.25012484045831 114.818 > 33 49 In 0.16666666666667 0.08333333333333 0.75012484045831 114.818 > 33 50 In 0.41666666666667 0.08333333333333 0.75012484045831 114.818 > 33 51 In 0.66666666666667 0.08333333333333 0.75012484045831 114.818 > 33 52 In 0.91666666666667 0.08333333333333 0.75012484045831 114.818 > 33 53 In 0.16666666666667 0.33333333333333 0.75012484045831 114.818 > 33 54 In 0.41666666666667 0.33333333333333 0.75012484045831 114.818 > 33 55 In 0.66666666666667 0.33333333333333 0.75012484045831 114.818 > 33 56 In 0.91666666666667 0.33333333333333 0.75012484045831 114.818 > 33 57 In 0.16666666666667 0.58333333333333 0.75012484045831 114.818 > 33 58 In 0.41666666666667 0.58333333333333 0.75012484045831 114.818 > 33 59 In 0.66666666666667 0.58333333333333 0.75012484045831 114.818 > 33 60 In 0.91666666666667 0.58333333333333 0.75012484045831 114.818 > 33 61 In 0.16666666666667 0.83333333333333 0.75012484045831 114.818 > 33 62 In 0.41666666666667 0.83333333333333 0.75012484045831 114.818 > 33 63 In 0.66666666666667 0.83333333333333 0.75012484045831 114.818 > 33 64 In 0.91666666666667 0.83333333333333 0.75012484045831 114.818 > 33 65 P 0.08333333333333 0.16666666666667 0.18737565954169 30.974 > 65 66 P 0.33333333333333 0.16666666666667 0.18737565954169 30.974 > 65 67 P 0.58333333333333 0.16666666666667 0.18737565954169 30.974 > 65 68 P 0.83333333333333 0.16666666666667 0.18737565954169 30.974 > 65 69 P 0.08333333333333 0.41666666666667 0.18737565954169 30.974 > 65 70 P 0.33333333333333 0.41666666666667 0.18737565954169 30.974 > 65 71 P 0.58333333333333 0.41666666666667 0.18737565954169 30.974 > 65 72 P 0.83333333333333 0.41666666666667 0.18737565954169 30.974 > 65 73 P 0.08333333333333 0.66666666666667 0.18737565954169 30.974 > 65 74 P 0.33333333333333 0.66666666666667 0.18737565954169 30.974 > 65 75 P 0.58333333333333 0.66666666666667 0.18737565954169 30.974 > 65 76 P 0.83333333333333 0.66666666666667 0.18737565954169 30.974 > 65 77 P 0.08333333333333 0.91666666666667 0.18737565954169 30.974 > 65 78 P 0.33333333333333 0.91666666666667 0.18737565954169 30.974 > 65 79 P 0.58333333333333 0.91666666666667 0.18737565954169 30.974 > 65 80 P 0.83333333333333 0.91666666666667 0.18737565954169 30.974 > 65 81 P 0.08333333333333 0.16666666666667 0.68737565954169 30.974 > 65 82 P 0.33333333333333 0.16666666666667 0.68737565954169 30.974 > 65 83 P 0.58333333333333 0.16666666666667 0.68737565954169 30.974 > 65 84 P 0.83333333333333 0.16666666666667 0.68737565954169 30.974 > 65 85 P 0.08333333333333 0.41666666666667 0.68737565954169 30.974 > 65 86 P 0.33333333333333 0.41666666666667 0.68737565954169 30.974 > 65 87 P 0.58333333333333 0.41666666666667 0.68737565954169 30.974 > 65 88 P 0.83333333333333 0.41666666666667 0.68737565954169 30.974 > 65 89 P 0.08333333333333 0.66666666666667 0.68737565954169 30.974 > 65 90 P 0.33333333333333 0.66666666666667 0.68737565954169 30.974 > 65 91 P 0.58333333333333 0.66666666666667 0.68737565954169 30.974 > 65 92 P 0.83333333333333 0.66666666666667 0.68737565954169 30.974 > 65 93 P 0.08333333333333 0.91666666666667 0.68737565954169 30.974 > 65 94 P 0.33333333333333 0.91666666666667 0.68737565954169 30.974 > 65 95 P 0.58333333333333 0.91666666666667 0.68737565954169 30.974 > 65 96 P 0.83333333333333 0.91666666666667 0.68737565954169 30.974 > 65 97 P 0.16666666666667 0.08333333333333 0.43737565954169 30.974 > 97 98 P 0.41666666666667 0.08333333333333 0.43737565954169 30.974 > 97 99 P 0.66666666666667 0.08333333333333 0.43737565954169 30.974 > 97 100 P 0.91666666666667 0.08333333333333 0.43737565954169 30.974 > 97 101 P 0.16666666666667 0.33333333333333 0.43737565954169 30.974 > 97 102 P 0.41666666666667 0.33333333333333 0.43737565954169 30.974 > 97 103 P 0.66666666666667 0.33333333333333 0.43737565954169 30.974 > 97 104 P 0.91666666666667 0.33333333333333 0.43737565954169 30.974 > 97 105 P 0.16666666666667 0.58333333333333 0.43737565954169 30.974 > 97 106 P 0.41666666666667 0.58333333333333 0.43737565954169 30.974 > 97 107 P 0.66666666666667 0.58333333333333 0.43737565954169 30.974 > 97 108 P 0.91666666666667 0.58333333333333 0.43737565954169 30.974 > 97 109 P 0.16666666666667 0.83333333333333 0.43737565954169 30.974 > 97 110 P 0.41666666666667 0.83333333333333 0.43737565954169 30.974 > 97 111 P 0.66666666666667 0.83333333333333 0.43737565954169 30.974 > 97 112 P 0.91666666666667 0.83333333333333 0.43737565954169 30.974 > 97 113 P 0.16666666666667 0.08333333333333 0.93737565954169 30.974 > 97 114 P 0.41666666666667 0.08333333333333 0.93737565954169 30.974 > 97 115 P 0.66666666666667 0.08333333333333 0.93737565954169 30.974 > 97 116 P 0.91666666666667 0.08333333333333 0.93737565954169 30.974 > 97 117 P 0.16666666666667 0.33333333333333 0.93737565954169 30.974 > 97 118 P 0.41666666666667 0.33333333333333 0.93737565954169 30.974 > 97 119 P 0.66666666666667 0.33333333333333 0.93737565954169 30.974 > 97 120 P 0.91666666666667 0.33333333333333 0.93737565954169 30.974 > 97 121 P 0.16666666666667 0.58333333333333 0.93737565954169 30.974 > 97 122 P 0.41666666666667 0.58333333333333 0.93737565954169 30.974 > 97 123 P 0.66666666666667 0.58333333333333 0.93737565954169 30.974 > 97 124 P 0.91666666666667 0.58333333333333 0.93737565954169 30.974 > 97 125 P 0.16666666666667 0.83333333333333 0.93737565954169 30.974 > 97 126 P 0.41666666666667 0.83333333333333 0.93737565954169 30.974 > 97 127 P 0.66666666666667 0.83333333333333 0.93737565954169 30.974 > 97 128 P 0.91666666666667 0.83333333333333 0.93737565954169 30.974 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 21.4794297 -0.0000000 0.0000000 0.0000000 21.4794296 0.0000000 0.0000000 0.0000000 18.8542021 -------------------------- Born effective charges -------------------------- 1 In 2.5737552 -0.0000000 0.0000000 0.0000000 2.5737552 0.0000000 0.0000000 0.0000000 2.8011319 2 In 2.5737552 -0.0000000 0.0000000 0.0000000 2.5737552 0.0000000 0.0000000 0.0000000 2.8011319 3 P -2.5737552 0.0000000 0.0000000 0.0000000 -2.5737552 0.0000000 0.0000000 0.0000000 -2.8011319 4 P -2.5737552 0.0000000 0.0000000 0.0000000 -2.5737552 0.0000000 0.0000000 0.0000000 -2.8011319 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000118 (zzz) 0.00000118 (zzz) 0.00000118 (zzz) Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 14 14 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.89, Number of G-points: 313, Lambda: 0.42 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/96) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 0.000 ( 0.000 0.000 0.000) 0.000 1.439 ( 0.000 0.000 0.000) 0.000 1.439 ( 0.000 0.000 0.000) 0.000 4.861 ( 0.000 -0.000 -0.000) 0.000 8.864 ( -0.000 0.000 -0.000) 0.000 8.934 ( 0.000 0.000 -0.000) 0.000 8.934 ( 0.000 0.000 0.000) 0.000 9.044 ( 0.000 -0.000 0.000) 0.000 9.044 ( 0.000 -0.000 0.000) 0.000 9.777 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/96) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.429 ( 18.282 10.555 0.000) 21.110 0.439 ( 19.508 11.263 0.000) 22.526 0.974 ( 41.137 23.751 0.000) 47.501 1.487 ( 4.089 2.361 0.000) 4.721 1.706 ( 20.272 11.704 0.000) 23.409 4.798 ( -5.141 -2.968 0.000) 5.937 8.871 ( 0.572 0.330 0.000) 0.660 8.945 ( 0.971 0.561 0.000) 1.121 8.987 ( 4.131 2.385 0.000) 4.770 9.061 ( 1.386 0.800 0.000) 1.600 9.704 ( 5.494 3.172 0.000) 6.344 9.796 ( 1.866 1.077 0.000) 2.155 ======================= Grid point 2 (3/96) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.829 ( 16.209 9.358 0.000) 18.717 0.876 ( 18.518 10.691 0.000) 21.382 1.617 ( 6.890 3.978 0.000) 7.955 1.874 ( 36.238 20.922 0.000) 41.844 2.222 ( 22.389 12.926 0.000) 25.853 4.654 ( -6.487 -3.745 0.000) 7.490 8.891 ( 1.100 0.635 0.000) 1.271 8.978 ( 1.833 1.059 0.000) 2.117 9.089 ( 3.914 2.260 0.000) 4.519 9.098 ( 1.638 0.946 0.000) 1.891 9.843 ( 5.201 3.003 0.000) 6.006 9.888 ( 6.736 3.889 0.000) 7.778 ======================= Grid point 3 (4/96) ======================= q-point: ( 0.21 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.159 ( 12.165 7.023 0.000) 14.047 1.287 ( 17.046 9.841 0.000) 19.683 1.782 ( 6.904 3.986 0.000) 7.972 2.509 ( 14.956 8.635 0.000) 17.269 2.774 ( 28.528 16.471 0.000) 32.941 4.548 ( -1.794 -1.036 0.000) 2.071 8.919 ( 1.309 0.756 0.000) 1.512 9.024 ( 1.964 1.134 0.000) 2.268 9.133 ( 1.342 0.775 0.000) 1.550 9.147 ( 1.180 0.681 0.000) 1.363 9.953 ( 4.144 2.392 0.000) 4.785 10.046 ( 5.942 3.431 0.000) 6.862 ======================= Grid point 4 (5/96) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.378 ( 6.855 3.958 0.000) 7.916 1.655 ( 14.892 8.598 0.000) 17.196 1.911 ( 3.914 2.260 0.000) 4.519 2.667 ( 1.280 0.739 0.000) 1.478 3.423 ( 26.258 15.160 0.000) 30.320 4.583 ( 4.283 2.473 0.000) 4.946 8.942 ( 0.475 0.274 0.000) 0.548 9.063 ( 1.352 0.780 0.000) 1.561 9.159 ( 0.344 0.199 0.000) 0.397 9.163 ( 1.347 0.778 0.000) 1.556 10.005 ( -0.075 -0.043 0.000) 0.087 10.141 ( 2.164 1.249 0.000) 2.498 ======================= Grid point 5 (6/96) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.480 ( 2.322 1.340 0.000) 2.681 1.953 ( -0.200 -0.116 0.000) 0.231 1.965 ( 11.942 6.894 0.000) 13.789 2.639 ( -2.642 -1.525 0.000) 3.051 3.958 ( 20.119 11.616 0.000) 23.231 4.702 ( 5.121 2.957 0.000) 5.913 8.938 ( -0.712 -0.411 0.000) 0.822 9.087 ( 0.811 0.468 0.000) 0.937 9.174 ( 0.863 0.498 0.000) 0.997 9.194 ( 1.322 0.763 0.000) 1.526 9.955 ( -3.674 -2.121 0.000) 4.242 10.146 ( -1.515 -0.875 0.000) 1.750 ======================= Grid point 6 (7/96) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.505 ( 0.291 0.168 0.000) 0.336 1.919 ( -2.093 -1.208 0.000) 2.417 2.195 ( 7.782 4.493 0.000) 8.986 2.573 ( -2.791 -1.612 0.000) 3.223 4.335 ( 12.259 7.078 0.000) 14.156 4.787 ( 1.948 1.125 0.000) 2.249 8.923 ( -0.400 -0.231 0.000) 0.462 9.102 ( 0.503 0.291 0.000) 0.581 9.187 ( 0.050 0.029 0.000) 0.057 9.219 ( 0.803 0.463 0.000) 0.927 9.872 ( -2.812 -1.623 0.000) 3.247 10.091 ( -2.689 -1.553 0.000) 3.105 ======================= Grid point 7 (8/96) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.507 ( 0.000 0.000 0.000) 0.000 1.889 ( 0.000 0.000 0.000) 0.000 2.296 ( 0.000 0.000 0.000) 0.000 2.532 ( -0.000 -0.000 0.000) 0.000 4.487 ( 0.000 0.000 0.000) 0.000 4.802 ( 0.000 0.000 0.000) 0.000 8.920 ( 0.000 0.000 0.000) 0.000 9.108 ( 0.000 0.000 0.000) 0.000 9.183 ( 0.000 0.000 0.000) 0.000 9.229 ( 0.000 0.000 0.000) 0.000 9.840 ( 0.000 0.000 0.000) 0.000 10.053 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 16 (9/96) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.761 ( 10.822 18.745 0.000) 21.645 0.784 ( 11.900 20.611 0.000) 23.799 1.548 ( 11.391 19.730 0.000) 22.782 1.690 ( 15.030 26.033 0.000) 30.060 2.047 ( 12.074 20.912 0.000) 24.147 4.694 ( -3.900 -6.755 0.000) 7.800 8.882 ( 0.468 0.810 0.000) 0.936 8.997 ( 1.897 3.286 0.000) 3.795 9.039 ( 2.185 3.785 0.000) 4.370 9.094 ( 0.974 1.688 0.000) 1.949 9.811 ( 3.949 6.840 0.000) 7.898 9.842 ( 2.117 3.667 0.000) 4.234 ======================= Grid point 17 (10/96) ======================= q-point: ( 0.14 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.142 ( 13.187 15.180 0.000) 20.108 1.168 ( 7.629 21.392 0.000) 22.711 1.787 ( 4.597 12.265 0.000) 13.098 2.346 ( 18.685 23.597 0.000) 30.099 2.479 ( 21.073 13.283 0.000) 24.910 4.568 ( -3.154 -4.089 0.000) 5.164 8.905 ( 1.086 0.948 0.000) 1.442 9.055 ( -0.275 5.916 0.000) 5.922 9.106 ( 0.272 -0.986 0.000) 1.023 9.137 ( 3.790 3.245 0.000) 4.989 9.922 ( 2.708 4.579 0.000) 5.320 9.964 ( 6.418 3.288 0.000) 7.211 ======================= Grid point 18 (11/96) ======================= q-point: ( 0.21 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.451 ( 1.555 19.738 0.000) 19.799 1.511 ( 13.128 12.339 0.000) 18.016 1.964 ( 0.700 13.184 0.000) 13.202 2.633 ( 4.733 3.977 0.000) 6.182 3.103 ( 27.512 16.124 0.000) 31.889 4.538 ( 2.146 0.139 0.000) 2.150 8.932 ( 0.787 0.629 0.000) 1.007 9.081 ( 0.630 -3.990 0.000) 4.039 9.106 ( -1.279 6.374 0.000) 6.501 9.226 ( 1.385 5.491 0.000) 5.663 9.991 ( 1.398 1.043 0.000) 1.744 10.085 ( 4.905 0.429 0.000) 4.924 ======================= Grid point 19 (12/96) ======================= q-point: ( 0.29 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.611 ( -4.220 17.271 0.000) 17.779 1.837 ( 11.915 9.545 0.000) 15.266 2.060 ( -4.043 11.716 0.000) 12.394 2.667 ( -1.473 -0.293 0.000) 1.502 3.685 ( 23.870 11.132 0.000) 26.338 4.624 ( 5.805 1.322 0.000) 5.953 8.929 ( -0.555 -1.631 0.000) 1.723 9.063 ( 2.685 -4.936 0.000) 5.619 9.141 ( -2.137 6.444 0.000) 6.789 9.296 ( -0.809 7.750 0.000) 7.792 9.981 ( -1.843 -2.080 0.000) 2.779 10.133 ( 1.504 -1.804 0.000) 2.349 ======================= Grid point 20 (13/96) ======================= q-point: ( 0.36 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.665 ( -7.688 15.357 0.000) 17.174 2.062 ( -7.435 10.218 0.000) 12.637 2.106 ( 8.837 7.210 0.000) 11.405 2.614 ( -3.346 -0.856 0.000) 3.454 4.136 ( 18.000 6.252 0.000) 19.055 4.725 ( 4.863 -0.839 0.000) 4.935 8.899 ( 0.281 -2.967 0.000) 2.981 9.059 ( 2.832 -6.257 0.000) 6.868 9.163 ( -2.653 6.328 0.000) 6.861 9.355 ( -2.287 9.260 0.000) 9.539 9.907 ( -3.344 -2.217 0.000) 4.013 10.109 ( -1.377 -2.427 0.000) 2.790 ======================= Grid point 21 (14/96) ======================= q-point: ( 0.43 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.671 ( -8.514 14.727 0.000) 17.011 2.024 ( -7.940 10.729 0.000) 13.348 2.282 ( 3.880 4.286 0.000) 5.782 2.550 ( -2.521 -0.697 0.000) 2.616 4.409 ( 8.838 0.396 0.000) 8.847 4.763 ( 2.312 -3.476 0.000) 4.175 8.888 ( 1.403 -2.349 0.000) 2.736 9.036 ( 3.378 -8.080 0.000) 8.757 9.176 ( -3.148 6.140 0.000) 6.900 9.395 ( -4.355 9.882 0.000) 10.799 9.847 ( -1.385 -0.835 0.000) 1.618 10.062 ( -1.347 -0.945 0.000) 1.646 ======================= Grid point 31 (15/96) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.451 ( 8.959 15.517 0.000) 17.917 1.577 ( 10.382 17.983 0.000) 20.764 2.046 ( 7.141 12.369 0.000) 14.282 2.617 ( 3.702 6.412 0.000) 7.404 2.902 ( 14.953 25.899 0.000) 29.906 4.519 ( -0.226 -0.392 0.000) 0.453 8.929 ( 0.808 1.400 0.000) 1.616 9.055 ( -2.103 -3.642 0.000) 4.205 9.183 ( 3.235 5.604 0.000) 6.471 9.219 ( 3.047 5.278 0.000) 6.094 10.004 ( 0.974 1.688 0.000) 1.949 10.013 ( 1.761 3.050 0.000) 3.522 ======================= Grid point 32 (16/96) ======================= q-point: ( 0.21 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.748 ( 9.165 12.112 0.000) 15.189 1.857 ( 1.354 18.172 0.000) 18.222 2.242 ( 0.297 13.177 0.000) 13.180 2.686 ( 0.047 1.720 0.000) 1.721 3.415 ( 19.882 15.496 0.000) 25.208 4.553 ( 3.179 1.255 0.000) 3.418 8.940 ( -1.866 -0.538 0.000) 1.942 8.987 ( 0.758 -4.495 0.000) 4.559 9.257 ( -1.805 7.931 0.000) 8.134 9.341 ( 3.820 5.538 0.000) 6.728 9.997 ( -1.915 0.046 0.000) 1.916 10.068 ( 3.162 -1.971 0.000) 3.726 ======================= Grid point 33 (17/96) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.955 ( -2.703 14.309 0.000) 14.562 2.046 ( 5.670 12.395 0.000) 13.631 2.316 ( -6.107 12.443 0.000) 13.861 2.668 ( -2.960 0.504 0.000) 3.003 3.874 ( 19.081 8.284 0.000) 20.802 4.628 ( 4.517 -0.996 0.000) 4.626 8.874 ( -1.399 -3.799 0.000) 4.048 8.977 ( 1.066 -3.679 0.000) 3.830 9.290 ( -3.218 7.864 0.000) 8.498 9.449 ( 1.453 6.961 0.000) 7.111 9.946 ( -3.322 -0.798 0.000) 3.417 10.082 ( 1.454 -2.766 0.000) 3.125 ======================= Grid point 34 (18/96) ======================= q-point: ( 0.36 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.996 ( -8.493 15.740 0.000) 17.885 2.258 ( 5.719 8.442 0.000) 10.197 2.302 ( -8.916 12.122 0.000) 15.047 2.604 ( -4.259 -0.288 0.000) 4.268 4.218 ( 13.872 2.435 0.000) 14.084 4.668 ( 3.244 -4.678 0.000) 5.692 8.838 ( 1.265 -2.857 0.000) 3.124 8.937 ( 0.084 -5.756 0.000) 5.757 9.308 ( -3.655 7.707 0.000) 8.530 9.531 ( -0.518 7.786 0.000) 7.803 9.877 ( -3.253 -0.965 0.000) 3.393 10.072 ( -0.449 -0.522 0.000) 0.689 ======================= Grid point 35 (19/96) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.999 ( -9.231 15.988 0.000) 18.461 2.280 ( -7.647 13.246 0.000) 15.295 2.370 ( -2.555 4.425 0.000) 5.109 2.547 ( -0.349 0.605 0.000) 0.699 4.362 ( 2.662 -4.610 0.000) 5.323 4.663 ( 3.425 -5.933 0.000) 6.851 8.848 ( 1.057 -1.831 0.000) 2.115 8.894 ( 3.292 -5.702 0.000) 6.584 9.316 ( -4.285 7.422 0.000) 8.570 9.571 ( -3.985 6.902 0.000) 7.969 9.838 ( -0.011 0.019 0.000) 0.022 10.065 ( -0.754 1.305 0.000) 1.507 ======================= Grid point 47 (20/96) ======================= q-point: ( 0.21 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.981 ( 6.422 11.122 0.000) 12.843 2.139 ( 5.793 10.034 0.000) 11.586 2.444 ( 3.816 6.610 0.000) 7.632 2.706 ( 0.086 0.149 0.000) 0.173 3.729 ( 9.227 15.982 0.000) 18.455 4.579 ( 0.619 1.072 0.000) 1.238 8.868 ( -2.633 -4.561 0.000) 5.267 8.948 ( -0.722 -1.251 0.000) 1.445 9.386 ( 2.317 4.014 0.000) 4.635 9.447 ( 3.344 5.792 0.000) 6.687 10.010 ( -1.605 -2.780 0.000) 3.211 10.015 ( 0.611 1.059 0.000) 1.223 ======================= Grid point 48 (21/96) ======================= q-point: ( 0.29 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.174 ( 2.979 9.488 0.000) 9.945 2.287 ( 1.833 9.613 0.000) 9.787 2.510 ( -3.061 6.430 0.000) 7.121 2.680 ( -3.484 0.113 0.000) 3.486 4.035 ( 11.520 8.177 0.000) 14.127 4.585 ( 1.075 -3.137 0.000) 3.317 8.800 ( -0.559 -3.094 0.000) 3.144 8.905 ( -1.380 -3.662 0.000) 3.913 9.433 ( -1.815 6.072 0.000) 6.338 9.561 ( 3.376 4.482 0.000) 5.611 9.930 ( -3.724 -2.083 0.000) 4.267 10.064 ( 1.912 1.981 0.000) 2.753 ======================= Grid point 49 (22/96) ======================= q-point: ( 0.36 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.251 ( -4.379 10.154 0.000) 11.058 2.431 ( 3.556 7.946 0.000) 8.705 2.499 ( -5.039 6.968 0.000) 8.600 2.594 ( -5.295 -1.024 0.000) 5.393 4.257 ( 6.987 2.034 0.000) 7.277 4.554 ( 1.394 -6.372 0.000) 6.522 8.792 ( 2.312 -1.734 0.000) 2.890 8.840 ( -1.428 -3.849 0.000) 4.105 9.454 ( -3.143 6.680 0.000) 7.383 9.643 ( 1.276 3.159 0.000) 3.407 9.855 ( -2.378 -1.223 0.000) 2.675 10.104 ( -0.611 3.284 0.000) 3.340 ======================= Grid point 62 (23/96) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.349 ( 4.315 7.473 0.000) 8.629 2.423 ( 2.912 5.044 0.000) 5.825 2.578 ( 0.598 1.036 0.000) 1.197 2.649 ( -2.028 -3.513 0.000) 4.056 4.200 ( 4.645 8.045 0.000) 9.290 4.508 ( -2.632 -4.559 0.000) 5.264 8.762 ( -0.325 -0.563 0.000) 0.650 8.834 ( -1.900 -3.292 0.000) 3.801 9.528 ( 2.074 3.593 0.000) 4.149 9.631 ( 1.357 2.350 0.000) 2.714 9.873 ( -1.599 -2.770 0.000) 3.199 10.129 ( 1.871 3.240 0.000) 3.741 ======================= Grid point 63 (24/96) ======================= q-point: ( 0.36 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.423 ( -4.164 7.212 0.000) 8.328 2.540 ( 0.721 -1.249 0.000) 1.443 2.544 ( -0.383 0.663 0.000) 0.765 2.584 ( -0.678 1.174 0.000) 1.355 4.311 ( -1.846 3.197 0.000) 3.692 4.434 ( 2.896 -5.016 0.000) 5.792 8.769 ( 0.312 -0.541 0.000) 0.625 8.788 ( 0.737 -1.277 0.000) 1.475 9.568 ( -2.654 4.596 0.000) 5.307 9.652 ( 1.067 -1.849 0.000) 2.135 9.842 ( 0.016 -0.028 0.000) 0.032 10.165 ( -1.117 1.935 0.000) 2.234 ======================= Grid point 211 (25/96) ======================= q-point: ( 0.00 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.475 ( -0.000 -0.000 21.554) 21.554 0.475 ( 0.000 0.000 21.554) 21.554 1.096 ( 0.000 0.000 51.075) 51.075 1.445 ( -0.000 -0.000 0.374) 0.374 1.445 ( 0.000 -0.000 0.374) 0.374 4.734 ( 0.000 -0.000 -11.973) 11.973 8.944 ( 0.000 0.000 0.933) 0.933 8.944 ( -0.000 -0.000 0.933) 0.933 9.043 ( 0.000 0.000 -0.097) 0.097 9.043 ( -0.000 0.000 -0.097) 0.097 9.700 ( -0.000 0.000 2.652) 2.652 9.800 ( 0.000 -0.000 1.969) 1.969 ======================= Grid point 212 (26/96) ======================= q-point: ( 0.07 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.639 ( 12.207 7.047 15.933) 21.273 0.825 ( 13.820 7.979 21.882) 27.083 1.339 ( 25.342 14.631 33.232) 44.279 1.495 ( 4.212 2.432 0.542) 4.894 1.704 ( 19.892 11.485 -0.021) 22.969 4.676 ( -4.673 -2.698 -11.497) 12.700 8.955 ( 0.968 0.559 0.936) 1.458 8.985 ( 2.712 1.566 1.524) 3.483 9.033 ( -0.163 -0.094 6.428) 6.431 9.060 ( 1.354 0.782 -0.142) 1.570 9.689 ( 2.658 1.535 0.228) 3.078 9.818 ( 1.777 1.026 1.892) 2.791 ======================= Grid point 213 (27/96) ======================= q-point: ( 0.14 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.953 ( 13.990 8.077 10.518) 19.277 1.086 ( 11.181 6.456 18.180) 22.298 1.627 ( 6.890 3.978 0.715) 7.988 2.045 ( 29.994 17.317 13.225) 37.073 2.225 ( 24.469 14.127 2.207) 28.341 4.550 ( -5.361 -3.095 -9.856) 11.639 8.988 ( 1.818 1.050 0.919) 2.292 8.998 ( -0.323 -0.186 7.889) 7.898 9.091 ( 3.441 1.986 0.696) 4.033 9.097 ( 1.644 0.949 -0.180) 1.907 9.811 ( 6.306 3.641 -3.028) 7.886 9.893 ( 4.951 2.858 1.122) 5.826 ======================= Grid point 214 (28/96) ======================= q-point: ( 0.21 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.248 ( 11.127 6.424 7.751) 15.005 1.355 ( 12.076 6.972 7.439) 15.804 1.790 ( 6.786 3.918 0.597) 7.858 2.501 ( 10.584 6.111 -1.990) 12.382 2.885 ( 30.362 17.529 10.798) 36.684 4.479 ( 0.123 0.071 -6.560) 6.562 8.997 ( 0.165 0.095 6.294) 6.297 9.033 ( 1.957 1.130 0.888) 2.428 9.132 ( 1.399 0.808 -0.107) 1.619 9.146 ( 1.306 0.754 0.251) 1.529 9.930 ( 3.729 2.153 -2.627) 5.044 10.027 ( 5.650 3.262 -0.961) 6.595 ======================= Grid point 215 (29/96) ======================= q-point: ( 0.29 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.450 ( 6.356 3.670 6.319) 9.685 1.630 ( 11.623 6.710 -0.989) 13.457 1.917 ( 3.851 2.223 0.402) 4.465 2.604 ( -0.310 -0.179 -5.950) 5.960 3.539 ( 26.023 15.024 10.193) 31.730 4.561 ( 6.261 3.615 -2.293) 7.584 8.997 ( -0.283 -0.163 4.504) 4.516 9.073 ( 1.386 0.800 0.912) 1.842 9.163 ( 1.389 0.802 0.008) 1.604 9.165 ( 0.746 0.431 0.882) 1.232 9.968 ( -0.665 -0.384 -3.375) 3.462 10.118 ( 1.933 1.116 -1.881) 2.919 ======================= Grid point 216 (30/96) ======================= q-point: ( 0.36 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.544 ( 2.133 1.232 5.668) 6.180 1.881 ( 9.977 5.760 -6.632) 13.294 1.959 ( -0.092 -0.053 0.453) 0.465 2.546 ( -3.740 -2.159 -8.476) 9.513 4.069 ( 19.921 11.502 9.600) 24.926 4.721 ( 6.665 3.848 1.342) 7.813 8.984 ( -0.708 -0.409 3.897) 3.981 9.098 ( 0.858 0.495 1.007) 1.413 9.186 ( 0.945 0.546 1.381) 1.760 9.195 ( 1.320 0.762 0.044) 1.525 9.908 ( -4.037 -2.331 -4.205) 6.278 10.115 ( -1.944 -1.122 -2.814) 3.600 ======================= Grid point 217 (31/96) ======================= q-point: ( 0.43 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.567 ( 0.227 0.131 5.428) 5.434 1.929 ( -1.899 -1.097 0.809) 2.337 2.079 ( 6.923 3.997 -9.786) 12.636 2.459 ( -3.519 -2.032 -10.229) 11.006 4.442 ( 12.035 6.949 9.111) 16.618 4.835 ( 2.965 1.712 3.981) 5.251 8.976 ( 0.106 0.061 4.717) 4.718 9.114 ( 0.522 0.301 1.078) 1.235 9.196 ( -0.291 -0.168 1.050) 1.103 9.220 ( 0.792 0.457 0.023) 0.915 9.819 ( -2.917 -1.684 -4.760) 5.831 10.050 ( -3.071 -1.773 -3.771) 5.176 ======================= Grid point 218 (32/96) ======================= q-point: (-0.50 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.568 ( 0.000 0.000 5.353) 5.353 1.902 ( 0.000 0.000 1.053) 1.053 2.170 ( 0.000 0.000 -10.751) 10.751 2.408 ( -0.000 -0.000 -10.960) 10.960 4.590 ( 0.000 0.000 8.784) 8.784 4.863 ( -0.000 -0.000 5.117) 5.117 8.981 ( 0.000 0.000 5.625) 5.625 9.120 ( 0.000 0.000 1.095) 1.095 9.186 ( 0.000 0.000 0.374) 0.374 9.229 ( 0.000 0.000 0.010) 0.010 9.787 ( 0.000 0.000 -4.827) 4.827 10.008 ( -0.000 -0.000 -4.256) 4.256 ======================= Grid point 227 (33/96) ======================= q-point: ( 0.07 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.914 ( 10.116 17.522 11.006) 23.032 1.009 ( 6.101 10.567 20.329) 23.709 1.598 ( 5.439 9.421 1.739) 11.016 1.870 ( 18.399 31.868 18.581) 41.223 2.039 ( 11.678 20.227 -0.894) 23.373 4.584 ( -3.377 -5.850 -10.448) 12.442 8.977 ( 0.482 0.834 4.057) 4.170 9.032 ( 0.814 1.410 5.581) 5.814 9.040 ( 2.003 3.470 0.245) 4.014 9.091 ( 0.962 1.666 -0.236) 1.938 9.770 ( 3.587 6.213 -2.822) 7.709 9.862 ( 1.996 3.457 1.679) 4.330 ======================= Grid point 228 (34/96) ======================= q-point: ( 0.14 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.230 ( 7.252 13.581 8.841) 17.754 1.272 ( 9.032 15.648 9.351) 20.344 1.796 ( 4.425 12.137 0.601) 12.932 2.391 ( 10.605 13.293 -0.484) 17.012 2.568 ( 26.731 21.735 11.973) 36.473 4.484 ( -1.901 -2.733 -8.000) 8.665 8.988 ( 0.439 0.260 6.424) 6.444 9.065 ( -0.541 4.868 1.050) 5.009 9.100 ( 0.705 -0.625 -0.145) 0.953 9.141 ( 3.356 3.271 0.408) 4.703 9.895 ( 3.563 4.244 -3.240) 6.419 9.958 ( 5.032 3.257 0.405) 6.008 ======================= Grid point 229 (35/96) ======================= q-point: ( 0.21 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.478 ( 6.248 12.020 2.234) 13.730 1.552 ( 5.495 15.096 5.170) 16.876 1.970 ( 0.649 13.081 0.430) 13.104 2.586 ( 3.308 2.413 -4.653) 6.198 3.229 ( 26.716 16.678 11.154) 33.411 4.495 ( 3.939 1.616 -4.134) 5.935 8.994 ( -0.291 -0.210 4.838) 4.851 9.077 ( 1.176 -2.994 0.133) 3.219 9.118 ( -0.933 5.406 1.135) 5.602 9.227 ( 1.280 5.475 0.095) 5.624 9.961 ( 1.063 0.635 -3.025) 3.269 10.067 ( 4.462 0.573 -1.228) 4.663 ======================= Grid point 230 (36/96) ======================= q-point: ( 0.29 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.650 ( -1.817 13.982 2.621) 14.342 1.785 ( 7.903 9.397 -2.684) 12.569 2.065 ( -3.997 11.680 0.352) 12.350 2.587 ( -2.646 -1.087 -7.381) 7.916 3.803 ( 23.324 11.414 10.219) 27.906 4.625 ( 7.438 2.520 -0.195) 7.856 8.965 ( -0.593 -3.061 2.605) 4.063 9.080 ( 2.527 -3.427 2.027) 4.716 9.152 ( -1.748 5.901 1.091) 6.251 9.295 ( -0.776 7.520 -0.092) 7.560 9.941 ( -2.323 -1.942 -3.510) 4.636 10.108 ( 1.056 -1.946 -2.248) 3.155 ======================= Grid point 231 (37/96) ======================= q-point: ( 0.36 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.709 ( -6.479 13.494 3.587) 15.392 2.001 ( 7.103 6.768 -7.916) 12.607 2.072 ( -7.182 10.207 0.364) 12.486 2.509 ( -4.286 -1.433 -9.448) 10.473 4.248 ( 17.566 6.447 9.626) 21.042 4.762 ( 6.082 0.111 2.949) 6.760 8.933 ( 1.228 -3.790 2.772) 4.853 9.077 ( 2.018 -5.498 1.992) 6.186 9.175 ( -2.343 5.885 1.181) 6.443 9.352 ( -2.125 8.896 -0.230) 9.150 9.862 ( -3.859 -1.690 -3.935) 5.765 10.074 ( -1.839 -2.532 -3.200) 4.476 ======================= Grid point 232 (38/96) ======================= q-point: ( 0.43 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.717 ( -7.600 13.198 3.817) 15.701 2.033 ( -7.651 10.586 0.709) 13.081 2.161 ( 3.264 4.419 -10.095) 11.493 2.428 ( -2.921 -1.090 -10.844) 11.283 4.516 ( 8.491 0.562 9.108) 12.465 4.823 ( 2.789 -2.902 4.995) 6.414 8.936 ( 2.503 -2.743 4.162) 5.578 9.043 ( 2.580 -8.100 0.872) 8.546 9.188 ( -2.963 5.826 1.180) 6.642 9.393 ( -4.140 9.529 -0.249) 10.392 9.800 ( -1.804 -0.188 -4.154) 4.533 10.020 ( -1.691 -0.771 -3.865) 4.289 ======================= Grid point 242 (39/96) ======================= q-point: ( 0.14 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.456 ( 6.484 11.231 2.027) 13.125 1.636 ( 9.723 16.840 5.217) 20.133 2.052 ( 7.093 12.285 0.479) 14.194 2.566 ( 2.964 5.134 -4.721) 7.578 3.050 ( 14.631 25.341 12.732) 31.912 4.468 ( 0.875 1.516 -4.931) 5.232 9.001 ( 0.488 0.846 6.098) 6.176 9.045 ( -2.028 -3.512 -0.761) 4.126 9.183 ( 2.872 4.974 0.020) 5.743 9.225 ( 3.099 5.367 0.596) 6.226 9.960 ( 1.066 1.847 -3.729) 4.296 10.015 ( 1.408 2.439 0.055) 2.817 ======================= Grid point 243 (40/96) ======================= q-point: ( 0.21 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.679 ( 7.157 9.518 -4.667) 12.790 1.901 ( 1.675 16.510 3.643) 16.990 2.247 ( 0.224 13.184 0.391) 13.192 2.612 ( -0.561 0.697 -6.684) 6.744 3.550 ( 19.031 15.840 11.609) 27.347 4.542 ( 4.549 2.858 -1.343) 5.537 8.959 ( -2.572 -4.133 0.331) 4.879 9.023 ( 1.308 -1.084 4.604) 4.908 9.253 ( -1.272 6.847 -0.512) 6.983 9.343 ( 3.288 5.806 0.238) 6.677 9.964 ( -1.621 0.210 -3.035) 3.447 10.051 ( 2.377 -2.079 -1.449) 3.474 ======================= Grid point 244 (41/96) ======================= q-point: ( 0.29 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.877 ( 4.435 8.703 -7.907) 12.567 2.044 ( -1.814 14.690 1.489) 14.877 2.321 ( -6.080 12.459 0.321) 13.867 2.575 ( -3.911 0.054 -8.241) 9.122 4.000 ( 18.352 8.693 10.719) 22.962 4.654 ( 5.791 0.322 2.003) 6.136 8.884 ( -0.621 -4.576 0.734) 4.676 9.018 ( 0.205 -3.122 3.959) 5.046 9.286 ( -2.500 6.699 -0.628) 7.178 9.444 ( 1.104 6.906 -0.337) 7.002 9.917 ( -3.634 0.001 -2.534) 4.431 10.055 ( 1.063 -2.772 -2.501) 3.882 ======================= Grid point 245 (42/96) ======================= q-point: ( 0.36 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.976 ( -3.983 11.084 -4.025) 12.447 2.173 ( 1.904 10.827 -4.921) 12.045 2.308 ( -8.732 12.049 0.389) 14.885 2.489 ( -5.309 -0.549 -9.835) 11.190 4.335 ( 13.262 2.719 9.981) 16.819 4.725 ( 4.067 -3.613 4.690) 7.182 8.857 ( 2.351 -3.505 1.498) 4.479 8.966 ( -0.949 -5.496 2.953) 6.311 9.306 ( -2.887 6.645 -0.446) 7.258 9.520 ( -0.537 7.229 -1.027) 7.322 9.849 ( -3.806 0.024 -2.347) 4.471 10.041 ( -0.851 0.044 -2.849) 2.974 ======================= Grid point 246 (43/96) ======================= q-point: ( 0.43 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.993 ( -6.901 11.952 -2.171) 13.971 2.272 ( -4.144 7.178 -6.541) 10.558 2.286 ( -7.539 13.057 0.151) 15.078 2.417 ( -0.183 0.317 -11.156) 11.162 4.475 ( 2.436 -4.219 9.531) 10.704 4.734 ( 3.146 -5.449 5.906) 8.629 8.890 ( 1.347 -2.333 3.845) 4.695 8.898 ( 3.233 -5.600 0.496) 6.485 9.315 ( -3.687 6.385 -0.291) 7.379 9.556 ( -3.477 6.023 -1.437) 7.101 9.810 ( -0.726 1.258 -2.273) 2.698 10.033 ( -1.217 2.109 -2.875) 3.768 ======================= Grid point 258 (44/96) ======================= q-point: ( 0.21 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.867 ( 5.359 9.282 -9.057) 14.032 2.154 ( 5.098 8.831 0.927) 10.239 2.450 ( 3.838 6.648 0.497) 7.693 2.619 ( -0.159 -0.275 -7.499) 7.506 3.866 ( 9.143 15.836 11.638) 21.675 4.600 ( 1.588 2.750 1.566) 3.540 8.866 ( -2.491 -4.315 -0.010) 4.982 9.001 ( -1.017 -1.761 4.692) 5.114 9.356 ( 1.852 3.207 -2.798) 4.641 9.458 ( 3.327 5.762 0.914) 6.716 9.988 ( 0.983 1.703 -2.369) 3.079 9.991 ( -1.763 -3.054 -1.706) 3.917 ======================= Grid point 259 (45/96) ======================= q-point: ( 0.29 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.040 ( 4.256 8.018 -11.472) 14.629 2.277 ( 0.573 7.977 -1.779) 8.193 2.516 ( -3.093 6.590 0.540) 7.299 2.585 ( -4.179 0.586 -7.191) 8.338 4.165 ( 10.954 8.203 11.017) 17.569 4.639 ( 2.013 -1.701 4.401) 5.129 8.804 ( -0.186 -3.048 0.536) 3.100 8.945 ( -1.897 -4.117 3.651) 5.821 9.400 ( -0.732 4.712 -3.200) 5.742 9.559 ( 2.139 4.506 -0.219) 4.993 9.912 ( -3.545 -1.671 -1.474) 4.187 10.046 ( 1.773 2.785 -1.639) 3.686 ======================= Grid point 260 (46/96) ======================= q-point: ( 0.36 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.150 ( -0.962 7.341 -10.634) 12.958 2.369 ( 2.481 6.214 -6.530) 9.349 2.490 ( -7.348 2.325 -5.682) 9.575 2.509 ( -5.062 6.534 0.739) 8.298 4.378 ( 6.523 2.028 10.162) 12.245 4.632 ( 1.698 -5.286 6.533) 8.573 8.802 ( 2.545 -1.927 0.961) 3.334 8.869 ( -1.602 -4.031 2.825) 5.177 9.428 ( -2.169 5.514 -2.677) 6.502 9.616 ( 0.607 2.048 -2.593) 3.360 9.850 ( -2.395 0.312 -0.100) 2.417 10.089 ( -0.961 4.048 -1.357) 4.377 ======================= Grid point 273 (47/96) ======================= q-point: ( 0.29 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.194 ( 4.060 7.031 -13.229) 15.522 2.373 ( 1.323 2.292 -7.547) 7.998 2.577 ( -0.920 -1.594 -3.203) 3.694 2.588 ( 0.718 1.244 0.990) 1.744 4.325 ( 4.349 7.533 10.505) 13.639 4.589 ( -1.902 -3.294 6.747) 7.745 8.771 ( -0.169 -0.292 0.894) 0.956 8.863 ( -2.223 -3.850 2.682) 5.192 9.480 ( 1.991 3.449 -4.400) 5.935 9.613 ( 0.392 0.679 -2.108) 2.249 9.875 ( -0.719 -1.246 0.583) 1.552 10.122 ( 2.024 3.505 -0.752) 4.116 ======================= Grid point 274 (48/96) ======================= q-point: ( 0.36 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.281 ( -3.147 5.451 -13.152) 14.580 2.388 ( 1.322 -2.289 -13.511) 13.767 2.550 ( -1.294 2.241 1.527) 3.004 2.595 ( -1.092 1.891 1.239) 2.511 4.427 ( -1.606 2.782 9.702) 10.220 4.532 ( 2.436 -4.219 8.165) 9.508 8.781 ( 0.135 -0.234 1.230) 1.259 8.811 ( 1.000 -1.731 2.103) 2.901 9.524 ( -2.281 3.951 -4.183) 6.190 9.607 ( 1.274 -2.206 -4.185) 4.900 9.865 ( -0.450 0.780 2.214) 2.390 10.159 ( -1.236 2.141 -0.560) 2.535 ======================= Grid point 422 (49/96) ======================= q-point: ( 0.00 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.887 ( 0.000 0.000 17.305) 17.305 0.887 ( 0.000 0.000 17.305) 17.305 1.439 ( 0.000 0.000 -1.360) 1.360 1.439 ( 0.000 -0.000 -1.360) 1.360 2.126 ( -0.000 -0.000 46.510) 46.510 4.367 ( 0.000 0.000 -22.662) 22.662 8.970 ( 0.000 -0.000 1.439) 1.439 8.970 ( 0.000 0.000 1.439) 1.439 9.038 ( 0.000 0.000 -0.394) 0.394 9.038 ( 0.000 -0.000 -0.394) 0.394 9.770 ( 0.000 0.000 3.691) 3.691 9.845 ( 0.000 -0.000 1.961) 1.961 ======================= Grid point 423 (50/96) ======================= q-point: ( 0.07 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.983 ( 7.818 4.514 15.544) 17.975 1.192 ( 19.097 11.026 14.575) 26.433 1.494 ( 4.543 2.623 -1.089) 5.358 1.684 ( 18.341 10.589 -1.561) 21.235 2.188 ( 7.603 4.390 42.298) 43.199 4.323 ( -3.436 -1.984 -21.870) 22.227 8.982 ( 0.980 0.566 1.468) 1.854 9.009 ( 3.074 1.775 1.375) 3.807 9.054 ( 1.266 0.731 -0.474) 1.537 9.104 ( 4.226 2.440 0.575) 4.914 9.736 ( -1.344 -0.776 3.810) 4.114 9.860 ( 1.409 0.813 1.600) 2.282 ======================= Grid point 424 (51/96) ======================= q-point: ( 0.14 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.207 ( 10.769 6.218 12.412) 17.570 1.506 ( 7.327 4.230 19.051) 20.845 1.632 ( 7.016 4.051 -0.668) 8.129 2.128 ( 17.680 10.208 -2.418) 20.558 2.568 ( 24.429 14.104 26.263) 38.541 4.245 ( -2.203 -1.272 -18.908) 19.078 9.014 ( 1.780 1.027 1.464) 2.523 9.081 ( 2.384 1.376 1.320) 3.053 9.089 ( 1.662 0.960 -0.537) 1.993 9.170 ( 1.310 0.756 3.948) 4.228 9.768 ( 3.944 2.277 0.179) 4.557 9.918 ( 3.635 2.098 0.618) 4.242 ======================= Grid point 425 (52/96) ======================= q-point: ( 0.21 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.444 ( 9.202 5.313 9.919) 14.536 1.607 ( 2.992 1.728 15.502) 15.883 1.794 ( 6.632 3.829 -0.554) 7.678 2.411 ( 6.085 3.513 -6.365) 9.481 3.191 ( 27.856 16.083 17.165) 36.459 4.274 ( 5.390 3.112 -12.800) 14.233 9.058 ( 1.927 1.113 1.398) 2.628 9.102 ( -0.317 -0.183 2.549) 2.575 9.127 ( 1.544 0.891 -0.417) 1.831 9.189 ( 0.737 0.425 3.619) 3.718 9.865 ( 3.302 1.906 -2.125) 4.365 10.011 ( 3.884 2.242 -1.104) 4.619 ======================= Grid point 426 (53/96) ======================= q-point: ( 0.29 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.613 ( 5.333 3.079 8.364) 10.386 1.677 ( 3.422 1.976 6.311) 7.446 1.917 ( 3.774 2.179 -0.620) 4.402 2.427 ( -3.428 -1.979 -10.576) 11.293 3.801 ( 24.653 14.234 13.664) 31.577 4.477 ( 11.154 6.440 -5.823) 14.134 9.085 ( -0.849 -0.490 3.262) 3.407 9.098 ( 1.466 0.847 1.412) 2.204 9.161 ( 1.480 0.854 -0.256) 1.728 9.209 ( 1.012 0.584 2.694) 2.936 9.889 ( -1.392 -0.804 -2.916) 3.329 10.067 ( 0.507 0.293 -3.137) 3.191 ======================= Grid point 427 (54/96) ======================= q-point: ( 0.36 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.691 ( 1.685 0.973 7.547) 7.794 1.766 ( 4.215 2.433 -2.565) 5.501 1.962 ( 0.206 0.119 -0.434) 0.495 2.311 ( -5.740 -3.314 -13.072) 14.656 4.308 ( 19.086 11.019 11.887) 25.040 4.734 ( 10.277 5.933 -0.566) 11.880 9.073 ( 0.013 0.008 4.156) 4.156 9.126 ( 0.979 0.565 1.556) 1.923 9.194 ( 1.303 0.752 -0.234) 1.522 9.235 ( 1.088 0.628 2.640) 2.924 9.809 ( -5.042 -2.911 -3.965) 7.044 10.032 ( -3.318 -1.916 -5.079) 6.362 ======================= Grid point 428 (55/96) ======================= q-point: ( 0.43 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.706 ( 0.024 0.014 7.128) 7.128 1.863 ( 3.923 2.265 -8.792) 9.891 1.942 ( -1.416 -0.818 0.131) 1.640 2.188 ( -4.611 -2.662 -14.308) 15.266 4.662 ( 11.322 6.537 10.644) 16.859 4.918 ( 5.297 3.058 3.089) 6.853 9.094 ( 1.922 1.110 5.795) 6.205 9.144 ( 0.572 0.330 1.672) 1.798 9.218 ( 0.760 0.439 -0.286) 0.923 9.242 ( -0.938 -0.542 3.355) 3.526 9.697 ( -3.686 -2.128 -5.790) 7.187 9.937 ( -4.172 -2.409 -6.810) 8.342 ======================= Grid point 429 (56/96) ======================= q-point: (-0.50 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.705 ( 0.000 0.000 6.958) 6.958 1.919 ( 0.000 0.000 -11.199) 11.199 1.921 ( 0.000 0.000 0.500) 0.500 2.125 ( -0.000 -0.000 -14.544) 14.544 4.800 ( 0.000 0.000 9.961) 9.961 4.976 ( 0.000 0.000 4.671) 4.671 9.139 ( -0.000 -0.000 8.742) 8.742 9.151 ( 0.000 0.000 1.702) 1.702 9.207 ( 0.000 0.000 1.801) 1.801 9.227 ( 0.000 0.000 -0.310) 0.310 9.658 ( 0.000 0.000 -6.658) 6.658 9.880 ( -0.000 -0.000 -7.720) 7.720 ======================= Grid point 438 (57/96) ======================= q-point: ( 0.07 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.178 ( 7.672 13.288 12.787) 19.973 1.432 ( 5.705 9.882 17.485) 20.879 1.612 ( 4.834 8.373 -0.277) 9.672 2.004 ( 10.762 18.640 -2.637) 21.685 2.417 ( 11.353 19.664 31.470) 38.806 4.262 ( -1.932 -3.346 -20.008) 20.377 9.017 ( 1.410 2.441 1.364) 3.132 9.053 ( 1.556 2.695 0.960) 3.257 9.085 ( 1.120 1.940 -0.239) 2.253 9.157 ( 1.066 1.846 3.195) 3.841 9.746 ( 1.502 2.601 1.353) 3.295 9.895 ( 1.609 2.787 0.983) 3.365 ======================= Grid point 439 (58/96) ======================= q-point: ( 0.14 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.445 ( 4.865 15.642 10.272) 19.336 1.565 ( 4.335 2.243 17.137) 17.819 1.800 ( 4.353 12.100 -0.483) 12.868 2.324 ( 8.893 8.069 -5.125) 13.056 2.927 ( 21.806 20.957 20.409) 36.486 4.236 ( 1.629 1.079 -15.433) 15.556 9.066 ( 0.677 3.611 1.284) 3.892 9.081 ( 0.201 -1.050 0.299) 1.110 9.140 ( 1.382 3.655 1.266) 4.107 9.177 ( 1.180 0.289 3.972) 4.153 9.828 ( 3.357 3.323 -1.602) 4.988 9.967 ( 3.277 2.506 -0.191) 4.130 ======================= Grid point 440 (59/96) ======================= q-point: ( 0.21 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.593 ( 2.480 3.764 8.977) 10.045 1.704 ( 2.111 12.573 8.496) 15.320 1.973 ( 0.509 13.099 -0.330) 13.113 2.440 ( 0.567 -0.311 -9.015) 9.038 3.520 ( 24.140 16.379 15.379) 32.978 4.360 ( 8.577 5.427 -8.708) 13.373 9.051 ( -0.289 -3.720 0.902) 3.839 9.108 ( 0.224 2.753 1.331) 3.066 9.174 ( 0.091 2.499 4.070) 4.777 9.231 ( 1.058 4.773 0.345) 4.901 9.889 ( 0.855 0.419 -2.566) 2.737 10.036 ( 2.651 -0.048 -2.030) 3.339 ======================= Grid point 441 (60/96) ======================= q-point: ( 0.29 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.670 ( 3.472 1.978 1.524) 4.277 1.829 ( -2.261 12.638 5.370) 13.917 2.067 ( -4.025 11.636 -0.483) 12.322 2.379 ( -4.960 -2.063 -11.649) 12.828 4.059 ( 21.805 11.482 12.959) 27.843 4.598 ( 11.336 5.353 -2.761) 12.837 9.010 ( 0.947 -5.128 1.839) 5.529 9.130 ( 0.410 0.711 1.202) 1.456 9.202 ( -0.245 3.081 3.953) 5.018 9.292 ( -0.769 6.925 -0.218) 6.971 9.861 ( -3.224 -1.481 -3.070) 4.691 10.043 ( -0.460 -2.543 -3.983) 4.748 ======================= Grid point 442 (61/96) ======================= q-point: ( 0.36 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.744 ( 1.752 2.945 -2.182) 4.062 1.895 ( -2.798 10.722 0.952) 11.122 2.074 ( -7.044 9.851 -0.710) 12.131 2.261 ( -6.353 -1.202 -12.692) 14.243 4.483 ( 16.428 6.525 11.470) 21.072 4.818 ( 8.904 2.250 1.628) 9.327 8.993 ( 3.381 -5.298 2.821) 6.889 9.133 ( 0.477 -2.355 1.976) 3.111 9.221 ( -1.199 3.040 3.504) 4.791 9.344 ( -1.752 7.910 -0.551) 8.121 9.770 ( -5.409 -0.273 -3.862) 6.652 9.979 ( -3.219 -2.718 -5.668) 7.062 ======================= Grid point 443 (62/96) ======================= q-point: ( 0.43 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.779 ( -1.984 4.778 -0.549) 5.202 1.961 ( -1.516 9.124 -4.756) 10.400 2.040 ( -6.990 8.670 -1.667) 11.261 2.156 ( -4.406 0.284 -12.412) 13.173 4.734 ( 7.697 0.696 10.356) 12.921 4.932 ( 3.919 -1.658 4.442) 6.151 9.024 ( 5.973 -4.796 3.242) 8.318 9.092 ( -0.026 -6.249 3.773) 7.300 9.226 ( -2.027 3.773 2.609) 5.015 9.383 ( -3.360 8.224 -0.693) 8.911 9.698 ( -3.432 2.080 -4.685) 6.169 9.906 ( -2.771 -0.049 -6.779) 7.323 ======================= Grid point 453 (63/96) ======================= q-point: ( 0.14 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.589 ( 1.081 1.872 10.614) 10.831 1.771 ( 8.182 14.171 6.999) 17.798 2.058 ( 7.140 12.367 -0.112) 14.281 2.425 ( 1.238 2.144 -8.659) 9.006 3.373 ( 13.478 23.345 16.916) 31.825 4.310 ( 3.777 6.541 -10.034) 12.560 9.035 ( -1.832 -3.174 0.137) 3.667 9.123 ( 0.697 1.207 2.932) 3.246 9.196 ( 1.504 2.604 2.335) 3.807 9.241 ( 3.229 5.592 0.845) 6.513 9.882 ( 1.030 1.784 -2.579) 3.301 10.005 ( 0.640 1.108 -1.258) 1.794 ======================= Grid point 454 (64/96) ======================= q-point: ( 0.21 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.637 ( 2.403 2.244 2.047) 3.873 1.985 ( 1.118 12.049 3.479) 12.591 2.252 ( -0.115 13.362 -0.097) 13.363 2.429 ( -1.485 0.102 -9.892) 10.004 3.838 ( 17.284 15.710 14.537) 27.511 4.483 ( 7.869 6.714 -4.471) 11.269 8.965 ( -1.427 -4.624 0.879) 4.919 9.130 ( -0.095 -0.242 2.730) 2.742 9.246 ( 0.010 3.722 1.341) 3.957 9.353 ( 2.115 6.703 0.619) 7.056 9.894 ( -1.115 0.663 -2.758) 3.048 10.008 ( 0.654 -2.642 -2.738) 3.861 ======================= Grid point 455 (65/96) ======================= q-point: ( 0.29 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.715 ( 3.512 2.811 -5.184) 6.863 2.065 ( -2.541 8.086 -1.878) 8.681 2.315 ( -6.721 11.561 -1.615) 13.470 2.389 ( -5.554 4.681 -6.214) 9.559 4.260 ( 16.793 8.966 12.782) 22.930 4.684 ( 8.743 3.274 0.234) 9.339 8.908 ( 0.526 -4.771 1.834) 5.139 9.112 ( -0.754 -2.476 2.867) 3.862 9.269 ( -0.783 2.816 0.307) 2.939 9.439 ( -0.123 7.028 -0.207) 7.032 9.860 ( -3.988 1.660 -2.330) 4.908 9.984 ( -0.159 -2.511 -4.115) 4.823 ======================= Grid point 456 (66/96) ======================= q-point: ( 0.36 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.801 ( 2.486 3.066 -9.302) 10.105 2.086 ( -1.626 4.284 -8.809) 9.930 2.272 ( -9.595 9.904 -2.453) 14.006 2.341 ( -8.398 10.250 -0.763) 13.272 4.573 ( 11.954 2.891 11.353) 16.738 4.825 ( 5.991 -1.349 3.864) 7.255 8.896 ( 3.342 -4.079 2.460) 5.819 9.051 ( -1.521 -5.182 4.011) 6.727 9.285 ( -0.771 3.084 -1.088) 3.360 9.489 ( -1.232 5.465 -2.006) 5.950 9.801 ( -4.782 2.730 -1.581) 5.729 9.961 ( -2.046 1.938 -4.520) 5.327 ======================= Grid point 457 (67/96) ======================= q-point: ( 0.43 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.842 ( -1.150 1.992 -10.068) 10.327 2.097 ( -0.017 0.030 -15.624) 15.624 2.228 ( -8.262 14.310 1.835) 16.626 2.317 ( -7.206 12.481 0.376) 14.416 4.700 ( 2.124 -3.678 10.536) 11.360 4.865 ( 2.610 -4.521 5.524) 7.601 8.922 ( 3.015 -5.222 1.987) 6.349 8.992 ( 1.948 -3.374 5.139) 6.449 9.295 ( -1.764 3.055 -1.393) 3.792 9.506 ( -1.962 3.399 -3.478) 5.244 9.770 ( -2.825 4.893 -0.876) 5.718 9.953 ( -2.598 4.500 -4.486) 6.864 ======================= Grid point 469 (68/96) ======================= q-point: ( 0.21 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.697 ( 2.226 3.856 -5.010) 6.703 2.123 ( 1.309 2.268 -5.852) 6.411 2.460 ( 3.964 6.866 0.366) 7.936 2.464 ( 1.653 2.863 -4.792) 5.822 4.147 ( 8.736 15.132 13.673) 22.186 4.618 ( 3.782 6.550 -0.430) 7.576 8.880 ( -2.099 -3.635 1.574) 4.483 9.106 ( -1.333 -2.308 3.042) 4.044 9.292 ( 0.669 1.159 -1.356) 1.905 9.476 ( 2.974 5.151 0.458) 5.965 9.931 ( 1.572 2.723 -2.641) 4.106 9.943 ( -1.795 -3.109 -2.686) 4.484 ======================= Grid point 470 (69/96) ======================= q-point: ( 0.29 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.787 ( 3.085 4.374 -10.330) 11.634 2.124 ( -1.208 -0.424 -12.866) 12.929 2.494 ( -3.930 6.782 -0.087) 7.838 2.535 ( -2.646 6.582 0.620) 7.121 4.426 ( 9.800 7.817 12.397) 17.630 4.730 ( 4.111 1.354 3.299) 5.442 8.831 ( 0.337 -2.639 2.312) 3.524 9.036 ( -2.400 -4.860 3.683) 6.553 9.314 ( 0.373 2.156 -3.674) 4.276 9.540 ( -0.147 2.345 -2.173) 3.200 9.885 ( -2.444 0.516 -0.580) 2.565 10.002 ( 1.270 4.496 -2.438) 5.270 ======================= Grid point 471 (70/96) ======================= q-point: ( 0.36 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.872 ( 1.052 4.068 -13.128) 13.784 2.093 ( -0.713 -1.585 -16.285) 16.377 2.480 ( -6.012 8.907 2.117) 10.952 2.532 ( -4.959 7.478 1.032) 9.032 4.615 ( 5.626 1.648 10.997) 12.462 4.777 ( 2.480 -3.097 6.095) 7.272 8.837 ( 2.418 -1.766 2.606) 3.969 8.948 ( -1.256 -4.726 4.112) 6.389 9.345 ( -0.343 3.052 -4.557) 5.495 9.532 ( -0.579 -0.874 -5.395) 5.496 9.868 ( -2.491 3.886 1.957) 5.014 10.050 ( -1.732 5.624 -2.309) 6.322 ======================= Grid point 484 (71/96) ======================= q-point: ( 0.29 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.884 ( 2.926 5.067 -14.010) 15.183 2.091 ( -1.465 -2.538 -16.354) 16.614 2.580 ( 1.219 2.111 1.917) 3.101 2.614 ( 0.995 1.723 1.239) 2.344 4.569 ( 3.623 6.276 11.272) 13.400 4.737 ( -0.377 -0.654 6.222) 6.267 8.808 ( 0.141 0.245 2.808) 2.823 8.936 ( -2.712 -4.698 3.717) 6.576 9.364 ( 1.625 2.815 -5.789) 6.640 9.530 ( -1.445 -2.503 -5.628) 6.327 9.911 ( 1.012 1.753 2.793) 3.449 10.097 ( 2.280 3.949 -1.671) 4.856 ======================= Grid point 485 (72/96) ======================= q-point: ( 0.36 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.955 ( -2.171 3.760 -15.626) 16.218 2.048 ( 1.520 -2.632 -17.009) 17.279 2.597 ( -1.660 2.875 2.499) 4.155 2.627 ( -1.280 2.217 1.523) 2.979 4.649 ( -1.039 1.800 10.060) 10.273 4.716 ( 1.516 -2.626 8.037) 8.590 8.825 ( -0.218 0.378 3.027) 3.059 8.873 ( 1.481 -2.564 3.635) 4.689 9.406 ( -1.689 2.926 -6.343) 7.186 9.487 ( 1.600 -2.772 -6.786) 7.503 9.932 ( -1.111 1.924 3.890) 4.479 10.138 ( -1.452 2.515 -1.533) 3.284 ======================= Grid point 633 (73/96) ======================= q-point: ( 0.00 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.191 ( 0.000 0.000 11.489) 11.489 1.191 ( -0.000 -0.000 11.489) 11.489 1.369 ( 0.000 0.000 -5.649) 5.649 1.369 ( 0.000 0.000 -5.649) 5.649 3.036 ( 0.000 0.000 39.867) 39.867 3.792 ( -0.000 -0.000 -31.868) 31.868 9.001 ( 0.000 -0.000 1.364) 1.364 9.001 ( 0.000 0.000 1.364) 1.364 9.025 ( 0.000 0.000 -0.897) 0.897 9.025 ( -0.000 0.000 -0.897) 0.897 9.839 ( 0.000 -0.000 2.546) 2.546 9.868 ( 0.000 -0.000 -0.067) 0.067 ======================= Grid point 634 (74/96) ======================= q-point: ( 0.07 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.263 ( 5.959 3.440 10.758) 12.770 1.431 ( 5.064 2.924 -5.261) 7.866 1.453 ( 18.613 10.746 10.122) 23.757 1.614 ( 18.154 10.481 -5.313) 21.626 3.041 ( 1.366 0.789 38.057) 38.089 3.768 ( -1.580 -0.912 -30.763) 30.817 9.013 ( 1.039 0.600 1.434) 1.869 9.039 ( 1.146 0.662 -0.990) 1.653 9.045 ( 3.354 1.937 1.950) 4.337 9.086 ( 4.462 2.576 -1.764) 5.446 9.821 ( -1.021 -0.589 3.589) 3.778 9.870 ( 0.440 0.254 -0.891) 1.026 ======================= Grid point 635 (75/96) ======================= q-point: ( 0.14 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.439 ( 8.668 5.005 9.144) 13.557 1.582 ( 7.521 4.342 -4.422) 9.746 1.837 ( 11.510 6.645 12.382) 18.164 2.034 ( 15.216 8.785 -6.635) 18.781 3.173 ( 11.666 6.735 29.722) 32.632 3.768 ( 3.051 1.762 -26.102) 26.339 9.046 ( 1.734 1.001 1.466) 2.481 9.073 ( 1.693 0.977 -1.052) 2.219 9.121 ( 2.592 1.496 2.624) 3.980 9.179 ( 2.832 1.635 -2.248) 3.968 9.826 ( 1.786 1.031 4.257) 4.731 9.901 ( 2.407 1.390 -2.346) 3.637 ======================= Grid point 636 (76/96) ======================= q-point: ( 0.21 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.633 ( 7.668 4.427 7.599) 11.668 1.752 ( 6.834 3.946 -3.711) 8.720 1.948 ( -0.484 -0.279 15.659) 15.669 2.231 ( 1.832 1.058 -10.997) 11.199 3.573 ( 21.412 12.362 18.641) 30.964 3.955 ( 13.043 7.530 -17.193) 22.856 9.088 ( 1.850 1.068 1.363) 2.533 9.113 ( 1.714 0.989 -0.924) 2.184 9.153 ( 0.359 0.207 2.912) 2.942 9.212 ( 0.465 0.268 -1.805) 1.883 9.882 ( 2.188 1.263 3.448) 4.274 9.964 ( 2.397 1.384 -3.456) 4.427 ======================= Grid point 637 (77/96) ======================= q-point: ( 0.29 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.774 ( 4.432 2.559 6.566) 8.325 1.878 ( 3.912 2.259 -3.305) 5.597 1.888 ( -3.613 -2.086 12.783) 13.446 2.172 ( -5.482 -3.165 -13.391) 14.812 4.081 ( 21.648 12.499 12.515) 27.955 4.315 ( 16.885 9.749 -9.562) 21.716 9.128 ( 1.539 0.888 1.318) 2.212 9.150 ( 1.547 0.893 -0.800) 1.957 9.153 ( -0.048 -0.027 3.531) 3.532 9.221 ( 0.554 0.320 -2.051) 2.148 9.890 ( -1.864 -1.076 3.106) 3.779 9.983 ( -1.170 -0.676 -4.716) 4.906 ======================= Grid point 638 (78/96) ======================= q-point: ( 0.36 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.812 ( -2.630 -1.518 6.932) 7.568 1.836 ( 1.153 0.665 5.901) 6.049 1.929 ( 0.639 0.369 -2.891) 2.983 2.030 ( -6.135 -3.542 -12.838) 14.662 4.534 ( 17.245 9.956 9.172) 21.923 4.681 ( 14.195 8.195 -4.715) 17.055 9.158 ( 1.125 0.649 1.439) 1.938 9.162 ( 1.056 0.610 4.463) 4.626 9.183 ( 1.241 0.717 -0.847) 1.665 9.244 ( 1.541 0.890 -2.409) 2.995 9.796 ( -5.969 -3.446 3.220) 7.607 9.909 ( -5.030 -2.904 -6.301) 8.570 ======================= Grid point 639 (79/96) ======================= q-point: ( 0.43 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.774 ( -0.718 -0.415 1.120) 1.393 1.841 ( -0.337 -0.194 5.354) 5.368 1.909 ( -4.099 -2.366 -10.997) 11.972 1.922 ( -0.842 -0.486 -2.280) 2.479 4.852 ( 9.956 5.748 6.886) 13.401 4.940 ( 7.856 4.536 -1.414) 9.181 9.179 ( 0.640 0.370 1.560) 1.726 9.206 ( 0.707 0.408 -0.938) 1.244 9.206 ( 2.767 1.597 4.985) 5.921 9.285 ( 1.576 0.910 -1.155) 2.156 9.646 ( -6.097 -3.520 2.602) 7.505 9.777 ( -5.660 -3.268 -7.970) 10.307 ======================= Grid point 640 (80/96) ======================= q-point: (-0.50 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.770 ( 0.000 0.000 -1.815) 1.815 1.835 ( 0.000 0.000 5.077) 5.077 1.857 ( -0.000 -0.000 -9.573) 9.573 1.909 ( 0.000 0.000 -1.918) 1.918 4.971 ( -0.000 -0.000 5.747) 5.747 5.030 ( 0.000 0.000 0.100) 0.100 9.186 ( 0.000 0.000 1.595) 1.595 9.214 ( 0.000 0.000 -0.971) 0.971 9.264 ( -0.000 -0.000 3.344) 3.344 9.292 ( 0.000 0.000 3.321) 3.321 9.568 ( 0.000 0.000 0.626) 0.626 9.698 ( -0.000 -0.000 -8.995) 8.995 ======================= Grid point 649 (81/96) ======================= q-point: ( 0.07 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.415 ( 6.165 10.677 9.305) 15.447 1.558 ( 5.412 9.374 -4.190) 11.607 1.739 ( 8.050 13.943 10.941) 19.465 1.916 ( 9.660 16.732 -6.037) 20.242 3.109 ( 4.396 7.613 32.872) 34.028 3.755 ( 0.428 0.741 -27.889) 27.902 9.044 ( 1.244 2.155 1.163) 2.746 9.063 ( 1.015 1.759 -0.504) 2.093 9.101 ( 1.792 3.104 2.231) 4.222 9.156 ( 2.022 3.502 -2.273) 4.638 9.816 ( 0.445 0.770 4.278) 4.370 9.887 ( 0.969 1.678 -1.958) 2.754 ======================= Grid point 650 (82/96) ======================= q-point: ( 0.14 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.638 ( 4.426 13.628 7.682) 16.259 1.757 ( 3.988 11.977 -3.579) 13.121 1.917 ( 3.753 1.249 14.887) 15.403 2.174 ( 6.247 4.216 -9.475) 12.107 3.386 ( 14.543 14.513 22.618) 30.556 3.850 ( 7.215 7.030 -20.818) 23.127 9.085 ( 0.606 1.579 0.686) 1.825 9.090 ( 0.741 0.722 0.292) 1.076 9.164 ( 1.686 3.461 1.669) 4.196 9.205 ( 1.164 1.218 -1.402) 2.192 9.856 ( 2.089 1.812 3.672) 4.596 9.935 ( 2.213 1.678 -2.949) 4.051 ======================= Grid point 651 (83/96) ======================= q-point: ( 0.21 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.796 ( -0.173 7.350 8.035) 10.891 1.876 ( 0.240 4.296 3.756) 5.711 1.978 ( -0.437 9.590 5.195) 10.915 2.214 ( -1.670 -1.682 -12.305) 12.531 3.843 ( 20.030 14.179 14.949) 28.736 4.135 ( 14.498 10.160 -12.559) 21.706 9.079 ( 0.610 -3.143 1.796) 3.671 9.112 ( 0.498 0.213 -0.903) 1.053 9.225 ( -0.023 5.572 1.020) 5.664 9.239 ( 0.440 3.505 -0.120) 3.534 9.893 ( 0.472 -0.343 2.884) 2.942 9.973 ( 1.042 -0.707 -3.882) 4.081 ======================= Grid point 652 (84/96) ======================= q-point: ( 0.29 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.779 ( -0.998 -3.436 8.337) 9.073 1.915 ( -2.109 4.454 -0.156) 4.930 2.046 ( -5.130 11.719 -0.150) 12.794 2.142 ( -6.217 3.396 -9.189) 11.603 4.313 ( 19.313 10.555 10.740) 24.490 4.499 ( 15.706 8.393 -6.777) 19.053 9.059 ( 1.936 -4.688 2.904) 5.845 9.118 ( 1.406 -1.842 -2.136) 3.152 9.270 ( -0.893 5.317 1.794) 5.682 9.287 ( -1.071 6.686 -0.258) 6.776 9.853 ( -3.509 -2.107 2.635) 4.867 9.944 ( -2.316 -2.903 -5.103) 6.311 ======================= Grid point 653 (85/96) ======================= q-point: ( 0.36 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.738 ( 0.194 -2.989 3.745) 4.796 1.885 ( -2.216 -1.642 -7.208) 7.717 2.035 ( -8.372 10.933 0.817) 13.794 2.089 ( -8.061 10.067 -3.267) 13.304 4.694 ( 14.604 5.921 8.032) 17.687 4.808 ( 12.013 4.429 -2.835) 13.113 9.059 ( 3.933 -5.215 3.540) 7.430 9.127 ( 2.361 -2.977 -2.278) 4.430 9.298 ( -1.138 3.634 2.766) 4.706 9.331 ( -1.735 6.816 -0.778) 7.076 9.747 ( -6.583 -0.423 2.298) 6.985 9.848 ( -5.167 -2.263 -6.344) 8.489 ======================= Grid point 654 (86/96) ======================= q-point: ( 0.43 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.726 ( 1.064 -1.991 -0.043) 2.258 1.836 ( -0.635 -2.150 -9.723) 9.978 2.003 ( -8.202 12.208 2.982) 15.006 2.057 ( -7.620 11.155 -1.835) 13.633 4.913 ( 6.611 0.407 6.140) 9.031 4.980 ( 5.287 -0.422 -0.232) 5.309 9.086 ( 4.920 -5.636 2.931) 8.035 9.132 ( 2.223 -4.190 -0.697) 4.794 9.304 ( -0.561 0.828 3.973) 4.097 9.361 ( -1.916 5.204 -1.587) 5.768 9.652 ( -5.766 4.386 1.364) 7.371 9.753 ( -4.646 1.951 -7.102) 8.708 ======================= Grid point 664 (87/96) ======================= q-point: ( 0.14 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.851 ( -1.752 -3.035 11.694) 12.208 1.901 ( 5.860 10.150 4.510) 12.558 2.050 ( 6.985 12.098 0.449) 13.977 2.208 ( 0.340 0.589 -11.583) 11.603 3.729 ( 11.103 19.231 16.497) 27.663 4.054 ( 7.570 13.111 -14.117) 20.700 9.059 ( -1.362 -2.358 2.097) 3.437 9.111 ( -0.201 -0.349 -2.339) 2.373 9.253 ( 2.082 3.606 1.404) 4.394 9.257 ( 3.120 5.404 0.537) 6.263 9.884 ( 0.371 0.642 2.597) 2.701 9.956 ( 0.095 0.165 -3.356) 3.361 ======================= Grid point 665 (88/96) ======================= q-point: ( 0.21 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.772 ( -2.176 -3.915 9.882) 10.850 1.975 ( -2.071 -0.410 -6.468) 6.804 2.240 ( -0.562 12.976 -0.436) 12.995 2.269 ( 0.767 9.149 -3.562) 9.848 4.126 ( 14.910 14.159 12.377) 24.000 4.347 ( 11.678 10.975 -8.541) 18.160 9.008 ( -0.738 -3.685 3.250) 4.969 9.092 ( -0.351 -1.760 -4.150) 4.522 9.321 ( 0.613 4.203 3.470) 5.484 9.359 ( 1.320 6.538 -0.367) 6.680 9.877 ( -1.055 -0.649 1.426) 1.889 9.938 ( -0.660 -2.543 -3.647) 4.495 ======================= Grid point 666 (89/96) ======================= q-point: ( 0.29 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.714 ( -0.624 -2.646 5.295) 5.952 1.904 ( -2.741 -3.313 -11.050) 11.857 2.304 ( -6.507 13.074 1.441) 14.674 2.338 ( -5.727 12.204 -1.384) 13.552 4.499 ( 14.719 8.230 9.416) 19.315 4.644 ( 11.986 6.244 -4.295) 14.181 8.970 ( 1.051 -3.992 4.190) 5.882 9.073 ( 0.648 -2.512 -4.852) 5.502 9.344 ( -1.267 1.451 4.983) 5.343 9.421 ( -1.122 5.149 -2.017) 5.643 9.841 ( -3.416 1.525 0.853) 3.837 9.895 ( -1.854 -0.927 -3.896) 4.414 ======================= Grid point 667 (90/96) ======================= q-point: ( 0.36 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.697 ( 0.593 -1.173 0.690) 1.485 1.833 ( -1.692 -2.673 -11.308) 11.742 2.296 ( -9.059 13.460 2.259) 16.381 2.331 ( -8.574 12.532 -0.943) 15.214 4.773 ( 10.271 2.319 7.085) 12.691 4.859 ( 8.219 0.847 -1.039) 8.328 8.968 ( 3.062 -3.787 4.331) 6.518 9.060 ( 1.268 -3.583 -3.404) 5.102 9.319 ( -1.267 -0.866 4.670) 4.916 9.424 ( -1.851 1.411 -4.345) 4.930 9.801 ( -4.808 5.338 1.693) 7.381 9.866 ( -3.674 4.689 -4.144) 7.257 ======================= Grid point 668 (91/96) ======================= q-point: ( 0.43 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.699 ( 0.463 -0.802 -1.543) 1.799 1.802 ( 0.693 -1.200 -10.768) 10.856 2.281 ( -8.269 14.322 2.448) 16.718 2.316 ( -7.753 13.428 -0.746) 15.524 4.878 ( 2.021 -3.500 5.850) 7.110 4.934 ( 2.186 -3.787 0.547) 4.407 8.984 ( 2.595 -4.494 3.674) 6.359 9.048 ( 2.197 -3.805 -1.047) 4.516 9.302 ( 0.204 -0.353 3.380) 3.404 9.409 ( -0.121 0.209 -5.498) 5.503 9.786 ( -4.557 7.892 2.348) 9.411 9.858 ( -4.242 7.348 -4.221) 9.477 ======================= Grid point 680 (92/96) ======================= q-point: ( 0.21 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.706 ( -1.374 -2.379 5.954) 6.557 1.911 ( -2.376 -4.116 -11.992) 12.899 2.440 ( 4.052 7.018 0.904) 8.154 2.459 ( 4.160 7.206 -0.583) 8.341 4.403 ( 7.762 13.445 10.197) 18.574 4.564 ( 6.063 10.502 -5.030) 13.129 8.942 ( -1.558 -2.699 4.407) 5.398 9.055 ( -1.036 -1.794 -5.660) 6.027 9.364 ( -0.145 -0.250 5.755) 5.763 9.458 ( 1.523 2.637 -2.627) 4.022 9.890 ( 1.150 1.993 -1.072) 2.538 9.894 ( -0.628 -1.088 -1.526) 1.977 ======================= Grid point 681 (93/96) ======================= q-point: ( 0.29 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.682 ( 0.128 -0.554 1.326) 1.443 1.837 ( -1.927 -3.095 -12.263) 12.794 2.519 ( -2.904 7.910 1.519) 8.562 2.540 ( -2.648 7.374 -0.327) 7.842 4.646 ( 8.307 6.589 7.991) 13.277 4.750 ( 6.377 4.377 -1.806) 7.943 8.908 ( 0.506 -2.110 4.944) 5.399 9.022 ( -0.348 -2.538 -4.823) 5.461 9.330 ( -1.516 -2.016 5.109) 5.697 9.455 ( -1.817 -1.561 -5.663) 6.149 9.899 ( -0.726 4.002 1.778) 4.439 9.948 ( 0.455 5.357 -2.603) 5.973 ======================= Grid point 682 (94/96) ======================= q-point: ( 0.36 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.691 ( 0.474 0.518 -2.549) 2.644 1.787 ( -0.666 -2.004 -11.016) 11.217 2.524 ( -5.420 8.862 1.702) 10.527 2.543 ( -5.106 8.179 -0.087) 9.642 4.799 ( 4.610 0.647 5.991) 7.587 4.854 ( 3.517 -0.974 0.770) 3.730 8.914 ( 1.846 -1.628 4.424) 5.062 8.991 ( 0.318 -3.267 -1.409) 3.572 9.298 ( -0.310 -0.657 1.755) 1.900 9.401 ( -0.886 -2.678 -6.535) 7.117 9.926 ( -2.597 6.211 3.255) 7.478 9.994 ( -2.368 6.626 -3.056) 7.671 ======================= Grid point 695 (95/96) ======================= q-point: ( 0.29 0.29 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.690 ( 0.683 1.182 -2.912) 3.217 1.783 ( -1.332 -2.307 -11.073) 11.389 2.619 ( 1.445 2.503 1.373) 3.200 2.631 ( 1.279 2.215 0.213) 2.567 4.758 ( 2.572 4.454 6.154) 8.020 4.816 ( 1.209 2.094 0.736) 2.528 8.889 ( 0.113 0.195 4.723) 4.728 8.972 ( -1.420 -2.460 -1.617) 3.269 9.298 ( -0.366 -0.633 1.101) 1.321 9.395 ( -2.062 -3.572 -6.581) 7.767 9.981 ( 1.997 3.459 3.551) 5.344 10.050 ( 2.335 4.044 -2.818) 5.454 ======================= Grid point 696 (96/96) ======================= q-point: ( 0.36 0.29 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.711 ( -0.743 1.287 -5.753) 5.942 1.752 ( 0.846 -1.465 -9.315) 9.467 2.640 ( -1.623 2.811 1.337) 3.511 2.649 ( -1.471 2.548 0.433) 2.974 4.810 ( -0.297 0.515 4.634) 4.672 4.832 ( 0.566 -0.981 2.493) 2.738 8.902 ( -0.166 0.287 3.707) 3.722 8.937 ( 1.112 -1.926 1.266) 2.559 9.301 ( -0.487 0.844 -2.185) 2.392 9.347 ( 1.311 -2.270 -5.675) 6.251 10.018 ( -1.469 2.544 4.011) 4.971 10.092 ( -1.562 2.705 -2.913) 4.271 =================== 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/16464 10.0 4355.533 4355.533 2623.612 -0.000 0.000 -0.000 3/16464 20.0 1826.732 1826.732 1724.047 -0.000 0.000 -0.000 3/16464 30.0 1183.866 1183.866 1144.529 -0.000 0.000 -0.000 3/16464 40.0 832.856 832.856 814.980 -0.000 0.000 -0.000 3/16464 50.0 629.940 629.940 621.155 -0.000 0.000 -0.000 3/16464 60.0 500.266 500.266 495.184 -0.000 0.000 -0.000 3/16464 70.0 410.483 410.483 406.836 -0.000 0.000 -0.000 3/16464 80.0 345.262 345.262 342.232 -0.000 0.000 -0.000 3/16464 90.0 296.458 296.458 293.799 -0.000 0.000 -0.000 3/16464 100.0 259.080 259.080 256.723 -0.000 0.000 -0.000 3/16464 110.0 229.836 229.836 227.748 -0.000 0.000 -0.000 3/16464 120.0 206.487 206.487 204.638 -0.000 0.000 -0.000 3/16464 130.0 187.487 187.487 185.846 -0.000 0.000 -0.000 3/16464 140.0 171.759 171.759 170.296 -0.000 0.000 -0.000 3/16464 150.0 158.537 158.537 157.226 -0.000 0.000 -0.000 3/16464 160.0 147.270 147.270 146.087 -0.000 0.000 -0.000 3/16464 170.0 137.553 137.553 136.480 -0.000 -0.000 -0.000 3/16464 180.0 129.084 129.084 128.104 -0.000 -0.000 -0.000 3/16464 190.0 121.635 121.635 120.735 -0.000 -0.000 -0.000 3/16464 200.0 115.029 115.029 114.197 -0.000 -0.000 -0.000 3/16464 210.0 109.127 109.127 108.356 -0.000 -0.000 -0.000 3/16464 220.0 103.822 103.822 103.102 -0.000 -0.000 -0.000 3/16464 230.0 99.024 99.024 98.350 -0.000 -0.000 -0.000 3/16464 240.0 94.664 94.664 94.030 -0.000 -0.000 -0.000 3/16464 250.0 90.682 90.682 90.084 -0.000 -0.000 -0.000 3/16464 260.0 87.030 87.030 86.464 -0.000 -0.000 -0.000 3/16464 270.0 83.669 83.669 83.132 -0.000 -0.000 -0.000 3/16464 280.0 80.564 80.564 80.052 -0.000 -0.000 -0.000 3/16464 290.0 77.686 77.686 77.198 -0.000 -0.000 -0.000 3/16464 300.0 75.011 75.011 74.544 -0.000 -0.000 -0.000 3/16464 310.0 72.518 72.518 72.071 -0.000 -0.000 -0.000 3/16464 320.0 70.189 70.189 69.759 -0.000 -0.000 -0.000 3/16464 330.0 68.007 68.007 67.594 -0.000 -0.000 -0.000 3/16464 340.0 65.959 65.959 65.561 -0.000 -0.000 -0.000 3/16464 350.0 64.033 64.033 63.649 -0.000 -0.000 -0.000 3/16464 360.0 62.218 62.218 61.846 -0.000 -0.000 -0.000 3/16464 370.0 60.505 60.505 60.145 -0.000 -0.000 -0.000 3/16464 380.0 58.884 58.884 58.536 -0.000 -0.000 -0.000 3/16464 390.0 57.350 57.350 57.012 -0.000 -0.000 -0.000 3/16464 400.0 55.894 55.894 55.566 -0.000 -0.000 -0.000 3/16464 410.0 54.511 54.511 54.192 -0.000 -0.000 -0.000 3/16464 420.0 53.196 53.196 52.886 -0.000 -0.000 -0.000 3/16464 430.0 51.943 51.943 51.642 -0.000 -0.000 -0.000 3/16464 440.0 50.749 50.749 50.455 -0.000 -0.000 -0.000 3/16464 450.0 49.609 49.609 49.322 -0.000 -0.000 -0.000 3/16464 460.0 48.520 48.520 48.240 -0.000 -0.000 -0.000 3/16464 470.0 47.478 47.478 47.204 -0.000 -0.000 -0.000 3/16464 480.0 46.480 46.480 46.213 -0.000 -0.000 -0.000 3/16464 490.0 45.523 45.523 45.262 -0.000 -0.000 -0.000 3/16464 500.0 44.606 44.606 44.350 -0.000 -0.000 -0.000 3/16464 510.0 43.725 43.725 43.475 -0.000 -0.000 -0.000 3/16464 520.0 42.878 42.878 42.633 -0.000 -0.000 -0.000 3/16464 530.0 42.064 42.064 41.824 -0.000 -0.000 -0.000 3/16464 540.0 41.280 41.280 41.045 -0.000 -0.000 -0.000 3/16464 550.0 40.525 40.525 40.295 -0.000 -0.000 -0.000 3/16464 560.0 39.798 39.798 39.572 -0.000 -0.000 -0.000 3/16464 570.0 39.096 39.096 38.874 -0.000 -0.000 -0.000 3/16464 580.0 38.419 38.419 38.201 -0.000 -0.000 -0.000 3/16464 590.0 37.765 37.765 37.551 -0.000 -0.000 -0.000 3/16464 600.0 37.133 37.133 36.923 -0.000 -0.000 -0.000 3/16464 610.0 36.522 36.522 36.315 -0.000 -0.000 -0.000 3/16464 620.0 35.931 35.931 35.727 -0.000 -0.000 -0.000 3/16464 630.0 35.358 35.358 35.158 -0.000 -0.000 -0.000 3/16464 640.0 34.804 34.804 34.607 -0.000 -0.000 -0.000 3/16464 650.0 34.267 34.267 34.074 -0.000 -0.000 -0.000 3/16464 660.0 33.746 33.746 33.556 -0.000 -0.000 -0.000 3/16464 670.0 33.241 33.241 33.054 -0.000 -0.000 -0.000 3/16464 680.0 32.751 32.751 32.567 -0.000 -0.000 -0.000 3/16464 690.0 32.276 32.276 32.094 -0.000 -0.000 -0.000 3/16464 700.0 31.814 31.814 31.634 -0.000 -0.000 -0.000 3/16464 710.0 31.365 31.365 31.188 -0.000 -0.000 -0.000 3/16464 720.0 30.928 30.928 30.754 -0.000 -0.000 -0.000 3/16464 730.0 30.504 30.504 30.332 -0.000 -0.000 -0.000 3/16464 740.0 30.091 30.091 29.922 -0.000 -0.000 -0.000 3/16464 750.0 29.689 29.689 29.522 -0.000 -0.000 -0.000 3/16464 760.0 29.298 29.298 29.133 -0.000 -0.000 -0.000 3/16464 770.0 28.917 28.917 28.755 -0.000 -0.000 -0.000 3/16464 780.0 28.546 28.546 28.386 -0.000 -0.000 -0.000 3/16464 790.0 28.185 28.185 28.026 -0.000 -0.000 -0.000 3/16464 800.0 27.832 27.832 27.675 -0.000 -0.000 -0.000 3/16464 810.0 27.488 27.488 27.333 -0.000 -0.000 -0.000 3/16464 820.0 27.153 27.153 27.000 -0.000 -0.000 -0.000 3/16464 830.0 26.825 26.825 26.674 -0.000 -0.000 -0.000 3/16464 840.0 26.506 26.506 26.357 -0.000 -0.000 -0.000 3/16464 850.0 26.194 26.194 26.047 -0.000 -0.000 -0.000 3/16464 860.0 25.889 25.889 25.744 -0.000 -0.000 -0.000 3/16464 870.0 25.591 25.591 25.448 -0.000 -0.000 -0.000 3/16464 880.0 25.301 25.301 25.158 -0.000 -0.000 -0.000 3/16464 890.0 25.016 25.016 24.876 -0.000 -0.000 -0.000 3/16464 900.0 24.738 24.738 24.599 -0.000 -0.000 -0.000 3/16464 910.0 24.466 24.466 24.329 -0.000 -0.000 -0.000 3/16464 920.0 24.200 24.200 24.064 -0.000 -0.000 -0.000 3/16464 930.0 23.940 23.940 23.806 -0.000 -0.000 -0.000 3/16464 940.0 23.686 23.686 23.552 -0.000 -0.000 -0.000 3/16464 950.0 23.436 23.436 23.305 -0.000 -0.000 -0.000 3/16464 960.0 23.192 23.192 23.062 -0.000 -0.000 -0.000 3/16464 970.0 22.953 22.953 22.824 -0.000 -0.000 -0.000 3/16464 980.0 22.719 22.719 22.591 -0.000 -0.000 -0.000 3/16464 990.0 22.490 22.490 22.363 -0.000 -0.000 -0.000 3/16464 1000.0 22.265 22.265 22.140 -0.000 -0.000 -0.000 3/16464 Thermal conductivity related properties were written into "kappa-m14147.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-09 04:30:06]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|