----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 07:20:04]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [1 1 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) Number of symmetry operations in supercell: 36 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 6.612214075966644 3.817563576698766 5.482577393333327 b -6.612214075966630 3.817563576698766 5.482577393333327 c 0.000000000000007 -7.635127153397544 5.482577393333327 Atomic positions (fractional): *1 Br 0.57259840926760 0.44021628734293 0.80474043540191 79.904 2 Br 0.44021628734293 0.80474043540191 0.57259840926760 79.904 3 Br 0.80474043540191 0.57259840926760 0.44021628734293 79.904 4 Br 0.69525956459809 0.05978371265707 0.92740159073240 79.904 5 Br 0.05978371265707 0.92740159073240 0.69525956459809 79.904 6 Br 0.92740159073240 0.69525956459809 0.05978371265707 79.904 7 Br 0.42740159073240 0.55978371265707 0.19525956459809 79.904 8 Br 0.55978371265707 0.19525956459809 0.42740159073240 79.904 9 Br 0.19525956459809 0.42740159073240 0.55978371265707 79.904 10 Br 0.30474043540191 0.94021628734293 0.07259840926760 79.904 11 Br 0.07259840926760 0.30474043540191 0.94021628734293 79.904 12 Br 0.94021628734293 0.07259840926760 0.30474043540191 79.904 *13 Rb 0.25000000000000 0.25000000000000 0.25000000000000 85.468 14 Rb 0.75000000000000 0.75000000000000 0.75000000000000 85.468 *15 Rb 0.13124390165606 0.75000000000000 0.36875609834394 85.468 16 Rb 0.75000000000000 0.36875609834394 0.13124390165606 85.468 17 Rb 0.36875609834394 0.13124390165606 0.75000000000000 85.468 18 Rb 0.86875609834394 0.25000000000000 0.63124390165606 85.468 19 Rb 0.25000000000000 0.63124390165606 0.86875609834394 85.468 20 Rb 0.63124390165606 0.86875609834394 0.25000000000000 85.468 *21 Pb 0.00000000000000 0.00000000000000 0.00000000000000 207.200 22 Pb 0.50000000000000 0.50000000000000 0.50000000000000 207.200 -------------------------------- unit cell --------------------------------- Lattice vectors: a 13.224428151933274 0.000000000000000 0.000000000000000 b -6.612214075966637 11.452690730096311 0.000000000000000 c 0.000000000000000 0.000000000000000 16.447732179999999 Atomic positions (fractional): *1 Br 0.30008003193012 0.46777794193557 0.27251837733748 79.904 > 1 2 Br 0.16769791000545 0.69991996806988 0.27251837733748 79.904 > 2 3 Br 0.53222205806443 0.83230208999455 0.27251837733748 79.904 > 3 4 Br 0.46777794193557 0.30008003193012 0.22748162266252 79.904 > 4 5 Br 0.83230208999455 0.53222205806443 0.22748162266252 79.904 > 5 6 Br 0.69991996806988 0.16769791000545 0.22748162266252 79.904 > 6 7 Br 0.36658663473655 0.86555539139777 0.06081495599585 79.904 > 7 8 Br 0.49896875666122 0.63341336526345 0.06081495599585 79.904 > 8 9 Br 0.13444460860224 0.50103124333878 0.06081495599585 79.904 > 9 10 Br 0.19888872473110 0.03325330140322 0.10585171067081 79.904 > 10 11 Br 0.96674669859678 0.16563542332788 0.10585171067081 79.904 > 11 12 Br 0.83436457667212 0.80111127526890 0.10585171067081 79.904 > 12 13 Br 0.96674669859678 0.80111127526890 0.60585171067081 79.904 > 1 14 Br 0.83436457667212 0.03325330140322 0.60585171067081 79.904 > 2 15 Br 0.19888872473110 0.16563542332788 0.60585171067081 79.904 > 3 16 Br 0.13444460860224 0.63341336526345 0.56081495599585 79.904 > 4 17 Br 0.49896875666122 0.86555539139777 0.56081495599585 79.904 > 5 18 Br 0.36658663473655 0.50103124333878 0.56081495599585 79.904 > 6 19 Br 0.03325330140322 0.19888872473110 0.39414828932919 79.904 > 7 20 Br 0.16563542332788 0.96674669859678 0.39414828932919 79.904 > 8 21 Br 0.80111127526890 0.83436457667212 0.39414828932919 79.904 > 9 22 Br 0.86555539139777 0.36658663473655 0.43918504400414 79.904 > 10 23 Br 0.63341336526345 0.49896875666122 0.43918504400414 79.904 > 11 24 Br 0.50103124333878 0.13444460860223 0.43918504400414 79.904 > 12 25 Br 0.63341336526345 0.13444460860224 0.93918504400415 79.904 > 1 26 Br 0.50103124333878 0.36658663473655 0.93918504400415 79.904 > 2 27 Br 0.86555539139777 0.49896875666122 0.93918504400415 79.904 > 3 28 Br 0.80111127526890 0.96674669859678 0.89414828932919 79.904 > 4 29 Br 0.16563542332788 0.19888872473110 0.89414828932919 79.904 > 5 30 Br 0.03325330140322 0.83436457667212 0.89414828932919 79.904 > 6 31 Br 0.69991996806988 0.53222205806443 0.72748162266252 79.904 > 7 32 Br 0.83230208999455 0.30008003193012 0.72748162266252 79.904 > 8 33 Br 0.46777794193557 0.16769791000545 0.72748162266252 79.904 > 9 34 Br 0.53222205806443 0.69991996806988 0.77251837733748 79.904 > 10 35 Br 0.30008003193012 0.83230208999455 0.77251837733748 79.904 > 11 36 Br 0.16769791000545 0.46777794193557 0.77251837733748 79.904 > 12 *37 Rb 0.33333333333333 0.66666666666667 0.91666666666667 85.468 > 13 38 Rb 0.33333333333333 0.66666666666667 0.41666666666667 85.468 > 14 *39 Rb 0.04791056832273 0.71457723498940 0.08333333333333 85.468 > 15 40 Rb 0.66666666666667 0.95208943167727 0.08333333333333 85.468 > 16 41 Rb 0.28542276501060 0.33333333333333 0.08333333333333 85.468 > 17 42 Rb 0.61875609834394 0.61875609834394 0.25000000000000 85.468 > 18 43 Rb 0.00000000000000 0.38124390165606 0.25000000000000 85.468 > 19 44 Rb 0.38124390165606 0.00000000000000 0.25000000000000 85.468 > 20 45 Rb 0.00000000000000 0.00000000000000 0.25000000000000 85.468 > 13 46 Rb 0.00000000000000 0.00000000000000 0.75000000000000 85.468 > 14 47 Rb 0.71457723498940 0.04791056832273 0.41666666666667 85.468 > 15 48 Rb 0.33333333333333 0.28542276501060 0.41666666666667 85.468 > 16 49 Rb 0.95208943167727 0.66666666666667 0.41666666666667 85.468 > 17 50 Rb 0.28542276501060 0.95208943167727 0.58333333333333 85.468 > 18 51 Rb 0.66666666666667 0.71457723498940 0.58333333333333 85.468 > 19 52 Rb 0.04791056832273 0.33333333333333 0.58333333333333 85.468 > 20 53 Rb 0.66666666666667 0.33333333333333 0.58333333333333 85.468 > 13 54 Rb 0.66666666666667 0.33333333333333 0.08333333333333 85.468 > 14 55 Rb 0.38124390165606 0.38124390165606 0.75000000000000 85.468 > 15 56 Rb 0.00000000000000 0.61875609834394 0.75000000000000 85.468 > 16 57 Rb 0.61875609834394 0.00000000000000 0.75000000000000 85.468 > 17 58 Rb 0.95208943167727 0.28542276501060 0.91666666666667 85.468 > 18 59 Rb 0.33333333333333 0.04791056832273 0.91666666666667 85.468 > 19 60 Rb 0.71457723498940 0.66666666666667 0.91666666666667 85.468 > 20 *61 Pb 0.00000000000000 0.00000000000000 0.00000000000000 207.200 > 21 62 Pb 0.33333333333333 0.66666666666667 0.16666666666667 207.200 > 22 63 Pb 0.66666666666667 0.33333333333333 0.33333333333333 207.200 > 21 64 Pb 0.00000000000000 0.00000000000000 0.50000000000000 207.200 > 22 65 Pb 0.33333333333333 0.66666666666667 0.66666666666667 207.200 > 21 66 Pb 0.66666666666667 0.33333333333333 0.83333333333333 207.200 > 22 -------------------------------- super cell -------------------------------- Lattice vectors: a 13.224428151933274 0.000000000000000 0.000000000000000 b -6.612214075966637 11.452690730096311 0.000000000000000 c 0.000000000000000 0.000000000000000 16.447732179999999 Atomic positions (fractional): *1 Br 0.30008003193012 0.46777794193557 0.27251837733748 79.904 > 1 2 Br 0.16769791000545 0.69991996806988 0.27251837733748 79.904 > 2 3 Br 0.53222205806443 0.83230208999455 0.27251837733748 79.904 > 3 4 Br 0.46777794193557 0.30008003193012 0.22748162266252 79.904 > 4 5 Br 0.83230208999455 0.53222205806443 0.22748162266252 79.904 > 5 6 Br 0.69991996806988 0.16769791000545 0.22748162266252 79.904 > 6 7 Br 0.36658663473655 0.86555539139777 0.06081495599585 79.904 > 7 8 Br 0.49896875666122 0.63341336526345 0.06081495599585 79.904 > 8 9 Br 0.13444460860224 0.50103124333878 0.06081495599585 79.904 > 9 10 Br 0.19888872473110 0.03325330140322 0.10585171067081 79.904 > 10 11 Br 0.96674669859678 0.16563542332788 0.10585171067081 79.904 > 11 12 Br 0.83436457667212 0.80111127526890 0.10585171067081 79.904 > 12 13 Br 0.96674669859678 0.80111127526890 0.60585171067081 79.904 > 1 14 Br 0.83436457667212 0.03325330140322 0.60585171067081 79.904 > 2 15 Br 0.19888872473110 0.16563542332788 0.60585171067081 79.904 > 3 16 Br 0.13444460860224 0.63341336526345 0.56081495599585 79.904 > 4 17 Br 0.49896875666122 0.86555539139777 0.56081495599585 79.904 > 5 18 Br 0.36658663473655 0.50103124333878 0.56081495599585 79.904 > 6 19 Br 0.03325330140322 0.19888872473110 0.39414828932919 79.904 > 7 20 Br 0.16563542332788 0.96674669859678 0.39414828932919 79.904 > 8 21 Br 0.80111127526890 0.83436457667212 0.39414828932919 79.904 > 9 22 Br 0.86555539139777 0.36658663473655 0.43918504400414 79.904 > 10 23 Br 0.63341336526345 0.49896875666122 0.43918504400414 79.904 > 11 24 Br 0.50103124333878 0.13444460860223 0.43918504400414 79.904 > 12 25 Br 0.63341336526345 0.13444460860224 0.93918504400415 79.904 > 1 26 Br 0.50103124333878 0.36658663473655 0.93918504400415 79.904 > 2 27 Br 0.86555539139777 0.49896875666122 0.93918504400415 79.904 > 3 28 Br 0.80111127526890 0.96674669859678 0.89414828932919 79.904 > 4 29 Br 0.16563542332788 0.19888872473110 0.89414828932919 79.904 > 5 30 Br 0.03325330140322 0.83436457667212 0.89414828932919 79.904 > 6 31 Br 0.69991996806988 0.53222205806443 0.72748162266252 79.904 > 7 32 Br 0.83230208999455 0.30008003193012 0.72748162266252 79.904 > 8 33 Br 0.46777794193557 0.16769791000545 0.72748162266252 79.904 > 9 34 Br 0.53222205806443 0.69991996806988 0.77251837733748 79.904 > 10 35 Br 0.30008003193012 0.83230208999455 0.77251837733748 79.904 > 11 36 Br 0.16769791000545 0.46777794193557 0.77251837733748 79.904 > 12 *37 Rb 0.33333333333333 0.66666666666667 0.91666666666667 85.468 > 13 38 Rb 0.33333333333333 0.66666666666667 0.41666666666667 85.468 > 14 *39 Rb 0.04791056832273 0.71457723498940 0.08333333333333 85.468 > 15 40 Rb 0.66666666666667 0.95208943167727 0.08333333333333 85.468 > 16 41 Rb 0.28542276501060 0.33333333333333 0.08333333333333 85.468 > 17 42 Rb 0.61875609834394 0.61875609834394 0.25000000000000 85.468 > 18 43 Rb 0.00000000000000 0.38124390165606 0.25000000000000 85.468 > 19 44 Rb 0.38124390165606 0.00000000000000 0.25000000000000 85.468 > 20 45 Rb 0.00000000000000 0.00000000000000 0.25000000000000 85.468 > 13 46 Rb 0.00000000000000 0.00000000000000 0.75000000000000 85.468 > 14 47 Rb 0.71457723498940 0.04791056832273 0.41666666666667 85.468 > 15 48 Rb 0.33333333333333 0.28542276501060 0.41666666666667 85.468 > 16 49 Rb 0.95208943167727 0.66666666666667 0.41666666666667 85.468 > 17 50 Rb 0.28542276501060 0.95208943167727 0.58333333333333 85.468 > 18 51 Rb 0.66666666666667 0.71457723498940 0.58333333333333 85.468 > 19 52 Rb 0.04791056832273 0.33333333333333 0.58333333333333 85.468 > 20 53 Rb 0.66666666666667 0.33333333333333 0.58333333333333 85.468 > 13 54 Rb 0.66666666666667 0.33333333333333 0.08333333333333 85.468 > 14 55 Rb 0.38124390165606 0.38124390165606 0.75000000000000 85.468 > 15 56 Rb 0.00000000000000 0.61875609834394 0.75000000000000 85.468 > 16 57 Rb 0.61875609834394 0.00000000000000 0.75000000000000 85.468 > 17 58 Rb 0.95208943167727 0.28542276501060 0.91666666666667 85.468 > 18 59 Rb 0.33333333333333 0.04791056832273 0.91666666666667 85.468 > 19 60 Rb 0.71457723498940 0.66666666666667 0.91666666666667 85.468 > 20 *61 Pb 0.00000000000000 0.00000000000000 0.00000000000000 207.200 > 21 62 Pb 0.33333333333333 0.66666666666667 0.16666666666667 207.200 > 22 63 Pb 0.66666666666667 0.33333333333333 0.33333333333333 207.200 > 21 64 Pb 0.00000000000000 0.00000000000000 0.50000000000000 207.200 > 22 65 Pb 0.33333333333333 0.66666666666667 0.66666666666667 207.200 > 21 66 Pb 0.66666666666667 0.33333333333333 0.83333333333333 207.200 > 22 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.1828860 0.0000000 0.0000000 -0.0000000 3.1828860 0.0000000 0.0000000 0.0000000 3.1968694 -------------------------- Born effective charges -------------------------- 1 Br -1.0375216 0.3229362 -0.2279541 0.2388177 -1.5573504 0.5110383 -0.2313857 0.5085382 -1.3620732 2 Br -1.6706398 0.1267132 0.5565492 0.0425947 -0.9242322 -0.0581051 0.5560999 -0.0538832 -1.3620732 3 Br -1.1841466 -0.3234717 -0.3285951 -0.4075902 -1.4107254 -0.4529332 -0.3247142 -0.4546550 -1.3620732 4 Br -1.6706398 -0.1267132 -0.5565492 -0.0425947 -0.9242322 -0.0581051 -0.5560999 -0.0538832 -1.3620732 5 Br -1.0375216 -0.3229362 0.2279541 -0.2388177 -1.5573504 0.5110383 0.2313857 0.5085382 -1.3620732 6 Br -1.1841466 0.3234717 0.3285951 0.4075902 -1.4107254 -0.4529332 0.3247142 -0.4546550 -1.3620732 7 Br -1.0375216 0.3229362 -0.2279541 0.2388177 -1.5573504 0.5110383 -0.2313857 0.5085382 -1.3620732 8 Br -1.6706398 0.1267132 0.5565492 0.0425947 -0.9242322 -0.0581051 0.5560999 -0.0538832 -1.3620732 9 Br -1.1841466 -0.3234717 -0.3285951 -0.4075902 -1.4107254 -0.4529332 -0.3247142 -0.4546550 -1.3620732 10 Br -1.6706398 -0.1267132 -0.5565492 -0.0425947 -0.9242322 -0.0581051 -0.5560999 -0.0538832 -1.3620732 11 Br -1.1841466 0.3234717 0.3285951 0.4075902 -1.4107254 -0.4529332 0.3247142 -0.4546550 -1.3620732 12 Br -1.0375216 -0.3229362 0.2279541 -0.2388177 -1.5573504 0.5110383 0.2313857 0.5085382 -1.3620732 13 Rb 0.9094608 0.0000000 -0.0000000 -0.0000000 0.9094608 0.0000000 0.0000000 -0.0000000 1.6152910 14 Rb 0.9094608 0.0000000 -0.0000000 -0.0000000 0.9094608 0.0000000 0.0000000 -0.0000000 1.6152910 15 Rb 1.3573660 -0.0688557 0.1709532 -0.0688557 1.2778583 -0.0986999 0.1363511 -0.0787224 1.1663259 16 Rb 1.3573660 0.0688557 -0.1709532 0.0688557 1.2778583 -0.0986999 -0.1363511 -0.0787224 1.1663259 17 Rb 1.2381044 -0.0000000 -0.0000000 0.0000000 1.3971199 0.1973997 0.0000000 0.1574447 1.1663259 18 Rb 1.3573660 -0.0688557 0.1709532 -0.0688557 1.2778583 -0.0986999 0.1363511 -0.0787224 1.1663259 19 Rb 1.3573660 0.0688557 -0.1709532 0.0688557 1.2778583 -0.0986999 -0.1363511 -0.0787224 1.1663259 20 Rb 1.2381044 -0.0000000 -0.0000000 0.0000000 1.3971199 0.1973997 0.0000000 0.1574447 1.1663259 21 Pb 2.9223190 0.0862955 -0.0000000 -0.0862955 2.9223190 0.0000000 0.0000000 -0.0000000 3.0581705 22 Pb 2.9223190 -0.0862955 -0.0000000 0.0862955 2.9223190 0.0000000 0.0000000 -0.0000000 3.0581705 ---------------------------------------------------------------------------- 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: 198/198 Permutation basis: 8409/8409 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 592 Number of blocks in projector: 300 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 275 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 234 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 83 Use standard eigh solver. Tree of FC basis block matrices: - (592, 573), data: False |-- (83, 79), data: True |-- (234, 228), data: True |-- (275, 266), data: True ----- Solver_atoms: 1 -- 66 / 66 Time (Solver_compr_matrix_reshape): 0.002 Solver_block: 100 / 240 - Time: 0.299 Solver_block: 200 / 240 - Time: 0.317 Solver_block: 240 / 240 - Time: 0.145 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.774 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 198/198 Permutation basis: 8409/8409 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 592 Number of blocks in projector: 300 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 275 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 234 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 83 Use standard eigh solver. Tree of FC basis block matrices: - (592, 573), data: False |-- (83, 79), data: True |-- (234, 228), data: True |-- (275, 266), data: True Max drift after symmetrization by symfc projector: 0.00000000 (yy) 0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 07:20:11]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 07:20:12]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [1 1 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 6.612214075966644 3.817563576698766 5.482577393333327 b -6.612214075966630 3.817563576698766 5.482577393333327 c 0.000000000000007 -7.635127153397544 5.482577393333327 Atomic positions (fractional): 1 Br 0.57259840926760 0.44021628734293 0.80474043540191 79.904 2 Br 0.44021628734293 0.80474043540191 0.57259840926760 79.904 3 Br 0.80474043540191 0.57259840926760 0.44021628734293 79.904 4 Br 0.69525956459809 0.05978371265707 0.92740159073240 79.904 5 Br 0.05978371265707 0.92740159073240 0.69525956459809 79.904 6 Br 0.92740159073240 0.69525956459809 0.05978371265707 79.904 7 Br 0.42740159073240 0.55978371265707 0.19525956459809 79.904 8 Br 0.55978371265707 0.19525956459809 0.42740159073240 79.904 9 Br 0.19525956459809 0.42740159073240 0.55978371265707 79.904 10 Br 0.30474043540191 0.94021628734293 0.07259840926760 79.904 11 Br 0.07259840926760 0.30474043540191 0.94021628734293 79.904 12 Br 0.94021628734293 0.07259840926760 0.30474043540191 79.904 13 Rb 0.25000000000000 0.25000000000000 0.25000000000000 85.468 14 Rb 0.75000000000000 0.75000000000000 0.75000000000000 85.468 15 Rb 0.13124390165606 0.75000000000000 0.36875609834394 85.468 16 Rb 0.75000000000000 0.36875609834394 0.13124390165606 85.468 17 Rb 0.36875609834394 0.13124390165606 0.75000000000000 85.468 18 Rb 0.86875609834394 0.25000000000000 0.63124390165606 85.468 19 Rb 0.25000000000000 0.63124390165606 0.86875609834394 85.468 20 Rb 0.63124390165606 0.86875609834394 0.25000000000000 85.468 21 Pb 0.00000000000000 0.00000000000000 0.00000000000000 207.200 22 Pb 0.50000000000000 0.50000000000000 0.50000000000000 207.200 -------------------------------- supercell --------------------------------- Lattice vectors: a 13.224428151933274 0.000000000000000 0.000000000000000 b -6.612214075966637 11.452690730096311 0.000000000000000 c 0.000000000000000 0.000000000000000 16.447732179999999 Atomic positions (fractional): 1 Br 0.30008003193012 0.46777794193557 0.27251837733748 79.904 > 1 2 Br 0.16769791000545 0.69991996806988 0.27251837733748 79.904 > 2 3 Br 0.53222205806443 0.83230208999455 0.27251837733748 79.904 > 3 4 Br 0.46777794193557 0.30008003193012 0.22748162266252 79.904 > 4 5 Br 0.83230208999455 0.53222205806443 0.22748162266252 79.904 > 5 6 Br 0.69991996806988 0.16769791000545 0.22748162266252 79.904 > 6 7 Br 0.36658663473655 0.86555539139777 0.06081495599585 79.904 > 7 8 Br 0.49896875666122 0.63341336526345 0.06081495599585 79.904 > 8 9 Br 0.13444460860224 0.50103124333878 0.06081495599585 79.904 > 9 10 Br 0.19888872473110 0.03325330140322 0.10585171067081 79.904 > 10 11 Br 0.96674669859678 0.16563542332788 0.10585171067081 79.904 > 11 12 Br 0.83436457667212 0.80111127526890 0.10585171067081 79.904 > 12 13 Br 0.96674669859678 0.80111127526890 0.60585171067081 79.904 > 1 14 Br 0.83436457667212 0.03325330140322 0.60585171067081 79.904 > 2 15 Br 0.19888872473110 0.16563542332788 0.60585171067081 79.904 > 3 16 Br 0.13444460860224 0.63341336526345 0.56081495599585 79.904 > 4 17 Br 0.49896875666122 0.86555539139777 0.56081495599585 79.904 > 5 18 Br 0.36658663473655 0.50103124333878 0.56081495599585 79.904 > 6 19 Br 0.03325330140322 0.19888872473110 0.39414828932919 79.904 > 7 20 Br 0.16563542332788 0.96674669859678 0.39414828932919 79.904 > 8 21 Br 0.80111127526890 0.83436457667212 0.39414828932919 79.904 > 9 22 Br 0.86555539139777 0.36658663473655 0.43918504400414 79.904 > 10 23 Br 0.63341336526345 0.49896875666122 0.43918504400414 79.904 > 11 24 Br 0.50103124333878 0.13444460860223 0.43918504400414 79.904 > 12 25 Br 0.63341336526345 0.13444460860224 0.93918504400415 79.904 > 1 26 Br 0.50103124333878 0.36658663473655 0.93918504400415 79.904 > 2 27 Br 0.86555539139777 0.49896875666122 0.93918504400415 79.904 > 3 28 Br 0.80111127526890 0.96674669859678 0.89414828932919 79.904 > 4 29 Br 0.16563542332788 0.19888872473110 0.89414828932919 79.904 > 5 30 Br 0.03325330140322 0.83436457667212 0.89414828932919 79.904 > 6 31 Br 0.69991996806988 0.53222205806443 0.72748162266252 79.904 > 7 32 Br 0.83230208999455 0.30008003193012 0.72748162266252 79.904 > 8 33 Br 0.46777794193557 0.16769791000545 0.72748162266252 79.904 > 9 34 Br 0.53222205806443 0.69991996806988 0.77251837733748 79.904 > 10 35 Br 0.30008003193012 0.83230208999455 0.77251837733748 79.904 > 11 36 Br 0.16769791000545 0.46777794193557 0.77251837733748 79.904 > 12 37 Rb 0.33333333333333 0.66666666666667 0.91666666666667 85.468 > 37 38 Rb 0.33333333333333 0.66666666666667 0.41666666666667 85.468 > 38 39 Rb 0.04791056832273 0.71457723498940 0.08333333333333 85.468 > 39 40 Rb 0.66666666666667 0.95208943167727 0.08333333333333 85.468 > 40 41 Rb 0.28542276501060 0.33333333333333 0.08333333333333 85.468 > 41 42 Rb 0.61875609834394 0.61875609834394 0.25000000000000 85.468 > 42 43 Rb 0.00000000000000 0.38124390165606 0.25000000000000 85.468 > 43 44 Rb 0.38124390165606 0.00000000000000 0.25000000000000 85.468 > 44 45 Rb 0.00000000000000 0.00000000000000 0.25000000000000 85.468 > 37 46 Rb 0.00000000000000 0.00000000000000 0.75000000000000 85.468 > 38 47 Rb 0.71457723498940 0.04791056832273 0.41666666666667 85.468 > 39 48 Rb 0.33333333333333 0.28542276501060 0.41666666666667 85.468 > 40 49 Rb 0.95208943167727 0.66666666666667 0.41666666666667 85.468 > 41 50 Rb 0.28542276501060 0.95208943167727 0.58333333333333 85.468 > 42 51 Rb 0.66666666666667 0.71457723498940 0.58333333333333 85.468 > 43 52 Rb 0.04791056832273 0.33333333333333 0.58333333333333 85.468 > 44 53 Rb 0.66666666666667 0.33333333333333 0.58333333333333 85.468 > 37 54 Rb 0.66666666666667 0.33333333333333 0.08333333333333 85.468 > 38 55 Rb 0.38124390165606 0.38124390165606 0.75000000000000 85.468 > 39 56 Rb 0.00000000000000 0.61875609834394 0.75000000000000 85.468 > 40 57 Rb 0.61875609834394 0.00000000000000 0.75000000000000 85.468 > 41 58 Rb 0.95208943167727 0.28542276501060 0.91666666666667 85.468 > 42 59 Rb 0.33333333333333 0.04791056832273 0.91666666666667 85.468 > 43 60 Rb 0.71457723498940 0.66666666666667 0.91666666666667 85.468 > 44 61 Pb 0.00000000000000 0.00000000000000 0.00000000000000 207.200 > 61 62 Pb 0.33333333333333 0.66666666666667 0.16666666666667 207.200 > 62 63 Pb 0.66666666666667 0.33333333333333 0.33333333333333 207.200 > 61 64 Pb 0.00000000000000 0.00000000000000 0.50000000000000 207.200 > 62 65 Pb 0.33333333333333 0.66666666666667 0.66666666666667 207.200 > 61 66 Pb 0.66666666666667 0.33333333333333 0.83333333333333 207.200 > 62 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.1828860 0.0000000 0.0000000 -0.0000000 3.1828860 0.0000000 0.0000000 0.0000000 3.1968694 -------------------------- Born effective charges -------------------------- 1 Br -1.0375216 0.3229362 -0.2279541 0.2388177 -1.5573504 0.5110383 -0.2313857 0.5085382 -1.3620732 2 Br -1.6706398 0.1267132 0.5565492 0.0425947 -0.9242322 -0.0581051 0.5560999 -0.0538832 -1.3620732 3 Br -1.1841466 -0.3234717 -0.3285951 -0.4075902 -1.4107254 -0.4529332 -0.3247142 -0.4546550 -1.3620732 4 Br -1.6706398 -0.1267132 -0.5565492 -0.0425947 -0.9242322 -0.0581051 -0.5560999 -0.0538832 -1.3620732 5 Br -1.0375216 -0.3229362 0.2279541 -0.2388177 -1.5573504 0.5110383 0.2313857 0.5085382 -1.3620732 6 Br -1.1841466 0.3234717 0.3285951 0.4075902 -1.4107254 -0.4529332 0.3247142 -0.4546550 -1.3620732 7 Br -1.0375216 0.3229362 -0.2279541 0.2388177 -1.5573504 0.5110383 -0.2313857 0.5085382 -1.3620732 8 Br -1.6706398 0.1267132 0.5565492 0.0425947 -0.9242322 -0.0581051 0.5560999 -0.0538832 -1.3620732 9 Br -1.1841466 -0.3234717 -0.3285951 -0.4075902 -1.4107254 -0.4529332 -0.3247142 -0.4546550 -1.3620732 10 Br -1.6706398 -0.1267132 -0.5565492 -0.0425947 -0.9242322 -0.0581051 -0.5560999 -0.0538832 -1.3620732 11 Br -1.1841466 0.3234717 0.3285951 0.4075902 -1.4107254 -0.4529332 0.3247142 -0.4546550 -1.3620732 12 Br -1.0375216 -0.3229362 0.2279541 -0.2388177 -1.5573504 0.5110383 0.2313857 0.5085382 -1.3620732 13 Rb 0.9094608 0.0000000 -0.0000000 -0.0000000 0.9094608 0.0000000 0.0000000 -0.0000000 1.6152910 14 Rb 0.9094608 0.0000000 -0.0000000 -0.0000000 0.9094608 0.0000000 0.0000000 -0.0000000 1.6152910 15 Rb 1.3573660 -0.0688557 0.1709532 -0.0688557 1.2778583 -0.0986999 0.1363511 -0.0787224 1.1663259 16 Rb 1.3573660 0.0688557 -0.1709532 0.0688557 1.2778583 -0.0986999 -0.1363511 -0.0787224 1.1663259 17 Rb 1.2381044 -0.0000000 -0.0000000 0.0000000 1.3971199 0.1973997 0.0000000 0.1574447 1.1663259 18 Rb 1.3573660 -0.0688557 0.1709532 -0.0688557 1.2778583 -0.0986999 0.1363511 -0.0787224 1.1663259 19 Rb 1.3573660 0.0688557 -0.1709532 0.0688557 1.2778583 -0.0986999 -0.1363511 -0.0787224 1.1663259 20 Rb 1.2381044 -0.0000000 -0.0000000 0.0000000 1.3971199 0.1973997 0.0000000 0.1574447 1.1663259 21 Pb 2.9223190 0.0862955 -0.0000000 -0.0862955 2.9223190 0.0000000 0.0000000 -0.0000000 3.0581705 22 Pb 2.9223190 -0.0862955 -0.0000000 0.0862955 2.9223190 0.0000000 0.0000000 -0.0000000 3.0581705 ---------------------------------------------------------------------------- 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.0050 0.0087 0.0000] [ 0.0050 -0.0087 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 37, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 39, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 61, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000000 (xxx) -0.00000000 (xxx) -0.00000000 (xxx) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 07:20:18]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 07:20:18]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [1 1 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 6.612214075966644 3.817563576698766 5.482577393333327 b -6.612214075966630 3.817563576698766 5.482577393333327 c 0.000000000000007 -7.635127153397544 5.482577393333327 Atomic positions (fractional): 1 Br 0.57259840926760 0.44021628734293 0.80474043540191 79.904 2 Br 0.44021628734293 0.80474043540191 0.57259840926760 79.904 3 Br 0.80474043540191 0.57259840926760 0.44021628734293 79.904 4 Br 0.69525956459809 0.05978371265707 0.92740159073240 79.904 5 Br 0.05978371265707 0.92740159073240 0.69525956459809 79.904 6 Br 0.92740159073240 0.69525956459809 0.05978371265707 79.904 7 Br 0.42740159073240 0.55978371265707 0.19525956459809 79.904 8 Br 0.55978371265707 0.19525956459809 0.42740159073240 79.904 9 Br 0.19525956459809 0.42740159073240 0.55978371265707 79.904 10 Br 0.30474043540191 0.94021628734293 0.07259840926760 79.904 11 Br 0.07259840926760 0.30474043540191 0.94021628734293 79.904 12 Br 0.94021628734293 0.07259840926760 0.30474043540191 79.904 13 Rb 0.25000000000000 0.25000000000000 0.25000000000000 85.468 14 Rb 0.75000000000000 0.75000000000000 0.75000000000000 85.468 15 Rb 0.13124390165606 0.75000000000000 0.36875609834394 85.468 16 Rb 0.75000000000000 0.36875609834394 0.13124390165606 85.468 17 Rb 0.36875609834394 0.13124390165606 0.75000000000000 85.468 18 Rb 0.86875609834394 0.25000000000000 0.63124390165606 85.468 19 Rb 0.25000000000000 0.63124390165606 0.86875609834394 85.468 20 Rb 0.63124390165606 0.86875609834394 0.25000000000000 85.468 21 Pb 0.00000000000000 0.00000000000000 0.00000000000000 207.200 22 Pb 0.50000000000000 0.50000000000000 0.50000000000000 207.200 -------------------------------- supercell --------------------------------- Lattice vectors: a 13.224428151933274 0.000000000000000 0.000000000000000 b -6.612214075966637 11.452690730096311 0.000000000000000 c 0.000000000000000 0.000000000000000 16.447732179999999 Atomic positions (fractional): 1 Br 0.30008003193012 0.46777794193557 0.27251837733748 79.904 > 1 2 Br 0.16769791000545 0.69991996806988 0.27251837733748 79.904 > 2 3 Br 0.53222205806443 0.83230208999455 0.27251837733748 79.904 > 3 4 Br 0.46777794193557 0.30008003193012 0.22748162266252 79.904 > 4 5 Br 0.83230208999455 0.53222205806443 0.22748162266252 79.904 > 5 6 Br 0.69991996806988 0.16769791000545 0.22748162266252 79.904 > 6 7 Br 0.36658663473655 0.86555539139777 0.06081495599585 79.904 > 7 8 Br 0.49896875666122 0.63341336526345 0.06081495599585 79.904 > 8 9 Br 0.13444460860224 0.50103124333878 0.06081495599585 79.904 > 9 10 Br 0.19888872473110 0.03325330140322 0.10585171067081 79.904 > 10 11 Br 0.96674669859678 0.16563542332788 0.10585171067081 79.904 > 11 12 Br 0.83436457667212 0.80111127526890 0.10585171067081 79.904 > 12 13 Br 0.96674669859678 0.80111127526890 0.60585171067081 79.904 > 1 14 Br 0.83436457667212 0.03325330140322 0.60585171067081 79.904 > 2 15 Br 0.19888872473110 0.16563542332788 0.60585171067081 79.904 > 3 16 Br 0.13444460860224 0.63341336526345 0.56081495599585 79.904 > 4 17 Br 0.49896875666122 0.86555539139777 0.56081495599585 79.904 > 5 18 Br 0.36658663473655 0.50103124333878 0.56081495599585 79.904 > 6 19 Br 0.03325330140322 0.19888872473110 0.39414828932919 79.904 > 7 20 Br 0.16563542332788 0.96674669859678 0.39414828932919 79.904 > 8 21 Br 0.80111127526890 0.83436457667212 0.39414828932919 79.904 > 9 22 Br 0.86555539139777 0.36658663473655 0.43918504400414 79.904 > 10 23 Br 0.63341336526345 0.49896875666122 0.43918504400414 79.904 > 11 24 Br 0.50103124333878 0.13444460860223 0.43918504400414 79.904 > 12 25 Br 0.63341336526345 0.13444460860224 0.93918504400415 79.904 > 1 26 Br 0.50103124333878 0.36658663473655 0.93918504400415 79.904 > 2 27 Br 0.86555539139777 0.49896875666122 0.93918504400415 79.904 > 3 28 Br 0.80111127526890 0.96674669859678 0.89414828932919 79.904 > 4 29 Br 0.16563542332788 0.19888872473110 0.89414828932919 79.904 > 5 30 Br 0.03325330140322 0.83436457667212 0.89414828932919 79.904 > 6 31 Br 0.69991996806988 0.53222205806443 0.72748162266252 79.904 > 7 32 Br 0.83230208999455 0.30008003193012 0.72748162266252 79.904 > 8 33 Br 0.46777794193557 0.16769791000545 0.72748162266252 79.904 > 9 34 Br 0.53222205806443 0.69991996806988 0.77251837733748 79.904 > 10 35 Br 0.30008003193012 0.83230208999455 0.77251837733748 79.904 > 11 36 Br 0.16769791000545 0.46777794193557 0.77251837733748 79.904 > 12 37 Rb 0.33333333333333 0.66666666666667 0.91666666666667 85.468 > 37 38 Rb 0.33333333333333 0.66666666666667 0.41666666666667 85.468 > 38 39 Rb 0.04791056832273 0.71457723498940 0.08333333333333 85.468 > 39 40 Rb 0.66666666666667 0.95208943167727 0.08333333333333 85.468 > 40 41 Rb 0.28542276501060 0.33333333333333 0.08333333333333 85.468 > 41 42 Rb 0.61875609834394 0.61875609834394 0.25000000000000 85.468 > 42 43 Rb 0.00000000000000 0.38124390165606 0.25000000000000 85.468 > 43 44 Rb 0.38124390165606 0.00000000000000 0.25000000000000 85.468 > 44 45 Rb 0.00000000000000 0.00000000000000 0.25000000000000 85.468 > 37 46 Rb 0.00000000000000 0.00000000000000 0.75000000000000 85.468 > 38 47 Rb 0.71457723498940 0.04791056832273 0.41666666666667 85.468 > 39 48 Rb 0.33333333333333 0.28542276501060 0.41666666666667 85.468 > 40 49 Rb 0.95208943167727 0.66666666666667 0.41666666666667 85.468 > 41 50 Rb 0.28542276501060 0.95208943167727 0.58333333333333 85.468 > 42 51 Rb 0.66666666666667 0.71457723498940 0.58333333333333 85.468 > 43 52 Rb 0.04791056832273 0.33333333333333 0.58333333333333 85.468 > 44 53 Rb 0.66666666666667 0.33333333333333 0.58333333333333 85.468 > 37 54 Rb 0.66666666666667 0.33333333333333 0.08333333333333 85.468 > 38 55 Rb 0.38124390165606 0.38124390165606 0.75000000000000 85.468 > 39 56 Rb 0.00000000000000 0.61875609834394 0.75000000000000 85.468 > 40 57 Rb 0.61875609834394 0.00000000000000 0.75000000000000 85.468 > 41 58 Rb 0.95208943167727 0.28542276501060 0.91666666666667 85.468 > 42 59 Rb 0.33333333333333 0.04791056832273 0.91666666666667 85.468 > 43 60 Rb 0.71457723498940 0.66666666666667 0.91666666666667 85.468 > 44 61 Pb 0.00000000000000 0.00000000000000 0.00000000000000 207.200 > 61 62 Pb 0.33333333333333 0.66666666666667 0.16666666666667 207.200 > 62 63 Pb 0.66666666666667 0.33333333333333 0.33333333333333 207.200 > 61 64 Pb 0.00000000000000 0.00000000000000 0.50000000000000 207.200 > 62 65 Pb 0.33333333333333 0.66666666666667 0.66666666666667 207.200 > 61 66 Pb 0.66666666666667 0.33333333333333 0.83333333333333 207.200 > 62 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 3.1828860 0.0000000 0.0000000 -0.0000000 3.1828860 0.0000000 0.0000000 0.0000000 3.1968694 -------------------------- Born effective charges -------------------------- 1 Br -1.0375216 0.3229362 -0.2279541 0.2388177 -1.5573504 0.5110383 -0.2313857 0.5085382 -1.3620732 2 Br -1.6706398 0.1267132 0.5565492 0.0425947 -0.9242322 -0.0581051 0.5560999 -0.0538832 -1.3620732 3 Br -1.1841466 -0.3234717 -0.3285951 -0.4075902 -1.4107254 -0.4529332 -0.3247142 -0.4546550 -1.3620732 4 Br -1.6706398 -0.1267132 -0.5565492 -0.0425947 -0.9242322 -0.0581051 -0.5560999 -0.0538832 -1.3620732 5 Br -1.0375216 -0.3229362 0.2279541 -0.2388177 -1.5573504 0.5110383 0.2313857 0.5085382 -1.3620732 6 Br -1.1841466 0.3234717 0.3285951 0.4075902 -1.4107254 -0.4529332 0.3247142 -0.4546550 -1.3620732 7 Br -1.0375216 0.3229362 -0.2279541 0.2388177 -1.5573504 0.5110383 -0.2313857 0.5085382 -1.3620732 8 Br -1.6706398 0.1267132 0.5565492 0.0425947 -0.9242322 -0.0581051 0.5560999 -0.0538832 -1.3620732 9 Br -1.1841466 -0.3234717 -0.3285951 -0.4075902 -1.4107254 -0.4529332 -0.3247142 -0.4546550 -1.3620732 10 Br -1.6706398 -0.1267132 -0.5565492 -0.0425947 -0.9242322 -0.0581051 -0.5560999 -0.0538832 -1.3620732 11 Br -1.1841466 0.3234717 0.3285951 0.4075902 -1.4107254 -0.4529332 0.3247142 -0.4546550 -1.3620732 12 Br -1.0375216 -0.3229362 0.2279541 -0.2388177 -1.5573504 0.5110383 0.2313857 0.5085382 -1.3620732 13 Rb 0.9094608 0.0000000 -0.0000000 -0.0000000 0.9094608 0.0000000 0.0000000 -0.0000000 1.6152910 14 Rb 0.9094608 0.0000000 -0.0000000 -0.0000000 0.9094608 0.0000000 0.0000000 -0.0000000 1.6152910 15 Rb 1.3573660 -0.0688557 0.1709532 -0.0688557 1.2778583 -0.0986999 0.1363511 -0.0787224 1.1663259 16 Rb 1.3573660 0.0688557 -0.1709532 0.0688557 1.2778583 -0.0986999 -0.1363511 -0.0787224 1.1663259 17 Rb 1.2381044 -0.0000000 -0.0000000 0.0000000 1.3971199 0.1973997 0.0000000 0.1574447 1.1663259 18 Rb 1.3573660 -0.0688557 0.1709532 -0.0688557 1.2778583 -0.0986999 0.1363511 -0.0787224 1.1663259 19 Rb 1.3573660 0.0688557 -0.1709532 0.0688557 1.2778583 -0.0986999 -0.1363511 -0.0787224 1.1663259 20 Rb 1.2381044 -0.0000000 -0.0000000 0.0000000 1.3971199 0.1973997 0.0000000 0.1574447 1.1663259 21 Pb 2.9223190 0.0862955 -0.0000000 -0.0862955 2.9223190 0.0000000 0.0000000 -0.0000000 3.0581705 22 Pb 2.9223190 -0.0862955 -0.0000000 0.0862955 2.9223190 0.0000000 0.0000000 -0.0000000 3.0581705 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (xxx) -0.00000000 (xxx) -0.00000000 (xxx) Max drift of fc2: 0.00000000 (yy) 0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 5 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: 305, Lambda: 0.08 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/19) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 19 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.702 ( 0.000 0.000 0.000) 0.000 0.702 ( 0.000 0.000 0.000) 0.000 1.010 ( 0.000 0.000 0.000) 0.000 1.010 ( 0.000 0.000 0.000) 0.000 1.118 ( 0.000 0.000 0.000) 0.000 1.118 ( 0.000 0.000 0.000) 0.000 1.232 ( 0.000 0.000 0.000) 0.000 1.284 ( 0.000 0.000 0.000) 0.000 1.321 ( 0.000 0.000 0.000) 0.000 1.321 ( 0.000 0.000 0.000) 0.000 1.381 ( 0.000 0.000 0.000) 0.000 1.381 ( 0.000 0.000 0.000) 0.000 1.466 ( 0.000 0.000 0.000) 0.000 1.483 ( 0.000 0.000 0.000) 0.000 1.521 ( 0.000 0.000 0.000) 0.000 1.521 ( 0.000 0.000 0.000) 0.000 1.544 ( 0.000 0.000 0.000) 0.000 1.544 ( 0.000 0.000 0.000) 0.000 1.573 ( 0.000 0.000 0.000) 0.000 1.573 ( 0.000 0.000 0.000) 0.000 1.596 ( 0.000 0.000 0.000) 0.000 1.645 ( 0.000 0.000 0.000) 0.000 1.650 ( 0.000 0.000 0.000) 0.000 1.717 ( 0.000 0.000 0.000) 0.000 1.756 ( 0.000 0.000 0.000) 0.000 1.756 ( 0.000 0.000 0.000) 0.000 1.805 ( 0.000 0.000 0.000) 0.000 1.805 ( 0.000 0.000 0.000) 0.000 1.844 ( 0.000 0.000 0.000) 0.000 1.844 ( 0.000 0.000 0.000) 0.000 1.874 ( 0.000 0.000 0.000) 0.000 1.968 ( 0.000 0.000 0.000) 0.000 1.968 ( 0.000 0.000 0.000) 0.000 1.987 ( 0.000 0.000 0.000) 0.000 2.010 ( 0.000 0.000 0.000) 0.000 2.010 ( 0.000 0.000 0.000) 0.000 2.127 ( 0.000 0.000 0.000) 0.000 2.127 ( 0.000 0.000 0.000) 0.000 2.163 ( 0.000 0.000 0.000) 0.000 2.163 ( 0.000 0.000 0.000) 0.000 2.230 ( 0.000 0.000 0.000) 0.000 2.230 ( 0.000 0.000 0.000) 0.000 2.363 ( 0.000 0.000 0.000) 0.000 2.501 ( 0.000 0.000 0.000) 0.000 2.532 ( 0.000 0.000 0.000) 0.000 2.532 ( 0.000 0.000 0.000) 0.000 2.544 ( 0.000 0.000 0.000) 0.000 2.544 ( 0.000 0.000 0.000) 0.000 2.556 ( 0.000 0.000 0.000) 0.000 2.721 ( 0.000 0.000 0.000) 0.000 2.945 ( 0.000 0.000 0.000) 0.000 2.945 ( 0.000 0.000 0.000) 0.000 3.141 ( 0.000 0.000 0.000) 0.000 3.219 ( 0.000 0.000 0.000) 0.000 3.291 ( 0.000 0.000 0.000) 0.000 3.291 ( 0.000 0.000 0.000) 0.000 3.305 ( 0.000 0.000 0.000) 0.000 3.401 ( 0.000 0.000 0.000) 0.000 3.401 ( 0.000 0.000 0.000) 0.000 3.486 ( 0.000 0.000 0.000) 0.000 3.486 ( 0.000 0.000 0.000) 0.000 3.801 ( 0.000 0.000 0.000) 0.000 3.856 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/19) ======================= q-point: ( 0.20 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.224 ( 7.641 4.412 5.861) 10.592 0.253 ( 6.432 3.714 9.535) 12.086 0.506 ( 13.559 7.828 16.139) 22.486 0.713 ( 0.685 0.395 -0.442) 0.906 0.737 ( 0.922 0.532 3.023) 3.205 0.968 ( -3.223 -1.861 -0.755) 3.798 1.038 ( -1.111 -0.641 3.521) 3.747 1.078 ( -0.230 -0.133 -2.804) 2.816 1.133 ( 3.438 1.985 -2.008) 4.449 1.229 ( -2.143 -1.237 1.922) 3.133 1.246 ( -0.882 -0.509 -3.749) 3.885 1.278 ( -0.907 -0.524 -5.393) 5.494 1.327 ( 0.751 0.434 -1.390) 1.638 1.353 ( -1.382 -0.798 -0.298) 1.624 1.408 ( 1.102 0.636 1.024) 1.633 1.463 ( -4.549 -2.626 -1.476) 5.456 1.466 ( -2.430 -1.403 1.465) 3.166 1.479 ( -1.338 -0.772 1.283) 2.008 1.531 ( 0.288 0.166 4.458) 4.471 1.556 ( -2.189 -1.264 0.715) 2.627 1.567 ( 0.423 0.244 0.581) 0.759 1.581 ( 1.781 1.028 0.242) 2.070 1.608 ( 2.953 1.705 0.380) 3.431 1.618 ( 1.639 0.947 3.137) 3.664 1.643 ( -1.763 -1.018 0.382) 2.072 1.666 ( 2.369 1.368 -0.796) 2.849 1.730 ( 0.653 0.377 0.624) 0.978 1.753 ( 0.895 0.517 -1.691) 1.982 1.784 ( 2.507 1.448 0.126) 2.898 1.789 ( -3.104 -1.792 -0.879) 3.690 1.830 ( -0.357 -0.206 -0.988) 1.070 1.839 ( -0.219 -0.126 -0.036) 0.255 1.850 ( -0.283 -0.163 2.343) 2.366 1.867 ( -0.376 -0.217 0.378) 0.575 1.947 ( -0.667 -0.385 -0.097) 0.776 1.967 ( -0.915 -0.528 -1.554) 1.880 1.977 ( -0.313 -0.181 -0.382) 0.526 1.992 ( -0.468 -0.270 -1.416) 1.516 2.055 ( -0.379 -0.219 -2.076) 2.122 2.109 ( -0.404 -0.234 -5.465) 5.485 2.130 ( 3.875 2.237 -6.845) 8.178 2.182 ( -1.974 -1.139 -2.017) 3.043 2.247 ( 4.538 2.620 -5.718) 7.756 2.255 ( 0.903 0.521 6.062) 6.151 2.273 ( 0.646 0.373 4.795) 4.852 2.328 ( 0.076 0.044 -4.001) 4.002 2.493 ( -0.144 -0.083 -0.215) 0.272 2.509 ( -0.309 -0.178 -2.184) 2.213 2.525 ( -0.013 -0.007 -2.642) 2.642 2.549 ( 1.190 0.687 -1.469) 2.012 2.550 ( -0.328 -0.189 0.091) 0.389 2.571 ( -1.617 -0.934 3.508) 3.974 2.870 ( -4.447 -2.567 11.152) 12.277 2.976 ( 2.969 1.714 0.485) 3.463 3.064 ( 0.074 0.043 -7.918) 7.919 3.225 ( -2.421 -1.398 -0.911) 2.940 3.276 ( -0.271 -0.157 4.689) 4.700 3.284 ( -4.895 -2.826 6.824) 8.861 3.292 ( 1.116 0.644 -1.771) 2.190 3.342 ( 2.220 1.281 -2.194) 3.374 3.383 ( 0.599 0.346 0.996) 1.212 3.490 ( 0.262 0.151 -0.193) 0.359 3.617 ( 5.074 2.929 6.933) 9.077 3.806 ( 0.368 0.212 0.361) 0.557 3.851 ( -0.258 -0.149 -0.521) 0.601 4.056 ( -2.008 -1.159 -1.207) 2.614 ======================= Grid point 2 (3/19) ======================= q-point: ( 0.40 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.445 ( 7.283 4.205 4.800) 9.682 0.481 ( 4.653 2.686 8.240) 9.837 0.662 ( -3.964 -2.288 -2.715) 5.322 0.703 ( -3.700 -2.136 -4.771) 6.405 0.892 ( 7.761 4.481 7.813) 11.889 0.985 ( 1.825 1.054 9.195) 9.434 1.046 ( -0.226 -0.131 -2.729) 2.742 1.071 ( 0.088 0.051 2.948) 2.950 1.144 ( -1.311 -0.757 -7.901) 8.045 1.177 ( -3.284 -1.896 -1.044) 3.933 1.194 ( 3.121 1.802 -1.502) 3.905 1.214 ( 0.320 0.185 -2.831) 2.855 1.317 ( -2.743 -1.584 1.460) 3.488 1.329 ( -1.296 -0.748 -0.048) 1.498 1.393 ( 2.097 1.211 0.521) 2.477 1.428 ( -1.631 -0.942 0.794) 2.044 1.438 ( 0.926 0.535 0.457) 1.163 1.488 ( -0.004 -0.002 2.108) 2.108 1.527 ( -1.468 -0.847 -0.501) 1.767 1.546 ( -4.078 -2.354 -2.319) 5.249 1.591 ( 1.762 1.017 0.337) 2.062 1.603 ( -0.025 -0.014 2.976) 2.976 1.623 ( 0.134 0.077 0.451) 0.477 1.672 ( 0.459 0.265 -4.334) 4.366 1.688 ( -0.659 -0.380 0.201) 0.787 1.696 ( 1.891 1.092 3.073) 3.770 1.717 ( -0.772 -0.446 1.409) 1.667 1.729 ( -0.536 -0.309 -1.866) 1.966 1.748 ( 1.200 0.693 0.637) 1.525 1.775 ( -1.298 -0.749 -1.208) 1.925 1.826 ( 0.838 0.484 0.330) 1.023 1.831 ( -0.321 -0.185 0.623) 0.725 1.851 ( -0.684 -0.395 -0.880) 1.182 1.860 ( 0.544 0.314 -0.339) 0.714 1.926 ( -1.730 -0.999 1.280) 2.372 1.937 ( 0.935 0.540 -3.122) 3.304 1.989 ( -0.017 -0.010 0.894) 0.895 2.012 ( 1.533 0.885 0.844) 1.961 2.029 ( -0.100 -0.058 -4.137) 4.138 2.047 ( -0.909 -0.525 1.019) 1.463 2.119 ( 2.836 1.637 -6.077) 6.903 2.122 ( -2.777 -1.603 1.697) 3.628 2.269 ( 1.226 0.708 -0.342) 1.456 2.290 ( -0.067 -0.038 -1.913) 1.915 2.347 ( 1.373 0.793 3.415) 3.765 2.348 ( 0.602 0.347 4.151) 4.209 2.472 ( -0.404 -0.233 -2.067) 2.119 2.482 ( -0.508 -0.293 -2.988) 3.045 2.501 ( 0.419 0.242 1.109) 1.210 2.530 ( -0.495 -0.286 -1.629) 1.726 2.551 ( 1.806 1.042 -2.839) 3.522 2.562 ( -3.023 -1.745 3.447) 4.906 2.920 ( -2.060 -1.189 7.658) 8.019 2.979 ( 0.997 0.576 -7.376) 7.465 3.069 ( 4.790 2.766 0.487) 5.553 3.161 ( -3.779 -2.182 0.093) 4.365 3.274 ( -2.561 -1.479 4.301) 5.220 3.282 ( 1.668 0.963 -3.609) 4.091 3.319 ( -0.491 -0.284 3.421) 3.468 3.344 ( 1.042 0.602 -2.934) 3.172 3.435 ( 0.973 0.562 3.255) 3.443 3.480 ( -0.263 -0.152 -1.768) 1.794 3.788 ( 2.907 1.678 3.826) 5.090 3.832 ( -0.624 -0.360 -1.021) 1.250 3.841 ( 2.774 1.602 3.163) 4.502 3.970 ( -3.520 -2.032 -3.798) 5.562 ======================= Grid point 6 (4/19) ======================= q-point: ( 0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.339 ( -0.000 2.003 11.691) 11.861 0.403 ( -0.000 5.054 11.444) 12.510 0.701 ( -0.000 4.034 13.015) 13.626 0.746 ( -0.000 6.143 10.754) 12.385 0.793 ( -0.000 2.982 -0.213) 2.990 0.986 ( -0.000 -4.297 1.584) 4.580 1.022 ( 0.000 -0.191 -1.801) 1.811 1.073 ( 0.000 1.681 -2.812) 3.276 1.122 ( -0.000 2.168 4.356) 4.865 1.187 ( 0.000 -3.622 -7.208) 8.067 1.187 ( 0.000 -0.625 -5.181) 5.219 1.275 ( -0.000 -1.618 2.491) 2.970 1.284 ( 0.000 0.573 -2.965) 3.020 1.365 ( 0.000 0.025 -1.042) 1.042 1.405 ( -0.000 0.218 2.170) 2.181 1.474 ( 0.000 -1.737 -1.062) 2.036 1.479 ( 0.000 -2.638 -0.388) 2.666 1.523 ( -0.000 -1.424 1.878) 2.357 1.556 ( -0.000 -1.935 2.136) 2.882 1.561 ( 0.000 0.166 -0.109) 0.199 1.573 ( -0.000 0.264 1.099) 1.130 1.574 ( 0.000 1.172 -0.596) 1.315 1.610 ( -0.000 -0.484 0.268) 0.553 1.628 ( -0.000 -0.852 0.652) 1.073 1.646 ( -0.000 2.397 0.811) 2.530 1.650 ( -0.000 -2.170 1.350) 2.556 1.712 ( 0.000 -1.456 -0.788) 1.655 1.762 ( -0.000 4.170 1.044) 4.298 1.804 ( -0.000 1.838 -0.348) 1.871 1.816 ( 0.000 -2.070 -2.094) 2.945 1.818 ( 0.000 -0.867 0.180) 0.885 1.860 ( -0.000 3.326 2.869) 4.392 1.865 ( 0.000 -1.347 0.336) 1.388 1.878 ( -0.000 -0.043 1.741) 1.742 1.916 ( 0.000 0.275 -3.027) 3.039 1.937 ( 0.000 -1.617 -2.304) 2.814 1.972 ( 0.000 0.454 -3.166) 3.198 2.008 ( 0.000 -0.762 -0.032) 0.762 2.012 ( 0.000 2.121 -4.650) 5.111 2.012 ( -0.000 -0.298 0.408) 0.506 2.073 ( 0.000 1.820 -5.058) 5.375 2.166 ( 0.000 -1.110 0.091) 1.114 2.208 ( -0.000 0.930 1.173) 1.497 2.275 ( 0.000 1.214 -3.736) 3.928 2.316 ( -0.000 -0.171 4.857) 4.860 2.349 ( -0.000 3.491 3.049) 4.635 2.476 ( 0.000 -0.524 -3.404) 3.444 2.489 ( 0.000 -0.344 -2.774) 2.795 2.489 ( 0.000 0.831 -0.466) 0.953 2.529 ( 0.000 -0.399 -1.048) 1.121 2.541 ( 0.000 0.309 -2.083) 2.106 2.624 ( -0.000 -2.124 3.807) 4.360 2.936 ( -0.000 -2.222 3.864) 4.457 2.960 ( 0.000 1.081 -8.967) 9.032 3.022 ( -0.000 4.400 2.031) 4.846 3.192 ( 0.000 -4.076 0.297) 4.087 3.260 ( 0.000 -3.654 -1.480) 3.943 3.284 ( 0.000 1.752 -1.563) 2.348 3.311 ( -0.000 -3.307 3.313) 4.681 3.401 ( -0.000 4.255 1.741) 4.597 3.431 ( -0.000 4.109 1.991) 4.566 3.535 ( -0.000 -2.824 3.270) 4.321 3.677 ( -0.000 3.142 7.895) 8.497 3.810 ( -0.000 0.314 0.573) 0.653 3.845 ( 0.000 0.082 -0.803) 0.807 4.030 ( 0.000 -3.427 -2.319) 4.138 ======================= Grid point 7 (5/19) ======================= q-point: ( 0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 63 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.504 ( 1.874 2.356 8.969) 9.460 0.592 ( 2.471 4.255 9.722) 10.896 0.715 ( -5.938 -1.652 -3.930) 7.310 0.728 ( -8.055 -1.604 2.424) 8.563 0.956 ( 2.446 0.710 4.345) 5.036 0.995 ( 0.941 -0.222 -1.104) 1.468 1.054 ( -0.918 1.410 -2.840) 3.300 1.088 ( 2.335 2.683 3.074) 4.701 1.111 ( 3.559 0.881 0.280) 3.677 1.144 ( 1.874 0.166 -0.943) 2.105 1.171 ( -2.787 0.855 3.981) 4.935 1.217 ( -4.076 0.132 -1.645) 4.397 1.303 ( 1.582 -2.651 0.593) 3.143 1.376 ( 0.044 -1.161 -0.153) 1.172 1.396 ( -2.107 -0.564 2.214) 3.108 1.443 ( 0.185 -0.369 2.208) 2.246 1.460 ( -1.303 -0.094 0.891) 1.582 1.493 ( -0.727 -1.240 -1.898) 2.381 1.527 ( -1.441 0.526 -1.743) 2.322 1.540 ( -1.635 -2.191 -1.985) 3.378 1.570 ( -0.043 -0.802 -1.729) 1.907 1.589 ( 1.640 0.635 -2.336) 2.924 1.613 ( 1.694 -1.683 -1.037) 2.604 1.639 ( 1.120 -0.590 -0.581) 1.393 1.660 ( 0.628 0.387 1.465) 1.640 1.689 ( 0.412 -2.837 -0.436) 2.900 1.709 ( -1.335 -0.031 -0.924) 1.624 1.728 ( 3.519 -0.396 2.028) 4.081 1.796 ( -2.162 0.903 1.321) 2.690 1.798 ( -2.705 1.985 -1.997) 3.904 1.824 ( -0.749 0.138 0.137) 0.774 1.843 ( -0.462 -0.300 -0.416) 0.690 1.869 ( -0.850 0.501 0.167) 1.001 1.889 ( -1.742 1.420 -0.233) 2.259 1.909 ( 1.082 0.915 2.488) 2.864 1.935 ( 2.099 1.099 -2.847) 3.704 1.979 ( 0.018 -0.615 -1.871) 1.969 1.998 ( 1.243 0.996 2.486) 2.952 2.006 ( 1.066 0.593 0.217) 1.239 2.024 ( 1.647 2.076 -3.021) 4.019 2.065 ( 0.433 -0.126 1.160) 1.245 2.153 ( -0.086 -1.663 1.321) 2.125 2.251 ( 1.308 0.459 2.211) 2.610 2.264 ( 0.681 0.397 -1.263) 1.489 2.360 ( -0.178 0.665 -1.837) 1.962 2.369 ( 0.079 -0.017 -1.901) 1.903 2.451 ( 0.896 -0.153 0.457) 1.017 2.462 ( 0.223 -0.368 0.448) 0.621 2.504 ( 0.563 1.757 -0.814) 2.017 2.510 ( -0.133 0.124 -1.729) 1.738 2.516 ( 0.350 0.107 -0.886) 0.959 2.623 ( -1.931 -2.288 2.444) 3.865 2.873 ( 0.756 1.238 -8.019) 8.150 2.981 ( 1.373 -0.772 3.684) 4.007 3.079 ( -0.117 4.271 1.482) 4.522 3.176 ( -0.413 -2.078 -0.023) 2.118 3.209 ( 0.882 -2.312 -1.858) 3.095 3.290 ( 0.786 1.092 -1.803) 2.250 3.322 ( -0.434 -2.457 3.069) 3.955 3.426 ( -4.047 1.962 3.418) 5.649 3.494 ( -0.144 3.461 2.595) 4.328 3.505 ( -2.908 -1.027 0.399) 3.110 3.804 ( 1.463 0.415 2.973) 3.340 3.823 ( -0.190 -0.333 -0.368) 0.531 3.854 ( 2.663 1.742 2.472) 4.029 3.946 ( -2.339 -3.226 -3.275) 5.158 ======================= Grid point 8 (6/19) ======================= q-point: (-0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 63 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.508 ( -8.854 -0.794 -3.315) 9.488 0.558 ( -10.801 -1.913 0.383) 10.976 0.676 ( 0.697 0.566 5.566) 5.638 0.772 ( 1.739 4.640 5.603) 7.479 0.879 ( 0.624 -3.239 -4.271) 5.396 0.934 ( -9.840 -2.293 -5.515) 11.511 1.041 ( -1.820 -0.557 -3.673) 4.137 1.118 ( -2.594 2.463 1.832) 4.019 1.151 ( 1.145 2.257 -0.181) 2.538 1.178 ( 3.340 -1.126 -0.460) 3.555 1.204 ( 2.099 -0.341 3.610) 4.190 1.250 ( 0.767 1.568 6.425) 6.658 1.304 ( 0.074 -0.795 -0.239) 0.834 1.334 ( -0.790 -0.416 0.302) 0.943 1.392 ( 0.241 -1.218 -0.048) 1.243 1.429 ( -2.470 -0.236 1.181) 2.748 1.448 ( -0.212 -0.914 -1.972) 2.184 1.465 ( -0.423 -0.589 -0.837) 1.107 1.489 ( 0.315 0.348 -3.205) 3.239 1.516 ( -0.269 -1.143 -1.639) 2.016 1.549 ( 2.859 -0.831 0.969) 3.131 1.583 ( 0.012 0.743 -1.775) 1.924 1.638 ( 1.972 -0.422 -2.247) 3.020 1.647 ( 0.631 1.932 -0.031) 2.033 1.680 ( -0.826 1.045 -3.117) 3.389 1.689 ( 0.474 -1.498 1.996) 2.540 1.715 ( 1.979 0.151 2.321) 3.054 1.743 ( -0.612 -0.393 -0.860) 1.126 1.771 ( 0.847 -1.732 2.861) 3.450 1.798 ( -1.579 1.420 0.592) 2.205 1.810 ( 2.074 -0.145 0.001) 2.079 1.840 ( -2.176 0.382 -0.794) 2.348 1.862 ( -0.912 0.380 0.897) 1.334 1.894 ( -0.811 2.303 0.100) 2.444 1.929 ( 0.542 0.480 0.052) 0.726 1.959 ( 1.519 0.058 0.503) 1.601 1.966 ( -0.585 -0.267 -0.638) 0.906 2.014 ( -0.152 -0.389 -2.401) 2.437 2.032 ( -0.351 0.214 4.574) 4.593 2.067 ( 3.599 1.337 -1.163) 4.011 2.099 ( 0.073 2.188 4.686) 5.172 2.193 ( 1.092 -2.117 5.599) 6.085 2.240 ( -0.157 -1.483 -3.938) 4.211 2.282 ( 0.099 -1.904 -2.392) 3.059 2.307 ( 0.884 0.589 -0.696) 1.270 2.328 ( -0.511 0.738 -3.630) 3.739 2.482 ( 1.816 0.119 -1.559) 2.396 2.488 ( 0.646 0.289 0.053) 0.709 2.499 ( 0.612 0.565 1.893) 2.068 2.513 ( 0.672 0.396 0.393) 0.874 2.543 ( 0.652 0.897 0.911) 1.435 2.608 ( -1.726 -1.097 2.729) 3.410 2.802 ( 0.683 0.591 -7.279) 7.335 3.045 ( 1.301 0.273 3.500) 3.744 3.119 ( -2.273 -0.482 0.157) 2.329 3.135 ( 0.144 1.659 0.346) 1.701 3.189 ( 1.469 -1.364 -2.981) 3.592 3.278 ( -0.246 0.408 -1.944) 2.002 3.327 ( -0.553 -1.305 2.014) 2.463 3.401 ( -3.647 -2.616 2.482) 5.129 3.462 ( -2.936 2.753 -0.882) 4.121 3.545 ( -1.679 2.923 3.162) 4.621 3.768 ( -3.703 -3.177 -3.868) 6.226 3.820 ( -0.103 -0.874 -0.424) 0.977 3.851 ( -0.891 0.210 -0.468) 1.028 3.969 ( 3.324 1.000 2.593) 4.332 ======================= Grid point 9 (7/19) ======================= q-point: (-0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.325 ( -10.522 -0.000 0.000) 10.522 0.390 ( -11.749 -0.000 0.000) 11.749 0.618 ( -18.781 -0.000 0.000) 18.781 0.701 ( -0.819 -0.000 0.000) 0.819 0.792 ( -4.008 -0.000 0.000) 4.008 0.939 ( 4.757 0.000 -0.000) 4.757 0.985 ( 1.039 0.000 -0.000) 1.039 1.104 ( -1.394 -0.000 0.000) 1.394 1.188 ( -3.551 -0.000 0.000) 3.551 1.188 ( 2.110 0.000 -0.000) 2.110 1.236 ( 3.216 0.000 -0.000) 3.216 1.311 ( 1.055 0.000 -0.000) 1.055 1.323 ( 1.233 0.000 -0.000) 1.233 1.352 ( -0.113 -0.000 0.000) 0.113 1.395 ( 1.206 0.000 -0.000) 1.206 1.437 ( 1.485 0.000 -0.000) 1.485 1.447 ( 0.937 0.000 -0.000) 0.937 1.471 ( 0.647 0.000 -0.000) 0.647 1.473 ( 0.404 0.000 -0.000) 0.404 1.501 ( 4.412 0.000 -0.000) 4.412 1.564 ( -0.870 -0.000 0.000) 0.870 1.580 ( 1.312 0.000 -0.000) 1.312 1.603 ( -0.491 -0.000 0.000) 0.491 1.633 ( -3.686 -0.000 0.000) 3.686 1.650 ( -0.340 -0.000 0.000) 0.340 1.687 ( -1.043 -0.000 0.000) 1.043 1.754 ( -0.429 -0.000 0.000) 0.429 1.775 ( 0.558 0.000 -0.000) 0.558 1.775 ( 0.032 0.000 -0.000) 0.032 1.817 ( 0.362 0.000 -0.000) 0.362 1.821 ( -0.401 -0.000 0.000) 0.401 1.829 ( 1.994 0.000 -0.000) 1.994 1.860 ( -1.364 -0.000 0.000) 1.364 1.869 ( -2.330 -0.000 0.000) 2.330 1.938 ( -0.751 -0.000 0.000) 0.751 1.955 ( 0.536 0.000 -0.000) 0.536 1.970 ( 1.395 0.000 -0.000) 1.395 1.999 ( 0.485 0.000 -0.000) 0.485 2.098 ( 1.302 0.000 -0.000) 1.302 2.118 ( 0.772 0.000 -0.000) 0.772 2.142 ( 2.319 0.000 -0.000) 2.319 2.200 ( 0.671 0.000 -0.000) 0.671 2.221 ( -1.499 -0.000 0.000) 1.499 2.229 ( -1.389 -0.000 0.000) 1.389 2.315 ( 2.676 0.000 -0.000) 2.676 2.347 ( 1.171 0.000 -0.000) 1.171 2.486 ( 0.721 0.000 -0.000) 0.721 2.511 ( 2.075 0.000 -0.000) 2.075 2.524 ( 0.542 0.000 -0.000) 0.542 2.533 ( 0.833 0.000 -0.000) 0.833 2.555 ( -0.230 -0.000 0.000) 0.230 2.611 ( -1.550 -0.000 0.000) 1.550 2.735 ( -0.875 -0.000 0.000) 0.875 3.035 ( -4.772 -0.000 0.000) 4.772 3.115 ( 1.066 0.000 -0.000) 1.066 3.178 ( 2.341 0.000 -0.000) 2.341 3.184 ( 3.195 0.000 -0.000) 3.195 3.248 ( -1.361 -0.000 0.000) 1.361 3.336 ( 0.693 0.000 -0.000) 0.693 3.351 ( -2.759 -0.000 0.000) 2.759 3.429 ( -2.860 -0.000 0.000) 2.860 3.561 ( -2.824 -0.000 0.000) 2.824 3.617 ( -6.648 -0.000 0.000) 6.648 3.808 ( -0.393 -0.000 0.000) 0.393 3.852 ( 0.132 0.000 -0.000) 0.132 4.034 ( 2.628 0.000 -0.000) 2.628 ======================= Grid point 12 (8/19) ======================= q-point: ( 0.40 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.616 ( -0.000 1.343 8.520) 8.625 0.685 ( -0.000 1.569 4.037) 4.331 0.713 ( 0.000 -0.869 -2.012) 2.192 0.865 ( 0.000 -1.248 -2.028) 2.381 0.868 ( 0.000 -3.114 -5.387) 6.222 0.951 ( -0.000 1.448 5.426) 5.616 1.048 ( -0.000 -0.105 0.346) 0.361 1.051 ( 0.000 -0.751 -3.877) 3.949 1.165 ( -0.000 4.479 1.941) 4.881 1.219 ( -0.000 2.351 5.696) 6.162 1.253 ( -0.000 0.016 2.449) 2.449 1.255 ( -0.000 1.160 1.767) 2.114 1.277 ( -0.000 -1.568 0.847) 1.782 1.356 ( 0.000 -3.242 -1.816) 3.716 1.408 ( 0.000 -0.336 -1.920) 1.950 1.462 ( -0.000 -0.874 2.867) 2.997 1.471 ( 0.000 -1.385 -0.997) 1.706 1.483 ( 0.000 1.463 -2.495) 2.892 1.488 ( 0.000 -0.477 -0.639) 0.797 1.515 ( -0.000 -1.737 1.626) 2.379 1.540 ( 0.000 1.426 -2.427) 2.815 1.591 ( -0.000 1.637 0.762) 1.806 1.613 ( -0.000 -0.684 1.476) 1.627 1.616 ( -0.000 0.748 0.791) 1.089 1.650 ( 0.000 -3.212 -0.389) 3.235 1.667 ( 0.000 -1.564 -0.704) 1.715 1.670 ( 0.000 0.616 -4.772) 4.812 1.733 ( -0.000 -0.089 1.952) 1.954 1.762 ( 0.000 -2.554 -2.737) 3.743 1.831 ( -0.000 -0.135 0.982) 0.992 1.849 ( 0.000 2.094 -1.251) 2.439 1.854 ( -0.000 0.672 0.572) 0.883 1.863 ( 0.000 -1.591 -0.428) 1.647 1.881 ( 0.000 2.473 -2.301) 3.378 1.892 ( 0.000 -0.988 -0.203) 1.008 1.945 ( 0.000 0.293 -1.897) 1.920 1.957 ( 0.000 1.004 -2.090) 2.319 1.988 ( -0.000 3.237 3.010) 4.420 2.022 ( -0.000 -0.061 1.677) 1.678 2.083 ( -0.000 3.023 2.949) 4.223 2.090 ( -0.000 0.358 4.988) 5.001 2.149 ( -0.000 -1.548 1.628) 2.247 2.246 ( 0.000 -0.005 -1.362) 1.362 2.281 ( 0.000 0.393 -5.727) 5.740 2.285 ( -0.000 -1.348 0.663) 1.503 2.334 ( 0.000 -1.220 -3.402) 3.614 2.450 ( 0.000 0.119 -0.476) 0.491 2.486 ( 0.000 0.311 -1.075) 1.119 2.486 ( -0.000 0.696 2.694) 2.782 2.498 ( 0.000 0.204 -0.764) 0.791 2.530 ( -0.000 1.262 0.999) 1.610 2.634 ( -0.000 -2.070 0.740) 2.199 2.798 ( 0.000 0.911 -5.037) 5.119 3.008 ( -0.000 0.986 2.768) 2.939 3.132 ( -0.000 2.868 0.641) 2.939 3.156 ( 0.000 -2.967 -1.857) 3.500 3.162 ( 0.000 -1.019 0.303) 1.063 3.283 ( 0.000 1.289 -0.823) 1.530 3.360 ( -0.000 0.088 3.399) 3.400 3.502 ( 0.000 -4.307 -0.223) 4.312 3.533 ( -0.000 2.031 0.867) 2.208 3.535 ( -0.000 4.179 2.376) 4.808 3.816 ( -0.000 -0.198 0.242) 0.313 3.832 ( -0.000 0.205 0.494) 0.535 3.838 ( 0.000 -2.219 -1.016) 2.440 3.928 ( -0.000 0.121 0.642) 0.654 ======================= Grid point 13 (9/19) ======================= q-point: (-0.40 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.615 ( -7.509 -0.000 0.000) 7.509 0.668 ( -5.041 -0.000 0.000) 5.041 0.745 ( 1.135 0.000 -0.000) 1.135 0.791 ( 1.847 0.000 -0.000) 1.847 0.821 ( -1.643 -0.000 0.000) 1.643 1.003 ( -2.931 -0.000 0.000) 2.931 1.053 ( -2.793 -0.000 0.000) 2.793 1.069 ( 6.192 0.000 -0.000) 6.192 1.193 ( -4.653 -0.000 0.000) 4.653 1.208 ( 2.869 0.000 -0.000) 2.869 1.228 ( -3.995 -0.000 0.000) 3.995 1.267 ( -0.837 -0.000 0.000) 0.837 1.297 ( 0.388 0.000 -0.000) 0.388 1.347 ( 0.630 0.000 -0.000) 0.630 1.386 ( -0.277 -0.000 0.000) 0.277 1.414 ( -1.067 -0.000 0.000) 1.067 1.443 ( 0.565 0.000 -0.000) 0.565 1.482 ( -1.548 -0.000 0.000) 1.548 1.496 ( -2.940 -0.000 0.000) 2.940 1.516 ( 0.460 0.000 -0.000) 0.460 1.548 ( -5.021 -0.000 0.000) 5.021 1.579 ( 3.355 0.000 -0.000) 3.355 1.605 ( -1.110 -0.000 0.000) 1.110 1.625 ( 1.475 0.000 -0.000) 1.475 1.664 ( 2.269 0.000 -0.000) 2.269 1.667 ( 1.748 0.000 -0.000) 1.748 1.685 ( 3.513 0.000 -0.000) 3.513 1.744 ( 1.544 0.000 -0.000) 1.544 1.760 ( 1.955 0.000 -0.000) 1.955 1.790 ( -3.877 -0.000 0.000) 3.877 1.837 ( -0.844 -0.000 0.000) 0.844 1.856 ( -0.397 -0.000 0.000) 0.397 1.889 ( 0.098 0.000 -0.000) 0.098 1.897 ( 1.483 0.000 -0.000) 1.483 1.903 ( 2.272 0.000 -0.000) 2.272 1.934 ( 0.946 0.000 -0.000) 0.946 1.959 ( 1.604 0.000 -0.000) 1.604 2.013 ( -1.609 -0.000 0.000) 1.609 2.029 ( 0.488 0.000 -0.000) 0.488 2.106 ( -1.850 -0.000 0.000) 1.850 2.162 ( 0.190 0.000 -0.000) 0.190 2.178 ( 0.738 0.000 -0.000) 0.738 2.191 ( 1.688 0.000 -0.000) 1.688 2.241 ( 0.596 0.000 -0.000) 0.596 2.292 ( -0.436 -0.000 0.000) 0.436 2.308 ( 0.633 0.000 -0.000) 0.633 2.449 ( 0.891 0.000 -0.000) 0.891 2.483 ( 0.062 0.000 -0.000) 0.062 2.503 ( 0.943 0.000 -0.000) 0.943 2.511 ( 0.238 0.000 -0.000) 0.238 2.546 ( 0.436 0.000 -0.000) 0.436 2.629 ( 0.105 0.000 -0.000) 0.105 2.764 ( -0.702 -0.000 0.000) 0.702 3.048 ( 1.320 0.000 -0.000) 1.320 3.129 ( 0.801 0.000 -0.000) 0.801 3.139 ( -0.622 -0.000 0.000) 0.622 3.165 ( -0.235 -0.000 0.000) 0.235 3.280 ( -0.600 -0.000 0.000) 0.600 3.367 ( -2.617 -0.000 0.000) 2.617 3.449 ( -2.740 -0.000 0.000) 2.740 3.514 ( -2.742 -0.000 0.000) 2.742 3.595 ( 0.963 0.000 -0.000) 0.963 3.777 ( -3.225 -0.000 0.000) 3.225 3.816 ( -0.053 -0.000 0.000) 0.053 3.850 ( 0.653 0.000 -0.000) 0.653 3.952 ( 1.666 0.000 -0.000) 1.666 ======================= Grid point 31 (10/19) ======================= q-point: ( 0.20 0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 19 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.516 ( -0.000 0.000 12.747) 12.747 0.516 ( -0.000 0.000 12.747) 12.747 0.809 ( -0.000 0.000 5.033) 5.033 0.809 ( -0.000 0.000 5.033) 5.033 0.857 ( -0.000 0.000 18.986) 18.986 0.967 ( 0.000 -0.000 -4.252) 4.252 1.019 ( 0.000 -0.000 -2.856) 2.856 1.107 ( -0.000 0.000 2.832) 2.832 1.107 ( -0.000 0.000 2.832) 2.832 1.118 ( 0.000 -0.000 -5.808) 5.808 1.197 ( 0.000 -0.000 -3.680) 3.680 1.197 ( 0.000 -0.000 -3.680) 3.680 1.327 ( -0.000 0.000 3.307) 3.307 1.368 ( 0.000 -0.000 -0.914) 0.914 1.368 ( 0.000 -0.000 -0.914) 0.914 1.506 ( 0.000 -0.000 -0.149) 0.149 1.506 ( 0.000 -0.000 -0.149) 0.149 1.539 ( 0.000 -0.000 -2.544) 2.544 1.539 ( 0.000 -0.000 -2.544) 2.544 1.548 ( 0.000 -0.000 -0.803) 0.803 1.569 ( -0.000 0.000 0.359) 0.359 1.604 ( 0.000 0.000 0.035) 0.035 1.604 ( 0.000 -0.000 0.035) 0.035 1.639 ( -0.000 0.000 1.268) 1.268 1.681 ( -0.000 0.000 5.761) 5.761 1.683 ( 0.000 -0.000 -0.513) 0.513 1.729 ( -0.000 0.000 1.928) 1.928 1.729 ( -0.000 0.000 1.928) 1.928 1.800 ( 0.000 -0.000 -4.283) 4.283 1.800 ( 0.000 -0.000 -4.283) 4.283 1.837 ( -0.000 0.000 1.785) 1.785 1.862 ( -0.000 -0.000 -0.232) 0.232 1.862 ( 0.000 -0.000 -0.232) 0.232 1.884 ( 0.000 -0.000 -2.442) 2.442 1.884 ( 0.000 -0.000 -2.442) 2.442 1.906 ( -0.000 0.000 2.747) 2.747 1.962 ( 0.000 -0.000 -1.502) 1.502 1.962 ( 0.000 -0.000 -1.502) 1.502 1.992 ( -0.000 0.000 1.791) 1.791 2.018 ( 0.000 -0.000 -0.903) 0.903 2.018 ( 0.000 -0.000 -0.903) 0.903 2.203 ( 0.000 -0.000 -4.065) 4.065 2.204 ( -0.000 0.000 2.308) 2.308 2.204 ( -0.000 0.000 2.308) 2.308 2.356 ( -0.000 0.000 1.370) 1.370 2.356 ( -0.000 0.000 1.370) 1.370 2.456 ( 0.000 -0.000 -1.950) 1.950 2.465 ( 0.000 -0.000 -0.788) 0.788 2.465 ( 0.000 -0.000 -0.788) 0.788 2.511 ( 0.000 -0.000 -2.225) 2.225 2.511 ( 0.000 -0.000 -2.225) 2.225 2.680 ( -0.000 0.000 1.646) 1.646 2.847 ( 0.000 -0.000 -7.412) 7.412 2.997 ( -0.000 0.000 2.241) 2.241 2.997 ( -0.000 0.000 2.241) 2.241 3.236 ( 0.000 -0.000 -1.862) 1.862 3.236 ( 0.000 -0.000 -1.862) 1.862 3.242 ( 0.000 -0.000 -2.553) 2.553 3.387 ( -0.000 0.000 2.489) 2.489 3.387 ( -0.000 0.000 2.489) 2.489 3.434 ( -0.000 0.000 5.366) 5.366 3.684 ( -0.000 0.000 5.389) 5.389 3.684 ( -0.000 0.000 5.389) 5.389 3.815 ( -0.000 0.000 0.713) 0.713 3.832 ( 0.000 -0.000 -1.111) 1.111 4.028 ( 0.000 -0.000 -3.839) 3.839 ======================= Grid point 32 (11/19) ======================= q-point: ( 0.40 0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.634 ( -1.492 -0.862 10.491) 10.631 0.689 ( 1.667 0.962 10.907) 11.076 0.750 ( -7.006 -4.045 -0.025) 8.089 0.811 ( -7.289 -4.208 2.465) 8.770 0.976 ( 2.866 1.654 -0.866) 3.420 1.010 ( 2.481 1.433 1.009) 3.038 1.012 ( -3.205 -1.850 -3.792) 5.299 1.069 ( 1.241 0.717 5.217) 5.410 1.071 ( 1.293 0.746 -4.526) 4.765 1.150 ( 1.295 0.748 -1.019) 1.810 1.194 ( -0.508 -0.293 3.085) 3.140 1.220 ( 0.321 0.185 5.132) 5.145 1.343 ( -0.132 -0.076 -0.996) 1.007 1.376 ( 0.177 0.102 -3.584) 3.590 1.396 ( 2.516 1.453 4.118) 5.039 1.465 ( -1.169 -0.675 -1.767) 2.223 1.492 ( -1.210 -0.698 1.525) 2.068 1.498 ( -0.586 -0.339 0.942) 1.160 1.511 ( -2.389 -1.379 -1.175) 2.999 1.520 ( -1.734 -1.001 -2.244) 3.007 1.543 ( 0.371 0.214 -2.319) 2.358 1.587 ( 2.223 1.284 0.621) 2.641 1.615 ( -6.892 -3.979 -3.945) 8.882 1.644 ( -0.943 -0.544 -1.197) 1.618 1.677 ( 0.636 0.367 -1.513) 1.682 1.681 ( 3.686 2.128 -1.389) 4.478 1.691 ( -1.390 -0.802 0.690) 1.747 1.756 ( -0.785 -0.453 -0.775) 1.193 1.774 ( 1.013 0.585 0.604) 1.316 1.796 ( -1.015 -0.586 0.988) 1.533 1.832 ( 0.035 0.020 1.729) 1.730 1.848 ( 1.923 1.110 -1.122) 2.488 1.870 ( 3.046 1.759 -0.459) 3.547 1.880 ( -1.443 -0.833 0.251) 1.685 1.906 ( -1.919 -1.108 -0.813) 2.360 1.938 ( 0.992 0.573 4.224) 4.377 1.948 ( 1.856 1.071 -1.436) 2.579 1.971 ( 0.092 0.053 1.118) 1.123 1.973 ( 0.273 0.158 -1.085) 1.130 2.054 ( 0.339 0.196 5.816) 5.830 2.065 ( 1.931 1.115 -0.033) 2.230 2.173 ( 1.549 0.894 -0.755) 1.942 2.233 ( 0.685 0.396 -1.604) 1.789 2.246 ( 1.898 1.096 0.828) 2.342 2.327 ( -0.988 -0.571 -2.842) 3.062 2.346 ( -0.485 -0.280 -3.457) 3.503 2.438 ( 0.790 0.456 -1.601) 1.843 2.479 ( 0.380 0.219 -0.164) 0.468 2.479 ( 0.538 0.311 0.734) 0.961 2.493 ( 0.410 0.237 1.264) 1.350 2.502 ( 1.101 0.636 -0.068) 1.273 2.667 ( -1.395 -0.805 -0.454) 1.673 2.782 ( 0.237 0.137 -3.942) 3.951 3.028 ( 1.371 0.791 1.587) 2.241 3.041 ( 0.506 0.292 2.493) 2.560 3.197 ( -0.676 -0.390 -1.752) 1.918 3.199 ( -0.412 -0.238 -2.054) 2.108 3.241 ( 1.098 0.634 -2.506) 2.808 3.390 ( -1.856 -1.072 2.268) 3.120 3.444 ( -0.356 -0.206 5.496) 5.511 3.508 ( 0.998 0.576 4.354) 4.504 3.606 ( -6.957 -4.017 2.872) 8.532 3.816 ( -0.286 -0.165 -1.104) 1.152 3.821 ( 0.963 0.556 1.659) 1.997 3.851 ( 4.246 2.452 2.422) 5.469 3.943 ( -3.007 -1.736 -4.619) 5.779 ======================= Grid point 33 (12/19) ======================= q-point: (-0.40 0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.566 ( -9.246 -5.338 1.813) 10.829 0.634 ( -8.979 -5.184 -0.078) 10.368 0.758 ( 1.283 0.741 7.008) 7.163 0.794 ( -3.981 -2.299 2.093) 5.051 0.895 ( -2.883 -1.665 0.300) 3.343 0.998 ( -3.622 -2.091 3.380) 5.377 1.029 ( -3.701 -2.137 -2.399) 4.901 1.077 ( 2.487 1.436 -2.331) 3.699 1.080 ( 2.177 1.257 -0.738) 2.620 1.150 ( 1.969 1.137 -2.151) 3.130 1.247 ( 2.015 1.163 4.841) 5.371 1.302 ( -0.400 -0.231 0.990) 1.092 1.312 ( -0.815 -0.471 -3.467) 3.592 1.340 ( 1.301 0.751 2.665) 3.060 1.409 ( -2.177 -1.257 -3.405) 4.232 1.443 ( -2.501 -1.444 1.895) 3.454 1.444 ( -1.410 -0.814 -0.559) 1.721 1.464 ( 0.349 0.201 -2.125) 2.163 1.470 ( 1.248 0.720 -1.187) 1.867 1.480 ( -1.911 -1.103 -0.248) 2.220 1.538 ( 1.235 0.713 -1.287) 1.922 1.603 ( 0.608 0.351 -0.590) 0.917 1.616 ( -0.131 -0.076 -2.032) 2.038 1.624 ( 2.474 1.428 4.127) 5.020 1.638 ( -1.968 -1.136 -1.192) 2.566 1.699 ( -1.374 -0.793 -0.409) 1.638 1.727 ( 2.994 1.729 1.797) 3.897 1.737 ( 3.281 1.894 -1.936) 4.254 1.776 ( -2.422 -1.398 -1.943) 3.405 1.790 ( 0.480 0.277 -1.044) 1.182 1.848 ( -0.545 -0.314 0.811) 1.026 1.852 ( -0.924 -0.533 3.001) 3.185 1.878 ( 1.671 0.965 1.154) 2.249 1.884 ( -0.387 -0.223 2.420) 2.461 1.939 ( -0.822 -0.474 -1.116) 1.465 1.964 ( 1.330 0.768 -0.152) 1.544 1.972 ( -1.668 -0.963 2.317) 3.012 1.999 ( 4.109 2.372 0.293) 4.753 1.999 ( -0.081 -0.047 0.529) 0.537 2.093 ( 1.786 1.031 1.845) 2.767 2.135 ( -0.050 -0.029 7.084) 7.084 2.172 ( -0.070 -0.041 -7.627) 7.627 2.237 ( 2.667 1.540 -2.205) 3.787 2.253 ( -0.859 -0.496 -5.132) 5.227 2.287 ( -0.543 -0.313 1.822) 1.927 2.304 ( 1.452 0.838 -0.819) 1.866 2.451 ( 2.097 1.211 -0.525) 2.478 2.479 ( 0.279 0.161 0.145) 0.354 2.500 ( 0.888 0.513 0.264) 1.059 2.530 ( 0.493 0.285 1.131) 1.267 2.536 ( 1.061 0.613 0.730) 1.427 2.635 ( -0.452 -0.261 -1.770) 1.846 2.757 ( -0.834 -0.481 1.340) 1.650 3.076 ( 2.147 1.240 0.169) 2.485 3.079 ( 0.639 0.369 1.443) 1.621 3.163 ( -2.268 -1.309 0.603) 2.687 3.172 ( 0.504 0.291 -1.904) 1.991 3.213 ( 0.178 0.103 -4.027) 4.033 3.373 ( -1.264 -0.730 0.666) 1.604 3.483 ( -2.294 -1.324 4.698) 5.394 3.491 ( -6.577 -3.797 1.575) 7.756 3.579 ( -0.005 -0.003 5.649) 5.649 3.770 ( -3.778 -2.181 -3.766) 5.763 3.829 ( -1.787 -1.032 -1.754) 2.708 3.844 ( 0.450 0.260 0.674) 0.851 3.951 ( 3.747 2.163 -0.356) 4.341 ======================= Grid point 34 (13/19) ======================= q-point: (-0.20 0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.429 ( -8.430 -4.867 6.090) 11.482 0.442 ( -10.419 -6.015 -0.471) 12.040 0.682 ( -4.041 -2.333 0.144) 4.668 0.783 ( -10.688 -6.170 10.381) 16.127 0.867 ( -0.633 -0.365 7.310) 7.347 0.906 ( -1.486 -0.858 -5.493) 5.754 1.065 ( -0.029 -0.017 7.317) 7.317 1.066 ( -0.370 -0.214 -3.307) 3.334 1.109 ( 2.495 1.440 -2.571) 3.862 1.196 ( 3.581 2.067 -2.409) 4.785 1.232 ( -1.209 -0.698 -3.358) 3.636 1.273 ( -0.906 -0.523 -3.714) 3.858 1.294 ( 1.481 0.855 -1.772) 2.462 1.374 ( 0.859 0.496 -0.691) 1.209 1.379 ( -0.290 -0.167 0.412) 0.531 1.436 ( 1.024 0.591 0.154) 1.192 1.454 ( 1.688 0.974 0.979) 2.181 1.471 ( -0.302 -0.174 2.825) 2.846 1.491 ( 0.997 0.575 1.768) 2.110 1.519 ( 2.987 1.724 2.697) 4.378 1.559 ( 0.503 0.290 1.429) 1.542 1.579 ( 0.182 0.105 -0.589) 0.625 1.600 ( 1.300 0.751 -1.871) 2.399 1.616 ( 0.252 0.146 0.420) 0.511 1.649 ( -3.162 -1.825 -0.127) 3.653 1.686 ( -1.716 -0.991 -1.385) 2.418 1.753 ( 3.214 1.855 1.800) 4.124 1.754 ( -1.336 -0.772 -3.039) 3.408 1.774 ( 2.469 1.426 -1.302) 3.134 1.799 ( 1.531 0.884 1.072) 2.067 1.852 ( -0.560 -0.323 1.351) 1.498 1.869 ( -1.774 -1.024 3.544) 4.094 1.880 ( -2.879 -1.662 1.042) 3.484 1.883 ( -2.969 -1.714 -0.743) 3.508 1.928 ( 1.607 0.928 -0.657) 1.968 1.953 ( -1.683 -0.972 0.239) 1.958 1.980 ( 1.496 0.864 0.301) 1.754 1.991 ( 0.379 0.219 0.229) 0.494 2.062 ( -0.783 -0.452 -7.976) 8.027 2.072 ( 1.683 0.971 -1.600) 2.517 2.155 ( 1.689 0.975 -2.552) 3.212 2.168 ( -1.114 -0.643 -2.494) 2.807 2.194 ( 0.055 0.032 -2.742) 2.742 2.232 ( 0.983 0.567 6.717) 6.812 2.316 ( 1.101 0.636 -2.685) 2.971 2.349 ( 0.346 0.200 3.425) 3.449 2.489 ( 0.228 0.132 0.327) 0.420 2.501 ( 2.207 1.274 -1.360) 2.889 2.518 ( 0.817 0.472 0.896) 1.301 2.529 ( -0.231 -0.133 -2.314) 2.329 2.555 ( 0.478 0.276 1.199) 1.320 2.583 ( -1.124 -0.649 -3.511) 3.744 2.816 ( 0.767 0.443 8.738) 8.783 3.064 ( -0.987 -0.570 -3.642) 3.816 3.099 ( -4.242 -2.449 -0.612) 4.936 3.121 ( 2.501 1.444 0.098) 2.890 3.190 ( 2.861 1.652 -0.329) 3.320 3.218 ( 2.426 1.401 1.853) 3.359 3.341 ( -1.919 -1.108 -1.670) 2.775 3.425 ( -1.855 -1.071 3.406) 4.024 3.463 ( -3.874 -2.237 2.376) 5.065 3.606 ( -5.866 -3.387 -3.977) 7.855 3.636 ( -1.987 -1.147 6.895) 7.267 3.811 ( -0.325 -0.188 0.234) 0.442 3.852 ( -0.041 -0.023 -0.024) 0.053 3.998 ( 3.059 1.766 -2.500) 4.327 ======================= Grid point 37 (14/19) ======================= q-point: ( 0.40 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.758 ( -0.000 -2.452 5.944) 6.430 0.769 ( -0.000 -3.954 8.831) 9.676 0.792 ( -0.000 -6.126 2.891) 6.774 0.891 ( -0.000 -2.162 0.760) 2.292 0.898 ( 0.000 -3.640 -2.114) 4.209 0.962 ( 0.000 -1.372 -1.268) 1.869 1.002 ( 0.000 0.467 -3.706) 3.736 1.014 ( -0.000 3.679 0.286) 3.690 1.137 ( 0.000 5.582 -3.843) 6.777 1.160 ( -0.000 1.840 4.449) 4.815 1.226 ( 0.000 0.719 -3.255) 3.334 1.298 ( -0.000 1.873 6.950) 7.198 1.318 ( 0.000 -0.419 -0.034) 0.421 1.333 ( 0.000 1.269 -5.374) 5.521 1.409 ( -0.000 -2.560 0.958) 2.733 1.450 ( 0.000 -0.742 -3.541) 3.618 1.460 ( 0.000 -0.360 -0.957) 1.022 1.508 ( 0.000 -1.962 -0.907) 2.162 1.526 ( -0.000 0.584 1.794) 1.886 1.529 ( -0.000 -0.452 2.615) 2.654 1.558 ( -0.000 -2.809 0.914) 2.954 1.583 ( -0.000 1.520 0.427) 1.579 1.600 ( 0.000 0.383 -1.057) 1.124 1.629 ( -0.000 -0.559 2.257) 2.325 1.643 ( 0.000 -1.304 -0.653) 1.458 1.699 ( 0.000 -3.285 -0.227) 3.292 1.701 ( 0.000 -1.871 -1.869) 2.644 1.729 ( 0.000 0.646 -1.251) 1.408 1.730 ( 0.000 1.146 -2.831) 3.054 1.808 ( -0.000 0.133 1.291) 1.298 1.828 ( 0.000 2.031 -1.337) 2.431 1.829 ( 0.000 1.337 -0.694) 1.507 1.867 ( 0.000 -0.104 -1.152) 1.157 1.892 ( -0.000 -2.059 2.271) 3.066 1.912 ( -0.000 1.498 0.203) 1.512 1.920 ( -0.000 1.452 0.023) 1.452 1.972 ( -0.000 -2.459 3.592) 4.353 1.993 ( -0.000 0.415 4.364) 4.384 2.025 ( 0.000 0.937 -0.323) 0.991 2.049 ( -0.000 2.620 0.306) 2.638 2.143 ( -0.000 1.376 6.743) 6.882 2.152 ( -0.000 1.598 -0.108) 1.601 2.199 ( 0.000 -0.357 -6.249) 6.259 2.223 ( 0.000 0.820 -1.322) 1.556 2.300 ( 0.000 -0.163 -1.253) 1.263 2.306 ( 0.000 -0.374 -3.583) 3.602 2.418 ( 0.000 1.367 -1.751) 2.221 2.474 ( -0.000 0.154 0.320) 0.355 2.483 ( -0.000 0.624 0.244) 0.670 2.506 ( -0.000 0.077 1.697) 1.699 2.524 ( -0.000 1.193 1.547) 1.953 2.640 ( 0.000 -0.469 -3.858) 3.887 2.772 ( -0.000 -0.238 1.941) 1.956 3.034 ( -0.000 -0.590 2.261) 2.337 3.091 ( -0.000 2.945 1.942) 3.527 3.149 ( 0.000 -0.689 -4.539) 4.591 3.171 ( 0.000 -1.012 -1.352) 1.689 3.229 ( 0.000 2.000 -1.668) 2.604 3.436 ( -0.000 -2.775 4.436) 5.233 3.536 ( -0.000 2.202 6.105) 6.490 3.548 ( -0.000 -0.744 3.603) 3.679 3.630 ( 0.000 -7.729 0.557) 7.749 3.804 ( 0.000 -0.074 -1.385) 1.387 3.833 ( -0.000 0.218 0.260) 0.339 3.855 ( -0.000 2.295 0.449) 2.339 3.896 ( 0.000 -0.092 -4.196) 4.197 ======================= Grid point 38 (15/19) ======================= q-point: (-0.40 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 63 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.691 ( -7.591 -4.432 3.426) 9.434 0.714 ( -5.703 -5.197 0.671) 7.745 0.797 ( 0.414 -3.672 2.350) 4.379 0.825 ( -4.103 -4.061 -1.723) 6.025 0.925 ( 2.387 -5.042 4.769) 7.339 0.976 ( -0.641 2.024 0.269) 2.140 1.029 ( 3.943 1.398 -2.020) 4.646 1.036 ( 3.755 0.662 0.656) 3.869 1.116 ( -3.618 5.248 -1.108) 6.470 1.182 ( -2.480 3.630 -2.652) 5.135 1.227 ( 1.742 0.399 -1.293) 2.206 1.264 ( 0.284 0.270 -6.243) 6.255 1.300 ( -2.082 -1.731 2.761) 3.868 1.373 ( 0.258 -2.327 0.952) 2.528 1.386 ( -0.291 -1.278 0.119) 1.316 1.435 ( 0.603 -0.386 0.215) 0.748 1.455 ( 0.358 -0.415 1.271) 1.384 1.474 ( -3.037 0.138 -1.016) 3.206 1.497 ( -0.539 -0.154 -0.796) 0.974 1.507 ( -0.350 -0.044 0.700) 0.783 1.532 ( -4.014 0.366 0.848) 4.119 1.593 ( -1.248 0.476 0.105) 1.340 1.616 ( -0.917 1.042 -1.218) 1.846 1.633 ( 2.240 -1.426 1.268) 2.942 1.645 ( -1.252 -1.621 -0.145) 2.053 1.684 ( 2.550 0.271 -0.083) 2.566 1.693 ( 3.005 0.195 1.565) 3.394 1.728 ( 0.362 0.277 0.198) 0.497 1.763 ( 2.551 0.629 -0.712) 2.723 1.784 ( -3.668 -0.335 -0.164) 3.687 1.836 ( 0.970 0.307 -1.234) 1.599 1.846 ( 0.063 -1.151 0.225) 1.174 1.885 ( 2.185 1.482 1.947) 3.281 1.898 ( 3.445 0.610 -1.473) 3.796 1.930 ( -0.574 0.572 1.480) 1.687 1.975 ( -0.563 -2.788 2.847) 4.024 1.990 ( 0.667 1.128 4.758) 4.935 2.011 ( -1.623 -0.725 -1.269) 2.184 2.026 ( -0.105 0.192 0.107) 0.244 2.073 ( 0.365 2.732 -4.227) 5.046 2.097 ( 0.167 -1.814 -2.041) 2.736 2.166 ( 1.544 -1.226 -3.557) 4.067 2.216 ( -0.804 1.661 -2.399) 3.026 2.232 ( -0.428 1.302 4.934) 5.121 2.283 ( 0.203 1.043 0.322) 1.110 2.306 ( 1.580 1.432 2.115) 3.003 2.428 ( 1.599 2.080 -1.135) 2.858 2.486 ( -0.098 -0.012 1.078) 1.082 2.499 ( 0.525 0.175 0.611) 0.825 2.515 ( 0.630 -0.755 -0.905) 1.336 2.543 ( 0.230 0.593 -0.073) 0.641 2.594 ( 0.287 0.884 -3.818) 3.930 2.816 ( -0.271 -0.558 6.283) 6.314 3.051 ( 0.810 -1.014 -1.711) 2.148 3.094 ( -0.451 1.392 -1.193) 1.888 3.137 ( 0.460 0.205 -0.295) 0.584 3.161 ( 1.299 1.284 -0.118) 1.831 3.217 ( -1.288 2.156 -1.706) 3.036 3.425 ( -3.982 -2.386 2.850) 5.447 3.509 ( -2.679 -7.001 -0.772) 7.536 3.573 ( -2.737 -2.287 3.575) 5.050 3.639 ( 1.442 -1.287 4.470) 4.870 3.762 ( -2.432 -0.067 -2.184) 3.269 3.822 ( -0.547 -0.691 -0.593) 1.062 3.845 ( 0.658 0.410 -0.129) 0.786 3.906 ( 2.017 2.913 -2.783) 4.506 ======================= Grid point 39 (16/19) ======================= q-point: (-0.20 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 63 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.572 ( -10.226 -2.845 1.640) 10.741 0.606 ( -8.437 -2.139 6.528) 10.880 0.697 ( 0.950 -5.752 0.392) 5.843 0.786 ( 2.873 -5.162 -5.758) 8.250 0.934 ( -1.845 -1.070 3.393) 4.008 0.992 ( -1.750 -0.803 4.137) 4.563 1.083 ( -1.037 1.931 2.877) 3.617 1.109 ( -1.595 4.582 0.388) 4.868 1.116 ( 2.380 0.630 -2.141) 3.263 1.148 ( -2.423 -1.155 -2.159) 3.444 1.156 ( -0.290 -0.139 -0.596) 0.677 1.223 ( 3.332 -1.676 -3.150) 4.882 1.299 ( -3.079 0.097 -0.420) 3.109 1.351 ( 2.810 -1.661 0.787) 3.358 1.385 ( -1.694 0.804 -0.883) 2.072 1.423 ( 2.635 -1.975 0.696) 3.366 1.436 ( -1.982 0.272 0.439) 2.048 1.481 ( 1.003 1.019 -0.438) 1.495 1.518 ( 0.552 -1.173 1.989) 2.374 1.522 ( 1.652 1.450 0.994) 2.413 1.545 ( 0.780 0.676 3.109) 3.276 1.599 ( 0.016 0.536 -1.535) 1.626 1.610 ( -0.661 1.571 0.391) 1.749 1.627 ( 0.658 2.133 -0.162) 2.238 1.667 ( -1.619 1.267 -2.141) 2.968 1.676 ( 0.796 -0.085 -1.759) 1.933 1.722 ( -0.416 0.785 0.164) 0.903 1.731 ( 1.548 -2.073 -0.215) 2.596 1.778 ( -0.079 0.330 1.568) 1.604 1.800 ( 0.477 1.146 2.075) 2.418 1.823 ( 1.669 -2.860 -2.268) 4.014 1.847 ( 1.302 -0.930 0.014) 1.600 1.909 ( -1.589 -1.254 0.754) 2.160 1.924 ( -2.315 -0.789 0.175) 2.452 1.943 ( -0.852 -0.538 -1.722) 1.996 1.952 ( 1.664 -0.676 -1.460) 2.314 1.978 ( -0.016 -0.693 2.231) 2.336 1.995 ( -0.646 -0.022 -1.582) 1.709 2.010 ( -2.651 1.143 0.956) 3.041 2.070 ( -0.267 -0.062 -2.387) 2.402 2.083 ( 1.012 0.647 -1.265) 1.744 2.122 ( 1.090 -0.943 -2.359) 2.764 2.219 ( -0.450 2.491 1.482) 2.934 2.261 ( 1.759 1.031 -1.437) 2.494 2.338 ( -1.066 0.946 5.314) 5.502 2.379 ( 1.161 0.380 2.639) 2.908 2.455 ( 1.823 1.179 -2.067) 2.997 2.462 ( 0.467 -0.019 -4.247) 4.272 2.509 ( -0.092 -0.360 1.228) 1.283 2.525 ( 0.258 0.254 -0.318) 0.482 2.555 ( 0.561 -0.062 -0.075) 0.570 2.558 ( 0.091 0.312 -1.079) 1.126 2.904 ( 0.462 -0.037 7.359) 7.374 2.987 ( -0.342 -0.366 -7.218) 7.235 3.109 ( -0.438 0.352 1.091) 1.227 3.142 ( 1.851 0.252 -0.640) 1.974 3.193 ( -2.260 2.800 -0.020) 3.599 3.223 ( 2.758 1.563 2.120) 3.813 3.362 ( -3.972 -1.138 -0.388) 4.150 3.417 ( 0.292 -6.280 0.896) 6.350 3.540 ( -3.876 -2.710 2.521) 5.360 3.589 ( -2.846 -4.825 -2.816) 6.269 3.766 ( -2.356 2.623 4.578) 5.778 3.817 ( 0.253 -0.115 0.187) 0.335 3.857 ( -0.587 0.357 0.752) 1.019 3.923 ( 3.585 1.391 -3.731) 5.358 ======================= Grid point 43 (17/19) ======================= q-point: (-0.40 -0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.663 ( 0.000 -9.756 -1.373) 9.852 0.716 ( 0.000 -6.620 -3.321) 7.406 0.789 ( -0.000 -2.828 4.235) 5.092 0.847 ( 0.000 -2.736 -7.601) 8.078 0.912 ( -0.000 -2.165 5.206) 5.638 0.980 ( -0.000 0.152 4.206) 4.209 1.005 ( -0.000 2.976 2.845) 4.117 1.056 ( -0.000 4.065 2.810) 4.942 1.150 ( 0.000 -0.410 -6.175) 6.188 1.177 ( 0.000 0.199 -4.526) 4.531 1.203 ( 0.000 4.387 -1.534) 4.647 1.212 ( 0.000 0.548 -2.399) 2.461 1.318 ( 0.000 -4.960 1.203) 5.104 1.360 ( -0.000 -1.996 1.382) 2.427 1.387 ( 0.000 0.708 -0.875) 1.125 1.452 ( -0.000 1.654 -0.149) 1.661 1.457 ( -0.000 0.034 0.254) 0.256 1.469 ( 0.000 -0.341 -1.660) 1.695 1.479 ( 0.000 0.496 -0.694) 0.853 1.520 ( -0.000 1.056 0.829) 1.343 1.584 ( -0.000 -0.426 0.598) 0.734 1.595 ( 0.000 -1.057 -0.620) 1.225 1.599 ( 0.000 -0.467 -1.240) 1.326 1.623 ( 0.000 -1.600 -0.042) 1.601 1.639 ( 0.000 -1.249 -0.016) 1.250 1.667 ( -0.000 0.688 1.666) 1.802 1.709 ( 0.000 -1.227 -1.001) 1.583 1.733 ( -0.000 2.076 0.705) 2.193 1.747 ( -0.000 1.054 3.288) 3.453 1.825 ( -0.000 1.539 0.807) 1.738 1.835 ( 0.000 0.550 -2.214) 2.281 1.849 ( -0.000 -0.087 0.297) 0.310 1.858 ( -0.000 0.654 4.091) 4.143 1.912 ( -0.000 -0.096 2.470) 2.472 1.954 ( 0.000 1.284 -0.402) 1.346 1.963 ( -0.000 -0.585 0.740) 0.943 1.977 ( 0.000 -0.206 -4.387) 4.391 2.028 ( 0.000 -2.224 -3.095) 3.811 2.032 ( -0.000 -2.367 3.991) 4.640 2.056 ( -0.000 -1.677 3.149) 3.567 2.079 ( 0.000 -0.947 -2.824) 2.979 2.102 ( 0.000 1.246 -0.953) 1.569 2.228 ( -0.000 2.564 0.178) 2.570 2.234 ( 0.000 0.814 -3.655) 3.745 2.345 ( -0.000 2.052 2.341) 3.113 2.355 ( -0.000 1.523 5.758) 5.956 2.419 ( -0.000 1.804 -0.542) 1.884 2.454 ( 0.000 -0.344 -4.334) 4.348 2.511 ( -0.000 -0.080 0.953) 0.957 2.513 ( -0.000 0.022 2.245) 2.245 2.535 ( 0.000 -0.333 -0.668) 0.747 2.573 ( 0.000 2.193 -2.207) 3.112 2.898 ( -0.000 -0.784 7.840) 7.880 2.983 ( 0.000 -1.013 -5.467) 5.560 3.118 ( -0.000 2.022 0.858) 2.197 3.127 ( 0.000 0.993 -1.733) 1.997 3.159 ( -0.000 -0.242 0.727) 0.766 3.227 ( 0.000 3.523 -2.158) 4.131 3.423 ( 0.000 -7.638 -1.431) 7.771 3.494 ( -0.000 -2.484 1.198) 2.758 3.606 ( 0.000 -5.177 -1.248) 5.326 3.611 ( -0.000 -3.867 3.454) 5.185 3.782 ( 0.000 1.046 -1.262) 1.639 3.828 ( 0.000 -0.163 -0.568) 0.591 3.837 ( -0.000 1.646 1.878) 2.497 3.881 ( 0.000 2.677 -2.564) 3.706 ======================= Grid point 62 (18/19) ======================= q-point: ( 0.40 0.40 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 19 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.885 ( -0.000 0.000 3.683) 3.683 0.885 ( -0.000 0.000 3.683) 3.683 0.919 ( -0.000 0.000 0.615) 0.615 0.929 ( 0.000 -0.000 -1.179) 1.179 0.934 ( -0.000 0.000 5.608) 5.608 0.934 ( -0.000 0.000 5.608) 5.608 0.982 ( -0.000 0.000 0.332) 0.332 0.988 ( 0.000 -0.000 -3.920) 3.920 0.988 ( 0.000 -0.000 -3.920) 3.920 1.128 ( 0.000 -0.000 -3.157) 3.157 1.208 ( 0.000 -0.000 -6.611) 6.611 1.208 ( 0.000 -0.000 -6.611) 6.611 1.358 ( -0.000 0.000 5.452) 5.452 1.358 ( -0.000 0.000 5.452) 5.452 1.443 ( 0.000 -0.000 -0.510) 0.510 1.443 ( 0.000 -0.000 -0.510) 0.510 1.466 ( -0.000 0.000 1.982) 1.982 1.522 ( 0.000 -0.000 -1.274) 1.274 1.533 ( 0.000 -0.000 -1.152) 1.152 1.533 ( 0.000 -0.000 -1.152) 1.152 1.543 ( 0.000 -0.000 -1.585) 1.585 1.543 ( 0.000 -0.000 -1.585) 1.585 1.653 ( 0.000 -0.000 -0.957) 0.957 1.654 ( 0.000 -0.000 -1.060) 1.060 1.685 ( -0.000 0.000 6.025) 6.025 1.687 ( -0.000 0.000 1.256) 1.256 1.687 ( 0.000 0.000 1.256) 1.256 1.748 ( 0.000 -0.000 -2.354) 2.354 1.748 ( 0.000 -0.000 -2.354) 2.354 1.752 ( 0.000 -0.000 -0.175) 0.175 1.821 ( -0.000 0.000 2.450) 2.450 1.834 ( -0.000 0.000 0.746) 0.746 1.834 ( -0.000 0.000 0.746) 0.746 1.920 ( -0.000 0.000 3.937) 3.937 1.920 ( -0.000 0.000 3.937) 3.937 1.939 ( 0.000 -0.000 -3.452) 3.452 1.939 ( 0.000 -0.000 -3.452) 3.452 2.006 ( 0.000 -0.000 -3.780) 3.780 2.053 ( -0.000 0.000 2.619) 2.619 2.053 ( -0.000 0.000 2.619) 2.619 2.161 ( 0.000 -0.000 -3.190) 3.190 2.161 ( 0.000 -0.000 -3.190) 3.190 2.168 ( -0.000 0.000 5.508) 5.508 2.224 ( -0.000 0.000 7.924) 7.924 2.264 ( 0.000 -0.000 -3.842) 3.842 2.264 ( 0.000 -0.000 -3.842) 3.842 2.375 ( 0.000 -0.000 -2.275) 2.275 2.489 ( -0.000 0.000 1.688) 1.688 2.489 ( -0.000 0.000 1.688) 1.688 2.526 ( -0.000 0.000 0.714) 0.714 2.526 ( -0.000 0.000 0.714) 0.714 2.571 ( 0.000 -0.000 -7.042) 7.042 2.827 ( -0.000 0.000 5.881) 5.881 3.081 ( 0.000 -0.000 -6.257) 6.257 3.083 ( -0.000 0.000 2.600) 2.600 3.083 ( -0.000 0.000 2.600) 2.600 3.175 ( 0.000 -0.000 -2.041) 2.041 3.175 ( 0.000 -0.000 -2.041) 2.041 3.567 ( -0.000 0.000 6.473) 6.473 3.567 ( -0.000 0.000 6.473) 6.473 3.611 ( -0.000 0.000 4.870) 4.870 3.751 ( 0.000 0.000 -2.190) 2.190 3.751 ( 0.000 -0.000 -2.190) 2.190 3.786 ( 0.000 -0.000 -2.518) 2.518 3.831 ( 0.000 -0.000 -0.082) 0.082 3.853 ( 0.000 -0.000 -3.760) 3.760 ======================= Grid point 63 (19/19) ======================= q-point: (-0.40 0.40 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.56e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.94e-05 1.94e-05 Number of triplets: 39 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.810 ( -6.754 -3.899 3.122) 8.400 0.830 ( -7.583 -4.378 2.017) 8.985 0.866 ( -4.439 -2.563 -0.855) 5.196 0.878 ( -3.586 -2.070 -1.224) 4.318 0.937 ( 1.965 1.135 -0.611) 2.350 0.978 ( 4.221 2.437 -0.236) 4.880 1.027 ( 0.910 0.525 3.086) 3.260 1.041 ( -2.529 -1.460 -3.891) 4.864 1.071 ( 2.125 1.227 2.962) 3.847 1.107 ( 3.869 2.234 4.317) 6.213 1.153 ( 2.689 1.552 -4.451) 5.427 1.166 ( 2.922 1.687 -6.500) 7.323 1.374 ( -2.884 -1.665 3.740) 5.008 1.411 ( -1.355 -0.782 0.674) 1.703 1.425 ( -0.216 -0.125 0.464) 0.526 1.445 ( -0.205 -0.118 -3.163) 3.172 1.449 ( -3.613 -2.086 0.385) 4.190 1.476 ( -0.460 -0.266 0.907) 1.051 1.491 ( -0.308 -0.178 -1.572) 1.611 1.517 ( -0.747 -0.431 -0.910) 1.254 1.541 ( 0.780 0.450 3.354) 3.473 1.586 ( 1.863 1.075 -0.350) 2.179 1.618 ( -2.335 -1.348 0.812) 2.816 1.658 ( 0.413 0.239 -1.994) 2.051 1.666 ( -1.775 -1.025 1.226) 2.388 1.687 ( -0.873 -0.504 0.185) 1.024 1.702 ( -0.866 -0.500 0.358) 1.062 1.720 ( 0.040 0.023 0.061) 0.077 1.749 ( 0.793 0.458 -0.760) 1.190 1.777 ( 2.353 1.358 0.581) 2.778 1.833 ( -2.314 -1.336 0.203) 2.680 1.864 ( 0.667 0.385 -1.815) 1.972 1.877 ( 1.830 1.057 2.641) 3.383 1.897 ( 1.640 0.947 4.983) 5.331 1.929 ( 0.583 0.337 -4.443) 4.494 1.957 ( 0.580 0.335 -4.675) 4.723 1.979 ( 1.096 0.633 1.898) 2.281 2.014 ( -0.933 -0.539 -3.452) 3.616 2.017 ( 2.295 1.325 4.842) 5.520 2.065 ( 0.374 0.216 -0.109) 0.445 2.138 ( -1.944 -1.122 -1.836) 2.900 2.148 ( -1.965 -1.135 1.498) 2.719 2.184 ( -0.516 -0.298 0.396) 0.716 2.208 ( 0.533 0.308 -2.796) 2.863 2.297 ( 3.477 2.007 2.520) 4.740 2.313 ( 1.713 0.989 1.860) 2.715 2.392 ( 1.168 0.675 3.166) 3.441 2.468 ( -0.468 -0.270 -6.482) 6.505 2.514 ( 0.234 0.135 1.350) 1.376 2.515 ( 0.234 0.135 1.028) 1.063 2.529 ( 0.596 0.344 -0.209) 0.720 2.538 ( 0.917 0.530 -1.867) 2.147 2.908 ( -0.145 -0.084 7.440) 7.442 2.998 ( -0.098 -0.057 -7.168) 7.169 3.098 ( -0.462 -0.266 2.319) 2.380 3.128 ( -0.479 -0.277 -2.153) 2.223 3.145 ( 1.639 0.946 2.389) 3.047 3.173 ( 1.169 0.675 -2.109) 2.505 3.529 ( -6.511 -3.759 4.272) 8.648 3.601 ( -5.284 -3.050 -3.081) 6.835 3.645 ( -2.461 -1.421 4.440) 5.272 3.717 ( 2.318 1.338 6.529) 7.056 3.721 ( -1.854 -1.071 -5.050) 5.486 3.802 ( 4.109 2.372 -4.556) 6.578 3.827 ( -0.256 -0.148 -0.034) 0.297 3.837 ( 0.917 0.529 -0.156) 1.070 =================== 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/8250 10.0 14.000 14.000 22.340 0.000 -0.000 -0.000 3/8250 20.0 2.325 2.325 3.831 0.000 -0.000 -0.000 3/8250 30.0 1.111 1.111 1.894 0.000 -0.000 -0.000 3/8250 40.0 0.735 0.735 1.285 0.000 -0.000 -0.000 3/8250 50.0 0.556 0.556 0.989 0.000 -0.000 -0.000 3/8250 60.0 0.451 0.451 0.810 0.000 -0.000 -0.000 3/8250 70.0 0.381 0.381 0.689 0.000 -0.000 -0.000 3/8250 80.0 0.330 0.330 0.601 0.000 -0.000 -0.000 3/8250 90.0 0.292 0.292 0.533 0.000 -0.000 -0.000 3/8250 100.0 0.262 0.262 0.480 0.000 -0.000 -0.000 3/8250 110.0 0.237 0.237 0.436 0.000 -0.000 -0.000 3/8250 120.0 0.217 0.217 0.400 0.000 -0.000 -0.000 3/8250 130.0 0.200 0.200 0.369 0.000 -0.000 -0.000 3/8250 140.0 0.186 0.186 0.342 0.000 -0.000 -0.000 3/8250 150.0 0.173 0.173 0.320 0.000 -0.000 -0.000 3/8250 160.0 0.162 0.162 0.300 0.000 -0.000 -0.000 3/8250 170.0 0.153 0.153 0.282 0.000 -0.000 -0.000 3/8250 180.0 0.144 0.144 0.267 0.000 -0.000 -0.000 3/8250 190.0 0.137 0.137 0.253 0.000 -0.000 -0.000 3/8250 200.0 0.130 0.130 0.240 0.000 -0.000 -0.000 3/8250 210.0 0.124 0.124 0.229 0.000 -0.000 -0.000 3/8250 220.0 0.118 0.118 0.218 0.000 -0.000 -0.000 3/8250 230.0 0.113 0.113 0.209 0.000 -0.000 -0.000 3/8250 240.0 0.108 0.108 0.200 0.000 -0.000 -0.000 3/8250 250.0 0.104 0.104 0.192 0.000 -0.000 -0.000 3/8250 260.0 0.100 0.100 0.185 0.000 -0.000 -0.000 3/8250 270.0 0.096 0.096 0.178 0.000 -0.000 -0.000 3/8250 280.0 0.093 0.093 0.172 0.000 -0.000 -0.000 3/8250 290.0 0.089 0.089 0.166 0.000 -0.000 -0.000 3/8250 300.0 0.086 0.086 0.160 0.000 -0.000 -0.000 3/8250 310.0 0.084 0.084 0.155 0.000 -0.000 -0.000 3/8250 320.0 0.081 0.081 0.150 0.000 -0.000 -0.000 3/8250 330.0 0.079 0.079 0.146 0.000 -0.000 -0.000 3/8250 340.0 0.076 0.076 0.141 0.000 -0.000 -0.000 3/8250 350.0 0.074 0.074 0.137 0.000 -0.000 -0.000 3/8250 360.0 0.072 0.072 0.133 0.000 -0.000 -0.000 3/8250 370.0 0.070 0.070 0.130 0.000 -0.000 -0.000 3/8250 380.0 0.068 0.068 0.126 0.000 -0.000 -0.000 3/8250 390.0 0.067 0.067 0.123 0.000 -0.000 -0.000 3/8250 400.0 0.065 0.065 0.120 0.000 -0.000 -0.000 3/8250 410.0 0.063 0.063 0.117 0.000 -0.000 -0.000 3/8250 420.0 0.062 0.062 0.114 0.000 -0.000 -0.000 3/8250 430.0 0.060 0.060 0.112 0.000 -0.000 -0.000 3/8250 440.0 0.059 0.059 0.109 0.000 -0.000 -0.000 3/8250 450.0 0.058 0.058 0.107 0.000 -0.000 -0.000 3/8250 460.0 0.056 0.056 0.104 0.000 -0.000 -0.000 3/8250 470.0 0.055 0.055 0.102 0.000 -0.000 -0.000 3/8250 480.0 0.054 0.054 0.100 0.000 -0.000 -0.000 3/8250 490.0 0.053 0.053 0.098 0.000 -0.000 -0.000 3/8250 500.0 0.052 0.052 0.096 0.000 -0.000 -0.000 3/8250 510.0 0.051 0.051 0.094 0.000 -0.000 -0.000 3/8250 520.0 0.050 0.050 0.092 0.000 -0.000 -0.000 3/8250 530.0 0.049 0.049 0.091 0.000 -0.000 -0.000 3/8250 540.0 0.048 0.048 0.089 0.000 -0.000 -0.000 3/8250 550.0 0.047 0.047 0.087 0.000 -0.000 -0.000 3/8250 560.0 0.046 0.046 0.086 0.000 -0.000 -0.000 3/8250 570.0 0.046 0.046 0.084 0.000 -0.000 -0.000 3/8250 580.0 0.045 0.045 0.083 0.000 -0.000 -0.000 3/8250 590.0 0.044 0.044 0.081 0.000 -0.000 -0.000 3/8250 600.0 0.043 0.043 0.080 0.000 -0.000 -0.000 3/8250 610.0 0.043 0.043 0.079 0.000 -0.000 -0.000 3/8250 620.0 0.042 0.042 0.078 0.000 -0.000 -0.000 3/8250 630.0 0.041 0.041 0.076 0.000 -0.000 -0.000 3/8250 640.0 0.041 0.041 0.075 0.000 -0.000 -0.000 3/8250 650.0 0.040 0.040 0.074 0.000 -0.000 -0.000 3/8250 660.0 0.039 0.039 0.073 0.000 -0.000 -0.000 3/8250 670.0 0.039 0.039 0.072 0.000 -0.000 -0.000 3/8250 680.0 0.038 0.038 0.071 0.000 -0.000 -0.000 3/8250 690.0 0.038 0.038 0.070 0.000 -0.000 -0.000 3/8250 700.0 0.037 0.037 0.069 0.000 -0.000 -0.000 3/8250 710.0 0.037 0.037 0.068 0.000 -0.000 -0.000 3/8250 720.0 0.036 0.036 0.067 0.000 -0.000 -0.000 3/8250 730.0 0.036 0.036 0.066 0.000 -0.000 -0.000 3/8250 740.0 0.035 0.035 0.065 0.000 -0.000 -0.000 3/8250 750.0 0.035 0.035 0.064 0.000 -0.000 -0.000 3/8250 760.0 0.034 0.034 0.063 0.000 -0.000 -0.000 3/8250 770.0 0.034 0.034 0.062 0.000 -0.000 -0.000 3/8250 780.0 0.033 0.033 0.062 0.000 -0.000 -0.000 3/8250 790.0 0.033 0.033 0.061 0.000 -0.000 -0.000 3/8250 800.0 0.032 0.032 0.060 0.000 -0.000 -0.000 3/8250 810.0 0.032 0.032 0.059 0.000 -0.000 -0.000 3/8250 820.0 0.032 0.032 0.059 0.000 -0.000 -0.000 3/8250 830.0 0.031 0.031 0.058 0.000 -0.000 -0.000 3/8250 840.0 0.031 0.031 0.057 0.000 -0.000 -0.000 3/8250 850.0 0.031 0.031 0.057 0.000 -0.000 -0.000 3/8250 860.0 0.030 0.030 0.056 0.000 -0.000 -0.000 3/8250 870.0 0.030 0.030 0.055 0.000 -0.000 -0.000 3/8250 880.0 0.029 0.029 0.055 0.000 -0.000 -0.000 3/8250 890.0 0.029 0.029 0.054 0.000 -0.000 -0.000 3/8250 900.0 0.029 0.029 0.053 0.000 -0.000 -0.000 3/8250 910.0 0.029 0.029 0.053 0.000 -0.000 -0.000 3/8250 920.0 0.028 0.028 0.052 0.000 -0.000 -0.000 3/8250 930.0 0.028 0.028 0.052 0.000 -0.000 -0.000 3/8250 940.0 0.028 0.028 0.051 0.000 -0.000 -0.000 3/8250 950.0 0.027 0.027 0.051 0.000 -0.000 -0.000 3/8250 960.0 0.027 0.027 0.050 0.000 -0.000 -0.000 3/8250 970.0 0.027 0.027 0.050 0.000 -0.000 -0.000 3/8250 980.0 0.026 0.026 0.049 0.000 -0.000 -0.000 3/8250 990.0 0.026 0.026 0.049 0.000 -0.000 -0.000 3/8250 1000.0 0.026 0.026 0.048 0.000 -0.000 -0.000 3/8250 Thermal conductivity related properties were written into "kappa-m555.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 07:21:22]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|