----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 18:53:03]------------------------- 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: P-6m2 (187) Number of symmetry operations in supercell: 192 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.747820461219762 0.000000000000000 0.000000000000000 b -1.873910230609881 3.245707728239425 0.000000000000000 c 0.000000000000000 0.000000000000000 16.122787270000000 Atomic positions (fractional): *1 Ga 0.00000000000000 0.00000000000000 0.92470939142131 69.723 2 Ga 0.00000000000000 0.00000000000000 0.07529060857869 69.723 *3 Ga 0.66666666666667 0.33333333333333 0.57525982534184 69.723 4 Ga 0.66666666666667 0.33333333333333 0.42474017465816 69.723 *5 Se 0.66666666666667 0.33333333333333 0.14798852090080 78.960 6 Se 0.66666666666667 0.33333333333333 0.85201147909920 78.960 *7 Se 0.33333333333333 0.66666666666667 0.35203925517484 78.960 8 Se 0.33333333333333 0.66666666666667 0.64796074482516 78.960 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.747820461219762 0.000000000000000 0.000000000000000 b -1.873910230609881 3.245707728239425 0.000000000000000 c 0.000000000000000 0.000000000000000 16.122787270000000 Atomic positions (fractional): *1 Ga 0.00000000000000 0.00000000000000 0.92470939142131 69.723 > 1 2 Ga 0.00000000000000 0.00000000000000 0.07529060857869 69.723 > 2 *3 Ga 0.66666666666667 0.33333333333333 0.57525982534184 69.723 > 3 4 Ga 0.66666666666667 0.33333333333333 0.42474017465816 69.723 > 4 *5 Se 0.66666666666667 0.33333333333333 0.14798852090080 78.960 > 5 6 Se 0.66666666666667 0.33333333333333 0.85201147909920 78.960 > 6 *7 Se 0.33333333333333 0.66666666666667 0.35203925517484 78.960 > 7 8 Se 0.33333333333333 0.66666666666667 0.64796074482516 78.960 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 14.991281844879047 0.000000000000000 0.000000000000000 b -7.495640922439524 12.982830912957700 0.000000000000000 c 0.000000000000000 0.000000000000000 16.122787270000000 Atomic positions (fractional): *1 Ga 0.00000000000000 0.00000000000000 0.92470939142131 69.723 > 1 2 Ga 0.25000000000000 0.00000000000000 0.92470939142131 69.723 > 1 3 Ga 0.50000000000000 0.00000000000000 0.92470939142131 69.723 > 1 4 Ga 0.75000000000000 0.00000000000000 0.92470939142131 69.723 > 1 5 Ga 0.00000000000000 0.25000000000000 0.92470939142131 69.723 > 1 6 Ga 0.25000000000000 0.25000000000000 0.92470939142131 69.723 > 1 7 Ga 0.50000000000000 0.25000000000000 0.92470939142131 69.723 > 1 8 Ga 0.75000000000000 0.25000000000000 0.92470939142131 69.723 > 1 9 Ga 0.00000000000000 0.50000000000000 0.92470939142131 69.723 > 1 10 Ga 0.25000000000000 0.50000000000000 0.92470939142131 69.723 > 1 11 Ga 0.50000000000000 0.50000000000000 0.92470939142131 69.723 > 1 12 Ga 0.75000000000000 0.50000000000000 0.92470939142131 69.723 > 1 13 Ga 0.00000000000000 0.75000000000000 0.92470939142131 69.723 > 1 14 Ga 0.25000000000000 0.75000000000000 0.92470939142131 69.723 > 1 15 Ga 0.50000000000000 0.75000000000000 0.92470939142131 69.723 > 1 16 Ga 0.75000000000000 0.75000000000000 0.92470939142131 69.723 > 1 17 Ga 0.00000000000000 0.00000000000000 0.07529060857869 69.723 > 2 18 Ga 0.25000000000000 0.00000000000000 0.07529060857869 69.723 > 2 19 Ga 0.50000000000000 0.00000000000000 0.07529060857869 69.723 > 2 20 Ga 0.75000000000000 0.00000000000000 0.07529060857869 69.723 > 2 21 Ga 0.00000000000000 0.25000000000000 0.07529060857869 69.723 > 2 22 Ga 0.25000000000000 0.25000000000000 0.07529060857869 69.723 > 2 23 Ga 0.50000000000000 0.25000000000000 0.07529060857869 69.723 > 2 24 Ga 0.75000000000000 0.25000000000000 0.07529060857869 69.723 > 2 25 Ga 0.00000000000000 0.50000000000000 0.07529060857869 69.723 > 2 26 Ga 0.25000000000000 0.50000000000000 0.07529060857869 69.723 > 2 27 Ga 0.50000000000000 0.50000000000000 0.07529060857869 69.723 > 2 28 Ga 0.75000000000000 0.50000000000000 0.07529060857869 69.723 > 2 29 Ga 0.00000000000000 0.75000000000000 0.07529060857869 69.723 > 2 30 Ga 0.25000000000000 0.75000000000000 0.07529060857869 69.723 > 2 31 Ga 0.50000000000000 0.75000000000000 0.07529060857869 69.723 > 2 32 Ga 0.75000000000000 0.75000000000000 0.07529060857869 69.723 > 2 *33 Ga 0.16666666666667 0.08333333333333 0.57525982534184 69.723 > 3 34 Ga 0.41666666666667 0.08333333333333 0.57525982534184 69.723 > 3 35 Ga 0.66666666666667 0.08333333333333 0.57525982534184 69.723 > 3 36 Ga 0.91666666666667 0.08333333333333 0.57525982534184 69.723 > 3 37 Ga 0.16666666666667 0.33333333333333 0.57525982534184 69.723 > 3 38 Ga 0.41666666666667 0.33333333333333 0.57525982534184 69.723 > 3 39 Ga 0.66666666666667 0.33333333333333 0.57525982534184 69.723 > 3 40 Ga 0.91666666666667 0.33333333333333 0.57525982534184 69.723 > 3 41 Ga 0.16666666666667 0.58333333333333 0.57525982534184 69.723 > 3 42 Ga 0.41666666666667 0.58333333333333 0.57525982534184 69.723 > 3 43 Ga 0.66666666666667 0.58333333333333 0.57525982534184 69.723 > 3 44 Ga 0.91666666666667 0.58333333333333 0.57525982534184 69.723 > 3 45 Ga 0.16666666666667 0.83333333333333 0.57525982534184 69.723 > 3 46 Ga 0.41666666666667 0.83333333333333 0.57525982534184 69.723 > 3 47 Ga 0.66666666666667 0.83333333333333 0.57525982534184 69.723 > 3 48 Ga 0.91666666666667 0.83333333333333 0.57525982534184 69.723 > 3 49 Ga 0.16666666666667 0.08333333333333 0.42474017465816 69.723 > 4 50 Ga 0.41666666666667 0.08333333333333 0.42474017465816 69.723 > 4 51 Ga 0.66666666666667 0.08333333333333 0.42474017465816 69.723 > 4 52 Ga 0.91666666666667 0.08333333333333 0.42474017465816 69.723 > 4 53 Ga 0.16666666666667 0.33333333333333 0.42474017465816 69.723 > 4 54 Ga 0.41666666666667 0.33333333333333 0.42474017465816 69.723 > 4 55 Ga 0.66666666666667 0.33333333333333 0.42474017465816 69.723 > 4 56 Ga 0.91666666666667 0.33333333333333 0.42474017465816 69.723 > 4 57 Ga 0.16666666666667 0.58333333333333 0.42474017465816 69.723 > 4 58 Ga 0.41666666666667 0.58333333333333 0.42474017465816 69.723 > 4 59 Ga 0.66666666666667 0.58333333333333 0.42474017465816 69.723 > 4 60 Ga 0.91666666666667 0.58333333333333 0.42474017465816 69.723 > 4 61 Ga 0.16666666666667 0.83333333333333 0.42474017465816 69.723 > 4 62 Ga 0.41666666666667 0.83333333333333 0.42474017465816 69.723 > 4 63 Ga 0.66666666666667 0.83333333333333 0.42474017465816 69.723 > 4 64 Ga 0.91666666666667 0.83333333333333 0.42474017465816 69.723 > 4 *65 Se 0.16666666666667 0.08333333333333 0.14798852090080 78.960 > 5 66 Se 0.41666666666667 0.08333333333333 0.14798852090080 78.960 > 5 67 Se 0.66666666666667 0.08333333333333 0.14798852090080 78.960 > 5 68 Se 0.91666666666667 0.08333333333333 0.14798852090080 78.960 > 5 69 Se 0.16666666666667 0.33333333333333 0.14798852090080 78.960 > 5 70 Se 0.41666666666667 0.33333333333333 0.14798852090080 78.960 > 5 71 Se 0.66666666666667 0.33333333333333 0.14798852090080 78.960 > 5 72 Se 0.91666666666667 0.33333333333333 0.14798852090080 78.960 > 5 73 Se 0.16666666666667 0.58333333333333 0.14798852090080 78.960 > 5 74 Se 0.41666666666667 0.58333333333333 0.14798852090080 78.960 > 5 75 Se 0.66666666666667 0.58333333333333 0.14798852090080 78.960 > 5 76 Se 0.91666666666667 0.58333333333333 0.14798852090080 78.960 > 5 77 Se 0.16666666666667 0.83333333333333 0.14798852090080 78.960 > 5 78 Se 0.41666666666667 0.83333333333333 0.14798852090080 78.960 > 5 79 Se 0.66666666666667 0.83333333333333 0.14798852090080 78.960 > 5 80 Se 0.91666666666667 0.83333333333333 0.14798852090080 78.960 > 5 81 Se 0.16666666666667 0.08333333333333 0.85201147909920 78.960 > 6 82 Se 0.41666666666667 0.08333333333333 0.85201147909920 78.960 > 6 83 Se 0.66666666666667 0.08333333333333 0.85201147909920 78.960 > 6 84 Se 0.91666666666667 0.08333333333333 0.85201147909920 78.960 > 6 85 Se 0.16666666666667 0.33333333333333 0.85201147909920 78.960 > 6 86 Se 0.41666666666667 0.33333333333333 0.85201147909920 78.960 > 6 87 Se 0.66666666666667 0.33333333333333 0.85201147909920 78.960 > 6 88 Se 0.91666666666667 0.33333333333333 0.85201147909920 78.960 > 6 89 Se 0.16666666666667 0.58333333333333 0.85201147909920 78.960 > 6 90 Se 0.41666666666667 0.58333333333333 0.85201147909920 78.960 > 6 91 Se 0.66666666666667 0.58333333333333 0.85201147909920 78.960 > 6 92 Se 0.91666666666667 0.58333333333333 0.85201147909920 78.960 > 6 93 Se 0.16666666666667 0.83333333333333 0.85201147909920 78.960 > 6 94 Se 0.41666666666667 0.83333333333333 0.85201147909920 78.960 > 6 95 Se 0.66666666666667 0.83333333333333 0.85201147909920 78.960 > 6 96 Se 0.91666666666667 0.83333333333333 0.85201147909920 78.960 > 6 *97 Se 0.08333333333333 0.16666666666667 0.35203925517484 78.960 > 7 98 Se 0.33333333333333 0.16666666666667 0.35203925517484 78.960 > 7 99 Se 0.58333333333333 0.16666666666667 0.35203925517484 78.960 > 7 100 Se 0.83333333333333 0.16666666666667 0.35203925517484 78.960 > 7 101 Se 0.08333333333333 0.41666666666667 0.35203925517484 78.960 > 7 102 Se 0.33333333333333 0.41666666666667 0.35203925517484 78.960 > 7 103 Se 0.58333333333333 0.41666666666667 0.35203925517484 78.960 > 7 104 Se 0.83333333333333 0.41666666666667 0.35203925517484 78.960 > 7 105 Se 0.08333333333333 0.66666666666667 0.35203925517484 78.960 > 7 106 Se 0.33333333333333 0.66666666666667 0.35203925517484 78.960 > 7 107 Se 0.58333333333333 0.66666666666667 0.35203925517484 78.960 > 7 108 Se 0.83333333333333 0.66666666666667 0.35203925517484 78.960 > 7 109 Se 0.08333333333333 0.91666666666667 0.35203925517484 78.960 > 7 110 Se 0.33333333333333 0.91666666666667 0.35203925517484 78.960 > 7 111 Se 0.58333333333333 0.91666666666667 0.35203925517484 78.960 > 7 112 Se 0.83333333333333 0.91666666666667 0.35203925517484 78.960 > 7 113 Se 0.08333333333333 0.16666666666667 0.64796074482516 78.960 > 8 114 Se 0.33333333333333 0.16666666666667 0.64796074482516 78.960 > 8 115 Se 0.58333333333333 0.16666666666667 0.64796074482516 78.960 > 8 116 Se 0.83333333333333 0.16666666666667 0.64796074482516 78.960 > 8 117 Se 0.08333333333333 0.41666666666667 0.64796074482516 78.960 > 8 118 Se 0.33333333333333 0.41666666666667 0.64796074482516 78.960 > 8 119 Se 0.58333333333333 0.41666666666667 0.64796074482516 78.960 > 8 120 Se 0.83333333333333 0.41666666666667 0.64796074482516 78.960 > 8 121 Se 0.08333333333333 0.66666666666667 0.64796074482516 78.960 > 8 122 Se 0.33333333333333 0.66666666666667 0.64796074482516 78.960 > 8 123 Se 0.58333333333333 0.66666666666667 0.64796074482516 78.960 > 8 124 Se 0.83333333333333 0.66666666666667 0.64796074482516 78.960 > 8 125 Se 0.08333333333333 0.91666666666667 0.64796074482516 78.960 > 8 126 Se 0.33333333333333 0.91666666666667 0.64796074482516 78.960 > 8 127 Se 0.58333333333333 0.91666666666667 0.64796074482516 78.960 > 8 128 Se 0.83333333333333 0.91666666666667 0.64796074482516 78.960 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.4174257 -0.0000000 0.0000000 0.0000000 9.4174257 0.0000000 0.0000000 0.0000000 9.7056307 -------------------------- Born effective charges -------------------------- 1 Ga 2.3421441 -0.0000000 0.0000000 0.0000000 2.3421441 0.0000000 0.0000000 0.0000000 0.8036420 2 Ga 2.3421441 -0.0000000 0.0000000 0.0000000 2.3421441 0.0000000 0.0000000 0.0000000 0.8036420 3 Ga 2.3263402 0.0000000 0.0000000 0.0000000 2.3263402 0.0000000 0.0000000 0.0000000 0.8121901 4 Ga 2.3263402 0.0000000 0.0000000 0.0000000 2.3263402 0.0000000 0.0000000 0.0000000 0.8121901 5 Se -2.3276685 -0.0000000 0.0000000 0.0000000 -2.3276685 0.0000000 0.0000000 0.0000000 -0.7843259 6 Se -2.3276685 -0.0000000 0.0000000 0.0000000 -2.3276685 0.0000000 0.0000000 0.0000000 -0.7843259 7 Se -2.3408158 0.0000000 0.0000000 0.0000000 -2.3408158 0.0000000 0.0000000 0.0000000 -0.8315062 8 Se -2.3408158 0.0000000 0.0000000 0.0000000 -2.3408158 0.0000000 0.0000000 0.0000000 -0.8315062 ---------------------------------------------------------------------------- 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: 448 Number of blocks in projector: 312 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: 196 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 160 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 92 Use standard eigh solver. Tree of FC basis block matrices: - (448, 440), data: False |-- (92, 88), data: True |-- (160, 160), data: True |-- (196, 192), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.010 Solver_block: 80 / 80 - Time: 0.336 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.349 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) 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: 448 Number of blocks in projector: 312 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: 196 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 160 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 92 Use standard eigh solver. Tree of FC basis block matrices: - (448, 440), data: False |-- (92, 88), data: True |-- (160, 160), data: True |-- (196, 192), data: True Max drift after symmetrization by symfc projector: 0.00000000 (zz) 0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 18:53:08]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:53:08]------------------------- 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: P-6m2 (187) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.747820461219762 0.000000000000000 0.000000000000000 b -1.873910230609881 3.245707728239425 0.000000000000000 c 0.000000000000000 0.000000000000000 16.122787270000000 Atomic positions (fractional): 1 Ga 0.00000000000000 0.00000000000000 0.92470939142131 69.723 2 Ga 0.00000000000000 0.00000000000000 0.07529060857869 69.723 3 Ga 0.66666666666667 0.33333333333333 0.57525982534184 69.723 4 Ga 0.66666666666667 0.33333333333333 0.42474017465816 69.723 5 Se 0.66666666666667 0.33333333333333 0.14798852090080 78.960 6 Se 0.66666666666667 0.33333333333333 0.85201147909920 78.960 7 Se 0.33333333333333 0.66666666666667 0.35203925517484 78.960 8 Se 0.33333333333333 0.66666666666667 0.64796074482516 78.960 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.991281844879047 0.000000000000000 0.000000000000000 b -7.495640922439524 12.982830912957700 0.000000000000000 c 0.000000000000000 0.000000000000000 16.122787270000000 Atomic positions (fractional): 1 Ga 0.00000000000000 0.00000000000000 0.92470939142131 69.723 > 1 2 Ga 0.25000000000000 0.00000000000000 0.92470939142131 69.723 > 1 3 Ga 0.50000000000000 0.00000000000000 0.92470939142131 69.723 > 1 4 Ga 0.75000000000000 0.00000000000000 0.92470939142131 69.723 > 1 5 Ga 0.00000000000000 0.25000000000000 0.92470939142131 69.723 > 1 6 Ga 0.25000000000000 0.25000000000000 0.92470939142131 69.723 > 1 7 Ga 0.50000000000000 0.25000000000000 0.92470939142131 69.723 > 1 8 Ga 0.75000000000000 0.25000000000000 0.92470939142131 69.723 > 1 9 Ga 0.00000000000000 0.50000000000000 0.92470939142131 69.723 > 1 10 Ga 0.25000000000000 0.50000000000000 0.92470939142131 69.723 > 1 11 Ga 0.50000000000000 0.50000000000000 0.92470939142131 69.723 > 1 12 Ga 0.75000000000000 0.50000000000000 0.92470939142131 69.723 > 1 13 Ga 0.00000000000000 0.75000000000000 0.92470939142131 69.723 > 1 14 Ga 0.25000000000000 0.75000000000000 0.92470939142131 69.723 > 1 15 Ga 0.50000000000000 0.75000000000000 0.92470939142131 69.723 > 1 16 Ga 0.75000000000000 0.75000000000000 0.92470939142131 69.723 > 1 17 Ga 0.00000000000000 0.00000000000000 0.07529060857869 69.723 > 17 18 Ga 0.25000000000000 0.00000000000000 0.07529060857869 69.723 > 17 19 Ga 0.50000000000000 0.00000000000000 0.07529060857869 69.723 > 17 20 Ga 0.75000000000000 0.00000000000000 0.07529060857869 69.723 > 17 21 Ga 0.00000000000000 0.25000000000000 0.07529060857869 69.723 > 17 22 Ga 0.25000000000000 0.25000000000000 0.07529060857869 69.723 > 17 23 Ga 0.50000000000000 0.25000000000000 0.07529060857869 69.723 > 17 24 Ga 0.75000000000000 0.25000000000000 0.07529060857869 69.723 > 17 25 Ga 0.00000000000000 0.50000000000000 0.07529060857869 69.723 > 17 26 Ga 0.25000000000000 0.50000000000000 0.07529060857869 69.723 > 17 27 Ga 0.50000000000000 0.50000000000000 0.07529060857869 69.723 > 17 28 Ga 0.75000000000000 0.50000000000000 0.07529060857869 69.723 > 17 29 Ga 0.00000000000000 0.75000000000000 0.07529060857869 69.723 > 17 30 Ga 0.25000000000000 0.75000000000000 0.07529060857869 69.723 > 17 31 Ga 0.50000000000000 0.75000000000000 0.07529060857869 69.723 > 17 32 Ga 0.75000000000000 0.75000000000000 0.07529060857869 69.723 > 17 33 Ga 0.16666666666667 0.08333333333333 0.57525982534184 69.723 > 33 34 Ga 0.41666666666667 0.08333333333333 0.57525982534184 69.723 > 33 35 Ga 0.66666666666667 0.08333333333333 0.57525982534184 69.723 > 33 36 Ga 0.91666666666667 0.08333333333333 0.57525982534184 69.723 > 33 37 Ga 0.16666666666667 0.33333333333333 0.57525982534184 69.723 > 33 38 Ga 0.41666666666667 0.33333333333333 0.57525982534184 69.723 > 33 39 Ga 0.66666666666667 0.33333333333333 0.57525982534184 69.723 > 33 40 Ga 0.91666666666667 0.33333333333333 0.57525982534184 69.723 > 33 41 Ga 0.16666666666667 0.58333333333333 0.57525982534184 69.723 > 33 42 Ga 0.41666666666667 0.58333333333333 0.57525982534184 69.723 > 33 43 Ga 0.66666666666667 0.58333333333333 0.57525982534184 69.723 > 33 44 Ga 0.91666666666667 0.58333333333333 0.57525982534184 69.723 > 33 45 Ga 0.16666666666667 0.83333333333333 0.57525982534184 69.723 > 33 46 Ga 0.41666666666667 0.83333333333333 0.57525982534184 69.723 > 33 47 Ga 0.66666666666667 0.83333333333333 0.57525982534184 69.723 > 33 48 Ga 0.91666666666667 0.83333333333333 0.57525982534184 69.723 > 33 49 Ga 0.16666666666667 0.08333333333333 0.42474017465816 69.723 > 49 50 Ga 0.41666666666667 0.08333333333333 0.42474017465816 69.723 > 49 51 Ga 0.66666666666667 0.08333333333333 0.42474017465816 69.723 > 49 52 Ga 0.91666666666667 0.08333333333333 0.42474017465816 69.723 > 49 53 Ga 0.16666666666667 0.33333333333333 0.42474017465816 69.723 > 49 54 Ga 0.41666666666667 0.33333333333333 0.42474017465816 69.723 > 49 55 Ga 0.66666666666667 0.33333333333333 0.42474017465816 69.723 > 49 56 Ga 0.91666666666667 0.33333333333333 0.42474017465816 69.723 > 49 57 Ga 0.16666666666667 0.58333333333333 0.42474017465816 69.723 > 49 58 Ga 0.41666666666667 0.58333333333333 0.42474017465816 69.723 > 49 59 Ga 0.66666666666667 0.58333333333333 0.42474017465816 69.723 > 49 60 Ga 0.91666666666667 0.58333333333333 0.42474017465816 69.723 > 49 61 Ga 0.16666666666667 0.83333333333333 0.42474017465816 69.723 > 49 62 Ga 0.41666666666667 0.83333333333333 0.42474017465816 69.723 > 49 63 Ga 0.66666666666667 0.83333333333333 0.42474017465816 69.723 > 49 64 Ga 0.91666666666667 0.83333333333333 0.42474017465816 69.723 > 49 65 Se 0.16666666666667 0.08333333333333 0.14798852090080 78.960 > 65 66 Se 0.41666666666667 0.08333333333333 0.14798852090080 78.960 > 65 67 Se 0.66666666666667 0.08333333333333 0.14798852090080 78.960 > 65 68 Se 0.91666666666667 0.08333333333333 0.14798852090080 78.960 > 65 69 Se 0.16666666666667 0.33333333333333 0.14798852090080 78.960 > 65 70 Se 0.41666666666667 0.33333333333333 0.14798852090080 78.960 > 65 71 Se 0.66666666666667 0.33333333333333 0.14798852090080 78.960 > 65 72 Se 0.91666666666667 0.33333333333333 0.14798852090080 78.960 > 65 73 Se 0.16666666666667 0.58333333333333 0.14798852090080 78.960 > 65 74 Se 0.41666666666667 0.58333333333333 0.14798852090080 78.960 > 65 75 Se 0.66666666666667 0.58333333333333 0.14798852090080 78.960 > 65 76 Se 0.91666666666667 0.58333333333333 0.14798852090080 78.960 > 65 77 Se 0.16666666666667 0.83333333333333 0.14798852090080 78.960 > 65 78 Se 0.41666666666667 0.83333333333333 0.14798852090080 78.960 > 65 79 Se 0.66666666666667 0.83333333333333 0.14798852090080 78.960 > 65 80 Se 0.91666666666667 0.83333333333333 0.14798852090080 78.960 > 65 81 Se 0.16666666666667 0.08333333333333 0.85201147909920 78.960 > 81 82 Se 0.41666666666667 0.08333333333333 0.85201147909920 78.960 > 81 83 Se 0.66666666666667 0.08333333333333 0.85201147909920 78.960 > 81 84 Se 0.91666666666667 0.08333333333333 0.85201147909920 78.960 > 81 85 Se 0.16666666666667 0.33333333333333 0.85201147909920 78.960 > 81 86 Se 0.41666666666667 0.33333333333333 0.85201147909920 78.960 > 81 87 Se 0.66666666666667 0.33333333333333 0.85201147909920 78.960 > 81 88 Se 0.91666666666667 0.33333333333333 0.85201147909920 78.960 > 81 89 Se 0.16666666666667 0.58333333333333 0.85201147909920 78.960 > 81 90 Se 0.41666666666667 0.58333333333333 0.85201147909920 78.960 > 81 91 Se 0.66666666666667 0.58333333333333 0.85201147909920 78.960 > 81 92 Se 0.91666666666667 0.58333333333333 0.85201147909920 78.960 > 81 93 Se 0.16666666666667 0.83333333333333 0.85201147909920 78.960 > 81 94 Se 0.41666666666667 0.83333333333333 0.85201147909920 78.960 > 81 95 Se 0.66666666666667 0.83333333333333 0.85201147909920 78.960 > 81 96 Se 0.91666666666667 0.83333333333333 0.85201147909920 78.960 > 81 97 Se 0.08333333333333 0.16666666666667 0.35203925517484 78.960 > 97 98 Se 0.33333333333333 0.16666666666667 0.35203925517484 78.960 > 97 99 Se 0.58333333333333 0.16666666666667 0.35203925517484 78.960 > 97 100 Se 0.83333333333333 0.16666666666667 0.35203925517484 78.960 > 97 101 Se 0.08333333333333 0.41666666666667 0.35203925517484 78.960 > 97 102 Se 0.33333333333333 0.41666666666667 0.35203925517484 78.960 > 97 103 Se 0.58333333333333 0.41666666666667 0.35203925517484 78.960 > 97 104 Se 0.83333333333333 0.41666666666667 0.35203925517484 78.960 > 97 105 Se 0.08333333333333 0.66666666666667 0.35203925517484 78.960 > 97 106 Se 0.33333333333333 0.66666666666667 0.35203925517484 78.960 > 97 107 Se 0.58333333333333 0.66666666666667 0.35203925517484 78.960 > 97 108 Se 0.83333333333333 0.66666666666667 0.35203925517484 78.960 > 97 109 Se 0.08333333333333 0.91666666666667 0.35203925517484 78.960 > 97 110 Se 0.33333333333333 0.91666666666667 0.35203925517484 78.960 > 97 111 Se 0.58333333333333 0.91666666666667 0.35203925517484 78.960 > 97 112 Se 0.83333333333333 0.91666666666667 0.35203925517484 78.960 > 97 113 Se 0.08333333333333 0.16666666666667 0.64796074482516 78.960 > 113 114 Se 0.33333333333333 0.16666666666667 0.64796074482516 78.960 > 113 115 Se 0.58333333333333 0.16666666666667 0.64796074482516 78.960 > 113 116 Se 0.83333333333333 0.16666666666667 0.64796074482516 78.960 > 113 117 Se 0.08333333333333 0.41666666666667 0.64796074482516 78.960 > 113 118 Se 0.33333333333333 0.41666666666667 0.64796074482516 78.960 > 113 119 Se 0.58333333333333 0.41666666666667 0.64796074482516 78.960 > 113 120 Se 0.83333333333333 0.41666666666667 0.64796074482516 78.960 > 113 121 Se 0.08333333333333 0.66666666666667 0.64796074482516 78.960 > 113 122 Se 0.33333333333333 0.66666666666667 0.64796074482516 78.960 > 113 123 Se 0.58333333333333 0.66666666666667 0.64796074482516 78.960 > 113 124 Se 0.83333333333333 0.66666666666667 0.64796074482516 78.960 > 113 125 Se 0.08333333333333 0.91666666666667 0.64796074482516 78.960 > 113 126 Se 0.33333333333333 0.91666666666667 0.64796074482516 78.960 > 113 127 Se 0.58333333333333 0.91666666666667 0.64796074482516 78.960 > 113 128 Se 0.83333333333333 0.91666666666667 0.64796074482516 78.960 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.4174257 -0.0000000 0.0000000 0.0000000 9.4174257 0.0000000 0.0000000 0.0000000 9.7056307 -------------------------- Born effective charges -------------------------- 1 Ga 2.3421441 -0.0000000 0.0000000 0.0000000 2.3421441 0.0000000 0.0000000 0.0000000 0.8036420 2 Ga 2.3421441 -0.0000000 0.0000000 0.0000000 2.3421441 0.0000000 0.0000000 0.0000000 0.8036420 3 Ga 2.3263402 0.0000000 0.0000000 0.0000000 2.3263402 0.0000000 0.0000000 0.0000000 0.8121901 4 Ga 2.3263402 0.0000000 0.0000000 0.0000000 2.3263402 0.0000000 0.0000000 0.0000000 0.8121901 5 Se -2.3276685 -0.0000000 0.0000000 0.0000000 -2.3276685 0.0000000 0.0000000 0.0000000 -0.7843259 6 Se -2.3276685 -0.0000000 0.0000000 0.0000000 -2.3276685 0.0000000 0.0000000 0.0000000 -0.7843259 7 Se -2.3408158 0.0000000 0.0000000 0.0000000 -2.3408158 0.0000000 0.0000000 0.0000000 -0.8315062 8 Se -2.3408158 0.0000000 0.0000000 0.0000000 -2.3408158 0.0000000 0.0000000 0.0000000 -0.8315062 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] 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.00000002 (xzx) 0.00000002 (xzx) 0.00000002 (xxz) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) 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-07 18:53:18]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:53:19]------------------------- 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: P-6m2 (187) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.747820461219762 0.000000000000000 0.000000000000000 b -1.873910230609881 3.245707728239425 0.000000000000000 c 0.000000000000000 0.000000000000000 16.122787270000000 Atomic positions (fractional): 1 Ga 0.00000000000000 0.00000000000000 0.92470939142131 69.723 2 Ga 0.00000000000000 0.00000000000000 0.07529060857869 69.723 3 Ga 0.66666666666667 0.33333333333333 0.57525982534184 69.723 4 Ga 0.66666666666667 0.33333333333333 0.42474017465816 69.723 5 Se 0.66666666666667 0.33333333333333 0.14798852090080 78.960 6 Se 0.66666666666667 0.33333333333333 0.85201147909920 78.960 7 Se 0.33333333333333 0.66666666666667 0.35203925517484 78.960 8 Se 0.33333333333333 0.66666666666667 0.64796074482516 78.960 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.991281844879047 0.000000000000000 0.000000000000000 b -7.495640922439524 12.982830912957700 0.000000000000000 c 0.000000000000000 0.000000000000000 16.122787270000000 Atomic positions (fractional): 1 Ga 0.00000000000000 0.00000000000000 0.92470939142131 69.723 > 1 2 Ga 0.25000000000000 0.00000000000000 0.92470939142131 69.723 > 1 3 Ga 0.50000000000000 0.00000000000000 0.92470939142131 69.723 > 1 4 Ga 0.75000000000000 0.00000000000000 0.92470939142131 69.723 > 1 5 Ga 0.00000000000000 0.25000000000000 0.92470939142131 69.723 > 1 6 Ga 0.25000000000000 0.25000000000000 0.92470939142131 69.723 > 1 7 Ga 0.50000000000000 0.25000000000000 0.92470939142131 69.723 > 1 8 Ga 0.75000000000000 0.25000000000000 0.92470939142131 69.723 > 1 9 Ga 0.00000000000000 0.50000000000000 0.92470939142131 69.723 > 1 10 Ga 0.25000000000000 0.50000000000000 0.92470939142131 69.723 > 1 11 Ga 0.50000000000000 0.50000000000000 0.92470939142131 69.723 > 1 12 Ga 0.75000000000000 0.50000000000000 0.92470939142131 69.723 > 1 13 Ga 0.00000000000000 0.75000000000000 0.92470939142131 69.723 > 1 14 Ga 0.25000000000000 0.75000000000000 0.92470939142131 69.723 > 1 15 Ga 0.50000000000000 0.75000000000000 0.92470939142131 69.723 > 1 16 Ga 0.75000000000000 0.75000000000000 0.92470939142131 69.723 > 1 17 Ga 0.00000000000000 0.00000000000000 0.07529060857869 69.723 > 17 18 Ga 0.25000000000000 0.00000000000000 0.07529060857869 69.723 > 17 19 Ga 0.50000000000000 0.00000000000000 0.07529060857869 69.723 > 17 20 Ga 0.75000000000000 0.00000000000000 0.07529060857869 69.723 > 17 21 Ga 0.00000000000000 0.25000000000000 0.07529060857869 69.723 > 17 22 Ga 0.25000000000000 0.25000000000000 0.07529060857869 69.723 > 17 23 Ga 0.50000000000000 0.25000000000000 0.07529060857869 69.723 > 17 24 Ga 0.75000000000000 0.25000000000000 0.07529060857869 69.723 > 17 25 Ga 0.00000000000000 0.50000000000000 0.07529060857869 69.723 > 17 26 Ga 0.25000000000000 0.50000000000000 0.07529060857869 69.723 > 17 27 Ga 0.50000000000000 0.50000000000000 0.07529060857869 69.723 > 17 28 Ga 0.75000000000000 0.50000000000000 0.07529060857869 69.723 > 17 29 Ga 0.00000000000000 0.75000000000000 0.07529060857869 69.723 > 17 30 Ga 0.25000000000000 0.75000000000000 0.07529060857869 69.723 > 17 31 Ga 0.50000000000000 0.75000000000000 0.07529060857869 69.723 > 17 32 Ga 0.75000000000000 0.75000000000000 0.07529060857869 69.723 > 17 33 Ga 0.16666666666667 0.08333333333333 0.57525982534184 69.723 > 33 34 Ga 0.41666666666667 0.08333333333333 0.57525982534184 69.723 > 33 35 Ga 0.66666666666667 0.08333333333333 0.57525982534184 69.723 > 33 36 Ga 0.91666666666667 0.08333333333333 0.57525982534184 69.723 > 33 37 Ga 0.16666666666667 0.33333333333333 0.57525982534184 69.723 > 33 38 Ga 0.41666666666667 0.33333333333333 0.57525982534184 69.723 > 33 39 Ga 0.66666666666667 0.33333333333333 0.57525982534184 69.723 > 33 40 Ga 0.91666666666667 0.33333333333333 0.57525982534184 69.723 > 33 41 Ga 0.16666666666667 0.58333333333333 0.57525982534184 69.723 > 33 42 Ga 0.41666666666667 0.58333333333333 0.57525982534184 69.723 > 33 43 Ga 0.66666666666667 0.58333333333333 0.57525982534184 69.723 > 33 44 Ga 0.91666666666667 0.58333333333333 0.57525982534184 69.723 > 33 45 Ga 0.16666666666667 0.83333333333333 0.57525982534184 69.723 > 33 46 Ga 0.41666666666667 0.83333333333333 0.57525982534184 69.723 > 33 47 Ga 0.66666666666667 0.83333333333333 0.57525982534184 69.723 > 33 48 Ga 0.91666666666667 0.83333333333333 0.57525982534184 69.723 > 33 49 Ga 0.16666666666667 0.08333333333333 0.42474017465816 69.723 > 49 50 Ga 0.41666666666667 0.08333333333333 0.42474017465816 69.723 > 49 51 Ga 0.66666666666667 0.08333333333333 0.42474017465816 69.723 > 49 52 Ga 0.91666666666667 0.08333333333333 0.42474017465816 69.723 > 49 53 Ga 0.16666666666667 0.33333333333333 0.42474017465816 69.723 > 49 54 Ga 0.41666666666667 0.33333333333333 0.42474017465816 69.723 > 49 55 Ga 0.66666666666667 0.33333333333333 0.42474017465816 69.723 > 49 56 Ga 0.91666666666667 0.33333333333333 0.42474017465816 69.723 > 49 57 Ga 0.16666666666667 0.58333333333333 0.42474017465816 69.723 > 49 58 Ga 0.41666666666667 0.58333333333333 0.42474017465816 69.723 > 49 59 Ga 0.66666666666667 0.58333333333333 0.42474017465816 69.723 > 49 60 Ga 0.91666666666667 0.58333333333333 0.42474017465816 69.723 > 49 61 Ga 0.16666666666667 0.83333333333333 0.42474017465816 69.723 > 49 62 Ga 0.41666666666667 0.83333333333333 0.42474017465816 69.723 > 49 63 Ga 0.66666666666667 0.83333333333333 0.42474017465816 69.723 > 49 64 Ga 0.91666666666667 0.83333333333333 0.42474017465816 69.723 > 49 65 Se 0.16666666666667 0.08333333333333 0.14798852090080 78.960 > 65 66 Se 0.41666666666667 0.08333333333333 0.14798852090080 78.960 > 65 67 Se 0.66666666666667 0.08333333333333 0.14798852090080 78.960 > 65 68 Se 0.91666666666667 0.08333333333333 0.14798852090080 78.960 > 65 69 Se 0.16666666666667 0.33333333333333 0.14798852090080 78.960 > 65 70 Se 0.41666666666667 0.33333333333333 0.14798852090080 78.960 > 65 71 Se 0.66666666666667 0.33333333333333 0.14798852090080 78.960 > 65 72 Se 0.91666666666667 0.33333333333333 0.14798852090080 78.960 > 65 73 Se 0.16666666666667 0.58333333333333 0.14798852090080 78.960 > 65 74 Se 0.41666666666667 0.58333333333333 0.14798852090080 78.960 > 65 75 Se 0.66666666666667 0.58333333333333 0.14798852090080 78.960 > 65 76 Se 0.91666666666667 0.58333333333333 0.14798852090080 78.960 > 65 77 Se 0.16666666666667 0.83333333333333 0.14798852090080 78.960 > 65 78 Se 0.41666666666667 0.83333333333333 0.14798852090080 78.960 > 65 79 Se 0.66666666666667 0.83333333333333 0.14798852090080 78.960 > 65 80 Se 0.91666666666667 0.83333333333333 0.14798852090080 78.960 > 65 81 Se 0.16666666666667 0.08333333333333 0.85201147909920 78.960 > 81 82 Se 0.41666666666667 0.08333333333333 0.85201147909920 78.960 > 81 83 Se 0.66666666666667 0.08333333333333 0.85201147909920 78.960 > 81 84 Se 0.91666666666667 0.08333333333333 0.85201147909920 78.960 > 81 85 Se 0.16666666666667 0.33333333333333 0.85201147909920 78.960 > 81 86 Se 0.41666666666667 0.33333333333333 0.85201147909920 78.960 > 81 87 Se 0.66666666666667 0.33333333333333 0.85201147909920 78.960 > 81 88 Se 0.91666666666667 0.33333333333333 0.85201147909920 78.960 > 81 89 Se 0.16666666666667 0.58333333333333 0.85201147909920 78.960 > 81 90 Se 0.41666666666667 0.58333333333333 0.85201147909920 78.960 > 81 91 Se 0.66666666666667 0.58333333333333 0.85201147909920 78.960 > 81 92 Se 0.91666666666667 0.58333333333333 0.85201147909920 78.960 > 81 93 Se 0.16666666666667 0.83333333333333 0.85201147909920 78.960 > 81 94 Se 0.41666666666667 0.83333333333333 0.85201147909920 78.960 > 81 95 Se 0.66666666666667 0.83333333333333 0.85201147909920 78.960 > 81 96 Se 0.91666666666667 0.83333333333333 0.85201147909920 78.960 > 81 97 Se 0.08333333333333 0.16666666666667 0.35203925517484 78.960 > 97 98 Se 0.33333333333333 0.16666666666667 0.35203925517484 78.960 > 97 99 Se 0.58333333333333 0.16666666666667 0.35203925517484 78.960 > 97 100 Se 0.83333333333333 0.16666666666667 0.35203925517484 78.960 > 97 101 Se 0.08333333333333 0.41666666666667 0.35203925517484 78.960 > 97 102 Se 0.33333333333333 0.41666666666667 0.35203925517484 78.960 > 97 103 Se 0.58333333333333 0.41666666666667 0.35203925517484 78.960 > 97 104 Se 0.83333333333333 0.41666666666667 0.35203925517484 78.960 > 97 105 Se 0.08333333333333 0.66666666666667 0.35203925517484 78.960 > 97 106 Se 0.33333333333333 0.66666666666667 0.35203925517484 78.960 > 97 107 Se 0.58333333333333 0.66666666666667 0.35203925517484 78.960 > 97 108 Se 0.83333333333333 0.66666666666667 0.35203925517484 78.960 > 97 109 Se 0.08333333333333 0.91666666666667 0.35203925517484 78.960 > 97 110 Se 0.33333333333333 0.91666666666667 0.35203925517484 78.960 > 97 111 Se 0.58333333333333 0.91666666666667 0.35203925517484 78.960 > 97 112 Se 0.83333333333333 0.91666666666667 0.35203925517484 78.960 > 97 113 Se 0.08333333333333 0.16666666666667 0.64796074482516 78.960 > 113 114 Se 0.33333333333333 0.16666666666667 0.64796074482516 78.960 > 113 115 Se 0.58333333333333 0.16666666666667 0.64796074482516 78.960 > 113 116 Se 0.83333333333333 0.16666666666667 0.64796074482516 78.960 > 113 117 Se 0.08333333333333 0.41666666666667 0.64796074482516 78.960 > 113 118 Se 0.33333333333333 0.41666666666667 0.64796074482516 78.960 > 113 119 Se 0.58333333333333 0.41666666666667 0.64796074482516 78.960 > 113 120 Se 0.83333333333333 0.41666666666667 0.64796074482516 78.960 > 113 121 Se 0.08333333333333 0.66666666666667 0.64796074482516 78.960 > 113 122 Se 0.33333333333333 0.66666666666667 0.64796074482516 78.960 > 113 123 Se 0.58333333333333 0.66666666666667 0.64796074482516 78.960 > 113 124 Se 0.83333333333333 0.66666666666667 0.64796074482516 78.960 > 113 125 Se 0.08333333333333 0.91666666666667 0.64796074482516 78.960 > 113 126 Se 0.33333333333333 0.91666666666667 0.64796074482516 78.960 > 113 127 Se 0.58333333333333 0.91666666666667 0.64796074482516 78.960 > 113 128 Se 0.83333333333333 0.91666666666667 0.64796074482516 78.960 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 9.4174257 -0.0000000 0.0000000 0.0000000 9.4174257 0.0000000 0.0000000 0.0000000 9.7056307 -------------------------- Born effective charges -------------------------- 1 Ga 2.3421441 -0.0000000 0.0000000 0.0000000 2.3421441 0.0000000 0.0000000 0.0000000 0.8036420 2 Ga 2.3421441 -0.0000000 0.0000000 0.0000000 2.3421441 0.0000000 0.0000000 0.0000000 0.8036420 3 Ga 2.3263402 0.0000000 0.0000000 0.0000000 2.3263402 0.0000000 0.0000000 0.0000000 0.8121901 4 Ga 2.3263402 0.0000000 0.0000000 0.0000000 2.3263402 0.0000000 0.0000000 0.0000000 0.8121901 5 Se -2.3276685 -0.0000000 0.0000000 0.0000000 -2.3276685 0.0000000 0.0000000 0.0000000 -0.7843259 6 Se -2.3276685 -0.0000000 0.0000000 0.0000000 -2.3276685 0.0000000 0.0000000 0.0000000 -0.7843259 7 Se -2.3408158 0.0000000 0.0000000 0.0000000 -2.3408158 0.0000000 0.0000000 0.0000000 -0.8315062 8 Se -2.3408158 0.0000000 0.0000000 0.0000000 -2.3408158 0.0000000 0.0000000 0.0000000 -0.8315062 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000002 (xzx) 0.00000002 (xzx) 0.00000002 (xxz) Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 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.71, Number of G-points: 305, Lambda: 0.23 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.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 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 0.493 ( 0.000 0.000 0.000) 0.000 0.493 ( 0.000 -0.000 0.000) 0.000 0.765 ( 0.000 -0.000 0.000) 0.000 1.609 ( 0.000 0.000 0.000) 0.000 1.609 ( 0.000 -0.000 0.000) 0.000 1.642 ( 0.000 0.000 0.000) 0.000 1.642 ( 0.000 -0.000 0.000) 0.000 3.768 ( 0.000 0.000 0.000) 0.000 3.854 ( 0.000 0.000 0.000) 0.000 6.094 ( 0.000 0.000 0.000) 0.000 6.094 ( 0.000 -0.000 0.000) 0.000 6.119 ( 0.000 0.000 0.000) 0.000 6.119 ( 0.000 -0.000 0.000) 0.000 6.266 ( 0.000 0.000 0.000) 0.000 6.266 ( 0.000 -0.000 0.000) 0.000 6.300 ( 0.000 0.000 0.000) 0.000 6.300 ( 0.000 -0.000 0.000) 0.000 7.066 ( 0.000 0.000 0.000) 0.000 7.264 ( 0.000 0.000 0.000) 0.000 9.121 ( 0.000 0.000 0.000) 0.000 9.134 ( 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.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.266 ( 13.918 8.035 0.000) 16.071 0.465 ( 19.282 11.133 0.000) 22.266 0.688 ( 13.760 7.944 0.000) 15.888 0.746 ( 0.784 0.453 0.000) 0.905 0.882 ( 35.806 20.672 0.000) 41.345 0.986 ( 29.810 17.211 0.000) 34.422 1.689 ( 6.413 3.703 0.000) 7.405 1.719 ( 6.160 3.556 0.000) 7.113 1.886 ( 18.843 10.879 0.000) 21.758 1.911 ( 21.372 12.339 0.000) 24.678 3.735 ( -2.373 -1.370 0.000) 2.740 3.828 ( -2.076 -1.198 0.000) 2.397 6.084 ( -0.890 -0.514 0.000) 1.028 6.109 ( -0.811 -0.468 0.000) 0.936 6.123 ( 1.276 0.737 0.000) 1.474 6.140 ( 1.221 0.705 0.000) 1.409 6.256 ( -0.767 -0.443 0.000) 0.885 6.291 ( -0.734 -0.424 0.000) 0.848 6.451 ( 11.470 6.622 0.000) 13.245 7.008 ( -4.113 -2.375 0.000) 4.749 7.199 ( -4.760 -2.748 0.000) 5.496 7.303 ( -0.265 -0.153 0.000) 0.306 9.088 ( -2.658 -1.534 0.000) 3.069 9.100 ( -2.767 -1.598 0.000) 3.195 ======================= Grid point 2 (3/54) ======================= q-point: ( 0.13 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.637 ( 16.676 9.628 0.000) 19.256 0.856 ( 8.283 4.782 0.000) 9.565 0.905 ( 17.636 10.182 0.000) 20.365 1.047 ( 15.388 8.884 0.000) 17.769 1.651 ( 26.954 15.562 0.000) 31.124 1.679 ( 27.918 16.118 0.000) 32.236 1.889 ( 9.846 5.684 0.000) 11.369 1.913 ( 9.668 5.582 0.000) 11.163 2.408 ( 23.331 13.470 0.000) 26.940 2.487 ( 24.933 14.395 0.000) 28.790 3.686 ( -0.877 -0.506 0.000) 1.012 3.771 ( -2.128 -1.228 0.000) 2.457 6.055 ( -1.483 -0.856 0.000) 1.713 6.082 ( -1.378 -0.796 0.000) 1.591 6.097 ( -4.219 -2.436 0.000) 4.872 6.115 ( -4.496 -2.596 0.000) 5.192 6.231 ( -1.290 -0.745 0.000) 1.490 6.266 ( -1.302 -0.752 0.000) 1.503 6.756 ( 12.629 7.291 0.000) 14.583 6.929 ( -1.333 -0.769 0.000) 1.539 7.099 ( -2.419 -1.396 0.000) 2.793 7.290 ( -0.918 -0.530 0.000) 1.060 9.005 ( -4.043 -2.334 0.000) 4.668 9.013 ( -4.229 -2.442 0.000) 4.883 ======================= Grid point 3 (4/54) ======================= q-point: ( 0.20 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.029 ( 15.696 9.062 0.000) 18.124 1.116 ( 13.287 7.671 0.000) 15.343 1.288 ( 14.391 8.309 0.000) 16.618 1.387 ( 12.929 7.464 0.000) 14.929 2.125 ( 12.842 7.414 0.000) 14.829 2.128 ( 9.835 5.678 0.000) 11.357 2.152 ( 9.981 5.763 0.000) 11.525 2.229 ( 16.532 9.545 0.000) 19.090 2.923 ( 18.977 10.957 0.000) 21.913 3.023 ( 19.255 11.117 0.000) 22.233 3.738 ( 5.617 3.243 0.000) 6.486 3.771 ( 3.664 2.115 0.000) 4.230 5.939 ( -7.673 -4.430 0.000) 8.860 5.953 ( -7.073 -4.083 0.000) 8.167 6.018 ( -1.524 -0.880 0.000) 1.760 6.047 ( -1.465 -0.846 0.000) 1.692 6.199 ( -1.346 -0.777 0.000) 1.554 6.233 ( -1.432 -0.827 0.000) 1.653 6.955 ( 2.730 1.576 0.000) 3.152 6.998 ( 7.257 4.190 0.000) 8.380 7.096 ( 1.476 0.852 0.000) 1.704 7.255 ( -2.146 -1.239 0.000) 2.478 8.909 ( -3.798 -2.193 0.000) 4.386 8.917 ( -3.499 -2.020 0.000) 4.040 ======================= Grid point 4 (5/54) ======================= q-point: ( 0.27 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.366 ( 13.046 7.532 0.000) 15.064 1.463 ( 14.859 8.579 0.000) 17.157 1.575 ( 9.587 5.535 0.000) 11.071 1.648 ( 8.932 5.157 0.000) 10.314 2.266 ( 0.071 0.041 0.000) 0.082 2.338 ( 7.563 4.367 0.000) 8.733 2.366 ( 7.786 4.495 0.000) 8.991 2.432 ( 1.253 0.723 0.000) 1.447 3.275 ( 10.528 6.079 0.000) 12.157 3.369 ( 9.815 5.667 0.000) 11.334 3.929 ( 11.065 6.388 0.000) 12.776 3.997 ( 13.429 7.753 0.000) 15.507 5.822 ( -0.926 -0.535 0.000) 1.069 5.856 ( -0.261 -0.151 0.000) 0.302 5.987 ( -1.014 -0.586 0.000) 1.171 6.018 ( -0.928 -0.536 0.000) 1.071 6.172 ( -0.863 -0.498 0.000) 0.997 6.204 ( -0.925 -0.534 0.000) 1.068 7.017 ( 1.944 1.122 0.000) 2.244 7.091 ( 0.640 0.370 0.000) 0.739 7.133 ( 1.082 0.624 0.000) 1.249 7.186 ( -3.640 -2.102 0.000) 4.203 8.838 ( -2.089 -1.206 0.000) 2.412 8.856 ( -1.580 -0.912 0.000) 1.825 ======================= Grid point 5 (6/54) ======================= q-point: ( 0.33 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.653 ( 10.986 6.343 0.000) 12.685 1.741 ( 4.509 2.603 0.000) 5.207 1.782 ( 11.629 6.714 0.000) 13.428 1.808 ( 4.570 2.638 0.000) 5.277 2.192 ( -5.147 -2.971 0.000) 5.943 2.356 ( -6.177 -3.566 0.000) 7.133 2.481 ( 4.473 2.583 0.000) 5.165 2.513 ( 4.542 2.622 0.000) 5.245 3.444 ( 4.579 2.644 0.000) 5.288 3.519 ( 3.852 2.224 0.000) 4.448 4.228 ( 13.043 7.531 0.000) 15.061 4.332 ( 13.742 7.934 0.000) 15.868 5.900 ( 6.718 3.879 0.000) 7.758 5.944 ( 6.996 4.039 0.000) 8.078 5.971 ( -0.290 -0.168 0.000) 0.335 6.006 ( -0.059 -0.034 0.000) 0.068 6.160 ( -0.105 -0.061 0.000) 0.121 6.192 ( -0.038 -0.022 0.000) 0.044 7.033 ( -0.660 -0.381 0.000) 0.762 7.047 ( -3.772 -2.178 0.000) 4.356 7.086 ( -4.631 -2.674 0.000) 5.347 7.134 ( -1.078 -0.622 0.000) 1.245 8.811 ( -0.288 -0.166 0.000) 0.332 8.839 ( -0.089 -0.051 0.000) 0.103 ======================= Grid point 6 (7/54) ======================= q-point: ( 0.40 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.804 ( 1.195 0.690 0.000) 1.380 1.872 ( 1.086 0.627 0.000) 1.254 1.878 ( 7.863 4.540 0.000) 9.079 2.001 ( 6.644 3.836 0.000) 7.672 2.064 ( -4.937 -2.850 0.000) 5.701 2.194 ( -6.597 -3.809 0.000) 7.618 2.555 ( 1.945 1.123 0.000) 2.245 2.586 ( 1.798 1.038 0.000) 2.076 3.527 ( 2.722 1.572 0.000) 3.143 3.589 ( 2.416 1.395 0.000) 2.789 4.511 ( 10.313 5.954 0.000) 11.908 4.612 ( 9.332 5.388 0.000) 10.775 5.971 ( 0.197 0.114 0.000) 0.228 6.012 ( 0.433 0.250 0.000) 0.500 6.091 ( 8.354 4.823 0.000) 9.646 6.138 ( 8.209 4.740 0.000) 9.479 6.164 ( 0.384 0.222 0.000) 0.444 6.199 ( 0.515 0.297 0.000) 0.595 6.942 ( -4.399 -2.540 0.000) 5.079 6.976 ( -4.323 -2.496 0.000) 4.992 6.994 ( -2.344 -1.353 0.000) 2.707 7.085 ( -2.858 -1.650 0.000) 3.300 8.817 ( 0.610 0.352 0.000) 0.704 8.843 ( 0.268 0.155 0.000) 0.309 ======================= Grid point 7 (8/54) ======================= q-point: ( 0.47 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.816 ( 0.108 0.062 0.000) 0.125 1.878 ( -0.121 -0.070 0.000) 0.139 1.980 ( -1.870 -1.080 0.000) 2.160 2.014 ( 3.230 1.865 0.000) 3.730 2.076 ( -2.830 -1.634 0.000) 3.268 2.098 ( 1.770 1.022 0.000) 2.044 2.582 ( 0.494 0.285 0.000) 0.570 2.609 ( 0.354 0.204 0.000) 0.409 3.574 ( 1.132 0.654 0.000) 1.308 3.632 ( 1.050 0.606 0.000) 1.213 4.693 ( 4.481 2.587 0.000) 5.174 4.757 ( 2.792 1.612 0.000) 3.224 5.976 ( 0.160 0.092 0.000) 0.185 6.021 ( 0.248 0.143 0.000) 0.287 6.173 ( 0.241 0.139 0.000) 0.278 6.210 ( 0.303 0.175 0.000) 0.350 6.253 ( 4.319 2.494 0.000) 4.987 6.285 ( 3.359 1.939 0.000) 3.879 6.864 ( -1.796 -1.037 0.000) 2.074 6.897 ( -1.944 -1.123 0.000) 2.245 6.944 ( -1.287 -0.743 0.000) 1.486 7.018 ( -2.055 -1.187 0.000) 2.373 8.832 ( 0.497 0.287 0.000) 0.573 8.846 ( -0.015 -0.008 0.000) 0.017 ======================= Grid point 16 (9/54) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.529 ( 9.239 16.002 0.000) 18.477 0.805 ( 3.455 5.985 0.000) 6.911 0.869 ( 12.220 21.166 0.000) 24.440 0.982 ( 10.525 18.230 0.000) 21.050 1.430 ( 17.203 29.796 0.000) 34.405 1.469 ( 15.056 26.077 0.000) 30.112 1.849 ( 6.413 11.108 0.000) 12.826 1.872 ( 6.161 10.670 0.000) 12.321 2.246 ( 12.701 21.999 0.000) 25.402 2.312 ( 13.831 23.956 0.000) 27.662 3.694 ( -1.104 -1.913 0.000) 2.208 3.785 ( -1.474 -2.553 0.000) 2.948 6.054 ( -1.455 -2.521 0.000) 2.911 6.078 ( -1.553 -2.690 0.000) 3.106 6.132 ( -0.348 -0.602 0.000) 0.695 6.151 ( -0.419 -0.725 0.000) 0.837 6.244 ( -0.421 -0.729 0.000) 0.841 6.278 ( -0.452 -0.783 0.000) 0.905 6.670 ( 7.464 12.928 0.000) 14.928 6.939 ( -1.820 -3.153 0.000) 3.641 7.116 ( -2.390 -4.140 0.000) 4.780 7.292 ( -0.530 -0.917 0.000) 1.059 9.030 ( -2.235 -3.871 0.000) 4.470 9.039 ( -2.373 -4.110 0.000) 4.746 ======================= Grid point 17 (10/54) ======================= q-point: ( 0.13 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.883 ( 12.708 13.595 0.000) 18.610 0.991 ( 8.902 8.149 0.000) 12.069 1.272 ( 8.017 22.030 0.000) 23.444 1.357 ( 8.634 19.329 0.000) 21.170 1.948 ( 15.248 11.964 0.000) 19.381 2.006 ( 19.130 14.564 0.000) 24.043 2.105 ( 6.583 14.508 0.000) 15.931 2.121 ( 6.731 13.965 0.000) 15.503 2.717 ( 17.026 15.285 0.000) 22.880 2.812 ( 17.598 15.919 0.000) 23.729 3.694 ( 2.239 1.869 0.000) 2.917 3.749 ( 0.230 -0.345 0.000) 0.415 5.967 ( -3.483 -6.316 0.000) 7.213 5.985 ( -3.749 -6.269 0.000) 7.305 6.088 ( -3.636 -0.026 0.000) 3.636 6.105 ( -3.156 0.244 0.000) 3.166 6.226 ( -1.285 0.286 0.000) 1.316 6.259 ( -1.260 0.095 0.000) 1.263 6.918 ( 7.407 7.501 0.000) 10.542 6.919 ( 1.512 -0.079 0.000) 1.514 7.070 ( 0.554 -1.247 0.000) 1.364 7.263 ( -1.110 -2.117 0.000) 2.390 8.941 ( -3.119 -3.719 0.000) 4.854 8.945 ( -3.008 -3.949 0.000) 4.964 ======================= Grid point 18 (11/54) ======================= q-point: ( 0.20 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.224 ( 11.693 9.514 0.000) 15.075 1.285 ( 14.080 8.784 0.000) 16.595 1.596 ( 3.959 19.649 0.000) 20.044 1.661 ( 3.904 17.545 0.000) 17.974 2.208 ( 6.623 0.727 0.000) 6.663 2.350 ( 9.008 2.034 0.000) 9.235 2.354 ( 3.576 15.409 0.000) 15.819 2.371 ( 4.114 14.830 0.000) 15.390 3.118 ( 14.046 7.581 0.000) 15.961 3.217 ( 13.728 7.338 0.000) 15.566 3.819 ( 7.165 4.696 0.000) 8.566 3.853 ( 9.824 5.762 0.000) 11.389 5.841 ( -2.902 -3.578 0.000) 4.607 5.865 ( -2.350 -3.201 0.000) 3.971 6.033 ( -3.175 1.333 0.000) 3.444 6.060 ( -2.835 1.333 0.000) 3.133 6.207 ( -1.859 1.540 0.000) 2.414 6.238 ( -1.840 1.320 0.000) 2.264 6.966 ( 3.297 -0.688 0.000) 3.368 7.061 ( 3.173 1.470 0.000) 3.497 7.089 ( 2.447 -1.689 0.000) 2.973 7.207 ( -2.096 -3.386 0.000) 3.982 8.855 ( -2.299 -2.989 0.000) 3.771 8.866 ( -1.758 -2.769 0.000) 3.280 ======================= Grid point 19 (12/54) ======================= q-point: ( 0.27 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.514 ( 11.506 6.932 0.000) 13.433 1.620 ( 13.590 6.831 0.000) 15.210 1.804 ( -0.705 14.670 0.000) 14.687 1.850 ( 0.003 12.617 0.000) 12.617 2.236 ( -1.591 -1.297 0.000) 2.052 2.391 ( -1.523 -3.895 0.000) 4.183 2.539 ( -0.082 14.299 0.000) 14.299 2.559 ( 0.489 13.508 0.000) 13.517 3.364 ( 7.924 2.918 0.000) 8.444 3.447 ( 7.142 2.427 0.000) 7.543 4.062 ( 12.717 6.265 0.000) 14.177 4.147 ( 14.482 6.416 0.000) 15.840 5.837 ( 3.485 1.865 0.000) 3.953 5.872 ( 3.931 1.674 0.000) 4.273 6.004 ( -1.878 2.142 0.000) 2.848 6.036 ( -1.702 2.182 0.000) 2.767 6.195 ( -1.822 2.699 0.000) 3.256 6.225 ( -1.735 2.522 0.000) 3.061 6.998 ( 2.105 -3.089 0.000) 3.738 7.065 ( -1.422 -2.870 0.000) 3.203 7.101 ( 1.706 -3.733 0.000) 4.104 7.119 ( -3.222 -4.288 0.000) 5.363 8.805 ( -0.434 -1.928 0.000) 1.976 8.828 ( 0.011 -1.794 0.000) 1.794 ======================= Grid point 20 (13/54) ======================= q-point: ( 0.33 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.768 ( 9.727 5.181 0.000) 11.021 1.875 ( -2.773 6.257 0.000) 6.844 1.889 ( 9.975 4.066 0.000) 10.772 1.925 ( -2.974 7.377 0.000) 7.954 2.173 ( -6.142 4.045 0.000) 7.354 2.279 ( -6.119 -2.219 0.000) 6.509 2.646 ( -3.242 12.615 0.000) 13.025 2.666 ( -2.744 11.692 0.000) 12.010 3.488 ( 3.977 1.923 0.000) 4.418 3.555 ( 3.459 1.657 0.000) 3.835 4.354 ( 12.817 4.619 0.000) 13.624 4.453 ( 13.155 3.729 0.000) 13.673 5.954 ( 4.270 0.811 0.000) 4.346 5.974 ( 3.630 0.654 0.000) 3.688 6.042 ( 3.068 5.295 0.000) 6.120 6.070 ( 3.510 5.773 0.000) 6.756 6.199 ( -1.491 3.722 0.000) 4.010 6.229 ( -1.295 3.527 0.000) 3.757 6.975 ( -2.691 -4.810 0.000) 5.512 6.975 ( -0.559 -5.077 0.000) 5.108 7.010 ( -3.660 -4.663 0.000) 5.928 7.068 ( 0.099 -5.669 0.000) 5.670 8.796 ( 1.038 -1.298 0.000) 1.662 8.823 ( 1.029 -1.475 0.000) 1.799 ======================= Grid point 21 (14/54) ======================= q-point: ( 0.40 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.865 ( -1.236 0.217 0.000) 1.255 1.911 ( -4.032 3.490 0.000) 5.333 1.955 ( 6.024 3.219 0.000) 6.830 2.050 ( 4.834 1.488 0.000) 5.057 2.119 ( -7.341 9.067 0.000) 11.666 2.164 ( -7.047 3.711 0.000) 7.964 2.691 ( -5.170 11.272 0.000) 12.401 2.710 ( -4.868 10.429 0.000) 11.510 3.559 ( 2.001 1.596 0.000) 2.560 3.618 ( 1.805 1.418 0.000) 2.295 4.586 ( 9.213 1.328 0.000) 9.308 4.665 ( 8.056 -0.150 0.000) 8.058 5.980 ( 0.183 0.610 0.000) 0.636 6.017 ( 0.901 0.249 0.000) 0.934 6.194 ( 5.938 4.819 0.000) 7.647 6.215 ( -1.678 4.537 0.000) 4.838 6.232 ( 5.738 4.246 0.000) 7.138 6.247 ( -1.424 4.197 0.000) 4.432 6.871 ( -2.185 -4.282 0.000) 4.807 6.906 ( -2.329 -4.301 0.000) 4.891 6.917 ( 0.223 -5.959 0.000) 5.963 6.997 ( -0.405 -6.678 0.000) 6.690 8.807 ( 1.542 -1.227 0.000) 1.971 8.827 ( 1.089 -1.666 0.000) 1.990 ======================= Grid point 22 (15/54) ======================= q-point: ( 0.47 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.855 ( 0.712 -1.234 0.000) 1.425 1.886 ( -2.129 3.687 0.000) 4.257 2.030 ( 0.054 -0.093 0.000) 0.108 2.097 ( 0.394 -0.683 0.000) 0.789 2.103 ( -6.462 11.192 0.000) 12.923 2.116 ( -4.563 7.904 0.000) 9.127 2.702 ( -6.100 10.566 0.000) 12.200 2.718 ( -5.706 9.883 0.000) 11.412 3.586 ( -0.267 0.463 0.000) 0.534 3.643 ( -0.195 0.338 0.000) 0.390 4.683 ( 2.009 -3.480 0.000) 4.018 4.732 ( 2.124 -3.679 0.000) 4.249 5.985 ( -0.320 0.554 0.000) 0.639 6.026 ( -0.079 0.137 0.000) 0.158 6.225 ( -2.633 4.560 0.000) 5.265 6.257 ( -2.364 4.095 0.000) 4.728 6.286 ( -0.390 0.675 0.000) 0.780 6.308 ( -0.227 0.394 0.000) 0.454 6.828 ( 1.389 -2.405 0.000) 2.777 6.859 ( 1.454 -2.518 0.000) 2.907 6.885 ( 2.787 -4.828 0.000) 5.574 6.949 ( 3.003 -5.201 0.000) 6.006 8.817 ( 0.941 -1.631 0.000) 1.883 8.826 ( 1.042 -1.805 0.000) 2.084 ======================= Grid point 33 (16/54) ======================= q-point: ( 0.13 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.176 ( 8.243 14.277 0.000) 16.486 1.186 ( 6.361 11.018 0.000) 12.722 1.709 ( 11.004 19.060 0.000) 22.008 1.757 ( 10.326 17.885 0.000) 20.652 2.117 ( 2.723 4.716 0.000) 5.446 2.243 ( 4.502 7.798 0.000) 9.005 2.429 ( 9.026 15.634 0.000) 18.052 2.436 ( 8.896 15.408 0.000) 17.791 3.007 ( 7.372 12.769 0.000) 14.744 3.109 ( 7.324 12.685 0.000) 14.647 3.782 ( 3.117 5.399 0.000) 6.234 3.788 ( 3.677 6.369 0.000) 7.354 5.847 ( -2.541 -4.400 0.000) 5.081 5.869 ( -2.381 -4.124 0.000) 4.762 6.079 ( -0.482 -0.834 0.000) 0.964 6.102 ( -0.298 -0.516 0.000) 0.596 6.238 ( 0.408 0.706 0.000) 0.815 6.267 ( 0.319 0.553 0.000) 0.639 6.925 ( 0.273 0.472 0.000) 0.545 7.035 ( 1.973 3.418 0.000) 3.947 7.050 ( -0.455 -0.788 0.000) 0.910 7.205 ( -2.019 -3.497 0.000) 4.038 8.860 ( -2.239 -3.877 0.000) 4.477 8.867 ( -1.977 -3.424 0.000) 3.954 ======================= Grid point 34 (17/54) ======================= q-point: ( 0.20 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.417 ( 8.245 9.423 0.000) 12.521 1.479 ( 10.804 10.014 0.000) 14.731 1.992 ( -2.233 11.335 0.000) 11.553 2.020 ( 0.528 17.509 0.000) 17.517 2.211 ( 5.561 4.182 0.000) 6.958 2.327 ( 3.515 -3.186 0.000) 4.744 2.706 ( 2.450 17.328 0.000) 17.500 2.712 ( 2.667 17.169 0.000) 17.375 3.240 ( 8.693 4.862 0.000) 9.961 3.332 ( 8.038 4.420 0.000) 9.173 3.939 ( 7.914 7.023 0.000) 10.581 3.991 ( 9.666 7.477 0.000) 12.220 5.812 ( 0.937 0.611 0.000) 1.119 5.838 ( 1.188 0.465 0.000) 1.276 6.062 ( -1.521 1.514 0.000) 2.146 6.091 ( -1.335 1.660 0.000) 2.131 6.253 ( -0.673 2.605 0.000) 2.691 6.279 ( -0.661 2.459 0.000) 2.547 6.927 ( 1.154 -2.912 0.000) 3.133 7.030 ( 0.671 -3.776 0.000) 3.835 7.052 ( -0.585 -2.294 0.000) 2.368 7.121 ( -2.570 -4.791 0.000) 5.437 8.793 ( -1.188 -2.789 0.000) 3.031 8.811 ( -0.784 -2.506 0.000) 2.626 ======================= Grid point 35 (18/54) ======================= q-point: ( 0.27 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.661 ( 9.039 7.718 0.000) 11.886 1.756 ( 10.330 6.742 0.000) 12.335 1.957 ( -3.771 -1.559 0.000) 4.080 2.088 ( -5.457 4.725 0.000) 7.218 2.342 ( -2.045 11.178 0.000) 11.363 2.368 ( -0.989 7.982 0.000) 8.043 2.871 ( -3.143 16.531 0.000) 16.827 2.880 ( -3.191 16.915 0.000) 17.213 3.416 ( 6.662 2.539 0.000) 7.130 3.490 ( 6.024 2.097 0.000) 6.379 4.182 ( 11.248 5.424 0.000) 12.487 4.256 ( 12.387 4.366 0.000) 13.134 5.879 ( 4.356 1.716 0.000) 4.682 5.906 ( 4.304 1.177 0.000) 4.462 6.072 ( -0.041 4.324 0.000) 4.324 6.104 ( -0.073 4.377 0.000) 4.378 6.270 ( -1.185 4.185 0.000) 4.350 6.296 ( -1.125 4.082 0.000) 4.234 6.899 ( 0.812 -6.151 0.000) 6.204 6.978 ( -2.828 -5.416 0.000) 6.110 6.989 ( 0.650 -6.725 0.000) 6.756 7.016 ( -3.164 -5.544 0.000) 6.383 8.762 ( 0.434 -2.136 0.000) 2.179 8.786 ( 0.617 -2.163 0.000) 2.250 ======================= Grid point 36 (19/54) ======================= q-point: ( 0.33 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.876 ( -2.068 -3.025 0.000) 3.664 1.882 ( 7.183 6.158 0.000) 9.461 1.970 ( 7.009 4.236 0.000) 8.190 1.997 ( -4.341 -2.611 0.000) 5.066 2.380 ( -7.572 14.031 0.000) 15.944 2.409 ( -7.591 14.120 0.000) 16.031 2.946 ( -6.448 15.262 0.000) 16.569 2.958 ( -6.621 15.788 0.000) 17.121 3.530 ( 3.594 2.257 0.000) 4.244 3.592 ( 3.230 1.869 0.000) 3.732 4.420 ( 10.513 2.025 0.000) 10.707 4.492 ( 10.628 0.395 0.000) 10.635 5.959 ( 2.895 -0.251 0.000) 2.906 5.978 ( 3.215 -0.358 0.000) 3.235 6.162 ( 3.685 5.797 0.000) 6.869 6.194 ( 3.810 5.948 0.000) 7.063 6.299 ( -1.205 5.556 0.000) 5.685 6.325 ( -1.181 5.399 0.000) 5.526 6.833 ( 0.348 -8.132 0.000) 8.140 6.861 ( -2.931 -6.115 0.000) 6.781 6.900 ( -3.053 -5.827 0.000) 6.579 6.915 ( 0.051 -8.762 0.000) 8.762 8.760 ( 1.574 -2.051 0.000) 2.585 8.782 ( 1.316 -2.388 0.000) 2.727 ======================= Grid point 37 (20/54) ======================= q-point: ( 0.40 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.832 ( -0.041 -1.901 0.000) 1.902 1.916 ( -1.370 -2.278 0.000) 2.658 2.023 ( 2.444 3.521 0.000) 4.286 2.082 ( 1.919 1.916 0.000) 2.712 2.377 ( -9.315 15.551 0.000) 18.127 2.404 ( -9.435 15.481 0.000) 18.130 2.967 ( -7.954 14.484 0.000) 16.524 2.981 ( -8.228 14.994 0.000) 17.103 3.591 ( 0.955 1.447 0.000) 1.734 3.644 ( 0.916 0.947 0.000) 1.317 4.573 ( 6.143 -2.320 0.000) 6.566 4.623 ( 5.593 -3.360 0.000) 6.525 5.985 ( 0.823 -0.343 0.000) 0.891 6.013 ( 1.456 -0.848 0.000) 1.684 6.283 ( 2.982 3.788 0.000) 4.821 6.308 ( 2.021 3.099 0.000) 3.700 6.333 ( -2.283 6.416 0.000) 6.810 6.356 ( -2.259 5.926 0.000) 6.343 6.766 ( 2.332 -8.066 0.000) 8.397 6.774 ( 0.124 -4.926 0.000) 4.928 6.806 ( -0.417 -5.227 0.000) 5.244 6.828 ( 1.544 -9.139 0.000) 9.269 8.770 ( 1.869 -2.292 0.000) 2.957 8.780 ( 1.358 -2.738 0.000) 3.056 ======================= Grid point 49 (21/54) ======================= q-point: ( 0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.620 ( 6.012 10.414 0.000) 12.025 1.692 ( 6.141 10.637 0.000) 12.282 1.982 ( -2.966 -5.137 0.000) 5.932 2.163 ( -3.178 -5.504 0.000) 6.355 2.411 ( 7.037 12.188 0.000) 14.074 2.447 ( 7.187 12.448 0.000) 14.374 3.013 ( 6.810 11.795 0.000) 13.620 3.026 ( 7.129 12.348 0.000) 14.259 3.345 ( 3.367 5.832 0.000) 6.734 3.425 ( 2.939 5.090 0.000) 5.878 4.104 ( 5.108 8.847 0.000) 10.216 4.156 ( 4.946 8.567 0.000) 9.893 5.848 ( 1.406 2.435 0.000) 2.812 5.870 ( 1.302 2.254 0.000) 2.603 6.100 ( 1.406 2.436 0.000) 2.813 6.132 ( 1.502 2.601 0.000) 3.003 6.305 ( 1.411 2.444 0.000) 2.822 6.331 ( 1.421 2.462 0.000) 2.842 6.854 ( -2.373 -4.111 0.000) 4.747 6.942 ( -2.675 -4.633 0.000) 5.350 6.973 ( -2.946 -5.102 0.000) 5.891 7.015 ( -3.157 -5.468 0.000) 6.314 8.746 ( -1.023 -1.772 0.000) 2.046 8.768 ( -0.927 -1.606 0.000) 1.854 ======================= Grid point 50 (22/54) ======================= q-point: ( 0.27 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.837 ( 6.580 9.261 0.000) 11.361 1.879 ( -2.488 -4.215 0.000) 4.895 1.905 ( 6.700 7.745 0.000) 10.241 2.026 ( -4.161 -6.034 0.000) 7.329 2.605 ( -0.238 12.671 0.000) 12.673 2.636 ( -0.534 12.441 0.000) 12.452 3.181 ( -1.757 12.991 0.000) 13.110 3.204 ( -1.752 13.509 0.000) 13.622 3.475 ( 5.227 3.270 0.000) 6.166 3.539 ( 4.748 2.637 0.000) 5.431 4.291 ( 7.236 5.275 0.000) 8.954 4.337 ( 7.763 3.785 0.000) 8.636 5.900 ( 2.599 0.449 0.000) 2.637 5.919 ( 2.467 0.242 0.000) 2.478 6.169 ( 2.186 4.797 0.000) 5.271 6.204 ( 2.176 4.992 0.000) 5.445 6.355 ( 0.616 3.894 0.000) 3.943 6.381 ( 0.607 4.015 0.000) 4.061 6.762 ( -1.576 -6.825 0.000) 7.004 6.842 ( -1.778 -7.272 0.000) 7.486 6.855 ( -3.479 -6.283 0.000) 7.182 6.898 ( -3.398 -5.850 0.000) 6.765 8.721 ( 0.266 -1.732 0.000) 1.752 8.744 ( 0.156 -1.901 0.000) 1.907 ======================= Grid point 51 (23/54) ======================= q-point: ( 0.33 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.817 ( -0.707 -2.187 0.000) 2.298 1.912 ( -2.271 -4.290 0.000) 4.854 2.021 ( 3.895 7.125 0.000) 8.121 2.069 ( 3.657 5.288 0.000) 6.429 2.676 ( -5.885 13.049 0.000) 14.314 2.701 ( -5.985 12.787 0.000) 14.118 3.241 ( -5.824 12.862 0.000) 14.119 3.266 ( -6.059 13.213 0.000) 14.536 3.576 ( 2.829 2.090 0.000) 3.517 3.625 ( 2.463 1.211 0.000) 2.744 4.448 ( 6.222 0.890 0.000) 6.285 4.488 ( 6.598 -0.330 0.000) 6.607 5.944 ( 2.520 -1.128 0.000) 2.761 5.960 ( 2.646 -1.302 0.000) 2.949 6.269 ( 3.233 4.514 0.000) 5.552 6.303 ( 2.847 4.539 0.000) 5.358 6.410 ( 0.444 4.991 0.000) 5.011 6.436 ( 0.090 5.079 0.000) 5.080 6.664 ( -0.519 -7.797 0.000) 7.814 6.731 ( -1.259 -8.557 0.000) 8.649 6.736 ( -2.292 -5.758 0.000) 6.197 6.779 ( -2.933 -5.701 0.000) 6.411 8.717 ( 1.257 -2.039 0.000) 2.395 8.731 ( 0.835 -2.362 0.000) 2.505 ======================= Grid point 52 (24/54) ======================= q-point: ( 0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.803 ( 0.521 -0.902 0.000) 1.042 1.867 ( 1.143 -1.980 0.000) 2.286 2.097 ( -2.022 3.503 0.000) 4.044 2.130 ( -1.465 2.538 0.000) 2.930 2.687 ( -7.625 13.208 0.000) 15.251 2.710 ( -7.542 13.063 0.000) 15.083 3.253 ( -7.246 12.550 0.000) 14.492 3.277 ( -7.423 12.857 0.000) 14.846 3.614 ( -0.382 0.661 0.000) 0.764 3.652 ( 0.114 -0.197 0.000) 0.228 4.511 ( 1.823 -3.157 0.000) 3.645 4.547 ( 2.101 -3.639 0.000) 4.202 5.961 ( 1.120 -1.940 0.000) 2.240 5.978 ( 1.389 -2.406 0.000) 2.778 6.342 ( -0.897 1.553 0.000) 1.794 6.356 ( -0.926 1.603 0.000) 1.851 6.453 ( -2.591 4.488 0.000) 5.182 6.466 ( -2.469 4.277 0.000) 4.938 6.615 ( 3.450 -5.976 0.000) 6.901 6.657 ( 3.799 -6.580 0.000) 7.598 6.683 ( 2.030 -3.516 0.000) 4.060 6.708 ( 2.099 -3.635 0.000) 4.198 8.720 ( 1.393 -2.413 0.000) 2.787 8.725 ( 1.338 -2.318 0.000) 2.677 ======================= Grid point 66 (25/54) ======================= q-point: ( 0.27 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.812 ( -1.199 -2.077 0.000) 2.398 1.910 ( -2.721 -4.713 0.000) 5.442 2.025 ( 4.831 8.367 0.000) 9.661 2.063 ( 4.112 7.123 0.000) 8.225 2.803 ( 3.724 6.450 0.000) 7.448 2.828 ( 3.569 6.181 0.000) 7.138 3.382 ( 3.646 6.314 0.000) 7.291 3.407 ( 3.527 6.109 0.000) 7.054 3.548 ( 2.132 3.692 0.000) 4.264 3.596 ( 1.599 2.769 0.000) 3.198 4.395 ( 2.628 4.552 0.000) 5.256 4.418 ( 2.365 4.097 0.000) 4.730 5.908 ( 0.317 0.550 0.000) 0.635 5.923 ( 0.195 0.338 0.000) 0.390 6.263 ( 2.510 4.348 0.000) 5.020 6.301 ( 2.536 4.392 0.000) 5.072 6.429 ( 1.975 3.420 0.000) 3.950 6.458 ( 2.034 3.523 0.000) 4.068 6.629 ( -3.420 -5.923 0.000) 6.840 6.699 ( -3.770 -6.529 0.000) 7.539 6.730 ( -3.229 -5.593 0.000) 6.459 6.780 ( -3.223 -5.583 0.000) 6.446 8.693 ( -0.557 -0.965 0.000) 1.114 8.710 ( -0.762 -1.320 0.000) 1.524 ======================= Grid point 67 (26/54) ======================= q-point: ( 0.33 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 226 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.792 ( 0.050 -0.300 0.000) 0.304 1.840 ( -0.773 -2.640 0.000) 2.751 2.152 ( 1.107 4.981 0.000) 5.102 2.173 ( 0.995 4.292 0.000) 4.406 2.881 ( -2.001 6.612 0.000) 6.908 2.902 ( -2.048 6.559 0.000) 6.872 3.453 ( -2.781 7.412 0.000) 7.917 3.475 ( -2.505 6.837 0.000) 7.282 3.612 ( 1.633 1.321 0.000) 2.100 3.637 ( 1.082 -0.081 0.000) 1.085 4.463 ( 1.747 0.655 0.000) 1.865 4.488 ( 2.210 0.268 0.000) 2.226 5.921 ( 1.089 -0.847 0.000) 1.379 5.931 ( 1.074 -1.223 0.000) 1.628 6.346 ( 2.240 3.026 0.000) 3.765 6.376 ( 1.057 2.396 0.000) 2.619 6.496 ( 1.503 3.293 0.000) 3.620 6.522 ( 0.596 3.067 0.000) 3.124 6.533 ( -0.392 -4.548 0.000) 4.565 6.580 ( -1.893 -5.893 0.000) 6.190 6.638 ( -1.075 -3.499 0.000) 3.661 6.676 ( -2.344 -4.178 0.000) 4.790 8.682 ( 0.484 -1.215 0.000) 1.308 8.691 ( 0.197 -1.495 0.000) 1.508 ======================= Grid point 82 (27/54) ======================= q-point: ( 0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 54 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.799 ( -0.000 -0.000 0.000) 0.000 1.805 ( -0.000 -0.000 0.000) 0.000 2.205 ( -0.000 -0.000 0.000) 0.000 2.222 ( -0.000 -0.000 0.000) 0.000 2.948 ( -0.000 -0.000 0.000) 0.000 2.969 ( -0.000 -0.000 0.000) 0.000 3.535 ( -0.000 -0.000 0.000) 0.000 3.543 ( 0.000 0.000 0.000) 0.000 3.628 ( -0.000 -0.000 0.000) 0.000 3.630 ( 0.000 0.000 0.000) 0.000 4.471 ( -0.000 -0.000 0.000) 0.000 4.492 ( -0.000 -0.000 0.000) 0.000 5.914 ( 0.000 0.000 0.000) 0.000 5.916 ( -0.000 -0.000 0.000) 0.000 6.390 ( -0.000 -0.000 0.000) 0.000 6.393 ( 0.000 0.000 0.000) 0.000 6.495 ( 0.000 0.000 0.000) 0.000 6.503 ( -0.000 -0.000 0.000) 0.000 6.539 ( -0.000 -0.000 0.000) 0.000 6.551 ( -0.000 -0.000 0.000) 0.000 6.607 ( 0.000 0.000 0.000) 0.000 6.620 ( -0.000 -0.000 0.000) 0.000 8.669 ( -0.000 -0.000 0.000) 0.000 8.674 ( -0.000 -0.000 0.000) 0.000 ======================= Grid point 241 (28/54) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 54 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.241 ( 0.000 0.000 10.710) 10.710 0.241 ( 0.000 0.000 10.710) 10.710 0.376 ( 0.000 -0.000 16.676) 16.676 0.424 ( 0.000 0.000 -6.447) 6.447 0.424 ( 0.000 0.000 -6.447) 6.447 0.659 ( 0.000 -0.000 -9.978) 9.978 1.617 ( 0.000 0.000 0.539) 0.539 1.617 ( 0.000 0.000 0.539) 0.539 1.637 ( 0.000 0.000 -0.443) 0.443 1.637 ( 0.000 0.000 -0.443) 0.443 3.790 ( 0.000 0.000 1.921) 1.921 3.833 ( 0.000 -0.000 -1.850) 1.850 6.096 ( 0.000 0.000 0.115) 0.115 6.096 ( 0.000 0.000 0.115) 0.115 6.115 ( 0.000 0.000 -0.193) 0.193 6.115 ( 0.000 0.000 -0.193) 0.193 6.273 ( -0.000 0.000 0.567) 0.567 6.273 ( 0.000 0.000 0.567) 0.567 6.294 ( 0.000 0.000 -0.489) 0.489 6.294 ( 0.000 0.000 -0.489) 0.489 7.201 ( -0.000 0.000 1.736) 1.736 7.245 ( 0.000 0.000 -1.724) 1.724 9.124 ( 0.000 0.000 0.270) 0.270 9.131 ( 0.000 0.000 -0.257) 0.257 ======================= Grid point 242 (29/54) ======================= q-point: ( 0.07 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.409 ( 7.769 4.485 10.618) 13.900 0.527 ( 17.235 9.950 5.190) 20.566 0.635 ( 3.152 1.820 -9.202) 9.895 0.637 ( 14.649 8.457 -4.539) 17.513 0.915 ( 33.588 19.392 2.662) 38.876 0.977 ( 30.091 17.373 -1.360) 34.773 1.696 ( 6.370 3.678 0.495) 7.372 1.714 ( 6.203 3.581 -0.396) 7.174 1.899 ( 19.863 11.468 0.526) 22.942 1.912 ( 20.576 11.879 0.043) 23.759 3.758 ( -2.344 -1.353 1.985) 3.356 3.803 ( -2.174 -1.255 -2.043) 3.237 6.086 ( -0.858 -0.496 0.145) 1.002 6.105 ( -0.860 -0.497 -0.237) 1.021 6.125 ( 1.281 0.739 0.122) 1.484 6.143 ( 1.164 0.672 0.161) 1.354 6.264 ( -0.737 -0.426 0.602) 1.042 6.285 ( -0.747 -0.431 -0.511) 1.002 6.493 ( 13.935 8.045 4.906) 16.822 6.721 ( 17.351 10.018 -18.909) 27.549 7.162 ( -1.999 -1.154 -0.770) 2.433 7.183 ( -4.368 -2.522 -1.325) 5.215 9.091 ( -2.690 -1.553 0.236) 3.115 9.097 ( -2.743 -1.584 -0.231) 3.176 ======================= Grid point 243 (30/54) ======================= q-point: ( 0.13 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.685 ( 14.313 8.264 4.274) 17.071 0.794 ( 10.154 5.862 -5.123) 12.795 0.941 ( 17.001 9.816 3.109) 19.876 1.012 ( 15.895 9.177 -3.049) 18.605 1.659 ( 27.353 15.792 0.524) 31.589 1.678 ( 27.332 15.780 -0.093) 31.560 1.895 ( 9.833 5.677 0.446) 11.363 1.908 ( 9.665 5.580 -0.358) 11.166 2.435 ( 23.638 13.647 1.992) 27.367 2.473 ( 24.482 14.135 -1.295) 28.299 3.706 ( -1.226 -0.708 1.782) 2.276 3.749 ( -1.838 -1.061 -1.928) 2.867 6.058 ( -1.437 -0.830 0.224) 1.674 6.077 ( -1.430 -0.826 -0.334) 1.685 6.100 ( -4.286 -2.474 0.181) 4.952 6.114 ( -4.526 -2.613 -0.046) 5.226 6.240 ( -1.277 -0.737 0.655) 1.614 6.260 ( -1.309 -0.756 -0.546) 1.607 6.805 ( 10.751 6.207 4.618) 13.246 6.907 ( 2.568 1.483 -2.977) 4.202 7.099 ( -1.427 -0.824 0.478) 1.716 7.172 ( 2.155 1.244 -7.184) 7.603 9.007 ( -4.093 -2.363 0.143) 4.729 9.011 ( -4.185 -2.416 -0.150) 4.835 ======================= Grid point 244 (31/54) ======================= q-point: ( 0.20 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.043 ( 15.181 8.765 1.428) 17.588 1.086 ( 14.001 8.084 -2.299) 16.330 1.313 ( 14.001 8.083 2.103) 16.303 1.363 ( 13.295 7.676 -2.125) 15.499 2.135 ( 9.880 5.704 0.533) 11.421 2.146 ( 9.923 5.729 -0.467) 11.467 2.154 ( 13.735 7.930 2.379) 16.037 2.205 ( 15.558 8.982 -2.084) 18.085 2.953 ( 18.985 10.961 2.438) 22.057 3.003 ( 19.124 11.041 -1.902) 22.164 3.744 ( 5.115 2.953 0.617) 5.939 3.761 ( 4.176 2.411 -0.840) 4.895 5.942 ( -7.438 -4.294 0.164) 8.590 5.952 ( -7.335 -4.235 -0.115) 8.471 6.022 ( -1.494 -0.863 0.303) 1.752 6.041 ( -1.485 -0.858 -0.408) 1.763 6.207 ( -1.365 -0.788 0.655) 1.707 6.226 ( -1.421 -0.821 -0.550) 1.731 6.961 ( 3.567 2.060 0.474) 4.146 7.003 ( 4.904 2.831 0.053) 5.663 7.119 ( 1.946 1.123 2.446) 3.321 7.204 ( 0.148 0.085 -3.967) 3.971 8.911 ( -3.714 -2.145 0.161) 4.292 8.915 ( -3.581 -2.068 -0.172) 4.139 ======================= Grid point 245 (32/54) ======================= q-point: ( 0.27 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.385 ( 13.582 7.842 1.823) 15.789 1.434 ( 14.493 8.368 -2.409) 16.908 1.592 ( 9.396 5.425 1.466) 10.948 1.630 ( 9.114 5.262 -1.509) 10.632 2.309 ( 0.320 0.185 3.708) 3.726 2.345 ( 7.594 4.385 0.612) 8.791 2.360 ( 7.755 4.477 -0.555) 8.972 2.392 ( 0.921 0.532 -3.564) 3.719 3.302 ( 10.271 5.930 2.220) 12.066 3.348 ( 9.937 5.737 -1.868) 11.625 3.945 ( 11.704 6.757 1.381) 13.585 3.980 ( 12.915 7.456 -1.466) 14.984 5.831 ( -0.741 -0.428 0.710) 1.112 5.848 ( -0.444 -0.256 -0.672) 0.846 5.991 ( -0.991 -0.572 0.347) 1.196 6.011 ( -0.950 -0.548 -0.443) 1.183 6.180 ( -0.874 -0.505 0.623) 1.186 6.197 ( -0.915 -0.528 -0.526) 1.181 7.029 ( 1.760 1.016 0.842) 2.199 7.079 ( 1.481 0.855 -0.919) 1.941 7.138 ( -0.802 -0.463 0.522) 1.063 7.178 ( -2.182 -1.260 -0.604) 2.591 8.842 ( -1.960 -1.132 0.410) 2.300 8.852 ( -1.704 -0.984 -0.415) 2.011 ======================= Grid point 246 (33/54) ======================= q-point: ( 0.33 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.681 ( 11.165 6.446 2.596) 13.151 1.746 ( 11.497 6.638 -3.026) 13.616 1.756 ( 4.505 2.601 1.270) 5.355 1.792 ( 4.575 2.642 -1.320) 5.446 2.234 ( -5.376 -3.104 3.655) 7.204 2.316 ( -5.880 -3.395 -3.524) 7.650 2.488 ( 4.462 2.576 0.606) 5.187 2.506 ( 4.553 2.628 -0.548) 5.285 3.464 ( 4.332 2.501 1.675) 5.275 3.502 ( 4.006 2.313 -1.479) 4.857 4.253 ( 13.207 7.625 2.159) 15.403 4.307 ( 13.598 7.851 -2.145) 15.847 5.911 ( 6.787 3.918 0.961) 7.895 5.934 ( 6.925 3.998 -0.922) 8.049 5.977 ( -0.238 -0.137 0.422) 0.503 5.999 ( -0.130 -0.075 -0.529) 0.550 6.168 ( -0.069 -0.040 0.649) 0.654 6.186 ( -0.056 -0.032 -0.541) 0.545 7.028 ( -2.084 -1.203 -0.131) 2.410 7.058 ( -3.398 -1.962 0.389) 3.943 7.090 ( -2.603 -1.503 0.515) 3.050 7.124 ( -2.008 -1.159 -0.775) 2.445 8.818 ( -0.241 -0.139 0.607) 0.668 8.832 ( -0.134 -0.077 -0.607) 0.626 ======================= Grid point 247 (34/54) ======================= q-point: ( 0.40 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.819 ( 1.153 0.666 1.234) 1.815 1.856 ( 1.137 0.656 -1.275) 1.830 1.904 ( 7.434 4.292 2.363) 8.904 1.965 ( 6.732 3.887 -2.856) 8.282 2.100 ( -5.083 -2.935 3.000) 6.592 2.165 ( -6.001 -3.465 -2.649) 7.418 2.561 ( 1.882 1.086 0.512) 2.232 2.580 ( 1.851 1.069 -0.461) 2.187 3.542 ( 2.610 1.507 1.332) 3.295 3.575 ( 2.462 1.421 -1.245) 3.103 4.535 ( 10.017 5.783 1.999) 11.738 4.589 ( 9.600 5.542 -1.998) 11.263 5.978 ( 0.250 0.144 0.540) 0.612 6.002 ( 0.355 0.205 -0.678) 0.792 6.103 ( 8.307 4.796 1.003) 9.645 6.126 ( 8.247 4.761 -0.975) 9.573 6.174 ( 0.441 0.255 0.751) 0.907 6.192 ( 0.484 0.279 -0.613) 0.829 6.944 ( -4.079 -2.355 0.201) 4.714 6.960 ( -4.142 -2.391 -0.695) 4.833 7.026 ( -2.822 -1.629 2.117) 3.886 7.067 ( -2.867 -1.655 -1.586) 3.671 8.823 ( 0.516 0.298 0.550) 0.810 8.837 ( 0.360 0.208 -0.551) 0.690 ======================= Grid point 248 (35/54) ======================= q-point: ( 0.47 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.829 ( 0.040 0.023 1.071) 1.072 1.864 ( -0.048 -0.028 -1.102) 1.104 1.993 ( -1.151 -0.664 1.271) 1.839 2.032 ( 0.656 0.379 0.188) 0.781 2.052 ( -0.236 -0.136 -0.563) 0.626 2.091 ( 1.082 0.625 -0.917) 1.550 2.587 ( 0.445 0.257 0.376) 0.637 2.604 ( 0.398 0.230 -0.337) 0.570 3.588 ( 1.099 0.634 1.206) 1.751 3.619 ( 1.055 0.609 -1.188) 1.702 4.706 ( 3.973 2.294 1.004) 4.696 4.744 ( 3.286 1.897 -1.038) 3.933 5.984 ( 0.179 0.103 0.619) 0.653 6.010 ( 0.217 0.125 -0.784) 0.823 6.184 ( 0.267 0.154 0.840) 0.894 6.203 ( 0.289 0.167 -0.675) 0.753 6.260 ( 4.064 2.346 0.653) 4.738 6.277 ( 3.627 2.094 -0.623) 4.234 6.870 ( -1.763 -1.018 0.519) 2.101 6.885 ( -1.796 -1.037 -0.850) 2.242 6.968 ( -1.628 -0.940 1.723) 2.550 7.003 ( -1.914 -1.105 -1.361) 2.596 8.835 ( 0.349 0.201 0.252) 0.475 8.843 ( 0.131 0.076 -0.256) 0.298 ======================= Grid point 257 (36/54) ======================= q-point: ( 0.07 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.589 ( 7.414 12.842 5.393) 15.779 0.726 ( 4.683 8.111 -6.327) 11.303 0.897 ( 11.778 20.400 2.352) 23.673 0.953 ( 10.970 19.000 -2.504) 22.082 1.443 ( 16.496 28.572 0.946) 33.006 1.469 ( 15.487 26.825 -0.152) 30.975 1.855 ( 6.369 11.032 0.399) 12.745 1.869 ( 6.225 10.781 -0.263) 12.452 2.269 ( 12.921 22.380 1.712) 25.898 2.301 ( 13.517 23.413 -1.015) 27.054 3.716 ( -1.213 -2.101 1.929) 3.100 3.761 ( -1.386 -2.401 -2.054) 3.450 6.057 ( -1.441 -2.495 0.195) 2.888 6.076 ( -1.501 -2.600 -0.204) 3.010 6.134 ( -0.403 -0.699 0.133) 0.818 6.149 ( -0.518 -0.897 -0.085) 1.040 6.253 ( -0.400 -0.693 0.642) 1.026 6.272 ( -0.445 -0.771 -0.503) 1.023 6.728 ( 6.994 12.115 5.820) 15.151 6.884 ( 2.001 3.466 -6.576) 7.698 7.110 ( -1.891 -3.275 -0.176) 3.785 7.156 ( 0.816 1.414 -6.559) 6.759 9.032 ( -2.272 -3.935 0.163) 4.547 9.037 ( -2.338 -4.049 -0.162) 4.679 ======================= Grid point 258 (37/54) ======================= q-point: ( 0.13 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.899 ( 11.931 12.187 1.643) 17.134 0.952 ( 10.041 9.461 -2.953) 14.109 1.294 ( 8.189 21.363 1.847) 22.953 1.337 ( 8.514 20.065 -1.761) 21.868 1.966 ( 16.162 12.637 1.361) 20.561 1.995 ( 17.918 13.800 -0.999) 22.639 2.110 ( 6.677 14.416 0.325) 15.890 2.120 ( 6.678 14.143 -0.151) 15.641 2.745 ( 17.139 15.358 2.332) 23.131 2.793 ( 17.431 15.686 -1.787) 23.517 3.706 ( 1.711 1.291 1.109) 2.414 3.734 ( 0.713 0.191 -1.313) 1.507 5.970 ( -3.504 -6.246 0.175) 7.164 5.983 ( -3.738 -6.370 -0.097) 7.386 6.089 ( -3.498 0.065 0.034) 3.499 6.102 ( -3.282 0.118 -0.185) 3.290 6.235 ( -1.308 0.275 0.659) 1.490 6.253 ( -1.286 0.138 -0.510) 1.390 6.914 ( 3.496 2.497 0.204) 4.301 6.944 ( 4.156 3.156 1.077) 5.329 7.090 ( 1.420 0.166 2.468) 2.853 7.190 ( 0.993 0.160 -5.583) 5.673 8.942 ( -3.082 -3.777 0.055) 4.876 8.944 ( -3.053 -3.888 -0.060) 4.944 ======================= Grid point 259 (38/54) ======================= q-point: ( 0.20 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.232 ( 12.468 9.354 0.860) 15.610 1.262 ( 13.674 9.006 -1.751) 16.467 1.613 ( 3.903 19.184 1.391) 19.626 1.647 ( 3.906 18.228 -1.235) 18.683 2.244 ( 7.098 0.985 3.096) 7.806 2.313 ( 8.137 2.143 -2.966) 8.922 2.360 ( 3.799 15.034 0.406) 15.512 2.369 ( 4.086 14.630 -0.289) 15.193 3.147 ( 13.928 7.449 2.368) 15.971 3.196 ( 13.770 7.335 -1.968) 15.725 3.826 ( 7.880 4.990 0.597) 9.346 3.843 ( 9.209 5.562 -0.756) 10.785 5.846 ( -2.720 -3.459 0.468) 4.425 5.860 ( -2.541 -3.322 -0.423) 4.204 6.036 ( -3.080 1.351 0.280) 3.375 6.054 ( -2.895 1.306 -0.399) 3.200 6.215 ( -1.877 1.502 0.615) 2.481 6.232 ( -1.857 1.362 -0.504) 2.357 6.979 ( 3.306 -0.473 1.074) 3.508 7.033 ( 3.349 -0.238 -1.586) 3.714 7.122 ( 1.247 -0.927 2.154) 2.656 7.182 ( -0.690 -2.303 -1.991) 3.121 8.858 ( -2.162 -2.930 0.256) 3.650 8.863 ( -1.895 -2.822 -0.262) 3.410 ======================= Grid point 260 (39/54) ======================= q-point: ( 0.27 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.536 ( 12.108 6.935 2.002) 14.097 1.589 ( 13.137 6.898 -2.568) 15.059 1.816 ( -0.583 14.399 0.968) 14.443 1.843 ( -0.244 13.508 -0.676) 13.527 2.273 ( -1.526 -2.276 3.284) 4.277 2.350 ( -1.467 -3.482 -3.494) 5.146 2.544 ( -0.006 14.052 0.431) 14.059 2.556 ( 0.352 13.751 -0.282) 13.759 3.386 ( 7.669 2.729 1.925) 8.365 3.429 ( 7.313 2.515 -1.682) 7.914 4.082 ( 13.177 6.299 1.725) 14.706 4.126 ( 14.082 6.419 -1.739) 15.573 5.846 ( 3.613 1.825 0.778) 4.122 5.864 ( 3.820 1.734 -0.727) 4.257 6.009 ( -1.855 2.157 0.378) 2.870 6.029 ( -1.739 2.144 -0.489) 2.804 6.203 ( -1.799 2.663 0.590) 3.267 6.219 ( -1.755 2.552 -0.496) 3.136 7.008 ( 1.574 -3.059 0.770) 3.525 7.055 ( 0.846 -3.220 -0.843) 3.434 7.095 ( -1.400 -3.656 -0.125) 3.917 7.124 ( -1.770 -3.986 0.154) 4.364 8.811 ( -0.324 -1.895 0.506) 1.988 8.822 ( -0.097 -1.823 -0.507) 1.895 ======================= Grid point 261 (40/54) ======================= q-point: ( 0.33 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.793 ( 9.656 4.875 2.293) 11.057 1.851 ( 8.571 4.402 -2.556) 9.968 1.895 ( -1.658 7.123 0.900) 7.368 1.920 ( -2.845 7.365 -0.524) 7.913 2.196 ( -5.964 2.038 2.094) 6.642 2.249 ( -5.929 -1.013 -2.495) 6.512 2.650 ( -3.153 12.324 0.309) 12.725 2.664 ( -2.897 12.018 -0.200) 12.363 3.505 ( 3.810 1.812 1.466) 4.467 3.540 ( 3.568 1.693 -1.347) 4.173 4.378 ( 12.879 4.364 1.993) 13.743 4.430 ( 13.086 3.974 -1.979) 13.818 5.960 ( 3.913 0.622 0.445) 3.987 5.973 ( 3.977 0.541 -0.189) 4.019 6.044 ( 3.292 5.603 0.233) 6.503 6.062 ( 3.254 5.753 -0.572) 6.635 6.207 ( -1.416 3.678 0.629) 3.991 6.223 ( -1.332 3.571 -0.518) 3.846 6.964 ( -1.857 -4.774 -0.231) 5.128 6.986 ( -2.964 -4.872 0.172) 5.705 7.022 ( -1.377 -5.067 1.132) 5.371 7.056 ( -0.586 -5.496 -1.053) 5.626 8.802 ( 1.030 -1.347 0.583) 1.793 8.816 ( 1.036 -1.427 -0.582) 1.857 ======================= Grid point 262 (41/54) ======================= q-point: ( 0.40 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.877 ( -1.826 0.930 0.897) 2.237 1.903 ( -2.802 2.292 -0.640) 3.676 1.979 ( 5.211 3.058 1.901) 6.334 2.026 ( 4.728 2.450 -1.925) 5.663 2.126 ( -7.027 7.459 0.720) 10.273 2.151 ( -6.976 5.053 -1.070) 8.681 2.694 ( -5.101 10.984 0.181) 12.112 2.708 ( -4.981 10.734 -0.114) 11.834 3.573 ( 1.936 1.525 1.236) 2.757 3.604 ( 1.835 1.438 -1.203) 2.623 4.604 ( 8.862 0.910 1.412) 9.020 4.647 ( 8.383 0.261 -1.426) 8.508 5.987 ( 0.274 0.602 0.560) 0.867 6.009 ( 0.666 0.354 -0.615) 0.973 6.201 ( 5.657 4.473 0.591) 7.236 6.218 ( 2.841 4.301 -0.176) 5.157 6.228 ( 1.397 4.549 0.097) 4.760 6.241 ( -1.170 4.374 -0.436) 4.549 6.874 ( -1.899 -4.338 0.306) 4.745 6.889 ( -1.895 -4.484 -0.734) 4.923 6.947 ( -0.477 -5.944 1.845) 6.242 6.982 ( -0.428 -6.461 -1.381) 6.621 8.812 ( 1.415 -1.348 0.381) 1.991 8.822 ( 1.213 -1.547 -0.384) 2.003 ======================= Grid point 263 (42/54) ======================= q-point: ( 0.47 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.862 ( -0.007 0.013 0.514) 0.514 1.883 ( -1.020 1.766 -0.259) 2.056 2.046 ( -0.172 0.298 1.281) 1.326 2.077 ( -1.363 2.361 -1.393) 3.062 2.105 ( -4.692 8.127 0.128) 9.385 2.113 ( -4.732 8.196 -0.311) 9.470 2.703 ( -5.962 10.326 0.100) 11.924 2.717 ( -5.855 10.141 -0.051) 11.710 3.600 ( -0.243 0.420 1.160) 1.258 3.629 ( -0.208 0.361 -1.164) 1.236 4.691 ( 2.036 -3.526 0.643) 4.122 4.723 ( 2.091 -3.621 -0.673) 4.235 5.993 ( -0.308 0.534 0.627) 0.880 6.016 ( -0.155 0.268 -0.716) 0.780 6.233 ( -2.502 4.333 0.690) 5.051 6.249 ( -2.329 4.035 -0.658) 4.705 6.291 ( -0.407 0.706 0.412) 0.913 6.304 ( -0.310 0.537 -0.321) 0.698 6.833 ( 1.441 -2.496 0.455) 2.918 6.846 ( 1.541 -2.669 -0.777) 3.178 6.906 ( 2.730 -4.728 1.489) 5.659 6.937 ( 2.925 -5.066 -1.150) 5.962 8.819 ( 0.973 -1.685 0.106) 1.949 8.825 ( 1.011 -1.751 -0.111) 2.024 ======================= Grid point 274 (43/54) ======================= q-point: ( 0.13 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.169 ( 7.805 13.519 -0.590) 15.622 1.174 ( 7.003 12.129 -0.496) 14.015 1.722 ( 10.851 18.794 1.042) 21.726 1.748 ( 10.570 18.309 -0.845) 21.158 2.150 ( 3.077 5.330 2.836) 6.776 2.212 ( 3.987 6.905 -2.635) 8.398 2.432 ( 8.998 15.586 0.199) 17.998 2.436 ( 8.887 15.393 -0.053) 17.774 3.036 ( 7.327 12.691 2.411) 14.852 3.087 ( 7.305 12.652 -2.023) 14.749 3.781 ( 3.292 5.702 -0.092) 6.585 3.786 ( 3.555 6.158 -0.088) 7.111 5.852 ( -2.476 -4.288 0.380) 4.966 5.865 ( -2.454 -4.250 -0.331) 4.918 6.082 ( -0.418 -0.723 0.205) 0.860 6.097 ( -0.347 -0.602 -0.357) 0.781 6.246 ( 0.380 0.659 0.615) 0.978 6.261 ( 0.328 0.568 -0.467) 0.806 6.939 ( 0.419 0.726 1.196) 1.460 6.993 ( 0.612 1.059 -1.947) 2.300 7.101 ( 0.170 0.295 2.951) 2.970 7.172 ( -1.093 -1.894 -2.667) 3.449 8.861 ( -2.169 -3.757 0.145) 4.341 8.865 ( -2.045 -3.542 -0.148) 4.093 ======================= Grid point 275 (44/54) ======================= q-point: ( 0.20 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.426 ( 8.995 9.686 0.982) 13.255 1.458 ( 10.261 10.007 -1.646) 14.427 2.007 ( -1.376 14.441 0.906) 14.535 2.024 ( -0.309 16.466 0.154) 16.469 2.232 ( 5.009 1.232 2.119) 5.577 2.291 ( 3.908 -2.376 -2.992) 5.465 2.708 ( 2.483 17.210 0.130) 17.389 2.713 ( 2.675 17.275 0.004) 17.481 3.265 ( 8.480 4.676 2.118) 9.913 3.311 ( 8.175 4.490 -1.863) 9.511 3.951 ( 8.382 7.134 0.981) 11.050 3.979 ( 9.259 7.431 -1.031) 11.917 5.818 ( 1.018 0.579 0.545) 1.292 5.833 ( 1.113 0.497 -0.501) 1.318 6.066 ( -1.475 1.568 0.363) 2.183 6.084 ( -1.363 1.609 -0.485) 2.164 6.260 ( -0.689 2.572 0.560) 2.721 6.274 ( -0.674 2.478 -0.455) 2.608 6.943 ( 1.040 -2.921 1.297) 3.361 6.994 ( 0.831 -3.216 -2.000) 3.877 7.079 ( -1.159 -3.292 1.226) 3.699 7.113 ( -1.945 -4.234 -0.614) 4.699 8.798 ( -1.088 -2.719 0.374) 2.952 8.806 ( -0.883 -2.573 -0.374) 2.746 ======================= Grid point 276 (45/54) ======================= q-point: ( 0.27 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.681 ( 9.383 7.530 1.805) 12.165 1.729 ( 10.014 7.082 -2.205) 12.462 1.994 ( -4.186 -0.587 3.124) 5.257 2.061 ( -5.067 2.411 -2.571) 6.172 2.343 ( -1.726 10.740 0.172) 10.879 2.359 ( -1.248 9.751 -0.644) 9.852 2.872 ( -3.232 16.656 0.096) 16.967 2.881 ( -3.119 16.834 0.076) 17.121 3.435 ( 6.471 2.351 1.623) 7.073 3.473 ( 6.171 2.159 -1.510) 6.710 4.199 ( 11.537 5.111 1.426) 12.699 4.239 ( 12.112 4.671 -1.425) 13.060 5.886 ( 4.373 1.578 0.581) 4.685 5.900 ( 4.350 1.320 -0.508) 4.574 6.077 ( -0.106 4.332 0.402) 4.352 6.097 ( -0.072 4.354 -0.542) 4.388 6.277 ( -1.173 4.167 0.542) 4.362 6.291 ( -1.132 4.099 -0.455) 4.277 6.910 ( 0.448 -6.099 0.916) 6.183 6.954 ( -0.043 -6.317 -1.352) 6.460 6.997 ( -2.398 -5.589 0.211) 6.086 7.020 ( -2.507 -5.814 0.223) 6.335 8.768 ( 0.476 -2.147 0.521) 2.259 8.780 ( 0.574 -2.152 -0.519) 2.287 ======================= Grid point 277 (46/54) ======================= q-point: ( 0.33 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.890 ( 2.920 2.111 1.662) 3.967 1.912 ( 1.993 1.222 1.638) 2.855 1.947 ( 4.466 2.660 -1.350) 5.371 1.975 ( -1.450 -1.425 -2.029) 2.873 2.384 ( -7.637 14.177 0.322) 16.106 2.403 ( -7.589 14.251 -0.420) 16.152 2.947 ( -6.583 15.454 0.076) 16.798 2.959 ( -6.564 15.631 0.068) 16.953 3.545 ( 3.496 2.110 1.277) 4.279 3.577 ( 3.315 1.925 -1.257) 4.034 4.436 ( 10.523 1.555 1.282) 10.714 4.476 ( 10.614 0.846 -1.288) 10.726 5.964 ( 2.913 -0.274 0.360) 2.948 5.975 ( 3.158 -0.381 -0.257) 3.191 6.166 ( 3.737 5.822 0.356) 6.927 6.186 ( 3.745 5.955 -0.558) 7.057 6.306 ( -1.178 5.519 0.534) 5.669 6.320 ( -1.158 5.441 -0.422) 5.579 6.837 ( -0.443 -7.512 0.253) 7.530 6.863 ( -1.864 -6.988 0.072) 7.232 6.896 ( -1.708 -6.836 -0.242) 7.051 6.914 ( -1.554 -7.500 -0.058) 7.659 8.765 ( 1.500 -2.143 0.443) 2.653 8.777 ( 1.387 -2.297 -0.443) 2.719 ======================= Grid point 278 (47/54) ======================= q-point: ( 0.40 0.13 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.851 ( -0.367 -1.950 1.690) 2.606 1.893 ( -0.985 -2.150 -1.858) 3.007 2.037 ( 2.299 3.014 1.139) 3.959 2.071 ( 2.108 2.294 -0.962) 3.260 2.380 ( -9.400 15.569 0.221) 18.188 2.400 ( -9.405 15.540 -0.297) 18.167 2.968 ( -8.086 14.673 0.048) 16.754 2.982 ( -8.146 14.832 0.041) 16.922 3.603 ( 0.948 1.292 1.059) 1.921 3.631 ( 0.924 1.049 -1.089) 1.773 4.582 ( 5.959 -2.627 0.694) 6.549 4.614 ( 5.773 -3.054 -0.712) 6.570 5.992 ( 0.927 -0.411 0.501) 1.131 6.007 ( 1.288 -0.737 -0.469) 1.557 6.285 ( 2.624 3.587 0.123) 4.446 6.303 ( 2.326 3.488 -0.282) 4.202 6.340 ( -2.224 6.164 0.429) 6.567 6.353 ( -2.178 5.956 -0.279) 6.348 6.767 ( 1.562 -6.688 0.042) 6.868 6.778 ( 0.568 -5.877 0.225) 5.909 6.808 ( 0.427 -6.609 0.269) 6.628 6.821 ( 1.004 -8.201 -0.515) 8.279 8.772 ( 1.708 -2.431 0.132) 2.974 8.778 ( 1.516 -2.600 -0.136) 3.013 ======================= Grid point 290 (48/54) ======================= q-point: ( 0.20 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.634 ( 6.088 10.545 1.318) 12.247 1.671 ( 6.173 10.691 -1.702) 12.462 2.027 ( -3.067 -5.312 3.958) 7.300 2.117 ( -3.189 -5.523 -3.930) 7.491 2.419 ( 7.105 12.307 0.677) 14.227 2.440 ( 7.193 12.458 -0.619) 14.398 3.017 ( 6.916 11.978 0.215) 13.833 3.026 ( 7.087 12.275 -0.018) 14.174 3.365 ( 3.198 5.539 1.742) 6.629 3.406 ( 3.007 5.208 -1.653) 6.237 4.115 ( 5.041 8.731 0.895) 10.121 4.145 ( 5.016 8.687 -0.908) 10.072 5.854 ( 1.385 2.398 0.461) 2.807 5.866 ( 1.328 2.300 -0.412) 2.688 6.106 ( 1.422 2.462 0.442) 2.877 6.125 ( 1.481 2.565 -0.565) 3.015 6.312 ( 1.413 2.448 0.543) 2.879 6.326 ( 1.416 2.453 -0.453) 2.869 6.868 ( -2.431 -4.211 1.205) 5.009 6.913 ( -2.605 -4.512 -1.815) 5.518 6.990 ( -2.988 -5.175 0.705) 6.016 7.013 ( -3.110 -5.386 -0.110) 6.220 8.751 ( -1.001 -1.735 0.476) 2.059 8.763 ( -0.948 -1.643 -0.474) 1.955 ======================= Grid point 291 (49/54) ======================= q-point: ( 0.27 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.850 ( 6.368 8.583 1.220) 10.757 1.886 ( 6.271 7.884 -1.370) 10.167 1.918 ( -2.298 -4.097 3.083) 5.619 1.989 ( -3.611 -5.504 -3.132) 7.290 2.611 ( -0.370 12.576 0.460) 12.590 2.630 ( -0.457 12.502 -0.451) 12.519 3.187 ( -1.836 13.195 0.423) 13.328 3.202 ( -1.770 13.447 -0.138) 13.564 3.490 ( 5.130 3.006 1.278) 6.082 3.523 ( 4.886 2.709 -1.354) 5.748 4.299 ( 7.377 4.801 0.654) 8.826 4.328 ( 7.623 4.235 -0.676) 8.746 5.905 ( 2.577 0.403 0.410) 2.641 5.915 ( 2.513 0.282 -0.343) 2.552 6.175 ( 2.143 4.823 0.436) 5.296 6.196 ( 2.182 4.945 -0.590) 5.438 6.362 ( 0.614 3.941 0.565) 4.028 6.376 ( 0.620 3.996 -0.464) 4.070 6.774 ( -1.697 -6.834 1.009) 7.114 6.812 ( -1.940 -7.072 -1.576) 7.501 6.874 ( -3.239 -6.192 0.742) 7.027 6.896 ( -3.339 -6.127 -0.166) 6.980 8.727 ( 0.232 -1.781 0.462) 1.855 8.738 ( 0.187 -1.853 -0.461) 1.918 ======================= Grid point 292 (50/54) ======================= q-point: ( 0.33 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.840 ( -1.101 -2.632 2.014) 3.492 1.888 ( -1.871 -3.666 -2.114) 4.627 2.032 ( 3.835 6.564 0.870) 7.651 2.059 ( 3.781 5.756 -0.821) 6.936 2.679 ( -5.948 12.944 0.218) 14.247 2.697 ( -5.938 12.846 -0.268) 14.155 3.247 ( -5.983 12.985 0.361) 14.302 3.264 ( -6.037 13.162 -0.154) 14.481 3.587 ( 2.769 1.816 0.939) 3.442 3.613 ( 2.583 1.386 -1.033) 3.108 4.454 ( 6.330 0.505 0.467) 6.367 4.481 ( 6.501 0.036 -0.496) 6.520 5.948 ( 2.536 -1.157 0.343) 2.808 5.957 ( 2.625 -1.284 -0.270) 2.935 6.274 ( 3.121 4.520 0.360) 5.504 6.296 ( 2.997 4.564 -0.535) 5.486 6.417 ( 0.330 5.015 0.539) 5.055 6.431 ( 0.166 5.031 -0.427) 5.051 6.674 ( -0.640 -7.718 0.862) 7.792 6.705 ( -1.068 -7.851 -1.301) 8.029 6.756 ( -2.463 -6.218 0.782) 6.734 6.775 ( -2.817 -6.017 -0.316) 6.651 8.720 ( 1.135 -2.133 0.263) 2.430 8.728 ( 0.955 -2.268 -0.265) 2.476 ======================= Grid point 293 (51/54) ======================= q-point: ( 0.40 0.20 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.818 ( 0.649 -1.124 1.368) 1.886 1.851 ( 0.965 -1.672 -1.427) 2.401 2.103 ( -1.844 3.193 0.458) 3.716 2.125 ( -1.607 2.784 -0.382) 3.237 2.689 ( -7.599 13.162 0.134) 15.198 2.706 ( -7.552 13.080 -0.207) 15.105 3.258 ( -7.317 12.674 0.305) 14.638 3.275 ( -7.398 12.814 -0.155) 14.797 3.622 ( -0.236 0.409 0.692) 0.838 3.642 ( 0.001 -0.002 -0.785) 0.785 4.516 ( 1.905 -3.300 0.364) 3.828 4.542 ( 2.028 -3.513 -0.391) 4.075 5.965 ( 1.170 -2.027 0.370) 2.370 5.975 ( 1.332 -2.307 -0.309) 2.682 6.341 ( -0.935 1.619 -0.089) 1.871 6.357 ( -0.939 1.627 0.084) 1.880 6.455 ( -2.534 4.389 0.213) 5.072 6.463 ( -2.398 4.154 -0.273) 4.804 6.624 ( 3.397 -5.884 0.714) 6.832 6.644 ( 3.568 -6.180 -0.892) 7.192 6.692 ( 2.170 -3.758 0.520) 4.371 6.705 ( 2.173 -3.764 -0.297) 4.356 8.720 ( 1.374 -2.380 0.015) 2.748 8.725 ( 1.358 -2.352 -0.019) 2.715 ======================= Grid point 307 (52/54) ======================= q-point: ( 0.27 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 184 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.836 ( -1.554 -2.691 2.121) 3.763 1.886 ( -2.314 -4.007 -2.166) 5.109 2.032 ( 4.614 7.992 0.609) 9.249 2.055 ( 4.338 7.514 -0.624) 8.698 2.807 ( 3.653 6.327 0.246) 7.310 2.824 ( 3.606 6.245 -0.307) 7.218 3.389 ( 3.590 6.218 0.473) 7.196 3.406 ( 3.590 6.218 -0.040) 7.181 3.557 ( 1.983 3.434 0.826) 4.051 3.582 ( 1.712 2.965 -1.107) 3.598 4.396 ( 2.531 4.384 0.053) 5.063 4.417 ( 2.455 4.251 -0.098) 4.910 5.912 ( 0.292 0.506 0.356) 0.684 5.920 ( 0.222 0.385 -0.289) 0.530 6.269 ( 2.507 4.343 0.454) 5.035 6.292 ( 2.529 4.381 -0.637) 5.098 6.437 ( 1.997 3.459 0.646) 4.046 6.453 ( 2.016 3.492 -0.512) 4.065 6.641 ( -3.440 -5.958 1.001) 6.952 6.674 ( -3.603 -6.241 -1.538) 7.368 6.750 ( -3.297 -5.711 1.107) 6.687 6.773 ( -3.293 -5.704 -0.555) 6.610 8.697 ( -0.616 -1.067 0.320) 1.273 8.707 ( -0.704 -1.220 -0.321) 1.445 ======================= Grid point 308 (53/54) ======================= q-point: ( 0.33 0.27 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 338 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.804 ( -0.160 -0.858 1.047) 1.363 1.828 ( -0.568 -2.030 -1.034) 2.348 2.155 ( 1.072 4.747 0.144) 4.869 2.171 ( 1.056 4.530 -0.117) 4.653 2.883 ( -2.025 6.581 0.110) 6.886 2.899 ( -2.022 6.552 -0.217) 6.860 3.458 ( -2.777 7.231 0.296) 7.752 3.475 ( -2.688 7.128 0.025) 7.618 3.616 ( 1.548 0.911 0.332) 1.827 3.630 ( 1.310 0.221 -0.554) 1.439 4.464 ( 1.909 0.528 0.082) 1.983 4.486 ( 2.056 0.374 -0.136) 2.095 5.924 ( 1.083 -0.934 0.257) 1.453 5.929 ( 1.078 -1.135 -0.189) 1.576 6.350 ( 1.992 2.976 0.313) 3.595 6.371 ( 1.574 2.680 -0.374) 3.130 6.500 ( 0.652 1.957 0.302) 2.085 6.512 ( 0.489 2.327 -0.554) 2.441 6.547 ( -0.196 -3.549 0.962) 3.682 6.568 ( -0.897 -4.540 -0.955) 4.726 6.650 ( -1.692 -3.967 0.896) 4.405 6.669 ( -2.225 -4.217 -0.619) 4.808 8.683 ( 0.391 -1.307 0.096) 1.368 8.689 ( 0.289 -1.403 -0.099) 1.436 ======================= Grid point 323 (54/54) ======================= q-point: ( 0.33 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.97e-04 1.97e-04 1.97e-04 1.97e-04 4.63e-04 4.63e-04 4.63e-04 4.63e-04 Number of triplets: 73 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.801 ( -0.000 -0.000 0.096) 0.096 1.804 ( -0.000 -0.000 -0.051) 0.051 2.206 ( -0.000 -0.000 0.057) 0.057 2.222 ( -0.000 -0.000 -0.027) 0.027 2.949 ( -0.000 -0.000 0.086) 0.086 2.965 ( -0.000 -0.000 -0.219) 0.219 3.535 ( -0.000 -0.000 -0.023) 0.023 3.550 ( -0.000 -0.000 0.468) 0.468 3.625 ( -0.000 -0.000 -0.322) 0.322 3.628 ( 0.000 0.000 -0.018) 0.018 4.471 ( -0.000 -0.000 0.003) 0.003 4.491 ( -0.000 -0.000 -0.071) 0.071 5.915 ( 0.000 0.000 0.066) 0.066 5.916 ( -0.000 -0.000 0.004) 0.004 6.390 ( -0.000 -0.000 -0.024) 0.024 6.399 ( -0.000 -0.000 0.439) 0.439 6.478 ( -0.000 -0.000 -1.181) 1.181 6.494 ( -0.000 -0.000 -0.003) 0.003 6.551 ( 0.000 0.000 -0.002) 0.002 6.565 ( -0.000 -0.000 1.339) 1.339 6.607 ( 0.000 0.000 0.004) 0.004 6.612 ( -0.000 -0.000 -0.600) 0.600 8.669 ( -0.000 -0.000 -0.002) 0.002 8.674 ( -0.000 -0.000 -0.001) 0.001 =================== 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 290.190 290.190 25.685 -0.000 -0.000 0.000 3/16200 20.0 233.749 233.749 12.432 -0.000 0.000 0.000 3/16200 30.0 188.703 188.703 7.981 -0.000 -0.000 0.000 3/16200 40.0 156.625 156.625 5.862 -0.000 -0.000 0.000 3/16200 50.0 130.663 130.663 4.608 -0.000 -0.000 0.000 3/16200 60.0 109.770 109.770 3.779 -0.000 -0.000 0.000 3/16200 70.0 93.390 93.390 3.195 -0.000 -0.000 0.000 3/16200 80.0 80.673 80.673 2.766 -0.000 -0.000 0.000 3/16200 90.0 70.746 70.746 2.439 -0.000 -0.000 0.000 3/16200 100.0 62.887 62.887 2.183 -0.000 -0.000 0.000 3/16200 110.0 56.562 56.562 1.976 -0.000 -0.000 0.000 3/16200 120.0 51.386 51.386 1.806 -0.000 -0.000 0.000 3/16200 130.0 47.083 47.083 1.663 -0.000 0.000 0.000 3/16200 140.0 43.455 43.455 1.542 -0.000 0.000 0.000 3/16200 150.0 40.356 40.356 1.438 -0.000 0.000 0.000 3/16200 160.0 37.680 37.680 1.348 -0.000 0.000 0.000 3/16200 170.0 35.345 35.345 1.268 -0.000 0.000 0.000 3/16200 180.0 33.290 33.290 1.197 -0.000 0.000 0.000 3/16200 190.0 31.467 31.467 1.134 -0.000 0.000 0.000 3/16200 200.0 29.839 29.839 1.078 -0.000 0.000 0.000 3/16200 210.0 28.376 28.376 1.027 -0.000 0.000 0.000 3/16200 220.0 27.053 27.053 0.980 -0.000 0.000 0.000 3/16200 230.0 25.851 25.851 0.938 -0.000 0.000 0.000 3/16200 240.0 24.754 24.754 0.899 -0.000 0.000 0.000 3/16200 250.0 23.748 23.748 0.864 -0.000 0.000 0.000 3/16200 260.0 22.823 22.823 0.831 -0.000 0.000 0.000 3/16200 270.0 21.969 21.969 0.800 -0.000 0.000 0.000 3/16200 280.0 21.178 21.178 0.772 -0.000 0.000 0.000 3/16200 290.0 20.443 20.443 0.746 -0.000 0.000 0.000 3/16200 300.0 19.758 19.758 0.721 -0.000 0.000 0.000 3/16200 310.0 19.119 19.119 0.698 -0.000 0.000 0.000 3/16200 320.0 18.520 18.520 0.676 -0.000 0.000 0.000 3/16200 330.0 17.959 17.959 0.656 -0.000 0.000 0.000 3/16200 340.0 17.431 17.431 0.637 -0.000 0.000 0.000 3/16200 350.0 16.933 16.933 0.619 -0.000 0.000 0.000 3/16200 360.0 16.464 16.464 0.602 -0.000 0.000 0.000 3/16200 370.0 16.021 16.021 0.586 -0.000 0.000 0.000 3/16200 380.0 15.601 15.601 0.571 -0.000 0.000 0.000 3/16200 390.0 15.202 15.202 0.556 -0.000 0.000 0.000 3/16200 400.0 14.824 14.824 0.543 -0.000 0.000 0.000 3/16200 410.0 14.465 14.465 0.529 -0.000 0.000 0.000 3/16200 420.0 14.122 14.122 0.517 -0.000 0.000 0.000 3/16200 430.0 13.796 13.796 0.505 -0.000 0.000 0.000 3/16200 440.0 13.485 13.485 0.494 -0.000 0.000 0.000 3/16200 450.0 13.187 13.187 0.483 -0.000 0.000 0.000 3/16200 460.0 12.903 12.903 0.473 -0.000 0.000 0.000 3/16200 470.0 12.630 12.630 0.463 -0.000 0.000 0.000 3/16200 480.0 12.369 12.369 0.453 -0.000 0.000 0.000 3/16200 490.0 12.119 12.119 0.444 -0.000 0.000 0.000 3/16200 500.0 11.879 11.879 0.435 -0.000 0.000 0.000 3/16200 510.0 11.648 11.648 0.427 -0.000 0.000 0.000 3/16200 520.0 11.426 11.426 0.419 -0.000 0.000 0.000 3/16200 530.0 11.212 11.212 0.411 -0.000 0.000 0.000 3/16200 540.0 11.006 11.006 0.403 -0.000 0.000 0.000 3/16200 550.0 10.808 10.808 0.396 -0.000 0.000 0.000 3/16200 560.0 10.617 10.617 0.389 -0.000 0.000 0.000 3/16200 570.0 10.432 10.432 0.382 -0.000 0.000 0.000 3/16200 580.0 10.254 10.254 0.376 -0.000 0.000 0.000 3/16200 590.0 10.082 10.082 0.369 -0.000 0.000 0.000 3/16200 600.0 9.916 9.916 0.363 -0.000 0.000 0.000 3/16200 610.0 9.755 9.755 0.357 -0.000 0.000 0.000 3/16200 620.0 9.599 9.599 0.352 -0.000 0.000 0.000 3/16200 630.0 9.448 9.448 0.346 -0.000 0.000 0.000 3/16200 640.0 9.302 9.302 0.341 -0.000 0.000 0.000 3/16200 650.0 9.160 9.160 0.336 -0.000 0.000 0.000 3/16200 660.0 9.023 9.023 0.331 -0.000 0.000 0.000 3/16200 670.0 8.890 8.890 0.326 -0.000 0.000 0.000 3/16200 680.0 8.760 8.760 0.321 -0.000 0.000 0.000 3/16200 690.0 8.634 8.634 0.316 -0.000 0.000 0.000 3/16200 700.0 8.512 8.512 0.312 -0.000 0.000 0.000 3/16200 710.0 8.394 8.394 0.308 -0.000 0.000 0.000 3/16200 720.0 8.278 8.278 0.303 -0.000 0.000 0.000 3/16200 730.0 8.166 8.166 0.299 -0.000 0.000 0.000 3/16200 740.0 8.057 8.057 0.295 -0.000 0.000 0.000 3/16200 750.0 7.950 7.950 0.291 -0.000 0.000 0.000 3/16200 760.0 7.847 7.847 0.287 -0.000 0.000 0.000 3/16200 770.0 7.746 7.746 0.284 -0.000 0.000 0.000 3/16200 780.0 7.648 7.648 0.280 -0.000 0.000 0.000 3/16200 790.0 7.552 7.552 0.277 -0.000 0.000 0.000 3/16200 800.0 7.458 7.458 0.273 -0.000 0.000 0.000 3/16200 810.0 7.367 7.367 0.270 -0.000 0.000 0.000 3/16200 820.0 7.278 7.278 0.267 -0.000 0.000 0.000 3/16200 830.0 7.191 7.191 0.263 -0.000 0.000 0.000 3/16200 840.0 7.107 7.107 0.260 -0.000 0.000 0.000 3/16200 850.0 7.024 7.024 0.257 -0.000 0.000 0.000 3/16200 860.0 6.943 6.943 0.254 -0.000 0.000 0.000 3/16200 870.0 6.864 6.864 0.251 -0.000 0.000 0.000 3/16200 880.0 6.787 6.787 0.249 -0.000 0.000 0.000 3/16200 890.0 6.711 6.711 0.246 -0.000 0.000 0.000 3/16200 900.0 6.637 6.637 0.243 -0.000 0.000 0.000 3/16200 910.0 6.565 6.565 0.240 -0.000 0.000 0.000 3/16200 920.0 6.494 6.494 0.238 -0.000 0.000 0.000 3/16200 930.0 6.425 6.425 0.235 -0.000 0.000 0.000 3/16200 940.0 6.357 6.357 0.233 -0.000 0.000 0.000 3/16200 950.0 6.291 6.291 0.230 -0.000 0.000 0.000 3/16200 960.0 6.226 6.226 0.228 -0.000 0.000 0.000 3/16200 970.0 6.163 6.163 0.226 -0.000 0.000 0.000 3/16200 980.0 6.100 6.100 0.223 -0.000 0.000 0.000 3/16200 990.0 6.039 6.039 0.221 -0.000 0.000 0.000 3/16200 1000.0 5.979 5.979 0.219 -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-07 18:53:33]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|