----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 02:21:34]------------------------- 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: [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) Number of symmetry operations in supercell: 384 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.447159740204642 0.000000000000000 0.000000000000000 b -2.223579870102321 3.851353309704624 0.000000000000000 c 0.000000000000000 0.000000000000000 16.881745339999998 Atomic positions (fractional): *1 Rb 0.33333333333333 0.66666666666667 0.25000000000000 85.468 2 Rb 0.66666666666667 0.33333333333333 0.75000000000000 85.468 *3 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 4 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 *5 Te 0.33333333333333 0.66666666666667 0.60024677538198 127.600 6 Te 0.66666666666667 0.33333333333333 0.10024677538198 127.600 7 Te 0.66666666666667 0.33333333333333 0.39975322461802 127.600 8 Te 0.33333333333333 0.66666666666667 0.89975322461802 127.600 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.447159740204642 0.000000000000000 0.000000000000000 b -2.223579870102321 3.851353309704624 0.000000000000000 c 0.000000000000000 0.000000000000000 16.881745339999998 Atomic positions (fractional): *1 Rb 0.33333333333333 0.66666666666667 0.25000000000000 85.468 > 1 2 Rb 0.66666666666667 0.33333333333333 0.75000000000000 85.468 > 2 *3 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 > 3 4 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 > 4 *5 Te 0.33333333333333 0.66666666666667 0.60024677538198 127.600 > 5 6 Te 0.66666666666667 0.33333333333333 0.10024677538198 127.600 > 6 7 Te 0.66666666666667 0.33333333333333 0.39975322461802 127.600 > 7 8 Te 0.33333333333333 0.66666666666667 0.89975322461802 127.600 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 17.788638960818567 0.000000000000000 0.000000000000000 b -8.894319480409283 15.405413238818497 0.000000000000000 c 0.000000000000000 0.000000000000000 16.881745339999998 Atomic positions (fractional): *1 Rb 0.08333333333333 0.16666666666667 0.25000000000000 85.468 > 1 2 Rb 0.33333333333333 0.16666666666667 0.25000000000000 85.468 > 1 3 Rb 0.58333333333333 0.16666666666667 0.25000000000000 85.468 > 1 4 Rb 0.83333333333333 0.16666666666667 0.25000000000000 85.468 > 1 5 Rb 0.08333333333333 0.41666666666667 0.25000000000000 85.468 > 1 6 Rb 0.33333333333333 0.41666666666667 0.25000000000000 85.468 > 1 7 Rb 0.58333333333333 0.41666666666667 0.25000000000000 85.468 > 1 8 Rb 0.83333333333333 0.41666666666667 0.25000000000000 85.468 > 1 9 Rb 0.08333333333333 0.66666666666667 0.25000000000000 85.468 > 1 10 Rb 0.33333333333333 0.66666666666667 0.25000000000000 85.468 > 1 11 Rb 0.58333333333333 0.66666666666667 0.25000000000000 85.468 > 1 12 Rb 0.83333333333333 0.66666666666667 0.25000000000000 85.468 > 1 13 Rb 0.08333333333333 0.91666666666667 0.25000000000000 85.468 > 1 14 Rb 0.33333333333333 0.91666666666667 0.25000000000000 85.468 > 1 15 Rb 0.58333333333333 0.91666666666667 0.25000000000000 85.468 > 1 16 Rb 0.83333333333333 0.91666666666667 0.25000000000000 85.468 > 1 17 Rb 0.16666666666667 0.08333333333333 0.75000000000000 85.468 > 2 18 Rb 0.41666666666667 0.08333333333333 0.75000000000000 85.468 > 2 19 Rb 0.66666666666667 0.08333333333333 0.75000000000000 85.468 > 2 20 Rb 0.91666666666667 0.08333333333333 0.75000000000000 85.468 > 2 21 Rb 0.16666666666667 0.33333333333333 0.75000000000000 85.468 > 2 22 Rb 0.41666666666667 0.33333333333333 0.75000000000000 85.468 > 2 23 Rb 0.66666666666667 0.33333333333333 0.75000000000000 85.468 > 2 24 Rb 0.91666666666667 0.33333333333333 0.75000000000000 85.468 > 2 25 Rb 0.16666666666667 0.58333333333333 0.75000000000000 85.468 > 2 26 Rb 0.41666666666667 0.58333333333333 0.75000000000000 85.468 > 2 27 Rb 0.66666666666667 0.58333333333333 0.75000000000000 85.468 > 2 28 Rb 0.91666666666667 0.58333333333333 0.75000000000000 85.468 > 2 29 Rb 0.16666666666667 0.83333333333333 0.75000000000000 85.468 > 2 30 Rb 0.41666666666667 0.83333333333333 0.75000000000000 85.468 > 2 31 Rb 0.66666666666667 0.83333333333333 0.75000000000000 85.468 > 2 32 Rb 0.91666666666667 0.83333333333333 0.75000000000000 85.468 > 2 *33 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 > 3 34 Y 0.25000000000000 0.00000000000000 0.00000000000000 88.906 > 3 35 Y 0.50000000000000 0.00000000000000 0.00000000000000 88.906 > 3 36 Y 0.75000000000000 0.00000000000000 0.00000000000000 88.906 > 3 37 Y 0.00000000000000 0.25000000000000 0.00000000000000 88.906 > 3 38 Y 0.25000000000000 0.25000000000000 0.00000000000000 88.906 > 3 39 Y 0.50000000000000 0.25000000000000 0.00000000000000 88.906 > 3 40 Y 0.75000000000000 0.25000000000000 0.00000000000000 88.906 > 3 41 Y 0.00000000000000 0.50000000000000 0.00000000000000 88.906 > 3 42 Y 0.25000000000000 0.50000000000000 0.00000000000000 88.906 > 3 43 Y 0.50000000000000 0.50000000000000 0.00000000000000 88.906 > 3 44 Y 0.75000000000000 0.50000000000000 0.00000000000000 88.906 > 3 45 Y 0.00000000000000 0.75000000000000 0.00000000000000 88.906 > 3 46 Y 0.25000000000000 0.75000000000000 0.00000000000000 88.906 > 3 47 Y 0.50000000000000 0.75000000000000 0.00000000000000 88.906 > 3 48 Y 0.75000000000000 0.75000000000000 0.00000000000000 88.906 > 3 49 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 > 4 50 Y 0.25000000000000 0.00000000000000 0.50000000000000 88.906 > 4 51 Y 0.50000000000000 0.00000000000000 0.50000000000000 88.906 > 4 52 Y 0.75000000000000 0.00000000000000 0.50000000000000 88.906 > 4 53 Y 0.00000000000000 0.25000000000000 0.50000000000000 88.906 > 4 54 Y 0.25000000000000 0.25000000000000 0.50000000000000 88.906 > 4 55 Y 0.50000000000000 0.25000000000000 0.50000000000000 88.906 > 4 56 Y 0.75000000000000 0.25000000000000 0.50000000000000 88.906 > 4 57 Y 0.00000000000000 0.50000000000000 0.50000000000000 88.906 > 4 58 Y 0.25000000000000 0.50000000000000 0.50000000000000 88.906 > 4 59 Y 0.50000000000000 0.50000000000000 0.50000000000000 88.906 > 4 60 Y 0.75000000000000 0.50000000000000 0.50000000000000 88.906 > 4 61 Y 0.00000000000000 0.75000000000000 0.50000000000000 88.906 > 4 62 Y 0.25000000000000 0.75000000000000 0.50000000000000 88.906 > 4 63 Y 0.50000000000000 0.75000000000000 0.50000000000000 88.906 > 4 64 Y 0.75000000000000 0.75000000000000 0.50000000000000 88.906 > 4 *65 Te 0.08333333333333 0.16666666666667 0.60024677538198 127.600 > 5 66 Te 0.33333333333333 0.16666666666667 0.60024677538198 127.600 > 5 67 Te 0.58333333333333 0.16666666666667 0.60024677538198 127.600 > 5 68 Te 0.83333333333333 0.16666666666667 0.60024677538198 127.600 > 5 69 Te 0.08333333333333 0.41666666666667 0.60024677538198 127.600 > 5 70 Te 0.33333333333333 0.41666666666667 0.60024677538198 127.600 > 5 71 Te 0.58333333333333 0.41666666666667 0.60024677538198 127.600 > 5 72 Te 0.83333333333333 0.41666666666667 0.60024677538198 127.600 > 5 73 Te 0.08333333333333 0.66666666666667 0.60024677538198 127.600 > 5 74 Te 0.33333333333333 0.66666666666667 0.60024677538198 127.600 > 5 75 Te 0.58333333333333 0.66666666666667 0.60024677538198 127.600 > 5 76 Te 0.83333333333333 0.66666666666667 0.60024677538198 127.600 > 5 77 Te 0.08333333333333 0.91666666666667 0.60024677538198 127.600 > 5 78 Te 0.33333333333333 0.91666666666667 0.60024677538198 127.600 > 5 79 Te 0.58333333333333 0.91666666666667 0.60024677538198 127.600 > 5 80 Te 0.83333333333333 0.91666666666667 0.60024677538198 127.600 > 5 81 Te 0.16666666666667 0.08333333333333 0.10024677538198 127.600 > 6 82 Te 0.41666666666667 0.08333333333333 0.10024677538198 127.600 > 6 83 Te 0.66666666666667 0.08333333333333 0.10024677538198 127.600 > 6 84 Te 0.91666666666667 0.08333333333333 0.10024677538198 127.600 > 6 85 Te 0.16666666666667 0.33333333333333 0.10024677538198 127.600 > 6 86 Te 0.41666666666667 0.33333333333333 0.10024677538198 127.600 > 6 87 Te 0.66666666666667 0.33333333333333 0.10024677538198 127.600 > 6 88 Te 0.91666666666667 0.33333333333333 0.10024677538198 127.600 > 6 89 Te 0.16666666666667 0.58333333333333 0.10024677538198 127.600 > 6 90 Te 0.41666666666667 0.58333333333333 0.10024677538198 127.600 > 6 91 Te 0.66666666666667 0.58333333333333 0.10024677538198 127.600 > 6 92 Te 0.91666666666667 0.58333333333333 0.10024677538198 127.600 > 6 93 Te 0.16666666666667 0.83333333333333 0.10024677538198 127.600 > 6 94 Te 0.41666666666667 0.83333333333333 0.10024677538198 127.600 > 6 95 Te 0.66666666666667 0.83333333333333 0.10024677538198 127.600 > 6 96 Te 0.91666666666667 0.83333333333333 0.10024677538198 127.600 > 6 97 Te 0.16666666666667 0.08333333333333 0.39975322461802 127.600 > 7 98 Te 0.41666666666667 0.08333333333333 0.39975322461802 127.600 > 7 99 Te 0.66666666666667 0.08333333333333 0.39975322461802 127.600 > 7 100 Te 0.91666666666667 0.08333333333333 0.39975322461802 127.600 > 7 101 Te 0.16666666666667 0.33333333333333 0.39975322461802 127.600 > 7 102 Te 0.41666666666667 0.33333333333333 0.39975322461802 127.600 > 7 103 Te 0.66666666666667 0.33333333333333 0.39975322461802 127.600 > 7 104 Te 0.91666666666667 0.33333333333333 0.39975322461802 127.600 > 7 105 Te 0.16666666666667 0.58333333333333 0.39975322461802 127.600 > 7 106 Te 0.41666666666667 0.58333333333333 0.39975322461802 127.600 > 7 107 Te 0.66666666666667 0.58333333333333 0.39975322461802 127.600 > 7 108 Te 0.91666666666667 0.58333333333333 0.39975322461802 127.600 > 7 109 Te 0.16666666666667 0.83333333333333 0.39975322461802 127.600 > 7 110 Te 0.41666666666667 0.83333333333333 0.39975322461802 127.600 > 7 111 Te 0.66666666666667 0.83333333333333 0.39975322461802 127.600 > 7 112 Te 0.91666666666667 0.83333333333333 0.39975322461802 127.600 > 7 113 Te 0.08333333333333 0.16666666666667 0.89975322461802 127.600 > 8 114 Te 0.33333333333333 0.16666666666667 0.89975322461802 127.600 > 8 115 Te 0.58333333333333 0.16666666666667 0.89975322461802 127.600 > 8 116 Te 0.83333333333333 0.16666666666667 0.89975322461802 127.600 > 8 117 Te 0.08333333333333 0.41666666666667 0.89975322461802 127.600 > 8 118 Te 0.33333333333333 0.41666666666667 0.89975322461802 127.600 > 8 119 Te 0.58333333333333 0.41666666666667 0.89975322461802 127.600 > 8 120 Te 0.83333333333333 0.41666666666667 0.89975322461802 127.600 > 8 121 Te 0.08333333333333 0.66666666666667 0.89975322461802 127.600 > 8 122 Te 0.33333333333333 0.66666666666667 0.89975322461802 127.600 > 8 123 Te 0.58333333333333 0.66666666666667 0.89975322461802 127.600 > 8 124 Te 0.83333333333333 0.66666666666667 0.89975322461802 127.600 > 8 125 Te 0.08333333333333 0.91666666666667 0.89975322461802 127.600 > 8 126 Te 0.33333333333333 0.91666666666667 0.89975322461802 127.600 > 8 127 Te 0.58333333333333 0.91666666666667 0.89975322461802 127.600 > 8 128 Te 0.83333333333333 0.91666666666667 0.89975322461802 127.600 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 8.5031623 0.0000000 0.0000000 -0.0000000 8.5031623 0.0000000 0.0000000 0.0000000 6.0008867 -------------------------- Born effective charges -------------------------- 1 Rb 1.1052389 0.0000000 0.0000000 0.0000000 1.1052389 0.0000000 0.0000000 0.0000000 2.1323423 2 Rb 1.1052389 0.0000000 0.0000000 0.0000000 1.1052389 0.0000000 0.0000000 0.0000000 2.1323423 3 Y 4.8553139 0.0000000 0.0000000 0.0000000 4.8553139 0.0000000 0.0000000 0.0000000 1.8077586 4 Y 4.8553139 0.0000000 0.0000000 0.0000000 4.8553139 0.0000000 0.0000000 0.0000000 1.8077586 5 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 6 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 7 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 8 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 ---------------------------------------------------------------------------- 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: 5136/5136 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 259 Number of blocks in projector: 213 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 123 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 76 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 60 Use standard eigh solver. Tree of FC basis block matrices: - (259, 253), data: False |-- (60, 57), data: True |-- (76, 76), data: True |-- (123, 120), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.009 Solver_block: 80 / 80 - Time: 0.199 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.210 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (yy) 0.00000000 (yy) Permutation basis: 384/384 Permutation basis: 5136/5136 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 259 Number of blocks in projector: 213 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 123 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 76 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 60 Use standard eigh solver. Tree of FC basis block matrices: - (259, 253), data: False |-- (60, 57), data: True |-- (76, 76), data: True |-- (123, 120), data: True Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 02:21:40]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 02:21:40]------------------------- 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): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.447159740204642 0.000000000000000 0.000000000000000 b -2.223579870102321 3.851353309704624 0.000000000000000 c 0.000000000000000 0.000000000000000 16.881745339999998 Atomic positions (fractional): 1 Rb 0.33333333333333 0.66666666666667 0.25000000000000 85.468 2 Rb 0.66666666666667 0.33333333333333 0.75000000000000 85.468 3 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 4 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 5 Te 0.33333333333333 0.66666666666667 0.60024677538198 127.600 6 Te 0.66666666666667 0.33333333333333 0.10024677538198 127.600 7 Te 0.66666666666667 0.33333333333333 0.39975322461802 127.600 8 Te 0.33333333333333 0.66666666666667 0.89975322461802 127.600 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.788638960818567 0.000000000000000 0.000000000000000 b -8.894319480409283 15.405413238818497 0.000000000000000 c 0.000000000000000 0.000000000000000 16.881745339999998 Atomic positions (fractional): 1 Rb 0.08333333333333 0.16666666666667 0.25000000000000 85.468 > 1 2 Rb 0.33333333333333 0.16666666666667 0.25000000000000 85.468 > 1 3 Rb 0.58333333333333 0.16666666666667 0.25000000000000 85.468 > 1 4 Rb 0.83333333333333 0.16666666666667 0.25000000000000 85.468 > 1 5 Rb 0.08333333333333 0.41666666666667 0.25000000000000 85.468 > 1 6 Rb 0.33333333333333 0.41666666666667 0.25000000000000 85.468 > 1 7 Rb 0.58333333333333 0.41666666666667 0.25000000000000 85.468 > 1 8 Rb 0.83333333333333 0.41666666666667 0.25000000000000 85.468 > 1 9 Rb 0.08333333333333 0.66666666666667 0.25000000000000 85.468 > 1 10 Rb 0.33333333333333 0.66666666666667 0.25000000000000 85.468 > 1 11 Rb 0.58333333333333 0.66666666666667 0.25000000000000 85.468 > 1 12 Rb 0.83333333333333 0.66666666666667 0.25000000000000 85.468 > 1 13 Rb 0.08333333333333 0.91666666666667 0.25000000000000 85.468 > 1 14 Rb 0.33333333333333 0.91666666666667 0.25000000000000 85.468 > 1 15 Rb 0.58333333333333 0.91666666666667 0.25000000000000 85.468 > 1 16 Rb 0.83333333333333 0.91666666666667 0.25000000000000 85.468 > 1 17 Rb 0.16666666666667 0.08333333333333 0.75000000000000 85.468 > 17 18 Rb 0.41666666666667 0.08333333333333 0.75000000000000 85.468 > 17 19 Rb 0.66666666666667 0.08333333333333 0.75000000000000 85.468 > 17 20 Rb 0.91666666666667 0.08333333333333 0.75000000000000 85.468 > 17 21 Rb 0.16666666666667 0.33333333333333 0.75000000000000 85.468 > 17 22 Rb 0.41666666666667 0.33333333333333 0.75000000000000 85.468 > 17 23 Rb 0.66666666666667 0.33333333333333 0.75000000000000 85.468 > 17 24 Rb 0.91666666666667 0.33333333333333 0.75000000000000 85.468 > 17 25 Rb 0.16666666666667 0.58333333333333 0.75000000000000 85.468 > 17 26 Rb 0.41666666666667 0.58333333333333 0.75000000000000 85.468 > 17 27 Rb 0.66666666666667 0.58333333333333 0.75000000000000 85.468 > 17 28 Rb 0.91666666666667 0.58333333333333 0.75000000000000 85.468 > 17 29 Rb 0.16666666666667 0.83333333333333 0.75000000000000 85.468 > 17 30 Rb 0.41666666666667 0.83333333333333 0.75000000000000 85.468 > 17 31 Rb 0.66666666666667 0.83333333333333 0.75000000000000 85.468 > 17 32 Rb 0.91666666666667 0.83333333333333 0.75000000000000 85.468 > 17 33 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 > 33 34 Y 0.25000000000000 0.00000000000000 0.00000000000000 88.906 > 33 35 Y 0.50000000000000 0.00000000000000 0.00000000000000 88.906 > 33 36 Y 0.75000000000000 0.00000000000000 0.00000000000000 88.906 > 33 37 Y 0.00000000000000 0.25000000000000 0.00000000000000 88.906 > 33 38 Y 0.25000000000000 0.25000000000000 0.00000000000000 88.906 > 33 39 Y 0.50000000000000 0.25000000000000 0.00000000000000 88.906 > 33 40 Y 0.75000000000000 0.25000000000000 0.00000000000000 88.906 > 33 41 Y 0.00000000000000 0.50000000000000 0.00000000000000 88.906 > 33 42 Y 0.25000000000000 0.50000000000000 0.00000000000000 88.906 > 33 43 Y 0.50000000000000 0.50000000000000 0.00000000000000 88.906 > 33 44 Y 0.75000000000000 0.50000000000000 0.00000000000000 88.906 > 33 45 Y 0.00000000000000 0.75000000000000 0.00000000000000 88.906 > 33 46 Y 0.25000000000000 0.75000000000000 0.00000000000000 88.906 > 33 47 Y 0.50000000000000 0.75000000000000 0.00000000000000 88.906 > 33 48 Y 0.75000000000000 0.75000000000000 0.00000000000000 88.906 > 33 49 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 > 49 50 Y 0.25000000000000 0.00000000000000 0.50000000000000 88.906 > 49 51 Y 0.50000000000000 0.00000000000000 0.50000000000000 88.906 > 49 52 Y 0.75000000000000 0.00000000000000 0.50000000000000 88.906 > 49 53 Y 0.00000000000000 0.25000000000000 0.50000000000000 88.906 > 49 54 Y 0.25000000000000 0.25000000000000 0.50000000000000 88.906 > 49 55 Y 0.50000000000000 0.25000000000000 0.50000000000000 88.906 > 49 56 Y 0.75000000000000 0.25000000000000 0.50000000000000 88.906 > 49 57 Y 0.00000000000000 0.50000000000000 0.50000000000000 88.906 > 49 58 Y 0.25000000000000 0.50000000000000 0.50000000000000 88.906 > 49 59 Y 0.50000000000000 0.50000000000000 0.50000000000000 88.906 > 49 60 Y 0.75000000000000 0.50000000000000 0.50000000000000 88.906 > 49 61 Y 0.00000000000000 0.75000000000000 0.50000000000000 88.906 > 49 62 Y 0.25000000000000 0.75000000000000 0.50000000000000 88.906 > 49 63 Y 0.50000000000000 0.75000000000000 0.50000000000000 88.906 > 49 64 Y 0.75000000000000 0.75000000000000 0.50000000000000 88.906 > 49 65 Te 0.08333333333333 0.16666666666667 0.60024677538198 127.600 > 65 66 Te 0.33333333333333 0.16666666666667 0.60024677538198 127.600 > 65 67 Te 0.58333333333333 0.16666666666667 0.60024677538198 127.600 > 65 68 Te 0.83333333333333 0.16666666666667 0.60024677538198 127.600 > 65 69 Te 0.08333333333333 0.41666666666667 0.60024677538198 127.600 > 65 70 Te 0.33333333333333 0.41666666666667 0.60024677538198 127.600 > 65 71 Te 0.58333333333333 0.41666666666667 0.60024677538198 127.600 > 65 72 Te 0.83333333333333 0.41666666666667 0.60024677538198 127.600 > 65 73 Te 0.08333333333333 0.66666666666667 0.60024677538198 127.600 > 65 74 Te 0.33333333333333 0.66666666666667 0.60024677538198 127.600 > 65 75 Te 0.58333333333333 0.66666666666667 0.60024677538198 127.600 > 65 76 Te 0.83333333333333 0.66666666666667 0.60024677538198 127.600 > 65 77 Te 0.08333333333333 0.91666666666667 0.60024677538198 127.600 > 65 78 Te 0.33333333333333 0.91666666666667 0.60024677538198 127.600 > 65 79 Te 0.58333333333333 0.91666666666667 0.60024677538198 127.600 > 65 80 Te 0.83333333333333 0.91666666666667 0.60024677538198 127.600 > 65 81 Te 0.16666666666667 0.08333333333333 0.10024677538198 127.600 > 81 82 Te 0.41666666666667 0.08333333333333 0.10024677538198 127.600 > 81 83 Te 0.66666666666667 0.08333333333333 0.10024677538198 127.600 > 81 84 Te 0.91666666666667 0.08333333333333 0.10024677538198 127.600 > 81 85 Te 0.16666666666667 0.33333333333333 0.10024677538198 127.600 > 81 86 Te 0.41666666666667 0.33333333333333 0.10024677538198 127.600 > 81 87 Te 0.66666666666667 0.33333333333333 0.10024677538198 127.600 > 81 88 Te 0.91666666666667 0.33333333333333 0.10024677538198 127.600 > 81 89 Te 0.16666666666667 0.58333333333333 0.10024677538198 127.600 > 81 90 Te 0.41666666666667 0.58333333333333 0.10024677538198 127.600 > 81 91 Te 0.66666666666667 0.58333333333333 0.10024677538198 127.600 > 81 92 Te 0.91666666666667 0.58333333333333 0.10024677538198 127.600 > 81 93 Te 0.16666666666667 0.83333333333333 0.10024677538198 127.600 > 81 94 Te 0.41666666666667 0.83333333333333 0.10024677538198 127.600 > 81 95 Te 0.66666666666667 0.83333333333333 0.10024677538198 127.600 > 81 96 Te 0.91666666666667 0.83333333333333 0.10024677538198 127.600 > 81 97 Te 0.16666666666667 0.08333333333333 0.39975322461802 127.600 > 97 98 Te 0.41666666666667 0.08333333333333 0.39975322461802 127.600 > 97 99 Te 0.66666666666667 0.08333333333333 0.39975322461802 127.600 > 97 100 Te 0.91666666666667 0.08333333333333 0.39975322461802 127.600 > 97 101 Te 0.16666666666667 0.33333333333333 0.39975322461802 127.600 > 97 102 Te 0.41666666666667 0.33333333333333 0.39975322461802 127.600 > 97 103 Te 0.66666666666667 0.33333333333333 0.39975322461802 127.600 > 97 104 Te 0.91666666666667 0.33333333333333 0.39975322461802 127.600 > 97 105 Te 0.16666666666667 0.58333333333333 0.39975322461802 127.600 > 97 106 Te 0.41666666666667 0.58333333333333 0.39975322461802 127.600 > 97 107 Te 0.66666666666667 0.58333333333333 0.39975322461802 127.600 > 97 108 Te 0.91666666666667 0.58333333333333 0.39975322461802 127.600 > 97 109 Te 0.16666666666667 0.83333333333333 0.39975322461802 127.600 > 97 110 Te 0.41666666666667 0.83333333333333 0.39975322461802 127.600 > 97 111 Te 0.66666666666667 0.83333333333333 0.39975322461802 127.600 > 97 112 Te 0.91666666666667 0.83333333333333 0.39975322461802 127.600 > 97 113 Te 0.08333333333333 0.16666666666667 0.89975322461802 127.600 > 113 114 Te 0.33333333333333 0.16666666666667 0.89975322461802 127.600 > 113 115 Te 0.58333333333333 0.16666666666667 0.89975322461802 127.600 > 113 116 Te 0.83333333333333 0.16666666666667 0.89975322461802 127.600 > 113 117 Te 0.08333333333333 0.41666666666667 0.89975322461802 127.600 > 113 118 Te 0.33333333333333 0.41666666666667 0.89975322461802 127.600 > 113 119 Te 0.58333333333333 0.41666666666667 0.89975322461802 127.600 > 113 120 Te 0.83333333333333 0.41666666666667 0.89975322461802 127.600 > 113 121 Te 0.08333333333333 0.66666666666667 0.89975322461802 127.600 > 113 122 Te 0.33333333333333 0.66666666666667 0.89975322461802 127.600 > 113 123 Te 0.58333333333333 0.66666666666667 0.89975322461802 127.600 > 113 124 Te 0.83333333333333 0.66666666666667 0.89975322461802 127.600 > 113 125 Te 0.08333333333333 0.91666666666667 0.89975322461802 127.600 > 113 126 Te 0.33333333333333 0.91666666666667 0.89975322461802 127.600 > 113 127 Te 0.58333333333333 0.91666666666667 0.89975322461802 127.600 > 113 128 Te 0.83333333333333 0.91666666666667 0.89975322461802 127.600 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 8.5031623 0.0000000 0.0000000 -0.0000000 8.5031623 0.0000000 0.0000000 0.0000000 6.0008867 -------------------------- Born effective charges -------------------------- 1 Rb 1.1052389 0.0000000 0.0000000 0.0000000 1.1052389 0.0000000 0.0000000 0.0000000 2.1323423 2 Rb 1.1052389 0.0000000 0.0000000 0.0000000 1.1052389 0.0000000 0.0000000 0.0000000 2.1323423 3 Y 4.8553139 0.0000000 0.0000000 0.0000000 4.8553139 0.0000000 0.0000000 0.0000000 1.8077586 4 Y 4.8553139 0.0000000 0.0000000 0.0000000 4.8553139 0.0000000 0.0000000 0.0000000 1.8077586 5 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 6 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 7 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 8 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 ---------------------------------------------------------------------------- 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.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.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.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.00000004 (yyy) 0.00000004 (yyy) 0.00000004 (yyy) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 02:21:49]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 02:21:49]------------------------- 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): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.447159740204642 0.000000000000000 0.000000000000000 b -2.223579870102321 3.851353309704624 0.000000000000000 c 0.000000000000000 0.000000000000000 16.881745339999998 Atomic positions (fractional): 1 Rb 0.33333333333333 0.66666666666667 0.25000000000000 85.468 2 Rb 0.66666666666667 0.33333333333333 0.75000000000000 85.468 3 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 4 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 5 Te 0.33333333333333 0.66666666666667 0.60024677538198 127.600 6 Te 0.66666666666667 0.33333333333333 0.10024677538198 127.600 7 Te 0.66666666666667 0.33333333333333 0.39975322461802 127.600 8 Te 0.33333333333333 0.66666666666667 0.89975322461802 127.600 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.788638960818567 0.000000000000000 0.000000000000000 b -8.894319480409283 15.405413238818497 0.000000000000000 c 0.000000000000000 0.000000000000000 16.881745339999998 Atomic positions (fractional): 1 Rb 0.08333333333333 0.16666666666667 0.25000000000000 85.468 > 1 2 Rb 0.33333333333333 0.16666666666667 0.25000000000000 85.468 > 1 3 Rb 0.58333333333333 0.16666666666667 0.25000000000000 85.468 > 1 4 Rb 0.83333333333333 0.16666666666667 0.25000000000000 85.468 > 1 5 Rb 0.08333333333333 0.41666666666667 0.25000000000000 85.468 > 1 6 Rb 0.33333333333333 0.41666666666667 0.25000000000000 85.468 > 1 7 Rb 0.58333333333333 0.41666666666667 0.25000000000000 85.468 > 1 8 Rb 0.83333333333333 0.41666666666667 0.25000000000000 85.468 > 1 9 Rb 0.08333333333333 0.66666666666667 0.25000000000000 85.468 > 1 10 Rb 0.33333333333333 0.66666666666667 0.25000000000000 85.468 > 1 11 Rb 0.58333333333333 0.66666666666667 0.25000000000000 85.468 > 1 12 Rb 0.83333333333333 0.66666666666667 0.25000000000000 85.468 > 1 13 Rb 0.08333333333333 0.91666666666667 0.25000000000000 85.468 > 1 14 Rb 0.33333333333333 0.91666666666667 0.25000000000000 85.468 > 1 15 Rb 0.58333333333333 0.91666666666667 0.25000000000000 85.468 > 1 16 Rb 0.83333333333333 0.91666666666667 0.25000000000000 85.468 > 1 17 Rb 0.16666666666667 0.08333333333333 0.75000000000000 85.468 > 17 18 Rb 0.41666666666667 0.08333333333333 0.75000000000000 85.468 > 17 19 Rb 0.66666666666667 0.08333333333333 0.75000000000000 85.468 > 17 20 Rb 0.91666666666667 0.08333333333333 0.75000000000000 85.468 > 17 21 Rb 0.16666666666667 0.33333333333333 0.75000000000000 85.468 > 17 22 Rb 0.41666666666667 0.33333333333333 0.75000000000000 85.468 > 17 23 Rb 0.66666666666667 0.33333333333333 0.75000000000000 85.468 > 17 24 Rb 0.91666666666667 0.33333333333333 0.75000000000000 85.468 > 17 25 Rb 0.16666666666667 0.58333333333333 0.75000000000000 85.468 > 17 26 Rb 0.41666666666667 0.58333333333333 0.75000000000000 85.468 > 17 27 Rb 0.66666666666667 0.58333333333333 0.75000000000000 85.468 > 17 28 Rb 0.91666666666667 0.58333333333333 0.75000000000000 85.468 > 17 29 Rb 0.16666666666667 0.83333333333333 0.75000000000000 85.468 > 17 30 Rb 0.41666666666667 0.83333333333333 0.75000000000000 85.468 > 17 31 Rb 0.66666666666667 0.83333333333333 0.75000000000000 85.468 > 17 32 Rb 0.91666666666667 0.83333333333333 0.75000000000000 85.468 > 17 33 Y 0.00000000000000 0.00000000000000 0.00000000000000 88.906 > 33 34 Y 0.25000000000000 0.00000000000000 0.00000000000000 88.906 > 33 35 Y 0.50000000000000 0.00000000000000 0.00000000000000 88.906 > 33 36 Y 0.75000000000000 0.00000000000000 0.00000000000000 88.906 > 33 37 Y 0.00000000000000 0.25000000000000 0.00000000000000 88.906 > 33 38 Y 0.25000000000000 0.25000000000000 0.00000000000000 88.906 > 33 39 Y 0.50000000000000 0.25000000000000 0.00000000000000 88.906 > 33 40 Y 0.75000000000000 0.25000000000000 0.00000000000000 88.906 > 33 41 Y 0.00000000000000 0.50000000000000 0.00000000000000 88.906 > 33 42 Y 0.25000000000000 0.50000000000000 0.00000000000000 88.906 > 33 43 Y 0.50000000000000 0.50000000000000 0.00000000000000 88.906 > 33 44 Y 0.75000000000000 0.50000000000000 0.00000000000000 88.906 > 33 45 Y 0.00000000000000 0.75000000000000 0.00000000000000 88.906 > 33 46 Y 0.25000000000000 0.75000000000000 0.00000000000000 88.906 > 33 47 Y 0.50000000000000 0.75000000000000 0.00000000000000 88.906 > 33 48 Y 0.75000000000000 0.75000000000000 0.00000000000000 88.906 > 33 49 Y 0.00000000000000 0.00000000000000 0.50000000000000 88.906 > 49 50 Y 0.25000000000000 0.00000000000000 0.50000000000000 88.906 > 49 51 Y 0.50000000000000 0.00000000000000 0.50000000000000 88.906 > 49 52 Y 0.75000000000000 0.00000000000000 0.50000000000000 88.906 > 49 53 Y 0.00000000000000 0.25000000000000 0.50000000000000 88.906 > 49 54 Y 0.25000000000000 0.25000000000000 0.50000000000000 88.906 > 49 55 Y 0.50000000000000 0.25000000000000 0.50000000000000 88.906 > 49 56 Y 0.75000000000000 0.25000000000000 0.50000000000000 88.906 > 49 57 Y 0.00000000000000 0.50000000000000 0.50000000000000 88.906 > 49 58 Y 0.25000000000000 0.50000000000000 0.50000000000000 88.906 > 49 59 Y 0.50000000000000 0.50000000000000 0.50000000000000 88.906 > 49 60 Y 0.75000000000000 0.50000000000000 0.50000000000000 88.906 > 49 61 Y 0.00000000000000 0.75000000000000 0.50000000000000 88.906 > 49 62 Y 0.25000000000000 0.75000000000000 0.50000000000000 88.906 > 49 63 Y 0.50000000000000 0.75000000000000 0.50000000000000 88.906 > 49 64 Y 0.75000000000000 0.75000000000000 0.50000000000000 88.906 > 49 65 Te 0.08333333333333 0.16666666666667 0.60024677538198 127.600 > 65 66 Te 0.33333333333333 0.16666666666667 0.60024677538198 127.600 > 65 67 Te 0.58333333333333 0.16666666666667 0.60024677538198 127.600 > 65 68 Te 0.83333333333333 0.16666666666667 0.60024677538198 127.600 > 65 69 Te 0.08333333333333 0.41666666666667 0.60024677538198 127.600 > 65 70 Te 0.33333333333333 0.41666666666667 0.60024677538198 127.600 > 65 71 Te 0.58333333333333 0.41666666666667 0.60024677538198 127.600 > 65 72 Te 0.83333333333333 0.41666666666667 0.60024677538198 127.600 > 65 73 Te 0.08333333333333 0.66666666666667 0.60024677538198 127.600 > 65 74 Te 0.33333333333333 0.66666666666667 0.60024677538198 127.600 > 65 75 Te 0.58333333333333 0.66666666666667 0.60024677538198 127.600 > 65 76 Te 0.83333333333333 0.66666666666667 0.60024677538198 127.600 > 65 77 Te 0.08333333333333 0.91666666666667 0.60024677538198 127.600 > 65 78 Te 0.33333333333333 0.91666666666667 0.60024677538198 127.600 > 65 79 Te 0.58333333333333 0.91666666666667 0.60024677538198 127.600 > 65 80 Te 0.83333333333333 0.91666666666667 0.60024677538198 127.600 > 65 81 Te 0.16666666666667 0.08333333333333 0.10024677538198 127.600 > 81 82 Te 0.41666666666667 0.08333333333333 0.10024677538198 127.600 > 81 83 Te 0.66666666666667 0.08333333333333 0.10024677538198 127.600 > 81 84 Te 0.91666666666667 0.08333333333333 0.10024677538198 127.600 > 81 85 Te 0.16666666666667 0.33333333333333 0.10024677538198 127.600 > 81 86 Te 0.41666666666667 0.33333333333333 0.10024677538198 127.600 > 81 87 Te 0.66666666666667 0.33333333333333 0.10024677538198 127.600 > 81 88 Te 0.91666666666667 0.33333333333333 0.10024677538198 127.600 > 81 89 Te 0.16666666666667 0.58333333333333 0.10024677538198 127.600 > 81 90 Te 0.41666666666667 0.58333333333333 0.10024677538198 127.600 > 81 91 Te 0.66666666666667 0.58333333333333 0.10024677538198 127.600 > 81 92 Te 0.91666666666667 0.58333333333333 0.10024677538198 127.600 > 81 93 Te 0.16666666666667 0.83333333333333 0.10024677538198 127.600 > 81 94 Te 0.41666666666667 0.83333333333333 0.10024677538198 127.600 > 81 95 Te 0.66666666666667 0.83333333333333 0.10024677538198 127.600 > 81 96 Te 0.91666666666667 0.83333333333333 0.10024677538198 127.600 > 81 97 Te 0.16666666666667 0.08333333333333 0.39975322461802 127.600 > 97 98 Te 0.41666666666667 0.08333333333333 0.39975322461802 127.600 > 97 99 Te 0.66666666666667 0.08333333333333 0.39975322461802 127.600 > 97 100 Te 0.91666666666667 0.08333333333333 0.39975322461802 127.600 > 97 101 Te 0.16666666666667 0.33333333333333 0.39975322461802 127.600 > 97 102 Te 0.41666666666667 0.33333333333333 0.39975322461802 127.600 > 97 103 Te 0.66666666666667 0.33333333333333 0.39975322461802 127.600 > 97 104 Te 0.91666666666667 0.33333333333333 0.39975322461802 127.600 > 97 105 Te 0.16666666666667 0.58333333333333 0.39975322461802 127.600 > 97 106 Te 0.41666666666667 0.58333333333333 0.39975322461802 127.600 > 97 107 Te 0.66666666666667 0.58333333333333 0.39975322461802 127.600 > 97 108 Te 0.91666666666667 0.58333333333333 0.39975322461802 127.600 > 97 109 Te 0.16666666666667 0.83333333333333 0.39975322461802 127.600 > 97 110 Te 0.41666666666667 0.83333333333333 0.39975322461802 127.600 > 97 111 Te 0.66666666666667 0.83333333333333 0.39975322461802 127.600 > 97 112 Te 0.91666666666667 0.83333333333333 0.39975322461802 127.600 > 97 113 Te 0.08333333333333 0.16666666666667 0.89975322461802 127.600 > 113 114 Te 0.33333333333333 0.16666666666667 0.89975322461802 127.600 > 113 115 Te 0.58333333333333 0.16666666666667 0.89975322461802 127.600 > 113 116 Te 0.83333333333333 0.16666666666667 0.89975322461802 127.600 > 113 117 Te 0.08333333333333 0.41666666666667 0.89975322461802 127.600 > 113 118 Te 0.33333333333333 0.41666666666667 0.89975322461802 127.600 > 113 119 Te 0.58333333333333 0.41666666666667 0.89975322461802 127.600 > 113 120 Te 0.83333333333333 0.41666666666667 0.89975322461802 127.600 > 113 121 Te 0.08333333333333 0.66666666666667 0.89975322461802 127.600 > 113 122 Te 0.33333333333333 0.66666666666667 0.89975322461802 127.600 > 113 123 Te 0.58333333333333 0.66666666666667 0.89975322461802 127.600 > 113 124 Te 0.83333333333333 0.66666666666667 0.89975322461802 127.600 > 113 125 Te 0.08333333333333 0.91666666666667 0.89975322461802 127.600 > 113 126 Te 0.33333333333333 0.91666666666667 0.89975322461802 127.600 > 113 127 Te 0.58333333333333 0.91666666666667 0.89975322461802 127.600 > 113 128 Te 0.83333333333333 0.91666666666667 0.89975322461802 127.600 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 8.5031623 0.0000000 0.0000000 -0.0000000 8.5031623 0.0000000 0.0000000 0.0000000 6.0008867 -------------------------- Born effective charges -------------------------- 1 Rb 1.1052389 0.0000000 0.0000000 0.0000000 1.1052389 0.0000000 0.0000000 0.0000000 2.1323423 2 Rb 1.1052389 0.0000000 0.0000000 0.0000000 1.1052389 0.0000000 0.0000000 0.0000000 2.1323423 3 Y 4.8553139 0.0000000 0.0000000 0.0000000 4.8553139 0.0000000 0.0000000 0.0000000 1.8077586 4 Y 4.8553139 0.0000000 0.0000000 0.0000000 4.8553139 0.0000000 0.0000000 0.0000000 1.8077586 5 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 6 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 7 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 8 Te -2.9802764 -0.0000000 0.0000000 -0.0000000 -2.9802764 0.0000000 0.0000000 0.0000000 -1.9700505 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000004 (yyy) 0.00000004 (yyy) 0.00000004 (yyy) Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 13 13 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.63, Number of G-points: 291, Lambda: 0.18 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/42) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 42 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.883 ( 0.000 0.000 0.000) 0.000 0.883 ( 0.000 0.000 0.000) 0.000 1.506 ( 0.000 0.000 0.000) 0.000 1.616 ( 0.000 0.000 0.000) 0.000 1.616 ( 0.000 0.000 0.000) 0.000 1.992 ( 0.000 0.000 0.000) 0.000 1.992 ( 0.000 0.000 0.000) 0.000 2.457 ( 0.000 0.000 0.000) 0.000 2.755 ( 0.000 0.000 0.000) 0.000 2.755 ( 0.000 0.000 0.000) 0.000 2.816 ( 0.000 0.000 0.000) 0.000 2.880 ( 0.000 0.000 0.000) 0.000 2.880 ( 0.000 0.000 0.000) 0.000 3.510 ( 0.000 0.000 0.000) 0.000 3.596 ( 0.000 0.000 0.000) 0.000 3.596 ( 0.000 0.000 0.000) 0.000 3.615 ( 0.000 0.000 0.000) 0.000 3.615 ( 0.000 0.000 0.000) 0.000 3.962 ( 0.000 0.000 0.000) 0.000 5.342 ( 0.000 0.000 0.000) 0.000 5.683 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/42) ======================= q-point: ( 0.08 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.379 ( 16.458 9.502 0.000) 19.005 0.382 ( 15.630 9.024 0.000) 18.048 0.662 ( 27.865 16.088 0.000) 32.176 0.968 ( 6.737 3.890 0.000) 7.779 1.013 ( 9.694 5.597 0.000) 11.194 1.301 ( -10.043 -5.798 0.000) 11.596 1.621 ( 0.494 0.285 0.000) 0.571 1.823 ( 11.238 6.488 0.000) 12.977 2.001 ( 0.811 0.468 0.000) 0.936 2.028 ( 1.738 1.004 0.000) 2.007 2.407 ( -1.472 -0.850 0.000) 1.699 2.645 ( -5.803 -3.350 0.000) 6.701 2.743 ( -0.920 -0.531 0.000) 1.063 2.870 ( -0.743 -0.429 0.000) 0.858 2.924 ( 5.700 3.291 0.000) 6.582 2.957 ( 5.219 3.013 0.000) 6.026 3.501 ( -0.729 -0.421 0.000) 0.842 3.606 ( 0.768 0.443 0.000) 0.886 3.625 ( 0.824 0.476 0.000) 0.952 3.690 ( 0.874 0.504 0.000) 1.009 4.157 ( 17.685 10.210 0.000) 20.421 5.247 ( -7.884 -4.552 0.000) 9.104 5.287 ( -1.686 -0.974 0.000) 1.947 5.545 ( -10.348 -5.974 0.000) 11.949 ======================= Grid point 2 (3/42) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.704 ( 12.077 6.972 0.000) 13.945 0.755 ( 15.836 9.143 0.000) 18.286 1.113 ( -3.345 -1.931 0.000) 3.863 1.154 ( 8.611 4.972 0.000) 9.944 1.257 ( 10.595 6.117 0.000) 12.234 1.285 ( 23.834 13.761 0.000) 27.522 1.638 ( 0.902 0.521 0.000) 1.042 2.029 ( 1.549 0.895 0.000) 1.789 2.042 ( 2.894 1.671 0.000) 3.342 2.122 ( 10.438 6.026 0.000) 12.052 2.429 ( 0.087 0.050 0.000) 0.100 2.531 ( -0.968 -0.559 0.000) 1.117 2.720 ( -0.922 -0.532 0.000) 1.065 2.855 ( -0.446 -0.258 0.000) 0.515 3.051 ( 4.974 2.872 0.000) 5.743 3.059 ( 2.960 1.709 0.000) 3.417 3.482 ( -0.807 -0.466 0.000) 0.932 3.617 ( -5.666 -3.271 0.000) 6.543 3.625 ( 0.796 0.459 0.000) 0.919 3.646 ( 0.808 0.467 0.000) 0.933 4.552 ( 11.612 6.704 0.000) 13.408 5.019 ( -9.642 -5.567 0.000) 11.133 5.094 ( -15.982 -9.227 0.000) 18.454 5.385 ( -1.690 -0.976 0.000) 1.952 ======================= Grid point 3 (4/42) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.938 ( 8.369 4.832 0.000) 9.664 1.090 ( 12.707 7.336 0.000) 14.672 1.103 ( 3.225 1.862 0.000) 3.724 1.342 ( 7.475 4.316 0.000) 8.631 1.500 ( 11.018 6.361 0.000) 12.722 1.660 ( 0.963 0.556 0.000) 1.112 1.817 ( 21.448 12.383 0.000) 24.765 2.046 ( -1.346 -0.777 0.000) 1.554 2.073 ( 2.231 1.288 0.000) 2.576 2.310 ( -6.959 -4.018 0.000) 8.036 2.380 ( 10.185 5.880 0.000) 11.761 2.666 ( 8.193 4.730 0.000) 9.460 2.705 ( -0.331 -0.191 0.000) 0.382 2.853 ( 0.241 0.139 0.000) 0.278 3.095 ( 0.550 0.317 0.000) 0.635 3.135 ( 1.866 1.078 0.000) 2.155 3.467 ( -0.458 -0.265 0.000) 0.529 3.471 ( -6.047 -3.491 0.000) 6.982 3.637 ( 0.128 0.074 0.000) 0.148 3.657 ( 0.105 0.061 0.000) 0.121 4.504 ( -12.603 -7.276 0.000) 14.553 4.677 ( -18.683 -10.787 0.000) 21.573 4.989 ( 5.056 2.919 0.000) 5.838 5.373 ( -0.249 -0.143 0.000) 0.287 ======================= Grid point 4 (5/42) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.102 ( 6.091 3.517 0.000) 7.034 1.250 ( 9.042 5.220 0.000) 10.441 1.327 ( 7.758 4.479 0.000) 8.958 1.497 ( 5.938 3.428 0.000) 6.856 1.680 ( 0.814 0.470 0.000) 0.939 1.773 ( 12.099 6.986 0.000) 13.971 1.982 ( -3.192 -1.843 0.000) 3.686 2.128 ( 2.418 1.396 0.000) 2.792 2.166 ( -4.530 -2.615 0.000) 5.231 2.280 ( 17.815 10.286 0.000) 20.571 2.553 ( 4.143 2.392 0.000) 4.784 2.702 ( 0.006 0.004 0.000) 0.007 2.797 ( 1.183 0.683 0.000) 1.366 2.861 ( 0.357 0.206 0.000) 0.412 3.102 ( 0.352 0.203 0.000) 0.406 3.128 ( -1.913 -1.104 0.000) 2.209 3.385 ( -0.746 -0.431 0.000) 0.862 3.459 ( -0.375 -0.217 0.000) 0.434 3.630 ( -0.742 -0.429 0.000) 0.857 3.650 ( -0.736 -0.425 0.000) 0.849 4.198 ( -12.007 -6.932 0.000) 13.865 4.275 ( -15.291 -8.828 0.000) 17.657 5.089 ( 2.930 1.691 0.000) 3.383 5.357 ( -1.241 -0.717 0.000) 1.433 ======================= Grid point 5 (6/42) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.224 ( 4.370 2.523 0.000) 5.046 1.439 ( 1.780 1.028 0.000) 2.055 1.488 ( 10.645 6.146 0.000) 12.292 1.617 ( 4.383 2.531 0.000) 5.061 1.699 ( 0.792 0.457 0.000) 0.915 1.916 ( -1.731 -1.000 0.000) 1.999 1.994 ( 6.002 3.465 0.000) 6.930 2.178 ( 1.808 1.044 0.000) 2.087 2.184 ( 7.278 4.202 0.000) 8.404 2.512 ( -3.949 -2.280 0.000) 4.560 2.673 ( 12.246 7.070 0.000) 14.140 2.702 ( -0.091 -0.053 0.000) 0.106 2.707 ( -7.216 -4.166 0.000) 8.332 2.865 ( -0.040 -0.023 0.000) 0.047 3.051 ( -5.322 -3.073 0.000) 6.145 3.116 ( 0.709 0.409 0.000) 0.819 3.435 ( 4.180 2.413 0.000) 4.827 3.447 ( -0.671 -0.387 0.000) 0.774 3.606 ( -1.138 -0.657 0.000) 1.314 3.627 ( -1.076 -0.621 0.000) 1.243 3.986 ( -6.317 -3.647 0.000) 7.294 4.003 ( -8.054 -4.650 0.000) 9.300 5.128 ( 0.706 0.408 0.000) 0.815 5.317 ( -2.016 -1.164 0.000) 2.328 ======================= Grid point 6 (7/42) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.295 ( 1.658 0.957 0.000) 1.914 1.436 ( -0.787 -0.454 0.000) 0.908 1.665 ( 3.262 1.884 0.000) 3.767 1.690 ( 1.749 1.010 0.000) 2.019 1.715 ( 0.444 0.257 0.000) 0.513 1.956 ( 5.215 3.011 0.000) 6.022 2.122 ( 4.936 2.850 0.000) 5.700 2.207 ( 0.659 0.381 0.000) 0.761 2.336 ( -8.237 -4.755 0.000) 9.511 2.357 ( 3.731 2.154 0.000) 4.308 2.560 ( -2.989 -1.726 0.000) 3.452 2.699 ( -0.110 -0.063 0.000) 0.127 2.861 ( -0.155 -0.090 0.000) 0.179 2.904 ( 6.223 3.593 0.000) 7.186 2.908 ( -5.383 -3.108 0.000) 6.216 3.119 ( -0.534 -0.308 0.000) 0.616 3.432 ( -0.414 -0.239 0.000) 0.478 3.522 ( 2.445 1.412 0.000) 2.823 3.585 ( -0.561 -0.324 0.000) 0.648 3.607 ( -0.515 -0.297 0.000) 0.595 3.893 ( -1.991 -1.149 0.000) 2.299 3.895 ( -1.883 -1.087 0.000) 2.174 5.133 ( -0.008 -0.004 0.000) 0.009 5.278 ( -1.077 -0.622 0.000) 1.243 ======================= Grid point 14 (8/42) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.654 ( 9.309 16.124 0.000) 18.618 0.668 ( 9.325 16.152 0.000) 18.651 1.086 ( 13.352 23.126 0.000) 26.704 1.120 ( 5.812 10.067 0.000) 11.624 1.149 ( -4.146 -7.180 0.000) 8.291 1.176 ( 5.589 9.680 0.000) 11.177 1.658 ( 1.908 3.304 0.000) 3.815 1.972 ( 3.124 5.411 0.000) 6.248 2.042 ( 0.301 0.521 0.000) 0.601 2.071 ( 4.688 8.121 0.000) 9.377 2.422 ( 1.780 3.082 0.000) 3.559 2.541 ( -3.556 -6.160 0.000) 7.113 2.737 ( -0.096 -0.166 0.000) 0.192 2.883 ( 0.121 0.209 0.000) 0.242 3.004 ( 1.272 2.202 0.000) 2.543 3.006 ( 2.463 4.266 0.000) 4.925 3.487 ( -0.506 -0.877 0.000) 1.013 3.632 ( 1.154 1.999 0.000) 2.308 3.649 ( -2.707 -4.689 0.000) 5.415 3.650 ( 1.021 1.769 0.000) 2.042 4.470 ( 9.866 17.089 0.000) 19.733 5.076 ( -6.917 -11.981 0.000) 13.834 5.276 ( -0.863 -1.494 0.000) 1.726 5.311 ( -8.645 -14.973 0.000) 17.289 ======================= Grid point 15 (9/42) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.957 ( 4.514 16.201 0.000) 16.818 0.982 ( 10.974 12.997 0.000) 17.010 1.064 ( 0.532 -1.829 0.000) 1.905 1.315 ( 4.227 10.043 0.000) 10.897 1.387 ( 9.357 6.738 0.000) 11.531 1.508 ( 14.869 7.078 0.000) 16.468 1.773 ( 3.919 13.318 0.000) 13.883 2.035 ( -0.327 -1.621 0.000) 1.653 2.049 ( 3.205 1.499 0.000) 3.538 2.286 ( 7.940 8.584 0.000) 11.693 2.396 ( -6.417 -2.843 0.000) 7.019 2.536 ( 6.659 0.293 0.000) 6.666 2.749 ( 0.134 3.812 0.000) 3.815 2.875 ( -1.458 1.552 0.000) 2.130 3.038 ( 2.746 -2.774 0.000) 3.903 3.079 ( 3.151 -0.088 0.000) 3.152 3.470 ( -0.465 -0.676 0.000) 0.820 3.522 ( -5.070 -5.970 0.000) 7.832 3.675 ( 0.004 4.178 0.000) 4.178 3.686 ( 0.082 3.265 0.000) 3.266 4.667 ( -8.357 2.724 0.000) 8.790 4.873 ( 4.667 -7.524 0.000) 8.854 4.928 ( -19.068 -8.253 0.000) 20.777 5.284 ( 4.640 -7.506 0.000) 8.825 ======================= Grid point 16 (10/42) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.114 ( 5.638 3.482 0.000) 6.626 1.170 ( 1.791 13.292 0.000) 13.412 1.266 ( 7.886 10.215 0.000) 12.905 1.467 ( 3.432 7.311 0.000) 8.076 1.631 ( 11.384 6.283 0.000) 13.003 1.693 ( 3.371 3.675 0.000) 4.987 1.988 ( -0.242 -4.147 0.000) 4.154 2.058 ( 7.207 7.579 0.000) 10.458 2.164 ( 8.243 9.887 0.000) 12.873 2.259 ( -6.892 -0.208 0.000) 6.895 2.473 ( 6.491 2.460 0.000) 6.941 2.623 ( 4.986 -3.823 0.000) 6.282 2.826 ( 1.084 6.484 0.000) 6.574 2.873 ( -0.587 1.558 0.000) 1.664 3.050 ( 2.556 -3.356 0.000) 4.218 3.097 ( 1.267 -4.159 0.000) 4.347 3.388 ( -3.102 -4.538 0.000) 5.497 3.458 ( -0.230 -0.722 0.000) 0.758 3.711 ( -2.088 6.370 0.000) 6.703 3.715 ( -1.643 5.053 0.000) 5.313 4.380 ( -15.769 -5.096 0.000) 16.572 4.513 ( -19.868 -5.840 0.000) 20.708 4.982 ( 7.972 -3.341 0.000) 8.644 5.284 ( 4.146 -8.043 0.000) 9.049 ======================= Grid point 17 (11/42) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.274 ( 3.889 8.826 0.000) 9.645 1.328 ( 5.759 7.106 0.000) 9.147 1.447 ( 2.293 7.461 0.000) 7.806 1.578 ( 3.257 3.902 0.000) 5.083 1.765 ( -0.487 7.071 0.000) 7.088 1.894 ( 11.317 5.268 0.000) 12.483 1.929 ( -1.585 -2.690 0.000) 3.122 2.155 ( 2.234 1.469 0.000) 2.674 2.190 ( -2.672 5.291 0.000) 5.928 2.470 ( 8.223 4.476 0.000) 9.362 2.563 ( 6.427 2.265 0.000) 6.815 2.643 ( 2.020 -4.511 0.000) 4.943 2.827 ( -5.558 1.042 0.000) 5.655 2.878 ( -0.692 1.036 0.000) 1.246 3.055 ( -0.771 -5.443 0.000) 5.497 3.069 ( 2.740 -2.075 0.000) 3.437 3.363 ( 3.445 -0.878 0.000) 3.555 3.445 ( -0.315 -1.120 0.000) 1.163 3.728 ( -4.423 8.427 0.000) 9.517 3.731 ( -3.501 7.239 0.000) 8.041 4.095 ( -11.437 -3.427 0.000) 11.940 4.156 ( -15.284 -3.242 0.000) 15.624 5.054 ( 5.233 -4.800 0.000) 7.101 5.264 ( 2.850 -8.074 0.000) 8.562 ======================= Grid point 18 (12/42) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.385 ( -3.029 12.572 0.000) 12.931 1.487 ( -2.934 3.758 0.000) 4.767 1.497 ( 8.303 -1.727 0.000) 8.480 1.670 ( 2.862 3.030 0.000) 4.168 1.846 ( 1.002 9.685 0.000) 9.736 1.900 ( -0.169 0.850 0.000) 0.867 2.102 ( 5.622 7.653 0.000) 9.497 2.175 ( 0.039 -2.693 0.000) 2.693 2.284 ( 5.660 4.410 0.000) 7.175 2.432 ( -4.819 -3.113 0.000) 5.736 2.573 ( -3.841 -6.294 0.000) 7.373 2.752 ( -2.859 2.626 0.000) 3.882 2.809 ( 9.586 5.753 0.000) 11.179 2.870 ( -0.769 0.216 0.000) 0.798 2.947 ( -4.436 -6.048 0.000) 7.501 3.091 ( 1.806 -2.136 0.000) 2.797 3.424 ( -0.486 -1.408 0.000) 1.489 3.443 ( 4.867 -0.483 0.000) 4.891 3.728 ( -6.306 10.338 0.000) 12.109 3.740 ( -5.455 9.955 0.000) 11.352 3.928 ( -4.650 -2.240 0.000) 5.161 3.940 ( -6.722 -1.672 0.000) 6.927 5.077 ( 3.509 -5.147 0.000) 6.229 5.228 ( 2.260 -7.344 0.000) 7.684 ======================= Grid point 19 (13/42) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.423 ( -6.381 11.052 0.000) 12.762 1.473 ( -2.340 4.053 0.000) 4.680 1.553 ( 3.104 -5.376 0.000) 6.208 1.714 ( -0.801 1.388 0.000) 1.602 1.890 ( -5.336 9.242 0.000) 10.672 1.991 ( 0.061 -0.106 0.000) 0.122 2.103 ( 4.016 -6.956 0.000) 8.032 2.223 ( -3.756 6.506 0.000) 7.513 2.353 ( 1.207 -2.090 0.000) 2.414 2.363 ( -2.755 4.772 0.000) 5.510 2.487 ( 2.401 -4.159 0.000) 4.802 2.765 ( -2.822 4.888 0.000) 5.644 2.842 ( 1.553 -2.690 0.000) 3.106 2.863 ( -0.079 0.137 0.000) 0.158 2.934 ( 0.473 -0.819 0.000) 0.946 3.089 ( 1.263 -2.187 0.000) 2.525 3.410 ( 0.509 -0.882 0.000) 1.018 3.486 ( 1.541 -2.670 0.000) 3.083 3.724 ( -6.355 11.008 0.000) 12.711 3.742 ( -6.784 11.751 0.000) 13.569 3.878 ( 0.267 -0.462 0.000) 0.534 3.879 ( 0.305 -0.528 0.000) 0.610 5.080 ( 2.954 -5.117 0.000) 5.908 5.209 ( 3.526 -6.108 0.000) 7.053 ======================= Grid point 29 (14/42) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.061 ( 1.445 2.503 0.000) 2.890 1.251 ( 7.624 13.206 0.000) 15.249 1.287 ( 8.609 14.912 0.000) 17.219 1.518 ( 4.883 8.457 0.000) 9.766 1.520 ( 4.484 7.767 0.000) 8.969 1.632 ( 3.853 6.673 0.000) 7.705 1.963 ( -1.449 -2.510 0.000) 2.898 2.095 ( 1.816 3.145 0.000) 3.631 2.102 ( 8.648 14.978 0.000) 17.295 2.346 ( -1.023 -1.771 0.000) 2.045 2.424 ( 2.632 4.559 0.000) 5.265 2.530 ( -0.210 -0.363 0.000) 0.419 2.855 ( 3.326 5.762 0.000) 6.653 2.895 ( 0.181 0.314 0.000) 0.362 2.998 ( -0.263 -0.456 0.000) 0.527 3.055 ( -1.337 -2.315 0.000) 2.673 3.395 ( -3.540 -6.131 0.000) 7.079 3.456 ( -0.453 -0.785 0.000) 0.907 3.758 ( 2.023 3.504 0.000) 4.046 3.766 ( 2.491 4.314 0.000) 4.981 4.555 ( -6.777 -11.738 0.000) 13.554 4.720 ( -7.528 -13.039 0.000) 15.056 4.834 ( 1.221 2.115 0.000) 2.443 5.150 ( -2.813 -4.872 0.000) 5.626 ======================= Grid point 30 (15/42) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.166 ( 6.642 4.038 0.000) 7.773 1.485 ( 5.006 10.990 0.000) 12.076 1.518 ( -1.374 18.737 0.000) 18.787 1.601 ( -0.735 5.402 0.000) 5.452 1.742 ( 11.133 5.468 0.000) 12.403 1.788 ( 6.051 5.114 0.000) 7.922 1.898 ( 0.025 -4.472 0.000) 4.472 2.150 ( 0.821 3.378 0.000) 3.476 2.327 ( -4.176 6.846 0.000) 8.019 2.407 ( 7.142 10.705 0.000) 12.869 2.471 ( 4.249 -0.682 0.000) 4.304 2.540 ( 4.050 -4.382 0.000) 5.967 2.897 ( -0.355 0.618 0.000) 0.713 2.917 ( -3.129 1.975 0.000) 3.700 2.999 ( 0.519 -5.240 0.000) 5.266 3.015 ( 1.810 0.130 0.000) 1.815 3.316 ( 0.254 -1.861 0.000) 1.878 3.439 ( -0.396 -1.118 0.000) 1.186 3.818 ( 0.276 4.347 0.000) 4.356 3.835 ( 0.103 4.891 0.000) 4.892 4.306 ( -12.792 -2.133 0.000) 12.969 4.429 ( -15.577 -2.896 0.000) 15.844 4.881 ( 6.025 -5.965 0.000) 8.478 5.089 ( 3.933 -10.421 0.000) 11.139 ======================= Grid point 31 (16/42) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.345 ( 8.351 3.733 0.000) 9.147 1.601 ( -3.714 6.328 0.000) 7.337 1.638 ( -5.424 18.000 0.000) 18.799 1.656 ( 3.531 5.601 0.000) 6.622 1.817 ( 0.727 -2.717 0.000) 2.812 1.990 ( 3.202 8.953 0.000) 9.508 2.000 ( 11.901 5.069 0.000) 12.935 2.160 ( -1.573 -2.095 0.000) 2.620 2.336 ( -3.121 7.346 0.000) 7.982 2.428 ( 0.884 -2.613 0.000) 2.758 2.558 ( 0.862 -2.645 0.000) 2.782 2.707 ( 9.462 7.649 0.000) 12.167 2.822 ( -4.555 -2.127 0.000) 5.027 2.889 ( -1.121 -0.015 0.000) 1.121 2.951 ( -1.799 -3.665 0.000) 4.083 3.043 ( 1.450 -0.854 0.000) 1.683 3.359 ( 4.185 0.176 0.000) 4.188 3.420 ( -0.412 -1.087 0.000) 1.163 3.852 ( -0.485 2.976 0.000) 3.016 3.867 ( -0.832 2.901 0.000) 3.018 4.097 ( -12.581 5.279 0.000) 13.644 4.168 ( -14.987 5.135 0.000) 15.843 4.918 ( 6.388 -8.167 0.000) 10.369 5.062 ( 5.306 -11.397 0.000) 12.571 ======================= Grid point 32 (17/42) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.477 ( 3.160 1.024 0.000) 3.322 1.587 ( -4.280 6.201 0.000) 7.535 1.696 ( -7.871 16.520 0.000) 18.300 1.746 ( 0.374 4.412 0.000) 4.428 1.870 ( 7.956 -0.832 0.000) 8.000 2.038 ( -5.022 -7.007 0.000) 8.621 2.101 ( -2.659 9.284 0.000) 9.657 2.228 ( 9.231 1.164 0.000) 9.304 2.362 ( -0.505 4.902 0.000) 4.928 2.439 ( -0.496 2.865 0.000) 2.907 2.481 ( -6.537 -1.845 0.000) 6.792 2.781 ( 1.295 -1.012 0.000) 1.643 2.864 ( -0.595 -1.131 0.000) 1.278 2.878 ( -4.124 0.871 0.000) 4.215 2.887 ( 4.889 2.018 0.000) 5.290 3.044 ( 0.860 -2.593 0.000) 2.732 3.406 ( -0.359 -0.148 0.000) 0.388 3.421 ( 3.184 -1.925 0.000) 3.721 3.863 ( 0.276 -0.311 0.000) 0.416 3.866 ( -0.122 -0.864 0.000) 0.872 3.993 ( -9.703 12.319 0.000) 15.681 4.032 ( -10.865 12.471 0.000) 16.541 4.937 ( 5.422 -8.469 0.000) 10.056 5.044 ( 5.457 -10.530 0.000) 11.861 ======================= Grid point 43 (18/42) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.289 ( 4.613 7.990 0.000) 9.226 1.645 ( 2.110 3.654 0.000) 4.220 1.708 ( 3.814 6.607 0.000) 7.629 1.796 ( -2.082 -3.606 0.000) 4.164 1.875 ( 8.252 14.292 0.000) 16.503 1.884 ( 5.105 8.842 0.000) 10.209 1.910 ( 4.103 7.106 0.000) 8.206 2.177 ( -1.180 -2.044 0.000) 2.360 2.402 ( -0.582 -1.009 0.000) 1.165 2.447 ( -1.858 -3.218 0.000) 3.716 2.500 ( 4.428 7.669 0.000) 8.856 2.656 ( 6.030 10.444 0.000) 12.060 2.855 ( -3.795 -6.572 0.000) 7.589 2.900 ( -0.273 -0.472 0.000) 0.545 2.945 ( -0.521 -0.902 0.000) 1.041 3.027 ( 0.201 0.347 0.000) 0.401 3.320 ( 0.972 1.684 0.000) 1.945 3.418 ( -0.479 -0.830 0.000) 0.958 3.867 ( 0.327 0.566 0.000) 0.653 3.882 ( 0.007 0.012 0.000) 0.014 4.285 ( -0.380 -0.658 0.000) 0.760 4.379 ( -1.519 -2.631 0.000) 3.038 4.780 ( -1.958 -3.391 0.000) 3.915 4.912 ( -3.707 -6.421 0.000) 7.414 ======================= Grid point 44 (19/42) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.456 ( 4.460 7.327 0.000) 8.578 1.698 ( -0.302 5.000 0.000) 5.009 1.754 ( 2.262 -2.826 0.000) 3.619 1.808 ( -0.869 7.619 0.000) 7.668 1.976 ( -4.974 13.682 0.000) 14.558 2.044 ( -4.585 -7.453 0.000) 8.750 2.066 ( 8.669 0.774 0.000) 8.703 2.186 ( 5.982 9.643 0.000) 11.348 2.389 ( -0.928 -0.535 0.000) 1.071 2.438 ( 2.674 3.091 0.000) 4.087 2.597 ( -3.397 2.702 0.000) 4.340 2.759 ( 0.152 -0.682 0.000) 0.698 2.827 ( 4.953 4.379 0.000) 6.611 2.866 ( -1.904 -2.833 0.000) 3.413 2.930 ( -0.925 0.527 0.000) 1.064 3.017 ( -0.425 -1.964 0.000) 2.010 3.356 ( 2.010 -0.745 0.000) 2.143 3.411 ( 0.005 0.222 0.000) 0.222 3.855 ( -0.563 -1.511 0.000) 1.613 3.858 ( -0.834 -1.995 0.000) 2.163 4.265 ( -6.720 9.308 0.000) 11.480 4.324 ( -7.902 8.777 0.000) 11.809 4.750 ( 4.625 -8.130 0.000) 9.353 4.838 ( 4.226 -10.350 0.000) 11.180 ======================= Grid point 45 (20/42) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.534 ( -2.601 4.504 0.000) 5.201 1.724 ( -3.908 6.769 0.000) 7.817 1.837 ( -2.690 4.660 0.000) 5.381 1.842 ( 3.062 -5.304 0.000) 6.124 1.868 ( 2.497 -4.324 0.000) 4.993 2.009 ( -7.819 13.543 0.000) 15.638 2.146 ( 4.012 -6.950 0.000) 8.025 2.290 ( -4.242 7.347 0.000) 8.484 2.415 ( -0.149 0.258 0.000) 0.298 2.504 ( -1.735 3.005 0.000) 3.470 2.517 ( -1.449 2.510 0.000) 2.898 2.786 ( -1.856 3.214 0.000) 3.711 2.814 ( 2.136 -3.700 0.000) 4.272 2.902 ( 0.044 -0.076 0.000) 0.088 2.924 ( -1.110 1.922 0.000) 2.220 2.991 ( 1.344 -2.328 0.000) 2.688 3.369 ( 1.651 -2.859 0.000) 3.302 3.414 ( -0.364 0.631 0.000) 0.729 3.837 ( 0.777 -1.345 0.000) 1.554 3.839 ( 0.724 -1.254 0.000) 1.448 4.251 ( -7.184 12.443 0.000) 14.368 4.297 ( -7.391 12.802 0.000) 14.782 4.753 ( 5.513 -9.549 0.000) 11.026 4.827 ( 6.246 -10.818 0.000) 12.491 ======================= Grid point 58 (21/42) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.619 ( 4.428 7.670 0.000) 8.857 1.721 ( 0.277 0.480 0.000) 0.554 1.802 ( 2.523 4.370 0.000) 5.047 1.852 ( -5.642 -9.773 0.000) 11.285 1.921 ( 1.881 3.258 0.000) 3.762 2.030 ( -1.808 -3.131 0.000) 3.616 2.166 ( 2.412 4.178 0.000) 4.825 2.362 ( 3.172 5.495 0.000) 6.345 2.415 ( 1.991 3.449 0.000) 3.982 2.534 ( -3.266 -5.657 0.000) 6.532 2.542 ( 4.410 7.638 0.000) 8.819 2.759 ( -4.290 -7.430 0.000) 8.580 2.850 ( 3.158 5.470 0.000) 6.316 2.897 ( 1.696 2.937 0.000) 3.391 2.934 ( -0.389 -0.674 0.000) 0.778 2.967 ( -1.609 -2.787 0.000) 3.218 3.334 ( -0.591 -1.024 0.000) 1.183 3.422 ( 0.302 0.524 0.000) 0.605 3.822 ( -0.673 -1.166 0.000) 1.346 3.825 ( -0.572 -0.991 0.000) 1.145 4.427 ( 3.627 6.283 0.000) 7.255 4.479 ( 3.519 6.095 0.000) 7.038 4.608 ( -3.292 -5.702 0.000) 6.584 4.666 ( -3.710 -6.427 0.000) 7.421 ======================= Grid point 181 (22/42) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 42 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.384 ( 0.000 -0.000 18.771) 18.771 0.384 ( -0.000 0.000 18.771) 18.771 0.632 ( 0.000 0.000 31.167) 31.167 0.720 ( 0.000 0.000 -14.348) 14.348 0.720 ( 0.000 -0.000 -14.348) 14.348 1.203 ( -0.000 0.000 -25.355) 25.355 1.730 ( 0.000 -0.000 9.379) 9.379 1.730 ( 0.000 0.000 9.379) 9.379 1.910 ( 0.000 0.000 -7.588) 7.588 1.910 ( 0.000 -0.000 -7.588) 7.588 2.701 ( 0.000 0.000 19.425) 19.425 2.791 ( 0.000 0.000 3.129) 3.129 2.791 ( 0.000 0.000 3.129) 3.129 2.853 ( 0.000 0.000 -2.607) 2.607 2.853 ( -0.000 0.000 -2.607) 2.607 3.098 ( 0.000 0.000 -19.532) 19.532 3.601 ( 0.000 0.000 0.446) 0.446 3.601 ( -0.000 0.000 0.446) 0.446 3.610 ( 0.000 -0.000 -0.419) 0.419 3.610 ( 0.000 0.000 -0.419) 0.419 3.736 ( 0.000 0.000 11.121) 11.121 3.904 ( 0.000 0.000 -5.823) 5.823 5.688 ( 0.000 -0.000 0.497) 0.497 5.699 ( 0.000 -0.000 -0.488) 0.488 ======================= Grid point 182 (23/42) ======================= q-point: ( 0.08 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.519 ( 9.656 5.575 13.057) 17.170 0.541 ( 11.020 6.362 13.309) 18.413 0.814 ( 5.807 3.353 -14.958) 16.392 0.817 ( 7.490 4.324 -12.900) 15.531 0.899 ( 19.416 11.210 18.974) 29.371 1.250 ( 5.188 2.995 -11.972) 13.388 1.739 ( 0.715 0.413 9.489) 9.525 1.792 ( 4.778 2.759 3.654) 6.618 1.920 ( 0.790 0.456 -7.608) 7.663 1.938 ( 1.313 0.758 -7.915) 8.058 2.556 ( -4.150 -2.396 9.789) 10.899 2.661 ( -6.076 -3.508 -0.399) 7.027 2.780 ( -0.857 -0.495 3.179) 3.329 2.842 ( -0.774 -0.447 -2.636) 2.783 2.916 ( 3.461 1.998 -0.413) 4.018 2.937 ( 3.702 2.137 -1.951) 4.699 3.578 ( -1.874 -1.082 6.061) 6.436 3.611 ( 0.782 0.452 0.466) 1.016 3.621 ( 0.811 0.468 -0.439) 1.034 3.671 ( -1.147 -0.662 -2.562) 2.884 4.251 ( 23.990 13.851 10.814) 29.737 4.645 ( 26.345 15.210 -29.675) 42.497 5.496 ( -13.400 -7.737 4.020) 15.987 5.536 ( -11.489 -6.633 -0.931) 13.298 ======================= Grid point 183 (24/42) ======================= q-point: ( 0.15 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.786 ( 12.854 7.421 4.183) 15.421 0.802 ( 10.668 6.159 9.086) 15.307 0.935 ( 5.610 3.239 -10.440) 12.286 1.015 ( 8.839 5.103 -11.066) 15.054 1.367 ( 16.695 9.639 9.366) 21.433 1.443 ( 12.191 7.038 5.426) 15.086 1.763 ( 1.378 0.795 9.777) 9.906 1.926 ( 4.175 2.410 -5.581) 7.375 1.946 ( 1.505 0.869 -7.651) 7.846 1.989 ( 6.377 3.682 -9.198) 11.782 2.545 ( -0.271 -0.157 6.756) 6.763 2.589 ( 0.487 0.281 2.555) 2.616 2.759 ( -0.762 -0.440 3.381) 3.494 2.825 ( -0.536 -0.309 -2.773) 2.841 3.020 ( 4.101 2.368 -1.754) 5.050 3.028 ( 4.311 2.489 -2.246) 5.462 3.524 ( -2.636 -1.522 3.513) 4.648 3.590 ( -4.965 -2.867 -2.638) 6.311 3.631 ( 0.800 0.462 0.488) 1.044 3.641 ( 0.806 0.465 -0.458) 1.037 4.662 ( 8.294 4.789 10.113) 13.928 4.885 ( -2.953 -1.705 -10.511) 11.050 5.176 ( -12.422 -7.172 5.605) 15.400 5.302 ( -4.987 -2.879 -6.681) 8.820 ======================= Grid point 184 (25/42) ======================= q-point: ( 0.23 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.016 ( 7.802 4.504 7.394) 11.654 1.072 ( 10.622 6.132 0.635) 12.281 1.107 ( 9.820 5.669 -1.697) 11.465 1.200 ( 7.021 4.053 -10.323) 13.125 1.570 ( 6.726 3.883 7.933) 11.102 1.772 ( 7.758 4.479 -9.647) 13.165 1.801 ( 1.832 1.058 10.259) 10.474 1.988 ( 2.099 1.212 -7.786) 8.155 2.032 ( 10.476 6.048 7.553) 14.261 2.243 ( 12.608 7.279 -11.069) 18.288 2.479 ( -4.073 -2.351 11.113) 12.067 2.647 ( 1.957 1.130 -4.454) 4.995 2.749 ( -0.121 -0.070 3.726) 3.729 2.821 ( 0.144 0.083 -3.000) 3.005 3.092 ( 2.795 1.614 0.966) 3.369 3.126 ( 3.095 1.787 -1.481) 3.869 3.467 ( -2.106 -1.216 0.031) 2.432 3.470 ( -4.867 -2.810 -0.192) 5.623 3.642 ( 0.123 0.071 0.484) 0.504 3.653 ( 0.111 0.064 -0.451) 0.469 4.542 ( -13.873 -8.010 3.512) 16.400 4.624 ( -16.343 -9.435 -4.221) 19.337 5.098 ( 3.037 1.754 9.264) 9.906 5.285 ( 0.888 0.512 -8.177) 8.241 ======================= Grid point 185 (26/42) ======================= q-point: ( 0.31 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.169 ( 5.701 3.292 6.499) 9.251 1.282 ( 8.187 4.727 2.816) 9.864 1.327 ( 7.660 4.423 -0.892) 8.890 1.339 ( 5.108 2.949 -10.058) 11.660 1.755 ( 9.089 5.247 0.953) 10.538 1.845 ( 2.000 1.155 10.847) 11.090 1.856 ( 4.191 2.419 -8.289) 9.598 2.040 ( 2.258 1.304 -8.053) 8.465 2.280 ( 0.745 0.430 9.528) 9.567 2.443 ( 7.368 4.254 11.008) 13.912 2.467 ( 4.297 2.481 -6.913) 8.509 2.677 ( 2.870 1.657 -10.396) 10.912 2.751 ( 0.164 0.095 4.045) 4.049 2.828 ( 0.316 0.182 -3.170) 3.191 3.132 ( -0.928 -0.536 1.678) 1.991 3.141 ( -1.389 -0.802 0.670) 1.738 3.401 ( -0.371 -0.214 1.571) 1.628 3.439 ( -0.170 -0.098 -1.817) 1.827 3.635 ( -0.740 -0.427 0.478) 0.979 3.645 ( -0.737 -0.425 -0.444) 0.960 4.216 ( -12.757 -7.365 1.662) 14.824 4.254 ( -14.367 -8.295 -1.856) 16.693 5.160 ( 1.710 0.987 6.416) 6.713 5.294 ( -0.336 -0.194 -5.903) 5.916 ======================= Grid point 186 (27/42) ======================= q-point: ( 0.38 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.283 ( 4.088 2.360 5.792) 7.471 1.436 ( 3.721 2.149 0.507) 4.327 1.439 ( 3.552 2.051 -9.645) 10.481 1.467 ( 6.416 3.704 -2.045) 7.685 1.891 ( 1.835 1.059 11.094) 11.295 1.938 ( 7.014 4.050 4.453) 9.242 2.051 ( 9.982 5.763 -3.631) 12.084 2.088 ( 1.784 1.030 -8.193) 8.448 2.185 ( -3.371 -1.946 7.417) 8.376 2.391 ( -6.398 -3.694 -10.449) 12.797 2.703 ( 7.780 4.492 4.389) 9.998 2.752 ( -0.041 -0.023 4.193) 4.194 2.760 ( 1.788 1.032 0.285) 2.084 2.832 ( -0.031 -0.018 -3.208) 3.208 3.051 ( -4.490 -2.592 0.793) 5.244 3.090 ( -1.908 -1.102 -2.354) 3.224 3.442 ( 2.877 1.661 0.505) 3.361 3.448 ( 0.594 0.343 -0.042) 0.687 3.612 ( -1.121 -0.647 0.497) 1.387 3.622 ( -1.091 -0.630 -0.462) 1.341 3.990 ( -6.723 -3.882 0.332) 7.770 3.998 ( -7.590 -4.382 -0.413) 8.773 5.177 ( -0.040 -0.023 4.475) 4.475 5.272 ( -1.394 -0.805 -4.244) 4.539 ======================= Grid point 187 (28/42) ======================= q-point: ( 0.46 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.350 ( 1.572 0.908 5.398) 5.695 1.486 ( 1.183 0.683 4.745) 4.938 1.498 ( 1.377 0.795 -9.237) 9.373 1.600 ( 3.776 2.180 -5.775) 7.236 1.924 ( 0.842 0.486 10.921) 10.964 1.998 ( -0.210 -0.121 5.477) 5.483 2.117 ( 0.705 0.407 -8.092) 8.133 2.159 ( -4.048 -2.337 -9.887) 10.936 2.223 ( 6.434 3.715 7.834) 10.796 2.364 ( 5.190 2.996 -5.049) 7.836 2.644 ( -5.506 -3.179 7.285) 9.669 2.750 ( -0.116 -0.067 4.193) 4.195 2.787 ( -3.321 -1.917 -6.704) 7.723 2.828 ( -0.142 -0.082 -3.176) 3.181 3.010 ( 1.786 1.031 4.891) 5.308 3.090 ( 0.507 0.293 -2.905) 2.964 3.455 ( 0.160 0.093 2.106) 2.114 3.500 ( 1.605 0.927 -2.034) 2.752 3.591 ( -0.549 -0.317 0.532) 0.827 3.602 ( -0.526 -0.304 -0.494) 0.783 3.894 ( -1.894 -1.094 0.074) 2.188 3.895 ( -1.869 -1.079 -0.018) 2.158 5.170 ( -0.292 -0.169 3.388) 3.405 5.243 ( -0.826 -0.477 -3.262) 3.399 ======================= Grid point 195 (29/42) ======================= q-point: ( 0.08 0.08 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.694 ( 6.686 11.581 5.423) 14.430 0.766 ( 7.693 13.325 8.330) 17.496 0.891 ( 2.627 4.550 -11.849) 12.962 0.987 ( 6.419 11.119 -10.770) 16.758 1.235 ( 11.411 19.765 11.854) 25.717 1.368 ( 4.610 7.986 3.187) 9.756 1.773 ( 1.569 2.718 8.697) 9.246 1.858 ( 2.081 3.604 -5.331) 6.763 1.944 ( 0.098 0.169 -7.130) 7.133 1.977 ( 3.232 5.598 -7.499) 9.900 2.537 ( 0.689 1.194 6.778) 6.917 2.585 ( -1.998 -3.460 1.972) 4.456 2.780 ( 0.104 0.180 3.632) 3.638 2.851 ( 0.202 0.349 -3.013) 3.040 2.975 ( 2.136 3.700 -1.934) 4.690 2.976 ( 1.332 2.306 -1.893) 3.267 3.538 ( -1.548 -2.680 4.241) 5.250 3.616 ( -2.438 -4.223 -3.024) 5.737 3.638 ( 1.140 1.975 0.505) 2.336 3.648 ( 0.906 1.570 -0.433) 1.864 4.600 ( 9.263 16.043 12.404) 22.295 4.900 ( 3.332 5.771 -15.374) 16.756 5.279 ( -6.731 -11.658 0.759) 13.483 5.300 ( -8.613 -14.919 -1.047) 17.258 ======================= Grid point 196 (30/42) ======================= q-point: ( 0.15 0.08 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 254 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.945 ( 7.987 9.805 -0.353) 12.651 0.985 ( 7.761 5.138 -3.688) 10.012 1.050 ( 5.419 13.205 5.671) 15.359 1.208 ( 4.828 12.041 -8.680) 15.609 1.513 ( 4.661 5.701 10.589) 12.897 1.619 ( 13.642 8.363 5.683) 16.981 1.823 ( -0.496 4.060 -0.417) 4.112 1.939 ( 3.722 0.182 -8.979) 9.722 1.977 ( 4.192 4.854 3.126) 7.135 2.151 ( 7.414 9.815 -10.208) 15.984 2.511 ( -2.755 -2.549 7.389) 8.287 2.583 ( 3.091 -1.162 1.097) 3.480 2.788 ( -0.875 2.743 2.884) 4.075 2.846 ( -1.462 1.468 -2.685) 3.392 3.027 ( 3.593 -0.738 0.122) 3.670 3.058 ( 4.165 0.575 -2.066) 4.685 3.481 ( -1.812 -2.489 1.066) 3.258 3.507 ( -4.055 -5.068 -1.333) 6.626 3.679 ( -0.008 3.963 0.336) 3.977 3.685 ( 0.020 3.498 -0.185) 3.503 4.725 ( -11.280 -0.276 4.899) 12.301 4.820 ( -10.403 -4.958 -4.140) 12.245 5.018 ( -0.737 -7.756 7.754) 10.992 5.192 ( 4.805 -6.927 -8.312) 11.839 ======================= Grid point 197 (31/42) ======================= q-point: ( 0.23 0.08 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 254 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.136 ( 7.173 5.184 2.276) 9.138 1.179 ( 5.165 7.137 -0.556) 8.827 1.284 ( 6.256 12.634 2.305) 14.285 1.379 ( 1.907 10.546 -6.103) 12.333 1.657 ( 6.880 4.878 5.262) 9.941 1.825 ( 3.118 1.877 9.283) 9.971 1.843 ( 1.622 3.813 -11.206) 11.947 2.005 ( 3.292 1.235 -6.560) 7.443 2.196 ( 10.077 7.121 7.524) 14.453 2.376 ( 7.992 3.986 -8.379) 12.247 2.444 ( -2.156 0.992 11.497) 11.740 2.617 ( 3.233 -1.744 -3.821) 5.300 2.794 ( -1.791 2.929 -0.276) 3.444 2.840 ( -0.689 1.227 -3.023) 3.335 3.080 ( 3.381 -2.177 2.077) 4.526 3.102 ( 2.083 -3.569 -0.071) 4.133 3.402 ( -2.195 -3.454 1.389) 4.322 3.437 ( -0.703 -1.449 -1.829) 2.437 3.712 ( -2.002 6.022 0.131) 6.348 3.714 ( -1.779 5.365 -0.061) 5.652 4.411 ( -16.638 -5.170 2.856) 17.655 4.476 ( -18.523 -5.480 -3.214) 19.582 5.065 ( 6.638 -4.722 7.292) 10.933 5.214 ( 4.889 -7.003 -6.532) 10.753 ======================= Grid point 198 (32/42) ======================= q-point: ( 0.31 0.08 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 254 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.294 ( 3.477 6.519 1.394) 7.518 1.312 ( 5.111 3.720 -0.761) 6.367 1.468 ( 0.595 9.649 1.133) 9.734 1.514 ( 3.102 10.541 -3.635) 11.574 1.826 ( 5.840 2.237 1.941) 6.548 1.886 ( 3.791 3.652 4.545) 6.955 1.950 ( 5.359 4.809 -1.890) 7.444 2.066 ( 2.269 1.295 -5.773) 6.337 2.263 ( -3.518 0.789 7.033) 7.903 2.429 ( -0.868 -3.619 -7.207) 8.112 2.588 ( 9.550 6.711 5.791) 13.029 2.681 ( 3.875 -0.864 -0.095) 3.972 2.797 ( -0.787 2.730 -0.322) 2.859 2.852 ( -0.144 2.315 -3.196) 3.949 3.052 ( -1.158 -5.635 0.955) 5.832 3.074 ( -0.263 -4.867 -0.587) 4.909 3.387 ( 2.590 -0.707 2.099) 3.408 3.428 ( 0.673 -0.911 -1.640) 1.993 3.729 ( -4.171 8.107 0.089) 9.118 3.731 ( -3.725 7.549 -0.044) 8.418 4.109 ( -12.390 -3.355 1.337) 12.905 4.140 ( -14.299 -3.256 -1.455) 14.737 5.109 ( 4.543 -5.696 4.987) 8.829 5.214 ( 3.368 -7.324 -4.669) 9.315 ======================= Grid point 199 (33/42) ======================= q-point: ( 0.38 0.08 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 254 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.383 ( -0.121 6.483 0.723) 6.524 1.432 ( 4.785 1.392 -4.063) 6.429 1.546 ( -2.037 9.383 5.133) 10.888 1.657 ( 0.628 10.343 -4.577) 11.328 1.885 ( 0.871 -0.637 6.407) 6.498 1.966 ( 3.627 -1.290 5.058) 6.356 2.077 ( 1.767 -1.948 -6.300) 6.826 2.142 ( 5.673 2.353 -4.312) 7.504 2.250 ( -1.472 7.234 3.645) 8.233 2.371 ( -1.579 3.266 -5.260) 6.390 2.667 ( -1.833 -5.445 5.751) 8.129 2.730 ( 2.226 -0.974 -0.825) 2.566 2.834 ( -1.018 2.699 0.062) 2.885 2.845 ( -1.607 0.669 -1.776) 2.487 2.977 ( 0.322 -2.981 3.421) 4.549 3.054 ( 2.190 -1.708 -3.456) 4.434 3.431 ( 0.808 -1.291 0.570) 1.626 3.440 ( 3.437 -0.825 -0.316) 3.548 3.731 ( -6.171 10.218 0.281) 11.940 3.737 ( -5.742 10.055 -0.270) 11.582 3.931 ( -5.103 -2.099 0.253) 5.524 3.937 ( -6.132 -1.813 -0.306) 6.401 5.116 ( 3.168 -5.728 3.544) 7.444 5.192 ( 2.545 -6.825 -3.396) 8.037 ======================= Grid point 200 (34/42) ======================= q-point: ( 0.46 0.08 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.417 ( -2.914 5.048 1.475) 6.013 1.484 ( 0.254 -0.439 -4.750) 4.777 1.577 ( -4.852 8.404 7.858) 12.487 1.719 ( -3.763 6.518 -3.764) 8.415 1.889 ( 0.369 -0.640 1.887) 2.026 2.006 ( 3.631 -6.290 -7.550) 10.476 2.027 ( -0.507 0.878 6.616) 6.693 2.192 ( -0.686 1.188 -8.966) 9.070 2.323 ( -2.993 5.185 6.890) 9.128 2.421 ( -1.114 1.929 -1.805) 2.867 2.580 ( 1.172 -2.030 7.958) 8.296 2.738 ( 1.153 -1.997 -7.448) 7.797 2.784 ( -1.684 2.916 2.135) 3.987 2.837 ( -0.446 0.772 -2.569) 2.719 3.003 ( 1.019 -1.764 4.161) 4.633 3.068 ( 1.220 -2.114 -2.192) 3.281 3.430 ( 0.780 -1.351 1.840) 2.413 3.469 ( 1.309 -2.268 -1.646) 3.093 3.727 ( -6.475 11.215 0.324) 12.954 3.736 ( -6.678 11.567 -0.500) 13.366 3.879 ( 0.270 -0.468 0.081) 0.547 3.879 ( 0.289 -0.501 0.016) 0.578 5.113 ( 3.103 -5.374 2.995) 6.890 5.177 ( 3.389 -5.869 -2.894) 7.369 ======================= Grid point 210 (35/42) ======================= q-point: ( 0.15 0.15 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.087 ( 3.724 6.450 2.855) 7.977 1.174 ( 6.635 11.493 -5.531) 14.377 1.317 ( 6.962 12.058 3.658) 14.396 1.428 ( 5.035 8.721 -6.717) 12.105 1.636 ( 4.245 7.353 10.405) 13.429 1.759 ( 3.548 6.146 9.653) 11.981 1.875 ( 1.382 2.393 -10.882) 11.227 1.977 ( 1.882 3.259 -9.942) 10.631 2.134 ( 5.009 8.675 8.616) 13.212 2.324 ( 3.560 6.165 -8.808) 11.325 2.480 ( 0.315 0.545 8.283) 8.307 2.565 ( 0.303 0.525 0.662) 0.898 2.831 ( 0.735 1.274 -1.063) 1.814 2.859 ( -0.132 -0.229 -2.920) 2.932 3.028 ( 0.530 0.918 2.217) 2.457 3.055 ( -0.545 -0.945 -0.350) 1.146 3.405 ( -2.637 -4.567 1.097) 5.387 3.436 ( -1.044 -1.809 -1.710) 2.700 3.760 ( 2.112 3.659 0.219) 4.230 3.764 ( 2.346 4.063 -0.148) 4.694 4.591 ( -6.715 -11.631 3.322) 13.835 4.668 ( -6.632 -11.487 -3.902) 13.826 4.928 ( -0.516 -0.895 7.617) 7.686 5.080 ( -2.066 -3.579 -6.576) 7.766 ======================= Grid point 211 (36/42) ======================= q-point: ( 0.23 0.15 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 254 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.228 ( 5.447 5.305 5.903) 9.626 1.371 ( 2.663 9.977 -7.534) 12.783 1.536 ( 5.286 10.786 -0.578) 12.026 1.556 ( 2.692 7.062 -2.959) 8.116 1.786 ( 2.363 9.152 7.292) 11.938 1.885 ( 4.981 4.051 7.090) 9.565 1.944 ( 1.487 5.036 -2.885) 5.991 2.039 ( 1.628 2.141 -6.709) 7.228 2.265 ( 4.142 2.498 0.857) 4.913 2.375 ( 1.642 -2.482 -6.512) 7.159 2.557 ( 2.159 7.330 7.696) 10.845 2.617 ( 4.018 1.075 3.318) 5.320 2.838 ( -1.605 0.569 -2.920) 3.380 2.868 ( 1.196 2.454 -4.374) 5.156 3.014 ( -0.944 -4.289 1.550) 4.658 3.032 ( 0.118 -2.598 0.076) 2.602 3.352 ( 0.317 -1.268 3.057) 3.325 3.412 ( -0.087 -1.040 -2.531) 2.738 3.822 ( 0.227 4.454 0.388) 4.476 3.831 ( 0.144 4.724 -0.402) 4.743 4.336 ( -13.427 -2.266 2.752) 13.893 4.397 ( -14.758 -2.629 -2.877) 15.263 4.937 ( 5.322 -7.209 4.989) 10.256 5.040 ( 4.339 -9.411 -4.557) 11.321 ======================= Grid point 212 (37/42) ======================= q-point: ( 0.31 0.15 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 254 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.373 ( 5.057 5.026 3.023) 7.744 1.470 ( -0.421 7.900 -6.722) 10.381 1.674 ( 3.061 8.027 3.029) 9.109 1.722 ( 1.213 10.830 0.383) 10.905 1.906 ( 3.168 5.884 0.676) 6.717 1.926 ( -2.168 0.650 7.174) 7.523 2.027 ( 1.356 2.795 -5.293) 6.137 2.089 ( 2.823 1.411 -3.280) 4.552 2.344 ( 0.505 5.867 0.761) 5.938 2.400 ( 1.882 2.459 -2.103) 3.743 2.607 ( 2.310 -3.901 4.325) 6.266 2.693 ( 0.712 0.269 -2.770) 2.873 2.842 ( 1.860 2.447 2.847) 4.189 2.897 ( -1.596 1.578 -1.018) 2.465 2.953 ( -1.249 -3.438 1.501) 3.954 3.008 ( 0.901 -1.327 -3.031) 3.429 3.378 ( 2.822 -0.353 1.554) 3.241 3.408 ( 0.587 -0.934 -1.220) 1.645 3.856 ( -0.578 2.953 0.321) 3.026 3.863 ( -0.747 2.910 -0.361) 3.026 4.115 ( -13.171 5.246 1.604) 14.268 4.150 ( -14.366 5.174 -1.618) 15.355 4.956 ( 6.067 -9.020 3.389) 11.387 5.027 ( 5.535 -10.632 -3.214) 12.410 ======================= Grid point 213 (38/42) ======================= q-point: ( 0.38 0.15 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 254 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.485 ( 1.822 4.432 1.122) 4.922 1.530 ( -2.143 7.308 -3.390) 8.336 1.769 ( 0.480 1.258 3.973) 4.195 1.790 ( -4.290 12.340 4.813) 13.923 1.897 ( -0.104 -1.897 -0.556) 1.979 1.987 ( -2.324 5.820 -7.416) 9.709 2.029 ( -0.080 1.147 3.110) 3.316 2.162 ( 3.382 -0.294 -7.103) 7.872 2.420 ( -0.496 7.746 3.934) 8.702 2.445 ( 0.127 1.707 1.534) 2.299 2.574 ( -1.644 -2.262 5.701) 6.350 2.693 ( 1.818 -1.613 -6.184) 6.644 2.855 ( -1.888 0.655 -0.389) 2.036 2.865 ( -3.198 1.095 -1.062) 3.543 2.953 ( 2.787 -0.542 4.204) 5.072 3.020 ( 1.467 -1.934 -2.386) 3.403 3.411 ( 0.582 -0.626 0.376) 0.934 3.418 ( 2.304 -1.504 -0.285) 2.766 3.863 ( 0.161 -0.399 0.067) 0.436 3.865 ( -0.027 -0.690 -0.097) 0.697 4.002 ( -9.982 12.304 0.890) 15.869 4.022 ( -10.569 12.390 -0.919) 16.311 4.965 ( 5.422 -9.001 2.506) 10.803 5.018 ( 5.441 -10.032 -2.419) 11.666 ======================= Grid point 224 (39/42) ======================= q-point: ( 0.23 0.23 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.354 ( 4.119 7.134 6.381) 10.420 1.525 ( 2.918 5.055 -10.466) 11.983 1.670 ( 3.796 6.575 2.474) 7.985 1.747 ( 5.633 9.758 -5.385) 12.488 1.909 ( -0.206 -0.356 6.433) 6.446 2.020 ( 2.746 4.757 7.040) 8.930 2.031 ( 5.892 10.205 6.117) 13.277 2.060 ( -0.845 -1.463 -8.440) 8.607 2.351 ( 1.007 1.745 -3.275) 3.845 2.383 ( 4.940 8.556 -6.870) 12.033 2.565 ( -1.950 -3.377 8.166) 9.049 2.695 ( 0.569 0.986 -2.847) 3.066 2.836 ( 0.921 1.595 3.088) 3.595 2.919 ( 0.873 1.512 -3.054) 3.518 2.947 ( -1.280 -2.218 2.606) 3.654 2.995 ( -0.428 -0.742 -2.461) 2.606 3.349 ( 0.423 0.734 2.428) 2.572 3.397 ( -0.250 -0.434 -2.003) 2.065 3.871 ( 0.255 0.441 0.335) 0.609 3.879 ( 0.096 0.166 -0.356) 0.404 4.308 ( -0.670 -1.160 2.114) 2.502 4.355 ( -1.230 -2.131 -2.187) 3.292 4.815 ( -2.437 -4.222 3.143) 5.800 4.881 ( -3.304 -5.723 -2.920) 7.225 ======================= Grid point 225 (40/42) ======================= q-point: ( 0.31 0.23 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 254 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.498 ( 3.570 7.198 4.115) 9.027 1.610 ( 1.450 5.952 -6.735) 9.105 1.771 ( 0.577 1.689 4.024) 4.402 1.854 ( -1.741 2.942 -1.055) 3.578 1.949 ( 4.306 -1.755 -1.320) 4.834 2.017 ( 1.534 2.289 -4.595) 5.358 2.079 ( -3.271 9.153 6.652) 11.778 2.192 ( 0.445 7.091 -2.417) 7.505 2.437 ( 2.248 1.116 3.121) 4.005 2.484 ( -0.351 3.603 4.410) 5.705 2.553 ( 2.494 -0.008 -2.192) 3.320 2.653 ( 1.161 -0.420 -7.053) 7.160 2.885 ( -1.308 0.159 2.160) 2.531 2.916 ( 0.299 -0.216 3.428) 3.448 2.920 ( -0.718 0.752 0.018) 1.040 2.991 ( 0.180 -0.992 -2.616) 2.804 3.370 ( 1.506 -0.549 1.274) 2.047 3.397 ( 0.496 -0.059 -1.251) 1.347 3.856 ( -0.633 -1.621 0.087) 1.743 3.858 ( -0.769 -1.865 -0.089) 2.019 4.279 ( -7.003 9.154 1.307) 11.599 4.308 ( -7.586 8.885 -1.382) 11.765 4.773 ( 4.491 -8.692 2.072) 10.001 4.817 ( 4.299 -9.805 -1.964) 10.885 ======================= Grid point 226 (41/42) ======================= q-point: ( 0.38 0.23 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.578 ( -2.595 4.494 4.187) 6.668 1.674 ( -3.489 6.044 -4.524) 8.317 1.781 ( -0.225 0.390 -1.516) 1.581 1.795 ( 1.762 -3.051 -2.783) 4.490 1.980 ( 0.643 -1.113 7.016) 7.133 2.078 ( -6.675 11.561 5.086) 14.286 2.112 ( 1.987 -3.441 -3.639) 5.388 2.216 ( -5.220 9.042 -6.352) 12.221 2.420 ( 1.033 -1.789 2.980) 3.626 2.493 ( 0.532 -0.922 -2.933) 3.120 2.597 ( -1.615 2.798 6.418) 7.185 2.716 ( -1.894 3.281 -5.285) 6.503 2.841 ( 1.143 -1.979 2.223) 3.189 2.893 ( -0.387 0.671 -2.624) 2.735 2.943 ( -0.050 0.087 2.770) 2.771 2.983 ( 0.649 -1.124 -1.107) 1.706 3.381 ( 1.141 -1.977 1.068) 2.521 3.403 ( 0.132 -0.228 -1.003) 1.037 3.838 ( 0.757 -1.310 0.039) 1.514 3.839 ( 0.730 -1.264 -0.035) 1.460 4.262 ( -7.227 12.517 1.014) 14.489 4.285 ( -7.328 12.693 -1.087) 14.697 4.772 ( 5.689 -9.854 1.719) 11.508 4.809 ( 6.058 -10.493 -1.646) 12.227 ======================= Grid point 239 (42/42) ======================= q-point: ( 0.31 0.31 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.10e-04 1.10e-04 0.00e+00 0.00e+00 2.84e-04 2.84e-04 2.84e-04 2.84e-04 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.644 ( 3.819 6.615 2.301) 7.978 1.694 ( 1.454 2.518 -2.258) 3.682 1.784 ( 0.399 0.691 -0.777) 1.114 1.806 ( -3.638 -6.301 -2.835) 7.809 1.973 ( 1.990 3.446 4.385) 5.922 2.041 ( -0.236 -0.409 -1.105) 1.202 2.226 ( 2.154 3.730 4.916) 6.536 2.314 ( 2.602 4.507 -3.381) 6.206 2.405 ( -0.468 -0.810 0.344) 0.996 2.469 ( -2.928 -5.071 -4.930) 7.655 2.625 ( 4.564 7.905 6.444) 11.173 2.755 ( 3.563 6.171 -6.222) 9.460 2.812 ( -3.265 -5.655 3.918) 7.615 2.889 ( -2.034 -3.524 -3.471) 5.348 2.953 ( 1.007 1.745 2.496) 3.208 2.977 ( -0.161 -0.279 -0.025) 0.323 3.357 ( -0.322 -0.558 2.076) 2.174 3.401 ( 0.117 0.203 -1.970) 1.984 3.823 ( -0.640 -1.108 0.072) 1.281 3.825 ( -0.590 -1.021 -0.043) 1.180 4.439 ( 3.568 6.179 1.096) 7.219 4.464 ( 3.487 6.040 -1.254) 7.087 4.624 ( -3.341 -5.787 1.412) 6.830 4.653 ( -3.577 -6.195 -1.260) 7.263 =================== 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/12168 10.0 429.855 429.855 217.442 0.000 0.000 -0.000 3/12168 20.0 78.523 78.523 38.490 0.000 0.000 0.000 3/12168 30.0 34.964 34.964 18.135 0.000 0.000 0.000 3/12168 40.0 21.735 21.735 12.049 0.000 0.000 0.000 3/12168 50.0 15.765 15.765 9.105 0.000 0.000 0.000 3/12168 60.0 12.449 12.449 7.354 0.000 0.000 0.000 3/12168 70.0 10.346 10.346 6.189 0.000 -0.000 0.000 3/12168 80.0 8.886 8.886 5.354 0.000 -0.000 0.000 3/12168 90.0 7.808 7.808 4.725 0.000 -0.000 0.000 3/12168 100.0 6.975 6.975 4.232 0.000 -0.000 0.000 3/12168 110.0 6.310 6.310 3.835 0.000 -0.000 0.000 3/12168 120.0 5.765 5.765 3.508 0.000 -0.000 0.000 3/12168 130.0 5.310 5.310 3.234 0.000 -0.000 0.000 3/12168 140.0 4.922 4.922 3.000 0.000 -0.000 0.000 3/12168 150.0 4.589 4.589 2.798 0.000 -0.000 0.000 3/12168 160.0 4.299 4.299 2.622 0.000 -0.000 0.000 3/12168 170.0 4.044 4.044 2.467 0.000 -0.000 0.000 3/12168 180.0 3.818 3.818 2.330 0.000 -0.000 0.000 3/12168 190.0 3.616 3.616 2.207 0.000 -0.000 0.000 3/12168 200.0 3.435 3.435 2.097 0.000 -0.000 0.000 3/12168 210.0 3.271 3.271 1.997 0.000 -0.000 0.000 3/12168 220.0 3.122 3.122 1.907 0.000 -0.000 0.000 3/12168 230.0 2.986 2.986 1.824 0.000 -0.000 0.000 3/12168 240.0 2.862 2.862 1.748 0.000 -0.000 0.000 3/12168 250.0 2.748 2.748 1.679 0.000 -0.000 0.000 3/12168 260.0 2.642 2.642 1.614 0.000 -0.000 0.000 3/12168 270.0 2.544 2.544 1.555 0.000 -0.000 0.000 3/12168 280.0 2.454 2.454 1.500 0.000 -0.000 0.000 3/12168 290.0 2.369 2.369 1.448 0.000 -0.000 0.000 3/12168 300.0 2.291 2.291 1.400 0.000 -0.000 0.000 3/12168 310.0 2.217 2.217 1.355 0.000 -0.000 0.000 3/12168 320.0 2.148 2.148 1.313 0.000 -0.000 0.000 3/12168 330.0 2.083 2.083 1.274 0.000 -0.000 0.000 3/12168 340.0 2.022 2.022 1.236 0.000 -0.000 0.000 3/12168 350.0 1.964 1.964 1.201 0.000 -0.000 0.000 3/12168 360.0 1.910 1.910 1.168 0.000 -0.000 0.000 3/12168 370.0 1.859 1.859 1.137 0.000 -0.000 0.000 3/12168 380.0 1.810 1.810 1.107 0.000 -0.000 0.000 3/12168 390.0 1.764 1.764 1.079 0.000 -0.000 0.000 3/12168 400.0 1.720 1.720 1.052 0.000 -0.000 0.000 3/12168 410.0 1.678 1.678 1.026 0.000 -0.000 0.000 3/12168 420.0 1.638 1.638 1.002 0.000 -0.000 0.000 3/12168 430.0 1.600 1.600 0.979 0.000 -0.000 0.000 3/12168 440.0 1.564 1.564 0.957 0.000 -0.000 0.000 3/12168 450.0 1.529 1.529 0.936 0.000 -0.000 0.000 3/12168 460.0 1.496 1.496 0.915 0.000 -0.000 0.000 3/12168 470.0 1.465 1.465 0.896 0.000 -0.000 0.000 3/12168 480.0 1.434 1.434 0.877 0.000 -0.000 0.000 3/12168 490.0 1.405 1.405 0.860 0.000 -0.000 0.000 3/12168 500.0 1.377 1.377 0.842 0.000 -0.000 0.000 3/12168 510.0 1.350 1.350 0.826 0.000 -0.000 0.000 3/12168 520.0 1.324 1.324 0.810 0.000 -0.000 0.000 3/12168 530.0 1.299 1.299 0.795 0.000 -0.000 0.000 3/12168 540.0 1.275 1.275 0.780 0.000 -0.000 0.000 3/12168 550.0 1.252 1.252 0.766 0.000 -0.000 0.000 3/12168 560.0 1.230 1.230 0.753 0.000 -0.000 0.000 3/12168 570.0 1.208 1.208 0.740 0.000 -0.000 0.000 3/12168 580.0 1.188 1.188 0.727 0.000 -0.000 0.000 3/12168 590.0 1.168 1.168 0.715 0.000 -0.000 0.000 3/12168 600.0 1.148 1.148 0.703 0.000 -0.000 0.000 3/12168 610.0 1.130 1.130 0.691 0.000 -0.000 0.000 3/12168 620.0 1.111 1.111 0.680 0.000 -0.000 0.000 3/12168 630.0 1.094 1.094 0.669 0.000 -0.000 0.000 3/12168 640.0 1.077 1.077 0.659 0.000 -0.000 0.000 3/12168 650.0 1.060 1.060 0.649 0.000 -0.000 0.000 3/12168 660.0 1.044 1.044 0.639 0.000 -0.000 0.000 3/12168 670.0 1.029 1.029 0.630 0.000 -0.000 0.000 3/12168 680.0 1.014 1.014 0.620 0.000 -0.000 0.000 3/12168 690.0 0.999 0.999 0.611 0.000 -0.000 0.000 3/12168 700.0 0.985 0.985 0.603 0.000 -0.000 0.000 3/12168 710.0 0.971 0.971 0.594 0.000 -0.000 0.000 3/12168 720.0 0.957 0.957 0.586 0.000 -0.000 0.000 3/12168 730.0 0.944 0.944 0.578 0.000 -0.000 0.000 3/12168 740.0 0.932 0.932 0.570 0.000 -0.000 0.000 3/12168 750.0 0.919 0.919 0.563 0.000 -0.000 0.000 3/12168 760.0 0.907 0.907 0.555 0.000 -0.000 0.000 3/12168 770.0 0.895 0.895 0.548 0.000 -0.000 0.000 3/12168 780.0 0.884 0.884 0.541 0.000 -0.000 0.000 3/12168 790.0 0.873 0.873 0.534 0.000 -0.000 0.000 3/12168 800.0 0.862 0.862 0.528 0.000 -0.000 0.000 3/12168 810.0 0.851 0.851 0.521 0.000 -0.000 0.000 3/12168 820.0 0.841 0.841 0.515 0.000 -0.000 0.000 3/12168 830.0 0.831 0.831 0.509 0.000 -0.000 0.000 3/12168 840.0 0.821 0.821 0.503 0.000 -0.000 0.000 3/12168 850.0 0.811 0.811 0.497 0.000 -0.000 0.000 3/12168 860.0 0.802 0.802 0.491 0.000 -0.000 0.000 3/12168 870.0 0.793 0.793 0.485 0.000 -0.000 0.000 3/12168 880.0 0.784 0.784 0.480 0.000 -0.000 0.000 3/12168 890.0 0.775 0.775 0.475 0.000 -0.000 0.000 3/12168 900.0 0.766 0.766 0.469 0.000 -0.000 0.000 3/12168 910.0 0.758 0.758 0.464 0.000 -0.000 0.000 3/12168 920.0 0.750 0.750 0.459 0.000 -0.000 0.000 3/12168 930.0 0.742 0.742 0.454 0.000 -0.000 0.000 3/12168 940.0 0.734 0.734 0.449 0.000 -0.000 0.000 3/12168 950.0 0.726 0.726 0.445 0.000 -0.000 0.000 3/12168 960.0 0.719 0.719 0.440 0.000 -0.000 0.000 3/12168 970.0 0.711 0.711 0.436 0.000 -0.000 0.000 3/12168 980.0 0.704 0.704 0.431 0.000 -0.000 0.000 3/12168 990.0 0.697 0.697 0.427 0.000 -0.000 0.000 3/12168 1000.0 0.690 0.690 0.423 0.000 -0.000 0.000 3/12168 Thermal conductivity related properties were written into "kappa-m13133.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 02:22:11]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|