----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 19:38:49]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) Number of symmetry operations in supercell: 1296 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.997036150000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.997036150000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.997036150000000 Atomic positions (fractional): *1 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 2 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 3 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 *4 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 *5 Ta 0.00000000000000 0.00000000000000 0.00000000000000 180.948 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.997036150000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.997036150000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.997036150000000 Atomic positions (fractional): *1 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 > 1 2 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 2 3 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 > 3 *4 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 4 *5 Ta 0.00000000000000 0.00000000000000 0.00000000000000 180.948 > 5 -------------------------------- super cell -------------------------------- Lattice vectors: a 11.991108450000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.991108450000000 0.000000000000000 c 0.000000000000000 0.000000000000000 11.991108450000000 Atomic positions (fractional): *1 O 0.00000000000000 0.16666666666667 0.00000000000000 15.999 > 1 2 O 0.33333333333333 0.16666666666667 0.00000000000000 15.999 > 1 3 O 0.66666666666667 0.16666666666667 0.00000000000000 15.999 > 1 4 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 > 1 5 O 0.33333333333333 0.50000000000000 0.00000000000000 15.999 > 1 6 O 0.66666666666667 0.50000000000000 0.00000000000000 15.999 > 1 7 O 0.00000000000000 0.83333333333333 0.00000000000000 15.999 > 1 8 O 0.33333333333333 0.83333333333333 0.00000000000000 15.999 > 1 9 O 0.66666666666667 0.83333333333333 0.00000000000000 15.999 > 1 10 O 0.00000000000000 0.16666666666667 0.33333333333333 15.999 > 1 11 O 0.33333333333333 0.16666666666667 0.33333333333333 15.999 > 1 12 O 0.66666666666667 0.16666666666667 0.33333333333333 15.999 > 1 13 O 0.00000000000000 0.50000000000000 0.33333333333333 15.999 > 1 14 O 0.33333333333333 0.50000000000000 0.33333333333333 15.999 > 1 15 O 0.66666666666667 0.50000000000000 0.33333333333333 15.999 > 1 16 O 0.00000000000000 0.83333333333333 0.33333333333333 15.999 > 1 17 O 0.33333333333333 0.83333333333333 0.33333333333333 15.999 > 1 18 O 0.66666666666667 0.83333333333333 0.33333333333333 15.999 > 1 19 O 0.00000000000000 0.16666666666667 0.66666666666667 15.999 > 1 20 O 0.33333333333333 0.16666666666667 0.66666666666667 15.999 > 1 21 O 0.66666666666667 0.16666666666667 0.66666666666667 15.999 > 1 22 O 0.00000000000000 0.50000000000000 0.66666666666667 15.999 > 1 23 O 0.33333333333333 0.50000000000000 0.66666666666667 15.999 > 1 24 O 0.66666666666667 0.50000000000000 0.66666666666667 15.999 > 1 25 O 0.00000000000000 0.83333333333333 0.66666666666667 15.999 > 1 26 O 0.33333333333333 0.83333333333333 0.66666666666667 15.999 > 1 27 O 0.66666666666667 0.83333333333333 0.66666666666667 15.999 > 1 28 O 0.00000000000000 0.00000000000000 0.16666666666667 15.999 > 2 29 O 0.33333333333333 0.00000000000000 0.16666666666667 15.999 > 2 30 O 0.66666666666667 0.00000000000000 0.16666666666667 15.999 > 2 31 O 0.00000000000000 0.33333333333333 0.16666666666667 15.999 > 2 32 O 0.33333333333333 0.33333333333333 0.16666666666667 15.999 > 2 33 O 0.66666666666667 0.33333333333333 0.16666666666667 15.999 > 2 34 O 0.00000000000000 0.66666666666667 0.16666666666667 15.999 > 2 35 O 0.33333333333333 0.66666666666667 0.16666666666667 15.999 > 2 36 O 0.66666666666667 0.66666666666667 0.16666666666667 15.999 > 2 37 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 2 38 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 2 39 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 2 40 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 2 41 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 2 42 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 2 43 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 2 44 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 2 45 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 2 46 O 0.00000000000000 0.00000000000000 0.83333333333333 15.999 > 2 47 O 0.33333333333333 0.00000000000000 0.83333333333333 15.999 > 2 48 O 0.66666666666667 0.00000000000000 0.83333333333333 15.999 > 2 49 O 0.00000000000000 0.33333333333333 0.83333333333333 15.999 > 2 50 O 0.33333333333333 0.33333333333333 0.83333333333333 15.999 > 2 51 O 0.66666666666667 0.33333333333333 0.83333333333333 15.999 > 2 52 O 0.00000000000000 0.66666666666667 0.83333333333333 15.999 > 2 53 O 0.33333333333333 0.66666666666667 0.83333333333333 15.999 > 2 54 O 0.66666666666667 0.66666666666667 0.83333333333333 15.999 > 2 55 O 0.16666666666667 0.00000000000000 0.00000000000000 15.999 > 3 56 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 > 3 57 O 0.83333333333333 0.00000000000000 0.00000000000000 15.999 > 3 58 O 0.16666666666667 0.33333333333333 0.00000000000000 15.999 > 3 59 O 0.50000000000000 0.33333333333333 0.00000000000000 15.999 > 3 60 O 0.83333333333333 0.33333333333333 0.00000000000000 15.999 > 3 61 O 0.16666666666667 0.66666666666667 0.00000000000000 15.999 > 3 62 O 0.50000000000000 0.66666666666667 0.00000000000000 15.999 > 3 63 O 0.83333333333333 0.66666666666667 0.00000000000000 15.999 > 3 64 O 0.16666666666667 0.00000000000000 0.33333333333333 15.999 > 3 65 O 0.50000000000000 0.00000000000000 0.33333333333333 15.999 > 3 66 O 0.83333333333333 0.00000000000000 0.33333333333333 15.999 > 3 67 O 0.16666666666667 0.33333333333333 0.33333333333333 15.999 > 3 68 O 0.50000000000000 0.33333333333333 0.33333333333333 15.999 > 3 69 O 0.83333333333333 0.33333333333333 0.33333333333333 15.999 > 3 70 O 0.16666666666667 0.66666666666667 0.33333333333333 15.999 > 3 71 O 0.50000000000000 0.66666666666667 0.33333333333333 15.999 > 3 72 O 0.83333333333333 0.66666666666667 0.33333333333333 15.999 > 3 73 O 0.16666666666667 0.00000000000000 0.66666666666667 15.999 > 3 74 O 0.50000000000000 0.00000000000000 0.66666666666667 15.999 > 3 75 O 0.83333333333333 0.00000000000000 0.66666666666667 15.999 > 3 76 O 0.16666666666667 0.33333333333333 0.66666666666667 15.999 > 3 77 O 0.50000000000000 0.33333333333333 0.66666666666667 15.999 > 3 78 O 0.83333333333333 0.33333333333333 0.66666666666667 15.999 > 3 79 O 0.16666666666667 0.66666666666667 0.66666666666667 15.999 > 3 80 O 0.50000000000000 0.66666666666667 0.66666666666667 15.999 > 3 81 O 0.83333333333333 0.66666666666667 0.66666666666667 15.999 > 3 *82 K 0.16666666666667 0.16666666666667 0.16666666666667 39.098 > 4 83 K 0.50000000000000 0.16666666666667 0.16666666666667 39.098 > 4 84 K 0.83333333333333 0.16666666666667 0.16666666666667 39.098 > 4 85 K 0.16666666666667 0.50000000000000 0.16666666666667 39.098 > 4 86 K 0.50000000000000 0.50000000000000 0.16666666666667 39.098 > 4 87 K 0.83333333333333 0.50000000000000 0.16666666666667 39.098 > 4 88 K 0.16666666666667 0.83333333333333 0.16666666666667 39.098 > 4 89 K 0.50000000000000 0.83333333333333 0.16666666666667 39.098 > 4 90 K 0.83333333333333 0.83333333333333 0.16666666666667 39.098 > 4 91 K 0.16666666666667 0.16666666666667 0.50000000000000 39.098 > 4 92 K 0.50000000000000 0.16666666666667 0.50000000000000 39.098 > 4 93 K 0.83333333333333 0.16666666666667 0.50000000000000 39.098 > 4 94 K 0.16666666666667 0.50000000000000 0.50000000000000 39.098 > 4 95 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 4 96 K 0.83333333333333 0.50000000000000 0.50000000000000 39.098 > 4 97 K 0.16666666666667 0.83333333333333 0.50000000000000 39.098 > 4 98 K 0.50000000000000 0.83333333333333 0.50000000000000 39.098 > 4 99 K 0.83333333333333 0.83333333333333 0.50000000000000 39.098 > 4 100 K 0.16666666666667 0.16666666666667 0.83333333333333 39.098 > 4 101 K 0.50000000000000 0.16666666666667 0.83333333333333 39.098 > 4 102 K 0.83333333333333 0.16666666666667 0.83333333333333 39.098 > 4 103 K 0.16666666666667 0.50000000000000 0.83333333333333 39.098 > 4 104 K 0.50000000000000 0.50000000000000 0.83333333333333 39.098 > 4 105 K 0.83333333333333 0.50000000000000 0.83333333333333 39.098 > 4 106 K 0.16666666666667 0.83333333333333 0.83333333333333 39.098 > 4 107 K 0.50000000000000 0.83333333333333 0.83333333333333 39.098 > 4 108 K 0.83333333333333 0.83333333333333 0.83333333333333 39.098 > 4 *109 Ta 0.00000000000000 0.00000000000000 0.00000000000000 180.948 > 5 110 Ta 0.33333333333333 0.00000000000000 0.00000000000000 180.948 > 5 111 Ta 0.66666666666667 0.00000000000000 0.00000000000000 180.948 > 5 112 Ta 0.00000000000000 0.33333333333333 0.00000000000000 180.948 > 5 113 Ta 0.33333333333333 0.33333333333333 0.00000000000000 180.948 > 5 114 Ta 0.66666666666667 0.33333333333333 0.00000000000000 180.948 > 5 115 Ta 0.00000000000000 0.66666666666667 0.00000000000000 180.948 > 5 116 Ta 0.33333333333333 0.66666666666667 0.00000000000000 180.948 > 5 117 Ta 0.66666666666667 0.66666666666667 0.00000000000000 180.948 > 5 118 Ta 0.00000000000000 0.00000000000000 0.33333333333333 180.948 > 5 119 Ta 0.33333333333333 0.00000000000000 0.33333333333333 180.948 > 5 120 Ta 0.66666666666667 0.00000000000000 0.33333333333333 180.948 > 5 121 Ta 0.00000000000000 0.33333333333333 0.33333333333333 180.948 > 5 122 Ta 0.33333333333333 0.33333333333333 0.33333333333333 180.948 > 5 123 Ta 0.66666666666667 0.33333333333333 0.33333333333333 180.948 > 5 124 Ta 0.00000000000000 0.66666666666667 0.33333333333333 180.948 > 5 125 Ta 0.33333333333333 0.66666666666667 0.33333333333333 180.948 > 5 126 Ta 0.66666666666667 0.66666666666667 0.33333333333333 180.948 > 5 127 Ta 0.00000000000000 0.00000000000000 0.66666666666667 180.948 > 5 128 Ta 0.33333333333333 0.00000000000000 0.66666666666667 180.948 > 5 129 Ta 0.66666666666667 0.00000000000000 0.66666666666667 180.948 > 5 130 Ta 0.00000000000000 0.33333333333333 0.66666666666667 180.948 > 5 131 Ta 0.33333333333333 0.33333333333333 0.66666666666667 180.948 > 5 132 Ta 0.66666666666667 0.33333333333333 0.66666666666667 180.948 > 5 133 Ta 0.00000000000000 0.66666666666667 0.66666666666667 180.948 > 5 134 Ta 0.33333333333333 0.66666666666667 0.66666666666667 180.948 > 5 135 Ta 0.66666666666667 0.66666666666667 0.66666666666667 180.948 > 5 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.7507916 0.0000000 0.0000000 0.0000000 5.7507916 0.0000000 0.0000000 0.0000000 5.7507916 -------------------------- Born effective charges -------------------------- 1 O -1.7009814 0.0000000 0.0000000 0.0000000 -6.7657002 0.0000000 0.0000000 0.0000000 -1.7009814 2 O -1.7009814 0.0000000 0.0000000 0.0000000 -1.7009814 0.0000000 0.0000000 0.0000000 -6.7657002 3 O -6.7657002 0.0000000 0.0000000 0.0000000 -1.7009814 0.0000000 0.0000000 0.0000000 -1.7009814 4 K 1.1716213 0.0000000 0.0000000 0.0000000 1.1716213 0.0000000 0.0000000 0.0000000 1.1716213 5 Ta 8.9960418 0.0000000 0.0000000 0.0000000 8.9960418 0.0000000 0.0000000 0.0000000 8.9960418 ---------------------------------------------------------------------------- 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: 405/405 Permutation basis: 3615/3615 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 104 Number of blocks in projector: 104 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 74 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 30 Use standard eigh solver. Tree of FC basis block matrices: - (104, 100), data: False |-- (30, 30), data: True |-- (74, 70), data: True ----- Solver_atoms: 1 -- 135 / 135 Time (Solver_compr_matrix_reshape): 0.004 Solver_block: 80 / 80 - Time: 0.054 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.058 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 405/405 Permutation basis: 3615/3615 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 104 Number of blocks in projector: 104 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 74 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 30 Use standard eigh solver. Tree of FC basis block matrices: - (104, 100), data: False |-- (30, 30), data: True |-- (74, 70), data: True Max drift after symmetrization by symfc projector: 0.00000000 (zz) 0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 19:38:54]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 19:38:54]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.997036150000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.997036150000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.997036150000000 Atomic positions (fractional): 1 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 2 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 3 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 4 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 5 Ta 0.00000000000000 0.00000000000000 0.00000000000000 180.948 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.991108450000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.991108450000000 0.000000000000000 c 0.000000000000000 0.000000000000000 11.991108450000000 Atomic positions (fractional): 1 O 0.00000000000000 0.16666666666667 0.00000000000000 15.999 > 1 2 O 0.33333333333333 0.16666666666667 0.00000000000000 15.999 > 1 3 O 0.66666666666667 0.16666666666667 0.00000000000000 15.999 > 1 4 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 > 1 5 O 0.33333333333333 0.50000000000000 0.00000000000000 15.999 > 1 6 O 0.66666666666667 0.50000000000000 0.00000000000000 15.999 > 1 7 O 0.00000000000000 0.83333333333333 0.00000000000000 15.999 > 1 8 O 0.33333333333333 0.83333333333333 0.00000000000000 15.999 > 1 9 O 0.66666666666667 0.83333333333333 0.00000000000000 15.999 > 1 10 O 0.00000000000000 0.16666666666667 0.33333333333333 15.999 > 1 11 O 0.33333333333333 0.16666666666667 0.33333333333333 15.999 > 1 12 O 0.66666666666667 0.16666666666667 0.33333333333333 15.999 > 1 13 O 0.00000000000000 0.50000000000000 0.33333333333333 15.999 > 1 14 O 0.33333333333333 0.50000000000000 0.33333333333333 15.999 > 1 15 O 0.66666666666667 0.50000000000000 0.33333333333333 15.999 > 1 16 O 0.00000000000000 0.83333333333333 0.33333333333333 15.999 > 1 17 O 0.33333333333333 0.83333333333333 0.33333333333333 15.999 > 1 18 O 0.66666666666667 0.83333333333333 0.33333333333333 15.999 > 1 19 O 0.00000000000000 0.16666666666667 0.66666666666667 15.999 > 1 20 O 0.33333333333333 0.16666666666667 0.66666666666667 15.999 > 1 21 O 0.66666666666667 0.16666666666667 0.66666666666667 15.999 > 1 22 O 0.00000000000000 0.50000000000000 0.66666666666667 15.999 > 1 23 O 0.33333333333333 0.50000000000000 0.66666666666667 15.999 > 1 24 O 0.66666666666667 0.50000000000000 0.66666666666667 15.999 > 1 25 O 0.00000000000000 0.83333333333333 0.66666666666667 15.999 > 1 26 O 0.33333333333333 0.83333333333333 0.66666666666667 15.999 > 1 27 O 0.66666666666667 0.83333333333333 0.66666666666667 15.999 > 1 28 O 0.00000000000000 0.00000000000000 0.16666666666667 15.999 > 28 29 O 0.33333333333333 0.00000000000000 0.16666666666667 15.999 > 28 30 O 0.66666666666667 0.00000000000000 0.16666666666667 15.999 > 28 31 O 0.00000000000000 0.33333333333333 0.16666666666667 15.999 > 28 32 O 0.33333333333333 0.33333333333333 0.16666666666667 15.999 > 28 33 O 0.66666666666667 0.33333333333333 0.16666666666667 15.999 > 28 34 O 0.00000000000000 0.66666666666667 0.16666666666667 15.999 > 28 35 O 0.33333333333333 0.66666666666667 0.16666666666667 15.999 > 28 36 O 0.66666666666667 0.66666666666667 0.16666666666667 15.999 > 28 37 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 28 38 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 28 39 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 28 40 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 28 41 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 28 42 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 28 43 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 28 44 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 28 45 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 28 46 O 0.00000000000000 0.00000000000000 0.83333333333333 15.999 > 28 47 O 0.33333333333333 0.00000000000000 0.83333333333333 15.999 > 28 48 O 0.66666666666667 0.00000000000000 0.83333333333333 15.999 > 28 49 O 0.00000000000000 0.33333333333333 0.83333333333333 15.999 > 28 50 O 0.33333333333333 0.33333333333333 0.83333333333333 15.999 > 28 51 O 0.66666666666667 0.33333333333333 0.83333333333333 15.999 > 28 52 O 0.00000000000000 0.66666666666667 0.83333333333333 15.999 > 28 53 O 0.33333333333333 0.66666666666667 0.83333333333333 15.999 > 28 54 O 0.66666666666667 0.66666666666667 0.83333333333333 15.999 > 28 55 O 0.16666666666667 0.00000000000000 0.00000000000000 15.999 > 55 56 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 > 55 57 O 0.83333333333333 0.00000000000000 0.00000000000000 15.999 > 55 58 O 0.16666666666667 0.33333333333333 0.00000000000000 15.999 > 55 59 O 0.50000000000000 0.33333333333333 0.00000000000000 15.999 > 55 60 O 0.83333333333333 0.33333333333333 0.00000000000000 15.999 > 55 61 O 0.16666666666667 0.66666666666667 0.00000000000000 15.999 > 55 62 O 0.50000000000000 0.66666666666667 0.00000000000000 15.999 > 55 63 O 0.83333333333333 0.66666666666667 0.00000000000000 15.999 > 55 64 O 0.16666666666667 0.00000000000000 0.33333333333333 15.999 > 55 65 O 0.50000000000000 0.00000000000000 0.33333333333333 15.999 > 55 66 O 0.83333333333333 0.00000000000000 0.33333333333333 15.999 > 55 67 O 0.16666666666667 0.33333333333333 0.33333333333333 15.999 > 55 68 O 0.50000000000000 0.33333333333333 0.33333333333333 15.999 > 55 69 O 0.83333333333333 0.33333333333333 0.33333333333333 15.999 > 55 70 O 0.16666666666667 0.66666666666667 0.33333333333333 15.999 > 55 71 O 0.50000000000000 0.66666666666667 0.33333333333333 15.999 > 55 72 O 0.83333333333333 0.66666666666667 0.33333333333333 15.999 > 55 73 O 0.16666666666667 0.00000000000000 0.66666666666667 15.999 > 55 74 O 0.50000000000000 0.00000000000000 0.66666666666667 15.999 > 55 75 O 0.83333333333333 0.00000000000000 0.66666666666667 15.999 > 55 76 O 0.16666666666667 0.33333333333333 0.66666666666667 15.999 > 55 77 O 0.50000000000000 0.33333333333333 0.66666666666667 15.999 > 55 78 O 0.83333333333333 0.33333333333333 0.66666666666667 15.999 > 55 79 O 0.16666666666667 0.66666666666667 0.66666666666667 15.999 > 55 80 O 0.50000000000000 0.66666666666667 0.66666666666667 15.999 > 55 81 O 0.83333333333333 0.66666666666667 0.66666666666667 15.999 > 55 82 K 0.16666666666667 0.16666666666667 0.16666666666667 39.098 > 82 83 K 0.50000000000000 0.16666666666667 0.16666666666667 39.098 > 82 84 K 0.83333333333333 0.16666666666667 0.16666666666667 39.098 > 82 85 K 0.16666666666667 0.50000000000000 0.16666666666667 39.098 > 82 86 K 0.50000000000000 0.50000000000000 0.16666666666667 39.098 > 82 87 K 0.83333333333333 0.50000000000000 0.16666666666667 39.098 > 82 88 K 0.16666666666667 0.83333333333333 0.16666666666667 39.098 > 82 89 K 0.50000000000000 0.83333333333333 0.16666666666667 39.098 > 82 90 K 0.83333333333333 0.83333333333333 0.16666666666667 39.098 > 82 91 K 0.16666666666667 0.16666666666667 0.50000000000000 39.098 > 82 92 K 0.50000000000000 0.16666666666667 0.50000000000000 39.098 > 82 93 K 0.83333333333333 0.16666666666667 0.50000000000000 39.098 > 82 94 K 0.16666666666667 0.50000000000000 0.50000000000000 39.098 > 82 95 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 82 96 K 0.83333333333333 0.50000000000000 0.50000000000000 39.098 > 82 97 K 0.16666666666667 0.83333333333333 0.50000000000000 39.098 > 82 98 K 0.50000000000000 0.83333333333333 0.50000000000000 39.098 > 82 99 K 0.83333333333333 0.83333333333333 0.50000000000000 39.098 > 82 100 K 0.16666666666667 0.16666666666667 0.83333333333333 39.098 > 82 101 K 0.50000000000000 0.16666666666667 0.83333333333333 39.098 > 82 102 K 0.83333333333333 0.16666666666667 0.83333333333333 39.098 > 82 103 K 0.16666666666667 0.50000000000000 0.83333333333333 39.098 > 82 104 K 0.50000000000000 0.50000000000000 0.83333333333333 39.098 > 82 105 K 0.83333333333333 0.50000000000000 0.83333333333333 39.098 > 82 106 K 0.16666666666667 0.83333333333333 0.83333333333333 39.098 > 82 107 K 0.50000000000000 0.83333333333333 0.83333333333333 39.098 > 82 108 K 0.83333333333333 0.83333333333333 0.83333333333333 39.098 > 82 109 Ta 0.00000000000000 0.00000000000000 0.00000000000000 180.948 > 109 110 Ta 0.33333333333333 0.00000000000000 0.00000000000000 180.948 > 109 111 Ta 0.66666666666667 0.00000000000000 0.00000000000000 180.948 > 109 112 Ta 0.00000000000000 0.33333333333333 0.00000000000000 180.948 > 109 113 Ta 0.33333333333333 0.33333333333333 0.00000000000000 180.948 > 109 114 Ta 0.66666666666667 0.33333333333333 0.00000000000000 180.948 > 109 115 Ta 0.00000000000000 0.66666666666667 0.00000000000000 180.948 > 109 116 Ta 0.33333333333333 0.66666666666667 0.00000000000000 180.948 > 109 117 Ta 0.66666666666667 0.66666666666667 0.00000000000000 180.948 > 109 118 Ta 0.00000000000000 0.00000000000000 0.33333333333333 180.948 > 109 119 Ta 0.33333333333333 0.00000000000000 0.33333333333333 180.948 > 109 120 Ta 0.66666666666667 0.00000000000000 0.33333333333333 180.948 > 109 121 Ta 0.00000000000000 0.33333333333333 0.33333333333333 180.948 > 109 122 Ta 0.33333333333333 0.33333333333333 0.33333333333333 180.948 > 109 123 Ta 0.66666666666667 0.33333333333333 0.33333333333333 180.948 > 109 124 Ta 0.00000000000000 0.66666666666667 0.33333333333333 180.948 > 109 125 Ta 0.33333333333333 0.66666666666667 0.33333333333333 180.948 > 109 126 Ta 0.66666666666667 0.66666666666667 0.33333333333333 180.948 > 109 127 Ta 0.00000000000000 0.00000000000000 0.66666666666667 180.948 > 109 128 Ta 0.33333333333333 0.00000000000000 0.66666666666667 180.948 > 109 129 Ta 0.66666666666667 0.00000000000000 0.66666666666667 180.948 > 109 130 Ta 0.00000000000000 0.33333333333333 0.66666666666667 180.948 > 109 131 Ta 0.33333333333333 0.33333333333333 0.66666666666667 180.948 > 109 132 Ta 0.66666666666667 0.33333333333333 0.66666666666667 180.948 > 109 133 Ta 0.00000000000000 0.66666666666667 0.66666666666667 180.948 > 109 134 Ta 0.33333333333333 0.66666666666667 0.66666666666667 180.948 > 109 135 Ta 0.66666666666667 0.66666666666667 0.66666666666667 180.948 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.7507916 0.0000000 0.0000000 0.0000000 5.7507916 0.0000000 0.0000000 0.0000000 5.7507916 -------------------------- Born effective charges -------------------------- 1 O -1.7009814 0.0000000 0.0000000 0.0000000 -6.7657002 0.0000000 0.0000000 0.0000000 -1.7009814 2 O -1.7009814 0.0000000 0.0000000 0.0000000 -1.7009814 0.0000000 0.0000000 0.0000000 -6.7657002 3 O -6.7657002 0.0000000 0.0000000 0.0000000 -1.7009814 0.0000000 0.0000000 0.0000000 -1.7009814 4 K 1.1716213 0.0000000 0.0000000 0.0000000 1.1716213 0.0000000 0.0000000 0.0000000 1.1716213 5 Ta 8.9960418 0.0000000 0.0000000 0.0000000 8.9960418 0.0000000 0.0000000 0.0000000 8.9960418 ---------------------------------------------------------------------------- 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.0100 0.0000] [ 0.0000 -0.0100 0.0000] Computing fc3[ 82, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 109, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 19:38:58]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 19:38:58]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.997036150000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.997036150000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.997036150000000 Atomic positions (fractional): 1 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 2 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 3 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 4 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 5 Ta 0.00000000000000 0.00000000000000 0.00000000000000 180.948 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.991108450000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.991108450000000 0.000000000000000 c 0.000000000000000 0.000000000000000 11.991108450000000 Atomic positions (fractional): 1 O 0.00000000000000 0.16666666666667 0.00000000000000 15.999 > 1 2 O 0.33333333333333 0.16666666666667 0.00000000000000 15.999 > 1 3 O 0.66666666666667 0.16666666666667 0.00000000000000 15.999 > 1 4 O 0.00000000000000 0.50000000000000 0.00000000000000 15.999 > 1 5 O 0.33333333333333 0.50000000000000 0.00000000000000 15.999 > 1 6 O 0.66666666666667 0.50000000000000 0.00000000000000 15.999 > 1 7 O 0.00000000000000 0.83333333333333 0.00000000000000 15.999 > 1 8 O 0.33333333333333 0.83333333333333 0.00000000000000 15.999 > 1 9 O 0.66666666666667 0.83333333333333 0.00000000000000 15.999 > 1 10 O 0.00000000000000 0.16666666666667 0.33333333333333 15.999 > 1 11 O 0.33333333333333 0.16666666666667 0.33333333333333 15.999 > 1 12 O 0.66666666666667 0.16666666666667 0.33333333333333 15.999 > 1 13 O 0.00000000000000 0.50000000000000 0.33333333333333 15.999 > 1 14 O 0.33333333333333 0.50000000000000 0.33333333333333 15.999 > 1 15 O 0.66666666666667 0.50000000000000 0.33333333333333 15.999 > 1 16 O 0.00000000000000 0.83333333333333 0.33333333333333 15.999 > 1 17 O 0.33333333333333 0.83333333333333 0.33333333333333 15.999 > 1 18 O 0.66666666666667 0.83333333333333 0.33333333333333 15.999 > 1 19 O 0.00000000000000 0.16666666666667 0.66666666666667 15.999 > 1 20 O 0.33333333333333 0.16666666666667 0.66666666666667 15.999 > 1 21 O 0.66666666666667 0.16666666666667 0.66666666666667 15.999 > 1 22 O 0.00000000000000 0.50000000000000 0.66666666666667 15.999 > 1 23 O 0.33333333333333 0.50000000000000 0.66666666666667 15.999 > 1 24 O 0.66666666666667 0.50000000000000 0.66666666666667 15.999 > 1 25 O 0.00000000000000 0.83333333333333 0.66666666666667 15.999 > 1 26 O 0.33333333333333 0.83333333333333 0.66666666666667 15.999 > 1 27 O 0.66666666666667 0.83333333333333 0.66666666666667 15.999 > 1 28 O 0.00000000000000 0.00000000000000 0.16666666666667 15.999 > 28 29 O 0.33333333333333 0.00000000000000 0.16666666666667 15.999 > 28 30 O 0.66666666666667 0.00000000000000 0.16666666666667 15.999 > 28 31 O 0.00000000000000 0.33333333333333 0.16666666666667 15.999 > 28 32 O 0.33333333333333 0.33333333333333 0.16666666666667 15.999 > 28 33 O 0.66666666666667 0.33333333333333 0.16666666666667 15.999 > 28 34 O 0.00000000000000 0.66666666666667 0.16666666666667 15.999 > 28 35 O 0.33333333333333 0.66666666666667 0.16666666666667 15.999 > 28 36 O 0.66666666666667 0.66666666666667 0.16666666666667 15.999 > 28 37 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 28 38 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 28 39 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 28 40 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 28 41 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 28 42 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 28 43 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 28 44 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 28 45 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 28 46 O 0.00000000000000 0.00000000000000 0.83333333333333 15.999 > 28 47 O 0.33333333333333 0.00000000000000 0.83333333333333 15.999 > 28 48 O 0.66666666666667 0.00000000000000 0.83333333333333 15.999 > 28 49 O 0.00000000000000 0.33333333333333 0.83333333333333 15.999 > 28 50 O 0.33333333333333 0.33333333333333 0.83333333333333 15.999 > 28 51 O 0.66666666666667 0.33333333333333 0.83333333333333 15.999 > 28 52 O 0.00000000000000 0.66666666666667 0.83333333333333 15.999 > 28 53 O 0.33333333333333 0.66666666666667 0.83333333333333 15.999 > 28 54 O 0.66666666666667 0.66666666666667 0.83333333333333 15.999 > 28 55 O 0.16666666666667 0.00000000000000 0.00000000000000 15.999 > 55 56 O 0.50000000000000 0.00000000000000 0.00000000000000 15.999 > 55 57 O 0.83333333333333 0.00000000000000 0.00000000000000 15.999 > 55 58 O 0.16666666666667 0.33333333333333 0.00000000000000 15.999 > 55 59 O 0.50000000000000 0.33333333333333 0.00000000000000 15.999 > 55 60 O 0.83333333333333 0.33333333333333 0.00000000000000 15.999 > 55 61 O 0.16666666666667 0.66666666666667 0.00000000000000 15.999 > 55 62 O 0.50000000000000 0.66666666666667 0.00000000000000 15.999 > 55 63 O 0.83333333333333 0.66666666666667 0.00000000000000 15.999 > 55 64 O 0.16666666666667 0.00000000000000 0.33333333333333 15.999 > 55 65 O 0.50000000000000 0.00000000000000 0.33333333333333 15.999 > 55 66 O 0.83333333333333 0.00000000000000 0.33333333333333 15.999 > 55 67 O 0.16666666666667 0.33333333333333 0.33333333333333 15.999 > 55 68 O 0.50000000000000 0.33333333333333 0.33333333333333 15.999 > 55 69 O 0.83333333333333 0.33333333333333 0.33333333333333 15.999 > 55 70 O 0.16666666666667 0.66666666666667 0.33333333333333 15.999 > 55 71 O 0.50000000000000 0.66666666666667 0.33333333333333 15.999 > 55 72 O 0.83333333333333 0.66666666666667 0.33333333333333 15.999 > 55 73 O 0.16666666666667 0.00000000000000 0.66666666666667 15.999 > 55 74 O 0.50000000000000 0.00000000000000 0.66666666666667 15.999 > 55 75 O 0.83333333333333 0.00000000000000 0.66666666666667 15.999 > 55 76 O 0.16666666666667 0.33333333333333 0.66666666666667 15.999 > 55 77 O 0.50000000000000 0.33333333333333 0.66666666666667 15.999 > 55 78 O 0.83333333333333 0.33333333333333 0.66666666666667 15.999 > 55 79 O 0.16666666666667 0.66666666666667 0.66666666666667 15.999 > 55 80 O 0.50000000000000 0.66666666666667 0.66666666666667 15.999 > 55 81 O 0.83333333333333 0.66666666666667 0.66666666666667 15.999 > 55 82 K 0.16666666666667 0.16666666666667 0.16666666666667 39.098 > 82 83 K 0.50000000000000 0.16666666666667 0.16666666666667 39.098 > 82 84 K 0.83333333333333 0.16666666666667 0.16666666666667 39.098 > 82 85 K 0.16666666666667 0.50000000000000 0.16666666666667 39.098 > 82 86 K 0.50000000000000 0.50000000000000 0.16666666666667 39.098 > 82 87 K 0.83333333333333 0.50000000000000 0.16666666666667 39.098 > 82 88 K 0.16666666666667 0.83333333333333 0.16666666666667 39.098 > 82 89 K 0.50000000000000 0.83333333333333 0.16666666666667 39.098 > 82 90 K 0.83333333333333 0.83333333333333 0.16666666666667 39.098 > 82 91 K 0.16666666666667 0.16666666666667 0.50000000000000 39.098 > 82 92 K 0.50000000000000 0.16666666666667 0.50000000000000 39.098 > 82 93 K 0.83333333333333 0.16666666666667 0.50000000000000 39.098 > 82 94 K 0.16666666666667 0.50000000000000 0.50000000000000 39.098 > 82 95 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 82 96 K 0.83333333333333 0.50000000000000 0.50000000000000 39.098 > 82 97 K 0.16666666666667 0.83333333333333 0.50000000000000 39.098 > 82 98 K 0.50000000000000 0.83333333333333 0.50000000000000 39.098 > 82 99 K 0.83333333333333 0.83333333333333 0.50000000000000 39.098 > 82 100 K 0.16666666666667 0.16666666666667 0.83333333333333 39.098 > 82 101 K 0.50000000000000 0.16666666666667 0.83333333333333 39.098 > 82 102 K 0.83333333333333 0.16666666666667 0.83333333333333 39.098 > 82 103 K 0.16666666666667 0.50000000000000 0.83333333333333 39.098 > 82 104 K 0.50000000000000 0.50000000000000 0.83333333333333 39.098 > 82 105 K 0.83333333333333 0.50000000000000 0.83333333333333 39.098 > 82 106 K 0.16666666666667 0.83333333333333 0.83333333333333 39.098 > 82 107 K 0.50000000000000 0.83333333333333 0.83333333333333 39.098 > 82 108 K 0.83333333333333 0.83333333333333 0.83333333333333 39.098 > 82 109 Ta 0.00000000000000 0.00000000000000 0.00000000000000 180.948 > 109 110 Ta 0.33333333333333 0.00000000000000 0.00000000000000 180.948 > 109 111 Ta 0.66666666666667 0.00000000000000 0.00000000000000 180.948 > 109 112 Ta 0.00000000000000 0.33333333333333 0.00000000000000 180.948 > 109 113 Ta 0.33333333333333 0.33333333333333 0.00000000000000 180.948 > 109 114 Ta 0.66666666666667 0.33333333333333 0.00000000000000 180.948 > 109 115 Ta 0.00000000000000 0.66666666666667 0.00000000000000 180.948 > 109 116 Ta 0.33333333333333 0.66666666666667 0.00000000000000 180.948 > 109 117 Ta 0.66666666666667 0.66666666666667 0.00000000000000 180.948 > 109 118 Ta 0.00000000000000 0.00000000000000 0.33333333333333 180.948 > 109 119 Ta 0.33333333333333 0.00000000000000 0.33333333333333 180.948 > 109 120 Ta 0.66666666666667 0.00000000000000 0.33333333333333 180.948 > 109 121 Ta 0.00000000000000 0.33333333333333 0.33333333333333 180.948 > 109 122 Ta 0.33333333333333 0.33333333333333 0.33333333333333 180.948 > 109 123 Ta 0.66666666666667 0.33333333333333 0.33333333333333 180.948 > 109 124 Ta 0.00000000000000 0.66666666666667 0.33333333333333 180.948 > 109 125 Ta 0.33333333333333 0.66666666666667 0.33333333333333 180.948 > 109 126 Ta 0.66666666666667 0.66666666666667 0.33333333333333 180.948 > 109 127 Ta 0.00000000000000 0.00000000000000 0.66666666666667 180.948 > 109 128 Ta 0.33333333333333 0.00000000000000 0.66666666666667 180.948 > 109 129 Ta 0.66666666666667 0.00000000000000 0.66666666666667 180.948 > 109 130 Ta 0.00000000000000 0.33333333333333 0.66666666666667 180.948 > 109 131 Ta 0.33333333333333 0.33333333333333 0.66666666666667 180.948 > 109 132 Ta 0.66666666666667 0.33333333333333 0.66666666666667 180.948 > 109 133 Ta 0.00000000000000 0.66666666666667 0.66666666666667 180.948 > 109 134 Ta 0.33333333333333 0.66666666666667 0.66666666666667 180.948 > 109 135 Ta 0.66666666666667 0.66666666666667 0.66666666666667 180.948 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 5.7507916 0.0000000 0.0000000 0.0000000 5.7507916 0.0000000 0.0000000 0.0000000 5.7507916 -------------------------- Born effective charges -------------------------- 1 O -1.7009814 0.0000000 0.0000000 0.0000000 -6.7657002 0.0000000 0.0000000 0.0000000 -1.7009814 2 O -1.7009814 0.0000000 0.0000000 0.0000000 -1.7009814 0.0000000 0.0000000 0.0000000 -6.7657002 3 O -6.7657002 0.0000000 0.0000000 0.0000000 -1.7009814 0.0000000 0.0000000 0.0000000 -1.7009814 4 K 1.1716213 0.0000000 0.0000000 0.0000000 1.1716213 0.0000000 0.0000000 0.0000000 1.1716213 5 Ta 8.9960418 0.0000000 0.0000000 0.0000000 8.9960418 0.0000000 0.0000000 0.0000000 8.9960418 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz) Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 13 13 13 ] 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: 1.04, Number of G-points: 305, Lambda: 0.26 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/84) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 84 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.103 ( 0.000 0.000 0.000) 0.000 2.103 ( 0.000 0.000 0.000) 0.000 2.103 ( 0.000 0.000 0.000) 0.000 5.808 ( 0.000 0.000 0.000) 0.000 5.808 ( 0.000 0.000 0.000) 0.000 5.808 ( 0.000 0.000 0.000) 0.000 7.910 ( 0.000 0.000 0.000) 0.000 7.910 ( 0.000 0.000 0.000) 0.000 7.910 ( 0.000 0.000 0.000) 0.000 15.778 ( 0.000 0.000 0.000) 0.000 15.778 ( 0.000 0.000 0.000) 0.000 15.778 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/84) ======================= q-point: ( 0.08 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.636 ( 25.786 0.000 0.000) 25.786 0.636 ( 25.786 0.000 0.000) 25.786 1.500 ( 75.738 0.000 0.000) 75.738 2.459 ( 34.474 0.000 0.000) 34.474 2.459 ( 34.474 0.000 0.000) 34.474 5.432 ( 6.362 0.000 0.000) 6.362 5.814 ( 0.538 0.000 0.000) 0.538 5.814 ( 0.538 0.000 0.000) 0.538 7.955 ( 4.512 0.000 0.000) 4.512 7.987 ( 8.183 0.000 0.000) 8.183 7.987 ( 8.183 0.000 0.000) 8.183 12.061 ( 3.010 0.000 0.000) 3.010 15.739 ( -3.984 0.000 0.000) 3.984 15.739 ( -3.984 0.000 0.000) 3.984 23.731 ( 7.478 0.000 0.000) 7.478 ======================= Grid point 2 (3/84) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.015 ( 15.355 0.000 0.000) 15.355 1.015 ( 15.355 0.000 0.000) 15.355 2.877 ( 66.553 0.000 0.000) 66.553 3.254 ( 43.529 0.000 0.000) 43.529 3.254 ( 43.529 0.000 0.000) 43.529 5.603 ( 11.020 0.000 0.000) 11.020 5.824 ( 0.419 0.000 0.000) 0.419 5.824 ( 0.419 0.000 0.000) 0.419 8.076 ( 7.864 0.000 0.000) 7.864 8.234 ( 17.628 0.000 0.000) 17.628 8.234 ( 17.628 0.000 0.000) 17.628 12.143 ( 5.284 0.000 0.000) 5.284 15.629 ( -7.368 0.000 0.000) 7.368 15.629 ( -7.368 0.000 0.000) 7.368 23.937 ( 13.499 0.000 0.000) 13.499 ======================= Grid point 3 (4/84) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.271 ( 11.723 0.000 0.000) 11.723 1.271 ( 11.723 0.000 0.000) 11.723 4.019 ( 34.570 0.000 0.000) 34.570 4.019 ( 34.570 0.000 0.000) 34.570 4.037 ( 53.589 0.000 0.000) 53.589 5.823 ( -0.687 0.000 0.000) 0.687 5.823 ( -0.687 0.000 0.000) 0.687 5.836 ( 12.686 0.000 0.000) 12.686 8.245 ( 9.282 0.000 0.000) 9.282 8.662 ( 26.266 0.000 0.000) 26.266 8.662 ( 26.266 0.000 0.000) 26.266 12.256 ( 6.294 0.000 0.000) 6.294 15.463 ( -9.579 0.000 0.000) 9.579 15.463 ( -9.579 0.000 0.000) 9.579 24.233 ( 16.683 0.000 0.000) 16.683 ======================= Grid point 4 (5/84) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.471 ( 9.019 0.000 0.000) 9.019 1.471 ( 9.019 0.000 0.000) 9.019 4.569 ( 22.860 0.000 0.000) 22.860 4.569 ( 22.860 0.000 0.000) 22.860 4.932 ( 39.358 0.000 0.000) 39.358 5.794 ( -2.331 0.000 0.000) 2.331 5.794 ( -2.331 0.000 0.000) 2.331 6.068 ( 10.714 0.000 0.000) 10.714 8.420 ( 8.551 0.000 0.000) 8.551 9.205 ( 28.787 0.000 0.000) 28.787 9.205 ( 28.787 0.000 0.000) 28.787 12.375 ( 5.850 0.000 0.000) 5.850 15.271 ( -10.009 0.000 0.000) 10.009 15.271 ( -10.009 0.000 0.000) 10.009 24.554 ( 16.072 0.000 0.000) 16.072 ======================= Grid point 5 (6/84) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.614 ( 5.789 0.000 0.000) 5.789 1.614 ( 5.789 0.000 0.000) 5.789 4.917 ( 13.638 0.000 0.000) 13.638 4.917 ( 13.638 0.000 0.000) 13.638 5.555 ( 25.417 0.000 0.000) 25.417 5.737 ( -3.370 0.000 0.000) 3.370 5.737 ( -3.370 0.000 0.000) 3.370 6.228 ( 5.541 0.000 0.000) 5.541 8.562 ( 5.958 0.000 0.000) 5.958 9.711 ( 22.304 0.000 0.000) 22.304 9.711 ( 22.304 0.000 0.000) 22.304 12.473 ( 4.104 0.000 0.000) 4.104 15.094 ( -7.947 0.000 0.000) 7.947 15.094 ( -7.947 0.000 0.000) 7.947 24.826 ( 11.587 0.000 0.000) 11.587 ======================= Grid point 6 (7/84) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.690 ( 2.002 0.000 0.000) 2.002 1.690 ( 2.002 0.000 0.000) 2.002 5.099 ( 5.037 0.000 0.000) 5.037 5.099 ( 5.037 0.000 0.000) 5.037 5.680 ( -1.975 0.000 0.000) 1.975 5.680 ( -1.975 0.000 0.000) 1.975 5.903 ( 10.133 0.000 0.000) 10.133 6.279 ( 0.318 0.000 0.000) 0.318 8.641 ( 2.129 0.000 0.000) 2.129 10.015 ( 8.405 0.000 0.000) 8.405 10.015 ( 8.405 0.000 0.000) 8.405 12.528 ( 1.472 0.000 0.000) 1.472 14.984 ( -3.104 0.000 0.000) 3.104 14.984 ( -3.104 0.000 0.000) 3.104 24.982 ( 4.216 0.000 0.000) 4.216 ======================= Grid point 14 (8/84) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.820 ( 14.117 14.117 0.000) 19.965 1.264 ( 29.512 29.512 0.000) 41.737 1.909 ( 47.624 47.624 0.000) 67.350 2.800 ( 32.539 32.539 0.000) 46.017 3.128 ( 44.929 44.929 0.000) 63.539 5.353 ( -0.875 -0.875 0.000) 1.238 5.827 ( 1.340 1.340 0.000) 1.895 5.907 ( 5.200 5.200 0.000) 7.355 7.886 ( -1.328 -1.328 0.000) 1.878 8.057 ( 7.360 7.360 0.000) 10.409 8.160 ( 12.265 12.265 0.000) 17.345 12.094 ( 3.155 3.155 0.000) 4.462 15.700 ( -4.025 -4.025 0.000) 5.693 15.801 ( 1.404 1.404 0.000) 1.986 23.718 ( 3.559 3.559 0.000) 5.034 ======================= Grid point 15 (9/84) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.079 ( 12.441 5.995 0.000) 13.810 1.631 ( 14.116 44.491 0.000) 46.677 3.037 ( 60.723 17.806 0.000) 63.280 3.543 ( 40.480 27.418 0.000) 48.892 3.937 ( 38.550 56.510 0.000) 68.407 5.441 ( 9.105 -10.089 0.000) 13.590 5.863 ( 2.303 4.165 0.000) 4.759 6.000 ( 5.011 12.020 0.000) 13.022 7.909 ( 3.463 -8.545 0.000) 9.220 8.281 ( 16.277 5.100 0.000) 17.057 8.455 ( 18.569 14.760 0.000) 23.721 12.177 ( 5.356 3.402 0.000) 6.345 15.588 ( -7.450 -4.149 0.000) 8.528 15.755 ( -5.414 10.231 0.000) 11.575 23.897 ( 13.964 -1.237 0.000) 14.019 ======================= Grid point 16 (10/84) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.301 ( 10.631 2.972 0.000) 11.039 1.874 ( 11.524 46.676 0.000) 48.078 4.110 ( 49.034 8.802 0.000) 49.817 4.247 ( 31.119 21.084 0.000) 37.589 4.570 ( 26.723 45.660 0.000) 52.905 5.678 ( 14.891 -3.641 0.000) 15.330 5.910 ( 2.372 9.661 0.000) 9.948 6.107 ( 6.036 18.241 0.000) 19.214 8.012 ( 6.933 -16.067 0.000) 17.499 8.684 ( 25.109 2.396 0.000) 25.223 8.872 ( 24.435 14.851 0.000) 28.594 12.293 ( 6.504 3.700 0.000) 7.483 15.421 ( -9.702 -4.356 0.000) 10.635 15.616 ( -8.654 14.473 0.000) 16.863 24.217 ( 18.407 -0.238 0.000) 18.409 ======================= Grid point 17 (11/84) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.485 ( 8.465 1.499 0.000) 8.597 2.078 ( 9.600 48.280 0.000) 49.225 4.724 ( 18.790 13.362 0.000) 23.057 4.791 ( 17.867 7.033 0.000) 19.201 5.076 ( 30.555 20.408 0.000) 36.744 5.948 ( 1.461 17.153 0.000) 17.215 5.975 ( 14.744 16.406 0.000) 22.057 6.223 ( 5.757 22.345 0.000) 23.075 8.159 ( 7.879 -21.608 0.000) 22.999 9.211 ( 28.254 0.625 0.000) 28.261 9.370 ( 26.319 12.762 0.000) 29.250 12.418 ( 6.230 4.272 0.000) 7.554 15.226 ( -10.160 -4.633 0.000) 11.166 15.438 ( -9.458 16.676 0.000) 19.172 24.573 ( 17.771 2.620 0.000) 17.963 ======================= Grid point 18 (12/84) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.621 ( 5.506 0.694 0.000) 5.549 2.236 ( 6.593 50.080 0.000) 50.512 4.975 ( 4.765 -1.298 0.000) 4.938 4.992 ( 9.610 5.734 0.000) 11.191 5.613 ( 22.572 8.703 0.000) 24.192 5.965 ( 0.414 24.430 0.000) 24.433 6.206 ( 8.709 38.932 0.000) 39.894 6.319 ( 4.008 17.934 0.000) 18.377 8.296 ( 5.966 -24.129 0.000) 24.855 9.711 ( 22.189 0.005 0.000) 22.189 9.836 ( 20.698 10.331 0.000) 23.133 12.523 ( 4.495 5.011 0.000) 6.731 15.046 ( -8.089 -4.936 0.000) 9.476 15.269 ( -7.624 17.953 0.000) 19.505 24.873 ( 12.726 5.190 0.000) 13.744 ======================= Grid point 19 (13/84) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.693 ( 1.914 0.326 0.000) 1.941 2.324 ( 2.367 51.370 0.000) 51.424 5.026 ( 1.101 -7.629 0.000) 7.708 5.111 ( 2.964 0.832 0.000) 3.079 5.914 ( 8.238 4.856 0.000) 9.563 5.968 ( 0.003 29.055 0.000) 29.055 6.310 ( 2.419 55.974 0.000) 56.026 6.371 ( 1.287 9.646 0.000) 9.732 8.377 ( 2.189 -24.847 0.000) 24.943 10.015 ( 8.407 -0.075 0.000) 8.407 10.120 ( 7.887 8.965 0.000) 11.940 12.584 ( 1.639 5.545 0.000) 5.782 14.934 ( -3.166 -5.153 0.000) 6.048 15.163 ( -2.992 18.586 0.000) 18.825 25.043 ( 4.610 6.608 0.000) 8.057 ======================= Grid point 28 (14/84) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.219 ( 8.093 8.093 0.000) 11.446 2.329 ( 25.944 25.944 0.000) 36.690 3.608 ( 39.987 39.987 0.000) 56.551 4.158 ( 32.416 32.416 0.000) 45.844 4.990 ( 45.925 45.925 0.000) 64.948 5.305 ( -1.366 -1.366 0.000) 1.932 5.995 ( 9.832 9.832 0.000) 13.905 6.239 ( 13.390 13.390 0.000) 18.936 7.807 ( -2.268 -2.268 0.000) 3.208 8.443 ( 12.169 12.169 0.000) 17.210 8.764 ( 17.727 17.727 0.000) 25.070 12.266 ( 5.727 5.727 0.000) 8.099 15.473 ( -7.700 -7.700 0.000) 10.889 15.915 ( 4.998 4.998 0.000) 7.068 23.982 ( 10.792 10.792 0.000) 15.262 ======================= Grid point 29 (15/84) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.378 ( 8.206 4.882 0.000) 9.548 2.705 ( 15.516 36.112 0.000) 39.305 4.402 ( 37.572 21.751 0.000) 43.414 4.679 ( 20.499 20.046 0.000) 28.672 5.222 ( 1.964 13.707 0.000) 13.847 5.806 ( 29.812 18.424 0.000) 35.046 6.234 ( 14.105 24.565 0.000) 28.327 6.555 ( 20.378 26.455 0.000) 33.394 7.799 ( 0.826 -2.418 0.000) 2.555 8.763 ( 21.127 6.290 0.000) 22.044 9.149 ( 22.272 13.995 0.000) 26.304 12.391 ( 7.074 6.294 0.000) 9.469 15.299 ( -10.073 -8.117 0.000) 12.936 15.941 ( -1.823 17.259 0.000) 17.355 24.289 ( 19.792 8.944 0.000) 21.719 ======================= Grid point 30 (16/84) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.526 ( 6.982 2.652 0.000) 7.468 2.971 ( 12.515 40.743 0.000) 42.621 4.859 ( 8.073 3.749 0.000) 8.901 4.957 ( 9.402 9.020 0.000) 13.029 5.427 ( 20.193 15.150 0.000) 25.244 6.221 ( 15.144 10.640 0.000) 18.508 6.489 ( 11.347 38.174 0.000) 39.825 6.994 ( 23.046 36.316 0.000) 43.011 7.822 ( 1.156 -1.015 0.000) 1.538 9.232 ( 26.305 1.777 0.000) 26.365 9.604 ( 24.123 10.905 0.000) 26.474 12.530 ( 7.042 7.023 0.000) 9.945 15.097 ( -10.619 -8.682 0.000) 13.717 15.868 ( -5.233 26.315 0.000) 26.830 24.691 ( 20.740 10.529 0.000) 23.260 ======================= Grid point 31 (17/84) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.640 ( 4.676 1.216 0.000) 4.832 3.184 ( 9.212 44.540 0.000) 45.483 4.888 ( -1.532 -4.876 0.000) 5.111 5.078 ( 3.892 3.004 0.000) 4.917 5.814 ( 16.809 12.391 0.000) 20.882 6.442 ( 8.376 8.524 0.000) 11.950 6.655 ( 5.995 45.077 0.000) 45.474 7.380 ( 16.374 41.923 0.000) 45.007 7.833 ( -0.154 -4.066 0.000) 4.069 9.712 ( 21.761 0.097 0.000) 21.762 10.034 ( 19.289 9.482 0.000) 21.494 12.651 ( 5.273 7.797 0.000) 9.412 14.908 ( -8.524 -9.308 0.000) 12.622 15.762 ( -5.170 32.138 0.000) 32.551 25.045 ( 15.112 13.162 0.000) 20.040 ======================= Grid point 32 (18/84) ======================= q-point: ( 0.46 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.702 ( 1.645 0.507 0.000) 1.721 3.310 ( 3.529 47.471 0.000) 47.602 4.856 ( -1.134 -8.759 0.000) 8.832 5.124 ( 1.070 0.548 0.000) 1.202 6.033 ( 5.733 8.897 0.000) 10.584 6.550 ( 2.815 8.566 0.000) 9.016 6.728 ( 1.777 47.578 0.000) 47.612 7.611 ( 7.165 60.362 0.000) 60.786 7.813 ( -1.656 -22.086 0.000) 22.148 10.013 ( 8.407 -0.131 0.000) 8.408 10.300 ( 7.427 9.029 0.000) 11.691 12.723 ( 1.966 8.324 0.000) 8.553 14.789 ( -3.358 -9.764 0.000) 10.325 15.688 ( -2.156 35.089 0.000) 35.155 25.247 ( 5.481 14.880 0.000) 15.857 ======================= Grid point 42 (19/84) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.480 ( 5.540 5.540 0.000) 7.834 3.241 ( 21.138 21.138 0.000) 29.894 4.836 ( 16.562 16.562 0.000) 23.422 4.956 ( 9.204 9.204 0.000) 13.016 5.387 ( 12.954 12.954 0.000) 18.320 6.159 ( 15.584 15.584 0.000) 22.040 6.799 ( 30.988 30.988 0.000) 43.824 6.911 ( 8.613 8.613 0.000) 12.181 8.023 ( 23.319 23.319 0.000) 32.978 8.947 ( 13.437 13.437 0.000) 19.002 9.433 ( 16.303 16.303 0.000) 23.056 12.528 ( 7.668 7.668 0.000) 10.844 15.116 ( -10.698 -10.698 0.000) 15.130 16.225 ( 11.510 11.510 0.000) 16.277 24.570 ( 19.522 19.522 0.000) 27.609 ======================= Grid point 43 (20/84) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.584 ( 5.002 3.200 0.000) 5.937 3.586 ( 15.816 23.475 0.000) 28.306 4.933 ( -2.611 2.552 0.000) 3.650 5.072 ( 3.626 3.585 0.000) 5.099 5.781 ( 22.409 23.017 0.000) 32.124 6.402 ( 9.765 7.239 0.000) 12.156 7.015 ( 5.193 -14.532 0.000) 15.432 7.315 ( 21.026 44.509 0.000) 49.225 8.475 ( 19.993 47.827 0.000) 51.838 9.290 ( 21.553 4.628 0.000) 22.045 9.795 ( 21.006 9.441 0.000) 23.031 12.677 ( 7.573 8.008 0.000) 11.021 14.900 ( -11.408 -11.562 0.000) 16.243 16.381 ( 5.067 25.443 0.000) 25.943 24.990 ( 22.590 20.124 0.000) 30.254 ======================= Grid point 44 (21/84) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.666 ( 3.457 1.495 0.000) 3.766 3.866 ( 13.183 25.499 0.000) 28.705 4.835 ( -6.456 -0.529 0.000) 6.477 5.117 ( 1.398 1.270 0.000) 1.889 6.129 ( 13.121 20.793 0.000) 24.587 6.534 ( 4.144 -2.894 0.000) 5.054 7.143 ( 7.315 -18.925 0.000) 20.290 7.602 ( 9.507 50.973 0.000) 51.852 8.757 ( 9.717 54.937 0.000) 55.789 9.716 ( 20.655 0.462 0.000) 20.660 10.196 ( 18.908 7.460 0.000) 20.326 12.808 ( 5.680 8.112 0.000) 9.904 14.695 ( -9.288 -12.544 0.000) 15.608 16.437 ( 1.272 35.920 0.000) 35.943 25.384 ( 17.151 21.729 0.000) 27.682 ======================= Grid point 45 (22/84) ======================= q-point: ( 0.46 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.712 ( 1.233 0.602 0.000) 1.372 4.067 ( 6.466 29.700 0.000) 30.396 4.721 ( -4.193 -4.766 0.000) 6.348 5.133 ( 0.366 0.414 0.000) 0.553 6.291 ( 4.084 18.310 0.000) 18.760 6.579 ( 0.945 -11.157 0.000) 11.197 7.254 ( 3.345 -16.176 0.000) 16.518 7.712 ( 2.578 52.738 0.000) 52.801 8.873 ( 2.803 55.738 0.000) 55.809 10.010 ( 8.390 -0.152 0.000) 8.391 10.461 ( 7.489 7.741 0.000) 10.771 12.885 ( 2.131 8.130 0.000) 8.404 14.565 ( -3.701 -13.278 0.000) 13.784 16.447 ( 0.070 41.469 0.000) 41.469 25.616 ( 6.302 22.974 0.000) 23.823 ======================= Grid point 56 (23/84) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.645 ( 3.020 3.020 0.000) 4.271 3.938 ( 14.955 14.955 0.000) 21.150 4.930 ( -2.727 -2.727 0.000) 3.857 5.116 ( 1.406 1.406 0.000) 1.989 6.252 ( 22.944 22.944 0.000) 32.448 6.417 ( -7.021 -7.021 0.000) 9.929 6.869 ( 2.561 2.561 0.000) 3.622 8.089 ( 33.247 33.247 0.000) 47.019 9.270 ( 32.015 32.015 0.000) 45.276 9.438 ( 11.622 11.622 0.000) 16.436 10.001 ( 13.341 13.341 0.000) 18.867 12.828 ( 7.412 7.412 0.000) 10.482 14.664 ( -12.529 -12.529 0.000) 17.718 16.805 ( 18.078 18.078 0.000) 25.566 25.426 ( 23.758 23.758 0.000) 33.599 ======================= Grid point 57 (24/84) ======================= q-point: ( 0.38 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.695 ( 2.148 1.459 0.000) 2.597 4.204 ( 13.114 11.873 0.000) 17.690 4.838 ( -6.622 -0.159 0.000) 6.624 5.134 ( 0.539 0.592 0.000) 0.801 6.238 ( -8.799 -22.421 0.000) 24.086 6.583 ( 10.992 23.893 0.000) 26.301 7.029 ( 10.765 3.090 0.000) 11.200 8.539 ( 14.102 44.612 0.000) 46.787 9.669 ( 11.076 39.793 0.000) 41.306 9.735 ( 17.520 1.831 0.000) 17.615 10.327 ( 18.456 6.385 0.000) 19.530 12.953 ( 5.330 6.773 0.000) 8.619 14.437 ( -10.418 -13.856 0.000) 17.336 17.073 ( 9.933 28.801 0.000) 30.466 25.846 ( 18.502 24.954 0.000) 31.064 ======================= Grid point 58 (25/84) ======================= q-point: ( 0.46 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.724 ( 0.777 0.606 0.000) 0.985 4.442 ( 10.952 10.452 0.000) 15.139 4.685 ( -8.709 0.923 0.000) 8.757 5.140 ( 0.139 0.326 0.000) 0.355 6.119 ( -3.257 -29.741 0.000) 29.919 6.708 ( 2.879 22.966 0.000) 23.146 7.192 ( 4.753 3.648 0.000) 5.992 8.692 ( 3.258 47.635 0.000) 47.747 9.786 ( 2.515 39.871 0.000) 39.950 10.007 ( 8.303 -0.126 0.000) 8.304 10.599 ( 7.915 6.608 0.000) 10.311 13.025 ( 1.967 6.227 0.000) 6.530 14.290 ( -4.226 -14.896 0.000) 15.484 17.196 ( 3.082 34.592 0.000) 34.729 26.097 ( 6.878 25.806 0.000) 26.707 ======================= Grid point 70 (26/84) ======================= q-point: ( 0.38 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.720 ( 1.069 1.069 0.000) 1.512 4.395 ( 8.991 8.991 0.000) 12.715 4.809 ( -3.076 -3.076 0.000) 4.350 5.142 ( 0.276 0.276 0.000) 0.391 5.830 ( -18.099 -18.099 0.000) 25.595 6.967 ( 13.777 13.777 0.000) 19.483 7.132 ( 7.141 7.141 0.000) 10.100 9.237 ( 25.310 25.310 0.000) 35.794 9.815 ( 7.722 7.722 0.000) 10.921 10.265 ( 19.371 19.371 0.000) 27.395 10.468 ( 10.761 10.761 0.000) 15.218 13.062 ( 4.467 4.467 0.000) 6.317 14.183 ( -11.835 -11.835 0.000) 16.738 17.518 ( 17.250 17.250 0.000) 24.396 26.288 ( 19.490 19.490 0.000) 27.563 ======================= Grid point 71 (27/84) ======================= q-point: ( 0.46 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.735 ( 0.393 0.464 0.000) 0.608 4.564 ( 8.459 4.623 0.000) 9.640 4.718 ( -6.302 1.339 0.000) 6.442 5.145 ( 0.082 0.216 0.000) 0.231 5.583 ( -6.809 -24.203 0.000) 25.142 7.101 ( 2.395 16.474 0.000) 16.647 7.272 ( 4.849 3.790 0.000) 6.155 9.493 ( 4.351 34.161 0.000) 34.437 10.006 ( 7.885 0.010 0.000) 7.885 10.422 ( 2.366 26.070 0.000) 26.177 10.711 ( 8.402 4.776 0.000) 9.665 13.121 ( 1.570 3.741 0.000) 4.057 14.014 ( -4.920 -13.021 0.000) 13.919 17.737 ( 5.617 21.196 0.000) 21.928 26.553 ( 7.278 20.087 0.000) 21.365 ======================= Grid point 84 (28/84) ======================= q-point: ( 0.46 0.46 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.741 ( 0.175 0.175 0.000) 0.248 4.656 ( 4.654 4.654 0.000) 6.582 4.720 ( -1.294 -1.294 0.000) 1.830 5.148 ( 0.076 0.076 0.000) 0.107 5.236 ( -10.438 -10.438 0.000) 14.762 7.312 ( 4.324 4.324 0.000) 6.115 7.327 ( 2.668 2.668 0.000) 3.773 9.945 ( 10.139 10.139 0.000) 14.338 10.018 ( 2.691 2.691 0.000) 3.806 10.764 ( 6.549 6.549 0.000) 9.261 10.781 ( 4.639 4.639 0.000) 6.560 13.169 ( 1.229 1.229 0.000) 1.739 13.825 ( -5.537 -5.537 0.000) 7.830 18.008 ( 6.973 6.973 0.000) 9.861 26.826 ( 7.506 7.506 0.000) 10.615 ======================= Grid point 183 (29/84) ======================= q-point: ( 0.08 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.368 ( 20.846 20.846 20.846) 36.106 1.368 ( 20.846 20.846 20.846) 36.106 2.242 ( 36.885 36.885 36.885) 63.886 3.371 ( 34.758 34.758 34.758) 60.202 3.371 ( 34.758 34.758 34.758) 60.202 5.281 ( -2.718 -2.718 -2.718) 4.708 5.912 ( 3.652 3.652 3.652) 6.326 5.912 ( 3.652 3.652 3.652) 6.326 7.820 ( -3.183 -3.183 -3.183) 5.513 8.242 ( 11.254 11.254 11.254) 19.493 8.242 ( 11.254 11.254 11.254) 19.493 12.122 ( 2.881 2.881 2.881) 4.990 15.765 ( -0.229 -0.229 -0.229) 0.396 15.765 ( -0.229 -0.229 -0.229) 0.396 23.713 ( 2.717 2.717 2.717) 4.706 ======================= Grid point 184 (30/84) ======================= q-point: ( 0.15 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.580 ( 11.836 23.673 23.673) 35.509 1.779 ( 14.084 25.470 25.470) 38.676 3.194 ( 54.600 17.273 17.273) 59.815 3.928 ( 31.241 30.479 30.479) 53.234 4.249 ( 41.212 33.052 33.052) 62.316 5.341 ( 9.780 -7.522 -7.522) 14.451 5.929 ( 1.217 5.536 5.536) 7.924 6.066 ( 9.579 6.381 6.381) 13.159 7.803 ( 1.496 -7.553 -7.553) 10.785 8.407 ( 10.548 12.188 12.188) 20.208 8.612 ( 22.733 12.924 12.924) 29.169 12.202 ( 5.291 2.603 2.603) 6.446 15.661 ( -6.915 1.868 1.868) 7.402 15.759 ( -3.788 3.839 3.839) 6.620 23.879 ( 14.024 0.211 0.211) 14.028 ======================= Grid point 185 (31/84) ======================= q-point: ( 0.23 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.808 ( 11.567 23.206 23.206) 34.797 2.000 ( 9.892 27.365 27.365) 39.944 4.168 ( 43.798 7.124 7.124) 44.942 4.482 ( 24.847 20.686 20.686) 38.382 4.846 ( 21.172 23.659 23.659) 39.595 5.677 ( 23.300 1.499 1.499) 23.396 5.958 ( 1.736 7.172 7.172) 10.291 6.292 ( 14.083 16.742 16.742) 27.548 7.872 ( 5.427 -12.155 -12.155) 18.026 8.667 ( 16.817 4.952 4.952) 18.217 9.104 ( 27.853 15.909 15.909) 35.805 12.319 ( 6.640 2.659 2.659) 7.631 15.506 ( -9.003 2.375 2.375) 9.609 15.638 ( -8.227 7.560 7.560) 13.490 24.213 ( 19.680 0.756 0.756) 19.709 ======================= Grid point 186 (32/84) ======================= q-point: ( 0.31 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.016 ( 9.835 23.838 23.838) 35.118 2.173 ( 8.126 27.196 27.196) 39.310 4.758 ( 15.836 -1.761 -1.761) 16.030 4.860 ( 14.552 12.206 12.206) 22.578 5.247 ( 24.355 13.326 13.326) 30.794 5.993 ( 1.807 10.507 10.507) 14.968 6.089 ( 17.179 6.353 6.353) 19.387 6.592 ( 15.857 36.732 36.732) 54.313 7.996 ( 6.892 -15.774 -15.774) 23.348 9.051 ( 22.362 -4.499 -4.499) 23.249 9.652 ( 28.005 17.119 17.119) 37.019 12.448 ( 6.537 3.074 3.074) 7.850 15.325 ( -9.425 2.954 2.954) 10.310 15.462 ( -9.625 9.243 9.243) 16.233 24.597 ( 19.187 3.023 3.023) 19.657 ======================= Grid point 187 (33/84) ======================= q-point: ( 0.38 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.177 ( 6.665 24.759 24.759) 35.643 2.308 ( 5.689 26.997 26.997) 38.601 4.909 ( 3.352 -5.390 -5.390) 8.327 5.059 ( 6.695 5.460 5.460) 10.220 5.697 ( 19.684 7.117 7.117) 22.108 6.024 ( 1.311 14.283 14.283) 20.242 6.332 ( 8.635 5.410 5.410) 11.537 6.850 ( 10.154 51.329 51.329) 73.297 8.119 ( 5.447 -17.865 -17.865) 25.845 9.474 ( 19.879 -10.290 -10.290) 24.636 10.136 ( 21.081 17.003 17.003) 31.978 12.560 ( 4.820 3.679 3.679) 7.092 15.158 ( -7.499 3.482 3.482) 8.971 15.287 ( -7.952 9.882 9.882) 16.079 24.920 ( 13.714 5.206 5.206) 15.566 ======================= Grid point 188 (34/84) ======================= q-point: ( 0.46 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.265 ( 2.367 25.402 25.402) 36.001 2.384 ( 2.071 26.989 26.989) 38.225 4.943 ( 0.656 -6.939 -6.939) 9.835 5.138 ( 1.851 1.676 1.676) 3.007 5.960 ( 7.085 3.604 3.604) 8.728 6.041 ( 0.470 16.637 16.637) 23.533 6.436 ( 2.491 6.378 6.378) 9.357 6.977 ( 3.215 56.432 56.432) 79.872 8.193 ( 2.033 -18.726 -18.726) 26.560 9.754 ( 7.927 -12.139 -12.139) 18.908 10.422 ( 7.887 16.609 16.609) 24.778 12.625 ( 1.780 4.129 4.129) 6.105 15.054 ( -2.933 3.807 3.807) 6.131 15.176 ( -3.151 10.058 10.058) 14.569 25.103 ( 4.955 6.441 6.441) 10.369 ======================= Grid point 197 (35/84) ======================= q-point: ( 0.15 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.802 ( 9.070 9.070 44.828) 46.627 2.390 ( 25.659 25.659 6.194) 36.813 3.717 ( 35.316 35.316 13.018) 51.613 4.489 ( 26.697 26.697 30.633) 48.620 5.057 ( 41.550 41.550 6.424) 59.111 5.238 ( 0.494 0.494 -5.181) 5.228 6.062 ( 8.191 8.191 6.890) 13.478 6.244 ( 13.370 13.370 0.718) 18.922 7.690 ( -3.635 -3.635 -8.974) 10.342 8.682 ( 13.755 13.755 21.443) 28.952 8.913 ( 19.844 19.844 13.373) 31.086 12.277 ( 5.178 5.178 1.323) 7.441 15.588 ( -6.756 -6.756 10.487) 14.187 15.884 ( 5.151 5.151 -3.139) 7.932 23.984 ( 11.494 11.494 1.465) 16.321 ======================= Grid point 198 (36/84) ======================= q-point: ( 0.23 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.976 ( 8.909 6.192 48.018) 49.229 2.762 ( 15.114 35.878 5.993) 39.390 4.389 ( 29.543 15.412 2.528) 33.417 4.917 ( 17.126 21.399 19.244) 33.489 5.223 ( 5.552 12.085 1.516) 13.386 5.822 ( 28.986 14.489 1.526) 32.441 6.232 ( 9.205 22.033 2.026) 23.965 6.662 ( 30.739 18.643 8.333) 36.903 7.669 ( 1.461 -5.955 -11.873) 13.363 8.903 ( 10.235 18.625 16.674) 27.013 9.387 ( 28.025 13.332 19.734) 36.777 12.396 ( 6.926 5.301 0.638) 8.745 15.431 ( -9.259 -7.154 12.467) 17.097 15.918 ( -1.472 16.660 -2.158) 16.864 24.304 ( 20.581 9.765 2.350) 22.901 ======================= Grid point 199 (37/84) ======================= q-point: ( 0.31 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.138 ( 7.660 3.795 50.597) 51.314 3.018 ( 11.886 41.187 4.837) 43.139 4.741 ( 8.188 1.022 -5.915) 10.152 5.109 ( 4.258 10.651 9.801) 15.088 5.503 ( 20.712 13.994 7.772) 26.177 6.226 ( 15.205 8.150 0.581) 17.261 6.407 ( 8.547 31.959 -1.913) 33.138 7.282 ( 27.732 12.527 19.707) 36.254 7.761 ( 9.090 9.153 -7.469) 14.906 9.115 ( 12.460 11.615 1.780) 17.127 9.938 ( 27.912 12.203 23.363) 38.390 12.534 ( 7.158 5.753 0.590) 9.202 15.244 ( -9.869 -7.787 14.213) 18.975 15.847 ( -5.256 25.644 -1.624) 26.227 24.723 ( 21.617 10.885 3.725) 24.487 ======================= Grid point 200 (38/84) ======================= q-point: ( 0.38 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.264 ( 5.244 2.042 52.560) 52.860 3.219 ( 8.708 45.402 3.603) 46.370 4.806 ( 0.493 -4.235 -5.993) 7.355 5.153 ( 0.990 3.591 5.732) 6.836 5.864 ( 15.019 10.794 5.187) 19.209 6.455 ( 8.944 6.737 1.335) 11.277 6.545 ( 5.524 37.532 -3.592) 38.106 7.608 ( 8.372 -21.657 0.598) 23.227 8.015 ( 13.829 52.100 29.082) 61.249 9.383 ( 14.411 0.992 -15.939) 21.511 10.416 ( 20.642 11.854 24.022) 33.818 12.659 ( 5.470 6.360 0.888) 8.435 15.067 ( -7.968 -8.484 15.720) 19.560 15.739 ( -5.363 31.545 -1.806) 32.049 25.091 ( 15.756 13.052 5.095) 21.085 ======================= Grid point 201 (39/84) ======================= q-point: ( 0.46 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.333 ( 1.875 1.101 53.678) 53.722 3.338 ( 3.327 48.465 2.852) 48.662 4.801 ( -0.422 -6.956 -4.612) 8.357 5.162 ( 0.122 0.847 3.098) 3.214 6.057 ( 5.009 7.661 3.164) 9.685 6.573 ( 3.244 6.915 2.181) 7.943 6.614 ( 1.694 39.378 -3.992) 39.616 7.701 ( 2.203 -29.339 -9.061) 30.785 8.195 ( 4.523 63.755 46.675) 79.144 9.604 ( 6.746 -3.857 -24.219) 25.435 10.695 ( 7.659 11.602 23.712) 27.487 12.733 ( 2.061 6.804 1.175) 7.206 14.956 ( -3.147 -8.977 16.634) 19.162 15.662 ( -2.258 34.553 -2.124) 34.692 25.302 ( 5.712 14.516 5.901) 16.678 ======================= Grid point 211 (40/84) ======================= q-point: ( 0.23 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.098 ( 6.313 6.313 51.262) 52.034 3.296 ( 21.026 21.026 5.711) 30.278 4.675 ( 11.367 11.367 -6.536) 17.353 5.214 ( 6.031 6.031 15.445) 17.644 5.440 ( 17.628 17.628 7.075) 25.915 6.131 ( 15.011 15.011 -2.910) 21.428 6.762 ( 27.544 27.544 -2.741) 39.050 6.916 ( 8.849 8.849 0.513) 12.525 7.746 ( 15.696 15.696 -21.657) 31.012 9.297 ( 18.085 18.085 29.041) 38.698 9.663 ( 18.252 18.252 21.450) 33.561 12.518 ( 7.189 7.189 -0.832) 10.201 15.263 ( -9.941 -9.941 14.547) 20.230 16.196 ( 11.404 11.404 -2.914) 16.389 24.597 ( 20.032 20.032 3.364) 28.529 ======================= Grid point 212 (41/84) ======================= q-point: ( 0.31 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.213 ( 5.548 4.005 53.467) 53.903 3.637 ( 15.549 23.511 5.322) 28.685 4.794 ( 1.731 3.472 -8.062) 8.947 5.231 ( -1.464 2.487 10.690) 11.073 5.842 ( 20.008 21.594 6.353) 30.116 6.375 ( 10.255 6.372 -2.698) 12.371 6.980 ( 4.907 -4.323 -3.883) 7.605 7.225 ( 14.624 29.291 -6.860) 33.450 8.254 ( 33.760 24.121 -21.110) 46.553 9.512 ( 5.969 28.267 24.985) 38.196 10.125 ( 26.778 7.714 26.718) 38.606 12.662 ( 7.495 7.312 -1.392) 10.563 15.061 ( -10.716 -10.860 16.131) 22.204 16.350 ( 4.860 24.885 -3.145) 25.550 25.026 ( 22.987 20.233 4.061) 30.891 ======================= Grid point 213 (42/84) ======================= q-point: ( 0.38 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.305 ( 3.835 2.244 54.961) 55.140 3.912 ( 12.892 25.707 4.772) 29.152 4.773 ( -3.284 0.782 -4.804) 5.872 5.196 ( -1.744 1.099 6.770) 7.077 6.146 ( 11.323 18.865 2.109) 22.103 6.518 ( 4.693 -3.628 -1.500) 6.119 7.095 ( 6.672 -15.464 -5.247) 17.640 7.412 ( 5.988 42.747 -14.321) 45.478 8.888 ( 28.979 28.420 -4.865) 40.880 9.583 ( 2.053 23.641 11.756) 26.482 10.596 ( 20.599 7.179 27.812) 35.346 12.793 ( 5.771 7.295 -1.380) 9.404 14.869 ( -8.758 -11.863 17.454) 22.849 16.401 ( 0.996 35.064 -3.602) 35.263 25.427 ( 17.431 21.430 4.631) 28.010 ======================= Grid point 214 (43/84) ======================= q-point: ( 0.46 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.356 ( 1.376 1.280 55.754) 55.786 4.108 ( 6.290 30.029 4.300) 30.980 4.699 ( -3.119 -3.212 -1.790) 4.822 5.173 ( -0.635 0.350 3.701) 3.771 6.286 ( 3.503 16.441 0.044) 16.810 6.571 ( 1.140 -12.143 -0.819) 12.224 7.199 ( 3.159 -15.342 -6.169) 16.835 7.482 ( 1.665 47.432 -17.096) 50.446 9.289 ( 11.401 44.789 27.358) 53.708 9.602 ( 0.106 6.673 -18.724) 19.878 10.875 ( 7.647 7.410 27.183) 29.194 12.872 ( 2.190 7.253 -1.229) 7.675 14.746 ( -3.492 -12.586 18.269) 22.458 16.407 ( -0.057 40.465 -4.010) 40.663 25.662 ( 6.404 22.436 4.968) 23.855 ======================= Grid point 225 (44/84) ======================= q-point: ( 0.31 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.288 ( 3.607 3.607 55.490) 55.724 3.987 ( 14.746 14.746 5.085) 21.465 4.838 ( 0.583 0.583 -6.172) 6.227 5.236 ( -1.316 -1.316 9.329) 9.512 6.271 ( 20.398 20.398 2.129) 28.926 6.395 ( -6.107 -6.107 -2.345) 8.949 6.834 ( 0.824 0.824 -3.619) 3.802 7.920 ( 30.042 30.042 -14.840) 45.004 8.664 ( 22.552 22.552 -34.272) 46.816 10.045 ( 19.834 19.834 34.604) 44.545 10.293 ( 14.532 14.532 26.545) 33.570 12.805 ( 7.210 7.210 -2.285) 10.449 14.839 ( -11.858 -11.858 17.553) 24.277 16.762 ( 17.429 17.429 -4.386) 25.036 25.463 ( 23.734 23.734 4.022) 33.805 ======================= Grid point 226 (45/84) ======================= q-point: ( 0.38 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.348 ( 2.529 2.123 56.844) 56.940 4.249 ( 12.858 11.652 4.631) 17.960 4.807 ( -3.845 1.693 -2.295) 4.787 5.202 ( -1.804 -0.365 6.276) 6.540 6.225 ( -8.706 -22.295 -1.393) 23.975 6.558 ( 9.289 21.899 -2.340) 23.903 6.974 ( 10.125 1.954 -5.860) 11.860 8.284 ( 9.969 43.131 -21.234) 49.097 9.168 ( 26.737 7.536 -31.158) 41.743 10.194 ( 1.733 33.174 31.256) 45.612 10.704 ( 21.011 4.276 28.660) 35.794 12.928 ( 5.351 6.499 -2.454) 8.769 14.624 ( -9.856 -13.147 18.796) 24.965 17.018 ( 9.376 27.771 -5.599) 29.841 25.882 ( 18.467 24.581 3.888) 30.990 ======================= Grid point 227 (46/84) ======================= q-point: ( 0.46 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.382 ( 0.913 1.316 57.543) 57.565 4.482 ( 10.826 9.988 4.125) 15.296 4.692 ( -7.711 2.438 1.008) 8.150 5.176 ( -0.699 0.066 3.584) 3.652 6.106 ( -3.265 -29.668 -1.347) 29.878 6.663 ( 2.434 21.006 -4.290) 21.577 7.127 ( 4.469 2.650 -6.958) 8.684 8.391 ( 2.317 45.087 -24.018) 51.138 9.563 ( 11.830 2.600 -30.165) 32.506 10.196 ( -0.673 36.287 31.894) 48.315 10.991 ( 7.920 4.876 27.638) 29.161 13.001 ( 2.000 5.917 -2.392) 6.688 14.485 ( -3.985 -14.113 19.623) 24.497 17.133 ( 2.861 33.365 -6.385) 34.090 26.133 ( 6.867 25.236 3.812) 26.430 ======================= Grid point 239 (47/84) ======================= q-point: ( 0.38 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.384 ( 1.542 1.542 58.221) 58.261 4.435 ( 8.770 8.770 4.176) 13.087 4.811 ( -1.548 -1.548 0.547) 2.256 5.189 ( -0.782 -0.782 4.661) 4.790 5.818 ( -18.085 -18.085 -1.212) 25.605 6.908 ( 12.368 12.368 -6.033) 18.502 7.058 ( 6.468 6.468 -7.711) 11.964 8.958 ( 22.827 22.827 -23.227) 39.770 9.328 ( 12.332 12.332 -33.032) 37.354 10.723 ( 14.359 14.359 30.271) 36.451 10.789 ( 11.126 11.126 27.163) 31.391 13.035 ( 4.449 4.449 -2.720) 6.855 14.383 ( -11.156 -11.156 20.098) 25.551 17.445 ( 16.444 16.444 -7.468) 24.426 26.317 ( 19.212 19.212 3.187) 27.356 ======================= Grid point 240 (48/84) ======================= q-point: ( 0.46 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.405 ( 0.567 1.014 58.966) 58.978 4.598 ( 7.874 4.606 3.502) 9.772 4.745 ( -5.173 1.862 2.943) 6.236 5.176 ( -0.439 -0.051 3.122) 3.153 5.572 ( -6.796 -24.191 -1.195) 25.155 7.025 ( 2.011 15.244 -7.877) 17.277 7.190 ( 4.560 3.032 -8.710) 10.288 9.157 ( 2.891 32.896 -26.832) 42.550 9.604 ( 10.105 1.779 -29.930) 31.639 10.783 ( 0.709 24.235 28.151) 37.152 11.067 ( 8.373 3.045 26.558) 28.013 13.094 ( 1.591 3.708 -2.701) 4.855 14.225 ( -4.595 -12.184 21.077) 24.775 17.653 ( 5.315 20.269 -8.585) 22.645 26.578 ( 7.180 19.675 2.791) 21.129 ======================= Grid point 253 (49/84) ======================= q-point: ( 0.46 0.46 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.419 ( 0.384 0.384 59.779) 59.781 4.691 ( 4.684 4.684 3.619) 7.549 4.754 ( -0.991 -0.991 3.533) 3.801 5.175 ( -0.091 -0.091 2.750) 2.753 5.223 ( -10.583 -10.583 -1.364) 15.029 7.220 ( 3.957 3.957 -9.612) 11.123 7.235 ( 2.396 2.396 -9.719) 10.293 9.591 ( 8.945 8.945 -27.927) 30.659 9.640 ( 4.017 4.017 -29.189) 29.737 11.105 ( 5.157 5.157 25.678) 26.694 11.109 ( 4.714 4.714 25.410) 26.270 13.142 ( 1.247 1.247 -2.694) 3.220 14.050 ( -5.099 -5.099 22.284) 23.422 17.911 ( 6.635 6.635 -9.860) 13.611 26.846 ( 7.360 7.360 2.218) 10.642 ======================= Grid point 366 (50/84) ======================= q-point: ( 0.15 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.557 ( 20.493 20.493 20.493) 35.495 2.557 ( 20.493 20.493 20.493) 35.495 4.102 ( 24.433 24.433 24.433) 42.319 5.185 ( 2.715 2.715 2.715) 4.702 5.201 ( 23.046 23.046 23.046) 39.918 5.201 ( 23.046 23.046 23.046) 39.918 6.278 ( 10.121 10.121 10.121) 17.529 6.278 ( 10.121 10.121 10.121) 17.529 7.532 ( -6.864 -6.864 -6.864) 11.889 9.182 ( 20.999 20.999 20.999) 36.371 9.182 ( 20.999 20.999 20.999) 36.371 12.326 ( 3.892 3.892 3.892) 6.741 15.797 ( 1.786 1.786 1.786) 3.093 15.797 ( 1.786 1.786 1.786) 3.093 24.099 ( 11.472 11.472 11.472) 19.870 ======================= Grid point 367 (51/84) ======================= q-point: ( 0.23 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.789 ( 12.955 25.854 25.854) 38.790 2.977 ( 15.331 21.917 21.917) 34.580 4.506 ( 14.572 7.656 7.656) 18.154 5.304 ( 2.922 19.344 19.344) 27.512 5.349 ( 15.643 12.102 12.102) 23.186 5.869 ( 26.478 4.294 4.294) 27.165 6.375 ( 6.676 13.874 13.874) 20.725 6.807 ( 39.147 6.358 6.358) 40.167 7.432 ( -3.590 -10.691 -10.691) 15.540 9.373 ( 7.061 29.587 29.587) 42.434 9.773 ( 31.381 18.805 18.805) 41.134 12.426 ( 6.387 2.828 2.828) 7.536 15.666 ( -7.577 6.237 6.237) 11.628 15.910 ( 2.135 5.283 5.283) 7.771 24.422 ( 21.168 10.853 10.853) 26.147 ======================= Grid point 368 (52/84) ======================= q-point: ( 0.31 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.039 ( 12.611 27.321 27.321) 40.644 3.205 ( 9.309 24.731 24.731) 36.192 4.669 ( 4.042 -1.507 -1.507) 4.570 5.314 ( -1.532 10.917 10.917) 15.515 5.707 ( 18.830 13.338 13.338) 26.653 6.250 ( 15.215 2.123 2.123) 15.509 6.516 ( 7.348 13.895 13.895) 20.979 7.391 ( -0.880 -14.943 -14.943) 21.151 7.673 ( 46.624 6.357 6.357) 47.483 9.460 ( 2.149 30.743 30.743) 43.530 10.356 ( 28.263 18.972 18.972) 38.970 12.561 ( 7.272 2.489 2.489) 8.079 15.515 ( -7.924 7.002 7.002) 12.682 15.891 ( -3.543 11.485 11.485) 16.624 24.861 ( 23.006 11.313 11.313) 28.023 ======================= Grid point 369 (53/84) ======================= q-point: ( 0.38 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.256 ( 9.489 29.408 29.408) 42.658 3.352 ( 5.987 25.686 25.686) 36.816 4.709 ( 0.795 -3.438 -3.438) 4.926 5.265 ( -2.972 5.529 5.529) 8.365 6.010 ( 12.062 9.774 9.774) 18.345 6.490 ( 9.802 1.825 1.825) 10.136 6.638 ( 4.975 14.163 14.163) 20.638 7.390 ( 0.578 -18.917 -18.917) 26.758 8.494 ( 36.956 9.554 9.554) 39.349 9.459 ( -2.356 28.206 28.206) 39.959 10.828 ( 19.912 18.109 18.109) 32.440 12.691 ( 5.829 2.657 2.657) 6.935 15.374 ( -6.285 7.733 7.733) 12.614 15.801 ( -4.988 15.374 15.374) 22.307 25.257 ( 17.015 12.486 12.486) 24.522 ======================= Grid point 370 (54/84) ======================= q-point: ( 0.46 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.386 ( 3.606 31.113 31.113) 44.148 3.432 ( 2.182 26.008 26.008) 36.846 4.714 ( 0.029 -3.798 -3.798) 5.371 5.219 ( -1.442 2.567 2.567) 3.906 6.165 ( 3.971 7.431 7.431) 11.234 6.620 ( 3.507 1.737 1.737) 4.282 6.702 ( 1.665 14.515 14.515) 20.595 7.402 ( 0.446 -21.284 -21.284) 30.104 9.061 ( 21.828 17.035 17.035) 32.509 9.353 ( -9.359 21.102 21.102) 31.276 11.094 ( 7.236 17.303 17.303) 25.518 12.772 ( 2.245 2.901 2.901) 4.677 15.287 ( -2.450 8.209 8.209) 11.865 15.726 ( -2.282 17.183 17.183) 24.408 25.485 ( 6.195 13.353 13.353) 19.874 ======================= Grid point 380 (55/84) ======================= q-point: ( 0.23 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.068 ( 11.304 11.304 43.923) 46.741 3.456 ( 20.795 20.795 10.752) 31.313 4.626 ( 4.467 4.467 0.574) 6.343 5.506 ( 2.026 2.026 18.180) 18.405 5.692 ( 19.525 19.525 17.028) 32.441 6.047 ( 13.663 13.663 -5.737) 20.155 6.748 ( 18.413 18.413 2.753) 26.186 6.931 ( 9.675 9.675 0.984) 13.718 7.336 ( 4.747 4.747 -19.602) 20.720 9.910 ( 20.500 20.500 33.064) 43.974 10.133 ( 22.610 22.610 24.412) 40.229 12.509 ( 5.655 5.655 0.311) 8.004 15.621 ( -6.847 -6.847 20.766) 22.912 16.115 ( 11.105 11.105 -5.404) 16.609 24.729 ( 20.497 20.497 11.056) 31.024 ======================= Grid point 381 (56/84) ======================= q-point: ( 0.31 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.270 ( 9.449 8.709 50.295) 51.911 3.787 ( 14.787 23.690 10.091) 29.693 4.684 ( 1.764 2.764 -2.946) 4.409 5.461 ( -5.072 3.192 13.903) 15.139 6.022 ( 14.201 18.476 12.304) 26.352 6.307 ( 11.947 3.273 -3.945) 13.000 6.880 ( 2.190 8.461 -5.473) 10.312 7.102 ( 4.122 3.923 -6.534) 8.665 7.819 ( 39.055 7.894 -20.582) 44.846 10.121 ( 5.117 34.614 34.588) 49.201 10.652 ( 27.137 11.891 25.955) 39.388 12.635 ( 7.047 5.092 -1.072) 8.760 15.468 ( -8.432 -7.634 24.060) 26.613 16.269 ( 4.776 22.173 -4.784) 23.180 25.166 ( 23.429 20.012 11.006) 32.719 ======================= Grid point 382 (57/84) ======================= q-point: ( 0.38 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.423 ( 6.250 6.386 55.223) 55.942 4.048 ( 12.268 26.674 9.284) 30.793 4.696 ( -0.508 2.400 -2.588) 3.566 5.347 ( -6.023 2.562 8.824) 10.986 6.234 ( 7.974 13.439 7.946) 17.531 6.482 ( 6.013 -5.411 -2.133) 8.366 6.931 ( 3.400 -9.839 -11.725) 15.679 7.158 ( 2.080 24.339 -11.780) 27.120 8.549 ( 33.259 1.243 -19.997) 38.828 10.167 ( 0.368 38.166 37.491) 53.501 11.104 ( 19.022 10.767 23.892) 32.382 12.763 ( 5.862 4.741 -1.405) 7.669 15.315 ( -6.993 -8.890 26.588) 28.894 16.314 ( 0.515 31.087 -4.561) 31.424 25.576 ( 17.852 20.306 11.173) 29.255 ======================= Grid point 383 (58/84) ======================= q-point: ( 0.46 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.505 ( 2.195 5.044 57.678) 57.940 4.235 ( 6.001 31.634 9.114) 33.464 4.673 ( -1.341 -0.130 -0.558) 1.458 5.256 ( -2.801 1.250 4.581) 5.514 6.335 ( 2.565 10.636 6.376) 12.663 6.549 ( 1.446 -14.410 -1.424) 14.552 6.991 ( 2.051 -13.110 -16.046) 20.822 7.184 ( 0.639 32.513 -13.104) 35.060 9.012 ( 13.081 -3.147 -22.194) 25.953 10.156 ( -0.834 41.523 40.910) 58.297 11.358 ( 6.849 10.286 22.217) 25.422 12.845 ( 2.291 4.561 -1.346) 5.279 15.217 ( -2.783 -9.810 28.229) 30.014 16.311 ( -0.333 35.955 -4.816) 36.278 25.817 ( 6.577 20.743 11.376) 24.555 ======================= Grid point 394 (59/84) ======================= q-point: ( 0.31 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.424 ( 7.125 7.125 56.542) 57.433 4.131 ( 14.121 14.121 9.670) 22.188 4.746 ( 2.943 2.943 -2.830) 5.033 5.451 ( -3.404 -3.404 13.699) 14.521 6.324 ( -3.373 -3.373 -5.081) 6.969 6.347 ( 13.371 13.371 6.369) 19.953 6.729 ( -4.321 -4.321 -7.371) 9.575 7.612 ( 24.466 24.466 -15.638) 37.970 8.045 ( 20.269 20.269 -29.305) 40.993 10.685 ( 18.646 18.646 31.530) 41.103 10.866 ( 15.498 15.498 30.061) 37.202 12.749 ( 6.397 6.397 -3.211) 9.600 15.293 ( -9.855 -9.855 27.867) 31.158 16.641 ( 15.634 15.634 -8.025) 23.521 25.596 ( 23.269 23.269 10.151) 34.437 ======================= Grid point 395 (60/84) ======================= q-point: ( 0.38 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.538 ( 4.634 5.453 61.015) 61.433 4.380 ( 12.251 10.787 8.981) 18.631 4.779 ( 0.170 4.820 -0.203) 4.827 5.358 ( -5.534 -1.255 10.191) 11.665 6.185 ( -8.266 -21.705 -2.805) 23.395 6.515 ( 4.895 14.323 -1.095) 15.176 6.795 ( 7.346 -0.639 -13.027) 14.969 7.853 ( 5.468 39.127 -22.283) 45.358 8.591 ( 28.985 3.093 -27.268) 39.916 10.828 ( 2.439 29.841 32.546) 44.223 11.254 ( 18.509 5.297 27.153) 33.286 12.865 ( 5.283 5.518 -3.864) 8.561 15.112 ( -8.329 -11.403 30.075) 33.226 16.865 ( 7.935 24.697 -9.978) 27.793 26.007 ( 18.137 23.297 9.370) 30.976 ======================= Grid point 396 (61/84) ======================= q-point: ( 0.46 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.599 ( 1.614 4.587 63.184) 63.370 4.597 ( 9.793 7.478 7.613) 14.484 4.741 ( -4.486 7.105 4.938) 9.746 5.267 ( -3.030 -0.052 5.620) 6.384 6.068 ( -3.279 -29.466 -2.604) 29.762 6.572 ( 1.379 13.204 -3.938) 13.848 6.907 ( 3.299 0.432 -16.311) 16.647 7.914 ( 1.391 39.976 -24.516) 46.916 8.998 ( 11.441 1.028 -26.616) 28.989 10.850 ( 0.372 30.647 33.518) 45.418 11.503 ( 6.741 5.209 24.625) 26.057 12.938 ( 2.051 4.868 -3.902) 6.567 14.995 ( -3.348 -12.463 31.421) 33.967 16.960 ( 2.284 29.663 -11.196) 31.788 26.253 ( 6.761 23.467 8.967) 26.016 ======================= Grid point 408 (62/84) ======================= q-point: ( 0.38 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.627 ( 3.708 3.708 65.576) 65.786 4.552 ( 8.086 8.086 7.815) 13.851 4.843 ( 1.549 1.549 3.090) 3.788 5.313 ( -3.004 -3.004 8.380) 9.396 5.785 ( -17.965 -17.965 -2.266) 25.507 6.743 ( 7.545 7.545 -10.626) 15.059 6.831 ( 4.320 4.320 -16.029) 17.154 8.491 ( 19.736 19.736 -23.615) 36.561 8.705 ( 13.609 13.609 -30.538) 36.096 11.283 ( 12.036 12.036 27.316) 32.185 11.345 ( 9.587 9.587 28.373) 31.445 12.963 ( 4.316 4.316 -4.630) 7.661 14.902 ( -9.708 -9.708 31.864) 34.696 17.240 ( 14.209 14.209 -13.472) 24.192 26.420 ( 18.281 18.281 7.761) 26.994 ======================= Grid point 409 (63/84) ======================= q-point: ( 0.46 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.677 ( 1.330 3.363 68.071) 68.167 4.694 ( 5.827 4.647 6.234) 9.717 4.839 ( -1.658 2.466 7.277) 7.860 5.258 ( -2.106 -0.680 5.242) 5.690 5.540 ( -6.670 -23.920 -2.062) 24.918 6.805 ( 0.953 10.119 -14.514) 17.719 6.928 ( 3.349 1.162 -18.696) 19.029 8.621 ( 1.609 31.266 -27.292) 41.533 9.029 ( 10.929 1.860 -28.125) 30.231 11.337 ( 0.621 19.795 27.823) 34.151 11.573 ( 6.846 2.319 24.891) 25.919 13.022 ( 1.625 3.562 -4.687) 6.107 14.765 ( -3.931 -10.657 33.001) 34.901 17.418 ( 4.486 17.646 -15.428) 23.865 26.669 ( 6.855 18.405 6.870) 20.807 ======================= Grid point 422 (64/84) ======================= q-point: ( 0.46 0.46 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.722 ( 1.269 1.269 71.163) 71.186 4.790 ( 4.749 4.749 6.408) 9.283 4.855 ( -0.413 -0.413 6.968) 6.992 5.188 ( -10.966 -10.966 -2.105) 15.650 5.249 ( -0.529 -0.529 4.737) 4.796 6.937 ( 2.543 2.543 -19.748) 20.073 6.947 ( 1.547 1.547 -20.256) 20.374 9.038 ( 7.718 7.718 -27.860) 29.922 9.067 ( 4.859 4.859 -28.795) 29.604 11.597 ( 4.154 4.154 24.354) 25.052 11.603 ( 3.561 3.561 24.603) 25.114 13.069 ( 1.286 1.286 -4.725) 5.063 14.614 ( -4.338 -4.338 34.093) 34.640 17.641 ( 5.687 5.687 -17.738) 19.476 26.920 ( 6.912 6.912 5.699) 11.315 ======================= Grid point 549 (65/84) ======================= q-point: ( 0.23 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.702 ( 18.677 18.677 18.677) 32.349 3.702 ( 18.677 18.677 18.677) 32.349 4.655 ( 1.712 1.712 1.712) 2.965 5.912 ( 5.382 5.382 5.382) 9.323 5.912 ( 5.382 5.382 5.382) 9.323 6.057 ( 20.984 20.984 20.984) 36.346 6.942 ( 9.085 9.085 9.085) 15.737 6.942 ( 9.085 9.085 9.085) 15.737 7.022 ( -10.790 -10.790 -10.790) 18.688 10.532 ( 23.651 23.651 23.651) 40.965 10.532 ( 23.651 23.651 23.651) 40.965 12.543 ( 3.388 3.388 3.388) 5.868 15.994 ( 4.776 4.776 4.776) 8.273 15.994 ( 4.776 4.776 4.776) 8.273 25.031 ( 19.856 19.856 19.856) 34.392 ======================= Grid point 550 (66/84) ======================= q-point: ( 0.31 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.002 ( 16.352 21.449 21.449) 34.461 4.024 ( 12.344 17.992 17.992) 28.281 4.681 ( 0.897 2.168 2.168) 3.195 5.732 ( -10.285 10.732 10.732) 18.335 6.195 ( 11.770 -0.675 -0.675) 11.809 6.374 ( 10.701 16.606 16.606) 25.808 6.775 ( -8.782 -8.244 -8.244) 14.596 6.996 ( 1.332 2.122 2.122) 3.283 7.548 ( 38.856 -7.113 -7.113) 40.137 10.760 ( 7.766 30.725 30.725) 44.140 11.080 ( 25.402 18.090 18.090) 36.052 12.636 ( 5.950 1.514 1.514) 6.323 15.854 ( -7.278 10.228 10.228) 16.193 16.232 ( 8.159 5.198 5.198) 10.982 25.456 ( 22.955 18.785 18.785) 35.110 ======================= Grid point 551 (67/84) ======================= q-point: ( 0.38 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.213 ( 7.535 18.581 18.581) 27.337 4.314 ( 15.394 23.611 23.611) 36.769 4.690 ( 0.074 2.507 2.507) 3.546 5.524 ( -11.028 8.335 8.335) 16.141 6.383 ( 7.776 -2.540 -2.540) 8.565 6.510 ( 4.316 11.045 11.045) 16.206 6.647 ( -4.645 -16.499 -16.499) 23.791 7.021 ( 1.121 2.761 2.761) 4.062 8.261 ( 31.983 -9.544 -9.544) 34.714 10.854 ( 2.640 32.601 32.601) 46.180 11.491 ( 16.904 15.980 15.980) 28.222 12.752 ( 5.599 0.603 0.603) 5.664 15.726 ( -5.684 9.888 9.888) 15.094 16.324 ( 2.026 9.988 9.988) 14.270 25.860 ( 17.770 18.127 18.127) 31.192 ======================= Grid point 552 (68/84) ======================= q-point: ( 0.46 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.314 ( 2.799 19.266 19.266) 27.390 4.556 ( 8.051 27.928 27.928) 40.308 4.687 ( -0.200 2.368 2.368) 3.355 5.340 ( -6.471 3.967 3.967) 8.564 6.493 ( 3.421 -6.265 -6.265) 9.497 6.554 ( 0.580 10.392 10.392) 14.707 6.590 ( -1.402 -20.249 -20.249) 28.670 7.036 ( 0.418 2.488 2.488) 3.543 8.697 ( 11.990 -10.500 -10.500) 19.085 10.881 ( 0.497 33.533 33.533) 47.426 11.715 ( 5.977 14.630 14.630) 21.537 12.832 ( 2.275 0.303 0.303) 2.315 15.647 ( -2.187 9.722 9.722) 13.921 16.338 ( 0.026 12.447 12.447) 17.603 26.101 ( 6.602 17.915 17.915) 26.181 ======================= Grid point 563 (69/84) ======================= q-point: ( 0.31 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.285 ( 13.480 13.480 30.251) 35.757 4.354 ( 13.088 13.088 13.346) 22.819 4.740 ( 3.395 3.395 2.307) 5.327 5.778 ( -4.950 -4.950 18.318) 19.610 6.197 ( 0.198 0.198 -8.148) 8.152 6.549 ( -10.901 -10.901 -11.265) 19.094 6.574 ( 3.373 3.373 20.063) 20.622 7.350 ( 18.684 18.684 -11.145) 28.677 7.561 ( 16.874 16.874 -20.507) 31.464 11.248 ( 16.437 16.437 26.734) 35.427 11.379 ( 16.555 16.555 22.210) 32.271 12.696 ( 4.518 4.518 -1.822) 6.643 15.834 ( -6.754 -6.754 26.506) 28.175 16.463 ( 13.110 13.110 -10.195) 21.159 25.858 ( 21.860 21.860 16.929) 35.246 ======================= Grid point 564 (70/84) ======================= q-point: ( 0.38 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.478 ( 6.915 12.753 28.980) 32.408 4.612 ( 12.601 7.323 17.976) 23.142 4.809 ( 4.705 8.290 3.538) 10.168 5.616 ( -10.764 0.631 15.904) 19.215 6.119 ( -7.277 -20.952 -3.830) 22.508 6.462 ( 0.865 -2.766 -19.157) 19.375 6.599 ( 0.270 -1.653 11.251) 11.375 7.446 ( 1.220 35.476 -19.851) 40.670 8.157 ( 30.643 -2.239 -17.441) 35.330 11.403 ( 3.544 24.242 26.766) 36.286 11.718 ( 15.252 7.965 20.564) 26.813 12.792 ( 4.833 3.343 -3.371) 6.775 15.695 ( -6.691 -7.969 28.141) 30.004 16.653 ( 6.564 18.699 -11.292) 22.809 26.247 ( 17.201 20.965 15.348) 31.160 ======================= Grid point 565 (71/84) ======================= q-point: ( 0.46 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.569 ( 2.519 11.386 30.085) 32.266 4.743 ( 1.847 4.063 7.016) 8.315 4.978 ( 10.256 12.825 21.845) 27.329 5.398 ( -10.041 1.730 8.855) 13.499 6.005 ( -3.363 -29.129 -3.970) 29.590 6.496 ( 1.445 -1.548 -24.737) 24.827 6.598 ( -0.211 -4.144 9.226) 10.116 7.461 ( 0.401 35.336 -22.465) 41.875 8.578 ( 11.632 -2.857 -16.743) 20.586 11.439 ( 0.708 24.517 27.192) 36.620 11.921 ( 5.446 7.109 18.516) 20.568 12.862 ( 2.026 2.672 -3.718) 5.007 15.601 ( -2.687 -9.174 28.838) 30.382 16.728 ( 1.717 22.180 -11.765) 25.166 26.481 ( 6.453 20.558 14.524) 25.984 ======================= Grid point 577 (72/84) ======================= q-point: ( 0.38 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.689 ( 7.394 7.394 36.104) 37.588 4.729 ( 6.836 6.836 10.412) 14.208 4.952 ( 7.150 7.150 9.923) 14.168 5.556 ( -6.276 -6.276 18.597) 20.606 5.732 ( -17.498 -17.498 -3.135) 24.943 6.449 ( 0.720 0.720 -22.958) 22.981 6.555 ( -2.212 -2.212 -5.221) 6.087 8.085 ( 16.939 16.939 -18.126) 30.040 8.178 ( 13.993 13.993 -23.966) 31.080 11.763 ( 10.150 10.150 22.382) 26.590 11.838 ( 7.995 7.995 22.431) 25.119 12.866 ( 3.850 3.850 -5.141) 7.488 15.526 ( -8.314 -8.314 30.955) 33.113 16.944 ( 11.057 11.057 -16.722) 22.894 26.620 ( 16.628 16.628 12.931) 26.836 ======================= Grid point 578 (73/84) ======================= q-point: ( 0.46 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.775 ( 2.396 9.780 31.992) 33.539 4.831 ( 2.400 4.846 8.165) 9.793 5.114 ( 8.873 2.930 26.819) 28.400 5.392 ( -9.472 -3.315 11.426) 15.207 5.499 ( -5.942 -21.311 -1.589) 22.181 6.473 ( 1.296 -0.835 -26.581) 26.626 6.519 ( -1.397 -3.149 -13.150) 13.593 8.135 ( 0.336 30.822 -22.975) 38.444 8.555 ( 11.505 -0.116 -20.871) 23.832 11.825 ( 0.714 15.494 22.563) 27.380 12.007 ( 5.221 2.363 19.979) 20.785 12.922 ( 1.625 3.148 -5.426) 6.480 15.408 ( -3.389 -9.462 31.741) 33.295 17.081 ( 3.376 13.689 -18.938) 23.610 26.847 ( 6.278 16.331 11.598) 20.991 ======================= Grid point 591 (74/84) ======================= q-point: ( 0.46 0.46 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.919 ( 3.850 3.850 6.128) 8.197 4.931 ( 5.668 5.668 26.394) 27.584 5.131 ( -0.353 -0.353 45.616) 45.618 5.156 ( -10.516 -10.516 -0.360) 14.876 5.356 ( -2.333 -2.333 7.713) 8.389 6.466 ( 0.156 0.156 -28.947) 28.948 6.479 ( -1.206 -1.206 -26.338) 26.393 8.552 ( 6.503 6.503 -22.413) 24.226 8.563 ( 5.336 5.336 -23.227) 24.422 12.025 ( 3.419 3.419 19.880) 20.459 12.034 ( 2.453 2.453 20.008) 20.307 12.968 ( 1.322 1.322 -5.620) 5.923 15.274 ( -3.858 -3.858 32.427) 32.882 17.252 ( 4.332 4.332 -21.897) 22.738 27.071 ( 6.182 6.182 9.967) 13.258 ======================= Grid point 732 (75/84) ======================= q-point: ( 0.31 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.637 ( 13.052 13.052 13.052) 22.606 4.637 ( 13.052 13.052 13.052) 22.606 4.818 ( 4.590 4.590 4.590) 7.951 6.017 ( -2.173 -2.173 -2.173) 3.763 6.017 ( -2.173 -2.173 -2.173) 3.763 6.304 ( -13.739 -13.739 -13.739) 23.797 7.100 ( 14.166 14.166 14.166) 24.537 7.259 ( 3.612 3.612 3.612) 6.256 7.259 ( 3.612 3.612 3.612) 6.256 11.706 ( 16.203 16.203 16.203) 28.064 11.706 ( 16.203 16.203 16.203) 28.064 12.693 ( 1.693 1.693 1.693) 2.932 16.262 ( 3.554 3.554 3.554) 6.156 16.262 ( 3.554 3.554 3.554) 6.156 26.222 ( 19.832 19.832 19.832) 34.349 ======================= Grid point 733 (76/84) ======================= q-point: ( 0.38 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.784 ( 5.803 11.458 11.458) 17.211 4.888 ( 2.694 7.376 7.376) 10.774 4.955 ( 16.987 9.848 9.848) 21.967 5.817 ( -10.511 2.612 2.612) 11.141 6.042 ( -12.785 -14.081 -14.081) 23.664 6.056 ( -0.399 -10.865 -10.865) 15.370 7.085 ( -2.039 15.203 15.203) 21.597 7.104 ( -5.260 5.035 5.035) 8.853 7.909 ( 31.624 -8.920 -8.920) 34.047 11.851 ( 4.051 19.628 19.628) 28.053 12.033 ( 13.173 11.931 11.931) 21.406 12.752 ( 3.762 -0.404 -0.404) 3.805 16.085 ( -7.738 8.278 8.278) 14.033 16.502 ( 7.999 -0.625 -0.625) 8.048 26.576 ( 15.823 18.056 18.056) 30.040 ======================= Grid point 734 (77/84) ======================= q-point: ( 0.46 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.857 ( 1.692 7.166 7.166) 10.275 4.929 ( 1.468 11.781 11.781) 16.725 5.280 ( 16.376 7.289 7.289) 19.350 5.589 ( -13.733 7.304 7.304) 17.184 5.860 ( -4.920 -17.191 -17.191) 24.804 6.038 ( -0.827 -14.421 -14.421) 20.411 7.041 ( -1.528 1.665 1.665) 2.807 7.063 ( -0.461 13.365 13.365) 18.906 8.340 ( 11.814 -8.639 -8.639) 16.995 11.891 ( 0.771 19.811 19.811) 28.027 12.205 ( 4.567 10.725 10.725) 15.841 12.811 ( 1.830 -1.127 -1.127) 2.426 15.980 ( -2.910 7.071 7.071) 10.414 16.593 ( 2.050 0.815 0.815) 2.352 26.793 ( 5.990 17.099 17.099) 24.912 ======================= Grid point 746 (78/84) ======================= q-point: ( 0.38 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.943 ( 4.721 4.721 11.459) 13.262 5.004 ( 4.714 4.714 8.771) 11.017 5.156 ( 12.262 12.262 4.990) 18.045 5.661 ( -16.315 -16.315 -4.533) 23.514 5.786 ( -4.963 -4.963 -2.113) 7.330 5.989 ( -3.563 -3.563 -23.325) 23.863 6.831 ( -11.003 -11.003 34.835) 38.153 7.782 ( 13.694 13.694 -17.333) 25.990 7.797 ( 14.064 14.064 -11.912) 23.183 12.140 ( 8.479 8.479 16.660) 20.527 12.199 ( 7.780 7.780 14.829) 18.465 12.777 ( 2.587 2.587 -3.675) 5.186 16.037 ( -7.165 -7.165 21.349) 23.632 16.620 ( 7.731 7.731 -16.306) 19.632 26.901 ( 14.567 14.567 15.523) 25.794 ======================= Grid point 747 (79/84) ======================= q-point: ( 0.46 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.976 ( 0.145 4.582 7.385) 8.693 5.099 ( 3.909 5.004 11.650) 13.268 5.356 ( 3.490 -1.832 1.139) 4.103 5.456 ( -0.929 -13.482 -1.064) 13.556 5.676 ( -5.810 4.006 -3.234) 7.763 5.923 ( -2.479 -2.853 -27.893) 28.148 6.686 ( -3.882 -17.331 34.775) 39.048 7.741 ( -1.145 31.689 -17.803) 36.365 8.223 ( 11.879 -4.005 -13.965) 18.766 12.201 ( 0.769 12.282 16.442) 20.537 12.335 ( 4.081 3.335 13.987) 14.947 12.823 ( 1.522 2.010 -4.522) 5.178 15.927 ( -3.295 -7.447 21.120) 22.636 16.720 ( 2.517 8.142 -17.541) 19.502 27.101 ( 5.557 13.867 14.113) 20.551 ======================= Grid point 760 (80/84) ======================= q-point: ( 0.46 0.46 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.976 ( -3.417 -3.417 -0.289) 4.841 5.156 ( 2.692 2.692 8.687) 9.485 5.211 ( -3.785 -3.785 5.356) 7.572 5.440 ( 0.667 0.667 4.807) 4.899 5.786 ( 6.116 6.116 -18.404) 20.336 5.884 ( -1.543 -1.543 -29.608) 29.688 6.412 ( -9.403 -9.403 44.321) 46.273 8.173 ( 5.542 5.542 -17.472) 19.149 8.176 ( 5.186 5.186 -16.715) 18.253 12.357 ( 2.821 2.821 14.580) 15.116 12.368 ( 1.790 1.790 14.566) 14.785 12.861 ( 1.315 1.315 -5.203) 5.525 15.808 ( -3.622 -3.622 22.375) 22.954 16.830 ( 2.875 2.875 -21.171) 21.558 27.292 ( 5.310 5.310 12.441) 14.532 ======================= Grid point 915 (81/84) ======================= q-point: ( 0.38 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.155 ( 3.961 3.961 3.961) 6.860 5.155 ( 3.961 3.961 3.961) 6.860 5.166 ( 6.880 6.880 6.880) 11.917 5.503 ( -13.123 -13.123 -13.123) 22.729 5.657 ( -8.042 -8.042 -8.042) 13.929 5.657 ( -8.042 -8.042 -8.042) 13.929 7.507 ( 4.849 4.849 4.849) 8.399 7.507 ( 4.849 4.849 4.849) 8.399 7.653 ( 5.089 5.089 5.089) 8.814 12.401 ( 8.220 8.220 8.220) 14.238 12.401 ( 8.220 8.220 8.220) 14.238 12.741 ( 0.137 0.137 0.137) 0.237 16.342 ( -0.645 -0.645 -0.645) 1.117 16.342 ( -0.645 -0.645 -0.645) 1.117 27.182 ( 12.749 12.749 12.749) 22.082 ======================= Grid point 916 (82/84) ======================= q-point: ( 0.46 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 5.087 ( -2.506 2.196 2.196) 3.991 5.240 ( -3.329 -10.590 -10.590) 15.342 5.301 ( 5.116 7.108 7.108) 11.279 5.378 ( 3.923 -9.271 -9.271) 13.686 5.541 ( -4.213 -7.584 -7.584) 11.523 5.542 ( -4.720 -2.972 -2.972) 6.320 7.357 ( -4.073 12.166 12.166) 17.681 7.510 ( -2.063 9.246 9.246) 13.237 8.009 ( 11.988 -8.372 -8.372) 16.849 12.456 ( 0.792 10.063 10.063) 14.253 12.535 ( 3.527 6.400 6.400) 9.714 12.764 ( 1.221 -1.133 -1.133) 2.014 16.181 ( -4.409 3.502 3.502) 6.631 16.468 ( 2.965 -5.840 -5.840) 8.775 27.358 ( 4.914 11.724 11.724) 17.293 ======================= Grid point 929 (83/84) ======================= q-point: ( 0.46 0.46 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.950 ( -7.467 -7.467 -1.685) 10.694 5.190 ( -1.823 -1.823 -12.577) 12.838 5.194 ( -2.494 -2.494 -13.771) 14.216 5.442 ( 4.771 4.771 -2.106) 7.068 5.502 ( -1.012 -1.012 -3.036) 3.357 5.519 ( -0.091 -0.091 2.545) 2.548 7.256 ( -4.576 -4.576 36.416) 36.987 7.890 ( 5.584 5.584 -11.745) 14.153 7.908 ( 3.807 3.807 -11.075) 12.314 12.584 ( 2.366 2.366 8.896) 9.504 12.593 ( 1.575 1.575 8.759) 9.037 12.775 ( 1.151 1.151 -3.559) 3.914 16.129 ( -3.486 -3.486 11.327) 12.353 16.470 ( 1.665 1.665 -15.365) 15.545 27.521 ( 4.548 4.548 10.512) 12.323 ======================= Grid point 1098 (84/84) ======================= q-point: ( 0.46 0.46 0.46) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.36e-05 3.36e-05 3.36e-05 1.64e-04 3.67e-09 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.921 ( -5.815 -5.815 -5.815) 10.072 4.928 ( -5.973 -5.973 -5.973) 10.346 4.928 ( -5.973 -5.973 -5.973) 10.346 5.508 ( 3.810 3.810 3.810) 6.600 5.546 ( 1.193 1.193 1.193) 2.066 5.546 ( 1.193 1.193 1.193) 2.066 7.732 ( 2.275 2.275 2.275) 3.941 7.732 ( 2.275 2.275 2.275) 3.941 7.767 ( 0.049 0.049 0.049) 0.084 12.698 ( 2.402 2.402 2.402) 4.160 12.698 ( 2.402 2.402 2.402) 4.160 12.735 ( -0.152 -0.152 -0.152) 0.264 16.264 ( -1.215 -1.215 -1.215) 2.105 16.264 ( -1.215 -1.215 -1.215) 2.105 27.667 ( 4.113 4.113 4.113) 7.125 =================== 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/32955 10.0 755.946 755.946 755.946 0.000 -0.000 0.000 3/32955 20.0 144.757 144.757 144.757 0.000 -0.000 0.000 3/32955 30.0 85.125 85.125 85.125 0.000 -0.000 0.000 3/32955 40.0 60.227 60.227 60.227 0.000 -0.000 0.000 3/32955 50.0 46.640 46.640 46.640 0.000 -0.000 0.000 3/32955 60.0 38.379 38.379 38.379 0.000 -0.000 0.000 3/32955 70.0 32.935 32.935 32.935 0.000 -0.000 0.000 3/32955 80.0 29.089 29.089 29.089 0.000 -0.000 0.000 3/32955 90.0 26.206 26.206 26.206 0.000 -0.000 0.000 3/32955 100.0 23.937 23.937 23.937 0.000 -0.000 0.000 3/32955 110.0 22.086 22.086 22.086 0.000 -0.000 0.000 3/32955 120.0 20.531 20.531 20.531 0.000 -0.000 0.000 3/32955 130.0 19.199 19.199 19.199 0.000 -0.000 0.000 3/32955 140.0 18.040 18.040 18.040 0.000 -0.000 0.000 3/32955 150.0 17.019 17.019 17.019 0.000 -0.000 0.000 3/32955 160.0 16.112 16.112 16.112 0.000 -0.000 0.000 3/32955 170.0 15.301 15.301 15.301 0.000 -0.000 0.000 3/32955 180.0 14.570 14.570 14.570 0.000 -0.000 0.000 3/32955 190.0 13.908 13.908 13.908 0.000 -0.000 0.000 3/32955 200.0 13.306 13.306 13.306 0.000 -0.000 0.000 3/32955 210.0 12.756 12.756 12.756 0.000 -0.000 0.000 3/32955 220.0 12.251 12.251 12.251 0.000 -0.000 0.000 3/32955 230.0 11.787 11.787 11.787 0.000 -0.000 0.000 3/32955 240.0 11.358 11.358 11.358 0.000 -0.000 0.000 3/32955 250.0 10.960 10.960 10.960 0.000 -0.000 0.000 3/32955 260.0 10.591 10.591 10.591 0.000 -0.000 0.000 3/32955 270.0 10.247 10.247 10.247 0.000 -0.000 0.000 3/32955 280.0 9.925 9.925 9.925 0.000 -0.000 0.000 3/32955 290.0 9.624 9.624 9.624 0.000 -0.000 0.000 3/32955 300.0 9.342 9.342 9.342 0.000 -0.000 0.000 3/32955 310.0 9.076 9.076 9.076 0.000 -0.000 0.000 3/32955 320.0 8.825 8.825 8.825 0.000 -0.000 0.000 3/32955 330.0 8.588 8.588 8.588 0.000 -0.000 0.000 3/32955 340.0 8.364 8.364 8.364 0.000 -0.000 0.000 3/32955 350.0 8.151 8.151 8.151 0.000 -0.000 0.000 3/32955 360.0 7.950 7.950 7.950 0.000 -0.000 0.000 3/32955 370.0 7.758 7.758 7.758 0.000 -0.000 0.000 3/32955 380.0 7.576 7.576 7.576 0.000 -0.000 0.000 3/32955 390.0 7.402 7.402 7.402 0.000 -0.000 0.000 3/32955 400.0 7.235 7.235 7.235 0.000 -0.000 0.000 3/32955 410.0 7.077 7.077 7.077 0.000 -0.000 0.000 3/32955 420.0 6.925 6.925 6.925 0.000 -0.000 0.000 3/32955 430.0 6.779 6.779 6.779 0.000 -0.000 0.000 3/32955 440.0 6.640 6.640 6.640 0.000 -0.000 0.000 3/32955 450.0 6.506 6.506 6.506 0.000 -0.000 0.000 3/32955 460.0 6.377 6.377 6.377 0.000 -0.000 0.000 3/32955 470.0 6.254 6.254 6.254 0.000 -0.000 0.000 3/32955 480.0 6.135 6.135 6.135 0.000 -0.000 0.000 3/32955 490.0 6.021 6.021 6.021 0.000 -0.000 0.000 3/32955 500.0 5.910 5.910 5.910 0.000 -0.000 0.000 3/32955 510.0 5.804 5.804 5.804 0.000 -0.000 0.000 3/32955 520.0 5.701 5.701 5.701 0.000 -0.000 0.000 3/32955 530.0 5.602 5.602 5.602 0.000 -0.000 0.000 3/32955 540.0 5.506 5.506 5.506 0.000 -0.000 0.000 3/32955 550.0 5.414 5.414 5.414 0.000 -0.000 0.000 3/32955 560.0 5.324 5.324 5.324 0.000 -0.000 0.000 3/32955 570.0 5.238 5.238 5.238 0.000 -0.000 0.000 3/32955 580.0 5.154 5.154 5.154 0.000 -0.000 0.000 3/32955 590.0 5.073 5.073 5.073 0.000 -0.000 0.000 3/32955 600.0 4.994 4.994 4.994 0.000 -0.000 0.000 3/32955 610.0 4.917 4.917 4.917 0.000 -0.000 0.000 3/32955 620.0 4.843 4.843 4.843 0.000 -0.000 0.000 3/32955 630.0 4.771 4.771 4.771 0.000 -0.000 0.000 3/32955 640.0 4.701 4.701 4.701 0.000 -0.000 0.000 3/32955 650.0 4.633 4.633 4.633 0.000 -0.000 0.000 3/32955 660.0 4.567 4.567 4.567 0.000 -0.000 0.000 3/32955 670.0 4.503 4.503 4.503 0.000 -0.000 0.000 3/32955 680.0 4.441 4.441 4.441 0.000 -0.000 0.000 3/32955 690.0 4.380 4.380 4.380 0.000 -0.000 0.000 3/32955 700.0 4.321 4.321 4.321 0.000 -0.000 0.000 3/32955 710.0 4.263 4.263 4.263 0.000 -0.000 0.000 3/32955 720.0 4.207 4.207 4.207 0.000 -0.000 0.000 3/32955 730.0 4.153 4.153 4.153 0.000 -0.000 0.000 3/32955 740.0 4.099 4.099 4.099 0.000 -0.000 0.000 3/32955 750.0 4.048 4.048 4.048 0.000 0.000 0.000 3/32955 760.0 3.997 3.997 3.997 0.000 0.000 0.000 3/32955 770.0 3.947 3.947 3.947 0.000 0.000 0.000 3/32955 780.0 3.899 3.899 3.899 0.000 0.000 0.000 3/32955 790.0 3.852 3.852 3.852 0.000 0.000 0.000 3/32955 800.0 3.806 3.806 3.806 0.000 0.000 0.000 3/32955 810.0 3.761 3.761 3.761 0.000 0.000 0.000 3/32955 820.0 3.717 3.717 3.717 0.000 0.000 0.000 3/32955 830.0 3.675 3.675 3.675 0.000 0.000 0.000 3/32955 840.0 3.633 3.633 3.633 0.000 0.000 0.000 3/32955 850.0 3.592 3.592 3.592 0.000 0.000 0.000 3/32955 860.0 3.552 3.552 3.552 0.000 0.000 0.000 3/32955 870.0 3.512 3.512 3.512 0.000 0.000 0.000 3/32955 880.0 3.474 3.474 3.474 0.000 0.000 0.000 3/32955 890.0 3.437 3.437 3.437 0.000 0.000 0.000 3/32955 900.0 3.400 3.400 3.400 0.000 0.000 0.000 3/32955 910.0 3.364 3.364 3.364 0.000 0.000 0.000 3/32955 920.0 3.329 3.329 3.329 0.000 0.000 0.000 3/32955 930.0 3.294 3.294 3.294 0.000 0.000 0.000 3/32955 940.0 3.260 3.260 3.260 0.000 0.000 0.000 3/32955 950.0 3.227 3.227 3.227 0.000 0.000 0.000 3/32955 960.0 3.195 3.195 3.195 0.000 0.000 0.000 3/32955 970.0 3.163 3.163 3.163 0.000 0.000 0.000 3/32955 980.0 3.132 3.132 3.132 0.000 0.000 0.000 3/32955 990.0 3.101 3.101 3.101 0.000 0.000 0.000 3/32955 1000.0 3.071 3.071 3.071 0.000 0.000 0.000 3/32955 Thermal conductivity related properties were written into "kappa-m131313.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 19:39:09]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|