----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 17:51: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 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) Number of symmetry operations in supercell: 96 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 7.752067007667781 0.000000000000000 0.000000000000000 b -3.876033503833891 6.713486960479515 0.000000000000000 c 0.000000000000000 0.000000000000000 9.349081210000000 Atomic positions (fractional): *1 Hg 0.49797400511094 0.50202599488906 0.24249200353310 200.590 2 Hg 0.50202599488906 0.00405198977813 0.74249200353310 200.590 3 Hg 0.50202599488906 0.49797400511094 0.74249200353310 200.590 4 Hg 0.00405198977813 0.50202599488906 0.24249200353310 200.590 5 Hg 0.49797400511094 0.99594801022187 0.24249200353310 200.590 6 Hg 0.99594801022187 0.49797400511094 0.74249200353310 200.590 *7 As 0.00000000000000 0.00000000000000 0.22397030328982 74.922 8 As 0.00000000000000 0.00000000000000 0.72397030328982 74.922 *9 Se 0.84975093076138 0.69950186152275 0.36848393203731 78.960 *10 Se 0.33333333333333 0.66666666666667 0.09937833036002 78.960 11 Se 0.15024906923862 0.84975093076138 0.86848393203731 78.960 12 Se 0.15024906923862 0.30049813847725 0.86848393203731 78.960 13 Se 0.69950186152275 0.84975093076138 0.86848393203731 78.960 14 Se 0.84975093076138 0.15024906923862 0.36848393203731 78.960 15 Se 0.66666666666667 0.33333333333333 0.59937833036002 78.960 16 Se 0.30049813847725 0.15024906923862 0.36848393203731 78.960 *17 Br 0.33333333333333 0.66666666666667 0.52782356963893 79.904 18 Br 0.66666666666667 0.33333333333333 0.02782356963893 79.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 7.752067007667781 0.000000000000000 0.000000000000000 b -3.876033503833891 6.713486960479515 0.000000000000000 c 0.000000000000000 0.000000000000000 9.349081210000000 Atomic positions (fractional): *1 Hg 0.49797400511094 0.50202599488906 0.24249200353310 200.590 > 1 2 Hg 0.50202599488906 0.00405198977813 0.74249200353310 200.590 > 2 3 Hg 0.50202599488906 0.49797400511094 0.74249200353310 200.590 > 3 4 Hg 0.00405198977813 0.50202599488906 0.24249200353310 200.590 > 4 5 Hg 0.49797400511094 0.99594801022187 0.24249200353310 200.590 > 5 6 Hg 0.99594801022187 0.49797400511094 0.74249200353310 200.590 > 6 *7 As 0.00000000000000 0.00000000000000 0.22397030328982 74.922 > 7 8 As 0.00000000000000 0.00000000000000 0.72397030328982 74.922 > 8 *9 Se 0.84975093076138 0.69950186152275 0.36848393203731 78.960 > 9 *10 Se 0.33333333333333 0.66666666666667 0.09937833036002 78.960 > 10 11 Se 0.15024906923862 0.84975093076138 0.86848393203731 78.960 > 11 12 Se 0.15024906923862 0.30049813847725 0.86848393203731 78.960 > 12 13 Se 0.69950186152275 0.84975093076138 0.86848393203731 78.960 > 13 14 Se 0.84975093076138 0.15024906923862 0.36848393203731 78.960 > 14 15 Se 0.66666666666667 0.33333333333333 0.59937833036002 78.960 > 15 16 Se 0.30049813847725 0.15024906923862 0.36848393203731 78.960 > 16 *17 Br 0.33333333333333 0.66666666666667 0.52782356963893 79.904 > 17 18 Br 0.66666666666667 0.33333333333333 0.02782356963893 79.904 > 18 -------------------------------- super cell -------------------------------- Lattice vectors: a 15.504134015335563 0.000000000000000 0.000000000000000 b -7.752067007667782 13.426973920959030 0.000000000000000 c 0.000000000000000 0.000000000000000 18.698162419999999 Atomic positions (fractional): *1 Hg 0.24898700255547 0.25101299744453 0.12124600176655 200.590 > 1 2 Hg 0.74898700255547 0.25101299744453 0.12124600176655 200.590 > 1 3 Hg 0.24898700255547 0.75101299744453 0.12124600176655 200.590 > 1 4 Hg 0.74898700255547 0.75101299744453 0.12124600176655 200.590 > 1 5 Hg 0.24898700255547 0.25101299744453 0.62124600176655 200.590 > 1 6 Hg 0.74898700255547 0.25101299744453 0.62124600176655 200.590 > 1 7 Hg 0.24898700255547 0.75101299744453 0.62124600176655 200.590 > 1 8 Hg 0.74898700255547 0.75101299744453 0.62124600176655 200.590 > 1 9 Hg 0.25101299744453 0.00202599488906 0.37124600176655 200.590 > 2 10 Hg 0.75101299744453 0.00202599488906 0.37124600176655 200.590 > 2 11 Hg 0.25101299744453 0.50202599488906 0.37124600176655 200.590 > 2 12 Hg 0.75101299744453 0.50202599488906 0.37124600176655 200.590 > 2 13 Hg 0.25101299744453 0.00202599488906 0.87124600176655 200.590 > 2 14 Hg 0.75101299744453 0.00202599488906 0.87124600176655 200.590 > 2 15 Hg 0.25101299744453 0.50202599488906 0.87124600176655 200.590 > 2 16 Hg 0.75101299744453 0.50202599488906 0.87124600176655 200.590 > 2 17 Hg 0.25101299744453 0.24898700255547 0.37124600176655 200.590 > 3 18 Hg 0.75101299744453 0.24898700255547 0.37124600176655 200.590 > 3 19 Hg 0.25101299744453 0.74898700255547 0.37124600176655 200.590 > 3 20 Hg 0.75101299744453 0.74898700255547 0.37124600176655 200.590 > 3 21 Hg 0.25101299744453 0.24898700255547 0.87124600176655 200.590 > 3 22 Hg 0.75101299744453 0.24898700255547 0.87124600176655 200.590 > 3 23 Hg 0.25101299744453 0.74898700255547 0.87124600176655 200.590 > 3 24 Hg 0.75101299744453 0.74898700255547 0.87124600176655 200.590 > 3 25 Hg 0.00202599488906 0.25101299744453 0.12124600176655 200.590 > 4 26 Hg 0.50202599488906 0.25101299744453 0.12124600176655 200.590 > 4 27 Hg 0.00202599488906 0.75101299744453 0.12124600176655 200.590 > 4 28 Hg 0.50202599488906 0.75101299744453 0.12124600176655 200.590 > 4 29 Hg 0.00202599488906 0.25101299744453 0.62124600176655 200.590 > 4 30 Hg 0.50202599488906 0.25101299744453 0.62124600176655 200.590 > 4 31 Hg 0.00202599488906 0.75101299744453 0.62124600176655 200.590 > 4 32 Hg 0.50202599488906 0.75101299744453 0.62124600176655 200.590 > 4 33 Hg 0.24898700255547 0.49797400511094 0.12124600176655 200.590 > 5 34 Hg 0.74898700255547 0.49797400511094 0.12124600176655 200.590 > 5 35 Hg 0.24898700255547 0.99797400511094 0.12124600176655 200.590 > 5 36 Hg 0.74898700255547 0.99797400511094 0.12124600176655 200.590 > 5 37 Hg 0.24898700255547 0.49797400511094 0.62124600176655 200.590 > 5 38 Hg 0.74898700255547 0.49797400511094 0.62124600176655 200.590 > 5 39 Hg 0.24898700255547 0.99797400511094 0.62124600176655 200.590 > 5 40 Hg 0.74898700255547 0.99797400511094 0.62124600176655 200.590 > 5 41 Hg 0.49797400511094 0.24898700255547 0.37124600176655 200.590 > 6 42 Hg 0.99797400511094 0.24898700255547 0.37124600176655 200.590 > 6 43 Hg 0.49797400511094 0.74898700255547 0.37124600176655 200.590 > 6 44 Hg 0.99797400511094 0.74898700255547 0.37124600176655 200.590 > 6 45 Hg 0.49797400511094 0.24898700255547 0.87124600176655 200.590 > 6 46 Hg 0.99797400511094 0.24898700255547 0.87124600176655 200.590 > 6 47 Hg 0.49797400511094 0.74898700255547 0.87124600176655 200.590 > 6 48 Hg 0.99797400511094 0.74898700255547 0.87124600176655 200.590 > 6 *49 As 0.00000000000000 0.00000000000000 0.11198515164491 74.922 > 7 50 As 0.50000000000000 0.00000000000000 0.11198515164491 74.922 > 7 51 As 0.00000000000000 0.50000000000000 0.11198515164491 74.922 > 7 52 As 0.50000000000000 0.50000000000000 0.11198515164491 74.922 > 7 53 As 0.00000000000000 0.00000000000000 0.61198515164491 74.922 > 7 54 As 0.50000000000000 0.00000000000000 0.61198515164491 74.922 > 7 55 As 0.00000000000000 0.50000000000000 0.61198515164491 74.922 > 7 56 As 0.50000000000000 0.50000000000000 0.61198515164491 74.922 > 7 57 As 0.00000000000000 0.00000000000000 0.36198515164491 74.922 > 8 58 As 0.50000000000000 0.00000000000000 0.36198515164491 74.922 > 8 59 As 0.00000000000000 0.50000000000000 0.36198515164491 74.922 > 8 60 As 0.50000000000000 0.50000000000000 0.36198515164491 74.922 > 8 61 As 0.00000000000000 0.00000000000000 0.86198515164491 74.922 > 8 62 As 0.50000000000000 0.00000000000000 0.86198515164491 74.922 > 8 63 As 0.00000000000000 0.50000000000000 0.86198515164491 74.922 > 8 64 As 0.50000000000000 0.50000000000000 0.86198515164491 74.922 > 8 *65 Se 0.42487546538069 0.34975093076138 0.18424196601866 78.960 > 9 66 Se 0.92487546538069 0.34975093076138 0.18424196601866 78.960 > 9 67 Se 0.42487546538069 0.84975093076138 0.18424196601866 78.960 > 9 68 Se 0.92487546538069 0.84975093076138 0.18424196601866 78.960 > 9 69 Se 0.42487546538069 0.34975093076138 0.68424196601866 78.960 > 9 70 Se 0.92487546538069 0.34975093076138 0.68424196601866 78.960 > 9 71 Se 0.42487546538069 0.84975093076138 0.68424196601866 78.960 > 9 72 Se 0.92487546538069 0.84975093076138 0.68424196601866 78.960 > 9 *73 Se 0.16666666666667 0.33333333333333 0.04968916518001 78.960 > 10 74 Se 0.66666666666667 0.33333333333333 0.04968916518001 78.960 > 10 75 Se 0.16666666666667 0.83333333333333 0.04968916518001 78.960 > 10 76 Se 0.66666666666667 0.83333333333333 0.04968916518001 78.960 > 10 77 Se 0.16666666666667 0.33333333333333 0.54968916518001 78.960 > 10 78 Se 0.66666666666667 0.33333333333333 0.54968916518001 78.960 > 10 79 Se 0.16666666666667 0.83333333333333 0.54968916518001 78.960 > 10 80 Se 0.66666666666667 0.83333333333333 0.54968916518001 78.960 > 10 81 Se 0.07512453461931 0.42487546538069 0.43424196601866 78.960 > 11 82 Se 0.57512453461931 0.42487546538069 0.43424196601866 78.960 > 11 83 Se 0.07512453461931 0.92487546538069 0.43424196601866 78.960 > 11 84 Se 0.57512453461931 0.92487546538069 0.43424196601866 78.960 > 11 85 Se 0.07512453461931 0.42487546538069 0.93424196601866 78.960 > 11 86 Se 0.57512453461931 0.42487546538069 0.93424196601866 78.960 > 11 87 Se 0.07512453461931 0.92487546538069 0.93424196601866 78.960 > 11 88 Se 0.57512453461931 0.92487546538069 0.93424196601866 78.960 > 11 89 Se 0.07512453461931 0.15024906923862 0.43424196601866 78.960 > 12 90 Se 0.57512453461931 0.15024906923862 0.43424196601866 78.960 > 12 91 Se 0.07512453461931 0.65024906923862 0.43424196601866 78.960 > 12 92 Se 0.57512453461931 0.65024906923862 0.43424196601866 78.960 > 12 93 Se 0.07512453461931 0.15024906923862 0.93424196601866 78.960 > 12 94 Se 0.57512453461931 0.15024906923862 0.93424196601866 78.960 > 12 95 Se 0.07512453461931 0.65024906923862 0.93424196601866 78.960 > 12 96 Se 0.57512453461931 0.65024906923862 0.93424196601866 78.960 > 12 97 Se 0.34975093076138 0.42487546538069 0.43424196601866 78.960 > 13 98 Se 0.84975093076138 0.42487546538069 0.43424196601866 78.960 > 13 99 Se 0.34975093076138 0.92487546538069 0.43424196601866 78.960 > 13 100 Se 0.84975093076138 0.92487546538069 0.43424196601866 78.960 > 13 101 Se 0.34975093076138 0.42487546538069 0.93424196601866 78.960 > 13 102 Se 0.84975093076138 0.42487546538069 0.93424196601866 78.960 > 13 103 Se 0.34975093076138 0.92487546538069 0.93424196601866 78.960 > 13 104 Se 0.84975093076138 0.92487546538069 0.93424196601866 78.960 > 13 105 Se 0.42487546538069 0.07512453461931 0.18424196601866 78.960 > 14 106 Se 0.92487546538069 0.07512453461931 0.18424196601866 78.960 > 14 107 Se 0.42487546538069 0.57512453461931 0.18424196601866 78.960 > 14 108 Se 0.92487546538069 0.57512453461931 0.18424196601866 78.960 > 14 109 Se 0.42487546538069 0.07512453461931 0.68424196601866 78.960 > 14 110 Se 0.92487546538069 0.07512453461931 0.68424196601866 78.960 > 14 111 Se 0.42487546538069 0.57512453461931 0.68424196601866 78.960 > 14 112 Se 0.92487546538069 0.57512453461931 0.68424196601866 78.960 > 14 113 Se 0.33333333333333 0.16666666666667 0.29968916518001 78.960 > 15 114 Se 0.83333333333333 0.16666666666667 0.29968916518001 78.960 > 15 115 Se 0.33333333333333 0.66666666666667 0.29968916518001 78.960 > 15 116 Se 0.83333333333333 0.66666666666667 0.29968916518001 78.960 > 15 117 Se 0.33333333333333 0.16666666666667 0.79968916518001 78.960 > 15 118 Se 0.83333333333333 0.16666666666667 0.79968916518001 78.960 > 15 119 Se 0.33333333333333 0.66666666666667 0.79968916518001 78.960 > 15 120 Se 0.83333333333333 0.66666666666667 0.79968916518001 78.960 > 15 121 Se 0.15024906923862 0.07512453461931 0.18424196601866 78.960 > 16 122 Se 0.65024906923862 0.07512453461931 0.18424196601866 78.960 > 16 123 Se 0.15024906923862 0.57512453461931 0.18424196601866 78.960 > 16 124 Se 0.65024906923862 0.57512453461931 0.18424196601866 78.960 > 16 125 Se 0.15024906923862 0.07512453461931 0.68424196601866 78.960 > 16 126 Se 0.65024906923862 0.07512453461931 0.68424196601866 78.960 > 16 127 Se 0.15024906923862 0.57512453461931 0.68424196601866 78.960 > 16 128 Se 0.65024906923862 0.57512453461931 0.68424196601866 78.960 > 16 *129 Br 0.16666666666667 0.33333333333333 0.26391178481947 79.904 > 17 130 Br 0.66666666666667 0.33333333333333 0.26391178481947 79.904 > 17 131 Br 0.16666666666667 0.83333333333333 0.26391178481947 79.904 > 17 132 Br 0.66666666666667 0.83333333333333 0.26391178481947 79.904 > 17 133 Br 0.16666666666667 0.33333333333333 0.76391178481947 79.904 > 17 134 Br 0.66666666666667 0.33333333333333 0.76391178481947 79.904 > 17 135 Br 0.16666666666667 0.83333333333333 0.76391178481947 79.904 > 17 136 Br 0.66666666666667 0.83333333333333 0.76391178481947 79.904 > 17 137 Br 0.33333333333333 0.16666666666667 0.01391178481947 79.904 > 18 138 Br 0.83333333333333 0.16666666666667 0.01391178481947 79.904 > 18 139 Br 0.33333333333333 0.66666666666667 0.01391178481947 79.904 > 18 140 Br 0.83333333333333 0.66666666666667 0.01391178481947 79.904 > 18 141 Br 0.33333333333333 0.16666666666667 0.51391178481947 79.904 > 18 142 Br 0.83333333333333 0.16666666666667 0.51391178481947 79.904 > 18 143 Br 0.33333333333333 0.66666666666667 0.51391178481947 79.904 > 18 144 Br 0.83333333333333 0.66666666666667 0.51391178481947 79.904 > 18 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 13.0621474 -0.0000000 0.0000000 -0.0000000 13.0621474 0.0000000 0.0000000 0.0000000 7.6492589 -------------------------- Born effective charges -------------------------- 1 Hg 3.0607681 0.1557450 0.5855930 0.1557450 3.2406069 -0.3380923 0.0924071 -0.0533513 1.7726225 2 Hg 3.3305263 -0.0000000 0.0000000 0.0000000 2.9708487 -0.6761845 -0.0000000 -0.1067026 1.7726225 3 Hg 3.0607681 0.1557450 -0.5855930 0.1557450 3.2406069 0.3380923 -0.0924071 0.0533513 1.7726225 4 Hg 3.0607681 -0.1557450 -0.5855930 -0.1557450 3.2406069 -0.3380923 -0.0924071 -0.0533513 1.7726225 5 Hg 3.3305263 -0.0000000 -0.0000000 0.0000000 2.9708487 0.6761845 0.0000000 0.1067026 1.7726225 6 Hg 3.0607681 -0.1557450 0.5855930 -0.1557450 3.2406069 0.3380923 0.0924071 0.0533513 1.7726225 7 As 2.3742169 -0.0000000 0.0000000 0.0000000 2.3742169 0.0000000 0.0000000 0.0000000 0.2905210 8 As 2.3742169 -0.0000000 0.0000000 0.0000000 2.3742169 0.0000000 0.0000000 0.0000000 0.2905210 9 Se -2.0887676 0.0000000 0.0000000 0.0000000 -2.9930734 -0.0350361 -0.0000000 -0.3000816 -0.6352551 10 Se -2.2945737 0.0000000 0.0000000 -0.0000000 -2.2945737 0.0000000 0.0000000 0.0000000 -1.2220514 11 Se -2.7669970 0.3915759 0.0303422 0.3915759 -2.3148441 -0.0175181 0.2598783 -0.1500408 -0.6352551 12 Se -2.0887676 0.0000000 -0.0000000 0.0000000 -2.9930734 0.0350361 0.0000000 0.3000816 -0.6352551 13 Se -2.7669970 -0.3915759 -0.0303422 -0.3915759 -2.3148441 -0.0175181 -0.2598783 -0.1500408 -0.6352551 14 Se -2.7669970 0.3915759 -0.0303422 0.3915759 -2.3148441 0.0175181 -0.2598783 0.1500408 -0.6352551 15 Se -2.2945737 0.0000000 0.0000000 -0.0000000 -2.2945737 0.0000000 0.0000000 0.0000000 -1.2220514 16 Se -2.7669970 -0.3915759 0.0303422 -0.3915759 -2.3148441 0.0175181 0.2598783 0.1500408 -0.6352551 17 Br -1.9089441 0.0000000 0.0000000 -0.0000000 -1.9089441 0.0000000 0.0000000 0.0000000 -2.4805719 18 Br -1.9089441 0.0000000 0.0000000 -0.0000000 -1.9089441 0.0000000 0.0000000 0.0000000 -2.4805719 ---------------------------------------------------------------------------- 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: 12204/12204 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 1067 Number of blocks in projector: 724 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: 467 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 420 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 180 Use standard eigh solver. Tree of FC basis block matrices: - (1067, 1051), data: False |-- (180, 175), data: True |-- (420, 416), data: True |-- (467, 460), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.018 Solver_block: 100 / 240 - Time: 1.432 Solver_block: 200 / 240 - Time: 1.473 Solver_block: 240 / 240 - Time: 0.658 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 3.606 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 432/432 Permutation basis: 12204/12204 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 1067 Number of blocks in projector: 724 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: 467 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 420 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 180 Use standard eigh solver. Tree of FC basis block matrices: - (1067, 1051), data: False |-- (180, 175), data: True |-- (420, 416), data: True |-- (467, 460), data: True Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 17:51:51]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 17:51:51]------------------------- 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 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 7.752067007667781 0.000000000000000 0.000000000000000 b -3.876033503833891 6.713486960479515 0.000000000000000 c 0.000000000000000 0.000000000000000 9.349081210000000 Atomic positions (fractional): 1 Hg 0.49797400511094 0.50202599488906 0.24249200353310 200.590 2 Hg 0.50202599488906 0.00405198977813 0.74249200353310 200.590 3 Hg 0.50202599488906 0.49797400511094 0.74249200353310 200.590 4 Hg 0.00405198977813 0.50202599488906 0.24249200353310 200.590 5 Hg 0.49797400511094 0.99594801022187 0.24249200353310 200.590 6 Hg 0.99594801022187 0.49797400511094 0.74249200353310 200.590 7 As 0.00000000000000 0.00000000000000 0.22397030328982 74.922 8 As 0.00000000000000 0.00000000000000 0.72397030328982 74.922 9 Se 0.84975093076138 0.69950186152275 0.36848393203731 78.960 10 Se 0.33333333333333 0.66666666666667 0.09937833036002 78.960 11 Se 0.15024906923862 0.84975093076138 0.86848393203731 78.960 12 Se 0.15024906923862 0.30049813847725 0.86848393203731 78.960 13 Se 0.69950186152275 0.84975093076138 0.86848393203731 78.960 14 Se 0.84975093076138 0.15024906923862 0.36848393203731 78.960 15 Se 0.66666666666667 0.33333333333333 0.59937833036002 78.960 16 Se 0.30049813847725 0.15024906923862 0.36848393203731 78.960 17 Br 0.33333333333333 0.66666666666667 0.52782356963893 79.904 18 Br 0.66666666666667 0.33333333333333 0.02782356963893 79.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.504134015335563 0.000000000000000 0.000000000000000 b -7.752067007667782 13.426973920959030 0.000000000000000 c 0.000000000000000 0.000000000000000 18.698162419999999 Atomic positions (fractional): 1 Hg 0.24898700255547 0.25101299744453 0.12124600176655 200.590 > 1 2 Hg 0.74898700255547 0.25101299744453 0.12124600176655 200.590 > 1 3 Hg 0.24898700255547 0.75101299744453 0.12124600176655 200.590 > 1 4 Hg 0.74898700255547 0.75101299744453 0.12124600176655 200.590 > 1 5 Hg 0.24898700255547 0.25101299744453 0.62124600176655 200.590 > 1 6 Hg 0.74898700255547 0.25101299744453 0.62124600176655 200.590 > 1 7 Hg 0.24898700255547 0.75101299744453 0.62124600176655 200.590 > 1 8 Hg 0.74898700255547 0.75101299744453 0.62124600176655 200.590 > 1 9 Hg 0.25101299744453 0.00202599488906 0.37124600176655 200.590 > 9 10 Hg 0.75101299744453 0.00202599488906 0.37124600176655 200.590 > 9 11 Hg 0.25101299744453 0.50202599488906 0.37124600176655 200.590 > 9 12 Hg 0.75101299744453 0.50202599488906 0.37124600176655 200.590 > 9 13 Hg 0.25101299744453 0.00202599488906 0.87124600176655 200.590 > 9 14 Hg 0.75101299744453 0.00202599488906 0.87124600176655 200.590 > 9 15 Hg 0.25101299744453 0.50202599488906 0.87124600176655 200.590 > 9 16 Hg 0.75101299744453 0.50202599488906 0.87124600176655 200.590 > 9 17 Hg 0.25101299744453 0.24898700255547 0.37124600176655 200.590 > 17 18 Hg 0.75101299744453 0.24898700255547 0.37124600176655 200.590 > 17 19 Hg 0.25101299744453 0.74898700255547 0.37124600176655 200.590 > 17 20 Hg 0.75101299744453 0.74898700255547 0.37124600176655 200.590 > 17 21 Hg 0.25101299744453 0.24898700255547 0.87124600176655 200.590 > 17 22 Hg 0.75101299744453 0.24898700255547 0.87124600176655 200.590 > 17 23 Hg 0.25101299744453 0.74898700255547 0.87124600176655 200.590 > 17 24 Hg 0.75101299744453 0.74898700255547 0.87124600176655 200.590 > 17 25 Hg 0.00202599488906 0.25101299744453 0.12124600176655 200.590 > 25 26 Hg 0.50202599488906 0.25101299744453 0.12124600176655 200.590 > 25 27 Hg 0.00202599488906 0.75101299744453 0.12124600176655 200.590 > 25 28 Hg 0.50202599488906 0.75101299744453 0.12124600176655 200.590 > 25 29 Hg 0.00202599488906 0.25101299744453 0.62124600176655 200.590 > 25 30 Hg 0.50202599488906 0.25101299744453 0.62124600176655 200.590 > 25 31 Hg 0.00202599488906 0.75101299744453 0.62124600176655 200.590 > 25 32 Hg 0.50202599488906 0.75101299744453 0.62124600176655 200.590 > 25 33 Hg 0.24898700255547 0.49797400511094 0.12124600176655 200.590 > 33 34 Hg 0.74898700255547 0.49797400511094 0.12124600176655 200.590 > 33 35 Hg 0.24898700255547 0.99797400511094 0.12124600176655 200.590 > 33 36 Hg 0.74898700255547 0.99797400511094 0.12124600176655 200.590 > 33 37 Hg 0.24898700255547 0.49797400511094 0.62124600176655 200.590 > 33 38 Hg 0.74898700255547 0.49797400511094 0.62124600176655 200.590 > 33 39 Hg 0.24898700255547 0.99797400511094 0.62124600176655 200.590 > 33 40 Hg 0.74898700255547 0.99797400511094 0.62124600176655 200.590 > 33 41 Hg 0.49797400511094 0.24898700255547 0.37124600176655 200.590 > 41 42 Hg 0.99797400511094 0.24898700255547 0.37124600176655 200.590 > 41 43 Hg 0.49797400511094 0.74898700255547 0.37124600176655 200.590 > 41 44 Hg 0.99797400511094 0.74898700255547 0.37124600176655 200.590 > 41 45 Hg 0.49797400511094 0.24898700255547 0.87124600176655 200.590 > 41 46 Hg 0.99797400511094 0.24898700255547 0.87124600176655 200.590 > 41 47 Hg 0.49797400511094 0.74898700255547 0.87124600176655 200.590 > 41 48 Hg 0.99797400511094 0.74898700255547 0.87124600176655 200.590 > 41 49 As 0.00000000000000 0.00000000000000 0.11198515164491 74.922 > 49 50 As 0.50000000000000 0.00000000000000 0.11198515164491 74.922 > 49 51 As 0.00000000000000 0.50000000000000 0.11198515164491 74.922 > 49 52 As 0.50000000000000 0.50000000000000 0.11198515164491 74.922 > 49 53 As 0.00000000000000 0.00000000000000 0.61198515164491 74.922 > 49 54 As 0.50000000000000 0.00000000000000 0.61198515164491 74.922 > 49 55 As 0.00000000000000 0.50000000000000 0.61198515164491 74.922 > 49 56 As 0.50000000000000 0.50000000000000 0.61198515164491 74.922 > 49 57 As 0.00000000000000 0.00000000000000 0.36198515164491 74.922 > 57 58 As 0.50000000000000 0.00000000000000 0.36198515164491 74.922 > 57 59 As 0.00000000000000 0.50000000000000 0.36198515164491 74.922 > 57 60 As 0.50000000000000 0.50000000000000 0.36198515164491 74.922 > 57 61 As 0.00000000000000 0.00000000000000 0.86198515164491 74.922 > 57 62 As 0.50000000000000 0.00000000000000 0.86198515164491 74.922 > 57 63 As 0.00000000000000 0.50000000000000 0.86198515164491 74.922 > 57 64 As 0.50000000000000 0.50000000000000 0.86198515164491 74.922 > 57 65 Se 0.42487546538069 0.34975093076138 0.18424196601866 78.960 > 65 66 Se 0.92487546538069 0.34975093076138 0.18424196601866 78.960 > 65 67 Se 0.42487546538069 0.84975093076138 0.18424196601866 78.960 > 65 68 Se 0.92487546538069 0.84975093076138 0.18424196601866 78.960 > 65 69 Se 0.42487546538069 0.34975093076138 0.68424196601866 78.960 > 65 70 Se 0.92487546538069 0.34975093076138 0.68424196601866 78.960 > 65 71 Se 0.42487546538069 0.84975093076138 0.68424196601866 78.960 > 65 72 Se 0.92487546538069 0.84975093076138 0.68424196601866 78.960 > 65 73 Se 0.16666666666667 0.33333333333333 0.04968916518001 78.960 > 73 74 Se 0.66666666666667 0.33333333333333 0.04968916518001 78.960 > 73 75 Se 0.16666666666667 0.83333333333333 0.04968916518001 78.960 > 73 76 Se 0.66666666666667 0.83333333333333 0.04968916518001 78.960 > 73 77 Se 0.16666666666667 0.33333333333333 0.54968916518001 78.960 > 73 78 Se 0.66666666666667 0.33333333333333 0.54968916518001 78.960 > 73 79 Se 0.16666666666667 0.83333333333333 0.54968916518001 78.960 > 73 80 Se 0.66666666666667 0.83333333333333 0.54968916518001 78.960 > 73 81 Se 0.07512453461931 0.42487546538069 0.43424196601866 78.960 > 81 82 Se 0.57512453461931 0.42487546538069 0.43424196601866 78.960 > 81 83 Se 0.07512453461931 0.92487546538069 0.43424196601866 78.960 > 81 84 Se 0.57512453461931 0.92487546538069 0.43424196601866 78.960 > 81 85 Se 0.07512453461931 0.42487546538069 0.93424196601866 78.960 > 81 86 Se 0.57512453461931 0.42487546538069 0.93424196601866 78.960 > 81 87 Se 0.07512453461931 0.92487546538069 0.93424196601866 78.960 > 81 88 Se 0.57512453461931 0.92487546538069 0.93424196601866 78.960 > 81 89 Se 0.07512453461931 0.15024906923862 0.43424196601866 78.960 > 89 90 Se 0.57512453461931 0.15024906923862 0.43424196601866 78.960 > 89 91 Se 0.07512453461931 0.65024906923862 0.43424196601866 78.960 > 89 92 Se 0.57512453461931 0.65024906923862 0.43424196601866 78.960 > 89 93 Se 0.07512453461931 0.15024906923862 0.93424196601866 78.960 > 89 94 Se 0.57512453461931 0.15024906923862 0.93424196601866 78.960 > 89 95 Se 0.07512453461931 0.65024906923862 0.93424196601866 78.960 > 89 96 Se 0.57512453461931 0.65024906923862 0.93424196601866 78.960 > 89 97 Se 0.34975093076138 0.42487546538069 0.43424196601866 78.960 > 97 98 Se 0.84975093076138 0.42487546538069 0.43424196601866 78.960 > 97 99 Se 0.34975093076138 0.92487546538069 0.43424196601866 78.960 > 97 100 Se 0.84975093076138 0.92487546538069 0.43424196601866 78.960 > 97 101 Se 0.34975093076138 0.42487546538069 0.93424196601866 78.960 > 97 102 Se 0.84975093076138 0.42487546538069 0.93424196601866 78.960 > 97 103 Se 0.34975093076138 0.92487546538069 0.93424196601866 78.960 > 97 104 Se 0.84975093076138 0.92487546538069 0.93424196601866 78.960 > 97 105 Se 0.42487546538069 0.07512453461931 0.18424196601866 78.960 > 105 106 Se 0.92487546538069 0.07512453461931 0.18424196601866 78.960 > 105 107 Se 0.42487546538069 0.57512453461931 0.18424196601866 78.960 > 105 108 Se 0.92487546538069 0.57512453461931 0.18424196601866 78.960 > 105 109 Se 0.42487546538069 0.07512453461931 0.68424196601866 78.960 > 105 110 Se 0.92487546538069 0.07512453461931 0.68424196601866 78.960 > 105 111 Se 0.42487546538069 0.57512453461931 0.68424196601866 78.960 > 105 112 Se 0.92487546538069 0.57512453461931 0.68424196601866 78.960 > 105 113 Se 0.33333333333333 0.16666666666667 0.29968916518001 78.960 > 113 114 Se 0.83333333333333 0.16666666666667 0.29968916518001 78.960 > 113 115 Se 0.33333333333333 0.66666666666667 0.29968916518001 78.960 > 113 116 Se 0.83333333333333 0.66666666666667 0.29968916518001 78.960 > 113 117 Se 0.33333333333333 0.16666666666667 0.79968916518001 78.960 > 113 118 Se 0.83333333333333 0.16666666666667 0.79968916518001 78.960 > 113 119 Se 0.33333333333333 0.66666666666667 0.79968916518001 78.960 > 113 120 Se 0.83333333333333 0.66666666666667 0.79968916518001 78.960 > 113 121 Se 0.15024906923862 0.07512453461931 0.18424196601866 78.960 > 121 122 Se 0.65024906923862 0.07512453461931 0.18424196601866 78.960 > 121 123 Se 0.15024906923862 0.57512453461931 0.18424196601866 78.960 > 121 124 Se 0.65024906923862 0.57512453461931 0.18424196601866 78.960 > 121 125 Se 0.15024906923862 0.07512453461931 0.68424196601866 78.960 > 121 126 Se 0.65024906923862 0.07512453461931 0.68424196601866 78.960 > 121 127 Se 0.15024906923862 0.57512453461931 0.68424196601866 78.960 > 121 128 Se 0.65024906923862 0.57512453461931 0.68424196601866 78.960 > 121 129 Br 0.16666666666667 0.33333333333333 0.26391178481947 79.904 > 129 130 Br 0.66666666666667 0.33333333333333 0.26391178481947 79.904 > 129 131 Br 0.16666666666667 0.83333333333333 0.26391178481947 79.904 > 129 132 Br 0.66666666666667 0.83333333333333 0.26391178481947 79.904 > 129 133 Br 0.16666666666667 0.33333333333333 0.76391178481947 79.904 > 129 134 Br 0.66666666666667 0.33333333333333 0.76391178481947 79.904 > 129 135 Br 0.16666666666667 0.83333333333333 0.76391178481947 79.904 > 129 136 Br 0.66666666666667 0.83333333333333 0.76391178481947 79.904 > 129 137 Br 0.33333333333333 0.16666666666667 0.01391178481947 79.904 > 137 138 Br 0.83333333333333 0.16666666666667 0.01391178481947 79.904 > 137 139 Br 0.33333333333333 0.66666666666667 0.01391178481947 79.904 > 137 140 Br 0.83333333333333 0.66666666666667 0.01391178481947 79.904 > 137 141 Br 0.33333333333333 0.16666666666667 0.51391178481947 79.904 > 137 142 Br 0.83333333333333 0.16666666666667 0.51391178481947 79.904 > 137 143 Br 0.33333333333333 0.66666666666667 0.51391178481947 79.904 > 137 144 Br 0.83333333333333 0.66666666666667 0.51391178481947 79.904 > 137 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 13.0621474 -0.0000000 0.0000000 -0.0000000 13.0621474 0.0000000 0.0000000 0.0000000 7.6492589 -------------------------- Born effective charges -------------------------- 1 Hg 3.0607681 0.1557450 0.5855930 0.1557450 3.2406069 -0.3380923 0.0924071 -0.0533513 1.7726225 2 Hg 3.3305263 -0.0000000 0.0000000 0.0000000 2.9708487 -0.6761845 -0.0000000 -0.1067026 1.7726225 3 Hg 3.0607681 0.1557450 -0.5855930 0.1557450 3.2406069 0.3380923 -0.0924071 0.0533513 1.7726225 4 Hg 3.0607681 -0.1557450 -0.5855930 -0.1557450 3.2406069 -0.3380923 -0.0924071 -0.0533513 1.7726225 5 Hg 3.3305263 -0.0000000 -0.0000000 0.0000000 2.9708487 0.6761845 0.0000000 0.1067026 1.7726225 6 Hg 3.0607681 -0.1557450 0.5855930 -0.1557450 3.2406069 0.3380923 0.0924071 0.0533513 1.7726225 7 As 2.3742169 -0.0000000 0.0000000 0.0000000 2.3742169 0.0000000 0.0000000 0.0000000 0.2905210 8 As 2.3742169 -0.0000000 0.0000000 0.0000000 2.3742169 0.0000000 0.0000000 0.0000000 0.2905210 9 Se -2.0887676 0.0000000 0.0000000 0.0000000 -2.9930734 -0.0350361 -0.0000000 -0.3000816 -0.6352551 10 Se -2.2945737 0.0000000 0.0000000 -0.0000000 -2.2945737 0.0000000 0.0000000 0.0000000 -1.2220514 11 Se -2.7669970 0.3915759 0.0303422 0.3915759 -2.3148441 -0.0175181 0.2598783 -0.1500408 -0.6352551 12 Se -2.0887676 0.0000000 -0.0000000 0.0000000 -2.9930734 0.0350361 0.0000000 0.3000816 -0.6352551 13 Se -2.7669970 -0.3915759 -0.0303422 -0.3915759 -2.3148441 -0.0175181 -0.2598783 -0.1500408 -0.6352551 14 Se -2.7669970 0.3915759 -0.0303422 0.3915759 -2.3148441 0.0175181 -0.2598783 0.1500408 -0.6352551 15 Se -2.2945737 0.0000000 0.0000000 -0.0000000 -2.2945737 0.0000000 0.0000000 0.0000000 -1.2220514 16 Se -2.7669970 -0.3915759 0.0303422 -0.3915759 -2.3148441 0.0175181 0.2598783 0.1500408 -0.6352551 17 Br -1.9089441 0.0000000 0.0000000 -0.0000000 -1.9089441 0.0000000 0.0000000 0.0000000 -2.4805719 18 Br -1.9089441 0.0000000 0.0000000 -0.0000000 -1.9089441 0.0000000 0.0000000 0.0000000 -2.4805719 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 49, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, 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[ 129, 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.00000083 (yyz) 0.00000083 (yyz) 0.00000083 (yzy) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000085 (yz) -0.00000085 (yz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 17:52:10]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 17:52:10]------------------------- 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 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 7.752067007667781 0.000000000000000 0.000000000000000 b -3.876033503833891 6.713486960479515 0.000000000000000 c 0.000000000000000 0.000000000000000 9.349081210000000 Atomic positions (fractional): 1 Hg 0.49797400511094 0.50202599488906 0.24249200353310 200.590 2 Hg 0.50202599488906 0.00405198977813 0.74249200353310 200.590 3 Hg 0.50202599488906 0.49797400511094 0.74249200353310 200.590 4 Hg 0.00405198977813 0.50202599488906 0.24249200353310 200.590 5 Hg 0.49797400511094 0.99594801022187 0.24249200353310 200.590 6 Hg 0.99594801022187 0.49797400511094 0.74249200353310 200.590 7 As 0.00000000000000 0.00000000000000 0.22397030328982 74.922 8 As 0.00000000000000 0.00000000000000 0.72397030328982 74.922 9 Se 0.84975093076138 0.69950186152275 0.36848393203731 78.960 10 Se 0.33333333333333 0.66666666666667 0.09937833036002 78.960 11 Se 0.15024906923862 0.84975093076138 0.86848393203731 78.960 12 Se 0.15024906923862 0.30049813847725 0.86848393203731 78.960 13 Se 0.69950186152275 0.84975093076138 0.86848393203731 78.960 14 Se 0.84975093076138 0.15024906923862 0.36848393203731 78.960 15 Se 0.66666666666667 0.33333333333333 0.59937833036002 78.960 16 Se 0.30049813847725 0.15024906923862 0.36848393203731 78.960 17 Br 0.33333333333333 0.66666666666667 0.52782356963893 79.904 18 Br 0.66666666666667 0.33333333333333 0.02782356963893 79.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.504134015335563 0.000000000000000 0.000000000000000 b -7.752067007667782 13.426973920959030 0.000000000000000 c 0.000000000000000 0.000000000000000 18.698162419999999 Atomic positions (fractional): 1 Hg 0.24898700255547 0.25101299744453 0.12124600176655 200.590 > 1 2 Hg 0.74898700255547 0.25101299744453 0.12124600176655 200.590 > 1 3 Hg 0.24898700255547 0.75101299744453 0.12124600176655 200.590 > 1 4 Hg 0.74898700255547 0.75101299744453 0.12124600176655 200.590 > 1 5 Hg 0.24898700255547 0.25101299744453 0.62124600176655 200.590 > 1 6 Hg 0.74898700255547 0.25101299744453 0.62124600176655 200.590 > 1 7 Hg 0.24898700255547 0.75101299744453 0.62124600176655 200.590 > 1 8 Hg 0.74898700255547 0.75101299744453 0.62124600176655 200.590 > 1 9 Hg 0.25101299744453 0.00202599488906 0.37124600176655 200.590 > 9 10 Hg 0.75101299744453 0.00202599488906 0.37124600176655 200.590 > 9 11 Hg 0.25101299744453 0.50202599488906 0.37124600176655 200.590 > 9 12 Hg 0.75101299744453 0.50202599488906 0.37124600176655 200.590 > 9 13 Hg 0.25101299744453 0.00202599488906 0.87124600176655 200.590 > 9 14 Hg 0.75101299744453 0.00202599488906 0.87124600176655 200.590 > 9 15 Hg 0.25101299744453 0.50202599488906 0.87124600176655 200.590 > 9 16 Hg 0.75101299744453 0.50202599488906 0.87124600176655 200.590 > 9 17 Hg 0.25101299744453 0.24898700255547 0.37124600176655 200.590 > 17 18 Hg 0.75101299744453 0.24898700255547 0.37124600176655 200.590 > 17 19 Hg 0.25101299744453 0.74898700255547 0.37124600176655 200.590 > 17 20 Hg 0.75101299744453 0.74898700255547 0.37124600176655 200.590 > 17 21 Hg 0.25101299744453 0.24898700255547 0.87124600176655 200.590 > 17 22 Hg 0.75101299744453 0.24898700255547 0.87124600176655 200.590 > 17 23 Hg 0.25101299744453 0.74898700255547 0.87124600176655 200.590 > 17 24 Hg 0.75101299744453 0.74898700255547 0.87124600176655 200.590 > 17 25 Hg 0.00202599488906 0.25101299744453 0.12124600176655 200.590 > 25 26 Hg 0.50202599488906 0.25101299744453 0.12124600176655 200.590 > 25 27 Hg 0.00202599488906 0.75101299744453 0.12124600176655 200.590 > 25 28 Hg 0.50202599488906 0.75101299744453 0.12124600176655 200.590 > 25 29 Hg 0.00202599488906 0.25101299744453 0.62124600176655 200.590 > 25 30 Hg 0.50202599488906 0.25101299744453 0.62124600176655 200.590 > 25 31 Hg 0.00202599488906 0.75101299744453 0.62124600176655 200.590 > 25 32 Hg 0.50202599488906 0.75101299744453 0.62124600176655 200.590 > 25 33 Hg 0.24898700255547 0.49797400511094 0.12124600176655 200.590 > 33 34 Hg 0.74898700255547 0.49797400511094 0.12124600176655 200.590 > 33 35 Hg 0.24898700255547 0.99797400511094 0.12124600176655 200.590 > 33 36 Hg 0.74898700255547 0.99797400511094 0.12124600176655 200.590 > 33 37 Hg 0.24898700255547 0.49797400511094 0.62124600176655 200.590 > 33 38 Hg 0.74898700255547 0.49797400511094 0.62124600176655 200.590 > 33 39 Hg 0.24898700255547 0.99797400511094 0.62124600176655 200.590 > 33 40 Hg 0.74898700255547 0.99797400511094 0.62124600176655 200.590 > 33 41 Hg 0.49797400511094 0.24898700255547 0.37124600176655 200.590 > 41 42 Hg 0.99797400511094 0.24898700255547 0.37124600176655 200.590 > 41 43 Hg 0.49797400511094 0.74898700255547 0.37124600176655 200.590 > 41 44 Hg 0.99797400511094 0.74898700255547 0.37124600176655 200.590 > 41 45 Hg 0.49797400511094 0.24898700255547 0.87124600176655 200.590 > 41 46 Hg 0.99797400511094 0.24898700255547 0.87124600176655 200.590 > 41 47 Hg 0.49797400511094 0.74898700255547 0.87124600176655 200.590 > 41 48 Hg 0.99797400511094 0.74898700255547 0.87124600176655 200.590 > 41 49 As 0.00000000000000 0.00000000000000 0.11198515164491 74.922 > 49 50 As 0.50000000000000 0.00000000000000 0.11198515164491 74.922 > 49 51 As 0.00000000000000 0.50000000000000 0.11198515164491 74.922 > 49 52 As 0.50000000000000 0.50000000000000 0.11198515164491 74.922 > 49 53 As 0.00000000000000 0.00000000000000 0.61198515164491 74.922 > 49 54 As 0.50000000000000 0.00000000000000 0.61198515164491 74.922 > 49 55 As 0.00000000000000 0.50000000000000 0.61198515164491 74.922 > 49 56 As 0.50000000000000 0.50000000000000 0.61198515164491 74.922 > 49 57 As 0.00000000000000 0.00000000000000 0.36198515164491 74.922 > 57 58 As 0.50000000000000 0.00000000000000 0.36198515164491 74.922 > 57 59 As 0.00000000000000 0.50000000000000 0.36198515164491 74.922 > 57 60 As 0.50000000000000 0.50000000000000 0.36198515164491 74.922 > 57 61 As 0.00000000000000 0.00000000000000 0.86198515164491 74.922 > 57 62 As 0.50000000000000 0.00000000000000 0.86198515164491 74.922 > 57 63 As 0.00000000000000 0.50000000000000 0.86198515164491 74.922 > 57 64 As 0.50000000000000 0.50000000000000 0.86198515164491 74.922 > 57 65 Se 0.42487546538069 0.34975093076138 0.18424196601866 78.960 > 65 66 Se 0.92487546538069 0.34975093076138 0.18424196601866 78.960 > 65 67 Se 0.42487546538069 0.84975093076138 0.18424196601866 78.960 > 65 68 Se 0.92487546538069 0.84975093076138 0.18424196601866 78.960 > 65 69 Se 0.42487546538069 0.34975093076138 0.68424196601866 78.960 > 65 70 Se 0.92487546538069 0.34975093076138 0.68424196601866 78.960 > 65 71 Se 0.42487546538069 0.84975093076138 0.68424196601866 78.960 > 65 72 Se 0.92487546538069 0.84975093076138 0.68424196601866 78.960 > 65 73 Se 0.16666666666667 0.33333333333333 0.04968916518001 78.960 > 73 74 Se 0.66666666666667 0.33333333333333 0.04968916518001 78.960 > 73 75 Se 0.16666666666667 0.83333333333333 0.04968916518001 78.960 > 73 76 Se 0.66666666666667 0.83333333333333 0.04968916518001 78.960 > 73 77 Se 0.16666666666667 0.33333333333333 0.54968916518001 78.960 > 73 78 Se 0.66666666666667 0.33333333333333 0.54968916518001 78.960 > 73 79 Se 0.16666666666667 0.83333333333333 0.54968916518001 78.960 > 73 80 Se 0.66666666666667 0.83333333333333 0.54968916518001 78.960 > 73 81 Se 0.07512453461931 0.42487546538069 0.43424196601866 78.960 > 81 82 Se 0.57512453461931 0.42487546538069 0.43424196601866 78.960 > 81 83 Se 0.07512453461931 0.92487546538069 0.43424196601866 78.960 > 81 84 Se 0.57512453461931 0.92487546538069 0.43424196601866 78.960 > 81 85 Se 0.07512453461931 0.42487546538069 0.93424196601866 78.960 > 81 86 Se 0.57512453461931 0.42487546538069 0.93424196601866 78.960 > 81 87 Se 0.07512453461931 0.92487546538069 0.93424196601866 78.960 > 81 88 Se 0.57512453461931 0.92487546538069 0.93424196601866 78.960 > 81 89 Se 0.07512453461931 0.15024906923862 0.43424196601866 78.960 > 89 90 Se 0.57512453461931 0.15024906923862 0.43424196601866 78.960 > 89 91 Se 0.07512453461931 0.65024906923862 0.43424196601866 78.960 > 89 92 Se 0.57512453461931 0.65024906923862 0.43424196601866 78.960 > 89 93 Se 0.07512453461931 0.15024906923862 0.93424196601866 78.960 > 89 94 Se 0.57512453461931 0.15024906923862 0.93424196601866 78.960 > 89 95 Se 0.07512453461931 0.65024906923862 0.93424196601866 78.960 > 89 96 Se 0.57512453461931 0.65024906923862 0.93424196601866 78.960 > 89 97 Se 0.34975093076138 0.42487546538069 0.43424196601866 78.960 > 97 98 Se 0.84975093076138 0.42487546538069 0.43424196601866 78.960 > 97 99 Se 0.34975093076138 0.92487546538069 0.43424196601866 78.960 > 97 100 Se 0.84975093076138 0.92487546538069 0.43424196601866 78.960 > 97 101 Se 0.34975093076138 0.42487546538069 0.93424196601866 78.960 > 97 102 Se 0.84975093076138 0.42487546538069 0.93424196601866 78.960 > 97 103 Se 0.34975093076138 0.92487546538069 0.93424196601866 78.960 > 97 104 Se 0.84975093076138 0.92487546538069 0.93424196601866 78.960 > 97 105 Se 0.42487546538069 0.07512453461931 0.18424196601866 78.960 > 105 106 Se 0.92487546538069 0.07512453461931 0.18424196601866 78.960 > 105 107 Se 0.42487546538069 0.57512453461931 0.18424196601866 78.960 > 105 108 Se 0.92487546538069 0.57512453461931 0.18424196601866 78.960 > 105 109 Se 0.42487546538069 0.07512453461931 0.68424196601866 78.960 > 105 110 Se 0.92487546538069 0.07512453461931 0.68424196601866 78.960 > 105 111 Se 0.42487546538069 0.57512453461931 0.68424196601866 78.960 > 105 112 Se 0.92487546538069 0.57512453461931 0.68424196601866 78.960 > 105 113 Se 0.33333333333333 0.16666666666667 0.29968916518001 78.960 > 113 114 Se 0.83333333333333 0.16666666666667 0.29968916518001 78.960 > 113 115 Se 0.33333333333333 0.66666666666667 0.29968916518001 78.960 > 113 116 Se 0.83333333333333 0.66666666666667 0.29968916518001 78.960 > 113 117 Se 0.33333333333333 0.16666666666667 0.79968916518001 78.960 > 113 118 Se 0.83333333333333 0.16666666666667 0.79968916518001 78.960 > 113 119 Se 0.33333333333333 0.66666666666667 0.79968916518001 78.960 > 113 120 Se 0.83333333333333 0.66666666666667 0.79968916518001 78.960 > 113 121 Se 0.15024906923862 0.07512453461931 0.18424196601866 78.960 > 121 122 Se 0.65024906923862 0.07512453461931 0.18424196601866 78.960 > 121 123 Se 0.15024906923862 0.57512453461931 0.18424196601866 78.960 > 121 124 Se 0.65024906923862 0.57512453461931 0.18424196601866 78.960 > 121 125 Se 0.15024906923862 0.07512453461931 0.68424196601866 78.960 > 121 126 Se 0.65024906923862 0.07512453461931 0.68424196601866 78.960 > 121 127 Se 0.15024906923862 0.57512453461931 0.68424196601866 78.960 > 121 128 Se 0.65024906923862 0.57512453461931 0.68424196601866 78.960 > 121 129 Br 0.16666666666667 0.33333333333333 0.26391178481947 79.904 > 129 130 Br 0.66666666666667 0.33333333333333 0.26391178481947 79.904 > 129 131 Br 0.16666666666667 0.83333333333333 0.26391178481947 79.904 > 129 132 Br 0.66666666666667 0.83333333333333 0.26391178481947 79.904 > 129 133 Br 0.16666666666667 0.33333333333333 0.76391178481947 79.904 > 129 134 Br 0.66666666666667 0.33333333333333 0.76391178481947 79.904 > 129 135 Br 0.16666666666667 0.83333333333333 0.76391178481947 79.904 > 129 136 Br 0.66666666666667 0.83333333333333 0.76391178481947 79.904 > 129 137 Br 0.33333333333333 0.16666666666667 0.01391178481947 79.904 > 137 138 Br 0.83333333333333 0.16666666666667 0.01391178481947 79.904 > 137 139 Br 0.33333333333333 0.66666666666667 0.01391178481947 79.904 > 137 140 Br 0.83333333333333 0.66666666666667 0.01391178481947 79.904 > 137 141 Br 0.33333333333333 0.16666666666667 0.51391178481947 79.904 > 137 142 Br 0.83333333333333 0.16666666666667 0.51391178481947 79.904 > 137 143 Br 0.33333333333333 0.66666666666667 0.51391178481947 79.904 > 137 144 Br 0.83333333333333 0.66666666666667 0.51391178481947 79.904 > 137 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 13.0621474 -0.0000000 0.0000000 -0.0000000 13.0621474 0.0000000 0.0000000 0.0000000 7.6492589 -------------------------- Born effective charges -------------------------- 1 Hg 3.0607681 0.1557450 0.5855930 0.1557450 3.2406069 -0.3380923 0.0924071 -0.0533513 1.7726225 2 Hg 3.3305263 -0.0000000 0.0000000 0.0000000 2.9708487 -0.6761845 -0.0000000 -0.1067026 1.7726225 3 Hg 3.0607681 0.1557450 -0.5855930 0.1557450 3.2406069 0.3380923 -0.0924071 0.0533513 1.7726225 4 Hg 3.0607681 -0.1557450 -0.5855930 -0.1557450 3.2406069 -0.3380923 -0.0924071 -0.0533513 1.7726225 5 Hg 3.3305263 -0.0000000 -0.0000000 0.0000000 2.9708487 0.6761845 0.0000000 0.1067026 1.7726225 6 Hg 3.0607681 -0.1557450 0.5855930 -0.1557450 3.2406069 0.3380923 0.0924071 0.0533513 1.7726225 7 As 2.3742169 -0.0000000 0.0000000 0.0000000 2.3742169 0.0000000 0.0000000 0.0000000 0.2905210 8 As 2.3742169 -0.0000000 0.0000000 0.0000000 2.3742169 0.0000000 0.0000000 0.0000000 0.2905210 9 Se -2.0887676 0.0000000 0.0000000 0.0000000 -2.9930734 -0.0350361 -0.0000000 -0.3000816 -0.6352551 10 Se -2.2945737 0.0000000 0.0000000 -0.0000000 -2.2945737 0.0000000 0.0000000 0.0000000 -1.2220514 11 Se -2.7669970 0.3915759 0.0303422 0.3915759 -2.3148441 -0.0175181 0.2598783 -0.1500408 -0.6352551 12 Se -2.0887676 0.0000000 -0.0000000 0.0000000 -2.9930734 0.0350361 0.0000000 0.3000816 -0.6352551 13 Se -2.7669970 -0.3915759 -0.0303422 -0.3915759 -2.3148441 -0.0175181 -0.2598783 -0.1500408 -0.6352551 14 Se -2.7669970 0.3915759 -0.0303422 0.3915759 -2.3148441 0.0175181 -0.2598783 0.1500408 -0.6352551 15 Se -2.2945737 0.0000000 0.0000000 -0.0000000 -2.2945737 0.0000000 0.0000000 0.0000000 -1.2220514 16 Se -2.7669970 -0.3915759 0.0303422 -0.3915759 -2.3148441 0.0175181 0.2598783 0.1500408 -0.6352551 17 Br -1.9089441 0.0000000 0.0000000 -0.0000000 -1.9089441 0.0000000 0.0000000 0.0000000 -2.4805719 18 Br -1.9089441 0.0000000 0.0000000 -0.0000000 -1.9089441 0.0000000 0.0000000 0.0000000 -2.4805719 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000083 (yyz) 0.00000083 (yyz) 0.00000083 (yzy) 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 5 ] 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.53, Number of G-points: 303, Lambda: 0.18 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/24) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 0.478 ( 0.000 -0.000 -0.000) 0.000 0.478 ( 0.000 -0.000 0.000) 0.000 0.639 ( 0.000 0.000 -0.000) 0.000 0.639 ( 0.000 0.000 -0.000) 0.000 0.775 ( 0.000 0.000 -0.000) 0.000 0.828 ( 0.000 0.000 -0.000) 0.000 0.990 ( 0.000 -0.000 -0.000) 0.000 0.990 ( 0.000 0.000 -0.000) 0.000 1.186 ( -0.000 -0.000 0.000) 0.000 1.432 ( 0.000 -0.000 0.000) 0.000 1.432 ( -0.000 0.000 0.000) 0.000 1.435 ( 0.000 0.000 0.000) 0.000 1.435 ( 0.000 -0.000 0.000) 0.000 1.682 ( 0.000 -0.000 0.000) 0.000 1.707 ( 0.000 -0.000 0.000) 0.000 1.707 ( 0.000 -0.000 0.000) 0.000 1.731 ( 0.000 0.000 0.000) 0.000 1.731 ( 0.000 0.000 0.000) 0.000 1.944 ( -0.000 -0.000 0.000) 0.000 2.293 ( 0.000 0.000 0.000) 0.000 2.293 ( 0.000 -0.000 0.000) 0.000 2.322 ( 0.000 0.000 0.000) 0.000 2.322 ( -0.000 -0.000 -0.000) 0.000 2.375 ( -0.000 0.000 0.000) 0.000 2.376 ( 0.000 0.000 -0.000) 0.000 2.737 ( 0.000 0.000 -0.000) 0.000 3.136 ( 0.000 0.000 0.000) 0.000 4.155 ( 0.000 -0.000 -0.000) 0.000 4.155 ( -0.000 -0.000 0.000) 0.000 4.180 ( -0.000 0.000 -0.000) 0.000 4.180 ( 0.000 -0.000 0.000) 0.000 4.308 ( 0.000 0.000 0.000) 0.000 4.329 ( 0.000 -0.000 -0.000) 0.000 4.764 ( 0.000 0.000 0.000) 0.000 4.845 ( -0.000 -0.000 -0.000) 0.000 4.969 ( 0.000 0.000 0.000) 0.000 4.969 ( 0.000 0.000 0.000) 0.000 5.006 ( 0.000 0.000 -0.000) 0.000 5.006 ( 0.000 0.000 -0.000) 0.000 5.323 ( -0.000 0.000 -0.000) 0.000 5.323 ( -0.000 0.000 -0.000) 0.000 5.374 ( -0.000 0.000 0.000) 0.000 5.374 ( 0.000 0.000 -0.000) 0.000 5.618 ( 0.000 0.000 -0.000) 0.000 5.666 ( 0.000 0.000 0.000) 0.000 6.814 ( -0.000 0.000 -0.000) 0.000 6.814 ( -0.000 0.000 -0.000) 0.000 6.866 ( 0.000 -0.000 -0.000) 0.000 6.866 ( -0.000 0.000 -0.000) 0.000 7.871 ( -0.000 0.000 0.000) 0.000 8.000 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/24) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.216 ( 8.654 4.996 0.000) 9.993 0.328 ( 12.274 7.087 0.000) 14.173 0.511 ( 2.405 1.389 0.000) 2.777 0.573 ( 6.051 3.494 0.000) 6.987 0.617 ( -1.329 -0.767 0.000) 1.535 0.644 ( 24.475 14.130 0.000) 28.261 0.729 ( 2.024 1.168 0.000) 2.337 0.839 ( 4.501 2.599 0.000) 5.197 0.882 ( 2.985 1.723 0.000) 3.447 0.968 ( -2.710 -1.565 0.000) 3.129 0.992 ( 0.737 0.426 0.000) 0.851 1.182 ( 0.976 0.563 0.000) 1.127 1.401 ( -1.893 -1.093 0.000) 2.186 1.407 ( -2.468 -1.425 0.000) 2.849 1.434 ( -0.603 -0.348 0.000) 0.696 1.502 ( 3.562 2.056 0.000) 4.113 1.699 ( 0.962 0.555 0.000) 1.111 1.701 ( 0.150 0.086 0.000) 0.173 1.722 ( 2.135 1.232 0.000) 2.465 1.760 ( 1.839 1.062 0.000) 2.123 1.842 ( 2.787 1.609 0.000) 3.218 1.967 ( 2.045 1.181 0.000) 2.362 2.261 ( -1.956 -1.129 0.000) 2.258 2.278 ( -3.673 -2.121 0.000) 4.241 2.307 ( -4.062 -2.345 0.000) 4.691 2.324 ( -0.348 -0.201 0.000) 0.402 2.358 ( -1.063 -0.614 0.000) 1.227 2.484 ( -0.094 -0.054 0.000) 0.108 2.738 ( 0.052 0.030 0.000) 0.060 3.070 ( -4.865 -2.809 0.000) 5.617 4.160 ( 0.384 0.222 0.000) 0.444 4.177 ( 0.225 0.130 0.000) 0.260 4.242 ( 4.081 2.356 0.000) 4.712 4.330 ( 0.312 0.180 0.000) 0.361 4.331 ( 1.046 0.604 0.000) 1.208 4.452 ( -3.292 -1.901 0.000) 3.801 4.810 ( 2.837 1.638 0.000) 3.276 4.858 ( 1.160 0.670 0.000) 1.339 4.980 ( -0.172 -0.100 0.000) 0.199 4.984 ( 1.012 0.584 0.000) 1.169 5.008 ( 0.123 0.071 0.000) 0.142 5.042 ( 1.923 1.110 0.000) 2.221 5.329 ( 0.514 0.297 0.000) 0.593 5.355 ( 2.222 1.283 0.000) 2.565 5.384 ( 0.724 0.418 0.000) 0.837 5.476 ( 0.307 0.178 0.000) 0.355 5.584 ( -2.574 -1.486 0.000) 2.972 5.641 ( -1.982 -1.145 0.000) 2.289 6.811 ( -0.158 -0.091 0.000) 0.182 6.857 ( 3.305 1.908 0.000) 3.816 6.860 ( -0.494 -0.285 0.000) 0.571 7.073 ( 1.238 0.715 0.000) 1.430 7.842 ( -2.219 -1.281 0.000) 2.562 7.953 ( -3.666 -2.117 0.000) 4.233 ======================= Grid point 2 (3/24) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.413 ( 7.121 4.111 0.000) 8.222 0.524 ( 1.625 0.938 0.000) 1.876 0.602 ( 0.345 0.199 0.000) 0.398 0.628 ( 8.611 4.971 0.000) 9.943 0.698 ( 3.211 1.854 0.000) 3.708 0.780 ( 0.638 0.368 0.000) 0.737 0.886 ( -1.084 -0.626 0.000) 1.251 0.901 ( 0.961 0.555 0.000) 1.110 0.974 ( 3.151 1.819 0.000) 3.639 1.032 ( 2.120 1.224 0.000) 2.448 1.156 ( 12.899 7.447 0.000) 14.895 1.230 ( 4.309 2.488 0.000) 4.975 1.339 ( -2.892 -1.670 0.000) 3.340 1.356 ( -0.499 -0.288 0.000) 0.576 1.418 ( -0.642 -0.370 0.000) 0.741 1.528 ( -1.053 -0.608 0.000) 1.216 1.723 ( 1.502 0.867 0.000) 1.734 1.740 ( 0.473 0.273 0.000) 0.546 1.800 ( 3.947 2.279 0.000) 4.558 1.816 ( 2.682 1.549 0.000) 3.097 1.924 ( 3.851 2.223 0.000) 4.447 2.056 ( 5.089 2.938 0.000) 5.876 2.170 ( -4.289 -2.476 0.000) 4.953 2.213 ( -2.953 -1.705 0.000) 3.410 2.213 ( -1.741 -1.005 0.000) 2.011 2.300 ( -1.485 -0.857 0.000) 1.715 2.333 ( -0.895 -0.516 0.000) 1.033 2.460 ( -1.637 -0.945 0.000) 1.890 2.732 ( -0.592 -0.342 0.000) 0.683 2.927 ( -6.203 -3.581 0.000) 7.163 4.182 ( 1.559 0.900 0.000) 1.800 4.199 ( 1.498 0.865 0.000) 1.730 4.327 ( -0.400 -0.231 0.000) 0.462 4.327 ( 2.276 1.314 0.000) 2.629 4.348 ( -0.253 -0.146 0.000) 0.292 4.367 ( -3.301 -1.906 0.000) 3.811 4.869 ( 1.723 0.995 0.000) 1.989 4.887 ( 0.835 0.482 0.000) 0.964 5.003 ( 0.134 0.077 0.000) 0.154 5.014 ( 0.571 0.330 0.000) 0.660 5.024 ( 3.823 2.207 0.000) 4.414 5.112 ( 3.936 2.273 0.000) 4.545 5.351 ( 1.190 0.687 0.000) 1.375 5.406 ( 1.046 0.604 0.000) 1.208 5.406 ( 1.053 0.608 0.000) 1.216 5.471 ( -1.264 -0.730 0.000) 1.459 5.516 ( -2.097 -1.211 0.000) 2.422 5.577 ( -2.788 -1.610 0.000) 3.219 6.807 ( -0.147 -0.085 0.000) 0.170 6.845 ( -0.599 -0.346 0.000) 0.691 6.957 ( 4.368 2.522 0.000) 5.044 7.113 ( 1.835 1.059 0.000) 2.118 7.773 ( -3.127 -1.805 0.000) 3.611 7.839 ( -5.058 -2.920 0.000) 5.841 ======================= Grid point 3 (4/24) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.537 ( 0.279 0.161 0.000) 0.322 0.545 ( 3.102 1.791 0.000) 3.582 0.608 ( -0.077 -0.044 0.000) 0.089 0.707 ( -1.135 -0.655 0.000) 1.310 0.769 ( -0.109 -0.063 0.000) 0.126 0.825 ( 6.759 3.902 0.000) 7.804 0.932 ( 0.817 0.472 0.000) 0.943 0.936 ( 2.606 1.505 0.000) 3.010 0.961 ( -2.849 -1.645 0.000) 3.290 1.082 ( 0.830 0.479 0.000) 0.958 1.175 ( -1.778 -1.027 0.000) 2.053 1.287 ( -0.883 -0.510 0.000) 1.020 1.354 ( 1.207 0.697 0.000) 1.393 1.399 ( -0.669 -0.386 0.000) 0.773 1.451 ( 1.890 1.091 0.000) 2.182 1.546 ( 3.920 2.263 0.000) 4.526 1.719 ( -0.830 -0.479 0.000) 0.959 1.758 ( 0.968 0.559 0.000) 1.118 1.834 ( -1.194 -0.690 0.000) 1.379 1.878 ( 1.752 1.011 0.000) 2.023 2.020 ( 2.849 1.645 0.000) 3.290 2.100 ( -0.788 -0.455 0.000) 0.910 2.148 ( 0.079 0.046 0.000) 0.092 2.181 ( -0.764 -0.441 0.000) 0.882 2.213 ( 3.303 1.907 0.000) 3.814 2.265 ( -1.003 -0.579 0.000) 1.158 2.321 ( 0.173 0.100 0.000) 0.200 2.413 ( -2.015 -1.163 0.000) 2.326 2.719 ( 0.016 0.009 0.000) 0.018 2.793 ( -4.018 -2.320 0.000) 4.639 4.227 ( 1.527 0.881 0.000) 1.763 4.250 ( 2.618 1.512 0.000) 3.023 4.305 ( -1.397 -0.806 0.000) 1.613 4.309 ( -0.850 -0.491 0.000) 0.982 4.315 ( -2.251 -1.300 0.000) 2.599 4.357 ( 0.398 0.230 0.000) 0.460 4.887 ( 0.180 0.104 0.000) 0.207 4.906 ( 0.461 0.266 0.000) 0.533 5.000 ( -0.128 -0.074 0.000) 0.148 5.026 ( 0.235 0.136 0.000) 0.271 5.132 ( 3.872 2.235 0.000) 4.471 5.213 ( 3.075 1.776 0.000) 3.551 5.377 ( -1.527 -0.882 0.000) 1.764 5.377 ( 0.696 0.402 0.000) 0.803 5.420 ( -3.233 -1.867 0.000) 3.733 5.426 ( 0.474 0.274 0.000) 0.548 5.497 ( 0.334 0.193 0.000) 0.385 5.535 ( -0.610 -0.352 0.000) 0.705 6.805 ( -0.043 -0.025 0.000) 0.049 6.834 ( -0.265 -0.153 0.000) 0.306 7.046 ( 2.277 1.315 0.000) 2.630 7.151 ( 0.923 0.533 0.000) 1.065 7.709 ( -1.607 -0.928 0.000) 1.856 7.737 ( -2.568 -1.482 0.000) 2.965 ======================= Grid point 8 (5/24) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.354 ( 3.891 6.740 0.000) 7.782 0.568 ( 1.639 2.839 0.000) 3.278 0.581 ( 1.853 3.210 0.000) 3.707 0.589 ( 4.423 7.661 0.000) 8.846 0.672 ( 2.086 3.613 0.000) 4.172 0.764 ( 1.475 2.555 0.000) 2.950 0.863 ( -2.657 -4.603 0.000) 5.315 0.894 ( 3.390 5.872 0.000) 6.780 0.961 ( 2.308 3.998 0.000) 4.616 1.029 ( 4.386 7.598 0.000) 8.773 1.071 ( 7.370 12.765 0.000) 14.740 1.197 ( 0.456 0.791 0.000) 0.913 1.362 ( -1.115 -1.931 0.000) 2.230 1.366 ( -1.274 -2.206 0.000) 2.547 1.421 ( -0.456 -0.790 0.000) 0.912 1.538 ( 0.244 0.423 0.000) 0.489 1.700 ( 0.362 0.628 0.000) 0.725 1.757 ( 0.966 1.673 0.000) 1.932 1.760 ( 2.332 4.039 0.000) 4.664 1.790 ( 1.038 1.798 0.000) 2.076 1.911 ( 2.570 4.451 0.000) 5.140 2.002 ( 1.365 2.363 0.000) 2.729 2.190 ( -2.507 -4.343 0.000) 5.014 2.223 ( -2.143 -3.712 0.000) 4.287 2.276 ( -0.513 -0.888 0.000) 1.026 2.289 ( -0.725 -1.256 0.000) 1.450 2.338 ( -0.885 -1.532 0.000) 1.769 2.466 ( -0.995 -1.724 0.000) 1.991 2.735 ( -0.291 -0.503 0.000) 0.581 2.967 ( -3.586 -6.211 0.000) 7.172 4.175 ( 0.562 0.973 0.000) 1.124 4.185 ( 0.669 1.159 0.000) 1.338 4.307 ( 1.104 1.911 0.000) 2.207 4.337 ( 0.732 1.267 0.000) 1.464 4.349 ( 0.376 0.652 0.000) 0.752 4.400 ( -1.620 -2.806 0.000) 3.240 4.863 ( 1.416 2.452 0.000) 2.831 4.880 ( 0.805 1.395 0.000) 1.610 4.966 ( 0.219 0.380 0.000) 0.438 4.980 ( -0.066 -0.115 0.000) 0.133 5.036 ( 1.836 3.181 0.000) 3.673 5.116 ( 1.984 3.436 0.000) 3.968 5.328 ( -0.164 -0.284 0.000) 0.328 5.375 ( -0.124 -0.215 0.000) 0.248 5.413 ( 1.484 2.570 0.000) 2.967 5.486 ( 0.173 0.300 0.000) 0.347 5.530 ( -1.754 -3.038 0.000) 3.508 5.597 ( -1.574 -2.726 0.000) 3.148 6.813 ( 0.136 0.236 0.000) 0.272 6.853 ( -0.157 -0.272 0.000) 0.314 6.927 ( 2.454 4.250 0.000) 4.908 7.097 ( 0.830 1.438 0.000) 1.660 7.793 ( -1.829 -3.167 0.000) 3.657 7.871 ( -2.967 -5.140 0.000) 5.935 ======================= Grid point 9 (6/24) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 75 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.475 ( 3.631 2.003 0.000) 4.147 0.554 ( -0.808 -0.518 0.000) 0.960 0.613 ( 0.368 1.114 0.000) 1.173 0.683 ( -0.205 -3.727 0.000) 3.733 0.734 ( 1.891 -0.701 0.000) 2.016 0.811 ( 1.626 4.473 0.000) 4.760 0.867 ( 3.070 2.039 0.000) 3.686 0.980 ( 0.154 5.141 0.000) 5.143 1.037 ( -1.642 3.927 0.000) 4.256 1.103 ( 0.445 3.819 0.000) 3.845 1.205 ( 0.290 1.214 0.000) 1.248 1.329 ( -0.084 1.373 0.000) 1.376 1.341 ( 0.971 2.536 0.000) 2.715 1.391 ( 3.009 2.210 0.000) 3.734 1.412 ( 3.380 3.634 0.000) 4.963 1.525 ( -0.696 0.194 0.000) 0.722 1.707 ( -0.116 -2.090 0.000) 2.093 1.779 ( -0.320 2.997 0.000) 3.014 1.815 ( 0.514 -1.692 0.000) 1.769 1.852 ( 3.592 2.161 0.000) 4.192 2.000 ( 1.239 4.297 0.000) 4.472 2.085 ( 4.698 0.279 0.000) 4.706 2.103 ( -0.517 -3.170 0.000) 3.212 2.170 ( -0.526 -0.989 0.000) 1.120 2.242 ( -2.693 1.569 0.000) 3.117 2.286 ( -0.407 0.038 0.000) 0.409 2.303 ( 0.250 -2.307 0.000) 2.320 2.415 ( -1.540 -3.419 0.000) 3.750 2.717 ( -0.396 -0.943 0.000) 1.023 2.832 ( -4.008 -4.826 0.000) 6.274 4.204 ( 1.460 1.037 0.000) 1.791 4.224 ( 1.777 1.396 0.000) 2.260 4.303 ( -0.722 -1.382 0.000) 1.559 4.340 ( -2.530 -0.747 0.000) 2.638 4.352 ( -0.912 0.933 0.000) 1.305 4.377 ( 0.358 2.138 0.000) 2.168 4.887 ( -0.083 0.474 0.000) 0.481 4.903 ( 0.448 0.713 0.000) 0.842 5.004 ( 1.484 1.336 0.000) 1.997 5.013 ( 1.950 1.555 0.000) 2.494 5.112 ( 1.984 4.327 0.000) 4.760 5.183 ( 2.080 2.984 0.000) 3.637 5.326 ( 1.464 -2.254 0.000) 2.687 5.358 ( 0.129 -3.594 0.000) 3.596 5.422 ( -1.219 -1.589 0.000) 2.003 5.454 ( -2.624 -2.698 0.000) 3.763 5.511 ( 0.479 2.031 0.000) 2.087 5.557 ( -0.765 0.361 0.000) 0.846 6.819 ( -0.327 1.182 0.000) 1.226 6.849 ( -0.478 0.605 0.000) 0.771 7.020 ( 2.749 3.104 0.000) 4.146 7.129 ( 1.369 0.304 0.000) 1.402 7.722 ( -1.898 -2.666 0.000) 3.273 7.759 ( -3.033 -4.192 0.000) 5.174 ======================= Grid point 10 (7/24) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.518 ( 0.563 -0.974 0.000) 1.125 0.546 ( 0.758 -1.313 0.000) 1.516 0.612 ( -0.001 0.003 0.000) 0.003 0.662 ( 1.612 -2.792 0.000) 3.224 0.738 ( 1.448 -2.508 0.000) 2.896 0.878 ( 1.338 -2.317 0.000) 2.675 0.907 ( 0.821 -1.422 0.000) 1.643 0.995 ( -3.319 5.748 0.000) 6.637 1.022 ( -3.273 5.669 0.000) 6.546 1.124 ( -1.509 2.614 0.000) 3.018 1.231 ( -2.962 5.130 0.000) 5.923 1.336 ( -2.272 3.936 0.000) 4.545 1.365 ( -0.113 0.195 0.000) 0.225 1.385 ( 0.384 -0.665 0.000) 0.768 1.484 ( -1.337 2.315 0.000) 2.673 1.565 ( 0.175 -0.302 0.000) 0.349 1.687 ( 1.363 -2.361 0.000) 2.727 1.779 ( 2.004 -3.471 0.000) 4.008 1.802 ( -1.850 3.204 0.000) 3.700 1.911 ( -0.969 1.678 0.000) 1.937 2.027 ( 1.130 -1.958 0.000) 2.261 2.078 ( -0.103 0.179 0.000) 0.207 2.142 ( 0.268 -0.464 0.000) 0.535 2.201 ( -0.655 1.134 0.000) 1.310 2.201 ( 0.910 -1.577 0.000) 1.821 2.288 ( -0.702 1.215 0.000) 1.403 2.305 ( 0.812 -1.407 0.000) 1.624 2.356 ( 2.044 -3.541 0.000) 4.089 2.715 ( 0.289 -0.500 0.000) 0.578 2.752 ( 0.611 -1.059 0.000) 1.223 4.236 ( 0.073 -0.126 0.000) 0.146 4.265 ( 0.229 -0.397 0.000) 0.459 4.271 ( 0.985 -1.706 0.000) 1.969 4.304 ( -0.951 1.648 0.000) 1.902 4.345 ( -1.531 2.651 0.000) 3.062 4.391 ( -0.903 1.565 0.000) 1.807 4.887 ( -0.073 0.126 0.000) 0.146 4.912 ( -0.046 0.079 0.000) 0.092 5.026 ( -1.325 2.295 0.000) 2.649 5.041 ( -0.868 1.503 0.000) 1.736 5.173 ( -0.566 0.981 0.000) 1.133 5.233 ( -0.136 0.236 0.000) 0.272 5.323 ( 1.995 -3.455 0.000) 3.989 5.337 ( 1.696 -2.938 0.000) 3.392 5.380 ( 0.873 -1.512 0.000) 1.746 5.402 ( 1.156 -2.003 0.000) 2.313 5.533 ( -1.229 2.128 0.000) 2.458 5.555 ( -0.829 1.436 0.000) 1.658 6.823 ( -0.876 1.516 0.000) 1.751 6.848 ( -0.640 1.108 0.000) 1.279 7.068 ( -0.445 0.771 0.000) 0.890 7.144 ( 0.610 -1.056 0.000) 1.220 7.687 ( 0.571 -0.989 0.000) 1.142 7.702 ( 0.857 -1.484 0.000) 1.714 ======================= Grid point 17 (8/24) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.516 ( 1.201 2.081 0.000) 2.402 0.528 ( -0.570 -0.988 0.000) 1.141 0.594 ( -1.919 -3.324 0.000) 3.838 0.632 ( 0.194 0.337 0.000) 0.389 0.699 ( -0.821 -1.422 0.000) 1.642 0.844 ( 0.630 1.091 0.000) 1.259 0.889 ( -1.147 -1.987 0.000) 2.294 1.081 ( 1.847 3.199 0.000) 3.694 1.090 ( 0.920 1.593 0.000) 1.839 1.158 ( 0.641 1.110 0.000) 1.282 1.306 ( 3.048 5.280 0.000) 6.097 1.360 ( 0.432 0.748 0.000) 0.863 1.377 ( 0.449 0.777 0.000) 0.897 1.383 ( -0.112 -0.195 0.000) 0.225 1.535 ( 1.004 1.738 0.000) 2.007 1.561 ( 3.388 5.869 0.000) 6.777 1.641 ( -1.808 -3.131 0.000) 3.615 1.731 ( -2.988 -5.176 0.000) 5.977 1.875 ( 2.439 4.224 0.000) 4.877 1.898 ( 1.594 2.760 0.000) 3.187 2.008 ( -2.590 -4.486 0.000) 5.180 2.100 ( 1.194 2.068 0.000) 2.388 2.116 ( 1.971 3.414 0.000) 3.943 2.168 ( 0.556 0.963 0.000) 1.112 2.227 ( -1.904 -3.297 0.000) 3.807 2.273 ( 0.063 0.109 0.000) 0.126 2.296 ( 0.224 0.388 0.000) 0.448 2.320 ( -2.550 -4.417 0.000) 5.100 2.706 ( 0.140 0.243 0.000) 0.280 2.745 ( -1.703 -2.949 0.000) 3.406 4.231 ( 0.822 1.423 0.000) 1.643 4.244 ( 0.247 0.427 0.000) 0.493 4.270 ( -0.948 -1.642 0.000) 1.896 4.331 ( -0.081 -0.140 0.000) 0.162 4.383 ( 0.798 1.382 0.000) 1.596 4.403 ( 0.171 0.296 0.000) 0.342 4.895 ( 0.167 0.289 0.000) 0.333 4.907 ( -0.125 -0.216 0.000) 0.250 5.060 ( 1.924 3.332 0.000) 3.847 5.073 ( 1.440 2.494 0.000) 2.880 5.176 ( 0.356 0.617 0.000) 0.712 5.229 ( 0.618 1.070 0.000) 1.235 5.287 ( -0.255 -0.442 0.000) 0.510 5.289 ( -0.796 -1.379 0.000) 1.592 5.367 ( -2.457 -4.255 0.000) 4.913 5.372 ( -1.250 -2.165 0.000) 2.500 5.558 ( 0.917 1.588 0.000) 1.834 5.572 ( 0.253 0.438 0.000) 0.505 6.847 ( 0.628 1.088 0.000) 1.257 6.861 ( 0.154 0.267 0.000) 0.308 7.080 ( 1.313 2.275 0.000) 2.627 7.124 ( -0.443 -0.767 0.000) 0.886 7.675 ( -0.802 -1.388 0.000) 1.603 7.685 ( -1.242 -2.151 0.000) 2.484 ======================= Grid point 55 (9/24) ======================= q-point: ( 0.00 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.244 ( 0.000 0.000 10.543) 10.543 0.244 ( -0.000 0.000 10.543) 10.543 0.359 ( 0.000 -0.000 16.423) 16.423 0.483 ( -0.000 0.000 0.379) 0.379 0.483 ( 0.000 -0.000 0.379) 0.379 0.650 ( 0.000 0.000 1.183) 1.183 0.650 ( 0.000 0.000 1.183) 1.183 0.779 ( -0.000 -0.000 0.375) 0.375 0.822 ( 0.000 0.000 -0.538) 0.538 0.945 ( -0.000 0.000 -4.069) 4.069 0.945 ( 0.000 0.000 -4.069) 4.069 1.134 ( 0.000 0.000 -5.051) 5.051 1.426 ( 0.000 0.000 -0.453) 0.453 1.426 ( 0.000 0.000 -0.453) 0.453 1.428 ( 0.000 0.000 -0.503) 0.503 1.428 ( -0.000 0.000 -0.503) 0.503 1.714 ( -0.000 -0.000 0.623) 0.623 1.714 ( 0.000 0.000 0.623) 0.623 1.733 ( -0.000 0.000 0.112) 0.112 1.733 ( 0.000 0.000 0.112) 0.112 1.905 ( -0.000 -0.000 4.227) 4.227 1.969 ( -0.000 0.000 2.309) 2.309 2.296 ( 0.000 0.000 0.201) 0.201 2.296 ( 0.000 -0.000 0.201) 0.201 2.318 ( -0.000 -0.000 -0.313) 0.313 2.318 ( 0.000 0.000 -0.313) 0.313 2.343 ( 0.000 -0.000 -3.024) 3.024 2.344 ( -0.000 -0.000 -2.905) 2.905 3.111 ( -0.000 0.000 -0.059) 0.059 3.130 ( 0.000 0.000 -0.482) 0.482 4.158 ( 0.000 -0.000 0.249) 0.249 4.158 ( 0.000 0.000 0.249) 0.249 4.178 ( 0.000 0.000 -0.162) 0.162 4.178 ( -0.000 0.000 -0.162) 0.162 4.309 ( 0.000 0.000 0.124) 0.124 4.326 ( 0.000 0.000 -0.244) 0.244 4.775 ( -0.000 -0.000 0.894) 0.894 4.840 ( 0.000 -0.000 -0.494) 0.494 4.971 ( 0.000 0.000 0.154) 0.154 4.971 ( -0.000 -0.000 0.154) 0.154 5.001 ( 0.000 -0.000 -0.486) 0.486 5.001 ( -0.000 0.000 -0.486) 0.486 5.327 ( 0.000 0.000 0.294) 0.294 5.327 ( -0.000 0.000 0.294) 0.294 5.368 ( -0.000 -0.000 -0.583) 0.583 5.368 ( 0.000 0.000 -0.583) 0.583 5.630 ( 0.000 -0.000 1.099) 1.099 5.738 ( -0.000 0.000 -1.185) 1.185 6.818 ( 0.000 -0.000 0.432) 0.432 6.818 ( 0.000 -0.000 0.432) 0.432 6.861 ( 0.000 0.000 -0.479) 0.479 6.861 ( 0.000 0.000 -0.479) 0.479 7.883 ( 0.000 -0.000 0.756) 0.756 7.984 ( 0.000 0.000 -1.405) 1.405 ======================= Grid point 56 (10/24) ======================= q-point: ( 0.14 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.306 ( 5.250 3.031 8.343) 10.313 0.401 ( 9.407 5.431 5.630) 12.235 0.519 ( 2.515 1.452 0.600) 2.965 0.581 ( 6.112 3.529 0.633) 7.086 0.629 ( -1.329 -0.767 1.283) 2.001 0.702 ( 12.414 7.167 1.471) 14.410 0.731 ( 10.904 6.295 3.844) 13.165 0.839 ( 3.827 2.210 -0.122) 4.421 0.874 ( 2.903 1.676 -0.713) 3.427 0.920 ( -2.291 -1.323 -4.182) 4.949 0.957 ( 1.902 1.098 -3.016) 3.731 1.148 ( 1.892 1.093 -3.285) 3.946 1.400 ( -1.595 -0.921 -0.020) 1.841 1.412 ( -1.477 -0.853 0.500) 1.778 1.427 ( -0.552 -0.319 -0.550) 0.842 1.489 ( 3.491 2.015 -1.153) 4.192 1.708 ( 0.102 0.059 0.656) 0.666 1.720 ( 1.376 0.794 0.664) 1.722 1.757 ( 1.571 0.907 -0.274) 1.835 1.759 ( 1.489 0.860 2.103) 2.716 1.873 ( 0.479 0.277 3.874) 3.914 1.993 ( 1.970 1.138 2.265) 3.211 2.239 ( -2.822 -1.629 -3.157) 4.537 2.267 ( -1.791 -1.034 0.518) 2.132 2.270 ( -3.516 -2.030 -1.804) 4.443 2.318 ( -0.473 -0.273 -0.576) 0.794 2.336 ( -0.730 -0.421 -1.626) 1.831 2.439 ( 1.981 1.144 -3.237) 3.963 2.868 ( -6.775 -3.911 7.473) 10.819 3.062 ( -4.883 -2.819 -0.708) 5.683 4.162 ( 0.396 0.229 0.239) 0.516 4.177 ( 0.263 0.152 -0.055) 0.309 4.245 ( 4.281 2.472 0.347) 4.956 4.329 ( 0.418 0.241 -0.031) 0.484 4.331 ( 0.996 0.575 -0.057) 1.151 4.384 ( 2.880 1.663 -4.826) 5.861 4.803 ( 2.820 1.628 -0.075) 3.257 4.852 ( 1.136 0.656 -0.575) 1.433 4.976 ( 0.169 0.098 -0.201) 0.281 4.984 ( 0.904 0.522 0.002) 1.044 5.003 ( 0.176 0.102 -0.422) 0.468 5.014 ( 1.823 1.053 -1.550) 2.615 5.333 ( 0.601 0.347 0.405) 0.804 5.359 ( 2.285 1.319 0.382) 2.666 5.378 ( 0.766 0.442 -0.543) 1.038 5.434 ( 2.631 1.519 -2.259) 3.786 5.594 ( -2.733 -1.578 0.921) 3.287 5.663 ( -3.854 -2.225 0.514) 4.480 6.816 ( -0.200 -0.115 0.383) 0.448 6.855 ( -0.472 -0.272 -0.451) 0.707 6.864 ( 3.458 1.997 0.659) 4.047 6.996 ( 5.048 2.914 -4.569) 7.406 7.851 ( -2.362 -1.364 0.786) 2.839 7.940 ( -3.454 -1.994 -1.165) 4.155 ======================= Grid point 57 (11/24) ======================= q-point: ( 0.29 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.461 ( 6.265 3.617 4.501) 8.520 0.546 ( 0.811 0.468 1.714) 1.953 0.601 ( 0.936 0.540 -0.054) 1.082 0.655 ( 7.512 4.337 2.484) 9.023 0.698 ( 2.477 1.430 0.105) 2.862 0.777 ( 1.405 0.811 -0.341) 1.658 0.860 ( -0.177 -0.102 -2.476) 2.484 0.878 ( 0.970 0.560 -2.037) 2.325 0.976 ( 3.378 1.950 0.239) 3.908 1.026 ( 2.703 1.561 -0.618) 3.182 1.160 ( 11.413 6.589 0.409) 13.185 1.211 ( 4.131 2.385 -1.312) 4.947 1.347 ( -2.521 -1.456 0.720) 2.999 1.382 ( 0.366 0.211 2.153) 2.194 1.411 ( -0.711 -0.410 -0.631) 1.035 1.526 ( -0.457 -0.264 -0.640) 0.829 1.730 ( 1.575 0.909 0.684) 1.943 1.772 ( -1.053 -0.608 2.601) 2.871 1.803 ( 3.980 2.298 0.453) 4.618 1.807 ( 2.460 1.420 -0.812) 2.954 1.935 ( 4.212 2.432 0.928) 4.951 2.064 ( 3.341 1.929 0.499) 3.890 2.161 ( -3.070 -1.772 -0.695) 3.612 2.189 ( -1.821 -1.051 -1.676) 2.689 2.221 ( -1.734 -1.001 0.724) 2.129 2.292 ( -1.508 -0.871 -0.746) 1.895 2.312 ( -1.016 -0.586 -1.511) 1.913 2.434 ( -1.496 -0.864 -2.434) 2.985 2.777 ( -2.053 -1.185 3.565) 4.281 2.925 ( -5.781 -3.337 -0.337) 6.683 4.185 ( 1.548 0.894 0.257) 1.806 4.199 ( 1.510 0.872 -0.034) 1.744 4.328 ( 1.946 1.123 0.137) 2.251 4.329 ( -0.365 -0.211 0.203) 0.468 4.346 ( -0.251 -0.145 -0.168) 0.335 4.361 ( -2.479 -1.431 -0.505) 2.906 4.868 ( 1.857 1.072 -0.071) 2.145 4.885 ( 0.999 0.577 -0.237) 1.178 5.001 ( 0.210 0.121 -0.097) 0.261 5.011 ( 0.592 0.342 -0.294) 0.744 5.023 ( 3.893 2.248 -0.013) 4.495 5.089 ( 4.266 2.463 -1.867) 5.268 5.357 ( 1.211 0.699 0.529) 1.495 5.401 ( 1.055 0.609 -0.435) 1.294 5.412 ( 1.239 0.715 0.548) 1.532 5.466 ( -0.326 -0.188 -0.502) 0.628 5.519 ( -2.479 -1.431 0.333) 2.882 5.571 ( -3.416 -1.972 -0.617) 3.992 6.810 ( -0.204 -0.118 0.270) 0.358 6.841 ( -0.569 -0.328 -0.387) 0.762 6.967 ( 4.439 2.563 0.929) 5.209 7.087 ( 2.899 1.674 -2.112) 3.958 7.779 ( -3.232 -1.866 0.560) 3.774 7.832 ( -4.775 -2.757 -0.602) 5.546 ======================= Grid point 58 (12/24) ======================= q-point: ( 0.43 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.550 ( 0.100 0.057 1.136) 1.141 0.573 ( 2.323 1.341 3.028) 4.045 0.607 ( -0.162 -0.094 -0.163) 0.248 0.711 ( -0.062 -0.036 0.601) 0.605 0.790 ( 0.129 0.074 0.948) 0.959 0.834 ( 6.069 3.504 0.805) 7.054 0.910 ( 2.241 1.294 -2.459) 3.570 0.912 ( 0.680 0.392 -1.302) 1.520 0.982 ( -1.683 -0.972 1.421) 2.407 1.074 ( 0.604 0.349 -0.796) 1.059 1.189 ( -1.056 -0.610 1.419) 1.871 1.302 ( -0.786 -0.454 1.290) 1.577 1.344 ( 1.397 0.807 -1.208) 2.015 1.392 ( -0.620 -0.358 -0.667) 0.979 1.452 ( 2.233 1.289 0.612) 2.650 1.552 ( 3.035 1.752 0.107) 3.507 1.719 ( -1.662 -0.960 0.452) 1.972 1.768 ( 1.066 0.615 0.904) 1.527 1.829 ( -1.733 -1.001 -0.652) 2.105 1.864 ( 1.645 0.949 -1.163) 2.227 2.040 ( 2.811 1.623 1.316) 3.503 2.098 ( 1.168 0.674 -0.915) 1.630 2.134 ( -0.837 -0.483 0.130) 0.975 2.188 ( -0.798 -0.461 0.660) 1.134 2.217 ( 2.658 1.535 0.301) 3.084 2.257 ( -0.991 -0.572 -0.742) 1.364 2.298 ( 0.396 0.228 -1.764) 1.822 2.390 ( -1.941 -1.121 -2.140) 3.099 2.745 ( -0.378 -0.218 2.087) 2.132 2.803 ( -3.482 -2.010 0.756) 4.091 4.230 ( 1.604 0.926 0.289) 1.875 4.250 ( 2.550 1.472 -0.047) 2.945 4.309 ( -1.263 -0.729 0.344) 1.499 4.310 ( -1.044 -0.603 0.060) 1.207 4.315 ( -2.082 -1.202 -0.038) 2.404 4.351 ( 0.203 0.117 -0.544) 0.593 4.888 ( 0.216 0.125 0.116) 0.275 4.904 ( 0.448 0.259 -0.177) 0.547 5.002 ( -0.040 -0.023 0.144) 0.151 5.022 ( 0.243 0.140 -0.294) 0.406 5.134 ( 3.930 2.269 0.157) 4.540 5.197 ( 3.387 1.955 -1.336) 4.132 5.381 ( -1.864 -1.076 0.422) 2.193 5.383 ( 0.679 0.392 0.509) 0.935 5.418 ( -3.351 -1.934 -0.205) 3.874 5.422 ( 0.504 0.291 -0.332) 0.670 5.498 ( 0.321 0.186 0.071) 0.378 5.528 ( -0.308 -0.178 -0.577) 0.677 6.807 ( -0.070 -0.041 0.171) 0.189 6.830 ( -0.250 -0.144 -0.330) 0.438 7.057 ( 2.259 1.304 0.941) 2.773 7.140 ( 1.227 0.708 -0.913) 1.685 7.713 ( -1.661 -0.959 0.361) 1.952 7.735 ( -2.431 -1.404 -0.116) 2.810 ======================= Grid point 63 (13/24) ======================= q-point: ( 0.14 0.14 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.409 ( 3.457 5.987 5.054) 8.564 0.575 ( 1.568 2.716 0.654) 3.204 0.585 ( 1.182 2.048 0.298) 2.383 0.629 ( 3.556 6.159 3.163) 7.783 0.681 ( 2.568 4.449 1.090) 5.251 0.758 ( 1.430 2.477 -0.530) 2.909 0.839 ( -1.474 -2.553 -2.277) 3.725 0.895 ( 2.771 4.799 -0.056) 5.542 0.963 ( 2.328 4.033 0.003) 4.657 1.003 ( 2.687 4.653 -2.447) 5.904 1.073 ( 8.419 14.583 0.595) 16.849 1.178 ( 0.870 1.507 -1.664) 2.408 1.373 ( -0.993 -1.720 0.598) 2.075 1.380 ( -0.706 -1.223 1.659) 2.178 1.415 ( -0.399 -0.692 -0.556) 0.973 1.524 ( 0.201 0.347 -1.314) 1.373 1.732 ( 0.796 1.378 1.411) 2.126 1.763 ( 2.118 3.669 0.285) 4.246 1.779 ( 0.094 0.163 2.685) 2.692 1.783 ( 1.112 1.926 -0.476) 2.274 1.914 ( 2.132 3.693 0.628) 4.311 2.027 ( 1.192 2.065 1.971) 3.094 2.183 ( -2.550 -4.416 -0.891) 5.177 2.190 ( -1.402 -2.428 -2.267) 3.605 2.272 ( -0.286 -0.495 -0.174) 0.598 2.285 ( -0.657 -1.138 -0.319) 1.352 2.317 ( -0.960 -1.663 -1.539) 2.461 2.440 ( -0.786 -1.361 -2.574) 3.016 2.791 ( -1.562 -2.705 4.364) 5.366 2.962 ( -3.414 -5.913 -0.542) 6.849 4.178 ( 0.565 0.979 0.204) 1.149 4.186 ( 0.656 1.136 0.042) 1.312 4.311 ( 1.121 1.941 0.419) 2.280 4.338 ( 0.675 1.170 0.089) 1.354 4.347 ( 0.383 0.664 -0.149) 0.781 4.386 ( -0.885 -1.533 -1.214) 2.146 4.859 ( 1.533 2.655 -0.330) 3.083 4.875 ( 0.866 1.499 -0.445) 1.787 4.964 ( 0.207 0.359 -0.087) 0.424 4.977 ( -0.041 -0.071 -0.253) 0.265 5.036 ( 1.870 3.239 -0.013) 3.740 5.092 ( 2.197 3.805 -1.818) 4.755 5.334 ( -0.118 -0.204 0.518) 0.569 5.373 ( -0.150 -0.261 -0.269) 0.404 5.417 ( 1.604 2.779 0.327) 3.225 5.474 ( 0.846 1.466 -0.996) 1.964 5.537 ( -1.908 -3.305 0.615) 3.866 5.595 ( -2.095 -3.629 -0.343) 4.204 6.816 ( 0.112 0.194 0.315) 0.386 6.849 ( -0.125 -0.216 -0.363) 0.440 6.936 ( 2.494 4.319 0.824) 5.055 7.063 ( 1.670 2.892 -2.661) 4.270 7.799 ( -1.896 -3.284 0.620) 3.843 7.863 ( -2.802 -4.854 -0.761) 5.656 ======================= Grid point 64 (14/24) ======================= q-point: ( 0.29 0.14 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 123 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.519 ( 2.829 2.179 3.952) 5.326 0.568 ( -0.748 -0.511 1.316) 1.598 0.614 ( 0.097 0.972 0.094) 0.981 0.671 ( 1.198 -2.657 0.035) 2.915 0.779 ( 1.156 0.419 2.242) 2.557 0.809 ( 2.249 3.548 -0.204) 4.206 0.868 ( 2.400 3.550 0.136) 4.287 0.972 ( 0.069 5.318 -0.808) 5.379 1.027 ( -1.786 3.063 -0.742) 3.622 1.084 ( 1.049 3.299 -1.512) 3.777 1.209 ( 1.032 1.830 0.520) 2.164 1.307 ( 3.034 3.238 -1.545) 4.698 1.351 ( -0.534 1.882 0.788) 2.109 1.390 ( 0.207 -0.741 -0.133) 0.781 1.428 ( 3.158 4.388 1.103) 5.517 1.517 ( 0.753 -0.322 -0.822) 1.161 1.732 ( -1.926 -2.696 2.212) 3.983 1.779 ( 0.093 2.875 0.156) 2.881 1.807 ( 0.829 -1.400 -0.621) 1.741 1.854 ( 2.788 2.173 0.074) 3.536 2.009 ( 2.579 3.960 0.445) 4.747 2.074 ( 1.537 -0.508 -0.996) 1.900 2.113 ( 0.776 -2.266 0.977) 2.587 2.162 ( 0.643 -0.629 -0.499) 1.028 2.241 ( -2.475 0.957 -0.028) 2.654 2.275 ( -0.686 -0.920 -1.339) 1.764 2.288 ( 0.379 -0.735 -0.602) 1.023 2.393 ( -1.469 -3.095 -2.032) 3.983 2.747 ( -0.831 -1.399 2.465) 2.953 2.837 ( -3.574 -4.380 0.309) 5.661 4.207 ( 1.485 1.061 0.254) 1.843 4.223 ( 1.773 1.361 -0.060) 2.236 4.307 ( -0.861 -1.313 0.360) 1.611 4.338 ( -2.298 -0.878 -0.270) 2.475 4.353 ( -0.889 1.067 0.105) 1.393 4.373 ( 0.193 2.004 -0.315) 2.038 4.887 ( 0.061 0.500 -0.005) 0.504 4.901 ( 0.496 0.652 -0.249) 0.857 5.003 ( 1.577 1.415 -0.104) 2.121 5.010 ( 1.991 1.534 -0.265) 2.527 5.113 ( 2.033 4.405 0.122) 4.853 5.168 ( 2.223 3.455 -1.253) 4.295 5.332 ( 1.390 -2.336 0.497) 2.763 5.358 ( 0.376 -3.471 0.006) 3.491 5.423 ( -1.579 -1.858 0.191) 2.445 5.453 ( -2.730 -2.889 -0.139) 3.977 5.512 ( 0.432 2.036 0.142) 2.086 5.548 ( -0.509 0.724 -0.741) 1.154 6.822 ( -0.365 1.186 0.238) 1.264 6.846 ( -0.496 0.730 -0.257) 0.919 7.029 ( 2.786 2.966 0.831) 4.153 7.116 ( 1.775 0.836 -1.145) 2.272 7.727 ( -1.960 -2.771 0.391) 3.417 7.756 ( -2.871 -3.993 -0.239) 4.924 ======================= Grid point 65 (15/24) ======================= q-point: ( 0.43 0.14 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.545 ( 0.046 -0.079 2.097) 2.099 0.565 ( 0.773 -1.339 2.646) 3.065 0.615 ( -0.333 0.576 0.074) 0.670 0.676 ( 1.582 -2.740 1.323) 3.429 0.783 ( 0.357 -0.618 2.535) 2.633 0.890 ( 0.788 -1.365 -1.699) 2.317 0.891 ( 1.538 -2.665 1.527) 3.435 0.989 ( -3.475 6.019 -0.499) 6.968 1.016 ( -3.023 5.236 -0.370) 6.058 1.104 ( -1.183 2.049 -1.706) 2.917 1.238 ( -2.583 4.473 0.768) 5.222 1.343 ( -1.891 3.276 0.621) 3.833 1.355 ( -0.245 0.424 -1.250) 1.342 1.381 ( 0.371 -0.642 -0.097) 0.748 1.487 ( -1.144 1.981 0.409) 2.324 1.560 ( 0.391 -0.678 -0.549) 0.956 1.691 ( 1.027 -1.779 0.524) 2.120 1.770 ( 1.890 -3.273 -0.820) 3.867 1.807 ( -1.691 2.928 0.554) 3.426 1.901 ( -1.039 1.799 -1.002) 2.306 2.051 ( 1.074 -1.859 1.443) 2.587 2.086 ( 0.165 -0.285 1.262) 1.305 2.115 ( 0.284 -0.492 -2.112) 2.187 2.195 ( 0.938 -1.625 -0.671) 1.993 2.205 ( -0.323 0.560 0.406) 0.763 2.283 ( -0.700 1.212 -0.362) 1.445 2.291 ( 0.441 -0.763 -0.983) 1.320 2.338 ( 1.759 -3.047 -1.702) 3.908 2.737 ( 0.312 -0.540 1.787) 1.893 2.767 ( 0.585 -1.013 1.140) 1.633 4.240 ( 0.085 -0.148 0.392) 0.428 4.269 ( 0.204 -0.354 0.332) 0.526 4.270 ( 0.884 -1.531 -0.109) 1.772 4.299 ( -0.810 1.403 -0.425) 1.674 4.348 ( -1.459 2.527 0.276) 2.931 4.385 ( -0.956 1.656 -0.522) 1.982 4.889 ( -0.041 0.071 0.112) 0.139 4.909 ( -0.009 0.015 -0.289) 0.290 5.027 ( -1.284 2.224 0.065) 2.568 5.039 ( -0.913 1.581 -0.205) 1.837 5.176 ( -0.619 1.072 0.256) 1.265 5.224 ( -0.297 0.514 -0.812) 1.006 5.326 ( 2.042 -3.537 0.247) 4.091 5.336 ( 1.761 -3.050 -0.106) 3.523 5.380 ( 0.916 -1.587 -0.009) 1.833 5.398 ( 1.156 -2.002 -0.309) 2.332 5.534 ( -1.230 2.130 0.076) 2.461 5.552 ( -0.890 1.542 -0.273) 1.801 6.825 ( -0.900 1.560 0.206) 1.813 6.846 ( -0.715 1.239 -0.209) 1.445 7.076 ( -0.333 0.576 0.742) 0.996 7.138 ( 0.509 -0.881 -0.586) 1.174 7.690 ( 0.604 -1.046 0.267) 1.237 7.702 ( 0.832 -1.441 0.004) 1.664 ======================= Grid point 72 (16/24) ======================= q-point: ( 0.29 0.29 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.543 ( -0.538 -0.933 1.481) 1.831 0.563 ( 1.009 1.747 4.017) 4.495 0.617 ( -1.002 -1.735 2.341) 3.082 0.630 ( 0.129 0.224 -0.256) 0.363 0.772 ( -0.184 -0.318 5.290) 5.303 0.837 ( 0.614 1.064 -0.831) 1.483 0.883 ( -1.516 -2.626 -0.322) 3.049 1.073 ( 1.911 3.310 -0.763) 3.898 1.083 ( 1.086 1.882 -0.564) 2.245 1.131 ( 0.442 0.766 -2.274) 2.440 1.306 ( 2.767 4.793 0.021) 5.534 1.357 ( 0.730 1.264 -0.398) 1.513 1.375 ( -0.134 -0.232 -0.591) 0.649 1.380 ( 0.017 0.029 0.228) 0.231 1.524 ( 1.321 2.288 -0.791) 2.758 1.550 ( 2.459 4.259 -0.941) 5.007 1.660 ( -1.563 -2.708 1.572) 3.500 1.727 ( -2.732 -4.732 -0.298) 5.472 1.873 ( 2.222 3.849 -0.113) 4.446 1.900 ( 1.488 2.578 0.035) 2.977 2.019 ( -2.278 -3.945 0.864) 4.637 2.093 ( 1.467 2.542 -1.757) 3.420 2.108 ( 1.061 1.838 0.866) 2.292 2.163 ( 0.518 0.897 -0.492) 1.147 2.210 ( -1.684 -2.916 -1.318) 3.616 2.267 ( -0.054 -0.094 -0.581) 0.591 2.295 ( 0.216 0.374 -0.075) 0.438 2.313 ( -1.943 -3.366 -0.678) 3.945 2.729 ( -0.032 -0.056 1.817) 1.818 2.759 ( -1.477 -2.558 1.110) 3.155 4.234 ( 0.828 1.435 0.295) 1.683 4.245 ( 0.323 0.560 0.081) 0.652 4.274 ( -0.946 -1.638 0.414) 1.936 4.323 ( -0.207 -0.358 -0.660) 0.779 4.383 ( 0.723 1.252 0.045) 1.447 4.400 ( 0.219 0.380 -0.324) 0.545 4.894 ( 0.136 0.235 -0.054) 0.277 4.904 ( -0.103 -0.179 -0.257) 0.330 5.060 ( 1.936 3.354 0.058) 3.873 5.071 ( 1.552 2.689 -0.163) 3.109 5.179 ( 0.471 0.815 0.336) 1.000 5.222 ( 0.782 1.354 -0.571) 1.664 5.288 ( -0.466 -0.807 0.072) 0.935 5.289 ( -0.893 -1.547 0.025) 1.787 5.365 ( -2.411 -4.176 -0.217) 4.827 5.369 ( -1.337 -2.316 -0.263) 2.687 5.559 ( 0.881 1.526 0.064) 1.764 5.570 ( 0.363 0.629 -0.189) 0.750 6.850 ( 0.616 1.067 0.222) 1.252 6.861 ( 0.227 0.392 0.002) 0.453 7.085 ( 1.169 2.024 0.422) 2.375 7.120 ( -0.239 -0.414 -0.341) 0.587 7.677 ( -0.839 -1.454 0.203) 1.691 7.685 ( -1.193 -2.066 0.044) 2.385 ======================= Grid point 110 (17/24) ======================= q-point: ( 0.00 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.424 ( -0.000 0.000 5.600) 5.600 0.424 ( 0.000 0.000 5.600) 5.600 0.484 ( 0.000 -0.000 -0.716) 0.716 0.484 ( 0.000 0.000 -0.716) 0.716 0.693 ( 0.000 0.000 14.540) 14.540 0.704 ( 0.000 0.000 4.335) 4.335 0.704 ( 0.000 -0.000 4.335) 4.335 0.790 ( 0.000 -0.000 0.694) 0.694 0.807 ( 0.000 0.000 -0.795) 0.795 0.828 ( 0.000 0.000 -6.365) 6.365 0.828 ( 0.000 0.000 -6.365) 6.365 0.966 ( 0.000 0.000 -10.601) 10.601 1.418 ( -0.000 -0.000 -0.247) 0.247 1.418 ( -0.000 -0.000 -0.247) 0.247 1.418 ( -0.000 -0.000 -0.318) 0.318 1.418 ( 0.000 0.000 -0.318) 0.318 1.728 ( -0.000 0.000 0.508) 0.508 1.728 ( 0.000 0.000 0.508) 0.508 1.734 ( 0.000 -0.000 -0.083) 0.083 1.734 ( 0.000 -0.000 -0.083) 0.083 2.009 ( 0.000 0.000 4.357) 4.357 2.030 ( 0.000 -0.000 2.430) 2.430 2.267 ( -0.000 -0.000 -3.078) 3.078 2.268 ( 0.000 -0.000 -3.189) 3.189 2.301 ( -0.000 0.000 0.344) 0.344 2.301 ( 0.000 -0.000 0.344) 0.344 2.310 ( -0.000 0.000 -0.417) 0.417 2.310 ( -0.000 0.000 -0.417) 0.417 3.111 ( 0.000 -0.000 0.171) 0.171 3.119 ( -0.000 0.000 -0.505) 0.505 4.165 ( 0.000 0.000 0.380) 0.380 4.165 ( 0.000 0.000 0.380) 0.380 4.173 ( 0.000 -0.000 -0.328) 0.328 4.173 ( -0.000 0.000 -0.328) 0.328 4.313 ( 0.000 0.000 0.260) 0.260 4.320 ( 0.000 -0.000 -0.335) 0.335 4.799 ( 0.000 0.000 1.250) 1.250 4.824 ( 0.000 0.000 -1.008) 1.008 4.977 ( -0.000 0.000 0.408) 0.408 4.977 ( 0.000 0.000 0.408) 0.408 4.988 ( 0.000 0.000 -0.613) 0.613 4.988 ( 0.000 0.000 -0.613) 0.613 5.336 ( -0.000 -0.000 0.621) 0.621 5.336 ( 0.000 0.000 0.621) 0.621 5.352 ( 0.000 -0.000 -0.799) 0.799 5.352 ( 0.000 0.000 -0.799) 0.799 5.663 ( 0.000 0.000 1.841) 1.841 5.704 ( 0.000 0.000 -1.893) 1.893 6.831 ( 0.000 0.000 0.724) 0.724 6.831 ( 0.000 0.000 0.724) 0.724 6.847 ( 0.000 0.000 -0.753) 0.753 6.847 ( -0.000 0.000 -0.753) 0.753 7.908 ( 0.000 -0.000 1.539) 1.539 7.947 ( 0.000 -0.000 -1.941) 1.941 ======================= Grid point 111 (18/24) ======================= q-point: ( 0.14 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.503 ( 3.281 1.894 8.149) 8.987 0.505 ( 5.218 3.013 2.937) 6.703 0.529 ( 3.020 1.744 0.091) 3.488 0.589 ( 6.073 3.506 -0.500) 7.030 0.690 ( -0.715 -0.413 4.752) 4.823 0.708 ( 2.467 1.424 1.351) 3.153 0.799 ( -1.728 -0.998 -6.107) 6.425 0.802 ( -0.289 -0.167 -4.271) 4.284 0.866 ( 4.811 2.778 0.446) 5.573 0.881 ( 11.973 6.913 7.686) 15.818 0.892 ( 5.067 2.926 -2.137) 6.229 1.041 ( 5.714 3.299 -6.586) 9.323 1.403 ( -1.032 -0.596 0.289) 1.226 1.413 ( -0.624 -0.360 -0.664) 0.980 1.429 ( 0.440 0.254 1.106) 1.217 1.458 ( 2.391 1.380 -1.535) 3.159 1.727 ( 0.304 0.176 1.000) 1.060 1.740 ( 1.612 0.931 1.180) 2.205 1.746 ( 0.924 0.533 -0.749) 1.303 1.764 ( 2.383 1.376 -0.877) 2.888 1.992 ( -1.179 -0.681 5.628) 5.790 2.044 ( 0.962 0.556 1.385) 1.776 2.190 ( -2.831 -1.634 -1.307) 3.521 2.209 ( -3.396 -1.961 -3.226) 5.078 2.282 ( -1.354 -0.782 0.868) 1.789 2.301 ( -0.840 -0.485 -0.905) 1.327 2.320 ( -0.103 -0.060 0.603) 0.615 2.364 ( 1.453 0.839 -3.287) 3.690 2.985 ( -6.720 -3.880 3.608) 8.558 3.038 ( -5.327 -3.076 -1.638) 6.366 4.169 ( 0.395 0.228 0.312) 0.553 4.174 ( 0.342 0.197 -0.201) 0.443 4.260 ( 4.929 2.846 1.087) 5.794 4.299 ( 5.617 3.243 -2.786) 7.059 4.329 ( 0.647 0.374 -0.001) 0.747 4.329 ( 0.858 0.496 -0.053) 0.993 4.814 ( 1.734 1.001 0.887) 2.190 4.835 ( 1.200 0.693 -0.978) 1.696 4.977 ( 0.386 0.223 0.314) 0.545 4.986 ( 0.639 0.369 0.197) 0.763 4.989 ( 0.764 0.441 -0.830) 1.211 4.993 ( 0.355 0.205 -0.455) 0.613 5.346 ( 0.757 0.437 0.719) 1.132 5.362 ( 0.816 0.471 -0.806) 1.240 5.372 ( 2.455 1.417 0.865) 2.963 5.396 ( 2.706 1.562 -1.376) 3.414 5.621 ( -3.112 -1.797 1.502) 3.895 5.652 ( -3.625 -2.093 -1.212) 4.358 6.827 ( -0.300 -0.173 0.657) 0.742 6.842 ( -0.404 -0.233 -0.699) 0.840 6.885 ( 3.934 2.271 1.315) 4.729 6.923 ( 4.818 2.782 -2.364) 6.045 7.875 ( -2.565 -1.481 1.373) 3.265 7.908 ( -2.980 -1.720 -1.642) 3.812 ======================= Grid point 112 (19/24) ======================= q-point: ( 0.29 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.580 ( 0.647 0.374 1.216) 1.427 0.590 ( 3.973 2.294 6.758) 8.168 0.597 ( 0.983 0.568 -0.429) 1.213 0.694 ( 2.146 1.239 -1.589) 2.943 0.730 ( 4.781 2.761 4.173) 6.921 0.765 ( 1.720 0.993 -0.534) 2.056 0.788 ( 1.489 0.860 -3.301) 3.722 0.817 ( 2.078 1.200 -3.647) 4.366 0.986 ( 3.628 2.094 0.758) 4.257 1.007 ( 3.475 2.006 -1.052) 4.148 1.171 ( 8.405 4.852 0.604) 9.724 1.187 ( 5.638 3.255 -0.888) 6.571 1.367 ( -1.682 -0.971 1.118) 2.241 1.392 ( -1.008 -0.582 -1.097) 1.599 1.437 ( 0.871 0.503 2.783) 2.959 1.493 ( 0.457 0.264 -2.293) 2.353 1.752 ( 1.730 0.999 1.352) 2.412 1.783 ( 2.031 1.172 -1.401) 2.731 1.810 ( 2.047 1.182 1.092) 2.604 1.824 ( 0.997 0.575 0.373) 1.210 1.976 ( 1.534 0.886 3.386) 3.821 2.051 ( -0.721 -0.416 -2.412) 2.552 2.149 ( 0.134 0.077 0.038) 0.159 2.163 ( 0.392 0.226 -0.813) 0.930 2.243 ( -1.680 -0.970 1.223) 2.293 2.270 ( -1.584 -0.914 -1.232) 2.205 2.302 ( -1.161 -0.670 1.282) 1.855 2.363 ( -1.076 -0.621 -3.753) 3.953 2.855 ( -3.868 -2.233 3.248) 5.522 2.906 ( -4.915 -2.838 -1.532) 5.878 4.191 ( 1.535 0.886 0.321) 1.801 4.197 ( 1.528 0.882 -0.187) 1.774 4.334 ( 0.965 0.557 0.445) 1.200 4.335 ( -0.299 -0.173 0.307) 0.462 4.341 ( -0.258 -0.149 -0.283) 0.411 4.347 ( -0.665 -0.384 -0.695) 1.036 4.869 ( 1.881 1.086 0.239) 2.185 4.877 ( 1.480 0.854 -0.437) 1.764 5.000 ( 0.398 0.230 0.047) 0.462 5.004 ( 0.559 0.323 -0.286) 0.706 5.027 ( 4.036 2.330 0.536) 4.691 5.049 ( 4.269 2.465 -1.534) 5.163 5.371 ( 1.228 0.709 0.796) 1.627 5.388 ( 1.170 0.676 -0.741) 1.541 5.428 ( 1.397 0.807 0.936) 1.866 5.450 ( 1.010 0.583 -0.948) 1.503 5.532 ( -3.408 -1.968 0.860) 4.028 5.553 ( -3.907 -2.256 -1.005) 4.623 6.819 ( -0.338 -0.195 0.498) 0.633 6.830 ( -0.477 -0.276 -0.570) 0.793 6.996 ( 4.509 2.603 1.677) 5.470 7.038 ( 4.159 2.401 -2.255) 5.305 7.795 ( -3.560 -2.055 0.901) 4.208 7.815 ( -4.140 -2.390 -0.929) 4.870 ======================= Grid point 113 (20/24) ======================= q-point: ( 0.43 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.578 ( -0.189 -0.109 1.251) 1.270 0.598 ( -0.274 -0.158 -0.656) 0.728 0.672 ( 2.030 1.172 5.356) 5.847 0.729 ( 0.942 0.544 0.488) 1.192 0.792 ( 0.470 0.272 -0.347) 0.645 0.834 ( 1.422 0.821 -4.236) 4.543 0.857 ( 4.260 2.460 1.221) 5.069 0.884 ( 2.160 1.247 -1.306) 2.816 1.015 ( -0.673 -0.389 1.611) 1.788 1.049 ( 0.097 0.056 -1.457) 1.461 1.237 ( 0.399 0.230 2.859) 2.896 1.299 ( 1.275 0.736 -2.739) 3.109 1.336 ( -0.614 -0.354 1.750) 1.888 1.371 ( -0.554 -0.320 -1.367) 1.509 1.485 ( 2.391 1.381 2.220) 3.543 1.533 ( 2.502 1.444 -1.879) 3.446 1.751 ( -3.096 -1.788 2.319) 4.262 1.795 ( 1.251 0.722 1.588) 2.147 1.801 ( -2.948 -1.702 -1.993) 3.945 1.832 ( 1.435 0.829 -1.743) 2.405 2.055 ( 2.178 1.257 0.228) 2.525 2.067 ( 1.638 0.946 -1.149) 2.213 2.174 ( 1.523 0.879 2.627) 3.161 2.208 ( -0.880 -0.508 1.168) 1.548 2.216 ( 2.224 1.284 -0.965) 2.743 2.235 ( -0.953 -0.550 -1.217) 1.641 2.276 ( -0.346 -0.200 0.700) 0.806 2.328 ( -1.666 -0.962 -3.302) 3.822 2.791 ( -1.291 -0.745 1.792) 2.331 2.812 ( -2.395 -1.383 -0.179) 2.771 4.238 ( 1.859 1.073 0.433) 2.189 4.246 ( 2.256 1.302 -0.299) 2.622 4.312 ( -1.399 -0.808 0.101) 1.619 4.314 ( -1.708 -0.986 -0.073) 1.973 4.320 ( -0.859 -0.496 0.658) 1.191 4.336 ( -0.290 -0.167 -0.780) 0.848 4.893 ( 0.307 0.177 0.266) 0.444 4.899 ( 0.395 0.228 -0.291) 0.541 5.007 ( 0.127 0.073 0.312) 0.345 5.015 ( 0.228 0.132 -0.404) 0.482 5.142 ( 3.999 2.309 0.697) 4.670 5.165 ( 3.848 2.222 -1.435) 4.669 5.394 ( -2.636 -1.522 0.707) 3.125 5.397 ( 0.635 0.367 0.734) 1.038 5.409 ( -3.281 -1.894 -0.634) 3.841 5.412 ( 0.571 0.330 -0.625) 0.909 5.503 ( 0.327 0.189 0.406) 0.554 5.514 ( 0.163 0.094 -0.662) 0.688 6.812 ( -0.136 -0.079 0.359) 0.392 6.822 ( -0.205 -0.118 -0.458) 0.515 7.083 ( 2.137 1.234 1.442) 2.859 7.115 ( 1.780 1.028 -1.418) 2.498 7.722 ( -1.834 -1.059 0.464) 2.168 7.731 ( -2.124 -1.226 -0.308) 2.472 ======================= Grid point 118 (21/24) ======================= q-point: ( 0.14 0.14 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.544 ( 1.833 3.174 5.970) 7.005 0.588 ( 1.225 2.122 -0.196) 2.459 0.596 ( 1.435 2.485 1.711) 3.341 0.661 ( 0.915 1.585 -1.098) 2.134 0.744 ( 2.117 3.667 2.446) 4.890 0.745 ( 1.240 2.148 -0.515) 2.533 0.782 ( 1.716 2.972 -0.246) 3.441 0.865 ( 2.080 3.602 -3.530) 5.455 0.955 ( 3.113 5.391 -0.908) 6.291 0.957 ( 2.656 4.601 -0.480) 5.334 1.098 ( 6.028 10.441 1.577) 12.159 1.136 ( 2.782 4.819 -1.973) 5.903 1.386 ( -0.553 -0.957 0.577) 1.246 1.399 ( -0.459 -0.795 -0.673) 1.138 1.429 ( -0.008 -0.014 2.567) 2.567 1.484 ( 0.179 0.310 -2.382) 2.408 1.748 ( 0.886 1.534 0.566) 1.859 1.761 ( 0.816 1.414 -0.559) 1.726 1.789 ( 2.041 3.536 1.230) 4.264 1.815 ( 1.346 2.332 -0.433) 2.727 1.972 ( -0.347 -0.602 4.821) 4.871 2.053 ( -0.101 -0.174 -0.954) 0.975 2.149 ( -1.405 -2.434 -1.425) 3.151 2.152 ( -0.721 -1.248 -1.156) 1.847 2.277 ( -0.320 -0.554 0.264) 0.692 2.279 ( -0.584 -1.012 -0.177) 1.182 2.302 ( -0.854 -1.480 1.144) 2.056 2.366 ( -0.512 -0.887 -3.878) 4.011 2.882 ( -2.698 -4.672 3.622) 6.498 2.939 ( -3.093 -5.357 -1.763) 6.432 4.183 ( 0.583 1.010 0.210) 1.185 4.186 ( 0.620 1.075 -0.062) 1.242 4.325 ( 0.969 1.678 0.942) 2.155 4.341 ( 0.595 1.030 0.235) 1.213 4.342 ( 0.404 0.699 -0.348) 0.879 4.356 ( 0.408 0.707 -1.379) 1.603 4.855 ( 1.426 2.470 0.116) 2.855 4.863 ( 1.084 1.877 -0.548) 2.236 4.965 ( 0.131 0.227 0.160) 0.308 4.971 ( 0.029 0.050 -0.325) 0.330 5.039 ( 1.981 3.431 0.393) 3.981 5.056 ( 2.162 3.744 -1.285) 4.511 5.348 ( -0.072 -0.125 0.748) 0.762 5.363 ( -0.111 -0.193 -0.631) 0.669 5.429 ( 1.694 2.934 0.809) 3.483 5.450 ( 1.506 2.608 -1.117) 3.212 5.556 ( -2.235 -3.871 1.097) 4.602 5.579 ( -2.421 -4.194 -1.023) 4.950 6.826 ( 0.046 0.079 0.549) 0.557 6.838 ( -0.043 -0.075 -0.584) 0.590 6.962 ( 2.556 4.427 1.612) 5.360 7.005 ( 2.435 4.217 -2.437) 5.445 7.818 ( -2.091 -3.621 1.038) 4.308 7.842 ( -2.431 -4.211 -1.132) 4.992 ======================= Grid point 119 (22/24) ======================= q-point: ( 0.29 0.14 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 123 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.592 ( -0.794 0.400 1.424) 1.679 0.608 ( 0.225 0.737 0.208) 0.798 0.636 ( 2.044 1.410 4.243) 4.916 0.684 ( 2.238 -0.623 0.427) 2.362 0.784 ( 0.811 1.103 -0.781) 1.576 0.796 ( 2.035 1.021 -1.668) 2.823 0.889 ( 1.025 7.330 2.012) 7.670 0.938 ( -0.140 6.680 -2.261) 7.053 1.023 ( -0.124 1.515 0.670) 1.661 1.049 ( 1.318 1.970 -1.572) 2.844 1.232 ( 2.391 2.815 1.584) 4.019 1.271 ( 3.515 3.578 -1.836) 5.341 1.369 ( -1.260 1.010 0.838) 1.819 1.385 ( -0.478 -0.051 -0.490) 0.686 1.455 ( 1.956 1.774 1.492) 3.033 1.490 ( 2.138 -0.095 -1.597) 2.670 1.771 ( -1.346 -1.898 0.874) 2.486 1.786 ( 0.142 -1.428 -0.446) 1.502 1.811 ( -0.587 1.599 2.065) 2.677 1.848 ( 0.511 1.637 -0.885) 1.929 2.015 ( 3.257 2.264 0.600) 4.012 2.044 ( 1.036 -0.026 -1.752) 2.036 2.138 ( 0.850 -0.766 1.045) 1.550 2.155 ( 2.059 -0.673 -0.443) 2.211 2.244 ( -1.421 -1.218 0.481) 1.933 2.252 ( -0.945 0.047 -0.378) 1.019 2.287 ( -0.507 -0.012 0.766) 0.918 2.334 ( -1.139 -1.956 -3.116) 3.851 2.803 ( -1.724 -2.370 2.284) 3.716 2.835 ( -2.672 -3.413 -0.706) 4.392 4.214 ( 1.581 1.144 0.346) 1.982 4.220 ( 1.710 1.265 -0.230) 2.140 4.317 ( -1.242 -1.167 0.559) 1.794 4.329 ( -1.767 -1.025 -0.523) 2.109 4.357 ( -0.678 1.307 0.290) 1.500 4.365 ( -0.224 1.667 -0.405) 1.730 4.888 ( 0.365 0.533 0.161) 0.666 4.894 ( 0.534 0.558 -0.326) 0.838 5.001 ( 1.791 1.522 0.018) 2.350 5.004 ( 1.966 1.528 -0.241) 2.502 5.120 ( 2.165 4.479 0.567) 5.007 5.139 ( 2.296 4.191 -1.224) 4.933 5.344 ( 1.227 -2.583 0.623) 2.926 5.355 ( 0.899 -3.034 -0.330) 3.182 5.432 ( -2.359 -2.585 0.588) 3.549 5.446 ( -2.802 -3.069 -0.595) 4.198 5.518 ( 0.312 1.936 0.427) 2.007 5.531 ( 0.009 1.476 -0.760) 1.660 6.829 ( -0.446 1.144 0.410) 1.295 6.839 ( -0.502 0.976 -0.425) 1.176 7.053 ( 2.766 2.553 1.387) 4.012 7.085 ( 2.440 1.816 -1.552) 3.415 7.737 ( -2.161 -3.076 0.554) 3.800 7.748 ( -2.504 -3.535 -0.458) 4.356 ======================= Grid point 120 (23/24) ======================= q-point: ( 0.43 0.14 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.583 ( -0.134 0.232 1.563) 1.586 0.608 ( -0.404 0.700 -0.757) 1.107 0.666 ( 0.473 -0.819 5.385) 5.467 0.711 ( 0.959 -1.661 1.354) 2.348 0.796 ( -0.162 0.280 -0.672) 0.746 0.829 ( 0.415 -0.719 -3.724) 3.815 0.945 ( -0.622 1.077 2.511) 2.802 0.976 ( -3.200 5.543 -0.773) 6.447 1.019 ( -1.685 2.918 0.990) 3.512 1.059 ( -0.621 1.076 -2.304) 2.618 1.266 ( -1.550 2.685 1.881) 3.626 1.314 ( -0.617 1.069 -2.328) 2.635 1.360 ( -1.071 1.854 0.938) 2.338 1.377 ( -0.216 0.375 -0.459) 0.631 1.506 ( -0.348 0.602 1.329) 1.500 1.538 ( 0.380 -0.659 -1.449) 1.637 1.713 ( 0.863 -1.495 1.482) 2.276 1.746 ( 1.401 -2.426 -1.418) 3.140 1.829 ( -1.364 2.362 1.546) 3.135 1.869 ( -1.169 2.024 -1.889) 3.005 2.063 ( 0.431 -0.746 0.166) 0.877 2.075 ( 0.255 -0.442 -1.206) 1.310 2.136 ( 1.074 -1.860 2.421) 3.236 2.174 ( 1.181 -2.045 -1.280) 2.686 2.226 ( 0.471 -0.816 1.922) 2.141 2.275 ( -0.534 0.925 -0.418) 1.147 2.276 ( -0.352 0.610 -0.330) 0.778 2.286 ( 0.993 -1.720 -2.881) 3.499 2.775 ( 0.385 -0.667 1.344) 1.549 2.786 ( 0.494 -0.856 0.315) 1.037 4.251 ( 0.178 -0.309 0.635) 0.728 4.264 ( 0.519 -0.898 -0.471) 1.139 4.278 ( 0.068 -0.118 0.480) 0.499 4.288 ( -0.353 0.611 -0.520) 0.876 4.357 ( -1.287 2.228 0.546) 2.630 4.371 ( -1.097 1.901 -0.703) 2.305 4.893 ( 0.026 -0.045 0.302) 0.307 4.901 ( 0.046 -0.079 -0.413) 0.423 5.029 ( -1.174 2.034 0.160) 2.354 5.034 ( -1.032 1.788 -0.245) 2.079 5.185 ( -0.689 1.193 0.641) 1.519 5.203 ( -0.581 1.006 -0.997) 1.531 5.331 ( 2.060 -3.568 0.169) 4.124 5.334 ( 1.921 -3.328 -0.087) 3.843 5.382 ( 1.007 -1.745 0.233) 2.028 5.390 ( 1.115 -1.931 -0.440) 2.273 5.537 ( -1.170 2.026 0.271) 2.355 5.545 ( -1.032 1.787 -0.369) 2.097 6.832 ( -0.921 1.596 0.347) 1.875 6.839 ( -0.854 1.479 -0.351) 1.743 7.097 ( -0.062 0.107 1.107) 1.114 7.120 ( 0.252 -0.437 -1.007) 1.126 7.696 ( 0.683 -1.182 0.307) 1.399 7.701 ( 0.768 -1.330 -0.136) 1.542 ======================= Grid point 127 (24/24) ======================= q-point: ( 0.29 0.29 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 6.53e-05 0.00e+00 0.00e+00 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 1.56e-04 1.56e-04 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.583 ( -0.351 -0.609 1.930) 2.054 0.616 ( -0.033 -0.057 -1.049) 1.051 0.664 ( 0.579 1.004 4.377) 4.528 0.681 ( 0.039 0.067 2.764) 2.765 0.805 ( 0.369 0.640 -1.802) 1.947 0.812 ( 0.348 0.603 -1.978) 2.097 0.937 ( -0.737 -1.276 5.332) 5.531 1.037 ( 1.365 2.364 -3.314) 4.294 1.073 ( 1.247 2.161 -0.157) 2.500 1.083 ( 0.544 0.943 -1.383) 1.760 1.313 ( 2.026 3.510 0.737) 4.119 1.337 ( 1.255 2.173 -1.324) 2.837 1.372 ( -0.313 -0.541 0.368) 0.725 1.381 ( -0.539 -0.934 -0.217) 1.100 1.514 ( 1.787 3.095 0.087) 3.575 1.527 ( 1.822 3.156 -1.005) 3.780 1.695 ( -1.881 -3.259 1.470) 4.040 1.718 ( -2.272 -3.935 -0.655) 4.590 1.876 ( 1.609 2.786 0.533) 3.261 1.892 ( 1.168 2.023 -0.738) 2.450 2.035 ( -0.616 -1.067 0.671) 1.403 2.056 ( 0.724 1.255 -1.430) 2.036 2.135 ( 0.508 0.880 1.081) 1.484 2.149 ( 0.232 0.402 -0.637) 0.788 2.198 ( -1.397 -2.420 0.984) 2.963 2.241 ( -0.820 -1.420 -2.027) 2.607 2.294 ( 0.142 0.245 -0.003) 0.283 2.298 ( -0.481 -0.832 -0.467) 1.069 2.767 ( -0.449 -0.777 1.341) 1.614 2.778 ( -0.976 -1.691 0.326) 1.980 4.241 ( 0.768 1.329 0.327) 1.569 4.246 ( 0.541 0.937 -0.067) 1.084 4.287 ( -0.843 -1.460 0.753) 1.846 4.305 ( -0.529 -0.916 -0.927) 1.406 4.385 ( 0.545 0.944 0.198) 1.108 4.392 ( 0.355 0.615 -0.378) 0.804 4.894 ( 0.057 0.099 0.070) 0.134 4.898 ( -0.036 -0.063 -0.266) 0.276 5.063 ( 1.914 3.316 0.146) 3.832 5.067 ( 1.772 3.070 -0.210) 3.551 5.190 ( 0.742 1.285 0.669) 1.628 5.207 ( 0.919 1.592 -0.831) 2.017 5.290 ( -1.027 -1.778 0.042) 2.054 5.290 ( -0.864 -1.496 0.064) 1.728 5.361 ( -2.101 -3.639 -0.084) 4.202 5.363 ( -1.627 -2.818 -0.249) 3.264 5.561 ( 0.769 1.332 0.156) 1.546 5.565 ( 0.582 1.008 -0.222) 1.185 6.856 ( 0.547 0.948 0.264) 1.126 6.860 ( 0.395 0.684 -0.129) 0.801 7.097 ( 0.775 1.342 0.638) 1.676 7.110 ( 0.245 0.425 -0.582) 0.761 7.682 ( -0.941 -1.629 0.220) 1.894 7.685 ( -1.073 -1.859 -0.066) 2.148 =================== 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/13230 10.0 14.816 14.816 8.909 0.000 -0.000 0.000 3/13230 20.0 5.096 5.096 2.577 0.000 -0.000 0.000 3/13230 30.0 2.868 2.868 1.392 0.000 -0.000 0.000 3/13230 40.0 1.974 1.974 0.947 0.000 -0.000 0.000 3/13230 50.0 1.504 1.504 0.720 0.000 -0.000 0.000 3/13230 60.0 1.217 1.217 0.582 0.000 -0.000 0.000 3/13230 70.0 1.025 1.025 0.489 0.000 -0.000 0.000 3/13230 80.0 0.887 0.887 0.423 0.000 -0.000 0.000 3/13230 90.0 0.784 0.784 0.373 0.000 -0.000 0.000 3/13230 100.0 0.703 0.703 0.334 0.000 -0.000 0.000 3/13230 110.0 0.638 0.638 0.303 0.000 -0.000 0.000 3/13230 120.0 0.584 0.584 0.277 0.000 -0.000 0.000 3/13230 130.0 0.539 0.539 0.255 0.000 -0.000 0.000 3/13230 140.0 0.501 0.501 0.236 0.000 -0.000 0.000 3/13230 150.0 0.468 0.468 0.220 0.000 -0.000 0.000 3/13230 160.0 0.439 0.439 0.207 0.000 -0.000 0.000 3/13230 170.0 0.413 0.413 0.194 0.000 -0.000 0.000 3/13230 180.0 0.391 0.391 0.183 0.000 -0.000 0.000 3/13230 190.0 0.370 0.370 0.174 0.000 -0.000 0.000 3/13230 200.0 0.352 0.352 0.165 0.000 -0.000 0.000 3/13230 210.0 0.335 0.335 0.157 0.000 -0.000 0.000 3/13230 220.0 0.320 0.320 0.150 0.000 -0.000 0.000 3/13230 230.0 0.307 0.307 0.144 0.000 -0.000 0.000 3/13230 240.0 0.294 0.294 0.138 0.000 -0.000 0.000 3/13230 250.0 0.283 0.283 0.132 0.000 -0.000 0.000 3/13230 260.0 0.272 0.272 0.127 0.000 -0.000 0.000 3/13230 270.0 0.262 0.262 0.122 0.000 -0.000 0.000 3/13230 280.0 0.253 0.253 0.118 0.000 -0.000 0.000 3/13230 290.0 0.244 0.244 0.114 0.000 -0.000 0.000 3/13230 300.0 0.236 0.236 0.110 0.000 -0.000 0.000 3/13230 310.0 0.229 0.229 0.107 0.000 -0.000 0.000 3/13230 320.0 0.222 0.222 0.103 0.000 -0.000 0.000 3/13230 330.0 0.215 0.215 0.100 0.000 -0.000 0.000 3/13230 340.0 0.209 0.209 0.097 0.000 -0.000 0.000 3/13230 350.0 0.203 0.203 0.094 0.000 -0.000 0.000 3/13230 360.0 0.197 0.197 0.092 0.000 -0.000 0.000 3/13230 370.0 0.192 0.192 0.089 0.000 -0.000 0.000 3/13230 380.0 0.187 0.187 0.087 0.000 -0.000 0.000 3/13230 390.0 0.182 0.182 0.085 0.000 -0.000 0.000 3/13230 400.0 0.178 0.178 0.083 0.000 -0.000 0.000 3/13230 410.0 0.173 0.173 0.081 0.000 -0.000 0.000 3/13230 420.0 0.169 0.169 0.079 0.000 -0.000 0.000 3/13230 430.0 0.165 0.165 0.077 0.000 -0.000 0.000 3/13230 440.0 0.162 0.162 0.075 0.000 -0.000 0.000 3/13230 450.0 0.158 0.158 0.074 0.000 -0.000 0.000 3/13230 460.0 0.155 0.155 0.072 0.000 -0.000 0.000 3/13230 470.0 0.151 0.151 0.070 0.000 -0.000 0.000 3/13230 480.0 0.148 0.148 0.069 0.000 -0.000 0.000 3/13230 490.0 0.145 0.145 0.068 0.000 -0.000 0.000 3/13230 500.0 0.142 0.142 0.066 0.000 -0.000 0.000 3/13230 510.0 0.140 0.140 0.065 0.000 -0.000 0.000 3/13230 520.0 0.137 0.137 0.064 0.000 -0.000 0.000 3/13230 530.0 0.134 0.134 0.062 0.000 -0.000 0.000 3/13230 540.0 0.132 0.132 0.061 0.000 -0.000 0.000 3/13230 550.0 0.130 0.130 0.060 0.000 -0.000 0.000 3/13230 560.0 0.127 0.127 0.059 0.000 -0.000 0.000 3/13230 570.0 0.125 0.125 0.058 0.000 -0.000 0.000 3/13230 580.0 0.123 0.123 0.057 0.000 -0.000 0.000 3/13230 590.0 0.121 0.121 0.056 0.000 -0.000 0.000 3/13230 600.0 0.119 0.119 0.055 0.000 -0.000 0.000 3/13230 610.0 0.117 0.117 0.054 0.000 -0.000 0.000 3/13230 620.0 0.115 0.115 0.053 0.000 -0.000 0.000 3/13230 630.0 0.113 0.113 0.053 0.000 -0.000 0.000 3/13230 640.0 0.112 0.112 0.052 0.000 -0.000 0.000 3/13230 650.0 0.110 0.110 0.051 0.000 -0.000 0.000 3/13230 660.0 0.108 0.108 0.050 0.000 -0.000 0.000 3/13230 670.0 0.107 0.107 0.049 0.000 -0.000 0.000 3/13230 680.0 0.105 0.105 0.049 0.000 -0.000 0.000 3/13230 690.0 0.103 0.103 0.048 0.000 -0.000 0.000 3/13230 700.0 0.102 0.102 0.047 0.000 -0.000 0.000 3/13230 710.0 0.101 0.101 0.047 0.000 -0.000 0.000 3/13230 720.0 0.099 0.099 0.046 0.000 -0.000 0.000 3/13230 730.0 0.098 0.098 0.045 0.000 -0.000 0.000 3/13230 740.0 0.097 0.097 0.045 0.000 -0.000 0.000 3/13230 750.0 0.095 0.095 0.044 0.000 -0.000 0.000 3/13230 760.0 0.094 0.094 0.044 0.000 -0.000 0.000 3/13230 770.0 0.093 0.093 0.043 0.000 -0.000 0.000 3/13230 780.0 0.092 0.092 0.042 0.000 -0.000 0.000 3/13230 790.0 0.090 0.090 0.042 0.000 -0.000 0.000 3/13230 800.0 0.089 0.089 0.041 0.000 -0.000 0.000 3/13230 810.0 0.088 0.088 0.041 0.000 -0.000 0.000 3/13230 820.0 0.087 0.087 0.040 0.000 -0.000 0.000 3/13230 830.0 0.086 0.086 0.040 0.000 -0.000 0.000 3/13230 840.0 0.085 0.085 0.039 0.000 -0.000 0.000 3/13230 850.0 0.084 0.084 0.039 0.000 -0.000 0.000 3/13230 860.0 0.083 0.083 0.039 0.000 -0.000 0.000 3/13230 870.0 0.082 0.082 0.038 0.000 -0.000 0.000 3/13230 880.0 0.081 0.081 0.038 0.000 -0.000 0.000 3/13230 890.0 0.080 0.080 0.037 0.000 -0.000 0.000 3/13230 900.0 0.079 0.079 0.037 0.000 -0.000 0.000 3/13230 910.0 0.079 0.079 0.036 0.000 -0.000 0.000 3/13230 920.0 0.078 0.078 0.036 0.000 -0.000 0.000 3/13230 930.0 0.077 0.077 0.036 0.000 -0.000 0.000 3/13230 940.0 0.076 0.076 0.035 0.000 -0.000 0.000 3/13230 950.0 0.075 0.075 0.035 0.000 -0.000 0.000 3/13230 960.0 0.075 0.075 0.035 0.000 -0.000 0.000 3/13230 970.0 0.074 0.074 0.034 0.000 -0.000 0.000 3/13230 980.0 0.073 0.073 0.034 0.000 -0.000 0.000 3/13230 990.0 0.072 0.072 0.033 0.000 -0.000 0.000 3/13230 1000.0 0.072 0.072 0.033 0.000 -0.000 0.000 3/13230 Thermal conductivity related properties were written into "kappa-m775.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 17:52:50]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|