----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 23:51:04]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [1 1 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4_1/acd (142) Number of symmetry operations in supercell: 32 ------------------------------ primitive cell ------------------------------ Lattice vectors: a -5.809994520000000 5.809994520000000 9.553604125000000 b 5.809994520000000 -5.809994520000000 9.553604125000000 c 5.809994520000000 5.809994520000000 -9.553604125000000 Atomic positions (fractional): *1 Rb 0.12500000000000 0.23605311599251 0.61105311599251 85.468 2 Rb 0.98605311599251 0.87500000000000 0.61105311599251 85.468 3 Rb 0.26394688400749 0.37500000000000 0.38894688400749 85.468 4 Rb 0.62500000000000 0.51394688400749 0.38894688400749 85.468 5 Rb 0.62500000000000 0.01394688400749 0.88894688400749 85.468 6 Rb 0.76394688400749 0.37500000000000 0.88894688400749 85.468 7 Rb 0.48605311599251 0.87500000000000 0.11105311599251 85.468 8 Rb 0.12500000000000 0.73605311599251 0.11105311599251 85.468 *9 Rb 0.62880753841490 0.62880753841490 0.75761507682980 85.468 10 Rb 0.37880753841490 0.12119246158510 0.50000000000000 85.468 11 Rb 0.62119246158510 0.87880753841490 0.50000000000000 85.468 12 Rb 0.87119246158510 0.87119246158510 0.24238492317020 85.468 13 Rb 0.12119246158510 0.62119246158510 0.74238492317020 85.468 14 Rb 0.37119246158510 0.12880753841490 0.00000000000000 85.468 15 Rb 0.12880753841490 0.37119246158510 0.00000000000000 85.468 16 Rb 0.87880753841490 0.37880753841490 0.25761507682980 85.468 *17 Sn 0.94450260302027 0.11781192881177 0.95257901663650 118.710 18 Sn 0.86781192881177 0.41523291217527 0.67330932579150 118.710 19 Sn 0.74192358638377 0.19450260302027 0.32669067420850 118.710 20 Sn 0.30549739697973 0.75807641361623 0.67330932579150 118.710 21 Sn 0.08476708782473 0.63218807118823 0.32669067420850 118.710 22 Sn 0.50807641361623 0.33476708782473 0.95257901663650 118.710 23 Sn 0.38218807118823 0.55549739697973 0.04742098336350 118.710 24 Sn 0.16523291217527 0.99192358638377 0.04742098336350 118.710 25 Sn 0.80549739697973 0.13218807118823 0.54742098336350 118.710 26 Sn 0.88218807118823 0.83476708782473 0.82669067420850 118.710 27 Sn 0.00807641361623 0.05549739697973 0.17330932579150 118.710 28 Sn 0.44450260302027 0.49192358638377 0.82669067420850 118.710 29 Sn 0.66523291217527 0.61781192881177 0.17330932579150 118.710 30 Sn 0.24192358638377 0.91523291217527 0.54742098336350 118.710 31 Sn 0.36781192881177 0.69450260302027 0.45257901663650 118.710 32 Sn 0.58476708782473 0.25807641361623 0.45257901663650 118.710 -------------------------------- unit cell --------------------------------- Lattice vectors: a 11.619989040000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.619989040000000 0.000000000000000 c 0.000000000000000 0.000000000000000 19.107208249999999 Atomic positions (fractional): *1 Rb 0.86105311599251 0.75000000000000 0.37500000000000 85.468 > 1 2 Rb 0.75000000000000 0.86105311599251 0.12500000000000 85.468 > 2 3 Rb 0.25000000000000 0.13894688400749 0.12500000000000 85.468 > 3 4 Rb 0.13894688400749 0.25000000000000 0.37500000000000 85.468 > 4 5 Rb 0.63894688400749 0.25000000000000 0.37500000000000 85.468 > 5 6 Rb 0.25000000000000 0.63894688400749 0.12500000000000 85.468 > 6 7 Rb 0.75000000000000 0.36105311599251 0.12500000000000 85.468 > 7 8 Rb 0.36105311599251 0.75000000000000 0.37500000000000 85.468 > 8 *9 Rb 0.37880753841490 0.37880753841490 0.25000000000000 85.468 > 9 10 Rb 0.12119246158510 0.37880753841490 0.00000000000000 85.468 > 10 11 Rb 0.87880753841490 0.62119246158510 0.00000000000000 85.468 > 11 12 Rb 0.62119246158510 0.62119246158510 0.25000000000000 85.468 > 12 13 Rb 0.62119246158510 0.12119246158510 0.00000000000000 85.468 > 13 14 Rb 0.87880753841490 0.12119246158510 0.25000000000000 85.468 > 14 15 Rb 0.12119246158510 0.87880753841490 0.25000000000000 85.468 > 15 16 Rb 0.37880753841490 0.87880753841490 0.00000000000000 85.468 > 16 17 Rb 0.36105311599251 0.25000000000000 0.87500000000000 85.468 > 1 18 Rb 0.25000000000000 0.36105311599251 0.62500000000000 85.468 > 2 19 Rb 0.75000000000000 0.63894688400749 0.62500000000000 85.468 > 3 20 Rb 0.63894688400749 0.75000000000000 0.87500000000000 85.468 > 4 21 Rb 0.13894688400749 0.75000000000000 0.87500000000000 85.468 > 5 22 Rb 0.75000000000000 0.13894688400749 0.62500000000000 85.468 > 6 23 Rb 0.25000000000000 0.86105311599251 0.62500000000000 85.468 > 7 24 Rb 0.86105311599251 0.25000000000000 0.87500000000000 85.468 > 8 25 Rb 0.87880753841490 0.87880753841490 0.75000000000000 85.468 > 9 26 Rb 0.62119246158510 0.87880753841490 0.50000000000000 85.468 > 10 27 Rb 0.37880753841490 0.12119246158510 0.50000000000000 85.468 > 11 28 Rb 0.12119246158510 0.12119246158510 0.75000000000000 85.468 > 12 29 Rb 0.12119246158510 0.62119246158510 0.50000000000000 85.468 > 13 30 Rb 0.37880753841490 0.62119246158510 0.75000000000000 85.468 > 14 31 Rb 0.62119246158510 0.37880753841490 0.75000000000000 85.468 > 15 32 Rb 0.87880753841490 0.37880753841490 0.50000000000000 85.468 > 16 *33 Sn 0.06294417121400 0.88963484542250 0.05486775759777 118.710 > 17 34 Sn 0.11036515457750 0.56294417121400 0.30486775759777 118.710 > 18 35 Sn 0.88963484542250 0.43705582878600 0.30486775759777 118.710 > 19 36 Sn 0.56294417121400 0.11036515457750 0.19513224240223 118.710 > 20 37 Sn 0.43705582878600 0.88963484542250 0.19513224240223 118.710 > 21 38 Sn 0.88963484542250 0.06294417121400 0.44513224240223 118.710 > 22 39 Sn 0.11036515457750 0.93705582878600 0.44513224240223 118.710 > 23 40 Sn 0.93705582878600 0.11036515457750 0.05486775759777 118.710 > 24 41 Sn 0.93705582878600 0.61036515457750 0.19513224240223 118.710 > 25 42 Sn 0.38963484542250 0.43705582878600 0.44513224240223 118.710 > 26 43 Sn 0.61036515457750 0.56294417121400 0.44513224240223 118.710 > 27 44 Sn 0.43705582878600 0.38963484542250 0.05486775759777 118.710 > 28 45 Sn 0.56294417121400 0.61036515457750 0.05486775759777 118.710 > 29 46 Sn 0.61036515457750 0.93705582878600 0.30486775759777 118.710 > 30 47 Sn 0.38963484542250 0.06294417121400 0.30486775759777 118.710 > 31 48 Sn 0.06294417121400 0.38963484542250 0.19513224240223 118.710 > 32 49 Sn 0.56294417121400 0.38963484542250 0.55486775759777 118.710 > 17 50 Sn 0.61036515457750 0.06294417121400 0.80486775759777 118.710 > 18 51 Sn 0.38963484542250 0.93705582878600 0.80486775759777 118.710 > 19 52 Sn 0.06294417121400 0.61036515457750 0.69513224240223 118.710 > 20 53 Sn 0.93705582878600 0.38963484542250 0.69513224240223 118.710 > 21 54 Sn 0.38963484542250 0.56294417121400 0.94513224240223 118.710 > 22 55 Sn 0.61036515457750 0.43705582878600 0.94513224240223 118.710 > 23 56 Sn 0.43705582878600 0.61036515457750 0.55486775759777 118.710 > 24 57 Sn 0.43705582878600 0.11036515457750 0.69513224240223 118.710 > 25 58 Sn 0.88963484542250 0.93705582878600 0.94513224240223 118.710 > 26 59 Sn 0.11036515457750 0.06294417121400 0.94513224240223 118.710 > 27 60 Sn 0.93705582878600 0.88963484542250 0.55486775759777 118.710 > 28 61 Sn 0.06294417121400 0.11036515457750 0.55486775759777 118.710 > 29 62 Sn 0.11036515457750 0.43705582878600 0.80486775759777 118.710 > 30 63 Sn 0.88963484542250 0.56294417121400 0.80486775759777 118.710 > 31 64 Sn 0.56294417121400 0.88963484542250 0.69513224240223 118.710 > 32 -------------------------------- super cell -------------------------------- Lattice vectors: a 11.619989040000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.619989040000000 0.000000000000000 c 0.000000000000000 0.000000000000000 19.107208249999999 Atomic positions (fractional): *1 Rb 0.86105311599251 0.75000000000000 0.37500000000000 85.468 > 1 2 Rb 0.75000000000000 0.86105311599251 0.12500000000000 85.468 > 2 3 Rb 0.25000000000000 0.13894688400749 0.12500000000000 85.468 > 3 4 Rb 0.13894688400749 0.25000000000000 0.37500000000000 85.468 > 4 5 Rb 0.63894688400749 0.25000000000000 0.37500000000000 85.468 > 5 6 Rb 0.25000000000000 0.63894688400749 0.12500000000000 85.468 > 6 7 Rb 0.75000000000000 0.36105311599251 0.12500000000000 85.468 > 7 8 Rb 0.36105311599251 0.75000000000000 0.37500000000000 85.468 > 8 *9 Rb 0.37880753841490 0.37880753841490 0.25000000000000 85.468 > 9 10 Rb 0.12119246158510 0.37880753841490 0.00000000000000 85.468 > 10 11 Rb 0.87880753841490 0.62119246158510 0.00000000000000 85.468 > 11 12 Rb 0.62119246158510 0.62119246158510 0.25000000000000 85.468 > 12 13 Rb 0.62119246158510 0.12119246158510 0.00000000000000 85.468 > 13 14 Rb 0.87880753841490 0.12119246158510 0.25000000000000 85.468 > 14 15 Rb 0.12119246158510 0.87880753841490 0.25000000000000 85.468 > 15 16 Rb 0.37880753841490 0.87880753841490 0.00000000000000 85.468 > 16 17 Rb 0.36105311599251 0.25000000000000 0.87500000000000 85.468 > 1 18 Rb 0.25000000000000 0.36105311599251 0.62500000000000 85.468 > 2 19 Rb 0.75000000000000 0.63894688400749 0.62500000000000 85.468 > 3 20 Rb 0.63894688400749 0.75000000000000 0.87500000000000 85.468 > 4 21 Rb 0.13894688400749 0.75000000000000 0.87500000000000 85.468 > 5 22 Rb 0.75000000000000 0.13894688400749 0.62500000000000 85.468 > 6 23 Rb 0.25000000000000 0.86105311599251 0.62500000000000 85.468 > 7 24 Rb 0.86105311599251 0.25000000000000 0.87500000000000 85.468 > 8 25 Rb 0.87880753841490 0.87880753841490 0.75000000000000 85.468 > 9 26 Rb 0.62119246158510 0.87880753841490 0.50000000000000 85.468 > 10 27 Rb 0.37880753841490 0.12119246158510 0.50000000000000 85.468 > 11 28 Rb 0.12119246158510 0.12119246158510 0.75000000000000 85.468 > 12 29 Rb 0.12119246158510 0.62119246158510 0.50000000000000 85.468 > 13 30 Rb 0.37880753841490 0.62119246158510 0.75000000000000 85.468 > 14 31 Rb 0.62119246158510 0.37880753841490 0.75000000000000 85.468 > 15 32 Rb 0.87880753841490 0.37880753841490 0.50000000000000 85.468 > 16 *33 Sn 0.06294417121400 0.88963484542250 0.05486775759777 118.710 > 17 34 Sn 0.11036515457750 0.56294417121400 0.30486775759777 118.710 > 18 35 Sn 0.88963484542250 0.43705582878600 0.30486775759777 118.710 > 19 36 Sn 0.56294417121400 0.11036515457750 0.19513224240223 118.710 > 20 37 Sn 0.43705582878600 0.88963484542250 0.19513224240223 118.710 > 21 38 Sn 0.88963484542250 0.06294417121400 0.44513224240223 118.710 > 22 39 Sn 0.11036515457750 0.93705582878600 0.44513224240223 118.710 > 23 40 Sn 0.93705582878600 0.11036515457750 0.05486775759777 118.710 > 24 41 Sn 0.93705582878600 0.61036515457750 0.19513224240223 118.710 > 25 42 Sn 0.38963484542250 0.43705582878600 0.44513224240223 118.710 > 26 43 Sn 0.61036515457750 0.56294417121400 0.44513224240223 118.710 > 27 44 Sn 0.43705582878600 0.38963484542250 0.05486775759777 118.710 > 28 45 Sn 0.56294417121400 0.61036515457750 0.05486775759777 118.710 > 29 46 Sn 0.61036515457750 0.93705582878600 0.30486775759777 118.710 > 30 47 Sn 0.38963484542250 0.06294417121400 0.30486775759777 118.710 > 31 48 Sn 0.06294417121400 0.38963484542250 0.19513224240223 118.710 > 32 49 Sn 0.56294417121400 0.38963484542250 0.55486775759777 118.710 > 17 50 Sn 0.61036515457750 0.06294417121400 0.80486775759777 118.710 > 18 51 Sn 0.38963484542250 0.93705582878600 0.80486775759777 118.710 > 19 52 Sn 0.06294417121400 0.61036515457750 0.69513224240223 118.710 > 20 53 Sn 0.93705582878600 0.38963484542250 0.69513224240223 118.710 > 21 54 Sn 0.38963484542250 0.56294417121400 0.94513224240223 118.710 > 22 55 Sn 0.61036515457750 0.43705582878600 0.94513224240223 118.710 > 23 56 Sn 0.43705582878600 0.61036515457750 0.55486775759777 118.710 > 24 57 Sn 0.43705582878600 0.11036515457750 0.69513224240223 118.710 > 25 58 Sn 0.88963484542250 0.93705582878600 0.94513224240223 118.710 > 26 59 Sn 0.11036515457750 0.06294417121400 0.94513224240223 118.710 > 27 60 Sn 0.93705582878600 0.88963484542250 0.55486775759777 118.710 > 28 61 Sn 0.06294417121400 0.11036515457750 0.55486775759777 118.710 > 29 62 Sn 0.11036515457750 0.43705582878600 0.80486775759777 118.710 > 30 63 Sn 0.88963484542250 0.56294417121400 0.80486775759777 118.710 > 31 64 Sn 0.56294417121400 0.88963484542250 0.69513224240223 118.710 > 32 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 10.2020420 0.0000000 0.0000000 0.0000000 10.2020420 0.0000000 0.0000000 0.0000000 10.5977093 -------------------------- Born effective charges -------------------------- 1 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 0.3078549 0.0000000 0.4050185 0.9832835 2 Rb 1.3881946 0.0000000 -0.3078549 0.0000000 1.1852511 0.0000000 -0.4050185 0.0000000 0.9832835 3 Rb 1.3881946 0.0000000 0.3078549 0.0000000 1.1852511 0.0000000 0.4050185 0.0000000 0.9832835 4 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 -0.3078549 0.0000000 -0.4050185 0.9832835 5 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 0.3078549 0.0000000 0.4050185 0.9832835 6 Rb 1.3881946 0.0000000 -0.3078549 0.0000000 1.1852511 0.0000000 -0.4050185 0.0000000 0.9832835 7 Rb 1.3881946 0.0000000 0.3078549 0.0000000 1.1852511 0.0000000 0.4050185 0.0000000 0.9832835 8 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 -0.3078549 0.0000000 -0.4050185 0.9832835 9 Rb 1.1296195 -0.3130094 0.0875629 -0.3130094 1.1296195 -0.0875629 -0.0168349 0.0168349 1.3706846 10 Rb 1.1296195 0.3130094 0.0875629 0.3130094 1.1296195 0.0875629 -0.0168349 -0.0168349 1.3706846 11 Rb 1.1296195 0.3130094 -0.0875629 0.3130094 1.1296195 -0.0875629 0.0168349 0.0168349 1.3706846 12 Rb 1.1296195 -0.3130094 -0.0875629 -0.3130094 1.1296195 0.0875629 0.0168349 -0.0168349 1.3706846 13 Rb 1.1296195 -0.3130094 0.0875629 -0.3130094 1.1296195 -0.0875629 -0.0168349 0.0168349 1.3706846 14 Rb 1.1296195 0.3130094 0.0875629 0.3130094 1.1296195 0.0875629 -0.0168349 -0.0168349 1.3706846 15 Rb 1.1296195 0.3130094 -0.0875629 0.3130094 1.1296195 -0.0875629 0.0168349 0.0168349 1.3706846 16 Rb 1.1296195 -0.3130094 -0.0875629 -0.3130094 1.1296195 0.0875629 0.0168349 -0.0168349 1.3706846 17 Sn -1.1168476 0.6426186 -0.6713674 0.5173347 -1.2994948 0.3948036 -0.6006502 0.4134729 -1.1769840 18 Sn -1.2994948 -0.5173347 -0.3948036 -0.6426186 -1.1168476 -0.6713674 -0.4134729 -0.6006502 -1.1769840 19 Sn -1.2994948 -0.5173347 0.3948036 -0.6426186 -1.1168476 0.6713674 0.4134729 0.6006502 -1.1769840 20 Sn -1.1168476 -0.6426186 0.6713674 -0.5173347 -1.2994948 0.3948036 0.6006502 0.4134729 -1.1769840 21 Sn -1.1168476 -0.6426186 -0.6713674 -0.5173347 -1.2994948 -0.3948036 -0.6006502 -0.4134729 -1.1769840 22 Sn -1.2994948 0.5173347 -0.3948036 0.6426186 -1.1168476 0.6713674 -0.4134729 0.6006502 -1.1769840 23 Sn -1.2994948 0.5173347 0.3948036 0.6426186 -1.1168476 -0.6713674 0.4134729 -0.6006502 -1.1769840 24 Sn -1.1168476 0.6426186 0.6713674 0.5173347 -1.2994948 -0.3948036 0.6006502 -0.4134729 -1.1769840 25 Sn -1.1168476 0.6426186 -0.6713674 0.5173347 -1.2994948 0.3948036 -0.6006502 0.4134729 -1.1769840 26 Sn -1.2994948 -0.5173347 -0.3948036 -0.6426186 -1.1168476 -0.6713674 -0.4134729 -0.6006502 -1.1769840 27 Sn -1.2994948 -0.5173347 0.3948036 -0.6426186 -1.1168476 0.6713674 0.4134729 0.6006502 -1.1769840 28 Sn -1.1168476 -0.6426186 0.6713674 -0.5173347 -1.2994948 0.3948036 0.6006502 0.4134729 -1.1769840 29 Sn -1.1168476 -0.6426186 -0.6713674 -0.5173347 -1.2994948 -0.3948036 -0.6006502 -0.4134729 -1.1769840 30 Sn -1.2994948 0.5173347 -0.3948036 0.6426186 -1.1168476 0.6713674 -0.4134729 0.6006502 -1.1769840 31 Sn -1.2994948 0.5173347 0.3948036 0.6426186 -1.1168476 -0.6713674 0.4134729 -0.6006502 -1.1769840 32 Sn -1.1168476 0.6426186 0.6713674 0.5173347 -1.2994948 -0.3948036 0.6006502 -0.4134729 -1.1769840 ---------------------------------------------------------------------------- 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: 192/192 Permutation basis: 11472/11472 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 639 Number of blocks in projector: 639 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 4 --- Eigsh_solver_block: 1 / 4 --- Block_size: 256 Use standard eigh solver. --- Eigsh_solver_block: 2 / 4 --- Block_size: 164 Use standard eigh solver. --- Eigsh_solver_block: 3 / 4 --- Block_size: 132 Use standard eigh solver. --- Eigsh_solver_block: 4 / 4 --- Block_size: 87 Use standard eigh solver. Tree of FC basis block matrices: - (639, 620), data: False |-- (87, 84), data: True |-- (132, 129), data: True |-- (164, 159), data: True |-- (256, 248), data: True ----- Solver_atoms: 1 -- 64 / 64 Time (Solver_compr_matrix_reshape): 0.001 Solver_block: 100 / 240 - Time: 0.385 Solver_block: 200 / 240 - Time: 0.385 Solver_block: 240 / 240 - Time: 0.220 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 1.058 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 192/192 Permutation basis: 11472/11472 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 639 Number of blocks in projector: 639 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 4 --- Eigsh_solver_block: 1 / 4 --- Block_size: 256 Use standard eigh solver. --- Eigsh_solver_block: 2 / 4 --- Block_size: 164 Use standard eigh solver. --- Eigsh_solver_block: 3 / 4 --- Block_size: 132 Use standard eigh solver. --- Eigsh_solver_block: 4 / 4 --- Block_size: 87 Use standard eigh solver. Tree of FC basis block matrices: - (639, 620), data: False |-- (87, 84), data: True |-- (132, 129), data: True |-- (164, 159), data: True |-- (256, 248), data: True Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 23:51:14]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 23:51:15]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [1 1 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4_1/acd (142) ------------------------------ primitive cell ------------------------------ Lattice vectors: a -5.809994520000000 5.809994520000000 9.553604125000000 b 5.809994520000000 -5.809994520000000 9.553604125000000 c 5.809994520000000 5.809994520000000 -9.553604125000000 Atomic positions (fractional): 1 Rb 0.12500000000000 0.23605311599251 0.61105311599251 85.468 2 Rb 0.98605311599251 0.87500000000000 0.61105311599251 85.468 3 Rb 0.26394688400749 0.37500000000000 0.38894688400749 85.468 4 Rb 0.62500000000000 0.51394688400749 0.38894688400749 85.468 5 Rb 0.62500000000000 0.01394688400749 0.88894688400749 85.468 6 Rb 0.76394688400749 0.37500000000000 0.88894688400749 85.468 7 Rb 0.48605311599251 0.87500000000000 0.11105311599251 85.468 8 Rb 0.12500000000000 0.73605311599251 0.11105311599251 85.468 9 Rb 0.62880753841490 0.62880753841490 0.75761507682980 85.468 10 Rb 0.37880753841490 0.12119246158510 0.50000000000000 85.468 11 Rb 0.62119246158510 0.87880753841490 0.50000000000000 85.468 12 Rb 0.87119246158510 0.87119246158510 0.24238492317020 85.468 13 Rb 0.12119246158510 0.62119246158510 0.74238492317020 85.468 14 Rb 0.37119246158510 0.12880753841490 0.00000000000000 85.468 15 Rb 0.12880753841490 0.37119246158510 0.00000000000000 85.468 16 Rb 0.87880753841490 0.37880753841490 0.25761507682980 85.468 17 Sn 0.94450260302027 0.11781192881177 0.95257901663650 118.710 18 Sn 0.86781192881177 0.41523291217527 0.67330932579150 118.710 19 Sn 0.74192358638377 0.19450260302027 0.32669067420850 118.710 20 Sn 0.30549739697973 0.75807641361623 0.67330932579150 118.710 21 Sn 0.08476708782473 0.63218807118823 0.32669067420850 118.710 22 Sn 0.50807641361623 0.33476708782473 0.95257901663650 118.710 23 Sn 0.38218807118823 0.55549739697973 0.04742098336350 118.710 24 Sn 0.16523291217527 0.99192358638377 0.04742098336350 118.710 25 Sn 0.80549739697973 0.13218807118823 0.54742098336350 118.710 26 Sn 0.88218807118823 0.83476708782473 0.82669067420850 118.710 27 Sn 0.00807641361623 0.05549739697973 0.17330932579150 118.710 28 Sn 0.44450260302027 0.49192358638377 0.82669067420850 118.710 29 Sn 0.66523291217527 0.61781192881177 0.17330932579150 118.710 30 Sn 0.24192358638377 0.91523291217527 0.54742098336350 118.710 31 Sn 0.36781192881177 0.69450260302027 0.45257901663650 118.710 32 Sn 0.58476708782473 0.25807641361623 0.45257901663650 118.710 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.619989040000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.619989040000000 0.000000000000000 c 0.000000000000000 0.000000000000000 19.107208249999999 Atomic positions (fractional): 1 Rb 0.86105311599251 0.75000000000000 0.37500000000000 85.468 > 1 2 Rb 0.75000000000000 0.86105311599251 0.12500000000000 85.468 > 2 3 Rb 0.25000000000000 0.13894688400749 0.12500000000000 85.468 > 3 4 Rb 0.13894688400749 0.25000000000000 0.37500000000000 85.468 > 4 5 Rb 0.63894688400749 0.25000000000000 0.37500000000000 85.468 > 5 6 Rb 0.25000000000000 0.63894688400749 0.12500000000000 85.468 > 6 7 Rb 0.75000000000000 0.36105311599251 0.12500000000000 85.468 > 7 8 Rb 0.36105311599251 0.75000000000000 0.37500000000000 85.468 > 8 9 Rb 0.37880753841490 0.37880753841490 0.25000000000000 85.468 > 9 10 Rb 0.12119246158510 0.37880753841490 0.00000000000000 85.468 > 10 11 Rb 0.87880753841490 0.62119246158510 0.00000000000000 85.468 > 11 12 Rb 0.62119246158510 0.62119246158510 0.25000000000000 85.468 > 12 13 Rb 0.62119246158510 0.12119246158510 0.00000000000000 85.468 > 13 14 Rb 0.87880753841490 0.12119246158510 0.25000000000000 85.468 > 14 15 Rb 0.12119246158510 0.87880753841490 0.25000000000000 85.468 > 15 16 Rb 0.37880753841490 0.87880753841490 0.00000000000000 85.468 > 16 17 Rb 0.36105311599251 0.25000000000000 0.87500000000000 85.468 > 1 18 Rb 0.25000000000000 0.36105311599251 0.62500000000000 85.468 > 2 19 Rb 0.75000000000000 0.63894688400749 0.62500000000000 85.468 > 3 20 Rb 0.63894688400749 0.75000000000000 0.87500000000000 85.468 > 4 21 Rb 0.13894688400749 0.75000000000000 0.87500000000000 85.468 > 5 22 Rb 0.75000000000000 0.13894688400749 0.62500000000000 85.468 > 6 23 Rb 0.25000000000000 0.86105311599251 0.62500000000000 85.468 > 7 24 Rb 0.86105311599251 0.25000000000000 0.87500000000000 85.468 > 8 25 Rb 0.87880753841490 0.87880753841490 0.75000000000000 85.468 > 9 26 Rb 0.62119246158510 0.87880753841490 0.50000000000000 85.468 > 10 27 Rb 0.37880753841490 0.12119246158510 0.50000000000000 85.468 > 11 28 Rb 0.12119246158510 0.12119246158510 0.75000000000000 85.468 > 12 29 Rb 0.12119246158510 0.62119246158510 0.50000000000000 85.468 > 13 30 Rb 0.37880753841490 0.62119246158510 0.75000000000000 85.468 > 14 31 Rb 0.62119246158510 0.37880753841490 0.75000000000000 85.468 > 15 32 Rb 0.87880753841490 0.37880753841490 0.50000000000000 85.468 > 16 33 Sn 0.06294417121400 0.88963484542250 0.05486775759777 118.710 > 33 34 Sn 0.11036515457750 0.56294417121400 0.30486775759777 118.710 > 34 35 Sn 0.88963484542250 0.43705582878600 0.30486775759777 118.710 > 35 36 Sn 0.56294417121400 0.11036515457750 0.19513224240223 118.710 > 36 37 Sn 0.43705582878600 0.88963484542250 0.19513224240223 118.710 > 37 38 Sn 0.88963484542250 0.06294417121400 0.44513224240223 118.710 > 38 39 Sn 0.11036515457750 0.93705582878600 0.44513224240223 118.710 > 39 40 Sn 0.93705582878600 0.11036515457750 0.05486775759777 118.710 > 40 41 Sn 0.93705582878600 0.61036515457750 0.19513224240223 118.710 > 41 42 Sn 0.38963484542250 0.43705582878600 0.44513224240223 118.710 > 42 43 Sn 0.61036515457750 0.56294417121400 0.44513224240223 118.710 > 43 44 Sn 0.43705582878600 0.38963484542250 0.05486775759777 118.710 > 44 45 Sn 0.56294417121400 0.61036515457750 0.05486775759777 118.710 > 45 46 Sn 0.61036515457750 0.93705582878600 0.30486775759777 118.710 > 46 47 Sn 0.38963484542250 0.06294417121400 0.30486775759777 118.710 > 47 48 Sn 0.06294417121400 0.38963484542250 0.19513224240223 118.710 > 48 49 Sn 0.56294417121400 0.38963484542250 0.55486775759777 118.710 > 33 50 Sn 0.61036515457750 0.06294417121400 0.80486775759777 118.710 > 34 51 Sn 0.38963484542250 0.93705582878600 0.80486775759777 118.710 > 35 52 Sn 0.06294417121400 0.61036515457750 0.69513224240223 118.710 > 36 53 Sn 0.93705582878600 0.38963484542250 0.69513224240223 118.710 > 37 54 Sn 0.38963484542250 0.56294417121400 0.94513224240223 118.710 > 38 55 Sn 0.61036515457750 0.43705582878600 0.94513224240223 118.710 > 39 56 Sn 0.43705582878600 0.61036515457750 0.55486775759777 118.710 > 40 57 Sn 0.43705582878600 0.11036515457750 0.69513224240223 118.710 > 41 58 Sn 0.88963484542250 0.93705582878600 0.94513224240223 118.710 > 42 59 Sn 0.11036515457750 0.06294417121400 0.94513224240223 118.710 > 43 60 Sn 0.93705582878600 0.88963484542250 0.55486775759777 118.710 > 44 61 Sn 0.06294417121400 0.11036515457750 0.55486775759777 118.710 > 45 62 Sn 0.11036515457750 0.43705582878600 0.80486775759777 118.710 > 46 63 Sn 0.88963484542250 0.56294417121400 0.80486775759777 118.710 > 47 64 Sn 0.56294417121400 0.88963484542250 0.69513224240223 118.710 > 48 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 10.2020420 0.0000000 0.0000000 0.0000000 10.2020420 0.0000000 0.0000000 0.0000000 10.5977093 -------------------------- Born effective charges -------------------------- 1 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 0.3078549 0.0000000 0.4050185 0.9832835 2 Rb 1.3881946 0.0000000 -0.3078549 0.0000000 1.1852511 0.0000000 -0.4050185 0.0000000 0.9832835 3 Rb 1.3881946 0.0000000 0.3078549 0.0000000 1.1852511 0.0000000 0.4050185 0.0000000 0.9832835 4 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 -0.3078549 0.0000000 -0.4050185 0.9832835 5 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 0.3078549 0.0000000 0.4050185 0.9832835 6 Rb 1.3881946 0.0000000 -0.3078549 0.0000000 1.1852511 0.0000000 -0.4050185 0.0000000 0.9832835 7 Rb 1.3881946 0.0000000 0.3078549 0.0000000 1.1852511 0.0000000 0.4050185 0.0000000 0.9832835 8 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 -0.3078549 0.0000000 -0.4050185 0.9832835 9 Rb 1.1296195 -0.3130094 0.0875629 -0.3130094 1.1296195 -0.0875629 -0.0168349 0.0168349 1.3706846 10 Rb 1.1296195 0.3130094 0.0875629 0.3130094 1.1296195 0.0875629 -0.0168349 -0.0168349 1.3706846 11 Rb 1.1296195 0.3130094 -0.0875629 0.3130094 1.1296195 -0.0875629 0.0168349 0.0168349 1.3706846 12 Rb 1.1296195 -0.3130094 -0.0875629 -0.3130094 1.1296195 0.0875629 0.0168349 -0.0168349 1.3706846 13 Rb 1.1296195 -0.3130094 0.0875629 -0.3130094 1.1296195 -0.0875629 -0.0168349 0.0168349 1.3706846 14 Rb 1.1296195 0.3130094 0.0875629 0.3130094 1.1296195 0.0875629 -0.0168349 -0.0168349 1.3706846 15 Rb 1.1296195 0.3130094 -0.0875629 0.3130094 1.1296195 -0.0875629 0.0168349 0.0168349 1.3706846 16 Rb 1.1296195 -0.3130094 -0.0875629 -0.3130094 1.1296195 0.0875629 0.0168349 -0.0168349 1.3706846 17 Sn -1.1168476 0.6426186 -0.6713674 0.5173347 -1.2994948 0.3948036 -0.6006502 0.4134729 -1.1769840 18 Sn -1.2994948 -0.5173347 -0.3948036 -0.6426186 -1.1168476 -0.6713674 -0.4134729 -0.6006502 -1.1769840 19 Sn -1.2994948 -0.5173347 0.3948036 -0.6426186 -1.1168476 0.6713674 0.4134729 0.6006502 -1.1769840 20 Sn -1.1168476 -0.6426186 0.6713674 -0.5173347 -1.2994948 0.3948036 0.6006502 0.4134729 -1.1769840 21 Sn -1.1168476 -0.6426186 -0.6713674 -0.5173347 -1.2994948 -0.3948036 -0.6006502 -0.4134729 -1.1769840 22 Sn -1.2994948 0.5173347 -0.3948036 0.6426186 -1.1168476 0.6713674 -0.4134729 0.6006502 -1.1769840 23 Sn -1.2994948 0.5173347 0.3948036 0.6426186 -1.1168476 -0.6713674 0.4134729 -0.6006502 -1.1769840 24 Sn -1.1168476 0.6426186 0.6713674 0.5173347 -1.2994948 -0.3948036 0.6006502 -0.4134729 -1.1769840 25 Sn -1.1168476 0.6426186 -0.6713674 0.5173347 -1.2994948 0.3948036 -0.6006502 0.4134729 -1.1769840 26 Sn -1.2994948 -0.5173347 -0.3948036 -0.6426186 -1.1168476 -0.6713674 -0.4134729 -0.6006502 -1.1769840 27 Sn -1.2994948 -0.5173347 0.3948036 -0.6426186 -1.1168476 0.6713674 0.4134729 0.6006502 -1.1769840 28 Sn -1.1168476 -0.6426186 0.6713674 -0.5173347 -1.2994948 0.3948036 0.6006502 0.4134729 -1.1769840 29 Sn -1.1168476 -0.6426186 -0.6713674 -0.5173347 -1.2994948 -0.3948036 -0.6006502 -0.4134729 -1.1769840 30 Sn -1.2994948 0.5173347 -0.3948036 0.6426186 -1.1168476 0.6713674 -0.4134729 0.6006502 -1.1769840 31 Sn -1.2994948 0.5173347 0.3948036 0.6426186 -1.1168476 -0.6713674 0.4134729 -0.6006502 -1.1769840 32 Sn -1.1168476 0.6426186 0.6713674 0.5173347 -1.2994948 -0.3948036 0.6006502 -0.4134729 -1.1769840 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 9, 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[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000000 (yyy) 0.00000000 (yyy) 0.00000000 (yyy) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (xz) -0.00000000 (xz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 23:51:22]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 23:51:22]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [1 1 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4_1/acd (142) ------------------------------ primitive cell ------------------------------ Lattice vectors: a -5.809994520000000 5.809994520000000 9.553604125000000 b 5.809994520000000 -5.809994520000000 9.553604125000000 c 5.809994520000000 5.809994520000000 -9.553604125000000 Atomic positions (fractional): 1 Rb 0.12500000000000 0.23605311599251 0.61105311599251 85.468 2 Rb 0.98605311599251 0.87500000000000 0.61105311599251 85.468 3 Rb 0.26394688400749 0.37500000000000 0.38894688400749 85.468 4 Rb 0.62500000000000 0.51394688400749 0.38894688400749 85.468 5 Rb 0.62500000000000 0.01394688400749 0.88894688400749 85.468 6 Rb 0.76394688400749 0.37500000000000 0.88894688400749 85.468 7 Rb 0.48605311599251 0.87500000000000 0.11105311599251 85.468 8 Rb 0.12500000000000 0.73605311599251 0.11105311599251 85.468 9 Rb 0.62880753841490 0.62880753841490 0.75761507682980 85.468 10 Rb 0.37880753841490 0.12119246158510 0.50000000000000 85.468 11 Rb 0.62119246158510 0.87880753841490 0.50000000000000 85.468 12 Rb 0.87119246158510 0.87119246158510 0.24238492317020 85.468 13 Rb 0.12119246158510 0.62119246158510 0.74238492317020 85.468 14 Rb 0.37119246158510 0.12880753841490 0.00000000000000 85.468 15 Rb 0.12880753841490 0.37119246158510 0.00000000000000 85.468 16 Rb 0.87880753841490 0.37880753841490 0.25761507682980 85.468 17 Sn 0.94450260302027 0.11781192881177 0.95257901663650 118.710 18 Sn 0.86781192881177 0.41523291217527 0.67330932579150 118.710 19 Sn 0.74192358638377 0.19450260302027 0.32669067420850 118.710 20 Sn 0.30549739697973 0.75807641361623 0.67330932579150 118.710 21 Sn 0.08476708782473 0.63218807118823 0.32669067420850 118.710 22 Sn 0.50807641361623 0.33476708782473 0.95257901663650 118.710 23 Sn 0.38218807118823 0.55549739697973 0.04742098336350 118.710 24 Sn 0.16523291217527 0.99192358638377 0.04742098336350 118.710 25 Sn 0.80549739697973 0.13218807118823 0.54742098336350 118.710 26 Sn 0.88218807118823 0.83476708782473 0.82669067420850 118.710 27 Sn 0.00807641361623 0.05549739697973 0.17330932579150 118.710 28 Sn 0.44450260302027 0.49192358638377 0.82669067420850 118.710 29 Sn 0.66523291217527 0.61781192881177 0.17330932579150 118.710 30 Sn 0.24192358638377 0.91523291217527 0.54742098336350 118.710 31 Sn 0.36781192881177 0.69450260302027 0.45257901663650 118.710 32 Sn 0.58476708782473 0.25807641361623 0.45257901663650 118.710 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.619989040000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.619989040000000 0.000000000000000 c 0.000000000000000 0.000000000000000 19.107208249999999 Atomic positions (fractional): 1 Rb 0.86105311599251 0.75000000000000 0.37500000000000 85.468 > 1 2 Rb 0.75000000000000 0.86105311599251 0.12500000000000 85.468 > 2 3 Rb 0.25000000000000 0.13894688400749 0.12500000000000 85.468 > 3 4 Rb 0.13894688400749 0.25000000000000 0.37500000000000 85.468 > 4 5 Rb 0.63894688400749 0.25000000000000 0.37500000000000 85.468 > 5 6 Rb 0.25000000000000 0.63894688400749 0.12500000000000 85.468 > 6 7 Rb 0.75000000000000 0.36105311599251 0.12500000000000 85.468 > 7 8 Rb 0.36105311599251 0.75000000000000 0.37500000000000 85.468 > 8 9 Rb 0.37880753841490 0.37880753841490 0.25000000000000 85.468 > 9 10 Rb 0.12119246158510 0.37880753841490 0.00000000000000 85.468 > 10 11 Rb 0.87880753841490 0.62119246158510 0.00000000000000 85.468 > 11 12 Rb 0.62119246158510 0.62119246158510 0.25000000000000 85.468 > 12 13 Rb 0.62119246158510 0.12119246158510 0.00000000000000 85.468 > 13 14 Rb 0.87880753841490 0.12119246158510 0.25000000000000 85.468 > 14 15 Rb 0.12119246158510 0.87880753841490 0.25000000000000 85.468 > 15 16 Rb 0.37880753841490 0.87880753841490 0.00000000000000 85.468 > 16 17 Rb 0.36105311599251 0.25000000000000 0.87500000000000 85.468 > 1 18 Rb 0.25000000000000 0.36105311599251 0.62500000000000 85.468 > 2 19 Rb 0.75000000000000 0.63894688400749 0.62500000000000 85.468 > 3 20 Rb 0.63894688400749 0.75000000000000 0.87500000000000 85.468 > 4 21 Rb 0.13894688400749 0.75000000000000 0.87500000000000 85.468 > 5 22 Rb 0.75000000000000 0.13894688400749 0.62500000000000 85.468 > 6 23 Rb 0.25000000000000 0.86105311599251 0.62500000000000 85.468 > 7 24 Rb 0.86105311599251 0.25000000000000 0.87500000000000 85.468 > 8 25 Rb 0.87880753841490 0.87880753841490 0.75000000000000 85.468 > 9 26 Rb 0.62119246158510 0.87880753841490 0.50000000000000 85.468 > 10 27 Rb 0.37880753841490 0.12119246158510 0.50000000000000 85.468 > 11 28 Rb 0.12119246158510 0.12119246158510 0.75000000000000 85.468 > 12 29 Rb 0.12119246158510 0.62119246158510 0.50000000000000 85.468 > 13 30 Rb 0.37880753841490 0.62119246158510 0.75000000000000 85.468 > 14 31 Rb 0.62119246158510 0.37880753841490 0.75000000000000 85.468 > 15 32 Rb 0.87880753841490 0.37880753841490 0.50000000000000 85.468 > 16 33 Sn 0.06294417121400 0.88963484542250 0.05486775759777 118.710 > 33 34 Sn 0.11036515457750 0.56294417121400 0.30486775759777 118.710 > 34 35 Sn 0.88963484542250 0.43705582878600 0.30486775759777 118.710 > 35 36 Sn 0.56294417121400 0.11036515457750 0.19513224240223 118.710 > 36 37 Sn 0.43705582878600 0.88963484542250 0.19513224240223 118.710 > 37 38 Sn 0.88963484542250 0.06294417121400 0.44513224240223 118.710 > 38 39 Sn 0.11036515457750 0.93705582878600 0.44513224240223 118.710 > 39 40 Sn 0.93705582878600 0.11036515457750 0.05486775759777 118.710 > 40 41 Sn 0.93705582878600 0.61036515457750 0.19513224240223 118.710 > 41 42 Sn 0.38963484542250 0.43705582878600 0.44513224240223 118.710 > 42 43 Sn 0.61036515457750 0.56294417121400 0.44513224240223 118.710 > 43 44 Sn 0.43705582878600 0.38963484542250 0.05486775759777 118.710 > 44 45 Sn 0.56294417121400 0.61036515457750 0.05486775759777 118.710 > 45 46 Sn 0.61036515457750 0.93705582878600 0.30486775759777 118.710 > 46 47 Sn 0.38963484542250 0.06294417121400 0.30486775759777 118.710 > 47 48 Sn 0.06294417121400 0.38963484542250 0.19513224240223 118.710 > 48 49 Sn 0.56294417121400 0.38963484542250 0.55486775759777 118.710 > 33 50 Sn 0.61036515457750 0.06294417121400 0.80486775759777 118.710 > 34 51 Sn 0.38963484542250 0.93705582878600 0.80486775759777 118.710 > 35 52 Sn 0.06294417121400 0.61036515457750 0.69513224240223 118.710 > 36 53 Sn 0.93705582878600 0.38963484542250 0.69513224240223 118.710 > 37 54 Sn 0.38963484542250 0.56294417121400 0.94513224240223 118.710 > 38 55 Sn 0.61036515457750 0.43705582878600 0.94513224240223 118.710 > 39 56 Sn 0.43705582878600 0.61036515457750 0.55486775759777 118.710 > 40 57 Sn 0.43705582878600 0.11036515457750 0.69513224240223 118.710 > 41 58 Sn 0.88963484542250 0.93705582878600 0.94513224240223 118.710 > 42 59 Sn 0.11036515457750 0.06294417121400 0.94513224240223 118.710 > 43 60 Sn 0.93705582878600 0.88963484542250 0.55486775759777 118.710 > 44 61 Sn 0.06294417121400 0.11036515457750 0.55486775759777 118.710 > 45 62 Sn 0.11036515457750 0.43705582878600 0.80486775759777 118.710 > 46 63 Sn 0.88963484542250 0.56294417121400 0.80486775759777 118.710 > 47 64 Sn 0.56294417121400 0.88963484542250 0.69513224240223 118.710 > 48 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 10.2020420 0.0000000 0.0000000 0.0000000 10.2020420 0.0000000 0.0000000 0.0000000 10.5977093 -------------------------- Born effective charges -------------------------- 1 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 0.3078549 0.0000000 0.4050185 0.9832835 2 Rb 1.3881946 0.0000000 -0.3078549 0.0000000 1.1852511 0.0000000 -0.4050185 0.0000000 0.9832835 3 Rb 1.3881946 0.0000000 0.3078549 0.0000000 1.1852511 0.0000000 0.4050185 0.0000000 0.9832835 4 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 -0.3078549 0.0000000 -0.4050185 0.9832835 5 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 0.3078549 0.0000000 0.4050185 0.9832835 6 Rb 1.3881946 0.0000000 -0.3078549 0.0000000 1.1852511 0.0000000 -0.4050185 0.0000000 0.9832835 7 Rb 1.3881946 0.0000000 0.3078549 0.0000000 1.1852511 0.0000000 0.4050185 0.0000000 0.9832835 8 Rb 1.1852511 0.0000000 -0.0000000 0.0000000 1.3881946 -0.3078549 0.0000000 -0.4050185 0.9832835 9 Rb 1.1296195 -0.3130094 0.0875629 -0.3130094 1.1296195 -0.0875629 -0.0168349 0.0168349 1.3706846 10 Rb 1.1296195 0.3130094 0.0875629 0.3130094 1.1296195 0.0875629 -0.0168349 -0.0168349 1.3706846 11 Rb 1.1296195 0.3130094 -0.0875629 0.3130094 1.1296195 -0.0875629 0.0168349 0.0168349 1.3706846 12 Rb 1.1296195 -0.3130094 -0.0875629 -0.3130094 1.1296195 0.0875629 0.0168349 -0.0168349 1.3706846 13 Rb 1.1296195 -0.3130094 0.0875629 -0.3130094 1.1296195 -0.0875629 -0.0168349 0.0168349 1.3706846 14 Rb 1.1296195 0.3130094 0.0875629 0.3130094 1.1296195 0.0875629 -0.0168349 -0.0168349 1.3706846 15 Rb 1.1296195 0.3130094 -0.0875629 0.3130094 1.1296195 -0.0875629 0.0168349 0.0168349 1.3706846 16 Rb 1.1296195 -0.3130094 -0.0875629 -0.3130094 1.1296195 0.0875629 0.0168349 -0.0168349 1.3706846 17 Sn -1.1168476 0.6426186 -0.6713674 0.5173347 -1.2994948 0.3948036 -0.6006502 0.4134729 -1.1769840 18 Sn -1.2994948 -0.5173347 -0.3948036 -0.6426186 -1.1168476 -0.6713674 -0.4134729 -0.6006502 -1.1769840 19 Sn -1.2994948 -0.5173347 0.3948036 -0.6426186 -1.1168476 0.6713674 0.4134729 0.6006502 -1.1769840 20 Sn -1.1168476 -0.6426186 0.6713674 -0.5173347 -1.2994948 0.3948036 0.6006502 0.4134729 -1.1769840 21 Sn -1.1168476 -0.6426186 -0.6713674 -0.5173347 -1.2994948 -0.3948036 -0.6006502 -0.4134729 -1.1769840 22 Sn -1.2994948 0.5173347 -0.3948036 0.6426186 -1.1168476 0.6713674 -0.4134729 0.6006502 -1.1769840 23 Sn -1.2994948 0.5173347 0.3948036 0.6426186 -1.1168476 -0.6713674 0.4134729 -0.6006502 -1.1769840 24 Sn -1.1168476 0.6426186 0.6713674 0.5173347 -1.2994948 -0.3948036 0.6006502 -0.4134729 -1.1769840 25 Sn -1.1168476 0.6426186 -0.6713674 0.5173347 -1.2994948 0.3948036 -0.6006502 0.4134729 -1.1769840 26 Sn -1.2994948 -0.5173347 -0.3948036 -0.6426186 -1.1168476 -0.6713674 -0.4134729 -0.6006502 -1.1769840 27 Sn -1.2994948 -0.5173347 0.3948036 -0.6426186 -1.1168476 0.6713674 0.4134729 0.6006502 -1.1769840 28 Sn -1.1168476 -0.6426186 0.6713674 -0.5173347 -1.2994948 0.3948036 0.6006502 0.4134729 -1.1769840 29 Sn -1.1168476 -0.6426186 -0.6713674 -0.5173347 -1.2994948 -0.3948036 -0.6006502 -0.4134729 -1.1769840 30 Sn -1.2994948 0.5173347 -0.3948036 0.6426186 -1.1168476 0.6713674 -0.4134729 0.6006502 -1.1769840 31 Sn -1.2994948 0.5173347 0.3948036 0.6426186 -1.1168476 -0.6713674 0.4134729 -0.6006502 -1.1769840 32 Sn -1.1168476 0.6426186 0.6713674 0.5173347 -1.2994948 -0.3948036 0.6006502 -0.4134729 -1.1769840 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000000 (yyy) 0.00000000 (yyy) 0.00000000 (yyy) Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... Generalized regular grid: [ 1 4 24 ] Grid generation matrix: [ 0 4 4 ] [ 4 0 4 ] [ 3 3 0 ] 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.38, Number of G-points: 299, Lambda: 0.13 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/18) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 18 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 0.000 ( 0.000 0.000 0.000) 0.000 0.907 ( 0.000 0.000 0.000) 0.000 0.907 ( 0.000 0.000 -0.000) 0.000 0.920 ( 0.000 0.000 0.000) 0.000 0.920 ( 0.000 0.000 -0.000) 0.000 1.104 ( 0.000 0.000 0.000) 0.000 1.104 ( -0.000 -0.000 0.000) 0.000 1.121 ( 0.000 0.000 0.000) 0.000 1.121 ( 0.000 -0.000 -0.000) 0.000 1.130 ( 0.000 -0.000 0.000) 0.000 1.135 ( -0.000 0.000 -0.000) 0.000 1.135 ( 0.000 0.000 0.000) 0.000 1.145 ( 0.000 -0.000 0.000) 0.000 1.206 ( -0.000 0.000 0.000) 0.000 1.206 ( 0.000 0.000 -0.000) 0.000 1.221 ( 0.000 0.000 -0.000) 0.000 1.281 ( -0.000 -0.000 0.000) 0.000 1.284 ( 0.000 -0.000 0.000) 0.000 1.284 ( 0.000 0.000 0.000) 0.000 1.338 ( 0.000 -0.000 -0.000) 0.000 1.343 ( 0.000 0.000 0.000) 0.000 1.360 ( 0.000 0.000 0.000) 0.000 1.394 ( 0.000 0.000 0.000) 0.000 1.394 ( 0.000 0.000 0.000) 0.000 1.400 ( 0.000 0.000 -0.000) 0.000 1.406 ( 0.000 0.000 0.000) 0.000 1.406 ( 0.000 -0.000 -0.000) 0.000 1.410 ( -0.000 0.000 0.000) 0.000 1.428 ( 0.000 -0.000 0.000) 0.000 1.431 ( -0.000 0.000 0.000) 0.000 1.557 ( 0.000 -0.000 0.000) 0.000 1.557 ( 0.000 -0.000 0.000) 0.000 1.616 ( 0.000 -0.000 0.000) 0.000 1.626 ( 0.000 0.000 0.000) 0.000 1.636 ( 0.000 0.000 -0.000) 0.000 1.764 ( 0.000 0.000 -0.000) 0.000 1.764 ( 0.000 0.000 0.000) 0.000 1.804 ( 0.000 -0.000 -0.000) 0.000 1.815 ( 0.000 -0.000 0.000) 0.000 1.818 ( 0.000 0.000 0.000) 0.000 1.839 ( 0.000 0.000 0.000) 0.000 1.839 ( 0.000 0.000 -0.000) 0.000 1.890 ( -0.000 0.000 0.000) 0.000 1.913 ( 0.000 0.000 0.000) 0.000 1.913 ( -0.000 0.000 0.000) 0.000 1.920 ( -0.000 0.000 0.000) 0.000 1.930 ( 0.000 0.000 0.000) 0.000 1.975 ( -0.000 0.000 0.000) 0.000 1.999 ( -0.000 0.000 0.000) 0.000 1.999 ( 0.000 -0.000 0.000) 0.000 2.052 ( 0.000 -0.000 0.000) 0.000 2.067 ( 0.000 0.000 0.000) 0.000 2.067 ( -0.000 -0.000 0.000) 0.000 2.072 ( -0.000 0.000 0.000) 0.000 2.125 ( -0.000 0.000 0.000) 0.000 2.143 ( 0.000 0.000 0.000) 0.000 2.209 ( 0.000 -0.000 -0.000) 0.000 2.209 ( -0.000 0.000 0.000) 0.000 2.211 ( 0.000 -0.000 0.000) 0.000 2.254 ( 0.000 -0.000 0.000) 0.000 2.265 ( 0.000 0.000 0.000) 0.000 2.265 ( 0.000 -0.000 -0.000) 0.000 2.382 ( -0.000 0.000 -0.000) 0.000 2.382 ( 0.000 0.000 -0.000) 0.000 2.385 ( 0.000 0.000 0.000) 0.000 2.385 ( 0.000 0.000 -0.000) 0.000 2.416 ( 0.000 0.000 0.000) 0.000 2.434 ( 0.000 -0.000 0.000) 0.000 2.465 ( -0.000 0.000 0.000) 0.000 2.468 ( -0.000 0.000 0.000) 0.000 3.113 ( -0.000 -0.000 0.000) 0.000 3.167 ( -0.000 -0.000 0.000) 0.000 3.178 ( 0.000 -0.000 0.000) 0.000 3.188 ( 0.000 0.000 0.000) 0.000 3.269 ( 0.000 0.000 0.000) 0.000 3.314 ( -0.000 0.000 -0.000) 0.000 3.354 ( 0.000 0.000 0.000) 0.000 3.369 ( -0.000 0.000 0.000) 0.000 3.945 ( -0.000 -0.000 0.000) 0.000 3.945 ( 0.000 -0.000 0.000) 0.000 3.948 ( 0.000 0.000 0.000) 0.000 3.948 ( 0.000 -0.000 0.000) 0.000 3.991 ( 0.000 0.000 0.000) 0.000 3.991 ( 0.000 0.000 0.000) 0.000 4.034 ( -0.000 0.000 0.000) 0.000 4.034 ( 0.000 0.000 0.000) 0.000 4.079 ( 0.000 0.000 -0.000) 0.000 4.089 ( 0.000 0.000 0.000) 0.000 4.097 ( -0.000 0.000 0.000) 0.000 4.098 ( 0.000 0.000 0.000) 0.000 5.551 ( 0.000 -0.000 0.000) 0.000 5.569 ( 0.000 0.000 0.000) 0.000 5.584 ( 0.000 0.000 0.000) 0.000 5.587 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/18) ======================= q-point: ( 0.00 0.00 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.407 ( 9.044 9.044 0.000) 12.790 0.428 ( 9.104 9.104 0.000) 12.875 0.731 ( 15.280 15.280 0.000) 21.609 0.875 ( -2.107 -2.107 0.000) 2.979 0.886 ( -1.838 -1.838 0.000) 2.599 0.939 ( 1.281 1.281 0.000) 1.812 0.978 ( 1.791 1.791 0.000) 2.533 0.990 ( -5.075 -5.075 0.000) 7.176 1.002 ( -4.695 -4.695 0.000) 6.640 1.078 ( -1.290 -1.290 0.000) 1.825 1.082 ( -0.448 -0.448 0.000) 0.634 1.089 ( -1.493 -1.493 0.000) 2.111 1.119 ( -1.045 -1.045 0.000) 1.478 1.137 ( -1.735 -1.735 0.000) 2.454 1.137 ( -0.039 -0.039 0.000) 0.055 1.171 ( 0.800 0.800 0.000) 1.131 1.185 ( -0.814 -0.814 0.000) 1.151 1.221 ( -1.173 -1.173 0.000) 1.659 1.228 ( -1.322 -1.322 0.000) 1.869 1.229 ( -2.465 -2.465 0.000) 3.486 1.256 ( 1.203 1.203 0.000) 1.701 1.264 ( -1.627 -1.627 0.000) 2.300 1.333 ( 1.148 1.148 0.000) 1.624 1.353 ( 0.614 0.614 0.000) 0.869 1.424 ( 0.853 0.853 0.000) 1.207 1.437 ( 0.138 0.138 0.000) 0.195 1.449 ( 2.915 2.915 0.000) 4.123 1.454 ( 1.343 1.343 0.000) 1.899 1.463 ( 1.591 1.591 0.000) 2.251 1.464 ( 1.154 1.154 0.000) 1.632 1.477 ( 0.914 0.914 0.000) 1.292 1.517 ( 0.675 0.675 0.000) 0.955 1.568 ( 1.696 1.696 0.000) 2.398 1.580 ( -0.240 -0.240 0.000) 0.339 1.607 ( -0.508 -0.508 0.000) 0.718 1.627 ( -0.723 -0.723 0.000) 1.022 1.640 ( 0.512 0.512 0.000) 0.725 1.671 ( -4.907 -4.907 0.000) 6.939 1.729 ( -0.684 -0.684 0.000) 0.968 1.757 ( 0.546 0.546 0.000) 0.772 1.763 ( -2.165 -2.165 0.000) 3.062 1.824 ( 0.359 0.359 0.000) 0.507 1.852 ( -1.693 -1.693 0.000) 2.395 1.861 ( 0.621 0.621 0.000) 0.878 1.864 ( 1.751 1.751 0.000) 2.476 1.883 ( -0.893 -0.893 0.000) 1.262 1.912 ( -1.030 -1.030 0.000) 1.457 1.925 ( -0.185 -0.185 0.000) 0.262 1.949 ( 0.061 0.061 0.000) 0.086 1.958 ( -2.894 -2.894 0.000) 4.093 1.972 ( 1.032 1.032 0.000) 1.460 2.019 ( 0.503 0.503 0.000) 0.711 2.036 ( -0.598 -0.598 0.000) 0.845 2.043 ( 0.023 0.023 0.000) 0.032 2.104 ( -3.453 -3.453 0.000) 4.883 2.125 ( -1.081 -1.081 0.000) 1.529 2.126 ( -2.317 -2.317 0.000) 3.276 2.164 ( 1.107 1.107 0.000) 1.565 2.165 ( 3.017 3.017 0.000) 4.267 2.231 ( 1.750 1.750 0.000) 2.475 2.260 ( -0.799 -0.799 0.000) 1.129 2.278 ( 1.606 1.606 0.000) 2.272 2.279 ( -2.119 -2.119 0.000) 2.997 2.302 ( 1.366 1.366 0.000) 1.932 2.367 ( 0.394 0.394 0.000) 0.557 2.376 ( 1.079 1.079 0.000) 1.526 2.406 ( 0.912 0.912 0.000) 1.289 2.409 ( 1.049 1.049 0.000) 1.483 2.440 ( 0.424 0.424 0.000) 0.600 2.462 ( -0.191 -0.191 0.000) 0.269 2.476 ( 0.037 0.037 0.000) 0.052 2.563 ( -1.214 -1.214 0.000) 1.717 3.131 ( 0.676 0.676 0.000) 0.956 3.168 ( -0.047 -0.047 0.000) 0.067 3.179 ( 0.066 0.066 0.000) 0.094 3.190 ( 0.015 0.015 0.000) 0.021 3.285 ( 0.722 0.722 0.000) 1.021 3.318 ( 0.177 0.177 0.000) 0.250 3.345 ( -0.354 -0.354 0.000) 0.501 3.357 ( -0.552 -0.552 0.000) 0.780 3.948 ( 0.116 0.116 0.000) 0.165 3.951 ( 0.136 0.136 0.000) 0.192 3.959 ( 0.419 0.419 0.000) 0.593 3.978 ( -0.054 -0.054 0.000) 0.076 3.984 ( -0.219 -0.219 0.000) 0.310 3.993 ( 0.142 0.142 0.000) 0.201 4.016 ( -0.724 -0.724 0.000) 1.024 4.026 ( -0.406 -0.406 0.000) 0.574 4.088 ( -0.033 -0.033 0.000) 0.046 4.093 ( 0.418 0.418 0.000) 0.591 4.094 ( -0.137 -0.137 0.000) 0.193 4.102 ( 0.157 0.157 0.000) 0.221 5.558 ( 0.315 0.315 0.000) 0.446 5.569 ( 0.014 0.014 0.000) 0.020 5.581 ( -0.215 -0.215 0.000) 0.304 5.582 ( -0.116 -0.116 0.000) 0.164 ======================= Grid point 2 (3/18) ======================= q-point: ( 0.00 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 16 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.716 ( -0.769 -0.769 0.000) 1.088 0.716 ( 0.769 0.769 0.000) 1.088 0.768 ( -0.790 -0.790 0.000) 1.117 0.768 ( 0.790 0.790 0.000) 1.117 0.834 ( -0.455 -0.455 0.000) 0.644 0.834 ( 0.455 0.455 0.000) 0.644 1.007 ( -2.119 -2.119 0.000) 2.997 1.007 ( 2.119 2.119 0.000) 2.997 1.044 ( -1.001 -1.001 0.000) 1.416 1.044 ( 1.001 1.001 0.000) 1.416 1.061 ( -0.830 -0.830 0.000) 1.174 1.061 ( 0.830 0.830 0.000) 1.174 1.080 ( -1.048 -1.048 0.000) 1.482 1.080 ( 1.048 1.048 0.000) 1.482 1.185 ( -0.692 -0.692 0.000) 0.978 1.185 ( 0.692 0.692 0.000) 0.978 1.197 ( -0.440 -0.440 0.000) 0.622 1.197 ( 0.440 0.440 0.000) 0.622 1.215 ( -0.425 -0.425 0.000) 0.601 1.215 ( 0.425 0.425 0.000) 0.601 1.291 ( -4.598 -4.598 0.000) 6.502 1.291 ( 4.598 4.598 0.000) 6.502 1.335 ( -1.090 -1.090 0.000) 1.542 1.335 ( 1.090 1.090 0.000) 1.542 1.447 ( -0.019 -0.019 0.000) 0.027 1.447 ( 0.019 0.019 0.000) 0.027 1.470 ( -0.449 -0.449 0.000) 0.635 1.470 ( 0.449 0.449 0.000) 0.635 1.476 ( -0.502 -0.502 0.000) 0.710 1.476 ( 0.502 0.502 0.000) 0.710 1.524 ( -0.201 -0.201 0.000) 0.284 1.524 ( 0.201 0.201 0.000) 0.284 1.583 ( -2.436 -2.436 0.000) 3.445 1.583 ( 2.436 2.436 0.000) 3.445 1.599 ( -0.412 -0.412 0.000) 0.582 1.599 ( 0.412 0.412 0.000) 0.582 1.612 ( -1.720 -1.720 0.000) 2.433 1.612 ( 1.720 1.720 0.000) 2.433 1.721 ( -0.081 -0.081 0.000) 0.115 1.721 ( 0.081 0.081 0.000) 0.115 1.824 ( -2.249 -2.249 0.000) 3.180 1.824 ( 2.249 2.249 0.000) 3.180 1.831 ( -0.741 -0.741 0.000) 1.048 1.831 ( 0.741 0.741 0.000) 1.048 1.881 ( -0.318 -0.318 0.000) 0.450 1.881 ( 0.318 0.318 0.000) 0.450 1.897 ( -0.850 -0.850 0.000) 1.202 1.897 ( 0.850 0.850 0.000) 1.202 1.943 ( -1.768 -1.768 0.000) 2.500 1.943 ( 1.768 1.768 0.000) 2.500 1.966 ( -1.754 -1.754 0.000) 2.480 1.966 ( 1.754 1.754 0.000) 2.480 2.038 ( -0.143 -0.143 0.000) 0.202 2.038 ( 0.143 0.143 0.000) 0.202 2.069 ( -1.214 -1.214 0.000) 1.717 2.069 ( 1.214 1.214 0.000) 1.717 2.109 ( -3.207 -3.207 0.000) 4.535 2.109 ( 3.207 3.207 0.000) 4.535 2.210 ( -1.088 -1.088 0.000) 1.539 2.210 ( 1.088 1.088 0.000) 1.539 2.273 ( -2.036 -2.036 0.000) 2.879 2.273 ( 2.036 2.036 0.000) 2.879 2.319 ( -0.924 -0.924 0.000) 1.306 2.319 ( 0.924 0.924 0.000) 1.306 2.416 ( -0.352 -0.352 0.000) 0.497 2.416 ( 0.352 0.352 0.000) 0.497 2.437 ( -0.127 -0.127 0.000) 0.180 2.437 ( 0.127 0.127 0.000) 0.180 2.461 ( -0.392 -0.392 0.000) 0.554 2.461 ( 0.392 0.392 0.000) 0.554 2.488 ( -1.937 -1.937 0.000) 2.739 2.488 ( 1.937 1.937 0.000) 2.739 3.157 ( -0.465 -0.465 0.000) 0.657 3.157 ( 0.465 0.465 0.000) 0.657 3.185 ( -0.201 -0.201 0.000) 0.284 3.185 ( 0.201 0.201 0.000) 0.284 3.323 ( -0.930 -0.930 0.000) 1.315 3.323 ( 0.930 0.930 0.000) 1.315 3.329 ( -0.329 -0.329 0.000) 0.465 3.329 ( 0.329 0.329 0.000) 0.465 3.953 ( -0.091 -0.091 0.000) 0.129 3.953 ( 0.091 0.091 0.000) 0.129 3.973 ( -0.190 -0.190 0.000) 0.269 3.973 ( 0.190 0.190 0.000) 0.269 3.987 ( -0.455 -0.455 0.000) 0.643 3.987 ( 0.455 0.455 0.000) 0.643 4.007 ( -0.452 -0.452 0.000) 0.639 4.007 ( 0.452 0.452 0.000) 0.639 4.089 ( -0.088 -0.088 0.000) 0.125 4.089 ( 0.088 0.088 0.000) 0.125 4.105 ( -0.093 -0.093 0.000) 0.131 4.105 ( 0.093 0.093 0.000) 0.131 5.572 ( -0.159 -0.159 0.000) 0.225 5.572 ( 0.159 0.159 0.000) 0.225 5.573 ( -0.302 -0.302 0.000) 0.427 5.573 ( 0.302 0.302 0.000) 0.427 ======================= Grid point 5 (4/18) ======================= q-point: (-0.29 -0.04 0.04) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.409 ( 0.000 -11.969 -7.769) 14.269 0.453 ( 0.000 -10.297 -11.645) 15.545 0.609 ( 0.000 -16.143 -11.913) 20.063 0.844 ( -0.000 -1.151 4.596) 4.738 0.845 ( -0.000 -1.687 4.097) 4.431 0.963 ( -0.000 2.381 -0.287) 2.398 0.968 ( -0.000 1.083 -0.646) 1.260 0.988 ( -0.000 5.598 2.523) 6.140 1.009 ( -0.000 2.406 6.657) 7.079 1.046 ( -0.000 -0.197 4.949) 4.953 1.069 ( -0.000 1.654 3.189) 3.592 1.093 ( -0.000 1.031 3.438) 3.590 1.108 ( -0.000 3.027 -1.789) 3.516 1.134 ( 0.000 0.159 -0.220) 0.271 1.161 ( 0.000 1.268 -4.246) 4.431 1.198 ( -0.000 1.859 0.015) 1.859 1.218 ( 0.000 0.760 -1.825) 1.977 1.245 ( 0.000 -1.151 -0.695) 1.345 1.252 ( -0.000 1.897 0.608) 1.992 1.277 ( -0.000 1.109 -0.615) 1.268 1.287 ( -0.000 0.500 3.297) 3.334 1.301 ( 0.000 1.355 -2.544) 2.883 1.335 ( 0.000 0.376 -2.358) 2.388 1.335 ( 0.000 -1.426 0.213) 1.442 1.381 ( 0.000 -1.098 -0.887) 1.411 1.413 ( 0.000 1.285 -2.926) 3.196 1.419 ( 0.000 -0.887 0.263) 0.926 1.426 ( -0.000 -0.430 1.428) 1.492 1.449 ( 0.000 -2.437 -0.862) 2.585 1.469 ( 0.000 0.012 -1.387) 1.387 1.478 ( 0.000 -2.752 -1.059) 2.949 1.536 ( 0.000 -4.711 -1.782) 5.037 1.539 ( -0.000 1.113 1.678) 2.014 1.545 ( -0.000 1.626 0.088) 1.629 1.632 ( -0.000 1.114 0.705) 1.318 1.643 ( 0.000 -0.431 -2.100) 2.143 1.659 ( 0.000 -2.100 -0.340) 2.127 1.717 ( -0.000 1.039 1.960) 2.218 1.729 ( -0.000 0.363 2.338) 2.366 1.741 ( -0.000 2.010 0.857) 2.185 1.771 ( -0.000 2.559 2.557) 3.618 1.801 ( -0.000 3.662 0.762) 3.740 1.822 ( -0.000 -0.042 1.919) 1.919 1.830 ( -0.000 1.544 0.814) 1.745 1.841 ( 0.000 -2.702 1.425) 3.055 1.893 ( -0.000 1.685 0.655) 1.807 1.904 ( -0.000 0.034 0.472) 0.473 1.939 ( -0.000 -1.818 1.608) 2.427 1.944 ( -0.000 0.508 1.036) 1.154 1.968 ( -0.000 -2.756 2.257) 3.562 2.004 ( 0.000 1.582 -3.742) 4.063 2.033 ( -0.000 3.542 -0.369) 3.561 2.042 ( 0.000 -0.734 0.101) 0.741 2.043 ( 0.000 0.905 -0.857) 1.246 2.089 ( -0.000 -1.911 2.781) 3.374 2.134 ( -0.000 1.071 1.016) 1.476 2.135 ( 0.000 0.238 -4.321) 4.328 2.173 ( 0.000 1.097 -2.929) 3.128 2.188 ( 0.000 -1.459 1.127) 1.844 2.188 ( -0.000 0.742 0.177) 0.762 2.241 ( 0.000 -2.012 1.270) 2.379 2.287 ( 0.000 0.298 -2.051) 2.073 2.313 ( 0.000 0.294 -2.121) 2.141 2.331 ( 0.000 -2.212 -2.899) 3.646 2.350 ( -0.000 1.165 -0.298) 1.202 2.381 ( 0.000 1.153 -0.973) 1.509 2.386 ( 0.000 0.148 -0.350) 0.380 2.438 ( -0.000 -0.052 1.107) 1.108 2.442 ( 0.000 -0.588 -0.455) 0.743 2.461 ( 0.000 -0.819 0.552) 0.988 2.483 ( 0.000 -0.634 -0.581) 0.860 2.537 ( -0.000 -0.728 2.682) 2.779 3.136 ( 0.000 -0.655 -1.395) 1.541 3.175 ( 0.000 -0.058 -0.655) 0.658 3.178 ( -0.000 0.363 0.306) 0.475 3.185 ( -0.000 -0.065 0.505) 0.509 3.284 ( 0.000 -0.617 -0.799) 1.009 3.316 ( -0.000 -0.044 0.053) 0.069 3.346 ( -0.000 0.300 0.310) 0.432 3.359 ( -0.000 0.429 0.464) 0.632 3.945 ( 0.000 -0.424 0.070) 0.430 3.950 ( 0.000 -0.304 0.145) 0.337 3.963 ( 0.000 -0.228 -0.803) 0.835 3.976 ( 0.000 -0.837 -0.051) 0.838 3.986 ( -0.000 0.158 0.319) 0.356 3.987 ( -0.000 -0.008 0.179) 0.180 4.021 ( -0.000 0.660 0.602) 0.894 4.022 ( -0.000 0.502 0.604) 0.786 4.082 ( 0.000 -0.301 0.133) 0.329 4.091 ( 0.000 0.094 -0.303) 0.318 4.096 ( -0.000 -0.219 0.382) 0.441 4.102 ( -0.000 0.587 -0.420) 0.722 5.556 ( 0.000 -0.324 -0.170) 0.366 5.569 ( 0.000 0.048 -0.041) 0.063 5.581 ( -0.000 0.046 0.249) 0.253 5.582 ( -0.000 0.213 0.129) 0.250 ======================= Grid point 7 (5/18) ======================= q-point: (-0.29 -0.04 0.54) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.689 ( 9.955 2.227 -3.667) 10.840 0.706 ( 7.132 3.988 -8.438) 11.746 0.756 ( -8.124 0.075 3.421) 8.816 0.831 ( -2.095 -1.824 -0.251) 2.789 0.871 ( -2.173 -0.208 4.259) 4.785 0.884 ( -2.075 0.073 4.141) 4.632 0.933 ( -2.999 -1.301 -2.528) 4.132 0.976 ( -0.520 -1.067 -0.037) 1.187 0.991 ( 0.853 2.915 -0.713) 3.120 1.021 ( -1.712 0.862 0.767) 2.064 1.058 ( -0.670 1.146 -0.878) 1.592 1.062 ( 0.702 0.087 -0.173) 0.728 1.105 ( -0.922 0.791 1.173) 1.688 1.117 ( 2.908 1.310 1.705) 3.616 1.141 ( -2.325 -0.071 -2.992) 3.790 1.178 ( 5.586 0.440 -3.115) 6.411 1.193 ( 1.299 -1.874 0.355) 2.308 1.214 ( -0.503 -1.247 1.475) 1.996 1.233 ( -1.831 -0.579 1.392) 2.372 1.252 ( -0.749 1.029 2.825) 3.098 1.264 ( 0.706 -0.258 -1.451) 1.634 1.285 ( 0.760 -1.632 1.438) 2.304 1.324 ( 1.146 -0.517 -2.047) 2.402 1.347 ( -0.105 -0.485 -1.742) 1.812 1.387 ( -1.458 0.199 1.707) 2.253 1.418 ( -1.000 1.309 0.438) 1.704 1.425 ( 1.455 -0.970 3.418) 3.839 1.449 ( 1.129 0.701 3.092) 3.365 1.494 ( 0.159 1.279 -1.247) 1.793 1.510 ( -0.150 1.405 0.483) 1.493 1.535 ( 0.812 0.744 -1.353) 1.745 1.552 ( -1.320 0.538 -0.636) 1.560 1.563 ( 2.276 1.012 0.385) 2.521 1.578 ( -1.950 0.248 1.596) 2.532 1.616 ( -0.354 -0.588 -3.677) 3.741 1.625 ( 1.843 -1.893 -2.717) 3.789 1.637 ( -1.144 -1.111 -1.754) 2.370 1.665 ( -0.035 -1.413 -0.581) 1.528 1.691 ( -0.596 -0.503 2.433) 2.555 1.717 ( 0.947 -1.418 2.611) 3.119 1.761 ( -0.780 2.066 0.073) 2.209 1.778 ( -2.240 1.496 1.649) 3.159 1.809 ( 0.152 1.850 -1.301) 2.267 1.838 ( -0.339 0.508 1.847) 1.945 1.883 ( -0.361 1.584 -0.853) 1.834 1.894 ( 0.630 0.488 0.320) 0.859 1.904 ( -2.022 -1.076 -1.810) 2.919 1.923 ( 0.518 -1.173 -0.466) 1.364 1.943 ( -1.001 -0.931 -0.964) 1.673 1.946 ( -0.099 -1.068 0.214) 1.093 1.976 ( 0.354 -0.818 -0.466) 1.006 1.989 ( -1.219 1.831 0.565) 2.271 2.035 ( -0.842 -1.272 -0.472) 1.596 2.047 ( 0.401 -0.660 0.800) 1.112 2.075 ( 0.704 -2.440 -0.369) 2.566 2.092 ( -0.436 -1.224 1.318) 1.851 2.124 ( -0.268 -0.576 -3.885) 3.937 2.165 ( 1.614 0.578 -1.397) 2.212 2.181 ( -0.691 -0.074 1.771) 1.902 2.211 ( 0.099 -2.163 3.298) 3.945 2.267 ( -0.091 -0.397 -0.227) 0.466 2.277 ( 1.459 0.753 1.096) 1.974 2.325 ( 0.359 -0.856 -0.778) 1.212 2.341 ( 2.040 -0.860 0.401) 2.250 2.396 ( 0.638 2.136 -0.339) 2.255 2.404 ( 0.329 2.021 -0.639) 2.145 2.415 ( -0.551 2.246 -0.542) 2.376 2.440 ( 0.310 0.892 0.987) 1.366 2.456 ( 0.127 0.176 -1.137) 1.158 2.463 ( 0.713 -0.443 -0.382) 0.922 2.478 ( -0.048 -1.589 -0.674) 1.727 2.501 ( -2.337 -1.019 1.489) 2.952 3.156 ( 0.686 0.308 -0.873) 1.153 3.172 ( -0.113 -0.290 -0.587) 0.664 3.173 ( -0.246 0.179 0.857) 0.910 3.184 ( 0.045 -0.184 0.563) 0.594 3.308 ( 0.869 0.604 -0.353) 1.116 3.317 ( -0.361 0.401 0.546) 0.768 3.339 ( 0.122 -0.383 -0.201) 0.449 3.341 ( -0.628 -0.571 0.056) 0.851 3.954 ( -0.047 0.286 -0.425) 0.514 3.955 ( -0.157 0.291 -0.308) 0.452 3.972 ( 0.377 -0.035 -0.062) 0.383 3.977 ( 0.062 -0.385 0.415) 0.569 3.987 ( 0.071 -0.114 -0.291) 0.320 3.991 ( 0.134 0.250 0.242) 0.373 4.000 ( -0.192 -0.371 -0.180) 0.454 4.008 ( -0.468 -0.278 0.586) 0.800 4.089 ( -0.207 0.197 -0.295) 0.411 4.093 ( -0.012 0.097 -0.203) 0.226 4.096 ( 0.070 0.252 0.132) 0.293 4.102 ( 0.110 0.013 0.350) 0.367 5.567 ( 0.341 0.221 -0.052) 0.410 5.568 ( -0.161 0.165 0.077) 0.243 5.576 ( -0.135 -0.128 -0.090) 0.206 5.578 ( -0.103 -0.168 0.063) 0.207 ======================= Grid point 9 (6/18) ======================= q-point: ( 0.42 -0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 33 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.652 ( -0.000 8.371 9.561) 12.708 0.674 ( -0.000 11.677 4.118) 12.382 0.767 ( 0.000 -9.525 -5.231) 10.868 0.826 ( 0.000 -5.391 -8.174) 9.792 0.859 ( 0.000 -5.147 1.321) 5.314 0.892 ( -0.000 2.974 -2.249) 3.729 0.931 ( -0.000 -0.201 4.276) 4.281 0.975 ( -0.000 6.027 1.429) 6.194 0.981 ( 0.000 -4.576 1.405) 4.787 1.002 ( 0.000 -0.260 -3.178) 3.189 1.042 ( -0.000 -1.696 3.341) 3.747 1.076 ( 0.000 -4.058 -2.282) 4.656 1.080 ( -0.000 4.048 -0.608) 4.093 1.112 ( 0.000 -0.575 -3.286) 3.336 1.143 ( -0.000 0.470 3.295) 3.328 1.188 ( 0.000 0.328 -0.600) 0.684 1.198 ( -0.000 -0.516 1.262) 1.363 1.210 ( -0.000 1.468 -0.671) 1.614 1.257 ( 0.000 -1.430 0.760) 1.619 1.271 ( -0.000 1.263 0.667) 1.428 1.281 ( -0.000 -0.330 1.936) 1.964 1.305 ( 0.000 1.605 -4.515) 4.792 1.331 ( 0.000 -1.096 0.498) 1.203 1.358 ( -0.000 0.353 2.177) 2.205 1.373 ( 0.000 -1.348 -0.281) 1.377 1.412 ( -0.000 0.867 0.047) 0.868 1.414 ( 0.000 1.200 -2.381) 2.666 1.443 ( 0.000 1.227 -3.516) 3.724 1.474 ( -0.000 0.428 1.095) 1.176 1.494 ( -0.000 -0.046 0.121) 0.129 1.507 ( -0.000 1.080 -0.791) 1.338 1.526 ( 0.000 -0.274 -2.099) 2.117 1.557 ( -0.000 2.571 1.108) 2.800 1.579 ( 0.000 -2.995 -1.431) 3.320 1.624 ( -0.000 1.910 2.992) 3.549 1.653 ( 0.000 -3.145 0.788) 3.242 1.656 ( -0.000 1.454 3.189) 3.505 1.694 ( 0.000 -2.031 -0.664) 2.137 1.715 ( 0.000 -1.034 0.072) 1.037 1.719 ( -0.000 0.784 -0.508) 0.934 1.742 ( 0.000 -1.647 -2.657) 3.126 1.749 ( 0.000 1.570 -1.964) 2.514 1.799 ( -0.000 -1.541 1.537) 2.177 1.822 ( 0.000 -0.355 -1.119) 1.174 1.849 ( -0.000 -1.756 2.744) 3.257 1.889 ( 0.000 0.001 -0.809) 0.809 1.923 ( -0.000 2.227 1.802) 2.865 1.933 ( 0.000 -3.780 1.303) 3.998 1.956 ( -0.000 0.251 0.421) 0.490 1.958 ( 0.000 -2.487 -0.187) 2.494 1.966 ( -0.000 0.386 0.234) 0.451 2.000 ( -0.000 3.608 -2.445) 4.359 2.040 ( -0.000 0.173 0.692) 0.713 2.066 ( 0.000 0.476 -1.928) 1.986 2.094 ( -0.000 -3.299 3.473) 4.790 2.112 ( -0.000 -2.997 2.789) 4.094 2.125 ( -0.000 1.178 -0.130) 1.185 2.164 ( -0.000 3.529 -2.865) 4.546 2.182 ( -0.000 -0.213 1.436) 1.452 2.237 ( 0.000 1.796 -4.484) 4.830 2.269 ( 0.000 -1.318 0.552) 1.429 2.289 ( -0.000 2.046 -1.589) 2.590 2.341 ( -0.000 -0.784 0.740) 1.078 2.353 ( -0.000 2.963 -0.416) 2.992 2.362 ( -0.000 0.365 0.159) 0.398 2.370 ( 0.000 0.212 -0.399) 0.452 2.377 ( -0.000 0.868 2.499) 2.645 2.424 ( 0.000 -1.103 -2.437) 2.675 2.453 ( -0.000 0.298 1.177) 1.215 2.468 ( 0.000 -0.547 -0.230) 0.593 2.500 ( -0.000 0.633 1.345) 1.487 2.515 ( 0.000 -2.002 -0.140) 2.006 3.153 ( -0.000 0.838 0.949) 1.267 3.171 ( 0.000 -0.337 -0.969) 1.026 3.176 ( -0.000 0.006 0.554) 0.554 3.186 ( 0.000 0.033 -0.468) 0.469 3.301 ( -0.000 0.905 0.438) 1.006 3.312 ( 0.000 -0.389 -0.712) 0.811 3.343 ( -0.000 0.060 0.349) 0.354 3.347 ( 0.000 -0.570 -0.141) 0.587 3.950 ( -0.000 -0.027 0.344) 0.345 3.951 ( -0.000 -0.223 0.224) 0.317 3.971 ( -0.000 0.489 0.258) 0.553 3.981 ( 0.000 -0.115 -0.729) 0.738 3.988 ( -0.000 0.080 0.146) 0.167 3.990 ( -0.000 0.289 0.062) 0.295 4.005 ( 0.000 -0.498 0.110) 0.510 4.012 ( 0.000 -0.392 -0.466) 0.609 4.086 ( 0.000 -0.308 0.166) 0.350 4.090 ( -0.000 0.422 -0.188) 0.462 4.094 ( -0.000 -0.281 0.433) 0.516 4.101 ( 0.000 0.197 -0.341) 0.394 5.564 ( -0.000 0.391 0.081) 0.399 5.566 ( 0.000 -0.222 -0.118) 0.252 5.578 ( 0.000 -0.195 0.117) 0.227 5.579 ( 0.000 -0.056 -0.083) 0.100 ======================= Grid point 10 (7/18) ======================= q-point: (-0.58 -0.08 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 32 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.598 ( 5.169 -5.169 -6.659) 9.888 0.726 ( 7.479 -7.479 -6.490) 12.409 0.763 ( 1.032 -1.032 2.854) 3.206 0.813 ( 0.693 -0.693 2.963) 3.121 0.873 ( -2.615 2.615 5.413) 6.556 0.911 ( -0.518 0.518 -0.798) 1.083 0.917 ( -2.501 2.501 3.223) 4.785 0.937 ( 0.205 -0.205 3.611) 3.622 0.967 ( 2.458 -2.458 -6.695) 7.544 0.998 ( 0.417 -0.417 0.404) 0.715 0.999 ( 0.564 -0.564 0.245) 0.835 1.047 ( 2.023 -2.023 -0.494) 2.904 1.083 ( -0.143 0.143 5.085) 5.089 1.124 ( -0.862 0.862 0.793) 1.454 1.135 ( 0.054 -0.054 -1.314) 1.316 1.224 ( -1.568 1.568 0.822) 2.365 1.225 ( -1.113 1.113 -0.750) 1.744 1.230 ( -0.031 0.031 -1.062) 1.063 1.248 ( -1.078 1.078 -3.291) 3.627 1.255 ( -1.352 1.352 -0.445) 1.964 1.306 ( -0.517 0.517 1.343) 1.529 1.323 ( 0.846 -0.846 -0.026) 1.197 1.344 ( -1.548 1.548 -2.266) 3.151 1.358 ( -0.048 0.048 -0.610) 0.614 1.372 ( -2.406 2.406 0.298) 3.415 1.399 ( 0.944 -0.944 -0.026) 1.335 1.442 ( 0.324 -0.324 -3.068) 3.102 1.488 ( 1.102 -1.102 0.618) 1.677 1.497 ( -0.779 0.779 -4.133) 4.277 1.498 ( 0.667 -0.667 0.604) 1.120 1.523 ( -0.638 0.638 -2.682) 2.829 1.532 ( 1.305 -1.305 0.656) 1.959 1.540 ( 0.505 -0.505 0.525) 0.887 1.592 ( 0.649 -0.649 -1.604) 1.848 1.613 ( -0.154 0.154 3.583) 3.589 1.617 ( -0.524 0.524 0.086) 0.746 1.660 ( -0.822 0.822 2.735) 2.972 1.678 ( -1.558 1.558 2.303) 3.188 1.693 ( -0.819 0.819 -1.912) 2.235 1.707 ( 2.064 -2.064 -0.209) 2.927 1.745 ( -1.403 1.403 -1.712) 2.620 1.770 ( 0.262 -0.262 1.456) 1.503 1.793 ( 1.101 -1.101 2.261) 2.745 1.830 ( 1.803 -1.803 1.421) 2.919 1.854 ( 1.506 -1.506 -0.090) 2.131 1.870 ( 0.422 -0.422 0.555) 0.815 1.892 ( -0.354 0.354 1.476) 1.558 1.915 ( 0.925 -0.925 1.107) 1.713 1.936 ( 0.630 -0.630 0.620) 1.086 1.950 ( 0.368 -0.368 2.208) 2.268 1.956 ( -0.253 0.253 0.404) 0.539 1.986 ( -1.725 1.725 2.213) 3.293 2.040 ( -0.904 0.904 -1.508) 1.977 2.086 ( -1.858 1.858 0.764) 2.737 2.106 ( -0.567 0.567 0.287) 0.852 2.121 ( -1.494 1.494 0.255) 2.129 2.132 ( 0.561 -0.561 -1.152) 1.399 2.148 ( -1.553 1.553 -1.752) 2.809 2.203 ( 0.516 -0.516 0.769) 1.060 2.205 ( -0.132 0.132 1.621) 1.631 2.255 ( -0.324 0.324 -5.610) 5.629 2.318 ( -0.822 0.822 -1.887) 2.217 2.341 ( -1.344 1.344 -0.513) 1.968 2.371 ( 0.141 -0.141 -2.449) 2.458 2.393 ( 0.735 -0.735 -0.481) 1.146 2.405 ( 1.842 -1.842 -0.551) 2.662 2.412 ( -0.370 0.370 1.587) 1.671 2.422 ( 0.204 -0.204 -0.280) 0.402 2.437 ( 1.126 -1.126 0.178) 1.602 2.456 ( 0.839 -0.839 -0.805) 1.434 2.476 ( -0.348 0.348 1.064) 1.172 2.488 ( -0.593 0.593 2.245) 2.396 3.167 ( 0.274 -0.274 -1.389) 1.442 3.169 ( -0.233 0.233 0.475) 0.578 3.171 ( -0.163 0.163 0.667) 0.706 3.185 ( -0.066 0.066 -0.438) 0.447 3.307 ( 0.378 -0.378 -1.131) 1.251 3.319 ( 0.075 -0.075 -0.203) 0.229 3.334 ( -0.230 0.230 0.517) 0.611 3.344 ( -0.239 0.239 0.765) 0.836 3.947 ( 0.379 -0.379 0.304) 0.616 3.948 ( 0.324 -0.324 0.415) 0.618 3.975 ( 0.138 -0.138 -0.504) 0.541 3.979 ( -0.127 0.127 -0.071) 0.193 3.988 ( -0.117 0.117 -0.181) 0.245 3.990 ( 0.250 -0.250 -0.433) 0.559 4.004 ( -0.315 0.315 0.334) 0.557 4.007 ( -0.326 0.326 0.575) 0.737 4.083 ( 0.271 -0.271 0.147) 0.410 4.090 ( 0.149 -0.149 0.476) 0.520 4.098 ( 0.028 -0.028 -0.404) 0.406 4.104 ( -0.158 0.158 -0.134) 0.261 5.563 ( 0.231 -0.231 -0.274) 0.426 5.571 ( 0.002 -0.002 -0.138) 0.138 5.575 ( 0.013 -0.013 0.336) 0.337 5.578 ( -0.122 0.122 0.194) 0.259 ======================= Grid point 12 (8/18) ======================= q-point: ( 0.42 -0.08 -0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.484 ( -7.690 7.690 5.705) 12.281 0.547 ( -6.788 6.788 12.305) 15.606 0.782 ( -12.900 12.900 6.096) 19.235 0.838 ( -0.232 0.232 -4.106) 4.119 0.860 ( 0.307 -0.307 -1.859) 1.909 0.939 ( 3.411 -3.411 -1.177) 4.966 0.958 ( 1.782 -1.782 -2.114) 3.289 0.985 ( 1.181 -1.181 -1.552) 2.280 0.996 ( 2.368 -2.368 -0.558) 3.394 1.025 ( 0.622 -0.622 -4.968) 5.045 1.028 ( -0.272 0.272 -4.262) 4.279 1.087 ( 0.546 -0.546 -2.080) 2.218 1.112 ( 0.927 -0.927 -0.702) 1.488 1.117 ( 0.390 -0.390 -0.565) 0.790 1.148 ( 1.080 -1.080 4.194) 4.464 1.185 ( 1.259 -1.259 -0.428) 1.832 1.197 ( 1.442 -1.442 2.814) 3.476 1.228 ( 0.667 -0.667 -0.483) 1.060 1.240 ( 1.554 -1.554 1.176) 2.493 1.271 ( -0.108 0.108 -4.774) 4.777 1.285 ( 0.797 -0.797 2.539) 2.778 1.287 ( 0.280 -0.280 2.917) 2.943 1.309 ( 0.284 -0.284 3.221) 3.246 1.344 ( -0.219 0.219 0.096) 0.324 1.402 ( 0.427 -0.427 -0.948) 1.124 1.412 ( -1.895 1.895 -1.232) 2.950 1.420 ( 0.610 -0.610 -2.240) 2.401 1.447 ( -0.475 0.475 0.719) 0.984 1.476 ( -1.956 1.956 2.044) 3.440 1.494 ( -1.457 1.457 0.178) 2.069 1.508 ( -1.631 1.631 1.808) 2.931 1.542 ( 0.852 -0.852 0.172) 1.217 1.553 ( -1.421 1.421 -1.889) 2.758 1.558 ( -2.400 2.400 1.236) 3.612 1.616 ( 1.328 -1.328 -0.161) 1.885 1.645 ( -0.063 0.063 3.127) 3.128 1.654 ( -0.268 0.268 0.091) 0.390 1.674 ( 3.515 -3.515 0.291) 4.979 1.705 ( 0.667 -0.667 -2.312) 2.497 1.715 ( 0.682 -0.682 -1.865) 2.100 1.784 ( 1.500 -1.500 -1.877) 2.833 1.804 ( -0.052 0.052 -2.079) 2.080 1.815 ( 0.235 -0.235 -1.319) 1.360 1.845 ( -1.364 1.364 -0.895) 2.126 1.862 ( -0.547 0.547 -0.520) 0.932 1.877 ( 0.323 -0.323 -0.282) 0.537 1.915 ( -0.216 0.216 -0.969) 1.016 1.920 ( 0.208 -0.208 0.424) 0.516 1.945 ( 0.107 -0.107 -0.356) 0.386 1.964 ( -0.131 0.131 -0.548) 0.579 1.999 ( 0.543 -0.543 0.268) 0.814 2.002 ( 1.348 -1.348 -1.113) 2.208 2.030 ( 0.751 -0.751 -0.383) 1.130 2.088 ( -0.008 0.008 2.816) 2.816 2.094 ( 2.711 -2.711 -0.441) 3.859 2.123 ( 1.363 -1.363 -1.005) 2.174 2.127 ( -0.563 0.563 -0.253) 0.835 2.140 ( -1.626 1.626 -0.748) 2.418 2.158 ( 0.869 -0.869 2.747) 3.009 2.225 ( -1.530 1.530 -0.616) 2.250 2.240 ( 0.515 -0.515 -2.280) 2.394 2.282 ( 0.498 -0.498 1.993) 2.114 2.313 ( 0.065 -0.065 2.434) 2.436 2.341 ( 0.676 -0.676 1.670) 1.925 2.372 ( -1.856 1.856 1.356) 2.955 2.380 ( -0.719 0.719 0.607) 1.184 2.418 ( -0.983 0.983 0.521) 1.484 2.424 ( -0.842 0.842 1.230) 1.712 2.442 ( -0.367 0.367 0.492) 0.715 2.449 ( 0.294 -0.294 -1.372) 1.434 2.481 ( 0.096 -0.096 0.322) 0.350 2.537 ( 0.776 -0.776 -2.706) 2.920 3.142 ( -0.551 0.551 1.222) 1.449 3.174 ( 0.064 -0.064 0.638) 0.644 3.177 ( 0.101 -0.101 -0.335) 0.365 3.184 ( 0.041 -0.041 -0.679) 0.681 3.291 ( -0.636 0.636 0.665) 1.118 3.318 ( -0.156 0.156 -0.024) 0.222 3.342 ( 0.325 -0.325 -0.327) 0.564 3.353 ( 0.467 -0.467 -0.361) 0.753 3.949 ( -0.347 0.347 0.111) 0.503 3.953 ( -0.216 0.216 0.001) 0.305 3.965 ( -0.196 0.196 0.598) 0.659 3.980 ( 0.005 -0.005 0.269) 0.269 3.982 ( -0.006 0.006 -0.102) 0.102 3.989 ( -0.085 0.085 -0.327) 0.348 4.012 ( 0.633 -0.633 -0.429) 0.992 4.019 ( 0.395 -0.395 -0.662) 0.866 4.087 ( -0.376 0.376 -0.317) 0.619 4.091 ( 0.009 -0.009 0.340) 0.340 4.098 ( -0.155 0.155 -0.377) 0.436 4.100 ( 0.217 -0.217 0.339) 0.458 5.559 ( -0.295 0.295 0.144) 0.441 5.569 ( -0.015 0.015 0.053) 0.057 5.580 ( 0.080 -0.080 -0.207) 0.236 5.580 ( 0.188 -0.188 -0.107) 0.287 ======================= Grid point 13 (9/18) ======================= q-point: ( 0.13 -0.12 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 20 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.289 ( -0.000 13.066 0.000) 13.066 0.331 ( -0.000 14.956 0.000) 14.956 0.515 ( -0.000 22.767 0.000) 22.767 0.897 ( 0.000 -0.955 -0.000) 0.955 0.912 ( 0.000 -0.620 -0.000) 0.620 0.919 ( -0.000 0.136 0.000) 0.136 0.949 ( -0.000 1.945 0.000) 1.945 1.053 ( 0.000 -5.147 -0.000) 5.147 1.057 ( 0.000 -5.255 -0.000) 5.255 1.086 ( 0.000 -4.700 -0.000) 4.700 1.101 ( 0.000 -1.860 -0.000) 1.860 1.114 ( 0.000 -2.167 -0.000) 2.167 1.126 ( 0.000 -0.238 -0.000) 0.238 1.126 ( 0.000 -0.844 -0.000) 0.844 1.130 ( 0.000 -0.982 -0.000) 0.982 1.179 ( 0.000 -0.100 -0.000) 0.100 1.200 ( 0.000 -0.484 -0.000) 0.484 1.229 ( 0.000 -3.488 -0.000) 3.488 1.242 ( -0.000 1.165 0.000) 1.165 1.256 ( 0.000 -0.281 -0.000) 0.281 1.259 ( 0.000 -3.528 -0.000) 3.528 1.323 ( -0.000 1.162 0.000) 1.162 1.337 ( -0.000 0.522 0.000) 0.522 1.353 ( -0.000 0.334 0.000) 0.334 1.380 ( 0.000 -0.565 -0.000) 0.565 1.385 ( -0.000 0.795 0.000) 0.795 1.418 ( -0.000 1.484 0.000) 1.484 1.424 ( -0.000 3.112 0.000) 3.112 1.441 ( -0.000 1.132 0.000) 1.132 1.457 ( -0.000 0.750 0.000) 0.750 1.482 ( -0.000 3.213 0.000) 3.213 1.519 ( -0.000 5.034 0.000) 5.034 1.556 ( 0.000 -0.271 -0.000) 0.271 1.561 ( -0.000 0.204 0.000) 0.204 1.614 ( 0.000 -0.230 -0.000) 0.230 1.640 ( 0.000 -0.368 -0.000) 0.368 1.646 ( -0.000 0.737 0.000) 0.737 1.743 ( 0.000 -0.384 -0.000) 0.384 1.745 ( 0.000 -4.014 -0.000) 4.014 1.758 ( 0.000 -0.352 -0.000) 0.352 1.799 ( 0.000 -3.947 -0.000) 3.947 1.810 ( 0.000 -1.041 -0.000) 1.041 1.836 ( 0.000 -3.540 -0.000) 3.540 1.846 ( -0.000 2.115 0.000) 2.115 1.856 ( 0.000 -1.287 -0.000) 1.287 1.887 ( 0.000 -1.318 -0.000) 1.318 1.906 ( 0.000 -4.320 -0.000) 4.320 1.948 ( -0.000 2.035 0.000) 2.035 1.954 ( -0.000 2.207 0.000) 2.207 1.980 ( 0.000 -1.140 -0.000) 1.140 1.996 ( 0.000 -4.272 -0.000) 4.272 2.010 ( -0.000 0.021 0.000) 0.021 2.041 ( 0.000 -1.642 -0.000) 1.642 2.048 ( -0.000 1.736 0.000) 1.736 2.078 ( 0.000 -0.250 -0.000) 0.250 2.121 ( -0.000 3.276 0.000) 3.276 2.141 ( 0.000 -0.316 -0.000) 0.316 2.148 ( 0.000 -0.970 -0.000) 0.970 2.192 ( 0.000 -1.153 -0.000) 1.153 2.203 ( -0.000 2.601 0.000) 2.601 2.261 ( -0.000 4.084 0.000) 4.084 2.267 ( -0.000 0.091 0.000) 0.091 2.289 ( -0.000 2.004 0.000) 2.004 2.291 ( -0.000 1.740 0.000) 1.740 2.347 ( 0.000 -1.627 -0.000) 1.627 2.371 ( 0.000 -0.690 -0.000) 0.690 2.380 ( 0.000 -0.191 -0.000) 0.191 2.436 ( -0.000 0.185 0.000) 0.185 2.450 ( -0.000 1.392 0.000) 1.392 2.464 ( 0.000 -0.057 -0.000) 0.057 2.476 ( -0.000 0.593 0.000) 0.593 2.567 ( 0.000 -1.799 -0.000) 1.799 3.122 ( -0.000 0.817 0.000) 0.817 3.168 ( -0.000 0.101 0.000) 0.101 3.179 ( -0.000 0.057 0.000) 0.057 3.190 ( -0.000 0.087 0.000) 0.087 3.277 ( -0.000 0.739 0.000) 0.739 3.316 ( -0.000 0.170 0.000) 0.170 3.349 ( 0.000 -0.428 -0.000) 0.428 3.363 ( 0.000 -0.520 -0.000) 0.520 3.946 ( -0.000 0.069 0.000) 0.069 3.949 ( -0.000 0.105 0.000) 0.105 3.954 ( -0.000 0.582 0.000) 0.582 3.981 ( -0.000 0.245 0.000) 0.245 3.989 ( 0.000 -0.138 -0.000) 0.138 3.990 ( 0.000 -0.112 -0.000) 0.112 4.028 ( 0.000 -0.491 -0.000) 0.491 4.029 ( 0.000 -0.833 -0.000) 0.833 4.083 ( -0.000 0.318 0.000) 0.318 4.088 ( 0.000 -0.087 -0.000) 0.087 4.095 ( 0.000 -0.210 -0.000) 0.210 4.100 ( -0.000 0.210 0.000) 0.210 5.555 ( -0.000 0.343 0.000) 0.343 5.569 ( 0.000 -0.032 -0.000) 0.032 5.583 ( 0.000 -0.090 -0.000) 0.090 5.584 ( 0.000 -0.262 -0.000) 0.262 ======================= Grid point 14 (10/18) ======================= q-point: ( 0.13 -0.12 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 32 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.614 ( 5.156 9.237 0.000) 10.579 0.658 ( 2.469 10.907 0.000) 11.183 0.793 ( -1.805 -5.235 -0.000) 5.537 0.815 ( -2.216 -5.383 -0.000) 5.821 0.893 ( -1.746 -3.674 -0.000) 4.068 0.906 ( -4.925 -1.118 0.000) 5.050 0.957 ( 2.150 -0.744 -0.000) 2.275 0.980 ( 3.429 -3.400 -0.000) 4.829 0.982 ( 4.537 1.990 -0.000) 4.954 1.014 ( 2.835 -4.088 -0.000) 4.975 1.043 ( -1.523 -2.310 -0.000) 2.766 1.059 ( 1.348 -4.137 -0.000) 4.351 1.112 ( -2.285 2.309 0.000) 3.249 1.120 ( -3.161 -0.135 0.000) 3.164 1.141 ( 3.616 8.203 0.000) 8.965 1.177 ( 1.474 0.524 -0.000) 1.565 1.193 ( 0.748 -0.293 -0.000) 0.803 1.217 ( -0.973 3.145 0.000) 3.292 1.229 ( -1.107 0.986 0.000) 1.483 1.239 ( -1.499 2.002 0.000) 2.501 1.271 ( -3.862 1.161 0.000) 4.033 1.274 ( 1.879 -0.419 -0.000) 1.925 1.333 ( 0.473 -2.087 -0.000) 2.140 1.353 ( -0.543 0.936 0.000) 1.082 1.418 ( 2.164 -1.208 -0.000) 2.478 1.428 ( 1.186 -1.579 -0.000) 1.975 1.467 ( 0.223 0.800 0.000) 0.831 1.476 ( -0.203 2.144 0.000) 2.154 1.494 ( -0.164 -0.321 -0.000) 0.360 1.506 ( -0.719 0.513 0.000) 0.883 1.514 ( 2.219 2.614 0.000) 3.429 1.547 ( 0.707 1.433 0.000) 1.598 1.564 ( -2.365 -2.911 0.000) 3.751 1.564 ( -0.875 0.388 0.000) 0.957 1.600 ( 0.733 -1.950 -0.000) 2.083 1.602 ( -0.361 -0.203 0.000) 0.414 1.606 ( 0.008 0.747 0.000) 0.747 1.656 ( -2.039 0.015 0.000) 2.039 1.728 ( -0.921 0.229 0.000) 0.949 1.730 ( -0.975 -1.064 -0.000) 1.443 1.768 ( 3.027 -3.592 -0.000) 4.697 1.781 ( 1.395 1.708 0.000) 2.205 1.825 ( 1.953 -0.659 -0.000) 2.061 1.850 ( 2.586 -2.690 -0.000) 3.731 1.866 ( -0.066 -1.138 -0.000) 1.140 1.885 ( -0.881 -0.511 0.000) 1.018 1.900 ( -0.141 -0.585 -0.000) 0.602 1.912 ( -0.327 1.388 0.000) 1.426 1.929 ( -1.771 -2.261 -0.000) 2.872 1.949 ( 0.254 -1.049 -0.000) 1.079 1.964 ( -1.591 0.345 0.000) 1.628 2.002 ( 0.603 -3.264 -0.000) 3.319 2.031 ( 0.523 0.573 0.000) 0.776 2.055 ( -0.101 0.443 0.000) 0.455 2.076 ( -0.191 -0.250 -0.000) 0.315 2.079 ( -4.779 2.160 0.000) 5.244 2.106 ( -1.269 -0.331 0.000) 1.311 2.150 ( 3.270 0.636 -0.000) 3.332 2.198 ( -2.321 -2.192 0.000) 3.193 2.240 ( -1.791 -1.282 0.000) 2.202 2.268 ( -1.591 3.576 0.000) 3.913 2.301 ( 0.447 3.284 0.000) 3.314 2.320 ( -1.338 1.879 0.000) 2.306 2.332 ( -0.766 1.067 0.000) 1.314 2.385 ( 2.364 0.758 -0.000) 2.482 2.399 ( 1.893 0.862 -0.000) 2.080 2.417 ( 2.233 0.559 -0.000) 2.302 2.441 ( 0.503 -0.688 -0.000) 0.852 2.442 ( -0.626 0.769 0.000) 0.992 2.462 ( 0.374 0.288 -0.000) 0.473 2.472 ( -0.855 -0.530 0.000) 1.006 2.515 ( -0.348 -3.068 -0.000) 3.088 3.148 ( 0.364 0.849 0.000) 0.924 3.166 ( -0.342 -0.098 0.000) 0.356 3.182 ( 0.117 0.155 0.000) 0.194 3.189 ( -0.102 -0.154 -0.000) 0.184 3.305 ( 0.662 1.038 0.000) 1.231 3.323 ( 0.221 0.231 0.000) 0.320 3.336 ( -0.216 -0.465 -0.000) 0.513 3.341 ( -0.636 -0.811 -0.000) 1.030 3.948 ( 0.230 -0.098 -0.000) 0.250 3.951 ( 0.208 -0.249 -0.000) 0.324 3.971 ( 0.076 0.661 0.000) 0.665 3.981 ( -0.509 0.315 0.000) 0.599 3.983 ( -0.122 0.143 0.000) 0.188 3.997 ( 0.466 0.222 -0.000) 0.516 3.998 ( -0.531 -0.857 -0.000) 1.008 4.014 ( -0.214 -0.578 -0.000) 0.617 4.086 ( 0.106 -0.101 -0.000) 0.147 4.089 ( 0.009 -0.271 -0.000) 0.272 4.099 ( 0.437 0.232 -0.000) 0.495 4.105 ( 0.031 0.061 0.000) 0.069 5.566 ( 0.243 0.407 0.000) 0.474 5.569 ( 0.127 -0.002 -0.000) 0.127 5.575 ( -0.099 -0.322 -0.000) 0.337 5.578 ( -0.184 -0.199 -0.000) 0.271 ======================= Grid point 17 (11/18) ======================= q-point: (-0.17 -0.17 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 15 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.242 ( 0.000 0.000 -13.407) 13.407 0.242 ( 0.000 0.000 -13.407) 13.407 0.435 ( 0.000 0.000 -24.285) 24.285 0.825 ( -0.000 -0.000 5.436) 5.436 0.825 ( -0.000 -0.000 5.436) 5.436 0.969 ( 0.000 0.000 -1.631) 1.631 0.969 ( 0.000 0.000 -1.631) 1.631 1.038 ( -0.000 -0.000 5.309) 5.309 1.038 ( -0.000 -0.000 5.309) 5.309 1.089 ( -0.000 -0.000 3.362) 3.362 1.089 ( -0.000 -0.000 3.362) 3.362 1.092 ( -0.000 -0.000 5.230) 5.230 1.141 ( 0.000 0.000 -0.566) 0.566 1.141 ( 0.000 0.000 -0.566) 0.566 1.176 ( 0.000 0.000 -4.274) 4.274 1.228 ( 0.000 0.000 -0.920) 0.920 1.235 ( 0.000 0.000 -2.709) 2.709 1.235 ( 0.000 0.000 -2.709) 2.709 1.252 ( -0.000 -0.000 8.874) 8.874 1.273 ( -0.000 -0.000 0.823) 0.823 1.324 ( -0.000 -0.000 1.502) 1.502 1.340 ( 0.000 0.000 -3.200) 3.200 1.340 ( 0.000 0.000 -3.200) 3.200 1.344 ( -0.000 -0.000 1.515) 1.515 1.370 ( 0.000 0.000 -0.803) 0.803 1.370 ( 0.000 0.000 -0.803) 0.803 1.391 ( 0.000 0.000 -0.837) 0.837 1.433 ( -0.000 -0.000 0.273) 0.273 1.439 ( 0.000 0.000 -1.642) 1.642 1.439 ( 0.000 0.000 -1.642) 1.642 1.458 ( 0.000 0.000 -1.252) 1.252 1.473 ( 0.000 0.000 -4.614) 4.614 1.570 ( 0.000 0.000 -1.749) 1.749 1.570 ( 0.000 0.000 -1.749) 1.749 1.615 ( -0.000 -0.000 0.923) 0.923 1.637 ( 0.000 0.000 -1.772) 1.772 1.646 ( 0.000 0.000 -0.664) 0.664 1.734 ( -0.000 -0.000 2.832) 2.832 1.734 ( -0.000 -0.000 2.832) 2.832 1.770 ( -0.000 -0.000 2.584) 2.584 1.790 ( -0.000 -0.000 2.152) 2.152 1.822 ( -0.000 -0.000 3.600) 3.600 1.823 ( -0.000 -0.000 2.358) 2.358 1.823 ( -0.000 -0.000 2.358) 2.358 1.857 ( 0.000 0.000 -1.845) 1.845 1.903 ( -0.000 -0.000 3.088) 3.088 1.903 ( -0.000 -0.000 1.227) 1.227 1.903 ( -0.000 -0.000 1.227) 1.227 1.969 ( -0.000 -0.000 0.635) 0.635 1.986 ( -0.000 -0.000 0.909) 0.909 1.986 ( -0.000 -0.000 0.909) 0.909 2.004 ( 0.000 0.000 -4.995) 4.995 2.056 ( 0.000 0.000 -0.895) 0.895 2.063 ( -0.000 -0.000 0.156) 0.156 2.063 ( -0.000 -0.000 0.156) 0.156 2.114 ( 0.000 0.000 -3.813) 3.813 2.169 ( -0.000 -0.000 1.956) 1.956 2.178 ( 0.000 0.000 -3.880) 3.880 2.202 ( -0.000 -0.000 0.812) 0.812 2.202 ( -0.000 -0.000 0.812) 0.812 2.232 ( -0.000 -0.000 2.520) 2.520 2.289 ( 0.000 0.000 -2.421) 2.421 2.289 ( 0.000 0.000 -2.421) 2.421 2.379 ( 0.000 0.000 -0.197) 0.197 2.379 ( 0.000 0.000 -0.197) 0.197 2.407 ( -0.000 -0.000 0.934) 0.934 2.409 ( 0.000 0.000 -1.783) 1.783 2.409 ( 0.000 0.000 -1.783) 1.783 2.435 ( 0.000 0.000 -0.074) 0.074 2.448 ( -0.000 -0.000 1.177) 1.177 2.457 ( -0.000 -0.000 0.896) 0.896 2.477 ( 0.000 0.000 -0.764) 0.764 3.128 ( 0.000 0.000 -1.588) 1.588 3.174 ( 0.000 0.000 -0.704) 0.704 3.184 ( -0.000 -0.000 0.485) 0.485 3.186 ( -0.000 -0.000 0.489) 0.489 3.278 ( 0.000 0.000 -0.961) 0.961 3.314 ( 0.000 0.000 -0.042) 0.042 3.350 ( -0.000 -0.000 0.411) 0.411 3.364 ( -0.000 -0.000 0.573) 0.573 3.939 ( -0.000 -0.000 0.188) 0.188 3.939 ( -0.000 -0.000 0.188) 0.188 3.961 ( 0.000 0.000 -0.891) 0.891 3.961 ( 0.000 0.000 -0.891) 0.891 3.987 ( -0.000 -0.000 0.291) 0.291 3.987 ( -0.000 -0.000 0.291) 0.291 4.029 ( -0.000 -0.000 0.594) 0.594 4.029 ( -0.000 -0.000 0.594) 0.594 4.078 ( -0.000 -0.000 0.107) 0.107 4.092 ( 0.000 0.000 -0.297) 0.297 4.094 ( -0.000 -0.000 0.393) 0.393 4.115 ( -0.000 -0.000 0.098) 0.098 5.552 ( 0.000 0.000 -0.194) 0.194 5.569 ( 0.000 0.000 -0.050) 0.050 5.581 ( -0.000 -0.000 0.302) 0.302 5.585 ( -0.000 -0.000 0.178) 0.178 ======================= Grid point 19 (12/18) ======================= q-point: (-0.17 -0.17 -0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 21 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.756 ( -3.001 -3.001 -2.814) 5.092 0.756 ( 3.001 3.001 -2.814) 5.092 0.778 ( -0.000 -0.000 -0.373) 0.373 0.778 ( 0.000 0.000 -0.373) 0.373 0.851 ( -3.425 -3.425 -3.285) 5.853 0.851 ( 3.425 3.425 -3.285) 5.853 0.947 ( -0.092 -0.092 5.884) 5.886 0.947 ( 0.092 0.092 5.884) 5.886 1.037 ( -1.047 -1.047 -0.994) 1.783 1.037 ( 1.047 1.047 -0.994) 1.783 1.062 ( -0.635 -0.635 0.887) 1.263 1.062 ( 0.635 0.635 0.887) 1.263 1.120 ( -0.272 -0.272 -1.428) 1.479 1.120 ( 0.272 0.272 -1.428) 1.479 1.140 ( -0.287 -0.287 1.262) 1.325 1.140 ( 0.287 0.287 1.262) 1.325 1.198 ( -0.171 -0.171 -1.832) 1.847 1.198 ( 0.171 0.171 -1.832) 1.847 1.265 ( -0.649 -0.649 4.161) 4.262 1.265 ( 0.649 0.649 4.161) 4.262 1.285 ( -1.575 -1.575 -3.744) 4.356 1.285 ( 1.575 1.575 -3.744) 4.356 1.322 ( -0.020 -0.020 -0.219) 0.221 1.322 ( 0.020 0.020 -0.219) 0.221 1.406 ( -1.366 -1.366 2.121) 2.869 1.406 ( 1.366 1.366 2.121) 2.869 1.422 ( -0.523 -0.523 2.670) 2.771 1.422 ( 0.523 0.523 2.670) 2.771 1.519 ( -0.000 -0.000 -1.901) 1.901 1.519 ( 0.000 0.000 -1.901) 1.901 1.544 ( -1.170 -1.170 0.212) 1.668 1.544 ( 1.170 1.170 0.212) 1.668 1.579 ( -0.586 -0.586 -1.417) 1.641 1.579 ( 0.586 0.586 -1.417) 1.641 1.602 ( -1.508 -1.508 0.274) 2.150 1.602 ( 1.508 1.508 0.274) 2.150 1.642 ( -1.378 -1.378 -2.658) 3.296 1.642 ( 1.378 1.378 -2.658) 3.296 1.693 ( -0.079 -0.079 2.785) 2.787 1.693 ( 0.079 0.079 2.785) 2.787 1.814 ( -1.963 -1.963 0.665) 2.854 1.814 ( 1.963 1.963 0.665) 2.854 1.814 ( -0.030 -0.030 0.712) 0.714 1.814 ( 0.030 0.030 0.712) 0.714 1.911 ( -0.533 -0.533 -1.238) 1.450 1.911 ( 0.533 0.533 -1.238) 1.450 1.916 ( -1.261 -1.261 -2.414) 3.001 1.916 ( 1.261 1.261 -2.414) 3.001 1.932 ( -1.504 -1.504 0.850) 2.291 1.932 ( 1.504 1.504 0.850) 2.291 1.954 ( -1.786 -1.786 1.412) 2.894 1.954 ( 1.786 1.786 1.412) 2.894 2.043 ( -0.376 -0.376 -0.497) 0.728 2.043 ( 0.376 0.376 -0.497) 0.728 2.056 ( -1.095 -1.095 0.985) 1.835 2.056 ( 1.095 1.095 0.985) 1.835 2.136 ( -1.736 -1.736 -2.621) 3.591 2.136 ( 1.736 1.736 -2.621) 3.591 2.185 ( -0.726 -0.726 2.572) 2.769 2.185 ( 0.726 0.726 2.572) 2.769 2.289 ( -1.972 -1.972 -1.534) 3.183 2.289 ( 1.972 1.972 -1.534) 3.183 2.313 ( -0.513 -0.513 0.878) 1.139 2.313 ( 0.513 0.513 0.878) 1.139 2.428 ( -0.387 -0.387 -1.288) 1.400 2.428 ( 0.387 0.387 -1.288) 1.400 2.446 ( -0.404 -0.404 -0.289) 0.641 2.446 ( 0.404 0.404 -0.289) 0.641 2.454 ( -0.413 -0.413 0.110) 0.595 2.454 ( 0.413 0.413 0.110) 0.595 2.474 ( -1.343 -1.343 1.363) 2.338 2.474 ( 1.343 1.343 1.363) 2.338 3.164 ( -0.401 -0.401 -0.737) 0.930 3.164 ( 0.401 0.401 -0.737) 0.930 3.178 ( -0.135 -0.135 0.727) 0.751 3.178 ( 0.135 0.135 0.727) 0.751 3.325 ( -0.819 -0.819 -0.149) 1.167 3.325 ( 0.819 0.819 -0.149) 1.167 3.328 ( -0.262 -0.262 0.164) 0.405 3.328 ( 0.262 0.262 0.164) 0.405 3.959 ( -0.039 -0.039 -0.626) 0.628 3.959 ( 0.039 0.039 -0.626) 0.628 3.970 ( -0.043 -0.043 0.434) 0.439 3.970 ( 0.043 0.043 0.434) 0.439 3.990 ( -0.192 -0.192 -0.460) 0.534 3.990 ( 0.192 0.192 -0.460) 0.534 4.001 ( -0.382 -0.382 0.632) 0.831 4.001 ( 0.382 0.382 0.632) 0.831 4.093 ( -0.082 -0.082 -0.417) 0.432 4.093 ( 0.082 0.082 -0.417) 0.432 4.101 ( -0.101 -0.101 0.411) 0.435 4.101 ( 0.101 0.101 0.411) 0.435 5.572 ( -0.134 -0.134 -0.012) 0.190 5.572 ( 0.134 0.134 -0.012) 0.190 5.573 ( -0.260 -0.260 0.010) 0.368 5.573 ( 0.260 0.260 0.010) 0.368 ======================= Grid point 22 (13/18) ======================= q-point: (-0.46 -0.21 0.21) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.561 ( 0.000 -8.437 -8.778) 12.175 0.606 ( 0.000 -10.902 -6.830) 12.865 0.759 ( -0.000 -0.407 4.462) 4.480 0.777 ( 0.000 -3.567 2.380) 4.288 0.881 ( 0.000 -4.031 -12.018) 12.676 0.882 ( -0.000 2.340 7.192) 7.563 0.915 ( -0.000 8.988 3.076) 9.499 0.961 ( 0.000 -2.964 -0.590) 3.022 0.970 ( -0.000 1.516 0.323) 1.550 1.002 ( 0.000 -0.543 0.019) 0.543 1.018 ( -0.000 -0.593 2.479) 2.549 1.039 ( -0.000 0.979 2.221) 2.427 1.088 ( -0.000 -0.766 6.747) 6.791 1.132 ( 0.000 0.948 -1.231) 1.553 1.145 ( -0.000 -0.907 0.907) 1.283 1.220 ( 0.000 1.291 -1.280) 1.818 1.242 ( -0.000 1.205 0.341) 1.253 1.242 ( 0.000 1.041 -0.975) 1.426 1.259 ( 0.000 -0.927 0.583) 1.096 1.290 ( 0.000 1.139 -1.356) 1.771 1.298 ( -0.000 -1.374 1.811) 2.274 1.301 ( 0.000 2.302 -2.057) 3.087 1.362 ( -0.000 0.392 0.704) 0.806 1.364 ( 0.000 2.423 -3.584) 4.327 1.399 ( -0.000 2.348 1.185) 2.630 1.401 ( 0.000 -1.871 -1.082) 2.161 1.414 ( 0.000 -0.362 -0.489) 0.608 1.457 ( 0.000 -2.302 -0.556) 2.368 1.485 ( 0.000 1.120 -4.078) 4.229 1.501 ( 0.000 -0.085 -1.983) 1.985 1.522 ( 0.000 -1.424 -1.623) 2.160 1.527 ( -0.000 2.514 -0.700) 2.609 1.548 ( 0.000 2.003 -2.272) 3.029 1.589 ( 0.000 -1.816 -2.818) 3.352 1.613 ( -0.000 2.394 2.425) 3.407 1.622 ( -0.000 -0.728 2.383) 2.492 1.658 ( 0.000 -0.637 0.383) 0.743 1.698 ( -0.000 1.220 2.632) 2.901 1.705 ( 0.000 -1.582 0.518) 1.664 1.712 ( -0.000 0.625 1.514) 1.638 1.734 ( -0.000 0.917 1.328) 1.614 1.777 ( -0.000 -0.644 2.696) 2.772 1.796 ( 0.000 -0.979 -0.526) 1.111 1.803 ( -0.000 -2.429 2.461) 3.458 1.809 ( -0.000 1.772 0.573) 1.863 1.879 ( -0.000 -1.409 1.699) 2.207 1.886 ( -0.000 0.490 0.268) 0.559 1.913 ( 0.000 -1.454 0.877) 1.698 1.931 ( -0.000 0.776 0.754) 1.082 1.935 ( 0.000 -2.322 1.364) 2.693 2.010 ( -0.000 -1.314 2.112) 2.487 2.043 ( -0.000 2.475 2.837) 3.765 2.052 ( -0.000 1.795 -1.170) 2.142 2.068 ( 0.000 0.709 -3.071) 3.152 2.076 ( 0.000 2.140 -2.880) 3.588 2.108 ( -0.000 1.244 1.714) 2.118 2.159 ( -0.000 1.520 1.352) 2.034 2.167 ( -0.000 -1.059 1.092) 1.521 2.201 ( 0.000 -1.030 -1.456) 1.784 2.225 ( 0.000 -2.767 0.672) 2.847 2.259 ( 0.000 0.431 -5.549) 5.566 2.324 ( 0.000 1.198 -1.782) 2.147 2.346 ( -0.000 1.875 -1.522) 2.415 2.359 ( 0.000 0.070 -0.644) 0.648 2.379 ( 0.000 0.070 -1.473) 1.475 2.388 ( 0.000 -1.464 -1.005) 1.776 2.399 ( -0.000 0.832 -0.558) 1.002 2.419 ( -0.000 1.853 0.667) 1.970 2.445 ( 0.000 -1.144 0.274) 1.177 2.446 ( -0.000 -0.825 0.805) 1.153 2.483 ( -0.000 -0.480 0.885) 1.007 2.495 ( -0.000 -1.348 1.746) 2.206 3.164 ( 0.000 -0.314 -1.576) 1.607 3.169 ( -0.000 0.595 0.570) 0.824 3.175 ( -0.000 -0.014 0.466) 0.467 3.186 ( 0.000 -0.010 -0.418) 0.418 3.303 ( 0.000 -0.340 -1.309) 1.352 3.316 ( -0.000 0.315 -0.146) 0.348 3.339 ( -0.000 -0.071 0.465) 0.471 3.346 ( -0.000 0.189 0.936) 0.955 3.943 ( 0.000 -0.429 0.190) 0.469 3.944 ( -0.000 -0.449 0.387) 0.593 3.976 ( 0.000 -0.202 0.137) 0.244 3.976 ( 0.000 -0.420 -0.576) 0.713 3.988 ( 0.000 -0.015 -0.319) 0.320 3.989 ( 0.000 -0.031 -0.556) 0.557 4.010 ( -0.000 0.475 0.768) 0.903 4.011 ( -0.000 0.302 0.520) 0.601 4.080 ( 0.000 -0.258 0.047) 0.263 4.088 ( -0.000 -0.230 0.485) 0.537 4.097 ( 0.000 0.171 -0.368) 0.406 4.106 ( -0.000 0.343 -0.050) 0.346 5.560 ( 0.000 -0.258 -0.323) 0.413 5.570 ( 0.000 0.126 -0.118) 0.173 5.575 ( -0.000 -0.072 0.401) 0.407 5.579 ( -0.000 0.070 0.204) 0.216 ======================= Grid point 26 (14/18) ======================= q-point: ( 0.25 -0.25 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 26 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.552 ( -0.000 10.957 0.000) 10.957 0.634 ( -0.000 12.961 0.000) 12.961 0.847 ( 0.000 -5.445 -0.000) 5.445 0.853 ( 0.000 -6.568 -0.000) 6.568 0.856 ( 0.000 -4.305 -0.000) 4.305 0.915 ( -0.000 7.217 0.000) 7.217 0.928 ( 0.000 -6.153 -0.000) 6.153 0.935 ( 0.000 -5.312 -0.000) 5.312 0.976 ( -0.000 0.194 0.000) 0.194 0.990 ( -0.000 0.942 0.000) 0.942 1.027 ( 0.000 -2.422 -0.000) 2.422 1.087 ( 0.000 -2.907 -0.000) 2.907 1.114 ( -0.000 6.573 0.000) 6.573 1.132 ( -0.000 1.161 0.000) 1.161 1.141 ( -0.000 0.604 0.000) 0.604 1.178 ( 0.000 -1.691 -0.000) 1.691 1.181 ( 0.000 -1.160 -0.000) 1.160 1.216 ( -0.000 2.277 0.000) 2.277 1.231 ( -0.000 0.156 0.000) 0.156 1.256 ( -0.000 0.810 0.000) 0.810 1.272 ( -0.000 1.760 0.000) 1.760 1.329 ( 0.000 -0.070 -0.000) 0.070 1.342 ( 0.000 -1.401 -0.000) 1.401 1.368 ( -0.000 1.478 0.000) 1.478 1.372 ( 0.000 -0.874 -0.000) 0.874 1.418 ( -0.000 1.056 0.000) 1.056 1.437 ( -0.000 0.445 0.000) 0.445 1.465 ( -0.000 0.583 0.000) 0.583 1.488 ( -0.000 1.717 0.000) 1.717 1.488 ( -0.000 2.662 0.000) 2.662 1.528 ( -0.000 1.449 0.000) 1.449 1.536 ( -0.000 0.126 0.000) 0.126 1.549 ( 0.000 -1.708 -0.000) 1.708 1.601 ( -0.000 1.854 0.000) 1.854 1.604 ( 0.000 -0.420 -0.000) 0.420 1.613 ( 0.000 -2.756 -0.000) 2.756 1.621 ( 0.000 -0.298 -0.000) 0.298 1.685 ( 0.000 -2.049 -0.000) 2.049 1.716 ( 0.000 -2.293 -0.000) 2.293 1.749 ( -0.000 0.008 0.000) 0.008 1.766 ( -0.000 1.636 0.000) 1.636 1.776 ( 0.000 -3.793 -0.000) 3.793 1.787 ( 0.000 -1.160 -0.000) 1.160 1.807 ( 0.000 -2.877 -0.000) 2.877 1.848 ( -0.000 0.990 0.000) 0.990 1.861 ( 0.000 -1.502 -0.000) 1.502 1.913 ( 0.000 -3.148 -0.000) 3.148 1.918 ( -0.000 3.157 0.000) 3.157 1.941 ( 0.000 -1.571 -0.000) 1.571 1.959 ( 0.000 -3.289 -0.000) 3.289 1.997 ( -0.000 0.949 0.000) 0.949 2.018 ( 0.000 -0.162 -0.000) 0.162 2.024 ( -0.000 4.196 0.000) 4.196 2.063 ( 0.000 -0.187 -0.000) 0.187 2.074 ( 0.000 -5.036 -0.000) 5.036 2.088 ( -0.000 1.902 0.000) 1.902 2.120 ( 0.000 -1.145 -0.000) 1.145 2.172 ( 0.000 -0.684 -0.000) 0.684 2.192 ( -0.000 3.108 0.000) 3.108 2.266 ( 0.000 -0.620 -0.000) 0.620 2.297 ( -0.000 4.057 0.000) 4.057 2.318 ( -0.000 1.199 0.000) 1.199 2.325 ( -0.000 1.370 0.000) 1.370 2.334 ( -0.000 0.443 0.000) 0.443 2.356 ( -0.000 3.250 0.000) 3.250 2.362 ( -0.000 0.003 0.000) 0.003 2.382 ( -0.000 1.062 0.000) 1.062 2.440 ( -0.000 0.101 0.000) 0.101 2.457 ( 0.000 -0.852 -0.000) 0.852 2.460 ( 0.000 -0.754 -0.000) 0.754 2.484 ( -0.000 0.245 0.000) 0.245 2.512 ( 0.000 -3.197 -0.000) 3.197 3.144 ( -0.000 1.044 0.000) 1.044 3.170 ( -0.000 0.057 0.000) 0.057 3.180 ( -0.000 0.109 0.000) 0.109 3.190 ( 0.000 -0.107 -0.000) 0.107 3.298 ( -0.000 1.079 0.000) 1.079 3.321 ( -0.000 0.204 0.000) 0.204 3.338 ( 0.000 -0.531 -0.000) 0.531 3.348 ( 0.000 -0.758 -0.000) 0.758 3.946 ( 0.000 -0.149 -0.000) 0.149 3.948 ( 0.000 -0.275 -0.000) 0.275 3.970 ( -0.000 0.739 0.000) 0.739 3.987 ( 0.000 -0.047 -0.000) 0.047 3.988 ( -0.000 0.411 0.000) 0.411 3.989 ( -0.000 0.226 0.000) 0.226 4.006 ( 0.000 -1.183 -0.000) 1.183 4.016 ( 0.000 -0.549 -0.000) 0.549 4.085 ( 0.000 -0.178 -0.000) 0.178 4.089 ( 0.000 -0.344 -0.000) 0.344 4.092 ( -0.000 0.456 0.000) 0.456 4.104 ( -0.000 0.151 0.000) 0.151 5.564 ( -0.000 0.451 0.000) 0.451 5.568 ( 0.000 -0.068 -0.000) 0.068 5.576 ( 0.000 -0.390 -0.000) 0.390 5.580 ( 0.000 -0.157 -0.000) 0.157 ======================= Grid point 27 (15/18) ======================= q-point: ( 0.25 0.75 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 20 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.692 ( 3.052 -3.074 0.000) 4.332 0.692 ( 3.074 -3.052 0.000) 4.332 0.795 ( -4.167 -1.034 0.000) 4.293 0.795 ( 1.034 4.167 0.000) 4.293 0.840 ( 5.296 -5.695 0.000) 7.777 0.840 ( 5.695 -5.296 0.000) 7.777 0.910 ( -6.175 0.034 0.000) 6.176 0.910 ( -0.034 6.175 0.000) 6.176 0.997 ( -1.008 -1.986 -0.000) 2.227 0.997 ( 1.986 1.008 -0.000) 2.227 1.008 ( 0.066 -0.609 0.000) 0.613 1.008 ( 0.609 -0.066 0.000) 0.613 1.127 ( -0.470 0.444 0.000) 0.647 1.127 ( -0.444 0.470 0.000) 0.647 1.182 ( 1.019 -1.116 0.000) 1.511 1.182 ( 1.116 -1.019 0.000) 1.511 1.244 ( -1.699 -0.312 0.000) 1.727 1.244 ( 0.312 1.699 0.000) 1.727 1.256 ( -1.034 -1.459 -0.000) 1.788 1.256 ( 1.459 1.034 -0.000) 1.788 1.292 ( -2.238 0.507 0.000) 2.294 1.292 ( -0.507 2.238 0.000) 2.294 1.366 ( -1.454 0.384 0.000) 1.504 1.366 ( -0.384 1.454 0.000) 1.504 1.402 ( -0.438 0.329 0.000) 0.548 1.402 ( -0.329 0.438 0.000) 0.548 1.475 ( 0.956 -1.096 0.000) 1.455 1.475 ( 1.096 -0.956 0.000) 1.455 1.507 ( -1.047 0.723 0.000) 1.272 1.507 ( -0.723 1.047 0.000) 1.272 1.531 ( 0.878 -1.313 -0.000) 1.580 1.531 ( 1.313 -0.878 -0.000) 1.580 1.590 ( -3.216 2.010 0.000) 3.793 1.590 ( -2.010 3.216 0.000) 3.793 1.603 ( 0.242 -0.265 0.000) 0.359 1.603 ( 0.265 -0.242 0.000) 0.359 1.620 ( -2.254 -1.794 0.000) 2.880 1.620 ( 1.794 2.254 0.000) 2.880 1.722 ( -0.378 -1.253 0.000) 1.309 1.722 ( 1.253 0.378 0.000) 1.309 1.760 ( -2.334 -4.078 -0.000) 4.699 1.760 ( 4.078 2.334 -0.000) 4.699 1.791 ( 2.126 -2.307 0.000) 3.137 1.791 ( 2.307 -2.126 0.000) 3.137 1.867 ( 0.796 -0.804 0.000) 1.131 1.867 ( 0.804 -0.796 0.000) 1.131 1.884 ( -1.000 -0.020 0.000) 1.000 1.884 ( 0.020 1.000 0.000) 1.000 1.930 ( 1.423 -1.915 -0.000) 2.386 1.930 ( 1.915 -1.423 -0.000) 2.386 1.956 ( -2.037 0.504 0.000) 2.099 1.956 ( -0.504 2.037 0.000) 2.099 2.069 ( -1.191 0.813 0.000) 1.442 2.069 ( -0.813 1.191 0.000) 1.442 2.110 ( -2.401 0.727 0.000) 2.508 2.110 ( -0.727 2.401 0.000) 2.508 2.156 ( -4.093 -1.648 0.000) 4.412 2.156 ( 1.648 4.093 0.000) 4.412 2.190 ( -2.064 -3.409 -0.000) 3.985 2.190 ( 3.409 2.064 -0.000) 3.985 2.341 ( -1.352 0.112 0.000) 1.357 2.341 ( -0.112 1.352 0.000) 1.357 2.347 ( -1.143 1.030 0.000) 1.538 2.347 ( -1.030 1.143 0.000) 1.538 2.396 ( -0.400 -1.742 -0.000) 1.788 2.396 ( 1.742 0.400 -0.000) 1.788 2.424 ( -0.405 -2.200 -0.000) 2.237 2.424 ( 2.200 0.405 -0.000) 2.237 2.447 ( -0.495 -0.508 0.000) 0.709 2.447 ( 0.508 0.495 0.000) 0.709 2.467 ( -0.335 -0.961 0.000) 1.017 2.467 ( 0.961 0.335 0.000) 1.017 3.164 ( -0.589 -0.086 0.000) 0.595 3.164 ( 0.086 0.589 0.000) 0.595 3.185 ( 0.001 -0.076 0.000) 0.076 3.185 ( 0.076 -0.001 0.000) 0.076 3.325 ( -0.602 -0.451 0.000) 0.752 3.325 ( 0.451 0.602 0.000) 0.752 3.327 ( 0.029 -0.106 0.000) 0.110 3.327 ( 0.106 -0.029 0.000) 0.110 3.944 ( 0.273 -0.304 -0.000) 0.409 3.944 ( 0.304 -0.273 -0.000) 0.409 3.981 ( -0.119 -0.523 -0.000) 0.536 3.981 ( 0.523 0.119 -0.000) 0.536 3.991 ( -0.509 0.392 0.000) 0.642 3.991 ( -0.392 0.509 0.000) 0.642 4.002 ( -0.260 -0.495 -0.000) 0.560 4.002 ( 0.495 0.260 -0.000) 0.560 4.084 ( 0.167 -0.191 0.000) 0.254 4.084 ( 0.191 -0.167 0.000) 0.254 4.104 ( -0.133 -0.093 0.000) 0.162 4.104 ( 0.093 0.133 0.000) 0.162 5.569 ( 0.095 -0.163 0.000) 0.189 5.569 ( 0.163 -0.095 0.000) 0.189 5.574 ( -0.201 -0.105 0.000) 0.227 5.574 ( 0.105 0.201 0.000) 0.227 ======================= Grid point 35 (16/18) ======================= q-point: (-0.33 -0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 18 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.458 ( 0.000 0.000 -10.935) 10.935 0.458 ( 0.000 0.000 -10.935) 10.935 0.733 ( -0.000 -0.000 4.872) 4.872 0.733 ( -0.000 -0.000 4.872) 4.872 0.834 ( 0.000 0.000 -20.950) 20.950 0.930 ( -0.000 -0.000 5.192) 5.192 0.930 ( -0.000 -0.000 5.192) 5.192 0.993 ( 0.000 0.000 -0.113) 0.113 0.993 ( 0.000 0.000 -0.113) 0.113 1.001 ( -0.000 -0.000 4.165) 4.165 1.038 ( -0.000 -0.000 2.355) 2.355 1.038 ( -0.000 -0.000 2.355) 2.355 1.085 ( -0.000 -0.000 8.901) 8.901 1.143 ( -0.000 -0.000 0.172) 0.172 1.143 ( -0.000 -0.000 0.172) 0.172 1.221 ( 0.000 0.000 -1.042) 1.042 1.255 ( -0.000 -0.000 0.121) 0.121 1.256 ( -0.000 -0.000 1.088) 1.088 1.268 ( 0.000 0.000 -0.775) 0.775 1.268 ( 0.000 0.000 -0.775) 0.775 1.303 ( -0.000 -0.000 0.651) 0.651 1.318 ( -0.000 -0.000 1.294) 1.294 1.374 ( -0.000 -0.000 0.945) 0.945 1.374 ( -0.000 -0.000 0.945) 0.945 1.397 ( 0.000 0.000 -4.250) 4.250 1.397 ( 0.000 0.000 -4.250) 4.250 1.407 ( 0.000 0.000 -0.693) 0.693 1.424 ( -0.000 -0.000 0.584) 0.584 1.488 ( 0.000 0.000 -2.140) 2.140 1.501 ( 0.000 0.000 -4.653) 4.653 1.501 ( 0.000 0.000 -4.653) 4.653 1.558 ( 0.000 0.000 -4.884) 4.884 1.590 ( -0.000 -0.000 2.426) 2.426 1.599 ( 0.000 0.000 -0.439) 0.439 1.599 ( 0.000 0.000 -0.439) 0.439 1.653 ( 0.000 0.000 -0.422) 0.422 1.658 ( 0.000 0.000 -0.272) 0.272 1.688 ( -0.000 -0.000 1.547) 1.547 1.688 ( -0.000 -0.000 1.547) 1.547 1.724 ( -0.000 -0.000 2.727) 2.727 1.758 ( -0.000 -0.000 1.337) 1.337 1.763 ( -0.000 -0.000 4.024) 4.024 1.763 ( -0.000 -0.000 4.024) 4.024 1.771 ( -0.000 -0.000 2.348) 2.348 1.856 ( -0.000 -0.000 0.429) 0.429 1.865 ( -0.000 -0.000 2.544) 2.544 1.865 ( -0.000 -0.000 2.544) 2.544 1.884 ( 0.000 0.000 -1.108) 1.108 1.961 ( -0.000 -0.000 0.001) 0.001 1.973 ( -0.000 -0.000 0.928) 0.928 1.973 ( -0.000 -0.000 0.928) 0.928 2.079 ( 0.000 0.000 -2.239) 2.239 2.079 ( 0.000 0.000 -2.239) 2.239 2.082 ( 0.000 0.000 -4.186) 4.186 2.086 ( 0.000 0.000 -2.349) 2.349 2.101 ( -0.000 -0.000 6.000) 6.000 2.166 ( 0.000 0.000 -1.223) 1.223 2.179 ( -0.000 -0.000 1.895) 1.895 2.179 ( -0.000 -0.000 1.895) 1.895 2.180 ( -0.000 -0.000 2.918) 2.918 2.262 ( 0.000 0.000 -5.163) 5.163 2.335 ( 0.000 0.000 -2.493) 2.493 2.335 ( 0.000 0.000 -2.493) 2.493 2.382 ( -0.000 -0.000 0.080) 0.080 2.382 ( -0.000 -0.000 0.080) 0.080 2.388 ( -0.000 -0.000 1.381) 1.381 2.432 ( -0.000 -0.000 0.481) 0.481 2.435 ( -0.000 -0.000 0.177) 0.177 2.438 ( 0.000 0.000 -1.358) 1.358 2.438 ( 0.000 0.000 -1.358) 1.358 2.441 ( -0.000 -0.000 0.296) 0.296 2.482 ( -0.000 -0.000 0.450) 0.450 3.160 ( 0.000 0.000 -1.796) 1.796 3.175 ( -0.000 -0.000 0.437) 0.437 3.176 ( -0.000 -0.000 0.510) 0.510 3.186 ( 0.000 0.000 -0.419) 0.419 3.300 ( 0.000 0.000 -1.458) 1.458 3.318 ( 0.000 0.000 -0.371) 0.371 3.340 ( -0.000 -0.000 0.699) 0.699 3.349 ( -0.000 -0.000 1.117) 1.117 3.937 ( -0.000 -0.000 0.093) 0.093 3.937 ( -0.000 -0.000 0.093) 0.093 3.971 ( 0.000 0.000 -0.148) 0.148 3.971 ( 0.000 0.000 -0.148) 0.148 3.989 ( 0.000 0.000 -0.564) 0.564 3.989 ( 0.000 0.000 -0.564) 0.564 4.016 ( -0.000 -0.000 0.820) 0.820 4.016 ( -0.000 -0.000 0.820) 0.820 4.077 ( -0.000 -0.000 0.019) 0.019 4.085 ( -0.000 -0.000 0.486) 0.486 4.098 ( 0.000 0.000 -0.445) 0.445 4.112 ( -0.000 -0.000 0.249) 0.249 5.557 ( 0.000 0.000 -0.380) 0.380 5.571 ( 0.000 0.000 -0.174) 0.174 5.574 ( -0.000 -0.000 0.487) 0.487 5.581 ( -0.000 -0.000 0.298) 0.298 ======================= Grid point 40 (17/18) ======================= q-point: ( 0.37 0.62 -0.62) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 20 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.661 ( 0.000 -5.100 -0.000) 5.100 0.675 ( 0.000 -7.379 -0.000) 7.379 0.728 ( 0.000 -5.199 -0.000) 5.199 0.758 ( 0.000 -7.455 -0.000) 7.455 0.784 ( -0.000 5.839 0.000) 5.839 0.870 ( -0.000 3.716 0.000) 3.716 0.903 ( -0.000 4.956 0.000) 4.956 0.968 ( -0.000 1.003 0.000) 1.003 1.000 ( 0.000 -0.290 -0.000) 0.290 1.002 ( 0.000 -0.573 -0.000) 0.573 1.012 ( 0.000 -1.599 -0.000) 1.599 1.019 ( -0.000 0.233 0.000) 0.233 1.137 ( 0.000 -0.360 -0.000) 0.360 1.138 ( 0.000 -1.262 -0.000) 1.262 1.147 ( 0.000 -0.160 -0.000) 0.160 1.170 ( 0.000 -0.983 -0.000) 0.983 1.237 ( -0.000 0.274 0.000) 0.274 1.238 ( -0.000 2.143 0.000) 2.143 1.258 ( -0.000 1.222 0.000) 1.222 1.275 ( 0.000 -3.110 -0.000) 3.110 1.311 ( -0.000 1.088 0.000) 1.088 1.317 ( 0.000 -0.003 -0.000) 0.003 1.344 ( 0.000 -0.683 -0.000) 0.683 1.395 ( -0.000 6.397 0.000) 6.397 1.408 ( 0.000 -2.929 -0.000) 2.929 1.416 ( -0.000 2.573 0.000) 2.573 1.422 ( -0.000 0.018 0.000) 0.018 1.463 ( 0.000 -2.245 -0.000) 2.245 1.506 ( -0.000 1.616 0.000) 1.616 1.521 ( -0.000 0.893 0.000) 0.893 1.525 ( -0.000 1.015 0.000) 1.015 1.546 ( 0.000 -0.048 -0.000) 0.048 1.577 ( -0.000 2.117 0.000) 2.117 1.588 ( -0.000 0.595 0.000) 0.595 1.615 ( -0.000 0.428 0.000) 0.428 1.630 ( -0.000 0.558 0.000) 0.558 1.650 ( 0.000 -0.499 -0.000) 0.499 1.652 ( -0.000 4.015 0.000) 4.015 1.678 ( 0.000 -3.121 -0.000) 3.121 1.694 ( 0.000 -0.506 -0.000) 0.506 1.734 ( 0.000 -1.595 -0.000) 1.595 1.757 ( 0.000 -1.012 -0.000) 1.012 1.772 ( 0.000 -1.752 -0.000) 1.752 1.792 ( 0.000 -2.134 -0.000) 2.134 1.831 ( 0.000 -0.574 -0.000) 0.574 1.861 ( 0.000 -1.801 -0.000) 1.801 1.880 ( -0.000 3.032 0.000) 3.032 1.910 ( 0.000 -1.408 -0.000) 1.408 1.922 ( 0.000 -2.623 -0.000) 2.623 1.922 ( -0.000 0.629 0.000) 0.629 1.995 ( 0.000 -0.105 0.000) 0.105 1.996 ( 0.000 -0.031 0.000) 0.031 2.073 ( -0.000 3.617 0.000) 3.617 2.084 ( -0.000 1.607 0.000) 1.607 2.104 ( -0.000 1.256 0.000) 1.256 2.118 ( -0.000 2.290 0.000) 2.290 2.128 ( -0.000 1.562 0.000) 1.562 2.157 ( 0.000 -0.954 -0.000) 0.954 2.209 ( 0.000 -2.200 -0.000) 2.200 2.218 ( 0.000 -3.820 -0.000) 3.820 2.335 ( -0.000 0.141 0.000) 0.141 2.339 ( -0.000 1.215 0.000) 1.215 2.356 ( -0.000 0.530 0.000) 0.530 2.363 ( -0.000 1.875 0.000) 1.875 2.375 ( -0.000 0.486 0.000) 0.486 2.384 ( -0.000 0.022 0.000) 0.022 2.424 ( 0.000 -1.811 -0.000) 1.811 2.430 ( 0.000 -1.007 -0.000) 1.007 2.433 ( -0.000 2.012 0.000) 2.012 2.445 ( -0.000 0.626 0.000) 0.626 2.458 ( 0.000 -0.801 -0.000) 0.801 2.481 ( 0.000 -0.876 -0.000) 0.876 3.163 ( -0.000 0.699 0.000) 0.699 3.171 ( -0.000 0.011 0.000) 0.011 3.183 ( -0.000 0.069 0.000) 0.069 3.186 ( 0.000 -0.206 -0.000) 0.206 3.319 ( -0.000 0.769 0.000) 0.769 3.324 ( -0.000 0.126 0.000) 0.126 3.329 ( 0.000 -0.290 -0.000) 0.290 3.334 ( 0.000 -0.535 -0.000) 0.535 3.940 ( 0.000 -0.351 -0.000) 0.351 3.940 ( 0.000 -0.309 -0.000) 0.309 3.975 ( 0.000 -0.202 -0.000) 0.202 3.981 ( 0.000 -0.859 -0.000) 0.859 3.995 ( -0.000 0.625 0.000) 0.625 3.996 ( -0.000 0.327 0.000) 0.327 3.998 ( -0.000 0.370 0.000) 0.370 4.006 ( 0.000 -0.435 -0.000) 0.435 4.081 ( 0.000 -0.182 -0.000) 0.182 4.082 ( 0.000 -0.282 -0.000) 0.282 4.101 ( -0.000 0.393 0.000) 0.393 4.106 ( -0.000 0.014 0.000) 0.014 5.566 ( 0.000 -0.066 -0.000) 0.066 5.568 ( 0.000 -0.290 -0.000) 0.290 5.572 ( -0.000 0.296 0.000) 0.296 5.577 ( 0.000 -0.131 -0.000) 0.131 ======================= Grid point 54 (18/18) ======================= q-point: ( 0.50 0.50 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 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.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 1.10e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 3.34e-04 Number of triplets: 12 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.586 ( 0.000 0.000 -0.000) 0.000 0.586 ( 0.000 0.000 -0.000) 0.000 0.668 ( 0.000 0.000 -0.000) 0.000 0.668 ( 0.000 0.000 -0.000) 0.000 0.873 ( 0.000 0.000 -0.000) 0.000 0.873 ( 0.000 0.000 -0.000) 0.000 0.980 ( 0.000 0.000 -2.428) 2.428 0.980 ( -0.000 -0.000 2.428) 2.428 0.996 ( 0.000 0.000 -0.000) 0.000 0.996 ( 0.000 0.000 -0.000) 0.000 1.014 ( 0.000 0.000 -0.000) 0.000 1.014 ( 0.000 0.000 -0.000) 0.000 1.141 ( 0.000 0.000 -13.116) 13.116 1.141 ( 0.000 0.000 0.010) 0.010 1.141 ( -0.000 -0.000 -0.010) 0.010 1.141 ( -0.000 -0.000 13.116) 13.116 1.238 ( 0.000 0.000 -0.951) 0.951 1.238 ( -0.000 -0.000 0.951) 0.951 1.271 ( 0.000 0.000 -0.000) 0.000 1.271 ( 0.000 0.000 -0.000) 0.000 1.302 ( 0.000 0.000 -0.477) 0.477 1.302 ( -0.000 -0.000 0.477) 0.477 1.365 ( 0.000 0.000 -0.000) 0.000 1.365 ( 0.000 0.000 -0.000) 0.000 1.415 ( 0.000 0.000 -0.373) 0.373 1.415 ( -0.000 -0.000 0.373) 0.373 1.450 ( 0.000 0.000 -0.000) 0.000 1.450 ( 0.000 0.000 -0.000) 0.000 1.534 ( 0.000 0.000 -3.158) 3.158 1.534 ( -0.000 -0.000 3.158) 3.158 1.546 ( 0.000 0.000 -0.000) 0.000 1.546 ( 0.000 0.000 -0.000) 0.000 1.591 ( 0.000 0.000 -0.000) 0.000 1.591 ( 0.000 0.000 -0.000) 0.000 1.633 ( 0.000 0.000 -3.254) 3.254 1.633 ( -0.000 -0.000 3.254) 3.254 1.678 ( 0.000 0.000 -2.315) 2.315 1.678 ( -0.000 -0.000 2.315) 2.315 1.690 ( 0.000 0.000 -0.000) 0.000 1.690 ( 0.000 0.000 -0.000) 0.000 1.709 ( -0.000 -0.000 0.000) 0.000 1.709 ( -0.000 -0.000 0.000) 0.000 1.747 ( 0.000 0.000 -0.237) 0.237 1.747 ( -0.000 -0.000 0.237) 0.237 1.839 ( 0.000 0.000 -0.000) 0.000 1.839 ( 0.000 0.000 -0.000) 0.000 1.885 ( 0.000 0.000 -1.499) 1.499 1.885 ( -0.000 -0.000 1.499) 1.499 1.961 ( 0.000 0.000 -0.000) 0.000 1.961 ( 0.000 0.000 -0.000) 0.000 1.993 ( 0.000 0.000 -4.704) 4.704 1.993 ( -0.000 -0.000 4.704) 4.704 2.117 ( 0.000 0.000 -0.000) 0.000 2.117 ( 0.000 0.000 -0.000) 0.000 2.131 ( 0.000 0.000 -2.749) 2.749 2.131 ( -0.000 -0.000 2.749) 2.749 2.147 ( 0.000 0.000 -0.000) 0.000 2.147 ( 0.000 0.000 -0.000) 0.000 2.147 ( 0.000 0.000 -2.907) 2.907 2.147 ( -0.000 -0.000 2.907) 2.907 2.343 ( 0.000 0.000 -3.838) 3.838 2.343 ( -0.000 -0.000 3.838) 3.838 2.365 ( -0.000 -0.000 0.000) 0.000 2.365 ( -0.000 -0.000 0.000) 0.000 2.374 ( 0.000 0.000 -0.000) 0.000 2.374 ( 0.000 0.000 -0.000) 0.000 2.431 ( 0.000 0.000 -0.196) 0.196 2.431 ( -0.000 -0.000 0.196) 0.196 2.451 ( 0.000 0.000 -0.000) 0.000 2.451 ( 0.000 0.000 -0.000) 0.000 2.459 ( 0.000 0.000 -1.820) 1.820 2.459 ( -0.000 -0.000 1.820) 1.820 3.171 ( 0.000 0.000 -0.064) 0.064 3.171 ( -0.000 -0.000 0.064) 0.064 3.183 ( 0.000 0.000 -0.768) 0.768 3.183 ( -0.000 -0.000 0.768) 0.768 3.326 ( 0.000 0.000 -1.469) 1.469 3.326 ( -0.000 -0.000 1.469) 1.469 3.327 ( 0.000 0.000 -0.714) 0.714 3.327 ( -0.000 -0.000 0.714) 0.714 3.936 ( 0.000 0.000 -0.000) 0.000 3.936 ( 0.000 0.000 -0.000) 0.000 3.971 ( 0.000 0.000 -0.000) 0.000 3.971 ( 0.000 0.000 -0.000) 0.000 4.000 ( 0.000 0.000 -0.000) 0.000 4.000 ( 0.000 0.000 -0.000) 0.000 4.003 ( 0.000 0.000 -0.000) 0.000 4.003 ( 0.000 0.000 -0.000) 0.000 4.079 ( 0.000 0.000 -0.261) 0.261 4.079 ( -0.000 -0.000 0.261) 0.261 4.106 ( 0.000 0.000 -0.403) 0.403 4.106 ( -0.000 -0.000 0.403) 0.403 5.565 ( 0.000 0.000 -0.503) 0.503 5.565 ( -0.000 -0.000 0.503) 0.503 5.575 ( 0.000 0.000 -0.294) 0.294 5.575 ( -0.000 -0.000 0.294) 0.294 =================== 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/9216 10.0 28.066 28.066 33.397 -0.000 -0.000 -0.000 3/9216 20.0 5.607 5.607 6.008 -0.000 -0.000 -0.000 3/9216 30.0 2.760 2.760 2.960 -0.000 -0.000 -0.000 3/9216 40.0 1.829 1.829 1.971 -0.000 -0.000 -0.000 3/9216 50.0 1.377 1.377 1.489 -0.000 -0.000 -0.000 3/9216 60.0 1.109 1.109 1.202 -0.000 -0.000 -0.000 3/9216 70.0 0.931 0.931 1.011 -0.000 -0.000 -0.000 3/9216 80.0 0.804 0.804 0.873 -0.000 -0.000 -0.000 3/9216 90.0 0.708 0.708 0.770 -0.000 -0.000 -0.000 3/9216 100.0 0.633 0.633 0.689 -0.000 -0.000 -0.000 3/9216 110.0 0.573 0.573 0.623 -0.000 -0.000 -0.000 3/9216 120.0 0.523 0.523 0.570 -0.000 -0.000 -0.000 3/9216 130.0 0.481 0.481 0.525 -0.000 -0.000 -0.000 3/9216 140.0 0.446 0.446 0.486 -0.000 -0.000 -0.000 3/9216 150.0 0.416 0.416 0.453 -0.000 -0.000 -0.000 3/9216 160.0 0.389 0.389 0.424 -0.000 -0.000 -0.000 3/9216 170.0 0.366 0.366 0.399 -0.000 -0.000 -0.000 3/9216 180.0 0.345 0.345 0.376 -0.000 -0.000 -0.000 3/9216 190.0 0.327 0.327 0.356 -0.000 -0.000 -0.000 3/9216 200.0 0.310 0.310 0.338 -0.000 -0.000 -0.000 3/9216 210.0 0.295 0.295 0.322 -0.000 -0.000 -0.000 3/9216 220.0 0.282 0.282 0.307 -0.000 -0.000 -0.000 3/9216 230.0 0.269 0.269 0.294 -0.000 -0.000 -0.000 3/9216 240.0 0.258 0.258 0.282 -0.000 -0.000 -0.000 3/9216 250.0 0.248 0.248 0.270 -0.000 -0.000 -0.000 3/9216 260.0 0.238 0.238 0.260 -0.000 -0.000 -0.000 3/9216 270.0 0.229 0.229 0.250 -0.000 -0.000 -0.000 3/9216 280.0 0.221 0.221 0.241 -0.000 -0.000 -0.000 3/9216 290.0 0.213 0.213 0.233 -0.000 -0.000 -0.000 3/9216 300.0 0.206 0.206 0.225 -0.000 -0.000 -0.000 3/9216 310.0 0.199 0.199 0.218 -0.000 -0.000 -0.000 3/9216 320.0 0.193 0.193 0.211 -0.000 -0.000 -0.000 3/9216 330.0 0.187 0.187 0.205 -0.000 -0.000 -0.000 3/9216 340.0 0.182 0.182 0.199 -0.000 -0.000 -0.000 3/9216 350.0 0.177 0.177 0.193 -0.000 -0.000 -0.000 3/9216 360.0 0.172 0.172 0.187 -0.000 -0.000 -0.000 3/9216 370.0 0.167 0.167 0.182 -0.000 -0.000 -0.000 3/9216 380.0 0.163 0.163 0.178 -0.000 -0.000 -0.000 3/9216 390.0 0.158 0.158 0.173 -0.000 -0.000 -0.000 3/9216 400.0 0.154 0.154 0.169 -0.000 -0.000 -0.000 3/9216 410.0 0.151 0.151 0.165 -0.000 -0.000 -0.000 3/9216 420.0 0.147 0.147 0.161 -0.000 -0.000 -0.000 3/9216 430.0 0.144 0.144 0.157 -0.000 -0.000 -0.000 3/9216 440.0 0.140 0.140 0.153 -0.000 -0.000 -0.000 3/9216 450.0 0.137 0.137 0.150 -0.000 -0.000 -0.000 3/9216 460.0 0.134 0.134 0.147 -0.000 -0.000 -0.000 3/9216 470.0 0.131 0.131 0.144 -0.000 -0.000 -0.000 3/9216 480.0 0.129 0.129 0.141 -0.000 -0.000 -0.000 3/9216 490.0 0.126 0.126 0.138 -0.000 -0.000 -0.000 3/9216 500.0 0.124 0.124 0.135 -0.000 -0.000 -0.000 3/9216 510.0 0.121 0.121 0.132 -0.000 -0.000 -0.000 3/9216 520.0 0.119 0.119 0.130 -0.000 -0.000 -0.000 3/9216 530.0 0.117 0.117 0.127 -0.000 -0.000 -0.000 3/9216 540.0 0.114 0.114 0.125 -0.000 -0.000 -0.000 3/9216 550.0 0.112 0.112 0.123 -0.000 -0.000 -0.000 3/9216 560.0 0.110 0.110 0.120 -0.000 -0.000 -0.000 3/9216 570.0 0.108 0.108 0.118 -0.000 -0.000 -0.000 3/9216 580.0 0.106 0.106 0.116 -0.000 -0.000 -0.000 3/9216 590.0 0.105 0.105 0.114 -0.000 -0.000 -0.000 3/9216 600.0 0.103 0.103 0.112 -0.000 -0.000 -0.000 3/9216 610.0 0.101 0.101 0.111 -0.000 -0.000 -0.000 3/9216 620.0 0.100 0.100 0.109 -0.000 -0.000 -0.000 3/9216 630.0 0.098 0.098 0.107 -0.000 -0.000 -0.000 3/9216 640.0 0.096 0.096 0.105 -0.000 -0.000 -0.000 3/9216 650.0 0.095 0.095 0.104 -0.000 -0.000 -0.000 3/9216 660.0 0.094 0.094 0.102 -0.000 -0.000 -0.000 3/9216 670.0 0.092 0.092 0.101 -0.000 -0.000 -0.000 3/9216 680.0 0.091 0.091 0.099 -0.000 -0.000 -0.000 3/9216 690.0 0.090 0.090 0.098 -0.000 -0.000 -0.000 3/9216 700.0 0.088 0.088 0.096 -0.000 -0.000 -0.000 3/9216 710.0 0.087 0.087 0.095 -0.000 -0.000 -0.000 3/9216 720.0 0.086 0.086 0.094 -0.000 -0.000 -0.000 3/9216 730.0 0.085 0.085 0.092 -0.000 -0.000 -0.000 3/9216 740.0 0.083 0.083 0.091 -0.000 -0.000 -0.000 3/9216 750.0 0.082 0.082 0.090 -0.000 -0.000 -0.000 3/9216 760.0 0.081 0.081 0.089 -0.000 -0.000 -0.000 3/9216 770.0 0.080 0.080 0.088 -0.000 -0.000 -0.000 3/9216 780.0 0.079 0.079 0.087 -0.000 -0.000 -0.000 3/9216 790.0 0.078 0.078 0.085 -0.000 -0.000 -0.000 3/9216 800.0 0.077 0.077 0.084 -0.000 -0.000 -0.000 3/9216 810.0 0.076 0.076 0.083 -0.000 -0.000 -0.000 3/9216 820.0 0.075 0.075 0.082 -0.000 -0.000 -0.000 3/9216 830.0 0.074 0.074 0.081 -0.000 -0.000 -0.000 3/9216 840.0 0.074 0.074 0.080 -0.000 -0.000 -0.000 3/9216 850.0 0.073 0.073 0.079 -0.000 -0.000 -0.000 3/9216 860.0 0.072 0.072 0.078 -0.000 -0.000 -0.000 3/9216 870.0 0.071 0.071 0.078 -0.000 -0.000 -0.000 3/9216 880.0 0.070 0.070 0.077 -0.000 -0.000 -0.000 3/9216 890.0 0.069 0.069 0.076 -0.000 -0.000 -0.000 3/9216 900.0 0.069 0.069 0.075 -0.000 -0.000 -0.000 3/9216 910.0 0.068 0.068 0.074 -0.000 -0.000 -0.000 3/9216 920.0 0.067 0.067 0.073 -0.000 -0.000 -0.000 3/9216 930.0 0.066 0.066 0.073 -0.000 -0.000 -0.000 3/9216 940.0 0.066 0.066 0.072 -0.000 -0.000 -0.000 3/9216 950.0 0.065 0.065 0.071 -0.000 -0.000 -0.000 3/9216 960.0 0.064 0.064 0.070 -0.000 -0.000 -0.000 3/9216 970.0 0.064 0.064 0.070 -0.000 -0.000 -0.000 3/9216 980.0 0.063 0.063 0.069 -0.000 -0.000 -0.000 3/9216 990.0 0.062 0.062 0.068 -0.000 -0.000 -0.000 3/9216 1000.0 0.062 0.062 0.067 -0.000 -0.000 -0.000 3/9216 Thermal conductivity related properties were written into "kappa-m1424.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 23:53:17]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|