----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 23:44: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 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-6m2 (187) Number of symmetry operations in supercell: 576 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.488603559793336 0.000000000000000 0.000000000000000 b -2.244301779896668 3.887244710298293 0.000000000000000 c 0.000000000000000 0.000000000000000 4.539274940000000 Atomic positions (fractional): *1 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 *2 As 0.33333333333333 0.66666666666667 0.50000000000000 74.922 *3 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.488603559793336 0.000000000000000 0.000000000000000 b -2.244301779896668 3.887244710298293 0.000000000000000 c 0.000000000000000 0.000000000000000 4.539274940000000 Atomic positions (fractional): *1 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 > 1 *2 As 0.33333333333333 0.66666666666667 0.50000000000000 74.922 > 2 *3 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 > 3 -------------------------------- super cell -------------------------------- Lattice vectors: a 17.954414239173342 0.000000000000000 0.000000000000000 b -8.977207119586671 15.548978841193172 0.000000000000000 c 0.000000000000000 0.000000000000000 13.617824820000001 Atomic positions (fractional): *1 Li 0.00000000000000 0.00000000000000 0.16666666666667 6.941 > 1 2 Li 0.25000000000000 0.00000000000000 0.16666666666667 6.941 > 1 3 Li 0.50000000000000 0.00000000000000 0.16666666666667 6.941 > 1 4 Li 0.75000000000000 0.00000000000000 0.16666666666667 6.941 > 1 5 Li 0.00000000000000 0.25000000000000 0.16666666666667 6.941 > 1 6 Li 0.25000000000000 0.25000000000000 0.16666666666667 6.941 > 1 7 Li 0.50000000000000 0.25000000000000 0.16666666666667 6.941 > 1 8 Li 0.75000000000000 0.25000000000000 0.16666666666667 6.941 > 1 9 Li 0.00000000000000 0.50000000000000 0.16666666666667 6.941 > 1 10 Li 0.25000000000000 0.50000000000000 0.16666666666667 6.941 > 1 11 Li 0.50000000000000 0.50000000000000 0.16666666666667 6.941 > 1 12 Li 0.75000000000000 0.50000000000000 0.16666666666667 6.941 > 1 13 Li 0.00000000000000 0.75000000000000 0.16666666666667 6.941 > 1 14 Li 0.25000000000000 0.75000000000000 0.16666666666667 6.941 > 1 15 Li 0.50000000000000 0.75000000000000 0.16666666666667 6.941 > 1 16 Li 0.75000000000000 0.75000000000000 0.16666666666667 6.941 > 1 17 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 > 1 18 Li 0.25000000000000 0.00000000000000 0.50000000000000 6.941 > 1 19 Li 0.50000000000000 0.00000000000000 0.50000000000000 6.941 > 1 20 Li 0.75000000000000 0.00000000000000 0.50000000000000 6.941 > 1 21 Li 0.00000000000000 0.25000000000000 0.50000000000000 6.941 > 1 22 Li 0.25000000000000 0.25000000000000 0.50000000000000 6.941 > 1 23 Li 0.50000000000000 0.25000000000000 0.50000000000000 6.941 > 1 24 Li 0.75000000000000 0.25000000000000 0.50000000000000 6.941 > 1 25 Li 0.00000000000000 0.50000000000000 0.50000000000000 6.941 > 1 26 Li 0.25000000000000 0.50000000000000 0.50000000000000 6.941 > 1 27 Li 0.50000000000000 0.50000000000000 0.50000000000000 6.941 > 1 28 Li 0.75000000000000 0.50000000000000 0.50000000000000 6.941 > 1 29 Li 0.00000000000000 0.75000000000000 0.50000000000000 6.941 > 1 30 Li 0.25000000000000 0.75000000000000 0.50000000000000 6.941 > 1 31 Li 0.50000000000000 0.75000000000000 0.50000000000000 6.941 > 1 32 Li 0.75000000000000 0.75000000000000 0.50000000000000 6.941 > 1 33 Li 0.00000000000000 0.00000000000000 0.83333333333333 6.941 > 1 34 Li 0.25000000000000 0.00000000000000 0.83333333333333 6.941 > 1 35 Li 0.50000000000000 0.00000000000000 0.83333333333333 6.941 > 1 36 Li 0.75000000000000 0.00000000000000 0.83333333333333 6.941 > 1 37 Li 0.00000000000000 0.25000000000000 0.83333333333333 6.941 > 1 38 Li 0.25000000000000 0.25000000000000 0.83333333333333 6.941 > 1 39 Li 0.50000000000000 0.25000000000000 0.83333333333333 6.941 > 1 40 Li 0.75000000000000 0.25000000000000 0.83333333333333 6.941 > 1 41 Li 0.00000000000000 0.50000000000000 0.83333333333333 6.941 > 1 42 Li 0.25000000000000 0.50000000000000 0.83333333333333 6.941 > 1 43 Li 0.50000000000000 0.50000000000000 0.83333333333333 6.941 > 1 44 Li 0.75000000000000 0.50000000000000 0.83333333333333 6.941 > 1 45 Li 0.00000000000000 0.75000000000000 0.83333333333333 6.941 > 1 46 Li 0.25000000000000 0.75000000000000 0.83333333333333 6.941 > 1 47 Li 0.50000000000000 0.75000000000000 0.83333333333333 6.941 > 1 48 Li 0.75000000000000 0.75000000000000 0.83333333333333 6.941 > 1 *49 As 0.08333333333333 0.16666666666667 0.16666666666667 74.922 > 2 50 As 0.33333333333333 0.16666666666667 0.16666666666667 74.922 > 2 51 As 0.58333333333333 0.16666666666667 0.16666666666667 74.922 > 2 52 As 0.83333333333333 0.16666666666667 0.16666666666667 74.922 > 2 53 As 0.08333333333333 0.41666666666667 0.16666666666667 74.922 > 2 54 As 0.33333333333333 0.41666666666667 0.16666666666667 74.922 > 2 55 As 0.58333333333333 0.41666666666667 0.16666666666667 74.922 > 2 56 As 0.83333333333333 0.41666666666667 0.16666666666667 74.922 > 2 57 As 0.08333333333333 0.66666666666667 0.16666666666667 74.922 > 2 58 As 0.33333333333333 0.66666666666667 0.16666666666667 74.922 > 2 59 As 0.58333333333333 0.66666666666667 0.16666666666667 74.922 > 2 60 As 0.83333333333333 0.66666666666667 0.16666666666667 74.922 > 2 61 As 0.08333333333333 0.91666666666667 0.16666666666667 74.922 > 2 62 As 0.33333333333333 0.91666666666667 0.16666666666667 74.922 > 2 63 As 0.58333333333333 0.91666666666667 0.16666666666667 74.922 > 2 64 As 0.83333333333333 0.91666666666667 0.16666666666667 74.922 > 2 65 As 0.08333333333333 0.16666666666667 0.50000000000000 74.922 > 2 66 As 0.33333333333333 0.16666666666667 0.50000000000000 74.922 > 2 67 As 0.58333333333333 0.16666666666667 0.50000000000000 74.922 > 2 68 As 0.83333333333333 0.16666666666667 0.50000000000000 74.922 > 2 69 As 0.08333333333333 0.41666666666667 0.50000000000000 74.922 > 2 70 As 0.33333333333333 0.41666666666667 0.50000000000000 74.922 > 2 71 As 0.58333333333333 0.41666666666667 0.50000000000000 74.922 > 2 72 As 0.83333333333333 0.41666666666667 0.50000000000000 74.922 > 2 73 As 0.08333333333333 0.66666666666667 0.50000000000000 74.922 > 2 74 As 0.33333333333333 0.66666666666667 0.50000000000000 74.922 > 2 75 As 0.58333333333333 0.66666666666667 0.50000000000000 74.922 > 2 76 As 0.83333333333333 0.66666666666667 0.50000000000000 74.922 > 2 77 As 0.08333333333333 0.91666666666667 0.50000000000000 74.922 > 2 78 As 0.33333333333333 0.91666666666667 0.50000000000000 74.922 > 2 79 As 0.58333333333333 0.91666666666667 0.50000000000000 74.922 > 2 80 As 0.83333333333333 0.91666666666667 0.50000000000000 74.922 > 2 81 As 0.08333333333333 0.16666666666667 0.83333333333333 74.922 > 2 82 As 0.33333333333333 0.16666666666667 0.83333333333333 74.922 > 2 83 As 0.58333333333333 0.16666666666667 0.83333333333333 74.922 > 2 84 As 0.83333333333333 0.16666666666667 0.83333333333333 74.922 > 2 85 As 0.08333333333333 0.41666666666667 0.83333333333333 74.922 > 2 86 As 0.33333333333333 0.41666666666667 0.83333333333333 74.922 > 2 87 As 0.58333333333333 0.41666666666667 0.83333333333333 74.922 > 2 88 As 0.83333333333333 0.41666666666667 0.83333333333333 74.922 > 2 89 As 0.08333333333333 0.66666666666667 0.83333333333333 74.922 > 2 90 As 0.33333333333333 0.66666666666667 0.83333333333333 74.922 > 2 91 As 0.58333333333333 0.66666666666667 0.83333333333333 74.922 > 2 92 As 0.83333333333333 0.66666666666667 0.83333333333333 74.922 > 2 93 As 0.08333333333333 0.91666666666667 0.83333333333333 74.922 > 2 94 As 0.33333333333333 0.91666666666667 0.83333333333333 74.922 > 2 95 As 0.58333333333333 0.91666666666667 0.83333333333333 74.922 > 2 96 As 0.83333333333333 0.91666666666667 0.83333333333333 74.922 > 2 *97 Ba 0.16666666666667 0.08333333333333 0.00000000000000 137.327 > 3 98 Ba 0.41666666666667 0.08333333333333 0.00000000000000 137.327 > 3 99 Ba 0.66666666666667 0.08333333333333 0.00000000000000 137.327 > 3 100 Ba 0.91666666666667 0.08333333333333 0.00000000000000 137.327 > 3 101 Ba 0.16666666666667 0.33333333333333 0.00000000000000 137.327 > 3 102 Ba 0.41666666666667 0.33333333333333 0.00000000000000 137.327 > 3 103 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 > 3 104 Ba 0.91666666666667 0.33333333333333 0.00000000000000 137.327 > 3 105 Ba 0.16666666666667 0.58333333333333 0.00000000000000 137.327 > 3 106 Ba 0.41666666666667 0.58333333333333 0.00000000000000 137.327 > 3 107 Ba 0.66666666666667 0.58333333333333 0.00000000000000 137.327 > 3 108 Ba 0.91666666666667 0.58333333333333 0.00000000000000 137.327 > 3 109 Ba 0.16666666666667 0.83333333333333 0.00000000000000 137.327 > 3 110 Ba 0.41666666666667 0.83333333333333 0.00000000000000 137.327 > 3 111 Ba 0.66666666666667 0.83333333333333 0.00000000000000 137.327 > 3 112 Ba 0.91666666666667 0.83333333333333 0.00000000000000 137.327 > 3 113 Ba 0.16666666666667 0.08333333333333 0.33333333333333 137.327 > 3 114 Ba 0.41666666666667 0.08333333333333 0.33333333333333 137.327 > 3 115 Ba 0.66666666666667 0.08333333333333 0.33333333333333 137.327 > 3 116 Ba 0.91666666666667 0.08333333333333 0.33333333333333 137.327 > 3 117 Ba 0.16666666666667 0.33333333333333 0.33333333333333 137.327 > 3 118 Ba 0.41666666666667 0.33333333333333 0.33333333333333 137.327 > 3 119 Ba 0.66666666666667 0.33333333333333 0.33333333333333 137.327 > 3 120 Ba 0.91666666666667 0.33333333333333 0.33333333333333 137.327 > 3 121 Ba 0.16666666666667 0.58333333333333 0.33333333333333 137.327 > 3 122 Ba 0.41666666666667 0.58333333333333 0.33333333333333 137.327 > 3 123 Ba 0.66666666666667 0.58333333333333 0.33333333333333 137.327 > 3 124 Ba 0.91666666666667 0.58333333333333 0.33333333333333 137.327 > 3 125 Ba 0.16666666666667 0.83333333333333 0.33333333333333 137.327 > 3 126 Ba 0.41666666666667 0.83333333333333 0.33333333333333 137.327 > 3 127 Ba 0.66666666666667 0.83333333333333 0.33333333333333 137.327 > 3 128 Ba 0.91666666666667 0.83333333333333 0.33333333333333 137.327 > 3 129 Ba 0.16666666666667 0.08333333333333 0.66666666666667 137.327 > 3 130 Ba 0.41666666666667 0.08333333333333 0.66666666666667 137.327 > 3 131 Ba 0.66666666666667 0.08333333333333 0.66666666666667 137.327 > 3 132 Ba 0.91666666666667 0.08333333333333 0.66666666666667 137.327 > 3 133 Ba 0.16666666666667 0.33333333333333 0.66666666666667 137.327 > 3 134 Ba 0.41666666666667 0.33333333333333 0.66666666666667 137.327 > 3 135 Ba 0.66666666666667 0.33333333333333 0.66666666666667 137.327 > 3 136 Ba 0.91666666666667 0.33333333333333 0.66666666666667 137.327 > 3 137 Ba 0.16666666666667 0.58333333333333 0.66666666666667 137.327 > 3 138 Ba 0.41666666666667 0.58333333333333 0.66666666666667 137.327 > 3 139 Ba 0.66666666666667 0.58333333333333 0.66666666666667 137.327 > 3 140 Ba 0.91666666666667 0.58333333333333 0.66666666666667 137.327 > 3 141 Ba 0.16666666666667 0.83333333333333 0.66666666666667 137.327 > 3 142 Ba 0.41666666666667 0.83333333333333 0.66666666666667 137.327 > 3 143 Ba 0.66666666666667 0.83333333333333 0.66666666666667 137.327 > 3 144 Ba 0.91666666666667 0.83333333333333 0.66666666666667 137.327 > 3 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.6531051 -0.0000000 0.0000000 0.0000000 9.6531051 0.0000000 0.0000000 0.0000000 45.3917432 -------------------------- Born effective charges -------------------------- 1 Li 0.9957698 0.0000000 0.0000000 0.0000000 0.9957698 0.0000000 0.0000000 0.0000000 -0.0040015 2 As -3.1893956 -0.0000000 0.0000000 -0.0000000 -3.1893956 0.0000000 0.0000000 0.0000000 -5.8823552 3 Ba 2.1936258 0.0000000 0.0000000 0.0000000 2.1936258 0.0000000 0.0000000 0.0000000 5.8863567 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 432/432 Permutation basis: 2574/2574 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 213 Number of blocks in projector: 147 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 114 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 54 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 45 Use standard eigh solver. Tree of FC basis block matrices: - (213, 207), data: False |-- (45, 45), data: True |-- (54, 51), data: True |-- (114, 111), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.009 Solver_block: 80 / 80 - Time: 0.192 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.223 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) Permutation basis: 432/432 Permutation basis: 2574/2574 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 213 Number of blocks in projector: 147 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 114 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 54 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 45 Use standard eigh solver. Tree of FC basis block matrices: - (213, 207), data: False |-- (45, 45), data: True |-- (54, 51), data: True |-- (114, 111), data: True Max drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) 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 23:44:39]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 23:44: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 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-6m2 (187) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.488603559793336 0.000000000000000 0.000000000000000 b -2.244301779896668 3.887244710298293 0.000000000000000 c 0.000000000000000 0.000000000000000 4.539274940000000 Atomic positions (fractional): 1 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 2 As 0.33333333333333 0.66666666666667 0.50000000000000 74.922 3 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.954414239173342 0.000000000000000 0.000000000000000 b -8.977207119586671 15.548978841193172 0.000000000000000 c 0.000000000000000 0.000000000000000 13.617824820000001 Atomic positions (fractional): 1 Li 0.00000000000000 0.00000000000000 0.16666666666667 6.941 > 1 2 Li 0.25000000000000 0.00000000000000 0.16666666666667 6.941 > 1 3 Li 0.50000000000000 0.00000000000000 0.16666666666667 6.941 > 1 4 Li 0.75000000000000 0.00000000000000 0.16666666666667 6.941 > 1 5 Li 0.00000000000000 0.25000000000000 0.16666666666667 6.941 > 1 6 Li 0.25000000000000 0.25000000000000 0.16666666666667 6.941 > 1 7 Li 0.50000000000000 0.25000000000000 0.16666666666667 6.941 > 1 8 Li 0.75000000000000 0.25000000000000 0.16666666666667 6.941 > 1 9 Li 0.00000000000000 0.50000000000000 0.16666666666667 6.941 > 1 10 Li 0.25000000000000 0.50000000000000 0.16666666666667 6.941 > 1 11 Li 0.50000000000000 0.50000000000000 0.16666666666667 6.941 > 1 12 Li 0.75000000000000 0.50000000000000 0.16666666666667 6.941 > 1 13 Li 0.00000000000000 0.75000000000000 0.16666666666667 6.941 > 1 14 Li 0.25000000000000 0.75000000000000 0.16666666666667 6.941 > 1 15 Li 0.50000000000000 0.75000000000000 0.16666666666667 6.941 > 1 16 Li 0.75000000000000 0.75000000000000 0.16666666666667 6.941 > 1 17 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 > 1 18 Li 0.25000000000000 0.00000000000000 0.50000000000000 6.941 > 1 19 Li 0.50000000000000 0.00000000000000 0.50000000000000 6.941 > 1 20 Li 0.75000000000000 0.00000000000000 0.50000000000000 6.941 > 1 21 Li 0.00000000000000 0.25000000000000 0.50000000000000 6.941 > 1 22 Li 0.25000000000000 0.25000000000000 0.50000000000000 6.941 > 1 23 Li 0.50000000000000 0.25000000000000 0.50000000000000 6.941 > 1 24 Li 0.75000000000000 0.25000000000000 0.50000000000000 6.941 > 1 25 Li 0.00000000000000 0.50000000000000 0.50000000000000 6.941 > 1 26 Li 0.25000000000000 0.50000000000000 0.50000000000000 6.941 > 1 27 Li 0.50000000000000 0.50000000000000 0.50000000000000 6.941 > 1 28 Li 0.75000000000000 0.50000000000000 0.50000000000000 6.941 > 1 29 Li 0.00000000000000 0.75000000000000 0.50000000000000 6.941 > 1 30 Li 0.25000000000000 0.75000000000000 0.50000000000000 6.941 > 1 31 Li 0.50000000000000 0.75000000000000 0.50000000000000 6.941 > 1 32 Li 0.75000000000000 0.75000000000000 0.50000000000000 6.941 > 1 33 Li 0.00000000000000 0.00000000000000 0.83333333333333 6.941 > 1 34 Li 0.25000000000000 0.00000000000000 0.83333333333333 6.941 > 1 35 Li 0.50000000000000 0.00000000000000 0.83333333333333 6.941 > 1 36 Li 0.75000000000000 0.00000000000000 0.83333333333333 6.941 > 1 37 Li 0.00000000000000 0.25000000000000 0.83333333333333 6.941 > 1 38 Li 0.25000000000000 0.25000000000000 0.83333333333333 6.941 > 1 39 Li 0.50000000000000 0.25000000000000 0.83333333333333 6.941 > 1 40 Li 0.75000000000000 0.25000000000000 0.83333333333333 6.941 > 1 41 Li 0.00000000000000 0.50000000000000 0.83333333333333 6.941 > 1 42 Li 0.25000000000000 0.50000000000000 0.83333333333333 6.941 > 1 43 Li 0.50000000000000 0.50000000000000 0.83333333333333 6.941 > 1 44 Li 0.75000000000000 0.50000000000000 0.83333333333333 6.941 > 1 45 Li 0.00000000000000 0.75000000000000 0.83333333333333 6.941 > 1 46 Li 0.25000000000000 0.75000000000000 0.83333333333333 6.941 > 1 47 Li 0.50000000000000 0.75000000000000 0.83333333333333 6.941 > 1 48 Li 0.75000000000000 0.75000000000000 0.83333333333333 6.941 > 1 49 As 0.08333333333333 0.16666666666667 0.16666666666667 74.922 > 49 50 As 0.33333333333333 0.16666666666667 0.16666666666667 74.922 > 49 51 As 0.58333333333333 0.16666666666667 0.16666666666667 74.922 > 49 52 As 0.83333333333333 0.16666666666667 0.16666666666667 74.922 > 49 53 As 0.08333333333333 0.41666666666667 0.16666666666667 74.922 > 49 54 As 0.33333333333333 0.41666666666667 0.16666666666667 74.922 > 49 55 As 0.58333333333333 0.41666666666667 0.16666666666667 74.922 > 49 56 As 0.83333333333333 0.41666666666667 0.16666666666667 74.922 > 49 57 As 0.08333333333333 0.66666666666667 0.16666666666667 74.922 > 49 58 As 0.33333333333333 0.66666666666667 0.16666666666667 74.922 > 49 59 As 0.58333333333333 0.66666666666667 0.16666666666667 74.922 > 49 60 As 0.83333333333333 0.66666666666667 0.16666666666667 74.922 > 49 61 As 0.08333333333333 0.91666666666667 0.16666666666667 74.922 > 49 62 As 0.33333333333333 0.91666666666667 0.16666666666667 74.922 > 49 63 As 0.58333333333333 0.91666666666667 0.16666666666667 74.922 > 49 64 As 0.83333333333333 0.91666666666667 0.16666666666667 74.922 > 49 65 As 0.08333333333333 0.16666666666667 0.50000000000000 74.922 > 49 66 As 0.33333333333333 0.16666666666667 0.50000000000000 74.922 > 49 67 As 0.58333333333333 0.16666666666667 0.50000000000000 74.922 > 49 68 As 0.83333333333333 0.16666666666667 0.50000000000000 74.922 > 49 69 As 0.08333333333333 0.41666666666667 0.50000000000000 74.922 > 49 70 As 0.33333333333333 0.41666666666667 0.50000000000000 74.922 > 49 71 As 0.58333333333333 0.41666666666667 0.50000000000000 74.922 > 49 72 As 0.83333333333333 0.41666666666667 0.50000000000000 74.922 > 49 73 As 0.08333333333333 0.66666666666667 0.50000000000000 74.922 > 49 74 As 0.33333333333333 0.66666666666667 0.50000000000000 74.922 > 49 75 As 0.58333333333333 0.66666666666667 0.50000000000000 74.922 > 49 76 As 0.83333333333333 0.66666666666667 0.50000000000000 74.922 > 49 77 As 0.08333333333333 0.91666666666667 0.50000000000000 74.922 > 49 78 As 0.33333333333333 0.91666666666667 0.50000000000000 74.922 > 49 79 As 0.58333333333333 0.91666666666667 0.50000000000000 74.922 > 49 80 As 0.83333333333333 0.91666666666667 0.50000000000000 74.922 > 49 81 As 0.08333333333333 0.16666666666667 0.83333333333333 74.922 > 49 82 As 0.33333333333333 0.16666666666667 0.83333333333333 74.922 > 49 83 As 0.58333333333333 0.16666666666667 0.83333333333333 74.922 > 49 84 As 0.83333333333333 0.16666666666667 0.83333333333333 74.922 > 49 85 As 0.08333333333333 0.41666666666667 0.83333333333333 74.922 > 49 86 As 0.33333333333333 0.41666666666667 0.83333333333333 74.922 > 49 87 As 0.58333333333333 0.41666666666667 0.83333333333333 74.922 > 49 88 As 0.83333333333333 0.41666666666667 0.83333333333333 74.922 > 49 89 As 0.08333333333333 0.66666666666667 0.83333333333333 74.922 > 49 90 As 0.33333333333333 0.66666666666667 0.83333333333333 74.922 > 49 91 As 0.58333333333333 0.66666666666667 0.83333333333333 74.922 > 49 92 As 0.83333333333333 0.66666666666667 0.83333333333333 74.922 > 49 93 As 0.08333333333333 0.91666666666667 0.83333333333333 74.922 > 49 94 As 0.33333333333333 0.91666666666667 0.83333333333333 74.922 > 49 95 As 0.58333333333333 0.91666666666667 0.83333333333333 74.922 > 49 96 As 0.83333333333333 0.91666666666667 0.83333333333333 74.922 > 49 97 Ba 0.16666666666667 0.08333333333333 0.00000000000000 137.327 > 97 98 Ba 0.41666666666667 0.08333333333333 0.00000000000000 137.327 > 97 99 Ba 0.66666666666667 0.08333333333333 0.00000000000000 137.327 > 97 100 Ba 0.91666666666667 0.08333333333333 0.00000000000000 137.327 > 97 101 Ba 0.16666666666667 0.33333333333333 0.00000000000000 137.327 > 97 102 Ba 0.41666666666667 0.33333333333333 0.00000000000000 137.327 > 97 103 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 > 97 104 Ba 0.91666666666667 0.33333333333333 0.00000000000000 137.327 > 97 105 Ba 0.16666666666667 0.58333333333333 0.00000000000000 137.327 > 97 106 Ba 0.41666666666667 0.58333333333333 0.00000000000000 137.327 > 97 107 Ba 0.66666666666667 0.58333333333333 0.00000000000000 137.327 > 97 108 Ba 0.91666666666667 0.58333333333333 0.00000000000000 137.327 > 97 109 Ba 0.16666666666667 0.83333333333333 0.00000000000000 137.327 > 97 110 Ba 0.41666666666667 0.83333333333333 0.00000000000000 137.327 > 97 111 Ba 0.66666666666667 0.83333333333333 0.00000000000000 137.327 > 97 112 Ba 0.91666666666667 0.83333333333333 0.00000000000000 137.327 > 97 113 Ba 0.16666666666667 0.08333333333333 0.33333333333333 137.327 > 97 114 Ba 0.41666666666667 0.08333333333333 0.33333333333333 137.327 > 97 115 Ba 0.66666666666667 0.08333333333333 0.33333333333333 137.327 > 97 116 Ba 0.91666666666667 0.08333333333333 0.33333333333333 137.327 > 97 117 Ba 0.16666666666667 0.33333333333333 0.33333333333333 137.327 > 97 118 Ba 0.41666666666667 0.33333333333333 0.33333333333333 137.327 > 97 119 Ba 0.66666666666667 0.33333333333333 0.33333333333333 137.327 > 97 120 Ba 0.91666666666667 0.33333333333333 0.33333333333333 137.327 > 97 121 Ba 0.16666666666667 0.58333333333333 0.33333333333333 137.327 > 97 122 Ba 0.41666666666667 0.58333333333333 0.33333333333333 137.327 > 97 123 Ba 0.66666666666667 0.58333333333333 0.33333333333333 137.327 > 97 124 Ba 0.91666666666667 0.58333333333333 0.33333333333333 137.327 > 97 125 Ba 0.16666666666667 0.83333333333333 0.33333333333333 137.327 > 97 126 Ba 0.41666666666667 0.83333333333333 0.33333333333333 137.327 > 97 127 Ba 0.66666666666667 0.83333333333333 0.33333333333333 137.327 > 97 128 Ba 0.91666666666667 0.83333333333333 0.33333333333333 137.327 > 97 129 Ba 0.16666666666667 0.08333333333333 0.66666666666667 137.327 > 97 130 Ba 0.41666666666667 0.08333333333333 0.66666666666667 137.327 > 97 131 Ba 0.66666666666667 0.08333333333333 0.66666666666667 137.327 > 97 132 Ba 0.91666666666667 0.08333333333333 0.66666666666667 137.327 > 97 133 Ba 0.16666666666667 0.33333333333333 0.66666666666667 137.327 > 97 134 Ba 0.41666666666667 0.33333333333333 0.66666666666667 137.327 > 97 135 Ba 0.66666666666667 0.33333333333333 0.66666666666667 137.327 > 97 136 Ba 0.91666666666667 0.33333333333333 0.66666666666667 137.327 > 97 137 Ba 0.16666666666667 0.58333333333333 0.66666666666667 137.327 > 97 138 Ba 0.41666666666667 0.58333333333333 0.66666666666667 137.327 > 97 139 Ba 0.66666666666667 0.58333333333333 0.66666666666667 137.327 > 97 140 Ba 0.91666666666667 0.58333333333333 0.66666666666667 137.327 > 97 141 Ba 0.16666666666667 0.83333333333333 0.66666666666667 137.327 > 97 142 Ba 0.41666666666667 0.83333333333333 0.66666666666667 137.327 > 97 143 Ba 0.66666666666667 0.83333333333333 0.66666666666667 137.327 > 97 144 Ba 0.91666666666667 0.83333333333333 0.66666666666667 137.327 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 9.6531051 -0.0000000 0.0000000 0.0000000 9.6531051 0.0000000 0.0000000 0.0000000 45.3917432 -------------------------- Born effective charges -------------------------- 1 Li 0.9957698 0.0000000 0.0000000 0.0000000 0.9957698 0.0000000 0.0000000 0.0000000 -0.0040015 2 As -3.1893956 -0.0000000 0.0000000 -0.0000000 -3.1893956 0.0000000 0.0000000 0.0000000 -5.8823552 3 Ba 2.1936258 0.0000000 0.0000000 0.0000000 2.1936258 0.0000000 0.0000000 0.0000000 5.8863567 ---------------------------------------------------------------------------- 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[ 49, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 97, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000012 (xxy) -0.00000012 (xxy) -0.00000012 (xyx) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) 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 23:44:45]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 23:44:46]------------------------- 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 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-6m2 (187) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.488603559793336 0.000000000000000 0.000000000000000 b -2.244301779896668 3.887244710298293 0.000000000000000 c 0.000000000000000 0.000000000000000 4.539274940000000 Atomic positions (fractional): 1 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 2 As 0.33333333333333 0.66666666666667 0.50000000000000 74.922 3 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.954414239173342 0.000000000000000 0.000000000000000 b -8.977207119586671 15.548978841193172 0.000000000000000 c 0.000000000000000 0.000000000000000 13.617824820000001 Atomic positions (fractional): 1 Li 0.00000000000000 0.00000000000000 0.16666666666667 6.941 > 1 2 Li 0.25000000000000 0.00000000000000 0.16666666666667 6.941 > 1 3 Li 0.50000000000000 0.00000000000000 0.16666666666667 6.941 > 1 4 Li 0.75000000000000 0.00000000000000 0.16666666666667 6.941 > 1 5 Li 0.00000000000000 0.25000000000000 0.16666666666667 6.941 > 1 6 Li 0.25000000000000 0.25000000000000 0.16666666666667 6.941 > 1 7 Li 0.50000000000000 0.25000000000000 0.16666666666667 6.941 > 1 8 Li 0.75000000000000 0.25000000000000 0.16666666666667 6.941 > 1 9 Li 0.00000000000000 0.50000000000000 0.16666666666667 6.941 > 1 10 Li 0.25000000000000 0.50000000000000 0.16666666666667 6.941 > 1 11 Li 0.50000000000000 0.50000000000000 0.16666666666667 6.941 > 1 12 Li 0.75000000000000 0.50000000000000 0.16666666666667 6.941 > 1 13 Li 0.00000000000000 0.75000000000000 0.16666666666667 6.941 > 1 14 Li 0.25000000000000 0.75000000000000 0.16666666666667 6.941 > 1 15 Li 0.50000000000000 0.75000000000000 0.16666666666667 6.941 > 1 16 Li 0.75000000000000 0.75000000000000 0.16666666666667 6.941 > 1 17 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 > 1 18 Li 0.25000000000000 0.00000000000000 0.50000000000000 6.941 > 1 19 Li 0.50000000000000 0.00000000000000 0.50000000000000 6.941 > 1 20 Li 0.75000000000000 0.00000000000000 0.50000000000000 6.941 > 1 21 Li 0.00000000000000 0.25000000000000 0.50000000000000 6.941 > 1 22 Li 0.25000000000000 0.25000000000000 0.50000000000000 6.941 > 1 23 Li 0.50000000000000 0.25000000000000 0.50000000000000 6.941 > 1 24 Li 0.75000000000000 0.25000000000000 0.50000000000000 6.941 > 1 25 Li 0.00000000000000 0.50000000000000 0.50000000000000 6.941 > 1 26 Li 0.25000000000000 0.50000000000000 0.50000000000000 6.941 > 1 27 Li 0.50000000000000 0.50000000000000 0.50000000000000 6.941 > 1 28 Li 0.75000000000000 0.50000000000000 0.50000000000000 6.941 > 1 29 Li 0.00000000000000 0.75000000000000 0.50000000000000 6.941 > 1 30 Li 0.25000000000000 0.75000000000000 0.50000000000000 6.941 > 1 31 Li 0.50000000000000 0.75000000000000 0.50000000000000 6.941 > 1 32 Li 0.75000000000000 0.75000000000000 0.50000000000000 6.941 > 1 33 Li 0.00000000000000 0.00000000000000 0.83333333333333 6.941 > 1 34 Li 0.25000000000000 0.00000000000000 0.83333333333333 6.941 > 1 35 Li 0.50000000000000 0.00000000000000 0.83333333333333 6.941 > 1 36 Li 0.75000000000000 0.00000000000000 0.83333333333333 6.941 > 1 37 Li 0.00000000000000 0.25000000000000 0.83333333333333 6.941 > 1 38 Li 0.25000000000000 0.25000000000000 0.83333333333333 6.941 > 1 39 Li 0.50000000000000 0.25000000000000 0.83333333333333 6.941 > 1 40 Li 0.75000000000000 0.25000000000000 0.83333333333333 6.941 > 1 41 Li 0.00000000000000 0.50000000000000 0.83333333333333 6.941 > 1 42 Li 0.25000000000000 0.50000000000000 0.83333333333333 6.941 > 1 43 Li 0.50000000000000 0.50000000000000 0.83333333333333 6.941 > 1 44 Li 0.75000000000000 0.50000000000000 0.83333333333333 6.941 > 1 45 Li 0.00000000000000 0.75000000000000 0.83333333333333 6.941 > 1 46 Li 0.25000000000000 0.75000000000000 0.83333333333333 6.941 > 1 47 Li 0.50000000000000 0.75000000000000 0.83333333333333 6.941 > 1 48 Li 0.75000000000000 0.75000000000000 0.83333333333333 6.941 > 1 49 As 0.08333333333333 0.16666666666667 0.16666666666667 74.922 > 49 50 As 0.33333333333333 0.16666666666667 0.16666666666667 74.922 > 49 51 As 0.58333333333333 0.16666666666667 0.16666666666667 74.922 > 49 52 As 0.83333333333333 0.16666666666667 0.16666666666667 74.922 > 49 53 As 0.08333333333333 0.41666666666667 0.16666666666667 74.922 > 49 54 As 0.33333333333333 0.41666666666667 0.16666666666667 74.922 > 49 55 As 0.58333333333333 0.41666666666667 0.16666666666667 74.922 > 49 56 As 0.83333333333333 0.41666666666667 0.16666666666667 74.922 > 49 57 As 0.08333333333333 0.66666666666667 0.16666666666667 74.922 > 49 58 As 0.33333333333333 0.66666666666667 0.16666666666667 74.922 > 49 59 As 0.58333333333333 0.66666666666667 0.16666666666667 74.922 > 49 60 As 0.83333333333333 0.66666666666667 0.16666666666667 74.922 > 49 61 As 0.08333333333333 0.91666666666667 0.16666666666667 74.922 > 49 62 As 0.33333333333333 0.91666666666667 0.16666666666667 74.922 > 49 63 As 0.58333333333333 0.91666666666667 0.16666666666667 74.922 > 49 64 As 0.83333333333333 0.91666666666667 0.16666666666667 74.922 > 49 65 As 0.08333333333333 0.16666666666667 0.50000000000000 74.922 > 49 66 As 0.33333333333333 0.16666666666667 0.50000000000000 74.922 > 49 67 As 0.58333333333333 0.16666666666667 0.50000000000000 74.922 > 49 68 As 0.83333333333333 0.16666666666667 0.50000000000000 74.922 > 49 69 As 0.08333333333333 0.41666666666667 0.50000000000000 74.922 > 49 70 As 0.33333333333333 0.41666666666667 0.50000000000000 74.922 > 49 71 As 0.58333333333333 0.41666666666667 0.50000000000000 74.922 > 49 72 As 0.83333333333333 0.41666666666667 0.50000000000000 74.922 > 49 73 As 0.08333333333333 0.66666666666667 0.50000000000000 74.922 > 49 74 As 0.33333333333333 0.66666666666667 0.50000000000000 74.922 > 49 75 As 0.58333333333333 0.66666666666667 0.50000000000000 74.922 > 49 76 As 0.83333333333333 0.66666666666667 0.50000000000000 74.922 > 49 77 As 0.08333333333333 0.91666666666667 0.50000000000000 74.922 > 49 78 As 0.33333333333333 0.91666666666667 0.50000000000000 74.922 > 49 79 As 0.58333333333333 0.91666666666667 0.50000000000000 74.922 > 49 80 As 0.83333333333333 0.91666666666667 0.50000000000000 74.922 > 49 81 As 0.08333333333333 0.16666666666667 0.83333333333333 74.922 > 49 82 As 0.33333333333333 0.16666666666667 0.83333333333333 74.922 > 49 83 As 0.58333333333333 0.16666666666667 0.83333333333333 74.922 > 49 84 As 0.83333333333333 0.16666666666667 0.83333333333333 74.922 > 49 85 As 0.08333333333333 0.41666666666667 0.83333333333333 74.922 > 49 86 As 0.33333333333333 0.41666666666667 0.83333333333333 74.922 > 49 87 As 0.58333333333333 0.41666666666667 0.83333333333333 74.922 > 49 88 As 0.83333333333333 0.41666666666667 0.83333333333333 74.922 > 49 89 As 0.08333333333333 0.66666666666667 0.83333333333333 74.922 > 49 90 As 0.33333333333333 0.66666666666667 0.83333333333333 74.922 > 49 91 As 0.58333333333333 0.66666666666667 0.83333333333333 74.922 > 49 92 As 0.83333333333333 0.66666666666667 0.83333333333333 74.922 > 49 93 As 0.08333333333333 0.91666666666667 0.83333333333333 74.922 > 49 94 As 0.33333333333333 0.91666666666667 0.83333333333333 74.922 > 49 95 As 0.58333333333333 0.91666666666667 0.83333333333333 74.922 > 49 96 As 0.83333333333333 0.91666666666667 0.83333333333333 74.922 > 49 97 Ba 0.16666666666667 0.08333333333333 0.00000000000000 137.327 > 97 98 Ba 0.41666666666667 0.08333333333333 0.00000000000000 137.327 > 97 99 Ba 0.66666666666667 0.08333333333333 0.00000000000000 137.327 > 97 100 Ba 0.91666666666667 0.08333333333333 0.00000000000000 137.327 > 97 101 Ba 0.16666666666667 0.33333333333333 0.00000000000000 137.327 > 97 102 Ba 0.41666666666667 0.33333333333333 0.00000000000000 137.327 > 97 103 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 > 97 104 Ba 0.91666666666667 0.33333333333333 0.00000000000000 137.327 > 97 105 Ba 0.16666666666667 0.58333333333333 0.00000000000000 137.327 > 97 106 Ba 0.41666666666667 0.58333333333333 0.00000000000000 137.327 > 97 107 Ba 0.66666666666667 0.58333333333333 0.00000000000000 137.327 > 97 108 Ba 0.91666666666667 0.58333333333333 0.00000000000000 137.327 > 97 109 Ba 0.16666666666667 0.83333333333333 0.00000000000000 137.327 > 97 110 Ba 0.41666666666667 0.83333333333333 0.00000000000000 137.327 > 97 111 Ba 0.66666666666667 0.83333333333333 0.00000000000000 137.327 > 97 112 Ba 0.91666666666667 0.83333333333333 0.00000000000000 137.327 > 97 113 Ba 0.16666666666667 0.08333333333333 0.33333333333333 137.327 > 97 114 Ba 0.41666666666667 0.08333333333333 0.33333333333333 137.327 > 97 115 Ba 0.66666666666667 0.08333333333333 0.33333333333333 137.327 > 97 116 Ba 0.91666666666667 0.08333333333333 0.33333333333333 137.327 > 97 117 Ba 0.16666666666667 0.33333333333333 0.33333333333333 137.327 > 97 118 Ba 0.41666666666667 0.33333333333333 0.33333333333333 137.327 > 97 119 Ba 0.66666666666667 0.33333333333333 0.33333333333333 137.327 > 97 120 Ba 0.91666666666667 0.33333333333333 0.33333333333333 137.327 > 97 121 Ba 0.16666666666667 0.58333333333333 0.33333333333333 137.327 > 97 122 Ba 0.41666666666667 0.58333333333333 0.33333333333333 137.327 > 97 123 Ba 0.66666666666667 0.58333333333333 0.33333333333333 137.327 > 97 124 Ba 0.91666666666667 0.58333333333333 0.33333333333333 137.327 > 97 125 Ba 0.16666666666667 0.83333333333333 0.33333333333333 137.327 > 97 126 Ba 0.41666666666667 0.83333333333333 0.33333333333333 137.327 > 97 127 Ba 0.66666666666667 0.83333333333333 0.33333333333333 137.327 > 97 128 Ba 0.91666666666667 0.83333333333333 0.33333333333333 137.327 > 97 129 Ba 0.16666666666667 0.08333333333333 0.66666666666667 137.327 > 97 130 Ba 0.41666666666667 0.08333333333333 0.66666666666667 137.327 > 97 131 Ba 0.66666666666667 0.08333333333333 0.66666666666667 137.327 > 97 132 Ba 0.91666666666667 0.08333333333333 0.66666666666667 137.327 > 97 133 Ba 0.16666666666667 0.33333333333333 0.66666666666667 137.327 > 97 134 Ba 0.41666666666667 0.33333333333333 0.66666666666667 137.327 > 97 135 Ba 0.66666666666667 0.33333333333333 0.66666666666667 137.327 > 97 136 Ba 0.91666666666667 0.33333333333333 0.66666666666667 137.327 > 97 137 Ba 0.16666666666667 0.58333333333333 0.66666666666667 137.327 > 97 138 Ba 0.41666666666667 0.58333333333333 0.66666666666667 137.327 > 97 139 Ba 0.66666666666667 0.58333333333333 0.66666666666667 137.327 > 97 140 Ba 0.91666666666667 0.58333333333333 0.66666666666667 137.327 > 97 141 Ba 0.16666666666667 0.83333333333333 0.66666666666667 137.327 > 97 142 Ba 0.41666666666667 0.83333333333333 0.66666666666667 137.327 > 97 143 Ba 0.66666666666667 0.83333333333333 0.66666666666667 137.327 > 97 144 Ba 0.91666666666667 0.83333333333333 0.66666666666667 137.327 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 9.6531051 -0.0000000 0.0000000 0.0000000 9.6531051 0.0000000 0.0000000 0.0000000 45.3917432 -------------------------- Born effective charges -------------------------- 1 Li 0.9957698 0.0000000 0.0000000 0.0000000 0.9957698 0.0000000 0.0000000 0.0000000 -0.0040015 2 As -3.1893956 -0.0000000 0.0000000 -0.0000000 -3.1893956 0.0000000 0.0000000 0.0000000 -5.8823552 3 Ba 2.1936258 0.0000000 0.0000000 0.0000000 2.1936258 0.0000000 0.0000000 0.0000000 5.8863567 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000012 (xxy) -0.00000012 (xxy) -0.00000012 (xyx) Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- 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 ... [ 13 13 11 ] 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.97, Number of G-points: 315, Lambda: 0.47 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/126) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126 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.927 ( 0.000 0.000 0.000) 0.000 3.032 ( 0.000 0.000 0.000) 0.000 3.032 ( 0.000 0.000 0.000) 0.000 9.134 ( 0.000 0.000 0.000) 0.000 10.763 ( 0.000 -0.000 0.000) 0.000 10.763 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/126) ======================= q-point: ( 0.08 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.391 ( 16.779 9.688 0.000) 19.375 0.485 ( 20.587 11.886 0.000) 23.771 0.812 ( 34.575 19.962 0.000) 39.924 1.951 ( 1.974 1.140 0.000) 2.280 3.035 ( 0.289 0.167 0.000) 0.334 3.947 ( 2.713 1.566 0.000) 3.133 9.134 ( 0.050 0.029 0.000) 0.057 10.759 ( -0.321 -0.185 0.000) 0.371 11.306 ( -0.151 -0.087 0.000) 0.175 ======================= Grid point 2 (3/126) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.760 ( 15.403 8.893 0.000) 17.786 0.928 ( 17.984 10.383 0.000) 20.766 1.558 ( 30.283 17.484 0.000) 34.968 2.007 ( 2.661 1.536 0.000) 3.073 3.047 ( 0.735 0.425 0.000) 0.849 4.032 ( 4.661 2.691 0.000) 5.382 9.137 ( 0.230 0.133 0.000) 0.266 10.748 ( -0.574 -0.332 0.000) 0.663 11.296 ( -0.863 -0.498 0.000) 0.996 ======================= Grid point 3 (4/126) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.089 ( 13.243 7.646 0.000) 15.292 1.303 ( 14.819 8.556 0.000) 17.111 2.061 ( 1.833 1.058 0.000) 2.117 2.177 ( 23.536 13.589 0.000) 27.178 3.063 ( 0.459 0.265 0.000) 0.530 4.155 ( 5.960 3.441 0.000) 6.882 9.146 ( 0.521 0.301 0.000) 0.602 10.735 ( -0.554 -0.320 0.000) 0.640 11.262 ( -2.147 -1.239 0.000) 2.479 ======================= Grid point 4 (5/126) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.360 ( 10.396 6.002 0.000) 12.004 1.609 ( 12.102 6.987 0.000) 13.974 2.085 ( 0.291 0.168 0.000) 0.336 2.625 ( 15.647 9.034 0.000) 18.068 3.057 ( -1.256 -0.725 0.000) 1.450 4.300 ( 6.583 3.801 0.000) 7.602 9.160 ( 0.736 0.425 0.000) 0.850 10.726 ( -0.220 -0.127 0.000) 0.254 11.199 ( -3.267 -1.886 0.000) 3.772 ======================= Grid point 5 (6/126) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.559 ( 6.897 3.982 0.000) 7.964 1.857 ( 9.401 5.428 0.000) 10.855 2.077 ( -0.827 -0.477 0.000) 0.954 2.898 ( 8.505 4.910 0.000) 9.820 3.004 ( -3.186 -1.839 0.000) 3.679 4.442 ( 5.432 3.136 0.000) 6.272 9.177 ( 0.667 0.385 0.000) 0.770 10.725 ( 0.108 0.063 0.000) 0.125 11.122 ( -3.141 -1.814 0.000) 3.627 ======================= Grid point 6 (7/126) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.669 ( 2.535 1.464 0.000) 2.928 2.019 ( 4.112 2.374 0.000) 4.748 2.058 ( -0.614 -0.355 0.000) 0.709 2.934 ( -2.129 -1.229 0.000) 2.458 3.024 ( 2.668 1.540 0.000) 3.081 4.531 ( 2.114 1.220 0.000) 2.440 9.188 ( 0.271 0.156 0.000) 0.312 10.728 ( 0.105 0.061 0.000) 0.122 11.068 ( -1.318 -0.761 0.000) 1.522 ======================= Grid point 14 (8/126) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.665 ( 9.174 15.889 0.000) 18.347 0.841 ( 11.841 20.509 0.000) 23.682 1.355 ( 17.819 30.864 0.000) 35.638 1.991 ( 1.536 2.660 0.000) 3.072 3.051 ( 0.824 1.427 0.000) 1.648 3.999 ( 2.066 3.578 0.000) 4.132 9.136 ( 0.104 0.180 0.000) 0.208 10.757 ( -0.046 -0.079 0.000) 0.092 11.295 ( -0.619 -1.072 0.000) 1.237 ======================= Grid point 15 (9/126) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.979 ( 10.653 12.469 0.000) 16.400 1.232 ( 10.632 18.354 0.000) 21.211 1.951 ( 20.566 20.806 0.000) 29.255 2.045 ( 1.725 2.117 0.000) 2.731 3.088 ( 0.293 3.719 0.000) 3.730 4.085 ( 4.481 2.558 0.000) 5.160 9.142 ( 0.323 0.451 0.000) 0.555 10.753 ( -0.770 0.689 0.000) 1.033 11.262 ( -1.000 -2.959 0.000) 3.124 ======================= Grid point 16 (10/126) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.264 ( 9.560 9.654 0.000) 13.587 1.561 ( 7.862 16.126 0.000) 17.941 2.083 ( 0.708 1.146 0.000) 1.347 2.442 ( 16.983 12.539 0.000) 21.111 3.122 ( -2.171 5.772 0.000) 6.167 4.203 ( 6.911 1.366 0.000) 7.045 9.156 ( 0.529 0.803 0.000) 0.962 10.744 ( -1.217 1.026 0.000) 1.592 11.202 ( -1.438 -4.679 0.000) 4.895 ======================= Grid point 17 (11/126) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.491 ( 7.149 6.928 0.000) 9.956 1.823 ( 5.704 13.421 0.000) 14.583 2.091 ( -0.635 0.479 0.000) 0.796 2.775 ( 11.383 5.777 0.000) 12.765 3.117 ( -5.854 6.975 0.000) 9.106 4.342 ( 7.889 0.518 0.000) 7.906 9.175 ( 0.502 1.065 0.000) 1.177 10.734 ( -0.888 0.671 0.000) 1.113 11.125 ( -1.537 -5.306 0.000) 5.524 ======================= Grid point 18 (12/126) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.640 ( 3.788 4.014 0.000) 5.519 2.016 ( 2.999 9.633 0.000) 10.089 2.078 ( -1.317 0.608 0.000) 1.451 2.955 ( 6.228 0.736 0.000) 6.272 3.069 ( -8.082 8.538 0.000) 11.756 4.464 ( 5.887 -0.827 0.000) 5.945 9.192 ( 0.127 1.077 0.000) 1.085 10.728 ( -0.223 0.068 0.000) 0.233 11.056 ( -0.474 -4.658 0.000) 4.682 ======================= Grid point 19 (13/126) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.694 ( -0.592 1.025 0.000) 1.184 2.068 ( -0.751 1.300 0.000) 1.502 2.094 ( -2.838 4.914 0.000) 5.674 3.010 ( 1.135 -1.967 0.000) 2.271 3.038 ( -5.925 10.260 0.000) 11.848 4.512 ( 1.787 -3.101 0.000) 3.579 9.198 ( -0.473 0.820 0.000) 0.947 10.727 ( 0.109 -0.189 0.000) 0.218 11.027 ( 1.865 -3.233 0.000) 3.732 ======================= Grid point 29 (14/126) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.229 ( 7.126 12.342 0.000) 14.252 1.590 ( 9.758 16.902 0.000) 19.517 2.083 ( 0.901 1.560 0.000) 1.802 2.346 ( 10.611 18.379 0.000) 21.222 3.183 ( 3.149 5.455 0.000) 6.299 4.138 ( 1.757 3.044 0.000) 3.515 9.156 ( 0.527 0.912 0.000) 1.053 10.766 ( 0.206 0.356 0.000) 0.411 11.184 ( -2.666 -4.617 0.000) 5.332 ======================= Grid point 30 (15/126) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.454 ( 6.389 9.179 0.000) 11.184 1.884 ( 5.733 15.448 0.000) 16.478 2.104 ( 0.086 0.874 0.000) 0.878 2.656 ( 9.563 8.768 0.000) 12.974 3.281 ( -1.185 9.667 0.000) 9.739 4.216 ( 5.536 0.339 0.000) 5.547 9.178 ( 0.645 1.338 0.000) 1.486 10.762 ( -1.142 0.412 0.000) 1.214 11.089 ( -1.845 -6.312 0.000) 6.576 ======================= Grid point 31 (16/126) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.621 ( 4.162 5.892 0.000) 7.214 2.093 ( 2.150 13.000 0.000) 13.177 2.106 ( -0.785 1.031 0.000) 1.296 2.844 ( 6.903 1.106 0.000) 6.991 3.323 ( -6.907 13.280 0.000) 14.969 4.326 ( 7.612 -1.979 0.000) 7.865 9.201 ( 0.359 1.553 0.000) 1.594 10.740 ( -1.214 -0.277 0.000) 1.246 11.004 ( -0.492 -6.403 0.000) 6.422 ======================= Grid point 32 (17/126) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.709 ( 1.010 2.735 0.000) 2.915 2.102 ( -1.172 1.729 0.000) 2.088 2.209 ( -1.928 9.554 0.000) 9.747 2.925 ( 4.286 -3.687 0.000) 5.654 3.322 ( -9.564 15.928 0.000) 18.578 4.412 ( 5.382 -4.400 0.000) 6.951 9.217 ( -0.311 1.412 0.000) 1.446 10.722 ( -0.174 -0.696 0.000) 0.717 10.953 ( 1.418 -5.300 0.000) 5.487 ======================= Grid point 43 (18/126) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.615 ( 3.897 6.751 0.000) 7.795 2.119 ( 0.392 0.678 0.000) 0.783 2.150 ( 6.477 11.218 0.000) 12.953 2.781 ( 2.249 3.896 0.000) 4.498 3.473 ( 5.134 8.892 0.000) 10.268 4.230 ( 0.777 1.345 0.000) 1.554 9.207 ( 0.859 1.488 0.000) 1.718 10.753 ( -0.738 -1.278 0.000) 1.476 10.971 ( -3.014 -5.220 0.000) 6.028 ======================= Grid point 44 (19/126) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 510 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.715 ( 1.835 3.467 0.000) 3.923 2.131 ( -0.178 1.398 0.000) 1.409 2.324 ( 1.268 10.346 0.000) 10.424 2.821 ( 2.503 -3.220 0.000) 4.078 3.605 ( -3.257 14.103 0.000) 14.475 4.269 ( 4.656 -3.771 0.000) 5.992 9.232 ( 0.323 1.422 0.000) 1.458 10.721 ( -0.923 -1.498 0.000) 1.760 10.891 ( -0.521 -4.629 0.000) 4.658 ======================= Grid point 45 (20/126) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.747 ( -0.692 1.199 0.000) 1.384 2.137 ( -0.896 1.551 0.000) 1.791 2.388 ( -4.828 8.359 0.000) 9.653 2.823 ( 3.505 -6.071 0.000) 7.010 3.645 ( -8.819 15.271 0.000) 17.635 4.294 ( 4.146 -7.186 0.000) 8.296 9.242 ( -0.579 1.004 0.000) 1.159 10.704 ( 0.542 -0.941 0.000) 1.086 10.863 ( 1.928 -3.343 0.000) 3.859 ======================= Grid point 58 (21/126) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.759 ( 0.603 1.044 0.000) 1.206 2.154 ( 0.404 0.700 0.000) 0.808 2.495 ( 3.589 6.216 0.000) 7.178 2.737 ( -2.405 -4.165 0.000) 4.809 3.836 ( 4.447 7.703 0.000) 8.895 4.187 ( -2.087 -3.615 0.000) 4.174 9.253 ( 0.334 0.578 0.000) 0.668 10.693 ( -0.537 -0.931 0.000) 1.075 10.830 ( -0.905 -1.568 0.000) 1.810 ======================= Grid point 181 (22/126) ======================= q-point: ( 0.00 0.00 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.360 ( 0.000 0.000 18.732) 18.732 0.360 ( 0.000 0.000 18.732) 18.732 0.792 ( 0.000 0.000 39.166) 39.166 3.002 ( 0.000 0.000 -2.938) 2.938 3.002 ( 0.000 0.000 -2.938) 2.938 3.650 ( 0.000 0.000 10.673) 10.673 9.133 ( 0.000 0.000 -0.084) 0.084 10.760 ( 0.000 0.000 -0.201) 0.201 10.760 ( 0.000 0.000 -0.201) 0.201 ======================= Grid point 182 (23/126) ======================= q-point: ( 0.08 0.00 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.605 ( 16.544 9.552 11.174) 22.131 0.680 ( 15.390 8.885 15.527) 23.598 1.039 ( 21.967 12.683 24.882) 35.532 2.844 ( -10.859 -6.269 18.103) 22.021 3.006 ( 0.387 0.224 -2.810) 2.845 3.724 ( 5.171 2.985 -1.758) 6.224 9.135 ( 0.200 0.115 0.089) 0.247 10.756 ( -0.357 -0.206 -0.243) 0.479 10.864 ( 7.357 4.247 -8.010) 11.676 ======================= Grid point 183 (24/126) ======================= q-point: ( 0.15 0.00 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.983 ( 12.335 7.121 17.113) 22.265 0.996 ( 16.777 9.686 6.800) 20.531 1.631 ( 26.440 15.265 9.121) 31.863 2.613 ( -8.343 -4.817 30.579) 32.061 3.021 ( 0.845 0.488 -2.541) 2.721 3.856 ( 6.270 3.620 -9.921) 12.282 9.143 ( 0.497 0.287 0.560) 0.802 10.745 ( -0.621 -0.358 -0.338) 0.793 11.021 ( 5.454 3.149 -12.743) 14.215 ======================= Grid point 184 (25/126) ======================= q-point: ( 0.23 0.00 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.252 ( 11.099 6.408 13.578) 18.671 1.352 ( 14.233 8.217 4.961) 17.167 2.169 ( 19.211 11.092 4.141) 22.567 2.477 ( -2.996 -1.730 25.899) 26.129 3.039 ( 0.563 0.325 -2.270) 2.362 4.009 ( 7.138 4.121 -11.148) 13.864 9.158 ( 0.853 0.492 1.230) 1.576 10.730 ( -0.597 -0.344 -0.435) 0.815 11.101 ( 1.652 0.954 -10.504) 10.676 ======================= Grid point 185 (26/126) ======================= q-point: ( 0.31 0.00 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.482 ( 8.893 5.134 10.336) 14.570 1.647 ( 11.589 6.691 3.750) 13.897 2.335 ( -1.591 -0.919 19.529) 19.615 2.623 ( 12.400 7.159 0.438) 14.325 3.036 ( -1.056 -0.610 -1.901) 2.258 4.178 ( 7.406 4.276 -10.197) 13.308 9.181 ( 1.060 0.612 1.966) 2.316 10.720 ( -0.262 -0.151 -0.526) 0.607 11.105 ( -1.035 -0.597 -7.084) 7.184 ======================= Grid point 186 (27/126) ======================= q-point: ( 0.38 0.00 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.653 ( 5.902 3.407 8.058) 10.553 1.882 ( 8.816 5.090 2.488) 10.479 2.290 ( -1.883 -1.087 17.765) 17.897 2.851 ( 7.270 4.197 -4.293) 9.428 2.989 ( -2.826 -1.631 -1.220) 3.484 4.333 ( 5.814 3.357 -9.296) 11.467 9.204 ( 0.895 0.517 2.584) 2.783 10.718 ( 0.070 0.041 -0.616) 0.621 11.068 ( -1.881 -1.086 -4.381) 4.890 ======================= Grid point 187 (28/126) ======================= q-point: ( 0.46 0.00 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.746 ( 2.127 1.228 6.852) 7.279 2.032 ( 3.765 2.174 1.328) 4.545 2.258 ( -0.782 -0.452 16.642) 16.667 2.928 ( -1.870 -1.080 -0.404) 2.197 2.959 ( 2.288 1.321 -5.923) 6.485 4.428 ( 2.202 1.271 -8.971) 9.324 9.218 ( 0.349 0.201 2.922) 2.949 10.720 ( 0.089 0.051 -0.679) 0.687 11.033 ( -0.892 -0.515 -2.885) 3.064 ======================= Grid point 195 (29/126) ======================= q-point: ( 0.08 0.08 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.904 ( 7.263 12.579 17.421) 22.682 0.914 ( 10.831 18.760 7.271) 22.850 1.457 ( 15.149 26.238 12.123) 32.633 2.663 ( -5.675 -9.829 29.016) 31.156 3.026 ( 1.005 1.741 -2.297) 3.053 3.813 ( 3.151 5.457 -8.581) 10.646 9.140 ( 0.237 0.411 0.404) 0.623 10.753 ( -0.077 -0.134 -0.321) 0.357 10.979 ( 3.431 5.943 -12.456) 14.222 ======================= Grid point 196 (30/126) ======================= q-point: ( 0.15 0.08 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.159 ( 8.892 10.194 14.563) 19.876 1.281 ( 10.336 17.270 5.013) 20.742 1.978 ( 17.940 18.366 5.151) 26.185 2.501 ( -3.821 -5.529 28.137) 28.929 3.069 ( 0.232 4.266 -1.661) 4.584 3.932 ( 5.377 3.779 -11.020) 12.831 9.152 ( 0.560 0.681 0.941) 1.289 10.748 ( -0.808 0.636 -0.430) 1.114 11.067 ( 2.491 1.066 -11.610) 11.922 ======================= Grid point 197 (31/126) ======================= q-point: ( 0.23 0.08 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.398 ( 8.225 7.967 11.173) 15.999 1.597 ( 7.749 15.217 3.704) 17.474 2.329 ( 3.393 2.279 14.638) 15.198 2.489 ( 9.211 5.559 8.279) 13.575 3.106 ( -2.385 6.319 -1.369) 6.892 4.074 ( 7.676 2.290 -10.539) 13.238 9.172 ( 0.834 0.969 1.587) 2.038 10.738 ( -1.238 0.980 -0.523) 1.663 11.087 ( 0.841 -2.426 -8.331) 8.717 ======================= Grid point 198 (32/126) ======================= q-point: ( 0.31 0.08 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.590 ( 6.277 5.529 8.429) 11.875 1.848 ( 5.531 12.560 2.569) 13.962 2.309 ( -1.603 -1.349 18.022) 18.143 2.743 ( 9.828 4.634 -2.871) 11.239 3.103 ( -5.878 7.377 -1.144) 9.501 4.232 ( 8.289 1.201 -9.392) 12.584 9.198 ( 0.807 1.114 2.202) 2.597 10.728 ( -0.904 0.662 -0.577) 1.260 11.058 ( -0.155 -4.006 -5.310) 6.653 ======================= Grid point 199 (33/126) ======================= q-point: ( 0.38 0.08 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.717 ( 3.454 2.882 6.671) 8.046 2.029 ( 2.775 8.889 1.386) 9.414 2.275 ( -1.339 -0.184 16.487) 16.543 2.897 ( 5.452 0.348 -5.325) 7.629 3.061 ( -7.830 8.777 -0.602) 11.777 4.362 ( 5.893 -0.362 -8.774) 10.576 9.219 ( 0.356 0.991 2.629) 2.832 10.722 ( -0.260 0.120 -0.599) 0.664 11.015 ( 0.196 -4.003 -3.289) 5.185 ======================= Grid point 200 (34/126) ======================= q-point: ( 0.46 0.08 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.761 ( -0.192 0.333 6.058) 6.070 2.101 ( -2.599 4.501 0.773) 5.255 2.262 ( -0.500 0.865 16.105) 16.136 2.944 ( 1.413 -2.449 -6.018) 6.649 3.033 ( -6.230 10.788 -0.258) 12.460 4.413 ( 1.588 -2.756 -8.641) 9.208 9.227 ( -0.364 0.631 2.775) 2.869 10.720 ( 0.060 -0.105 -0.606) 0.618 10.996 ( 1.786 -3.096 -2.569) 4.401 ======================= Grid point 210 (35/126) ======================= q-point: ( 0.15 0.15 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.366 ( 5.957 10.317 11.406) 16.493 1.621 ( 9.320 16.144 3.168) 18.908 2.295 ( 5.607 9.712 8.482) 14.061 2.441 ( 1.643 2.846 15.318) 15.667 3.173 ( 3.310 5.733 -0.725) 6.659 4.008 ( 2.381 4.124 -10.640) 11.657 9.170 ( 0.666 1.153 1.398) 1.930 10.760 ( 0.194 0.336 -0.489) 0.624 11.060 ( -0.937 -1.623 -8.821) 9.018 ======================= Grid point 211 (36/126) ======================= q-point: ( 0.23 0.15 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.554 ( 5.412 7.478 8.414) 12.490 1.902 ( 5.580 14.547 1.854) 15.691 2.318 ( -0.979 -1.223 17.655) 17.725 2.637 ( 8.138 7.023 -1.736) 10.888 3.272 ( -1.368 9.821 -0.587) 9.933 4.106 ( 6.101 1.210 -9.407) 11.278 9.196 ( 0.833 1.422 1.827) 2.461 10.756 ( -1.178 0.481 -0.494) 1.365 11.013 ( -0.536 -4.589 -5.882) 7.480 ======================= Grid point 212 (37/126) ======================= q-point: ( 0.31 0.15 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.692 ( 3.614 4.519 6.166) 8.456 2.099 ( 2.090 11.999 0.752) 12.203 2.295 ( -0.950 0.089 15.951) 15.980 2.796 ( 6.091 0.545 -4.489) 7.586 3.313 ( -6.894 13.182 -0.762) 14.895 4.230 ( 7.769 -1.268 -8.310) 11.447 9.224 ( 0.565 1.452 2.144) 2.650 10.735 ( -1.243 -0.125 -0.434) 1.323 10.960 ( 0.286 -5.470 -3.647) 6.581 ======================= Grid point 213 (38/126) ======================= q-point: ( 0.38 0.15 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.764 ( 1.047 1.853 5.043) 5.474 2.205 ( -1.742 8.359 0.297) 8.544 2.287 ( -0.937 1.554 14.945) 15.055 2.864 ( 3.895 -3.585 -5.529) 7.655 3.312 ( -9.306 15.553 -0.735) 18.140 4.321 ( 5.164 -3.791 -7.894) 10.166 9.241 ( -0.128 1.176 2.299) 2.585 10.718 ( -0.250 -0.497 -0.346) 0.655 10.923 ( 1.635 -4.931 -2.468) 5.751 ======================= Grid point 224 (39/126) ======================= q-point: ( 0.23 0.23 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.683 ( 3.126 5.414 5.962) 8.639 2.150 ( 5.962 10.326 0.292) 11.927 2.305 ( -0.010 -0.017 15.741) 15.741 2.737 ( 1.776 3.075 -4.110) 5.431 3.463 ( 4.977 8.621 -0.713) 9.980 4.135 ( 1.191 2.063 -8.209) 8.548 9.226 ( 0.863 1.495 1.906) 2.572 10.750 ( -0.655 -1.135 -0.239) 1.333 10.923 ( -2.368 -4.101 -3.901) 6.135 ======================= Grid point 225 (40/126) ======================= q-point: ( 0.31 0.23 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.764 ( 1.558 2.708 4.513) 5.488 2.285 ( 0.365 3.614 7.548) 8.377 2.335 ( 0.598 6.912 6.342) 9.400 2.767 ( 2.234 -3.271 -4.936) 6.329 3.589 ( -3.344 13.663 -1.294) 14.126 4.185 ( 4.804 -3.184 -7.280) 9.285 9.252 ( 0.385 1.267 1.878) 2.298 10.720 ( -0.941 -1.241 0.028) 1.557 10.860 ( -0.155 -4.159 -2.617) 4.917 ======================= Grid point 226 (41/126) ======================= q-point: ( 0.38 0.23 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.790 ( -0.504 0.872 4.110) 4.232 2.298 ( -0.829 1.436 11.890) 12.005 2.385 ( -4.340 7.516 1.192) 8.760 2.766 ( 3.356 -5.814 -5.045) 8.397 3.626 ( -8.546 14.797 -1.579) 17.160 4.215 ( 3.859 -6.689 -6.943) 10.385 9.262 ( -0.462 0.801 1.848) 2.067 10.705 ( 0.408 -0.708 0.137) 0.829 10.838 ( 1.878 -3.256 -2.167) 4.339 ======================= Grid point 239 (42/126) ======================= q-point: ( 0.31 0.31 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.799 ( 0.501 0.867 3.863) 3.991 2.312 ( 0.309 0.536 12.558) 12.573 2.487 ( 3.578 6.197 -0.203) 7.158 2.682 ( -2.545 -4.408 -4.906) 7.070 3.813 ( 4.391 7.605 -1.769) 8.958 4.112 ( -2.024 -3.505 -6.770) 7.887 9.270 ( 0.291 0.504 1.639) 1.739 10.696 ( -0.477 -0.827 0.429) 1.047 10.804 ( -0.816 -1.413 -2.169) 2.714 ======================= Grid point 362 (43/126) ======================= q-point: ( 0.00 0.00 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.757 ( 0.000 0.000 20.887) 20.887 0.757 ( 0.000 0.000 20.887) 20.887 1.563 ( 0.000 -0.000 37.709) 37.709 2.918 ( 0.000 0.000 -5.435) 5.435 2.918 ( 0.000 -0.000 -5.435) 5.435 3.908 ( 0.000 0.000 13.412) 13.412 9.132 ( 0.000 0.000 -0.010) 0.010 10.756 ( 0.000 0.000 -0.158) 0.158 10.756 ( 0.000 0.000 -0.158) 0.158 ======================= Grid point 363 (44/126) ======================= q-point: ( 0.08 0.00 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.899 ( 11.074 6.393 17.497) 21.672 0.972 ( 14.867 8.584 15.824) 23.348 1.662 ( 8.781 5.070 33.623) 35.119 2.926 ( 0.723 0.417 -5.143) 5.210 2.952 ( 2.534 1.463 -1.709) 3.388 3.850 ( -4.438 -2.562 10.312) 11.515 9.138 ( 0.599 0.346 0.268) 0.742 10.751 ( -0.447 -0.258 -0.221) 0.561 10.793 ( 2.958 1.708 -1.416) 3.697 ======================= Grid point 364 (45/126) ======================= q-point: ( 0.15 0.00 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.196 ( 13.842 7.992 12.881) 20.528 1.291 ( 11.860 6.847 13.271) 19.069 1.954 ( 15.748 9.092 22.032) 28.567 2.950 ( 1.301 0.751 -4.436) 4.684 2.993 ( -0.053 -0.031 9.479) 9.479 3.759 ( -1.953 -1.128 -0.121) 2.258 9.159 ( 1.196 0.691 1.032) 1.724 10.737 ( -0.737 -0.425 -0.363) 0.925 10.874 ( 3.777 2.181 -3.633) 5.677 ======================= Grid point 365 (46/126) ======================= q-point: ( 0.23 0.00 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.501 ( 12.461 7.194 9.720) 17.364 1.519 ( 8.529 4.924 11.515) 15.152 2.313 ( 14.402 8.315 10.608) 19.726 2.921 ( -5.826 -3.363 17.466) 18.717 2.979 ( 1.088 0.628 -3.530) 3.747 3.804 ( 5.683 3.281 -8.626) 10.838 9.192 ( 1.712 0.989 2.112) 2.893 10.720 ( -0.704 -0.406 -0.512) 0.960 10.949 ( 2.559 1.477 -4.649) 5.508 ======================= Grid point 366 (47/126) ======================= q-point: ( 0.31 0.00 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.691 ( 6.651 3.840 9.062) 11.879 1.758 ( 10.005 5.777 7.232) 13.630 2.579 ( 8.193 4.730 2.710) 9.841 2.777 ( -5.321 -3.072 15.493) 16.667 2.991 ( -0.253 -0.146 -2.342) 2.360 3.975 ( 8.243 4.759 -8.708) 12.900 9.235 ( 1.886 1.089 3.273) 3.932 10.707 ( -0.367 -0.212 -0.653) 0.778 10.987 ( 0.777 0.449 -4.059) 4.157 ======================= Grid point 367 (48/126) ======================= q-point: ( 0.38 0.00 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.822 ( 4.676 2.700 7.546) 9.279 1.955 ( 7.137 4.121 4.777) 9.526 2.605 ( -3.527 -2.036 6.942) 8.049 2.786 ( 4.088 2.360 2.831) 5.504 2.967 ( -1.659 -0.958 -0.694) 2.038 4.146 ( 6.239 3.602 -7.804) 10.621 9.274 ( 1.471 0.849 4.232) 4.560 10.703 ( -0.026 -0.015 -0.790) 0.791 10.991 ( -0.268 -0.155 -2.658) 2.676 ======================= Grid point 368 (49/126) ======================= q-point: ( 0.46 0.00 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.898 ( 1.783 1.029 7.141) 7.431 2.073 ( 2.828 1.633 2.751) 4.270 2.542 ( -1.531 -0.884 6.262) 6.507 2.853 ( 1.505 0.869 -0.850) 1.935 2.930 ( -1.125 -0.649 0.909) 1.585 4.245 ( 2.247 1.297 -7.637) 8.066 9.298 ( 0.548 0.317 4.760) 4.802 10.704 ( 0.048 0.027 -0.887) 0.889 10.983 ( -0.246 -0.142 -1.551) 1.577 ======================= Grid point 376 (50/126) ======================= q-point: ( 0.08 0.08 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.126 ( 8.558 14.823 13.615) 21.871 1.208 ( 7.434 12.876 13.544) 20.112 1.858 ( 8.214 14.227 25.696) 30.499 2.930 ( -0.161 -0.280 -1.673) 1.704 3.004 ( 1.889 3.271 2.908) 4.767 3.774 ( -2.265 -3.923 3.368) 5.645 9.152 ( 0.591 1.023 0.781) 1.416 10.746 ( -0.156 -0.269 -0.338) 0.459 10.846 ( 1.936 3.353 -3.077) 4.945 ======================= Grid point 377 (51/126) ======================= q-point: ( 0.15 0.08 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.424 ( 8.119 12.397 10.694) 18.275 1.446 ( 8.532 9.912 11.396) 17.347 2.178 ( 12.065 12.533 14.692) 22.771 2.895 ( -1.392 -5.063 8.270) 9.796 3.060 ( -0.683 5.406 3.888) 6.693 3.749 ( 2.433 0.620 -6.752) 7.203 9.179 ( 1.175 1.290 1.658) 2.407 10.738 ( -0.891 0.514 -0.494) 1.141 10.912 ( 2.655 1.302 -4.439) 5.334 ======================= Grid point 378 (52/126) ======================= q-point: ( 0.23 0.08 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.620 ( 6.406 5.545 9.631) 12.827 1.710 ( 7.120 12.559 7.435) 16.238 2.474 ( 10.075 7.449 5.574) 13.714 2.789 ( -3.121 -7.011 13.843) 15.828 3.088 ( -3.740 7.905 0.562) 8.763 3.866 ( 8.163 2.758 -9.170) 12.583 9.216 ( 1.607 1.404 2.690) 3.434 10.726 ( -1.296 0.905 -0.611) 1.695 10.955 ( 2.145 -1.031 -4.378) 4.983 ======================= Grid point 379 (53/126) ======================= q-point: ( 0.31 0.08 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.764 ( 5.174 3.512 7.691) 9.913 1.926 ( 4.834 10.140 5.035) 12.310 2.619 ( -0.169 -1.574 3.733) 4.055 2.722 ( 3.787 -0.319 8.167) 9.008 3.085 ( -6.081 8.615 -0.261) 10.548 4.044 ( 8.812 2.157 -7.913) 12.038 9.258 ( 1.574 1.237 3.645) 4.158 10.714 ( -0.980 0.697 -0.683) 1.383 10.966 ( 1.347 -2.537 -3.211) 4.309 ======================= Grid point 380 (54/126) ======================= q-point: ( 0.38 0.08 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.863 ( 3.303 1.453 6.855) 7.747 2.072 ( 2.143 6.907 2.841) 7.770 2.560 ( -1.971 -1.893 6.054) 6.642 2.806 ( 3.858 -0.232 0.444) 3.891 3.057 ( -7.212 9.347 0.530) 11.818 4.185 ( 5.736 0.341 -7.371) 9.347 9.290 ( 0.936 0.766 4.298) 4.465 10.707 ( -0.392 0.289 -0.739) 0.885 10.958 ( 1.166 -3.043 -1.863) 3.754 ======================= Grid point 381 (55/126) ======================= q-point: ( 0.46 0.08 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.899 ( 0.380 -0.658 6.733) 6.776 2.126 ( -2.049 3.548 1.803) 4.476 2.533 ( 0.003 -0.005 5.383) 5.383 2.837 ( 1.168 -2.024 -0.653) 2.427 3.039 ( -6.128 10.612 1.063) 12.300 4.237 ( 1.235 -2.145 -7.343) 7.749 9.302 ( -0.083 0.143 4.524) 4.527 10.706 ( -0.076 0.132 -0.768) 0.783 10.952 ( 1.655 -2.871 -1.273) 3.550 ======================= Grid point 391 (56/126) ======================= q-point: ( 0.15 0.15 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.592 ( 4.287 7.425 9.801) 13.022 1.720 ( 7.988 13.836 6.703) 17.326 2.417 ( 6.341 10.983 7.137) 14.553 2.778 ( -3.585 -6.209 13.102) 14.935 3.177 ( 3.394 5.878 1.981) 7.071 3.796 ( 2.458 4.258 -9.521) 10.715 9.209 ( 1.026 1.777 2.406) 3.162 10.749 ( 0.195 0.337 -0.550) 0.674 10.923 ( -0.064 -0.111 -4.469) 4.471 ======================= Grid point 392 (57/126) ======================= q-point: ( 0.23 0.15 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.726 ( 4.059 4.991 7.568) 9.933 1.961 ( 4.989 11.945 4.053) 13.564 2.590 ( 3.817 3.472 -0.814) 5.224 2.676 ( -0.909 -3.523 14.304) 14.760 3.271 ( -2.070 10.047 0.674) 10.280 3.915 ( 7.142 2.386 -8.218) 11.146 9.247 ( 1.302 1.631 3.076) 3.717 10.746 ( -1.275 0.762 -0.481) 1.562 10.913 ( 0.816 -2.977 -3.545) 4.701 ======================= Grid point 393 (58/126) ======================= q-point: ( 0.31 0.15 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.827 ( 3.134 2.746 6.330) 7.578 2.125 ( 1.808 9.458 1.864) 9.808 2.572 ( -1.752 -2.115 4.062) 4.903 2.729 ( 4.337 -0.027 3.869) 5.812 3.305 ( -6.895 12.955 0.104) 14.675 4.063 ( 7.991 -0.184 -6.878) 10.545 9.282 ( 1.076 1.183 3.543) 3.887 10.726 ( -1.398 0.344 -0.395) 1.493 10.895 ( 1.426 -4.283 -2.182) 5.014 ======================= Grid point 394 (59/126) ======================= q-point: ( 0.38 0.15 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.883 ( 1.321 0.709 6.032) 6.216 2.209 ( -1.425 6.732 0.390) 6.892 2.537 ( -0.654 -0.633 3.292) 3.415 2.779 ( 2.809 -2.335 1.999) 4.164 3.304 ( -8.735 14.700 0.143) 17.101 4.161 ( 4.714 -2.796 -6.575) 8.560 9.304 ( 0.337 0.561 3.764) 3.821 10.711 ( -0.512 0.043 -0.351) 0.622 10.881 ( 2.034 -4.402 -1.190) 4.993 ======================= Grid point 405 (60/126) ======================= q-point: ( 0.23 0.23 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.813 ( 2.176 3.769 6.087) 7.483 2.164 ( 4.824 8.355 1.116) 9.712 2.577 ( -1.357 -2.351 1.273) 2.998 2.685 ( 1.160 2.010 7.491) 7.842 3.457 ( 4.604 7.974 0.131) 9.209 3.969 ( 1.869 3.237 -6.868) 7.819 9.279 ( 0.863 1.495 3.178) 3.616 10.748 ( -0.375 -0.649 0.081) 0.754 10.852 ( -1.592 -2.757 -2.616) 4.121 ======================= Grid point 406 (61/126) ======================= q-point: ( 0.31 0.23 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.872 ( 1.387 1.852 5.575) 6.036 2.289 ( 0.909 7.009 -1.069) 7.148 2.540 ( -0.810 -1.354 -0.279) 1.602 2.714 ( 1.847 -1.263 6.608) 6.976 3.568 ( -3.613 12.809 -0.724) 13.328 4.038 ( 5.095 -2.347 -5.891) 8.134 9.303 ( 0.535 0.856 3.099) 3.259 10.725 ( -1.209 -0.393 0.421) 1.339 10.812 ( 0.665 -3.708 -1.605) 4.095 ======================= Grid point 407 (62/126) ======================= q-point: ( 0.38 0.23 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.894 ( -0.231 0.399 5.546) 5.565 2.332 ( -3.152 5.457 -1.939) 6.593 2.525 ( 0.438 -0.760 -1.408) 1.659 2.724 ( 1.613 -2.795 7.095) 7.794 3.598 ( -8.029 13.902 -1.086) 16.090 4.074 ( 3.427 -5.940 -5.601) 8.855 9.312 ( -0.159 0.275 3.038) 3.054 10.711 ( -0.064 0.111 0.436) 0.454 10.800 ( 1.969 -3.414 -1.077) 4.086 ======================= Grid point 420 (63/126) ======================= q-point: ( 0.31 0.31 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.898 ( 0.399 0.691 5.410) 5.468 2.398 ( 2.113 3.660 -2.774) 5.055 2.511 ( -0.861 -1.492 -5.639) 5.896 2.691 ( -0.438 -0.758 11.556) 11.589 3.783 ( 4.651 8.055 -1.497) 9.421 3.971 ( -2.342 -4.057 -5.346) 7.108 9.315 ( 0.177 0.307 2.701) 2.725 10.715 ( -0.271 -0.469 1.402) 1.503 10.761 ( -0.737 -1.276 -1.687) 2.240 ======================= Grid point 543 (64/126) ======================= q-point: ( 0.00 0.00 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.189 ( 0.000 0.000 21.914) 21.914 1.189 ( 0.000 0.000 21.914) 21.914 2.298 ( 0.000 0.000 35.560) 35.560 2.786 ( 0.000 0.000 -7.769) 7.769 2.786 ( 0.000 0.000 -7.769) 7.769 4.128 ( 0.000 0.000 7.442) 7.442 9.133 ( 0.000 0.000 0.207) 0.207 10.755 ( 0.000 0.000 0.058) 0.058 10.755 ( 0.000 0.000 0.058) 0.058 ======================= Grid point 544 (65/126) ======================= q-point: ( 0.08 0.00 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.281 ( 7.461 4.308 20.043) 21.817 1.327 ( 10.218 5.900 19.203) 22.538 2.328 ( 2.837 1.638 32.387) 32.552 2.802 ( 1.377 0.795 -7.277) 7.449 2.867 ( 6.827 3.942 -6.004) 9.909 4.031 ( -8.167 -4.715 6.230) 11.302 9.146 ( 1.093 0.631 0.487) 1.353 10.749 ( -0.548 -0.316 0.001) 0.632 10.778 ( 1.881 1.086 -0.287) 2.191 ======================= Grid point 545 (66/126) ======================= q-point: ( 0.15 0.00 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.494 ( 10.387 5.997 16.234) 20.184 1.561 ( 8.929 5.155 14.458) 17.757 2.430 ( 5.716 3.300 23.951) 24.844 2.846 ( 2.385 1.377 -5.922) 6.531 3.070 ( 9.790 5.652 -0.291) 11.308 3.794 ( -11.071 -6.392 1.913) 12.926 9.182 ( 2.068 1.194 1.265) 2.703 10.732 ( -0.873 -0.504 -0.132) 1.017 10.833 ( 2.758 1.592 -0.974) 3.330 ======================= Grid point 546 (67/126) ======================= q-point: ( 0.23 0.00 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.723 ( 5.545 3.201 9.867) 11.763 1.727 ( 9.632 5.561 12.349) 16.619 2.556 ( 4.505 2.601 12.059) 13.133 2.904 ( 2.486 1.435 -3.925) 4.863 3.191 ( -2.628 -1.517 10.085) 10.532 3.662 ( 3.177 1.834 -5.853) 6.907 9.238 ( 2.767 1.598 2.360) 3.972 10.712 ( -0.833 -0.481 -0.280) 1.002 10.895 ( 2.454 1.417 -1.297) 3.117 ======================= Grid point 547 (68/126) ======================= q-point: ( 0.31 0.00 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.833 ( 4.396 2.538 5.766) 7.682 1.923 ( 7.402 4.274 8.753) 12.234 2.606 ( -0.623 -0.360 -0.103) 0.727 2.951 ( 1.507 0.870 -1.479) 2.284 3.048 ( -6.024 -3.478 10.890) 12.922 3.854 ( 10.049 5.802 -3.330) 12.072 9.304 ( 2.880 1.663 3.515) 4.839 10.696 ( -0.495 -0.286 -0.423) 0.711 10.941 ( 1.571 0.907 -0.852) 2.005 ======================= Grid point 548 (69/126) ======================= q-point: ( 0.38 0.00 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.928 ( 3.897 2.250 3.231) 5.539 2.063 ( 4.757 2.746 5.569) 7.822 2.539 ( -4.218 -2.435 -7.254) 8.737 2.971 ( 0.264 0.152 1.106) 1.148 2.972 ( -1.177 -0.679 9.511) 9.607 4.051 ( 6.761 3.904 -1.701) 7.990 9.364 ( 2.167 1.251 4.468) 5.121 10.689 ( -0.140 -0.081 -0.562) 0.584 10.967 ( 0.770 0.444 0.119) 0.897 ======================= Grid point 549 (70/126) ======================= q-point: ( 0.46 0.00 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.999 ( 1.856 1.072 2.669) 3.423 2.137 ( 1.689 0.975 3.439) 3.954 2.459 ( -2.130 -1.230 -9.486) 9.800 2.963 ( -0.038 -0.022 7.291) 7.291 2.970 ( -0.133 -0.077 3.000) 3.004 4.155 ( 2.256 1.302 -1.442) 2.977 9.398 ( 0.797 0.460 5.008) 5.092 10.688 ( -0.001 -0.001 -0.657) 0.657 10.978 ( 0.224 0.129 0.908) 0.944 ======================= Grid point 557 (71/126) ======================= q-point: ( 0.08 0.08 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.439 ( 6.199 10.737 17.071) 21.098 1.498 ( 5.560 9.631 15.843) 19.356 2.396 ( 2.932 5.078 26.719) 27.354 2.836 ( 1.442 2.497 -6.142) 6.785 3.008 ( 5.589 9.681 -2.638) 11.485 3.862 ( -6.580 -11.397 3.557) 13.632 9.170 ( 1.035 1.793 1.012) 2.305 10.741 ( -0.249 -0.432 -0.116) 0.512 10.812 ( 1.290 2.235 -0.773) 2.694 ======================= Grid point 558 (72/126) ======================= q-point: ( 0.15 0.08 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.659 ( 5.632 7.597 12.315) 15.527 1.675 ( 6.800 8.198 12.351) 16.310 2.510 ( 4.151 4.166 16.656) 17.663 2.892 ( 1.651 2.828 -3.601) 4.867 3.199 ( 4.056 6.238 5.633) 9.332 3.651 ( -3.817 -6.156 -3.429) 8.015 9.215 ( 1.935 2.069 1.917) 3.420 10.731 ( -0.969 0.362 -0.272) 1.069 10.860 ( 2.281 0.962 -1.211) 2.756 ======================= Grid point 559 (73/126) ======================= q-point: ( 0.23 0.08 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.782 ( 4.267 3.031 7.309) 8.990 1.881 ( 5.693 9.049 9.288) 14.162 2.589 ( 1.759 0.132 4.656) 4.979 2.932 ( 0.744 0.287 1.620) 1.805 3.176 ( -7.389 2.033 7.208) 10.520 3.725 ( 10.091 4.110 -4.730) 11.878 9.275 ( 2.527 1.960 2.961) 4.358 10.716 ( -1.362 0.846 -0.378) 1.647 10.904 ( 2.503 -0.565 -1.072) 2.781 ======================= Grid point 560 (74/126) ======================= q-point: ( 0.31 0.08 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.878 ( 4.111 2.140 4.034) 6.145 2.040 ( 3.561 6.815 6.000) 9.753 2.561 ( -2.212 -3.881 -4.947) 6.665 2.931 ( 0.714 -1.384 7.612) 7.770 3.108 ( -6.383 7.978 2.808) 10.596 3.945 ( 9.669 3.394 -1.879) 10.419 9.336 ( 2.474 1.393 3.896) 4.821 10.703 ( -1.105 0.775 -0.448) 1.422 10.933 ( 2.210 -1.689 -0.275) 2.795 ======================= Grid point 561 (75/126) ======================= q-point: ( 0.38 0.08 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.961 ( 3.348 1.065 2.759) 4.467 2.137 ( 1.338 4.304 3.424) 5.660 2.477 ( -2.423 -3.311 -8.708) 9.626 2.939 ( 1.083 -0.526 7.648) 7.742 3.088 ( -6.151 9.106 2.520) 11.274 4.099 ( 5.546 1.006 -1.174) 5.757 9.381 ( 1.638 0.484 4.536) 4.847 10.694 ( -0.574 0.483 -0.515) 0.910 10.948 ( 1.856 -2.340 0.692) 3.066 ======================= Grid point 562 (76/126) ======================= q-point: ( 0.46 0.08 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.998 ( 0.528 -0.916 2.768) 2.963 2.170 ( -1.267 2.193 2.386) 3.479 2.436 ( 0.544 -0.943 -9.876) 9.936 2.947 ( 0.388 -0.674 7.124) 7.166 3.079 ( -5.571 9.647 2.861) 11.502 4.152 ( 0.901 -1.566 -1.134) 2.133 9.398 ( 0.269 -0.467 4.764) 4.794 10.692 ( -0.227 0.393 -0.548) 0.712 10.952 ( 1.535 -2.663 1.126) 3.273 ======================= Grid point 572 (77/126) ======================= q-point: ( 0.15 0.15 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.759 ( 2.529 4.381 7.770) 9.271 1.882 ( 5.906 10.230 9.056) 14.885 2.571 ( 0.802 1.389 6.696) 6.886 2.928 ( 0.112 0.194 3.128) 3.136 3.270 ( 0.671 1.162 6.279) 6.421 3.639 ( 3.274 5.670 -5.906) 8.817 9.262 ( 1.475 2.555 2.703) 4.001 10.740 ( 0.204 0.354 -0.379) 0.557 10.871 ( 0.131 0.227 -1.090) 1.121 ======================= Grid point 573 (78/126) ======================= q-point: ( 0.23 0.15 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.842 ( 2.925 3.036 4.462) 6.139 2.058 ( 3.732 8.205 5.296) 10.455 2.552 ( -1.485 -3.584 -2.831) 4.803 2.914 ( -0.137 -1.100 9.016) 9.084 3.294 ( -3.253 8.497 1.125) 9.168 3.811 ( 8.841 4.470 -1.831) 10.074 9.313 ( 1.841 1.879 3.367) 4.273 10.737 ( -1.352 1.074 -0.337) 1.759 10.875 ( 1.546 -2.198 -0.456) 2.726 ======================= Grid point 574 (79/126) ======================= q-point: ( 0.31 0.15 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.919 ( 3.036 2.028 2.815) 4.610 2.170 ( 1.218 5.904 2.343) 6.468 2.476 ( -2.001 -4.217 -7.629) 8.943 2.916 ( 0.982 -0.351 9.069) 9.129 3.315 ( -6.307 12.112 0.708) 13.675 3.989 ( 8.018 1.121 -0.366) 8.105 9.358 ( 1.666 0.832 3.789) 4.222 10.720 ( -1.595 0.793 -0.253) 1.799 10.880 ( 2.257 -3.443 0.512) 4.148 ======================= Grid point 575 (80/126) ======================= q-point: ( 0.38 0.15 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.975 ( 1.630 0.501 2.749) 3.235 2.220 ( -0.865 3.870 0.563) 4.005 2.419 ( -0.445 -2.428 -9.598) 9.910 2.930 ( 0.890 -0.524 8.129) 8.194 3.319 ( -7.750 13.376 1.246) 15.509 4.092 ( 4.115 -1.717 -0.203) 4.464 9.383 ( 0.902 -0.213 3.979) 4.086 10.705 ( -0.792 0.530 -0.219) 0.978 10.883 ( 2.333 -3.902 1.260) 4.718 ======================= Grid point 586 (81/126) ======================= q-point: ( 0.23 0.23 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.903 ( 1.781 3.086 2.793) 4.527 2.192 ( 3.031 5.249 1.554) 6.258 2.470 ( -2.357 -4.082 -6.975) 8.418 2.904 ( 0.062 0.108 9.664) 9.665 3.457 ( 4.045 7.006 -0.479) 8.104 3.902 ( 2.774 4.805 0.497) 5.571 9.347 ( 0.836 1.449 3.436) 3.822 10.749 ( -0.115 -0.199 -0.118) 0.259 10.830 ( -1.066 -1.847 0.434) 2.177 ======================= Grid point 587 (82/126) ======================= q-point: ( 0.31 0.23 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.958 ( 1.599 1.855 2.592) 3.566 2.265 ( 0.414 3.671 -1.367) 3.940 2.405 ( -1.070 -2.939 -8.748) 9.290 2.912 ( 0.625 -0.005 9.033) 9.054 3.553 ( -3.170 11.119 -0.820) 11.591 3.992 ( 4.787 -0.758 1.459) 5.062 9.369 ( 0.698 0.318 3.304) 3.392 10.730 ( -1.550 0.276 0.051) 1.575 10.810 ( 1.431 -3.234 1.327) 3.777 ======================= Grid point 588 (83/126) ======================= q-point: ( 0.38 0.23 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.983 ( -0.182 0.315 2.837) 2.860 2.285 ( -1.544 2.672 -2.625) 4.052 2.380 ( 0.916 -1.587 -9.202) 9.383 2.919 ( 0.280 -0.486 8.767) 8.785 3.579 ( -6.846 11.853 -0.772) 13.710 4.032 ( 2.417 -4.192 1.570) 5.087 9.376 ( 0.220 -0.382 3.221) 3.251 10.717 ( -0.431 0.745 0.138) 0.871 10.807 ( 1.958 -3.394 1.679) 4.263 ======================= Grid point 601 (84/126) ======================= q-point: ( 0.31 0.31 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.987 ( 0.473 0.820 2.843) 2.997 2.317 ( 0.980 1.697 -4.787) 5.172 2.360 ( -0.903 -1.564 -8.070) 8.270 2.914 ( 0.039 0.068 9.133) 9.134 3.724 ( 3.053 5.288 -3.562) 7.069 3.966 ( -0.818 -1.417 4.242) 4.546 9.372 ( 0.031 0.053 2.874) 2.874 10.731 ( -0.101 -0.175 -0.091) 0.221 10.765 ( -0.642 -1.112 2.219) 2.564 ======================= Grid point 724 (85/126) ======================= q-point: ( 0.00 0.00 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.620 ( 0.000 0.000 20.660) 20.660 1.620 ( -0.000 0.000 20.660) 20.660 2.604 ( 0.000 0.000 -10.468) 10.468 2.604 ( 0.000 0.000 -10.468) 10.468 2.977 ( 0.000 0.000 31.745) 31.745 4.173 ( 0.000 0.000 -3.470) 3.470 9.139 ( -0.000 -0.000 0.335) 0.335 10.758 ( 0.000 0.000 0.188) 0.188 10.758 ( 0.000 0.000 0.188) 0.188 ======================= Grid point 725 (86/126) ======================= q-point: ( 0.08 0.00 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.680 ( 4.927 2.845 19.277) 20.099 1.716 ( 6.920 3.995 18.882) 20.503 2.632 ( 2.460 1.420 -9.683) 10.091 2.712 ( 9.541 5.509 -8.274) 13.778 2.951 ( -2.202 -1.271 28.273) 28.387 4.064 ( -9.128 -5.270 -3.322) 11.051 9.157 ( 1.511 0.872 0.528) 1.823 10.751 ( -0.626 -0.361 0.152) 0.738 10.776 ( 1.556 0.899 0.077) 1.799 ======================= Grid point 726 (87/126) ======================= q-point: ( 0.15 0.00 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.820 ( 6.793 3.922 15.673) 17.526 1.862 ( 4.922 2.842 14.698) 15.759 2.713 ( 4.430 2.558 -7.322) 8.932 2.874 ( -4.523 -2.611 19.153) 19.853 3.024 ( 16.552 9.556 -3.070) 19.358 3.791 ( -13.528 -7.810 -2.597) 15.835 9.207 ( 2.814 1.625 1.073) 3.422 10.731 ( -0.985 -0.569 0.063) 1.139 10.824 ( 2.476 1.430 -0.110) 2.861 ======================= Grid point 727 (88/126) ======================= q-point: ( 0.23 0.00 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.934 ( 1.707 0.986 11.006) 11.181 1.969 ( 5.966 3.444 11.108) 13.071 2.739 ( -7.296 -4.213 5.318) 9.963 2.825 ( 5.087 2.937 -3.868) 7.033 3.344 ( 3.030 1.749 5.692) 6.681 3.567 ( 2.228 1.286 -3.640) 4.458 9.282 ( 3.659 2.112 1.842) 4.609 10.709 ( -0.946 -0.546 -0.043) 1.093 10.884 ( 2.621 1.513 -0.052) 3.027 ======================= Grid point 728 (89/126) ======================= q-point: ( 0.31 0.00 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.957 ( 0.653 0.377 7.106) 7.146 2.087 ( 4.246 2.451 7.018) 8.561 2.552 ( -8.535 -4.928 -5.361) 11.219 2.932 ( 4.097 2.365 -0.453) 4.752 3.217 ( -6.201 -3.580 6.146) 9.437 3.830 ( 12.471 7.200 0.434) 14.407 9.368 ( 3.711 2.142 2.646) 5.036 10.691 ( -0.607 -0.350 -0.149) 0.716 10.940 ( 2.241 1.294 0.445) 2.625 ======================= Grid point 729 (90/126) ======================= q-point: ( 0.38 0.00 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.979 ( 1.476 0.852 2.698) 3.191 2.164 ( 2.505 1.446 4.186) 5.088 2.375 ( -6.565 -3.791 -8.802) 11.617 3.006 ( 2.303 1.330 2.115) 3.398 3.116 ( -2.902 -1.675 4.960) 5.986 4.058 ( 7.306 4.218 1.646) 8.596 9.444 ( 2.766 1.597 3.319) 4.607 10.681 ( -0.238 -0.137 -0.251) 0.372 10.984 ( 1.542 0.890 1.194) 2.144 ======================= Grid point 730 (91/126) ======================= q-point: ( 0.46 0.00 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.017 ( 1.366 0.789 -0.126) 1.583 2.202 ( 0.830 0.479 2.759) 2.920 2.264 ( -2.831 -1.634 -9.069) 9.640 3.039 ( 0.685 0.395 3.501) 3.589 3.075 ( -0.841 -0.485 3.798) 3.920 4.167 ( 2.291 1.323 1.861) 3.235 9.488 ( 1.020 0.589 3.718) 3.900 10.678 ( -0.043 -0.025 -0.319) 0.323 11.008 ( 0.561 0.324 1.745) 1.861 ======================= Grid point 738 (92/126) ======================= q-point: ( 0.08 0.08 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.784 ( 4.024 6.969 16.647) 18.490 1.826 ( 3.315 5.742 15.954) 17.277 2.693 ( 2.643 4.577 -7.979) 9.571 2.899 ( -1.524 -2.639 22.866) 23.068 2.921 ( 7.930 13.735 -5.484) 16.782 3.873 ( -7.583 -13.134 -2.845) 15.431 9.190 ( 1.419 2.459 0.897) 2.977 10.741 ( -0.335 -0.580 0.065) 0.672 10.805 ( 1.117 1.934 -0.047) 2.234 ======================= Grid point 739 (93/126) ======================= q-point: ( 0.15 0.08 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.907 ( 2.389 2.212 12.203) 12.630 1.934 ( 4.372 6.201 12.489) 14.613 2.774 ( -1.937 -2.753 3.791) 5.070 2.825 ( 0.664 3.649 2.035) 4.230 3.259 ( 12.033 13.092 1.169) 17.820 3.598 ( -8.177 -10.348 -2.165) 13.365 9.251 ( 2.579 2.764 1.546) 4.084 10.727 ( -1.012 0.182 -0.060) 1.030 10.850 ( 2.226 1.029 -0.069) 2.453 ======================= Grid point 740 (94/126) ======================= q-point: ( 0.23 0.08 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.941 ( 1.031 -0.154 8.766) 8.828 2.061 ( 3.378 5.184 7.994) 10.109 2.618 ( -6.584 -7.574 -2.004) 10.234 2.936 ( 2.235 7.537 -0.366) 7.870 3.302 ( -7.433 -3.390 4.795) 9.473 3.676 ( 13.480 7.724 -0.405) 15.542 9.329 ( 3.275 2.466 2.284) 4.693 10.711 ( -1.386 0.745 -0.144) 1.580 10.899 ( 2.854 -0.096 0.242) 2.866 ======================= Grid point 741 (95/126) ======================= q-point: ( 0.31 0.08 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.957 ( 1.399 -0.177 4.544) 4.758 2.149 ( 1.975 3.473 4.560) 6.063 2.434 ( -6.123 -6.500 -7.775) 11.840 3.027 ( 0.562 5.420 3.052) 6.246 3.180 ( -5.238 1.313 3.400) 6.381 3.951 ( 10.277 4.896 1.694) 11.509 9.406 ( 3.206 1.542 2.928) 4.608 10.696 ( -1.188 0.813 -0.187) 1.451 10.943 ( 2.870 -1.031 0.904) 3.181 ======================= Grid point 742 (96/126) ======================= q-point: ( 0.38 0.08 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.990 ( 2.180 0.338 0.917) 2.389 2.199 ( 0.769 1.806 2.486) 3.168 2.295 ( -3.802 -3.957 -8.901) 10.457 3.060 ( -0.120 0.956 4.687) 4.785 3.142 ( -4.587 7.049 2.272) 8.712 4.117 ( 5.294 1.765 2.132) 5.973 9.463 ( 2.237 0.241 3.375) 4.056 10.687 ( -0.716 0.642 -0.216) 0.985 10.974 ( 2.356 -1.820 1.593) 3.377 ======================= Grid point 743 (97/126) ======================= q-point: ( 0.46 0.08 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.017 ( 0.418 -0.725 -0.359) 0.910 2.214 ( -0.384 0.664 1.738) 1.899 2.238 ( 0.493 -0.857 -8.905) 8.960 3.060 ( 0.018 -0.032 4.126) 4.126 3.141 ( -4.845 8.389 2.905) 10.114 4.170 ( 0.534 -0.931 2.167) 2.418 9.483 ( 0.580 -1.005 3.540) 3.725 10.684 ( -0.356 0.616 -0.230) 0.747 10.985 ( 1.444 -2.504 1.885) 3.451 ======================= Grid point 753 (98/126) ======================= q-point: ( 0.15 0.15 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.930 ( 0.246 0.426 9.285) 9.298 2.060 ( 3.387 5.867 8.086) 10.549 2.636 ( -5.032 -8.716 -0.556) 10.079 2.962 ( 4.105 7.110 0.856) 8.254 3.370 ( -3.251 -5.630 2.883) 7.112 3.574 ( 7.060 12.229 -0.388) 14.126 9.312 ( 1.868 3.235 2.129) 4.300 10.734 ( 0.161 0.278 -0.219) 0.389 10.865 ( 0.396 0.686 0.258) 0.833 ======================= Grid point 754 (99/126) ======================= q-point: ( 0.23 0.15 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.936 ( 0.655 -0.119 5.566) 5.606 2.156 ( 1.968 4.155 4.152) 6.195 2.459 ( -5.294 -8.007 -6.607) 11.653 3.058 ( 0.724 2.720 5.890) 6.528 3.300 ( -3.094 4.654 -0.610) 5.622 3.827 ( 9.595 7.251 2.550) 12.294 9.375 ( 2.262 2.088 2.593) 4.025 10.732 ( -1.271 1.124 -0.218) 1.711 10.882 ( 1.995 -1.401 0.841) 2.579 ======================= Grid point 755 (100/126) ======================= q-point: ( 0.31 0.15 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.959 ( 1.781 0.558 1.910) 2.671 2.209 ( 0.548 2.292 1.396) 2.739 2.313 ( -3.744 -5.329 -8.405) 10.633 3.064 ( -0.251 -0.397 5.968) 5.986 3.325 ( -4.395 11.160 0.153) 11.996 4.026 ( 7.430 2.764 3.086) 8.507 9.426 ( 2.130 0.512 2.856) 3.599 10.716 ( -1.607 1.037 -0.080) 1.915 10.904 ( 2.756 -2.697 1.500) 4.137 ======================= Grid point 756 (101/126) ======================= q-point: ( 0.38 0.15 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.997 ( 1.488 0.384 -0.068) 1.539 2.226 ( -0.213 0.810 -0.459) 0.955 2.230 ( -1.268 -2.393 -8.189) 8.625 3.058 ( 0.016 -0.288 4.697) 4.706 3.350 ( -6.294 12.125 1.651) 13.761 4.130 ( 3.301 -0.381 3.084) 4.533 9.455 ( 1.384 -0.901 2.966) 3.395 10.703 ( -0.957 0.885 0.020) 1.304 10.919 ( 2.495 -3.479 1.964) 4.710 ======================= Grid point 767 (102/126) ======================= q-point: ( 0.23 0.23 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.945 ( 0.704 1.219 2.272) 2.673 2.217 ( 1.104 1.913 0.748) 2.332 2.319 ( -3.389 -5.871 -8.094) 10.558 3.063 ( -0.364 -0.630 6.497) 6.537 3.437 ( 4.078 7.064 -1.300) 8.260 3.961 ( 3.571 6.185 4.145) 8.258 9.409 ( 0.785 1.360 2.618) 3.053 10.745 ( -0.023 -0.040 -0.156) 0.162 10.855 ( -0.547 -0.947 1.602) 1.940 ======================= Grid point 768 (103/126) ======================= q-point: ( 0.31 0.23 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.980 ( 1.403 1.464 0.152) 2.033 2.230 ( -1.491 -2.494 -6.633) 7.241 2.236 ( -0.401 -0.121 -3.188) 3.215 3.055 ( -0.047 -0.232 5.361) 5.367 3.542 ( -1.679 9.673 -0.147) 9.819 4.064 ( 3.847 1.258 4.596) 6.124 9.428 ( 0.814 -0.178 2.483) 2.619 10.732 ( -1.323 0.460 0.159) 1.409 10.848 ( 1.533 -2.522 2.059) 3.598 ======================= Grid point 769 (104/126) ======================= q-point: ( 0.38 0.23 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.004 ( -0.135 0.234 -0.450) 0.525 2.200 ( 0.320 -0.556 -7.938) 7.963 2.233 ( 0.083 -0.145 -2.704) 2.709 3.054 ( 0.034 -0.060 4.899) 4.899 3.574 ( -5.490 9.504 0.392) 10.983 4.105 ( 1.080 -1.874 4.602) 5.084 9.434 ( 0.554 -0.959 2.408) 2.650 10.721 ( -0.550 0.951 0.331) 1.147 10.850 ( 1.777 -3.081 2.199) 4.181 ======================= Grid point 782 (105/126) ======================= q-point: ( 0.31 0.31 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.006 ( 0.471 0.816 -0.603) 1.119 2.198 ( -0.299 -0.517 -7.189) 7.214 2.229 ( -0.332 -0.576 -4.202) 4.255 3.054 ( 0.017 0.029 5.074) 5.075 3.679 ( 2.147 3.719 -0.830) 4.374 4.079 ( 0.183 0.316 5.888) 5.900 9.424 ( -0.103 -0.178 2.155) 2.165 10.734 ( -0.066 -0.114 0.407) 0.427 10.815 ( -0.470 -0.814 2.340) 2.522 ======================= Grid point 905 (106/126) ======================= q-point: ( 0.00 0.00 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.997 ( 0.000 0.000 15.929) 15.929 1.997 ( -0.000 0.000 15.929) 15.929 2.370 ( 0.000 0.000 -12.200) 12.200 2.370 ( 0.000 0.000 -12.200) 12.200 3.546 ( 0.000 0.000 23.970) 23.970 3.989 ( 0.000 0.000 -14.145) 14.145 9.145 ( 0.000 0.000 0.168) 0.168 10.761 ( 0.000 0.000 0.099) 0.099 10.761 ( 0.000 0.000 0.099) 0.099 ======================= Grid point 906 (107/126) ======================= q-point: ( 0.08 0.00 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.026 ( 2.245 1.296 13.794) 14.036 2.044 ( 2.903 1.676 11.981) 12.441 2.422 ( 4.555 2.630 -10.235) 11.507 2.529 ( 13.857 8.000 -8.174) 17.967 3.457 ( -7.743 -4.471 20.404) 22.277 3.902 ( -7.064 -4.078 -11.590) 14.172 9.165 ( 1.738 1.004 0.234) 2.021 10.753 ( -0.666 -0.385 0.087) 0.774 10.778 ( 1.460 0.843 0.072) 1.688 ======================= Grid point 907 (108/126) ======================= q-point: ( 0.15 0.00 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.084 ( 2.586 1.493 8.944) 9.429 2.098 ( 1.755 1.013 7.262) 7.540 2.569 ( 7.919 4.572 -5.758) 10.806 2.929 ( 16.707 9.646 -3.003) 19.524 3.213 ( -10.686 -6.170 10.711) 16.340 3.702 ( -9.646 -5.569 -4.822) 12.137 9.222 ( 3.227 1.863 0.424) 3.750 10.733 ( -1.046 -0.604 0.057) 1.209 10.824 ( 2.439 1.408 0.039) 2.817 ======================= Grid point 908 (109/126) ======================= q-point: ( 0.23 0.00 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.123 ( 0.305 0.176 6.150) 6.160 2.139 ( 2.270 1.310 4.930) 5.584 2.760 ( 8.262 4.770 -2.098) 9.769 2.773 ( -16.911 -9.763 -0.588) 19.536 3.438 ( 9.025 5.211 3.839) 11.106 3.507 ( 0.946 0.546 -2.707) 2.919 9.308 ( 4.151 2.396 0.695) 4.843 10.709 ( -1.012 -0.584 0.019) 1.168 10.885 ( 2.824 1.630 0.102) 3.262 ======================= Grid point 909 (110/126) ======================= q-point: ( 0.31 0.00 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.109 ( -1.524 -0.880 6.797) 7.021 2.186 ( 1.835 1.059 2.647) 3.391 2.417 ( -13.089 -7.557 -6.496) 16.451 2.929 ( 6.265 3.617 -0.023) 7.235 3.298 ( -6.648 -3.838 1.997) 7.932 3.846 ( 13.786 7.959 0.653) 15.932 9.405 ( 4.169 2.407 0.977) 4.912 10.689 ( -0.672 -0.388 -0.021) 0.777 10.949 ( 2.684 1.550 0.318) 3.115 ======================= Grid point 910 (111/126) ======================= q-point: ( 0.38 0.00 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.060 ( -2.400 -1.386 5.557) 6.210 2.201 ( -5.984 -3.455 -8.216) 10.736 2.222 ( 1.219 0.704 1.505) 2.061 3.041 ( 3.547 2.048 1.059) 4.231 3.178 ( -3.828 -2.210 1.463) 4.656 4.090 ( 7.611 4.394 1.046) 8.850 9.490 ( 3.108 1.794 1.217) 3.789 10.678 ( -0.293 -0.169 -0.058) 0.344 11.004 ( 1.994 1.151 0.609) 2.381 ======================= Grid point 911 (112/126) ======================= q-point: ( 0.46 0.00 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.021 ( -0.816 -0.471 0.776) 1.220 2.119 ( -1.621 -0.936 -4.788) 5.141 2.241 ( 0.429 0.248 1.040) 1.152 3.093 ( 1.097 0.634 1.536) 1.991 3.121 ( -1.215 -0.701 1.061) 1.759 4.202 ( 2.324 1.342 1.129) 2.911 9.540 ( 1.152 0.665 1.363) 1.904 10.674 ( -0.066 -0.038 -0.083) 0.112 11.036 ( 0.748 0.432 0.811) 1.185 ======================= Grid point 919 (113/126) ======================= q-point: ( 0.08 0.08 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.068 ( 1.395 2.416 9.822) 10.211 2.087 ( 1.344 2.328 8.495) 8.910 2.531 ( 4.669 8.086 -6.860) 11.587 2.810 ( 10.709 18.549 -4.176) 21.822 3.285 ( -7.093 -12.286 13.666) 19.699 3.761 ( -5.438 -9.419 -6.727) 12.789 9.203 ( 1.634 2.830 0.364) 3.288 10.743 ( -0.387 -0.671 0.055) 0.776 10.806 ( 1.091 1.890 0.053) 2.183 ======================= Grid point 920 (114/126) ======================= q-point: ( 0.15 0.08 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.103 ( 1.023 0.288 5.900) 5.995 2.136 ( 1.402 2.592 6.319) 6.972 2.713 ( 5.158 8.265 -2.950) 10.179 2.916 ( -10.702 -11.129 1.562) 15.519 3.284 ( 14.361 14.290 1.434) 20.310 3.561 ( -6.961 -8.195 -1.348) 10.836 9.273 ( 2.935 3.167 0.595) 4.359 10.727 ( -1.027 0.045 0.005) 1.028 10.851 ( 2.274 1.190 0.085) 2.568 ======================= Grid point 921 (115/126) ======================= q-point: ( 0.23 0.08 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.107 ( -0.025 -1.440 6.049) 6.218 2.181 ( 0.978 2.239 3.646) 4.388 2.536 ( -12.066 -12.797 -4.470) 18.147 2.928 ( 3.379 11.150 -0.281) 11.655 3.365 ( -6.480 -5.776 1.542) 8.817 3.684 ( 14.887 10.464 0.556) 18.205 9.361 ( 3.673 2.768 0.856) 4.678 10.710 ( -1.376 0.646 -0.027) 1.521 10.905 ( 3.107 0.255 0.236) 3.126 ======================= Grid point 922 (116/126) ======================= q-point: ( 0.31 0.08 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.073 ( -1.454 -2.518 6.736) 7.337 2.212 ( 0.573 0.861 1.465) 1.794 2.265 ( -7.086 -6.807 -8.129) 12.753 3.071 ( 0.757 9.977 1.275) 10.087 3.225 ( -4.558 -3.101 1.030) 5.608 3.984 ( 10.389 5.863 1.107) 11.981 9.447 ( 3.597 1.639 1.082) 4.098 10.694 ( -1.211 0.817 -0.040) 1.461 10.959 ( 3.264 -0.630 0.499) 3.361 ======================= Grid point 923 (117/126) ======================= q-point: ( 0.38 0.08 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.025 ( -1.018 -1.866 2.819) 3.530 2.136 ( -2.459 -2.815 -6.308) 7.333 2.234 ( 0.551 0.202 0.848) 1.031 3.132 ( -1.516 0.715 2.449) 2.967 3.166 ( -2.898 7.049 0.065) 7.622 4.156 ( 5.087 2.298 1.241) 5.719 9.510 ( 2.574 0.105 1.239) 2.858 10.685 ( -0.787 0.741 -0.045) 1.082 11.000 ( 2.637 -1.535 0.754) 3.143 ======================= Grid point 924 (118/126) ======================= q-point: ( 0.46 0.08 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.008 ( 0.407 -0.706 -0.157) 0.829 2.102 ( 0.400 -0.697 -3.969) 4.049 2.239 ( 0.262 -0.455 0.648) 0.834 3.114 ( -0.045 0.077 1.281) 1.284 3.185 ( -4.610 7.982 1.259) 9.304 4.210 ( 0.277 -0.484 1.255) 1.374 9.533 ( 0.760 -1.316 1.298) 1.999 10.682 ( -0.438 0.759 -0.047) 0.877 11.015 ( 1.400 -2.428 0.858) 2.931 ======================= Grid point 934 (119/126) ======================= q-point: ( 0.15 0.15 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.098 ( -0.491 -0.851 5.731) 5.815 2.185 ( 1.241 2.149 3.991) 4.700 2.575 ( -8.956 -15.512 -3.723) 18.294 2.970 ( 6.303 10.916 0.073) 12.605 3.396 ( -4.504 -7.801 0.388) 9.016 3.594 ( 8.942 15.489 1.101) 17.918 9.342 ( 2.094 3.626 0.809) 4.265 10.731 ( 0.111 0.192 -0.070) 0.232 10.872 ( 0.613 1.062 0.246) 1.250 ======================= Grid point 935 (120/126) ======================= q-point: ( 0.23 0.15 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.068 ( -1.094 -2.445 6.775) 7.285 2.215 ( 0.270 0.847 1.737) 1.951 2.304 ( -6.790 -9.702 -7.805) 14.183 3.158 ( 2.022 9.219 3.754) 10.157 3.271 ( -3.147 -1.521 -1.917) 3.987 3.875 ( 9.454 8.721 1.527) 12.952 9.411 ( 2.477 2.206 0.973) 3.457 10.729 ( -1.163 1.101 -0.070) 1.603 10.897 ( 2.229 -0.866 0.482) 2.440 ======================= Grid point 936 (121/126) ======================= q-point: ( 0.31 0.15 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.021 ( -1.072 -2.378 4.579) 5.270 2.147 ( -2.820 -4.068 -7.423) 8.922 2.229 ( -0.025 -0.282 0.285) 0.402 3.154 ( -1.920 -1.753 2.723) 3.765 3.320 ( -2.070 12.119 -0.465) 12.303 4.079 ( 6.899 3.800 1.632) 8.043 9.466 ( 2.371 0.323 1.057) 2.616 10.716 ( -1.554 1.187 -0.000) 1.956 10.928 ( 3.003 -2.274 0.720) 3.835 ======================= Grid point 937 (122/126) ======================= q-point: ( 0.38 0.15 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.996 ( 0.003 -0.891 0.340) 0.953 2.094 ( -0.625 -1.423 -4.268) 4.542 2.224 ( 0.446 -1.170 -0.152) 1.261 3.120 ( -0.757 -0.364 1.568) 1.779 3.377 ( -5.157 11.732 0.810) 12.841 4.183 ( 2.750 0.503 1.619) 3.230 9.496 ( 1.654 -1.300 1.091) 2.369 10.704 ( -1.044 1.128 0.053) 1.537 10.949 ( 2.583 -3.278 0.878) 4.265 ======================= Grid point 948 (123/126) ======================= q-point: ( 0.23 0.23 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.018 ( -1.301 -2.253 5.193) 5.808 2.151 ( -2.679 -4.641 -7.765) 9.435 2.228 ( -0.368 -0.638 -0.007) 0.737 3.162 ( -1.405 -2.434 2.953) 4.077 3.412 ( 5.150 8.920 -0.984) 10.346 4.029 ( 3.920 6.790 2.010) 8.094 9.446 ( 0.742 1.285 0.976) 1.776 10.743 ( 0.078 0.134 -0.029) 0.158 10.881 ( -0.285 -0.494 0.753) 0.945 ======================= Grid point 949 (124/126) ======================= q-point: ( 0.31 0.23 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 930 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.990 ( -0.217 -0.670 1.163) 1.360 2.094 ( -0.918 -1.414 -4.870) 5.154 2.207 ( -0.394 -1.782 -0.834) 2.006 3.129 ( -0.729 -0.739 1.895) 2.161 3.544 ( -0.507 9.606 0.174) 9.621 4.138 ( 3.342 2.249 2.144) 4.564 9.463 ( 0.867 -0.475 0.918) 1.349 10.735 ( -1.165 0.691 0.118) 1.360 10.880 ( 1.570 -2.239 0.896) 2.878 ======================= Grid point 950 (125/126) ======================= q-point: ( 0.38 0.23 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.987 ( 0.078 -0.136 -0.664) 0.683 2.081 ( -0.055 0.093 -3.337) 3.339 2.192 ( 1.041 -1.806 -1.170) 2.390 3.119 ( -0.081 0.139 1.617) 1.625 3.587 ( -4.917 8.512 0.504) 9.844 4.179 ( 0.408 -0.711 2.148) 2.299 9.468 ( 0.745 -1.290 0.887) 1.734 10.727 ( -0.668 1.155 0.194) 1.348 10.883 ( 1.726 -2.992 0.937) 3.579 ======================= Grid point 963 (126/126) ======================= q-point: ( 0.31 0.31 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.46e-03 0.00e+00 6.24e-05 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.986 ( 0.020 0.035 -0.737) 0.738 2.086 ( 0.186 0.322 -3.261) 3.282 2.169 ( -0.894 -1.549 -1.587) 2.391 3.121 ( -0.063 -0.109 1.681) 1.686 3.678 ( 2.036 3.526 0.211) 4.077 4.169 ( 0.460 0.796 2.500) 2.664 9.454 ( -0.181 -0.314 0.796) 0.874 10.742 ( 0.032 0.055 0.256) 0.264 10.849 ( -0.450 -0.779 0.953) 1.310 =================== 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/16731 10.0 711.432 711.432 465.950 0.000 0.000 0.000 3/16731 20.0 82.945 82.945 76.669 0.000 0.000 -0.000 3/16731 30.0 36.588 36.588 39.702 0.000 0.000 -0.000 3/16731 40.0 24.536 24.536 28.518 -0.000 0.000 -0.000 3/16731 50.0 19.100 19.100 22.764 -0.000 0.000 -0.000 3/16731 60.0 15.878 15.878 19.087 -0.000 0.000 -0.000 3/16731 70.0 13.670 13.670 16.475 -0.000 0.000 -0.000 3/16731 80.0 12.031 12.031 14.502 -0.000 0.000 -0.000 3/16731 90.0 10.752 10.752 12.954 -0.000 0.000 -0.000 3/16731 100.0 9.722 9.722 11.703 -0.000 0.000 -0.000 3/16731 110.0 8.873 8.873 10.672 -0.000 0.000 -0.000 3/16731 120.0 8.160 8.160 9.807 -0.000 0.000 -0.000 3/16731 130.0 7.553 7.553 9.071 -0.000 0.000 -0.000 3/16731 140.0 7.029 7.029 8.437 -0.000 0.000 -0.000 3/16731 150.0 6.573 6.573 7.886 -0.000 0.000 -0.000 3/16731 160.0 6.172 6.172 7.402 -0.000 0.000 -0.000 3/16731 170.0 5.818 5.818 6.974 -0.000 0.000 -0.000 3/16731 180.0 5.501 5.501 6.593 -0.000 0.000 -0.000 3/16731 190.0 5.217 5.217 6.251 -0.000 0.000 -0.000 3/16731 200.0 4.961 4.961 5.943 -0.000 0.000 -0.000 3/16731 210.0 4.729 4.729 5.663 -0.000 0.000 -0.000 3/16731 220.0 4.518 4.518 5.409 -0.000 0.000 -0.000 3/16731 230.0 4.324 4.324 5.177 -0.000 0.000 -0.000 3/16731 240.0 4.147 4.147 4.964 -0.000 0.000 -0.000 3/16731 250.0 3.983 3.983 4.768 -0.000 0.000 -0.000 3/16731 260.0 3.832 3.832 4.586 -0.000 0.000 -0.000 3/16731 270.0 3.692 3.692 4.418 -0.000 0.000 -0.000 3/16731 280.0 3.562 3.562 4.262 -0.000 0.000 -0.000 3/16731 290.0 3.440 3.440 4.117 -0.000 0.000 -0.000 3/16731 300.0 3.327 3.327 3.981 -0.000 0.000 -0.000 3/16731 310.0 3.221 3.221 3.854 -0.000 0.000 -0.000 3/16731 320.0 3.121 3.121 3.735 -0.000 0.000 -0.000 3/16731 330.0 3.028 3.028 3.622 -0.000 0.000 -0.000 3/16731 340.0 2.939 2.939 3.517 -0.000 0.000 -0.000 3/16731 350.0 2.856 2.856 3.417 -0.000 0.000 -0.000 3/16731 360.0 2.778 2.778 3.323 -0.000 0.000 -0.000 3/16731 370.0 2.703 2.703 3.234 -0.000 0.000 -0.000 3/16731 380.0 2.633 2.633 3.150 -0.000 0.000 -0.000 3/16731 390.0 2.566 2.566 3.070 -0.000 0.000 -0.000 3/16731 400.0 2.502 2.502 2.994 -0.000 0.000 -0.000 3/16731 410.0 2.442 2.442 2.921 -0.000 0.000 -0.000 3/16731 420.0 2.384 2.384 2.852 -0.000 0.000 -0.000 3/16731 430.0 2.329 2.329 2.786 -0.000 0.000 -0.000 3/16731 440.0 2.277 2.277 2.724 -0.000 0.000 -0.000 3/16731 450.0 2.226 2.226 2.664 -0.000 0.000 -0.000 3/16731 460.0 2.178 2.178 2.606 -0.000 0.000 -0.000 3/16731 470.0 2.132 2.132 2.551 -0.000 0.000 -0.000 3/16731 480.0 2.088 2.088 2.498 -0.000 0.000 -0.000 3/16731 490.0 2.046 2.046 2.448 -0.000 0.000 -0.000 3/16731 500.0 2.005 2.005 2.399 -0.000 0.000 -0.000 3/16731 510.0 1.966 1.966 2.352 -0.000 0.000 -0.000 3/16731 520.0 1.929 1.929 2.307 -0.000 0.000 -0.000 3/16731 530.0 1.893 1.893 2.264 -0.000 0.000 -0.000 3/16731 540.0 1.858 1.858 2.223 -0.000 0.000 -0.000 3/16731 550.0 1.824 1.824 2.182 -0.000 0.000 -0.000 3/16731 560.0 1.792 1.792 2.144 -0.000 0.000 -0.000 3/16731 570.0 1.761 1.761 2.106 -0.000 0.000 -0.000 3/16731 580.0 1.730 1.730 2.070 -0.000 0.000 -0.000 3/16731 590.0 1.701 1.701 2.035 -0.000 0.000 -0.000 3/16731 600.0 1.673 1.673 2.002 -0.000 0.000 -0.000 3/16731 610.0 1.646 1.646 1.969 -0.000 0.000 -0.000 3/16731 620.0 1.619 1.619 1.938 -0.000 0.000 -0.000 3/16731 630.0 1.594 1.594 1.907 -0.000 0.000 -0.000 3/16731 640.0 1.569 1.569 1.877 -0.000 0.000 -0.000 3/16731 650.0 1.545 1.545 1.849 -0.000 0.000 -0.000 3/16731 660.0 1.522 1.522 1.821 -0.000 0.000 -0.000 3/16731 670.0 1.499 1.499 1.794 -0.000 0.000 -0.000 3/16731 680.0 1.477 1.477 1.768 -0.000 0.000 -0.000 3/16731 690.0 1.456 1.456 1.742 -0.000 0.000 -0.000 3/16731 700.0 1.435 1.435 1.717 -0.000 0.000 -0.000 3/16731 710.0 1.415 1.415 1.693 -0.000 0.000 -0.000 3/16731 720.0 1.396 1.396 1.670 -0.000 0.000 -0.000 3/16731 730.0 1.377 1.377 1.647 -0.000 0.000 -0.000 3/16731 740.0 1.358 1.358 1.625 -0.000 0.000 -0.000 3/16731 750.0 1.340 1.340 1.604 -0.000 0.000 -0.000 3/16731 760.0 1.323 1.323 1.583 -0.000 0.000 -0.000 3/16731 770.0 1.306 1.306 1.562 -0.000 0.000 -0.000 3/16731 780.0 1.289 1.289 1.542 -0.000 0.000 -0.000 3/16731 790.0 1.273 1.273 1.523 -0.000 0.000 -0.000 3/16731 800.0 1.257 1.257 1.504 -0.000 0.000 -0.000 3/16731 810.0 1.241 1.241 1.485 -0.000 0.000 -0.000 3/16731 820.0 1.226 1.226 1.467 -0.000 0.000 -0.000 3/16731 830.0 1.212 1.212 1.450 -0.000 0.000 -0.000 3/16731 840.0 1.197 1.197 1.433 -0.000 0.000 -0.000 3/16731 850.0 1.183 1.183 1.416 -0.000 0.000 -0.000 3/16731 860.0 1.170 1.170 1.400 -0.000 0.000 -0.000 3/16731 870.0 1.156 1.156 1.384 -0.000 0.000 -0.000 3/16731 880.0 1.143 1.143 1.368 -0.000 0.000 -0.000 3/16731 890.0 1.130 1.130 1.353 -0.000 0.000 -0.000 3/16731 900.0 1.118 1.118 1.338 -0.000 0.000 -0.000 3/16731 910.0 1.106 1.106 1.323 -0.000 0.000 -0.000 3/16731 920.0 1.094 1.094 1.309 -0.000 0.000 -0.000 3/16731 930.0 1.082 1.082 1.295 -0.000 0.000 -0.000 3/16731 940.0 1.070 1.070 1.281 -0.000 0.000 -0.000 3/16731 950.0 1.059 1.059 1.268 -0.000 0.000 -0.000 3/16731 960.0 1.048 1.048 1.254 -0.000 0.000 -0.000 3/16731 970.0 1.037 1.037 1.242 -0.000 0.000 -0.000 3/16731 980.0 1.027 1.027 1.229 -0.000 0.000 -0.000 3/16731 990.0 1.017 1.017 1.217 -0.000 0.000 -0.000 3/16731 1000.0 1.006 1.006 1.205 -0.000 0.000 -0.000 3/16731 Thermal conductivity related properties were written into "kappa-m131311.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 23:45:42]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|