----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 18:35:33]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [4 4 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) Number of symmetry operations in supercell: 384 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.123353369962359 0.000000000000000 0.000000000000000 b -2.061676684981180 3.570928767167578 0.000000000000000 c 0.000000000000000 0.000000000000000 6.710167980000000 Atomic positions (fractional): *1 S 0.66666666666667 0.33333333333333 0.87672146406862 32.065 2 S 0.33333333333333 0.66666666666667 0.37672146406862 32.065 *3 Cd 0.66666666666667 0.33333333333333 0.50027853593138 112.411 4 Cd 0.33333333333333 0.66666666666667 0.00027853593138 112.411 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.123353369962359 0.000000000000000 0.000000000000000 b -2.061676684981180 3.570928767167578 0.000000000000000 c 0.000000000000000 0.000000000000000 6.710167980000000 Atomic positions (fractional): *1 S 0.66666666666667 0.33333333333333 0.87672146406862 32.065 > 1 2 S 0.33333333333333 0.66666666666667 0.37672146406862 32.065 > 2 *3 Cd 0.66666666666667 0.33333333333333 0.50027853593138 112.411 > 3 4 Cd 0.33333333333333 0.66666666666667 0.00027853593138 112.411 > 4 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.493413479849437 0.000000000000000 0.000000000000000 b -8.246706739924720 14.283715068670311 0.000000000000000 c 0.000000000000000 0.000000000000000 13.420335959999999 Atomic positions (fractional): *1 S 0.16666666666667 0.08333333333333 0.43836073203431 32.065 > 1 2 S 0.41666666666667 0.08333333333333 0.43836073203431 32.065 > 1 3 S 0.66666666666667 0.08333333333333 0.43836073203431 32.065 > 1 4 S 0.91666666666667 0.08333333333333 0.43836073203431 32.065 > 1 5 S 0.16666666666667 0.33333333333333 0.43836073203431 32.065 > 1 6 S 0.41666666666667 0.33333333333333 0.43836073203431 32.065 > 1 7 S 0.66666666666667 0.33333333333333 0.43836073203431 32.065 > 1 8 S 0.91666666666667 0.33333333333333 0.43836073203431 32.065 > 1 9 S 0.16666666666667 0.58333333333333 0.43836073203431 32.065 > 1 10 S 0.41666666666667 0.58333333333333 0.43836073203431 32.065 > 1 11 S 0.66666666666667 0.58333333333333 0.43836073203431 32.065 > 1 12 S 0.91666666666667 0.58333333333333 0.43836073203431 32.065 > 1 13 S 0.16666666666667 0.83333333333333 0.43836073203431 32.065 > 1 14 S 0.41666666666667 0.83333333333333 0.43836073203431 32.065 > 1 15 S 0.66666666666667 0.83333333333333 0.43836073203431 32.065 > 1 16 S 0.91666666666667 0.83333333333333 0.43836073203431 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.93836073203431 32.065 > 1 18 S 0.41666666666667 0.08333333333333 0.93836073203431 32.065 > 1 19 S 0.66666666666667 0.08333333333333 0.93836073203431 32.065 > 1 20 S 0.91666666666667 0.08333333333333 0.93836073203431 32.065 > 1 21 S 0.16666666666667 0.33333333333333 0.93836073203431 32.065 > 1 22 S 0.41666666666667 0.33333333333333 0.93836073203431 32.065 > 1 23 S 0.66666666666667 0.33333333333333 0.93836073203431 32.065 > 1 24 S 0.91666666666667 0.33333333333333 0.93836073203431 32.065 > 1 25 S 0.16666666666667 0.58333333333333 0.93836073203431 32.065 > 1 26 S 0.41666666666667 0.58333333333333 0.93836073203431 32.065 > 1 27 S 0.66666666666667 0.58333333333333 0.93836073203431 32.065 > 1 28 S 0.91666666666667 0.58333333333333 0.93836073203431 32.065 > 1 29 S 0.16666666666667 0.83333333333333 0.93836073203431 32.065 > 1 30 S 0.41666666666667 0.83333333333333 0.93836073203431 32.065 > 1 31 S 0.66666666666667 0.83333333333333 0.93836073203431 32.065 > 1 32 S 0.91666666666667 0.83333333333333 0.93836073203431 32.065 > 1 33 S 0.08333333333333 0.16666666666667 0.18836073203431 32.065 > 2 34 S 0.33333333333333 0.16666666666667 0.18836073203431 32.065 > 2 35 S 0.58333333333333 0.16666666666667 0.18836073203431 32.065 > 2 36 S 0.83333333333333 0.16666666666667 0.18836073203431 32.065 > 2 37 S 0.08333333333333 0.41666666666667 0.18836073203431 32.065 > 2 38 S 0.33333333333333 0.41666666666667 0.18836073203431 32.065 > 2 39 S 0.58333333333333 0.41666666666667 0.18836073203431 32.065 > 2 40 S 0.83333333333333 0.41666666666667 0.18836073203431 32.065 > 2 41 S 0.08333333333333 0.66666666666667 0.18836073203431 32.065 > 2 42 S 0.33333333333333 0.66666666666667 0.18836073203431 32.065 > 2 43 S 0.58333333333333 0.66666666666667 0.18836073203431 32.065 > 2 44 S 0.83333333333333 0.66666666666667 0.18836073203431 32.065 > 2 45 S 0.08333333333333 0.91666666666667 0.18836073203431 32.065 > 2 46 S 0.33333333333333 0.91666666666667 0.18836073203431 32.065 > 2 47 S 0.58333333333333 0.91666666666667 0.18836073203431 32.065 > 2 48 S 0.83333333333333 0.91666666666667 0.18836073203431 32.065 > 2 49 S 0.08333333333333 0.16666666666667 0.68836073203431 32.065 > 2 50 S 0.33333333333333 0.16666666666667 0.68836073203431 32.065 > 2 51 S 0.58333333333333 0.16666666666667 0.68836073203431 32.065 > 2 52 S 0.83333333333333 0.16666666666667 0.68836073203431 32.065 > 2 53 S 0.08333333333333 0.41666666666667 0.68836073203431 32.065 > 2 54 S 0.33333333333333 0.41666666666667 0.68836073203431 32.065 > 2 55 S 0.58333333333333 0.41666666666667 0.68836073203431 32.065 > 2 56 S 0.83333333333333 0.41666666666667 0.68836073203431 32.065 > 2 57 S 0.08333333333333 0.66666666666667 0.68836073203431 32.065 > 2 58 S 0.33333333333333 0.66666666666667 0.68836073203431 32.065 > 2 59 S 0.58333333333333 0.66666666666667 0.68836073203431 32.065 > 2 60 S 0.83333333333333 0.66666666666667 0.68836073203431 32.065 > 2 61 S 0.08333333333333 0.91666666666667 0.68836073203431 32.065 > 2 62 S 0.33333333333333 0.91666666666667 0.68836073203431 32.065 > 2 63 S 0.58333333333333 0.91666666666667 0.68836073203431 32.065 > 2 64 S 0.83333333333333 0.91666666666667 0.68836073203431 32.065 > 2 *65 Cd 0.16666666666667 0.08333333333333 0.25013926796569 112.411 > 3 66 Cd 0.41666666666667 0.08333333333333 0.25013926796569 112.411 > 3 67 Cd 0.66666666666667 0.08333333333333 0.25013926796569 112.411 > 3 68 Cd 0.91666666666667 0.08333333333333 0.25013926796569 112.411 > 3 69 Cd 0.16666666666667 0.33333333333333 0.25013926796569 112.411 > 3 70 Cd 0.41666666666667 0.33333333333333 0.25013926796569 112.411 > 3 71 Cd 0.66666666666667 0.33333333333333 0.25013926796569 112.411 > 3 72 Cd 0.91666666666667 0.33333333333333 0.25013926796569 112.411 > 3 73 Cd 0.16666666666667 0.58333333333333 0.25013926796569 112.411 > 3 74 Cd 0.41666666666667 0.58333333333333 0.25013926796569 112.411 > 3 75 Cd 0.66666666666667 0.58333333333333 0.25013926796569 112.411 > 3 76 Cd 0.91666666666667 0.58333333333333 0.25013926796569 112.411 > 3 77 Cd 0.16666666666667 0.83333333333333 0.25013926796569 112.411 > 3 78 Cd 0.41666666666667 0.83333333333333 0.25013926796569 112.411 > 3 79 Cd 0.66666666666667 0.83333333333333 0.25013926796569 112.411 > 3 80 Cd 0.91666666666667 0.83333333333333 0.25013926796569 112.411 > 3 81 Cd 0.16666666666667 0.08333333333333 0.75013926796569 112.411 > 3 82 Cd 0.41666666666667 0.08333333333333 0.75013926796569 112.411 > 3 83 Cd 0.66666666666667 0.08333333333333 0.75013926796569 112.411 > 3 84 Cd 0.91666666666667 0.08333333333333 0.75013926796569 112.411 > 3 85 Cd 0.16666666666667 0.33333333333333 0.75013926796569 112.411 > 3 86 Cd 0.41666666666667 0.33333333333333 0.75013926796569 112.411 > 3 87 Cd 0.66666666666667 0.33333333333333 0.75013926796569 112.411 > 3 88 Cd 0.91666666666667 0.33333333333333 0.75013926796569 112.411 > 3 89 Cd 0.16666666666667 0.58333333333333 0.75013926796569 112.411 > 3 90 Cd 0.41666666666667 0.58333333333333 0.75013926796569 112.411 > 3 91 Cd 0.66666666666667 0.58333333333333 0.75013926796569 112.411 > 3 92 Cd 0.91666666666667 0.58333333333333 0.75013926796569 112.411 > 3 93 Cd 0.16666666666667 0.83333333333333 0.75013926796569 112.411 > 3 94 Cd 0.41666666666667 0.83333333333333 0.75013926796569 112.411 > 3 95 Cd 0.66666666666667 0.83333333333333 0.75013926796569 112.411 > 3 96 Cd 0.91666666666667 0.83333333333333 0.75013926796569 112.411 > 3 97 Cd 0.08333333333333 0.16666666666667 0.00013926796569 112.411 > 4 98 Cd 0.33333333333333 0.16666666666667 0.00013926796569 112.411 > 4 99 Cd 0.58333333333333 0.16666666666667 0.00013926796569 112.411 > 4 100 Cd 0.83333333333333 0.16666666666667 0.00013926796569 112.411 > 4 101 Cd 0.08333333333333 0.41666666666667 0.00013926796569 112.411 > 4 102 Cd 0.33333333333333 0.41666666666667 0.00013926796569 112.411 > 4 103 Cd 0.58333333333333 0.41666666666667 0.00013926796569 112.411 > 4 104 Cd 0.83333333333333 0.41666666666667 0.00013926796569 112.411 > 4 105 Cd 0.08333333333333 0.66666666666667 0.00013926796569 112.411 > 4 106 Cd 0.33333333333333 0.66666666666667 0.00013926796569 112.411 > 4 107 Cd 0.58333333333333 0.66666666666667 0.00013926796569 112.411 > 4 108 Cd 0.83333333333333 0.66666666666667 0.00013926796569 112.411 > 4 109 Cd 0.08333333333333 0.91666666666667 0.00013926796569 112.411 > 4 110 Cd 0.33333333333333 0.91666666666667 0.00013926796569 112.411 > 4 111 Cd 0.58333333333333 0.91666666666667 0.00013926796569 112.411 > 4 112 Cd 0.83333333333333 0.91666666666667 0.00013926796569 112.411 > 4 113 Cd 0.08333333333333 0.16666666666667 0.50013926796569 112.411 > 4 114 Cd 0.33333333333333 0.16666666666667 0.50013926796569 112.411 > 4 115 Cd 0.58333333333333 0.16666666666667 0.50013926796569 112.411 > 4 116 Cd 0.83333333333333 0.16666666666667 0.50013926796569 112.411 > 4 117 Cd 0.08333333333333 0.41666666666667 0.50013926796569 112.411 > 4 118 Cd 0.33333333333333 0.41666666666667 0.50013926796569 112.411 > 4 119 Cd 0.58333333333333 0.41666666666667 0.50013926796569 112.411 > 4 120 Cd 0.83333333333333 0.41666666666667 0.50013926796569 112.411 > 4 121 Cd 0.08333333333333 0.66666666666667 0.50013926796569 112.411 > 4 122 Cd 0.33333333333333 0.66666666666667 0.50013926796569 112.411 > 4 123 Cd 0.58333333333333 0.66666666666667 0.50013926796569 112.411 > 4 124 Cd 0.83333333333333 0.66666666666667 0.50013926796569 112.411 > 4 125 Cd 0.08333333333333 0.91666666666667 0.50013926796569 112.411 > 4 126 Cd 0.33333333333333 0.91666666666667 0.50013926796569 112.411 > 4 127 Cd 0.58333333333333 0.91666666666667 0.50013926796569 112.411 > 4 128 Cd 0.83333333333333 0.91666666666667 0.50013926796569 112.411 > 4 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.4947083 -0.0000000 0.0000000 0.0000000 7.4947083 0.0000000 0.0000000 0.0000000 7.5360942 -------------------------- Born effective charges -------------------------- 1 S -2.1930538 -0.0000000 0.0000000 0.0000000 -2.1930538 0.0000000 0.0000000 0.0000000 -2.2856256 2 S -2.1930538 -0.0000000 0.0000000 0.0000000 -2.1930538 0.0000000 0.0000000 0.0000000 -2.2856256 3 Cd 2.1930538 0.0000000 0.0000000 0.0000000 2.1930538 0.0000000 0.0000000 0.0000000 2.2856256 4 Cd 2.1930538 0.0000000 0.0000000 0.0000000 2.1930538 0.0000000 0.0000000 0.0000000 2.2856256 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 384/384 Permutation basis: 2856/2856 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 224 Number of blocks in projector: 188 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 11 --- Eigsh_solver_block: 1 / 11 --- Block_size: 98 Use standard eigh solver. --- Eigsh_solver_block: 2 / 11 --- Block_size: 72 Use standard eigh solver. --- Eigsh_solver_block: 3 / 11 --- Block_size: 46 Use standard eigh solver. --- Eigsh_solver_block: 4 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 5 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 6 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 7 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 8 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 9 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 10 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 11 / 11 --- Block_size: 1 Use standard eigh solver. Tree of FC basis block matrices: - (224, 220), data: False |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (46, 44), data: True |-- (72, 72), data: True |-- (98, 96), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.022 Solver_block: 80 / 80 - Time: 0.186 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.241 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (xx) -0.00000000 (xx) Permutation basis: 384/384 Permutation basis: 2856/2856 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 224 Number of blocks in projector: 188 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 11 --- Eigsh_solver_block: 1 / 11 --- Block_size: 98 Use standard eigh solver. --- Eigsh_solver_block: 2 / 11 --- Block_size: 72 Use standard eigh solver. --- Eigsh_solver_block: 3 / 11 --- Block_size: 46 Use standard eigh solver. --- Eigsh_solver_block: 4 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 5 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 6 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 7 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 8 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 9 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 10 / 11 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 11 / 11 --- Block_size: 1 Use standard eigh solver. Tree of FC basis block matrices: - (224, 220), data: False |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (1, 1), data: True |-- (46, 44), data: True |-- (72, 72), data: True |-- (98, 96), data: True Max drift after symmetrization by symfc projector: -0.00000000 (xx) -0.00000000 (xx) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 18:35:38]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:35: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): [4 4 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.123353369962359 0.000000000000000 0.000000000000000 b -2.061676684981180 3.570928767167578 0.000000000000000 c 0.000000000000000 0.000000000000000 6.710167980000000 Atomic positions (fractional): 1 S 0.66666666666667 0.33333333333333 0.87672146406862 32.065 2 S 0.33333333333333 0.66666666666667 0.37672146406862 32.065 3 Cd 0.66666666666667 0.33333333333333 0.50027853593138 112.411 4 Cd 0.33333333333333 0.66666666666667 0.00027853593138 112.411 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.493413479849437 0.000000000000000 0.000000000000000 b -8.246706739924720 14.283715068670311 0.000000000000000 c 0.000000000000000 0.000000000000000 13.420335959999999 Atomic positions (fractional): 1 S 0.16666666666667 0.08333333333333 0.43836073203431 32.065 > 1 2 S 0.41666666666667 0.08333333333333 0.43836073203431 32.065 > 1 3 S 0.66666666666667 0.08333333333333 0.43836073203431 32.065 > 1 4 S 0.91666666666667 0.08333333333333 0.43836073203431 32.065 > 1 5 S 0.16666666666667 0.33333333333333 0.43836073203431 32.065 > 1 6 S 0.41666666666667 0.33333333333333 0.43836073203431 32.065 > 1 7 S 0.66666666666667 0.33333333333333 0.43836073203431 32.065 > 1 8 S 0.91666666666667 0.33333333333333 0.43836073203431 32.065 > 1 9 S 0.16666666666667 0.58333333333333 0.43836073203431 32.065 > 1 10 S 0.41666666666667 0.58333333333333 0.43836073203431 32.065 > 1 11 S 0.66666666666667 0.58333333333333 0.43836073203431 32.065 > 1 12 S 0.91666666666667 0.58333333333333 0.43836073203431 32.065 > 1 13 S 0.16666666666667 0.83333333333333 0.43836073203431 32.065 > 1 14 S 0.41666666666667 0.83333333333333 0.43836073203431 32.065 > 1 15 S 0.66666666666667 0.83333333333333 0.43836073203431 32.065 > 1 16 S 0.91666666666667 0.83333333333333 0.43836073203431 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.93836073203431 32.065 > 1 18 S 0.41666666666667 0.08333333333333 0.93836073203431 32.065 > 1 19 S 0.66666666666667 0.08333333333333 0.93836073203431 32.065 > 1 20 S 0.91666666666667 0.08333333333333 0.93836073203431 32.065 > 1 21 S 0.16666666666667 0.33333333333333 0.93836073203431 32.065 > 1 22 S 0.41666666666667 0.33333333333333 0.93836073203431 32.065 > 1 23 S 0.66666666666667 0.33333333333333 0.93836073203431 32.065 > 1 24 S 0.91666666666667 0.33333333333333 0.93836073203431 32.065 > 1 25 S 0.16666666666667 0.58333333333333 0.93836073203431 32.065 > 1 26 S 0.41666666666667 0.58333333333333 0.93836073203431 32.065 > 1 27 S 0.66666666666667 0.58333333333333 0.93836073203431 32.065 > 1 28 S 0.91666666666667 0.58333333333333 0.93836073203431 32.065 > 1 29 S 0.16666666666667 0.83333333333333 0.93836073203431 32.065 > 1 30 S 0.41666666666667 0.83333333333333 0.93836073203431 32.065 > 1 31 S 0.66666666666667 0.83333333333333 0.93836073203431 32.065 > 1 32 S 0.91666666666667 0.83333333333333 0.93836073203431 32.065 > 1 33 S 0.08333333333333 0.16666666666667 0.18836073203431 32.065 > 33 34 S 0.33333333333333 0.16666666666667 0.18836073203431 32.065 > 33 35 S 0.58333333333333 0.16666666666667 0.18836073203431 32.065 > 33 36 S 0.83333333333333 0.16666666666667 0.18836073203431 32.065 > 33 37 S 0.08333333333333 0.41666666666667 0.18836073203431 32.065 > 33 38 S 0.33333333333333 0.41666666666667 0.18836073203431 32.065 > 33 39 S 0.58333333333333 0.41666666666667 0.18836073203431 32.065 > 33 40 S 0.83333333333333 0.41666666666667 0.18836073203431 32.065 > 33 41 S 0.08333333333333 0.66666666666667 0.18836073203431 32.065 > 33 42 S 0.33333333333333 0.66666666666667 0.18836073203431 32.065 > 33 43 S 0.58333333333333 0.66666666666667 0.18836073203431 32.065 > 33 44 S 0.83333333333333 0.66666666666667 0.18836073203431 32.065 > 33 45 S 0.08333333333333 0.91666666666667 0.18836073203431 32.065 > 33 46 S 0.33333333333333 0.91666666666667 0.18836073203431 32.065 > 33 47 S 0.58333333333333 0.91666666666667 0.18836073203431 32.065 > 33 48 S 0.83333333333333 0.91666666666667 0.18836073203431 32.065 > 33 49 S 0.08333333333333 0.16666666666667 0.68836073203431 32.065 > 33 50 S 0.33333333333333 0.16666666666667 0.68836073203431 32.065 > 33 51 S 0.58333333333333 0.16666666666667 0.68836073203431 32.065 > 33 52 S 0.83333333333333 0.16666666666667 0.68836073203431 32.065 > 33 53 S 0.08333333333333 0.41666666666667 0.68836073203431 32.065 > 33 54 S 0.33333333333333 0.41666666666667 0.68836073203431 32.065 > 33 55 S 0.58333333333333 0.41666666666667 0.68836073203431 32.065 > 33 56 S 0.83333333333333 0.41666666666667 0.68836073203431 32.065 > 33 57 S 0.08333333333333 0.66666666666667 0.68836073203431 32.065 > 33 58 S 0.33333333333333 0.66666666666667 0.68836073203431 32.065 > 33 59 S 0.58333333333333 0.66666666666667 0.68836073203431 32.065 > 33 60 S 0.83333333333333 0.66666666666667 0.68836073203431 32.065 > 33 61 S 0.08333333333333 0.91666666666667 0.68836073203431 32.065 > 33 62 S 0.33333333333333 0.91666666666667 0.68836073203431 32.065 > 33 63 S 0.58333333333333 0.91666666666667 0.68836073203431 32.065 > 33 64 S 0.83333333333333 0.91666666666667 0.68836073203431 32.065 > 33 65 Cd 0.16666666666667 0.08333333333333 0.25013926796569 112.411 > 65 66 Cd 0.41666666666667 0.08333333333333 0.25013926796569 112.411 > 65 67 Cd 0.66666666666667 0.08333333333333 0.25013926796569 112.411 > 65 68 Cd 0.91666666666667 0.08333333333333 0.25013926796569 112.411 > 65 69 Cd 0.16666666666667 0.33333333333333 0.25013926796569 112.411 > 65 70 Cd 0.41666666666667 0.33333333333333 0.25013926796569 112.411 > 65 71 Cd 0.66666666666667 0.33333333333333 0.25013926796569 112.411 > 65 72 Cd 0.91666666666667 0.33333333333333 0.25013926796569 112.411 > 65 73 Cd 0.16666666666667 0.58333333333333 0.25013926796569 112.411 > 65 74 Cd 0.41666666666667 0.58333333333333 0.25013926796569 112.411 > 65 75 Cd 0.66666666666667 0.58333333333333 0.25013926796569 112.411 > 65 76 Cd 0.91666666666667 0.58333333333333 0.25013926796569 112.411 > 65 77 Cd 0.16666666666667 0.83333333333333 0.25013926796569 112.411 > 65 78 Cd 0.41666666666667 0.83333333333333 0.25013926796569 112.411 > 65 79 Cd 0.66666666666667 0.83333333333333 0.25013926796569 112.411 > 65 80 Cd 0.91666666666667 0.83333333333333 0.25013926796569 112.411 > 65 81 Cd 0.16666666666667 0.08333333333333 0.75013926796569 112.411 > 65 82 Cd 0.41666666666667 0.08333333333333 0.75013926796569 112.411 > 65 83 Cd 0.66666666666667 0.08333333333333 0.75013926796569 112.411 > 65 84 Cd 0.91666666666667 0.08333333333333 0.75013926796569 112.411 > 65 85 Cd 0.16666666666667 0.33333333333333 0.75013926796569 112.411 > 65 86 Cd 0.41666666666667 0.33333333333333 0.75013926796569 112.411 > 65 87 Cd 0.66666666666667 0.33333333333333 0.75013926796569 112.411 > 65 88 Cd 0.91666666666667 0.33333333333333 0.75013926796569 112.411 > 65 89 Cd 0.16666666666667 0.58333333333333 0.75013926796569 112.411 > 65 90 Cd 0.41666666666667 0.58333333333333 0.75013926796569 112.411 > 65 91 Cd 0.66666666666667 0.58333333333333 0.75013926796569 112.411 > 65 92 Cd 0.91666666666667 0.58333333333333 0.75013926796569 112.411 > 65 93 Cd 0.16666666666667 0.83333333333333 0.75013926796569 112.411 > 65 94 Cd 0.41666666666667 0.83333333333333 0.75013926796569 112.411 > 65 95 Cd 0.66666666666667 0.83333333333333 0.75013926796569 112.411 > 65 96 Cd 0.91666666666667 0.83333333333333 0.75013926796569 112.411 > 65 97 Cd 0.08333333333333 0.16666666666667 0.00013926796569 112.411 > 97 98 Cd 0.33333333333333 0.16666666666667 0.00013926796569 112.411 > 97 99 Cd 0.58333333333333 0.16666666666667 0.00013926796569 112.411 > 97 100 Cd 0.83333333333333 0.16666666666667 0.00013926796569 112.411 > 97 101 Cd 0.08333333333333 0.41666666666667 0.00013926796569 112.411 > 97 102 Cd 0.33333333333333 0.41666666666667 0.00013926796569 112.411 > 97 103 Cd 0.58333333333333 0.41666666666667 0.00013926796569 112.411 > 97 104 Cd 0.83333333333333 0.41666666666667 0.00013926796569 112.411 > 97 105 Cd 0.08333333333333 0.66666666666667 0.00013926796569 112.411 > 97 106 Cd 0.33333333333333 0.66666666666667 0.00013926796569 112.411 > 97 107 Cd 0.58333333333333 0.66666666666667 0.00013926796569 112.411 > 97 108 Cd 0.83333333333333 0.66666666666667 0.00013926796569 112.411 > 97 109 Cd 0.08333333333333 0.91666666666667 0.00013926796569 112.411 > 97 110 Cd 0.33333333333333 0.91666666666667 0.00013926796569 112.411 > 97 111 Cd 0.58333333333333 0.91666666666667 0.00013926796569 112.411 > 97 112 Cd 0.83333333333333 0.91666666666667 0.00013926796569 112.411 > 97 113 Cd 0.08333333333333 0.16666666666667 0.50013926796569 112.411 > 97 114 Cd 0.33333333333333 0.16666666666667 0.50013926796569 112.411 > 97 115 Cd 0.58333333333333 0.16666666666667 0.50013926796569 112.411 > 97 116 Cd 0.83333333333333 0.16666666666667 0.50013926796569 112.411 > 97 117 Cd 0.08333333333333 0.41666666666667 0.50013926796569 112.411 > 97 118 Cd 0.33333333333333 0.41666666666667 0.50013926796569 112.411 > 97 119 Cd 0.58333333333333 0.41666666666667 0.50013926796569 112.411 > 97 120 Cd 0.83333333333333 0.41666666666667 0.50013926796569 112.411 > 97 121 Cd 0.08333333333333 0.66666666666667 0.50013926796569 112.411 > 97 122 Cd 0.33333333333333 0.66666666666667 0.50013926796569 112.411 > 97 123 Cd 0.58333333333333 0.66666666666667 0.50013926796569 112.411 > 97 124 Cd 0.83333333333333 0.66666666666667 0.50013926796569 112.411 > 97 125 Cd 0.08333333333333 0.91666666666667 0.50013926796569 112.411 > 97 126 Cd 0.33333333333333 0.91666666666667 0.50013926796569 112.411 > 97 127 Cd 0.58333333333333 0.91666666666667 0.50013926796569 112.411 > 97 128 Cd 0.83333333333333 0.91666666666667 0.50013926796569 112.411 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.4947083 -0.0000000 0.0000000 0.0000000 7.4947083 0.0000000 0.0000000 0.0000000 7.5360942 -------------------------- Born effective charges -------------------------- 1 S -2.1930538 -0.0000000 0.0000000 0.0000000 -2.1930538 0.0000000 0.0000000 0.0000000 -2.2856256 2 S -2.1930538 -0.0000000 0.0000000 0.0000000 -2.1930538 0.0000000 0.0000000 0.0000000 -2.2856256 3 Cd 2.1930538 0.0000000 0.0000000 0.0000000 2.1930538 0.0000000 0.0000000 0.0000000 2.2856256 4 Cd 2.1930538 0.0000000 0.0000000 0.0000000 2.1930538 0.0000000 0.0000000 0.0000000 2.2856256 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000031 (zzz) 0.00000031 (zzz) 0.00000031 (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-07 18:35:44]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:35:45]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.123353369962359 0.000000000000000 0.000000000000000 b -2.061676684981180 3.570928767167578 0.000000000000000 c 0.000000000000000 0.000000000000000 6.710167980000000 Atomic positions (fractional): 1 S 0.66666666666667 0.33333333333333 0.87672146406862 32.065 2 S 0.33333333333333 0.66666666666667 0.37672146406862 32.065 3 Cd 0.66666666666667 0.33333333333333 0.50027853593138 112.411 4 Cd 0.33333333333333 0.66666666666667 0.00027853593138 112.411 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.493413479849437 0.000000000000000 0.000000000000000 b -8.246706739924720 14.283715068670311 0.000000000000000 c 0.000000000000000 0.000000000000000 13.420335959999999 Atomic positions (fractional): 1 S 0.16666666666667 0.08333333333333 0.43836073203431 32.065 > 1 2 S 0.41666666666667 0.08333333333333 0.43836073203431 32.065 > 1 3 S 0.66666666666667 0.08333333333333 0.43836073203431 32.065 > 1 4 S 0.91666666666667 0.08333333333333 0.43836073203431 32.065 > 1 5 S 0.16666666666667 0.33333333333333 0.43836073203431 32.065 > 1 6 S 0.41666666666667 0.33333333333333 0.43836073203431 32.065 > 1 7 S 0.66666666666667 0.33333333333333 0.43836073203431 32.065 > 1 8 S 0.91666666666667 0.33333333333333 0.43836073203431 32.065 > 1 9 S 0.16666666666667 0.58333333333333 0.43836073203431 32.065 > 1 10 S 0.41666666666667 0.58333333333333 0.43836073203431 32.065 > 1 11 S 0.66666666666667 0.58333333333333 0.43836073203431 32.065 > 1 12 S 0.91666666666667 0.58333333333333 0.43836073203431 32.065 > 1 13 S 0.16666666666667 0.83333333333333 0.43836073203431 32.065 > 1 14 S 0.41666666666667 0.83333333333333 0.43836073203431 32.065 > 1 15 S 0.66666666666667 0.83333333333333 0.43836073203431 32.065 > 1 16 S 0.91666666666667 0.83333333333333 0.43836073203431 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.93836073203431 32.065 > 1 18 S 0.41666666666667 0.08333333333333 0.93836073203431 32.065 > 1 19 S 0.66666666666667 0.08333333333333 0.93836073203431 32.065 > 1 20 S 0.91666666666667 0.08333333333333 0.93836073203431 32.065 > 1 21 S 0.16666666666667 0.33333333333333 0.93836073203431 32.065 > 1 22 S 0.41666666666667 0.33333333333333 0.93836073203431 32.065 > 1 23 S 0.66666666666667 0.33333333333333 0.93836073203431 32.065 > 1 24 S 0.91666666666667 0.33333333333333 0.93836073203431 32.065 > 1 25 S 0.16666666666667 0.58333333333333 0.93836073203431 32.065 > 1 26 S 0.41666666666667 0.58333333333333 0.93836073203431 32.065 > 1 27 S 0.66666666666667 0.58333333333333 0.93836073203431 32.065 > 1 28 S 0.91666666666667 0.58333333333333 0.93836073203431 32.065 > 1 29 S 0.16666666666667 0.83333333333333 0.93836073203431 32.065 > 1 30 S 0.41666666666667 0.83333333333333 0.93836073203431 32.065 > 1 31 S 0.66666666666667 0.83333333333333 0.93836073203431 32.065 > 1 32 S 0.91666666666667 0.83333333333333 0.93836073203431 32.065 > 1 33 S 0.08333333333333 0.16666666666667 0.18836073203431 32.065 > 33 34 S 0.33333333333333 0.16666666666667 0.18836073203431 32.065 > 33 35 S 0.58333333333333 0.16666666666667 0.18836073203431 32.065 > 33 36 S 0.83333333333333 0.16666666666667 0.18836073203431 32.065 > 33 37 S 0.08333333333333 0.41666666666667 0.18836073203431 32.065 > 33 38 S 0.33333333333333 0.41666666666667 0.18836073203431 32.065 > 33 39 S 0.58333333333333 0.41666666666667 0.18836073203431 32.065 > 33 40 S 0.83333333333333 0.41666666666667 0.18836073203431 32.065 > 33 41 S 0.08333333333333 0.66666666666667 0.18836073203431 32.065 > 33 42 S 0.33333333333333 0.66666666666667 0.18836073203431 32.065 > 33 43 S 0.58333333333333 0.66666666666667 0.18836073203431 32.065 > 33 44 S 0.83333333333333 0.66666666666667 0.18836073203431 32.065 > 33 45 S 0.08333333333333 0.91666666666667 0.18836073203431 32.065 > 33 46 S 0.33333333333333 0.91666666666667 0.18836073203431 32.065 > 33 47 S 0.58333333333333 0.91666666666667 0.18836073203431 32.065 > 33 48 S 0.83333333333333 0.91666666666667 0.18836073203431 32.065 > 33 49 S 0.08333333333333 0.16666666666667 0.68836073203431 32.065 > 33 50 S 0.33333333333333 0.16666666666667 0.68836073203431 32.065 > 33 51 S 0.58333333333333 0.16666666666667 0.68836073203431 32.065 > 33 52 S 0.83333333333333 0.16666666666667 0.68836073203431 32.065 > 33 53 S 0.08333333333333 0.41666666666667 0.68836073203431 32.065 > 33 54 S 0.33333333333333 0.41666666666667 0.68836073203431 32.065 > 33 55 S 0.58333333333333 0.41666666666667 0.68836073203431 32.065 > 33 56 S 0.83333333333333 0.41666666666667 0.68836073203431 32.065 > 33 57 S 0.08333333333333 0.66666666666667 0.68836073203431 32.065 > 33 58 S 0.33333333333333 0.66666666666667 0.68836073203431 32.065 > 33 59 S 0.58333333333333 0.66666666666667 0.68836073203431 32.065 > 33 60 S 0.83333333333333 0.66666666666667 0.68836073203431 32.065 > 33 61 S 0.08333333333333 0.91666666666667 0.68836073203431 32.065 > 33 62 S 0.33333333333333 0.91666666666667 0.68836073203431 32.065 > 33 63 S 0.58333333333333 0.91666666666667 0.68836073203431 32.065 > 33 64 S 0.83333333333333 0.91666666666667 0.68836073203431 32.065 > 33 65 Cd 0.16666666666667 0.08333333333333 0.25013926796569 112.411 > 65 66 Cd 0.41666666666667 0.08333333333333 0.25013926796569 112.411 > 65 67 Cd 0.66666666666667 0.08333333333333 0.25013926796569 112.411 > 65 68 Cd 0.91666666666667 0.08333333333333 0.25013926796569 112.411 > 65 69 Cd 0.16666666666667 0.33333333333333 0.25013926796569 112.411 > 65 70 Cd 0.41666666666667 0.33333333333333 0.25013926796569 112.411 > 65 71 Cd 0.66666666666667 0.33333333333333 0.25013926796569 112.411 > 65 72 Cd 0.91666666666667 0.33333333333333 0.25013926796569 112.411 > 65 73 Cd 0.16666666666667 0.58333333333333 0.25013926796569 112.411 > 65 74 Cd 0.41666666666667 0.58333333333333 0.25013926796569 112.411 > 65 75 Cd 0.66666666666667 0.58333333333333 0.25013926796569 112.411 > 65 76 Cd 0.91666666666667 0.58333333333333 0.25013926796569 112.411 > 65 77 Cd 0.16666666666667 0.83333333333333 0.25013926796569 112.411 > 65 78 Cd 0.41666666666667 0.83333333333333 0.25013926796569 112.411 > 65 79 Cd 0.66666666666667 0.83333333333333 0.25013926796569 112.411 > 65 80 Cd 0.91666666666667 0.83333333333333 0.25013926796569 112.411 > 65 81 Cd 0.16666666666667 0.08333333333333 0.75013926796569 112.411 > 65 82 Cd 0.41666666666667 0.08333333333333 0.75013926796569 112.411 > 65 83 Cd 0.66666666666667 0.08333333333333 0.75013926796569 112.411 > 65 84 Cd 0.91666666666667 0.08333333333333 0.75013926796569 112.411 > 65 85 Cd 0.16666666666667 0.33333333333333 0.75013926796569 112.411 > 65 86 Cd 0.41666666666667 0.33333333333333 0.75013926796569 112.411 > 65 87 Cd 0.66666666666667 0.33333333333333 0.75013926796569 112.411 > 65 88 Cd 0.91666666666667 0.33333333333333 0.75013926796569 112.411 > 65 89 Cd 0.16666666666667 0.58333333333333 0.75013926796569 112.411 > 65 90 Cd 0.41666666666667 0.58333333333333 0.75013926796569 112.411 > 65 91 Cd 0.66666666666667 0.58333333333333 0.75013926796569 112.411 > 65 92 Cd 0.91666666666667 0.58333333333333 0.75013926796569 112.411 > 65 93 Cd 0.16666666666667 0.83333333333333 0.75013926796569 112.411 > 65 94 Cd 0.41666666666667 0.83333333333333 0.75013926796569 112.411 > 65 95 Cd 0.66666666666667 0.83333333333333 0.75013926796569 112.411 > 65 96 Cd 0.91666666666667 0.83333333333333 0.75013926796569 112.411 > 65 97 Cd 0.08333333333333 0.16666666666667 0.00013926796569 112.411 > 97 98 Cd 0.33333333333333 0.16666666666667 0.00013926796569 112.411 > 97 99 Cd 0.58333333333333 0.16666666666667 0.00013926796569 112.411 > 97 100 Cd 0.83333333333333 0.16666666666667 0.00013926796569 112.411 > 97 101 Cd 0.08333333333333 0.41666666666667 0.00013926796569 112.411 > 97 102 Cd 0.33333333333333 0.41666666666667 0.00013926796569 112.411 > 97 103 Cd 0.58333333333333 0.41666666666667 0.00013926796569 112.411 > 97 104 Cd 0.83333333333333 0.41666666666667 0.00013926796569 112.411 > 97 105 Cd 0.08333333333333 0.66666666666667 0.00013926796569 112.411 > 97 106 Cd 0.33333333333333 0.66666666666667 0.00013926796569 112.411 > 97 107 Cd 0.58333333333333 0.66666666666667 0.00013926796569 112.411 > 97 108 Cd 0.83333333333333 0.66666666666667 0.00013926796569 112.411 > 97 109 Cd 0.08333333333333 0.91666666666667 0.00013926796569 112.411 > 97 110 Cd 0.33333333333333 0.91666666666667 0.00013926796569 112.411 > 97 111 Cd 0.58333333333333 0.91666666666667 0.00013926796569 112.411 > 97 112 Cd 0.83333333333333 0.91666666666667 0.00013926796569 112.411 > 97 113 Cd 0.08333333333333 0.16666666666667 0.50013926796569 112.411 > 97 114 Cd 0.33333333333333 0.16666666666667 0.50013926796569 112.411 > 97 115 Cd 0.58333333333333 0.16666666666667 0.50013926796569 112.411 > 97 116 Cd 0.83333333333333 0.16666666666667 0.50013926796569 112.411 > 97 117 Cd 0.08333333333333 0.41666666666667 0.50013926796569 112.411 > 97 118 Cd 0.33333333333333 0.41666666666667 0.50013926796569 112.411 > 97 119 Cd 0.58333333333333 0.41666666666667 0.50013926796569 112.411 > 97 120 Cd 0.83333333333333 0.41666666666667 0.50013926796569 112.411 > 97 121 Cd 0.08333333333333 0.66666666666667 0.50013926796569 112.411 > 97 122 Cd 0.33333333333333 0.66666666666667 0.50013926796569 112.411 > 97 123 Cd 0.58333333333333 0.66666666666667 0.50013926796569 112.411 > 97 124 Cd 0.83333333333333 0.66666666666667 0.50013926796569 112.411 > 97 125 Cd 0.08333333333333 0.91666666666667 0.50013926796569 112.411 > 97 126 Cd 0.33333333333333 0.91666666666667 0.50013926796569 112.411 > 97 127 Cd 0.58333333333333 0.91666666666667 0.50013926796569 112.411 > 97 128 Cd 0.83333333333333 0.91666666666667 0.50013926796569 112.411 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 7.4947083 -0.0000000 0.0000000 0.0000000 7.4947083 0.0000000 0.0000000 0.0000000 7.5360942 -------------------------- Born effective charges -------------------------- 1 S -2.1930538 -0.0000000 0.0000000 0.0000000 -2.1930538 0.0000000 0.0000000 0.0000000 -2.2856256 2 S -2.1930538 -0.0000000 0.0000000 0.0000000 -2.1930538 0.0000000 0.0000000 0.0000000 -2.2856256 3 Cd 2.1930538 0.0000000 0.0000000 0.0000000 2.1930538 0.0000000 0.0000000 0.0000000 2.2856256 4 Cd 2.1930538 0.0000000 0.0000000 0.0000000 2.1930538 0.0000000 0.0000000 0.0000000 2.2856256 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000031 (zzz) 0.00000031 (zzz) 0.00000031 (zzz) Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 14 14 7 ] 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.90, Number of G-points: 313, Lambda: 0.26 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/96) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 1.229 ( 0.000 -0.000 0.000) 0.000 1.229 ( -0.000 0.000 0.000) 0.000 4.025 ( 0.000 0.000 0.000) 0.000 6.975 ( 0.000 0.000 0.000) 0.000 7.200 ( 0.000 0.000 -0.000) 0.000 7.200 ( 0.000 -0.000 -0.000) 0.000 7.574 ( 0.000 -0.000 0.000) 0.000 7.574 ( 0.000 -0.000 -0.000) 0.000 8.939 ( 0.000 0.000 -0.000) 0.000 ======================= Grid point 1 (2/96) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.328 ( 13.839 7.990 0.000) 15.980 0.341 ( 14.791 8.540 0.000) 17.080 0.852 ( 36.043 20.810 0.000) 41.619 1.252 ( 2.147 1.240 0.000) 2.480 1.443 ( 15.956 9.212 0.000) 18.424 3.996 ( -2.243 -1.295 0.000) 2.590 7.004 ( 2.480 1.432 0.000) 2.864 7.229 ( 2.471 1.427 0.000) 2.853 7.616 ( 3.480 2.009 0.000) 4.018 7.673 ( 8.093 4.673 0.000) 9.345 8.687 ( 2.849 1.645 0.000) 3.290 8.924 ( -1.178 -0.680 0.000) 1.360 ======================= Grid point 2 (3/96) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.630 ( 12.152 7.016 0.000) 14.032 0.676 ( 14.161 8.176 0.000) 16.352 1.331 ( 4.582 2.645 0.000) 5.291 1.646 ( 31.788 18.353 0.000) 36.706 1.834 ( 16.156 9.328 0.000) 18.655 3.943 ( -1.802 -1.041 0.000) 2.081 7.089 ( 4.801 2.772 0.000) 5.543 7.307 ( 4.012 2.316 0.000) 4.632 7.726 ( 5.836 3.369 0.000) 6.739 7.916 ( 12.037 6.949 0.000) 13.899 8.770 ( 3.736 2.157 0.000) 4.313 8.890 ( -1.388 -0.801 0.000) 1.602 ======================= Grid point 3 (4/96) ======================= q-point: ( 0.21 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.878 ( 9.134 5.274 0.000) 10.547 0.992 ( 13.110 7.569 0.000) 15.139 1.450 ( 5.336 3.081 0.000) 6.161 2.071 ( 7.277 4.201 0.000) 8.402 2.371 ( 27.619 15.946 0.000) 31.892 3.936 ( 1.524 0.880 0.000) 1.760 7.223 ( 6.691 3.863 0.000) 7.727 7.407 ( 4.513 2.606 0.000) 5.211 7.873 ( 6.618 3.821 0.000) 7.642 8.183 ( 10.176 5.875 0.000) 11.750 8.801 ( -1.301 -0.751 0.000) 1.502 8.909 ( 3.151 1.819 0.000) 3.639 ======================= Grid point 4 (5/96) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.047 ( 5.451 3.147 0.000) 6.294 1.278 ( 11.600 6.698 0.000) 13.395 1.556 ( 3.462 1.999 0.000) 3.997 2.149 ( 0.396 0.229 0.000) 0.457 2.953 ( 22.544 13.016 0.000) 26.031 4.007 ( 4.210 2.430 0.000) 4.861 7.391 ( 7.755 4.477 0.000) 8.954 7.510 ( 4.307 2.487 0.000) 4.974 8.021 ( 6.017 3.474 0.000) 6.948 8.351 ( 3.906 2.255 0.000) 4.510 8.766 ( -1.227 -0.708 0.000) 1.416 8.978 ( 2.273 1.313 0.000) 2.625 ======================= Grid point 5 (6/96) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.136 ( 2.464 1.423 0.000) 2.845 1.524 ( 9.557 5.518 0.000) 11.036 1.601 ( 0.411 0.237 0.000) 0.474 2.126 ( -1.848 -1.067 0.000) 2.133 3.406 ( 16.493 9.522 0.000) 19.045 4.109 ( 4.143 2.392 0.000) 4.784 7.571 ( 7.511 4.336 0.000) 8.673 7.601 ( 3.523 2.034 0.000) 4.067 8.144 ( 4.468 2.579 0.000) 5.159 8.373 ( -1.275 -0.736 0.000) 1.473 8.742 ( -1.214 -0.701 0.000) 1.402 9.009 ( 0.499 0.288 0.000) 0.576 ======================= Grid point 6 (7/96) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.172 ( 0.896 0.517 0.000) 1.034 1.587 ( -1.180 -0.681 0.000) 1.363 1.711 ( 6.334 3.657 0.000) 7.314 2.078 ( -2.049 -1.183 0.000) 2.366 3.706 ( 9.187 5.304 0.000) 10.609 4.184 ( 2.202 1.271 0.000) 2.543 7.667 ( 2.047 1.182 0.000) 2.364 7.723 ( 5.204 3.005 0.000) 6.009 8.224 ( 2.363 1.364 0.000) 2.729 8.326 ( -2.156 -1.245 0.000) 2.489 8.706 ( -1.640 -0.947 0.000) 1.894 9.009 ( -0.296 -0.171 0.000) 0.342 ======================= Grid point 7 (8/96) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.182 ( 0.000 0.000 0.000) 0.000 1.569 ( 0.000 0.000 0.000) 0.000 1.792 ( 0.000 0.000 0.000) 0.000 2.048 ( 0.000 0.000 0.000) 0.000 3.815 ( 0.000 0.000 0.000) 0.000 4.209 ( 0.000 0.000 0.000) 0.000 7.691 ( 0.000 0.000 0.000) 0.000 7.788 ( 0.000 0.000 0.000) 0.000 8.251 ( 0.000 0.000 0.000) 0.000 8.297 ( 0.000 0.000 0.000) 0.000 8.683 ( 0.000 0.000 0.000) 0.000 9.004 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 16 (9/96) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.586 ( 8.182 14.171 0.000) 16.364 0.600 ( 9.055 15.684 0.000) 18.111 1.344 ( 6.256 10.836 0.000) 12.512 1.448 ( 16.334 28.290 0.000) 32.667 1.680 ( 8.705 15.078 0.000) 17.411 3.955 ( -1.358 -2.352 0.000) 2.716 7.061 ( 2.425 4.200 0.000) 4.850 7.299 ( 2.906 5.033 0.000) 5.811 7.684 ( 3.059 5.298 0.000) 6.118 7.842 ( 6.361 11.017 0.000) 12.721 8.742 ( 2.118 3.669 0.000) 4.236 8.898 ( -1.053 -1.824 0.000) 2.106 ======================= Grid point 17 (10/96) ======================= q-point: ( 0.14 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.876 ( 10.309 11.277 0.000) 15.279 0.891 ( 5.462 16.445 0.000) 17.328 1.482 ( 1.208 10.285 0.000) 10.356 1.962 ( 9.904 7.571 0.000) 12.466 2.097 ( 23.694 21.843 0.000) 32.227 3.925 ( 0.125 -0.289 0.000) 0.314 7.170 ( 4.609 5.296 0.000) 7.021 7.411 ( 2.410 7.556 0.000) 7.931 7.819 ( 5.411 6.075 0.000) 8.135 8.079 ( 8.796 8.414 0.000) 12.172 8.812 ( 0.778 1.942 0.000) 2.092 8.867 ( 1.090 -1.411 0.000) 1.783 ======================= Grid point 18 (11/96) ======================= q-point: ( 0.21 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.107 ( 0.909 15.464 0.000) 15.491 1.160 ( 10.490 9.105 0.000) 13.890 1.599 ( -0.008 10.323 0.000) 10.323 2.118 ( 3.730 1.049 0.000) 3.875 2.696 ( 22.758 16.145 0.000) 27.903 3.962 ( 3.123 1.738 0.000) 3.574 7.322 ( 6.234 5.985 0.000) 8.642 7.525 ( 1.252 8.701 0.000) 8.791 7.976 ( 5.083 6.322 0.000) 8.112 8.281 ( 6.843 3.482 0.000) 7.678 8.787 ( -2.144 -0.620 0.000) 2.231 8.925 ( 4.095 -0.231 0.000) 4.101 ======================= Grid point 19 (12/96) ======================= q-point: ( 0.29 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.233 ( -2.980 13.685 0.000) 14.006 1.417 ( 9.577 7.287 0.000) 12.034 1.675 ( -2.530 8.794 0.000) 9.150 2.139 ( -0.608 -0.789 0.000) 0.996 3.196 ( 18.985 11.098 0.000) 21.991 4.052 ( 4.512 1.817 0.000) 4.864 7.498 ( 6.842 6.099 0.000) 9.165 7.626 ( 0.277 8.638 0.000) 8.642 8.109 ( 3.580 4.895 0.000) 6.064 8.372 ( 1.772 0.306 0.000) 1.798 8.749 ( -1.254 -1.157 0.000) 1.706 8.978 ( 2.648 -1.004 0.000) 2.832 ======================= Grid point 20 (13/96) ======================= q-point: ( 0.36 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.288 ( -5.258 12.168 0.000) 13.256 1.628 ( 6.502 5.639 0.000) 8.607 1.691 ( -4.051 7.465 0.000) 8.493 2.103 ( -2.091 -1.060 0.000) 2.344 3.567 ( 13.325 6.574 0.000) 14.859 4.137 ( 3.763 0.250 0.000) 3.772 7.669 ( 5.788 5.373 0.000) 7.898 7.705 ( -0.833 7.926 0.000) 7.970 8.189 ( 1.924 1.086 0.000) 2.209 8.375 ( -1.673 1.734 0.000) 2.409 8.712 ( -1.120 -2.197 0.000) 2.466 8.998 ( 1.001 -1.016 0.000) 1.426 ======================= Grid point 21 (14/96) ======================= q-point: ( 0.43 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.305 ( -6.097 11.302 0.000) 12.842 1.661 ( -5.220 7.510 0.000) 9.146 1.781 ( 2.960 3.317 0.000) 4.446 2.057 ( -1.548 -0.878 0.000) 1.780 3.775 ( 5.409 1.678 0.000) 5.663 4.180 ( 2.103 -1.674 0.000) 2.688 7.749 ( -2.388 6.535 0.000) 6.958 7.788 ( 1.878 3.613 0.000) 4.072 8.203 ( 1.192 -2.487 0.000) 2.758 8.369 ( -2.602 4.479 0.000) 5.180 8.675 ( -0.237 -1.990 0.000) 2.004 9.000 ( 0.244 -0.557 0.000) 0.608 ======================= Grid point 31 (15/96) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.111 ( 6.888 11.931 0.000) 13.777 1.203 ( 7.852 13.600 0.000) 15.704 1.676 ( 4.829 8.365 0.000) 9.659 2.077 ( 2.500 4.330 0.000) 4.999 2.560 ( 13.523 23.422 0.000) 27.045 3.943 ( 1.214 2.103 0.000) 2.429 7.297 ( 4.241 7.346 0.000) 8.483 7.576 ( 4.735 8.200 0.000) 9.469 7.964 ( 4.598 7.964 0.000) 9.195 8.227 ( 3.498 6.059 0.000) 6.996 8.820 ( -1.010 -1.750 0.000) 2.020 8.865 ( 1.092 1.892 0.000) 2.185 ======================= Grid point 32 (16/96) ======================= q-point: ( 0.21 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.344 ( 7.304 9.530 0.000) 12.007 1.419 ( 1.321 13.854 0.000) 13.917 1.807 ( -0.083 9.400 0.000) 9.401 2.131 ( 0.999 0.387 0.000) 1.071 3.013 ( 15.733 15.608 0.000) 22.162 4.004 ( 2.941 2.275 0.000) 3.718 7.458 ( 5.540 7.565 0.000) 9.377 7.722 ( 1.400 10.123 0.000) 10.220 8.115 ( 3.252 7.236 0.000) 7.933 8.318 ( 3.607 0.365 0.000) 3.625 8.775 ( -2.512 -0.691 0.000) 2.605 8.911 ( 3.161 -0.765 0.000) 3.252 ======================= Grid point 33 (17/96) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.506 ( -2.279 11.531 0.000) 11.754 1.579 ( 5.495 9.272 0.000) 10.778 1.862 ( -3.979 8.905 0.000) 9.754 2.125 ( -1.267 -0.607 0.000) 1.405 3.401 ( 14.044 9.640 0.000) 17.035 4.075 ( 3.291 0.408 0.000) 3.316 7.632 ( 5.641 7.181 0.000) 9.132 7.827 ( -1.202 10.727 0.000) 10.794 8.194 ( -0.376 3.129 0.000) 3.151 8.379 ( 3.267 0.577 0.000) 3.317 8.720 ( -2.508 -1.708 0.000) 3.035 8.954 ( 2.708 -0.965 0.000) 2.875 ======================= Grid point 34 (18/96) ======================= q-point: ( 0.36 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.544 ( -5.971 11.898 0.000) 13.312 1.753 ( 4.661 6.623 0.000) 8.099 1.859 ( -5.857 8.536 0.000) 10.352 2.082 ( -2.394 -1.149 0.000) 2.656 3.677 ( 8.938 4.711 0.000) 10.103 4.117 ( 2.506 -2.228 0.000) 3.353 7.782 ( 3.503 5.818 0.000) 6.791 7.898 ( -3.042 10.758 0.000) 11.180 8.182 ( -1.276 -1.248 0.000) 1.784 8.436 ( 0.631 3.469 0.000) 3.526 8.660 ( -1.410 -2.940 0.000) 3.261 8.982 ( 1.269 -0.257 0.000) 1.295 ======================= Grid point 35 (19/96) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.551 ( -6.744 11.681 0.000) 13.488 1.843 ( -5.440 9.422 0.000) 10.879 1.843 ( -1.677 2.905 0.000) 3.354 2.045 ( 0.155 -0.269 0.000) 0.311 3.783 ( 0.305 -0.529 0.000) 0.611 4.124 ( 2.087 -3.615 0.000) 4.175 7.845 ( -1.785 3.092 0.000) 3.571 7.930 ( -5.734 9.932 0.000) 11.469 8.157 ( 1.071 -1.856 0.000) 2.143 8.461 ( -2.318 4.014 0.000) 4.635 8.629 ( 1.506 -2.608 0.000) 3.012 8.992 ( 0.071 -0.122 0.000) 0.141 ======================= Grid point 47 (20/96) ======================= q-point: ( 0.21 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.532 ( 5.230 9.059 0.000) 10.461 1.636 ( 4.528 7.843 0.000) 9.056 1.951 ( 2.795 4.840 0.000) 5.589 2.133 ( -0.209 -0.362 0.000) 0.418 3.318 ( 8.500 14.722 0.000) 17.000 4.046 ( 0.998 1.729 0.000) 1.996 7.619 ( 4.774 8.269 0.000) 9.548 7.902 ( 4.328 7.497 0.000) 8.657 8.252 ( 3.256 5.640 0.000) 6.513 8.306 ( -0.384 -0.665 0.000) 0.768 8.745 ( -1.478 -2.560 0.000) 2.957 8.914 ( 0.802 1.390 0.000) 1.605 ======================= Grid point 48 (21/96) ======================= q-point: ( 0.29 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.688 ( 1.726 7.881 0.000) 8.068 1.757 ( 2.443 7.063 0.000) 7.474 2.004 ( -1.957 5.079 0.000) 5.443 2.109 ( -1.714 -1.281 0.000) 2.140 3.588 ( 8.360 9.080 0.000) 12.343 4.066 ( 0.914 -1.327 0.000) 1.611 7.779 ( 4.348 7.378 0.000) 8.564 8.020 ( 0.814 7.441 0.000) 7.485 8.244 ( -5.118 2.672 0.000) 5.773 8.389 ( 6.295 0.234 0.000) 6.299 8.676 ( -3.049 -2.856 0.000) 4.178 8.957 ( 1.947 1.419 0.000) 2.410 ======================= Grid point 49 (22/96) ======================= q-point: ( 0.36 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.739 ( -3.629 7.841 0.000) 8.640 1.879 ( 2.883 5.329 0.000) 6.058 2.002 ( -3.774 5.083 0.000) 6.331 2.053 ( -2.219 -1.638 0.000) 2.758 3.767 ( 3.512 4.503 0.000) 5.711 4.054 ( 1.162 -3.860 0.000) 4.032 7.896 ( 0.891 5.505 0.000) 5.576 8.094 ( 0.785 5.472 0.000) 5.528 8.172 ( -6.344 3.235 0.000) 7.121 8.488 ( 3.906 1.511 0.000) 4.188 8.593 ( -2.260 -3.868 0.000) 4.480 8.994 ( 0.477 1.310 0.000) 1.394 ======================= Grid point 62 (23/96) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.834 ( 3.462 5.996 0.000) 6.924 1.853 ( 2.038 3.529 0.000) 4.075 2.065 ( -1.839 -3.185 0.000) 3.678 2.066 ( 0.815 1.411 0.000) 1.629 3.758 ( 4.376 7.580 0.000) 8.752 4.022 ( -1.715 -2.970 0.000) 3.429 7.918 ( 3.530 6.115 0.000) 7.060 8.096 ( 0.305 0.529 0.000) 0.611 8.309 ( 1.451 2.514 0.000) 2.903 8.408 ( 1.638 2.836 0.000) 3.275 8.602 ( -2.589 -4.483 0.000) 5.177 8.998 ( 1.186 2.055 0.000) 2.373 ======================= Grid point 63 (24/96) ======================= q-point: ( 0.36 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.871 ( -3.095 5.361 0.000) 6.190 1.944 ( -0.290 0.502 0.000) 0.579 1.998 ( 1.218 -2.110 0.000) 2.436 2.078 ( -0.894 1.548 0.000) 1.788 3.856 ( -2.376 4.116 0.000) 4.753 3.971 ( 2.275 -3.940 0.000) 4.550 7.995 ( -2.307 3.996 0.000) 4.615 8.084 ( 1.392 -2.410 0.000) 2.783 8.326 ( -4.395 7.613 0.000) 8.791 8.495 ( 1.903 -3.296 0.000) 3.806 8.508 ( 0.929 -1.610 0.000) 1.859 9.021 ( -0.518 0.898 0.000) 1.037 ======================= Grid point 211 (25/96) ======================= q-point: ( 0.00 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.366 ( -0.000 0.000 16.460) 16.460 0.366 ( -0.000 -0.000 16.460) 16.460 0.957 ( 0.000 0.000 43.842) 43.842 1.220 ( 0.000 0.000 -0.911) 0.911 1.220 ( 0.000 -0.000 -0.911) 0.911 3.934 ( 0.000 -0.000 -8.474) 8.474 7.219 ( -0.000 0.000 1.762) 1.762 7.219 ( 0.000 -0.000 1.762) 1.762 7.557 ( 0.000 0.000 -1.636) 1.636 7.557 ( 0.000 0.000 -1.636) 1.636 8.607 ( -0.000 0.000 2.899) 2.899 8.936 ( 0.000 -0.000 -0.330) 0.330 ======================= Grid point 212 (26/96) ======================= q-point: ( 0.07 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.490 ( 9.192 5.307 12.200) 16.171 0.623 ( 3.741 2.160 23.858) 24.246 1.154 ( 25.354 14.638 24.905) 38.437 1.247 ( 2.393 1.382 -0.597) 2.827 1.421 ( 15.385 8.882 -1.216) 17.806 3.910 ( -1.851 -1.069 -8.087) 8.365 7.208 ( 0.300 0.173 9.937) 9.943 7.249 ( 2.515 1.452 1.806) 3.420 7.598 ( 3.433 1.982 -1.689) 4.309 7.655 ( 8.018 4.629 -1.627) 9.400 8.632 ( 2.689 1.552 -1.165) 3.316 8.921 ( -1.233 -0.712 -0.394) 1.477 ======================= Grid point 213 (27/96) ======================= q-point: ( 0.14 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.728 ( 10.553 6.093 8.194) 14.684 0.756 ( 9.151 5.283 13.702) 17.303 1.330 ( 4.652 2.686 -0.217) 5.376 1.764 ( 17.571 10.145 0.540) 20.297 1.861 ( 26.886 15.523 10.817) 32.876 3.871 ( -0.902 -0.521 -6.753) 6.833 7.248 ( 2.986 1.724 11.966) 12.453 7.328 ( 4.127 2.383 1.975) 5.159 7.707 ( 5.768 3.330 -1.820) 6.905 7.894 ( 11.736 6.776 -1.973) 13.695 8.723 ( 4.513 2.605 -3.182) 6.106 8.885 ( -1.636 -0.945 -0.544) 1.966 ======================= Grid point 214 (28/96) ======================= q-point: ( 0.21 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.951 ( 8.401 4.850 6.260) 11.546 0.998 ( 10.942 6.317 3.731) 13.174 1.449 ( 5.245 3.028 -0.252) 6.061 2.022 ( 6.244 3.605 -4.450) 8.472 2.501 ( 26.399 15.241 10.968) 32.396 3.891 ( 2.786 1.608 -4.404) 5.454 7.342 ( 5.055 2.918 10.010) 11.587 7.431 ( 4.671 2.697 2.272) 5.853 7.852 ( 6.563 3.789 -1.959) 7.827 8.152 ( 9.689 5.594 -2.795) 11.531 8.789 ( 0.019 0.011 -1.185) 1.185 8.883 ( 2.548 1.471 -1.962) 3.536 ======================= Grid point 215 (29/96) ======================= q-point: ( 0.29 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.107 ( 5.063 2.923 5.230) 7.844 1.242 ( 9.950 5.744 -1.699) 11.614 1.553 ( 3.422 1.976 -0.398) 3.971 2.082 ( -0.278 -0.160 -5.981) 5.989 3.059 ( 21.677 12.515 8.984) 26.594 3.990 ( 5.368 3.099 -1.795) 6.453 7.475 ( 6.311 3.644 7.302) 10.316 7.538 ( 4.444 2.566 2.592) 5.749 7.999 ( 5.983 3.455 -2.059) 7.209 8.309 ( 3.596 2.076 -3.710) 5.568 8.769 ( -0.948 -0.547 0.200) 1.112 8.952 ( 2.463 1.422 -2.237) 3.618 ======================= Grid point 216 (30/96) ======================= q-point: ( 0.36 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.190 ( 2.295 1.325 4.708) 5.402 1.454 ( 8.318 4.802 -5.404) 11.021 1.598 ( 0.497 0.287 -0.348) 0.671 2.043 ( -2.506 -1.447 -7.331) 7.882 3.496 ( 15.976 9.223 7.699) 19.989 4.115 ( 5.016 2.896 0.335) 5.801 7.624 ( 6.320 3.649 4.717) 8.690 7.632 ( 3.597 2.077 2.816) 5.018 8.121 ( 4.448 2.568 -2.118) 5.556 8.330 ( -1.135 -0.655 -3.824) 4.042 8.746 ( -1.247 -0.720 0.424) 1.501 8.987 ( 0.580 0.335 -2.032) 2.139 ======================= Grid point 217 (31/96) ======================= q-point: ( 0.43 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.223 ( 0.826 0.477 4.477) 4.578 1.587 ( -1.040 -0.600 -0.081) 1.204 1.619 ( 5.751 3.320 -7.674) 10.148 1.981 ( -2.611 -1.508 -8.592) 9.106 3.786 ( 8.886 5.130 6.892) 12.360 4.207 ( 2.735 1.579 1.809) 3.640 7.699 ( 2.069 1.195 2.911) 3.766 7.753 ( 4.452 2.571 2.712) 5.813 8.200 ( 2.353 1.359 -2.151) 3.466 8.287 ( -1.985 -1.146 -3.450) 4.142 8.710 ( -1.614 -0.932 0.457) 1.918 8.987 ( -0.290 -0.167 -1.996) 2.024 ======================= Grid point 218 (32/96) ======================= q-point: (-0.50 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.233 ( 0.000 0.000 4.397) 4.397 1.571 ( 0.000 0.000 0.095) 0.095 1.694 ( 0.000 0.000 -8.346) 8.346 1.944 ( -0.000 -0.000 -9.242) 9.242 3.891 ( 0.000 0.000 6.566) 6.566 4.239 ( -0.000 -0.000 2.390) 2.390 7.723 ( 0.000 0.000 2.930) 2.930 7.808 ( 0.000 0.000 1.896) 1.896 8.228 ( 0.000 0.000 -2.162) 2.162 8.261 ( -0.000 -0.000 -3.244) 3.244 8.688 ( 0.000 0.000 0.524) 0.524 8.982 ( 0.000 0.000 -2.008) 2.008 ======================= Grid point 227 (33/96) ======================= q-point: ( 0.07 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.697 ( 4.170 7.222 17.266) 19.174 0.702 ( 7.741 13.408 8.536) 17.679 1.339 ( 3.878 6.717 -1.096) 7.834 1.644 ( 16.774 29.054 16.035) 37.183 1.660 ( 8.354 14.470 -1.877) 16.813 3.878 ( -0.927 -1.605 -7.224) 7.458 7.228 ( 1.267 2.194 11.797) 12.066 7.319 ( 2.921 5.059 1.840) 6.124 7.669 ( 3.063 5.305 -1.441) 6.293 7.820 ( 6.242 10.812 -1.962) 12.638 8.692 ( 2.389 4.138 -3.007) 5.645 8.894 ( -1.066 -1.847 -0.455) 2.181 ======================= Grid point 228 (34/96) ======================= q-point: ( 0.14 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.894 ( 7.838 9.152 5.900) 13.417 0.967 ( 5.731 14.683 6.401) 17.013 1.472 ( 1.827 9.884 -0.831) 10.086 1.926 ( 8.785 6.257 -3.434) 11.319 2.246 ( 22.071 21.529 12.381) 33.226 3.869 ( 1.017 0.683 -5.366) 5.504 7.298 ( 3.432 3.576 10.442) 11.559 7.433 ( 2.653 7.487 2.049) 8.203 7.802 ( 5.167 6.278 -1.558) 8.278 8.052 ( 8.546 8.116 -2.445) 12.036 8.779 ( 2.501 2.823 -2.747) 4.666 8.858 ( -0.186 -1.672 -0.672) 1.811 ======================= Grid point 229 (35/96) ======================= q-point: ( 0.21 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.118 ( 6.795 8.827 0.389) 11.147 1.182 ( 3.246 13.031 4.420) 14.138 1.593 ( 0.209 10.153 -0.568) 10.171 2.060 ( 2.945 0.640 -5.196) 6.007 2.814 ( 21.580 15.609 9.989) 28.445 3.933 ( 4.198 2.630 -2.885) 5.732 7.413 ( 5.076 4.432 7.826) 10.327 7.552 ( 1.584 8.655 2.497) 9.146 7.957 ( 4.844 6.527 -1.760) 8.316 8.246 ( 6.425 3.366 -3.193) 7.925 8.784 ( -1.639 -0.406 -0.325) 1.719 8.899 ( 4.100 -0.037 -2.159) 4.634 ======================= Grid point 230 (36/96) ======================= q-point: ( 0.29 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.269 ( -1.110 11.312 2.347) 11.606 1.369 ( 6.806 7.471 -2.398) 10.387 1.670 ( -2.455 8.755 -0.492) 9.106 2.065 ( -1.299 -1.123 -6.603) 6.823 3.294 ( 18.239 10.838 8.377) 22.810 4.048 ( 5.453 2.491 -0.567) 6.022 7.555 ( 5.703 4.606 4.901) 8.818 7.658 ( 0.672 8.696 3.081) 9.250 8.088 ( 3.380 4.914 -1.996) 6.289 8.332 ( 1.606 0.668 -3.573) 3.974 8.752 ( -1.065 -1.224 0.281) 1.647 8.956 ( 2.702 -0.856 -2.025) 3.484 ======================= Grid point 231 (37/96) ======================= q-point: ( 0.36 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.327 ( -4.468 10.873 3.210) 12.185 1.551 ( 6.347 5.254 -6.342) 10.397 1.684 ( -4.871 7.484 -0.432) 8.940 2.014 ( -2.773 -1.395 -7.937) 8.522 3.652 ( 12.843 6.437 7.332) 16.129 4.154 ( 4.442 0.759 1.240) 4.674 7.688 ( 3.989 4.122 1.979) 6.067 7.750 ( 0.439 7.964 3.880) 8.870 8.162 ( 1.765 0.632 -2.476) 3.106 8.341 ( -1.587 2.544 -3.040) 4.270 8.717 ( -1.080 -2.159 0.449) 2.455 8.977 ( 0.998 -0.979 -1.968) 2.414 ======================= Grid point 232 (38/96) ======================= q-point: ( 0.43 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.345 ( -5.481 10.225 3.314) 12.065 1.659 ( -3.718 6.977 -1.075) 7.978 1.687 ( 1.310 3.725 -7.075) 8.102 1.955 ( -1.936 -1.073 -9.022) 9.289 3.853 ( 5.186 1.643 6.710) 8.638 4.209 ( 2.339 -1.374 2.347) 3.587 7.758 ( -0.411 3.929 1.131) 4.109 7.827 ( -0.134 5.440 3.382) 6.407 8.170 ( 1.249 -2.960 -3.007) 4.400 8.343 ( -2.679 5.080 -2.354) 6.207 8.681 ( -0.252 -1.846 0.606) 1.959 8.978 ( 0.234 -0.551 -1.993) 2.081 ======================= Grid point 242 (39/96) ======================= q-point: ( 0.14 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.088 ( 5.727 9.919 0.221) 11.456 1.254 ( 7.398 12.814 4.455) 15.452 1.666 ( 4.985 8.635 -0.828) 10.005 2.025 ( 2.035 3.525 -4.707) 6.223 2.689 ( 12.869 22.290 10.772) 27.901 3.908 ( 1.892 3.277 -3.357) 5.058 7.392 ( 3.265 5.655 8.038) 10.356 7.596 ( 4.728 8.188 1.879) 9.640 7.950 ( 4.675 8.097 -1.199) 9.426 8.195 ( 3.359 5.817 -2.944) 7.334 8.818 ( -0.969 -1.678 -0.265) 1.956 8.836 ( 1.336 2.314 -2.410) 3.598 ======================= Grid point 243 (40/96) ======================= q-point: ( 0.21 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.283 ( 6.207 8.000 -4.103) 10.925 1.459 ( 1.439 12.788 3.388) 13.308 1.801 ( -0.022 9.599 -0.488) 9.612 2.064 ( 0.513 -0.103 -5.982) 6.004 3.122 ( 14.993 15.144 9.216) 23.218 3.993 ( 3.720 3.183 -1.237) 5.050 7.519 ( 4.517 6.005 5.276) 9.182 7.744 ( 1.794 9.764 2.001) 10.127 8.102 ( 2.800 7.619 -1.189) 8.204 8.284 ( 3.667 0.590 -3.047) 4.804 8.774 ( -2.390 -0.708 -0.165) 2.498 8.888 ( 3.282 -0.657 -1.951) 3.875 ======================= Grid point 244 (41/96) ======================= q-point: ( 0.29 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.452 ( 4.089 6.980 -6.238) 10.215 1.575 ( -1.168 11.573 1.760) 11.764 1.857 ( -4.023 8.995 -0.462) 9.864 2.045 ( -1.903 -0.867 -7.013) 7.318 3.495 ( 13.445 9.427 7.987) 18.261 4.085 ( 4.012 1.130 0.686) 4.224 7.659 ( 4.599 5.746 2.459) 7.760 7.854 ( -0.411 10.068 2.373) 10.352 8.171 ( -1.051 3.086 -2.030) 3.840 8.355 ( 3.312 1.724 -2.142) 4.305 8.721 ( -2.240 -1.804 0.131) 2.879 8.933 ( 2.707 -0.961 -1.915) 3.452 ======================= Grid point 245 (42/96) ======================= q-point: ( 0.36 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.539 ( -2.764 8.618 -2.515) 9.394 1.683 ( 2.000 8.168 -3.954) 9.293 1.854 ( -5.945 8.349 -0.737) 10.276 1.989 ( -3.123 -1.091 -7.953) 8.614 3.760 ( 8.537 4.557 7.079) 11.990 4.144 ( 2.958 -1.640 2.213) 4.042 7.778 ( 2.534 4.744 -0.009) 5.378 7.936 ( -1.674 9.972 3.082) 10.571 8.149 ( -1.871 -1.167 -2.772) 3.543 8.420 ( 0.352 4.325 -1.542) 4.606 8.666 ( -1.157 -2.825 0.606) 3.112 8.961 ( 1.264 -0.296 -1.984) 2.371 ======================= Grid point 246 (43/96) ======================= q-point: ( 0.43 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.555 ( -5.175 8.963 -0.988) 10.396 1.761 ( -2.669 4.622 -5.688) 7.800 1.836 ( -4.851 8.401 -1.513) 9.818 1.944 ( -0.294 0.509 -8.128) 8.149 3.861 ( 0.303 -0.524 6.665) 6.693 4.159 ( 1.920 -3.326 2.878) 4.799 7.826 ( -1.627 2.818 -1.020) 3.410 7.981 ( -5.367 9.296 4.019) 11.462 8.115 ( 1.258 -2.179 -3.610) 4.400 8.446 ( -2.610 4.520 -1.495) 5.429 8.638 ( 1.375 -2.382 0.955) 2.912 8.971 ( 0.090 -0.155 -2.020) 2.027 ======================= Grid point 258 (44/96) ======================= q-point: ( 0.21 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.440 ( 4.450 7.708 -7.296) 11.509 1.660 ( 4.159 7.204 1.831) 8.518 1.950 ( 2.888 5.002 -0.124) 5.777 2.056 ( -0.470 -0.814 -6.668) 6.734 3.417 ( 8.198 14.199 8.390) 18.418 4.053 ( 1.535 2.659 0.429) 3.100 7.648 ( 3.880 6.720 2.600) 8.184 7.915 ( 3.978 6.890 0.961) 8.014 8.249 ( 3.459 5.991 -0.276) 6.924 8.279 ( -0.022 -0.037 -2.300) 2.300 8.743 ( -1.462 -2.532 -0.191) 2.930 8.893 ( 0.795 1.377 -1.942) 2.510 ======================= Grid point 259 (45/96) ======================= q-point: ( 0.29 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.584 ( 3.612 6.469 -9.242) 11.845 1.762 ( 0.449 6.629 0.379) 6.655 1.999 ( -2.857 4.112 -2.051) 5.411 2.029 ( -1.431 -0.281 -5.080) 5.286 3.676 ( 7.953 8.683 7.450) 13.933 4.092 ( 1.450 -0.510 2.064) 2.573 7.779 ( 3.509 6.009 0.024) 6.959 8.018 ( 1.235 5.084 -0.647) 5.272 8.235 ( -5.226 4.431 -0.250) 6.856 8.388 ( 6.133 1.211 -0.099) 6.252 8.678 ( -2.804 -2.732 0.097) 3.916 8.935 ( 1.954 1.330 -2.053) 3.131 ======================= Grid point 260 (46/96) ======================= q-point: ( 0.36 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.674 ( -0.693 5.703 -8.115) 9.942 1.838 ( 1.425 6.048 -2.152) 6.576 1.939 ( -4.342 -0.961 -7.475) 8.698 2.016 ( -3.063 4.608 -0.303) 5.541 3.845 ( 3.329 4.141 6.613) 8.483 4.094 ( 1.316 -3.193 3.326) 4.795 7.873 ( 0.583 4.568 -1.860) 4.966 8.068 ( 1.516 0.324 -3.142) 3.504 8.190 ( -6.549 7.811 2.449) 10.483 8.487 ( 3.116 2.238 -0.180) 3.840 8.602 ( -1.700 -3.676 0.883) 4.145 8.972 ( 0.505 1.238 -2.121) 2.507 ======================= Grid point 273 (47/96) ======================= q-point: ( 0.29 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.709 ( 3.270 5.664 -10.825) 12.647 1.848 ( 1.401 2.427 -2.813) 3.971 1.980 ( -1.599 -2.770 -5.237) 6.136 2.069 ( 0.852 1.476 0.236) 1.720 3.836 ( 4.044 7.004 6.570) 10.420 4.064 ( -1.254 -2.173 3.513) 4.316 7.892 ( 2.901 5.024 -2.132) 6.181 8.054 ( -0.300 -0.519 -3.691) 3.739 8.336 ( 2.088 3.616 2.476) 4.855 8.417 ( 1.750 3.031 0.767) 3.582 8.608 ( -2.412 -4.177 0.486) 4.848 8.974 ( 1.169 2.026 -2.192) 3.206 ======================= Grid point 274 (48/96) ======================= q-point: ( 0.36 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.779 ( -2.591 4.487 -10.407) 11.626 1.875 ( 1.410 -2.442 -10.924) 11.282 1.934 ( -0.723 1.252 1.482) 2.070 2.082 ( -0.923 1.598 0.297) 1.869 3.925 ( -2.093 3.625 5.812) 7.163 4.025 ( 1.934 -3.350 4.466) 5.908 7.955 ( -1.985 3.438 -3.251) 5.131 8.035 ( 1.353 -2.344 -4.106) 4.918 8.360 ( -4.409 7.637 2.926) 9.291 8.513 ( 1.997 -3.458 1.586) 4.297 8.513 ( 0.698 -1.209 0.276) 1.423 8.997 ( -0.506 0.877 -2.221) 2.441 ======================= Grid point 422 (49/96) ======================= q-point: ( 0.00 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.688 ( 0.000 -0.000 13.522) 13.522 0.688 ( 0.000 0.000 13.522) 13.522 1.183 ( 0.000 -0.000 -2.694) 2.694 1.183 ( 0.000 0.000 -2.694) 2.694 1.846 ( 0.000 -0.000 39.178) 39.178 3.665 ( -0.000 -0.000 -16.761) 16.761 7.272 ( -0.000 0.000 3.165) 3.165 7.272 ( 0.000 0.000 3.165) 3.165 7.507 ( 0.000 -0.000 -3.005) 3.005 7.507 ( 0.000 0.000 -3.005) 3.005 8.690 ( 0.000 -0.000 4.579) 4.579 8.919 ( -0.000 0.000 -1.396) 1.396 ======================= Grid point 423 (50/96) ======================= q-point: ( 0.07 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.760 ( 5.857 3.381 12.202) 13.951 0.967 ( 14.995 8.657 11.004) 20.515 1.219 ( 2.974 1.717 -2.273) 4.118 1.369 ( 13.699 7.909 -2.619) 16.033 1.878 ( 6.478 3.740 36.758) 37.511 3.652 ( -0.872 -0.504 -16.202) 16.233 7.304 ( 2.648 1.529 3.255) 4.466 7.355 ( 5.897 3.405 4.901) 8.389 7.546 ( 3.296 1.903 -3.106) 4.912 7.605 ( 7.954 4.592 -3.043) 9.675 8.669 ( -0.874 -0.505 3.784) 3.916 8.902 ( -1.388 -0.802 -1.544) 2.226 ======================= Grid point 424 (51/96) ======================= q-point: ( 0.14 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.930 ( 8.217 4.744 9.867) 13.689 1.152 ( 2.819 1.628 17.993) 18.285 1.313 ( 4.929 2.846 -1.624) 5.919 1.700 ( 12.768 7.372 -4.129) 15.311 2.241 ( 22.485 12.982 22.610) 34.428 3.652 ( 1.546 0.892 -13.946) 14.060 7.388 ( 4.436 2.561 3.555) 6.234 7.487 ( 4.974 2.872 9.426) 11.038 7.651 ( 5.567 3.214 -3.356) 7.251 7.835 ( 11.053 6.381 -3.431) 13.216 8.687 ( 2.617 1.511 0.430) 3.052 8.863 ( -1.750 -1.010 -1.669) 2.620 ======================= Grid point 425 (52/96) ======================= q-point: ( 0.21 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.112 ( 7.050 4.070 8.075) 11.466 1.206 ( 2.665 1.539 14.042) 14.375 1.433 ( 5.132 2.963 -1.427) 6.095 1.898 ( 4.160 2.402 -6.857) 8.372 2.788 ( 23.420 13.521 15.323) 31.083 3.742 ( 6.237 3.601 -9.697) 12.079 7.500 ( 5.083 2.934 4.039) 7.124 7.584 ( 3.681 2.125 11.461) 12.224 7.792 ( 6.383 3.685 -3.643) 8.221 8.071 ( 8.462 4.885 -4.723) 10.852 8.759 ( 2.278 1.315 -1.118) 2.859 8.842 ( 0.785 0.453 -1.904) 2.109 ======================= Grid point 426 (53/96) ======================= q-point: ( 0.29 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.244 ( 4.309 2.488 6.995) 8.584 1.279 ( 3.505 2.024 5.800) 7.073 1.534 ( 3.359 1.939 -1.499) 4.158 1.914 ( -1.956 -1.129 -9.284) 9.555 3.291 ( 19.828 11.448 12.005) 25.852 3.919 ( 8.448 4.878 -5.156) 11.034 7.615 ( 4.799 2.771 4.550) 7.169 7.668 ( 3.724 2.150 10.048) 10.930 7.935 ( 5.862 3.384 -3.879) 7.801 8.203 ( 2.913 1.682 -6.077) 6.946 8.774 ( -0.129 -0.075 0.419) 0.445 8.894 ( 2.415 1.394 -2.962) 4.068 ======================= Grid point 427 (54/96) ======================= q-point: ( 0.36 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.314 ( 1.909 1.102 6.389) 6.759 1.362 ( 3.676 2.122 -1.722) 4.581 1.581 ( 0.733 0.423 -1.368) 1.609 1.840 ( -3.977 -2.296 -11.080) 11.994 3.693 ( 14.751 8.516 9.956) 19.729 4.106 ( 7.273 4.199 -1.628) 8.554 7.716 ( 3.791 2.189 4.907) 6.576 7.755 ( 3.713 2.144 7.239) 8.414 8.055 ( 4.380 2.529 -4.030) 6.467 8.221 ( -0.781 -0.451 -6.152) 6.218 8.761 ( -1.070 -0.618 1.107) 1.659 8.928 ( 0.552 0.319 -3.242) 3.304 ======================= Grid point 428 (55/96) ======================= q-point: ( 0.43 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.341 ( 0.618 0.357 6.006) 6.048 1.445 ( 3.275 1.891 -7.121) 8.063 1.578 ( -0.671 -0.387 -0.920) 1.203 1.747 ( -3.684 -2.127 -12.374) 13.085 3.961 ( 8.156 4.709 8.632) 12.775 4.239 ( 4.050 2.338 0.691) 4.727 7.785 ( 2.127 1.228 5.062) 5.626 7.832 ( 2.682 1.549 4.593) 5.539 8.133 ( 2.323 1.341 -4.105) 4.904 8.190 ( -1.487 -0.859 -5.437) 5.702 8.729 ( -1.430 -0.825 1.464) 2.206 8.928 ( -0.376 -0.217 -3.471) 3.498 ======================= Grid point 429 (56/96) ======================= q-point: (-0.50 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.347 ( 0.000 0.000 5.836) 5.836 1.491 ( 0.000 0.000 -9.183) 9.183 1.567 ( 0.000 0.000 -0.631) 0.631 1.696 ( -0.000 -0.000 -12.837) 12.837 4.057 ( 0.000 0.000 8.117) 8.117 4.287 ( 0.000 0.000 1.556) 1.556 7.811 ( 0.000 0.000 5.094) 5.094 7.866 ( 0.000 0.000 3.442) 3.442 8.160 ( 0.000 0.000 -4.127) 4.127 8.170 ( -0.000 -0.000 -5.010) 5.010 8.710 ( 0.000 0.000 1.647) 1.647 8.922 ( 0.000 0.000 -3.579) 3.579 ======================= Grid point 438 (57/96) ======================= q-point: ( 0.07 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.911 ( 5.939 10.287 10.146) 15.622 1.114 ( 2.333 4.040 16.394) 17.045 1.308 ( 3.786 6.557 -1.826) 7.788 1.602 ( 7.575 13.120 -3.448) 15.537 2.101 ( 11.086 19.201 26.578) 34.611 3.645 ( 0.228 0.394 -14.781) 14.788 7.374 ( 2.959 5.125 3.315) 6.783 7.444 ( 2.822 4.888 7.784) 9.615 7.624 ( 3.238 5.608 -2.654) 6.998 7.762 ( 5.965 10.332 -3.368) 12.397 8.672 ( 0.852 1.475 1.455) 2.240 8.873 ( -1.094 -1.894 -1.662) 2.747 ======================= Grid point 439 (58/96) ======================= q-point: ( 0.14 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.113 ( 2.971 11.143 9.151) 14.722 1.176 ( 3.631 2.376 15.006) 15.621 1.447 ( 2.552 9.482 -1.661) 9.959 1.827 ( 6.712 4.774 -5.622) 9.972 2.570 ( 18.625 18.877 17.649) 31.855 3.691 ( 3.476 3.370 -11.438) 12.421 7.488 ( 2.603 6.941 3.937) 8.393 7.541 ( 4.054 2.693 9.857) 10.993 7.757 ( 4.394 7.065 -2.485) 8.683 7.980 ( 7.913 7.368 -4.123) 11.572 8.729 ( 2.780 2.576 -1.167) 3.965 8.836 ( -0.746 -1.535 -1.507) 2.277 ======================= Grid point 440 (59/96) ======================= q-point: ( 0.21 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.208 ( 2.800 2.494 8.190) 9.007 1.313 ( 1.428 10.988 7.172) 13.198 1.573 ( 0.369 10.065 -1.395) 10.168 1.914 ( 1.153 -0.424 -8.169) 8.261 3.072 ( 19.267 14.415 13.445) 27.564 3.830 ( 7.085 5.052 -7.017) 11.179 7.590 ( 2.697 4.373 6.376) 8.188 7.647 ( 3.355 5.587 7.565) 9.985 7.905 ( 3.985 7.328 -2.913) 8.835 8.153 ( 5.415 3.173 -5.287) 8.207 8.773 ( -0.037 0.167 -0.496) 0.524 8.848 ( 3.241 -0.087 -2.384) 4.024 ======================= Grid point 441 (60/96) ======================= q-point: ( 0.29 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.284 ( 3.483 1.649 1.472) 4.125 1.416 ( -1.756 10.447 5.043) 11.733 1.650 ( -2.417 8.682 -1.483) 9.133 1.883 ( -2.931 -1.533 -9.861) 10.401 3.509 ( 16.621 10.147 10.944) 22.338 4.013 ( 7.905 4.286 -3.088) 9.507 7.676 ( 3.265 2.639 5.823) 7.179 7.759 ( 1.999 7.132 6.222) 9.673 8.028 ( 2.629 5.207 -3.452) 6.778 8.231 ( 1.212 1.762 -5.702) 6.090 8.762 ( -0.402 -1.243 0.707) 1.485 8.900 ( 2.529 -0.706 -3.016) 3.998 ======================= Grid point 442 (61/96) ======================= q-point: ( 0.36 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.355 ( 2.385 1.952 -2.426) 3.922 1.474 ( -2.375 9.079 2.034) 9.602 1.665 ( -4.541 6.813 -1.884) 8.402 1.801 ( -4.535 -0.987 -10.699) 11.662 3.838 ( 11.748 5.986 9.253) 16.108 4.170 ( 6.139 2.065 -0.176) 6.479 7.752 ( 2.796 1.763 3.948) 5.149 7.855 ( 0.776 7.070 5.585) 9.043 8.088 ( 1.163 0.097 -4.176) 4.336 8.255 ( -1.465 4.480 -4.853) 6.765 8.734 ( -0.790 -1.974 1.317) 2.501 8.919 ( 0.851 -0.999 -3.336) 3.585 ======================= Grid point 443 (62/96) ======================= q-point: ( 0.43 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.396 ( -0.585 2.938 -1.529) 3.363 1.526 ( -1.600 7.562 -2.628) 8.164 1.641 ( -4.245 4.429 -3.852) 7.244 1.726 ( -3.684 1.458 -9.560) 10.349 4.022 ( 4.722 1.448 8.253) 9.618 4.256 ( 2.894 -0.630 1.507) 3.323 7.800 ( 1.080 0.716 2.779) 3.067 7.921 ( -1.473 6.656 5.156) 8.547 8.084 ( 1.155 -3.466 -4.690) 5.945 8.273 ( -2.866 6.151 -4.174) 7.966 8.705 ( -0.219 -1.485 1.731) 2.291 8.918 ( 0.172 -0.587 -3.555) 3.607 ======================= Grid point 453 (63/96) ======================= q-point: ( 0.14 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.196 ( 1.274 2.206 9.679) 10.009 1.371 ( 6.318 10.943 5.970) 13.975 1.643 ( 5.261 9.112 -1.335) 10.606 1.890 ( 0.953 1.650 -7.659) 7.892 2.965 ( 11.595 20.084 14.405) 27.300 3.791 ( 3.738 6.474 -7.837) 10.831 7.590 ( 1.783 3.088 9.356) 10.012 7.652 ( 4.617 7.997 3.141) 9.754 7.918 ( 4.921 8.524 -1.630) 9.977 8.110 ( 3.085 5.344 -4.799) 7.817 8.783 ( 1.423 2.465 -2.087) 3.530 8.805 ( -0.773 -1.338 -0.941) 1.810 ======================= Grid point 454 (64/96) ======================= q-point: ( 0.21 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.250 ( 2.368 2.168 2.202) 3.893 1.541 ( 1.116 9.874 3.595) 10.567 1.785 ( -0.332 9.939 -1.209) 10.018 1.901 ( -0.136 -0.368 -8.379) 8.388 3.357 ( 13.484 14.027 11.988) 22.854 3.938 ( 5.779 5.613 -4.183) 9.078 7.656 ( 2.349 3.009 6.846) 7.838 7.800 ( 2.930 8.003 2.914) 9.007 8.069 ( 1.492 8.704 -1.639) 8.982 8.200 ( 3.956 1.799 -4.464) 6.230 8.768 ( -1.727 -0.703 -0.293) 1.888 8.838 ( 3.097 -0.692 -2.526) 4.056 ======================= Grid point 455 (65/96) ======================= q-point: ( 0.29 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.319 ( 3.119 1.995 -4.179) 5.584 1.617 ( -1.901 7.500 0.093) 7.738 1.821 ( -4.733 6.208 -4.063) 8.801 1.889 ( -3.124 4.135 -4.495) 6.860 3.697 ( 12.141 8.744 10.062) 18.031 4.086 ( 5.841 2.993 -1.029) 6.644 7.728 ( 2.562 2.731 3.817) 5.347 7.910 ( 1.522 6.373 1.989) 6.847 8.123 ( -2.452 4.882 -1.638) 5.704 8.296 ( 2.803 4.649 -3.185) 6.294 8.728 ( -1.296 -1.998 0.651) 2.469 8.877 ( 2.478 -1.109 -3.118) 4.134 ======================= Grid point 456 (66/96) ======================= q-point: ( 0.36 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.391 ( 2.581 1.848 -8.214) 8.806 1.642 ( -1.473 4.100 -5.648) 7.133 1.771 ( -6.699 5.281 -3.606) 9.261 1.873 ( -5.597 7.797 -1.407) 9.700 3.936 ( 7.679 4.015 8.614) 12.218 4.188 ( 4.046 -0.197 1.328) 4.263 7.791 ( 1.619 2.108 1.328) 2.971 7.983 ( 1.645 2.977 -0.641) 3.461 8.107 ( -4.439 4.695 0.989) 6.537 8.371 ( -0.617 6.081 -3.056) 6.833 8.689 ( -0.416 -2.466 1.678) 3.011 8.901 ( 1.178 -0.500 -3.540) 3.764 ======================= Grid point 457 (67/96) ======================= q-point: ( 0.43 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.427 ( -0.400 0.692 -9.169) 9.204 1.653 ( 0.402 -0.696 -13.434) 13.458 1.730 ( -6.039 10.460 2.277) 12.290 1.861 ( -5.186 8.983 -0.849) 10.408 4.026 ( 0.393 -0.681 7.975) 8.013 4.221 ( 1.512 -2.619 2.305) 3.803 7.817 ( -0.548 0.949 0.365) 1.155 8.019 ( 1.485 -2.572 -4.995) 5.811 8.084 ( -5.129 8.883 5.396) 11.590 8.396 ( -3.242 5.615 -3.203) 7.232 8.672 ( 1.073 -1.859 2.229) 3.094 8.909 ( 0.167 -0.290 -3.690) 3.705 ======================= Grid point 469 (68/96) ======================= q-point: ( 0.21 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.301 ( 1.717 2.974 -4.001) 5.273 1.668 ( 1.568 2.716 -3.302) 4.554 1.906 ( 0.623 1.080 -4.848) 5.005 1.944 ( 3.067 5.312 -0.477) 6.152 3.629 ( 7.459 12.919 10.520) 18.254 4.047 ( 2.917 5.052 -1.432) 6.006 7.720 ( 1.939 3.358 3.879) 5.485 7.926 ( 2.424 4.199 -0.415) 4.867 8.231 ( 1.639 2.839 -1.574) 3.637 8.243 ( 3.948 6.838 -0.254) 7.899 8.738 ( -1.351 -2.340 -0.193) 2.708 8.838 ( 0.592 1.025 -2.966) 3.193 ======================= Grid point 470 (69/96) ======================= q-point: ( 0.29 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.368 ( 2.513 2.995 -9.028) 9.838 1.674 ( -1.129 -0.753 -10.465) 10.552 1.917 ( -2.628 4.276 0.128) 5.020 2.008 ( -1.456 5.109 -0.320) 5.322 3.861 ( 7.035 7.627 8.984) 13.725 4.131 ( 2.756 1.483 1.106) 3.319 7.786 ( 1.893 3.002 0.834) 3.646 7.972 ( 0.732 0.386 -3.521) 3.617 8.258 ( -3.704 8.055 2.385) 9.181 8.384 ( 4.929 3.408 -0.354) 6.003 8.683 ( -1.838 -2.507 0.660) 3.178 8.873 ( 1.847 0.899 -3.638) 4.178 ======================= Grid point 471 (70/96) ======================= q-point: ( 0.36 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.436 ( 1.253 2.708 -11.820) 12.191 1.637 ( -0.631 -2.318 -14.220) 14.422 1.908 ( -4.015 6.076 2.477) 7.693 2.015 ( -3.331 5.552 -0.175) 6.477 4.007 ( 2.965 3.187 7.702) 8.847 4.168 ( 1.678 -1.631 2.956) 3.770 7.836 ( 0.400 2.309 -1.231) 2.647 7.975 ( 0.844 -1.830 -4.877) 5.277 8.267 ( -4.989 8.939 4.250) 11.084 8.472 ( 0.817 3.674 -1.435) 4.029 8.635 ( -0.138 -2.980 2.392) 3.823 8.906 ( 0.582 0.929 -3.981) 4.129 ======================= Grid point 484 (71/96) ======================= q-point: ( 0.29 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.439 ( 2.231 3.864 -12.668) 13.431 1.631 ( -1.772 -3.070 -14.352) 14.783 1.975 ( 0.859 1.488 2.427) 2.973 2.073 ( 0.941 1.629 -0.039) 1.882 3.996 ( 3.256 5.639 7.569) 9.984 4.143 ( -0.183 -0.317 3.214) 3.234 7.845 ( 1.580 2.736 -1.830) 3.651 7.958 ( -0.795 -1.378 -4.815) 5.071 8.403 ( 2.713 4.698 3.264) 6.331 8.435 ( 1.806 3.129 0.787) 3.697 8.626 ( -1.883 -3.262 1.512) 4.059 8.905 ( 1.060 1.835 -4.188) 4.693 ======================= Grid point 485 (72/96) ======================= q-point: ( 0.36 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.496 ( -1.716 2.972 -14.170) 14.580 1.581 ( 1.630 -2.823 -15.202) 15.547 1.987 ( -1.092 1.891 2.910) 3.638 2.087 ( -0.980 1.697 0.045) 1.960 4.065 ( -1.427 2.472 6.458) 7.061 4.129 ( 1.158 -2.005 4.557) 5.111 7.881 ( -1.185 2.052 -3.199) 3.981 7.935 ( 1.112 -1.926 -4.663) 5.166 8.428 ( -3.960 6.859 2.834) 8.412 8.511 ( 0.117 -0.203 -0.498) 0.550 8.569 ( 1.826 -3.163 3.985) 5.406 8.927 ( -0.443 0.768 -4.341) 4.431 ======================= Grid point 633 (73/96) ======================= q-point: ( 0.00 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.934 ( 0.000 0.000 9.587) 9.587 0.934 ( 0.000 -0.000 9.587) 9.587 1.096 ( -0.000 0.000 -5.696) 5.696 1.096 ( -0.000 0.000 -5.696) 5.696 2.611 ( -0.000 -0.000 32.500) 32.500 3.222 ( 0.000 -0.000 -24.824) 24.824 7.349 ( 0.000 0.000 3.905) 3.905 7.349 ( 0.000 -0.000 3.905) 3.905 7.433 ( 0.000 -0.000 -3.833) 3.833 7.433 ( -0.000 0.000 -3.833) 3.833 8.789 ( -0.000 0.000 4.505) 4.505 8.872 ( -0.000 -0.000 -3.096) 3.096 ======================= Grid point 634 (74/96) ======================= q-point: ( 0.07 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.988 ( 4.435 2.560 9.025) 10.377 1.140 ( 3.642 2.103 -5.292) 6.759 1.160 ( 15.113 8.726 7.471) 18.983 1.290 ( 13.949 8.054 -4.847) 16.821 2.619 ( 1.931 1.115 31.948) 32.026 3.220 ( 0.345 0.199 -24.316) 24.319 7.383 ( 2.853 1.647 4.032) 5.207 7.447 ( 7.630 4.405 4.047) 9.696 7.469 ( 3.088 1.783 -3.963) 5.332 7.531 ( 7.877 4.548 -3.804) 9.859 8.764 ( -1.791 -1.034 4.630) 5.071 8.851 ( -1.646 -0.950 -3.327) 3.832 ======================= Grid point 635 (75/96) ======================= q-point: ( 0.14 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.121 ( 6.586 3.802 7.709) 10.828 1.251 ( 5.539 3.198 -4.413) 7.771 1.441 ( 7.220 4.168 9.684) 12.778 1.598 ( 10.495 6.059 -5.668) 13.379 2.772 ( 12.028 6.944 25.652) 29.171 3.278 ( 5.446 3.144 -21.058) 21.977 7.474 ( 4.843 2.796 4.390) 7.110 7.568 ( 5.246 3.029 -4.299) 7.428 7.645 ( 8.363 4.828 5.816) 11.273 7.752 ( 10.136 5.852 -4.402) 12.505 8.735 ( -0.192 -0.111 3.460) 3.468 8.810 ( -1.460 -0.843 -3.225) 3.639 ======================= Grid point 636 (76/96) ======================= q-point: ( 0.21 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.271 ( 5.943 3.431 6.477) 9.437 1.380 ( 5.303 3.062 -3.738) 7.174 1.504 ( -0.450 -0.260 12.721) 12.732 1.732 ( 1.412 0.815 -8.919) 9.066 3.136 ( 17.939 10.357 16.901) 26.734 3.480 ( 11.476 6.626 -14.700) 19.791 7.597 ( 5.595 3.230 4.907) 8.113 7.701 ( 6.056 3.497 -4.729) 8.442 7.799 ( 4.803 2.773 8.576) 10.213 7.954 ( 6.777 3.913 -6.233) 10.005 8.759 ( 1.835 1.059 1.172) 2.422 8.798 ( 0.658 0.380 -2.183) 2.311 ======================= Grid point 637 (77/96) ======================= q-point: ( 0.29 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.383 ( 3.705 2.139 5.786) 7.196 1.466 ( -2.356 -1.360 10.703) 11.044 1.483 ( 3.408 1.967 -3.489) 5.259 1.699 ( -3.479 -2.009 -10.764) 11.489 3.546 ( 16.815 9.708 11.540) 22.587 3.768 ( 12.640 7.297 -8.980) 17.136 7.724 ( 5.235 3.022 5.429) 8.125 7.838 ( 5.613 3.241 -5.128) 8.265 7.877 ( 2.282 1.317 9.221) 9.590 8.056 ( 2.269 1.310 -7.636) 8.073 8.792 ( 0.903 0.521 1.396) 1.743 8.835 ( 1.759 1.015 -2.495) 3.217 ======================= Grid point 638 (78/96) ======================= q-point: ( 0.36 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.411 ( -2.172 -1.254 6.178) 6.668 1.442 ( 1.477 0.853 5.336) 5.602 1.534 ( 1.073 0.620 -3.224) 3.454 1.600 ( -4.674 -2.699 -10.804) 12.077 3.892 ( 12.827 7.405 8.319) 16.988 4.036 ( 10.153 5.862 -5.047) 12.764 7.832 ( 4.031 2.327 5.792) 7.430 7.917 ( 1.338 0.772 7.695) 7.849 7.953 ( 4.241 2.448 -5.399) 7.289 8.077 ( -0.088 -0.051 -7.246) 7.247 8.797 ( -0.522 -0.302 2.285) 2.364 8.857 ( 0.132 0.076 -3.252) 3.256 ======================= Grid point 639 (79/96) ======================= q-point: ( 0.43 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.374 ( -0.892 -0.515 1.289) 1.650 1.460 ( 0.268 0.155 4.856) 4.866 1.500 ( -3.657 -2.111 -9.799) 10.670 1.541 ( -0.192 -0.111 -2.711) 2.720 4.125 ( 7.076 4.085 6.356) 10.352 4.221 ( 5.636 3.254 -2.610) 7.011 7.905 ( 2.201 1.271 5.953) 6.472 7.942 ( 0.818 0.472 5.609) 5.688 8.029 ( 2.271 1.311 -5.524) 6.114 8.066 ( -0.591 -0.341 -5.972) 6.010 8.775 ( -1.065 -0.615 2.871) 3.124 8.848 ( -0.660 -0.381 -3.782) 3.858 ======================= Grid point 640 (80/96) ======================= q-point: (-0.50 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.367 ( 0.000 0.000 -1.328) 1.328 1.452 ( -0.000 -0.000 -8.797) 8.797 1.462 ( 0.000 0.000 4.596) 4.596 1.536 ( 0.000 0.000 -2.398) 2.398 4.208 ( 0.000 0.000 5.641) 5.641 4.287 ( 0.000 0.000 -1.733) 1.733 7.931 ( 0.000 0.000 5.988) 5.988 7.952 ( 0.000 0.000 4.595) 4.595 8.056 ( 0.000 0.000 -5.554) 5.554 8.058 ( -0.000 -0.000 -5.260) 5.260 8.761 ( 0.000 0.000 3.112) 3.112 8.838 ( 0.000 0.000 -3.994) 3.994 ======================= Grid point 649 (81/96) ======================= q-point: ( 0.07 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.107 ( 4.780 8.280 7.829) 12.357 1.236 ( 4.110 7.119 -4.160) 9.213 1.375 ( 5.401 9.355 8.115) 13.510 1.512 ( 6.747 11.686 -5.119) 14.433 2.704 ( 5.019 8.693 28.160) 29.896 3.248 ( 1.995 3.455 -22.367) 22.720 7.454 ( 2.973 5.150 4.011) 7.173 7.531 ( 2.902 5.026 -1.966) 6.127 7.597 ( 5.003 8.666 3.054) 10.462 7.682 ( 5.635 9.761 -4.035) 11.972 8.737 ( -0.594 -1.028 4.020) 4.192 8.820 ( -1.047 -1.813 -3.369) 3.967 ======================= Grid point 650 (82/96) ======================= q-point: ( 0.14 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.279 ( 3.145 10.718 6.557) 12.952 1.387 ( 2.798 9.296 -3.452) 10.303 1.484 ( 2.373 0.338 11.686) 11.929 1.688 ( 4.488 2.409 -7.705) 9.237 2.979 ( 12.893 13.328 19.993) 27.268 3.384 ( 7.378 7.575 -17.239) 20.223 7.573 ( 3.495 6.765 4.142) 8.668 7.650 ( 4.708 4.821 -1.977) 7.022 7.762 ( 4.649 6.530 5.294) 9.606 7.881 ( 6.837 6.494 -5.175) 10.756 8.740 ( 1.134 0.710 1.904) 2.327 8.792 ( -0.271 -0.848 -2.559) 2.710 ======================= Grid point 651 (83/96) ======================= q-point: ( 0.21 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.395 ( -0.145 5.013 7.261) 8.825 1.467 ( 0.301 2.969 4.462) 5.368 1.545 ( -0.004 8.818 1.818) 9.003 1.722 ( -0.504 -0.950 -9.708) 9.768 3.364 ( 15.709 11.991 13.625) 24.004 3.634 ( 11.210 8.466 -11.250) 17.997 7.698 ( 3.096 6.154 4.640) 8.306 7.774 ( 4.695 3.048 -0.463) 5.617 7.886 ( 2.313 6.578 4.675) 8.395 8.026 ( 4.037 3.373 -6.600) 8.440 8.773 ( 1.379 0.250 0.695) 1.564 8.802 ( 2.052 -0.138 -1.875) 2.783 ======================= Grid point 652 (84/96) ======================= q-point: ( 0.29 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.384 ( -0.517 -2.829 7.316) 7.861 1.503 ( -1.368 3.949 0.506) 4.210 1.606 ( -3.449 8.401 -1.252) 9.168 1.684 ( -4.191 3.118 -7.427) 9.080 3.733 ( 14.255 8.772 9.695) 19.343 3.909 ( 11.127 6.631 -6.662) 14.566 7.794 ( 2.762 2.637 5.216) 6.464 7.870 ( 3.734 1.755 0.717) 4.187 7.978 ( 0.289 7.103 1.912) 7.362 8.098 ( 0.701 3.743 -6.620) 7.638 8.787 ( 0.689 -1.079 1.743) 2.163 8.836 ( 1.793 -0.836 -2.788) 3.419 ======================= Grid point 653 (85/96) ======================= q-point: ( 0.36 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.349 ( 0.071 -2.980 3.454) 4.562 1.484 ( -1.680 -1.889 -6.474) 6.950 1.599 ( -6.048 8.237 0.982) 10.266 1.659 ( -5.614 8.123 -3.242) 10.393 4.018 ( 10.214 5.076 7.149) 13.461 4.132 ( 8.241 3.674 -3.526) 9.688 7.849 ( 2.491 -0.458 4.912) 5.526 7.938 ( 2.882 0.109 -1.277) 3.154 8.032 ( -1.644 5.956 2.344) 6.608 8.141 ( -1.492 6.398 -5.747) 8.728 8.776 ( -0.209 -1.622 2.626) 3.093 8.843 ( 0.418 -1.236 -3.547) 3.780 ======================= Grid point 654 (86/96) ======================= q-point: ( 0.43 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.333 ( 1.035 -2.261 -0.177) 2.493 1.434 ( -0.578 -2.671 -8.971) 9.378 1.582 ( -5.937 9.383 3.307) 11.586 1.644 ( -5.509 8.558 -2.334) 10.442 4.175 ( 4.173 0.985 5.736) 7.161 4.256 ( 3.546 0.247 -1.796) 3.982 7.874 ( 1.801 -1.822 4.135) 4.864 7.970 ( 2.409 -2.711 -4.399) 5.701 8.049 ( -3.347 6.586 5.323) 9.106 8.169 ( -3.060 6.898 -5.435) 9.300 8.757 ( -0.106 -1.076 3.151) 3.331 8.835 ( 0.043 -0.756 -3.988) 4.059 ======================= Grid point 664 (87/96) ======================= q-point: ( 0.14 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.432 ( -1.431 -2.478 10.357) 10.745 1.485 ( 4.531 7.847 4.015) 9.911 1.610 ( 5.521 9.563 -1.226) 11.110 1.709 ( 0.458 0.794 -8.987) 9.034 3.281 ( 9.353 16.201 14.784) 23.843 3.574 ( 6.411 11.104 -12.340) 17.795 7.720 ( 4.058 7.028 2.869) 8.607 7.736 ( 2.052 3.553 2.930) 5.042 7.907 ( 4.519 7.828 1.925) 9.241 7.999 ( 3.079 5.334 -5.373) 8.173 8.761 ( 0.732 1.269 0.089) 1.468 8.780 ( -0.211 -0.366 -1.306) 1.373 ======================= Grid point 665 (88/96) ======================= q-point: ( 0.21 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.374 ( -1.550 -3.081 8.498) 9.172 1.547 ( -1.298 -0.081 -5.217) 5.377 1.746 ( -0.562 8.950 -0.778) 9.002 1.776 ( 1.313 7.842 -2.701) 8.398 3.605 ( 11.335 11.949 10.928) 19.765 3.809 ( 8.556 8.862 -7.954) 14.663 7.785 ( 1.130 1.778 4.825) 5.265 7.841 ( 2.926 3.755 -0.009) 4.760 8.052 ( 1.444 9.206 0.731) 9.347 8.107 ( 3.758 5.272 -3.960) 7.589 8.768 ( -0.112 -0.789 0.421) 0.901 8.791 ( 1.831 -0.901 -1.732) 2.676 ======================= Grid point 666 (89/96) ======================= q-point: ( 0.29 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.324 ( -0.653 -2.920 4.807) 5.662 1.493 ( -2.195 -3.179 -9.559) 10.310 1.796 ( -4.445 9.722 1.863) 10.851 1.840 ( -3.746 9.245 -1.884) 10.151 3.894 ( 10.376 7.387 8.055) 15.070 4.029 ( 8.157 5.321 -4.479) 10.720 7.816 ( 1.343 0.269 4.246) 4.462 7.898 ( 1.662 0.287 -2.996) 3.438 8.140 ( -1.372 8.990 3.107) 9.610 8.221 ( 1.102 7.718 -3.886) 8.711 8.754 ( 0.242 -1.966 1.923) 2.761 8.809 ( 1.670 -1.546 -3.082) 3.831 ======================= Grid point 667 (90/96) ======================= q-point: ( 0.36 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.296 ( 0.489 -2.215 0.651) 2.359 1.423 ( -1.609 -3.458 -10.234) 10.922 1.798 ( -6.300 10.007 2.667) 12.122 1.845 ( -5.907 9.360 -1.672) 11.193 4.097 ( 6.613 2.991 6.022) 9.431 4.183 ( 5.356 1.508 -2.018) 5.919 7.837 ( 1.239 -0.495 3.002) 3.285 7.915 ( 1.369 -1.700 -4.113) 4.656 8.188 ( -3.341 8.735 4.968) 10.590 8.290 ( -1.960 7.737 -4.434) 9.130 8.740 ( 0.395 -1.803 3.143) 3.645 8.818 ( 0.916 -1.029 -4.001) 4.231 ======================= Grid point 668 (91/96) ======================= q-point: ( 0.43 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.292 ( 1.017 -1.761 -1.588) 2.580 1.387 ( 1.162 -2.013 -9.836) 10.107 1.791 ( -6.016 10.421 2.824) 12.360 1.838 ( -5.659 9.802 -1.558) 11.426 4.171 ( 0.645 -1.118 5.181) 5.339 4.238 ( 1.041 -1.803 -1.023) 2.320 7.846 ( 0.599 -1.037 2.385) 2.669 7.918 ( 1.317 -2.281 -4.259) 5.008 8.202 ( -4.673 8.094 5.503) 10.846 8.311 ( -3.976 6.887 -4.647) 9.211 8.734 ( 0.712 -1.233 3.624) 3.894 8.821 ( 0.373 -0.646 -4.319) 4.383 ======================= Grid point 680 (92/96) ======================= q-point: ( 0.21 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.315 ( -1.521 -2.635 5.369) 6.171 1.496 ( -2.120 -3.672 -10.304) 11.143 1.892 ( 3.068 5.314 1.605) 6.342 1.926 ( 3.218 5.573 -1.292) 6.564 3.836 ( 6.283 10.883 8.525) 15.186 3.981 ( 4.665 8.080 -4.946) 10.560 7.805 ( 0.429 0.743 4.071) 4.160 7.886 ( 0.475 0.823 -3.194) 3.332 8.223 ( 3.819 6.614 0.513) 7.654 8.236 ( 4.339 7.516 -0.486) 8.692 8.743 ( -0.897 -1.554 0.822) 1.973 8.778 ( -0.063 -0.109 -2.433) 2.436 ======================= Grid point 681 (93/96) ======================= q-point: ( 0.29 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.275 ( -0.398 -1.872 1.294) 2.310 1.419 ( -2.140 -3.885 -11.111) 11.964 1.955 ( -1.776 5.764 2.208) 6.423 1.993 ( -1.560 5.544 -1.208) 5.885 4.029 ( 5.818 5.998 6.342) 10.491 4.121 ( 4.333 3.853 -2.294) 6.236 7.820 ( 0.657 0.285 2.422) 2.526 7.889 ( 0.223 -0.937 -3.882) 4.000 8.318 ( -1.095 8.184 2.818) 8.724 8.366 ( 2.107 6.026 -1.511) 6.560 8.716 ( -0.067 -1.728 2.542) 3.074 8.788 ( 1.312 -0.209 -3.997) 4.212 ======================= Grid point 682 (94/96) ======================= q-point: ( 0.36 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.264 ( 0.613 -0.890 -2.517) 2.739 1.356 ( -0.594 -3.153 -10.329) 10.816 1.964 ( -3.634 6.235 2.383) 7.600 2.002 ( -3.453 5.915 -1.147) 6.944 4.144 ( 2.558 1.824 4.707) 5.660 4.198 ( 2.116 0.167 -0.403) 2.161 7.832 ( 0.458 -0.010 1.041) 1.137 7.882 ( 0.588 -1.541 -3.568) 3.931 8.355 ( -3.440 7.490 3.839) 9.092 8.426 ( -1.531 5.373 -2.797) 6.247 8.707 ( 0.756 -1.488 4.277) 4.591 8.808 ( 0.730 0.072 -4.939) 4.993 ======================= Grid point 695 (95/96) ======================= q-point: ( 0.29 0.29 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.254 ( -0.043 -0.075 -2.842) 2.844 1.345 ( -1.919 -3.323 -10.518) 11.196 2.028 ( 1.071 1.856 2.230) 3.093 2.063 ( 1.030 1.784 -1.000) 2.290 4.129 ( 2.233 3.868 4.512) 6.349 4.179 ( 1.057 1.831 -0.164) 2.121 7.828 ( 0.309 0.535 0.483) 0.785 7.869 ( -0.550 -0.952 -3.237) 3.419 8.449 ( 1.787 3.095 0.498) 3.608 8.449 ( 2.204 3.818 0.809) 4.482 8.690 ( -0.479 -0.830 4.542) 4.642 8.801 ( 0.649 1.125 -5.414) 5.567 ======================= Grid point 696 (96/96) ======================= q-point: ( 0.36 0.29 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 2.72e-04 2.72e-04 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.263 ( -0.350 0.606 -5.714) 5.756 1.302 ( 1.163 -2.014 -9.054) 9.348 2.044 ( -1.098 1.902 2.243) 3.139 2.078 ( -1.046 1.811 -0.938) 2.292 4.173 ( -0.610 1.057 3.234) 3.456 4.194 ( 0.265 -0.459 1.214) 1.324 7.837 ( -0.254 0.440 -0.697) 0.863 7.856 ( 0.561 -0.972 -2.485) 2.727 8.470 ( -2.114 3.662 1.328) 4.432 8.494 ( -0.758 1.313 -1.002) 1.818 8.682 ( 0.496 -0.859 6.237) 6.315 8.815 ( -0.213 0.368 -5.963) 5.978 =================== 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/16464 10.0 451.887 451.887 491.089 -0.000 -0.000 0.000 3/16464 20.0 217.649 217.649 267.531 -0.000 -0.000 0.000 3/16464 30.0 168.834 168.834 191.151 -0.000 -0.000 0.000 3/16464 40.0 139.062 139.062 152.453 -0.000 -0.000 0.000 3/16464 50.0 117.170 117.170 126.355 -0.000 -0.000 0.000 3/16464 60.0 99.670 99.670 106.275 -0.000 -0.000 0.000 3/16464 70.0 85.389 85.389 90.181 -0.000 -0.000 0.000 3/16464 80.0 73.826 73.826 77.291 -0.000 -0.000 0.000 3/16464 90.0 64.550 64.550 67.038 -0.000 -0.000 0.000 3/16464 100.0 57.118 57.118 58.887 -0.000 -0.000 0.000 3/16464 110.0 51.127 51.127 52.367 -0.000 -0.000 0.000 3/16464 120.0 46.247 46.247 47.093 -0.000 -0.000 0.000 3/16464 130.0 42.221 42.221 42.772 -0.000 -0.000 0.000 3/16464 140.0 38.857 38.857 39.183 -0.000 -0.000 0.000 3/16464 150.0 36.010 36.010 36.164 -0.000 -0.000 0.000 3/16464 160.0 33.572 33.572 33.592 -0.000 -0.000 0.000 3/16464 170.0 31.462 31.462 31.377 -0.000 -0.000 0.000 3/16464 180.0 29.618 29.618 29.449 -0.000 -0.000 0.000 3/16464 190.0 27.992 27.992 27.757 -0.000 -0.000 0.000 3/16464 200.0 26.546 26.546 26.259 -0.000 -0.000 0.000 3/16464 210.0 25.253 25.253 24.923 -0.000 -0.000 0.000 3/16464 220.0 24.087 24.087 23.724 -0.000 -0.000 0.000 3/16464 230.0 23.032 23.032 22.642 -0.000 -0.000 0.000 3/16464 240.0 22.071 22.071 21.659 -0.000 -0.000 0.000 3/16464 250.0 21.192 21.192 20.762 -0.000 -0.000 0.000 3/16464 260.0 20.384 20.384 19.941 -0.000 -0.000 0.000 3/16464 270.0 19.639 19.639 19.185 -0.000 -0.000 0.000 3/16464 280.0 18.950 18.950 18.487 -0.000 -0.000 0.000 3/16464 290.0 18.310 18.310 17.841 -0.000 -0.000 0.000 3/16464 300.0 17.714 17.714 17.240 -0.000 -0.000 0.000 3/16464 310.0 17.157 17.157 16.680 -0.000 -0.000 0.000 3/16464 320.0 16.636 16.636 16.157 -0.000 -0.000 0.000 3/16464 330.0 16.147 16.147 15.667 -0.000 -0.000 0.000 3/16464 340.0 15.688 15.688 15.207 -0.000 -0.000 0.000 3/16464 350.0 15.254 15.254 14.775 -0.000 -0.000 0.000 3/16464 360.0 14.846 14.846 14.367 -0.000 -0.000 0.000 3/16464 370.0 14.459 14.459 13.982 -0.000 -0.000 0.000 3/16464 380.0 14.092 14.092 13.617 -0.000 -0.000 0.000 3/16464 390.0 13.745 13.745 13.272 -0.000 -0.000 0.000 3/16464 400.0 13.415 13.415 12.945 -0.000 -0.000 0.000 3/16464 410.0 13.100 13.100 12.633 -0.000 -0.000 0.000 3/16464 420.0 12.801 12.801 12.337 -0.000 -0.000 0.000 3/16464 430.0 12.515 12.515 12.055 -0.000 -0.000 0.000 3/16464 440.0 12.243 12.243 11.785 -0.000 -0.000 0.000 3/16464 450.0 11.982 11.982 11.528 -0.000 -0.000 0.000 3/16464 460.0 11.732 11.732 11.282 -0.000 -0.000 0.000 3/16464 470.0 11.493 11.493 11.047 -0.000 -0.000 0.000 3/16464 480.0 11.264 11.264 10.821 -0.000 -0.000 0.000 3/16464 490.0 11.043 11.043 10.605 -0.000 -0.000 0.000 3/16464 500.0 10.832 10.832 10.397 -0.000 -0.000 0.000 3/16464 510.0 10.629 10.629 10.197 -0.000 -0.000 0.000 3/16464 520.0 10.433 10.433 10.006 -0.000 -0.000 0.000 3/16464 530.0 10.244 10.244 9.821 -0.000 -0.000 0.000 3/16464 540.0 10.063 10.063 9.643 -0.000 -0.000 0.000 3/16464 550.0 9.887 9.887 9.472 -0.000 -0.000 0.000 3/16464 560.0 9.718 9.718 9.306 -0.000 -0.000 0.000 3/16464 570.0 9.555 9.555 9.147 -0.000 -0.000 0.000 3/16464 580.0 9.397 9.397 8.993 -0.000 -0.000 0.000 3/16464 590.0 9.245 9.245 8.844 -0.000 -0.000 0.000 3/16464 600.0 9.097 9.097 8.700 -0.000 -0.000 0.000 3/16464 610.0 8.954 8.954 8.561 -0.000 -0.000 0.000 3/16464 620.0 8.816 8.816 8.426 -0.000 -0.000 0.000 3/16464 630.0 8.682 8.682 8.295 -0.000 -0.000 0.000 3/16464 640.0 8.552 8.552 8.169 -0.000 -0.000 0.000 3/16464 650.0 8.425 8.425 8.046 -0.000 -0.000 0.000 3/16464 660.0 8.303 8.303 7.927 -0.000 -0.000 0.000 3/16464 670.0 8.184 8.184 7.812 -0.000 -0.000 0.000 3/16464 680.0 8.069 8.069 7.700 -0.000 -0.000 0.000 3/16464 690.0 7.956 7.956 7.591 -0.000 -0.000 0.000 3/16464 700.0 7.847 7.847 7.485 -0.000 -0.000 0.000 3/16464 710.0 7.741 7.741 7.382 -0.000 -0.000 0.000 3/16464 720.0 7.638 7.638 7.282 -0.000 -0.000 0.000 3/16464 730.0 7.538 7.538 7.185 -0.000 -0.000 0.000 3/16464 740.0 7.440 7.440 7.090 -0.000 -0.000 0.000 3/16464 750.0 7.344 7.344 6.998 -0.000 -0.000 0.000 3/16464 760.0 7.252 7.252 6.908 -0.000 -0.000 0.000 3/16464 770.0 7.161 7.161 6.820 -0.000 -0.000 0.000 3/16464 780.0 7.073 7.073 6.735 -0.000 -0.000 0.000 3/16464 790.0 6.987 6.987 6.652 -0.000 -0.000 0.000 3/16464 800.0 6.903 6.903 6.571 -0.000 -0.000 0.000 3/16464 810.0 6.821 6.821 6.491 -0.000 -0.000 0.000 3/16464 820.0 6.741 6.741 6.414 -0.000 -0.000 0.000 3/16464 830.0 6.662 6.662 6.339 -0.000 -0.000 0.000 3/16464 840.0 6.586 6.586 6.265 -0.000 -0.000 0.000 3/16464 850.0 6.511 6.511 6.193 -0.000 -0.000 0.000 3/16464 860.0 6.438 6.438 6.123 -0.000 -0.000 0.000 3/16464 870.0 6.367 6.367 6.054 -0.000 -0.000 0.000 3/16464 880.0 6.297 6.297 5.987 -0.000 -0.000 0.000 3/16464 890.0 6.229 6.229 5.921 -0.000 -0.000 0.000 3/16464 900.0 6.163 6.163 5.857 -0.000 -0.000 0.000 3/16464 910.0 6.097 6.097 5.794 -0.000 -0.000 0.000 3/16464 920.0 6.033 6.033 5.732 -0.000 -0.000 0.000 3/16464 930.0 5.971 5.971 5.672 -0.000 -0.000 0.000 3/16464 940.0 5.909 5.909 5.613 -0.000 -0.000 0.000 3/16464 950.0 5.849 5.849 5.555 -0.000 -0.000 0.000 3/16464 960.0 5.791 5.791 5.499 -0.000 -0.000 0.000 3/16464 970.0 5.733 5.733 5.443 -0.000 -0.000 0.000 3/16464 980.0 5.676 5.676 5.389 -0.000 -0.000 0.000 3/16464 990.0 5.621 5.621 5.336 -0.000 -0.000 0.000 3/16464 1000.0 5.567 5.567 5.284 -0.000 -0.000 0.000 3/16464 Thermal conductivity related properties were written into "kappa-m14147.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 18:36:23]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|