----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 19:34: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 4.004049210000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.004049210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.004049210000000 Atomic positions (fractional): *1 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 2 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 3 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 *4 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 *5 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.004049210000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.004049210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.004049210000000 Atomic positions (fractional): *1 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 1 2 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 2 3 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 3 *4 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 > 4 *5 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 5 -------------------------------- super cell -------------------------------- Lattice vectors: a 12.012147629999999 0.000000000000000 0.000000000000000 b 0.000000000000000 12.012147629999999 0.000000000000000 c 0.000000000000000 0.000000000000000 12.012147629999999 Atomic positions (fractional): *1 F 0.00000000000000 0.16666666666667 0.00000000000000 18.998 > 1 2 F 0.33333333333333 0.16666666666667 0.00000000000000 18.998 > 1 3 F 0.66666666666667 0.16666666666667 0.00000000000000 18.998 > 1 4 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 1 5 F 0.33333333333333 0.50000000000000 0.00000000000000 18.998 > 1 6 F 0.66666666666667 0.50000000000000 0.00000000000000 18.998 > 1 7 F 0.00000000000000 0.83333333333333 0.00000000000000 18.998 > 1 8 F 0.33333333333333 0.83333333333333 0.00000000000000 18.998 > 1 9 F 0.66666666666667 0.83333333333333 0.00000000000000 18.998 > 1 10 F 0.00000000000000 0.16666666666667 0.33333333333333 18.998 > 1 11 F 0.33333333333333 0.16666666666667 0.33333333333333 18.998 > 1 12 F 0.66666666666667 0.16666666666667 0.33333333333333 18.998 > 1 13 F 0.00000000000000 0.50000000000000 0.33333333333333 18.998 > 1 14 F 0.33333333333333 0.50000000000000 0.33333333333333 18.998 > 1 15 F 0.66666666666667 0.50000000000000 0.33333333333333 18.998 > 1 16 F 0.00000000000000 0.83333333333333 0.33333333333333 18.998 > 1 17 F 0.33333333333333 0.83333333333333 0.33333333333333 18.998 > 1 18 F 0.66666666666667 0.83333333333333 0.33333333333333 18.998 > 1 19 F 0.00000000000000 0.16666666666667 0.66666666666667 18.998 > 1 20 F 0.33333333333333 0.16666666666667 0.66666666666667 18.998 > 1 21 F 0.66666666666667 0.16666666666667 0.66666666666667 18.998 > 1 22 F 0.00000000000000 0.50000000000000 0.66666666666667 18.998 > 1 23 F 0.33333333333333 0.50000000000000 0.66666666666667 18.998 > 1 24 F 0.66666666666667 0.50000000000000 0.66666666666667 18.998 > 1 25 F 0.00000000000000 0.83333333333333 0.66666666666667 18.998 > 1 26 F 0.33333333333333 0.83333333333333 0.66666666666667 18.998 > 1 27 F 0.66666666666667 0.83333333333333 0.66666666666667 18.998 > 1 28 F 0.00000000000000 0.00000000000000 0.16666666666667 18.998 > 2 29 F 0.33333333333333 0.00000000000000 0.16666666666667 18.998 > 2 30 F 0.66666666666667 0.00000000000000 0.16666666666667 18.998 > 2 31 F 0.00000000000000 0.33333333333333 0.16666666666667 18.998 > 2 32 F 0.33333333333333 0.33333333333333 0.16666666666667 18.998 > 2 33 F 0.66666666666667 0.33333333333333 0.16666666666667 18.998 > 2 34 F 0.00000000000000 0.66666666666667 0.16666666666667 18.998 > 2 35 F 0.33333333333333 0.66666666666667 0.16666666666667 18.998 > 2 36 F 0.66666666666667 0.66666666666667 0.16666666666667 18.998 > 2 37 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 2 38 F 0.33333333333333 0.00000000000000 0.50000000000000 18.998 > 2 39 F 0.66666666666667 0.00000000000000 0.50000000000000 18.998 > 2 40 F 0.00000000000000 0.33333333333333 0.50000000000000 18.998 > 2 41 F 0.33333333333333 0.33333333333333 0.50000000000000 18.998 > 2 42 F 0.66666666666667 0.33333333333333 0.50000000000000 18.998 > 2 43 F 0.00000000000000 0.66666666666667 0.50000000000000 18.998 > 2 44 F 0.33333333333333 0.66666666666667 0.50000000000000 18.998 > 2 45 F 0.66666666666667 0.66666666666667 0.50000000000000 18.998 > 2 46 F 0.00000000000000 0.00000000000000 0.83333333333333 18.998 > 2 47 F 0.33333333333333 0.00000000000000 0.83333333333333 18.998 > 2 48 F 0.66666666666667 0.00000000000000 0.83333333333333 18.998 > 2 49 F 0.00000000000000 0.33333333333333 0.83333333333333 18.998 > 2 50 F 0.33333333333333 0.33333333333333 0.83333333333333 18.998 > 2 51 F 0.66666666666667 0.33333333333333 0.83333333333333 18.998 > 2 52 F 0.00000000000000 0.66666666666667 0.83333333333333 18.998 > 2 53 F 0.33333333333333 0.66666666666667 0.83333333333333 18.998 > 2 54 F 0.66666666666667 0.66666666666667 0.83333333333333 18.998 > 2 55 F 0.16666666666667 0.00000000000000 0.00000000000000 18.998 > 3 56 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 3 57 F 0.83333333333333 0.00000000000000 0.00000000000000 18.998 > 3 58 F 0.16666666666667 0.33333333333333 0.00000000000000 18.998 > 3 59 F 0.50000000000000 0.33333333333333 0.00000000000000 18.998 > 3 60 F 0.83333333333333 0.33333333333333 0.00000000000000 18.998 > 3 61 F 0.16666666666667 0.66666666666667 0.00000000000000 18.998 > 3 62 F 0.50000000000000 0.66666666666667 0.00000000000000 18.998 > 3 63 F 0.83333333333333 0.66666666666667 0.00000000000000 18.998 > 3 64 F 0.16666666666667 0.00000000000000 0.33333333333333 18.998 > 3 65 F 0.50000000000000 0.00000000000000 0.33333333333333 18.998 > 3 66 F 0.83333333333333 0.00000000000000 0.33333333333333 18.998 > 3 67 F 0.16666666666667 0.33333333333333 0.33333333333333 18.998 > 3 68 F 0.50000000000000 0.33333333333333 0.33333333333333 18.998 > 3 69 F 0.83333333333333 0.33333333333333 0.33333333333333 18.998 > 3 70 F 0.16666666666667 0.66666666666667 0.33333333333333 18.998 > 3 71 F 0.50000000000000 0.66666666666667 0.33333333333333 18.998 > 3 72 F 0.83333333333333 0.66666666666667 0.33333333333333 18.998 > 3 73 F 0.16666666666667 0.00000000000000 0.66666666666667 18.998 > 3 74 F 0.50000000000000 0.00000000000000 0.66666666666667 18.998 > 3 75 F 0.83333333333333 0.00000000000000 0.66666666666667 18.998 > 3 76 F 0.16666666666667 0.33333333333333 0.66666666666667 18.998 > 3 77 F 0.50000000000000 0.33333333333333 0.66666666666667 18.998 > 3 78 F 0.83333333333333 0.33333333333333 0.66666666666667 18.998 > 3 79 F 0.16666666666667 0.66666666666667 0.66666666666667 18.998 > 3 80 F 0.50000000000000 0.66666666666667 0.66666666666667 18.998 > 3 81 F 0.83333333333333 0.66666666666667 0.66666666666667 18.998 > 3 *82 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 > 4 83 Mg 0.33333333333333 0.00000000000000 0.00000000000000 24.305 > 4 84 Mg 0.66666666666667 0.00000000000000 0.00000000000000 24.305 > 4 85 Mg 0.00000000000000 0.33333333333333 0.00000000000000 24.305 > 4 86 Mg 0.33333333333333 0.33333333333333 0.00000000000000 24.305 > 4 87 Mg 0.66666666666667 0.33333333333333 0.00000000000000 24.305 > 4 88 Mg 0.00000000000000 0.66666666666667 0.00000000000000 24.305 > 4 89 Mg 0.33333333333333 0.66666666666667 0.00000000000000 24.305 > 4 90 Mg 0.66666666666667 0.66666666666667 0.00000000000000 24.305 > 4 91 Mg 0.00000000000000 0.00000000000000 0.33333333333333 24.305 > 4 92 Mg 0.33333333333333 0.00000000000000 0.33333333333333 24.305 > 4 93 Mg 0.66666666666667 0.00000000000000 0.33333333333333 24.305 > 4 94 Mg 0.00000000000000 0.33333333333333 0.33333333333333 24.305 > 4 95 Mg 0.33333333333333 0.33333333333333 0.33333333333333 24.305 > 4 96 Mg 0.66666666666667 0.33333333333333 0.33333333333333 24.305 > 4 97 Mg 0.00000000000000 0.66666666666667 0.33333333333333 24.305 > 4 98 Mg 0.33333333333333 0.66666666666667 0.33333333333333 24.305 > 4 99 Mg 0.66666666666667 0.66666666666667 0.33333333333333 24.305 > 4 100 Mg 0.00000000000000 0.00000000000000 0.66666666666667 24.305 > 4 101 Mg 0.33333333333333 0.00000000000000 0.66666666666667 24.305 > 4 102 Mg 0.66666666666667 0.00000000000000 0.66666666666667 24.305 > 4 103 Mg 0.00000000000000 0.33333333333333 0.66666666666667 24.305 > 4 104 Mg 0.33333333333333 0.33333333333333 0.66666666666667 24.305 > 4 105 Mg 0.66666666666667 0.33333333333333 0.66666666666667 24.305 > 4 106 Mg 0.00000000000000 0.66666666666667 0.66666666666667 24.305 > 4 107 Mg 0.33333333333333 0.66666666666667 0.66666666666667 24.305 > 4 108 Mg 0.66666666666667 0.66666666666667 0.66666666666667 24.305 > 4 *109 K 0.16666666666667 0.16666666666667 0.16666666666667 39.098 > 5 110 K 0.50000000000000 0.16666666666667 0.16666666666667 39.098 > 5 111 K 0.83333333333333 0.16666666666667 0.16666666666667 39.098 > 5 112 K 0.16666666666667 0.50000000000000 0.16666666666667 39.098 > 5 113 K 0.50000000000000 0.50000000000000 0.16666666666667 39.098 > 5 114 K 0.83333333333333 0.50000000000000 0.16666666666667 39.098 > 5 115 K 0.16666666666667 0.83333333333333 0.16666666666667 39.098 > 5 116 K 0.50000000000000 0.83333333333333 0.16666666666667 39.098 > 5 117 K 0.83333333333333 0.83333333333333 0.16666666666667 39.098 > 5 118 K 0.16666666666667 0.16666666666667 0.50000000000000 39.098 > 5 119 K 0.50000000000000 0.16666666666667 0.50000000000000 39.098 > 5 120 K 0.83333333333333 0.16666666666667 0.50000000000000 39.098 > 5 121 K 0.16666666666667 0.50000000000000 0.50000000000000 39.098 > 5 122 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 5 123 K 0.83333333333333 0.50000000000000 0.50000000000000 39.098 > 5 124 K 0.16666666666667 0.83333333333333 0.50000000000000 39.098 > 5 125 K 0.50000000000000 0.83333333333333 0.50000000000000 39.098 > 5 126 K 0.83333333333333 0.83333333333333 0.50000000000000 39.098 > 5 127 K 0.16666666666667 0.16666666666667 0.83333333333333 39.098 > 5 128 K 0.50000000000000 0.16666666666667 0.83333333333333 39.098 > 5 129 K 0.83333333333333 0.16666666666667 0.83333333333333 39.098 > 5 130 K 0.16666666666667 0.50000000000000 0.83333333333333 39.098 > 5 131 K 0.50000000000000 0.50000000000000 0.83333333333333 39.098 > 5 132 K 0.83333333333333 0.50000000000000 0.83333333333333 39.098 > 5 133 K 0.16666666666667 0.83333333333333 0.83333333333333 39.098 > 5 134 K 0.50000000000000 0.83333333333333 0.83333333333333 39.098 > 5 135 K 0.83333333333333 0.83333333333333 0.83333333333333 39.098 > 5 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 2.1091791 0.0000000 0.0000000 0.0000000 2.1091791 0.0000000 0.0000000 0.0000000 2.1091791 -------------------------- Born effective charges -------------------------- 1 F -0.9584263 0.0000000 0.0000000 0.0000000 -1.2734104 0.0000000 0.0000000 0.0000000 -0.9584263 2 F -0.9584263 0.0000000 0.0000000 0.0000000 -0.9584263 0.0000000 0.0000000 0.0000000 -1.2734104 3 F -1.2734104 0.0000000 0.0000000 0.0000000 -0.9584263 0.0000000 0.0000000 0.0000000 -0.9584263 4 Mg 2.0017747 0.0000000 0.0000000 0.0000000 2.0017747 0.0000000 0.0000000 0.0000000 2.0017747 5 K 1.1884883 0.0000000 0.0000000 0.0000000 1.1884883 0.0000000 0.0000000 0.0000000 1.1884883 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 405/405 Permutation basis: 3615/3615 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 104 Number of blocks in projector: 104 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 74 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 30 Use standard eigh solver. Tree of FC basis block matrices: - (104, 100), data: False |-- (30, 30), data: True |-- (74, 70), data: True ----- Solver_atoms: 1 -- 135 / 135 Time (Solver_compr_matrix_reshape): 0.004 Solver_block: 80 / 80 - Time: 0.057 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.061 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (zz) 0.00000000 (zz) Permutation basis: 405/405 Permutation basis: 3615/3615 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 104 Number of blocks in projector: 104 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 74 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 30 Use standard eigh solver. Tree of FC basis block matrices: - (104, 100), data: False |-- (30, 30), data: True |-- (74, 70), data: True Max drift after symmetrization by symfc projector: 0.00000000 (zz) 0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 19:34:08]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 19:34: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 4.004049210000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.004049210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.004049210000000 Atomic positions (fractional): 1 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 2 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 3 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 4 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 5 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 -------------------------------- supercell --------------------------------- Lattice vectors: a 12.012147629999999 0.000000000000000 0.000000000000000 b 0.000000000000000 12.012147629999999 0.000000000000000 c 0.000000000000000 0.000000000000000 12.012147629999999 Atomic positions (fractional): 1 F 0.00000000000000 0.16666666666667 0.00000000000000 18.998 > 1 2 F 0.33333333333333 0.16666666666667 0.00000000000000 18.998 > 1 3 F 0.66666666666667 0.16666666666667 0.00000000000000 18.998 > 1 4 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 1 5 F 0.33333333333333 0.50000000000000 0.00000000000000 18.998 > 1 6 F 0.66666666666667 0.50000000000000 0.00000000000000 18.998 > 1 7 F 0.00000000000000 0.83333333333333 0.00000000000000 18.998 > 1 8 F 0.33333333333333 0.83333333333333 0.00000000000000 18.998 > 1 9 F 0.66666666666667 0.83333333333333 0.00000000000000 18.998 > 1 10 F 0.00000000000000 0.16666666666667 0.33333333333333 18.998 > 1 11 F 0.33333333333333 0.16666666666667 0.33333333333333 18.998 > 1 12 F 0.66666666666667 0.16666666666667 0.33333333333333 18.998 > 1 13 F 0.00000000000000 0.50000000000000 0.33333333333333 18.998 > 1 14 F 0.33333333333333 0.50000000000000 0.33333333333333 18.998 > 1 15 F 0.66666666666667 0.50000000000000 0.33333333333333 18.998 > 1 16 F 0.00000000000000 0.83333333333333 0.33333333333333 18.998 > 1 17 F 0.33333333333333 0.83333333333333 0.33333333333333 18.998 > 1 18 F 0.66666666666667 0.83333333333333 0.33333333333333 18.998 > 1 19 F 0.00000000000000 0.16666666666667 0.66666666666667 18.998 > 1 20 F 0.33333333333333 0.16666666666667 0.66666666666667 18.998 > 1 21 F 0.66666666666667 0.16666666666667 0.66666666666667 18.998 > 1 22 F 0.00000000000000 0.50000000000000 0.66666666666667 18.998 > 1 23 F 0.33333333333333 0.50000000000000 0.66666666666667 18.998 > 1 24 F 0.66666666666667 0.50000000000000 0.66666666666667 18.998 > 1 25 F 0.00000000000000 0.83333333333333 0.66666666666667 18.998 > 1 26 F 0.33333333333333 0.83333333333333 0.66666666666667 18.998 > 1 27 F 0.66666666666667 0.83333333333333 0.66666666666667 18.998 > 1 28 F 0.00000000000000 0.00000000000000 0.16666666666667 18.998 > 28 29 F 0.33333333333333 0.00000000000000 0.16666666666667 18.998 > 28 30 F 0.66666666666667 0.00000000000000 0.16666666666667 18.998 > 28 31 F 0.00000000000000 0.33333333333333 0.16666666666667 18.998 > 28 32 F 0.33333333333333 0.33333333333333 0.16666666666667 18.998 > 28 33 F 0.66666666666667 0.33333333333333 0.16666666666667 18.998 > 28 34 F 0.00000000000000 0.66666666666667 0.16666666666667 18.998 > 28 35 F 0.33333333333333 0.66666666666667 0.16666666666667 18.998 > 28 36 F 0.66666666666667 0.66666666666667 0.16666666666667 18.998 > 28 37 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 28 38 F 0.33333333333333 0.00000000000000 0.50000000000000 18.998 > 28 39 F 0.66666666666667 0.00000000000000 0.50000000000000 18.998 > 28 40 F 0.00000000000000 0.33333333333333 0.50000000000000 18.998 > 28 41 F 0.33333333333333 0.33333333333333 0.50000000000000 18.998 > 28 42 F 0.66666666666667 0.33333333333333 0.50000000000000 18.998 > 28 43 F 0.00000000000000 0.66666666666667 0.50000000000000 18.998 > 28 44 F 0.33333333333333 0.66666666666667 0.50000000000000 18.998 > 28 45 F 0.66666666666667 0.66666666666667 0.50000000000000 18.998 > 28 46 F 0.00000000000000 0.00000000000000 0.83333333333333 18.998 > 28 47 F 0.33333333333333 0.00000000000000 0.83333333333333 18.998 > 28 48 F 0.66666666666667 0.00000000000000 0.83333333333333 18.998 > 28 49 F 0.00000000000000 0.33333333333333 0.83333333333333 18.998 > 28 50 F 0.33333333333333 0.33333333333333 0.83333333333333 18.998 > 28 51 F 0.66666666666667 0.33333333333333 0.83333333333333 18.998 > 28 52 F 0.00000000000000 0.66666666666667 0.83333333333333 18.998 > 28 53 F 0.33333333333333 0.66666666666667 0.83333333333333 18.998 > 28 54 F 0.66666666666667 0.66666666666667 0.83333333333333 18.998 > 28 55 F 0.16666666666667 0.00000000000000 0.00000000000000 18.998 > 55 56 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 55 57 F 0.83333333333333 0.00000000000000 0.00000000000000 18.998 > 55 58 F 0.16666666666667 0.33333333333333 0.00000000000000 18.998 > 55 59 F 0.50000000000000 0.33333333333333 0.00000000000000 18.998 > 55 60 F 0.83333333333333 0.33333333333333 0.00000000000000 18.998 > 55 61 F 0.16666666666667 0.66666666666667 0.00000000000000 18.998 > 55 62 F 0.50000000000000 0.66666666666667 0.00000000000000 18.998 > 55 63 F 0.83333333333333 0.66666666666667 0.00000000000000 18.998 > 55 64 F 0.16666666666667 0.00000000000000 0.33333333333333 18.998 > 55 65 F 0.50000000000000 0.00000000000000 0.33333333333333 18.998 > 55 66 F 0.83333333333333 0.00000000000000 0.33333333333333 18.998 > 55 67 F 0.16666666666667 0.33333333333333 0.33333333333333 18.998 > 55 68 F 0.50000000000000 0.33333333333333 0.33333333333333 18.998 > 55 69 F 0.83333333333333 0.33333333333333 0.33333333333333 18.998 > 55 70 F 0.16666666666667 0.66666666666667 0.33333333333333 18.998 > 55 71 F 0.50000000000000 0.66666666666667 0.33333333333333 18.998 > 55 72 F 0.83333333333333 0.66666666666667 0.33333333333333 18.998 > 55 73 F 0.16666666666667 0.00000000000000 0.66666666666667 18.998 > 55 74 F 0.50000000000000 0.00000000000000 0.66666666666667 18.998 > 55 75 F 0.83333333333333 0.00000000000000 0.66666666666667 18.998 > 55 76 F 0.16666666666667 0.33333333333333 0.66666666666667 18.998 > 55 77 F 0.50000000000000 0.33333333333333 0.66666666666667 18.998 > 55 78 F 0.83333333333333 0.33333333333333 0.66666666666667 18.998 > 55 79 F 0.16666666666667 0.66666666666667 0.66666666666667 18.998 > 55 80 F 0.50000000000000 0.66666666666667 0.66666666666667 18.998 > 55 81 F 0.83333333333333 0.66666666666667 0.66666666666667 18.998 > 55 82 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 > 82 83 Mg 0.33333333333333 0.00000000000000 0.00000000000000 24.305 > 82 84 Mg 0.66666666666667 0.00000000000000 0.00000000000000 24.305 > 82 85 Mg 0.00000000000000 0.33333333333333 0.00000000000000 24.305 > 82 86 Mg 0.33333333333333 0.33333333333333 0.00000000000000 24.305 > 82 87 Mg 0.66666666666667 0.33333333333333 0.00000000000000 24.305 > 82 88 Mg 0.00000000000000 0.66666666666667 0.00000000000000 24.305 > 82 89 Mg 0.33333333333333 0.66666666666667 0.00000000000000 24.305 > 82 90 Mg 0.66666666666667 0.66666666666667 0.00000000000000 24.305 > 82 91 Mg 0.00000000000000 0.00000000000000 0.33333333333333 24.305 > 82 92 Mg 0.33333333333333 0.00000000000000 0.33333333333333 24.305 > 82 93 Mg 0.66666666666667 0.00000000000000 0.33333333333333 24.305 > 82 94 Mg 0.00000000000000 0.33333333333333 0.33333333333333 24.305 > 82 95 Mg 0.33333333333333 0.33333333333333 0.33333333333333 24.305 > 82 96 Mg 0.66666666666667 0.33333333333333 0.33333333333333 24.305 > 82 97 Mg 0.00000000000000 0.66666666666667 0.33333333333333 24.305 > 82 98 Mg 0.33333333333333 0.66666666666667 0.33333333333333 24.305 > 82 99 Mg 0.66666666666667 0.66666666666667 0.33333333333333 24.305 > 82 100 Mg 0.00000000000000 0.00000000000000 0.66666666666667 24.305 > 82 101 Mg 0.33333333333333 0.00000000000000 0.66666666666667 24.305 > 82 102 Mg 0.66666666666667 0.00000000000000 0.66666666666667 24.305 > 82 103 Mg 0.00000000000000 0.33333333333333 0.66666666666667 24.305 > 82 104 Mg 0.33333333333333 0.33333333333333 0.66666666666667 24.305 > 82 105 Mg 0.66666666666667 0.33333333333333 0.66666666666667 24.305 > 82 106 Mg 0.00000000000000 0.66666666666667 0.66666666666667 24.305 > 82 107 Mg 0.33333333333333 0.66666666666667 0.66666666666667 24.305 > 82 108 Mg 0.66666666666667 0.66666666666667 0.66666666666667 24.305 > 82 109 K 0.16666666666667 0.16666666666667 0.16666666666667 39.098 > 109 110 K 0.50000000000000 0.16666666666667 0.16666666666667 39.098 > 109 111 K 0.83333333333333 0.16666666666667 0.16666666666667 39.098 > 109 112 K 0.16666666666667 0.50000000000000 0.16666666666667 39.098 > 109 113 K 0.50000000000000 0.50000000000000 0.16666666666667 39.098 > 109 114 K 0.83333333333333 0.50000000000000 0.16666666666667 39.098 > 109 115 K 0.16666666666667 0.83333333333333 0.16666666666667 39.098 > 109 116 K 0.50000000000000 0.83333333333333 0.16666666666667 39.098 > 109 117 K 0.83333333333333 0.83333333333333 0.16666666666667 39.098 > 109 118 K 0.16666666666667 0.16666666666667 0.50000000000000 39.098 > 109 119 K 0.50000000000000 0.16666666666667 0.50000000000000 39.098 > 109 120 K 0.83333333333333 0.16666666666667 0.50000000000000 39.098 > 109 121 K 0.16666666666667 0.50000000000000 0.50000000000000 39.098 > 109 122 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 109 123 K 0.83333333333333 0.50000000000000 0.50000000000000 39.098 > 109 124 K 0.16666666666667 0.83333333333333 0.50000000000000 39.098 > 109 125 K 0.50000000000000 0.83333333333333 0.50000000000000 39.098 > 109 126 K 0.83333333333333 0.83333333333333 0.50000000000000 39.098 > 109 127 K 0.16666666666667 0.16666666666667 0.83333333333333 39.098 > 109 128 K 0.50000000000000 0.16666666666667 0.83333333333333 39.098 > 109 129 K 0.83333333333333 0.16666666666667 0.83333333333333 39.098 > 109 130 K 0.16666666666667 0.50000000000000 0.83333333333333 39.098 > 109 131 K 0.50000000000000 0.50000000000000 0.83333333333333 39.098 > 109 132 K 0.83333333333333 0.50000000000000 0.83333333333333 39.098 > 109 133 K 0.16666666666667 0.83333333333333 0.83333333333333 39.098 > 109 134 K 0.50000000000000 0.83333333333333 0.83333333333333 39.098 > 109 135 K 0.83333333333333 0.83333333333333 0.83333333333333 39.098 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 2.1091791 0.0000000 0.0000000 0.0000000 2.1091791 0.0000000 0.0000000 0.0000000 2.1091791 -------------------------- Born effective charges -------------------------- 1 F -0.9584263 0.0000000 0.0000000 0.0000000 -1.2734104 0.0000000 0.0000000 0.0000000 -0.9584263 2 F -0.9584263 0.0000000 0.0000000 0.0000000 -0.9584263 0.0000000 0.0000000 0.0000000 -1.2734104 3 F -1.2734104 0.0000000 0.0000000 0.0000000 -0.9584263 0.0000000 0.0000000 0.0000000 -0.9584263 4 Mg 2.0017747 0.0000000 0.0000000 0.0000000 2.0017747 0.0000000 0.0000000 0.0000000 2.0017747 5 K 1.1884883 0.0000000 0.0000000 0.0000000 1.1884883 0.0000000 0.0000000 0.0000000 1.1884883 ---------------------------------------------------------------------------- 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 (zz) -0.00000000 (zz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 19:34:13]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 19:34: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 4.004049210000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.004049210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.004049210000000 Atomic positions (fractional): 1 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 2 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 3 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 4 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 5 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 -------------------------------- supercell --------------------------------- Lattice vectors: a 12.012147629999999 0.000000000000000 0.000000000000000 b 0.000000000000000 12.012147629999999 0.000000000000000 c 0.000000000000000 0.000000000000000 12.012147629999999 Atomic positions (fractional): 1 F 0.00000000000000 0.16666666666667 0.00000000000000 18.998 > 1 2 F 0.33333333333333 0.16666666666667 0.00000000000000 18.998 > 1 3 F 0.66666666666667 0.16666666666667 0.00000000000000 18.998 > 1 4 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 1 5 F 0.33333333333333 0.50000000000000 0.00000000000000 18.998 > 1 6 F 0.66666666666667 0.50000000000000 0.00000000000000 18.998 > 1 7 F 0.00000000000000 0.83333333333333 0.00000000000000 18.998 > 1 8 F 0.33333333333333 0.83333333333333 0.00000000000000 18.998 > 1 9 F 0.66666666666667 0.83333333333333 0.00000000000000 18.998 > 1 10 F 0.00000000000000 0.16666666666667 0.33333333333333 18.998 > 1 11 F 0.33333333333333 0.16666666666667 0.33333333333333 18.998 > 1 12 F 0.66666666666667 0.16666666666667 0.33333333333333 18.998 > 1 13 F 0.00000000000000 0.50000000000000 0.33333333333333 18.998 > 1 14 F 0.33333333333333 0.50000000000000 0.33333333333333 18.998 > 1 15 F 0.66666666666667 0.50000000000000 0.33333333333333 18.998 > 1 16 F 0.00000000000000 0.83333333333333 0.33333333333333 18.998 > 1 17 F 0.33333333333333 0.83333333333333 0.33333333333333 18.998 > 1 18 F 0.66666666666667 0.83333333333333 0.33333333333333 18.998 > 1 19 F 0.00000000000000 0.16666666666667 0.66666666666667 18.998 > 1 20 F 0.33333333333333 0.16666666666667 0.66666666666667 18.998 > 1 21 F 0.66666666666667 0.16666666666667 0.66666666666667 18.998 > 1 22 F 0.00000000000000 0.50000000000000 0.66666666666667 18.998 > 1 23 F 0.33333333333333 0.50000000000000 0.66666666666667 18.998 > 1 24 F 0.66666666666667 0.50000000000000 0.66666666666667 18.998 > 1 25 F 0.00000000000000 0.83333333333333 0.66666666666667 18.998 > 1 26 F 0.33333333333333 0.83333333333333 0.66666666666667 18.998 > 1 27 F 0.66666666666667 0.83333333333333 0.66666666666667 18.998 > 1 28 F 0.00000000000000 0.00000000000000 0.16666666666667 18.998 > 28 29 F 0.33333333333333 0.00000000000000 0.16666666666667 18.998 > 28 30 F 0.66666666666667 0.00000000000000 0.16666666666667 18.998 > 28 31 F 0.00000000000000 0.33333333333333 0.16666666666667 18.998 > 28 32 F 0.33333333333333 0.33333333333333 0.16666666666667 18.998 > 28 33 F 0.66666666666667 0.33333333333333 0.16666666666667 18.998 > 28 34 F 0.00000000000000 0.66666666666667 0.16666666666667 18.998 > 28 35 F 0.33333333333333 0.66666666666667 0.16666666666667 18.998 > 28 36 F 0.66666666666667 0.66666666666667 0.16666666666667 18.998 > 28 37 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 28 38 F 0.33333333333333 0.00000000000000 0.50000000000000 18.998 > 28 39 F 0.66666666666667 0.00000000000000 0.50000000000000 18.998 > 28 40 F 0.00000000000000 0.33333333333333 0.50000000000000 18.998 > 28 41 F 0.33333333333333 0.33333333333333 0.50000000000000 18.998 > 28 42 F 0.66666666666667 0.33333333333333 0.50000000000000 18.998 > 28 43 F 0.00000000000000 0.66666666666667 0.50000000000000 18.998 > 28 44 F 0.33333333333333 0.66666666666667 0.50000000000000 18.998 > 28 45 F 0.66666666666667 0.66666666666667 0.50000000000000 18.998 > 28 46 F 0.00000000000000 0.00000000000000 0.83333333333333 18.998 > 28 47 F 0.33333333333333 0.00000000000000 0.83333333333333 18.998 > 28 48 F 0.66666666666667 0.00000000000000 0.83333333333333 18.998 > 28 49 F 0.00000000000000 0.33333333333333 0.83333333333333 18.998 > 28 50 F 0.33333333333333 0.33333333333333 0.83333333333333 18.998 > 28 51 F 0.66666666666667 0.33333333333333 0.83333333333333 18.998 > 28 52 F 0.00000000000000 0.66666666666667 0.83333333333333 18.998 > 28 53 F 0.33333333333333 0.66666666666667 0.83333333333333 18.998 > 28 54 F 0.66666666666667 0.66666666666667 0.83333333333333 18.998 > 28 55 F 0.16666666666667 0.00000000000000 0.00000000000000 18.998 > 55 56 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 55 57 F 0.83333333333333 0.00000000000000 0.00000000000000 18.998 > 55 58 F 0.16666666666667 0.33333333333333 0.00000000000000 18.998 > 55 59 F 0.50000000000000 0.33333333333333 0.00000000000000 18.998 > 55 60 F 0.83333333333333 0.33333333333333 0.00000000000000 18.998 > 55 61 F 0.16666666666667 0.66666666666667 0.00000000000000 18.998 > 55 62 F 0.50000000000000 0.66666666666667 0.00000000000000 18.998 > 55 63 F 0.83333333333333 0.66666666666667 0.00000000000000 18.998 > 55 64 F 0.16666666666667 0.00000000000000 0.33333333333333 18.998 > 55 65 F 0.50000000000000 0.00000000000000 0.33333333333333 18.998 > 55 66 F 0.83333333333333 0.00000000000000 0.33333333333333 18.998 > 55 67 F 0.16666666666667 0.33333333333333 0.33333333333333 18.998 > 55 68 F 0.50000000000000 0.33333333333333 0.33333333333333 18.998 > 55 69 F 0.83333333333333 0.33333333333333 0.33333333333333 18.998 > 55 70 F 0.16666666666667 0.66666666666667 0.33333333333333 18.998 > 55 71 F 0.50000000000000 0.66666666666667 0.33333333333333 18.998 > 55 72 F 0.83333333333333 0.66666666666667 0.33333333333333 18.998 > 55 73 F 0.16666666666667 0.00000000000000 0.66666666666667 18.998 > 55 74 F 0.50000000000000 0.00000000000000 0.66666666666667 18.998 > 55 75 F 0.83333333333333 0.00000000000000 0.66666666666667 18.998 > 55 76 F 0.16666666666667 0.33333333333333 0.66666666666667 18.998 > 55 77 F 0.50000000000000 0.33333333333333 0.66666666666667 18.998 > 55 78 F 0.83333333333333 0.33333333333333 0.66666666666667 18.998 > 55 79 F 0.16666666666667 0.66666666666667 0.66666666666667 18.998 > 55 80 F 0.50000000000000 0.66666666666667 0.66666666666667 18.998 > 55 81 F 0.83333333333333 0.66666666666667 0.66666666666667 18.998 > 55 82 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 > 82 83 Mg 0.33333333333333 0.00000000000000 0.00000000000000 24.305 > 82 84 Mg 0.66666666666667 0.00000000000000 0.00000000000000 24.305 > 82 85 Mg 0.00000000000000 0.33333333333333 0.00000000000000 24.305 > 82 86 Mg 0.33333333333333 0.33333333333333 0.00000000000000 24.305 > 82 87 Mg 0.66666666666667 0.33333333333333 0.00000000000000 24.305 > 82 88 Mg 0.00000000000000 0.66666666666667 0.00000000000000 24.305 > 82 89 Mg 0.33333333333333 0.66666666666667 0.00000000000000 24.305 > 82 90 Mg 0.66666666666667 0.66666666666667 0.00000000000000 24.305 > 82 91 Mg 0.00000000000000 0.00000000000000 0.33333333333333 24.305 > 82 92 Mg 0.33333333333333 0.00000000000000 0.33333333333333 24.305 > 82 93 Mg 0.66666666666667 0.00000000000000 0.33333333333333 24.305 > 82 94 Mg 0.00000000000000 0.33333333333333 0.33333333333333 24.305 > 82 95 Mg 0.33333333333333 0.33333333333333 0.33333333333333 24.305 > 82 96 Mg 0.66666666666667 0.33333333333333 0.33333333333333 24.305 > 82 97 Mg 0.00000000000000 0.66666666666667 0.33333333333333 24.305 > 82 98 Mg 0.33333333333333 0.66666666666667 0.33333333333333 24.305 > 82 99 Mg 0.66666666666667 0.66666666666667 0.33333333333333 24.305 > 82 100 Mg 0.00000000000000 0.00000000000000 0.66666666666667 24.305 > 82 101 Mg 0.33333333333333 0.00000000000000 0.66666666666667 24.305 > 82 102 Mg 0.66666666666667 0.00000000000000 0.66666666666667 24.305 > 82 103 Mg 0.00000000000000 0.33333333333333 0.66666666666667 24.305 > 82 104 Mg 0.33333333333333 0.33333333333333 0.66666666666667 24.305 > 82 105 Mg 0.66666666666667 0.33333333333333 0.66666666666667 24.305 > 82 106 Mg 0.00000000000000 0.66666666666667 0.66666666666667 24.305 > 82 107 Mg 0.33333333333333 0.66666666666667 0.66666666666667 24.305 > 82 108 Mg 0.66666666666667 0.66666666666667 0.66666666666667 24.305 > 82 109 K 0.16666666666667 0.16666666666667 0.16666666666667 39.098 > 109 110 K 0.50000000000000 0.16666666666667 0.16666666666667 39.098 > 109 111 K 0.83333333333333 0.16666666666667 0.16666666666667 39.098 > 109 112 K 0.16666666666667 0.50000000000000 0.16666666666667 39.098 > 109 113 K 0.50000000000000 0.50000000000000 0.16666666666667 39.098 > 109 114 K 0.83333333333333 0.50000000000000 0.16666666666667 39.098 > 109 115 K 0.16666666666667 0.83333333333333 0.16666666666667 39.098 > 109 116 K 0.50000000000000 0.83333333333333 0.16666666666667 39.098 > 109 117 K 0.83333333333333 0.83333333333333 0.16666666666667 39.098 > 109 118 K 0.16666666666667 0.16666666666667 0.50000000000000 39.098 > 109 119 K 0.50000000000000 0.16666666666667 0.50000000000000 39.098 > 109 120 K 0.83333333333333 0.16666666666667 0.50000000000000 39.098 > 109 121 K 0.16666666666667 0.50000000000000 0.50000000000000 39.098 > 109 122 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 109 123 K 0.83333333333333 0.50000000000000 0.50000000000000 39.098 > 109 124 K 0.16666666666667 0.83333333333333 0.50000000000000 39.098 > 109 125 K 0.50000000000000 0.83333333333333 0.50000000000000 39.098 > 109 126 K 0.83333333333333 0.83333333333333 0.50000000000000 39.098 > 109 127 K 0.16666666666667 0.16666666666667 0.83333333333333 39.098 > 109 128 K 0.50000000000000 0.16666666666667 0.83333333333333 39.098 > 109 129 K 0.83333333333333 0.16666666666667 0.83333333333333 39.098 > 109 130 K 0.16666666666667 0.50000000000000 0.83333333333333 39.098 > 109 131 K 0.50000000000000 0.50000000000000 0.83333333333333 39.098 > 109 132 K 0.83333333333333 0.50000000000000 0.83333333333333 39.098 > 109 133 K 0.16666666666667 0.83333333333333 0.83333333333333 39.098 > 109 134 K 0.50000000000000 0.83333333333333 0.83333333333333 39.098 > 109 135 K 0.83333333333333 0.83333333333333 0.83333333333333 39.098 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 2.1091791 0.0000000 0.0000000 0.0000000 2.1091791 0.0000000 0.0000000 0.0000000 2.1091791 -------------------------- Born effective charges -------------------------- 1 F -0.9584263 0.0000000 0.0000000 0.0000000 -1.2734104 0.0000000 0.0000000 0.0000000 -0.9584263 2 F -0.9584263 0.0000000 0.0000000 0.0000000 -0.9584263 0.0000000 0.0000000 0.0000000 -1.2734104 3 F -1.2734104 0.0000000 0.0000000 0.0000000 -0.9584263 0.0000000 0.0000000 0.0000000 -0.9584263 4 Mg 2.0017747 0.0000000 0.0000000 0.0000000 2.0017747 0.0000000 0.0000000 0.0000000 2.0017747 5 K 1.1884883 0.0000000 0.0000000 0.0000000 1.1884883 0.0000000 0.0000000 0.0000000 1.1884883 ---------------------------------------------------------------------------- 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 (zz) 0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 12 12 12 ] 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.16 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: 0.00e+00 0.00e+00 0.00e+00 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 4.656 ( 0.000 0.000 0.000) 0.000 4.656 ( 0.000 0.000 0.000) 0.000 4.656 ( 0.000 0.000 0.000) 0.000 5.587 ( 0.000 0.000 0.000) 0.000 5.587 ( 0.000 0.000 0.000) 0.000 5.587 ( 0.000 0.000 0.000) 0.000 8.467 ( 0.000 0.000 0.000) 0.000 8.467 ( 0.000 0.000 0.000) 0.000 8.467 ( 0.000 0.000 0.000) 0.000 13.784 ( 0.000 0.000 0.000) 0.000 13.784 ( 0.000 0.000 0.000) 0.000 13.784 ( 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: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.798 ( 37.792 0.000 0.000) 37.792 0.798 ( 37.792 0.000 0.000) 37.792 1.374 ( 63.949 0.000 0.000) 63.949 4.646 ( -1.042 0.000 0.000) 1.042 4.646 ( -1.042 0.000 0.000) 1.042 5.608 ( 1.967 0.000 0.000) 1.967 5.609 ( 2.201 0.000 0.000) 2.201 5.609 ( 2.201 0.000 0.000) 2.201 5.745 ( 11.555 0.000 0.000) 11.555 8.437 ( -2.885 0.000 0.000) 2.885 8.437 ( -2.885 0.000 0.000) 2.885 10.249 ( 4.831 0.000 0.000) 4.831 13.781 ( -0.355 0.000 0.000) 0.355 13.781 ( -0.355 0.000 0.000) 0.355 16.423 ( -9.084 0.000 0.000) 9.084 ======================= Grid point 2 (3/84) ======================= q-point: ( 0.17 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 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.562 ( 35.461 0.000 0.000) 35.461 1.562 ( 35.461 0.000 0.000) 35.461 2.625 ( 55.518 0.000 0.000) 55.518 4.607 ( -2.889 0.000 0.000) 2.889 4.607 ( -2.889 0.000 0.000) 2.889 5.664 ( 3.365 0.000 0.000) 3.365 5.680 ( 4.671 0.000 0.000) 4.671 5.680 ( 4.671 0.000 0.000) 4.671 6.067 ( 18.144 0.000 0.000) 18.144 8.349 ( -5.444 0.000 0.000) 5.444 8.349 ( -5.444 0.000 0.000) 5.444 10.423 ( 12.551 0.000 0.000) 12.551 13.770 ( -0.618 0.000 0.000) 0.618 13.770 ( -0.618 0.000 0.000) 0.618 16.140 ( -18.106 0.000 0.000) 18.106 ======================= Grid point 3 (4/84) ======================= q-point: ( 0.25 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.264 ( 31.797 0.000 0.000) 31.797 2.264 ( 31.797 0.000 0.000) 31.797 3.666 ( 44.357 0.000 0.000) 44.357 4.519 ( -5.683 0.000 0.000) 5.683 4.519 ( -5.683 0.000 0.000) 5.683 5.741 ( 3.822 0.000 0.000) 3.822 5.802 ( 6.894 0.000 0.000) 6.894 5.802 ( 6.894 0.000 0.000) 6.894 6.438 ( 15.898 0.000 0.000) 15.898 8.216 ( -7.226 0.000 0.000) 7.226 8.216 ( -7.226 0.000 0.000) 7.226 10.802 ( 24.506 0.000 0.000) 24.506 13.756 ( -0.719 0.000 0.000) 0.719 13.756 ( -0.719 0.000 0.000) 0.719 15.671 ( -26.858 0.000 0.000) 26.858 ======================= Grid point 4 (5/84) ======================= q-point: ( 0.33 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 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.879 ( 27.173 0.000 0.000) 27.173 2.879 ( 27.173 0.000 0.000) 27.173 4.366 ( -9.096 0.000 0.000) 9.096 4.366 ( -9.096 0.000 0.000) 9.096 4.467 ( 32.511 0.000 0.000) 32.511 5.816 ( 3.257 0.000 0.000) 3.257 5.957 ( 7.596 0.000 0.000) 7.596 5.957 ( 7.596 0.000 0.000) 7.596 6.682 ( 7.206 0.000 0.000) 7.206 8.059 ( -7.478 0.000 0.000) 7.478 8.059 ( -7.478 0.000 0.000) 7.478 11.448 ( 37.018 0.000 0.000) 37.018 13.742 ( -0.626 0.000 0.000) 0.626 13.742 ( -0.626 0.000 0.000) 0.626 15.027 ( -34.811 0.000 0.000) 34.811 ======================= Grid point 5 (6/84) ======================= q-point: ( 0.42 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.392 ( 21.955 0.000 0.000) 21.955 3.392 ( 21.955 0.000 0.000) 21.955 4.139 ( -12.754 0.000 0.000) 12.754 4.139 ( -12.754 0.000 0.000) 12.754 5.006 ( 18.637 0.000 0.000) 18.637 5.871 ( 1.858 0.000 0.000) 1.858 6.098 ( 5.349 0.000 0.000) 5.349 6.098 ( 5.349 0.000 0.000) 5.349 6.751 ( 0.319 0.000 0.000) 0.319 7.923 ( -5.132 0.000 0.000) 5.132 7.923 ( -5.132 0.000 0.000) 5.132 12.300 ( 43.414 0.000 0.000) 43.414 13.731 ( -0.363 0.000 0.000) 0.363 13.731 ( -0.363 0.000 0.000) 0.363 14.238 ( -40.403 0.000 0.000) 40.403 ======================= Grid point 6 (7/84) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.738 ( -0.000 0.000 0.000) 0.000 3.738 ( -0.000 0.000 0.000) 0.000 3.888 ( 0.000 0.000 0.000) 0.000 3.888 ( 0.000 0.000 0.000) 0.000 5.207 ( -0.000 0.000 0.000) 0.000 5.891 ( 0.000 0.000 0.000) 0.000 6.157 ( 0.000 0.000 0.000) 0.000 6.157 ( 0.000 0.000 0.000) 0.000 6.745 ( 0.000 0.000 0.000) 0.000 7.866 ( 0.000 0.000 0.000) 0.000 7.866 ( 0.000 0.000 0.000) 0.000 13.037 ( 0.000 0.000 0.000) 0.000 13.536 ( -0.000 0.000 0.000) 0.000 13.727 ( 0.000 0.000 0.000) 0.000 13.727 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 14 (8/84) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.127 ( 26.632 26.632 0.000) 37.663 1.142 ( 26.333 26.333 0.000) 37.240 1.912 ( 43.711 43.711 0.000) 61.816 4.639 ( -0.755 -0.755 0.000) 1.068 4.745 ( 3.875 3.875 0.000) 5.481 5.531 ( -2.687 -2.687 0.000) 3.800 5.616 ( 0.029 0.029 0.000) 0.041 5.630 ( 2.006 2.006 0.000) 2.837 5.752 ( 7.576 7.576 0.000) 10.715 8.406 ( -2.925 -2.925 0.000) 4.137 8.609 ( 6.881 6.881 0.000) 9.732 10.252 ( 2.512 2.512 0.000) 3.552 13.638 ( -6.901 -6.901 0.000) 9.760 13.777 ( -0.355 -0.355 0.000) 0.502 16.418 ( -4.774 -4.774 0.000) 6.752 ======================= Grid point 15 (9/84) ======================= q-point: ( 0.17 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.750 ( 30.693 16.012 0.000) 34.618 1.751 ( 31.424 16.834 0.000) 35.649 2.900 ( 47.409 24.567 0.000) 53.396 4.609 ( -2.338 0.102 0.000) 2.341 4.737 ( -3.454 10.603 0.000) 11.152 5.484 ( -1.578 -7.894 0.000) 8.050 5.695 ( 4.348 1.471 0.000) 4.590 5.778 ( 11.280 -7.625 0.000) 13.616 5.976 ( 14.776 5.228 0.000) 15.674 8.317 ( -5.531 -3.053 0.000) 6.317 8.643 ( -2.421 21.797 0.000) 21.931 10.398 ( 11.958 -1.570 0.000) 12.061 13.547 ( -2.617 -16.627 0.000) 16.832 13.766 ( -0.618 -0.354 0.000) 0.712 16.195 ( -16.090 1.947 0.000) 16.208 ======================= Grid point 16 (10/84) ======================= q-point: ( 0.25 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.384 ( 29.678 10.290 0.000) 31.411 2.393 ( 29.704 11.735 0.000) 31.938 3.804 ( 38.552 12.412 0.000) 40.501 4.535 ( -4.937 1.456 0.000) 5.148 4.636 ( -6.027 10.582 0.000) 12.178 5.466 ( -0.177 -14.563 0.000) 14.564 5.810 ( 6.594 0.808 0.000) 6.643 5.990 ( 9.007 1.559 0.000) 9.141 6.342 ( 18.437 0.144 0.000) 18.438 8.181 ( -7.365 -3.281 0.000) 8.063 8.541 ( -6.887 26.802 0.000) 27.672 10.769 ( 24.062 -2.960 0.000) 24.243 13.508 ( -1.408 -20.379 0.000) 20.428 13.752 ( -0.718 -0.353 0.000) 0.800 15.764 ( -24.999 6.107 0.000) 25.735 ======================= Grid point 17 (11/84) ======================= q-point: ( 0.33 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.968 ( 26.179 7.464 0.000) 27.222 2.971 ( 25.673 8.412 0.000) 27.016 4.399 ( -8.162 3.221 0.000) 8.775 4.420 ( 4.904 2.755 0.000) 5.624 4.555 ( 13.875 10.560 0.000) 17.436 5.471 ( 0.468 -20.740 0.000) 20.745 5.960 ( 7.448 0.327 0.000) 7.455 6.151 ( 6.213 3.517 0.000) 7.140 6.686 ( 14.073 6.428 0.000) 15.471 8.021 ( -7.663 -3.616 0.000) 8.473 8.373 ( -8.860 26.734 0.000) 28.164 11.393 ( 35.086 -5.028 0.000) 35.444 13.484 ( -0.941 -21.125 0.000) 21.146 13.738 ( -0.625 -0.351 0.000) 0.717 15.168 ( -31.877 10.017 0.000) 33.414 ======================= Grid point 18 (12/84) ======================= q-point: ( 0.42 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.452 ( 20.214 5.261 0.000) 20.888 3.461 ( 20.761 5.504 0.000) 21.478 4.195 ( -11.233 5.722 0.000) 12.606 4.247 ( -11.540 9.854 0.000) 15.175 4.967 ( 15.531 -2.721 0.000) 15.768 5.477 ( -0.049 -25.466 0.000) 25.467 6.099 ( 5.342 0.175 0.000) 5.345 6.240 ( 2.348 2.475 0.000) 3.411 6.922 ( 8.434 16.704 0.000) 18.713 7.881 ( -5.293 -3.990 0.000) 6.628 8.197 ( -7.244 24.447 0.000) 25.498 12.164 ( 36.323 -12.369 0.000) 38.371 13.466 ( -0.784 -18.366 0.000) 18.382 13.727 ( -0.362 -0.349 0.000) 0.503 14.476 ( -32.570 15.476 0.000) 36.060 ======================= Grid point 19 (13/84) ======================= q-point: (-0.50 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.739 ( -0.000 0.094 0.000) 0.094 3.751 ( -0.000 1.222 0.000) 1.222 4.002 ( 0.000 10.705 0.000) 10.705 4.053 ( 0.000 14.373 0.000) 14.373 5.138 ( -0.000 -5.539 0.000) 5.539 5.473 ( 0.000 -27.893 0.000) 27.893 6.159 ( 0.000 0.186 0.000) 0.186 6.261 ( 0.000 2.095 0.000) 2.095 7.016 ( 0.000 21.898 0.000) 21.898 7.822 ( 0.000 -4.177 0.000) 4.177 8.113 ( 0.000 22.753 0.000) 22.753 12.654 ( 0.000 -27.147 0.000) 27.147 13.451 ( 0.000 -7.923 0.000) 7.923 13.723 ( 0.000 -0.348 0.000) 0.348 14.037 ( -0.000 20.091 0.000) 20.091 ======================= Grid point 28 (14/84) ======================= q-point: ( 0.17 0.17 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 301 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.157 ( 22.101 22.101 0.000) 31.256 2.195 ( 24.320 24.320 0.000) 34.393 3.518 ( 31.534 31.534 0.000) 44.596 4.606 ( -0.561 -0.561 0.000) 0.794 4.904 ( 2.008 2.008 0.000) 2.840 5.384 ( -3.245 -3.245 0.000) 4.589 5.618 ( 0.425 0.425 0.000) 0.600 5.745 ( 3.383 3.383 0.000) 4.784 6.226 ( 15.366 15.366 0.000) 21.730 8.224 ( -5.804 -5.804 0.000) 8.208 9.040 ( 13.588 13.588 0.000) 19.216 10.420 ( 5.688 5.688 0.000) 8.044 13.235 ( -11.995 -11.995 0.000) 16.964 13.756 ( -0.616 -0.616 0.000) 0.871 16.125 ( -9.190 -9.190 0.000) 12.997 ======================= Grid point 29 (15/84) ======================= q-point: ( 0.25 0.17 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.643 ( 24.139 14.005 0.000) 27.907 2.714 ( 24.727 18.073 0.000) 30.628 4.113 ( 24.083 15.261 0.000) 28.511 4.578 ( -2.300 2.617 0.000) 3.484 4.860 ( -3.378 4.155 0.000) 5.355 5.247 ( -8.368 -0.480 0.000) 8.382 5.826 ( 13.875 -10.619 0.000) 17.472 5.838 ( 5.591 1.993 0.000) 5.936 6.593 ( 19.572 18.695 0.000) 27.066 8.081 ( -7.809 -6.298 0.000) 10.032 9.163 ( -1.292 30.161 0.000) 30.189 10.699 ( 21.302 -2.803 0.000) 21.486 13.054 ( -5.593 -20.815 0.000) 21.553 13.742 ( -0.715 -0.612 0.000) 0.941 15.820 ( -19.735 -2.099 0.000) 19.846 ======================= Grid point 30 (16/84) ======================= q-point: ( 0.33 0.17 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 518 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.139 ( 22.961 8.201 0.000) 24.381 3.197 ( 21.189 12.297 0.000) 24.499 4.455 ( 7.057 0.778 0.000) 7.100 4.506 ( -4.613 7.226 0.000) 8.573 4.853 ( 4.443 8.258 0.000) 9.377 5.071 ( -8.223 -9.052 0.000) 12.230 5.971 ( 6.888 0.801 0.000) 6.934 6.076 ( 9.418 -6.860 0.000) 11.651 7.011 ( 19.539 21.411 0.000) 28.986 7.910 ( -8.260 -7.033 0.000) 10.849 9.043 ( -8.767 35.564 0.000) 36.628 11.255 ( 30.038 -7.591 0.000) 30.983 12.996 ( 0.104 -23.759 0.000) 23.760 13.728 ( -0.622 -0.608 0.000) 0.870 15.326 ( -26.996 3.733 0.000) 27.253 ======================= Grid point 31 (17/84) ======================= q-point: ( 0.42 0.17 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.570 ( 17.456 4.296 0.000) 17.977 3.574 ( 14.361 5.631 0.000) 15.425 4.393 ( -5.723 13.538 0.000) 14.698 4.448 ( -4.788 4.936 0.000) 6.876 4.924 ( -5.200 -23.548 0.000) 24.116 4.997 ( 6.554 6.304 0.000) 9.094 6.105 ( 5.282 0.334 0.000) 5.293 6.215 ( 4.147 -3.091 0.000) 5.173 7.364 ( 13.350 23.680 0.000) 27.184 7.757 ( -5.832 -7.888 0.000) 9.810 8.851 ( -8.427 36.576 0.000) 37.534 11.836 ( 22.456 -17.611 0.000) 28.538 13.077 ( 7.976 -18.546 0.000) 20.189 13.717 ( -0.360 -0.604 0.000) 0.703 14.744 ( -26.875 9.009 0.000) 28.344 ======================= Grid point 32 (18/84) ======================= q-point: (-0.50 0.17 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 301 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.741 ( -0.000 0.152 0.000) 0.152 3.783 ( 0.000 1.502 0.000) 1.502 4.316 ( 0.000 18.936 0.000) 18.936 4.367 ( 0.000 11.350 0.000) 11.350 4.869 ( 0.000 -29.035 0.000) 29.035 5.070 ( 0.000 2.562 0.000) 2.562 6.164 ( 0.000 0.331 0.000) 0.331 6.257 ( 0.000 -1.616 0.000) 1.616 7.515 ( 0.000 24.460 0.000) 24.460 7.692 ( 0.000 -8.333 0.000) 8.333 8.751 ( 0.000 36.835 0.000) 36.835 12.079 ( 0.000 -26.484 0.000) 26.484 13.224 ( 0.000 -13.410 0.000) 13.410 13.713 ( 0.000 -0.603 0.000) 0.603 14.395 ( -0.000 12.984 0.000) 12.984 ======================= Grid point 42 (19/84) ======================= q-point: ( 0.25 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.957 ( 15.994 15.994 0.000) 22.619 3.102 ( 18.278 18.278 0.000) 25.849 4.377 ( 9.084 9.084 0.000) 12.847 4.639 ( 3.159 3.159 0.000) 4.467 4.728 ( -10.012 -10.012 0.000) 14.159 5.424 ( 4.886 4.886 0.000) 6.909 5.707 ( 4.331 4.331 0.000) 6.125 5.897 ( 3.775 3.775 0.000) 5.339 7.030 ( 22.420 22.420 0.000) 31.706 7.924 ( -8.626 -8.626 0.000) 12.200 9.678 ( 15.812 15.812 0.000) 22.361 10.739 ( 9.819 9.819 0.000) 13.886 12.696 ( -12.934 -12.934 0.000) 18.291 13.728 ( -0.709 -0.709 0.000) 1.003 15.661 ( -13.027 -13.027 0.000) 18.423 ======================= Grid point 43 (20/84) ======================= q-point: ( 0.33 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.299 ( 16.515 7.135 0.000) 17.990 3.444 ( 13.993 10.627 0.000) 17.571 4.443 ( -3.229 -3.042 0.000) 4.436 4.581 ( -3.269 -15.613 0.000) 15.951 4.703 ( 2.797 11.563 0.000) 11.896 5.377 ( -3.938 19.711 0.000) 20.100 5.946 ( 12.404 -4.937 0.000) 13.350 5.996 ( 5.570 1.656 0.000) 5.811 7.490 ( 20.818 23.470 0.000) 31.372 7.732 ( -9.409 -9.899 0.000) 13.658 9.773 ( -3.900 31.950 0.000) 32.187 11.125 ( 23.040 -2.877 0.000) 23.219 12.546 ( -0.454 -18.590 0.000) 18.595 13.714 ( -0.616 -0.703 0.000) 0.935 15.296 ( -21.316 -6.771 0.000) 22.366 ======================= Grid point 44 (21/84) ======================= q-point: ( 0.42 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.618 ( 13.144 -0.543 0.000) 13.155 3.668 ( 7.346 3.331 0.000) 8.066 4.354 ( -3.355 -14.154 0.000) 14.546 4.498 ( -5.489 -11.251 0.000) 12.519 4.750 ( 1.589 20.076 0.000) 20.139 5.327 ( -0.967 19.851 0.000) 19.874 6.113 ( 5.028 0.487 0.000) 5.051 6.145 ( 6.419 -3.411 0.000) 7.269 7.555 ( -6.949 -11.500 0.000) 13.436 7.856 ( 13.352 22.800 0.000) 26.422 9.642 ( -6.594 36.956 0.000) 37.539 11.505 ( 11.135 -11.898 0.000) 16.296 12.705 ( 14.202 -16.987 0.000) 22.142 13.703 ( -0.357 -0.698 0.000) 0.784 14.825 ( -21.794 -1.498 0.000) 21.845 ======================= Grid point 45 (22/84) ======================= q-point: (-0.50 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.745 ( 0.000 0.157 0.000) 0.157 3.776 ( 0.000 -4.314 0.000) 4.314 4.323 ( 0.000 -20.773 0.000) 20.773 4.414 ( 0.000 -4.629 0.000) 4.629 4.766 ( 0.000 23.743 0.000) 23.743 5.320 ( 0.000 18.546 0.000) 18.546 6.172 ( 0.000 0.394 0.000) 0.394 6.211 ( 0.000 -2.570 0.000) 2.570 7.476 ( 0.000 -12.409 0.000) 12.409 8.003 ( 0.000 22.032 0.000) 22.032 9.563 ( 0.000 38.585 0.000) 38.585 11.608 ( 0.000 -16.995 0.000) 16.995 12.919 ( 0.000 -15.251 0.000) 15.251 13.700 ( 0.000 -0.696 0.000) 0.696 14.549 ( -0.000 1.563 0.000) 1.563 ======================= Grid point 56 (23/84) ======================= q-point: ( 0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 301 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.431 ( 5.210 5.210 0.000) 7.368 3.623 ( 6.500 6.500 0.000) 9.192 4.247 ( -10.403 -10.403 0.000) 14.712 4.428 ( -3.777 -3.777 0.000) 5.342 4.966 ( 12.843 12.843 0.000) 18.162 5.688 ( 7.173 7.173 0.000) 10.145 5.919 ( 4.711 4.711 0.000) 6.663 6.045 ( 3.167 3.167 0.000) 4.479 7.507 ( -11.362 -11.362 0.000) 16.069 7.950 ( 19.913 19.913 0.000) 28.161 10.213 ( 8.030 8.030 0.000) 11.356 11.246 ( 14.453 14.453 0.000) 20.439 12.270 ( -5.993 -5.993 0.000) 8.476 13.700 ( -0.610 -0.610 0.000) 0.863 15.054 ( -15.855 -15.855 0.000) 22.422 ======================= Grid point 57 (24/84) ======================= q-point: ( 0.42 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.502 ( 1.626 -10.428 0.000) 10.554 3.718 ( 2.881 1.677 0.000) 3.333 4.116 ( -2.162 -5.326 0.000) 5.748 4.341 ( -4.270 -6.128 0.000) 7.469 5.195 ( 8.219 21.828 0.000) 23.325 5.719 ( -0.651 16.605 0.000) 16.617 6.081 ( 7.082 -2.510 0.000) 7.513 6.126 ( 4.176 0.768 0.000) 4.247 7.284 ( -9.163 -14.248 0.000) 16.940 8.287 ( 11.831 18.128 0.000) 21.647 10.225 ( -2.923 14.730 0.000) 15.018 11.483 ( 4.371 13.039 0.000) 13.752 12.386 ( 15.382 -13.033 0.000) 20.161 13.689 ( -0.353 -0.604 0.000) 0.700 14.684 ( -17.663 -11.598 0.000) 21.130 ======================= Grid point 58 (25/84) ======================= q-point: (-0.50 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 301 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.514 ( 0.000 -18.228 0.000) 18.228 3.748 ( 0.000 0.121 0.000) 0.121 4.111 ( 0.000 -1.741 0.000) 1.741 4.284 ( 0.000 -6.475 0.000) 6.475 5.283 ( 0.000 25.418 0.000) 25.418 5.710 ( 0.000 17.251 0.000) 17.251 6.158 ( 0.000 -2.405 0.000) 2.405 6.180 ( 0.000 0.352 0.000) 0.352 7.177 ( 0.000 -16.244 0.000) 16.244 8.415 ( 0.000 17.156 0.000) 17.156 10.185 ( 0.000 16.265 0.000) 16.265 11.505 ( 0.000 10.927 0.000) 10.927 12.615 ( 0.000 -13.334 0.000) 13.334 13.686 ( 0.000 -0.602 0.000) 0.602 14.461 ( 0.000 -9.669 0.000) 9.669 ======================= Grid point 70 (26/84) ======================= q-point: ( 0.42 0.42 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.274 ( -9.205 -9.205 0.000) 13.018 3.742 ( 0.707 0.707 0.000) 1.000 4.091 ( 1.190 1.190 0.000) 1.683 4.231 ( -4.374 -4.374 0.000) 6.186 5.613 ( 16.778 16.778 0.000) 23.728 5.972 ( 5.714 5.714 0.000) 8.081 6.060 ( 2.078 2.078 0.000) 2.939 6.150 ( 1.789 1.789 0.000) 2.530 6.982 ( -13.740 -13.740 0.000) 19.431 8.589 ( 10.316 10.316 0.000) 14.589 10.305 ( -1.994 -1.994 0.000) 2.820 11.897 ( 15.516 15.516 0.000) 21.943 12.228 ( 2.400 2.400 0.000) 3.394 13.679 ( -0.349 -0.349 0.000) 0.494 14.384 ( -15.255 -15.255 0.000) 21.574 ======================= Grid point 71 (27/84) ======================= q-point: (-0.50 0.42 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.173 ( 0.000 -12.554 0.000) 12.554 3.750 ( 0.000 0.063 0.000) 0.063 4.119 ( 0.000 1.073 0.000) 1.073 4.173 ( -0.000 -3.863 0.000) 3.863 5.805 ( 0.000 24.242 0.000) 24.242 5.998 ( 0.000 9.781 0.000) 9.781 6.117 ( 0.000 -1.430 0.000) 1.430 6.186 ( 0.000 0.207 0.000) 0.207 6.805 ( 0.000 -19.177 0.000) 19.177 8.699 ( 0.000 9.671 0.000) 9.671 10.272 ( 0.000 -2.381 0.000) 2.381 11.972 ( 0.000 26.272 0.000) 26.272 12.388 ( 0.000 -7.925 0.000) 7.925 13.676 ( 0.000 -0.348 0.000) 0.348 14.187 ( 0.000 -14.768 0.000) 14.768 ======================= Grid point 90 (28/84) ======================= q-point: (-0.50 -0.50 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.037 ( 0.000 0.000 0.000) 0.000 3.750 ( 0.000 0.000 0.000) 0.000 4.132 ( 0.000 0.000 0.000) 0.000 4.132 ( 0.000 0.000 0.000) 0.000 6.102 ( 0.000 0.000 0.000) 0.000 6.102 ( 0.000 0.000 0.000) 0.000 6.188 ( -0.000 -0.000 0.000) 0.000 6.188 ( 0.000 0.000 0.000) 0.000 6.476 ( 0.000 0.000 0.000) 0.000 8.802 ( 0.000 0.000 0.000) 0.000 10.233 ( 0.000 0.000 0.000) 0.000 12.303 ( 0.000 0.000 0.000) 0.000 12.303 ( 0.000 0.000 0.000) 0.000 13.672 ( 0.000 0.000 0.000) 0.000 13.990 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 183 (29/84) ======================= q-point: ( 0.08 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 182 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.377 ( 20.851 20.851 20.851) 36.116 1.377 ( 20.851 20.851 20.851) 36.116 2.326 ( 35.072 35.072 35.072) 60.746 4.730 ( 1.908 1.908 1.908) 3.305 4.730 ( 1.908 1.908 1.908) 3.305 5.456 ( -4.148 -4.148 -4.148) 7.184 5.548 ( -1.828 -1.828 -1.828) 3.166 5.749 ( 5.145 5.145 5.145) 8.911 5.749 ( 5.145 5.145 5.145) 8.911 8.577 ( 3.550 3.550 3.550) 6.149 8.577 ( 3.550 3.550 3.550) 6.149 10.253 ( 1.551 1.551 1.551) 2.686 13.634 ( -4.689 -4.689 -4.689) 8.122 13.634 ( -4.689 -4.689 -4.689) 8.122 16.415 ( -3.144 -3.144 -3.144) 5.445 ======================= Grid point 184 (30/84) ======================= q-point: ( 0.17 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.893 ( 27.463 13.840 13.840) 33.724 1.909 ( 28.183 13.487 13.487) 34.031 3.157 ( 40.937 22.520 22.520) 51.866 4.682 ( -3.198 3.124 3.124) 5.454 4.762 ( -1.827 4.819 4.819) 7.056 5.385 ( -2.628 -7.062 -7.062) 10.328 5.618 ( 5.154 -5.779 -5.779) 9.662 5.893 ( 13.711 2.001 2.001) 14.000 5.945 ( 10.332 6.419 6.419) 13.753 8.486 ( -5.686 6.647 6.647) 10.987 8.691 ( 0.700 9.756 9.756) 13.815 10.375 ( 11.125 -1.487 -1.487) 11.322 13.491 ( -4.794 -8.639 -8.639) 13.124 13.625 ( -0.537 -6.820 -6.820) 9.659 16.234 ( -14.153 1.340 1.340) 14.279 ======================= Grid point 185 (31/84) ======================= q-point: ( 0.25 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.473 ( 27.572 9.075 9.075) 30.413 2.504 ( 28.282 9.247 9.247) 31.159 3.933 ( 32.239 11.058 11.058) 35.832 4.592 ( -5.581 2.984 2.984) 6.997 4.700 ( -3.291 8.243 8.243) 12.114 5.347 ( -1.079 -10.227 -10.227) 14.503 5.696 ( 2.805 -6.436 -6.436) 9.524 6.171 ( 11.141 9.122 9.122) 17.046 6.277 ( 20.636 -0.367 -0.367) 20.643 8.346 ( -7.645 6.334 6.334) 11.777 8.629 ( -6.108 14.631 14.631) 21.574 10.736 ( 23.621 -3.049 -3.049) 24.011 13.419 ( -2.444 -12.374 -12.374) 17.669 13.613 ( -0.623 -6.736 -6.736) 9.547 15.841 ( -23.200 5.159 5.159) 24.321 ======================= Grid point 186 (32/84) ======================= q-point: ( 0.33 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.021 ( 24.712 6.147 6.147) 26.197 3.062 ( 24.767 6.634 6.634) 26.485 4.398 ( 7.273 -0.668 -0.668) 7.334 4.445 ( -8.641 3.289 3.289) 9.813 4.711 ( 9.482 13.887 13.887) 21.808 5.337 ( 0.039 -13.113 -13.113) 18.545 5.741 ( 1.411 -10.213 -10.213) 14.512 6.398 ( 10.373 11.401 11.401) 19.172 6.683 ( 17.600 3.934 3.934) 18.458 8.178 ( -8.139 5.838 5.838) 11.593 8.461 ( -9.412 15.198 15.198) 23.464 11.341 ( 33.495 -4.824 -4.824) 34.183 13.381 ( -1.253 -13.521 -13.521) 19.163 13.601 ( -0.541 -6.651 -6.651) 9.422 15.287 ( -29.484 8.774 8.774) 31.989 ======================= Grid point 187 (33/84) ======================= q-point: ( 0.42 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.490 ( 20.041 4.168 4.168) 20.890 3.513 ( 17.783 3.622 3.622) 18.506 4.232 ( -11.656 4.031 4.031) 12.975 4.343 ( -6.831 7.148 7.148) 12.200 4.990 ( 11.814 3.587 3.587) 12.857 5.343 ( 0.344 -15.102 -15.102) 21.360 5.749 ( -0.550 -13.486 -13.486) 19.080 6.584 ( 6.911 13.150 13.150) 19.839 6.991 ( 11.427 9.601 9.601) 17.746 8.026 ( -5.852 5.167 5.167) 9.362 8.266 ( -8.348 13.954 13.954) 21.428 12.051 ( 31.647 -10.435 -10.435) 34.919 13.366 ( -0.159 -11.745 -11.745) 16.610 13.591 ( -0.313 -6.588 -6.588) 9.321 14.662 ( -28.282 13.302 13.302) 33.968 ======================= Grid point 188 (34/84) ======================= q-point: (-0.50 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 276 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.729 ( 0.000 -0.220 -0.220) 0.311 3.780 ( -0.000 1.939 1.939) 2.742 4.030 ( 0.000 5.899 5.899) 8.343 4.239 ( 0.000 12.689 12.689) 17.946 5.122 ( -0.000 -0.117 -0.117) 0.166 5.347 ( 0.000 -15.803 -15.803) 22.349 5.737 ( 0.000 -14.592 -14.592) 20.637 6.661 ( 0.000 13.909 13.909) 19.670 7.120 ( 0.000 12.389 12.389) 17.521 7.960 ( 0.000 4.764 4.764) 6.738 8.167 ( 0.000 12.899 12.899) 18.242 12.445 ( 0.000 -18.211 -18.211) 25.754 13.370 ( 0.000 -7.619 -7.619) 10.775 13.588 ( 0.000 -6.564 -6.564) 9.283 14.301 ( -0.000 17.181 17.181) 24.297 ======================= Grid point 197 (35/84) ======================= q-point: ( 0.17 0.17 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.262 ( 20.247 20.247 9.832) 30.275 2.285 ( 21.918 21.918 8.574) 32.161 3.706 ( 27.008 27.008 16.060) 41.434 4.686 ( -0.957 -0.957 6.516) 6.656 4.877 ( 2.081 2.081 -2.558) 3.899 5.255 ( -5.077 -5.077 -8.919) 11.450 5.567 ( 2.410 2.410 -2.127) 4.017 5.919 ( 4.794 4.794 13.146) 14.792 6.208 ( 15.260 15.260 -1.122) 21.611 8.459 ( -5.230 -5.230 19.746) 21.085 9.005 ( 13.477 13.477 -3.345) 19.351 10.394 ( 5.236 5.236 -2.070) 7.689 13.232 ( -11.974 -11.974 -0.214) 16.935 13.570 ( -1.078 -1.078 -15.415) 15.490 16.169 ( -8.326 -8.326 2.391) 12.015 ======================= Grid point 198 (36/84) ======================= q-point: ( 0.25 0.17 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.715 ( 22.719 12.991 6.853) 27.053 2.768 ( 23.601 16.039 5.071) 28.982 4.176 ( 16.493 10.242 5.108) 20.075 4.651 ( -2.925 3.709 4.393) 6.450 4.880 ( 1.107 0.910 2.494) 2.876 5.154 ( -4.706 -2.764 -8.422) 10.036 5.630 ( 1.433 2.104 -3.865) 4.628 6.140 ( 15.094 -2.953 16.588) 22.621 6.569 ( 19.398 19.453 -1.989) 27.544 8.319 ( -7.973 -5.454 20.386) 22.559 9.139 ( -0.700 28.337 -1.907) 28.409 10.663 ( 20.880 -2.995 -3.321) 21.354 13.041 ( -6.373 -19.480 -1.503) 20.551 13.551 ( -0.843 -1.384 -15.999) 16.081 15.888 ( -18.298 -1.961 4.839) 19.029 ======================= Grid point 199 (37/84) ======================= q-point: ( 0.33 0.17 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.184 ( 21.773 7.891 4.413) 23.576 3.235 ( 20.607 10.511 3.371) 23.377 4.378 ( 4.224 -1.000 -3.282) 5.441 4.551 ( -6.586 6.311 0.727) 9.150 4.939 ( 1.193 -13.951 0.461) 14.010 5.095 ( 1.278 11.014 6.077) 12.644 5.630 ( -0.847 0.338 -13.475) 13.506 6.453 ( 13.776 -2.210 20.641) 24.914 6.992 ( 20.170 21.615 -1.471) 29.601 8.140 ( -8.870 -6.064 19.876) 22.595 9.025 ( -8.740 33.644 -1.143) 34.780 11.208 ( 29.271 -7.378 -4.399) 30.505 12.967 ( -0.583 -22.328 -3.126) 22.553 13.534 ( -0.849 -1.704 -15.993) 16.106 15.430 ( -24.931 3.530 7.897) 26.389 ======================= Grid point 200 (38/84) ======================= q-point: ( 0.42 0.17 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.587 ( 14.085 3.269 1.672) 14.555 3.604 ( 16.053 5.927 2.754) 17.332 4.373 ( -7.748 9.486 -0.621) 12.264 4.432 ( -0.348 0.140 -2.104) 2.137 4.890 ( -3.618 -24.665 -3.920) 25.236 5.187 ( 4.629 14.941 15.932) 22.327 5.609 ( -0.920 -1.957 -24.199) 24.296 6.691 ( 8.597 0.168 23.829) 25.333 7.361 ( 14.096 22.567 0.141) 26.608 7.970 ( -6.814 -6.928 18.811) 21.173 8.831 ( -8.627 35.006 -1.360) 36.079 11.766 ( 21.270 -15.661 -6.576) 27.220 13.030 ( 6.697 -17.965 -4.651) 19.729 13.514 ( -1.154 -2.353 -15.822) 16.037 14.904 ( -23.398 8.472 11.974) 27.616 ======================= Grid point 201 (39/84) ======================= q-point: (-0.50 0.17 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.732 ( 0.000 0.500 -0.128) 0.516 3.822 ( -0.000 1.807 2.911) 3.426 4.244 ( 0.000 13.840 -2.143) 14.005 4.436 ( 0.000 3.671 1.485) 3.960 4.844 ( 0.000 -27.783 -2.701) 27.914 5.238 ( 0.000 13.269 15.886) 20.699 5.598 ( 0.000 -2.740 -28.839) 28.969 6.786 ( 0.000 1.205 25.339) 25.368 7.521 ( 0.000 22.772 1.021) 22.795 7.891 ( 0.000 -7.572 18.039) 19.563 8.728 ( 0.000 35.556 -1.670) 35.596 11.996 ( 0.000 -22.806 -7.996) 24.167 13.151 ( 0.000 -12.926 -6.836) 14.622 13.493 ( -0.000 -3.195 -15.444) 15.771 14.614 ( -0.000 11.632 15.397) 19.297 ======================= Grid point 211 (40/84) ======================= q-point: ( 0.25 0.25 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.009 ( 15.127 15.127 4.890) 21.945 3.130 ( 17.717 17.717 2.681) 25.199 4.322 ( 3.986 3.986 -2.522) 6.175 4.710 ( -10.028 -10.028 -1.748) 14.290 4.742 ( 2.923 2.923 4.272) 5.945 5.190 ( 2.026 2.026 -16.230) 16.481 5.784 ( 7.761 7.761 10.601) 15.259 6.149 ( 6.096 6.096 17.920) 19.886 7.011 ( 22.454 22.454 -1.688) 31.799 8.169 ( -8.574 -8.574 21.823) 24.965 9.639 ( 15.774 15.774 -3.712) 22.614 10.697 ( 9.490 9.490 -3.881) 13.971 12.693 ( -12.989 -12.989 -0.231) 18.371 13.527 ( -1.119 -1.119 -17.417) 17.489 15.741 ( -12.042 -12.042 6.070) 18.080 ======================= Grid point 212 (41/84) ======================= q-point: ( 0.33 0.25 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.334 ( 15.722 6.681 3.330) 17.404 3.469 ( 14.177 10.838 2.277) 17.990 4.361 ( -0.080 -0.825 -5.571) 5.632 4.531 ( -6.019 -15.146 -2.590) 16.502 4.740 ( -3.663 6.492 -1.703) 7.646 5.272 ( 6.344 2.537 -14.249) 15.802 5.819 ( -2.658 16.073 10.282) 19.265 6.398 ( 14.553 -2.396 21.612) 26.166 7.470 ( 20.732 23.209 -1.602) 31.161 7.975 ( -9.680 -9.374 21.672) 25.520 9.736 ( -3.854 31.724 -3.534) 32.152 11.079 ( 23.095 -3.043 -4.284) 23.685 12.534 ( -1.361 -18.204 -1.166) 18.292 13.502 ( -1.284 -1.489 -18.094) 18.201 15.405 ( -19.490 -6.074 8.435) 22.088 ======================= Grid point 213 (42/84) ======================= q-point: ( 0.42 0.25 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.630 ( 11.222 -0.176 1.432) 11.314 3.704 ( 8.682 3.331 3.055) 9.788 4.330 ( -1.845 -11.614 -2.270) 11.977 4.447 ( -3.631 -12.085 -3.903) 13.209 4.643 ( -3.860 15.157 -8.003) 17.570 5.427 ( 7.830 3.364 -13.033) 15.572 5.730 ( -5.458 15.991 7.664) 18.554 6.667 ( 10.202 -1.887 25.031) 27.096 7.770 ( -4.582 -6.477 17.246) 18.983 7.851 ( 10.123 18.018 2.038) 20.767 9.605 ( -6.586 36.732 -3.473) 37.479 11.466 ( 11.566 -11.059 -3.785) 16.444 12.663 ( 12.526 -16.805 -3.832) 21.308 13.469 ( -1.955 -2.094 -19.022) 19.236 14.984 ( -18.755 -0.987 12.015) 22.296 ======================= Grid point 214 (43/84) ======================= q-point: (-0.50 0.25 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.740 ( 0.000 -0.676 0.394) 0.782 3.823 ( 0.000 -3.280 3.485) 4.786 4.313 ( -0.000 -20.119 -1.036) 20.146 4.383 ( 0.000 -4.462 -2.535) 5.132 4.609 ( 0.000 18.253 -10.869) 21.244 5.558 ( 0.000 -1.550 -30.938) 30.976 5.625 ( -0.000 21.348 24.042) 32.152 6.782 ( 0.000 -0.967 26.809) 26.826 7.688 ( 0.000 -11.690 19.807) 23.000 7.989 ( 0.000 21.362 -1.314) 21.402 9.526 ( 0.000 38.406 -3.466) 38.562 11.574 ( 0.000 -15.721 -3.356) 16.075 12.856 ( 0.000 -14.729 -5.845) 15.846 13.436 ( -0.000 -2.507 -20.038) 20.194 14.756 ( 0.000 1.727 14.799) 14.899 ======================= Grid point 225 (44/84) ======================= q-point: ( 0.33 0.33 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.455 ( 4.589 4.589 2.334) 6.898 3.657 ( 6.874 6.874 3.023) 10.181 4.235 ( -9.974 -9.974 -1.118) 14.150 4.355 ( -1.533 -1.533 -5.431) 5.847 4.899 ( 6.779 6.779 -8.141) 12.577 5.312 ( 3.383 3.383 -23.682) 24.161 6.108 ( 7.049 7.049 18.975) 21.434 6.408 ( 6.111 6.111 22.505) 24.107 7.747 ( -10.022 -10.022 20.935) 25.282 7.935 ( 18.592 18.592 -0.672) 26.302 10.180 ( 8.421 8.421 -3.144) 12.317 11.192 ( 14.247 14.247 -4.999) 20.760 12.259 ( -6.404 -6.404 -1.026) 9.115 13.468 ( -1.830 -1.830 -19.508) 19.679 15.188 ( -14.122 -14.122 10.261) 22.453 ======================= Grid point 226 (45/84) ======================= q-point: ( 0.42 0.33 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.512 ( 0.944 -11.046 0.911) 11.123 3.756 ( 2.779 1.913 3.379) 4.775 4.116 ( -1.532 -4.523 0.011) 4.776 4.312 ( -2.725 -4.617 -2.342) 5.850 4.989 ( 2.119 16.925 -16.493) 23.727 5.435 ( 7.234 -0.513 -26.954) 27.912 6.109 ( -3.378 17.526 22.052) 28.370 6.625 ( 10.909 -1.760 24.999) 27.333 7.529 ( -9.539 -12.931 22.059) 27.292 8.264 ( 11.591 17.666 -2.222) 21.245 10.196 ( -2.834 15.735 -2.719) 16.217 11.443 ( 5.421 12.146 -3.825) 13.840 12.350 ( 13.832 -12.735 -3.361) 19.099 13.422 ( -2.729 -2.656 -21.399) 21.735 14.868 ( -14.631 -9.664 13.446) 22.096 ======================= Grid point 227 (46/84) ======================= q-point: (-0.50 0.33 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.516 ( 0.000 -18.541 0.206) 18.542 3.782 ( 0.000 0.393 3.174) 3.198 4.120 ( 0.000 -1.179 0.828) 1.440 4.270 ( 0.000 -5.390 -1.024) 5.487 5.008 ( 0.000 19.038 -19.267) 27.086 5.529 ( 0.000 -1.149 -31.485) 31.506 6.059 ( -0.000 19.595 25.972) 32.535 6.758 ( 0.000 -1.190 26.710) 26.736 7.412 ( 0.000 -14.633 21.578) 26.072 8.389 ( 0.000 16.659 -2.475) 16.842 10.157 ( 0.000 17.395 -2.640) 17.594 11.476 ( 0.000 10.186 -2.784) 10.559 12.563 ( 0.000 -12.879 -4.886) 13.775 13.379 ( -0.000 -3.282 -23.121) 23.353 14.690 ( 0.000 -7.524 15.829) 17.526 ======================= Grid point 239 (47/84) ======================= q-point: ( 0.42 0.42 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.274 ( -9.467 -9.467 -0.007) 13.389 3.784 ( 0.843 0.843 3.878) 4.057 4.100 ( 1.411 1.411 0.913) 2.195 4.228 ( -3.464 -3.464 -0.160) 4.901 5.287 ( 9.336 9.336 -24.267) 27.626 5.445 ( 2.632 2.632 -29.592) 29.826 6.386 ( 7.046 7.046 25.592) 27.463 6.633 ( 4.299 4.299 25.400) 26.118 7.258 ( -12.372 -12.372 23.859) 29.587 8.557 ( 10.024 10.024 -2.987) 14.487 10.293 ( -1.515 -1.515 -1.089) 2.403 11.841 ( 15.655 15.655 -5.370) 22.782 12.196 ( 1.963 1.963 -2.995) 4.084 13.352 ( -4.080 -4.080 -24.219) 24.896 14.629 ( -11.443 -11.443 16.409) 23.046 ======================= Grid point 240 (48/84) ======================= q-point: (-0.50 0.42 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.171 ( 0.000 -12.676 -0.237) 12.679 3.793 ( 0.000 0.415 3.993) 4.015 4.131 ( 0.000 1.090 1.200) 1.621 4.178 ( -0.000 -3.183 0.575) 3.234 5.357 ( -0.000 13.191 -27.371) 30.384 5.511 ( 0.000 -0.640 -31.605) 31.612 6.430 ( -0.000 16.251 27.421) 31.875 6.737 ( 0.000 -0.740 26.374) 26.384 7.090 ( 0.000 -15.945 24.563) 29.285 8.664 ( 0.000 9.360 -3.288) 9.921 10.265 ( 0.000 -1.768 -0.655) 1.885 11.928 ( 0.000 25.988 -4.240) 26.332 12.344 ( 0.000 -7.642 -4.155) 8.699 13.289 ( 0.000 -5.477 -26.481) 27.041 14.489 ( 0.000 -9.981 18.547) 21.062 ======================= Grid point 259 (49/84) ======================= q-point: (-0.50 -0.50 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.033 ( 0.000 0.000 -0.349) 0.349 3.797 ( 0.000 0.000 4.359) 4.359 4.144 ( 0.000 0.000 1.182) 1.182 4.144 ( 0.000 0.000 1.182) 1.182 5.504 ( 0.000 0.000 -31.608) 31.608 5.504 ( 0.000 0.000 -31.608) 31.608 6.729 ( -0.000 -0.000 26.216) 26.216 6.729 ( 0.000 0.000 26.216) 26.216 6.787 ( 0.000 0.000 28.706) 28.706 8.764 ( 0.000 0.000 -3.605) 3.605 10.232 ( 0.000 0.000 -0.093) 0.093 12.262 ( 0.000 0.000 -3.870) 3.870 12.262 ( 0.000 0.000 -3.870) 3.870 13.201 ( 0.000 0.000 -29.292) 29.292 14.367 ( 0.000 0.000 20.429) 20.429 ======================= Grid point 366 (50/84) ======================= q-point: ( 0.17 0.17 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 189 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.543 ( 16.136 16.136 16.136) 27.949 2.543 ( 16.136 16.136 16.136) 27.949 4.064 ( 15.509 15.509 15.509) 26.862 4.801 ( -0.395 -0.395 -0.395) 0.683 4.801 ( -0.395 -0.395 -0.395) 0.683 5.076 ( -7.915 -7.915 -7.915) 13.709 5.653 ( 10.120 10.120 10.120) 17.529 6.222 ( 9.850 9.850 9.850) 17.061 6.222 ( 9.850 9.850 9.850) 17.061 8.903 ( 6.658 6.658 6.658) 11.532 8.903 ( 6.658 6.658 6.658) 11.532 10.377 ( 2.317 2.317 2.317) 4.014 13.226 ( -8.066 -8.066 -8.066) 13.971 13.226 ( -8.066 -8.066 -8.066) 13.971 16.151 ( -5.118 -5.118 -5.118) 8.865 ======================= Grid point 367 (51/84) ======================= q-point: ( 0.25 0.17 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.920 ( 19.641 10.564 10.564) 24.677 2.920 ( 19.250 11.294 11.294) 25.014 4.276 ( 4.819 3.437 3.437) 6.844 4.692 ( -5.682 0.819 0.819) 5.799 4.901 ( 6.701 -4.546 -4.546) 9.286 4.925 ( -6.733 -10.018 -10.018) 15.686 5.801 ( 1.360 19.211 19.211) 27.202 6.499 ( 22.105 8.779 8.779) 25.353 6.544 ( 15.394 9.022 9.022) 19.995 8.746 ( -8.627 12.701 12.701) 19.926 9.162 ( 3.818 9.677 9.677) 14.208 10.581 ( 18.372 -3.611 -3.611) 19.069 12.955 ( -9.941 -8.953 -8.953) 16.098 13.218 ( -0.425 -11.832 -11.832) 16.739 15.944 ( -14.538 -0.670 -0.670) 14.569 ======================= Grid point 368 (52/84) ======================= q-point: ( 0.33 0.17 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 476 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.320 ( 18.830 7.513 7.513) 21.620 3.326 ( 18.335 6.096 6.096) 20.260 4.323 ( 1.207 -2.599 -2.599) 3.869 4.562 ( -7.086 1.122 1.122) 7.262 4.798 ( -5.189 -12.511 -12.511) 18.439 5.045 ( 5.699 -9.409 -9.409) 14.475 5.761 ( -3.614 22.177 22.177) 31.571 6.851 ( 13.797 10.487 10.487) 20.256 7.006 ( 24.036 10.009 10.009) 27.895 8.554 ( -9.467 12.141 12.141) 19.606 9.099 ( -7.899 14.298 14.298) 21.708 11.084 ( 27.039 -6.953 -6.953) 28.772 12.823 ( -2.435 -12.218 -12.218) 17.450 13.209 ( -0.366 -11.748 -11.748) 16.619 15.567 ( -20.922 3.833 3.833) 21.612 ======================= Grid point 369 (53/84) ======================= q-point: ( 0.42 0.17 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.636 ( 10.487 2.869 2.869) 11.244 3.684 ( 15.402 4.595 4.595) 16.717 4.366 ( 3.114 -3.694 -3.694) 6.081 4.392 ( -8.943 2.364 2.364) 9.548 4.717 ( -2.529 -14.642 -14.642) 20.861 5.104 ( 0.017 -14.525 -14.525) 20.541 5.695 ( -2.143 24.175 24.175) 34.256 7.099 ( 9.333 11.543 11.543) 18.804 7.435 ( 16.030 10.815 10.815) 22.156 8.372 ( -7.265 11.553 11.553) 17.881 8.907 ( -8.939 15.656 15.656) 23.878 11.576 ( 17.511 -11.063 -11.063) 23.483 12.861 ( 5.558 -12.257 -12.257) 18.203 13.203 ( -0.211 -11.686 -11.686) 16.528 15.129 ( -19.033 8.167 8.167) 22.263 ======================= Grid point 370 (54/84) ======================= q-point: (-0.50 0.17 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 277 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.745 ( 0.000 1.231 1.231) 1.741 3.886 ( -0.000 2.844 2.844) 4.022 4.258 ( 0.000 3.469 3.469) 4.905 4.413 ( 0.000 -2.894 -2.894) 4.093 4.692 ( 0.000 -15.286 -15.286) 21.618 5.090 ( 0.000 -16.082 -16.082) 22.743 5.675 ( 0.000 24.195 24.195) 34.217 7.204 ( 0.000 12.009 12.009) 16.983 7.616 ( 0.000 10.834 10.834) 15.322 8.289 ( 0.000 11.229 11.229) 15.881 8.798 ( 0.000 16.186 16.186) 22.891 11.758 ( 0.000 -14.059 -14.059) 19.882 12.953 ( 0.000 -11.664 -11.664) 16.495 13.201 ( 0.000 -11.663 -11.663) 16.494 14.900 ( 0.000 10.722 10.722) 15.163 ======================= Grid point 380 (55/84) ======================= q-point: ( 0.25 0.25 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.153 ( 12.767 12.767 8.803) 20.086 3.214 ( 15.252 15.252 5.486) 22.257 4.297 ( -1.026 -1.026 -0.435) 1.515 4.650 ( -10.401 -10.401 -4.393) 15.352 4.721 ( -1.402 -1.402 -4.669) 5.072 4.858 ( -1.027 -1.027 -14.590) 14.662 6.209 ( 12.684 12.684 28.030) 33.278 6.540 ( 8.048 8.048 19.034) 22.177 6.964 ( 22.315 22.315 -2.705) 31.674 8.742 ( -7.151 -7.151 31.151) 32.752 9.527 ( 15.684 15.684 -6.979) 23.253 10.595 ( 8.140 8.140 -5.195) 12.630 12.687 ( -13.136 -13.136 -0.395) 18.581 13.098 ( -2.378 -2.378 -21.797) 22.055 15.837 ( -9.612 -9.612 1.876) 13.723 ======================= Grid point 381 (56/84) ======================= q-point: ( 0.33 0.25 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.434 ( 13.616 5.279 6.115) 15.832 3.527 ( 14.073 10.460 3.043) 17.796 4.268 ( -1.030 -2.276 -3.208) 4.066 4.459 ( -7.206 -14.909 -5.030) 17.306 4.628 ( -6.793 2.051 -7.125) 10.055 4.894 ( 2.645 -5.340 -18.463) 19.401 6.263 ( -2.244 23.464 26.957) 35.809 6.860 ( 17.251 -2.219 22.478) 28.422 7.428 ( 21.217 22.691 -2.350) 31.154 8.556 ( -9.951 -7.037 32.138) 34.371 9.636 ( -3.382 29.890 -5.540) 30.586 10.958 ( 22.877 -3.495 -6.949) 24.164 12.497 ( -3.890 -16.351 -2.537) 16.998 13.064 ( -1.343 -4.033 -21.998) 22.405 15.560 ( -16.278 -4.660 5.061) 17.672 ======================= Grid point 382 (57/84) ======================= q-point: ( 0.42 0.25 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.682 ( 8.765 -0.058 3.438) 9.415 3.776 ( 9.787 3.737 3.320) 10.989 4.270 ( 0.829 -7.972 -3.179) 8.622 4.356 ( -2.702 -12.536 -5.367) 13.901 4.475 ( -7.096 5.269 -7.367) 11.506 4.926 ( 0.351 -4.861 -25.055) 25.525 6.210 ( -2.225 24.022 27.413) 36.517 7.167 ( 11.469 -1.021 22.616) 25.379 7.803 ( 13.794 20.155 -1.422) 24.465 8.357 ( -8.345 -7.776 32.456) 34.402 9.510 ( -6.506 34.808 -4.978) 35.759 11.352 ( 12.250 -8.770 -6.942) 16.587 12.552 ( 8.287 -15.707 -6.706) 18.983 13.034 ( -1.589 -5.724 -20.878) 21.707 15.213 ( -15.141 -0.351 8.463) 17.349 ======================= Grid point 383 (58/84) ======================= q-point: (-0.50 0.25 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.768 ( 0.000 -0.416 2.272) 2.310 3.901 ( -0.000 -3.307 3.541) 4.845 4.276 ( 0.000 -17.140 -2.682) 17.349 4.326 ( 0.000 -2.706 -2.915) 3.978 4.396 ( 0.000 6.944 -8.766) 11.184 4.924 ( 0.000 -3.595 -29.032) 29.254 6.185 ( 0.000 23.925 27.889) 36.745 7.296 ( 0.000 0.038 22.137) 22.137 7.957 ( 0.000 18.634 -1.275) 18.678 8.257 ( 0.000 -8.621 33.081) 34.186 9.431 ( 0.000 36.627 -4.901) 36.953 11.468 ( 0.000 -11.861 -6.716) 13.631 12.687 ( 0.000 -13.356 -10.043) 16.711 13.008 ( -0.000 -7.341 -19.493) 20.829 15.033 ( 0.000 1.788 10.408) 10.560 ======================= Grid point 394 (59/84) ======================= q-point: ( 0.33 0.33 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 476 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.524 ( 2.760 2.760 4.125) 5.679 3.728 ( 8.144 8.144 3.182) 11.948 4.194 ( -8.596 -8.596 -3.193) 12.569 4.248 ( -0.016 -0.016 -4.437) 4.437 4.654 ( -0.971 -0.971 -13.225) 13.296 4.827 ( -0.419 -0.419 -21.923) 21.931 6.646 ( 8.872 8.872 29.750) 32.287 6.868 ( 7.403 7.403 21.301) 23.735 7.869 ( 18.895 18.895 -2.953) 26.884 8.365 ( -10.341 -10.341 34.367) 37.349 10.085 ( 9.496 9.496 -5.940) 14.684 11.050 ( 13.576 13.576 -8.186) 20.872 12.230 ( -7.488 -7.488 -1.690) 10.723 12.997 ( -2.884 -2.884 -24.034) 24.377 15.383 ( -11.710 -11.710 7.037) 17.993 ======================= Grid point 395 (60/84) ======================= q-point: ( 0.42 0.33 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.538 ( -0.908 -13.027 1.487) 13.143 3.844 ( 2.856 3.023 4.525) 6.146 4.113 ( 0.552 -2.199 -0.650) 2.359 4.263 ( 0.353 -1.719 -2.183) 2.801 4.606 ( -2.270 6.965 -18.641) 20.029 4.853 ( 1.699 -2.420 -27.438) 27.596 6.668 ( -1.546 19.673 28.391) 34.575 7.130 ( 12.295 -1.938 23.260) 26.381 8.128 ( -4.494 -5.837 28.194) 29.140 8.212 ( 6.403 10.907 3.669) 13.169 10.114 ( -2.561 18.361 -5.099) 19.227 11.329 ( 8.056 9.824 -6.838) 14.428 12.255 ( 9.548 -12.002 -5.402) 16.260 12.931 ( -3.583 -4.691 -24.036) 24.750 15.121 ( -11.743 -7.940 9.454) 17.039 ======================= Grid point 396 (61/84) ======================= q-point: (-0.50 0.33 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 506 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.522 ( 0.000 -19.655 0.296) 19.657 3.868 ( 0.000 1.646 4.548) 4.837 4.145 ( 0.000 0.222 1.634) 1.649 4.253 ( 0.000 -3.330 -0.640) 3.391 4.584 ( 0.000 9.728 -20.441) 22.638 4.873 ( 0.000 -1.725 -31.092) 31.139 6.649 ( 0.000 20.384 28.550) 35.080 7.275 ( 0.000 -1.510 22.826) 22.876 8.031 ( 0.000 -12.308 36.839) 38.841 8.317 ( 0.000 15.202 -4.230) 15.779 10.077 ( 0.000 20.333 -4.924) 20.920 11.392 ( 0.000 8.374 -5.141) 9.826 12.421 ( 0.000 -11.683 -8.464) 14.426 12.873 ( -0.000 -6.108 -23.802) 24.573 14.981 ( 0.000 -6.240 10.850) 12.516 ======================= Grid point 408 (62/84) ======================= q-point: ( 0.42 0.42 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.274 ( -10.130 -10.130 -0.080) 14.326 3.892 ( 1.393 1.393 5.987) 6.303 4.123 ( 2.086 2.086 0.954) 3.101 4.232 ( -1.679 -1.679 0.659) 2.464 4.733 ( 4.043 4.043 -27.536) 28.123 4.828 ( 0.343 0.343 -29.009) 29.013 6.987 ( 7.842 7.842 29.213) 31.247 7.129 ( 4.780 4.780 22.156) 23.164 7.905 ( -10.624 -10.624 37.758) 40.637 8.470 ( 9.199 9.199 -5.162) 13.996 10.260 ( -0.095 -0.095 -2.153) 2.158 11.681 ( 15.707 15.707 -9.613) 24.203 12.109 ( 0.705 0.705 -5.130) 5.226 12.824 ( -5.836 -5.836 -25.056) 26.380 14.926 ( -9.168 -9.168 11.014) 17.013 ======================= Grid point 409 (63/84) ======================= q-point: (-0.50 0.42 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.163 ( 0.000 -12.987 -0.433) 12.994 3.905 ( 0.000 1.481 6.369) 6.539 4.166 ( 0.000 1.008 1.950) 2.195 4.197 ( -0.000 -1.900 1.051) 2.171 4.764 ( -0.000 6.772 -28.657) 29.446 4.847 ( 0.000 -0.793 -31.660) 31.670 7.026 ( -0.000 15.619 27.311) 31.462 7.246 ( 0.000 -1.069 22.443) 22.469 7.767 ( 0.000 -12.419 40.162) 42.038 8.568 ( 0.000 8.518 -5.712) 10.256 10.244 ( 0.000 0.086 -1.303) 1.306 11.799 ( 0.000 24.913 -8.015) 26.171 12.222 ( 0.000 -6.905 -7.268) 10.025 12.730 ( -0.000 -8.125 -25.821) 27.070 14.816 ( 0.000 -8.089 11.968) 14.445 ======================= Grid point 428 (64/84) ======================= q-point: (-0.50 -0.50 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.023 ( 0.000 0.000 -0.606) 0.606 3.921 ( 0.000 0.000 7.195) 7.195 4.177 ( 0.000 0.000 1.683) 1.683 4.177 ( 0.000 0.000 1.683) 1.683 4.839 ( 0.000 0.000 -31.766) 31.766 4.839 ( 0.000 0.000 -31.766) 31.766 7.235 ( 0.000 0.000 22.243) 22.243 7.235 ( 0.000 0.000 22.243) 22.243 7.590 ( 0.000 0.000 46.109) 46.109 8.659 ( 0.000 0.000 -6.248) 6.248 10.229 ( 0.000 0.000 -0.161) 0.161 12.148 ( 0.000 0.000 -6.807) 6.807 12.148 ( 0.000 0.000 -6.807) 6.807 12.586 ( 0.000 0.000 -28.849) 28.849 14.717 ( 0.000 0.000 12.704) 12.704 ======================= Grid point 549 (65/84) ======================= q-point: ( 0.25 0.25 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 182 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.363 ( 9.977 9.977 9.977) 17.280 3.363 ( 9.977 9.977 9.977) 17.280 4.272 ( -2.378 -2.378 -2.378) 4.119 4.492 ( -10.533 -10.533 -10.533) 18.244 4.633 ( -4.889 -4.889 -4.889) 8.467 4.633 ( -4.889 -4.889 -4.889) 8.467 6.811 ( 20.624 20.624 20.624) 35.721 6.926 ( 12.056 12.056 12.056) 20.882 6.926 ( 12.056 12.056 12.056) 20.882 9.355 ( 7.306 7.306 7.306) 12.654 9.355 ( 7.306 7.306 7.306) 12.654 10.547 ( 3.410 3.410 3.410) 5.907 12.678 ( -9.032 -9.032 -9.032) 15.643 12.678 ( -9.032 -9.032 -9.032) 15.643 15.792 ( -6.344 -6.344 -6.344) 10.988 ======================= Grid point 550 (66/84) ======================= q-point: ( 0.33 0.25 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.572 ( 10.506 5.339 5.339) 12.937 3.601 ( 11.320 5.439 5.439) 13.686 4.217 ( -2.027 -2.063 -2.063) 3.553 4.285 ( -8.325 -11.493 -11.493) 18.262 4.502 ( -6.435 -4.542 -4.542) 9.092 4.592 ( -2.025 -9.776 -9.776) 13.972 6.815 ( -3.164 24.587 24.587) 34.914 7.273 ( 15.180 9.613 9.613) 20.378 7.427 ( 24.148 9.630 9.630) 27.723 9.151 ( -9.987 15.694 15.694) 24.337 9.599 ( 0.318 9.057 9.057) 12.813 10.825 ( 20.855 -4.321 -4.321) 21.732 12.403 ( -7.802 -7.839 -7.839) 13.556 12.669 ( -0.378 -13.501 -13.501) 19.097 15.585 ( -12.832 -2.828 -2.828) 13.441 ======================= Grid point 551 (67/84) ======================= q-point: ( 0.42 0.25 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.760 ( 6.476 2.297 2.297) 7.245 3.813 ( 8.314 -0.123 -0.123) 8.316 4.167 ( -2.434 -10.562 -10.562) 15.133 4.222 ( 3.081 -2.727 -2.727) 4.936 4.358 ( -7.283 -3.819 -3.819) 9.067 4.524 ( -4.670 -11.691 -11.691) 17.180 6.758 ( -2.049 24.445 24.445) 34.631 7.540 ( 9.835 9.507 9.507) 16.658 7.834 ( 14.246 7.905 7.905) 18.109 8.962 ( -7.437 16.014 16.014) 23.837 9.514 ( -5.572 13.026 13.026) 19.246 11.202 ( 11.988 -6.256 -6.256) 14.899 12.377 ( 4.392 -10.574 -10.574) 15.585 12.663 ( -0.215 -13.623 -13.623) 19.267 15.307 ( -12.188 0.366 0.366) 12.199 ======================= Grid point 552 (68/84) ======================= q-point: (-0.50 0.25 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 276 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.824 ( 0.000 1.773 1.773) 2.508 3.904 ( -0.000 -6.496 -6.496) 9.186 4.192 ( 0.000 -4.081 -4.081) 5.771 4.228 ( -0.000 -5.322 -5.322) 7.527 4.297 ( -0.000 -1.610 -1.610) 2.277 4.439 ( 0.000 -14.186 -14.186) 20.062 6.736 ( 0.000 24.452 24.452) 34.580 7.649 ( 0.000 9.285 9.285) 13.132 7.989 ( 0.000 6.840 6.840) 9.673 8.877 ( 0.000 16.290 16.290) 23.038 9.445 ( 0.000 14.266 14.266) 20.175 11.315 ( 0.000 -6.827 -6.827) 9.655 12.457 ( 0.000 -11.574 -11.574) 16.367 12.660 ( 0.000 -13.666 -13.666) 19.327 15.163 ( 0.000 1.918 1.918) 2.712 ======================= Grid point 563 (69/84) ======================= q-point: ( 0.33 0.33 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.620 ( -0.110 -0.110 4.851) 4.853 3.759 ( 7.884 7.884 -0.779) 11.177 4.068 ( -6.859 -6.859 -10.059) 13.974 4.184 ( -0.885 -0.885 -1.477) 1.936 4.419 ( -4.690 -4.690 -8.550) 10.821 4.453 ( -4.072 -4.072 -12.402) 13.674 7.246 ( 7.860 7.860 26.219) 28.478 7.285 ( 7.761 7.761 18.548) 21.552 7.800 ( 18.231 18.231 -3.229) 25.984 9.133 ( -8.157 -8.157 37.168) 38.917 9.939 ( 11.026 11.026 -7.881) 17.471 10.884 ( 12.014 12.014 -6.516) 18.197 12.192 ( -8.875 -8.875 -1.824) 12.682 12.519 ( -3.402 -3.402 -20.947) 21.493 15.451 ( -9.485 -9.485 -0.712) 13.433 ======================= Grid point 564 (70/84) ======================= q-point: ( 0.42 0.33 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.569 ( -3.572 -15.945 1.283) 16.391 3.888 ( 3.719 5.092 -3.549) 7.236 4.044 ( 3.822 -0.837 -8.394) 9.261 4.200 ( 2.787 0.186 -4.748) 5.509 4.302 ( -6.090 -2.192 -7.442) 9.863 4.382 ( -3.763 -3.848 -13.508) 14.541 7.211 ( -1.913 18.925 23.030) 29.869 7.578 ( 11.063 -1.714 19.267) 22.283 8.109 ( 10.716 14.172 -3.940) 18.199 8.959 ( -7.343 -7.830 39.916) 41.334 9.993 ( -1.937 20.871 -6.041) 21.813 11.179 ( 10.665 6.808 -6.683) 14.309 12.136 ( 3.857 -10.881 -5.840) 12.938 12.466 ( -2.259 -6.552 -19.931) 21.102 15.235 ( -9.711 -6.733 1.355) 11.895 ======================= Grid point 565 (71/84) ======================= q-point: (-0.50 0.33 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.525 ( 0.000 -20.870 -0.228) 20.871 3.921 ( 0.000 4.073 -3.282) 5.231 4.148 ( 0.000 -2.538 -15.606) 15.811 4.185 ( -0.000 1.144 1.914) 2.230 4.283 ( -0.000 0.111 -1.872) 1.875 4.294 ( 0.000 -2.759 -20.133) 20.321 7.190 ( 0.000 19.058 22.990) 29.862 7.700 ( 0.000 -1.202 17.699) 17.739 8.224 ( 0.000 12.456 -4.312) 13.182 8.874 ( 0.000 -7.931 41.458) 42.210 9.962 ( 0.000 23.096 -5.571) 23.758 11.276 ( 0.000 6.227 -5.360) 8.216 12.226 ( 0.000 -10.137 -9.840) 14.127 12.430 ( -0.000 -8.999 -18.210) 20.313 15.120 ( 0.000 -5.528 2.379) 6.018 ======================= Grid point 577 (72/84) ======================= q-point: ( 0.42 0.42 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.271 ( -10.902 -10.902 -0.208) 15.419 3.981 ( 3.105 3.105 -2.917) 5.271 4.060 ( 2.118 2.118 -13.507) 13.835 4.215 ( -0.705 -0.705 -16.245) 16.276 4.255 ( -0.629 -0.629 1.284) 1.562 4.332 ( -1.625 -1.625 -11.828) 12.050 7.524 ( 5.365 5.365 21.861) 23.141 7.549 ( 4.669 4.669 18.069) 19.238 8.349 ( 8.264 8.264 -5.932) 13.106 8.785 ( -7.495 -7.495 44.178) 45.432 10.205 ( 2.146 2.146 -3.030) 4.288 11.463 ( 14.669 14.669 -10.456) 23.231 11.992 ( -1.173 -1.173 -5.809) 6.042 12.347 ( -5.490 -5.490 -20.175) 21.618 15.068 ( -7.886 -7.886 2.509) 11.431 ======================= Grid point 578 (73/84) ======================= q-point: (-0.50 0.42 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.153 ( 0.000 -13.329 -0.546) 13.340 4.008 ( 0.000 4.146 -3.231) 5.256 4.113 ( 0.000 -0.948 -20.010) 20.032 4.165 ( -0.000 -2.806 -7.044) 7.582 4.289 ( -0.000 1.242 -10.141) 10.217 4.295 ( 0.000 0.801 -5.042) 5.105 7.518 ( 0.000 11.886 20.010) 23.275 7.667 ( 0.000 -1.445 17.834) 17.893 8.437 ( 0.000 7.359 -6.604) 9.888 8.697 ( 0.000 -7.659 46.471) 47.098 10.211 ( 0.000 3.107 -1.843) 3.612 11.605 ( -0.000 21.693 -9.989) 23.882 12.054 ( 0.000 -5.969 -8.551) 10.428 12.255 ( -0.000 -8.799 -19.202) 21.122 14.973 ( 0.000 -7.152 3.126) 7.805 ======================= Grid point 597 (74/84) ======================= q-point: (-0.50 -0.50 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.009 ( 0.000 0.000 -0.703) 0.703 4.081 ( 0.000 0.000 7.797) 7.797 4.104 ( 0.000 0.000 -20.753) 20.753 4.104 ( 0.000 0.000 -20.753) 20.753 4.305 ( 0.000 0.000 -7.060) 7.060 4.305 ( 0.000 0.000 -7.060) 7.060 7.650 ( 0.000 0.000 17.618) 17.618 7.650 ( 0.000 0.000 17.618) 17.618 8.515 ( 0.000 0.000 -7.218) 7.218 8.606 ( 0.000 0.000 49.190) 49.190 10.226 ( 0.000 0.000 -0.186) 0.186 11.990 ( -0.000 -0.000 -8.064) 8.064 11.990 ( -0.000 -0.000 -8.064) 8.064 12.035 ( 0.000 0.000 -23.191) 23.191 14.886 ( 0.000 0.000 3.568) 3.568 ======================= Grid point 732 (75/84) ======================= q-point: ( 0.33 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 189 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.707 ( -0.860 -0.860 -0.860) 1.490 3.707 ( -0.860 -0.860 -0.860) 1.490 3.810 ( -10.788 -10.788 -10.788) 18.685 4.178 ( 0.245 0.245 0.245) 0.425 4.305 ( -3.172 -3.172 -3.172) 5.494 4.305 ( -3.172 -3.172 -3.172) 5.494 7.632 ( 9.825 9.825 9.825) 17.017 7.632 ( 9.825 9.825 9.825) 17.017 7.817 ( 10.628 10.628 10.628) 18.408 9.768 ( 5.780 5.780 5.780) 10.011 9.768 ( 5.780 5.780 5.780) 10.011 10.866 ( 7.150 7.150 7.150) 12.383 12.157 ( -7.261 -7.261 -7.261) 12.577 12.157 ( -7.261 -7.261 -7.261) 12.577 15.361 ( -7.385 -7.385 -7.385) 12.791 ======================= Grid point 733 (76/84) ======================= q-point: ( 0.42 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.550 ( -8.024 -9.063 -9.063) 15.122 3.652 ( -1.752 -8.174 -8.174) 11.691 3.807 ( 2.470 -11.085 -11.085) 15.870 4.176 ( -0.367 1.699 1.699) 2.431 4.257 ( -1.348 0.728 0.728) 1.696 4.285 ( 0.130 -0.532 -0.532) 0.763 7.620 ( -1.727 16.281 16.281) 23.089 7.879 ( 8.629 6.666 6.666) 12.780 8.076 ( 9.798 3.783 3.783) 11.163 9.618 ( -5.672 14.365 14.365) 21.093 9.940 ( 0.768 6.719 6.719) 9.534 11.115 ( 11.298 2.693 2.693) 11.923 12.009 ( -1.069 -6.936 -6.936) 9.867 12.133 ( -0.811 -11.024 -11.024) 15.611 15.185 ( -8.111 -5.642 -5.642) 11.377 ======================= Grid point 734 (77/84) ======================= q-point: (-0.50 0.33 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 277 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.464 ( 0.000 -10.650 -10.650) 15.061 3.635 ( 0.000 -11.531 -11.531) 16.307 3.834 ( 0.000 -12.825 -12.825) 18.137 4.171 ( 0.000 1.983 1.983) 2.805 4.243 ( 0.000 2.666 2.666) 3.770 4.289 ( 0.000 0.921 0.921) 1.302 7.601 ( 0.000 16.392 16.392) 23.182 7.972 ( 0.000 6.199 6.199) 8.766 8.181 ( 0.000 2.547 2.547) 3.601 9.555 ( 0.000 15.113 15.113) 21.373 9.934 ( 0.000 8.627 8.627) 12.201 11.227 ( 0.000 2.806 2.806) 3.969 12.029 ( 0.000 -8.662 -8.662) 12.251 12.124 ( 0.000 -11.328 -11.328) 16.020 15.088 ( 0.000 -4.881 -4.881) 6.902 ======================= Grid point 746 (78/84) ======================= q-point: ( 0.42 0.42 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.265 ( -11.419 -11.419 -0.428) 16.155 3.626 ( -2.848 -2.848 -21.357) 21.734 3.629 ( -2.746 -2.746 -19.767) 20.145 4.206 ( 0.892 0.892 4.630) 4.799 4.288 ( 1.120 1.120 1.675) 2.305 4.289 ( 0.970 0.970 1.481) 2.019 7.876 ( 4.120 4.120 13.166) 14.398 7.899 ( 3.255 3.255 14.205) 14.932 8.231 ( 6.941 6.941 -4.961) 10.998 9.626 ( -3.832 -3.832 33.385) 33.822 10.138 ( 4.828 4.828 -3.278) 7.574 11.305 ( 11.774 11.774 -2.961) 16.912 11.878 ( -3.246 -3.246 -4.904) 6.717 11.997 ( -3.946 -3.946 -13.369) 14.487 15.039 ( -6.973 -6.973 -4.702) 10.926 ======================= Grid point 747 (79/84) ======================= q-point: (-0.50 0.42 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.141 ( 0.000 -13.397 -0.608) 13.411 3.566 ( 0.000 -0.526 -24.124) 24.130 3.631 ( 0.000 -6.587 -21.502) 22.488 4.214 ( 0.000 1.601 5.392) 5.625 4.294 ( 0.000 1.976 2.297) 3.031 4.308 ( 0.000 0.744 1.882) 2.024 7.868 ( 0.000 9.018 13.672) 16.379 7.983 ( 0.000 -1.587 12.447) 12.548 8.306 ( 0.000 5.973 -5.608) 8.193 9.580 ( 0.000 -3.430 35.334) 35.500 10.171 ( 0.000 6.535 -1.832) 6.788 11.436 ( -0.000 15.067 -4.348) 15.682 11.882 ( 0.000 -5.103 -7.603) 9.156 11.939 ( -0.000 -7.102 -11.338) 13.378 14.955 ( 0.000 -6.533 -4.283) 7.811 ======================= Grid point 766 (80/84) ======================= q-point: (-0.50 -0.50 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.995 ( 0.000 0.000 -0.612) 0.612 3.562 ( 0.000 0.000 -24.846) 24.846 3.562 ( 0.000 0.000 -24.846) 24.846 4.231 ( -0.000 -0.000 6.348) 6.348 4.316 ( 0.000 0.000 2.205) 2.205 4.316 ( 0.000 0.000 2.205) 2.205 7.964 ( 0.000 0.000 12.372) 12.372 7.964 ( 0.000 0.000 12.372) 12.372 8.372 ( 0.000 0.000 -6.250) 6.250 9.537 ( 0.000 0.000 37.668) 37.668 10.222 ( 0.000 0.000 -0.161) 0.161 11.664 ( -0.000 -0.000 -11.279) 11.279 11.828 ( 0.000 0.000 -7.213) 7.213 11.828 ( 0.000 0.000 -7.213) 7.213 14.876 ( 0.000 0.000 -4.006) 4.006 ======================= Grid point 915 (81/84) ======================= q-point: ( 0.42 0.42 0.42) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 182 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.223 ( -7.300 -7.300 -7.300) 12.645 3.296 ( -8.517 -8.517 -8.517) 14.752 3.296 ( -8.517 -8.517 -8.517) 14.752 4.283 ( 2.354 2.354 2.354) 4.077 4.322 ( 1.467 1.467 1.467) 2.541 4.322 ( 1.467 1.467 1.467) 2.541 8.089 ( 4.649 4.649 4.649) 8.052 8.089 ( 4.649 4.649 4.649) 8.052 8.180 ( 2.190 2.190 2.190) 3.793 10.077 ( 4.019 4.019 4.019) 6.962 10.077 ( 4.019 4.019 4.019) 6.962 11.377 ( 7.873 7.873 7.873) 13.637 11.797 ( -4.179 -4.179 -4.179) 7.239 11.797 ( -4.179 -4.179 -4.179) 7.239 14.910 ( -6.373 -6.373 -6.373) 11.038 ======================= Grid point 916 (82/84) ======================= q-point: (-0.50 0.42 0.42) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 276 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.109 ( 0.000 -6.028 -6.028) 8.526 3.175 ( 0.000 -8.799 -8.799) 12.444 3.271 ( 0.000 -12.332 -12.332) 17.440 4.305 ( -0.000 3.011 3.011) 4.258 4.337 ( 0.000 1.658 1.658) 2.345 4.343 ( 0.000 1.285 1.285) 1.818 8.085 ( 0.000 7.157 7.157) 10.121 8.165 ( 0.000 3.064 3.064) 4.334 8.229 ( 0.000 0.194 0.194) 0.275 10.052 ( 0.000 8.128 8.128) 11.495 10.162 ( 0.000 2.924 2.924) 4.136 11.478 ( -0.000 7.266 7.266) 10.276 11.753 ( 0.000 -4.503 -4.503) 6.368 11.767 ( 0.000 -5.512 -5.512) 7.794 14.833 ( 0.000 -6.158 -6.158) 8.709 ======================= Grid point 935 (83/84) ======================= q-point: (-0.50 -0.50 0.42) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.985 ( 0.000 0.000 -0.354) 0.354 3.140 ( 0.000 0.000 -14.750) 14.750 3.140 ( 0.000 0.000 -14.750) 14.750 4.335 ( -0.000 -0.000 3.504) 3.504 4.356 ( 0.000 0.000 1.481) 1.481 4.356 ( 0.000 0.000 1.481) 1.481 8.161 ( 0.000 0.000 6.448) 6.448 8.161 ( 0.000 0.000 6.448) 6.448 8.267 ( 0.000 0.000 -3.606) 3.606 10.087 ( 0.000 0.000 14.685) 14.685 10.219 ( 0.000 0.000 -0.093) 0.093 11.592 ( -0.000 -0.000 3.183) 3.183 11.705 ( 0.000 0.000 -4.288) 4.288 11.705 ( 0.000 0.000 -4.288) 4.288 14.758 ( 0.000 0.000 -6.040) 6.040 ======================= Grid point 1194 (84/84) ======================= q-point: (-0.50 -0.50 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 7.40e-04 1.64e-04 Number of triplets: 44 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.981 ( 0.000 0.000 0.000) 0.000 2.981 ( 0.000 0.000 0.000) 0.000 2.981 ( 0.000 0.000 0.000) 0.000 4.372 ( 0.000 0.000 0.000) 0.000 4.372 ( 0.000 0.000 0.000) 0.000 4.372 ( 0.000 0.000 0.000) 0.000 8.228 ( 0.000 0.000 0.000) 0.000 8.228 ( 0.000 0.000 0.000) 0.000 8.228 ( 0.000 0.000 0.000) 0.000 10.218 ( 0.000 0.000 0.000) 0.000 10.218 ( 0.000 0.000 0.000) 0.000 11.659 ( 0.000 0.000 0.000) 0.000 11.659 ( 0.000 0.000 0.000) 0.000 11.659 ( 0.000 0.000 0.000) 0.000 14.684 ( 0.000 0.000 0.000) 0.000 =================== 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/25920 10.0 3833.649 3833.649 3833.649 -0.000 -0.000 0.000 3/25920 20.0 771.133 771.133 771.133 -0.000 0.000 0.000 3/25920 30.0 239.088 239.088 239.088 -0.000 0.000 0.000 3/25920 40.0 115.680 115.680 115.680 -0.000 0.000 0.000 3/25920 50.0 72.265 72.265 72.265 -0.000 0.000 0.000 3/25920 60.0 52.334 52.334 52.334 -0.000 0.000 0.000 3/25920 70.0 41.265 41.265 41.265 -0.000 0.000 0.000 3/25920 80.0 34.250 34.250 34.250 -0.000 0.000 0.000 3/25920 90.0 29.389 29.389 29.389 -0.000 0.000 0.000 3/25920 100.0 25.807 25.807 25.807 -0.000 0.000 0.000 3/25920 110.0 23.051 23.051 23.051 -0.000 0.000 0.000 3/25920 120.0 20.858 20.858 20.858 -0.000 0.000 0.000 3/25920 130.0 19.070 19.070 19.070 -0.000 0.000 0.000 3/25920 140.0 17.582 17.582 17.582 -0.000 0.000 0.000 3/25920 150.0 16.322 16.322 16.322 -0.000 0.000 0.000 3/25920 160.0 15.241 15.241 15.241 -0.000 0.000 0.000 3/25920 170.0 14.301 14.301 14.301 -0.000 0.000 0.000 3/25920 180.0 13.476 13.476 13.476 -0.000 0.000 0.000 3/25920 190.0 12.745 12.745 12.745 -0.000 0.000 0.000 3/25920 200.0 12.093 12.093 12.093 -0.000 0.000 0.000 3/25920 210.0 11.507 11.507 11.507 -0.000 0.000 0.000 3/25920 220.0 10.977 10.977 10.977 -0.000 0.000 0.000 3/25920 230.0 10.495 10.495 10.495 -0.000 0.000 0.000 3/25920 240.0 10.055 10.055 10.055 -0.000 0.000 0.000 3/25920 250.0 9.651 9.651 9.651 -0.000 0.000 0.000 3/25920 260.0 9.279 9.279 9.279 -0.000 0.000 0.000 3/25920 270.0 8.936 8.936 8.936 -0.000 0.000 0.000 3/25920 280.0 8.617 8.617 8.617 -0.000 0.000 0.000 3/25920 290.0 8.321 8.321 8.321 -0.000 0.000 0.000 3/25920 300.0 8.044 8.044 8.044 -0.000 0.000 0.000 3/25920 310.0 7.786 7.786 7.786 -0.000 0.000 0.000 3/25920 320.0 7.544 7.544 7.544 -0.000 0.000 0.000 3/25920 330.0 7.317 7.317 7.317 -0.000 0.000 0.000 3/25920 340.0 7.103 7.103 7.103 -0.000 0.000 0.000 3/25920 350.0 6.901 6.901 6.901 -0.000 0.000 0.000 3/25920 360.0 6.711 6.711 6.711 -0.000 0.000 0.000 3/25920 370.0 6.531 6.531 6.531 -0.000 0.000 0.000 3/25920 380.0 6.360 6.360 6.360 -0.000 0.000 0.000 3/25920 390.0 6.199 6.199 6.199 -0.000 0.000 0.000 3/25920 400.0 6.045 6.045 6.045 -0.000 0.000 0.000 3/25920 410.0 5.899 5.899 5.899 -0.000 0.000 0.000 3/25920 420.0 5.759 5.759 5.759 -0.000 0.000 0.000 3/25920 430.0 5.626 5.626 5.626 -0.000 0.000 0.000 3/25920 440.0 5.499 5.499 5.499 -0.000 0.000 0.000 3/25920 450.0 5.378 5.378 5.378 -0.000 0.000 0.000 3/25920 460.0 5.262 5.262 5.262 -0.000 0.000 0.000 3/25920 470.0 5.151 5.151 5.151 -0.000 0.000 0.000 3/25920 480.0 5.045 5.045 5.045 -0.000 0.000 0.000 3/25920 490.0 4.942 4.942 4.942 -0.000 0.000 0.000 3/25920 500.0 4.844 4.844 4.844 -0.000 0.000 0.000 3/25920 510.0 4.750 4.750 4.750 -0.000 0.000 0.000 3/25920 520.0 4.659 4.659 4.659 -0.000 0.000 0.000 3/25920 530.0 4.572 4.572 4.572 -0.000 0.000 0.000 3/25920 540.0 4.488 4.488 4.488 -0.000 0.000 0.000 3/25920 550.0 4.407 4.407 4.407 -0.000 0.000 0.000 3/25920 560.0 4.329 4.329 4.329 -0.000 0.000 0.000 3/25920 570.0 4.253 4.253 4.253 -0.000 0.000 0.000 3/25920 580.0 4.181 4.181 4.181 -0.000 0.000 0.000 3/25920 590.0 4.110 4.110 4.110 -0.000 0.000 0.000 3/25920 600.0 4.042 4.042 4.042 -0.000 0.000 0.000 3/25920 610.0 3.976 3.976 3.976 -0.000 0.000 0.000 3/25920 620.0 3.913 3.913 3.913 -0.000 0.000 0.000 3/25920 630.0 3.851 3.851 3.851 -0.000 0.000 0.000 3/25920 640.0 3.791 3.791 3.791 -0.000 0.000 0.000 3/25920 650.0 3.733 3.733 3.733 -0.000 0.000 0.000 3/25920 660.0 3.677 3.677 3.677 -0.000 0.000 0.000 3/25920 670.0 3.622 3.622 3.622 -0.000 0.000 0.000 3/25920 680.0 3.569 3.569 3.569 -0.000 0.000 0.000 3/25920 690.0 3.518 3.518 3.518 -0.000 0.000 0.000 3/25920 700.0 3.468 3.468 3.468 -0.000 0.000 0.000 3/25920 710.0 3.420 3.420 3.420 -0.000 0.000 0.000 3/25920 720.0 3.372 3.372 3.372 -0.000 0.000 0.000 3/25920 730.0 3.326 3.326 3.326 -0.000 0.000 0.000 3/25920 740.0 3.282 3.282 3.282 -0.000 0.000 0.000 3/25920 750.0 3.238 3.238 3.238 -0.000 0.000 0.000 3/25920 760.0 3.196 3.196 3.196 -0.000 0.000 0.000 3/25920 770.0 3.154 3.154 3.154 -0.000 0.000 0.000 3/25920 780.0 3.114 3.114 3.114 -0.000 0.000 0.000 3/25920 790.0 3.075 3.075 3.075 -0.000 0.000 0.000 3/25920 800.0 3.037 3.037 3.037 -0.000 0.000 0.000 3/25920 810.0 3.000 3.000 3.000 -0.000 0.000 0.000 3/25920 820.0 2.963 2.963 2.963 -0.000 0.000 0.000 3/25920 830.0 2.928 2.928 2.928 -0.000 0.000 0.000 3/25920 840.0 2.893 2.893 2.893 -0.000 0.000 0.000 3/25920 850.0 2.859 2.859 2.859 -0.000 0.000 0.000 3/25920 860.0 2.826 2.826 2.826 -0.000 0.000 0.000 3/25920 870.0 2.794 2.794 2.794 -0.000 0.000 0.000 3/25920 880.0 2.762 2.762 2.762 -0.000 0.000 0.000 3/25920 890.0 2.731 2.731 2.731 -0.000 0.000 0.000 3/25920 900.0 2.701 2.701 2.701 -0.000 0.000 0.000 3/25920 910.0 2.671 2.671 2.671 -0.000 0.000 0.000 3/25920 920.0 2.642 2.642 2.642 -0.000 0.000 0.000 3/25920 930.0 2.614 2.614 2.614 -0.000 0.000 0.000 3/25920 940.0 2.586 2.586 2.586 -0.000 0.000 0.000 3/25920 950.0 2.559 2.559 2.559 -0.000 0.000 0.000 3/25920 960.0 2.533 2.533 2.533 -0.000 0.000 0.000 3/25920 970.0 2.507 2.507 2.507 -0.000 0.000 0.000 3/25920 980.0 2.481 2.481 2.481 -0.000 0.000 0.000 3/25920 990.0 2.456 2.456 2.456 -0.000 0.000 0.000 3/25920 1000.0 2.432 2.432 2.432 -0.000 0.000 0.000 3/25920 Thermal conductivity related properties were written into "kappa-m121212.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 19:34:22]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|