----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 01:36:34]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [2 2 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) Number of symmetry operations in supercell: 144 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.452592450260334 2.570705449749494 6.035744523333326 b -4.452592450260325 2.570705449749494 6.035744523333326 c 0.000000000000004 -5.141410899498996 6.035744523333326 Atomic positions (fractional): *1 K 0.95977788716489 0.54022211283511 0.25000000000000 39.098 2 K 0.25000000000000 0.95977788716489 0.54022211283511 39.098 3 K 0.54022211283511 0.25000000000000 0.95977788716489 39.098 4 K 0.04022211283511 0.45977788716489 0.75000000000000 39.098 5 K 0.75000000000000 0.04022211283511 0.45977788716489 39.098 6 K 0.45977788716489 0.75000000000000 0.04022211283511 39.098 *7 Se 0.63349403243070 0.63349403243070 0.63349403243070 78.960 8 Se 0.86650596756930 0.86650596756930 0.86650596756930 78.960 9 Se 0.13349403243070 0.13349403243071 0.13349403243070 78.960 10 Se 0.36650596756930 0.36650596756930 0.36650596756930 78.960 *11 Au 0.50000000000000 0.50000000000000 0.50000000000000 196.967 12 Au 0.00000000000000 0.00000000000000 0.00000000000000 196.967 -------------------------------- unit cell --------------------------------- Lattice vectors: a 8.905184900520659 0.000000000000000 0.000000000000000 b -4.452592450260330 7.712116349248490 0.000000000000000 c 0.000000000000000 0.000000000000000 18.107233569999998 Atomic positions (fractional): *1 K 0.70977788716489 0.00000000000000 0.25000000000000 39.098 > 1 2 K 0.00000000000000 0.70977788716489 0.25000000000000 39.098 > 2 3 K 0.29022211283511 0.29022211283511 0.25000000000000 39.098 > 3 4 K 0.95688877950177 0.33333333333333 0.08333333333333 39.098 > 4 5 K 0.66666666666667 0.62355544616844 0.08333333333333 39.098 > 5 6 K 0.37644455383156 0.04311122049823 0.08333333333333 39.098 > 6 7 K 0.37644455383156 0.33333333333333 0.58333333333333 39.098 > 1 8 K 0.66666666666667 0.04311122049823 0.58333333333333 39.098 > 2 9 K 0.95688877950177 0.62355544616844 0.58333333333333 39.098 > 3 10 K 0.62355544616844 0.66666666666667 0.41666666666667 39.098 > 4 11 K 0.33333333333333 0.95688877950177 0.41666666666667 39.098 > 5 12 K 0.04311122049823 0.37644455383156 0.41666666666667 39.098 > 6 13 K 0.04311122049823 0.66666666666667 0.91666666666667 39.098 > 1 14 K 0.33333333333333 0.37644455383156 0.91666666666667 39.098 > 2 15 K 0.62355544616844 0.95688877950177 0.91666666666667 39.098 > 3 16 K 0.29022211283511 0.00000000000000 0.75000000000000 39.098 > 4 17 K 0.00000000000000 0.29022211283511 0.75000000000000 39.098 > 5 18 K 0.70977788716489 0.70977788716489 0.75000000000000 39.098 > 6 *19 Se 0.33333333333333 0.66666666666667 0.30016069909737 78.960 > 7 20 Se 0.33333333333333 0.66666666666667 0.53317263423596 78.960 > 8 21 Se 0.33333333333333 0.66666666666667 0.80016069909737 78.960 > 9 22 Se 0.33333333333333 0.66666666666667 0.03317263423596 78.960 > 10 23 Se 0.00000000000000 0.00000000000000 0.63349403243070 78.960 > 7 24 Se 0.00000000000000 0.00000000000000 0.86650596756930 78.960 > 8 25 Se 0.00000000000000 0.00000000000000 0.13349403243070 78.960 > 9 26 Se 0.00000000000000 0.00000000000000 0.36650596756930 78.960 > 10 27 Se 0.66666666666667 0.33333333333333 0.96682736576404 78.960 > 7 28 Se 0.66666666666667 0.33333333333333 0.19983930090263 78.960 > 8 29 Se 0.66666666666667 0.33333333333333 0.46682736576403 78.960 > 9 30 Se 0.66666666666667 0.33333333333333 0.69983930090263 78.960 > 10 *31 Au 0.33333333333333 0.66666666666667 0.16666666666667 196.967 > 11 32 Au 0.00000000000000 0.00000000000000 0.00000000000000 196.967 > 12 33 Au 0.00000000000000 0.00000000000000 0.50000000000000 196.967 > 11 34 Au 0.66666666666667 0.33333333333333 0.33333333333333 196.967 > 12 35 Au 0.66666666666667 0.33333333333333 0.83333333333333 196.967 > 11 36 Au 0.33333333333333 0.66666666666667 0.66666666666667 196.967 > 12 -------------------------------- super cell -------------------------------- Lattice vectors: a 17.810369801041318 0.000000000000000 0.000000000000000 b -8.905184900520659 15.424232698496979 0.000000000000000 c 0.000000000000000 0.000000000000000 18.107233569999998 Atomic positions (fractional): *1 K 0.35488894358245 0.00000000000000 0.25000000000000 39.098 > 1 2 K 0.85488894358245 0.00000000000000 0.25000000000000 39.098 > 1 3 K 0.35488894358245 0.50000000000000 0.25000000000000 39.098 > 1 4 K 0.85488894358245 0.50000000000000 0.25000000000000 39.098 > 1 5 K 0.00000000000000 0.35488894358245 0.25000000000000 39.098 > 2 6 K 0.50000000000000 0.35488894358245 0.25000000000000 39.098 > 2 7 K 0.00000000000000 0.85488894358245 0.25000000000000 39.098 > 2 8 K 0.50000000000000 0.85488894358245 0.25000000000000 39.098 > 2 9 K 0.14511105641755 0.14511105641755 0.25000000000000 39.098 > 3 10 K 0.64511105641755 0.14511105641755 0.25000000000000 39.098 > 3 11 K 0.14511105641755 0.64511105641755 0.25000000000000 39.098 > 3 12 K 0.64511105641755 0.64511105641755 0.25000000000000 39.098 > 3 13 K 0.47844438975089 0.16666666666667 0.08333333333333 39.098 > 4 14 K 0.97844438975089 0.16666666666667 0.08333333333333 39.098 > 4 15 K 0.47844438975089 0.66666666666667 0.08333333333333 39.098 > 4 16 K 0.97844438975089 0.66666666666667 0.08333333333333 39.098 > 4 17 K 0.33333333333333 0.31177772308422 0.08333333333333 39.098 > 5 18 K 0.83333333333333 0.31177772308422 0.08333333333333 39.098 > 5 19 K 0.33333333333333 0.81177772308422 0.08333333333333 39.098 > 5 20 K 0.83333333333333 0.81177772308422 0.08333333333333 39.098 > 5 21 K 0.18822227691578 0.02155561024911 0.08333333333333 39.098 > 6 22 K 0.68822227691578 0.02155561024911 0.08333333333333 39.098 > 6 23 K 0.18822227691578 0.52155561024911 0.08333333333333 39.098 > 6 24 K 0.68822227691578 0.52155561024911 0.08333333333333 39.098 > 6 25 K 0.18822227691578 0.16666666666667 0.58333333333333 39.098 > 1 26 K 0.68822227691578 0.16666666666667 0.58333333333333 39.098 > 1 27 K 0.18822227691578 0.66666666666667 0.58333333333333 39.098 > 1 28 K 0.68822227691578 0.66666666666667 0.58333333333333 39.098 > 1 29 K 0.33333333333333 0.02155561024911 0.58333333333333 39.098 > 2 30 K 0.83333333333333 0.02155561024911 0.58333333333333 39.098 > 2 31 K 0.33333333333333 0.52155561024911 0.58333333333333 39.098 > 2 32 K 0.83333333333333 0.52155561024911 0.58333333333333 39.098 > 2 33 K 0.47844438975089 0.31177772308422 0.58333333333333 39.098 > 3 34 K 0.97844438975089 0.31177772308422 0.58333333333333 39.098 > 3 35 K 0.47844438975089 0.81177772308422 0.58333333333333 39.098 > 3 36 K 0.97844438975089 0.81177772308422 0.58333333333333 39.098 > 3 37 K 0.31177772308422 0.33333333333333 0.41666666666667 39.098 > 4 38 K 0.81177772308422 0.33333333333333 0.41666666666667 39.098 > 4 39 K 0.31177772308422 0.83333333333333 0.41666666666667 39.098 > 4 40 K 0.81177772308422 0.83333333333333 0.41666666666667 39.098 > 4 41 K 0.16666666666667 0.47844438975089 0.41666666666667 39.098 > 5 42 K 0.66666666666667 0.47844438975089 0.41666666666667 39.098 > 5 43 K 0.16666666666667 0.97844438975089 0.41666666666667 39.098 > 5 44 K 0.66666666666667 0.97844438975089 0.41666666666667 39.098 > 5 45 K 0.02155561024911 0.18822227691578 0.41666666666667 39.098 > 6 46 K 0.52155561024911 0.18822227691578 0.41666666666667 39.098 > 6 47 K 0.02155561024911 0.68822227691578 0.41666666666667 39.098 > 6 48 K 0.52155561024911 0.68822227691578 0.41666666666667 39.098 > 6 49 K 0.02155561024911 0.33333333333333 0.91666666666667 39.098 > 1 50 K 0.52155561024911 0.33333333333333 0.91666666666667 39.098 > 1 51 K 0.02155561024911 0.83333333333333 0.91666666666667 39.098 > 1 52 K 0.52155561024911 0.83333333333333 0.91666666666667 39.098 > 1 53 K 0.16666666666667 0.18822227691578 0.91666666666667 39.098 > 2 54 K 0.66666666666667 0.18822227691578 0.91666666666667 39.098 > 2 55 K 0.16666666666667 0.68822227691578 0.91666666666667 39.098 > 2 56 K 0.66666666666667 0.68822227691578 0.91666666666667 39.098 > 2 57 K 0.31177772308422 0.47844438975089 0.91666666666667 39.098 > 3 58 K 0.81177772308422 0.47844438975089 0.91666666666667 39.098 > 3 59 K 0.31177772308422 0.97844438975089 0.91666666666667 39.098 > 3 60 K 0.81177772308422 0.97844438975089 0.91666666666667 39.098 > 3 61 K 0.14511105641755 0.00000000000000 0.75000000000000 39.098 > 4 62 K 0.64511105641755 0.00000000000000 0.75000000000000 39.098 > 4 63 K 0.14511105641755 0.50000000000000 0.75000000000000 39.098 > 4 64 K 0.64511105641755 0.50000000000000 0.75000000000000 39.098 > 4 65 K 0.00000000000000 0.14511105641755 0.75000000000000 39.098 > 5 66 K 0.50000000000000 0.14511105641755 0.75000000000000 39.098 > 5 67 K 0.00000000000000 0.64511105641755 0.75000000000000 39.098 > 5 68 K 0.50000000000000 0.64511105641755 0.75000000000000 39.098 > 5 69 K 0.35488894358245 0.35488894358245 0.75000000000000 39.098 > 6 70 K 0.85488894358245 0.35488894358245 0.75000000000000 39.098 > 6 71 K 0.35488894358245 0.85488894358245 0.75000000000000 39.098 > 6 72 K 0.85488894358245 0.85488894358245 0.75000000000000 39.098 > 6 *73 Se 0.16666666666667 0.33333333333333 0.30016069909737 78.960 > 7 74 Se 0.66666666666667 0.33333333333333 0.30016069909737 78.960 > 7 75 Se 0.16666666666667 0.83333333333333 0.30016069909737 78.960 > 7 76 Se 0.66666666666667 0.83333333333333 0.30016069909737 78.960 > 7 77 Se 0.16666666666667 0.33333333333333 0.53317263423596 78.960 > 8 78 Se 0.66666666666667 0.33333333333333 0.53317263423596 78.960 > 8 79 Se 0.16666666666667 0.83333333333333 0.53317263423596 78.960 > 8 80 Se 0.66666666666667 0.83333333333333 0.53317263423596 78.960 > 8 81 Se 0.16666666666667 0.33333333333333 0.80016069909737 78.960 > 9 82 Se 0.66666666666667 0.33333333333333 0.80016069909737 78.960 > 9 83 Se 0.16666666666667 0.83333333333333 0.80016069909737 78.960 > 9 84 Se 0.66666666666667 0.83333333333333 0.80016069909737 78.960 > 9 85 Se 0.16666666666667 0.33333333333333 0.03317263423596 78.960 > 10 86 Se 0.66666666666667 0.33333333333333 0.03317263423596 78.960 > 10 87 Se 0.16666666666667 0.83333333333333 0.03317263423596 78.960 > 10 88 Se 0.66666666666667 0.83333333333333 0.03317263423596 78.960 > 10 89 Se 0.00000000000000 0.00000000000000 0.63349403243070 78.960 > 7 90 Se 0.50000000000000 0.00000000000000 0.63349403243070 78.960 > 7 91 Se 0.00000000000000 0.50000000000000 0.63349403243070 78.960 > 7 92 Se 0.50000000000000 0.50000000000000 0.63349403243070 78.960 > 7 93 Se 0.00000000000000 0.00000000000000 0.86650596756930 78.960 > 8 94 Se 0.50000000000000 0.00000000000000 0.86650596756930 78.960 > 8 95 Se 0.00000000000000 0.50000000000000 0.86650596756930 78.960 > 8 96 Se 0.50000000000000 0.50000000000000 0.86650596756930 78.960 > 8 97 Se 0.00000000000000 0.00000000000000 0.13349403243070 78.960 > 9 98 Se 0.50000000000000 0.00000000000000 0.13349403243070 78.960 > 9 99 Se 0.00000000000000 0.50000000000000 0.13349403243070 78.960 > 9 100 Se 0.50000000000000 0.50000000000000 0.13349403243070 78.960 > 9 101 Se 0.00000000000000 0.00000000000000 0.36650596756930 78.960 > 10 102 Se 0.50000000000000 0.00000000000000 0.36650596756930 78.960 > 10 103 Se 0.00000000000000 0.50000000000000 0.36650596756930 78.960 > 10 104 Se 0.50000000000000 0.50000000000000 0.36650596756930 78.960 > 10 105 Se 0.33333333333333 0.16666666666667 0.96682736576404 78.960 > 7 106 Se 0.83333333333333 0.16666666666667 0.96682736576404 78.960 > 7 107 Se 0.33333333333333 0.66666666666667 0.96682736576404 78.960 > 7 108 Se 0.83333333333333 0.66666666666667 0.96682736576404 78.960 > 7 109 Se 0.33333333333333 0.16666666666667 0.19983930090263 78.960 > 8 110 Se 0.83333333333333 0.16666666666667 0.19983930090263 78.960 > 8 111 Se 0.33333333333333 0.66666666666667 0.19983930090263 78.960 > 8 112 Se 0.83333333333333 0.66666666666667 0.19983930090263 78.960 > 8 113 Se 0.33333333333333 0.16666666666667 0.46682736576403 78.960 > 9 114 Se 0.83333333333333 0.16666666666667 0.46682736576403 78.960 > 9 115 Se 0.33333333333333 0.66666666666667 0.46682736576403 78.960 > 9 116 Se 0.83333333333333 0.66666666666667 0.46682736576403 78.960 > 9 117 Se 0.33333333333333 0.16666666666667 0.69983930090263 78.960 > 10 118 Se 0.83333333333333 0.16666666666667 0.69983930090263 78.960 > 10 119 Se 0.33333333333333 0.66666666666667 0.69983930090263 78.960 > 10 120 Se 0.83333333333333 0.66666666666667 0.69983930090263 78.960 > 10 *121 Au 0.16666666666667 0.33333333333333 0.16666666666667 196.967 > 11 122 Au 0.66666666666667 0.33333333333333 0.16666666666667 196.967 > 11 123 Au 0.16666666666667 0.83333333333333 0.16666666666667 196.967 > 11 124 Au 0.66666666666667 0.83333333333333 0.16666666666667 196.967 > 11 125 Au 0.00000000000000 0.00000000000000 0.00000000000000 196.967 > 12 126 Au 0.50000000000000 0.00000000000000 0.00000000000000 196.967 > 12 127 Au 0.00000000000000 0.50000000000000 0.00000000000000 196.967 > 12 128 Au 0.50000000000000 0.50000000000000 0.00000000000000 196.967 > 12 129 Au 0.00000000000000 0.00000000000000 0.50000000000000 196.967 > 11 130 Au 0.50000000000000 0.00000000000000 0.50000000000000 196.967 > 11 131 Au 0.00000000000000 0.50000000000000 0.50000000000000 196.967 > 11 132 Au 0.50000000000000 0.50000000000000 0.50000000000000 196.967 > 11 133 Au 0.33333333333333 0.16666666666667 0.33333333333333 196.967 > 12 134 Au 0.83333333333333 0.16666666666667 0.33333333333333 196.967 > 12 135 Au 0.33333333333333 0.66666666666667 0.33333333333333 196.967 > 12 136 Au 0.83333333333333 0.66666666666667 0.33333333333333 196.967 > 12 137 Au 0.33333333333333 0.16666666666667 0.83333333333333 196.967 > 11 138 Au 0.83333333333333 0.16666666666667 0.83333333333333 196.967 > 11 139 Au 0.33333333333333 0.66666666666667 0.83333333333333 196.967 > 11 140 Au 0.83333333333333 0.66666666666667 0.83333333333333 196.967 > 11 141 Au 0.16666666666667 0.33333333333333 0.66666666666667 196.967 > 12 142 Au 0.66666666666667 0.33333333333333 0.66666666666667 196.967 > 12 143 Au 0.16666666666667 0.83333333333333 0.66666666666667 196.967 > 12 144 Au 0.66666666666667 0.83333333333333 0.66666666666667 196.967 > 12 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.2908547 0.0000000 0.0000000 0.0000000 3.2908547 0.0000000 0.0000000 0.0000000 4.9458754 -------------------------- Born effective charges -------------------------- 1 K 1.1537931 -0.0000000 0.0000000 -0.0000000 1.0852954 0.0387301 0.0000000 0.1821122 1.0559230 2 K 1.1024198 -0.0296604 -0.0335413 -0.0296604 1.1366687 -0.0193651 -0.1577138 -0.0910561 1.0559230 3 K 1.1024198 0.0296604 0.0335413 0.0296604 1.1366687 -0.0193651 0.1577138 -0.0910561 1.0559230 4 K 1.1537931 -0.0000000 0.0000000 -0.0000000 1.0852954 0.0387301 0.0000000 0.1821122 1.0559230 5 K 1.1024198 -0.0296604 -0.0335413 -0.0296604 1.1366687 -0.0193651 -0.1577138 -0.0910561 1.0559230 6 K 1.1024198 0.0296604 0.0335413 0.0296604 1.1366687 -0.0193651 0.1577138 -0.0910561 1.0559230 7 Se -1.2765823 -0.0119592 0.0000000 0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 8 Se -1.2765823 0.0119592 0.0000000 -0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 9 Se -1.2765823 0.0119592 0.0000000 -0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 10 Se -1.2765823 -0.0119592 0.0000000 0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 11 Au -0.8054682 -0.0085713 0.0000000 0.0085713 -0.8054682 -0.0000000 0.0000000 -0.0000000 1.6296806 12 Au -0.8054682 0.0085713 0.0000000 -0.0085713 -0.8054682 -0.0000000 0.0000000 -0.0000000 1.6296806 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 432/432 Permutation basis: 9936/9936 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 715 Number of blocks in projector: 376 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: 338 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 272 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 105 Use standard eigh solver. Tree of FC basis block matrices: - (715, 704), data: False |-- (105, 102), data: True |-- (272, 270), data: True |-- (338, 332), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.022 Solver_block: 100 / 160 - Time: 0.780 Solver_block: 160 / 160 - Time: 0.526 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 1.335 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) Permutation basis: 432/432 Permutation basis: 9936/9936 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 715 Number of blocks in projector: 376 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: 338 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 272 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 105 Use standard eigh solver. Tree of FC basis block matrices: - (715, 704), data: False |-- (105, 102), data: True |-- (272, 270), data: True |-- (338, 332), data: True Max drift after symmetrization by symfc projector: 0.00000000 (zz) 0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 01:36:44]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 01:36:45]------------------------- 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): [2 2 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.452592450260334 2.570705449749494 6.035744523333326 b -4.452592450260325 2.570705449749494 6.035744523333326 c 0.000000000000004 -5.141410899498996 6.035744523333326 Atomic positions (fractional): 1 K 0.95977788716489 0.54022211283511 0.25000000000000 39.098 2 K 0.25000000000000 0.95977788716489 0.54022211283511 39.098 3 K 0.54022211283511 0.25000000000000 0.95977788716489 39.098 4 K 0.04022211283511 0.45977788716489 0.75000000000000 39.098 5 K 0.75000000000000 0.04022211283511 0.45977788716489 39.098 6 K 0.45977788716489 0.75000000000000 0.04022211283511 39.098 7 Se 0.63349403243070 0.63349403243070 0.63349403243070 78.960 8 Se 0.86650596756930 0.86650596756930 0.86650596756930 78.960 9 Se 0.13349403243070 0.13349403243071 0.13349403243070 78.960 10 Se 0.36650596756930 0.36650596756930 0.36650596756930 78.960 11 Au 0.50000000000000 0.50000000000000 0.50000000000000 196.967 12 Au 0.00000000000000 0.00000000000000 0.00000000000000 196.967 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.810369801041318 0.000000000000000 0.000000000000000 b -8.905184900520659 15.424232698496979 0.000000000000000 c 0.000000000000000 0.000000000000000 18.107233569999998 Atomic positions (fractional): 1 K 0.35488894358245 0.00000000000000 0.25000000000000 39.098 > 1 2 K 0.85488894358245 0.00000000000000 0.25000000000000 39.098 > 1 3 K 0.35488894358245 0.50000000000000 0.25000000000000 39.098 > 1 4 K 0.85488894358245 0.50000000000000 0.25000000000000 39.098 > 1 5 K 0.00000000000000 0.35488894358245 0.25000000000000 39.098 > 5 6 K 0.50000000000000 0.35488894358245 0.25000000000000 39.098 > 5 7 K 0.00000000000000 0.85488894358245 0.25000000000000 39.098 > 5 8 K 0.50000000000000 0.85488894358245 0.25000000000000 39.098 > 5 9 K 0.14511105641755 0.14511105641755 0.25000000000000 39.098 > 9 10 K 0.64511105641755 0.14511105641755 0.25000000000000 39.098 > 9 11 K 0.14511105641755 0.64511105641755 0.25000000000000 39.098 > 9 12 K 0.64511105641755 0.64511105641755 0.25000000000000 39.098 > 9 13 K 0.47844438975089 0.16666666666667 0.08333333333333 39.098 > 13 14 K 0.97844438975089 0.16666666666667 0.08333333333333 39.098 > 13 15 K 0.47844438975089 0.66666666666667 0.08333333333333 39.098 > 13 16 K 0.97844438975089 0.66666666666667 0.08333333333333 39.098 > 13 17 K 0.33333333333333 0.31177772308422 0.08333333333333 39.098 > 17 18 K 0.83333333333333 0.31177772308422 0.08333333333333 39.098 > 17 19 K 0.33333333333333 0.81177772308422 0.08333333333333 39.098 > 17 20 K 0.83333333333333 0.81177772308422 0.08333333333333 39.098 > 17 21 K 0.18822227691578 0.02155561024911 0.08333333333333 39.098 > 21 22 K 0.68822227691578 0.02155561024911 0.08333333333333 39.098 > 21 23 K 0.18822227691578 0.52155561024911 0.08333333333333 39.098 > 21 24 K 0.68822227691578 0.52155561024911 0.08333333333333 39.098 > 21 25 K 0.18822227691578 0.16666666666667 0.58333333333333 39.098 > 1 26 K 0.68822227691578 0.16666666666667 0.58333333333333 39.098 > 1 27 K 0.18822227691578 0.66666666666667 0.58333333333333 39.098 > 1 28 K 0.68822227691578 0.66666666666667 0.58333333333333 39.098 > 1 29 K 0.33333333333333 0.02155561024911 0.58333333333333 39.098 > 5 30 K 0.83333333333333 0.02155561024911 0.58333333333333 39.098 > 5 31 K 0.33333333333333 0.52155561024911 0.58333333333333 39.098 > 5 32 K 0.83333333333333 0.52155561024911 0.58333333333333 39.098 > 5 33 K 0.47844438975089 0.31177772308422 0.58333333333333 39.098 > 9 34 K 0.97844438975089 0.31177772308422 0.58333333333333 39.098 > 9 35 K 0.47844438975089 0.81177772308422 0.58333333333333 39.098 > 9 36 K 0.97844438975089 0.81177772308422 0.58333333333333 39.098 > 9 37 K 0.31177772308422 0.33333333333333 0.41666666666667 39.098 > 13 38 K 0.81177772308422 0.33333333333333 0.41666666666667 39.098 > 13 39 K 0.31177772308422 0.83333333333333 0.41666666666667 39.098 > 13 40 K 0.81177772308422 0.83333333333333 0.41666666666667 39.098 > 13 41 K 0.16666666666667 0.47844438975089 0.41666666666667 39.098 > 17 42 K 0.66666666666667 0.47844438975089 0.41666666666667 39.098 > 17 43 K 0.16666666666667 0.97844438975089 0.41666666666667 39.098 > 17 44 K 0.66666666666667 0.97844438975089 0.41666666666667 39.098 > 17 45 K 0.02155561024911 0.18822227691578 0.41666666666667 39.098 > 21 46 K 0.52155561024911 0.18822227691578 0.41666666666667 39.098 > 21 47 K 0.02155561024911 0.68822227691578 0.41666666666667 39.098 > 21 48 K 0.52155561024911 0.68822227691578 0.41666666666667 39.098 > 21 49 K 0.02155561024911 0.33333333333333 0.91666666666667 39.098 > 1 50 K 0.52155561024911 0.33333333333333 0.91666666666667 39.098 > 1 51 K 0.02155561024911 0.83333333333333 0.91666666666667 39.098 > 1 52 K 0.52155561024911 0.83333333333333 0.91666666666667 39.098 > 1 53 K 0.16666666666667 0.18822227691578 0.91666666666667 39.098 > 5 54 K 0.66666666666667 0.18822227691578 0.91666666666667 39.098 > 5 55 K 0.16666666666667 0.68822227691578 0.91666666666667 39.098 > 5 56 K 0.66666666666667 0.68822227691578 0.91666666666667 39.098 > 5 57 K 0.31177772308422 0.47844438975089 0.91666666666667 39.098 > 9 58 K 0.81177772308422 0.47844438975089 0.91666666666667 39.098 > 9 59 K 0.31177772308422 0.97844438975089 0.91666666666667 39.098 > 9 60 K 0.81177772308422 0.97844438975089 0.91666666666667 39.098 > 9 61 K 0.14511105641755 0.00000000000000 0.75000000000000 39.098 > 13 62 K 0.64511105641755 0.00000000000000 0.75000000000000 39.098 > 13 63 K 0.14511105641755 0.50000000000000 0.75000000000000 39.098 > 13 64 K 0.64511105641755 0.50000000000000 0.75000000000000 39.098 > 13 65 K 0.00000000000000 0.14511105641755 0.75000000000000 39.098 > 17 66 K 0.50000000000000 0.14511105641755 0.75000000000000 39.098 > 17 67 K 0.00000000000000 0.64511105641755 0.75000000000000 39.098 > 17 68 K 0.50000000000000 0.64511105641755 0.75000000000000 39.098 > 17 69 K 0.35488894358245 0.35488894358245 0.75000000000000 39.098 > 21 70 K 0.85488894358245 0.35488894358245 0.75000000000000 39.098 > 21 71 K 0.35488894358245 0.85488894358245 0.75000000000000 39.098 > 21 72 K 0.85488894358245 0.85488894358245 0.75000000000000 39.098 > 21 73 Se 0.16666666666667 0.33333333333333 0.30016069909737 78.960 > 73 74 Se 0.66666666666667 0.33333333333333 0.30016069909737 78.960 > 73 75 Se 0.16666666666667 0.83333333333333 0.30016069909737 78.960 > 73 76 Se 0.66666666666667 0.83333333333333 0.30016069909737 78.960 > 73 77 Se 0.16666666666667 0.33333333333333 0.53317263423596 78.960 > 77 78 Se 0.66666666666667 0.33333333333333 0.53317263423596 78.960 > 77 79 Se 0.16666666666667 0.83333333333333 0.53317263423596 78.960 > 77 80 Se 0.66666666666667 0.83333333333333 0.53317263423596 78.960 > 77 81 Se 0.16666666666667 0.33333333333333 0.80016069909737 78.960 > 81 82 Se 0.66666666666667 0.33333333333333 0.80016069909737 78.960 > 81 83 Se 0.16666666666667 0.83333333333333 0.80016069909737 78.960 > 81 84 Se 0.66666666666667 0.83333333333333 0.80016069909737 78.960 > 81 85 Se 0.16666666666667 0.33333333333333 0.03317263423596 78.960 > 85 86 Se 0.66666666666667 0.33333333333333 0.03317263423596 78.960 > 85 87 Se 0.16666666666667 0.83333333333333 0.03317263423596 78.960 > 85 88 Se 0.66666666666667 0.83333333333333 0.03317263423596 78.960 > 85 89 Se 0.00000000000000 0.00000000000000 0.63349403243070 78.960 > 73 90 Se 0.50000000000000 0.00000000000000 0.63349403243070 78.960 > 73 91 Se 0.00000000000000 0.50000000000000 0.63349403243070 78.960 > 73 92 Se 0.50000000000000 0.50000000000000 0.63349403243070 78.960 > 73 93 Se 0.00000000000000 0.00000000000000 0.86650596756930 78.960 > 77 94 Se 0.50000000000000 0.00000000000000 0.86650596756930 78.960 > 77 95 Se 0.00000000000000 0.50000000000000 0.86650596756930 78.960 > 77 96 Se 0.50000000000000 0.50000000000000 0.86650596756930 78.960 > 77 97 Se 0.00000000000000 0.00000000000000 0.13349403243070 78.960 > 81 98 Se 0.50000000000000 0.00000000000000 0.13349403243070 78.960 > 81 99 Se 0.00000000000000 0.50000000000000 0.13349403243070 78.960 > 81 100 Se 0.50000000000000 0.50000000000000 0.13349403243070 78.960 > 81 101 Se 0.00000000000000 0.00000000000000 0.36650596756930 78.960 > 85 102 Se 0.50000000000000 0.00000000000000 0.36650596756930 78.960 > 85 103 Se 0.00000000000000 0.50000000000000 0.36650596756930 78.960 > 85 104 Se 0.50000000000000 0.50000000000000 0.36650596756930 78.960 > 85 105 Se 0.33333333333333 0.16666666666667 0.96682736576404 78.960 > 73 106 Se 0.83333333333333 0.16666666666667 0.96682736576404 78.960 > 73 107 Se 0.33333333333333 0.66666666666667 0.96682736576404 78.960 > 73 108 Se 0.83333333333333 0.66666666666667 0.96682736576404 78.960 > 73 109 Se 0.33333333333333 0.16666666666667 0.19983930090263 78.960 > 77 110 Se 0.83333333333333 0.16666666666667 0.19983930090263 78.960 > 77 111 Se 0.33333333333333 0.66666666666667 0.19983930090263 78.960 > 77 112 Se 0.83333333333333 0.66666666666667 0.19983930090263 78.960 > 77 113 Se 0.33333333333333 0.16666666666667 0.46682736576403 78.960 > 81 114 Se 0.83333333333333 0.16666666666667 0.46682736576403 78.960 > 81 115 Se 0.33333333333333 0.66666666666667 0.46682736576403 78.960 > 81 116 Se 0.83333333333333 0.66666666666667 0.46682736576403 78.960 > 81 117 Se 0.33333333333333 0.16666666666667 0.69983930090263 78.960 > 85 118 Se 0.83333333333333 0.16666666666667 0.69983930090263 78.960 > 85 119 Se 0.33333333333333 0.66666666666667 0.69983930090263 78.960 > 85 120 Se 0.83333333333333 0.66666666666667 0.69983930090263 78.960 > 85 121 Au 0.16666666666667 0.33333333333333 0.16666666666667 196.967 > 121 122 Au 0.66666666666667 0.33333333333333 0.16666666666667 196.967 > 121 123 Au 0.16666666666667 0.83333333333333 0.16666666666667 196.967 > 121 124 Au 0.66666666666667 0.83333333333333 0.16666666666667 196.967 > 121 125 Au 0.00000000000000 0.00000000000000 0.00000000000000 196.967 > 125 126 Au 0.50000000000000 0.00000000000000 0.00000000000000 196.967 > 125 127 Au 0.00000000000000 0.50000000000000 0.00000000000000 196.967 > 125 128 Au 0.50000000000000 0.50000000000000 0.00000000000000 196.967 > 125 129 Au 0.00000000000000 0.00000000000000 0.50000000000000 196.967 > 121 130 Au 0.50000000000000 0.00000000000000 0.50000000000000 196.967 > 121 131 Au 0.00000000000000 0.50000000000000 0.50000000000000 196.967 > 121 132 Au 0.50000000000000 0.50000000000000 0.50000000000000 196.967 > 121 133 Au 0.33333333333333 0.16666666666667 0.33333333333333 196.967 > 125 134 Au 0.83333333333333 0.16666666666667 0.33333333333333 196.967 > 125 135 Au 0.33333333333333 0.66666666666667 0.33333333333333 196.967 > 125 136 Au 0.83333333333333 0.66666666666667 0.33333333333333 196.967 > 125 137 Au 0.33333333333333 0.16666666666667 0.83333333333333 196.967 > 121 138 Au 0.83333333333333 0.16666666666667 0.83333333333333 196.967 > 121 139 Au 0.33333333333333 0.66666666666667 0.83333333333333 196.967 > 121 140 Au 0.83333333333333 0.66666666666667 0.83333333333333 196.967 > 121 141 Au 0.16666666666667 0.33333333333333 0.66666666666667 196.967 > 125 142 Au 0.66666666666667 0.33333333333333 0.66666666666667 196.967 > 125 143 Au 0.16666666666667 0.83333333333333 0.66666666666667 196.967 > 125 144 Au 0.66666666666667 0.83333333333333 0.66666666666667 196.967 > 125 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.2908547 0.0000000 0.0000000 0.0000000 3.2908547 0.0000000 0.0000000 0.0000000 4.9458754 -------------------------- Born effective charges -------------------------- 1 K 1.1537931 -0.0000000 0.0000000 -0.0000000 1.0852954 0.0387301 0.0000000 0.1821122 1.0559230 2 K 1.1024198 -0.0296604 -0.0335413 -0.0296604 1.1366687 -0.0193651 -0.1577138 -0.0910561 1.0559230 3 K 1.1024198 0.0296604 0.0335413 0.0296604 1.1366687 -0.0193651 0.1577138 -0.0910561 1.0559230 4 K 1.1537931 -0.0000000 0.0000000 -0.0000000 1.0852954 0.0387301 0.0000000 0.1821122 1.0559230 5 K 1.1024198 -0.0296604 -0.0335413 -0.0296604 1.1366687 -0.0193651 -0.1577138 -0.0910561 1.0559230 6 K 1.1024198 0.0296604 0.0335413 0.0296604 1.1366687 -0.0193651 0.1577138 -0.0910561 1.0559230 7 Se -1.2765823 -0.0119592 0.0000000 0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 8 Se -1.2765823 0.0119592 0.0000000 -0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 9 Se -1.2765823 0.0119592 0.0000000 -0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 10 Se -1.2765823 -0.0119592 0.0000000 0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 11 Au -0.8054682 -0.0085713 0.0000000 0.0085713 -0.8054682 -0.0000000 0.0000000 -0.0000000 1.6296806 12 Au -0.8054682 0.0085713 0.0000000 -0.0085713 -0.8054682 -0.0000000 0.0000000 -0.0000000 1.6296806 ---------------------------------------------------------------------------- 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.0050 0.0087 0.0000] [ 0.0050 -0.0087 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 73, 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[ 121, 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.00000000 (zxz) -0.00000000 (zxz) -0.00000000 (zzx) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000001 (yx) 0.00000001 (yx) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 01:36:57]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 01:36:58]------------------------- 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): [2 2 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.452592450260334 2.570705449749494 6.035744523333326 b -4.452592450260325 2.570705449749494 6.035744523333326 c 0.000000000000004 -5.141410899498996 6.035744523333326 Atomic positions (fractional): 1 K 0.95977788716489 0.54022211283511 0.25000000000000 39.098 2 K 0.25000000000000 0.95977788716489 0.54022211283511 39.098 3 K 0.54022211283511 0.25000000000000 0.95977788716489 39.098 4 K 0.04022211283511 0.45977788716489 0.75000000000000 39.098 5 K 0.75000000000000 0.04022211283511 0.45977788716489 39.098 6 K 0.45977788716489 0.75000000000000 0.04022211283511 39.098 7 Se 0.63349403243070 0.63349403243070 0.63349403243070 78.960 8 Se 0.86650596756930 0.86650596756930 0.86650596756930 78.960 9 Se 0.13349403243070 0.13349403243071 0.13349403243070 78.960 10 Se 0.36650596756930 0.36650596756930 0.36650596756930 78.960 11 Au 0.50000000000000 0.50000000000000 0.50000000000000 196.967 12 Au 0.00000000000000 0.00000000000000 0.00000000000000 196.967 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.810369801041318 0.000000000000000 0.000000000000000 b -8.905184900520659 15.424232698496979 0.000000000000000 c 0.000000000000000 0.000000000000000 18.107233569999998 Atomic positions (fractional): 1 K 0.35488894358245 0.00000000000000 0.25000000000000 39.098 > 1 2 K 0.85488894358245 0.00000000000000 0.25000000000000 39.098 > 1 3 K 0.35488894358245 0.50000000000000 0.25000000000000 39.098 > 1 4 K 0.85488894358245 0.50000000000000 0.25000000000000 39.098 > 1 5 K 0.00000000000000 0.35488894358245 0.25000000000000 39.098 > 5 6 K 0.50000000000000 0.35488894358245 0.25000000000000 39.098 > 5 7 K 0.00000000000000 0.85488894358245 0.25000000000000 39.098 > 5 8 K 0.50000000000000 0.85488894358245 0.25000000000000 39.098 > 5 9 K 0.14511105641755 0.14511105641755 0.25000000000000 39.098 > 9 10 K 0.64511105641755 0.14511105641755 0.25000000000000 39.098 > 9 11 K 0.14511105641755 0.64511105641755 0.25000000000000 39.098 > 9 12 K 0.64511105641755 0.64511105641755 0.25000000000000 39.098 > 9 13 K 0.47844438975089 0.16666666666667 0.08333333333333 39.098 > 13 14 K 0.97844438975089 0.16666666666667 0.08333333333333 39.098 > 13 15 K 0.47844438975089 0.66666666666667 0.08333333333333 39.098 > 13 16 K 0.97844438975089 0.66666666666667 0.08333333333333 39.098 > 13 17 K 0.33333333333333 0.31177772308422 0.08333333333333 39.098 > 17 18 K 0.83333333333333 0.31177772308422 0.08333333333333 39.098 > 17 19 K 0.33333333333333 0.81177772308422 0.08333333333333 39.098 > 17 20 K 0.83333333333333 0.81177772308422 0.08333333333333 39.098 > 17 21 K 0.18822227691578 0.02155561024911 0.08333333333333 39.098 > 21 22 K 0.68822227691578 0.02155561024911 0.08333333333333 39.098 > 21 23 K 0.18822227691578 0.52155561024911 0.08333333333333 39.098 > 21 24 K 0.68822227691578 0.52155561024911 0.08333333333333 39.098 > 21 25 K 0.18822227691578 0.16666666666667 0.58333333333333 39.098 > 1 26 K 0.68822227691578 0.16666666666667 0.58333333333333 39.098 > 1 27 K 0.18822227691578 0.66666666666667 0.58333333333333 39.098 > 1 28 K 0.68822227691578 0.66666666666667 0.58333333333333 39.098 > 1 29 K 0.33333333333333 0.02155561024911 0.58333333333333 39.098 > 5 30 K 0.83333333333333 0.02155561024911 0.58333333333333 39.098 > 5 31 K 0.33333333333333 0.52155561024911 0.58333333333333 39.098 > 5 32 K 0.83333333333333 0.52155561024911 0.58333333333333 39.098 > 5 33 K 0.47844438975089 0.31177772308422 0.58333333333333 39.098 > 9 34 K 0.97844438975089 0.31177772308422 0.58333333333333 39.098 > 9 35 K 0.47844438975089 0.81177772308422 0.58333333333333 39.098 > 9 36 K 0.97844438975089 0.81177772308422 0.58333333333333 39.098 > 9 37 K 0.31177772308422 0.33333333333333 0.41666666666667 39.098 > 13 38 K 0.81177772308422 0.33333333333333 0.41666666666667 39.098 > 13 39 K 0.31177772308422 0.83333333333333 0.41666666666667 39.098 > 13 40 K 0.81177772308422 0.83333333333333 0.41666666666667 39.098 > 13 41 K 0.16666666666667 0.47844438975089 0.41666666666667 39.098 > 17 42 K 0.66666666666667 0.47844438975089 0.41666666666667 39.098 > 17 43 K 0.16666666666667 0.97844438975089 0.41666666666667 39.098 > 17 44 K 0.66666666666667 0.97844438975089 0.41666666666667 39.098 > 17 45 K 0.02155561024911 0.18822227691578 0.41666666666667 39.098 > 21 46 K 0.52155561024911 0.18822227691578 0.41666666666667 39.098 > 21 47 K 0.02155561024911 0.68822227691578 0.41666666666667 39.098 > 21 48 K 0.52155561024911 0.68822227691578 0.41666666666667 39.098 > 21 49 K 0.02155561024911 0.33333333333333 0.91666666666667 39.098 > 1 50 K 0.52155561024911 0.33333333333333 0.91666666666667 39.098 > 1 51 K 0.02155561024911 0.83333333333333 0.91666666666667 39.098 > 1 52 K 0.52155561024911 0.83333333333333 0.91666666666667 39.098 > 1 53 K 0.16666666666667 0.18822227691578 0.91666666666667 39.098 > 5 54 K 0.66666666666667 0.18822227691578 0.91666666666667 39.098 > 5 55 K 0.16666666666667 0.68822227691578 0.91666666666667 39.098 > 5 56 K 0.66666666666667 0.68822227691578 0.91666666666667 39.098 > 5 57 K 0.31177772308422 0.47844438975089 0.91666666666667 39.098 > 9 58 K 0.81177772308422 0.47844438975089 0.91666666666667 39.098 > 9 59 K 0.31177772308422 0.97844438975089 0.91666666666667 39.098 > 9 60 K 0.81177772308422 0.97844438975089 0.91666666666667 39.098 > 9 61 K 0.14511105641755 0.00000000000000 0.75000000000000 39.098 > 13 62 K 0.64511105641755 0.00000000000000 0.75000000000000 39.098 > 13 63 K 0.14511105641755 0.50000000000000 0.75000000000000 39.098 > 13 64 K 0.64511105641755 0.50000000000000 0.75000000000000 39.098 > 13 65 K 0.00000000000000 0.14511105641755 0.75000000000000 39.098 > 17 66 K 0.50000000000000 0.14511105641755 0.75000000000000 39.098 > 17 67 K 0.00000000000000 0.64511105641755 0.75000000000000 39.098 > 17 68 K 0.50000000000000 0.64511105641755 0.75000000000000 39.098 > 17 69 K 0.35488894358245 0.35488894358245 0.75000000000000 39.098 > 21 70 K 0.85488894358245 0.35488894358245 0.75000000000000 39.098 > 21 71 K 0.35488894358245 0.85488894358245 0.75000000000000 39.098 > 21 72 K 0.85488894358245 0.85488894358245 0.75000000000000 39.098 > 21 73 Se 0.16666666666667 0.33333333333333 0.30016069909737 78.960 > 73 74 Se 0.66666666666667 0.33333333333333 0.30016069909737 78.960 > 73 75 Se 0.16666666666667 0.83333333333333 0.30016069909737 78.960 > 73 76 Se 0.66666666666667 0.83333333333333 0.30016069909737 78.960 > 73 77 Se 0.16666666666667 0.33333333333333 0.53317263423596 78.960 > 77 78 Se 0.66666666666667 0.33333333333333 0.53317263423596 78.960 > 77 79 Se 0.16666666666667 0.83333333333333 0.53317263423596 78.960 > 77 80 Se 0.66666666666667 0.83333333333333 0.53317263423596 78.960 > 77 81 Se 0.16666666666667 0.33333333333333 0.80016069909737 78.960 > 81 82 Se 0.66666666666667 0.33333333333333 0.80016069909737 78.960 > 81 83 Se 0.16666666666667 0.83333333333333 0.80016069909737 78.960 > 81 84 Se 0.66666666666667 0.83333333333333 0.80016069909737 78.960 > 81 85 Se 0.16666666666667 0.33333333333333 0.03317263423596 78.960 > 85 86 Se 0.66666666666667 0.33333333333333 0.03317263423596 78.960 > 85 87 Se 0.16666666666667 0.83333333333333 0.03317263423596 78.960 > 85 88 Se 0.66666666666667 0.83333333333333 0.03317263423596 78.960 > 85 89 Se 0.00000000000000 0.00000000000000 0.63349403243070 78.960 > 73 90 Se 0.50000000000000 0.00000000000000 0.63349403243070 78.960 > 73 91 Se 0.00000000000000 0.50000000000000 0.63349403243070 78.960 > 73 92 Se 0.50000000000000 0.50000000000000 0.63349403243070 78.960 > 73 93 Se 0.00000000000000 0.00000000000000 0.86650596756930 78.960 > 77 94 Se 0.50000000000000 0.00000000000000 0.86650596756930 78.960 > 77 95 Se 0.00000000000000 0.50000000000000 0.86650596756930 78.960 > 77 96 Se 0.50000000000000 0.50000000000000 0.86650596756930 78.960 > 77 97 Se 0.00000000000000 0.00000000000000 0.13349403243070 78.960 > 81 98 Se 0.50000000000000 0.00000000000000 0.13349403243070 78.960 > 81 99 Se 0.00000000000000 0.50000000000000 0.13349403243070 78.960 > 81 100 Se 0.50000000000000 0.50000000000000 0.13349403243070 78.960 > 81 101 Se 0.00000000000000 0.00000000000000 0.36650596756930 78.960 > 85 102 Se 0.50000000000000 0.00000000000000 0.36650596756930 78.960 > 85 103 Se 0.00000000000000 0.50000000000000 0.36650596756930 78.960 > 85 104 Se 0.50000000000000 0.50000000000000 0.36650596756930 78.960 > 85 105 Se 0.33333333333333 0.16666666666667 0.96682736576404 78.960 > 73 106 Se 0.83333333333333 0.16666666666667 0.96682736576404 78.960 > 73 107 Se 0.33333333333333 0.66666666666667 0.96682736576404 78.960 > 73 108 Se 0.83333333333333 0.66666666666667 0.96682736576404 78.960 > 73 109 Se 0.33333333333333 0.16666666666667 0.19983930090263 78.960 > 77 110 Se 0.83333333333333 0.16666666666667 0.19983930090263 78.960 > 77 111 Se 0.33333333333333 0.66666666666667 0.19983930090263 78.960 > 77 112 Se 0.83333333333333 0.66666666666667 0.19983930090263 78.960 > 77 113 Se 0.33333333333333 0.16666666666667 0.46682736576403 78.960 > 81 114 Se 0.83333333333333 0.16666666666667 0.46682736576403 78.960 > 81 115 Se 0.33333333333333 0.66666666666667 0.46682736576403 78.960 > 81 116 Se 0.83333333333333 0.66666666666667 0.46682736576403 78.960 > 81 117 Se 0.33333333333333 0.16666666666667 0.69983930090263 78.960 > 85 118 Se 0.83333333333333 0.16666666666667 0.69983930090263 78.960 > 85 119 Se 0.33333333333333 0.66666666666667 0.69983930090263 78.960 > 85 120 Se 0.83333333333333 0.66666666666667 0.69983930090263 78.960 > 85 121 Au 0.16666666666667 0.33333333333333 0.16666666666667 196.967 > 121 122 Au 0.66666666666667 0.33333333333333 0.16666666666667 196.967 > 121 123 Au 0.16666666666667 0.83333333333333 0.16666666666667 196.967 > 121 124 Au 0.66666666666667 0.83333333333333 0.16666666666667 196.967 > 121 125 Au 0.00000000000000 0.00000000000000 0.00000000000000 196.967 > 125 126 Au 0.50000000000000 0.00000000000000 0.00000000000000 196.967 > 125 127 Au 0.00000000000000 0.50000000000000 0.00000000000000 196.967 > 125 128 Au 0.50000000000000 0.50000000000000 0.00000000000000 196.967 > 125 129 Au 0.00000000000000 0.00000000000000 0.50000000000000 196.967 > 121 130 Au 0.50000000000000 0.00000000000000 0.50000000000000 196.967 > 121 131 Au 0.00000000000000 0.50000000000000 0.50000000000000 196.967 > 121 132 Au 0.50000000000000 0.50000000000000 0.50000000000000 196.967 > 121 133 Au 0.33333333333333 0.16666666666667 0.33333333333333 196.967 > 125 134 Au 0.83333333333333 0.16666666666667 0.33333333333333 196.967 > 125 135 Au 0.33333333333333 0.66666666666667 0.33333333333333 196.967 > 125 136 Au 0.83333333333333 0.66666666666667 0.33333333333333 196.967 > 125 137 Au 0.33333333333333 0.16666666666667 0.83333333333333 196.967 > 121 138 Au 0.83333333333333 0.16666666666667 0.83333333333333 196.967 > 121 139 Au 0.33333333333333 0.66666666666667 0.83333333333333 196.967 > 121 140 Au 0.83333333333333 0.66666666666667 0.83333333333333 196.967 > 121 141 Au 0.16666666666667 0.33333333333333 0.66666666666667 196.967 > 125 142 Au 0.66666666666667 0.33333333333333 0.66666666666667 196.967 > 125 143 Au 0.16666666666667 0.83333333333333 0.66666666666667 196.967 > 125 144 Au 0.66666666666667 0.83333333333333 0.66666666666667 196.967 > 125 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 3.2908547 0.0000000 0.0000000 0.0000000 3.2908547 0.0000000 0.0000000 0.0000000 4.9458754 -------------------------- Born effective charges -------------------------- 1 K 1.1537931 -0.0000000 0.0000000 -0.0000000 1.0852954 0.0387301 0.0000000 0.1821122 1.0559230 2 K 1.1024198 -0.0296604 -0.0335413 -0.0296604 1.1366687 -0.0193651 -0.1577138 -0.0910561 1.0559230 3 K 1.1024198 0.0296604 0.0335413 0.0296604 1.1366687 -0.0193651 0.1577138 -0.0910561 1.0559230 4 K 1.1537931 -0.0000000 0.0000000 -0.0000000 1.0852954 0.0387301 0.0000000 0.1821122 1.0559230 5 K 1.1024198 -0.0296604 -0.0335413 -0.0296604 1.1366687 -0.0193651 -0.1577138 -0.0910561 1.0559230 6 K 1.1024198 0.0296604 0.0335413 0.0296604 1.1366687 -0.0193651 0.1577138 -0.0910561 1.0559230 7 Se -1.2765823 -0.0119592 0.0000000 0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 8 Se -1.2765823 0.0119592 0.0000000 -0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 9 Se -1.2765823 0.0119592 0.0000000 -0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 10 Se -1.2765823 -0.0119592 0.0000000 0.0119592 -1.2765823 -0.0000000 0.0000000 -0.0000000 -2.3987248 11 Au -0.8054682 -0.0085713 0.0000000 0.0085713 -0.8054682 -0.0000000 0.0000000 -0.0000000 1.6296806 12 Au -0.8054682 0.0085713 0.0000000 -0.0085713 -0.8054682 -0.0000000 0.0000000 -0.0000000 1.6296806 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (zxz) -0.00000000 (zxz) -0.00000000 (zzx) 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 ... [ 7 7 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.56, Number of G-points: 289, Lambda: 0.11 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/44) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 44 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.411 ( 0.000 0.000 0.000) 0.000 1.411 ( 0.000 0.000 0.000) 0.000 1.462 ( 0.000 0.000 0.000) 0.000 1.987 ( 0.000 0.000 0.000) 0.000 1.987 ( 0.000 0.000 0.000) 0.000 2.159 ( 0.000 0.000 0.000) 0.000 2.159 ( 0.000 0.000 0.000) 0.000 2.225 ( 0.000 0.000 0.000) 0.000 2.225 ( 0.000 0.000 0.000) 0.000 2.408 ( 0.000 0.000 0.000) 0.000 2.542 ( 0.000 0.000 0.000) 0.000 2.542 ( 0.000 0.000 0.000) 0.000 2.660 ( 0.000 0.000 0.000) 0.000 2.660 ( 0.000 0.000 0.000) 0.000 2.734 ( 0.000 0.000 0.000) 0.000 2.734 ( 0.000 0.000 0.000) 0.000 2.976 ( 0.000 0.000 0.000) 0.000 3.023 ( 0.000 0.000 0.000) 0.000 3.240 ( 0.000 0.000 0.000) 0.000 3.327 ( 0.000 0.000 0.000) 0.000 3.327 ( 0.000 0.000 0.000) 0.000 3.630 ( 0.000 0.000 0.000) 0.000 3.630 ( 0.000 0.000 0.000) 0.000 3.717 ( 0.000 0.000 0.000) 0.000 3.865 ( 0.000 0.000 0.000) 0.000 3.865 ( 0.000 0.000 0.000) 0.000 3.937 ( 0.000 0.000 0.000) 0.000 4.150 ( 0.000 0.000 0.000) 0.000 4.150 ( 0.000 0.000 0.000) 0.000 6.219 ( 0.000 0.000 0.000) 0.000 6.308 ( 0.000 0.000 0.000) 0.000 7.133 ( 0.000 0.000 0.000) 0.000 7.269 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/44) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.272 ( 10.752 6.208 5.293) 13.496 0.293 ( 10.880 6.281 7.184) 14.472 0.543 ( 19.938 11.511 12.133) 26.024 1.363 ( -4.879 -2.817 2.011) 5.983 1.369 ( -3.640 -2.102 -0.678) 4.258 1.424 ( -2.281 -1.317 -3.386) 4.290 1.886 ( -6.997 -4.040 -4.522) 9.259 1.987 ( -0.200 -0.116 -0.848) 0.879 2.143 ( -2.163 -1.249 -1.169) 2.758 2.152 ( -1.810 -1.045 0.345) 2.118 2.208 ( -0.065 -0.037 1.004) 1.007 2.264 ( 1.761 1.017 4.658) 5.083 2.454 ( -2.956 -1.707 -2.450) 4.202 2.517 ( -3.418 -1.973 3.996) 5.617 2.550 ( -4.497 -2.596 10.817) 11.998 2.592 ( -1.308 -0.755 -7.337) 7.491 2.690 ( 1.304 0.753 -8.220) 8.356 2.777 ( -0.631 -0.364 6.306) 6.348 2.798 ( 0.379 0.219 8.993) 9.003 3.019 ( -0.270 -0.156 -1.471) 1.503 3.048 ( 5.958 3.440 -0.362) 6.889 3.270 ( 0.676 0.390 6.352) 6.400 3.289 ( -2.737 -1.580 0.411) 3.187 3.352 ( 1.288 0.744 2.702) 3.084 3.623 ( -0.245 -0.142 -1.505) 1.531 3.709 ( 0.375 0.217 1.264) 1.336 3.757 ( 3.236 1.869 -1.188) 3.921 3.843 ( -1.122 -0.648 -1.978) 2.364 3.992 ( 3.840 2.217 1.888) 4.819 4.132 ( -1.556 -0.899 -0.403) 1.842 4.139 ( -0.153 -0.088 -1.239) 1.252 4.603 ( 11.588 6.690 -36.022) 38.427 6.226 ( 0.216 0.125 1.166) 1.192 6.303 ( -0.105 -0.060 -0.846) 0.855 7.247 ( -7.631 -4.406 21.002) 22.775 7.266 ( 0.419 0.242 -1.691) 1.759 ======================= Grid point 2 (3/44) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.541 ( 10.677 6.164 4.706) 13.196 0.576 ( 10.105 5.834 7.533) 13.888 1.016 ( 14.623 8.442 10.072) 19.661 1.250 ( -6.652 -3.840 0.617) 7.706 1.252 ( -5.922 -3.419 -0.898) 6.896 1.327 ( -3.672 -2.120 -2.505) 4.924 1.670 ( -8.103 -4.678 -5.305) 10.756 1.923 ( -3.474 -2.006 -6.602) 7.725 2.096 ( -0.056 -0.032 -1.079) 1.081 2.148 ( 0.240 0.139 2.337) 2.354 2.251 ( 0.817 0.472 6.003) 6.076 2.317 ( -0.602 -0.347 -1.458) 1.615 2.380 ( -0.006 -0.004 0.584) 0.584 2.474 ( -4.767 -2.752 7.730) 9.489 2.511 ( -0.952 -0.550 -6.532) 6.624 2.560 ( -0.513 -0.296 2.941) 3.001 2.624 ( 0.347 0.200 -10.847) 10.855 2.779 ( -3.003 -1.734 3.310) 4.794 2.872 ( 1.020 0.589 5.648) 5.770 2.983 ( -1.064 -0.614 -2.377) 2.676 3.158 ( 4.470 2.581 -2.381) 5.685 3.249 ( -2.068 -1.194 2.718) 3.618 3.372 ( 1.702 0.983 9.526) 9.727 3.429 ( 3.722 2.149 5.077) 6.651 3.589 ( -0.834 -0.482 -4.292) 4.399 3.734 ( 0.389 0.225 3.205) 3.237 3.775 ( -1.406 -0.812 -3.144) 3.539 3.822 ( 2.368 1.367 -3.476) 4.423 4.044 ( -0.435 -0.251 -2.116) 2.175 4.080 ( -2.633 -1.520 -0.695) 3.118 4.193 ( 5.045 2.913 4.273) 7.225 4.552 ( 3.281 1.894 -16.942) 17.360 6.244 ( 0.294 0.170 1.890) 1.920 6.292 ( -0.026 -0.015 -1.359) 1.360 7.249 ( -3.698 -2.135 10.236) 11.091 7.260 ( 0.869 0.502 -3.232) 3.384 ======================= Grid point 3 (4/44) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.799 ( 10.207 5.893 3.703) 12.354 0.841 ( 8.978 5.183 7.560) 12.831 1.033 ( -9.110 -5.259 -2.098) 10.726 1.075 ( -7.586 -4.380 -6.032) 10.635 1.283 ( -0.482 -0.278 4.043) 4.081 1.320 ( 2.053 1.185 0.177) 2.377 1.564 ( 3.247 1.875 3.645) 5.229 1.756 ( -5.998 -3.463 -9.124) 11.455 2.115 ( 1.779 1.027 0.210) 2.065 2.158 ( -1.810 -1.045 1.352) 2.489 2.284 ( 1.494 0.863 -2.267) 2.848 2.304 ( -0.157 -0.090 2.379) 2.386 2.424 ( 2.060 1.190 -0.581) 2.449 2.444 ( -4.532 -2.617 11.375) 12.521 2.465 ( 0.064 0.037 -1.637) 1.638 2.518 ( -0.276 -0.160 -10.391) 10.396 2.613 ( 2.610 1.507 0.883) 3.140 2.703 ( -4.635 -2.676 -0.254) 5.358 2.927 ( 0.599 0.346 3.641) 3.706 2.955 ( -0.053 -0.031 -2.302) 2.303 3.210 ( 2.683 1.549 -2.878) 4.229 3.233 ( -1.964 -1.134 3.410) 4.095 3.456 ( 0.445 0.257 7.581) 7.598 3.525 ( -0.832 -0.480 -6.682) 6.751 3.579 ( 4.741 2.737 5.852) 8.013 3.699 ( -2.993 -1.728 -3.795) 5.133 3.778 ( 0.483 0.279 5.515) 5.543 3.826 ( 0.996 0.575 -6.141) 6.248 4.007 ( -0.357 -0.206 -1.625) 1.677 4.022 ( -2.066 -1.193 -0.122) 2.389 4.348 ( 4.221 2.437 7.248) 8.735 4.472 ( -0.864 -0.499 -10.325) 10.373 6.265 ( 0.151 0.087 1.970) 1.978 6.281 ( 0.012 0.007 -1.687) 1.688 7.250 ( -2.315 -1.336 6.539) 7.064 7.254 ( 1.466 0.847 -4.662) 4.960 ======================= Grid point 8 (5/44) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.343 ( 0.000 10.149 9.847) 14.141 0.360 ( 0.000 11.237 9.168) 14.503 0.661 ( 0.000 17.035 19.079) 25.577 1.334 ( -0.000 -7.484 -0.935) 7.542 1.389 ( -0.000 -2.212 -5.588) 6.010 1.412 ( 0.000 -2.252 2.912) 3.681 1.862 ( -0.000 -3.635 -7.976) 8.766 1.932 ( -0.000 -3.029 -2.331) 3.822 2.127 ( -0.000 -6.215 -3.697) 7.231 2.191 ( 0.000 -1.803 2.003) 2.695 2.209 ( 0.000 0.299 3.646) 3.658 2.284 ( 0.000 1.077 3.072) 3.255 2.404 ( -0.000 -4.902 -3.299) 5.909 2.535 ( -0.000 1.995 -2.332) 3.069 2.586 ( -0.000 -1.325 -0.240) 1.346 2.614 ( -0.000 -1.992 -6.395) 6.699 2.688 ( 0.000 -10.951 8.397) 13.800 2.781 ( 0.000 -0.818 7.035) 7.082 2.875 ( 0.000 2.169 6.043) 6.420 3.001 ( -0.000 -0.479 -2.301) 2.351 3.050 ( -0.000 6.466 0.503) 6.485 3.311 ( 0.000 -1.142 1.009) 1.524 3.335 ( 0.000 -0.078 8.384) 8.385 3.340 ( 0.000 -1.366 4.152) 4.371 3.658 ( 0.000 1.313 1.516) 2.006 3.667 ( -0.000 -1.118 -1.923) 2.224 3.754 ( -0.000 4.195 -1.600) 4.490 3.857 ( -0.000 -0.020 -0.820) 0.820 4.014 ( 0.000 4.169 3.513) 5.452 4.099 ( -0.000 -2.294 -3.283) 4.005 4.156 ( 0.000 0.749 0.099) 0.756 4.317 ( -0.000 21.688 -29.892) 36.931 6.239 ( 0.000 0.168 2.189) 2.196 6.299 ( -0.000 -0.693 -0.186) 0.717 7.247 ( -0.000 0.418 -3.128) 3.156 7.384 ( 0.000 -11.306 12.855) 17.120 ======================= Grid point 9 (6/44) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.533 ( 3.691 8.317 8.681) 12.576 0.588 ( 5.739 9.405 8.908) 14.169 1.028 ( 8.450 12.744 14.553) 21.109 1.222 ( -2.322 -8.619 -2.441) 9.254 1.295 ( -2.885 -3.281 -5.569) 7.078 1.383 ( -3.730 -2.368 4.445) 6.267 1.691 ( -6.543 -5.329 -7.440) 11.249 1.889 ( -0.926 -1.905 -7.477) 7.772 2.031 ( 0.959 -6.030 -1.599) 6.312 2.163 ( 0.404 -2.876 -2.082) 3.573 2.257 ( 1.389 -1.535 1.430) 2.516 2.279 ( 0.795 -2.409 1.711) 3.060 2.407 ( -0.013 0.818 5.970) 6.026 2.504 ( -2.082 2.153 -5.871) 6.591 2.536 ( -0.384 -1.464 -8.226) 8.364 2.569 ( -1.785 -1.551 3.622) 4.325 2.620 ( -3.584 -1.788 6.201) 7.382 2.852 ( 0.915 -0.670 4.559) 4.698 2.909 ( -2.484 1.624 4.052) 5.023 2.968 ( -1.154 -0.996 -2.498) 2.926 3.129 ( 1.687 4.921 -0.343) 5.214 3.289 ( -1.063 -1.207 1.413) 2.141 3.402 ( 3.531 -1.389 5.451) 6.642 3.428 ( 2.218 0.729 9.076) 9.371 3.600 ( -3.287 0.338 -1.770) 3.749 3.719 ( 2.664 -0.042 0.649) 2.743 3.780 ( 0.242 3.913 -3.499) 5.255 3.816 ( -3.203 0.808 -0.932) 3.432 4.022 ( -0.546 -1.770 -2.145) 2.834 4.125 ( -3.641 1.977 -2.076) 4.634 4.160 ( 6.189 2.864 3.252) 7.555 4.354 ( 7.028 7.204 -17.403) 20.103 6.261 ( 0.081 0.113 2.677) 2.680 6.287 ( -0.348 -0.525 -0.399) 0.746 7.225 ( 0.390 0.523 -4.155) 4.206 7.351 ( -3.150 -5.466 8.005) 10.193 ======================= Grid point 10 (7/44) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.741 ( 6.410 5.182 7.194) 10.941 0.835 ( 6.870 7.416 7.478) 12.574 1.047 ( -5.048 -8.619 -2.995) 10.428 1.130 ( -6.920 -3.557 -5.493) 9.524 1.324 ( 0.852 2.399 4.240) 4.946 1.351 ( -3.054 -0.271 4.113) 5.130 1.549 ( 2.826 1.080 1.081) 3.213 1.745 ( -4.910 -2.227 -9.112) 10.587 2.030 ( 4.588 -4.094 -0.897) 6.214 2.107 ( 1.751 -6.259 -4.735) 8.041 2.271 ( 1.000 2.243 0.299) 2.474 2.316 ( 0.592 -1.263 2.466) 2.833 2.407 ( -1.599 -2.443 -4.118) 5.048 2.462 ( 0.083 0.892 -1.421) 1.680 2.490 ( -0.705 2.719 0.382) 2.835 2.562 ( -1.673 -1.550 5.753) 6.188 2.641 ( -1.004 3.287 5.641) 6.606 2.811 ( -8.951 1.110 1.074) 9.083 2.925 ( 0.357 -1.012 -2.534) 2.752 2.930 ( 1.130 -0.986 2.001) 2.501 3.177 ( 1.245 2.099 -1.128) 2.689 3.280 ( -1.252 -1.571 3.282) 3.848 3.473 ( 1.049 -1.243 0.310) 1.655 3.504 ( 1.326 -0.426 0.577) 1.508 3.591 ( 1.284 1.835 7.063) 7.410 3.719 ( -4.245 0.084 -1.815) 4.618 3.759 ( 1.558 1.191 -5.809) 6.131 3.818 ( 0.255 1.309 3.069) 3.346 3.993 ( 0.637 -0.070 -1.972) 2.074 4.060 ( -3.322 1.276 -2.663) 4.444 4.329 ( 6.937 2.195 4.922) 8.785 4.368 ( 3.182 1.375 -10.020) 10.602 6.273 ( -0.410 -0.227 -0.721) 0.860 6.283 ( -0.051 -0.012 2.510) 2.511 7.204 ( 0.744 0.532 -4.783) 4.870 7.320 ( -2.468 -2.671 5.204) 6.349 ======================= Grid point 11 (8/44) ======================= q-point: (-0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.843 ( -7.267 -7.219 -3.114) 10.706 0.914 ( -9.610 -3.345 -5.600) 11.615 0.941 ( 7.531 1.905 5.480) 9.506 1.056 ( 6.205 5.328 3.854) 9.042 1.257 ( -1.283 -0.889 -3.885) 4.187 1.374 ( -1.876 1.224 4.603) 5.119 1.572 ( -4.625 -1.841 -5.683) 7.555 1.754 ( 7.381 4.468 7.668) 11.543 2.045 ( 2.433 -5.525 -3.395) 6.927 2.071 ( 5.265 -3.036 -3.128) 6.836 2.289 ( -1.198 -1.358 -1.314) 2.238 2.311 ( -0.601 -0.713 0.756) 1.200 2.342 ( 1.849 3.217 -7.430) 8.305 2.473 ( -1.460 -0.467 2.761) 3.158 2.527 ( -1.258 2.819 5.488) 6.296 2.587 ( -0.508 -1.945 7.327) 7.598 2.662 ( -6.886 0.355 -2.538) 7.347 2.718 ( 0.233 5.540 3.603) 6.613 2.889 ( -0.137 -1.287 -4.470) 4.653 2.952 ( 1.114 -1.247 1.023) 1.960 3.199 ( 1.082 1.136 -1.547) 2.203 3.264 ( -2.096 -1.839 3.576) 4.535 3.414 ( -1.390 -0.110 -6.246) 6.400 3.560 ( 0.807 -1.341 2.625) 3.057 3.605 ( -5.679 -1.133 -0.280) 5.798 3.689 ( 3.071 1.377 2.603) 4.255 3.739 ( -0.035 1.163 -4.877) 5.014 3.870 ( 0.993 -0.311 4.757) 4.870 4.001 ( 0.321 1.301 -2.032) 2.434 4.033 ( 0.446 2.072 0.584) 2.199 4.318 ( -1.296 -0.867 -6.336) 6.525 4.490 ( 5.740 0.605 6.242) 8.501 6.257 ( -0.541 -0.135 -0.987) 1.134 6.298 ( -0.230 -0.105 2.016) 2.032 7.183 ( 0.961 0.420 -5.253) 5.357 7.303 ( -1.708 -1.226 3.341) 3.948 ======================= Grid point 12 (9/44) ======================= q-point: (-0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.628 ( -9.852 -4.741 -2.056) 11.124 0.677 ( -11.492 -2.398 -4.009) 12.405 1.096 ( 5.757 -0.465 2.168) 6.170 1.152 ( -3.460 0.243 -4.444) 5.637 1.244 ( -3.787 1.043 -1.256) 4.124 1.302 ( 3.794 1.452 2.096) 4.572 1.612 ( 7.078 2.254 3.353) 8.150 1.928 ( 3.593 2.350 3.228) 5.371 2.026 ( 3.495 -4.132 0.551) 5.440 2.099 ( 4.863 -0.987 -5.534) 7.433 2.236 ( -1.694 -0.102 -1.644) 2.363 2.330 ( 1.614 -1.756 5.675) 6.156 2.370 ( 2.945 2.778 -5.506) 6.834 2.404 ( -6.289 -0.388 -0.396) 6.314 2.536 ( -0.859 -0.581 -3.152) 3.318 2.588 ( -0.250 2.117 5.778) 6.159 2.684 ( 1.391 -0.522 7.148) 7.301 2.782 ( -0.625 4.008 0.315) 4.069 2.824 ( -0.675 -1.740 -6.501) 6.764 2.971 ( 1.881 -0.510 0.782) 2.100 3.212 ( -0.055 0.215 0.392) 0.451 3.244 ( -1.579 0.346 2.228) 2.753 3.342 ( -2.005 -0.063 -5.238) 5.609 3.423 ( -6.273 -3.569 -4.689) 8.607 3.631 ( -0.224 0.575 5.032) 5.070 3.705 ( 0.039 0.341 -4.159) 4.173 3.767 ( 3.250 0.809 1.374) 3.620 3.887 ( -2.775 -0.077 2.691) 3.866 4.028 ( 3.079 -0.172 1.108) 3.277 4.074 ( 0.943 1.647 -0.821) 2.068 4.222 ( -4.481 -0.727 -2.651) 5.257 4.627 ( 4.630 0.000 6.944) 8.346 6.239 ( -0.611 -0.110 -0.818) 1.027 6.307 ( -0.226 -0.079 1.402) 1.423 7.158 ( 0.724 0.172 -5.288) 5.340 7.293 ( -1.186 -0.417 1.715) 2.126 ======================= Grid point 13 (10/44) ======================= q-point: (-0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.418 ( -11.995 0.000 0.000) 11.995 0.448 ( -13.152 0.000 0.000) 13.152 0.826 ( -22.946 0.000 0.000) 22.946 1.243 ( 8.780 -0.000 -0.000) 8.780 1.316 ( 5.312 -0.000 -0.000) 5.312 1.393 ( 4.050 -0.000 -0.000) 4.050 1.782 ( 9.553 -0.000 -0.000) 9.553 1.966 ( 1.187 -0.000 -0.000) 1.187 2.082 ( 3.640 -0.000 -0.000) 3.640 2.149 ( 4.086 -0.000 -0.000) 4.086 2.211 ( -1.056 0.000 0.000) 1.056 2.283 ( -5.108 0.000 0.000) 5.108 2.418 ( 4.279 -0.000 -0.000) 4.279 2.427 ( 5.618 -0.000 -0.000) 5.618 2.499 ( -1.657 0.000 0.000) 1.657 2.644 ( 1.671 -0.000 -0.000) 1.671 2.737 ( -1.821 0.000 0.000) 1.821 2.740 ( -0.386 0.000 0.000) 0.386 2.780 ( -1.661 0.000 0.000) 1.661 3.008 ( 1.893 -0.000 -0.000) 1.893 3.158 ( -7.630 0.000 0.000) 7.630 3.251 ( -0.548 0.000 0.000) 0.548 3.303 ( 0.021 -0.000 -0.000) 0.021 3.307 ( -0.553 0.000 0.000) 0.553 3.646 ( -1.033 0.000 0.000) 1.033 3.691 ( 0.582 -0.000 -0.000) 0.582 3.827 ( 2.621 -0.000 -0.000) 2.621 3.830 ( -4.986 0.000 0.000) 4.986 4.045 ( -3.821 0.000 0.000) 3.821 4.094 ( 3.390 -0.000 -0.000) 3.390 4.157 ( -1.769 0.000 0.000) 1.769 4.738 ( 4.136 -0.000 -0.000) 4.136 6.226 ( -0.465 0.000 0.000) 0.465 6.310 ( -0.122 0.000 0.000) 0.122 7.135 ( -0.077 0.000 0.000) 0.077 7.282 ( -0.811 0.000 0.000) 0.811 ======================= Grid point 16 (11/44) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.684 ( 0.000 9.541 10.023) 13.838 0.690 ( 0.000 9.459 7.563) 12.110 1.133 ( -0.000 -9.051 -3.197) 9.599 1.193 ( 0.000 11.513 13.698) 17.894 1.238 ( -0.000 -3.938 -5.456) 6.729 1.426 ( 0.000 -3.313 5.531) 6.448 1.638 ( -0.000 -5.465 -9.492) 10.953 1.814 ( -0.000 -1.740 -7.202) 7.409 1.943 ( -0.000 -7.131 -1.557) 7.299 2.055 ( -0.000 -7.158 -5.991) 9.335 2.265 ( 0.000 -1.324 2.168) 2.541 2.278 ( 0.000 0.795 3.154) 3.252 2.443 ( -0.000 -4.086 -0.525) 4.119 2.487 ( -0.000 0.980 -7.114) 7.181 2.497 ( 0.000 3.833 1.134) 3.997 2.590 ( 0.000 -0.246 1.678) 1.696 2.687 ( 0.000 -2.051 7.636) 7.906 2.845 ( 0.000 -1.206 2.807) 3.056 2.965 ( -0.000 -0.102 -0.809) 0.815 2.996 ( 0.000 3.163 2.683) 4.148 3.165 ( -0.000 4.047 -0.208) 4.053 3.281 ( -0.000 -1.505 -0.460) 1.573 3.407 ( 0.000 -1.715 7.319) 7.517 3.478 ( 0.000 -0.213 5.674) 5.678 3.630 ( 0.000 -0.557 2.663) 2.720 3.700 ( -0.000 2.645 -2.973) 3.980 3.800 ( 0.000 2.984 0.485) 3.023 3.847 ( -0.000 1.025 -1.587) 1.889 3.987 ( -0.000 -2.911 -3.687) 4.698 4.143 ( 0.000 3.630 1.917) 4.105 4.191 ( 0.000 1.212 2.385) 2.676 4.219 ( -0.000 7.651 -16.884) 18.536 6.283 ( 0.000 -0.009 2.823) 2.823 6.283 ( 0.000 -0.612 -0.029) 0.613 7.193 ( -0.000 0.511 -4.338) 4.368 7.373 ( 0.000 -4.935 4.827) 6.903 ======================= Grid point 17 (12/44) ======================= q-point: ( 0.43 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.819 ( 0.084 7.295 7.920) 10.768 0.873 ( 3.402 6.564 6.960) 10.154 1.012 ( -2.629 -6.861 -2.202) 7.670 1.136 ( -3.379 -3.507 -4.732) 6.790 1.379 ( -2.176 1.601 1.608) 3.144 1.426 ( -2.161 -2.463 5.575) 6.466 1.509 ( 3.502 0.518 -0.491) 3.574 1.720 ( -2.410 -1.023 -8.786) 9.168 1.911 ( 3.684 -5.214 -0.029) 6.384 1.970 ( 2.078 -5.922 -4.774) 7.885 2.278 ( 0.717 0.205 0.957) 1.213 2.311 ( 0.546 -0.450 1.777) 1.913 2.374 ( -1.275 -2.320 -3.926) 4.735 2.437 ( 0.272 0.115 -4.491) 4.500 2.558 ( -0.596 2.813 4.793) 5.590 2.599 ( -0.535 0.556 1.843) 1.998 2.734 ( -0.483 1.277 6.997) 7.129 2.876 ( 2.608 -1.199 1.270) 3.139 2.928 ( -3.058 -0.600 -1.122) 3.312 2.993 ( -6.327 2.837 1.239) 7.043 3.183 ( -0.155 1.553 -0.782) 1.746 3.277 ( 1.293 -1.883 0.621) 2.367 3.457 ( -0.421 -2.214 2.258) 3.190 3.482 ( 1.098 -1.005 2.057) 2.539 3.662 ( -1.496 0.529 6.495) 6.686 3.702 ( 1.492 2.258 -4.831) 5.537 3.808 ( -4.346 2.046 1.044) 4.916 3.842 ( 0.281 0.760 0.909) 1.218 3.952 ( 1.729 -0.147 -1.997) 2.646 4.095 ( -3.761 2.299 -7.754) 8.919 4.248 ( 5.565 2.710 -2.799) 6.793 4.276 ( 5.136 0.926 0.706) 5.266 6.275 ( -0.429 -0.302 -0.049) 0.527 6.303 ( -0.071 -0.101 2.262) 2.265 7.167 ( 0.188 0.390 -3.800) 3.825 7.355 ( -1.159 -2.843 2.972) 4.273 ======================= Grid point 18 (13/44) ======================= q-point: ( 0.57 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.865 ( -4.407 -4.236 -1.222) 6.234 0.952 ( 0.927 3.618 7.066) 7.992 0.985 ( -6.404 -1.612 -3.278) 7.373 1.035 ( 2.803 2.111 1.886) 3.984 1.239 ( -2.706 -2.397 -5.216) 6.346 1.449 ( -0.230 0.392 6.527) 6.542 1.584 ( -4.265 -0.855 -6.306) 7.661 1.714 ( 8.395 3.707 4.269) 10.122 1.929 ( 2.575 -4.248 -2.211) 5.438 1.937 ( 4.357 -2.741 -2.224) 5.607 2.291 ( 0.511 -0.591 -1.762) 1.928 2.298 ( -1.333 -0.938 0.645) 1.753 2.323 ( -0.550 1.074 -2.383) 2.671 2.456 ( 3.599 -0.163 0.122) 3.605 2.590 ( -0.341 0.239 2.456) 2.492 2.634 ( 0.162 2.198 4.586) 5.088 2.788 ( -1.538 1.509 1.770) 2.788 2.828 ( -8.594 0.420 0.516) 8.620 2.880 ( -5.552 0.590 -1.354) 5.745 2.919 ( 1.721 -1.275 0.092) 2.144 3.191 ( 0.483 1.102 -0.750) 1.418 3.281 ( 0.202 -2.381 1.586) 2.868 3.400 ( -1.870 -0.353 -3.130) 3.663 3.540 ( 2.272 -1.490 3.016) 4.060 3.670 ( 0.450 1.021 -3.279) 3.463 3.683 ( -3.606 -0.112 1.000) 3.744 3.778 ( -4.018 2.006 3.734) 5.840 3.858 ( 1.810 -0.482 0.054) 1.874 3.990 ( 1.454 1.960 -1.556) 2.894 4.028 ( 0.562 0.054 -2.502) 2.565 4.287 ( 1.000 0.850 -2.465) 2.793 4.411 ( 7.958 0.032 2.244) 8.269 6.263 ( -0.766 -0.099 -0.209) 0.800 6.314 ( -0.186 -0.091 1.261) 1.278 7.147 ( 0.137 0.172 -2.553) 2.562 7.335 ( -1.273 -1.072 1.537) 2.266 ======================= Grid point 19 (14/44) ======================= q-point: (-0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.739 ( -7.623 0.000 0.000) 7.623 0.822 ( -9.920 0.000 0.000) 9.920 1.023 ( 2.522 -0.000 -0.000) 2.522 1.101 ( 4.915 -0.000 -0.000) 4.915 1.240 ( 3.895 -0.000 -0.000) 3.895 1.413 ( -10.619 0.000 0.000) 10.619 1.538 ( 4.674 -0.000 -0.000) 4.674 1.885 ( 6.306 -0.000 -0.000) 6.306 1.949 ( 4.975 -0.000 -0.000) 4.975 1.972 ( 3.531 -0.000 -0.000) 3.531 2.257 ( -1.836 0.000 0.000) 1.836 2.309 ( 1.967 -0.000 -0.000) 1.967 2.324 ( 0.914 -0.000 -0.000) 0.914 2.506 ( 1.723 -0.000 -0.000) 1.723 2.596 ( -13.761 0.000 0.000) 13.761 2.609 ( 1.053 -0.000 -0.000) 1.053 2.696 ( 2.158 -0.000 -0.000) 2.158 2.796 ( -2.258 0.000 0.000) 2.258 2.826 ( -1.336 0.000 0.000) 1.336 2.941 ( 1.660 -0.000 -0.000) 1.660 3.207 ( 0.926 -0.000 -0.000) 0.926 3.269 ( -0.688 0.000 0.000) 0.688 3.353 ( -2.194 0.000 0.000) 2.194 3.509 ( -6.450 0.000 0.000) 6.450 3.668 ( 0.758 0.000 -0.000) 0.758 3.668 ( 0.758 0.000 -0.000) 0.758 3.750 ( -1.193 0.000 0.000) 1.193 3.904 ( 2.655 -0.000 -0.000) 2.655 3.995 ( 0.383 -0.000 -0.000) 0.383 4.056 ( 1.863 -0.000 -0.000) 1.863 4.268 ( -2.377 0.000 0.000) 2.377 4.548 ( 7.509 -0.000 -0.000) 7.509 6.248 ( -0.846 0.000 0.000) 0.846 6.315 ( -0.199 0.000 0.000) 0.199 7.137 ( -0.046 0.000 0.000) 0.046 7.317 ( -1.256 0.000 0.000) 1.256 ======================= Grid point 24 (15/44) ======================= q-point: ( 0.43 0.43 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.913 ( 0.000 4.175 3.204) 5.263 0.925 ( 0.000 -2.635 0.895) 2.783 1.028 ( 0.000 4.934 5.028) 7.044 1.114 ( -0.000 -0.819 -0.889) 1.208 1.370 ( -0.000 -6.777 -9.398) 11.587 1.458 ( 0.000 -1.998 4.588) 5.005 1.494 ( 0.000 3.855 5.407) 6.640 1.671 ( -0.000 0.000 -6.528) 6.528 1.841 ( 0.000 -3.194 0.401) 3.219 1.884 ( -0.000 -3.125 -2.109) 3.770 2.270 ( 0.000 -0.667 0.041) 0.669 2.321 ( 0.000 -0.065 1.422) 1.423 2.346 ( -0.000 -1.077 -1.565) 1.900 2.397 ( -0.000 -0.028 -4.005) 4.005 2.617 ( 0.000 0.483 0.750) 0.892 2.617 ( 0.000 1.656 3.079) 3.496 2.792 ( 0.000 0.593 4.768) 4.805 2.853 ( 0.000 -0.162 0.772) 0.789 2.952 ( -0.000 -0.512 -0.597) 0.786 3.077 ( -0.000 3.051 0.294) 3.065 3.185 ( -0.000 -0.074 -0.710) 0.714 3.244 ( -0.000 -1.652 -0.398) 1.699 3.467 ( -0.000 -1.307 -0.826) 1.546 3.486 ( 0.000 -1.633 5.794) 6.020 3.658 ( -0.000 1.642 -5.213) 5.465 3.750 ( 0.000 0.193 6.415) 6.418 3.838 ( -0.000 0.352 -0.819) 0.892 3.880 ( 0.000 1.808 1.664) 2.457 3.927 ( 0.000 0.256 1.760) 1.779 4.072 ( -0.000 1.420 -11.863) 11.948 4.210 ( -0.000 2.827 -0.273) 2.841 4.240 ( 0.000 -0.990 2.062) 2.287 6.278 ( 0.000 -0.149 0.105) 0.182 6.317 ( 0.000 -0.096 1.312) 1.316 7.145 ( -0.000 0.224 -2.129) 2.141 7.359 ( 0.000 -1.618 1.308) 2.081 ======================= Grid point 25 (16/44) ======================= q-point: (-0.43 -0.57 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.892 ( -2.076 0.000 0.000) 2.076 0.969 ( 2.769 -0.000 -0.000) 2.769 1.037 ( -2.986 0.000 0.000) 2.986 1.109 ( -5.207 0.000 0.000) 5.207 1.232 ( -3.107 0.000 0.000) 3.107 1.472 ( 0.459 -0.000 -0.000) 0.459 1.588 ( 2.998 -0.000 -0.000) 2.998 1.638 ( 1.801 -0.000 -0.000) 1.801 1.855 ( 3.083 -0.000 -0.000) 3.083 1.871 ( 1.244 -0.000 -0.000) 1.244 2.272 ( 0.631 -0.000 -0.000) 0.631 2.323 ( -1.479 0.000 0.000) 1.479 2.324 ( -0.314 0.000 0.000) 0.314 2.399 ( 2.486 -0.000 -0.000) 2.486 2.610 ( -0.952 0.000 0.000) 0.952 2.652 ( 1.133 -0.000 -0.000) 1.133 2.813 ( -0.368 0.000 0.000) 0.368 2.876 ( 2.128 -0.000 -0.000) 2.128 2.908 ( -3.821 0.000 0.000) 3.821 3.039 ( -7.657 0.000 0.000) 7.657 3.184 ( 0.455 -0.000 -0.000) 0.455 3.252 ( 1.985 -0.000 -0.000) 1.985 3.425 ( -2.967 0.000 0.000) 2.967 3.515 ( 1.483 -0.000 -0.000) 1.483 3.651 ( 0.613 -0.000 -0.000) 0.613 3.746 ( -3.056 0.000 0.000) 3.056 3.832 ( 0.013 -0.000 -0.000) 0.013 3.889 ( -2.892 0.000 0.000) 2.892 3.973 ( 1.509 -0.000 -0.000) 1.509 3.989 ( -0.252 0.000 0.000) 0.252 4.254 ( 2.574 -0.000 -0.000) 2.574 4.297 ( 6.380 -0.000 -0.000) 6.380 6.273 ( -0.536 0.000 0.000) 0.536 6.321 ( -0.110 0.000 0.000) 0.110 7.137 ( -0.003 0.000 0.000) 0.003 7.351 ( -0.662 0.000 0.000) 0.662 ======================= Grid point 59 (17/44) ======================= q-point: ( 0.14 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 44 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.343 ( 0.000 -0.000 14.886) 14.886 0.343 ( 0.000 -0.000 14.886) 14.886 0.690 ( 0.000 -0.000 26.354) 26.354 1.343 ( -0.000 0.000 -8.951) 8.951 1.434 ( 0.000 -0.000 1.847) 1.847 1.434 ( 0.000 -0.000 1.847) 1.847 1.876 ( -0.000 0.000 -6.719) 6.719 1.876 ( -0.000 0.000 -6.719) 6.719 2.189 ( -0.000 0.000 -5.494) 5.494 2.189 ( -0.000 0.000 -5.494) 5.494 2.266 ( 0.000 -0.000 5.985) 5.985 2.266 ( 0.000 -0.000 5.985) 5.985 2.446 ( -0.000 0.000 -3.878) 3.878 2.446 ( -0.000 0.000 -3.878) 3.878 2.622 ( 0.000 -0.000 0.533) 0.533 2.622 ( 0.000 -0.000 0.533) 0.533 2.862 ( 0.000 -0.000 5.539) 5.539 2.862 ( 0.000 -0.000 5.539) 5.539 2.901 ( -0.000 0.000 -7.935) 7.935 2.986 ( -0.000 0.000 -2.360) 2.360 3.025 ( 0.000 -0.000 0.516) 0.516 3.355 ( 0.000 -0.000 2.025) 2.025 3.355 ( 0.000 -0.000 2.025) 2.025 3.457 ( 0.000 -0.000 9.728) 9.728 3.647 ( 0.000 -0.000 1.543) 1.543 3.647 ( 0.000 -0.000 1.543) 1.543 3.694 ( -0.000 0.000 -2.598) 2.598 3.828 ( -0.000 0.000 -3.055) 3.055 3.828 ( -0.000 0.000 -3.055) 3.055 4.005 ( 0.000 -0.000 4.908) 4.908 4.127 ( -0.000 0.000 -1.963) 1.963 4.127 ( -0.000 0.000 -1.963) 1.963 6.259 ( 0.000 -0.000 3.074) 3.074 6.311 ( 0.000 -0.000 0.128) 0.128 7.215 ( -0.000 0.000 -3.917) 3.917 7.561 ( -0.000 0.000 -1.641) 1.641 ======================= Grid point 60 (18/44) ======================= q-point: ( 0.29 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.509 ( 4.776 2.758 12.843) 13.977 0.546 ( 6.559 3.787 12.831) 14.899 0.992 ( 9.123 5.267 19.637) 22.284 1.243 ( -2.482 -1.433 -9.319) 9.750 1.377 ( -4.861 -2.806 1.131) 5.726 1.443 ( -2.483 -1.433 3.792) 4.753 1.717 ( -6.351 -3.667 -7.642) 10.592 1.851 ( -0.947 -0.547 -8.456) 8.526 2.059 ( -1.780 -1.028 -5.997) 6.339 2.162 ( -0.070 -0.041 -7.894) 7.895 2.294 ( -0.913 -0.527 0.683) 1.256 2.335 ( 0.914 0.528 5.946) 6.039 2.424 ( 0.576 0.333 3.257) 3.325 2.440 ( 0.244 0.141 -2.110) 2.129 2.552 ( -4.215 -2.434 -1.487) 5.090 2.643 ( -0.789 -0.456 4.613) 4.702 2.697 ( -7.916 -4.570 0.812) 9.177 2.898 ( -1.645 -0.950 3.796) 4.245 2.933 ( 0.521 0.301 2.759) 2.824 2.963 ( -0.899 -0.519 -2.121) 2.362 3.060 ( 3.307 1.909 1.589) 4.136 3.331 ( -1.405 -0.811 2.158) 2.699 3.427 ( 3.122 1.802 3.802) 5.240 3.518 ( 0.295 0.171 10.818) 10.823 3.590 ( -3.586 -2.070 -1.215) 4.315 3.695 ( 1.648 0.951 -3.941) 4.376 3.738 ( 2.742 1.583 1.746) 3.616 3.770 ( -1.436 -0.829 -3.925) 4.261 3.986 ( 9.323 5.383 -10.135) 14.786 4.037 ( -1.387 -0.801 -0.954) 1.864 4.136 ( 4.742 2.738 2.428) 5.990 4.155 ( 4.465 2.578 -3.431) 6.193 6.284 ( -0.034 -0.020 3.215) 3.215 6.300 ( -0.924 -0.533 0.113) 1.073 7.186 ( 0.279 0.161 -4.066) 4.078 7.480 ( -5.067 -2.925 1.399) 6.015 ======================= Grid point 61 (19/44) ======================= q-point: ( 0.43 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.723 ( 6.399 3.695 10.017) 12.448 0.792 ( 7.497 4.328 10.186) 13.367 1.095 ( -4.997 -2.885 -7.511) 9.472 1.219 ( -7.151 -4.128 2.042) 8.505 1.317 ( 3.120 1.801 6.030) 7.025 1.413 ( -3.142 -1.814 5.773) 6.818 1.542 ( -1.218 -0.703 -2.810) 3.142 1.713 ( -4.012 -2.316 -9.712) 10.760 2.032 ( 2.059 1.189 -4.607) 5.185 2.090 ( 0.684 0.395 -11.538) 11.565 2.277 ( -1.187 -0.685 2.717) 3.044 2.359 ( -1.394 -0.805 2.212) 2.736 2.419 ( -3.183 -1.838 -4.204) 5.585 2.464 ( 1.504 0.868 2.953) 3.426 2.494 ( -0.892 -0.515 5.928) 6.017 2.608 ( -1.545 -0.892 1.662) 2.438 2.669 ( -0.798 -0.461 6.700) 6.763 2.852 ( -6.282 -3.627 4.078) 8.322 2.911 ( -1.001 -0.578 -3.530) 3.714 2.958 ( 0.220 0.127 1.855) 1.872 3.135 ( 2.642 1.525 0.105) 3.052 3.329 ( -0.854 -0.493 2.608) 2.788 3.484 ( -4.088 -2.360 -3.677) 5.984 3.529 ( 2.987 1.724 4.393) 5.585 3.624 ( 1.390 0.803 9.590) 9.724 3.682 ( 0.359 0.207 -5.770) 5.785 3.707 ( -1.204 -0.695 -4.391) 4.606 3.815 ( 2.723 1.572 3.949) 5.048 3.983 ( -1.030 -0.595 -3.102) 3.322 4.013 ( 0.982 0.567 -4.825) 4.957 4.239 ( 5.825 3.363 -7.659) 10.193 4.304 ( 7.434 4.292 3.656) 9.330 6.279 ( -0.957 -0.552 -0.040) 1.106 6.305 ( -0.205 -0.118 2.533) 2.544 7.163 ( 0.396 0.229 -3.591) 3.620 7.398 ( -3.749 -2.164 1.872) 4.716 ======================= Grid point 62 (20/44) ======================= q-point: ( 0.57 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.904 ( -3.161 -1.825 -1.922) 4.125 0.948 ( 2.740 1.582 5.549) 6.388 0.989 ( -2.544 -1.469 3.135) 4.296 1.023 ( -2.980 -1.721 1.983) 3.972 1.254 ( -2.968 -1.714 -3.650) 5.007 1.411 ( -2.483 -1.433 5.197) 5.935 1.560 ( -3.058 -1.765 -5.799) 6.789 1.714 ( 8.990 5.190 5.051) 11.544 2.018 ( 1.780 1.028 -12.560) 12.727 2.045 ( 2.552 1.474 -5.426) 6.175 2.277 ( -0.316 -0.182 -0.796) 0.875 2.318 ( -3.300 -1.905 -0.005) 3.811 2.354 ( -0.733 -0.423 3.579) 3.677 2.512 ( 0.170 0.098 3.412) 3.417 2.542 ( 0.449 0.259 7.281) 7.299 2.626 ( 1.963 1.133 -0.200) 2.275 2.669 ( -6.287 -3.630 -1.637) 7.442 2.746 ( -3.268 -1.887 11.127) 11.749 2.865 ( -0.672 -0.388 -4.839) 4.901 2.970 ( 0.101 0.058 0.472) 0.486 3.176 ( 1.307 0.755 0.047) 1.510 3.322 ( -1.569 -0.906 1.651) 2.451 3.378 ( -2.890 -1.668 -3.409) 4.770 3.599 ( -4.168 -2.406 -4.315) 6.464 3.619 ( 1.386 0.800 6.004) 6.213 3.666 ( 1.291 0.745 -6.570) 6.737 3.707 ( 1.167 0.674 4.807) 4.992 3.904 ( 1.448 0.836 3.905) 4.248 3.965 ( 1.414 0.816 -2.192) 2.733 4.009 ( 1.321 0.762 0.308) 1.556 4.260 ( 1.410 0.814 -6.601) 6.799 4.483 ( 6.726 3.883 2.375) 8.122 6.259 ( -0.838 -0.484 -0.213) 0.991 6.315 ( -0.308 -0.178 1.388) 1.433 7.146 ( 0.247 0.143 -2.527) 2.543 7.345 ( -2.301 -1.329 0.923) 2.813 ======================= Grid point 63 (21/44) ======================= q-point: (-0.29 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.724 ( -9.070 -5.237 0.105) 10.474 0.761 ( -11.652 -6.728 0.436) 13.462 1.086 ( 3.287 1.898 -3.192) 4.959 1.125 ( 4.420 2.552 5.680) 7.636 1.268 ( -9.599 -5.542 -1.378) 11.169 1.276 ( 4.303 2.485 1.085) 5.086 1.601 ( 5.190 2.996 5.115) 7.879 1.879 ( 3.905 2.255 -3.498) 5.707 2.004 ( 5.817 3.358 -6.290) 9.202 2.050 ( 2.341 1.351 -6.251) 6.810 2.235 ( -1.336 -0.771 -0.759) 1.719 2.308 ( 0.668 0.386 0.958) 1.230 2.391 ( 0.937 0.541 3.902) 4.049 2.463 ( -8.650 -4.994 2.256) 10.240 2.539 ( -0.553 -0.319 -0.668) 0.924 2.610 ( 0.179 0.104 6.029) 6.033 2.677 ( 3.049 1.760 -0.836) 3.618 2.775 ( -0.766 -0.442 -3.817) 3.918 2.829 ( -1.077 -0.622 6.713) 6.828 2.977 ( 0.676 0.390 -0.173) 0.800 3.204 ( 1.120 0.647 1.413) 1.915 3.289 ( -1.306 -0.754 -2.374) 2.813 3.306 ( -3.033 -1.751 1.062) 3.660 3.453 ( -5.983 -3.454 -3.554) 7.769 3.655 ( 1.753 1.012 -4.832) 5.239 3.671 ( -0.945 -0.545 5.572) 5.678 3.769 ( 2.355 1.359 1.990) 3.369 3.914 ( -1.655 -0.956 1.074) 2.192 4.003 ( 3.216 1.857 -1.679) 4.075 4.073 ( 1.623 0.937 4.056) 4.467 4.204 ( -2.218 -1.280 -4.127) 4.857 4.621 ( 5.393 3.114 -1.144) 6.331 6.241 ( -0.700 -0.404 0.058) 0.810 6.314 ( -0.234 -0.135 -0.031) 0.272 7.136 ( -0.063 -0.036 0.171) 0.186 7.310 ( -1.333 -0.770 -0.394) 1.589 ======================= Grid point 64 (22/44) ======================= q-point: (-0.14 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.522 ( -11.629 -6.714 3.215) 13.808 0.530 ( -10.592 -6.115 4.733) 13.115 0.939 ( -17.008 -9.819 4.602) 20.171 1.183 ( 7.651 4.417 -2.384) 9.150 1.310 ( 4.756 2.746 3.886) 6.727 1.363 ( 3.630 2.096 -2.316) 4.788 1.789 ( 7.288 4.208 5.886) 10.269 1.873 ( 2.396 1.383 -9.343) 9.744 2.058 ( 3.389 1.957 -6.049) 7.204 2.122 ( 4.710 2.719 0.763) 5.492 2.236 ( 1.295 0.748 1.963) 2.467 2.287 ( -5.299 -3.059 -3.020) 6.824 2.344 ( 2.432 1.404 -0.401) 2.837 2.489 ( 2.315 1.337 8.073) 8.504 2.528 ( -0.948 -0.547 2.521) 2.748 2.633 ( 1.222 0.705 -4.052) 4.291 2.692 ( -1.451 -0.838 -8.813) 8.971 2.727 ( 0.146 0.084 4.989) 4.992 2.867 ( -1.972 -1.138 8.198) 8.509 2.997 ( 1.323 0.764 -1.221) 1.956 3.194 ( -5.809 -3.354 -2.635) 7.207 3.261 ( 2.235 1.291 4.072) 4.821 3.283 ( -1.430 -0.826 5.012) 5.277 3.336 ( -2.119 -1.223 -1.546) 2.894 3.659 ( 1.450 0.837 -2.971) 3.410 3.674 ( -1.632 -0.942 2.865) 3.429 3.821 ( 0.510 0.294 0.295) 0.658 3.865 ( -2.638 -1.523 0.347) 3.066 4.060 ( 3.250 1.877 -2.178) 4.339 4.081 ( -3.692 -2.132 1.610) 4.557 4.155 ( -0.098 -0.056 0.883) 0.891 4.658 ( 1.551 0.895 -13.781) 13.897 6.233 ( -0.429 -0.248 0.988) 1.105 6.304 ( -0.152 -0.088 -1.440) 1.450 7.171 ( 1.371 0.791 8.352) 8.500 7.279 ( -0.768 -0.443 -1.851) 2.052 ======================= Grid point 67 (23/44) ======================= q-point: ( 0.29 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.627 ( 0.000 5.191 12.881) 13.887 0.642 ( 0.000 5.339 14.408) 15.366 1.129 ( 0.000 7.668 16.529) 18.221 1.176 ( -0.000 -2.304 -7.850) 8.181 1.362 ( 0.000 -8.201 1.173) 8.284 1.497 ( 0.000 0.021 4.051) 4.052 1.677 ( -0.000 -5.557 -7.379) 9.237 1.779 ( -0.000 -1.336 -9.180) 9.277 1.993 ( -0.000 -2.470 -7.565) 7.958 2.056 ( -0.000 -1.060 -12.720) 12.764 2.320 ( 0.000 -1.978 3.140) 3.710 2.352 ( 0.000 -2.237 5.638) 6.066 2.437 ( -0.000 0.586 -0.768) 0.966 2.454 ( 0.000 1.128 1.403) 1.800 2.554 ( 0.000 -5.513 3.467) 6.512 2.634 ( 0.000 -3.422 1.164) 3.614 2.776 ( 0.000 -1.180 2.846) 3.080 2.891 ( 0.000 -2.002 1.557) 2.536 2.957 ( -0.000 0.302 -0.637) 0.705 2.966 ( 0.000 1.509 1.908) 2.433 3.087 ( 0.000 4.064 2.061) 4.557 3.329 ( 0.000 -3.204 1.451) 3.517 3.454 ( 0.000 2.497 5.119) 5.696 3.568 ( 0.000 -3.012 7.556) 8.134 3.629 ( -0.000 0.627 -5.759) 5.793 3.652 ( 0.000 -0.065 2.278) 2.279 3.747 ( 0.000 3.585 2.980) 4.662 3.805 ( -0.000 1.908 -3.486) 3.974 3.887 ( -0.000 10.618 -8.773) 13.774 4.037 ( -0.000 -2.397 -4.596) 5.183 4.092 ( 0.000 2.319 0.644) 2.407 4.158 ( 0.000 3.766 0.960) 3.887 6.301 ( 0.000 -0.921 0.085) 0.925 6.309 ( 0.000 -0.179 2.877) 2.883 7.157 ( -0.000 0.191 -3.123) 3.129 7.480 ( -0.000 -4.039 -1.047) 4.173 ======================= Grid point 68 (24/44) ======================= q-point: ( 0.43 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.783 ( 1.882 4.614 11.616) 12.639 0.835 ( 3.493 5.177 12.040) 13.564 1.081 ( -2.661 -3.196 -4.279) 5.967 1.255 ( -3.337 -8.868 1.175) 9.548 1.347 ( 2.823 4.573 6.900) 8.746 1.500 ( -3.599 -1.619 3.211) 5.088 1.525 ( -1.998 -3.255 -3.113) 4.927 1.668 ( -2.111 -1.527 -10.104) 10.435 1.945 ( 2.436 -1.197 -6.823) 7.343 1.962 ( 1.862 -0.731 -11.724) 11.894 2.313 ( -1.077 -1.862 3.306) 3.944 2.366 ( 0.119 -2.893 2.283) 3.687 2.416 ( -1.838 -3.541 -1.512) 4.267 2.450 ( 1.591 -1.636 2.385) 3.300 2.566 ( -1.154 1.121 6.131) 6.339 2.613 ( -0.067 -2.595 1.624) 3.062 2.772 ( -3.277 1.737 5.949) 7.010 2.891 ( -0.323 -0.710 -1.920) 2.072 2.945 ( 0.304 -1.448 0.566) 1.584 2.966 ( -3.786 0.499 1.621) 4.149 3.139 ( 0.402 3.193 1.025) 3.378 3.324 ( 1.277 -2.375 1.416) 3.046 3.504 ( -4.424 -3.353 -2.034) 5.912 3.532 ( 1.814 0.795 4.921) 5.305 3.618 ( 1.870 0.829 -5.891) 6.235 3.700 ( -1.347 1.415 6.441) 6.731 3.757 ( -3.481 3.714 0.940) 5.177 3.824 ( 1.690 2.488 2.456) 3.883 3.938 ( 2.434 1.557 -4.599) 5.431 3.970 ( 0.151 1.266 -6.631) 6.752 4.153 ( 3.968 3.606 -3.100) 6.194 4.248 ( 5.590 3.586 1.410) 6.789 6.287 ( -0.693 -0.778 0.154) 1.053 6.323 ( -0.172 -0.326 1.710) 1.749 7.141 ( 0.080 0.119 -1.500) 1.507 7.421 ( -1.752 -2.931 -1.076) 3.580 ======================= Grid point 69 (25/44) ======================= q-point: ( 0.57 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.933 ( -0.048 0.744 5.608) 5.657 0.982 ( -1.258 -1.671 4.484) 4.948 1.008 ( 2.283 1.684 6.546) 7.134 1.086 ( -7.127 -8.230 0.104) 10.887 1.256 ( -4.705 -3.786 -4.836) 7.737 1.497 ( -1.232 -0.182 6.046) 6.173 1.542 ( -3.323 -0.510 -7.169) 7.919 1.651 ( 8.276 5.477 2.410) 10.212 1.882 ( 2.101 0.659 -12.787) 12.975 1.953 ( 3.450 -0.608 -5.422) 6.455 2.293 ( -1.850 -1.538 0.685) 2.502 2.342 ( -1.170 -2.753 1.488) 3.341 2.374 ( -0.319 -2.552 4.133) 4.868 2.488 ( 2.972 -2.020 1.970) 4.098 2.603 ( 0.552 0.871 -0.304) 1.075 2.624 ( -0.530 0.622 6.451) 6.502 2.777 ( -5.559 1.953 2.310) 6.329 2.831 ( -2.923 0.293 -3.946) 4.919 2.887 ( -5.690 -2.324 5.442) 8.209 2.952 ( 1.506 -1.892 -0.106) 2.421 3.173 ( 0.510 1.138 1.103) 1.665 3.334 ( -0.181 0.024 0.007) 0.183 3.392 ( -2.943 -3.538 -1.529) 4.849 3.558 ( 1.763 -0.914 -4.681) 5.084 3.634 ( -0.074 0.960 0.034) 0.963 3.678 ( -2.836 2.421 1.392) 3.980 3.778 ( -3.169 1.567 5.863) 6.847 3.870 ( 2.862 -1.559 -1.255) 3.493 3.944 ( 0.971 1.749 -1.847) 2.723 3.997 ( 3.123 2.221 1.053) 3.975 4.213 ( 1.981 3.771 -3.553) 5.547 4.394 ( 8.006 2.485 -1.060) 8.450 6.269 ( -0.944 -0.528 0.236) 1.107 6.324 ( -0.303 -0.329 0.198) 0.489 7.138 ( -0.027 -0.037 0.615) 0.617 7.362 ( -1.668 -1.567 -1.511) 2.742 ======================= Grid point 70 (26/44) ======================= q-point: (-0.29 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.833 ( -7.426 -5.389 3.770) 9.920 0.869 ( -9.504 -5.618 1.146) 11.100 1.028 ( 0.339 -2.016 -2.501) 3.230 1.118 ( 4.934 1.315 6.922) 8.602 1.230 ( 3.288 1.035 0.145) 3.450 1.406 ( -9.168 0.102 -0.876) 9.211 1.561 ( 4.536 -1.034 2.932) 5.500 1.818 ( 3.292 4.996 -4.366) 7.407 1.867 ( 7.546 3.327 -7.015) 10.827 1.960 ( 3.246 0.369 -5.408) 6.318 2.252 ( -1.783 -0.508 0.170) 1.862 2.328 ( -0.558 -1.837 3.373) 3.881 2.400 ( 0.889 -1.502 4.684) 4.999 2.513 ( 0.417 -1.795 -1.629) 2.459 2.591 ( -2.229 -1.206 -1.282) 2.840 2.639 ( -7.225 -1.253 2.313) 7.689 2.682 ( 0.751 0.800 0.974) 1.467 2.764 ( -2.721 2.320 -2.007) 4.100 2.890 ( -0.006 -2.394 7.365) 7.745 2.954 ( 1.225 -1.400 -0.605) 1.956 3.201 ( 0.897 0.113 1.675) 1.904 3.311 ( -2.316 2.483 -1.646) 3.774 3.323 ( -1.328 -3.458 -0.230) 3.711 3.523 ( -4.946 -1.139 0.996) 5.172 3.589 ( 1.250 1.330 -5.557) 5.849 3.731 ( 0.279 0.551 5.507) 5.542 3.769 ( -1.325 -0.498 2.124) 2.553 3.895 ( 3.509 -1.205 -3.050) 4.803 3.967 ( 0.344 2.160 -0.179) 2.194 4.078 ( 2.504 -0.009 3.546) 4.341 4.226 ( -1.704 3.208 -1.499) 3.929 4.512 ( 6.857 1.596 -5.306) 8.816 6.254 ( -0.856 -0.328 0.663) 1.131 6.312 ( -0.330 -0.211 -1.166) 1.231 7.153 ( 0.137 -0.059 3.615) 3.618 7.313 ( -1.277 -0.588 -2.419) 2.798 ======================= Grid point 71 (27/44) ======================= q-point: (-0.14 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.664 ( -9.501 -5.451 4.373) 11.795 0.714 ( -10.640 -1.695 6.315) 12.488 1.062 ( 6.959 1.218 -2.387) 7.457 1.201 ( -10.192 -0.895 6.503) 12.123 1.271 ( 3.326 -0.397 -0.974) 3.488 1.298 ( 1.914 0.881 1.441) 2.553 1.641 ( 7.770 -3.700 -3.525) 9.300 1.874 ( 1.286 7.452 -1.441) 7.698 1.937 ( 5.470 -0.261 -7.356) 9.171 2.050 ( 7.247 5.549 -1.810) 9.305 2.240 ( -0.364 0.774 1.359) 1.606 2.312 ( -2.982 -2.451 -1.108) 4.016 2.382 ( -2.807 0.697 2.086) 3.566 2.501 ( 2.470 -2.843 6.292) 7.333 2.561 ( -0.465 -0.576 0.773) 1.071 2.584 ( -1.342 -1.357 -1.273) 2.294 2.625 ( -0.969 -1.180 -8.257) 8.397 2.800 ( -0.631 3.305 4.738) 5.811 2.913 ( 0.071 -4.416 5.694) 7.206 2.962 ( 1.770 -0.713 -1.258) 2.285 3.206 ( -1.288 -1.154 -0.993) 1.994 3.279 ( 0.479 -2.201 -0.703) 2.360 3.324 ( -0.869 4.451 4.208) 6.187 3.432 ( -6.828 0.059 5.482) 8.756 3.603 ( 2.717 -0.024 -3.501) 4.432 3.729 ( -2.489 -0.380 2.673) 3.672 3.801 ( 2.112 -0.818 1.554) 2.746 3.876 ( -1.071 -2.129 -3.715) 4.413 3.996 ( 2.037 1.347 -1.538) 2.886 4.127 ( 1.210 -1.509 3.264) 3.794 4.200 ( -4.361 2.667 0.586) 5.145 4.531 ( 1.735 1.020 -13.017) 13.172 6.249 ( -0.500 -0.121 1.622) 1.701 6.294 ( -0.228 -0.084 -1.600) 1.618 7.213 ( 2.087 0.332 7.720) 8.004 7.269 ( -0.973 -0.103 -3.471) 3.606 ======================= Grid point 75 (28/44) ======================= q-point: ( 0.43 0.43 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.889 ( 0.000 3.213 9.052) 9.606 0.963 ( 0.000 4.509 12.328) 13.127 1.055 ( 0.000 -3.385 0.790) 3.476 1.199 ( 0.000 -8.972 0.752) 9.003 1.425 ( 0.000 5.843 3.499) 6.810 1.442 ( -0.000 -7.553 -6.811) 10.171 1.527 ( -0.000 -2.989 -2.521) 3.911 1.622 ( -0.000 1.343 -1.885) 2.315 1.851 ( -0.000 0.670 -10.749) 10.770 1.865 ( -0.000 -1.229 -5.465) 5.601 2.315 ( 0.000 -2.834 1.829) 3.373 2.359 ( 0.000 -2.686 2.213) 3.480 2.386 ( 0.000 -3.746 -0.072) 3.747 2.442 ( 0.000 -3.122 4.500) 5.477 2.606 ( 0.000 -0.114 0.036) 0.119 2.631 ( 0.000 0.293 5.173) 5.181 2.851 ( 0.000 -1.480 5.676) 5.866 2.860 ( -0.000 -0.724 -0.759) 1.049 2.955 ( -0.000 -0.333 -0.582) 0.671 3.002 ( -0.000 3.278 -1.716) 3.700 3.174 ( 0.000 2.208 1.437) 2.635 3.297 ( 0.000 -3.148 1.826) 3.639 3.508 ( -0.000 -0.310 -6.492) 6.500 3.518 ( 0.000 -4.507 -0.284) 4.516 3.604 ( 0.000 2.995 3.885) 4.905 3.793 ( -0.000 2.555 -2.753) 3.756 3.801 ( 0.000 0.523 6.625) 6.646 3.868 ( 0.000 2.343 3.705) 4.384 3.916 ( -0.000 1.402 -7.342) 7.475 3.923 ( -0.000 4.827 -0.687) 4.876 4.128 ( -0.000 4.071 -2.517) 4.786 4.235 ( -0.000 2.430 -0.371) 2.459 6.288 ( 0.000 -0.687 0.236) 0.726 6.330 ( 0.000 -0.473 0.383) 0.609 7.138 ( 0.000 -0.052 0.858) 0.860 7.393 ( -0.000 -2.315 -2.429) 3.356 ======================= Grid point 76 (29/44) ======================= q-point: ( 0.57 0.43 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.980 ( -0.468 -2.805 4.828) 5.604 0.993 ( 0.919 -3.364 4.446) 5.650 1.073 ( -2.446 -1.748 4.372) 5.306 1.136 ( -4.228 -1.382 5.111) 6.775 1.263 ( -5.266 -4.534 -3.505) 7.782 1.453 ( 1.180 -1.359 -7.338) 7.556 1.551 ( 5.349 3.403 -0.806) 6.390 1.641 ( -1.133 1.594 4.741) 5.128 1.791 ( 0.530 3.025 -8.697) 9.224 1.849 ( 3.359 -0.353 -5.286) 6.273 2.306 ( -0.146 -2.284 1.785) 2.903 2.335 ( -0.800 -2.274 0.927) 2.582 2.378 ( 0.906 -2.592 3.316) 4.306 2.456 ( 0.649 -3.079 5.195) 6.074 2.599 ( -0.621 0.211 -1.261) 1.421 2.661 ( 0.531 -0.718 1.791) 2.001 2.828 ( -1.786 -0.442 0.026) 1.840 2.890 ( 2.446 -2.463 -0.548) 3.514 2.912 ( -3.266 0.662 0.855) 3.440 2.996 ( -4.884 2.459 -0.444) 5.487 3.198 ( -0.559 0.533 2.722) 2.830 3.302 ( 2.525 -3.575 1.592) 4.657 3.410 ( -3.873 -0.174 -3.919) 5.513 3.488 ( 0.905 -0.703 -4.952) 5.083 3.666 ( 0.159 1.321 4.646) 4.832 3.765 ( -3.416 -0.377 1.008) 3.582 3.788 ( 0.044 0.675 -4.400) 4.452 3.877 ( -0.808 0.839 -3.926) 4.095 3.963 ( 4.139 1.071 2.443) 4.924 3.976 ( -2.394 4.483 6.740) 8.442 4.191 ( 3.694 4.688 -1.669) 6.197 4.285 ( 5.160 1.379 -2.469) 5.884 6.280 ( -0.555 -0.511 0.357) 0.835 6.321 ( -0.182 -0.445 -0.908) 1.028 7.152 ( -0.020 -0.257 3.010) 3.021 7.349 ( -0.711 -1.032 -3.184) 3.421 ======================= Grid point 77 (30/44) ======================= q-point: (-0.29 0.43 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.898 ( -4.506 -6.024 3.328) 8.227 0.990 ( -3.054 -4.923 -3.175) 6.606 1.007 ( -7.040 -0.817 3.294) 7.816 1.078 ( 4.276 -2.469 1.067) 5.052 1.243 ( -1.735 2.051 5.933) 6.513 1.411 ( 1.809 -1.261 -6.929) 7.271 1.605 ( -3.869 1.718 4.753) 6.365 1.711 ( 7.158 0.481 -4.169) 8.298 1.813 ( 2.321 6.736 -2.055) 7.415 1.875 ( 3.970 2.459 -4.348) 6.381 2.288 ( -1.987 -0.258 2.100) 2.903 2.327 ( 0.772 -1.876 0.727) 2.155 2.407 ( 1.348 -2.842 5.878) 6.667 2.482 ( 0.701 -1.432 2.177) 2.699 2.569 ( -0.806 -0.438 -2.937) 3.077 2.634 ( -0.849 -2.509 -4.897) 5.568 2.805 ( -2.920 0.954 2.758) 4.129 2.857 ( -6.158 -0.016 -1.287) 6.291 2.890 ( -0.596 -2.648 1.929) 3.330 2.916 ( -0.374 -1.059 -0.278) 1.157 3.215 ( -0.386 -0.061 2.825) 2.852 3.284 ( -0.099 -3.716 -1.856) 4.155 3.381 ( -1.560 3.294 1.927) 4.123 3.476 ( 3.243 1.228 -5.352) 6.377 3.668 ( -5.857 1.719 2.143) 6.469 3.730 ( 2.365 -3.968 -0.197) 4.623 3.783 ( 0.192 0.412 -0.101) 0.466 3.840 ( -0.767 0.401 -2.067) 2.241 4.020 ( -2.463 3.141 3.839) 5.538 4.042 ( 3.860 -1.320 1.704) 4.421 4.270 ( 0.762 4.006 1.383) 4.306 4.370 ( 6.841 1.327 -5.101) 8.636 6.270 ( -0.664 -0.340 0.833) 1.118 6.303 ( -0.268 -0.306 -1.689) 1.738 7.182 ( 0.294 -0.320 4.839) 4.858 7.303 ( -0.964 0.040 -4.010) 4.124 ======================= Grid point 78 (31/44) ======================= q-point: (-0.14 0.43 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.801 ( -6.814 -0.889 4.592) 8.265 0.920 ( -6.898 0.718 4.020) 8.016 0.939 ( 3.035 -1.310 -1.400) 3.590 1.118 ( 6.084 -4.402 -4.590) 8.801 1.286 ( -0.165 0.360 6.136) 6.149 1.372 ( 0.883 -2.230 -4.203) 4.839 1.575 ( -5.645 2.257 4.119) 7.344 1.761 ( 6.068 -2.993 -7.664) 10.223 1.924 ( 2.333 8.302 1.054) 8.687 1.978 ( 6.449 6.075 -1.993) 9.081 2.283 ( -1.351 1.114 2.723) 3.237 2.320 ( -0.336 -1.253 -0.603) 1.431 2.437 ( -2.819 0.561 -1.135) 3.090 2.473 ( 2.130 -2.936 3.515) 5.051 2.530 ( -0.700 -1.820 -1.395) 2.397 2.548 ( -0.713 -2.444 -5.007) 5.617 2.700 ( -7.742 -0.995 3.711) 8.642 2.833 ( 2.174 -5.809 0.852) 6.261 2.888 ( 0.280 0.383 2.625) 2.667 2.924 ( 1.376 0.077 -1.409) 1.971 3.218 ( -1.385 0.163 -0.109) 1.398 3.247 ( 0.853 -1.768 -0.203) 1.973 3.417 ( -1.218 4.259 4.648) 6.421 3.515 ( 3.840 0.334 -4.935) 6.262 3.592 ( -6.699 1.460 4.483) 8.191 3.729 ( 2.833 -4.638 -1.141) 5.553 3.795 ( -0.641 -0.997 3.528) 3.722 3.822 ( 0.563 0.807 -4.447) 4.555 4.014 ( -2.394 1.322 0.418) 2.767 4.082 ( 2.484 -3.214 -0.010) 4.062 4.312 ( -2.608 4.188 4.089) 6.408 4.420 ( 3.424 1.616 -8.222) 9.052 6.269 ( -0.344 -0.147 1.738) 1.777 6.284 ( -0.172 -0.193 -1.612) 1.633 7.236 ( 1.601 -0.282 6.048) 6.263 7.254 ( -1.051 0.453 -4.771) 4.906 ======================= Grid point 83 (32/44) ======================= q-point: (-0.43 -0.43 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.959 ( 0.000 -6.268 2.351) 6.695 0.981 ( 0.000 -3.607 1.742) 4.006 1.110 ( 0.000 -2.803 3.832) 4.748 1.152 ( -0.000 -6.190 -3.308) 7.019 1.348 ( 0.000 4.186 11.316) 12.065 1.373 ( -0.000 1.173 -8.562) 8.642 1.486 ( -0.000 0.671 -4.790) 4.837 1.716 ( 0.000 1.127 6.931) 7.022 1.770 ( -0.000 1.181 -5.620) 5.742 1.782 ( -0.000 5.835 -4.865) 7.598 2.288 ( 0.000 -2.048 0.564) 2.124 2.350 ( 0.000 -2.599 3.021) 3.986 2.364 ( 0.000 -0.882 2.507) 2.657 2.481 ( 0.000 -2.779 7.651) 8.139 2.599 ( -0.000 0.581 -2.053) 2.134 2.640 ( -0.000 -2.505 -1.195) 2.776 2.813 ( -0.000 -4.384 -3.092) 5.365 2.839 ( -0.000 0.009 -2.306) 2.306 2.947 ( 0.000 0.417 0.488) 0.642 3.060 ( -0.000 2.987 -1.308) 3.261 3.232 ( 0.000 -1.279 4.425) 4.606 3.254 ( 0.000 -2.773 0.785) 2.882 3.393 ( -0.000 1.698 -8.616) 8.782 3.466 ( 0.000 -0.273 0.979) 1.016 3.712 ( 0.000 0.607 4.659) 4.698 3.722 ( -0.000 -2.626 -6.740) 7.233 3.816 ( -0.000 2.380 -0.528) 2.438 3.895 ( -0.000 2.209 -1.360) 2.594 3.909 ( 0.000 -1.371 0.191) 1.384 4.142 ( 0.000 4.118 12.215) 12.890 4.167 ( -0.000 4.061 -2.764) 4.912 4.229 ( -0.000 0.278 -2.875) 2.889 6.282 ( 0.000 -0.481 0.414) 0.635 6.308 ( -0.000 -0.459 -1.715) 1.776 7.179 ( 0.000 -0.526 4.490) 4.521 7.320 ( -0.000 0.123 -4.286) 4.288 ======================= Grid point 84 (33/44) ======================= q-point: (-0.29 -0.43 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.909 ( -1.544 -5.890 -1.259) 6.218 0.948 ( 1.365 -5.988 -1.349) 6.288 0.989 ( -2.850 -5.308 -2.941) 6.705 1.140 ( -3.572 1.582 5.386) 6.653 1.331 ( 0.395 1.622 -7.824) 8.000 1.386 ( -4.534 2.012 4.777) 6.886 1.539 ( 5.816 -1.940 -0.295) 6.138 1.751 ( -3.698 3.338 4.963) 7.033 1.792 ( 4.806 4.339 -2.585) 6.972 1.853 ( 2.282 7.307 0.685) 7.686 2.283 ( 1.458 -1.461 -1.157) 2.367 2.336 ( -1.439 -1.263 2.572) 3.206 2.398 ( 1.501 -1.522 3.842) 4.396 2.510 ( -1.215 -1.701 6.449) 6.779 2.569 ( -1.597 -0.225 -2.084) 2.636 2.593 ( -0.991 -2.639 -3.219) 4.279 2.750 ( 2.003 -5.283 -4.577) 7.271 2.840 ( 2.294 0.695 -3.084) 3.906 2.929 ( -2.756 0.435 1.347) 3.099 3.010 ( -5.956 -0.196 -2.763) 6.569 3.229 ( -0.059 -2.904 1.120) 3.113 3.260 ( -0.487 -0.001 2.038) 2.095 3.372 ( 3.130 2.109 -5.705) 6.841 3.477 ( -1.112 2.001 2.960) 3.742 3.635 ( 0.279 -3.070 -6.599) 7.284 3.763 ( 1.182 0.237 4.404) 4.566 3.813 ( -3.663 0.163 -0.619) 3.718 3.835 ( -1.858 0.600 -1.617) 2.535 3.969 ( 3.989 -0.622 0.364) 4.053 4.135 ( -4.306 1.940 1.179) 4.868 4.247 ( 5.270 0.529 -2.718) 5.953 4.307 ( 2.434 3.063 7.948) 8.859 6.280 ( -0.258 -0.388 0.833) 0.954 6.288 ( -0.064 -0.407 -2.026) 2.067 7.213 ( 0.209 -0.669 5.089) 5.137 7.284 ( -0.597 1.130 -5.322) 5.473 ======================= Grid point 91 (34/44) ======================= q-point: (-0.29 -0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.821 ( -0.000 -9.303 -5.135) 10.626 0.822 ( -0.000 -10.985 -4.768) 11.975 1.017 ( 0.000 3.074 2.211) 3.786 1.230 ( 0.000 2.807 5.204) 5.913 1.295 ( -0.000 0.598 -8.266) 8.288 1.318 ( -0.000 -5.618 -5.847) 8.109 1.632 ( 0.000 4.415 12.160) 12.937 1.748 ( -0.000 6.145 -5.842) 8.479 1.875 ( 0.000 2.675 8.215) 8.639 1.918 ( 0.000 10.052 2.408) 10.336 2.249 ( -0.000 -1.048 -1.457) 1.795 2.341 ( 0.000 -4.443 3.364) 5.573 2.395 ( -0.000 0.641 -1.912) 2.016 2.563 ( -0.000 0.943 -2.557) 2.725 2.567 ( 0.000 -2.593 1.993) 3.270 2.577 ( 0.000 -1.884 7.719) 7.946 2.639 ( -0.000 -3.827 -6.758) 7.767 2.795 ( -0.000 1.284 -4.629) 4.803 2.966 ( 0.000 0.363 0.560) 0.668 3.012 ( -0.000 -1.490 -5.372) 5.575 3.224 ( 0.000 -1.779 1.714) 2.470 3.286 ( 0.000 -0.825 3.021) 3.131 3.307 ( -0.000 2.274 -6.974) 7.336 3.488 ( -0.000 -3.922 -1.625) 4.245 3.587 ( -0.000 1.696 -2.436) 2.968 3.760 ( 0.000 -2.291 0.576) 2.362 3.802 ( -0.000 1.047 -2.559) 2.765 3.895 ( 0.000 -1.484 2.151) 2.613 3.953 ( 0.000 3.164 2.235) 3.874 4.163 ( -0.000 0.434 -3.388) 3.416 4.182 ( -0.000 -0.745 -1.851) 1.995 4.404 ( 0.000 1.495 14.301) 14.379 6.268 ( -0.000 -0.474 -2.061) 2.115 6.286 ( 0.000 -0.330 1.107) 1.156 7.242 ( 0.000 -0.884 4.633) 4.717 7.257 ( -0.000 2.127 -6.769) 7.095 ======================= Grid point 116 (35/44) ======================= q-point: ( 0.29 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 44 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.707 ( 0.000 -0.000 15.789) 15.789 0.707 ( 0.000 -0.000 15.789) 15.789 1.146 ( -0.000 0.000 -4.850) 4.850 1.178 ( 0.000 -0.000 13.546) 13.546 1.488 ( 0.000 -0.000 2.537) 2.537 1.488 ( 0.000 -0.000 2.537) 2.537 1.703 ( -0.000 0.000 -8.103) 8.103 1.703 ( -0.000 0.000 -8.103) 8.103 1.938 ( -0.000 0.000 -12.996) 12.996 1.938 ( -0.000 0.000 -12.996) 12.996 2.388 ( 0.000 -0.000 4.177) 4.177 2.388 ( 0.000 -0.000 4.177) 4.177 2.454 ( 0.000 -0.000 2.218) 2.218 2.454 ( 0.000 -0.000 2.218) 2.218 2.674 ( 0.000 -0.000 3.823) 3.823 2.674 ( 0.000 -0.000 3.823) 3.823 2.778 ( -0.000 0.000 -1.157) 1.157 2.938 ( -0.000 0.000 -1.644) 1.644 2.946 ( 0.000 -0.000 1.171) 1.171 2.946 ( 0.000 -0.000 1.171) 1.171 3.079 ( 0.000 -0.000 3.311) 3.311 3.416 ( 0.000 -0.000 3.177) 3.177 3.416 ( 0.000 -0.000 3.177) 3.177 3.581 ( -0.000 0.000 -6.810) 6.810 3.672 ( 0.000 -0.000 8.919) 8.919 3.697 ( 0.000 -0.000 0.602) 0.602 3.697 ( 0.000 -0.000 0.602) 0.602 3.742 ( -0.000 0.000 -1.796) 1.796 3.742 ( -0.000 0.000 -1.796) 1.796 4.058 ( -0.000 0.000 -3.781) 3.781 4.058 ( -0.000 0.000 -3.781) 3.781 4.103 ( 0.000 -0.000 2.021) 2.021 6.310 ( -0.000 0.000 -0.173) 0.173 6.331 ( 0.000 -0.000 2.038) 2.038 7.138 ( -0.000 0.000 -1.085) 1.085 7.494 ( -0.000 0.000 -4.193) 4.193 ======================= Grid point 117 (36/44) ======================= q-point: ( 0.43 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.838 ( 1.526 0.881 14.528) 14.634 0.876 ( 3.055 1.764 14.566) 14.987 1.101 ( -2.196 -1.268 -0.427) 2.572 1.316 ( 3.962 2.288 6.260) 7.754 1.410 ( -5.589 -3.227 1.344) 6.592 1.540 ( -0.455 -0.263 4.331) 4.362 1.556 ( -5.470 -3.158 -5.563) 8.416 1.613 ( -2.149 -1.241 -11.462) 11.728 1.837 ( 0.921 0.532 -15.659) 15.695 1.885 ( 2.744 1.584 -7.670) 8.299 2.374 ( -2.068 -1.194 3.544) 4.273 2.429 ( -0.905 -0.522 3.039) 3.214 2.449 ( -2.313 -1.335 -0.372) 2.697 2.479 ( 1.374 0.793 2.809) 3.226 2.620 ( -3.494 -2.017 5.304) 6.665 2.739 ( -3.973 -2.294 3.917) 6.032 2.743 ( 0.599 0.346 3.987) 4.047 2.892 ( -2.062 -1.191 -2.580) 3.511 2.958 ( 1.309 0.756 -0.630) 1.637 2.982 ( 0.247 0.142 1.065) 1.102 3.114 ( 0.936 0.540 3.234) 3.410 3.380 ( -2.181 -1.259 2.180) 3.331 3.511 ( -1.489 -0.860 -6.177) 6.412 3.517 ( 3.493 2.017 3.702) 5.474 3.606 ( -3.211 -1.854 -3.900) 5.381 3.710 ( -1.784 -1.030 0.782) 2.203 3.755 ( 1.166 0.673 7.999) 8.112 3.781 ( 4.920 2.840 -0.910) 5.753 3.841 ( 5.352 3.090 1.157) 6.287 3.961 ( -2.577 -1.488 -4.447) 5.351 4.082 ( 2.763 1.595 -3.502) 4.737 4.177 ( 5.680 3.280 -0.796) 6.607 6.300 ( -0.728 -0.420 0.032) 0.841 6.338 ( -0.313 -0.181 0.602) 0.702 7.139 ( -0.028 -0.016 1.278) 1.278 7.436 ( -1.758 -1.015 -4.291) 4.747 ======================= Grid point 118 (37/44) ======================= q-point: ( 0.57 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.978 ( 0.802 0.463 9.953) 9.996 1.040 ( 2.538 1.465 9.957) 10.379 1.053 ( -2.718 -1.569 5.425) 6.267 1.233 ( -10.367 -5.986 -0.878) 12.003 1.297 ( -5.861 -3.384 -3.205) 7.488 1.494 ( -1.191 -0.688 -8.212) 8.327 1.556 ( -1.414 -0.816 6.271) 6.480 1.566 ( 7.257 4.190 0.479) 8.393 1.742 ( 1.959 1.131 -15.867) 16.027 1.892 ( 2.789 1.610 -6.631) 7.372 2.340 ( -3.414 -1.971 1.346) 4.166 2.396 ( -1.660 -0.958 1.101) 2.210 2.431 ( -1.053 -0.608 3.525) 3.729 2.521 ( 1.009 0.583 0.521) 1.277 2.609 ( -1.367 -0.789 4.245) 4.529 2.705 ( -2.531 -1.461 5.882) 6.568 2.760 ( -1.565 -0.904 1.488) 2.341 2.813 ( -3.036 -1.753 -3.840) 5.199 2.970 ( 0.313 0.181 -0.717) 0.803 2.971 ( -2.441 -1.410 3.983) 4.880 3.166 ( 1.459 0.842 3.080) 3.511 3.336 ( -2.639 -1.524 -0.785) 3.147 3.436 ( -2.273 -1.312 -3.017) 3.999 3.524 ( -0.235 -0.136 -4.907) 4.915 3.624 ( 3.756 2.168 3.646) 5.666 3.663 ( -3.407 -1.967 0.741) 4.003 3.819 ( -1.276 -0.737 6.577) 6.740 3.857 ( 3.998 2.308 -3.750) 5.948 3.901 ( 0.241 0.139 -3.495) 3.506 3.976 ( 4.545 2.624 5.603) 7.677 4.126 ( 3.528 2.037 -3.146) 5.148 4.307 ( 7.514 4.338 -4.078) 9.587 6.283 ( -0.988 -0.570 0.464) 1.231 6.328 ( -0.530 -0.306 -0.754) 0.971 7.156 ( -0.157 -0.090 3.443) 3.448 7.364 ( -1.741 -1.005 -4.158) 4.618 ======================= Grid point 119 (38/44) ======================= q-point: (-0.29 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.938 ( -6.348 -3.665 4.635) 8.673 1.002 ( -9.207 -5.315 -0.754) 10.657 1.073 ( -5.099 -2.944 1.297) 6.030 1.157 ( 3.215 1.856 9.648) 10.338 1.207 ( 1.154 0.666 -1.259) 1.833 1.415 ( -3.506 -2.024 -0.923) 4.153 1.610 ( 2.095 1.210 6.828) 7.244 1.664 ( 2.890 1.668 -13.355) 13.766 1.746 ( 8.725 5.038 -2.698) 10.430 1.892 ( 2.217 1.280 -7.194) 7.636 2.278 ( -2.630 -1.519 0.816) 3.145 2.383 ( -1.004 -0.580 2.418) 2.682 2.427 ( -1.503 -0.868 -0.290) 1.760 2.533 ( 0.282 0.163 0.368) 0.491 2.610 ( -0.862 -0.498 0.090) 1.000 2.687 ( -2.456 -1.418 -7.132) 7.675 2.704 ( -6.317 -3.647 3.811) 8.230 2.730 ( 0.135 0.078 6.845) 6.847 2.966 ( 0.034 0.020 -0.702) 0.704 2.979 ( -0.657 -0.379 6.293) 6.339 3.225 ( 1.398 0.807 4.614) 4.888 3.272 ( -2.098 -1.211 -1.422) 2.809 3.354 ( -2.442 -1.410 -4.243) 5.094 3.526 ( 2.651 1.530 -3.674) 4.782 3.582 ( -4.943 -2.854 3.114) 6.502 3.727 ( 3.290 1.900 2.734) 4.680 3.807 ( -2.604 -1.503 2.853) 4.145 3.904 ( 1.748 1.009 -1.233) 2.365 3.912 ( 2.349 1.356 -1.933) 3.330 4.098 ( 2.719 1.570 4.598) 5.568 4.168 ( 1.623 0.937 -0.703) 2.002 4.408 ( 6.027 3.480 -8.499) 10.985 6.269 ( -0.845 -0.488 1.127) 1.491 6.307 ( -0.553 -0.319 -1.472) 1.605 7.189 ( 0.167 0.096 5.389) 5.393 7.299 ( -1.206 -0.696 -4.380) 4.596 ======================= Grid point 124 (39/44) ======================= q-point: ( 0.43 0.43 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.947 ( 0.000 0.187 13.637) 13.638 1.000 ( 0.000 3.156 15.413) 15.733 1.119 ( 0.000 -1.789 3.673) 4.085 1.348 ( 0.000 4.164 1.634) 4.473 1.393 ( 0.000 -9.301 0.952) 9.350 1.515 ( -0.000 -2.422 -12.594) 12.825 1.520 ( -0.000 -5.140 -5.744) 7.708 1.600 ( 0.000 2.765 4.477) 5.262 1.726 ( -0.000 0.899 -13.510) 13.540 1.805 ( -0.000 1.861 -7.469) 7.698 2.390 ( 0.000 -3.318 2.413) 4.102 2.440 ( 0.000 -2.654 3.074) 4.061 2.454 ( 0.000 -2.866 -0.282) 2.880 2.494 ( 0.000 1.107 2.562) 2.791 2.677 ( 0.000 -5.743 3.605) 6.781 2.695 ( 0.000 -2.462 6.138) 6.614 2.872 ( -0.000 -1.881 -2.058) 2.788 2.875 ( 0.000 4.421 4.579) 6.365 2.949 ( -0.000 0.557 -1.075) 1.210 2.957 ( -0.000 -0.064 -2.721) 2.722 3.150 ( 0.000 1.405 4.061) 4.297 3.390 ( 0.000 -4.373 3.757) 5.766 3.443 ( -0.000 -0.538 -9.042) 9.057 3.570 ( 0.000 5.083 4.427) 6.740 3.593 ( -0.000 -3.482 -2.922) 4.545 3.706 ( -0.000 3.479 -4.486) 5.678 3.776 ( 0.000 2.845 0.844) 2.967 3.847 ( 0.000 2.816 8.744) 9.186 3.849 ( 0.000 3.325 5.049) 6.046 3.922 ( -0.000 -3.168 -4.515) 5.516 4.040 ( -0.000 2.575 -3.870) 4.648 4.152 ( -0.000 4.442 -3.031) 5.377 6.302 ( 0.000 -0.585 0.218) 0.624 6.338 ( -0.000 -0.437 -0.565) 0.714 7.159 ( 0.000 -0.173 3.727) 3.731 7.398 ( -0.000 -1.405 -5.474) 5.651 ======================= Grid point 125 (40/44) ======================= q-point: ( 0.57 0.43 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.050 ( 1.132 -0.958 9.984) 10.094 1.103 ( -2.053 -1.826 6.358) 6.926 1.143 ( -1.472 0.917 10.335) 10.479 1.276 ( -4.282 -9.334 0.407) 10.278 1.319 ( -5.676 -3.901 -4.567) 8.263 1.414 ( 1.337 -1.430 -9.555) 9.754 1.486 ( 5.772 1.602 -0.131) 5.991 1.629 ( -1.345 3.069 -9.116) 9.712 1.653 ( -1.058 3.260 2.327) 4.143 1.798 ( 3.390 0.652 -6.676) 7.515 2.372 ( -1.117 -3.053 1.669) 3.654 2.393 ( -1.456 -3.016 -0.010) 3.349 2.457 ( 0.760 -2.073 3.669) 4.282 2.514 ( -0.381 0.345 0.302) 0.596 2.643 ( -1.303 -4.417 3.995) 6.097 2.712 ( 0.041 -3.209 2.128) 3.851 2.829 ( -2.527 -0.385 -0.738) 2.661 2.870 ( -4.199 2.651 -0.467) 4.987 2.941 ( 1.031 -0.732 1.197) 1.741 2.981 ( 1.053 -0.665 -0.864) 1.516 3.200 ( 0.061 1.684 4.760) 5.049 3.351 ( -1.124 0.149 -4.683) 4.819 3.403 ( 1.185 -4.152 0.541) 4.351 3.493 ( -3.418 -2.750 -4.192) 6.068 3.655 ( 0.495 4.157 4.341) 6.031 3.721 ( -0.640 1.237 0.097) 1.396 3.791 ( 0.071 2.415 -2.600) 3.550 3.830 ( 1.215 -2.860 -3.777) 4.891 3.933 ( -1.537 4.286 4.869) 6.666 3.975 ( 3.170 1.743 5.588) 6.657 4.085 ( 3.777 4.854 -1.239) 6.274 4.200 ( 4.766 4.101 -5.465) 8.330 6.295 ( -0.563 -0.631 0.676) 1.083 6.322 ( -0.321 -0.583 -1.310) 1.469 7.192 ( -0.136 -0.381 5.280) 5.296 7.334 ( -0.752 -1.066 -5.634) 5.783 ======================= Grid point 126 (41/44) ======================= q-point: (-0.29 0.43 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.044 ( -5.576 -1.508 3.426) 6.716 1.070 ( -8.405 -6.492 -2.438) 10.897 1.105 ( -3.590 -9.242 2.677) 10.270 1.135 ( 2.961 -3.173 0.518) 4.371 1.255 ( -1.890 2.330 7.184) 7.785 1.386 ( 3.814 -1.787 -6.062) 7.381 1.569 ( -1.972 4.471 -4.975) 6.974 1.636 ( 7.001 3.391 0.340) 7.787 1.708 ( -0.077 4.609 4.931) 6.750 1.803 ( 3.048 0.630 -6.732) 7.417 2.318 ( -3.028 -1.173 1.128) 3.437 2.367 ( 0.731 -2.708 -1.037) 2.991 2.464 ( -2.859 0.735 -0.512) 2.997 2.510 ( 1.960 -1.783 1.059) 2.854 2.617 ( -1.231 -3.423 3.330) 4.932 2.644 ( -2.182 -2.945 -5.804) 6.865 2.809 ( -2.065 1.587 4.335) 5.057 2.826 ( -4.957 0.778 0.971) 5.111 2.938 ( 0.148 -1.609 -0.713) 1.766 2.982 ( 2.264 -3.134 1.595) 4.182 3.253 ( -0.267 -0.156 5.625) 5.633 3.302 ( -0.614 -0.160 -4.287) 4.333 3.352 ( -4.231 0.476 -0.907) 4.354 3.447 ( 3.832 -2.410 -1.940) 4.925 3.670 ( -4.505 -0.232 3.480) 5.697 3.745 ( 0.866 3.340 0.759) 3.533 3.800 ( 3.731 -2.334 -0.710) 4.458 3.820 ( 0.964 -0.785 -2.374) 2.679 3.968 ( -2.862 4.604 2.832) 6.116 4.071 ( 4.805 -0.255 1.271) 4.977 4.190 ( 1.074 5.371 2.824) 6.163 4.274 ( 5.704 3.657 -7.722) 10.273 6.287 ( -0.664 -0.528 1.328) 1.576 6.300 ( -0.431 -0.564 -1.398) 1.568 7.231 ( 0.084 -0.366 5.842) 5.855 7.270 ( -0.826 -0.506 -5.530) 5.615 ======================= Grid point 132 (42/44) ======================= q-point: ( 0.57 0.57 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.072 ( 0.000 -3.878 3.928) 5.520 1.144 ( 0.000 -6.181 5.768) 8.454 1.181 ( 0.000 -5.278 0.606) 5.313 1.270 ( 0.000 -6.207 0.991) 6.285 1.283 ( -0.000 0.157 -11.012) 11.013 1.397 ( 0.000 -0.481 9.857) 9.869 1.410 ( -0.000 3.590 -2.784) 4.543 1.617 ( -0.000 6.219 -6.812) 9.224 1.711 ( -0.000 -0.058 -6.461) 6.461 1.736 ( 0.000 2.559 5.992) 6.516 2.347 ( 0.000 -3.720 0.591) 3.767 2.404 ( 0.000 -1.806 0.735) 1.950 2.438 ( 0.000 -3.017 3.568) 4.673 2.532 ( -0.000 1.026 -0.474) 1.130 2.648 ( 0.000 -4.874 5.723) 7.517 2.687 ( -0.000 -3.450 -1.033) 3.601 2.822 ( -0.000 -0.692 -3.120) 3.196 2.868 ( -0.000 -3.357 -3.835) 5.097 2.944 ( 0.000 0.457 0.519) 0.691 2.999 ( 0.000 2.854 1.164) 3.083 3.257 ( 0.000 1.532 5.482) 5.692 3.301 ( -0.000 0.822 -8.181) 8.222 3.356 ( 0.000 -5.641 2.533) 6.183 3.493 ( -0.000 -2.208 -1.103) 2.468 3.694 ( -0.000 3.062 -0.915) 3.195 3.740 ( -0.000 -4.277 -5.221) 6.750 3.754 ( 0.000 5.008 2.656) 5.669 3.843 ( 0.000 2.444 0.916) 2.610 3.928 ( -0.000 -0.811 -0.570) 0.991 4.061 ( -0.000 5.084 -1.478) 5.294 4.124 ( 0.000 6.946 10.367) 12.479 4.156 ( -0.000 3.802 -5.548) 6.726 6.300 ( 0.000 -0.709 0.970) 1.202 6.307 ( -0.000 -0.791 -1.504) 1.699 7.233 ( 0.000 -0.769 5.985) 6.035 7.285 ( -0.000 -0.643 -6.062) 6.096 ======================= Grid point 174 (43/44) ======================= q-point: ( 0.43 0.43 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 44 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.084 ( 0.000 -0.000 15.867) 15.867 1.084 ( 0.000 -0.000 15.867) 15.867 1.177 ( 0.000 -0.000 6.181) 6.181 1.306 ( -0.000 0.000 -1.943) 1.943 1.440 ( -0.000 0.000 -13.258) 13.258 1.440 ( -0.000 0.000 -13.258) 13.258 1.567 ( 0.000 -0.000 3.076) 3.076 1.567 ( 0.000 -0.000 3.076) 3.076 1.685 ( -0.000 0.000 -7.311) 7.311 1.685 ( -0.000 0.000 -7.311) 7.311 2.463 ( 0.000 -0.000 2.130) 2.130 2.463 ( 0.000 -0.000 2.130) 2.130 2.486 ( 0.000 -0.000 0.231) 0.231 2.486 ( 0.000 -0.000 0.231) 0.231 2.793 ( 0.000 -0.000 5.679) 5.679 2.793 ( 0.000 -0.000 5.679) 5.679 2.828 ( 0.000 -0.000 3.744) 3.744 2.898 ( -0.000 0.000 -2.043) 2.043 2.912 ( -0.000 0.000 -3.763) 3.763 2.912 ( -0.000 0.000 -3.763) 3.763 3.189 ( 0.000 -0.000 6.562) 6.562 3.386 ( -0.000 0.000 -9.157) 9.157 3.503 ( 0.000 -0.000 4.046) 4.046 3.503 ( 0.000 -0.000 4.046) 4.046 3.606 ( -0.000 0.000 -4.662) 4.662 3.606 ( -0.000 0.000 -4.662) 4.662 3.824 ( 0.000 -0.000 5.186) 5.186 3.824 ( 0.000 -0.000 5.186) 5.186 3.889 ( 0.000 -0.000 9.076) 9.076 3.951 ( -0.000 0.000 -5.164) 5.164 3.951 ( -0.000 0.000 -5.164) 5.164 4.067 ( -0.000 0.000 -5.080) 5.080 6.311 ( 0.000 -0.000 0.505) 0.505 6.335 ( -0.000 0.000 -1.165) 1.165 7.200 ( 0.000 -0.000 5.780) 5.780 7.361 ( -0.000 0.000 -6.815) 6.815 ======================= Grid point 175 (44/44) ======================= q-point: ( 0.57 0.43 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 1.64e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 0.00e+00 0.00e+00 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.165 ( -3.250 -1.877 6.476) 7.485 1.186 ( -1.592 -0.919 7.309) 7.537 1.198 ( -1.086 -0.627 8.750) 8.839 1.318 ( -0.488 -0.282 -3.614) 3.657 1.350 ( 2.430 1.403 -6.882) 7.433 1.351 ( 0.291 0.168 -7.182) 7.190 1.513 ( -2.293 -1.324 2.307) 3.512 1.523 ( -4.333 -2.502 -6.495) 8.199 1.653 ( 1.216 0.702 5.438) 5.617 1.722 ( 3.573 2.063 -5.955) 7.244 2.428 ( -3.400 -1.963 -0.203) 3.931 2.439 ( -1.927 -1.112 0.616) 2.308 2.488 ( -0.579 -0.334 0.607) 0.903 2.507 ( 0.767 0.443 0.180) 0.904 2.748 ( -4.194 -2.421 3.775) 6.140 2.787 ( -4.250 -2.454 -3.989) 6.324 2.826 ( -1.859 -1.073 1.727) 2.755 2.861 ( -0.026 -0.015 3.649) 3.649 2.925 ( 2.097 1.211 -0.586) 2.492 2.963 ( 3.036 1.753 -2.253) 4.167 3.237 ( -0.413 -0.238 7.428) 7.443 3.299 ( -0.688 -0.397 -7.898) 7.938 3.480 ( -4.078 -2.354 -0.007) 4.709 3.510 ( -1.217 -0.702 -0.526) 1.500 3.604 ( 3.301 1.906 3.824) 5.399 3.672 ( 5.749 3.319 -4.381) 7.954 3.788 ( -3.888 -2.245 4.670) 6.478 3.850 ( -1.417 -0.818 -4.329) 4.628 3.918 ( 0.671 0.387 4.332) 4.401 3.971 ( 3.365 1.943 -0.756) 3.958 4.018 ( 4.631 2.674 5.233) 7.481 4.073 ( 4.751 2.743 -7.658) 9.420 6.311 ( -0.514 -0.297 1.015) 1.175 6.320 ( -0.432 -0.249 -1.279) 1.372 7.245 ( -0.385 -0.222 6.594) 6.609 7.298 ( -0.639 -0.369 -6.793) 6.833 =================== 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/12348 10.0 217.938 217.938 226.511 0.000 -0.000 -0.000 3/12348 20.0 44.399 44.399 50.886 0.000 -0.000 -0.000 3/12348 30.0 16.581 16.581 19.857 0.000 -0.000 -0.000 3/12348 40.0 9.046 9.046 11.229 0.000 -0.000 -0.000 3/12348 50.0 6.056 6.056 7.726 0.000 -0.000 -0.000 3/12348 60.0 4.539 4.539 5.910 0.000 -0.000 -0.000 3/12348 70.0 3.639 3.639 4.810 0.000 -0.000 -0.000 3/12348 80.0 3.045 3.045 4.072 0.000 -0.000 -0.000 3/12348 90.0 2.624 2.624 3.541 0.000 -0.000 -0.000 3/12348 100.0 2.309 2.309 3.138 0.000 -0.000 -0.000 3/12348 110.0 2.065 2.065 2.822 0.000 -0.000 -0.000 3/12348 120.0 1.869 1.869 2.567 0.000 -0.000 -0.000 3/12348 130.0 1.708 1.708 2.355 0.000 -0.000 -0.000 3/12348 140.0 1.574 1.574 2.177 0.000 -0.000 -0.000 3/12348 150.0 1.460 1.460 2.025 0.000 -0.000 -0.000 3/12348 160.0 1.362 1.362 1.893 0.000 -0.000 -0.000 3/12348 170.0 1.276 1.276 1.778 0.000 -0.000 -0.000 3/12348 180.0 1.201 1.201 1.676 0.000 -0.000 -0.000 3/12348 190.0 1.134 1.134 1.586 0.000 -0.000 -0.000 3/12348 200.0 1.075 1.075 1.505 0.000 -0.000 -0.000 3/12348 210.0 1.022 1.022 1.432 0.000 -0.000 -0.000 3/12348 220.0 0.973 0.973 1.366 0.000 -0.000 -0.000 3/12348 230.0 0.930 0.930 1.305 0.000 -0.000 -0.000 3/12348 240.0 0.890 0.890 1.250 0.000 -0.000 -0.000 3/12348 250.0 0.853 0.853 1.200 0.000 -0.000 -0.000 3/12348 260.0 0.819 0.819 1.153 0.000 -0.000 -0.000 3/12348 270.0 0.788 0.788 1.110 0.000 -0.000 -0.000 3/12348 280.0 0.759 0.759 1.070 0.000 -0.000 -0.000 3/12348 290.0 0.733 0.733 1.033 0.000 -0.000 -0.000 3/12348 300.0 0.708 0.708 0.998 0.000 -0.000 -0.000 3/12348 310.0 0.684 0.684 0.966 0.000 -0.000 -0.000 3/12348 320.0 0.663 0.663 0.936 0.000 -0.000 -0.000 3/12348 330.0 0.642 0.642 0.907 0.000 -0.000 -0.000 3/12348 340.0 0.623 0.623 0.880 0.000 -0.000 -0.000 3/12348 350.0 0.605 0.605 0.855 0.000 -0.000 -0.000 3/12348 360.0 0.588 0.588 0.831 0.000 -0.000 -0.000 3/12348 370.0 0.572 0.572 0.809 0.000 -0.000 -0.000 3/12348 380.0 0.557 0.557 0.787 0.000 -0.000 -0.000 3/12348 390.0 0.542 0.542 0.767 0.000 -0.000 -0.000 3/12348 400.0 0.529 0.529 0.748 0.000 -0.000 -0.000 3/12348 410.0 0.516 0.516 0.730 0.000 -0.000 -0.000 3/12348 420.0 0.503 0.503 0.712 0.000 -0.000 -0.000 3/12348 430.0 0.491 0.491 0.696 0.000 -0.000 -0.000 3/12348 440.0 0.480 0.480 0.680 0.000 -0.000 -0.000 3/12348 450.0 0.469 0.469 0.665 0.000 -0.000 -0.000 3/12348 460.0 0.459 0.459 0.650 0.000 -0.000 -0.000 3/12348 470.0 0.449 0.449 0.637 0.000 -0.000 -0.000 3/12348 480.0 0.440 0.440 0.623 0.000 -0.000 -0.000 3/12348 490.0 0.431 0.431 0.611 0.000 -0.000 -0.000 3/12348 500.0 0.422 0.422 0.598 0.000 -0.000 -0.000 3/12348 510.0 0.414 0.414 0.587 0.000 -0.000 -0.000 3/12348 520.0 0.406 0.406 0.575 0.000 -0.000 -0.000 3/12348 530.0 0.398 0.398 0.565 0.000 -0.000 -0.000 3/12348 540.0 0.391 0.391 0.554 0.000 -0.000 -0.000 3/12348 550.0 0.384 0.384 0.544 0.000 -0.000 -0.000 3/12348 560.0 0.377 0.377 0.534 0.000 -0.000 -0.000 3/12348 570.0 0.370 0.370 0.525 0.000 -0.000 -0.000 3/12348 580.0 0.364 0.364 0.516 0.000 -0.000 -0.000 3/12348 590.0 0.357 0.357 0.507 0.000 -0.000 -0.000 3/12348 600.0 0.351 0.351 0.499 0.000 -0.000 -0.000 3/12348 610.0 0.346 0.346 0.491 0.000 -0.000 -0.000 3/12348 620.0 0.340 0.340 0.483 0.000 -0.000 -0.000 3/12348 630.0 0.335 0.335 0.475 0.000 -0.000 -0.000 3/12348 640.0 0.329 0.329 0.468 0.000 -0.000 -0.000 3/12348 650.0 0.324 0.324 0.460 0.000 -0.000 -0.000 3/12348 660.0 0.319 0.319 0.453 0.000 -0.000 -0.000 3/12348 670.0 0.315 0.315 0.447 0.000 -0.000 -0.000 3/12348 680.0 0.310 0.310 0.440 0.000 -0.000 -0.000 3/12348 690.0 0.305 0.305 0.434 0.000 -0.000 -0.000 3/12348 700.0 0.301 0.301 0.428 0.000 -0.000 -0.000 3/12348 710.0 0.297 0.297 0.421 0.000 -0.000 -0.000 3/12348 720.0 0.293 0.293 0.416 0.000 -0.000 -0.000 3/12348 730.0 0.289 0.289 0.410 0.000 -0.000 -0.000 3/12348 740.0 0.285 0.285 0.404 0.000 -0.000 -0.000 3/12348 750.0 0.281 0.281 0.399 0.000 -0.000 -0.000 3/12348 760.0 0.277 0.277 0.394 0.000 -0.000 -0.000 3/12348 770.0 0.274 0.274 0.389 0.000 -0.000 -0.000 3/12348 780.0 0.270 0.270 0.384 0.000 -0.000 -0.000 3/12348 790.0 0.267 0.267 0.379 0.000 -0.000 -0.000 3/12348 800.0 0.263 0.263 0.374 0.000 -0.000 -0.000 3/12348 810.0 0.260 0.260 0.370 0.000 -0.000 -0.000 3/12348 820.0 0.257 0.257 0.365 0.000 -0.000 -0.000 3/12348 830.0 0.254 0.254 0.361 0.000 -0.000 -0.000 3/12348 840.0 0.251 0.251 0.356 0.000 -0.000 -0.000 3/12348 850.0 0.248 0.248 0.352 0.000 -0.000 -0.000 3/12348 860.0 0.245 0.245 0.348 0.000 -0.000 -0.000 3/12348 870.0 0.242 0.242 0.344 0.000 -0.000 -0.000 3/12348 880.0 0.239 0.239 0.340 0.000 -0.000 -0.000 3/12348 890.0 0.237 0.237 0.336 0.000 -0.000 -0.000 3/12348 900.0 0.234 0.234 0.333 0.000 -0.000 -0.000 3/12348 910.0 0.231 0.231 0.329 0.000 -0.000 -0.000 3/12348 920.0 0.229 0.229 0.325 0.000 -0.000 -0.000 3/12348 930.0 0.227 0.227 0.322 0.000 -0.000 -0.000 3/12348 940.0 0.224 0.224 0.318 0.000 -0.000 -0.000 3/12348 950.0 0.222 0.222 0.315 0.000 -0.000 -0.000 3/12348 960.0 0.219 0.219 0.312 0.000 -0.000 -0.000 3/12348 970.0 0.217 0.217 0.309 0.000 -0.000 -0.000 3/12348 980.0 0.215 0.215 0.305 0.000 -0.000 -0.000 3/12348 990.0 0.213 0.213 0.302 0.000 -0.000 -0.000 3/12348 1000.0 0.211 0.211 0.299 0.000 -0.000 -0.000 3/12348 Thermal conductivity related properties were written into "kappa-m777.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 01:37:18]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|