----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 17:55:33]------------------------- Compiled with OpenMP support (max 48 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [2 2 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pnna (52) Number of symmetry operations in supercell: 32 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 7.824470680000000 0.000000000000000 0.000000000000000 b 0.000000000000000 10.136594210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.287178179999998 Atomic positions (fractional): *1 In 0.25000000000000 0.00000000000000 0.32046657352454 114.818 *2 In 0.82856942391417 0.75000000000000 0.25000000000000 114.818 3 In 0.75000000000000 0.00000000000000 0.67953342647546 114.818 4 In 0.17143057608583 0.25000000000000 0.75000000000000 114.818 5 In 0.67143057608583 0.25000000000000 0.25000000000000 114.818 6 In 0.75000000000000 0.50000000000000 0.82046657352454 114.818 7 In 0.32856942391417 0.75000000000000 0.75000000000000 114.818 8 In 0.25000000000000 0.50000000000000 0.17953342647546 114.818 *9 Br 0.99262223094321 0.45258576902811 0.32037942549330 79.904 10 Br 0.00737776905679 0.54741423097189 0.67962057450670 79.904 11 Br 0.49262223094321 0.04741423097189 0.82037942549330 79.904 *12 Br 0.34647419754853 0.30613211532536 0.04422525730086 79.904 13 Br 0.00737776905679 0.95258576902811 0.82037942549330 79.904 14 Br 0.84647419754853 0.30613211532536 0.95577474269914 79.904 15 Br 0.34647419754853 0.19386788467464 0.45577474269914 79.904 16 Br 0.84647419754853 0.19386788467464 0.54422525730086 79.904 17 Br 0.15352580245147 0.80613211532536 0.45577474269914 79.904 18 Br 0.49262223094321 0.45258576902811 0.67962057450670 79.904 19 Br 0.65352580245147 0.80613211532536 0.54422525730086 79.904 20 Br 0.99262223094321 0.04741423097189 0.17962057450670 79.904 21 Br 0.15352580245147 0.69386788467464 0.04422525730086 79.904 22 Br 0.50737776905679 0.54741423097189 0.32037942549330 79.904 23 Br 0.65352580245147 0.69386788467464 0.95577474269914 79.904 24 Br 0.50737776905679 0.95258576902811 0.17962057450670 79.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 7.824470680000000 0.000000000000000 0.000000000000000 b 0.000000000000000 10.136594210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.287178179999998 Atomic positions (fractional): *1 In 0.25000000000000 0.00000000000000 0.32046657352454 114.818 > 1 *2 In 0.82856942391417 0.75000000000000 0.25000000000000 114.818 > 2 3 In 0.75000000000000 0.00000000000000 0.67953342647546 114.818 > 3 4 In 0.17143057608583 0.25000000000000 0.75000000000000 114.818 > 4 5 In 0.67143057608583 0.25000000000000 0.25000000000000 114.818 > 5 6 In 0.75000000000000 0.50000000000000 0.82046657352454 114.818 > 6 7 In 0.32856942391417 0.75000000000000 0.75000000000000 114.818 > 7 8 In 0.25000000000000 0.50000000000000 0.17953342647546 114.818 > 8 *9 Br 0.99262223094321 0.45258576902811 0.32037942549330 79.904 > 9 10 Br 0.00737776905679 0.54741423097189 0.67962057450670 79.904 > 10 11 Br 0.49262223094321 0.04741423097189 0.82037942549330 79.904 > 11 *12 Br 0.34647419754853 0.30613211532536 0.04422525730086 79.904 > 12 13 Br 0.00737776905679 0.95258576902811 0.82037942549330 79.904 > 13 14 Br 0.84647419754853 0.30613211532536 0.95577474269914 79.904 > 14 15 Br 0.34647419754853 0.19386788467464 0.45577474269914 79.904 > 15 16 Br 0.84647419754853 0.19386788467464 0.54422525730086 79.904 > 16 17 Br 0.15352580245147 0.80613211532536 0.45577474269914 79.904 > 17 18 Br 0.49262223094321 0.45258576902811 0.67962057450670 79.904 > 18 19 Br 0.65352580245147 0.80613211532536 0.54422525730086 79.904 > 19 20 Br 0.99262223094321 0.04741423097189 0.17962057450670 79.904 > 20 21 Br 0.15352580245147 0.69386788467464 0.04422525730086 79.904 > 21 22 Br 0.50737776905679 0.54741423097189 0.32037942549330 79.904 > 22 23 Br 0.65352580245147 0.69386788467464 0.95577474269914 79.904 > 23 24 Br 0.50737776905679 0.95258576902811 0.17962057450670 79.904 > 24 -------------------------------- super cell -------------------------------- Lattice vectors: a 15.648941360000000 0.000000000000000 0.000000000000000 b 0.000000000000000 20.273188420000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.287178179999998 Atomic positions (fractional): *1 In 0.12500000000000 0.00000000000000 0.32046657352454 114.818 > 1 2 In 0.62500000000000 0.00000000000000 0.32046657352454 114.818 > 1 3 In 0.12500000000000 0.50000000000000 0.32046657352454 114.818 > 1 4 In 0.62500000000000 0.50000000000000 0.32046657352454 114.818 > 1 *5 In 0.41428471195709 0.37500000000000 0.25000000000000 114.818 > 2 6 In 0.91428471195709 0.37500000000000 0.25000000000000 114.818 > 2 7 In 0.41428471195709 0.87500000000000 0.25000000000000 114.818 > 2 8 In 0.91428471195709 0.87500000000000 0.25000000000000 114.818 > 2 9 In 0.37500000000000 0.00000000000000 0.67953342647546 114.818 > 3 10 In 0.87500000000000 0.00000000000000 0.67953342647546 114.818 > 3 11 In 0.37500000000000 0.50000000000000 0.67953342647546 114.818 > 3 12 In 0.87500000000000 0.50000000000000 0.67953342647546 114.818 > 3 13 In 0.08571528804291 0.12500000000000 0.75000000000000 114.818 > 4 14 In 0.58571528804291 0.12500000000000 0.75000000000000 114.818 > 4 15 In 0.08571528804291 0.62500000000000 0.75000000000000 114.818 > 4 16 In 0.58571528804291 0.62500000000000 0.75000000000000 114.818 > 4 17 In 0.33571528804291 0.12500000000000 0.25000000000000 114.818 > 5 18 In 0.83571528804291 0.12500000000000 0.25000000000000 114.818 > 5 19 In 0.33571528804291 0.62500000000000 0.25000000000000 114.818 > 5 20 In 0.83571528804291 0.62500000000000 0.25000000000000 114.818 > 5 21 In 0.37500000000000 0.25000000000000 0.82046657352454 114.818 > 6 22 In 0.87500000000000 0.25000000000000 0.82046657352454 114.818 > 6 23 In 0.37500000000000 0.75000000000000 0.82046657352454 114.818 > 6 24 In 0.87500000000000 0.75000000000000 0.82046657352454 114.818 > 6 25 In 0.16428471195709 0.37500000000000 0.75000000000000 114.818 > 7 26 In 0.66428471195709 0.37500000000000 0.75000000000000 114.818 > 7 27 In 0.16428471195709 0.87500000000000 0.75000000000000 114.818 > 7 28 In 0.66428471195709 0.87500000000000 0.75000000000000 114.818 > 7 29 In 0.12500000000000 0.25000000000000 0.17953342647546 114.818 > 8 30 In 0.62500000000000 0.25000000000000 0.17953342647546 114.818 > 8 31 In 0.12500000000000 0.75000000000000 0.17953342647546 114.818 > 8 32 In 0.62500000000000 0.75000000000000 0.17953342647546 114.818 > 8 *33 Br 0.49631111547161 0.22629288451405 0.32037942549330 79.904 > 9 34 Br 0.99631111547161 0.22629288451405 0.32037942549330 79.904 > 9 35 Br 0.49631111547161 0.72629288451405 0.32037942549330 79.904 > 9 36 Br 0.99631111547161 0.72629288451405 0.32037942549330 79.904 > 9 37 Br 0.00368888452839 0.27370711548595 0.67962057450670 79.904 > 10 38 Br 0.50368888452839 0.27370711548595 0.67962057450670 79.904 > 10 39 Br 0.00368888452839 0.77370711548595 0.67962057450670 79.904 > 10 40 Br 0.50368888452839 0.77370711548595 0.67962057450670 79.904 > 10 41 Br 0.24631111547161 0.02370711548595 0.82037942549330 79.904 > 11 42 Br 0.74631111547161 0.02370711548595 0.82037942549330 79.904 > 11 43 Br 0.24631111547161 0.52370711548595 0.82037942549330 79.904 > 11 44 Br 0.74631111547161 0.52370711548595 0.82037942549330 79.904 > 11 *45 Br 0.17323709877426 0.15306605766268 0.04422525730086 79.904 > 12 46 Br 0.67323709877426 0.15306605766268 0.04422525730086 79.904 > 12 47 Br 0.17323709877426 0.65306605766268 0.04422525730086 79.904 > 12 48 Br 0.67323709877426 0.65306605766268 0.04422525730086 79.904 > 12 49 Br 0.00368888452839 0.47629288451405 0.82037942549330 79.904 > 13 50 Br 0.50368888452839 0.47629288451405 0.82037942549330 79.904 > 13 51 Br 0.00368888452839 0.97629288451405 0.82037942549330 79.904 > 13 52 Br 0.50368888452839 0.97629288451405 0.82037942549330 79.904 > 13 53 Br 0.42323709877426 0.15306605766268 0.95577474269914 79.904 > 14 54 Br 0.92323709877426 0.15306605766268 0.95577474269914 79.904 > 14 55 Br 0.42323709877426 0.65306605766268 0.95577474269914 79.904 > 14 56 Br 0.92323709877426 0.65306605766268 0.95577474269914 79.904 > 14 57 Br 0.17323709877426 0.09693394233732 0.45577474269914 79.904 > 15 58 Br 0.67323709877426 0.09693394233732 0.45577474269914 79.904 > 15 59 Br 0.17323709877426 0.59693394233732 0.45577474269914 79.904 > 15 60 Br 0.67323709877426 0.59693394233732 0.45577474269914 79.904 > 15 61 Br 0.42323709877426 0.09693394233732 0.54422525730086 79.904 > 16 62 Br 0.92323709877426 0.09693394233732 0.54422525730086 79.904 > 16 63 Br 0.42323709877426 0.59693394233732 0.54422525730086 79.904 > 16 64 Br 0.92323709877426 0.59693394233732 0.54422525730086 79.904 > 16 65 Br 0.07676290122574 0.40306605766268 0.45577474269914 79.904 > 17 66 Br 0.57676290122574 0.40306605766268 0.45577474269914 79.904 > 17 67 Br 0.07676290122574 0.90306605766268 0.45577474269914 79.904 > 17 68 Br 0.57676290122574 0.90306605766268 0.45577474269914 79.904 > 17 69 Br 0.24631111547161 0.22629288451405 0.67962057450670 79.904 > 18 70 Br 0.74631111547161 0.22629288451405 0.67962057450670 79.904 > 18 71 Br 0.24631111547161 0.72629288451405 0.67962057450670 79.904 > 18 72 Br 0.74631111547161 0.72629288451405 0.67962057450670 79.904 > 18 73 Br 0.32676290122574 0.40306605766268 0.54422525730086 79.904 > 19 74 Br 0.82676290122574 0.40306605766268 0.54422525730086 79.904 > 19 75 Br 0.32676290122574 0.90306605766268 0.54422525730086 79.904 > 19 76 Br 0.82676290122574 0.90306605766268 0.54422525730086 79.904 > 19 77 Br 0.49631111547161 0.02370711548595 0.17962057450670 79.904 > 20 78 Br 0.99631111547161 0.02370711548595 0.17962057450670 79.904 > 20 79 Br 0.49631111547161 0.52370711548595 0.17962057450670 79.904 > 20 80 Br 0.99631111547161 0.52370711548595 0.17962057450670 79.904 > 20 81 Br 0.07676290122574 0.34693394233732 0.04422525730086 79.904 > 21 82 Br 0.57676290122574 0.34693394233732 0.04422525730086 79.904 > 21 83 Br 0.07676290122574 0.84693394233732 0.04422525730086 79.904 > 21 84 Br 0.57676290122574 0.84693394233732 0.04422525730086 79.904 > 21 85 Br 0.25368888452839 0.27370711548595 0.32037942549330 79.904 > 22 86 Br 0.75368888452839 0.27370711548595 0.32037942549330 79.904 > 22 87 Br 0.25368888452839 0.77370711548595 0.32037942549330 79.904 > 22 88 Br 0.75368888452839 0.77370711548595 0.32037942549330 79.904 > 22 89 Br 0.32676290122574 0.34693394233732 0.95577474269914 79.904 > 23 90 Br 0.82676290122574 0.34693394233732 0.95577474269914 79.904 > 23 91 Br 0.32676290122574 0.84693394233732 0.95577474269914 79.904 > 23 92 Br 0.82676290122574 0.84693394233732 0.95577474269914 79.904 > 23 93 Br 0.25368888452839 0.47629288451405 0.17962057450670 79.904 > 24 94 Br 0.75368888452839 0.47629288451405 0.17962057450670 79.904 > 24 95 Br 0.25368888452839 0.97629288451405 0.17962057450670 79.904 > 24 96 Br 0.75368888452839 0.97629288451405 0.17962057450670 79.904 > 24 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 4.7068177 0.0000000 0.0000000 0.0000000 4.3610801 0.0000000 0.0000000 0.0000000 4.2926020 -------------------------- Born effective charges -------------------------- 1 In 3.2009179 0.3083122 -0.0000000 0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 2 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 0.2056455 -0.0000000 0.3485131 1.9142448 3 In 3.2009179 0.3083122 -0.0000000 0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 4 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 0.2056455 -0.0000000 0.3485131 1.9142448 5 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 -0.2056455 -0.0000000 -0.3485131 1.9142448 6 In 3.2009179 -0.3083122 -0.0000000 -0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 7 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 -0.2056455 -0.0000000 -0.3485131 1.9142448 8 In 3.2009179 -0.3083122 -0.0000000 -0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 9 Br -1.7506695 -0.1377095 0.5111143 -0.2532621 -1.2356452 0.2630911 0.4907325 0.1799490 -0.7513838 10 Br -1.7506695 -0.1377095 0.5111143 -0.2532621 -1.2356452 0.2630911 0.4907325 0.1799490 -0.7513838 11 Br -1.7506695 0.1377095 0.5111143 0.2532621 -1.2356452 -0.2630911 0.4907325 -0.1799490 -0.7513838 12 Br -1.0555407 0.2666999 -0.3476411 0.2653963 -1.1258863 -0.3852760 -0.2679068 -0.3901344 -1.3658623 13 Br -1.7506695 0.1377095 -0.5111143 0.2532621 -1.2356452 0.2630911 -0.4907325 0.1799490 -0.7513838 14 Br -1.0555407 0.2666999 0.3476411 0.2653963 -1.1258863 0.3852760 0.2679068 0.3901344 -1.3658623 15 Br -1.0555407 -0.2666999 0.3476411 -0.2653963 -1.1258863 -0.3852760 0.2679068 -0.3901344 -1.3658623 16 Br -1.0555407 -0.2666999 -0.3476411 -0.2653963 -1.1258863 0.3852760 -0.2679068 0.3901344 -1.3658623 17 Br -1.0555407 -0.2666999 -0.3476411 -0.2653963 -1.1258863 0.3852760 -0.2679068 0.3901344 -1.3658623 18 Br -1.7506695 -0.1377095 -0.5111143 -0.2532621 -1.2356452 -0.2630911 -0.4907325 -0.1799490 -0.7513838 19 Br -1.0555407 -0.2666999 0.3476411 -0.2653963 -1.1258863 -0.3852760 0.2679068 -0.3901344 -1.3658623 20 Br -1.7506695 0.1377095 -0.5111143 0.2532621 -1.2356452 0.2630911 -0.4907325 0.1799490 -0.7513838 21 Br -1.0555407 0.2666999 0.3476411 0.2653963 -1.1258863 0.3852760 0.2679068 0.3901344 -1.3658623 22 Br -1.7506695 -0.1377095 -0.5111143 -0.2532621 -1.2356452 -0.2630911 -0.4907325 -0.1799490 -0.7513838 23 Br -1.0555407 0.2666999 -0.3476411 0.2653963 -1.1258863 -0.3852760 -0.2679068 -0.3901344 -1.3658623 24 Br -1.7506695 0.1377095 0.5111143 0.2532621 -1.2356452 -0.2630911 0.4907325 -0.1799490 -0.7513838 ---------------------------------------------------------------------------- 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: 288/288 Permutation basis: 10656/10656 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 1371 Number of blocks in projector: 1371 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 6 --- Eigsh_solver_block: 1 / 6 --- Block_size: 290 Use standard eigh solver. --- Eigsh_solver_block: 2 / 6 --- Block_size: 286 Use standard eigh solver. --- Eigsh_solver_block: 3 / 6 --- Block_size: 282 Use standard eigh solver. --- Eigsh_solver_block: 4 / 6 --- Block_size: 171 Use standard eigh solver. --- Eigsh_solver_block: 5 / 6 --- Block_size: 171 Use standard eigh solver. --- Eigsh_solver_block: 6 / 6 --- Block_size: 171 Use standard eigh solver. Tree of FC basis block matrices: - (1371, 1343), data: False |-- (171, 167), data: True |-- (171, 167), data: True |-- (171, 167), data: True |-- (282, 278), data: True |-- (286, 280), data: True |-- (290, 284), data: True ----- Solver_atoms: 1 -- 96 / 96 Time (Solver_compr_matrix_reshape): 0.002 Solver_block: 100 / 400 - Time: 0.305 Solver_block: 200 / 400 - Time: 0.284 Solver_block: 300 / 400 - Time: 0.277 Solver_block: 400 / 400 - Time: 0.280 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 1.176 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) Permutation basis: 288/288 Permutation basis: 10656/10656 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 1371 Number of blocks in projector: 1371 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 6 --- Eigsh_solver_block: 1 / 6 --- Block_size: 290 Use standard eigh solver. --- Eigsh_solver_block: 2 / 6 --- Block_size: 286 Use standard eigh solver. --- Eigsh_solver_block: 3 / 6 --- Block_size: 282 Use standard eigh solver. --- Eigsh_solver_block: 4 / 6 --- Block_size: 171 Use standard eigh solver. --- Eigsh_solver_block: 5 / 6 --- Block_size: 171 Use standard eigh solver. --- Eigsh_solver_block: 6 / 6 --- Block_size: 171 Use standard eigh solver. Tree of FC basis block matrices: - (1371, 1343), data: False |-- (171, 167), data: True |-- (171, 167), data: True |-- (171, 167), data: True |-- (282, 278), data: True |-- (286, 280), data: True |-- (290, 284), data: True Max drift after symmetrization by symfc projector: -0.00000000 (zx) -0.00000000 (zx) 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:55:47]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 17:55:47]------------------------- Compiled with OpenMP support (max 48 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [2 2 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pnna (52) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 7.824470680000000 0.000000000000000 0.000000000000000 b 0.000000000000000 10.136594210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.287178179999998 Atomic positions (fractional): 1 In 0.25000000000000 0.00000000000000 0.32046657352454 114.818 2 In 0.82856942391417 0.75000000000000 0.25000000000000 114.818 3 In 0.75000000000000 0.00000000000000 0.67953342647546 114.818 4 In 0.17143057608583 0.25000000000000 0.75000000000000 114.818 5 In 0.67143057608583 0.25000000000000 0.25000000000000 114.818 6 In 0.75000000000000 0.50000000000000 0.82046657352454 114.818 7 In 0.32856942391417 0.75000000000000 0.75000000000000 114.818 8 In 0.25000000000000 0.50000000000000 0.17953342647546 114.818 9 Br 0.99262223094321 0.45258576902811 0.32037942549330 79.904 10 Br 0.00737776905679 0.54741423097189 0.67962057450670 79.904 11 Br 0.49262223094321 0.04741423097189 0.82037942549330 79.904 12 Br 0.34647419754853 0.30613211532536 0.04422525730086 79.904 13 Br 0.00737776905679 0.95258576902811 0.82037942549330 79.904 14 Br 0.84647419754853 0.30613211532536 0.95577474269914 79.904 15 Br 0.34647419754853 0.19386788467464 0.45577474269914 79.904 16 Br 0.84647419754853 0.19386788467464 0.54422525730086 79.904 17 Br 0.15352580245147 0.80613211532536 0.45577474269914 79.904 18 Br 0.49262223094321 0.45258576902811 0.67962057450670 79.904 19 Br 0.65352580245147 0.80613211532536 0.54422525730086 79.904 20 Br 0.99262223094321 0.04741423097189 0.17962057450670 79.904 21 Br 0.15352580245147 0.69386788467464 0.04422525730086 79.904 22 Br 0.50737776905679 0.54741423097189 0.32037942549330 79.904 23 Br 0.65352580245147 0.69386788467464 0.95577474269914 79.904 24 Br 0.50737776905679 0.95258576902811 0.17962057450670 79.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.648941360000000 0.000000000000000 0.000000000000000 b 0.000000000000000 20.273188420000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.287178179999998 Atomic positions (fractional): 1 In 0.12500000000000 0.00000000000000 0.32046657352454 114.818 > 1 2 In 0.62500000000000 0.00000000000000 0.32046657352454 114.818 > 1 3 In 0.12500000000000 0.50000000000000 0.32046657352454 114.818 > 1 4 In 0.62500000000000 0.50000000000000 0.32046657352454 114.818 > 1 5 In 0.41428471195709 0.37500000000000 0.25000000000000 114.818 > 5 6 In 0.91428471195709 0.37500000000000 0.25000000000000 114.818 > 5 7 In 0.41428471195709 0.87500000000000 0.25000000000000 114.818 > 5 8 In 0.91428471195709 0.87500000000000 0.25000000000000 114.818 > 5 9 In 0.37500000000000 0.00000000000000 0.67953342647546 114.818 > 9 10 In 0.87500000000000 0.00000000000000 0.67953342647546 114.818 > 9 11 In 0.37500000000000 0.50000000000000 0.67953342647546 114.818 > 9 12 In 0.87500000000000 0.50000000000000 0.67953342647546 114.818 > 9 13 In 0.08571528804291 0.12500000000000 0.75000000000000 114.818 > 13 14 In 0.58571528804291 0.12500000000000 0.75000000000000 114.818 > 13 15 In 0.08571528804291 0.62500000000000 0.75000000000000 114.818 > 13 16 In 0.58571528804291 0.62500000000000 0.75000000000000 114.818 > 13 17 In 0.33571528804291 0.12500000000000 0.25000000000000 114.818 > 17 18 In 0.83571528804291 0.12500000000000 0.25000000000000 114.818 > 17 19 In 0.33571528804291 0.62500000000000 0.25000000000000 114.818 > 17 20 In 0.83571528804291 0.62500000000000 0.25000000000000 114.818 > 17 21 In 0.37500000000000 0.25000000000000 0.82046657352454 114.818 > 21 22 In 0.87500000000000 0.25000000000000 0.82046657352454 114.818 > 21 23 In 0.37500000000000 0.75000000000000 0.82046657352454 114.818 > 21 24 In 0.87500000000000 0.75000000000000 0.82046657352454 114.818 > 21 25 In 0.16428471195709 0.37500000000000 0.75000000000000 114.818 > 25 26 In 0.66428471195709 0.37500000000000 0.75000000000000 114.818 > 25 27 In 0.16428471195709 0.87500000000000 0.75000000000000 114.818 > 25 28 In 0.66428471195709 0.87500000000000 0.75000000000000 114.818 > 25 29 In 0.12500000000000 0.25000000000000 0.17953342647546 114.818 > 29 30 In 0.62500000000000 0.25000000000000 0.17953342647546 114.818 > 29 31 In 0.12500000000000 0.75000000000000 0.17953342647546 114.818 > 29 32 In 0.62500000000000 0.75000000000000 0.17953342647546 114.818 > 29 33 Br 0.49631111547161 0.22629288451405 0.32037942549330 79.904 > 33 34 Br 0.99631111547161 0.22629288451405 0.32037942549330 79.904 > 33 35 Br 0.49631111547161 0.72629288451405 0.32037942549330 79.904 > 33 36 Br 0.99631111547161 0.72629288451405 0.32037942549330 79.904 > 33 37 Br 0.00368888452839 0.27370711548595 0.67962057450670 79.904 > 37 38 Br 0.50368888452839 0.27370711548595 0.67962057450670 79.904 > 37 39 Br 0.00368888452839 0.77370711548595 0.67962057450670 79.904 > 37 40 Br 0.50368888452839 0.77370711548595 0.67962057450670 79.904 > 37 41 Br 0.24631111547161 0.02370711548595 0.82037942549330 79.904 > 41 42 Br 0.74631111547161 0.02370711548595 0.82037942549330 79.904 > 41 43 Br 0.24631111547161 0.52370711548595 0.82037942549330 79.904 > 41 44 Br 0.74631111547161 0.52370711548595 0.82037942549330 79.904 > 41 45 Br 0.17323709877426 0.15306605766268 0.04422525730086 79.904 > 45 46 Br 0.67323709877426 0.15306605766268 0.04422525730086 79.904 > 45 47 Br 0.17323709877426 0.65306605766268 0.04422525730086 79.904 > 45 48 Br 0.67323709877426 0.65306605766268 0.04422525730086 79.904 > 45 49 Br 0.00368888452839 0.47629288451405 0.82037942549330 79.904 > 49 50 Br 0.50368888452839 0.47629288451405 0.82037942549330 79.904 > 49 51 Br 0.00368888452839 0.97629288451405 0.82037942549330 79.904 > 49 52 Br 0.50368888452839 0.97629288451405 0.82037942549330 79.904 > 49 53 Br 0.42323709877426 0.15306605766268 0.95577474269914 79.904 > 53 54 Br 0.92323709877426 0.15306605766268 0.95577474269914 79.904 > 53 55 Br 0.42323709877426 0.65306605766268 0.95577474269914 79.904 > 53 56 Br 0.92323709877426 0.65306605766268 0.95577474269914 79.904 > 53 57 Br 0.17323709877426 0.09693394233732 0.45577474269914 79.904 > 57 58 Br 0.67323709877426 0.09693394233732 0.45577474269914 79.904 > 57 59 Br 0.17323709877426 0.59693394233732 0.45577474269914 79.904 > 57 60 Br 0.67323709877426 0.59693394233732 0.45577474269914 79.904 > 57 61 Br 0.42323709877426 0.09693394233732 0.54422525730086 79.904 > 61 62 Br 0.92323709877426 0.09693394233732 0.54422525730086 79.904 > 61 63 Br 0.42323709877426 0.59693394233732 0.54422525730086 79.904 > 61 64 Br 0.92323709877426 0.59693394233732 0.54422525730086 79.904 > 61 65 Br 0.07676290122574 0.40306605766268 0.45577474269914 79.904 > 65 66 Br 0.57676290122574 0.40306605766268 0.45577474269914 79.904 > 65 67 Br 0.07676290122574 0.90306605766268 0.45577474269914 79.904 > 65 68 Br 0.57676290122574 0.90306605766268 0.45577474269914 79.904 > 65 69 Br 0.24631111547161 0.22629288451405 0.67962057450670 79.904 > 69 70 Br 0.74631111547161 0.22629288451405 0.67962057450670 79.904 > 69 71 Br 0.24631111547161 0.72629288451405 0.67962057450670 79.904 > 69 72 Br 0.74631111547161 0.72629288451405 0.67962057450670 79.904 > 69 73 Br 0.32676290122574 0.40306605766268 0.54422525730086 79.904 > 73 74 Br 0.82676290122574 0.40306605766268 0.54422525730086 79.904 > 73 75 Br 0.32676290122574 0.90306605766268 0.54422525730086 79.904 > 73 76 Br 0.82676290122574 0.90306605766268 0.54422525730086 79.904 > 73 77 Br 0.49631111547161 0.02370711548595 0.17962057450670 79.904 > 77 78 Br 0.99631111547161 0.02370711548595 0.17962057450670 79.904 > 77 79 Br 0.49631111547161 0.52370711548595 0.17962057450670 79.904 > 77 80 Br 0.99631111547161 0.52370711548595 0.17962057450670 79.904 > 77 81 Br 0.07676290122574 0.34693394233732 0.04422525730086 79.904 > 81 82 Br 0.57676290122574 0.34693394233732 0.04422525730086 79.904 > 81 83 Br 0.07676290122574 0.84693394233732 0.04422525730086 79.904 > 81 84 Br 0.57676290122574 0.84693394233732 0.04422525730086 79.904 > 81 85 Br 0.25368888452839 0.27370711548595 0.32037942549330 79.904 > 85 86 Br 0.75368888452839 0.27370711548595 0.32037942549330 79.904 > 85 87 Br 0.25368888452839 0.77370711548595 0.32037942549330 79.904 > 85 88 Br 0.75368888452839 0.77370711548595 0.32037942549330 79.904 > 85 89 Br 0.32676290122574 0.34693394233732 0.95577474269914 79.904 > 89 90 Br 0.82676290122574 0.34693394233732 0.95577474269914 79.904 > 89 91 Br 0.32676290122574 0.84693394233732 0.95577474269914 79.904 > 89 92 Br 0.82676290122574 0.84693394233732 0.95577474269914 79.904 > 89 93 Br 0.25368888452839 0.47629288451405 0.17962057450670 79.904 > 93 94 Br 0.75368888452839 0.47629288451405 0.17962057450670 79.904 > 93 95 Br 0.25368888452839 0.97629288451405 0.17962057450670 79.904 > 93 96 Br 0.75368888452839 0.97629288451405 0.17962057450670 79.904 > 93 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 4.7068177 0.0000000 0.0000000 0.0000000 4.3610801 0.0000000 0.0000000 0.0000000 4.2926020 -------------------------- Born effective charges -------------------------- 1 In 3.2009179 0.3083122 -0.0000000 0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 2 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 0.2056455 -0.0000000 0.3485131 1.9142448 3 In 3.2009179 0.3083122 -0.0000000 0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 4 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 0.2056455 -0.0000000 0.3485131 1.9142448 5 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 -0.2056455 -0.0000000 -0.3485131 1.9142448 6 In 3.2009179 -0.3083122 -0.0000000 -0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 7 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 -0.2056455 -0.0000000 -0.3485131 1.9142448 8 In 3.2009179 -0.3083122 -0.0000000 -0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 9 Br -1.7506695 -0.1377095 0.5111143 -0.2532621 -1.2356452 0.2630911 0.4907325 0.1799490 -0.7513838 10 Br -1.7506695 -0.1377095 0.5111143 -0.2532621 -1.2356452 0.2630911 0.4907325 0.1799490 -0.7513838 11 Br -1.7506695 0.1377095 0.5111143 0.2532621 -1.2356452 -0.2630911 0.4907325 -0.1799490 -0.7513838 12 Br -1.0555407 0.2666999 -0.3476411 0.2653963 -1.1258863 -0.3852760 -0.2679068 -0.3901344 -1.3658623 13 Br -1.7506695 0.1377095 -0.5111143 0.2532621 -1.2356452 0.2630911 -0.4907325 0.1799490 -0.7513838 14 Br -1.0555407 0.2666999 0.3476411 0.2653963 -1.1258863 0.3852760 0.2679068 0.3901344 -1.3658623 15 Br -1.0555407 -0.2666999 0.3476411 -0.2653963 -1.1258863 -0.3852760 0.2679068 -0.3901344 -1.3658623 16 Br -1.0555407 -0.2666999 -0.3476411 -0.2653963 -1.1258863 0.3852760 -0.2679068 0.3901344 -1.3658623 17 Br -1.0555407 -0.2666999 -0.3476411 -0.2653963 -1.1258863 0.3852760 -0.2679068 0.3901344 -1.3658623 18 Br -1.7506695 -0.1377095 -0.5111143 -0.2532621 -1.2356452 -0.2630911 -0.4907325 -0.1799490 -0.7513838 19 Br -1.0555407 -0.2666999 0.3476411 -0.2653963 -1.1258863 -0.3852760 0.2679068 -0.3901344 -1.3658623 20 Br -1.7506695 0.1377095 -0.5111143 0.2532621 -1.2356452 0.2630911 -0.4907325 0.1799490 -0.7513838 21 Br -1.0555407 0.2666999 0.3476411 0.2653963 -1.1258863 0.3852760 0.2679068 0.3901344 -1.3658623 22 Br -1.7506695 -0.1377095 -0.5111143 -0.2532621 -1.2356452 -0.2630911 -0.4907325 -0.1799490 -0.7513838 23 Br -1.0555407 0.2666999 -0.3476411 0.2653963 -1.1258863 -0.3852760 -0.2679068 -0.3901344 -1.3658623 24 Br -1.7506695 0.1377095 0.5111143 0.2532621 -1.2356452 -0.2630911 0.4907325 -0.1799490 -0.7513838 ---------------------------------------------------------------------------- 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.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 5, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 45, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 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.00000001 (zzz) -0.00000001 (zzz) -0.00000001 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (yz) -0.00000000 (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:56:00]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 17:56:00]------------------------- Compiled with OpenMP support (max 48 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [2 2 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pnna (52) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 7.824470680000000 0.000000000000000 0.000000000000000 b 0.000000000000000 10.136594210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.287178179999998 Atomic positions (fractional): 1 In 0.25000000000000 0.00000000000000 0.32046657352454 114.818 2 In 0.82856942391417 0.75000000000000 0.25000000000000 114.818 3 In 0.75000000000000 0.00000000000000 0.67953342647546 114.818 4 In 0.17143057608583 0.25000000000000 0.75000000000000 114.818 5 In 0.67143057608583 0.25000000000000 0.25000000000000 114.818 6 In 0.75000000000000 0.50000000000000 0.82046657352454 114.818 7 In 0.32856942391417 0.75000000000000 0.75000000000000 114.818 8 In 0.25000000000000 0.50000000000000 0.17953342647546 114.818 9 Br 0.99262223094321 0.45258576902811 0.32037942549330 79.904 10 Br 0.00737776905679 0.54741423097189 0.67962057450670 79.904 11 Br 0.49262223094321 0.04741423097189 0.82037942549330 79.904 12 Br 0.34647419754853 0.30613211532536 0.04422525730086 79.904 13 Br 0.00737776905679 0.95258576902811 0.82037942549330 79.904 14 Br 0.84647419754853 0.30613211532536 0.95577474269914 79.904 15 Br 0.34647419754853 0.19386788467464 0.45577474269914 79.904 16 Br 0.84647419754853 0.19386788467464 0.54422525730086 79.904 17 Br 0.15352580245147 0.80613211532536 0.45577474269914 79.904 18 Br 0.49262223094321 0.45258576902811 0.67962057450670 79.904 19 Br 0.65352580245147 0.80613211532536 0.54422525730086 79.904 20 Br 0.99262223094321 0.04741423097189 0.17962057450670 79.904 21 Br 0.15352580245147 0.69386788467464 0.04422525730086 79.904 22 Br 0.50737776905679 0.54741423097189 0.32037942549330 79.904 23 Br 0.65352580245147 0.69386788467464 0.95577474269914 79.904 24 Br 0.50737776905679 0.95258576902811 0.17962057450670 79.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.648941360000000 0.000000000000000 0.000000000000000 b 0.000000000000000 20.273188420000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.287178179999998 Atomic positions (fractional): 1 In 0.12500000000000 0.00000000000000 0.32046657352454 114.818 > 1 2 In 0.62500000000000 0.00000000000000 0.32046657352454 114.818 > 1 3 In 0.12500000000000 0.50000000000000 0.32046657352454 114.818 > 1 4 In 0.62500000000000 0.50000000000000 0.32046657352454 114.818 > 1 5 In 0.41428471195709 0.37500000000000 0.25000000000000 114.818 > 5 6 In 0.91428471195709 0.37500000000000 0.25000000000000 114.818 > 5 7 In 0.41428471195709 0.87500000000000 0.25000000000000 114.818 > 5 8 In 0.91428471195709 0.87500000000000 0.25000000000000 114.818 > 5 9 In 0.37500000000000 0.00000000000000 0.67953342647546 114.818 > 9 10 In 0.87500000000000 0.00000000000000 0.67953342647546 114.818 > 9 11 In 0.37500000000000 0.50000000000000 0.67953342647546 114.818 > 9 12 In 0.87500000000000 0.50000000000000 0.67953342647546 114.818 > 9 13 In 0.08571528804291 0.12500000000000 0.75000000000000 114.818 > 13 14 In 0.58571528804291 0.12500000000000 0.75000000000000 114.818 > 13 15 In 0.08571528804291 0.62500000000000 0.75000000000000 114.818 > 13 16 In 0.58571528804291 0.62500000000000 0.75000000000000 114.818 > 13 17 In 0.33571528804291 0.12500000000000 0.25000000000000 114.818 > 17 18 In 0.83571528804291 0.12500000000000 0.25000000000000 114.818 > 17 19 In 0.33571528804291 0.62500000000000 0.25000000000000 114.818 > 17 20 In 0.83571528804291 0.62500000000000 0.25000000000000 114.818 > 17 21 In 0.37500000000000 0.25000000000000 0.82046657352454 114.818 > 21 22 In 0.87500000000000 0.25000000000000 0.82046657352454 114.818 > 21 23 In 0.37500000000000 0.75000000000000 0.82046657352454 114.818 > 21 24 In 0.87500000000000 0.75000000000000 0.82046657352454 114.818 > 21 25 In 0.16428471195709 0.37500000000000 0.75000000000000 114.818 > 25 26 In 0.66428471195709 0.37500000000000 0.75000000000000 114.818 > 25 27 In 0.16428471195709 0.87500000000000 0.75000000000000 114.818 > 25 28 In 0.66428471195709 0.87500000000000 0.75000000000000 114.818 > 25 29 In 0.12500000000000 0.25000000000000 0.17953342647546 114.818 > 29 30 In 0.62500000000000 0.25000000000000 0.17953342647546 114.818 > 29 31 In 0.12500000000000 0.75000000000000 0.17953342647546 114.818 > 29 32 In 0.62500000000000 0.75000000000000 0.17953342647546 114.818 > 29 33 Br 0.49631111547161 0.22629288451405 0.32037942549330 79.904 > 33 34 Br 0.99631111547161 0.22629288451405 0.32037942549330 79.904 > 33 35 Br 0.49631111547161 0.72629288451405 0.32037942549330 79.904 > 33 36 Br 0.99631111547161 0.72629288451405 0.32037942549330 79.904 > 33 37 Br 0.00368888452839 0.27370711548595 0.67962057450670 79.904 > 37 38 Br 0.50368888452839 0.27370711548595 0.67962057450670 79.904 > 37 39 Br 0.00368888452839 0.77370711548595 0.67962057450670 79.904 > 37 40 Br 0.50368888452839 0.77370711548595 0.67962057450670 79.904 > 37 41 Br 0.24631111547161 0.02370711548595 0.82037942549330 79.904 > 41 42 Br 0.74631111547161 0.02370711548595 0.82037942549330 79.904 > 41 43 Br 0.24631111547161 0.52370711548595 0.82037942549330 79.904 > 41 44 Br 0.74631111547161 0.52370711548595 0.82037942549330 79.904 > 41 45 Br 0.17323709877426 0.15306605766268 0.04422525730086 79.904 > 45 46 Br 0.67323709877426 0.15306605766268 0.04422525730086 79.904 > 45 47 Br 0.17323709877426 0.65306605766268 0.04422525730086 79.904 > 45 48 Br 0.67323709877426 0.65306605766268 0.04422525730086 79.904 > 45 49 Br 0.00368888452839 0.47629288451405 0.82037942549330 79.904 > 49 50 Br 0.50368888452839 0.47629288451405 0.82037942549330 79.904 > 49 51 Br 0.00368888452839 0.97629288451405 0.82037942549330 79.904 > 49 52 Br 0.50368888452839 0.97629288451405 0.82037942549330 79.904 > 49 53 Br 0.42323709877426 0.15306605766268 0.95577474269914 79.904 > 53 54 Br 0.92323709877426 0.15306605766268 0.95577474269914 79.904 > 53 55 Br 0.42323709877426 0.65306605766268 0.95577474269914 79.904 > 53 56 Br 0.92323709877426 0.65306605766268 0.95577474269914 79.904 > 53 57 Br 0.17323709877426 0.09693394233732 0.45577474269914 79.904 > 57 58 Br 0.67323709877426 0.09693394233732 0.45577474269914 79.904 > 57 59 Br 0.17323709877426 0.59693394233732 0.45577474269914 79.904 > 57 60 Br 0.67323709877426 0.59693394233732 0.45577474269914 79.904 > 57 61 Br 0.42323709877426 0.09693394233732 0.54422525730086 79.904 > 61 62 Br 0.92323709877426 0.09693394233732 0.54422525730086 79.904 > 61 63 Br 0.42323709877426 0.59693394233732 0.54422525730086 79.904 > 61 64 Br 0.92323709877426 0.59693394233732 0.54422525730086 79.904 > 61 65 Br 0.07676290122574 0.40306605766268 0.45577474269914 79.904 > 65 66 Br 0.57676290122574 0.40306605766268 0.45577474269914 79.904 > 65 67 Br 0.07676290122574 0.90306605766268 0.45577474269914 79.904 > 65 68 Br 0.57676290122574 0.90306605766268 0.45577474269914 79.904 > 65 69 Br 0.24631111547161 0.22629288451405 0.67962057450670 79.904 > 69 70 Br 0.74631111547161 0.22629288451405 0.67962057450670 79.904 > 69 71 Br 0.24631111547161 0.72629288451405 0.67962057450670 79.904 > 69 72 Br 0.74631111547161 0.72629288451405 0.67962057450670 79.904 > 69 73 Br 0.32676290122574 0.40306605766268 0.54422525730086 79.904 > 73 74 Br 0.82676290122574 0.40306605766268 0.54422525730086 79.904 > 73 75 Br 0.32676290122574 0.90306605766268 0.54422525730086 79.904 > 73 76 Br 0.82676290122574 0.90306605766268 0.54422525730086 79.904 > 73 77 Br 0.49631111547161 0.02370711548595 0.17962057450670 79.904 > 77 78 Br 0.99631111547161 0.02370711548595 0.17962057450670 79.904 > 77 79 Br 0.49631111547161 0.52370711548595 0.17962057450670 79.904 > 77 80 Br 0.99631111547161 0.52370711548595 0.17962057450670 79.904 > 77 81 Br 0.07676290122574 0.34693394233732 0.04422525730086 79.904 > 81 82 Br 0.57676290122574 0.34693394233732 0.04422525730086 79.904 > 81 83 Br 0.07676290122574 0.84693394233732 0.04422525730086 79.904 > 81 84 Br 0.57676290122574 0.84693394233732 0.04422525730086 79.904 > 81 85 Br 0.25368888452839 0.27370711548595 0.32037942549330 79.904 > 85 86 Br 0.75368888452839 0.27370711548595 0.32037942549330 79.904 > 85 87 Br 0.25368888452839 0.77370711548595 0.32037942549330 79.904 > 85 88 Br 0.75368888452839 0.77370711548595 0.32037942549330 79.904 > 85 89 Br 0.32676290122574 0.34693394233732 0.95577474269914 79.904 > 89 90 Br 0.82676290122574 0.34693394233732 0.95577474269914 79.904 > 89 91 Br 0.32676290122574 0.84693394233732 0.95577474269914 79.904 > 89 92 Br 0.82676290122574 0.84693394233732 0.95577474269914 79.904 > 89 93 Br 0.25368888452839 0.47629288451405 0.17962057450670 79.904 > 93 94 Br 0.75368888452839 0.47629288451405 0.17962057450670 79.904 > 93 95 Br 0.25368888452839 0.97629288451405 0.17962057450670 79.904 > 93 96 Br 0.75368888452839 0.97629288451405 0.17962057450670 79.904 > 93 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 4.7068177 0.0000000 0.0000000 0.0000000 4.3610801 0.0000000 0.0000000 0.0000000 4.2926020 -------------------------- Born effective charges -------------------------- 1 In 3.2009179 0.3083122 -0.0000000 0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 2 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 0.2056455 -0.0000000 0.3485131 1.9142448 3 In 3.2009179 0.3083122 -0.0000000 0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 4 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 0.2056455 -0.0000000 0.3485131 1.9142448 5 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 -0.2056455 -0.0000000 -0.3485131 1.9142448 6 In 3.2009179 -0.3083122 -0.0000000 -0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 7 In 2.4115026 -0.0000000 -0.0000000 0.0000000 2.0893035 -0.2056455 -0.0000000 -0.3485131 1.9142448 8 In 3.2009179 -0.3083122 -0.0000000 -0.2417425 2.6337595 0.0000000 -0.0000000 -0.0000000 2.3202475 9 Br -1.7506695 -0.1377095 0.5111143 -0.2532621 -1.2356452 0.2630911 0.4907325 0.1799490 -0.7513838 10 Br -1.7506695 -0.1377095 0.5111143 -0.2532621 -1.2356452 0.2630911 0.4907325 0.1799490 -0.7513838 11 Br -1.7506695 0.1377095 0.5111143 0.2532621 -1.2356452 -0.2630911 0.4907325 -0.1799490 -0.7513838 12 Br -1.0555407 0.2666999 -0.3476411 0.2653963 -1.1258863 -0.3852760 -0.2679068 -0.3901344 -1.3658623 13 Br -1.7506695 0.1377095 -0.5111143 0.2532621 -1.2356452 0.2630911 -0.4907325 0.1799490 -0.7513838 14 Br -1.0555407 0.2666999 0.3476411 0.2653963 -1.1258863 0.3852760 0.2679068 0.3901344 -1.3658623 15 Br -1.0555407 -0.2666999 0.3476411 -0.2653963 -1.1258863 -0.3852760 0.2679068 -0.3901344 -1.3658623 16 Br -1.0555407 -0.2666999 -0.3476411 -0.2653963 -1.1258863 0.3852760 -0.2679068 0.3901344 -1.3658623 17 Br -1.0555407 -0.2666999 -0.3476411 -0.2653963 -1.1258863 0.3852760 -0.2679068 0.3901344 -1.3658623 18 Br -1.7506695 -0.1377095 -0.5111143 -0.2532621 -1.2356452 -0.2630911 -0.4907325 -0.1799490 -0.7513838 19 Br -1.0555407 -0.2666999 0.3476411 -0.2653963 -1.1258863 -0.3852760 0.2679068 -0.3901344 -1.3658623 20 Br -1.7506695 0.1377095 -0.5111143 0.2532621 -1.2356452 0.2630911 -0.4907325 0.1799490 -0.7513838 21 Br -1.0555407 0.2666999 0.3476411 0.2653963 -1.1258863 0.3852760 0.2679068 0.3901344 -1.3658623 22 Br -1.7506695 -0.1377095 -0.5111143 -0.2532621 -1.2356452 -0.2630911 -0.4907325 -0.1799490 -0.7513838 23 Br -1.0555407 0.2666999 -0.3476411 0.2653963 -1.1258863 -0.3852760 -0.2679068 -0.3901344 -1.3658623 24 Br -1.7506695 0.1377095 0.5111143 0.2532621 -1.2356452 -0.2630911 0.4907325 -0.1799490 -0.7513838 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000001 (zzz) -0.00000001 (zzz) -0.00000001 (zzz) Max drift of fc2: -0.00000000 (zx) -0.00000000 (zx) --------------------------- 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 ... [ 6 5 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.44, Number of G-points: 323, Lambda: 0.10 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/36) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 36 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.561 ( 0.000 0.000 0.000) 0.000 0.610 ( 0.000 0.000 0.000) 0.000 0.625 ( 0.000 0.000 0.000) 0.000 0.685 ( 0.000 0.000 0.000) 0.000 0.763 ( 0.000 0.000 0.000) 0.000 0.836 ( 0.000 0.000 0.000) 0.000 0.878 ( 0.000 0.000 0.000) 0.000 0.883 ( 0.000 0.000 0.000) 0.000 0.894 ( 0.000 0.000 0.000) 0.000 0.947 ( 0.000 0.000 0.000) 0.000 1.090 ( 0.000 0.000 0.000) 0.000 1.147 ( 0.000 0.000 0.000) 0.000 1.157 ( 0.000 0.000 0.000) 0.000 1.159 ( 0.000 0.000 0.000) 0.000 1.269 ( 0.000 0.000 0.000) 0.000 1.270 ( 0.000 0.000 0.000) 0.000 1.281 ( 0.000 0.000 0.000) 0.000 1.306 ( 0.000 0.000 0.000) 0.000 1.322 ( 0.000 0.000 0.000) 0.000 1.340 ( 0.000 0.000 0.000) 0.000 1.347 ( 0.000 0.000 0.000) 0.000 1.380 ( 0.000 0.000 0.000) 0.000 1.382 ( 0.000 0.000 0.000) 0.000 1.412 ( 0.000 0.000 0.000) 0.000 1.415 ( 0.000 0.000 0.000) 0.000 1.475 ( 0.000 0.000 0.000) 0.000 1.492 ( 0.000 0.000 0.000) 0.000 1.549 ( 0.000 0.000 0.000) 0.000 1.587 ( 0.000 0.000 0.000) 0.000 1.615 ( 0.000 0.000 0.000) 0.000 1.689 ( 0.000 0.000 0.000) 0.000 1.726 ( 0.000 0.000 0.000) 0.000 1.739 ( 0.000 0.000 0.000) 0.000 1.811 ( 0.000 0.000 0.000) 0.000 1.838 ( 0.000 0.000 0.000) 0.000 1.862 ( 0.000 0.000 0.000) 0.000 1.915 ( 0.000 0.000 0.000) 0.000 1.953 ( 0.000 0.000 0.000) 0.000 1.979 ( 0.000 0.000 0.000) 0.000 2.091 ( 0.000 0.000 0.000) 0.000 2.095 ( 0.000 0.000 0.000) 0.000 2.139 ( 0.000 0.000 0.000) 0.000 2.158 ( 0.000 0.000 0.000) 0.000 2.242 ( 0.000 0.000 0.000) 0.000 2.299 ( 0.000 0.000 0.000) 0.000 2.321 ( 0.000 0.000 0.000) 0.000 2.361 ( 0.000 0.000 0.000) 0.000 2.394 ( 0.000 0.000 0.000) 0.000 2.446 ( 0.000 0.000 0.000) 0.000 2.488 ( 0.000 0.000 0.000) 0.000 2.489 ( 0.000 0.000 0.000) 0.000 2.538 ( 0.000 0.000 0.000) 0.000 2.746 ( 0.000 0.000 0.000) 0.000 5.462 ( 0.000 0.000 0.000) 0.000 5.493 ( 0.000 0.000 0.000) 0.000 5.508 ( 0.000 0.000 0.000) 0.000 5.518 ( 0.000 0.000 0.000) 0.000 6.246 ( 0.000 0.000 0.000) 0.000 6.294 ( 0.000 0.000 0.000) 0.000 6.342 ( 0.000 0.000 0.000) 0.000 6.443 ( 0.000 0.000 0.000) 0.000 6.445 ( 0.000 0.000 0.000) 0.000 6.471 ( 0.000 0.000 0.000) 0.000 6.548 ( 0.000 0.000 0.000) 0.000 6.673 ( 0.000 0.000 0.000) 0.000 6.840 ( 0.000 0.000 0.000) 0.000 6.870 ( 0.000 0.000 0.000) 0.000 6.978 ( 0.000 0.000 0.000) 0.000 6.995 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/36) ======================= q-point: ( 0.17 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 27 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.189 ( 8.390 0.000 0.000) 8.390 0.272 ( 11.980 0.000 0.000) 11.980 0.481 ( 20.006 0.000 0.000) 20.006 0.592 ( 2.806 0.000 0.000) 2.806 0.624 ( 0.893 0.000 0.000) 0.893 0.684 ( 5.096 0.000 0.000) 5.096 0.717 ( -3.555 0.000 0.000) 3.555 0.793 ( 8.305 0.000 0.000) 8.305 0.847 ( 0.585 0.000 0.000) 0.585 0.878 ( 0.215 0.000 0.000) 0.215 0.905 ( -0.193 0.000 0.000) 0.193 0.909 ( 1.249 0.000 0.000) 1.249 0.952 ( -0.056 0.000 0.000) 0.056 1.058 ( -2.292 0.000 0.000) 2.292 1.145 ( 0.514 0.000 0.000) 0.514 1.157 ( -5.829 0.000 0.000) 5.829 1.166 ( -3.225 0.000 0.000) 3.225 1.173 ( -0.230 0.000 0.000) 0.230 1.181 ( -3.653 0.000 0.000) 3.653 1.224 ( -1.548 0.000 0.000) 1.548 1.275 ( -4.691 0.000 0.000) 4.691 1.314 ( 3.463 0.000 0.000) 3.463 1.334 ( -3.722 0.000 0.000) 3.722 1.343 ( -1.421 0.000 0.000) 1.421 1.382 ( 3.224 0.000 0.000) 3.224 1.416 ( -0.705 0.000 0.000) 0.705 1.462 ( -0.110 0.000 0.000) 0.110 1.475 ( 6.367 0.000 0.000) 6.367 1.493 ( 2.082 0.000 0.000) 2.082 1.497 ( 1.182 0.000 0.000) 1.182 1.580 ( 0.906 0.000 0.000) 0.906 1.614 ( 2.485 0.000 0.000) 2.485 1.671 ( 3.050 0.000 0.000) 3.050 1.692 ( 0.519 0.000 0.000) 0.519 1.769 ( 2.318 0.000 0.000) 2.318 1.776 ( -3.368 0.000 0.000) 3.368 1.848 ( -0.845 0.000 0.000) 0.845 1.852 ( 0.861 0.000 0.000) 0.861 1.893 ( -2.450 0.000 0.000) 2.450 1.935 ( -1.203 0.000 0.000) 1.203 1.987 ( 0.396 0.000 0.000) 0.396 2.050 ( -3.119 0.000 0.000) 3.119 2.076 ( -1.312 0.000 0.000) 1.312 2.146 ( -0.092 0.000 0.000) 0.092 2.181 ( 3.292 0.000 0.000) 3.292 2.217 ( -2.262 0.000 0.000) 2.262 2.256 ( -1.580 0.000 0.000) 1.580 2.302 ( 0.170 0.000 0.000) 0.170 2.330 ( -2.641 0.000 0.000) 2.641 2.382 ( -1.535 0.000 0.000) 1.535 2.400 ( 1.212 0.000 0.000) 1.212 2.437 ( -0.799 0.000 0.000) 0.799 2.464 ( -1.768 0.000 0.000) 1.768 2.529 ( -0.921 0.000 0.000) 0.921 2.725 ( -1.664 0.000 0.000) 1.664 2.811 ( -1.168 0.000 0.000) 1.168 5.459 ( -0.216 0.000 0.000) 0.216 5.488 ( -0.427 0.000 0.000) 0.427 5.506 ( -0.874 0.000 0.000) 0.874 5.506 ( -0.340 0.000 0.000) 0.340 6.257 ( 0.946 0.000 0.000) 0.946 6.308 ( 1.240 0.000 0.000) 1.240 6.336 ( -0.462 0.000 0.000) 0.462 6.457 ( -1.282 0.000 0.000) 1.282 6.458 ( 1.194 0.000 0.000) 1.194 6.543 ( -0.247 0.000 0.000) 0.247 6.653 ( -1.848 0.000 0.000) 1.848 6.707 ( -2.230 0.000 0.000) 2.230 6.848 ( 0.672 0.000 0.000) 0.672 6.873 ( 0.230 0.000 0.000) 0.230 7.012 ( 2.998 0.000 0.000) 2.998 7.243 ( -0.826 0.000 0.000) 0.826 ======================= Grid point 2 (3/36) ======================= q-point: ( 0.33 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.348 ( 6.430 0.000 0.000) 6.430 0.481 ( 5.745 0.000 0.000) 5.745 0.648 ( 1.364 0.000 0.000) 1.364 0.651 ( -0.881 0.000 0.000) 0.881 0.671 ( 4.091 0.000 0.000) 4.091 0.708 ( 1.324 0.000 0.000) 1.324 0.812 ( 6.131 0.000 0.000) 6.131 0.853 ( -0.020 0.000 0.000) 0.020 0.859 ( -2.127 0.000 0.000) 2.127 0.927 ( -2.112 0.000 0.000) 2.112 0.934 ( 0.293 0.000 0.000) 0.293 0.944 ( 3.144 0.000 0.000) 3.144 0.970 ( 2.889 0.000 0.000) 2.889 1.031 ( -4.185 0.000 0.000) 4.185 1.050 ( -6.372 0.000 0.000) 6.372 1.091 ( -2.892 0.000 0.000) 2.892 1.122 ( 11.475 0.000 0.000) 11.475 1.174 ( -1.031 0.000 0.000) 1.031 1.187 ( 3.967 0.000 0.000) 3.967 1.226 ( -1.444 0.000 0.000) 1.444 1.241 ( 2.458 0.000 0.000) 2.458 1.280 ( 0.103 0.000 0.000) 0.103 1.318 ( 0.274 0.000 0.000) 0.274 1.363 ( -3.002 0.000 0.000) 3.002 1.406 ( 4.917 0.000 0.000) 4.917 1.457 ( 0.344 0.000 0.000) 0.344 1.483 ( 5.757 0.000 0.000) 5.757 1.513 ( -0.138 0.000 0.000) 0.138 1.519 ( 0.311 0.000 0.000) 0.311 1.586 ( 2.958 0.000 0.000) 2.958 1.601 ( 1.023 0.000 0.000) 1.023 1.661 ( -7.534 0.000 0.000) 7.534 1.672 ( -2.795 0.000 0.000) 2.795 1.683 ( 3.448 0.000 0.000) 3.448 1.717 ( 1.722 0.000 0.000) 1.722 1.752 ( -8.734 0.000 0.000) 8.734 1.826 ( 2.459 0.000 0.000) 2.459 1.833 ( -2.395 0.000 0.000) 2.395 1.869 ( 0.338 0.000 0.000) 0.338 1.949 ( 2.497 0.000 0.000) 2.497 1.954 ( -4.582 0.000 0.000) 4.582 1.977 ( -0.937 0.000 0.000) 0.937 2.047 ( -1.338 0.000 0.000) 1.338 2.146 ( -4.506 0.000 0.000) 4.506 2.194 ( 4.464 0.000 0.000) 4.464 2.237 ( 2.051 0.000 0.000) 2.051 2.274 ( -1.928 0.000 0.000) 1.928 2.298 ( 4.064 0.000 0.000) 4.064 2.300 ( -0.489 0.000 0.000) 0.489 2.320 ( -3.538 0.000 0.000) 3.538 2.390 ( -1.897 0.000 0.000) 1.897 2.421 ( -0.491 0.000 0.000) 0.491 2.437 ( -0.893 0.000 0.000) 0.893 2.490 ( -2.946 0.000 0.000) 2.946 2.700 ( -0.040 0.000 0.000) 0.040 2.775 ( -2.166 0.000 0.000) 2.166 5.459 ( 0.375 0.000 0.000) 0.375 5.478 ( -0.477 0.000 0.000) 0.477 5.484 ( -0.925 0.000 0.000) 0.925 5.496 ( -0.802 0.000 0.000) 0.802 6.281 ( 1.160 0.000 0.000) 1.160 6.323 ( -0.781 0.000 0.000) 0.781 6.341 ( 1.734 0.000 0.000) 1.734 6.421 ( -1.907 0.000 0.000) 1.907 6.494 ( 2.112 0.000 0.000) 2.112 6.547 ( 0.856 0.000 0.000) 0.856 6.602 ( -2.673 0.000 0.000) 2.673 6.647 ( -3.181 0.000 0.000) 3.181 6.864 ( 0.736 0.000 0.000) 0.736 6.877 ( 0.113 0.000 0.000) 0.113 7.089 ( 3.855 0.000 0.000) 3.855 7.215 ( -1.790 0.000 0.000) 1.790 ======================= Grid point 3 (4/36) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 18 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.464 ( -4.399 0.000 0.000) 4.399 0.464 ( 4.399 0.000 0.000) 4.399 0.683 ( -1.967 0.000 0.000) 1.967 0.683 ( 1.967 0.000 0.000) 1.967 0.711 ( -2.165 0.000 0.000) 2.165 0.711 ( 2.165 0.000 0.000) 2.165 0.849 ( -0.432 0.000 0.000) 0.432 0.849 ( 0.432 0.000 0.000) 0.432 0.899 ( -0.466 0.000 0.000) 0.466 0.899 ( 0.466 0.000 0.000) 0.466 0.926 ( -4.658 0.000 0.000) 4.658 0.926 ( 4.658 0.000 0.000) 4.658 0.966 ( -0.902 0.000 0.000) 0.902 0.966 ( 0.902 0.000 0.000) 0.902 1.033 ( -2.839 0.000 0.000) 2.839 1.033 ( 2.839 0.000 0.000) 2.839 1.178 ( -1.943 0.000 0.000) 1.943 1.178 ( 1.943 0.000 0.000) 1.943 1.232 ( -2.557 0.000 0.000) 2.557 1.232 ( 2.557 0.000 0.000) 2.557 1.311 ( -1.963 0.000 0.000) 1.963 1.311 ( 1.963 0.000 0.000) 1.963 1.452 ( -11.540 0.000 0.000) 11.540 1.452 ( 11.540 0.000 0.000) 11.540 1.464 ( -10.122 0.000 0.000) 10.122 1.464 ( 10.122 0.000 0.000) 10.122 1.481 ( -1.758 0.000 0.000) 1.758 1.481 ( 1.758 0.000 0.000) 1.758 1.516 ( -6.651 0.000 0.000) 6.651 1.516 ( 6.651 0.000 0.000) 6.651 1.591 ( -2.637 0.000 0.000) 2.637 1.591 ( 2.637 0.000 0.000) 2.637 1.631 ( -2.038 0.000 0.000) 2.038 1.631 ( 2.038 0.000 0.000) 2.038 1.737 ( -1.181 0.000 0.000) 1.181 1.737 ( 1.181 0.000 0.000) 1.181 1.815 ( -2.000 0.000 0.000) 2.000 1.815 ( 2.000 0.000 0.000) 2.000 1.886 ( -1.740 0.000 0.000) 1.740 1.886 ( 1.740 0.000 0.000) 1.740 2.005 ( -1.914 0.000 0.000) 1.914 2.005 ( 1.914 0.000 0.000) 1.914 2.032 ( -5.001 0.000 0.000) 5.001 2.032 ( 5.001 0.000 0.000) 5.001 2.266 ( -1.218 0.000 0.000) 1.218 2.266 ( 1.218 0.000 0.000) 1.218 2.278 ( -1.588 0.000 0.000) 1.588 2.278 ( 1.588 0.000 0.000) 1.588 2.298 ( -4.350 0.000 0.000) 4.350 2.298 ( 4.350 0.000 0.000) 4.350 2.410 ( -4.093 0.000 0.000) 4.093 2.410 ( 4.093 0.000 0.000) 4.093 2.421 ( -0.505 0.000 0.000) 0.505 2.421 ( 0.505 0.000 0.000) 0.505 2.726 ( -2.127 0.000 0.000) 2.127 2.726 ( 2.127 0.000 0.000) 2.127 5.473 ( -0.111 0.000 0.000) 0.111 5.473 ( 0.111 0.000 0.000) 0.111 5.475 ( -0.990 0.000 0.000) 0.990 5.475 ( 0.990 0.000 0.000) 0.990 6.304 ( -0.985 0.000 0.000) 0.985 6.304 ( 0.985 0.000 0.000) 0.985 6.380 ( -1.894 0.000 0.000) 1.894 6.380 ( 1.894 0.000 0.000) 1.894 6.545 ( -2.621 0.000 0.000) 2.621 6.545 ( 2.621 0.000 0.000) 2.621 6.584 ( -2.493 0.000 0.000) 2.493 6.584 ( 2.493 0.000 0.000) 2.493 6.876 ( -0.322 0.000 0.000) 0.322 6.876 ( 0.322 0.000 0.000) 0.322 7.164 ( -3.008 0.000 0.000) 3.008 7.164 ( 3.008 0.000 0.000) 3.008 ======================= Grid point 7 (5/36) ======================= q-point: ( 0.00 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.177 ( 0.000 8.701 0.000) 8.701 0.252 ( 0.000 12.278 0.000) 12.278 0.378 ( 0.000 18.653 0.000) 18.653 0.555 ( 0.000 -2.175 0.000) 2.175 0.629 ( 0.000 1.887 0.000) 1.887 0.630 ( 0.000 1.922 0.000) 1.922 0.719 ( 0.000 3.357 0.000) 3.357 0.761 ( 0.000 -5.107 0.000) 5.107 0.807 ( 0.000 -1.810 0.000) 1.810 0.874 ( 0.000 -0.909 0.000) 0.909 0.907 ( 0.000 2.004 0.000) 2.004 0.925 ( 0.000 2.424 0.000) 2.424 0.944 ( 0.000 -0.617 0.000) 0.617 1.019 ( 0.000 -4.984 0.000) 4.984 1.132 ( 0.000 -4.060 0.000) 4.060 1.142 ( 0.000 -1.488 0.000) 1.488 1.177 ( 0.000 -7.085 0.000) 7.085 1.188 ( 0.000 -7.177 0.000) 7.177 1.191 ( 0.000 -2.278 0.000) 2.278 1.222 ( 0.000 2.888 0.000) 2.888 1.282 ( 0.000 -2.813 0.000) 2.813 1.324 ( 0.000 2.837 0.000) 2.837 1.335 ( 0.000 0.797 0.000) 0.797 1.378 ( 0.000 -0.407 0.000) 0.407 1.380 ( 0.000 -0.652 0.000) 0.652 1.430 ( 0.000 3.953 0.000) 3.953 1.448 ( 0.000 1.817 0.000) 1.817 1.483 ( 0.000 -5.091 0.000) 5.091 1.489 ( 0.000 2.707 0.000) 2.707 1.514 ( 0.000 4.915 0.000) 4.915 1.592 ( 0.000 0.626 0.000) 0.626 1.615 ( 0.000 -0.116 0.000) 0.116 1.698 ( 0.000 -1.301 0.000) 1.301 1.741 ( 0.000 2.311 0.000) 2.311 1.751 ( 0.000 1.100 0.000) 1.100 1.771 ( 0.000 -1.069 0.000) 1.069 1.824 ( 0.000 -1.083 0.000) 1.083 1.832 ( 0.000 -2.682 0.000) 2.682 1.839 ( 0.000 1.775 0.000) 1.775 1.844 ( 0.000 -2.865 0.000) 2.865 2.072 ( 0.000 3.396 0.000) 3.396 2.073 ( 0.000 -1.727 0.000) 1.727 2.094 ( 0.000 -0.079 0.000) 0.079 2.115 ( 0.000 -2.302 0.000) 2.302 2.169 ( 0.000 2.675 0.000) 2.675 2.191 ( 0.000 0.611 0.000) 0.611 2.237 ( 0.000 -0.406 0.000) 0.406 2.286 ( 0.000 -0.693 0.000) 0.693 2.335 ( 0.000 1.256 0.000) 1.256 2.391 ( 0.000 2.597 0.000) 2.597 2.443 ( 0.000 -2.105 0.000) 2.105 2.460 ( 0.000 0.830 0.000) 0.830 2.537 ( 0.000 0.128 0.000) 0.128 2.551 ( 0.000 3.429 0.000) 3.429 2.696 ( 0.000 -0.963 0.000) 0.963 2.719 ( 0.000 -2.691 0.000) 2.691 5.467 ( 0.000 0.421 0.000) 0.421 5.495 ( 0.000 0.153 0.000) 0.153 5.507 ( 0.000 -0.123 0.000) 0.123 5.512 ( 0.000 -0.609 0.000) 0.609 6.297 ( 0.000 0.272 0.000) 0.272 6.336 ( 0.000 -0.580 0.000) 0.580 6.441 ( 0.000 -0.123 0.000) 0.123 6.447 ( 0.000 -0.330 0.000) 0.330 6.455 ( 0.000 0.966 0.000) 0.966 6.489 ( 0.000 1.675 0.000) 1.675 6.537 ( 0.000 -1.038 0.000) 1.038 6.648 ( 0.000 -2.339 0.000) 2.339 6.855 ( 0.000 1.423 0.000) 1.423 6.966 ( 0.000 -1.205 0.000) 1.205 7.006 ( 0.000 0.942 0.000) 0.942 7.053 ( 0.000 -0.379 0.000) 0.379 ======================= Grid point 8 (6/36) ======================= q-point: ( 0.17 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 45 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.266 ( 6.538 6.377 0.000) 9.133 0.307 ( 8.661 5.446 0.000) 10.231 0.580 ( 2.273 -1.615 0.000) 2.789 0.587 ( 9.492 6.967 0.000) 11.774 0.616 ( 0.258 -0.185 0.000) 0.317 0.678 ( 3.449 -0.412 0.000) 3.473 0.721 ( -1.199 -1.845 0.000) 2.200 0.820 ( 6.441 2.225 0.000) 6.815 0.826 ( 1.740 1.746 0.000) 2.465 0.862 ( -2.077 -3.375 0.000) 3.963 0.863 ( -0.580 -1.243 0.000) 1.371 0.940 ( 0.989 1.750 0.000) 2.010 0.968 ( 0.484 -0.717 0.000) 0.865 1.019 ( -0.228 -2.332 0.000) 2.343 1.086 ( -0.717 -3.300 0.000) 3.377 1.118 ( 0.441 -1.915 0.000) 1.965 1.151 ( 1.194 -1.926 0.000) 2.266 1.167 ( -2.339 -1.097 0.000) 2.584 1.216 ( -4.698 3.009 0.000) 5.579 1.236 ( 1.489 2.645 0.000) 3.035 1.304 ( -0.129 -2.770 0.000) 2.773 1.310 ( -2.971 -1.078 0.000) 3.160 1.336 ( 2.694 0.826 0.000) 2.818 1.347 ( -1.446 -0.015 0.000) 1.446 1.392 ( 2.958 0.217 0.000) 2.965 1.417 ( 0.763 1.308 0.000) 1.514 1.455 ( 0.468 2.807 0.000) 2.845 1.460 ( -1.455 0.017 0.000) 1.455 1.523 ( 2.610 -1.919 0.000) 3.239 1.566 ( -3.422 3.710 0.000) 5.047 1.586 ( 2.294 3.858 0.000) 4.488 1.620 ( 2.621 0.535 0.000) 2.676 1.659 ( 3.099 -0.868 0.000) 3.219 1.718 ( -1.845 1.570 0.000) 2.423 1.743 ( -2.432 -2.488 0.000) 3.479 1.772 ( 1.377 0.154 0.000) 1.386 1.814 ( -0.600 -3.502 0.000) 3.553 1.827 ( -0.356 -1.933 0.000) 1.966 1.851 ( 0.468 -1.777 0.000) 1.837 1.936 ( 1.624 0.361 0.000) 1.663 2.022 ( -1.557 -3.298 0.000) 3.647 2.040 ( -2.592 2.167 0.000) 3.379 2.079 ( -1.483 -0.130 0.000) 1.488 2.138 ( -0.059 -3.503 0.000) 3.503 2.167 ( 0.177 2.803 0.000) 2.808 2.213 ( 0.414 0.486 0.000) 0.638 2.259 ( 1.381 -1.085 0.000) 1.756 2.283 ( -0.227 -0.670 0.000) 0.707 2.311 ( -0.497 0.777 0.000) 0.922 2.370 ( -2.037 2.837 0.000) 3.492 2.422 ( -2.395 -2.212 0.000) 3.260 2.439 ( -1.404 0.071 0.000) 1.406 2.532 ( -1.286 4.893 0.000) 5.059 2.534 ( -0.537 0.542 0.000) 0.763 2.697 ( -1.744 -2.726 0.000) 3.236 2.747 ( 1.047 -4.147 0.000) 4.278 5.462 ( -0.272 0.364 0.000) 0.454 5.490 ( -0.401 0.197 0.000) 0.446 5.501 ( -0.873 -0.502 0.000) 1.008 5.505 ( -0.254 -0.144 0.000) 0.292 6.310 ( 1.118 0.232 0.000) 1.142 6.321 ( -1.403 4.654 0.000) 4.861 6.338 ( 0.233 -0.002 0.000) 0.233 6.462 ( 0.652 0.433 0.000) 0.783 6.472 ( -1.324 1.531 0.000) 2.024 6.531 ( -0.365 -1.127 0.000) 1.185 6.604 ( 1.674 -6.019 0.000) 6.247 6.632 ( -1.088 -2.565 0.000) 2.787 6.866 ( 0.989 1.753 0.000) 2.013 6.941 ( -1.680 4.126 0.000) 4.455 6.998 ( 2.791 -1.429 0.000) 3.136 7.155 ( 2.887 -5.411 0.000) 6.133 ======================= Grid point 9 (7/36) ======================= q-point: ( 0.33 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-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.402 ( 5.616 5.039 0.000) 7.545 0.510 ( 7.554 2.490 0.000) 7.954 0.619 ( -1.597 -0.741 0.000) 1.761 0.630 ( 1.574 -1.286 0.000) 2.032 0.651 ( 3.882 -1.792 0.000) 4.275 0.725 ( 0.239 -0.282 0.000) 0.370 0.771 ( 4.052 -3.796 0.000) 5.552 0.790 ( -2.035 -4.860 0.000) 5.268 0.832 ( -0.938 -1.645 0.000) 1.893 0.872 ( -1.575 -4.448 0.000) 4.718 0.882 ( -2.016 -2.224 0.000) 3.002 0.956 ( -1.250 -1.744 0.000) 2.145 0.999 ( 4.314 4.460 0.000) 6.205 1.062 ( -1.999 0.482 0.000) 2.056 1.094 ( -1.717 1.180 0.000) 2.083 1.132 ( -1.808 5.522 0.000) 5.811 1.168 ( -0.374 -0.570 0.000) 0.682 1.174 ( 6.020 3.423 0.000) 6.926 1.207 ( 5.319 1.925 0.000) 5.657 1.242 ( 3.032 -0.291 0.000) 3.046 1.266 ( 2.130 4.447 0.000) 4.931 1.275 ( 0.374 0.125 0.000) 0.395 1.343 ( 0.196 1.762 0.000) 1.773 1.378 ( -4.163 1.331 0.000) 4.371 1.423 ( 4.187 1.145 0.000) 4.341 1.443 ( -0.046 -1.291 0.000) 1.292 1.476 ( 2.562 -2.471 0.000) 3.559 1.477 ( 3.787 -0.373 0.000) 3.805 1.516 ( -1.915 0.274 0.000) 1.935 1.600 ( 3.096 -0.270 0.000) 3.108 1.611 ( -0.353 1.130 0.000) 1.184 1.654 ( -5.216 -0.550 0.000) 5.245 1.662 ( -0.539 -0.579 0.000) 0.791 1.706 ( 2.300 1.548 0.000) 2.772 1.717 ( 1.255 -0.003 0.000) 1.255 1.755 ( -4.672 -0.117 0.000) 4.673 1.799 ( -0.595 -2.354 0.000) 2.428 1.815 ( 0.422 -0.184 0.000) 0.460 1.867 ( 0.486 -0.350 0.000) 0.599 1.950 ( -3.601 -0.481 0.000) 3.633 1.965 ( 0.866 0.674 0.000) 1.097 1.992 ( -0.954 1.317 0.000) 1.626 2.045 ( -1.323 -0.344 0.000) 1.367 2.121 ( -2.050 -2.266 0.000) 3.056 2.192 ( 2.365 0.361 0.000) 2.393 2.231 ( 0.966 -0.697 0.000) 1.192 2.281 ( 0.233 0.278 0.000) 0.362 2.281 ( 1.233 -0.498 0.000) 1.329 2.315 ( -2.557 1.483 0.000) 2.956 2.322 ( 0.663 -0.035 0.000) 0.664 2.358 ( -2.837 -2.072 0.000) 3.513 2.410 ( -1.435 -0.760 0.000) 1.624 2.492 ( -2.115 3.896 0.000) 4.433 2.505 ( -2.081 1.324 0.000) 2.466 2.669 ( -0.258 -2.944 0.000) 2.955 2.732 ( -1.797 -3.764 0.000) 4.171 5.461 ( 0.274 0.191 0.000) 0.334 5.480 ( -0.518 0.201 0.000) 0.556 5.482 ( -0.781 -0.243 0.000) 0.818 5.494 ( -0.770 -0.142 0.000) 0.783 6.318 ( 0.398 3.383 0.000) 3.406 6.328 ( 0.215 0.631 0.000) 0.666 6.355 ( 1.542 1.025 0.000) 1.852 6.436 ( -1.973 1.421 0.000) 2.431 6.484 ( 1.443 -0.839 0.000) 1.670 6.531 ( 0.635 -1.424 0.000) 1.559 6.575 ( -2.460 -2.769 0.000) 3.704 6.606 ( -1.953 -3.194 0.000) 3.744 6.891 ( 1.154 2.485 0.000) 2.740 6.922 ( -0.465 3.749 0.000) 3.778 7.069 ( 3.527 -1.979 0.000) 4.044 7.172 ( -0.742 -3.826 0.000) 3.898 ======================= Grid point 10 (8/36) ======================= q-point: (-0.50 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.501 ( -1.626 3.636 0.000) 3.983 0.501 ( 1.626 3.636 0.000) 3.983 0.675 ( -0.282 -0.602 0.000) 0.665 0.675 ( 0.282 -0.602 0.000) 0.665 0.703 ( -0.269 -0.738 0.000) 0.785 0.703 ( 0.269 -0.738 0.000) 0.785 0.788 ( -0.013 -4.073 0.000) 4.073 0.788 ( 0.013 -4.073 0.000) 4.073 0.823 ( -1.168 -4.610 0.000) 4.755 0.823 ( 1.168 -4.610 0.000) 4.755 0.906 ( -3.075 -2.791 0.000) 4.152 0.906 ( 3.075 -2.791 0.000) 4.152 1.029 ( -0.923 3.571 0.000) 3.688 1.029 ( 0.923 3.571 0.000) 3.688 1.082 ( -0.007 3.464 0.000) 3.464 1.082 ( 0.007 3.464 0.000) 3.464 1.187 ( -1.742 1.190 0.000) 2.110 1.187 ( 1.742 1.190 0.000) 2.110 1.258 ( -1.576 2.271 0.000) 2.764 1.258 ( 1.576 2.271 0.000) 2.764 1.302 ( -2.252 -0.813 0.000) 2.394 1.302 ( 2.252 -0.813 0.000) 2.394 1.455 ( -0.105 -1.130 0.000) 1.135 1.455 ( 0.105 -1.130 0.000) 1.135 1.466 ( -0.140 -0.878 0.000) 0.889 1.466 ( 0.140 -0.878 0.000) 0.889 1.478 ( -3.249 1.005 0.000) 3.401 1.478 ( 3.249 1.005 0.000) 3.401 1.538 ( -5.073 2.186 0.000) 5.524 1.538 ( 5.073 2.186 0.000) 5.524 1.577 ( -1.127 0.838 0.000) 1.404 1.577 ( 1.127 0.838 0.000) 1.404 1.657 ( -1.712 1.196 0.000) 2.088 1.657 ( 1.712 1.196 0.000) 2.088 1.761 ( -0.184 2.047 0.000) 2.055 1.761 ( 0.184 2.047 0.000) 2.055 1.799 ( -1.029 -1.541 0.000) 1.853 1.799 ( 1.029 -1.541 0.000) 1.853 1.871 ( -0.787 -1.222 0.000) 1.454 1.871 ( 0.787 -1.222 0.000) 1.454 2.011 ( -1.717 0.570 0.000) 1.809 2.011 ( 1.717 0.570 0.000) 1.809 2.041 ( -3.945 0.629 0.000) 3.995 2.041 ( 3.945 0.629 0.000) 3.995 2.249 ( -1.133 -1.158 0.000) 1.620 2.249 ( 1.133 -1.158 0.000) 1.620 2.271 ( -0.135 -0.887 0.000) 0.897 2.271 ( 0.135 -0.887 0.000) 0.897 2.307 ( -1.798 0.574 0.000) 1.888 2.307 ( 1.798 0.574 0.000) 1.888 2.369 ( -2.513 -1.978 0.000) 3.198 2.369 ( 2.513 -1.978 0.000) 3.198 2.470 ( -0.040 2.848 0.000) 2.848 2.470 ( 0.040 2.848 0.000) 2.848 2.689 ( -1.641 -3.423 0.000) 3.796 2.689 ( 1.641 -3.423 0.000) 3.796 5.473 ( -0.013 0.024 0.000) 0.027 5.473 ( 0.013 0.024 0.000) 0.027 5.475 ( -0.563 -0.005 0.000) 0.563 5.475 ( 0.563 -0.005 0.000) 0.563 6.326 ( -0.027 2.082 0.000) 2.083 6.326 ( 0.027 2.082 0.000) 2.083 6.393 ( -1.060 1.256 0.000) 1.644 6.393 ( 1.060 1.256 0.000) 1.644 6.523 ( -1.484 -2.043 0.000) 2.525 6.523 ( 1.484 -2.043 0.000) 2.525 6.560 ( -0.480 -2.029 0.000) 2.085 6.560 ( 0.480 -2.029 0.000) 2.085 6.911 ( -0.137 3.199 0.000) 3.202 6.911 ( 0.137 3.199 0.000) 3.202 7.136 ( -1.413 -2.722 0.000) 3.067 7.136 ( 1.413 -2.722 0.000) 3.067 ======================= Grid point 14 (9/36) ======================= q-point: ( 0.00 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.339 ( 0.000 7.573 0.000) 7.573 0.471 ( 0.000 -5.749 0.000) 5.749 0.473 ( 0.000 9.810 0.000) 9.810 0.638 ( 0.000 -6.950 0.000) 6.950 0.678 ( 0.000 2.631 0.000) 2.631 0.688 ( 0.000 2.413 0.000) 2.413 0.718 ( 0.000 14.356 0.000) 14.356 0.803 ( 0.000 4.369 0.000) 4.369 0.805 ( 0.000 2.066 0.000) 2.066 0.848 ( 0.000 -1.182 0.000) 1.182 0.906 ( 0.000 -5.676 0.000) 5.676 0.910 ( 0.000 -3.473 0.000) 3.473 0.953 ( 0.000 -0.783 0.000) 0.783 0.957 ( 0.000 3.332 0.000) 3.332 0.989 ( 0.000 -6.597 0.000) 6.597 1.075 ( 0.000 -6.692 0.000) 6.692 1.080 ( 0.000 -0.996 0.000) 0.996 1.114 ( 0.000 -2.751 0.000) 2.751 1.152 ( 0.000 -0.506 0.000) 0.506 1.180 ( 0.000 -4.094 0.000) 4.094 1.300 ( 0.000 4.154 0.000) 4.154 1.364 ( 0.000 0.305 0.000) 0.305 1.365 ( 0.000 -0.409 0.000) 0.409 1.368 ( 0.000 -0.512 0.000) 0.512 1.370 ( 0.000 1.382 0.000) 1.382 1.390 ( 0.000 -3.522 0.000) 3.522 1.475 ( 0.000 1.469 0.000) 1.469 1.502 ( 0.000 -0.872 0.000) 0.872 1.527 ( 0.000 5.719 0.000) 5.719 1.588 ( 0.000 1.324 0.000) 1.324 1.607 ( 0.000 -0.724 0.000) 0.724 1.637 ( 0.000 -4.798 0.000) 4.798 1.647 ( 0.000 5.291 0.000) 5.291 1.743 ( 0.000 -3.681 0.000) 3.681 1.756 ( 0.000 -0.337 0.000) 0.337 1.759 ( 0.000 -0.499 0.000) 0.499 1.796 ( 0.000 0.569 0.000) 0.569 1.813 ( 0.000 -0.882 0.000) 0.882 1.837 ( 0.000 2.286 0.000) 2.286 1.863 ( 0.000 0.162 0.000) 0.162 2.046 ( 0.000 -0.563 0.000) 0.563 2.064 ( 0.000 -2.276 0.000) 2.276 2.084 ( 0.000 -0.368 0.000) 0.368 2.085 ( 0.000 -0.503 0.000) 0.503 2.209 ( 0.000 1.022 0.000) 1.022 2.226 ( 0.000 -0.734 0.000) 0.734 2.251 ( 0.000 3.429 0.000) 3.429 2.281 ( 0.000 -0.184 0.000) 0.184 2.365 ( 0.000 1.739 0.000) 1.739 2.402 ( 0.000 -1.992 0.000) 1.992 2.441 ( 0.000 1.910 0.000) 1.910 2.464 ( 0.000 -0.509 0.000) 0.509 2.564 ( 0.000 3.074 0.000) 3.074 2.616 ( 0.000 3.054 0.000) 3.054 2.645 ( 0.000 -4.465 0.000) 4.465 2.667 ( 0.000 -2.035 0.000) 2.035 5.479 ( 0.000 0.792 0.000) 0.792 5.496 ( 0.000 -0.905 0.000) 0.905 5.499 ( 0.000 0.247 0.000) 0.247 5.504 ( 0.000 -0.228 0.000) 0.228 6.306 ( 0.000 0.644 0.000) 0.644 6.321 ( 0.000 -0.833 0.000) 0.833 6.439 ( 0.000 -0.025 0.000) 0.025 6.441 ( 0.000 -0.209 0.000) 0.209 6.480 ( 0.000 1.489 0.000) 1.489 6.511 ( 0.000 -1.516 0.000) 1.516 6.533 ( 0.000 2.717 0.000) 2.717 6.592 ( 0.000 -3.106 0.000) 3.106 6.891 ( 0.000 2.119 0.000) 2.119 6.933 ( 0.000 -2.003 0.000) 2.003 7.026 ( 0.000 0.985 0.000) 0.985 7.043 ( 0.000 -0.703 0.000) 0.703 ======================= Grid point 15 (10/36) ======================= q-point: ( 0.17 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 45 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.400 ( 5.221 6.646 0.000) 8.452 0.466 ( 1.591 9.096 0.000) 9.234 0.515 ( 3.824 -4.803 0.000) 6.140 0.621 ( -0.236 -5.684 0.000) 5.689 0.637 ( -2.005 1.953 0.000) 2.799 0.669 ( -1.083 -0.812 0.000) 1.353 0.696 ( -1.212 3.016 0.000) 3.251 0.782 ( -4.006 -4.795 0.000) 6.249 0.831 ( -1.358 -1.089 0.000) 1.740 0.837 ( -0.066 -1.201 0.000) 1.203 0.910 ( 0.771 2.982 0.000) 3.080 0.934 ( -0.055 -1.207 0.000) 1.208 0.978 ( 4.574 3.703 0.000) 5.886 1.021 ( 5.541 3.874 0.000) 6.761 1.063 ( 4.811 -3.094 0.000) 5.720 1.090 ( 2.251 -1.443 0.000) 2.674 1.125 ( 4.389 2.088 0.000) 4.860 1.156 ( 0.172 0.442 0.000) 0.474 1.202 ( 6.452 -4.824 0.000) 8.056 1.211 ( 1.881 -3.669 0.000) 4.123 1.318 ( 0.761 4.196 0.000) 4.264 1.320 ( -3.694 0.812 0.000) 3.782 1.339 ( -2.423 -0.842 0.000) 2.565 1.349 ( 1.092 1.382 0.000) 1.761 1.390 ( 1.790 -0.618 0.000) 1.894 1.417 ( 2.863 -4.209 0.000) 5.090 1.447 ( -1.391 1.163 0.000) 1.814 1.463 ( -2.441 -0.217 0.000) 2.451 1.567 ( 2.073 4.929 0.000) 5.348 1.625 ( 1.381 0.343 0.000) 1.423 1.628 ( 1.316 -0.082 0.000) 1.318 1.647 ( -0.033 0.308 0.000) 0.310 1.676 ( 2.774 3.226 0.000) 4.254 1.679 ( -7.041 -0.203 0.000) 7.044 1.737 ( -4.833 -2.853 0.000) 5.612 1.746 ( 0.249 -3.112 0.000) 3.122 1.782 ( 0.324 3.117 0.000) 3.134 1.812 ( 1.361 0.720 0.000) 1.539 1.829 ( 1.335 -0.803 0.000) 1.558 1.876 ( 2.068 -5.646 0.000) 6.013 2.009 ( -1.543 2.256 0.000) 2.733 2.057 ( -2.148 0.300 0.000) 2.169 2.067 ( -1.724 -0.747 0.000) 1.879 2.068 ( -0.421 -3.413 0.000) 3.439 2.234 ( 0.889 1.269 0.000) 1.549 2.236 ( -1.108 3.169 0.000) 3.357 2.250 ( 1.682 -0.306 0.000) 1.710 2.276 ( -0.472 -0.779 0.000) 0.911 2.338 ( -1.684 1.472 0.000) 2.236 2.374 ( -2.267 -2.291 0.000) 3.223 2.423 ( -2.158 2.003 0.000) 2.944 2.439 ( -1.725 -0.066 0.000) 1.726 2.558 ( -0.591 2.246 0.000) 2.323 2.613 ( -0.929 3.483 0.000) 3.605 2.626 ( -1.579 -4.100 0.000) 4.394 2.677 ( -0.056 -3.165 0.000) 3.166 5.473 ( -0.425 0.692 0.000) 0.812 5.488 ( -0.648 -0.775 0.000) 1.010 5.495 ( -0.336 0.285 0.000) 0.441 5.501 ( -0.290 -0.252 0.000) 0.384 6.319 ( 1.146 0.638 0.000) 1.312 6.332 ( 1.052 -0.573 0.000) 1.198 6.408 ( -0.570 4.160 0.000) 4.199 6.477 ( -0.198 1.148 0.000) 1.165 6.489 ( 0.821 -3.462 0.000) 3.558 6.504 ( -0.467 -1.468 0.000) 1.541 6.523 ( -0.623 1.449 0.000) 1.577 6.577 ( -1.334 -3.034 0.000) 3.315 6.910 ( 1.638 2.518 0.000) 3.004 6.959 ( 2.283 -2.350 0.000) 3.277 7.011 ( -0.495 3.128 0.000) 3.167 7.072 ( 1.464 -3.318 0.000) 3.626 ======================= Grid point 16 (11/36) ======================= q-point: ( 0.33 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-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.525 ( 5.306 6.729 0.000) 8.569 0.571 ( 6.840 3.884 0.000) 7.865 0.611 ( 3.397 -1.444 0.000) 3.691 0.618 ( 2.126 -0.337 0.000) 2.153 0.648 ( -1.294 2.818 0.000) 3.101 0.663 ( 1.998 -4.416 0.000) 4.848 0.676 ( 1.749 -5.449 0.000) 5.723 0.722 ( -2.000 -3.122 0.000) 3.708 0.804 ( -1.347 -0.889 0.000) 1.614 0.810 ( -1.366 -1.740 0.000) 2.212 0.857 ( -4.123 0.353 0.000) 4.139 0.898 ( -3.043 -3.230 0.000) 4.437 1.081 ( -0.389 1.220 0.000) 1.281 1.089 ( 4.870 4.059 0.000) 6.339 1.100 ( -1.134 -0.514 0.000) 1.245 1.147 ( 2.034 -1.784 0.000) 2.706 1.205 ( -2.858 1.731 0.000) 3.341 1.226 ( -0.804 -0.835 0.000) 1.159 1.233 ( 3.126 2.134 0.000) 3.785 1.244 ( 2.860 1.162 0.000) 3.087 1.305 ( 4.417 2.652 0.000) 5.152 1.351 ( 2.234 -1.375 0.000) 2.623 1.359 ( 3.523 3.351 0.000) 4.862 1.411 ( 3.229 -1.544 0.000) 3.579 1.427 ( -0.273 0.192 0.000) 0.334 1.434 ( -0.416 3.758 0.000) 3.781 1.443 ( 1.866 -1.173 0.000) 2.205 1.478 ( 1.068 -0.057 0.000) 1.070 1.524 ( -4.344 0.773 0.000) 4.412 1.570 ( 2.344 -3.077 0.000) 3.868 1.628 ( -0.988 0.838 0.000) 1.295 1.647 ( 0.283 -0.989 0.000) 1.029 1.662 ( 0.661 1.680 0.000) 1.805 1.700 ( 1.349 -1.820 0.000) 2.265 1.728 ( 1.202 0.747 0.000) 1.415 1.742 ( -0.992 -0.775 0.000) 1.259 1.796 ( 0.430 2.861 0.000) 2.893 1.832 ( -0.116 1.317 0.000) 1.322 1.854 ( 0.363 -0.902 0.000) 0.972 1.887 ( -1.946 -5.061 0.000) 5.422 2.004 ( 1.066 2.768 0.000) 2.966 2.017 ( -0.706 0.980 0.000) 1.207 2.033 ( -0.782 -0.744 0.000) 1.079 2.064 ( 0.145 -3.181 0.000) 3.185 2.220 ( -0.132 2.212 0.000) 2.216 2.235 ( -0.184 1.282 0.000) 1.296 2.265 ( -0.578 -1.927 0.000) 2.012 2.268 ( -0.514 -1.475 0.000) 1.562 2.320 ( -0.408 0.009 0.000) 0.408 2.327 ( -1.164 -0.777 0.000) 1.399 2.358 ( -3.130 2.458 0.000) 3.980 2.395 ( -2.650 -1.109 0.000) 2.873 2.539 ( -0.995 2.229 0.000) 2.441 2.570 ( -2.573 3.963 0.000) 4.725 2.599 ( -0.512 -3.725 0.000) 3.760 2.650 ( -2.066 -4.179 0.000) 4.662 5.467 ( -0.016 0.377 0.000) 0.377 5.475 ( -0.421 -0.406 0.000) 0.585 5.485 ( -0.604 0.287 0.000) 0.669 5.490 ( -0.700 -0.245 0.000) 0.741 6.348 ( 1.425 1.228 0.000) 1.881 6.366 ( 2.064 -0.245 0.000) 2.078 6.398 ( -0.640 4.434 0.000) 4.480 6.466 ( -2.134 0.074 0.000) 2.135 6.475 ( 0.104 0.391 0.000) 0.405 6.499 ( 0.172 -1.680 0.000) 1.689 6.509 ( -1.075 -2.000 0.000) 2.271 6.544 ( -1.575 -3.022 0.000) 3.408 6.951 ( 1.962 3.388 0.000) 3.915 7.003 ( -0.444 4.229 0.000) 4.252 7.017 ( 2.790 -3.148 0.000) 4.206 7.088 ( -0.108 -4.366 0.000) 4.367 ======================= Grid point 17 (12/36) ======================= q-point: (-0.50 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.594 ( -0.727 5.060 0.000) 5.112 0.594 ( 0.727 5.060 0.000) 5.112 0.663 ( -0.249 -1.208 0.000) 1.233 0.663 ( 0.249 -1.208 0.000) 1.233 0.691 ( -0.228 -0.141 0.000) 0.268 0.691 ( 0.228 -0.141 0.000) 0.268 0.714 ( -0.000 -3.068 0.000) 3.068 0.714 ( 0.000 -3.068 0.000) 3.068 0.771 ( -0.381 0.135 0.000) 0.404 0.771 ( 0.381 0.135 0.000) 0.404 0.827 ( -2.759 -4.504 0.000) 5.282 0.827 ( 2.759 -4.504 0.000) 5.282 1.087 ( -0.926 2.273 0.000) 2.454 1.087 ( 0.926 2.273 0.000) 2.454 1.114 ( -0.590 -0.270 0.000) 0.649 1.114 ( 0.590 -0.270 0.000) 0.649 1.236 ( -3.810 3.012 0.000) 4.857 1.236 ( 3.810 3.012 0.000) 4.857 1.278 ( -1.190 -1.394 0.000) 1.833 1.278 ( 1.190 -1.394 0.000) 1.833 1.313 ( -0.048 3.230 0.000) 3.230 1.313 ( 0.048 3.230 0.000) 3.230 1.391 ( -2.274 -4.205 0.000) 4.781 1.391 ( 2.274 -4.205 0.000) 4.781 1.468 ( -0.027 0.534 0.000) 0.535 1.468 ( 0.027 0.534 0.000) 0.535 1.480 ( -0.381 -0.526 0.000) 0.649 1.480 ( 0.381 -0.526 0.000) 0.649 1.592 ( -1.509 2.890 0.000) 3.260 1.592 ( 1.509 2.890 0.000) 3.260 1.624 ( -0.435 1.073 0.000) 1.158 1.624 ( 0.435 1.073 0.000) 1.158 1.675 ( -1.200 2.188 0.000) 2.495 1.675 ( 1.200 2.188 0.000) 2.495 1.746 ( -0.826 -3.454 0.000) 3.551 1.746 ( 0.826 -3.454 0.000) 3.551 1.810 ( -0.281 2.182 0.000) 2.200 1.810 ( 0.281 2.182 0.000) 2.200 1.845 ( -0.003 -1.413 0.000) 1.413 1.845 ( 0.003 -1.413 0.000) 1.413 2.027 ( -0.119 1.006 0.000) 1.013 2.027 ( 0.119 1.006 0.000) 1.013 2.044 ( -1.248 -0.452 0.000) 1.327 2.044 ( 1.248 -0.452 0.000) 1.327 2.234 ( -1.049 -0.079 0.000) 1.052 2.234 ( 1.049 -0.079 0.000) 1.052 2.247 ( -0.108 -1.218 0.000) 1.223 2.247 ( 0.108 -1.218 0.000) 1.223 2.319 ( -0.176 0.753 0.000) 0.773 2.319 ( 0.176 0.753 0.000) 0.773 2.338 ( -1.970 -1.211 0.000) 2.312 2.338 ( 1.970 -1.211 0.000) 2.312 2.531 ( -0.084 3.506 0.000) 3.507 2.531 ( 0.084 3.506 0.000) 3.507 2.609 ( -0.621 -4.271 0.000) 4.316 2.609 ( 0.621 -4.271 0.000) 4.316 5.474 ( -0.099 0.053 0.000) 0.112 5.474 ( 0.099 0.053 0.000) 0.112 5.475 ( -0.215 -0.037 0.000) 0.218 5.475 ( 0.215 -0.037 0.000) 0.218 6.377 ( -0.396 2.818 0.000) 2.845 6.377 ( 0.396 2.818 0.000) 2.845 6.415 ( -1.014 -0.036 0.000) 1.014 6.415 ( 1.014 -0.036 0.000) 1.014 6.485 ( -0.763 -0.642 0.000) 0.997 6.485 ( 0.763 -0.642 0.000) 0.997 6.514 ( -0.170 -2.474 0.000) 2.480 6.514 ( 0.170 -2.474 0.000) 2.480 6.986 ( -0.266 4.135 0.000) 4.144 6.986 ( 0.266 4.135 0.000) 4.144 7.067 ( -0.509 -3.967 0.000) 4.000 7.067 ( 0.509 -3.967 0.000) 4.000 ======================= Grid point 35 (13/36) ======================= q-point: ( 0.00 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.157 ( 0.000 0.000 7.521) 7.521 0.175 ( 0.000 0.000 9.219) 9.219 0.354 ( 0.000 0.000 17.506) 17.506 0.511 ( 0.000 0.000 -4.756) 4.756 0.557 ( 0.000 0.000 -3.758) 3.758 0.624 ( 0.000 0.000 -4.387) 4.387 0.670 ( 0.000 0.000 2.397) 2.397 0.786 ( 0.000 0.000 -1.830) 1.830 0.820 ( 0.000 0.000 -3.763) 3.763 0.829 ( 0.000 0.000 4.018) 4.018 0.888 ( 0.000 0.000 0.245) 0.245 0.920 ( 0.000 0.000 1.402) 1.402 1.002 ( 0.000 0.000 2.674) 2.674 1.041 ( 0.000 0.000 -2.733) 2.733 1.146 ( 0.000 0.000 -1.710) 1.710 1.207 ( 0.000 0.000 4.399) 4.399 1.243 ( 0.000 0.000 -5.771) 5.771 1.246 ( 0.000 0.000 -1.028) 1.028 1.259 ( 0.000 0.000 -2.558) 2.558 1.267 ( 0.000 0.000 0.148) 0.148 1.311 ( 0.000 0.000 0.564) 0.564 1.314 ( 0.000 0.000 0.167) 0.167 1.376 ( 0.000 0.000 -0.799) 0.799 1.387 ( 0.000 0.000 -0.536) 0.536 1.395 ( 0.000 0.000 1.665) 1.665 1.397 ( 0.000 0.000 -0.190) 0.190 1.405 ( 0.000 0.000 1.319) 1.319 1.446 ( 0.000 0.000 -2.513) 2.513 1.503 ( 0.000 0.000 -3.093) 3.093 1.519 ( 0.000 0.000 4.692) 4.692 1.536 ( 0.000 0.000 1.777) 1.777 1.590 ( 0.000 0.000 -1.777) 1.777 1.617 ( 0.000 0.000 1.638) 1.638 1.671 ( 0.000 0.000 -1.662) 1.662 1.733 ( 0.000 0.000 -2.409) 2.409 1.751 ( 0.000 0.000 1.867) 1.867 1.817 ( 0.000 0.000 -2.021) 2.021 1.834 ( 0.000 0.000 -1.796) 1.796 1.853 ( 0.000 0.000 2.464) 2.464 1.917 ( 0.000 0.000 0.146) 0.146 1.985 ( 0.000 0.000 2.144) 2.144 1.989 ( 0.000 0.000 0.976) 0.976 2.119 ( 0.000 0.000 2.444) 2.444 2.135 ( 0.000 0.000 -0.334) 0.334 2.141 ( 0.000 0.000 -1.175) 1.175 2.223 ( 0.000 0.000 1.556) 1.556 2.261 ( 0.000 0.000 1.368) 1.368 2.290 ( 0.000 0.000 -1.020) 1.020 2.320 ( 0.000 0.000 -0.031) 0.031 2.343 ( 0.000 0.000 -1.291) 1.291 2.374 ( 0.000 0.000 -1.863) 1.863 2.428 ( 0.000 0.000 -3.923) 3.923 2.531 ( 0.000 0.000 3.662) 3.662 2.596 ( 0.000 0.000 3.633) 3.633 2.729 ( 0.000 0.000 -1.712) 1.712 2.742 ( 0.000 0.000 -1.090) 1.090 5.469 ( 0.000 0.000 0.599) 0.599 5.496 ( 0.000 0.000 0.269) 0.269 5.507 ( 0.000 0.000 -0.213) 0.213 5.515 ( 0.000 0.000 -0.272) 0.272 6.247 ( 0.000 0.000 0.268) 0.268 6.318 ( 0.000 0.000 2.155) 2.155 6.325 ( 0.000 0.000 -1.522) 1.522 6.438 ( 0.000 0.000 -0.594) 0.594 6.475 ( 0.000 0.000 0.373) 0.373 6.532 ( 0.000 0.000 -1.300) 1.300 6.693 ( 0.000 0.000 1.769) 1.769 6.796 ( 0.000 0.000 -0.969) 0.969 6.846 ( 0.000 0.000 0.609) 0.609 6.874 ( 0.000 0.000 0.302) 0.302 6.970 ( 0.000 0.000 -0.711) 0.711 6.986 ( 0.000 0.000 -0.864) 0.864 ======================= Grid point 36 (14/36) ======================= q-point: ( 0.17 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 45 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.292 ( 4.780 0.000 7.880) 9.216 0.309 ( 9.928 0.000 3.216) 10.436 0.521 ( 7.203 0.000 0.168) 7.205 0.569 ( 0.945 0.000 -3.289) 3.422 0.582 ( 10.415 0.000 2.844) 10.797 0.687 ( 2.087 0.000 -3.760) 4.300 0.723 ( 5.143 0.000 2.350) 5.655 0.801 ( 1.640 0.000 -2.531) 3.016 0.821 ( 3.770 0.000 0.683) 3.832 0.854 ( 0.599 0.000 0.114) 0.610 0.887 ( 0.349 0.000 -0.054) 0.353 0.951 ( 1.690 0.000 0.342) 1.724 0.987 ( -1.461 0.000 2.199) 2.640 1.009 ( -3.064 0.000 -0.948) 3.208 1.129 ( -2.732 0.000 -4.376) 5.159 1.136 ( -3.150 0.000 0.433) 3.180 1.160 ( -3.951 0.000 -1.683) 4.294 1.177 ( -5.357 0.000 1.731) 5.630 1.201 ( -1.824 0.000 -3.180) 3.666 1.284 ( -2.913 0.000 1.398) 3.231 1.309 ( -2.533 0.000 7.399) 7.821 1.311 ( -1.441 0.000 -1.811) 2.314 1.326 ( 1.976 0.000 1.071) 2.248 1.336 ( -3.280 0.000 0.194) 3.285 1.370 ( -0.512 0.000 -2.294) 2.350 1.425 ( 1.882 0.000 2.072) 2.799 1.449 ( 2.653 0.000 -2.101) 3.384 1.463 ( 0.694 0.000 -0.571) 0.899 1.481 ( 5.540 0.000 0.851) 5.605 1.554 ( 3.231 0.000 -0.789) 3.326 1.557 ( 0.033 0.000 2.735) 2.735 1.635 ( 3.009 0.000 -3.190) 4.385 1.647 ( 1.630 0.000 1.716) 2.367 1.683 ( 1.051 0.000 -1.240) 1.626 1.728 ( 0.990 0.000 -1.107) 1.485 1.814 ( 0.603 0.000 -1.836) 1.933 1.825 ( -2.103 0.000 2.415) 3.202 1.871 ( 1.660 0.000 0.598) 1.764 1.894 ( -2.007 0.000 0.273) 2.025 1.958 ( -2.184 0.000 1.709) 2.773 1.975 ( 1.932 0.000 -1.933) 2.733 2.006 ( 0.761 0.000 -2.546) 2.658 2.091 ( -2.372 0.000 1.780) 2.966 2.146 ( 0.929 0.000 0.209) 0.952 2.200 ( 1.732 0.000 1.793) 2.494 2.221 ( -2.919 0.000 0.077) 2.920 2.222 ( 2.963 0.000 -1.925) 3.534 2.286 ( -0.200 0.000 -1.593) 1.605 2.290 ( -1.050 0.000 -1.090) 1.513 2.372 ( -0.620 0.000 -0.823) 1.030 2.411 ( 2.870 0.000 -0.685) 2.951 2.435 ( 0.297 0.000 -2.249) 2.269 2.493 ( -3.154 0.000 4.438) 5.445 2.560 ( -2.949 0.000 2.301) 3.741 2.710 ( -1.559 0.000 -1.413) 2.105 2.750 ( 0.080 0.000 -3.859) 3.860 5.465 ( -0.209 0.000 0.641) 0.674 5.490 ( -0.498 0.000 0.164) 0.524 5.503 ( -0.690 0.000 -0.240) 0.731 5.506 ( -0.474 0.000 -0.211) 0.519 6.257 ( 0.906 0.000 0.239) 0.937 6.334 ( 1.317 0.000 2.264) 2.619 6.337 ( 0.232 0.000 -0.737) 0.773 6.441 ( 0.336 0.000 -0.598) 0.686 6.465 ( -0.762 0.000 0.292) 0.816 6.527 ( -0.348 0.000 -1.393) 1.436 6.663 ( -2.459 0.000 0.838) 2.598 6.692 ( -3.352 0.000 -0.918) 3.476 6.863 ( 0.202 0.000 1.026) 1.046 6.876 ( 0.162 0.000 0.266) 0.311 7.009 ( 3.248 0.000 -0.239) 3.257 7.137 ( 4.924 0.000 -6.554) 8.198 ======================= Grid point 37 (15/36) ======================= q-point: ( 0.33 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-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.388 ( 4.386 0.000 4.012) 5.945 0.483 ( 4.401 0.000 0.289) 4.410 0.600 ( 2.631 0.000 -1.222) 2.900 0.614 ( 3.879 0.000 -3.670) 5.340 0.684 ( -0.211 0.000 -0.466) 0.512 0.713 ( 0.325 0.000 -0.804) 0.867 0.829 ( 4.067 0.000 1.429) 4.311 0.849 ( -0.368 0.000 0.068) 0.374 0.880 ( 1.041 0.000 1.393) 1.739 0.904 ( -0.312 0.000 -4.176) 4.188 0.905 ( -1.207 0.000 -4.868) 5.015 0.939 ( -2.022 0.000 0.090) 2.024 0.967 ( -1.877 0.000 1.927) 2.690 1.003 ( 4.046 0.000 0.251) 4.053 1.039 ( -4.691 0.000 -1.742) 5.004 1.113 ( 6.525 0.000 1.191) 6.633 1.124 ( -0.622 0.000 2.170) 2.257 1.154 ( 2.069 0.000 0.146) 2.074 1.203 ( -3.050 0.000 -1.695) 3.489 1.223 ( 3.111 0.000 -1.544) 3.473 1.233 ( -1.778 0.000 0.845) 1.969 1.282 ( -0.356 0.000 -3.247) 3.267 1.306 ( 0.215 0.000 1.084) 1.105 1.351 ( -0.523 0.000 0.060) 0.526 1.406 ( 4.188 0.000 1.449) 4.431 1.459 ( 0.115 0.000 -3.861) 3.863 1.465 ( 0.818 0.000 -0.050) 0.819 1.529 ( 4.928 0.000 2.801) 5.668 1.551 ( -0.851 0.000 2.012) 2.184 1.585 ( 2.407 0.000 -0.032) 2.408 1.616 ( 1.570 0.000 0.853) 1.786 1.647 ( -1.898 0.000 -2.323) 3.000 1.681 ( -4.331 0.000 2.237) 4.875 1.704 ( 2.550 0.000 2.073) 3.286 1.715 ( 1.812 0.000 -0.012) 1.812 1.781 ( -6.334 0.000 2.176) 6.697 1.826 ( 2.320 0.000 0.126) 2.323 1.835 ( -1.543 0.000 -0.120) 1.547 1.873 ( -0.490 0.000 0.587) 0.764 1.933 ( -4.124 0.000 -1.339) 4.336 1.956 ( 2.208 0.000 0.712) 2.320 1.993 ( -0.700 0.000 0.512) 0.867 2.051 ( -1.648 0.000 0.589) 1.750 2.134 ( -4.388 0.000 -1.139) 4.534 2.210 ( 3.854 0.000 1.421) 4.108 2.214 ( 0.759 0.000 -0.691) 1.026 2.281 ( 0.408 0.000 -0.131) 0.429 2.288 ( -0.191 0.000 -1.098) 1.115 2.310 ( 3.691 0.000 0.988) 3.821 2.332 ( -2.831 0.000 0.957) 2.989 2.402 ( -3.168 0.000 0.882) 3.288 2.432 ( -0.420 0.000 -0.348) 0.545 2.451 ( -0.264 0.000 2.781) 2.794 2.495 ( -2.891 0.000 0.459) 2.927 2.682 ( -0.546 0.000 -1.747) 1.830 2.734 ( -1.498 0.000 -3.732) 4.022 5.465 ( 0.229 0.000 0.502) 0.552 5.478 ( -0.481 0.000 0.079) 0.488 5.483 ( -0.854 0.000 -0.122) 0.862 5.495 ( -0.708 0.000 -0.106) 0.716 6.281 ( 1.254 0.000 0.261) 1.281 6.329 ( -0.800 0.000 0.407) 0.897 6.366 ( 1.505 0.000 2.142) 2.618 6.429 ( -1.468 0.000 0.639) 1.601 6.473 ( 1.535 0.000 -1.667) 2.266 6.524 ( -0.000 0.000 -2.052) 2.052 6.607 ( -2.308 0.000 0.234) 2.320 6.634 ( -2.502 0.000 -1.140) 2.749 6.869 ( 0.390 0.000 0.451) 0.596 6.879 ( 0.048 0.000 0.112) 0.122 7.087 ( 3.668 0.000 -0.135) 3.670 7.179 ( -0.049 0.000 -3.088) 3.089 ======================= Grid point 38 (16/36) ======================= q-point: (-0.50 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.473 ( -1.554 0.000 1.258) 1.999 0.473 ( 1.554 0.000 1.258) 1.999 0.663 ( -0.513 0.000 -1.976) 2.042 0.663 ( 0.513 0.000 -1.976) 2.042 0.696 ( -0.109 0.000 -1.162) 1.167 0.696 ( 0.109 0.000 -1.162) 1.167 0.859 ( -0.450 0.000 -3.786) 3.812 0.859 ( 0.450 0.000 -3.786) 3.812 0.869 ( -0.038 0.000 0.634) 0.635 0.869 ( 0.038 0.000 0.634) 0.635 0.921 ( -0.015 0.000 -2.338) 2.338 0.921 ( 0.015 0.000 -2.338) 2.338 0.949 ( -2.343 0.000 1.371) 2.715 0.949 ( 2.343 0.000 1.371) 2.715 1.071 ( -1.399 0.000 3.339) 3.621 1.071 ( 1.399 0.000 3.339) 3.621 1.158 ( -0.352 0.000 -1.064) 1.121 1.158 ( 0.352 0.000 -1.064) 1.121 1.204 ( -1.980 0.000 -2.519) 3.204 1.204 ( 1.980 0.000 -2.519) 3.204 1.307 ( -0.142 0.000 -0.397) 0.422 1.307 ( 0.142 0.000 -0.397) 0.422 1.400 ( -5.410 0.000 -2.909) 6.142 1.400 ( 5.410 0.000 -2.909) 6.142 1.460 ( -3.189 0.000 -0.141) 3.192 1.460 ( 3.189 0.000 -0.141) 3.192 1.506 ( -0.309 0.000 0.148) 0.343 1.506 ( 0.309 0.000 0.148) 0.343 1.534 ( -4.916 0.000 1.117) 5.041 1.534 ( 4.916 0.000 1.117) 5.041 1.614 ( -0.210 0.000 1.369) 1.385 1.614 ( 0.210 0.000 1.369) 1.385 1.673 ( -1.868 0.000 4.254) 4.646 1.673 ( 1.868 0.000 4.254) 4.646 1.752 ( -0.385 0.000 1.260) 1.317 1.752 ( 0.385 0.000 1.260) 1.317 1.813 ( -1.912 0.000 0.003) 1.912 1.813 ( 1.912 0.000 0.003) 1.912 1.892 ( -0.674 0.000 0.498) 0.837 1.892 ( 0.674 0.000 0.498) 0.837 2.009 ( -2.335 0.000 0.230) 2.347 2.009 ( 2.335 0.000 0.230) 2.347 2.025 ( -4.440 0.000 -0.499) 4.468 2.025 ( 4.440 0.000 -0.499) 4.468 2.246 ( -2.304 0.000 -1.371) 2.681 2.246 ( 2.304 0.000 -1.371) 2.681 2.294 ( -0.467 0.000 0.984) 1.090 2.294 ( 0.467 0.000 0.984) 1.090 2.305 ( -2.239 0.000 0.723) 2.353 2.305 ( 2.239 0.000 0.723) 2.353 2.406 ( -1.521 0.000 0.297) 1.550 2.406 ( 1.521 0.000 0.297) 1.550 2.450 ( -0.052 0.000 2.081) 2.082 2.450 ( 0.052 0.000 2.081) 2.082 2.696 ( -1.512 0.000 -2.860) 3.235 2.696 ( 1.512 0.000 -2.860) 3.235 5.472 ( -0.040 0.000 0.031) 0.051 5.472 ( 0.040 0.000 0.031) 0.051 5.478 ( -0.614 0.000 0.155) 0.634 5.478 ( 0.614 0.000 0.155) 0.634 6.308 ( -0.548 0.000 0.424) 0.693 6.308 ( 0.548 0.000 0.424) 0.693 6.397 ( -0.945 0.000 1.481) 1.757 6.397 ( 0.945 0.000 1.481) 1.757 6.511 ( -1.542 0.000 -2.729) 3.134 6.511 ( 1.542 0.000 -2.729) 3.134 6.591 ( -0.283 0.000 0.180) 0.335 6.591 ( 0.283 0.000 0.180) 0.335 6.877 ( -0.025 0.000 0.104) 0.107 6.877 ( 0.025 0.000 0.104) 0.107 7.153 ( -1.529 0.000 -1.036) 1.847 7.153 ( 1.529 0.000 -1.036) 1.847 ======================= Grid point 42 (17/36) ======================= q-point: ( 0.00 0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 52 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.196 ( 0.000 5.433 4.779) 7.236 0.289 ( 0.000 8.338 4.043) 9.267 0.462 ( 0.000 6.990 2.539) 7.437 0.528 ( 0.000 1.294 -3.819) 4.032 0.614 ( 0.000 6.413 2.597) 6.919 0.658 ( 0.000 3.278 -3.822) 5.035 0.669 ( 0.000 0.156 4.198) 4.201 0.735 ( 0.000 -2.713 -5.435) 6.074 0.788 ( 0.000 -3.427 2.084) 4.011 0.857 ( 0.000 1.670 -2.621) 3.108 0.897 ( 0.000 -2.114 0.303) 2.136 0.916 ( 0.000 0.364 -1.532) 1.575 0.992 ( 0.000 -1.867 0.963) 2.101 1.002 ( 0.000 -0.262 2.866) 2.878 1.089 ( 0.000 -4.434 -3.781) 5.827 1.144 ( 0.000 -4.894 -2.045) 5.304 1.175 ( 0.000 -9.042 0.948) 9.091 1.202 ( 0.000 -2.592 3.389) 4.266 1.243 ( 0.000 -0.333 2.903) 2.922 1.258 ( 0.000 -1.512 -1.799) 2.350 1.286 ( 0.000 1.027 -1.453) 1.779 1.309 ( 0.000 -0.395 -1.753) 1.797 1.357 ( 0.000 -1.023 5.260) 5.359 1.394 ( 0.000 -0.694 0.468) 0.837 1.396 ( 0.000 0.194 2.267) 2.275 1.423 ( 0.000 3.099 -0.184) 3.104 1.434 ( 0.000 0.961 -1.486) 1.769 1.454 ( 0.000 -1.226 -1.656) 2.061 1.496 ( 0.000 3.374 -1.617) 3.741 1.506 ( 0.000 2.870 1.477) 3.228 1.565 ( 0.000 0.024 -4.677) 4.677 1.601 ( 0.000 -0.759 0.909) 1.184 1.694 ( 0.000 2.509 -3.613) 4.398 1.726 ( 0.000 2.987 0.796) 3.092 1.737 ( 0.000 -0.694 0.421) 0.812 1.743 ( 0.000 1.064 -0.959) 1.433 1.796 ( 0.000 -1.609 -1.817) 2.427 1.822 ( 0.000 -0.809 0.982) 1.273 1.852 ( 0.000 -2.705 0.413) 2.736 1.880 ( 0.000 0.069 -0.048) 0.084 2.012 ( 0.000 2.517 -1.965) 3.193 2.059 ( 0.000 2.230 -0.924) 2.414 2.091 ( 0.000 -2.126 1.170) 2.427 2.148 ( 0.000 -1.625 -1.482) 2.199 2.168 ( 0.000 -0.383 2.165) 2.199 2.175 ( 0.000 3.967 1.013) 4.095 2.256 ( 0.000 -0.627 1.423) 1.555 2.304 ( 0.000 -0.470 0.240) 0.528 2.336 ( 0.000 1.244 0.272) 1.273 2.360 ( 0.000 1.282 -1.959) 2.341 2.424 ( 0.000 0.017 -2.168) 2.168 2.426 ( 0.000 3.178 -2.548) 4.073 2.568 ( 0.000 3.129 3.273) 4.528 2.591 ( 0.000 -0.229 2.999) 3.007 2.666 ( 0.000 -2.754 -0.938) 2.909 2.766 ( 0.000 -0.840 0.415) 0.937 5.472 ( 0.000 0.339 0.516) 0.617 5.497 ( 0.000 0.090 0.204) 0.223 5.505 ( 0.000 -0.148 -0.226) 0.270 5.510 ( 0.000 -0.474 -0.157) 0.500 6.290 ( 0.000 1.881 -1.785) 2.593 6.322 ( 0.000 -0.298 -1.270) 1.304 6.331 ( 0.000 1.823 0.241) 1.839 6.432 ( 0.000 -0.692 -0.942) 1.169 6.490 ( 0.000 1.527 0.184) 1.538 6.530 ( 0.000 -0.163 -0.665) 0.685 6.633 ( 0.000 -6.192 5.876) 8.536 6.664 ( 0.000 -2.664 1.405) 3.012 6.858 ( 0.000 1.196 0.382) 1.255 6.942 ( 0.000 -0.227 -2.276) 2.288 6.998 ( 0.000 1.048 -0.660) 1.239 7.010 ( 0.000 2.462 -2.169) 3.281 ======================= Grid point 43 (18/36) ======================= q-point: ( 0.17 0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-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.294 ( 7.380 2.009 3.822) 8.550 0.363 ( 6.029 4.793 4.165) 8.756 0.525 ( 3.296 -0.117 -2.734) 4.284 0.550 ( 1.369 0.353 -3.641) 3.906 0.666 ( 2.984 4.334 2.170) 5.691 0.678 ( 2.379 -1.266 -0.571) 2.754 0.751 ( 3.479 0.804 1.568) 3.900 0.765 ( 2.934 -1.379 -2.803) 4.285 0.807 ( 2.698 -1.582 -0.042) 3.128 0.843 ( -0.261 -0.029 -1.104) 1.135 0.906 ( 1.833 0.147 1.319) 2.263 0.925 ( -0.154 -0.675 -3.230) 3.303 0.962 ( -0.202 -0.755 2.126) 2.265 1.015 ( -0.592 -1.363 0.827) 1.701 1.078 ( -3.020 -4.533 -0.777) 5.502 1.090 ( 1.738 -0.305 -1.501) 2.317 1.151 ( -1.756 -0.695 -0.976) 2.126 1.191 ( -3.001 2.167 -1.176) 3.884 1.218 ( -0.730 -2.226 1.036) 2.562 1.278 ( -1.944 -1.822 0.990) 2.842 1.289 ( 0.077 -1.477 -0.877) 1.719 1.304 ( -2.304 1.710 1.497) 3.236 1.321 ( 0.007 -0.921 -0.121) 0.929 1.359 ( -1.584 1.763 0.141) 2.374 1.398 ( 3.210 0.510 0.515) 3.290 1.433 ( 0.314 0.786 1.113) 1.398 1.450 ( 0.341 0.282 -0.690) 0.820 1.463 ( 0.724 1.984 0.600) 2.195 1.518 ( 3.536 -1.803 -0.288) 3.980 1.529 ( 0.112 3.052 -1.526) 3.414 1.608 ( 0.192 3.119 0.257) 3.136 1.618 ( 1.919 -0.644 -1.545) 2.547 1.655 ( 2.649 -0.031 0.738) 2.750 1.714 ( -0.008 2.143 -0.524) 2.206 1.731 ( -1.360 -0.074 0.173) 1.373 1.777 ( 0.375 -2.017 0.082) 2.053 1.793 ( 0.805 -1.629 -0.604) 1.915 1.848 ( 0.161 -2.192 0.211) 2.208 1.864 ( 0.605 -2.107 1.510) 2.662 1.945 ( 1.975 -1.457 0.379) 2.483 1.985 ( -0.858 -0.666 -1.866) 2.159 2.032 ( -1.992 2.260 -0.920) 3.150 2.090 ( -0.820 -1.099 0.684) 1.532 2.153 ( -1.560 -2.033 1.114) 2.794 2.175 ( 0.168 1.707 0.930) 1.951 2.222 ( -0.119 1.069 -0.437) 1.161 2.227 ( 2.977 0.677 0.710) 3.135 2.291 ( -0.381 1.258 -0.848) 1.564 2.326 ( 0.773 -1.588 1.315) 2.202 2.363 ( -2.688 2.733 -1.145) 4.001 2.408 ( -1.542 -0.874 -1.355) 2.231 2.426 ( 0.431 -0.013 -1.067) 1.151 2.546 ( -3.018 1.974 1.341) 3.848 2.567 ( -0.941 2.209 2.870) 3.742 2.667 ( -0.528 -3.060 -1.848) 3.614 2.747 ( -1.499 -1.165 -0.760) 2.045 5.469 ( -0.253 0.296 0.576) 0.695 5.490 ( -0.627 -0.037 -0.059) 0.631 5.501 ( -0.688 -0.179 -0.017) 0.711 5.503 ( -0.239 -0.188 -0.178) 0.352 6.288 ( 0.226 2.223 -1.261) 2.566 6.324 ( 0.648 -0.287 0.550) 0.897 6.356 ( 1.328 1.255 0.668) 1.945 6.442 ( 0.435 -0.078 -1.102) 1.187 6.481 ( -0.776 1.507 0.206) 1.707 6.518 ( -0.812 -0.677 -0.639) 1.235 6.614 ( -1.680 -3.961 0.524) 4.334 6.655 ( -1.028 -2.672 1.829) 3.397 6.864 ( 0.507 1.207 -0.122) 1.315 6.939 ( -0.889 3.829 -0.149) 3.934 6.987 ( 2.855 -1.939 -1.084) 3.617 7.109 ( 3.774 -2.203 -3.228) 5.433 ======================= Grid point 44 (19/36) ======================= q-point: ( 0.33 0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 76 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.414 ( 3.973 3.080 1.522) 5.253 0.521 ( 5.851 2.958 0.773) 6.602 0.587 ( 3.169 -0.852 -1.407) 3.570 0.601 ( 3.740 -1.194 -2.089) 4.447 0.661 ( -0.884 -0.531 0.258) 1.063 0.737 ( 0.207 0.665 -1.461) 1.619 0.777 ( 0.928 -3.890 1.433) 4.248 0.790 ( -0.168 -5.155 -0.179) 5.160 0.817 ( 0.050 -3.603 -1.735) 3.999 0.848 ( -0.063 -2.801 -1.517) 3.186 0.908 ( -1.102 -1.417 -1.334) 2.237 0.931 ( -2.110 -2.300 1.193) 3.342 0.995 ( 0.869 3.361 -1.137) 3.653 1.035 ( -1.953 -0.877 -1.338) 2.524 1.055 ( -2.780 2.036 -3.358) 4.811 1.110 ( 4.429 1.806 -0.228) 4.788 1.161 ( 1.048 1.021 0.489) 1.543 1.191 ( 0.496 2.293 -0.000) 2.346 1.214 ( 4.740 1.436 -0.456) 4.974 1.227 ( -0.367 1.736 -0.708) 1.910 1.267 ( -1.227 1.890 -0.149) 2.258 1.298 ( -0.506 0.504 -0.572) 0.915 1.308 ( 3.296 2.267 -0.265) 4.009 1.359 ( -0.924 0.279 -0.521) 1.097 1.402 ( 0.146 -0.552 -1.675) 1.770 1.445 ( 2.539 -0.103 0.259) 2.554 1.493 ( 0.944 0.960 -1.032) 1.696 1.502 ( 2.198 -0.861 1.914) 3.039 1.549 ( 0.632 -0.101 2.934) 3.003 1.583 ( 2.022 -1.426 -0.662) 2.561 1.634 ( 1.957 1.867 1.409) 3.050 1.644 ( 0.665 0.066 -1.150) 1.330 1.667 ( -3.770 -0.722 1.238) 4.033 1.716 ( 2.303 0.364 1.399) 2.719 1.732 ( 1.044 1.308 1.261) 2.096 1.782 ( -1.814 -0.338 1.738) 2.535 1.812 ( -1.350 -1.683 0.995) 2.376 1.819 ( -0.247 0.008 0.353) 0.431 1.863 ( 0.267 -0.719 -0.313) 0.829 1.938 ( -3.552 -0.754 -1.108) 3.796 1.967 ( 1.243 1.199 0.686) 1.858 2.003 ( -0.590 1.091 0.481) 1.331 2.053 ( -1.666 0.054 0.809) 1.853 2.113 ( -2.334 -1.919 -0.790) 3.123 2.189 ( 1.383 -0.282 0.293) 1.442 2.240 ( 2.140 1.000 0.489) 2.412 2.284 ( 0.483 -1.714 0.982) 2.034 2.293 ( 0.255 0.958 0.176) 1.007 2.307 ( -1.606 1.542 -0.735) 2.344 2.334 ( -0.185 0.069 0.914) 0.935 2.367 ( -1.615 -3.003 0.906) 3.528 2.412 ( -1.395 -0.858 0.159) 1.646 2.478 ( -2.799 2.133 -0.381) 3.540 2.536 ( -1.686 2.873 2.055) 3.914 2.650 ( -0.385 -2.964 -1.925) 3.555 2.710 ( -1.971 -2.210 -1.978) 3.561 5.466 ( 0.127 0.123 0.440) 0.474 5.474 ( -0.690 -0.245 -0.289) 0.787 5.487 ( -0.583 0.214 0.209) 0.655 5.494 ( -0.636 -0.103 -0.030) 0.645 6.299 ( 0.625 1.633 -0.906) 1.969 6.338 ( 0.540 0.803 0.743) 1.220 6.382 ( 0.926 1.412 1.795) 2.464 6.440 ( -0.743 0.702 0.034) 1.023 6.473 ( 0.517 0.364 -0.737) 0.971 6.506 ( -0.240 -1.516 -1.668) 2.267 6.569 ( -2.324 -2.937 -0.636) 3.799 6.624 ( -1.635 -1.486 1.141) 2.487 6.877 ( 0.644 1.685 -0.953) 2.040 6.931 ( -0.102 3.783 0.547) 3.824 7.057 ( 3.231 -2.770 -1.120) 4.401 7.151 ( 0.435 -2.567 -1.673) 3.095 ======================= Grid point 45 (20/36) ======================= q-point: (-0.50 0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 50 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.461 ( 0.000 2.280 -0.655) 2.372 0.558 ( 0.000 5.047 2.558) 5.658 0.634 ( 0.000 -1.945 -2.006) 2.794 0.670 ( -0.000 0.204 -1.596) 1.609 0.670 ( 0.000 -0.971 -2.021) 2.242 0.716 ( -0.000 0.578 -0.045) 0.579 0.779 ( -0.000 -4.025 0.796) 4.103 0.785 ( 0.000 -4.048 -2.039) 4.533 0.818 ( 0.000 -4.876 -0.664) 4.921 0.884 ( -0.000 -4.178 3.090) 5.196 0.892 ( -0.000 -2.020 1.757) 2.677 0.899 ( 0.000 -1.373 -3.155) 3.441 0.969 ( 0.000 2.050 -4.011) 4.505 0.984 ( 0.000 2.647 -2.652) 3.747 1.042 ( -0.000 0.312 -0.742) 0.805 1.138 ( 0.000 -1.357 -2.071) 2.477 1.170 ( -0.000 5.043 3.739) 6.278 1.176 ( 0.000 0.736 -2.547) 2.651 1.216 ( -0.000 1.965 -2.495) 3.176 1.229 ( 0.000 1.653 -2.350) 2.874 1.286 ( 0.000 -0.867 -1.076) 1.382 1.328 ( 0.000 0.683 1.618) 1.757 1.376 ( -0.000 -2.055 -4.365) 4.824 1.407 ( -0.000 -0.676 -0.694) 0.969 1.469 ( 0.000 1.783 0.064) 1.784 1.479 ( -0.000 0.969 -1.154) 1.506 1.512 ( -0.000 0.611 1.433) 1.558 1.517 ( 0.000 1.335 2.044) 2.441 1.537 ( -0.000 1.330 -1.217) 1.803 1.568 ( -0.000 0.144 2.009) 2.015 1.573 ( 0.000 1.705 0.903) 1.929 1.665 ( 0.000 -1.187 1.290) 1.753 1.671 ( -0.000 2.225 2.629) 3.444 1.681 ( 0.000 1.196 5.388) 5.519 1.762 ( 0.000 1.683 0.766) 1.849 1.780 ( 0.000 1.240 1.100) 1.657 1.803 ( -0.000 -0.956 0.180) 0.973 1.807 ( -0.000 -0.337 1.053) 1.106 1.856 ( 0.000 -1.922 -0.494) 1.985 1.893 ( 0.000 -1.020 1.125) 1.519 2.009 ( -0.000 0.323 0.181) 0.370 2.017 ( 0.000 0.858 0.183) 0.878 2.023 ( -0.000 0.442 -1.063) 1.151 2.053 ( 0.000 1.329 0.468) 1.409 2.228 ( -0.000 -0.885 -1.500) 1.742 2.238 ( 0.000 -0.761 -0.300) 0.818 2.280 ( 0.000 -1.197 1.592) 1.992 2.306 ( -0.000 0.189 -0.093) 0.210 2.311 ( -0.000 0.251 1.759) 1.777 2.314 ( 0.000 0.651 0.523) 0.835 2.369 ( 0.000 -1.590 -0.009) 1.590 2.373 ( -0.000 -3.375 0.564) 3.421 2.448 ( 0.000 2.176 0.182) 2.184 2.514 ( 0.000 2.844 2.254) 3.630 2.658 ( 0.000 -3.286 -2.837) 4.341 2.675 ( -0.000 -2.240 -1.509) 2.700 5.466 ( 0.000 -0.229 -0.292) 0.370 5.473 ( -0.000 -0.255 -0.020) 0.256 5.478 ( 0.000 0.274 0.295) 0.403 5.483 ( 0.000 0.226 0.369) 0.432 6.307 ( 0.000 1.052 -0.792) 1.317 6.345 ( -0.000 2.388 0.999) 2.589 6.391 ( 0.000 0.150 0.598) 0.617 6.428 ( 0.000 2.180 2.158) 3.068 6.495 ( 0.000 -1.253 -1.780) 2.177 6.497 ( -0.000 -1.329 -3.004) 3.285 6.536 ( -0.000 -3.110 -0.106) 3.112 6.602 ( -0.000 -0.789 1.549) 1.739 6.885 ( 0.000 1.981 -1.313) 2.376 6.930 ( 0.000 3.773 0.736) 3.844 7.101 ( 0.000 -3.699 -2.092) 4.250 7.147 ( -0.000 -1.713 0.033) 1.714 ======================= Grid point 49 (21/36) ======================= q-point: ( 0.00 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 52 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.326 ( 0.000 6.220 0.184) 6.223 0.398 ( 0.000 0.425 -1.505) 1.564 0.568 ( 0.000 4.452 -6.811) 8.137 0.572 ( 0.000 4.563 2.539) 5.222 0.643 ( 0.000 -2.977 -2.847) 4.120 0.684 ( 0.000 -5.220 3.504) 6.287 0.744 ( 0.000 3.919 -5.142) 6.465 0.750 ( 0.000 3.737 -3.819) 5.344 0.789 ( 0.000 7.556 3.424) 8.296 0.848 ( 0.000 -2.314 -1.850) 2.963 0.870 ( 0.000 -0.590 -1.501) 1.613 0.895 ( 0.000 -4.627 0.062) 4.627 0.961 ( 0.000 2.607 -1.018) 2.798 0.967 ( 0.000 -2.347 1.097) 2.591 1.004 ( 0.000 -5.376 0.196) 5.380 1.076 ( 0.000 1.067 0.729) 1.293 1.088 ( 0.000 -2.729 0.694) 2.816 1.110 ( 0.000 -5.515 2.040) 5.881 1.204 ( 0.000 -1.754 2.200) 2.814 1.213 ( 0.000 -1.898 3.812) 4.259 1.302 ( 0.000 0.473 -2.844) 2.883 1.307 ( 0.000 -0.134 -3.341) 3.344 1.363 ( 0.000 0.816 2.563) 2.690 1.383 ( 0.000 -0.926 0.423) 1.018 1.406 ( 0.000 0.582 2.144) 2.222 1.415 ( 0.000 -0.756 1.197) 1.415 1.463 ( 0.000 1.416 -0.611) 1.542 1.477 ( 0.000 0.948 -1.327) 1.631 1.557 ( 0.000 3.853 2.019) 4.350 1.568 ( 0.000 1.099 -1.776) 2.088 1.569 ( 0.000 0.255 -3.343) 3.353 1.611 ( 0.000 2.357 -1.377) 2.730 1.668 ( 0.000 -4.070 0.912) 4.171 1.689 ( 0.000 -2.221 -3.637) 4.261 1.764 ( 0.000 -0.324 0.199) 0.380 1.766 ( 0.000 0.137 0.283) 0.314 1.798 ( 0.000 -0.414 0.577) 0.710 1.800 ( 0.000 0.613 -0.547) 0.821 1.845 ( 0.000 2.237 0.168) 2.243 1.879 ( 0.000 -0.500 0.333) 0.601 2.042 ( 0.000 0.465 -0.995) 1.098 2.054 ( 0.000 -1.210 0.622) 1.361 2.074 ( 0.000 0.247 -0.756) 0.795 2.094 ( 0.000 -2.070 -0.046) 2.071 2.199 ( 0.000 1.923 0.222) 1.935 2.234 ( 0.000 -1.416 0.721) 1.589 2.257 ( 0.000 2.903 0.459) 2.939 2.289 ( 0.000 -0.848 0.149) 0.861 2.361 ( 0.000 1.129 -0.127) 1.136 2.374 ( 0.000 -0.043 -1.664) 1.665 2.440 ( 0.000 1.320 -0.594) 1.447 2.457 ( 0.000 0.021 -0.926) 0.927 2.592 ( 0.000 0.884 2.096) 2.275 2.622 ( 0.000 -1.902 -0.408) 1.946 2.643 ( 0.000 3.973 2.533) 4.712 2.719 ( 0.000 -3.453 2.375) 4.191 5.482 ( 0.000 0.654 0.293) 0.716 5.496 ( 0.000 -0.730 -0.023) 0.731 5.499 ( 0.000 0.170 0.062) 0.181 5.503 ( 0.000 -0.202 -0.046) 0.207 6.306 ( 0.000 0.405 -0.221) 0.461 6.314 ( 0.000 -0.425 -0.665) 0.790 6.377 ( 0.000 2.148 -2.277) 3.130 6.411 ( 0.000 -1.380 -1.753) 2.231 6.527 ( 0.000 1.586 -0.124) 1.591 6.530 ( 0.000 0.221 0.699) 0.733 6.566 ( 0.000 -0.809 4.108) 4.187 6.603 ( 0.000 -3.156 1.245) 3.392 6.889 ( 0.000 1.762 -0.184) 1.772 6.922 ( 0.000 -1.493 -1.019) 1.808 7.018 ( 0.000 0.847 -0.634) 1.058 7.029 ( 0.000 -0.066 -1.131) 1.133 ======================= Grid point 50 (22/36) ======================= q-point: ( 0.17 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-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.381 ( 5.165 5.473 -0.045) 7.525 0.447 ( 4.439 0.811 -0.591) 4.551 0.533 ( 0.182 3.211 -1.923) 3.747 0.580 ( -0.955 2.363 -2.471) 3.550 0.629 ( 0.189 -3.054 -1.217) 3.293 0.662 ( -1.068 -6.392 0.908) 6.544 0.727 ( -0.104 2.320 0.545) 2.386 0.771 ( -1.039 0.386 0.226) 1.131 0.797 ( 0.358 -2.079 -2.808) 3.512 0.821 ( -1.596 -1.077 -0.965) 2.154 0.899 ( 0.672 0.034 -3.783) 3.843 0.916 ( -0.452 -1.527 0.546) 1.683 0.998 ( 3.500 2.209 0.888) 4.233 1.033 ( 1.595 0.957 1.109) 2.166 1.040 ( 6.653 5.000 1.415) 8.442 1.084 ( 1.410 -0.668 -3.877) 4.179 1.120 ( 2.037 -3.009 0.625) 3.687 1.168 ( 3.124 -1.156 2.107) 3.942 1.200 ( 1.107 -2.044 -1.276) 2.651 1.222 ( 2.258 -2.448 1.545) 3.671 1.281 ( -1.894 0.447 -2.149) 2.899 1.298 ( -0.179 -1.028 -1.444) 1.781 1.357 ( -1.945 2.694 1.267) 3.557 1.384 ( -0.892 0.353 1.183) 1.523 1.404 ( 2.081 0.371 0.651) 2.212 1.421 ( 1.802 -1.879 0.209) 2.612 1.446 ( -0.420 0.352 -0.339) 0.644 1.464 ( -1.257 -0.874 0.420) 1.588 1.561 ( 0.349 3.486 -1.023) 3.650 1.580 ( -0.675 1.605 -2.714) 3.224 1.618 ( 1.582 0.905 -0.520) 1.895 1.648 ( 0.524 0.522 -1.631) 1.791 1.655 ( 1.552 0.054 0.572) 1.655 1.687 ( -0.070 -2.980 1.095) 3.175 1.753 ( -0.949 0.131 1.529) 1.804 1.765 ( -0.021 0.675 0.072) 0.680 1.795 ( 0.292 1.875 1.151) 2.219 1.819 ( 0.729 -0.219 0.407) 0.863 1.827 ( 0.960 -1.005 0.028) 1.390 1.874 ( 1.057 -4.505 -0.228) 4.632 2.011 ( -1.168 2.142 0.260) 2.453 2.048 ( -1.149 -1.048 -0.546) 1.648 2.065 ( -1.015 0.186 0.376) 1.098 2.097 ( -1.005 -2.625 1.464) 3.169 2.221 ( 0.169 1.608 -0.470) 1.684 2.236 ( 0.397 0.127 0.224) 0.474 2.261 ( 0.830 1.865 0.826) 2.202 2.289 ( -0.070 -1.047 0.689) 1.255 2.330 ( -1.090 1.779 -0.550) 2.157 2.370 ( -1.572 -1.909 -0.358) 2.499 2.421 ( -2.642 1.566 -0.483) 3.109 2.431 ( -1.077 0.320 -0.764) 1.359 2.565 ( -2.003 1.277 0.733) 2.485 2.606 ( -1.428 -2.661 -0.889) 3.148 2.638 ( -1.039 3.494 2.097) 4.205 2.704 ( -1.578 -2.995 1.409) 3.667 5.477 ( -0.433 0.497 0.316) 0.731 5.486 ( -0.748 -0.387 -0.148) 0.855 5.498 ( -0.246 -0.019 0.206) 0.321 5.500 ( -0.273 -0.152 -0.090) 0.325 6.318 ( 0.904 0.691 -0.970) 1.495 6.322 ( 0.874 0.004 -0.168) 0.890 6.388 ( 1.015 2.030 -0.889) 2.438 6.427 ( 0.821 -1.586 -2.202) 2.835 6.508 ( -1.280 0.997 0.202) 1.635 6.527 ( -0.578 1.297 0.315) 1.454 6.551 ( -1.313 -1.361 2.497) 3.133 6.595 ( -0.694 -3.239 2.064) 3.903 6.899 ( 1.005 2.144 -0.957) 2.554 6.943 ( 1.897 -2.257 -1.492) 3.304 7.008 ( -0.046 3.041 -0.277) 3.054 7.062 ( 1.753 -2.539 -0.768) 3.179 ======================= Grid point 51 (23/36) ======================= q-point: ( 0.33 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 76 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.509 ( 5.119 5.781 -0.970) 7.782 0.567 ( 2.203 0.944 -2.941) 3.794 0.576 ( 5.183 -1.127 -1.304) 5.462 0.592 ( 4.471 2.163 0.546) 4.996 0.652 ( 1.621 -3.686 -0.583) 4.069 0.675 ( 1.673 -4.301 0.178) 4.619 0.712 ( -1.168 0.685 2.956) 3.251 0.751 ( -0.711 -1.354 1.538) 2.169 0.779 ( -1.686 -0.293 -2.040) 2.663 0.795 ( -0.109 -2.119 -1.036) 2.361 0.870 ( -2.718 -0.959 1.104) 3.086 0.889 ( -2.032 -2.063 -1.392) 3.213 1.029 ( -0.886 0.725 -3.174) 3.374 1.051 ( -0.857 -0.707 -2.778) 2.992 1.070 ( -0.276 2.812 -3.639) 4.608 1.114 ( -0.582 -1.246 -2.730) 3.057 1.185 ( -0.156 0.701 -1.096) 1.311 1.200 ( 0.650 -1.073 -1.436) 1.907 1.264 ( 6.251 3.195 1.165) 7.117 1.268 ( -0.567 2.343 -0.105) 2.413 1.288 ( -0.662 0.026 0.086) 0.668 1.320 ( 3.690 -1.042 0.499) 3.866 1.357 ( -1.532 1.659 -2.949) 3.714 1.367 ( 0.522 1.192 -1.966) 2.358 1.399 ( 0.741 1.688 -0.459) 1.899 1.442 ( 3.869 -0.661 1.414) 4.172 1.498 ( 2.891 0.235 0.165) 2.905 1.511 ( 1.551 -0.120 1.616) 2.243 1.529 ( 1.239 -0.998 0.386) 1.637 1.549 ( 0.768 -1.881 -0.231) 2.044 1.647 ( -0.817 1.461 1.340) 2.144 1.666 ( 0.831 0.242 1.667) 1.879 1.680 ( 1.521 1.344 1.509) 2.529 1.703 ( 1.160 -1.470 0.903) 2.079 1.755 ( -0.039 1.084 2.012) 2.285 1.775 ( -0.323 -0.760 2.315) 2.458 1.803 ( 0.577 2.267 0.643) 2.426 1.832 ( 0.095 0.872 -0.022) 0.877 1.848 ( 0.264 -0.783 -0.556) 0.996 1.882 ( -1.058 -4.426 -0.202) 4.555 2.002 ( 0.301 2.226 0.128) 2.250 2.032 ( -0.592 1.516 0.978) 1.899 2.042 ( -0.314 -1.294 0.667) 1.490 2.068 ( -0.533 -2.323 0.417) 2.419 2.212 ( -0.089 2.215 -0.303) 2.237 2.242 ( 0.274 -0.512 0.114) 0.592 2.263 ( -1.023 -0.311 0.457) 1.162 2.286 ( -0.437 -2.253 1.118) 2.552 2.326 ( -1.101 -0.055 0.402) 1.174 2.331 ( -0.104 -0.168 0.656) 0.685 2.355 ( -2.830 2.551 -0.305) 3.822 2.395 ( -2.416 -1.198 0.009) 2.697 2.522 ( -1.698 2.350 -0.871) 3.027 2.577 ( -0.881 -2.486 -1.274) 2.929 2.599 ( -2.031 2.020 1.431) 3.202 2.656 ( -2.491 -2.989 0.379) 3.909 5.468 ( -0.210 0.112 0.110) 0.262 5.471 ( -0.583 -0.117 -0.300) 0.666 5.490 ( -0.456 0.113 0.296) 0.555 5.492 ( -0.478 -0.104 0.198) 0.528 6.335 ( 0.716 1.746 -1.568) 2.453 6.353 ( 1.739 0.239 -0.153) 1.762 6.408 ( 0.437 1.128 0.376) 1.267 6.433 ( -0.005 -1.237 -1.689) 2.093 6.481 ( -1.032 0.220 -0.568) 1.198 6.510 ( -1.075 -0.225 -0.170) 1.112 6.515 ( -1.567 0.695 2.119) 2.725 6.577 ( -0.903 -3.089 2.453) 4.046 6.926 ( 1.298 3.036 -1.994) 3.857 6.991 ( 2.182 -3.206 -2.101) 4.411 7.012 ( 0.236 3.813 0.411) 3.842 7.088 ( 0.605 -3.634 -0.023) 3.684 ======================= Grid point 52 (24/36) ======================= q-point: (-0.50 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 50 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.555 ( 0.000 6.403 -0.951) 6.473 0.609 ( 0.000 -0.138 -2.542) 2.545 0.621 ( -0.000 -0.370 -1.819) 1.856 0.664 ( 0.000 -2.563 -0.466) 2.605 0.687 ( 0.000 3.305 1.911) 3.818 0.712 ( 0.000 -1.801 -0.422) 1.850 0.718 ( -0.000 -0.515 0.371) 0.634 0.726 ( 0.000 -1.730 1.096) 2.048 0.750 ( -0.000 -0.030 -0.642) 0.643 0.821 ( 0.000 -1.433 2.789) 3.136 0.830 ( -0.000 -5.123 -0.241) 5.129 0.836 ( 0.000 -2.569 1.508) 2.979 1.004 ( 0.000 1.603 -5.286) 5.524 1.034 ( -0.000 -1.251 -4.089) 4.276 1.047 ( 0.000 3.260 -4.960) 5.935 1.102 ( 0.000 -2.281 -4.044) 4.643 1.206 ( -0.000 2.222 -3.141) 3.847 1.240 ( -0.000 -0.505 -3.171) 3.211 1.248 ( 0.000 2.350 0.651) 2.438 1.257 ( 0.000 0.933 -1.727) 1.963 1.286 ( 0.000 1.275 -1.367) 1.869 1.326 ( 0.000 -2.605 -3.477) 4.345 1.337 ( -0.000 0.882 3.038) 3.164 1.373 ( -0.000 -2.258 1.275) 2.593 1.503 ( 0.000 1.137 -0.815) 1.399 1.505 ( -0.000 1.604 2.222) 2.740 1.514 ( 0.000 0.049 1.322) 1.323 1.532 ( 0.000 -0.722 3.060) 3.144 1.565 ( 0.000 1.405 0.004) 1.405 1.596 ( 0.000 -0.610 0.573) 0.837 1.620 ( 0.000 3.569 0.086) 3.570 1.663 ( -0.000 3.025 -0.234) 3.034 1.679 ( 0.000 -1.626 1.916) 2.513 1.726 ( 0.000 2.918 3.828) 4.814 1.744 ( 0.000 -3.883 0.505) 3.916 1.775 ( -0.000 -1.836 2.389) 3.014 1.811 ( -0.000 2.725 0.861) 2.858 1.825 ( 0.000 1.005 0.819) 1.297 1.836 ( 0.000 -0.358 -0.295) 0.464 1.860 ( 0.000 -2.161 0.990) 2.377 2.019 ( -0.000 0.543 -0.189) 0.575 2.027 ( -0.000 -0.152 -0.680) 0.697 2.043 ( 0.000 1.544 0.806) 1.742 2.064 ( 0.000 -0.307 0.924) 0.974 2.223 ( -0.000 0.523 -0.713) 0.885 2.235 ( 0.000 0.848 0.354) 0.919 2.248 ( 0.000 -1.764 0.467) 1.825 2.274 ( 0.000 -2.841 1.388) 3.162 2.316 ( 0.000 0.939 0.014) 0.939 2.328 ( 0.000 0.771 0.520) 0.930 2.338 ( -0.000 -1.067 0.259) 1.098 2.345 ( 0.000 -0.938 0.517) 1.071 2.503 ( 0.000 3.449 -0.990) 3.588 2.566 ( 0.000 2.635 1.699) 3.135 2.580 ( -0.000 -4.220 -1.994) 4.667 2.621 ( -0.000 -2.907 0.394) 2.934 5.464 ( 0.000 0.066 -0.410) 0.415 5.468 ( 0.000 -0.264 -0.258) 0.370 5.483 ( 0.000 0.219 0.473) 0.522 5.485 ( 0.000 0.003 0.521) 0.521 6.344 ( 0.000 2.480 -1.752) 3.036 6.387 ( -0.000 -1.218 -1.200) 1.710 6.394 ( 0.000 2.537 0.622) 2.612 6.440 ( 0.000 -1.814 0.218) 1.827 6.470 ( 0.000 -0.613 -1.473) 1.595 6.486 ( -0.000 -1.823 -0.196) 1.834 6.504 ( 0.000 2.915 2.004) 3.537 6.565 ( 0.000 -2.782 2.536) 3.765 6.942 ( 0.000 3.551 -2.454) 4.316 7.013 ( 0.000 4.325 0.888) 4.416 7.020 ( 0.000 -4.200 -2.705) 4.995 7.093 ( -0.000 -3.606 0.713) 3.676 ======================= Grid point 70 (25/36) ======================= q-point: ( 0.00 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.285 ( 0.000 0.000 5.787) 5.787 0.351 ( 0.000 0.000 8.512) 8.512 0.400 ( 0.000 0.000 -6.028) 6.028 0.496 ( 0.000 0.000 -5.174) 5.174 0.549 ( 0.000 0.000 0.066) 0.066 0.619 ( 0.000 0.000 -5.103) 5.103 0.654 ( 0.000 0.000 11.844) 11.844 0.764 ( 0.000 0.000 -2.248) 2.248 0.858 ( 0.000 0.000 5.711) 5.711 0.876 ( 0.000 0.000 0.762) 0.762 0.887 ( 0.000 0.000 1.659) 1.659 0.944 ( 0.000 0.000 1.842) 1.842 0.982 ( 0.000 0.000 -6.065) 6.065 1.047 ( 0.000 0.000 -7.054) 7.054 1.049 ( 0.000 0.000 2.941) 2.941 1.123 ( 0.000 0.000 -4.906) 4.906 1.232 ( 0.000 0.000 0.144) 0.144 1.255 ( 0.000 0.000 -1.447) 1.447 1.268 ( 0.000 0.000 2.116) 2.116 1.282 ( 0.000 0.000 1.698) 1.698 1.332 ( 0.000 0.000 1.800) 1.800 1.337 ( 0.000 0.000 2.122) 2.122 1.356 ( 0.000 0.000 -0.798) 0.798 1.370 ( 0.000 0.000 -1.797) 1.797 1.403 ( 0.000 0.000 1.641) 1.641 1.406 ( 0.000 0.000 -0.121) 0.121 1.423 ( 0.000 0.000 1.136) 1.136 1.431 ( 0.000 0.000 0.518) 0.518 1.447 ( 0.000 0.000 -2.512) 2.512 1.480 ( 0.000 0.000 -5.218) 5.218 1.601 ( 0.000 0.000 2.646) 2.646 1.611 ( 0.000 0.000 1.974) 1.974 1.633 ( 0.000 0.000 -0.126) 0.126 1.635 ( 0.000 0.000 -1.738) 1.738 1.715 ( 0.000 0.000 1.249) 1.249 1.766 ( 0.000 0.000 -3.030) 3.030 1.783 ( 0.000 0.000 1.347) 1.347 1.806 ( 0.000 0.000 -1.169) 1.169 1.893 ( 0.000 0.000 1.571) 1.571 1.914 ( 0.000 0.000 -0.598) 0.598 2.010 ( 0.000 0.000 0.912) 0.912 2.014 ( 0.000 0.000 0.662) 0.662 2.126 ( 0.000 0.000 -0.265) 0.265 2.128 ( 0.000 0.000 -0.333) 0.333 2.183 ( 0.000 0.000 3.896) 3.896 2.254 ( 0.000 0.000 -2.872) 2.872 2.261 ( 0.000 0.000 2.041) 2.041 2.280 ( 0.000 0.000 0.296) 0.296 2.323 ( 0.000 0.000 0.497) 0.497 2.326 ( 0.000 0.000 -0.328) 0.328 2.336 ( 0.000 0.000 -1.482) 1.482 2.357 ( 0.000 0.000 -3.094) 3.094 2.609 ( 0.000 0.000 4.081) 4.081 2.661 ( 0.000 0.000 2.985) 2.985 2.681 ( 0.000 0.000 -3.166) 3.166 2.711 ( 0.000 0.000 -2.093) 2.093 5.483 ( 0.000 0.000 0.850) 0.850 5.499 ( 0.000 0.000 -0.634) 0.634 5.502 ( 0.000 0.000 0.355) 0.355 5.509 ( 0.000 0.000 -0.333) 0.333 6.261 ( 0.000 0.000 1.198) 1.198 6.291 ( 0.000 0.000 -1.783) 1.783 6.367 ( 0.000 0.000 2.682) 2.682 6.414 ( 0.000 0.000 -1.933) 1.933 6.486 ( 0.000 0.000 0.813) 0.813 6.506 ( 0.000 0.000 -1.257) 1.257 6.732 ( 0.000 0.000 2.069) 2.069 6.770 ( 0.000 0.000 -1.705) 1.705 6.862 ( 0.000 0.000 0.990) 0.990 6.876 ( 0.000 0.000 -0.258) 0.258 6.958 ( 0.000 0.000 -0.196) 0.196 6.967 ( 0.000 0.000 -0.921) 0.921 ======================= Grid point 71 (26/36) ======================= q-point: ( 0.17 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 45 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.376 ( 6.747 0.000 3.581) 7.639 0.438 ( 4.205 0.000 6.578) 7.807 0.455 ( 4.186 0.000 -4.271) 5.980 0.533 ( 2.088 0.000 -1.831) 2.777 0.611 ( 3.919 0.000 -1.955) 4.380 0.651 ( 2.249 0.000 -3.116) 3.843 0.752 ( 9.674 0.000 6.640) 11.734 0.801 ( 3.744 0.000 -0.862) 3.842 0.846 ( -1.097 0.000 3.265) 3.444 0.864 ( -0.680 0.000 0.557) 0.879 0.883 ( -0.849 0.000 0.637) 1.062 0.957 ( 0.114 0.000 1.820) 1.824 0.974 ( -0.564 0.000 -6.486) 6.510 1.018 ( -2.529 0.000 -4.330) 5.015 1.050 ( 0.286 0.000 3.226) 3.239 1.118 ( -0.966 0.000 -3.688) 3.813 1.136 ( -2.162 0.000 -0.738) 2.284 1.139 ( -4.066 0.000 -0.197) 4.070 1.224 ( -4.730 0.000 2.909) 5.553 1.276 ( -4.135 0.000 -2.091) 4.634 1.305 ( -3.601 0.000 0.397) 3.622 1.320 ( -3.453 0.000 -1.104) 3.626 1.355 ( -1.558 0.000 1.700) 2.306 1.370 ( 1.344 0.000 2.336) 2.695 1.398 ( -0.595 0.000 -2.520) 2.590 1.415 ( -0.283 0.000 0.714) 0.767 1.445 ( 2.821 0.000 1.903) 3.403 1.465 ( 4.141 0.000 1.374) 4.362 1.490 ( 4.708 0.000 0.100) 4.709 1.550 ( 2.821 0.000 -5.051) 5.786 1.567 ( 1.864 0.000 1.441) 2.356 1.586 ( -0.310 0.000 -0.142) 0.341 1.665 ( 2.220 0.000 -0.161) 2.226 1.665 ( 1.957 0.000 0.279) 1.977 1.738 ( 1.881 0.000 1.740) 2.563 1.778 ( 1.781 0.000 -1.990) 2.670 1.862 ( -2.361 0.000 1.347) 2.718 1.877 ( -0.015 0.000 -0.136) 0.136 1.903 ( 1.992 0.000 0.652) 2.096 1.925 ( 5.584 0.000 -1.872) 5.890 1.984 ( -2.728 0.000 0.847) 2.856 1.992 ( -2.056 0.000 -0.126) 2.060 2.136 ( -3.120 0.000 2.583) 4.051 2.157 ( 1.825 0.000 0.644) 1.935 2.197 ( 4.977 0.000 -0.889) 5.056 2.206 ( -2.686 0.000 -1.687) 3.172 2.233 ( -0.310 0.000 1.136) 1.177 2.253 ( -0.084 0.000 -1.345) 1.348 2.298 ( -1.020 0.000 2.087) 2.323 2.348 ( -0.431 0.000 -2.130) 2.174 2.397 ( 3.267 0.000 -0.556) 3.314 2.398 ( 2.650 0.000 -0.979) 2.825 2.570 ( -3.529 0.000 3.230) 4.784 2.597 ( -4.589 0.000 0.867) 4.670 2.679 ( -1.047 0.000 -1.144) 1.551 2.692 ( -1.414 0.000 -2.135) 2.560 5.481 ( -0.205 0.000 0.875) 0.899 5.494 ( -0.658 0.000 0.241) 0.701 5.496 ( -0.325 0.000 -0.570) 0.655 5.499 ( -0.713 0.000 -0.323) 0.783 6.272 ( 0.958 0.000 1.327) 1.637 6.306 ( 1.087 0.000 -1.984) 2.262 6.384 ( 1.342 0.000 2.676) 2.994 6.427 ( 0.861 0.000 -1.425) 1.665 6.474 ( -1.009 0.000 0.824) 1.303 6.498 ( -0.712 0.000 -1.488) 1.649 6.677 ( -3.822 0.000 0.458) 3.849 6.682 ( -4.583 0.000 -0.222) 4.589 6.877 ( 0.148 0.000 0.434) 0.458 6.880 ( 0.096 0.000 0.092) 0.133 7.011 ( 3.941 0.000 0.730) 4.008 7.046 ( 4.875 0.000 -3.008) 5.729 ======================= Grid point 72 (27/36) ======================= q-point: ( 0.33 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-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.488 ( 1.593 0.000 5.285) 5.520 0.501 ( 3.748 0.000 1.697) 4.114 0.555 ( 4.469 0.000 -2.153) 4.960 0.586 ( 2.962 0.000 -2.126) 3.646 0.653 ( 0.558 0.000 -0.348) 0.658 0.686 ( 0.530 0.000 -1.955) 2.026 0.796 ( -3.360 0.000 -3.076) 4.555 0.811 ( -3.392 0.000 -3.501) 4.875 0.862 ( 1.139 0.000 0.764) 1.372 0.868 ( 0.118 0.000 0.185) 0.220 0.912 ( 0.609 0.000 2.179) 2.262 0.942 ( -3.432 0.000 -0.070) 3.432 0.967 ( 1.383 0.000 -0.397) 1.438 0.986 ( 0.441 0.000 -3.011) 3.043 1.068 ( 1.176 0.000 3.659) 3.844 1.124 ( 0.978 0.000 -2.009) 2.235 1.159 ( 6.997 0.000 2.400) 7.397 1.169 ( -0.882 0.000 1.745) 1.956 1.189 ( 4.055 0.000 -1.284) 4.253 1.200 ( -3.243 0.000 -0.122) 3.245 1.234 ( 0.252 0.000 -1.058) 1.088 1.237 ( -2.112 0.000 -0.322) 2.136 1.313 ( -1.162 0.000 1.058) 1.572 1.369 ( -0.402 0.000 1.722) 1.769 1.372 ( -0.734 0.000 -4.478) 4.538 1.411 ( -0.198 0.000 -1.962) 1.972 1.497 ( 2.693 0.000 3.742) 4.611 1.560 ( 1.208 0.000 -1.138) 1.659 1.586 ( 1.612 0.000 0.106) 1.616 1.588 ( 2.054 0.000 1.920) 2.812 1.596 ( 0.386 0.000 -1.786) 1.827 1.616 ( 1.279 0.000 -0.825) 1.522 1.716 ( 2.042 0.000 0.854) 2.213 1.719 ( 2.737 0.000 0.274) 2.751 1.763 ( -0.384 0.000 3.267) 3.289 1.806 ( -3.743 0.000 -0.603) 3.791 1.827 ( 1.322 0.000 0.031) 1.322 1.830 ( 1.473 0.000 0.107) 1.477 1.890 ( -3.151 0.000 0.995) 3.304 1.910 ( -4.192 0.000 -1.049) 4.321 1.973 ( 2.027 0.000 0.956) 2.241 1.990 ( 1.228 0.000 -0.619) 1.375 2.071 ( -2.944 0.000 1.463) 3.287 2.104 ( -4.219 0.000 -1.784) 4.581 2.226 ( -0.871 0.000 1.568) 1.793 2.241 ( 2.823 0.000 1.585) 3.238 2.260 ( -0.595 0.000 -1.647) 1.751 2.269 ( 1.777 0.000 -1.025) 2.051 2.331 ( 2.903 0.000 1.059) 3.090 2.344 ( 0.078 0.000 -0.065) 0.102 2.420 ( -1.496 0.000 0.848) 1.720 2.429 ( 0.295 0.000 -0.082) 0.307 2.501 ( -1.456 0.000 1.710) 2.246 2.506 ( -2.727 0.000 0.650) 2.804 2.646 ( -1.560 0.000 -1.341) 2.058 2.659 ( -1.545 0.000 -2.976) 3.353 5.476 ( -0.139 0.000 0.563) 0.579 5.480 ( -0.621 0.000 -0.003) 0.621 5.481 ( -0.555 0.000 0.120) 0.567 5.490 ( -0.480 0.000 -0.460) 0.665 6.297 ( 1.220 0.000 1.315) 1.794 6.324 ( 0.410 0.000 -1.120) 1.192 6.411 ( 0.966 0.000 2.219) 2.420 6.439 ( 0.065 0.000 -0.097) 0.117 6.451 ( -0.699 0.000 0.029) 0.700 6.480 ( -0.925 0.000 -2.189) 2.377 6.608 ( -2.215 0.000 -0.075) 2.216 6.613 ( -2.036 0.000 -0.742) 2.167 6.877 ( 0.007 0.000 0.318) 0.318 6.880 ( -0.056 0.000 0.002) 0.056 7.091 ( 3.125 0.000 0.743) 3.213 7.121 ( 2.050 0.000 -2.411) 3.164 ======================= Grid point 73 (28/36) ======================= q-point: (-0.50 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.520 ( -0.154 0.000 3.109) 3.113 0.520 ( 0.154 0.000 3.109) 3.113 0.627 ( -1.320 0.000 -1.088) 1.710 0.627 ( 1.320 0.000 -1.088) 1.710 0.662 ( -1.643 0.000 -2.691) 3.153 0.662 ( 1.643 0.000 -2.691) 3.153 0.757 ( -0.521 0.000 -3.466) 3.505 0.757 ( 0.521 0.000 -3.466) 3.505 0.888 ( -0.004 0.000 -0.941) 0.941 0.888 ( 0.004 0.000 -0.941) 0.941 0.893 ( -0.004 0.000 0.667) 0.667 0.893 ( 0.004 0.000 0.667) 0.667 0.980 ( -0.190 0.000 1.363) 1.376 0.980 ( 0.190 0.000 1.363) 1.376 1.112 ( -1.308 0.000 -0.139) 1.315 1.112 ( 1.308 0.000 -0.139) 1.315 1.158 ( -0.001 0.000 0.032) 0.032 1.158 ( 0.001 0.000 0.032) 0.032 1.171 ( -0.002 0.000 0.747) 0.747 1.171 ( 0.002 0.000 0.747) 0.747 1.303 ( -0.000 0.000 0.493) 0.493 1.303 ( 0.000 0.000 0.493) 0.493 1.346 ( -0.767 0.000 -2.697) 2.804 1.346 ( 0.767 0.000 -2.697) 2.804 1.456 ( -0.382 0.000 -0.167) 0.417 1.456 ( 0.382 0.000 -0.167) 0.417 1.502 ( -0.011 0.000 -0.186) 0.186 1.502 ( 0.011 0.000 -0.186) 0.186 1.549 ( -0.911 0.000 0.395) 0.993 1.549 ( 0.911 0.000 0.395) 0.993 1.632 ( -0.023 0.000 0.503) 0.504 1.632 ( 0.023 0.000 0.503) 0.504 1.744 ( -0.610 0.000 2.144) 2.229 1.744 ( 0.610 0.000 2.144) 2.229 1.770 ( -0.996 0.000 0.244) 1.026 1.770 ( 0.996 0.000 0.244) 1.026 1.822 ( -0.133 0.000 0.795) 0.806 1.822 ( 0.133 0.000 0.795) 0.806 1.901 ( -0.172 0.000 0.340) 0.381 1.901 ( 0.172 0.000 0.340) 0.381 1.999 ( -0.746 0.000 -0.768) 1.070 1.999 ( 0.746 0.000 -0.768) 1.070 2.029 ( -0.007 0.000 0.480) 0.480 2.029 ( 0.007 0.000 0.480) 0.480 2.223 ( -0.773 0.000 -0.741) 1.070 2.223 ( 0.773 0.000 -0.741) 1.070 2.310 ( -2.739 0.000 0.665) 2.818 2.310 ( 2.739 0.000 0.665) 2.818 2.317 ( -1.412 0.000 0.161) 1.421 2.317 ( 1.412 0.000 0.161) 1.421 2.421 ( -0.237 0.000 0.856) 0.889 2.421 ( 0.237 0.000 0.856) 0.889 2.490 ( -0.010 0.000 1.702) 1.702 2.490 ( 0.010 0.000 1.702) 1.702 2.633 ( -0.105 0.000 -2.738) 2.740 2.633 ( 0.105 0.000 -2.738) 2.740 5.474 ( -0.002 0.000 0.123) 0.123 5.474 ( 0.002 0.000 0.123) 0.123 5.479 ( -0.218 0.000 -0.008) 0.219 5.479 ( 0.218 0.000 -0.008) 0.219 6.319 ( -0.153 0.000 0.579) 0.599 6.319 ( 0.153 0.000 0.579) 0.599 6.427 ( -0.188 0.000 1.431) 1.443 6.427 ( 0.188 0.000 1.431) 1.443 6.459 ( -0.625 0.000 -2.174) 2.262 6.459 ( 0.625 0.000 -2.174) 2.262 6.588 ( -0.013 0.000 -0.243) 0.243 6.588 ( 0.013 0.000 -0.243) 0.243 6.879 ( -0.002 0.000 0.072) 0.072 6.879 ( 0.002 0.000 0.072) 0.072 7.134 ( -0.148 0.000 -0.604) 0.622 7.134 ( 0.148 0.000 -0.604) 0.622 ======================= Grid point 77 (29/36) ======================= q-point: ( 0.00 0.20 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 52 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.324 ( 0.000 3.427 6.047) 6.950 0.390 ( 0.000 2.096 3.392) 3.987 0.443 ( 0.000 4.997 -1.797) 5.310 0.465 ( 0.000 -2.902 2.155) 3.615 0.549 ( 0.000 3.169 -5.930) 6.724 0.589 ( 0.000 -2.313 -7.526) 7.874 0.769 ( 0.000 7.900 6.674) 10.341 0.776 ( 0.000 0.262 3.335) 3.345 0.837 ( 0.000 -1.408 1.466) 2.033 0.837 ( 0.000 -2.740 2.945) 4.022 0.865 ( 0.000 -2.808 -1.451) 3.161 0.896 ( 0.000 -3.660 -0.456) 3.688 1.001 ( 0.000 -3.933 -3.806) 5.473 1.016 ( 0.000 2.597 -3.613) 4.449 1.045 ( 0.000 -0.355 1.977) 2.009 1.074 ( 0.000 -4.316 -2.267) 4.875 1.197 ( 0.000 -5.653 1.240) 5.787 1.217 ( 0.000 -1.676 -0.832) 1.871 1.245 ( 0.000 -1.710 -1.092) 2.029 1.251 ( 0.000 -2.528 -1.214) 2.805 1.313 ( 0.000 1.730 3.198) 3.636 1.330 ( 0.000 -2.231 5.062) 5.532 1.369 ( 0.000 0.785 -2.536) 2.655 1.394 ( 0.000 -0.828 0.510) 0.973 1.410 ( 0.000 0.379 -0.593) 0.704 1.423 ( 0.000 1.477 0.314) 1.510 1.441 ( 0.000 -1.663 -2.069) 2.654 1.445 ( 0.000 0.273 -0.744) 0.792 1.462 ( 0.000 3.223 -2.200) 3.902 1.464 ( 0.000 3.382 -0.858) 3.489 1.575 ( 0.000 0.244 4.303) 4.310 1.616 ( 0.000 -1.339 0.231) 1.358 1.652 ( 0.000 0.866 -0.857) 1.218 1.682 ( 0.000 2.571 -3.457) 4.308 1.741 ( 0.000 2.172 0.493) 2.228 1.752 ( 0.000 0.384 -0.436) 0.581 1.797 ( 0.000 0.304 1.836) 1.861 1.828 ( 0.000 1.323 -0.815) 1.554 1.859 ( 0.000 -3.445 0.252) 3.454 1.865 ( 0.000 -2.521 -0.604) 2.593 2.014 ( 0.000 0.905 0.990) 1.341 2.037 ( 0.000 1.348 -1.274) 1.854 2.112 ( 0.000 -2.147 0.947) 2.346 2.128 ( 0.000 -1.591 -0.743) 1.756 2.207 ( 0.000 -0.234 1.464) 1.482 2.207 ( 0.000 2.871 1.564) 3.269 2.281 ( 0.000 -0.607 1.081) 1.240 2.297 ( 0.000 -0.124 -0.665) 0.676 2.334 ( 0.000 1.136 -0.408) 1.208 2.343 ( 0.000 1.090 0.007) 1.090 2.376 ( 0.000 3.935 -1.891) 4.365 2.382 ( 0.000 2.113 -1.712) 2.719 2.620 ( 0.000 0.671 1.910) 2.024 2.648 ( 0.000 -1.742 -1.062) 2.040 2.665 ( 0.000 0.727 4.193) 4.255 2.737 ( 0.000 0.957 -2.920) 3.073 5.485 ( 0.000 0.125 0.721) 0.732 5.497 ( 0.000 -0.103 -0.563) 0.573 5.502 ( 0.000 -0.057 0.265) 0.271 5.506 ( 0.000 -0.231 -0.214) 0.315 6.275 ( 0.000 1.205 0.334) 1.250 6.294 ( 0.000 0.278 -1.383) 1.410 6.365 ( 0.000 -0.115 2.239) 2.242 6.405 ( 0.000 -0.810 -1.745) 1.924 6.497 ( 0.000 1.097 0.585) 1.243 6.513 ( 0.000 0.656 -0.972) 1.172 6.691 ( 0.000 -3.358 1.100) 3.534 6.695 ( 0.000 -4.834 1.096) 4.957 6.871 ( 0.000 0.763 1.000) 1.258 6.898 ( 0.000 0.667 -1.858) 1.974 6.986 ( 0.000 1.428 -0.454) 1.498 6.986 ( 0.000 1.944 -0.537) 2.017 ======================= Grid point 78 (30/36) ======================= q-point: ( 0.17 0.20 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-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.388 ( 5.458 1.309 4.372) 7.115 0.431 ( 2.828 -0.287 2.013) 3.483 0.468 ( 2.811 0.956 -2.448) 3.849 0.529 ( 3.919 0.305 2.158) 4.484 0.603 ( 3.277 -0.331 -4.281) 5.402 0.646 ( 4.183 -0.972 -4.686) 6.356 0.754 ( 0.735 1.010 2.064) 2.412 0.808 ( -1.164 -3.107 0.863) 3.428 0.829 ( -0.552 -1.520 3.867) 4.192 0.843 ( -1.677 -1.614 -1.676) 2.868 0.868 ( 1.430 -2.048 1.902) 3.139 0.937 ( 5.756 2.809 -0.298) 6.412 0.975 ( -0.454 -0.609 -5.350) 5.404 1.001 ( -1.549 -2.051 1.011) 2.762 1.058 ( -0.434 -1.288 0.788) 1.571 1.075 ( 2.807 1.517 -0.901) 3.316 1.154 ( -4.003 -1.759 1.453) 4.608 1.181 ( -1.481 -0.204 0.950) 1.772 1.196 ( -2.908 2.041 -2.022) 4.088 1.241 ( -2.573 -1.540 -2.792) 4.097 1.315 ( -0.417 -0.768 1.607) 1.829 1.325 ( -1.462 -1.825 0.444) 2.380 1.344 ( -1.779 0.713 -0.580) 2.002 1.348 ( -1.199 -1.539 1.012) 2.197 1.420 ( -1.831 0.739 -0.888) 2.165 1.431 ( 1.842 0.075 2.242) 2.902 1.434 ( 1.247 0.142 -1.043) 1.632 1.487 ( 3.120 1.456 0.855) 3.548 1.502 ( 1.653 1.158 -0.998) 2.252 1.521 ( 3.771 0.244 -0.355) 3.795 1.573 ( 1.360 0.462 -1.622) 2.166 1.611 ( 0.852 1.338 0.147) 1.593 1.678 ( 1.302 0.924 1.248) 2.026 1.711 ( 2.205 2.843 -0.340) 3.614 1.740 ( 0.475 -0.064 0.475) 0.675 1.776 ( 1.309 0.002 -0.475) 1.392 1.782 ( 0.895 -1.890 -0.017) 2.091 1.853 ( -0.228 -2.627 0.236) 2.647 1.884 ( 1.594 -2.709 0.455) 3.177 1.935 ( 3.274 -2.730 -0.516) 4.294 1.987 ( -1.563 1.723 1.095) 2.571 2.011 ( -1.699 2.031 -1.186) 2.901 2.097 ( -0.649 -2.298 0.103) 2.390 2.171 ( -1.771 -2.418 0.493) 3.037 2.201 ( -0.488 1.641 1.468) 2.256 2.226 ( 0.377 2.347 0.408) 2.412 2.248 ( 2.724 1.531 0.782) 3.221 2.274 ( 0.265 1.076 -0.591) 1.256 2.332 ( -2.195 1.826 -1.237) 3.112 2.348 ( 1.097 -0.174 0.992) 1.489 2.374 ( -0.729 0.326 -1.874) 2.037 2.409 ( 1.744 0.110 -0.506) 1.819 2.570 ( -4.258 0.326 0.835) 4.351 2.619 ( -2.003 0.712 1.233) 2.457 2.651 ( -1.802 -0.549 1.068) 2.166 2.715 ( -1.796 0.705 -2.212) 2.936 5.482 ( -0.259 0.054 0.691) 0.739 5.489 ( -0.889 -0.332 -0.047) 0.950 5.496 ( -0.265 0.028 -0.453) 0.526 5.502 ( -0.405 0.091 0.007) 0.415 6.279 ( 0.438 0.702 0.152) 0.842 6.314 ( 1.490 0.562 -0.982) 1.871 6.384 ( 1.451 0.004 1.960) 2.438 6.420 ( 0.934 -0.765 -1.351) 1.812 6.486 ( -0.884 1.128 0.565) 1.540 6.503 ( -0.841 0.656 -0.984) 1.451 6.625 ( -3.689 -2.962 0.417) 4.749 6.682 ( -2.143 -1.530 0.694) 2.724 6.865 ( -0.251 0.324 0.141) 0.433 6.935 ( 0.878 2.912 -0.289) 3.055 6.965 ( 2.331 -2.356 -0.727) 3.393 7.066 ( 3.678 0.337 -1.092) 3.851 ======================= Grid point 79 (31/36) ======================= q-point: ( 0.33 0.20 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 76 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.461 ( 1.528 0.322 2.706) 3.124 0.525 ( 4.595 -0.154 -0.128) 4.599 0.551 ( 4.263 -0.074 -1.902) 4.669 0.608 ( 3.473 1.530 1.606) 4.121 0.661 ( 1.522 0.886 0.451) 1.818 0.692 ( 0.256 0.206 -2.337) 2.360 0.763 ( -1.592 -3.931 -2.448) 4.897 0.768 ( -0.769 -3.479 -2.358) 4.273 0.800 ( -0.748 -3.387 1.001) 3.610 0.823 ( -0.570 -1.669 -1.322) 2.204 0.874 ( 0.319 -1.794 1.276) 2.225 0.941 ( -2.170 -1.206 -3.899) 4.623 0.977 ( -0.492 -0.188 -0.252) 0.584 1.035 ( -2.288 -1.826 3.448) 4.523 1.046 ( -0.927 0.293 0.894) 1.321 1.100 ( 0.800 -0.466 -1.462) 1.731 1.162 ( -0.132 0.137 -0.203) 0.278 1.175 ( -0.615 0.938 -0.898) 1.437 1.205 ( 5.416 1.954 -0.029) 5.757 1.232 ( 1.473 1.959 0.606) 2.525 1.250 ( -1.695 1.354 -0.995) 2.386 1.282 ( -1.013 1.961 -0.324) 2.231 1.317 ( 3.928 1.838 -0.093) 4.337 1.353 ( -1.565 -0.001 -2.164) 2.671 1.382 ( 0.322 0.657 2.663) 2.762 1.447 ( 0.148 1.641 -3.017) 3.437 1.477 ( 2.184 -1.307 2.015) 3.246 1.536 ( 2.142 -1.317 1.329) 2.844 1.558 ( 0.656 -0.446 -0.842) 1.157 1.570 ( 0.817 -1.567 -1.317) 2.203 1.645 ( 3.032 2.334 1.142) 3.993 1.654 ( 2.636 2.864 0.327) 3.906 1.690 ( -0.231 -0.993 0.523) 1.146 1.750 ( 1.132 0.432 1.006) 1.575 1.757 ( 1.321 0.780 1.453) 2.113 1.803 ( 1.345 -1.086 0.474) 1.793 1.824 ( 0.272 0.260 0.054) 0.380 1.826 ( -0.777 -0.199 0.243) 0.838 1.859 ( -1.944 -1.740 -0.038) 2.609 1.922 ( -2.683 -0.291 -0.413) 2.730 1.984 ( 0.403 0.938 0.702) 1.239 2.007 ( 0.474 1.623 0.005) 1.691 2.070 ( -1.677 -0.361 0.708) 1.856 2.095 ( -3.007 -0.736 -0.869) 3.215 2.205 ( 0.310 -0.935 0.863) 1.310 2.252 ( 1.098 0.393 0.823) 1.427 2.280 ( -0.825 1.710 -1.235) 2.265 2.296 ( 0.152 2.227 -0.113) 2.235 2.309 ( 1.162 -1.819 1.117) 2.430 2.337 ( -0.139 -0.707 -0.619) 0.950 2.395 ( -0.730 -1.849 1.470) 2.472 2.416 ( -0.177 -2.449 0.171) 2.461 2.480 ( -2.384 2.042 0.369) 3.161 2.572 ( -2.015 2.513 1.215) 3.443 2.606 ( -1.755 -2.371 -1.820) 3.467 2.672 ( -1.998 0.071 -1.440) 2.464 5.470 ( -0.698 -0.418 -0.102) 0.820 5.475 ( -0.264 -0.171 0.382) 0.495 5.488 ( -0.379 0.253 -0.040) 0.457 5.493 ( -0.415 0.292 -0.035) 0.509 6.291 ( 0.597 0.228 -0.026) 0.640 6.347 ( 1.406 1.207 0.015) 1.853 6.412 ( 0.659 0.034 1.230) 1.396 6.429 ( 0.124 -0.867 -0.663) 1.098 6.468 ( -0.765 0.652 -0.169) 1.019 6.484 ( -0.726 1.340 -0.500) 1.604 6.558 ( -2.356 -2.496 -0.290) 3.444 6.638 ( -1.589 0.177 0.252) 1.619 6.863 ( 0.112 0.448 -0.378) 0.596 6.938 ( 0.090 3.739 0.201) 3.745 7.035 ( 2.847 -4.135 -0.839) 5.090 7.125 ( 1.807 -0.790 -0.723) 2.101 ======================= Grid point 80 (32/36) ======================= q-point: (-0.50 0.20 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 50 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.479 ( 0.000 1.004 1.733) 2.003 0.587 ( 0.000 -0.427 0.165) 0.458 0.610 ( -0.000 -0.357 -0.728) 0.811 0.641 ( 0.000 1.025 -0.401) 1.100 0.655 ( -0.000 2.150 1.508) 2.626 0.679 ( 0.000 0.680 -3.524) 3.589 0.751 ( -0.000 -0.856 -0.687) 1.098 0.768 ( 0.000 0.521 -4.494) 4.524 0.803 ( 0.000 -4.324 0.390) 4.342 0.821 ( 0.000 -3.481 -1.165) 3.671 0.881 ( -0.000 -2.293 0.152) 2.298 0.902 ( 0.000 -0.917 -2.149) 2.337 0.976 ( 0.000 0.126 0.433) 0.451 1.004 ( 0.000 -2.681 5.477) 6.099 1.054 ( -0.000 -2.698 1.281) 2.987 1.092 ( 0.000 -3.127 -1.650) 3.536 1.107 ( 0.000 -0.555 -4.177) 4.214 1.168 ( 0.000 0.940 -2.051) 2.256 1.201 ( 0.000 2.405 -0.574) 2.473 1.228 ( 0.000 2.061 1.390) 2.485 1.275 ( 0.000 -0.693 0.080) 0.698 1.309 ( 0.000 -1.233 -2.286) 2.598 1.359 ( 0.000 2.960 1.955) 3.547 1.409 ( -0.000 2.843 -0.475) 2.882 1.466 ( -0.000 -1.064 0.432) 1.148 1.468 ( 0.000 2.430 -0.112) 2.433 1.496 ( -0.000 0.264 -2.191) 2.207 1.527 ( 0.000 1.567 0.235) 1.584 1.545 ( 0.000 2.096 0.812) 2.248 1.578 ( 0.000 1.327 -0.132) 1.334 1.601 ( -0.000 0.070 0.903) 0.906 1.686 ( -0.000 -3.091 0.492) 3.130 1.696 ( 0.000 3.800 0.709) 3.865 1.770 ( -0.000 -3.003 1.542) 3.376 1.786 ( 0.000 2.372 1.419) 2.764 1.803 ( 0.000 2.061 0.920) 2.257 1.810 ( -0.000 1.657 0.741) 1.815 1.828 ( 0.000 -0.339 0.707) 0.784 1.851 ( 0.000 -3.040 -0.049) 3.040 1.911 ( 0.000 -0.569 0.640) 0.856 2.000 ( 0.000 0.450 -0.899) 1.006 2.017 ( -0.000 1.513 -0.123) 1.518 2.018 ( 0.000 0.023 0.358) 0.359 2.060 ( 0.000 1.561 0.234) 1.578 2.208 ( 0.000 -0.352 -0.483) 0.598 2.236 ( 0.000 0.317 0.052) 0.321 2.298 ( 0.000 -0.303 -0.495) 0.580 2.309 ( 0.000 -1.138 1.215) 1.664 2.321 ( 0.000 -0.307 -0.010) 0.307 2.333 ( -0.000 1.372 0.527) 1.469 2.368 ( 0.000 -1.375 -0.042) 1.375 2.389 ( 0.000 -4.589 0.652) 4.635 2.481 ( 0.000 1.205 2.706) 2.963 2.549 ( 0.000 2.508 1.114) 2.744 2.589 ( 0.000 -2.786 -3.611) 4.561 2.645 ( 0.000 -0.272 -1.130) 1.162 5.462 ( 0.000 -0.442 -0.076) 0.449 5.473 ( 0.000 -0.387 -0.013) 0.387 5.483 ( 0.000 0.367 0.232) 0.434 5.487 ( 0.000 0.451 0.077) 0.457 6.298 ( 0.000 0.018 -0.119) 0.120 6.365 ( 0.000 2.321 0.922) 2.498 6.407 ( 0.000 -1.235 0.841) 1.494 6.447 ( -0.000 -0.437 -1.812) 1.864 6.448 ( 0.000 0.061 -0.718) 0.720 6.477 ( 0.000 2.354 0.540) 2.416 6.532 ( 0.000 -2.898 -0.236) 2.908 6.620 ( 0.000 0.494 0.349) 0.605 6.866 ( 0.000 0.683 -0.586) 0.899 6.939 ( 0.000 3.913 0.243) 3.920 7.067 ( 0.000 -4.414 -1.053) 4.538 7.144 ( 0.000 -1.040 -0.183) 1.056 ======================= Grid point 84 (33/36) ======================= q-point: ( 0.00 0.40 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 52 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.390 ( 0.000 2.193 6.236) 6.610 0.408 ( 0.000 -0.025 1.586) 1.586 0.456 ( 0.000 2.544 -3.867) 4.629 0.526 ( 0.000 -2.624 -4.983) 5.632 0.582 ( 0.000 6.201 -2.908) 6.849 0.631 ( 0.000 3.909 -5.578) 6.811 0.746 ( 0.000 -4.916 2.259) 5.410 0.747 ( 0.000 -2.743 2.312) 3.587 0.816 ( 0.000 -1.622 -0.580) 1.722 0.817 ( 0.000 -1.397 -0.197) 1.411 0.859 ( 0.000 1.244 0.276) 1.274 0.901 ( 0.000 -1.242 -3.163) 3.398 0.953 ( 0.000 4.291 4.119) 5.948 0.985 ( 0.000 -3.020 0.524) 3.065 1.005 ( 0.000 -2.060 0.475) 2.114 1.062 ( 0.000 -0.178 -3.757) 3.761 1.142 ( 0.000 -0.203 4.305) 4.310 1.152 ( 0.000 -3.180 2.528) 4.062 1.213 ( 0.000 -1.193 -3.069) 3.293 1.219 ( 0.000 -1.158 -2.428) 2.690 1.285 ( 0.000 -0.889 2.156) 2.332 1.303 ( 0.000 -1.953 1.702) 2.590 1.381 ( 0.000 0.135 -0.392) 0.415 1.382 ( 0.000 -0.073 -0.286) 0.295 1.416 ( 0.000 0.266 -0.217) 0.344 1.419 ( 0.000 -0.183 -0.766) 0.787 1.468 ( 0.000 1.603 0.109) 1.606 1.472 ( 0.000 1.708 -0.361) 1.745 1.523 ( 0.000 2.778 0.213) 2.786 1.539 ( 0.000 2.803 -0.808) 2.917 1.592 ( 0.000 -0.335 -0.333) 0.473 1.595 ( 0.000 1.917 1.000) 2.162 1.637 ( 0.000 -2.135 -0.594) 2.216 1.666 ( 0.000 -3.411 -1.448) 3.706 1.765 ( 0.000 -0.613 0.213) 0.649 1.769 ( 0.000 -0.189 -0.032) 0.192 1.803 ( 0.000 0.775 0.779) 1.099 1.810 ( 0.000 -0.830 0.389) 0.917 1.843 ( 0.000 1.749 -0.016) 1.749 1.861 ( 0.000 0.523 -1.704) 1.783 2.038 ( 0.000 1.074 0.460) 1.169 2.051 ( 0.000 -0.012 -0.766) 0.767 2.076 ( 0.000 -0.781 0.755) 1.087 2.088 ( 0.000 -1.627 -0.452) 1.688 2.223 ( 0.000 1.294 1.781) 2.201 2.249 ( 0.000 0.053 0.061) 0.081 2.259 ( 0.000 -0.237 0.223) 0.325 2.279 ( 0.000 -1.204 -1.034) 1.587 2.358 ( 0.000 0.967 0.144) 0.977 2.363 ( 0.000 0.559 0.073) 0.564 2.426 ( 0.000 1.755 -0.508) 1.827 2.432 ( 0.000 1.332 -1.202) 1.794 2.623 ( 0.000 -0.157 0.859) 0.873 2.625 ( 0.000 -0.562 0.383) 0.680 2.697 ( 0.000 2.148 2.593) 3.367 2.733 ( 0.000 -1.233 -0.849) 1.497 5.488 ( 0.000 0.271 0.351) 0.444 5.494 ( 0.000 -0.252 -0.200) 0.321 5.501 ( 0.000 -0.006 0.073) 0.073 5.502 ( 0.000 -0.151 -0.059) 0.162 6.296 ( 0.000 0.643 -0.376) 0.745 6.299 ( 0.000 0.176 -0.651) 0.674 6.368 ( 0.000 0.537 0.715) 0.894 6.386 ( 0.000 -1.081 -0.985) 1.462 6.524 ( 0.000 1.271 -0.026) 1.271 6.529 ( 0.000 0.710 -0.436) 0.833 6.623 ( 0.000 -2.051 1.592) 2.596 6.629 ( 0.000 -2.416 0.834) 2.556 6.888 ( 0.000 0.843 0.241) 0.877 6.900 ( 0.000 -0.334 -1.004) 1.058 7.008 ( 0.000 0.654 -0.308) 0.723 7.010 ( 0.000 0.460 -0.588) 0.747 ======================= Grid point 85 (34/36) ======================= q-point: ( 0.17 0.40 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-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.424 ( 2.805 2.001 3.269) 4.749 0.459 ( 4.845 0.692 1.931) 5.262 0.465 ( 1.605 1.703 -2.735) 3.600 0.553 ( 2.527 -3.716 -1.299) 4.678 0.575 ( -0.056 2.692 -2.957) 3.999 0.628 ( -0.091 -1.562 -3.843) 4.149 0.737 ( 0.058 -0.625 0.797) 1.015 0.779 ( 0.105 -3.940 0.033) 3.941 0.801 ( -0.804 -0.424 0.201) 0.931 0.811 ( -0.517 -0.541 1.772) 1.924 0.832 ( 0.030 2.471 -0.918) 2.637 0.905 ( -0.605 -3.780 -2.788) 4.735 0.987 ( 0.805 1.012 -1.074) 1.681 1.033 ( -0.308 -2.460 -0.218) 2.489 1.073 ( 7.604 4.547 1.828) 9.046 1.087 ( -0.054 1.499 1.853) 2.384 1.129 ( 1.005 1.488 -0.165) 1.803 1.165 ( 0.626 -1.672 -0.294) 1.809 1.194 ( -1.589 -1.591 -0.219) 2.259 1.230 ( 3.146 -2.354 -0.857) 4.021 1.262 ( 0.858 -0.696 0.196) 1.122 1.297 ( -0.582 -2.095 0.799) 2.316 1.367 ( -2.724 1.352 -1.574) 3.424 1.373 ( -1.890 2.842 0.786) 3.502 1.404 ( -0.287 -0.888 -2.074) 2.274 1.423 ( 0.807 -0.706 1.006) 1.471 1.435 ( 0.131 0.656 0.018) 0.669 1.476 ( -0.998 -2.037 0.088) 2.270 1.545 ( -0.177 1.852 -0.097) 1.863 1.556 ( 0.016 1.171 -0.322) 1.215 1.602 ( 3.374 2.309 -0.751) 4.157 1.627 ( 2.546 0.118 -0.467) 2.591 1.678 ( 2.999 0.194 1.346) 3.293 1.711 ( 2.824 -2.340 0.310) 3.681 1.766 ( 0.656 1.106 0.269) 1.313 1.771 ( 0.683 -0.086 0.757) 1.023 1.803 ( -1.100 2.743 -0.306) 2.972 1.822 ( 0.568 -0.262 -0.052) 0.627 1.832 ( 1.236 -1.580 0.264) 2.023 1.870 ( 0.867 -3.461 -0.160) 3.571 2.021 ( -1.018 1.308 0.681) 1.792 2.036 ( -0.892 0.048 -0.724) 1.151 2.073 ( -0.509 0.714 0.255) 0.913 2.115 ( 0.157 -2.781 0.401) 2.814 2.229 ( -0.407 0.695 1.032) 1.309 2.241 ( -0.883 -0.431 0.088) 0.986 2.277 ( 1.696 0.946 0.510) 2.008 2.286 ( 1.119 0.000 -0.686) 1.313 2.328 ( -0.573 1.390 0.249) 1.524 2.368 ( -1.256 -1.469 0.256) 1.950 2.401 ( -2.306 2.052 -1.127) 3.286 2.416 ( -0.857 0.512 -0.568) 1.148 2.580 ( -3.237 0.660 0.594) 3.357 2.599 ( -2.324 -1.354 0.028) 2.690 2.679 ( -1.844 2.003 1.878) 3.308 2.711 ( -2.011 -1.037 -0.613) 2.344 5.482 ( -0.549 0.022 0.200) 0.585 5.484 ( -0.835 -0.146 -0.063) 0.850 5.498 ( -0.149 0.169 -0.071) 0.236 5.501 ( -0.156 -0.124 0.104) 0.225 6.298 ( 0.462 1.179 -0.589) 1.396 6.317 ( 1.276 -0.460 -0.376) 1.407 6.387 ( 1.423 0.394 0.560) 1.579 6.401 ( 1.000 -0.965 -0.790) 1.599 6.508 ( -1.325 0.984 -0.001) 1.651 6.520 ( -0.828 0.273 -0.595) 1.056 6.591 ( -1.839 0.359 1.243) 2.249 6.633 ( -0.462 -2.929 1.215) 3.204 6.883 ( -0.026 1.335 -0.518) 1.432 6.914 ( 1.318 -1.896 -1.061) 2.541 7.003 ( 0.536 2.921 -0.157) 2.974 7.050 ( 2.107 -1.750 -0.347) 2.761 ======================= Grid point 86 (35/36) ======================= q-point: ( 0.33 0.40 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 76 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.501 ( 3.862 2.062 -0.333) 4.390 0.525 ( 3.153 -0.064 -0.637) 3.217 0.552 ( 2.276 1.929 -0.666) 3.057 0.619 ( 5.957 -1.106 1.777) 6.314 0.645 ( 4.610 -2.275 -0.342) 5.153 0.662 ( 3.338 -2.664 -1.415) 4.499 0.717 ( -1.573 0.125 -1.803) 2.396 0.734 ( -2.812 -0.548 -1.546) 3.255 0.776 ( -1.744 1.164 1.400) 2.521 0.803 ( -0.572 -1.046 0.467) 1.280 0.850 ( 0.467 0.450 0.076) 0.653 0.896 ( -0.442 -3.215 -0.922) 3.373 0.971 ( -1.958 0.382 -2.231) 2.993 0.995 ( -1.997 -1.936 -0.806) 2.896 1.046 ( -3.431 0.624 0.279) 3.499 1.072 ( -2.530 -1.754 -1.337) 3.356 1.159 ( -1.739 -0.041 -1.057) 2.036 1.176 ( -0.601 -1.191 -1.011) 1.674 1.248 ( -0.624 1.324 -0.614) 1.587 1.256 ( -1.112 0.448 -1.135) 1.651 1.302 ( 6.060 3.805 0.532) 7.176 1.312 ( -2.650 1.453 -1.041) 3.197 1.347 ( -0.937 -0.201 -0.014) 0.959 1.373 ( 5.281 -1.187 2.334) 5.895 1.398 ( 4.123 3.427 0.563) 5.390 1.448 ( 2.134 -1.614 -1.525) 3.080 1.471 ( 3.188 1.392 -0.556) 3.523 1.528 ( 3.992 -3.116 -0.320) 5.074 1.535 ( 0.573 0.759 0.662) 1.159 1.549 ( 0.905 -0.533 -0.169) 1.064 1.674 ( 0.833 1.487 0.977) 1.965 1.688 ( 1.126 0.696 0.408) 1.385 1.708 ( 1.723 1.256 1.113) 2.405 1.731 ( -1.256 -1.643 1.187) 2.384 1.784 ( -0.255 1.598 0.813) 1.811 1.807 ( 1.832 0.713 0.960) 2.188 1.816 ( 1.316 0.525 0.333) 1.455 1.830 ( 0.427 0.128 -0.184) 0.482 1.837 ( -0.239 -0.438 -0.451) 0.673 1.884 ( 0.605 -3.361 0.228) 3.423 2.009 ( -0.062 1.680 0.359) 1.719 2.043 ( -1.897 -0.040 0.094) 1.900 2.049 ( 1.764 0.291 0.238) 1.804 2.080 ( -2.080 -1.092 0.516) 2.406 2.211 ( -0.264 1.402 0.078) 1.429 2.241 ( 0.265 -1.300 -0.094) 1.330 2.279 ( -1.206 -0.556 1.000) 1.663 2.296 ( -0.581 -1.923 -0.199) 2.018 2.331 ( -1.028 1.814 0.294) 2.105 2.343 ( -0.730 2.149 -0.488) 2.322 2.358 ( -0.450 -0.706 1.379) 1.613 2.391 ( -1.735 -1.577 -0.511) 2.400 2.515 ( -2.003 1.450 0.042) 2.473 2.548 ( -1.816 -2.209 -1.137) 3.077 2.626 ( -2.642 2.053 0.934) 3.474 2.659 ( -2.472 -1.261 -0.069) 2.776 5.466 ( -0.669 0.022 -0.127) 0.682 5.470 ( -0.507 -0.303 0.027) 0.591 5.494 ( -0.244 0.262 0.127) 0.380 5.496 ( -0.274 0.012 0.164) 0.319 6.312 ( 0.631 1.728 -0.659) 1.954 6.347 ( 1.441 -1.343 -0.334) 1.998 6.410 ( 0.064 0.121 -0.026) 0.139 6.417 ( 0.972 -0.352 -0.136) 1.043 6.475 ( -1.525 0.678 -0.099) 1.672 6.497 ( -1.137 -1.258 -0.573) 1.790 6.557 ( -1.231 2.921 1.365) 3.451 6.615 ( -0.922 -2.359 1.158) 2.785 6.892 ( 0.732 2.329 -1.137) 2.693 6.953 ( 1.878 -3.749 -1.327) 4.398 7.017 ( 0.549 3.927 0.123) 3.967 7.087 ( 1.184 -2.937 -0.074) 3.168 ======================= Grid point 87 (36/36) ======================= q-point: (-0.50 0.40 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.25e-05 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 Number of triplets: 50 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.562 ( -0.000 6.446 1.185) 6.554 0.562 ( 0.000 -0.455 -1.536) 1.602 0.580 ( -0.000 -1.906 -1.633) 2.510 0.665 ( -0.000 -0.177 0.651) 0.675 0.683 ( 0.000 -0.536 -1.930) 2.003 0.707 ( -0.000 -2.798 -1.299) 3.085 0.736 ( 0.000 3.796 2.136) 4.355 0.742 ( -0.000 -0.612 0.297) 0.680 0.749 ( -0.000 -1.563 0.870) 1.789 0.798 ( 0.000 -1.050 -2.910) 3.094 0.846 ( 0.000 -0.247 0.908) 0.941 0.873 ( 0.000 -2.125 -0.700) 2.237 0.944 ( 0.000 -0.728 1.146) 1.358 0.979 ( -0.000 0.785 -0.996) 1.268 0.982 ( 0.000 -3.653 -1.598) 3.987 1.022 ( 0.000 -3.265 -2.994) 4.430 1.139 ( -0.000 2.135 -2.679) 3.426 1.171 ( 0.000 -0.908 -2.650) 2.801 1.238 ( 0.000 1.206 -0.421) 1.277 1.249 ( 0.000 0.106 -0.246) 0.268 1.275 ( 0.000 0.527 -0.093) 0.535 1.289 ( 0.000 -0.857 -0.831) 1.193 1.404 ( 0.000 1.502 3.295) 3.622 1.425 ( -0.000 -0.831 3.055) 3.166 1.467 ( 0.000 1.348 -2.258) 2.630 1.488 ( 0.000 -0.700 -2.712) 2.801 1.533 ( -0.000 3.136 0.694) 3.212 1.568 ( -0.000 -0.230 0.806) 0.838 1.574 ( 0.000 2.270 0.444) 2.313 1.617 ( 0.000 3.121 -0.121) 3.123 1.622 ( -0.000 -2.335 0.904) 2.503 1.660 ( 0.000 4.298 -0.083) 4.299 1.693 ( 0.000 -4.078 0.100) 4.079 1.758 ( 0.000 -4.847 0.600) 4.884 1.772 ( 0.000 3.223 0.795) 3.320 1.819 ( 0.000 -1.211 1.281) 1.763 1.835 ( 0.000 0.361 0.175) 0.401 1.837 ( 0.000 0.514 0.339) 0.615 1.839 ( 0.000 2.847 1.274) 3.119 1.887 ( 0.000 -1.730 1.375) 2.210 2.013 ( 0.000 0.691 -0.438) 0.818 2.021 ( -0.000 -0.029 -0.038) 0.048 2.054 ( 0.000 1.931 0.279) 1.951 2.077 ( 0.000 -0.102 0.343) 0.358 2.214 ( -0.000 1.043 -0.183) 1.058 2.245 ( -0.000 1.050 0.490) 1.159 2.253 ( 0.000 -2.866 0.086) 2.868 2.288 ( 0.000 -2.726 0.026) 2.726 2.321 ( 0.000 0.550 0.413) 0.688 2.336 ( 0.000 1.093 0.235) 1.118 2.347 ( -0.000 -0.681 0.449) 0.815 2.355 ( 0.000 -0.704 0.364) 0.792 2.500 ( 0.000 1.174 0.861) 1.456 2.536 ( 0.000 -2.367 -2.337) 3.327 2.592 ( 0.000 1.961 0.760) 2.103 2.626 ( 0.000 -1.473 0.148) 1.480 5.459 ( 0.000 0.078 -0.148) 0.167 5.464 ( 0.000 -0.430 -0.096) 0.441 5.490 ( 0.000 0.287 0.237) 0.372 5.493 ( 0.000 0.080 0.215) 0.229 6.319 ( 0.000 1.975 -0.664) 2.083 6.368 ( 0.000 -2.589 -0.639) 2.667 6.405 ( -0.000 1.692 0.427) 1.745 6.433 ( 0.000 -1.194 -0.497) 1.293 6.453 ( 0.000 0.995 -0.360) 1.058 6.483 ( 0.000 -2.002 -0.121) 2.005 6.543 ( 0.000 3.673 1.378) 3.923 6.604 ( 0.000 -2.132 1.165) 2.430 6.902 ( 0.000 2.895 -1.321) 3.182 6.976 ( 0.000 -4.449 -1.417) 4.669 7.023 ( 0.000 4.277 0.230) 4.283 7.099 ( 0.000 -3.277 0.108) 3.279 =================== 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/10800 10.0 2.437 2.139 1.480 0.000 0.000 0.000 3/10800 20.0 0.849 0.781 0.530 0.000 0.000 0.000 3/10800 30.0 0.525 0.486 0.327 0.000 0.000 0.000 3/10800 40.0 0.386 0.359 0.239 0.000 0.000 0.000 3/10800 50.0 0.308 0.289 0.191 0.000 0.000 0.000 3/10800 60.0 0.258 0.244 0.160 0.000 0.000 0.000 3/10800 70.0 0.223 0.213 0.138 0.000 0.000 0.000 3/10800 80.0 0.196 0.190 0.122 0.000 0.000 0.000 3/10800 90.0 0.176 0.172 0.110 0.000 0.000 0.000 3/10800 100.0 0.160 0.158 0.100 0.000 0.000 0.000 3/10800 110.0 0.146 0.146 0.091 0.000 0.000 0.000 3/10800 120.0 0.135 0.135 0.084 0.000 0.000 0.000 3/10800 130.0 0.125 0.126 0.078 0.000 0.000 0.000 3/10800 140.0 0.116 0.119 0.073 0.000 0.000 0.000 3/10800 150.0 0.109 0.112 0.069 0.000 0.000 0.000 3/10800 160.0 0.103 0.105 0.065 0.000 0.000 0.000 3/10800 170.0 0.097 0.100 0.061 0.000 0.000 0.000 3/10800 180.0 0.092 0.095 0.058 0.000 0.000 0.000 3/10800 190.0 0.087 0.090 0.055 0.000 0.000 0.000 3/10800 200.0 0.083 0.086 0.052 0.000 0.000 0.000 3/10800 210.0 0.079 0.082 0.050 0.000 0.000 0.000 3/10800 220.0 0.076 0.079 0.048 0.000 0.000 0.000 3/10800 230.0 0.073 0.076 0.046 0.000 0.000 0.000 3/10800 240.0 0.070 0.073 0.044 0.000 0.000 0.000 3/10800 250.0 0.067 0.070 0.042 0.000 0.000 0.000 3/10800 260.0 0.064 0.068 0.041 0.000 0.000 0.000 3/10800 270.0 0.062 0.065 0.039 0.000 0.000 0.000 3/10800 280.0 0.060 0.063 0.038 0.000 0.000 0.000 3/10800 290.0 0.058 0.061 0.037 0.000 0.000 0.000 3/10800 300.0 0.056 0.059 0.035 0.000 0.000 0.000 3/10800 310.0 0.054 0.057 0.034 0.000 0.000 0.000 3/10800 320.0 0.053 0.055 0.033 0.000 0.000 0.000 3/10800 330.0 0.051 0.054 0.032 0.000 0.000 0.000 3/10800 340.0 0.050 0.052 0.031 0.000 0.000 0.000 3/10800 350.0 0.048 0.051 0.030 0.000 0.000 0.000 3/10800 360.0 0.047 0.050 0.030 0.000 0.000 0.000 3/10800 370.0 0.046 0.048 0.029 0.000 0.000 0.000 3/10800 380.0 0.044 0.047 0.028 0.000 0.000 0.000 3/10800 390.0 0.043 0.046 0.027 0.000 0.000 0.000 3/10800 400.0 0.042 0.045 0.027 0.000 0.000 0.000 3/10800 410.0 0.041 0.044 0.026 0.000 0.000 0.000 3/10800 420.0 0.040 0.043 0.025 0.000 0.000 0.000 3/10800 430.0 0.039 0.042 0.025 0.000 0.000 0.000 3/10800 440.0 0.038 0.041 0.024 0.000 0.000 0.000 3/10800 450.0 0.038 0.040 0.024 0.000 0.000 0.000 3/10800 460.0 0.037 0.039 0.023 0.000 0.000 0.000 3/10800 470.0 0.036 0.038 0.023 0.000 0.000 0.000 3/10800 480.0 0.035 0.037 0.022 0.000 0.000 0.000 3/10800 490.0 0.035 0.037 0.022 0.000 0.000 0.000 3/10800 500.0 0.034 0.036 0.021 0.000 0.000 0.000 3/10800 510.0 0.033 0.035 0.021 0.000 0.000 0.000 3/10800 520.0 0.033 0.035 0.021 0.000 0.000 0.000 3/10800 530.0 0.032 0.034 0.020 0.000 0.000 0.000 3/10800 540.0 0.031 0.033 0.020 0.000 0.000 0.000 3/10800 550.0 0.031 0.033 0.020 0.000 0.000 0.000 3/10800 560.0 0.030 0.032 0.019 0.000 0.000 0.000 3/10800 570.0 0.030 0.032 0.019 0.000 0.000 0.000 3/10800 580.0 0.029 0.031 0.019 0.000 0.000 0.000 3/10800 590.0 0.029 0.031 0.018 0.000 0.000 0.000 3/10800 600.0 0.028 0.030 0.018 0.000 0.000 0.000 3/10800 610.0 0.028 0.030 0.018 0.000 0.000 0.000 3/10800 620.0 0.027 0.029 0.017 0.000 0.000 0.000 3/10800 630.0 0.027 0.029 0.017 0.000 0.000 0.000 3/10800 640.0 0.027 0.028 0.017 0.000 0.000 0.000 3/10800 650.0 0.026 0.028 0.017 0.000 0.000 0.000 3/10800 660.0 0.026 0.027 0.016 0.000 0.000 0.000 3/10800 670.0 0.025 0.027 0.016 0.000 0.000 0.000 3/10800 680.0 0.025 0.027 0.016 0.000 0.000 0.000 3/10800 690.0 0.025 0.026 0.016 0.000 0.000 0.000 3/10800 700.0 0.024 0.026 0.015 0.000 0.000 0.000 3/10800 710.0 0.024 0.025 0.015 0.000 0.000 0.000 3/10800 720.0 0.024 0.025 0.015 0.000 0.000 0.000 3/10800 730.0 0.023 0.025 0.015 0.000 0.000 0.000 3/10800 740.0 0.023 0.024 0.015 0.000 0.000 0.000 3/10800 750.0 0.023 0.024 0.014 0.000 0.000 0.000 3/10800 760.0 0.022 0.024 0.014 0.000 0.000 0.000 3/10800 770.0 0.022 0.024 0.014 0.000 0.000 0.000 3/10800 780.0 0.022 0.023 0.014 0.000 0.000 0.000 3/10800 790.0 0.022 0.023 0.014 0.000 0.000 0.000 3/10800 800.0 0.021 0.023 0.013 0.000 0.000 0.000 3/10800 810.0 0.021 0.022 0.013 0.000 0.000 0.000 3/10800 820.0 0.021 0.022 0.013 0.000 0.000 0.000 3/10800 830.0 0.021 0.022 0.013 0.000 0.000 0.000 3/10800 840.0 0.020 0.022 0.013 0.000 0.000 0.000 3/10800 850.0 0.020 0.021 0.013 0.000 0.000 0.000 3/10800 860.0 0.020 0.021 0.013 0.000 0.000 0.000 3/10800 870.0 0.020 0.021 0.012 0.000 0.000 0.000 3/10800 880.0 0.019 0.021 0.012 0.000 0.000 0.000 3/10800 890.0 0.019 0.020 0.012 0.000 0.000 0.000 3/10800 900.0 0.019 0.020 0.012 0.000 0.000 0.000 3/10800 910.0 0.019 0.020 0.012 0.000 0.000 0.000 3/10800 920.0 0.019 0.020 0.012 0.000 0.000 0.000 3/10800 930.0 0.018 0.020 0.012 0.000 0.000 0.000 3/10800 940.0 0.018 0.019 0.011 0.000 0.000 0.000 3/10800 950.0 0.018 0.019 0.011 0.000 0.000 0.000 3/10800 960.0 0.018 0.019 0.011 0.000 0.000 0.000 3/10800 970.0 0.018 0.019 0.011 0.000 0.000 0.000 3/10800 980.0 0.017 0.019 0.011 0.000 0.000 0.000 3/10800 990.0 0.017 0.018 0.011 0.000 0.000 0.000 3/10800 1000.0 0.017 0.018 0.011 0.000 0.000 0.000 3/10800 Thermal conductivity related properties were written into "kappa-m655.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 17:57:20]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|