----------------------------- ------- 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.378137629250306 0.000000000000000 0.000000000000000 b -2.189068814625153 3.791578408195341 0.000000000000000 c 0.000000000000000 0.000000000000000 4.456480789999999 Atomic positions (fractional): *1 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 *2 P 0.33333333333333 0.66666666666667 0.50000000000000 30.974 *3 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.378137629250306 0.000000000000000 0.000000000000000 b -2.189068814625153 3.791578408195341 0.000000000000000 c 0.000000000000000 0.000000000000000 4.456480789999999 Atomic positions (fractional): *1 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 > 1 *2 P 0.33333333333333 0.66666666666667 0.50000000000000 30.974 > 2 *3 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 > 3 -------------------------------- super cell -------------------------------- Lattice vectors: a 17.512550517001223 0.000000000000000 0.000000000000000 b -8.756275258500612 15.166313632781364 0.000000000000000 c 0.000000000000000 0.000000000000000 13.369442369999996 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 P 0.08333333333333 0.16666666666667 0.16666666666667 30.974 > 2 50 P 0.33333333333333 0.16666666666667 0.16666666666667 30.974 > 2 51 P 0.58333333333333 0.16666666666667 0.16666666666667 30.974 > 2 52 P 0.83333333333333 0.16666666666667 0.16666666666667 30.974 > 2 53 P 0.08333333333333 0.41666666666667 0.16666666666667 30.974 > 2 54 P 0.33333333333333 0.41666666666667 0.16666666666667 30.974 > 2 55 P 0.58333333333333 0.41666666666667 0.16666666666667 30.974 > 2 56 P 0.83333333333333 0.41666666666667 0.16666666666667 30.974 > 2 57 P 0.08333333333333 0.66666666666667 0.16666666666667 30.974 > 2 58 P 0.33333333333333 0.66666666666667 0.16666666666667 30.974 > 2 59 P 0.58333333333333 0.66666666666667 0.16666666666667 30.974 > 2 60 P 0.83333333333333 0.66666666666667 0.16666666666667 30.974 > 2 61 P 0.08333333333333 0.91666666666667 0.16666666666667 30.974 > 2 62 P 0.33333333333333 0.91666666666667 0.16666666666667 30.974 > 2 63 P 0.58333333333333 0.91666666666667 0.16666666666667 30.974 > 2 64 P 0.83333333333333 0.91666666666667 0.16666666666667 30.974 > 2 65 P 0.08333333333333 0.16666666666667 0.50000000000000 30.974 > 2 66 P 0.33333333333333 0.16666666666667 0.50000000000000 30.974 > 2 67 P 0.58333333333333 0.16666666666667 0.50000000000000 30.974 > 2 68 P 0.83333333333333 0.16666666666667 0.50000000000000 30.974 > 2 69 P 0.08333333333333 0.41666666666667 0.50000000000000 30.974 > 2 70 P 0.33333333333333 0.41666666666667 0.50000000000000 30.974 > 2 71 P 0.58333333333333 0.41666666666667 0.50000000000000 30.974 > 2 72 P 0.83333333333333 0.41666666666667 0.50000000000000 30.974 > 2 73 P 0.08333333333333 0.66666666666667 0.50000000000000 30.974 > 2 74 P 0.33333333333333 0.66666666666667 0.50000000000000 30.974 > 2 75 P 0.58333333333333 0.66666666666667 0.50000000000000 30.974 > 2 76 P 0.83333333333333 0.66666666666667 0.50000000000000 30.974 > 2 77 P 0.08333333333333 0.91666666666667 0.50000000000000 30.974 > 2 78 P 0.33333333333333 0.91666666666667 0.50000000000000 30.974 > 2 79 P 0.58333333333333 0.91666666666667 0.50000000000000 30.974 > 2 80 P 0.83333333333333 0.91666666666667 0.50000000000000 30.974 > 2 81 P 0.08333333333333 0.16666666666667 0.83333333333333 30.974 > 2 82 P 0.33333333333333 0.16666666666667 0.83333333333333 30.974 > 2 83 P 0.58333333333333 0.16666666666667 0.83333333333333 30.974 > 2 84 P 0.83333333333333 0.16666666666667 0.83333333333333 30.974 > 2 85 P 0.08333333333333 0.41666666666667 0.83333333333333 30.974 > 2 86 P 0.33333333333333 0.41666666666667 0.83333333333333 30.974 > 2 87 P 0.58333333333333 0.41666666666667 0.83333333333333 30.974 > 2 88 P 0.83333333333333 0.41666666666667 0.83333333333333 30.974 > 2 89 P 0.08333333333333 0.66666666666667 0.83333333333333 30.974 > 2 90 P 0.33333333333333 0.66666666666667 0.83333333333333 30.974 > 2 91 P 0.58333333333333 0.66666666666667 0.83333333333333 30.974 > 2 92 P 0.83333333333333 0.66666666666667 0.83333333333333 30.974 > 2 93 P 0.08333333333333 0.91666666666667 0.83333333333333 30.974 > 2 94 P 0.33333333333333 0.91666666666667 0.83333333333333 30.974 > 2 95 P 0.58333333333333 0.91666666666667 0.83333333333333 30.974 > 2 96 P 0.83333333333333 0.91666666666667 0.83333333333333 30.974 > 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 ---------------------------- 8.9317197 0.0000000 0.0000000 0.0000000 8.9317197 0.0000000 0.0000000 0.0000000 30.7747607 -------------------------- Born effective charges -------------------------- 1 Li 1.0053597 0.0000000 0.0000000 0.0000000 1.0053597 0.0000000 0.0000000 0.0000000 0.2180664 2 P -3.2075593 -0.0000000 0.0000000 0.0000000 -3.2075593 0.0000000 0.0000000 0.0000000 -5.5391043 3 Ba 2.2021996 0.0000000 0.0000000 -0.0000000 2.2021996 0.0000000 0.0000000 0.0000000 5.3210379 ---------------------------------------------------------------------------- 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: 174 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.014 Solver_block: 80 / 80 - Time: 0.168 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.183 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (xx) -0.00000000 (xx) 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: 174 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 (xx) -0.00000000 (xx) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 23:44:38]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ 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.378137629250306 0.000000000000000 0.000000000000000 b -2.189068814625153 3.791578408195341 0.000000000000000 c 0.000000000000000 0.000000000000000 4.456480789999999 Atomic positions (fractional): 1 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 2 P 0.33333333333333 0.66666666666667 0.50000000000000 30.974 3 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.512550517001223 0.000000000000000 0.000000000000000 b -8.756275258500612 15.166313632781364 0.000000000000000 c 0.000000000000000 0.000000000000000 13.369442369999996 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 P 0.08333333333333 0.16666666666667 0.16666666666667 30.974 > 49 50 P 0.33333333333333 0.16666666666667 0.16666666666667 30.974 > 49 51 P 0.58333333333333 0.16666666666667 0.16666666666667 30.974 > 49 52 P 0.83333333333333 0.16666666666667 0.16666666666667 30.974 > 49 53 P 0.08333333333333 0.41666666666667 0.16666666666667 30.974 > 49 54 P 0.33333333333333 0.41666666666667 0.16666666666667 30.974 > 49 55 P 0.58333333333333 0.41666666666667 0.16666666666667 30.974 > 49 56 P 0.83333333333333 0.41666666666667 0.16666666666667 30.974 > 49 57 P 0.08333333333333 0.66666666666667 0.16666666666667 30.974 > 49 58 P 0.33333333333333 0.66666666666667 0.16666666666667 30.974 > 49 59 P 0.58333333333333 0.66666666666667 0.16666666666667 30.974 > 49 60 P 0.83333333333333 0.66666666666667 0.16666666666667 30.974 > 49 61 P 0.08333333333333 0.91666666666667 0.16666666666667 30.974 > 49 62 P 0.33333333333333 0.91666666666667 0.16666666666667 30.974 > 49 63 P 0.58333333333333 0.91666666666667 0.16666666666667 30.974 > 49 64 P 0.83333333333333 0.91666666666667 0.16666666666667 30.974 > 49 65 P 0.08333333333333 0.16666666666667 0.50000000000000 30.974 > 49 66 P 0.33333333333333 0.16666666666667 0.50000000000000 30.974 > 49 67 P 0.58333333333333 0.16666666666667 0.50000000000000 30.974 > 49 68 P 0.83333333333333 0.16666666666667 0.50000000000000 30.974 > 49 69 P 0.08333333333333 0.41666666666667 0.50000000000000 30.974 > 49 70 P 0.33333333333333 0.41666666666667 0.50000000000000 30.974 > 49 71 P 0.58333333333333 0.41666666666667 0.50000000000000 30.974 > 49 72 P 0.83333333333333 0.41666666666667 0.50000000000000 30.974 > 49 73 P 0.08333333333333 0.66666666666667 0.50000000000000 30.974 > 49 74 P 0.33333333333333 0.66666666666667 0.50000000000000 30.974 > 49 75 P 0.58333333333333 0.66666666666667 0.50000000000000 30.974 > 49 76 P 0.83333333333333 0.66666666666667 0.50000000000000 30.974 > 49 77 P 0.08333333333333 0.91666666666667 0.50000000000000 30.974 > 49 78 P 0.33333333333333 0.91666666666667 0.50000000000000 30.974 > 49 79 P 0.58333333333333 0.91666666666667 0.50000000000000 30.974 > 49 80 P 0.83333333333333 0.91666666666667 0.50000000000000 30.974 > 49 81 P 0.08333333333333 0.16666666666667 0.83333333333333 30.974 > 49 82 P 0.33333333333333 0.16666666666667 0.83333333333333 30.974 > 49 83 P 0.58333333333333 0.16666666666667 0.83333333333333 30.974 > 49 84 P 0.83333333333333 0.16666666666667 0.83333333333333 30.974 > 49 85 P 0.08333333333333 0.41666666666667 0.83333333333333 30.974 > 49 86 P 0.33333333333333 0.41666666666667 0.83333333333333 30.974 > 49 87 P 0.58333333333333 0.41666666666667 0.83333333333333 30.974 > 49 88 P 0.83333333333333 0.41666666666667 0.83333333333333 30.974 > 49 89 P 0.08333333333333 0.66666666666667 0.83333333333333 30.974 > 49 90 P 0.33333333333333 0.66666666666667 0.83333333333333 30.974 > 49 91 P 0.58333333333333 0.66666666666667 0.83333333333333 30.974 > 49 92 P 0.83333333333333 0.66666666666667 0.83333333333333 30.974 > 49 93 P 0.08333333333333 0.91666666666667 0.83333333333333 30.974 > 49 94 P 0.33333333333333 0.91666666666667 0.83333333333333 30.974 > 49 95 P 0.58333333333333 0.91666666666667 0.83333333333333 30.974 > 49 96 P 0.83333333333333 0.91666666666667 0.83333333333333 30.974 > 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 ---------------------------- 8.9317197 0.0000000 0.0000000 0.0000000 8.9317197 0.0000000 0.0000000 0.0000000 30.7747607 -------------------------- Born effective charges -------------------------- 1 Li 1.0053597 0.0000000 0.0000000 0.0000000 1.0053597 0.0000000 0.0000000 0.0000000 0.2180664 2 P -3.2075593 -0.0000000 0.0000000 0.0000000 -3.2075593 0.0000000 0.0000000 0.0000000 -5.5391043 3 Ba 2.2021996 0.0000000 0.0000000 -0.0000000 2.2021996 0.0000000 0.0000000 0.0000000 5.3210379 ---------------------------------------------------------------------------- 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.00000013 (xyx) -0.00000013 (xyx) -0.00000013 (xxy) 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:45]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 3] Primitive matrix: [ 1. -0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P-6m2 (187) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.378137629250306 0.000000000000000 0.000000000000000 b -2.189068814625153 3.791578408195341 0.000000000000000 c 0.000000000000000 0.000000000000000 4.456480789999999 Atomic positions (fractional): 1 Li 0.00000000000000 0.00000000000000 0.50000000000000 6.941 2 P 0.33333333333333 0.66666666666667 0.50000000000000 30.974 3 Ba 0.66666666666667 0.33333333333333 0.00000000000000 137.327 -------------------------------- supercell --------------------------------- Lattice vectors: a 17.512550517001223 0.000000000000000 0.000000000000000 b -8.756275258500612 15.166313632781364 0.000000000000000 c 0.000000000000000 0.000000000000000 13.369442369999996 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 P 0.08333333333333 0.16666666666667 0.16666666666667 30.974 > 49 50 P 0.33333333333333 0.16666666666667 0.16666666666667 30.974 > 49 51 P 0.58333333333333 0.16666666666667 0.16666666666667 30.974 > 49 52 P 0.83333333333333 0.16666666666667 0.16666666666667 30.974 > 49 53 P 0.08333333333333 0.41666666666667 0.16666666666667 30.974 > 49 54 P 0.33333333333333 0.41666666666667 0.16666666666667 30.974 > 49 55 P 0.58333333333333 0.41666666666667 0.16666666666667 30.974 > 49 56 P 0.83333333333333 0.41666666666667 0.16666666666667 30.974 > 49 57 P 0.08333333333333 0.66666666666667 0.16666666666667 30.974 > 49 58 P 0.33333333333333 0.66666666666667 0.16666666666667 30.974 > 49 59 P 0.58333333333333 0.66666666666667 0.16666666666667 30.974 > 49 60 P 0.83333333333333 0.66666666666667 0.16666666666667 30.974 > 49 61 P 0.08333333333333 0.91666666666667 0.16666666666667 30.974 > 49 62 P 0.33333333333333 0.91666666666667 0.16666666666667 30.974 > 49 63 P 0.58333333333333 0.91666666666667 0.16666666666667 30.974 > 49 64 P 0.83333333333333 0.91666666666667 0.16666666666667 30.974 > 49 65 P 0.08333333333333 0.16666666666667 0.50000000000000 30.974 > 49 66 P 0.33333333333333 0.16666666666667 0.50000000000000 30.974 > 49 67 P 0.58333333333333 0.16666666666667 0.50000000000000 30.974 > 49 68 P 0.83333333333333 0.16666666666667 0.50000000000000 30.974 > 49 69 P 0.08333333333333 0.41666666666667 0.50000000000000 30.974 > 49 70 P 0.33333333333333 0.41666666666667 0.50000000000000 30.974 > 49 71 P 0.58333333333333 0.41666666666667 0.50000000000000 30.974 > 49 72 P 0.83333333333333 0.41666666666667 0.50000000000000 30.974 > 49 73 P 0.08333333333333 0.66666666666667 0.50000000000000 30.974 > 49 74 P 0.33333333333333 0.66666666666667 0.50000000000000 30.974 > 49 75 P 0.58333333333333 0.66666666666667 0.50000000000000 30.974 > 49 76 P 0.83333333333333 0.66666666666667 0.50000000000000 30.974 > 49 77 P 0.08333333333333 0.91666666666667 0.50000000000000 30.974 > 49 78 P 0.33333333333333 0.91666666666667 0.50000000000000 30.974 > 49 79 P 0.58333333333333 0.91666666666667 0.50000000000000 30.974 > 49 80 P 0.83333333333333 0.91666666666667 0.50000000000000 30.974 > 49 81 P 0.08333333333333 0.16666666666667 0.83333333333333 30.974 > 49 82 P 0.33333333333333 0.16666666666667 0.83333333333333 30.974 > 49 83 P 0.58333333333333 0.16666666666667 0.83333333333333 30.974 > 49 84 P 0.83333333333333 0.16666666666667 0.83333333333333 30.974 > 49 85 P 0.08333333333333 0.41666666666667 0.83333333333333 30.974 > 49 86 P 0.33333333333333 0.41666666666667 0.83333333333333 30.974 > 49 87 P 0.58333333333333 0.41666666666667 0.83333333333333 30.974 > 49 88 P 0.83333333333333 0.41666666666667 0.83333333333333 30.974 > 49 89 P 0.08333333333333 0.66666666666667 0.83333333333333 30.974 > 49 90 P 0.33333333333333 0.66666666666667 0.83333333333333 30.974 > 49 91 P 0.58333333333333 0.66666666666667 0.83333333333333 30.974 > 49 92 P 0.83333333333333 0.66666666666667 0.83333333333333 30.974 > 49 93 P 0.08333333333333 0.91666666666667 0.83333333333333 30.974 > 49 94 P 0.33333333333333 0.91666666666667 0.83333333333333 30.974 > 49 95 P 0.58333333333333 0.91666666666667 0.83333333333333 30.974 > 49 96 P 0.83333333333333 0.91666666666667 0.83333333333333 30.974 > 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 ---------------------------- 8.9317197 0.0000000 0.0000000 0.0000000 8.9317197 0.0000000 0.0000000 0.0000000 30.7747607 -------------------------- Born effective charges -------------------------- 1 Li 1.0053597 0.0000000 0.0000000 0.0000000 1.0053597 0.0000000 0.0000000 0.0000000 0.2180664 2 P -3.2075593 -0.0000000 0.0000000 0.0000000 -3.2075593 0.0000000 0.0000000 0.0000000 -5.5391043 3 Ba 2.2021996 0.0000000 0.0000000 -0.0000000 2.2021996 0.0000000 0.0000000 0.0000000 5.3210379 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000013 (xyx) -0.00000013 (xyx) -0.00000013 (xxy) Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 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.99, Number of G-points: 315, Lambda: 0.42 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 2.603 ( -0.000 0.000 0.000) 0.000 4.183 ( 0.000 -0.000 0.000) 0.000 4.183 ( 0.000 -0.000 0.000) 0.000 9.798 ( -0.000 0.000 0.000) 0.000 11.419 ( 0.000 -0.000 0.000) 0.000 11.419 ( 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.445 ( 18.420 10.635 0.000) 21.269 0.551 ( 22.568 13.029 0.000) 26.059 0.928 ( 37.993 21.935 0.000) 43.871 2.662 ( 4.862 2.807 0.000) 5.615 4.214 ( 2.590 1.495 0.000) 2.991 5.468 ( 7.871 4.544 0.000) 9.089 9.792 ( -0.502 -0.290 0.000) 0.580 11.407 ( -1.021 -0.589 0.000) 1.179 12.353 ( -2.655 -1.533 0.000) 3.066 ======================= 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.852 ( 16.097 9.294 0.000) 18.588 1.040 ( 18.925 10.926 0.000) 21.853 1.746 ( 31.260 18.048 0.000) 36.095 2.813 ( 7.551 4.360 0.000) 8.719 4.302 ( 4.803 2.773 0.000) 5.546 5.737 ( 14.812 8.551 0.000) 17.103 9.777 ( -0.720 -0.416 0.000) 0.832 11.373 ( -1.886 -1.089 0.000) 2.177 12.255 ( -5.798 -3.347 0.000) 6.695 ======================= 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.192 ( 12.832 7.409 0.000) 14.818 1.434 ( 14.769 8.527 0.000) 17.054 2.375 ( 22.346 12.901 0.000) 25.802 2.995 ( 7.571 4.371 0.000) 8.742 4.424 ( 5.226 3.017 0.000) 6.034 6.145 ( 19.560 11.293 0.000) 22.586 9.761 ( -0.584 -0.337 0.000) 0.675 11.324 ( -2.174 -1.255 0.000) 2.511 12.079 ( -9.207 -5.316 0.000) 10.631 ======================= 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.449 ( 9.017 5.206 0.000) 10.412 1.738 ( 11.266 6.504 0.000) 13.009 2.799 ( 14.093 8.137 0.000) 16.273 3.155 ( 5.881 3.396 0.000) 6.791 4.527 ( 3.218 1.858 0.000) 3.716 6.623 ( 20.422 11.790 0.000) 23.581 9.751 ( -0.269 -0.155 0.000) 0.311 11.277 ( -1.661 -0.959 0.000) 1.918 11.832 ( -11.561 -6.675 0.000) 13.349 ======================= 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.614 ( 5.131 2.963 0.000) 5.925 1.963 ( 7.817 4.513 0.000) 9.026 3.051 ( 7.699 4.445 0.000) 8.890 3.266 ( 3.586 2.070 0.000) 4.140 4.570 ( 0.561 0.324 0.000) 0.647 7.061 ( 16.100 9.295 0.000) 18.590 9.748 ( -0.026 -0.015 0.000) 0.030 11.249 ( -0.783 -0.452 0.000) 0.904 11.563 ( -10.583 -6.110 0.000) 12.221 ======================= 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.692 ( 1.623 0.937 0.000) 1.874 2.092 ( 2.963 1.711 0.000) 3.421 3.169 ( 2.471 1.426 0.000) 2.853 3.322 ( 1.196 0.690 0.000) 1.381 4.568 ( -0.286 -0.165 0.000) 0.330 7.333 ( 6.292 3.633 0.000) 7.265 9.748 ( 0.024 0.014 0.000) 0.028 11.238 ( -0.177 -0.102 0.000) 0.204 11.377 ( -4.492 -2.593 0.000) 5.186 ======================= 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.749 ( 9.741 16.872 0.000) 19.482 0.946 ( 12.593 21.811 0.000) 25.185 1.529 ( 18.814 32.588 0.000) 37.629 2.767 ( 4.124 7.143 0.000) 8.248 4.291 ( 3.491 6.046 0.000) 6.982 5.635 ( 6.802 11.782 0.000) 13.605 9.781 ( -0.397 -0.687 0.000) 0.793 11.389 ( -0.752 -1.303 0.000) 1.505 12.284 ( -3.049 -5.281 0.000) 6.098 ======================= 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| 1.082 ( 10.636 12.471 0.000) 16.391 1.361 ( 10.713 18.511 0.000) 21.387 2.152 ( 20.126 20.450 0.000) 28.693 2.933 ( 5.938 7.067 0.000) 9.231 4.444 ( 3.352 11.011 0.000) 11.510 5.942 ( 14.493 11.505 0.000) 18.504 9.767 ( -0.496 -0.503 0.000) 0.707 11.353 ( -1.931 -0.930 0.000) 2.143 12.138 ( -5.835 -8.112 0.000) 9.992 ======================= 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.362 ( 8.690 8.868 0.000) 12.416 1.692 ( 7.446 15.391 0.000) 17.098 2.628 ( 15.522 11.368 0.000) 19.240 3.103 ( 5.462 6.027 0.000) 8.133 4.607 ( -0.349 14.858 0.000) 14.862 6.361 ( 20.165 9.687 0.000) 22.371 9.756 ( -0.361 -0.095 0.000) 0.373 11.306 ( -2.367 -0.770 0.000) 2.489 11.919 ( -8.308 -10.182 0.000) 13.141 ======================= 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.562 ( 5.635 5.737 0.000) 8.042 1.939 ( 4.775 12.228 0.000) 13.127 2.935 ( 10.319 4.925 0.000) 11.434 3.238 ( 3.619 4.773 0.000) 5.990 4.716 ( -5.558 16.579 0.000) 17.486 6.813 ( 20.600 6.723 0.000) 21.669 9.752 ( -0.216 0.339 0.000) 0.402 11.261 ( -1.564 -1.015 0.000) 1.865 11.656 ( -9.200 -10.171 0.000) 13.714 ======================= 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.677 ( 2.402 3.148 0.000) 3.959 2.104 ( 1.673 8.702 0.000) 8.861 3.102 ( 5.873 0.430 0.000) 5.889 3.324 ( 1.229 3.547 0.000) 3.754 4.754 ( -9.328 17.439 0.000) 19.777 7.175 ( 14.448 1.769 0.000) 14.556 9.753 ( -0.229 0.601 0.000) 0.643 11.233 ( -0.199 -1.339 0.000) 1.353 11.422 ( -6.235 -7.347 0.000) 9.636 ======================= 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.714 ( -0.664 1.150 0.000) 1.327 2.165 ( -2.999 5.195 0.000) 5.999 3.155 ( 1.568 -2.715 0.000) 3.135 3.353 ( -1.277 2.211 0.000) 2.554 4.758 ( -10.368 17.958 0.000) 20.736 7.317 ( 3.013 -5.219 0.000) 6.027 9.755 ( -0.361 0.625 0.000) 0.722 11.225 ( 0.712 -1.234 0.000) 1.425 11.323 ( 1.413 -2.447 0.000) 2.826 ======================= 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.330 ( 6.630 11.484 0.000) 13.260 1.721 ( 9.288 16.088 0.000) 18.577 2.537 ( 9.719 16.833 0.000) 19.437 3.085 ( 4.397 7.617 0.000) 8.795 4.716 ( 8.476 14.682 0.000) 16.953 6.207 ( 8.532 14.777 0.000) 17.063 9.759 ( -0.126 -0.217 0.000) 0.251 11.327 ( -1.045 -1.810 0.000) 2.090 11.936 ( -6.656 -11.528 0.000) 13.312 ======================= 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.535 ( 5.313 7.813 0.000) 9.448 2.002 ( 5.142 14.283 0.000) 15.180 2.821 ( 8.509 7.424 0.000) 11.292 3.231 ( 3.832 6.380 0.000) 7.442 4.990 ( 0.768 21.469 0.000) 21.483 6.541 ( 15.695 8.307 0.000) 17.758 9.759 ( -0.013 0.378 0.000) 0.378 11.279 ( -2.008 -2.077 0.000) 2.889 11.686 ( -7.549 -12.134 0.000) 14.291 ======================= 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.671 ( 2.930 4.751 0.000) 5.581 2.193 ( 1.325 12.224 0.000) 12.296 2.991 ( 6.468 0.541 0.000) 6.490 3.342 ( 1.997 5.256 0.000) 5.622 5.161 ( -8.545 25.773 0.000) 27.152 6.894 ( 17.530 1.287 0.000) 17.577 9.764 ( -0.080 0.835 0.000) 0.838 11.223 ( -1.211 -2.830 0.000) 3.078 11.445 ( -6.447 -9.891 0.000) 11.807 ======================= 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.736 ( 0.354 2.498 0.000) 2.523 2.294 ( -2.655 9.682 0.000) 10.040 3.067 ( 4.243 -3.731 0.000) 5.650 3.403 ( -0.537 4.007 0.000) 4.043 5.227 ( -14.188 27.327 0.000) 30.791 7.131 ( 10.979 -6.051 0.000) 12.536 9.770 ( -0.364 0.971 0.000) 1.037 11.189 ( 0.780 -2.886 0.000) 2.990 11.286 ( -1.919 -5.614 0.000) 5.933 ======================= 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.669 ( 3.085 5.343 0.000) 6.169 2.255 ( 6.073 10.519 0.000) 12.147 2.927 ( 1.824 3.159 0.000) 3.648 3.353 ( 3.121 5.406 0.000) 6.242 5.411 ( 10.637 18.423 0.000) 21.273 6.708 ( 4.772 8.265 0.000) 9.544 9.770 ( 0.392 0.678 0.000) 0.783 11.220 ( -2.061 -3.569 0.000) 4.121 11.451 ( -6.087 -10.543 0.000) 12.173 ======================= 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.750 ( 1.274 2.932 0.000) 3.197 2.424 ( 0.868 10.342 0.000) 10.379 2.960 ( 2.487 -3.309 0.000) 4.139 3.441 ( 1.196 4.208 0.000) 4.374 5.694 ( -4.016 25.109 0.000) 25.428 6.866 ( 9.906 -3.867 0.000) 10.634 9.783 ( 0.095 0.888 0.000) 0.893 11.154 ( -0.759 -3.416 0.000) 3.500 11.273 ( -3.720 -6.967 0.000) 7.898 ======================= 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.775 ( -0.773 1.339 0.000) 1.546 2.486 ( -5.077 8.794 0.000) 10.155 2.963 ( 3.482 -6.032 0.000) 6.965 3.474 ( -1.549 2.683 0.000) 3.098 5.784 ( -14.685 25.435 0.000) 29.370 6.941 ( 6.913 -11.975 0.000) 13.827 9.789 ( -0.425 0.735 0.000) 0.849 11.133 ( 1.202 -2.082 0.000) 2.404 11.197 ( 1.853 -3.211 0.000) 3.707 ======================= 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.789 ( 0.556 0.963 0.000) 1.112 2.599 ( 3.515 6.088 0.000) 7.030 2.876 ( -2.251 -3.899 0.000) 4.503 3.502 ( 0.935 1.619 0.000) 1.869 6.107 ( 7.603 13.168 0.000) 15.205 6.759 ( -3.002 -5.200 0.000) 6.005 9.797 ( 0.219 0.380 0.000) 0.439 11.113 ( -0.110 -0.191 0.000) 0.220 11.162 ( -2.220 -3.846 0.000) 4.441 ======================= 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.438 ( 0.000 0.000 21.886) 21.886 0.438 ( 0.000 0.000 21.886) 21.886 0.920 ( 0.000 0.000 44.777) 44.777 4.145 ( 0.000 0.000 -3.680) 3.680 4.145 ( 0.000 0.000 -3.680) 3.680 5.538 ( -0.000 -0.000 10.315) 10.315 9.843 ( 0.000 0.000 0.881) 0.881 11.420 ( 0.000 -0.000 0.072) 0.072 11.420 ( -0.000 0.000 0.072) 0.072 ======================= 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.704 ( 17.679 10.207 13.588) 24.522 0.782 ( 15.858 9.155 18.343) 25.918 1.198 ( 23.894 13.795 28.427) 39.615 3.994 ( -11.298 -6.523 26.284) 29.343 4.179 ( 2.847 1.644 -3.365) 4.704 5.469 ( -0.737 -0.426 5.944) 6.005 9.824 ( -1.147 -0.662 1.782) 2.220 11.407 ( -1.068 -0.616 0.015) 1.233 11.647 ( 13.636 7.873 -17.028) 23.192 ======================= 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| 1.093 ( 11.927 6.886 18.449) 23.022 1.128 ( 17.392 10.041 8.367) 21.756 1.836 ( 26.875 15.516 10.866) 32.880 3.737 ( -8.745 -5.049 49.140) 50.167 4.275 ( 5.188 2.995 -2.586) 6.525 5.616 ( 12.731 7.350 -7.492) 16.500 9.801 ( -0.816 -0.471 2.035) 2.242 11.371 ( -1.944 -1.122 -0.109) 2.247 11.874 ( 4.721 2.726 -20.828) 21.529 ======================= 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.350 ( 9.880 5.704 13.231) 17.470 1.497 ( 13.963 8.062 5.987) 17.199 2.389 ( 20.005 11.550 3.486) 23.361 3.603 ( -3.070 -1.772 43.025) 43.170 4.407 ( 5.650 3.262 -1.626) 6.724 6.012 ( 19.881 11.479 -10.904) 25.415 9.786 ( -0.434 -0.251 2.325) 2.379 11.321 ( -2.224 -1.284 -0.230) 2.578 11.878 ( -3.799 -2.193 -14.390) 15.043 ======================= 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.550 ( 7.210 4.163 8.909) 12.193 1.783 ( 10.579 6.108 4.348) 12.966 2.772 ( 12.850 7.419 -0.995) 14.871 3.570 ( -0.188 -0.109 32.538) 32.539 4.520 ( 3.694 2.133 -0.574) 4.304 6.499 ( 20.736 11.972 -10.597) 26.184 9.780 ( -0.080 -0.046 2.755) 2.757 11.274 ( -1.709 -0.987 -0.338) 2.002 11.725 ( -8.580 -4.954 -8.509) 13.060 ======================= 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.686 ( 4.316 2.492 6.416) 8.124 1.992 ( 7.175 4.143 2.837) 8.757 3.002 ( 6.966 4.022 -4.043) 9.002 3.579 ( 0.691 0.399 25.850) 25.862 4.574 ( 1.039 0.600 0.547) 1.318 6.939 ( 15.962 9.215 -10.614) 21.268 9.781 ( 0.080 0.046 3.099) 3.100 11.244 ( -0.830 -0.479 -0.449) 1.058 11.511 ( -8.820 -5.092 -4.352) 11.076 ======================= 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.753 ( 1.426 0.823 5.476) 5.718 2.110 ( 2.669 1.541 1.741) 3.540 3.107 ( 2.178 1.257 -5.756) 6.281 3.594 ( 0.389 0.225 23.040) 23.045 4.582 ( -0.061 -0.035 1.395) 1.396 7.207 ( 6.130 3.539 -11.094) 13.160 9.782 ( 0.039 0.023 3.227) 3.227 11.233 ( -0.198 -0.114 -0.531) 0.578 11.353 ( -3.831 -2.212 -1.956) 4.836 ======================= 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| 1.013 ( 7.267 12.587 18.950) 23.882 1.040 ( 11.240 19.468 9.241) 24.305 1.654 ( 15.717 27.223 14.277) 34.525 3.787 ( -6.309 -10.928 46.033) 47.731 4.266 ( 3.878 6.716 -2.156) 8.049 5.535 ( 4.850 8.401 -4.572) 10.724 9.806 ( -0.530 -0.918 1.984) 2.250 11.387 ( -0.795 -1.376 -0.092) 1.592 11.831 ( 4.122 7.139 -21.678) 23.192 ======================= 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.263 ( 8.190 9.381 14.793) 19.336 1.424 ( 10.250 17.176 6.108) 20.913 2.193 ( 17.722 18.210 5.935) 26.094 3.620 ( -3.043 -5.068 45.051) 45.438 4.431 ( 3.418 11.887 -0.877) 12.399 5.811 ( 14.234 11.567 -10.252) 21.012 9.790 ( -0.423 -0.482 2.148) 2.242 11.351 ( -1.959 -0.999 -0.225) 2.211 11.887 ( -0.667 -2.104 -16.782) 16.926 ======================= 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.479 ( 6.930 6.785 10.111) 14.010 1.737 ( 7.201 14.303 4.369) 16.599 2.618 ( 13.980 10.335 0.731) 17.401 3.566 ( -0.117 -1.179 35.104) 35.124 4.602 ( -0.381 15.523 -0.224) 15.529 6.236 ( 20.433 10.284 -10.563) 25.197 9.782 ( -0.168 -0.066 2.473) 2.480 11.302 ( -2.382 -0.798 -0.306) 2.531 11.782 ( -5.245 -7.164 -10.567) 13.802 ======================= 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.639 ( 4.765 4.421 6.851) 9.444 1.970 ( 4.518 11.292 2.953) 12.516 2.896 ( 9.382 4.435 -2.782) 10.744 3.574 ( 0.754 0.861 27.180) 27.205 4.718 ( -5.353 16.882 0.310) 17.713 6.695 ( 20.473 7.105 -10.191) 23.948 9.782 ( -0.009 0.278 2.815) 2.829 11.257 ( -1.602 -0.963 -0.337) 1.899 11.584 ( -7.311 -8.526 -5.894) 12.684 ======================= 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.736 ( 2.272 2.351 5.277) 6.208 2.122 ( 1.481 8.032 1.756) 8.355 3.046 ( 5.298 0.232 -5.040) 7.316 3.598 ( 0.298 1.628 23.119) 23.178 4.762 ( -8.960 17.368 0.905) 19.564 7.055 ( 14.021 2.028 -10.508) 17.639 9.786 ( -0.083 0.440 3.008) 3.041 11.230 ( -0.270 -1.155 -0.326) 1.230 11.388 ( -5.174 -6.600 -2.829) 8.851 ======================= 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.767 ( -0.360 0.623 4.875) 4.928 2.177 ( -2.810 4.867 1.241) 5.755 3.093 ( 1.490 -2.581 -5.846) 6.562 3.610 ( -0.788 1.365 21.950) 22.007 4.769 ( -10.159 17.596 1.194) 20.353 7.194 ( 2.779 -4.815 -10.819) 12.163 9.787 ( -0.248 0.430 3.045) 3.085 11.223 ( 0.541 -0.936 -0.254) 1.111 11.303 ( 1.471 -2.549 -1.733) 3.416 ======================= 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.452 ( 5.129 8.883 10.466) 14.655 1.761 ( 8.749 15.154 3.877) 17.923 2.537 ( 8.731 15.123 1.617) 17.537 3.559 ( -0.552 -0.956 35.414) 35.431 4.716 ( 8.663 15.005 0.354) 17.330 6.084 ( 8.860 15.346 -10.454) 20.574 9.784 ( -0.082 -0.142 2.308) 2.314 11.324 ( -1.051 -1.820 -0.303) 2.123 11.786 ( -4.312 -7.468 -11.349) 14.254 ======================= 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.612 ( 4.310 6.084 6.883) 10.147 2.026 ( 4.913 13.349 2.344) 14.416 2.792 ( 7.671 6.627 -1.804) 10.296 3.570 ( 1.008 1.521 27.295) 27.356 4.992 ( 0.699 21.476 0.444) 21.491 6.428 ( 15.862 8.808 -9.775) 20.609 9.785 ( 0.097 0.332 2.472) 2.496 11.277 ( -1.986 -1.982 -0.234) 2.815 11.603 ( -5.902 -9.906 -6.813) 13.393 ======================= 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.723 ( 2.605 3.709 4.769) 6.580 2.206 ( 1.236 11.403 1.207) 11.533 2.944 ( 5.844 0.233 -4.138) 7.165 3.611 ( 0.860 2.634 22.526) 22.696 5.162 ( -8.394 25.431 0.251) 26.781 6.784 ( 17.246 1.692 -9.560) 19.791 9.792 ( 0.059 0.642 2.572) 2.652 11.222 ( -1.196 -2.617 -0.049) 2.878 11.402 ( -5.396 -8.727 -3.636) 10.886 ======================= 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.779 ( 0.504 1.851 3.974) 4.413 2.299 ( -2.536 9.083 0.505) 9.444 3.010 ( 3.905 -3.657 -5.316) 7.542 3.643 ( -0.440 2.547 20.485) 20.648 5.227 ( -13.836 26.735 0.204) 30.104 7.019 ( 10.542 -5.528 -9.821) 15.432 9.798 ( -0.215 0.699 2.586) 2.687 11.191 ( 0.843 -2.607 0.241) 2.750 11.263 ( -1.666 -5.255 -2.017) 5.870 ======================= 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.719 ( 2.508 4.344 4.544) 6.768 2.263 ( 5.708 9.887 0.823) 11.446 2.884 ( 1.505 2.606 -3.691) 4.762 3.617 ( 1.602 2.774 22.197) 22.427 5.409 ( 10.455 18.109 0.050) 20.910 6.602 ( 4.911 8.507 -9.224) 13.475 9.795 ( 0.341 0.590 2.348) 2.444 11.220 ( -1.971 -3.414 0.050) 3.942 11.405 ( -5.241 -9.077 -3.929) 11.194 ======================= 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.786 ( 1.161 2.420 3.428) 4.354 2.422 ( 0.788 9.792 -0.168) 9.825 2.909 ( 2.224 -3.399 -4.697) 6.209 3.669 ( 0.737 2.732 19.444) 19.649 5.686 ( -4.083 24.709 -0.571) 25.050 6.764 ( 9.845 -3.632 -9.013) 13.833 9.806 ( 0.124 0.686 2.188) 2.296 11.156 ( -0.726 -3.425 0.178) 3.506 11.250 ( -3.247 -6.081 -1.910) 7.153 ======================= 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.809 ( -0.603 1.045 3.154) 3.377 2.480 ( -4.861 8.420 -0.535) 9.737 2.909 ( 3.389 -5.870 -4.987) 8.415 3.690 ( -1.069 1.851 18.574) 18.696 5.773 ( -14.473 25.066 -0.845) 28.957 6.838 ( 6.730 -11.659 -8.980) 16.182 9.811 ( -0.304 0.527 2.111) 2.197 11.135 ( 1.376 -2.384 -0.020) 2.753 11.184 ( 1.510 -2.616 -0.931) 3.160 ======================= 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.819 ( 0.483 0.836 2.905) 3.061 2.592 ( 3.568 6.180 -0.797) 7.181 2.823 ( -2.420 -4.191 -4.918) 6.899 3.710 ( 0.651 1.127 17.855) 17.903 6.098 ( 7.825 13.554 -0.590) 15.661 6.654 ( -3.270 -5.664 -9.305) 11.374 9.817 ( 0.169 0.293 1.903) 1.932 11.110 ( -0.447 -0.774 -0.625) 1.091 11.156 ( -1.712 -2.965 -0.066) 3.424 ======================= 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.897 ( 0.000 0.000 23.146) 23.146 0.897 ( 0.000 0.000 23.146) 23.146 1.818 ( 0.000 0.000 43.030) 43.030 4.037 ( 0.000 0.000 -6.848) 6.848 4.037 ( 0.000 -0.000 -6.848) 6.848 5.796 ( 0.000 0.000 13.467) 13.467 9.872 ( 0.000 0.000 1.938) 1.938 11.423 ( 0.000 0.000 0.180) 0.180 11.423 ( 0.000 0.000 0.180) 0.180 ======================= 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| 1.050 ( 11.595 6.694 19.438) 23.603 1.128 ( 15.172 8.759 17.637) 24.859 1.917 ( 8.560 4.942 37.932) 39.199 4.080 ( 3.655 2.110 -6.183) 7.486 4.157 ( 9.371 5.410 -1.534) 10.928 5.686 ( -7.897 -4.559 12.808) 15.722 9.864 ( -0.566 -0.327 2.276) 2.368 11.408 ( -1.190 -0.687 0.090) 1.377 11.489 ( 4.914 2.837 -2.987) 6.412 ======================= 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.363 ( 13.858 8.001 14.101) 21.329 1.441 ( 10.498 6.061 14.979) 19.269 2.202 ( 14.448 8.342 23.516) 28.833 4.201 ( 6.442 3.719 -4.487) 8.687 4.378 ( 6.292 3.632 16.366) 17.906 5.556 ( 0.490 0.283 2.313) 2.381 9.852 ( -0.398 -0.230 2.872) 2.908 11.369 ( -2.099 -1.212 -0.111) 2.427 11.603 ( 3.792 2.189 -7.381) 8.582 ======================= 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.633 ( 6.354 3.669 13.019) 14.944 1.666 ( 11.679 6.743 10.187) 16.901 2.527 ( 12.332 7.120 9.458) 17.094 4.365 ( 7.043 4.066 -2.337) 8.462 4.373 ( -5.281 -3.049 30.489) 31.093 5.792 ( 17.906 10.338 -9.296) 22.670 9.848 ( 0.072 0.042 3.523) 3.524 11.315 ( -2.361 -1.363 -0.320) 2.745 11.638 ( -0.965 -0.557 -8.312) 8.386 ======================= 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.755 ( 4.335 2.503 10.033) 11.212 1.904 ( 8.657 4.998 7.160) 12.296 2.761 ( 7.651 4.417 -0.401) 8.843 4.230 ( -5.605 -3.236 29.142) 29.853 4.511 ( 5.153 2.975 -0.096) 5.951 6.270 ( 20.828 12.025 -10.341) 26.179 9.854 ( 0.407 0.235 4.201) 4.227 11.265 ( -1.844 -1.065 -0.513) 2.190 11.562 ( -5.154 -2.976 -6.423) 8.756 ======================= 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.840 ( 2.852 1.646 7.794) 8.461 2.071 ( 5.536 3.196 4.637) 7.897 2.892 ( 3.791 2.189 -6.562) 7.888 4.132 ( -2.766 -1.597 25.281) 25.482 4.599 ( 2.368 1.367 2.014) 3.396 6.705 ( 15.404 8.893 -10.723) 20.769 9.864 ( 0.365 0.211 4.731) 4.750 11.232 ( -0.959 -0.553 -0.705) 1.312 11.424 ( -5.970 -3.447 -3.518) 7.739 ======================= 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.886 ( 1.050 0.606 6.829) 6.935 2.160 ( 1.961 1.132 3.005) 3.762 2.948 ( 1.113 0.643 -9.436) 9.523 4.093 ( -0.736 -0.425 22.874) 22.890 4.630 ( 0.525 0.303 3.435) 3.488 6.958 ( 5.698 3.290 -11.522) 13.268 9.869 ( 0.107 0.062 4.985) 4.987 11.218 ( -0.255 -0.147 -0.844) 0.894 11.316 ( -2.637 -1.522 -1.316) 3.317 ======================= 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.290 ( 8.660 14.999 14.998) 22.911 1.364 ( 6.887 11.929 15.258) 20.556 2.110 ( 7.702 13.341 28.003) 31.960 4.161 ( 3.126 5.414 -3.330) 7.084 4.337 ( 6.031 10.446 8.551) 14.786 5.563 ( -2.697 -4.670 6.783) 8.665 9.854 ( -0.301 -0.522 2.700) 2.767 11.386 ( -0.903 -1.563 -0.076) 1.807 11.572 ( 2.482 4.300 -6.457) 8.145 ======================= 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.565 ( 6.019 6.898 13.338) 16.178 1.602 ( 8.714 13.056 11.006) 19.171 2.407 ( 10.527 10.840 14.355) 20.842 4.256 ( 2.680 -0.236 9.621) 9.991 4.504 ( 1.008 10.098 14.903) 18.030 5.625 ( 10.022 7.362 -6.619) 14.087 9.847 ( -0.069 -0.193 3.232) 3.238 11.345 ( -2.038 -1.161 -0.281) 2.362 11.628 ( 0.859 -0.260 -8.357) 8.405 ======================= 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.706 ( 4.470 3.604 10.823) 12.252 1.862 ( 6.284 11.439 7.525) 15.066 2.663 ( 8.612 5.846 3.327) 10.927 4.221 ( -1.346 -6.255 24.645) 25.462 4.629 ( -2.287 16.810 4.675) 17.597 6.010 ( 20.370 10.612 -10.200) 25.132 9.849 ( 0.329 0.072 3.803) 3.818 11.294 ( -2.446 -0.823 -0.406) 2.612 11.589 ( -2.485 -4.439 -7.378) 8.962 ======================= 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.802 ( 3.326 2.117 8.190) 9.089 2.053 ( 3.706 8.909 4.943) 10.842 2.826 ( 5.665 1.802 -4.151) 7.250 4.139 ( -1.493 -3.754 24.993) 25.318 4.741 ( -5.324 17.730 2.356) 18.662 6.471 ( 19.990 7.603 -10.110) 23.656 9.858 ( 0.496 0.166 4.319) 4.351 11.248 ( -1.752 -0.770 -0.474) 1.971 11.466 ( -4.541 -6.256 -4.826) 9.113 ======================= 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.865 ( 2.017 0.927 6.701) 7.059 2.172 ( 0.996 6.313 2.978) 7.051 2.908 ( 3.255 -0.708 -8.254) 8.901 4.097 ( -0.622 -1.163 22.766) 22.804 4.798 ( -8.139 17.161 2.733) 19.189 6.821 ( 12.990 2.445 -10.723) 17.020 9.866 ( 0.295 0.087 4.639) 4.649 11.220 ( -0.530 -0.714 -0.566) 1.054 11.332 ( -3.221 -5.379 -2.185) 6.640 ======================= 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.888 ( 0.166 -0.287 6.295) 6.303 2.213 ( -2.287 3.961 2.188) 5.070 2.932 ( 1.257 -2.177 -9.519) 9.846 4.087 ( 0.040 -0.069 21.793) 21.793 4.812 ( -9.609 16.643 3.079) 19.463 6.951 ( 2.294 -3.975 -11.196) 12.100 9.869 ( 0.021 -0.036 4.734) 4.734 11.214 ( 0.178 -0.309 -0.707) 0.792 11.271 ( 1.598 -2.769 -0.934) 3.331 ======================= 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.683 ( 2.964 5.133 10.956) 12.457 1.874 ( 7.238 12.536 6.959) 16.061 2.610 ( 5.055 8.755 5.139) 11.341 4.199 ( -2.350 -4.071 23.347) 23.816 4.766 ( 8.136 14.092 6.036) 17.356 5.859 ( 8.800 15.242 -10.214) 20.349 9.846 ( 0.070 0.121 3.582) 3.585 11.317 ( -1.033 -1.789 -0.336) 2.093 11.582 ( -2.429 -4.207 -7.676) 9.084 ======================= 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.775 ( 2.686 3.144 8.111) 9.104 2.094 ( 4.198 10.832 4.174) 12.344 2.753 ( 4.460 2.951 -2.184) 5.777 4.127 ( -1.314 -2.980 24.219) 24.437 5.023 ( 0.014 20.984 2.888) 21.182 6.213 ( 16.168 9.467 -9.730) 21.112 9.852 ( 0.364 0.255 3.852) 3.878 11.272 ( -1.952 -1.606 -0.127) 2.531 11.467 ( -3.726 -7.152 -5.684) 9.866 ======================= 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.842 ( 2.070 1.826 6.217) 6.803 2.240 ( 0.927 9.155 2.100) 9.439 2.832 ( 3.673 -1.146 -6.629) 7.665 4.095 ( -0.364 -0.828 22.047) 22.066 5.180 ( -8.205 24.401 1.676) 25.798 6.573 ( 16.590 2.380 -9.469) 19.250 9.861 ( 0.399 0.207 4.008) 4.033 11.223 ( -1.308 -2.008 0.087) 2.397 11.326 ( -3.496 -6.947 -3.258) 8.432 ======================= 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.881 ( 0.789 0.676 5.489) 5.586 2.313 ( -2.177 7.331 0.803) 7.689 2.865 ( 2.838 -3.377 -8.495) 9.572 4.089 ( -0.147 0.187 20.433) 20.434 5.242 ( -12.990 25.173 1.377) 28.360 6.801 ( 9.575 -4.460 -9.862) 14.451 9.868 ( 0.156 0.059 4.054) 4.058 11.195 ( 1.003 -2.395 -0.379) 2.624 11.224 ( -1.281 -4.213 -0.958) 4.506 ======================= 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.833 ( 1.486 2.574 5.968) 6.668 2.287 ( 4.644 8.044 1.468) 9.403 2.786 ( 0.296 0.513 -5.771) 5.802 4.094 ( -0.332 -0.576 21.717) 21.727 5.423 ( 9.840 17.044 1.303) 19.723 6.399 ( 5.152 8.923 -9.039) 13.706 9.859 ( 0.227 0.394 3.692) 3.720 11.224 ( -1.783 -3.088 0.493) 3.599 11.320 ( -3.841 -6.652 -3.804) 8.572 ======================= 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.876 ( 0.985 1.494 4.951) 5.264 2.415 ( 0.505 7.944 -0.602) 7.983 2.784 ( 1.439 -3.430 -7.171) 8.078 4.093 ( 0.074 0.361 19.437) 19.441 5.679 ( -4.349 23.544 -0.139) 23.942 6.566 ( 9.741 -3.159 -8.563) 13.348 9.866 ( 0.199 0.209 3.469) 3.481 11.155 ( -1.138 -4.408 -0.484) 4.578 11.217 ( -1.898 -3.407 -0.732) 3.968 ======================= 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.893 ( -0.287 0.497 4.726) 4.761 2.461 ( -4.022 6.967 -1.525) 8.188 2.777 ( 2.876 -4.981 -7.350) 9.333 4.097 ( -0.223 0.386 18.568) 18.573 5.758 ( -13.867 24.018 -0.710) 27.743 6.642 ( 6.311 -10.933 -8.435) 15.183 9.869 ( -0.015 0.026 3.363) 3.363 11.120 ( 1.904 -3.297 -1.366) 4.045 11.183 ( 0.622 -1.078 0.978) 1.583 ======================= 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.899 ( 0.353 0.612 4.492) 4.547 2.558 ( 3.319 5.749 -3.049) 7.305 2.701 ( -2.497 -4.324 -6.157) 7.927 4.101 ( 0.155 0.268 17.827) 17.830 6.088 ( 8.708 15.082 -1.218) 17.458 6.449 ( -4.323 -7.487 -8.166) 11.893 9.869 ( 0.049 0.085 3.055) 3.056 11.085 ( -1.126 -1.950 -1.500) 2.705 11.172 ( -0.684 -1.185 1.435) 1.983 ======================= 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.374 ( 0.000 0.000 23.174) 23.174 1.374 ( -0.000 0.000 23.174) 23.174 2.658 ( 0.000 -0.000 38.705) 38.705 3.871 ( 0.000 0.000 -9.223) 9.223 3.871 ( 0.000 0.000 -9.223) 9.223 6.034 ( 0.000 -0.000 8.935) 8.935 9.921 ( 0.000 0.000 2.815) 2.815 11.427 ( 0.000 0.000 0.276) 0.276 11.427 ( 0.000 0.000 0.276) 0.276 ======================= 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.468 ( 7.378 4.260 20.809) 22.485 1.512 ( 9.782 5.648 19.510) 22.543 2.669 ( 0.997 0.576 34.792) 34.811 3.932 ( 5.096 2.942 -8.157) 10.058 4.064 ( 15.803 9.124 -6.447) 19.354 5.920 ( -9.097 -5.252 9.170) 13.943 9.918 ( -0.216 -0.125 2.953) 2.964 11.411 ( -1.339 -0.773 0.187) 1.557 11.457 ( 2.247 1.297 -0.717) 2.692 ======================= 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.677 ( 9.619 5.553 15.945) 19.432 1.734 ( 7.940 4.584 14.207) 16.909 2.702 ( 1.591 0.919 23.654) 23.726 4.098 ( 8.729 5.040 -5.425) 11.446 4.538 ( 21.868 12.625 2.127) 25.340 5.669 ( -10.108 -5.836 6.933) 13.576 9.915 ( 0.064 0.037 3.281) 3.282 11.368 ( -2.301 -1.329 -0.022) 2.657 11.511 ( 1.838 1.061 -2.536) 3.307 ======================= 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.872 ( 4.061 2.345 10.910) 11.875 1.890 ( 8.293 4.788 11.220) 14.750 2.732 ( 0.842 0.486 9.412) 9.462 4.318 ( 9.508 5.490 -2.066) 11.172 4.860 ( 0.892 0.515 18.172) 18.202 5.651 ( 12.871 7.431 -4.608) 15.560 9.923 ( 0.611 0.353 3.725) 3.791 11.309 ( -2.552 -1.473 -0.256) 2.958 11.525 ( -0.841 -0.485 -3.392) 3.529 ======================= 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.940 ( 2.072 1.196 8.161) 8.505 2.058 ( 5.977 3.451 7.469) 10.170 2.742 ( 0.066 0.038 -1.908) 1.909 4.521 ( 7.444 4.298 1.086) 8.664 4.723 ( -8.187 -4.727 19.344) 21.530 6.107 ( 21.128 12.198 -5.576) 25.026 9.942 ( 0.934 0.539 4.245) 4.380 11.254 ( -2.032 -1.173 -0.484) 2.396 11.472 ( -3.504 -2.023 -2.680) 4.853 ======================= 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.980 ( 1.455 0.840 5.783) 6.022 2.170 ( 3.598 2.077 4.796) 6.345 2.741 ( -0.058 -0.033 -8.109) 8.110 4.557 ( -5.533 -3.194 16.232) 17.444 4.658 ( 4.159 2.401 3.614) 6.011 6.539 ( 14.830 8.562 -5.418) 17.961 9.963 ( 0.739 0.426 4.742) 4.818 11.217 ( -1.130 -0.652 -0.703) 1.482 11.378 ( -3.974 -2.294 -0.980) 4.692 ======================= 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| 2.006 ( 0.648 0.374 4.552) 4.613 2.227 ( 1.209 0.698 3.337) 3.617 2.740 ( -0.010 -0.006 -10.493) 10.493 4.470 ( -1.856 -1.071 13.929) 14.092 4.719 ( 1.240 0.716 5.096) 5.293 6.777 ( 5.246 3.029 -5.991) 8.520 9.974 ( 0.242 0.140 5.069) 5.077 11.200 ( -0.330 -0.190 -0.855) 0.936 11.307 ( -1.720 -0.993 0.443) 2.035 ======================= 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.624 ( 5.804 10.052 17.019) 20.601 1.676 ( 5.062 8.768 15.564) 18.568 2.691 ( 0.836 1.448 27.247) 27.299 4.056 ( 5.133 8.890 -6.043) 11.912 4.393 ( 12.504 21.658 -1.178) 25.036 5.735 ( -6.694 -11.594 8.327) 15.766 9.915 ( -0.056 -0.097 3.178) 3.180 11.385 ( -1.030 -1.785 0.025) 2.061 11.494 ( 1.109 1.921 -2.072) 3.036 ======================= 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.816 ( 4.589 5.056 11.888) 13.709 1.846 ( 6.041 8.372 12.289) 16.050 2.718 ( 0.899 0.428 14.572) 14.606 4.268 ( 6.678 10.873 -2.219) 12.951 4.822 ( 10.010 12.951 11.194) 19.830 5.562 ( 0.820 -2.075 0.260) 2.246 9.917 ( 0.312 0.156 3.530) 3.547 11.341 ( -2.156 -1.316 -0.182) 2.533 11.517 ( 0.376 -0.793 -3.263) 3.379 ======================= 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.905 ( 2.723 1.235 9.005) 9.489 2.028 ( 4.545 7.871 8.279) 12.294 2.728 ( 0.610 -1.019 2.215) 2.513 4.473 ( 4.811 7.308 3.831) 9.551 4.874 ( -7.752 3.747 16.045) 18.209 5.846 ( 20.593 10.934 -5.800) 24.026 9.930 ( 0.832 0.251 3.948) 4.042 11.287 ( -2.596 -0.776 -0.334) 2.730 11.484 ( -1.521 -3.636 -3.182) 5.065 ======================= 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.952 ( 2.050 0.221 6.442) 6.764 2.160 ( 2.436 6.039 5.214) 8.342 2.725 ( 0.749 -1.599 -5.788) 6.051 4.533 ( -0.161 -3.436 13.457) 13.889 4.826 ( -7.197 15.037 5.960) 17.704 6.317 ( 19.540 8.317 -4.820) 21.777 9.949 ( 1.042 0.078 4.366) 4.490 11.238 ( -2.048 -0.542 -0.498) 2.176 11.399 ( -2.678 -4.968 -1.813) 5.928 ======================= 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.985 ( 1.612 -0.149 4.765) 5.032 2.237 ( 0.470 4.145 3.145) 5.224 2.724 ( 0.913 -1.555 -9.413) 9.584 4.472 ( -1.281 -3.290 13.962) 14.401 4.871 ( -7.127 16.312 4.105) 18.268 6.657 ( 11.822 3.058 -5.126) 13.243 9.964 ( 0.775 -0.274 4.700) 4.771 11.206 ( -0.914 -0.652 -0.797) 1.377 11.310 ( -1.623 -4.172 0.037) 4.476 ======================= 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| 2.000 ( 0.467 -0.809 4.257) 4.358 2.262 ( -1.515 2.623 2.378) 3.851 2.724 ( 0.870 -1.508 -10.357) 10.502 4.443 ( 0.713 -1.234 12.994) 13.072 4.891 ( -8.741 15.140 4.453) 18.041 6.779 ( 1.682 -2.915 -5.486) 6.436 9.970 ( 0.336 -0.582 4.837) 4.883 11.194 ( 0.310 -0.537 -1.027) 1.200 11.273 ( 1.248 -2.162 0.974) 2.679 ======================= 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.886 ( 1.405 2.433 9.195) 9.615 2.032 ( 5.095 8.825 8.148) 13.047 2.714 ( -0.439 -0.761 4.120) 4.212 4.473 ( 4.122 7.139 6.067) 10.236 4.988 ( 2.670 4.624 13.618) 14.628 5.692 ( 8.627 14.943 -6.026) 18.276 9.923 ( 0.240 0.416 3.799) 3.829 11.311 ( -0.952 -1.650 -0.192) 1.914 11.473 ( -2.023 -3.504 -3.399) 5.285 ======================= 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.924 ( 1.385 0.800 6.616) 6.806 2.187 ( 2.975 7.462 4.714) 9.314 2.695 ( 0.088 -1.993 -3.764) 4.260 4.516 ( -0.852 -2.101 14.355) 14.533 5.102 ( -2.086 16.871 3.994) 17.463 6.070 ( 16.814 10.824 -3.990) 20.391 9.934 ( 0.630 0.181 3.995) 4.049 11.272 ( -2.084 -0.869 0.009) 2.258 11.382 ( -2.361 -6.110 -2.610) 7.051 ======================= 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.956 ( 1.561 0.349 4.745) 5.007 2.285 ( 0.396 5.974 2.072) 6.335 2.683 ( 1.100 -2.378 -7.726) 8.158 4.461 ( -1.316 -3.210 13.875) 14.302 5.220 ( -7.613 22.110 1.733) 23.448 6.440 ( 15.651 3.745 -3.268) 16.422 9.946 ( 0.776 -0.266 4.103) 4.184 11.221 ( -2.172 -1.737 -0.473) 2.821 11.285 ( -1.347 -5.426 -0.526) 5.615 ======================= 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.981 ( 0.930 -0.106 3.981) 4.090 2.327 ( -1.532 4.516 0.383) 4.784 2.682 ( 1.464 -2.406 -8.937) 9.370 4.421 ( -0.064 -1.916 11.980) 12.132 5.278 ( -11.503 22.607 1.868) 25.434 6.661 ( 8.204 -2.692 -3.538) 9.331 9.954 ( 0.605 -0.674 4.154) 4.252 11.171 ( 0.103 -3.369 -1.233) 3.589 11.236 ( -0.001 -2.165 1.446) 2.603 ======================= 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.943 ( 0.664 1.150 4.682) 4.867 2.318 ( 2.976 5.155 1.320) 6.096 2.656 ( -0.923 -1.600 -6.771) 7.018 4.456 ( -1.737 -3.009 13.840) 14.269 5.443 ( 8.556 14.819 -0.050) 17.112 6.283 ( 5.901 10.221 -1.748) 11.931 9.937 ( 0.091 0.158 3.822) 3.827 11.240 ( -2.332 -4.039 0.785) 4.729 11.260 ( -2.121 -3.674 -1.867) 4.635 ======================= 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.969 ( 0.861 0.847 3.768) 3.956 2.392 ( -0.028 4.323 -1.903) 4.724 2.635 ( 0.691 -2.169 -6.938) 7.302 4.408 ( -0.718 -1.981 11.263) 11.458 5.664 ( -3.848 20.353 -1.703) 20.783 6.470 ( 8.982 -0.747 -0.380) 9.021 9.940 ( 0.284 -0.343 3.580) 3.608 11.138 ( -1.729 -5.102 -0.893) 5.461 11.226 ( -0.560 -1.294 1.233) 1.873 ======================= 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.982 ( -0.086 0.150 3.586) 3.590 2.411 ( -2.062 3.572 -3.463) 5.385 2.632 ( 1.381 -2.392 -6.442) 7.010 4.391 ( 0.589 -1.020 10.119) 10.188 5.731 ( -11.801 20.439 -2.094) 23.693 6.551 ( 4.518 -7.826 -0.086) 9.037 9.940 ( 0.325 -0.562 3.470) 3.530 11.095 ( 1.855 -3.213 -0.864) 3.809 11.214 ( 0.334 -0.578 1.784) 1.905 ======================= 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.984 ( 0.285 0.494 3.476) 3.522 2.454 ( 0.960 1.663 -6.590) 6.864 2.602 ( -0.564 -0.977 -3.756) 3.922 4.381 ( -0.384 -0.665 9.371) 9.402 5.992 ( 5.938 10.286 -7.375) 13.981 6.422 ( -1.730 -2.997 5.022) 6.099 9.934 ( -0.092 -0.160 3.164) 3.169 11.060 ( -1.216 -2.106 -0.774) 2.553 11.208 ( -0.285 -0.494 1.945) 2.027 ======================= 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.818 ( 0.000 0.000 19.605) 19.605 1.818 ( -0.000 0.000 19.605) 19.605 3.363 ( 0.000 0.000 29.266) 29.266 3.671 ( 0.000 0.000 -9.864) 9.864 3.671 ( 0.000 0.000 -9.864) 9.864 6.143 ( 0.000 0.000 1.961) 1.961 9.980 ( 0.000 0.000 2.656) 2.656 11.433 ( 0.000 0.000 0.263) 0.263 11.433 ( 0.000 0.000 0.263) 0.263 ======================= 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.869 ( 3.957 2.285 17.606) 18.190 1.893 ( 5.226 3.017 16.777) 17.830 3.297 ( -5.526 -3.190 25.656) 26.437 3.758 ( 7.203 4.159 -8.418) 11.834 3.921 ( 20.360 11.755 -6.974) 24.522 6.043 ( -8.228 -4.750 2.972) 9.955 9.978 ( -0.092 -0.053 2.688) 2.690 11.415 ( -1.465 -0.846 0.201) 1.703 11.449 ( 1.227 0.709 -0.166) 1.427 ======================= 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.982 ( 5.231 3.020 13.114) 14.438 2.013 ( 4.440 2.563 12.428) 13.444 3.115 ( -9.656 -5.575 15.970) 19.477 3.989 ( 11.891 6.866 -4.987) 14.609 4.540 ( 29.996 17.318 -0.917) 34.649 5.789 ( -12.380 -7.147 4.392) 14.954 9.979 ( 0.260 0.150 2.756) 2.773 11.368 ( -2.484 -1.434 0.049) 2.869 11.478 ( 0.924 0.534 -0.947) 1.426 ======================= 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| 2.089 ( 1.974 1.140 10.177) 10.429 2.099 ( 4.598 2.655 8.650) 10.149 2.872 ( -10.285 -5.938 3.944) 12.514 4.283 ( 12.558 7.251 -1.311) 14.560 5.147 ( 12.383 7.149 10.550) 17.769 5.588 ( 3.806 2.198 -2.126) 4.882 9.992 ( 0.889 0.513 2.870) 3.048 11.305 ( -2.739 -1.581 -0.134) 3.165 11.481 ( -0.799 -0.461 -1.223) 1.532 ======================= 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| 2.105 ( -0.429 -0.248 8.434) 8.449 2.193 ( 3.383 1.953 5.417) 6.678 2.668 ( -6.576 -3.797 -5.400) 9.318 4.551 ( 9.875 5.701 1.624) 11.518 5.034 ( -10.737 -6.199 11.461) 16.884 6.030 ( 22.440 12.956 -2.345) 26.017 10.019 ( 1.288 0.744 3.086) 3.426 11.246 ( -2.217 -1.280 -0.320) 2.580 11.441 ( -2.515 -1.452 -0.654) 2.977 ======================= 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| 2.085 ( -0.885 -0.511 4.979) 5.083 2.257 ( 2.049 1.183 3.470) 4.200 2.567 ( -2.387 -1.378 -8.718) 9.144 4.735 ( 5.791 3.344 3.644) 7.616 4.806 ( -8.024 -4.633 8.581) 12.629 6.472 ( 14.644 8.455 -1.531) 16.978 10.048 ( 1.092 0.630 3.411) 3.637 11.204 ( -1.291 -0.745 -0.496) 1.571 11.375 ( -2.737 -1.580 0.407) 3.187 ======================= 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.073 ( -0.186 -0.107 2.305) 2.315 2.289 ( 0.686 0.396 2.573) 2.692 2.537 ( -0.481 -0.278 -8.925) 8.943 4.676 ( -2.894 -1.671 6.721) 7.506 4.824 ( 1.833 1.058 4.702) 5.156 6.703 ( 4.962 2.865 -1.747) 5.990 10.066 ( 0.409 0.236 3.702) 3.732 11.185 ( -0.397 -0.229 -0.614) 0.766 11.327 ( -1.163 -0.671 1.223) 1.816 ======================= 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.951 ( 3.010 5.213 14.178) 15.403 1.982 ( 2.875 4.979 13.594) 14.760 3.172 ( -5.131 -8.887 18.945) 21.546 3.931 ( 7.083 12.268 -5.801) 15.307 4.344 ( 16.204 28.066 -2.826) 32.531 5.865 ( -6.934 -12.010 4.238) 14.501 9.977 ( 0.046 0.080 2.739) 2.740 11.386 ( -1.143 -1.980 0.078) 2.288 11.468 ( 0.506 0.877 -0.742) 1.255 ======================= 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| 2.047 ( 2.679 2.252 10.126) 10.714 2.084 ( 3.206 4.868 10.521) 12.028 2.955 ( -8.005 -9.658 8.176) 14.973 4.221 ( 8.991 15.246 -2.045) 17.818 4.963 ( 20.988 22.395 3.777) 30.923 5.605 ( -7.245 -9.665 2.966) 12.438 9.984 ( 0.490 0.343 2.828) 2.890 11.338 ( -2.276 -1.481 -0.106) 2.717 11.474 ( 0.101 -1.110 -1.199) 1.637 ======================= 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| 2.086 ( 1.136 -0.863 8.716) 8.832 2.183 ( 2.237 4.722 6.711) 8.505 2.729 ( -6.755 -7.558 -2.232) 10.379 4.520 ( 6.723 14.860 1.641) 16.392 5.157 ( -8.896 -4.058 11.075) 14.774 5.761 ( 21.984 13.158 -2.885) 25.783 10.002 ( 1.106 0.380 2.968) 3.190 11.281 ( -2.749 -0.834 -0.278) 2.886 11.445 ( -0.948 -3.193 -0.960) 3.466 ======================= 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| 2.078 ( 0.047 -1.967 6.363) 6.661 2.254 ( 1.030 3.437 3.757) 5.195 2.587 ( -2.811 -3.788 -7.662) 8.997 4.723 ( 2.056 6.535 6.359) 9.347 4.952 ( -9.338 2.473 5.290) 11.013 6.264 ( 19.169 9.851 -0.926) 21.572 10.027 ( 1.432 0.046 3.175) 3.484 11.228 ( -2.274 -0.605 -0.457) 2.397 11.384 ( -1.509 -4.031 -0.017) 4.304 ======================= 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| 2.063 ( 0.350 -1.357 3.193) 3.487 2.291 ( 0.256 1.809 1.909) 2.643 2.538 ( -0.468 -1.261 -8.316) 8.424 4.686 ( -2.610 -4.784 7.557) 9.317 4.947 ( -5.259 14.676 2.970) 15.870 6.601 ( 10.680 4.189 -0.922) 11.509 10.049 ( 1.212 -0.532 3.420) 3.667 11.191 ( -1.050 -0.877 -0.617) 1.500 11.324 ( -0.768 -3.165 1.044) 3.421 ======================= 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.061 ( 0.558 -0.967 1.917) 2.218 2.302 ( -0.465 0.806 1.335) 1.627 2.530 ( 0.234 -0.406 -8.121) 8.135 4.633 ( 1.203 -2.085 6.227) 6.676 4.979 ( -7.633 13.221 3.772) 15.726 6.717 ( 0.885 -1.534 -1.114) 2.092 10.057 ( 0.608 -1.053 3.546) 3.748 11.177 ( 0.412 -0.713 -0.659) 1.055 11.300 ( 0.990 -1.716 1.477) 2.471 ======================= 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| 2.070 ( 0.081 0.140 8.711) 8.712 2.190 ( 2.901 5.024 7.001) 9.092 2.750 ( -5.614 -9.724 -0.779) 11.255 4.551 ( 8.897 15.410 2.580) 17.980 5.233 ( -1.728 -2.993 9.502) 10.111 5.606 ( 9.967 17.264 -2.730) 20.121 9.993 ( 0.337 0.583 2.931) 3.008 11.307 ( -0.943 -1.633 -0.275) 1.906 11.430 ( -1.824 -3.160 -1.049) 3.797 ======================= 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| 2.058 ( -0.224 -1.697 6.818) 7.029 2.275 ( 1.343 3.934 3.616) 5.509 2.591 ( -3.150 -5.486 -6.429) 9.019 4.752 ( 0.646 2.926 9.675) 10.129 5.148 ( -4.134 7.894 -0.059) 8.911 6.044 ( 17.070 13.847 0.806) 21.995 10.007 ( 0.775 0.112 2.997) 3.097 11.267 ( -2.409 -0.699 -0.499) 2.557 11.356 ( -1.293 -5.250 -0.196) 5.410 ======================= 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| 2.041 ( 0.334 -1.405 3.926) 4.183 2.313 ( -0.276 2.029 0.231) 2.061 2.529 ( -0.301 -1.709 -6.855) 7.072 4.682 ( -3.115 -5.822 8.304) 10.609 5.234 ( -5.190 19.490 -0.448) 20.175 6.427 ( 14.085 6.328 1.284) 15.494 10.021 ( 1.053 -0.646 3.021) 3.264 11.206 ( -2.472 -2.321 -0.759) 3.474 11.296 ( -0.240 -3.720 1.042) 3.870 ======================= 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| 2.040 ( 0.610 -0.764 1.987) 2.214 2.318 ( -0.432 0.777 -1.493) 1.738 2.523 ( 0.185 -0.236 -6.087) 6.094 4.594 ( -0.406 -4.118 5.537) 6.913 5.309 ( -9.141 19.628 1.110) 21.680 6.643 ( 6.349 0.064 1.062) 6.437 10.029 ( 0.994 -1.295 3.049) 3.459 11.153 ( 0.018 -3.061 -0.575) 3.114 11.269 ( 0.300 -1.764 1.568) 2.379 ======================= 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| 2.030 ( -0.466 -0.807 4.163) 4.266 2.328 ( 0.656 1.137 -0.790) 1.532 2.521 ( -0.792 -1.371 -5.930) 6.138 4.679 ( -3.584 -6.208 8.562) 11.167 5.409 ( 7.702 13.340 -3.008) 15.695 6.311 ( 7.109 12.313 3.356) 14.609 10.007 ( -0.045 -0.078 2.849) 2.851 11.221 ( -2.814 -4.874 -1.165) 5.747 11.273 ( -0.896 -1.551 1.386) 2.265 ======================= 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| 2.027 ( 0.395 -0.024 2.116) 2.152 2.330 ( -0.455 -0.100 -4.045) 4.072 2.519 ( 0.339 0.399 -4.175) 4.208 4.569 ( -1.965 -4.796 5.086) 7.262 5.615 ( -1.825 16.497 -2.591) 16.799 6.525 ( 7.338 3.359 4.438) 9.210 10.005 ( 0.340 -0.827 2.640) 2.788 11.125 ( -1.666 -4.647 -0.370) 4.950 11.255 ( -0.145 -0.822 1.444) 1.667 ======================= 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.032 ( 0.059 -0.102 1.415) 1.420 2.324 ( 0.011 -0.020 -4.597) 4.597 2.525 ( -0.185 0.320 -3.980) 3.998 4.526 ( 1.421 -2.462 3.726) 4.686 5.682 ( -8.885 15.388 -2.287) 17.915 6.610 ( 1.660 -2.876 4.625) 5.694 10.003 ( 0.618 -1.071 2.550) 2.834 11.085 ( 1.707 -2.957 -0.193) 3.420 11.250 ( 0.252 -0.437 1.607) 1.685 ======================= 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.032 ( 0.161 0.280 1.258) 1.299 2.323 ( -0.174 -0.302 -5.608) 5.619 2.529 ( 0.216 0.374 -3.380) 3.407 4.500 ( -1.033 -1.789 3.031) 3.668 5.856 ( 3.682 6.377 -5.184) 9.005 6.568 ( 0.583 1.010 7.594) 7.683 9.992 ( -0.217 -0.375 2.326) 2.366 11.052 ( -1.167 -2.021 -0.078) 2.335 11.246 ( -0.133 -0.230 1.638) 1.660 ======================= 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| 2.124 ( 0.000 0.000 8.737) 8.737 2.124 ( 0.000 0.000 8.737) 8.737 3.504 ( 0.000 0.000 -5.210) 5.210 3.504 ( 0.000 0.000 -5.210) 5.210 3.795 ( -0.000 0.000 11.585) 11.585 6.144 ( 0.000 -0.000 -0.776) 0.776 10.020 ( 0.000 -0.000 1.098) 1.098 11.437 ( 0.000 0.000 0.109) 0.109 11.437 ( 0.000 0.000 0.109) 0.109 ======================= 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.139 ( 1.239 0.715 7.580) 7.714 2.149 ( 1.915 1.105 7.121) 7.456 3.619 ( 9.395 5.424 -4.208) 11.637 3.659 ( -8.627 -4.981 8.216) 12.913 3.819 ( 22.458 12.966 -1.767) 25.992 6.064 ( -6.659 -3.844 -0.033) 7.689 10.019 ( -0.074 -0.042 1.102) 1.105 11.419 ( -1.535 -0.886 0.087) 1.774 11.448 ( 0.802 0.463 -0.029) 0.927 ======================= 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.177 ( 1.890 1.091 5.255) 5.690 2.202 ( 2.429 1.402 5.196) 5.904 3.334 ( -17.195 -9.927 5.359) 20.565 3.911 ( 14.666 8.467 -2.202) 17.077 4.523 ( 34.027 19.645 -0.445) 39.293 5.845 ( -11.675 -6.741 1.229) 13.537 10.020 ( 0.272 0.157 1.095) 1.139 11.369 ( -2.592 -1.496 0.033) 2.993 11.466 ( 0.587 0.339 -0.263) 0.727 ======================= 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.224 ( 2.037 1.176 3.243) 4.006 2.255 ( 1.951 1.127 4.980) 5.466 2.905 ( -18.298 -10.564 0.078) 21.129 4.266 ( 14.813 8.552 -0.442) 17.110 5.300 ( 28.190 16.275 4.251) 32.827 5.552 ( -10.178 -5.876 -1.277) 11.821 10.034 ( 0.941 0.544 1.084) 1.535 11.303 ( -2.855 -1.648 -0.036) 3.297 11.467 ( -0.660 -0.381 -0.301) 0.819 ======================= 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.269 ( 1.803 1.041 1.949) 2.852 2.279 ( -0.614 -0.354 7.671) 7.704 2.533 ( -11.948 -6.898 -6.760) 15.363 4.578 ( 11.438 6.603 0.774) 13.230 5.192 ( -13.194 -7.618 3.974) 15.745 6.000 ( 24.211 13.978 -0.706) 27.965 10.063 ( 1.435 0.829 1.121) 2.001 11.241 ( -2.334 -1.347 -0.109) 2.697 11.435 ( -1.941 -1.121 -0.053) 2.242 ======================= 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.191 ( -5.257 -3.035 5.010) 7.871 2.306 ( 1.240 0.716 1.252) 1.902 2.401 ( -0.972 -0.561 -6.721) 6.814 4.791 ( 6.733 3.888 1.534) 7.925 4.919 ( -9.637 -5.564 2.651) 11.439 6.459 ( 14.709 8.492 -0.159) 16.985 10.096 ( 1.309 0.756 1.230) 1.948 11.197 ( -1.389 -0.802 -0.177) 1.614 11.385 ( -2.109 -1.218 0.331) 2.458 ======================= 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.109 ( -1.661 -0.959 1.128) 2.225 2.326 ( 0.442 0.255 0.963) 1.089 2.398 ( 0.143 0.082 -4.017) 4.021 4.761 ( -3.533 -2.040 1.925) 4.510 4.894 ( 2.154 1.244 1.910) 3.136 6.687 ( 4.846 2.798 -0.191) 5.599 10.119 ( 0.526 0.304 1.346) 1.477 11.176 ( -0.438 -0.253 -0.221) 0.552 11.347 ( -0.899 -0.519 0.602) 1.200 ======================= 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.164 ( 0.887 1.537 5.747) 6.015 2.189 ( 1.523 2.638 5.693) 6.457 3.435 ( -9.160 -15.866 6.472) 19.430 3.840 ( 8.772 15.193 -2.551) 17.728 4.300 ( 18.353 31.789 -1.142) 36.725 5.915 ( -6.125 -10.609 0.964) 12.288 10.018 ( 0.057 0.099 1.101) 1.107 11.387 ( -1.219 -2.111 0.037) 2.438 11.459 ( 0.279 0.483 -0.200) 0.592 ======================= 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.197 ( 1.448 0.777 3.953) 4.281 2.248 ( 1.772 2.966 4.687) 5.823 3.054 ( -14.031 -16.675 2.015) 21.885 4.191 ( 10.690 18.091 -0.801) 21.029 5.008 ( 26.024 26.479 0.982) 37.140 5.652 ( -10.014 -11.835 1.304) 15.558 10.025 ( 0.506 0.388 1.095) 1.267 11.336 ( -2.340 -1.659 -0.046) 2.869 11.460 ( 0.046 -1.106 -0.298) 1.147 ======================= 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.226 ( 1.662 -0.614 3.685) 4.089 2.303 ( 0.688 2.857 4.725) 5.565 2.655 ( -12.677 -13.032 -3.743) 18.562 4.545 ( 7.861 18.081 0.685) 19.728 5.316 ( -11.160 -10.525 4.172) 15.898 5.719 ( 24.679 17.806 -1.186) 30.455 10.045 ( 1.171 0.441 1.104) 1.669 11.276 ( -2.791 -1.023 -0.127) 2.975 11.435 ( -0.617 -2.784 -0.149) 2.855 ======================= 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.219 ( -1.675 -4.020 7.101) 8.330 2.305 ( -1.626 -0.358 1.022) 1.954 2.424 ( -2.740 -2.133 -7.533) 8.295 4.813 ( 2.899 15.143 2.412) 15.605 5.021 ( -9.957 -8.238 1.503) 13.010 6.260 ( 18.838 11.359 0.124) 21.998 10.073 ( 1.608 0.053 1.154) 1.980 11.221 ( -2.303 -0.784 -0.189) 2.440 11.388 ( -0.968 -3.419 0.204) 3.559 ======================= 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.123 ( -2.359 -3.281 2.333) 4.666 2.301 ( 1.211 -1.191 -1.031) 1.987 2.415 ( -0.455 1.651 -2.919) 3.384 4.791 ( -4.790 -6.330 2.660) 8.371 4.985 ( -2.644 14.056 0.802) 14.325 6.598 ( 9.889 5.210 0.166) 11.179 10.098 ( 1.470 -0.654 1.242) 2.032 11.182 ( -1.053 -0.958 -0.218) 1.440 11.342 ( -0.419 -2.707 0.542) 2.792 ======================= 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.087 ( 0.602 -1.043 0.739) 1.413 2.306 ( 0.969 -1.678 -0.712) 2.064 2.420 ( -0.897 1.553 -2.302) 2.919 4.713 ( 1.450 -2.512 1.818) 3.423 5.034 ( -6.801 11.778 1.438) 13.677 6.711 ( 0.254 -0.440 0.107) 0.519 10.108 ( 0.763 -1.321 1.293) 2.000 11.168 ( 0.402 -0.696 -0.219) 0.832 11.324 ( 0.918 -1.591 0.669) 1.955 ======================= 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.208 ( 0.156 0.270 3.746) 3.759 2.311 ( 1.673 2.898 4.339) 5.479 2.702 ( -9.746 -16.880 -2.746) 19.684 4.588 ( 11.295 19.563 0.992) 22.611 5.364 ( -7.169 -12.417 3.300) 14.713 5.571 ( 14.244 24.671 -0.822) 28.500 10.036 ( 0.363 0.630 1.110) 1.327 11.301 ( -1.055 -1.828 -0.185) 2.119 11.420 ( -1.575 -2.728 -0.136) 3.153 ======================= 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.200 ( -0.530 -2.328 6.126) 6.575 2.332 ( -2.254 -2.146 2.298) 3.869 2.444 ( -3.266 -4.580 -7.230) 9.161 4.922 ( 3.422 13.142 6.414) 15.019 5.103 ( -6.614 -5.202 -3.531) 9.125 6.070 ( 16.688 16.278 1.041) 23.336 10.050 ( 0.809 0.070 1.114) 1.379 11.257 ( -2.413 -1.159 -0.321) 2.696 11.360 ( -0.873 -4.253 0.243) 4.349 ======================= 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.125 ( -2.423 -4.349 4.048) 6.416 2.274 ( 0.048 -1.661 -3.680) 4.037 2.437 ( 0.042 1.714 -1.810) 2.493 4.804 ( -5.742 -9.221 3.343) 11.366 5.213 ( -1.790 19.375 -1.015) 19.484 6.459 ( 12.791 8.370 1.145) 15.328 10.064 ( 1.181 -0.854 1.108) 1.831 11.195 ( -2.252 -2.371 -0.282) 3.282 11.314 ( -0.038 -3.108 0.552) 3.157 ======================= 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| 2.070 ( -0.352 -1.794 0.928) 2.051 2.272 ( 0.977 -1.412 -2.114) 2.723 2.449 ( -0.518 1.435 -1.443) 2.100 4.664 ( -0.884 -5.786 1.593) 6.066 5.324 ( -7.220 17.857 0.315) 19.264 6.672 ( 4.978 2.170 1.054) 5.532 10.073 ( 1.216 -1.648 1.114) 2.331 11.146 ( 0.001 -2.770 -0.156) 2.775 11.293 ( 0.418 -1.700 0.669) 1.874 ======================= 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.124 ( -2.472 -4.281 4.841) 6.919 2.265 ( -1.072 -1.857 -4.822) 5.278 2.445 ( 0.797 1.381 -1.417) 2.133 4.807 ( -5.993 -10.380 3.567) 12.505 5.348 ( 8.591 14.880 -2.141) 17.315 6.374 ( 7.814 13.534 1.971) 15.752 10.048 ( -0.137 -0.237 1.055) 1.089 11.206 ( -2.392 -4.143 -0.383) 4.800 11.295 ( -0.933 -1.616 0.628) 1.968 ======================= 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| 2.063 ( -0.851 -1.799 1.296) 2.375 2.249 ( 0.251 -0.767 -2.938) 3.047 2.465 ( 0.056 0.912 -1.253) 1.551 4.633 ( -3.115 -7.089 1.458) 7.879 5.576 ( -0.166 15.028 -1.078) 15.067 6.602 ( 6.343 5.732 2.329) 8.861 10.042 ( 0.359 -1.115 0.968) 1.520 11.121 ( -1.535 -4.243 -0.080) 4.512 11.278 ( -0.034 -0.764 0.619) 0.984 ======================= 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| 2.048 ( 0.284 -0.492 0.325) 0.654 2.249 ( 0.456 -0.790 -2.233) 2.412 2.469 ( -0.334 0.578 -1.458) 1.604 4.568 ( 2.003 -3.470 0.817) 4.089 5.649 ( -7.376 12.775 -0.806) 14.774 6.691 ( 0.038 -0.066 2.397) 2.398 10.040 ( 0.788 -1.365 0.932) 1.831 11.084 ( 1.593 -2.760 0.006) 3.186 11.274 ( 0.254 -0.440 0.647) 0.823 ======================= 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| 2.044 ( -0.167 -0.290 0.176) 0.378 2.241 ( -0.079 -0.137 -2.214) 2.220 2.474 ( 0.081 0.140 -1.612) 1.620 4.531 ( -1.523 -2.637 0.534) 3.092 5.788 ( 3.142 5.442 -1.573) 6.477 6.686 ( 1.297 2.247 3.246) 4.155 10.025 ( -0.291 -0.504 0.849) 1.030 11.053 ( -1.112 -1.925 0.059) 2.224 11.270 ( -0.090 -0.156 0.642) 0.666 =================== 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 1156.911 1156.911 530.732 -0.000 0.000 0.000 3/16731 20.0 262.569 262.569 145.412 -0.000 0.000 0.000 3/16731 30.0 101.325 101.325 61.297 -0.000 0.000 0.000 3/16731 40.0 55.423 55.423 36.608 -0.000 0.000 0.000 3/16731 50.0 36.468 36.468 25.714 -0.000 0.000 0.000 3/16731 60.0 26.808 26.808 19.779 -0.000 0.000 0.000 3/16731 70.0 21.141 21.141 16.092 -0.000 0.000 0.000 3/16731 80.0 17.461 17.461 13.588 -0.000 0.000 0.000 3/16731 90.0 14.888 14.888 11.777 -0.000 0.000 0.000 3/16731 100.0 12.988 12.988 10.404 -0.000 0.000 0.000 3/16731 110.0 11.525 11.525 9.327 -0.000 0.000 0.000 3/16731 120.0 10.365 10.365 8.458 -0.000 0.000 0.000 3/16731 130.0 9.421 9.421 7.741 -0.000 0.000 0.000 3/16731 140.0 8.638 8.638 7.140 -0.000 0.000 0.000 3/16731 150.0 7.978 7.978 6.628 -0.000 0.000 0.000 3/16731 160.0 7.414 7.414 6.186 -0.000 0.000 0.000 3/16731 170.0 6.926 6.926 5.801 -0.000 0.000 0.000 3/16731 180.0 6.499 6.499 5.462 -0.000 0.000 0.000 3/16731 190.0 6.124 6.124 5.161 -0.000 0.000 0.000 3/16731 200.0 5.790 5.790 4.892 -0.000 0.000 0.000 3/16731 210.0 5.492 5.492 4.651 -0.000 0.000 0.000 3/16731 220.0 5.223 5.223 4.432 -0.000 0.000 0.000 3/16731 230.0 4.981 4.981 4.234 -0.000 0.000 0.000 3/16731 240.0 4.760 4.760 4.053 -0.000 0.000 0.000 3/16731 250.0 4.558 4.558 3.887 -0.000 0.000 0.000 3/16731 260.0 4.373 4.373 3.734 -0.000 0.000 0.000 3/16731 270.0 4.203 4.203 3.593 -0.000 0.000 0.000 3/16731 280.0 4.046 4.046 3.462 -0.000 0.000 0.000 3/16731 290.0 3.900 3.900 3.341 -0.000 0.000 0.000 3/16731 300.0 3.765 3.765 3.228 -0.000 0.000 0.000 3/16731 310.0 3.639 3.639 3.123 -0.000 0.000 0.000 3/16731 320.0 3.521 3.521 3.024 -0.000 0.000 0.000 3/16731 330.0 3.411 3.411 2.931 -0.000 0.000 0.000 3/16731 340.0 3.307 3.307 2.844 -0.000 0.000 0.000 3/16731 350.0 3.210 3.210 2.762 -0.000 0.000 0.000 3/16731 360.0 3.118 3.118 2.685 -0.000 0.000 0.000 3/16731 370.0 3.032 3.032 2.611 -0.000 0.000 0.000 3/16731 380.0 2.950 2.950 2.542 -0.000 0.000 0.000 3/16731 390.0 2.873 2.873 2.477 -0.000 0.000 0.000 3/16731 400.0 2.799 2.799 2.414 -0.000 0.000 0.000 3/16731 410.0 2.729 2.729 2.355 -0.000 0.000 0.000 3/16731 420.0 2.663 2.663 2.299 -0.000 0.000 0.000 3/16731 430.0 2.600 2.600 2.245 -0.000 0.000 0.000 3/16731 440.0 2.540 2.540 2.194 -0.000 0.000 0.000 3/16731 450.0 2.483 2.483 2.145 -0.000 0.000 0.000 3/16731 460.0 2.428 2.428 2.098 -0.000 0.000 0.000 3/16731 470.0 2.375 2.375 2.053 -0.000 0.000 0.000 3/16731 480.0 2.325 2.325 2.011 -0.000 0.000 0.000 3/16731 490.0 2.277 2.277 1.969 -0.000 0.000 0.000 3/16731 500.0 2.231 2.231 1.930 -0.000 0.000 0.000 3/16731 510.0 2.186 2.186 1.892 -0.000 0.000 0.000 3/16731 520.0 2.144 2.144 1.856 -0.000 0.000 0.000 3/16731 530.0 2.103 2.103 1.821 -0.000 0.000 0.000 3/16731 540.0 2.064 2.064 1.787 -0.000 0.000 0.000 3/16731 550.0 2.026 2.026 1.754 -0.000 0.000 0.000 3/16731 560.0 1.989 1.989 1.723 -0.000 0.000 0.000 3/16731 570.0 1.954 1.954 1.693 -0.000 0.000 0.000 3/16731 580.0 1.920 1.920 1.664 -0.000 0.000 0.000 3/16731 590.0 1.887 1.887 1.635 -0.000 0.000 0.000 3/16731 600.0 1.855 1.855 1.608 -0.000 0.000 0.000 3/16731 610.0 1.825 1.825 1.582 -0.000 0.000 0.000 3/16731 620.0 1.795 1.795 1.556 -0.000 0.000 0.000 3/16731 630.0 1.766 1.766 1.532 -0.000 0.000 0.000 3/16731 640.0 1.739 1.739 1.508 -0.000 0.000 0.000 3/16731 650.0 1.712 1.712 1.484 -0.000 0.000 0.000 3/16731 660.0 1.686 1.686 1.462 -0.000 0.000 0.000 3/16731 670.0 1.660 1.660 1.440 -0.000 0.000 0.000 3/16731 680.0 1.636 1.636 1.419 -0.000 0.000 0.000 3/16731 690.0 1.612 1.612 1.398 -0.000 0.000 0.000 3/16731 700.0 1.589 1.589 1.379 -0.000 0.000 0.000 3/16731 710.0 1.566 1.566 1.359 -0.000 0.000 0.000 3/16731 720.0 1.544 1.544 1.340 -0.000 0.000 0.000 3/16731 730.0 1.523 1.523 1.322 -0.000 0.000 0.000 3/16731 740.0 1.503 1.503 1.304 -0.000 0.000 0.000 3/16731 750.0 1.482 1.482 1.287 -0.000 0.000 0.000 3/16731 760.0 1.463 1.463 1.270 -0.000 0.000 0.000 3/16731 770.0 1.444 1.444 1.253 -0.000 0.000 0.000 3/16731 780.0 1.425 1.425 1.237 -0.000 0.000 0.000 3/16731 790.0 1.407 1.407 1.222 -0.000 0.000 0.000 3/16731 800.0 1.389 1.389 1.207 -0.000 0.000 0.000 3/16731 810.0 1.372 1.372 1.192 -0.000 0.000 0.000 3/16731 820.0 1.355 1.355 1.177 -0.000 0.000 0.000 3/16731 830.0 1.339 1.339 1.163 -0.000 0.000 0.000 3/16731 840.0 1.323 1.323 1.149 -0.000 0.000 0.000 3/16731 850.0 1.307 1.307 1.136 -0.000 0.000 0.000 3/16731 860.0 1.292 1.292 1.123 -0.000 0.000 0.000 3/16731 870.0 1.277 1.277 1.110 -0.000 0.000 0.000 3/16731 880.0 1.263 1.263 1.097 -0.000 0.000 0.000 3/16731 890.0 1.249 1.249 1.085 -0.000 0.000 0.000 3/16731 900.0 1.235 1.235 1.073 -0.000 0.000 0.000 3/16731 910.0 1.221 1.221 1.061 -0.000 0.000 0.000 3/16731 920.0 1.208 1.208 1.050 -0.000 0.000 0.000 3/16731 930.0 1.195 1.195 1.038 -0.000 0.000 0.000 3/16731 940.0 1.182 1.182 1.027 -0.000 0.000 0.000 3/16731 950.0 1.170 1.170 1.016 -0.000 0.000 0.000 3/16731 960.0 1.157 1.157 1.006 -0.000 0.000 0.000 3/16731 970.0 1.145 1.145 0.996 -0.000 0.000 0.000 3/16731 980.0 1.134 1.134 0.985 -0.000 0.000 0.000 3/16731 990.0 1.122 1.122 0.976 -0.000 0.000 0.000 3/16731 1000.0 1.111 1.111 0.966 -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:44:52]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|