----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 19:00:19]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [4 4 1] Primitive matrix: [ 1. -0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) Number of symmetry operations in supercell: 384 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.744748761320569 0.000000000000000 0.000000000000000 b -1.872374380660284 3.243047558093922 0.000000000000000 c 0.000000000000000 0.000000000000000 16.369318050000000 Atomic positions (fractional): *1 Ga 0.33333333333333 0.66666666666667 0.82414740655805 69.723 2 Ga 0.66666666666667 0.33333333333333 0.32414740655805 69.723 3 Ga 0.66666666666667 0.33333333333333 0.17585259344195 69.723 4 Ga 0.33333333333333 0.66666666666667 0.67585259344195 69.723 *5 Se 0.33333333333333 0.66666666666667 0.39585370278844 78.960 6 Se 0.66666666666667 0.33333333333333 0.89585370278844 78.960 7 Se 0.66666666666667 0.33333333333333 0.60414629721156 78.960 8 Se 0.33333333333333 0.66666666666667 0.10414629721156 78.960 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.744748761320569 0.000000000000000 0.000000000000000 b -1.872374380660284 3.243047558093922 0.000000000000000 c 0.000000000000000 0.000000000000000 16.369318050000000 Atomic positions (fractional): *1 Ga 0.33333333333333 0.66666666666667 0.82414740655805 69.723 > 1 2 Ga 0.66666666666667 0.33333333333333 0.32414740655805 69.723 > 2 3 Ga 0.66666666666667 0.33333333333333 0.17585259344195 69.723 > 3 4 Ga 0.33333333333333 0.66666666666667 0.67585259344195 69.723 > 4 *5 Se 0.33333333333333 0.66666666666667 0.39585370278844 78.960 > 5 6 Se 0.66666666666667 0.33333333333333 0.89585370278844 78.960 > 6 7 Se 0.66666666666667 0.33333333333333 0.60414629721156 78.960 > 7 8 Se 0.33333333333333 0.66666666666667 0.10414629721156 78.960 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 14.978995045282277 0.000000000000000 0.000000000000000 b -7.489497522641136 12.972190232375688 0.000000000000000 c 0.000000000000000 0.000000000000000 16.369318050000000 Atomic positions (fractional): *1 Ga 0.08333333333333 0.16666666666667 0.82414740655805 69.723 > 1 2 Ga 0.33333333333333 0.16666666666667 0.82414740655805 69.723 > 1 3 Ga 0.58333333333333 0.16666666666667 0.82414740655805 69.723 > 1 4 Ga 0.83333333333333 0.16666666666667 0.82414740655805 69.723 > 1 5 Ga 0.08333333333333 0.41666666666667 0.82414740655805 69.723 > 1 6 Ga 0.33333333333333 0.41666666666667 0.82414740655805 69.723 > 1 7 Ga 0.58333333333333 0.41666666666667 0.82414740655805 69.723 > 1 8 Ga 0.83333333333333 0.41666666666667 0.82414740655805 69.723 > 1 9 Ga 0.08333333333333 0.66666666666667 0.82414740655805 69.723 > 1 10 Ga 0.33333333333333 0.66666666666667 0.82414740655805 69.723 > 1 11 Ga 0.58333333333333 0.66666666666667 0.82414740655805 69.723 > 1 12 Ga 0.83333333333333 0.66666666666667 0.82414740655805 69.723 > 1 13 Ga 0.08333333333333 0.91666666666667 0.82414740655805 69.723 > 1 14 Ga 0.33333333333333 0.91666666666667 0.82414740655805 69.723 > 1 15 Ga 0.58333333333333 0.91666666666667 0.82414740655805 69.723 > 1 16 Ga 0.83333333333333 0.91666666666667 0.82414740655805 69.723 > 1 17 Ga 0.16666666666667 0.08333333333333 0.32414740655805 69.723 > 2 18 Ga 0.41666666666667 0.08333333333333 0.32414740655805 69.723 > 2 19 Ga 0.66666666666667 0.08333333333333 0.32414740655805 69.723 > 2 20 Ga 0.91666666666667 0.08333333333333 0.32414740655805 69.723 > 2 21 Ga 0.16666666666667 0.33333333333333 0.32414740655805 69.723 > 2 22 Ga 0.41666666666667 0.33333333333333 0.32414740655805 69.723 > 2 23 Ga 0.66666666666667 0.33333333333333 0.32414740655805 69.723 > 2 24 Ga 0.91666666666667 0.33333333333333 0.32414740655805 69.723 > 2 25 Ga 0.16666666666667 0.58333333333333 0.32414740655805 69.723 > 2 26 Ga 0.41666666666667 0.58333333333333 0.32414740655805 69.723 > 2 27 Ga 0.66666666666667 0.58333333333333 0.32414740655805 69.723 > 2 28 Ga 0.91666666666667 0.58333333333333 0.32414740655805 69.723 > 2 29 Ga 0.16666666666667 0.83333333333333 0.32414740655805 69.723 > 2 30 Ga 0.41666666666667 0.83333333333333 0.32414740655805 69.723 > 2 31 Ga 0.66666666666667 0.83333333333333 0.32414740655805 69.723 > 2 32 Ga 0.91666666666667 0.83333333333333 0.32414740655805 69.723 > 2 33 Ga 0.16666666666667 0.08333333333333 0.17585259344195 69.723 > 3 34 Ga 0.41666666666667 0.08333333333333 0.17585259344195 69.723 > 3 35 Ga 0.66666666666667 0.08333333333333 0.17585259344195 69.723 > 3 36 Ga 0.91666666666667 0.08333333333333 0.17585259344195 69.723 > 3 37 Ga 0.16666666666667 0.33333333333333 0.17585259344195 69.723 > 3 38 Ga 0.41666666666667 0.33333333333333 0.17585259344195 69.723 > 3 39 Ga 0.66666666666667 0.33333333333333 0.17585259344195 69.723 > 3 40 Ga 0.91666666666667 0.33333333333333 0.17585259344195 69.723 > 3 41 Ga 0.16666666666667 0.58333333333333 0.17585259344195 69.723 > 3 42 Ga 0.41666666666667 0.58333333333333 0.17585259344195 69.723 > 3 43 Ga 0.66666666666667 0.58333333333333 0.17585259344195 69.723 > 3 44 Ga 0.91666666666667 0.58333333333333 0.17585259344195 69.723 > 3 45 Ga 0.16666666666667 0.83333333333333 0.17585259344195 69.723 > 3 46 Ga 0.41666666666667 0.83333333333333 0.17585259344195 69.723 > 3 47 Ga 0.66666666666667 0.83333333333333 0.17585259344195 69.723 > 3 48 Ga 0.91666666666667 0.83333333333333 0.17585259344195 69.723 > 3 49 Ga 0.08333333333333 0.16666666666667 0.67585259344195 69.723 > 4 50 Ga 0.33333333333333 0.16666666666667 0.67585259344195 69.723 > 4 51 Ga 0.58333333333333 0.16666666666667 0.67585259344195 69.723 > 4 52 Ga 0.83333333333333 0.16666666666667 0.67585259344195 69.723 > 4 53 Ga 0.08333333333333 0.41666666666667 0.67585259344195 69.723 > 4 54 Ga 0.33333333333333 0.41666666666667 0.67585259344195 69.723 > 4 55 Ga 0.58333333333333 0.41666666666667 0.67585259344195 69.723 > 4 56 Ga 0.83333333333333 0.41666666666667 0.67585259344195 69.723 > 4 57 Ga 0.08333333333333 0.66666666666667 0.67585259344195 69.723 > 4 58 Ga 0.33333333333333 0.66666666666667 0.67585259344195 69.723 > 4 59 Ga 0.58333333333333 0.66666666666667 0.67585259344195 69.723 > 4 60 Ga 0.83333333333333 0.66666666666667 0.67585259344195 69.723 > 4 61 Ga 0.08333333333333 0.91666666666667 0.67585259344195 69.723 > 4 62 Ga 0.33333333333333 0.91666666666667 0.67585259344195 69.723 > 4 63 Ga 0.58333333333333 0.91666666666667 0.67585259344195 69.723 > 4 64 Ga 0.83333333333333 0.91666666666667 0.67585259344195 69.723 > 4 *65 Se 0.08333333333333 0.16666666666667 0.39585370278844 78.960 > 5 66 Se 0.33333333333333 0.16666666666667 0.39585370278844 78.960 > 5 67 Se 0.58333333333333 0.16666666666667 0.39585370278844 78.960 > 5 68 Se 0.83333333333333 0.16666666666667 0.39585370278844 78.960 > 5 69 Se 0.08333333333333 0.41666666666667 0.39585370278844 78.960 > 5 70 Se 0.33333333333333 0.41666666666667 0.39585370278844 78.960 > 5 71 Se 0.58333333333333 0.41666666666667 0.39585370278844 78.960 > 5 72 Se 0.83333333333333 0.41666666666667 0.39585370278844 78.960 > 5 73 Se 0.08333333333333 0.66666666666667 0.39585370278844 78.960 > 5 74 Se 0.33333333333333 0.66666666666667 0.39585370278844 78.960 > 5 75 Se 0.58333333333333 0.66666666666667 0.39585370278844 78.960 > 5 76 Se 0.83333333333333 0.66666666666667 0.39585370278844 78.960 > 5 77 Se 0.08333333333333 0.91666666666667 0.39585370278844 78.960 > 5 78 Se 0.33333333333333 0.91666666666667 0.39585370278844 78.960 > 5 79 Se 0.58333333333333 0.91666666666667 0.39585370278844 78.960 > 5 80 Se 0.83333333333333 0.91666666666667 0.39585370278844 78.960 > 5 81 Se 0.16666666666667 0.08333333333333 0.89585370278844 78.960 > 6 82 Se 0.41666666666667 0.08333333333333 0.89585370278844 78.960 > 6 83 Se 0.66666666666667 0.08333333333333 0.89585370278844 78.960 > 6 84 Se 0.91666666666667 0.08333333333333 0.89585370278844 78.960 > 6 85 Se 0.16666666666667 0.33333333333333 0.89585370278844 78.960 > 6 86 Se 0.41666666666667 0.33333333333333 0.89585370278844 78.960 > 6 87 Se 0.66666666666667 0.33333333333333 0.89585370278844 78.960 > 6 88 Se 0.91666666666667 0.33333333333333 0.89585370278844 78.960 > 6 89 Se 0.16666666666667 0.58333333333333 0.89585370278844 78.960 > 6 90 Se 0.41666666666667 0.58333333333333 0.89585370278844 78.960 > 6 91 Se 0.66666666666667 0.58333333333333 0.89585370278844 78.960 > 6 92 Se 0.91666666666667 0.58333333333333 0.89585370278844 78.960 > 6 93 Se 0.16666666666667 0.83333333333333 0.89585370278844 78.960 > 6 94 Se 0.41666666666667 0.83333333333333 0.89585370278844 78.960 > 6 95 Se 0.66666666666667 0.83333333333333 0.89585370278844 78.960 > 6 96 Se 0.91666666666667 0.83333333333333 0.89585370278844 78.960 > 6 97 Se 0.16666666666667 0.08333333333333 0.60414629721156 78.960 > 7 98 Se 0.41666666666667 0.08333333333333 0.60414629721156 78.960 > 7 99 Se 0.66666666666667 0.08333333333333 0.60414629721156 78.960 > 7 100 Se 0.91666666666667 0.08333333333333 0.60414629721156 78.960 > 7 101 Se 0.16666666666667 0.33333333333333 0.60414629721156 78.960 > 7 102 Se 0.41666666666667 0.33333333333333 0.60414629721156 78.960 > 7 103 Se 0.66666666666667 0.33333333333333 0.60414629721156 78.960 > 7 104 Se 0.91666666666667 0.33333333333333 0.60414629721156 78.960 > 7 105 Se 0.16666666666667 0.58333333333333 0.60414629721156 78.960 > 7 106 Se 0.41666666666667 0.58333333333333 0.60414629721156 78.960 > 7 107 Se 0.66666666666667 0.58333333333333 0.60414629721156 78.960 > 7 108 Se 0.91666666666667 0.58333333333333 0.60414629721156 78.960 > 7 109 Se 0.16666666666667 0.83333333333333 0.60414629721156 78.960 > 7 110 Se 0.41666666666667 0.83333333333333 0.60414629721156 78.960 > 7 111 Se 0.66666666666667 0.83333333333333 0.60414629721156 78.960 > 7 112 Se 0.91666666666667 0.83333333333333 0.60414629721156 78.960 > 7 113 Se 0.08333333333333 0.16666666666667 0.10414629721156 78.960 > 8 114 Se 0.33333333333333 0.16666666666667 0.10414629721156 78.960 > 8 115 Se 0.58333333333333 0.16666666666667 0.10414629721156 78.960 > 8 116 Se 0.83333333333333 0.16666666666667 0.10414629721156 78.960 > 8 117 Se 0.08333333333333 0.41666666666667 0.10414629721156 78.960 > 8 118 Se 0.33333333333333 0.41666666666667 0.10414629721156 78.960 > 8 119 Se 0.58333333333333 0.41666666666667 0.10414629721156 78.960 > 8 120 Se 0.83333333333333 0.41666666666667 0.10414629721156 78.960 > 8 121 Se 0.08333333333333 0.66666666666667 0.10414629721156 78.960 > 8 122 Se 0.33333333333333 0.66666666666667 0.10414629721156 78.960 > 8 123 Se 0.58333333333333 0.66666666666667 0.10414629721156 78.960 > 8 124 Se 0.83333333333333 0.66666666666667 0.10414629721156 78.960 > 8 125 Se 0.08333333333333 0.91666666666667 0.10414629721156 78.960 > 8 126 Se 0.33333333333333 0.91666666666667 0.10414629721156 78.960 > 8 127 Se 0.58333333333333 0.91666666666667 0.10414629721156 78.960 > 8 128 Se 0.83333333333333 0.91666666666667 0.10414629721156 78.960 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.2187833 -0.0000000 0.0000000 0.0000000 9.2187833 0.0000000 0.0000000 0.0000000 8.3707462 -------------------------- Born effective charges -------------------------- 1 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 2 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 3 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 4 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 5 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 6 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 7 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 8 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 ---------------------------------------------------------------------------- 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: 252 Number of blocks in projector: 192 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: 116 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 56 Use standard eigh solver. Tree of FC basis block matrices: - (252, 248), data: False |-- (56, 54), data: True |-- (80, 80), data: True |-- (116, 114), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.023 Solver_block: 80 / 80 - Time: 0.211 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.245 --------------------------------- 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: 252 Number of blocks in projector: 192 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: 116 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 56 Use standard eigh solver. Tree of FC basis block matrices: - (252, 248), data: False |-- (56, 54), data: True |-- (80, 80), data: True |-- (116, 114), 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 19:00:24]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 19:00:25]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 1] Primitive matrix: [ 1. -0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.744748761320569 0.000000000000000 0.000000000000000 b -1.872374380660284 3.243047558093922 0.000000000000000 c 0.000000000000000 0.000000000000000 16.369318050000000 Atomic positions (fractional): 1 Ga 0.33333333333333 0.66666666666667 0.82414740655805 69.723 2 Ga 0.66666666666667 0.33333333333333 0.32414740655805 69.723 3 Ga 0.66666666666667 0.33333333333333 0.17585259344195 69.723 4 Ga 0.33333333333333 0.66666666666667 0.67585259344195 69.723 5 Se 0.33333333333333 0.66666666666667 0.39585370278844 78.960 6 Se 0.66666666666667 0.33333333333333 0.89585370278844 78.960 7 Se 0.66666666666667 0.33333333333333 0.60414629721156 78.960 8 Se 0.33333333333333 0.66666666666667 0.10414629721156 78.960 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.978995045282277 0.000000000000000 0.000000000000000 b -7.489497522641136 12.972190232375688 0.000000000000000 c 0.000000000000000 0.000000000000000 16.369318050000000 Atomic positions (fractional): 1 Ga 0.08333333333333 0.16666666666667 0.82414740655805 69.723 > 1 2 Ga 0.33333333333333 0.16666666666667 0.82414740655805 69.723 > 1 3 Ga 0.58333333333333 0.16666666666667 0.82414740655805 69.723 > 1 4 Ga 0.83333333333333 0.16666666666667 0.82414740655805 69.723 > 1 5 Ga 0.08333333333333 0.41666666666667 0.82414740655805 69.723 > 1 6 Ga 0.33333333333333 0.41666666666667 0.82414740655805 69.723 > 1 7 Ga 0.58333333333333 0.41666666666667 0.82414740655805 69.723 > 1 8 Ga 0.83333333333333 0.41666666666667 0.82414740655805 69.723 > 1 9 Ga 0.08333333333333 0.66666666666667 0.82414740655805 69.723 > 1 10 Ga 0.33333333333333 0.66666666666667 0.82414740655805 69.723 > 1 11 Ga 0.58333333333333 0.66666666666667 0.82414740655805 69.723 > 1 12 Ga 0.83333333333333 0.66666666666667 0.82414740655805 69.723 > 1 13 Ga 0.08333333333333 0.91666666666667 0.82414740655805 69.723 > 1 14 Ga 0.33333333333333 0.91666666666667 0.82414740655805 69.723 > 1 15 Ga 0.58333333333333 0.91666666666667 0.82414740655805 69.723 > 1 16 Ga 0.83333333333333 0.91666666666667 0.82414740655805 69.723 > 1 17 Ga 0.16666666666667 0.08333333333333 0.32414740655805 69.723 > 17 18 Ga 0.41666666666667 0.08333333333333 0.32414740655805 69.723 > 17 19 Ga 0.66666666666667 0.08333333333333 0.32414740655805 69.723 > 17 20 Ga 0.91666666666667 0.08333333333333 0.32414740655805 69.723 > 17 21 Ga 0.16666666666667 0.33333333333333 0.32414740655805 69.723 > 17 22 Ga 0.41666666666667 0.33333333333333 0.32414740655805 69.723 > 17 23 Ga 0.66666666666667 0.33333333333333 0.32414740655805 69.723 > 17 24 Ga 0.91666666666667 0.33333333333333 0.32414740655805 69.723 > 17 25 Ga 0.16666666666667 0.58333333333333 0.32414740655805 69.723 > 17 26 Ga 0.41666666666667 0.58333333333333 0.32414740655805 69.723 > 17 27 Ga 0.66666666666667 0.58333333333333 0.32414740655805 69.723 > 17 28 Ga 0.91666666666667 0.58333333333333 0.32414740655805 69.723 > 17 29 Ga 0.16666666666667 0.83333333333333 0.32414740655805 69.723 > 17 30 Ga 0.41666666666667 0.83333333333333 0.32414740655805 69.723 > 17 31 Ga 0.66666666666667 0.83333333333333 0.32414740655805 69.723 > 17 32 Ga 0.91666666666667 0.83333333333333 0.32414740655805 69.723 > 17 33 Ga 0.16666666666667 0.08333333333333 0.17585259344195 69.723 > 33 34 Ga 0.41666666666667 0.08333333333333 0.17585259344195 69.723 > 33 35 Ga 0.66666666666667 0.08333333333333 0.17585259344195 69.723 > 33 36 Ga 0.91666666666667 0.08333333333333 0.17585259344195 69.723 > 33 37 Ga 0.16666666666667 0.33333333333333 0.17585259344195 69.723 > 33 38 Ga 0.41666666666667 0.33333333333333 0.17585259344195 69.723 > 33 39 Ga 0.66666666666667 0.33333333333333 0.17585259344195 69.723 > 33 40 Ga 0.91666666666667 0.33333333333333 0.17585259344195 69.723 > 33 41 Ga 0.16666666666667 0.58333333333333 0.17585259344195 69.723 > 33 42 Ga 0.41666666666667 0.58333333333333 0.17585259344195 69.723 > 33 43 Ga 0.66666666666667 0.58333333333333 0.17585259344195 69.723 > 33 44 Ga 0.91666666666667 0.58333333333333 0.17585259344195 69.723 > 33 45 Ga 0.16666666666667 0.83333333333333 0.17585259344195 69.723 > 33 46 Ga 0.41666666666667 0.83333333333333 0.17585259344195 69.723 > 33 47 Ga 0.66666666666667 0.83333333333333 0.17585259344195 69.723 > 33 48 Ga 0.91666666666667 0.83333333333333 0.17585259344195 69.723 > 33 49 Ga 0.08333333333333 0.16666666666667 0.67585259344195 69.723 > 49 50 Ga 0.33333333333333 0.16666666666667 0.67585259344195 69.723 > 49 51 Ga 0.58333333333333 0.16666666666667 0.67585259344195 69.723 > 49 52 Ga 0.83333333333333 0.16666666666667 0.67585259344195 69.723 > 49 53 Ga 0.08333333333333 0.41666666666667 0.67585259344195 69.723 > 49 54 Ga 0.33333333333333 0.41666666666667 0.67585259344195 69.723 > 49 55 Ga 0.58333333333333 0.41666666666667 0.67585259344195 69.723 > 49 56 Ga 0.83333333333333 0.41666666666667 0.67585259344195 69.723 > 49 57 Ga 0.08333333333333 0.66666666666667 0.67585259344195 69.723 > 49 58 Ga 0.33333333333333 0.66666666666667 0.67585259344195 69.723 > 49 59 Ga 0.58333333333333 0.66666666666667 0.67585259344195 69.723 > 49 60 Ga 0.83333333333333 0.66666666666667 0.67585259344195 69.723 > 49 61 Ga 0.08333333333333 0.91666666666667 0.67585259344195 69.723 > 49 62 Ga 0.33333333333333 0.91666666666667 0.67585259344195 69.723 > 49 63 Ga 0.58333333333333 0.91666666666667 0.67585259344195 69.723 > 49 64 Ga 0.83333333333333 0.91666666666667 0.67585259344195 69.723 > 49 65 Se 0.08333333333333 0.16666666666667 0.39585370278844 78.960 > 65 66 Se 0.33333333333333 0.16666666666667 0.39585370278844 78.960 > 65 67 Se 0.58333333333333 0.16666666666667 0.39585370278844 78.960 > 65 68 Se 0.83333333333333 0.16666666666667 0.39585370278844 78.960 > 65 69 Se 0.08333333333333 0.41666666666667 0.39585370278844 78.960 > 65 70 Se 0.33333333333333 0.41666666666667 0.39585370278844 78.960 > 65 71 Se 0.58333333333333 0.41666666666667 0.39585370278844 78.960 > 65 72 Se 0.83333333333333 0.41666666666667 0.39585370278844 78.960 > 65 73 Se 0.08333333333333 0.66666666666667 0.39585370278844 78.960 > 65 74 Se 0.33333333333333 0.66666666666667 0.39585370278844 78.960 > 65 75 Se 0.58333333333333 0.66666666666667 0.39585370278844 78.960 > 65 76 Se 0.83333333333333 0.66666666666667 0.39585370278844 78.960 > 65 77 Se 0.08333333333333 0.91666666666667 0.39585370278844 78.960 > 65 78 Se 0.33333333333333 0.91666666666667 0.39585370278844 78.960 > 65 79 Se 0.58333333333333 0.91666666666667 0.39585370278844 78.960 > 65 80 Se 0.83333333333333 0.91666666666667 0.39585370278844 78.960 > 65 81 Se 0.16666666666667 0.08333333333333 0.89585370278844 78.960 > 81 82 Se 0.41666666666667 0.08333333333333 0.89585370278844 78.960 > 81 83 Se 0.66666666666667 0.08333333333333 0.89585370278844 78.960 > 81 84 Se 0.91666666666667 0.08333333333333 0.89585370278844 78.960 > 81 85 Se 0.16666666666667 0.33333333333333 0.89585370278844 78.960 > 81 86 Se 0.41666666666667 0.33333333333333 0.89585370278844 78.960 > 81 87 Se 0.66666666666667 0.33333333333333 0.89585370278844 78.960 > 81 88 Se 0.91666666666667 0.33333333333333 0.89585370278844 78.960 > 81 89 Se 0.16666666666667 0.58333333333333 0.89585370278844 78.960 > 81 90 Se 0.41666666666667 0.58333333333333 0.89585370278844 78.960 > 81 91 Se 0.66666666666667 0.58333333333333 0.89585370278844 78.960 > 81 92 Se 0.91666666666667 0.58333333333333 0.89585370278844 78.960 > 81 93 Se 0.16666666666667 0.83333333333333 0.89585370278844 78.960 > 81 94 Se 0.41666666666667 0.83333333333333 0.89585370278844 78.960 > 81 95 Se 0.66666666666667 0.83333333333333 0.89585370278844 78.960 > 81 96 Se 0.91666666666667 0.83333333333333 0.89585370278844 78.960 > 81 97 Se 0.16666666666667 0.08333333333333 0.60414629721156 78.960 > 97 98 Se 0.41666666666667 0.08333333333333 0.60414629721156 78.960 > 97 99 Se 0.66666666666667 0.08333333333333 0.60414629721156 78.960 > 97 100 Se 0.91666666666667 0.08333333333333 0.60414629721156 78.960 > 97 101 Se 0.16666666666667 0.33333333333333 0.60414629721156 78.960 > 97 102 Se 0.41666666666667 0.33333333333333 0.60414629721156 78.960 > 97 103 Se 0.66666666666667 0.33333333333333 0.60414629721156 78.960 > 97 104 Se 0.91666666666667 0.33333333333333 0.60414629721156 78.960 > 97 105 Se 0.16666666666667 0.58333333333333 0.60414629721156 78.960 > 97 106 Se 0.41666666666667 0.58333333333333 0.60414629721156 78.960 > 97 107 Se 0.66666666666667 0.58333333333333 0.60414629721156 78.960 > 97 108 Se 0.91666666666667 0.58333333333333 0.60414629721156 78.960 > 97 109 Se 0.16666666666667 0.83333333333333 0.60414629721156 78.960 > 97 110 Se 0.41666666666667 0.83333333333333 0.60414629721156 78.960 > 97 111 Se 0.66666666666667 0.83333333333333 0.60414629721156 78.960 > 97 112 Se 0.91666666666667 0.83333333333333 0.60414629721156 78.960 > 97 113 Se 0.08333333333333 0.16666666666667 0.10414629721156 78.960 > 113 114 Se 0.33333333333333 0.16666666666667 0.10414629721156 78.960 > 113 115 Se 0.58333333333333 0.16666666666667 0.10414629721156 78.960 > 113 116 Se 0.83333333333333 0.16666666666667 0.10414629721156 78.960 > 113 117 Se 0.08333333333333 0.41666666666667 0.10414629721156 78.960 > 113 118 Se 0.33333333333333 0.41666666666667 0.10414629721156 78.960 > 113 119 Se 0.58333333333333 0.41666666666667 0.10414629721156 78.960 > 113 120 Se 0.83333333333333 0.41666666666667 0.10414629721156 78.960 > 113 121 Se 0.08333333333333 0.66666666666667 0.10414629721156 78.960 > 113 122 Se 0.33333333333333 0.66666666666667 0.10414629721156 78.960 > 113 123 Se 0.58333333333333 0.66666666666667 0.10414629721156 78.960 > 113 124 Se 0.83333333333333 0.66666666666667 0.10414629721156 78.960 > 113 125 Se 0.08333333333333 0.91666666666667 0.10414629721156 78.960 > 113 126 Se 0.33333333333333 0.91666666666667 0.10414629721156 78.960 > 113 127 Se 0.58333333333333 0.91666666666667 0.10414629721156 78.960 > 113 128 Se 0.83333333333333 0.91666666666667 0.10414629721156 78.960 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.2187833 -0.0000000 0.0000000 0.0000000 9.2187833 0.0000000 0.0000000 0.0000000 8.3707462 -------------------------- Born effective charges -------------------------- 1 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 2 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 3 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 4 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 5 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 6 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 7 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 8 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000002 (zzz) -0.00000002 (zzz) -0.00000002 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (yy) 0.00000000 (yy) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 19:00:31]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 19:00:32]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 1] Primitive matrix: [ 1. -0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.744748761320569 0.000000000000000 0.000000000000000 b -1.872374380660284 3.243047558093922 0.000000000000000 c 0.000000000000000 0.000000000000000 16.369318050000000 Atomic positions (fractional): 1 Ga 0.33333333333333 0.66666666666667 0.82414740655805 69.723 2 Ga 0.66666666666667 0.33333333333333 0.32414740655805 69.723 3 Ga 0.66666666666667 0.33333333333333 0.17585259344195 69.723 4 Ga 0.33333333333333 0.66666666666667 0.67585259344195 69.723 5 Se 0.33333333333333 0.66666666666667 0.39585370278844 78.960 6 Se 0.66666666666667 0.33333333333333 0.89585370278844 78.960 7 Se 0.66666666666667 0.33333333333333 0.60414629721156 78.960 8 Se 0.33333333333333 0.66666666666667 0.10414629721156 78.960 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.978995045282277 0.000000000000000 0.000000000000000 b -7.489497522641136 12.972190232375688 0.000000000000000 c 0.000000000000000 0.000000000000000 16.369318050000000 Atomic positions (fractional): 1 Ga 0.08333333333333 0.16666666666667 0.82414740655805 69.723 > 1 2 Ga 0.33333333333333 0.16666666666667 0.82414740655805 69.723 > 1 3 Ga 0.58333333333333 0.16666666666667 0.82414740655805 69.723 > 1 4 Ga 0.83333333333333 0.16666666666667 0.82414740655805 69.723 > 1 5 Ga 0.08333333333333 0.41666666666667 0.82414740655805 69.723 > 1 6 Ga 0.33333333333333 0.41666666666667 0.82414740655805 69.723 > 1 7 Ga 0.58333333333333 0.41666666666667 0.82414740655805 69.723 > 1 8 Ga 0.83333333333333 0.41666666666667 0.82414740655805 69.723 > 1 9 Ga 0.08333333333333 0.66666666666667 0.82414740655805 69.723 > 1 10 Ga 0.33333333333333 0.66666666666667 0.82414740655805 69.723 > 1 11 Ga 0.58333333333333 0.66666666666667 0.82414740655805 69.723 > 1 12 Ga 0.83333333333333 0.66666666666667 0.82414740655805 69.723 > 1 13 Ga 0.08333333333333 0.91666666666667 0.82414740655805 69.723 > 1 14 Ga 0.33333333333333 0.91666666666667 0.82414740655805 69.723 > 1 15 Ga 0.58333333333333 0.91666666666667 0.82414740655805 69.723 > 1 16 Ga 0.83333333333333 0.91666666666667 0.82414740655805 69.723 > 1 17 Ga 0.16666666666667 0.08333333333333 0.32414740655805 69.723 > 17 18 Ga 0.41666666666667 0.08333333333333 0.32414740655805 69.723 > 17 19 Ga 0.66666666666667 0.08333333333333 0.32414740655805 69.723 > 17 20 Ga 0.91666666666667 0.08333333333333 0.32414740655805 69.723 > 17 21 Ga 0.16666666666667 0.33333333333333 0.32414740655805 69.723 > 17 22 Ga 0.41666666666667 0.33333333333333 0.32414740655805 69.723 > 17 23 Ga 0.66666666666667 0.33333333333333 0.32414740655805 69.723 > 17 24 Ga 0.91666666666667 0.33333333333333 0.32414740655805 69.723 > 17 25 Ga 0.16666666666667 0.58333333333333 0.32414740655805 69.723 > 17 26 Ga 0.41666666666667 0.58333333333333 0.32414740655805 69.723 > 17 27 Ga 0.66666666666667 0.58333333333333 0.32414740655805 69.723 > 17 28 Ga 0.91666666666667 0.58333333333333 0.32414740655805 69.723 > 17 29 Ga 0.16666666666667 0.83333333333333 0.32414740655805 69.723 > 17 30 Ga 0.41666666666667 0.83333333333333 0.32414740655805 69.723 > 17 31 Ga 0.66666666666667 0.83333333333333 0.32414740655805 69.723 > 17 32 Ga 0.91666666666667 0.83333333333333 0.32414740655805 69.723 > 17 33 Ga 0.16666666666667 0.08333333333333 0.17585259344195 69.723 > 33 34 Ga 0.41666666666667 0.08333333333333 0.17585259344195 69.723 > 33 35 Ga 0.66666666666667 0.08333333333333 0.17585259344195 69.723 > 33 36 Ga 0.91666666666667 0.08333333333333 0.17585259344195 69.723 > 33 37 Ga 0.16666666666667 0.33333333333333 0.17585259344195 69.723 > 33 38 Ga 0.41666666666667 0.33333333333333 0.17585259344195 69.723 > 33 39 Ga 0.66666666666667 0.33333333333333 0.17585259344195 69.723 > 33 40 Ga 0.91666666666667 0.33333333333333 0.17585259344195 69.723 > 33 41 Ga 0.16666666666667 0.58333333333333 0.17585259344195 69.723 > 33 42 Ga 0.41666666666667 0.58333333333333 0.17585259344195 69.723 > 33 43 Ga 0.66666666666667 0.58333333333333 0.17585259344195 69.723 > 33 44 Ga 0.91666666666667 0.58333333333333 0.17585259344195 69.723 > 33 45 Ga 0.16666666666667 0.83333333333333 0.17585259344195 69.723 > 33 46 Ga 0.41666666666667 0.83333333333333 0.17585259344195 69.723 > 33 47 Ga 0.66666666666667 0.83333333333333 0.17585259344195 69.723 > 33 48 Ga 0.91666666666667 0.83333333333333 0.17585259344195 69.723 > 33 49 Ga 0.08333333333333 0.16666666666667 0.67585259344195 69.723 > 49 50 Ga 0.33333333333333 0.16666666666667 0.67585259344195 69.723 > 49 51 Ga 0.58333333333333 0.16666666666667 0.67585259344195 69.723 > 49 52 Ga 0.83333333333333 0.16666666666667 0.67585259344195 69.723 > 49 53 Ga 0.08333333333333 0.41666666666667 0.67585259344195 69.723 > 49 54 Ga 0.33333333333333 0.41666666666667 0.67585259344195 69.723 > 49 55 Ga 0.58333333333333 0.41666666666667 0.67585259344195 69.723 > 49 56 Ga 0.83333333333333 0.41666666666667 0.67585259344195 69.723 > 49 57 Ga 0.08333333333333 0.66666666666667 0.67585259344195 69.723 > 49 58 Ga 0.33333333333333 0.66666666666667 0.67585259344195 69.723 > 49 59 Ga 0.58333333333333 0.66666666666667 0.67585259344195 69.723 > 49 60 Ga 0.83333333333333 0.66666666666667 0.67585259344195 69.723 > 49 61 Ga 0.08333333333333 0.91666666666667 0.67585259344195 69.723 > 49 62 Ga 0.33333333333333 0.91666666666667 0.67585259344195 69.723 > 49 63 Ga 0.58333333333333 0.91666666666667 0.67585259344195 69.723 > 49 64 Ga 0.83333333333333 0.91666666666667 0.67585259344195 69.723 > 49 65 Se 0.08333333333333 0.16666666666667 0.39585370278844 78.960 > 65 66 Se 0.33333333333333 0.16666666666667 0.39585370278844 78.960 > 65 67 Se 0.58333333333333 0.16666666666667 0.39585370278844 78.960 > 65 68 Se 0.83333333333333 0.16666666666667 0.39585370278844 78.960 > 65 69 Se 0.08333333333333 0.41666666666667 0.39585370278844 78.960 > 65 70 Se 0.33333333333333 0.41666666666667 0.39585370278844 78.960 > 65 71 Se 0.58333333333333 0.41666666666667 0.39585370278844 78.960 > 65 72 Se 0.83333333333333 0.41666666666667 0.39585370278844 78.960 > 65 73 Se 0.08333333333333 0.66666666666667 0.39585370278844 78.960 > 65 74 Se 0.33333333333333 0.66666666666667 0.39585370278844 78.960 > 65 75 Se 0.58333333333333 0.66666666666667 0.39585370278844 78.960 > 65 76 Se 0.83333333333333 0.66666666666667 0.39585370278844 78.960 > 65 77 Se 0.08333333333333 0.91666666666667 0.39585370278844 78.960 > 65 78 Se 0.33333333333333 0.91666666666667 0.39585370278844 78.960 > 65 79 Se 0.58333333333333 0.91666666666667 0.39585370278844 78.960 > 65 80 Se 0.83333333333333 0.91666666666667 0.39585370278844 78.960 > 65 81 Se 0.16666666666667 0.08333333333333 0.89585370278844 78.960 > 81 82 Se 0.41666666666667 0.08333333333333 0.89585370278844 78.960 > 81 83 Se 0.66666666666667 0.08333333333333 0.89585370278844 78.960 > 81 84 Se 0.91666666666667 0.08333333333333 0.89585370278844 78.960 > 81 85 Se 0.16666666666667 0.33333333333333 0.89585370278844 78.960 > 81 86 Se 0.41666666666667 0.33333333333333 0.89585370278844 78.960 > 81 87 Se 0.66666666666667 0.33333333333333 0.89585370278844 78.960 > 81 88 Se 0.91666666666667 0.33333333333333 0.89585370278844 78.960 > 81 89 Se 0.16666666666667 0.58333333333333 0.89585370278844 78.960 > 81 90 Se 0.41666666666667 0.58333333333333 0.89585370278844 78.960 > 81 91 Se 0.66666666666667 0.58333333333333 0.89585370278844 78.960 > 81 92 Se 0.91666666666667 0.58333333333333 0.89585370278844 78.960 > 81 93 Se 0.16666666666667 0.83333333333333 0.89585370278844 78.960 > 81 94 Se 0.41666666666667 0.83333333333333 0.89585370278844 78.960 > 81 95 Se 0.66666666666667 0.83333333333333 0.89585370278844 78.960 > 81 96 Se 0.91666666666667 0.83333333333333 0.89585370278844 78.960 > 81 97 Se 0.16666666666667 0.08333333333333 0.60414629721156 78.960 > 97 98 Se 0.41666666666667 0.08333333333333 0.60414629721156 78.960 > 97 99 Se 0.66666666666667 0.08333333333333 0.60414629721156 78.960 > 97 100 Se 0.91666666666667 0.08333333333333 0.60414629721156 78.960 > 97 101 Se 0.16666666666667 0.33333333333333 0.60414629721156 78.960 > 97 102 Se 0.41666666666667 0.33333333333333 0.60414629721156 78.960 > 97 103 Se 0.66666666666667 0.33333333333333 0.60414629721156 78.960 > 97 104 Se 0.91666666666667 0.33333333333333 0.60414629721156 78.960 > 97 105 Se 0.16666666666667 0.58333333333333 0.60414629721156 78.960 > 97 106 Se 0.41666666666667 0.58333333333333 0.60414629721156 78.960 > 97 107 Se 0.66666666666667 0.58333333333333 0.60414629721156 78.960 > 97 108 Se 0.91666666666667 0.58333333333333 0.60414629721156 78.960 > 97 109 Se 0.16666666666667 0.83333333333333 0.60414629721156 78.960 > 97 110 Se 0.41666666666667 0.83333333333333 0.60414629721156 78.960 > 97 111 Se 0.66666666666667 0.83333333333333 0.60414629721156 78.960 > 97 112 Se 0.91666666666667 0.83333333333333 0.60414629721156 78.960 > 97 113 Se 0.08333333333333 0.16666666666667 0.10414629721156 78.960 > 113 114 Se 0.33333333333333 0.16666666666667 0.10414629721156 78.960 > 113 115 Se 0.58333333333333 0.16666666666667 0.10414629721156 78.960 > 113 116 Se 0.83333333333333 0.16666666666667 0.10414629721156 78.960 > 113 117 Se 0.08333333333333 0.41666666666667 0.10414629721156 78.960 > 113 118 Se 0.33333333333333 0.41666666666667 0.10414629721156 78.960 > 113 119 Se 0.58333333333333 0.41666666666667 0.10414629721156 78.960 > 113 120 Se 0.83333333333333 0.41666666666667 0.10414629721156 78.960 > 113 121 Se 0.08333333333333 0.66666666666667 0.10414629721156 78.960 > 113 122 Se 0.33333333333333 0.66666666666667 0.10414629721156 78.960 > 113 123 Se 0.58333333333333 0.66666666666667 0.10414629721156 78.960 > 113 124 Se 0.83333333333333 0.66666666666667 0.10414629721156 78.960 > 113 125 Se 0.08333333333333 0.91666666666667 0.10414629721156 78.960 > 113 126 Se 0.33333333333333 0.91666666666667 0.10414629721156 78.960 > 113 127 Se 0.58333333333333 0.91666666666667 0.10414629721156 78.960 > 113 128 Se 0.83333333333333 0.91666666666667 0.10414629721156 78.960 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 9.2187833 -0.0000000 0.0000000 0.0000000 9.2187833 0.0000000 0.0000000 0.0000000 8.3707462 -------------------------- Born effective charges -------------------------- 1 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 2 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 3 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 4 Ga 2.3282432 0.0000000 0.0000000 0.0000000 2.3282432 0.0000000 0.0000000 0.0000000 0.7877927 5 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 6 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 7 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 8 Se -2.3282432 -0.0000000 0.0000000 0.0000000 -2.3282432 0.0000000 0.0000000 0.0000000 -0.7877927 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000002 (zzz) -0.00000002 (zzz) -0.00000002 (zzz) 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.22 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.395 ( 0.000 0.000 0.000) 0.000 0.395 ( 0.000 0.000 0.000) 0.000 0.664 ( 0.000 0.000 0.000) 0.000 1.616 ( 0.000 0.000 0.000) 0.000 1.616 ( 0.000 0.000 0.000) 0.000 1.642 ( 0.000 0.000 0.000) 0.000 1.642 ( 0.000 0.000 0.000) 0.000 3.776 ( 0.000 0.000 0.000) 0.000 3.853 ( 0.000 0.000 0.000) 0.000 6.100 ( 0.000 0.000 0.000) 0.000 6.100 ( 0.000 0.000 0.000) 0.000 6.110 ( 0.000 0.000 0.000) 0.000 6.110 ( 0.000 0.000 0.000) 0.000 6.271 ( 0.000 0.000 0.000) 0.000 6.271 ( 0.000 0.000 0.000) 0.000 6.272 ( 0.000 0.000 0.000) 0.000 6.272 ( 0.000 0.000 0.000) 0.000 7.100 ( 0.000 0.000 0.000) 0.000 7.284 ( 0.000 0.000 0.000) 0.000 9.147 ( 0.000 0.000 0.000) 0.000 9.153 ( 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.247 ( 13.739 7.932 0.000) 15.865 0.470 ( 19.464 11.238 0.000) 22.475 0.620 ( 15.109 8.723 0.000) 17.446 0.667 ( 2.432 1.404 0.000) 2.808 0.866 ( 35.195 20.320 0.000) 40.640 0.966 ( 32.299 18.648 0.000) 37.295 1.694 ( 6.321 3.650 0.000) 7.299 1.721 ( 6.266 3.618 0.000) 7.235 1.894 ( 19.080 11.016 0.000) 22.032 1.912 ( 21.303 12.299 0.000) 24.598 3.742 ( -2.435 -1.406 0.000) 2.811 3.827 ( -2.070 -1.195 0.000) 2.391 6.087 ( -0.989 -0.571 0.000) 1.142 6.101 ( -0.762 -0.440 0.000) 0.880 6.121 ( 1.159 0.669 0.000) 1.338 6.147 ( 1.881 1.086 0.000) 2.173 6.257 ( -0.906 -0.523 0.000) 1.046 6.268 ( -0.561 -0.324 0.000) 0.648 6.441 ( 12.636 7.295 0.000) 14.590 7.040 ( -4.358 -2.516 0.000) 5.032 7.219 ( -4.701 -2.714 0.000) 5.428 7.308 ( -0.388 -0.224 0.000) 0.448 9.114 ( -2.688 -1.552 0.000) 3.104 9.119 ( -2.757 -1.592 0.000) 3.184 ======================= 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.620 ( 16.912 9.764 0.000) 19.528 0.818 ( 9.703 5.602 0.000) 11.204 0.914 ( 17.676 10.205 0.000) 20.410 1.000 ( 15.972 9.221 0.000) 18.443 1.634 ( 28.455 16.428 0.000) 32.857 1.686 ( 26.872 15.515 0.000) 31.030 1.893 ( 9.836 5.679 0.000) 11.357 1.916 ( 9.627 5.558 0.000) 11.116 2.423 ( 23.536 13.589 0.000) 27.177 2.484 ( 24.734 14.280 0.000) 28.560 3.690 ( -1.128 -0.651 0.000) 1.302 3.771 ( -2.043 -1.179 0.000) 2.359 6.057 ( -1.417 -0.818 0.000) 1.636 6.074 ( -1.505 -0.869 0.000) 1.738 6.101 ( -3.834 -2.214 0.000) 4.427 6.128 ( -4.422 -2.553 0.000) 5.106 6.231 ( -1.234 -0.712 0.000) 1.425 6.246 ( -1.280 -0.739 0.000) 1.479 6.768 ( 13.160 7.598 0.000) 15.196 6.953 ( -1.639 -0.946 0.000) 1.893 7.119 ( -2.596 -1.499 0.000) 2.997 7.292 ( -1.053 -0.608 0.000) 1.215 9.031 ( -4.033 -2.329 0.000) 4.657 9.032 ( -4.255 -2.457 0.000) 4.913 ======================= 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.026 ( 16.831 9.718 0.000) 19.435 1.094 ( 13.009 7.511 0.000) 15.021 1.296 ( 14.239 8.221 0.000) 16.442 1.352 ( 13.282 7.668 0.000) 15.337 2.134 ( 9.934 5.735 0.000) 11.471 2.141 ( 13.038 7.527 0.000) 15.055 2.152 ( 9.738 5.622 0.000) 11.244 2.206 ( 16.058 9.271 0.000) 18.542 2.941 ( 18.920 10.924 0.000) 21.847 3.017 ( 19.143 11.052 0.000) 22.105 3.737 ( 5.948 3.434 0.000) 6.868 3.772 ( 3.215 1.856 0.000) 3.712 5.952 ( -7.186 -4.149 0.000) 8.298 5.965 ( -7.389 -4.266 0.000) 8.532 6.023 ( -1.388 -0.801 0.000) 1.603 6.035 ( -1.632 -0.942 0.000) 1.884 6.201 ( -1.235 -0.713 0.000) 1.426 6.212 ( -1.418 -0.819 0.000) 1.637 6.974 ( 2.588 1.494 0.000) 2.989 7.015 ( 7.162 4.135 0.000) 8.270 7.109 ( 1.149 0.664 0.000) 1.327 7.255 ( -2.118 -1.223 0.000) 2.446 8.932 ( -3.846 -2.220 0.000) 4.441 8.938 ( -3.510 -2.027 0.000) 4.053 ======================= 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.385 ( 13.240 7.644 0.000) 15.289 1.432 ( 15.044 8.685 0.000) 17.371 1.579 ( 9.443 5.452 0.000) 10.904 1.619 ( 9.094 5.250 0.000) 10.501 2.272 ( -0.684 -0.395 0.000) 0.790 2.345 ( 7.556 4.362 0.000) 8.725 2.361 ( 7.612 4.395 0.000) 8.789 2.409 ( 1.504 0.868 0.000) 1.737 3.288 ( 10.192 5.884 0.000) 11.769 3.363 ( 10.019 5.784 0.000) 11.569 3.947 ( 11.235 6.486 0.000) 12.973 3.977 ( 13.116 7.573 0.000) 15.145 5.840 ( -0.981 -0.566 0.000) 1.132 5.863 ( -0.291 -0.168 0.000) 0.336 5.993 ( -1.047 -0.605 0.000) 1.209 6.004 ( -0.892 -0.515 0.000) 1.030 6.175 ( -0.885 -0.511 0.000) 1.021 6.185 ( -0.751 -0.434 0.000) 0.868 7.033 ( 1.719 0.993 0.000) 1.985 7.102 ( 0.260 0.150 0.000) 0.301 7.141 ( 1.062 0.613 0.000) 1.227 7.190 ( -3.289 -1.899 0.000) 3.798 8.861 ( -2.073 -1.197 0.000) 2.393 8.874 ( -1.736 -1.002 0.000) 2.005 ======================= 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.672 ( 10.786 6.227 0.000) 12.454 1.742 ( 4.497 2.596 0.000) 5.192 1.759 ( 12.008 6.933 0.000) 13.866 1.780 ( 4.540 2.621 0.000) 5.242 2.180 ( -5.826 -3.363 0.000) 6.727 2.340 ( -5.842 -3.373 0.000) 6.746 2.486 ( 4.343 2.508 0.000) 5.015 2.505 ( 4.533 2.617 0.000) 5.234 3.448 ( 4.240 2.448 0.000) 4.896 3.521 ( 4.207 2.429 0.000) 4.858 4.247 ( 12.964 7.485 0.000) 14.969 4.306 ( 13.604 7.854 0.000) 15.708 5.914 ( 6.463 3.731 0.000) 7.463 5.951 ( 6.979 4.030 0.000) 8.059 5.977 ( -0.279 -0.161 0.000) 0.322 5.993 ( -0.059 -0.034 0.000) 0.068 6.163 ( -0.091 -0.052 0.000) 0.105 6.177 ( 0.087 0.050 0.000) 0.100 7.043 ( -0.842 -0.486 0.000) 0.972 7.050 ( -3.875 -2.237 0.000) 4.474 7.099 ( -4.428 -2.556 0.000) 5.113 7.144 ( -0.882 -0.509 0.000) 1.018 8.833 ( -0.351 -0.202 0.000) 0.405 8.853 ( -0.191 -0.110 0.000) 0.220 ======================= 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.805 ( 1.182 0.682 0.000) 1.364 1.843 ( 1.057 0.610 0.000) 1.220 1.889 ( 7.387 4.265 0.000) 8.529 1.989 ( 7.216 4.166 0.000) 8.333 2.041 ( -5.177 -2.989 0.000) 5.978 2.184 ( -6.505 -3.756 0.000) 7.512 2.557 ( 1.909 1.102 0.000) 2.204 2.578 ( 1.768 1.021 0.000) 2.041 3.525 ( 2.598 1.500 0.000) 3.000 3.598 ( 2.617 1.511 0.000) 3.022 4.525 ( 9.898 5.715 0.000) 11.430 4.588 ( 9.570 5.525 0.000) 11.050 5.978 ( 0.254 0.147 0.000) 0.293 5.998 ( 0.379 0.219 0.000) 0.438 6.099 ( 8.143 4.702 0.000) 9.403 6.143 ( 8.077 4.663 0.000) 9.326 6.168 ( 0.433 0.250 0.000) 0.499 6.186 ( 0.559 0.323 0.000) 0.645 6.947 ( -4.258 -2.459 0.000) 4.917 6.990 ( -4.332 -2.501 0.000) 5.002 7.001 ( -2.410 -1.392 0.000) 2.783 7.101 ( -2.582 -1.491 0.000) 2.982 8.837 ( 0.446 0.258 0.000) 0.515 8.857 ( 0.311 0.180 0.000) 0.359 ======================= 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.817 ( 0.107 0.062 0.000) 0.123 1.849 ( -0.137 -0.079 0.000) 0.158 1.955 ( -1.836 -1.060 0.000) 2.120 2.013 ( 2.766 1.597 0.000) 3.194 2.065 ( -2.895 -1.672 0.000) 3.343 2.101 ( 2.268 1.309 0.000) 2.619 2.586 ( 0.609 0.352 0.000) 0.703 2.599 ( 0.224 0.129 0.000) 0.259 3.572 ( 1.119 0.646 0.000) 1.292 3.644 ( 1.102 0.636 0.000) 1.272 4.695 ( 3.974 2.294 0.000) 4.589 4.742 ( 3.234 1.867 0.000) 3.735 5.984 ( 0.184 0.106 0.000) 0.212 6.006 ( 0.224 0.129 0.000) 0.258 6.178 ( 0.262 0.151 0.000) 0.303 6.198 ( 0.313 0.180 0.000) 0.361 6.256 ( 4.149 2.395 0.000) 4.790 6.288 ( 3.326 1.920 0.000) 3.841 6.871 ( -1.811 -1.046 0.000) 2.092 6.912 ( -1.873 -1.081 0.000) 2.162 6.950 ( -1.393 -0.804 0.000) 1.609 7.042 ( -1.783 -1.029 0.000) 2.059 8.847 ( 0.302 0.174 0.000) 0.349 8.863 ( 0.137 0.079 0.000) 0.158 ======================= 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.509 ( 9.273 16.062 0.000) 18.547 0.757 ( 4.448 7.704 0.000) 8.896 0.846 ( 12.265 21.243 0.000) 24.530 0.963 ( 11.119 19.258 0.000) 22.238 1.427 ( 17.073 29.571 0.000) 34.145 1.456 ( 15.519 26.880 0.000) 31.038 1.852 ( 6.335 10.972 0.000) 12.670 1.877 ( 6.265 10.852 0.000) 12.530 2.260 ( 12.843 22.244 0.000) 25.685 2.308 ( 13.694 23.718 0.000) 27.387 3.699 ( -1.228 -2.128 0.000) 2.457 3.785 ( -1.435 -2.485 0.000) 2.870 6.060 ( -1.663 -2.880 0.000) 3.326 6.068 ( -1.108 -1.920 0.000) 2.216 6.142 ( -0.056 -0.096 0.000) 0.111 6.155 ( -0.684 -1.185 0.000) 1.369 6.250 ( -0.439 -0.761 0.000) 0.879 6.252 ( -0.372 -0.645 0.000) 0.744 6.678 ( 7.889 13.664 0.000) 15.778 6.965 ( -2.020 -3.498 0.000) 4.039 7.136 ( -2.431 -4.210 0.000) 4.861 7.295 ( -0.615 -1.065 0.000) 1.229 9.055 ( -2.262 -3.918 0.000) 4.524 9.059 ( -2.356 -4.080 0.000) 4.711 ======================= 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.868 ( 13.173 13.621 0.000) 18.949 0.971 ( 9.322 9.016 0.000) 12.969 1.265 ( 9.225 21.431 0.000) 23.333 1.338 ( 7.480 20.561 0.000) 21.879 1.954 ( 16.111 12.275 0.000) 20.255 1.992 ( 18.393 14.007 0.000) 23.119 2.107 ( 6.705 14.387 0.000) 15.872 2.127 ( 6.461 14.151 0.000) 15.556 2.734 ( 17.019 15.341 0.000) 22.912 2.804 ( 17.568 15.612 0.000) 23.502 3.692 ( 2.022 1.475 0.000) 2.503 3.752 ( 0.321 -0.142 0.000) 0.350 5.972 ( -3.106 -5.943 0.000) 6.706 5.992 ( -3.418 -6.018 0.000) 6.921 6.094 ( -3.801 -0.141 0.000) 3.804 6.104 ( -3.525 -0.018 0.000) 3.525 6.230 ( -1.409 0.470 0.000) 1.485 6.238 ( -1.081 0.116 0.000) 1.087 6.935 ( 7.489 7.532 0.000) 10.621 6.939 ( 1.399 -0.222 0.000) 1.417 7.085 ( 0.274 -1.573 0.000) 1.596 7.263 ( -1.158 -2.192 0.000) 2.479 8.964 ( -3.208 -3.955 0.000) 5.093 8.966 ( -2.944 -3.741 0.000) 4.761 ======================= 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.233 ( 12.934 9.765 0.000) 16.206 1.260 ( 13.337 8.865 0.000) 16.014 1.596 ( 4.200 19.215 0.000) 19.669 1.639 ( 3.501 18.418 0.000) 18.748 2.219 ( 6.396 0.243 0.000) 6.400 2.325 ( 8.943 1.942 0.000) 9.151 2.358 ( 3.742 15.271 0.000) 15.723 2.372 ( 3.682 14.995 0.000) 15.440 3.134 ( 13.839 7.415 0.000) 15.700 3.209 ( 13.894 7.233 0.000) 15.664 3.831 ( 8.464 5.279 0.000) 9.975 3.838 ( 8.475 4.996 0.000) 9.838 5.859 ( -2.742 -3.338 0.000) 4.320 5.870 ( -2.282 -3.414 0.000) 4.106 6.041 ( -3.149 1.495 0.000) 3.486 6.047 ( -3.113 1.225 0.000) 3.345 6.211 ( -1.882 1.588 0.000) 2.462 6.219 ( -1.767 1.547 0.000) 2.349 6.984 ( 3.145 -0.739 0.000) 3.231 7.075 ( 2.969 1.282 0.000) 3.234 7.096 ( 2.325 -1.956 0.000) 3.038 7.207 ( -1.912 -3.346 0.000) 3.854 8.877 ( -2.241 -3.043 0.000) 3.779 8.887 ( -1.909 -2.791 0.000) 3.382 ======================= 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.530 ( 11.758 6.493 0.000) 13.432 1.598 ( 13.605 7.386 0.000) 15.480 1.801 ( -0.489 14.022 0.000) 14.031 1.832 ( -0.495 13.743 0.000) 13.752 2.229 ( -2.248 -2.267 0.000) 3.192 2.372 ( -1.109 -3.546 0.000) 3.715 2.542 ( -0.058 13.937 0.000) 13.938 2.556 ( 0.219 13.884 0.000) 13.885 3.372 ( 7.574 2.681 0.000) 8.035 3.444 ( 7.534 2.550 0.000) 7.954 4.080 ( 12.854 6.128 0.000) 14.240 4.123 ( 14.232 6.297 0.000) 15.563 5.853 ( 3.225 1.736 0.000) 3.662 5.879 ( 4.042 1.689 0.000) 4.381 6.010 ( -2.056 1.971 0.000) 2.848 6.024 ( -1.595 2.428 0.000) 2.905 6.198 ( -1.868 2.687 0.000) 3.272 6.211 ( -1.655 2.830 0.000) 3.279 7.011 ( 1.840 -3.119 0.000) 3.622 7.072 ( -1.759 -2.989 0.000) 3.468 7.108 ( 1.918 -3.836 0.000) 4.289 7.127 ( -2.782 -4.213 0.000) 5.049 8.827 ( -0.411 -2.016 0.000) 2.058 8.845 ( -0.163 -1.817 0.000) 1.825 ======================= 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.780 ( 9.631 4.725 0.000) 10.728 1.872 ( -2.675 5.887 0.000) 6.467 1.875 ( 10.224 4.591 0.000) 11.208 1.905 ( -3.362 7.576 0.000) 8.289 2.150 ( -6.305 3.599 0.000) 7.260 2.270 ( -6.033 -1.622 0.000) 6.247 2.646 ( -3.164 12.245 0.000) 12.647 2.664 ( -2.989 12.114 0.000) 12.477 3.489 ( 3.749 1.827 0.000) 4.170 3.560 ( 3.793 1.764 0.000) 4.183 4.370 ( 12.681 4.301 0.000) 13.390 4.427 ( 13.169 3.819 0.000) 13.712 5.960 ( 3.952 0.691 0.000) 4.012 5.980 ( 2.572 0.321 0.000) 2.591 6.039 ( 3.262 5.474 0.000) 6.372 6.072 ( 4.338 5.941 0.000) 7.356 6.201 ( -1.471 3.627 0.000) 3.914 6.219 ( -1.322 3.899 0.000) 4.117 6.980 ( -3.507 -4.536 0.000) 5.733 6.984 ( 0.129 -5.132 0.000) 5.134 7.022 ( -3.537 -4.798 0.000) 5.961 7.079 ( 0.440 -5.698 0.000) 5.715 8.816 ( 0.970 -1.447 0.000) 1.743 8.837 ( 0.958 -1.430 0.000) 1.721 ======================= 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.853 ( -1.449 -0.539 0.000) 1.546 1.892 ( -3.642 3.484 0.000) 5.040 1.958 ( 5.621 2.835 0.000) 6.295 2.047 ( 5.305 1.944 0.000) 5.650 2.095 ( -7.642 9.849 0.000) 12.466 2.158 ( -7.120 3.783 0.000) 8.062 2.691 ( -4.968 11.072 0.000) 12.136 2.706 ( -5.168 10.708 0.000) 11.890 3.557 ( 1.922 1.613 0.000) 2.509 3.628 ( 1.962 1.446 0.000) 2.438 4.593 ( 8.789 0.960 0.000) 8.841 4.644 ( 8.415 0.079 0.000) 8.415 5.988 ( 0.213 0.659 0.000) 0.693 6.006 ( 0.514 0.519 0.000) 0.731 6.196 ( 5.976 4.471 0.000) 7.464 6.216 ( -1.578 4.291 0.000) 4.572 6.236 ( 5.691 4.165 0.000) 7.052 6.239 ( -1.583 4.678 0.000) 4.939 6.878 ( -2.195 -4.221 0.000) 4.758 6.919 ( -2.189 -4.331 0.000) 4.853 6.924 ( 0.069 -5.909 0.000) 5.910 7.015 ( 0.017 -6.624 0.000) 6.624 8.824 ( 1.384 -1.427 0.000) 1.988 8.843 ( 1.172 -1.533 0.000) 1.930 ======================= 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.835 ( 0.594 -1.029 0.000) 1.188 1.873 ( -1.611 2.791 0.000) 3.223 2.025 ( 0.044 -0.077 0.000) 0.089 2.078 ( -6.969 12.070 0.000) 13.937 2.104 ( 0.377 -0.653 0.000) 0.754 2.111 ( -4.546 7.875 0.000) 9.093 2.703 ( -5.972 10.343 0.000) 11.944 2.712 ( -5.888 10.198 0.000) 11.776 3.584 ( -0.290 0.502 0.000) 0.580 3.655 ( -0.180 0.312 0.000) 0.360 4.681 ( 1.984 -3.436 0.000) 3.967 4.719 ( 2.183 -3.781 0.000) 4.366 5.994 ( -0.338 0.586 0.000) 0.676 6.013 ( -0.236 0.408 0.000) 0.471 6.226 ( -2.447 4.238 0.000) 4.894 6.251 ( -2.663 4.613 0.000) 5.327 6.285 ( -0.194 0.336 0.000) 0.388 6.312 ( -0.286 0.495 0.000) 0.571 6.834 ( 1.399 -2.423 0.000) 2.798 6.874 ( 1.459 -2.528 0.000) 2.919 6.890 ( 2.755 -4.772 0.000) 5.511 6.974 ( 3.049 -5.280 0.000) 6.097 8.830 ( 1.002 -1.736 0.000) 2.005 8.845 ( 0.992 -1.718 0.000) 1.984 ======================= 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.162 ( 8.489 14.703 0.000) 16.977 1.180 ( 6.434 11.143 0.000) 12.867 1.701 ( 11.110 19.243 0.000) 22.220 1.750 ( 10.419 18.046 0.000) 20.838 2.125 ( 2.757 4.775 0.000) 5.514 2.217 ( 4.130 7.153 0.000) 8.260 2.431 ( 9.060 15.692 0.000) 18.120 2.443 ( 8.823 15.282 0.000) 17.647 3.024 ( 7.272 12.596 0.000) 14.544 3.097 ( 7.274 12.599 0.000) 14.549 3.775 ( 3.986 6.904 0.000) 7.973 3.792 ( 2.700 4.677 0.000) 5.400 5.864 ( -2.175 -3.768 0.000) 4.350 5.873 ( -2.584 -4.475 0.000) 5.167 6.088 ( -0.293 -0.508 0.000) 0.587 6.093 ( -0.589 -1.021 0.000) 1.179 6.243 ( 0.405 0.701 0.000) 0.809 6.249 ( 0.467 0.810 0.000) 0.935 6.943 ( 0.255 0.441 0.000) 0.509 7.049 ( 1.882 3.260 0.000) 3.764 7.057 ( -0.678 -1.175 0.000) 1.356 7.204 ( -2.012 -3.485 0.000) 4.024 8.881 ( -2.237 -3.874 0.000) 4.473 8.889 ( -2.027 -3.511 0.000) 4.054 ======================= 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.420 ( 8.720 9.062 0.000) 12.576 1.465 ( 10.576 10.596 0.000) 14.971 1.979 ( -1.330 13.067 0.000) 13.135 2.029 ( -0.907 17.415 0.000) 17.439 2.201 ( 4.746 1.846 0.000) 5.093 2.301 ( 4.477 -2.758 0.000) 5.258 2.707 ( 2.242 17.317 0.000) 17.462 2.717 ( 2.628 17.208 0.000) 17.408 3.251 ( 8.414 4.606 0.000) 9.593 3.323 ( 8.395 4.474 0.000) 9.512 3.951 ( 8.121 6.781 0.000) 10.579 3.972 ( 9.360 7.439 0.000) 11.956 5.832 ( 0.610 0.560 0.000) 0.828 5.842 ( 1.499 0.485 0.000) 1.576 6.070 ( -1.772 1.288 0.000) 2.191 6.079 ( -1.173 1.908 0.000) 2.240 6.257 ( -0.740 2.626 0.000) 2.728 6.267 ( -0.518 2.766 0.000) 2.814 6.945 ( 0.993 -2.895 0.000) 3.060 7.032 ( 0.725 -4.017 0.000) 4.082 7.063 ( -0.782 -2.380 0.000) 2.505 7.122 ( -2.326 -4.770 0.000) 5.307 8.814 ( -1.163 -2.865 0.000) 3.092 8.830 ( -0.920 -2.563 0.000) 2.723 ======================= 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.667 ( 9.301 7.249 0.000) 11.792 1.746 ( 10.258 7.353 0.000) 12.621 1.959 ( -4.105 -1.366 0.000) 4.327 2.060 ( -5.552 2.991 0.000) 6.307 2.322 ( -1.995 13.024 0.000) 13.176 2.365 ( -0.993 7.572 0.000) 7.636 2.868 ( -3.365 16.673 0.000) 17.009 2.885 ( -3.222 16.885 0.000) 17.190 3.419 ( 6.429 2.348 0.000) 6.844 3.489 ( 6.450 2.212 0.000) 6.819 4.193 ( 11.385 4.954 0.000) 12.416 4.233 ( 12.174 4.558 0.000) 13.000 5.891 ( 3.970 1.483 0.000) 4.238 5.913 ( 4.397 1.179 0.000) 4.552 6.073 ( -0.314 4.152 0.000) 4.163 6.100 ( 0.216 4.716 0.000) 4.721 6.272 ( -1.309 4.159 0.000) 4.360 6.289 ( -1.000 4.456 0.000) 4.567 6.913 ( 0.569 -6.053 0.000) 6.080 6.985 ( -2.971 -5.249 0.000) 6.031 6.993 ( 0.986 -6.880 0.000) 6.950 7.023 ( -2.856 -5.686 0.000) 6.363 8.782 ( 0.417 -2.274 0.000) 2.312 8.803 ( 0.499 -2.165 0.000) 2.222 ======================= 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.869 ( -2.546 -3.247 0.000) 4.126 1.885 ( 7.157 5.792 0.000) 9.207 1.966 ( -3.930 -2.708 0.000) 4.773 1.967 ( 7.155 4.774 0.000) 8.601 2.377 ( -7.739 15.001 0.000) 16.880 2.395 ( -7.612 13.517 0.000) 15.513 2.944 ( -6.503 15.632 0.000) 16.931 2.960 ( -6.758 15.546 0.000) 16.951 3.531 ( 3.493 2.246 0.000) 4.153 3.599 ( 3.512 1.895 0.000) 3.991 4.428 ( 10.438 1.579 0.000) 10.557 4.470 ( 10.667 0.664 0.000) 10.687 5.965 ( 2.954 -0.219 0.000) 2.962 5.981 ( 2.605 -0.357 0.000) 2.629 6.160 ( 3.694 5.795 0.000) 6.872 6.196 ( 3.958 5.846 0.000) 7.060 6.298 ( -1.339 5.415 0.000) 5.578 6.324 ( -1.125 5.863 0.000) 5.970 6.845 ( 0.152 -7.968 0.000) 7.969 6.870 ( -2.961 -5.916 0.000) 6.615 6.909 ( -2.844 -5.983 0.000) 6.625 6.924 ( 0.518 -8.878 0.000) 8.893 8.777 ( 1.470 -2.224 0.000) 2.666 8.798 ( 1.323 -2.322 0.000) 2.673 ======================= 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.816 ( -0.263 -1.981 0.000) 1.999 1.891 ( -1.091 -2.348 0.000) 2.589 2.021 ( 2.202 3.407 0.000) 4.057 2.087 ( 2.155 2.202 0.000) 3.081 2.373 ( -9.201 14.976 0.000) 17.577 2.390 ( -9.744 16.345 0.000) 19.029 2.969 ( -8.002 14.887 0.000) 16.901 2.978 ( -8.294 14.713 0.000) 16.890 3.590 ( 0.891 1.480 0.000) 1.728 3.654 ( 1.049 0.938 0.000) 1.407 4.574 ( 5.927 -2.471 0.000) 6.421 4.606 ( 5.812 -3.290 0.000) 6.678 5.994 ( 0.895 -0.351 0.000) 0.961 6.009 ( 1.091 -0.540 0.000) 1.217 6.278 ( 2.717 3.372 0.000) 4.330 6.312 ( 2.283 3.214 0.000) 3.942 6.328 ( -2.314 6.059 0.000) 6.486 6.361 ( -2.361 6.586 0.000) 6.996 6.774 ( 2.011 -8.014 0.000) 8.263 6.781 ( -0.087 -4.950 0.000) 4.951 6.818 ( -0.119 -5.221 0.000) 5.222 6.847 ( 2.098 -9.152 0.000) 9.389 8.782 ( 1.672 -2.484 0.000) 2.995 8.798 ( 1.531 -2.605 0.000) 3.022 ======================= 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.619 ( 5.940 10.289 0.000) 11.880 1.688 ( 6.348 10.994 0.000) 12.695 1.988 ( -3.095 -5.361 0.000) 6.191 2.125 ( -3.412 -5.910 0.000) 6.824 2.419 ( 7.239 12.537 0.000) 14.477 2.425 ( 7.186 12.447 0.000) 14.372 3.020 ( 6.872 11.902 0.000) 13.743 3.027 ( 7.075 12.254 0.000) 14.149 3.350 ( 3.172 5.494 0.000) 6.344 3.419 ( 3.118 5.401 0.000) 6.237 4.109 ( 4.913 8.509 0.000) 9.825 4.138 ( 5.002 8.664 0.000) 10.004 5.863 ( 1.189 2.059 0.000) 2.377 5.876 ( 1.382 2.394 0.000) 2.764 6.103 ( 1.286 2.228 0.000) 2.573 6.127 ( 1.679 2.907 0.000) 3.357 6.309 ( 1.389 2.405 0.000) 2.777 6.324 ( 1.598 2.768 0.000) 3.196 6.871 ( -2.393 -4.144 0.000) 4.786 6.942 ( -2.681 -4.643 0.000) 5.362 6.984 ( -2.939 -5.091 0.000) 5.879 7.017 ( -3.109 -5.385 0.000) 6.218 8.765 ( -1.071 -1.856 0.000) 2.143 8.786 ( -0.979 -1.695 0.000) 1.957 ======================= 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.835 ( 6.749 8.984 0.000) 11.236 1.877 ( -2.865 -4.574 0.000) 5.397 1.907 ( 6.689 8.282 0.000) 10.646 1.989 ( -3.911 -5.896 0.000) 7.075 2.606 ( -0.614 12.054 0.000) 12.069 2.622 ( -0.222 13.193 0.000) 13.195 3.186 ( -1.974 13.287 0.000) 13.433 3.205 ( -1.812 13.438 0.000) 13.560 3.475 ( 5.249 3.029 0.000) 6.060 3.541 ( 5.037 2.765 0.000) 5.746 4.291 ( 7.395 4.763 0.000) 8.796 4.320 ( 7.570 4.139 0.000) 8.628 5.909 ( 2.402 0.366 0.000) 2.430 5.925 ( 2.438 0.141 0.000) 2.442 6.170 ( 2.086 4.898 0.000) 5.324 6.203 ( 2.265 4.973 0.000) 5.464 6.357 ( 0.446 3.891 0.000) 3.917 6.381 ( 0.812 4.291 0.000) 4.368 6.778 ( -1.696 -6.713 0.000) 6.924 6.843 ( -1.513 -7.348 0.000) 7.502 6.866 ( -3.483 -6.105 0.000) 7.029 6.902 ( -3.259 -5.925 0.000) 6.762 8.739 ( 0.202 -1.838 0.000) 1.850 8.760 ( 0.123 -1.936 0.000) 1.940 ======================= 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.805 ( -1.038 -2.459 0.000) 2.669 1.883 ( -1.957 -4.150 0.000) 4.588 2.019 ( 3.982 7.006 0.000) 8.058 2.076 ( 3.642 5.673 0.000) 6.742 2.664 ( -6.025 12.134 0.000) 13.547 2.701 ( -5.956 13.887 0.000) 15.110 3.247 ( -5.820 13.022 0.000) 14.264 3.263 ( -6.314 13.206 0.000) 14.637 3.576 ( 2.854 2.050 0.000) 3.514 3.632 ( 2.697 1.231 0.000) 2.965 4.448 ( 6.370 0.672 0.000) 6.406 4.471 ( 6.462 -0.184 0.000) 6.465 5.951 ( 2.538 -1.096 0.000) 2.764 5.963 ( 2.430 -1.290 0.000) 2.751 6.268 ( 2.902 4.565 0.000) 5.410 6.303 ( 3.010 4.419 0.000) 5.347 6.408 ( 0.066 4.975 0.000) 4.975 6.443 ( 0.411 5.380 0.000) 5.395 6.679 ( -0.742 -7.607 0.000) 7.643 6.738 ( -0.733 -8.712 0.000) 8.743 6.748 ( -2.458 -5.587 0.000) 6.104 6.785 ( -2.661 -5.847 0.000) 6.424 8.731 ( 1.110 -2.115 0.000) 2.389 8.748 ( 0.934 -2.375 0.000) 2.552 ======================= 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.785 ( 0.562 -0.973 0.000) 1.124 1.842 ( 1.158 -2.006 0.000) 2.317 2.095 ( -2.076 3.596 0.000) 4.152 2.139 ( -1.520 2.632 0.000) 3.040 2.668 ( -7.209 12.487 0.000) 14.419 2.716 ( -8.080 13.995 0.000) 16.160 3.260 ( -7.236 12.533 0.000) 14.471 3.272 ( -7.526 13.036 0.000) 15.053 3.613 ( -0.381 0.660 0.000) 0.762 3.662 ( 0.149 -0.257 0.000) 0.297 4.511 ( 1.812 -3.138 0.000) 3.624 4.530 ( 2.134 -3.696 0.000) 4.267 5.969 ( 1.149 -1.991 0.000) 2.299 5.979 ( 1.262 -2.187 0.000) 2.525 6.331 ( -0.931 1.612 0.000) 1.862 6.363 ( -0.919 1.592 0.000) 1.838 6.442 ( -2.550 4.417 0.000) 5.100 6.484 ( -2.708 4.690 0.000) 5.415 6.625 ( 3.323 -5.756 0.000) 6.647 6.672 ( 4.023 -6.968 0.000) 8.046 6.690 ( 1.952 -3.382 0.000) 3.905 6.719 ( 2.254 -3.905 0.000) 4.509 8.731 ( 1.330 -2.303 0.000) 2.660 8.744 ( 1.425 -2.468 0.000) 2.850 ======================= 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.802 ( -1.440 -2.495 0.000) 2.881 1.879 ( -2.553 -4.422 0.000) 5.106 2.021 ( 4.883 8.458 0.000) 9.767 2.073 ( 4.217 7.304 0.000) 8.434 2.788 ( 3.201 5.544 0.000) 6.402 2.834 ( 4.139 7.169 0.000) 8.278 3.385 ( 3.466 6.003 0.000) 6.932 3.415 ( 3.746 6.488 0.000) 7.491 3.542 ( 2.051 3.553 0.000) 4.103 3.601 ( 1.703 2.949 0.000) 3.405 4.389 ( 2.638 4.569 0.000) 5.276 4.406 ( 2.334 4.043 0.000) 4.668 5.916 ( 0.289 0.500 0.000) 0.578 5.928 ( 0.176 0.304 0.000) 0.351 6.266 ( 2.520 4.365 0.000) 5.040 6.298 ( 2.488 4.309 0.000) 4.975 6.431 ( 1.942 3.364 0.000) 3.884 6.464 ( 2.173 3.765 0.000) 4.347 6.647 ( -3.369 -5.835 0.000) 6.737 6.700 ( -3.753 -6.500 0.000) 7.506 6.745 ( -3.150 -5.456 0.000) 6.299 6.782 ( -3.254 -5.636 0.000) 6.508 8.709 ( -0.596 -1.033 0.000) 1.193 8.726 ( -0.782 -1.354 0.000) 1.563 ======================= 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.774 ( -0.102 -0.543 0.000) 0.553 1.814 ( -0.620 -2.475 0.000) 2.552 2.153 ( 1.215 5.420 0.000) 5.554 2.182 ( 0.910 4.048 0.000) 4.149 2.851 ( -2.037 5.871 0.000) 6.214 2.924 ( -2.049 7.428 0.000) 7.706 3.453 ( -2.446 6.741 0.000) 7.171 3.482 ( -3.224 7.733 0.000) 8.378 3.610 ( 1.926 1.147 0.000) 2.242 3.645 ( 1.176 -0.056 0.000) 1.177 4.462 ( 1.952 0.728 0.000) 2.084 4.471 ( 2.018 0.169 0.000) 2.025 5.928 ( 1.096 -0.874 0.000) 1.402 5.936 ( 1.016 -1.164 0.000) 1.545 6.344 ( 1.600 2.748 0.000) 3.180 6.376 ( 1.628 2.793 0.000) 3.233 6.494 ( 0.934 3.265 0.000) 3.396 6.536 ( 1.312 3.500 0.000) 3.738 6.551 ( -0.723 -4.553 0.000) 4.610 6.583 ( -1.460 -6.028 0.000) 6.203 6.653 ( -1.363 -3.488 0.000) 3.744 6.677 ( -2.190 -4.422 0.000) 4.934 8.696 ( 0.388 -1.167 0.000) 1.230 8.706 ( 0.280 -1.592 0.000) 1.617 ======================= Grid point 82 (27/54) ======================= q-point: ( 0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.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.779 ( -0.481 -0.832 0.000) 0.961 1.779 ( 0.481 0.832 0.000) 0.961 2.220 ( -0.353 -0.612 0.000) 0.706 2.220 ( 0.353 0.612 0.000) 0.706 2.909 ( -0.000 -0.000 0.000) 0.000 3.000 ( -0.000 -0.000 0.000) 0.000 3.523 ( -0.000 -0.000 0.000) 0.000 3.566 ( -0.000 -0.000 0.000) 0.000 3.630 ( -0.574 -0.994 0.000) 1.147 3.630 ( 0.574 0.994 0.000) 1.147 4.473 ( -0.089 -0.154 0.000) 0.177 4.473 ( 0.089 0.154 0.000) 0.177 5.921 ( -0.103 -0.178 0.000) 0.206 5.921 ( 0.103 0.178 0.000) 0.206 6.374 ( -0.000 -0.000 0.000) 0.000 6.413 ( -0.000 -0.000 0.000) 0.000 6.506 ( -0.398 -0.689 0.000) 0.796 6.506 ( 0.398 0.689 0.000) 0.796 6.533 ( -0.000 -0.000 0.000) 0.000 6.583 ( -0.000 -0.000 0.000) 0.000 6.615 ( -0.272 -0.471 0.000) 0.544 6.615 ( 0.272 0.471 0.000) 0.544 8.686 ( -0.124 -0.215 0.000) 0.248 8.686 ( 0.124 0.215 0.000) 0.248 ======================= 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.195 ( -0.000 0.000 8.745) 8.745 0.195 ( 0.000 0.000 8.745) 8.745 0.328 ( 0.000 0.000 14.721) 14.721 0.340 ( -0.000 0.000 -5.190) 5.190 0.340 ( -0.000 -0.000 -5.190) 5.190 0.573 ( 0.000 -0.000 -8.722) 8.722 1.623 ( -0.000 0.000 0.617) 0.617 1.623 ( 0.000 -0.000 0.617) 0.617 1.636 ( 0.000 0.000 -0.574) 0.574 1.636 ( 0.000 0.000 -0.574) 0.574 3.795 ( 0.000 0.000 1.745) 1.745 3.834 ( -0.000 0.000 -1.703) 1.703 6.102 ( 0.000 0.000 0.223) 0.223 6.102 ( 0.000 0.000 0.223) 0.223 6.107 ( -0.000 0.000 -0.227) 0.227 6.107 ( 0.000 0.000 -0.227) 0.227 6.271 ( -0.000 0.000 0.033) 0.033 6.271 ( -0.000 0.000 0.033) 0.033 6.272 ( -0.000 0.000 -0.029) 0.029 6.272 ( 0.000 0.000 -0.029) 0.029 7.239 ( 0.000 0.000 1.327) 1.327 7.269 ( 0.000 -0.000 -1.319) 1.319 9.149 ( -0.000 -0.000 0.132) 0.132 9.152 ( -0.000 0.000 -0.127) 0.127 ======================= 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.383 ( 7.997 4.617 10.155) 13.726 0.510 ( 18.023 10.405 3.501) 21.103 0.579 ( 3.985 2.301 -7.935) 9.173 0.585 ( 15.913 9.188 -3.167) 18.646 0.895 ( 34.087 19.680 2.432) 39.435 0.944 ( 32.543 18.789 -1.966) 37.629 1.701 ( 6.307 3.641 0.610) 7.308 1.714 ( 6.279 3.625 -0.567) 7.272 1.907 ( 19.740 11.397 0.627) 22.802 1.913 ( 20.801 12.010 -0.001) 24.019 3.763 ( -2.363 -1.364 1.851) 3.297 3.805 ( -2.161 -1.247 -1.892) 3.131 6.090 ( -0.940 -0.543 0.315) 1.130 6.097 ( -0.827 -0.478 -0.335) 1.012 6.128 ( 1.308 0.755 0.624) 1.634 6.141 ( 1.677 0.968 -0.538) 2.010 6.260 ( -0.811 -0.468 0.264) 0.973 6.265 ( -0.640 -0.369 -0.244) 0.778 6.499 ( 15.072 8.702 6.330) 18.519 6.752 ( 17.094 9.869 -19.708) 27.893 7.196 ( -2.337 -1.349 -0.896) 2.844 7.208 ( -4.379 -2.528 -0.986) 5.152 9.115 ( -2.710 -1.565 0.114) 3.132 9.118 ( -2.742 -1.583 -0.117) 3.169 ======================= 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.669 ( 14.726 8.502 4.310) 17.542 0.767 ( 11.208 6.471 -4.471) 13.692 0.935 ( 17.220 9.942 1.944) 19.979 0.979 ( 16.371 9.452 -1.921) 19.000 1.648 ( 28.083 16.214 1.202) 32.449 1.674 ( 27.295 15.759 -1.112) 31.538 1.899 ( 9.779 5.646 0.535) 11.305 1.911 ( 9.675 5.586 -0.499) 11.183 2.444 ( 23.678 13.671 1.619) 27.389 2.473 ( 24.340 14.053 -1.035) 28.125 3.710 ( -1.376 -0.794 1.772) 2.380 3.750 ( -1.826 -1.054 -1.837) 2.796 6.061 ( -1.437 -0.830 0.353) 1.696 6.069 ( -1.481 -0.855 -0.383) 1.752 6.108 ( -3.978 -2.297 0.636) 4.637 6.122 ( -4.275 -2.468 -0.564) 4.968 6.235 ( -1.248 -0.721 0.348) 1.483 6.242 ( -1.271 -0.734 -0.318) 1.502 6.823 ( 10.901 6.294 5.075) 13.572 6.927 ( 2.129 1.229 -3.340) 4.147 7.123 ( -1.564 -0.903 0.834) 1.990 7.191 ( 1.391 0.803 -6.479) 6.675 9.031 ( -4.090 -2.362 0.033) 4.723 9.032 ( -4.202 -2.426 -0.044) 4.852 ======================= 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.039 ( 15.896 9.177 1.283) 18.399 1.073 ( 13.995 8.080 -1.724) 16.252 1.310 ( 13.996 8.081 1.244) 16.209 1.338 ( 13.518 7.805 -1.254) 15.660 2.138 ( 9.881 5.705 0.411) 11.417 2.147 ( 9.783 5.648 -0.387) 11.303 2.159 ( 13.809 7.972 1.531) 16.018 2.191 ( 15.328 8.849 -1.333) 17.749 2.963 ( 18.905 10.915 1.850) 21.907 3.000 ( 19.027 10.986 -1.534) 22.025 3.744 ( 5.258 3.036 0.716) 6.114 3.762 ( 3.889 2.245 -0.846) 4.570 5.956 ( -7.248 -4.185 0.309) 8.376 5.962 ( -7.351 -4.244 -0.261) 8.492 6.026 ( -1.444 -0.834 0.252) 1.687 6.032 ( -1.566 -0.904 -0.268) 1.828 6.204 ( -1.286 -0.742 0.258) 1.507 6.209 ( -1.377 -0.795 -0.242) 1.608 6.986 ( 3.083 1.780 1.056) 3.713 7.008 ( 5.053 2.917 -0.817) 5.892 7.137 ( 1.713 0.989 2.813) 3.439 7.210 ( -0.196 -0.113 -3.721) 3.728 8.934 ( -3.761 -2.171 0.114) 4.344 8.936 ( -3.592 -2.074 -0.128) 4.150 ======================= 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.395 ( 13.747 7.937 0.896) 15.899 1.418 ( 14.648 8.457 -1.194) 16.956 1.589 ( 9.355 5.401 0.898) 10.840 1.609 ( 9.181 5.301 -0.911) 10.640 2.307 ( -0.191 -0.110 3.091) 3.099 2.349 ( 7.569 4.370 0.363) 8.748 2.357 ( 7.597 4.386 -0.345) 8.779 2.375 ( 0.906 0.523 -2.974) 3.153 3.309 ( 10.138 5.853 1.789) 11.842 3.346 ( 10.052 5.803 -1.565) 11.712 3.954 ( 11.743 6.780 0.648) 13.575 3.969 ( 12.683 7.323 -0.710) 14.662 5.846 ( -0.814 -0.470 0.517) 1.073 5.857 ( -0.469 -0.271 -0.494) 0.733 5.996 ( -1.011 -0.584 0.219) 1.188 6.001 ( -0.934 -0.539 -0.231) 1.103 6.177 ( -0.848 -0.490 0.240) 1.008 6.183 ( -0.782 -0.451 -0.227) 0.931 7.048 ( 1.582 0.914 1.368) 2.282 7.081 ( 1.079 0.623 -1.625) 2.048 7.156 ( -0.413 -0.239 1.134) 1.231 7.179 ( -2.359 -1.362 -0.985) 2.896 8.864 ( -1.985 -1.146 0.301) 2.311 8.871 ( -1.816 -1.049 -0.306) 2.119 ======================= 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.692 ( 11.116 6.418 1.860) 12.969 1.736 ( 11.724 6.769 -2.043) 13.691 1.752 ( 4.506 2.602 0.837) 5.270 1.771 ( 4.528 2.614 -0.853) 5.297 2.220 ( -5.841 -3.372 3.582) 7.636 2.300 ( -5.848 -3.376 -3.559) 7.633 2.491 ( 4.392 2.536 0.447) 5.091 2.501 ( 4.487 2.590 -0.428) 5.198 3.468 ( 4.216 2.434 1.710) 5.160 3.504 ( 4.201 2.426 -1.516) 5.083 4.262 ( 13.140 7.587 1.327) 15.231 4.292 ( 13.460 7.771 -1.313) 15.597 5.923 ( 6.588 3.804 0.846) 7.654 5.942 ( 6.846 3.953 -0.835) 7.950 5.981 ( -0.229 -0.132 0.329) 0.422 5.989 ( -0.120 -0.069 -0.355) 0.381 6.167 ( -0.040 -0.023 0.340) 0.343 6.174 ( 0.048 0.028 -0.314) 0.319 7.046 ( -1.990 -1.149 0.217) 2.308 7.049 ( -3.406 -1.967 -0.096) 3.935 7.109 ( -3.168 -1.829 0.955) 3.781 7.132 ( -1.439 -0.831 -1.095) 1.990 8.838 ( -0.309 -0.178 0.447) 0.572 8.848 ( -0.229 -0.132 -0.443) 0.516 ======================= 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.815 ( 1.152 0.665 0.836) 1.571 1.834 ( 1.090 0.629 -0.852) 1.519 1.913 ( 7.302 4.216 2.127) 8.696 1.962 ( 7.205 4.160 -2.306) 8.633 2.077 ( -5.434 -3.138 3.196) 7.042 2.148 ( -6.130 -3.539 -3.142) 7.744 2.563 ( 1.872 1.081 0.476) 2.213 2.573 ( 1.801 1.040 -0.456) 2.129 3.545 ( 2.583 1.491 1.692) 3.429 3.581 ( 2.595 1.498 -1.548) 3.373 4.541 ( 9.819 5.669 1.424) 11.427 4.572 ( 9.655 5.574 -1.386) 11.234 5.983 ( 0.281 0.162 0.413) 0.526 5.992 ( 0.343 0.198 -0.454) 0.602 6.110 ( 8.126 4.691 0.982) 9.434 6.132 ( 8.093 4.672 -0.980) 9.396 6.173 ( 0.469 0.271 0.427) 0.690 6.182 ( 0.531 0.307 -0.387) 0.725 6.951 ( -4.171 -2.408 0.412) 4.834 6.966 ( -3.937 -2.273 -1.019) 4.658 7.039 ( -2.830 -1.634 2.301) 3.997 7.082 ( -2.645 -1.527 -1.698) 3.494 8.842 ( 0.412 0.238 0.452) 0.657 8.852 ( 0.345 0.199 -0.446) 0.598 ======================= 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.825 ( 0.048 0.028 0.714) 0.716 1.841 ( -0.074 -0.042 -0.722) 0.727 1.980 ( -1.923 -1.110 2.293) 3.192 2.033 ( -0.214 -0.124 -0.782) 0.821 2.038 ( 0.314 0.181 0.265) 0.450 2.080 ( 2.134 1.232 -1.880) 3.099 2.589 ( 0.510 0.294 0.303) 0.662 2.596 ( 0.318 0.183 -0.291) 0.468 3.591 ( 1.106 0.638 1.670) 2.102 3.627 ( 1.098 0.634 -1.567) 2.016 4.707 ( 3.789 2.188 1.069) 4.504 4.731 ( 3.420 1.974 -1.029) 4.081 5.989 ( 0.192 0.111 0.452) 0.503 6.000 ( 0.212 0.122 -0.501) 0.558 6.183 ( 0.277 0.160 0.480) 0.577 6.193 ( 0.302 0.174 -0.431) 0.554 6.264 ( 3.952 2.282 0.719) 4.619 6.280 ( 3.541 2.044 -0.704) 4.149 6.876 ( -1.802 -1.041 0.510) 2.143 6.892 ( -1.784 -1.030 -1.108) 2.339 6.982 ( -1.575 -0.910 2.229) 2.877 7.024 ( -1.721 -0.993 -1.654) 2.585 8.851 ( 0.260 0.150 0.352) 0.463 8.859 ( 0.177 0.102 -0.347) 0.403 ======================= 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.573 ( 7.655 13.259 5.526) 16.277 0.695 ( 5.374 9.308 -5.441) 12.047 0.874 ( 11.965 20.723 2.487) 24.058 0.932 ( 11.434 19.804 -2.702) 23.027 1.437 ( 16.655 28.847 0.828) 33.320 1.452 ( 15.862 27.473 -0.494) 31.727 1.859 ( 6.326 10.958 0.592) 12.667 1.871 ( 6.290 10.895 -0.492) 12.591 2.278 ( 12.951 22.431 1.358) 25.937 2.301 ( 13.429 23.260 -0.747) 26.869 3.720 ( -1.284 -2.224 1.880) 3.182 3.763 ( -1.380 -2.391 -1.935) 3.371 6.063 ( -1.542 -2.670 0.233) 3.092 6.067 ( -1.259 -2.180 -0.123) 2.520 6.144 ( -0.191 -0.331 0.257) 0.460 6.151 ( -0.509 -0.882 -0.340) 1.074 6.251 ( -0.427 -0.740 0.055) 0.857 6.251 ( -0.391 -0.678 -0.020) 0.783 6.745 ( 7.121 12.333 6.541) 15.672 6.907 ( 1.755 3.039 -6.795) 7.648 7.134 ( -1.945 -3.369 0.136) 3.892 7.181 ( 0.363 0.629 -5.937) 5.981 9.056 ( -2.285 -3.959 0.073) 4.572 9.058 ( -2.332 -4.039 -0.077) 4.664 ======================= 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.890 ( 12.206 12.431 2.051) 17.542 0.941 ( 10.287 10.135 -2.543) 14.664 1.282 ( 8.853 21.187 1.595) 23.018 1.319 ( 7.987 20.755 -1.671) 22.302 1.965 ( 16.626 12.756 0.952) 20.977 1.984 ( 17.770 13.621 -0.732) 22.401 2.113 ( 6.632 14.352 0.486) 15.817 2.123 ( 6.507 14.229 -0.389) 15.651 2.755 ( 17.107 15.304 1.729) 23.018 2.790 ( 17.392 15.459 -1.379) 23.310 3.706 ( 1.574 1.057 1.281) 2.288 3.736 ( 0.724 0.248 -1.380) 1.578 5.977 ( -3.181 -5.971 0.447) 6.780 5.987 ( -3.337 -6.009 -0.436) 6.887 6.097 ( -3.740 -0.100 0.240) 3.749 6.102 ( -3.603 -0.040 -0.205) 3.610 6.232 ( -1.328 0.369 0.181) 1.390 6.236 ( -1.163 0.193 -0.171) 1.192 6.946 ( 2.807 1.593 0.506) 3.267 6.946 ( 4.611 3.866 0.572) 6.045 7.109 ( 1.209 -0.175 2.794) 3.050 7.199 ( 0.625 -0.314 -5.150) 5.197 8.964 ( -3.145 -3.900 0.039) 5.010 8.965 ( -3.014 -3.792 -0.046) 4.844 ======================= 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.236 ( 13.081 9.577 0.425) 16.217 1.250 ( 13.304 9.140 -0.808) 16.161 1.607 ( 4.046 19.035 0.952) 19.484 1.628 ( 3.694 18.638 -0.946) 19.025 2.246 ( 6.972 0.638 2.403) 7.402 2.299 ( 8.223 1.577 -2.277) 8.677 2.362 ( 3.729 15.138 0.337) 15.594 2.369 ( 3.705 15.016 -0.256) 15.468 3.154 ( 13.843 7.334 1.778) 15.766 3.192 ( 13.871 7.247 -1.557) 15.728 3.832 ( 8.487 5.234 0.111) 9.972 3.835 ( 8.495 5.094 -0.209) 9.908 5.862 ( -2.626 -3.355 0.276) 4.270 5.868 ( -2.396 -3.393 -0.240) 4.161 6.043 ( -3.150 1.421 0.124) 3.458 6.045 ( -3.133 1.285 -0.139) 3.389 6.213 ( -1.849 1.576 0.199) 2.438 6.217 ( -1.792 1.556 -0.193) 2.381 7.000 ( 3.150 -0.633 1.519) 3.554 7.039 ( 3.116 -0.212 -2.128) 3.779 7.135 ( 1.282 -1.204 2.446) 3.013 7.184 ( -0.777 -2.482 -2.060) 3.318 8.879 ( -2.157 -2.977 0.216) 3.682 8.884 ( -1.991 -2.851 -0.223) 3.484 ======================= 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.545 ( 12.251 6.743 1.396) 14.054 1.579 ( 13.168 7.187 -1.634) 15.090 1.810 ( -0.489 14.074 0.750) 14.103 1.825 ( -0.494 13.930 -0.657) 13.955 2.264 ( -1.980 -2.729 3.140) 4.608 2.335 ( -1.407 -3.339 -3.195) 4.831 2.546 ( -0.006 13.929 0.341) 13.933 2.553 ( 0.139 13.899 -0.290) 13.903 3.392 ( 7.565 2.626 1.710) 8.188 3.428 ( 7.545 2.563 -1.519) 8.112 4.091 ( 13.221 6.187 0.955) 14.629 4.112 ( 13.910 6.271 -0.967) 15.289 5.860 ( 3.422 1.727 0.597) 3.879 5.873 ( 3.831 1.704 -0.568) 4.231 6.013 ( -1.942 2.073 0.296) 2.856 6.020 ( -1.713 2.300 -0.332) 2.887 6.202 ( -1.808 2.730 0.289) 3.287 6.208 ( -1.702 2.800 -0.269) 3.288 7.025 ( 1.485 -3.149 1.253) 3.700 7.054 ( 0.146 -3.135 -1.433) 3.450 7.114 ( -0.242 -3.801 0.472) 3.838 7.123 ( -2.129 -4.040 -0.333) 4.579 8.832 ( -0.347 -1.965 0.395) 2.034 8.841 ( -0.223 -1.866 -0.395) 1.920 ======================= 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.802 ( 9.683 4.658 1.958) 10.922 1.848 ( 9.406 4.637 -2.130) 10.701 1.885 ( -2.232 6.623 0.829) 7.037 1.900 ( -3.092 7.490 -0.521) 8.120 2.177 ( -6.204 1.973 2.530) 6.985 2.237 ( -6.060 -0.603 -2.829) 6.715 2.651 ( -3.129 12.222 0.428) 12.624 2.660 ( -3.039 12.156 -0.387) 12.536 3.508 ( 3.750 1.785 1.677) 4.479 3.544 ( 3.773 1.757 -1.524) 4.432 4.384 ( 12.811 4.183 1.299) 13.539 4.413 ( 13.055 3.942 -1.270) 13.696 5.966 ( 3.601 0.548 0.469) 3.672 5.976 ( 2.920 0.354 -0.398) 2.969 6.046 ( 3.512 5.633 0.660) 6.671 6.063 ( 4.065 5.845 -0.777) 7.162 6.206 ( -1.422 3.718 0.427) 4.004 6.215 ( -1.351 3.849 -0.377) 4.097 6.980 ( -3.019 -4.625 0.041) 5.524 6.982 ( -1.664 -4.870 -0.144) 5.149 7.037 ( -1.717 -5.129 1.321) 5.568 7.066 ( -0.063 -5.535 -1.227) 5.670 8.822 ( 0.968 -1.443 0.474) 1.802 8.832 ( 0.962 -1.435 -0.469) 1.790 ======================= 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.866 ( -1.907 0.356 1.052) 2.207 1.885 ( -3.010 2.420 -0.692) 3.924 1.980 ( 5.399 2.692 1.948) 6.340 2.024 ( 5.261 2.263 -2.008) 6.069 2.107 ( -7.472 8.374 1.191) 11.286 2.139 ( -7.176 5.262 -1.621) 9.045 2.695 ( -5.029 10.991 0.362) 12.093 2.703 ( -5.128 10.810 -0.327) 11.969 3.576 ( 1.924 1.549 1.645) 2.968 3.611 ( 1.945 1.468 -1.540) 2.883 4.606 ( 8.697 0.737 1.154) 8.804 4.632 ( 8.510 0.297 -1.118) 8.588 5.992 ( 0.270 0.631 0.387) 0.788 6.001 ( 0.418 0.562 -0.412) 0.812 6.205 ( 5.599 4.266 0.717) 7.076 6.217 ( 1.441 4.026 -0.206) 4.281 6.231 ( 2.635 4.599 -0.140) 5.302 6.235 ( -1.117 4.703 -0.343) 4.846 6.882 ( -2.073 -4.311 0.376) 4.798 6.895 ( -1.717 -4.574 -0.943) 4.975 6.961 ( -0.425 -5.822 2.064) 6.192 6.999 ( -0.098 -6.387 -1.506) 6.563 8.829 ( 1.331 -1.454 0.408) 2.013 8.838 ( 1.225 -1.508 -0.403) 1.984 ======================= 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.847 ( 0.169 -0.293 0.992) 1.048 1.866 ( -0.960 1.662 -0.668) 2.032 2.044 ( -0.082 0.142 1.628) 1.636 2.081 ( -0.544 0.943 -1.771) 2.079 2.084 ( -6.279 10.875 0.617) 12.572 2.102 ( -4.413 7.643 -0.912) 8.873 2.706 ( -5.959 10.320 0.211) 11.919 2.710 ( -5.917 10.249 -0.182) 11.836 3.602 ( -0.255 0.442 1.617) 1.695 3.637 ( -0.201 0.347 -1.535) 1.587 4.691 ( 2.035 -3.525 0.869) 4.162 4.710 ( 2.135 -3.698 -0.833) 4.350 5.998 ( -0.318 0.550 0.422) 0.763 6.008 ( -0.267 0.463 -0.460) 0.705 6.232 ( -2.443 4.232 0.530) 4.915 6.244 ( -2.507 4.343 -0.571) 5.047 6.293 ( -0.314 0.544 0.654) 0.907 6.306 ( -0.315 0.545 -0.552) 0.838 6.838 ( 1.464 -2.536 0.416) 2.958 6.852 ( 1.628 -2.820 -1.006) 3.408 6.922 ( 2.645 -4.581 1.978) 5.648 6.958 ( 2.925 -5.067 -1.421) 6.021 8.834 ( 1.000 -1.732 0.333) 2.028 8.841 ( 0.995 -1.723 -0.330) 2.017 ======================= 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.163 ( 8.014 13.880 0.186) 16.028 1.172 ( 6.987 12.102 -0.604) 13.987 1.713 ( 10.963 18.988 1.069) 21.952 1.737 ( 10.616 18.388 -1.125) 21.262 2.150 ( 3.059 5.298 2.156) 6.486 2.196 ( 3.754 6.503 -1.941) 7.755 2.435 ( 8.993 15.576 0.330) 17.988 2.441 ( 8.877 15.376 -0.235) 17.757 3.044 ( 7.251 12.560 1.715) 14.604 3.080 ( 7.255 12.566 -1.517) 14.589 3.778 ( 3.676 6.367 0.344) 7.360 3.787 ( 3.032 5.252 -0.443) 6.080 5.867 ( -2.272 -3.934 0.218) 4.548 5.871 ( -2.476 -4.288 -0.187) 4.955 6.089 ( -0.379 -0.657 0.122) 0.768 6.092 ( -0.527 -0.913 -0.121) 1.062 6.244 ( 0.422 0.730 0.138) 0.854 6.247 ( 0.453 0.785 -0.140) 0.917 6.960 ( 0.329 0.569 1.577) 1.709 7.002 ( 0.566 0.981 -2.353) 2.611 7.112 ( 0.003 0.006 3.110) 3.110 7.174 ( -1.234 -2.137 -2.628) 3.605 8.883 ( -2.184 -3.782 0.175) 4.371 8.887 ( -2.079 -3.601 -0.179) 4.162 ======================= 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.429 ( 9.218 9.498 0.876) 13.265 1.452 ( 10.143 10.263 -1.131) 14.474 1.994 ( -1.157 14.604 1.272) 14.705 2.019 ( -0.940 16.740 -0.946) 16.793 2.224 ( 4.582 0.219 2.112) 5.050 2.274 ( 4.448 -2.042 -2.368) 5.437 2.710 ( 2.329 17.291 0.255) 17.449 2.715 ( 2.523 17.239 -0.180) 17.424 3.271 ( 8.414 4.543 1.663) 9.706 3.306 ( 8.403 4.480 -1.513) 9.643 3.956 ( 8.449 6.969 0.460) 10.962 3.967 ( 9.067 7.298 -0.491) 11.650 5.835 ( 0.830 0.540 0.233) 1.017 5.840 ( 1.275 0.503 -0.193) 1.384 6.072 ( -1.630 1.432 0.188) 2.178 6.077 ( -1.330 1.743 -0.227) 2.204 6.259 ( -0.677 2.669 0.228) 2.762 6.264 ( -0.566 2.738 -0.214) 2.804 6.961 ( 0.916 -2.989 1.518) 3.475 7.000 ( 0.720 -3.269 -2.118) 3.961 7.088 ( -1.072 -3.414 1.453) 3.862 7.112 ( -1.905 -4.327 -0.914) 4.815 8.818 ( -1.101 -2.788 0.351) 3.018 8.826 ( -0.980 -2.637 -0.352) 2.835 ======================= 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.685 ( 9.532 7.282 1.664) 12.111 1.725 ( 10.006 7.335 -1.850) 12.544 1.987 ( -4.404 -0.550 2.429) 5.060 2.038 ( -5.137 1.624 -2.077) 5.774 2.330 ( -1.795 11.950 0.800) 12.111 2.351 ( -1.287 9.225 -1.116) 9.381 2.873 ( -3.337 16.736 0.411) 17.070 2.882 ( -3.264 16.842 -0.346) 17.158 3.438 ( 6.438 2.278 1.624) 7.020 3.473 ( 6.448 2.214 -1.500) 6.980 4.203 ( 11.592 4.859 0.902) 12.601 4.223 ( 11.986 4.661 -0.888) 12.891 5.897 ( 4.071 1.397 0.493) 4.332 5.908 ( 4.284 1.244 -0.475) 4.486 6.080 ( -0.189 4.286 0.562) 4.326 6.093 ( 0.073 4.565 -0.621) 4.608 6.277 ( -1.217 4.250 0.404) 4.440 6.286 ( -1.066 4.397 -0.357) 4.538 6.926 ( 0.360 -6.098 1.229) 6.231 6.957 ( -0.360 -6.084 -1.669) 6.319 7.010 ( -1.695 -5.943 0.708) 6.220 7.020 ( -2.564 -5.734 -0.285) 6.288 8.787 ( 0.438 -2.247 0.472) 2.337 8.798 ( 0.479 -2.192 -0.469) 2.292 ======================= 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.887 ( 0.081 -0.774 1.761) 1.925 1.911 ( 4.044 3.043 1.850) 5.388 1.935 ( 1.396 0.482 -2.009) 2.493 1.956 ( 2.382 1.597 -1.586) 3.277 2.380 ( -7.740 14.768 0.325) 16.677 2.389 ( -7.676 14.044 -0.470) 16.012 2.948 ( -6.575 15.624 0.400) 16.955 2.957 ( -6.699 15.580 -0.333) 16.963 3.548 ( 3.496 2.122 1.558) 4.376 3.582 ( 3.505 1.950 -1.477) 4.274 4.438 ( 10.499 1.345 0.957) 10.628 4.460 ( 10.613 0.888 -0.932) 10.691 5.969 ( 2.864 -0.264 0.355) 2.898 5.977 ( 2.690 -0.333 -0.343) 2.732 6.168 ( 3.730 5.789 0.741) 6.927 6.186 ( 3.846 5.803 -0.872) 7.016 6.306 ( -1.237 5.565 0.631) 5.736 6.319 ( -1.148 5.779 -0.512) 5.915 6.852 ( -0.490 -7.580 0.618) 7.621 6.865 ( -2.089 -6.523 -0.509) 6.869 6.912 ( -2.138 -6.581 0.286) 6.925 6.919 ( -0.415 -8.059 -0.402) 8.080 8.782 ( 1.434 -2.250 0.456) 2.706 8.793 ( 1.360 -2.299 -0.454) 2.709 ======================= 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.835 ( -0.437 -2.117 1.649) 2.719 1.872 ( -0.845 -2.307 -1.685) 2.979 2.037 ( 2.170 3.100 1.477) 4.062 2.070 ( 2.146 2.504 -1.482) 3.615 2.377 ( -9.352 15.382 0.346) 18.005 2.386 ( -9.623 16.059 -0.416) 18.726 2.972 ( -8.085 14.859 0.237) 16.918 2.976 ( -8.230 14.774 -0.170) 16.912 3.606 ( 0.937 1.314 1.459) 2.176 3.638 ( 1.015 1.046 -1.416) 2.032 4.582 ( 5.901 -2.682 0.734) 6.524 4.598 ( 5.844 -3.091 -0.708) 6.649 5.998 ( 0.935 -0.394 0.319) 1.063 6.005 ( 1.032 -0.489 -0.331) 1.189 6.284 ( 2.506 3.360 0.594) 4.234 6.299 ( 2.004 3.415 -0.831) 4.046 6.341 ( -1.917 6.031 0.812) 6.380 6.355 ( -2.226 6.429 -0.549) 6.826 6.776 ( 1.281 -7.022 0.149) 7.139 6.779 ( 0.272 -5.537 -0.148) 5.545 6.825 ( 0.625 -6.429 0.617) 6.489 6.839 ( 1.701 -8.360 -0.639) 8.556 8.786 ( 1.637 -2.515 0.355) 3.022 8.794 ( 1.567 -2.576 -0.353) 3.036 ======================= 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.048 10.475 1.461) 12.183 1.669 ( 6.251 10.828 -1.634) 12.609 2.023 ( -3.174 -5.498 3.072) 7.053 2.091 ( -3.332 -5.771 -3.037) 7.324 2.421 ( 7.228 12.520 0.121) 14.457 2.423 ( 7.202 12.475 -0.115) 14.405 3.022 ( 6.929 12.001 0.220) 13.859 3.026 ( 7.030 12.175 -0.123) 14.060 3.368 ( 3.141 5.440 1.583) 6.478 3.403 ( 3.116 5.398 -1.506) 6.412 4.116 ( 4.940 8.556 0.658) 9.901 4.131 ( 4.984 8.633 -0.651) 9.990 5.867 ( 1.229 2.129 0.297) 2.477 5.873 ( 1.326 2.297 -0.277) 2.666 6.109 ( 1.382 2.394 0.511) 2.811 6.121 ( 1.577 2.731 -0.567) 3.204 6.314 ( 1.451 2.513 0.353) 2.923 6.321 ( 1.554 2.692 -0.312) 3.124 6.885 ( -2.461 -4.262 1.291) 5.088 6.918 ( -2.604 -4.510 -1.738) 5.490 6.998 ( -2.982 -5.166 0.887) 6.030 7.012 ( -3.066 -5.311 -0.461) 6.150 8.770 ( -1.048 -1.815 0.466) 2.147 8.781 ( -1.002 -1.735 -0.465) 2.056 ======================= 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.849 ( 6.238 8.098 1.444) 10.323 1.886 ( 6.303 7.961 -1.772) 10.308 1.908 ( -2.498 -4.096 2.581) 5.448 1.963 ( -3.347 -5.168 -2.348) 6.590 2.610 ( -0.533 12.355 0.347) 12.372 2.618 ( -0.336 12.922 -0.358) 12.931 3.192 ( -1.946 13.341 0.485) 13.490 3.201 ( -1.860 13.415 -0.384) 13.549 3.492 ( 5.201 2.925 1.483) 6.149 3.525 ( 5.093 2.798 -1.450) 5.989 4.298 ( 7.440 4.604 0.656) 8.773 4.313 ( 7.527 4.292 -0.641) 8.689 5.913 ( 2.412 0.306 0.360) 2.457 5.921 ( 2.430 0.194 -0.360) 2.464 6.177 ( 2.111 4.899 0.687) 5.379 6.194 ( 2.195 4.932 -0.780) 5.454 6.364 ( 0.561 4.015 0.587) 4.096 6.376 ( 0.739 4.210 -0.494) 4.303 6.790 ( -1.748 -6.798 1.106) 7.106 6.819 ( -1.845 -6.987 -1.587) 7.398 6.883 ( -3.138 -6.259 0.936) 7.064 6.897 ( -3.214 -6.040 -0.463) 6.858 8.744 ( 0.182 -1.864 0.477) 1.932 8.755 ( 0.142 -1.913 -0.476) 1.976 ======================= 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.824 ( -1.242 -2.846 1.690) 3.535 1.863 ( -1.698 -3.687 -1.779) 4.432 2.034 ( 3.874 6.647 1.271) 7.798 2.062 ( 3.701 5.980 -1.248) 7.143 2.673 ( -6.013 12.570 0.804) 13.957 2.691 ( -5.979 13.444 -0.833) 14.737 3.252 ( -5.951 13.091 0.404) 14.386 3.260 ( -6.195 13.187 -0.291) 14.573 3.590 ( 2.818 1.807 1.248) 3.573 3.618 ( 2.738 1.401 -1.257) 3.322 4.454 ( 6.395 0.452 0.530) 6.433 4.465 ( 6.441 0.024 -0.518) 6.462 5.954 ( 2.510 -1.146 0.273) 2.773 5.960 ( 2.456 -1.242 -0.273) 2.766 6.275 ( 2.902 4.520 0.692) 5.416 6.292 ( 2.937 4.445 -0.865) 5.397 6.419 ( 0.203 5.089 0.860) 5.165 6.436 ( 0.357 5.289 -0.686) 5.345 6.689 ( -0.828 -7.651 0.943) 7.753 6.714 ( -1.031 -7.733 -1.395) 7.925 6.767 ( -2.218 -6.353 0.893) 6.788 6.780 ( -2.527 -6.012 -0.445) 6.536 8.735 ( 1.066 -2.181 0.375) 2.456 8.744 ( 0.978 -2.311 -0.374) 2.537 ======================= 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.799 ( 0.705 -1.220 1.241) 1.878 1.828 ( 1.003 -1.737 -1.294) 2.387 2.106 ( -1.939 3.359 0.969) 3.997 2.128 ( -1.661 2.878 -0.971) 3.462 2.680 ( -7.420 12.852 1.036) 14.877 2.704 ( -7.854 13.604 -1.075) 15.746 3.264 ( -7.325 12.687 0.340) 14.654 3.270 ( -7.470 12.939 -0.218) 14.942 3.625 ( -0.231 0.401 1.079) 1.174 3.650 ( 0.032 -0.056 -1.109) 1.111 4.516 ( 1.896 -3.283 0.420) 3.815 4.525 ( 2.057 -3.562 -0.410) 4.133 5.971 ( 1.176 -2.036 0.230) 2.362 5.977 ( 1.232 -2.134 -0.233) 2.475 6.337 ( -0.941 1.629 0.610) 1.978 6.353 ( -0.950 1.645 -0.810) 2.065 6.455 ( -2.553 4.422 1.035) 5.210 6.475 ( -2.647 4.585 -0.813) 5.357 6.633 ( 3.349 -5.800 0.759) 6.741 6.654 ( 3.469 -6.009 -1.161) 7.035 6.704 ( 2.399 -4.155 0.757) 4.857 6.716 ( 2.320 -4.018 -0.375) 4.654 8.734 ( 1.354 -2.345 0.289) 2.724 8.741 ( 1.402 -2.428 -0.289) 2.818 ======================= 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.820 ( -1.687 -2.922 1.665) 3.763 1.859 ( -2.239 -3.879 -1.758) 4.812 2.034 ( 4.697 8.135 1.186) 9.468 2.060 ( 4.360 7.552 -1.153) 8.797 2.799 ( 3.424 5.930 1.023) 6.924 2.822 ( 3.892 6.741 -1.051) 7.854 3.393 ( 3.528 6.111 0.701) 7.091 3.408 ( 3.670 6.357 -0.614) 7.366 3.557 ( 1.970 3.412 1.331) 4.158 3.587 ( 1.793 3.105 -1.318) 3.820 4.394 ( 2.560 4.434 0.378) 5.134 4.402 ( 2.408 4.171 -0.372) 4.831 5.919 ( 0.261 0.453 0.268) 0.588 5.925 ( 0.205 0.355 -0.270) 0.490 6.273 ( 2.501 4.332 0.644) 5.043 6.289 ( 2.481 4.297 -0.787) 5.024 6.440 ( 2.014 3.489 0.798) 4.107 6.457 ( 2.126 3.682 -0.656) 4.302 6.657 ( -3.421 -5.926 0.945) 6.908 6.682 ( -3.570 -6.183 -1.316) 7.260 6.759 ( -3.261 -5.648 0.982) 6.595 6.776 ( -3.269 -5.663 -0.609) 6.567 8.713 ( -0.643 -1.114 0.378) 1.341 8.722 ( -0.736 -1.274 -0.378) 1.519 ======================= 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.784 ( -0.230 -1.007 0.871) 1.352 1.804 ( -0.489 -1.973 -0.896) 2.221 2.161 ( 1.137 5.080 0.644) 5.245 2.175 ( 0.984 4.394 -0.636) 4.547 2.868 ( -2.040 6.228 1.589) 6.744 2.905 ( -2.047 7.002 -1.674) 7.485 3.462 ( -2.603 6.999 0.732) 7.503 3.476 ( -2.993 7.489 -0.568) 8.085 3.618 ( 1.700 0.851 0.760) 2.047 3.636 ( 1.325 0.250 -0.819) 1.577 4.464 ( 1.969 0.586 0.197) 2.064 4.469 ( 2.002 0.307 -0.196) 2.035 5.930 ( 1.076 -0.943 0.169) 1.441 5.934 ( 1.036 -1.088 -0.165) 1.511 6.350 ( 1.608 2.752 0.633) 3.250 6.366 ( 1.620 2.773 -0.817) 3.314 6.506 ( 0.983 3.260 0.988) 3.545 6.526 ( 1.159 3.315 -0.835) 3.610 6.558 ( -0.803 -4.669 0.664) 4.784 6.574 ( -1.144 -5.379 -0.780) 5.554 6.660 ( -1.649 -3.918 0.620) 4.296 6.672 ( -2.035 -4.318 -0.471) 4.796 8.699 ( 0.361 -1.274 0.215) 1.342 8.703 ( 0.307 -1.487 -0.217) 1.533 ======================= 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.779 ( -0.240 -0.416 0.000) 0.481 1.779 ( 0.240 0.416 0.000) 0.481 2.220 ( -0.177 -0.306 0.007) 0.353 2.220 ( 0.177 0.306 0.007) 0.353 2.931 ( -0.000 -0.000 1.954) 1.954 2.976 ( -0.000 -0.000 -2.092) 2.092 3.535 ( -0.000 -0.000 1.045) 1.045 3.557 ( -0.000 -0.000 -0.897) 0.897 3.630 ( -0.287 -0.497 -0.010) 0.574 3.630 ( 0.287 0.497 -0.010) 0.574 4.473 ( -0.044 -0.077 -0.001) 0.089 4.473 ( 0.044 0.077 -0.001) 0.089 5.921 ( -0.051 -0.089 0.005) 0.103 5.921 ( 0.051 0.089 0.005) 0.103 6.382 ( -0.000 -0.000 0.720) 0.720 6.400 ( -0.000 -0.000 -0.970) 0.970 6.506 ( -0.199 -0.345 0.001) 0.398 6.506 ( 0.199 0.345 0.001) 0.398 6.548 ( -0.000 -0.000 1.236) 1.236 6.572 ( -0.000 -0.000 -0.985) 0.985 6.615 ( -0.136 -0.236 0.000) 0.272 6.615 ( 0.136 0.236 0.000) 0.272 8.686 ( -0.062 -0.107 -0.001) 0.124 8.686 ( 0.062 0.107 -0.001) 0.124 =================== 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 311.984 311.984 18.670 0.000 -0.000 -0.000 3/16200 20.0 289.757 289.757 10.399 0.000 0.000 -0.000 3/16200 30.0 235.485 235.485 6.768 0.000 0.000 -0.000 3/16200 40.0 194.345 194.345 4.988 0.000 0.000 -0.000 3/16200 50.0 161.564 161.564 3.933 0.000 0.000 -0.000 3/16200 60.0 135.517 135.517 3.235 0.000 0.000 -0.000 3/16200 70.0 115.228 115.228 2.743 0.000 0.000 -0.000 3/16200 80.0 99.527 99.527 2.382 0.000 0.000 -0.000 3/16200 90.0 87.285 87.285 2.106 0.000 0.000 -0.000 3/16200 100.0 77.601 77.601 1.888 0.000 0.000 -0.000 3/16200 110.0 69.809 69.809 1.713 0.000 0.000 -0.000 3/16200 120.0 63.432 63.432 1.569 0.000 0.000 -0.000 3/16200 130.0 58.131 58.131 1.447 0.000 0.000 -0.000 3/16200 140.0 53.660 53.660 1.344 0.000 0.000 -0.000 3/16200 150.0 49.841 49.841 1.255 0.000 0.000 -0.000 3/16200 160.0 46.543 46.543 1.177 0.000 0.000 -0.000 3/16200 170.0 43.665 43.665 1.108 0.000 0.000 -0.000 3/16200 180.0 41.132 41.132 1.048 0.000 0.000 -0.000 3/16200 190.0 38.885 38.885 0.993 0.000 0.000 -0.000 3/16200 200.0 36.878 36.878 0.945 0.000 0.000 -0.000 3/16200 210.0 35.073 35.073 0.900 0.000 0.000 -0.000 3/16200 220.0 33.442 33.442 0.860 0.000 0.000 -0.000 3/16200 230.0 31.960 31.960 0.824 0.000 0.000 -0.000 3/16200 240.0 30.606 30.606 0.790 0.000 0.000 -0.000 3/16200 250.0 29.366 29.366 0.759 0.000 0.000 -0.000 3/16200 260.0 28.225 28.225 0.730 0.000 0.000 -0.000 3/16200 270.0 27.171 27.171 0.704 0.000 0.000 -0.000 3/16200 280.0 26.195 26.195 0.679 0.000 0.000 -0.000 3/16200 290.0 25.288 25.288 0.656 0.000 0.000 -0.000 3/16200 300.0 24.443 24.443 0.635 0.000 0.000 -0.000 3/16200 310.0 23.653 23.653 0.615 0.000 0.000 -0.000 3/16200 320.0 22.914 22.914 0.596 0.000 0.000 -0.000 3/16200 330.0 22.221 22.221 0.578 0.000 0.000 -0.000 3/16200 340.0 21.569 21.569 0.562 0.000 0.000 -0.000 3/16200 350.0 20.955 20.955 0.546 0.000 0.000 -0.000 3/16200 360.0 20.376 20.376 0.531 0.000 0.000 -0.000 3/16200 370.0 19.828 19.828 0.517 0.000 0.000 -0.000 3/16200 380.0 19.309 19.309 0.504 0.000 0.000 -0.000 3/16200 390.0 18.817 18.817 0.491 0.000 0.000 -0.000 3/16200 400.0 18.350 18.350 0.479 0.000 0.000 -0.000 3/16200 410.0 17.906 17.906 0.468 0.000 0.000 -0.000 3/16200 420.0 17.483 17.483 0.457 0.000 0.000 -0.000 3/16200 430.0 17.080 17.080 0.446 0.000 0.000 -0.000 3/16200 440.0 16.695 16.695 0.436 0.000 0.000 -0.000 3/16200 450.0 16.327 16.327 0.427 0.000 0.000 -0.000 3/16200 460.0 15.976 15.976 0.418 0.000 0.000 -0.000 3/16200 470.0 15.639 15.639 0.409 0.000 0.000 -0.000 3/16200 480.0 15.317 15.317 0.400 0.000 0.000 -0.000 3/16200 490.0 15.007 15.007 0.392 0.000 0.000 -0.000 3/16200 500.0 14.710 14.710 0.385 0.000 0.000 -0.000 3/16200 510.0 14.425 14.425 0.377 0.000 0.000 -0.000 3/16200 520.0 14.150 14.150 0.370 0.000 0.000 -0.000 3/16200 530.0 13.886 13.886 0.363 0.000 0.000 -0.000 3/16200 540.0 13.632 13.632 0.357 0.000 0.000 -0.000 3/16200 550.0 13.387 13.387 0.350 0.000 0.000 -0.000 3/16200 560.0 13.150 13.150 0.344 0.000 0.000 -0.000 3/16200 570.0 12.922 12.922 0.338 0.000 0.000 -0.000 3/16200 580.0 12.702 12.702 0.332 0.000 0.000 -0.000 3/16200 590.0 12.489 12.489 0.327 0.000 0.000 -0.000 3/16200 600.0 12.283 12.283 0.322 0.000 0.000 -0.000 3/16200 610.0 12.084 12.084 0.316 0.000 0.000 -0.000 3/16200 620.0 11.891 11.891 0.311 0.000 0.000 -0.000 3/16200 630.0 11.705 11.705 0.306 0.000 0.000 -0.000 3/16200 640.0 11.524 11.524 0.302 0.000 0.000 -0.000 3/16200 650.0 11.349 11.349 0.297 0.000 0.000 -0.000 3/16200 660.0 11.179 11.179 0.293 0.000 0.000 -0.000 3/16200 670.0 11.014 11.014 0.288 0.000 0.000 -0.000 3/16200 680.0 10.854 10.854 0.284 0.000 0.000 -0.000 3/16200 690.0 10.698 10.698 0.280 0.000 0.000 -0.000 3/16200 700.0 10.547 10.547 0.276 0.000 0.000 -0.000 3/16200 710.0 10.400 10.400 0.272 0.000 0.000 -0.000 3/16200 720.0 10.257 10.257 0.269 0.000 0.000 -0.000 3/16200 730.0 10.119 10.119 0.265 0.000 0.000 -0.000 3/16200 740.0 9.983 9.983 0.262 0.000 0.000 -0.000 3/16200 750.0 9.852 9.852 0.258 0.000 0.000 -0.000 3/16200 760.0 9.724 9.724 0.255 0.000 0.000 -0.000 3/16200 770.0 9.599 9.599 0.251 0.000 0.000 -0.000 3/16200 780.0 9.477 9.477 0.248 0.000 0.000 -0.000 3/16200 790.0 9.358 9.358 0.245 0.000 0.000 -0.000 3/16200 800.0 9.243 9.243 0.242 0.000 0.000 -0.000 3/16200 810.0 9.130 9.130 0.239 0.000 0.000 -0.000 3/16200 820.0 9.020 9.020 0.236 0.000 0.000 -0.000 3/16200 830.0 8.912 8.912 0.233 0.000 0.000 -0.000 3/16200 840.0 8.807 8.807 0.231 0.000 0.000 -0.000 3/16200 850.0 8.705 8.705 0.228 0.000 0.000 -0.000 3/16200 860.0 8.605 8.605 0.225 0.000 0.000 -0.000 3/16200 870.0 8.507 8.507 0.223 0.000 0.000 -0.000 3/16200 880.0 8.411 8.411 0.220 0.000 0.000 -0.000 3/16200 890.0 8.318 8.318 0.218 0.000 0.000 -0.000 3/16200 900.0 8.226 8.226 0.216 0.000 0.000 -0.000 3/16200 910.0 8.137 8.137 0.213 0.000 0.000 -0.000 3/16200 920.0 8.049 8.049 0.211 0.000 0.000 -0.000 3/16200 930.0 7.964 7.964 0.209 0.000 0.000 -0.000 3/16200 940.0 7.880 7.880 0.206 0.000 0.000 -0.000 3/16200 950.0 7.798 7.798 0.204 0.000 0.000 -0.000 3/16200 960.0 7.718 7.718 0.202 0.000 0.000 -0.000 3/16200 970.0 7.639 7.639 0.200 0.000 0.000 -0.000 3/16200 980.0 7.562 7.562 0.198 0.000 0.000 -0.000 3/16200 990.0 7.486 7.486 0.196 0.000 0.000 -0.000 3/16200 1000.0 7.412 7.412 0.194 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 19:01:04]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|