----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 04:57:19]------------------------- 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.898739870000000 5.898739870000000 9.688779625000000 b 5.898739870000000 -5.898739870000000 9.688779625000000 c 5.898739870000000 5.898739870000000 -9.688779625000000 Atomic positions (fractional): *1 Rb 0.12500000000000 0.23946474961314 0.61446474961314 85.468 2 Rb 0.98946474961314 0.87500000000000 0.61446474961314 85.468 3 Rb 0.26053525038686 0.37500000000000 0.38553525038686 85.468 4 Rb 0.62500000000000 0.51053525038686 0.38553525038686 85.468 5 Rb 0.62500000000000 0.01053525038686 0.88553525038686 85.468 6 Rb 0.76053525038686 0.37500000000000 0.88553525038686 85.468 7 Rb 0.48946474961314 0.87500000000000 0.11446474961314 85.468 8 Rb 0.12500000000000 0.73946474961314 0.11446474961314 85.468 *9 Rb 0.62883857854968 0.62883857854968 0.75767715709935 85.468 10 Rb 0.37883857854968 0.12116142145032 0.50000000000000 85.468 11 Rb 0.62116142145032 0.87883857854968 0.50000000000000 85.468 12 Rb 0.87116142145032 0.87116142145032 0.24232284290065 85.468 13 Rb 0.12116142145032 0.62116142145032 0.74232284290065 85.468 14 Rb 0.37116142145032 0.12883857854968 1.00000000000000 85.468 15 Rb 0.12883857854968 0.37116142145032 1.00000000000000 85.468 16 Rb 0.87883857854968 0.37883857854968 0.25767715709935 85.468 *17 Pb 0.94117099749323 0.11785508293667 0.94471279471231 207.200 18 Pb 0.86785508293667 0.42314228822437 0.67668408544345 207.200 19 Pb 0.74645820278092 0.19117099749323 0.32331591455655 207.200 20 Pb 0.30882900250677 0.75354179721908 0.67668408544345 207.200 21 Pb 0.07685771177564 0.63214491706333 0.32331591455655 207.200 22 Pb 0.50354179721908 0.32685771177564 0.94471279471231 207.200 23 Pb 0.38214491706333 0.55882900250677 0.05528720528769 207.200 24 Pb 0.17314228822437 0.99645820278092 0.05528720528769 207.200 25 Pb 0.80882900250677 0.13214491706333 0.55528720528769 207.200 26 Pb 0.88214491706333 0.82685771177564 0.82331591455655 207.200 27 Pb 0.00354179721908 0.05882900250677 0.17668408544345 207.200 28 Pb 0.44117099749323 0.49645820278092 0.82331591455655 207.200 29 Pb 0.67314228822436 0.61785508293667 0.17668408544345 207.200 30 Pb 0.24645820278092 0.92314228822436 0.55528720528769 207.200 31 Pb 0.36785508293667 0.69117099749323 0.44471279471231 207.200 32 Pb 0.57685771177564 0.25354179721908 0.44471279471231 207.200 -------------------------------- unit cell --------------------------------- Lattice vectors: a 11.797479740000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.797479740000000 0.000000000000000 c 0.000000000000000 0.000000000000000 19.377559250000001 Atomic positions (fractional): *1 Rb 0.86446474961314 0.75000000000000 0.37500000000000 85.468 > 1 2 Rb 0.75000000000000 0.86446474961314 0.12500000000000 85.468 > 2 3 Rb 0.25000000000000 0.13553525038686 0.12500000000000 85.468 > 3 4 Rb 0.13553525038686 0.25000000000000 0.37500000000000 85.468 > 4 5 Rb 0.63553525038686 0.25000000000000 0.37500000000000 85.468 > 5 6 Rb 0.25000000000000 0.63553525038686 0.12500000000000 85.468 > 6 7 Rb 0.75000000000000 0.36446474961314 0.12500000000000 85.468 > 7 8 Rb 0.36446474961314 0.75000000000000 0.37500000000000 85.468 > 8 *9 Rb 0.37883857854968 0.37883857854968 0.25000000000000 85.468 > 9 10 Rb 0.12116142145032 0.37883857854968 0.00000000000000 85.468 > 10 11 Rb 0.87883857854968 0.62116142145032 0.00000000000000 85.468 > 11 12 Rb 0.62116142145032 0.62116142145032 0.25000000000000 85.468 > 12 13 Rb 0.62116142145032 0.12116142145032 0.00000000000000 85.468 > 13 14 Rb 0.87883857854968 0.12116142145032 0.25000000000000 85.468 > 14 15 Rb 0.12116142145032 0.87883857854968 0.25000000000000 85.468 > 15 16 Rb 0.37883857854968 0.87883857854968 0.00000000000000 85.468 > 16 17 Rb 0.36446474961314 0.25000000000000 0.87500000000000 85.468 > 1 18 Rb 0.25000000000000 0.36446474961314 0.62500000000000 85.468 > 2 19 Rb 0.75000000000000 0.63553525038686 0.62500000000000 85.468 > 3 20 Rb 0.63553525038686 0.75000000000000 0.87500000000000 85.468 > 4 21 Rb 0.13553525038686 0.75000000000000 0.87500000000000 85.468 > 5 22 Rb 0.75000000000000 0.13553525038686 0.62500000000000 85.468 > 6 23 Rb 0.25000000000000 0.86446474961314 0.62500000000000 85.468 > 7 24 Rb 0.86446474961314 0.25000000000000 0.87500000000000 85.468 > 8 25 Rb 0.87883857854968 0.87883857854968 0.75000000000000 85.468 > 9 26 Rb 0.62116142145032 0.87883857854968 0.50000000000000 85.468 > 10 27 Rb 0.37883857854968 0.12116142145032 0.50000000000000 85.468 > 11 28 Rb 0.12116142145032 0.12116142145032 0.75000000000000 85.468 > 12 29 Rb 0.12116142145032 0.62116142145032 0.50000000000000 85.468 > 13 30 Rb 0.37883857854968 0.62116142145032 0.75000000000000 85.468 > 14 31 Rb 0.62116142145032 0.37883857854968 0.75000000000000 85.468 > 15 32 Rb 0.87883857854968 0.37883857854968 0.50000000000000 85.468 > 16 *33 Pb 0.06069844007788 0.88401435463443 0.05715664285880 207.200 > 17 34 Pb 0.11598564536557 0.56069844007788 0.30715664285880 207.200 > 18 35 Pb 0.88401435463443 0.43930155992212 0.30715664285880 207.200 > 19 36 Pb 0.56069844007788 0.11598564536557 0.19284335714120 207.200 > 20 37 Pb 0.43930155992212 0.88401435463443 0.19284335714120 207.200 > 21 38 Pb 0.88401435463443 0.06069844007788 0.44284335714120 207.200 > 22 39 Pb 0.11598564536557 0.93930155992212 0.44284335714120 207.200 > 23 40 Pb 0.93930155992212 0.11598564536557 0.05715664285880 207.200 > 24 41 Pb 0.93930155992212 0.61598564536557 0.19284335714120 207.200 > 25 42 Pb 0.38401435463443 0.43930155992212 0.44284335714120 207.200 > 26 43 Pb 0.61598564536557 0.56069844007788 0.44284335714120 207.200 > 27 44 Pb 0.43930155992212 0.38401435463443 0.05715664285880 207.200 > 28 45 Pb 0.56069844007788 0.61598564536557 0.05715664285880 207.200 > 29 46 Pb 0.61598564536557 0.93930155992212 0.30715664285880 207.200 > 30 47 Pb 0.38401435463443 0.06069844007788 0.30715664285880 207.200 > 31 48 Pb 0.06069844007788 0.38401435463443 0.19284335714120 207.200 > 32 49 Pb 0.56069844007788 0.38401435463443 0.55715664285880 207.200 > 17 50 Pb 0.61598564536557 0.06069844007788 0.80715664285880 207.200 > 18 51 Pb 0.38401435463443 0.93930155992212 0.80715664285880 207.200 > 19 52 Pb 0.06069844007788 0.61598564536557 0.69284335714120 207.200 > 20 53 Pb 0.93930155992212 0.38401435463443 0.69284335714120 207.200 > 21 54 Pb 0.38401435463443 0.56069844007788 0.94284335714120 207.200 > 22 55 Pb 0.61598564536557 0.43930155992212 0.94284335714120 207.200 > 23 56 Pb 0.43930155992212 0.61598564536557 0.55715664285880 207.200 > 24 57 Pb 0.43930155992212 0.11598564536557 0.69284335714120 207.200 > 25 58 Pb 0.88401435463443 0.93930155992212 0.94284335714120 207.200 > 26 59 Pb 0.11598564536557 0.06069844007788 0.94284335714120 207.200 > 27 60 Pb 0.93930155992212 0.88401435463443 0.55715664285880 207.200 > 28 61 Pb 0.06069844007788 0.11598564536557 0.55715664285880 207.200 > 29 62 Pb 0.11598564536557 0.43930155992212 0.80715664285880 207.200 > 30 63 Pb 0.88401435463443 0.56069844007788 0.80715664285880 207.200 > 31 64 Pb 0.56069844007788 0.88401435463443 0.69284335714120 207.200 > 32 -------------------------------- super cell -------------------------------- Lattice vectors: a 11.797479740000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.797479740000000 0.000000000000000 c 0.000000000000000 0.000000000000000 19.377559250000001 Atomic positions (fractional): *1 Rb 0.86446474961314 0.75000000000000 0.37500000000000 85.468 > 1 2 Rb 0.75000000000000 0.86446474961314 0.12500000000000 85.468 > 2 3 Rb 0.25000000000000 0.13553525038686 0.12500000000000 85.468 > 3 4 Rb 0.13553525038686 0.25000000000000 0.37500000000000 85.468 > 4 5 Rb 0.63553525038686 0.25000000000000 0.37500000000000 85.468 > 5 6 Rb 0.25000000000000 0.63553525038686 0.12500000000000 85.468 > 6 7 Rb 0.75000000000000 0.36446474961314 0.12500000000000 85.468 > 7 8 Rb 0.36446474961314 0.75000000000000 0.37500000000000 85.468 > 8 *9 Rb 0.37883857854968 0.37883857854968 0.25000000000000 85.468 > 9 10 Rb 0.12116142145032 0.37883857854968 0.00000000000000 85.468 > 10 11 Rb 0.87883857854968 0.62116142145032 0.00000000000000 85.468 > 11 12 Rb 0.62116142145032 0.62116142145032 0.25000000000000 85.468 > 12 13 Rb 0.62116142145032 0.12116142145032 0.00000000000000 85.468 > 13 14 Rb 0.87883857854968 0.12116142145032 0.25000000000000 85.468 > 14 15 Rb 0.12116142145032 0.87883857854968 0.25000000000000 85.468 > 15 16 Rb 0.37883857854968 0.87883857854968 0.00000000000000 85.468 > 16 17 Rb 0.36446474961314 0.25000000000000 0.87500000000000 85.468 > 1 18 Rb 0.25000000000000 0.36446474961314 0.62500000000000 85.468 > 2 19 Rb 0.75000000000000 0.63553525038686 0.62500000000000 85.468 > 3 20 Rb 0.63553525038686 0.75000000000000 0.87500000000000 85.468 > 4 21 Rb 0.13553525038686 0.75000000000000 0.87500000000000 85.468 > 5 22 Rb 0.75000000000000 0.13553525038686 0.62500000000000 85.468 > 6 23 Rb 0.25000000000000 0.86446474961314 0.62500000000000 85.468 > 7 24 Rb 0.86446474961314 0.25000000000000 0.87500000000000 85.468 > 8 25 Rb 0.87883857854968 0.87883857854968 0.75000000000000 85.468 > 9 26 Rb 0.62116142145032 0.87883857854968 0.50000000000000 85.468 > 10 27 Rb 0.37883857854968 0.12116142145032 0.50000000000000 85.468 > 11 28 Rb 0.12116142145032 0.12116142145032 0.75000000000000 85.468 > 12 29 Rb 0.12116142145032 0.62116142145032 0.50000000000000 85.468 > 13 30 Rb 0.37883857854968 0.62116142145032 0.75000000000000 85.468 > 14 31 Rb 0.62116142145032 0.37883857854968 0.75000000000000 85.468 > 15 32 Rb 0.87883857854968 0.37883857854968 0.50000000000000 85.468 > 16 *33 Pb 0.06069844007788 0.88401435463443 0.05715664285880 207.200 > 17 34 Pb 0.11598564536557 0.56069844007788 0.30715664285880 207.200 > 18 35 Pb 0.88401435463443 0.43930155992212 0.30715664285880 207.200 > 19 36 Pb 0.56069844007788 0.11598564536557 0.19284335714120 207.200 > 20 37 Pb 0.43930155992212 0.88401435463443 0.19284335714120 207.200 > 21 38 Pb 0.88401435463443 0.06069844007788 0.44284335714120 207.200 > 22 39 Pb 0.11598564536557 0.93930155992212 0.44284335714120 207.200 > 23 40 Pb 0.93930155992212 0.11598564536557 0.05715664285880 207.200 > 24 41 Pb 0.93930155992212 0.61598564536557 0.19284335714120 207.200 > 25 42 Pb 0.38401435463443 0.43930155992212 0.44284335714120 207.200 > 26 43 Pb 0.61598564536557 0.56069844007788 0.44284335714120 207.200 > 27 44 Pb 0.43930155992212 0.38401435463443 0.05715664285880 207.200 > 28 45 Pb 0.56069844007788 0.61598564536557 0.05715664285880 207.200 > 29 46 Pb 0.61598564536557 0.93930155992212 0.30715664285880 207.200 > 30 47 Pb 0.38401435463443 0.06069844007788 0.30715664285880 207.200 > 31 48 Pb 0.06069844007788 0.38401435463443 0.19284335714120 207.200 > 32 49 Pb 0.56069844007788 0.38401435463443 0.55715664285880 207.200 > 17 50 Pb 0.61598564536557 0.06069844007788 0.80715664285880 207.200 > 18 51 Pb 0.38401435463443 0.93930155992212 0.80715664285880 207.200 > 19 52 Pb 0.06069844007788 0.61598564536557 0.69284335714120 207.200 > 20 53 Pb 0.93930155992212 0.38401435463443 0.69284335714120 207.200 > 21 54 Pb 0.38401435463443 0.56069844007788 0.94284335714120 207.200 > 22 55 Pb 0.61598564536557 0.43930155992212 0.94284335714120 207.200 > 23 56 Pb 0.43930155992212 0.61598564536557 0.55715664285880 207.200 > 24 57 Pb 0.43930155992212 0.11598564536557 0.69284335714120 207.200 > 25 58 Pb 0.88401435463443 0.93930155992212 0.94284335714120 207.200 > 26 59 Pb 0.11598564536557 0.06069844007788 0.94284335714120 207.200 > 27 60 Pb 0.93930155992212 0.88401435463443 0.55715664285880 207.200 > 28 61 Pb 0.06069844007788 0.11598564536557 0.55715664285880 207.200 > 29 62 Pb 0.11598564536557 0.43930155992212 0.80715664285880 207.200 > 30 63 Pb 0.88401435463443 0.56069844007788 0.80715664285880 207.200 > 31 64 Pb 0.56069844007788 0.88401435463443 0.69284335714120 207.200 > 32 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 14.0257037 0.0000000 0.0000000 0.0000000 14.0257037 0.0000000 0.0000000 0.0000000 15.1745652 -------------------------- Born effective charges -------------------------- 1 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 0.3420357 0.0000000 0.4232832 1.0357203 2 Rb 1.4554383 0.0000000 -0.3420357 0.0000000 1.2627091 0.0000000 -0.4232832 0.0000000 1.0357203 3 Rb 1.4554383 0.0000000 0.3420357 0.0000000 1.2627091 0.0000000 0.4232832 0.0000000 1.0357203 4 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 -0.3420357 0.0000000 -0.4232832 1.0357203 5 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 0.3420357 0.0000000 0.4232832 1.0357203 6 Rb 1.4554383 0.0000000 -0.3420357 0.0000000 1.2627091 0.0000000 -0.4232832 0.0000000 1.0357203 7 Rb 1.4554383 0.0000000 0.3420357 0.0000000 1.2627091 0.0000000 0.4232832 0.0000000 1.0357203 8 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 -0.3420357 0.0000000 -0.4232832 1.0357203 9 Rb 1.1189246 -0.3312070 0.1224567 -0.3312070 1.1189246 -0.1224567 -0.0314029 0.0314029 1.4439728 10 Rb 1.1189246 0.3312070 0.1224567 0.3312070 1.1189246 0.1224567 -0.0314029 -0.0314029 1.4439728 11 Rb 1.1189246 0.3312070 -0.1224567 0.3312070 1.1189246 -0.1224567 0.0314029 0.0314029 1.4439728 12 Rb 1.1189246 -0.3312070 -0.1224567 -0.3312070 1.1189246 0.1224567 0.0314029 -0.0314029 1.4439728 13 Rb 1.1189246 -0.3312070 0.1224567 -0.3312070 1.1189246 -0.1224567 -0.0314029 0.0314029 1.4439728 14 Rb 1.1189246 0.3312070 0.1224567 0.3312070 1.1189246 0.1224567 -0.0314029 -0.0314029 1.4439728 15 Rb 1.1189246 0.3312070 -0.1224567 0.3312070 1.1189246 -0.1224567 0.0314029 0.0314029 1.4439728 16 Rb 1.1189246 -0.3312070 -0.1224567 -0.3312070 1.1189246 0.1224567 0.0314029 -0.0314029 1.4439728 17 Pb -1.0571990 0.9568053 -0.9396816 0.8536039 -1.4207993 0.6328732 -0.8773152 0.6957805 -1.2398466 18 Pb -1.4207993 -0.8536039 -0.6328732 -0.9568053 -1.0571990 -0.9396816 -0.6957805 -0.8773152 -1.2398466 19 Pb -1.4207993 -0.8536039 0.6328732 -0.9568053 -1.0571990 0.9396816 0.6957805 0.8773152 -1.2398466 20 Pb -1.0571990 -0.9568053 0.9396816 -0.8536039 -1.4207993 0.6328732 0.8773152 0.6957805 -1.2398466 21 Pb -1.0571990 -0.9568053 -0.9396816 -0.8536039 -1.4207993 -0.6328732 -0.8773152 -0.6957805 -1.2398466 22 Pb -1.4207993 0.8536039 -0.6328732 0.9568053 -1.0571990 0.9396816 -0.6957805 0.8773152 -1.2398466 23 Pb -1.4207993 0.8536039 0.6328732 0.9568053 -1.0571990 -0.9396816 0.6957805 -0.8773152 -1.2398466 24 Pb -1.0571990 0.9568053 0.9396816 0.8536039 -1.4207993 -0.6328732 0.8773152 -0.6957805 -1.2398466 25 Pb -1.0571990 0.9568053 -0.9396816 0.8536039 -1.4207993 0.6328732 -0.8773152 0.6957805 -1.2398466 26 Pb -1.4207993 -0.8536039 -0.6328732 -0.9568053 -1.0571990 -0.9396816 -0.6957805 -0.8773152 -1.2398466 27 Pb -1.4207993 -0.8536039 0.6328732 -0.9568053 -1.0571990 0.9396816 0.6957805 0.8773152 -1.2398466 28 Pb -1.0571990 -0.9568053 0.9396816 -0.8536039 -1.4207993 0.6328732 0.8773152 0.6957805 -1.2398466 29 Pb -1.0571990 -0.9568053 -0.9396816 -0.8536039 -1.4207993 -0.6328732 -0.8773152 -0.6957805 -1.2398466 30 Pb -1.4207993 0.8536039 -0.6328732 0.9568053 -1.0571990 0.9396816 -0.6957805 0.8773152 -1.2398466 31 Pb -1.4207993 0.8536039 0.6328732 0.9568053 -1.0571990 -0.9396816 0.6957805 -0.8773152 -1.2398466 32 Pb -1.0571990 0.9568053 0.9396816 0.8536039 -1.4207993 -0.6328732 0.8773152 -0.6957805 -1.2398466 ---------------------------------------------------------------------------- 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.345 Solver_block: 200 / 240 - Time: 0.349 Solver_block: 240 / 240 - Time: 0.161 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.868 --------------------------------- 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 (xx) 0.00000000 (xx) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 04:57:26]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 04:57:26]------------------------- 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.898739870000000 5.898739870000000 9.688779625000000 b 5.898739870000000 -5.898739870000000 9.688779625000000 c 5.898739870000000 5.898739870000000 -9.688779625000000 Atomic positions (fractional): 1 Rb 0.12500000000000 0.23946474961314 0.61446474961314 85.468 2 Rb 0.98946474961314 0.87500000000000 0.61446474961314 85.468 3 Rb 0.26053525038686 0.37500000000000 0.38553525038686 85.468 4 Rb 0.62500000000000 0.51053525038686 0.38553525038686 85.468 5 Rb 0.62500000000000 0.01053525038686 0.88553525038686 85.468 6 Rb 0.76053525038686 0.37500000000000 0.88553525038686 85.468 7 Rb 0.48946474961314 0.87500000000000 0.11446474961314 85.468 8 Rb 0.12500000000000 0.73946474961314 0.11446474961314 85.468 9 Rb 0.62883857854968 0.62883857854968 0.75767715709935 85.468 10 Rb 0.37883857854968 0.12116142145032 0.50000000000000 85.468 11 Rb 0.62116142145032 0.87883857854968 0.50000000000000 85.468 12 Rb 0.87116142145032 0.87116142145032 0.24232284290065 85.468 13 Rb 0.12116142145032 0.62116142145032 0.74232284290065 85.468 14 Rb 0.37116142145032 0.12883857854968 1.00000000000000 85.468 15 Rb 0.12883857854968 0.37116142145032 1.00000000000000 85.468 16 Rb 0.87883857854968 0.37883857854968 0.25767715709935 85.468 17 Pb 0.94117099749323 0.11785508293667 0.94471279471231 207.200 18 Pb 0.86785508293667 0.42314228822437 0.67668408544345 207.200 19 Pb 0.74645820278092 0.19117099749323 0.32331591455655 207.200 20 Pb 0.30882900250677 0.75354179721908 0.67668408544345 207.200 21 Pb 0.07685771177564 0.63214491706333 0.32331591455655 207.200 22 Pb 0.50354179721908 0.32685771177564 0.94471279471231 207.200 23 Pb 0.38214491706333 0.55882900250677 0.05528720528769 207.200 24 Pb 0.17314228822437 0.99645820278092 0.05528720528769 207.200 25 Pb 0.80882900250677 0.13214491706333 0.55528720528769 207.200 26 Pb 0.88214491706333 0.82685771177564 0.82331591455655 207.200 27 Pb 0.00354179721908 0.05882900250677 0.17668408544345 207.200 28 Pb 0.44117099749323 0.49645820278092 0.82331591455655 207.200 29 Pb 0.67314228822436 0.61785508293667 0.17668408544345 207.200 30 Pb 0.24645820278092 0.92314228822436 0.55528720528769 207.200 31 Pb 0.36785508293667 0.69117099749323 0.44471279471231 207.200 32 Pb 0.57685771177564 0.25354179721908 0.44471279471231 207.200 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.797479740000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.797479740000000 0.000000000000000 c 0.000000000000000 0.000000000000000 19.377559250000001 Atomic positions (fractional): 1 Rb 0.86446474961314 0.75000000000000 0.37500000000000 85.468 > 1 2 Rb 0.75000000000000 0.86446474961314 0.12500000000000 85.468 > 2 3 Rb 0.25000000000000 0.13553525038686 0.12500000000000 85.468 > 3 4 Rb 0.13553525038686 0.25000000000000 0.37500000000000 85.468 > 4 5 Rb 0.63553525038686 0.25000000000000 0.37500000000000 85.468 > 5 6 Rb 0.25000000000000 0.63553525038686 0.12500000000000 85.468 > 6 7 Rb 0.75000000000000 0.36446474961314 0.12500000000000 85.468 > 7 8 Rb 0.36446474961314 0.75000000000000 0.37500000000000 85.468 > 8 9 Rb 0.37883857854968 0.37883857854968 0.25000000000000 85.468 > 9 10 Rb 0.12116142145032 0.37883857854968 0.00000000000000 85.468 > 10 11 Rb 0.87883857854968 0.62116142145032 0.00000000000000 85.468 > 11 12 Rb 0.62116142145032 0.62116142145032 0.25000000000000 85.468 > 12 13 Rb 0.62116142145032 0.12116142145032 0.00000000000000 85.468 > 13 14 Rb 0.87883857854968 0.12116142145032 0.25000000000000 85.468 > 14 15 Rb 0.12116142145032 0.87883857854968 0.25000000000000 85.468 > 15 16 Rb 0.37883857854968 0.87883857854968 0.00000000000000 85.468 > 16 17 Rb 0.36446474961314 0.25000000000000 0.87500000000000 85.468 > 1 18 Rb 0.25000000000000 0.36446474961314 0.62500000000000 85.468 > 2 19 Rb 0.75000000000000 0.63553525038686 0.62500000000000 85.468 > 3 20 Rb 0.63553525038686 0.75000000000000 0.87500000000000 85.468 > 4 21 Rb 0.13553525038686 0.75000000000000 0.87500000000000 85.468 > 5 22 Rb 0.75000000000000 0.13553525038686 0.62500000000000 85.468 > 6 23 Rb 0.25000000000000 0.86446474961314 0.62500000000000 85.468 > 7 24 Rb 0.86446474961314 0.25000000000000 0.87500000000000 85.468 > 8 25 Rb 0.87883857854968 0.87883857854968 0.75000000000000 85.468 > 9 26 Rb 0.62116142145032 0.87883857854968 0.50000000000000 85.468 > 10 27 Rb 0.37883857854968 0.12116142145032 0.50000000000000 85.468 > 11 28 Rb 0.12116142145032 0.12116142145032 0.75000000000000 85.468 > 12 29 Rb 0.12116142145032 0.62116142145032 0.50000000000000 85.468 > 13 30 Rb 0.37883857854968 0.62116142145032 0.75000000000000 85.468 > 14 31 Rb 0.62116142145032 0.37883857854968 0.75000000000000 85.468 > 15 32 Rb 0.87883857854968 0.37883857854968 0.50000000000000 85.468 > 16 33 Pb 0.06069844007788 0.88401435463443 0.05715664285880 207.200 > 33 34 Pb 0.11598564536557 0.56069844007788 0.30715664285880 207.200 > 34 35 Pb 0.88401435463443 0.43930155992212 0.30715664285880 207.200 > 35 36 Pb 0.56069844007788 0.11598564536557 0.19284335714120 207.200 > 36 37 Pb 0.43930155992212 0.88401435463443 0.19284335714120 207.200 > 37 38 Pb 0.88401435463443 0.06069844007788 0.44284335714120 207.200 > 38 39 Pb 0.11598564536557 0.93930155992212 0.44284335714120 207.200 > 39 40 Pb 0.93930155992212 0.11598564536557 0.05715664285880 207.200 > 40 41 Pb 0.93930155992212 0.61598564536557 0.19284335714120 207.200 > 41 42 Pb 0.38401435463443 0.43930155992212 0.44284335714120 207.200 > 42 43 Pb 0.61598564536557 0.56069844007788 0.44284335714120 207.200 > 43 44 Pb 0.43930155992212 0.38401435463443 0.05715664285880 207.200 > 44 45 Pb 0.56069844007788 0.61598564536557 0.05715664285880 207.200 > 45 46 Pb 0.61598564536557 0.93930155992212 0.30715664285880 207.200 > 46 47 Pb 0.38401435463443 0.06069844007788 0.30715664285880 207.200 > 47 48 Pb 0.06069844007788 0.38401435463443 0.19284335714120 207.200 > 48 49 Pb 0.56069844007788 0.38401435463443 0.55715664285880 207.200 > 33 50 Pb 0.61598564536557 0.06069844007788 0.80715664285880 207.200 > 34 51 Pb 0.38401435463443 0.93930155992212 0.80715664285880 207.200 > 35 52 Pb 0.06069844007788 0.61598564536557 0.69284335714120 207.200 > 36 53 Pb 0.93930155992212 0.38401435463443 0.69284335714120 207.200 > 37 54 Pb 0.38401435463443 0.56069844007788 0.94284335714120 207.200 > 38 55 Pb 0.61598564536557 0.43930155992212 0.94284335714120 207.200 > 39 56 Pb 0.43930155992212 0.61598564536557 0.55715664285880 207.200 > 40 57 Pb 0.43930155992212 0.11598564536557 0.69284335714120 207.200 > 41 58 Pb 0.88401435463443 0.93930155992212 0.94284335714120 207.200 > 42 59 Pb 0.11598564536557 0.06069844007788 0.94284335714120 207.200 > 43 60 Pb 0.93930155992212 0.88401435463443 0.55715664285880 207.200 > 44 61 Pb 0.06069844007788 0.11598564536557 0.55715664285880 207.200 > 45 62 Pb 0.11598564536557 0.43930155992212 0.80715664285880 207.200 > 46 63 Pb 0.88401435463443 0.56069844007788 0.80715664285880 207.200 > 47 64 Pb 0.56069844007788 0.88401435463443 0.69284335714120 207.200 > 48 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 14.0257037 0.0000000 0.0000000 0.0000000 14.0257037 0.0000000 0.0000000 0.0000000 15.1745652 -------------------------- Born effective charges -------------------------- 1 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 0.3420357 0.0000000 0.4232832 1.0357203 2 Rb 1.4554383 0.0000000 -0.3420357 0.0000000 1.2627091 0.0000000 -0.4232832 0.0000000 1.0357203 3 Rb 1.4554383 0.0000000 0.3420357 0.0000000 1.2627091 0.0000000 0.4232832 0.0000000 1.0357203 4 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 -0.3420357 0.0000000 -0.4232832 1.0357203 5 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 0.3420357 0.0000000 0.4232832 1.0357203 6 Rb 1.4554383 0.0000000 -0.3420357 0.0000000 1.2627091 0.0000000 -0.4232832 0.0000000 1.0357203 7 Rb 1.4554383 0.0000000 0.3420357 0.0000000 1.2627091 0.0000000 0.4232832 0.0000000 1.0357203 8 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 -0.3420357 0.0000000 -0.4232832 1.0357203 9 Rb 1.1189246 -0.3312070 0.1224567 -0.3312070 1.1189246 -0.1224567 -0.0314029 0.0314029 1.4439728 10 Rb 1.1189246 0.3312070 0.1224567 0.3312070 1.1189246 0.1224567 -0.0314029 -0.0314029 1.4439728 11 Rb 1.1189246 0.3312070 -0.1224567 0.3312070 1.1189246 -0.1224567 0.0314029 0.0314029 1.4439728 12 Rb 1.1189246 -0.3312070 -0.1224567 -0.3312070 1.1189246 0.1224567 0.0314029 -0.0314029 1.4439728 13 Rb 1.1189246 -0.3312070 0.1224567 -0.3312070 1.1189246 -0.1224567 -0.0314029 0.0314029 1.4439728 14 Rb 1.1189246 0.3312070 0.1224567 0.3312070 1.1189246 0.1224567 -0.0314029 -0.0314029 1.4439728 15 Rb 1.1189246 0.3312070 -0.1224567 0.3312070 1.1189246 -0.1224567 0.0314029 0.0314029 1.4439728 16 Rb 1.1189246 -0.3312070 -0.1224567 -0.3312070 1.1189246 0.1224567 0.0314029 -0.0314029 1.4439728 17 Pb -1.0571990 0.9568053 -0.9396816 0.8536039 -1.4207993 0.6328732 -0.8773152 0.6957805 -1.2398466 18 Pb -1.4207993 -0.8536039 -0.6328732 -0.9568053 -1.0571990 -0.9396816 -0.6957805 -0.8773152 -1.2398466 19 Pb -1.4207993 -0.8536039 0.6328732 -0.9568053 -1.0571990 0.9396816 0.6957805 0.8773152 -1.2398466 20 Pb -1.0571990 -0.9568053 0.9396816 -0.8536039 -1.4207993 0.6328732 0.8773152 0.6957805 -1.2398466 21 Pb -1.0571990 -0.9568053 -0.9396816 -0.8536039 -1.4207993 -0.6328732 -0.8773152 -0.6957805 -1.2398466 22 Pb -1.4207993 0.8536039 -0.6328732 0.9568053 -1.0571990 0.9396816 -0.6957805 0.8773152 -1.2398466 23 Pb -1.4207993 0.8536039 0.6328732 0.9568053 -1.0571990 -0.9396816 0.6957805 -0.8773152 -1.2398466 24 Pb -1.0571990 0.9568053 0.9396816 0.8536039 -1.4207993 -0.6328732 0.8773152 -0.6957805 -1.2398466 25 Pb -1.0571990 0.9568053 -0.9396816 0.8536039 -1.4207993 0.6328732 -0.8773152 0.6957805 -1.2398466 26 Pb -1.4207993 -0.8536039 -0.6328732 -0.9568053 -1.0571990 -0.9396816 -0.6957805 -0.8773152 -1.2398466 27 Pb -1.4207993 -0.8536039 0.6328732 -0.9568053 -1.0571990 0.9396816 0.6957805 0.8773152 -1.2398466 28 Pb -1.0571990 -0.9568053 0.9396816 -0.8536039 -1.4207993 0.6328732 0.8773152 0.6957805 -1.2398466 29 Pb -1.0571990 -0.9568053 -0.9396816 -0.8536039 -1.4207993 -0.6328732 -0.8773152 -0.6957805 -1.2398466 30 Pb -1.4207993 0.8536039 -0.6328732 0.9568053 -1.0571990 0.9396816 -0.6957805 0.8773152 -1.2398466 31 Pb -1.4207993 0.8536039 0.6328732 0.9568053 -1.0571990 -0.9396816 0.6957805 -0.8773152 -1.2398466 32 Pb -1.0571990 0.9568053 0.9396816 0.8536039 -1.4207993 -0.6328732 0.8773152 -0.6957805 -1.2398466 ---------------------------------------------------------------------------- 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 (zy) 0.00000000 (zy) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 04:57:33]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 04:57:33]------------------------- 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.898739870000000 5.898739870000000 9.688779625000000 b 5.898739870000000 -5.898739870000000 9.688779625000000 c 5.898739870000000 5.898739870000000 -9.688779625000000 Atomic positions (fractional): 1 Rb 0.12500000000000 0.23946474961314 0.61446474961314 85.468 2 Rb 0.98946474961314 0.87500000000000 0.61446474961314 85.468 3 Rb 0.26053525038686 0.37500000000000 0.38553525038686 85.468 4 Rb 0.62500000000000 0.51053525038686 0.38553525038686 85.468 5 Rb 0.62500000000000 0.01053525038686 0.88553525038686 85.468 6 Rb 0.76053525038686 0.37500000000000 0.88553525038686 85.468 7 Rb 0.48946474961314 0.87500000000000 0.11446474961314 85.468 8 Rb 0.12500000000000 0.73946474961314 0.11446474961314 85.468 9 Rb 0.62883857854968 0.62883857854968 0.75767715709935 85.468 10 Rb 0.37883857854968 0.12116142145032 0.50000000000000 85.468 11 Rb 0.62116142145032 0.87883857854968 0.50000000000000 85.468 12 Rb 0.87116142145032 0.87116142145032 0.24232284290065 85.468 13 Rb 0.12116142145032 0.62116142145032 0.74232284290065 85.468 14 Rb 0.37116142145032 0.12883857854968 1.00000000000000 85.468 15 Rb 0.12883857854968 0.37116142145032 1.00000000000000 85.468 16 Rb 0.87883857854968 0.37883857854968 0.25767715709935 85.468 17 Pb 0.94117099749323 0.11785508293667 0.94471279471231 207.200 18 Pb 0.86785508293667 0.42314228822437 0.67668408544345 207.200 19 Pb 0.74645820278092 0.19117099749323 0.32331591455655 207.200 20 Pb 0.30882900250677 0.75354179721908 0.67668408544345 207.200 21 Pb 0.07685771177564 0.63214491706333 0.32331591455655 207.200 22 Pb 0.50354179721908 0.32685771177564 0.94471279471231 207.200 23 Pb 0.38214491706333 0.55882900250677 0.05528720528769 207.200 24 Pb 0.17314228822437 0.99645820278092 0.05528720528769 207.200 25 Pb 0.80882900250677 0.13214491706333 0.55528720528769 207.200 26 Pb 0.88214491706333 0.82685771177564 0.82331591455655 207.200 27 Pb 0.00354179721908 0.05882900250677 0.17668408544345 207.200 28 Pb 0.44117099749323 0.49645820278092 0.82331591455655 207.200 29 Pb 0.67314228822436 0.61785508293667 0.17668408544345 207.200 30 Pb 0.24645820278092 0.92314228822436 0.55528720528769 207.200 31 Pb 0.36785508293667 0.69117099749323 0.44471279471231 207.200 32 Pb 0.57685771177564 0.25354179721908 0.44471279471231 207.200 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.797479740000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.797479740000000 0.000000000000000 c 0.000000000000000 0.000000000000000 19.377559250000001 Atomic positions (fractional): 1 Rb 0.86446474961314 0.75000000000000 0.37500000000000 85.468 > 1 2 Rb 0.75000000000000 0.86446474961314 0.12500000000000 85.468 > 2 3 Rb 0.25000000000000 0.13553525038686 0.12500000000000 85.468 > 3 4 Rb 0.13553525038686 0.25000000000000 0.37500000000000 85.468 > 4 5 Rb 0.63553525038686 0.25000000000000 0.37500000000000 85.468 > 5 6 Rb 0.25000000000000 0.63553525038686 0.12500000000000 85.468 > 6 7 Rb 0.75000000000000 0.36446474961314 0.12500000000000 85.468 > 7 8 Rb 0.36446474961314 0.75000000000000 0.37500000000000 85.468 > 8 9 Rb 0.37883857854968 0.37883857854968 0.25000000000000 85.468 > 9 10 Rb 0.12116142145032 0.37883857854968 0.00000000000000 85.468 > 10 11 Rb 0.87883857854968 0.62116142145032 0.00000000000000 85.468 > 11 12 Rb 0.62116142145032 0.62116142145032 0.25000000000000 85.468 > 12 13 Rb 0.62116142145032 0.12116142145032 0.00000000000000 85.468 > 13 14 Rb 0.87883857854968 0.12116142145032 0.25000000000000 85.468 > 14 15 Rb 0.12116142145032 0.87883857854968 0.25000000000000 85.468 > 15 16 Rb 0.37883857854968 0.87883857854968 0.00000000000000 85.468 > 16 17 Rb 0.36446474961314 0.25000000000000 0.87500000000000 85.468 > 1 18 Rb 0.25000000000000 0.36446474961314 0.62500000000000 85.468 > 2 19 Rb 0.75000000000000 0.63553525038686 0.62500000000000 85.468 > 3 20 Rb 0.63553525038686 0.75000000000000 0.87500000000000 85.468 > 4 21 Rb 0.13553525038686 0.75000000000000 0.87500000000000 85.468 > 5 22 Rb 0.75000000000000 0.13553525038686 0.62500000000000 85.468 > 6 23 Rb 0.25000000000000 0.86446474961314 0.62500000000000 85.468 > 7 24 Rb 0.86446474961314 0.25000000000000 0.87500000000000 85.468 > 8 25 Rb 0.87883857854968 0.87883857854968 0.75000000000000 85.468 > 9 26 Rb 0.62116142145032 0.87883857854968 0.50000000000000 85.468 > 10 27 Rb 0.37883857854968 0.12116142145032 0.50000000000000 85.468 > 11 28 Rb 0.12116142145032 0.12116142145032 0.75000000000000 85.468 > 12 29 Rb 0.12116142145032 0.62116142145032 0.50000000000000 85.468 > 13 30 Rb 0.37883857854968 0.62116142145032 0.75000000000000 85.468 > 14 31 Rb 0.62116142145032 0.37883857854968 0.75000000000000 85.468 > 15 32 Rb 0.87883857854968 0.37883857854968 0.50000000000000 85.468 > 16 33 Pb 0.06069844007788 0.88401435463443 0.05715664285880 207.200 > 33 34 Pb 0.11598564536557 0.56069844007788 0.30715664285880 207.200 > 34 35 Pb 0.88401435463443 0.43930155992212 0.30715664285880 207.200 > 35 36 Pb 0.56069844007788 0.11598564536557 0.19284335714120 207.200 > 36 37 Pb 0.43930155992212 0.88401435463443 0.19284335714120 207.200 > 37 38 Pb 0.88401435463443 0.06069844007788 0.44284335714120 207.200 > 38 39 Pb 0.11598564536557 0.93930155992212 0.44284335714120 207.200 > 39 40 Pb 0.93930155992212 0.11598564536557 0.05715664285880 207.200 > 40 41 Pb 0.93930155992212 0.61598564536557 0.19284335714120 207.200 > 41 42 Pb 0.38401435463443 0.43930155992212 0.44284335714120 207.200 > 42 43 Pb 0.61598564536557 0.56069844007788 0.44284335714120 207.200 > 43 44 Pb 0.43930155992212 0.38401435463443 0.05715664285880 207.200 > 44 45 Pb 0.56069844007788 0.61598564536557 0.05715664285880 207.200 > 45 46 Pb 0.61598564536557 0.93930155992212 0.30715664285880 207.200 > 46 47 Pb 0.38401435463443 0.06069844007788 0.30715664285880 207.200 > 47 48 Pb 0.06069844007788 0.38401435463443 0.19284335714120 207.200 > 48 49 Pb 0.56069844007788 0.38401435463443 0.55715664285880 207.200 > 33 50 Pb 0.61598564536557 0.06069844007788 0.80715664285880 207.200 > 34 51 Pb 0.38401435463443 0.93930155992212 0.80715664285880 207.200 > 35 52 Pb 0.06069844007788 0.61598564536557 0.69284335714120 207.200 > 36 53 Pb 0.93930155992212 0.38401435463443 0.69284335714120 207.200 > 37 54 Pb 0.38401435463443 0.56069844007788 0.94284335714120 207.200 > 38 55 Pb 0.61598564536557 0.43930155992212 0.94284335714120 207.200 > 39 56 Pb 0.43930155992212 0.61598564536557 0.55715664285880 207.200 > 40 57 Pb 0.43930155992212 0.11598564536557 0.69284335714120 207.200 > 41 58 Pb 0.88401435463443 0.93930155992212 0.94284335714120 207.200 > 42 59 Pb 0.11598564536557 0.06069844007788 0.94284335714120 207.200 > 43 60 Pb 0.93930155992212 0.88401435463443 0.55715664285880 207.200 > 44 61 Pb 0.06069844007788 0.11598564536557 0.55715664285880 207.200 > 45 62 Pb 0.11598564536557 0.43930155992212 0.80715664285880 207.200 > 46 63 Pb 0.88401435463443 0.56069844007788 0.80715664285880 207.200 > 47 64 Pb 0.56069844007788 0.88401435463443 0.69284335714120 207.200 > 48 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 14.0257037 0.0000000 0.0000000 0.0000000 14.0257037 0.0000000 0.0000000 0.0000000 15.1745652 -------------------------- Born effective charges -------------------------- 1 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 0.3420357 0.0000000 0.4232832 1.0357203 2 Rb 1.4554383 0.0000000 -0.3420357 0.0000000 1.2627091 0.0000000 -0.4232832 0.0000000 1.0357203 3 Rb 1.4554383 0.0000000 0.3420357 0.0000000 1.2627091 0.0000000 0.4232832 0.0000000 1.0357203 4 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 -0.3420357 0.0000000 -0.4232832 1.0357203 5 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 0.3420357 0.0000000 0.4232832 1.0357203 6 Rb 1.4554383 0.0000000 -0.3420357 0.0000000 1.2627091 0.0000000 -0.4232832 0.0000000 1.0357203 7 Rb 1.4554383 0.0000000 0.3420357 0.0000000 1.2627091 0.0000000 0.4232832 0.0000000 1.0357203 8 Rb 1.2627091 0.0000000 0.0000000 0.0000000 1.4554383 -0.3420357 0.0000000 -0.4232832 1.0357203 9 Rb 1.1189246 -0.3312070 0.1224567 -0.3312070 1.1189246 -0.1224567 -0.0314029 0.0314029 1.4439728 10 Rb 1.1189246 0.3312070 0.1224567 0.3312070 1.1189246 0.1224567 -0.0314029 -0.0314029 1.4439728 11 Rb 1.1189246 0.3312070 -0.1224567 0.3312070 1.1189246 -0.1224567 0.0314029 0.0314029 1.4439728 12 Rb 1.1189246 -0.3312070 -0.1224567 -0.3312070 1.1189246 0.1224567 0.0314029 -0.0314029 1.4439728 13 Rb 1.1189246 -0.3312070 0.1224567 -0.3312070 1.1189246 -0.1224567 -0.0314029 0.0314029 1.4439728 14 Rb 1.1189246 0.3312070 0.1224567 0.3312070 1.1189246 0.1224567 -0.0314029 -0.0314029 1.4439728 15 Rb 1.1189246 0.3312070 -0.1224567 0.3312070 1.1189246 -0.1224567 0.0314029 0.0314029 1.4439728 16 Rb 1.1189246 -0.3312070 -0.1224567 -0.3312070 1.1189246 0.1224567 0.0314029 -0.0314029 1.4439728 17 Pb -1.0571990 0.9568053 -0.9396816 0.8536039 -1.4207993 0.6328732 -0.8773152 0.6957805 -1.2398466 18 Pb -1.4207993 -0.8536039 -0.6328732 -0.9568053 -1.0571990 -0.9396816 -0.6957805 -0.8773152 -1.2398466 19 Pb -1.4207993 -0.8536039 0.6328732 -0.9568053 -1.0571990 0.9396816 0.6957805 0.8773152 -1.2398466 20 Pb -1.0571990 -0.9568053 0.9396816 -0.8536039 -1.4207993 0.6328732 0.8773152 0.6957805 -1.2398466 21 Pb -1.0571990 -0.9568053 -0.9396816 -0.8536039 -1.4207993 -0.6328732 -0.8773152 -0.6957805 -1.2398466 22 Pb -1.4207993 0.8536039 -0.6328732 0.9568053 -1.0571990 0.9396816 -0.6957805 0.8773152 -1.2398466 23 Pb -1.4207993 0.8536039 0.6328732 0.9568053 -1.0571990 -0.9396816 0.6957805 -0.8773152 -1.2398466 24 Pb -1.0571990 0.9568053 0.9396816 0.8536039 -1.4207993 -0.6328732 0.8773152 -0.6957805 -1.2398466 25 Pb -1.0571990 0.9568053 -0.9396816 0.8536039 -1.4207993 0.6328732 -0.8773152 0.6957805 -1.2398466 26 Pb -1.4207993 -0.8536039 -0.6328732 -0.9568053 -1.0571990 -0.9396816 -0.6957805 -0.8773152 -1.2398466 27 Pb -1.4207993 -0.8536039 0.6328732 -0.9568053 -1.0571990 0.9396816 0.6957805 0.8773152 -1.2398466 28 Pb -1.0571990 -0.9568053 0.9396816 -0.8536039 -1.4207993 0.6328732 0.8773152 0.6957805 -1.2398466 29 Pb -1.0571990 -0.9568053 -0.9396816 -0.8536039 -1.4207993 -0.6328732 -0.8773152 -0.6957805 -1.2398466 30 Pb -1.4207993 0.8536039 -0.6328732 0.9568053 -1.0571990 0.9396816 -0.6957805 0.8773152 -1.2398466 31 Pb -1.4207993 0.8536039 0.6328732 0.9568053 -1.0571990 -0.9396816 0.6957805 -0.8773152 -1.2398466 32 Pb -1.0571990 0.9568053 0.9396816 0.8536039 -1.4207993 -0.6328732 0.8773152 -0.6957805 -1.2398466 ---------------------------------------------------------------------------- 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 (xx) 0.00000000 (xx) --------------------------- 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.15 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 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.671 ( -0.000 0.000 -0.000) 0.000 0.671 ( 0.000 0.000 0.000) 0.000 0.706 ( -0.000 0.000 -0.000) 0.000 0.706 ( 0.000 0.000 0.000) 0.000 0.814 ( -0.000 0.000 0.000) 0.000 0.819 ( 0.000 -0.000 -0.000) 0.000 0.819 ( 0.000 0.000 0.000) 0.000 0.852 ( 0.000 -0.000 0.000) 0.000 0.862 ( 0.000 -0.000 -0.000) 0.000 0.862 ( 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.899 ( 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.915 ( 0.000 -0.000 0.000) 0.000 0.948 ( 0.000 0.000 0.000) 0.000 0.952 ( -0.000 0.000 -0.000) 0.000 0.952 ( 0.000 0.000 0.000) 0.000 0.962 ( 0.000 0.000 -0.000) 0.000 0.997 ( 0.000 -0.000 0.000) 0.000 1.203 ( 0.000 0.000 0.000) 0.000 1.203 ( 0.000 0.000 -0.000) 0.000 1.252 ( 0.000 -0.000 0.000) 0.000 1.252 ( 0.000 0.000 -0.000) 0.000 1.284 ( 0.000 0.000 0.000) 0.000 1.291 ( 0.000 0.000 0.000) 0.000 1.302 ( 0.000 0.000 -0.000) 0.000 1.305 ( 0.000 -0.000 -0.000) 0.000 1.313 ( 0.000 -0.000 0.000) 0.000 1.313 ( -0.000 0.000 0.000) 0.000 1.465 ( -0.000 0.000 -0.000) 0.000 1.508 ( 0.000 -0.000 0.000) 0.000 1.561 ( -0.000 0.000 -0.000) 0.000 1.561 ( 0.000 0.000 0.000) 0.000 1.580 ( 0.000 -0.000 0.000) 0.000 1.612 ( 0.000 0.000 -0.000) 0.000 1.623 ( 0.000 0.000 0.000) 0.000 1.623 ( 0.000 0.000 0.000) 0.000 1.663 ( 0.000 0.000 0.000) 0.000 1.707 ( 0.000 0.000 -0.000) 0.000 1.728 ( -0.000 0.000 -0.000) 0.000 1.731 ( 0.000 0.000 0.000) 0.000 1.761 ( 0.000 -0.000 0.000) 0.000 1.800 ( 0.000 -0.000 0.000) 0.000 1.801 ( 0.000 -0.000 0.000) 0.000 1.801 ( 0.000 0.000 0.000) 0.000 1.849 ( 0.000 0.000 0.000) 0.000 1.849 ( 0.000 -0.000 -0.000) 0.000 1.891 ( 0.000 0.000 0.000) 0.000 1.933 ( 0.000 -0.000 0.000) 0.000 1.933 ( -0.000 0.000 -0.000) 0.000 1.940 ( -0.000 0.000 0.000) 0.000 1.957 ( -0.000 0.000 -0.000) 0.000 1.977 ( 0.000 0.000 0.000) 0.000 2.053 ( 0.000 0.000 0.000) 0.000 2.060 ( -0.000 0.000 0.000) 0.000 2.060 ( -0.000 0.000 -0.000) 0.000 2.085 ( 0.000 -0.000 -0.000) 0.000 2.096 ( 0.000 0.000 0.000) 0.000 2.096 ( 0.000 0.000 -0.000) 0.000 2.126 ( 0.000 0.000 0.000) 0.000 2.167 ( 0.000 0.000 -0.000) 0.000 2.167 ( 0.000 -0.000 -0.000) 0.000 2.200 ( 0.000 -0.000 0.000) 0.000 2.200 ( 0.000 0.000 0.000) 0.000 2.295 ( 0.000 -0.000 0.000) 0.000 2.313 ( -0.000 0.000 -0.000) 0.000 2.319 ( -0.000 0.000 0.000) 0.000 2.330 ( 0.000 0.000 0.000) 0.000 2.330 ( 0.000 0.000 -0.000) 0.000 2.339 ( 0.000 0.000 0.000) 0.000 2.356 ( 0.000 -0.000 0.000) 0.000 2.429 ( 0.000 -0.000 -0.000) 0.000 2.442 ( 0.000 0.000 0.000) 0.000 2.457 ( 0.000 0.000 0.000) 0.000 2.478 ( 0.000 0.000 -0.000) 0.000 2.699 ( 0.000 -0.000 -0.000) 0.000 2.699 ( -0.000 0.000 0.000) 0.000 2.721 ( 0.000 0.000 0.000) 0.000 2.721 ( 0.000 0.000 -0.000) 0.000 2.856 ( 0.000 0.000 0.000) 0.000 2.856 ( -0.000 0.000 -0.000) 0.000 2.861 ( 0.000 0.000 0.000) 0.000 2.861 ( 0.000 0.000 0.000) 0.000 2.874 ( 0.000 0.000 0.000) 0.000 2.878 ( 0.000 0.000 0.000) 0.000 2.887 ( 0.000 0.000 -0.000) 0.000 2.917 ( 0.000 -0.000 -0.000) 0.000 3.850 ( 0.000 0.000 0.000) 0.000 3.889 ( 0.000 0.000 0.000) 0.000 3.890 ( 0.000 0.000 0.000) 0.000 3.903 ( -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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.306 ( 6.615 6.615 0.000) 9.355 0.314 ( 6.530 6.530 0.000) 9.235 0.559 ( 11.025 11.025 0.000) 15.592 0.665 ( -0.525 -0.525 0.000) 0.742 0.669 ( -1.644 -1.644 0.000) 2.325 0.698 ( 1.061 1.061 0.000) 1.501 0.726 ( -4.815 -4.815 0.000) 6.809 0.729 ( -3.708 -3.708 0.000) 5.244 0.744 ( 1.024 1.024 0.000) 1.448 0.813 ( -0.533 -0.533 0.000) 0.753 0.817 ( -1.288 -1.288 0.000) 1.822 0.829 ( 0.416 0.416 0.000) 0.589 0.841 ( -1.698 -1.698 0.000) 2.401 0.853 ( 0.259 0.259 0.000) 0.366 0.867 ( 0.031 0.031 0.000) 0.043 0.876 ( 0.211 0.211 0.000) 0.298 0.903 ( -0.413 -0.413 0.000) 0.585 0.904 ( -0.529 -0.529 0.000) 0.747 0.904 ( -0.861 -0.861 0.000) 1.217 0.928 ( 1.199 1.199 0.000) 1.696 0.940 ( 0.305 0.305 0.000) 0.432 0.971 ( -0.519 -0.519 0.000) 0.734 0.974 ( -0.240 -0.240 0.000) 0.340 1.004 ( -0.496 -0.496 0.000) 0.702 1.149 ( -0.395 -0.395 0.000) 0.559 1.201 ( -1.042 -1.042 0.000) 1.474 1.241 ( 1.157 1.157 0.000) 1.636 1.273 ( 0.531 0.531 0.000) 0.751 1.307 ( 0.220 0.220 0.000) 0.311 1.322 ( 0.927 0.927 0.000) 1.311 1.327 ( 0.287 0.287 0.000) 0.407 1.331 ( 0.263 0.263 0.000) 0.372 1.343 ( 1.549 1.549 0.000) 2.191 1.403 ( 2.612 2.612 0.000) 3.694 1.468 ( -3.133 -3.133 0.000) 4.431 1.472 ( -0.162 -0.162 0.000) 0.229 1.485 ( -1.115 -1.115 0.000) 1.577 1.575 ( 0.426 0.426 0.000) 0.602 1.591 ( -1.574 -1.574 0.000) 2.226 1.593 ( -0.750 -0.750 0.000) 1.061 1.614 ( 0.624 0.624 0.000) 0.883 1.688 ( 1.254 1.254 0.000) 1.773 1.690 ( -0.093 -0.093 0.000) 0.132 1.690 ( -1.295 -1.295 0.000) 1.832 1.724 ( 1.167 1.167 0.000) 1.650 1.760 ( 0.393 0.393 0.000) 0.556 1.768 ( -1.492 -1.492 0.000) 2.110 1.788 ( -1.421 -1.421 0.000) 2.010 1.792 ( 1.707 1.707 0.000) 2.415 1.815 ( 2.267 2.267 0.000) 3.206 1.836 ( -1.217 -1.217 0.000) 1.721 1.838 ( 0.423 0.423 0.000) 0.598 1.860 ( -0.314 -0.314 0.000) 0.443 1.865 ( 0.404 0.404 0.000) 0.571 1.913 ( -4.286 -4.286 0.000) 6.062 1.981 ( -3.285 -3.285 0.000) 4.645 1.992 ( 0.836 0.836 0.000) 1.182 1.995 ( 0.268 0.268 0.000) 0.379 1.998 ( -1.797 -1.797 0.000) 2.541 2.034 ( 0.453 0.453 0.000) 0.641 2.041 ( -2.312 -2.312 0.000) 3.270 2.058 ( -0.838 -0.838 0.000) 1.186 2.099 ( 0.943 0.943 0.000) 1.333 2.111 ( 1.410 1.410 0.000) 1.994 2.182 ( 1.103 1.103 0.000) 1.560 2.189 ( 1.571 1.571 0.000) 2.222 2.209 ( 2.190 2.190 0.000) 3.097 2.212 ( 1.101 1.101 0.000) 1.557 2.276 ( -0.392 -0.392 0.000) 0.554 2.297 ( -0.952 -0.952 0.000) 1.347 2.314 ( 1.309 1.309 0.000) 1.852 2.325 ( 0.139 0.139 0.000) 0.197 2.351 ( -0.010 -0.010 0.000) 0.014 2.352 ( 0.368 0.368 0.000) 0.521 2.361 ( 0.642 0.642 0.000) 0.907 2.370 ( 0.818 0.818 0.000) 1.157 2.436 ( 0.337 0.337 0.000) 0.476 2.445 ( -0.732 -0.732 0.000) 1.035 2.462 ( 0.546 0.546 0.000) 0.773 2.470 ( -0.147 -0.147 0.000) 0.208 2.708 ( 0.399 0.399 0.000) 0.564 2.729 ( 0.226 0.226 0.000) 0.319 2.742 ( 0.648 0.648 0.000) 0.917 2.744 ( 0.192 0.192 0.000) 0.272 2.825 ( -1.291 -1.291 0.000) 1.825 2.825 ( -1.378 -1.378 0.000) 1.949 2.844 ( -0.355 -0.355 0.000) 0.503 2.854 ( -0.385 -0.385 0.000) 0.545 2.876 ( 0.123 0.123 0.000) 0.174 2.895 ( 0.566 0.566 0.000) 0.800 2.901 ( 0.406 0.406 0.000) 0.574 2.906 ( -0.403 -0.403 0.000) 0.570 3.859 ( 0.394 0.394 0.000) 0.558 3.887 ( -0.070 -0.070 0.000) 0.099 3.888 ( -0.123 -0.123 0.000) 0.174 3.896 ( -0.257 -0.257 0.000) 0.364 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 16 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.490 ( -0.096 -0.096 0.000) 0.136 0.490 ( 0.096 0.096 0.000) 0.136 0.595 ( -0.552 -0.552 0.000) 0.780 0.595 ( 0.552 0.552 0.000) 0.780 0.632 ( -0.587 -0.587 0.000) 0.830 0.632 ( 0.587 0.587 0.000) 0.830 0.752 ( -1.102 -1.102 0.000) 1.559 0.752 ( 1.102 1.102 0.000) 1.559 0.775 ( -0.828 -0.828 0.000) 1.171 0.775 ( 0.828 0.828 0.000) 1.171 0.788 ( -0.401 -0.401 0.000) 0.567 0.788 ( 0.401 0.401 0.000) 0.567 0.846 ( -0.731 -0.731 0.000) 1.034 0.846 ( 0.731 0.731 0.000) 1.034 0.884 ( -0.243 -0.243 0.000) 0.344 0.884 ( 0.243 0.243 0.000) 0.344 0.888 ( -0.184 -0.184 0.000) 0.261 0.888 ( 0.184 0.184 0.000) 0.261 0.953 ( -0.696 -0.696 0.000) 0.985 0.953 ( 0.696 0.696 0.000) 0.985 0.962 ( -0.909 -0.909 0.000) 1.285 0.962 ( 0.909 0.909 0.000) 1.285 0.976 ( -0.623 -0.623 0.000) 0.881 0.976 ( 0.623 0.623 0.000) 0.881 1.164 ( -0.727 -0.727 0.000) 1.028 1.164 ( 0.727 0.727 0.000) 1.028 1.248 ( -0.706 -0.706 0.000) 0.999 1.248 ( 0.706 0.706 0.000) 0.999 1.322 ( -0.372 -0.372 0.000) 0.526 1.322 ( 0.372 0.372 0.000) 0.526 1.341 ( -0.220 -0.220 0.000) 0.311 1.341 ( 0.220 0.220 0.000) 0.311 1.424 ( -1.833 -1.833 0.000) 2.592 1.424 ( 1.833 1.833 0.000) 2.592 1.446 ( -1.276 -1.276 0.000) 1.804 1.446 ( 1.276 1.276 0.000) 1.804 1.497 ( -2.141 -2.141 0.000) 3.028 1.497 ( 2.141 2.141 0.000) 3.028 1.580 ( -0.113 -0.113 0.000) 0.160 1.580 ( 0.113 0.113 0.000) 0.160 1.681 ( -2.402 -2.402 0.000) 3.396 1.681 ( 2.402 2.402 0.000) 3.396 1.701 ( -1.119 -1.119 0.000) 1.583 1.701 ( 1.119 1.119 0.000) 1.583 1.720 ( -1.516 -1.516 0.000) 2.144 1.720 ( 1.516 1.516 0.000) 2.144 1.778 ( -0.997 -0.997 0.000) 1.410 1.778 ( 0.997 0.997 0.000) 1.410 1.794 ( -2.157 -2.157 0.000) 3.050 1.794 ( 2.157 2.157 0.000) 3.050 1.847 ( -0.559 -0.559 0.000) 0.790 1.847 ( 0.559 0.559 0.000) 0.790 1.860 ( -0.873 -0.873 0.000) 1.234 1.860 ( 0.873 0.873 0.000) 1.234 1.900 ( -2.354 -2.354 0.000) 3.330 1.900 ( 2.354 2.354 0.000) 3.330 1.913 ( -2.106 -2.106 0.000) 2.979 1.913 ( 2.106 2.106 0.000) 2.979 2.028 ( -0.891 -0.891 0.000) 1.260 2.028 ( 0.891 0.891 0.000) 1.260 2.045 ( -2.485 -2.485 0.000) 3.514 2.045 ( 2.485 2.485 0.000) 3.514 2.105 ( -1.708 -1.708 0.000) 2.416 2.105 ( 1.708 1.708 0.000) 2.416 2.228 ( -0.590 -0.590 0.000) 0.835 2.228 ( 0.590 0.590 0.000) 0.835 2.271 ( -0.331 -0.331 0.000) 0.468 2.271 ( 0.331 0.331 0.000) 0.468 2.299 ( -1.102 -1.102 0.000) 1.558 2.299 ( 1.102 1.102 0.000) 1.558 2.306 ( -1.385 -1.385 0.000) 1.959 2.306 ( 1.385 1.385 0.000) 1.959 2.349 ( -0.077 -0.077 0.000) 0.109 2.349 ( 0.077 0.077 0.000) 0.109 2.395 ( -1.173 -1.173 0.000) 1.659 2.395 ( 1.173 1.173 0.000) 1.659 2.425 ( -1.067 -1.067 0.000) 1.509 2.425 ( 1.067 1.067 0.000) 1.509 2.473 ( -0.075 -0.075 0.000) 0.106 2.473 ( 0.075 0.075 0.000) 0.106 2.729 ( -0.383 -0.383 0.000) 0.541 2.729 ( 0.383 0.383 0.000) 0.541 2.753 ( -0.136 -0.136 0.000) 0.193 2.753 ( 0.136 0.136 0.000) 0.193 2.783 ( -0.007 -0.007 0.000) 0.010 2.783 ( 0.007 0.007 0.000) 0.010 2.841 ( -0.209 -0.209 0.000) 0.295 2.841 ( 0.209 0.209 0.000) 0.295 2.888 ( -0.420 -0.420 0.000) 0.594 2.888 ( 0.420 0.420 0.000) 0.594 2.910 ( -0.046 -0.046 0.000) 0.065 2.910 ( 0.046 0.046 0.000) 0.065 3.877 ( -0.384 -0.384 0.000) 0.543 3.877 ( 0.384 0.384 0.000) 0.543 3.887 ( -0.107 -0.107 0.000) 0.151 3.887 ( 0.107 0.107 0.000) 0.151 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.315 ( 0.000 -9.264 -6.005) 11.040 0.359 ( 0.000 -8.123 -9.190) 12.266 0.466 ( 0.000 -11.143 -9.903) 14.908 0.646 ( 0.000 -0.768 2.268) 2.395 0.649 ( 0.000 -1.204 2.126) 2.443 0.723 ( -0.000 1.198 -0.587) 1.334 0.729 ( -0.000 1.249 0.862) 1.517 0.743 ( -0.000 3.962 2.453) 4.661 0.757 ( -0.000 -0.286 4.284) 4.294 0.775 ( -0.000 1.957 2.499) 3.174 0.796 ( -0.000 0.607 3.445) 3.498 0.819 ( -0.000 -0.426 2.915) 2.946 0.838 ( -0.000 2.189 -0.478) 2.240 0.840 ( -0.000 1.807 0.282) 1.829 0.871 ( -0.000 1.086 -2.412) 2.645 0.899 ( -0.000 0.527 0.192) 0.561 0.911 ( 0.000 0.083 -0.347) 0.357 0.913 ( 0.000 0.107 -0.814) 0.821 0.915 ( -0.000 1.608 1.623) 2.284 0.928 ( -0.000 0.217 0.770) 0.800 0.952 ( -0.000 0.943 -1.991) 2.203 0.961 ( 0.000 0.056 -0.682) 0.684 0.978 ( 0.000 -0.057 -0.983) 0.985 1.014 ( 0.000 -0.137 -0.019) 0.138 1.174 ( 0.000 -0.386 0.639) 0.746 1.208 ( 0.000 -1.151 0.082) 1.154 1.241 ( -0.000 1.438 -1.103) 1.812 1.248 ( -0.000 0.414 -0.759) 0.865 1.292 ( 0.000 -1.452 -0.009) 1.452 1.315 ( 0.000 -2.249 0.953) 2.442 1.331 ( 0.000 -0.191 0.441) 0.481 1.338 ( 0.000 0.266 -1.518) 1.541 1.347 ( 0.000 -0.878 -2.232) 2.399 1.382 ( 0.000 -1.716 -2.388) 2.941 1.486 ( -0.000 2.419 1.247) 2.721 1.500 ( 0.000 -2.029 -1.035) 2.278 1.517 ( -0.000 1.325 1.412) 1.936 1.544 ( -0.000 0.766 0.976) 1.241 1.557 ( -0.000 1.501 2.257) 2.711 1.586 ( -0.000 0.034 1.325) 1.325 1.627 ( 0.000 -0.452 -0.810) 0.927 1.643 ( 0.000 -0.866 0.844) 1.209 1.674 ( -0.000 0.500 0.196) 0.537 1.679 ( -0.000 0.432 1.693) 1.747 1.736 ( -0.000 0.552 1.520) 1.617 1.737 ( 0.000 -1.816 1.610) 2.427 1.771 ( 0.000 -2.971 0.871) 3.095 1.787 ( -0.000 0.364 1.766) 1.803 1.794 ( -0.000 1.369 0.376) 1.420 1.816 ( 0.000 -0.768 0.394) 0.863 1.844 ( 0.000 -0.759 1.730) 1.890 1.847 ( 0.000 0.194 -4.489) 4.493 1.861 ( -0.000 1.216 -2.640) 2.907 1.883 ( -0.000 4.821 -1.306) 4.995 1.935 ( -0.000 -0.600 3.864) 3.910 1.951 ( -0.000 0.356 -0.334) 0.488 1.999 ( -0.000 3.946 0.171) 3.949 2.002 ( 0.000 -1.807 1.713) 2.490 2.020 ( -0.000 2.964 0.299) 2.979 2.028 ( -0.000 1.428 -1.572) 2.124 2.076 ( 0.000 0.290 -2.675) 2.691 2.084 ( 0.000 -0.836 1.056) 1.346 2.097 ( 0.000 -0.606 -2.028) 2.117 2.144 ( 0.000 -1.449 -0.410) 1.506 2.149 ( 0.000 -1.752 -1.757) 2.482 2.188 ( -0.000 0.943 -2.120) 2.321 2.218 ( 0.000 -1.542 -0.729) 1.706 2.229 ( 0.000 -1.906 -0.043) 1.907 2.247 ( -0.000 0.653 1.948) 2.054 2.282 ( -0.000 0.135 -0.263) 0.296 2.308 ( 0.000 -0.915 0.497) 1.041 2.308 ( -0.000 0.417 0.085) 0.426 2.352 ( 0.000 -0.036 -0.048) 0.059 2.372 ( 0.000 0.109 -2.775) 2.777 2.376 ( 0.000 -0.548 -1.736) 1.821 2.389 ( 0.000 -0.877 -2.826) 2.959 2.427 ( 0.000 -0.386 0.705) 0.804 2.439 ( -0.000 0.347 1.740) 1.774 2.450 ( 0.000 -0.599 0.513) 0.789 2.462 ( -0.000 0.459 1.092) 1.184 2.710 ( 0.000 -0.697 -0.811) 1.069 2.731 ( 0.000 -1.222 0.331) 1.266 2.744 ( 0.000 -0.733 -1.532) 1.699 2.749 ( 0.000 -1.014 -1.005) 1.428 2.826 ( -0.000 1.050 1.845) 2.123 2.826 ( -0.000 0.975 1.808) 2.054 2.842 ( -0.000 0.903 0.871) 1.254 2.844 ( -0.000 0.655 0.966) 1.167 2.880 ( 0.000 -0.017 -0.567) 0.567 2.883 ( 0.000 -0.349 0.450) 0.569 2.887 ( 0.000 -0.568 0.218) 0.608 2.917 ( -0.000 0.820 -0.268) 0.863 3.859 ( 0.000 -0.416 -0.484) 0.639 3.885 ( -0.000 0.151 0.336) 0.368 3.889 ( -0.000 0.016 0.046) 0.048 3.896 ( -0.000 0.266 0.363) 0.450 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.505 ( 6.804 0.634 -3.176) 7.535 0.525 ( 2.447 1.433 -5.025) 5.770 0.578 ( -3.826 -0.376 0.591) 3.890 0.629 ( -2.271 -2.121 -0.270) 3.119 0.657 ( 1.720 -1.758 1.133) 2.707 0.675 ( -1.833 1.205 1.565) 2.694 0.700 ( -2.461 0.312 0.520) 2.535 0.721 ( 0.658 0.875 -0.208) 1.115 0.745 ( -0.048 2.318 0.145) 2.323 0.771 ( -0.167 -0.202 1.489) 1.512 0.779 ( 1.060 -0.155 -0.090) 1.075 0.800 ( -1.052 -0.279 1.594) 1.931 0.831 ( -0.812 0.612 -0.447) 1.111 0.849 ( -1.023 1.796 -1.084) 2.334 0.861 ( -1.511 0.240 -1.141) 1.908 0.879 ( 2.163 1.534 0.292) 2.668 0.901 ( 0.169 -0.393 1.479) 1.540 0.907 ( -0.255 -0.588 -1.219) 1.377 0.914 ( -0.092 0.260 2.380) 2.396 0.938 ( 0.562 0.176 0.748) 0.952 0.941 ( 1.015 -0.605 0.199) 1.198 0.956 ( 1.100 0.073 1.110) 1.565 0.970 ( -0.189 -0.557 -2.131) 2.211 0.993 ( -0.657 0.095 0.001) 0.664 1.185 ( 1.609 -1.388 -0.778) 2.263 1.197 ( -0.492 -0.843 -0.351) 1.037 1.222 ( -0.906 -0.952 0.840) 1.559 1.238 ( 0.008 -0.319 1.570) 1.602 1.312 ( -1.104 1.031 0.050) 1.511 1.322 ( 2.221 -0.063 1.195) 2.523 1.359 ( -0.309 1.119 -1.803) 2.144 1.382 ( 2.080 0.081 -0.724) 2.204 1.398 ( -0.150 2.089 0.551) 2.166 1.416 ( -0.764 0.027 1.709) 1.872 1.433 ( 1.498 -1.392 -2.364) 3.126 1.459 ( -0.924 1.422 -0.732) 1.847 1.494 ( -0.543 -1.378 -1.460) 2.080 1.532 ( -0.783 -0.076 2.463) 2.586 1.551 ( -2.008 -1.159 -1.026) 2.535 1.580 ( 0.629 -1.312 1.151) 1.855 1.639 ( -0.094 1.673 -0.458) 1.737 1.652 ( 1.454 0.873 0.119) 1.700 1.679 ( -1.038 1.794 1.069) 2.332 1.688 ( -0.227 1.456 1.306) 1.969 1.742 ( 0.260 0.158 -0.802) 0.858 1.751 ( -1.120 -0.467 0.167) 1.225 1.773 ( -1.013 0.875 -1.323) 1.882 1.791 ( 1.767 -0.036 -0.351) 1.802 1.806 ( -2.991 -1.269 -0.961) 3.389 1.821 ( -0.493 0.236 0.227) 0.591 1.827 ( -0.079 0.756 -0.213) 0.789 1.832 ( -0.810 0.196 0.063) 0.836 1.854 ( -0.593 -0.667 -0.990) 1.333 1.865 ( 1.133 0.492 -0.048) 1.236 1.904 ( 1.151 -2.617 -0.989) 3.025 1.928 ( -0.848 -1.984 -1.637) 2.708 1.936 ( -0.383 -1.502 1.533) 2.180 1.981 ( 1.921 -0.605 1.607) 2.576 1.984 ( -1.397 -0.999 0.270) 1.738 1.993 ( 1.433 0.131 2.289) 2.704 2.068 ( 0.076 -0.360 -0.463) 0.591 2.072 ( -0.067 -1.219 0.122) 1.227 2.119 ( 0.930 -1.430 -1.532) 2.293 2.140 ( 0.004 -1.485 0.592) 1.599 2.212 ( 0.123 1.878 0.400) 1.924 2.225 ( 0.814 1.188 -0.355) 1.483 2.235 ( 0.218 1.659 0.362) 1.712 2.260 ( 0.716 1.283 0.619) 1.594 2.267 ( 0.953 0.490 -1.259) 1.653 2.283 ( -0.066 1.041 -0.320) 1.091 2.309 ( -0.250 -0.700 0.967) 1.219 2.323 ( -0.695 0.727 1.411) 1.733 2.351 ( -0.285 -0.452 -1.103) 1.225 2.365 ( 0.893 0.121 0.310) 0.953 2.377 ( 0.339 -1.245 -1.485) 1.967 2.406 ( -0.589 -1.012 1.501) 1.904 2.418 ( 0.685 0.877 -1.990) 2.280 2.448 ( -0.407 0.322 -1.003) 1.129 2.453 ( -0.501 0.539 1.465) 1.639 2.464 ( 0.149 0.362 0.676) 0.781 2.730 ( 0.611 0.323 -0.757) 1.025 2.741 ( 0.339 -0.223 0.196) 0.451 2.756 ( 0.402 -0.504 -0.696) 0.949 2.767 ( 0.013 -0.487 -0.735) 0.882 2.787 ( -0.734 -0.259 -0.003) 0.778 2.808 ( -0.038 -0.521 -0.427) 0.675 2.815 ( -0.768 0.205 1.226) 1.461 2.831 ( -0.539 -0.159 1.095) 1.230 2.885 ( 0.041 0.502 -0.618) 0.797 2.896 ( -0.018 0.422 0.387) 0.573 2.900 ( -0.575 0.083 -0.201) 0.615 2.906 ( 0.470 0.172 0.227) 0.549 3.873 ( 0.482 0.262 -0.372) 0.663 3.880 ( -0.335 0.166 0.437) 0.575 3.886 ( -0.122 -0.217 -0.176) 0.305 3.889 ( -0.115 -0.139 0.112) 0.212 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 33 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.495 ( -0.000 4.856 7.061) 8.569 0.514 ( -0.000 8.575 2.617) 8.965 0.596 ( 0.000 -6.276 -2.474) 6.746 0.643 ( 0.000 -3.677 -4.652) 5.930 0.653 ( 0.000 -3.296 0.730) 3.376 0.665 ( -0.000 4.438 0.013) 4.438 0.693 ( 0.000 -0.960 1.691) 1.945 0.724 ( -0.000 1.591 1.713) 2.338 0.724 ( 0.000 -1.685 0.037) 1.685 0.768 ( 0.000 -0.870 -3.817) 3.915 0.789 ( -0.000 -0.234 1.831) 1.846 0.794 ( 0.000 -2.687 -1.261) 2.969 0.803 ( 0.000 -0.354 -0.563) 0.665 0.825 ( -0.000 0.607 -1.680) 1.787 0.856 ( -0.000 0.766 1.244) 1.461 0.879 ( -0.000 0.231 0.386) 0.449 0.900 ( -0.000 -0.053 0.193) 0.200 0.911 ( 0.000 0.368 -2.654) 2.680 0.921 ( 0.000 -0.878 0.463) 0.993 0.943 ( -0.000 2.084 1.328) 2.472 0.946 ( -0.000 1.726 -0.879) 1.937 0.959 ( 0.000 -0.472 -0.471) 0.667 0.980 ( -0.000 0.083 1.388) 1.390 0.988 ( 0.000 -1.538 0.185) 1.549 1.201 ( 0.000 -2.168 0.298) 2.189 1.216 ( -0.000 2.492 -1.421) 2.869 1.234 ( -0.000 1.193 0.121) 1.199 1.237 ( 0.000 -0.540 -0.174) 0.567 1.308 ( -0.000 0.466 1.163) 1.253 1.326 ( 0.000 -0.989 0.624) 1.169 1.348 ( 0.000 0.643 -2.972) 3.040 1.349 ( -0.000 1.417 -1.721) 2.229 1.388 ( -0.000 2.714 2.303) 3.559 1.416 ( 0.000 -2.917 -0.381) 2.942 1.429 ( -0.000 2.255 1.422) 2.666 1.452 ( 0.000 -2.723 -1.282) 3.009 1.526 ( 0.000 -0.790 -0.381) 0.877 1.529 ( 0.000 -0.737 -1.062) 1.293 1.571 ( -0.000 1.201 2.720) 2.973 1.604 ( 0.000 -0.246 -0.376) 0.450 1.615 ( -0.000 -0.056 2.577) 2.578 1.625 ( 0.000 -0.999 1.932) 2.175 1.665 ( 0.000 -0.790 -3.010) 3.112 1.667 ( 0.000 -0.276 -2.651) 2.665 1.724 ( -0.000 -0.037 1.205) 1.206 1.748 ( 0.000 -1.538 -1.330) 2.034 1.785 ( -0.000 1.974 0.927) 2.180 1.789 ( 0.000 -2.851 0.576) 2.908 1.807 ( -0.000 1.310 1.267) 1.822 1.827 ( 0.000 -2.548 -0.775) 2.664 1.827 ( 0.000 -2.500 1.955) 3.174 1.842 ( 0.000 -0.787 1.533) 1.723 1.855 ( -0.000 2.734 -1.198) 2.985 1.864 ( -0.000 3.716 -1.059) 3.864 1.930 ( 0.000 -2.976 1.329) 3.259 1.955 ( -0.000 -0.312 1.561) 1.592 1.956 ( -0.000 0.657 -1.445) 1.588 1.963 ( -0.000 0.514 0.866) 1.007 1.975 ( 0.000 -0.239 -1.254) 1.277 2.016 ( 0.000 0.788 -4.128) 4.203 2.086 ( -0.000 -0.014 1.474) 1.474 2.104 ( -0.000 0.171 -0.280) 0.328 2.147 ( 0.000 -1.046 0.428) 1.130 2.156 ( -0.000 3.765 0.282) 3.775 2.190 ( -0.000 1.545 0.220) 1.561 2.200 ( -0.000 1.536 1.388) 2.070 2.218 ( 0.000 -1.319 -2.485) 2.813 2.227 ( 0.000 -1.096 -1.731) 2.049 2.259 ( 0.000 -0.817 0.693) 1.071 2.265 ( -0.000 0.900 -0.039) 0.901 2.307 ( -0.000 0.293 1.104) 1.143 2.325 ( -0.000 0.482 -1.242) 1.332 2.362 ( 0.000 -0.700 0.868) 1.115 2.364 ( -0.000 1.036 0.302) 1.079 2.395 ( -0.000 1.060 1.502) 1.838 2.399 ( -0.000 0.019 2.151) 2.151 2.429 ( 0.000 -0.609 -1.485) 1.605 2.432 ( 0.000 0.026 -0.396) 0.397 2.454 ( 0.000 -0.364 -0.264) 0.450 2.456 ( 0.000 -0.273 -0.955) 0.993 2.726 ( -0.000 0.612 0.903) 1.091 2.745 ( 0.000 0.016 -0.926) 0.926 2.764 ( -0.000 1.094 1.808) 2.113 2.769 ( -0.000 0.830 1.272) 1.519 2.800 ( 0.000 -1.231 -1.469) 1.916 2.802 ( 0.000 -1.292 -1.602) 2.058 2.823 ( 0.000 -0.622 0.321) 0.700 2.830 ( 0.000 -0.512 -0.296) 0.591 2.879 ( -0.000 -0.162 0.524) 0.549 2.889 ( -0.000 0.284 -0.545) 0.614 2.899 ( 0.000 -0.863 0.330) 0.924 2.903 ( -0.000 0.732 -0.106) 0.740 3.870 ( -0.000 0.550 0.419) 0.691 3.878 ( 0.000 -0.422 -0.512) 0.663 3.888 ( -0.000 -0.043 0.140) 0.147 3.890 ( 0.000 -0.212 -0.052) 0.218 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 32 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.451 ( 3.442 -3.442 -5.334) 7.221 0.556 ( 4.362 -4.362 -2.259) 6.570 0.592 ( 0.252 -0.252 1.882) 1.915 0.623 ( -0.910 0.910 3.139) 3.393 0.633 ( -0.548 0.548 1.305) 1.518 0.679 ( -0.920 0.920 2.100) 2.470 0.690 ( -1.515 1.515 1.562) 2.651 0.716 ( 0.333 -0.333 3.205) 3.240 0.729 ( 1.976 -1.976 -4.432) 5.239 0.749 ( 0.924 -0.924 -3.558) 3.791 0.771 ( 0.735 -0.735 -0.304) 1.083 0.802 ( 1.620 -1.621 0.438) 2.334 0.802 ( 0.400 -0.399 3.116) 3.167 0.828 ( -1.320 1.320 0.459) 1.922 0.853 ( -0.092 0.092 -0.036) 0.135 0.897 ( -0.481 0.481 -0.060) 0.682 0.901 ( -0.472 0.472 -1.004) 1.206 0.917 ( -0.067 0.067 0.088) 0.129 0.921 ( -0.776 0.776 -0.347) 1.151 0.935 ( -0.280 0.280 -1.386) 1.442 0.963 ( 0.207 -0.207 0.946) 0.991 0.972 ( 0.081 -0.081 -0.120) 0.166 0.991 ( -0.855 0.855 -0.876) 1.493 1.006 ( -1.303 1.303 -0.162) 1.849 1.168 ( 0.225 -0.225 -1.011) 1.059 1.207 ( 0.545 -0.545 0.777) 1.095 1.241 ( -1.860 1.860 -0.089) 2.632 1.247 ( -0.831 0.831 -0.061) 1.177 1.299 ( 0.183 -0.183 -2.275) 2.290 1.347 ( 2.165 -2.165 -1.689) 3.497 1.352 ( -0.478 0.478 -0.733) 0.997 1.362 ( 0.859 -0.859 -0.377) 1.272 1.394 ( 0.666 -0.666 -1.029) 1.395 1.437 ( 0.124 -0.124 -2.606) 2.612 1.443 ( -1.582 1.582 1.444) 2.663 1.470 ( -0.799 0.799 3.767) 3.933 1.492 ( -0.528 0.528 2.957) 3.050 1.513 ( 0.166 -0.166 1.118) 1.142 1.529 ( 0.110 -0.110 -2.140) 2.145 1.582 ( 0.189 -0.189 0.472) 0.543 1.637 ( -0.428 0.428 -0.510) 0.791 1.653 ( 2.234 -2.234 0.769) 3.252 1.667 ( 0.856 -0.856 -1.753) 2.130 1.698 ( 0.275 -0.275 0.262) 0.469 1.718 ( 0.938 -0.938 0.353) 1.373 1.721 ( 0.907 -0.907 0.528) 1.388 1.744 ( 0.377 -0.377 2.727) 2.779 1.766 ( 1.637 -1.637 1.710) 2.878 1.791 ( 0.782 -0.782 1.284) 1.695 1.794 ( 0.583 -0.583 3.655) 3.746 1.811 ( 0.416 -0.416 0.432) 0.730 1.838 ( -1.941 1.941 0.761) 2.848 1.872 ( -0.379 0.379 -1.346) 1.449 1.916 ( -2.904 2.904 0.627) 4.154 1.925 ( -0.966 0.966 -1.368) 1.933 1.948 ( -1.968 1.968 0.625) 2.853 1.954 ( -0.286 0.286 -1.736) 1.783 2.001 ( -2.614 2.614 -1.480) 3.982 2.010 ( -0.099 0.099 -2.042) 2.047 2.035 ( 0.133 -0.133 1.052) 1.069 2.051 ( -0.163 0.163 0.873) 0.903 2.076 ( -0.599 0.599 0.647) 1.066 2.103 ( -0.858 0.858 -3.277) 3.495 2.182 ( -1.207 1.207 1.770) 2.460 2.201 ( 0.525 -0.525 -2.532) 2.639 2.228 ( 0.557 -0.557 0.741) 1.082 2.238 ( 2.512 -2.512 -0.950) 3.678 2.244 ( 1.516 -1.516 -0.465) 2.193 2.249 ( -0.331 0.331 -2.493) 2.537 2.287 ( 0.442 -0.442 -0.062) 0.628 2.309 ( -0.471 0.471 -0.514) 0.841 2.317 ( -0.173 0.173 -0.744) 0.784 2.352 ( -0.204 0.204 0.017) 0.288 2.383 ( -0.573 0.573 2.381) 2.515 2.395 ( -0.555 0.555 -1.779) 1.944 2.396 ( -0.472 0.472 0.241) 0.709 2.430 ( 0.249 -0.249 0.743) 0.822 2.432 ( 0.224 -0.224 1.484) 1.518 2.441 ( 0.682 -0.682 -1.328) 1.642 2.459 ( 0.619 -0.619 -0.397) 0.961 2.735 ( 0.383 -0.383 -0.408) 0.678 2.736 ( 0.266 -0.266 -0.712) 0.805 2.762 ( 0.007 -0.007 -1.042) 1.042 2.774 ( -0.826 0.826 1.161) 1.647 2.786 ( -0.371 0.371 2.055) 2.121 2.796 ( 0.866 -0.866 -1.234) 1.738 2.816 ( -0.454 0.454 1.646) 1.767 2.825 ( 0.045 -0.045 0.048) 0.079 2.883 ( 0.463 -0.463 0.499) 0.823 2.884 ( 0.389 -0.389 0.459) 0.717 2.895 ( 0.137 -0.137 -0.802) 0.825 2.914 ( -0.237 0.237 0.126) 0.359 3.873 ( 0.303 -0.303 -0.613) 0.748 3.879 ( -0.044 0.044 0.193) 0.203 3.885 ( -0.166 0.166 0.453) 0.510 3.887 ( 0.044 -0.044 0.145) 0.158 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.360 ( -5.454 5.454 4.708) 9.037 0.430 ( -4.711 4.711 10.767) 12.661 0.585 ( -8.010 8.010 2.906) 11.695 0.646 ( -0.552 0.552 -1.432) 1.631 0.658 ( 0.020 -0.020 -1.041) 1.041 0.694 ( 2.640 -2.640 -1.947) 4.211 0.715 ( 3.287 -3.287 -1.613) 4.921 0.729 ( 0.067 -0.067 0.100) 0.138 0.752 ( 0.217 -0.217 0.836) 0.890 0.757 ( 0.436 -0.436 -4.244) 4.288 0.784 ( 0.184 -0.184 -2.946) 2.958 0.815 ( 0.536 -0.536 -0.948) 1.213 0.832 ( -0.490 0.490 -2.252) 2.356 0.850 ( 0.182 -0.182 0.010) 0.257 0.866 ( 0.677 -0.677 1.889) 2.117 0.891 ( 0.766 -0.766 1.742) 2.051 0.898 ( 0.142 -0.142 -3.100) 3.107 0.908 ( 0.175 -0.175 0.434) 0.500 0.921 ( 0.397 -0.397 -0.751) 0.937 0.933 ( -0.110 0.110 -1.875) 1.882 0.949 ( 0.384 -0.384 1.902) 1.978 0.950 ( 0.510 -0.510 2.228) 2.342 0.980 ( 0.209 -0.209 0.615) 0.683 1.010 ( 0.417 -0.417 -0.133) 0.605 1.154 ( 0.250 -0.250 0.640) 0.731 1.209 ( 0.205 -0.205 0.550) 0.621 1.242 ( -0.093 0.093 0.353) 0.377 1.250 ( 0.788 -0.788 -1.400) 1.790 1.283 ( -0.476 0.476 -0.715) 0.982 1.320 ( -0.667 0.667 0.640) 1.140 1.347 ( -0.473 0.473 0.716) 0.980 1.360 ( -1.957 1.957 1.040) 2.956 1.370 ( -1.094 1.094 2.360) 2.822 1.400 ( -1.385 1.385 0.496) 2.021 1.458 ( 2.145 -2.145 -0.708) 3.115 1.479 ( 1.582 -1.582 1.454) 2.668 1.500 ( 0.011 -0.011 0.919) 0.920 1.539 ( -0.656 0.656 -3.171) 3.304 1.561 ( 0.717 -0.717 -3.048) 3.212 1.605 ( 0.107 -0.107 1.016) 1.027 1.626 ( 1.022 -1.022 -1.066) 1.796 1.652 ( 0.554 -0.554 -0.546) 0.955 1.671 ( -1.605 1.605 -1.150) 2.544 1.699 ( -0.777 0.777 0.161) 1.111 1.725 ( -0.989 0.989 -0.385) 1.450 1.739 ( -0.008 0.008 -1.293) 1.293 1.780 ( -0.108 0.108 0.137) 0.205 1.781 ( 0.313 -0.313 -1.050) 1.140 1.800 ( -0.987 0.987 -0.550) 1.501 1.828 ( -0.749 0.749 -0.362) 1.119 1.848 ( -0.360 0.360 -1.967) 2.031 1.850 ( 1.873 -1.873 -0.151) 2.654 1.851 ( -0.248 0.248 1.434) 1.476 1.871 ( -0.096 0.096 3.190) 3.193 1.916 ( 4.094 -4.094 0.488) 5.810 1.961 ( -0.192 0.192 -1.309) 1.337 1.962 ( -0.601 0.601 -1.879) 2.062 1.981 ( 2.334 -2.334 -0.609) 3.357 1.989 ( 1.875 -1.875 0.400) 2.682 2.045 ( 2.041 -2.041 -0.018) 2.886 2.055 ( -0.192 0.192 2.024) 2.042 2.062 ( 0.632 -0.632 0.037) 0.895 2.090 ( -0.255 0.255 -1.023) 1.084 2.148 ( -1.019 1.019 3.097) 3.415 2.200 ( -1.678 1.678 1.422) 2.766 2.202 ( -0.589 0.589 1.887) 2.063 2.224 ( -1.763 1.763 1.213) 2.772 2.233 ( -0.117 0.117 -2.648) 2.653 2.238 ( -0.872 0.872 1.586) 2.009 2.289 ( 0.630 -0.630 -0.302) 0.941 2.314 ( -0.151 0.151 -0.685) 0.718 2.319 ( -0.980 0.980 -1.231) 1.853 2.353 ( -0.158 0.158 -0.121) 0.255 2.363 ( 0.229 -0.229 1.645) 1.676 2.377 ( 0.506 -0.506 1.942) 2.069 2.398 ( -0.882 0.882 2.501) 2.795 2.423 ( 0.562 -0.562 -1.808) 1.975 2.446 ( -0.511 0.511 -0.097) 0.729 2.454 ( -0.372 0.372 -0.475) 0.708 2.462 ( 0.063 -0.063 -1.062) 1.065 2.717 ( -0.500 0.500 1.038) 1.256 2.729 ( -0.315 0.315 0.146) 0.469 2.749 ( -0.422 0.422 0.672) 0.899 2.759 ( -0.567 0.567 1.588) 1.779 2.810 ( 1.113 -1.113 -1.546) 2.207 2.820 ( 0.780 -0.780 -1.130) 1.579 2.829 ( 0.631 -0.631 -0.923) 1.284 2.843 ( 0.446 -0.446 -1.306) 1.450 2.882 ( -0.170 0.170 0.641) 0.685 2.892 ( -0.538 0.538 -0.408) 0.863 2.892 ( -0.417 0.417 -0.499) 0.772 2.912 ( 0.396 -0.396 0.273) 0.623 3.863 ( -0.373 0.373 0.456) 0.697 3.884 ( 0.054 -0.054 -0.296) 0.306 3.888 ( 0.076 -0.076 0.011) 0.107 3.893 ( 0.233 -0.233 -0.354) 0.483 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 20 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.220 ( -0.000 9.981 0.000) 9.981 0.256 ( -0.000 11.723 0.000) 11.723 0.393 ( -0.000 17.007 0.000) 17.007 0.667 ( 0.000 -0.320 -0.000) 0.320 0.686 ( -0.000 0.825 0.000) 0.825 0.693 ( 0.000 -1.087 -0.000) 1.087 0.725 ( -0.000 1.247 0.000) 1.247 0.778 ( 0.000 -3.626 -0.000) 3.626 0.803 ( 0.000 -4.028 -0.000) 4.028 0.815 ( 0.000 -0.788 -0.000) 0.788 0.821 ( 0.000 -3.276 -0.000) 3.276 0.829 ( 0.000 -0.097 -0.000) 0.097 0.844 ( 0.000 -1.795 -0.000) 1.795 0.853 ( 0.000 -1.661 -0.000) 1.661 0.870 ( 0.000 -0.827 -0.000) 0.827 0.888 ( -0.000 0.543 0.000) 0.543 0.896 ( 0.000 -2.182 -0.000) 2.182 0.898 ( 0.000 -0.282 -0.000) 0.282 0.912 ( -0.000 0.487 0.000) 0.487 0.923 ( 0.000 -0.465 -0.000) 0.465 0.935 ( -0.000 0.526 0.000) 0.526 0.970 ( 0.000 -0.144 -0.000) 0.144 0.980 ( 0.000 -1.150 -0.000) 1.150 1.010 ( -0.000 0.732 0.000) 0.732 1.183 ( -0.000 0.835 0.000) 0.835 1.211 ( -0.000 0.677 0.000) 0.677 1.229 ( 0.000 -1.232 -0.000) 1.232 1.241 ( 0.000 -0.421 -0.000) 0.421 1.294 ( -0.000 0.249 0.000) 0.249 1.315 ( -0.000 0.959 0.000) 0.959 1.322 ( -0.000 0.722 0.000) 0.722 1.339 ( -0.000 2.159 0.000) 2.159 1.340 ( -0.000 3.004 0.000) 3.004 1.357 ( -0.000 2.925 0.000) 2.925 1.486 ( -0.000 0.383 0.000) 0.383 1.501 ( 0.000 -1.098 -0.000) 1.098 1.523 ( 0.000 -2.905 -0.000) 2.905 1.552 ( 0.000 -0.769 -0.000) 0.769 1.602 ( 0.000 -1.237 -0.000) 1.237 1.605 ( 0.000 -1.025 -0.000) 1.025 1.608 ( 0.000 -0.479 -0.000) 0.479 1.643 ( 0.000 -3.849 -0.000) 3.849 1.677 ( -0.000 1.484 0.000) 1.484 1.712 ( -0.000 0.809 0.000) 0.809 1.755 ( 0.000 -2.951 -0.000) 2.951 1.756 ( -0.000 2.058 0.000) 2.058 1.765 ( 0.000 -1.375 -0.000) 1.375 1.783 ( -0.000 3.953 0.000) 3.953 1.795 ( -0.000 1.051 0.000) 1.051 1.816 ( -0.000 1.249 0.000) 1.249 1.825 ( 0.000 -3.282 -0.000) 3.282 1.838 ( 0.000 -0.925 -0.000) 0.925 1.851 ( 0.000 -3.795 -0.000) 3.795 1.866 ( -0.000 0.700 0.000) 0.700 1.955 ( 0.000 -1.566 -0.000) 1.566 1.980 ( 0.000 -0.080 -0.000) 0.080 1.998 ( -0.000 2.142 0.000) 2.142 2.008 ( 0.000 -4.650 -0.000) 4.650 2.019 ( 0.000 -0.423 -0.000) 0.423 2.025 ( -0.000 1.512 0.000) 1.512 2.051 ( 0.000 -1.477 -0.000) 1.477 2.074 ( -0.000 3.264 0.000) 3.264 2.101 ( -0.000 1.456 0.000) 1.456 2.125 ( -0.000 0.947 0.000) 0.947 2.138 ( -0.000 0.402 0.000) 0.402 2.167 ( 0.000 -0.833 -0.000) 0.833 2.203 ( -0.000 0.947 0.000) 0.947 2.226 ( -0.000 2.982 0.000) 2.982 2.282 ( 0.000 -0.096 -0.000) 0.096 2.282 ( 0.000 -1.133 -0.000) 1.133 2.309 ( 0.000 -1.836 -0.000) 1.836 2.323 ( -0.000 0.785 0.000) 0.785 2.335 ( -0.000 1.206 0.000) 1.206 2.351 ( -0.000 1.169 0.000) 1.169 2.353 ( -0.000 0.941 0.000) 0.941 2.358 ( -0.000 0.305 0.000) 0.305 2.430 ( -0.000 0.070 0.000) 0.070 2.455 ( 0.000 -0.235 -0.000) 0.235 2.459 ( -0.000 1.037 0.000) 1.037 2.472 ( 0.000 -0.590 -0.000) 0.590 2.706 ( -0.000 0.554 0.000) 0.554 2.727 ( -0.000 0.516 0.000) 0.516 2.734 ( -0.000 0.985 0.000) 0.985 2.744 ( -0.000 0.420 0.000) 0.420 2.844 ( 0.000 -1.088 -0.000) 1.088 2.845 ( 0.000 -0.910 -0.000) 0.910 2.851 ( 0.000 -0.862 -0.000) 0.862 2.852 ( 0.000 -1.238 -0.000) 1.238 2.876 ( 0.000 -0.025 -0.000) 0.025 2.888 ( -0.000 0.188 0.000) 0.188 2.888 ( -0.000 0.707 0.000) 0.707 2.911 ( 0.000 -0.569 -0.000) 0.569 3.854 ( -0.000 0.431 0.000) 0.431 3.888 ( 0.000 -0.129 -0.000) 0.129 3.889 ( 0.000 -0.070 -0.000) 0.070 3.899 ( 0.000 -0.337 -0.000) 0.337 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 32 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.451 ( 2.779 6.421 0.000) 6.996 0.478 ( 0.128 7.535 0.000) 7.536 0.608 ( -1.737 -4.352 -0.000) 4.686 0.611 ( -1.291 -4.158 -0.000) 4.354 0.660 ( -1.069 -1.237 -0.000) 1.635 0.673 ( -3.199 -0.524 0.000) 3.242 0.713 ( 2.718 0.681 -0.000) 2.802 0.725 ( 2.826 1.535 -0.000) 3.216 0.739 ( 1.572 -2.074 -0.000) 2.602 0.786 ( -0.232 -1.710 -0.000) 1.726 0.792 ( 1.740 -1.359 -0.000) 2.207 0.824 ( -0.674 0.120 0.000) 0.685 0.827 ( -2.338 0.299 0.000) 2.358 0.831 ( 1.551 -3.100 -0.000) 3.467 0.878 ( 0.885 -0.257 -0.000) 0.921 0.882 ( -0.067 0.289 0.000) 0.296 0.886 ( 1.042 1.718 0.000) 2.009 0.905 ( -1.303 0.744 0.000) 1.501 0.927 ( 0.789 3.728 0.000) 3.810 0.945 ( 1.226 0.192 -0.000) 1.241 0.950 ( 0.588 -0.251 -0.000) 0.639 0.957 ( 0.123 -0.920 -0.000) 0.928 0.975 ( 0.171 1.580 0.000) 1.590 0.982 ( -0.742 -1.413 -0.000) 1.596 1.180 ( -2.162 2.156 0.000) 3.054 1.188 ( -1.014 -0.156 0.000) 1.026 1.245 ( 0.725 -1.801 -0.000) 1.942 1.247 ( -0.453 -0.883 -0.000) 0.992 1.320 ( 0.837 1.387 0.000) 1.620 1.336 ( 0.420 0.550 0.000) 0.692 1.351 ( -2.419 2.433 0.000) 3.431 1.366 ( -1.594 1.319 0.000) 2.069 1.378 ( 2.398 1.952 -0.000) 3.092 1.414 ( 1.077 -1.675 -0.000) 1.991 1.446 ( 2.304 1.195 -0.000) 2.596 1.456 ( 0.046 -1.773 -0.000) 1.773 1.466 ( -2.577 -0.847 0.000) 2.712 1.536 ( -1.559 -3.391 -0.000) 3.732 1.579 ( 0.457 -0.253 -0.000) 0.523 1.579 ( -0.128 -0.308 -0.000) 0.334 1.635 ( 2.009 1.547 -0.000) 2.536 1.639 ( 3.736 -0.080 -0.000) 3.737 1.710 ( 0.221 -0.962 -0.000) 0.987 1.719 ( 0.050 -1.892 -0.000) 1.893 1.725 ( 0.812 0.888 0.000) 1.203 1.743 ( -1.844 -1.301 0.000) 2.257 1.760 ( 0.003 0.213 0.000) 0.213 1.771 ( 0.101 0.299 0.000) 0.316 1.800 ( -1.593 -3.746 -0.000) 4.070 1.819 ( 1.009 -1.133 -0.000) 1.517 1.822 ( 1.602 -2.353 -0.000) 2.846 1.835 ( -0.568 -1.505 -0.000) 1.608 1.846 ( 0.055 0.649 0.000) 0.652 1.859 ( -0.543 2.440 0.000) 2.500 1.905 ( -1.204 -1.008 0.000) 1.570 1.908 ( -1.556 2.609 0.000) 3.038 1.963 ( -1.964 -0.707 0.000) 2.088 1.985 ( -2.419 -0.514 0.000) 2.473 1.990 ( 2.672 0.564 -0.000) 2.731 2.025 ( -1.503 1.192 0.000) 1.918 2.050 ( -1.391 0.049 0.000) 1.392 2.092 ( -1.276 -1.100 0.000) 1.684 2.097 ( -1.348 2.883 0.000) 3.182 2.147 ( -2.088 1.101 0.000) 2.360 2.199 ( 1.961 0.611 -0.000) 2.054 2.221 ( 1.432 1.353 -0.000) 1.970 2.245 ( 1.028 0.471 -0.000) 1.131 2.259 ( 0.665 1.425 0.000) 1.573 2.265 ( 1.319 0.575 -0.000) 1.439 2.291 ( 1.320 0.495 -0.000) 1.409 2.324 ( -1.214 -0.303 0.000) 1.251 2.332 ( 1.973 -1.113 -0.000) 2.265 2.340 ( -1.237 0.310 0.000) 1.275 2.351 ( -0.644 0.057 0.000) 0.647 2.379 ( 0.467 1.059 0.000) 1.158 2.392 ( 1.606 0.959 -0.000) 1.871 2.426 ( -1.516 -0.885 0.000) 1.755 2.433 ( 0.050 -0.738 -0.000) 0.740 2.466 ( 0.615 -0.237 -0.000) 0.659 2.472 ( 0.135 0.037 -0.000) 0.139 2.723 ( 0.235 0.892 0.000) 0.923 2.738 ( -0.404 0.628 0.000) 0.747 2.751 ( -0.160 0.469 0.000) 0.496 2.754 ( 0.032 0.288 0.000) 0.290 2.796 ( -1.002 -1.143 -0.000) 1.520 2.798 ( -1.166 -1.038 0.000) 1.561 2.833 ( 0.496 -0.648 -0.000) 0.816 2.840 ( 0.303 -0.896 -0.000) 0.946 2.878 ( 0.433 0.098 -0.000) 0.444 2.894 ( -0.117 -0.713 -0.000) 0.722 2.905 ( 0.580 0.109 -0.000) 0.590 2.906 ( 0.234 0.338 0.000) 0.411 3.869 ( 0.298 0.536 0.000) 0.613 3.884 ( -0.245 -0.220 0.000) 0.329 3.885 ( 0.084 -0.131 -0.000) 0.156 3.889 ( -0.068 -0.363 -0.000) 0.369 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 15 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.188 ( 0.000 0.000 -10.604) 10.604 0.188 ( 0.000 0.000 -10.604) 10.604 0.352 ( 0.000 0.000 -19.769) 19.769 0.634 ( -0.000 -0.000 2.984) 2.984 0.634 ( -0.000 -0.000 2.984) 2.984 0.728 ( 0.000 0.000 -1.047) 1.047 0.728 ( 0.000 0.000 -1.047) 1.047 0.769 ( -0.000 -0.000 4.167) 4.167 0.790 ( -0.000 -0.000 2.931) 2.931 0.790 ( -0.000 -0.000 2.931) 2.931 0.807 ( -0.000 -0.000 3.654) 3.654 0.807 ( -0.000 -0.000 3.654) 3.654 0.869 ( -0.000 -0.000 0.500) 0.500 0.869 ( -0.000 -0.000 0.500) 0.500 0.881 ( 0.000 0.000 -1.997) 1.997 0.907 ( 0.000 0.000 -1.683) 1.683 0.910 ( -0.000 -0.000 6.258) 6.258 0.931 ( 0.000 0.000 -1.090) 1.090 0.932 ( 0.000 0.000 -1.220) 1.220 0.932 ( 0.000 0.000 -1.220) 1.220 0.940 ( -0.000 -0.000 0.779) 0.779 0.981 ( -0.000 -0.000 0.933) 0.933 0.986 ( 0.000 0.000 -2.704) 2.704 0.986 ( 0.000 0.000 -2.704) 2.704 1.194 ( -0.000 -0.000 0.694) 0.694 1.194 ( -0.000 -0.000 0.694) 0.694 1.259 ( 0.000 0.000 -0.963) 0.963 1.259 ( 0.000 0.000 -0.963) 0.963 1.262 ( -0.000 -0.000 1.185) 1.185 1.268 ( -0.000 -0.000 1.435) 1.435 1.341 ( 0.000 0.000 -2.692) 2.692 1.341 ( 0.000 0.000 -2.692) 2.692 1.345 ( 0.000 0.000 -1.072) 1.072 1.366 ( 0.000 0.000 -4.732) 4.732 1.459 ( -0.000 -0.000 0.751) 0.751 1.506 ( -0.000 -0.000 0.181) 0.181 1.552 ( -0.000 -0.000 0.967) 0.967 1.552 ( -0.000 -0.000 0.967) 0.967 1.564 ( -0.000 -0.000 1.556) 1.556 1.601 ( -0.000 -0.000 2.380) 2.380 1.601 ( -0.000 -0.000 2.380) 2.380 1.610 ( -0.000 -0.000 0.116) 0.116 1.674 ( 0.000 0.000 -0.706) 0.706 1.682 ( -0.000 -0.000 1.141) 1.141 1.715 ( -0.000 -0.000 1.425) 1.425 1.725 ( -0.000 -0.000 0.267) 0.267 1.787 ( -0.000 -0.000 2.469) 2.469 1.787 ( -0.000 -0.000 2.469) 2.469 1.803 ( -0.000 -0.000 0.446) 0.446 1.827 ( -0.000 -0.000 1.194) 1.194 1.827 ( -0.000 -0.000 1.194) 1.194 1.840 ( 0.000 0.000 -6.343) 6.343 1.908 ( 0.000 0.000 -1.846) 1.846 1.920 ( -0.000 -0.000 1.280) 1.280 1.927 ( -0.000 -0.000 0.290) 0.290 1.927 ( -0.000 -0.000 0.290) 0.290 2.003 ( -0.000 -0.000 4.975) 4.975 2.015 ( 0.000 0.000 -3.570) 3.570 2.072 ( -0.000 -0.000 1.441) 1.441 2.082 ( 0.000 0.000 -1.392) 1.392 2.082 ( 0.000 0.000 -1.392) 1.392 2.097 ( 0.000 0.000 -0.708) 0.708 2.097 ( 0.000 0.000 -0.708) 0.708 2.107 ( -0.000 -0.000 1.955) 1.955 2.181 ( 0.000 0.000 -1.085) 1.085 2.181 ( 0.000 0.000 -1.085) 1.085 2.184 ( 0.000 0.000 -1.556) 1.556 2.223 ( 0.000 0.000 -2.558) 2.558 2.223 ( 0.000 0.000 -2.558) 2.558 2.254 ( -0.000 -0.000 2.244) 2.244 2.295 ( -0.000 -0.000 0.049) 0.049 2.319 ( -0.000 -0.000 1.173) 1.173 2.357 ( 0.000 0.000 -0.157) 0.157 2.370 ( 0.000 0.000 -3.493) 3.493 2.372 ( 0.000 0.000 -1.377) 1.377 2.376 ( 0.000 0.000 -3.237) 3.237 2.422 ( -0.000 -0.000 0.932) 0.932 2.443 ( 0.000 0.000 -0.108) 0.108 2.443 ( -0.000 -0.000 1.642) 1.642 2.469 ( -0.000 -0.000 1.161) 1.161 2.702 ( 0.000 0.000 -0.586) 0.586 2.702 ( 0.000 0.000 -0.586) 0.586 2.735 ( 0.000 0.000 -1.295) 1.295 2.735 ( 0.000 0.000 -1.295) 1.295 2.838 ( -0.000 -0.000 1.835) 1.835 2.838 ( -0.000 -0.000 1.835) 1.835 2.852 ( -0.000 -0.000 1.090) 1.090 2.852 ( -0.000 -0.000 1.090) 1.090 2.875 ( -0.000 -0.000 0.329) 0.329 2.880 ( 0.000 0.000 -0.625) 0.625 2.883 ( -0.000 -0.000 0.370) 0.370 2.930 ( -0.000 -0.000 0.259) 0.259 3.854 ( 0.000 0.000 -0.503) 0.503 3.886 ( -0.000 -0.000 0.333) 0.333 3.889 ( -0.000 -0.000 0.111) 0.111 3.899 ( -0.000 -0.000 0.426) 0.426 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 21 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.534 ( -1.418 -1.418 -3.813) 4.308 0.534 ( 1.418 1.418 -3.813) 4.308 0.582 ( -0.307 -0.307 1.629) 1.686 0.582 ( 0.307 0.307 1.629) 1.686 0.656 ( -1.823 -1.823 -3.111) 4.040 0.656 ( 1.823 1.823 -3.111) 4.040 0.729 ( -0.031 -0.031 4.126) 4.126 0.729 ( 0.031 0.031 4.126) 4.126 0.762 ( -0.979 -0.979 -1.118) 1.780 0.762 ( 0.979 0.979 -1.118) 1.780 0.792 ( -0.789 -0.789 2.247) 2.509 0.792 ( 0.789 0.789 2.247) 2.509 0.842 ( -0.378 -0.378 -2.461) 2.518 0.842 ( 0.378 0.378 -2.461) 2.518 0.878 ( -0.354 -0.354 -0.645) 0.816 0.878 ( 0.354 0.354 -0.645) 0.816 0.910 ( -0.235 -0.235 1.748) 1.779 0.910 ( 0.235 0.235 1.748) 1.779 0.927 ( -0.898 -0.898 2.265) 2.596 0.927 ( 0.898 0.898 2.265) 2.596 0.954 ( -0.679 -0.679 -1.791) 2.032 0.954 ( 0.679 0.679 -1.791) 2.032 0.980 ( -1.479 -1.479 0.582) 2.172 0.980 ( 1.479 1.479 0.582) 2.172 1.179 ( -0.964 -0.964 -1.616) 2.114 1.179 ( 0.964 0.964 -1.616) 2.114 1.216 ( -0.084 -0.084 2.523) 2.526 1.216 ( 0.084 0.084 2.523) 2.526 1.347 ( -0.907 -0.907 -0.876) 1.554 1.347 ( 0.907 0.907 -0.876) 1.554 1.351 ( -0.379 -0.379 -1.013) 1.146 1.351 ( 0.379 0.379 -1.013) 1.146 1.416 ( -0.429 -0.429 -0.747) 0.962 1.416 ( 0.429 0.429 -0.747) 0.962 1.448 ( -1.899 -1.899 1.510) 3.081 1.448 ( 1.899 1.899 1.510) 3.081 1.508 ( -1.591 -1.591 -2.094) 3.073 1.508 ( 1.591 1.591 -2.094) 3.073 1.554 ( -0.409 -0.409 2.611) 2.674 1.554 ( 0.409 0.409 2.611) 2.674 1.677 ( -1.676 -1.676 0.080) 2.371 1.677 ( 1.676 1.676 0.080) 2.371 1.682 ( -1.357 -1.357 0.598) 2.010 1.682 ( 1.357 1.357 0.598) 2.010 1.760 ( -0.708 -0.708 -1.973) 2.213 1.760 ( 0.708 0.708 -1.973) 2.213 1.776 ( -1.273 -1.273 0.276) 1.821 1.776 ( 1.273 1.273 0.276) 1.821 1.813 ( -1.741 -1.741 -1.896) 3.107 1.813 ( 1.741 1.741 -1.896) 3.107 1.839 ( -0.180 -0.180 0.854) 0.891 1.839 ( 0.180 0.180 0.854) 0.891 1.858 ( -0.123 -0.123 -0.188) 0.256 1.858 ( 0.123 0.123 -0.188) 0.256 1.880 ( -1.327 -1.327 1.862) 2.644 1.880 ( 1.327 1.327 1.862) 2.644 1.934 ( -2.054 -2.054 -2.132) 3.602 1.934 ( 2.054 2.054 -2.132) 3.602 1.980 ( -0.061 -0.061 3.080) 3.081 1.980 ( 0.061 0.061 3.080) 3.081 2.084 ( -0.700 -0.700 -1.882) 2.127 2.084 ( 0.700 0.700 -1.882) 2.127 2.101 ( -1.855 -1.855 0.371) 2.650 2.101 ( 1.855 1.855 0.371) 2.650 2.241 ( -0.246 -0.246 -0.760) 0.836 2.241 ( 0.246 0.246 -0.760) 0.836 2.262 ( -0.680 -0.680 1.980) 2.201 2.262 ( 0.680 0.680 1.980) 2.201 2.282 ( -0.399 -0.399 -1.808) 1.894 2.282 ( 0.399 0.399 -1.808) 1.894 2.317 ( -0.801 -0.801 1.274) 1.705 2.317 ( 0.801 0.801 1.274) 1.705 2.354 ( -0.261 -0.261 -1.147) 1.205 2.354 ( 0.261 0.261 -1.147) 1.205 2.378 ( -1.103 -1.103 1.405) 2.100 2.378 ( 1.103 1.103 1.405) 2.100 2.443 ( -0.859 -0.859 -1.496) 1.927 2.443 ( 0.859 0.859 -1.496) 1.927 2.465 ( -0.001 -0.001 0.922) 0.922 2.465 ( 0.001 0.001 0.922) 0.922 2.738 ( -0.123 -0.123 -0.881) 0.898 2.738 ( 0.123 0.123 -0.881) 0.898 2.752 ( -0.315 -0.315 0.491) 0.663 2.752 ( 0.315 0.315 0.491) 0.663 2.793 ( -0.372 -0.372 -1.395) 1.491 2.793 ( 0.372 0.372 -1.395) 1.491 2.824 ( -0.325 -0.325 1.717) 1.777 2.824 ( 0.325 0.325 1.717) 1.777 2.895 ( -0.355 -0.355 -0.623) 0.800 2.895 ( 0.355 0.355 -0.623) 0.800 2.906 ( -0.013 -0.013 0.500) 0.500 2.906 ( 0.013 0.013 0.500) 0.500 3.879 ( -0.330 -0.330 -0.268) 0.538 3.879 ( 0.330 0.330 -0.268) 0.538 3.885 ( -0.090 -0.090 0.268) 0.296 3.885 ( 0.090 0.090 0.268) 0.296 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.432 ( 0.000 -6.133 -6.900) 9.231 0.470 ( 0.000 -8.386 -4.911) 9.719 0.594 ( -0.000 0.723 3.391) 3.467 0.606 ( 0.000 -1.961 0.603) 2.052 0.656 ( -0.000 0.594 4.944) 4.979 0.686 ( -0.000 4.107 0.178) 4.111 0.700 ( 0.000 -0.986 -1.480) 1.778 0.711 ( 0.000 -1.175 -5.672) 5.792 0.721 ( -0.000 0.445 2.334) 2.376 0.736 ( -0.000 0.587 1.716) 1.814 0.772 ( 0.000 -0.133 -0.393) 0.415 0.780 ( 0.000 -0.555 1.236) 1.355 0.795 ( -0.000 -0.427 4.049) 4.072 0.844 ( -0.000 1.097 -0.276) 1.131 0.868 ( 0.000 -0.302 0.001) 0.302 0.900 ( -0.000 0.548 -0.303) 0.626 0.911 ( 0.000 -0.196 -0.555) 0.589 0.913 ( -0.000 0.996 0.744) 1.243 0.934 ( -0.000 1.415 -1.119) 1.804 0.939 ( 0.000 0.264 -1.176) 1.205 0.960 ( 0.000 -1.307 1.346) 1.876 0.963 ( -0.000 -0.097 0.783) 0.789 1.003 ( -0.000 2.195 -1.697) 2.774 1.016 ( -0.000 1.918 -0.334) 1.947 1.169 ( 0.000 -0.258 0.041) 0.261 1.208 ( 0.000 -1.870 -0.011) 1.870 1.256 ( -0.000 0.378 -0.162) 0.411 1.266 ( -0.000 2.104 -1.347) 2.498 1.295 ( 0.000 -2.146 0.728) 2.266 1.308 ( 0.000 -2.610 -1.304) 2.917 1.338 ( -0.000 2.162 -1.671) 2.732 1.370 ( 0.000 0.570 -2.163) 2.237 1.395 ( 0.000 -0.887 -2.722) 2.863 1.431 ( -0.000 1.510 -3.058) 3.411 1.454 ( -0.000 3.175 2.643) 4.131 1.464 ( -0.000 0.285 3.438) 3.449 1.511 ( 0.000 -0.083 -0.040) 0.092 1.525 ( -0.000 0.444 0.975) 1.071 1.552 ( 0.000 -2.026 -0.119) 2.030 1.561 ( 0.000 -0.950 1.270) 1.586 1.606 ( -0.000 -0.769 2.833) 2.936 1.630 ( 0.000 -2.792 -0.066) 2.793 1.670 ( -0.000 2.594 0.020) 2.594 1.679 ( 0.000 0.224 -1.045) 1.069 1.712 ( 0.000 -1.188 1.253) 1.727 1.717 ( 0.000 -0.601 0.747) 0.959 1.747 ( 0.000 -3.363 1.834) 3.830 1.750 ( 0.000 -1.443 1.881) 2.370 1.783 ( 0.000 -2.159 2.776) 3.516 1.783 ( -0.000 0.190 0.920) 0.940 1.830 ( -0.000 -0.266 0.884) 0.924 1.868 ( -0.000 3.165 2.116) 3.808 1.896 ( -0.000 1.132 0.115) 1.138 1.899 ( -0.000 1.712 -2.978) 3.435 1.935 ( -0.000 2.068 -3.695) 4.234 1.978 ( 0.000 -0.561 -2.623) 2.682 1.983 ( -0.000 0.332 0.631) 0.713 1.987 ( -0.000 3.590 1.208) 3.788 2.002 ( -0.000 1.593 1.329) 2.074 2.054 ( -0.000 1.346 -0.878) 1.607 2.064 ( 0.000 -0.832 1.146) 1.416 2.096 ( 0.000 -0.555 0.767) 0.947 2.140 ( -0.000 0.926 -2.579) 2.740 2.161 ( 0.000 -2.400 -1.034) 2.614 2.189 ( 0.000 -1.022 -2.809) 2.989 2.210 ( 0.000 -1.059 0.078) 1.062 2.216 ( 0.000 -0.818 0.869) 1.194 2.226 ( -0.000 2.043 -2.153) 2.968 2.260 ( -0.000 0.838 -2.071) 2.234 2.275 ( -0.000 0.621 1.280) 1.423 2.316 ( 0.000 0.335 -1.065) 1.116 2.319 ( 0.000 -1.420 -1.198) 1.858 2.356 ( -0.000 0.332 -0.458) 0.566 2.388 ( -0.000 1.006 2.068) 2.299 2.398 ( -0.000 0.077 1.097) 1.100 2.412 ( 0.000 -0.921 -0.978) 1.344 2.417 ( -0.000 0.819 0.207) 0.845 2.429 ( 0.000 -0.196 -1.051) 1.069 2.441 ( 0.000 -0.811 0.053) 0.812 2.446 ( 0.000 -0.348 0.313) 0.469 2.729 ( 0.000 -0.401 -0.074) 0.408 2.732 ( 0.000 -0.929 -1.405) 1.684 2.770 ( 0.000 -0.828 -1.342) 1.577 2.772 ( 0.000 -1.030 -1.510) 1.828 2.790 ( -0.000 0.912 2.001) 2.199 2.792 ( -0.000 0.650 1.791) 1.905 2.821 ( -0.000 0.309 1.490) 1.522 2.825 ( -0.000 0.025 1.069) 1.070 2.875 ( 0.000 -0.238 0.426) 0.488 2.882 ( 0.000 -0.668 0.360) 0.758 2.893 ( 0.000 0.096 -0.829) 0.835 2.916 ( -0.000 0.516 0.286) 0.590 3.869 ( 0.000 -0.362 -0.673) 0.764 3.879 ( -0.000 0.244 0.244) 0.346 3.887 ( 0.000 -0.121 0.233) 0.263 3.888 ( -0.000 0.074 0.446) 0.452 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 26 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.415 ( -0.000 8.229 0.000) 8.229 0.488 ( -0.000 9.835 0.000) 9.835 0.636 ( 0.000 -5.273 -0.000) 5.273 0.646 ( -0.000 3.295 0.000) 3.295 0.646 ( 0.000 -3.513 -0.000) 3.513 0.682 ( 0.000 -0.985 -0.000) 0.985 0.686 ( 0.000 -3.954 -0.000) 3.954 0.711 ( 0.000 -3.569 -0.000) 3.569 0.717 ( 0.000 -0.420 -0.000) 0.420 0.761 ( 0.000 -0.534 -0.000) 0.534 0.787 ( 0.000 -1.261 -0.000) 1.261 0.805 ( 0.000 -2.677 -0.000) 2.677 0.842 ( -0.000 0.475 0.000) 0.475 0.859 ( -0.000 3.772 0.000) 3.772 0.863 ( -0.000 0.864 0.000) 0.864 0.866 ( 0.000 -0.381 -0.000) 0.381 0.891 ( 0.000 -0.496 -0.000) 0.496 0.914 ( -0.000 0.885 0.000) 0.885 0.924 ( -0.000 0.773 0.000) 0.773 0.926 ( -0.000 1.291 0.000) 1.291 0.941 ( 0.000 -0.139 -0.000) 0.139 0.963 ( -0.000 0.101 0.000) 0.101 0.974 ( -0.000 1.260 0.000) 1.260 0.990 ( 0.000 -2.023 -0.000) 2.023 1.199 ( 0.000 -1.800 -0.000) 1.800 1.232 ( -0.000 1.344 0.000) 1.344 1.233 ( -0.000 1.309 0.000) 1.309 1.239 ( -0.000 0.537 0.000) 0.537 1.306 ( -0.000 2.166 0.000) 2.166 1.321 ( 0.000 -0.832 -0.000) 0.832 1.351 ( -0.000 2.232 0.000) 2.232 1.386 ( -0.000 1.383 0.000) 1.383 1.392 ( -0.000 2.942 0.000) 2.942 1.403 ( -0.000 1.324 0.000) 1.324 1.406 ( 0.000 -4.820 -0.000) 4.820 1.463 ( 0.000 -2.089 -0.000) 2.089 1.531 ( 0.000 -1.174 -0.000) 1.174 1.535 ( -0.000 0.621 0.000) 0.621 1.569 ( 0.000 -2.721 -0.000) 2.721 1.586 ( 0.000 -1.615 -0.000) 1.615 1.586 ( 0.000 -0.424 -0.000) 0.424 1.606 ( -0.000 1.265 0.000) 1.265 1.698 ( 0.000 -1.378 -0.000) 1.378 1.707 ( -0.000 0.360 0.000) 0.360 1.711 ( 0.000 -2.477 -0.000) 2.477 1.743 ( 0.000 -0.029 -0.000) 0.029 1.779 ( -0.000 1.276 0.000) 1.276 1.781 ( 0.000 -2.350 -0.000) 2.350 1.798 ( -0.000 0.630 0.000) 0.630 1.802 ( 0.000 -2.912 -0.000) 2.912 1.825 ( -0.000 0.270 0.000) 0.270 1.832 ( 0.000 -3.628 -0.000) 3.628 1.867 ( -0.000 3.530 0.000) 3.530 1.873 ( -0.000 4.008 0.000) 4.008 1.923 ( 0.000 -5.361 -0.000) 5.361 1.936 ( -0.000 1.655 0.000) 1.655 1.959 ( 0.000 -0.300 -0.000) 0.300 1.979 ( 0.000 -1.912 -0.000) 1.912 1.982 ( -0.000 1.558 0.000) 1.558 2.065 ( -0.000 0.351 0.000) 0.351 2.071 ( -0.000 0.991 0.000) 0.991 2.100 ( 0.000 -0.011 -0.000) 0.011 2.150 ( -0.000 0.851 0.000) 0.851 2.151 ( -0.000 1.782 0.000) 1.782 2.182 ( -0.000 1.144 0.000) 1.144 2.188 ( -0.000 4.712 0.000) 4.712 2.242 ( -0.000 2.157 0.000) 2.157 2.254 ( 0.000 -1.231 -0.000) 1.231 2.260 ( 0.000 -0.023 -0.000) 0.023 2.261 ( 0.000 -2.260 -0.000) 2.260 2.286 ( -0.000 0.331 0.000) 0.331 2.340 ( -0.000 0.659 0.000) 0.659 2.355 ( -0.000 0.614 0.000) 0.614 2.366 ( -0.000 0.279 0.000) 0.279 2.371 ( -0.000 0.733 0.000) 0.733 2.371 ( -0.000 0.893 0.000) 0.893 2.433 ( -0.000 0.187 0.000) 0.187 2.447 ( 0.000 -0.507 -0.000) 0.507 2.455 ( 0.000 -0.979 -0.000) 0.979 2.468 ( 0.000 -0.289 -0.000) 0.289 2.720 ( -0.000 0.685 0.000) 0.685 2.744 ( -0.000 1.051 0.000) 1.051 2.753 ( -0.000 0.526 0.000) 0.526 2.756 ( -0.000 0.724 0.000) 0.724 2.815 ( 0.000 -1.375 -0.000) 1.375 2.821 ( 0.000 -1.564 -0.000) 1.564 2.821 ( 0.000 -1.260 -0.000) 1.260 2.829 ( 0.000 -1.095 -0.000) 1.095 2.874 ( 0.000 -0.132 -0.000) 0.132 2.894 ( -0.000 0.448 0.000) 0.448 2.895 ( 0.000 -0.842 -0.000) 0.842 2.903 ( -0.000 0.560 0.000) 0.560 3.866 ( -0.000 0.595 0.000) 0.595 3.884 ( 0.000 -0.225 -0.000) 0.225 3.887 ( 0.000 -0.149 -0.000) 0.149 3.890 ( 0.000 -0.474 -0.000) 0.474 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 20 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.530 ( 1.788 -1.788 0.000) 2.528 0.530 ( 1.788 -1.788 0.000) 2.528 0.576 ( -4.297 0.143 0.000) 4.299 0.576 ( -0.143 4.297 0.000) 4.299 0.654 ( 3.162 -3.986 0.000) 5.088 0.654 ( 3.986 -3.162 0.000) 5.088 0.678 ( -1.798 -0.414 0.000) 1.845 0.678 ( 0.414 1.798 0.000) 1.845 0.746 ( -0.347 -1.725 0.000) 1.759 0.746 ( 1.725 0.347 0.000) 1.759 0.779 ( 0.041 -1.053 -0.000) 1.053 0.779 ( 1.053 -0.041 -0.000) 1.053 0.841 ( -1.141 0.738 0.000) 1.359 0.841 ( -0.738 1.141 0.000) 1.359 0.883 ( -0.107 0.068 0.000) 0.126 0.883 ( -0.068 0.107 0.000) 0.126 0.912 ( -1.094 -0.140 0.000) 1.103 0.912 ( 0.140 1.094 0.000) 1.103 0.935 ( 0.365 -0.940 0.000) 1.008 0.935 ( 0.940 -0.365 0.000) 1.008 0.954 ( 0.328 -0.636 0.000) 0.716 0.954 ( 0.636 -0.328 0.000) 0.716 1.006 ( -1.806 1.246 0.000) 2.194 1.006 ( -1.246 1.806 0.000) 2.194 1.188 ( -0.153 -0.999 0.000) 1.011 1.188 ( 0.999 0.153 0.000) 1.011 1.245 ( -1.470 1.373 0.000) 2.011 1.245 ( -1.373 1.470 0.000) 2.011 1.335 ( 1.083 -1.597 0.000) 1.929 1.335 ( 1.597 -1.083 0.000) 1.929 1.391 ( -1.368 1.209 0.000) 1.826 1.391 ( -1.209 1.368 0.000) 1.826 1.405 ( -0.285 -1.370 0.000) 1.400 1.405 ( 1.370 0.285 0.000) 1.400 1.437 ( -2.996 -0.973 0.000) 3.150 1.437 ( 0.973 2.996 0.000) 3.150 1.478 ( -1.501 -1.119 0.000) 1.872 1.478 ( 1.119 1.501 0.000) 1.872 1.573 ( 0.440 -0.573 0.000) 0.722 1.573 ( 0.573 -0.440 0.000) 0.722 1.643 ( 0.612 -3.849 -0.000) 3.897 1.643 ( 3.849 -0.612 -0.000) 3.897 1.696 ( -1.024 -0.599 0.000) 1.186 1.696 ( 0.599 1.024 0.000) 1.186 1.716 ( 0.600 -0.806 0.000) 1.004 1.716 ( 0.806 -0.600 0.000) 1.004 1.724 ( 0.823 -1.734 -0.000) 1.919 1.724 ( 1.734 -0.823 -0.000) 1.919 1.765 ( 1.209 -1.777 0.000) 2.149 1.765 ( 1.777 -1.209 0.000) 2.149 1.821 ( -0.732 0.661 0.000) 0.987 1.821 ( -0.661 0.732 0.000) 0.987 1.897 ( -2.490 0.796 0.000) 2.614 1.897 ( -0.796 2.490 0.000) 2.614 1.939 ( -2.991 1.486 0.000) 3.339 1.939 ( -1.486 2.991 0.000) 3.339 1.993 ( -1.405 0.626 0.000) 1.538 1.993 ( -0.626 1.405 0.000) 1.538 2.028 ( -1.548 -1.059 0.000) 1.875 2.028 ( 1.059 1.548 0.000) 1.875 2.056 ( -1.581 -1.642 0.000) 2.279 2.056 ( 1.642 1.581 0.000) 2.279 2.156 ( -1.835 0.645 0.000) 1.945 2.156 ( -0.645 1.835 0.000) 1.945 2.209 ( -0.171 -1.278 0.000) 1.289 2.209 ( 1.278 0.171 0.000) 1.289 2.255 ( 0.371 -1.060 0.000) 1.123 2.255 ( 1.060 -0.371 0.000) 1.123 2.294 ( 0.453 -0.698 0.000) 0.833 2.294 ( 0.698 -0.453 0.000) 0.833 2.308 ( -1.871 -1.178 0.000) 2.212 2.308 ( 1.178 1.871 0.000) 2.212 2.360 ( -0.830 0.817 0.000) 1.164 2.360 ( -0.817 0.830 0.000) 1.164 2.395 ( 0.439 -0.511 0.000) 0.674 2.395 ( 0.511 -0.439 0.000) 0.674 2.419 ( -0.435 0.432 0.000) 0.613 2.419 ( -0.432 0.435 0.000) 0.613 2.461 ( 0.443 -0.869 0.000) 0.975 2.461 ( 0.869 -0.443 0.000) 0.975 2.740 ( -0.602 -0.673 0.000) 0.903 2.740 ( 0.673 0.602 0.000) 0.903 2.761 ( -0.687 0.552 0.000) 0.881 2.761 ( -0.552 0.687 0.000) 0.881 2.787 ( -0.097 0.091 0.000) 0.133 2.787 ( -0.091 0.097 0.000) 0.133 2.819 ( 0.826 -0.923 0.000) 1.238 2.819 ( 0.923 -0.826 0.000) 1.238 2.879 ( 0.065 -0.651 -0.000) 0.654 2.879 ( 0.651 -0.065 -0.000) 0.654 2.908 ( -0.044 -0.039 0.000) 0.059 2.908 ( 0.039 0.044 0.000) 0.059 3.879 ( -0.312 -0.138 0.000) 0.341 3.879 ( 0.138 0.312 0.000) 0.341 3.882 ( 0.176 -0.198 0.000) 0.264 3.882 ( 0.198 -0.176 0.000) 0.264 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 18 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.357 ( 0.000 0.000 -8.726) 8.726 0.357 ( 0.000 0.000 -8.726) 8.726 0.581 ( -0.000 -0.000 2.995) 2.995 0.581 ( -0.000 -0.000 2.995) 2.995 0.666 ( 0.000 0.000 -16.181) 16.181 0.699 ( -0.000 -0.000 3.811) 3.811 0.699 ( -0.000 -0.000 3.811) 3.811 0.702 ( -0.000 -0.000 2.880) 2.880 0.755 ( -0.000 -0.000 2.185) 2.185 0.755 ( -0.000 -0.000 2.185) 2.185 0.772 ( 0.000 0.000 -0.740) 0.740 0.772 ( 0.000 0.000 -0.740) 0.740 0.782 ( -0.000 -0.000 7.303) 7.303 0.863 ( -0.000 -0.000 0.338) 0.338 0.863 ( -0.000 -0.000 0.338) 0.338 0.895 ( 0.000 0.000 -0.457) 0.457 0.924 ( -0.000 -0.000 1.054) 1.054 0.931 ( -0.000 -0.000 0.925) 0.925 0.931 ( -0.000 -0.000 0.925) 0.925 0.945 ( 0.000 0.000 -0.600) 0.600 0.947 ( 0.000 0.000 -0.811) 0.811 0.967 ( -0.000 -0.000 0.783) 0.783 1.029 ( 0.000 0.000 -2.075) 2.075 1.029 ( 0.000 0.000 -2.075) 2.075 1.183 ( -0.000 -0.000 0.481) 0.481 1.183 ( -0.000 -0.000 0.481) 0.481 1.248 ( -0.000 -0.000 0.711) 0.711 1.249 ( -0.000 -0.000 0.465) 0.465 1.289 ( 0.000 0.000 -2.243) 2.243 1.289 ( 0.000 0.000 -2.243) 2.243 1.371 ( 0.000 0.000 -1.826) 1.826 1.388 ( 0.000 0.000 -2.436) 2.436 1.388 ( 0.000 0.000 -2.436) 2.436 1.438 ( -0.000 -0.000 1.661) 1.661 1.455 ( 0.000 0.000 -5.309) 5.309 1.505 ( 0.000 0.000 -0.107) 0.107 1.533 ( -0.000 -0.000 1.166) 1.166 1.533 ( -0.000 -0.000 1.166) 1.166 1.535 ( -0.000 -0.000 1.672) 1.672 1.554 ( -0.000 -0.000 2.576) 2.576 1.554 ( -0.000 -0.000 2.576) 2.576 1.595 ( -0.000 -0.000 2.031) 2.031 1.676 ( 0.000 0.000 -0.234) 0.234 1.685 ( 0.000 0.000 -0.505) 0.505 1.693 ( -0.000 -0.000 0.946) 0.946 1.717 ( -0.000 -0.000 1.042) 1.042 1.727 ( -0.000 -0.000 3.328) 3.328 1.727 ( -0.000 -0.000 3.328) 3.328 1.778 ( -0.000 -0.000 0.871) 0.871 1.811 ( -0.000 -0.000 1.098) 1.098 1.811 ( -0.000 -0.000 1.098) 1.098 1.902 ( -0.000 -0.000 6.225) 6.225 1.930 ( 0.000 0.000 -2.260) 2.260 1.942 ( 0.000 0.000 -5.137) 5.137 1.946 ( 0.000 0.000 -2.633) 2.633 1.946 ( 0.000 0.000 -2.633) 2.633 1.952 ( 0.000 0.000 -3.257) 3.257 2.034 ( -0.000 -0.000 2.796) 2.796 2.053 ( -0.000 -0.000 0.586) 0.586 2.067 ( -0.000 -0.000 2.195) 2.195 2.067 ( -0.000 -0.000 2.195) 2.195 2.068 ( -0.000 -0.000 2.177) 2.177 2.139 ( 0.000 0.000 -2.735) 2.735 2.139 ( 0.000 0.000 -2.735) 2.735 2.191 ( -0.000 -0.000 0.066) 0.066 2.191 ( -0.000 -0.000 0.066) 0.066 2.216 ( 0.000 0.000 -2.104) 2.104 2.247 ( 0.000 0.000 -2.205) 2.205 2.273 ( 0.000 0.000 -2.711) 2.711 2.273 ( 0.000 0.000 -2.711) 2.711 2.292 ( -0.000 -0.000 1.958) 1.958 2.318 ( 0.000 0.000 -2.475) 2.475 2.363 ( 0.000 0.000 -0.536) 0.536 2.398 ( -0.000 -0.000 3.065) 3.065 2.400 ( 0.000 0.000 -1.705) 1.705 2.400 ( -0.000 -0.000 1.461) 1.461 2.421 ( 0.000 0.000 -2.155) 2.155 2.426 ( 0.000 0.000 -1.930) 1.930 2.438 ( -0.000 -0.000 2.357) 2.357 2.442 ( -0.000 -0.000 0.439) 0.439 2.720 ( 0.000 0.000 -1.217) 1.217 2.720 ( 0.000 0.000 -1.217) 1.217 2.760 ( 0.000 0.000 -1.472) 1.472 2.760 ( 0.000 0.000 -1.472) 1.472 2.801 ( -0.000 -0.000 2.092) 2.092 2.801 ( -0.000 -0.000 2.092) 2.092 2.826 ( -0.000 -0.000 1.772) 1.772 2.826 ( -0.000 -0.000 1.772) 1.772 2.869 ( -0.000 -0.000 0.241) 0.241 2.875 ( -0.000 -0.000 0.483) 0.483 2.894 ( 0.000 0.000 -0.942) 0.942 2.923 ( -0.000 -0.000 0.581) 0.581 3.865 ( 0.000 0.000 -0.724) 0.724 3.881 ( -0.000 -0.000 0.179) 0.179 3.885 ( -0.000 -0.000 0.332) 0.332 3.890 ( -0.000 -0.000 0.569) 0.569 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 20 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.514 ( 0.000 -5.113 -0.000) 5.113 0.518 ( 0.000 -5.431 -0.000) 5.431 0.567 ( -0.000 6.084 0.000) 6.084 0.585 ( 0.000 -3.594 -0.000) 3.594 0.604 ( 0.000 -5.062 -0.000) 5.062 0.657 ( -0.000 1.471 0.000) 1.471 0.670 ( -0.000 4.194 0.000) 4.194 0.693 ( -0.000 0.118 0.000) 0.118 0.720 ( 0.000 -0.824 -0.000) 0.824 0.756 ( 0.000 -1.353 -0.000) 1.353 0.774 ( -0.000 0.779 0.000) 0.779 0.777 ( -0.000 0.186 0.000) 0.186 0.844 ( -0.000 0.332 0.000) 0.332 0.872 ( 0.000 -1.175 -0.000) 1.175 0.880 ( -0.000 0.945 0.000) 0.945 0.887 ( -0.000 0.589 0.000) 0.589 0.905 ( -0.000 1.107 0.000) 1.107 0.906 ( 0.000 -0.605 -0.000) 0.605 0.934 ( 0.000 -0.765 -0.000) 0.765 0.938 ( 0.000 -2.548 -0.000) 2.548 0.944 ( -0.000 0.914 0.000) 0.914 0.952 ( -0.000 0.283 0.000) 0.283 1.020 ( -0.000 2.325 0.000) 2.325 1.020 ( -0.000 3.168 0.000) 3.168 1.169 ( 0.000 -0.391 -0.000) 0.391 1.208 ( 0.000 -2.192 -0.000) 2.192 1.257 ( -0.000 0.142 0.000) 0.142 1.276 ( -0.000 2.685 0.000) 2.685 1.287 ( 0.000 -2.448 -0.000) 2.448 1.324 ( 0.000 -2.072 -0.000) 2.072 1.381 ( -0.000 2.945 0.000) 2.945 1.396 ( -0.000 1.466 0.000) 1.466 1.406 ( -0.000 0.659 0.000) 0.659 1.422 ( 0.000 -0.939 -0.000) 0.939 1.426 ( -0.000 1.749 0.000) 1.749 1.467 ( -0.000 3.487 0.000) 3.487 1.510 ( 0.000 -0.395 -0.000) 0.395 1.515 ( 0.000 -0.785 -0.000) 0.785 1.546 ( 0.000 -0.441 -0.000) 0.441 1.548 ( 0.000 -1.559 -0.000) 1.559 1.573 ( 0.000 -0.990 -0.000) 0.990 1.634 ( 0.000 -1.737 -0.000) 1.737 1.686 ( -0.000 0.048 0.000) 0.048 1.692 ( 0.000 -0.991 -0.000) 0.991 1.697 ( 0.000 -0.223 -0.000) 0.223 1.699 ( 0.000 -1.619 -0.000) 1.619 1.732 ( 0.000 -2.891 -0.000) 2.891 1.732 ( 0.000 -2.526 -0.000) 2.526 1.748 ( 0.000 -1.308 -0.000) 1.308 1.773 ( 0.000 -0.207 -0.000) 0.207 1.818 ( 0.000 -0.617 -0.000) 0.617 1.846 ( -0.000 1.047 0.000) 1.047 1.904 ( -0.000 3.466 0.000) 3.466 1.948 ( -0.000 2.924 0.000) 2.924 1.952 ( -0.000 3.377 0.000) 3.377 1.966 ( -0.000 0.810 0.000) 0.810 1.977 ( -0.000 2.245 0.000) 2.245 2.001 ( -0.000 2.810 0.000) 2.810 2.010 ( -0.000 0.611 0.000) 0.611 2.048 ( 0.000 -2.789 -0.000) 2.789 2.054 ( 0.000 -1.385 -0.000) 1.385 2.081 ( 0.000 -2.012 -0.000) 2.012 2.171 ( -0.000 0.158 0.000) 0.158 2.175 ( -0.000 1.119 0.000) 1.119 2.194 ( 0.000 -0.506 -0.000) 0.506 2.203 ( 0.000 -1.806 -0.000) 1.806 2.240 ( 0.000 -0.736 -0.000) 0.736 2.242 ( -0.000 0.369 0.000) 0.369 2.281 ( -0.000 2.463 0.000) 2.463 2.285 ( -0.000 1.850 0.000) 1.850 2.301 ( -0.000 0.993 0.000) 0.993 2.332 ( 0.000 -1.449 -0.000) 1.449 2.370 ( -0.000 0.694 0.000) 0.694 2.371 ( -0.000 0.367 0.000) 0.367 2.383 ( -0.000 0.362 0.000) 0.362 2.386 ( -0.000 0.302 0.000) 0.302 2.434 ( 0.000 -0.624 -0.000) 0.624 2.435 ( 0.000 -0.750 -0.000) 0.750 2.438 ( -0.000 0.349 0.000) 0.349 2.453 ( 0.000 -0.852 -0.000) 0.852 2.731 ( -0.000 0.275 0.000) 0.275 2.747 ( 0.000 -0.952 -0.000) 0.952 2.768 ( -0.000 1.046 0.000) 1.046 2.772 ( -0.000 0.676 0.000) 0.676 2.792 ( 0.000 -1.009 -0.000) 1.009 2.795 ( 0.000 -1.123 -0.000) 1.123 2.799 ( 0.000 -0.044 -0.000) 0.044 2.809 ( 0.000 -0.698 -0.000) 0.698 2.871 ( 0.000 -0.134 -0.000) 0.134 2.878 ( 0.000 -0.713 -0.000) 0.713 2.905 ( -0.000 0.443 0.000) 0.443 2.910 ( -0.000 0.079 0.000) 0.079 3.877 ( -0.000 0.410 0.000) 0.410 3.879 ( 0.000 -0.192 -0.000) 0.192 3.881 ( 0.000 -0.333 -0.000) 0.333 3.883 ( 0.000 -0.140 -0.000) 0.140 ======================= 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 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 1.94e-05 Number of triplets: 12 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.454 ( 0.000 0.000 -0.000) 0.000 0.454 ( 0.000 0.000 -0.000) 0.000 0.543 ( 0.000 0.000 -0.000) 0.000 0.543 ( 0.000 0.000 -0.000) 0.000 0.655 ( 0.000 0.000 -0.000) 0.000 0.655 ( 0.000 0.000 -0.000) 0.000 0.694 ( 0.000 0.000 -2.400) 2.400 0.694 ( -0.000 -0.000 2.400) 2.400 0.735 ( 0.000 0.000 -0.000) 0.000 0.735 ( 0.000 0.000 -0.000) 0.000 0.783 ( 0.000 0.000 -0.000) 0.000 0.783 ( 0.000 0.000 -0.000) 0.000 0.857 ( 0.000 0.000 -0.000) 0.000 0.857 ( 0.000 0.000 -0.000) 0.000 0.883 ( 0.000 0.000 -7.832) 7.832 0.883 ( -0.000 -0.000 7.832) 7.832 0.909 ( 0.000 0.000 -0.156) 0.156 0.909 ( -0.000 -0.000 0.156) 0.156 0.920 ( 0.000 0.000 -0.000) 0.000 0.920 ( 0.000 0.000 -0.000) 0.000 0.955 ( 0.000 0.000 -0.629) 0.629 0.955 ( -0.000 -0.000 0.629) 0.629 1.048 ( 0.000 0.000 -0.000) 0.000 1.048 ( 0.000 0.000 -0.000) 0.000 1.179 ( 0.000 0.000 -0.000) 0.000 1.179 ( 0.000 0.000 -0.000) 0.000 1.244 ( 0.000 0.000 -0.131) 0.131 1.244 ( -0.000 -0.000 0.131) 0.131 1.313 ( 0.000 0.000 -0.000) 0.000 1.313 ( 0.000 0.000 -0.000) 0.000 1.405 ( 0.000 0.000 -2.075) 2.075 1.405 ( -0.000 -0.000 2.075) 2.075 1.413 ( -0.000 -0.000 0.000) 0.000 1.413 ( -0.000 -0.000 0.000) 0.000 1.512 ( 0.000 0.000 -0.848) 0.848 1.512 ( -0.000 -0.000 0.848) 0.848 1.514 ( 0.000 0.000 -0.000) 0.000 1.514 ( 0.000 0.000 -0.000) 0.000 1.534 ( 0.000 0.000 -0.000) 0.000 1.534 ( 0.000 0.000 -0.000) 0.000 1.539 ( 0.000 0.000 -4.279) 4.279 1.539 ( -0.000 -0.000 4.279) 4.279 1.682 ( 0.000 0.000 -0.449) 0.449 1.682 ( -0.000 -0.000 0.449) 0.449 1.695 ( 0.000 0.000 -0.884) 0.884 1.695 ( -0.000 -0.000 0.884) 0.884 1.696 ( 0.000 0.000 -0.000) 0.000 1.696 ( 0.000 0.000 -0.000) 0.000 1.799 ( 0.000 0.000 -0.000) 0.000 1.799 ( 0.000 0.000 -0.000) 0.000 1.808 ( 0.000 0.000 -4.071) 4.071 1.808 ( -0.000 -0.000 4.071) 4.071 1.980 ( 0.000 0.000 -3.252) 3.252 1.980 ( -0.000 -0.000 3.252) 3.252 1.991 ( 0.000 0.000 -0.000) 0.000 1.991 ( 0.000 0.000 -0.000) 0.000 2.010 ( 0.000 0.000 -3.324) 3.324 2.010 ( -0.000 -0.000 3.324) 3.324 2.020 ( 0.000 0.000 -3.763) 3.763 2.020 ( -0.000 -0.000 3.763) 3.763 2.027 ( 0.000 0.000 -0.000) 0.000 2.027 ( 0.000 0.000 -0.000) 0.000 2.174 ( -0.000 -0.000 0.000) 0.000 2.174 ( -0.000 -0.000 0.000) 0.000 2.181 ( 0.000 0.000 -0.000) 0.000 2.181 ( 0.000 0.000 -0.000) 0.000 2.254 ( 0.000 0.000 -2.276) 2.276 2.254 ( -0.000 -0.000 2.276) 2.276 2.299 ( 0.000 0.000 -0.000) 0.000 2.299 ( 0.000 0.000 -0.000) 0.000 2.314 ( 0.000 0.000 -4.281) 4.281 2.314 ( -0.000 -0.000 4.281) 4.281 2.377 ( 0.000 0.000 -1.139) 1.139 2.377 ( -0.000 -0.000 1.139) 1.139 2.387 ( 0.000 0.000 -2.957) 2.957 2.387 ( -0.000 -0.000 2.957) 2.957 2.427 ( 0.000 0.000 -1.271) 1.271 2.427 ( -0.000 -0.000 1.271) 1.271 2.443 ( 0.000 0.000 -0.100) 0.100 2.443 ( -0.000 -0.000 0.100) 0.100 2.733 ( 0.000 0.000 -0.000) 0.000 2.733 ( 0.000 0.000 -0.000) 0.000 2.780 ( 0.000 0.000 -0.000) 0.000 2.780 ( 0.000 0.000 -0.000) 0.000 2.781 ( 0.000 0.000 -0.000) 0.000 2.781 ( 0.000 0.000 -0.000) 0.000 2.801 ( 0.000 0.000 -0.000) 0.000 2.801 ( 0.000 0.000 -0.000) 0.000 2.869 ( 0.000 0.000 -0.206) 0.206 2.869 ( -0.000 -0.000 0.206) 0.206 2.910 ( 0.000 0.000 -0.877) 0.877 2.910 ( -0.000 -0.000 0.877) 0.877 3.877 ( 0.000 0.000 -0.607) 0.607 3.877 ( -0.000 -0.000 0.607) 0.607 3.882 ( 0.000 0.000 -0.292) 0.292 3.882 ( -0.000 -0.000 0.292) 0.292 =================== 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 15.637 15.637 20.217 0.000 -0.000 -0.000 3/9216 20.0 3.340 3.340 4.071 0.000 -0.000 -0.000 3/9216 30.0 1.684 1.684 2.012 0.000 -0.000 -0.000 3/9216 40.0 1.127 1.127 1.337 0.000 -0.000 -0.000 3/9216 50.0 0.853 0.853 1.009 0.000 -0.000 -0.000 3/9216 60.0 0.689 0.689 0.815 0.000 -0.000 -0.000 3/9216 70.0 0.580 0.580 0.685 0.000 -0.000 -0.000 3/9216 80.0 0.501 0.501 0.592 0.000 -0.000 -0.000 3/9216 90.0 0.442 0.442 0.522 0.000 -0.000 -0.000 3/9216 100.0 0.395 0.395 0.467 0.000 -0.000 -0.000 3/9216 110.0 0.358 0.358 0.423 0.000 -0.000 -0.000 3/9216 120.0 0.327 0.327 0.386 0.000 -0.000 -0.000 3/9216 130.0 0.301 0.301 0.356 0.000 -0.000 -0.000 3/9216 140.0 0.279 0.279 0.330 0.000 -0.000 -0.000 3/9216 150.0 0.260 0.260 0.307 0.000 -0.000 -0.000 3/9216 160.0 0.243 0.243 0.287 0.000 -0.000 -0.000 3/9216 170.0 0.229 0.229 0.270 0.000 -0.000 -0.000 3/9216 180.0 0.216 0.216 0.255 0.000 -0.000 -0.000 3/9216 190.0 0.204 0.204 0.241 0.000 -0.000 -0.000 3/9216 200.0 0.194 0.194 0.229 0.000 -0.000 -0.000 3/9216 210.0 0.185 0.185 0.218 0.000 -0.000 -0.000 3/9216 220.0 0.176 0.176 0.208 0.000 -0.000 -0.000 3/9216 230.0 0.168 0.168 0.199 0.000 -0.000 -0.000 3/9216 240.0 0.161 0.161 0.191 0.000 -0.000 -0.000 3/9216 250.0 0.155 0.155 0.183 0.000 -0.000 -0.000 3/9216 260.0 0.149 0.149 0.176 0.000 -0.000 -0.000 3/9216 270.0 0.143 0.143 0.169 0.000 -0.000 -0.000 3/9216 280.0 0.138 0.138 0.163 0.000 -0.000 -0.000 3/9216 290.0 0.133 0.133 0.158 0.000 -0.000 -0.000 3/9216 300.0 0.129 0.129 0.152 0.000 -0.000 -0.000 3/9216 310.0 0.125 0.125 0.147 0.000 -0.000 -0.000 3/9216 320.0 0.121 0.121 0.143 0.000 -0.000 -0.000 3/9216 330.0 0.117 0.117 0.138 0.000 -0.000 -0.000 3/9216 340.0 0.114 0.114 0.134 0.000 -0.000 -0.000 3/9216 350.0 0.110 0.110 0.130 0.000 -0.000 -0.000 3/9216 360.0 0.107 0.107 0.127 0.000 -0.000 -0.000 3/9216 370.0 0.104 0.104 0.123 0.000 -0.000 -0.000 3/9216 380.0 0.102 0.102 0.120 0.000 -0.000 -0.000 3/9216 390.0 0.099 0.099 0.117 0.000 -0.000 -0.000 3/9216 400.0 0.097 0.097 0.114 0.000 -0.000 -0.000 3/9216 410.0 0.094 0.094 0.111 0.000 -0.000 -0.000 3/9216 420.0 0.092 0.092 0.109 0.000 -0.000 -0.000 3/9216 430.0 0.090 0.090 0.106 0.000 -0.000 -0.000 3/9216 440.0 0.088 0.088 0.104 0.000 -0.000 -0.000 3/9216 450.0 0.086 0.086 0.101 0.000 -0.000 -0.000 3/9216 460.0 0.084 0.084 0.099 0.000 -0.000 -0.000 3/9216 470.0 0.082 0.082 0.097 0.000 -0.000 -0.000 3/9216 480.0 0.080 0.080 0.095 0.000 -0.000 -0.000 3/9216 490.0 0.079 0.079 0.093 0.000 -0.000 -0.000 3/9216 500.0 0.077 0.077 0.091 0.000 -0.000 -0.000 3/9216 510.0 0.076 0.076 0.089 0.000 -0.000 -0.000 3/9216 520.0 0.074 0.074 0.088 0.000 -0.000 -0.000 3/9216 530.0 0.073 0.073 0.086 0.000 -0.000 -0.000 3/9216 540.0 0.071 0.071 0.084 0.000 -0.000 -0.000 3/9216 550.0 0.070 0.070 0.083 0.000 -0.000 -0.000 3/9216 560.0 0.069 0.069 0.081 0.000 -0.000 -0.000 3/9216 570.0 0.068 0.068 0.080 0.000 -0.000 -0.000 3/9216 580.0 0.067 0.067 0.079 0.000 -0.000 -0.000 3/9216 590.0 0.065 0.065 0.077 0.000 -0.000 -0.000 3/9216 600.0 0.064 0.064 0.076 0.000 -0.000 -0.000 3/9216 610.0 0.063 0.063 0.075 0.000 -0.000 -0.000 3/9216 620.0 0.062 0.062 0.074 0.000 -0.000 -0.000 3/9216 630.0 0.061 0.061 0.072 0.000 -0.000 -0.000 3/9216 640.0 0.060 0.060 0.071 0.000 -0.000 -0.000 3/9216 650.0 0.059 0.059 0.070 0.000 -0.000 -0.000 3/9216 660.0 0.058 0.058 0.069 0.000 -0.000 -0.000 3/9216 670.0 0.058 0.058 0.068 0.000 -0.000 -0.000 3/9216 680.0 0.057 0.057 0.067 0.000 -0.000 -0.000 3/9216 690.0 0.056 0.056 0.066 0.000 -0.000 -0.000 3/9216 700.0 0.055 0.055 0.065 0.000 -0.000 -0.000 3/9216 710.0 0.054 0.054 0.064 0.000 -0.000 -0.000 3/9216 720.0 0.054 0.054 0.063 0.000 -0.000 -0.000 3/9216 730.0 0.053 0.053 0.062 0.000 -0.000 -0.000 3/9216 740.0 0.052 0.052 0.062 0.000 -0.000 -0.000 3/9216 750.0 0.051 0.051 0.061 0.000 -0.000 -0.000 3/9216 760.0 0.051 0.051 0.060 0.000 -0.000 -0.000 3/9216 770.0 0.050 0.050 0.059 0.000 -0.000 -0.000 3/9216 780.0 0.049 0.049 0.058 0.000 -0.000 -0.000 3/9216 790.0 0.049 0.049 0.058 0.000 -0.000 -0.000 3/9216 800.0 0.048 0.048 0.057 0.000 -0.000 -0.000 3/9216 810.0 0.048 0.048 0.056 0.000 -0.000 -0.000 3/9216 820.0 0.047 0.047 0.056 0.000 -0.000 -0.000 3/9216 830.0 0.046 0.046 0.055 0.000 -0.000 -0.000 3/9216 840.0 0.046 0.046 0.054 0.000 -0.000 -0.000 3/9216 850.0 0.045 0.045 0.054 0.000 -0.000 -0.000 3/9216 860.0 0.045 0.045 0.053 0.000 -0.000 -0.000 3/9216 870.0 0.044 0.044 0.052 0.000 -0.000 -0.000 3/9216 880.0 0.044 0.044 0.052 0.000 -0.000 -0.000 3/9216 890.0 0.043 0.043 0.051 0.000 -0.000 -0.000 3/9216 900.0 0.043 0.043 0.051 0.000 -0.000 -0.000 3/9216 910.0 0.042 0.042 0.050 0.000 -0.000 -0.000 3/9216 920.0 0.042 0.042 0.050 0.000 -0.000 -0.000 3/9216 930.0 0.041 0.041 0.049 0.000 -0.000 -0.000 3/9216 940.0 0.041 0.041 0.049 0.000 -0.000 -0.000 3/9216 950.0 0.041 0.041 0.048 0.000 -0.000 -0.000 3/9216 960.0 0.040 0.040 0.048 0.000 -0.000 -0.000 3/9216 970.0 0.040 0.040 0.047 0.000 -0.000 -0.000 3/9216 980.0 0.039 0.039 0.047 0.000 -0.000 -0.000 3/9216 990.0 0.039 0.039 0.046 0.000 -0.000 -0.000 3/9216 1000.0 0.039 0.039 0.046 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-08 04:59:01]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|