# Fileset

[LTC-calc.log](https://mdr.nims.go.jp/filesets/8ed3ffd9-d492-4633-9057-e16efbf94eba/download)

## Creator

[Atsushi Togo](https://orcid.org/0000-0001-8393-9766)

## Rights

Creative Commons Attribution 4.0 International[Creative Commons BY Attribution 4.0 International](https://creativecommons.org/licenses/by/4.0/)

## Other metadata

[First-principles lattice thermal conductivity calculation for Rb4PbBr6 / R-3c (167) / materials id 28564](https://mdr.nims.go.jp/datasets/d056e2b6-b857-43b5-9c1d-3e65170119b1)

## Fulltext

------------------------------------ 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.2Spglib version 2.6.1Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz".Unit of length: angstromSettings:  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.482577393333327Atomic 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.447732179999999Atomic 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.447732179999999Atomic 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/198Permutation basis: 8409/8409Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 592Number of blocks in projector: 300Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 3--- Eigsh_solver_block: 1 / 3 ---Block_size: 275Use standard eigh solver.--- Eigsh_solver_block: 2 / 3 ---Block_size: 234Use standard eigh solver.--- Eigsh_solver_block: 3 / 3 ---Block_size: 83Use 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 / 66Time (Solver_compr_matrix_reshape): 0.002Solver_block: 100 / 240 - Time: 0.299Solver_block: 200 / 240 - Time: 0.317Solver_block: 240 / 240 - Time: 0.145Solver: 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/198Permutation basis: 8409/8409Construct permutation basis matrix.Finding block diagonal structure in projector.Using scipy connected_components.Rank of projector: 592Number of blocks in projector: 300Finding block diagonal structure in projector.Using scipy connected_components.Number of blocks in projector (Sum rule): 3--- Eigsh_solver_block: 1 / 3 ---Block_size: 275Use standard eigh solver.--- Eigsh_solver_block: 2 / 3 ---Block_size: 234Use standard eigh solver.--- Eigsh_solver_block: 3 / 3 ---Block_size: 83Use standard eigh solver.Tree of FC basis block matrices:- (592, 573), data: False|-- (83, 79), data: True|-- (234, 228), data: True|-- (275, 266), data: TrueMax 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.2Spglib version 2.6.1----------------------------- General settings -----------------------------Run mode: force constantsHDF5 data compression filter: gzipCrystal 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.482577393333327Atomic 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.447732179999999Atomic 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): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperatures: 0.0  300.0 Cutoff frequency: 0.01Frequency 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.2Spglib version 2.6.1----------------------------- General settings -----------------------------Run mode: conductivity-RTAHDF5 data compression filter: gzipCrystal 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.482577393333327Atomic 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.447732179999999Atomic 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): TrueNAC unit conversion factor:  14.39965BZ integration: Tetrahedron-methodTemperature:  0.0  10.0  20.0  30.0  40.0 ... 1000.0Cutoff frequency: 0.01Frequency conversion factor to THz:  15.63330Length for sampling mesh generation: 50.00Generating 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.08Running 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.000Mass 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: 19Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 63Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 63Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 19Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 63Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 63Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 19Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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.000Mass 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: 39Calculating Gamma of ph-isotope with tetrahedron methodFrequency     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/8250Thermal conductivity related properties were written into "kappa-m555.hdf5".Summary of calculation was written in "phono3py.yaml".-------------------------[time 2026-01-08 07:21:22]-------------------------                 _   ___ _ __   __| |  / _ \ '_ \ / _` | |  __/ | | | (_| |  \___|_| |_|\__,_|