----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-09 04:33:49]------------------------- 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.050675551774594 0.000000000000000 0.000000000000000 b -2.025337775887297 3.507987930325347 0.000000000000000 c 0.000000000000000 0.000000000000000 15.840398120000000 Atomic positions (fractional): *1 Tb 0.00000000000000 0.00000000000000 0.00000000000000 158.925 2 Tb 0.00000000000000 0.00000000000000 0.50000000000000 158.925 *3 Cs 0.33333333333333 0.66666666666667 0.75000000000000 132.905 4 Cs 0.66666666666667 0.33333333333333 0.25000000000000 132.905 *5 S 0.33333333333333 0.66666666666667 0.40770196528995 32.065 6 S 0.66666666666667 0.33333333333333 0.90770196528995 32.065 7 S 0.66666666666667 0.33333333333333 0.59229803471005 32.065 8 S 0.33333333333333 0.66666666666667 0.09229803471005 32.065 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.050675551774594 0.000000000000000 0.000000000000000 b -2.025337775887297 3.507987930325347 0.000000000000000 c 0.000000000000000 0.000000000000000 15.840398120000000 Atomic positions (fractional): *1 Tb 0.00000000000000 0.00000000000000 0.00000000000000 158.925 > 1 2 Tb 0.00000000000000 0.00000000000000 0.50000000000000 158.925 > 2 *3 Cs 0.33333333333333 0.66666666666667 0.75000000000000 132.905 > 3 4 Cs 0.66666666666667 0.33333333333333 0.25000000000000 132.905 > 4 *5 S 0.33333333333333 0.66666666666667 0.40770196528995 32.065 > 5 6 S 0.66666666666667 0.33333333333333 0.90770196528995 32.065 > 6 7 S 0.66666666666667 0.33333333333333 0.59229803471005 32.065 > 7 8 S 0.33333333333333 0.66666666666667 0.09229803471005 32.065 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.202702207098376 0.000000000000000 0.000000000000000 b -8.101351103549188 14.031951721301388 0.000000000000000 c 0.000000000000000 0.000000000000000 15.840398120000000 Atomic positions (fractional): *1 Tb 0.00000000000000 0.00000000000000 0.00000000000000 158.925 > 1 2 Tb 0.25000000000000 0.00000000000000 0.00000000000000 158.925 > 1 3 Tb 0.50000000000000 0.00000000000000 0.00000000000000 158.925 > 1 4 Tb 0.75000000000000 0.00000000000000 0.00000000000000 158.925 > 1 5 Tb 0.00000000000000 0.25000000000000 0.00000000000000 158.925 > 1 6 Tb 0.25000000000000 0.25000000000000 0.00000000000000 158.925 > 1 7 Tb 0.50000000000000 0.25000000000000 0.00000000000000 158.925 > 1 8 Tb 0.75000000000000 0.25000000000000 0.00000000000000 158.925 > 1 9 Tb 0.00000000000000 0.50000000000000 0.00000000000000 158.925 > 1 10 Tb 0.25000000000000 0.50000000000000 0.00000000000000 158.925 > 1 11 Tb 0.50000000000000 0.50000000000000 0.00000000000000 158.925 > 1 12 Tb 0.75000000000000 0.50000000000000 0.00000000000000 158.925 > 1 13 Tb 0.00000000000000 0.75000000000000 0.00000000000000 158.925 > 1 14 Tb 0.25000000000000 0.75000000000000 0.00000000000000 158.925 > 1 15 Tb 0.50000000000000 0.75000000000000 0.00000000000000 158.925 > 1 16 Tb 0.75000000000000 0.75000000000000 0.00000000000000 158.925 > 1 17 Tb 0.00000000000000 0.00000000000000 0.50000000000000 158.925 > 2 18 Tb 0.25000000000000 0.00000000000000 0.50000000000000 158.925 > 2 19 Tb 0.50000000000000 0.00000000000000 0.50000000000000 158.925 > 2 20 Tb 0.75000000000000 0.00000000000000 0.50000000000000 158.925 > 2 21 Tb 0.00000000000000 0.25000000000000 0.50000000000000 158.925 > 2 22 Tb 0.25000000000000 0.25000000000000 0.50000000000000 158.925 > 2 23 Tb 0.50000000000000 0.25000000000000 0.50000000000000 158.925 > 2 24 Tb 0.75000000000000 0.25000000000000 0.50000000000000 158.925 > 2 25 Tb 0.00000000000000 0.50000000000000 0.50000000000000 158.925 > 2 26 Tb 0.25000000000000 0.50000000000000 0.50000000000000 158.925 > 2 27 Tb 0.50000000000000 0.50000000000000 0.50000000000000 158.925 > 2 28 Tb 0.75000000000000 0.50000000000000 0.50000000000000 158.925 > 2 29 Tb 0.00000000000000 0.75000000000000 0.50000000000000 158.925 > 2 30 Tb 0.25000000000000 0.75000000000000 0.50000000000000 158.925 > 2 31 Tb 0.50000000000000 0.75000000000000 0.50000000000000 158.925 > 2 32 Tb 0.75000000000000 0.75000000000000 0.50000000000000 158.925 > 2 *33 Cs 0.08333333333333 0.16666666666667 0.75000000000000 132.905 > 3 34 Cs 0.33333333333333 0.16666666666667 0.75000000000000 132.905 > 3 35 Cs 0.58333333333333 0.16666666666667 0.75000000000000 132.905 > 3 36 Cs 0.83333333333333 0.16666666666667 0.75000000000000 132.905 > 3 37 Cs 0.08333333333333 0.41666666666667 0.75000000000000 132.905 > 3 38 Cs 0.33333333333333 0.41666666666667 0.75000000000000 132.905 > 3 39 Cs 0.58333333333333 0.41666666666667 0.75000000000000 132.905 > 3 40 Cs 0.83333333333333 0.41666666666667 0.75000000000000 132.905 > 3 41 Cs 0.08333333333333 0.66666666666667 0.75000000000000 132.905 > 3 42 Cs 0.33333333333333 0.66666666666667 0.75000000000000 132.905 > 3 43 Cs 0.58333333333333 0.66666666666667 0.75000000000000 132.905 > 3 44 Cs 0.83333333333333 0.66666666666667 0.75000000000000 132.905 > 3 45 Cs 0.08333333333333 0.91666666666667 0.75000000000000 132.905 > 3 46 Cs 0.33333333333333 0.91666666666667 0.75000000000000 132.905 > 3 47 Cs 0.58333333333333 0.91666666666667 0.75000000000000 132.905 > 3 48 Cs 0.83333333333333 0.91666666666667 0.75000000000000 132.905 > 3 49 Cs 0.16666666666667 0.08333333333333 0.25000000000000 132.905 > 4 50 Cs 0.41666666666667 0.08333333333333 0.25000000000000 132.905 > 4 51 Cs 0.66666666666667 0.08333333333333 0.25000000000000 132.905 > 4 52 Cs 0.91666666666667 0.08333333333333 0.25000000000000 132.905 > 4 53 Cs 0.16666666666667 0.33333333333333 0.25000000000000 132.905 > 4 54 Cs 0.41666666666667 0.33333333333333 0.25000000000000 132.905 > 4 55 Cs 0.66666666666667 0.33333333333333 0.25000000000000 132.905 > 4 56 Cs 0.91666666666667 0.33333333333333 0.25000000000000 132.905 > 4 57 Cs 0.16666666666667 0.58333333333333 0.25000000000000 132.905 > 4 58 Cs 0.41666666666667 0.58333333333333 0.25000000000000 132.905 > 4 59 Cs 0.66666666666667 0.58333333333333 0.25000000000000 132.905 > 4 60 Cs 0.91666666666667 0.58333333333333 0.25000000000000 132.905 > 4 61 Cs 0.16666666666667 0.83333333333333 0.25000000000000 132.905 > 4 62 Cs 0.41666666666667 0.83333333333333 0.25000000000000 132.905 > 4 63 Cs 0.66666666666667 0.83333333333333 0.25000000000000 132.905 > 4 64 Cs 0.91666666666667 0.83333333333333 0.25000000000000 132.905 > 4 *65 S 0.08333333333333 0.16666666666667 0.40770196528995 32.065 > 5 66 S 0.33333333333333 0.16666666666667 0.40770196528995 32.065 > 5 67 S 0.58333333333333 0.16666666666667 0.40770196528995 32.065 > 5 68 S 0.83333333333333 0.16666666666667 0.40770196528995 32.065 > 5 69 S 0.08333333333333 0.41666666666667 0.40770196528995 32.065 > 5 70 S 0.33333333333333 0.41666666666667 0.40770196528995 32.065 > 5 71 S 0.58333333333333 0.41666666666667 0.40770196528995 32.065 > 5 72 S 0.83333333333333 0.41666666666667 0.40770196528995 32.065 > 5 73 S 0.08333333333333 0.66666666666667 0.40770196528995 32.065 > 5 74 S 0.33333333333333 0.66666666666667 0.40770196528995 32.065 > 5 75 S 0.58333333333333 0.66666666666667 0.40770196528995 32.065 > 5 76 S 0.83333333333333 0.66666666666667 0.40770196528995 32.065 > 5 77 S 0.08333333333333 0.91666666666667 0.40770196528995 32.065 > 5 78 S 0.33333333333333 0.91666666666667 0.40770196528995 32.065 > 5 79 S 0.58333333333333 0.91666666666667 0.40770196528995 32.065 > 5 80 S 0.83333333333333 0.91666666666667 0.40770196528995 32.065 > 5 81 S 0.16666666666667 0.08333333333333 0.90770196528995 32.065 > 6 82 S 0.41666666666667 0.08333333333333 0.90770196528995 32.065 > 6 83 S 0.66666666666667 0.08333333333333 0.90770196528995 32.065 > 6 84 S 0.91666666666667 0.08333333333333 0.90770196528995 32.065 > 6 85 S 0.16666666666667 0.33333333333333 0.90770196528995 32.065 > 6 86 S 0.41666666666667 0.33333333333333 0.90770196528995 32.065 > 6 87 S 0.66666666666667 0.33333333333333 0.90770196528995 32.065 > 6 88 S 0.91666666666667 0.33333333333333 0.90770196528995 32.065 > 6 89 S 0.16666666666667 0.58333333333333 0.90770196528995 32.065 > 6 90 S 0.41666666666667 0.58333333333333 0.90770196528995 32.065 > 6 91 S 0.66666666666667 0.58333333333333 0.90770196528995 32.065 > 6 92 S 0.91666666666667 0.58333333333333 0.90770196528995 32.065 > 6 93 S 0.16666666666667 0.83333333333333 0.90770196528995 32.065 > 6 94 S 0.41666666666667 0.83333333333333 0.90770196528995 32.065 > 6 95 S 0.66666666666667 0.83333333333333 0.90770196528995 32.065 > 6 96 S 0.91666666666667 0.83333333333333 0.90770196528995 32.065 > 6 97 S 0.16666666666667 0.08333333333333 0.59229803471005 32.065 > 7 98 S 0.41666666666667 0.08333333333333 0.59229803471005 32.065 > 7 99 S 0.66666666666667 0.08333333333333 0.59229803471005 32.065 > 7 100 S 0.91666666666667 0.08333333333333 0.59229803471005 32.065 > 7 101 S 0.16666666666667 0.33333333333333 0.59229803471005 32.065 > 7 102 S 0.41666666666667 0.33333333333333 0.59229803471005 32.065 > 7 103 S 0.66666666666667 0.33333333333333 0.59229803471005 32.065 > 7 104 S 0.91666666666667 0.33333333333333 0.59229803471005 32.065 > 7 105 S 0.16666666666667 0.58333333333333 0.59229803471005 32.065 > 7 106 S 0.41666666666667 0.58333333333333 0.59229803471005 32.065 > 7 107 S 0.66666666666667 0.58333333333333 0.59229803471005 32.065 > 7 108 S 0.91666666666667 0.58333333333333 0.59229803471005 32.065 > 7 109 S 0.16666666666667 0.83333333333333 0.59229803471005 32.065 > 7 110 S 0.41666666666667 0.83333333333333 0.59229803471005 32.065 > 7 111 S 0.66666666666667 0.83333333333333 0.59229803471005 32.065 > 7 112 S 0.91666666666667 0.83333333333333 0.59229803471005 32.065 > 7 113 S 0.08333333333333 0.16666666666667 0.09229803471005 32.065 > 8 114 S 0.33333333333333 0.16666666666667 0.09229803471005 32.065 > 8 115 S 0.58333333333333 0.16666666666667 0.09229803471005 32.065 > 8 116 S 0.83333333333333 0.16666666666667 0.09229803471005 32.065 > 8 117 S 0.08333333333333 0.41666666666667 0.09229803471005 32.065 > 8 118 S 0.33333333333333 0.41666666666667 0.09229803471005 32.065 > 8 119 S 0.58333333333333 0.41666666666667 0.09229803471005 32.065 > 8 120 S 0.83333333333333 0.41666666666667 0.09229803471005 32.065 > 8 121 S 0.08333333333333 0.66666666666667 0.09229803471005 32.065 > 8 122 S 0.33333333333333 0.66666666666667 0.09229803471005 32.065 > 8 123 S 0.58333333333333 0.66666666666667 0.09229803471005 32.065 > 8 124 S 0.83333333333333 0.66666666666667 0.09229803471005 32.065 > 8 125 S 0.08333333333333 0.91666666666667 0.09229803471005 32.065 > 8 126 S 0.33333333333333 0.91666666666667 0.09229803471005 32.065 > 8 127 S 0.58333333333333 0.91666666666667 0.09229803471005 32.065 > 8 128 S 0.83333333333333 0.91666666666667 0.09229803471005 32.065 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.6752548 -0.0000000 0.0000000 0.0000000 5.6752548 0.0000000 0.0000000 0.0000000 5.0261496 -------------------------- Born effective charges -------------------------- 1 Tb 4.4777939 -0.0000000 0.0000000 0.0000000 4.4777939 0.0000000 0.0000000 0.0000000 2.1265888 2 Tb 4.4777939 -0.0000000 0.0000000 0.0000000 4.4777939 0.0000000 0.0000000 0.0000000 2.1265888 3 Cs 1.0692823 -0.0000000 0.0000000 0.0000000 1.0692823 0.0000000 0.0000000 0.0000000 2.2402605 4 Cs 1.0692823 -0.0000000 0.0000000 0.0000000 1.0692823 0.0000000 0.0000000 0.0000000 2.2402605 5 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 6 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 7 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 8 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 ---------------------------------------------------------------------------- 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: 198 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.174 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.185 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) 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: 198 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 (xx) 0.00000000 (xx) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-09 04:33:53]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:33:54]------------------------- 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.050675551774594 0.000000000000000 0.000000000000000 b -2.025337775887297 3.507987930325347 0.000000000000000 c 0.000000000000000 0.000000000000000 15.840398120000000 Atomic positions (fractional): 1 Tb 0.00000000000000 0.00000000000000 0.00000000000000 158.925 2 Tb 0.00000000000000 0.00000000000000 0.50000000000000 158.925 3 Cs 0.33333333333333 0.66666666666667 0.75000000000000 132.905 4 Cs 0.66666666666667 0.33333333333333 0.25000000000000 132.905 5 S 0.33333333333333 0.66666666666667 0.40770196528995 32.065 6 S 0.66666666666667 0.33333333333333 0.90770196528995 32.065 7 S 0.66666666666667 0.33333333333333 0.59229803471005 32.065 8 S 0.33333333333333 0.66666666666667 0.09229803471005 32.065 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.202702207098376 0.000000000000000 0.000000000000000 b -8.101351103549188 14.031951721301388 0.000000000000000 c 0.000000000000000 0.000000000000000 15.840398120000000 Atomic positions (fractional): 1 Tb 0.00000000000000 0.00000000000000 0.00000000000000 158.925 > 1 2 Tb 0.25000000000000 0.00000000000000 0.00000000000000 158.925 > 1 3 Tb 0.50000000000000 0.00000000000000 0.00000000000000 158.925 > 1 4 Tb 0.75000000000000 0.00000000000000 0.00000000000000 158.925 > 1 5 Tb 0.00000000000000 0.25000000000000 0.00000000000000 158.925 > 1 6 Tb 0.25000000000000 0.25000000000000 0.00000000000000 158.925 > 1 7 Tb 0.50000000000000 0.25000000000000 0.00000000000000 158.925 > 1 8 Tb 0.75000000000000 0.25000000000000 0.00000000000000 158.925 > 1 9 Tb 0.00000000000000 0.50000000000000 0.00000000000000 158.925 > 1 10 Tb 0.25000000000000 0.50000000000000 0.00000000000000 158.925 > 1 11 Tb 0.50000000000000 0.50000000000000 0.00000000000000 158.925 > 1 12 Tb 0.75000000000000 0.50000000000000 0.00000000000000 158.925 > 1 13 Tb 0.00000000000000 0.75000000000000 0.00000000000000 158.925 > 1 14 Tb 0.25000000000000 0.75000000000000 0.00000000000000 158.925 > 1 15 Tb 0.50000000000000 0.75000000000000 0.00000000000000 158.925 > 1 16 Tb 0.75000000000000 0.75000000000000 0.00000000000000 158.925 > 1 17 Tb 0.00000000000000 0.00000000000000 0.50000000000000 158.925 > 17 18 Tb 0.25000000000000 0.00000000000000 0.50000000000000 158.925 > 17 19 Tb 0.50000000000000 0.00000000000000 0.50000000000000 158.925 > 17 20 Tb 0.75000000000000 0.00000000000000 0.50000000000000 158.925 > 17 21 Tb 0.00000000000000 0.25000000000000 0.50000000000000 158.925 > 17 22 Tb 0.25000000000000 0.25000000000000 0.50000000000000 158.925 > 17 23 Tb 0.50000000000000 0.25000000000000 0.50000000000000 158.925 > 17 24 Tb 0.75000000000000 0.25000000000000 0.50000000000000 158.925 > 17 25 Tb 0.00000000000000 0.50000000000000 0.50000000000000 158.925 > 17 26 Tb 0.25000000000000 0.50000000000000 0.50000000000000 158.925 > 17 27 Tb 0.50000000000000 0.50000000000000 0.50000000000000 158.925 > 17 28 Tb 0.75000000000000 0.50000000000000 0.50000000000000 158.925 > 17 29 Tb 0.00000000000000 0.75000000000000 0.50000000000000 158.925 > 17 30 Tb 0.25000000000000 0.75000000000000 0.50000000000000 158.925 > 17 31 Tb 0.50000000000000 0.75000000000000 0.50000000000000 158.925 > 17 32 Tb 0.75000000000000 0.75000000000000 0.50000000000000 158.925 > 17 33 Cs 0.08333333333333 0.16666666666667 0.75000000000000 132.905 > 33 34 Cs 0.33333333333333 0.16666666666667 0.75000000000000 132.905 > 33 35 Cs 0.58333333333333 0.16666666666667 0.75000000000000 132.905 > 33 36 Cs 0.83333333333333 0.16666666666667 0.75000000000000 132.905 > 33 37 Cs 0.08333333333333 0.41666666666667 0.75000000000000 132.905 > 33 38 Cs 0.33333333333333 0.41666666666667 0.75000000000000 132.905 > 33 39 Cs 0.58333333333333 0.41666666666667 0.75000000000000 132.905 > 33 40 Cs 0.83333333333333 0.41666666666667 0.75000000000000 132.905 > 33 41 Cs 0.08333333333333 0.66666666666667 0.75000000000000 132.905 > 33 42 Cs 0.33333333333333 0.66666666666667 0.75000000000000 132.905 > 33 43 Cs 0.58333333333333 0.66666666666667 0.75000000000000 132.905 > 33 44 Cs 0.83333333333333 0.66666666666667 0.75000000000000 132.905 > 33 45 Cs 0.08333333333333 0.91666666666667 0.75000000000000 132.905 > 33 46 Cs 0.33333333333333 0.91666666666667 0.75000000000000 132.905 > 33 47 Cs 0.58333333333333 0.91666666666667 0.75000000000000 132.905 > 33 48 Cs 0.83333333333333 0.91666666666667 0.75000000000000 132.905 > 33 49 Cs 0.16666666666667 0.08333333333333 0.25000000000000 132.905 > 49 50 Cs 0.41666666666667 0.08333333333333 0.25000000000000 132.905 > 49 51 Cs 0.66666666666667 0.08333333333333 0.25000000000000 132.905 > 49 52 Cs 0.91666666666667 0.08333333333333 0.25000000000000 132.905 > 49 53 Cs 0.16666666666667 0.33333333333333 0.25000000000000 132.905 > 49 54 Cs 0.41666666666667 0.33333333333333 0.25000000000000 132.905 > 49 55 Cs 0.66666666666667 0.33333333333333 0.25000000000000 132.905 > 49 56 Cs 0.91666666666667 0.33333333333333 0.25000000000000 132.905 > 49 57 Cs 0.16666666666667 0.58333333333333 0.25000000000000 132.905 > 49 58 Cs 0.41666666666667 0.58333333333333 0.25000000000000 132.905 > 49 59 Cs 0.66666666666667 0.58333333333333 0.25000000000000 132.905 > 49 60 Cs 0.91666666666667 0.58333333333333 0.25000000000000 132.905 > 49 61 Cs 0.16666666666667 0.83333333333333 0.25000000000000 132.905 > 49 62 Cs 0.41666666666667 0.83333333333333 0.25000000000000 132.905 > 49 63 Cs 0.66666666666667 0.83333333333333 0.25000000000000 132.905 > 49 64 Cs 0.91666666666667 0.83333333333333 0.25000000000000 132.905 > 49 65 S 0.08333333333333 0.16666666666667 0.40770196528995 32.065 > 65 66 S 0.33333333333333 0.16666666666667 0.40770196528995 32.065 > 65 67 S 0.58333333333333 0.16666666666667 0.40770196528995 32.065 > 65 68 S 0.83333333333333 0.16666666666667 0.40770196528995 32.065 > 65 69 S 0.08333333333333 0.41666666666667 0.40770196528995 32.065 > 65 70 S 0.33333333333333 0.41666666666667 0.40770196528995 32.065 > 65 71 S 0.58333333333333 0.41666666666667 0.40770196528995 32.065 > 65 72 S 0.83333333333333 0.41666666666667 0.40770196528995 32.065 > 65 73 S 0.08333333333333 0.66666666666667 0.40770196528995 32.065 > 65 74 S 0.33333333333333 0.66666666666667 0.40770196528995 32.065 > 65 75 S 0.58333333333333 0.66666666666667 0.40770196528995 32.065 > 65 76 S 0.83333333333333 0.66666666666667 0.40770196528995 32.065 > 65 77 S 0.08333333333333 0.91666666666667 0.40770196528995 32.065 > 65 78 S 0.33333333333333 0.91666666666667 0.40770196528995 32.065 > 65 79 S 0.58333333333333 0.91666666666667 0.40770196528995 32.065 > 65 80 S 0.83333333333333 0.91666666666667 0.40770196528995 32.065 > 65 81 S 0.16666666666667 0.08333333333333 0.90770196528995 32.065 > 81 82 S 0.41666666666667 0.08333333333333 0.90770196528995 32.065 > 81 83 S 0.66666666666667 0.08333333333333 0.90770196528995 32.065 > 81 84 S 0.91666666666667 0.08333333333333 0.90770196528995 32.065 > 81 85 S 0.16666666666667 0.33333333333333 0.90770196528995 32.065 > 81 86 S 0.41666666666667 0.33333333333333 0.90770196528995 32.065 > 81 87 S 0.66666666666667 0.33333333333333 0.90770196528995 32.065 > 81 88 S 0.91666666666667 0.33333333333333 0.90770196528995 32.065 > 81 89 S 0.16666666666667 0.58333333333333 0.90770196528995 32.065 > 81 90 S 0.41666666666667 0.58333333333333 0.90770196528995 32.065 > 81 91 S 0.66666666666667 0.58333333333333 0.90770196528995 32.065 > 81 92 S 0.91666666666667 0.58333333333333 0.90770196528995 32.065 > 81 93 S 0.16666666666667 0.83333333333333 0.90770196528995 32.065 > 81 94 S 0.41666666666667 0.83333333333333 0.90770196528995 32.065 > 81 95 S 0.66666666666667 0.83333333333333 0.90770196528995 32.065 > 81 96 S 0.91666666666667 0.83333333333333 0.90770196528995 32.065 > 81 97 S 0.16666666666667 0.08333333333333 0.59229803471005 32.065 > 97 98 S 0.41666666666667 0.08333333333333 0.59229803471005 32.065 > 97 99 S 0.66666666666667 0.08333333333333 0.59229803471005 32.065 > 97 100 S 0.91666666666667 0.08333333333333 0.59229803471005 32.065 > 97 101 S 0.16666666666667 0.33333333333333 0.59229803471005 32.065 > 97 102 S 0.41666666666667 0.33333333333333 0.59229803471005 32.065 > 97 103 S 0.66666666666667 0.33333333333333 0.59229803471005 32.065 > 97 104 S 0.91666666666667 0.33333333333333 0.59229803471005 32.065 > 97 105 S 0.16666666666667 0.58333333333333 0.59229803471005 32.065 > 97 106 S 0.41666666666667 0.58333333333333 0.59229803471005 32.065 > 97 107 S 0.66666666666667 0.58333333333333 0.59229803471005 32.065 > 97 108 S 0.91666666666667 0.58333333333333 0.59229803471005 32.065 > 97 109 S 0.16666666666667 0.83333333333333 0.59229803471005 32.065 > 97 110 S 0.41666666666667 0.83333333333333 0.59229803471005 32.065 > 97 111 S 0.66666666666667 0.83333333333333 0.59229803471005 32.065 > 97 112 S 0.91666666666667 0.83333333333333 0.59229803471005 32.065 > 97 113 S 0.08333333333333 0.16666666666667 0.09229803471005 32.065 > 113 114 S 0.33333333333333 0.16666666666667 0.09229803471005 32.065 > 113 115 S 0.58333333333333 0.16666666666667 0.09229803471005 32.065 > 113 116 S 0.83333333333333 0.16666666666667 0.09229803471005 32.065 > 113 117 S 0.08333333333333 0.41666666666667 0.09229803471005 32.065 > 113 118 S 0.33333333333333 0.41666666666667 0.09229803471005 32.065 > 113 119 S 0.58333333333333 0.41666666666667 0.09229803471005 32.065 > 113 120 S 0.83333333333333 0.41666666666667 0.09229803471005 32.065 > 113 121 S 0.08333333333333 0.66666666666667 0.09229803471005 32.065 > 113 122 S 0.33333333333333 0.66666666666667 0.09229803471005 32.065 > 113 123 S 0.58333333333333 0.66666666666667 0.09229803471005 32.065 > 113 124 S 0.83333333333333 0.66666666666667 0.09229803471005 32.065 > 113 125 S 0.08333333333333 0.91666666666667 0.09229803471005 32.065 > 113 126 S 0.33333333333333 0.91666666666667 0.09229803471005 32.065 > 113 127 S 0.58333333333333 0.91666666666667 0.09229803471005 32.065 > 113 128 S 0.83333333333333 0.91666666666667 0.09229803471005 32.065 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.6752548 -0.0000000 0.0000000 0.0000000 5.6752548 0.0000000 0.0000000 0.0000000 5.0261496 -------------------------- Born effective charges -------------------------- 1 Tb 4.4777939 -0.0000000 0.0000000 0.0000000 4.4777939 0.0000000 0.0000000 0.0000000 2.1265888 2 Tb 4.4777939 -0.0000000 0.0000000 0.0000000 4.4777939 0.0000000 0.0000000 0.0000000 2.1265888 3 Cs 1.0692823 -0.0000000 0.0000000 0.0000000 1.0692823 0.0000000 0.0000000 0.0000000 2.2402605 4 Cs 1.0692823 -0.0000000 0.0000000 0.0000000 1.0692823 0.0000000 0.0000000 0.0000000 2.2402605 5 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 6 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 7 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 8 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 ---------------------------------------------------------------------------- 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.00000010 (xxy) 0.00000010 (xxy) 0.00000010 (xyx) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) 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-09 04:34:01]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:34:01]------------------------- 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.050675551774594 0.000000000000000 0.000000000000000 b -2.025337775887297 3.507987930325347 0.000000000000000 c 0.000000000000000 0.000000000000000 15.840398120000000 Atomic positions (fractional): 1 Tb 0.00000000000000 0.00000000000000 0.00000000000000 158.925 2 Tb 0.00000000000000 0.00000000000000 0.50000000000000 158.925 3 Cs 0.33333333333333 0.66666666666667 0.75000000000000 132.905 4 Cs 0.66666666666667 0.33333333333333 0.25000000000000 132.905 5 S 0.33333333333333 0.66666666666667 0.40770196528995 32.065 6 S 0.66666666666667 0.33333333333333 0.90770196528995 32.065 7 S 0.66666666666667 0.33333333333333 0.59229803471005 32.065 8 S 0.33333333333333 0.66666666666667 0.09229803471005 32.065 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.202702207098376 0.000000000000000 0.000000000000000 b -8.101351103549188 14.031951721301388 0.000000000000000 c 0.000000000000000 0.000000000000000 15.840398120000000 Atomic positions (fractional): 1 Tb 0.00000000000000 0.00000000000000 0.00000000000000 158.925 > 1 2 Tb 0.25000000000000 0.00000000000000 0.00000000000000 158.925 > 1 3 Tb 0.50000000000000 0.00000000000000 0.00000000000000 158.925 > 1 4 Tb 0.75000000000000 0.00000000000000 0.00000000000000 158.925 > 1 5 Tb 0.00000000000000 0.25000000000000 0.00000000000000 158.925 > 1 6 Tb 0.25000000000000 0.25000000000000 0.00000000000000 158.925 > 1 7 Tb 0.50000000000000 0.25000000000000 0.00000000000000 158.925 > 1 8 Tb 0.75000000000000 0.25000000000000 0.00000000000000 158.925 > 1 9 Tb 0.00000000000000 0.50000000000000 0.00000000000000 158.925 > 1 10 Tb 0.25000000000000 0.50000000000000 0.00000000000000 158.925 > 1 11 Tb 0.50000000000000 0.50000000000000 0.00000000000000 158.925 > 1 12 Tb 0.75000000000000 0.50000000000000 0.00000000000000 158.925 > 1 13 Tb 0.00000000000000 0.75000000000000 0.00000000000000 158.925 > 1 14 Tb 0.25000000000000 0.75000000000000 0.00000000000000 158.925 > 1 15 Tb 0.50000000000000 0.75000000000000 0.00000000000000 158.925 > 1 16 Tb 0.75000000000000 0.75000000000000 0.00000000000000 158.925 > 1 17 Tb 0.00000000000000 0.00000000000000 0.50000000000000 158.925 > 17 18 Tb 0.25000000000000 0.00000000000000 0.50000000000000 158.925 > 17 19 Tb 0.50000000000000 0.00000000000000 0.50000000000000 158.925 > 17 20 Tb 0.75000000000000 0.00000000000000 0.50000000000000 158.925 > 17 21 Tb 0.00000000000000 0.25000000000000 0.50000000000000 158.925 > 17 22 Tb 0.25000000000000 0.25000000000000 0.50000000000000 158.925 > 17 23 Tb 0.50000000000000 0.25000000000000 0.50000000000000 158.925 > 17 24 Tb 0.75000000000000 0.25000000000000 0.50000000000000 158.925 > 17 25 Tb 0.00000000000000 0.50000000000000 0.50000000000000 158.925 > 17 26 Tb 0.25000000000000 0.50000000000000 0.50000000000000 158.925 > 17 27 Tb 0.50000000000000 0.50000000000000 0.50000000000000 158.925 > 17 28 Tb 0.75000000000000 0.50000000000000 0.50000000000000 158.925 > 17 29 Tb 0.00000000000000 0.75000000000000 0.50000000000000 158.925 > 17 30 Tb 0.25000000000000 0.75000000000000 0.50000000000000 158.925 > 17 31 Tb 0.50000000000000 0.75000000000000 0.50000000000000 158.925 > 17 32 Tb 0.75000000000000 0.75000000000000 0.50000000000000 158.925 > 17 33 Cs 0.08333333333333 0.16666666666667 0.75000000000000 132.905 > 33 34 Cs 0.33333333333333 0.16666666666667 0.75000000000000 132.905 > 33 35 Cs 0.58333333333333 0.16666666666667 0.75000000000000 132.905 > 33 36 Cs 0.83333333333333 0.16666666666667 0.75000000000000 132.905 > 33 37 Cs 0.08333333333333 0.41666666666667 0.75000000000000 132.905 > 33 38 Cs 0.33333333333333 0.41666666666667 0.75000000000000 132.905 > 33 39 Cs 0.58333333333333 0.41666666666667 0.75000000000000 132.905 > 33 40 Cs 0.83333333333333 0.41666666666667 0.75000000000000 132.905 > 33 41 Cs 0.08333333333333 0.66666666666667 0.75000000000000 132.905 > 33 42 Cs 0.33333333333333 0.66666666666667 0.75000000000000 132.905 > 33 43 Cs 0.58333333333333 0.66666666666667 0.75000000000000 132.905 > 33 44 Cs 0.83333333333333 0.66666666666667 0.75000000000000 132.905 > 33 45 Cs 0.08333333333333 0.91666666666667 0.75000000000000 132.905 > 33 46 Cs 0.33333333333333 0.91666666666667 0.75000000000000 132.905 > 33 47 Cs 0.58333333333333 0.91666666666667 0.75000000000000 132.905 > 33 48 Cs 0.83333333333333 0.91666666666667 0.75000000000000 132.905 > 33 49 Cs 0.16666666666667 0.08333333333333 0.25000000000000 132.905 > 49 50 Cs 0.41666666666667 0.08333333333333 0.25000000000000 132.905 > 49 51 Cs 0.66666666666667 0.08333333333333 0.25000000000000 132.905 > 49 52 Cs 0.91666666666667 0.08333333333333 0.25000000000000 132.905 > 49 53 Cs 0.16666666666667 0.33333333333333 0.25000000000000 132.905 > 49 54 Cs 0.41666666666667 0.33333333333333 0.25000000000000 132.905 > 49 55 Cs 0.66666666666667 0.33333333333333 0.25000000000000 132.905 > 49 56 Cs 0.91666666666667 0.33333333333333 0.25000000000000 132.905 > 49 57 Cs 0.16666666666667 0.58333333333333 0.25000000000000 132.905 > 49 58 Cs 0.41666666666667 0.58333333333333 0.25000000000000 132.905 > 49 59 Cs 0.66666666666667 0.58333333333333 0.25000000000000 132.905 > 49 60 Cs 0.91666666666667 0.58333333333333 0.25000000000000 132.905 > 49 61 Cs 0.16666666666667 0.83333333333333 0.25000000000000 132.905 > 49 62 Cs 0.41666666666667 0.83333333333333 0.25000000000000 132.905 > 49 63 Cs 0.66666666666667 0.83333333333333 0.25000000000000 132.905 > 49 64 Cs 0.91666666666667 0.83333333333333 0.25000000000000 132.905 > 49 65 S 0.08333333333333 0.16666666666667 0.40770196528995 32.065 > 65 66 S 0.33333333333333 0.16666666666667 0.40770196528995 32.065 > 65 67 S 0.58333333333333 0.16666666666667 0.40770196528995 32.065 > 65 68 S 0.83333333333333 0.16666666666667 0.40770196528995 32.065 > 65 69 S 0.08333333333333 0.41666666666667 0.40770196528995 32.065 > 65 70 S 0.33333333333333 0.41666666666667 0.40770196528995 32.065 > 65 71 S 0.58333333333333 0.41666666666667 0.40770196528995 32.065 > 65 72 S 0.83333333333333 0.41666666666667 0.40770196528995 32.065 > 65 73 S 0.08333333333333 0.66666666666667 0.40770196528995 32.065 > 65 74 S 0.33333333333333 0.66666666666667 0.40770196528995 32.065 > 65 75 S 0.58333333333333 0.66666666666667 0.40770196528995 32.065 > 65 76 S 0.83333333333333 0.66666666666667 0.40770196528995 32.065 > 65 77 S 0.08333333333333 0.91666666666667 0.40770196528995 32.065 > 65 78 S 0.33333333333333 0.91666666666667 0.40770196528995 32.065 > 65 79 S 0.58333333333333 0.91666666666667 0.40770196528995 32.065 > 65 80 S 0.83333333333333 0.91666666666667 0.40770196528995 32.065 > 65 81 S 0.16666666666667 0.08333333333333 0.90770196528995 32.065 > 81 82 S 0.41666666666667 0.08333333333333 0.90770196528995 32.065 > 81 83 S 0.66666666666667 0.08333333333333 0.90770196528995 32.065 > 81 84 S 0.91666666666667 0.08333333333333 0.90770196528995 32.065 > 81 85 S 0.16666666666667 0.33333333333333 0.90770196528995 32.065 > 81 86 S 0.41666666666667 0.33333333333333 0.90770196528995 32.065 > 81 87 S 0.66666666666667 0.33333333333333 0.90770196528995 32.065 > 81 88 S 0.91666666666667 0.33333333333333 0.90770196528995 32.065 > 81 89 S 0.16666666666667 0.58333333333333 0.90770196528995 32.065 > 81 90 S 0.41666666666667 0.58333333333333 0.90770196528995 32.065 > 81 91 S 0.66666666666667 0.58333333333333 0.90770196528995 32.065 > 81 92 S 0.91666666666667 0.58333333333333 0.90770196528995 32.065 > 81 93 S 0.16666666666667 0.83333333333333 0.90770196528995 32.065 > 81 94 S 0.41666666666667 0.83333333333333 0.90770196528995 32.065 > 81 95 S 0.66666666666667 0.83333333333333 0.90770196528995 32.065 > 81 96 S 0.91666666666667 0.83333333333333 0.90770196528995 32.065 > 81 97 S 0.16666666666667 0.08333333333333 0.59229803471005 32.065 > 97 98 S 0.41666666666667 0.08333333333333 0.59229803471005 32.065 > 97 99 S 0.66666666666667 0.08333333333333 0.59229803471005 32.065 > 97 100 S 0.91666666666667 0.08333333333333 0.59229803471005 32.065 > 97 101 S 0.16666666666667 0.33333333333333 0.59229803471005 32.065 > 97 102 S 0.41666666666667 0.33333333333333 0.59229803471005 32.065 > 97 103 S 0.66666666666667 0.33333333333333 0.59229803471005 32.065 > 97 104 S 0.91666666666667 0.33333333333333 0.59229803471005 32.065 > 97 105 S 0.16666666666667 0.58333333333333 0.59229803471005 32.065 > 97 106 S 0.41666666666667 0.58333333333333 0.59229803471005 32.065 > 97 107 S 0.66666666666667 0.58333333333333 0.59229803471005 32.065 > 97 108 S 0.91666666666667 0.58333333333333 0.59229803471005 32.065 > 97 109 S 0.16666666666667 0.83333333333333 0.59229803471005 32.065 > 97 110 S 0.41666666666667 0.83333333333333 0.59229803471005 32.065 > 97 111 S 0.66666666666667 0.83333333333333 0.59229803471005 32.065 > 97 112 S 0.91666666666667 0.83333333333333 0.59229803471005 32.065 > 97 113 S 0.08333333333333 0.16666666666667 0.09229803471005 32.065 > 113 114 S 0.33333333333333 0.16666666666667 0.09229803471005 32.065 > 113 115 S 0.58333333333333 0.16666666666667 0.09229803471005 32.065 > 113 116 S 0.83333333333333 0.16666666666667 0.09229803471005 32.065 > 113 117 S 0.08333333333333 0.41666666666667 0.09229803471005 32.065 > 113 118 S 0.33333333333333 0.41666666666667 0.09229803471005 32.065 > 113 119 S 0.58333333333333 0.41666666666667 0.09229803471005 32.065 > 113 120 S 0.83333333333333 0.41666666666667 0.09229803471005 32.065 > 113 121 S 0.08333333333333 0.66666666666667 0.09229803471005 32.065 > 113 122 S 0.33333333333333 0.66666666666667 0.09229803471005 32.065 > 113 123 S 0.58333333333333 0.66666666666667 0.09229803471005 32.065 > 113 124 S 0.83333333333333 0.66666666666667 0.09229803471005 32.065 > 113 125 S 0.08333333333333 0.91666666666667 0.09229803471005 32.065 > 113 126 S 0.33333333333333 0.91666666666667 0.09229803471005 32.065 > 113 127 S 0.58333333333333 0.91666666666667 0.09229803471005 32.065 > 113 128 S 0.83333333333333 0.91666666666667 0.09229803471005 32.065 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 5.6752548 -0.0000000 0.0000000 0.0000000 5.6752548 0.0000000 0.0000000 0.0000000 5.0261496 -------------------------- Born effective charges -------------------------- 1 Tb 4.4777939 -0.0000000 0.0000000 0.0000000 4.4777939 0.0000000 0.0000000 0.0000000 2.1265888 2 Tb 4.4777939 -0.0000000 0.0000000 0.0000000 4.4777939 0.0000000 0.0000000 0.0000000 2.1265888 3 Cs 1.0692823 -0.0000000 0.0000000 0.0000000 1.0692823 0.0000000 0.0000000 0.0000000 2.2402605 4 Cs 1.0692823 -0.0000000 0.0000000 0.0000000 1.0692823 0.0000000 0.0000000 0.0000000 2.2402605 5 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 6 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 7 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 8 S -2.7735381 0.0000000 0.0000000 0.0000000 -2.7735381 0.0000000 0.0000000 0.0000000 -2.1834247 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000010 (xxy) 0.00000010 (xxy) 0.00000010 (xyx) Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 14 14 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.68, Number of G-points: 291, Lambda: 0.17 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/48) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 48 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 1.108 ( 0.000 0.000 0.000) 0.000 1.108 ( 0.000 0.000 0.000) 0.000 1.505 ( 0.000 0.000 0.000) 0.000 1.505 ( 0.000 0.000 0.000) 0.000 1.850 ( 0.000 0.000 0.000) 0.000 1.850 ( 0.000 0.000 0.000) 0.000 2.202 ( 0.000 0.000 0.000) 0.000 2.691 ( 0.000 0.000 0.000) 0.000 3.090 ( 0.000 0.000 0.000) 0.000 4.904 ( 0.000 0.000 0.000) 0.000 4.904 ( 0.000 0.000 0.000) 0.000 4.936 ( 0.000 0.000 0.000) 0.000 4.936 ( 0.000 0.000 0.000) 0.000 6.309 ( 0.000 0.000 0.000) 0.000 6.309 ( 0.000 0.000 0.000) 0.000 6.329 ( 0.000 0.000 0.000) 0.000 6.329 ( 0.000 0.000 0.000) 0.000 6.744 ( 0.000 0.000 0.000) 0.000 8.046 ( 0.000 0.000 0.000) 0.000 8.167 ( 0.000 0.000 0.000) 0.000 8.367 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/48) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.446 ( 18.772 10.838 0.000) 21.676 0.485 ( 19.531 11.276 0.000) 22.552 0.814 ( 32.969 19.035 0.000) 38.069 1.225 ( 9.092 5.250 0.000) 10.499 1.327 ( 15.734 9.084 0.000) 18.168 1.512 ( -0.139 -0.080 0.000) 0.161 1.516 ( 0.892 0.515 0.000) 1.030 1.868 ( 1.496 0.864 0.000) 1.728 1.931 ( 6.865 3.963 0.000) 7.927 2.226 ( 2.994 1.728 0.000) 3.457 2.627 ( -2.655 -1.533 0.000) 3.065 3.068 ( -1.880 -1.085 0.000) 2.171 4.957 ( 4.321 2.495 0.000) 4.990 4.983 ( 3.948 2.279 0.000) 4.559 5.384 ( 33.339 19.248 0.000) 38.497 6.253 ( -4.460 -2.575 0.000) 5.149 6.281 ( -3.768 -2.175 0.000) 4.351 6.288 ( -1.673 -0.966 0.000) 1.932 6.308 ( -1.713 -0.989 0.000) 1.978 6.753 ( 0.677 0.391 0.000) 0.782 7.933 ( -7.702 -4.447 0.000) 8.893 8.128 ( -3.018 -1.742 0.000) 3.484 8.210 ( -6.863 -3.962 0.000) 7.925 8.244 ( -2.139 -1.235 0.000) 2.470 ======================= Grid point 2 (3/48) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.877 ( 17.733 10.238 0.000) 20.477 0.898 ( 15.331 8.851 0.000) 17.702 1.466 ( 7.493 4.326 0.000) 8.652 1.483 ( 11.857 6.846 0.000) 13.692 1.519 ( 15.039 8.683 0.000) 17.366 1.546 ( 1.635 0.944 0.000) 1.888 1.695 ( 12.950 7.477 0.000) 14.953 1.921 ( 2.968 1.714 0.000) 3.428 2.185 ( 15.352 8.864 0.000) 17.727 2.370 ( 8.792 5.076 0.000) 10.152 2.664 ( 6.739 3.890 0.000) 7.781 3.003 ( -3.553 -2.051 0.000) 4.103 5.094 ( 6.929 4.000 0.000) 8.001 5.109 ( 6.393 3.691 0.000) 7.382 6.102 ( 2.605 1.504 0.000) 3.008 6.146 ( 16.311 9.417 0.000) 18.834 6.170 ( -5.313 -3.068 0.000) 6.135 6.236 ( -2.534 -1.463 0.000) 2.926 6.254 ( -2.655 -1.533 0.000) 3.065 6.775 ( 1.299 0.750 0.000) 1.500 7.713 ( -10.980 -6.339 0.000) 12.679 7.875 ( -16.571 -9.568 0.000) 19.135 8.056 ( -2.240 -1.293 0.000) 2.587 8.274 ( 0.838 0.484 0.000) 0.968 ======================= Grid point 3 (4/48) ======================= q-point: ( 0.21 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.202 ( 10.560 6.097 0.000) 12.193 1.268 ( 15.193 8.772 0.000) 17.544 1.488 ( 2.354 1.359 0.000) 2.718 1.589 ( 1.938 1.119 0.000) 2.238 1.754 ( 10.880 6.282 0.000) 12.563 1.892 ( 4.501 2.599 0.000) 5.198 1.910 ( 11.686 6.747 0.000) 13.494 2.007 ( 4.346 2.509 0.000) 5.019 2.589 ( 8.938 5.160 0.000) 10.320 2.649 ( 22.520 13.002 0.000) 26.004 2.854 ( -0.942 -0.544 0.000) 1.088 2.983 ( 10.812 6.242 0.000) 12.484 5.258 ( 6.502 3.754 0.000) 7.507 5.262 ( 6.124 3.536 0.000) 7.071 5.954 ( -6.180 -3.568 0.000) 7.136 6.050 ( -4.462 -2.576 0.000) 5.153 6.181 ( -1.906 -1.101 0.000) 2.201 6.194 ( -2.212 -1.277 0.000) 2.554 6.508 ( 6.771 3.909 0.000) 7.818 6.819 ( 2.129 1.229 0.000) 2.459 7.405 ( -15.122 -8.731 0.000) 17.462 7.507 ( -13.823 -7.980 0.000) 15.961 8.047 ( 1.303 0.752 0.000) 1.505 8.275 ( -0.746 -0.431 0.000) 0.862 ======================= Grid point 4 (5/48) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.404 ( 6.910 3.989 0.000) 7.979 1.580 ( 11.195 6.464 0.000) 12.927 1.597 ( 6.633 3.830 0.000) 7.659 1.632 ( 1.648 0.951 0.000) 1.903 1.952 ( 1.322 0.763 0.000) 1.527 1.987 ( 8.896 5.136 0.000) 10.273 2.107 ( 5.769 3.331 0.000) 6.661 2.120 ( 5.067 2.925 0.000) 5.850 2.728 ( -6.760 -3.903 0.000) 7.806 2.746 ( 3.907 2.256 0.000) 4.512 3.180 ( 21.665 12.508 0.000) 25.016 3.342 ( 15.969 9.219 0.000) 18.439 5.383 ( 4.083 2.357 0.000) 4.714 5.385 ( 4.213 2.432 0.000) 4.864 5.847 ( -2.520 -1.455 0.000) 2.910 5.978 ( -1.334 -0.770 0.000) 1.540 6.152 ( -0.674 -0.389 0.000) 0.778 6.155 ( -1.140 -0.658 0.000) 1.317 6.551 ( -1.335 -0.771 0.000) 1.542 6.845 ( -0.020 -0.012 0.000) 0.023 7.060 ( -13.344 -7.704 0.000) 15.409 7.206 ( -11.958 -6.904 0.000) 13.808 8.082 ( 0.814 0.470 0.000) 0.940 8.238 ( -2.362 -1.364 0.000) 2.728 ======================= Grid point 5 (6/48) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.536 ( 4.363 2.519 0.000) 5.038 1.665 ( 1.151 0.664 0.000) 1.329 1.776 ( 7.982 4.609 0.000) 9.217 1.791 ( 6.715 3.877 0.000) 7.754 1.973 ( 0.699 0.403 0.000) 0.807 2.170 ( 6.544 3.778 0.000) 7.556 2.204 ( 2.892 1.670 0.000) 3.340 2.236 ( 4.563 2.635 0.000) 5.269 2.564 ( -6.819 -3.937 0.000) 7.874 2.765 ( -2.152 -1.243 0.000) 2.485 3.640 ( 16.972 9.799 0.000) 19.598 3.688 ( 12.952 7.478 0.000) 14.956 5.460 ( 2.642 1.525 0.000) 3.051 5.464 ( 2.778 1.604 0.000) 3.208 5.845 ( 2.226 1.285 0.000) 2.570 6.000 ( 3.312 1.912 0.000) 3.824 6.132 ( -1.020 -0.589 0.000) 1.178 6.139 ( -0.674 -0.389 0.000) 0.778 6.489 ( -3.477 -2.008 0.000) 4.015 6.780 ( -10.807 -6.239 0.000) 12.478 6.833 ( -0.733 -0.423 0.000) 0.846 6.959 ( -8.590 -4.960 0.000) 9.919 8.065 ( -2.158 -1.246 0.000) 2.492 8.167 ( -3.551 -2.050 0.000) 4.101 ======================= Grid point 6 (7/48) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.610 ( 2.007 1.159 0.000) 2.317 1.688 ( 0.772 0.446 0.000) 0.891 1.896 ( 2.302 1.329 0.000) 2.658 1.940 ( 5.252 3.033 0.000) 6.065 1.995 ( 1.217 0.703 0.000) 1.406 2.264 ( 2.670 1.542 0.000) 3.083 2.290 ( 3.579 2.066 0.000) 4.132 2.324 ( 2.761 1.594 0.000) 3.188 2.426 ( -4.562 -2.634 0.000) 5.267 2.668 ( -5.182 -2.992 0.000) 5.984 3.932 ( 7.428 4.289 0.000) 8.577 3.956 ( 9.509 5.490 0.000) 10.980 5.512 ( 1.731 1.000 0.000) 1.999 5.520 ( 1.919 1.108 0.000) 2.216 5.932 ( 4.441 2.564 0.000) 5.128 6.105 ( -1.128 -0.651 0.000) 1.302 6.117 ( -1.020 -0.589 0.000) 1.178 6.127 ( 7.205 4.160 0.000) 8.320 6.402 ( -3.635 -2.099 0.000) 4.197 6.540 ( -9.745 -5.626 0.000) 11.252 6.812 ( -3.995 -2.307 0.000) 4.613 6.816 ( -0.563 -0.325 0.000) 0.650 8.000 ( -2.742 -1.583 0.000) 3.166 8.086 ( -2.918 -1.685 0.000) 3.369 ======================= Grid point 7 (8/48) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 64 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.633 ( 0.000 0.000 0.000) 0.000 1.698 ( 0.000 0.000 0.000) 0.000 1.918 ( 0.000 0.000 0.000) 0.000 2.002 ( 0.000 0.000 0.000) 0.000 2.014 ( 0.000 0.000 0.000) 0.000 2.312 ( 0.000 0.000 0.000) 0.000 2.333 ( 0.000 0.000 0.000) 0.000 2.358 ( 0.000 0.000 0.000) 0.000 2.369 ( 0.000 0.000 0.000) 0.000 2.586 ( 0.000 0.000 0.000) 0.000 4.022 ( 0.000 0.000 0.000) 0.000 4.071 ( 0.000 0.000 0.000) 0.000 5.534 ( 0.000 0.000 0.000) 0.000 5.545 ( 0.000 0.000 0.000) 0.000 5.999 ( 0.000 0.000 0.000) 0.000 6.090 ( 0.000 0.000 0.000) 0.000 6.103 ( 0.000 0.000 0.000) 0.000 6.319 ( 0.000 0.000 0.000) 0.001 6.322 ( -0.000 -0.000 0.000) 0.000 6.348 ( 0.000 0.000 0.000) 0.000 6.766 ( 0.000 0.000 0.000) 0.000 6.809 ( 0.000 0.000 0.000) 0.000 7.964 ( 0.000 0.000 0.000) 0.000 8.048 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 16 (9/48) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.762 ( 10.384 17.986 0.000) 20.768 0.827 ( 10.691 18.518 0.000) 21.383 1.309 ( 13.767 23.845 0.000) 27.534 1.424 ( 7.429 12.868 0.000) 14.858 1.482 ( -1.343 -2.326 0.000) 2.686 1.563 ( 2.728 4.725 0.000) 5.456 1.596 ( 8.396 14.542 0.000) 16.791 1.914 ( 2.020 3.499 0.000) 4.041 2.088 ( 6.900 11.951 0.000) 13.800 2.313 ( 4.154 7.195 0.000) 8.308 2.629 ( 2.044 3.541 0.000) 4.089 3.024 ( -1.854 -3.212 0.000) 3.708 5.061 ( 4.264 7.386 0.000) 8.529 5.077 ( 3.860 6.686 0.000) 7.720 5.952 ( 17.353 30.056 0.000) 34.706 6.155 ( -3.830 -6.634 0.000) 7.661 6.201 ( -3.088 -5.349 0.000) 6.177 6.255 ( -1.220 -2.112 0.000) 2.439 6.272 ( -1.377 -2.385 0.000) 2.754 6.765 ( 0.435 0.753 0.000) 0.870 7.781 ( -5.712 -9.893 0.000) 11.424 7.988 ( -8.971 -15.538 0.000) 17.942 8.071 ( -1.819 -3.150 0.000) 3.637 8.252 ( 0.230 0.399 0.000) 0.460 ======================= Grid point 17 (10/48) ======================= q-point: ( 0.14 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.133 ( 12.460 14.035 0.000) 18.768 1.164 ( 7.402 14.884 0.000) 16.623 1.450 ( 0.791 -1.142 0.000) 1.389 1.566 ( 4.014 2.700 0.000) 4.838 1.698 ( 7.632 13.156 0.000) 15.210 1.807 ( 9.667 12.389 0.000) 15.714 1.833 ( 6.047 5.402 0.000) 8.109 2.001 ( 2.387 5.945 0.000) 6.407 2.442 ( 16.522 16.413 0.000) 23.289 2.500 ( 7.560 7.281 0.000) 10.496 2.806 ( 8.929 10.035 0.000) 13.433 2.943 ( -2.414 -3.708 0.000) 4.425 5.236 ( 5.036 9.529 0.000) 10.778 5.237 ( 4.796 8.616 0.000) 9.861 6.004 ( -5.105 -6.347 0.000) 8.145 6.082 ( -3.793 -5.223 0.000) 6.455 6.205 ( -1.815 -1.434 0.000) 2.313 6.219 ( -1.928 -1.728 0.000) 2.589 6.434 ( 10.054 12.309 0.000) 15.893 6.786 ( 1.652 0.250 0.000) 1.671 7.540 ( -10.116 -10.217 0.000) 14.377 7.638 ( -12.525 -13.157 0.000) 18.166 8.028 ( 0.035 -1.301 0.000) 1.301 8.246 ( 0.567 -3.110 0.000) 3.161 ======================= Grid point 18 (11/48) ======================= q-point: ( 0.21 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.383 ( 4.446 8.721 0.000) 9.789 1.467 ( 10.892 10.384 0.000) 15.049 1.506 ( 3.796 3.275 0.000) 5.013 1.656 ( 2.550 4.745 0.000) 5.387 1.913 ( 3.130 3.668 0.000) 4.822 1.968 ( 5.167 9.763 0.000) 11.046 2.027 ( 5.974 4.741 0.000) 7.627 2.121 ( 2.984 8.072 0.000) 8.605 2.684 ( 5.948 3.818 0.000) 7.068 2.802 ( -5.793 -3.562 0.000) 6.800 2.946 ( 20.214 15.659 0.000) 25.570 3.165 ( 14.904 10.573 0.000) 18.273 5.399 ( 2.462 9.675 0.000) 9.984 5.406 ( 2.262 10.428 0.000) 10.670 5.870 ( -3.307 -4.252 0.000) 5.387 5.978 ( -2.024 -4.276 0.000) 4.731 6.167 ( -1.299 -0.353 0.000) 1.346 6.179 ( -1.408 -0.083 0.000) 1.411 6.572 ( -0.378 2.171 0.000) 2.204 6.823 ( 1.981 -0.747 0.000) 2.117 7.220 ( -13.632 -9.079 0.000) 16.378 7.320 ( -10.483 -10.019 0.000) 14.501 8.040 ( 2.268 -1.769 0.000) 2.876 8.211 ( 0.596 -5.555 0.000) 5.586 ======================= Grid point 19 (12/48) ======================= q-point: ( 0.29 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.530 ( 3.911 6.025 0.000) 7.182 1.607 ( 3.288 2.946 0.000) 4.415 1.718 ( 7.488 6.849 0.000) 10.148 1.763 ( 2.725 7.580 0.000) 8.055 1.955 ( 1.184 0.050 0.000) 1.185 2.140 ( 3.890 1.284 0.000) 4.097 2.165 ( 3.762 11.055 0.000) 11.677 2.257 ( 2.141 9.064 0.000) 9.314 2.652 ( -6.973 -3.092 0.000) 7.628 2.765 ( 0.979 -0.267 0.000) 1.014 3.434 ( 18.396 12.041 0.000) 21.987 3.528 ( 14.191 8.810 0.000) 16.703 5.514 ( -0.702 10.165 0.000) 10.189 5.522 ( -1.059 10.458 0.000) 10.512 5.818 ( 0.766 -1.077 0.000) 1.322 5.939 ( 2.042 -2.815 0.000) 3.478 6.149 ( -0.782 -0.122 0.000) 0.792 6.166 ( -1.219 1.766 0.000) 2.146 6.536 ( -3.496 -0.620 0.000) 3.550 6.831 ( 0.238 -1.153 0.000) 1.178 6.911 ( -11.813 -7.352 0.000) 13.914 7.056 ( -9.140 -7.115 0.000) 11.583 8.044 ( 1.153 -4.118 0.000) 4.277 8.155 ( -0.071 -6.397 0.000) 6.397 ======================= Grid point 20 (13/48) ======================= q-point: ( 0.36 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.635 ( 0.801 6.074 0.000) 6.127 1.674 ( 1.631 1.781 0.000) 2.415 1.870 ( 3.576 3.720 0.000) 5.160 1.902 ( 4.141 6.459 0.000) 7.672 1.972 ( 1.153 -0.215 0.000) 1.173 2.201 ( 3.014 -0.573 0.000) 3.068 2.314 ( 1.454 8.557 0.000) 8.679 2.376 ( 0.423 9.627 0.000) 9.637 2.510 ( -6.791 -0.018 0.000) 6.791 2.718 ( -3.927 -2.799 0.000) 4.822 3.817 ( 13.296 7.263 0.000) 15.151 3.823 ( 10.423 5.503 0.000) 11.787 5.589 ( -2.614 10.483 0.000) 10.803 5.592 ( -2.625 10.156 0.000) 10.490 5.870 ( 4.411 1.276 0.000) 4.591 5.993 ( 6.182 -1.852 0.000) 6.453 6.128 ( -1.036 -0.786 0.000) 1.300 6.167 ( -1.368 3.391 0.000) 3.656 6.456 ( -4.341 -0.911 0.000) 4.436 6.647 ( -9.918 -6.943 0.000) 12.107 6.818 ( -0.271 -0.982 0.000) 1.018 6.869 ( -5.805 -3.615 0.000) 6.839 7.999 ( -0.614 -4.800 0.000) 4.839 8.084 ( -0.751 -5.519 0.000) 5.569 ======================= Grid point 21 (14/48) ======================= q-point: ( 0.43 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.681 ( -1.893 5.256 0.000) 5.586 1.706 ( 0.080 1.628 0.000) 1.630 1.923 ( 0.178 1.366 0.000) 1.378 1.995 ( 1.417 -0.438 0.000) 1.483 2.006 ( 1.368 3.216 0.000) 3.495 2.245 ( 2.979 -1.651 0.000) 3.406 2.371 ( -0.430 0.668 0.000) 0.794 2.443 ( -5.606 9.013 0.000) 10.614 2.466 ( -1.674 10.199 0.000) 10.335 2.603 ( -4.476 -2.984 0.000) 5.379 3.990 ( 4.289 1.173 0.000) 4.447 4.030 ( 5.369 1.596 0.000) 5.601 5.634 ( -4.367 10.102 0.000) 11.005 5.635 ( -3.957 9.276 0.000) 10.084 5.964 ( 3.750 0.387 0.000) 3.770 6.064 ( 2.702 -2.005 0.000) 3.364 6.104 ( -0.498 -0.702 0.000) 0.861 6.231 ( 5.266 4.391 0.000) 6.856 6.379 ( -3.294 0.254 0.000) 3.303 6.420 ( -8.615 -5.527 0.000) 10.235 6.777 ( -1.934 -1.142 0.000) 2.246 6.806 ( -0.001 -0.682 0.000) 0.682 7.947 ( 0.057 -3.318 0.000) 3.319 8.030 ( 0.228 -3.509 0.000) 3.516 ======================= Grid point 31 (15/48) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.379 ( 2.724 4.718 0.000) 5.448 1.422 ( 7.935 13.744 0.000) 15.870 1.528 ( 5.604 9.706 0.000) 11.208 1.634 ( 2.345 4.062 0.000) 4.690 1.909 ( 1.414 2.449 0.000) 2.828 1.975 ( 7.541 13.061 0.000) 15.082 2.001 ( 4.039 6.996 0.000) 8.079 2.145 ( 4.584 7.941 0.000) 9.169 2.645 ( 3.692 6.395 0.000) 7.384 2.835 ( 12.260 21.235 0.000) 24.520 2.843 ( -3.000 -5.196 0.000) 6.000 3.074 ( 8.631 14.950 0.000) 17.263 5.437 ( 5.889 10.199 0.000) 11.777 5.452 ( 6.232 10.794 0.000) 12.463 5.877 ( -3.190 -5.525 0.000) 6.380 5.973 ( -2.943 -5.097 0.000) 5.886 6.179 ( -0.609 -1.054 0.000) 1.218 6.192 ( -0.409 -0.708 0.000) 0.818 6.596 ( 2.001 3.466 0.000) 4.002 6.789 ( 0.160 0.277 0.000) 0.320 7.305 ( -7.482 -12.959 0.000) 14.963 7.365 ( -7.404 -12.824 0.000) 14.808 8.004 ( -0.683 -1.183 0.000) 1.366 8.162 ( -2.680 -4.642 0.000) 5.360 ======================= Grid point 32 (16/48) ======================= q-point: ( 0.21 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.476 ( 4.662 3.611 0.000) 5.897 1.673 ( 6.822 9.448 0.000) 11.653 1.678 ( -0.630 9.366 0.000) 9.387 1.763 ( 5.918 5.992 0.000) 8.422 1.942 ( 0.775 0.781 0.000) 1.100 2.104 ( 2.450 2.726 0.000) 3.665 2.215 ( 4.869 11.830 0.000) 12.793 2.315 ( 3.573 10.545 0.000) 11.134 2.733 ( -4.597 -2.555 0.000) 5.260 2.742 ( 2.124 2.107 0.000) 2.992 3.270 ( 14.863 15.845 0.000) 21.725 3.394 ( 11.532 11.587 0.000) 16.348 5.625 ( 2.447 11.403 0.000) 11.663 5.644 ( 2.026 11.554 0.000) 11.731 5.800 ( -0.643 -2.119 0.000) 2.214 5.891 ( -0.716 -3.878 0.000) 3.944 6.161 ( -0.747 -0.422 0.000) 0.858 6.196 ( 0.139 1.868 0.000) 1.873 6.599 ( -2.821 0.726 0.000) 2.913 6.800 ( 1.325 -1.245 0.000) 1.819 7.021 ( -10.104 -11.026 0.000) 14.956 7.131 ( -7.445 -7.612 0.000) 10.648 7.974 ( 0.453 -4.330 0.000) 4.354 8.072 ( -0.321 -7.411 0.000) 7.418 ======================= Grid point 33 (17/48) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 200 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.608 ( 5.228 4.034 0.000) 6.604 1.744 ( -1.294 7.468 0.000) 7.579 1.847 ( 4.200 5.520 0.000) 6.936 1.911 ( 3.368 6.717 0.000) 7.514 1.958 ( 0.525 0.276 0.000) 0.593 2.155 ( 1.918 0.161 0.000) 1.924 2.367 ( 1.330 5.867 0.000) 6.016 2.467 ( 1.357 9.692 0.000) 9.786 2.656 ( -5.061 6.090 0.000) 7.918 2.761 ( -1.810 1.672 0.000) 2.464 3.656 ( 13.632 9.419 0.000) 16.570 3.691 ( 10.865 6.831 0.000) 12.834 5.749 ( -1.134 11.026 0.000) 11.084 5.751 ( -1.151 8.852 0.000) 8.927 5.823 ( 2.126 3.318 0.000) 3.941 5.884 ( 3.465 -1.830 0.000) 3.919 6.137 ( -1.102 -1.209 0.000) 1.636 6.223 ( -0.673 3.644 0.000) 3.706 6.541 ( -4.706 1.700 0.000) 5.004 6.745 ( -9.516 -8.593 0.000) 12.822 6.803 ( 0.613 -1.539 0.000) 1.657 6.946 ( -7.040 -3.562 0.000) 7.890 7.933 ( 0.992 -6.040 0.000) 6.121 8.004 ( 1.363 -7.644 0.000) 7.765 ======================= Grid point 34 (18/48) ======================= q-point: ( 0.36 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.717 ( 1.963 4.085 0.000) 4.532 1.785 ( -2.159 6.484 0.000) 6.834 1.934 ( 1.036 2.438 0.000) 2.649 1.970 ( 0.930 0.057 0.000) 0.932 2.019 ( 1.934 4.884 0.000) 5.253 2.186 ( 2.248 -0.942 0.000) 2.437 2.385 ( -0.951 -0.367 0.000) 1.020 2.534 ( -0.110 0.663 0.000) 0.672 2.680 ( -5.943 15.122 0.000) 16.248 2.744 ( -6.153 10.451 0.000) 12.127 3.898 ( 7.381 1.428 0.000) 7.518 3.922 ( 8.996 2.665 0.000) 9.382 5.810 ( -2.326 8.868 0.000) 9.168 5.820 ( -2.633 9.778 0.000) 10.126 5.905 ( 2.841 2.687 0.000) 3.910 5.961 ( 5.560 -0.699 0.000) 5.604 6.103 ( -0.722 -1.597 0.000) 1.753 6.238 ( -0.947 2.341 0.000) 2.526 6.474 ( -5.532 3.285 0.000) 6.433 6.525 ( -8.350 -2.893 0.000) 8.837 6.795 ( 0.235 -1.222 0.000) 1.244 6.816 ( -4.435 -1.594 0.000) 4.713 7.893 ( 1.074 -5.068 0.000) 5.180 7.963 ( 1.813 -5.944 0.000) 6.215 ======================= Grid point 35 (19/48) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.755 ( -1.779 3.082 0.000) 3.559 1.799 ( -3.294 5.705 0.000) 6.588 1.949 ( -0.634 1.099 0.000) 1.268 1.984 ( 0.346 -0.599 0.000) 0.692 2.063 ( -1.299 2.249 0.000) 2.597 2.207 ( 1.215 -2.104 0.000) 2.429 2.370 ( -0.085 0.147 0.000) 0.170 2.520 ( 1.744 -3.020 0.000) 3.487 2.710 ( -8.902 15.419 0.000) 17.804 2.750 ( -8.768 15.187 0.000) 17.536 3.971 ( 1.836 -3.180 0.000) 3.672 4.016 ( 1.776 -3.076 0.000) 3.552 5.838 ( -5.142 8.906 0.000) 10.283 5.850 ( -5.121 8.870 0.000) 10.242 5.948 ( 0.908 -1.573 0.000) 1.816 6.015 ( 1.845 -3.195 0.000) 3.690 6.090 ( 0.420 -0.727 0.000) 0.839 6.247 ( 1.760 -3.048 0.000) 3.520 6.421 ( -4.365 7.560 0.000) 8.729 6.435 ( -3.310 5.733 0.000) 6.620 6.767 ( -0.165 0.285 0.000) 0.329 6.790 ( 0.490 -0.849 0.000) 0.980 7.875 ( 2.053 -3.556 0.000) 4.106 7.949 ( 2.569 -4.450 0.000) 5.139 ======================= Grid point 47 (20/48) ======================= q-point: ( 0.21 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.579 ( 3.695 6.400 0.000) 7.390 1.808 ( 2.368 4.101 0.000) 4.735 1.838 ( 3.730 6.460 0.000) 7.459 1.907 ( 4.136 7.164 0.000) 8.272 1.955 ( 0.344 0.597 0.000) 0.689 2.140 ( 0.539 0.934 0.000) 1.078 2.384 ( 2.096 3.631 0.000) 4.193 2.504 ( 3.524 6.104 0.000) 7.048 2.746 ( 2.876 4.982 0.000) 5.752 2.796 ( 2.581 4.470 0.000) 5.162 3.573 ( 7.744 13.413 0.000) 15.488 3.620 ( 5.816 10.074 0.000) 11.633 5.772 ( 0.820 1.420 0.000) 1.640 5.809 ( 2.040 3.533 0.000) 4.080 5.851 ( 1.731 2.998 0.000) 3.462 5.859 ( 4.141 7.173 0.000) 8.282 6.141 ( -0.955 -1.655 0.000) 1.911 6.251 ( 1.730 2.996 0.000) 3.459 6.612 ( 0.366 0.634 0.000) 0.732 6.778 ( -0.466 -0.808 0.000) 0.933 6.787 ( -6.358 -11.013 0.000) 12.717 7.006 ( -3.078 -5.331 0.000) 6.156 7.884 ( -2.309 -3.999 0.000) 4.618 7.938 ( -3.009 -5.212 0.000) 6.019 ======================= Grid point 48 (21/48) ======================= q-point: ( 0.29 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.718 ( 3.899 6.555 0.000) 7.627 1.869 ( -0.309 4.665 0.000) 4.675 1.934 ( 1.864 2.870 0.000) 3.422 1.965 ( 0.169 0.443 0.000) 0.474 2.029 ( 2.666 4.690 0.000) 5.395 2.152 ( 0.738 -0.298 0.000) 0.796 2.394 ( -0.741 -0.756 0.000) 1.058 2.529 ( -1.186 -2.816 0.000) 3.056 2.899 ( 0.702 14.994 0.000) 15.011 2.939 ( 0.953 14.562 0.000) 14.594 3.787 ( 5.768 2.119 0.000) 6.145 3.799 ( 7.337 4.147 0.000) 8.427 5.824 ( 2.718 2.601 0.000) 3.762 5.841 ( 2.740 0.543 0.000) 2.793 5.947 ( 1.057 5.108 0.000) 5.216 5.963 ( 0.113 5.675 0.000) 5.677 6.098 ( -1.138 -2.453 0.000) 2.704 6.274 ( -1.689 -0.862 0.000) 1.896 6.614 ( -3.664 4.915 0.000) 6.131 6.623 ( -3.773 -1.117 0.000) 3.935 6.770 ( 0.865 -1.621 0.000) 1.837 6.899 ( -5.260 -1.523 0.000) 5.476 7.826 ( 0.368 -4.004 0.000) 4.020 7.868 ( 1.028 -5.280 0.000) 5.380 ======================= Grid point 49 (22/48) ======================= q-point: ( 0.36 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.821 ( -0.125 5.795 0.000) 5.796 1.897 ( -1.636 4.319 0.000) 4.618 1.967 ( -0.170 0.842 0.000) 0.859 1.973 ( 0.541 0.179 0.000) 0.570 2.094 ( 0.531 2.305 0.000) 2.365 2.165 ( 1.403 -0.990 0.000) 1.717 2.377 ( -0.503 -0.070 0.000) 0.508 2.465 ( -0.248 -4.683 0.000) 4.689 3.025 ( -6.176 17.701 0.000) 18.748 3.067 ( -6.334 18.147 0.000) 19.220 3.865 ( 4.862 -4.913 0.000) 6.913 3.911 ( 5.369 -3.956 0.000) 6.669 5.905 ( 2.855 2.391 0.000) 3.724 5.918 ( 4.277 1.247 0.000) 4.455 5.976 ( -3.075 2.391 0.000) 3.895 5.978 ( -1.501 1.340 0.000) 2.012 6.070 ( 0.490 -1.559 0.000) 1.634 6.197 ( -1.098 -5.106 0.000) 5.223 6.592 ( -5.251 7.350 0.000) 9.033 6.619 ( -3.872 10.052 0.000) 10.772 6.771 ( 0.688 -1.040 0.000) 1.247 6.799 ( -4.121 -0.233 0.000) 4.127 7.809 ( 1.448 -3.016 0.000) 3.346 7.853 ( 2.548 -4.731 0.000) 5.373 ======================= Grid point 62 (23/48) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.857 ( 3.649 6.319 0.000) 7.297 1.934 ( 1.097 1.900 0.000) 2.194 1.968 ( 0.463 0.801 0.000) 0.925 1.973 ( 0.145 0.251 0.000) 0.290 2.097 ( 1.093 1.893 0.000) 2.186 2.148 ( -0.020 -0.035 0.000) 0.040 2.382 ( -0.181 -0.313 0.000) 0.361 2.442 ( -2.598 -4.500 0.000) 5.196 3.199 ( 7.731 13.390 0.000) 15.462 3.247 ( 8.094 14.019 0.000) 16.188 3.769 ( -2.053 -3.556 0.000) 4.107 3.815 ( -1.334 -2.311 0.000) 2.669 5.877 ( 1.704 2.952 0.000) 3.408 5.890 ( 2.108 3.652 0.000) 4.217 6.002 ( 0.485 0.840 0.000) 0.970 6.022 ( 0.255 0.442 0.000) 0.510 6.049 ( -1.162 -2.013 0.000) 2.325 6.180 ( -3.592 -6.221 0.000) 7.183 6.693 ( 3.266 5.657 0.000) 6.532 6.703 ( 1.601 2.772 0.000) 3.201 6.744 ( -0.253 -0.438 0.000) 0.506 6.866 ( -1.230 -2.130 0.000) 2.459 7.777 ( -0.557 -0.965 0.000) 1.115 7.797 ( -1.092 -1.892 0.000) 2.184 ======================= Grid point 63 (24/48) ======================= q-point: ( 0.36 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.931 ( -2.350 4.070 0.000) 4.700 1.956 ( -0.778 1.347 0.000) 1.556 1.975 ( -0.090 0.155 0.000) 0.179 1.975 ( -0.024 0.041 0.000) 0.047 2.110 ( 0.250 -0.433 0.000) 0.500 2.149 ( 0.268 -0.464 0.000) 0.535 2.379 ( -0.089 0.154 0.000) 0.178 2.388 ( 1.297 -2.246 0.000) 2.594 3.368 ( -8.520 14.757 0.000) 17.040 3.423 ( -8.844 15.318 0.000) 17.688 3.702 ( 5.876 -10.177 0.000) 11.751 3.762 ( 5.733 -9.931 0.000) 11.467 5.925 ( 0.180 -0.312 0.000) 0.360 5.946 ( -0.831 1.440 0.000) 1.663 6.001 ( -0.547 0.947 0.000) 1.093 6.006 ( -1.230 2.131 0.000) 2.461 6.035 ( 1.135 -1.966 0.000) 2.271 6.096 ( 2.388 -4.136 0.000) 4.776 6.717 ( -2.235 3.872 0.000) 4.470 6.754 ( 0.434 -0.751 0.000) 0.867 6.800 ( -0.732 1.268 0.000) 1.464 6.800 ( -2.538 4.395 0.000) 5.075 7.771 ( 0.480 -0.832 0.000) 0.960 7.780 ( 1.205 -2.087 0.000) 2.410 ======================= Grid point 211 (25/48) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.463 ( 0.000 -0.000 21.304) 21.304 0.463 ( -0.000 0.000 21.304) 21.304 0.861 ( 0.000 0.000 40.095) 40.095 0.877 ( 0.000 -0.000 -17.354) 17.354 0.877 ( 0.000 0.000 -17.354) 17.354 1.645 ( -0.000 -0.000 8.801) 8.801 1.645 ( -0.000 0.000 8.801) 8.801 1.666 ( -0.000 -0.000 -35.622) 35.622 1.796 ( 0.000 0.000 -5.073) 5.073 1.796 ( -0.000 0.000 -5.073) 5.073 3.088 ( 0.000 -0.000 22.343) 22.343 3.468 ( 0.000 -0.000 -12.853) 12.853 4.912 ( 0.000 0.000 0.665) 0.665 4.912 ( -0.000 0.000 0.665) 0.665 4.928 ( 0.000 -0.000 -0.677) 0.677 4.928 ( 0.000 0.000 -0.677) 0.677 6.314 ( 0.000 0.000 0.436) 0.436 6.314 ( 0.000 -0.000 0.436) 0.436 6.324 ( 0.000 -0.000 -0.427) 0.427 6.324 ( 0.000 -0.000 -0.427) 0.427 8.074 ( -0.000 0.000 2.490) 2.490 8.134 ( 0.000 -0.000 -2.729) 2.729 8.377 ( 0.000 0.000 0.716) 0.716 8.389 ( 0.000 0.000 -0.416) 0.416 ======================= Grid point 212 (26/48) ======================= q-point: ( 0.07 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.669 ( 14.017 8.092 14.639) 21.823 0.708 ( 12.387 7.152 17.077) 22.276 1.003 ( 9.402 5.428 -15.543) 18.959 1.066 ( 11.394 6.579 -16.383) 21.013 1.126 ( 21.682 12.518 26.438) 36.411 1.666 ( 1.806 1.043 8.878) 9.120 1.686 ( 6.374 3.680 2.170) 7.674 1.776 ( 7.627 4.403 -22.518) 24.179 1.815 ( 1.621 0.936 -4.954) 5.296 1.871 ( 6.276 3.624 -5.450) 9.068 2.933 ( -7.374 -4.257 15.194) 17.418 3.110 ( -9.611 -5.549 -1.031) 11.145 4.963 ( 4.229 2.442 0.567) 4.916 4.977 ( 4.043 2.334 -0.576) 4.704 5.509 ( 37.335 21.555 12.894) 44.997 6.013 ( 39.264 22.669 -37.178) 58.632 6.266 ( -3.386 -1.955 0.059) 3.910 6.276 ( -3.799 -2.193 -0.507) 4.415 6.293 ( -1.684 -0.973 0.426) 1.991 6.303 ( -1.704 -0.984 -0.418) 2.012 8.012 ( -6.825 -3.940 3.853) 8.772 8.086 ( -4.706 -2.717 -3.158) 6.285 8.241 ( -9.094 -5.251 1.446) 10.600 8.256 ( -5.835 -3.369 0.146) 6.739 ======================= Grid point 213 (27/48) ======================= q-point: ( 0.14 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.993 ( 12.551 7.246 10.117) 17.674 1.005 ( 13.397 7.735 9.375) 18.089 1.248 ( 10.218 5.900 -12.706) 17.340 1.258 ( 5.587 3.226 -13.959) 15.378 1.665 ( 19.472 11.242 14.441) 26.723 1.734 ( 4.149 2.395 8.560) 9.810 1.873 ( 3.289 1.899 -4.490) 5.881 1.874 ( 6.501 3.753 7.333) 10.494 2.084 ( 18.287 10.558 -14.037) 25.356 2.115 ( 15.965 9.217 -8.214) 20.182 2.852 ( 0.600 0.346 9.584) 9.609 2.982 ( -3.057 -1.765 -2.865) 4.546 5.098 ( 6.797 3.924 0.327) 7.855 5.106 ( 6.529 3.769 -0.331) 7.546 6.127 ( -5.876 -3.392 1.404) 6.929 6.152 ( -4.176 -2.411 -0.822) 4.892 6.241 ( -2.566 -1.481 0.386) 2.988 6.250 ( -2.626 -1.516 -0.383) 3.056 6.263 ( 21.067 12.163 10.861) 26.640 6.569 ( 11.867 6.851 -15.839) 20.943 7.764 ( -12.697 -7.331 3.770) 15.139 7.839 ( -14.869 -8.585 -2.993) 17.429 8.103 ( -2.165 -1.250 4.001) 4.718 8.206 ( -0.079 -0.046 -5.187) 5.187 ======================= Grid point 214 (28/48) ======================= q-point: ( 0.21 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.275 ( 9.428 5.443 6.540) 12.699 1.294 ( 12.686 7.324 2.922) 14.938 1.396 ( 6.851 3.955 -6.244) 10.078 1.448 ( 6.574 3.796 -9.082) 11.837 1.864 ( 6.684 3.859 6.306) 9.967 1.958 ( 6.066 3.502 4.919) 8.559 1.970 ( 4.864 2.808 -3.522) 6.629 1.987 ( 4.369 2.522 3.412) 6.090 2.522 ( 14.886 8.594 -3.561) 17.554 2.549 ( 17.477 10.090 -5.517) 20.921 2.936 ( 4.321 2.495 4.442) 6.680 2.988 ( 8.134 4.696 -0.457) 9.403 5.259 ( 6.408 3.700 0.085) 7.400 5.261 ( 6.219 3.591 -0.086) 7.182 5.979 ( -5.810 -3.354 2.134) 7.040 6.027 ( -4.945 -2.855 -1.998) 6.049 6.184 ( -1.983 -1.145 0.282) 2.307 6.191 ( -2.136 -1.233 -0.281) 2.482 6.577 ( 5.623 3.246 6.125) 8.925 6.730 ( 3.222 1.860 -7.181) 8.087 7.441 ( -14.497 -8.370 2.661) 16.950 7.489 ( -13.952 -8.055 -1.657) 16.195 8.098 ( 1.122 0.648 4.560) 4.741 8.212 ( 0.142 0.082 -5.265) 5.268 ======================= Grid point 215 (29/48) ======================= q-point: ( 0.29 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.454 ( 5.994 3.461 4.394) 8.198 1.565 ( 3.636 2.099 -5.323) 6.779 1.569 ( 10.091 5.826 -0.182) 11.654 1.583 ( 9.008 5.201 -1.193) 10.470 2.011 ( 0.807 0.466 4.555) 4.650 2.032 ( 7.234 4.177 3.308) 8.985 2.090 ( 2.987 1.725 -2.138) 4.058 2.095 ( 5.590 3.227 -2.283) 6.846 2.710 ( -1.364 -0.788 -0.642) 1.701 2.724 ( 3.364 1.942 -1.561) 4.186 3.233 ( 18.635 10.759 4.039) 21.894 3.311 ( 16.314 9.419 -2.869) 19.055 5.383 ( 4.115 2.376 0.037) 4.752 5.384 ( 4.180 2.414 -0.037) 4.827 5.879 ( -2.291 -1.323 2.769) 3.830 5.944 ( -1.716 -0.991 -2.856) 3.477 6.153 ( -0.791 -0.456 0.077) 0.916 6.155 ( -1.024 -0.591 -0.077) 1.185 6.611 ( -1.370 -0.791 5.311) 5.542 6.749 ( -1.286 -0.743 -6.925) 7.083 7.120 ( -12.076 -6.972 3.815) 14.457 7.184 ( -11.932 -6.889 -2.085) 13.935 8.120 ( 0.159 0.092 3.331) 3.336 8.199 ( -1.440 -0.831 -3.408) 3.792 ======================= Grid point 216 (30/48) ======================= q-point: ( 0.36 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.566 ( 3.675 2.122 2.664) 5.010 1.630 ( 2.104 1.215 -2.906) 3.788 1.773 ( 7.646 4.414 0.232) 8.832 1.786 ( 7.277 4.202 -0.687) 8.431 2.024 ( 0.532 0.307 4.459) 4.501 2.135 ( 1.212 0.700 -5.313) 5.495 2.189 ( 5.885 3.398 1.550) 6.970 2.221 ( 4.929 2.846 -1.273) 5.833 2.630 ( -4.555 -2.630 4.916) 7.199 2.726 ( -2.644 -1.526 -3.554) 4.685 3.652 ( 15.898 9.179 1.056) 18.388 3.676 ( 13.888 8.018 -1.014) 16.069 5.461 ( 2.675 1.545 0.094) 3.091 5.463 ( 2.744 1.584 -0.094) 3.169 5.883 ( 2.580 1.490 3.287) 4.436 5.960 ( 3.125 1.804 -3.383) 4.947 6.134 ( -0.933 -0.539 0.141) 1.087 6.137 ( -0.760 -0.439 -0.140) 0.889 6.544 ( -4.083 -2.357 4.862) 6.773 6.671 ( -5.653 -3.264 -6.594) 9.278 6.900 ( -6.180 -3.568 3.020) 7.749 6.945 ( -7.990 -4.613 -1.365) 9.327 8.092 ( -2.437 -1.407 2.268) 3.614 8.143 ( -3.141 -1.814 -2.100) 4.191 ======================= Grid point 217 (31/48) ======================= q-point: ( 0.43 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.629 ( 1.723 0.995 1.631) 2.573 1.668 ( 1.110 0.641 -1.697) 2.126 1.901 ( 3.082 1.780 0.688) 3.625 1.923 ( 4.050 2.338 -1.179) 4.823 2.049 ( 1.675 0.967 4.496) 4.895 2.166 ( 1.660 0.958 -6.260) 6.547 2.299 ( 3.342 1.930 0.769) 3.935 2.316 ( 2.938 1.696 -0.689) 3.461 2.523 ( -3.972 -2.293 6.047) 7.590 2.629 ( -4.762 -2.749 -3.667) 6.609 3.938 ( 7.956 4.593 0.529) 9.201 3.950 ( 8.997 5.195 -0.520) 10.402 5.514 ( 1.777 1.026 0.181) 2.060 5.518 ( 1.871 1.080 -0.185) 2.168 5.984 ( 5.311 3.067 4.405) 7.551 6.082 ( 6.817 3.936 -4.026) 8.842 6.108 ( -1.100 -0.635 0.255) 1.295 6.114 ( -1.046 -0.604 -0.252) 1.234 6.429 ( -5.258 -3.036 2.511) 6.570 6.497 ( -8.254 -4.765 -3.375) 10.111 6.814 ( -3.018 -1.742 0.215) 3.492 6.818 ( -1.972 -1.139 -0.044) 2.277 8.025 ( -2.744 -1.584 1.990) 3.742 8.067 ( -2.838 -1.639 -1.683) 3.684 ======================= Grid point 218 (32/48) ======================= q-point: (-0.50 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 89 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.649 ( 0.000 0.000 1.383) 1.383 1.681 ( 0.000 0.000 -1.419) 1.419 1.933 ( 0.000 0.000 1.373) 1.373 1.968 ( 0.000 0.000 -1.753) 1.753 2.078 ( 0.000 0.000 4.490) 4.490 2.195 ( 0.000 0.000 -6.364) 6.364 2.340 ( 0.000 0.000 0.554) 0.554 2.352 ( 0.000 0.000 -0.507) 0.507 2.472 ( 0.000 0.000 5.163) 5.163 2.557 ( 0.000 0.000 -2.748) 2.748 4.034 ( 0.000 0.000 1.063) 1.063 4.058 ( 0.000 0.000 -1.041) 1.041 5.537 ( 0.000 0.000 0.237) 0.237 5.542 ( 0.000 0.000 -0.243) 0.243 6.070 ( 0.000 0.000 6.194) 6.194 6.093 ( 0.000 0.000 0.277) 0.277 6.099 ( 0.000 0.000 -0.271) 0.271 6.221 ( 0.000 0.000 -7.079) 7.079 6.334 ( -0.000 -0.000 1.042) 1.042 6.355 ( -0.000 -0.000 -0.600) 0.600 6.777 ( 0.000 0.000 0.971) 0.971 6.799 ( 0.000 0.000 -0.900) 0.900 7.988 ( 0.000 0.000 1.999) 1.999 8.030 ( 0.000 0.000 -1.637) 1.637 ======================= Grid point 227 (33/48) ======================= q-point: ( 0.07 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.905 ( 6.733 11.661 11.903) 17.972 0.942 ( 9.086 15.738 9.865) 20.677 1.203 ( 7.062 12.232 -12.913) 19.137 1.211 ( 3.457 5.988 -15.193) 16.692 1.513 ( 11.828 20.487 17.182) 29.238 1.731 ( 3.391 5.874 8.746) 11.068 1.811 ( 3.341 5.787 7.599) 10.119 1.867 ( 2.976 5.155 -3.019) 6.675 1.975 ( 9.231 15.989 -20.472) 27.568 2.015 ( 6.523 11.298 -6.511) 14.581 2.853 ( -1.001 -1.734 11.238) 11.415 3.003 ( -2.349 -4.069 -3.136) 5.649 5.065 ( 4.167 7.217 0.358) 8.341 5.073 ( 3.964 6.866 -0.359) 7.936 6.085 ( 15.277 26.460 12.016) 32.831 6.165 ( -3.447 -5.971 0.771) 6.937 6.197 ( -2.290 -3.966 0.237) 4.586 6.259 ( -1.258 -2.178 0.357) 2.541 6.268 ( -1.336 -2.314 -0.345) 2.694 6.479 ( 9.295 16.100 -20.863) 27.944 7.844 ( -6.784 -11.751 4.485) 14.291 7.935 ( -7.521 -13.027 -3.903) 15.540 8.119 ( -2.458 -4.258 3.402) 5.979 8.199 ( -0.856 -1.482 -3.956) 4.311 ======================= Grid point 228 (34/48) ======================= q-point: ( 0.14 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.166 ( 9.593 10.092 3.828) 14.440 1.231 ( 7.226 10.570 3.965) 13.404 1.334 ( 4.776 6.807 -5.583) 10.016 1.437 ( 5.426 10.360 -8.887) 14.688 1.792 ( 6.993 5.908 7.320) 11.721 1.893 ( 3.278 4.572 2.373) 6.106 1.932 ( 6.081 7.139 4.490) 10.398 1.985 ( 2.689 6.381 0.516) 6.943 2.365 ( 15.212 16.270 -5.699) 22.992 2.375 ( 14.943 14.955 -6.184) 22.027 2.891 ( 4.162 4.138 4.725) 7.534 2.947 ( -0.224 -1.204 -0.540) 1.339 5.236 ( 4.974 9.311 0.043) 10.556 5.237 ( 4.855 8.853 -0.032) 10.097 6.025 ( -4.774 -6.061 1.765) 7.915 6.064 ( -4.103 -5.471 -1.581) 7.018 6.209 ( -1.917 -1.486 0.333) 2.448 6.216 ( -1.947 -1.623 -0.260) 2.548 6.512 ( 8.424 9.942 7.011) 14.797 6.687 ( 4.245 4.050 -8.148) 10.041 7.571 ( -10.562 -11.113 2.370) 15.514 7.618 ( -11.759 -12.481 -1.794) 17.241 8.074 ( 0.307 -1.544 4.160) 4.448 8.182 ( 0.706 -2.328 -5.192) 5.734 ======================= Grid point 229 (35/48) ======================= q-point: ( 0.21 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.400 ( 6.431 6.204 1.303) 9.031 1.440 ( 8.775 5.163 -2.335) 10.446 1.517 ( 4.822 10.750 2.628) 12.071 1.604 ( 3.001 8.904 -4.430) 10.388 1.959 ( 3.955 3.146 2.776) 5.766 1.978 ( 4.936 2.575 1.144) 5.684 2.075 ( 2.524 6.525 2.349) 7.380 2.106 ( 2.843 8.221 -1.139) 8.773 2.666 ( 6.557 5.019 -0.696) 8.287 2.704 ( 3.661 3.675 -4.737) 7.025 3.048 ( 13.211 9.841 5.741) 17.446 3.138 ( 13.726 9.769 -2.685) 17.060 5.401 ( 2.407 9.863 0.163) 10.154 5.405 ( 2.306 10.240 -0.157) 10.498 5.898 ( -3.058 -4.290 2.354) 5.770 5.951 ( -2.426 -4.304 -2.290) 5.446 6.170 ( -1.372 -0.235 0.272) 1.418 6.176 ( -1.427 -0.099 -0.233) 1.449 6.626 ( -0.107 1.415 4.833) 5.037 6.749 ( 0.791 -0.053 -5.879) 5.932 7.257 ( -12.204 -9.301 2.670) 15.575 7.304 ( -10.905 -9.776 -1.555) 14.728 8.080 ( 2.043 -2.588 3.545) 4.842 8.166 ( 1.206 -4.487 -3.854) 6.036 ======================= Grid point 230 (36/48) ======================= q-point: ( 0.29 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.551 ( 4.166 4.996 1.830) 6.758 1.592 ( 3.712 3.552 -1.612) 5.384 1.716 ( 6.128 7.925 0.670) 10.040 1.747 ( 3.670 8.146 -1.600) 9.077 2.000 ( 1.502 -0.540 3.747) 4.073 2.084 ( 3.328 0.318 -3.731) 5.010 2.202 ( 2.312 10.078 1.954) 10.523 2.240 ( 1.949 9.518 -1.479) 9.828 2.681 ( -3.741 -1.702 2.478) 4.800 2.738 ( -0.045 -0.440 -2.403) 2.443 3.459 ( 17.069 11.021 2.137) 20.430 3.506 ( 15.004 9.429 -1.916) 17.824 5.516 ( -0.791 10.235 0.155) 10.266 5.520 ( -0.969 10.381 -0.157) 10.428 5.849 ( 1.109 -1.431 2.590) 3.160 5.909 ( 1.749 -2.289 -2.600) 3.880 6.153 ( -0.974 0.361 0.354) 1.098 6.161 ( -1.199 1.291 -0.391) 1.805 6.590 ( -3.342 -1.223 4.857) 6.022 6.720 ( -2.982 -2.147 -6.831) 7.757 6.984 ( -8.905 -6.268 3.631) 11.479 7.040 ( -9.045 -6.809 -1.648) 11.440 8.072 ( 0.946 -4.657 2.411) 5.329 8.127 ( 0.327 -5.800 -2.369) 6.274 ======================= Grid point 231 (37/48) ======================= q-point: ( 0.36 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.646 ( 0.875 5.014 0.912) 5.171 1.665 ( 1.207 2.894 -0.775) 3.229 1.870 ( 3.819 4.437 0.390) 5.867 1.888 ( 3.857 5.159 -1.028) 6.523 2.021 ( 1.633 0.118 4.114) 4.427 2.126 ( 2.283 -0.484 -5.334) 5.822 2.334 ( 0.832 9.450 1.262) 9.570 2.359 ( 0.255 9.671 -1.118) 9.739 2.582 ( -5.106 -1.269 5.205) 7.400 2.680 ( -4.143 -2.418 -3.541) 5.962 3.818 ( 12.579 6.815 0.118) 14.307 3.821 ( 11.144 5.934 -0.126) 12.626 5.590 ( -2.623 10.386 0.067) 10.713 5.591 ( -2.628 10.228 -0.079) 10.560 5.906 ( 5.245 0.839 2.957) 6.080 5.968 ( 6.246 -0.747 -2.350) 6.715 6.134 ( -1.404 -0.002 0.583) 1.520 6.153 ( -1.653 2.130 -1.058) 2.896 6.495 ( -5.228 -2.134 3.484) 6.636 6.587 ( -7.485 -4.825 -4.581) 10.015 6.840 ( -2.895 -2.328 1.411) 3.974 6.862 ( -5.097 -3.318 -0.642) 6.115 8.022 ( -0.581 -4.983 1.938) 5.378 8.065 ( -0.656 -5.343 -1.728) 5.654 ======================= Grid point 232 (38/48) ======================= q-point: ( 0.43 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.687 ( -1.409 4.320 0.523) 4.574 1.700 ( -0.432 2.518 -0.553) 2.614 1.936 ( 0.574 1.507 1.118) 1.962 1.965 ( 1.113 1.381 -1.499) 2.322 2.063 ( 1.604 0.908 3.943) 4.352 2.164 ( 2.124 -0.456 -5.323) 5.749 2.412 ( -0.878 5.523 1.719) 5.850 2.434 ( -2.955 8.805 -0.663) 9.312 2.505 ( -3.600 3.471 3.408) 6.052 2.575 ( -4.136 -1.768 -2.567) 5.179 4.000 ( 4.568 1.276 0.881) 4.824 4.020 ( 5.108 1.488 -0.871) 5.392 5.634 ( -4.249 9.835 -0.001) 10.714 5.635 ( -4.051 9.442 -0.048) 10.274 6.015 ( 5.010 -0.059 3.910) 6.356 6.062 ( 3.574 -1.806 -0.471) 4.032 6.111 ( -0.851 -0.001 1.125) 1.410 6.180 ( 3.058 2.989 -4.381) 6.122 6.383 ( -4.842 -1.002 0.650) 4.987 6.405 ( -6.713 -3.170 -1.179) 7.517 6.785 ( -1.462 -1.006 0.677) 1.899 6.799 ( -0.513 -0.798 -0.582) 1.113 7.971 ( 0.145 -3.399 1.931) 3.911 8.012 ( 0.225 -3.490 -1.633) 3.860 ======================= Grid point 242 (39/48) ======================= q-point: ( 0.14 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.371 ( 3.435 5.950 -0.438) 6.884 1.386 ( 6.279 10.875 -1.982) 12.713 1.552 ( 6.428 11.134 2.769) 13.151 1.622 ( 4.398 7.617 -2.477) 9.137 1.938 ( 4.197 7.270 0.550) 8.412 1.945 ( 1.036 1.794 1.497) 2.556 2.089 ( 4.522 7.833 3.688) 9.768 2.132 ( 4.592 7.954 -1.022) 9.241 2.627 ( 5.124 8.875 -0.491) 10.260 2.681 ( 5.551 9.614 -5.902) 12.573 2.985 ( 5.093 8.821 5.346) 11.503 3.058 ( 7.182 12.439 -1.790) 14.474 5.441 ( 5.970 10.340 0.334) 11.944 5.448 ( 6.141 10.637 -0.319) 12.286 5.902 ( -3.166 -5.484 2.110) 6.675 5.950 ( -3.046 -5.276 -2.002) 6.413 6.184 ( -0.521 -0.902 0.362) 1.103 6.190 ( -0.424 -0.735 -0.179) 0.868 6.639 ( 1.478 2.560 3.840) 4.846 6.734 ( 0.540 0.936 -4.457) 4.586 7.326 ( -7.376 -12.775 1.559) 14.833 7.355 ( -7.358 -12.744 -0.979) 14.748 8.040 ( -1.067 -1.848 3.189) 3.837 8.119 ( -2.061 -3.570 -3.634) 5.495 ======================= Grid point 243 (40/48) ======================= q-point: ( 0.21 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.505 ( 5.026 5.273 2.740) 7.783 1.585 ( 4.870 7.788 -4.538) 10.245 1.735 ( 2.673 9.553 2.369) 10.199 1.768 ( 4.017 7.065 -0.423) 8.139 1.974 ( 1.146 0.371 2.732) 2.985 2.044 ( 2.418 2.133 -3.794) 4.979 2.262 ( 3.878 10.615 2.658) 11.610 2.302 ( 3.587 10.592 -1.214) 11.249 2.724 ( -1.860 -0.169 -0.307) 1.893 2.731 ( 1.297 1.888 -0.787) 2.422 3.307 ( 13.475 14.136 2.930) 19.748 3.368 ( 11.952 12.163 -2.371) 17.216 5.629 ( 2.344 11.436 0.399) 11.681 5.639 ( 2.134 11.511 -0.413) 11.714 5.823 ( -0.664 -2.604 1.932) 3.310 5.868 ( -0.704 -3.491 -1.961) 4.066 6.173 ( -0.612 0.317 0.900) 1.133 6.190 ( -0.158 1.453 -0.592) 1.577 6.636 ( -2.273 -0.246 3.473) 4.158 6.732 ( -0.711 -1.573 -4.981) 5.271 7.065 ( -8.431 -9.510 2.940) 13.045 7.115 ( -7.614 -8.140 -1.584) 11.258 7.998 ( 0.342 -5.033 2.049) 5.445 8.047 ( -0.048 -6.579 -2.150) 6.921 ======================= Grid point 244 (41/48) ======================= q-point: ( 0.29 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 296 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.641 ( 4.075 4.809 2.812) 6.903 1.708 ( 1.165 6.359 -3.012) 7.133 1.863 ( 2.812 5.926 1.346) 6.696 1.893 ( 2.405 5.683 -1.276) 6.302 2.000 ( 1.423 0.924 3.386) 3.787 2.089 ( 1.753 0.456 -4.680) 5.019 2.408 ( 1.333 7.249 2.647) 7.832 2.452 ( 1.266 9.131 -1.424) 9.328 2.685 ( -4.023 4.580 2.389) 6.547 2.737 ( -2.421 2.436 -2.128) 4.040 3.665 ( 12.911 8.744 0.747) 15.611 3.682 ( 11.528 7.450 -0.749) 13.746 5.751 ( -1.253 10.747 0.088) 10.820 5.752 ( -1.296 9.788 0.023) 9.873 5.838 ( 2.820 1.672 1.317) 3.533 5.869 ( 3.443 -0.798 -1.305) 3.767 6.159 ( -1.312 -0.114 1.856) 2.276 6.202 ( -1.098 2.355 -1.855) 3.192 6.573 ( -4.552 0.155 2.940) 5.421 6.658 ( -5.068 -3.569 -4.874) 7.886 6.872 ( -4.509 -4.431 3.499) 7.225 6.928 ( -6.417 -4.000 -1.749) 7.761 7.952 ( 1.145 -6.441 1.568) 6.728 7.987 ( 1.327 -7.244 -1.497) 7.515 ======================= Grid point 245 (42/48) ======================= q-point: ( 0.36 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.734 ( 0.902 4.581 1.430) 4.883 1.768 ( -1.171 5.744 -1.476) 6.046 1.942 ( 1.100 2.401 0.692) 2.730 1.959 ( 1.137 1.520 -0.808) 2.063 2.052 ( 1.815 2.947 2.811) 4.458 2.128 ( 2.015 0.673 -4.032) 4.557 2.439 ( -0.717 -0.706 3.730) 3.863 2.507 ( -0.427 0.153 -2.445) 2.487 2.698 ( -5.852 13.915 1.438) 15.164 2.729 ( -6.014 11.610 -1.301) 13.139 3.903 ( 7.793 1.745 0.516) 8.003 3.916 ( 8.601 2.364 -0.525) 8.935 5.813 ( -2.576 9.131 0.182) 9.489 5.817 ( -2.744 9.542 -0.221) 9.932 5.928 ( 4.371 1.840 1.662) 5.026 5.954 ( 5.565 0.052 -0.707) 5.610 6.126 ( -1.382 -1.032 2.283) 2.861 6.194 ( -1.380 0.886 -3.515) 3.879 6.487 ( -6.241 2.386 1.102) 6.772 6.512 ( -7.549 -0.594 -1.081) 7.649 6.802 ( -1.090 -1.519 0.539) 1.946 6.812 ( -3.383 -1.653 -0.387) 3.785 7.912 ( 1.311 -5.332 1.590) 5.716 7.947 ( 1.676 -5.766 -1.437) 6.174 ======================= Grid point 246 (43/48) ======================= q-point: ( 0.43 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.765 ( -2.123 3.678 0.860) 4.333 1.786 ( -2.870 4.970 -1.019) 5.829 1.960 ( -0.475 0.823 0.841) 1.269 1.977 ( 0.001 -0.001 -0.638) 0.638 2.089 ( -0.797 1.381 2.387) 2.871 2.156 ( 0.261 -0.452 -3.581) 3.619 2.421 ( 0.636 -1.102 3.572) 3.792 2.490 ( 1.311 -2.270 -2.658) 3.733 2.722 ( -8.686 15.044 0.981) 17.399 2.742 ( -8.684 15.042 -0.718) 17.384 3.982 ( 1.819 -3.150 0.975) 3.766 4.005 ( 1.789 -3.098 -0.966) 3.706 5.839 ( -5.064 8.771 0.112) 10.129 5.843 ( -4.965 8.599 -0.354) 9.936 5.987 ( 1.379 -2.388 2.528) 3.741 6.015 ( 1.832 -3.172 -0.277) 3.674 6.105 ( 0.241 -0.418 2.051) 2.107 6.187 ( 0.934 -1.618 -4.704) 5.061 6.428 ( -4.010 6.946 0.501) 8.036 6.436 ( -3.402 5.892 -0.156) 6.805 6.774 ( -0.016 0.027 0.549) 0.550 6.785 ( 0.315 -0.546 -0.421) 0.758 7.895 ( 2.222 -3.849 1.664) 4.746 7.932 ( 2.477 -4.289 -1.486) 5.171 ======================= Grid point 258 (44/48) ======================= q-point: ( 0.21 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.626 ( 3.729 6.459 4.146) 8.533 1.729 ( 3.531 6.116 -5.011) 8.659 1.878 ( 2.680 4.642 1.628) 5.602 1.900 ( 3.067 5.312 -0.635) 6.167 1.990 ( 0.972 1.684 3.130) 3.685 2.076 ( 0.806 1.395 -4.565) 4.841 2.431 ( 2.443 4.232 3.175) 5.828 2.486 ( 3.184 5.515 -1.770) 6.610 2.758 ( 2.854 4.943 1.056) 5.805 2.783 ( 2.695 4.669 -1.091) 5.500 3.585 ( 7.228 12.520 1.011) 14.492 3.608 ( 6.266 10.853 -0.991) 12.571 5.783 ( 0.761 1.318 0.887) 1.762 5.802 ( 1.222 2.116 -0.723) 2.549 5.850 ( 2.860 4.955 0.052) 5.721 5.855 ( 3.912 6.776 -0.313) 7.831 6.171 ( -0.461 -0.799 2.531) 2.694 6.226 ( 0.899 1.558 -2.198) 2.840 6.627 ( -0.268 -0.464 1.492) 1.585 6.682 ( -1.909 -3.306 -3.812) 5.395 6.894 ( -3.761 -6.515 4.848) 8.949 6.976 ( -3.333 -5.773 -2.765) 7.216 7.898 ( -2.464 -4.267 1.166) 5.064 7.925 ( -2.815 -4.876 -1.151) 5.747 ======================= Grid point 259 (45/48) ======================= q-point: ( 0.29 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.754 ( 3.000 6.123 3.105) 7.492 1.828 ( 1.068 5.247 -3.325) 6.303 1.944 ( 1.195 2.177 0.711) 2.583 1.958 ( 0.617 1.209 -0.566) 1.471 2.053 ( 2.103 3.566 2.076) 4.631 2.109 ( 1.275 1.446 -3.047) 3.606 2.440 ( -1.028 -1.807 3.365) 3.955 2.504 ( -1.201 -2.631 -2.309) 3.701 2.909 ( 0.777 14.871 0.896) 14.918 2.930 ( 0.900 14.655 -0.852) 14.708 3.790 ( 6.165 2.644 0.260) 6.713 3.796 ( 6.949 3.658 -0.263) 7.857 5.829 ( 2.829 2.058 0.387) 3.519 5.837 ( 2.827 1.035 -0.313) 3.027 5.951 ( 0.866 5.125 0.340) 5.209 5.959 ( 0.463 5.444 -0.384) 5.477 6.136 ( -1.524 -2.130 3.390) 4.284 6.222 ( -1.792 -1.460 -4.133) 4.735 6.616 ( -2.849 3.476 0.282) 4.503 6.622 ( -2.961 1.620 -0.187) 3.380 6.807 ( -1.587 -2.075 3.048) 4.014 6.871 ( -4.318 -1.810 -2.480) 5.298 7.837 ( 0.542 -4.335 0.946) 4.470 7.858 ( 0.872 -4.973 -0.892) 5.127 ======================= Grid point 260 (46/48) ======================= q-point: ( 0.36 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.838 ( -0.552 5.417 1.523) 5.654 1.876 ( -1.328 4.668 -1.720) 5.149 1.970 ( 0.129 0.651 0.225) 0.701 1.973 ( 0.443 0.312 -0.037) 0.544 2.108 ( 0.675 1.872 1.294) 2.374 2.143 ( 1.093 0.398 -1.743) 2.095 2.403 ( -0.409 -1.718 2.090) 2.736 2.446 ( -0.304 -3.857 -1.672) 4.215 3.036 ( -6.193 17.760 0.933) 18.832 3.057 ( -6.273 17.985 -0.905) 19.069 3.876 ( 4.987 -4.649 1.021) 6.894 3.900 ( 5.240 -4.172 -0.991) 6.771 5.908 ( 3.029 1.947 0.329) 3.616 5.916 ( 3.763 1.594 -0.288) 4.097 5.977 ( -1.568 1.249 0.277) 2.023 5.983 ( -1.802 1.512 -0.036) 2.352 6.096 ( -0.303 -1.523 2.484) 2.929 6.161 ( -0.952 -3.274 -3.017) 4.553 6.598 ( -5.052 7.170 0.575) 8.789 6.613 ( -4.605 8.612 -0.629) 9.786 6.780 ( -0.384 -0.382 0.749) 0.925 6.795 ( -2.767 0.010 -0.483) 2.809 7.820 ( 1.741 -3.488 0.973) 4.018 7.842 ( 2.290 -4.343 -0.931) 4.998 ======================= Grid point 273 (47/48) ======================= q-point: ( 0.29 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.876 ( 3.039 5.264 1.645) 6.297 1.914 ( 1.765 3.057 -1.710) 3.922 1.969 ( 0.389 0.673 0.106) 0.784 1.972 ( 0.231 0.401 -0.117) 0.477 2.109 ( 0.916 1.586 1.021) 2.097 2.134 ( 0.401 0.695 -1.162) 1.412 2.399 ( -0.966 -1.673 1.402) 2.387 2.429 ( -2.133 -3.695 -1.187) 4.429 3.211 ( 7.799 13.509 1.025) 15.632 3.235 ( 7.978 13.819 -1.058) 15.992 3.781 ( -1.852 -3.208 1.026) 3.843 3.804 ( -1.495 -2.590 -0.967) 3.143 5.879 ( 1.750 3.032 0.249) 3.510 5.886 ( 1.951 3.379 -0.330) 3.915 5.999 ( 0.010 0.017 0.205) 0.206 6.012 ( 0.055 0.095 -0.808) 0.816 6.089 ( -1.070 -1.854 3.035) 3.714 6.152 ( -2.392 -4.144 -2.488) 5.393 6.682 ( 1.843 3.192 -0.071) 3.686 6.693 ( 1.467 2.541 -0.772) 3.034 6.788 ( 0.241 0.418 2.974) 3.013 6.844 ( -0.653 -1.130 -2.013) 2.399 7.782 ( -0.700 -1.212 0.425) 1.463 7.792 ( -0.967 -1.675 -0.405) 1.976 ======================= Grid point 274 (48/48) ======================= q-point: ( 0.36 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.935 ( -1.921 3.328 0.423) 3.866 1.947 ( -1.140 1.975 -0.619) 2.363 1.977 ( -0.079 0.137 0.110) 0.193 1.978 ( -0.107 0.185 0.088) 0.231 2.120 ( 0.247 -0.428 0.826) 0.963 2.139 ( 0.254 -0.440 -0.843) 0.984 2.381 ( 0.267 -0.463 0.187) 0.566 2.386 ( 0.959 -1.660 -0.183) 1.926 3.381 ( -8.547 14.803 1.134) 17.131 3.408 ( -8.690 15.051 -1.250) 17.424 3.718 ( 5.771 -9.996 1.354) 11.621 3.748 ( 5.719 -9.906 -1.218) 11.503 5.931 ( -0.035 0.061 0.412) 0.418 5.939 ( -0.597 1.034 -0.430) 1.269 5.987 ( -0.058 0.099 -0.723) 0.732 5.987 ( 0.010 -0.016 -0.733) 0.733 6.073 ( 0.308 -0.534 2.139) 2.226 6.100 ( 1.298 -2.248 -0.243) 2.608 6.707 ( -1.514 2.623 -0.190) 3.034 6.724 ( -0.571 0.990 -1.660) 2.016 6.815 ( -1.157 2.004 0.788) 2.445 6.817 ( -1.657 2.870 0.643) 3.376 7.773 ( 0.659 -1.142 0.198) 1.333 7.778 ( 1.022 -1.770 -0.199) 2.054 =================== 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/14112 10.0 290.785 290.785 323.996 0.000 -0.000 -0.000 3/14112 20.0 70.152 70.152 57.478 0.000 -0.000 -0.000 3/14112 30.0 43.293 43.293 28.579 0.000 -0.000 -0.000 3/14112 40.0 31.471 31.471 18.714 0.000 -0.000 -0.000 3/14112 50.0 24.122 24.122 13.654 0.000 -0.000 -0.000 3/14112 60.0 19.172 19.172 10.605 0.000 -0.000 -0.000 3/14112 70.0 15.733 15.733 8.604 0.000 -0.000 -0.000 3/14112 80.0 13.268 13.268 7.213 0.000 -0.000 -0.000 3/14112 90.0 11.445 11.445 6.200 0.000 -0.000 -0.000 3/14112 100.0 10.055 10.055 5.435 0.000 -0.000 -0.000 3/14112 110.0 8.966 8.966 4.839 0.000 -0.000 -0.000 3/14112 120.0 8.091 8.091 4.362 0.000 -0.000 -0.000 3/14112 130.0 7.375 7.375 3.973 0.000 -0.000 -0.000 3/14112 140.0 6.778 6.778 3.648 0.000 -0.000 -0.000 3/14112 150.0 6.273 6.273 3.374 0.000 -0.000 -0.000 3/14112 160.0 5.839 5.839 3.139 0.000 -0.000 -0.000 3/14112 170.0 5.463 5.463 2.936 0.000 -0.000 -0.000 3/14112 180.0 5.134 5.134 2.758 0.000 -0.000 -0.000 3/14112 190.0 4.843 4.843 2.601 0.000 -0.000 -0.000 3/14112 200.0 4.584 4.584 2.461 0.000 -0.000 -0.000 3/14112 210.0 4.352 4.352 2.336 0.000 -0.000 -0.000 3/14112 220.0 4.143 4.143 2.223 0.000 -0.000 -0.000 3/14112 230.0 3.954 3.954 2.121 0.000 -0.000 -0.000 3/14112 240.0 3.781 3.781 2.028 0.000 -0.000 -0.000 3/14112 250.0 3.623 3.623 1.943 0.000 -0.000 -0.000 3/14112 260.0 3.478 3.478 1.865 0.000 -0.000 -0.000 3/14112 270.0 3.345 3.345 1.793 0.000 -0.000 -0.000 3/14112 280.0 3.221 3.221 1.727 0.000 -0.000 -0.000 3/14112 290.0 3.106 3.106 1.665 0.000 -0.000 -0.000 3/14112 300.0 3.000 3.000 1.608 0.000 -0.000 -0.000 3/14112 310.0 2.900 2.900 1.554 0.000 -0.000 -0.000 3/14112 320.0 2.807 2.807 1.504 0.000 -0.000 -0.000 3/14112 330.0 2.720 2.720 1.457 0.000 -0.000 -0.000 3/14112 340.0 2.638 2.638 1.413 0.000 -0.000 -0.000 3/14112 350.0 2.561 2.561 1.372 0.000 -0.000 -0.000 3/14112 360.0 2.489 2.489 1.333 0.000 -0.000 -0.000 3/14112 370.0 2.420 2.420 1.296 0.000 -0.000 -0.000 3/14112 380.0 2.355 2.355 1.261 0.000 -0.000 -0.000 3/14112 390.0 2.294 2.294 1.228 0.000 -0.000 -0.000 3/14112 400.0 2.236 2.236 1.197 0.000 -0.000 -0.000 3/14112 410.0 2.180 2.180 1.167 0.000 -0.000 -0.000 3/14112 420.0 2.128 2.128 1.139 0.000 -0.000 -0.000 3/14112 430.0 2.078 2.078 1.112 0.000 -0.000 -0.000 3/14112 440.0 2.030 2.030 1.087 0.000 -0.000 -0.000 3/14112 450.0 1.984 1.984 1.062 0.000 -0.000 -0.000 3/14112 460.0 1.941 1.941 1.039 0.000 -0.000 -0.000 3/14112 470.0 1.899 1.899 1.016 0.000 -0.000 -0.000 3/14112 480.0 1.859 1.859 0.995 0.000 -0.000 -0.000 3/14112 490.0 1.821 1.821 0.974 0.000 -0.000 -0.000 3/14112 500.0 1.784 1.784 0.955 0.000 -0.000 -0.000 3/14112 510.0 1.748 1.748 0.936 0.000 -0.000 -0.000 3/14112 520.0 1.715 1.715 0.918 0.000 -0.000 -0.000 3/14112 530.0 1.682 1.682 0.900 0.000 -0.000 -0.000 3/14112 540.0 1.651 1.651 0.883 0.000 -0.000 -0.000 3/14112 550.0 1.620 1.620 0.867 0.000 -0.000 -0.000 3/14112 560.0 1.591 1.591 0.851 0.000 -0.000 -0.000 3/14112 570.0 1.563 1.563 0.836 0.000 -0.000 -0.000 3/14112 580.0 1.536 1.536 0.822 0.000 -0.000 -0.000 3/14112 590.0 1.510 1.510 0.808 0.000 -0.000 -0.000 3/14112 600.0 1.484 1.484 0.794 0.000 -0.000 -0.000 3/14112 610.0 1.460 1.460 0.781 0.000 -0.000 -0.000 3/14112 620.0 1.436 1.436 0.768 0.000 -0.000 -0.000 3/14112 630.0 1.413 1.413 0.756 0.000 -0.000 -0.000 3/14112 640.0 1.391 1.391 0.744 0.000 -0.000 -0.000 3/14112 650.0 1.370 1.370 0.733 0.000 -0.000 -0.000 3/14112 660.0 1.349 1.349 0.722 0.000 -0.000 -0.000 3/14112 670.0 1.329 1.329 0.711 0.000 -0.000 -0.000 3/14112 680.0 1.309 1.309 0.700 0.000 -0.000 -0.000 3/14112 690.0 1.290 1.290 0.690 0.000 -0.000 -0.000 3/14112 700.0 1.271 1.271 0.680 0.000 -0.000 -0.000 3/14112 710.0 1.253 1.253 0.671 0.000 -0.000 -0.000 3/14112 720.0 1.236 1.236 0.661 0.000 -0.000 -0.000 3/14112 730.0 1.219 1.219 0.652 0.000 -0.000 -0.000 3/14112 740.0 1.202 1.202 0.643 0.000 -0.000 -0.000 3/14112 750.0 1.186 1.186 0.635 0.000 -0.000 -0.000 3/14112 760.0 1.171 1.171 0.626 0.000 -0.000 -0.000 3/14112 770.0 1.155 1.155 0.618 0.000 -0.000 -0.000 3/14112 780.0 1.141 1.141 0.610 0.000 -0.000 -0.000 3/14112 790.0 1.126 1.126 0.602 0.000 -0.000 -0.000 3/14112 800.0 1.112 1.112 0.595 0.000 -0.000 -0.000 3/14112 810.0 1.098 1.098 0.587 0.000 -0.000 -0.000 3/14112 820.0 1.085 1.085 0.580 0.000 -0.000 -0.000 3/14112 830.0 1.072 1.072 0.573 0.000 -0.000 -0.000 3/14112 840.0 1.059 1.059 0.566 0.000 -0.000 -0.000 3/14112 850.0 1.046 1.046 0.560 0.000 -0.000 -0.000 3/14112 860.0 1.034 1.034 0.553 0.000 -0.000 -0.000 3/14112 870.0 1.022 1.022 0.547 0.000 -0.000 -0.000 3/14112 880.0 1.011 1.011 0.541 0.000 -0.000 -0.000 3/14112 890.0 0.999 0.999 0.534 0.000 -0.000 -0.000 3/14112 900.0 0.988 0.988 0.528 0.000 -0.000 -0.000 3/14112 910.0 0.977 0.977 0.523 0.000 -0.000 -0.000 3/14112 920.0 0.967 0.967 0.517 0.000 -0.000 -0.000 3/14112 930.0 0.956 0.956 0.511 0.000 -0.000 -0.000 3/14112 940.0 0.946 0.946 0.506 0.000 -0.000 -0.000 3/14112 950.0 0.936 0.936 0.501 0.000 -0.000 -0.000 3/14112 960.0 0.926 0.926 0.495 0.000 -0.000 -0.000 3/14112 970.0 0.917 0.917 0.490 0.000 -0.000 -0.000 3/14112 980.0 0.907 0.907 0.485 0.000 -0.000 -0.000 3/14112 990.0 0.898 0.898 0.480 0.000 -0.000 -0.000 3/14112 1000.0 0.889 0.889 0.476 0.000 -0.000 -0.000 3/14112 Thermal conductivity related properties were written into "kappa-m14143.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-09 04:34:13]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|