----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-09 04:43:19]------------------------- 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.442207480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.442207480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.442207480000000 Atomic positions (fractional): *1 Mg 0.50000000000000 0.50000000000000 0.00000000000000 24.305 2 Mg 0.50000000000000 0.00000000000000 0.50000000000000 24.305 3 Mg 0.00000000000000 0.50000000000000 0.50000000000000 24.305 *4 Sn 0.00000000000000 0.00000000000000 0.00000000000000 118.710 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.442207480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.442207480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.442207480000000 Atomic positions (fractional): *1 Mg 0.50000000000000 0.50000000000000 0.00000000000000 24.305 > 1 2 Mg 0.50000000000000 0.00000000000000 0.50000000000000 24.305 > 2 3 Mg 0.00000000000000 0.50000000000000 0.50000000000000 24.305 > 3 *4 Sn 0.00000000000000 0.00000000000000 0.00000000000000 118.710 > 4 -------------------------------- super cell -------------------------------- Lattice vectors: a 13.326622440000001 0.000000000000000 0.000000000000000 b 0.000000000000000 13.326622440000001 0.000000000000000 c 0.000000000000000 0.000000000000000 13.326622440000001 Atomic positions (fractional): *1 Mg 0.16666666666667 0.16666666666667 0.00000000000000 24.305 > 1 2 Mg 0.50000000000000 0.16666666666667 0.00000000000000 24.305 > 1 3 Mg 0.83333333333333 0.16666666666667 0.00000000000000 24.305 > 1 4 Mg 0.16666666666667 0.50000000000000 0.00000000000000 24.305 > 1 5 Mg 0.50000000000000 0.50000000000000 0.00000000000000 24.305 > 1 6 Mg 0.83333333333333 0.50000000000000 0.00000000000000 24.305 > 1 7 Mg 0.16666666666667 0.83333333333333 0.00000000000000 24.305 > 1 8 Mg 0.50000000000000 0.83333333333333 0.00000000000000 24.305 > 1 9 Mg 0.83333333333333 0.83333333333333 0.00000000000000 24.305 > 1 10 Mg 0.16666666666667 0.16666666666667 0.33333333333333 24.305 > 1 11 Mg 0.50000000000000 0.16666666666667 0.33333333333333 24.305 > 1 12 Mg 0.83333333333333 0.16666666666667 0.33333333333333 24.305 > 1 13 Mg 0.16666666666667 0.50000000000000 0.33333333333333 24.305 > 1 14 Mg 0.50000000000000 0.50000000000000 0.33333333333333 24.305 > 1 15 Mg 0.83333333333333 0.50000000000000 0.33333333333333 24.305 > 1 16 Mg 0.16666666666667 0.83333333333333 0.33333333333333 24.305 > 1 17 Mg 0.50000000000000 0.83333333333333 0.33333333333333 24.305 > 1 18 Mg 0.83333333333333 0.83333333333333 0.33333333333333 24.305 > 1 19 Mg 0.16666666666667 0.16666666666667 0.66666666666667 24.305 > 1 20 Mg 0.50000000000000 0.16666666666667 0.66666666666667 24.305 > 1 21 Mg 0.83333333333333 0.16666666666667 0.66666666666667 24.305 > 1 22 Mg 0.16666666666667 0.50000000000000 0.66666666666667 24.305 > 1 23 Mg 0.50000000000000 0.50000000000000 0.66666666666667 24.305 > 1 24 Mg 0.83333333333333 0.50000000000000 0.66666666666667 24.305 > 1 25 Mg 0.16666666666667 0.83333333333333 0.66666666666667 24.305 > 1 26 Mg 0.50000000000000 0.83333333333333 0.66666666666667 24.305 > 1 27 Mg 0.83333333333333 0.83333333333333 0.66666666666667 24.305 > 1 28 Mg 0.16666666666667 0.00000000000000 0.16666666666667 24.305 > 2 29 Mg 0.50000000000000 0.00000000000000 0.16666666666667 24.305 > 2 30 Mg 0.83333333333333 0.00000000000000 0.16666666666667 24.305 > 2 31 Mg 0.16666666666667 0.33333333333333 0.16666666666667 24.305 > 2 32 Mg 0.50000000000000 0.33333333333333 0.16666666666667 24.305 > 2 33 Mg 0.83333333333333 0.33333333333333 0.16666666666667 24.305 > 2 34 Mg 0.16666666666667 0.66666666666667 0.16666666666667 24.305 > 2 35 Mg 0.50000000000000 0.66666666666667 0.16666666666667 24.305 > 2 36 Mg 0.83333333333333 0.66666666666667 0.16666666666667 24.305 > 2 37 Mg 0.16666666666667 0.00000000000000 0.50000000000000 24.305 > 2 38 Mg 0.50000000000000 0.00000000000000 0.50000000000000 24.305 > 2 39 Mg 0.83333333333333 0.00000000000000 0.50000000000000 24.305 > 2 40 Mg 0.16666666666667 0.33333333333333 0.50000000000000 24.305 > 2 41 Mg 0.50000000000000 0.33333333333333 0.50000000000000 24.305 > 2 42 Mg 0.83333333333333 0.33333333333333 0.50000000000000 24.305 > 2 43 Mg 0.16666666666667 0.66666666666667 0.50000000000000 24.305 > 2 44 Mg 0.50000000000000 0.66666666666667 0.50000000000000 24.305 > 2 45 Mg 0.83333333333333 0.66666666666667 0.50000000000000 24.305 > 2 46 Mg 0.16666666666667 0.00000000000000 0.83333333333333 24.305 > 2 47 Mg 0.50000000000000 0.00000000000000 0.83333333333333 24.305 > 2 48 Mg 0.83333333333333 0.00000000000000 0.83333333333333 24.305 > 2 49 Mg 0.16666666666667 0.33333333333333 0.83333333333333 24.305 > 2 50 Mg 0.50000000000000 0.33333333333333 0.83333333333333 24.305 > 2 51 Mg 0.83333333333333 0.33333333333333 0.83333333333333 24.305 > 2 52 Mg 0.16666666666667 0.66666666666667 0.83333333333333 24.305 > 2 53 Mg 0.50000000000000 0.66666666666667 0.83333333333333 24.305 > 2 54 Mg 0.83333333333333 0.66666666666667 0.83333333333333 24.305 > 2 55 Mg 0.00000000000000 0.16666666666667 0.16666666666667 24.305 > 3 56 Mg 0.33333333333333 0.16666666666667 0.16666666666667 24.305 > 3 57 Mg 0.66666666666667 0.16666666666667 0.16666666666667 24.305 > 3 58 Mg 0.00000000000000 0.50000000000000 0.16666666666667 24.305 > 3 59 Mg 0.33333333333333 0.50000000000000 0.16666666666667 24.305 > 3 60 Mg 0.66666666666667 0.50000000000000 0.16666666666667 24.305 > 3 61 Mg 0.00000000000000 0.83333333333333 0.16666666666667 24.305 > 3 62 Mg 0.33333333333333 0.83333333333333 0.16666666666667 24.305 > 3 63 Mg 0.66666666666667 0.83333333333333 0.16666666666667 24.305 > 3 64 Mg 0.00000000000000 0.16666666666667 0.50000000000000 24.305 > 3 65 Mg 0.33333333333333 0.16666666666667 0.50000000000000 24.305 > 3 66 Mg 0.66666666666667 0.16666666666667 0.50000000000000 24.305 > 3 67 Mg 0.00000000000000 0.50000000000000 0.50000000000000 24.305 > 3 68 Mg 0.33333333333333 0.50000000000000 0.50000000000000 24.305 > 3 69 Mg 0.66666666666667 0.50000000000000 0.50000000000000 24.305 > 3 70 Mg 0.00000000000000 0.83333333333333 0.50000000000000 24.305 > 3 71 Mg 0.33333333333333 0.83333333333333 0.50000000000000 24.305 > 3 72 Mg 0.66666666666667 0.83333333333333 0.50000000000000 24.305 > 3 73 Mg 0.00000000000000 0.16666666666667 0.83333333333333 24.305 > 3 74 Mg 0.33333333333333 0.16666666666667 0.83333333333333 24.305 > 3 75 Mg 0.66666666666667 0.16666666666667 0.83333333333333 24.305 > 3 76 Mg 0.00000000000000 0.50000000000000 0.83333333333333 24.305 > 3 77 Mg 0.33333333333333 0.50000000000000 0.83333333333333 24.305 > 3 78 Mg 0.66666666666667 0.50000000000000 0.83333333333333 24.305 > 3 79 Mg 0.00000000000000 0.83333333333333 0.83333333333333 24.305 > 3 80 Mg 0.33333333333333 0.83333333333333 0.83333333333333 24.305 > 3 81 Mg 0.66666666666667 0.83333333333333 0.83333333333333 24.305 > 3 *82 Sn 0.00000000000000 0.00000000000000 0.00000000000000 118.710 > 4 83 Sn 0.33333333333333 0.00000000000000 0.00000000000000 118.710 > 4 84 Sn 0.66666666666667 0.00000000000000 0.00000000000000 118.710 > 4 85 Sn 0.00000000000000 0.33333333333333 0.00000000000000 118.710 > 4 86 Sn 0.33333333333333 0.33333333333333 0.00000000000000 118.710 > 4 87 Sn 0.66666666666667 0.33333333333333 0.00000000000000 118.710 > 4 88 Sn 0.00000000000000 0.66666666666667 0.00000000000000 118.710 > 4 89 Sn 0.33333333333333 0.66666666666667 0.00000000000000 118.710 > 4 90 Sn 0.66666666666667 0.66666666666667 0.00000000000000 118.710 > 4 91 Sn 0.00000000000000 0.00000000000000 0.33333333333333 118.710 > 4 92 Sn 0.33333333333333 0.00000000000000 0.33333333333333 118.710 > 4 93 Sn 0.66666666666667 0.00000000000000 0.33333333333333 118.710 > 4 94 Sn 0.00000000000000 0.33333333333333 0.33333333333333 118.710 > 4 95 Sn 0.33333333333333 0.33333333333333 0.33333333333333 118.710 > 4 96 Sn 0.66666666666667 0.33333333333333 0.33333333333333 118.710 > 4 97 Sn 0.00000000000000 0.66666666666667 0.33333333333333 118.710 > 4 98 Sn 0.33333333333333 0.66666666666667 0.33333333333333 118.710 > 4 99 Sn 0.66666666666667 0.66666666666667 0.33333333333333 118.710 > 4 100 Sn 0.00000000000000 0.00000000000000 0.66666666666667 118.710 > 4 101 Sn 0.33333333333333 0.00000000000000 0.66666666666667 118.710 > 4 102 Sn 0.66666666666667 0.00000000000000 0.66666666666667 118.710 > 4 103 Sn 0.00000000000000 0.33333333333333 0.66666666666667 118.710 > 4 104 Sn 0.33333333333333 0.33333333333333 0.66666666666667 118.710 > 4 105 Sn 0.66666666666667 0.33333333333333 0.66666666666667 118.710 > 4 106 Sn 0.00000000000000 0.66666666666667 0.66666666666667 118.710 > 4 107 Sn 0.33333333333333 0.66666666666667 0.66666666666667 118.710 > 4 108 Sn 0.66666666666667 0.66666666666667 0.66666666666667 118.710 > 4 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 51.2236657 0.0000000 0.0000000 0.0000000 51.2236657 0.0000000 0.0000000 0.0000000 51.2236657 -------------------------- Born effective charges -------------------------- 1 Mg 2.1308425 0.0000000 0.0000000 0.0000000 2.1308425 0.0000000 0.0000000 0.0000000 0.9460822 2 Mg 2.1308425 0.0000000 0.0000000 0.0000000 0.9460822 0.0000000 0.0000000 0.0000000 2.1308425 3 Mg 0.9460822 0.0000000 0.0000000 0.0000000 2.1308425 0.0000000 0.0000000 0.0000000 2.1308425 4 Sn -5.2077672 0.0000000 0.0000000 0.0000000 -5.2077672 0.0000000 0.0000000 0.0000000 -5.2077672 ---------------------------------------------------------------------------- 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: 324/324 Permutation basis: 2406/2406 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 68 Number of blocks in projector: 68 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: 48 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 20 Use standard eigh solver. Tree of FC basis block matrices: - (68, 65), data: False |-- (20, 20), data: True |-- (48, 45), data: True ----- Solver_atoms: 1 -- 108 / 108 Time (Solver_compr_matrix_reshape): 0.003 Solver_block: 80 / 80 - Time: 0.021 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.025 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 324/324 Permutation basis: 2406/2406 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 68 Number of blocks in projector: 68 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: 48 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 20 Use standard eigh solver. Tree of FC basis block matrices: - (68, 65), data: False |-- (20, 20), data: True |-- (48, 45), data: True Max drift after symmetrization by symfc projector: 0.00000000 (xx) 0.00000000 (xx) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-09 04:43:23]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:43:23]------------------------- 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.442207480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.442207480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.442207480000000 Atomic positions (fractional): 1 Mg 0.50000000000000 0.50000000000000 0.00000000000000 24.305 2 Mg 0.50000000000000 0.00000000000000 0.50000000000000 24.305 3 Mg 0.00000000000000 0.50000000000000 0.50000000000000 24.305 4 Sn 0.00000000000000 0.00000000000000 0.00000000000000 118.710 -------------------------------- supercell --------------------------------- Lattice vectors: a 13.326622440000001 0.000000000000000 0.000000000000000 b 0.000000000000000 13.326622440000001 0.000000000000000 c 0.000000000000000 0.000000000000000 13.326622440000001 Atomic positions (fractional): 1 Mg 0.16666666666667 0.16666666666667 0.00000000000000 24.305 > 1 2 Mg 0.50000000000000 0.16666666666667 0.00000000000000 24.305 > 1 3 Mg 0.83333333333333 0.16666666666667 0.00000000000000 24.305 > 1 4 Mg 0.16666666666667 0.50000000000000 0.00000000000000 24.305 > 1 5 Mg 0.50000000000000 0.50000000000000 0.00000000000000 24.305 > 1 6 Mg 0.83333333333333 0.50000000000000 0.00000000000000 24.305 > 1 7 Mg 0.16666666666667 0.83333333333333 0.00000000000000 24.305 > 1 8 Mg 0.50000000000000 0.83333333333333 0.00000000000000 24.305 > 1 9 Mg 0.83333333333333 0.83333333333333 0.00000000000000 24.305 > 1 10 Mg 0.16666666666667 0.16666666666667 0.33333333333333 24.305 > 1 11 Mg 0.50000000000000 0.16666666666667 0.33333333333333 24.305 > 1 12 Mg 0.83333333333333 0.16666666666667 0.33333333333333 24.305 > 1 13 Mg 0.16666666666667 0.50000000000000 0.33333333333333 24.305 > 1 14 Mg 0.50000000000000 0.50000000000000 0.33333333333333 24.305 > 1 15 Mg 0.83333333333333 0.50000000000000 0.33333333333333 24.305 > 1 16 Mg 0.16666666666667 0.83333333333333 0.33333333333333 24.305 > 1 17 Mg 0.50000000000000 0.83333333333333 0.33333333333333 24.305 > 1 18 Mg 0.83333333333333 0.83333333333333 0.33333333333333 24.305 > 1 19 Mg 0.16666666666667 0.16666666666667 0.66666666666667 24.305 > 1 20 Mg 0.50000000000000 0.16666666666667 0.66666666666667 24.305 > 1 21 Mg 0.83333333333333 0.16666666666667 0.66666666666667 24.305 > 1 22 Mg 0.16666666666667 0.50000000000000 0.66666666666667 24.305 > 1 23 Mg 0.50000000000000 0.50000000000000 0.66666666666667 24.305 > 1 24 Mg 0.83333333333333 0.50000000000000 0.66666666666667 24.305 > 1 25 Mg 0.16666666666667 0.83333333333333 0.66666666666667 24.305 > 1 26 Mg 0.50000000000000 0.83333333333333 0.66666666666667 24.305 > 1 27 Mg 0.83333333333333 0.83333333333333 0.66666666666667 24.305 > 1 28 Mg 0.16666666666667 0.00000000000000 0.16666666666667 24.305 > 28 29 Mg 0.50000000000000 0.00000000000000 0.16666666666667 24.305 > 28 30 Mg 0.83333333333333 0.00000000000000 0.16666666666667 24.305 > 28 31 Mg 0.16666666666667 0.33333333333333 0.16666666666667 24.305 > 28 32 Mg 0.50000000000000 0.33333333333333 0.16666666666667 24.305 > 28 33 Mg 0.83333333333333 0.33333333333333 0.16666666666667 24.305 > 28 34 Mg 0.16666666666667 0.66666666666667 0.16666666666667 24.305 > 28 35 Mg 0.50000000000000 0.66666666666667 0.16666666666667 24.305 > 28 36 Mg 0.83333333333333 0.66666666666667 0.16666666666667 24.305 > 28 37 Mg 0.16666666666667 0.00000000000000 0.50000000000000 24.305 > 28 38 Mg 0.50000000000000 0.00000000000000 0.50000000000000 24.305 > 28 39 Mg 0.83333333333333 0.00000000000000 0.50000000000000 24.305 > 28 40 Mg 0.16666666666667 0.33333333333333 0.50000000000000 24.305 > 28 41 Mg 0.50000000000000 0.33333333333333 0.50000000000000 24.305 > 28 42 Mg 0.83333333333333 0.33333333333333 0.50000000000000 24.305 > 28 43 Mg 0.16666666666667 0.66666666666667 0.50000000000000 24.305 > 28 44 Mg 0.50000000000000 0.66666666666667 0.50000000000000 24.305 > 28 45 Mg 0.83333333333333 0.66666666666667 0.50000000000000 24.305 > 28 46 Mg 0.16666666666667 0.00000000000000 0.83333333333333 24.305 > 28 47 Mg 0.50000000000000 0.00000000000000 0.83333333333333 24.305 > 28 48 Mg 0.83333333333333 0.00000000000000 0.83333333333333 24.305 > 28 49 Mg 0.16666666666667 0.33333333333333 0.83333333333333 24.305 > 28 50 Mg 0.50000000000000 0.33333333333333 0.83333333333333 24.305 > 28 51 Mg 0.83333333333333 0.33333333333333 0.83333333333333 24.305 > 28 52 Mg 0.16666666666667 0.66666666666667 0.83333333333333 24.305 > 28 53 Mg 0.50000000000000 0.66666666666667 0.83333333333333 24.305 > 28 54 Mg 0.83333333333333 0.66666666666667 0.83333333333333 24.305 > 28 55 Mg 0.00000000000000 0.16666666666667 0.16666666666667 24.305 > 55 56 Mg 0.33333333333333 0.16666666666667 0.16666666666667 24.305 > 55 57 Mg 0.66666666666667 0.16666666666667 0.16666666666667 24.305 > 55 58 Mg 0.00000000000000 0.50000000000000 0.16666666666667 24.305 > 55 59 Mg 0.33333333333333 0.50000000000000 0.16666666666667 24.305 > 55 60 Mg 0.66666666666667 0.50000000000000 0.16666666666667 24.305 > 55 61 Mg 0.00000000000000 0.83333333333333 0.16666666666667 24.305 > 55 62 Mg 0.33333333333333 0.83333333333333 0.16666666666667 24.305 > 55 63 Mg 0.66666666666667 0.83333333333333 0.16666666666667 24.305 > 55 64 Mg 0.00000000000000 0.16666666666667 0.50000000000000 24.305 > 55 65 Mg 0.33333333333333 0.16666666666667 0.50000000000000 24.305 > 55 66 Mg 0.66666666666667 0.16666666666667 0.50000000000000 24.305 > 55 67 Mg 0.00000000000000 0.50000000000000 0.50000000000000 24.305 > 55 68 Mg 0.33333333333333 0.50000000000000 0.50000000000000 24.305 > 55 69 Mg 0.66666666666667 0.50000000000000 0.50000000000000 24.305 > 55 70 Mg 0.00000000000000 0.83333333333333 0.50000000000000 24.305 > 55 71 Mg 0.33333333333333 0.83333333333333 0.50000000000000 24.305 > 55 72 Mg 0.66666666666667 0.83333333333333 0.50000000000000 24.305 > 55 73 Mg 0.00000000000000 0.16666666666667 0.83333333333333 24.305 > 55 74 Mg 0.33333333333333 0.16666666666667 0.83333333333333 24.305 > 55 75 Mg 0.66666666666667 0.16666666666667 0.83333333333333 24.305 > 55 76 Mg 0.00000000000000 0.50000000000000 0.83333333333333 24.305 > 55 77 Mg 0.33333333333333 0.50000000000000 0.83333333333333 24.305 > 55 78 Mg 0.66666666666667 0.50000000000000 0.83333333333333 24.305 > 55 79 Mg 0.00000000000000 0.83333333333333 0.83333333333333 24.305 > 55 80 Mg 0.33333333333333 0.83333333333333 0.83333333333333 24.305 > 55 81 Mg 0.66666666666667 0.83333333333333 0.83333333333333 24.305 > 55 82 Sn 0.00000000000000 0.00000000000000 0.00000000000000 118.710 > 82 83 Sn 0.33333333333333 0.00000000000000 0.00000000000000 118.710 > 82 84 Sn 0.66666666666667 0.00000000000000 0.00000000000000 118.710 > 82 85 Sn 0.00000000000000 0.33333333333333 0.00000000000000 118.710 > 82 86 Sn 0.33333333333333 0.33333333333333 0.00000000000000 118.710 > 82 87 Sn 0.66666666666667 0.33333333333333 0.00000000000000 118.710 > 82 88 Sn 0.00000000000000 0.66666666666667 0.00000000000000 118.710 > 82 89 Sn 0.33333333333333 0.66666666666667 0.00000000000000 118.710 > 82 90 Sn 0.66666666666667 0.66666666666667 0.00000000000000 118.710 > 82 91 Sn 0.00000000000000 0.00000000000000 0.33333333333333 118.710 > 82 92 Sn 0.33333333333333 0.00000000000000 0.33333333333333 118.710 > 82 93 Sn 0.66666666666667 0.00000000000000 0.33333333333333 118.710 > 82 94 Sn 0.00000000000000 0.33333333333333 0.33333333333333 118.710 > 82 95 Sn 0.33333333333333 0.33333333333333 0.33333333333333 118.710 > 82 96 Sn 0.66666666666667 0.33333333333333 0.33333333333333 118.710 > 82 97 Sn 0.00000000000000 0.66666666666667 0.33333333333333 118.710 > 82 98 Sn 0.33333333333333 0.66666666666667 0.33333333333333 118.710 > 82 99 Sn 0.66666666666667 0.66666666666667 0.33333333333333 118.710 > 82 100 Sn 0.00000000000000 0.00000000000000 0.66666666666667 118.710 > 82 101 Sn 0.33333333333333 0.00000000000000 0.66666666666667 118.710 > 82 102 Sn 0.66666666666667 0.00000000000000 0.66666666666667 118.710 > 82 103 Sn 0.00000000000000 0.33333333333333 0.66666666666667 118.710 > 82 104 Sn 0.33333333333333 0.33333333333333 0.66666666666667 118.710 > 82 105 Sn 0.66666666666667 0.33333333333333 0.66666666666667 118.710 > 82 106 Sn 0.00000000000000 0.66666666666667 0.66666666666667 118.710 > 82 107 Sn 0.33333333333333 0.66666666666667 0.66666666666667 118.710 > 82 108 Sn 0.66666666666667 0.66666666666667 0.66666666666667 118.710 > 82 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 51.2236657 0.0000000 0.0000000 0.0000000 51.2236657 0.0000000 0.0000000 0.0000000 51.2236657 -------------------------- Born effective charges -------------------------- 1 Mg 2.1308425 0.0000000 0.0000000 0.0000000 2.1308425 0.0000000 0.0000000 0.0000000 0.9460822 2 Mg 2.1308425 0.0000000 0.0000000 0.0000000 0.9460822 0.0000000 0.0000000 0.0000000 2.1308425 3 Mg 0.9460822 0.0000000 0.0000000 0.0000000 2.1308425 0.0000000 0.0000000 0.0000000 2.1308425 4 Sn -5.2077672 0.0000000 0.0000000 0.0000000 -5.2077672 0.0000000 0.0000000 0.0000000 -5.2077672 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 82, 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 (yyz) 0.00000000 (yyz) 0.00000000 (yzy) 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-09 04:43:25]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:43:26]------------------------- 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.442207480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.442207480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.442207480000000 Atomic positions (fractional): 1 Mg 0.50000000000000 0.50000000000000 0.00000000000000 24.305 2 Mg 0.50000000000000 0.00000000000000 0.50000000000000 24.305 3 Mg 0.00000000000000 0.50000000000000 0.50000000000000 24.305 4 Sn 0.00000000000000 0.00000000000000 0.00000000000000 118.710 -------------------------------- supercell --------------------------------- Lattice vectors: a 13.326622440000001 0.000000000000000 0.000000000000000 b 0.000000000000000 13.326622440000001 0.000000000000000 c 0.000000000000000 0.000000000000000 13.326622440000001 Atomic positions (fractional): 1 Mg 0.16666666666667 0.16666666666667 0.00000000000000 24.305 > 1 2 Mg 0.50000000000000 0.16666666666667 0.00000000000000 24.305 > 1 3 Mg 0.83333333333333 0.16666666666667 0.00000000000000 24.305 > 1 4 Mg 0.16666666666667 0.50000000000000 0.00000000000000 24.305 > 1 5 Mg 0.50000000000000 0.50000000000000 0.00000000000000 24.305 > 1 6 Mg 0.83333333333333 0.50000000000000 0.00000000000000 24.305 > 1 7 Mg 0.16666666666667 0.83333333333333 0.00000000000000 24.305 > 1 8 Mg 0.50000000000000 0.83333333333333 0.00000000000000 24.305 > 1 9 Mg 0.83333333333333 0.83333333333333 0.00000000000000 24.305 > 1 10 Mg 0.16666666666667 0.16666666666667 0.33333333333333 24.305 > 1 11 Mg 0.50000000000000 0.16666666666667 0.33333333333333 24.305 > 1 12 Mg 0.83333333333333 0.16666666666667 0.33333333333333 24.305 > 1 13 Mg 0.16666666666667 0.50000000000000 0.33333333333333 24.305 > 1 14 Mg 0.50000000000000 0.50000000000000 0.33333333333333 24.305 > 1 15 Mg 0.83333333333333 0.50000000000000 0.33333333333333 24.305 > 1 16 Mg 0.16666666666667 0.83333333333333 0.33333333333333 24.305 > 1 17 Mg 0.50000000000000 0.83333333333333 0.33333333333333 24.305 > 1 18 Mg 0.83333333333333 0.83333333333333 0.33333333333333 24.305 > 1 19 Mg 0.16666666666667 0.16666666666667 0.66666666666667 24.305 > 1 20 Mg 0.50000000000000 0.16666666666667 0.66666666666667 24.305 > 1 21 Mg 0.83333333333333 0.16666666666667 0.66666666666667 24.305 > 1 22 Mg 0.16666666666667 0.50000000000000 0.66666666666667 24.305 > 1 23 Mg 0.50000000000000 0.50000000000000 0.66666666666667 24.305 > 1 24 Mg 0.83333333333333 0.50000000000000 0.66666666666667 24.305 > 1 25 Mg 0.16666666666667 0.83333333333333 0.66666666666667 24.305 > 1 26 Mg 0.50000000000000 0.83333333333333 0.66666666666667 24.305 > 1 27 Mg 0.83333333333333 0.83333333333333 0.66666666666667 24.305 > 1 28 Mg 0.16666666666667 0.00000000000000 0.16666666666667 24.305 > 28 29 Mg 0.50000000000000 0.00000000000000 0.16666666666667 24.305 > 28 30 Mg 0.83333333333333 0.00000000000000 0.16666666666667 24.305 > 28 31 Mg 0.16666666666667 0.33333333333333 0.16666666666667 24.305 > 28 32 Mg 0.50000000000000 0.33333333333333 0.16666666666667 24.305 > 28 33 Mg 0.83333333333333 0.33333333333333 0.16666666666667 24.305 > 28 34 Mg 0.16666666666667 0.66666666666667 0.16666666666667 24.305 > 28 35 Mg 0.50000000000000 0.66666666666667 0.16666666666667 24.305 > 28 36 Mg 0.83333333333333 0.66666666666667 0.16666666666667 24.305 > 28 37 Mg 0.16666666666667 0.00000000000000 0.50000000000000 24.305 > 28 38 Mg 0.50000000000000 0.00000000000000 0.50000000000000 24.305 > 28 39 Mg 0.83333333333333 0.00000000000000 0.50000000000000 24.305 > 28 40 Mg 0.16666666666667 0.33333333333333 0.50000000000000 24.305 > 28 41 Mg 0.50000000000000 0.33333333333333 0.50000000000000 24.305 > 28 42 Mg 0.83333333333333 0.33333333333333 0.50000000000000 24.305 > 28 43 Mg 0.16666666666667 0.66666666666667 0.50000000000000 24.305 > 28 44 Mg 0.50000000000000 0.66666666666667 0.50000000000000 24.305 > 28 45 Mg 0.83333333333333 0.66666666666667 0.50000000000000 24.305 > 28 46 Mg 0.16666666666667 0.00000000000000 0.83333333333333 24.305 > 28 47 Mg 0.50000000000000 0.00000000000000 0.83333333333333 24.305 > 28 48 Mg 0.83333333333333 0.00000000000000 0.83333333333333 24.305 > 28 49 Mg 0.16666666666667 0.33333333333333 0.83333333333333 24.305 > 28 50 Mg 0.50000000000000 0.33333333333333 0.83333333333333 24.305 > 28 51 Mg 0.83333333333333 0.33333333333333 0.83333333333333 24.305 > 28 52 Mg 0.16666666666667 0.66666666666667 0.83333333333333 24.305 > 28 53 Mg 0.50000000000000 0.66666666666667 0.83333333333333 24.305 > 28 54 Mg 0.83333333333333 0.66666666666667 0.83333333333333 24.305 > 28 55 Mg 0.00000000000000 0.16666666666667 0.16666666666667 24.305 > 55 56 Mg 0.33333333333333 0.16666666666667 0.16666666666667 24.305 > 55 57 Mg 0.66666666666667 0.16666666666667 0.16666666666667 24.305 > 55 58 Mg 0.00000000000000 0.50000000000000 0.16666666666667 24.305 > 55 59 Mg 0.33333333333333 0.50000000000000 0.16666666666667 24.305 > 55 60 Mg 0.66666666666667 0.50000000000000 0.16666666666667 24.305 > 55 61 Mg 0.00000000000000 0.83333333333333 0.16666666666667 24.305 > 55 62 Mg 0.33333333333333 0.83333333333333 0.16666666666667 24.305 > 55 63 Mg 0.66666666666667 0.83333333333333 0.16666666666667 24.305 > 55 64 Mg 0.00000000000000 0.16666666666667 0.50000000000000 24.305 > 55 65 Mg 0.33333333333333 0.16666666666667 0.50000000000000 24.305 > 55 66 Mg 0.66666666666667 0.16666666666667 0.50000000000000 24.305 > 55 67 Mg 0.00000000000000 0.50000000000000 0.50000000000000 24.305 > 55 68 Mg 0.33333333333333 0.50000000000000 0.50000000000000 24.305 > 55 69 Mg 0.66666666666667 0.50000000000000 0.50000000000000 24.305 > 55 70 Mg 0.00000000000000 0.83333333333333 0.50000000000000 24.305 > 55 71 Mg 0.33333333333333 0.83333333333333 0.50000000000000 24.305 > 55 72 Mg 0.66666666666667 0.83333333333333 0.50000000000000 24.305 > 55 73 Mg 0.00000000000000 0.16666666666667 0.83333333333333 24.305 > 55 74 Mg 0.33333333333333 0.16666666666667 0.83333333333333 24.305 > 55 75 Mg 0.66666666666667 0.16666666666667 0.83333333333333 24.305 > 55 76 Mg 0.00000000000000 0.50000000000000 0.83333333333333 24.305 > 55 77 Mg 0.33333333333333 0.50000000000000 0.83333333333333 24.305 > 55 78 Mg 0.66666666666667 0.50000000000000 0.83333333333333 24.305 > 55 79 Mg 0.00000000000000 0.83333333333333 0.83333333333333 24.305 > 55 80 Mg 0.33333333333333 0.83333333333333 0.83333333333333 24.305 > 55 81 Mg 0.66666666666667 0.83333333333333 0.83333333333333 24.305 > 55 82 Sn 0.00000000000000 0.00000000000000 0.00000000000000 118.710 > 82 83 Sn 0.33333333333333 0.00000000000000 0.00000000000000 118.710 > 82 84 Sn 0.66666666666667 0.00000000000000 0.00000000000000 118.710 > 82 85 Sn 0.00000000000000 0.33333333333333 0.00000000000000 118.710 > 82 86 Sn 0.33333333333333 0.33333333333333 0.00000000000000 118.710 > 82 87 Sn 0.66666666666667 0.33333333333333 0.00000000000000 118.710 > 82 88 Sn 0.00000000000000 0.66666666666667 0.00000000000000 118.710 > 82 89 Sn 0.33333333333333 0.66666666666667 0.00000000000000 118.710 > 82 90 Sn 0.66666666666667 0.66666666666667 0.00000000000000 118.710 > 82 91 Sn 0.00000000000000 0.00000000000000 0.33333333333333 118.710 > 82 92 Sn 0.33333333333333 0.00000000000000 0.33333333333333 118.710 > 82 93 Sn 0.66666666666667 0.00000000000000 0.33333333333333 118.710 > 82 94 Sn 0.00000000000000 0.33333333333333 0.33333333333333 118.710 > 82 95 Sn 0.33333333333333 0.33333333333333 0.33333333333333 118.710 > 82 96 Sn 0.66666666666667 0.33333333333333 0.33333333333333 118.710 > 82 97 Sn 0.00000000000000 0.66666666666667 0.33333333333333 118.710 > 82 98 Sn 0.33333333333333 0.66666666666667 0.33333333333333 118.710 > 82 99 Sn 0.66666666666667 0.66666666666667 0.33333333333333 118.710 > 82 100 Sn 0.00000000000000 0.00000000000000 0.66666666666667 118.710 > 82 101 Sn 0.33333333333333 0.00000000000000 0.66666666666667 118.710 > 82 102 Sn 0.66666666666667 0.00000000000000 0.66666666666667 118.710 > 82 103 Sn 0.00000000000000 0.33333333333333 0.66666666666667 118.710 > 82 104 Sn 0.33333333333333 0.33333333333333 0.66666666666667 118.710 > 82 105 Sn 0.66666666666667 0.33333333333333 0.66666666666667 118.710 > 82 106 Sn 0.00000000000000 0.66666666666667 0.66666666666667 118.710 > 82 107 Sn 0.33333333333333 0.66666666666667 0.66666666666667 118.710 > 82 108 Sn 0.66666666666667 0.66666666666667 0.66666666666667 118.710 > 82 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 51.2236657 0.0000000 0.0000000 0.0000000 51.2236657 0.0000000 0.0000000 0.0000000 51.2236657 -------------------------- Born effective charges -------------------------- 1 Mg 2.1308425 0.0000000 0.0000000 0.0000000 2.1308425 0.0000000 0.0000000 0.0000000 0.9460822 2 Mg 2.1308425 0.0000000 0.0000000 0.0000000 0.9460822 0.0000000 0.0000000 0.0000000 2.1308425 3 Mg 0.9460822 0.0000000 0.0000000 0.0000000 2.1308425 0.0000000 0.0000000 0.0000000 2.1308425 4 Sn -5.2077672 0.0000000 0.0000000 0.0000000 -5.2077672 0.0000000 0.0000000 0.0000000 -5.2077672 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000000 (yyz) 0.00000000 (yyz) 0.00000000 (yzy) Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 11 11 11 ] fc3-r2q-transformation over three atoms: True --------------------------- Phonon calculations ---------------------------- Use NAC by Gonze et al. (no real space sum in current implementation) PRB 50, 13035(R) (1994), PRB 55, 10355 (1997) G-cutoff distance: 0.93, Number of G-points: 305, Lambda: 0.70 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/56) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 56 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 3.875 ( 0.000 0.000 0.000) 0.000 3.875 ( 0.000 0.000 0.000) 0.000 3.875 ( 0.000 0.000 0.000) 0.000 5.321 ( 0.000 0.000 0.000) 0.000 5.321 ( 0.000 0.000 0.000) 0.000 5.321 ( 0.000 0.000 0.000) 0.000 7.404 ( 0.000 0.000 0.000) 0.000 7.404 ( 0.000 0.000 0.000) 0.000 7.404 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/56) ======================= q-point: ( 0.09 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.551 ( 26.593 0.000 0.000) 26.593 0.551 ( 26.593 0.000 0.000) 26.593 0.912 ( 43.543 0.000 0.000) 43.543 3.886 ( 0.913 0.000 0.000) 0.913 3.886 ( 0.913 0.000 0.000) 0.913 4.446 ( -0.509 0.000 0.000) 0.509 5.323 ( 0.201 0.000 0.000) 0.201 5.330 ( 0.871 0.000 0.000) 0.871 5.330 ( 0.871 0.000 0.000) 0.871 7.369 ( -3.292 0.000 0.000) 3.292 7.369 ( -3.292 0.000 0.000) 3.292 7.455 ( -4.750 0.000 0.000) 4.750 ======================= Grid point 2 (3/56) ======================= q-point: ( 0.18 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.082 ( 25.159 0.000 0.000) 25.159 1.082 ( 25.159 0.000 0.000) 25.159 1.756 ( 38.079 0.000 0.000) 38.079 3.904 ( 0.563 0.000 0.000) 0.563 3.904 ( 0.563 0.000 0.000) 0.563 4.455 ( 2.184 0.000 0.000) 2.184 5.329 ( 0.337 0.000 0.000) 0.337 5.359 ( 2.100 0.000 0.000) 2.100 5.359 ( 2.100 0.000 0.000) 2.100 7.273 ( -5.921 0.000 0.000) 5.921 7.273 ( -5.921 0.000 0.000) 5.921 7.305 ( -10.014 0.000 0.000) 10.014 ======================= Grid point 3 (4/56) ======================= q-point: ( 0.27 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.573 ( 22.506 0.000 0.000) 22.506 1.573 ( 22.506 0.000 0.000) 22.506 2.429 ( 26.638 0.000 0.000) 26.638 3.896 ( -1.565 0.000 0.000) 1.565 3.896 ( -1.565 0.000 0.000) 1.565 4.579 ( 11.023 0.000 0.000) 11.023 5.336 ( 0.365 0.000 0.000) 0.365 5.418 ( 3.673 0.000 0.000) 3.673 5.418 ( 3.673 0.000 0.000) 3.673 7.041 ( -15.828 0.000 0.000) 15.828 7.136 ( -7.257 0.000 0.000) 7.257 7.136 ( -7.257 0.000 0.000) 7.257 ======================= Grid point 4 (5/56) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.986 ( 17.293 0.000 0.000) 17.293 1.986 ( 17.293 0.000 0.000) 17.293 2.823 ( 12.204 0.000 0.000) 12.204 3.839 ( -3.810 0.000 0.000) 3.810 3.839 ( -3.810 0.000 0.000) 3.810 4.925 ( 22.244 0.000 0.000) 22.244 5.343 ( 0.277 0.000 0.000) 0.277 5.504 ( 4.319 0.000 0.000) 4.319 5.504 ( 4.319 0.000 0.000) 4.319 6.659 ( -21.199 0.000 0.000) 21.199 6.991 ( -6.503 0.000 0.000) 6.503 6.991 ( -6.503 0.000 0.000) 6.503 ======================= Grid point 5 (6/56) ======================= q-point: ( 0.45 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.245 ( 7.039 0.000 0.000) 7.039 2.245 ( 7.039 0.000 0.000) 7.039 2.968 ( 2.962 0.000 0.000) 2.962 3.766 ( -2.425 0.000 0.000) 2.425 3.766 ( -2.425 0.000 0.000) 2.425 5.347 ( 0.103 0.000 0.000) 0.103 5.411 ( 21.267 0.000 0.000) 21.267 5.575 ( 2.101 0.000 0.000) 2.101 5.575 ( 2.101 0.000 0.000) 2.101 6.211 ( -19.774 0.000 0.000) 19.774 6.890 ( -2.775 0.000 0.000) 2.775 6.890 ( -2.775 0.000 0.000) 2.775 ======================= Grid point 12 (7/56) ======================= q-point: ( 0.09 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.536 ( 13.969 13.969 0.000) 19.755 0.785 ( 19.170 19.170 0.000) 27.110 1.391 ( 31.834 31.834 0.000) 45.021 3.848 ( -1.296 -1.296 0.000) 1.832 3.899 ( 1.195 1.195 0.000) 1.690 4.485 ( 1.597 1.597 0.000) 2.259 5.250 ( -3.681 -3.681 0.000) 5.205 5.336 ( 0.674 0.674 0.000) 0.953 5.412 ( 4.528 4.528 0.000) 6.404 7.242 ( -7.609 -7.609 0.000) 10.761 7.335 ( -3.331 -3.331 0.000) 4.711 7.499 ( -0.190 -0.190 0.000) 0.269 ======================= Grid point 13 (8/56) ======================= q-point: ( 0.18 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.995 ( 26.068 -1.651 0.000) 26.120 1.225 ( 22.474 12.845 0.000) 25.886 2.039 ( 29.892 21.369 0.000) 36.745 3.844 ( 0.506 -4.202 0.000) 4.232 3.927 ( 1.186 2.102 0.000) 2.414 4.515 ( 2.349 4.683 0.000) 5.239 5.177 ( -3.430 -8.429 0.000) 9.101 5.360 ( 1.767 0.145 0.000) 1.773 5.514 ( 5.526 8.593 0.000) 10.216 7.066 ( -9.780 -11.515 0.000) 15.107 7.237 ( -6.019 -3.472 0.000) 6.948 7.441 ( -5.200 4.798 0.000) 7.076 ======================= Grid point 14 (9/56) ======================= q-point: ( 0.27 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.519 ( 24.329 -2.828 0.000) 24.493 1.673 ( 20.795 9.188 0.000) 22.734 2.560 ( 19.934 10.341 0.000) 22.457 3.854 ( 0.132 -3.733 0.000) 3.735 3.935 ( -0.554 3.755 0.000) 3.795 4.624 ( 9.055 4.311 0.000) 10.029 5.116 ( -2.242 -12.534 0.000) 12.733 5.412 ( 3.308 -0.553 0.000) 3.354 5.637 ( 6.396 12.466 0.000) 14.011 6.840 ( -12.259 -12.279 0.000) 17.351 7.097 ( -7.408 -3.735 0.000) 8.296 7.293 ( -9.101 7.990 0.000) 12.111 ======================= Grid point 15 (10/56) ======================= q-point: ( 0.36 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.967 ( 18.749 -1.315 0.000) 18.795 2.054 ( 15.795 6.160 0.000) 16.954 2.836 ( 7.493 0.808 0.000) 7.537 3.840 ( -1.507 -0.251 0.000) 1.527 3.901 ( -2.570 6.066 0.000) 6.588 4.869 ( 11.721 -8.568 0.000) 14.519 5.121 ( 5.372 -7.516 0.000) 9.238 5.491 ( 4.049 -1.186 0.000) 4.219 5.758 ( 4.391 14.030 0.000) 14.701 6.577 ( -12.726 -6.317 0.000) 14.208 6.949 ( -6.651 -4.062 0.000) 7.793 7.085 ( -10.570 7.519 0.000) 12.972 ======================= Grid point 16 (11/56) ======================= q-point: ( 0.45 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.251 ( 7.901 0.548 0.000) 7.920 2.288 ( 6.268 3.825 0.000) 7.343 2.908 ( 0.806 -5.273 0.000) 5.334 3.807 ( -1.218 3.281 0.000) 3.500 3.851 ( -1.724 8.155 0.000) 8.335 4.969 ( 0.759 -25.119 0.000) 25.130 5.411 ( 18.138 0.988 0.000) 18.165 5.558 ( 2.006 -1.537 0.000) 2.527 5.740 ( -8.010 5.089 0.000) 9.490 6.380 ( -4.859 14.638 0.000) 15.423 6.846 ( -2.839 -4.307 0.000) 5.158 6.904 ( -5.663 1.630 0.000) 5.893 ======================= Grid point 24 (12/56) ======================= q-point: ( 0.18 0.18 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.144 ( 15.374 15.374 0.000) 21.742 1.550 ( 17.714 17.714 0.000) 25.051 2.469 ( 18.702 18.702 0.000) 26.449 3.777 ( -2.006 -2.006 0.000) 2.836 3.983 ( 3.156 3.156 0.000) 4.463 4.625 ( 6.212 6.212 0.000) 8.785 4.997 ( -8.869 -8.869 0.000) 12.542 5.369 ( 0.868 0.868 0.000) 1.227 5.700 ( 9.597 9.597 0.000) 13.573 6.824 ( -12.128 -12.128 0.000) 17.152 7.135 ( -6.348 -6.348 0.000) 8.978 7.484 ( -0.553 -0.553 0.000) 0.783 ======================= Grid point 25 (13/56) ======================= q-point: ( 0.27 0.18 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.540 ( 21.379 5.686 0.000) 22.122 1.910 ( 16.687 12.941 0.000) 21.116 2.755 ( 8.732 7.229 0.000) 11.336 3.766 ( 0.828 -4.380 0.000) 4.457 4.044 ( 2.562 6.653 0.000) 7.129 4.725 ( -2.392 1.365 0.000) 2.754 4.885 ( 4.181 -5.953 0.000) 7.275 5.401 ( 2.318 -0.342 0.000) 2.343 5.901 ( 9.507 13.493 0.000) 16.506 6.579 ( -11.129 -13.105 0.000) 17.193 6.986 ( -7.930 -6.968 0.000) 10.556 7.419 ( -5.792 4.143 0.000) 7.121 ======================= Grid point 26 (14/56) ======================= q-point: ( 0.36 0.18 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.962 ( 18.845 1.443 0.000) 18.900 2.208 ( 11.961 8.136 0.000) 14.466 2.833 ( -0.467 -1.519 0.000) 1.589 3.806 ( 2.878 -3.467 0.000) 4.506 4.081 ( 0.998 11.244 0.000) 11.288 4.576 ( -8.195 -16.499 0.000) 18.422 5.122 ( 14.944 4.614 0.000) 15.640 5.462 ( 3.354 -1.427 0.000) 3.645 5.997 ( -3.867 8.080 0.000) 8.958 6.463 ( 3.441 -3.127 0.000) 4.649 6.826 ( -7.222 -7.823 0.000) 10.647 7.249 ( -10.722 7.844 0.000) 13.284 ======================= Grid point 27 (15/56) ======================= q-point: ( 0.45 0.18 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.267 ( 9.379 1.008 0.000) 9.433 2.380 ( 4.413 4.575 0.000) 6.356 2.778 ( -3.283 -6.916 0.000) 7.656 3.864 ( 2.004 1.447 0.000) 2.472 4.090 ( 0.082 14.899 0.000) 14.899 4.442 ( -3.939 -25.334 0.000) 25.638 5.444 ( 15.826 2.179 0.000) 15.975 5.520 ( 1.797 -1.945 0.000) 2.649 5.778 ( -14.286 -0.167 0.000) 14.287 6.678 ( 14.013 11.475 0.000) 18.112 6.714 ( -3.119 -8.564 0.000) 9.114 6.992 ( -13.869 8.273 0.000) 16.149 ======================= Grid point 36 (16/56) ======================= q-point: ( 0.27 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.742 ( 13.124 13.124 0.000) 18.561 2.166 ( 11.414 11.414 0.000) 16.142 2.829 ( 0.005 0.005 0.000) 0.006 3.696 ( -2.021 -2.021 0.000) 2.859 4.196 ( 7.868 7.868 0.000) 11.127 4.534 ( -13.218 -13.218 0.000) 18.694 5.010 ( 11.084 11.084 0.000) 15.676 5.405 ( 0.900 0.900 0.000) 1.273 6.192 ( 14.025 14.025 0.000) 19.834 6.310 ( -12.231 -12.231 0.000) 17.298 6.820 ( -9.004 -9.004 0.000) 12.734 7.455 ( -0.772 -0.772 0.000) 1.092 ======================= Grid point 37 (17/56) ======================= q-point: ( 0.36 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.034 ( 14.451 5.489 0.000) 15.458 2.359 ( 7.182 6.106 0.000) 9.427 2.770 ( -5.316 -4.567 0.000) 7.008 3.698 ( 2.565 -6.695 0.000) 7.169 4.261 ( -13.133 -13.672 0.000) 18.957 4.348 ( 6.504 14.447 0.000) 15.843 5.250 ( 12.047 7.295 0.000) 14.083 5.440 ( 2.359 -0.556 0.000) 2.424 6.026 ( -13.382 -4.706 0.000) 14.185 6.525 ( 16.081 8.755 0.000) 18.310 6.634 ( -8.613 -10.748 0.000) 13.774 7.385 ( -6.007 5.097 0.000) 7.878 ======================= Grid point 38 (18/56) ======================= q-point: ( 0.45 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.293 ( 9.345 1.538 0.000) 9.471 2.456 ( 2.371 2.656 0.000) 3.560 2.639 ( -6.139 -6.458 0.000) 8.910 3.805 ( 7.786 -11.200 0.000) 13.641 4.006 ( -11.478 -12.630 0.000) 17.066 4.444 ( 2.518 19.292 0.000) 19.455 5.487 ( 1.604 -1.211 0.000) 2.010 5.498 ( 12.102 2.935 0.000) 12.452 5.757 ( -12.697 -1.597 0.000) 12.797 6.497 ( -3.920 -12.577 0.000) 13.174 6.827 ( 11.672 4.703 0.000) 12.583 7.227 ( -7.989 12.599 0.000) 14.918 ======================= Grid point 48 (19/56) ======================= q-point: ( 0.36 0.36 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.174 ( 7.724 7.724 0.000) 10.924 2.451 ( 3.000 3.000 0.000) 4.243 2.655 ( -6.362 -6.362 0.000) 8.997 3.559 ( -6.149 -6.149 0.000) 8.696 4.027 ( -9.108 -9.108 0.000) 12.880 4.645 ( 13.640 13.640 0.000) 19.290 5.405 ( 7.608 7.608 0.000) 10.760 5.443 ( 0.905 0.905 0.000) 1.280 5.883 ( -8.013 -8.013 0.000) 11.332 6.401 ( -11.472 -11.472 0.000) 16.223 6.738 ( 10.498 10.498 0.000) 14.846 7.450 ( 1.260 1.260 0.000) 1.782 ======================= Grid point 49 (20/56) ======================= q-point: ( 0.45 0.36 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.331 ( 6.856 2.138 0.000) 7.181 2.488 ( 0.793 0.626 0.000) 1.010 2.521 ( -6.208 -4.970 0.000) 7.952 3.470 ( -2.247 -15.908 0.000) 16.066 3.896 ( -3.436 -2.794 0.000) 4.429 4.859 ( 6.014 20.659 0.000) 21.517 5.472 ( 1.384 -0.298 0.000) 1.416 5.560 ( 7.577 3.042 0.000) 8.165 5.719 ( -7.905 -2.104 0.000) 8.181 6.205 ( -6.114 -15.685 0.000) 16.835 6.907 ( 4.932 3.244 0.000) 5.903 7.450 ( -0.587 8.616 0.000) 8.636 ======================= Grid point 60 (21/56) ======================= q-point: ( 0.45 0.45 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.379 ( 2.475 2.475 0.000) 3.500 2.438 ( -3.204 -3.204 0.000) 4.531 2.492 ( -0.040 -0.040 0.000) 0.057 3.235 ( -6.212 -6.212 0.000) 8.785 3.854 ( -1.133 -1.133 0.000) 1.602 5.247 ( 14.736 14.736 0.000) 20.839 5.472 ( 0.427 0.427 0.000) 0.603 5.620 ( 2.764 2.764 0.000) 3.909 5.672 ( -2.410 -2.410 0.000) 3.409 5.886 ( -13.326 -13.326 0.000) 18.845 6.955 ( 1.279 1.279 0.000) 1.809 7.570 ( 2.943 2.943 0.000) 4.163 ======================= Grid point 133 (22/56) ======================= q-point: ( 0.09 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.785 ( 13.449 13.449 13.449) 23.295 0.785 ( 13.449 13.449 13.449) 23.295 1.713 ( 24.621 24.621 24.621) 42.645 3.851 ( -1.028 -1.028 -1.028) 1.781 3.851 ( -1.028 -1.028 -1.028) 1.781 4.554 ( 4.167 4.167 4.167) 7.217 5.182 ( -4.623 -4.623 -4.623) 8.007 5.410 ( 2.696 2.696 2.696) 4.670 5.410 ( 2.696 2.696 2.696) 4.670 7.212 ( -5.913 -5.913 -5.913) 10.241 7.212 ( -5.913 -5.913 -5.913) 10.241 7.540 ( 1.086 1.086 1.086) 1.882 ======================= Grid point 134 (23/56) ======================= q-point: ( 0.18 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.101 ( 21.398 6.802 6.802) 23.461 1.241 ( 22.981 7.096 7.096) 25.077 2.234 ( 24.231 15.758 15.758) 32.921 3.844 ( 1.339 -2.641 -2.641) 3.968 3.852 ( -0.217 -2.366 -2.366) 3.353 4.629 ( 3.786 8.812 8.812) 13.024 5.083 ( -4.980 -7.108 -7.108) 11.219 5.425 ( 1.402 2.881 2.881) 4.309 5.533 ( 6.690 3.351 3.351) 8.198 7.009 ( -10.812 -6.462 -6.462) 14.157 7.124 ( -5.545 -6.943 -6.943) 11.277 7.513 ( -3.708 3.857 3.857) 6.596 ======================= Grid point 135 (24/56) ======================= q-point: ( 0.27 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.569 ( 23.009 2.392 2.392) 23.256 1.694 ( 20.772 5.321 5.321) 22.093 2.647 ( 15.205 7.107 7.107) 18.226 3.835 ( -1.528 -2.753 -2.753) 4.183 3.894 ( 3.276 -1.875 -1.875) 4.214 4.738 ( 7.321 9.814 9.814) 15.691 4.984 ( -4.552 -10.072 -10.072) 14.953 5.472 ( 3.300 2.172 2.172) 4.509 5.684 ( 7.964 5.626 5.626) 11.257 6.777 ( -11.436 -6.443 -6.443) 14.622 6.989 ( -7.445 -6.676 -6.676) 12.024 7.391 ( -8.071 6.139 6.139) 11.854 ======================= Grid point 136 (25/56) ======================= q-point: ( 0.36 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.002 ( 18.458 1.441 1.441) 18.570 2.071 ( 15.484 3.660 3.660) 16.326 2.844 ( 4.546 0.441 0.441) 4.589 3.789 ( -2.868 -2.250 -2.250) 4.284 3.960 ( 2.645 2.805 2.805) 4.768 4.904 ( -3.121 -12.920 -12.920) 18.536 4.930 ( 11.001 2.262 2.262) 11.457 5.550 ( 3.600 2.133 2.133) 4.697 5.851 ( 8.036 8.290 8.290) 14.213 6.565 ( -8.476 -1.638 -1.638) 8.787 6.830 ( -7.763 -7.287 -7.287) 12.902 7.196 ( -10.471 7.001 7.001) 14.411 ======================= Grid point 137 (26/56) ======================= q-point: ( 0.45 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.286 ( 8.017 1.848 1.848) 8.432 2.299 ( 6.091 2.320 2.320) 6.918 2.873 ( -0.404 -3.433 -3.433) 4.872 3.735 ( -1.744 -1.348 -1.348) 2.583 3.992 ( 0.654 8.026 8.026) 11.370 4.860 ( -1.077 -14.722 -14.722) 20.848 5.138 ( 7.206 -10.780 -10.780) 16.862 5.585 ( -0.255 2.469 2.469) 3.501 5.986 ( 4.146 11.264 11.264) 16.460 6.462 ( -1.623 8.946 8.946) 12.755 6.699 ( -4.048 -9.140 -9.140) 13.544 7.007 ( -6.215 5.240 5.240) 9.672 ======================= Grid point 145 (27/56) ======================= q-point: ( 0.18 0.18 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.280 ( 13.237 13.237 12.323) 22.412 1.535 ( 17.897 17.897 -0.256) 25.312 2.564 ( 14.619 14.619 7.966) 22.156 3.766 ( -2.237 -2.237 -1.033) 3.329 3.862 ( 1.108 1.108 -8.190) 8.339 4.820 ( 9.434 9.434 15.430) 20.398 4.920 ( -8.499 -8.499 -6.628) 13.726 5.438 ( -0.185 -0.185 5.846) 5.852 5.669 ( 8.698 8.698 -2.997) 12.661 6.812 ( -11.292 -11.292 -1.228) 16.017 7.003 ( -5.764 -5.764 -10.550) 13.332 7.551 ( -0.469 -0.469 4.293) 4.344 ======================= Grid point 146 (28/56) ======================= q-point: ( 0.27 0.18 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.627 ( 19.234 5.289 7.912) 21.460 1.902 ( 17.073 12.727 -0.263) 21.296 2.786 ( 6.606 5.283 2.629) 8.858 3.750 ( 0.229 -4.319 -1.514) 4.582 3.903 ( 3.306 2.386 -9.927) 10.731 4.738 ( -8.911 -11.975 -5.358) 15.860 4.999 ( 6.895 13.443 14.384) 20.860 5.457 ( 2.800 -2.978 5.095) 6.532 5.863 ( 9.792 10.832 -3.858) 15.103 6.589 ( -9.466 -10.869 0.937) 14.444 6.872 ( -6.987 -5.108 -9.078) 12.542 7.490 ( -5.446 3.123 5.145) 8.117 ======================= Grid point 147 (29/56) ======================= q-point: ( 0.36 0.18 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.015 ( 17.669 1.232 4.834) 18.360 2.208 ( 12.279 8.045 0.081) 14.680 2.837 ( -1.142 -1.506 0.236) 1.905 3.759 ( 0.438 -1.562 -4.104) 4.413 4.006 ( 6.528 2.770 -5.550) 9.005 4.563 ( -8.004 -16.577 -3.899) 18.817 5.070 ( 0.213 7.754 -3.041) 8.332 5.561 ( 6.312 -0.468 9.382) 11.318 6.032 ( 5.525 8.166 0.641) 9.881 6.477 ( -0.496 -7.135 2.610) 7.614 6.733 ( -5.475 -1.214 -6.933) 8.917 7.331 ( -9.931 5.940 6.194) 13.125 ======================= Grid point 148 (30/56) ======================= q-point: ( 0.45 0.18 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.303 ( 8.824 0.761 3.319) 9.458 2.386 ( 4.772 4.657 0.505) 6.687 2.775 ( -3.458 -5.814 -0.523) 6.785 3.764 ( 0.079 3.267 -6.619) 7.382 4.142 ( 5.381 6.022 2.629) 8.493 4.423 ( -4.957 -22.777 -2.594) 23.454 5.048 ( -1.241 -1.504 -22.621) 22.705 5.656 ( 2.238 4.132 9.075) 10.219 6.100 ( 2.096 1.647 18.340) 18.533 6.474 ( -2.062 -10.594 -12.615) 16.602 6.748 ( 6.873 14.366 0.800) 15.945 7.108 ( -10.162 5.151 8.159) 14.013 ======================= Grid point 157 (31/56) ======================= q-point: ( 0.27 0.27 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.806 ( 11.782 11.782 6.070) 17.734 2.161 ( 11.762 11.762 -0.497) 16.641 2.838 ( -0.375 -0.375 0.764) 0.930 3.683 ( -1.840 -1.840 -1.185) 2.859 3.989 ( 5.894 5.894 -14.867) 17.044 4.484 ( -12.417 -12.417 -4.803) 18.205 5.275 ( 11.031 11.031 20.322) 25.619 5.373 ( -3.380 -3.380 -2.748) 5.513 6.092 ( 10.215 10.215 -9.373) 17.220 6.375 ( -8.811 -8.811 5.936) 13.802 6.776 ( -3.812 -3.812 -3.463) 6.407 7.509 ( -1.410 -1.410 4.423) 4.852 ======================= Grid point 158 (32/56) ======================= q-point: ( 0.36 0.27 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.076 ( 13.636 4.772 3.959) 14.980 2.362 ( 7.621 6.514 0.222) 10.028 2.778 ( -5.078 -4.157 0.785) 6.610 3.686 ( 2.292 -5.627 -1.296) 6.213 4.110 ( 2.014 6.009 -14.139) 15.495 4.264 ( -5.412 -10.593 -4.762) 12.813 5.197 ( -8.723 3.353 -10.067) 13.736 5.555 ( 11.142 0.674 14.571) 18.355 6.152 ( -3.591 2.933 2.350) 5.198 6.313 ( 0.313 -7.744 -4.435) 8.930 6.775 ( 5.676 5.102 5.572) 9.449 7.429 ( -6.496 3.377 3.995) 8.340 ======================= Grid point 159 (33/56) ======================= q-point: ( 0.45 0.27 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.323 ( 8.995 1.277 2.870) 9.527 2.468 ( 2.744 3.085 0.992) 4.246 2.652 ( -6.036 -5.794 1.219) 8.455 3.781 ( 6.833 -7.135 -2.669) 10.234 3.964 ( -9.813 -12.581 -3.889) 16.422 4.341 ( 5.186 9.854 -8.090) 13.763 5.030 ( -5.865 -0.478 -24.580) 25.275 5.729 ( 4.930 2.070 12.400) 13.504 6.108 ( 0.043 -0.204 18.621) 18.622 6.256 ( -4.011 -9.730 -12.936) 16.676 6.985 ( 12.412 7.597 7.989) 16.601 7.249 ( -10.632 8.910 3.138) 14.222 ======================= Grid point 169 (34/56) ======================= q-point: ( 0.36 0.36 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.204 ( 7.416 7.416 2.936) 10.891 2.464 ( 3.483 3.483 1.131) 5.054 2.674 ( -5.850 -5.850 1.755) 8.457 3.558 ( -5.956 -5.956 -0.130) 8.423 4.006 ( -8.494 -8.494 -2.164) 12.206 4.361 ( 11.726 11.726 -21.003) 26.761 5.171 ( -5.533 -5.533 -18.683) 20.255 5.610 ( 4.335 4.335 9.175) 11.034 6.147 ( -2.635 -2.635 18.024) 18.405 6.190 ( -4.250 -4.250 -9.003) 10.825 6.925 ( 8.698 8.698 10.795) 16.366 7.461 ( -0.198 -0.198 1.510) 1.536 ======================= Grid point 170 (35/56) ======================= q-point: ( 0.45 0.36 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.359 ( 6.837 2.221 2.684) 7.674 2.508 ( 0.876 0.807 1.838) 2.190 2.548 ( -5.820 -4.079 2.566) 7.556 3.469 ( -2.331 -15.747 -0.102) 15.919 3.881 ( -3.286 -1.694 -1.575) 4.018 4.591 ( 8.783 13.648 -20.182) 25.898 5.023 ( -6.882 -0.372 -25.815) 26.719 5.694 ( 2.774 -6.030 -3.736) 7.617 6.102 ( -1.004 -0.465 14.680) 14.721 6.119 ( -2.552 -2.999 11.440) 12.098 7.083 ( 5.264 2.945 13.031) 14.359 7.433 ( -1.795 7.825 -1.448) 8.158 ======================= Grid point 181 (36/56) ======================= q-point: ( 0.45 0.45 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.411 ( 2.686 2.686 3.048) 4.870 2.469 ( -2.865 -2.865 2.988) 5.034 2.524 ( 0.086 0.086 2.998) 3.001 3.234 ( -6.298 -6.298 -0.096) 8.907 3.851 ( -0.948 -0.948 -0.339) 1.383 4.854 ( 10.366 10.366 -26.330) 30.136 5.002 ( -2.199 -2.199 -27.407) 27.583 5.540 ( -6.287 -6.287 -17.570) 19.692 6.094 ( 0.140 0.140 20.839) 20.840 6.101 ( -0.253 -0.253 23.393) 23.395 7.123 ( 0.963 0.963 13.678) 13.746 7.544 ( 2.809 2.809 -2.418) 4.650 ======================= Grid point 266 (37/56) ======================= q-point: ( 0.18 0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.595 ( 12.059 12.059 12.059) 20.886 1.595 ( 12.059 12.059 12.059) 20.886 2.735 ( 7.511 7.511 7.511) 13.009 3.738 ( -2.417 -2.417 -2.417) 4.187 3.738 ( -2.417 -2.417 -2.417) 4.187 4.756 ( -8.968 -8.968 -8.968) 15.533 5.158 ( 16.031 16.031 16.031) 27.767 5.582 ( 1.432 1.432 1.432) 2.480 5.582 ( 1.432 1.432 1.432) 2.480 6.772 ( -6.799 -6.799 -6.799) 11.776 6.772 ( -6.799 -6.799 -6.799) 11.776 7.609 ( 0.914 0.914 0.914) 1.583 ======================= Grid point 267 (38/56) ======================= q-point: ( 0.27 0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.816 ( 13.930 8.624 8.624) 18.514 1.945 ( 16.296 6.359 6.359) 18.613 2.840 ( 2.466 2.127 2.127) 3.889 3.703 ( -1.575 -3.222 -3.222) 4.821 3.743 ( 3.823 -4.353 -4.353) 7.247 4.566 ( -9.374 -10.207 -10.207) 17.211 5.297 ( -7.284 10.941 10.941) 17.101 5.603 ( 10.878 7.799 7.799) 15.491 5.754 ( 8.768 -4.299 -4.299) 10.669 6.599 ( -4.827 -1.539 -1.539) 5.295 6.701 ( -4.264 -4.918 -4.918) 8.158 7.579 ( -3.852 2.982 2.982) 5.711 ======================= Grid point 268 (39/56) ======================= q-point: ( 0.36 0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.123 ( 14.990 4.201 4.201) 16.125 2.237 ( 11.707 4.051 4.051) 13.033 2.834 ( -2.866 -0.786 -0.786) 3.075 3.678 ( -0.796 -2.848 -2.848) 4.106 3.895 ( 10.861 -5.008 -5.008) 12.966 4.386 ( -7.759 -11.903 -11.903) 18.535 5.030 ( -16.466 -1.590 -1.590) 16.618 5.814 ( 7.100 12.989 12.989) 19.693 5.933 ( 8.402 -8.780 -8.780) 14.993 6.601 ( -5.349 0.864 0.864) 5.486 6.615 ( 6.597 3.932 3.932) 8.627 7.453 ( -8.327 5.162 5.162) 11.074 ======================= Grid point 269 (40/56) ======================= q-point: ( 0.45 0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.382 ( 8.737 2.724 2.724) 9.548 2.405 ( 4.325 2.642 2.642) 5.715 2.747 ( -4.385 -2.368 -2.368) 5.518 3.670 ( -0.087 -1.718 -1.718) 2.431 4.158 ( 12.646 -1.722 -1.722) 12.879 4.268 ( -3.228 -13.576 -13.576) 19.469 4.685 ( -14.827 -10.476 -10.476) 20.960 5.891 ( 1.292 12.649 12.649) 17.935 6.075 ( 4.555 -16.090 -16.090) 23.206 6.501 ( -3.472 8.909 8.909) 13.068 6.810 ( 9.271 7.589 7.589) 14.182 7.265 ( -8.311 6.753 6.753) 12.660 ======================= Grid point 278 (41/56) ======================= q-point: ( 0.27 0.27 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.975 ( 8.741 8.741 9.959) 15.874 2.154 ( 11.588 11.588 0.250) 16.390 2.855 ( -0.973 -0.973 0.728) 1.556 3.648 ( -1.540 -1.540 -2.154) 3.063 3.721 ( 1.944 1.944 -10.034) 10.404 4.345 ( -10.991 -10.991 -8.511) 17.721 5.314 ( -11.136 -11.136 -2.634) 15.968 5.600 ( 0.835 0.835 -0.617) 1.333 5.933 ( 13.529 13.529 5.030) 19.784 6.525 ( -2.543 -2.543 8.101) 8.863 6.710 ( 3.483 3.483 -2.594) 5.567 7.600 ( -1.161 -1.161 3.750) 4.094 ======================= Grid point 279 (42/56) ======================= q-point: ( 0.36 0.27 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.188 ( 11.465 3.505 6.674) 13.722 2.366 ( 8.532 7.009 0.272) 11.045 2.798 ( -4.521 -2.775 0.960) 5.390 3.640 ( 0.904 -2.374 -3.268) 4.139 3.828 ( 7.743 0.091 -11.443) 13.817 4.135 ( -8.515 -11.684 -8.106) 16.575 4.972 ( -19.675 -4.870 -10.367) 22.767 5.653 ( 2.797 -13.574 -13.901) 19.630 6.129 ( 4.968 13.558 12.911) 19.370 6.532 ( 2.262 -2.470 14.650) 15.028 6.834 ( 9.107 11.572 0.657) 14.741 7.529 ( -5.684 2.110 5.220) 8.000 ======================= Grid point 280 (43/56) ======================= q-point: ( 0.45 0.27 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.404 ( 8.005 0.927 4.856) 9.409 2.491 ( 3.796 4.131 1.165) 5.730 2.682 ( -5.885 -3.548 1.487) 7.031 3.686 ( 2.993 1.026 -6.023) 6.803 3.880 ( -4.047 -16.913 -3.626) 17.764 4.155 ( 11.622 -1.007 -10.927) 15.984 4.569 ( -17.908 -3.114 -18.694) 26.074 5.694 ( 1.081 -18.542 -19.458) 26.900 6.175 ( 0.601 12.844 14.234) 19.181 6.569 ( 0.936 -0.066 19.523) 19.545 7.052 ( 9.991 13.202 -0.173) 16.557 7.381 ( -7.633 4.607 8.410) 12.257 ======================= Grid point 290 (44/56) ======================= q-point: ( 0.36 0.36 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.291 ( 6.552 6.552 5.465) 10.758 2.490 ( 4.779 4.779 1.218) 6.868 2.722 ( -4.498 -4.498 2.824) 6.960 3.550 ( -5.318 -5.318 -0.735) 7.556 3.930 ( -6.903 -6.903 -5.567) 11.238 3.935 ( 8.248 8.248 -19.158) 22.430 4.802 ( -11.924 -11.924 -16.253) 23.421 5.426 ( -7.703 -7.703 -21.934) 24.490 6.345 ( 6.604 6.604 18.336) 20.578 6.510 ( 0.828 0.828 16.803) 16.843 7.071 ( 11.351 11.351 3.026) 16.335 7.533 ( -1.817 -1.817 5.568) 6.132 ======================= Grid point 291 (45/56) ======================= q-point: ( 0.45 0.36 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.436 ( 6.606 2.356 4.778) 8.486 2.556 ( 1.751 1.931 2.638) 3.708 2.618 ( -5.349 -2.356 3.993) 7.079 3.464 ( -2.429 -15.209 -0.434) 15.408 3.814 ( -3.006 2.585 -5.715) 6.956 4.174 ( 12.930 3.721 -19.116) 23.377 4.515 ( -14.447 -2.565 -22.555) 26.908 5.334 ( -1.875 -15.956 -24.164) 29.017 6.403 ( 0.538 8.984 18.269) 20.366 6.554 ( 1.915 -0.900 19.980) 20.092 7.286 ( 9.275 9.121 2.493) 13.245 7.454 ( -5.777 2.558 7.141) 9.534 ======================= Grid point 302 (46/56) ======================= q-point: ( 0.45 0.45 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.496 ( 3.097 3.097 5.096) 6.719 2.552 ( -1.841 -1.841 4.758) 5.424 2.604 ( 0.284 0.284 4.446) 4.464 3.231 ( -6.453 -6.453 -0.235) 9.129 3.836 ( 0.039 0.039 -1.320) 1.321 4.305 ( 7.818 7.818 -26.353) 28.578 4.446 ( -4.735 -4.735 -26.177) 27.020 5.082 ( -7.296 -7.296 -24.895) 26.949 6.526 ( 2.291 2.291 19.931) 20.192 6.544 ( 0.529 0.529 19.650) 19.664 7.415 ( 3.353 3.353 11.304) 12.258 7.486 ( 0.630 0.630 -0.583) 1.064 ======================= Grid point 399 (47/56) ======================= q-point: ( 0.27 0.27 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.193 ( 7.429 7.429 7.429) 12.867 2.193 ( 7.429 7.429 7.429) 12.867 2.856 ( -0.965 -0.965 -0.965) 1.671 3.601 ( -1.855 -1.855 -1.855) 3.213 3.601 ( -1.855 -1.855 -1.855) 3.213 4.150 ( -10.106 -10.106 -10.106) 17.504 5.271 ( -12.001 -12.001 -12.001) 20.787 5.271 ( -12.001 -12.001 -12.001) 20.787 6.259 ( 17.084 17.084 17.084) 29.590 6.687 ( 4.358 4.358 4.358) 7.549 6.687 ( 4.358 4.358 4.358) 7.549 7.647 ( 0.614 0.614 0.614) 1.063 ======================= Grid point 400 (48/56) ======================= q-point: ( 0.36 0.27 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.324 ( 8.425 5.521 5.521) 11.486 2.395 ( 8.430 3.818 3.818) 10.011 2.807 ( -3.832 -0.297 -0.297) 3.855 3.570 ( -0.480 -2.423 -2.423) 3.460 3.661 ( 7.268 -5.656 -5.656) 10.808 3.946 ( -9.510 -9.679 -9.679) 16.668 4.783 ( -24.334 -8.804 -8.804) 27.334 5.249 ( -0.681 -20.915 -20.915) 29.585 6.472 ( 3.485 16.455 16.455) 23.530 6.816 ( 4.990 6.563 6.563) 10.537 6.831 ( 12.234 5.804 5.804) 14.732 7.621 ( -3.001 3.434 3.434) 5.709 ======================= Grid point 401 (49/56) ======================= q-point: ( 0.45 0.27 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.507 ( 8.247 3.415 3.415) 9.557 2.518 ( 3.201 3.033 3.033) 5.352 2.703 ( -5.912 0.352 0.352) 5.933 3.588 ( 1.806 -2.641 -2.641) 4.149 3.781 ( -5.661 -10.033 -10.033) 15.276 3.891 ( 14.578 -9.678 -9.678) 19.997 4.309 ( -20.851 -7.976 -7.976) 23.707 5.241 ( -0.186 -22.605 -22.605) 31.969 6.495 ( 0.123 14.922 14.922) 21.104 6.886 ( 1.755 8.436 8.436) 12.058 7.064 ( 7.837 4.714 4.714) 10.289 7.542 ( -3.564 6.851 6.851) 10.324 ======================= Grid point 411 (50/56) ======================= q-point: ( 0.36 0.36 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.421 ( 5.362 5.362 6.946) 10.283 2.514 ( 6.026 6.026 1.368) 8.631 2.778 ( -2.616 -2.616 2.282) 4.346 3.524 ( -4.129 -4.129 -1.918) 6.145 3.615 ( 2.974 2.974 -11.561) 12.303 3.772 ( -6.338 -6.338 -9.521) 13.076 4.537 ( -16.690 -16.690 -9.362) 25.391 4.918 ( -9.487 -9.487 -25.937) 29.202 6.729 ( 6.277 6.277 17.588) 19.702 6.819 ( 2.848 2.848 13.114) 13.719 7.076 ( 10.988 10.988 -1.490) 15.610 7.662 ( 0.536 0.536 6.191) 6.237 ======================= Grid point 412 (51/56) ======================= q-point: ( 0.45 0.36 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.546 ( 5.770 2.413 5.718) 8.474 2.609 ( 3.288 3.744 2.427) 5.542 2.700 ( -4.934 -0.405 3.709) 6.186 3.449 ( -2.220 -13.923 -1.171) 14.147 3.636 ( -1.883 3.436 -10.323) 11.041 3.808 ( 9.652 -2.029 -17.815) 20.363 4.169 ( -16.351 -6.881 -10.339) 20.533 4.819 ( -1.784 -17.450 -24.950) 30.499 6.753 ( -0.281 9.920 15.184) 18.139 6.916 ( 3.134 -1.247 15.200) 15.570 7.250 ( 5.041 9.300 -3.252) 11.067 7.653 ( -0.839 3.935 9.585) 10.395 ======================= Grid point 423 (52/56) ======================= q-point: ( 0.45 0.45 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.609 ( 3.192 3.192 5.615) 7.205 2.651 ( 0.034 0.034 4.587) 4.587 2.695 ( -0.046 -0.046 4.231) 4.231 3.224 ( -6.424 -6.424 -0.443) 9.096 3.748 ( 3.232 3.232 -10.788) 11.716 3.801 ( 4.474 4.474 -22.327) 23.206 4.001 ( -9.338 -9.338 -13.374) 18.795 4.572 ( -6.075 -6.075 -24.006) 25.498 6.888 ( 1.850 1.850 15.182) 15.405 6.897 ( 0.936 0.936 14.718) 14.777 7.384 ( 3.463 3.463 -4.042) 6.350 7.705 ( 1.186 1.186 11.552) 11.674 ======================= Grid point 532 (53/56) ======================= q-point: ( 0.36 0.36 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.561 ( 4.930 4.930 4.930) 8.540 2.561 ( 4.930 4.930 4.930) 8.540 2.801 ( -0.344 -0.344 -0.344) 0.595 3.470 ( -3.109 -3.109 -3.109) 5.385 3.470 ( -3.109 -3.109 -3.109) 5.385 3.578 ( -8.193 -8.193 -8.193) 14.191 4.407 ( -13.454 -13.454 -13.454) 23.304 4.407 ( -13.454 -13.454 -13.454) 23.304 7.034 ( 7.416 7.416 7.416) 12.845 7.039 ( 5.439 5.439 5.439) 9.421 7.039 ( 5.439 5.439 5.439) 9.421 7.766 ( 3.767 3.767 3.767) 6.524 ======================= Grid point 533 (54/56) ======================= q-point: ( 0.45 0.36 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.658 ( 2.638 3.597 3.597) 5.730 2.658 ( 5.454 3.752 3.752) 7.609 2.758 ( -3.784 1.798 1.798) 4.559 3.397 ( -2.508 -6.330 -6.330) 9.296 3.471 ( -2.366 -5.895 -5.895) 8.666 3.478 ( 2.221 -9.410 -9.410) 13.491 4.030 ( -11.708 -5.471 -5.471) 14.034 4.360 ( -1.040 -18.420 -18.420) 26.070 7.009 ( -0.766 9.604 9.604) 13.604 7.154 ( 2.833 4.547 4.547) 7.026 7.204 ( 4.262 2.365 2.365) 5.418 7.815 ( 1.131 6.090 6.090) 8.686 ======================= Grid point 544 (55/56) ======================= q-point: ( 0.45 0.45 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.714 ( 2.627 2.627 4.412) 5.768 2.729 ( 1.925 1.925 3.015) 4.062 2.770 ( -0.650 -0.650 2.874) 3.018 3.213 ( -6.009 -6.009 -0.604) 8.520 3.415 ( 0.161 0.161 -14.966) 14.968 3.427 ( -0.919 -0.919 -15.496) 15.551 3.893 ( -6.289 -6.289 -1.538) 9.026 4.135 ( -3.921 -3.921 -17.838) 18.680 7.140 ( 1.240 1.240 9.397) 9.559 7.142 ( 1.129 1.129 9.094) 9.233 7.310 ( 3.183 3.183 -2.964) 5.390 7.899 ( 2.059 2.059 7.295) 7.855 ======================= Grid point 665 (56/56) ======================= q-point: ( 0.45 0.45 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 7.40e-04 7.40e-04 7.40e-04 3.34e-04 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.778 ( 1.730 1.730 1.730) 2.996 2.778 ( 1.730 1.730 1.730) 2.996 2.801 ( 0.104 0.104 0.104) 0.180 3.198 ( -3.562 -3.562 -3.562) 6.169 3.206 ( -3.880 -3.880 -3.880) 6.721 3.206 ( -3.880 -3.880 -3.880) 6.721 3.876 ( -3.524 -3.524 -3.524) 6.104 3.876 ( -3.524 -3.524 -3.524) 6.104 7.267 ( 1.827 1.827 1.827) 3.164 7.267 ( 1.827 1.827 1.827) 3.164 7.271 ( 1.531 1.531 1.531) 2.652 8.000 ( 2.479 2.479 2.479) 4.293 =================== 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/15972 10.0 355.978 355.978 355.978 0.000 0.000 0.000 3/15972 20.0 125.778 125.778 125.778 -0.000 0.000 0.000 3/15972 30.0 50.531 50.531 50.531 -0.000 0.000 0.000 3/15972 40.0 27.383 27.383 27.383 -0.000 0.000 0.000 3/15972 50.0 17.854 17.854 17.854 0.000 0.000 0.000 3/15972 60.0 13.105 13.105 13.105 0.000 0.000 0.000 3/15972 70.0 10.390 10.390 10.390 0.000 0.000 0.000 3/15972 80.0 8.667 8.667 8.667 0.000 0.000 0.000 3/15972 90.0 7.481 7.481 7.481 0.000 0.000 0.000 3/15972 100.0 6.611 6.611 6.611 0.000 0.000 0.000 3/15972 110.0 5.944 5.944 5.944 0.000 0.000 0.000 3/15972 120.0 5.411 5.411 5.411 0.000 0.000 0.000 3/15972 130.0 4.975 4.975 4.975 0.000 0.000 0.000 3/15972 140.0 4.609 4.609 4.609 0.000 0.000 0.000 3/15972 150.0 4.297 4.297 4.297 0.000 0.000 0.000 3/15972 160.0 4.027 4.027 4.027 0.000 0.000 0.000 3/15972 170.0 3.791 3.791 3.791 0.000 0.000 0.000 3/15972 180.0 3.582 3.582 3.582 0.000 0.000 0.000 3/15972 190.0 3.395 3.395 3.395 0.000 0.000 0.000 3/15972 200.0 3.228 3.228 3.228 0.000 0.000 0.000 3/15972 210.0 3.077 3.077 3.077 0.000 0.000 0.000 3/15972 220.0 2.940 2.940 2.940 0.000 0.000 0.000 3/15972 230.0 2.814 2.814 2.814 0.000 0.000 0.000 3/15972 240.0 2.699 2.699 2.699 0.000 0.000 0.000 3/15972 250.0 2.594 2.594 2.594 0.000 0.000 0.000 3/15972 260.0 2.496 2.496 2.496 0.000 0.000 0.000 3/15972 270.0 2.405 2.405 2.405 0.000 0.000 0.000 3/15972 280.0 2.321 2.321 2.321 0.000 0.000 0.000 3/15972 290.0 2.243 2.243 2.243 0.000 0.000 0.000 3/15972 300.0 2.170 2.170 2.170 0.000 0.000 0.000 3/15972 310.0 2.101 2.101 2.101 0.000 0.000 0.000 3/15972 320.0 2.037 2.037 2.037 0.000 0.000 0.000 3/15972 330.0 1.976 1.976 1.976 0.000 0.000 0.000 3/15972 340.0 1.919 1.919 1.919 0.000 0.000 0.000 3/15972 350.0 1.866 1.866 1.866 0.000 0.000 0.000 3/15972 360.0 1.815 1.815 1.815 0.000 0.000 0.000 3/15972 370.0 1.767 1.767 1.767 0.000 0.000 0.000 3/15972 380.0 1.721 1.721 1.721 0.000 0.000 0.000 3/15972 390.0 1.678 1.678 1.678 0.000 0.000 0.000 3/15972 400.0 1.637 1.637 1.637 0.000 0.000 0.000 3/15972 410.0 1.597 1.597 1.597 0.000 0.000 0.000 3/15972 420.0 1.560 1.560 1.560 0.000 0.000 0.000 3/15972 430.0 1.524 1.524 1.524 0.000 0.000 0.000 3/15972 440.0 1.490 1.490 1.490 0.000 0.000 0.000 3/15972 450.0 1.458 1.458 1.458 0.000 0.000 0.000 3/15972 460.0 1.427 1.427 1.427 0.000 0.000 0.000 3/15972 470.0 1.397 1.397 1.397 0.000 0.000 0.000 3/15972 480.0 1.368 1.368 1.368 0.000 0.000 0.000 3/15972 490.0 1.341 1.341 1.341 0.000 0.000 0.000 3/15972 500.0 1.314 1.314 1.314 0.000 0.000 0.000 3/15972 510.0 1.289 1.289 1.289 0.000 0.000 0.000 3/15972 520.0 1.264 1.264 1.264 0.000 0.000 0.000 3/15972 530.0 1.241 1.241 1.241 0.000 0.000 0.000 3/15972 540.0 1.218 1.218 1.218 0.000 0.000 0.000 3/15972 550.0 1.196 1.196 1.196 0.000 0.000 0.000 3/15972 560.0 1.175 1.175 1.175 0.000 0.000 0.000 3/15972 570.0 1.155 1.155 1.155 0.000 0.000 0.000 3/15972 580.0 1.135 1.135 1.135 0.000 0.000 0.000 3/15972 590.0 1.116 1.116 1.116 0.000 0.000 0.000 3/15972 600.0 1.098 1.098 1.098 0.000 0.000 0.000 3/15972 610.0 1.080 1.080 1.080 0.000 0.000 0.000 3/15972 620.0 1.063 1.063 1.063 0.000 0.000 0.000 3/15972 630.0 1.046 1.046 1.046 0.000 0.000 0.000 3/15972 640.0 1.030 1.030 1.030 0.000 0.000 0.000 3/15972 650.0 1.015 1.015 1.015 0.000 0.000 0.000 3/15972 660.0 0.999 0.999 0.999 0.000 0.000 0.000 3/15972 670.0 0.985 0.985 0.985 0.000 0.000 0.000 3/15972 680.0 0.970 0.970 0.970 0.000 0.000 0.000 3/15972 690.0 0.956 0.956 0.956 0.000 0.000 0.000 3/15972 700.0 0.943 0.943 0.943 0.000 0.000 0.000 3/15972 710.0 0.930 0.930 0.930 0.000 0.000 0.000 3/15972 720.0 0.917 0.917 0.917 0.000 0.000 0.000 3/15972 730.0 0.905 0.905 0.905 0.000 0.000 0.000 3/15972 740.0 0.892 0.892 0.892 0.000 0.000 0.000 3/15972 750.0 0.881 0.881 0.881 0.000 0.000 0.000 3/15972 760.0 0.869 0.869 0.869 0.000 0.000 0.000 3/15972 770.0 0.858 0.858 0.858 0.000 0.000 0.000 3/15972 780.0 0.847 0.847 0.847 0.000 0.000 0.000 3/15972 790.0 0.837 0.837 0.837 0.000 0.000 0.000 3/15972 800.0 0.826 0.826 0.826 0.000 0.000 0.000 3/15972 810.0 0.816 0.816 0.816 0.000 0.000 0.000 3/15972 820.0 0.806 0.806 0.806 0.000 0.000 0.000 3/15972 830.0 0.797 0.797 0.797 0.000 0.000 0.000 3/15972 840.0 0.787 0.787 0.787 0.000 0.000 0.000 3/15972 850.0 0.778 0.778 0.778 0.000 0.000 0.000 3/15972 860.0 0.769 0.769 0.769 0.000 0.000 0.000 3/15972 870.0 0.760 0.760 0.760 0.000 0.000 0.000 3/15972 880.0 0.752 0.752 0.752 0.000 0.000 0.000 3/15972 890.0 0.743 0.743 0.743 0.000 0.000 0.000 3/15972 900.0 0.735 0.735 0.735 0.000 0.000 0.000 3/15972 910.0 0.727 0.727 0.727 0.000 0.000 0.000 3/15972 920.0 0.719 0.719 0.719 0.000 0.000 0.000 3/15972 930.0 0.712 0.712 0.712 0.000 0.000 0.000 3/15972 940.0 0.704 0.704 0.704 0.000 0.000 0.000 3/15972 950.0 0.697 0.697 0.697 0.000 0.000 0.000 3/15972 960.0 0.690 0.690 0.690 0.000 0.000 0.000 3/15972 970.0 0.683 0.683 0.683 0.000 0.000 0.000 3/15972 980.0 0.676 0.676 0.676 0.000 0.000 0.000 3/15972 990.0 0.669 0.669 0.669 0.000 0.000 0.000 3/15972 1000.0 0.662 0.662 0.662 0.000 0.000 0.000 3/15972 Thermal conductivity related properties were written into "kappa-m111111.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-09 04:43:30]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|