----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 05:10:34]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [3 3 2] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4/mmm (139) Number of symmetry operations in supercell: 576 ------------------------------ primitive cell ------------------------------ Lattice vectors: a -2.387467235000000 2.387467235000000 5.871517955000000 b 2.387467235000000 -2.387467235000000 5.871517955000000 c 2.387467235000000 2.387467235000000 -5.871517955000000 Atomic positions (fractional): *1 I 0.65165977896405 0.65165977896405 0.00000000000000 126.904 2 I 0.34834022103595 0.34834022103595 0.00000000000000 126.904 *3 Hg 0.88802910916933 0.88802910916933 0.00000000000000 200.590 4 Hg 0.11197089083067 0.11197089083067 0.00000000000000 200.590 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.774934470000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.774934470000000 0.000000000000000 c 0.000000000000000 0.000000000000000 11.743035910000000 Atomic positions (fractional): *1 I 0.50000000000000 0.50000000000000 0.15165977896405 126.904 > 1 2 I 0.00000000000000 0.00000000000000 0.34834022103595 126.904 > 2 3 I 0.00000000000000 0.00000000000000 0.65165977896405 126.904 > 1 4 I 0.50000000000000 0.50000000000000 0.84834022103595 126.904 > 2 *5 Hg 0.50000000000000 0.50000000000000 0.38802910916933 200.590 > 3 6 Hg 0.00000000000000 0.00000000000000 0.11197089083067 200.590 > 4 7 Hg 0.00000000000000 0.00000000000000 0.88802910916933 200.590 > 3 8 Hg 0.50000000000000 0.50000000000000 0.61197089083067 200.590 > 4 -------------------------------- super cell -------------------------------- Lattice vectors: a 14.324803409999999 0.000000000000000 0.000000000000000 b 0.000000000000000 14.324803409999999 0.000000000000000 c 0.000000000000000 0.000000000000000 23.486071819999999 Atomic positions (fractional): *1 I 0.16666666666667 0.16666666666667 0.07582988948203 126.904 > 1 2 I 0.50000000000000 0.16666666666667 0.07582988948203 126.904 > 1 3 I 0.83333333333333 0.16666666666667 0.07582988948203 126.904 > 1 4 I 0.16666666666667 0.50000000000000 0.07582988948203 126.904 > 1 5 I 0.50000000000000 0.50000000000000 0.07582988948203 126.904 > 1 6 I 0.83333333333333 0.50000000000000 0.07582988948203 126.904 > 1 7 I 0.16666666666667 0.83333333333333 0.07582988948203 126.904 > 1 8 I 0.50000000000000 0.83333333333333 0.07582988948203 126.904 > 1 9 I 0.83333333333333 0.83333333333333 0.07582988948203 126.904 > 1 10 I 0.16666666666667 0.16666666666667 0.57582988948202 126.904 > 1 11 I 0.50000000000000 0.16666666666667 0.57582988948202 126.904 > 1 12 I 0.83333333333333 0.16666666666667 0.57582988948202 126.904 > 1 13 I 0.16666666666667 0.50000000000000 0.57582988948202 126.904 > 1 14 I 0.50000000000000 0.50000000000000 0.57582988948202 126.904 > 1 15 I 0.83333333333333 0.50000000000000 0.57582988948202 126.904 > 1 16 I 0.16666666666667 0.83333333333333 0.57582988948202 126.904 > 1 17 I 0.50000000000000 0.83333333333333 0.57582988948202 126.904 > 1 18 I 0.83333333333333 0.83333333333333 0.57582988948202 126.904 > 1 19 I 0.00000000000000 0.00000000000000 0.17417011051798 126.904 > 2 20 I 0.33333333333333 0.00000000000000 0.17417011051798 126.904 > 2 21 I 0.66666666666667 0.00000000000000 0.17417011051798 126.904 > 2 22 I 0.00000000000000 0.33333333333333 0.17417011051798 126.904 > 2 23 I 0.33333333333333 0.33333333333333 0.17417011051798 126.904 > 2 24 I 0.66666666666667 0.33333333333333 0.17417011051798 126.904 > 2 25 I 0.00000000000000 0.66666666666667 0.17417011051798 126.904 > 2 26 I 0.33333333333333 0.66666666666667 0.17417011051798 126.904 > 2 27 I 0.66666666666667 0.66666666666667 0.17417011051798 126.904 > 2 28 I 0.00000000000000 0.00000000000000 0.67417011051798 126.904 > 2 29 I 0.33333333333333 0.00000000000000 0.67417011051798 126.904 > 2 30 I 0.66666666666667 0.00000000000000 0.67417011051798 126.904 > 2 31 I 0.00000000000000 0.33333333333333 0.67417011051798 126.904 > 2 32 I 0.33333333333333 0.33333333333333 0.67417011051798 126.904 > 2 33 I 0.66666666666667 0.33333333333333 0.67417011051798 126.904 > 2 34 I 0.00000000000000 0.66666666666667 0.67417011051798 126.904 > 2 35 I 0.33333333333333 0.66666666666667 0.67417011051798 126.904 > 2 36 I 0.66666666666667 0.66666666666667 0.67417011051798 126.904 > 2 37 I 0.00000000000000 0.00000000000000 0.32582988948203 126.904 > 1 38 I 0.33333333333333 0.00000000000000 0.32582988948203 126.904 > 1 39 I 0.66666666666667 0.00000000000000 0.32582988948203 126.904 > 1 40 I 0.00000000000000 0.33333333333333 0.32582988948203 126.904 > 1 41 I 0.33333333333333 0.33333333333333 0.32582988948203 126.904 > 1 42 I 0.66666666666667 0.33333333333333 0.32582988948203 126.904 > 1 43 I 0.00000000000000 0.66666666666667 0.32582988948203 126.904 > 1 44 I 0.33333333333333 0.66666666666667 0.32582988948203 126.904 > 1 45 I 0.66666666666667 0.66666666666667 0.32582988948203 126.904 > 1 46 I 0.00000000000000 0.00000000000000 0.82582988948203 126.904 > 1 47 I 0.33333333333333 0.00000000000000 0.82582988948203 126.904 > 1 48 I 0.66666666666667 0.00000000000000 0.82582988948203 126.904 > 1 49 I 0.00000000000000 0.33333333333333 0.82582988948203 126.904 > 1 50 I 0.33333333333333 0.33333333333333 0.82582988948203 126.904 > 1 51 I 0.66666666666667 0.33333333333333 0.82582988948203 126.904 > 1 52 I 0.00000000000000 0.66666666666667 0.82582988948203 126.904 > 1 53 I 0.33333333333333 0.66666666666667 0.82582988948203 126.904 > 1 54 I 0.66666666666667 0.66666666666667 0.82582988948203 126.904 > 1 55 I 0.16666666666667 0.16666666666667 0.42417011051797 126.904 > 2 56 I 0.50000000000000 0.16666666666667 0.42417011051797 126.904 > 2 57 I 0.83333333333333 0.16666666666667 0.42417011051797 126.904 > 2 58 I 0.16666666666667 0.50000000000000 0.42417011051797 126.904 > 2 59 I 0.50000000000000 0.50000000000000 0.42417011051797 126.904 > 2 60 I 0.83333333333333 0.50000000000000 0.42417011051797 126.904 > 2 61 I 0.16666666666667 0.83333333333333 0.42417011051797 126.904 > 2 62 I 0.50000000000000 0.83333333333333 0.42417011051797 126.904 > 2 63 I 0.83333333333333 0.83333333333333 0.42417011051797 126.904 > 2 64 I 0.16666666666667 0.16666666666667 0.92417011051797 126.904 > 2 65 I 0.50000000000000 0.16666666666667 0.92417011051797 126.904 > 2 66 I 0.83333333333333 0.16666666666667 0.92417011051797 126.904 > 2 67 I 0.16666666666667 0.50000000000000 0.92417011051797 126.904 > 2 68 I 0.50000000000000 0.50000000000000 0.92417011051797 126.904 > 2 69 I 0.83333333333333 0.50000000000000 0.92417011051797 126.904 > 2 70 I 0.16666666666667 0.83333333333333 0.92417011051797 126.904 > 2 71 I 0.50000000000000 0.83333333333333 0.92417011051797 126.904 > 2 72 I 0.83333333333333 0.83333333333333 0.92417011051797 126.904 > 2 *73 Hg 0.16666666666667 0.16666666666667 0.19401455458467 200.590 > 3 74 Hg 0.50000000000000 0.16666666666667 0.19401455458467 200.590 > 3 75 Hg 0.83333333333333 0.16666666666667 0.19401455458467 200.590 > 3 76 Hg 0.16666666666667 0.50000000000000 0.19401455458467 200.590 > 3 77 Hg 0.50000000000000 0.50000000000000 0.19401455458467 200.590 > 3 78 Hg 0.83333333333333 0.50000000000000 0.19401455458467 200.590 > 3 79 Hg 0.16666666666667 0.83333333333333 0.19401455458467 200.590 > 3 80 Hg 0.50000000000000 0.83333333333333 0.19401455458467 200.590 > 3 81 Hg 0.83333333333333 0.83333333333333 0.19401455458467 200.590 > 3 82 Hg 0.16666666666667 0.16666666666667 0.69401455458467 200.590 > 3 83 Hg 0.50000000000000 0.16666666666667 0.69401455458467 200.590 > 3 84 Hg 0.83333333333333 0.16666666666667 0.69401455458467 200.590 > 3 85 Hg 0.16666666666667 0.50000000000000 0.69401455458467 200.590 > 3 86 Hg 0.50000000000000 0.50000000000000 0.69401455458467 200.590 > 3 87 Hg 0.83333333333333 0.50000000000000 0.69401455458467 200.590 > 3 88 Hg 0.16666666666667 0.83333333333333 0.69401455458467 200.590 > 3 89 Hg 0.50000000000000 0.83333333333333 0.69401455458467 200.590 > 3 90 Hg 0.83333333333333 0.83333333333333 0.69401455458467 200.590 > 3 91 Hg 0.00000000000000 0.00000000000000 0.05598544541533 200.590 > 4 92 Hg 0.33333333333333 0.00000000000000 0.05598544541533 200.590 > 4 93 Hg 0.66666666666667 0.00000000000000 0.05598544541533 200.590 > 4 94 Hg 0.00000000000000 0.33333333333333 0.05598544541533 200.590 > 4 95 Hg 0.33333333333333 0.33333333333333 0.05598544541533 200.590 > 4 96 Hg 0.66666666666667 0.33333333333333 0.05598544541533 200.590 > 4 97 Hg 0.00000000000000 0.66666666666667 0.05598544541533 200.590 > 4 98 Hg 0.33333333333333 0.66666666666667 0.05598544541533 200.590 > 4 99 Hg 0.66666666666667 0.66666666666667 0.05598544541533 200.590 > 4 100 Hg 0.00000000000000 0.00000000000000 0.55598544541533 200.590 > 4 101 Hg 0.33333333333333 0.00000000000000 0.55598544541533 200.590 > 4 102 Hg 0.66666666666667 0.00000000000000 0.55598544541533 200.590 > 4 103 Hg 0.00000000000000 0.33333333333333 0.55598544541533 200.590 > 4 104 Hg 0.33333333333333 0.33333333333333 0.55598544541533 200.590 > 4 105 Hg 0.66666666666667 0.33333333333333 0.55598544541533 200.590 > 4 106 Hg 0.00000000000000 0.66666666666667 0.55598544541533 200.590 > 4 107 Hg 0.33333333333333 0.66666666666667 0.55598544541533 200.590 > 4 108 Hg 0.66666666666667 0.66666666666667 0.55598544541533 200.590 > 4 109 Hg 0.00000000000000 0.00000000000000 0.44401455458467 200.590 > 3 110 Hg 0.33333333333333 0.00000000000000 0.44401455458467 200.590 > 3 111 Hg 0.66666666666667 0.00000000000000 0.44401455458467 200.590 > 3 112 Hg 0.00000000000000 0.33333333333333 0.44401455458467 200.590 > 3 113 Hg 0.33333333333333 0.33333333333333 0.44401455458467 200.590 > 3 114 Hg 0.66666666666667 0.33333333333333 0.44401455458467 200.590 > 3 115 Hg 0.00000000000000 0.66666666666667 0.44401455458467 200.590 > 3 116 Hg 0.33333333333333 0.66666666666667 0.44401455458467 200.590 > 3 117 Hg 0.66666666666667 0.66666666666667 0.44401455458467 200.590 > 3 118 Hg 0.00000000000000 0.00000000000000 0.94401455458467 200.590 > 3 119 Hg 0.33333333333333 0.00000000000000 0.94401455458467 200.590 > 3 120 Hg 0.66666666666667 0.00000000000000 0.94401455458467 200.590 > 3 121 Hg 0.00000000000000 0.33333333333333 0.94401455458467 200.590 > 3 122 Hg 0.33333333333333 0.33333333333333 0.94401455458467 200.590 > 3 123 Hg 0.66666666666667 0.33333333333333 0.94401455458467 200.590 > 3 124 Hg 0.00000000000000 0.66666666666667 0.94401455458467 200.590 > 3 125 Hg 0.33333333333333 0.66666666666667 0.94401455458467 200.590 > 3 126 Hg 0.66666666666667 0.66666666666667 0.94401455458467 200.590 > 3 127 Hg 0.16666666666667 0.16666666666667 0.30598544541533 200.590 > 4 128 Hg 0.50000000000000 0.16666666666667 0.30598544541533 200.590 > 4 129 Hg 0.83333333333333 0.16666666666667 0.30598544541533 200.590 > 4 130 Hg 0.16666666666667 0.50000000000000 0.30598544541533 200.590 > 4 131 Hg 0.50000000000000 0.50000000000000 0.30598544541533 200.590 > 4 132 Hg 0.83333333333333 0.50000000000000 0.30598544541533 200.590 > 4 133 Hg 0.16666666666667 0.83333333333333 0.30598544541533 200.590 > 4 134 Hg 0.50000000000000 0.83333333333333 0.30598544541533 200.590 > 4 135 Hg 0.83333333333333 0.83333333333333 0.30598544541533 200.590 > 4 136 Hg 0.16666666666667 0.16666666666667 0.80598544541533 200.590 > 4 137 Hg 0.50000000000000 0.16666666666667 0.80598544541533 200.590 > 4 138 Hg 0.83333333333333 0.16666666666667 0.80598544541533 200.590 > 4 139 Hg 0.16666666666667 0.50000000000000 0.80598544541533 200.590 > 4 140 Hg 0.50000000000000 0.50000000000000 0.80598544541533 200.590 > 4 141 Hg 0.83333333333333 0.50000000000000 0.80598544541533 200.590 > 4 142 Hg 0.16666666666667 0.83333333333333 0.80598544541533 200.590 > 4 143 Hg 0.50000000000000 0.83333333333333 0.80598544541533 200.590 > 4 144 Hg 0.83333333333333 0.83333333333333 0.80598544541533 200.590 > 4 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.0906709 0.0000000 0.0000000 0.0000000 7.0906709 0.0000000 0.0000000 0.0000000 15.7187253 -------------------------- Born effective charges -------------------------- 1 I -1.7083236 0.0000000 0.0000000 0.0000000 -1.7083236 0.0000000 0.0000000 0.0000000 -2.9213064 2 I -1.7083236 0.0000000 0.0000000 0.0000000 -1.7083236 0.0000000 0.0000000 0.0000000 -2.9213064 3 Hg 1.7083236 0.0000000 0.0000000 0.0000000 1.7083236 0.0000000 0.0000000 0.0000000 2.9213064 4 Hg 1.7083236 0.0000000 0.0000000 0.0000000 1.7083236 0.0000000 0.0000000 0.0000000 2.9213064 ---------------------------------------------------------------------------- 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: 432/432 Permutation basis: 3540/3540 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 240 Number of blocks in projector: 240 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 4 --- Eigsh_solver_block: 1 / 4 --- Block_size: 88 Use standard eigh solver. --- Eigsh_solver_block: 2 / 4 --- Block_size: 66 Use standard eigh solver. --- Eigsh_solver_block: 3 / 4 --- Block_size: 64 Use standard eigh solver. --- Eigsh_solver_block: 4 / 4 --- Block_size: 22 Use standard eigh solver. Tree of FC basis block matrices: - (240, 236), data: False |-- (22, 22), data: True |-- (64, 64), data: True |-- (66, 64), data: True |-- (88, 86), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.007 Solver_block: 80 / 80 - Time: 0.204 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.213 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 432/432 Permutation basis: 3540/3540 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 240 Number of blocks in projector: 240 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 4 --- Eigsh_solver_block: 1 / 4 --- Block_size: 88 Use standard eigh solver. --- Eigsh_solver_block: 2 / 4 --- Block_size: 66 Use standard eigh solver. --- Eigsh_solver_block: 3 / 4 --- Block_size: 64 Use standard eigh solver. --- Eigsh_solver_block: 4 / 4 --- Block_size: 22 Use standard eigh solver. Tree of FC basis block matrices: - (240, 236), data: False |-- (22, 22), data: True |-- (64, 64), data: True |-- (66, 64), data: True |-- (88, 86), data: True Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 05:10:39]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 05:10:39]------------------------- 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): [3 3 2] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4/mmm (139) ------------------------------ primitive cell ------------------------------ Lattice vectors: a -2.387467235000000 2.387467235000000 5.871517955000000 b 2.387467235000000 -2.387467235000000 5.871517955000000 c 2.387467235000000 2.387467235000000 -5.871517955000000 Atomic positions (fractional): 1 I 0.65165977896405 0.65165977896405 0.00000000000000 126.904 2 I 0.34834022103595 0.34834022103595 0.00000000000000 126.904 3 Hg 0.88802910916933 0.88802910916933 0.00000000000000 200.590 4 Hg 0.11197089083067 0.11197089083067 0.00000000000000 200.590 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.324803409999999 0.000000000000000 0.000000000000000 b 0.000000000000000 14.324803409999999 0.000000000000000 c 0.000000000000000 0.000000000000000 23.486071819999999 Atomic positions (fractional): 1 I 0.16666666666667 0.16666666666667 0.07582988948203 126.904 > 1 2 I 0.50000000000000 0.16666666666667 0.07582988948203 126.904 > 1 3 I 0.83333333333333 0.16666666666667 0.07582988948203 126.904 > 1 4 I 0.16666666666667 0.50000000000000 0.07582988948203 126.904 > 1 5 I 0.50000000000000 0.50000000000000 0.07582988948203 126.904 > 1 6 I 0.83333333333333 0.50000000000000 0.07582988948203 126.904 > 1 7 I 0.16666666666667 0.83333333333333 0.07582988948203 126.904 > 1 8 I 0.50000000000000 0.83333333333333 0.07582988948203 126.904 > 1 9 I 0.83333333333333 0.83333333333333 0.07582988948203 126.904 > 1 10 I 0.16666666666667 0.16666666666667 0.57582988948202 126.904 > 1 11 I 0.50000000000000 0.16666666666667 0.57582988948202 126.904 > 1 12 I 0.83333333333333 0.16666666666667 0.57582988948202 126.904 > 1 13 I 0.16666666666667 0.50000000000000 0.57582988948202 126.904 > 1 14 I 0.50000000000000 0.50000000000000 0.57582988948202 126.904 > 1 15 I 0.83333333333333 0.50000000000000 0.57582988948202 126.904 > 1 16 I 0.16666666666667 0.83333333333333 0.57582988948202 126.904 > 1 17 I 0.50000000000000 0.83333333333333 0.57582988948202 126.904 > 1 18 I 0.83333333333333 0.83333333333333 0.57582988948202 126.904 > 1 19 I 0.00000000000000 0.00000000000000 0.17417011051798 126.904 > 19 20 I 0.33333333333333 0.00000000000000 0.17417011051798 126.904 > 19 21 I 0.66666666666667 0.00000000000000 0.17417011051798 126.904 > 19 22 I 0.00000000000000 0.33333333333333 0.17417011051798 126.904 > 19 23 I 0.33333333333333 0.33333333333333 0.17417011051798 126.904 > 19 24 I 0.66666666666667 0.33333333333333 0.17417011051798 126.904 > 19 25 I 0.00000000000000 0.66666666666667 0.17417011051798 126.904 > 19 26 I 0.33333333333333 0.66666666666667 0.17417011051798 126.904 > 19 27 I 0.66666666666667 0.66666666666667 0.17417011051798 126.904 > 19 28 I 0.00000000000000 0.00000000000000 0.67417011051798 126.904 > 19 29 I 0.33333333333333 0.00000000000000 0.67417011051798 126.904 > 19 30 I 0.66666666666667 0.00000000000000 0.67417011051798 126.904 > 19 31 I 0.00000000000000 0.33333333333333 0.67417011051798 126.904 > 19 32 I 0.33333333333333 0.33333333333333 0.67417011051798 126.904 > 19 33 I 0.66666666666667 0.33333333333333 0.67417011051798 126.904 > 19 34 I 0.00000000000000 0.66666666666667 0.67417011051798 126.904 > 19 35 I 0.33333333333333 0.66666666666667 0.67417011051798 126.904 > 19 36 I 0.66666666666667 0.66666666666667 0.67417011051798 126.904 > 19 37 I 0.00000000000000 0.00000000000000 0.32582988948203 126.904 > 1 38 I 0.33333333333333 0.00000000000000 0.32582988948203 126.904 > 1 39 I 0.66666666666667 0.00000000000000 0.32582988948203 126.904 > 1 40 I 0.00000000000000 0.33333333333333 0.32582988948203 126.904 > 1 41 I 0.33333333333333 0.33333333333333 0.32582988948203 126.904 > 1 42 I 0.66666666666667 0.33333333333333 0.32582988948203 126.904 > 1 43 I 0.00000000000000 0.66666666666667 0.32582988948203 126.904 > 1 44 I 0.33333333333333 0.66666666666667 0.32582988948203 126.904 > 1 45 I 0.66666666666667 0.66666666666667 0.32582988948203 126.904 > 1 46 I 0.00000000000000 0.00000000000000 0.82582988948203 126.904 > 1 47 I 0.33333333333333 0.00000000000000 0.82582988948203 126.904 > 1 48 I 0.66666666666667 0.00000000000000 0.82582988948203 126.904 > 1 49 I 0.00000000000000 0.33333333333333 0.82582988948203 126.904 > 1 50 I 0.33333333333333 0.33333333333333 0.82582988948203 126.904 > 1 51 I 0.66666666666667 0.33333333333333 0.82582988948203 126.904 > 1 52 I 0.00000000000000 0.66666666666667 0.82582988948203 126.904 > 1 53 I 0.33333333333333 0.66666666666667 0.82582988948203 126.904 > 1 54 I 0.66666666666667 0.66666666666667 0.82582988948203 126.904 > 1 55 I 0.16666666666667 0.16666666666667 0.42417011051797 126.904 > 19 56 I 0.50000000000000 0.16666666666667 0.42417011051797 126.904 > 19 57 I 0.83333333333333 0.16666666666667 0.42417011051797 126.904 > 19 58 I 0.16666666666667 0.50000000000000 0.42417011051797 126.904 > 19 59 I 0.50000000000000 0.50000000000000 0.42417011051797 126.904 > 19 60 I 0.83333333333333 0.50000000000000 0.42417011051797 126.904 > 19 61 I 0.16666666666667 0.83333333333333 0.42417011051797 126.904 > 19 62 I 0.50000000000000 0.83333333333333 0.42417011051797 126.904 > 19 63 I 0.83333333333333 0.83333333333333 0.42417011051797 126.904 > 19 64 I 0.16666666666667 0.16666666666667 0.92417011051797 126.904 > 19 65 I 0.50000000000000 0.16666666666667 0.92417011051797 126.904 > 19 66 I 0.83333333333333 0.16666666666667 0.92417011051797 126.904 > 19 67 I 0.16666666666667 0.50000000000000 0.92417011051797 126.904 > 19 68 I 0.50000000000000 0.50000000000000 0.92417011051797 126.904 > 19 69 I 0.83333333333333 0.50000000000000 0.92417011051797 126.904 > 19 70 I 0.16666666666667 0.83333333333333 0.92417011051797 126.904 > 19 71 I 0.50000000000000 0.83333333333333 0.92417011051797 126.904 > 19 72 I 0.83333333333333 0.83333333333333 0.92417011051797 126.904 > 19 73 Hg 0.16666666666667 0.16666666666667 0.19401455458467 200.590 > 73 74 Hg 0.50000000000000 0.16666666666667 0.19401455458467 200.590 > 73 75 Hg 0.83333333333333 0.16666666666667 0.19401455458467 200.590 > 73 76 Hg 0.16666666666667 0.50000000000000 0.19401455458467 200.590 > 73 77 Hg 0.50000000000000 0.50000000000000 0.19401455458467 200.590 > 73 78 Hg 0.83333333333333 0.50000000000000 0.19401455458467 200.590 > 73 79 Hg 0.16666666666667 0.83333333333333 0.19401455458467 200.590 > 73 80 Hg 0.50000000000000 0.83333333333333 0.19401455458467 200.590 > 73 81 Hg 0.83333333333333 0.83333333333333 0.19401455458467 200.590 > 73 82 Hg 0.16666666666667 0.16666666666667 0.69401455458467 200.590 > 73 83 Hg 0.50000000000000 0.16666666666667 0.69401455458467 200.590 > 73 84 Hg 0.83333333333333 0.16666666666667 0.69401455458467 200.590 > 73 85 Hg 0.16666666666667 0.50000000000000 0.69401455458467 200.590 > 73 86 Hg 0.50000000000000 0.50000000000000 0.69401455458467 200.590 > 73 87 Hg 0.83333333333333 0.50000000000000 0.69401455458467 200.590 > 73 88 Hg 0.16666666666667 0.83333333333333 0.69401455458467 200.590 > 73 89 Hg 0.50000000000000 0.83333333333333 0.69401455458467 200.590 > 73 90 Hg 0.83333333333333 0.83333333333333 0.69401455458467 200.590 > 73 91 Hg 0.00000000000000 0.00000000000000 0.05598544541533 200.590 > 91 92 Hg 0.33333333333333 0.00000000000000 0.05598544541533 200.590 > 91 93 Hg 0.66666666666667 0.00000000000000 0.05598544541533 200.590 > 91 94 Hg 0.00000000000000 0.33333333333333 0.05598544541533 200.590 > 91 95 Hg 0.33333333333333 0.33333333333333 0.05598544541533 200.590 > 91 96 Hg 0.66666666666667 0.33333333333333 0.05598544541533 200.590 > 91 97 Hg 0.00000000000000 0.66666666666667 0.05598544541533 200.590 > 91 98 Hg 0.33333333333333 0.66666666666667 0.05598544541533 200.590 > 91 99 Hg 0.66666666666667 0.66666666666667 0.05598544541533 200.590 > 91 100 Hg 0.00000000000000 0.00000000000000 0.55598544541533 200.590 > 91 101 Hg 0.33333333333333 0.00000000000000 0.55598544541533 200.590 > 91 102 Hg 0.66666666666667 0.00000000000000 0.55598544541533 200.590 > 91 103 Hg 0.00000000000000 0.33333333333333 0.55598544541533 200.590 > 91 104 Hg 0.33333333333333 0.33333333333333 0.55598544541533 200.590 > 91 105 Hg 0.66666666666667 0.33333333333333 0.55598544541533 200.590 > 91 106 Hg 0.00000000000000 0.66666666666667 0.55598544541533 200.590 > 91 107 Hg 0.33333333333333 0.66666666666667 0.55598544541533 200.590 > 91 108 Hg 0.66666666666667 0.66666666666667 0.55598544541533 200.590 > 91 109 Hg 0.00000000000000 0.00000000000000 0.44401455458467 200.590 > 73 110 Hg 0.33333333333333 0.00000000000000 0.44401455458467 200.590 > 73 111 Hg 0.66666666666667 0.00000000000000 0.44401455458467 200.590 > 73 112 Hg 0.00000000000000 0.33333333333333 0.44401455458467 200.590 > 73 113 Hg 0.33333333333333 0.33333333333333 0.44401455458467 200.590 > 73 114 Hg 0.66666666666667 0.33333333333333 0.44401455458467 200.590 > 73 115 Hg 0.00000000000000 0.66666666666667 0.44401455458467 200.590 > 73 116 Hg 0.33333333333333 0.66666666666667 0.44401455458467 200.590 > 73 117 Hg 0.66666666666667 0.66666666666667 0.44401455458467 200.590 > 73 118 Hg 0.00000000000000 0.00000000000000 0.94401455458467 200.590 > 73 119 Hg 0.33333333333333 0.00000000000000 0.94401455458467 200.590 > 73 120 Hg 0.66666666666667 0.00000000000000 0.94401455458467 200.590 > 73 121 Hg 0.00000000000000 0.33333333333333 0.94401455458467 200.590 > 73 122 Hg 0.33333333333333 0.33333333333333 0.94401455458467 200.590 > 73 123 Hg 0.66666666666667 0.33333333333333 0.94401455458467 200.590 > 73 124 Hg 0.00000000000000 0.66666666666667 0.94401455458467 200.590 > 73 125 Hg 0.33333333333333 0.66666666666667 0.94401455458467 200.590 > 73 126 Hg 0.66666666666667 0.66666666666667 0.94401455458467 200.590 > 73 127 Hg 0.16666666666667 0.16666666666667 0.30598544541533 200.590 > 91 128 Hg 0.50000000000000 0.16666666666667 0.30598544541533 200.590 > 91 129 Hg 0.83333333333333 0.16666666666667 0.30598544541533 200.590 > 91 130 Hg 0.16666666666667 0.50000000000000 0.30598544541533 200.590 > 91 131 Hg 0.50000000000000 0.50000000000000 0.30598544541533 200.590 > 91 132 Hg 0.83333333333333 0.50000000000000 0.30598544541533 200.590 > 91 133 Hg 0.16666666666667 0.83333333333333 0.30598544541533 200.590 > 91 134 Hg 0.50000000000000 0.83333333333333 0.30598544541533 200.590 > 91 135 Hg 0.83333333333333 0.83333333333333 0.30598544541533 200.590 > 91 136 Hg 0.16666666666667 0.16666666666667 0.80598544541533 200.590 > 91 137 Hg 0.50000000000000 0.16666666666667 0.80598544541533 200.590 > 91 138 Hg 0.83333333333333 0.16666666666667 0.80598544541533 200.590 > 91 139 Hg 0.16666666666667 0.50000000000000 0.80598544541533 200.590 > 91 140 Hg 0.50000000000000 0.50000000000000 0.80598544541533 200.590 > 91 141 Hg 0.83333333333333 0.50000000000000 0.80598544541533 200.590 > 91 142 Hg 0.16666666666667 0.83333333333333 0.80598544541533 200.590 > 91 143 Hg 0.50000000000000 0.83333333333333 0.80598544541533 200.590 > 91 144 Hg 0.83333333333333 0.83333333333333 0.80598544541533 200.590 > 91 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.0906709 0.0000000 0.0000000 0.0000000 7.0906709 0.0000000 0.0000000 0.0000000 15.7187253 -------------------------- Born effective charges -------------------------- 1 I -1.7083236 0.0000000 0.0000000 0.0000000 -1.7083236 0.0000000 0.0000000 0.0000000 -2.9213064 2 I -1.7083236 0.0000000 0.0000000 0.0000000 -1.7083236 0.0000000 0.0000000 0.0000000 -2.9213064 3 Hg 1.7083236 0.0000000 0.0000000 0.0000000 1.7083236 0.0000000 0.0000000 0.0000000 2.9213064 4 Hg 1.7083236 0.0000000 0.0000000 0.0000000 1.7083236 0.0000000 0.0000000 0.0000000 2.9213064 ---------------------------------------------------------------------------- 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[ 73, 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.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 05:10:44]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 05:10:44]------------------------- 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): [3 3 2] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4/mmm (139) ------------------------------ primitive cell ------------------------------ Lattice vectors: a -2.387467235000000 2.387467235000000 5.871517955000000 b 2.387467235000000 -2.387467235000000 5.871517955000000 c 2.387467235000000 2.387467235000000 -5.871517955000000 Atomic positions (fractional): 1 I 0.65165977896405 0.65165977896405 0.00000000000000 126.904 2 I 0.34834022103595 0.34834022103595 0.00000000000000 126.904 3 Hg 0.88802910916933 0.88802910916933 0.00000000000000 200.590 4 Hg 0.11197089083067 0.11197089083067 0.00000000000000 200.590 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.324803409999999 0.000000000000000 0.000000000000000 b 0.000000000000000 14.324803409999999 0.000000000000000 c 0.000000000000000 0.000000000000000 23.486071819999999 Atomic positions (fractional): 1 I 0.16666666666667 0.16666666666667 0.07582988948203 126.904 > 1 2 I 0.50000000000000 0.16666666666667 0.07582988948203 126.904 > 1 3 I 0.83333333333333 0.16666666666667 0.07582988948203 126.904 > 1 4 I 0.16666666666667 0.50000000000000 0.07582988948203 126.904 > 1 5 I 0.50000000000000 0.50000000000000 0.07582988948203 126.904 > 1 6 I 0.83333333333333 0.50000000000000 0.07582988948203 126.904 > 1 7 I 0.16666666666667 0.83333333333333 0.07582988948203 126.904 > 1 8 I 0.50000000000000 0.83333333333333 0.07582988948203 126.904 > 1 9 I 0.83333333333333 0.83333333333333 0.07582988948203 126.904 > 1 10 I 0.16666666666667 0.16666666666667 0.57582988948202 126.904 > 1 11 I 0.50000000000000 0.16666666666667 0.57582988948202 126.904 > 1 12 I 0.83333333333333 0.16666666666667 0.57582988948202 126.904 > 1 13 I 0.16666666666667 0.50000000000000 0.57582988948202 126.904 > 1 14 I 0.50000000000000 0.50000000000000 0.57582988948202 126.904 > 1 15 I 0.83333333333333 0.50000000000000 0.57582988948202 126.904 > 1 16 I 0.16666666666667 0.83333333333333 0.57582988948202 126.904 > 1 17 I 0.50000000000000 0.83333333333333 0.57582988948202 126.904 > 1 18 I 0.83333333333333 0.83333333333333 0.57582988948202 126.904 > 1 19 I 0.00000000000000 0.00000000000000 0.17417011051798 126.904 > 19 20 I 0.33333333333333 0.00000000000000 0.17417011051798 126.904 > 19 21 I 0.66666666666667 0.00000000000000 0.17417011051798 126.904 > 19 22 I 0.00000000000000 0.33333333333333 0.17417011051798 126.904 > 19 23 I 0.33333333333333 0.33333333333333 0.17417011051798 126.904 > 19 24 I 0.66666666666667 0.33333333333333 0.17417011051798 126.904 > 19 25 I 0.00000000000000 0.66666666666667 0.17417011051798 126.904 > 19 26 I 0.33333333333333 0.66666666666667 0.17417011051798 126.904 > 19 27 I 0.66666666666667 0.66666666666667 0.17417011051798 126.904 > 19 28 I 0.00000000000000 0.00000000000000 0.67417011051798 126.904 > 19 29 I 0.33333333333333 0.00000000000000 0.67417011051798 126.904 > 19 30 I 0.66666666666667 0.00000000000000 0.67417011051798 126.904 > 19 31 I 0.00000000000000 0.33333333333333 0.67417011051798 126.904 > 19 32 I 0.33333333333333 0.33333333333333 0.67417011051798 126.904 > 19 33 I 0.66666666666667 0.33333333333333 0.67417011051798 126.904 > 19 34 I 0.00000000000000 0.66666666666667 0.67417011051798 126.904 > 19 35 I 0.33333333333333 0.66666666666667 0.67417011051798 126.904 > 19 36 I 0.66666666666667 0.66666666666667 0.67417011051798 126.904 > 19 37 I 0.00000000000000 0.00000000000000 0.32582988948203 126.904 > 1 38 I 0.33333333333333 0.00000000000000 0.32582988948203 126.904 > 1 39 I 0.66666666666667 0.00000000000000 0.32582988948203 126.904 > 1 40 I 0.00000000000000 0.33333333333333 0.32582988948203 126.904 > 1 41 I 0.33333333333333 0.33333333333333 0.32582988948203 126.904 > 1 42 I 0.66666666666667 0.33333333333333 0.32582988948203 126.904 > 1 43 I 0.00000000000000 0.66666666666667 0.32582988948203 126.904 > 1 44 I 0.33333333333333 0.66666666666667 0.32582988948203 126.904 > 1 45 I 0.66666666666667 0.66666666666667 0.32582988948203 126.904 > 1 46 I 0.00000000000000 0.00000000000000 0.82582988948203 126.904 > 1 47 I 0.33333333333333 0.00000000000000 0.82582988948203 126.904 > 1 48 I 0.66666666666667 0.00000000000000 0.82582988948203 126.904 > 1 49 I 0.00000000000000 0.33333333333333 0.82582988948203 126.904 > 1 50 I 0.33333333333333 0.33333333333333 0.82582988948203 126.904 > 1 51 I 0.66666666666667 0.33333333333333 0.82582988948203 126.904 > 1 52 I 0.00000000000000 0.66666666666667 0.82582988948203 126.904 > 1 53 I 0.33333333333333 0.66666666666667 0.82582988948203 126.904 > 1 54 I 0.66666666666667 0.66666666666667 0.82582988948203 126.904 > 1 55 I 0.16666666666667 0.16666666666667 0.42417011051797 126.904 > 19 56 I 0.50000000000000 0.16666666666667 0.42417011051797 126.904 > 19 57 I 0.83333333333333 0.16666666666667 0.42417011051797 126.904 > 19 58 I 0.16666666666667 0.50000000000000 0.42417011051797 126.904 > 19 59 I 0.50000000000000 0.50000000000000 0.42417011051797 126.904 > 19 60 I 0.83333333333333 0.50000000000000 0.42417011051797 126.904 > 19 61 I 0.16666666666667 0.83333333333333 0.42417011051797 126.904 > 19 62 I 0.50000000000000 0.83333333333333 0.42417011051797 126.904 > 19 63 I 0.83333333333333 0.83333333333333 0.42417011051797 126.904 > 19 64 I 0.16666666666667 0.16666666666667 0.92417011051797 126.904 > 19 65 I 0.50000000000000 0.16666666666667 0.92417011051797 126.904 > 19 66 I 0.83333333333333 0.16666666666667 0.92417011051797 126.904 > 19 67 I 0.16666666666667 0.50000000000000 0.92417011051797 126.904 > 19 68 I 0.50000000000000 0.50000000000000 0.92417011051797 126.904 > 19 69 I 0.83333333333333 0.50000000000000 0.92417011051797 126.904 > 19 70 I 0.16666666666667 0.83333333333333 0.92417011051797 126.904 > 19 71 I 0.50000000000000 0.83333333333333 0.92417011051797 126.904 > 19 72 I 0.83333333333333 0.83333333333333 0.92417011051797 126.904 > 19 73 Hg 0.16666666666667 0.16666666666667 0.19401455458467 200.590 > 73 74 Hg 0.50000000000000 0.16666666666667 0.19401455458467 200.590 > 73 75 Hg 0.83333333333333 0.16666666666667 0.19401455458467 200.590 > 73 76 Hg 0.16666666666667 0.50000000000000 0.19401455458467 200.590 > 73 77 Hg 0.50000000000000 0.50000000000000 0.19401455458467 200.590 > 73 78 Hg 0.83333333333333 0.50000000000000 0.19401455458467 200.590 > 73 79 Hg 0.16666666666667 0.83333333333333 0.19401455458467 200.590 > 73 80 Hg 0.50000000000000 0.83333333333333 0.19401455458467 200.590 > 73 81 Hg 0.83333333333333 0.83333333333333 0.19401455458467 200.590 > 73 82 Hg 0.16666666666667 0.16666666666667 0.69401455458467 200.590 > 73 83 Hg 0.50000000000000 0.16666666666667 0.69401455458467 200.590 > 73 84 Hg 0.83333333333333 0.16666666666667 0.69401455458467 200.590 > 73 85 Hg 0.16666666666667 0.50000000000000 0.69401455458467 200.590 > 73 86 Hg 0.50000000000000 0.50000000000000 0.69401455458467 200.590 > 73 87 Hg 0.83333333333333 0.50000000000000 0.69401455458467 200.590 > 73 88 Hg 0.16666666666667 0.83333333333333 0.69401455458467 200.590 > 73 89 Hg 0.50000000000000 0.83333333333333 0.69401455458467 200.590 > 73 90 Hg 0.83333333333333 0.83333333333333 0.69401455458467 200.590 > 73 91 Hg 0.00000000000000 0.00000000000000 0.05598544541533 200.590 > 91 92 Hg 0.33333333333333 0.00000000000000 0.05598544541533 200.590 > 91 93 Hg 0.66666666666667 0.00000000000000 0.05598544541533 200.590 > 91 94 Hg 0.00000000000000 0.33333333333333 0.05598544541533 200.590 > 91 95 Hg 0.33333333333333 0.33333333333333 0.05598544541533 200.590 > 91 96 Hg 0.66666666666667 0.33333333333333 0.05598544541533 200.590 > 91 97 Hg 0.00000000000000 0.66666666666667 0.05598544541533 200.590 > 91 98 Hg 0.33333333333333 0.66666666666667 0.05598544541533 200.590 > 91 99 Hg 0.66666666666667 0.66666666666667 0.05598544541533 200.590 > 91 100 Hg 0.00000000000000 0.00000000000000 0.55598544541533 200.590 > 91 101 Hg 0.33333333333333 0.00000000000000 0.55598544541533 200.590 > 91 102 Hg 0.66666666666667 0.00000000000000 0.55598544541533 200.590 > 91 103 Hg 0.00000000000000 0.33333333333333 0.55598544541533 200.590 > 91 104 Hg 0.33333333333333 0.33333333333333 0.55598544541533 200.590 > 91 105 Hg 0.66666666666667 0.33333333333333 0.55598544541533 200.590 > 91 106 Hg 0.00000000000000 0.66666666666667 0.55598544541533 200.590 > 91 107 Hg 0.33333333333333 0.66666666666667 0.55598544541533 200.590 > 91 108 Hg 0.66666666666667 0.66666666666667 0.55598544541533 200.590 > 91 109 Hg 0.00000000000000 0.00000000000000 0.44401455458467 200.590 > 73 110 Hg 0.33333333333333 0.00000000000000 0.44401455458467 200.590 > 73 111 Hg 0.66666666666667 0.00000000000000 0.44401455458467 200.590 > 73 112 Hg 0.00000000000000 0.33333333333333 0.44401455458467 200.590 > 73 113 Hg 0.33333333333333 0.33333333333333 0.44401455458467 200.590 > 73 114 Hg 0.66666666666667 0.33333333333333 0.44401455458467 200.590 > 73 115 Hg 0.00000000000000 0.66666666666667 0.44401455458467 200.590 > 73 116 Hg 0.33333333333333 0.66666666666667 0.44401455458467 200.590 > 73 117 Hg 0.66666666666667 0.66666666666667 0.44401455458467 200.590 > 73 118 Hg 0.00000000000000 0.00000000000000 0.94401455458467 200.590 > 73 119 Hg 0.33333333333333 0.00000000000000 0.94401455458467 200.590 > 73 120 Hg 0.66666666666667 0.00000000000000 0.94401455458467 200.590 > 73 121 Hg 0.00000000000000 0.33333333333333 0.94401455458467 200.590 > 73 122 Hg 0.33333333333333 0.33333333333333 0.94401455458467 200.590 > 73 123 Hg 0.66666666666667 0.33333333333333 0.94401455458467 200.590 > 73 124 Hg 0.00000000000000 0.66666666666667 0.94401455458467 200.590 > 73 125 Hg 0.33333333333333 0.66666666666667 0.94401455458467 200.590 > 73 126 Hg 0.66666666666667 0.66666666666667 0.94401455458467 200.590 > 73 127 Hg 0.16666666666667 0.16666666666667 0.30598544541533 200.590 > 91 128 Hg 0.50000000000000 0.16666666666667 0.30598544541533 200.590 > 91 129 Hg 0.83333333333333 0.16666666666667 0.30598544541533 200.590 > 91 130 Hg 0.16666666666667 0.50000000000000 0.30598544541533 200.590 > 91 131 Hg 0.50000000000000 0.50000000000000 0.30598544541533 200.590 > 91 132 Hg 0.83333333333333 0.50000000000000 0.30598544541533 200.590 > 91 133 Hg 0.16666666666667 0.83333333333333 0.30598544541533 200.590 > 91 134 Hg 0.50000000000000 0.83333333333333 0.30598544541533 200.590 > 91 135 Hg 0.83333333333333 0.83333333333333 0.30598544541533 200.590 > 91 136 Hg 0.16666666666667 0.16666666666667 0.80598544541533 200.590 > 91 137 Hg 0.50000000000000 0.16666666666667 0.80598544541533 200.590 > 91 138 Hg 0.83333333333333 0.16666666666667 0.80598544541533 200.590 > 91 139 Hg 0.16666666666667 0.50000000000000 0.80598544541533 200.590 > 91 140 Hg 0.50000000000000 0.50000000000000 0.80598544541533 200.590 > 91 141 Hg 0.83333333333333 0.50000000000000 0.80598544541533 200.590 > 91 142 Hg 0.16666666666667 0.83333333333333 0.80598544541533 200.590 > 91 143 Hg 0.50000000000000 0.83333333333333 0.80598544541533 200.590 > 91 144 Hg 0.83333333333333 0.83333333333333 0.80598544541533 200.590 > 91 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 7.0906709 0.0000000 0.0000000 0.0000000 7.0906709 0.0000000 0.0000000 0.0000000 15.7187253 -------------------------- Born effective charges -------------------------- 1 I -1.7083236 0.0000000 0.0000000 0.0000000 -1.7083236 0.0000000 0.0000000 0.0000000 -2.9213064 2 I -1.7083236 0.0000000 0.0000000 0.0000000 -1.7083236 0.0000000 0.0000000 0.0000000 -2.9213064 3 Hg 1.7083236 0.0000000 0.0000000 0.0000000 1.7083236 0.0000000 0.0000000 0.0000000 2.9213064 4 Hg 1.7083236 0.0000000 0.0000000 0.0000000 1.7083236 0.0000000 0.0000000 0.0000000 2.9213064 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (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 ... Generalized regular grid: [ 2 10 40 ] Grid generation matrix: [ 0 10 10 ] [ 10 0 10 ] [ 4 4 0 ] 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.81, Number of G-points: 309, Lambda: 0.27 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/93) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 93 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.651 ( 0.000 0.000 0.000) 0.000 0.651 ( 0.000 0.000 0.000) 0.000 1.014 ( 0.000 0.000 0.000) 0.000 1.014 ( 0.000 0.000 0.000) 0.000 2.222 ( 0.000 0.000 0.000) 0.000 2.222 ( 0.000 0.000 0.000) 0.000 3.088 ( 0.000 0.000 0.000) 0.000 3.413 ( 0.000 0.000 0.000) 0.000 5.379 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/93) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 116 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.716 ( -0.000 0.000 0.000) 0.000 0.774 ( -0.000 0.000 -0.000) 0.000 0.919 ( -0.000 0.000 0.000) 0.000 1.075 ( 0.000 0.000 -0.000) 0.000 1.403 ( -0.000 0.000 0.000) 0.000 1.836 ( -0.000 0.000 0.000) 0.000 1.866 ( 0.000 0.000 -0.000) 0.000 1.887 ( -0.000 0.000 0.000) 0.000 1.937 ( -0.000 0.000 0.000) 0.000 2.771 ( 0.000 0.000 -0.000) 0.000 3.994 ( -0.000 0.000 0.000) 0.000 5.333 ( 0.000 0.000 -0.000) 0.000 ======================= Grid point 3 (3/93) ======================= q-point: (-0.50 0.00 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.673 ( -1.738 -2.793 -1.483) 3.609 0.768 ( -0.211 -0.559 -0.639) 0.875 0.867 ( -0.631 -3.021 -2.186) 3.782 1.047 ( -1.259 -1.611 -0.495) 2.103 1.322 ( -4.205 -3.658 -0.250) 5.579 1.586 ( -7.646 -7.095 -1.481) 10.536 1.817 ( -5.752 4.074 1.245) 7.158 1.935 ( 4.598 -4.567 -1.130) 6.578 2.083 ( 5.561 3.028 1.188) 6.442 2.764 ( -0.645 -0.043 0.169) 0.669 3.999 ( -0.288 0.651 -1.402) 1.572 5.337 ( 0.026 0.272 -0.437) 0.515 ======================= Grid point 4 (4/93) ======================= q-point: ( 0.00 0.00 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 150 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.084 ( 1.383 1.383 0.000) 1.956 0.323 ( 7.299 7.299 0.000) 10.322 0.590 ( 12.381 12.381 0.000) 17.510 0.620 ( -1.003 -1.003 0.000) 1.418 0.793 ( 5.322 5.322 0.000) 7.526 1.043 ( 0.968 0.968 0.000) 1.368 2.129 ( -0.191 -0.191 0.000) 0.270 2.141 ( -3.500 -3.500 0.000) 4.950 2.221 ( 0.148 0.148 0.000) 0.209 3.082 ( -0.322 -0.322 0.000) 0.455 3.451 ( 1.659 1.659 0.000) 2.346 5.375 ( -0.145 -0.145 0.000) 0.206 ======================= Grid point 5 (5/93) ======================= q-point: ( 0.00 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 153 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.117 ( 0.576 0.576 0.000) 0.814 0.586 ( 5.108 5.108 0.000) 7.223 0.616 ( 1.040 1.040 0.000) 1.471 0.980 ( 5.645 5.645 0.000) 7.983 0.987 ( 3.211 3.211 0.000) 4.541 1.063 ( -0.048 -0.048 0.000) 0.067 1.957 ( -5.208 -5.208 0.000) 7.365 2.096 ( -1.205 -1.205 0.000) 1.704 2.250 ( 1.304 1.304 0.000) 1.844 3.063 ( -0.527 -0.527 0.000) 0.745 3.534 ( 2.129 2.129 0.000) 3.010 5.370 ( -0.068 -0.068 0.000) 0.096 ======================= Grid point 6 (6/93) ======================= q-point: (-0.50 0.00 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.502 ( -4.206 -6.082 -6.054) 9.557 0.736 ( -0.054 -1.033 -2.246) 2.473 0.826 ( 1.053 -1.460 -1.337) 2.242 0.988 ( -0.474 -0.982 0.290) 1.128 1.095 ( -5.398 -4.561 1.026) 7.141 1.371 ( -2.266 -2.362 -3.870) 5.068 1.805 ( -5.251 5.405 0.979) 7.599 1.921 ( 5.236 -5.986 -0.554) 7.972 2.228 ( 3.479 0.808 0.655) 3.631 2.746 ( -0.890 0.128 0.197) 0.921 4.005 ( -0.723 0.897 -2.070) 2.369 5.343 ( -0.127 0.316 -0.712) 0.789 ======================= Grid point 8 (7/93) ======================= q-point: ( 0.00 0.00 0.30) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 150 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.181 ( 2.709 2.709 0.000) 3.832 0.673 ( 1.057 1.057 0.000) 1.494 0.756 ( 3.122 3.122 0.000) 4.416 1.020 ( -3.856 -3.856 0.000) 5.454 1.045 ( -0.283 -0.283 0.000) 0.401 1.066 ( 0.661 0.661 0.000) 0.935 1.681 ( -8.355 -8.355 0.000) 11.815 2.083 ( 1.287 1.287 0.000) 1.820 2.315 ( 1.475 1.475 0.000) 2.085 3.040 ( -0.566 -0.566 0.000) 0.800 3.618 ( 1.779 1.779 0.000) 2.516 5.371 ( 0.121 0.121 0.000) 0.172 ======================= Grid point 9 (8/93) ======================= q-point: ( 0.00 0.00 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 153 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.306 ( 2.581 2.581 0.000) 3.650 0.678 ( -0.708 -0.708 0.000) 1.001 0.723 ( -8.772 -8.772 0.000) 12.406 0.844 ( 0.579 0.579 0.000) 0.819 0.990 ( -1.849 -1.849 0.000) 2.615 1.129 ( 1.987 1.987 0.000) 2.810 1.277 ( -9.512 -9.512 0.000) 13.453 2.195 ( 3.277 3.277 0.000) 4.634 2.356 ( 0.458 0.458 0.000) 0.647 3.014 ( -0.655 -0.655 0.000) 0.926 3.678 ( 0.999 0.999 0.000) 1.413 5.378 ( 0.155 0.155 0.000) 0.220 ======================= Grid point 13 (9/93) ======================= q-point: ( 0.00 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 84 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.363 ( 0.000 0.000 0.000) 0.000 0.477 ( 0.000 0.000 0.000) 0.000 0.653 ( 0.000 0.000 0.000) 0.000 0.786 ( 0.000 0.000 0.000) 0.000 0.956 ( -0.000 -0.000 0.000) 0.000 1.070 ( 0.000 0.000 0.000) 0.000 1.181 ( 0.000 0.000 0.000) 0.000 2.279 ( 0.000 0.000 0.000) 0.000 2.363 ( 0.000 0.000 0.000) 0.000 2.996 ( 0.000 0.000 0.000) 0.000 3.699 ( 0.000 0.000 0.000) 0.000 5.381 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 23 (10/93) ======================= q-point: (-0.23 -0.03 0.03) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.498 ( 0.000 -6.437 -1.715) 6.662 0.504 ( 0.000 -8.538 -4.172) 9.503 0.695 ( 0.000 -1.491 -1.033) 1.813 0.832 ( 0.000 -7.921 -3.258) 8.565 0.931 ( 0.000 -7.179 -16.452) 17.950 1.254 ( 0.000 -10.633 0.228) 10.635 1.656 ( -0.000 -10.894 20.101) 22.863 2.111 ( -0.000 3.724 1.611) 4.057 2.146 ( -0.000 2.849 0.609) 2.914 2.900 ( -0.000 0.526 9.695) 9.709 3.874 ( 0.000 3.549 -17.342) 17.702 5.356 ( -0.000 0.392 0.761) 0.856 ======================= Grid point 24 (11/93) ======================= q-point: ( 0.28 -0.03 0.03) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.609 ( -0.000 4.344 0.522) 4.375 0.675 ( -0.000 7.706 2.589) 8.129 0.727 ( -0.000 1.536 0.997) 1.832 0.986 ( -0.000 6.103 -0.326) 6.111 1.085 ( -0.000 7.368 13.939) 15.766 1.485 ( -0.000 10.861 0.067) 10.861 1.807 ( 0.000 4.133 -11.061) 11.808 2.038 ( 0.000 -2.988 -1.804) 3.490 2.077 ( 0.000 -3.741 -0.458) 3.769 2.893 ( 0.000 -0.186 -9.372) 9.373 3.834 ( -0.000 -0.503 15.223) 15.231 5.349 ( 0.000 -0.220 -1.019) 1.043 ======================= Grid point 25 (12/93) ======================= q-point: ( 0.28 -0.03 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.651 ( -2.975 4.379 0.421) 5.311 0.752 ( -0.300 1.347 0.718) 1.555 0.829 ( 0.640 4.294 1.325) 4.539 1.068 ( -1.099 2.436 -1.686) 3.160 1.190 ( -2.715 4.209 10.407) 11.549 1.559 ( -9.523 7.957 0.009) 12.410 1.835 ( -6.088 -4.552 -0.960) 7.662 1.956 ( 4.046 1.881 -2.521) 5.125 2.112 ( 5.689 -2.572 -1.167) 6.351 2.885 ( -0.464 -0.103 -9.486) 9.498 3.842 ( -0.127 1.190 11.990) 12.049 5.348 ( 0.115 0.030 -1.154) 1.160 ======================= Grid point 26 (13/93) ======================= q-point: (-0.23 -0.03 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.215 ( 0.305 -0.305 -8.307) 8.318 0.490 ( 6.859 -6.859 -4.778) 10.812 0.641 ( -0.513 0.513 -1.905) 2.039 0.779 ( 5.847 -5.847 -5.023) 9.675 0.859 ( 5.271 -5.271 -17.436) 18.962 1.034 ( 1.146 -1.146 0.863) 1.837 1.538 ( 13.389 -13.389 22.034) 29.053 2.127 ( -3.553 3.553 1.260) 5.180 2.211 ( 0.174 -0.174 0.789) 0.827 2.907 ( -0.633 0.633 10.362) 10.400 3.931 ( -4.252 4.252 -16.272) 17.347 5.361 ( -0.334 0.334 0.482) 0.674 ======================= Grid point 28 (14/93) ======================= q-point: ( 0.28 -0.03 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.670 ( -2.859 2.728 -0.258) 3.960 0.761 ( -0.527 0.380 0.098) 0.657 0.851 ( -3.938 0.840 -1.460) 4.283 1.046 ( -3.439 -0.183 -1.096) 3.614 1.158 ( -4.319 0.960 8.723) 9.781 1.507 ( -10.382 3.770 -0.612) 11.062 1.520 ( -9.408 -3.669 1.991) 10.292 2.112 ( 4.655 3.504 -0.145) 5.828 2.164 ( 5.001 -3.530 -1.517) 6.306 2.870 ( -0.792 -0.070 -9.792) 9.824 3.868 ( -0.359 1.611 9.453) 9.596 5.355 ( 0.169 0.237 -1.314) 1.346 ======================= Grid point 29 (15/93) ======================= q-point: (-0.72 -0.03 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.598 ( -4.261 -1.902 -2.292) 5.199 0.760 ( -0.691 0.441 -0.364) 0.897 0.773 ( -5.392 4.238 -1.652) 7.054 0.932 ( -5.795 -0.107 1.030) 5.887 1.053 ( -2.901 -3.657 1.254) 4.833 1.219 ( -5.746 0.544 7.665) 9.595 1.519 ( -7.614 8.608 -0.988) 11.534 2.175 ( 4.529 -4.575 -1.777) 6.679 2.261 ( 2.731 1.893 0.190) 3.329 2.850 ( -0.866 0.067 -9.989) 10.027 3.887 ( -0.872 1.328 8.111) 8.265 5.362 ( 0.013 0.205 -1.439) 1.454 ======================= Grid point 30 (16/93) ======================= q-point: (-0.23 -0.03 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.526 ( 7.509 -6.028 -1.896) 9.815 0.709 ( 2.625 -1.236 -1.255) 3.161 0.731 ( 5.003 3.173 -2.817) 6.560 1.000 ( 4.062 0.994 0.868) 4.271 1.082 ( 5.650 -0.368 -11.864) 13.146 1.368 ( 10.016 -8.317 -0.474) 13.028 1.875 ( 2.459 3.595 8.017) 9.124 1.956 ( -3.354 -3.795 2.041) 5.460 2.162 ( -2.138 4.619 0.956) 5.179 2.888 ( -0.195 -0.480 9.630) 9.644 3.830 ( -0.201 -0.364 -14.895) 14.901 5.350 ( -0.165 0.047 1.071) 1.085 ======================= Grid point 31 (17/93) ======================= q-point: (-0.23 -0.03 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.561 ( 7.361 -5.665 -1.766) 9.455 0.741 ( 1.899 -0.852 -0.924) 2.278 0.832 ( 1.184 -0.688 -0.300) 1.402 1.038 ( 0.869 -1.785 1.702) 2.615 1.119 ( 2.512 -3.523 -9.011) 9.996 1.350 ( 7.115 -9.904 -1.459) 12.282 1.682 ( -7.510 -8.186 -0.072) 11.109 2.041 ( 3.145 3.818 1.672) 5.221 2.223 ( -2.036 4.738 1.316) 5.322 2.871 ( -0.069 -0.799 10.156) 10.187 3.839 ( 1.152 -0.097 -11.635) 11.692 5.352 ( 0.070 0.189 1.245) 1.261 ======================= Grid point 34 (18/93) ======================= q-point: (-0.23 -0.03 0.53) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.588 ( 5.253 -5.030 -0.577) 7.296 0.724 ( -1.845 -7.551 2.002) 8.026 0.747 ( 1.103 -1.240 -0.721) 1.809 0.952 ( -1.637 -5.536 -0.418) 5.788 1.088 ( 0.683 -2.177 -6.565) 6.950 1.288 ( -3.972 -9.475 -0.532) 10.287 1.340 ( 4.926 -7.804 -1.728) 9.389 2.200 ( 3.608 3.694 -0.148) 5.166 2.258 ( -3.175 3.862 1.589) 5.246 2.851 ( -0.076 -1.026 10.540) 10.590 3.860 ( 1.185 -0.384 -9.549) 9.630 5.358 ( 0.145 0.173 1.384) 1.402 ======================= Grid point 35 (19/93) ======================= q-point: ( 0.78 -0.03 -0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.522 ( 1.992 -1.992 2.145) 3.541 0.633 ( 7.369 -7.369 0.252) 10.424 0.766 ( -0.494 0.494 0.317) 0.767 0.809 ( 4.511 -4.511 -2.187) 6.744 1.036 ( -1.125 1.125 -0.243) 1.610 1.134 ( 2.478 -2.478 -8.262) 8.974 1.355 ( 7.860 -7.860 1.203) 11.181 2.262 ( -3.657 3.657 1.751) 5.461 2.290 ( 0.006 -0.006 -0.447) 0.447 2.835 ( 0.519 -0.519 10.729) 10.754 3.869 ( 0.825 -0.825 -8.818) 8.895 5.362 ( 0.007 -0.007 1.441) 1.441 ======================= Grid point 40 (20/93) ======================= q-point: ( 0.28 -0.03 -0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.345 ( -9.162 -8.122 4.498) 13.044 0.619 ( -5.638 5.028 4.351) 8.718 0.658 ( -1.959 -1.804 2.164) 3.432 0.900 ( -5.387 4.447 0.562) 7.008 0.969 ( -5.061 3.773 14.112) 15.460 1.156 ( -9.409 -6.079 0.506) 11.213 1.766 ( -7.892 7.275 -17.089) 20.180 2.043 ( 4.315 -4.139 -1.443) 6.151 2.198 ( 1.269 2.621 -0.808) 3.022 2.895 ( 0.488 -0.457 -9.766) 9.789 3.856 ( 2.464 -1.391 17.191) 17.422 5.355 ( 0.288 -0.091 -0.868) 0.919 ======================= Grid point 43 (21/93) ======================= q-point: (-0.45 -0.05 0.05) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 232 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.691 ( 0.000 -2.195 -0.780) 2.329 0.770 ( 0.000 -0.471 -0.662) 0.812 0.873 ( 0.000 -4.098 -2.965) 5.058 1.061 ( 0.000 -1.723 -0.624) 1.833 1.367 ( 0.000 -3.354 -0.568) 3.402 1.699 ( 0.000 -8.750 -0.499) 8.764 1.835 ( 0.000 -3.449 -1.647) 3.822 1.942 ( -0.000 1.056 2.358) 2.584 1.981 ( -0.000 4.705 0.287) 4.713 2.771 ( -0.000 -0.033 0.260) 0.262 4.002 ( 0.000 0.749 -1.449) 1.631 5.337 ( 0.000 0.317 -0.475) 0.572 ======================= Grid point 44 (22/93) ======================= q-point: ( 0.05 -0.05 0.05) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.233 ( -0.000 11.059 0.000) 11.059 0.278 ( -0.000 11.775 0.000) 11.775 0.334 ( -0.000 15.479 0.000) 15.479 0.659 ( -0.000 0.785 0.000) 0.785 0.708 ( -0.000 5.062 0.000) 5.062 1.078 ( -0.000 6.080 0.000) 6.080 2.110 ( 0.000 -1.654 -0.000) 1.654 2.195 ( 0.000 -2.454 -0.000) 2.454 2.204 ( 0.000 -1.688 -0.000) 1.688 3.085 ( 0.000 -0.330 -0.000) 0.330 3.434 ( -0.000 1.895 0.000) 1.895 5.377 ( 0.000 -0.185 -0.000) 0.185 ======================= Grid point 45 (23/93) ======================= q-point: ( 0.05 -0.05 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.308 ( -11.456 11.974 0.000) 16.572 0.494 ( 3.302 8.488 0.000) 9.108 0.636 ( -2.612 2.277 0.000) 3.465 0.795 ( 9.810 7.399 -0.000) 12.288 0.905 ( 4.353 5.244 0.000) 6.815 1.159 ( -6.753 9.107 0.000) 11.337 2.056 ( -4.292 -4.501 -0.000) 6.219 2.094 ( 1.575 -2.640 -0.000) 3.074 2.209 ( 2.613 -1.237 -0.000) 2.891 3.072 ( -0.271 -0.641 -0.000) 0.696 3.500 ( 1.116 2.924 0.000) 3.130 5.371 ( -0.043 -0.228 -0.000) 0.232 ======================= Grid point 46 (24/93) ======================= q-point: (-0.45 -0.05 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.591 ( -3.560 -5.163 -3.456) 7.161 0.739 ( -0.352 -2.122 -1.570) 2.663 0.792 ( 1.820 -3.768 -3.338) 5.353 0.994 ( -0.041 -3.393 -0.893) 3.509 1.219 ( -5.188 -5.875 -0.087) 7.839 1.428 ( -3.495 -7.692 -2.715) 8.874 1.795 ( 5.660 -7.385 -0.579) 9.323 1.917 ( -4.522 4.042 1.603) 6.274 2.142 ( 4.537 2.487 0.686) 5.220 2.762 ( -0.581 -0.172 0.491) 0.780 4.019 ( -0.535 1.262 -2.141) 2.542 5.345 ( -0.105 0.535 -0.914) 1.064 ======================= Grid point 47 (25/93) ======================= q-point: (-0.45 -0.05 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 227 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.413 ( 1.104 -1.104 -7.859) 8.013 0.720 ( 0.648 -0.648 -3.570) 3.686 0.773 ( 3.388 -3.388 -2.033) 5.205 0.961 ( 1.969 -1.969 0.033) 2.785 1.041 ( 0.028 -0.028 1.325) 1.326 1.322 ( 2.258 -2.258 -5.213) 6.114 1.771 ( 7.594 -7.594 0.108) 10.740 1.917 ( -4.619 4.619 1.330) 6.666 2.225 ( 1.144 -1.144 0.397) 1.665 2.750 ( -0.294 0.294 0.641) 0.764 4.027 ( -1.162 1.162 -2.192) 2.740 5.351 ( -0.412 0.412 -1.090) 1.236 ======================= Grid point 48 (26/93) ======================= q-point: ( 0.05 -0.05 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.324 ( -10.849 12.265 0.000) 16.375 0.665 ( -0.229 3.283 0.000) 3.291 0.702 ( 1.797 5.604 0.000) 5.885 0.992 ( 2.541 -0.533 -0.000) 2.597 1.038 ( 0.935 1.786 0.000) 2.016 1.185 ( -7.909 7.416 0.000) 10.842 1.835 ( -6.429 -6.535 -0.000) 9.167 2.066 ( 0.732 -1.057 -0.000) 1.286 2.261 ( 3.543 -0.288 -0.000) 3.555 3.050 ( -0.373 -0.699 -0.000) 0.792 3.588 ( 1.149 2.889 0.000) 3.109 5.369 ( 0.117 -0.046 -0.000) 0.126 ======================= Grid point 49 (27/93) ======================= q-point: ( 0.05 -0.05 0.35) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.380 ( -7.893 11.567 0.000) 14.003 0.704 ( -1.403 1.665 0.000) 2.177 0.826 ( 0.761 3.387 0.000) 3.471 0.888 ( -6.008 -6.709 -0.000) 9.006 1.031 ( -1.919 -0.844 0.000) 2.096 1.164 ( -5.727 6.992 0.000) 9.038 1.486 ( -9.944 -10.187 -0.000) 14.236 2.128 ( 3.025 2.993 -0.000) 4.256 2.320 ( 2.907 -0.969 -0.000) 3.064 3.027 ( -0.499 -0.718 -0.000) 0.874 3.661 ( 0.571 2.279 0.000) 2.350 5.374 ( 0.232 0.117 -0.000) 0.260 ======================= Grid point 52 (28/93) ======================= q-point: ( 0.05 -0.05 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.448 ( -6.558 8.813 0.000) 10.986 0.573 ( -8.271 -4.339 0.000) 9.340 0.691 ( -2.696 2.331 0.000) 3.564 0.825 ( -3.183 -2.626 0.000) 4.127 0.973 ( -1.642 0.697 0.000) 1.783 1.112 ( -4.440 -5.792 -0.000) 7.298 1.222 ( -3.877 6.317 0.000) 7.412 2.257 ( 1.995 2.400 0.000) 3.121 2.340 ( 2.126 -1.894 -0.000) 2.847 3.001 ( -0.394 -0.558 -0.000) 0.683 3.703 ( -0.335 1.362 0.000) 1.403 5.380 ( 0.130 0.060 -0.000) 0.144 ======================= Grid point 54 (29/93) ======================= q-point: ( 0.55 -0.05 -0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 210 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.705 ( 0.000 -1.108 -0.000) 1.108 0.770 ( 0.000 -0.295 -0.000) 0.295 0.902 ( 0.000 -1.745 -0.000) 1.745 1.063 ( 0.000 -1.187 0.000) 1.187 1.361 ( 0.000 -3.759 -0.000) 3.759 1.711 ( -0.000 -8.829 0.000) 8.829 1.725 ( 0.000 -7.466 -0.000) 7.466 2.014 ( -0.000 6.175 0.000) 6.175 2.019 ( 0.000 4.854 -0.000) 4.854 2.764 ( 0.000 -0.628 0.000) 0.628 3.992 ( 0.000 -0.195 -0.000) 0.195 5.334 ( 0.000 0.069 0.000) 0.069 ======================= Grid point 56 (30/93) ======================= q-point: ( 0.55 -0.05 -0.35) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 404 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.617 ( -4.310 -3.642 3.037) 6.408 0.760 ( -0.848 -0.575 0.726) 1.255 0.837 ( 0.004 -2.436 0.620) 2.514 1.011 ( -1.066 -2.158 -0.056) 2.407 1.204 ( -5.008 -6.395 -0.702) 8.153 1.425 ( -3.138 -7.388 2.400) 8.378 1.678 ( 6.587 -7.247 -0.842) 9.829 2.035 ( -4.813 4.640 1.467) 6.844 2.193 ( 2.521 4.740 -1.167) 5.494 2.745 ( -0.016 -1.145 0.011) 1.145 3.991 ( 0.444 -0.428 1.232) 1.378 5.338 ( 0.161 0.057 0.337) 0.378 ======================= Grid point 59 (31/93) ======================= q-point: ( 0.55 -0.05 -0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 225 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.456 ( -0.633 0.633 7.100) 7.156 0.736 ( -0.122 0.122 2.132) 2.139 0.823 ( 1.554 -1.554 -0.114) 2.201 0.975 ( 1.006 -1.006 -1.448) 2.030 1.032 ( -0.030 0.030 -1.470) 1.471 1.325 ( 2.083 -2.083 3.721) 4.746 1.677 ( 6.980 -6.980 -0.482) 9.883 2.029 ( -4.979 4.979 1.369) 7.173 2.279 ( -1.163 1.163 -0.757) 1.810 2.728 ( 0.696 -0.696 0.047) 0.985 3.991 ( 0.570 -0.570 1.659) 1.844 5.340 ( 0.088 -0.088 0.457) 0.474 ======================= Grid point 65 (32/93) ======================= q-point: ( 0.32 -0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.683 ( -0.000 2.787 0.061) 2.788 0.755 ( -0.000 1.078 0.598) 1.233 0.819 ( -0.000 5.851 1.590) 6.063 1.080 ( 0.000 3.007 -1.622) 3.416 1.220 ( -0.000 4.971 10.930) 12.007 1.693 ( -0.000 8.858 0.041) 8.858 1.852 ( 0.000 0.672 -4.224) 4.277 1.989 ( 0.000 -4.695 -0.395) 4.711 1.991 ( 0.000 -1.608 -1.337) 2.091 2.890 ( 0.000 -0.135 -9.215) 9.216 3.843 ( -0.000 1.160 12.120) 12.175 5.347 ( 0.000 0.014 -1.116) 1.116 ======================= Grid point 66 (33/93) ======================= q-point: (-0.18 -0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.323 ( 0.000 -8.430 -6.212) 10.472 0.335 ( 0.000 -8.734 -4.708) 9.922 0.669 ( 0.000 -0.898 -0.854) 1.239 0.694 ( 0.000 -4.363 -7.165) 8.389 0.792 ( 0.000 -5.702 -20.177) 20.967 1.069 ( -0.000 -6.403 0.876) 6.463 1.325 ( -0.000 -20.508 20.280) 28.842 2.181 ( -0.000 2.661 1.145) 2.897 2.194 ( -0.000 1.600 0.783) 1.781 2.915 ( -0.000 0.847 10.908) 10.941 3.989 ( 0.000 7.058 -13.153) 14.927 5.366 ( -0.000 0.559 0.257) 0.615 ======================= Grid point 68 (34/93) ======================= q-point: ( 0.32 -0.08 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.711 ( -1.124 1.259 -0.412) 1.737 0.769 ( -0.299 0.514 0.052) 0.597 0.900 ( -1.176 2.250 -0.189) 2.546 1.100 ( -1.698 0.504 -2.473) 3.041 1.245 ( -3.376 0.931 9.238) 9.880 1.697 ( -8.390 4.674 0.049) 9.604 1.725 ( -8.989 -5.271 0.027) 10.420 2.003 ( 5.799 2.625 -0.258) 6.371 2.052 ( 5.580 -3.175 -1.097) 6.513 2.883 ( -0.406 -0.107 -9.242) 9.252 3.875 ( -0.213 1.872 9.388) 9.575 5.352 ( 0.102 0.315 -1.252) 1.295 ======================= Grid point 69 (35/93) ======================= q-point: (-0.68 -0.08 0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.672 ( -2.537 -1.938 -1.841) 3.685 0.770 ( -0.423 0.239 -0.219) 0.533 0.870 ( -3.407 0.873 -2.727) 4.450 1.029 ( -3.618 -1.325 -0.834) 3.942 1.144 ( -4.958 -2.644 4.637) 7.286 1.377 ( -9.099 -4.754 4.734) 11.305 1.669 ( -6.716 6.951 -0.503) 9.679 2.074 ( 5.078 -5.107 -1.432) 7.343 2.181 ( 4.730 2.868 0.059) 5.532 2.868 ( -0.750 -0.145 -9.388) 9.419 3.904 ( -0.735 1.756 7.613) 7.847 5.362 ( 0.005 0.444 -1.436) 1.503 ======================= Grid point 71 (36/93) ======================= q-point: (-0.18 -0.08 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.376 ( 9.521 -7.559 -4.931) 13.119 0.689 ( 2.637 -0.462 -1.901) 3.283 0.790 ( 2.890 2.326 -3.430) 5.052 1.016 ( 1.252 0.432 1.683) 2.141 1.068 ( 1.858 -0.967 -5.240) 5.643 1.214 ( 6.122 -5.742 -5.030) 9.786 1.824 ( -5.948 -5.817 0.864) 8.364 1.980 ( 2.811 3.199 5.501) 6.957 2.249 ( -0.419 3.494 1.013) 3.662 2.875 ( -0.278 -0.710 10.248) 10.276 3.823 ( 0.270 -0.164 -13.964) 13.967 5.351 ( -0.034 0.088 1.182) 1.185 ======================= Grid point 72 (37/93) ======================= q-point: (-0.68 -0.08 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.577 ( -0.341 0.341 -4.092) 4.120 0.771 ( -0.647 0.647 -0.966) 1.330 0.845 ( -2.297 2.297 -2.423) 4.053 0.953 ( -2.060 2.060 0.878) 3.042 1.004 ( 0.183 -0.183 0.959) 0.993 1.259 ( -2.567 2.567 7.665) 8.481 1.679 ( -6.625 6.625 -0.906) 9.412 2.069 ( 5.406 -5.406 -1.554) 7.802 2.272 ( 0.923 -0.923 -0.045) 1.306 2.855 ( -0.376 0.376 -9.464) 9.479 3.915 ( -1.309 1.309 6.974) 7.216 5.367 ( -0.272 0.272 -1.518) 1.566 ======================= Grid point 74 (38/93) ======================= q-point: (-0.18 -0.08 0.48) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.425 ( 8.841 -6.394 -4.019) 11.628 0.718 ( 1.484 -1.242 -1.173) 2.263 0.791 ( -2.982 -3.423 1.580) 4.807 0.992 ( -1.880 -2.793 1.050) 3.526 1.067 ( 0.641 -1.092 -7.601) 7.706 1.166 ( 5.503 -6.812 -2.093) 9.004 1.496 ( -9.484 -9.402 -1.079) 13.398 2.119 ( 3.840 3.724 0.493) 5.372 2.304 ( -1.095 2.893 1.351) 3.376 2.854 ( -0.216 -0.854 10.905) 10.940 3.839 ( 0.871 0.087 -11.110) 11.144 5.356 ( 0.117 0.195 1.341) 1.360 ======================= Grid point 75 (39/93) ======================= q-point: (-0.18 -0.08 0.57) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.464 ( 7.368 -6.306 -1.547) 9.820 0.555 ( -4.209 -7.410 1.433) 8.642 0.716 ( 2.256 -1.945 -2.052) 3.617 0.833 ( -3.548 -4.782 0.374) 5.966 1.063 ( -0.069 -0.313 -4.322) 4.334 1.140 ( -3.928 -4.457 -3.671) 6.983 1.205 ( 5.567 -4.761 0.066) 7.326 2.263 ( 2.427 2.191 -0.559) 3.317 2.322 ( -1.981 2.180 1.531) 3.320 2.831 ( -0.085 -0.812 11.300) 11.330 3.854 ( 0.630 -0.253 -9.639) 9.662 5.361 ( 0.065 0.108 1.435) 1.441 ======================= Grid point 81 (40/93) ======================= q-point: ( 0.32 -0.08 -0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.237 ( -0.862 0.862 8.225) 8.315 0.648 ( -0.916 0.916 2.769) 3.057 0.715 ( -4.169 4.169 5.440) 8.022 0.976 ( -2.703 2.703 -1.037) 3.961 1.060 ( 0.025 -0.025 -0.412) 0.413 1.089 ( -4.243 4.243 11.683) 13.134 1.887 ( -3.186 3.186 -10.049) 11.013 1.955 ( 3.438 -3.438 -2.803) 5.613 2.240 ( -1.250 1.250 -0.826) 1.951 2.884 ( 0.537 -0.537 -10.014) 10.043 3.830 ( 0.894 -0.894 16.049) 16.099 5.353 ( 0.090 -0.090 -1.058) 1.065 ======================= Grid point 85 (41/93) ======================= q-point: ( 0.10 -0.10 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 158 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.458 ( -0.000 10.288 0.000) 10.288 0.479 ( -0.000 7.518 0.000) 7.518 0.639 ( -0.000 13.420 0.000) 13.420 0.683 ( -0.000 1.409 0.000) 1.409 0.841 ( -0.000 7.175 0.000) 7.175 1.257 ( -0.000 10.375 0.000) 10.375 2.062 ( 0.000 -2.827 -0.000) 2.827 2.128 ( 0.000 -3.655 -0.000) 3.655 2.154 ( 0.000 -2.963 -0.000) 2.963 3.074 ( 0.000 -0.684 -0.000) 0.684 3.488 ( -0.000 3.233 0.000) 3.233 5.372 ( 0.000 -0.290 -0.000) 0.290 ======================= Grid point 86 (42/93) ======================= q-point: (-0.40 -0.10 0.10) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.631 ( 0.000 -3.470 -1.619) 3.829 0.751 ( 0.000 -1.264 -1.083) 1.665 0.762 ( 0.000 -6.202 -5.048) 7.997 0.993 ( 0.000 -4.688 -1.520) 4.928 1.274 ( 0.000 -5.103 -1.333) 5.274 1.490 ( 0.000 -11.060 -0.394) 11.067 1.712 ( 0.000 -8.184 -1.418) 8.306 1.989 ( -0.000 3.339 2.218) 4.008 2.068 ( -0.000 3.677 0.234) 3.684 2.768 ( -0.000 -0.256 0.618) 0.669 4.025 ( 0.000 1.428 -2.107) 2.545 5.347 ( 0.000 0.630 -0.991) 1.174 ======================= Grid point 87 (43/93) ======================= q-point: (-0.40 -0.10 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.465 ( -4.942 -6.135 -5.889) 9.835 0.681 ( 3.736 -3.796 -2.944) 6.085 0.722 ( 0.436 -2.504 -3.293) 4.160 0.909 ( 2.630 -4.577 -0.903) 5.356 1.090 ( -4.151 -5.862 1.058) 7.261 1.279 ( 1.834 -6.302 -4.613) 8.022 1.605 ( 7.628 -10.236 -0.834) 12.793 2.010 ( -4.188 4.548 1.317) 6.321 2.183 ( 2.667 1.407 0.408) 3.043 2.756 ( -0.245 -0.339 1.074) 1.153 4.050 ( -0.884 1.617 -1.697) 2.505 5.359 ( -0.325 0.717 -1.359) 1.571 ======================= Grid point 88 (44/93) ======================= q-point: ( 0.10 -0.10 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.522 ( -7.253 8.211 0.000) 10.956 0.663 ( 0.619 7.103 0.000) 7.130 0.698 ( -1.200 3.430 0.000) 3.634 0.930 ( 3.137 5.884 0.000) 6.668 1.010 ( 1.523 4.703 0.000) 4.944 1.370 ( -8.902 10.191 0.000) 13.532 1.957 ( -5.363 -4.930 0.000) 7.285 2.035 ( 2.353 -2.661 -0.000) 3.552 2.173 ( 4.697 -2.062 -0.000) 5.129 3.056 ( -0.162 -0.846 -0.000) 0.862 3.570 ( 0.545 3.648 0.000) 3.688 5.367 ( 0.075 -0.188 -0.000) 0.202 ======================= Grid point 89 (45/93) ======================= q-point: ( 0.10 -0.10 0.30) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.549 ( -6.997 8.694 0.000) 11.159 0.728 ( -0.907 2.373 0.000) 2.540 0.815 ( 0.327 5.034 0.000) 5.044 0.981 ( -3.435 -0.744 0.000) 3.515 1.069 ( -1.630 1.220 0.000) 2.036 1.344 ( -10.395 7.974 0.000) 13.101 1.682 ( -8.681 -8.093 0.000) 11.868 2.067 ( 2.940 1.300 -0.000) 3.214 2.238 ( 4.813 -1.910 -0.000) 5.178 3.035 ( -0.256 -0.764 -0.000) 0.806 3.653 ( 0.251 3.279 0.000) 3.288 5.369 ( 0.236 0.058 -0.000) 0.243 ======================= Grid point 92 (46/93) ======================= q-point: ( 0.10 -0.10 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.584 ( -5.721 6.656 0.000) 8.777 0.737 ( -1.573 1.458 0.000) 2.145 0.774 ( -9.979 -2.780 0.000) 10.359 0.884 ( -1.942 1.762 0.000) 2.623 1.022 ( -2.936 0.344 0.000) 2.956 1.269 ( -10.127 -9.895 0.000) 14.158 1.341 ( -7.353 9.394 0.000) 11.930 2.199 ( 3.452 3.479 0.000) 4.901 2.277 ( 3.873 -3.121 -0.000) 4.974 3.011 ( -0.462 -0.789 -0.000) 0.915 3.712 ( -0.506 2.494 0.000) 2.545 5.376 ( 0.212 0.147 -0.000) 0.258 ======================= Grid point 93 (47/93) ======================= q-point: ( 0.10 0.90 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 143 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.544 ( -2.454 2.454 0.000) 3.471 0.647 ( -8.597 8.597 0.000) 12.157 0.767 ( 0.495 -0.495 -0.000) 0.700 0.782 ( -4.589 4.589 0.000) 6.490 1.011 ( -2.500 2.500 0.000) 3.536 1.037 ( 1.293 -1.293 -0.000) 1.828 1.374 ( -7.800 7.800 0.000) 11.031 2.281 ( 3.605 -3.605 -0.000) 5.098 2.286 ( -0.153 0.153 0.000) 0.216 2.993 ( 0.131 -0.131 -0.000) 0.186 3.733 ( -1.498 1.498 0.000) 2.118 5.381 ( 0.005 -0.005 -0.000) 0.007 ======================= Grid point 96 (48/93) ======================= q-point: ( 0.60 -0.10 -0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 210 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.668 ( 0.000 -2.501 -0.000) 2.501 0.762 ( 0.000 -0.433 -0.000) 0.433 0.838 ( 0.000 -4.551 -0.000) 4.551 1.024 ( 0.000 -2.651 0.000) 2.651 1.263 ( 0.000 -5.072 -0.000) 5.072 1.498 ( 0.000 -11.086 0.000) 11.086 1.563 ( 0.000 -7.916 0.000) 7.916 2.117 ( -0.000 4.511 0.000) 4.511 2.134 ( 0.000 5.174 -0.000) 5.174 2.745 ( 0.000 -1.181 0.000) 1.181 3.986 ( 0.000 -0.324 -0.000) 0.324 5.336 ( 0.000 0.120 0.000) 0.120 ======================= Grid point 100 (49/93) ======================= q-point: ( 0.60 -0.10 -0.30) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.524 ( -4.838 -4.901 4.535) 8.245 0.744 ( -0.401 -1.103 1.027) 1.560 0.766 ( 3.327 -4.030 -0.572) 5.257 0.953 ( 0.824 -3.567 -1.454) 3.939 1.088 ( -4.667 -3.855 -1.793) 6.313 1.288 ( 1.358 -5.462 3.009) 6.382 1.517 ( 7.951 -8.007 -0.030) 11.284 2.129 ( -4.581 4.430 1.781) 6.616 2.274 ( 1.528 2.897 -1.148) 3.470 2.718 ( 0.249 -1.295 0.099) 1.323 3.982 ( 0.284 -0.391 0.915) 1.035 5.339 ( 0.040 0.072 0.215) 0.230 ======================= Grid point 108 (50/93) ======================= q-point: (-0.13 -0.12 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.207 ( 0.000 0.000 -9.241) 9.241 0.207 ( 0.000 0.000 -9.241) 9.241 0.659 ( 0.000 0.000 -0.743) 0.743 0.659 ( 0.000 0.000 -0.743) 0.743 0.719 ( 0.000 0.000 -31.087) 31.087 1.001 ( -0.000 -0.000 1.241) 1.241 1.001 ( -0.000 -0.000 1.241) 1.241 2.211 ( -0.000 -0.000 0.860) 0.860 2.211 ( -0.000 -0.000 0.860) 0.860 2.927 ( -0.000 -0.000 12.027) 12.027 4.092 ( 0.000 0.000 -3.437) 3.437 5.374 ( -0.000 -0.000 0.111) 0.111 ======================= Grid point 109 (51/93) ======================= q-point: ( 0.38 -0.12 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.722 ( 0.000 0.747 -0.385) 0.840 0.773 ( -0.000 0.611 0.073) 0.615 0.911 ( -0.000 2.758 0.333) 2.778 1.117 ( 0.000 0.664 -2.907) 2.981 1.282 ( -0.000 0.954 9.314) 9.363 1.845 ( 0.000 2.759 -0.747) 2.858 1.854 ( 0.000 -0.109 -1.165) 1.170 1.884 ( -0.000 -2.524 0.236) 2.535 1.958 ( 0.000 -1.793 -0.065) 1.794 2.887 ( 0.000 -0.131 -9.035) 9.036 3.877 ( -0.000 1.959 9.366) 9.568 5.351 ( 0.000 0.346 -1.227) 1.275 ======================= Grid point 110 (52/93) ======================= q-point: (-0.62 -0.12 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.706 ( -0.885 -1.447 -1.406) 2.203 0.777 ( -0.242 0.023 -0.158) 0.290 0.917 ( -1.300 -0.652 -2.682) 3.051 1.088 ( -1.951 -1.652 -2.381) 3.494 1.232 ( -3.094 -1.968 7.357) 8.220 1.577 ( -9.267 -7.731 1.429) 12.152 1.799 ( -5.557 4.233 -0.186) 6.988 1.962 ( 5.447 -5.223 -0.517) 7.564 2.071 ( 5.510 3.292 0.075) 6.419 2.880 ( -0.391 -0.180 -9.020) 9.030 3.917 ( -0.423 2.065 7.338) 7.635 5.362 ( -0.004 0.645 -1.430) 1.568 ======================= Grid point 113 (53/93) ======================= q-point: (-0.62 -0.12 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.610 ( -2.367 -3.468 -3.743) 5.625 0.776 ( 0.278 0.541 -1.350) 1.480 0.873 ( -0.553 -0.891 -3.273) 3.437 0.998 ( -1.867 -1.516 -0.465) 2.450 1.062 ( -5.209 -4.215 1.831) 6.947 1.319 ( -3.311 -1.396 6.371) 7.315 1.799 ( -4.839 5.369 -0.560) 7.249 1.950 ( 5.948 -6.692 -0.919) 9.000 2.224 ( 3.482 1.114 -0.139) 3.659 2.865 ( -0.524 -0.125 -9.075) 9.091 3.940 ( -1.085 1.666 6.244) 6.553 5.373 ( -0.287 0.616 -1.598) 1.737 ======================= Grid point 115 (54/93) ======================= q-point: (-0.13 -0.12 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.290 ( 1.633 1.633 -7.857) 8.189 0.692 ( 0.778 0.778 -1.670) 2.000 0.820 ( 0.269 0.269 -1.775) 1.815 1.017 ( -0.428 -0.428 2.210) 2.292 1.058 ( 0.413 0.413 0.816) 1.004 1.130 ( -2.124 -2.124 -10.806) 11.216 1.680 ( -7.878 -7.878 0.013) 11.142 2.045 ( 3.202 3.202 2.448) 5.148 2.301 ( 1.387 1.387 1.120) 2.258 2.861 ( -0.526 -0.526 10.807) 10.832 3.825 ( 0.375 0.375 -12.631) 12.642 5.353 ( 0.107 0.107 1.280) 1.289 ======================= Grid point 116 (55/93) ======================= q-point: (-0.13 -0.12 0.53) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.356 ( 1.229 1.229 -4.150) 4.500 0.682 ( -6.588 -6.588 1.910) 9.511 0.698 ( -0.367 -0.367 -1.789) 1.863 0.915 ( -4.154 -4.154 0.330) 5.883 1.062 ( -0.178 -0.178 -8.623) 8.627 1.108 ( 1.683 1.683 1.890) 3.040 1.299 ( -8.933 -8.933 -2.488) 12.876 2.197 ( 3.613 3.613 -0.293) 5.118 2.339 ( 0.424 0.424 1.380) 1.505 2.838 ( -0.593 -0.593 11.375) 11.406 3.843 ( 0.355 0.355 -10.435) 10.447 5.359 ( 0.140 0.140 1.407) 1.421 ======================= Grid point 119 (56/93) ======================= q-point: (-0.13 -0.12 -0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 122 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.380 ( 0.000 0.000 -1.603) 1.603 0.465 ( 0.000 0.000 1.206) 1.206 0.686 ( 0.000 0.000 -2.765) 2.765 0.784 ( 0.000 0.000 0.493) 0.493 1.066 ( 0.000 0.000 -0.059) 0.059 1.092 ( -0.000 -0.000 -9.972) 9.972 1.153 ( 0.000 0.000 2.478) 2.478 2.286 ( 0.000 0.000 -0.745) 0.745 2.346 ( 0.000 0.000 1.448) 1.448 2.821 ( 0.000 0.000 11.664) 11.664 3.851 ( 0.000 0.000 -9.676) 9.676 5.363 ( 0.000 0.000 1.454) 1.454 ======================= Grid point 129 (57/93) ======================= q-point: (-0.35 -0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 232 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.548 ( 0.000 -4.498 -2.794) 5.295 0.618 ( 0.000 -7.429 -6.018) 9.560 0.721 ( 0.000 -1.523 -1.275) 1.986 0.874 ( 0.000 -6.280 -2.167) 6.644 1.162 ( 0.000 -5.580 -1.081) 5.684 1.249 ( 0.000 -11.383 0.174) 11.384 1.507 ( 0.000 -10.716 -3.735) 11.348 2.073 ( -0.000 4.364 1.477) 4.607 2.135 ( -0.000 2.703 0.265) 2.716 2.759 ( -0.000 -0.599 1.244) 1.380 4.060 ( 0.000 1.826 -1.346) 2.269 5.362 ( 0.000 0.859 -1.447) 1.683 ======================= Grid point 130 (58/93) ======================= q-point: ( 0.15 -0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.604 ( -0.000 4.717 0.000) 4.717 0.655 ( -0.000 8.392 0.000) 8.392 0.716 ( -0.000 1.650 0.000) 1.650 0.888 ( -0.000 10.197 0.000) 10.197 0.990 ( -0.000 6.679 0.000) 6.679 1.485 ( -0.000 10.818 0.000) 10.818 1.997 ( 0.000 -3.236 -0.000) 3.236 2.053 ( 0.000 -3.298 -0.000) 3.298 2.083 ( 0.000 -3.810 -0.000) 3.810 3.057 ( 0.000 -0.914 -0.000) 0.914 3.564 ( -0.000 3.913 0.000) 3.913 5.366 ( 0.000 -0.243 -0.000) 0.243 ======================= Grid point 131 (59/93) ======================= q-point: ( 0.15 -0.15 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.653 ( -3.039 4.512 0.000) 5.440 0.743 ( -0.504 1.899 0.000) 1.965 0.810 ( 0.242 5.737 0.000) 5.743 1.039 ( -1.984 4.168 0.000) 4.616 1.097 ( -1.005 3.408 0.000) 3.553 1.564 ( -9.802 8.104 0.000) 12.719 1.849 ( -7.359 -5.470 0.000) 9.169 1.996 ( 3.953 -0.750 -0.000) 4.024 2.124 ( 5.914 -2.578 -0.000) 6.451 3.038 ( -0.064 -0.823 -0.000) 0.825 3.650 ( 0.081 3.888 0.000) 3.889 5.364 ( 0.149 0.011 -0.000) 0.150 ======================= Grid point 132 (60/93) ======================= q-point: (-0.35 -0.15 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 225 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.385 ( 0.204 -0.204 -7.327) 7.333 0.577 ( 5.724 -5.724 -4.183) 9.112 0.692 ( 0.530 -0.530 -2.636) 2.740 0.810 ( 4.604 -4.604 -0.904) 6.573 1.005 ( 1.489 -1.489 1.741) 2.732 1.159 ( 5.655 -5.655 -2.533) 8.389 1.399 ( 8.146 -8.146 -6.112) 13.041 2.100 ( -3.912 3.912 0.908) 5.607 2.197 ( 0.208 -0.208 0.339) 0.448 2.749 ( 0.333 -0.333 1.774) 1.835 4.082 ( -1.259 1.259 -0.256) 1.799 5.373 ( -0.592 0.592 -1.587) 1.794 ======================= Grid point 134 (61/93) ======================= q-point: ( 0.15 -0.15 0.35) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.676 ( -3.037 3.270 0.000) 4.463 0.760 ( -0.661 0.782 0.000) 1.024 0.907 ( -0.497 3.600 0.000) 3.635 0.965 ( -7.851 -0.392 0.000) 7.861 1.089 ( -3.202 0.637 0.000) 3.265 1.501 ( -11.245 -8.996 0.000) 14.401 1.517 ( -9.078 8.265 0.000) 12.277 2.116 ( 4.493 3.039 -0.000) 5.424 2.181 ( 5.166 -3.535 -0.000) 6.260 3.019 ( -0.281 -0.752 -0.000) 0.803 3.722 ( -0.469 3.297 0.000) 3.331 5.372 ( 0.212 0.256 0.000) 0.333 ======================= Grid point 135 (62/93) ======================= q-point: (-0.85 -0.15 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.622 ( -4.096 -0.716 0.000) 4.158 0.761 ( -0.474 0.811 0.000) 0.939 0.820 ( -7.809 5.016 0.000) 9.281 0.900 ( -4.614 0.289 0.000) 4.623 1.034 ( -1.633 -0.987 0.000) 1.908 1.115 ( -7.436 -2.753 0.000) 7.929 1.535 ( -7.362 8.668 0.000) 11.372 2.194 ( 4.629 -4.640 -0.000) 6.554 2.260 ( 2.607 2.017 -0.000) 3.296 2.996 ( -0.310 -0.575 -0.000) 0.653 3.763 ( -1.370 2.394 0.000) 2.759 5.380 ( 0.025 0.223 0.000) 0.224 ======================= Grid point 142 (63/93) ======================= q-point: ( 0.65 -0.15 -0.35) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 210 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.600 ( 0.000 -4.047 -0.000) 4.047 0.702 ( 0.000 -7.882 0.000) 7.882 0.758 ( 0.000 -0.478 -0.000) 0.478 0.944 ( 0.000 -5.246 0.000) 5.246 1.169 ( 0.000 -3.649 -0.000) 3.649 1.273 ( 0.000 -9.513 0.000) 9.513 1.402 ( 0.000 -7.228 0.000) 7.228 2.207 ( -0.000 4.064 0.000) 4.064 2.228 ( 0.000 3.793 -0.000) 3.793 2.716 ( 0.000 -1.533 0.000) 1.533 3.979 ( 0.000 -0.340 -0.000) 0.340 5.339 ( 0.000 0.129 -0.000) 0.129 ======================= Grid point 144 (64/93) ======================= q-point: ( 0.65 -0.15 -0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 227 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.452 ( 0.750 -0.750 4.264) 4.394 0.648 ( 6.390 -6.390 -1.356) 9.138 0.742 ( -0.262 0.262 1.922) 1.957 0.866 ( 4.026 -4.026 -2.964) 6.419 1.061 ( -1.154 1.154 -2.066) 2.633 1.199 ( 3.372 -3.372 2.801) 5.530 1.356 ( 6.753 -6.753 0.527) 9.565 2.218 ( -3.884 3.884 1.834) 5.791 2.308 ( -0.266 0.266 -1.288) 1.342 2.694 ( 0.875 -0.875 0.058) 1.239 3.975 ( 0.226 -0.226 0.485) 0.581 5.340 ( -0.050 0.050 0.101) 0.123 ======================= Grid point 153 (65/93) ======================= q-point: (-0.58 -0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.714 ( 0.000 -1.310 -1.202) 1.778 0.779 ( 0.000 -0.088 -0.152) 0.175 0.930 ( 0.000 -0.954 -2.563) 2.735 1.109 ( 0.000 -1.655 -2.924) 3.360 1.265 ( -0.000 -2.223 8.110) 8.409 1.713 ( 0.000 -8.797 -0.680) 8.824 1.838 ( 0.000 -2.817 0.163) 2.822 1.923 ( -0.000 -0.311 0.726) 0.790 1.978 ( -0.000 4.864 0.021) 4.864 2.884 ( 0.000 -0.179 -8.896) 8.897 3.921 ( -0.000 2.180 7.253) 7.574 5.362 ( 0.000 0.724 -1.426) 1.600 ======================= Grid point 156 (66/93) ======================= q-point: (-0.58 -0.18 0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.655 ( -1.574 -3.386 -2.339) 4.406 0.768 ( 0.197 -0.739 -0.956) 1.224 0.875 ( 0.169 -3.227 -4.507) 5.546 1.029 ( -0.997 -3.999 -1.282) 4.316 1.169 ( -3.886 -4.115 3.859) 6.850 1.411 ( -5.547 -7.561 2.756) 9.774 1.807 ( 6.016 -7.525 0.041) 9.634 1.904 ( -3.857 3.855 -0.029) 5.454 2.137 ( 4.529 2.781 -0.095) 5.316 2.874 ( -0.339 -0.383 -8.853) 8.867 3.959 ( -0.626 1.909 5.886) 6.219 5.378 ( -0.192 0.904 -1.649) 1.890 ======================= Grid point 157 (67/93) ======================= q-point: (-0.58 -0.18 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.556 ( 1.079 -1.079 -5.175) 5.395 0.784 ( 0.349 -0.349 -2.201) 2.256 0.834 ( 2.399 -2.399 -3.660) 4.991 0.964 ( 1.879 -1.879 0.276) 2.672 1.008 ( 0.128 -0.128 1.571) 1.582 1.299 ( 0.743 -0.743 4.743) 4.858 1.785 ( 7.917 -7.917 -0.203) 11.198 1.910 ( -4.407 4.407 -0.388) 6.244 2.224 ( 1.078 -1.078 -0.252) 1.545 2.863 ( 0.017 -0.017 -8.864) 8.864 3.972 ( -1.336 1.336 5.393) 5.715 5.386 ( -0.637 0.637 -1.739) 1.959 ======================= Grid point 173 (68/93) ======================= q-point: ( 0.20 -0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 158 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.684 ( -0.000 2.934 0.000) 2.934 0.748 ( -0.000 1.397 0.000) 1.397 0.807 ( -0.000 6.059 0.000) 6.059 1.060 ( -0.000 6.068 0.000) 6.068 1.108 ( -0.000 4.321 0.000) 4.321 1.694 ( -0.000 8.946 0.000) 8.946 1.933 ( 0.000 -2.694 -0.000) 2.694 1.994 ( 0.000 -2.427 -0.000) 2.427 1.994 ( 0.000 -4.729 -0.000) 4.729 3.038 ( 0.000 -0.848 -0.000) 0.848 3.649 ( -0.000 4.090 0.000) 4.090 5.363 ( 0.000 -0.006 -0.000) 0.006 ======================= Grid point 174 (69/93) ======================= q-point: (-0.30 -0.20 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.446 ( 0.000 -4.871 -5.094) 7.048 0.464 ( 0.000 -6.647 -6.460) 9.269 0.693 ( 0.000 -1.072 -1.267) 1.660 0.747 ( 0.000 -5.399 -1.826) 5.699 1.038 ( -0.000 -5.935 1.894) 6.230 1.044 ( -0.000 -7.404 1.356) 7.527 1.311 ( 0.000 -6.716 -11.887) 13.653 2.158 ( -0.000 3.323 0.719) 3.400 2.180 ( -0.000 1.478 0.333) 1.516 2.745 ( -0.000 -0.665 2.113) 2.216 4.096 ( -0.000 1.444 0.638) 1.579 5.380 ( 0.000 0.732 -1.606) 1.765 ======================= Grid point 176 (70/93) ======================= q-point: ( 0.20 -0.20 0.30) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.717 ( -1.065 1.623 0.000) 1.941 0.769 ( -0.262 0.712 0.000) 0.759 0.912 ( -0.030 3.979 0.000) 3.979 1.093 ( -4.216 0.886 0.000) 4.308 1.146 ( -1.984 1.156 0.000) 2.297 1.706 ( -8.657 5.517 0.000) 10.266 1.722 ( -9.674 -6.807 0.000) 11.829 2.010 ( 5.624 2.056 -0.000) 5.987 2.061 ( 6.256 -3.612 -0.000) 7.224 3.022 ( -0.093 -0.637 -0.000) 0.644 3.730 ( -0.274 3.770 0.000) 3.780 5.368 ( 0.131 0.352 0.000) 0.376 ======================= Grid point 177 (71/93) ======================= q-point: (-0.80 -0.20 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 256 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.692 ( -2.392 -1.194 0.000) 2.673 0.772 ( -0.374 0.596 0.000) 0.704 0.959 ( -1.408 1.124 0.000) 1.802 0.970 ( -6.301 0.556 0.000) 6.326 1.088 ( -3.135 -0.966 0.000) 3.280 1.322 ( -11.439 -7.321 0.000) 13.581 1.679 ( -6.439 7.028 0.000) 9.532 2.089 ( 5.447 -5.328 -0.000) 7.620 2.182 ( 4.719 2.925 -0.000) 5.552 3.004 ( -0.400 -0.648 -0.000) 0.762 3.789 ( -1.034 3.025 0.000) 3.197 5.380 ( 0.017 0.502 0.000) 0.502 ======================= Grid point 181 (72/93) ======================= q-point: ( 0.20 0.80 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.622 ( -0.894 0.894 0.000) 1.265 0.779 ( -1.012 1.012 0.000) 1.431 0.915 ( -3.655 3.655 0.000) 5.170 0.916 ( -1.382 1.382 0.000) 1.954 0.994 ( 0.491 -0.491 -0.000) 0.695 1.143 ( -3.148 3.148 0.000) 4.452 1.694 ( -6.466 6.466 0.000) 9.144 2.086 ( 5.673 -5.673 -0.000) 8.023 2.274 ( 0.821 -0.821 -0.000) 1.160 2.989 ( 0.051 -0.051 -0.000) 0.072 3.810 ( -2.038 2.038 0.000) 2.882 5.386 ( -0.303 0.303 0.000) 0.429 ======================= Grid point 186 (73/93) ======================= q-point: ( 0.70 -0.20 -0.30) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 210 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.499 ( 0.000 -5.316 -0.000) 5.316 0.523 ( -0.000 -8.008 0.000) 8.008 0.752 ( 0.000 -0.165 -0.000) 0.165 0.816 ( 0.000 -5.782 0.000) 5.782 1.113 ( 0.000 -1.806 0.000) 1.806 1.132 ( 0.000 -3.931 0.000) 3.931 1.276 ( 0.000 -4.390 0.000) 4.390 2.281 ( -0.000 2.726 0.000) 2.726 2.290 ( 0.000 2.030 -0.000) 2.030 2.685 ( 0.000 -1.291 0.000) 1.291 3.973 ( 0.000 -0.221 -0.000) 0.221 5.341 ( 0.000 0.083 -0.000) 0.083 ======================= Grid point 199 (74/93) ======================= q-point: (-0.53 -0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.671 ( 0.000 -2.748 -1.832) 3.302 0.768 ( 0.000 -0.957 -0.389) 1.033 0.872 ( 0.000 -4.582 -5.230) 6.954 1.039 ( 0.000 -5.021 -1.546) 5.253 1.208 ( -0.000 -2.962 5.278) 6.052 1.500 ( 0.000 -11.370 -0.518) 11.381 1.714 ( -0.000 -8.168 1.711) 8.345 1.967 ( -0.000 3.824 0.131) 3.826 2.068 ( -0.000 3.755 -0.119) 3.757 2.878 ( 0.000 -0.447 -8.782) 8.793 3.965 ( -0.000 1.996 5.787) 6.121 5.380 ( 0.000 1.017 -1.665) 1.951 ======================= Grid point 200 (75/93) ======================= q-point: (-0.53 -0.23 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.572 ( -2.135 -4.111 -3.757) 5.964 0.747 ( 2.063 -1.542 -2.662) 3.704 0.792 ( 2.116 -4.238 -3.565) 5.929 0.927 ( 1.365 -5.444 -0.044) 5.613 1.057 ( -4.289 -5.993 1.524) 7.526 1.280 ( 1.390 -5.186 2.205) 5.805 1.619 ( 7.078 -9.281 0.984) 11.714 2.001 ( -4.144 4.842 -0.308) 6.381 2.183 ( 2.706 1.536 -0.240) 3.121 2.864 ( -0.039 -0.611 -8.737) 8.759 3.995 ( -0.753 1.473 4.985) 5.252 5.398 ( -0.419 0.932 -1.807) 2.076 ======================= Grid point 218 (76/93) ======================= q-point: (-0.25 -0.25 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 93 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.383 ( 0.000 0.000 -7.062) 7.062 0.383 ( 0.000 0.000 -7.062) 7.062 0.681 ( 0.000 0.000 -1.186) 1.186 0.681 ( 0.000 0.000 -1.186) 1.186 0.964 ( -0.000 -0.000 2.103) 2.103 0.964 ( -0.000 -0.000 2.103) 2.103 1.243 ( 0.000 0.000 -16.087) 16.087 2.195 ( -0.000 -0.000 0.367) 0.367 2.195 ( -0.000 -0.000 0.367) 0.367 2.737 ( -0.000 -0.000 2.598) 2.598 4.112 ( -0.000 -0.000 1.863) 1.863 5.389 ( 0.000 0.000 -1.553) 1.553 ======================= Grid point 219 (77/93) ======================= q-point: ( 0.25 -0.25 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.727 ( -0.000 1.156 0.000) 1.156 0.772 ( -0.000 0.759 0.000) 0.759 0.912 ( -0.000 4.036 0.000) 4.036 1.138 ( -0.000 1.289 0.000) 1.289 1.167 ( -0.000 1.361 0.000) 1.361 1.857 ( -0.000 6.699 0.000) 6.699 1.879 ( 0.000 -6.365 -0.000) 6.365 1.893 ( 0.000 -0.930 -0.000) 0.930 1.940 ( 0.000 -3.044 -0.000) 3.044 3.023 ( 0.000 -0.581 -0.000) 0.581 3.733 ( -0.000 3.927 0.000) 3.927 5.367 ( -0.000 0.390 0.000) 0.390 ======================= Grid point 220 (78/93) ======================= q-point: (-0.75 -0.25 0.35) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.723 ( -0.777 -0.935 -0.000) 1.215 0.777 ( -0.180 0.163 0.000) 0.243 0.974 ( -0.281 1.832 0.000) 1.854 1.075 ( -3.555 -2.525 0.000) 4.360 1.147 ( -2.146 -0.853 0.000) 2.309 1.561 ( -10.730 -8.550 0.000) 13.720 1.807 ( -5.811 4.473 0.000) 7.333 1.964 ( 6.512 -5.940 -0.000) 8.814 2.071 ( 5.591 3.368 -0.000) 6.527 3.011 ( -0.228 -0.520 -0.000) 0.568 3.806 ( -0.556 3.404 0.000) 3.449 5.379 ( 0.004 0.735 0.000) 0.735 ======================= Grid point 224 (79/93) ======================= q-point: (-0.75 -0.25 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.651 ( -1.627 -2.186 -0.000) 2.726 0.790 ( 0.220 1.141 0.000) 1.162 0.951 ( -1.597 -1.673 -0.000) 2.313 0.970 ( -1.726 -0.800 0.000) 1.902 1.038 ( -4.190 -3.512 0.000) 5.468 1.229 ( -5.312 -1.461 0.000) 5.509 1.810 ( -4.804 5.514 0.000) 7.313 1.957 ( 6.646 -7.276 -0.000) 9.855 2.227 ( 3.492 1.217 -0.000) 3.698 2.992 ( -0.316 -0.422 -0.000) 0.528 3.847 ( -1.487 2.530 0.000) 2.935 5.392 ( -0.313 0.687 0.000) 0.755 ======================= Grid point 233 (80/93) ======================= q-point: ( 0.75 -0.25 -0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 69 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.425 ( -0.000 -0.000 -2.317) 2.317 0.425 ( 0.000 0.000 2.317) 2.317 0.750 ( -0.000 -0.000 -2.725) 2.725 0.750 ( 0.000 0.000 2.725) 2.725 1.093 ( -0.000 -0.000 -2.505) 2.505 1.093 ( 0.000 0.000 2.505) 2.505 1.229 ( 0.000 0.000 -0.000) 0.000 2.311 ( 0.000 0.000 -1.544) 1.544 2.311 ( 0.000 0.000 1.544) 1.544 2.670 ( 0.000 0.000 0.000) 0.000 3.971 ( -0.000 -0.000 -0.000) 0.000 5.342 ( 0.000 0.000 -0.000) 0.000 ======================= Grid point 245 (81/93) ======================= q-point: (-0.48 -0.28 0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.604 ( 0.000 -3.479 -2.193) 4.113 0.742 ( 0.000 -7.603 -5.275) 9.254 0.743 ( 0.000 -1.306 -0.700) 1.482 0.913 ( 0.000 -6.554 -0.784) 6.601 1.131 ( -0.000 -5.011 1.852) 5.343 1.248 ( 0.000 -12.256 -0.055) 12.256 1.534 ( -0.000 -7.868 2.044) 8.129 2.063 ( -0.000 4.923 -0.257) 4.929 2.136 ( -0.000 2.733 -0.218) 2.741 2.864 ( 0.000 -0.859 -8.694) 8.736 4.003 ( -0.000 1.514 4.873) 5.103 5.403 ( 0.000 1.049 -1.820) 2.101 ======================= Grid point 248 (82/93) ======================= q-point: (-0.48 -0.28 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.514 ( 0.720 -0.720 -4.430) 4.545 0.665 ( 5.369 -5.369 -3.742) 8.465 0.742 ( 1.381 -1.381 -1.825) 2.673 0.819 ( 4.146 -4.146 0.278) 5.870 0.967 ( 1.632 -1.632 1.637) 2.829 1.155 ( 7.266 -7.266 0.565) 10.292 1.465 ( 4.572 -4.572 0.955) 6.536 2.097 ( -4.085 4.085 -0.373) 5.790 2.198 ( 0.203 -0.203 -0.289) 0.407 2.851 ( 0.535 -0.535 -8.694) 8.727 4.019 ( -0.806 0.806 4.529) 4.670 5.415 ( -0.611 0.611 -1.845) 2.037 ======================= Grid point 265 (83/93) ======================= q-point: (-0.70 -0.30 0.30) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 158 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.731 ( 0.000 -0.823 -0.000) 0.823 0.779 ( 0.000 -0.042 -0.000) 0.042 0.976 ( -0.000 2.046 0.000) 2.046 1.113 ( 0.000 -3.559 -0.000) 3.559 1.170 ( 0.000 -0.823 -0.000) 0.823 1.721 ( 0.000 -8.869 -0.000) 8.869 1.850 ( 0.000 -5.911 -0.000) 5.911 1.901 ( -0.000 1.812 0.000) 1.812 1.978 ( -0.000 4.983 0.000) 4.983 3.013 ( 0.000 -0.447 -0.000) 0.447 3.812 ( -0.000 3.531 0.000) 3.531 5.379 ( -0.000 0.827 0.000) 0.827 ======================= Grid point 269 (84/93) ======================= q-point: (-0.70 -0.30 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.682 ( -1.118 -2.884 -0.000) 3.093 0.778 ( 0.644 -0.056 -0.000) 0.647 0.976 ( -0.774 -1.906 -0.000) 2.057 0.992 ( -0.453 -5.304 -0.000) 5.323 1.119 ( -2.762 -1.874 0.000) 3.338 1.376 ( -8.314 -8.470 -0.000) 11.869 1.808 ( 7.473 -8.943 -0.000) 11.654 1.903 ( -4.198 4.788 0.000) 6.367 2.139 ( 4.570 2.892 -0.000) 5.408 2.999 ( -0.268 -0.607 -0.000) 0.663 3.872 ( -0.788 2.839 0.000) 2.946 5.398 ( -0.204 0.997 0.000) 1.018 ======================= Grid point 271 (85/93) ======================= q-point: ( 0.30 0.70 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 143 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.615 ( 1.036 -1.036 -0.000) 1.466 0.807 ( 0.150 -0.150 -0.000) 0.211 0.910 ( 1.507 -1.507 -0.000) 2.131 0.932 ( 2.854 -2.854 -0.000) 4.036 0.988 ( 0.229 -0.229 -0.000) 0.324 1.230 ( -0.843 0.843 0.000) 1.192 1.787 ( 8.811 -8.811 -0.000) 12.460 1.915 ( -4.618 4.618 0.000) 6.531 2.228 ( 1.063 -1.063 -0.000) 1.503 2.986 ( 0.161 -0.161 -0.000) 0.228 3.894 ( -1.840 1.840 0.000) 2.602 5.407 ( -0.684 0.684 0.000) 0.968 ======================= Grid point 293 (86/93) ======================= q-point: (-0.43 -0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 220 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.534 ( 0.000 -2.943 -2.923) 4.148 0.580 ( 0.000 -6.881 -4.139) 8.030 0.717 ( 0.000 -1.040 -0.924) 1.391 0.779 ( 0.000 -5.885 -0.939) 5.960 0.998 ( -0.000 -6.513 1.310) 6.644 1.018 ( -0.000 -8.723 0.905) 8.770 1.422 ( -0.000 -2.808 0.296) 2.824 2.156 ( 0.000 3.577 -0.362) 3.595 2.181 ( 0.000 1.485 -0.289) 1.512 2.845 ( 0.000 -0.878 -8.695) 8.740 4.027 ( -0.000 0.818 4.396) 4.472 5.421 ( 0.000 0.687 -1.834) 1.959 ======================= Grid point 313 (87/93) ======================= q-point: (-0.65 -0.35 0.35) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.694 ( 0.000 -2.577 -0.000) 2.577 0.770 ( 0.000 -0.753 -0.000) 0.753 0.984 ( 0.000 -1.728 -0.000) 1.728 0.995 ( 0.000 -7.431 -0.000) 7.431 1.148 ( 0.000 -1.041 -0.000) 1.041 1.507 ( 0.000 -11.521 -0.000) 11.521 1.690 ( 0.000 -9.144 -0.000) 9.144 1.966 ( -0.000 4.238 0.000) 4.238 2.070 ( -0.000 3.826 0.000) 3.826 3.002 ( 0.000 -0.648 -0.000) 0.648 3.880 ( -0.000 2.945 0.000) 2.945 5.400 ( -0.000 1.117 0.000) 1.117 ======================= Grid point 315 (88/93) ======================= q-point: (-0.65 -0.35 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 250 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.613 ( -0.811 -3.274 -0.000) 3.373 0.775 ( 2.079 -0.205 -0.000) 2.089 0.861 ( 3.501 -6.970 -0.000) 7.800 0.902 ( -0.268 -4.704 -0.000) 4.711 1.039 ( -4.472 -5.748 -0.000) 7.283 1.247 ( -0.535 -3.773 -0.000) 3.811 1.602 ( 8.175 -9.978 -0.000) 12.899 2.006 ( -4.101 4.937 0.000) 6.418 2.187 ( 2.732 1.591 -0.000) 3.161 2.985 ( 0.015 -0.741 -0.000) 0.741 3.923 ( -0.975 2.073 0.000) 2.291 5.419 ( -0.433 0.976 0.000) 1.068 ======================= Grid point 340 (89/93) ======================= q-point: (-0.38 -0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 72 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.499 ( 0.000 0.000 -3.692) 3.692 0.499 ( 0.000 0.000 -3.692) 3.692 0.705 ( 0.000 0.000 -0.956) 0.956 0.705 ( 0.000 0.000 -0.956) 0.956 0.921 ( -0.000 -0.000 1.631) 1.631 0.921 ( -0.000 -0.000 1.631) 1.631 1.397 ( 0.000 0.000 -0.847) 0.847 2.196 ( 0.000 0.000 -0.315) 0.315 2.196 ( 0.000 0.000 -0.315) 0.315 2.834 ( 0.000 0.000 -8.741) 8.741 4.036 ( -0.000 -0.000 4.253) 4.253 5.429 ( 0.000 0.000 -1.808) 1.808 ======================= Grid point 361 (90/93) ======================= q-point: (-0.60 -0.40 0.40) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 158 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.629 ( 0.000 -3.351 -0.000) 3.351 0.750 ( 0.000 -1.104 -0.000) 1.104 0.815 ( 0.000 -9.387 -0.000) 9.387 0.906 ( 0.000 -5.229 -0.000) 5.229 1.110 ( 0.000 -3.631 -0.000) 3.631 1.249 ( 0.000 -12.630 -0.000) 12.630 1.498 ( 0.000 -7.788 -0.000) 7.788 2.067 ( -0.000 5.062 0.000) 5.062 2.139 ( -0.000 2.793 0.000) 2.793 2.985 ( 0.000 -0.990 -0.000) 0.990 3.934 ( -0.000 2.155 0.000) 2.155 5.424 ( -0.000 1.087 0.000) 1.087 ======================= Grid point 365 (91/93) ======================= q-point: ( 0.40 0.60 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.563 ( 1.121 -1.121 -0.000) 1.585 0.715 ( 6.439 -6.439 -0.000) 9.106 0.763 ( 1.661 -1.661 -0.000) 2.349 0.806 ( 3.029 -3.029 -0.000) 4.284 0.947 ( 1.584 -1.584 -0.000) 2.241 1.149 ( 6.837 -6.837 -0.000) 9.669 1.440 ( 4.373 -4.373 -0.000) 6.185 2.103 ( -4.073 4.073 0.000) 5.760 2.202 ( 0.191 -0.191 -0.000) 0.269 2.970 ( 0.578 -0.578 -0.000) 0.817 3.957 ( -1.104 1.104 0.000) 1.561 5.436 ( -0.609 0.609 0.000) 0.861 ======================= Grid point 411 (92/93) ======================= q-point: (-0.55 -0.45 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.565 ( 0.000 -2.458 -0.000) 2.458 0.627 ( 0.000 -7.662 -0.000) 7.662 0.728 ( 0.000 -0.946 -0.000) 0.946 0.787 ( 0.000 -5.545 -0.000) 5.545 0.984 ( 0.000 -6.850 -0.000) 6.850 1.008 ( 0.000 -9.349 -0.000) 9.349 1.403 ( 0.000 -1.700 -0.000) 1.700 2.162 ( -0.000 3.548 0.000) 3.548 2.185 ( -0.000 1.522 0.000) 1.522 2.964 ( 0.000 -0.901 -0.000) 0.901 3.969 ( -0.000 1.151 0.000) 1.151 5.443 ( -0.000 0.675 0.000) 0.675 ======================= Grid point 460 (93/93) ======================= q-point: ( 0.50 0.50 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 6.53e-05 6.53e-05 Number of triplets: 63 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.538 ( 0.000 0.000 0.000) 0.000 0.538 ( 0.000 -0.000 0.000) 0.000 0.716 ( 0.000 -0.000 0.000) 0.000 0.716 ( 0.000 -0.000 0.000) 0.000 0.903 ( 0.000 0.000 0.000) 0.000 0.903 ( 0.000 -0.000 0.000) 0.000 1.391 ( 0.000 -0.000 0.000) 0.000 2.201 ( 0.000 -0.000 0.000) 0.000 2.201 ( 0.000 -0.000 0.000) 0.000 2.953 ( 0.000 -0.000 0.000) 0.000 3.981 ( 0.000 0.000 0.000) 0.000 5.450 ( 0.000 0.000 0.000) 0.000 =================== 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/9600 10.0 12.248 12.248 8.909 -0.000 -0.000 -0.000 3/9600 20.0 4.270 4.270 4.163 -0.000 -0.000 -0.000 3/9600 30.0 2.628 2.628 2.736 -0.000 -0.000 -0.000 3/9600 40.0 1.922 1.922 2.075 -0.000 -0.000 -0.000 3/9600 50.0 1.523 1.523 1.696 -0.000 -0.000 -0.000 3/9600 60.0 1.263 1.263 1.444 -0.000 -0.000 -0.000 3/9600 70.0 1.080 1.080 1.261 -0.000 -0.000 -0.000 3/9600 80.0 0.944 0.944 1.120 -0.000 -0.000 -0.000 3/9600 90.0 0.839 0.839 1.007 -0.000 -0.000 -0.000 3/9600 100.0 0.754 0.754 0.915 -0.000 -0.000 -0.000 3/9600 110.0 0.686 0.686 0.839 -0.000 -0.000 -0.000 3/9600 120.0 0.628 0.628 0.773 -0.000 -0.000 -0.000 3/9600 130.0 0.580 0.580 0.718 -0.000 -0.000 -0.000 3/9600 140.0 0.539 0.539 0.669 -0.000 -0.000 -0.000 3/9600 150.0 0.503 0.503 0.627 -0.000 -0.000 -0.000 3/9600 160.0 0.471 0.471 0.589 -0.000 -0.000 -0.000 3/9600 170.0 0.444 0.444 0.556 -0.000 -0.000 -0.000 3/9600 180.0 0.419 0.419 0.527 -0.000 -0.000 -0.000 3/9600 190.0 0.397 0.397 0.500 -0.000 -0.000 -0.000 3/9600 200.0 0.377 0.377 0.476 -0.000 -0.000 -0.000 3/9600 210.0 0.359 0.359 0.454 -0.000 -0.000 -0.000 3/9600 220.0 0.343 0.343 0.434 -0.000 -0.000 -0.000 3/9600 230.0 0.328 0.328 0.415 -0.000 -0.000 -0.000 3/9600 240.0 0.314 0.314 0.398 -0.000 -0.000 -0.000 3/9600 250.0 0.302 0.302 0.383 -0.000 -0.000 -0.000 3/9600 260.0 0.290 0.290 0.368 -0.000 -0.000 -0.000 3/9600 270.0 0.279 0.279 0.355 -0.000 -0.000 -0.000 3/9600 280.0 0.269 0.269 0.343 -0.000 -0.000 -0.000 3/9600 290.0 0.260 0.260 0.331 -0.000 -0.000 -0.000 3/9600 300.0 0.251 0.251 0.320 -0.000 -0.000 -0.000 3/9600 310.0 0.243 0.243 0.310 -0.000 -0.000 -0.000 3/9600 320.0 0.236 0.236 0.301 -0.000 -0.000 -0.000 3/9600 330.0 0.229 0.229 0.292 -0.000 -0.000 -0.000 3/9600 340.0 0.222 0.222 0.283 -0.000 -0.000 -0.000 3/9600 350.0 0.216 0.216 0.275 -0.000 -0.000 -0.000 3/9600 360.0 0.210 0.210 0.268 -0.000 -0.000 -0.000 3/9600 370.0 0.204 0.204 0.260 -0.000 -0.000 -0.000 3/9600 380.0 0.199 0.199 0.254 -0.000 -0.000 -0.000 3/9600 390.0 0.193 0.193 0.247 -0.000 -0.000 -0.000 3/9600 400.0 0.189 0.189 0.241 -0.000 -0.000 -0.000 3/9600 410.0 0.184 0.184 0.235 -0.000 -0.000 -0.000 3/9600 420.0 0.180 0.180 0.230 -0.000 -0.000 -0.000 3/9600 430.0 0.176 0.176 0.225 -0.000 -0.000 -0.000 3/9600 440.0 0.172 0.172 0.219 -0.000 -0.000 -0.000 3/9600 450.0 0.168 0.168 0.215 -0.000 -0.000 -0.000 3/9600 460.0 0.164 0.164 0.210 -0.000 -0.000 -0.000 3/9600 470.0 0.161 0.161 0.206 -0.000 -0.000 -0.000 3/9600 480.0 0.157 0.157 0.201 -0.000 -0.000 -0.000 3/9600 490.0 0.154 0.154 0.197 -0.000 -0.000 -0.000 3/9600 500.0 0.151 0.151 0.193 -0.000 -0.000 -0.000 3/9600 510.0 0.148 0.148 0.190 -0.000 -0.000 -0.000 3/9600 520.0 0.145 0.145 0.186 -0.000 -0.000 -0.000 3/9600 530.0 0.142 0.142 0.183 -0.000 -0.000 -0.000 3/9600 540.0 0.140 0.140 0.179 -0.000 -0.000 -0.000 3/9600 550.0 0.137 0.137 0.176 -0.000 -0.000 -0.000 3/9600 560.0 0.135 0.135 0.173 -0.000 -0.000 -0.000 3/9600 570.0 0.132 0.132 0.170 -0.000 -0.000 -0.000 3/9600 580.0 0.130 0.130 0.167 -0.000 -0.000 -0.000 3/9600 590.0 0.128 0.128 0.164 -0.000 -0.000 -0.000 3/9600 600.0 0.126 0.126 0.161 -0.000 -0.000 -0.000 3/9600 610.0 0.124 0.124 0.159 -0.000 -0.000 -0.000 3/9600 620.0 0.122 0.122 0.156 -0.000 -0.000 -0.000 3/9600 630.0 0.120 0.120 0.154 -0.000 -0.000 -0.000 3/9600 640.0 0.118 0.118 0.151 -0.000 -0.000 -0.000 3/9600 650.0 0.116 0.116 0.149 -0.000 -0.000 -0.000 3/9600 660.0 0.114 0.114 0.147 -0.000 -0.000 -0.000 3/9600 670.0 0.113 0.113 0.145 -0.000 -0.000 -0.000 3/9600 680.0 0.111 0.111 0.143 -0.000 -0.000 -0.000 3/9600 690.0 0.109 0.109 0.140 -0.000 -0.000 -0.000 3/9600 700.0 0.108 0.108 0.138 -0.000 -0.000 -0.000 3/9600 710.0 0.106 0.106 0.137 -0.000 -0.000 -0.000 3/9600 720.0 0.105 0.105 0.135 -0.000 -0.000 -0.000 3/9600 730.0 0.103 0.103 0.133 -0.000 -0.000 -0.000 3/9600 740.0 0.102 0.102 0.131 -0.000 -0.000 -0.000 3/9600 750.0 0.101 0.101 0.129 -0.000 -0.000 -0.000 3/9600 760.0 0.099 0.099 0.128 -0.000 -0.000 -0.000 3/9600 770.0 0.098 0.098 0.126 -0.000 -0.000 -0.000 3/9600 780.0 0.097 0.097 0.124 -0.000 -0.000 -0.000 3/9600 790.0 0.096 0.096 0.123 -0.000 -0.000 -0.000 3/9600 800.0 0.094 0.094 0.121 -0.000 -0.000 -0.000 3/9600 810.0 0.093 0.093 0.120 -0.000 -0.000 -0.000 3/9600 820.0 0.092 0.092 0.118 -0.000 -0.000 -0.000 3/9600 830.0 0.091 0.091 0.117 -0.000 -0.000 -0.000 3/9600 840.0 0.090 0.090 0.116 -0.000 -0.000 -0.000 3/9600 850.0 0.089 0.089 0.114 -0.000 -0.000 -0.000 3/9600 860.0 0.088 0.088 0.113 -0.000 -0.000 -0.000 3/9600 870.0 0.087 0.087 0.112 -0.000 -0.000 -0.000 3/9600 880.0 0.086 0.086 0.110 -0.000 -0.000 -0.000 3/9600 890.0 0.085 0.085 0.109 -0.000 -0.000 -0.000 3/9600 900.0 0.084 0.084 0.108 -0.000 -0.000 -0.000 3/9600 910.0 0.083 0.083 0.107 -0.000 -0.000 -0.000 3/9600 920.0 0.082 0.082 0.106 -0.000 -0.000 -0.000 3/9600 930.0 0.081 0.081 0.104 -0.000 -0.000 -0.000 3/9600 940.0 0.080 0.080 0.103 -0.000 -0.000 -0.000 3/9600 950.0 0.079 0.079 0.102 -0.000 -0.000 -0.000 3/9600 960.0 0.079 0.079 0.101 -0.000 -0.000 -0.000 3/9600 970.0 0.078 0.078 0.100 -0.000 -0.000 -0.000 3/9600 980.0 0.077 0.077 0.099 -0.000 -0.000 -0.000 3/9600 990.0 0.076 0.076 0.098 -0.000 -0.000 -0.000 3/9600 1000.0 0.076 0.076 0.097 -0.000 -0.000 -0.000 3/9600 Thermal conductivity related properties were written into "kappa-m21040.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 05:10:51]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|