----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 04:50:48]------------------------- 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 3.821510150938196 0.000000000000000 0.000000000000000 b -1.910755075469098 3.309524871532582 0.000000000000000 c 0.000000000000000 0.000000000000000 14.714207790000000 Atomic positions (fractional): *1 S 0.33333333333333 0.66666666666667 0.09803873074900 32.065 2 S 0.66666666666667 0.33333333333333 0.59803873074900 32.065 3 S 0.33333333333333 0.66666666666667 0.40196126925100 32.065 4 S 0.66666666666667 0.33333333333333 0.90196126925100 32.065 *5 In 0.00000000000000 0.00000000000000 0.50000000000000 114.818 6 In 0.00000000000000 0.00000000000000 0.00000000000000 114.818 *7 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 8 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.821510150938196 0.000000000000000 0.000000000000000 b -1.910755075469098 3.309524871532582 0.000000000000000 c 0.000000000000000 0.000000000000000 14.714207790000000 Atomic positions (fractional): *1 S 0.33333333333333 0.66666666666667 0.09803873074900 32.065 > 1 2 S 0.66666666666667 0.33333333333333 0.59803873074900 32.065 > 2 3 S 0.33333333333333 0.66666666666667 0.40196126925100 32.065 > 3 4 S 0.66666666666667 0.33333333333333 0.90196126925100 32.065 > 4 *5 In 0.00000000000000 0.00000000000000 0.50000000000000 114.818 > 5 6 In 0.00000000000000 0.00000000000000 0.00000000000000 114.818 > 6 *7 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 > 7 8 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 15.286040603752785 0.000000000000000 0.000000000000000 b -7.643020301876392 13.238099486130327 0.000000000000000 c 0.000000000000000 0.000000000000000 14.714207790000000 Atomic positions (fractional): *1 S 0.08333333333333 0.16666666666667 0.09803873074900 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.09803873074900 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.09803873074900 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.09803873074900 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.09803873074900 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.09803873074900 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.09803873074900 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.09803873074900 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.09803873074900 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.09803873074900 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.09803873074900 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.09803873074900 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.09803873074900 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.09803873074900 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.09803873074900 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.09803873074900 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.59803873074900 32.065 > 2 18 S 0.41666666666667 0.08333333333333 0.59803873074900 32.065 > 2 19 S 0.66666666666667 0.08333333333333 0.59803873074900 32.065 > 2 20 S 0.91666666666667 0.08333333333333 0.59803873074900 32.065 > 2 21 S 0.16666666666667 0.33333333333333 0.59803873074900 32.065 > 2 22 S 0.41666666666667 0.33333333333333 0.59803873074900 32.065 > 2 23 S 0.66666666666667 0.33333333333333 0.59803873074900 32.065 > 2 24 S 0.91666666666667 0.33333333333333 0.59803873074900 32.065 > 2 25 S 0.16666666666667 0.58333333333333 0.59803873074900 32.065 > 2 26 S 0.41666666666667 0.58333333333333 0.59803873074900 32.065 > 2 27 S 0.66666666666667 0.58333333333333 0.59803873074900 32.065 > 2 28 S 0.91666666666667 0.58333333333333 0.59803873074900 32.065 > 2 29 S 0.16666666666667 0.83333333333333 0.59803873074900 32.065 > 2 30 S 0.41666666666667 0.83333333333333 0.59803873074900 32.065 > 2 31 S 0.66666666666667 0.83333333333333 0.59803873074900 32.065 > 2 32 S 0.91666666666667 0.83333333333333 0.59803873074900 32.065 > 2 33 S 0.08333333333333 0.16666666666667 0.40196126925100 32.065 > 3 34 S 0.33333333333333 0.16666666666667 0.40196126925100 32.065 > 3 35 S 0.58333333333333 0.16666666666667 0.40196126925100 32.065 > 3 36 S 0.83333333333333 0.16666666666667 0.40196126925100 32.065 > 3 37 S 0.08333333333333 0.41666666666667 0.40196126925100 32.065 > 3 38 S 0.33333333333333 0.41666666666667 0.40196126925100 32.065 > 3 39 S 0.58333333333333 0.41666666666667 0.40196126925100 32.065 > 3 40 S 0.83333333333333 0.41666666666667 0.40196126925100 32.065 > 3 41 S 0.08333333333333 0.66666666666667 0.40196126925100 32.065 > 3 42 S 0.33333333333333 0.66666666666667 0.40196126925100 32.065 > 3 43 S 0.58333333333333 0.66666666666667 0.40196126925100 32.065 > 3 44 S 0.83333333333333 0.66666666666667 0.40196126925100 32.065 > 3 45 S 0.08333333333333 0.91666666666667 0.40196126925100 32.065 > 3 46 S 0.33333333333333 0.91666666666667 0.40196126925100 32.065 > 3 47 S 0.58333333333333 0.91666666666667 0.40196126925100 32.065 > 3 48 S 0.83333333333333 0.91666666666667 0.40196126925100 32.065 > 3 49 S 0.16666666666667 0.08333333333333 0.90196126925100 32.065 > 4 50 S 0.41666666666667 0.08333333333333 0.90196126925100 32.065 > 4 51 S 0.66666666666667 0.08333333333333 0.90196126925100 32.065 > 4 52 S 0.91666666666667 0.08333333333333 0.90196126925100 32.065 > 4 53 S 0.16666666666667 0.33333333333333 0.90196126925100 32.065 > 4 54 S 0.41666666666667 0.33333333333333 0.90196126925100 32.065 > 4 55 S 0.66666666666667 0.33333333333333 0.90196126925100 32.065 > 4 56 S 0.91666666666667 0.33333333333333 0.90196126925100 32.065 > 4 57 S 0.16666666666667 0.58333333333333 0.90196126925100 32.065 > 4 58 S 0.41666666666667 0.58333333333333 0.90196126925100 32.065 > 4 59 S 0.66666666666667 0.58333333333333 0.90196126925100 32.065 > 4 60 S 0.91666666666667 0.58333333333333 0.90196126925100 32.065 > 4 61 S 0.16666666666667 0.83333333333333 0.90196126925100 32.065 > 4 62 S 0.41666666666667 0.83333333333333 0.90196126925100 32.065 > 4 63 S 0.66666666666667 0.83333333333333 0.90196126925100 32.065 > 4 64 S 0.91666666666667 0.83333333333333 0.90196126925100 32.065 > 4 *65 In 0.00000000000000 0.00000000000000 0.50000000000000 114.818 > 5 66 In 0.25000000000000 0.00000000000000 0.50000000000000 114.818 > 5 67 In 0.50000000000000 0.00000000000000 0.50000000000000 114.818 > 5 68 In 0.75000000000000 0.00000000000000 0.50000000000000 114.818 > 5 69 In 0.00000000000000 0.25000000000000 0.50000000000000 114.818 > 5 70 In 0.25000000000000 0.25000000000000 0.50000000000000 114.818 > 5 71 In 0.50000000000000 0.25000000000000 0.50000000000000 114.818 > 5 72 In 0.75000000000000 0.25000000000000 0.50000000000000 114.818 > 5 73 In 0.00000000000000 0.50000000000000 0.50000000000000 114.818 > 5 74 In 0.25000000000000 0.50000000000000 0.50000000000000 114.818 > 5 75 In 0.50000000000000 0.50000000000000 0.50000000000000 114.818 > 5 76 In 0.75000000000000 0.50000000000000 0.50000000000000 114.818 > 5 77 In 0.00000000000000 0.75000000000000 0.50000000000000 114.818 > 5 78 In 0.25000000000000 0.75000000000000 0.50000000000000 114.818 > 5 79 In 0.50000000000000 0.75000000000000 0.50000000000000 114.818 > 5 80 In 0.75000000000000 0.75000000000000 0.50000000000000 114.818 > 5 81 In 0.00000000000000 0.00000000000000 0.00000000000000 114.818 > 6 82 In 0.25000000000000 0.00000000000000 0.00000000000000 114.818 > 6 83 In 0.50000000000000 0.00000000000000 0.00000000000000 114.818 > 6 84 In 0.75000000000000 0.00000000000000 0.00000000000000 114.818 > 6 85 In 0.00000000000000 0.25000000000000 0.00000000000000 114.818 > 6 86 In 0.25000000000000 0.25000000000000 0.00000000000000 114.818 > 6 87 In 0.50000000000000 0.25000000000000 0.00000000000000 114.818 > 6 88 In 0.75000000000000 0.25000000000000 0.00000000000000 114.818 > 6 89 In 0.00000000000000 0.50000000000000 0.00000000000000 114.818 > 6 90 In 0.25000000000000 0.50000000000000 0.00000000000000 114.818 > 6 91 In 0.50000000000000 0.50000000000000 0.00000000000000 114.818 > 6 92 In 0.75000000000000 0.50000000000000 0.00000000000000 114.818 > 6 93 In 0.00000000000000 0.75000000000000 0.00000000000000 114.818 > 6 94 In 0.25000000000000 0.75000000000000 0.00000000000000 114.818 > 6 95 In 0.50000000000000 0.75000000000000 0.00000000000000 114.818 > 6 96 In 0.75000000000000 0.75000000000000 0.00000000000000 114.818 > 6 *97 Tl 0.16666666666667 0.08333333333333 0.25000000000000 204.383 > 7 98 Tl 0.41666666666667 0.08333333333333 0.25000000000000 204.383 > 7 99 Tl 0.66666666666667 0.08333333333333 0.25000000000000 204.383 > 7 100 Tl 0.91666666666667 0.08333333333333 0.25000000000000 204.383 > 7 101 Tl 0.16666666666667 0.33333333333333 0.25000000000000 204.383 > 7 102 Tl 0.41666666666667 0.33333333333333 0.25000000000000 204.383 > 7 103 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 > 7 104 Tl 0.91666666666667 0.33333333333333 0.25000000000000 204.383 > 7 105 Tl 0.16666666666667 0.58333333333333 0.25000000000000 204.383 > 7 106 Tl 0.41666666666667 0.58333333333333 0.25000000000000 204.383 > 7 107 Tl 0.66666666666667 0.58333333333333 0.25000000000000 204.383 > 7 108 Tl 0.91666666666667 0.58333333333333 0.25000000000000 204.383 > 7 109 Tl 0.16666666666667 0.83333333333333 0.25000000000000 204.383 > 7 110 Tl 0.41666666666667 0.83333333333333 0.25000000000000 204.383 > 7 111 Tl 0.66666666666667 0.83333333333333 0.25000000000000 204.383 > 7 112 Tl 0.91666666666667 0.83333333333333 0.25000000000000 204.383 > 7 113 Tl 0.08333333333333 0.16666666666667 0.75000000000000 204.383 > 8 114 Tl 0.33333333333333 0.16666666666667 0.75000000000000 204.383 > 8 115 Tl 0.58333333333333 0.16666666666667 0.75000000000000 204.383 > 8 116 Tl 0.83333333333333 0.16666666666667 0.75000000000000 204.383 > 8 117 Tl 0.08333333333333 0.41666666666667 0.75000000000000 204.383 > 8 118 Tl 0.33333333333333 0.41666666666667 0.75000000000000 204.383 > 8 119 Tl 0.58333333333333 0.41666666666667 0.75000000000000 204.383 > 8 120 Tl 0.83333333333333 0.41666666666667 0.75000000000000 204.383 > 8 121 Tl 0.08333333333333 0.66666666666667 0.75000000000000 204.383 > 8 122 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 > 8 123 Tl 0.58333333333333 0.66666666666667 0.75000000000000 204.383 > 8 124 Tl 0.83333333333333 0.66666666666667 0.75000000000000 204.383 > 8 125 Tl 0.08333333333333 0.91666666666667 0.75000000000000 204.383 > 8 126 Tl 0.33333333333333 0.91666666666667 0.75000000000000 204.383 > 8 127 Tl 0.58333333333333 0.91666666666667 0.75000000000000 204.383 > 8 128 Tl 0.83333333333333 0.91666666666667 0.75000000000000 204.383 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.9878453 -0.0000000 0.0000000 0.0000000 9.9878453 0.0000000 0.0000000 0.0000000 19.1018958 -------------------------- Born effective charges -------------------------- 1 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 2 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 3 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 4 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 5 In 4.3954468 -0.0000000 0.0000000 0.0000000 4.3954468 0.0000000 0.0000000 0.0000000 4.2850499 6 In 4.3954468 -0.0000000 0.0000000 0.0000000 4.3954468 0.0000000 0.0000000 0.0000000 4.2850499 7 Tl 1.6709507 -0.0000000 0.0000000 0.0000000 1.6709507 0.0000000 0.0000000 0.0000000 3.7130069 8 Tl 1.6709507 -0.0000000 0.0000000 0.0000000 1.6709507 0.0000000 0.0000000 0.0000000 3.7130069 ---------------------------------------------------------------------------- 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: 182 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.171 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.182 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) Permutation basis: 384/384 Permutation basis: 5136/5136 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 259 Number of blocks in projector: 182 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-08 04:50:53]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 04:50:53]------------------------- 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 3.821510150938196 0.000000000000000 0.000000000000000 b -1.910755075469098 3.309524871532582 0.000000000000000 c 0.000000000000000 0.000000000000000 14.714207790000000 Atomic positions (fractional): 1 S 0.33333333333333 0.66666666666667 0.09803873074900 32.065 2 S 0.66666666666667 0.33333333333333 0.59803873074900 32.065 3 S 0.33333333333333 0.66666666666667 0.40196126925100 32.065 4 S 0.66666666666667 0.33333333333333 0.90196126925100 32.065 5 In 0.00000000000000 0.00000000000000 0.50000000000000 114.818 6 In 0.00000000000000 0.00000000000000 0.00000000000000 114.818 7 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 8 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.286040603752785 0.000000000000000 0.000000000000000 b -7.643020301876392 13.238099486130327 0.000000000000000 c 0.000000000000000 0.000000000000000 14.714207790000000 Atomic positions (fractional): 1 S 0.08333333333333 0.16666666666667 0.09803873074900 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.09803873074900 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.09803873074900 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.09803873074900 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.09803873074900 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.09803873074900 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.09803873074900 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.09803873074900 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.09803873074900 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.09803873074900 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.09803873074900 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.09803873074900 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.09803873074900 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.09803873074900 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.09803873074900 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.09803873074900 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.59803873074900 32.065 > 17 18 S 0.41666666666667 0.08333333333333 0.59803873074900 32.065 > 17 19 S 0.66666666666667 0.08333333333333 0.59803873074900 32.065 > 17 20 S 0.91666666666667 0.08333333333333 0.59803873074900 32.065 > 17 21 S 0.16666666666667 0.33333333333333 0.59803873074900 32.065 > 17 22 S 0.41666666666667 0.33333333333333 0.59803873074900 32.065 > 17 23 S 0.66666666666667 0.33333333333333 0.59803873074900 32.065 > 17 24 S 0.91666666666667 0.33333333333333 0.59803873074900 32.065 > 17 25 S 0.16666666666667 0.58333333333333 0.59803873074900 32.065 > 17 26 S 0.41666666666667 0.58333333333333 0.59803873074900 32.065 > 17 27 S 0.66666666666667 0.58333333333333 0.59803873074900 32.065 > 17 28 S 0.91666666666667 0.58333333333333 0.59803873074900 32.065 > 17 29 S 0.16666666666667 0.83333333333333 0.59803873074900 32.065 > 17 30 S 0.41666666666667 0.83333333333333 0.59803873074900 32.065 > 17 31 S 0.66666666666667 0.83333333333333 0.59803873074900 32.065 > 17 32 S 0.91666666666667 0.83333333333333 0.59803873074900 32.065 > 17 33 S 0.08333333333333 0.16666666666667 0.40196126925100 32.065 > 33 34 S 0.33333333333333 0.16666666666667 0.40196126925100 32.065 > 33 35 S 0.58333333333333 0.16666666666667 0.40196126925100 32.065 > 33 36 S 0.83333333333333 0.16666666666667 0.40196126925100 32.065 > 33 37 S 0.08333333333333 0.41666666666667 0.40196126925100 32.065 > 33 38 S 0.33333333333333 0.41666666666667 0.40196126925100 32.065 > 33 39 S 0.58333333333333 0.41666666666667 0.40196126925100 32.065 > 33 40 S 0.83333333333333 0.41666666666667 0.40196126925100 32.065 > 33 41 S 0.08333333333333 0.66666666666667 0.40196126925100 32.065 > 33 42 S 0.33333333333333 0.66666666666667 0.40196126925100 32.065 > 33 43 S 0.58333333333333 0.66666666666667 0.40196126925100 32.065 > 33 44 S 0.83333333333333 0.66666666666667 0.40196126925100 32.065 > 33 45 S 0.08333333333333 0.91666666666667 0.40196126925100 32.065 > 33 46 S 0.33333333333333 0.91666666666667 0.40196126925100 32.065 > 33 47 S 0.58333333333333 0.91666666666667 0.40196126925100 32.065 > 33 48 S 0.83333333333333 0.91666666666667 0.40196126925100 32.065 > 33 49 S 0.16666666666667 0.08333333333333 0.90196126925100 32.065 > 49 50 S 0.41666666666667 0.08333333333333 0.90196126925100 32.065 > 49 51 S 0.66666666666667 0.08333333333333 0.90196126925100 32.065 > 49 52 S 0.91666666666667 0.08333333333333 0.90196126925100 32.065 > 49 53 S 0.16666666666667 0.33333333333333 0.90196126925100 32.065 > 49 54 S 0.41666666666667 0.33333333333333 0.90196126925100 32.065 > 49 55 S 0.66666666666667 0.33333333333333 0.90196126925100 32.065 > 49 56 S 0.91666666666667 0.33333333333333 0.90196126925100 32.065 > 49 57 S 0.16666666666667 0.58333333333333 0.90196126925100 32.065 > 49 58 S 0.41666666666667 0.58333333333333 0.90196126925100 32.065 > 49 59 S 0.66666666666667 0.58333333333333 0.90196126925100 32.065 > 49 60 S 0.91666666666667 0.58333333333333 0.90196126925100 32.065 > 49 61 S 0.16666666666667 0.83333333333333 0.90196126925100 32.065 > 49 62 S 0.41666666666667 0.83333333333333 0.90196126925100 32.065 > 49 63 S 0.66666666666667 0.83333333333333 0.90196126925100 32.065 > 49 64 S 0.91666666666667 0.83333333333333 0.90196126925100 32.065 > 49 65 In 0.00000000000000 0.00000000000000 0.50000000000000 114.818 > 65 66 In 0.25000000000000 0.00000000000000 0.50000000000000 114.818 > 65 67 In 0.50000000000000 0.00000000000000 0.50000000000000 114.818 > 65 68 In 0.75000000000000 0.00000000000000 0.50000000000000 114.818 > 65 69 In 0.00000000000000 0.25000000000000 0.50000000000000 114.818 > 65 70 In 0.25000000000000 0.25000000000000 0.50000000000000 114.818 > 65 71 In 0.50000000000000 0.25000000000000 0.50000000000000 114.818 > 65 72 In 0.75000000000000 0.25000000000000 0.50000000000000 114.818 > 65 73 In 0.00000000000000 0.50000000000000 0.50000000000000 114.818 > 65 74 In 0.25000000000000 0.50000000000000 0.50000000000000 114.818 > 65 75 In 0.50000000000000 0.50000000000000 0.50000000000000 114.818 > 65 76 In 0.75000000000000 0.50000000000000 0.50000000000000 114.818 > 65 77 In 0.00000000000000 0.75000000000000 0.50000000000000 114.818 > 65 78 In 0.25000000000000 0.75000000000000 0.50000000000000 114.818 > 65 79 In 0.50000000000000 0.75000000000000 0.50000000000000 114.818 > 65 80 In 0.75000000000000 0.75000000000000 0.50000000000000 114.818 > 65 81 In 0.00000000000000 0.00000000000000 0.00000000000000 114.818 > 81 82 In 0.25000000000000 0.00000000000000 0.00000000000000 114.818 > 81 83 In 0.50000000000000 0.00000000000000 0.00000000000000 114.818 > 81 84 In 0.75000000000000 0.00000000000000 0.00000000000000 114.818 > 81 85 In 0.00000000000000 0.25000000000000 0.00000000000000 114.818 > 81 86 In 0.25000000000000 0.25000000000000 0.00000000000000 114.818 > 81 87 In 0.50000000000000 0.25000000000000 0.00000000000000 114.818 > 81 88 In 0.75000000000000 0.25000000000000 0.00000000000000 114.818 > 81 89 In 0.00000000000000 0.50000000000000 0.00000000000000 114.818 > 81 90 In 0.25000000000000 0.50000000000000 0.00000000000000 114.818 > 81 91 In 0.50000000000000 0.50000000000000 0.00000000000000 114.818 > 81 92 In 0.75000000000000 0.50000000000000 0.00000000000000 114.818 > 81 93 In 0.00000000000000 0.75000000000000 0.00000000000000 114.818 > 81 94 In 0.25000000000000 0.75000000000000 0.00000000000000 114.818 > 81 95 In 0.50000000000000 0.75000000000000 0.00000000000000 114.818 > 81 96 In 0.75000000000000 0.75000000000000 0.00000000000000 114.818 > 81 97 Tl 0.16666666666667 0.08333333333333 0.25000000000000 204.383 > 97 98 Tl 0.41666666666667 0.08333333333333 0.25000000000000 204.383 > 97 99 Tl 0.66666666666667 0.08333333333333 0.25000000000000 204.383 > 97 100 Tl 0.91666666666667 0.08333333333333 0.25000000000000 204.383 > 97 101 Tl 0.16666666666667 0.33333333333333 0.25000000000000 204.383 > 97 102 Tl 0.41666666666667 0.33333333333333 0.25000000000000 204.383 > 97 103 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 > 97 104 Tl 0.91666666666667 0.33333333333333 0.25000000000000 204.383 > 97 105 Tl 0.16666666666667 0.58333333333333 0.25000000000000 204.383 > 97 106 Tl 0.41666666666667 0.58333333333333 0.25000000000000 204.383 > 97 107 Tl 0.66666666666667 0.58333333333333 0.25000000000000 204.383 > 97 108 Tl 0.91666666666667 0.58333333333333 0.25000000000000 204.383 > 97 109 Tl 0.16666666666667 0.83333333333333 0.25000000000000 204.383 > 97 110 Tl 0.41666666666667 0.83333333333333 0.25000000000000 204.383 > 97 111 Tl 0.66666666666667 0.83333333333333 0.25000000000000 204.383 > 97 112 Tl 0.91666666666667 0.83333333333333 0.25000000000000 204.383 > 97 113 Tl 0.08333333333333 0.16666666666667 0.75000000000000 204.383 > 113 114 Tl 0.33333333333333 0.16666666666667 0.75000000000000 204.383 > 113 115 Tl 0.58333333333333 0.16666666666667 0.75000000000000 204.383 > 113 116 Tl 0.83333333333333 0.16666666666667 0.75000000000000 204.383 > 113 117 Tl 0.08333333333333 0.41666666666667 0.75000000000000 204.383 > 113 118 Tl 0.33333333333333 0.41666666666667 0.75000000000000 204.383 > 113 119 Tl 0.58333333333333 0.41666666666667 0.75000000000000 204.383 > 113 120 Tl 0.83333333333333 0.41666666666667 0.75000000000000 204.383 > 113 121 Tl 0.08333333333333 0.66666666666667 0.75000000000000 204.383 > 113 122 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 > 113 123 Tl 0.58333333333333 0.66666666666667 0.75000000000000 204.383 > 113 124 Tl 0.83333333333333 0.66666666666667 0.75000000000000 204.383 > 113 125 Tl 0.08333333333333 0.91666666666667 0.75000000000000 204.383 > 113 126 Tl 0.33333333333333 0.91666666666667 0.75000000000000 204.383 > 113 127 Tl 0.58333333333333 0.91666666666667 0.75000000000000 204.383 > 113 128 Tl 0.83333333333333 0.91666666666667 0.75000000000000 204.383 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.9878453 -0.0000000 0.0000000 0.0000000 9.9878453 0.0000000 0.0000000 0.0000000 19.1018958 -------------------------- Born effective charges -------------------------- 1 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 2 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 3 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 4 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 5 In 4.3954468 -0.0000000 0.0000000 0.0000000 4.3954468 0.0000000 0.0000000 0.0000000 4.2850499 6 In 4.3954468 -0.0000000 0.0000000 0.0000000 4.3954468 0.0000000 0.0000000 0.0000000 4.2850499 7 Tl 1.6709507 -0.0000000 0.0000000 0.0000000 1.6709507 0.0000000 0.0000000 0.0000000 3.7130069 8 Tl 1.6709507 -0.0000000 0.0000000 0.0000000 1.6709507 0.0000000 0.0000000 0.0000000 3.7130069 ---------------------------------------------------------------------------- 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[ 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] Computing fc3[ 97, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.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.00000007 (xxy) -0.00000007 (xxy) -0.00000007 (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-08 04:51:00]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 04:51:00]------------------------- 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 3.821510150938196 0.000000000000000 0.000000000000000 b -1.910755075469098 3.309524871532582 0.000000000000000 c 0.000000000000000 0.000000000000000 14.714207790000000 Atomic positions (fractional): 1 S 0.33333333333333 0.66666666666667 0.09803873074900 32.065 2 S 0.66666666666667 0.33333333333333 0.59803873074900 32.065 3 S 0.33333333333333 0.66666666666667 0.40196126925100 32.065 4 S 0.66666666666667 0.33333333333333 0.90196126925100 32.065 5 In 0.00000000000000 0.00000000000000 0.50000000000000 114.818 6 In 0.00000000000000 0.00000000000000 0.00000000000000 114.818 7 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 8 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.286040603752785 0.000000000000000 0.000000000000000 b -7.643020301876392 13.238099486130327 0.000000000000000 c 0.000000000000000 0.000000000000000 14.714207790000000 Atomic positions (fractional): 1 S 0.08333333333333 0.16666666666667 0.09803873074900 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.09803873074900 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.09803873074900 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.09803873074900 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.09803873074900 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.09803873074900 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.09803873074900 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.09803873074900 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.09803873074900 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.09803873074900 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.09803873074900 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.09803873074900 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.09803873074900 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.09803873074900 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.09803873074900 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.09803873074900 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.59803873074900 32.065 > 17 18 S 0.41666666666667 0.08333333333333 0.59803873074900 32.065 > 17 19 S 0.66666666666667 0.08333333333333 0.59803873074900 32.065 > 17 20 S 0.91666666666667 0.08333333333333 0.59803873074900 32.065 > 17 21 S 0.16666666666667 0.33333333333333 0.59803873074900 32.065 > 17 22 S 0.41666666666667 0.33333333333333 0.59803873074900 32.065 > 17 23 S 0.66666666666667 0.33333333333333 0.59803873074900 32.065 > 17 24 S 0.91666666666667 0.33333333333333 0.59803873074900 32.065 > 17 25 S 0.16666666666667 0.58333333333333 0.59803873074900 32.065 > 17 26 S 0.41666666666667 0.58333333333333 0.59803873074900 32.065 > 17 27 S 0.66666666666667 0.58333333333333 0.59803873074900 32.065 > 17 28 S 0.91666666666667 0.58333333333333 0.59803873074900 32.065 > 17 29 S 0.16666666666667 0.83333333333333 0.59803873074900 32.065 > 17 30 S 0.41666666666667 0.83333333333333 0.59803873074900 32.065 > 17 31 S 0.66666666666667 0.83333333333333 0.59803873074900 32.065 > 17 32 S 0.91666666666667 0.83333333333333 0.59803873074900 32.065 > 17 33 S 0.08333333333333 0.16666666666667 0.40196126925100 32.065 > 33 34 S 0.33333333333333 0.16666666666667 0.40196126925100 32.065 > 33 35 S 0.58333333333333 0.16666666666667 0.40196126925100 32.065 > 33 36 S 0.83333333333333 0.16666666666667 0.40196126925100 32.065 > 33 37 S 0.08333333333333 0.41666666666667 0.40196126925100 32.065 > 33 38 S 0.33333333333333 0.41666666666667 0.40196126925100 32.065 > 33 39 S 0.58333333333333 0.41666666666667 0.40196126925100 32.065 > 33 40 S 0.83333333333333 0.41666666666667 0.40196126925100 32.065 > 33 41 S 0.08333333333333 0.66666666666667 0.40196126925100 32.065 > 33 42 S 0.33333333333333 0.66666666666667 0.40196126925100 32.065 > 33 43 S 0.58333333333333 0.66666666666667 0.40196126925100 32.065 > 33 44 S 0.83333333333333 0.66666666666667 0.40196126925100 32.065 > 33 45 S 0.08333333333333 0.91666666666667 0.40196126925100 32.065 > 33 46 S 0.33333333333333 0.91666666666667 0.40196126925100 32.065 > 33 47 S 0.58333333333333 0.91666666666667 0.40196126925100 32.065 > 33 48 S 0.83333333333333 0.91666666666667 0.40196126925100 32.065 > 33 49 S 0.16666666666667 0.08333333333333 0.90196126925100 32.065 > 49 50 S 0.41666666666667 0.08333333333333 0.90196126925100 32.065 > 49 51 S 0.66666666666667 0.08333333333333 0.90196126925100 32.065 > 49 52 S 0.91666666666667 0.08333333333333 0.90196126925100 32.065 > 49 53 S 0.16666666666667 0.33333333333333 0.90196126925100 32.065 > 49 54 S 0.41666666666667 0.33333333333333 0.90196126925100 32.065 > 49 55 S 0.66666666666667 0.33333333333333 0.90196126925100 32.065 > 49 56 S 0.91666666666667 0.33333333333333 0.90196126925100 32.065 > 49 57 S 0.16666666666667 0.58333333333333 0.90196126925100 32.065 > 49 58 S 0.41666666666667 0.58333333333333 0.90196126925100 32.065 > 49 59 S 0.66666666666667 0.58333333333333 0.90196126925100 32.065 > 49 60 S 0.91666666666667 0.58333333333333 0.90196126925100 32.065 > 49 61 S 0.16666666666667 0.83333333333333 0.90196126925100 32.065 > 49 62 S 0.41666666666667 0.83333333333333 0.90196126925100 32.065 > 49 63 S 0.66666666666667 0.83333333333333 0.90196126925100 32.065 > 49 64 S 0.91666666666667 0.83333333333333 0.90196126925100 32.065 > 49 65 In 0.00000000000000 0.00000000000000 0.50000000000000 114.818 > 65 66 In 0.25000000000000 0.00000000000000 0.50000000000000 114.818 > 65 67 In 0.50000000000000 0.00000000000000 0.50000000000000 114.818 > 65 68 In 0.75000000000000 0.00000000000000 0.50000000000000 114.818 > 65 69 In 0.00000000000000 0.25000000000000 0.50000000000000 114.818 > 65 70 In 0.25000000000000 0.25000000000000 0.50000000000000 114.818 > 65 71 In 0.50000000000000 0.25000000000000 0.50000000000000 114.818 > 65 72 In 0.75000000000000 0.25000000000000 0.50000000000000 114.818 > 65 73 In 0.00000000000000 0.50000000000000 0.50000000000000 114.818 > 65 74 In 0.25000000000000 0.50000000000000 0.50000000000000 114.818 > 65 75 In 0.50000000000000 0.50000000000000 0.50000000000000 114.818 > 65 76 In 0.75000000000000 0.50000000000000 0.50000000000000 114.818 > 65 77 In 0.00000000000000 0.75000000000000 0.50000000000000 114.818 > 65 78 In 0.25000000000000 0.75000000000000 0.50000000000000 114.818 > 65 79 In 0.50000000000000 0.75000000000000 0.50000000000000 114.818 > 65 80 In 0.75000000000000 0.75000000000000 0.50000000000000 114.818 > 65 81 In 0.00000000000000 0.00000000000000 0.00000000000000 114.818 > 81 82 In 0.25000000000000 0.00000000000000 0.00000000000000 114.818 > 81 83 In 0.50000000000000 0.00000000000000 0.00000000000000 114.818 > 81 84 In 0.75000000000000 0.00000000000000 0.00000000000000 114.818 > 81 85 In 0.00000000000000 0.25000000000000 0.00000000000000 114.818 > 81 86 In 0.25000000000000 0.25000000000000 0.00000000000000 114.818 > 81 87 In 0.50000000000000 0.25000000000000 0.00000000000000 114.818 > 81 88 In 0.75000000000000 0.25000000000000 0.00000000000000 114.818 > 81 89 In 0.00000000000000 0.50000000000000 0.00000000000000 114.818 > 81 90 In 0.25000000000000 0.50000000000000 0.00000000000000 114.818 > 81 91 In 0.50000000000000 0.50000000000000 0.00000000000000 114.818 > 81 92 In 0.75000000000000 0.50000000000000 0.00000000000000 114.818 > 81 93 In 0.00000000000000 0.75000000000000 0.00000000000000 114.818 > 81 94 In 0.25000000000000 0.75000000000000 0.00000000000000 114.818 > 81 95 In 0.50000000000000 0.75000000000000 0.00000000000000 114.818 > 81 96 In 0.75000000000000 0.75000000000000 0.00000000000000 114.818 > 81 97 Tl 0.16666666666667 0.08333333333333 0.25000000000000 204.383 > 97 98 Tl 0.41666666666667 0.08333333333333 0.25000000000000 204.383 > 97 99 Tl 0.66666666666667 0.08333333333333 0.25000000000000 204.383 > 97 100 Tl 0.91666666666667 0.08333333333333 0.25000000000000 204.383 > 97 101 Tl 0.16666666666667 0.33333333333333 0.25000000000000 204.383 > 97 102 Tl 0.41666666666667 0.33333333333333 0.25000000000000 204.383 > 97 103 Tl 0.66666666666667 0.33333333333333 0.25000000000000 204.383 > 97 104 Tl 0.91666666666667 0.33333333333333 0.25000000000000 204.383 > 97 105 Tl 0.16666666666667 0.58333333333333 0.25000000000000 204.383 > 97 106 Tl 0.41666666666667 0.58333333333333 0.25000000000000 204.383 > 97 107 Tl 0.66666666666667 0.58333333333333 0.25000000000000 204.383 > 97 108 Tl 0.91666666666667 0.58333333333333 0.25000000000000 204.383 > 97 109 Tl 0.16666666666667 0.83333333333333 0.25000000000000 204.383 > 97 110 Tl 0.41666666666667 0.83333333333333 0.25000000000000 204.383 > 97 111 Tl 0.66666666666667 0.83333333333333 0.25000000000000 204.383 > 97 112 Tl 0.91666666666667 0.83333333333333 0.25000000000000 204.383 > 97 113 Tl 0.08333333333333 0.16666666666667 0.75000000000000 204.383 > 113 114 Tl 0.33333333333333 0.16666666666667 0.75000000000000 204.383 > 113 115 Tl 0.58333333333333 0.16666666666667 0.75000000000000 204.383 > 113 116 Tl 0.83333333333333 0.16666666666667 0.75000000000000 204.383 > 113 117 Tl 0.08333333333333 0.41666666666667 0.75000000000000 204.383 > 113 118 Tl 0.33333333333333 0.41666666666667 0.75000000000000 204.383 > 113 119 Tl 0.58333333333333 0.41666666666667 0.75000000000000 204.383 > 113 120 Tl 0.83333333333333 0.41666666666667 0.75000000000000 204.383 > 113 121 Tl 0.08333333333333 0.66666666666667 0.75000000000000 204.383 > 113 122 Tl 0.33333333333333 0.66666666666667 0.75000000000000 204.383 > 113 123 Tl 0.58333333333333 0.66666666666667 0.75000000000000 204.383 > 113 124 Tl 0.83333333333333 0.66666666666667 0.75000000000000 204.383 > 113 125 Tl 0.08333333333333 0.91666666666667 0.75000000000000 204.383 > 113 126 Tl 0.33333333333333 0.91666666666667 0.75000000000000 204.383 > 113 127 Tl 0.58333333333333 0.91666666666667 0.75000000000000 204.383 > 113 128 Tl 0.83333333333333 0.91666666666667 0.75000000000000 204.383 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 9.9878453 -0.0000000 0.0000000 0.0000000 9.9878453 0.0000000 0.0000000 0.0000000 19.1018958 -------------------------- Born effective charges -------------------------- 1 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 2 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 3 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 4 S -3.0331988 0.0000000 0.0000000 0.0000000 -3.0331988 0.0000000 0.0000000 0.0000000 -3.9990284 5 In 4.3954468 -0.0000000 0.0000000 0.0000000 4.3954468 0.0000000 0.0000000 0.0000000 4.2850499 6 In 4.3954468 -0.0000000 0.0000000 0.0000000 4.3954468 0.0000000 0.0000000 0.0000000 4.2850499 7 Tl 1.6709507 -0.0000000 0.0000000 0.0000000 1.6709507 0.0000000 0.0000000 0.0000000 3.7130069 8 Tl 1.6709507 -0.0000000 0.0000000 0.0000000 1.6709507 0.0000000 0.0000000 0.0000000 3.7130069 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000007 (xxy) -0.00000007 (xxy) -0.00000007 (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 ... [ 15 15 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.73, Number of G-points: 291, Lambda: 0.27 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/54) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 54 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.014 ( 0.000 0.000 0.000) 0.000 1.014 ( 0.000 0.000 0.000) 0.000 1.314 ( 0.000 0.000 0.000) 0.000 1.314 ( 0.000 0.000 0.000) 0.000 1.508 ( 0.000 0.000 0.000) 0.000 1.508 ( 0.000 0.000 0.000) 0.000 1.845 ( 0.000 0.000 0.000) 0.000 2.139 ( 0.000 0.000 0.000) 0.000 2.301 ( 0.000 0.000 0.000) 0.000 4.634 ( 0.000 0.000 0.000) 0.000 4.634 ( 0.000 0.000 0.000) 0.000 5.095 ( 0.000 0.000 0.000) 0.000 5.095 ( 0.000 0.000 0.000) 0.000 5.330 ( 0.000 0.000 0.000) 0.000 5.330 ( 0.000 0.000 0.000) 0.000 5.601 ( 0.000 0.000 0.000) 0.000 5.601 ( 0.000 0.000 0.000) 0.000 6.266 ( 0.000 0.000 0.000) 0.000 8.152 ( 0.000 0.000 0.000) 0.000 8.248 ( 0.000 0.000 0.000) 0.000 8.557 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/54) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.355 ( 14.581 8.418 0.000) 16.836 0.368 ( 14.944 8.628 0.000) 17.255 0.711 ( 26.196 15.125 0.000) 30.249 1.018 ( 0.293 0.169 0.000) 0.339 1.041 ( 1.893 1.093 0.000) 2.186 1.407 ( 7.424 4.286 0.000) 8.573 1.442 ( 7.381 4.262 0.000) 8.523 1.537 ( 2.598 1.500 0.000) 3.000 1.701 ( 9.927 5.732 0.000) 11.463 1.874 ( 2.503 1.445 0.000) 2.890 2.152 ( 3.402 1.964 0.000) 3.928 2.211 ( -6.012 -3.471 0.000) 6.942 4.653 ( 1.148 0.663 0.000) 1.325 4.897 ( 20.318 11.731 0.000) 23.462 5.069 ( -1.799 -1.039 0.000) 2.078 5.160 ( 5.739 3.314 0.000) 6.627 5.360 ( 2.813 1.624 0.000) 3.249 5.685 ( 6.639 3.833 0.000) 7.666 5.782 ( 14.179 8.186 0.000) 16.373 6.215 ( -4.143 -2.392 0.000) 4.783 7.991 ( -2.014 -1.163 0.000) 2.326 8.129 ( -1.974 -1.140 0.000) 2.279 8.223 ( -2.086 -1.204 0.000) 2.409 8.384 ( -13.677 -7.897 0.000) 15.793 ======================= Grid point 2 (3/54) ======================= q-point: ( 0.13 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.661 ( 11.338 6.546 0.000) 13.092 0.675 ( 11.111 6.415 0.000) 12.830 1.027 ( 0.465 0.268 0.000) 0.537 1.055 ( -0.076 -0.044 0.000) 0.087 1.152 ( 10.752 6.208 0.000) 12.415 1.479 ( -4.622 -2.669 0.000) 5.337 1.625 ( 10.599 6.119 0.000) 12.239 1.639 ( 6.359 3.671 0.000) 7.343 1.949 ( 3.881 2.241 0.000) 4.481 2.028 ( 7.044 4.067 0.000) 8.133 2.177 ( 14.443 8.339 0.000) 16.677 2.420 ( 19.028 10.986 0.000) 21.972 4.655 ( -1.608 -0.928 0.000) 1.856 5.033 ( -0.999 -0.577 0.000) 1.153 5.347 ( 9.457 5.460 0.000) 10.921 5.439 ( 23.120 13.348 0.000) 26.696 5.476 ( 7.363 4.251 0.000) 8.502 5.878 ( 9.330 5.386 0.000) 10.773 6.065 ( -3.855 -2.226 0.000) 4.452 6.214 ( 17.351 10.018 0.000) 20.036 7.774 ( -19.402 -11.202 0.000) 22.404 8.052 ( -5.310 -3.066 0.000) 6.131 8.120 ( -3.658 -2.112 0.000) 4.224 8.158 ( -3.321 -1.917 0.000) 3.834 ======================= Grid point 3 (4/54) ======================= q-point: ( 0.20 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.869 ( 6.398 3.694 0.000) 7.388 0.881 ( 6.692 3.863 0.000) 7.727 0.985 ( -5.130 -2.962 0.000) 5.924 1.039 ( 0.565 0.326 0.000) 0.653 1.308 ( -8.385 -4.841 0.000) 9.683 1.328 ( 5.301 3.060 0.000) 6.121 1.828 ( 9.469 5.467 0.000) 10.934 1.871 ( 10.071 5.815 0.000) 11.629 2.053 ( 4.692 2.709 0.000) 5.418 2.088 ( 3.778 2.181 0.000) 4.363 2.740 ( 26.275 15.170 0.000) 30.340 2.923 ( 22.455 12.964 0.000) 25.929 4.569 ( -5.505 -3.178 0.000) 6.357 5.026 ( 0.327 0.189 0.000) 0.378 5.551 ( 6.955 4.015 0.000) 8.031 5.691 ( 10.450 6.033 0.000) 12.066 5.802 ( 3.787 2.186 0.000) 4.372 5.950 ( -3.374 -1.948 0.000) 3.896 6.099 ( 9.363 5.406 0.000) 10.812 6.667 ( 19.253 11.115 0.000) 22.231 7.324 ( -13.478 -7.781 0.000) 15.563 7.911 ( -6.694 -3.865 0.000) 7.729 8.077 ( -3.452 -1.993 0.000) 3.986 8.093 ( -0.354 -0.204 0.000) 0.409 ======================= Grid point 4 (5/54) ======================= q-point: ( 0.27 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.862 ( -4.133 -2.386 0.000) 4.772 0.950 ( 0.312 0.180 0.000) 0.360 0.998 ( 3.508 2.025 0.000) 4.050 1.056 ( 1.112 0.642 0.000) 1.284 1.147 ( -4.412 -2.547 0.000) 5.095 1.424 ( 3.160 1.824 0.000) 3.648 2.053 ( 9.333 5.389 0.000) 10.777 2.081 ( 7.878 4.549 0.000) 9.097 2.143 ( 2.579 1.489 0.000) 2.978 2.204 ( 5.428 3.134 0.000) 6.268 3.321 ( 23.189 13.388 0.000) 26.776 3.420 ( 19.785 11.423 0.000) 22.845 4.423 ( -6.419 -3.706 0.000) 7.412 5.044 ( 1.052 0.608 0.000) 1.215 5.556 ( -21.044 -12.150 0.000) 24.299 5.641 ( 1.501 0.867 0.000) 1.733 5.930 ( 9.500 5.485 0.000) 10.970 5.953 ( 3.087 1.782 0.000) 3.565 6.302 ( 7.869 4.543 0.000) 9.087 7.078 ( 15.798 9.121 0.000) 18.242 7.231 ( 1.649 0.952 0.000) 1.904 7.761 ( -5.187 -2.995 0.000) 5.989 8.002 ( -2.997 -1.731 0.000) 3.461 8.085 ( -0.393 -0.227 0.000) 0.454 ======================= Grid point 5 (6/54) ======================= q-point: ( 0.33 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.857 ( 4.606 2.659 0.000) 5.318 0.911 ( -2.187 -1.263 0.000) 2.525 1.048 ( 1.256 0.725 0.000) 1.450 1.097 ( 2.061 1.190 0.000) 2.380 1.103 ( -0.521 -0.301 0.000) 0.602 1.477 ( 1.380 0.797 0.000) 1.593 2.162 ( -0.893 -0.515 0.000) 1.031 2.235 ( 5.379 3.106 0.000) 6.212 2.240 ( 6.443 3.720 0.000) 7.440 2.308 ( 2.884 1.665 0.000) 3.330 3.804 ( 18.118 10.461 0.000) 20.921 3.830 ( 15.272 8.817 0.000) 17.634 4.296 ( -4.131 -2.385 0.000) 4.770 5.025 ( -21.914 -12.652 0.000) 25.304 5.069 ( 1.048 0.605 0.000) 1.210 5.679 ( 2.394 1.382 0.000) 2.765 6.040 ( 3.674 2.121 0.000) 4.243 6.114 ( 6.110 3.528 0.000) 7.055 6.460 ( 5.658 3.267 0.000) 6.533 7.207 ( -4.828 -2.787 0.000) 5.575 7.393 ( 11.238 6.488 0.000) 12.976 7.717 ( 1.147 0.662 0.000) 1.324 7.937 ( -2.631 -1.519 0.000) 3.038 8.074 ( -0.498 -0.287 0.000) 0.575 ======================= Grid point 6 (7/54) ======================= q-point: ( 0.40 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.879 ( -0.631 -0.364 0.000) 0.729 1.065 ( 11.978 6.915 0.000) 13.831 1.073 ( 0.988 0.570 0.000) 1.141 1.093 ( -0.401 -0.232 0.000) 0.463 1.137 ( 1.283 0.741 0.000) 1.481 1.491 ( -0.151 -0.087 0.000) 0.174 2.116 ( -2.521 -1.455 0.000) 2.911 2.317 ( -2.210 -1.276 0.000) 2.552 2.333 ( 3.079 1.778 0.000) 3.555 2.349 ( 3.000 1.732 0.000) 3.465 4.124 ( 9.826 5.673 0.000) 11.346 4.153 ( 11.674 6.740 0.000) 13.480 4.237 ( -1.163 -0.671 0.000) 1.342 4.601 ( -13.874 -8.010 0.000) 16.021 5.089 ( 0.661 0.382 0.000) 0.763 5.750 ( 3.315 1.914 0.000) 3.828 6.110 ( 2.237 1.292 0.000) 2.584 6.215 ( 2.724 1.573 0.000) 3.145 6.564 ( 3.302 1.907 0.000) 3.813 7.060 ( -6.431 -3.713 0.000) 7.426 7.604 ( 7.064 4.079 0.000) 8.157 7.758 ( 1.536 0.887 0.000) 1.774 7.877 ( -2.545 -1.469 0.000) 2.938 8.063 ( -0.400 -0.231 0.000) 0.462 ======================= Grid point 7 (8/54) ======================= q-point: ( 0.47 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.873 ( -0.099 -0.057 0.000) 0.114 1.086 ( -0.179 -0.103 0.000) 0.207 1.091 ( 0.460 0.265 0.000) 0.531 1.156 ( 0.386 0.223 0.000) 0.446 1.333 ( 8.985 5.187 0.000) 10.375 1.471 ( -1.663 -0.960 0.000) 1.921 2.068 ( -1.262 -0.728 0.000) 1.457 2.233 ( -3.518 -2.031 0.000) 4.062 2.380 ( 0.999 0.577 0.000) 1.153 2.390 ( 0.745 0.430 0.000) 0.860 4.227 ( 0.011 0.007 0.000) 0.013 4.280 ( 3.442 1.987 0.000) 3.974 4.338 ( 4.050 2.339 0.000) 4.677 4.389 ( -4.456 -2.573 0.000) 5.145 5.100 ( 0.214 0.124 0.000) 0.247 5.812 ( 1.565 0.904 0.000) 1.807 6.143 ( 0.687 0.397 0.000) 0.793 6.251 ( 0.654 0.377 0.000) 0.755 6.614 ( 1.072 0.619 0.000) 1.237 6.949 ( -2.623 -1.515 0.000) 3.029 7.731 ( 4.007 2.314 0.000) 4.627 7.779 ( 0.384 0.222 0.000) 0.444 7.817 ( -2.674 -1.544 0.000) 3.088 8.057 ( -0.161 -0.093 0.000) 0.186 ======================= Grid point 16 (9/54) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.586 ( 7.114 12.322 0.000) 14.228 0.613 ( 7.464 12.928 0.000) 14.928 1.027 ( 2.399 4.154 0.000) 4.797 1.050 ( -0.246 -0.426 0.000) 0.492 1.067 ( 6.700 11.605 0.000) 13.400 1.496 ( -1.115 -1.931 0.000) 2.230 1.563 ( 5.969 10.339 0.000) 11.938 1.613 ( 3.607 6.248 0.000) 7.215 1.942 ( 3.059 5.299 0.000) 6.118 1.958 ( 11.672 20.217 0.000) 23.345 2.113 ( -2.514 -4.354 0.000) 5.028 2.298 ( 8.246 14.283 0.000) 16.493 4.776 ( 4.318 7.478 0.000) 8.635 5.088 ( 0.799 1.384 0.000) 1.598 5.137 ( 6.593 11.419 0.000) 13.186 5.216 ( 2.395 4.148 0.000) 4.790 5.476 ( 6.421 11.122 0.000) 12.842 5.886 ( 7.646 13.243 0.000) 15.292 6.029 ( 9.318 16.140 0.000) 18.637 6.130 ( -2.994 -5.186 0.000) 5.988 7.966 ( -1.618 -2.802 0.000) 3.236 8.075 ( -11.590 -20.074 0.000) 23.179 8.086 ( -1.669 -2.891 0.000) 3.338 8.178 ( -1.814 -3.141 0.000) 3.627 ======================= Grid point 17 (10/54) ======================= q-point: ( 0.13 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.806 ( 6.334 7.011 0.000) 9.448 0.836 ( 5.440 8.030 0.000) 9.699 0.997 ( -2.650 -4.626 0.000) 5.332 1.063 ( 0.257 2.194 0.000) 2.209 1.267 ( 5.348 4.831 0.000) 7.207 1.367 ( -6.016 -7.931 0.000) 9.955 1.785 ( 6.036 10.364 0.000) 11.993 1.788 ( 7.775 9.184 0.000) 12.033 2.068 ( 2.559 9.301 0.000) 9.646 2.075 ( 0.466 2.463 0.000) 2.506 2.477 ( 21.711 19.608 0.000) 29.255 2.699 ( 17.842 15.978 0.000) 23.950 4.828 ( -12.935 16.720 0.000) 21.139 5.089 ( -3.508 5.146 0.000) 6.228 5.365 ( 11.325 -1.250 0.000) 11.394 5.399 ( 20.794 -5.773 0.000) 21.580 5.728 ( 6.772 7.075 0.000) 9.794 5.982 ( -3.306 -4.073 0.000) 5.246 6.164 ( 1.222 19.204 0.000) 19.243 6.432 ( 19.163 10.767 0.000) 21.980 7.596 ( -21.073 -5.896 0.000) 21.882 7.920 ( -0.330 -10.223 0.000) 10.228 8.041 ( 1.006 -3.884 0.000) 4.012 8.104 ( -2.669 -3.292 0.000) 4.238 ======================= Grid point 18 (11/54) ======================= q-point: ( 0.20 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.870 ( -2.963 -6.282 0.000) 6.946 0.936 ( 2.720 2.745 0.000) 3.864 0.982 ( 2.890 4.812 0.000) 5.613 1.087 ( -0.819 3.627 0.000) 3.719 1.182 ( -5.240 -7.064 0.000) 8.796 1.378 ( 3.871 1.992 0.000) 4.354 1.994 ( 6.665 5.454 0.000) 8.612 2.007 ( 6.419 10.539 0.000) 12.340 2.163 ( 3.813 6.200 0.000) 7.278 2.205 ( 0.450 12.008 0.000) 12.016 3.053 ( 22.923 15.515 0.000) 27.680 3.179 ( 19.951 11.958 0.000) 23.260 4.702 ( -17.984 13.989 0.000) 22.784 5.088 ( -2.446 5.187 0.000) 5.735 5.526 ( 8.618 -3.492 0.000) 9.299 5.604 ( 0.488 -14.357 0.000) 14.365 5.809 ( 9.482 1.402 0.000) 9.585 5.921 ( 0.584 -0.676 0.000) 0.893 6.366 ( -0.060 18.050 0.000) 18.050 6.866 ( 18.492 8.551 0.000) 20.373 7.324 ( -9.546 5.546 0.000) 11.040 7.787 ( -3.030 -8.142 0.000) 8.687 8.025 ( -2.542 -3.096 0.000) 4.006 8.033 ( 2.707 -4.721 0.000) 5.442 ======================= Grid point 19 (12/54) ======================= q-point: ( 0.27 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.809 ( 2.005 -0.866 0.000) 2.184 0.944 ( -2.440 -0.605 0.000) 2.514 1.061 ( 1.367 3.466 0.000) 3.726 1.086 ( 0.162 -3.181 0.000) 3.185 1.110 ( -0.987 4.239 0.000) 4.352 1.447 ( 2.557 0.574 0.000) 2.620 2.112 ( 3.252 0.457 0.000) 3.284 2.188 ( 5.198 4.161 0.000) 6.658 2.308 ( 2.629 10.146 0.000) 10.481 2.322 ( -1.661 13.487 0.000) 13.589 3.568 ( 20.497 10.868 0.000) 23.200 3.612 ( 18.083 7.301 0.000) 19.501 4.516 ( -15.731 8.981 0.000) 18.114 5.105 ( -1.491 4.836 0.000) 5.061 5.234 ( -17.598 -14.411 0.000) 22.746 5.599 ( 4.863 -3.285 0.000) 5.869 5.986 ( 4.625 2.384 0.000) 5.204 6.063 ( 9.470 7.485 0.000) 12.071 6.523 ( -1.107 15.069 0.000) 15.110 7.231 ( 14.616 6.054 0.000) 15.820 7.247 ( -5.096 -2.165 0.000) 5.537 7.704 ( 0.121 -1.581 0.000) 1.586 7.956 ( -2.059 -2.830 0.000) 3.500 8.036 ( 2.569 -3.996 0.000) 4.751 ======================= Grid point 20 (13/54) ======================= q-point: ( 0.33 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.899 ( -1.664 0.160 0.000) 1.672 0.946 ( 9.354 6.916 0.000) 11.633 1.074 ( 0.894 -2.321 0.000) 2.487 1.103 ( -0.164 3.227 0.000) 3.231 1.138 ( -0.323 3.406 0.000) 3.421 1.478 ( 1.019 -0.507 0.000) 1.138 2.125 ( -0.483 -1.543 0.000) 1.616 2.263 ( 2.400 -2.661 0.000) 3.584 2.416 ( -2.950 12.788 0.000) 13.124 2.432 ( -1.334 11.284 0.000) 11.362 3.947 ( 14.370 2.720 0.000) 14.625 3.974 ( 15.874 6.188 0.000) 17.038 4.361 ( -8.946 4.291 0.000) 9.922 4.822 ( -18.432 -2.797 0.000) 18.643 5.130 ( -1.354 4.749 0.000) 4.938 5.657 ( 5.403 -2.946 0.000) 6.155 6.090 ( 3.097 3.722 0.000) 4.842 6.254 ( 2.072 10.219 0.000) 10.427 6.633 ( -2.084 12.388 0.000) 12.562 7.087 ( -5.083 -9.305 0.000) 10.603 7.495 ( 10.080 3.705 0.000) 10.739 7.730 ( 1.997 0.524 0.000) 2.064 7.895 ( -1.816 -2.576 0.000) 3.152 8.042 ( 1.768 -2.374 0.000) 2.960 ======================= Grid point 21 (14/54) ======================= q-point: ( 0.40 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.887 ( -0.821 1.056 0.000) 1.337 1.066 ( 0.904 -2.235 0.000) 2.411 1.128 ( -0.388 2.817 0.000) 2.844 1.144 ( 0.564 1.834 0.000) 1.919 1.228 ( 9.428 7.996 0.000) 12.362 1.476 ( -0.245 -1.178 0.000) 1.203 2.089 ( -1.713 -0.668 0.000) 1.838 2.226 ( -0.886 -5.908 0.000) 5.974 2.474 ( -4.302 11.293 0.000) 12.084 2.487 ( -4.567 10.536 0.000) 11.483 4.161 ( 9.424 -1.731 0.000) 9.581 4.233 ( 9.014 1.153 0.000) 9.088 4.311 ( -1.793 5.066 0.000) 5.374 4.541 ( -15.070 4.408 0.000) 15.702 5.151 ( -1.918 4.967 0.000) 5.324 5.717 ( 5.512 -4.667 0.000) 7.223 6.160 ( 0.220 4.228 0.000) 4.233 6.348 ( -3.437 11.423 0.000) 11.929 6.700 ( -3.431 10.513 0.000) 11.059 6.928 ( -1.478 -9.399 0.000) 9.515 7.663 ( 6.259 1.890 0.000) 6.538 7.760 ( 1.342 -0.534 0.000) 1.445 7.838 ( -1.957 -2.409 0.000) 3.104 8.047 ( 0.838 -1.069 0.000) 1.358 ======================= Grid point 22 (15/54) ======================= q-point: ( 0.47 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.885 ( -0.716 1.240 0.000) 1.431 1.063 ( 1.162 -2.013 0.000) 2.324 1.145 ( -2.631 4.557 0.000) 5.262 1.146 ( 0.799 -1.383 0.000) 1.597 1.397 ( -0.042 0.073 0.000) 0.084 1.449 ( 0.221 -0.383 0.000) 0.442 2.069 ( -0.480 0.831 0.000) 0.960 2.171 ( 1.877 -3.251 0.000) 3.754 2.493 ( -5.832 10.101 0.000) 11.664 2.498 ( -5.731 9.927 0.000) 11.462 4.235 ( 3.355 -5.812 0.000) 6.711 4.320 ( 1.244 -2.155 0.000) 2.489 4.374 ( -0.925 1.602 0.000) 1.850 4.380 ( -6.942 12.023 0.000) 13.884 5.160 ( -2.798 4.846 0.000) 5.596 5.743 ( 3.874 -6.709 0.000) 7.747 6.183 ( -2.150 3.725 0.000) 4.301 6.374 ( -6.457 11.184 0.000) 12.914 6.723 ( -5.297 9.175 0.000) 10.595 6.865 ( 3.927 -6.801 0.000) 7.854 7.749 ( 0.754 -1.306 0.000) 1.508 7.769 ( 0.625 -1.082 0.000) 1.249 7.786 ( 0.412 -0.713 0.000) 0.824 8.049 ( 0.382 -0.662 0.000) 0.764 ======================= Grid point 33 (16/54) ======================= q-point: ( 0.13 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.864 ( -3.656 -6.332 0.000) 7.311 0.924 ( 2.579 4.467 0.000) 5.158 0.983 ( 2.998 5.193 0.000) 5.997 1.116 ( 1.611 2.790 0.000) 3.221 1.192 ( -5.026 -8.705 0.000) 10.052 1.350 ( 2.073 3.590 0.000) 4.146 1.962 ( 4.337 7.512 0.000) 8.674 2.004 ( 5.491 9.510 0.000) 10.982 2.180 ( 5.005 8.669 0.000) 10.010 2.285 ( 6.441 11.156 0.000) 12.882 2.910 ( 12.971 22.467 0.000) 25.943 3.043 ( 10.272 17.791 0.000) 20.544 5.175 ( 1.145 1.983 0.000) 2.290 5.198 ( -4.118 -7.133 0.000) 8.236 5.262 ( 9.430 16.334 0.000) 18.861 5.387 ( 2.452 4.246 0.000) 4.903 5.727 ( -3.043 -5.271 0.000) 6.086 5.914 ( -0.945 -1.637 0.000) 1.890 6.558 ( 9.125 15.805 0.000) 18.250 6.654 ( 7.949 13.769 0.000) 15.899 7.548 ( -0.038 -0.067 0.000) 0.077 7.692 ( -6.224 -10.780 0.000) 12.448 7.972 ( -1.666 -2.886 0.000) 3.333 8.031 ( -2.185 -3.784 0.000) 4.370 ======================= Grid point 34 (17/54) ======================= q-point: ( 0.20 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.777 ( 0.185 -1.359 0.000) 1.372 0.977 ( -0.647 1.181 0.000) 1.347 1.050 ( -1.602 -5.483 0.000) 5.712 1.069 ( 1.737 3.595 0.000) 3.993 1.168 ( 0.287 4.158 0.000) 4.168 1.414 ( 2.136 1.615 0.000) 2.678 2.076 ( 2.665 2.766 0.000) 3.841 2.142 ( 3.587 3.176 0.000) 4.791 2.399 ( 4.218 15.426 0.000) 15.992 2.488 ( 1.474 14.834 0.000) 14.907 3.357 ( 16.286 14.528 0.000) 21.825 3.399 ( 13.988 9.970 0.000) 17.177 4.873 ( -17.493 -7.039 0.000) 18.856 5.187 ( -1.795 3.469 0.000) 3.906 5.438 ( -13.604 7.656 0.000) 15.611 5.486 ( 4.793 0.099 0.000) 4.794 5.830 ( 22.116 1.144 0.000) 22.146 5.944 ( 3.410 3.659 0.000) 5.002 6.746 ( -3.604 18.653 0.000) 18.998 7.010 ( 17.170 5.740 0.000) 18.104 7.387 ( -12.876 -4.611 0.000) 13.677 7.657 ( 4.575 -0.322 0.000) 4.586 7.934 ( 1.973 -4.728 0.000) 5.123 7.957 ( -1.988 -3.547 0.000) 4.066 ======================= Grid point 35 (18/54) ======================= q-point: ( 0.27 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.877 ( 7.049 7.640 0.000) 10.395 0.935 ( -2.211 -0.127 0.000) 2.215 1.028 ( 1.265 -2.529 0.000) 2.828 1.117 ( 0.633 1.849 0.000) 1.954 1.209 ( -0.953 5.145 0.000) 5.233 1.454 ( 1.338 0.217 0.000) 1.355 2.115 ( 0.539 -0.016 0.000) 0.539 2.197 ( 1.962 -2.017 0.000) 2.814 2.597 ( -0.302 16.462 0.000) 16.465 2.631 ( -3.005 16.318 0.000) 16.593 3.715 ( 14.234 3.235 0.000) 14.597 3.753 ( 15.993 7.369 0.000) 17.609 4.529 ( -10.982 -8.267 0.000) 13.746 5.197 ( -0.950 3.275 0.000) 3.410 5.225 ( -20.837 13.127 0.000) 24.627 5.548 ( 3.579 -1.722 0.000) 3.972 6.079 ( 4.473 7.306 0.000) 8.566 6.219 ( 15.067 7.476 0.000) 16.820 6.845 ( -4.806 16.116 0.000) 16.817 7.108 ( -9.759 -10.360 0.000) 14.233 7.325 ( 13.468 3.401 0.000) 13.891 7.718 ( 0.656 2.449 0.000) 2.536 7.893 ( -1.478 -3.282 0.000) 3.599 7.951 ( 4.651 -3.918 0.000) 6.081 ======================= Grid point 36 (19/54) ======================= q-point: ( 0.33 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.916 ( -0.973 1.311 0.000) 1.632 1.021 ( 1.086 -2.732 0.000) 2.940 1.089 ( 4.334 2.159 0.000) 4.842 1.160 ( 4.735 6.459 0.000) 8.008 1.240 ( -1.438 5.700 0.000) 5.879 1.465 ( 0.237 -0.720 0.000) 0.758 2.109 ( -0.882 0.016 0.000) 0.882 2.173 ( -0.021 -5.306 0.000) 5.306 2.715 ( -4.894 16.162 0.000) 16.887 2.716 ( -6.014 16.333 0.000) 17.405 3.950 ( 11.703 -2.048 0.000) 11.881 4.047 ( 12.516 0.969 0.000) 12.553 4.369 ( -2.305 0.036 0.000) 2.305 4.987 ( -22.373 14.571 0.000) 26.699 5.220 ( -0.326 3.320 0.000) 3.336 5.581 ( 3.766 -4.447 0.000) 5.828 6.217 ( 1.268 9.139 0.000) 9.226 6.485 ( 5.663 9.430 0.000) 11.000 6.862 ( -6.033 -10.044 0.000) 11.717 6.908 ( -5.382 14.073 0.000) 15.067 7.547 ( 9.195 1.588 0.000) 9.331 7.739 ( 0.417 0.583 0.000) 0.717 7.837 ( -1.282 -3.060 0.000) 3.318 7.996 ( 3.464 -1.950 0.000) 3.975 ======================= Grid point 37 (20/54) ======================= q-point: ( 0.40 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.916 ( -0.904 1.606 0.000) 1.843 1.015 ( 1.354 -2.542 0.000) 2.880 1.096 ( 1.629 -4.010 0.000) 4.329 1.248 ( -2.567 5.605 0.000) 6.165 1.360 ( 5.153 4.457 0.000) 6.813 1.451 ( -0.736 -1.197 0.000) 1.405 2.095 ( -1.184 1.202 0.000) 1.687 2.119 ( -0.029 -3.643 0.000) 3.643 2.754 ( -8.050 15.757 0.000) 17.694 2.761 ( -8.139 15.980 0.000) 17.933 4.075 ( 7.310 -6.262 0.000) 9.625 4.202 ( 7.072 -4.315 0.000) 8.285 4.417 ( 1.027 7.622 0.000) 7.691 4.759 ( -18.469 14.096 0.000) 23.233 5.244 ( -1.045 3.381 0.000) 3.539 5.591 ( 4.305 -7.198 0.000) 8.388 6.296 ( -2.980 9.219 0.000) 9.689 6.620 ( 1.934 5.643 0.000) 5.965 6.707 ( -4.014 -3.551 0.000) 5.359 6.941 ( -6.253 12.550 0.000) 14.021 7.683 ( 5.439 0.039 0.000) 5.439 7.749 ( 0.445 -0.229 0.000) 0.501 7.783 ( -1.445 -2.996 0.000) 3.326 8.024 ( 1.581 -1.340 0.000) 2.072 ======================= Grid point 49 (21/54) ======================= q-point: ( 0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.854 ( 5.020 8.696 0.000) 10.041 0.963 ( -1.255 -2.174 0.000) 2.511 0.994 ( -0.182 -0.315 0.000) 0.364 1.121 ( 0.873 1.512 0.000) 1.746 1.246 ( 1.977 3.425 0.000) 3.954 1.443 ( 0.683 1.183 0.000) 1.366 2.110 ( 0.460 0.797 0.000) 0.921 2.168 ( -0.272 -0.472 0.000) 0.545 2.701 ( 7.874 13.639 0.000) 15.748 2.762 ( 6.763 11.714 0.000) 13.526 3.587 ( 5.106 8.844 0.000) 10.212 3.633 ( 7.361 12.750 0.000) 14.722 4.580 ( -8.487 -14.700 0.000) 16.974 5.221 ( 0.298 0.515 0.000) 0.595 5.508 ( 0.624 1.080 0.000) 1.247 5.620 ( 0.203 0.351 0.000) 0.405 5.984 ( 9.594 16.618 0.000) 19.189 6.078 ( 5.288 9.159 0.000) 10.576 7.087 ( 6.004 10.399 0.000) 12.008 7.109 ( 5.039 8.728 0.000) 10.078 7.163 ( -7.811 -13.530 0.000) 15.623 7.747 ( 3.105 5.378 0.000) 6.210 7.856 ( -1.364 -2.363 0.000) 2.728 7.888 ( -1.880 -3.256 0.000) 3.759 ======================= Grid point 50 (22/54) ======================= q-point: ( 0.27 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.941 ( -0.592 0.612 0.000) 0.852 0.984 ( 0.276 -1.853 0.000) 1.874 1.050 ( 3.882 4.700 0.000) 6.096 1.150 ( 2.776 4.414 0.000) 5.215 1.304 ( 0.366 4.127 0.000) 4.143 1.458 ( 0.199 0.087 0.000) 0.217 2.119 ( 0.046 0.564 0.000) 0.565 2.132 ( -0.771 -3.678 0.000) 3.758 2.929 ( 1.040 15.789 0.000) 15.824 2.950 ( -0.470 14.692 0.000) 14.700 3.754 ( 8.149 0.905 0.000) 8.199 3.862 ( 9.215 3.449 0.000) 9.839 4.398 ( -1.860 -1.651 0.000) 2.487 5.231 ( 0.574 0.596 0.000) 0.827 5.509 ( -17.665 11.321 0.000) 20.982 5.511 ( 0.430 -2.469 0.000) 2.506 6.278 ( 4.508 11.733 0.000) 12.569 6.292 ( 13.375 -3.730 0.000) 13.885 6.946 ( -5.701 -1.558 0.000) 5.910 7.150 ( -5.797 13.717 0.000) 14.892 7.373 ( 12.454 1.605 0.000) 12.557 7.778 ( -2.968 3.063 0.000) 4.265 7.827 ( -1.378 -3.155 0.000) 3.443 7.891 ( 5.955 -1.850 0.000) 6.236 ======================= Grid point 51 (23/54) ======================= q-point: ( 0.33 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.943 ( -0.492 1.100 0.000) 1.205 0.973 ( 0.781 -1.856 0.000) 2.014 1.026 ( -0.432 -6.562 0.000) 6.576 1.310 ( 3.437 5.752 0.000) 6.700 1.347 ( 0.459 4.553 0.000) 4.576 1.452 ( -0.505 -0.577 0.000) 0.767 2.095 ( -0.035 -1.547 0.000) 1.547 2.126 ( -0.521 1.545 0.000) 1.631 3.048 ( -5.797 16.120 0.000) 17.131 3.057 ( -5.213 17.114 0.000) 17.891 3.878 ( 7.983 -4.839 0.000) 9.335 4.013 ( 8.414 -4.270 0.000) 9.435 4.476 ( 3.770 11.427 0.000) 12.033 5.256 ( 1.413 0.591 0.000) 1.532 5.259 ( -22.757 10.928 0.000) 25.245 5.478 ( 0.902 -5.599 0.000) 5.671 6.345 ( 9.269 -19.614 0.000) 21.694 6.442 ( -0.187 12.200 0.000) 12.202 6.951 ( -4.187 14.459 0.000) 15.053 7.176 ( -6.063 12.148 0.000) 13.577 7.564 ( 8.359 0.285 0.000) 8.364 7.763 ( -1.460 1.926 0.000) 2.417 7.772 ( -1.249 -3.469 0.000) 3.687 7.961 ( 3.839 -1.710 0.000) 4.203 ======================= Grid point 52 (24/54) ======================= q-point: ( 0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.943 ( -0.572 0.991 0.000) 1.145 0.972 ( 0.968 -1.677 0.000) 1.936 0.980 ( 3.932 -6.810 0.000) 7.864 1.345 ( -2.140 3.706 0.000) 4.280 1.408 ( -0.403 0.698 0.000) 0.806 1.433 ( 0.170 -0.294 0.000) 0.340 2.091 ( -0.196 0.339 0.000) 0.391 2.130 ( -1.114 1.930 0.000) 2.228 3.078 ( -9.040 15.658 0.000) 18.081 3.097 ( -9.518 16.486 0.000) 19.037 3.925 ( 4.822 -8.352 0.000) 9.644 4.066 ( 5.075 -8.790 0.000) 10.149 4.631 ( -7.027 12.171 0.000) 14.054 5.052 ( -8.804 15.249 0.000) 17.608 5.276 ( -0.017 0.029 0.000) 0.033 5.452 ( 3.493 -6.050 0.000) 6.986 6.317 ( 11.506 -19.928 0.000) 23.011 6.503 ( -5.921 10.255 0.000) 11.841 6.997 ( -8.892 15.401 0.000) 17.784 7.184 ( -6.516 11.287 0.000) 13.033 7.658 ( 1.652 -2.861 0.000) 3.303 7.724 ( 1.368 -2.370 0.000) 2.736 7.762 ( -0.937 1.623 0.000) 1.874 7.986 ( 1.490 -2.581 0.000) 2.981 ======================= Grid point 66 (25/54) ======================= q-point: ( 0.27 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.953 ( 0.240 0.416 0.000) 0.480 0.956 ( -0.480 -0.831 0.000) 0.960 0.995 ( -3.790 -6.565 0.000) 7.581 1.303 ( 3.873 6.708 0.000) 7.746 1.372 ( 1.487 2.575 0.000) 2.973 1.455 ( -0.244 -0.422 0.000) 0.487 2.087 ( -0.132 -0.228 0.000) 0.264 2.139 ( 0.734 1.271 0.000) 1.467 3.205 ( 5.926 10.264 0.000) 11.852 3.212 ( 6.826 11.823 0.000) 13.652 3.761 ( -0.012 -0.021 0.000) 0.025 3.895 ( 0.020 0.035 0.000) 0.040 4.502 ( 6.804 11.785 0.000) 13.608 5.240 ( 0.374 0.648 0.000) 0.748 5.444 ( -2.278 -3.945 0.000) 4.556 5.624 ( -0.624 -1.081 0.000) 1.248 6.058 ( -8.629 -14.946 0.000) 17.258 6.510 ( 5.967 10.335 0.000) 11.933 7.066 ( 5.667 9.815 0.000) 11.334 7.366 ( 2.522 4.369 0.000) 5.045 7.418 ( 3.592 6.221 0.000) 7.184 7.764 ( -1.858 -3.218 0.000) 3.716 7.817 ( 0.225 0.389 0.000) 0.450 7.887 ( 0.881 1.525 0.000) 1.761 ======================= Grid point 67 (26/54) ======================= q-point: ( 0.33 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.865 ( -0.596 -7.993 0.000) 8.016 0.949 ( 0.329 -0.420 0.000) 0.534 0.955 ( -0.196 0.070 0.000) 0.208 1.387 ( 0.727 2.310 0.000) 2.421 1.412 ( 0.856 1.844 0.000) 2.033 1.440 ( -0.912 -0.503 0.000) 1.041 2.106 ( 0.275 1.661 0.000) 1.684 2.161 ( -0.141 1.623 0.000) 1.629 3.350 ( -3.423 13.525 0.000) 13.952 3.385 ( -3.223 15.208 0.000) 15.546 3.758 ( 3.915 -7.095 0.000) 8.103 3.879 ( 4.023 -8.852 0.000) 9.723 4.821 ( 6.760 21.983 0.000) 22.999 5.262 ( 1.252 0.514 0.000) 1.354 5.368 ( -0.509 -4.974 0.000) 5.000 5.411 ( -21.214 4.135 0.000) 21.613 5.882 ( 11.509 -24.891 0.000) 27.423 6.662 ( -0.237 8.559 0.000) 8.563 7.228 ( -0.963 11.128 0.000) 11.169 7.386 ( -4.712 8.378 0.000) 9.612 7.566 ( 6.651 0.254 0.000) 6.656 7.694 ( -1.650 -4.325 0.000) 4.629 7.812 ( -1.972 2.704 0.000) 3.347 7.920 ( 2.684 -2.424 0.000) 3.617 ======================= Grid point 82 (27/54) ======================= q-point: ( 0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 54 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.772 ( -0.000 -0.000 0.000) 0.000 0.951 ( -0.174 -0.301 0.000) 0.348 0.951 ( 0.174 0.301 0.000) 0.348 1.417 ( -0.439 -0.760 0.000) 0.877 1.417 ( 0.438 0.759 0.000) 0.877 1.440 ( 0.000 0.000 0.000) 0.000 2.127 ( -0.000 -0.000 0.000) 0.000 2.179 ( -0.000 -0.000 0.000) 0.000 3.587 ( -5.014 -8.685 0.000) 10.028 3.587 ( 5.014 8.685 0.000) 10.028 3.664 ( -6.182 -10.707 0.000) 12.364 3.664 ( 6.182 10.707 0.000) 12.363 5.291 ( -1.274 -2.207 0.000) 2.548 5.291 ( 1.274 2.207 0.000) 2.548 5.303 ( -11.259 -19.502 0.000) 22.519 5.303 ( 11.259 19.501 0.000) 22.518 5.522 ( 0.000 0.000 0.000) 0.000 6.753 ( -0.000 -0.000 0.000) 0.000 7.348 ( -0.000 -0.000 0.000) 0.000 7.482 ( -0.000 -0.000 0.000) 0.000 7.604 ( -2.010 -3.481 0.000) 4.020 7.604 ( 2.010 3.481 0.000) 4.020 7.867 ( -1.359 -2.354 0.000) 2.718 7.867 ( 1.359 2.354 0.000) 2.718 ======================= Grid point 241 (28/54) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 54 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.438 ( -0.000 0.000 18.547) 18.547 0.438 ( -0.000 0.000 18.547) 18.547 0.817 ( -0.000 -0.000 -14.261) 14.261 0.817 ( 0.000 0.000 -14.261) 14.261 0.859 ( 0.000 -0.000 36.271) 36.271 1.399 ( -0.000 0.000 4.685) 4.685 1.399 ( -0.000 0.000 4.685) 4.685 1.481 ( 0.000 0.000 -2.435) 2.435 1.481 ( 0.000 0.000 -2.435) 2.435 1.614 ( 0.000 0.000 -29.913) 29.913 2.654 ( 0.000 -0.000 14.878) 14.878 2.909 ( 0.000 -0.000 -7.553) 7.553 4.725 ( 0.000 0.000 7.456) 7.456 4.725 ( -0.000 0.000 7.456) 7.456 4.935 ( 0.000 0.000 -10.075) 10.075 4.935 ( 0.000 0.000 -10.075) 10.075 5.445 ( 0.000 -0.000 6.379) 6.379 5.445 ( 0.000 0.000 6.379) 6.379 5.561 ( -0.000 0.000 -3.535) 3.535 5.561 ( 0.000 0.000 -3.535) 3.535 8.035 ( -0.000 0.000 3.691) 3.691 8.164 ( -0.000 0.000 -6.299) 6.299 8.290 ( 0.000 0.000 9.232) 9.232 8.482 ( -0.000 0.000 -6.502) 6.502 ======================= Grid point 242 (29/54) ======================= q-point: ( 0.07 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.563 ( 8.902 5.140 13.529) 16.991 0.602 ( 8.301 4.793 14.223) 17.152 0.866 ( 3.502 2.022 -11.632) 12.315 0.898 ( 4.712 2.720 -10.848) 12.136 1.018 ( 11.624 6.711 22.146) 25.896 1.416 ( 0.195 0.112 -6.564) 6.568 1.457 ( 5.061 2.922 3.150) 6.639 1.516 ( 3.148 1.818 -1.823) 4.067 1.647 ( 14.617 8.439 -0.066) 16.878 1.768 ( 14.114 8.149 -14.961) 22.123 2.497 ( -12.141 -7.010 10.110) 17.285 2.522 ( -19.896 -11.487 11.610) 25.740 4.735 ( 0.551 0.318 6.774) 6.804 4.926 ( -0.716 -0.414 -9.155) 9.192 4.958 ( 17.954 10.366 5.534) 21.457 5.116 ( 11.629 6.714 -6.205) 14.792 5.482 ( 3.427 1.979 7.460) 8.445 5.628 ( 5.547 3.203 -4.842) 8.029 5.730 ( 14.328 8.272 -5.098) 17.312 5.750 ( 15.448 8.919 -7.147) 19.216 8.014 ( -2.166 -1.251 3.639) 4.416 8.139 ( -2.576 -1.487 -6.126) 6.810 8.204 ( -5.918 -3.417 6.173) 9.209 8.335 ( -11.233 -6.485 -4.326) 13.673 ======================= Grid point 243 (30/54) ======================= q-point: ( 0.13 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.766 ( 6.227 3.595 7.898) 10.681 0.773 ( 8.134 4.696 7.465) 11.997 0.947 ( 3.052 1.762 -6.554) 7.442 0.962 ( 0.479 0.277 -7.812) 7.831 1.246 ( 6.695 3.866 7.154) 10.533 1.405 ( -1.261 -0.728 -5.867) 6.045 1.627 ( 9.134 5.273 0.217) 10.549 1.635 ( 7.099 4.098 -0.349) 8.204 2.073 ( 10.173 5.873 3.482) 12.252 2.123 ( 2.454 1.417 1.770) 3.341 2.219 ( 9.353 5.400 8.634) 13.827 2.381 ( 10.723 6.191 -3.897) 12.981 4.728 ( -1.502 -0.867 6.110) 6.352 4.904 ( -1.089 -0.629 -8.520) 8.612 5.387 ( 14.878 8.590 3.940) 17.626 5.477 ( 19.944 11.515 -0.796) 23.043 5.613 ( 7.776 4.490 9.095) 12.781 5.801 ( 8.866 5.119 -6.465) 12.108 5.959 ( 2.896 1.672 -1.534) 3.679 6.088 ( 14.419 8.325 -8.601) 18.739 7.847 ( -14.145 -8.167 5.381) 17.197 7.988 ( -8.183 -4.725 -6.491) 11.464 8.097 ( -2.344 -1.353 1.713) 3.203 8.145 ( -3.519 -2.031 -1.760) 4.428 ======================= Grid point 244 (31/54) ======================= q-point: ( 0.20 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.892 ( 4.376 2.526 1.604) 5.301 0.924 ( 4.952 2.859 3.390) 6.648 0.943 ( -1.367 -0.789 -2.777) 3.195 1.003 ( 1.898 1.096 -2.939) 3.667 1.320 ( -3.022 -1.745 1.290) 3.720 1.341 ( -1.107 -0.639 -0.026) 1.279 1.835 ( 9.654 5.574 0.648) 11.166 1.856 ( 9.964 5.753 -1.058) 11.554 2.085 ( 1.425 0.822 1.815) 2.450 2.099 ( 1.787 1.032 0.511) 2.125 2.773 ( 25.761 14.873 3.119) 29.909 2.868 ( 23.580 13.614 -4.294) 27.564 4.661 ( -4.028 -2.326 7.684) 8.982 4.879 ( -1.044 -0.603 -10.209) 10.279 5.638 ( 4.924 2.843 6.330) 8.508 5.817 ( 2.701 1.559 -8.626) 9.172 5.828 ( 10.104 5.834 9.334) 14.942 5.971 ( 3.529 2.037 14.456) 15.019 6.021 ( 9.631 5.561 -6.613) 12.939 6.385 ( 11.694 6.751 -19.061) 23.359 7.546 ( -7.779 -4.491 13.931) 16.576 7.822 ( -6.146 -3.548 -8.878) 11.366 8.066 ( -1.180 -0.681 0.404) 1.421 8.085 ( -1.779 -1.027 -0.717) 2.176 ======================= Grid point 245 (32/54) ======================= q-point: ( 0.27 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.923 ( 0.193 0.111 3.375) 3.382 0.958 ( 0.893 0.515 0.135) 1.040 1.012 ( 2.786 1.608 1.189) 3.430 1.042 ( 1.583 0.914 -1.156) 2.163 1.194 ( -4.645 -2.682 4.928) 7.284 1.346 ( 0.482 0.278 -6.452) 6.476 2.058 ( 8.999 5.195 0.407) 10.399 2.071 ( 8.274 4.777 -0.694) 9.580 2.155 ( 3.290 1.899 1.115) 3.959 2.187 ( 4.679 2.702 -1.370) 5.574 3.339 ( 22.495 12.988 1.700) 26.031 3.391 ( 20.867 12.048 -2.329) 24.208 4.558 ( -4.445 -2.566 11.052) 12.186 4.858 ( -0.693 -0.400 -13.439) 13.463 5.548 ( -12.349 -7.130 3.624) 14.713 5.683 ( -7.843 -4.528 -5.117) 10.401 6.056 ( 9.072 5.237 8.637) 13.576 6.114 ( 6.116 3.531 18.739) 20.026 6.233 ( 8.284 4.783 -5.903) 11.241 6.643 ( 9.195 5.309 -24.442) 26.649 7.516 ( 3.219 1.858 12.075) 12.634 7.715 ( -2.709 -1.564 -5.117) 5.997 8.026 ( -1.828 -1.055 1.774) 2.757 8.065 ( -0.484 -0.279 -1.459) 1.563 ======================= Grid point 246 (33/54) ======================= q-point: ( 0.33 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.939 ( -1.006 -0.581 2.171) 2.462 0.960 ( 1.939 1.119 2.815) 3.597 1.059 ( 1.423 0.822 0.897) 1.873 1.083 ( 1.822 1.052 -1.059) 2.356 1.150 ( 0.991 0.572 7.379) 7.468 1.374 ( 1.875 1.082 -8.678) 8.944 2.201 ( 0.297 0.171 3.039) 3.058 2.233 ( 5.674 3.276 -0.014) 6.552 2.237 ( 6.129 3.539 -0.255) 7.082 2.274 ( 2.134 1.232 -2.810) 3.738 3.807 ( 17.510 10.109 0.547) 20.226 3.825 ( 16.332 9.429 -0.666) 18.870 4.471 ( -2.753 -1.589 14.228) 14.578 4.848 ( -0.185 -0.107 -16.382) 16.384 5.189 ( -15.640 -9.030 13.594) 22.604 5.541 ( -4.585 -2.647 -13.858) 14.835 6.231 ( 3.930 2.269 18.864) 19.402 6.234 ( 6.043 3.489 8.164) 10.740 6.398 ( 5.810 3.355 -5.361) 8.588 6.773 ( 2.131 1.230 -25.425) 25.544 7.604 ( 3.235 1.868 7.518) 8.395 7.702 ( 1.119 0.646 -1.835) 2.244 7.992 ( -0.994 -0.574 3.624) 3.801 8.056 ( -0.225 -0.130 -1.738) 1.757 ======================= Grid point 247 (34/54) ======================= q-point: ( 0.40 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.921 ( -0.127 -0.073 3.402) 3.405 1.008 ( 1.531 0.884 -3.569) 3.983 1.087 ( 1.073 0.619 1.183) 1.713 1.119 ( 1.221 0.705 -1.363) 1.961 1.241 ( 6.450 3.724 9.353) 11.956 1.424 ( 2.183 1.260 -5.821) 6.343 2.171 ( -2.490 -1.438 4.204) 5.093 2.271 ( -2.338 -1.350 -3.822) 4.679 2.334 ( 3.057 1.765 0.186) 3.535 2.343 ( 3.018 1.742 -0.471) 3.516 4.135 ( 10.561 6.097 0.809) 12.222 4.150 ( 11.308 6.529 -0.394) 13.064 4.432 ( -0.722 -0.417 15.846) 15.867 4.848 ( 0.165 0.095 -17.947) 17.948 4.885 ( -9.947 -5.743 22.811) 25.540 5.464 ( -2.153 -1.243 -23.571) 23.701 6.302 ( 2.237 1.292 17.263) 17.456 6.337 ( 2.936 1.695 8.288) 8.955 6.502 ( 3.206 1.851 -5.358) 6.513 6.776 ( -0.972 -0.561 -20.966) 20.996 7.654 ( 1.160 0.670 3.910) 4.133 7.731 ( 0.973 0.561 -2.398) 2.648 7.980 ( -0.054 -0.031 5.480) 5.480 8.055 ( 0.113 0.065 -1.249) 1.256 ======================= Grid point 248 (35/54) ======================= q-point: ( 0.47 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.926 ( 0.242 0.140 4.195) 4.204 1.030 ( 0.467 0.269 -4.263) 4.297 1.106 ( 0.451 0.260 1.222) 1.328 1.138 ( 0.415 0.240 -1.367) 1.449 1.397 ( 5.342 3.084 3.962) 7.332 1.460 ( 0.415 0.239 -1.379) 1.460 2.114 ( -1.714 -0.990 3.539) 4.055 2.196 ( -2.810 -1.622 -3.051) 4.454 2.379 ( 0.935 0.540 0.061) 1.082 2.385 ( 0.808 0.467 -0.347) 0.995 4.301 ( 3.618 2.089 1.546) 4.454 4.330 ( 3.981 2.298 -0.839) 4.673 4.427 ( 0.035 0.020 16.184) 16.184 4.733 ( -3.174 -1.832 27.517) 27.760 4.852 ( 0.127 0.073 -18.372) 18.372 5.433 ( -0.605 -0.349 -29.106) 29.114 6.337 ( 0.780 0.451 16.091) 16.116 6.379 ( 0.791 0.456 8.612) 8.660 6.550 ( 0.990 0.571 -5.580) 5.696 6.755 ( -0.590 -0.341 -16.880) 16.894 7.665 ( 0.080 0.046 1.759) 1.761 7.744 ( 0.215 0.124 -3.512) 3.521 7.983 ( 0.120 0.069 6.545) 6.546 8.058 ( 0.103 0.060 -0.696) 0.706 ======================= Grid point 257 (36/54) ======================= q-point: ( 0.07 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.724 ( 3.742 6.482 9.832) 12.357 0.730 ( 5.503 9.531 8.860) 14.128 0.937 ( 2.403 4.161 -7.857) 9.210 0.949 ( 0.722 1.250 -8.365) 8.489 1.194 ( 4.637 8.031 9.706) 13.424 1.404 ( -0.779 -1.349 -7.498) 7.659 1.580 ( 5.162 8.940 1.074) 10.379 1.604 ( 4.015 6.954 -0.865) 8.077 1.974 ( 13.155 22.786 2.451) 26.425 2.051 ( 8.958 15.516 -3.932) 18.342 2.239 ( -6.256 -10.836 12.314) 17.555 2.342 ( -2.260 -3.914 1.261) 4.692 4.814 ( 2.801 4.852 3.897) 6.825 4.951 ( 1.118 1.937 -7.532) 7.857 5.200 ( 7.091 12.282 3.433) 14.591 5.255 ( 4.676 8.099 0.380) 9.360 5.609 ( 5.977 10.352 8.919) 14.914 5.797 ( 6.642 11.504 -7.161) 15.092 5.951 ( 4.715 8.167 -5.922) 11.135 5.967 ( 8.076 13.989 -4.316) 16.719 7.954 ( -3.103 -5.375 1.008) 6.288 8.019 ( -7.678 -13.299 -3.808) 15.822 8.113 ( -2.379 -4.121 2.375) 5.319 8.165 ( -3.187 -5.520 -1.565) 6.563 ======================= Grid point 258 (37/54) ======================= q-point: ( 0.13 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.850 ( 4.068 4.436 3.164) 6.799 0.897 ( 3.770 6.157 4.536) 8.527 0.942 ( -0.500 -1.334 -3.958) 4.207 1.011 ( 1.279 3.746 -4.214) 5.781 1.299 ( 2.277 0.926 2.436) 3.461 1.349 ( -2.470 -4.956 -1.614) 5.768 1.787 ( 6.374 10.354 0.144) 12.160 1.790 ( 7.123 9.875 -0.076) 12.176 2.114 ( -1.504 1.175 2.392) 3.060 2.122 ( -1.551 2.872 2.432) 4.070 2.520 ( 21.237 19.757 3.986) 29.278 2.637 ( 18.699 17.159 -5.010) 25.869 4.847 ( -8.631 11.144 2.625) 14.338 4.963 ( -3.669 5.207 -7.170) 9.591 5.438 ( 13.261 -0.716 4.962) 14.177 5.499 ( 18.047 -1.535 2.262) 18.253 5.820 ( 4.829 6.418 6.455) 10.304 5.920 ( -0.954 -0.489 -2.145) 2.398 6.125 ( 2.656 18.019 -0.683) 18.226 6.282 ( 12.346 11.476 -10.941) 20.096 7.683 ( -12.791 -7.054 6.484) 15.982 7.839 ( -2.569 -9.073 -6.338) 11.362 8.058 ( -0.496 -2.469 1.563) 2.965 8.093 ( -1.722 -2.927 -1.127) 3.579 ======================= Grid point 259 (38/54) ======================= q-point: ( 0.20 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.912 ( 0.141 -1.280 2.521) 2.831 0.942 ( 2.206 2.126 -0.077) 3.065 1.006 ( 1.877 4.483 2.034) 5.269 1.060 ( 0.026 3.834 -2.165) 4.403 1.222 ( -3.818 -6.402 3.499) 8.235 1.323 ( 1.263 -1.357 -4.399) 4.774 1.997 ( 6.482 6.674 0.358) 9.310 2.006 ( 6.483 8.969 -0.263) 11.070 2.171 ( 2.342 7.531 0.798) 7.927 2.194 ( 0.784 10.289 -0.923) 10.360 3.081 ( 22.077 15.206 2.381) 26.913 3.145 ( 20.649 13.366 -2.700) 24.745 4.762 ( -12.660 10.566 5.580) 17.409 4.943 ( -4.824 6.102 -9.384) 12.189 5.583 ( 2.412 -5.074 4.416) 7.146 5.682 ( 0.991 -9.848 -0.657) 9.919 5.896 ( 5.782 -2.294 3.207) 6.998 5.919 ( 10.881 -0.889 5.779) 12.352 6.351 ( 0.593 18.830 1.584) 18.906 6.594 ( 8.333 12.558 -16.931) 22.667 7.513 ( -2.652 0.536 11.024) 11.351 7.720 ( -1.656 -6.079 -6.372) 8.962 8.031 ( -0.562 -2.612 0.666) 2.753 8.040 ( 0.058 -3.349 0.128) 3.352 ======================= Grid point 260 (39/54) ======================= q-point: ( 0.27 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.929 ( 1.928 1.090 5.226) 5.676 0.957 ( -1.105 -0.369 0.868) 1.453 1.064 ( 0.920 2.670 1.313) 3.115 1.096 ( -0.056 3.750 -0.935) 3.865 1.136 ( -0.381 -1.799 5.679) 5.969 1.342 ( 2.133 -0.181 -8.730) 8.989 2.131 ( 4.107 1.151 1.534) 4.532 2.170 ( 5.029 2.830 -1.547) 5.975 2.308 ( 0.994 11.911 0.257) 11.955 2.318 ( -0.401 12.778 -0.439) 12.791 3.578 ( 19.722 10.325 0.954) 22.282 3.603 ( 18.741 8.742 -0.904) 20.699 4.644 ( -11.377 8.639 10.583) 17.778 4.926 ( -4.276 6.530 -12.174) 14.461 5.360 ( -14.360 -8.553 7.811) 18.449 5.551 ( -2.612 -7.385 -6.354) 10.086 6.033 ( 10.402 -2.578 8.015) 13.382 6.160 ( 11.937 2.399 2.697) 12.471 6.510 ( -2.510 16.205 2.069) 16.529 6.786 ( 1.203 7.930 -19.967) 21.518 7.544 ( 3.882 -0.017 9.313) 10.090 7.678 ( 1.264 -1.486 -2.975) 3.557 7.989 ( -0.593 -2.514 2.239) 3.419 8.027 ( 1.595 -3.088 -0.905) 3.592 ======================= Grid point 261 (40/54) ======================= q-point: ( 0.33 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.929 ( -0.860 -0.129 2.355) 2.510 0.982 ( 1.874 0.737 -1.236) 2.363 1.110 ( 0.124 3.772 1.232) 3.970 1.130 ( -0.259 3.123 -0.542) 3.181 1.171 ( 4.723 2.240 9.337) 10.701 1.389 ( 2.690 0.670 -7.588) 8.078 2.164 ( 0.509 -1.980 2.939) 3.580 2.232 ( 1.941 -2.590 -2.549) 4.120 2.417 ( -2.629 12.606 0.238) 12.879 2.426 ( -2.049 11.975 -0.455) 12.157 3.956 ( 14.912 4.109 0.798) 15.488 3.973 ( 15.330 5.394 -0.390) 16.256 4.558 ( -7.404 7.512 15.534) 18.777 4.924 ( -4.154 6.701 -10.000) 12.734 5.046 ( -13.772 -2.173 11.230) 17.902 5.451 ( -1.140 -5.631 -17.055) 17.997 6.169 ( 7.631 -2.025 9.989) 12.733 6.354 ( 6.249 4.622 1.771) 7.972 6.594 ( -4.328 12.707 -0.142) 13.425 6.810 ( -3.337 2.547 -16.490) 17.016 7.611 ( 3.828 -0.875 5.903) 7.090 7.705 ( 1.972 -0.342 -2.439) 3.155 7.966 ( 0.545 -1.818 4.336) 4.734 8.031 ( 1.740 -1.851 -1.291) 2.850 ======================= Grid point 262 (41/54) ======================= q-point: ( 0.40 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.925 ( 0.018 0.430 3.096) 3.126 1.007 ( 1.368 -0.838 -3.717) 4.048 1.137 ( -0.801 3.249 0.566) 3.394 1.148 ( -0.627 2.065 -0.275) 2.175 1.322 ( 6.860 4.189 6.231) 10.170 1.438 ( 1.985 0.194 -3.445) 3.981 2.130 ( -1.457 -1.763 3.054) 3.816 2.197 ( -1.006 -4.396 -2.396) 5.106 2.475 ( -4.386 11.224 0.170) 12.052 2.482 ( -4.551 10.890 -0.388) 11.809 4.186 ( 9.434 -0.888 1.832) 9.651 4.222 ( 9.431 0.613 -1.087) 9.513 4.543 ( -4.329 9.524 17.425) 20.324 4.804 ( -8.979 -0.847 22.286) 24.042 4.955 ( -5.008 9.598 -16.040) 19.352 5.398 ( 0.904 -4.788 -24.972) 25.443 6.256 ( 4.385 -1.521 11.110) 12.041 6.476 ( 0.683 7.902 4.010) 8.888 6.619 ( -4.360 7.533 -4.891) 9.984 6.762 ( -2.921 0.298 -11.374) 11.747 7.648 ( 1.538 -1.038 2.408) 3.040 7.728 ( 1.279 -0.880 -3.296) 3.643 7.967 ( 1.031 -1.074 6.249) 6.424 8.045 ( 1.163 -0.745 -0.958) 1.681 ======================= Grid point 263 (42/54) ======================= q-point: ( 0.47 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.929 ( -0.130 0.226 3.442) 3.452 1.015 ( 0.870 -1.507 -3.568) 3.970 1.148 ( -1.995 3.455 0.228) 3.996 1.151 ( -0.375 0.650 0.113) 0.759 1.426 ( 0.341 -0.591 1.764) 1.892 1.449 ( 0.413 -0.715 -0.227) 0.857 2.104 ( -0.066 0.114 2.557) 2.561 2.155 ( 1.120 -1.940 -1.543) 2.720 2.493 ( -5.858 10.146 0.033) 11.716 2.496 ( -5.808 10.059 -0.181) 11.616 4.264 ( 3.015 -5.222 2.199) 6.419 4.308 ( 2.292 -3.970 -1.331) 4.773 4.568 ( -7.274 12.599 15.353) 21.151 4.708 ( 0.190 -0.330 26.298) 26.301 4.962 ( -5.655 9.795 -15.736) 19.379 5.382 ( 2.369 -4.104 -27.796) 28.197 6.288 ( 1.304 -2.259 11.674) 11.962 6.528 ( -6.290 10.894 9.482) 15.753 6.611 ( -1.274 2.207 -11.424) 11.705 6.735 ( -0.568 0.985 -9.385) 9.453 7.654 ( 0.651 -1.128 1.356) 1.880 7.734 ( 0.655 -1.134 -3.925) 4.138 7.972 ( 0.575 -0.996 6.923) 7.018 8.052 ( 0.358 -0.620 -0.708) 1.007 ======================= Grid point 274 (43/54) ======================= q-point: ( 0.13 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.903 ( -0.572 -0.991 2.410) 2.667 0.932 ( 1.955 3.387 0.026) 3.911 1.017 ( 2.672 4.628 2.792) 6.029 1.085 ( 1.954 3.384 -2.566) 4.674 1.217 ( -4.152 -7.191 2.505) 8.673 1.299 ( -0.210 -0.365 -3.893) 3.915 1.972 ( 4.353 7.540 0.850) 8.747 1.994 ( 4.970 8.609 -0.838) 9.976 2.209 ( 5.073 8.787 2.294) 10.402 2.262 ( 5.759 9.975 -1.931) 11.679 2.945 ( 12.557 21.749 2.768) 25.266 3.012 ( 11.169 19.346 -2.556) 22.484 5.018 ( -0.411 -0.712 -2.736) 2.857 5.071 ( 3.735 6.469 -7.098) 10.304 5.440 ( 1.432 2.480 9.366) 9.794 5.554 ( 4.611 7.987 1.962) 9.429 5.802 ( -3.406 -5.899 5.608) 8.823 5.868 ( -1.811 -3.136 -1.028) 3.765 6.492 ( 8.795 15.234 -3.192) 17.878 6.548 ( 8.637 14.960 -7.060) 18.661 7.591 ( -1.324 -2.294 3.158) 4.122 7.661 ( -4.346 -7.527 -2.541) 9.056 8.005 ( -1.501 -2.599 1.938) 3.573 8.031 ( -1.807 -3.129 -0.322) 3.627 ======================= Grid point 275 (44/54) ======================= q-point: ( 0.20 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.914 ( 1.454 1.811 6.871) 7.253 0.973 ( -0.100 0.916 -0.503) 1.050 1.072 ( -0.671 -1.023 4.066) 4.246 1.113 ( 0.143 0.452 3.336) 3.369 1.149 ( 0.846 3.489 -1.314) 3.823 1.307 ( 1.498 0.194 -8.889) 9.016 2.089 ( 2.953 2.844 1.131) 4.252 2.122 ( 3.391 3.051 -1.510) 4.805 2.422 ( 3.558 15.355 1.865) 15.872 2.467 ( 2.201 15.043 -1.730) 15.302 3.374 ( 15.529 13.725 1.226) 20.761 3.396 ( 14.584 11.518 -0.495) 18.590 4.903 ( -9.352 -1.713 5.232) 10.852 5.112 ( -6.120 9.996 -10.754) 15.906 5.488 ( 1.728 -5.511 9.246) 10.902 5.561 ( -12.705 3.756 3.282) 13.649 5.859 ( 9.090 0.087 -3.510) 9.745 5.890 ( 14.914 -2.504 -0.549) 15.133 6.713 ( -0.247 16.328 -0.025) 16.330 6.836 ( 6.478 11.101 -9.688) 16.095 7.511 ( -2.680 -3.037 6.640) 7.778 7.623 ( 3.023 -1.393 -3.152) 4.584 7.963 ( 0.462 -3.762 1.609) 4.117 7.973 ( -0.959 -3.009 0.708) 3.236 ======================= Grid point 276 (45/54) ======================= q-point: ( 0.27 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.946 ( -0.958 -0.134 1.194) 1.537 0.969 ( 0.480 1.217 0.457) 1.386 1.096 ( 4.587 2.511 10.418) 11.657 1.146 ( 0.157 2.613 2.043) 3.321 1.192 ( -0.483 4.338 -1.263) 4.544 1.352 ( 2.495 1.266 -8.761) 9.197 2.134 ( 1.149 -0.615 1.540) 2.017 2.175 ( 1.784 -1.611 -1.747) 2.971 2.607 ( -1.088 16.490 0.780) 16.544 2.625 ( -2.376 16.401 -0.595) 16.582 3.732 ( 14.583 4.704 1.197) 15.370 3.752 ( 15.428 6.360 -0.355) 16.691 4.725 ( -6.958 -2.538 16.024) 17.653 5.115 ( -5.136 10.803 -12.629) 17.395 5.337 ( -16.223 6.898 6.449) 18.771 5.430 ( -3.569 -3.204 -2.924) 5.617 6.016 ( 8.865 1.122 -0.016) 8.935 6.173 ( 15.677 -1.523 -8.194) 17.754 6.826 ( -4.517 14.468 1.969) 15.284 6.950 ( -4.459 8.639 -7.445) 12.245 7.506 ( 5.320 -3.057 8.974) 10.871 7.664 ( 2.084 -0.377 -4.835) 5.279 7.932 ( 0.447 -2.775 2.680) 3.883 7.964 ( 2.718 -2.492 0.178) 3.692 ======================= Grid point 277 (46/54) ======================= q-point: ( 0.33 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.935 ( -0.431 0.553 1.548) 1.699 0.979 ( 0.691 -0.869 -2.199) 2.464 1.157 ( -0.150 0.482 3.104) 3.145 1.211 ( -0.923 4.515 -0.887) 4.693 1.258 ( 5.607 5.673 6.275) 10.149 1.406 ( 2.093 1.014 -5.037) 5.548 2.129 ( -0.258 -1.232 1.432) 1.907 2.160 ( 0.143 -3.870 -1.099) 4.025 2.714 ( -5.520 16.295 0.124) 17.205 2.718 ( -5.589 16.242 -0.095) 17.177 3.982 ( 11.850 -1.239 2.410) 12.156 4.031 ( 12.289 0.296 -1.499) 12.383 4.645 ( -1.370 1.601 21.899) 22.000 5.116 ( -18.395 12.190 10.029) 24.240 5.152 ( -3.873 9.198 -14.298) 17.437 5.350 ( -1.887 -2.953 -11.786) 12.296 6.162 ( 6.094 1.299 1.995) 6.543 6.385 ( 11.058 -2.724 -12.357) 16.804 6.841 ( -8.245 10.732 3.419) 13.959 6.920 ( -7.696 9.596 -1.719) 12.420 7.581 ( 4.726 -1.802 3.721) 6.279 7.685 ( 1.804 -1.723 -4.426) 5.081 7.931 ( 1.560 -1.243 5.377) 5.734 8.000 ( 2.911 -0.909 -0.675) 3.124 ======================= Grid point 278 (47/54) ======================= q-point: ( 0.40 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.936 ( -0.195 0.497 1.635) 1.720 0.981 ( 0.968 -1.602 -2.128) 2.834 1.154 ( 0.133 -0.665 3.651) 3.714 1.224 ( -1.944 3.900 -2.197) 4.880 1.392 ( 3.637 2.312 2.220) 4.848 1.436 ( 0.642 -0.347 -1.366) 1.549 2.113 ( -0.637 0.140 1.077) 1.260 2.125 ( -0.068 -2.230 0.135) 2.235 2.755 ( -8.057 15.809 0.103) 17.744 2.759 ( -8.139 15.867 -0.175) 17.834 4.114 ( 7.214 -5.829 3.045) 9.762 4.179 ( 7.103 -4.814 -2.096) 8.833 4.694 ( 1.880 4.772 21.509) 22.112 4.922 ( -15.489 11.559 11.960) 22.728 5.189 ( -5.542 12.108 -8.639) 15.873 5.306 ( 0.406 -3.093 -19.948) 20.190 6.247 ( 2.349 0.505 3.367) 4.136 6.501 ( 7.620 -7.088 -13.553) 17.087 6.797 ( -8.979 12.562 6.048) 16.584 6.889 ( -7.263 10.305 -2.580) 12.869 7.619 ( 1.849 -1.769 0.981) 2.740 7.697 ( 1.714 -2.211 -4.415) 5.227 7.949 ( 1.061 -0.452 7.202) 7.294 8.033 ( 1.309 -0.435 -0.567) 1.492 ======================= Grid point 290 (48/54) ======================= q-point: ( 0.20 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.952 ( 0.159 0.275 1.529) 1.561 0.979 ( -0.245 -0.425 -1.112) 1.216 1.061 ( 2.420 4.191 12.559) 13.459 1.160 ( 1.289 2.233 2.818) 3.819 1.220 ( 1.713 2.966 -2.150) 4.044 1.335 ( 1.431 2.478 -9.274) 9.705 2.119 ( 0.257 0.445 0.944) 1.074 2.149 ( -0.121 -0.210 -1.405) 1.425 2.719 ( 7.610 13.181 1.342) 15.279 2.749 ( 7.058 12.224 -1.089) 14.157 3.610 ( 5.683 9.843 1.450) 11.458 3.631 ( 6.791 11.763 -0.339) 13.587 4.766 ( -4.654 -8.062 15.856) 18.387 5.215 ( 0.389 0.674 -20.645) 20.659 5.445 ( -0.244 -0.422 17.715) 17.722 5.711 ( 1.433 2.482 7.672) 8.190 5.855 ( 2.595 4.495 -13.502) 14.465 5.937 ( 4.399 7.619 -11.102) 14.165 6.963 ( 4.505 7.803 0.551) 9.027 7.051 ( 5.589 9.680 -4.625) 12.096 7.422 ( -1.814 -3.141 11.542) 12.099 7.640 ( 0.641 1.110 -7.673) 7.779 7.905 ( -0.929 -1.609 2.784) 3.347 7.926 ( -0.656 -1.136 1.195) 1.775 ======================= Grid point 291 (49/54) ======================= q-point: ( 0.27 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.946 ( -0.338 0.116 0.647) 0.739 0.967 ( -0.225 -0.832 -1.027) 1.341 1.153 ( 0.833 0.568 6.248) 6.329 1.230 ( 3.811 6.460 6.458) 9.898 1.265 ( 0.343 2.773 -3.224) 4.266 1.390 ( 1.625 2.201 -5.765) 6.382 2.120 ( 0.151 -0.364 0.136) 0.417 2.127 ( -0.298 -2.471 -0.382) 2.518 2.935 ( 0.552 15.409 0.507) 15.427 2.946 ( -0.144 14.888 -0.344) 14.892 3.790 ( 8.300 1.379 2.625) 8.814 3.843 ( 8.898 2.709 -1.683) 9.452 4.665 ( -0.926 -0.812 21.422) 21.457 5.210 ( -0.361 -0.615 -22.326) 22.337 5.422 ( -3.353 2.376 13.909) 14.504 5.596 ( -18.432 14.149 4.227) 23.618 6.047 ( 12.881 -3.350 -11.321) 17.473 6.123 ( 9.655 1.719 -13.423) 16.624 7.047 ( -3.049 7.352 6.634) 10.362 7.155 ( -5.191 10.994 -1.880) 12.302 7.451 ( 7.028 -1.881 7.019) 10.110 7.642 ( 0.489 -2.030 -8.665) 8.913 7.899 ( 1.062 -0.047 4.322) 4.451 7.942 ( 2.631 0.469 0.988) 2.849 ======================= Grid point 292 (50/54) ======================= q-point: ( 0.33 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.944 ( -0.243 0.209 0.248) 0.405 0.956 ( 0.159 -1.146 -0.805) 1.409 1.136 ( -0.309 -2.600 7.366) 7.818 1.283 ( -0.864 2.630 -4.607) 5.375 1.361 ( 3.389 3.833 2.878) 5.870 1.425 ( 0.644 0.701 -2.257) 2.449 2.109 ( 0.234 -0.582 0.972) 1.156 2.125 ( 0.017 0.889 -0.280) 0.932 3.049 ( -5.668 16.271 0.114) 17.231 3.054 ( -5.441 16.609 -0.245) 17.479 3.918 ( 8.092 -4.830 3.138) 9.932 3.987 ( 8.351 -4.517 -2.324) 9.774 4.713 ( 2.360 6.020 18.556) 19.650 5.172 ( -3.119 -2.470 -17.478) 17.925 5.385 ( -6.065 6.066 7.662) 11.502 5.435 ( -13.067 12.455 6.733) 19.267 6.201 ( 8.104 -2.552 -9.273) 12.576 6.260 ( 6.816 -4.145 -10.827) 13.449 7.048 ( -5.512 10.218 5.610) 12.894 7.141 ( -8.457 10.590 -2.209) 13.731 7.546 ( 4.731 -1.640 1.086) 5.124 7.633 ( 1.927 -3.464 -6.131) 7.301 7.932 ( 1.368 1.477 6.668) 6.965 8.000 ( 2.463 0.746 0.680) 2.662 ======================= Grid point 293 (51/54) ======================= q-point: ( 0.40 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.942 ( -0.066 0.114 0.072) 0.150 0.951 ( 0.662 -1.147 -0.659) 1.479 1.117 ( 1.558 -2.698 8.735) 9.274 1.288 ( -1.376 2.383 -5.351) 6.017 1.415 ( -0.090 0.156 0.675) 0.698 1.429 ( 0.060 -0.105 -0.407) 0.424 2.112 ( -0.349 0.605 1.364) 1.533 2.131 ( -0.789 1.366 -0.182) 1.588 3.080 ( -9.076 15.719 0.225) 18.153 3.089 ( -9.296 16.101 -0.510) 18.599 3.966 ( 4.911 -8.507 3.229) 10.340 4.038 ( 5.057 -8.760 -2.459) 10.410 4.794 ( -3.044 5.272 12.995) 14.351 5.097 ( -1.726 2.990 -6.714) 7.550 5.351 ( -5.653 9.792 2.514) 11.583 5.401 ( -4.531 7.848 6.173) 10.965 6.277 ( -1.330 2.303 -10.820) 11.142 6.279 ( 8.010 -13.873 -5.580) 16.963 7.071 ( -7.510 13.008 2.210) 15.182 7.098 ( -5.825 10.090 -3.139) 12.066 7.576 ( 1.441 -2.496 -0.289) 2.897 7.636 ( 2.250 -3.898 -4.895) 6.650 7.955 ( -0.568 0.985 7.555) 7.640 8.027 ( 0.146 -0.254 0.861) 0.909 ======================= Grid point 307 (52/54) ======================= q-point: ( 0.27 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.947 ( -0.059 -0.102 -0.184) 0.219 0.950 ( -0.420 -0.727 -0.372) 0.918 1.124 ( -1.709 -2.959 8.947) 9.577 1.304 ( 0.747 1.294 -5.887) 6.074 1.350 ( 2.634 4.563 3.457) 6.302 1.424 ( 0.575 0.995 -2.571) 2.817 2.104 ( 0.283 0.490 1.212) 1.338 2.129 ( 0.705 1.222 -0.846) 1.645 3.205 ( 6.016 10.420 0.021) 12.032 3.207 ( 6.442 11.157 -0.222) 12.885 3.800 ( -0.037 -0.064 3.038) 3.039 3.868 ( -0.013 -0.023 -2.344) 2.344 4.724 ( 3.643 6.310 17.426) 18.888 5.174 ( -1.407 -2.438 -19.537) 19.739 5.453 ( 0.559 0.968 17.571) 17.607 5.813 ( 0.787 1.363 15.600) 15.679 5.894 ( -2.133 -3.695 -16.376) 16.923 6.235 ( 3.491 6.047 -18.985) 20.228 7.159 ( 2.523 4.369 7.071) 8.687 7.321 ( 2.270 3.932 -4.854) 6.646 7.450 ( 1.353 2.344 3.945) 4.784 7.589 ( -1.583 -2.741 -8.658) 9.219 7.930 ( 1.520 2.633 5.808) 6.556 7.977 ( 1.442 2.497 2.092) 3.563 ======================= Grid point 308 (53/54) ======================= q-point: ( 0.33 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.939 ( -0.144 -0.471 -0.294) 0.574 0.942 ( -0.118 -0.314 -0.658) 0.738 1.066 ( -0.094 -3.637 13.097) 13.593 1.320 ( -0.257 1.126 -7.901) 7.985 1.410 ( 1.080 1.383 1.158) 2.102 1.431 ( -0.196 -0.059 -0.551) 0.588 2.128 ( 0.293 1.731 1.563) 2.351 2.156 ( 0.094 1.695 -0.650) 1.818 3.352 ( -3.343 13.473 0.304) 13.885 3.368 ( -3.187 14.175 -1.074) 14.569 3.796 ( 3.879 -7.121 2.815) 8.583 3.857 ( 3.978 -8.180 -2.028) 9.319 4.878 ( 2.681 9.338 6.271) 11.563 5.114 ( -1.987 -2.977 -13.414) 13.883 5.489 ( 0.670 5.308 15.924) 16.799 5.633 ( -14.445 6.187 9.764) 18.500 5.966 ( 11.032 -13.846 4.971) 18.388 6.324 ( 0.413 3.979 -22.460) 22.814 7.241 ( -0.391 8.680 -1.017) 8.749 7.318 ( -7.039 6.703 -4.736) 10.812 7.512 ( 3.879 -1.952 -0.104) 4.344 7.549 ( 1.295 -4.538 -4.526) 6.538 7.979 ( 0.343 2.681 6.759) 7.279 8.020 ( 0.972 0.860 3.151) 3.408 ======================= Grid point 323 (54/54) ======================= q-point: ( 0.33 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 1.25e-05 1.25e-05 2.00e-05 2.00e-05 Number of triplets: 73 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.937 ( -0.060 -0.104 -0.738) 0.747 0.937 ( 0.060 0.104 -0.738) 0.747 1.023 ( -0.000 -0.000 16.270) 16.270 1.331 ( -0.000 -0.000 -9.585) 9.585 1.427 ( -0.203 -0.351 0.523) 0.662 1.427 ( 0.203 0.351 0.523) 0.662 2.150 ( -0.000 -0.000 1.546) 1.546 2.175 ( -0.000 -0.000 -0.535) 0.535 3.558 ( -1.332 -2.308 -1.212) 2.927 3.558 ( 1.332 2.307 -1.212) 2.926 3.690 ( -1.933 -3.348 1.079) 4.014 3.690 ( 1.933 3.349 1.079) 4.014 5.036 ( -2.773 -4.804 -6.329) 8.416 5.036 ( 2.773 4.804 -6.329) 8.415 5.632 ( -3.831 -6.635 10.612) 13.089 5.632 ( 3.831 6.635 10.612) 13.088 5.781 ( 0.000 0.000 21.437) 21.437 6.368 ( -0.000 -0.000 -26.607) 26.607 7.378 ( -1.517 -2.627 -9.725) 10.187 7.378 ( 1.516 2.627 -9.725) 10.187 7.460 ( 0.000 0.000 4.666) 4.666 7.497 ( 0.000 0.000 0.604) 0.604 8.020 ( -0.521 -0.903 5.506) 5.604 8.020 ( 0.521 0.903 5.506) 5.604 =================== 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/16200 10.0 171.946 171.946 148.719 0.000 0.000 -0.000 3/16200 20.0 35.503 35.503 36.848 0.000 0.000 -0.000 3/16200 30.0 17.970 17.970 18.218 0.000 0.000 -0.000 3/16200 40.0 12.124 12.124 11.723 0.000 0.000 -0.000 3/16200 50.0 9.052 9.052 8.514 0.000 0.000 -0.000 3/16200 60.0 7.147 7.147 6.639 -0.000 0.000 -0.000 3/16200 70.0 5.876 5.876 5.430 -0.000 0.000 -0.000 3/16200 80.0 4.982 4.982 4.596 -0.000 0.000 -0.000 3/16200 90.0 4.324 4.324 3.988 -0.000 0.000 -0.000 3/16200 100.0 3.823 3.823 3.527 -0.000 0.000 -0.000 3/16200 110.0 3.428 3.428 3.164 -0.000 0.000 -0.000 3/16200 120.0 3.109 3.109 2.872 -0.000 0.000 -0.000 3/16200 130.0 2.846 2.846 2.631 -0.000 0.000 -0.000 3/16200 140.0 2.625 2.625 2.428 -0.000 0.000 -0.000 3/16200 150.0 2.437 2.437 2.255 -0.000 0.000 -0.000 3/16200 160.0 2.274 2.274 2.106 -0.000 0.000 -0.000 3/16200 170.0 2.133 2.133 1.976 -0.000 0.000 -0.000 3/16200 180.0 2.008 2.008 1.861 -0.000 0.000 -0.000 3/16200 190.0 1.898 1.898 1.760 -0.000 0.000 -0.000 3/16200 200.0 1.799 1.799 1.668 -0.000 0.000 -0.000 3/16200 210.0 1.710 1.710 1.587 -0.000 0.000 -0.000 3/16200 220.0 1.630 1.630 1.512 -0.000 0.000 -0.000 3/16200 230.0 1.557 1.557 1.445 -0.000 0.000 -0.000 3/16200 240.0 1.490 1.490 1.383 -0.000 0.000 -0.000 3/16200 250.0 1.429 1.429 1.327 -0.000 0.000 -0.000 3/16200 260.0 1.373 1.373 1.275 -0.000 0.000 -0.000 3/16200 270.0 1.321 1.321 1.227 -0.000 0.000 -0.000 3/16200 280.0 1.273 1.273 1.182 -0.000 0.000 -0.000 3/16200 290.0 1.228 1.228 1.141 -0.000 0.000 -0.000 3/16200 300.0 1.186 1.186 1.102 -0.000 0.000 -0.000 3/16200 310.0 1.147 1.147 1.066 -0.000 0.000 -0.000 3/16200 320.0 1.111 1.111 1.033 -0.000 0.000 -0.000 3/16200 330.0 1.077 1.077 1.001 -0.000 0.000 -0.000 3/16200 340.0 1.045 1.045 0.971 -0.000 0.000 -0.000 3/16200 350.0 1.015 1.015 0.943 -0.000 0.000 -0.000 3/16200 360.0 0.986 0.986 0.917 -0.000 0.000 -0.000 3/16200 370.0 0.959 0.959 0.892 -0.000 0.000 -0.000 3/16200 380.0 0.934 0.934 0.868 -0.000 0.000 -0.000 3/16200 390.0 0.909 0.909 0.846 -0.000 0.000 -0.000 3/16200 400.0 0.887 0.887 0.825 -0.000 0.000 -0.000 3/16200 410.0 0.865 0.865 0.804 -0.000 0.000 -0.000 3/16200 420.0 0.844 0.844 0.785 -0.000 0.000 -0.000 3/16200 430.0 0.824 0.824 0.767 -0.000 0.000 -0.000 3/16200 440.0 0.805 0.805 0.749 -0.000 0.000 -0.000 3/16200 450.0 0.787 0.787 0.732 -0.000 0.000 -0.000 3/16200 460.0 0.770 0.770 0.716 -0.000 0.000 -0.000 3/16200 470.0 0.754 0.754 0.701 -0.000 0.000 -0.000 3/16200 480.0 0.738 0.738 0.686 -0.000 0.000 -0.000 3/16200 490.0 0.723 0.723 0.672 -0.000 0.000 -0.000 3/16200 500.0 0.708 0.708 0.659 -0.000 0.000 -0.000 3/16200 510.0 0.694 0.694 0.646 -0.000 0.000 -0.000 3/16200 520.0 0.681 0.681 0.633 -0.000 0.000 -0.000 3/16200 530.0 0.668 0.668 0.621 -0.000 0.000 -0.000 3/16200 540.0 0.655 0.655 0.610 -0.000 0.000 -0.000 3/16200 550.0 0.643 0.643 0.599 -0.000 0.000 -0.000 3/16200 560.0 0.632 0.632 0.588 -0.000 0.000 -0.000 3/16200 570.0 0.621 0.621 0.578 -0.000 0.000 -0.000 3/16200 580.0 0.610 0.610 0.568 -0.000 0.000 -0.000 3/16200 590.0 0.600 0.600 0.558 -0.000 0.000 -0.000 3/16200 600.0 0.590 0.590 0.549 -0.000 0.000 -0.000 3/16200 610.0 0.580 0.580 0.540 -0.000 0.000 -0.000 3/16200 620.0 0.570 0.570 0.531 -0.000 0.000 -0.000 3/16200 630.0 0.561 0.561 0.523 -0.000 0.000 -0.000 3/16200 640.0 0.553 0.553 0.514 -0.000 0.000 -0.000 3/16200 650.0 0.544 0.544 0.506 -0.000 0.000 -0.000 3/16200 660.0 0.536 0.536 0.499 -0.000 0.000 -0.000 3/16200 670.0 0.528 0.528 0.491 -0.000 0.000 -0.000 3/16200 680.0 0.520 0.520 0.484 -0.000 0.000 -0.000 3/16200 690.0 0.512 0.512 0.477 -0.000 0.000 -0.000 3/16200 700.0 0.505 0.505 0.470 -0.000 0.000 -0.000 3/16200 710.0 0.498 0.498 0.464 -0.000 0.000 -0.000 3/16200 720.0 0.491 0.491 0.457 -0.000 0.000 -0.000 3/16200 730.0 0.484 0.484 0.451 -0.000 0.000 -0.000 3/16200 740.0 0.478 0.478 0.445 -0.000 0.000 -0.000 3/16200 750.0 0.471 0.471 0.439 -0.000 0.000 -0.000 3/16200 760.0 0.465 0.465 0.433 -0.000 0.000 -0.000 3/16200 770.0 0.459 0.459 0.427 -0.000 0.000 -0.000 3/16200 780.0 0.453 0.453 0.422 -0.000 0.000 -0.000 3/16200 790.0 0.447 0.447 0.417 -0.000 0.000 -0.000 3/16200 800.0 0.442 0.442 0.411 -0.000 0.000 -0.000 3/16200 810.0 0.436 0.436 0.406 -0.000 0.000 -0.000 3/16200 820.0 0.431 0.431 0.401 -0.000 0.000 -0.000 3/16200 830.0 0.426 0.426 0.396 -0.000 0.000 -0.000 3/16200 840.0 0.421 0.421 0.392 -0.000 0.000 -0.000 3/16200 850.0 0.416 0.416 0.387 -0.000 0.000 -0.000 3/16200 860.0 0.411 0.411 0.383 -0.000 0.000 -0.000 3/16200 870.0 0.406 0.406 0.378 -0.000 0.000 -0.000 3/16200 880.0 0.402 0.402 0.374 -0.000 0.000 -0.000 3/16200 890.0 0.397 0.397 0.370 -0.000 0.000 -0.000 3/16200 900.0 0.393 0.393 0.366 -0.000 0.000 -0.000 3/16200 910.0 0.388 0.388 0.362 -0.000 0.000 -0.000 3/16200 920.0 0.384 0.384 0.358 -0.000 0.000 -0.000 3/16200 930.0 0.380 0.380 0.354 -0.000 0.000 -0.000 3/16200 940.0 0.376 0.376 0.350 -0.000 0.000 -0.000 3/16200 950.0 0.372 0.372 0.346 -0.000 0.000 -0.000 3/16200 960.0 0.368 0.368 0.343 -0.000 0.000 -0.000 3/16200 970.0 0.364 0.364 0.339 -0.000 0.000 -0.000 3/16200 980.0 0.361 0.361 0.336 -0.000 0.000 -0.000 3/16200 990.0 0.357 0.357 0.332 -0.000 0.000 -0.000 3/16200 1000.0 0.353 0.353 0.329 -0.000 0.000 -0.000 3/16200 Thermal conductivity related properties were written into "kappa-m15153.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 04:51:15]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|