----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 22:53:03]------------------------- 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.074999130000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.074999130000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.074999130000000 Atomic positions (fractional): *1 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 2 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 3 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 *4 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 *5 Rb 0.50000000000000 0.50000000000000 0.50000000000000 85.468 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.074999130000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.074999130000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.074999130000000 Atomic positions (fractional): *1 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 1 2 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 2 3 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 3 *4 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 > 4 *5 Rb 0.50000000000000 0.50000000000000 0.50000000000000 85.468 > 5 -------------------------------- super cell -------------------------------- Lattice vectors: a 12.224997390000000 0.000000000000000 0.000000000000000 b 0.000000000000000 12.224997390000000 0.000000000000000 c 0.000000000000000 0.000000000000000 12.224997390000000 Atomic positions (fractional): *1 F 0.16666666666667 0.00000000000000 0.00000000000000 18.998 > 1 2 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 1 3 F 0.83333333333333 0.00000000000000 0.00000000000000 18.998 > 1 4 F 0.16666666666667 0.33333333333333 0.00000000000000 18.998 > 1 5 F 0.50000000000000 0.33333333333333 0.00000000000000 18.998 > 1 6 F 0.83333333333333 0.33333333333333 0.00000000000000 18.998 > 1 7 F 0.16666666666667 0.66666666666667 0.00000000000000 18.998 > 1 8 F 0.50000000000000 0.66666666666667 0.00000000000000 18.998 > 1 9 F 0.83333333333333 0.66666666666667 0.00000000000000 18.998 > 1 10 F 0.16666666666667 0.00000000000000 0.33333333333333 18.998 > 1 11 F 0.50000000000000 0.00000000000000 0.33333333333333 18.998 > 1 12 F 0.83333333333333 0.00000000000000 0.33333333333333 18.998 > 1 13 F 0.16666666666667 0.33333333333333 0.33333333333333 18.998 > 1 14 F 0.50000000000000 0.33333333333333 0.33333333333333 18.998 > 1 15 F 0.83333333333333 0.33333333333333 0.33333333333333 18.998 > 1 16 F 0.16666666666667 0.66666666666667 0.33333333333333 18.998 > 1 17 F 0.50000000000000 0.66666666666667 0.33333333333333 18.998 > 1 18 F 0.83333333333333 0.66666666666667 0.33333333333333 18.998 > 1 19 F 0.16666666666667 0.00000000000000 0.66666666666667 18.998 > 1 20 F 0.50000000000000 0.00000000000000 0.66666666666667 18.998 > 1 21 F 0.83333333333333 0.00000000000000 0.66666666666667 18.998 > 1 22 F 0.16666666666667 0.33333333333333 0.66666666666667 18.998 > 1 23 F 0.50000000000000 0.33333333333333 0.66666666666667 18.998 > 1 24 F 0.83333333333333 0.33333333333333 0.66666666666667 18.998 > 1 25 F 0.16666666666667 0.66666666666667 0.66666666666667 18.998 > 1 26 F 0.50000000000000 0.66666666666667 0.66666666666667 18.998 > 1 27 F 0.83333333333333 0.66666666666667 0.66666666666667 18.998 > 1 28 F 0.00000000000000 0.16666666666667 0.00000000000000 18.998 > 2 29 F 0.33333333333333 0.16666666666667 0.00000000000000 18.998 > 2 30 F 0.66666666666667 0.16666666666667 0.00000000000000 18.998 > 2 31 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 2 32 F 0.33333333333333 0.50000000000000 0.00000000000000 18.998 > 2 33 F 0.66666666666667 0.50000000000000 0.00000000000000 18.998 > 2 34 F 0.00000000000000 0.83333333333333 0.00000000000000 18.998 > 2 35 F 0.33333333333333 0.83333333333333 0.00000000000000 18.998 > 2 36 F 0.66666666666667 0.83333333333333 0.00000000000000 18.998 > 2 37 F 0.00000000000000 0.16666666666667 0.33333333333333 18.998 > 2 38 F 0.33333333333333 0.16666666666667 0.33333333333333 18.998 > 2 39 F 0.66666666666667 0.16666666666667 0.33333333333333 18.998 > 2 40 F 0.00000000000000 0.50000000000000 0.33333333333333 18.998 > 2 41 F 0.33333333333333 0.50000000000000 0.33333333333333 18.998 > 2 42 F 0.66666666666667 0.50000000000000 0.33333333333333 18.998 > 2 43 F 0.00000000000000 0.83333333333333 0.33333333333333 18.998 > 2 44 F 0.33333333333333 0.83333333333333 0.33333333333333 18.998 > 2 45 F 0.66666666666667 0.83333333333333 0.33333333333333 18.998 > 2 46 F 0.00000000000000 0.16666666666667 0.66666666666667 18.998 > 2 47 F 0.33333333333333 0.16666666666667 0.66666666666667 18.998 > 2 48 F 0.66666666666667 0.16666666666667 0.66666666666667 18.998 > 2 49 F 0.00000000000000 0.50000000000000 0.66666666666667 18.998 > 2 50 F 0.33333333333333 0.50000000000000 0.66666666666667 18.998 > 2 51 F 0.66666666666667 0.50000000000000 0.66666666666667 18.998 > 2 52 F 0.00000000000000 0.83333333333333 0.66666666666667 18.998 > 2 53 F 0.33333333333333 0.83333333333333 0.66666666666667 18.998 > 2 54 F 0.66666666666667 0.83333333333333 0.66666666666667 18.998 > 2 55 F 0.00000000000000 0.00000000000000 0.16666666666667 18.998 > 3 56 F 0.33333333333333 0.00000000000000 0.16666666666667 18.998 > 3 57 F 0.66666666666667 0.00000000000000 0.16666666666667 18.998 > 3 58 F 0.00000000000000 0.33333333333333 0.16666666666667 18.998 > 3 59 F 0.33333333333333 0.33333333333333 0.16666666666667 18.998 > 3 60 F 0.66666666666667 0.33333333333333 0.16666666666667 18.998 > 3 61 F 0.00000000000000 0.66666666666667 0.16666666666667 18.998 > 3 62 F 0.33333333333333 0.66666666666667 0.16666666666667 18.998 > 3 63 F 0.66666666666667 0.66666666666667 0.16666666666667 18.998 > 3 64 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 3 65 F 0.33333333333333 0.00000000000000 0.50000000000000 18.998 > 3 66 F 0.66666666666667 0.00000000000000 0.50000000000000 18.998 > 3 67 F 0.00000000000000 0.33333333333333 0.50000000000000 18.998 > 3 68 F 0.33333333333333 0.33333333333333 0.50000000000000 18.998 > 3 69 F 0.66666666666667 0.33333333333333 0.50000000000000 18.998 > 3 70 F 0.00000000000000 0.66666666666667 0.50000000000000 18.998 > 3 71 F 0.33333333333333 0.66666666666667 0.50000000000000 18.998 > 3 72 F 0.66666666666667 0.66666666666667 0.50000000000000 18.998 > 3 73 F 0.00000000000000 0.00000000000000 0.83333333333333 18.998 > 3 74 F 0.33333333333333 0.00000000000000 0.83333333333333 18.998 > 3 75 F 0.66666666666667 0.00000000000000 0.83333333333333 18.998 > 3 76 F 0.00000000000000 0.33333333333333 0.83333333333333 18.998 > 3 77 F 0.33333333333333 0.33333333333333 0.83333333333333 18.998 > 3 78 F 0.66666666666667 0.33333333333333 0.83333333333333 18.998 > 3 79 F 0.00000000000000 0.66666666666667 0.83333333333333 18.998 > 3 80 F 0.33333333333333 0.66666666666667 0.83333333333333 18.998 > 3 81 F 0.66666666666667 0.66666666666667 0.83333333333333 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 Rb 0.16666666666667 0.16666666666667 0.16666666666667 85.468 > 5 110 Rb 0.50000000000000 0.16666666666667 0.16666666666667 85.468 > 5 111 Rb 0.83333333333333 0.16666666666667 0.16666666666667 85.468 > 5 112 Rb 0.16666666666667 0.50000000000000 0.16666666666667 85.468 > 5 113 Rb 0.50000000000000 0.50000000000000 0.16666666666667 85.468 > 5 114 Rb 0.83333333333333 0.50000000000000 0.16666666666667 85.468 > 5 115 Rb 0.16666666666667 0.83333333333333 0.16666666666667 85.468 > 5 116 Rb 0.50000000000000 0.83333333333333 0.16666666666667 85.468 > 5 117 Rb 0.83333333333333 0.83333333333333 0.16666666666667 85.468 > 5 118 Rb 0.16666666666667 0.16666666666667 0.50000000000000 85.468 > 5 119 Rb 0.50000000000000 0.16666666666667 0.50000000000000 85.468 > 5 120 Rb 0.83333333333333 0.16666666666667 0.50000000000000 85.468 > 5 121 Rb 0.16666666666667 0.50000000000000 0.50000000000000 85.468 > 5 122 Rb 0.50000000000000 0.50000000000000 0.50000000000000 85.468 > 5 123 Rb 0.83333333333333 0.50000000000000 0.50000000000000 85.468 > 5 124 Rb 0.16666666666667 0.83333333333333 0.50000000000000 85.468 > 5 125 Rb 0.50000000000000 0.83333333333333 0.50000000000000 85.468 > 5 126 Rb 0.83333333333333 0.83333333333333 0.50000000000000 85.468 > 5 127 Rb 0.16666666666667 0.16666666666667 0.83333333333333 85.468 > 5 128 Rb 0.50000000000000 0.16666666666667 0.83333333333333 85.468 > 5 129 Rb 0.83333333333333 0.16666666666667 0.83333333333333 85.468 > 5 130 Rb 0.16666666666667 0.50000000000000 0.83333333333333 85.468 > 5 131 Rb 0.50000000000000 0.50000000000000 0.83333333333333 85.468 > 5 132 Rb 0.83333333333333 0.50000000000000 0.83333333333333 85.468 > 5 133 Rb 0.16666666666667 0.83333333333333 0.83333333333333 85.468 > 5 134 Rb 0.50000000000000 0.83333333333333 0.83333333333333 85.468 > 5 135 Rb 0.83333333333333 0.83333333333333 0.83333333333333 85.468 > 5 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 2.2517795 0.0000000 0.0000000 0.0000000 2.2517795 0.0000000 0.0000000 0.0000000 2.2517795 -------------------------- Born effective charges -------------------------- 1 F -1.2627655 0.0000000 0.0000000 0.0000000 -1.0294943 0.0000000 0.0000000 0.0000000 -1.0294943 2 F -1.0294943 0.0000000 0.0000000 0.0000000 -1.2627655 0.0000000 0.0000000 0.0000000 -1.0294943 3 F -1.0294943 0.0000000 0.0000000 0.0000000 -1.0294943 0.0000000 0.0000000 0.0000000 -1.2627655 4 Mg 2.0696769 0.0000000 0.0000000 0.0000000 2.0696769 0.0000000 0.0000000 0.0000000 2.0696769 5 Rb 1.2520772 0.0000000 0.0000000 0.0000000 1.2520772 0.0000000 0.0000000 0.0000000 1.2520772 ---------------------------------------------------------------------------- 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.099 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.103 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) Permutation basis: 405/405 Permutation basis: 3615/3615 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 104 Number of blocks in projector: 104 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 74 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 30 Use standard eigh solver. Tree of FC basis block matrices: - (104, 100), data: False |-- (30, 30), data: True |-- (74, 70), data: True Max drift after symmetrization by symfc projector: -0.00000000 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 22:53:09]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 22:53: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.074999130000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.074999130000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.074999130000000 Atomic positions (fractional): 1 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 2 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 3 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 4 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 5 Rb 0.50000000000000 0.50000000000000 0.50000000000000 85.468 -------------------------------- supercell --------------------------------- Lattice vectors: a 12.224997390000000 0.000000000000000 0.000000000000000 b 0.000000000000000 12.224997390000000 0.000000000000000 c 0.000000000000000 0.000000000000000 12.224997390000000 Atomic positions (fractional): 1 F 0.16666666666667 0.00000000000000 0.00000000000000 18.998 > 1 2 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 1 3 F 0.83333333333333 0.00000000000000 0.00000000000000 18.998 > 1 4 F 0.16666666666667 0.33333333333333 0.00000000000000 18.998 > 1 5 F 0.50000000000000 0.33333333333333 0.00000000000000 18.998 > 1 6 F 0.83333333333333 0.33333333333333 0.00000000000000 18.998 > 1 7 F 0.16666666666667 0.66666666666667 0.00000000000000 18.998 > 1 8 F 0.50000000000000 0.66666666666667 0.00000000000000 18.998 > 1 9 F 0.83333333333333 0.66666666666667 0.00000000000000 18.998 > 1 10 F 0.16666666666667 0.00000000000000 0.33333333333333 18.998 > 1 11 F 0.50000000000000 0.00000000000000 0.33333333333333 18.998 > 1 12 F 0.83333333333333 0.00000000000000 0.33333333333333 18.998 > 1 13 F 0.16666666666667 0.33333333333333 0.33333333333333 18.998 > 1 14 F 0.50000000000000 0.33333333333333 0.33333333333333 18.998 > 1 15 F 0.83333333333333 0.33333333333333 0.33333333333333 18.998 > 1 16 F 0.16666666666667 0.66666666666667 0.33333333333333 18.998 > 1 17 F 0.50000000000000 0.66666666666667 0.33333333333333 18.998 > 1 18 F 0.83333333333333 0.66666666666667 0.33333333333333 18.998 > 1 19 F 0.16666666666667 0.00000000000000 0.66666666666667 18.998 > 1 20 F 0.50000000000000 0.00000000000000 0.66666666666667 18.998 > 1 21 F 0.83333333333333 0.00000000000000 0.66666666666667 18.998 > 1 22 F 0.16666666666667 0.33333333333333 0.66666666666667 18.998 > 1 23 F 0.50000000000000 0.33333333333333 0.66666666666667 18.998 > 1 24 F 0.83333333333333 0.33333333333333 0.66666666666667 18.998 > 1 25 F 0.16666666666667 0.66666666666667 0.66666666666667 18.998 > 1 26 F 0.50000000000000 0.66666666666667 0.66666666666667 18.998 > 1 27 F 0.83333333333333 0.66666666666667 0.66666666666667 18.998 > 1 28 F 0.00000000000000 0.16666666666667 0.00000000000000 18.998 > 28 29 F 0.33333333333333 0.16666666666667 0.00000000000000 18.998 > 28 30 F 0.66666666666667 0.16666666666667 0.00000000000000 18.998 > 28 31 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 28 32 F 0.33333333333333 0.50000000000000 0.00000000000000 18.998 > 28 33 F 0.66666666666667 0.50000000000000 0.00000000000000 18.998 > 28 34 F 0.00000000000000 0.83333333333333 0.00000000000000 18.998 > 28 35 F 0.33333333333333 0.83333333333333 0.00000000000000 18.998 > 28 36 F 0.66666666666667 0.83333333333333 0.00000000000000 18.998 > 28 37 F 0.00000000000000 0.16666666666667 0.33333333333333 18.998 > 28 38 F 0.33333333333333 0.16666666666667 0.33333333333333 18.998 > 28 39 F 0.66666666666667 0.16666666666667 0.33333333333333 18.998 > 28 40 F 0.00000000000000 0.50000000000000 0.33333333333333 18.998 > 28 41 F 0.33333333333333 0.50000000000000 0.33333333333333 18.998 > 28 42 F 0.66666666666667 0.50000000000000 0.33333333333333 18.998 > 28 43 F 0.00000000000000 0.83333333333333 0.33333333333333 18.998 > 28 44 F 0.33333333333333 0.83333333333333 0.33333333333333 18.998 > 28 45 F 0.66666666666667 0.83333333333333 0.33333333333333 18.998 > 28 46 F 0.00000000000000 0.16666666666667 0.66666666666667 18.998 > 28 47 F 0.33333333333333 0.16666666666667 0.66666666666667 18.998 > 28 48 F 0.66666666666667 0.16666666666667 0.66666666666667 18.998 > 28 49 F 0.00000000000000 0.50000000000000 0.66666666666667 18.998 > 28 50 F 0.33333333333333 0.50000000000000 0.66666666666667 18.998 > 28 51 F 0.66666666666667 0.50000000000000 0.66666666666667 18.998 > 28 52 F 0.00000000000000 0.83333333333333 0.66666666666667 18.998 > 28 53 F 0.33333333333333 0.83333333333333 0.66666666666667 18.998 > 28 54 F 0.66666666666667 0.83333333333333 0.66666666666667 18.998 > 28 55 F 0.00000000000000 0.00000000000000 0.16666666666667 18.998 > 55 56 F 0.33333333333333 0.00000000000000 0.16666666666667 18.998 > 55 57 F 0.66666666666667 0.00000000000000 0.16666666666667 18.998 > 55 58 F 0.00000000000000 0.33333333333333 0.16666666666667 18.998 > 55 59 F 0.33333333333333 0.33333333333333 0.16666666666667 18.998 > 55 60 F 0.66666666666667 0.33333333333333 0.16666666666667 18.998 > 55 61 F 0.00000000000000 0.66666666666667 0.16666666666667 18.998 > 55 62 F 0.33333333333333 0.66666666666667 0.16666666666667 18.998 > 55 63 F 0.66666666666667 0.66666666666667 0.16666666666667 18.998 > 55 64 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 55 65 F 0.33333333333333 0.00000000000000 0.50000000000000 18.998 > 55 66 F 0.66666666666667 0.00000000000000 0.50000000000000 18.998 > 55 67 F 0.00000000000000 0.33333333333333 0.50000000000000 18.998 > 55 68 F 0.33333333333333 0.33333333333333 0.50000000000000 18.998 > 55 69 F 0.66666666666667 0.33333333333333 0.50000000000000 18.998 > 55 70 F 0.00000000000000 0.66666666666667 0.50000000000000 18.998 > 55 71 F 0.33333333333333 0.66666666666667 0.50000000000000 18.998 > 55 72 F 0.66666666666667 0.66666666666667 0.50000000000000 18.998 > 55 73 F 0.00000000000000 0.00000000000000 0.83333333333333 18.998 > 55 74 F 0.33333333333333 0.00000000000000 0.83333333333333 18.998 > 55 75 F 0.66666666666667 0.00000000000000 0.83333333333333 18.998 > 55 76 F 0.00000000000000 0.33333333333333 0.83333333333333 18.998 > 55 77 F 0.33333333333333 0.33333333333333 0.83333333333333 18.998 > 55 78 F 0.66666666666667 0.33333333333333 0.83333333333333 18.998 > 55 79 F 0.00000000000000 0.66666666666667 0.83333333333333 18.998 > 55 80 F 0.33333333333333 0.66666666666667 0.83333333333333 18.998 > 55 81 F 0.66666666666667 0.66666666666667 0.83333333333333 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 Rb 0.16666666666667 0.16666666666667 0.16666666666667 85.468 > 109 110 Rb 0.50000000000000 0.16666666666667 0.16666666666667 85.468 > 109 111 Rb 0.83333333333333 0.16666666666667 0.16666666666667 85.468 > 109 112 Rb 0.16666666666667 0.50000000000000 0.16666666666667 85.468 > 109 113 Rb 0.50000000000000 0.50000000000000 0.16666666666667 85.468 > 109 114 Rb 0.83333333333333 0.50000000000000 0.16666666666667 85.468 > 109 115 Rb 0.16666666666667 0.83333333333333 0.16666666666667 85.468 > 109 116 Rb 0.50000000000000 0.83333333333333 0.16666666666667 85.468 > 109 117 Rb 0.83333333333333 0.83333333333333 0.16666666666667 85.468 > 109 118 Rb 0.16666666666667 0.16666666666667 0.50000000000000 85.468 > 109 119 Rb 0.50000000000000 0.16666666666667 0.50000000000000 85.468 > 109 120 Rb 0.83333333333333 0.16666666666667 0.50000000000000 85.468 > 109 121 Rb 0.16666666666667 0.50000000000000 0.50000000000000 85.468 > 109 122 Rb 0.50000000000000 0.50000000000000 0.50000000000000 85.468 > 109 123 Rb 0.83333333333333 0.50000000000000 0.50000000000000 85.468 > 109 124 Rb 0.16666666666667 0.83333333333333 0.50000000000000 85.468 > 109 125 Rb 0.50000000000000 0.83333333333333 0.50000000000000 85.468 > 109 126 Rb 0.83333333333333 0.83333333333333 0.50000000000000 85.468 > 109 127 Rb 0.16666666666667 0.16666666666667 0.83333333333333 85.468 > 109 128 Rb 0.50000000000000 0.16666666666667 0.83333333333333 85.468 > 109 129 Rb 0.83333333333333 0.16666666666667 0.83333333333333 85.468 > 109 130 Rb 0.16666666666667 0.50000000000000 0.83333333333333 85.468 > 109 131 Rb 0.50000000000000 0.50000000000000 0.83333333333333 85.468 > 109 132 Rb 0.83333333333333 0.50000000000000 0.83333333333333 85.468 > 109 133 Rb 0.16666666666667 0.83333333333333 0.83333333333333 85.468 > 109 134 Rb 0.50000000000000 0.83333333333333 0.83333333333333 85.468 > 109 135 Rb 0.83333333333333 0.83333333333333 0.83333333333333 85.468 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 2.2517795 0.0000000 0.0000000 0.0000000 2.2517795 0.0000000 0.0000000 0.0000000 2.2517795 -------------------------- Born effective charges -------------------------- 1 F -1.2627655 0.0000000 0.0000000 0.0000000 -1.0294943 0.0000000 0.0000000 0.0000000 -1.0294943 2 F -1.0294943 0.0000000 0.0000000 0.0000000 -1.2627655 0.0000000 0.0000000 0.0000000 -1.0294943 3 F -1.0294943 0.0000000 0.0000000 0.0000000 -1.0294943 0.0000000 0.0000000 0.0000000 -1.2627655 4 Mg 2.0696769 0.0000000 0.0000000 0.0000000 2.0696769 0.0000000 0.0000000 0.0000000 2.0696769 5 Rb 1.2520772 0.0000000 0.0000000 0.0000000 1.2520772 0.0000000 0.0000000 0.0000000 1.2520772 ---------------------------------------------------------------------------- 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.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 82, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 109, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (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 22:53:13]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 22:53:14]------------------------- 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.074999130000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.074999130000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.074999130000000 Atomic positions (fractional): 1 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 2 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 3 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 4 Mg 0.00000000000000 0.00000000000000 0.00000000000000 24.305 5 Rb 0.50000000000000 0.50000000000000 0.50000000000000 85.468 -------------------------------- supercell --------------------------------- Lattice vectors: a 12.224997390000000 0.000000000000000 0.000000000000000 b 0.000000000000000 12.224997390000000 0.000000000000000 c 0.000000000000000 0.000000000000000 12.224997390000000 Atomic positions (fractional): 1 F 0.16666666666667 0.00000000000000 0.00000000000000 18.998 > 1 2 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 1 3 F 0.83333333333333 0.00000000000000 0.00000000000000 18.998 > 1 4 F 0.16666666666667 0.33333333333333 0.00000000000000 18.998 > 1 5 F 0.50000000000000 0.33333333333333 0.00000000000000 18.998 > 1 6 F 0.83333333333333 0.33333333333333 0.00000000000000 18.998 > 1 7 F 0.16666666666667 0.66666666666667 0.00000000000000 18.998 > 1 8 F 0.50000000000000 0.66666666666667 0.00000000000000 18.998 > 1 9 F 0.83333333333333 0.66666666666667 0.00000000000000 18.998 > 1 10 F 0.16666666666667 0.00000000000000 0.33333333333333 18.998 > 1 11 F 0.50000000000000 0.00000000000000 0.33333333333333 18.998 > 1 12 F 0.83333333333333 0.00000000000000 0.33333333333333 18.998 > 1 13 F 0.16666666666667 0.33333333333333 0.33333333333333 18.998 > 1 14 F 0.50000000000000 0.33333333333333 0.33333333333333 18.998 > 1 15 F 0.83333333333333 0.33333333333333 0.33333333333333 18.998 > 1 16 F 0.16666666666667 0.66666666666667 0.33333333333333 18.998 > 1 17 F 0.50000000000000 0.66666666666667 0.33333333333333 18.998 > 1 18 F 0.83333333333333 0.66666666666667 0.33333333333333 18.998 > 1 19 F 0.16666666666667 0.00000000000000 0.66666666666667 18.998 > 1 20 F 0.50000000000000 0.00000000000000 0.66666666666667 18.998 > 1 21 F 0.83333333333333 0.00000000000000 0.66666666666667 18.998 > 1 22 F 0.16666666666667 0.33333333333333 0.66666666666667 18.998 > 1 23 F 0.50000000000000 0.33333333333333 0.66666666666667 18.998 > 1 24 F 0.83333333333333 0.33333333333333 0.66666666666667 18.998 > 1 25 F 0.16666666666667 0.66666666666667 0.66666666666667 18.998 > 1 26 F 0.50000000000000 0.66666666666667 0.66666666666667 18.998 > 1 27 F 0.83333333333333 0.66666666666667 0.66666666666667 18.998 > 1 28 F 0.00000000000000 0.16666666666667 0.00000000000000 18.998 > 28 29 F 0.33333333333333 0.16666666666667 0.00000000000000 18.998 > 28 30 F 0.66666666666667 0.16666666666667 0.00000000000000 18.998 > 28 31 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 28 32 F 0.33333333333333 0.50000000000000 0.00000000000000 18.998 > 28 33 F 0.66666666666667 0.50000000000000 0.00000000000000 18.998 > 28 34 F 0.00000000000000 0.83333333333333 0.00000000000000 18.998 > 28 35 F 0.33333333333333 0.83333333333333 0.00000000000000 18.998 > 28 36 F 0.66666666666667 0.83333333333333 0.00000000000000 18.998 > 28 37 F 0.00000000000000 0.16666666666667 0.33333333333333 18.998 > 28 38 F 0.33333333333333 0.16666666666667 0.33333333333333 18.998 > 28 39 F 0.66666666666667 0.16666666666667 0.33333333333333 18.998 > 28 40 F 0.00000000000000 0.50000000000000 0.33333333333333 18.998 > 28 41 F 0.33333333333333 0.50000000000000 0.33333333333333 18.998 > 28 42 F 0.66666666666667 0.50000000000000 0.33333333333333 18.998 > 28 43 F 0.00000000000000 0.83333333333333 0.33333333333333 18.998 > 28 44 F 0.33333333333333 0.83333333333333 0.33333333333333 18.998 > 28 45 F 0.66666666666667 0.83333333333333 0.33333333333333 18.998 > 28 46 F 0.00000000000000 0.16666666666667 0.66666666666667 18.998 > 28 47 F 0.33333333333333 0.16666666666667 0.66666666666667 18.998 > 28 48 F 0.66666666666667 0.16666666666667 0.66666666666667 18.998 > 28 49 F 0.00000000000000 0.50000000000000 0.66666666666667 18.998 > 28 50 F 0.33333333333333 0.50000000000000 0.66666666666667 18.998 > 28 51 F 0.66666666666667 0.50000000000000 0.66666666666667 18.998 > 28 52 F 0.00000000000000 0.83333333333333 0.66666666666667 18.998 > 28 53 F 0.33333333333333 0.83333333333333 0.66666666666667 18.998 > 28 54 F 0.66666666666667 0.83333333333333 0.66666666666667 18.998 > 28 55 F 0.00000000000000 0.00000000000000 0.16666666666667 18.998 > 55 56 F 0.33333333333333 0.00000000000000 0.16666666666667 18.998 > 55 57 F 0.66666666666667 0.00000000000000 0.16666666666667 18.998 > 55 58 F 0.00000000000000 0.33333333333333 0.16666666666667 18.998 > 55 59 F 0.33333333333333 0.33333333333333 0.16666666666667 18.998 > 55 60 F 0.66666666666667 0.33333333333333 0.16666666666667 18.998 > 55 61 F 0.00000000000000 0.66666666666667 0.16666666666667 18.998 > 55 62 F 0.33333333333333 0.66666666666667 0.16666666666667 18.998 > 55 63 F 0.66666666666667 0.66666666666667 0.16666666666667 18.998 > 55 64 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 55 65 F 0.33333333333333 0.00000000000000 0.50000000000000 18.998 > 55 66 F 0.66666666666667 0.00000000000000 0.50000000000000 18.998 > 55 67 F 0.00000000000000 0.33333333333333 0.50000000000000 18.998 > 55 68 F 0.33333333333333 0.33333333333333 0.50000000000000 18.998 > 55 69 F 0.66666666666667 0.33333333333333 0.50000000000000 18.998 > 55 70 F 0.00000000000000 0.66666666666667 0.50000000000000 18.998 > 55 71 F 0.33333333333333 0.66666666666667 0.50000000000000 18.998 > 55 72 F 0.66666666666667 0.66666666666667 0.50000000000000 18.998 > 55 73 F 0.00000000000000 0.00000000000000 0.83333333333333 18.998 > 55 74 F 0.33333333333333 0.00000000000000 0.83333333333333 18.998 > 55 75 F 0.66666666666667 0.00000000000000 0.83333333333333 18.998 > 55 76 F 0.00000000000000 0.33333333333333 0.83333333333333 18.998 > 55 77 F 0.33333333333333 0.33333333333333 0.83333333333333 18.998 > 55 78 F 0.66666666666667 0.33333333333333 0.83333333333333 18.998 > 55 79 F 0.00000000000000 0.66666666666667 0.83333333333333 18.998 > 55 80 F 0.33333333333333 0.66666666666667 0.83333333333333 18.998 > 55 81 F 0.66666666666667 0.66666666666667 0.83333333333333 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 Rb 0.16666666666667 0.16666666666667 0.16666666666667 85.468 > 109 110 Rb 0.50000000000000 0.16666666666667 0.16666666666667 85.468 > 109 111 Rb 0.83333333333333 0.16666666666667 0.16666666666667 85.468 > 109 112 Rb 0.16666666666667 0.50000000000000 0.16666666666667 85.468 > 109 113 Rb 0.50000000000000 0.50000000000000 0.16666666666667 85.468 > 109 114 Rb 0.83333333333333 0.50000000000000 0.16666666666667 85.468 > 109 115 Rb 0.16666666666667 0.83333333333333 0.16666666666667 85.468 > 109 116 Rb 0.50000000000000 0.83333333333333 0.16666666666667 85.468 > 109 117 Rb 0.83333333333333 0.83333333333333 0.16666666666667 85.468 > 109 118 Rb 0.16666666666667 0.16666666666667 0.50000000000000 85.468 > 109 119 Rb 0.50000000000000 0.16666666666667 0.50000000000000 85.468 > 109 120 Rb 0.83333333333333 0.16666666666667 0.50000000000000 85.468 > 109 121 Rb 0.16666666666667 0.50000000000000 0.50000000000000 85.468 > 109 122 Rb 0.50000000000000 0.50000000000000 0.50000000000000 85.468 > 109 123 Rb 0.83333333333333 0.50000000000000 0.50000000000000 85.468 > 109 124 Rb 0.16666666666667 0.83333333333333 0.50000000000000 85.468 > 109 125 Rb 0.50000000000000 0.83333333333333 0.50000000000000 85.468 > 109 126 Rb 0.83333333333333 0.83333333333333 0.50000000000000 85.468 > 109 127 Rb 0.16666666666667 0.16666666666667 0.83333333333333 85.468 > 109 128 Rb 0.50000000000000 0.16666666666667 0.83333333333333 85.468 > 109 129 Rb 0.83333333333333 0.16666666666667 0.83333333333333 85.468 > 109 130 Rb 0.16666666666667 0.50000000000000 0.83333333333333 85.468 > 109 131 Rb 0.50000000000000 0.50000000000000 0.83333333333333 85.468 > 109 132 Rb 0.83333333333333 0.50000000000000 0.83333333333333 85.468 > 109 133 Rb 0.16666666666667 0.83333333333333 0.83333333333333 85.468 > 109 134 Rb 0.50000000000000 0.83333333333333 0.83333333333333 85.468 > 109 135 Rb 0.83333333333333 0.83333333333333 0.83333333333333 85.468 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 2.2517795 0.0000000 0.0000000 0.0000000 2.2517795 0.0000000 0.0000000 0.0000000 2.2517795 -------------------------- Born effective charges -------------------------- 1 F -1.2627655 0.0000000 0.0000000 0.0000000 -1.0294943 0.0000000 0.0000000 0.0000000 -1.0294943 2 F -1.0294943 0.0000000 0.0000000 0.0000000 -1.2627655 0.0000000 0.0000000 0.0000000 -1.0294943 3 F -1.0294943 0.0000000 0.0000000 0.0000000 -1.0294943 0.0000000 0.0000000 0.0000000 -1.2627655 4 Mg 2.0696769 0.0000000 0.0000000 0.0000000 2.0696769 0.0000000 0.0000000 0.0000000 2.0696769 5 Rb 1.2520772 0.0000000 0.0000000 0.0000000 1.2520772 0.0000000 0.0000000 0.0000000 1.2520772 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (zzz) -0.00000000 (zzz) -0.00000000 (zzz) Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 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.02, 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.10e-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.523 ( 0.000 0.000 0.000) 0.000 4.523 ( 0.000 0.000 0.000) 0.000 4.523 ( 0.000 0.000 0.000) 0.000 6.212 ( 0.000 0.000 0.000) 0.000 6.212 ( 0.000 0.000 0.000) 0.000 6.212 ( 0.000 0.000 0.000) 0.000 8.408 ( 0.000 0.000 0.000) 0.000 8.408 ( 0.000 0.000 0.000) 0.000 8.408 ( 0.000 0.000 0.000) 0.000 11.929 ( 0.000 0.000 0.000) 0.000 11.929 ( 0.000 0.000 0.000) 0.000 11.929 ( 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.10e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.711 ( 34.204 0.000 0.000) 34.204 0.711 ( 34.204 0.000 0.000) 34.204 1.121 ( 52.702 0.000 0.000) 52.702 4.524 ( 0.129 0.000 0.000) 0.129 4.524 ( 0.129 0.000 0.000) 0.129 5.126 ( 12.978 0.000 0.000) 12.978 6.232 ( 1.884 0.000 0.000) 1.884 6.232 ( 1.884 0.000 0.000) 1.884 6.234 ( 2.054 0.000 0.000) 2.054 8.386 ( -2.125 0.000 0.000) 2.125 8.386 ( -2.125 0.000 0.000) 2.125 10.035 ( 3.151 0.000 0.000) 3.151 11.924 ( -0.452 0.000 0.000) 0.452 11.924 ( -0.452 0.000 0.000) 0.452 14.820 ( -6.695 0.000 0.000) 6.695 ======================= 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.10e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.387 ( 31.688 0.000 0.000) 31.688 1.387 ( 31.688 0.000 0.000) 31.688 2.119 ( 44.269 0.000 0.000) 44.269 4.524 ( -0.338 0.000 0.000) 0.338 4.524 ( -0.338 0.000 0.000) 0.338 5.495 ( 22.047 0.000 0.000) 22.047 6.289 ( 3.745 0.000 0.000) 3.745 6.289 ( 3.745 0.000 0.000) 3.745 6.292 ( 3.517 0.000 0.000) 3.517 8.323 ( -4.002 0.000 0.000) 4.002 8.323 ( -4.002 0.000 0.000) 4.002 10.145 ( 7.949 0.000 0.000) 7.949 11.911 ( -0.785 0.000 0.000) 0.785 11.911 ( -0.785 0.000 0.000) 0.785 14.615 ( -13.358 0.000 0.000) 13.358 ======================= 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.10e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.995 ( 27.442 0.000 0.000) 27.442 1.995 ( 27.442 0.000 0.000) 27.442 2.916 ( 33.559 0.000 0.000) 33.559 4.506 ( -1.523 0.000 0.000) 1.523 4.506 ( -1.523 0.000 0.000) 1.523 5.974 ( 23.366 0.000 0.000) 23.366 6.371 ( 3.999 0.000 0.000) 3.999 6.382 ( 5.245 0.000 0.000) 5.245 6.382 ( 5.245 0.000 0.000) 5.245 8.227 ( -5.287 0.000 0.000) 5.287 8.227 ( -5.287 0.000 0.000) 5.287 10.379 ( 15.299 0.000 0.000) 15.299 11.893 ( -0.908 0.000 0.000) 0.908 11.893 ( -0.908 0.000 0.000) 0.908 14.276 ( -19.650 0.000 0.000) 19.650 ======================= 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.10e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.496 ( 21.197 0.000 0.000) 21.197 2.496 ( 21.197 0.000 0.000) 21.197 3.491 ( 22.700 0.000 0.000) 22.700 4.461 ( -2.780 0.000 0.000) 2.780 4.461 ( -2.780 0.000 0.000) 2.780 6.397 ( 17.008 0.000 0.000) 17.008 6.448 ( 3.412 0.000 0.000) 3.412 6.496 ( 5.602 0.000 0.000) 5.602 6.496 ( 5.602 0.000 0.000) 5.602 8.115 ( -5.421 0.000 0.000) 5.421 8.115 ( -5.421 0.000 0.000) 5.421 10.775 ( 23.065 0.000 0.000) 23.065 11.875 ( -0.788 0.000 0.000) 0.788 11.875 ( -0.788 0.000 0.000) 0.788 13.821 ( -24.347 0.000 0.000) 24.347 ======================= 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.10e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.843 ( 12.179 0.000 0.000) 12.179 2.843 ( 12.179 0.000 0.000) 12.179 3.843 ( 11.630 0.000 0.000) 11.630 4.401 ( -2.643 0.000 0.000) 2.643 4.401 ( -2.643 0.000 0.000) 2.643 6.504 ( 1.949 0.000 0.000) 1.949 6.597 ( 3.865 0.000 0.000) 3.865 6.597 ( 3.865 0.000 0.000) 3.865 6.653 ( 8.048 0.000 0.000) 8.048 8.018 ( -3.672 0.000 0.000) 3.672 8.018 ( -3.672 0.000 0.000) 3.672 11.274 ( 23.484 0.000 0.000) 23.484 11.862 ( -0.456 0.000 0.000) 0.456 11.862 ( -0.456 0.000 0.000) 0.456 13.319 ( -22.922 0.000 0.000) 22.922 ======================= 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.10e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.972 ( -0.000 0.000 0.000) 0.000 2.972 ( -0.000 0.000 0.000) 0.000 3.963 ( -0.000 0.000 0.000) 0.000 4.371 ( 0.000 0.000 0.000) 0.000 4.371 ( 0.000 0.000 0.000) 0.000 6.525 ( 0.000 0.000 0.000) 0.000 6.639 ( 0.000 0.000 0.000) 0.000 6.639 ( 0.000 0.000 0.000) 0.000 6.733 ( 0.000 0.000 0.000) 0.000 7.978 ( 0.000 0.000 0.000) 0.000 7.978 ( 0.000 0.000 0.000) 0.000 11.567 ( 0.000 0.000 0.000) 0.000 11.858 ( 0.000 0.000 0.000) 0.000 11.858 ( 0.000 0.000 0.000) 0.000 13.035 ( -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.10e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.889 ( 20.771 20.771 0.000) 29.374 1.000 ( 23.760 23.760 0.000) 33.602 1.626 ( 37.174 37.174 0.000) 52.572 4.530 ( 0.482 0.482 0.000) 0.682 4.648 ( 5.817 5.817 0.000) 8.226 5.067 ( 3.844 3.844 0.000) 5.436 6.175 ( -1.897 -1.897 0.000) 2.683 6.250 ( 1.778 1.778 0.000) 2.515 6.329 ( 5.674 5.674 0.000) 8.025 8.364 ( -2.196 -2.196 0.000) 3.105 8.517 ( 5.309 5.309 0.000) 7.508 10.062 ( 2.762 2.762 0.000) 3.906 11.790 ( -6.601 -6.601 0.000) 9.336 11.919 ( -0.453 -0.453 0.000) 0.641 14.800 ( -4.313 -4.313 0.000) 6.099 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.422 ( 28.412 6.139 0.000) 29.067 1.543 ( 27.650 14.189 0.000) 31.078 2.403 ( 36.398 23.167 0.000) 43.145 4.540 ( 0.377 1.557 0.000) 1.602 4.673 ( -1.083 13.417 0.000) 13.460 5.370 ( 21.774 -8.443 0.000) 23.353 6.155 ( 0.023 -5.814 0.000) 5.814 6.304 ( 3.561 1.483 0.000) 3.858 6.464 ( 7.404 9.312 0.000) 11.897 8.298 ( -4.145 -2.410 0.000) 4.794 8.541 ( -2.107 16.584 0.000) 16.718 10.159 ( 7.020 1.596 0.000) 7.199 11.713 ( -1.607 -15.647 0.000) 15.729 11.906 ( -0.786 -0.455 0.000) 0.908 14.626 ( -12.447 -0.593 0.000) 12.461 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.992 ( 26.570 0.983 0.000) 26.588 2.089 ( 25.109 8.759 0.000) 26.592 3.061 ( 27.318 11.959 0.000) 29.821 4.540 ( -0.470 3.348 0.000) 3.381 4.650 ( -0.912 13.639 0.000) 13.670 5.829 ( 20.182 -10.995 0.000) 22.983 6.188 ( 4.097 -7.814 0.000) 8.823 6.394 ( 5.054 1.089 0.000) 5.170 6.631 ( 9.029 13.235 0.000) 16.021 8.198 ( -5.495 -2.764 0.000) 6.151 8.458 ( -5.551 19.962 0.000) 20.720 10.360 ( 12.797 -1.436 0.000) 12.878 11.711 ( 1.465 -16.155 0.000) 16.221 11.888 ( -0.909 -0.457 0.000) 1.018 14.299 ( -19.315 1.247 0.000) 19.355 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.486 ( 21.321 -0.444 0.000) 21.325 2.550 ( 19.442 4.947 0.000) 20.061 3.508 ( 16.296 1.316 0.000) 16.349 4.519 ( -1.512 5.703 0.000) 5.900 4.639 ( -0.153 16.270 0.000) 16.270 6.073 ( 4.584 -18.250 0.000) 18.817 6.376 ( 12.222 -3.718 0.000) 12.775 6.504 ( 5.493 0.768 0.000) 5.546 6.834 ( 10.648 17.639 0.000) 20.604 8.081 ( -5.667 -3.231 0.000) 6.524 8.327 ( -6.923 19.200 0.000) 20.410 10.667 ( 16.163 -9.648 0.000) 18.824 11.789 ( 6.965 -8.344 0.000) 10.869 11.871 ( -0.788 -0.459 0.000) 0.912 13.844 ( -24.914 1.748 0.000) 24.975 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.841 ( 12.681 -0.107 0.000) 12.682 2.865 ( 10.965 2.054 0.000) 11.156 3.733 ( 6.114 -7.718 0.000) 9.846 4.485 ( -1.586 8.079 0.000) 8.233 4.645 ( 0.653 19.957 0.000) 19.968 6.123 ( 1.256 -21.565 0.000) 21.602 6.569 ( 6.023 -1.762 0.000) 6.276 6.604 ( 3.855 0.649 0.000) 3.909 7.040 ( 8.374 21.534 0.000) 23.105 7.980 ( -3.868 -3.702 0.000) 5.354 8.194 ( -5.503 16.510 0.000) 17.403 10.953 ( 10.099 -23.754 0.000) 25.812 11.858 ( -0.456 -0.460 0.000) 0.647 12.031 ( 16.783 8.414 0.000) 18.774 13.305 ( -26.664 -0.677 0.000) 26.673 ======================= 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.10e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.976 ( 0.000 0.453 0.000) 0.453 2.980 ( -0.000 0.816 0.000) 0.816 3.789 ( 0.000 -11.326 0.000) 11.326 4.467 ( 0.000 9.196 0.000) 9.196 4.654 ( 0.000 21.410 0.000) 21.410 6.135 ( 0.000 -22.700 0.000) 22.700 6.628 ( 0.000 -0.959 0.000) 0.959 6.646 ( 0.000 0.653 0.000) 0.653 7.133 ( 0.000 23.605 0.000) 23.605 7.938 ( 0.000 -3.922 0.000) 3.922 8.132 ( 0.000 14.700 0.000) 14.700 11.054 ( 0.000 -30.411 0.000) 30.411 11.853 ( 0.000 -0.460 0.000) 0.460 12.318 ( 0.000 25.298 0.000) 25.298 12.921 ( -0.000 -10.848 0.000) 10.848 ======================= 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.10e-04 Number of triplets: 301 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.672 ( 17.340 17.340 0.000) 24.523 1.909 ( 20.159 20.159 0.000) 28.509 2.906 ( 23.554 23.554 0.000) 33.310 4.584 ( 2.569 2.569 0.000) 3.633 4.949 ( 8.038 8.038 0.000) 11.368 5.331 ( 9.535 9.535 0.000) 13.485 6.053 ( -3.961 -3.961 0.000) 5.602 6.350 ( 3.021 3.021 0.000) 4.272 6.671 ( 11.044 11.044 0.000) 15.619 8.226 ( -4.576 -4.576 0.000) 6.472 8.829 ( 9.486 9.486 0.000) 13.415 10.225 ( 5.088 5.088 0.000) 7.195 11.418 ( -11.004 -11.004 0.000) 15.562 11.893 ( -0.789 -0.789 0.000) 1.116 14.538 ( -8.468 -8.468 0.000) 11.976 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.078 ( 21.245 7.808 0.000) 22.635 2.323 ( 19.523 13.233 0.000) 23.585 3.310 ( 15.367 10.832 0.000) 18.801 4.639 ( 2.729 6.214 0.000) 6.787 5.001 ( 1.023 17.551 0.000) 17.580 5.660 ( 13.933 -2.514 0.000) 14.158 6.035 ( 5.057 -7.287 0.000) 8.869 6.427 ( 4.446 2.227 0.000) 4.972 6.917 ( 13.072 14.882 0.000) 19.808 8.115 ( -6.134 -5.300 0.000) 8.107 8.907 ( -1.338 20.862 0.000) 20.905 10.352 ( 6.337 1.909 0.000) 6.618 11.320 ( 2.518 -20.153 0.000) 20.310 11.875 ( -0.911 -0.792 0.000) 1.207 14.279 ( -16.679 -4.085 0.000) 17.172 ======================= 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.10e-04 Number of triplets: 518 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.503 ( 19.557 2.563 0.000) 19.724 2.681 ( 14.977 7.381 0.000) 16.697 3.529 ( 6.266 0.426 0.000) 6.280 4.690 ( 2.114 10.786 0.000) 10.991 5.042 ( 2.882 20.030 0.000) 20.237 5.747 ( -1.452 -11.860 0.000) 11.949 6.290 ( 14.753 -4.298 0.000) 15.367 6.527 ( 5.106 1.487 0.000) 5.318 7.204 ( 14.639 18.829 0.000) 23.851 7.984 ( -6.442 -6.282 0.000) 8.998 8.808 ( -7.428 25.303 0.000) 26.371 10.442 ( 2.223 -10.212 0.000) 10.451 11.545 ( 19.117 -15.038 0.000) 24.323 11.857 ( -0.790 -0.793 0.000) 1.119 13.863 ( -23.763 -0.668 0.000) 23.772 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.847 ( 13.122 0.896 0.000) 13.153 2.920 ( 8.082 3.044 0.000) 8.636 3.586 ( 0.035 -6.452 0.000) 6.452 4.723 ( 1.094 14.874 0.000) 14.914 5.104 ( 2.687 22.756 0.000) 22.914 5.703 ( -2.084 -18.543 0.000) 18.660 6.523 ( 7.532 -2.577 0.000) 7.960 6.623 ( 3.800 1.156 0.000) 3.972 7.485 ( 11.645 22.076 0.000) 24.959 7.867 ( -4.505 -7.320 0.000) 8.595 8.642 ( -7.702 25.562 0.000) 26.697 10.461 ( 0.131 -22.302 0.000) 22.302 11.844 ( -0.456 -0.793 0.000) 0.915 12.051 ( 28.871 -6.978 0.000) 29.702 13.319 ( -29.013 1.920 0.000) 29.077 ======================= 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.10e-04 Number of triplets: 301 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.993 ( 0.000 1.231 0.000) 1.231 3.003 ( 0.000 1.366 0.000) 1.366 3.575 ( 0.000 -9.367 0.000) 9.367 4.734 ( 0.000 16.561 0.000) 16.561 5.134 ( 0.000 24.340 0.000) 24.340 5.679 ( 0.000 -21.577 0.000) 21.577 6.599 ( 0.000 -1.913 0.000) 1.913 6.665 ( 0.000 1.146 0.000) 1.146 7.618 ( 0.000 23.645 0.000) 23.645 7.818 ( 0.000 -7.828 0.000) 7.828 8.550 ( 0.000 24.922 0.000) 24.922 10.461 ( 0.000 -26.398 0.000) 26.398 11.840 ( 0.000 -0.793 0.000) 0.793 12.612 ( 0.000 -18.681 0.000) 18.681 12.761 ( -0.000 17.260 0.000) 17.260 ======================= 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.10e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.301 ( 13.390 13.390 0.000) 18.936 2.598 ( 12.952 12.952 0.000) 18.317 3.474 ( 4.990 4.990 0.000) 7.057 4.786 ( 7.903 7.903 0.000) 11.177 5.210 ( 2.927 2.927 0.000) 4.140 5.821 ( 12.675 12.675 0.000) 17.925 5.902 ( -1.843 -1.843 0.000) 2.606 6.485 ( 3.400 3.400 0.000) 4.808 7.244 ( 17.088 17.088 0.000) 24.167 7.985 ( -7.240 -7.240 0.000) 10.238 9.207 ( 7.121 7.121 0.000) 10.070 10.467 ( 6.624 6.624 0.000) 9.368 10.978 ( -8.706 -8.706 0.000) 12.312 11.857 ( -0.913 -0.913 0.000) 1.291 14.115 ( -12.023 -12.023 0.000) 17.003 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.599 ( 15.008 6.675 0.000) 16.425 2.832 ( 9.618 7.015 0.000) 11.904 3.517 ( -0.543 -1.716 0.000) 1.800 4.948 ( 7.513 14.051 0.000) 15.934 5.234 ( 0.760 -3.863 0.000) 3.937 5.751 ( -5.717 13.413 0.000) 14.581 6.214 ( 16.230 -2.733 0.000) 16.459 6.565 ( 4.274 2.167 0.000) 4.792 7.607 ( 17.882 20.445 0.000) 27.162 7.828 ( -7.856 -8.821 0.000) 11.812 9.223 ( -3.955 11.067 0.000) 11.753 10.382 ( -11.687 8.302 0.000) 14.336 11.198 ( 25.508 -18.262 0.000) 31.372 11.840 ( -0.790 -0.913 0.000) 1.207 13.790 ( -19.530 -6.894 0.000) 20.711 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.884 ( 12.003 2.820 0.000) 12.330 2.982 ( 4.977 2.941 0.000) 5.781 3.466 ( -3.967 -5.454 0.000) 6.744 5.076 ( 4.610 19.176 0.000) 19.722 5.252 ( 0.859 -13.679 0.000) 13.706 5.654 ( -3.509 18.659 0.000) 18.986 6.467 ( 8.297 -2.772 0.000) 8.748 6.649 ( 3.608 1.419 0.000) 3.877 7.683 ( -5.771 -10.652 0.000) 12.115 7.940 ( 13.436 22.181 0.000) 25.933 9.107 ( -6.095 14.039 0.000) 15.305 10.162 ( -8.231 -1.583 0.000) 8.382 11.790 ( 29.195 -16.944 0.000) 33.756 11.827 ( -0.455 -0.912 0.000) 1.019 13.340 ( -23.326 -1.152 0.000) 23.354 ======================= 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.10e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.028 ( 0.000 2.268 0.000) 2.268 3.033 ( 0.000 1.504 0.000) 1.504 3.405 ( 0.000 -7.297 0.000) 7.297 5.125 ( 0.000 21.177 0.000) 21.177 5.262 ( 0.000 -19.016 0.000) 19.016 5.617 ( 0.000 22.312 0.000) 22.312 6.552 ( 0.000 -2.602 0.000) 2.602 6.691 ( 0.000 1.345 0.000) 1.345 7.618 ( 0.000 -11.644 0.000) 11.644 8.092 ( 0.000 22.517 0.000) 22.517 9.031 ( 0.000 16.086 0.000) 16.086 10.075 ( 0.000 -6.089 0.000) 6.089 11.822 ( 0.000 -0.911 0.000) 0.911 12.188 ( 0.000 -21.997 0.000) 21.997 13.005 ( -0.000 6.378 0.000) 6.378 ======================= 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.10e-04 Number of triplets: 301 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.766 ( 9.345 9.345 0.000) 13.216 2.957 ( 5.064 5.064 0.000) 7.162 3.459 ( -3.955 -3.955 0.000) 5.593 5.087 ( -7.623 -7.623 0.000) 10.780 5.241 ( 13.904 13.904 0.000) 19.663 6.011 ( 6.212 6.212 0.000) 8.786 6.230 ( 6.813 6.813 0.000) 9.635 6.616 ( 2.873 2.873 0.000) 4.063 7.631 ( -10.078 -10.078 0.000) 14.252 8.024 ( 19.661 19.661 0.000) 27.805 9.236 ( -6.183 -6.183 0.000) 8.744 10.757 ( 7.445 7.445 0.000) 10.529 10.867 ( 3.616 3.616 0.000) 5.113 11.822 ( -0.788 -0.788 0.000) 1.114 13.575 ( -14.041 -14.041 0.000) 19.857 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.961 ( 9.070 4.569 0.000) 10.155 3.036 ( 2.583 2.200 0.000) 3.393 3.360 ( -5.437 -4.865 0.000) 7.296 4.960 ( -4.507 -12.958 0.000) 13.719 5.484 ( 9.026 20.088 0.000) 22.023 6.037 ( -0.532 16.660 0.000) 16.668 6.415 ( 7.763 -2.173 0.000) 8.062 6.679 ( 3.042 1.463 0.000) 3.375 7.437 ( -8.116 -13.101 0.000) 15.412 8.376 ( 13.630 19.876 0.000) 24.100 9.093 ( -6.145 -10.455 0.000) 12.127 10.493 ( -12.202 28.254 0.000) 30.776 11.417 ( 26.959 -18.297 0.000) 32.582 11.809 ( -0.454 -0.786 0.000) 0.907 13.243 ( -16.853 -8.403 0.000) 18.832 ======================= 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.10e-04 Number of triplets: 301 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.062 ( 0.000 1.242 0.000) 1.242 3.083 ( 0.000 2.979 0.000) 2.979 3.277 ( 0.000 -5.173 0.000) 5.173 4.912 ( 0.000 -14.852 0.000) 14.852 5.581 ( 0.000 22.975 0.000) 22.975 6.030 ( 0.000 17.722 0.000) 17.722 6.497 ( 0.000 -2.582 0.000) 2.582 6.717 ( 0.000 1.184 0.000) 1.184 7.344 ( 0.000 -15.127 0.000) 15.127 8.527 ( 0.000 19.528 0.000) 19.528 9.025 ( 0.000 -11.809 0.000) 11.809 10.358 ( 0.000 27.938 0.000) 27.938 11.745 ( 0.000 -20.500 0.000) 20.500 11.804 ( 0.000 -0.784 0.000) 0.784 13.022 ( 0.000 -4.322 0.000) 4.322 ======================= 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.10e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.062 ( 5.034 5.034 0.000) 7.120 3.070 ( 1.158 1.158 0.000) 1.638 3.266 ( -4.363 -4.363 0.000) 6.170 4.747 ( -7.431 -7.431 0.000) 10.508 5.862 ( 15.628 15.628 0.000) 22.101 6.289 ( 6.277 6.277 0.000) 8.877 6.400 ( 2.036 2.036 0.000) 2.879 6.710 ( 1.630 1.630 0.000) 2.306 7.166 ( -12.457 -12.457 0.000) 17.617 8.722 ( 13.036 13.036 0.000) 18.435 8.870 ( -9.077 -9.077 0.000) 12.837 11.063 ( 7.108 7.108 0.000) 10.053 11.109 ( 5.925 5.925 0.000) 8.380 11.796 ( -0.452 -0.452 0.000) 0.639 13.019 ( -12.079 -12.079 0.000) 17.082 ======================= 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.10e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.082 ( 0.000 0.692 0.000) 0.692 3.141 ( 0.000 2.372 0.000) 2.372 3.196 ( -0.000 -2.719 0.000) 2.719 4.669 ( 0.000 -8.515 0.000) 8.515 6.044 ( 0.000 21.574 0.000) 21.574 6.324 ( 0.000 10.501 0.000) 10.501 6.453 ( 0.000 -1.625 0.000) 1.625 6.737 ( 0.000 0.691 0.000) 0.691 7.008 ( 0.000 -17.196 0.000) 17.196 8.771 ( 0.000 -10.350 0.000) 10.350 8.864 ( 0.000 12.590 0.000) 12.590 10.947 ( 0.000 25.443 0.000) 25.443 11.387 ( 0.000 -13.391 0.000) 13.391 11.791 ( 0.000 -0.451 0.000) 0.451 12.862 ( 0.000 -9.814 0.000) 9.814 ======================= 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.10e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.090 ( 0.000 0.000 0.000) 0.000 3.168 ( 0.000 0.000 0.000) 0.000 3.168 ( 0.000 0.000 0.000) 0.000 4.580 ( 0.000 0.000 0.000) 0.000 6.354 ( 0.000 0.000 0.000) 0.000 6.436 ( 0.000 0.000 0.000) 0.000 6.436 ( 0.000 0.000 0.000) 0.000 6.744 ( 0.000 0.000 0.000) 0.000 6.744 ( 0.000 0.000 0.000) 0.000 8.659 ( 0.000 0.000 0.000) 0.000 9.001 ( 0.000 0.000 0.000) 0.000 11.242 ( 0.000 0.000 0.000) 0.000 11.242 ( 0.000 0.000 0.000) 0.000 11.786 ( 0.000 0.000 0.000) 0.000 12.730 ( 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.10e-04 Number of triplets: 182 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.126 ( 17.342 17.342 17.342) 30.037 1.126 ( 17.342 17.342 17.342) 30.037 1.989 ( 29.221 29.221 29.221) 50.613 4.648 ( 3.751 3.751 3.751) 6.496 4.648 ( 3.751 3.751 3.751) 6.496 5.027 ( 2.046 2.046 2.046) 3.544 6.119 ( -3.042 -3.042 -3.042) 5.268 6.345 ( 4.160 4.160 4.160) 7.205 6.345 ( 4.160 4.160 4.160) 7.205 8.493 ( 2.704 2.704 2.704) 4.684 8.493 ( 2.704 2.704 2.704) 4.684 10.083 ( 2.308 2.308 2.308) 3.997 11.786 ( -4.525 -4.525 -4.525) 7.838 11.786 ( -4.525 -4.525 -4.525) 7.838 14.784 ( -3.252 -3.252 -3.252) 5.633 ======================= 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.10e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.541 ( 24.096 9.394 9.394) 27.516 1.615 ( 25.686 9.857 9.857) 29.225 2.626 ( 30.382 18.657 18.657) 40.240 4.642 ( -0.497 5.404 5.404) 7.659 4.713 ( 0.637 7.616 7.616) 10.790 5.295 ( 21.089 -4.727 -4.727) 22.123 6.071 ( -1.640 -5.407 -5.407) 7.821 6.398 ( 3.928 5.105 5.105) 8.219 6.502 ( 8.459 4.954 4.954) 10.984 8.424 ( -4.389 4.939 4.939) 8.250 8.575 ( 0.129 7.345 7.345) 10.389 10.167 ( 6.175 1.094 1.094) 6.366 11.664 ( -3.037 -8.327 -8.327) 12.161 11.774 ( -0.716 -6.527 -6.527) 9.259 14.634 ( -11.376 -0.496 -0.496) 11.397 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.047 ( 24.358 3.876 3.876) 24.967 2.132 ( 24.116 6.133 6.133) 25.628 3.168 ( 21.840 8.806 8.806) 25.141 4.624 ( -1.330 5.340 5.340) 7.668 4.738 ( 2.123 10.787 10.787) 15.403 5.747 ( 20.524 -6.585 -6.585) 22.538 6.051 ( -0.422 -7.939 -7.939) 11.235 6.523 ( 8.930 5.613 5.613) 11.948 6.684 ( 9.119 7.022 7.022) 13.482 8.318 ( -5.889 4.480 4.480) 8.650 8.518 ( -5.163 10.931 10.931) 16.298 10.341 ( 10.692 -1.544 -1.544) 10.912 11.660 ( 2.788 -9.762 -9.762) 14.084 11.758 ( -0.830 -6.454 -6.454) 9.164 14.326 ( -18.666 1.565 1.565) 18.796 ======================= 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.10e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.507 ( 20.028 1.224 1.224) 20.103 2.578 ( 19.030 3.722 3.722) 19.744 3.509 ( 11.625 0.042 0.042) 11.625 4.588 ( -2.189 5.673 5.673) 8.316 4.798 ( 3.499 15.201 15.201) 21.780 6.050 ( 0.313 -10.267 -10.267) 14.522 6.063 ( 10.182 -8.652 -8.652) 15.919 6.772 ( 14.554 7.770 7.770) 18.236 6.866 ( 8.426 8.908 8.908) 15.156 8.191 ( -6.266 3.812 3.812) 8.266 8.382 ( -7.648 10.971 10.971) 17.298 10.581 ( 11.736 -7.864 -7.864) 16.169 11.742 ( -0.721 -6.380 -6.380) 9.052 11.797 ( 11.442 -4.844 -4.844) 13.336 13.878 ( -24.829 2.755 2.755) 25.133 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.841 ( 11.967 0.171 0.171) 11.969 2.889 ( 10.868 2.134 2.134) 11.279 3.660 ( 3.764 -5.737 -5.737) 8.944 4.542 ( -1.983 6.268 6.268) 9.083 4.866 ( 2.758 18.668 18.668) 26.544 6.059 ( 0.404 -11.931 -11.931) 16.878 6.195 ( 3.598 -10.238 -10.238) 14.919 7.014 ( 5.573 10.506 10.506) 15.868 7.051 ( 11.151 10.535 10.535) 18.609 8.076 ( -4.520 2.972 2.972) 6.172 8.228 ( -6.612 9.206 9.206) 14.602 10.774 ( 6.443 -15.979 -15.979) 23.498 11.730 ( -0.417 -6.327 -6.327) 8.958 12.148 ( 22.588 3.701 3.701) 23.187 13.325 ( -28.752 2.438 2.438) 28.958 ======================= 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.10e-04 Number of triplets: 276 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.968 ( 0.000 -0.011 -0.011) 0.015 3.003 ( -0.000 1.511 1.511) 2.137 3.693 ( 0.000 -7.653 -7.653) 10.823 4.519 ( 0.000 6.602 6.602) 9.337 4.897 ( 0.000 19.784 19.784) 27.979 6.064 ( 0.000 -12.537 -12.537) 17.729 6.230 ( 0.000 -10.871 -10.871) 15.374 7.074 ( 0.000 11.251 11.251) 15.911 7.175 ( 0.000 12.242 12.242) 17.313 8.026 ( 0.000 2.484 2.484) 3.513 8.152 ( 0.000 7.853 7.853) 11.105 10.839 ( 0.000 -19.197 -19.197) 27.149 11.725 ( 0.000 -6.308 -6.308) 8.921 12.601 ( 0.000 19.740 19.740) 27.917 12.812 ( -0.000 -10.348 -10.348) 14.634 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.788 ( 15.586 15.586 10.826) 24.557 1.938 ( 19.163 19.163 3.630) 27.343 3.034 ( 18.909 18.909 10.699) 28.802 4.696 ( 2.097 2.097 9.871) 10.307 4.937 ( 7.665 7.665 -1.066) 10.892 5.329 ( 11.778 11.778 1.830) 16.757 5.969 ( -4.405 -4.405 -6.181) 8.776 6.472 ( 3.459 3.459 9.731) 10.892 6.675 ( 10.709 10.709 0.452) 15.151 8.399 ( -4.299 -4.299 14.769) 15.971 8.800 ( 9.350 9.350 -2.802) 13.516 10.219 ( 4.232 4.232 -0.308) 5.993 11.414 ( -10.972 -10.972 -0.305) 15.520 11.725 ( -0.953 -0.953 -14.460) 14.523 14.557 ( -7.667 -7.667 0.784) 10.870 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.155 ( 19.334 7.715 7.145) 22.009 2.340 ( 19.254 12.758 2.133) 23.196 3.349 ( 11.483 7.699 3.137) 14.176 4.734 ( 1.317 6.079 7.549) 9.781 5.035 ( 4.724 13.752 3.384) 14.929 5.669 ( 14.575 2.311 2.498) 14.967 5.909 ( 0.294 -6.031 -8.969) 10.813 6.578 ( 7.678 1.599 11.597) 14.000 6.917 ( 12.787 14.095 0.068) 19.031 8.286 ( -6.467 -4.794 14.999) 17.023 8.886 ( -0.914 19.436 -1.823) 19.543 10.327 ( 5.229 1.315 -2.201) 5.824 11.309 ( 2.015 -19.577 -1.135) 19.713 11.711 ( -0.386 -1.090 -14.583) 14.628 14.316 ( -15.737 -3.337 2.692) 16.311 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.544 ( 17.895 3.033 3.715) 18.527 2.698 ( 15.202 7.238 1.880) 16.942 3.506 ( 4.241 -0.456 -1.905) 4.671 4.746 ( -0.034 9.825 4.315) 10.731 5.160 ( 7.097 15.961 10.511) 20.386 5.748 ( -2.145 -10.346 -0.028) 10.566 6.041 ( 7.481 0.457 -13.481) 15.424 6.792 ( 12.195 -0.688 16.752) 20.732 7.196 ( 14.330 18.477 -0.770) 23.396 8.144 ( -7.199 -5.659 14.235) 16.926 8.791 ( -7.492 23.846 -1.305) 25.029 10.392 ( 1.002 -8.832 -4.770) 10.088 11.502 ( 13.978 -15.303 -3.937) 21.097 11.731 ( 4.647 -2.521 -11.435) 12.598 13.921 ( -22.705 0.470 4.717) 23.194 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.855 ( 11.605 1.340 0.893) 11.716 2.945 ( 8.645 3.110 2.286) 9.468 3.541 ( -0.245 -5.487 -3.767) 6.660 4.740 ( -0.361 12.800 2.236) 12.999 5.297 ( 5.587 20.100 16.349) 26.505 5.695 ( -2.504 -16.732 -0.876) 16.941 6.137 ( 2.402 -0.447 -19.420) 19.573 7.016 ( 8.529 -0.945 20.997) 22.683 7.477 ( 11.934 21.822 -0.758) 24.883 8.007 ( -5.567 -6.738 12.847) 15.539 8.620 ( -8.034 24.507 -1.769) 25.850 10.390 ( -0.537 -19.120 -6.898) 20.333 11.634 ( 1.014 -3.614 -14.183) 14.672 12.086 ( 26.093 -9.381 -1.776) 27.785 13.405 ( -26.932 4.539 7.591) 28.347 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.979 ( 0.000 1.094 -0.194) 1.111 3.039 ( 0.000 1.881 2.302) 2.972 3.530 ( 0.000 -7.675 -3.856) 8.589 4.735 ( 0.000 14.082 1.500) 14.162 5.360 ( 0.000 24.028 19.435) 30.904 5.665 ( 0.000 -21.127 -1.650) 21.191 6.159 ( 0.000 -1.241 -22.186) 22.221 7.108 ( 0.000 -0.356 22.637) 22.640 7.615 ( 0.000 23.118 -0.135) 23.118 7.944 ( 0.000 -7.428 11.813) 13.954 8.523 ( 0.000 24.155 -2.292) 24.263 10.383 ( 0.000 -22.692 -7.649) 23.946 11.640 ( 0.000 -2.671 -14.744) 14.984 12.517 ( 0.000 -17.841 -9.013) 19.988 12.973 ( -0.000 15.192 16.214) 22.219 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.366 ( 12.496 12.496 6.170) 18.718 2.611 ( 13.012 13.012 1.440) 18.458 3.461 ( 3.231 3.231 -1.065) 4.692 4.869 ( 6.643 6.643 6.321) 11.324 5.189 ( 3.093 3.093 -1.989) 4.805 5.786 ( -3.277 -3.277 -9.123) 10.233 5.901 ( 14.020 14.020 9.367) 21.928 6.629 ( 4.183 4.183 11.905) 13.294 7.235 ( 16.837 16.837 -0.913) 23.828 8.159 ( -7.348 -7.348 15.849) 18.952 9.175 ( 7.151 7.151 -3.139) 10.589 10.421 ( 5.574 5.574 -4.204) 8.933 10.975 ( -8.778 -8.778 -0.290) 12.417 11.692 ( -0.764 -0.764 -14.950) 14.989 14.171 ( -11.002 -11.002 4.431) 16.178 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.642 ( 13.833 6.565 4.146) 15.864 2.850 ( 9.910 7.142 1.699) 12.332 3.485 ( -0.797 -1.752 -2.739) 3.347 4.982 ( 4.176 12.899 1.797) 13.677 5.244 ( 2.434 -6.380 0.633) 6.858 5.733 ( -0.560 5.248 -5.056) 7.309 6.131 ( 5.964 9.377 2.001) 11.291 6.778 ( 10.182 -0.362 15.425) 18.486 7.594 ( 17.720 20.196 -1.324) 26.900 7.995 ( -8.362 -8.541 15.443) 19.528 9.192 ( -3.946 11.226 -2.977) 12.266 10.342 ( -11.101 7.818 -3.783) 14.094 11.154 ( 22.241 -17.983 -3.977) 28.877 11.687 ( 0.714 -1.857 -13.727) 13.871 13.870 ( -18.069 -5.775 6.671) 20.108 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.900 ( 10.410 3.165 1.648) 11.004 3.009 ( 5.653 2.954 2.402) 6.816 3.438 ( -3.362 -4.625 -2.314) 6.168 5.048 ( 2.448 16.274 -2.369) 16.627 5.265 ( -0.480 -14.745 0.228) 14.755 5.777 ( 4.169 15.076 7.481) 17.339 6.164 ( -1.339 4.398 -11.545) 12.426 6.992 ( 9.135 -1.406 19.822) 21.871 7.830 ( -6.513 -9.578 13.661) 17.910 7.927 ( 13.028 21.224 -0.980) 24.923 9.075 ( -6.184 14.323 -3.075) 15.901 10.128 ( -8.164 -1.191 -3.359) 8.908 11.551 ( 8.036 -5.833 -16.513) 19.268 11.828 ( 17.924 -13.495 -4.612) 22.906 13.462 ( -20.351 0.065 10.211) 22.769 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.013 ( 0.000 2.185 0.047) 2.185 3.078 ( 0.000 1.981 2.824) 3.449 3.387 ( 0.000 -6.204 -1.447) 6.371 5.077 ( 0.000 18.819 -3.648) 19.170 5.252 ( 0.000 -18.884 -0.996) 18.910 5.837 ( 0.000 21.992 19.520) 29.406 6.136 ( 0.000 -1.131 -23.040) 23.068 7.096 ( 0.000 -0.751 22.022) 22.035 7.751 ( 0.000 -11.258 12.765) 17.021 8.078 ( 0.000 21.888 -1.347) 21.929 8.998 ( 0.000 16.442 -3.196) 16.749 10.041 ( 0.000 -5.384 -3.346) 6.339 11.598 ( 0.000 -1.738 -18.142) 18.225 12.112 ( 0.000 -21.000 -7.196) 22.199 13.191 ( -0.000 5.794 14.663) 15.766 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.803 ( 8.834 8.834 3.544) 12.987 2.978 ( 5.199 5.199 1.921) 7.599 3.433 ( -3.341 -3.341 -2.212) 5.217 5.077 ( -7.500 -7.500 -0.992) 10.653 5.252 ( 11.869 11.869 -0.798) 16.804 5.818 ( 4.090 4.090 -14.913) 15.995 6.340 ( 7.178 7.178 12.470) 16.079 6.806 ( 4.370 4.370 14.719) 15.964 7.802 ( -9.835 -9.835 15.541) 20.856 8.006 ( 19.240 19.240 -1.620) 27.257 9.215 ( -5.755 -5.755 -2.116) 8.409 10.672 ( 6.742 6.742 -7.988) 12.438 10.854 ( 3.270 3.270 -1.230) 4.785 11.656 ( -1.139 -1.139 -14.631) 14.720 13.681 ( -12.543 -12.543 8.826) 19.813 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.982 ( 7.841 4.709 2.082) 9.380 3.061 ( 3.171 2.102 2.355) 4.475 3.350 ( -4.617 -4.051 -0.839) 6.199 4.951 ( -4.419 -12.783 -0.824) 13.550 5.435 ( 5.346 18.537 -5.928) 20.183 5.950 ( 9.055 3.039 -14.259) 17.162 6.356 ( -4.632 12.110 10.603) 16.749 6.962 ( 9.049 -1.203 18.145) 20.311 7.604 ( -8.900 -11.858 15.249) 21.269 8.351 ( 13.303 19.254 -2.406) 23.526 9.077 ( -6.017 -9.725 -1.600) 11.547 10.448 ( -11.629 27.264 -4.520) 29.983 11.329 ( 23.386 -14.904 -8.799) 29.094 11.644 ( 0.384 -4.423 -13.097) 13.829 13.389 ( -14.327 -7.071 11.937) 19.943 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.063 ( 0.000 2.462 0.747) 2.573 3.119 ( 0.000 2.047 2.786) 3.457 3.278 ( 0.000 -4.436 0.213) 4.441 4.905 ( 0.000 -14.725 -0.725) 14.743 5.484 ( 0.000 20.464 -8.802) 22.277 6.112 ( 0.000 -1.118 -23.066) 23.093 6.242 ( -0.000 17.168 19.599) 26.055 7.080 ( 0.000 -0.759 20.949) 20.963 7.491 ( 0.000 -13.840 13.877) 19.598 8.498 ( 0.000 18.735 -2.781) 18.940 9.010 ( 0.000 -10.974 -1.489) 11.075 10.319 ( 0.000 27.297 -3.928) 27.579 11.564 ( 0.000 -1.686 -19.613) 19.686 11.690 ( 0.000 -19.481 -5.205) 20.165 13.210 ( 0.000 -3.631 14.700) 15.142 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.084 ( 4.781 4.781 2.075) 7.073 3.093 ( 1.207 1.207 2.169) 2.760 3.271 ( -3.725 -3.725 0.525) 5.294 4.741 ( -7.338 -7.338 -0.586) 10.394 5.775 ( 12.368 12.368 -10.869) 20.593 6.000 ( 3.856 3.856 -20.215) 20.938 6.530 ( 2.883 2.883 16.633) 17.126 6.976 ( 3.644 3.644 17.912) 18.639 7.366 ( -10.860 -10.860 16.895) 22.832 8.684 ( 12.446 12.446 -3.659) 17.978 8.863 ( -8.835 -8.835 -0.682) 12.513 10.975 ( 7.886 7.886 -8.979) 14.318 11.086 ( 5.778 5.778 -2.233) 8.471 11.581 ( -2.804 -2.804 -16.927) 17.386 13.201 ( -9.862 -9.862 14.383) 20.034 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.103 ( 0.000 1.286 1.973) 2.355 3.162 ( 0.000 1.912 1.956) 2.735 3.208 ( -0.000 -2.339 1.228) 2.642 4.664 ( 0.000 -8.414 -0.485) 8.428 5.880 ( 0.000 16.989 -15.544) 23.027 6.093 ( 0.000 -0.730 -22.807) 22.818 6.520 ( -0.000 9.440 20.855) 22.892 7.067 ( 0.000 -0.463 19.952) 19.958 7.210 ( 0.000 -12.469 16.694) 20.837 8.767 ( 0.000 -10.094 -0.459) 10.104 8.819 ( 0.000 11.884 -4.320) 12.644 10.899 ( 0.000 25.539 -5.109) 26.045 11.352 ( 0.000 -12.592 -3.370) 13.035 11.521 ( -0.000 -2.807 -21.069) 21.255 13.078 ( 0.000 -7.781 16.383) 18.137 ======================= 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.10e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.116 ( 0.000 0.000 2.475) 2.475 3.184 ( 0.000 0.000 1.573) 1.573 3.184 ( 0.000 0.000 1.573) 1.573 4.576 ( 0.000 0.000 -0.376) 0.376 6.085 ( 0.000 0.000 -22.657) 22.657 6.085 ( 0.000 0.000 -22.657) 22.657 6.612 ( -0.000 -0.000 24.226) 24.226 7.062 ( 0.000 0.000 19.542) 19.542 7.062 ( 0.000 0.000 19.542) 19.542 8.656 ( 0.000 0.000 -0.233) 0.233 8.948 ( 0.000 0.000 -5.097) 5.097 11.216 ( 0.000 0.000 -2.507) 2.507 11.216 ( 0.000 0.000 -2.507) 2.507 11.461 ( 0.000 0.000 -25.516) 25.516 12.978 ( 0.000 0.000 18.036) 18.036 ======================= 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.10e-04 Number of triplets: 189 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.080 ( 13.505 13.505 13.505) 23.392 2.080 ( 13.505 13.505 13.505) 23.392 3.259 ( 9.846 9.846 9.846) 17.054 4.909 ( 3.996 3.996 3.996) 6.921 4.909 ( 3.996 3.996 3.996) 6.921 5.502 ( 15.589 15.589 15.589) 27.000 5.843 ( -5.880 -5.880 -5.880) 10.185 6.700 ( 7.240 7.240 7.240) 12.540 6.700 ( 7.240 7.240 7.240) 12.540 8.717 ( 4.202 4.202 4.202) 7.278 8.717 ( 4.202 4.202 4.202) 7.278 10.230 ( 1.897 1.897 1.897) 3.286 11.406 ( -7.434 -7.434 -7.434) 12.876 11.406 ( -7.434 -7.434 -7.434) 12.876 14.522 ( -4.946 -4.946 -4.946) 8.567 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.340 ( 14.434 9.767 9.767) 19.979 2.438 ( 17.522 8.222 8.222) 21.029 3.410 ( 4.810 2.323 2.323) 5.825 4.876 ( -1.542 6.015 6.015) 8.645 5.103 ( 10.162 1.773 1.773) 10.467 5.732 ( -4.908 -7.643 -7.643) 11.871 5.857 ( 15.219 16.358 16.358) 27.691 6.813 ( 9.970 8.371 8.371) 15.478 6.952 ( 12.440 6.243 6.243) 15.254 8.587 ( -7.232 8.359 8.359) 13.858 8.887 ( 2.346 5.610 5.610) 8.273 10.287 ( 2.193 -0.789 -0.789) 2.461 11.237 ( 0.245 -9.616 -9.616) 13.601 11.394 ( -0.616 -10.766 -10.766) 15.238 14.341 ( -12.750 -1.205 -1.205) 12.863 ======================= 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.10e-04 Number of triplets: 476 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.639 ( 14.127 5.287 5.287) 15.984 2.770 ( 14.349 5.306 5.306) 16.193 3.465 ( 0.875 -2.014 -2.014) 2.979 4.848 ( -1.206 5.908 5.908) 8.441 5.331 ( 11.951 2.213 2.213) 12.354 5.643 ( -3.717 -9.503 -9.503) 13.944 6.029 ( 1.766 14.209 14.209) 20.171 7.121 ( 18.658 9.119 9.119) 22.681 7.199 ( 11.459 7.474 7.474) 15.590 8.428 ( -8.067 7.575 7.575) 13.410 8.830 ( -6.922 9.154 9.154) 14.680 10.269 ( -3.484 -6.024 -6.024) 9.203 11.382 ( -0.536 -10.645 -10.645) 15.064 11.458 ( 19.891 -11.004 -11.004) 25.255 14.008 ( -19.468 2.761 2.761) 19.855 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.888 ( 9.398 2.384 2.384) 9.985 3.004 ( 8.143 3.328 3.328) 9.405 3.458 ( -1.096 -3.906 -3.906) 5.632 4.828 ( -0.709 6.446 6.446) 9.144 5.578 ( 11.821 2.498 2.498) 12.338 5.583 ( -2.073 -11.124 -11.124) 15.867 5.976 ( -6.004 11.857 11.857) 17.811 7.406 ( 8.137 8.609 8.609) 14.644 7.490 ( 15.500 10.652 10.652) 21.615 8.273 ( -6.504 6.628 6.628) 11.409 8.654 ( -8.934 10.124 10.124) 16.876 10.192 ( -3.209 -11.478 -11.478) 16.546 11.373 ( -0.310 -10.557 -10.557) 14.933 11.941 ( 24.345 -11.548 -11.548) 29.315 13.576 ( -21.300 7.680 7.680) 23.910 ======================= 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.10e-04 Number of triplets: 277 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.991 ( 0.000 1.359 1.359) 1.921 3.090 ( 0.000 2.594 2.594) 3.668 3.442 ( 0.000 -4.389 -4.389) 6.207 4.820 ( 0.000 6.823 6.823) 9.649 5.562 ( 0.000 -11.805 -11.805) 16.695 5.782 ( 0.000 -11.074 -11.074) 15.661 5.830 ( -0.000 24.716 24.716) 34.953 7.497 ( 0.000 9.281 9.281) 13.126 7.670 ( 0.000 11.460 11.460) 16.207 8.197 ( 0.000 6.002 6.002) 8.488 8.541 ( 0.000 10.068 10.068) 14.239 10.157 ( 0.000 -13.494 -13.494) 19.083 11.370 ( 0.000 -10.524 -10.524) 14.883 12.260 ( 0.000 -15.566 -15.566) 22.014 13.292 ( -0.000 13.371 13.371) 18.909 ======================= 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.10e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.539 ( 10.268 10.268 10.325) 17.817 2.663 ( 12.334 12.334 3.793) 17.851 3.436 ( 0.287 0.287 -1.346) 1.406 4.977 ( 2.888 2.888 3.944) 5.678 5.132 ( 3.185 3.185 -3.486) 5.696 5.592 ( -5.811 -5.811 -9.175) 12.317 6.242 ( 16.940 16.940 23.050) 33.245 6.907 ( 5.320 5.320 14.478) 16.316 7.209 ( 15.773 15.773 -1.467) 22.355 8.563 ( -6.484 -6.484 21.494) 23.368 9.081 ( 7.219 7.219 -5.956) 11.820 10.319 ( 2.540 2.540 -4.786) 5.984 10.967 ( -8.960 -8.960 -0.491) 12.682 11.316 ( -0.839 -0.839 -20.028) 20.063 14.245 ( -8.359 -8.359 1.764) 11.952 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.759 ( 10.735 6.525 6.931) 14.348 2.901 ( 10.356 6.904 3.380) 12.897 3.425 ( -1.299 -1.921 -2.797) 3.633 4.993 ( -0.076 8.107 -0.250) 8.111 5.235 ( 4.389 -8.757 -1.761) 9.952 5.516 ( -0.427 -3.051 -13.066) 13.424 6.425 ( 1.205 20.580 22.771) 30.717 7.117 ( 14.092 -1.803 17.404) 22.466 7.553 ( 17.356 19.542 -2.687) 26.275 8.398 ( -8.946 -6.841 21.996) 24.711 9.106 ( -3.738 10.759 -5.081) 12.471 10.237 ( -10.250 5.834 -5.843) 13.162 11.056 ( 13.607 -17.002 -5.131) 22.373 11.342 ( 5.605 -3.417 -18.520) 19.649 14.004 ( -14.805 -3.508 5.177) 16.071 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.953 ( 7.457 3.925 3.418) 9.094 3.072 ( 6.169 3.099 3.611) 7.791 3.388 ( -2.139 -2.999 -2.324) 4.356 4.999 ( 0.651 9.588 -1.978) 9.812 5.247 ( -1.957 -15.593 -1.787) 15.816 5.598 ( 6.440 -0.942 -17.663) 18.824 6.376 ( -4.106 21.527 23.239) 31.942 7.393 ( 11.331 -2.574 19.233) 22.471 7.884 ( 13.693 20.231 -2.260) 24.533 8.218 ( -8.064 -7.978 21.706) 24.492 8.988 ( -6.402 13.990 -4.975) 16.169 10.027 ( -8.331 0.046 -6.175) 10.370 11.202 ( 1.659 -6.610 -15.976) 17.369 11.658 ( 20.118 -14.982 -11.137) 27.445 13.677 ( -15.516 1.629 9.299) 18.162 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.034 ( 0.000 2.758 1.942) 3.373 3.141 ( 0.000 2.392 3.192) 3.988 3.358 ( 0.000 -3.750 -1.283) 3.963 5.009 ( 0.000 11.346 -2.442) 11.606 5.219 ( 0.000 -17.962 -2.433) 18.126 5.680 ( 0.000 -2.275 -20.927) 21.050 6.321 ( -0.000 22.607 25.682) 34.215 7.524 ( 0.000 -1.420 19.309) 19.361 8.043 ( 0.000 19.635 -1.726) 19.711 8.117 ( 0.000 -9.433 21.740) 23.698 8.908 ( 0.000 16.384 -5.169) 17.180 9.937 ( 0.000 -3.034 -6.564) 7.232 11.214 ( 0.000 -5.147 -17.337) 18.085 11.907 ( 0.000 -18.269 -12.410) 22.085 13.485 ( 0.000 5.186 12.491) 13.524 ======================= 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.10e-04 Number of triplets: 476 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.905 ( 7.482 7.482 6.137) 12.232 3.031 ( 5.501 5.501 3.218) 8.419 3.384 ( -2.141 -2.141 -2.265) 3.782 5.048 ( -7.120 -7.120 -1.824) 10.233 5.170 ( 6.558 6.558 -6.416) 11.278 5.494 ( 1.260 1.260 -15.807) 15.907 6.753 ( 8.597 8.597 26.038) 28.736 7.129 ( 5.415 5.415 16.010) 17.747 7.947 ( 18.189 18.189 -3.329) 25.938 8.223 ( -9.512 -9.512 23.397) 26.988 9.150 ( -4.544 -4.544 -4.188) 7.670 10.455 ( 4.376 4.376 -12.174) 13.656 10.820 ( 2.313 2.313 -2.055) 3.863 11.300 ( -0.394 -0.394 -18.653) 18.661 13.864 ( -9.857 -9.857 7.715) 15.932 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.044 ( 5.373 4.909 3.895) 8.255 3.125 ( 3.967 1.980 3.725) 5.790 3.334 ( -2.741 -2.296 -0.544) 3.617 4.927 ( -4.185 -12.263 -1.518) 13.046 5.255 ( 1.569 11.668 -10.422) 15.723 5.575 ( 5.521 -0.989 -19.458) 20.251 6.779 ( -2.576 17.224 25.755) 31.090 7.336 ( 11.644 -2.521 18.206) 21.758 8.011 ( -9.343 -10.701 23.342) 27.325 8.282 ( 12.226 17.277 -4.079) 21.555 9.027 ( -5.654 -7.589 -3.239) 10.003 10.307 ( -9.729 23.229 -9.074) 26.769 11.076 ( 9.239 -6.159 -14.060) 17.916 11.381 ( 10.154 -10.709 -11.909) 18.963 13.636 ( -10.903 -5.338 10.623) 16.132 ======================= 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.10e-04 Number of triplets: 506 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.098 ( 0.000 3.338 2.627) 4.248 3.185 ( 0.000 1.856 3.469) 3.935 3.290 ( 0.000 -2.804 1.010) 2.980 4.884 ( 0.000 -14.307 -1.345) 14.370 5.266 ( 0.000 13.108 -11.502) 17.439 5.645 ( 0.000 -1.365 -22.179) 22.221 6.742 ( 0.000 18.471 26.674) 32.445 7.488 ( 0.000 -1.839 18.617) 18.708 7.885 ( 0.000 -12.604 23.954) 27.067 8.417 ( 0.000 16.474 -4.910) 17.190 8.963 ( 0.000 -8.519 -3.034) 9.044 10.196 ( 0.000 24.685 -7.953) 25.934 11.142 ( 0.000 -2.149 -19.366) 19.485 11.542 ( 0.000 -16.701 -8.934) 18.941 13.504 ( 0.000 -2.962 12.519) 12.864 ======================= 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.10e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.143 ( 4.094 4.094 3.654) 6.846 3.154 ( 1.325 1.325 3.587) 4.047 3.289 ( -2.252 -2.252 1.270) 3.429 4.724 ( -7.070 -7.070 -1.042) 10.053 5.463 ( 6.632 6.632 -17.412) 19.777 5.574 ( 1.857 1.857 -20.600) 20.767 7.036 ( 6.273 6.273 28.989) 30.316 7.337 ( 4.450 4.450 16.867) 18.003 7.786 ( -10.709 -10.709 23.776) 28.189 8.578 ( 10.868 10.868 -6.517) 16.694 8.842 ( -8.083 -8.083 -1.396) 11.516 10.690 ( 6.883 6.883 -17.836) 20.320 11.022 ( 5.346 5.346 -3.913) 8.513 11.250 ( -2.382 -2.382 -13.914) 14.316 13.490 ( -7.549 -7.549 12.309) 16.294 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.159 ( 0.000 2.327 3.408) 4.127 3.216 ( 0.000 1.166 3.207) 3.413 3.245 ( -0.000 -1.525 2.351) 2.802 4.650 ( 0.000 -8.124 -0.853) 8.169 5.509 ( -0.000 9.643 -19.133) 21.425 5.623 ( 0.000 -0.733 -22.539) 22.551 7.074 ( -0.000 13.942 29.351) 32.494 7.455 ( 0.000 -1.195 17.668) 17.709 7.625 ( 0.000 -12.236 24.086) 27.016 8.695 ( 0.000 10.024 -7.598) 12.578 8.753 ( 0.000 -9.268 -0.922) 9.314 10.699 ( 0.000 20.073 -15.261) 25.216 11.122 ( 0.000 0.354 -14.745) 14.749 11.256 ( -0.000 -10.489 -5.827) 11.999 13.398 ( 0.000 -6.099 13.425) 14.746 ======================= 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.10e-04 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.184 ( -0.000 -0.000 3.965) 3.965 3.230 ( 0.000 0.000 2.787) 2.787 3.230 ( 0.000 0.000 2.787) 2.787 4.565 ( 0.000 0.000 -0.652) 0.652 5.616 ( 0.000 0.000 -22.625) 22.625 5.616 ( 0.000 0.000 -22.625) 22.625 7.271 ( -0.000 -0.000 37.955) 37.955 7.443 ( 0.000 0.000 17.281) 17.281 7.443 ( 0.000 0.000 17.281) 17.281 8.649 ( 0.000 0.000 -0.404) 0.404 8.802 ( 0.000 0.000 -8.839) 8.839 10.906 ( -0.000 -0.000 -26.672) 26.672 11.144 ( 0.000 0.000 -4.441) 4.441 11.144 ( 0.000 0.000 -4.441) 4.441 13.322 ( 0.000 0.000 14.194) 14.194 ======================= 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.10e-04 Number of triplets: 182 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.768 ( 8.894 8.894 8.894) 15.405 2.768 ( 8.894 8.894 8.894) 15.405 3.407 ( -1.494 -1.494 -1.494) 2.587 5.056 ( 0.420 0.420 0.420) 0.728 5.056 ( 0.420 0.420 0.420) 0.728 5.421 ( -7.535 -7.535 -7.535) 13.051 6.739 ( 21.337 21.337 21.337) 36.957 7.204 ( 8.843 8.843 8.843) 15.317 7.204 ( 8.843 8.843 8.843) 15.317 8.936 ( 2.150 2.150 2.150) 3.723 8.936 ( 2.150 2.150 2.150) 3.723 10.265 ( -0.708 -0.708 -0.708) 1.226 10.956 ( -6.303 -6.303 -6.303) 10.918 10.956 ( -6.303 -6.303 -6.303) 10.918 14.216 ( -4.846 -4.846 -4.846) 8.393 ======================= 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.10e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.909 ( 7.054 7.267 7.267) 12.465 2.989 ( 9.673 5.177 5.177) 12.132 3.376 ( -1.470 -1.911 -1.911) 3.076 4.998 ( -1.936 1.030 1.030) 2.423 5.164 ( 5.892 -5.907 -5.907) 10.223 5.282 ( -6.180 -8.146 -8.146) 13.074 6.923 ( -0.897 23.692 23.692) 33.517 7.485 ( 12.914 6.509 6.509) 15.859 7.486 ( 21.061 8.345 8.345) 24.143 8.757 ( -9.127 7.315 7.315) 13.796 9.029 ( -2.112 0.209 0.209) 2.132 10.165 ( -10.548 0.969 0.969) 10.636 10.946 ( -0.462 -9.369 -9.369) 13.257 10.994 ( 15.827 -10.331 -10.331) 21.540 14.050 ( -10.984 -1.041 -1.041) 11.082 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.036 ( 4.898 4.620 4.620) 8.166 3.147 ( 5.476 3.475 3.475) 7.358 3.349 ( -1.141 -1.437 -1.437) 2.330 4.982 ( 0.303 0.302 0.302) 0.524 5.170 ( -4.545 -8.644 -8.644) 13.043 5.279 ( 4.678 -9.359 -9.359) 14.038 6.864 ( -3.236 22.792 22.792) 32.394 7.718 ( 9.314 6.672 6.672) 13.258 7.886 ( 15.927 8.176 8.176) 19.682 8.579 ( -7.591 7.337 7.337) 12.856 8.929 ( -5.909 1.916 1.916) 6.501 9.935 ( -9.624 0.193 0.193) 9.628 10.938 ( -0.262 -9.488 -9.488) 13.421 11.403 ( 19.295 -13.112 -13.112) 26.761 13.805 ( -11.480 2.761 2.761) 12.126 ======================= 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.10e-04 Number of triplets: 276 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.090 ( 0.000 3.404 3.404) 4.814 3.204 ( 0.000 2.850 2.850) 4.031 3.335 ( 0.000 -0.962 -0.962) 1.360 4.992 ( 0.000 1.029 1.029) 1.456 5.116 ( 0.000 -9.767 -9.767) 13.813 5.331 ( 0.000 -10.807 -10.807) 15.283 6.826 ( 0.000 22.913 22.913) 32.404 7.825 ( -0.000 6.693 6.693) 9.465 8.065 ( 0.000 7.345 7.345) 10.388 8.490 ( 0.000 7.361 7.361) 10.410 8.853 ( 0.000 3.235 3.235) 4.575 9.829 ( 0.000 -0.762 -0.762) 1.077 10.935 ( 0.000 -9.528 -9.528) 13.475 11.627 ( 0.000 -14.433 -14.433) 20.411 13.667 ( 0.000 4.939 4.939) 6.985 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.043 ( 5.714 5.714 7.044) 10.719 3.104 ( 5.653 5.653 3.784) 8.845 3.346 ( -1.072 -1.072 -1.362) 2.037 5.005 ( -6.414 -6.414 -2.333) 9.366 5.038 ( 0.897 0.897 -5.743) 5.881 5.194 ( -1.665 -1.665 -13.230) 13.438 7.308 ( 8.079 8.079 26.320) 28.693 7.446 ( 5.624 5.624 14.779) 16.783 7.864 ( 17.753 17.753 -4.563) 25.518 8.691 ( -8.266 -8.266 19.539) 22.769 9.045 ( -2.768 -2.768 -6.029) 7.188 10.246 ( 0.649 0.649 -5.827) 5.899 10.775 ( 0.987 0.987 -2.250) 2.648 10.937 ( 1.478 1.478 -15.906) 16.043 13.961 ( -7.246 -7.246 1.450) 10.350 ======================= 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.10e-04 Number of triplets: 864 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.135 ( 3.080 4.883 4.776) 7.493 3.205 ( 4.049 2.080 3.934) 6.016 3.329 ( -0.761 -0.565 -0.015) 0.948 4.890 ( -3.892 -11.025 -2.164) 11.890 5.053 ( -0.156 2.454 -8.590) 8.936 5.199 ( 1.640 -1.881 -16.766) 16.951 7.283 ( -3.092 18.128 22.002) 28.676 7.696 ( 11.618 -3.505 16.455) 20.446 8.177 ( 11.878 14.678 -4.761) 19.473 8.512 ( -8.261 -8.641 21.890) 24.941 8.945 ( -5.114 -4.465 -4.719) 8.268 10.127 ( -7.864 15.788 -5.920) 18.605 10.826 ( 0.387 -0.657 -9.629) 9.660 11.152 ( 15.694 -10.727 -10.546) 21.739 13.788 ( -8.363 -4.041 3.857) 10.057 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.166 ( 0.000 3.836 3.836) 5.425 3.255 ( 0.000 2.022 3.181) 3.769 3.316 ( 0.000 -0.921 1.390) 1.667 4.850 ( 0.000 -13.206 -2.016) 13.359 5.048 ( 0.000 4.057 -9.203) 10.058 5.222 ( 0.000 -2.480 -18.702) 18.866 7.248 ( 0.000 18.405 21.821) 28.547 7.831 ( 0.000 -2.424 14.608) 14.808 8.308 ( 0.000 12.660 -5.310) 13.729 8.409 ( 0.000 -9.575 24.087) 25.920 8.886 ( 0.000 -4.921 -4.437) 6.625 10.035 ( 0.000 18.875 -5.350) 19.619 10.826 ( 0.000 -0.956 -11.083) 11.124 11.341 ( 0.000 -12.884 -10.270) 16.476 13.688 ( 0.000 -2.518 5.132) 5.716 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.226 ( 3.159 3.159 4.232) 6.153 3.233 ( 1.412 1.412 3.986) 4.458 3.320 ( -0.601 -0.601 1.562) 1.778 4.700 ( -6.660 -6.660 -1.262) 9.502 5.118 ( 2.440 2.440 -15.646) 16.022 5.177 ( 0.000 0.000 -17.831) 17.831 7.577 ( 5.900 5.900 22.124) 23.645 7.656 ( 4.417 4.417 14.194) 15.508 8.313 ( -8.628 -8.628 25.204) 28.002 8.430 ( 8.321 8.321 -7.214) 13.803 8.806 ( -6.779 -6.779 -2.136) 9.822 10.348 ( 3.697 3.697 -12.411) 13.467 10.933 ( 4.673 4.673 -4.602) 8.053 11.008 ( 1.529 1.529 -10.211) 10.438 13.671 ( -6.144 -6.144 5.044) 10.046 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.236 ( 0.000 2.792 3.972) 4.855 3.285 ( 0.000 0.910 3.307) 3.430 3.300 ( -0.000 -0.639 2.787) 2.860 4.631 ( 0.000 -7.679 -1.016) 7.746 5.135 ( -0.000 3.799 -16.843) 17.266 5.188 ( 0.000 -1.016 -19.597) 19.624 7.580 ( 0.000 13.901 19.232) 23.730 7.782 ( 0.000 -1.935 14.150) 14.281 8.196 ( 0.000 -10.295 29.010) 30.783 8.522 ( 0.000 7.619 -8.919) 11.730 8.729 ( 0.000 -7.732 -1.388) 7.856 10.365 ( 0.000 9.575 -13.109) 16.234 10.915 ( 0.000 8.907 -7.242) 11.479 11.124 ( 0.000 -7.794 -6.771) 10.325 13.596 ( 0.000 -5.232 5.706) 7.742 ======================= 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.10e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.268 ( -0.000 -0.000 4.067) 4.067 3.293 ( 0.000 0.000 3.222) 3.222 3.293 ( 0.000 0.000 3.222) 3.222 4.550 ( 0.000 0.000 -0.755) 0.755 5.178 ( 0.000 0.000 -19.751) 19.751 5.178 ( 0.000 0.000 -19.751) 19.751 7.761 ( -0.000 -0.000 13.757) 13.757 7.761 ( 0.000 0.000 13.757) 13.757 8.052 ( 0.000 0.000 35.185) 35.185 8.602 ( 0.000 0.000 -10.213) 10.213 8.640 ( 0.000 0.000 -0.466) 0.466 10.446 ( 0.000 0.000 -16.090) 16.090 11.042 ( 0.000 0.000 -5.322) 5.322 11.042 ( 0.000 0.000 -5.322) 5.322 13.532 ( 0.000 0.000 6.107) 6.107 ======================= 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.10e-04 Number of triplets: 189 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.179 ( 4.654 4.654 4.654) 8.062 3.179 ( 4.654 4.654 4.654) 8.062 3.329 ( -0.294 -0.294 -0.294) 0.509 4.951 ( -3.630 -3.630 -3.630) 6.288 4.951 ( -3.630 -3.630 -3.630) 6.288 4.969 ( -6.774 -6.774 -6.774) 11.734 7.727 ( 7.711 7.711 7.711) 13.355 7.727 ( 7.711 7.711 7.711) 13.355 7.826 ( 12.381 12.381 12.381) 21.445 8.908 ( -2.885 -2.885 -2.885) 4.997 8.908 ( -2.885 -2.885 -2.885) 4.997 10.249 ( 1.712 1.712 1.712) 2.966 10.732 ( -0.845 -0.845 -0.845) 1.463 10.732 ( -0.845 -0.845 -0.845) 1.463 13.923 ( -4.894 -4.894 -4.894) 8.477 ======================= 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.10e-04 Number of triplets: 468 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.231 ( 1.620 4.424 4.424) 6.462 3.276 ( 3.332 2.726 2.726) 5.096 3.333 ( 0.533 0.463 0.463) 0.844 4.824 ( -4.404 -5.699 -5.699) 9.185 4.918 ( -1.494 -5.425 -5.425) 7.815 4.922 ( -0.535 -7.708 -7.708) 10.913 7.664 ( -3.222 15.486 15.486) 22.137 7.949 ( 8.580 4.580 4.580) 10.750 8.127 ( 11.258 3.473 3.473) 12.283 8.768 ( -6.028 1.028 1.028) 6.201 8.848 ( -4.463 -4.029 -4.029) 7.237 10.170 ( -6.268 9.439 9.439) 14.747 10.702 ( -1.016 -1.319 -1.319) 2.124 10.959 ( 15.206 -7.958 -7.958) 18.918 13.793 ( -6.501 -2.964 -2.964) 7.735 ======================= 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.10e-04 Number of triplets: 277 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.247 ( 0.000 3.848 3.848) 5.442 3.311 ( 0.000 2.251 2.251) 3.184 3.341 ( 0.000 0.900 0.900) 1.273 4.779 ( 0.000 -6.561 -6.561) 9.278 4.902 ( 0.000 -4.973 -4.973) 7.032 4.918 ( 0.000 -9.043 -9.043) 12.789 7.629 ( 0.000 15.681 15.681) 22.176 8.045 ( -0.000 4.105 4.105) 5.805 8.246 ( 0.000 1.580 1.580) 2.235 8.697 ( 0.000 1.892 1.892) 2.676 8.795 ( 0.000 -3.963 -3.963) 5.604 10.096 ( 0.000 11.423 11.423) 16.154 10.692 ( 0.000 -1.672 -1.672) 2.364 11.140 ( 0.000 -8.973 -8.973) 12.689 13.715 ( 0.000 -2.092 -2.092) 2.959 ======================= 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.10e-04 Number of triplets: 474 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.308 ( 2.208 2.208 3.513) 4.701 3.310 ( 1.372 1.372 3.356) 3.876 3.348 ( 0.763 0.763 1.166) 1.588 4.674 ( -6.119 -6.119 -1.250) 8.744 4.843 ( -1.241 -1.241 -10.906) 11.047 4.860 ( -1.743 -1.743 -12.813) 13.048 7.914 ( 3.734 3.734 10.923) 12.132 7.933 ( 3.096 3.096 13.041) 13.756 8.275 ( 6.756 6.756 -6.602) 11.614 8.688 ( -5.027 -5.027 7.902) 10.629 8.755 ( -4.959 -4.959 -2.725) 7.524 10.306 ( 1.511 1.511 9.311) 9.553 10.828 ( 4.133 4.133 -7.077) 9.178 10.842 ( 3.880 3.880 -4.072) 6.833 13.698 ( -5.235 -5.235 -2.081) 7.690 ======================= 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.10e-04 Number of triplets: 504 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.314 ( 0.000 2.531 3.436) 4.268 3.345 ( 0.000 1.042 2.435) 2.648 3.352 ( -0.000 0.246 2.173) 2.187 4.610 ( 0.000 -7.043 -0.977) 7.111 4.838 ( -0.000 -1.425 -14.353) 14.424 4.838 ( 0.000 -1.590 -11.883) 11.989 7.892 ( 0.000 9.804 12.021) 15.512 8.029 ( 0.000 -2.378 9.853) 10.136 8.348 ( 0.000 5.120 -7.688) 9.237 8.627 ( 0.000 -4.827 10.207) 11.291 8.696 ( 0.000 -5.379 -1.777) 5.665 10.307 ( 0.000 5.267 8.744) 10.208 10.794 ( 0.000 10.431 -4.710) 11.445 10.991 ( 0.000 -5.280 -5.973) 7.972 13.634 ( 0.000 -4.743 -1.684) 5.033 ======================= 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.10e-04 Number of triplets: 172 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.343 ( -0.000 -0.000 3.064) 3.064 3.354 ( 0.000 0.000 2.596) 2.596 3.354 ( 0.000 0.000 2.596) 2.596 4.536 ( 0.000 0.000 -0.656) 0.656 4.824 ( 0.000 0.000 -14.488) 14.488 4.824 ( 0.000 0.000 -14.488) 14.488 8.001 ( 0.000 0.000 9.652) 9.652 8.001 ( 0.000 0.000 9.652) 9.652 8.403 ( 0.000 0.000 -8.834) 8.834 8.569 ( 0.000 0.000 13.079) 13.079 8.631 ( 0.000 0.000 -0.404) 0.404 10.346 ( 0.000 0.000 7.418) 7.418 10.935 ( 0.000 0.000 -4.826) 4.826 10.935 ( 0.000 0.000 -4.826) 4.826 13.575 ( 0.000 0.000 -1.488) 1.488 ======================= 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.10e-04 Number of triplets: 182 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.363 ( 1.595 1.595 1.595) 2.763 3.363 ( 1.595 1.595 1.595) 2.763 3.367 ( 1.054 1.054 1.054) 1.826 4.637 ( -3.772 -3.772 -3.772) 6.533 4.677 ( -4.442 -4.442 -4.442) 7.693 4.677 ( -4.442 -4.442 -4.442) 7.693 8.095 ( 3.924 3.924 3.924) 6.796 8.095 ( 3.924 3.924 3.924) 6.796 8.200 ( 1.239 1.239 1.239) 2.146 8.700 ( -2.786 -2.786 -2.786) 4.826 8.700 ( -2.786 -2.786 -2.786) 4.826 10.556 ( 6.952 6.952 6.952) 12.041 10.774 ( 1.338 1.338 1.338) 2.317 10.774 ( 1.338 1.338 1.338) 2.317 13.617 ( -4.696 -4.696 -4.696) 8.134 ======================= 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.10e-04 Number of triplets: 276 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.370 ( 0.000 1.931 1.931) 2.730 3.382 ( 0.000 1.184 1.184) 1.674 3.384 ( 0.000 0.978 0.978) 1.383 4.583 ( 0.000 -3.086 -3.086) 4.364 4.623 ( 0.000 -5.018 -5.018) 7.097 4.657 ( 0.000 -5.847 -5.847) 8.269 8.083 ( 0.000 6.571 6.571) 9.293 8.166 ( 0.000 1.890 1.890) 2.673 8.237 ( 0.000 -1.212 -1.212) 1.714 8.657 ( 0.000 -2.109 -2.109) 2.983 8.672 ( 0.000 -1.926 -1.926) 2.724 10.583 ( 0.000 10.562 10.562) 14.937 10.749 ( 0.000 2.931 2.931) 4.145 10.891 ( 0.000 -3.529 -3.529) 4.990 13.559 ( 0.000 -4.518 -4.518) 6.390 ======================= 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.10e-04 Number of triplets: 168 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.390 ( -0.000 -0.000 1.565) 1.565 3.394 ( 0.000 0.000 1.259) 1.259 3.394 ( 0.000 0.000 1.259) 1.259 4.525 ( 0.000 0.000 -0.380) 0.380 4.598 ( 0.000 0.000 -7.506) 7.506 4.598 ( 0.000 0.000 -7.506) 7.506 8.152 ( 0.000 0.000 5.024) 5.024 8.152 ( 0.000 0.000 5.024) 5.024 8.257 ( 0.000 0.000 -5.088) 5.088 8.625 ( 0.000 0.000 -0.234) 0.234 8.649 ( 0.000 0.000 -1.416) 1.416 10.637 ( 0.000 0.000 15.879) 15.879 10.854 ( 0.000 0.000 -2.902) 2.902 10.854 ( 0.000 0.000 -2.902) 2.902 13.502 ( 0.000 0.000 -4.476) 4.476 ======================= 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.10e-04 Number of triplets: 44 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.406 ( 0.000 0.000 0.000) 0.000 3.406 ( 0.000 0.000 0.000) 0.000 3.406 ( 0.000 0.000 0.000) 0.000 4.521 ( 0.000 0.000 0.000) 0.000 4.521 ( 0.000 0.000 0.000) 0.000 4.521 ( 0.000 0.000 0.000) 0.000 8.204 ( 0.000 0.000 0.000) 0.000 8.204 ( 0.000 0.000 0.000) 0.000 8.204 ( 0.000 0.000 0.000) 0.000 8.622 ( 0.000 0.000 0.000) 0.000 8.622 ( 0.000 0.000 0.000) 0.000 10.823 ( 0.000 0.000 0.000) 0.000 10.823 ( 0.000 0.000 0.000) 0.000 10.823 ( 0.000 0.000 0.000) 0.000 13.446 ( 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 2634.317 2634.317 2634.317 -0.000 -0.000 0.000 3/25920 20.0 668.857 668.857 668.857 -0.000 -0.000 0.000 3/25920 30.0 247.357 247.357 247.357 -0.000 -0.000 0.000 3/25920 40.0 114.737 114.737 114.737 -0.000 -0.000 0.000 3/25920 50.0 63.973 63.973 63.973 -0.000 -0.000 0.000 3/25920 60.0 41.241 41.241 41.241 -0.000 -0.000 0.000 3/25920 70.0 29.574 29.574 29.574 -0.000 -0.000 0.000 3/25920 80.0 22.883 22.883 22.883 -0.000 -0.000 0.000 3/25920 90.0 18.689 18.689 18.689 -0.000 -0.000 0.000 3/25920 100.0 15.865 15.865 15.865 -0.000 -0.000 0.000 3/25920 110.0 13.850 13.850 13.850 -0.000 -0.000 0.000 3/25920 120.0 12.342 12.342 12.342 -0.000 -0.000 0.000 3/25920 130.0 11.170 11.170 11.170 -0.000 -0.000 0.000 3/25920 140.0 10.229 10.229 10.229 -0.000 -0.000 0.000 3/25920 150.0 9.454 9.454 9.454 -0.000 -0.000 0.000 3/25920 160.0 8.803 8.803 8.803 -0.000 -0.000 0.000 3/25920 170.0 8.246 8.246 8.246 -0.000 -0.000 0.000 3/25920 180.0 7.763 7.763 7.763 -0.000 -0.000 0.000 3/25920 190.0 7.338 7.338 7.338 -0.000 -0.000 0.000 3/25920 200.0 6.962 6.962 6.962 -0.000 -0.000 0.000 3/25920 210.0 6.625 6.625 6.625 -0.000 -0.000 0.000 3/25920 220.0 6.321 6.321 6.321 -0.000 -0.000 0.000 3/25920 230.0 6.045 6.045 6.045 -0.000 -0.000 0.000 3/25920 240.0 5.794 5.794 5.794 -0.000 -0.000 0.000 3/25920 250.0 5.563 5.563 5.563 -0.000 -0.000 0.000 3/25920 260.0 5.351 5.351 5.351 -0.000 -0.000 0.000 3/25920 270.0 5.155 5.155 5.155 -0.000 -0.000 0.000 3/25920 280.0 4.974 4.974 4.974 -0.000 -0.000 0.000 3/25920 290.0 4.805 4.805 4.805 -0.000 -0.000 0.000 3/25920 300.0 4.647 4.647 4.647 -0.000 -0.000 0.000 3/25920 310.0 4.499 4.499 4.499 -0.000 -0.000 0.000 3/25920 320.0 4.361 4.361 4.361 -0.000 -0.000 0.000 3/25920 330.0 4.231 4.231 4.231 -0.000 -0.000 0.000 3/25920 340.0 4.109 4.109 4.109 -0.000 -0.000 0.000 3/25920 350.0 3.994 3.994 3.994 -0.000 -0.000 0.000 3/25920 360.0 3.885 3.885 3.885 -0.000 -0.000 0.000 3/25920 370.0 3.782 3.782 3.782 -0.000 -0.000 0.000 3/25920 380.0 3.684 3.684 3.684 -0.000 -0.000 0.000 3/25920 390.0 3.591 3.591 3.591 -0.000 -0.000 0.000 3/25920 400.0 3.503 3.503 3.503 -0.000 -0.000 0.000 3/25920 410.0 3.419 3.419 3.419 -0.000 -0.000 0.000 3/25920 420.0 3.339 3.339 3.339 -0.000 -0.000 0.000 3/25920 430.0 3.262 3.262 3.262 -0.000 -0.000 0.000 3/25920 440.0 3.189 3.189 3.189 -0.000 -0.000 0.000 3/25920 450.0 3.120 3.120 3.120 -0.000 -0.000 0.000 3/25920 460.0 3.053 3.053 3.053 -0.000 -0.000 0.000 3/25920 470.0 2.989 2.989 2.989 -0.000 -0.000 0.000 3/25920 480.0 2.928 2.928 2.928 -0.000 -0.000 0.000 3/25920 490.0 2.869 2.869 2.869 -0.000 -0.000 0.000 3/25920 500.0 2.812 2.812 2.812 -0.000 -0.000 0.000 3/25920 510.0 2.758 2.758 2.758 -0.000 -0.000 0.000 3/25920 520.0 2.706 2.706 2.706 -0.000 -0.000 0.000 3/25920 530.0 2.655 2.655 2.655 -0.000 -0.000 0.000 3/25920 540.0 2.607 2.607 2.607 -0.000 -0.000 0.000 3/25920 550.0 2.560 2.560 2.560 -0.000 -0.000 0.000 3/25920 560.0 2.515 2.515 2.515 -0.000 -0.000 0.000 3/25920 570.0 2.472 2.472 2.472 -0.000 -0.000 0.000 3/25920 580.0 2.429 2.429 2.429 -0.000 -0.000 0.000 3/25920 590.0 2.389 2.389 2.389 -0.000 -0.000 0.000 3/25920 600.0 2.349 2.349 2.349 -0.000 -0.000 0.000 3/25920 610.0 2.311 2.311 2.311 -0.000 -0.000 0.000 3/25920 620.0 2.275 2.275 2.275 -0.000 -0.000 0.000 3/25920 630.0 2.239 2.239 2.239 -0.000 -0.000 0.000 3/25920 640.0 2.204 2.204 2.204 -0.000 -0.000 0.000 3/25920 650.0 2.171 2.171 2.171 -0.000 -0.000 0.000 3/25920 660.0 2.138 2.138 2.138 -0.000 -0.000 0.000 3/25920 670.0 2.107 2.107 2.107 -0.000 -0.000 0.000 3/25920 680.0 2.076 2.076 2.076 -0.000 -0.000 0.000 3/25920 690.0 2.046 2.046 2.046 -0.000 -0.000 0.000 3/25920 700.0 2.017 2.017 2.017 -0.000 -0.000 0.000 3/25920 710.0 1.989 1.989 1.989 -0.000 -0.000 0.000 3/25920 720.0 1.962 1.962 1.962 -0.000 -0.000 0.000 3/25920 730.0 1.935 1.935 1.935 -0.000 -0.000 0.000 3/25920 740.0 1.909 1.909 1.909 -0.000 -0.000 0.000 3/25920 750.0 1.884 1.884 1.884 -0.000 -0.000 0.000 3/25920 760.0 1.860 1.860 1.860 -0.000 -0.000 0.000 3/25920 770.0 1.836 1.836 1.836 -0.000 -0.000 0.000 3/25920 780.0 1.812 1.812 1.812 -0.000 -0.000 0.000 3/25920 790.0 1.790 1.790 1.790 -0.000 -0.000 0.000 3/25920 800.0 1.767 1.767 1.767 -0.000 -0.000 0.000 3/25920 810.0 1.746 1.746 1.746 -0.000 -0.000 0.000 3/25920 820.0 1.725 1.725 1.725 -0.000 -0.000 0.000 3/25920 830.0 1.704 1.704 1.704 -0.000 -0.000 0.000 3/25920 840.0 1.684 1.684 1.684 -0.000 -0.000 0.000 3/25920 850.0 1.664 1.664 1.664 -0.000 -0.000 0.000 3/25920 860.0 1.645 1.645 1.645 -0.000 -0.000 0.000 3/25920 870.0 1.626 1.626 1.626 -0.000 -0.000 0.000 3/25920 880.0 1.608 1.608 1.608 -0.000 -0.000 0.000 3/25920 890.0 1.590 1.590 1.590 -0.000 -0.000 0.000 3/25920 900.0 1.573 1.573 1.573 -0.000 -0.000 0.000 3/25920 910.0 1.555 1.555 1.555 -0.000 -0.000 0.000 3/25920 920.0 1.539 1.539 1.539 -0.000 -0.000 0.000 3/25920 930.0 1.522 1.522 1.522 -0.000 -0.000 0.000 3/25920 940.0 1.506 1.506 1.506 -0.000 -0.000 0.000 3/25920 950.0 1.490 1.490 1.490 -0.000 -0.000 0.000 3/25920 960.0 1.475 1.475 1.475 -0.000 -0.000 0.000 3/25920 970.0 1.460 1.460 1.460 -0.000 -0.000 0.000 3/25920 980.0 1.445 1.445 1.445 -0.000 -0.000 0.000 3/25920 990.0 1.431 1.431 1.431 -0.000 -0.000 0.000 3/25920 1000.0 1.416 1.416 1.416 -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 22:53:53]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|