----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 09:58:04]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [2 2 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pnma (62) Number of symmetry operations in supercell: 32 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 9.434285200000000 0.000000000000000 0.000000000000000 b 0.000000000000000 7.381102380000000 0.000000000000000 c 0.000000000000000 0.000000000000000 15.111277439999999 Atomic positions (fractional): *1 Cs 0.96992290148160 0.25000000000000 0.20911211981328 132.905 2 Cs 0.03007709851840 0.75000000000000 0.79088788018672 132.905 3 Cs 0.46992290148160 0.25000000000000 0.29088788018672 132.905 4 Cs 0.53007709851840 0.75000000000000 0.70911211981328 132.905 *5 Cs 0.83096926700783 0.75000000000000 0.41534229816506 132.905 6 Cs 0.16903073299217 0.25000000000000 0.58465770183495 132.905 7 Cs 0.33096926700783 0.75000000000000 0.08465770183495 132.905 8 Cs 0.66903073299217 0.25000000000000 0.91534229816505 132.905 *9 Dy 0.81742696229724 0.75000000000000 0.06528498301047 162.500 10 Dy 0.18257303770276 0.25000000000000 0.93471501698953 162.500 11 Dy 0.31742696229724 0.75000000000000 0.43471501698953 162.500 12 Dy 0.68257303770276 0.25000000000000 0.56528498301047 162.500 *13 Cl 0.00000000000000 0.00000000000000 0.00000000000000 35.453 14 Cl 0.00000000000000 0.50000000000000 0.00000000000000 35.453 15 Cl 0.50000000000000 0.50000000000000 0.50000000000000 35.453 16 Cl 0.50000000000000 0.00000000000000 0.50000000000000 35.453 *17 Cl 0.65910274759175 0.49225769358265 0.11782660982414 35.453 18 Cl 0.34089725240825 0.99225769358265 0.88217339017586 35.453 19 Cl 0.15910274759175 0.00774230641735 0.38217339017586 35.453 20 Cl 0.84089725240825 0.50774230641735 0.61782660982414 35.453 21 Cl 0.34089725240825 0.50774230641735 0.88217339017586 35.453 22 Cl 0.65910274759175 0.00774230641735 0.11782660982414 35.453 23 Cl 0.84089725240825 0.99225769358265 0.61782660982414 35.453 24 Cl 0.15910274759175 0.49225769358265 0.38217339017586 35.453 *25 Cl 0.78847484992322 0.25000000000000 0.40531260520702 35.453 26 Cl 0.21152515007678 0.75000000000000 0.59468739479298 35.453 27 Cl 0.28847484992322 0.25000000000000 0.09468739479298 35.453 28 Cl 0.71152515007678 0.75000000000000 0.90531260520702 35.453 *29 Cl 0.97982617717568 0.75000000000000 0.20226779441939 35.453 30 Cl 0.02017382282432 0.25000000000000 0.79773220558061 35.453 31 Cl 0.47982617717568 0.75000000000000 0.29773220558061 35.453 32 Cl 0.52017382282432 0.25000000000000 0.70226779441939 35.453 -------------------------------- unit cell --------------------------------- Lattice vectors: a 9.434285200000000 0.000000000000000 0.000000000000000 b 0.000000000000000 7.381102380000000 0.000000000000000 c 0.000000000000000 0.000000000000000 15.111277439999999 Atomic positions (fractional): *1 Cs 0.96992290148160 0.25000000000000 0.20911211981328 132.905 > 1 2 Cs 0.03007709851840 0.75000000000000 0.79088788018672 132.905 > 2 3 Cs 0.46992290148160 0.25000000000000 0.29088788018672 132.905 > 3 4 Cs 0.53007709851840 0.75000000000000 0.70911211981328 132.905 > 4 *5 Cs 0.83096926700783 0.75000000000000 0.41534229816506 132.905 > 5 6 Cs 0.16903073299217 0.25000000000000 0.58465770183495 132.905 > 6 7 Cs 0.33096926700783 0.75000000000000 0.08465770183495 132.905 > 7 8 Cs 0.66903073299217 0.25000000000000 0.91534229816505 132.905 > 8 *9 Dy 0.81742696229724 0.75000000000000 0.06528498301047 162.500 > 9 10 Dy 0.18257303770276 0.25000000000000 0.93471501698953 162.500 > 10 11 Dy 0.31742696229724 0.75000000000000 0.43471501698953 162.500 > 11 12 Dy 0.68257303770276 0.25000000000000 0.56528498301047 162.500 > 12 *13 Cl 0.00000000000000 0.00000000000000 0.00000000000000 35.453 > 13 14 Cl 0.00000000000000 0.50000000000000 0.00000000000000 35.453 > 14 15 Cl 0.50000000000000 0.50000000000000 0.50000000000000 35.453 > 15 16 Cl 0.50000000000000 0.00000000000000 0.50000000000000 35.453 > 16 *17 Cl 0.65910274759175 0.49225769358265 0.11782660982414 35.453 > 17 18 Cl 0.34089725240825 0.99225769358265 0.88217339017586 35.453 > 18 19 Cl 0.15910274759175 0.00774230641735 0.38217339017586 35.453 > 19 20 Cl 0.84089725240825 0.50774230641735 0.61782660982414 35.453 > 20 21 Cl 0.34089725240825 0.50774230641735 0.88217339017586 35.453 > 21 22 Cl 0.65910274759175 0.00774230641735 0.11782660982414 35.453 > 22 23 Cl 0.84089725240825 0.99225769358265 0.61782660982414 35.453 > 23 24 Cl 0.15910274759175 0.49225769358265 0.38217339017586 35.453 > 24 *25 Cl 0.78847484992322 0.25000000000000 0.40531260520702 35.453 > 25 26 Cl 0.21152515007678 0.75000000000000 0.59468739479298 35.453 > 26 27 Cl 0.28847484992322 0.25000000000000 0.09468739479298 35.453 > 27 28 Cl 0.71152515007678 0.75000000000000 0.90531260520702 35.453 > 28 *29 Cl 0.97982617717568 0.75000000000000 0.20226779441939 35.453 > 29 30 Cl 0.02017382282432 0.25000000000000 0.79773220558061 35.453 > 30 31 Cl 0.47982617717568 0.75000000000000 0.29773220558061 35.453 > 31 32 Cl 0.52017382282432 0.25000000000000 0.70226779441939 35.453 > 32 -------------------------------- super cell -------------------------------- Lattice vectors: a 18.868570399999999 0.000000000000000 0.000000000000000 b 0.000000000000000 14.762204759999999 0.000000000000000 c 0.000000000000000 0.000000000000000 15.111277439999999 Atomic positions (fractional): *1 Cs 0.48496145074080 0.12500000000000 0.20911211981328 132.905 > 1 2 Cs 0.98496145074080 0.12500000000000 0.20911211981328 132.905 > 1 3 Cs 0.48496145074080 0.62500000000000 0.20911211981328 132.905 > 1 4 Cs 0.98496145074080 0.62500000000000 0.20911211981328 132.905 > 1 5 Cs 0.01503854925920 0.37500000000000 0.79088788018672 132.905 > 2 6 Cs 0.51503854925920 0.37500000000000 0.79088788018672 132.905 > 2 7 Cs 0.01503854925920 0.87500000000000 0.79088788018672 132.905 > 2 8 Cs 0.51503854925920 0.87500000000000 0.79088788018672 132.905 > 2 9 Cs 0.23496145074080 0.12500000000000 0.29088788018672 132.905 > 3 10 Cs 0.73496145074080 0.12500000000000 0.29088788018672 132.905 > 3 11 Cs 0.23496145074080 0.62500000000000 0.29088788018672 132.905 > 3 12 Cs 0.73496145074080 0.62500000000000 0.29088788018672 132.905 > 3 13 Cs 0.26503854925920 0.37500000000000 0.70911211981328 132.905 > 4 14 Cs 0.76503854925920 0.37500000000000 0.70911211981328 132.905 > 4 15 Cs 0.26503854925920 0.87500000000000 0.70911211981328 132.905 > 4 16 Cs 0.76503854925920 0.87500000000000 0.70911211981328 132.905 > 4 *17 Cs 0.41548463350392 0.37500000000000 0.41534229816506 132.905 > 5 18 Cs 0.91548463350392 0.37500000000000 0.41534229816506 132.905 > 5 19 Cs 0.41548463350392 0.87500000000000 0.41534229816506 132.905 > 5 20 Cs 0.91548463350392 0.87500000000000 0.41534229816506 132.905 > 5 21 Cs 0.08451536649608 0.12500000000000 0.58465770183495 132.905 > 6 22 Cs 0.58451536649608 0.12500000000000 0.58465770183495 132.905 > 6 23 Cs 0.08451536649608 0.62500000000000 0.58465770183495 132.905 > 6 24 Cs 0.58451536649608 0.62500000000000 0.58465770183495 132.905 > 6 25 Cs 0.16548463350392 0.37500000000000 0.08465770183495 132.905 > 7 26 Cs 0.66548463350392 0.37500000000000 0.08465770183495 132.905 > 7 27 Cs 0.16548463350392 0.87500000000000 0.08465770183495 132.905 > 7 28 Cs 0.66548463350392 0.87500000000000 0.08465770183495 132.905 > 7 29 Cs 0.33451536649608 0.12500000000000 0.91534229816505 132.905 > 8 30 Cs 0.83451536649608 0.12500000000000 0.91534229816505 132.905 > 8 31 Cs 0.33451536649608 0.62500000000000 0.91534229816505 132.905 > 8 32 Cs 0.83451536649608 0.62500000000000 0.91534229816505 132.905 > 8 *33 Dy 0.40871348114862 0.37500000000000 0.06528498301047 162.500 > 9 34 Dy 0.90871348114862 0.37500000000000 0.06528498301047 162.500 > 9 35 Dy 0.40871348114862 0.87500000000000 0.06528498301047 162.500 > 9 36 Dy 0.90871348114862 0.87500000000000 0.06528498301047 162.500 > 9 37 Dy 0.09128651885138 0.12500000000000 0.93471501698953 162.500 > 10 38 Dy 0.59128651885138 0.12500000000000 0.93471501698953 162.500 > 10 39 Dy 0.09128651885138 0.62500000000000 0.93471501698953 162.500 > 10 40 Dy 0.59128651885138 0.62500000000000 0.93471501698953 162.500 > 10 41 Dy 0.15871348114862 0.37500000000000 0.43471501698953 162.500 > 11 42 Dy 0.65871348114862 0.37500000000000 0.43471501698953 162.500 > 11 43 Dy 0.15871348114862 0.87500000000000 0.43471501698953 162.500 > 11 44 Dy 0.65871348114862 0.87500000000000 0.43471501698953 162.500 > 11 45 Dy 0.34128651885138 0.12500000000000 0.56528498301047 162.500 > 12 46 Dy 0.84128651885138 0.12500000000000 0.56528498301047 162.500 > 12 47 Dy 0.34128651885138 0.62500000000000 0.56528498301047 162.500 > 12 48 Dy 0.84128651885138 0.62500000000000 0.56528498301047 162.500 > 12 *49 Cl 0.00000000000000 0.00000000000000 0.00000000000000 35.453 > 13 50 Cl 0.50000000000000 0.00000000000000 0.00000000000000 35.453 > 13 51 Cl 0.00000000000000 0.50000000000000 0.00000000000000 35.453 > 13 52 Cl 0.50000000000000 0.50000000000000 0.00000000000000 35.453 > 13 53 Cl 0.00000000000000 0.25000000000000 0.00000000000000 35.453 > 14 54 Cl 0.50000000000000 0.25000000000000 0.00000000000000 35.453 > 14 55 Cl 0.00000000000000 0.75000000000000 0.00000000000000 35.453 > 14 56 Cl 0.50000000000000 0.75000000000000 0.00000000000000 35.453 > 14 57 Cl 0.25000000000000 0.25000000000000 0.50000000000000 35.453 > 15 58 Cl 0.75000000000000 0.25000000000000 0.50000000000000 35.453 > 15 59 Cl 0.25000000000000 0.75000000000000 0.50000000000000 35.453 > 15 60 Cl 0.75000000000000 0.75000000000000 0.50000000000000 35.453 > 15 61 Cl 0.25000000000000 0.00000000000000 0.50000000000000 35.453 > 16 62 Cl 0.75000000000000 0.00000000000000 0.50000000000000 35.453 > 16 63 Cl 0.25000000000000 0.50000000000000 0.50000000000000 35.453 > 16 64 Cl 0.75000000000000 0.50000000000000 0.50000000000000 35.453 > 16 *65 Cl 0.32955137379588 0.24612884679132 0.11782660982414 35.453 > 17 66 Cl 0.82955137379588 0.24612884679132 0.11782660982414 35.453 > 17 67 Cl 0.32955137379588 0.74612884679132 0.11782660982414 35.453 > 17 68 Cl 0.82955137379588 0.74612884679132 0.11782660982414 35.453 > 17 69 Cl 0.17044862620412 0.49612884679132 0.88217339017586 35.453 > 18 70 Cl 0.67044862620412 0.49612884679132 0.88217339017586 35.453 > 18 71 Cl 0.17044862620412 0.99612884679132 0.88217339017586 35.453 > 18 72 Cl 0.67044862620412 0.99612884679132 0.88217339017586 35.453 > 18 73 Cl 0.07955137379588 0.00387115320868 0.38217339017586 35.453 > 19 74 Cl 0.57955137379588 0.00387115320868 0.38217339017586 35.453 > 19 75 Cl 0.07955137379588 0.50387115320868 0.38217339017586 35.453 > 19 76 Cl 0.57955137379588 0.50387115320868 0.38217339017586 35.453 > 19 77 Cl 0.42044862620412 0.25387115320868 0.61782660982414 35.453 > 20 78 Cl 0.92044862620412 0.25387115320868 0.61782660982414 35.453 > 20 79 Cl 0.42044862620412 0.75387115320868 0.61782660982414 35.453 > 20 80 Cl 0.92044862620412 0.75387115320868 0.61782660982414 35.453 > 20 81 Cl 0.17044862620412 0.25387115320868 0.88217339017586 35.453 > 21 82 Cl 0.67044862620412 0.25387115320868 0.88217339017586 35.453 > 21 83 Cl 0.17044862620412 0.75387115320868 0.88217339017586 35.453 > 21 84 Cl 0.67044862620412 0.75387115320868 0.88217339017586 35.453 > 21 85 Cl 0.32955137379588 0.00387115320868 0.11782660982414 35.453 > 22 86 Cl 0.82955137379588 0.00387115320868 0.11782660982414 35.453 > 22 87 Cl 0.32955137379588 0.50387115320868 0.11782660982414 35.453 > 22 88 Cl 0.82955137379588 0.50387115320868 0.11782660982414 35.453 > 22 89 Cl 0.42044862620412 0.49612884679132 0.61782660982414 35.453 > 23 90 Cl 0.92044862620412 0.49612884679132 0.61782660982414 35.453 > 23 91 Cl 0.42044862620412 0.99612884679132 0.61782660982414 35.453 > 23 92 Cl 0.92044862620412 0.99612884679132 0.61782660982414 35.453 > 23 93 Cl 0.07955137379588 0.24612884679132 0.38217339017586 35.453 > 24 94 Cl 0.57955137379588 0.24612884679132 0.38217339017586 35.453 > 24 95 Cl 0.07955137379588 0.74612884679132 0.38217339017586 35.453 > 24 96 Cl 0.57955137379588 0.74612884679132 0.38217339017586 35.453 > 24 *97 Cl 0.39423742496161 0.12500000000000 0.40531260520702 35.453 > 25 98 Cl 0.89423742496161 0.12500000000000 0.40531260520702 35.453 > 25 99 Cl 0.39423742496161 0.62500000000000 0.40531260520702 35.453 > 25 100 Cl 0.89423742496161 0.62500000000000 0.40531260520702 35.453 > 25 101 Cl 0.10576257503839 0.37500000000000 0.59468739479298 35.453 > 26 102 Cl 0.60576257503839 0.37500000000000 0.59468739479298 35.453 > 26 103 Cl 0.10576257503839 0.87500000000000 0.59468739479298 35.453 > 26 104 Cl 0.60576257503839 0.87500000000000 0.59468739479298 35.453 > 26 105 Cl 0.14423742496161 0.12500000000000 0.09468739479298 35.453 > 27 106 Cl 0.64423742496161 0.12500000000000 0.09468739479298 35.453 > 27 107 Cl 0.14423742496161 0.62500000000000 0.09468739479298 35.453 > 27 108 Cl 0.64423742496161 0.62500000000000 0.09468739479298 35.453 > 27 109 Cl 0.35576257503839 0.37500000000000 0.90531260520702 35.453 > 28 110 Cl 0.85576257503839 0.37500000000000 0.90531260520702 35.453 > 28 111 Cl 0.35576257503839 0.87500000000000 0.90531260520702 35.453 > 28 112 Cl 0.85576257503839 0.87500000000000 0.90531260520702 35.453 > 28 *113 Cl 0.48991308858784 0.37500000000000 0.20226779441939 35.453 > 29 114 Cl 0.98991308858784 0.37500000000000 0.20226779441939 35.453 > 29 115 Cl 0.48991308858784 0.87500000000000 0.20226779441939 35.453 > 29 116 Cl 0.98991308858784 0.87500000000000 0.20226779441939 35.453 > 29 117 Cl 0.01008691141216 0.12500000000000 0.79773220558061 35.453 > 30 118 Cl 0.51008691141216 0.12500000000000 0.79773220558061 35.453 > 30 119 Cl 0.01008691141216 0.62500000000000 0.79773220558061 35.453 > 30 120 Cl 0.51008691141216 0.62500000000000 0.79773220558061 35.453 > 30 121 Cl 0.23991308858784 0.37500000000000 0.29773220558061 35.453 > 31 122 Cl 0.73991308858784 0.37500000000000 0.29773220558061 35.453 > 31 123 Cl 0.23991308858784 0.87500000000000 0.29773220558061 35.453 > 31 124 Cl 0.73991308858784 0.87500000000000 0.29773220558061 35.453 > 31 125 Cl 0.26008691141216 0.12500000000000 0.70226779441939 35.453 > 32 126 Cl 0.76008691141216 0.12500000000000 0.70226779441939 35.453 > 32 127 Cl 0.26008691141216 0.62500000000000 0.70226779441939 35.453 > 32 128 Cl 0.76008691141216 0.62500000000000 0.70226779441939 35.453 > 32 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 2.8609084 0.0000000 0.0000000 0.0000000 2.9279013 0.0000000 0.0000000 0.0000000 2.8719783 -------------------------- Born effective charges -------------------------- 1 Cs 1.2061394 0.0000000 0.0228863 0.0000000 1.3701435 0.0000000 -0.0198765 0.0000000 1.5198233 2 Cs 1.2061394 0.0000000 0.0228863 0.0000000 1.3701435 0.0000000 -0.0198765 0.0000000 1.5198233 3 Cs 1.2061394 0.0000000 -0.0228863 0.0000000 1.3701435 0.0000000 0.0198765 0.0000000 1.5198233 4 Cs 1.2061394 0.0000000 -0.0228863 0.0000000 1.3701435 0.0000000 0.0198765 0.0000000 1.5198233 5 Cs 1.4024399 0.0000000 -0.1171344 0.0000000 1.2964748 0.0000000 -0.1117394 0.0000000 1.2870865 6 Cs 1.4024399 0.0000000 -0.1171344 0.0000000 1.2964748 0.0000000 -0.1117394 0.0000000 1.2870865 7 Cs 1.4024399 0.0000000 0.1171344 0.0000000 1.2964748 0.0000000 0.1117394 0.0000000 1.2870865 8 Cs 1.4024399 0.0000000 0.1171344 0.0000000 1.2964748 0.0000000 0.1117394 0.0000000 1.2870865 9 Dy 3.4787294 0.0000000 -0.0268019 0.0000000 3.7567839 0.0000000 -0.0207386 0.0000000 3.3159223 10 Dy 3.4787294 0.0000000 -0.0268019 0.0000000 3.7567839 0.0000000 -0.0207386 0.0000000 3.3159223 11 Dy 3.4787294 0.0000000 0.0268019 0.0000000 3.7567839 0.0000000 0.0207386 0.0000000 3.3159223 12 Dy 3.4787294 0.0000000 0.0268019 0.0000000 3.7567839 0.0000000 0.0207386 0.0000000 3.3159223 13 Cl -1.4202677 -0.7175615 0.4011783 -0.8319935 -1.6733687 0.4384569 0.4046592 0.3664735 -0.9603484 14 Cl -1.4202677 0.7175615 0.4011783 0.8319935 -1.6733687 -0.4384569 0.4046592 -0.3664735 -0.9603484 15 Cl -1.4202677 0.7175615 -0.4011783 0.8319935 -1.6733687 0.4384569 -0.4046592 0.3664735 -0.9603484 16 Cl -1.4202677 -0.7175615 -0.4011783 -0.8319935 -1.6733687 -0.4384569 -0.4046592 -0.3664735 -0.9603484 17 Cl -1.2569846 -0.4347298 0.2330427 -0.4134532 -1.4341400 0.2805469 0.2291280 0.2751855 -0.9278093 18 Cl -1.2569846 0.4347298 0.2330427 0.4134532 -1.4341400 -0.2805469 0.2291280 -0.2751855 -0.9278093 19 Cl -1.2569846 0.4347298 -0.2330427 0.4134532 -1.4341400 0.2805469 -0.2291280 0.2751855 -0.9278093 20 Cl -1.2569846 -0.4347298 -0.2330427 -0.4134532 -1.4341400 -0.2805469 -0.2291280 -0.2751855 -0.9278093 21 Cl -1.2569846 -0.4347298 0.2330427 -0.4134532 -1.4341400 0.2805469 0.2291280 0.2751855 -0.9278093 22 Cl -1.2569846 0.4347298 0.2330427 0.4134532 -1.4341400 -0.2805469 0.2291280 -0.2751855 -0.9278093 23 Cl -1.2569846 0.4347298 -0.2330427 0.4134532 -1.4341400 0.2805469 -0.2291280 0.2751855 -0.9278093 24 Cl -1.2569846 -0.4347298 -0.2330427 -0.4134532 -1.4341400 -0.2805469 -0.2291280 -0.2751855 -0.9278093 25 Cl -1.0553548 0.0000000 0.4681014 0.0000000 -0.9481055 0.0000000 0.2972869 0.0000000 -1.7296133 26 Cl -1.0553548 0.0000000 0.4681014 0.0000000 -0.9481055 0.0000000 0.2972869 0.0000000 -1.7296133 27 Cl -1.0553548 0.0000000 -0.4681014 0.0000000 -0.9481055 0.0000000 -0.2972869 0.0000000 -1.7296133 28 Cl -1.0553548 0.0000000 -0.4681014 0.0000000 -0.9481055 0.0000000 -0.2972869 0.0000000 -1.7296133 29 Cl -1.0977170 0.0000000 -0.3435436 0.0000000 -0.9336480 0.0000000 -0.4166316 0.0000000 -1.5772518 30 Cl -1.0977170 0.0000000 -0.3435436 0.0000000 -0.9336480 0.0000000 -0.4166316 0.0000000 -1.5772518 31 Cl -1.0977170 0.0000000 0.3435436 0.0000000 -0.9336480 0.0000000 0.4166316 0.0000000 -1.5772518 32 Cl -1.0977170 0.0000000 0.3435436 0.0000000 -0.9336480 0.0000000 0.4166316 0.0000000 -1.5772518 ---------------------------------------------------------------------------- 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: 384/384 Permutation basis: 18816/18816 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 2459 Number of blocks in projector: 2459 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 6 --- Eigsh_solver_block: 1 / 6 --- Block_size: 614 Use submatrix version of eigh solver. (1) Block: 500 / 614 (2) Block: 250 / 500 (2) 236 eigenvectors found. (2) Block: 500 / 500 (2) 236 eigenvectors found. (2) Complementary block size: 28 (2) Compute compressed projector. (2) Use standard solver. (2) 14 eigenvectors found. (1) 486 eigenvectors found. (1) Block: 614 / 614 (1) 104 eigenvectors found. (1) Complementary block size: 24 (1) Compute compressed projector. (1) Use standard solver. (1) 10 eigenvectors found. --- Eigsh_solver_block: 2 / 6 --- Block_size: 414 Use standard eigh solver. --- Eigsh_solver_block: 3 / 6 --- Block_size: 414 Use standard eigh solver. --- Eigsh_solver_block: 4 / 6 --- Block_size: 339 Use standard eigh solver. --- Eigsh_solver_block: 5 / 6 --- Block_size: 339 Use standard eigh solver. --- Eigsh_solver_block: 6 / 6 --- Block_size: 339 Use standard eigh solver. Tree of FC basis block matrices: - (2459, 2416), data: False |-- (339, 332), data: True |-- (339, 332), data: True |-- (339, 332), data: True |-- (414, 410), data: True |-- (414, 410), data: True |-- (614, 600), data: False |-- (614, 10), data: (614, 24) @ (24, 10) |-- (114, 104), data: True |-- (500, 486), data: False |-- (500, 14), data: (500, 28) @ (28, 14) |-- (250, 236), data: True |-- (250, 236), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.003 Solver_block: 100 / 480 - Time: 0.881 Solver_block: 200 / 480 - Time: 1.302 Solver_block: 300 / 480 - Time: 0.645 Solver_block: 400 / 480 - Time: 0.651 Solver_block: 480 / 480 - Time: 0.625 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 4.472 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 384/384 Permutation basis: 18816/18816 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 2459 Number of blocks in projector: 2459 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 6 --- Eigsh_solver_block: 1 / 6 --- Block_size: 614 Use submatrix version of eigh solver. (1) Block: 500 / 614 (2) Block: 250 / 500 (2) 236 eigenvectors found. (2) Block: 500 / 500 (2) 236 eigenvectors found. (2) Complementary block size: 28 (2) Compute compressed projector. (2) Use standard solver. (2) 14 eigenvectors found. (1) 486 eigenvectors found. (1) Block: 614 / 614 (1) 104 eigenvectors found. (1) Complementary block size: 24 (1) Compute compressed projector. (1) Use standard solver. (1) 10 eigenvectors found. --- Eigsh_solver_block: 2 / 6 --- Block_size: 414 Use standard eigh solver. --- Eigsh_solver_block: 3 / 6 --- Block_size: 414 Use standard eigh solver. --- Eigsh_solver_block: 4 / 6 --- Block_size: 339 Use standard eigh solver. --- Eigsh_solver_block: 5 / 6 --- Block_size: 339 Use standard eigh solver. --- Eigsh_solver_block: 6 / 6 --- Block_size: 339 Use standard eigh solver. Tree of FC basis block matrices: - (2459, 2416), data: False |-- (339, 332), data: True |-- (339, 332), data: True |-- (339, 332), data: True |-- (414, 410), data: True |-- (414, 410), data: True |-- (614, 600), data: False |-- (614, 10), data: (614, 24) @ (24, 10) |-- (114, 104), data: True |-- (500, 486), data: False |-- (500, 14), data: (500, 28) @ (28, 14) |-- (250, 236), data: True |-- (250, 236), 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 09:58:36]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 09:58:37]------------------------- 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): [2 2 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pnma (62) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 9.434285200000000 0.000000000000000 0.000000000000000 b 0.000000000000000 7.381102380000000 0.000000000000000 c 0.000000000000000 0.000000000000000 15.111277439999999 Atomic positions (fractional): 1 Cs 0.96992290148160 0.25000000000000 0.20911211981328 132.905 2 Cs 0.03007709851840 0.75000000000000 0.79088788018672 132.905 3 Cs 0.46992290148160 0.25000000000000 0.29088788018672 132.905 4 Cs 0.53007709851840 0.75000000000000 0.70911211981328 132.905 5 Cs 0.83096926700783 0.75000000000000 0.41534229816506 132.905 6 Cs 0.16903073299217 0.25000000000000 0.58465770183495 132.905 7 Cs 0.33096926700783 0.75000000000000 0.08465770183495 132.905 8 Cs 0.66903073299217 0.25000000000000 0.91534229816505 132.905 9 Dy 0.81742696229724 0.75000000000000 0.06528498301047 162.500 10 Dy 0.18257303770276 0.25000000000000 0.93471501698953 162.500 11 Dy 0.31742696229724 0.75000000000000 0.43471501698953 162.500 12 Dy 0.68257303770276 0.25000000000000 0.56528498301047 162.500 13 Cl 0.00000000000000 0.00000000000000 0.00000000000000 35.453 14 Cl 0.00000000000000 0.50000000000000 0.00000000000000 35.453 15 Cl 0.50000000000000 0.50000000000000 0.50000000000000 35.453 16 Cl 0.50000000000000 0.00000000000000 0.50000000000000 35.453 17 Cl 0.65910274759175 0.49225769358265 0.11782660982414 35.453 18 Cl 0.34089725240825 0.99225769358265 0.88217339017586 35.453 19 Cl 0.15910274759175 0.00774230641735 0.38217339017586 35.453 20 Cl 0.84089725240825 0.50774230641735 0.61782660982414 35.453 21 Cl 0.34089725240825 0.50774230641735 0.88217339017586 35.453 22 Cl 0.65910274759175 0.00774230641735 0.11782660982414 35.453 23 Cl 0.84089725240825 0.99225769358265 0.61782660982414 35.453 24 Cl 0.15910274759175 0.49225769358265 0.38217339017586 35.453 25 Cl 0.78847484992322 0.25000000000000 0.40531260520702 35.453 26 Cl 0.21152515007678 0.75000000000000 0.59468739479298 35.453 27 Cl 0.28847484992322 0.25000000000000 0.09468739479298 35.453 28 Cl 0.71152515007678 0.75000000000000 0.90531260520702 35.453 29 Cl 0.97982617717568 0.75000000000000 0.20226779441939 35.453 30 Cl 0.02017382282432 0.25000000000000 0.79773220558061 35.453 31 Cl 0.47982617717568 0.75000000000000 0.29773220558061 35.453 32 Cl 0.52017382282432 0.25000000000000 0.70226779441939 35.453 -------------------------------- supercell --------------------------------- Lattice vectors: a 18.868570399999999 0.000000000000000 0.000000000000000 b 0.000000000000000 14.762204759999999 0.000000000000000 c 0.000000000000000 0.000000000000000 15.111277439999999 Atomic positions (fractional): 1 Cs 0.48496145074080 0.12500000000000 0.20911211981328 132.905 > 1 2 Cs 0.98496145074080 0.12500000000000 0.20911211981328 132.905 > 1 3 Cs 0.48496145074080 0.62500000000000 0.20911211981328 132.905 > 1 4 Cs 0.98496145074080 0.62500000000000 0.20911211981328 132.905 > 1 5 Cs 0.01503854925920 0.37500000000000 0.79088788018672 132.905 > 5 6 Cs 0.51503854925920 0.37500000000000 0.79088788018672 132.905 > 5 7 Cs 0.01503854925920 0.87500000000000 0.79088788018672 132.905 > 5 8 Cs 0.51503854925920 0.87500000000000 0.79088788018672 132.905 > 5 9 Cs 0.23496145074080 0.12500000000000 0.29088788018672 132.905 > 9 10 Cs 0.73496145074080 0.12500000000000 0.29088788018672 132.905 > 9 11 Cs 0.23496145074080 0.62500000000000 0.29088788018672 132.905 > 9 12 Cs 0.73496145074080 0.62500000000000 0.29088788018672 132.905 > 9 13 Cs 0.26503854925920 0.37500000000000 0.70911211981328 132.905 > 13 14 Cs 0.76503854925920 0.37500000000000 0.70911211981328 132.905 > 13 15 Cs 0.26503854925920 0.87500000000000 0.70911211981328 132.905 > 13 16 Cs 0.76503854925920 0.87500000000000 0.70911211981328 132.905 > 13 17 Cs 0.41548463350392 0.37500000000000 0.41534229816506 132.905 > 17 18 Cs 0.91548463350392 0.37500000000000 0.41534229816506 132.905 > 17 19 Cs 0.41548463350392 0.87500000000000 0.41534229816506 132.905 > 17 20 Cs 0.91548463350392 0.87500000000000 0.41534229816506 132.905 > 17 21 Cs 0.08451536649608 0.12500000000000 0.58465770183495 132.905 > 21 22 Cs 0.58451536649608 0.12500000000000 0.58465770183495 132.905 > 21 23 Cs 0.08451536649608 0.62500000000000 0.58465770183495 132.905 > 21 24 Cs 0.58451536649608 0.62500000000000 0.58465770183495 132.905 > 21 25 Cs 0.16548463350392 0.37500000000000 0.08465770183495 132.905 > 25 26 Cs 0.66548463350392 0.37500000000000 0.08465770183495 132.905 > 25 27 Cs 0.16548463350392 0.87500000000000 0.08465770183495 132.905 > 25 28 Cs 0.66548463350392 0.87500000000000 0.08465770183495 132.905 > 25 29 Cs 0.33451536649608 0.12500000000000 0.91534229816505 132.905 > 29 30 Cs 0.83451536649608 0.12500000000000 0.91534229816505 132.905 > 29 31 Cs 0.33451536649608 0.62500000000000 0.91534229816505 132.905 > 29 32 Cs 0.83451536649608 0.62500000000000 0.91534229816505 132.905 > 29 33 Dy 0.40871348114862 0.37500000000000 0.06528498301047 162.500 > 33 34 Dy 0.90871348114862 0.37500000000000 0.06528498301047 162.500 > 33 35 Dy 0.40871348114862 0.87500000000000 0.06528498301047 162.500 > 33 36 Dy 0.90871348114862 0.87500000000000 0.06528498301047 162.500 > 33 37 Dy 0.09128651885138 0.12500000000000 0.93471501698953 162.500 > 37 38 Dy 0.59128651885138 0.12500000000000 0.93471501698953 162.500 > 37 39 Dy 0.09128651885138 0.62500000000000 0.93471501698953 162.500 > 37 40 Dy 0.59128651885138 0.62500000000000 0.93471501698953 162.500 > 37 41 Dy 0.15871348114862 0.37500000000000 0.43471501698953 162.500 > 41 42 Dy 0.65871348114862 0.37500000000000 0.43471501698953 162.500 > 41 43 Dy 0.15871348114862 0.87500000000000 0.43471501698953 162.500 > 41 44 Dy 0.65871348114862 0.87500000000000 0.43471501698953 162.500 > 41 45 Dy 0.34128651885138 0.12500000000000 0.56528498301047 162.500 > 45 46 Dy 0.84128651885138 0.12500000000000 0.56528498301047 162.500 > 45 47 Dy 0.34128651885138 0.62500000000000 0.56528498301047 162.500 > 45 48 Dy 0.84128651885138 0.62500000000000 0.56528498301047 162.500 > 45 49 Cl 0.00000000000000 0.00000000000000 0.00000000000000 35.453 > 49 50 Cl 0.50000000000000 0.00000000000000 0.00000000000000 35.453 > 49 51 Cl 0.00000000000000 0.50000000000000 0.00000000000000 35.453 > 49 52 Cl 0.50000000000000 0.50000000000000 0.00000000000000 35.453 > 49 53 Cl 0.00000000000000 0.25000000000000 0.00000000000000 35.453 > 53 54 Cl 0.50000000000000 0.25000000000000 0.00000000000000 35.453 > 53 55 Cl 0.00000000000000 0.75000000000000 0.00000000000000 35.453 > 53 56 Cl 0.50000000000000 0.75000000000000 0.00000000000000 35.453 > 53 57 Cl 0.25000000000000 0.25000000000000 0.50000000000000 35.453 > 57 58 Cl 0.75000000000000 0.25000000000000 0.50000000000000 35.453 > 57 59 Cl 0.25000000000000 0.75000000000000 0.50000000000000 35.453 > 57 60 Cl 0.75000000000000 0.75000000000000 0.50000000000000 35.453 > 57 61 Cl 0.25000000000000 0.00000000000000 0.50000000000000 35.453 > 61 62 Cl 0.75000000000000 0.00000000000000 0.50000000000000 35.453 > 61 63 Cl 0.25000000000000 0.50000000000000 0.50000000000000 35.453 > 61 64 Cl 0.75000000000000 0.50000000000000 0.50000000000000 35.453 > 61 65 Cl 0.32955137379588 0.24612884679132 0.11782660982414 35.453 > 65 66 Cl 0.82955137379588 0.24612884679132 0.11782660982414 35.453 > 65 67 Cl 0.32955137379588 0.74612884679132 0.11782660982414 35.453 > 65 68 Cl 0.82955137379588 0.74612884679132 0.11782660982414 35.453 > 65 69 Cl 0.17044862620412 0.49612884679132 0.88217339017586 35.453 > 69 70 Cl 0.67044862620412 0.49612884679132 0.88217339017586 35.453 > 69 71 Cl 0.17044862620412 0.99612884679132 0.88217339017586 35.453 > 69 72 Cl 0.67044862620412 0.99612884679132 0.88217339017586 35.453 > 69 73 Cl 0.07955137379588 0.00387115320868 0.38217339017586 35.453 > 73 74 Cl 0.57955137379588 0.00387115320868 0.38217339017586 35.453 > 73 75 Cl 0.07955137379588 0.50387115320868 0.38217339017586 35.453 > 73 76 Cl 0.57955137379588 0.50387115320868 0.38217339017586 35.453 > 73 77 Cl 0.42044862620412 0.25387115320868 0.61782660982414 35.453 > 77 78 Cl 0.92044862620412 0.25387115320868 0.61782660982414 35.453 > 77 79 Cl 0.42044862620412 0.75387115320868 0.61782660982414 35.453 > 77 80 Cl 0.92044862620412 0.75387115320868 0.61782660982414 35.453 > 77 81 Cl 0.17044862620412 0.25387115320868 0.88217339017586 35.453 > 81 82 Cl 0.67044862620412 0.25387115320868 0.88217339017586 35.453 > 81 83 Cl 0.17044862620412 0.75387115320868 0.88217339017586 35.453 > 81 84 Cl 0.67044862620412 0.75387115320868 0.88217339017586 35.453 > 81 85 Cl 0.32955137379588 0.00387115320868 0.11782660982414 35.453 > 85 86 Cl 0.82955137379588 0.00387115320868 0.11782660982414 35.453 > 85 87 Cl 0.32955137379588 0.50387115320868 0.11782660982414 35.453 > 85 88 Cl 0.82955137379588 0.50387115320868 0.11782660982414 35.453 > 85 89 Cl 0.42044862620412 0.49612884679132 0.61782660982414 35.453 > 89 90 Cl 0.92044862620412 0.49612884679132 0.61782660982414 35.453 > 89 91 Cl 0.42044862620412 0.99612884679132 0.61782660982414 35.453 > 89 92 Cl 0.92044862620412 0.99612884679132 0.61782660982414 35.453 > 89 93 Cl 0.07955137379588 0.24612884679132 0.38217339017586 35.453 > 93 94 Cl 0.57955137379588 0.24612884679132 0.38217339017586 35.453 > 93 95 Cl 0.07955137379588 0.74612884679132 0.38217339017586 35.453 > 93 96 Cl 0.57955137379588 0.74612884679132 0.38217339017586 35.453 > 93 97 Cl 0.39423742496161 0.12500000000000 0.40531260520702 35.453 > 97 98 Cl 0.89423742496161 0.12500000000000 0.40531260520702 35.453 > 97 99 Cl 0.39423742496161 0.62500000000000 0.40531260520702 35.453 > 97 100 Cl 0.89423742496161 0.62500000000000 0.40531260520702 35.453 > 97 101 Cl 0.10576257503839 0.37500000000000 0.59468739479298 35.453 > 101 102 Cl 0.60576257503839 0.37500000000000 0.59468739479298 35.453 > 101 103 Cl 0.10576257503839 0.87500000000000 0.59468739479298 35.453 > 101 104 Cl 0.60576257503839 0.87500000000000 0.59468739479298 35.453 > 101 105 Cl 0.14423742496161 0.12500000000000 0.09468739479298 35.453 > 105 106 Cl 0.64423742496161 0.12500000000000 0.09468739479298 35.453 > 105 107 Cl 0.14423742496161 0.62500000000000 0.09468739479298 35.453 > 105 108 Cl 0.64423742496161 0.62500000000000 0.09468739479298 35.453 > 105 109 Cl 0.35576257503839 0.37500000000000 0.90531260520702 35.453 > 109 110 Cl 0.85576257503839 0.37500000000000 0.90531260520702 35.453 > 109 111 Cl 0.35576257503839 0.87500000000000 0.90531260520702 35.453 > 109 112 Cl 0.85576257503839 0.87500000000000 0.90531260520702 35.453 > 109 113 Cl 0.48991308858784 0.37500000000000 0.20226779441939 35.453 > 113 114 Cl 0.98991308858784 0.37500000000000 0.20226779441939 35.453 > 113 115 Cl 0.48991308858784 0.87500000000000 0.20226779441939 35.453 > 113 116 Cl 0.98991308858784 0.87500000000000 0.20226779441939 35.453 > 113 117 Cl 0.01008691141216 0.12500000000000 0.79773220558061 35.453 > 117 118 Cl 0.51008691141216 0.12500000000000 0.79773220558061 35.453 > 117 119 Cl 0.01008691141216 0.62500000000000 0.79773220558061 35.453 > 117 120 Cl 0.51008691141216 0.62500000000000 0.79773220558061 35.453 > 117 121 Cl 0.23991308858784 0.37500000000000 0.29773220558061 35.453 > 121 122 Cl 0.73991308858784 0.37500000000000 0.29773220558061 35.453 > 121 123 Cl 0.23991308858784 0.87500000000000 0.29773220558061 35.453 > 121 124 Cl 0.73991308858784 0.87500000000000 0.29773220558061 35.453 > 121 125 Cl 0.26008691141216 0.12500000000000 0.70226779441939 35.453 > 125 126 Cl 0.76008691141216 0.12500000000000 0.70226779441939 35.453 > 125 127 Cl 0.26008691141216 0.62500000000000 0.70226779441939 35.453 > 125 128 Cl 0.76008691141216 0.62500000000000 0.70226779441939 35.453 > 125 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 2.8609084 0.0000000 0.0000000 0.0000000 2.9279013 0.0000000 0.0000000 0.0000000 2.8719783 -------------------------- Born effective charges -------------------------- 1 Cs 1.2061394 0.0000000 0.0228863 0.0000000 1.3701435 0.0000000 -0.0198765 0.0000000 1.5198233 2 Cs 1.2061394 0.0000000 0.0228863 0.0000000 1.3701435 0.0000000 -0.0198765 0.0000000 1.5198233 3 Cs 1.2061394 0.0000000 -0.0228863 0.0000000 1.3701435 0.0000000 0.0198765 0.0000000 1.5198233 4 Cs 1.2061394 0.0000000 -0.0228863 0.0000000 1.3701435 0.0000000 0.0198765 0.0000000 1.5198233 5 Cs 1.4024399 0.0000000 -0.1171344 0.0000000 1.2964748 0.0000000 -0.1117394 0.0000000 1.2870865 6 Cs 1.4024399 0.0000000 -0.1171344 0.0000000 1.2964748 0.0000000 -0.1117394 0.0000000 1.2870865 7 Cs 1.4024399 0.0000000 0.1171344 0.0000000 1.2964748 0.0000000 0.1117394 0.0000000 1.2870865 8 Cs 1.4024399 0.0000000 0.1171344 0.0000000 1.2964748 0.0000000 0.1117394 0.0000000 1.2870865 9 Dy 3.4787294 0.0000000 -0.0268019 0.0000000 3.7567839 0.0000000 -0.0207386 0.0000000 3.3159223 10 Dy 3.4787294 0.0000000 -0.0268019 0.0000000 3.7567839 0.0000000 -0.0207386 0.0000000 3.3159223 11 Dy 3.4787294 0.0000000 0.0268019 0.0000000 3.7567839 0.0000000 0.0207386 0.0000000 3.3159223 12 Dy 3.4787294 0.0000000 0.0268019 0.0000000 3.7567839 0.0000000 0.0207386 0.0000000 3.3159223 13 Cl -1.4202677 -0.7175615 0.4011783 -0.8319935 -1.6733687 0.4384569 0.4046592 0.3664735 -0.9603484 14 Cl -1.4202677 0.7175615 0.4011783 0.8319935 -1.6733687 -0.4384569 0.4046592 -0.3664735 -0.9603484 15 Cl -1.4202677 0.7175615 -0.4011783 0.8319935 -1.6733687 0.4384569 -0.4046592 0.3664735 -0.9603484 16 Cl -1.4202677 -0.7175615 -0.4011783 -0.8319935 -1.6733687 -0.4384569 -0.4046592 -0.3664735 -0.9603484 17 Cl -1.2569846 -0.4347298 0.2330427 -0.4134532 -1.4341400 0.2805469 0.2291280 0.2751855 -0.9278093 18 Cl -1.2569846 0.4347298 0.2330427 0.4134532 -1.4341400 -0.2805469 0.2291280 -0.2751855 -0.9278093 19 Cl -1.2569846 0.4347298 -0.2330427 0.4134532 -1.4341400 0.2805469 -0.2291280 0.2751855 -0.9278093 20 Cl -1.2569846 -0.4347298 -0.2330427 -0.4134532 -1.4341400 -0.2805469 -0.2291280 -0.2751855 -0.9278093 21 Cl -1.2569846 -0.4347298 0.2330427 -0.4134532 -1.4341400 0.2805469 0.2291280 0.2751855 -0.9278093 22 Cl -1.2569846 0.4347298 0.2330427 0.4134532 -1.4341400 -0.2805469 0.2291280 -0.2751855 -0.9278093 23 Cl -1.2569846 0.4347298 -0.2330427 0.4134532 -1.4341400 0.2805469 -0.2291280 0.2751855 -0.9278093 24 Cl -1.2569846 -0.4347298 -0.2330427 -0.4134532 -1.4341400 -0.2805469 -0.2291280 -0.2751855 -0.9278093 25 Cl -1.0553548 0.0000000 0.4681014 0.0000000 -0.9481055 0.0000000 0.2972869 0.0000000 -1.7296133 26 Cl -1.0553548 0.0000000 0.4681014 0.0000000 -0.9481055 0.0000000 0.2972869 0.0000000 -1.7296133 27 Cl -1.0553548 0.0000000 -0.4681014 0.0000000 -0.9481055 0.0000000 -0.2972869 0.0000000 -1.7296133 28 Cl -1.0553548 0.0000000 -0.4681014 0.0000000 -0.9481055 0.0000000 -0.2972869 0.0000000 -1.7296133 29 Cl -1.0977170 0.0000000 -0.3435436 0.0000000 -0.9336480 0.0000000 -0.4166316 0.0000000 -1.5772518 30 Cl -1.0977170 0.0000000 -0.3435436 0.0000000 -0.9336480 0.0000000 -0.4166316 0.0000000 -1.5772518 31 Cl -1.0977170 0.0000000 0.3435436 0.0000000 -0.9336480 0.0000000 0.4166316 0.0000000 -1.5772518 32 Cl -1.0977170 0.0000000 0.3435436 0.0000000 -0.9336480 0.0000000 0.4166316 0.0000000 -1.5772518 ---------------------------------------------------------------------------- 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[ 17, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 49, 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[ 65, 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[ 97, 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[ 113, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000001 (xzx) 0.00000001 (xzx) 0.00000001 (xxz) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (xz) 0.00000000 (xz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 09:59:10]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 09:59:11]------------------------- 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): [2 2 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pnma (62) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 9.434285200000000 0.000000000000000 0.000000000000000 b 0.000000000000000 7.381102380000000 0.000000000000000 c 0.000000000000000 0.000000000000000 15.111277439999999 Atomic positions (fractional): 1 Cs 0.96992290148160 0.25000000000000 0.20911211981328 132.905 2 Cs 0.03007709851840 0.75000000000000 0.79088788018672 132.905 3 Cs 0.46992290148160 0.25000000000000 0.29088788018672 132.905 4 Cs 0.53007709851840 0.75000000000000 0.70911211981328 132.905 5 Cs 0.83096926700783 0.75000000000000 0.41534229816506 132.905 6 Cs 0.16903073299217 0.25000000000000 0.58465770183495 132.905 7 Cs 0.33096926700783 0.75000000000000 0.08465770183495 132.905 8 Cs 0.66903073299217 0.25000000000000 0.91534229816505 132.905 9 Dy 0.81742696229724 0.75000000000000 0.06528498301047 162.500 10 Dy 0.18257303770276 0.25000000000000 0.93471501698953 162.500 11 Dy 0.31742696229724 0.75000000000000 0.43471501698953 162.500 12 Dy 0.68257303770276 0.25000000000000 0.56528498301047 162.500 13 Cl 0.00000000000000 0.00000000000000 0.00000000000000 35.453 14 Cl 0.00000000000000 0.50000000000000 0.00000000000000 35.453 15 Cl 0.50000000000000 0.50000000000000 0.50000000000000 35.453 16 Cl 0.50000000000000 0.00000000000000 0.50000000000000 35.453 17 Cl 0.65910274759175 0.49225769358265 0.11782660982414 35.453 18 Cl 0.34089725240825 0.99225769358265 0.88217339017586 35.453 19 Cl 0.15910274759175 0.00774230641735 0.38217339017586 35.453 20 Cl 0.84089725240825 0.50774230641735 0.61782660982414 35.453 21 Cl 0.34089725240825 0.50774230641735 0.88217339017586 35.453 22 Cl 0.65910274759175 0.00774230641735 0.11782660982414 35.453 23 Cl 0.84089725240825 0.99225769358265 0.61782660982414 35.453 24 Cl 0.15910274759175 0.49225769358265 0.38217339017586 35.453 25 Cl 0.78847484992322 0.25000000000000 0.40531260520702 35.453 26 Cl 0.21152515007678 0.75000000000000 0.59468739479298 35.453 27 Cl 0.28847484992322 0.25000000000000 0.09468739479298 35.453 28 Cl 0.71152515007678 0.75000000000000 0.90531260520702 35.453 29 Cl 0.97982617717568 0.75000000000000 0.20226779441939 35.453 30 Cl 0.02017382282432 0.25000000000000 0.79773220558061 35.453 31 Cl 0.47982617717568 0.75000000000000 0.29773220558061 35.453 32 Cl 0.52017382282432 0.25000000000000 0.70226779441939 35.453 -------------------------------- supercell --------------------------------- Lattice vectors: a 18.868570399999999 0.000000000000000 0.000000000000000 b 0.000000000000000 14.762204759999999 0.000000000000000 c 0.000000000000000 0.000000000000000 15.111277439999999 Atomic positions (fractional): 1 Cs 0.48496145074080 0.12500000000000 0.20911211981328 132.905 > 1 2 Cs 0.98496145074080 0.12500000000000 0.20911211981328 132.905 > 1 3 Cs 0.48496145074080 0.62500000000000 0.20911211981328 132.905 > 1 4 Cs 0.98496145074080 0.62500000000000 0.20911211981328 132.905 > 1 5 Cs 0.01503854925920 0.37500000000000 0.79088788018672 132.905 > 5 6 Cs 0.51503854925920 0.37500000000000 0.79088788018672 132.905 > 5 7 Cs 0.01503854925920 0.87500000000000 0.79088788018672 132.905 > 5 8 Cs 0.51503854925920 0.87500000000000 0.79088788018672 132.905 > 5 9 Cs 0.23496145074080 0.12500000000000 0.29088788018672 132.905 > 9 10 Cs 0.73496145074080 0.12500000000000 0.29088788018672 132.905 > 9 11 Cs 0.23496145074080 0.62500000000000 0.29088788018672 132.905 > 9 12 Cs 0.73496145074080 0.62500000000000 0.29088788018672 132.905 > 9 13 Cs 0.26503854925920 0.37500000000000 0.70911211981328 132.905 > 13 14 Cs 0.76503854925920 0.37500000000000 0.70911211981328 132.905 > 13 15 Cs 0.26503854925920 0.87500000000000 0.70911211981328 132.905 > 13 16 Cs 0.76503854925920 0.87500000000000 0.70911211981328 132.905 > 13 17 Cs 0.41548463350392 0.37500000000000 0.41534229816506 132.905 > 17 18 Cs 0.91548463350392 0.37500000000000 0.41534229816506 132.905 > 17 19 Cs 0.41548463350392 0.87500000000000 0.41534229816506 132.905 > 17 20 Cs 0.91548463350392 0.87500000000000 0.41534229816506 132.905 > 17 21 Cs 0.08451536649608 0.12500000000000 0.58465770183495 132.905 > 21 22 Cs 0.58451536649608 0.12500000000000 0.58465770183495 132.905 > 21 23 Cs 0.08451536649608 0.62500000000000 0.58465770183495 132.905 > 21 24 Cs 0.58451536649608 0.62500000000000 0.58465770183495 132.905 > 21 25 Cs 0.16548463350392 0.37500000000000 0.08465770183495 132.905 > 25 26 Cs 0.66548463350392 0.37500000000000 0.08465770183495 132.905 > 25 27 Cs 0.16548463350392 0.87500000000000 0.08465770183495 132.905 > 25 28 Cs 0.66548463350392 0.87500000000000 0.08465770183495 132.905 > 25 29 Cs 0.33451536649608 0.12500000000000 0.91534229816505 132.905 > 29 30 Cs 0.83451536649608 0.12500000000000 0.91534229816505 132.905 > 29 31 Cs 0.33451536649608 0.62500000000000 0.91534229816505 132.905 > 29 32 Cs 0.83451536649608 0.62500000000000 0.91534229816505 132.905 > 29 33 Dy 0.40871348114862 0.37500000000000 0.06528498301047 162.500 > 33 34 Dy 0.90871348114862 0.37500000000000 0.06528498301047 162.500 > 33 35 Dy 0.40871348114862 0.87500000000000 0.06528498301047 162.500 > 33 36 Dy 0.90871348114862 0.87500000000000 0.06528498301047 162.500 > 33 37 Dy 0.09128651885138 0.12500000000000 0.93471501698953 162.500 > 37 38 Dy 0.59128651885138 0.12500000000000 0.93471501698953 162.500 > 37 39 Dy 0.09128651885138 0.62500000000000 0.93471501698953 162.500 > 37 40 Dy 0.59128651885138 0.62500000000000 0.93471501698953 162.500 > 37 41 Dy 0.15871348114862 0.37500000000000 0.43471501698953 162.500 > 41 42 Dy 0.65871348114862 0.37500000000000 0.43471501698953 162.500 > 41 43 Dy 0.15871348114862 0.87500000000000 0.43471501698953 162.500 > 41 44 Dy 0.65871348114862 0.87500000000000 0.43471501698953 162.500 > 41 45 Dy 0.34128651885138 0.12500000000000 0.56528498301047 162.500 > 45 46 Dy 0.84128651885138 0.12500000000000 0.56528498301047 162.500 > 45 47 Dy 0.34128651885138 0.62500000000000 0.56528498301047 162.500 > 45 48 Dy 0.84128651885138 0.62500000000000 0.56528498301047 162.500 > 45 49 Cl 0.00000000000000 0.00000000000000 0.00000000000000 35.453 > 49 50 Cl 0.50000000000000 0.00000000000000 0.00000000000000 35.453 > 49 51 Cl 0.00000000000000 0.50000000000000 0.00000000000000 35.453 > 49 52 Cl 0.50000000000000 0.50000000000000 0.00000000000000 35.453 > 49 53 Cl 0.00000000000000 0.25000000000000 0.00000000000000 35.453 > 53 54 Cl 0.50000000000000 0.25000000000000 0.00000000000000 35.453 > 53 55 Cl 0.00000000000000 0.75000000000000 0.00000000000000 35.453 > 53 56 Cl 0.50000000000000 0.75000000000000 0.00000000000000 35.453 > 53 57 Cl 0.25000000000000 0.25000000000000 0.50000000000000 35.453 > 57 58 Cl 0.75000000000000 0.25000000000000 0.50000000000000 35.453 > 57 59 Cl 0.25000000000000 0.75000000000000 0.50000000000000 35.453 > 57 60 Cl 0.75000000000000 0.75000000000000 0.50000000000000 35.453 > 57 61 Cl 0.25000000000000 0.00000000000000 0.50000000000000 35.453 > 61 62 Cl 0.75000000000000 0.00000000000000 0.50000000000000 35.453 > 61 63 Cl 0.25000000000000 0.50000000000000 0.50000000000000 35.453 > 61 64 Cl 0.75000000000000 0.50000000000000 0.50000000000000 35.453 > 61 65 Cl 0.32955137379588 0.24612884679132 0.11782660982414 35.453 > 65 66 Cl 0.82955137379588 0.24612884679132 0.11782660982414 35.453 > 65 67 Cl 0.32955137379588 0.74612884679132 0.11782660982414 35.453 > 65 68 Cl 0.82955137379588 0.74612884679132 0.11782660982414 35.453 > 65 69 Cl 0.17044862620412 0.49612884679132 0.88217339017586 35.453 > 69 70 Cl 0.67044862620412 0.49612884679132 0.88217339017586 35.453 > 69 71 Cl 0.17044862620412 0.99612884679132 0.88217339017586 35.453 > 69 72 Cl 0.67044862620412 0.99612884679132 0.88217339017586 35.453 > 69 73 Cl 0.07955137379588 0.00387115320868 0.38217339017586 35.453 > 73 74 Cl 0.57955137379588 0.00387115320868 0.38217339017586 35.453 > 73 75 Cl 0.07955137379588 0.50387115320868 0.38217339017586 35.453 > 73 76 Cl 0.57955137379588 0.50387115320868 0.38217339017586 35.453 > 73 77 Cl 0.42044862620412 0.25387115320868 0.61782660982414 35.453 > 77 78 Cl 0.92044862620412 0.25387115320868 0.61782660982414 35.453 > 77 79 Cl 0.42044862620412 0.75387115320868 0.61782660982414 35.453 > 77 80 Cl 0.92044862620412 0.75387115320868 0.61782660982414 35.453 > 77 81 Cl 0.17044862620412 0.25387115320868 0.88217339017586 35.453 > 81 82 Cl 0.67044862620412 0.25387115320868 0.88217339017586 35.453 > 81 83 Cl 0.17044862620412 0.75387115320868 0.88217339017586 35.453 > 81 84 Cl 0.67044862620412 0.75387115320868 0.88217339017586 35.453 > 81 85 Cl 0.32955137379588 0.00387115320868 0.11782660982414 35.453 > 85 86 Cl 0.82955137379588 0.00387115320868 0.11782660982414 35.453 > 85 87 Cl 0.32955137379588 0.50387115320868 0.11782660982414 35.453 > 85 88 Cl 0.82955137379588 0.50387115320868 0.11782660982414 35.453 > 85 89 Cl 0.42044862620412 0.49612884679132 0.61782660982414 35.453 > 89 90 Cl 0.92044862620412 0.49612884679132 0.61782660982414 35.453 > 89 91 Cl 0.42044862620412 0.99612884679132 0.61782660982414 35.453 > 89 92 Cl 0.92044862620412 0.99612884679132 0.61782660982414 35.453 > 89 93 Cl 0.07955137379588 0.24612884679132 0.38217339017586 35.453 > 93 94 Cl 0.57955137379588 0.24612884679132 0.38217339017586 35.453 > 93 95 Cl 0.07955137379588 0.74612884679132 0.38217339017586 35.453 > 93 96 Cl 0.57955137379588 0.74612884679132 0.38217339017586 35.453 > 93 97 Cl 0.39423742496161 0.12500000000000 0.40531260520702 35.453 > 97 98 Cl 0.89423742496161 0.12500000000000 0.40531260520702 35.453 > 97 99 Cl 0.39423742496161 0.62500000000000 0.40531260520702 35.453 > 97 100 Cl 0.89423742496161 0.62500000000000 0.40531260520702 35.453 > 97 101 Cl 0.10576257503839 0.37500000000000 0.59468739479298 35.453 > 101 102 Cl 0.60576257503839 0.37500000000000 0.59468739479298 35.453 > 101 103 Cl 0.10576257503839 0.87500000000000 0.59468739479298 35.453 > 101 104 Cl 0.60576257503839 0.87500000000000 0.59468739479298 35.453 > 101 105 Cl 0.14423742496161 0.12500000000000 0.09468739479298 35.453 > 105 106 Cl 0.64423742496161 0.12500000000000 0.09468739479298 35.453 > 105 107 Cl 0.14423742496161 0.62500000000000 0.09468739479298 35.453 > 105 108 Cl 0.64423742496161 0.62500000000000 0.09468739479298 35.453 > 105 109 Cl 0.35576257503839 0.37500000000000 0.90531260520702 35.453 > 109 110 Cl 0.85576257503839 0.37500000000000 0.90531260520702 35.453 > 109 111 Cl 0.35576257503839 0.87500000000000 0.90531260520702 35.453 > 109 112 Cl 0.85576257503839 0.87500000000000 0.90531260520702 35.453 > 109 113 Cl 0.48991308858784 0.37500000000000 0.20226779441939 35.453 > 113 114 Cl 0.98991308858784 0.37500000000000 0.20226779441939 35.453 > 113 115 Cl 0.48991308858784 0.87500000000000 0.20226779441939 35.453 > 113 116 Cl 0.98991308858784 0.87500000000000 0.20226779441939 35.453 > 113 117 Cl 0.01008691141216 0.12500000000000 0.79773220558061 35.453 > 117 118 Cl 0.51008691141216 0.12500000000000 0.79773220558061 35.453 > 117 119 Cl 0.01008691141216 0.62500000000000 0.79773220558061 35.453 > 117 120 Cl 0.51008691141216 0.62500000000000 0.79773220558061 35.453 > 117 121 Cl 0.23991308858784 0.37500000000000 0.29773220558061 35.453 > 121 122 Cl 0.73991308858784 0.37500000000000 0.29773220558061 35.453 > 121 123 Cl 0.23991308858784 0.87500000000000 0.29773220558061 35.453 > 121 124 Cl 0.73991308858784 0.87500000000000 0.29773220558061 35.453 > 121 125 Cl 0.26008691141216 0.12500000000000 0.70226779441939 35.453 > 125 126 Cl 0.76008691141216 0.12500000000000 0.70226779441939 35.453 > 125 127 Cl 0.26008691141216 0.62500000000000 0.70226779441939 35.453 > 125 128 Cl 0.76008691141216 0.62500000000000 0.70226779441939 35.453 > 125 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 2.8609084 0.0000000 0.0000000 0.0000000 2.9279013 0.0000000 0.0000000 0.0000000 2.8719783 -------------------------- Born effective charges -------------------------- 1 Cs 1.2061394 0.0000000 0.0228863 0.0000000 1.3701435 0.0000000 -0.0198765 0.0000000 1.5198233 2 Cs 1.2061394 0.0000000 0.0228863 0.0000000 1.3701435 0.0000000 -0.0198765 0.0000000 1.5198233 3 Cs 1.2061394 0.0000000 -0.0228863 0.0000000 1.3701435 0.0000000 0.0198765 0.0000000 1.5198233 4 Cs 1.2061394 0.0000000 -0.0228863 0.0000000 1.3701435 0.0000000 0.0198765 0.0000000 1.5198233 5 Cs 1.4024399 0.0000000 -0.1171344 0.0000000 1.2964748 0.0000000 -0.1117394 0.0000000 1.2870865 6 Cs 1.4024399 0.0000000 -0.1171344 0.0000000 1.2964748 0.0000000 -0.1117394 0.0000000 1.2870865 7 Cs 1.4024399 0.0000000 0.1171344 0.0000000 1.2964748 0.0000000 0.1117394 0.0000000 1.2870865 8 Cs 1.4024399 0.0000000 0.1171344 0.0000000 1.2964748 0.0000000 0.1117394 0.0000000 1.2870865 9 Dy 3.4787294 0.0000000 -0.0268019 0.0000000 3.7567839 0.0000000 -0.0207386 0.0000000 3.3159223 10 Dy 3.4787294 0.0000000 -0.0268019 0.0000000 3.7567839 0.0000000 -0.0207386 0.0000000 3.3159223 11 Dy 3.4787294 0.0000000 0.0268019 0.0000000 3.7567839 0.0000000 0.0207386 0.0000000 3.3159223 12 Dy 3.4787294 0.0000000 0.0268019 0.0000000 3.7567839 0.0000000 0.0207386 0.0000000 3.3159223 13 Cl -1.4202677 -0.7175615 0.4011783 -0.8319935 -1.6733687 0.4384569 0.4046592 0.3664735 -0.9603484 14 Cl -1.4202677 0.7175615 0.4011783 0.8319935 -1.6733687 -0.4384569 0.4046592 -0.3664735 -0.9603484 15 Cl -1.4202677 0.7175615 -0.4011783 0.8319935 -1.6733687 0.4384569 -0.4046592 0.3664735 -0.9603484 16 Cl -1.4202677 -0.7175615 -0.4011783 -0.8319935 -1.6733687 -0.4384569 -0.4046592 -0.3664735 -0.9603484 17 Cl -1.2569846 -0.4347298 0.2330427 -0.4134532 -1.4341400 0.2805469 0.2291280 0.2751855 -0.9278093 18 Cl -1.2569846 0.4347298 0.2330427 0.4134532 -1.4341400 -0.2805469 0.2291280 -0.2751855 -0.9278093 19 Cl -1.2569846 0.4347298 -0.2330427 0.4134532 -1.4341400 0.2805469 -0.2291280 0.2751855 -0.9278093 20 Cl -1.2569846 -0.4347298 -0.2330427 -0.4134532 -1.4341400 -0.2805469 -0.2291280 -0.2751855 -0.9278093 21 Cl -1.2569846 -0.4347298 0.2330427 -0.4134532 -1.4341400 0.2805469 0.2291280 0.2751855 -0.9278093 22 Cl -1.2569846 0.4347298 0.2330427 0.4134532 -1.4341400 -0.2805469 0.2291280 -0.2751855 -0.9278093 23 Cl -1.2569846 0.4347298 -0.2330427 0.4134532 -1.4341400 0.2805469 -0.2291280 0.2751855 -0.9278093 24 Cl -1.2569846 -0.4347298 -0.2330427 -0.4134532 -1.4341400 -0.2805469 -0.2291280 -0.2751855 -0.9278093 25 Cl -1.0553548 0.0000000 0.4681014 0.0000000 -0.9481055 0.0000000 0.2972869 0.0000000 -1.7296133 26 Cl -1.0553548 0.0000000 0.4681014 0.0000000 -0.9481055 0.0000000 0.2972869 0.0000000 -1.7296133 27 Cl -1.0553548 0.0000000 -0.4681014 0.0000000 -0.9481055 0.0000000 -0.2972869 0.0000000 -1.7296133 28 Cl -1.0553548 0.0000000 -0.4681014 0.0000000 -0.9481055 0.0000000 -0.2972869 0.0000000 -1.7296133 29 Cl -1.0977170 0.0000000 -0.3435436 0.0000000 -0.9336480 0.0000000 -0.4166316 0.0000000 -1.5772518 30 Cl -1.0977170 0.0000000 -0.3435436 0.0000000 -0.9336480 0.0000000 -0.4166316 0.0000000 -1.5772518 31 Cl -1.0977170 0.0000000 0.3435436 0.0000000 -0.9336480 0.0000000 0.4166316 0.0000000 -1.5772518 32 Cl -1.0977170 0.0000000 0.3435436 0.0000000 -0.9336480 0.0000000 0.4166316 0.0000000 -1.5772518 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000001 (xzx) 0.00000001 (xzx) 0.00000001 (xxz) 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 ... [ 5 7 3 ] 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.41, Number of G-points: 285, Lambda: 0.07 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/24) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 24 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.453 ( 0.000 0.000 0.000) 0.000 0.631 ( 0.000 0.000 0.000) 0.000 0.712 ( 0.000 0.000 0.000) 0.000 0.752 ( 0.000 0.000 0.000) 0.000 0.820 ( 0.000 0.000 0.000) 0.000 0.960 ( 0.000 0.000 0.000) 0.000 1.011 ( 0.000 0.000 0.000) 0.000 1.028 ( 0.000 0.000 0.000) 0.000 1.099 ( 0.000 0.000 0.000) 0.000 1.123 ( 0.000 0.000 0.000) 0.000 1.137 ( 0.000 0.000 0.000) 0.000 1.197 ( 0.000 0.000 0.000) 0.000 1.200 ( 0.000 0.000 0.000) 0.000 1.215 ( 0.000 0.000 0.000) 0.000 1.235 ( 0.000 0.000 0.000) 0.000 1.285 ( 0.000 0.000 0.000) 0.000 1.321 ( 0.000 0.000 0.000) 0.000 1.325 ( 0.000 0.000 0.000) 0.000 1.398 ( 0.000 0.000 0.000) 0.000 1.417 ( 0.000 0.000 0.000) 0.000 1.425 ( 0.000 0.000 0.000) 0.000 1.475 ( 0.000 0.000 0.000) 0.000 1.508 ( 0.000 0.000 0.000) 0.000 1.508 ( 0.000 0.000 0.000) 0.000 1.513 ( 0.000 0.000 0.000) 0.000 1.522 ( 0.000 0.000 0.000) 0.000 1.576 ( 0.000 0.000 0.000) 0.000 1.599 ( 0.000 0.000 0.000) 0.000 1.604 ( 0.000 0.000 0.000) 0.000 1.628 ( 0.000 0.000 0.000) 0.000 1.638 ( 0.000 0.000 0.000) 0.000 1.648 ( 0.000 0.000 0.000) 0.000 1.747 ( 0.000 0.000 0.000) 0.000 1.764 ( 0.000 0.000 0.000) 0.000 1.906 ( 0.000 0.000 0.000) 0.000 1.939 ( 0.000 0.000 0.000) 0.000 1.964 ( 0.000 0.000 0.000) 0.000 1.992 ( 0.000 0.000 0.000) 0.000 2.060 ( 0.000 0.000 0.000) 0.000 2.082 ( 0.000 0.000 0.000) 0.000 2.097 ( 0.000 0.000 0.000) 0.000 2.300 ( 0.000 0.000 0.000) 0.000 2.403 ( 0.000 0.000 0.000) 0.000 2.421 ( 0.000 0.000 0.000) 0.000 2.476 ( 0.000 0.000 0.000) 0.000 2.482 ( 0.000 0.000 0.000) 0.000 2.543 ( 0.000 0.000 0.000) 0.000 2.565 ( 0.000 0.000 0.000) 0.000 2.568 ( 0.000 0.000 0.000) 0.000 2.709 ( 0.000 0.000 0.000) 0.000 2.735 ( 0.000 0.000 0.000) 0.000 2.807 ( 0.000 0.000 0.000) 0.000 2.888 ( 0.000 0.000 0.000) 0.000 2.995 ( 0.000 0.000 0.000) 0.000 3.005 ( 0.000 0.000 0.000) 0.000 3.118 ( 0.000 0.000 0.000) 0.000 3.124 ( 0.000 0.000 0.000) 0.000 3.133 ( 0.000 0.000 0.000) 0.000 3.152 ( 0.000 0.000 0.000) 0.000 3.244 ( 0.000 0.000 0.000) 0.000 3.311 ( 0.000 0.000 0.000) 0.000 3.353 ( 0.000 0.000 0.000) 0.000 3.355 ( 0.000 0.000 0.000) 0.000 3.473 ( 0.000 0.000 0.000) 0.000 3.664 ( 0.000 0.000 0.000) 0.000 3.677 ( 0.000 0.000 0.000) 0.000 3.714 ( 0.000 0.000 0.000) 0.000 3.786 ( 0.000 0.000 0.000) 0.000 3.816 ( 0.000 0.000 0.000) 0.000 3.842 ( 0.000 0.000 0.000) 0.000 3.852 ( 0.000 0.000 0.000) 0.000 3.948 ( 0.000 0.000 0.000) 0.000 4.024 ( 0.000 0.000 0.000) 0.000 6.599 ( 0.000 0.000 0.000) 0.000 6.618 ( 0.000 0.000 0.000) 0.000 6.819 ( 0.000 0.000 0.000) 0.000 6.820 ( 0.000 0.000 0.000) 0.000 6.878 ( 0.000 0.000 0.000) 0.000 6.891 ( 0.000 0.000 0.000) 0.000 7.229 ( 0.000 0.000 0.000) 0.000 7.445 ( 0.000 0.000 0.000) 0.000 7.471 ( 0.000 0.000 0.000) 0.000 7.475 ( 0.000 0.000 0.000) 0.000 7.498 ( 0.000 0.000 0.000) 0.000 7.566 ( 0.000 0.000 0.000) 0.000 7.572 ( 0.000 0.000 0.000) 0.000 7.773 ( 0.000 0.000 0.000) 0.000 7.792 ( 0.000 0.000 0.000) 0.000 8.106 ( 0.000 0.000 0.000) 0.000 8.423 ( 0.000 0.000 0.000) 0.000 8.493 ( 0.000 0.000 0.000) 0.000 8.720 ( 0.000 0.000 0.000) 0.000 8.725 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/24) ======================= q-point: ( 0.20 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.304 ( 14.049 0.000 0.000) 14.049 0.308 ( 13.424 0.000 0.000) 13.424 0.485 ( 21.942 0.000 0.000) 21.942 0.557 ( 8.755 0.000 0.000) 8.755 0.630 ( -0.604 0.000 0.000) 0.604 0.765 ( 1.170 0.000 0.000) 1.170 0.795 ( 7.030 0.000 0.000) 7.030 0.838 ( 1.973 0.000 0.000) 1.973 0.964 ( -2.434 0.000 0.000) 2.434 0.975 ( -4.578 0.000 0.000) 4.578 1.050 ( 0.621 0.000 0.000) 0.621 1.108 ( 0.646 0.000 0.000) 0.646 1.126 ( -3.795 0.000 0.000) 3.795 1.136 ( 0.709 0.000 0.000) 0.709 1.136 ( -0.106 0.000 0.000) 0.106 1.189 ( -0.537 0.000 0.000) 0.537 1.190 ( -3.981 0.000 0.000) 3.981 1.242 ( -5.340 0.000 0.000) 5.340 1.265 ( -1.715 0.000 0.000) 1.715 1.283 ( 2.083 0.000 0.000) 2.083 1.308 ( -0.941 0.000 0.000) 0.941 1.361 ( -4.591 0.000 0.000) 4.591 1.421 ( 0.548 0.000 0.000) 0.548 1.459 ( -2.084 0.000 0.000) 2.084 1.460 ( -1.095 0.000 0.000) 1.095 1.471 ( -1.990 0.000 0.000) 1.990 1.542 ( 0.402 0.000 0.000) 0.402 1.545 ( 0.416 0.000 0.000) 0.416 1.558 ( 0.131 0.000 0.000) 0.131 1.564 ( 2.239 0.000 0.000) 2.239 1.605 ( -2.334 0.000 0.000) 2.334 1.632 ( 0.392 0.000 0.000) 0.392 1.649 ( 2.206 0.000 0.000) 2.206 1.652 ( 2.867 0.000 0.000) 2.867 1.737 ( -1.136 0.000 0.000) 1.136 1.778 ( 1.824 0.000 0.000) 1.824 1.791 ( 2.870 0.000 0.000) 2.870 1.906 ( -0.429 0.000 0.000) 0.429 1.966 ( -1.764 0.000 0.000) 1.764 1.980 ( 1.555 0.000 0.000) 1.555 2.000 ( -2.504 0.000 0.000) 2.504 2.042 ( 0.534 0.000 0.000) 0.534 2.124 ( 0.874 0.000 0.000) 0.874 2.161 ( 5.230 0.000 0.000) 5.230 2.308 ( 1.589 0.000 0.000) 1.589 2.361 ( -3.715 0.000 0.000) 3.715 2.449 ( 1.922 0.000 0.000) 1.922 2.460 ( -1.245 0.000 0.000) 1.245 2.471 ( -0.924 0.000 0.000) 0.924 2.527 ( -1.170 0.000 0.000) 1.170 2.554 ( -1.047 0.000 0.000) 1.047 2.555 ( -1.128 0.000 0.000) 1.128 2.731 ( 1.636 0.000 0.000) 1.636 2.745 ( 0.742 0.000 0.000) 0.742 2.805 ( -0.124 0.000 0.000) 0.124 2.875 ( -1.465 0.000 0.000) 1.465 2.993 ( -0.491 0.000 0.000) 0.491 3.012 ( 1.440 0.000 0.000) 1.440 3.076 ( -1.631 0.000 0.000) 1.631 3.110 ( -1.997 0.000 0.000) 1.997 3.230 ( 3.148 0.000 0.000) 3.148 3.269 ( 3.268 0.000 0.000) 3.268 3.301 ( 1.383 0.000 0.000) 1.383 3.317 ( 0.535 0.000 0.000) 0.535 3.352 ( -0.155 0.000 0.000) 0.155 3.438 ( -2.947 0.000 0.000) 2.947 3.645 ( -2.423 0.000 0.000) 2.423 3.676 ( 0.952 0.000 0.000) 0.952 3.696 ( 0.840 0.000 0.000) 0.840 3.715 ( 0.023 0.000 0.000) 0.023 3.745 ( -3.254 0.000 0.000) 3.254 3.803 ( -1.878 0.000 0.000) 1.878 3.904 ( 2.773 0.000 0.000) 2.773 3.957 ( 0.834 0.000 0.000) 0.834 3.988 ( 0.273 0.000 0.000) 0.273 4.022 ( -0.345 0.000 0.000) 0.345 6.619 ( 1.724 0.000 0.000) 1.724 6.760 ( -2.365 0.000 0.000) 2.365 6.820 ( 0.084 0.000 0.000) 0.084 6.821 ( 0.070 0.000 0.000) 0.070 6.925 ( 3.298 0.000 0.000) 3.298 6.943 ( 4.408 0.000 0.000) 4.408 7.221 ( -0.518 0.000 0.000) 0.518 7.428 ( -1.739 0.000 0.000) 1.739 7.477 ( 0.484 0.000 0.000) 0.484 7.481 ( 0.489 0.000 0.000) 0.489 7.562 ( -0.416 0.000 0.000) 0.416 7.565 ( -0.600 0.000 0.000) 0.600 7.643 ( -4.843 0.000 0.000) 4.843 7.872 ( -8.256 0.000 0.000) 8.256 7.946 ( 7.309 0.000 0.000) 7.309 8.286 ( -5.397 0.000 0.000) 5.397 8.507 ( 1.249 0.000 0.000) 1.249 8.533 ( 3.220 0.000 0.000) 3.220 8.723 ( 0.263 0.000 0.000) 0.263 8.729 ( 0.039 0.000 0.000) 0.039 ======================= Grid point 2 (3/24) ======================= q-point: ( 0.40 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.536 ( 7.030 0.000 0.000) 7.030 0.581 ( 11.497 0.000 0.000) 11.497 0.607 ( -1.316 0.000 0.000) 1.316 0.755 ( -3.757 0.000 0.000) 3.757 0.770 ( 10.348 0.000 0.000) 10.348 0.883 ( 0.903 0.000 0.000) 0.903 0.888 ( -2.012 0.000 0.000) 2.012 0.907 ( 16.038 0.000 0.000) 16.038 0.928 ( -0.445 0.000 0.000) 0.445 0.961 ( 7.154 0.000 0.000) 7.154 0.990 ( -5.546 0.000 0.000) 5.546 1.033 ( -4.752 0.000 0.000) 4.752 1.097 ( -1.473 0.000 0.000) 1.473 1.107 ( -2.281 0.000 0.000) 2.281 1.123 ( -1.869 0.000 0.000) 1.869 1.134 ( 0.223 0.000 0.000) 0.223 1.152 ( -1.993 0.000 0.000) 1.993 1.218 ( 2.337 0.000 0.000) 2.337 1.220 ( -8.777 0.000 0.000) 8.777 1.260 ( 1.409 0.000 0.000) 1.409 1.278 ( -2.823 0.000 0.000) 2.823 1.290 ( -1.021 0.000 0.000) 1.021 1.425 ( -0.927 0.000 0.000) 0.927 1.430 ( -1.391 0.000 0.000) 1.391 1.436 ( 0.353 0.000 0.000) 0.353 1.441 ( -0.503 0.000 0.000) 0.503 1.510 ( -1.546 0.000 0.000) 1.546 1.512 ( -1.659 0.000 0.000) 1.659 1.573 ( 0.884 0.000 0.000) 0.884 1.599 ( -1.904 0.000 0.000) 1.904 1.604 ( 0.963 0.000 0.000) 0.963 1.607 ( 1.025 0.000 0.000) 1.025 1.707 ( 2.223 0.000 0.000) 2.223 1.726 ( 4.189 0.000 0.000) 4.189 1.736 ( -0.245 0.000 0.000) 0.245 1.793 ( -1.128 0.000 0.000) 1.128 1.848 ( 1.870 0.000 0.000) 1.870 1.873 ( -1.666 0.000 0.000) 1.666 1.945 ( -0.486 0.000 0.000) 0.486 1.949 ( -1.503 0.000 0.000) 1.503 2.019 ( 1.812 0.000 0.000) 1.812 2.051 ( -0.860 0.000 0.000) 0.860 2.153 ( 3.460 0.000 0.000) 3.460 2.249 ( -3.385 0.000 0.000) 3.385 2.291 ( 3.066 0.000 0.000) 3.066 2.343 ( 0.142 0.000 0.000) 0.142 2.439 ( -0.369 0.000 0.000) 0.369 2.448 ( -1.025 0.000 0.000) 1.025 2.486 ( 1.643 0.000 0.000) 1.643 2.512 ( 0.055 0.000 0.000) 0.055 2.524 ( -1.771 0.000 0.000) 1.771 2.527 ( -1.195 0.000 0.000) 1.195 2.758 ( 0.660 0.000 0.000) 0.660 2.759 ( 0.432 0.000 0.000) 0.432 2.805 ( 0.323 0.000 0.000) 0.323 2.830 ( -2.107 0.000 0.000) 2.107 3.006 ( 1.717 0.000 0.000) 1.717 3.043 ( 1.203 0.000 0.000) 1.203 3.057 ( -2.772 0.000 0.000) 2.772 3.059 ( -0.394 0.000 0.000) 0.394 3.280 ( 2.523 0.000 0.000) 2.523 3.326 ( -0.204 0.000 0.000) 0.204 3.332 ( 0.735 0.000 0.000) 0.735 3.345 ( -0.495 0.000 0.000) 0.495 3.349 ( 2.449 0.000 0.000) 2.449 3.373 ( -1.996 0.000 0.000) 1.996 3.610 ( 0.062 0.000 0.000) 0.062 3.658 ( -3.715 0.000 0.000) 3.715 3.697 ( 0.947 0.000 0.000) 0.947 3.712 ( -0.388 0.000 0.000) 0.388 3.728 ( 1.396 0.000 0.000) 1.396 3.763 ( -1.851 0.000 0.000) 1.851 3.956 ( 2.001 0.000 0.000) 2.001 3.980 ( 1.272 0.000 0.000) 1.272 3.984 ( -0.683 0.000 0.000) 0.683 4.007 ( -1.081 0.000 0.000) 1.081 6.662 ( 2.226 0.000 0.000) 2.226 6.710 ( -2.294 0.000 0.000) 2.294 6.822 ( 0.065 0.000 0.000) 0.065 6.822 ( 0.030 0.000 0.000) 0.030 6.995 ( 3.013 0.000 0.000) 3.013 7.025 ( 1.389 0.000 0.000) 1.389 7.243 ( 3.971 0.000 0.000) 3.971 7.359 ( -5.001 0.000 0.000) 5.001 7.489 ( 0.516 0.000 0.000) 0.516 7.492 ( 0.315 0.000 0.000) 0.315 7.552 ( -0.390 0.000 0.000) 0.390 7.552 ( -0.466 0.000 0.000) 0.466 7.595 ( 0.646 0.000 0.000) 0.646 7.682 ( -7.640 0.000 0.000) 7.640 8.084 ( 5.642 0.000 0.000) 5.642 8.188 ( -4.401 0.000 0.000) 4.401 8.538 ( 1.500 0.000 0.000) 1.500 8.560 ( -0.249 0.000 0.000) 0.249 8.729 ( 0.208 0.000 0.000) 0.208 8.731 ( 0.024 0.000 0.000) 0.024 ======================= Grid point 5 (4/24) ======================= q-point: ( 0.00 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.276 ( 0.000 14.151 0.000) 14.151 0.293 ( 0.000 15.315 0.000) 15.315 0.480 ( 0.000 3.186 0.000) 3.186 0.565 ( 0.000 27.199 0.000) 27.199 0.734 ( 0.000 2.314 0.000) 2.314 0.752 ( 0.000 0.193 0.000) 0.193 0.806 ( 0.000 13.766 0.000) 13.766 0.924 ( 0.000 8.673 0.000) 8.673 0.931 ( 0.000 -2.464 0.000) 2.464 1.000 ( 0.000 -0.981 0.000) 0.981 1.089 ( 0.000 4.216 0.000) 4.216 1.121 ( 0.000 0.370 0.000) 0.370 1.129 ( 0.000 2.799 0.000) 2.799 1.156 ( 0.000 1.641 0.000) 1.641 1.180 ( 0.000 -0.589 0.000) 0.589 1.188 ( 0.000 -1.141 0.000) 1.141 1.204 ( 0.000 -0.922 0.000) 0.922 1.249 ( 0.000 0.037 0.000) 0.037 1.281 ( 0.000 -0.310 0.000) 0.310 1.317 ( 0.000 -0.772 0.000) 0.772 1.343 ( 0.000 -0.265 0.000) 0.265 1.379 ( 0.000 -1.764 0.000) 1.764 1.391 ( 0.000 -0.687 0.000) 0.687 1.426 ( 0.000 -0.945 0.000) 0.945 1.467 ( 0.000 -2.151 0.000) 2.151 1.500 ( 0.000 1.171 0.000) 1.171 1.504 ( 0.000 -0.354 0.000) 0.354 1.533 ( 0.000 1.966 0.000) 1.966 1.535 ( 0.000 1.347 0.000) 1.347 1.599 ( 0.000 0.024 0.000) 0.024 1.603 ( 0.000 -1.879 0.000) 1.879 1.644 ( 0.000 0.896 0.000) 0.896 1.645 ( 0.000 1.506 0.000) 1.506 1.657 ( 0.000 2.311 0.000) 2.311 1.798 ( 0.000 5.199 0.000) 5.199 1.804 ( 0.000 2.434 0.000) 2.434 1.859 ( 0.000 -1.731 0.000) 1.731 1.911 ( 0.000 1.273 0.000) 1.273 1.973 ( 0.000 -1.340 0.000) 1.340 1.975 ( 0.000 3.574 0.000) 3.574 2.027 ( 0.000 5.482 0.000) 5.482 2.103 ( 0.000 1.591 0.000) 1.591 2.115 ( 0.000 5.177 0.000) 5.177 2.152 ( 0.000 4.155 0.000) 4.155 2.335 ( 0.000 3.531 0.000) 3.531 2.420 ( 0.000 0.425 0.000) 0.425 2.436 ( 0.000 3.442 0.000) 3.442 2.467 ( 0.000 -1.181 0.000) 1.181 2.473 ( 0.000 -0.609 0.000) 0.609 2.527 ( 0.000 -1.584 0.000) 1.584 2.560 ( 0.000 -0.377 0.000) 0.377 2.609 ( 0.000 -0.126 0.000) 0.126 2.729 ( 0.000 2.181 0.000) 2.181 2.791 ( 0.000 4.595 0.000) 4.595 2.933 ( 0.000 3.244 0.000) 3.244 2.980 ( 0.000 -1.262 0.000) 1.262 3.003 ( 0.000 -0.459 0.000) 0.459 3.089 ( 0.000 -1.901 0.000) 1.901 3.119 ( 0.000 -2.760 0.000) 2.760 3.130 ( 0.000 -4.374 0.000) 4.374 3.137 ( 0.000 1.706 0.000) 1.706 3.184 ( 0.000 -1.070 0.000) 1.070 3.241 ( 0.000 -3.087 0.000) 3.087 3.408 ( 0.000 2.723 0.000) 2.723 3.417 ( 0.000 -4.793 0.000) 4.793 3.424 ( 0.000 3.585 0.000) 3.585 3.561 ( 0.000 -5.945 0.000) 5.945 3.618 ( 0.000 0.386 0.000) 0.386 3.727 ( 0.000 -1.956 0.000) 1.956 3.756 ( 0.000 3.311 0.000) 3.311 3.788 ( 0.000 2.716 0.000) 2.716 3.794 ( 0.000 -2.191 0.000) 2.191 3.828 ( 0.000 -0.737 0.000) 0.737 3.843 ( 0.000 -0.731 0.000) 0.731 3.944 ( 0.000 -0.396 0.000) 0.396 4.023 ( 0.000 -0.083 0.000) 0.083 6.603 ( 0.000 0.459 0.000) 0.459 6.627 ( 0.000 0.912 0.000) 0.912 6.763 ( 0.000 -3.293 0.000) 3.293 6.858 ( 0.000 -1.986 0.000) 1.986 6.956 ( 0.000 3.635 0.000) 3.635 7.068 ( 0.000 1.763 0.000) 1.763 7.190 ( 0.000 -3.554 0.000) 3.554 7.425 ( 0.000 -5.609 0.000) 5.609 7.435 ( 0.000 -0.899 0.000) 0.899 7.514 ( 0.000 3.146 0.000) 3.146 7.573 ( 0.000 0.104 0.000) 0.104 7.580 ( 0.000 -0.258 0.000) 0.258 7.817 ( 0.000 4.747 0.000) 4.747 7.823 ( 0.000 2.844 0.000) 2.844 8.105 ( 0.000 0.122 0.000) 0.122 8.309 ( 0.000 0.329 0.000) 0.329 8.412 ( 0.000 -0.354 0.000) 0.354 8.453 ( 0.000 -3.618 0.000) 3.618 8.708 ( 0.000 -1.331 0.000) 1.331 8.712 ( 0.000 -1.358 0.000) 1.358 ======================= Grid point 6 (5/24) ======================= q-point: ( 0.20 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.405 ( 10.036 9.173 0.000) 13.596 0.435 ( 11.255 8.169 0.000) 13.907 0.581 ( 8.151 2.604 0.000) 8.557 0.721 ( 11.632 19.984 0.000) 23.123 0.778 ( 2.480 2.024 0.000) 3.201 0.791 ( 1.925 6.598 0.000) 6.874 0.806 ( 2.794 8.064 0.000) 8.535 0.943 ( -0.649 -1.937 0.000) 2.043 0.948 ( 2.437 8.716 0.000) 9.050 0.997 ( 0.156 -1.811 0.000) 1.818 1.037 ( -4.022 4.414 0.000) 5.972 1.118 ( -3.035 -0.260 0.000) 3.047 1.125 ( -1.361 -0.826 0.000) 1.592 1.144 ( 0.377 0.341 0.000) 0.509 1.144 ( -0.225 2.481 0.000) 2.491 1.167 ( -3.689 -1.570 0.000) 4.009 1.190 ( 0.080 -0.476 0.000) 0.483 1.260 ( 0.853 -1.562 0.000) 1.780 1.260 ( -2.821 2.097 0.000) 3.515 1.275 ( -0.949 0.937 0.000) 1.334 1.343 ( -1.158 2.078 0.000) 2.379 1.358 ( -2.023 0.155 0.000) 2.029 1.436 ( 0.766 -1.397 0.000) 1.593 1.439 ( -0.051 -0.597 0.000) 0.600 1.467 ( -0.803 -0.965 0.000) 1.256 1.481 ( -2.764 1.347 0.000) 3.075 1.490 ( 0.723 -3.104 0.000) 3.187 1.521 ( 2.366 -1.308 0.000) 2.703 1.571 ( 0.125 2.465 0.000) 2.468 1.577 ( -1.182 0.474 0.000) 1.273 1.580 ( -0.043 0.022 0.000) 0.048 1.641 ( 2.616 1.813 0.000) 3.183 1.662 ( -0.990 1.009 0.000) 1.414 1.691 ( 1.612 0.640 0.000) 1.734 1.791 ( -0.330 3.983 0.000) 3.997 1.823 ( 1.254 4.421 0.000) 4.596 1.850 ( 1.174 2.641 0.000) 2.890 1.901 ( -1.059 0.098 0.000) 1.064 1.950 ( -1.375 -0.936 0.000) 1.664 1.969 ( -1.129 1.783 0.000) 2.111 2.060 ( 1.089 4.277 0.000) 4.413 2.082 ( -0.006 3.023 0.000) 3.023 2.155 ( 1.325 3.278 0.000) 3.535 2.198 ( 4.333 3.274 0.000) 5.431 2.347 ( 1.351 3.586 0.000) 3.832 2.395 ( -3.368 3.431 0.000) 4.808 2.431 ( 0.501 -0.463 0.000) 0.682 2.458 ( -0.308 -0.852 0.000) 0.906 2.461 ( -1.149 -0.518 0.000) 1.260 2.505 ( -1.772 -1.914 0.000) 2.608 2.551 ( -0.694 0.046 0.000) 0.696 2.581 ( -1.788 1.173 0.000) 2.139 2.756 ( 2.033 2.455 0.000) 3.188 2.796 ( 0.390 4.455 0.000) 4.472 2.913 ( -1.901 2.939 0.000) 3.500 2.942 ( -4.464 5.247 0.000) 6.889 2.994 ( -0.327 0.030 0.000) 0.329 3.020 ( -0.454 2.786 0.000) 2.822 3.079 ( -2.079 0.645 0.000) 2.177 3.092 ( -1.828 -1.659 0.000) 2.469 3.163 ( 1.507 -7.352 0.000) 7.505 3.237 ( -0.197 -5.540 0.000) 5.544 3.247 ( 3.821 1.495 0.000) 4.103 3.329 ( -0.469 1.449 0.000) 1.523 3.406 ( -1.544 -3.047 0.000) 3.416 3.408 ( -0.042 2.492 0.000) 2.492 3.558 ( -0.162 -5.442 0.000) 5.445 3.688 ( 0.022 -1.923 0.000) 1.923 3.706 ( -0.630 -1.561 0.000) 1.683 3.744 ( -0.688 3.375 0.000) 3.444 3.765 ( -0.990 3.440 0.000) 3.580 3.801 ( -1.617 0.159 0.000) 1.624 3.881 ( 2.233 -2.015 0.000) 3.008 3.900 ( 3.923 -5.890 0.000) 7.077 3.951 ( 0.693 -0.571 0.000) 0.898 4.018 ( -0.551 -0.374 0.000) 0.666 6.619 ( 1.425 0.065 0.000) 1.427 6.670 ( 2.113 -2.498 0.000) 3.272 6.769 ( 0.426 -3.484 0.000) 3.510 6.828 ( -2.240 -1.408 0.000) 2.646 6.994 ( 3.360 3.083 0.000) 4.561 7.029 ( -0.739 4.845 0.000) 4.901 7.185 ( -0.242 -3.198 0.000) 3.207 7.386 ( -2.126 -3.816 0.000) 4.369 7.500 ( 2.560 -3.066 0.000) 3.994 7.516 ( 0.175 2.641 0.000) 2.647 7.554 ( -1.623 2.166 0.000) 2.707 7.564 ( -0.700 -0.024 0.000) 0.700 7.706 ( -3.767 6.627 0.000) 7.623 7.930 ( -9.701 5.881 0.000) 11.345 7.961 ( 6.438 1.378 0.000) 6.584 8.272 ( -3.789 -1.096 0.000) 3.944 8.463 ( 0.890 -4.031 0.000) 4.128 8.483 ( 2.917 -3.711 0.000) 4.720 8.711 ( 0.259 -1.339 0.000) 1.364 8.714 ( 0.147 -1.526 0.000) 1.533 ======================= Grid point 7 (6/24) ======================= q-point: ( 0.40 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.631 ( 10.055 5.392 0.000) 11.409 0.670 ( 8.715 7.745 0.000) 11.659 0.766 ( 6.776 3.326 0.000) 7.549 0.778 ( 0.787 9.780 0.000) 9.812 0.796 ( -2.564 3.872 0.000) 4.644 0.908 ( -0.245 1.167 0.000) 1.193 0.927 ( 3.892 1.450 0.000) 4.153 0.927 ( -2.397 2.741 0.000) 3.641 0.972 ( -0.923 0.443 0.000) 1.023 0.983 ( -0.812 2.608 0.000) 2.731 0.992 ( 4.492 5.738 0.000) 7.287 1.045 ( -3.222 1.481 0.000) 3.546 1.093 ( -0.036 0.667 0.000) 0.668 1.097 ( 1.320 0.733 0.000) 1.510 1.161 ( 0.775 2.996 0.000) 3.095 1.162 ( -0.024 2.735 0.000) 2.735 1.178 ( -0.739 1.244 0.000) 1.447 1.235 ( 0.861 2.200 0.000) 2.363 1.258 ( -6.050 1.877 0.000) 6.334 1.260 ( -0.695 -0.946 0.000) 1.174 1.263 ( 0.384 0.624 0.000) 0.732 1.311 ( -3.379 2.911 0.000) 4.460 1.430 ( 0.214 -0.028 0.000) 0.215 1.432 ( -0.455 0.465 0.000) 0.651 1.448 ( -0.119 -0.698 0.000) 0.708 1.462 ( 0.805 -0.185 0.000) 0.825 1.492 ( -1.001 -2.127 0.000) 2.351 1.504 ( -3.747 0.419 0.000) 3.770 1.558 ( 1.582 0.091 0.000) 1.584 1.590 ( 0.872 1.413 0.000) 1.660 1.615 ( -1.822 0.996 0.000) 2.076 1.620 ( 0.432 0.133 0.000) 0.452 1.711 ( 4.198 1.651 0.000) 4.511 1.745 ( 3.172 0.927 0.000) 3.305 1.778 ( -1.540 3.133 0.000) 3.491 1.816 ( -2.901 1.484 0.000) 3.259 1.885 ( 0.864 2.019 0.000) 2.196 1.893 ( 1.967 2.770 0.000) 3.398 1.929 ( -2.050 0.521 0.000) 2.115 1.932 ( -0.845 -0.754 0.000) 1.133 2.075 ( 0.557 4.225 0.000) 4.261 2.083 ( -0.220 2.986 0.000) 2.994 2.196 ( 3.224 3.976 0.000) 5.118 2.268 ( -1.277 1.610 0.000) 2.055 2.339 ( 0.985 5.029 0.000) 5.124 2.368 ( -0.007 2.120 0.000) 2.120 2.430 ( -0.461 -0.541 0.000) 0.711 2.433 ( -1.089 -1.193 0.000) 1.616 2.464 ( 0.460 -1.416 0.000) 1.489 2.474 ( -0.829 -2.375 0.000) 2.516 2.542 ( -0.022 1.296 0.000) 1.296 2.552 ( -0.941 1.506 0.000) 1.775 2.791 ( 0.883 3.070 0.000) 3.194 2.799 ( -0.161 3.762 0.000) 3.765 2.871 ( -1.148 3.323 0.000) 3.515 2.878 ( -1.526 5.667 0.000) 5.869 3.002 ( 1.065 -0.137 0.000) 1.074 3.030 ( -0.298 -0.663 0.000) 0.727 3.050 ( -1.103 -0.540 0.000) 1.228 3.054 ( -0.524 -0.174 0.000) 0.552 3.199 ( 1.735 -7.425 0.000) 7.625 3.228 ( -0.828 -7.108 0.000) 7.156 3.314 ( 2.484 0.165 0.000) 2.490 3.351 ( 1.163 -0.901 0.000) 1.471 3.366 ( -1.367 0.039 0.000) 1.368 3.396 ( -1.388 1.037 0.000) 1.733 3.573 ( 1.974 -3.196 0.000) 3.757 3.634 ( -3.279 -2.798 0.000) 4.311 3.728 ( 1.372 0.416 0.000) 1.433 3.738 ( 0.049 3.142 0.000) 3.143 3.747 ( -0.598 3.013 0.000) 3.071 3.764 ( -1.785 0.196 0.000) 1.796 3.923 ( 1.510 -3.052 0.000) 3.405 3.937 ( 0.243 -4.132 0.000) 4.139 3.972 ( 1.235 -0.784 0.000) 1.463 3.999 ( -1.189 -0.757 0.000) 1.409 6.656 ( 1.839 -0.580 0.000) 1.928 6.686 ( -0.611 -1.807 0.000) 1.907 6.777 ( 0.337 -3.490 0.000) 3.506 6.789 ( -1.122 -2.385 0.000) 2.636 7.044 ( 1.447 3.064 0.000) 3.389 7.058 ( 1.166 2.069 0.000) 2.374 7.209 ( 3.566 -3.058 0.000) 4.697 7.314 ( -4.828 -4.247 0.000) 6.430 7.518 ( 0.048 1.539 0.000) 1.540 7.523 ( -0.378 0.283 0.000) 0.472 7.542 ( 0.703 0.863 0.000) 1.113 7.553 ( -0.339 -0.017 0.000) 0.339 7.673 ( 0.855 8.554 0.000) 8.597 7.752 ( -6.524 7.767 0.000) 10.144 8.085 ( 5.239 -0.043 0.000) 5.239 8.185 ( -4.299 -0.512 0.000) 4.330 8.487 ( 1.190 -4.773 0.000) 4.919 8.506 ( -0.339 -4.929 0.000) 4.940 8.716 ( 0.175 -1.378 0.000) 1.389 8.717 ( 0.118 -1.465 0.000) 1.470 ======================= Grid point 10 (7/24) ======================= q-point: ( 0.00 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.544 ( 0.000 13.316 0.000) 13.316 0.578 ( 0.000 6.631 0.000) 6.631 0.582 ( 0.000 14.013 0.000) 14.013 0.766 ( 0.000 1.401 0.000) 1.401 0.793 ( 0.000 3.141 0.000) 3.141 0.891 ( 0.000 -1.239 0.000) 1.239 0.976 ( 0.000 -1.310 0.000) 1.310 0.996 ( 0.000 12.456 0.000) 12.456 1.018 ( 0.000 5.641 0.000) 5.641 1.067 ( 0.000 5.523 0.000) 5.523 1.109 ( 0.000 0.126 0.000) 0.126 1.125 ( 0.000 0.358 0.000) 0.358 1.170 ( 0.000 -4.272 0.000) 4.272 1.179 ( 0.000 -2.225 0.000) 2.225 1.192 ( 0.000 3.183 0.000) 3.183 1.201 ( 0.000 2.857 0.000) 2.857 1.212 ( 0.000 4.165 0.000) 4.165 1.252 ( 0.000 -4.629 0.000) 4.629 1.259 ( 0.000 4.781 0.000) 4.781 1.278 ( 0.000 0.426 0.000) 0.426 1.302 ( 0.000 -0.222 0.000) 0.222 1.373 ( 0.000 -1.180 0.000) 1.180 1.389 ( 0.000 -2.423 0.000) 2.423 1.410 ( 0.000 -1.690 0.000) 1.690 1.451 ( 0.000 3.560 0.000) 3.560 1.456 ( 0.000 -2.125 0.000) 2.125 1.527 ( 0.000 0.546 0.000) 0.546 1.560 ( 0.000 0.501 0.000) 0.501 1.575 ( 0.000 1.243 0.000) 1.243 1.602 ( 0.000 0.573 0.000) 0.573 1.659 ( 0.000 4.956 0.000) 4.956 1.671 ( 0.000 0.714 0.000) 0.714 1.679 ( 0.000 0.016 0.000) 0.016 1.706 ( 0.000 5.012 0.000) 5.012 1.810 ( 0.000 -1.345 0.000) 1.345 1.814 ( 0.000 -0.948 0.000) 0.948 1.862 ( 0.000 -0.288 0.000) 0.288 1.970 ( 0.000 4.027 0.000) 4.027 2.005 ( 0.000 5.815 0.000) 5.815 2.061 ( 0.000 4.768 0.000) 4.768 2.125 ( 0.000 3.602 0.000) 3.602 2.126 ( 0.000 0.470 0.000) 0.470 2.213 ( 0.000 3.259 0.000) 3.259 2.232 ( 0.000 3.060 0.000) 3.060 2.415 ( 0.000 3.235 0.000) 3.235 2.431 ( 0.000 -2.030 0.000) 2.030 2.464 ( 0.000 -0.508 0.000) 0.508 2.477 ( 0.000 6.192 0.000) 6.192 2.480 ( 0.000 -2.484 0.000) 2.484 2.543 ( 0.000 6.814 0.000) 6.814 2.569 ( 0.000 2.302 0.000) 2.302 2.639 ( 0.000 3.800 0.000) 3.800 2.788 ( 0.000 3.506 0.000) 3.506 2.857 ( 0.000 0.477 0.000) 0.477 2.963 ( 0.000 -4.085 0.000) 4.085 2.971 ( 0.000 0.610 0.000) 0.610 2.978 ( 0.000 2.364 0.000) 2.364 3.019 ( 0.000 -5.366 0.000) 5.366 3.043 ( 0.000 -2.993 0.000) 2.993 3.073 ( 0.000 -1.010 0.000) 1.010 3.187 ( 0.000 2.718 0.000) 2.718 3.203 ( 0.000 2.151 0.000) 2.151 3.207 ( 0.000 -0.812 0.000) 0.812 3.318 ( 0.000 -4.872 0.000) 4.872 3.417 ( 0.000 -2.034 0.000) 2.034 3.446 ( 0.000 -5.814 0.000) 5.814 3.467 ( 0.000 0.556 0.000) 0.556 3.578 ( 0.000 -5.196 0.000) 5.196 3.695 ( 0.000 0.312 0.000) 0.312 3.746 ( 0.000 -2.184 0.000) 2.184 3.788 ( 0.000 0.941 0.000) 0.941 3.834 ( 0.000 -0.047 0.000) 0.047 3.846 ( 0.000 2.146 0.000) 2.146 3.865 ( 0.000 4.255 0.000) 4.255 3.932 ( 0.000 -0.804 0.000) 0.804 4.017 ( 0.000 -0.723 0.000) 0.723 6.617 ( 0.000 1.004 0.000) 1.004 6.653 ( 0.000 1.810 0.000) 1.810 6.711 ( 0.000 -2.154 0.000) 2.154 6.809 ( 0.000 -2.899 0.000) 2.899 7.003 ( 0.000 1.380 0.000) 1.380 7.113 ( 0.000 2.672 0.000) 2.672 7.113 ( 0.000 -3.878 0.000) 3.878 7.315 ( 0.000 -5.295 0.000) 5.295 7.418 ( 0.000 -0.583 0.000) 0.583 7.544 ( 0.000 -3.228 0.000) 3.228 7.564 ( 0.000 1.743 0.000) 1.743 7.577 ( 0.000 0.269 0.000) 0.269 7.884 ( 0.000 3.105 0.000) 3.105 7.934 ( 0.000 5.577 0.000) 5.577 8.118 ( 0.000 1.210 0.000) 1.210 8.275 ( 0.000 -3.452 0.000) 3.452 8.392 ( 0.000 -1.188 0.000) 1.188 8.448 ( 0.000 3.718 0.000) 3.718 8.663 ( 0.000 -3.484 0.000) 3.484 8.672 ( 0.000 -2.683 0.000) 2.683 ======================= Grid point 11 (8/24) ======================= q-point: ( 0.20 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.605 ( 5.671 9.686 0.000) 11.224 0.621 ( 4.433 10.640 0.000) 11.527 0.664 ( 6.711 6.894 0.000) 9.621 0.813 ( 4.021 2.617 0.000) 4.797 0.837 ( 3.813 2.297 0.000) 4.452 0.908 ( 0.697 -0.867 0.000) 1.112 0.977 ( 0.442 -0.513 0.000) 0.677 1.014 ( -0.431 6.276 0.000) 6.291 1.050 ( 2.668 2.514 0.000) 3.666 1.090 ( -4.509 -3.028 0.000) 5.431 1.093 ( 2.176 5.369 0.000) 5.793 1.100 ( -0.758 0.693 0.000) 1.027 1.140 ( -1.880 6.532 0.000) 6.797 1.145 ( -3.250 0.313 0.000) 3.265 1.166 ( 0.786 3.310 0.000) 3.402 1.190 ( -0.372 2.649 0.000) 2.675 1.193 ( 0.494 0.966 0.000) 1.085 1.267 ( 0.431 3.614 0.000) 3.640 1.296 ( 0.769 1.750 0.000) 1.912 1.305 ( 0.459 2.810 0.000) 2.847 1.313 ( 2.474 -4.629 0.000) 5.249 1.375 ( -0.525 -0.659 0.000) 0.843 1.398 ( 1.207 -2.494 0.000) 2.770 1.412 ( -0.683 -1.500 0.000) 1.649 1.423 ( -1.938 -1.953 0.000) 2.752 1.459 ( 0.284 1.929 0.000) 1.949 1.476 ( -3.129 0.180 0.000) 3.134 1.550 ( -0.731 0.184 0.000) 0.754 1.562 ( -0.780 -0.937 0.000) 1.219 1.623 ( 0.725 -0.049 0.000) 0.727 1.644 ( -0.990 2.347 0.000) 2.548 1.679 ( -0.371 1.682 0.000) 1.722 1.691 ( 2.426 2.426 0.000) 3.431 1.721 ( 1.817 6.176 0.000) 6.438 1.829 ( 0.670 0.178 0.000) 0.693 1.853 ( 2.603 -1.104 0.000) 2.828 1.874 ( 0.926 -0.217 0.000) 0.951 1.944 ( -2.215 3.998 0.000) 4.571 1.995 ( -1.120 6.130 0.000) 6.231 2.053 ( -0.992 5.534 0.000) 5.623 2.125 ( 0.590 1.204 0.000) 1.341 2.131 ( 0.388 2.120 0.000) 2.155 2.222 ( 0.760 1.922 0.000) 2.067 2.253 ( 1.979 1.467 0.000) 2.464 2.398 ( -1.005 0.231 0.000) 1.031 2.439 ( -1.896 -0.059 0.000) 1.897 2.442 ( 0.385 -0.530 0.000) 0.655 2.462 ( -1.369 -1.973 0.000) 2.402 2.476 ( -0.045 5.884 0.000) 5.884 2.521 ( -1.863 7.296 0.000) 7.530 2.583 ( 1.196 4.280 0.000) 4.443 2.621 ( -1.437 3.715 0.000) 3.983 2.813 ( 1.981 2.969 0.000) 3.569 2.870 ( 0.899 1.271 0.000) 1.557 2.943 ( -1.492 -0.780 0.000) 1.684 2.966 ( -1.083 1.734 0.000) 2.044 2.991 ( 0.967 -0.368 0.000) 1.035 3.031 ( -0.245 -5.105 0.000) 5.111 3.036 ( 1.016 -4.727 0.000) 4.835 3.063 ( -0.855 -0.952 0.000) 1.279 3.114 ( -2.812 2.876 0.000) 4.022 3.174 ( -3.081 0.952 0.000) 3.225 3.279 ( 2.982 1.558 0.000) 3.364 3.335 ( 1.143 -3.202 0.000) 3.400 3.360 ( -4.364 0.523 0.000) 4.396 3.418 ( -0.198 -1.537 0.000) 1.549 3.447 ( 0.601 -6.037 0.000) 6.067 3.593 ( -1.392 -6.462 0.000) 6.611 3.713 ( 1.485 1.060 0.000) 1.825 3.785 ( -0.270 1.335 0.000) 1.362 3.795 ( 2.982 -2.613 0.000) 3.965 3.828 ( -1.265 2.103 0.000) 2.454 3.847 ( 0.724 -0.919 0.000) 1.170 3.855 ( -0.749 3.293 0.000) 3.377 3.936 ( 0.381 -1.011 0.000) 1.080 4.005 ( -1.032 -1.035 0.000) 1.461 6.623 ( 0.668 0.448 0.000) 0.804 6.658 ( 0.325 0.663 0.000) 0.738 6.710 ( -0.049 -2.526 0.000) 2.527 6.786 ( -1.938 -2.680 0.000) 3.307 7.033 ( 2.541 1.017 0.000) 2.737 7.095 ( -0.945 2.591 0.000) 2.758 7.118 ( 0.466 -3.322 0.000) 3.354 7.289 ( -2.391 -5.355 0.000) 5.865 7.460 ( 3.001 -0.962 0.000) 3.151 7.548 ( 0.190 -2.089 0.000) 2.097 7.556 ( -0.704 1.373 0.000) 1.543 7.566 ( -0.820 0.216 0.000) 0.848 7.862 ( -0.568 7.530 0.000) 7.551 7.990 ( 3.753 1.302 0.000) 3.972 8.053 ( -5.042 5.457 0.000) 7.430 8.241 ( -3.027 -2.493 0.000) 3.922 8.392 ( 0.089 -1.697 0.000) 1.700 8.439 ( -0.778 0.334 0.000) 0.846 8.665 ( 0.191 -3.574 0.000) 3.579 8.670 ( -0.099 -2.992 0.000) 2.993 ======================= Grid point 12 (9/24) ======================= q-point: ( 0.40 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.757 ( 7.869 5.461 0.000) 9.578 0.769 ( 8.421 4.269 0.000) 9.441 0.804 ( 5.250 4.896 0.000) 7.179 0.881 ( 0.228 4.825 0.000) 4.831 0.891 ( -3.101 -0.117 0.000) 3.103 0.929 ( 3.571 -0.597 0.000) 3.621 0.957 ( -5.594 3.093 0.000) 6.393 0.988 ( 0.684 0.507 0.000) 0.852 1.024 ( -2.686 2.458 0.000) 3.641 1.032 ( 0.109 4.738 0.000) 4.740 1.069 ( -1.637 2.552 0.000) 3.032 1.070 ( -2.912 0.723 0.000) 3.001 1.140 ( 2.494 0.575 0.000) 2.560 1.150 ( 1.902 2.936 0.000) 3.499 1.172 ( 0.426 -0.082 0.000) 0.434 1.198 ( -0.607 -0.584 0.000) 0.843 1.247 ( 2.235 4.374 0.000) 4.912 1.266 ( -0.017 2.838 0.000) 2.838 1.287 ( -1.551 2.817 0.000) 3.216 1.295 ( -3.299 2.219 0.000) 3.976 1.330 ( 1.897 5.840 0.000) 6.140 1.365 ( -1.172 -1.788 0.000) 2.138 1.404 ( -0.084 -2.676 0.000) 2.678 1.406 ( -0.003 -2.289 0.000) 2.289 1.420 ( 2.871 0.164 0.000) 2.876 1.444 ( 0.409 0.329 0.000) 0.525 1.461 ( -0.228 1.785 0.000) 1.799 1.505 ( -3.762 -0.379 0.000) 3.781 1.569 ( 1.568 -2.213 0.000) 2.712 1.613 ( -2.073 -1.309 0.000) 2.452 1.643 ( 1.104 2.127 0.000) 2.396 1.663 ( -0.924 3.260 0.000) 3.388 1.758 ( 3.726 2.626 0.000) 4.558 1.779 ( 3.295 4.701 0.000) 5.741 1.818 ( -1.684 0.842 0.000) 1.883 1.862 ( -3.669 4.138 0.000) 5.530 1.895 ( 1.253 -0.688 0.000) 1.430 1.905 ( -0.045 0.516 0.000) 0.518 1.983 ( 0.567 5.271 0.000) 5.302 2.016 ( -2.136 5.746 0.000) 6.130 2.141 ( 0.477 1.695 0.000) 1.761 2.142 ( 0.520 2.188 0.000) 2.249 2.253 ( 2.091 0.508 0.000) 2.152 2.286 ( -0.090 -0.607 0.000) 0.613 2.394 ( 0.846 0.390 0.000) 0.932 2.413 ( -0.179 -1.066 0.000) 1.081 2.433 ( -0.756 -1.854 0.000) 2.002 2.433 ( -1.293 -0.643 0.000) 1.444 2.478 ( -0.076 6.403 0.000) 6.403 2.482 ( -1.135 8.618 0.000) 8.693 2.595 ( -0.523 3.430 0.000) 3.470 2.599 ( -0.110 4.788 0.000) 4.789 2.853 ( 1.503 2.505 0.000) 2.921 2.876 ( -0.565 2.762 0.000) 2.819 2.914 ( -0.301 0.832 0.000) 0.885 2.935 ( -1.644 0.705 0.000) 1.789 3.007 ( 0.865 0.255 0.000) 0.902 3.031 ( -1.265 -0.152 0.000) 1.274 3.044 ( 0.848 -6.216 0.000) 6.273 3.055 ( -0.104 -4.337 0.000) 4.338 3.084 ( 0.020 2.683 0.000) 2.683 3.108 ( -2.212 1.898 0.000) 2.915 3.317 ( 0.832 -0.219 0.000) 0.860 3.325 ( -0.183 -1.501 0.000) 1.512 3.359 ( 1.308 -0.959 0.000) 1.622 3.394 ( -1.756 -0.990 0.000) 2.015 3.488 ( 2.801 -5.651 0.000) 6.307 3.547 ( -2.567 -6.073 0.000) 6.593 3.742 ( 1.381 1.117 0.000) 1.777 3.772 ( -1.201 0.940 0.000) 1.525 3.804 ( -0.284 3.417 0.000) 3.429 3.814 ( -0.758 3.258 0.000) 3.345 3.859 ( 0.796 -3.044 0.000) 3.147 3.862 ( 0.464 -2.677 0.000) 2.717 3.951 ( 1.036 -1.424 0.000) 1.760 3.978 ( -1.420 -1.484 0.000) 2.054 6.645 ( 1.153 -0.456 0.000) 1.240 6.661 ( -0.236 -0.630 0.000) 0.673 6.713 ( 0.472 -2.918 0.000) 2.955 6.740 ( -2.047 -2.671 0.000) 3.365 7.082 ( 1.113 0.489 0.000) 1.216 7.087 ( 0.048 1.486 0.000) 1.486 7.146 ( 2.839 -2.969 0.000) 4.108 7.222 ( -3.671 -4.533 0.000) 5.833 7.514 ( 1.865 -0.348 0.000) 1.897 7.540 ( -0.750 0.816 0.000) 1.108 7.547 ( -0.180 -0.630 0.000) 0.655 7.550 ( -0.577 -0.219 0.000) 0.617 7.876 ( 1.939 10.125 0.000) 10.309 7.945 ( -4.509 9.764 0.000) 10.755 8.074 ( 4.094 -1.209 0.000) 4.269 8.162 ( -4.110 -1.939 0.000) 4.545 8.399 ( 0.569 -2.633 0.000) 2.694 8.417 ( -1.121 -2.571 0.000) 2.805 8.669 ( 0.132 -3.688 0.000) 3.690 8.670 ( 0.024 -3.509 0.000) 3.509 ======================= Grid point 15 (10/24) ======================= q-point: ( 0.00 0.43 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.716 ( 0.000 6.794 0.000) 6.794 0.773 ( 0.000 7.181 0.000) 7.181 0.802 ( 0.000 6.411 0.000) 6.411 0.803 ( 0.000 2.214 0.000) 2.214 0.832 ( 0.000 3.013 0.000) 3.013 0.887 ( 0.000 2.088 0.000) 2.088 0.972 ( 0.000 1.562 0.000) 1.562 0.994 ( 0.000 -5.208 0.000) 5.208 1.002 ( 0.000 -8.190 0.000) 8.190 1.032 ( 0.000 -4.156 0.000) 4.156 1.080 ( 0.000 -0.539 0.000) 0.539 1.108 ( 0.000 -3.332 0.000) 3.332 1.159 ( 0.000 3.603 0.000) 3.603 1.203 ( 0.000 -1.081 0.000) 1.081 1.238 ( 0.000 1.980 0.000) 1.980 1.239 ( 0.000 0.872 0.000) 0.872 1.257 ( 0.000 2.682 0.000) 2.682 1.270 ( 0.000 1.336 0.000) 1.336 1.300 ( 0.000 -0.239 0.000) 0.239 1.309 ( 0.000 -2.142 0.000) 2.142 1.345 ( 0.000 -0.560 0.000) 0.560 1.358 ( 0.000 -2.127 0.000) 2.127 1.377 ( 0.000 3.780 0.000) 3.780 1.403 ( 0.000 9.579 0.000) 9.579 1.485 ( 0.000 -4.940 0.000) 4.940 1.507 ( 0.000 5.843 0.000) 5.843 1.534 ( 0.000 -3.488 0.000) 3.488 1.542 ( 0.000 4.045 0.000) 4.045 1.598 ( 0.000 1.548 0.000) 1.548 1.617 ( 0.000 0.260 0.000) 0.260 1.652 ( 0.000 -2.263 0.000) 2.263 1.666 ( 0.000 -0.599 0.000) 0.599 1.679 ( 0.000 -3.175 0.000) 3.175 1.687 ( 0.000 -2.836 0.000) 2.836 1.822 ( 0.000 1.893 0.000) 1.893 1.846 ( 0.000 -0.717 0.000) 0.717 1.942 ( 0.000 10.038 0.000) 10.038 2.045 ( 0.000 3.765 0.000) 3.765 2.090 ( 0.000 -3.892 0.000) 3.892 2.105 ( 0.000 -1.321 0.000) 1.321 2.134 ( 0.000 5.057 0.000) 5.057 2.185 ( 0.000 3.918 0.000) 3.918 2.218 ( 0.000 -2.493 0.000) 2.493 2.246 ( 0.000 -1.616 0.000) 1.616 2.409 ( 0.000 0.148 0.000) 0.148 2.427 ( 0.000 -1.082 0.000) 1.082 2.442 ( 0.000 -0.763 0.000) 0.763 2.449 ( 0.000 -0.922 0.000) 0.922 2.622 ( 0.000 6.444 0.000) 6.444 2.651 ( 0.000 3.749 0.000) 3.749 2.685 ( 0.000 6.925 0.000) 6.925 2.745 ( 0.000 4.833 0.000) 4.833 2.822 ( 0.000 -1.706 0.000) 1.706 2.847 ( 0.000 2.032 0.000) 2.032 2.881 ( 0.000 -3.178 0.000) 3.178 2.908 ( 0.000 -6.244 0.000) 6.244 2.967 ( 0.000 -0.291 0.000) 0.291 2.994 ( 0.000 -2.492 0.000) 2.492 3.041 ( 0.000 2.348 0.000) 2.348 3.060 ( 0.000 -0.428 0.000) 0.428 3.213 ( 0.000 1.355 0.000) 1.355 3.242 ( 0.000 3.710 0.000) 3.710 3.247 ( 0.000 -2.314 0.000) 2.314 3.259 ( 0.000 3.547 0.000) 3.547 3.339 ( 0.000 -5.325 0.000) 5.325 3.342 ( 0.000 -4.772 0.000) 4.772 3.450 ( 0.000 -1.337 0.000) 1.337 3.472 ( 0.000 -4.207 0.000) 4.207 3.712 ( 0.000 0.784 0.000) 0.784 3.722 ( 0.000 -0.493 0.000) 0.493 3.818 ( 0.000 1.705 0.000) 1.705 3.836 ( 0.000 -0.109 0.000) 0.109 3.884 ( 0.000 1.634 0.000) 1.634 3.912 ( 0.000 -1.282 0.000) 1.282 3.940 ( 0.000 3.340 0.000) 3.340 3.991 ( 0.000 -1.962 0.000) 1.962 6.642 ( 0.000 1.528 0.000) 1.528 6.674 ( 0.000 -1.763 0.000) 1.763 6.696 ( 0.000 2.568 0.000) 2.568 6.751 ( 0.000 -2.993 0.000) 2.993 7.022 ( 0.000 0.958 0.000) 0.958 7.052 ( 0.000 -2.315 0.000) 2.315 7.167 ( 0.000 2.887 0.000) 2.887 7.228 ( 0.000 -3.662 0.000) 3.662 7.424 ( 0.000 1.423 0.000) 1.423 7.474 ( 0.000 -3.500 0.000) 3.500 7.583 ( 0.000 0.290 0.000) 0.290 7.583 ( 0.000 0.371 0.000) 0.371 7.938 ( 0.000 2.485 0.000) 2.485 7.976 ( 0.000 -0.938 0.000) 0.938 8.151 ( 0.000 2.226 0.000) 2.226 8.205 ( 0.000 -3.345 0.000) 3.345 8.445 ( 0.000 5.996 0.000) 5.996 8.530 ( 0.000 4.380 0.000) 4.380 8.570 ( 0.000 -6.007 0.000) 6.007 8.609 ( 0.000 -3.761 0.000) 3.761 ======================= Grid point 16 (11/24) ======================= q-point: ( 0.20 0.43 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.742 ( 2.813 6.075 0.000) 6.695 0.813 ( 1.215 7.896 0.000) 7.989 0.831 ( 3.471 0.595 0.000) 3.522 0.863 ( 4.502 1.778 0.000) 4.841 0.865 ( 3.442 6.840 0.000) 7.657 0.913 ( 2.045 1.375 0.000) 2.464 0.971 ( 0.101 -0.273 0.000) 0.291 1.003 ( 1.337 -3.556 0.000) 3.799 1.006 ( 0.303 -4.437 0.000) 4.447 1.038 ( 0.923 -3.961 0.000) 4.068 1.075 ( 0.078 1.078 0.000) 1.081 1.121 ( 0.322 -3.700 0.000) 3.714 1.165 ( 0.076 1.459 0.000) 1.461 1.181 ( -1.794 -0.539 0.000) 1.873 1.198 ( -3.904 2.493 0.000) 4.632 1.228 ( -2.156 0.041 0.000) 2.156 1.264 ( 1.421 2.591 0.000) 2.955 1.270 ( 1.653 1.285 0.000) 2.094 1.312 ( 0.920 -1.764 0.000) 1.989 1.323 ( 1.453 -2.473 0.000) 2.868 1.340 ( -0.636 0.342 0.000) 0.722 1.364 ( 0.240 -2.213 0.000) 2.226 1.374 ( -0.352 4.122 0.000) 4.137 1.442 ( 1.271 7.208 0.000) 7.320 1.476 ( -0.299 -4.204 0.000) 4.215 1.476 ( -2.363 5.063 0.000) 5.587 1.521 ( -3.102 1.155 0.000) 3.310 1.522 ( -1.386 2.190 0.000) 2.592 1.528 ( -1.785 -1.590 0.000) 2.390 1.596 ( 1.448 -2.765 0.000) 3.122 1.633 ( -1.437 -0.847 0.000) 1.668 1.675 ( -0.453 -3.114 0.000) 3.147 1.729 ( 4.022 2.238 0.000) 4.602 1.777 ( 4.441 -2.315 0.000) 5.008 1.848 ( 1.345 1.693 0.000) 2.162 1.869 ( 1.521 -0.344 0.000) 1.559 1.941 ( 0.053 9.827 0.000) 9.827 2.031 ( -1.259 4.724 0.000) 4.889 2.096 ( 0.523 -4.773 0.000) 4.801 2.097 ( -1.847 0.254 0.000) 1.864 2.138 ( 1.336 3.161 0.000) 3.432 2.167 ( -1.631 4.596 0.000) 4.876 2.212 ( -0.098 -2.548 0.000) 2.549 2.241 ( -0.285 -2.538 0.000) 2.554 2.394 ( -0.759 0.303 0.000) 0.817 2.416 ( -0.998 -1.716 0.000) 1.985 2.428 ( -1.133 -0.049 0.000) 1.134 2.440 ( -0.576 -0.866 0.000) 1.040 2.621 ( 0.273 7.082 0.000) 7.088 2.674 ( -1.215 7.215 0.000) 7.317 2.677 ( 1.428 2.564 0.000) 2.935 2.719 ( -2.149 4.506 0.000) 4.992 2.842 ( 1.546 -1.888 0.000) 2.440 2.857 ( 0.623 1.365 0.000) 1.501 2.885 ( 0.198 -2.957 0.000) 2.964 2.910 ( 0.097 -6.799 0.000) 6.800 2.982 ( 0.608 -0.809 0.000) 1.012 2.983 ( -0.038 -0.774 0.000) 0.775 3.025 ( -1.190 2.309 0.000) 2.598 3.048 ( -0.973 -0.507 0.000) 1.098 3.191 ( -1.849 5.282 0.000) 5.596 3.227 ( 1.283 3.374 0.000) 3.610 3.281 ( 2.469 -2.259 0.000) 3.346 3.304 ( -0.146 -1.313 0.000) 1.321 3.323 ( 1.956 -3.588 0.000) 4.086 3.367 ( 0.684 -3.255 0.000) 3.326 3.368 ( -4.925 1.542 0.000) 5.160 3.458 ( -2.192 -6.863 0.000) 7.204 3.735 ( 1.751 1.173 0.000) 2.107 3.757 ( 2.495 -1.224 0.000) 2.779 3.814 ( -0.439 1.431 0.000) 1.496 3.831 ( -0.340 -0.445 0.000) 0.560 3.875 ( -0.654 2.116 0.000) 2.215 3.911 ( -0.015 -1.559 0.000) 1.559 3.924 ( -1.434 3.258 0.000) 3.560 3.975 ( -1.392 -2.090 0.000) 2.511 6.639 ( -0.106 1.179 0.000) 1.184 6.669 ( -0.381 -1.846 0.000) 1.885 6.685 ( -0.933 1.998 0.000) 2.206 6.732 ( -1.627 -2.733 0.000) 3.180 7.045 ( 1.961 0.629 0.000) 2.060 7.067 ( 1.295 -1.824 0.000) 2.237 7.142 ( -1.680 2.497 0.000) 3.009 7.200 ( -2.206 -3.696 0.000) 4.304 7.461 ( 2.777 1.111 0.000) 2.991 7.499 ( 1.776 -2.551 0.000) 3.108 7.571 ( -0.971 0.313 0.000) 1.020 7.572 ( -0.945 0.285 0.000) 0.987 7.955 ( 1.567 2.624 0.000) 3.056 7.992 ( 1.528 -1.247 0.000) 1.973 8.127 ( -2.191 2.675 0.000) 3.458 8.181 ( -2.277 -3.173 0.000) 3.906 8.441 ( -0.363 6.049 0.000) 6.060 8.500 ( -2.355 4.941 0.000) 5.473 8.570 ( -0.055 -6.240 0.000) 6.240 8.597 ( -0.987 -4.590 0.000) 4.695 ======================= Grid point 17 (12/24) ======================= q-point: ( 0.40 0.43 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 36 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.822 ( 4.074 2.943 0.000) 5.026 0.856 ( 2.113 4.132 0.000) 4.641 0.880 ( 1.312 -2.637 0.000) 2.946 0.901 ( -1.392 0.011 0.000) 1.392 0.956 ( 3.714 6.968 0.000) 7.896 0.966 ( 1.864 2.851 0.000) 3.406 0.984 ( 1.663 0.061 0.000) 1.664 1.008 ( -0.774 1.973 0.000) 2.120 1.057 ( 2.140 0.293 0.000) 2.160 1.078 ( 1.721 -1.030 0.000) 2.006 1.080 ( 0.204 0.815 0.000) 0.840 1.085 ( -2.518 0.810 0.000) 2.645 1.118 ( -1.538 1.150 0.000) 1.921 1.147 ( -1.226 -2.906 0.000) 3.154 1.168 ( -0.171 -2.991 0.000) 2.996 1.179 ( -1.829 -2.309 0.000) 2.945 1.294 ( 1.327 1.876 0.000) 2.298 1.315 ( 1.298 1.249 0.000) 1.802 1.320 ( -0.293 1.551 0.000) 1.578 1.339 ( 0.444 0.384 0.000) 0.587 1.339 ( -1.067 -1.119 0.000) 1.546 1.353 ( -0.099 1.745 0.000) 1.748 1.362 ( -0.317 -2.373 0.000) 2.394 1.387 ( -3.487 1.041 0.000) 3.640 1.447 ( 0.844 3.297 0.000) 3.403 1.481 ( -0.750 -1.492 0.000) 1.670 1.503 ( 0.711 -2.706 0.000) 2.798 1.510 ( 0.275 -0.922 0.000) 0.962 1.547 ( 2.206 5.784 0.000) 6.190 1.594 ( -2.237 1.729 0.000) 2.828 1.631 ( 1.330 -2.647 0.000) 2.962 1.658 ( -1.164 -3.060 0.000) 3.274 1.802 ( 2.822 2.558 0.000) 3.808 1.847 ( -1.333 2.131 0.000) 2.514 1.853 ( 2.613 -2.103 0.000) 3.354 1.885 ( -0.362 -1.276 0.000) 1.326 1.954 ( 1.360 9.111 0.000) 9.212 1.995 ( -2.094 6.740 0.000) 7.058 2.072 ( 0.086 1.396 0.000) 1.398 2.094 ( -1.292 -2.950 0.000) 3.221 2.143 ( 0.193 3.504 0.000) 3.509 2.153 ( -0.143 1.042 0.000) 1.052 2.223 ( 0.809 -3.487 0.000) 3.580 2.236 ( -0.337 -4.309 0.000) 4.322 2.395 ( 0.360 -0.055 0.000) 0.364 2.402 ( 0.582 -0.733 0.000) 0.936 2.406 ( -0.822 0.858 0.000) 1.188 2.426 ( -0.911 -0.761 0.000) 1.187 2.629 ( 0.437 7.023 0.000) 7.037 2.646 ( -1.071 6.633 0.000) 6.719 2.674 ( -1.245 3.419 0.000) 3.639 2.682 ( -0.965 2.071 0.000) 2.285 2.866 ( 0.304 -0.428 0.000) 0.525 2.867 ( 0.478 -1.470 0.000) 1.546 2.891 ( 0.585 -3.375 0.000) 3.426 2.907 ( -0.560 -5.298 0.000) 5.327 2.983 ( 0.655 1.255 0.000) 1.416 3.001 ( -0.969 2.076 0.000) 2.291 3.009 ( 1.091 -0.688 0.000) 1.290 3.028 ( -0.834 -0.468 0.000) 0.956 3.183 ( 1.122 6.651 0.000) 6.745 3.226 ( -2.145 7.263 0.000) 7.573 3.300 ( 0.161 -2.620 0.000) 2.625 3.304 ( -0.460 -0.441 0.000) 0.637 3.331 ( 2.036 -1.308 0.000) 2.420 3.363 ( -0.418 -2.236 0.000) 2.275 3.370 ( 1.259 -4.212 0.000) 4.396 3.407 ( -2.139 -7.833 0.000) 8.119 3.770 ( 1.433 1.658 0.000) 2.191 3.797 ( -1.091 1.466 0.000) 1.827 3.803 ( 1.622 -1.801 0.000) 2.424 3.824 ( -0.468 -1.426 0.000) 1.501 3.867 ( 0.188 2.953 0.000) 2.959 3.888 ( -1.614 3.129 0.000) 3.520 3.917 ( 0.654 -2.131 0.000) 2.229 3.941 ( -1.590 -2.312 0.000) 2.806 6.644 ( 0.632 0.514 0.000) 0.815 6.663 ( -1.053 0.832 0.000) 1.342 6.667 ( 0.420 -1.804 0.000) 1.853 6.692 ( -1.840 -2.116 0.000) 2.804 7.087 ( 1.514 0.333 0.000) 1.550 7.103 ( 2.158 -1.452 0.000) 2.601 7.112 ( -1.116 1.437 0.000) 1.820 7.151 ( -2.279 -2.756 0.000) 3.576 7.515 ( 2.008 0.470 0.000) 2.062 7.530 ( 1.069 -0.967 0.000) 1.441 7.548 ( -1.240 0.157 0.000) 1.250 7.550 ( -0.975 0.062 0.000) 0.977 8.001 ( 2.712 2.990 0.000) 4.036 8.041 ( 3.067 -1.917 0.000) 3.617 8.066 ( -3.247 3.102 0.000) 4.490 8.115 ( -3.610 -2.576 0.000) 4.435 8.436 ( 0.130 6.071 0.000) 6.073 8.454 ( -1.621 5.747 0.000) 5.972 8.570 ( 0.138 -6.510 0.000) 6.512 8.578 ( -0.655 -5.996 0.000) 6.031 ======================= Grid point 35 (13/24) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 24 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.310 ( 0.000 0.000 13.289) 13.289 0.315 ( 0.000 0.000 13.498) 13.498 0.526 ( 0.000 0.000 7.143) 7.143 0.559 ( 0.000 0.000 -8.084) 8.084 0.577 ( 0.000 0.000 -6.863) 6.863 0.583 ( 0.000 0.000 21.718) 21.718 0.707 ( 0.000 0.000 -6.122) 6.122 0.831 ( 0.000 0.000 3.420) 3.420 0.946 ( 0.000 0.000 -5.659) 5.659 0.960 ( 0.000 0.000 -3.721) 3.721 1.042 ( 0.000 0.000 2.300) 2.300 1.099 ( 0.000 0.000 -2.571) 2.571 1.110 ( 0.000 0.000 -3.570) 3.570 1.125 ( 0.000 0.000 1.532) 1.532 1.167 ( 0.000 0.000 -8.752) 8.752 1.167 ( 0.000 0.000 -1.474) 1.474 1.207 ( 0.000 0.000 2.258) 2.258 1.234 ( 0.000 0.000 1.814) 1.814 1.299 ( 0.000 0.000 -3.224) 3.224 1.310 ( 0.000 0.000 -0.313) 0.313 1.326 ( 0.000 0.000 0.641) 0.641 1.327 ( 0.000 0.000 -1.680) 1.680 1.416 ( 0.000 0.000 0.196) 0.196 1.430 ( 0.000 0.000 1.876) 1.876 1.452 ( 0.000 0.000 -1.943) 1.943 1.472 ( 0.000 0.000 -0.448) 0.448 1.512 ( 0.000 0.000 3.532) 3.532 1.538 ( 0.000 0.000 2.139) 2.139 1.569 ( 0.000 0.000 -0.089) 0.089 1.571 ( 0.000 0.000 0.307) 0.307 1.581 ( 0.000 0.000 -1.115) 1.115 1.595 ( 0.000 0.000 -2.917) 2.917 1.636 ( 0.000 0.000 2.338) 2.338 1.640 ( 0.000 0.000 0.687) 0.687 1.674 ( 0.000 0.000 -1.347) 1.347 1.702 ( 0.000 0.000 -3.603) 3.603 1.808 ( 0.000 0.000 3.296) 3.296 1.885 ( 0.000 0.000 -2.641) 2.641 2.038 ( 0.000 0.000 3.598) 3.598 2.042 ( 0.000 0.000 7.075) 7.075 2.074 ( 0.000 0.000 -0.745) 0.745 2.078 ( 0.000 0.000 2.162) 2.162 2.115 ( 0.000 0.000 0.615) 0.615 2.174 ( 0.000 0.000 -4.500) 4.500 2.269 ( 0.000 0.000 -1.535) 1.535 2.303 ( 0.000 0.000 -4.202) 4.202 2.474 ( 0.000 0.000 -0.006) 0.006 2.479 ( 0.000 0.000 -0.285) 0.285 2.514 ( 0.000 0.000 3.717) 3.717 2.545 ( 0.000 0.000 0.325) 0.325 2.561 ( 0.000 0.000 -1.036) 1.036 2.566 ( 0.000 0.000 -0.861) 0.861 2.723 ( 0.000 0.000 0.961) 0.961 2.736 ( 0.000 0.000 -0.135) 0.135 2.818 ( 0.000 0.000 1.329) 1.329 2.862 ( 0.000 0.000 -2.110) 2.110 3.015 ( 0.000 0.000 1.761) 1.761 3.038 ( 0.000 0.000 2.647) 2.647 3.067 ( 0.000 0.000 -2.678) 2.678 3.101 ( 0.000 0.000 -2.580) 2.580 3.151 ( 0.000 0.000 0.761) 0.761 3.156 ( 0.000 0.000 0.159) 0.159 3.322 ( 0.000 0.000 0.883) 0.883 3.343 ( 0.000 0.000 -0.838) 0.838 3.388 ( 0.000 0.000 2.580) 2.580 3.446 ( 0.000 0.000 -2.231) 2.231 3.679 ( 0.000 0.000 1.135) 1.135 3.681 ( 0.000 0.000 0.299) 0.299 3.683 ( 0.000 0.000 0.330) 0.330 3.704 ( 0.000 0.000 -0.910) 0.910 3.790 ( 0.000 0.000 0.696) 0.696 3.820 ( 0.000 0.000 -1.655) 1.655 3.832 ( 0.000 0.000 1.140) 1.140 3.855 ( 0.000 0.000 -0.780) 0.780 4.027 ( 0.000 0.000 0.264) 0.264 4.033 ( 0.000 0.000 -0.237) 0.237 6.614 ( 0.000 0.000 0.162) 0.162 6.617 ( 0.000 0.000 -0.110) 0.110 6.819 ( 0.000 0.000 0.022) 0.022 6.820 ( 0.000 0.000 -0.021) 0.021 6.882 ( 0.000 0.000 0.279) 0.279 6.888 ( 0.000 0.000 -0.261) 0.261 7.457 ( 0.000 0.000 1.033) 1.033 7.472 ( 0.000 0.000 0.086) 0.086 7.474 ( 0.000 0.000 -0.084) 0.084 7.484 ( 0.000 0.000 -1.160) 1.160 7.568 ( 0.000 0.000 0.121) 0.121 7.571 ( 0.000 0.000 -0.123) 0.123 7.776 ( 0.000 0.000 0.231) 0.231 7.782 ( 0.000 0.000 -0.265) 0.265 8.135 ( 0.000 0.000 2.366) 2.366 8.191 ( 0.000 0.000 -2.249) 2.249 8.442 ( 0.000 0.000 1.458) 1.458 8.476 ( 0.000 0.000 -1.382) 1.382 8.724 ( 0.000 0.000 0.004) 0.004 8.724 ( 0.000 0.000 -0.037) 0.037 ======================= Grid point 36 (14/24) ======================= q-point: ( 0.20 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 32 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.404 ( 7.517 0.000 10.454) 12.876 0.420 ( 8.047 0.000 8.283) 11.548 0.590 ( 5.639 0.000 4.372) 7.135 0.591 ( 1.968 0.000 -5.317) 5.670 0.653 ( 5.836 0.000 -9.117) 10.825 0.727 ( 4.015 0.000 -2.340) 4.647 0.760 ( 11.366 0.000 12.734) 17.069 0.824 ( 0.124 0.000 1.696) 1.701 0.914 ( -2.758 0.000 -4.764) 5.505 1.005 ( 2.018 0.000 -1.959) 2.812 1.013 ( -2.171 0.000 1.227) 2.494 1.075 ( -2.467 0.000 -3.290) 4.113 1.104 ( -0.661 0.000 -1.637) 1.766 1.108 ( -1.406 0.000 0.608) 1.532 1.138 ( -2.957 0.000 -1.703) 3.413 1.164 ( 0.567 0.000 -4.516) 4.551 1.203 ( -0.106 0.000 3.142) 3.144 1.233 ( -0.285 0.000 2.110) 2.129 1.288 ( -1.645 0.000 -1.173) 2.020 1.290 ( -2.558 0.000 -0.003) 2.558 1.304 ( -0.802 0.000 -3.444) 3.536 1.343 ( 1.627 0.000 -3.642) 3.989 1.418 ( 0.683 0.000 0.898) 1.128 1.423 ( -2.122 0.000 2.736) 3.463 1.451 ( 0.274 0.000 1.213) 1.243 1.453 ( 0.182 0.000 -1.145) 1.159 1.524 ( -1.029 0.000 -1.158) 1.549 1.535 ( -2.232 0.000 -0.912) 2.411 1.563 ( 0.367 0.000 0.676) 0.770 1.583 ( 2.717 0.000 2.917) 3.986 1.598 ( 0.966 0.000 -2.310) 2.504 1.631 ( 0.623 0.000 0.278) 0.682 1.658 ( 0.922 0.000 0.803) 1.223 1.674 ( 1.816 0.000 1.745) 2.518 1.707 ( 5.352 0.000 -3.420) 6.351 1.715 ( 0.805 0.000 -1.736) 1.914 1.814 ( 0.493 0.000 1.899) 1.962 1.869 ( -1.382 0.000 -2.453) 2.816 2.006 ( -2.941 0.000 4.293) 5.203 2.022 ( -1.936 0.000 3.138) 3.687 2.065 ( 0.192 0.000 0.923) 0.943 2.075 ( -6.479 0.000 -0.693) 6.516 2.104 ( 1.690 0.000 -0.686) 1.823 2.143 ( 3.673 0.000 -2.270) 4.318 2.296 ( 2.660 0.000 0.717) 2.755 2.338 ( 1.380 0.000 -2.207) 2.603 2.471 ( -0.671 0.000 0.462) 0.815 2.479 ( -0.222 0.000 0.108) 0.247 2.493 ( -1.502 0.000 3.005) 3.360 2.533 ( -1.235 0.000 0.766) 1.453 2.546 ( -1.786 0.000 -1.364) 2.247 2.546 ( -1.165 0.000 -0.492) 1.264 2.736 ( 1.112 0.000 0.429) 1.192 2.744 ( 0.757 0.000 -0.148) 0.771 2.819 ( 0.162 0.000 1.556) 1.564 2.862 ( 0.020 0.000 -1.609) 1.609 3.015 ( -1.771 0.000 2.222) 2.842 3.027 ( 1.312 0.000 1.328) 1.867 3.061 ( -0.139 0.000 -1.435) 1.441 3.076 ( -2.139 0.000 -2.688) 3.436 3.216 ( 4.196 0.000 0.044) 4.196 3.236 ( 4.971 0.000 -1.358) 5.153 3.329 ( 0.543 0.000 0.860) 1.017 3.346 ( 0.210 0.000 -0.568) 0.606 3.349 ( -1.748 0.000 2.717) 3.230 3.409 ( -2.863 0.000 -2.396) 3.733 3.661 ( -2.148 0.000 1.286) 2.503 3.687 ( 0.644 0.000 0.873) 1.085 3.692 ( -0.777 0.000 -1.152) 1.389 3.706 ( 0.192 0.000 -0.742) 0.766 3.742 ( -1.650 0.000 1.139) 2.005 3.781 ( -1.519 0.000 -1.753) 2.319 3.896 ( 3.146 0.000 -0.256) 3.157 3.905 ( 2.811 0.000 -1.080) 3.012 4.023 ( -0.595 0.000 0.539) 0.803 4.025 ( -0.397 0.000 0.203) 0.446 6.628 ( 1.600 0.000 0.882) 1.827 6.667 ( 1.909 0.000 -3.189) 3.717 6.820 ( 0.080 0.000 0.019) 0.083 6.820 ( 0.074 0.000 -0.019) 0.076 6.931 ( 4.208 0.000 1.030) 4.332 6.953 ( 5.260 0.000 -0.603) 5.294 7.378 ( -6.374 0.000 5.555) 8.455 7.430 ( -2.814 0.000 -0.315) 2.832 7.478 ( 0.484 0.000 0.091) 0.493 7.480 ( 0.486 0.000 -0.090) 0.495 7.562 ( -0.460 0.000 0.075) 0.466 7.564 ( -0.552 0.000 -0.076) 0.558 7.660 ( -5.761 0.000 3.091) 6.537 7.769 ( -4.937 0.000 -5.655) 7.506 8.060 ( 0.339 0.000 6.129) 6.139 8.192 ( -1.442 0.000 -5.962) 6.134 8.507 ( 2.798 0.000 0.703) 2.885 8.529 ( 2.565 0.000 -0.822) 2.693 8.724 ( 0.177 0.000 0.093) 0.200 8.727 ( 0.151 0.000 -0.141) 0.206 ======================= Grid point 37 (15/24) ======================= q-point: ( 0.40 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 32 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.559 ( 3.780 0.000 1.952) 4.254 0.585 ( 8.566 0.000 1.971) 8.790 0.600 ( -0.928 0.000 -1.082) 1.425 0.717 ( 0.348 0.000 -2.610) 2.633 0.736 ( 2.792 0.000 -5.193) 5.896 0.777 ( 3.015 0.000 -3.288) 4.461 0.859 ( 2.399 0.000 -0.241) 2.411 0.882 ( 0.699 0.000 -1.439) 1.600 0.949 ( -1.491 0.000 1.827) 2.358 0.993 ( -3.188 0.000 -0.875) 3.306 1.037 ( 3.990 0.000 2.336) 4.624 1.054 ( -3.400 0.000 -1.309) 3.643 1.056 ( 6.886 0.000 0.932) 6.948 1.063 ( -2.946 0.000 -1.978) 3.549 1.090 ( 0.023 0.000 2.119) 2.120 1.185 ( -0.482 0.000 -3.553) 3.586 1.203 ( 0.117 0.000 3.313) 3.315 1.214 ( -1.051 0.000 3.110) 3.283 1.261 ( -0.067 0.000 -0.827) 0.830 1.267 ( -0.489 0.000 -2.099) 2.155 1.279 ( 0.697 0.000 3.843) 3.906 1.374 ( -4.296 0.000 5.615) 7.070 1.387 ( 2.280 0.000 -3.082) 3.834 1.433 ( -2.322 0.000 1.300) 2.661 1.441 ( 1.292 0.000 1.187) 1.754 1.458 ( -0.037 0.000 0.253) 0.256 1.467 ( -1.034 0.000 -3.019) 3.191 1.488 ( -2.443 0.000 -1.530) 2.883 1.582 ( 1.113 0.000 0.941) 1.457 1.604 ( -0.638 0.000 -0.232) 0.679 1.609 ( 0.864 0.000 0.666) 1.091 1.630 ( -0.687 0.000 0.090) 0.693 1.693 ( 0.878 0.000 -0.203) 0.901 1.718 ( -0.864 0.000 -1.572) 1.794 1.742 ( 5.214 0.000 2.251) 5.679 1.803 ( 0.535 0.000 -0.925) 1.068 1.837 ( 3.540 0.000 -1.134) 3.718 1.850 ( 0.806 0.000 -1.305) 1.533 1.932 ( -3.966 0.000 -0.314) 3.979 1.938 ( -6.055 0.000 -0.539) 6.079 2.010 ( 0.538 0.000 0.165) 0.562 2.033 ( -1.623 0.000 -1.537) 2.235 2.187 ( 4.946 0.000 2.723) 5.647 2.242 ( 3.271 0.000 -1.689) 3.681 2.312 ( -1.539 0.000 1.817) 2.382 2.339 ( -0.867 0.000 -0.496) 0.999 2.458 ( -0.219 0.000 1.145) 1.166 2.466 ( -0.821 0.000 0.920) 1.233 2.480 ( 0.430 0.000 0.065) 0.435 2.503 ( -1.764 0.000 -1.597) 2.380 2.512 ( -0.341 0.000 0.055) 0.345 2.520 ( -1.017 0.000 -0.443) 1.109 2.757 ( 0.569 0.000 -0.028) 0.570 2.758 ( 0.404 0.000 -0.113) 0.419 2.830 ( 1.041 0.000 2.129) 2.370 2.856 ( -0.890 0.000 0.812) 1.205 2.990 ( 0.053 0.000 -0.395) 0.398 3.022 ( -2.553 0.000 -2.525) 3.591 3.057 ( 1.129 0.000 0.890) 1.437 3.066 ( 0.223 0.000 0.249) 0.334 3.273 ( 1.319 0.000 0.040) 1.320 3.298 ( -0.299 0.000 -1.915) 1.938 3.341 ( 0.539 0.000 0.615) 0.818 3.348 ( -0.089 0.000 0.049) 0.101 3.357 ( 1.865 0.000 0.591) 1.957 3.369 ( 0.132 0.000 -0.406) 0.427 3.612 ( -1.297 0.000 0.445) 1.371 3.636 ( -3.080 0.000 -1.500) 3.425 3.701 ( 0.576 0.000 0.293) 0.646 3.708 ( -0.093 0.000 -0.302) 0.316 3.743 ( 0.785 0.000 1.038) 1.302 3.760 ( -0.820 0.000 -0.409) 0.916 3.949 ( 1.715 0.000 -0.347) 1.750 3.952 ( 1.423 0.000 -0.775) 1.620 4.008 ( -0.601 0.000 0.634) 0.874 4.012 ( -0.769 0.000 0.215) 0.799 6.662 ( 1.390 0.000 0.215) 1.406 6.681 ( -0.301 0.000 -1.564) 1.593 6.822 ( 0.057 0.000 0.009) 0.057 6.822 ( 0.039 0.000 -0.008) 0.040 7.031 ( 4.430 0.000 2.578) 5.126 7.062 ( 3.116 0.000 0.738) 3.203 7.281 ( -1.303 0.000 3.401) 3.642 7.346 ( -5.018 0.000 -1.611) 5.271 7.489 ( 0.464 0.000 0.054) 0.468 7.491 ( 0.364 0.000 -0.058) 0.368 7.552 ( -0.408 0.000 0.008) 0.408 7.552 ( -0.446 0.000 -0.005) 0.446 7.588 ( -0.777 0.000 0.298) 0.832 7.635 ( -5.251 0.000 -3.343) 6.225 8.114 ( 3.108 0.000 2.271) 3.849 8.164 ( -1.587 0.000 -1.970) 2.530 8.541 ( 0.900 0.000 0.327) 0.958 8.553 ( -0.060 0.000 -0.615) 0.618 8.729 ( 0.170 0.000 0.026) 0.172 8.730 ( 0.079 0.000 -0.051) 0.094 ======================= Grid point 40 (16/24) ======================= q-point: ( 0.00 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 33 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.421 ( 0.000 9.677 9.234) 13.376 0.460 ( 0.000 8.412 11.510) 14.256 0.551 ( 0.000 2.581 6.381) 6.884 0.623 ( 0.000 4.644 -7.968) 9.223 0.702 ( 0.000 8.990 -9.155) 12.831 0.710 ( 0.000 1.082 -5.475) 5.581 0.763 ( 0.000 16.866 14.648) 22.339 0.917 ( 0.000 7.792 2.787) 8.275 0.962 ( 0.000 1.924 -0.495) 1.987 1.014 ( 0.000 5.224 -4.817) 7.106 1.037 ( 0.000 -0.600 2.574) 2.643 1.094 ( 0.000 -0.402 -2.189) 2.225 1.109 ( 0.000 0.508 -3.737) 3.772 1.149 ( 0.000 1.228 1.028) 1.602 1.157 ( 0.000 0.055 -3.018) 3.018 1.181 ( 0.000 -0.134 -2.129) 2.134 1.223 ( 0.000 -0.580 1.878) 1.966 1.238 ( 0.000 2.590 1.206) 2.857 1.313 ( 0.000 -0.627 0.277) 0.686 1.315 ( 0.000 -1.152 0.013) 1.152 1.333 ( 0.000 3.760 -6.384) 7.408 1.342 ( 0.000 1.588 -2.210) 2.722 1.404 ( 0.000 -2.029 0.796) 2.179 1.416 ( 0.000 -4.448 1.172) 4.600 1.454 ( 0.000 1.772 2.148) 2.784 1.479 ( 0.000 0.940 -0.122) 0.948 1.527 ( 0.000 0.428 2.118) 2.161 1.536 ( 0.000 1.458 2.705) 3.073 1.561 ( 0.000 1.614 1.894) 2.489 1.583 ( 0.000 -1.905 -0.263) 1.923 1.610 ( 0.000 -1.345 -1.986) 2.399 1.614 ( 0.000 1.806 -2.726) 3.270 1.667 ( 0.000 0.488 1.246) 1.339 1.669 ( 0.000 4.408 0.371) 4.423 1.698 ( 0.000 3.961 -4.373) 5.900 1.734 ( 0.000 3.341 -4.655) 5.730 1.834 ( 0.000 3.292 2.529) 4.152 1.891 ( 0.000 0.967 -2.004) 2.225 2.018 ( 0.000 0.808 2.625) 2.747 2.038 ( 0.000 -3.763 1.280) 3.975 2.097 ( 0.000 3.285 3.632) 4.897 2.123 ( 0.000 4.120 -0.196) 4.125 2.123 ( 0.000 4.017 1.201) 4.193 2.159 ( 0.000 0.296 0.199) 0.356 2.321 ( 0.000 4.047 -1.332) 4.261 2.340 ( 0.000 3.567 -3.339) 4.886 2.455 ( 0.000 -0.986 0.666) 1.190 2.461 ( 0.000 -0.702 -0.118) 0.712 2.518 ( 0.000 -0.576 1.521) 1.626 2.541 ( 0.000 0.372 1.899) 1.935 2.544 ( 0.000 -0.965 -1.167) 1.514 2.586 ( 0.000 0.883 -1.901) 2.096 2.746 ( 0.000 2.351 1.176) 2.628 2.768 ( 0.000 3.205 -1.002) 3.358 2.912 ( 0.000 2.898 -1.050) 3.082 2.933 ( 0.000 2.437 -3.493) 4.259 3.036 ( 0.000 -0.082 2.141) 2.143 3.046 ( 0.000 2.004 -0.274) 2.022 3.067 ( 0.000 -0.328 -1.650) 1.682 3.093 ( 0.000 -0.959 -2.020) 2.236 3.152 ( 0.000 0.880 0.787) 1.181 3.164 ( 0.000 1.574 -0.331) 1.609 3.263 ( 0.000 -3.050 1.931) 3.610 3.319 ( 0.000 -0.522 -3.135) 3.178 3.431 ( 0.000 -2.445 0.712) 2.547 3.434 ( 0.000 1.779 0.445) 1.834 3.581 ( 0.000 -4.731 1.862) 5.084 3.634 ( 0.000 -2.245 -2.225) 3.161 3.748 ( 0.000 -0.465 2.449) 2.493 3.765 ( 0.000 3.115 0.645) 3.181 3.779 ( 0.000 1.673 -0.451) 1.732 3.803 ( 0.000 -0.948 -2.091) 2.296 3.810 ( 0.000 -1.868 1.274) 2.261 3.838 ( 0.000 -1.262 -0.951) 1.580 4.002 ( 0.000 -1.918 1.927) 2.718 4.022 ( 0.000 -0.494 -0.254) 0.556 6.611 ( 0.000 0.184 0.546) 0.576 6.622 ( 0.000 0.588 -0.405) 0.714 6.767 ( 0.000 -2.852 0.456) 2.889 6.794 ( 0.000 -1.454 -2.535) 2.922 6.960 ( 0.000 4.132 0.353) 4.147 6.975 ( 0.000 5.394 -1.450) 5.585 7.355 ( 0.000 -6.866 3.882) 7.888 7.406 ( 0.000 -5.319 -1.540) 5.538 7.491 ( 0.000 -1.839 3.932) 4.340 7.525 ( 0.000 2.738 0.830) 2.862 7.556 ( 0.000 0.772 -1.186) 1.415 7.561 ( 0.000 -1.279 -1.714) 2.139 7.815 ( 0.000 4.150 -0.080) 4.150 7.817 ( 0.000 3.374 -0.284) 3.386 8.141 ( 0.000 0.850 3.001) 3.118 8.226 ( 0.000 2.555 -4.587) 5.251 8.419 ( 0.000 -1.756 0.871) 1.961 8.442 ( 0.000 -3.077 -0.940) 3.217 8.711 ( 0.000 -1.365 0.078) 1.367 8.712 ( 0.000 -1.338 -0.010) 1.338 ======================= Grid point 41 (17/24) ======================= q-point: ( 0.20 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 53 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.481 ( 5.526 7.452 8.082) 12.304 0.554 ( 7.242 7.649 7.613) 12.997 0.618 ( 5.830 2.590 3.869) 7.461 0.680 ( 4.244 3.773 -7.326) 9.269 0.728 ( 2.979 5.650 -6.197) 8.899 0.756 ( 3.066 4.596 -5.993) 8.151 0.862 ( 6.279 11.863 8.315) 15.789 0.910 ( 0.802 7.946 2.522) 8.375 0.973 ( -1.256 3.809 -2.568) 4.763 0.999 ( 0.513 1.589 1.747) 2.417 1.021 ( -0.872 1.598 0.023) 1.821 1.074 ( -2.186 -0.146 -2.696) 3.473 1.090 ( -1.320 -0.392 -1.378) 1.948 1.131 ( -1.344 0.649 -0.382) 1.541 1.141 ( -2.088 1.046 -1.449) 2.748 1.162 ( 0.111 -0.394 -0.170) 0.443 1.224 ( 0.008 -0.289 1.565) 1.592 1.243 ( 0.326 2.801 0.636) 2.891 1.290 ( -1.917 0.285 -0.040) 1.938 1.299 ( -1.849 0.583 -0.269) 1.957 1.322 ( -0.639 1.320 -1.562) 2.143 1.370 ( 2.058 3.036 -4.373) 5.707 1.398 ( -1.059 -2.137 1.632) 2.890 1.422 ( 0.014 -2.007 -0.555) 2.082 1.463 ( 0.446 2.535 1.391) 2.925 1.488 ( 0.605 1.780 -0.397) 1.921 1.501 ( -1.174 -0.944 0.081) 1.508 1.522 ( -2.387 -0.427 -0.920) 2.594 1.547 ( 0.589 -2.047 1.985) 2.912 1.590 ( -0.628 -1.587 -0.281) 1.730 1.617 ( 1.964 2.184 3.054) 4.237 1.652 ( 1.219 0.669 1.099) 1.772 1.656 ( 0.837 1.826 0.340) 2.038 1.686 ( 1.240 2.659 -0.243) 2.944 1.733 ( 2.495 1.607 -3.426) 4.533 1.753 ( 1.542 3.386 -5.013) 6.243 1.854 ( 1.265 3.606 2.000) 4.314 1.896 ( 0.024 2.502 -1.679) 3.013 1.978 ( -2.787 -0.071 2.315) 3.624 2.010 ( -2.709 -1.128 0.742) 3.027 2.088 ( -1.484 2.962 1.942) 3.840 2.108 ( -1.872 2.089 0.793) 2.915 2.149 ( 2.360 4.123 0.208) 4.755 2.181 ( 3.616 3.258 -1.549) 5.107 2.334 ( 1.184 3.103 -0.506) 3.360 2.355 ( 0.499 1.499 -2.145) 2.664 2.453 ( -0.324 -1.034 1.209) 1.624 2.465 ( -0.141 -0.693 0.175) 0.729 2.499 ( -1.379 0.350 1.791) 2.287 2.518 ( -1.404 -0.118 1.440) 2.015 2.534 ( -1.033 -0.472 -0.819) 1.400 2.567 ( -1.712 0.978 -1.474) 2.462 2.765 ( 1.472 2.744 0.689) 3.189 2.781 ( 0.854 3.487 -0.787) 3.675 2.910 ( -1.200 4.413 -1.235) 4.737 2.912 ( -0.257 3.601 -0.871) 3.714 3.014 ( -2.108 0.055 2.006) 2.910 3.033 ( 0.269 1.180 1.264) 1.750 3.063 ( -0.995 0.291 -1.664) 1.961 3.069 ( -1.815 -0.375 -1.777) 2.568 3.164 ( 1.438 -5.081 0.828) 5.345 3.193 ( 1.483 -3.465 -1.625) 4.104 3.288 ( 2.153 -0.938 1.835) 2.980 3.327 ( 0.786 -0.359 -1.398) 1.643 3.394 ( -1.656 0.085 1.053) 1.964 3.418 ( -1.737 0.576 -0.587) 1.922 3.582 ( -0.063 -4.592 2.224) 5.102 3.645 ( -0.315 -2.858 -3.135) 4.254 3.729 ( -0.515 -0.535 1.774) 1.923 3.749 ( -1.141 2.114 0.195) 2.410 3.762 ( -0.868 3.317 -0.064) 3.429 3.781 ( -0.940 0.685 -1.552) 1.939 3.867 ( 3.121 -2.622 -0.093) 4.078 3.878 ( 2.441 -2.230 -0.792) 3.400 3.999 ( -0.275 -1.757 1.673) 2.442 4.017 ( -0.480 -0.723 -0.226) 0.897 6.627 ( 1.365 -0.041 0.702) 1.535 6.650 ( 1.680 -0.873 -1.334) 2.316 6.772 ( 0.379 -3.118 0.405) 3.167 6.792 ( -0.263 -1.885 -1.549) 2.454 6.996 ( 2.999 3.747 0.360) 4.813 7.008 ( 3.098 3.707 -0.697) 4.881 7.302 ( -3.784 -5.549 5.377) 8.604 7.375 ( -2.594 -4.892 -1.581) 5.759 7.501 ( 0.557 0.036 1.092) 1.227 7.531 ( 0.013 2.240 -0.375) 2.271 7.540 ( -0.386 -0.092 0.154) 0.426 7.559 ( -0.485 -0.380 -0.693) 0.928 7.734 ( -4.014 6.716 2.635) 8.256 7.821 ( -3.582 5.347 -4.977) 8.136 8.079 ( -0.768 1.932 5.803) 6.164 8.201 ( -1.977 0.751 -5.117) 5.537 8.466 ( 1.888 -3.789 0.474) 4.260 8.479 ( 2.185 -4.163 -0.531) 4.732 8.712 ( 0.166 -1.363 0.052) 1.374 8.713 ( 0.141 -1.439 -0.045) 1.447 ======================= Grid point 42 (18/24) ======================= q-point: ( 0.40 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 53 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.625 ( 7.099 4.303 1.352) 8.411 0.687 ( 3.552 6.864 1.424) 7.858 0.721 ( 1.075 6.302 -3.537) 7.306 0.754 ( -2.091 2.896 -6.746) 7.633 0.781 ( 6.194 6.061 2.342) 8.977 0.840 ( 3.232 7.867 -1.859) 8.706 0.908 ( 2.089 4.407 0.389) 4.892 0.930 ( 1.779 3.118 -1.187) 3.781 0.959 ( -1.686 1.208 1.425) 2.516 0.999 ( -3.758 0.804 -1.210) 4.029 1.028 ( -1.854 0.530 0.405) 1.970 1.036 ( -1.484 -0.345 -0.575) 1.629 1.069 ( 0.680 2.013 0.227) 2.137 1.087 ( 0.275 1.750 -0.975) 2.022 1.138 ( 2.765 2.976 1.102) 4.209 1.191 ( 1.066 0.565 -1.492) 1.919 1.214 ( -0.658 1.247 1.783) 2.273 1.236 ( -1.489 2.217 1.500) 3.063 1.261 ( -0.489 0.313 -1.337) 1.458 1.266 ( -0.633 0.357 -0.317) 0.793 1.305 ( 1.222 2.371 4.542) 5.267 1.378 ( -1.544 0.447 2.468) 2.945 1.414 ( 1.480 2.201 -0.790) 2.768 1.418 ( -0.700 -1.926 1.025) 2.291 1.439 ( -2.076 -2.012 -1.298) 3.169 1.467 ( -1.982 -0.218 -1.129) 2.292 1.478 ( 0.806 2.336 1.468) 2.875 1.506 ( -1.173 2.829 -0.695) 3.140 1.559 ( 0.649 -2.120 -0.713) 2.329 1.568 ( -0.352 -1.307 -1.503) 2.023 1.651 ( 1.297 2.040 3.879) 4.571 1.667 ( -0.037 2.481 1.876) 3.110 1.698 ( 1.038 1.709 -0.827) 2.164 1.728 ( 0.214 0.826 -2.916) 3.038 1.755 ( 3.720 1.867 1.410) 4.394 1.812 ( 1.417 0.640 -1.875) 2.436 1.873 ( 1.128 2.814 -0.355) 3.053 1.888 ( -0.081 3.437 -0.912) 3.557 1.921 ( -2.858 0.439 0.307) 2.908 1.934 ( -3.740 0.773 -0.346) 3.835 2.063 ( -0.495 4.033 -0.465) 4.090 2.068 ( -1.088 2.992 -0.943) 3.321 2.224 ( 3.928 3.320 2.181) 5.587 2.264 ( 2.006 1.790 -0.925) 2.844 2.345 ( 0.446 2.880 0.540) 2.964 2.361 ( 0.185 1.974 -0.741) 2.117 2.429 ( -1.462 -1.689 0.712) 2.345 2.440 ( -1.878 -1.763 0.154) 2.581 2.485 ( 0.126 0.220 1.310) 1.334 2.496 ( -0.891 -0.808 1.435) 1.872 2.518 ( -0.180 1.169 -1.672) 2.049 2.529 ( -1.350 1.115 -1.611) 2.380 2.792 ( 0.870 3.188 0.183) 3.310 2.797 ( 0.519 3.576 -0.171) 3.617 2.894 ( -0.017 4.915 1.427) 5.118 2.906 ( -0.769 4.173 1.662) 4.557 2.983 ( -0.014 -0.561 -0.816) 0.990 3.013 ( -2.369 -0.580 -2.167) 3.263 3.049 ( 0.716 -0.451 0.649) 1.066 3.058 ( -0.150 -0.388 0.095) 0.426 3.192 ( 1.037 -6.810 -0.145) 6.890 3.207 ( -0.227 -6.789 -1.268) 6.910 3.336 ( 1.846 -0.448 1.224) 2.260 3.351 ( 1.115 -0.903 -0.158) 1.444 3.369 ( -0.763 0.304 0.348) 0.892 3.384 ( -1.229 0.869 -0.867) 1.737 3.582 ( 0.486 -2.805 0.929) 2.995 3.613 ( -2.125 -2.520 -1.587) 3.658 3.731 ( 0.289 1.690 0.156) 1.721 3.736 ( -0.171 2.390 -0.234) 2.407 3.750 ( -0.049 1.395 0.563) 1.505 3.762 ( -0.896 0.354 -0.316) 1.014 3.916 ( 1.447 -2.909 -0.273) 3.260 3.920 ( 1.149 -2.931 -0.656) 3.216 3.993 ( -0.101 -1.287 0.916) 1.583 4.002 ( -0.790 -0.953 0.047) 1.239 6.657 ( 1.245 -0.461 0.237) 1.348 6.671 ( 0.132 -0.891 -0.937) 1.300 6.778 ( 0.213 -3.382 0.184) 3.394 6.784 ( -0.365 -2.907 -0.318) 2.947 7.064 ( 3.062 2.130 1.663) 4.084 7.085 ( 2.420 1.371 0.496) 2.826 7.241 ( -0.641 -3.521 2.853) 4.576 7.299 ( -4.319 -4.237 -1.632) 6.267 7.510 ( 0.366 0.830 0.010) 0.907 7.522 ( -0.645 1.279 -0.608) 1.556 7.543 ( 0.349 -0.147 0.259) 0.459 7.550 ( -0.268 -0.189 -0.270) 0.424 7.679 ( -0.664 9.351 0.830) 9.412 7.717 ( -4.212 8.802 -2.374) 10.043 8.114 ( 2.629 -0.102 2.211) 3.436 8.163 ( -1.859 -0.324 -1.852) 2.644 8.489 ( 0.668 -4.751 0.285) 4.806 8.499 ( -0.124 -4.881 -0.526) 4.911 8.716 ( 0.168 -1.395 0.004) 1.405 8.717 ( 0.132 -1.430 -0.021) 1.436 ======================= Grid point 45 (19/24) ======================= q-point: ( 0.00 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 33 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.614 ( 0.000 7.726 5.146) 9.283 0.615 ( 0.000 6.223 3.355) 7.069 0.644 ( 0.000 9.819 4.563) 10.828 0.731 ( 0.000 4.834 -2.658) 5.516 0.748 ( 0.000 4.221 -3.526) 5.501 0.800 ( 0.000 3.001 -6.684) 7.327 0.997 ( 0.000 -1.254 1.188) 1.727 1.022 ( 0.000 1.304 -0.636) 1.451 1.051 ( 0.000 3.258 2.506) 4.110 1.057 ( 0.000 4.433 0.402) 4.452 1.096 ( 0.000 1.092 -1.513) 1.866 1.114 ( 0.000 0.448 -0.485) 0.660 1.128 ( 0.000 4.452 -0.930) 4.548 1.131 ( 0.000 -3.128 -1.343) 3.405 1.142 ( 0.000 3.231 -1.235) 3.459 1.190 ( 0.000 1.457 -1.087) 1.818 1.256 ( 0.000 2.249 2.328) 3.237 1.279 ( 0.000 -1.464 0.931) 1.735 1.286 ( 0.000 -1.497 2.310) 2.753 1.328 ( 0.000 5.399 0.955) 5.482 1.335 ( 0.000 -0.644 2.517) 2.598 1.360 ( 0.000 -1.904 -0.109) 1.907 1.408 ( 0.000 4.082 -2.673) 4.879 1.429 ( 0.000 2.369 -2.958) 3.790 1.452 ( 0.000 0.956 2.421) 2.603 1.461 ( 0.000 -1.803 1.862) 2.592 1.481 ( 0.000 -1.774 -1.742) 2.486 1.539 ( 0.000 -1.702 -0.377) 1.743 1.575 ( 0.000 0.437 0.293) 0.526 1.604 ( 0.000 1.297 0.275) 1.326 1.648 ( 0.000 1.271 -1.458) 1.934 1.658 ( 0.000 1.843 -1.554) 2.411 1.705 ( 0.000 3.186 2.530) 4.069 1.713 ( 0.000 5.841 0.606) 5.872 1.756 ( 0.000 -0.328 -0.891) 0.949 1.787 ( 0.000 1.217 -2.705) 2.966 1.903 ( 0.000 3.018 1.522) 3.380 1.938 ( 0.000 3.984 -1.964) 4.442 2.038 ( 0.000 3.425 2.074) 4.004 2.067 ( 0.000 3.386 -0.559) 3.432 2.132 ( 0.000 0.955 1.568) 1.836 2.147 ( 0.000 -0.461 1.021) 1.120 2.216 ( 0.000 3.520 -0.214) 3.527 2.218 ( 0.000 3.290 -0.949) 3.424 2.387 ( 0.000 0.152 -1.966) 1.972 2.389 ( 0.000 1.691 -1.747) 2.431 2.469 ( 0.000 0.572 0.526) 0.777 2.474 ( 0.000 -0.122 0.026) 0.125 2.512 ( 0.000 4.446 2.074) 4.905 2.549 ( 0.000 3.256 -1.323) 3.514 2.582 ( 0.000 4.270 1.886) 4.668 2.621 ( 0.000 3.397 -1.742) 3.818 2.805 ( 0.000 3.357 1.269) 3.589 2.829 ( 0.000 1.990 -1.211) 2.330 2.917 ( 0.000 -2.266 -2.188) 3.150 2.919 ( 0.000 -2.365 -1.915) 3.043 3.011 ( 0.000 -2.183 1.081) 2.436 3.030 ( 0.000 -3.011 -0.399) 3.038 3.056 ( 0.000 0.421 1.362) 1.426 3.078 ( 0.000 -0.206 -0.188) 0.279 3.192 ( 0.000 1.988 0.330) 2.016 3.200 ( 0.000 1.824 -0.321) 1.852 3.238 ( 0.000 0.342 2.959) 2.979 3.318 ( 0.000 -1.512 -2.761) 3.148 3.368 ( 0.000 -2.971 0.091) 2.973 3.407 ( 0.000 -4.245 -2.124) 4.746 3.519 ( 0.000 -2.171 3.046) 3.740 3.575 ( 0.000 -4.549 -1.430) 4.768 3.729 ( 0.000 0.022 2.407) 2.407 3.763 ( 0.000 -1.458 1.272) 1.935 3.777 ( 0.000 -0.034 -1.313) 1.313 3.801 ( 0.000 -0.064 -1.945) 1.946 3.844 ( 0.000 2.786 0.190) 2.793 3.856 ( 0.000 2.921 -0.634) 2.989 3.975 ( 0.000 -0.997 2.291) 2.499 4.009 ( 0.000 -0.848 -0.823) 1.181 6.623 ( 0.000 1.011 0.571) 1.161 6.641 ( 0.000 1.279 -0.894) 1.560 6.723 ( 0.000 -1.800 1.044) 2.081 6.762 ( 0.000 -1.894 -2.600) 3.217 7.018 ( 0.000 2.028 1.228) 2.371 7.051 ( 0.000 2.627 -1.764) 3.164 7.231 ( 0.000 -5.242 4.656) 7.011 7.297 ( 0.000 -5.324 -1.679) 5.582 7.450 ( 0.000 -1.776 2.742) 3.267 7.515 ( 0.000 -3.133 -2.479) 3.995 7.566 ( 0.000 1.359 0.203) 1.374 7.572 ( 0.000 0.641 -0.314) 0.713 7.890 ( 0.000 3.593 0.632) 3.649 7.913 ( 0.000 4.708 -1.359) 4.901 8.165 ( 0.000 0.840 3.526) 3.625 8.242 ( 0.000 -1.682 -2.883) 3.337 8.397 ( 0.000 0.021 0.545) 0.545 8.420 ( 0.000 2.855 -1.602) 3.274 8.667 ( 0.000 -3.246 0.257) 3.256 8.671 ( 0.000 -2.844 -0.107) 2.846 ======================= Grid point 46 (20/24) ======================= q-point: ( 0.20 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 53 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.652 ( 2.844 8.565 4.375) 10.030 0.685 ( 5.729 5.773 3.698) 8.935 0.689 ( 4.925 6.505 3.044) 8.708 0.777 ( 3.562 5.098 -2.620) 6.748 0.785 ( 3.135 3.615 -4.167) 6.345 0.830 ( 2.431 2.221 -5.832) 6.698 0.993 ( -0.588 -0.685 0.971) 1.326 1.020 ( -0.460 -0.680 -1.960) 2.125 1.046 ( -0.418 3.109 1.839) 3.636 1.058 ( -0.489 3.833 -0.780) 3.942 1.080 ( -1.613 0.605 -1.966) 2.614 1.092 ( -0.985 -0.260 -1.589) 1.887 1.112 ( -1.106 -0.517 -0.067) 1.223 1.129 ( 0.061 3.989 -1.001) 4.113 1.168 ( -0.264 2.864 1.515) 3.251 1.183 ( 1.739 2.923 0.201) 3.408 1.241 ( -0.308 2.008 2.562) 3.270 1.291 ( -0.048 -0.366 -0.140) 0.395 1.304 ( 0.902 0.167 0.776) 1.202 1.320 ( -0.212 2.362 0.100) 2.374 1.345 ( 0.099 0.952 2.270) 2.463 1.363 ( -0.903 0.031 -0.450) 1.009 1.381 ( -0.191 0.882 -0.156) 0.916 1.427 ( -0.335 1.194 -0.936) 1.553 1.454 ( 0.097 -1.159 2.384) 2.652 1.470 ( -0.586 -1.997 0.743) 2.210 1.494 ( -0.398 0.394 0.289) 0.631 1.530 ( 0.146 0.783 -2.397) 2.526 1.546 ( 0.303 0.507 0.058) 0.593 1.581 ( -1.463 0.124 -2.204) 2.648 1.668 ( 0.021 1.329 2.210) 2.578 1.684 ( 2.641 2.835 0.331) 3.889 1.720 ( 0.256 5.298 1.595) 5.538 1.729 ( 1.724 2.160 1.195) 3.011 1.759 ( 1.927 0.646 -1.784) 2.704 1.818 ( 1.508 2.392 -4.087) 4.970 1.909 ( 0.265 1.851 2.193) 2.882 1.950 ( 0.462 3.515 -1.385) 3.806 2.022 ( -1.414 3.874 1.815) 4.506 2.046 ( -2.106 4.101 -0.452) 4.632 2.130 ( -0.524 0.834 0.551) 1.129 2.136 ( -0.538 1.187 0.351) 1.349 2.222 ( 1.314 1.918 0.300) 2.345 2.241 ( 2.055 1.705 -0.988) 2.847 2.371 ( -0.789 0.475 -1.339) 1.626 2.383 ( -0.649 0.891 -2.191) 2.452 2.449 ( -1.458 -0.317 0.372) 1.538 2.459 ( -1.129 -0.302 -0.216) 1.188 2.524 ( 0.335 3.970 2.775) 4.856 2.559 ( -1.200 4.452 1.038) 4.726 2.572 ( 1.010 5.298 -0.100) 5.395 2.603 ( -1.494 3.546 -1.758) 4.231 2.827 ( 1.651 3.121 1.101) 3.698 2.849 ( 1.357 2.475 -0.987) 2.990 2.917 ( 0.401 -2.393 -0.944) 2.604 2.933 ( 1.137 -1.176 -2.502) 2.990 3.019 ( -1.629 -2.008 0.375) 2.613 3.022 ( -0.254 -1.055 1.678) 1.998 3.045 ( 0.548 -2.926 -0.109) 2.979 3.066 ( -1.045 -0.571 -0.446) 1.272 3.130 ( -2.974 2.908 1.557) 4.441 3.163 ( -2.974 2.162 -1.055) 3.825 3.289 ( 2.380 0.385 0.781) 2.534 3.318 ( 0.255 -1.306 -1.854) 2.283 3.363 ( -0.170 -2.245 1.687) 2.813 3.403 ( -0.623 -2.465 -1.648) 3.030 3.495 ( -0.783 -4.555 3.561) 5.834 3.567 ( -1.344 -5.402 -2.423) 6.071 3.733 ( 0.734 1.005 1.635) 2.055 3.769 ( -0.163 1.038 -1.277) 1.654 3.791 ( 0.898 -1.384 -0.118) 1.654 3.806 ( -0.135 0.455 -1.546) 1.617 3.843 ( 0.520 1.540 1.028) 1.923 3.861 ( 0.524 1.485 -0.469) 1.643 3.970 ( -0.403 -1.305 1.939) 2.371 3.999 ( -0.888 -1.153 -0.697) 1.614 6.630 ( 0.635 0.437 0.578) 0.964 6.647 ( 0.493 0.468 -0.836) 1.078 6.721 ( -0.233 -2.247 0.909) 2.435 6.753 ( -0.879 -2.169 -1.993) 3.074 7.044 ( 2.279 1.488 0.857) 2.854 7.065 ( 1.195 2.080 -1.020) 2.606 7.201 ( -2.035 -4.395 4.270) 6.456 7.270 ( -2.331 -5.274 -1.883) 6.066 7.478 ( 1.909 -1.332 1.687) 2.875 7.524 ( 0.515 -2.024 -1.976) 2.875 7.559 ( -0.660 1.007 0.208) 1.222 7.564 ( -0.747 0.477 -0.186) 0.905 7.879 ( -0.250 6.737 1.566) 6.921 7.931 ( 0.691 4.637 -3.042) 5.589 8.125 ( -2.554 2.105 4.152) 5.310 8.207 ( -2.787 -0.807 -2.922) 4.118 8.400 ( 0.127 -1.278 0.714) 1.469 8.422 ( -0.170 -0.195 -1.163) 1.192 8.667 ( 0.053 -3.392 0.136) 3.395 8.670 ( -0.065 -3.106 -0.070) 3.108 ======================= Grid point 47 (21/24) ======================= q-point: ( 0.40 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 53 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.733 ( 4.369 6.074 0.758) 7.520 0.791 ( 1.911 5.118 -1.585) 5.689 0.813 ( 2.861 4.470 -1.139) 5.428 0.828 ( 1.956 4.245 -1.467) 4.898 0.856 ( 4.504 3.330 0.573) 5.630 0.909 ( 4.670 0.028 -2.208) 5.165 0.980 ( -3.297 -0.227 0.569) 3.354 0.981 ( -0.439 0.470 -0.202) 0.674 1.022 ( -2.469 0.679 -0.158) 2.566 1.036 ( -2.691 2.143 -1.433) 3.727 1.051 ( 0.967 5.572 1.237) 5.789 1.075 ( -0.876 3.093 -0.413) 3.241 1.112 ( 1.125 1.013 -0.736) 1.684 1.124 ( 0.156 2.054 -1.790) 2.729 1.166 ( 0.636 0.408 0.746) 1.062 1.195 ( -1.201 -1.018 -1.065) 1.900 1.260 ( 1.053 2.395 0.871) 2.758 1.278 ( -0.449 1.240 -0.383) 1.373 1.291 ( -1.666 2.916 0.660) 3.423 1.305 ( -2.095 3.698 -0.509) 4.281 1.340 ( 0.999 -0.920 1.201) 1.813 1.356 ( 0.537 -2.802 -0.610) 2.917 1.399 ( 0.667 2.102 2.076) 3.028 1.423 ( -0.771 1.082 -0.258) 1.353 1.438 ( 0.019 -1.454 0.636) 1.587 1.449 ( -0.768 -2.515 1.012) 2.818 1.482 ( -0.229 1.181 1.589) 1.993 1.502 ( -1.573 0.873 -0.451) 1.854 1.585 ( 1.969 1.607 -0.498) 2.590 1.615 ( 1.605 3.196 -2.771) 4.524 1.652 ( 0.745 -0.782 3.899) 4.046 1.701 ( -1.804 0.344 0.715) 1.971 1.755 ( 2.346 3.636 0.588) 4.367 1.768 ( 1.763 4.399 1.492) 4.968 1.809 ( 1.715 2.424 -0.374) 2.993 1.838 ( 0.421 2.686 -2.813) 3.912 1.901 ( -0.476 1.098 1.507) 1.925 1.931 ( -2.245 1.565 0.214) 2.745 1.993 ( -0.775 5.043 0.304) 5.111 2.000 ( -1.478 5.157 -0.996) 5.456 2.124 ( 0.185 2.057 -0.725) 2.189 2.130 ( -0.200 2.230 -0.897) 2.412 2.265 ( 1.879 -0.465 0.864) 2.119 2.281 ( 0.675 -0.923 -0.445) 1.227 2.375 ( 1.383 0.543 -0.810) 1.692 2.382 ( 0.906 0.114 -1.262) 1.558 2.421 ( -1.183 -0.692 -0.162) 1.380 2.428 ( -1.820 -0.131 -0.529) 1.899 2.512 ( -0.324 4.749 2.415) 5.338 2.530 ( -1.329 5.486 1.921) 5.962 2.574 ( -0.595 3.927 -1.390) 4.208 2.583 ( -0.617 5.247 -1.301) 5.441 2.859 ( 1.326 3.155 0.581) 3.471 2.871 ( 0.571 3.119 -0.280) 3.183 2.931 ( 0.632 -1.301 0.904) 1.705 2.958 ( -2.483 0.402 1.843) 3.118 2.967 ( 1.822 -0.921 -2.170) 2.980 3.006 ( -1.869 -1.241 -1.762) 2.852 3.055 ( 0.769 -4.530 0.452) 4.617 3.064 ( 0.281 -4.353 0.118) 4.364 3.091 ( -0.687 4.385 0.773) 4.505 3.107 ( -1.911 3.271 -0.378) 3.807 3.323 ( 0.996 -1.070 0.494) 1.543 3.330 ( 0.379 -1.071 0.111) 1.142 3.363 ( 0.357 -0.978 0.524) 1.165 3.382 ( -1.194 -1.007 -0.951) 1.829 3.501 ( 1.082 -5.666 1.150) 5.882 3.532 ( -1.670 -5.850 -1.289) 6.219 3.752 ( 0.985 1.162 0.809) 1.725 3.768 ( -0.177 0.956 -0.362) 1.038 3.786 ( -0.535 1.662 -0.779) 1.912 3.793 ( -0.772 1.759 -1.173) 2.251 3.866 ( 1.012 -1.311 0.656) 1.781 3.873 ( 0.324 -0.910 0.120) 0.973 3.965 ( 0.107 -1.669 0.886) 1.892 3.977 ( -0.975 -1.599 -0.195) 1.883 6.648 ( 0.842 -0.412 0.271) 0.976 6.656 ( 0.158 -0.494 -0.390) 0.649 6.716 ( -0.000 -2.855 0.349) 2.876 6.729 ( -1.138 -2.706 -0.736) 3.027 7.092 ( 1.577 0.793 0.764) 1.924 7.100 ( 1.555 0.348 0.313) 1.624 7.171 ( -0.096 -3.277 2.063) 3.874 7.211 ( -2.921 -4.208 -1.173) 5.255 7.513 ( 1.331 -0.236 0.204) 1.367 7.529 ( -0.012 -0.250 -1.032) 1.062 7.542 ( -0.684 0.044 -0.022) 0.686 7.546 ( -0.798 -0.152 -0.275) 0.858 7.889 ( 1.104 10.037 1.082) 10.155 7.921 ( -1.545 9.700 -1.671) 9.964 8.104 ( 0.833 -1.016 2.066) 2.448 8.145 ( -2.704 -1.558 -1.451) 3.441 8.402 ( 0.158 -2.542 0.306) 2.566 8.411 ( -0.667 -2.515 -0.460) 2.642 8.669 ( 0.103 -3.612 0.030) 3.613 8.670 ( 0.048 -3.529 -0.009) 3.529 ======================= Grid point 50 (22/24) ======================= q-point: ( 0.00 0.43 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 33 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.711 ( 0.000 3.832 -0.839) 3.922 0.725 ( 0.000 2.821 -3.100) 4.192 0.809 ( 0.000 4.174 0.902) 4.270 0.818 ( 0.000 2.844 2.195) 3.593 0.846 ( 0.000 6.257 0.203) 6.260 0.874 ( 0.000 5.361 -0.829) 5.424 0.971 ( 0.000 -0.222 0.361) 0.424 0.986 ( 0.000 -2.665 -0.744) 2.767 1.007 ( 0.000 -6.442 0.447) 6.458 1.020 ( 0.000 -4.826 -0.770) 4.887 1.093 ( 0.000 0.360 1.401) 1.447 1.114 ( 0.000 -1.178 -0.015) 1.178 1.135 ( 0.000 2.131 -1.252) 2.471 1.158 ( 0.000 -1.417 -2.340) 2.736 1.197 ( 0.000 1.942 -2.487) 3.155 1.205 ( 0.000 0.593 -2.250) 2.327 1.264 ( 0.000 1.102 0.781) 1.351 1.273 ( 0.000 0.251 -0.112) 0.275 1.295 ( 0.000 0.599 1.742) 1.842 1.323 ( 0.000 -1.510 -0.322) 1.544 1.365 ( 0.000 -1.381 1.587) 2.103 1.370 ( 0.000 -0.602 1.770) 1.870 1.418 ( 0.000 0.340 2.156) 2.182 1.444 ( 0.000 -1.555 -2.248) 2.734 1.503 ( 0.000 5.906 3.725) 6.982 1.512 ( 0.000 4.718 -0.222) 4.723 1.538 ( 0.000 3.378 -1.090) 3.549 1.546 ( 0.000 3.294 -2.664) 4.236 1.583 ( 0.000 -0.576 0.461) 0.738 1.600 ( 0.000 -1.715 1.871) 2.538 1.659 ( 0.000 -0.205 -0.389) 0.440 1.660 ( 0.000 -0.106 -0.486) 0.498 1.701 ( 0.000 -2.821 1.287) 3.101 1.702 ( 0.000 -3.430 1.388) 3.700 1.827 ( 0.000 2.302 0.189) 2.310 1.829 ( 0.000 2.891 -0.311) 2.908 1.977 ( 0.000 4.837 1.923) 5.206 2.011 ( 0.000 2.190 -1.395) 2.597 2.116 ( 0.000 3.129 -1.139) 3.330 2.126 ( 0.000 2.058 -2.143) 2.971 2.135 ( 0.000 -0.110 1.634) 1.637 2.141 ( 0.000 0.174 1.012) 1.027 2.224 ( 0.000 -1.868 -0.077) 1.869 2.231 ( 0.000 -1.731 -0.888) 1.946 2.384 ( 0.000 0.262 -0.994) 1.028 2.394 ( 0.000 -0.471 -1.658) 1.723 2.453 ( 0.000 -0.602 0.485) 0.773 2.458 ( 0.000 -0.917 0.077) 0.920 2.644 ( 0.000 6.114 1.905) 6.404 2.658 ( 0.000 5.090 0.883) 5.166 2.695 ( 0.000 5.980 0.012) 5.980 2.718 ( 0.000 5.259 -2.595) 5.864 2.812 ( 0.000 -2.932 0.017) 2.932 2.840 ( 0.000 -3.203 -2.660) 4.164 2.864 ( 0.000 0.234 0.407) 0.469 2.880 ( 0.000 -1.488 -0.452) 1.555 2.970 ( 0.000 -1.667 -0.423) 1.720 2.971 ( 0.000 -2.325 -0.096) 2.327 3.071 ( 0.000 0.630 1.189) 1.346 3.077 ( 0.000 -0.147 0.654) 0.670 3.220 ( 0.000 1.606 0.678) 1.743 3.241 ( 0.000 2.035 -1.002) 2.268 3.271 ( 0.000 2.649 2.169) 3.423 3.274 ( 0.000 -2.107 1.731) 2.727 3.317 ( 0.000 -1.816 -1.722) 2.503 3.318 ( 0.000 -4.050 -1.904) 4.475 3.465 ( 0.000 -2.747 0.966) 2.912 3.476 ( 0.000 -4.313 0.059) 4.314 3.740 ( 0.000 0.673 2.243) 2.342 3.748 ( 0.000 -0.229 2.083) 2.096 3.796 ( 0.000 1.223 -1.937) 2.290 3.806 ( 0.000 0.218 -2.389) 2.399 3.887 ( 0.000 1.626 0.423) 1.680 3.899 ( 0.000 0.855 -0.561) 1.022 3.966 ( 0.000 0.394 1.386) 1.441 3.988 ( 0.000 -1.333 -0.496) 1.422 6.647 ( 0.000 1.384 0.529) 1.482 6.664 ( 0.000 0.466 -0.826) 0.949 6.698 ( 0.000 -0.213 0.487) 0.532 6.723 ( 0.000 -1.990 -1.692) 2.612 7.049 ( 0.000 1.462 1.935) 2.426 7.074 ( 0.000 -0.536 0.454) 0.702 7.175 ( 0.000 -0.073 1.199) 1.201 7.212 ( 0.000 -3.361 -1.472) 3.669 7.432 ( 0.000 0.105 0.803) 0.810 7.457 ( 0.000 -2.410 -1.260) 2.720 7.581 ( 0.000 0.283 -0.076) 0.293 7.582 ( 0.000 0.324 -0.096) 0.338 7.946 ( 0.000 2.072 0.689) 2.183 7.965 ( 0.000 0.504 -0.885) 1.019 8.165 ( 0.000 0.115 1.111) 1.117 8.191 ( 0.000 -2.516 -1.093) 2.743 8.464 ( 0.000 5.760 1.571) 5.970 8.504 ( 0.000 4.695 -1.808) 5.031 8.585 ( 0.000 -5.032 0.902) 5.113 8.602 ( 0.000 -4.187 -0.618) 4.232 ======================= Grid point 51 (23/24) ======================= q-point: ( 0.20 0.43 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 53 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.762 ( 3.900 4.507 1.143) 6.069 0.797 ( 5.535 2.654 -2.853) 6.769 0.824 ( 1.634 5.014 0.834) 5.339 0.842 ( 2.272 1.525 -0.041) 2.737 0.878 ( 2.086 5.323 0.096) 5.718 0.887 ( 1.635 3.471 -0.851) 3.930 0.970 ( -0.049 -0.686 0.135) 0.701 0.981 ( -0.065 -1.373 -0.903) 1.644 1.022 ( 1.461 -4.354 0.342) 4.605 1.028 ( 1.101 -3.719 -0.437) 3.903 1.091 ( -0.259 0.589 1.036) 1.219 1.110 ( -0.582 -0.538 -0.412) 0.893 1.133 ( -0.930 2.313 -2.284) 3.381 1.159 ( -0.529 -1.898 -2.071) 2.859 1.179 ( -1.501 0.639 -0.454) 1.694 1.193 ( -1.056 0.090 -1.614) 1.931 1.273 ( -0.025 1.734 1.157) 2.085 1.284 ( 1.161 0.892 0.278) 1.490 1.303 ( 0.847 -0.330 0.069) 0.912 1.314 ( -0.097 -1.409 -0.672) 1.565 1.341 ( -1.155 -1.303 0.751) 1.896 1.357 ( -1.326 -2.330 0.028) 2.681 1.428 ( 0.512 0.860 1.406) 1.726 1.453 ( 1.009 -0.068 -1.272) 1.625 1.476 ( -2.122 2.225 1.270) 3.327 1.490 ( -2.119 2.696 1.342) 3.683 1.525 ( -1.320 1.707 -1.018) 2.386 1.529 ( -1.806 0.828 -0.350) 2.017 1.561 ( 0.864 2.806 3.894) 4.877 1.607 ( 0.838 0.516 -0.007) 0.985 1.672 ( -0.572 -0.690 1.756) 1.972 1.684 ( 0.614 -1.915 0.714) 2.134 1.727 ( 3.203 0.092 0.473) 3.239 1.745 ( 3.852 -1.517 -1.466) 4.392 1.851 ( 1.534 3.482 0.137) 3.807 1.862 ( 2.470 1.845 -0.433) 3.113 1.965 ( -0.697 4.828 1.871) 5.225 2.016 ( 0.178 1.978 -2.133) 2.915 2.093 ( -2.259 0.920 0.649) 2.524 2.110 ( -1.014 2.209 -1.191) 2.707 2.131 ( -0.370 0.327 1.591) 1.666 2.162 ( 0.831 2.432 -0.939) 2.736 2.209 ( -0.558 -2.539 0.306) 2.617 2.229 ( -0.040 -2.456 -1.069) 2.679 2.377 ( -0.440 0.416 -0.707) 0.931 2.385 ( -0.622 -0.331 -1.397) 1.564 2.437 ( -1.165 -0.282 0.381) 1.257 2.445 ( -0.946 -0.910 -0.155) 1.321 2.638 ( -0.497 5.807 1.744) 6.084 2.672 ( 0.598 3.781 0.430) 3.852 2.677 ( -1.665 5.461 -0.788) 5.764 2.702 ( -0.914 4.643 -1.781) 5.056 2.838 ( 2.045 -2.729 0.035) 3.411 2.854 ( 1.279 -3.259 -1.513) 3.814 2.872 ( 0.652 1.035 1.089) 1.638 2.898 ( 0.842 -1.936 -1.124) 2.392 2.965 ( -0.048 -1.510 -0.167) 1.520 2.970 ( 0.113 -2.090 -0.450) 2.141 3.056 ( -1.433 1.261 1.307) 2.313 3.066 ( -1.040 0.072 0.488) 1.151 3.209 ( -1.325 4.699 1.415) 5.083 3.232 ( 0.517 4.079 -0.124) 4.113 3.285 ( 1.759 -1.406 0.112) 2.255 3.294 ( -0.459 -1.926 -0.825) 2.145 3.314 ( 1.695 -0.668 0.551) 1.904 3.343 ( 0.905 -2.951 -1.635) 3.492 3.409 ( -3.488 -3.066 2.444) 5.248 3.450 ( -2.621 -5.859 -1.020) 6.499 3.755 ( 1.132 1.002 1.611) 2.209 3.769 ( 1.424 -0.635 1.080) 1.897 3.795 ( -0.118 1.207 -1.493) 1.924 3.807 ( 0.051 -0.068 -1.893) 1.895 3.882 ( -0.318 1.854 0.680) 2.000 3.897 ( -0.142 0.940 -0.580) 1.114 3.952 ( -1.185 -0.082 1.357) 1.804 3.972 ( -1.322 -1.531 -0.440) 2.070 6.645 ( -0.123 1.085 0.544) 1.220 6.661 ( -0.261 0.326 -0.721) 0.833 6.688 ( -0.845 -0.581 0.516) 1.148 6.711 ( -1.115 -2.038 -1.407) 2.716 7.067 ( 1.474 1.057 1.486) 2.345 7.085 ( 0.921 -0.313 0.395) 1.050 7.151 ( -1.555 -0.146 1.225) 1.985 7.186 ( -1.909 -3.220 -1.292) 3.959 7.464 ( 2.350 0.117 0.449) 2.395 7.484 ( 1.848 -1.765 -1.116) 2.789 7.571 ( -0.945 0.276 -0.009) 0.985 7.571 ( -0.932 0.255 -0.062) 0.968 7.962 ( 1.535 2.178 0.616) 2.735 7.980 ( 1.486 0.464 -0.890) 1.793 8.142 ( -2.095 0.393 1.167) 2.430 8.168 ( -2.165 -2.304 -1.035) 3.327 8.455 ( -0.771 5.908 1.169) 6.071 8.484 ( -1.700 5.298 -1.256) 5.704 8.578 ( -0.414 -5.687 0.635) 5.737 8.591 ( -0.819 -4.927 -0.463) 5.016 ======================= Grid point 52 (24/24) ======================= q-point: ( 0.40 0.43 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 5.20e-05 5.20e-05 5.20e-05 5.20e-05 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 5.83e-04 Number of triplets: 53 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.831 ( 2.372 3.782 1.052) 4.587 0.859 ( 0.651 2.707 -0.535) 2.835 0.886 ( 1.543 -1.502 0.074) 2.154 0.889 ( 0.843 0.865 -0.357) 1.259 0.933 ( 3.494 5.243 -0.906) 6.365 0.943 ( 3.162 3.808 -1.648) 5.217 0.974 ( 0.678 0.819 0.005) 1.063 0.992 ( 0.312 1.341 -1.198) 1.825 1.056 ( 0.210 1.207 -0.179) 1.238 1.066 ( 1.360 -0.900 -0.506) 1.708 1.073 ( -0.836 1.240 -0.096) 1.498 1.092 ( -1.715 1.165 -1.080) 2.338 1.105 ( 0.910 -1.737 0.607) 2.052 1.137 ( -0.877 -2.364 -1.940) 3.182 1.146 ( -1.575 -2.407 -0.399) 2.904 1.153 ( -2.020 -2.473 -1.009) 3.349 1.283 ( 0.715 0.894 -0.071) 1.147 1.293 ( 0.206 0.754 -0.918) 1.205 1.315 ( 0.360 -1.235 -0.681) 1.455 1.321 ( 0.168 -1.071 -0.931) 1.429 1.337 ( 1.049 1.033 0.677) 1.620 1.349 ( 0.450 -0.722 -0.412) 0.945 1.404 ( -2.927 0.482 1.711) 3.424 1.408 ( -3.858 0.826 1.263) 4.143 1.463 ( 0.486 3.310 0.742) 3.427 1.471 ( 0.201 1.706 -0.481) 1.784 1.513 ( 0.429 -0.456 0.466) 0.781 1.518 ( 0.341 -0.737 0.048) 0.814 1.602 ( 2.426 2.032 3.145) 4.462 1.636 ( 1.712 -0.458 0.515) 1.846 1.648 ( -1.651 0.252 2.374) 2.902 1.672 ( -1.411 -2.264 0.934) 2.826 1.807 ( 3.616 1.776 0.697) 4.089 1.828 ( 3.306 0.102 -0.569) 3.356 1.867 ( -0.399 2.145 0.421) 2.222 1.889 ( -0.724 0.235 -0.563) 0.947 1.958 ( 0.594 6.702 1.128) 6.822 1.995 ( -2.221 5.425 -0.717) 5.905 2.062 ( -0.135 -0.516 -0.581) 0.789 2.078 ( -1.486 0.095 -1.563) 2.159 2.134 ( 0.831 1.804 0.828) 2.152 2.156 ( -0.941 2.213 -0.489) 2.454 2.216 ( 0.666 -4.139 -0.123) 4.194 2.227 ( -0.320 -3.868 -0.762) 3.955 2.380 ( 0.906 0.143 -0.468) 1.030 2.386 ( 1.136 -0.244 -0.849) 1.439 2.411 ( -1.456 0.422 0.220) 1.532 2.421 ( -1.510 -0.348 -0.555) 1.646 2.630 ( 0.476 5.311 0.470) 5.353 2.649 ( -0.123 4.022 -0.958) 4.136 2.661 ( -1.445 4.344 0.567) 4.613 2.680 ( -1.728 3.182 -0.714) 3.691 2.875 ( 1.080 -1.860 0.770) 2.284 2.883 ( 0.332 -0.470 0.830) 1.010 2.886 ( 1.397 -2.826 0.051) 3.153 2.904 ( -0.134 -3.008 -0.725) 3.097 2.976 ( 1.159 0.547 0.060) 1.283 2.990 ( 1.421 0.025 -1.068) 1.778 3.017 ( -1.875 1.316 1.039) 2.515 3.033 ( -1.911 -0.069 -0.187) 1.921 3.203 ( 0.758 6.720 1.466) 6.920 3.227 ( -1.085 7.199 -0.374) 7.290 3.304 ( 0.886 -0.597 0.566) 1.209 3.314 ( 0.208 -0.365 0.022) 0.420 3.333 ( 0.667 -1.836 0.482) 2.012 3.350 ( -0.426 -2.056 -1.016) 2.333 3.371 ( 0.268 -6.713 0.546) 6.741 3.395 ( -2.018 -7.671 -1.105) 8.008 3.777 ( 0.827 1.184 0.601) 1.564 3.789 ( -0.213 0.851 -0.443) 0.983 3.794 ( 0.585 -0.171 -0.395) 0.727 3.804 ( -0.333 -0.301 -1.237) 1.315 3.881 ( 0.349 2.050 0.876) 2.256 3.892 ( -0.408 1.872 -0.052) 1.916 3.930 ( -0.455 -1.790 0.847) 2.031 3.941 ( -1.309 -2.039 -0.153) 2.428 6.647 ( 0.357 0.474 0.275) 0.654 6.655 ( -0.262 0.427 -0.437) 0.664 6.672 ( -0.345 -1.554 0.364) 1.633 6.683 ( -1.219 -1.884 -0.615) 2.327 7.100 ( 1.299 0.351 0.755) 1.543 7.107 ( 1.043 0.203 0.128) 1.070 7.127 ( -0.529 -0.815 1.145) 1.502 7.147 ( -1.651 -2.324 -0.487) 2.892 7.510 ( 1.840 0.094 -0.136) 1.847 7.518 ( 1.351 -0.649 -0.775) 1.687 7.545 ( -1.456 0.096 -0.158) 1.468 7.546 ( -1.321 0.052 -0.249) 1.345 8.008 ( 2.613 2.596 0.589) 3.730 8.025 ( 2.326 1.148 -0.919) 2.752 8.086 ( -2.572 -0.170 1.156) 2.825 8.107 ( -3.219 -2.005 -0.695) 3.856 8.441 ( -0.283 6.060 0.387) 6.079 8.450 ( -1.156 5.898 -0.329) 6.019 8.573 ( -0.070 -6.351 0.207) 6.355 8.577 ( -0.464 -6.095 -0.119) 6.114 =================== 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/10080 10.0 3.026 4.962 3.027 0.000 0.000 -0.000 3/10080 20.0 0.751 1.314 0.756 0.000 0.000 -0.000 3/10080 30.0 0.413 0.743 0.411 0.000 0.000 -0.000 3/10080 40.0 0.289 0.533 0.285 0.000 0.000 -0.000 3/10080 50.0 0.225 0.424 0.220 0.000 0.000 -0.000 3/10080 60.0 0.186 0.357 0.180 0.000 0.000 -0.000 3/10080 70.0 0.160 0.311 0.154 0.000 0.000 -0.000 3/10080 80.0 0.141 0.278 0.134 0.000 0.000 -0.000 3/10080 90.0 0.127 0.253 0.120 0.000 0.000 -0.000 3/10080 100.0 0.116 0.233 0.108 0.000 0.000 -0.000 3/10080 110.0 0.107 0.216 0.099 0.000 0.000 -0.000 3/10080 120.0 0.100 0.202 0.091 0.000 0.000 -0.000 3/10080 130.0 0.093 0.190 0.085 0.000 0.000 -0.000 3/10080 140.0 0.088 0.179 0.079 0.000 0.000 -0.000 3/10080 150.0 0.083 0.169 0.074 0.000 0.000 -0.000 3/10080 160.0 0.078 0.161 0.070 0.000 0.000 -0.000 3/10080 170.0 0.074 0.153 0.066 0.000 0.000 -0.000 3/10080 180.0 0.071 0.146 0.063 0.000 0.000 -0.000 3/10080 190.0 0.067 0.140 0.060 0.000 0.000 -0.000 3/10080 200.0 0.064 0.134 0.057 0.000 0.000 -0.000 3/10080 210.0 0.062 0.128 0.054 0.000 0.000 -0.000 3/10080 220.0 0.059 0.123 0.052 0.000 0.000 -0.000 3/10080 230.0 0.057 0.119 0.050 0.000 0.000 -0.000 3/10080 240.0 0.055 0.114 0.048 0.000 0.000 -0.000 3/10080 250.0 0.053 0.110 0.046 0.000 0.000 -0.000 3/10080 260.0 0.051 0.107 0.044 0.000 0.000 -0.000 3/10080 270.0 0.049 0.103 0.043 0.000 0.000 -0.000 3/10080 280.0 0.048 0.100 0.041 0.000 0.000 -0.000 3/10080 290.0 0.046 0.097 0.040 0.000 0.000 -0.000 3/10080 300.0 0.045 0.094 0.039 0.000 0.000 -0.000 3/10080 310.0 0.043 0.091 0.038 0.000 0.000 -0.000 3/10080 320.0 0.042 0.089 0.036 0.000 0.000 -0.000 3/10080 330.0 0.041 0.086 0.035 0.000 0.000 -0.000 3/10080 340.0 0.040 0.084 0.034 0.000 0.000 -0.000 3/10080 350.0 0.039 0.082 0.033 0.000 0.000 -0.000 3/10080 360.0 0.038 0.080 0.033 0.000 0.000 -0.000 3/10080 370.0 0.037 0.078 0.032 0.000 0.000 -0.000 3/10080 380.0 0.036 0.076 0.031 0.000 0.000 -0.000 3/10080 390.0 0.035 0.074 0.030 0.000 0.000 -0.000 3/10080 400.0 0.034 0.072 0.029 0.000 0.000 -0.000 3/10080 410.0 0.034 0.071 0.029 0.000 0.000 -0.000 3/10080 420.0 0.033 0.069 0.028 0.000 0.000 -0.000 3/10080 430.0 0.032 0.068 0.027 0.000 0.000 -0.000 3/10080 440.0 0.031 0.066 0.027 0.000 0.000 -0.000 3/10080 450.0 0.031 0.065 0.026 0.000 0.000 -0.000 3/10080 460.0 0.030 0.063 0.026 0.000 0.000 -0.000 3/10080 470.0 0.029 0.062 0.025 0.000 0.000 -0.000 3/10080 480.0 0.029 0.061 0.025 0.000 0.000 -0.000 3/10080 490.0 0.028 0.060 0.024 0.000 0.000 -0.000 3/10080 500.0 0.028 0.059 0.024 0.000 0.000 -0.000 3/10080 510.0 0.027 0.058 0.023 0.000 0.000 -0.000 3/10080 520.0 0.027 0.056 0.023 0.000 0.000 -0.000 3/10080 530.0 0.026 0.055 0.022 0.000 0.000 -0.000 3/10080 540.0 0.026 0.054 0.022 0.000 0.000 -0.000 3/10080 550.0 0.025 0.054 0.022 0.000 0.000 -0.000 3/10080 560.0 0.025 0.053 0.021 0.000 0.000 -0.000 3/10080 570.0 0.025 0.052 0.021 0.000 0.000 -0.000 3/10080 580.0 0.024 0.051 0.020 0.000 0.000 -0.000 3/10080 590.0 0.024 0.050 0.020 0.000 0.000 -0.000 3/10080 600.0 0.023 0.049 0.020 0.000 0.000 -0.000 3/10080 610.0 0.023 0.049 0.019 0.000 0.000 -0.000 3/10080 620.0 0.023 0.048 0.019 0.000 0.000 -0.000 3/10080 630.0 0.022 0.047 0.019 0.000 0.000 -0.000 3/10080 640.0 0.022 0.046 0.019 0.000 0.000 -0.000 3/10080 650.0 0.022 0.046 0.018 0.000 0.000 -0.000 3/10080 660.0 0.021 0.045 0.018 0.000 0.000 -0.000 3/10080 670.0 0.021 0.044 0.018 0.000 0.000 -0.000 3/10080 680.0 0.021 0.044 0.017 0.000 0.000 -0.000 3/10080 690.0 0.020 0.043 0.017 0.000 0.000 -0.000 3/10080 700.0 0.020 0.043 0.017 0.000 0.000 -0.000 3/10080 710.0 0.020 0.042 0.017 0.000 0.000 -0.000 3/10080 720.0 0.020 0.041 0.017 0.000 0.000 -0.000 3/10080 730.0 0.019 0.041 0.016 0.000 0.000 -0.000 3/10080 740.0 0.019 0.040 0.016 0.000 0.000 -0.000 3/10080 750.0 0.019 0.040 0.016 0.000 0.000 -0.000 3/10080 760.0 0.019 0.039 0.016 0.000 0.000 -0.000 3/10080 770.0 0.018 0.039 0.015 0.000 0.000 -0.000 3/10080 780.0 0.018 0.038 0.015 0.000 0.000 -0.000 3/10080 790.0 0.018 0.038 0.015 0.000 0.000 -0.000 3/10080 800.0 0.018 0.037 0.015 0.000 0.000 -0.000 3/10080 810.0 0.017 0.037 0.015 0.000 0.000 -0.000 3/10080 820.0 0.017 0.036 0.015 0.000 0.000 -0.000 3/10080 830.0 0.017 0.036 0.014 0.000 0.000 -0.000 3/10080 840.0 0.017 0.036 0.014 0.000 0.000 -0.000 3/10080 850.0 0.017 0.035 0.014 0.000 0.000 -0.000 3/10080 860.0 0.016 0.035 0.014 0.000 0.000 -0.000 3/10080 870.0 0.016 0.034 0.014 0.000 0.000 -0.000 3/10080 880.0 0.016 0.034 0.014 0.000 0.000 -0.000 3/10080 890.0 0.016 0.034 0.013 0.000 0.000 -0.000 3/10080 900.0 0.016 0.033 0.013 0.000 0.000 -0.000 3/10080 910.0 0.016 0.033 0.013 0.000 0.000 -0.000 3/10080 920.0 0.015 0.033 0.013 0.000 0.000 -0.000 3/10080 930.0 0.015 0.032 0.013 0.000 0.000 -0.000 3/10080 940.0 0.015 0.032 0.013 0.000 0.000 -0.000 3/10080 950.0 0.015 0.032 0.013 0.000 0.000 -0.000 3/10080 960.0 0.015 0.031 0.012 0.000 0.000 -0.000 3/10080 970.0 0.015 0.031 0.012 0.000 0.000 -0.000 3/10080 980.0 0.014 0.031 0.012 0.000 0.000 -0.000 3/10080 990.0 0.014 0.030 0.012 0.000 0.000 -0.000 3/10080 1000.0 0.014 0.030 0.012 0.000 0.000 -0.000 3/10080 Thermal conductivity related properties were written into "kappa-m573.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 10:02:20]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|