----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 05:35:04]------------------------- 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.989339610000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.989339610000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.989339610000000 Atomic positions (fractional): *1 H 0.50000000000000 0.00000000000000 0.50000000000000 1.008 2 H 0.50000000000000 0.50000000000000 0.00000000000000 1.008 3 H 0.00000000000000 0.50000000000000 0.50000000000000 1.008 *4 Mg 0.50000000000000 0.50000000000000 0.50000000000000 24.305 *5 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.989339610000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.989339610000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.989339610000000 Atomic positions (fractional): *1 H 0.50000000000000 0.00000000000000 0.50000000000000 1.008 > 1 2 H 0.50000000000000 0.50000000000000 0.00000000000000 1.008 > 2 3 H 0.00000000000000 0.50000000000000 0.50000000000000 1.008 > 3 *4 Mg 0.50000000000000 0.50000000000000 0.50000000000000 24.305 > 4 *5 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 > 5 -------------------------------- super cell -------------------------------- Lattice vectors: a 11.968018830000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.968018830000000 0.000000000000000 c 0.000000000000000 0.000000000000000 11.968018830000000 Atomic positions (fractional): *1 H 0.16666666666667 0.00000000000000 0.16666666666667 1.008 > 1 2 H 0.50000000000000 0.00000000000000 0.16666666666667 1.008 > 1 3 H 0.83333333333333 0.00000000000000 0.16666666666667 1.008 > 1 4 H 0.16666666666667 0.33333333333333 0.16666666666667 1.008 > 1 5 H 0.50000000000000 0.33333333333333 0.16666666666667 1.008 > 1 6 H 0.83333333333333 0.33333333333333 0.16666666666667 1.008 > 1 7 H 0.16666666666667 0.66666666666667 0.16666666666667 1.008 > 1 8 H 0.50000000000000 0.66666666666667 0.16666666666667 1.008 > 1 9 H 0.83333333333333 0.66666666666667 0.16666666666667 1.008 > 1 10 H 0.16666666666667 0.00000000000000 0.50000000000000 1.008 > 1 11 H 0.50000000000000 0.00000000000000 0.50000000000000 1.008 > 1 12 H 0.83333333333333 0.00000000000000 0.50000000000000 1.008 > 1 13 H 0.16666666666667 0.33333333333333 0.50000000000000 1.008 > 1 14 H 0.50000000000000 0.33333333333333 0.50000000000000 1.008 > 1 15 H 0.83333333333333 0.33333333333333 0.50000000000000 1.008 > 1 16 H 0.16666666666667 0.66666666666667 0.50000000000000 1.008 > 1 17 H 0.50000000000000 0.66666666666667 0.50000000000000 1.008 > 1 18 H 0.83333333333333 0.66666666666667 0.50000000000000 1.008 > 1 19 H 0.16666666666667 0.00000000000000 0.83333333333333 1.008 > 1 20 H 0.50000000000000 0.00000000000000 0.83333333333333 1.008 > 1 21 H 0.83333333333333 0.00000000000000 0.83333333333333 1.008 > 1 22 H 0.16666666666667 0.33333333333333 0.83333333333333 1.008 > 1 23 H 0.50000000000000 0.33333333333333 0.83333333333333 1.008 > 1 24 H 0.83333333333333 0.33333333333333 0.83333333333333 1.008 > 1 25 H 0.16666666666667 0.66666666666667 0.83333333333333 1.008 > 1 26 H 0.50000000000000 0.66666666666667 0.83333333333333 1.008 > 1 27 H 0.83333333333333 0.66666666666667 0.83333333333333 1.008 > 1 28 H 0.16666666666667 0.16666666666667 0.00000000000000 1.008 > 2 29 H 0.50000000000000 0.16666666666667 0.00000000000000 1.008 > 2 30 H 0.83333333333333 0.16666666666667 0.00000000000000 1.008 > 2 31 H 0.16666666666667 0.50000000000000 0.00000000000000 1.008 > 2 32 H 0.50000000000000 0.50000000000000 0.00000000000000 1.008 > 2 33 H 0.83333333333333 0.50000000000000 0.00000000000000 1.008 > 2 34 H 0.16666666666667 0.83333333333333 0.00000000000000 1.008 > 2 35 H 0.50000000000000 0.83333333333333 0.00000000000000 1.008 > 2 36 H 0.83333333333333 0.83333333333333 0.00000000000000 1.008 > 2 37 H 0.16666666666667 0.16666666666667 0.33333333333333 1.008 > 2 38 H 0.50000000000000 0.16666666666667 0.33333333333333 1.008 > 2 39 H 0.83333333333333 0.16666666666667 0.33333333333333 1.008 > 2 40 H 0.16666666666667 0.50000000000000 0.33333333333333 1.008 > 2 41 H 0.50000000000000 0.50000000000000 0.33333333333333 1.008 > 2 42 H 0.83333333333333 0.50000000000000 0.33333333333333 1.008 > 2 43 H 0.16666666666667 0.83333333333333 0.33333333333333 1.008 > 2 44 H 0.50000000000000 0.83333333333333 0.33333333333333 1.008 > 2 45 H 0.83333333333333 0.83333333333333 0.33333333333333 1.008 > 2 46 H 0.16666666666667 0.16666666666667 0.66666666666667 1.008 > 2 47 H 0.50000000000000 0.16666666666667 0.66666666666667 1.008 > 2 48 H 0.83333333333333 0.16666666666667 0.66666666666667 1.008 > 2 49 H 0.16666666666667 0.50000000000000 0.66666666666667 1.008 > 2 50 H 0.50000000000000 0.50000000000000 0.66666666666667 1.008 > 2 51 H 0.83333333333333 0.50000000000000 0.66666666666667 1.008 > 2 52 H 0.16666666666667 0.83333333333333 0.66666666666667 1.008 > 2 53 H 0.50000000000000 0.83333333333333 0.66666666666667 1.008 > 2 54 H 0.83333333333333 0.83333333333333 0.66666666666667 1.008 > 2 55 H 0.00000000000000 0.16666666666667 0.16666666666667 1.008 > 3 56 H 0.33333333333333 0.16666666666667 0.16666666666667 1.008 > 3 57 H 0.66666666666667 0.16666666666667 0.16666666666667 1.008 > 3 58 H 0.00000000000000 0.50000000000000 0.16666666666667 1.008 > 3 59 H 0.33333333333333 0.50000000000000 0.16666666666667 1.008 > 3 60 H 0.66666666666667 0.50000000000000 0.16666666666667 1.008 > 3 61 H 0.00000000000000 0.83333333333333 0.16666666666667 1.008 > 3 62 H 0.33333333333333 0.83333333333333 0.16666666666667 1.008 > 3 63 H 0.66666666666667 0.83333333333333 0.16666666666667 1.008 > 3 64 H 0.00000000000000 0.16666666666667 0.50000000000000 1.008 > 3 65 H 0.33333333333333 0.16666666666667 0.50000000000000 1.008 > 3 66 H 0.66666666666667 0.16666666666667 0.50000000000000 1.008 > 3 67 H 0.00000000000000 0.50000000000000 0.50000000000000 1.008 > 3 68 H 0.33333333333333 0.50000000000000 0.50000000000000 1.008 > 3 69 H 0.66666666666667 0.50000000000000 0.50000000000000 1.008 > 3 70 H 0.00000000000000 0.83333333333333 0.50000000000000 1.008 > 3 71 H 0.33333333333333 0.83333333333333 0.50000000000000 1.008 > 3 72 H 0.66666666666667 0.83333333333333 0.50000000000000 1.008 > 3 73 H 0.00000000000000 0.16666666666667 0.83333333333333 1.008 > 3 74 H 0.33333333333333 0.16666666666667 0.83333333333333 1.008 > 3 75 H 0.66666666666667 0.16666666666667 0.83333333333333 1.008 > 3 76 H 0.00000000000000 0.50000000000000 0.83333333333333 1.008 > 3 77 H 0.33333333333333 0.50000000000000 0.83333333333333 1.008 > 3 78 H 0.66666666666667 0.50000000000000 0.83333333333333 1.008 > 3 79 H 0.00000000000000 0.83333333333333 0.83333333333333 1.008 > 3 80 H 0.33333333333333 0.83333333333333 0.83333333333333 1.008 > 3 81 H 0.66666666666667 0.83333333333333 0.83333333333333 1.008 > 3 *82 Mg 0.16666666666667 0.16666666666667 0.16666666666667 24.305 > 4 83 Mg 0.50000000000000 0.16666666666667 0.16666666666667 24.305 > 4 84 Mg 0.83333333333333 0.16666666666667 0.16666666666667 24.305 > 4 85 Mg 0.16666666666667 0.50000000000000 0.16666666666667 24.305 > 4 86 Mg 0.50000000000000 0.50000000000000 0.16666666666667 24.305 > 4 87 Mg 0.83333333333333 0.50000000000000 0.16666666666667 24.305 > 4 88 Mg 0.16666666666667 0.83333333333333 0.16666666666667 24.305 > 4 89 Mg 0.50000000000000 0.83333333333333 0.16666666666667 24.305 > 4 90 Mg 0.83333333333333 0.83333333333333 0.16666666666667 24.305 > 4 91 Mg 0.16666666666667 0.16666666666667 0.50000000000000 24.305 > 4 92 Mg 0.50000000000000 0.16666666666667 0.50000000000000 24.305 > 4 93 Mg 0.83333333333333 0.16666666666667 0.50000000000000 24.305 > 4 94 Mg 0.16666666666667 0.50000000000000 0.50000000000000 24.305 > 4 95 Mg 0.50000000000000 0.50000000000000 0.50000000000000 24.305 > 4 96 Mg 0.83333333333333 0.50000000000000 0.50000000000000 24.305 > 4 97 Mg 0.16666666666667 0.83333333333333 0.50000000000000 24.305 > 4 98 Mg 0.50000000000000 0.83333333333333 0.50000000000000 24.305 > 4 99 Mg 0.83333333333333 0.83333333333333 0.50000000000000 24.305 > 4 100 Mg 0.16666666666667 0.16666666666667 0.83333333333333 24.305 > 4 101 Mg 0.50000000000000 0.16666666666667 0.83333333333333 24.305 > 4 102 Mg 0.83333333333333 0.16666666666667 0.83333333333333 24.305 > 4 103 Mg 0.16666666666667 0.50000000000000 0.83333333333333 24.305 > 4 104 Mg 0.50000000000000 0.50000000000000 0.83333333333333 24.305 > 4 105 Mg 0.83333333333333 0.50000000000000 0.83333333333333 24.305 > 4 106 Mg 0.16666666666667 0.83333333333333 0.83333333333333 24.305 > 4 107 Mg 0.50000000000000 0.83333333333333 0.83333333333333 24.305 > 4 108 Mg 0.83333333333333 0.83333333333333 0.83333333333333 24.305 > 4 *109 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 > 5 110 K 0.33333333333333 0.00000000000000 0.00000000000000 39.098 > 5 111 K 0.66666666666667 0.00000000000000 0.00000000000000 39.098 > 5 112 K 0.00000000000000 0.33333333333333 0.00000000000000 39.098 > 5 113 K 0.33333333333333 0.33333333333333 0.00000000000000 39.098 > 5 114 K 0.66666666666667 0.33333333333333 0.00000000000000 39.098 > 5 115 K 0.00000000000000 0.66666666666667 0.00000000000000 39.098 > 5 116 K 0.33333333333333 0.66666666666667 0.00000000000000 39.098 > 5 117 K 0.66666666666667 0.66666666666667 0.00000000000000 39.098 > 5 118 K 0.00000000000000 0.00000000000000 0.33333333333333 39.098 > 5 119 K 0.33333333333333 0.00000000000000 0.33333333333333 39.098 > 5 120 K 0.66666666666667 0.00000000000000 0.33333333333333 39.098 > 5 121 K 0.00000000000000 0.33333333333333 0.33333333333333 39.098 > 5 122 K 0.33333333333333 0.33333333333333 0.33333333333333 39.098 > 5 123 K 0.66666666666667 0.33333333333333 0.33333333333333 39.098 > 5 124 K 0.00000000000000 0.66666666666667 0.33333333333333 39.098 > 5 125 K 0.33333333333333 0.66666666666667 0.33333333333333 39.098 > 5 126 K 0.66666666666667 0.66666666666667 0.33333333333333 39.098 > 5 127 K 0.00000000000000 0.00000000000000 0.66666666666667 39.098 > 5 128 K 0.33333333333333 0.00000000000000 0.66666666666667 39.098 > 5 129 K 0.66666666666667 0.00000000000000 0.66666666666667 39.098 > 5 130 K 0.00000000000000 0.33333333333333 0.66666666666667 39.098 > 5 131 K 0.33333333333333 0.33333333333333 0.66666666666667 39.098 > 5 132 K 0.66666666666667 0.33333333333333 0.66666666666667 39.098 > 5 133 K 0.00000000000000 0.66666666666667 0.66666666666667 39.098 > 5 134 K 0.33333333333333 0.66666666666667 0.66666666666667 39.098 > 5 135 K 0.66666666666667 0.66666666666667 0.66666666666667 39.098 > 5 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.7023609 0.0000000 0.0000000 0.0000000 3.7023609 0.0000000 0.0000000 0.0000000 3.7023609 -------------------------- Born effective charges -------------------------- 1 H -0.8075801 0.0000000 0.0000000 0.0000000 -1.2289746 0.0000000 0.0000000 0.0000000 -0.8075801 2 H -0.8075801 0.0000000 0.0000000 0.0000000 -0.8075801 0.0000000 0.0000000 0.0000000 -1.2289746 3 H -1.2289746 0.0000000 0.0000000 0.0000000 -0.8075801 0.0000000 0.0000000 0.0000000 -0.8075801 4 Mg 1.8670214 0.0000000 0.0000000 0.0000000 1.8670214 0.0000000 0.0000000 0.0000000 1.8670214 5 K 0.9771135 0.0000000 0.0000000 0.0000000 0.9771135 0.0000000 0.0000000 0.0000000 0.9771135 ---------------------------------------------------------------------------- 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.005 Solver_block: 80 / 80 - Time: 0.063 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.068 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 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 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 05:35:09]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 05:35:09]------------------------- 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.989339610000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.989339610000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.989339610000000 Atomic positions (fractional): 1 H 0.50000000000000 0.00000000000000 0.50000000000000 1.008 2 H 0.50000000000000 0.50000000000000 0.00000000000000 1.008 3 H 0.00000000000000 0.50000000000000 0.50000000000000 1.008 4 Mg 0.50000000000000 0.50000000000000 0.50000000000000 24.305 5 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.968018830000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.968018830000000 0.000000000000000 c 0.000000000000000 0.000000000000000 11.968018830000000 Atomic positions (fractional): 1 H 0.16666666666667 0.00000000000000 0.16666666666667 1.008 > 1 2 H 0.50000000000000 0.00000000000000 0.16666666666667 1.008 > 1 3 H 0.83333333333333 0.00000000000000 0.16666666666667 1.008 > 1 4 H 0.16666666666667 0.33333333333333 0.16666666666667 1.008 > 1 5 H 0.50000000000000 0.33333333333333 0.16666666666667 1.008 > 1 6 H 0.83333333333333 0.33333333333333 0.16666666666667 1.008 > 1 7 H 0.16666666666667 0.66666666666667 0.16666666666667 1.008 > 1 8 H 0.50000000000000 0.66666666666667 0.16666666666667 1.008 > 1 9 H 0.83333333333333 0.66666666666667 0.16666666666667 1.008 > 1 10 H 0.16666666666667 0.00000000000000 0.50000000000000 1.008 > 1 11 H 0.50000000000000 0.00000000000000 0.50000000000000 1.008 > 1 12 H 0.83333333333333 0.00000000000000 0.50000000000000 1.008 > 1 13 H 0.16666666666667 0.33333333333333 0.50000000000000 1.008 > 1 14 H 0.50000000000000 0.33333333333333 0.50000000000000 1.008 > 1 15 H 0.83333333333333 0.33333333333333 0.50000000000000 1.008 > 1 16 H 0.16666666666667 0.66666666666667 0.50000000000000 1.008 > 1 17 H 0.50000000000000 0.66666666666667 0.50000000000000 1.008 > 1 18 H 0.83333333333333 0.66666666666667 0.50000000000000 1.008 > 1 19 H 0.16666666666667 0.00000000000000 0.83333333333333 1.008 > 1 20 H 0.50000000000000 0.00000000000000 0.83333333333333 1.008 > 1 21 H 0.83333333333333 0.00000000000000 0.83333333333333 1.008 > 1 22 H 0.16666666666667 0.33333333333333 0.83333333333333 1.008 > 1 23 H 0.50000000000000 0.33333333333333 0.83333333333333 1.008 > 1 24 H 0.83333333333333 0.33333333333333 0.83333333333333 1.008 > 1 25 H 0.16666666666667 0.66666666666667 0.83333333333333 1.008 > 1 26 H 0.50000000000000 0.66666666666667 0.83333333333333 1.008 > 1 27 H 0.83333333333333 0.66666666666667 0.83333333333333 1.008 > 1 28 H 0.16666666666667 0.16666666666667 0.00000000000000 1.008 > 28 29 H 0.50000000000000 0.16666666666667 0.00000000000000 1.008 > 28 30 H 0.83333333333333 0.16666666666667 0.00000000000000 1.008 > 28 31 H 0.16666666666667 0.50000000000000 0.00000000000000 1.008 > 28 32 H 0.50000000000000 0.50000000000000 0.00000000000000 1.008 > 28 33 H 0.83333333333333 0.50000000000000 0.00000000000000 1.008 > 28 34 H 0.16666666666667 0.83333333333333 0.00000000000000 1.008 > 28 35 H 0.50000000000000 0.83333333333333 0.00000000000000 1.008 > 28 36 H 0.83333333333333 0.83333333333333 0.00000000000000 1.008 > 28 37 H 0.16666666666667 0.16666666666667 0.33333333333333 1.008 > 28 38 H 0.50000000000000 0.16666666666667 0.33333333333333 1.008 > 28 39 H 0.83333333333333 0.16666666666667 0.33333333333333 1.008 > 28 40 H 0.16666666666667 0.50000000000000 0.33333333333333 1.008 > 28 41 H 0.50000000000000 0.50000000000000 0.33333333333333 1.008 > 28 42 H 0.83333333333333 0.50000000000000 0.33333333333333 1.008 > 28 43 H 0.16666666666667 0.83333333333333 0.33333333333333 1.008 > 28 44 H 0.50000000000000 0.83333333333333 0.33333333333333 1.008 > 28 45 H 0.83333333333333 0.83333333333333 0.33333333333333 1.008 > 28 46 H 0.16666666666667 0.16666666666667 0.66666666666667 1.008 > 28 47 H 0.50000000000000 0.16666666666667 0.66666666666667 1.008 > 28 48 H 0.83333333333333 0.16666666666667 0.66666666666667 1.008 > 28 49 H 0.16666666666667 0.50000000000000 0.66666666666667 1.008 > 28 50 H 0.50000000000000 0.50000000000000 0.66666666666667 1.008 > 28 51 H 0.83333333333333 0.50000000000000 0.66666666666667 1.008 > 28 52 H 0.16666666666667 0.83333333333333 0.66666666666667 1.008 > 28 53 H 0.50000000000000 0.83333333333333 0.66666666666667 1.008 > 28 54 H 0.83333333333333 0.83333333333333 0.66666666666667 1.008 > 28 55 H 0.00000000000000 0.16666666666667 0.16666666666667 1.008 > 55 56 H 0.33333333333333 0.16666666666667 0.16666666666667 1.008 > 55 57 H 0.66666666666667 0.16666666666667 0.16666666666667 1.008 > 55 58 H 0.00000000000000 0.50000000000000 0.16666666666667 1.008 > 55 59 H 0.33333333333333 0.50000000000000 0.16666666666667 1.008 > 55 60 H 0.66666666666667 0.50000000000000 0.16666666666667 1.008 > 55 61 H 0.00000000000000 0.83333333333333 0.16666666666667 1.008 > 55 62 H 0.33333333333333 0.83333333333333 0.16666666666667 1.008 > 55 63 H 0.66666666666667 0.83333333333333 0.16666666666667 1.008 > 55 64 H 0.00000000000000 0.16666666666667 0.50000000000000 1.008 > 55 65 H 0.33333333333333 0.16666666666667 0.50000000000000 1.008 > 55 66 H 0.66666666666667 0.16666666666667 0.50000000000000 1.008 > 55 67 H 0.00000000000000 0.50000000000000 0.50000000000000 1.008 > 55 68 H 0.33333333333333 0.50000000000000 0.50000000000000 1.008 > 55 69 H 0.66666666666667 0.50000000000000 0.50000000000000 1.008 > 55 70 H 0.00000000000000 0.83333333333333 0.50000000000000 1.008 > 55 71 H 0.33333333333333 0.83333333333333 0.50000000000000 1.008 > 55 72 H 0.66666666666667 0.83333333333333 0.50000000000000 1.008 > 55 73 H 0.00000000000000 0.16666666666667 0.83333333333333 1.008 > 55 74 H 0.33333333333333 0.16666666666667 0.83333333333333 1.008 > 55 75 H 0.66666666666667 0.16666666666667 0.83333333333333 1.008 > 55 76 H 0.00000000000000 0.50000000000000 0.83333333333333 1.008 > 55 77 H 0.33333333333333 0.50000000000000 0.83333333333333 1.008 > 55 78 H 0.66666666666667 0.50000000000000 0.83333333333333 1.008 > 55 79 H 0.00000000000000 0.83333333333333 0.83333333333333 1.008 > 55 80 H 0.33333333333333 0.83333333333333 0.83333333333333 1.008 > 55 81 H 0.66666666666667 0.83333333333333 0.83333333333333 1.008 > 55 82 Mg 0.16666666666667 0.16666666666667 0.16666666666667 24.305 > 82 83 Mg 0.50000000000000 0.16666666666667 0.16666666666667 24.305 > 82 84 Mg 0.83333333333333 0.16666666666667 0.16666666666667 24.305 > 82 85 Mg 0.16666666666667 0.50000000000000 0.16666666666667 24.305 > 82 86 Mg 0.50000000000000 0.50000000000000 0.16666666666667 24.305 > 82 87 Mg 0.83333333333333 0.50000000000000 0.16666666666667 24.305 > 82 88 Mg 0.16666666666667 0.83333333333333 0.16666666666667 24.305 > 82 89 Mg 0.50000000000000 0.83333333333333 0.16666666666667 24.305 > 82 90 Mg 0.83333333333333 0.83333333333333 0.16666666666667 24.305 > 82 91 Mg 0.16666666666667 0.16666666666667 0.50000000000000 24.305 > 82 92 Mg 0.50000000000000 0.16666666666667 0.50000000000000 24.305 > 82 93 Mg 0.83333333333333 0.16666666666667 0.50000000000000 24.305 > 82 94 Mg 0.16666666666667 0.50000000000000 0.50000000000000 24.305 > 82 95 Mg 0.50000000000000 0.50000000000000 0.50000000000000 24.305 > 82 96 Mg 0.83333333333333 0.50000000000000 0.50000000000000 24.305 > 82 97 Mg 0.16666666666667 0.83333333333333 0.50000000000000 24.305 > 82 98 Mg 0.50000000000000 0.83333333333333 0.50000000000000 24.305 > 82 99 Mg 0.83333333333333 0.83333333333333 0.50000000000000 24.305 > 82 100 Mg 0.16666666666667 0.16666666666667 0.83333333333333 24.305 > 82 101 Mg 0.50000000000000 0.16666666666667 0.83333333333333 24.305 > 82 102 Mg 0.83333333333333 0.16666666666667 0.83333333333333 24.305 > 82 103 Mg 0.16666666666667 0.50000000000000 0.83333333333333 24.305 > 82 104 Mg 0.50000000000000 0.50000000000000 0.83333333333333 24.305 > 82 105 Mg 0.83333333333333 0.50000000000000 0.83333333333333 24.305 > 82 106 Mg 0.16666666666667 0.83333333333333 0.83333333333333 24.305 > 82 107 Mg 0.50000000000000 0.83333333333333 0.83333333333333 24.305 > 82 108 Mg 0.83333333333333 0.83333333333333 0.83333333333333 24.305 > 82 109 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 > 109 110 K 0.33333333333333 0.00000000000000 0.00000000000000 39.098 > 109 111 K 0.66666666666667 0.00000000000000 0.00000000000000 39.098 > 109 112 K 0.00000000000000 0.33333333333333 0.00000000000000 39.098 > 109 113 K 0.33333333333333 0.33333333333333 0.00000000000000 39.098 > 109 114 K 0.66666666666667 0.33333333333333 0.00000000000000 39.098 > 109 115 K 0.00000000000000 0.66666666666667 0.00000000000000 39.098 > 109 116 K 0.33333333333333 0.66666666666667 0.00000000000000 39.098 > 109 117 K 0.66666666666667 0.66666666666667 0.00000000000000 39.098 > 109 118 K 0.00000000000000 0.00000000000000 0.33333333333333 39.098 > 109 119 K 0.33333333333333 0.00000000000000 0.33333333333333 39.098 > 109 120 K 0.66666666666667 0.00000000000000 0.33333333333333 39.098 > 109 121 K 0.00000000000000 0.33333333333333 0.33333333333333 39.098 > 109 122 K 0.33333333333333 0.33333333333333 0.33333333333333 39.098 > 109 123 K 0.66666666666667 0.33333333333333 0.33333333333333 39.098 > 109 124 K 0.00000000000000 0.66666666666667 0.33333333333333 39.098 > 109 125 K 0.33333333333333 0.66666666666667 0.33333333333333 39.098 > 109 126 K 0.66666666666667 0.66666666666667 0.33333333333333 39.098 > 109 127 K 0.00000000000000 0.00000000000000 0.66666666666667 39.098 > 109 128 K 0.33333333333333 0.00000000000000 0.66666666666667 39.098 > 109 129 K 0.66666666666667 0.00000000000000 0.66666666666667 39.098 > 109 130 K 0.00000000000000 0.33333333333333 0.66666666666667 39.098 > 109 131 K 0.33333333333333 0.33333333333333 0.66666666666667 39.098 > 109 132 K 0.66666666666667 0.33333333333333 0.66666666666667 39.098 > 109 133 K 0.00000000000000 0.66666666666667 0.66666666666667 39.098 > 109 134 K 0.33333333333333 0.66666666666667 0.66666666666667 39.098 > 109 135 K 0.66666666666667 0.66666666666667 0.66666666666667 39.098 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.7023609 0.0000000 0.0000000 0.0000000 3.7023609 0.0000000 0.0000000 0.0000000 3.7023609 -------------------------- Born effective charges -------------------------- 1 H -0.8075801 0.0000000 0.0000000 0.0000000 -1.2289746 0.0000000 0.0000000 0.0000000 -0.8075801 2 H -0.8075801 0.0000000 0.0000000 0.0000000 -0.8075801 0.0000000 0.0000000 0.0000000 -1.2289746 3 H -1.2289746 0.0000000 0.0000000 0.0000000 -0.8075801 0.0000000 0.0000000 0.0000000 -0.8075801 4 Mg 1.8670214 0.0000000 0.0000000 0.0000000 1.8670214 0.0000000 0.0000000 0.0000000 1.8670214 5 K 0.9771135 0.0000000 0.0000000 0.0000000 0.9771135 0.0000000 0.0000000 0.0000000 0.9771135 ---------------------------------------------------------------------------- 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 (yyy) -0.00000000 (yyy) -0.00000000 (yyy) 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-08 05:35:13]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 05:35:13]------------------------- 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.989339610000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.989339610000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.989339610000000 Atomic positions (fractional): 1 H 0.50000000000000 0.00000000000000 0.50000000000000 1.008 2 H 0.50000000000000 0.50000000000000 0.00000000000000 1.008 3 H 0.00000000000000 0.50000000000000 0.50000000000000 1.008 4 Mg 0.50000000000000 0.50000000000000 0.50000000000000 24.305 5 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.968018830000000 0.000000000000000 0.000000000000000 b 0.000000000000000 11.968018830000000 0.000000000000000 c 0.000000000000000 0.000000000000000 11.968018830000000 Atomic positions (fractional): 1 H 0.16666666666667 0.00000000000000 0.16666666666667 1.008 > 1 2 H 0.50000000000000 0.00000000000000 0.16666666666667 1.008 > 1 3 H 0.83333333333333 0.00000000000000 0.16666666666667 1.008 > 1 4 H 0.16666666666667 0.33333333333333 0.16666666666667 1.008 > 1 5 H 0.50000000000000 0.33333333333333 0.16666666666667 1.008 > 1 6 H 0.83333333333333 0.33333333333333 0.16666666666667 1.008 > 1 7 H 0.16666666666667 0.66666666666667 0.16666666666667 1.008 > 1 8 H 0.50000000000000 0.66666666666667 0.16666666666667 1.008 > 1 9 H 0.83333333333333 0.66666666666667 0.16666666666667 1.008 > 1 10 H 0.16666666666667 0.00000000000000 0.50000000000000 1.008 > 1 11 H 0.50000000000000 0.00000000000000 0.50000000000000 1.008 > 1 12 H 0.83333333333333 0.00000000000000 0.50000000000000 1.008 > 1 13 H 0.16666666666667 0.33333333333333 0.50000000000000 1.008 > 1 14 H 0.50000000000000 0.33333333333333 0.50000000000000 1.008 > 1 15 H 0.83333333333333 0.33333333333333 0.50000000000000 1.008 > 1 16 H 0.16666666666667 0.66666666666667 0.50000000000000 1.008 > 1 17 H 0.50000000000000 0.66666666666667 0.50000000000000 1.008 > 1 18 H 0.83333333333333 0.66666666666667 0.50000000000000 1.008 > 1 19 H 0.16666666666667 0.00000000000000 0.83333333333333 1.008 > 1 20 H 0.50000000000000 0.00000000000000 0.83333333333333 1.008 > 1 21 H 0.83333333333333 0.00000000000000 0.83333333333333 1.008 > 1 22 H 0.16666666666667 0.33333333333333 0.83333333333333 1.008 > 1 23 H 0.50000000000000 0.33333333333333 0.83333333333333 1.008 > 1 24 H 0.83333333333333 0.33333333333333 0.83333333333333 1.008 > 1 25 H 0.16666666666667 0.66666666666667 0.83333333333333 1.008 > 1 26 H 0.50000000000000 0.66666666666667 0.83333333333333 1.008 > 1 27 H 0.83333333333333 0.66666666666667 0.83333333333333 1.008 > 1 28 H 0.16666666666667 0.16666666666667 0.00000000000000 1.008 > 28 29 H 0.50000000000000 0.16666666666667 0.00000000000000 1.008 > 28 30 H 0.83333333333333 0.16666666666667 0.00000000000000 1.008 > 28 31 H 0.16666666666667 0.50000000000000 0.00000000000000 1.008 > 28 32 H 0.50000000000000 0.50000000000000 0.00000000000000 1.008 > 28 33 H 0.83333333333333 0.50000000000000 0.00000000000000 1.008 > 28 34 H 0.16666666666667 0.83333333333333 0.00000000000000 1.008 > 28 35 H 0.50000000000000 0.83333333333333 0.00000000000000 1.008 > 28 36 H 0.83333333333333 0.83333333333333 0.00000000000000 1.008 > 28 37 H 0.16666666666667 0.16666666666667 0.33333333333333 1.008 > 28 38 H 0.50000000000000 0.16666666666667 0.33333333333333 1.008 > 28 39 H 0.83333333333333 0.16666666666667 0.33333333333333 1.008 > 28 40 H 0.16666666666667 0.50000000000000 0.33333333333333 1.008 > 28 41 H 0.50000000000000 0.50000000000000 0.33333333333333 1.008 > 28 42 H 0.83333333333333 0.50000000000000 0.33333333333333 1.008 > 28 43 H 0.16666666666667 0.83333333333333 0.33333333333333 1.008 > 28 44 H 0.50000000000000 0.83333333333333 0.33333333333333 1.008 > 28 45 H 0.83333333333333 0.83333333333333 0.33333333333333 1.008 > 28 46 H 0.16666666666667 0.16666666666667 0.66666666666667 1.008 > 28 47 H 0.50000000000000 0.16666666666667 0.66666666666667 1.008 > 28 48 H 0.83333333333333 0.16666666666667 0.66666666666667 1.008 > 28 49 H 0.16666666666667 0.50000000000000 0.66666666666667 1.008 > 28 50 H 0.50000000000000 0.50000000000000 0.66666666666667 1.008 > 28 51 H 0.83333333333333 0.50000000000000 0.66666666666667 1.008 > 28 52 H 0.16666666666667 0.83333333333333 0.66666666666667 1.008 > 28 53 H 0.50000000000000 0.83333333333333 0.66666666666667 1.008 > 28 54 H 0.83333333333333 0.83333333333333 0.66666666666667 1.008 > 28 55 H 0.00000000000000 0.16666666666667 0.16666666666667 1.008 > 55 56 H 0.33333333333333 0.16666666666667 0.16666666666667 1.008 > 55 57 H 0.66666666666667 0.16666666666667 0.16666666666667 1.008 > 55 58 H 0.00000000000000 0.50000000000000 0.16666666666667 1.008 > 55 59 H 0.33333333333333 0.50000000000000 0.16666666666667 1.008 > 55 60 H 0.66666666666667 0.50000000000000 0.16666666666667 1.008 > 55 61 H 0.00000000000000 0.83333333333333 0.16666666666667 1.008 > 55 62 H 0.33333333333333 0.83333333333333 0.16666666666667 1.008 > 55 63 H 0.66666666666667 0.83333333333333 0.16666666666667 1.008 > 55 64 H 0.00000000000000 0.16666666666667 0.50000000000000 1.008 > 55 65 H 0.33333333333333 0.16666666666667 0.50000000000000 1.008 > 55 66 H 0.66666666666667 0.16666666666667 0.50000000000000 1.008 > 55 67 H 0.00000000000000 0.50000000000000 0.50000000000000 1.008 > 55 68 H 0.33333333333333 0.50000000000000 0.50000000000000 1.008 > 55 69 H 0.66666666666667 0.50000000000000 0.50000000000000 1.008 > 55 70 H 0.00000000000000 0.83333333333333 0.50000000000000 1.008 > 55 71 H 0.33333333333333 0.83333333333333 0.50000000000000 1.008 > 55 72 H 0.66666666666667 0.83333333333333 0.50000000000000 1.008 > 55 73 H 0.00000000000000 0.16666666666667 0.83333333333333 1.008 > 55 74 H 0.33333333333333 0.16666666666667 0.83333333333333 1.008 > 55 75 H 0.66666666666667 0.16666666666667 0.83333333333333 1.008 > 55 76 H 0.00000000000000 0.50000000000000 0.83333333333333 1.008 > 55 77 H 0.33333333333333 0.50000000000000 0.83333333333333 1.008 > 55 78 H 0.66666666666667 0.50000000000000 0.83333333333333 1.008 > 55 79 H 0.00000000000000 0.83333333333333 0.83333333333333 1.008 > 55 80 H 0.33333333333333 0.83333333333333 0.83333333333333 1.008 > 55 81 H 0.66666666666667 0.83333333333333 0.83333333333333 1.008 > 55 82 Mg 0.16666666666667 0.16666666666667 0.16666666666667 24.305 > 82 83 Mg 0.50000000000000 0.16666666666667 0.16666666666667 24.305 > 82 84 Mg 0.83333333333333 0.16666666666667 0.16666666666667 24.305 > 82 85 Mg 0.16666666666667 0.50000000000000 0.16666666666667 24.305 > 82 86 Mg 0.50000000000000 0.50000000000000 0.16666666666667 24.305 > 82 87 Mg 0.83333333333333 0.50000000000000 0.16666666666667 24.305 > 82 88 Mg 0.16666666666667 0.83333333333333 0.16666666666667 24.305 > 82 89 Mg 0.50000000000000 0.83333333333333 0.16666666666667 24.305 > 82 90 Mg 0.83333333333333 0.83333333333333 0.16666666666667 24.305 > 82 91 Mg 0.16666666666667 0.16666666666667 0.50000000000000 24.305 > 82 92 Mg 0.50000000000000 0.16666666666667 0.50000000000000 24.305 > 82 93 Mg 0.83333333333333 0.16666666666667 0.50000000000000 24.305 > 82 94 Mg 0.16666666666667 0.50000000000000 0.50000000000000 24.305 > 82 95 Mg 0.50000000000000 0.50000000000000 0.50000000000000 24.305 > 82 96 Mg 0.83333333333333 0.50000000000000 0.50000000000000 24.305 > 82 97 Mg 0.16666666666667 0.83333333333333 0.50000000000000 24.305 > 82 98 Mg 0.50000000000000 0.83333333333333 0.50000000000000 24.305 > 82 99 Mg 0.83333333333333 0.83333333333333 0.50000000000000 24.305 > 82 100 Mg 0.16666666666667 0.16666666666667 0.83333333333333 24.305 > 82 101 Mg 0.50000000000000 0.16666666666667 0.83333333333333 24.305 > 82 102 Mg 0.83333333333333 0.16666666666667 0.83333333333333 24.305 > 82 103 Mg 0.16666666666667 0.50000000000000 0.83333333333333 24.305 > 82 104 Mg 0.50000000000000 0.50000000000000 0.83333333333333 24.305 > 82 105 Mg 0.83333333333333 0.50000000000000 0.83333333333333 24.305 > 82 106 Mg 0.16666666666667 0.83333333333333 0.83333333333333 24.305 > 82 107 Mg 0.50000000000000 0.83333333333333 0.83333333333333 24.305 > 82 108 Mg 0.83333333333333 0.83333333333333 0.83333333333333 24.305 > 82 109 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 > 109 110 K 0.33333333333333 0.00000000000000 0.00000000000000 39.098 > 109 111 K 0.66666666666667 0.00000000000000 0.00000000000000 39.098 > 109 112 K 0.00000000000000 0.33333333333333 0.00000000000000 39.098 > 109 113 K 0.33333333333333 0.33333333333333 0.00000000000000 39.098 > 109 114 K 0.66666666666667 0.33333333333333 0.00000000000000 39.098 > 109 115 K 0.00000000000000 0.66666666666667 0.00000000000000 39.098 > 109 116 K 0.33333333333333 0.66666666666667 0.00000000000000 39.098 > 109 117 K 0.66666666666667 0.66666666666667 0.00000000000000 39.098 > 109 118 K 0.00000000000000 0.00000000000000 0.33333333333333 39.098 > 109 119 K 0.33333333333333 0.00000000000000 0.33333333333333 39.098 > 109 120 K 0.66666666666667 0.00000000000000 0.33333333333333 39.098 > 109 121 K 0.00000000000000 0.33333333333333 0.33333333333333 39.098 > 109 122 K 0.33333333333333 0.33333333333333 0.33333333333333 39.098 > 109 123 K 0.66666666666667 0.33333333333333 0.33333333333333 39.098 > 109 124 K 0.00000000000000 0.66666666666667 0.33333333333333 39.098 > 109 125 K 0.33333333333333 0.66666666666667 0.33333333333333 39.098 > 109 126 K 0.66666666666667 0.66666666666667 0.33333333333333 39.098 > 109 127 K 0.00000000000000 0.00000000000000 0.66666666666667 39.098 > 109 128 K 0.33333333333333 0.00000000000000 0.66666666666667 39.098 > 109 129 K 0.66666666666667 0.00000000000000 0.66666666666667 39.098 > 109 130 K 0.00000000000000 0.33333333333333 0.66666666666667 39.098 > 109 131 K 0.33333333333333 0.33333333333333 0.66666666666667 39.098 > 109 132 K 0.66666666666667 0.33333333333333 0.66666666666667 39.098 > 109 133 K 0.00000000000000 0.66666666666667 0.66666666666667 39.098 > 109 134 K 0.33333333333333 0.66666666666667 0.66666666666667 39.098 > 109 135 K 0.66666666666667 0.66666666666667 0.66666666666667 39.098 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 3.7023609 0.0000000 0.0000000 0.0000000 3.7023609 0.0000000 0.0000000 0.0000000 3.7023609 -------------------------- Born effective charges -------------------------- 1 H -0.8075801 0.0000000 0.0000000 0.0000000 -1.2289746 0.0000000 0.0000000 0.0000000 -0.8075801 2 H -0.8075801 0.0000000 0.0000000 0.0000000 -0.8075801 0.0000000 0.0000000 0.0000000 -1.2289746 3 H -1.2289746 0.0000000 0.0000000 0.0000000 -0.8075801 0.0000000 0.0000000 0.0000000 -0.8075801 4 Mg 1.8670214 0.0000000 0.0000000 0.0000000 1.8670214 0.0000000 0.0000000 0.0000000 1.8670214 5 K 0.9771135 0.0000000 0.0000000 0.0000000 0.9771135 0.0000000 0.0000000 0.0000000 0.9771135 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (yyy) -0.00000000 (yyy) -0.00000000 (yyy) Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 13 13 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.21 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: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 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 5.684 ( 0.000 0.000 0.000) 0.000 5.684 ( 0.000 0.000 0.000) 0.000 5.684 ( 0.000 0.000 0.000) 0.000 17.712 ( 0.000 0.000 0.000) 0.000 17.712 ( 0.000 0.000 0.000) 0.000 17.712 ( 0.000 0.000 0.000) 0.000 19.211 ( 0.000 0.000 0.000) 0.000 19.211 ( 0.000 0.000 0.000) 0.000 19.211 ( 0.000 0.000 0.000) 0.000 35.814 ( 0.000 0.000 0.000) 0.000 35.814 ( 0.000 0.000 0.000) 0.000 35.814 ( 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: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.902 ( 45.799 0.000 0.000) 45.799 0.902 ( 45.799 0.000 0.000) 45.799 1.225 ( 61.896 0.000 0.000) 61.896 5.670 ( -1.413 0.000 0.000) 1.413 5.670 ( -1.413 0.000 0.000) 1.413 5.808 ( 8.582 0.000 0.000) 8.582 17.767 ( 5.718 0.000 0.000) 5.718 17.767 ( 5.718 0.000 0.000) 5.718 19.238 ( 2.727 0.000 0.000) 2.727 19.238 ( 2.727 0.000 0.000) 2.727 19.273 ( 6.307 0.000 0.000) 6.307 25.432 ( 0.145 0.000 0.000) 0.145 35.797 ( -1.697 0.000 0.000) 1.697 35.797 ( -1.697 0.000 0.000) 1.697 38.795 ( -1.362 0.000 0.000) 1.362 ======================= Grid point 2 (3/84) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.748 ( 41.519 0.000 0.000) 41.519 1.748 ( 41.519 0.000 0.000) 41.519 2.357 ( 54.736 0.000 0.000) 54.736 5.630 ( -2.681 0.000 0.000) 2.681 5.630 ( -2.681 0.000 0.000) 2.681 6.046 ( 15.762 0.000 0.000) 15.762 17.926 ( 10.464 0.000 0.000) 10.464 17.926 ( 10.464 0.000 0.000) 10.464 19.311 ( 4.576 0.000 0.000) 4.576 19.311 ( 4.576 0.000 0.000) 4.576 19.444 ( 11.056 0.000 0.000) 11.056 25.436 ( 0.243 0.000 0.000) 0.243 35.751 ( -3.045 0.000 0.000) 3.045 35.751 ( -3.045 0.000 0.000) 3.045 38.758 ( -2.459 0.000 0.000) 2.459 ======================= Grid point 3 (4/84) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.489 ( 35.006 0.000 0.000) 35.006 2.489 ( 35.006 0.000 0.000) 35.006 3.311 ( 43.725 0.000 0.000) 43.725 5.569 ( -3.647 0.000 0.000) 3.647 5.569 ( -3.647 0.000 0.000) 3.647 6.398 ( 20.136 0.000 0.000) 20.136 18.159 ( 13.340 0.000 0.000) 13.340 18.159 ( 13.340 0.000 0.000) 13.340 19.405 ( 4.945 0.000 0.000) 4.945 19.405 ( 4.945 0.000 0.000) 4.945 19.682 ( 13.150 0.000 0.000) 13.150 25.441 ( 0.269 0.000 0.000) 0.269 35.684 ( -3.744 0.000 0.000) 3.744 35.684 ( -3.744 0.000 0.000) 3.744 38.704 ( -3.050 0.000 0.000) 3.050 ======================= Grid point 4 (5/84) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.088 ( 26.842 0.000 0.000) 26.842 3.088 ( 26.842 0.000 0.000) 26.842 4.030 ( 30.783 0.000 0.000) 30.783 5.493 ( -4.079 0.000 0.000) 4.079 5.493 ( -4.079 0.000 0.000) 4.079 6.795 ( 20.130 0.000 0.000) 20.130 18.423 ( 13.592 0.000 0.000) 13.592 18.423 ( 13.592 0.000 0.000) 13.592 19.491 ( 3.773 0.000 0.000) 3.773 19.491 ( 3.773 0.000 0.000) 3.773 19.932 ( 12.207 0.000 0.000) 12.207 25.446 ( 0.230 0.000 0.000) 0.230 35.612 ( -3.600 0.000 0.000) 3.600 35.612 ( -3.600 0.000 0.000) 3.600 38.645 ( -2.959 0.000 0.000) 2.959 ======================= Grid point 5 (6/84) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.515 ( 17.246 0.000 0.000) 17.246 3.515 ( 17.246 0.000 0.000) 17.246 4.499 ( 17.931 0.000 0.000) 17.931 5.417 ( -3.536 0.000 0.000) 3.536 5.417 ( -3.536 0.000 0.000) 3.536 7.140 ( 14.877 0.000 0.000) 14.877 18.662 ( 10.624 0.000 0.000) 10.624 18.662 ( 10.624 0.000 0.000) 10.624 19.545 ( 1.774 0.000 0.000) 1.774 19.545 ( 1.774 0.000 0.000) 1.774 20.136 ( 8.556 0.000 0.000) 8.556 25.450 ( 0.150 0.000 0.000) 0.150 35.551 ( -2.598 0.000 0.000) 2.598 35.551 ( -2.598 0.000 0.000) 2.598 38.594 ( -2.151 0.000 0.000) 2.151 ======================= Grid point 6 (7/84) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.742 ( 6.055 0.000 0.000) 6.055 3.742 ( 6.055 0.000 0.000) 6.055 4.727 ( 5.859 0.000 0.000) 5.859 5.366 ( -1.487 0.000 0.000) 1.487 5.366 ( -1.487 0.000 0.000) 1.487 7.341 ( 5.469 0.000 0.000) 5.469 18.810 ( 4.168 0.000 0.000) 4.168 18.810 ( 4.168 0.000 0.000) 4.168 19.563 ( 0.323 0.000 0.000) 0.323 19.563 ( 0.323 0.000 0.000) 0.323 20.250 ( 3.068 0.000 0.000) 3.068 25.452 ( 0.051 0.000 0.000) 0.051 35.516 ( -0.947 0.000 0.000) 0.947 35.516 ( -0.947 0.000 0.000) 0.947 38.565 ( -0.787 0.000 0.000) 0.787 ======================= Grid point 14 (8/84) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.005 ( 25.182 25.182 0.000) 35.613 1.265 ( 31.493 31.493 0.000) 44.538 1.873 ( 45.981 45.981 0.000) 65.027 5.658 ( -1.206 -1.206 0.000) 1.705 5.723 ( 0.209 0.209 0.000) 0.296 5.825 ( 7.126 7.126 0.000) 10.077 17.772 ( 2.504 2.504 0.000) 3.541 17.825 ( 5.908 5.908 0.000) 8.356 19.103 ( -5.104 -5.104 0.000) 7.219 19.265 ( 2.732 2.732 0.000) 3.864 19.494 ( 13.865 13.865 0.000) 19.608 25.432 ( -0.027 -0.027 0.000) 0.038 35.616 ( -10.352 -10.352 0.000) 14.640 35.781 ( -1.688 -1.688 0.000) 2.388 38.946 ( 6.817 6.817 0.000) 9.641 ======================= Grid point 15 (9/84) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.674 ( 38.794 1.658 0.000) 38.830 1.938 ( 35.897 18.223 0.000) 40.257 2.771 ( 45.076 32.898 0.000) 55.804 5.624 ( -2.260 -0.577 0.000) 2.332 5.690 ( -2.135 5.501 0.000) 5.901 6.067 ( 16.285 3.030 0.000) 16.565 17.878 ( 8.091 -4.910 0.000) 9.464 17.989 ( 10.772 6.406 0.000) 12.533 19.045 ( -1.205 -11.107 0.000) 11.172 19.338 ( 4.606 2.766 0.000) 5.373 19.788 ( 16.366 20.373 0.000) 26.132 25.431 ( -0.049 -0.486 0.000) 0.489 35.455 ( -7.273 -23.917 0.000) 24.999 35.734 ( -3.027 -1.662 0.000) 3.453 39.003 ( 0.269 18.511 0.000) 18.513 ======================= Grid point 16 (10/84) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.402 ( 35.553 -4.524 0.000) 35.840 2.598 ( 31.817 10.700 0.000) 33.568 3.559 ( 35.666 20.564 0.000) 41.170 5.573 ( -3.009 0.479 0.000) 3.046 5.654 ( -1.470 8.993 0.000) 9.113 6.428 ( 20.623 3.263 0.000) 20.879 18.071 ( 11.638 -12.522 0.000) 17.096 18.228 ( 13.653 7.023 0.000) 15.353 19.044 ( 0.717 -14.280 0.000) 14.298 19.433 ( 5.024 2.876 0.000) 5.789 20.109 ( 16.512 25.153 0.000) 30.088 25.430 ( -0.100 -1.134 0.000) 1.139 35.318 ( -7.039 -32.398 0.000) 33.153 35.668 ( -3.718 -1.617 0.000) 4.054 38.985 ( -1.818 23.525 0.000) 23.595 ======================= Grid point 17 (11/84) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.020 ( 28.120 -5.053 0.000) 28.571 3.145 ( 24.610 5.674 0.000) 25.255 4.122 ( 22.490 7.190 0.000) 23.611 5.512 ( -3.258 1.929 0.000) 3.786 5.640 ( 0.164 15.296 0.000) 15.296 6.835 ( 20.670 4.115 0.000) 21.076 18.313 ( 13.083 -20.583 0.000) 24.389 18.498 ( 13.814 7.548 0.000) 15.742 19.059 ( 0.432 -14.246 0.000) 14.252 19.522 ( 3.903 3.089 0.000) 4.977 20.408 ( 14.056 28.276 0.000) 31.577 25.427 ( -0.137 -1.851 0.000) 1.856 35.186 ( -6.423 -38.814 0.000) 39.342 35.596 ( -3.571 -1.560 0.000) 3.896 38.943 ( -2.263 25.783 0.000) 25.882 ======================= Grid point 18 (12/84) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.473 ( 18.558 -3.674 0.000) 18.918 3.536 ( 15.678 2.030 0.000) 15.809 4.430 ( 10.016 -5.765 0.000) 11.556 5.453 ( -2.709 3.567 0.000) 4.479 5.660 ( 1.707 22.913 0.000) 22.976 7.190 ( 15.297 5.035 0.000) 16.105 18.559 ( 12.001 -30.467 0.000) 32.745 18.740 ( 10.732 7.861 0.000) 13.303 19.049 ( -1.658 -9.479 0.000) 9.623 19.578 ( 1.905 3.360 0.000) 3.863 20.637 ( 9.400 30.011 0.000) 31.449 25.425 ( -0.123 -2.469 0.000) 2.472 35.078 ( -4.597 -43.652 0.000) 43.893 35.536 ( -2.575 -1.505 0.000) 2.982 38.903 ( -1.768 26.867 0.000) 26.925 ======================= Grid point 19 (13/84) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.721 ( 6.711 -2.027 0.000) 7.010 3.741 ( 5.434 -0.114 0.000) 5.436 4.540 ( 2.263 -13.704 0.000) 13.890 5.414 ( -1.099 4.792 0.000) 4.917 5.691 ( 1.047 27.519 0.000) 27.538 7.397 ( 5.623 5.588 0.000) 7.927 18.752 ( 6.904 -48.654 0.000) 49.141 18.889 ( 4.198 7.985 0.000) 9.021 18.997 ( -3.042 6.298 0.000) 6.994 19.598 ( 0.377 3.554 0.000) 3.574 20.761 ( 3.300 30.749 0.000) 30.925 25.423 ( -0.049 -2.827 0.000) 2.827 35.016 ( -1.677 -46.343 0.000) 46.373 35.501 ( -0.938 -1.470 0.000) 1.744 38.879 ( -0.663 27.315 0.000) 27.323 ======================= Grid point 28 (14/84) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.917 ( 21.848 21.848 0.000) 30.898 2.377 ( 25.469 25.469 0.000) 36.019 3.423 ( 32.142 32.142 0.000) 45.456 5.609 ( -1.002 -1.002 0.000) 1.417 5.804 ( 5.663 5.663 0.000) 8.008 6.209 ( 12.425 12.425 0.000) 17.572 17.756 ( -6.842 -6.842 0.000) 9.676 18.165 ( 11.574 11.574 0.000) 16.368 18.951 ( 0.555 0.555 0.000) 0.785 19.412 ( 4.720 4.720 0.000) 6.675 20.176 ( 19.835 19.835 0.000) 28.052 25.418 ( -0.828 -0.828 0.000) 1.171 35.005 ( -21.444 -21.444 0.000) 30.327 35.689 ( -2.972 -2.972 0.000) 4.203 39.294 ( 10.524 10.524 0.000) 14.883 ======================= Grid point 29 (15/84) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.424 ( 28.686 7.547 0.000) 29.663 2.866 ( 24.305 16.053 0.000) 29.128 3.946 ( 21.134 17.160 0.000) 27.223 5.588 ( -1.164 1.037 0.000) 1.559 5.918 ( 6.721 18.169 0.000) 19.372 6.534 ( 20.321 8.011 0.000) 21.843 17.641 ( -4.662 -28.282 0.000) 28.664 18.420 ( 14.461 12.560 0.000) 19.154 19.026 ( 6.252 9.244 0.000) 11.160 19.511 ( 5.274 4.976 0.000) 7.250 20.547 ( 18.370 20.281 0.000) 27.364 25.400 ( -1.034 -1.867 0.000) 2.134 34.629 ( -17.884 -36.334 0.000) 40.497 35.624 ( -3.636 -2.878 0.000) 4.637 39.422 ( 3.506 19.370 0.000) 19.684 ======================= Grid point 30 (16/84) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.964 ( 26.262 0.302 0.000) 26.264 3.287 ( 18.944 8.472 0.000) 20.752 4.229 ( 8.701 3.018 0.000) 9.210 5.566 ( -1.030 3.638 0.000) 3.781 6.065 ( 8.178 27.689 0.000) 28.871 6.951 ( 21.734 7.811 0.000) 23.095 17.581 ( -1.708 -46.695 0.000) 46.727 18.703 ( 14.379 13.383 0.000) 19.643 19.164 ( 7.519 17.369 0.000) 18.926 19.605 ( 4.283 5.403 0.000) 6.895 20.871 ( 14.916 19.613 0.000) 24.641 25.380 ( -1.027 -2.933 0.000) 3.108 34.313 ( -14.743 -48.469 0.000) 50.662 35.554 ( -3.479 -2.757 0.000) 4.439 39.453 ( 0.242 24.284 0.000) 24.285 ======================= Grid point 31 (17/84) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.408 ( 19.175 -2.309 0.000) 19.314 3.585 ( 11.768 2.824 0.000) 12.102 4.311 ( 0.705 -6.343 0.000) 6.382 5.549 ( -0.674 6.245 0.000) 6.281 6.213 ( 6.590 32.657 0.000) 33.315 7.327 ( 16.284 8.960 0.000) 18.587 17.567 ( 0.033 -59.400 0.000) 59.400 18.953 ( 10.997 13.855 0.000) 17.688 19.297 ( 5.849 23.095 0.000) 23.825 19.669 ( 2.279 5.914 0.000) 6.338 21.111 ( 9.702 18.796 0.000) 21.152 25.362 ( -0.759 -3.794 0.000) 3.869 34.070 ( -10.132 -57.353 0.000) 58.241 35.495 ( -2.501 -2.640 0.000) 3.636 39.446 ( -0.708 26.559 0.000) 26.569 ======================= Grid point 32 (18/84) ======================= q-point: ( 0.46 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.677 ( 7.817 -2.124 0.000) 8.101 3.737 ( 3.964 -0.322 0.000) 3.977 4.293 ( -1.392 -11.405 0.000) 11.489 5.541 ( -0.228 7.944 0.000) 7.948 6.303 ( 2.475 34.521 0.000) 34.609 7.547 ( 5.985 9.727 0.000) 11.421 17.572 ( 0.235 -65.809 0.000) 65.809 19.106 ( 4.271 14.035 0.000) 14.670 19.377 ( 2.195 26.061 0.000) 26.153 19.695 ( 0.528 6.266 0.000) 6.288 21.238 ( 3.363 18.294 0.000) 18.601 25.352 ( -0.279 -4.270 0.000) 4.279 33.934 ( -3.652 -62.185 0.000) 62.292 35.462 ( -0.910 -2.567 0.000) 2.723 39.434 ( -0.378 27.397 0.000) 27.399 ======================= Grid point 42 (19/84) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.675 ( 17.289 17.289 0.000) 24.450 3.180 ( 15.721 15.721 0.000) 22.233 4.179 ( 7.043 7.043 0.000) 9.961 5.613 ( 1.607 1.607 0.000) 2.273 6.313 ( 20.704 20.704 0.000) 29.280 6.755 ( 15.509 15.509 0.000) 21.934 17.107 ( -25.535 -25.535 0.000) 36.111 18.695 ( 15.435 15.435 0.000) 21.829 19.279 ( 15.229 15.229 0.000) 21.537 19.616 ( 5.718 5.718 0.000) 8.087 20.898 ( 16.261 16.261 0.000) 22.997 25.361 ( -2.156 -2.156 0.000) 3.049 33.962 ( -32.400 -32.400 0.000) 45.820 35.561 ( -3.497 -3.497 0.000) 4.945 39.704 ( 10.140 10.140 0.000) 14.340 ======================= Grid point 43 (20/84) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.043 ( 19.813 7.683 0.000) 21.251 3.451 ( 12.039 8.047 0.000) 14.481 4.237 ( -0.344 -1.845 0.000) 1.876 5.648 ( 1.909 4.701 0.000) 5.074 6.650 ( 13.975 31.316 0.000) 34.293 7.131 ( 21.981 10.740 0.000) 24.465 16.685 ( -17.909 -44.291 0.000) 47.775 18.994 ( 15.028 16.211 0.000) 22.105 19.571 ( 14.323 23.097 0.000) 27.178 19.721 ( 4.889 6.354 0.000) 8.017 21.175 ( 12.397 12.327 0.000) 17.483 25.320 ( -1.991 -3.158 0.000) 3.733 33.390 ( -26.598 -46.182 0.000) 53.294 35.494 ( -3.321 -3.317 0.000) 4.694 39.833 ( 3.895 14.976 0.000) 15.474 ======================= Grid point 44 (21/84) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.402 ( 16.609 1.842 0.000) 16.710 3.637 ( 7.208 2.411 0.000) 7.600 4.189 ( -4.153 -6.311 0.000) 7.555 5.683 ( 1.601 7.301 0.000) 7.474 6.866 ( 8.608 33.874 0.000) 34.950 7.525 ( 17.311 11.232 0.000) 20.635 16.416 ( -10.105 -58.323 0.000) 59.192 19.253 ( 11.273 16.605 0.000) 20.070 19.797 ( 2.843 7.055 0.000) 7.607 19.810 ( 10.033 28.801 0.000) 30.498 21.373 ( 7.961 8.895 0.000) 11.938 25.287 ( -1.393 -3.873 0.000) 4.116 32.954 ( -18.106 -56.892 0.000) 59.704 35.438 ( -2.373 -3.141 0.000) 3.937 39.875 ( 0.913 17.605 0.000) 17.628 ======================= Grid point 45 (22/84) ======================= q-point: ( 0.46 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.652 ( 8.016 -0.313 0.000) 8.023 3.729 ( 2.361 -0.533 0.000) 2.421 4.102 ( -3.745 -8.583 0.000) 9.365 5.705 ( 0.629 8.746 0.000) 8.768 6.980 ( 3.048 34.667 0.000) 34.800 7.759 ( 6.357 11.870 0.000) 13.465 16.290 ( -3.205 -65.812 0.000) 65.890 19.408 ( 4.349 16.737 0.000) 17.293 19.830 ( 0.743 7.511 0.000) 7.548 19.944 ( 3.602 31.713 0.000) 31.917 21.477 ( 2.777 7.119 0.000) 7.641 25.268 ( -0.498 -4.238 0.000) 4.267 32.712 ( -6.505 -62.846 0.000) 63.181 35.406 ( -0.860 -3.032 0.000) 3.151 39.882 ( 0.060 18.661 0.000) 18.661 ======================= Grid point 56 (23/84) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.239 ( 11.921 11.921 0.000) 16.859 3.586 ( 5.872 5.872 0.000) 8.304 4.172 ( -4.594 -4.594 0.000) 6.497 5.740 ( 4.649 4.649 0.000) 6.575 7.200 ( 22.670 22.670 0.000) 32.060 7.375 ( 16.226 16.226 0.000) 22.947 15.917 ( -34.625 -34.625 0.000) 48.967 19.305 ( 15.475 15.475 0.000) 21.885 19.840 ( 5.673 5.673 0.000) 8.023 20.009 ( 21.375 21.375 0.000) 30.229 21.363 ( 7.592 7.592 0.000) 10.737 25.262 ( -2.732 -2.732 0.000) 3.864 32.563 ( -38.832 -38.832 0.000) 54.916 35.431 ( -3.116 -3.116 0.000) 4.406 40.041 ( 7.112 7.112 0.000) 10.057 ======================= Grid point 57 (24/84) ======================= q-point: ( 0.38 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.475 ( 11.828 5.413 0.000) 13.008 3.675 ( 3.361 1.551 0.000) 3.702 4.069 ( -5.836 -6.010 0.000) 8.378 5.819 ( 3.376 6.566 0.000) 7.383 7.480 ( 9.480 28.906 0.000) 30.421 7.748 ( 18.060 11.576 0.000) 21.452 15.370 ( -21.764 -48.996 0.000) 53.612 19.569 ( 11.388 15.587 0.000) 19.304 19.930 ( 3.543 6.503 0.000) 7.405 20.364 ( 14.745 27.540 0.000) 31.239 21.479 ( 4.632 2.896 0.000) 5.463 25.217 ( -1.835 -3.205 0.000) 3.693 31.924 ( -26.641 -48.827 0.000) 55.623 35.379 ( -2.205 -2.915 0.000) 3.655 40.129 ( 2.523 9.133 0.000) 9.475 ======================= Grid point 58 (25/84) ======================= q-point: ( 0.46 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.667 ( 6.962 1.831 0.000) 7.198 3.718 ( 1.073 -0.603 0.000) 1.231 3.960 ( -4.678 -6.193 0.000) 7.761 5.865 ( 1.228 7.499 0.000) 7.598 7.605 ( 3.387 29.364 0.000) 29.558 7.992 ( 6.605 12.002 0.000) 13.699 15.088 ( -7.407 -57.448 0.000) 57.924 19.726 ( 4.393 15.564 0.000) 16.172 19.973 ( 0.981 7.026 0.000) 7.094 20.558 ( 5.156 30.894 0.000) 31.321 21.542 ( 1.744 0.449 0.000) 1.800 25.193 ( -0.644 -3.438 0.000) 3.498 31.568 ( -9.609 -54.532 0.000) 55.372 35.349 ( -0.795 -2.788 0.000) 2.899 40.156 ( 0.522 10.025 0.000) 10.039 ======================= Grid point 70 (26/84) ======================= q-point: ( 0.38 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.595 ( 6.640 6.640 0.000) 9.391 3.698 ( 0.784 0.784 0.000) 1.109 3.959 ( -5.446 -5.446 0.000) 7.702 5.927 ( 4.453 4.453 0.000) 6.297 7.945 ( 15.234 15.234 0.000) 21.544 7.960 ( 13.337 13.337 0.000) 18.861 14.570 ( -33.041 -33.041 0.000) 46.727 19.832 ( 11.215 11.215 0.000) 15.860 20.037 ( 4.331 4.331 0.000) 6.125 20.830 ( 19.776 19.776 0.000) 27.967 21.510 ( 0.798 0.798 0.000) 1.129 25.166 ( -2.117 -2.117 0.000) 2.994 31.115 ( -33.842 -33.842 0.000) 47.859 35.330 ( -2.040 -2.040 0.000) 2.885 40.247 ( 3.645 3.645 0.000) 5.156 ======================= Grid point 71 (27/84) ======================= q-point: ( 0.46 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.707 ( 0.230 -0.473 0.000) 0.526 3.714 ( 4.985 2.704 0.000) 5.671 3.862 ( -4.265 -3.992 0.000) 5.842 5.986 ( 1.564 4.933 0.000) 5.175 8.083 ( 3.684 19.599 0.000) 19.942 8.206 ( 6.747 9.655 0.000) 11.779 14.131 ( -11.771 -40.309 0.000) 41.992 19.986 ( 4.354 10.929 0.000) 11.764 20.091 ( 1.231 4.940 0.000) 5.091 21.091 ( 6.799 23.193 0.000) 24.169 21.524 ( 0.640 -1.702 0.000) 1.818 25.138 ( -0.743 -2.272 0.000) 2.390 30.662 ( -12.280 -38.060 0.000) 39.992 35.303 ( -0.731 -1.937 0.000) 2.070 40.287 ( 0.864 4.177 0.000) 4.266 ======================= Grid point 84 (28/84) ======================= q-point: ( 0.46 0.46 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.701 ( -0.178 -0.178 0.000) 0.252 3.761 ( 2.089 2.089 0.000) 2.954 3.803 ( -2.216 -2.216 0.000) 3.134 6.051 ( 1.708 1.708 0.000) 2.416 8.340 ( 5.526 5.526 0.000) 7.815 8.344 ( 5.245 5.245 0.000) 7.417 13.588 ( -14.761 -14.761 0.000) 20.875 20.133 ( 4.050 4.050 0.000) 5.728 20.156 ( 1.645 1.645 0.000) 2.326 21.407 ( 8.504 8.504 0.000) 12.027 21.497 ( -0.614 -0.614 0.000) 0.868 25.108 ( -0.804 -0.804 0.000) 1.137 30.150 ( -13.881 -13.881 0.000) 19.631 35.277 ( -0.690 -0.690 0.000) 0.976 40.334 ( 1.033 1.033 0.000) 1.461 ======================= Grid point 183 (29/84) ======================= q-point: ( 0.08 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.331 ( 21.750 21.750 21.750) 37.671 1.331 ( 21.750 21.750 21.750) 37.671 2.328 ( 36.893 36.893 36.893) 63.901 5.649 ( -1.804 -1.804 -1.804) 3.124 5.811 ( 4.263 4.263 4.263) 7.384 5.811 ( 4.263 4.263 4.263) 7.384 17.807 ( 2.261 2.261 2.261) 3.916 17.807 ( 2.261 2.261 2.261) 3.916 18.915 ( -10.227 -10.227 -10.227) 17.713 19.542 ( 11.543 11.543 11.543) 19.993 19.542 ( 11.543 11.543 11.543) 19.993 25.431 ( -0.107 -0.107 -0.107) 0.185 35.598 ( -7.546 -7.546 -7.546) 13.070 35.598 ( -7.546 -7.546 -7.546) 13.070 39.093 ( 9.305 9.305 9.305) 16.117 ======================= Grid point 184 (30/84) ======================= q-point: ( 0.15 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.811 ( 31.694 10.419 10.419) 34.951 1.971 ( 34.794 10.573 10.573) 37.871 3.078 ( 38.174 26.180 26.180) 53.179 5.621 ( -0.619 0.331 0.331) 0.775 5.771 ( -2.669 7.123 7.123) 10.421 6.076 ( 17.022 1.658 1.658) 17.183 17.905 ( 6.512 -1.182 -1.182) 6.724 17.914 ( 9.922 -2.088 -2.088) 10.352 18.764 ( -5.488 -16.976 -16.976) 24.627 19.613 ( 3.351 17.151 17.151) 24.486 19.940 ( 21.404 14.281 14.281) 29.428 25.427 ( -0.326 -0.463 -0.463) 0.731 35.366 ( -10.814 -13.208 -13.208) 21.583 35.548 ( -3.299 -10.579 -10.579) 15.320 39.207 ( 2.741 16.230 16.230) 23.116 ======================= Grid point 185 (31/84) ======================= q-point: ( 0.23 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.445 ( 32.274 2.315 2.315) 32.439 2.619 ( 31.654 6.417 6.417) 32.929 3.734 ( 28.709 14.859 14.859) 35.577 5.634 ( 2.125 5.002 5.002) 7.386 5.710 ( -3.487 7.232 7.232) 10.805 6.450 ( 21.269 2.490 2.490) 21.558 18.051 ( 8.418 -6.219 -6.219) 12.174 18.169 ( 15.944 -3.928 -3.928) 16.884 18.697 ( -1.715 -21.260 -21.260) 30.116 19.666 ( 1.946 17.124 17.124) 24.295 20.349 ( 20.487 19.594 19.594) 34.461 25.419 ( -0.493 -1.108 -1.108) 1.643 35.169 ( -9.813 -20.383 -20.383) 30.451 35.475 ( -4.109 -10.861 -10.861) 15.900 39.223 ( -0.563 20.519 20.519) 29.024 ======================= Grid point 186 (32/84) ======================= q-point: ( 0.31 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.015 ( 26.234 -1.160 -1.160) 26.285 3.169 ( 24.936 4.057 4.057) 25.587 4.163 ( 15.788 3.014 3.014) 16.354 5.640 ( -3.699 7.526 7.526) 11.268 5.703 ( 4.907 11.936 11.936) 17.579 6.869 ( 21.282 3.590 3.590) 21.879 18.219 ( 8.779 -13.532 -13.532) 21.055 18.508 ( 18.718 -3.329 -3.329) 19.301 18.685 ( 0.091 -21.849 -21.849) 30.900 19.678 ( -0.970 16.038 16.038) 22.702 20.713 ( 16.881 23.408 23.408) 37.160 25.409 ( -0.542 -1.877 -1.877) 2.709 34.990 ( -8.626 -25.988 -25.988) 37.751 35.395 ( -4.007 -11.219 -11.219) 16.364 39.199 ( -1.628 22.603 22.603) 32.007 ======================= Grid point 187 (33/84) ======================= q-point: ( 0.38 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.439 ( 17.347 -2.630 -2.630) 17.741 3.567 ( 16.084 2.622 2.622) 16.506 4.362 ( 5.783 -5.821 -5.821) 10.061 5.573 ( -3.023 7.996 7.996) 11.706 5.806 ( 5.079 17.579 17.579) 25.374 7.235 ( 15.723 4.573 4.573) 17.001 18.377 ( 7.246 -23.019 -23.019) 33.350 18.688 ( 0.021 -18.047 -18.047) 25.522 18.870 ( 18.323 -1.781 -1.781) 18.496 19.619 ( -5.320 14.754 14.754) 21.533 20.986 ( 11.086 25.837 25.837) 38.183 25.399 ( -0.426 -2.558 -2.558) 3.643 34.845 ( -6.085 -30.147 -30.147) 43.067 35.327 ( -2.926 -11.563 -11.563) 16.613 39.168 ( -1.464 23.563 23.563) 33.355 ======================= Grid point 188 (34/84) ======================= q-point: ( 0.46 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.670 ( 6.256 -3.106 -3.106) 7.644 3.779 ( 5.610 1.843 1.843) 6.186 4.420 ( 1.026 -9.893 -9.893) 14.029 5.531 ( -1.216 8.413 8.413) 11.959 5.878 ( 2.084 20.049 20.049) 28.430 7.447 ( 5.771 5.112 5.112) 9.251 18.483 ( 3.215 -33.825 -33.825) 47.943 18.683 ( -0.433 -10.203 -10.203) 14.436 19.193 ( 14.274 -3.678 -3.678) 15.192 19.472 ( -9.502 17.012 17.012) 25.867 21.131 ( 3.861 27.015 27.015) 38.399 25.393 ( -0.162 -2.958 -2.958) 4.186 34.764 ( -2.210 -32.435 -32.435) 45.923 35.287 ( -1.074 -11.778 -11.778) 16.692 39.147 ( -0.573 23.934 23.934) 33.852 ======================= Grid point 197 (35/84) ======================= q-point: ( 0.15 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.077 ( 19.290 19.290 15.701) 31.476 2.360 ( 25.292 25.292 0.650) 35.774 3.606 ( 25.798 25.798 15.885) 39.792 5.653 ( 3.133 3.133 5.092) 6.749 5.896 ( 4.329 4.329 9.409) 11.225 6.196 ( 12.475 12.475 -1.367) 17.696 17.779 ( -6.746 -6.746 2.233) 9.798 18.014 ( 7.855 7.855 -12.264) 16.547 18.526 ( -6.848 -6.848 -27.964) 29.594 19.878 ( 10.062 10.062 32.313) 35.307 20.300 ( 21.993 21.993 10.640) 32.872 25.412 ( -1.033 -1.033 -0.655) 1.601 34.985 ( -21.482 -21.482 -2.047) 30.449 35.436 ( -4.191 -4.191 -23.013) 23.763 39.480 ( 10.736 10.736 16.031) 22.079 ======================= Grid point 198 (36/84) ======================= q-point: ( 0.23 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.526 ( 25.599 8.090 9.935) 28.626 2.858 ( 25.144 16.053 0.423) 29.835 4.007 ( 15.162 11.548 4.984) 19.700 5.720 ( 2.611 4.030 13.745) 14.560 5.974 ( 5.628 19.292 5.480) 20.830 6.537 ( 21.442 6.939 0.496) 22.542 17.652 ( -6.121 -29.281 0.809) 29.925 18.212 ( 12.132 6.606 -18.248) 22.886 18.496 ( 3.122 -1.168 -35.464) 35.620 20.019 ( 4.923 18.308 34.133) 39.044 20.739 ( 22.455 19.926 16.671) 34.339 25.388 ( -1.444 -2.034 -1.256) 2.793 34.591 ( -19.220 -34.131 -4.198) 39.395 35.353 ( -4.452 -4.474 -24.931) 25.718 39.618 ( 4.106 18.148 17.563) 25.587 ======================= Grid point 199 (37/84) ======================= q-point: ( 0.31 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.009 ( 23.489 1.820 4.383) 23.964 3.300 ( 20.201 8.559 1.835) 22.016 4.198 ( 5.346 0.304 -2.891) 6.085 5.734 ( -0.692 4.794 15.896) 16.617 6.129 ( 9.427 28.068 6.883) 30.399 6.973 ( 22.537 7.116 2.386) 23.754 17.554 ( -3.898 -46.517 -3.852) 46.839 18.431 ( 8.996 3.197 -21.340) 23.379 18.645 ( 12.585 4.308 -38.321) 40.564 20.075 ( 0.972 23.474 33.489) 40.908 21.138 ( 18.480 19.438 22.218) 34.828 25.359 ( -1.499 -3.133 -2.108) 4.063 34.249 ( -16.096 -45.087 -7.392) 48.441 35.268 ( -4.224 -4.530 -26.160) 26.883 39.660 ( 0.637 22.505 18.767) 29.310 ======================= Grid point 200 (38/84) ======================= q-point: ( 0.38 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.400 ( 16.351 -0.873 -0.383) 16.379 3.624 ( 13.248 3.050 3.610) 14.066 4.244 ( 0.109 -6.011 -5.881) 8.409 5.713 ( -1.179 7.277 14.885) 16.610 6.298 ( 7.438 30.745 9.285) 32.966 7.361 ( 16.691 8.325 3.487) 18.975 17.501 ( -1.705 -57.073 -8.297) 57.699 18.557 ( 4.693 2.194 -25.719) 26.235 18.930 ( 14.970 4.321 -37.316) 40.438 20.059 ( -2.390 27.698 33.997) 43.917 21.436 ( 12.046 19.451 25.857) 34.526 25.333 ( -1.129 -4.062 -2.922) 5.130 33.983 ( -11.108 -53.310 -10.065) 55.377 35.196 ( -3.062 -4.517 -27.165) 27.708 39.658 ( -0.515 24.572 19.403) 31.313 ======================= Grid point 201 (39/84) ======================= q-point: ( 0.46 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.619 ( 5.909 -1.956 -3.419) 7.101 3.804 ( 5.045 0.764 4.300) 6.673 4.228 ( -1.015 -9.288 -5.734) 10.963 5.696 ( -0.490 8.983 14.186) 16.798 6.399 ( 2.722 31.583 10.270) 33.323 7.586 ( 6.103 9.025 3.942) 11.586 17.482 ( -0.415 -62.067 -10.808) 63.003 18.619 ( 1.623 2.184 -27.347) 27.482 19.158 ( 7.111 1.385 -40.011) 40.661 20.005 ( -2.289 31.929 37.338) 49.182 21.594 ( 4.178 19.526 27.644) 34.102 25.318 ( -0.419 -4.588 -3.416) 5.736 33.834 ( -4.009 -57.831 -11.584) 59.116 35.155 ( -1.121 -4.495 -27.763) 28.147 39.648 ( -0.322 25.337 19.654) 32.068 ======================= Grid point 211 (40/84) ======================= q-point: ( 0.23 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.766 ( 16.060 16.060 9.180) 24.497 3.181 ( 16.495 16.495 0.515) 23.334 4.153 ( 3.846 3.846 -2.510) 5.990 5.787 ( 2.870 2.870 17.113) 17.588 6.388 ( 21.049 21.049 7.777) 30.767 6.744 ( 15.633 15.633 -0.994) 22.131 17.096 ( -26.596 -26.596 -1.392) 37.638 18.391 ( 11.167 11.167 -27.093) 31.360 18.552 ( 6.001 6.001 -43.616) 44.434 20.332 ( 13.527 13.527 42.721) 46.808 21.112 ( 18.829 18.829 18.549) 32.452 25.343 ( -2.564 -2.564 -1.830) 4.062 33.943 ( -32.287 -32.287 -1.951) 45.702 35.268 ( -4.421 -4.421 -28.075) 28.763 39.889 ( 9.937 9.937 16.855) 21.945 ======================= Grid point 212 (41/84) ======================= q-point: ( 0.31 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.104 ( 17.916 7.894 6.120) 20.513 3.469 ( 12.973 8.378 1.945) 15.565 4.177 ( -0.837 -2.255 -5.388) 5.900 5.832 ( 1.814 5.170 17.643) 18.474 6.714 ( 13.124 31.150 6.474) 34.417 7.140 ( 23.125 10.081 1.170) 25.254 16.656 ( -18.773 -44.552 -3.320) 48.459 18.579 ( 7.182 11.196 -35.974) 38.355 18.741 ( 13.689 5.318 -46.752) 49.004 20.531 ( 7.172 22.691 44.039) 50.058 21.456 ( 16.398 13.849 23.762) 32.021 25.293 ( -2.506 -3.582 -2.739) 5.158 33.366 ( -27.063 -44.986 -2.512) 52.559 35.185 ( -4.104 -4.237 -29.815) 30.393 40.017 ( 3.946 14.289 16.970) 22.533 ======================= Grid point 213 (42/84) ======================= q-point: ( 0.38 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.418 ( 13.746 2.906 1.821) 14.167 3.676 ( 8.677 2.173 3.850) 9.738 4.137 ( -3.080 -5.111 -4.666) 7.575 5.860 ( 1.162 7.807 17.022) 18.764 6.921 ( 8.425 32.535 5.732) 34.094 7.547 ( 17.615 10.722 2.387) 20.759 16.373 ( -10.702 -57.846 -4.805) 59.023 18.685 ( 4.118 10.147 -41.857) 43.266 19.017 ( 13.314 4.524 -48.489) 50.487 20.619 ( 2.317 29.192 43.441) 52.390 21.727 ( 11.230 11.040 28.550) 32.605 25.250 ( -1.815 -4.328 -3.660) 5.952 32.920 ( -18.580 -54.984 -3.611) 58.151 35.115 ( -2.946 -4.032 -31.084) 31.483 40.060 ( 0.975 16.700 17.167) 23.970 ======================= Grid point 214 (43/84) ======================= q-point: ( 0.46 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.602 ( 4.845 0.399 -2.556) 5.493 3.809 ( 4.663 -0.210 5.732) 7.392 4.072 ( -2.896 -6.916 -2.586) 7.931 5.876 ( 0.433 9.424 16.479) 18.989 7.033 ( 3.030 33.015 5.700) 33.640 7.784 ( 6.401 11.206 2.593) 13.163 16.238 ( -3.422 -64.894 -5.637) 65.228 18.739 ( 1.439 9.633 -43.967) 45.033 19.208 ( 5.529 3.170 -50.527) 50.927 20.638 ( 0.188 32.884 43.595) 54.607 21.876 ( 3.991 10.086 31.047) 32.887 25.226 ( -0.660 -4.715 -4.235) 6.372 32.672 ( -6.698 -60.577 -4.369) 61.103 35.076 ( -1.075 -3.900 -31.803) 32.060 40.068 ( 0.082 17.680 17.272) 24.716 ======================= Grid point 225 (44/84) ======================= q-point: ( 0.31 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.296 ( 11.269 11.269 5.708) 16.927 3.609 ( 6.026 6.026 2.300) 8.827 4.120 ( -3.463 -3.463 -4.633) 6.741 5.927 ( 4.619 4.619 18.190) 19.327 7.263 ( 21.835 21.835 6.515) 31.560 7.371 ( 16.350 16.350 -0.466) 23.127 15.883 ( -34.905 -34.905 -3.704) 49.502 18.825 ( 10.663 10.663 -42.463) 45.061 18.848 ( 8.432 8.432 -56.292) 57.541 20.914 ( 16.372 16.372 47.689) 53.012 21.690 ( 10.891 10.891 27.190) 31.249 25.226 ( -3.258 -3.258 -3.691) 5.903 32.553 ( -38.446 -38.446 -0.994) 54.380 35.106 ( -3.850 -3.850 -31.673) 32.138 40.217 ( 6.893 6.893 16.350) 19.035 ======================= Grid point 226 (45/84) ======================= q-point: ( 0.38 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.507 ( 9.746 6.033 3.317) 11.932 3.705 ( 4.350 0.845 2.976) 5.338 4.046 ( -4.194 -4.278 -1.915) 6.289 6.004 ( 3.231 6.901 18.141) 19.676 7.512 ( 8.882 27.881 3.236) 29.440 7.763 ( 18.231 11.308 1.730) 21.523 15.331 ( -21.962 -48.949 -4.134) 53.809 18.918 ( 3.285 13.126 -53.646) 55.326 19.098 ( 12.049 3.818 -54.858) 56.295 21.149 ( 8.079 24.957 46.522) 53.408 21.881 ( 8.666 5.653 31.983) 33.615 25.172 ( -2.263 -3.719 -4.644) 6.366 31.918 ( -26.560 -47.759 -0.608) 54.650 35.041 ( -2.732 -3.586 -32.942) 33.249 40.303 ( 2.470 8.769 16.183) 18.571 ======================= Grid point 227 (46/84) ======================= q-point: ( 0.46 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.633 ( 2.909 2.625 -1.097) 4.069 3.798 ( 4.926 -0.888 5.716) 7.598 3.961 ( -4.275 -4.669 0.394) 6.343 6.048 ( 1.167 8.105 18.021) 19.794 7.633 ( 3.323 28.277 2.955) 28.624 8.006 ( 6.526 11.526 1.478) 13.327 15.045 ( -7.477 -57.236 -4.386) 57.888 18.961 ( 1.147 12.524 -55.743) 57.144 19.267 ( 4.771 2.805 -57.281) 57.548 21.243 ( 2.130 29.060 44.900) 53.526 22.004 ( 3.501 3.869 35.110) 35.496 25.142 ( -0.808 -3.939 -5.235) 6.602 31.562 ( -9.604 -53.071 -0.568) 53.936 35.004 ( -0.991 -3.418 -33.640) 33.828 40.329 ( 0.515 9.604 16.154) 18.800 ======================= Grid point 239 (47/84) ======================= q-point: ( 0.38 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.633 ( 6.325 6.325 3.884) 9.752 3.715 ( 0.658 0.658 1.780) 2.008 3.968 ( -3.976 -3.976 1.178) 5.745 6.118 ( 4.665 4.665 18.846) 19.968 7.959 ( 13.418 13.418 -0.054) 18.976 7.973 ( 14.402 14.402 3.001) 20.588 14.530 ( -33.080 -33.080 -4.084) 46.960 19.149 ( 6.654 6.654 -63.274) 63.970 19.168 ( 6.707 6.707 -56.206) 57.001 21.547 ( 15.433 15.433 45.797) 50.732 21.970 ( 4.211 4.211 35.478) 35.974 25.111 ( -2.517 -2.517 -5.567) 6.608 31.125 ( -33.241 -33.241 1.028) 47.021 34.981 ( -2.512 -2.512 -34.141) 34.325 40.416 ( 3.522 3.522 15.807) 16.573 ======================= Grid point 240 (48/84) ======================= q-point: ( 0.46 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.693 ( 0.341 3.129 0.737) 3.232 3.780 ( 5.313 -0.630 4.528) 7.010 3.888 ( -4.360 -3.006 2.840) 6.010 6.179 ( 1.639 5.365 19.210) 20.013 8.094 ( 3.681 18.944 1.181) 19.335 8.213 ( 6.506 9.408 0.787) 11.466 14.090 ( -11.791 -40.294 -4.132) 42.187 19.189 ( 0.835 10.336 -62.341) 63.197 19.311 ( 4.235 1.712 -62.106) 62.274 21.729 ( 4.009 20.598 43.508) 48.305 22.047 ( 2.884 1.063 38.492) 38.614 25.077 ( -0.894 -2.646 -6.121) 6.729 30.680 ( -12.075 -37.128 1.816) 39.084 34.948 ( -0.905 -2.373 -34.782) 34.875 40.455 ( 0.838 4.025 15.693) 16.223 ======================= Grid point 253 (49/84) ======================= q-point: ( 0.46 0.46 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.730 ( 0.538 0.538 3.058) 3.151 3.792 ( 1.998 1.998 3.122) 4.211 3.837 ( -2.422 -2.422 3.497) 4.895 6.250 ( 1.864 1.864 19.874) 20.048 8.342 ( 5.565 5.565 0.200) 7.873 8.349 ( 4.916 4.916 0.578) 6.976 13.548 ( -14.773 -14.773 -4.131) 21.297 19.330 ( 2.476 2.476 -65.270) 65.364 19.335 ( 2.062 2.062 -64.361) 64.427 22.003 ( 6.818 6.818 41.765) 42.864 22.061 ( 0.958 0.958 40.425) 40.448 25.042 ( -0.942 -0.942 -6.638) 6.770 30.181 ( -13.539 -13.539 3.093) 19.395 34.916 ( -0.850 -0.850 -35.382) 35.402 40.501 ( 0.998 0.998 15.543) 15.607 ======================= Grid point 366 (50/84) ======================= q-point: ( 0.15 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.470 ( 17.399 17.399 17.399) 30.136 2.470 ( 17.399 17.399 17.399) 30.136 3.908 ( 13.207 13.207 13.207) 22.875 5.822 ( 12.350 12.350 12.350) 21.390 6.159 ( 7.659 7.659 7.659) 13.266 6.159 ( 7.659 7.659 7.659) 13.266 17.823 ( -2.824 -2.824 -2.824) 4.892 17.823 ( -2.824 -2.824 -2.824) 4.892 18.036 ( -19.952 -19.952 -19.952) 34.558 20.506 ( 20.785 20.785 20.785) 36.001 20.506 ( 20.785 20.785 20.785) 36.001 25.391 ( -1.625 -1.625 -1.625) 2.814 34.929 ( -15.657 -15.657 -15.657) 27.119 34.929 ( -15.657 -15.657 -15.657) 27.119 39.784 ( 13.508 13.508 13.508) 23.396 ======================= Grid point 367 (51/84) ======================= q-point: ( 0.23 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.763 ( 18.033 13.063 13.063) 25.816 2.943 ( 24.369 9.760 9.760) 28.006 4.083 ( 5.138 2.139 2.139) 5.963 6.099 ( 15.319 19.026 19.026) 30.961 6.109 ( -2.727 13.107 13.107) 18.736 6.569 ( 24.182 3.498 3.498) 24.683 17.573 ( -17.401 -19.502 -19.502) 32.610 17.825 ( 5.589 -11.520 -11.520) 17.224 17.967 ( 4.256 -15.710 -15.710) 22.621 20.656 ( 5.874 30.361 30.361) 43.336 21.088 ( 27.984 18.725 18.725) 38.527 25.351 ( -2.432 -2.599 -2.599) 4.408 34.412 ( -25.605 -17.815 -17.815) 35.921 34.845 ( -4.796 -21.817 -21.817) 31.224 39.977 ( 6.713 17.383 17.383) 25.483 ======================= Grid point 368 (52/84) ======================= q-point: ( 0.31 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.117 ( 17.640 6.605 6.605) 19.961 3.379 ( 20.257 6.499 6.499) 22.244 4.131 ( 0.552 -3.776 -3.776) 5.369 6.058 ( -2.420 13.851 13.851) 19.738 6.380 ( 13.271 21.008 21.008) 32.539 7.049 ( 24.314 5.681 5.681) 25.607 17.288 ( -11.875 -29.470 -29.470) 43.335 18.010 ( 12.629 -19.199 -19.199) 29.945 18.056 ( 4.576 -15.815 -15.815) 22.829 20.735 ( 2.354 33.539 33.539) 47.490 21.571 ( 21.766 20.919 20.919) 36.728 25.301 ( -2.608 -3.811 -3.811) 5.987 33.954 ( -21.582 -25.762 -25.762) 42.345 34.750 ( -4.788 -22.106 -22.106) 31.628 40.058 ( 2.156 20.185 20.185) 28.628 ======================= Grid point 369 (53/84) ======================= q-point: ( 0.38 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.411 ( 12.212 1.869 1.869) 12.495 3.705 ( 13.192 4.329 4.329) 14.543 4.125 ( -0.765 -5.745 -5.745) 8.160 6.019 ( -1.604 14.767 14.767) 20.946 6.596 ( 8.859 21.748 21.748) 32.007 7.461 ( 17.483 6.857 6.857) 19.993 17.113 ( -6.469 -35.773 -35.773) 51.002 18.133 ( 3.208 -17.174 -17.174) 24.499 18.267 ( 12.688 -24.169 -24.169) 36.459 20.754 ( -0.116 36.247 36.247) 51.261 21.916 ( 13.828 22.159 22.159) 34.253 25.255 ( -1.995 -4.927 -4.927) 7.248 33.597 ( -14.887 -32.068 -32.068) 47.732 34.668 ( -3.566 -22.399 -22.399) 31.877 40.077 ( 0.168 21.664 21.664) 30.638 ======================= Grid point 370 (54/84) ======================= q-point: ( 0.46 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.573 ( 4.316 -1.093 -1.093) 4.584 3.878 ( 4.556 3.186 3.186) 6.408 4.112 ( -0.430 -5.680 -5.680) 8.044 5.998 ( -0.551 15.443 15.443) 21.846 6.713 ( 3.109 22.161 22.161) 31.494 7.695 ( 6.314 7.240 7.240) 12.029 17.033 ( -2.008 -38.961 -38.961) 55.136 18.176 ( 1.140 -18.118 -18.118) 25.648 18.452 ( 5.401 -26.834 -26.834) 38.331 20.744 ( -0.501 38.189 38.189) 54.009 22.097 ( 4.746 22.715 22.715) 32.472 25.228 ( -0.746 -5.594 -5.594) 7.947 33.398 ( -5.374 -35.616 -35.616) 50.655 34.619 ( -1.323 -22.587 -22.587) 31.970 40.075 ( -0.138 22.241 22.241) 31.453 ======================= Grid point 380 (55/84) ======================= q-point: ( 0.23 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.008 ( 13.092 13.092 15.295) 24.015 3.219 ( 16.728 16.728 3.766) 23.955 4.091 ( -0.936 -0.936 -3.737) 3.965 6.214 ( 3.583 3.583 25.354) 25.855 6.608 ( 22.358 22.358 14.899) 34.953 6.720 ( 16.028 16.028 -1.480) 22.715 17.019 ( -29.139 -29.139 -8.664) 42.110 17.782 ( -1.296 -1.296 -33.355) 33.406 17.827 ( 1.429 1.429 -28.065) 28.138 21.118 ( 16.422 16.422 38.180) 44.689 21.486 ( 23.084 23.084 18.378) 37.464 25.290 ( -3.617 -3.617 -3.644) 6.280 33.889 ( -31.973 -31.973 -3.627) 45.362 34.591 ( -7.936 -7.936 -39.648) 41.206 40.247 ( 10.426 10.426 18.147) 23.382 ======================= Grid point 381 (56/84) ======================= q-point: ( 0.31 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.266 ( 12.928 8.872 10.192) 18.701 3.530 ( 14.795 8.303 4.529) 17.559 4.063 ( -1.551 -3.011 -5.810) 6.726 6.253 ( 0.954 7.671 24.020) 25.233 6.882 ( 11.975 28.121 10.090) 32.187 7.197 ( 25.609 9.739 5.824) 28.011 16.529 ( -21.279 -42.687 -11.321) 49.022 17.818 ( 4.496 -5.534 -40.823) 41.441 17.891 ( 4.551 -1.827 -37.760) 38.078 21.330 ( 7.052 27.370 38.346) 47.637 21.927 ( 21.127 16.318 23.079) 35.288 25.217 ( -3.774 -4.753 -5.037) 7.887 33.286 ( -29.206 -39.802 -6.233) 49.760 34.464 ( -5.643 -9.698 -41.734) 43.216 40.387 ( 4.512 13.485 19.216) 23.905 ======================= Grid point 382 (57/84) ======================= q-point: ( 0.38 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.477 ( 8.384 5.002 4.358) 10.691 3.773 ( 10.205 2.218 5.834) 11.962 4.040 ( -0.850 -2.832 -4.817) 5.652 6.265 ( 0.421 11.149 22.923) 25.494 7.089 ( 8.921 27.430 11.727) 31.137 7.623 ( 17.727 9.795 5.681) 21.035 16.204 ( -12.491 -53.108 -14.128) 56.357 17.893 ( 2.583 -7.820 -38.719) 39.586 18.017 ( 7.840 -5.555 -50.956) 51.854 21.420 ( 2.745 32.248 38.978) 50.664 22.268 ( 13.801 14.417 25.549) 32.420 25.151 ( -2.851 -5.646 -6.394) 8.994 32.799 ( -20.518 -47.010 -9.832) 52.227 34.372 ( -3.888 -10.217 -42.355) 43.744 40.439 ( 1.298 15.295 19.933) 25.159 ======================= Grid point 383 (58/84) ======================= q-point: ( 0.46 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.581 ( 2.513 1.961 0.501) 3.227 3.919 ( 4.683 0.840 5.439) 7.226 4.023 ( -0.996 -3.232 -2.069) 3.965 6.271 ( 0.179 12.989 22.541) 26.017 7.209 ( 3.273 27.663 12.830) 30.669 7.858 ( 6.313 9.484 5.112) 12.487 16.045 ( -4.071 -58.671 -15.781) 60.893 17.924 ( 0.751 -8.141 -39.360) 40.200 18.140 ( 3.791 -8.689 -55.759) 56.559 21.450 ( 0.614 34.478 39.754) 52.626 22.449 ( 4.765 13.912 26.428) 30.244 25.113 ( -1.059 -6.136 -7.238) 9.548 32.524 ( -7.454 -51.247 -12.034) 53.166 34.319 ( -1.412 -10.352 -42.666) 43.927 40.451 ( 0.178 16.062 20.247) 25.845 ======================= Grid point 394 (59/84) ======================= q-point: ( 0.31 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.450 ( 9.572 9.572 9.946) 16.798 3.673 ( 6.164 6.164 4.213) 9.681 4.023 ( -1.033 -1.033 -4.903) 5.116 6.375 ( 5.190 5.190 26.542) 27.538 7.361 ( 16.728 16.728 -0.374) 23.660 7.449 ( 19.645 19.645 12.669) 30.535 15.766 ( -35.081 -35.081 -8.807) 50.388 17.783 ( 1.272 1.272 -52.902) 52.933 17.899 ( 1.822 1.822 -50.022) 50.088 21.759 ( 16.514 16.514 39.593) 45.967 22.222 ( 14.738 14.738 25.745) 33.124 25.125 ( -4.640 -4.640 -6.640) 9.335 32.526 ( -37.350 -37.350 -1.872) 52.854 34.317 ( -6.165 -6.165 -47.431) 48.226 40.579 ( 6.724 6.724 19.070) 21.309 ======================= Grid point 395 (60/84) ======================= q-point: ( 0.38 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.598 ( 4.883 7.264 5.911) 10.562 3.784 ( 6.170 -0.815 5.003) 7.985 4.015 ( -0.322 -0.185 -0.996) 1.063 6.457 ( 3.262 8.732 26.983) 28.548 7.594 ( 9.591 23.719 4.797) 26.031 7.832 ( 16.575 11.702 6.101) 21.188 15.207 ( -22.350 -47.991 -8.834) 53.672 17.825 ( 2.033 0.128 -57.379) 57.415 17.957 ( 4.112 -1.587 -60.928) 61.087 21.976 ( 6.887 24.849 38.692) 46.497 22.486 ( 11.801 8.707 28.593) 32.134 25.046 ( -3.380 -5.150 -8.185) 10.244 31.898 ( -26.527 -44.320 -1.500) 51.674 34.220 ( -3.962 -6.188 -49.392) 49.935 40.663 ( 2.471 8.192 19.154) 20.978 ======================= Grid point 396 (61/84) ======================= q-point: ( 0.46 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.644 ( 0.739 4.396 2.310) 5.020 3.907 ( 5.812 -1.655 5.081) 7.895 3.990 ( -2.670 -0.774 3.000) 4.090 6.500 ( 1.156 10.630 27.055) 29.091 7.726 ( 3.664 24.804 6.933) 26.014 8.052 ( 5.874 10.324 3.384) 12.351 14.916 ( -7.640 -55.735 -9.118) 56.991 17.848 ( 0.519 -0.265 -57.842) 57.845 18.029 ( 2.425 -3.848 -67.983) 68.135 22.055 ( 1.798 27.801 38.366) 47.414 22.647 ( 4.415 7.150 29.690) 30.856 25.000 ( -1.236 -5.378 -9.141) 10.677 31.542 ( -9.678 -48.349 -1.768) 49.340 34.167 ( -1.401 -6.101 -50.154) 50.543 40.689 ( 0.527 8.866 19.224) 21.177 ======================= Grid point 408 (62/84) ======================= q-point: ( 0.38 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.740 ( 5.489 5.489 6.957) 10.424 3.764 ( 0.737 0.737 3.305) 3.466 4.009 ( -0.894 -0.894 3.016) 3.270 6.597 ( 5.727 5.727 29.184) 30.287 7.961 ( 13.672 13.672 0.481) 19.341 8.065 ( 12.016 12.016 6.678) 18.258 14.415 ( -32.970 -32.970 -7.743) 47.265 17.861 ( 2.386 2.386 -68.699) 68.782 17.945 ( 0.508 0.508 -67.434) 67.437 22.354 ( 13.668 13.668 37.087) 41.822 22.629 ( 6.764 6.764 30.803) 32.254 24.960 ( -3.613 -3.613 -9.804) 11.056 31.153 ( -31.560 -31.560 1.807) 44.670 34.124 ( -3.810 -3.810 -51.963) 52.242 40.770 ( 3.323 3.323 18.860) 19.436 ======================= Grid point 409 (63/84) ======================= q-point: ( 0.46 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.739 ( -0.499 5.026 3.990) 6.437 3.877 ( 5.628 -1.092 5.470) 7.923 3.971 ( -3.051 -1.415 5.555) 6.494 6.673 ( 2.004 7.094 30.267) 31.152 8.132 ( 4.465 16.731 2.927) 17.562 8.242 ( 5.029 8.869 2.471) 10.491 13.976 ( -11.807 -40.043 -7.599) 42.434 17.884 ( 0.265 3.136 -70.786) 70.856 17.973 ( 1.459 -2.128 -74.045) 74.090 22.505 ( 3.085 18.432 35.751) 40.341 22.743 ( 3.842 3.146 31.921) 32.305 24.911 ( -1.303 -3.695 -10.762) 11.452 30.729 ( -11.531 -34.450 3.155) 36.465 34.074 ( -1.322 -3.665 -52.926) 53.069 40.807 ( 0.793 3.739 18.771) 19.157 ======================= Grid point 422 (64/84) ======================= q-point: ( 0.46 0.46 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.815 ( 1.997 1.997 5.601) 6.273 3.878 ( 1.755 1.755 5.687) 6.205 3.932 ( -2.627 -2.627 6.065) 7.112 6.766 ( 2.488 2.488 32.091) 32.283 8.352 ( 5.695 5.695 1.054) 8.123 8.371 ( 3.887 3.887 1.878) 5.809 13.434 ( -14.782 -14.782 -7.398) 22.175 17.935 ( 1.137 1.137 -77.118) 77.135 17.948 ( -0.119 -0.119 -76.964) 76.964 22.745 ( 5.670 5.670 33.761) 34.701 22.784 ( 1.751 1.751 32.964) 33.057 24.862 ( -1.325 -1.325 -11.681) 11.830 30.265 ( -12.599 -12.599 5.468) 18.637 34.026 ( -1.256 -1.256 -53.943) 53.972 40.849 ( 0.929 0.929 18.612) 18.658 ======================= Grid point 549 (65/84) ======================= q-point: ( 0.23 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.328 ( 12.172 12.172 12.172) 21.082 3.328 ( 12.172 12.172 12.172) 21.082 4.019 ( -3.450 -3.450 -3.450) 5.975 6.693 ( 10.768 10.768 10.768) 18.650 6.693 ( 10.768 10.768 10.768) 18.650 6.965 ( 21.920 21.920 21.920) 37.966 16.663 ( -26.730 -26.730 -26.730) 46.298 17.385 ( -12.895 -12.895 -12.895) 22.334 17.385 ( -12.895 -12.895 -12.895) 22.334 21.791 ( 22.206 22.206 22.206) 38.462 21.791 ( 22.206 22.206 22.206) 38.462 25.205 ( -4.991 -4.991 -4.991) 8.645 33.807 ( -22.576 -22.576 -22.576) 39.103 33.807 ( -22.576 -22.576 -22.576) 39.103 40.541 ( 11.891 11.891 11.891) 20.597 ======================= Grid point 550 (66/84) ======================= q-point: ( 0.31 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.476 ( 7.324 10.947 10.947) 17.126 3.642 ( 14.909 6.794 6.794) 17.737 3.968 ( -1.530 -3.704 -3.704) 5.457 6.677 ( -0.335 17.783 17.783) 25.152 7.078 ( 17.517 11.163 11.163) 23.581 7.406 ( 22.226 15.649 15.649) 31.365 16.144 ( -24.675 -29.062 -29.062) 47.938 17.182 ( -4.698 -22.861 -22.861) 32.670 17.298 ( -2.230 -23.629 -23.629) 33.490 21.991 ( 7.941 29.633 29.633) 42.654 22.320 ( 23.308 17.810 17.810) 34.317 25.105 ( -5.222 -6.415 -6.415) 10.468 33.068 ( -35.308 -19.848 -19.848) 45.106 33.714 ( -4.747 -31.030 -31.030) 44.139 40.707 ( 5.693 13.425 13.425) 19.820 ======================= Grid point 551 (67/84) ======================= q-point: ( 0.38 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.585 ( 3.664 6.636 6.636) 10.075 3.884 ( 9.860 4.738 4.738) 11.922 3.966 ( 1.191 -2.059 -2.059) 3.146 6.679 ( 0.357 18.925 18.925) 26.766 7.366 ( 11.980 16.363 16.363) 26.058 7.776 ( 15.372 10.403 10.403) 21.278 15.752 ( -15.572 -33.890 -33.890) 50.393 17.168 ( 2.112 -31.879 -31.879) 45.133 17.268 ( -1.094 -27.797 -27.797) 39.327 22.104 ( 4.001 31.684 31.684) 44.986 22.684 ( 14.376 17.548 17.548) 28.679 25.014 ( -3.980 -7.624 -7.624) 11.493 32.481 ( -24.668 -25.054 -25.054) 43.173 33.632 ( -3.576 -30.611 -30.611) 43.437 40.777 ( 1.938 14.455 14.455) 20.534 ======================= Grid point 552 (68/84) ======================= q-point: ( 0.46 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.623 ( 0.739 3.775 3.775) 5.390 3.994 ( 1.104 -0.467 -0.467) 1.286 4.013 ( 3.404 3.682 3.682) 6.221 6.686 ( 0.273 19.706 19.706) 27.870 7.525 ( 4.293 19.264 19.264) 27.580 7.980 ( 5.459 7.477 7.477) 11.901 15.551 ( -5.240 -36.977 -36.977) 52.555 17.216 ( 1.811 -36.606 -36.606) 51.801 17.253 ( -0.394 -29.783 -29.783) 42.121 22.153 ( 1.167 32.846 32.846) 46.466 22.870 ( 4.864 17.287 17.287) 24.927 24.959 ( -1.489 -8.331 -8.331) 11.876 32.150 ( -8.956 -28.039 -28.039) 40.653 33.583 ( -1.336 -30.369 -30.369) 42.968 40.796 ( 0.368 14.924 14.924) 21.108 ======================= Grid point 563 (69/84) ======================= q-point: ( 0.31 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.663 ( 7.432 7.432 11.684) 15.716 3.769 ( 5.847 5.847 5.646) 10.013 3.941 ( 1.634 1.634 -3.516) 4.208 6.898 ( 6.380 6.380 26.697) 28.180 7.362 ( 17.386 17.386 0.569) 24.594 7.738 ( 17.000 17.000 16.688) 29.266 15.516 ( -32.541 -32.541 -18.033) 49.428 16.871 ( -11.056 -11.056 -39.847) 42.805 16.989 ( -9.677 -9.677 -42.742) 44.879 22.435 ( 15.598 15.598 30.548) 37.680 22.649 ( 17.359 17.359 18.413) 30.687 24.979 ( -6.404 -6.404 -8.237) 12.242 32.483 ( -35.731 -35.731 -2.453) 50.591 33.359 ( -9.565 -9.565 -50.260) 52.048 40.902 ( 6.961 6.961 13.775) 16.931 ======================= Grid point 564 (70/84) ======================= q-point: ( 0.38 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.725 ( 0.212 7.645 6.976) 10.352 3.892 ( 6.582 -1.635 5.949) 9.022 4.013 ( 4.264 4.318 0.814) 6.123 6.983 ( 3.120 12.777 26.273) 29.381 7.687 ( 15.105 16.330 5.046) 22.809 8.013 ( 10.761 13.234 12.396) 21.085 14.975 ( -22.358 -42.639 -15.951) 50.719 16.741 ( -2.923 -14.648 -53.618) 55.659 16.859 ( -4.744 -15.149 -50.905) 53.323 22.632 ( 6.328 22.810 29.324) 37.686 22.955 ( 12.983 11.015 19.975) 26.247 24.868 ( -4.764 -7.164 -9.837) 13.069 31.852 ( -27.492 -36.792 -3.797) 46.086 33.234 ( -4.406 -12.595 -50.493) 52.226 40.990 ( 2.667 7.876 14.125) 16.391 ======================= Grid point 565 (71/84) ======================= q-point: ( 0.46 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.720 ( -0.303 5.992 5.355) 8.041 3.988 ( 2.850 -0.043 3.719) 4.686 4.071 ( 1.480 1.961 4.689) 5.294 7.024 ( 1.141 15.466 25.954) 30.234 7.904 ( 6.435 18.679 11.362) 22.791 8.144 ( 2.994 9.326 6.401) 11.700 14.681 ( -7.820 -49.146 -16.064) 52.293 16.727 ( 0.538 -16.549 -62.470) 64.627 16.792 ( -2.162 -18.124 -52.824) 55.889 22.708 ( 1.837 24.586 29.313) 38.302 23.126 ( 4.520 9.647 19.963) 22.628 24.804 ( -1.766 -7.523 -10.798) 13.278 31.480 ( -10.159 -38.305 -5.397) 39.995 33.178 ( -1.450 -13.309 -49.773) 51.542 41.019 ( 0.594 8.318 14.312) 16.565 ======================= Grid point 577 (72/84) ======================= q-point: ( 0.38 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.843 ( 1.576 1.576 4.872) 5.358 3.892 ( 4.402 4.402 8.467) 10.509 4.076 ( 1.943 1.943 3.641) 4.562 7.179 ( 7.688 7.688 29.716) 31.643 7.984 ( 14.162 14.162 2.024) 20.131 8.236 ( 8.694 8.694 10.991) 16.491 14.235 ( -31.917 -31.917 -10.859) 46.425 16.551 ( -5.640 -5.640 -65.821) 66.303 16.648 ( -7.690 -7.690 -65.471) 66.369 22.969 ( 11.638 11.638 26.768) 31.423 23.133 ( 8.086 8.086 21.406) 24.269 24.748 ( -5.113 -5.113 -11.742) 13.790 31.192 ( -29.147 -29.147 2.163) 41.277 33.071 ( -5.414 -5.414 -55.217) 55.745 41.093 ( 3.231 3.231 14.045) 14.770 ======================= Grid point 578 (73/84) ======================= q-point: ( 0.46 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.841 ( -0.270 6.242 6.433) 8.968 3.991 ( 3.773 0.277 6.208) 7.270 4.087 ( -0.553 -0.328 6.009) 6.043 7.279 ( 2.652 10.696 31.008) 32.908 8.209 ( 8.282 12.643 5.037) 15.932 8.323 ( 0.340 8.433 6.425) 10.607 13.804 ( -11.702 -38.783 -9.978) 41.721 16.502 ( 0.002 -7.631 -72.959) 73.356 16.541 ( -3.649 -8.682 -71.285) 71.905 23.095 ( 2.600 15.364 25.483) 29.870 23.259 ( 3.901 4.545 21.673) 22.486 24.679 ( -1.858 -5.203 -12.845) 13.983 30.795 ( -10.889 -30.217 3.551) 32.315 33.005 ( -1.640 -5.709 -55.785) 56.101 41.129 ( 0.777 3.519 14.041) 14.496 ======================= Grid point 591 (74/84) ======================= q-point: ( 0.46 0.46 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.940 ( 3.144 3.144 7.165) 8.432 4.004 ( 1.436 1.436 7.072) 7.358 4.060 ( -2.237 -2.237 6.881) 7.573 7.423 ( 3.911 3.911 34.457) 34.899 8.390 ( 5.976 5.976 3.170) 9.026 8.431 ( 2.016 2.016 4.644) 5.449 13.272 ( -14.711 -14.711 -9.079) 22.700 16.416 ( -1.611 -1.611 -79.108) 79.141 16.431 ( -3.089 -3.089 -78.962) 79.083 23.292 ( 4.531 4.531 23.041) 23.915 23.317 ( 2.040 2.040 22.374) 22.559 24.610 ( -1.861 -1.861 -13.966) 14.212 30.383 ( -11.296 -11.296 6.582) 17.277 32.936 ( -1.686 -1.686 -56.944) 56.994 41.169 ( 0.873 0.873 13.948) 14.003 ======================= Grid point 732 (75/84) ======================= q-point: ( 0.31 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.882 ( 7.160 7.160 7.160) 12.402 3.882 ( 7.160 7.160 7.160) 12.402 3.889 ( 0.168 0.168 0.168) 0.290 7.386 ( 12.905 12.905 12.905) 22.353 7.386 ( 12.905 12.905 12.905) 22.353 8.055 ( 14.941 14.941 14.941) 25.879 15.064 ( -27.321 -27.321 -27.321) 47.322 16.291 ( -24.896 -24.896 -24.896) 43.121 16.291 ( -24.896 -24.896 -24.896) 43.121 22.936 ( 16.525 16.525 16.525) 28.622 22.936 ( 16.525 16.525 16.525) 28.622 24.818 ( -8.054 -8.054 -8.054) 13.950 32.434 ( -23.444 -23.444 -23.444) 40.607 32.434 ( -23.444 -23.444 -23.444) 40.607 41.104 ( 7.416 7.416 7.416) 12.844 ======================= Grid point 733 (76/84) ======================= q-point: ( 0.38 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.856 ( -1.125 6.197 6.197) 8.835 3.991 ( 4.733 3.683 3.683) 7.037 4.056 ( 7.158 4.008 4.008) 9.130 7.433 ( 2.579 19.583 19.583) 27.815 7.810 ( 19.815 8.596 8.596) 23.247 8.274 ( 6.976 13.263 13.263) 20.011 14.559 ( -22.686 -27.393 -27.393) 44.893 15.869 ( -15.088 -34.623 -34.623) 51.236 16.014 ( -9.472 -36.429 -36.429) 52.382 23.107 ( 6.294 19.887 19.887) 28.820 23.265 ( 12.655 12.531 12.531) 21.776 24.680 ( -5.939 -9.275 -9.275) 14.399 31.696 ( -30.961 -15.533 -15.533) 37.962 32.391 ( -1.907 -32.311 -32.311) 45.734 41.200 ( 2.979 7.817 7.817) 11.449 ======================= Grid point 734 (77/84) ======================= q-point: ( 0.46 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.842 ( -0.345 7.006 7.006) 9.915 4.056 ( 1.787 3.236 3.236) 4.913 4.150 ( 2.493 3.270 3.270) 5.253 7.473 ( 1.232 20.529 20.529) 29.059 8.131 ( 12.731 11.129 11.129) 20.244 8.309 ( -3.160 10.524 10.524) 15.215 14.250 ( -8.458 -29.659 -29.659) 42.788 15.687 ( -4.356 -41.659 -41.659) 59.076 15.897 ( -2.922 -39.915 -39.915) 56.523 23.186 ( 1.997 20.299 20.299) 28.777 23.427 ( 4.177 11.584 11.584) 16.906 24.600 ( -2.201 -9.924 -9.924) 14.206 31.282 ( -11.204 -16.837 -16.837) 26.316 32.365 ( -0.714 -31.370 -31.370) 44.369 41.233 ( 0.699 8.032 8.032) 11.381 ======================= Grid point 746 (78/84) ======================= q-point: ( 0.38 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.947 ( 2.920 2.920 5.581) 6.943 4.053 ( 3.365 3.365 7.923) 9.243 4.141 ( 3.867 3.867 2.909) 6.194 7.708 ( 9.192 9.192 24.884) 28.075 8.044 ( 15.072 15.072 4.289) 21.742 8.473 ( 6.380 6.380 12.671) 15.555 13.992 ( -28.001 -28.001 -14.698) 42.240 15.386 ( -18.036 -18.036 -53.024) 58.840 15.478 ( -18.397 -18.397 -54.515) 60.405 23.393 ( 9.404 9.404 17.401) 21.901 23.462 ( 8.714 8.714 13.075) 17.967 24.524 ( -6.623 -6.623 -10.988) 14.439 31.233 ( -26.528 -26.528 2.027) 37.571 32.060 ( -6.636 -6.636 -48.117) 49.023 41.304 ( 3.260 3.260 7.873) 9.123 ======================= Grid point 747 (79/84) ======================= q-point: ( 0.46 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.977 ( 0.609 6.640 7.311) 9.895 4.109 ( 1.803 2.175 5.792) 6.444 4.191 ( 1.283 0.887 4.578) 4.836 7.821 ( 2.867 15.396 24.521) 29.096 8.326 ( 12.575 9.271 7.078) 17.152 8.497 ( -3.944 7.347 11.131) 13.908 13.594 ( -11.290 -34.107 -11.969) 37.868 15.142 ( -6.541 -17.095 -64.977) 67.506 15.272 ( -4.630 -24.520 -59.869) 64.861 23.495 ( 2.259 11.859 16.270) 20.260 23.588 ( 3.553 5.594 12.862) 14.469 24.435 ( -2.409 -6.861 -11.944) 13.984 30.854 ( -10.701 -23.679 2.111) 26.070 31.996 ( -1.238 -9.423 -46.956) 47.908 41.340 ( 0.796 3.417 7.927) 8.669 ======================= Grid point 760 (80/84) ======================= q-point: ( 0.46 0.46 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.082 ( 3.586 3.586 7.214) 8.818 4.140 ( 1.129 1.129 6.747) 6.933 4.185 ( -1.210 -1.210 5.853) 6.098 8.042 ( 6.504 6.504 28.543) 29.988 8.483 ( 6.576 6.576 6.572) 11.387 8.561 ( -0.748 -0.748 9.015) 9.077 13.097 ( -14.434 -14.434 -8.764) 22.215 14.941 ( -5.433 -5.433 -72.242) 72.649 14.961 ( -7.333 -7.333 -71.879) 72.623 23.644 ( 3.401 3.401 13.882) 14.692 23.657 ( 2.100 2.100 13.354) 13.680 24.344 ( -2.429 -2.429 -13.034) 13.479 30.509 ( -9.944 -9.944 6.187) 15.363 31.901 ( -1.900 -1.900 -48.839) 48.913 41.379 ( 0.852 0.852 7.905) 7.996 ======================= Grid point 915 (81/84) ======================= q-point: ( 0.38 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.047 ( 4.554 4.554 4.554) 7.888 4.185 ( 3.556 3.556 3.556) 6.159 4.185 ( 3.556 3.556 3.556) 6.159 8.142 ( 12.851 12.851 12.851) 22.259 8.142 ( 12.851 12.851 12.851) 22.259 8.673 ( 6.703 6.703 6.703) 11.611 13.650 ( -20.365 -20.365 -20.365) 35.273 14.593 ( -31.951 -31.951 -31.951) 55.341 14.593 ( -31.951 -31.951 -31.951) 55.341 23.650 ( 8.171 8.171 8.171) 14.153 23.650 ( 8.171 8.171 8.171) 14.153 24.340 ( -7.764 -7.764 -7.764) 13.447 31.267 ( -15.730 -15.730 -15.730) 27.245 31.267 ( -15.730 -15.730 -15.730) 27.245 41.408 ( 3.341 3.341 3.341) 5.788 ======================= Grid point 916 (82/84) ======================= q-point: ( 0.46 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.111 ( 1.800 6.282 6.282) 9.065 4.204 ( 0.458 3.903 3.903) 5.539 4.257 ( 1.923 2.194 2.194) 3.651 8.228 ( 2.706 18.067 18.067) 25.694 8.480 ( 13.703 8.544 8.544) 18.269 8.700 ( -3.790 8.132 8.132) 12.108 13.313 ( -10.870 -18.570 -18.570) 28.423 14.116 ( -12.673 -37.749 -37.749) 54.869 14.280 ( -6.882 -41.532 -41.532) 59.137 23.734 ( 2.073 8.790 8.790) 12.603 23.770 ( 3.021 6.358 6.358) 9.485 24.236 ( -2.801 -8.318 -8.318) 12.092 30.830 ( -11.728 -7.192 -7.192) 15.524 31.287 ( 0.517 -23.043 -23.043) 32.592 41.446 ( 0.831 3.402 3.402) 4.882 ======================= Grid point 929 (83/84) ======================= q-point: ( 0.46 0.46 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.208 ( 3.369 3.369 5.595) 7.348 4.253 ( 0.895 0.895 4.710) 4.877 4.277 ( 0.132 0.132 3.465) 3.470 8.500 ( 8.419 8.419 19.119) 22.523 8.639 ( 7.763 7.763 9.073) 14.242 8.762 ( -1.967 -1.967 10.477) 10.840 12.946 ( -13.460 -13.460 -6.707) 20.182 13.704 ( -10.372 -10.372 -53.705) 55.673 13.731 ( -12.607 -12.607 -53.475) 56.369 23.842 ( 2.488 2.488 7.049) 7.878 23.847 ( 2.064 2.064 6.704) 7.313 24.126 ( -2.912 -2.912 -9.166) 10.048 30.613 ( -8.849 -8.849 4.381) 13.259 31.104 ( -1.840 -1.840 -32.739) 32.843 41.484 ( 0.856 0.856 3.415) 3.624 ======================= Grid point 1098 (84/84) ======================= q-point: ( 0.46 0.46 0.46) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 7.40e-04 1.64e-04 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.289 ( 2.731 2.731 2.731) 4.730 4.315 ( 1.064 1.064 1.064) 1.843 4.315 ( 1.064 1.064 1.064) 1.843 8.790 ( 7.854 7.854 7.854) 13.603 8.790 ( 7.854 7.854 7.854) 13.603 8.891 ( 1.480 1.480 1.480) 2.564 12.823 ( -7.512 -7.512 -7.512) 13.011 12.992 ( -18.333 -18.333 -18.333) 31.753 12.992 ( -18.333 -18.333 -18.333) 31.753 23.928 ( 2.023 2.023 2.023) 3.504 23.928 ( 2.023 2.023 2.023) 3.504 24.004 ( -3.238 -3.238 -3.238) 5.609 30.671 ( -4.967 -4.967 -4.967) 8.604 30.671 ( -4.967 -4.967 -4.967) 8.604 41.523 ( 0.866 0.866 0.866) 1.500 =================== 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 949.208 949.208 949.208 0.000 0.000 0.000 3/32955 20.0 794.758 794.758 794.758 0.000 0.000 0.000 3/32955 30.0 416.147 416.147 416.147 0.000 0.000 0.000 3/32955 40.0 249.925 249.925 249.925 -0.000 0.000 0.000 3/32955 50.0 160.923 160.923 160.923 -0.000 0.000 0.000 3/32955 60.0 112.544 112.544 112.544 -0.000 0.000 0.000 3/32955 70.0 84.681 84.681 84.681 -0.000 0.000 0.000 3/32955 80.0 67.365 67.365 67.365 -0.000 0.000 0.000 3/32955 90.0 55.776 55.776 55.776 -0.000 0.000 0.000 3/32955 100.0 47.507 47.507 47.507 -0.000 0.000 0.000 3/32955 110.0 41.292 41.292 41.292 -0.000 0.000 0.000 3/32955 120.0 36.432 36.432 36.432 -0.000 0.000 0.000 3/32955 130.0 32.518 32.518 32.518 -0.000 0.000 0.000 3/32955 140.0 29.297 29.297 29.297 -0.000 0.000 0.000 3/32955 150.0 26.605 26.605 26.605 -0.000 0.000 0.000 3/32955 160.0 24.325 24.325 24.325 -0.000 0.000 0.000 3/32955 170.0 22.376 22.376 22.376 -0.000 0.000 0.000 3/32955 180.0 20.695 20.695 20.695 -0.000 0.000 0.000 3/32955 190.0 19.232 19.232 19.232 -0.000 0.000 0.000 3/32955 200.0 17.950 17.950 17.950 -0.000 0.000 0.000 3/32955 210.0 16.819 16.819 16.819 -0.000 0.000 0.000 3/32955 220.0 15.815 15.815 15.815 -0.000 0.000 0.000 3/32955 230.0 14.917 14.917 14.917 -0.000 0.000 0.000 3/32955 240.0 14.111 14.111 14.111 -0.000 0.000 0.000 3/32955 250.0 13.383 13.383 13.383 -0.000 0.000 0.000 3/32955 260.0 12.722 12.722 12.722 -0.000 0.000 0.000 3/32955 270.0 12.120 12.120 12.120 -0.000 0.000 0.000 3/32955 280.0 11.569 11.569 11.569 -0.000 0.000 0.000 3/32955 290.0 11.063 11.063 11.063 -0.000 0.000 0.000 3/32955 300.0 10.597 10.597 10.597 -0.000 0.000 0.000 3/32955 310.0 10.166 10.166 10.166 -0.000 0.000 0.000 3/32955 320.0 9.767 9.767 9.767 -0.000 0.000 0.000 3/32955 330.0 9.396 9.396 9.396 -0.000 0.000 0.000 3/32955 340.0 9.050 9.050 9.050 -0.000 0.000 0.000 3/32955 350.0 8.728 8.728 8.728 -0.000 0.000 0.000 3/32955 360.0 8.427 8.427 8.427 -0.000 0.000 0.000 3/32955 370.0 8.144 8.144 8.144 -0.000 0.000 0.000 3/32955 380.0 7.879 7.879 7.879 -0.000 0.000 0.000 3/32955 390.0 7.630 7.630 7.630 -0.000 0.000 0.000 3/32955 400.0 7.396 7.396 7.396 -0.000 0.000 0.000 3/32955 410.0 7.174 7.174 7.174 -0.000 0.000 0.000 3/32955 420.0 6.966 6.966 6.966 -0.000 0.000 0.000 3/32955 430.0 6.768 6.768 6.768 -0.000 0.000 0.000 3/32955 440.0 6.581 6.581 6.581 -0.000 0.000 0.000 3/32955 450.0 6.404 6.404 6.404 -0.000 0.000 0.000 3/32955 460.0 6.236 6.236 6.236 -0.000 0.000 0.000 3/32955 470.0 6.077 6.077 6.077 -0.000 0.000 0.000 3/32955 480.0 5.925 5.925 5.925 -0.000 0.000 0.000 3/32955 490.0 5.780 5.780 5.780 -0.000 0.000 0.000 3/32955 500.0 5.643 5.643 5.643 -0.000 0.000 0.000 3/32955 510.0 5.511 5.511 5.511 -0.000 0.000 0.000 3/32955 520.0 5.386 5.386 5.386 -0.000 0.000 0.000 3/32955 530.0 5.266 5.266 5.266 -0.000 0.000 0.000 3/32955 540.0 5.151 5.151 5.151 -0.000 0.000 0.000 3/32955 550.0 5.042 5.042 5.042 -0.000 0.000 0.000 3/32955 560.0 4.936 4.936 4.936 -0.000 0.000 0.000 3/32955 570.0 4.836 4.836 4.836 -0.000 0.000 0.000 3/32955 580.0 4.739 4.739 4.739 -0.000 0.000 0.000 3/32955 590.0 4.646 4.646 4.646 -0.000 0.000 0.000 3/32955 600.0 4.557 4.557 4.557 -0.000 0.000 0.000 3/32955 610.0 4.471 4.471 4.471 -0.000 0.000 0.000 3/32955 620.0 4.388 4.388 4.388 -0.000 0.000 0.000 3/32955 630.0 4.308 4.308 4.308 -0.000 0.000 0.000 3/32955 640.0 4.232 4.232 4.232 -0.000 0.000 0.000 3/32955 650.0 4.158 4.158 4.158 -0.000 0.000 0.000 3/32955 660.0 4.086 4.086 4.086 -0.000 0.000 0.000 3/32955 670.0 4.017 4.017 4.017 -0.000 0.000 0.000 3/32955 680.0 3.951 3.951 3.951 -0.000 0.000 0.000 3/32955 690.0 3.886 3.886 3.886 -0.000 0.000 0.000 3/32955 700.0 3.824 3.824 3.824 -0.000 0.000 0.000 3/32955 710.0 3.764 3.764 3.764 -0.000 0.000 0.000 3/32955 720.0 3.705 3.705 3.705 -0.000 0.000 0.000 3/32955 730.0 3.649 3.649 3.649 -0.000 0.000 0.000 3/32955 740.0 3.594 3.594 3.594 -0.000 0.000 0.000 3/32955 750.0 3.541 3.541 3.541 -0.000 0.000 0.000 3/32955 760.0 3.489 3.489 3.489 -0.000 0.000 0.000 3/32955 770.0 3.439 3.439 3.439 -0.000 0.000 0.000 3/32955 780.0 3.391 3.391 3.391 -0.000 0.000 0.000 3/32955 790.0 3.344 3.344 3.344 -0.000 0.000 0.000 3/32955 800.0 3.298 3.298 3.298 -0.000 0.000 0.000 3/32955 810.0 3.253 3.253 3.253 -0.000 0.000 0.000 3/32955 820.0 3.210 3.210 3.210 -0.000 0.000 0.000 3/32955 830.0 3.168 3.168 3.168 -0.000 0.000 0.000 3/32955 840.0 3.127 3.127 3.127 -0.000 0.000 0.000 3/32955 850.0 3.087 3.087 3.087 -0.000 0.000 0.000 3/32955 860.0 3.048 3.048 3.048 -0.000 0.000 0.000 3/32955 870.0 3.010 3.010 3.010 -0.000 0.000 0.000 3/32955 880.0 2.973 2.973 2.973 -0.000 0.000 0.000 3/32955 890.0 2.937 2.937 2.937 -0.000 0.000 0.000 3/32955 900.0 2.902 2.902 2.902 -0.000 0.000 0.000 3/32955 910.0 2.867 2.867 2.867 -0.000 0.000 0.000 3/32955 920.0 2.834 2.834 2.834 -0.000 0.000 0.000 3/32955 930.0 2.801 2.801 2.801 -0.000 0.000 0.000 3/32955 940.0 2.769 2.769 2.769 -0.000 0.000 0.000 3/32955 950.0 2.738 2.738 2.738 -0.000 0.000 0.000 3/32955 960.0 2.708 2.708 2.708 -0.000 0.000 0.000 3/32955 970.0 2.678 2.678 2.678 -0.000 0.000 0.000 3/32955 980.0 2.649 2.649 2.649 -0.000 0.000 0.000 3/32955 990.0 2.620 2.620 2.620 -0.000 0.000 0.000 3/32955 1000.0 2.593 2.593 2.593 -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-08 05:35:24]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|