----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 18:31:48]------------------------- 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: [2 2 2] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fm-3m (225) Number of symmetry operations in supercell: 1536 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 2.684118995000000 2.684118995000000 b 2.684118995000000 0.000000000000000 2.684118995000000 c 2.684118995000000 2.684118995000000 0.000000000000000 Atomic positions (fractional): *1 F 0.50000000000000 0.50000000000000 0.50000000000000 18.998 *2 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 -------------------------------- unit cell --------------------------------- Lattice vectors: a 5.368237990000000 0.000000000000000 0.000000000000000 b 0.000000000000000 5.368237990000000 0.000000000000000 c 0.000000000000000 0.000000000000000 5.368237990000000 Atomic positions (fractional): *1 F 0.50000000000000 0.00000000000000 0.00000000000000 18.998 > 1 2 F 0.50000000000000 0.50000000000000 0.50000000000000 18.998 > 1 3 F 0.00000000000000 0.00000000000000 0.50000000000000 18.998 > 1 4 F 0.00000000000000 0.50000000000000 0.00000000000000 18.998 > 1 *5 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 > 2 6 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 > 2 7 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 > 2 8 K 0.50000000000000 0.50000000000000 0.00000000000000 39.098 > 2 -------------------------------- super cell -------------------------------- Lattice vectors: a 10.736475980000000 0.000000000000000 0.000000000000000 b 0.000000000000000 10.736475980000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.736475980000000 Atomic positions (fractional): *1 F 0.25000000000000 0.00000000000000 0.00000000000000 18.998 > 1 2 F 0.75000000000000 0.00000000000000 0.00000000000000 18.998 > 1 3 F 0.25000000000000 0.50000000000000 0.00000000000000 18.998 > 1 4 F 0.75000000000000 0.50000000000000 0.00000000000000 18.998 > 1 5 F 0.25000000000000 0.00000000000000 0.50000000000000 18.998 > 1 6 F 0.75000000000000 0.00000000000000 0.50000000000000 18.998 > 1 7 F 0.25000000000000 0.50000000000000 0.50000000000000 18.998 > 1 8 F 0.75000000000000 0.50000000000000 0.50000000000000 18.998 > 1 9 F 0.25000000000000 0.25000000000000 0.25000000000000 18.998 > 1 10 F 0.75000000000000 0.25000000000000 0.25000000000000 18.998 > 1 11 F 0.25000000000000 0.75000000000000 0.25000000000000 18.998 > 1 12 F 0.75000000000000 0.75000000000000 0.25000000000000 18.998 > 1 13 F 0.25000000000000 0.25000000000000 0.75000000000000 18.998 > 1 14 F 0.75000000000000 0.25000000000000 0.75000000000000 18.998 > 1 15 F 0.25000000000000 0.75000000000000 0.75000000000000 18.998 > 1 16 F 0.75000000000000 0.75000000000000 0.75000000000000 18.998 > 1 17 F 0.00000000000000 0.00000000000000 0.25000000000000 18.998 > 1 18 F 0.50000000000000 0.00000000000000 0.25000000000000 18.998 > 1 19 F 0.00000000000000 0.50000000000000 0.25000000000000 18.998 > 1 20 F 0.50000000000000 0.50000000000000 0.25000000000000 18.998 > 1 21 F 0.00000000000000 0.00000000000000 0.75000000000000 18.998 > 1 22 F 0.50000000000000 0.00000000000000 0.75000000000000 18.998 > 1 23 F 0.00000000000000 0.50000000000000 0.75000000000000 18.998 > 1 24 F 0.50000000000000 0.50000000000000 0.75000000000000 18.998 > 1 25 F 0.00000000000000 0.25000000000000 0.00000000000000 18.998 > 1 26 F 0.50000000000000 0.25000000000000 0.00000000000000 18.998 > 1 27 F 0.00000000000000 0.75000000000000 0.00000000000000 18.998 > 1 28 F 0.50000000000000 0.75000000000000 0.00000000000000 18.998 > 1 29 F 0.00000000000000 0.25000000000000 0.50000000000000 18.998 > 1 30 F 0.50000000000000 0.25000000000000 0.50000000000000 18.998 > 1 31 F 0.00000000000000 0.75000000000000 0.50000000000000 18.998 > 1 32 F 0.50000000000000 0.75000000000000 0.50000000000000 18.998 > 1 *33 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 > 2 34 K 0.50000000000000 0.00000000000000 0.00000000000000 39.098 > 2 35 K 0.00000000000000 0.50000000000000 0.00000000000000 39.098 > 2 36 K 0.50000000000000 0.50000000000000 0.00000000000000 39.098 > 2 37 K 0.00000000000000 0.00000000000000 0.50000000000000 39.098 > 2 38 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 > 2 39 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 > 2 40 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 2 41 K 0.00000000000000 0.25000000000000 0.25000000000000 39.098 > 2 42 K 0.50000000000000 0.25000000000000 0.25000000000000 39.098 > 2 43 K 0.00000000000000 0.75000000000000 0.25000000000000 39.098 > 2 44 K 0.50000000000000 0.75000000000000 0.25000000000000 39.098 > 2 45 K 0.00000000000000 0.25000000000000 0.75000000000000 39.098 > 2 46 K 0.50000000000000 0.25000000000000 0.75000000000000 39.098 > 2 47 K 0.00000000000000 0.75000000000000 0.75000000000000 39.098 > 2 48 K 0.50000000000000 0.75000000000000 0.75000000000000 39.098 > 2 49 K 0.25000000000000 0.00000000000000 0.25000000000000 39.098 > 2 50 K 0.75000000000000 0.00000000000000 0.25000000000000 39.098 > 2 51 K 0.25000000000000 0.50000000000000 0.25000000000000 39.098 > 2 52 K 0.75000000000000 0.50000000000000 0.25000000000000 39.098 > 2 53 K 0.25000000000000 0.00000000000000 0.75000000000000 39.098 > 2 54 K 0.75000000000000 0.00000000000000 0.75000000000000 39.098 > 2 55 K 0.25000000000000 0.50000000000000 0.75000000000000 39.098 > 2 56 K 0.75000000000000 0.50000000000000 0.75000000000000 39.098 > 2 57 K 0.25000000000000 0.25000000000000 0.00000000000000 39.098 > 2 58 K 0.75000000000000 0.25000000000000 0.00000000000000 39.098 > 2 59 K 0.25000000000000 0.75000000000000 0.00000000000000 39.098 > 2 60 K 0.75000000000000 0.75000000000000 0.00000000000000 39.098 > 2 61 K 0.25000000000000 0.25000000000000 0.50000000000000 39.098 > 2 62 K 0.75000000000000 0.25000000000000 0.50000000000000 39.098 > 2 63 K 0.25000000000000 0.75000000000000 0.50000000000000 39.098 > 2 64 K 0.75000000000000 0.75000000000000 0.50000000000000 39.098 > 2 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 2.0006396 0.0000000 0.0000000 0.0000000 2.0006396 0.0000000 0.0000000 0.0000000 2.0006396 -------------------------- Born effective charges -------------------------- 1 F -1.1410513 0.0000000 0.0000000 0.0000000 -1.1410513 0.0000000 0.0000000 0.0000000 -1.1410513 2 K 1.1410513 0.0000000 0.0000000 0.0000000 1.1410513 0.0000000 0.0000000 0.0000000 1.1410513 ---------------------------------------------------------------------------- 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: 192/192 Permutation basis: 852/852 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 33 Number of blocks in projector: 33 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 28 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 3 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 2 Use standard eigh solver. Tree of FC basis block matrices: - (33, 31), data: False |-- (2, 2), data: True |-- (3, 3), data: True |-- (28, 26), data: True ----- Solver_atoms: 1 -- 64 / 64 Time (Solver_compr_matrix_reshape): 0.001 Solver_block: 80 / 80 - Time: 0.006 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.007 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) Permutation basis: 192/192 Permutation basis: 852/852 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 33 Number of blocks in projector: 33 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 28 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 3 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 2 Use standard eigh solver. Tree of FC basis block matrices: - (33, 31), data: False |-- (2, 2), data: True |-- (3, 3), data: True |-- (28, 26), 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 18:31:51]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:31:51]------------------------- 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): [2 2 2] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fm-3m (225) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 2.684118995000000 2.684118995000000 b 2.684118995000000 0.000000000000000 2.684118995000000 c 2.684118995000000 2.684118995000000 0.000000000000000 Atomic positions (fractional): 1 F 0.50000000000000 0.50000000000000 0.50000000000000 18.998 2 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 -------------------------------- supercell --------------------------------- Lattice vectors: a 10.736475980000000 0.000000000000000 0.000000000000000 b 0.000000000000000 10.736475980000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.736475980000000 Atomic positions (fractional): 1 F 0.25000000000000 0.00000000000000 0.00000000000000 18.998 > 1 2 F 0.75000000000000 0.00000000000000 0.00000000000000 18.998 > 1 3 F 0.25000000000000 0.50000000000000 0.00000000000000 18.998 > 1 4 F 0.75000000000000 0.50000000000000 0.00000000000000 18.998 > 1 5 F 0.25000000000000 0.00000000000000 0.50000000000000 18.998 > 1 6 F 0.75000000000000 0.00000000000000 0.50000000000000 18.998 > 1 7 F 0.25000000000000 0.50000000000000 0.50000000000000 18.998 > 1 8 F 0.75000000000000 0.50000000000000 0.50000000000000 18.998 > 1 9 F 0.25000000000000 0.25000000000000 0.25000000000000 18.998 > 1 10 F 0.75000000000000 0.25000000000000 0.25000000000000 18.998 > 1 11 F 0.25000000000000 0.75000000000000 0.25000000000000 18.998 > 1 12 F 0.75000000000000 0.75000000000000 0.25000000000000 18.998 > 1 13 F 0.25000000000000 0.25000000000000 0.75000000000000 18.998 > 1 14 F 0.75000000000000 0.25000000000000 0.75000000000000 18.998 > 1 15 F 0.25000000000000 0.75000000000000 0.75000000000000 18.998 > 1 16 F 0.75000000000000 0.75000000000000 0.75000000000000 18.998 > 1 17 F 0.00000000000000 0.00000000000000 0.25000000000000 18.998 > 1 18 F 0.50000000000000 0.00000000000000 0.25000000000000 18.998 > 1 19 F 0.00000000000000 0.50000000000000 0.25000000000000 18.998 > 1 20 F 0.50000000000000 0.50000000000000 0.25000000000000 18.998 > 1 21 F 0.00000000000000 0.00000000000000 0.75000000000000 18.998 > 1 22 F 0.50000000000000 0.00000000000000 0.75000000000000 18.998 > 1 23 F 0.00000000000000 0.50000000000000 0.75000000000000 18.998 > 1 24 F 0.50000000000000 0.50000000000000 0.75000000000000 18.998 > 1 25 F 0.00000000000000 0.25000000000000 0.00000000000000 18.998 > 1 26 F 0.50000000000000 0.25000000000000 0.00000000000000 18.998 > 1 27 F 0.00000000000000 0.75000000000000 0.00000000000000 18.998 > 1 28 F 0.50000000000000 0.75000000000000 0.00000000000000 18.998 > 1 29 F 0.00000000000000 0.25000000000000 0.50000000000000 18.998 > 1 30 F 0.50000000000000 0.25000000000000 0.50000000000000 18.998 > 1 31 F 0.00000000000000 0.75000000000000 0.50000000000000 18.998 > 1 32 F 0.50000000000000 0.75000000000000 0.50000000000000 18.998 > 1 33 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 > 33 34 K 0.50000000000000 0.00000000000000 0.00000000000000 39.098 > 33 35 K 0.00000000000000 0.50000000000000 0.00000000000000 39.098 > 33 36 K 0.50000000000000 0.50000000000000 0.00000000000000 39.098 > 33 37 K 0.00000000000000 0.00000000000000 0.50000000000000 39.098 > 33 38 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 > 33 39 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 > 33 40 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 33 41 K 0.00000000000000 0.25000000000000 0.25000000000000 39.098 > 33 42 K 0.50000000000000 0.25000000000000 0.25000000000000 39.098 > 33 43 K 0.00000000000000 0.75000000000000 0.25000000000000 39.098 > 33 44 K 0.50000000000000 0.75000000000000 0.25000000000000 39.098 > 33 45 K 0.00000000000000 0.25000000000000 0.75000000000000 39.098 > 33 46 K 0.50000000000000 0.25000000000000 0.75000000000000 39.098 > 33 47 K 0.00000000000000 0.75000000000000 0.75000000000000 39.098 > 33 48 K 0.50000000000000 0.75000000000000 0.75000000000000 39.098 > 33 49 K 0.25000000000000 0.00000000000000 0.25000000000000 39.098 > 33 50 K 0.75000000000000 0.00000000000000 0.25000000000000 39.098 > 33 51 K 0.25000000000000 0.50000000000000 0.25000000000000 39.098 > 33 52 K 0.75000000000000 0.50000000000000 0.25000000000000 39.098 > 33 53 K 0.25000000000000 0.00000000000000 0.75000000000000 39.098 > 33 54 K 0.75000000000000 0.00000000000000 0.75000000000000 39.098 > 33 55 K 0.25000000000000 0.50000000000000 0.75000000000000 39.098 > 33 56 K 0.75000000000000 0.50000000000000 0.75000000000000 39.098 > 33 57 K 0.25000000000000 0.25000000000000 0.00000000000000 39.098 > 33 58 K 0.75000000000000 0.25000000000000 0.00000000000000 39.098 > 33 59 K 0.25000000000000 0.75000000000000 0.00000000000000 39.098 > 33 60 K 0.75000000000000 0.75000000000000 0.00000000000000 39.098 > 33 61 K 0.25000000000000 0.25000000000000 0.50000000000000 39.098 > 33 62 K 0.75000000000000 0.25000000000000 0.50000000000000 39.098 > 33 63 K 0.25000000000000 0.75000000000000 0.50000000000000 39.098 > 33 64 K 0.75000000000000 0.75000000000000 0.50000000000000 39.098 > 33 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 2.0006396 0.0000000 0.0000000 0.0000000 2.0006396 0.0000000 0.0000000 0.0000000 2.0006396 -------------------------- Born effective charges -------------------------- 1 F -1.1410513 0.0000000 0.0000000 0.0000000 -1.1410513 0.0000000 0.0000000 0.0000000 -1.1410513 2 K 1.1410513 0.0000000 0.0000000 0.0000000 1.1410513 0.0000000 0.0000000 0.0000000 1.1410513 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 33, 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.00000001 (xxy) -0.00000001 (xxy) -0.00000001 (xyx) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) 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 18:31:52]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 18:31:53]------------------------- 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): [2 2 2] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fm-3m (225) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 2.684118995000000 2.684118995000000 b 2.684118995000000 0.000000000000000 2.684118995000000 c 2.684118995000000 2.684118995000000 0.000000000000000 Atomic positions (fractional): 1 F 0.50000000000000 0.50000000000000 0.50000000000000 18.998 2 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 -------------------------------- supercell --------------------------------- Lattice vectors: a 10.736475980000000 0.000000000000000 0.000000000000000 b 0.000000000000000 10.736475980000000 0.000000000000000 c 0.000000000000000 0.000000000000000 10.736475980000000 Atomic positions (fractional): 1 F 0.25000000000000 0.00000000000000 0.00000000000000 18.998 > 1 2 F 0.75000000000000 0.00000000000000 0.00000000000000 18.998 > 1 3 F 0.25000000000000 0.50000000000000 0.00000000000000 18.998 > 1 4 F 0.75000000000000 0.50000000000000 0.00000000000000 18.998 > 1 5 F 0.25000000000000 0.00000000000000 0.50000000000000 18.998 > 1 6 F 0.75000000000000 0.00000000000000 0.50000000000000 18.998 > 1 7 F 0.25000000000000 0.50000000000000 0.50000000000000 18.998 > 1 8 F 0.75000000000000 0.50000000000000 0.50000000000000 18.998 > 1 9 F 0.25000000000000 0.25000000000000 0.25000000000000 18.998 > 1 10 F 0.75000000000000 0.25000000000000 0.25000000000000 18.998 > 1 11 F 0.25000000000000 0.75000000000000 0.25000000000000 18.998 > 1 12 F 0.75000000000000 0.75000000000000 0.25000000000000 18.998 > 1 13 F 0.25000000000000 0.25000000000000 0.75000000000000 18.998 > 1 14 F 0.75000000000000 0.25000000000000 0.75000000000000 18.998 > 1 15 F 0.25000000000000 0.75000000000000 0.75000000000000 18.998 > 1 16 F 0.75000000000000 0.75000000000000 0.75000000000000 18.998 > 1 17 F 0.00000000000000 0.00000000000000 0.25000000000000 18.998 > 1 18 F 0.50000000000000 0.00000000000000 0.25000000000000 18.998 > 1 19 F 0.00000000000000 0.50000000000000 0.25000000000000 18.998 > 1 20 F 0.50000000000000 0.50000000000000 0.25000000000000 18.998 > 1 21 F 0.00000000000000 0.00000000000000 0.75000000000000 18.998 > 1 22 F 0.50000000000000 0.00000000000000 0.75000000000000 18.998 > 1 23 F 0.00000000000000 0.50000000000000 0.75000000000000 18.998 > 1 24 F 0.50000000000000 0.50000000000000 0.75000000000000 18.998 > 1 25 F 0.00000000000000 0.25000000000000 0.00000000000000 18.998 > 1 26 F 0.50000000000000 0.25000000000000 0.00000000000000 18.998 > 1 27 F 0.00000000000000 0.75000000000000 0.00000000000000 18.998 > 1 28 F 0.50000000000000 0.75000000000000 0.00000000000000 18.998 > 1 29 F 0.00000000000000 0.25000000000000 0.50000000000000 18.998 > 1 30 F 0.50000000000000 0.25000000000000 0.50000000000000 18.998 > 1 31 F 0.00000000000000 0.75000000000000 0.50000000000000 18.998 > 1 32 F 0.50000000000000 0.75000000000000 0.50000000000000 18.998 > 1 33 K 0.00000000000000 0.00000000000000 0.00000000000000 39.098 > 33 34 K 0.50000000000000 0.00000000000000 0.00000000000000 39.098 > 33 35 K 0.00000000000000 0.50000000000000 0.00000000000000 39.098 > 33 36 K 0.50000000000000 0.50000000000000 0.00000000000000 39.098 > 33 37 K 0.00000000000000 0.00000000000000 0.50000000000000 39.098 > 33 38 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 > 33 39 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 > 33 40 K 0.50000000000000 0.50000000000000 0.50000000000000 39.098 > 33 41 K 0.00000000000000 0.25000000000000 0.25000000000000 39.098 > 33 42 K 0.50000000000000 0.25000000000000 0.25000000000000 39.098 > 33 43 K 0.00000000000000 0.75000000000000 0.25000000000000 39.098 > 33 44 K 0.50000000000000 0.75000000000000 0.25000000000000 39.098 > 33 45 K 0.00000000000000 0.25000000000000 0.75000000000000 39.098 > 33 46 K 0.50000000000000 0.25000000000000 0.75000000000000 39.098 > 33 47 K 0.00000000000000 0.75000000000000 0.75000000000000 39.098 > 33 48 K 0.50000000000000 0.75000000000000 0.75000000000000 39.098 > 33 49 K 0.25000000000000 0.00000000000000 0.25000000000000 39.098 > 33 50 K 0.75000000000000 0.00000000000000 0.25000000000000 39.098 > 33 51 K 0.25000000000000 0.50000000000000 0.25000000000000 39.098 > 33 52 K 0.75000000000000 0.50000000000000 0.25000000000000 39.098 > 33 53 K 0.25000000000000 0.00000000000000 0.75000000000000 39.098 > 33 54 K 0.75000000000000 0.00000000000000 0.75000000000000 39.098 > 33 55 K 0.25000000000000 0.50000000000000 0.75000000000000 39.098 > 33 56 K 0.75000000000000 0.50000000000000 0.75000000000000 39.098 > 33 57 K 0.25000000000000 0.25000000000000 0.00000000000000 39.098 > 33 58 K 0.75000000000000 0.25000000000000 0.00000000000000 39.098 > 33 59 K 0.25000000000000 0.75000000000000 0.00000000000000 39.098 > 33 60 K 0.75000000000000 0.75000000000000 0.00000000000000 39.098 > 33 61 K 0.25000000000000 0.25000000000000 0.50000000000000 39.098 > 33 62 K 0.75000000000000 0.25000000000000 0.50000000000000 39.098 > 33 63 K 0.25000000000000 0.75000000000000 0.50000000000000 39.098 > 33 64 K 0.75000000000000 0.75000000000000 0.50000000000000 39.098 > 33 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 2.0006396 0.0000000 0.0000000 0.0000000 2.0006396 0.0000000 0.0000000 0.0000000 2.0006396 -------------------------- Born effective charges -------------------------- 1 F -1.1410513 0.0000000 0.0000000 0.0000000 -1.1410513 0.0000000 0.0000000 0.0000000 -1.1410513 2 K 1.1410513 0.0000000 0.0000000 0.0000000 1.1410513 0.0000000 0.0000000 0.0000000 1.1410513 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000001 (xxy) -0.00000001 (xxy) -0.00000001 (xyx) 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 ... [ 16 16 16 ] 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.23, Number of G-points: 307, Lambda: 0.18 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/145) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 145 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| -0.000 ( 0.000 0.000 0.000) 0.000 -0.000 ( 0.000 0.000 0.000) 0.000 0.000 ( 0.000 0.000 0.000) 0.000 5.617 ( 0.000 0.000 0.000) 0.000 5.617 ( 0.000 0.000 0.000) 0.000 5.617 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/145) ======================= q-point: ( 0.06 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 408 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.608 ( -17.270 17.270 17.270) 29.913 0.608 ( -17.270 17.270 17.270) 29.913 0.938 ( -26.721 26.721 26.721) 46.283 5.596 ( 1.178 -1.178 -1.178) 2.040 5.596 ( 1.178 -1.178 -1.178) 2.040 9.425 ( 2.743 -2.743 -2.743) 4.751 ======================= Grid point 2 (3/145) ======================= q-point: ( 0.12 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 417 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.202 ( -16.687 16.687 16.687) 28.903 1.202 ( -16.687 16.687 16.687) 28.903 1.862 ( -26.158 26.158 26.158) 45.308 5.539 ( 1.963 -1.963 -1.963) 3.399 5.539 ( 1.963 -1.963 -1.963) 3.399 9.286 ( 5.125 -5.125 -5.125) 8.876 ======================= Grid point 3 (4/145) ======================= q-point: ( 0.19 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 408 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.768 ( -15.659 15.659 15.659) 27.122 1.768 ( -15.659 15.659 15.659) 27.122 2.760 ( -25.168 25.168 25.168) 43.592 5.465 ( 2.230 -2.230 -2.230) 3.862 5.465 ( 2.230 -2.230 -2.230) 3.862 9.074 ( 6.973 -6.973 -6.973) 12.077 ======================= Grid point 4 (5/145) ======================= q-point: ( 0.25 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 417 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.290 ( -14.126 14.126 14.126) 24.467 2.290 ( -14.126 14.126 14.126) 24.467 3.614 ( -23.663 23.663 23.663) 40.985 5.387 ( 2.172 -2.172 -2.172) 3.762 5.387 ( 2.172 -2.172 -2.172) 3.762 8.805 ( 8.335 -8.335 -8.335) 14.436 ======================= Grid point 5 (6/145) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 408 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.748 ( -12.005 12.005 12.005) 20.793 2.748 ( -12.005 12.005 12.005) 20.793 4.405 ( -21.456 21.456 21.456) 37.162 5.313 ( 2.059 -2.059 -2.059) 3.567 5.313 ( 2.059 -2.059 -2.059) 3.567 8.496 ( 9.284 -9.284 -9.284) 16.080 ======================= Grid point 6 (7/145) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 417 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.119 ( -9.097 9.097 9.097) 15.756 3.119 ( -9.097 9.097 9.097) 15.756 5.099 ( -18.010 18.010 18.010) 31.194 5.244 ( 1.894 -1.894 -1.894) 3.280 5.244 ( 1.894 -1.894 -1.894) 3.280 8.165 ( 9.474 -9.474 -9.474) 16.410 ======================= Grid point 7 (8/145) ======================= q-point: ( 0.44 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 408 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.371 ( -5.081 5.081 5.081) 8.801 3.371 ( -5.081 5.081 5.081) 8.801 5.186 ( 1.287 -1.287 -1.287) 2.230 5.186 ( 1.287 -1.287 -1.287) 2.230 5.630 ( -11.685 11.685 11.685) 20.238 7.861 ( 7.269 -7.269 -7.269) 12.591 ======================= Grid point 8 (9/145) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 213 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.461 ( -0.000 0.000 0.000) 0.000 3.461 ( -0.000 0.000 0.000) 0.000 5.162 ( -0.000 0.000 0.000) 0.000 5.162 ( -0.000 0.000 0.000) 0.000 5.849 ( -0.000 0.000 0.000) 0.000 7.720 ( -0.000 0.000 0.000) 0.000 ======================= Grid point 18 (10/145) ======================= q-point: ( 0.06 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.530 ( 0.000 -0.000 22.396) 22.396 0.530 ( 0.000 -0.000 22.396) 22.396 1.251 ( 0.000 -0.000 52.601) 52.601 5.626 ( 0.000 -0.000 0.788) 0.788 5.626 ( 0.000 -0.000 0.788) 0.788 9.363 ( -0.000 0.000 -9.490) 9.490 ======================= Grid point 19 (11/145) ======================= q-point: ( 0.12 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.954 ( -13.029 13.029 17.712) 25.558 1.039 ( -17.768 17.768 17.722) 30.748 1.936 ( -7.638 7.638 47.625) 48.835 5.576 ( 4.131 -4.131 0.622) 5.875 5.616 ( 1.838 -1.838 1.236) 2.878 9.233 ( 0.017 -0.017 -13.388) 13.388 ======================= Grid point 20 (12/145) ======================= q-point: ( 0.19 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.481 ( -15.891 15.891 13.842) 26.394 1.664 ( -17.588 17.588 17.456) 30.387 2.690 ( -12.896 12.896 39.975) 43.939 5.475 ( 4.227 -4.227 -0.655) 6.013 5.574 ( 3.166 -3.166 1.849) 4.844 9.045 ( 1.698 -1.698 -15.478) 15.663 ======================= Grid point 21 (13/145) ======================= q-point: ( 0.25 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.001 ( -16.090 16.090 11.122) 25.327 2.247 ( -14.227 14.227 18.738) 27.494 3.455 ( -16.568 16.568 32.001) 39.661 5.373 ( 3.272 -3.272 -1.616) 4.902 5.515 ( 3.964 -3.964 2.418) 6.105 8.798 ( 3.352 -3.352 -16.572) 17.236 ======================= Grid point 22 (14/145) ======================= q-point: ( 0.31 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.481 ( -14.991 14.991 8.852) 22.974 2.754 ( -10.264 10.264 18.983) 23.897 4.195 ( -18.476 18.476 24.490) 35.812 5.288 ( 2.296 -2.296 -1.705) 3.667 5.447 ( 4.433 -4.433 2.619) 6.794 8.507 ( 4.701 -4.701 -17.099) 18.346 ======================= Grid point 23 (15/145) ======================= q-point: ( 0.38 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.898 ( -12.804 12.804 6.735) 19.319 3.158 ( -6.191 6.191 17.031) 19.150 4.867 ( -18.094 18.094 16.190) 30.280 5.232 ( 1.382 -1.382 0.439) 2.003 5.369 ( 4.685 -4.685 2.280) 7.007 8.187 ( 5.533 -5.533 -16.981) 18.697 ======================= Grid point 24 (16/145) ======================= q-point: ( 0.44 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.224 ( -9.319 9.319 4.566) 13.948 3.426 ( -1.920 1.920 12.189) 12.488 5.165 ( -0.801 0.801 -0.504) 1.240 5.284 ( 4.456 -4.456 1.533) 6.486 5.462 ( -13.270 13.270 11.757) 22.146 7.872 ( 4.777 -4.777 -15.163) 16.600 ======================= Grid point 25 (17/145) ======================= q-point: ( 0.50 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.424 ( -4.186 4.186 2.227) 6.325 3.520 ( 2.479 -2.479 4.742) 5.897 5.185 ( -0.020 0.020 3.445) 3.445 5.205 ( 3.181 -3.181 0.676) 4.548 5.781 ( -5.897 5.897 0.705) 8.369 7.663 ( -0.327 0.327 -9.134) 9.146 ======================= Grid point 26 (18/145) ======================= q-point: (-0.44 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.424 ( 6.337 -6.337 -3.087) 9.479 3.460 ( 2.277 -2.277 -0.114) 3.222 5.160 ( 0.815 -0.815 -0.152) 1.163 5.247 ( 0.068 -0.068 6.763) 6.764 5.710 ( 6.352 -6.352 -11.009) 14.209 7.712 ( -8.468 8.468 -0.916) 12.011 ======================= Grid point 27 (19/145) ======================= q-point: (-0.38 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.167 ( 9.192 -9.192 -9.375) 16.027 3.319 ( 8.524 -8.524 -2.033) 12.225 5.163 ( -1.518 1.518 -1.006) 2.371 5.250 ( 10.733 -10.733 -12.882) 19.908 5.368 ( 4.621 -4.621 4.741) 8.074 7.978 ( -12.303 12.303 2.481) 17.575 ======================= Grid point 28 (20/145) ======================= q-point: (-0.31 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.790 ( 11.200 -11.200 -14.033) 21.162 3.029 ( 13.399 -13.399 -3.342) 19.242 4.656 ( 19.646 -19.646 -18.569) 33.418 5.202 ( -3.030 3.030 -1.696) 4.608 5.433 ( 1.470 -1.470 9.897) 10.113 8.305 ( -12.765 12.765 2.679) 18.250 ======================= Grid point 29 (21/145) ======================= q-point: (-0.25 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.328 ( 12.669 -12.669 -17.375) 24.958 2.630 ( 16.972 -16.972 -4.093) 24.349 3.941 ( 23.714 -23.714 -17.013) 37.605 5.262 ( -3.948 3.948 -1.861) 5.886 5.510 ( 1.668 -1.668 9.531) 9.819 8.622 ( -12.031 12.031 1.975) 17.129 ======================= Grid point 30 (22/145) ======================= q-point: (-0.19 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.804 ( 13.886 -13.886 -19.004) 27.328 2.153 ( 19.618 -19.618 -4.199) 28.060 3.167 ( 27.225 -27.225 -13.555) 40.819 5.342 ( -4.481 4.481 -1.427) 6.495 5.576 ( 1.566 -1.566 8.265) 8.557 8.906 ( -10.747 10.747 1.092) 15.238 ======================= Grid point 31 (23/145) ======================= q-point: (-0.12 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.253 ( 15.246 -15.246 -16.507) 27.154 1.627 ( 21.580 -21.580 -3.238) 30.689 2.367 ( 30.379 -30.379 -8.213) 43.740 5.435 ( -4.449 4.449 -0.680) 6.328 5.625 ( 1.228 -1.228 5.423) 5.695 9.144 ( -8.882 8.882 0.299) 12.564 ======================= Grid point 32 (24/145) ======================= q-point: (-0.06 0.06 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 620 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.760 ( 16.550 -16.550 0.000) 23.405 1.086 ( 22.862 -22.862 0.000) 32.332 1.577 ( 32.809 -32.809 0.000) 46.398 5.526 ( -3.606 3.606 0.000) 5.100 5.636 ( 0.804 -0.804 0.000) 1.137 9.323 ( -6.311 6.311 0.000) 8.925 ======================= Grid point 35 (25/145) ======================= q-point: ( 0.12 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 349 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.035 ( 0.000 -0.000 20.886) 20.886 1.035 ( 0.000 -0.000 20.886) 20.886 2.425 ( 0.000 -0.000 47.715) 47.715 5.652 ( 0.000 -0.000 1.399) 1.399 5.652 ( 0.000 -0.000 1.399) 1.399 9.036 ( -0.000 0.000 -18.461) 18.461 ======================= Grid point 36 (26/145) ======================= q-point: ( 0.19 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.377 ( -8.628 8.628 18.072) 21.806 1.461 ( -14.658 14.658 18.109) 27.526 3.004 ( -3.814 3.814 43.114) 43.450 5.605 ( 4.957 -4.957 1.695) 7.213 5.651 ( 1.597 -1.597 1.730) 2.845 8.818 ( -1.410 1.410 -22.030) 22.120 ======================= Grid point 37 (27/145) ======================= q-point: ( 0.25 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.817 ( -12.246 12.246 14.542) 22.614 2.059 ( -18.256 18.256 16.571) 30.679 3.592 ( -6.873 6.873 36.275) 37.555 5.490 ( 6.069 -6.069 1.649) 8.740 5.622 ( 2.793 -2.793 2.145) 4.495 8.582 ( -1.201 1.201 -23.868) 23.928 ======================= Grid point 38 (28/145) ======================= q-point: ( 0.31 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.268 ( -13.230 13.230 11.508) 21.966 2.661 ( -16.096 16.096 16.960) 28.387 4.161 ( -9.981 9.981 27.373) 30.798 5.378 ( 5.045 -5.045 1.842) 7.369 5.572 ( 3.698 -3.698 2.381) 5.747 8.317 ( -0.138 0.138 -24.316) 24.317 ======================= Grid point 39 (29/145) ======================= q-point: ( 0.38 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.690 ( -12.619 12.619 8.897) 19.940 3.178 ( -10.572 10.572 17.422) 22.958 4.683 ( -12.644 12.644 15.934) 23.950 5.313 ( 3.061 -3.061 4.166) 6.008 5.504 ( 4.467 -4.467 2.222) 6.696 8.025 ( 1.016 -1.016 -23.763) 23.807 ======================= Grid point 40 (30/145) ======================= q-point: ( 0.44 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.054 ( -10.658 10.658 6.497) 16.414 3.535 ( -3.357 3.357 15.156) 15.883 5.061 ( -9.067 9.067 2.027) 12.981 5.386 ( -3.509 3.509 10.518) 11.630 5.416 ( 5.029 -5.029 1.690) 7.310 7.728 ( 1.323 -1.323 -21.776) 21.856 ======================= Grid point 41 (31/145) ======================= q-point: ( 0.50 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.327 ( -7.196 7.196 4.129) 10.982 3.673 ( 3.494 -3.494 8.750) 10.049 5.207 ( -3.149 3.149 3.513) 5.672 5.314 ( 4.969 -4.969 1.012) 7.099 5.631 ( -7.044 7.044 2.941) 10.386 7.492 ( -1.406 1.406 -16.666) 16.784 ======================= Grid point 42 (32/145) ======================= q-point: (-0.44 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.472 ( -2.054 2.054 1.798) 3.416 3.587 ( 8.412 -8.412 1.146) 11.951 5.218 ( 3.756 -3.756 0.386) 5.326 5.323 ( -1.372 1.372 6.770) 7.042 5.680 ( 1.172 -1.172 -7.313) 7.499 7.457 ( -8.555 8.555 -8.447) 14.756 ======================= Grid point 43 (33/145) ======================= q-point: (-0.38 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.328 ( 11.605 -11.605 -4.800) 17.099 3.457 ( 4.191 -4.191 -0.163) 5.929 5.156 ( 1.522 -1.522 -0.213) 2.163 5.316 ( 6.461 -6.461 -10.815) 14.158 5.531 ( 3.755 -3.755 6.414) 8.327 7.664 ( -13.917 13.917 -3.350) 19.965 ======================= Grid point 44 (34/145) ======================= q-point: (-0.31 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.950 ( 13.822 -13.822 -8.721) 21.404 3.279 ( 10.076 -10.076 -1.388) 14.318 4.872 ( 15.456 -15.456 -15.065) 26.547 5.142 ( -0.737 0.737 -0.787) 1.306 5.595 ( 1.905 -1.905 10.636) 10.972 7.980 ( -15.213 15.213 -1.817) 21.590 ======================= Grid point 45 (35/145) ======================= q-point: (-0.25 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.494 ( 15.495 -15.495 -10.382) 24.249 2.968 ( 14.710 -14.710 -1.828) 20.883 4.294 ( 20.500 -20.500 -11.481) 31.181 5.169 ( -2.467 2.467 -1.025) 3.637 5.663 ( 2.027 -2.027 8.609) 9.074 8.312 ( -14.600 14.600 -1.264) 20.686 ======================= Grid point 46 (36/145) ======================= q-point: (-0.19 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.004 ( 16.772 -16.772 -8.441) 25.176 2.564 ( 18.112 -18.112 -1.466) 25.656 3.662 ( 24.670 -24.670 -6.257) 35.446 5.232 ( -3.745 3.745 -0.702) 5.342 5.698 ( 1.855 -1.855 5.106) 5.740 8.625 ( -13.201 13.201 -0.685) 18.682 ======================= Grid point 47 (37/145) ======================= q-point: (-0.12 0.12 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 621 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.548 ( 17.271 -17.271 0.000) 24.425 2.103 ( 20.445 -20.445 0.000) 28.913 3.008 ( 28.162 -28.162 0.000) 39.827 5.325 ( -4.497 4.497 0.000) 6.360 5.690 ( 1.513 -1.513 0.000) 2.139 8.907 ( -11.302 11.302 0.000) 15.984 ======================= Grid point 52 (38/145) ======================= q-point: ( 0.19 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.496 ( 0.000 -0.000 18.523) 18.523 1.496 ( 0.000 -0.000 18.523) 18.523 3.454 ( 0.000 -0.000 40.285) 40.285 5.689 ( 0.000 -0.000 1.716) 1.716 5.689 ( 0.000 -0.000 1.716) 1.716 8.513 ( -0.000 0.000 -26.153) 26.153 ======================= Grid point 53 (39/145) ======================= q-point: ( 0.25 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.779 ( -6.333 6.333 16.183) 18.496 1.866 ( -13.071 13.071 16.432) 24.733 3.918 ( -1.867 1.867 35.031) 35.130 5.649 ( 4.792 -4.792 2.000) 7.065 5.693 ( 1.365 -1.365 1.856) 2.678 8.223 ( -2.485 2.485 -28.655) 28.870 ======================= Grid point 54 (40/145) ======================= q-point: ( 0.31 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.142 ( -9.778 9.778 13.154) 19.086 2.429 ( -18.498 18.498 15.019) 30.165 4.342 ( -2.968 2.968 27.524) 27.842 5.542 ( 6.111 -6.111 2.724) 9.062 5.670 ( 2.505 -2.505 1.981) 4.059 7.956 ( -3.415 3.415 -29.328) 29.723 ======================= Grid point 55 (41/145) ======================= q-point: ( 0.38 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.525 ( -11.069 11.069 10.297) 18.737 3.036 ( -17.723 17.723 15.087) 29.255 4.681 ( -3.817 3.817 16.357) 17.225 5.459 ( 4.377 -4.377 5.346) 8.179 5.623 ( 3.556 -3.556 1.879) 5.368 7.694 ( -2.948 2.948 -28.408) 28.712 ======================= Grid point 56 (42/145) ======================= q-point: ( 0.44 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.886 ( -10.722 10.722 7.728) 17.019 3.561 ( -11.042 11.042 15.374) 21.914 4.894 ( -5.542 5.542 2.644) 8.272 5.486 ( 0.914 -0.914 9.760) 9.845 5.548 ( 4.550 -4.550 1.499) 6.607 7.435 ( -2.253 2.253 -25.970) 26.164 ======================= Grid point 57 (43/145) ======================= q-point: ( 0.50 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.193 ( -8.870 8.870 5.346) 13.635 3.852 ( 0.451 -0.451 11.779) 11.797 5.058 ( -9.247 9.247 -2.215) 13.264 5.447 ( 5.269 -5.269 0.983) 7.516 5.604 ( -1.004 1.004 7.997) 8.122 7.215 ( -3.172 3.172 -21.300) 21.768 ======================= Grid point 58 (44/145) ======================= q-point: (-0.44 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.412 ( -5.419 5.419 3.093) 8.264 3.826 ( 8.884 -8.884 4.400) 13.312 5.253 ( -8.044 8.044 -2.074) 11.564 5.331 ( 5.259 -5.259 0.525) 7.455 5.675 ( 1.117 -1.117 3.935) 4.240 7.129 ( -7.998 7.998 -14.100) 18.076 ======================= Grid point 60 (45/145) ======================= q-point: (-0.38 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.506 ( -0.407 0.407 1.147) 1.283 3.583 ( 12.926 -12.926 -1.168) 18.318 5.224 ( 4.108 -4.108 0.169) 5.812 5.283 ( 0.879 -0.879 -10.105) 10.181 5.694 ( 2.410 -2.410 9.091) 9.709 7.267 ( -14.197 14.197 -7.779) 21.532 ======================= Grid point 61 (46/145) ======================= q-point: (-0.31 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.220 ( 15.320 -15.320 -4.011) 22.034 3.453 ( 5.454 -5.454 -0.117) 7.714 5.019 ( 9.878 -9.878 -10.614) 17.544 5.152 ( 2.009 -2.009 -0.150) 2.845 5.758 ( 2.035 -2.035 9.118) 9.561 7.569 ( -16.844 16.844 -4.141) 24.179 ======================= Grid point 62 (47/145) ======================= q-point: (-0.25 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.794 ( 16.904 -16.904 -3.917) 24.225 3.256 ( 10.906 -10.906 -0.492) 15.431 4.618 ( 15.849 -15.849 -5.888) 23.174 5.129 ( -0.275 0.275 -0.299) 0.490 5.793 ( 2.166 -2.166 5.064) 5.918 7.928 ( -16.587 16.587 -1.751) 23.523 ======================= Grid point 63 (48/145) ======================= q-point: (-0.19 0.19 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 620 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.363 ( 17.597 -17.597 0.000) 24.886 2.946 ( 15.178 -15.178 0.000) 21.466 4.156 ( 20.675 -20.675 0.000) 29.239 5.157 ( -2.256 2.256 0.000) 3.190 5.774 ( 2.035 -2.035 0.000) 2.878 8.289 ( -15.128 15.128 0.000) 21.394 ======================= Grid point 70 (49/145) ======================= q-point: ( 0.25 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 349 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.893 ( 0.000 -0.000 15.492) 15.492 1.893 ( 0.000 -0.000 15.492) 15.492 4.287 ( 0.000 -0.000 31.013) 31.013 5.730 ( 0.000 -0.000 1.714) 1.714 5.730 ( 0.000 -0.000 1.714) 1.714 7.838 ( -0.000 0.000 -31.227) 31.227 ======================= Grid point 71 (50/145) ======================= q-point: ( 0.31 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.123 ( -5.046 5.046 13.284) 15.079 2.218 ( -12.558 12.558 13.630) 22.387 4.619 ( -0.162 0.162 24.831) 24.832 5.696 ( 4.187 -4.187 2.070) 6.272 5.735 ( 1.219 -1.219 1.658) 2.392 7.511 ( -3.420 3.420 -31.613) 31.981 ======================= Grid point 72 (51/145) ======================= q-point: ( 0.38 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.421 ( -8.215 8.215 10.614) 15.736 2.750 ( -19.177 19.177 12.380) 29.813 4.847 ( 1.253 -1.253 15.155) 15.258 5.623 ( 4.379 -4.379 4.529) 7.672 5.711 ( 2.378 -2.378 1.518) 3.690 7.254 ( -4.981 4.981 -29.865) 30.684 ======================= Grid point 73 (52/145) ======================= q-point: ( 0.44 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.740 ( -9.606 9.606 8.027) 15.779 3.356 ( -19.618 19.618 12.273) 30.337 4.899 ( 2.675 -2.675 2.717) 4.658 5.630 ( 1.789 -1.789 8.784) 9.141 5.658 ( 3.554 -3.554 1.205) 5.168 7.040 ( -4.813 4.813 -26.591) 27.448 ======================= Grid point 74 (53/145) ======================= q-point: ( 0.50 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.043 ( -9.387 9.387 5.655) 14.429 3.887 ( -11.527 11.527 12.530) 20.561 4.850 ( -2.019 2.019 -5.778) 6.445 5.574 ( 4.655 -4.655 0.792) 6.631 5.712 ( 0.818 -0.818 9.153) 9.226 6.865 ( -4.762 4.762 -22.065) 23.070 ======================= Grid point 75 (54/145) ======================= q-point: ( 0.56 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.297 ( -7.596 7.596 3.495) 11.296 4.076 ( 5.692 -5.692 7.383) 10.923 4.947 ( -12.166 12.166 -6.987) 18.570 5.463 ( 5.398 -5.398 0.432) 7.647 5.788 ( 1.449 -1.449 8.398) 8.644 6.771 ( -7.118 7.118 -16.472) 19.305 ======================= Grid point 76 (55/145) ======================= q-point: (-0.38 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.467 ( -4.237 4.237 1.646) 6.214 3.887 ( 12.852 -12.852 1.184) 18.214 5.114 ( -6.855 6.855 -7.182) 12.065 5.339 ( 5.376 -5.376 0.198) 7.606 5.851 ( 1.861 -1.861 9.110) 9.483 6.842 ( -12.764 12.764 -10.288) 20.777 ======================= Grid point 77 (56/145) ======================= q-point: (-0.31 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.524 ( 0.489 -0.489 0.391) 0.794 3.553 ( 15.411 -15.411 -0.958) 21.816 5.078 ( 2.962 -2.962 -5.265) 6.728 5.226 ( 4.268 -4.268 0.044) 6.036 5.895 ( 2.178 -2.178 5.745) 6.519 7.117 ( -17.310 17.310 -4.116) 24.824 ======================= Grid point 78 (57/145) ======================= q-point: (-0.25 0.25 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 621 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.168 ( 16.744 -16.744 0.000) 23.680 3.452 ( 5.894 -5.894 0.000) 8.335 4.888 ( 10.185 -10.185 0.000) 14.403 5.150 ( 2.184 -2.184 0.000) 3.088 5.876 ( 2.308 -2.308 0.000) 3.263 7.512 ( -17.926 17.926 0.000) 25.352 ======================= Grid point 87 (58/145) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.213 ( 0.000 -0.000 11.978) 11.978 2.213 ( 0.000 -0.000 11.978) 11.978 4.884 ( 0.000 -0.000 19.840) 19.840 5.767 ( 0.000 -0.000 1.452) 1.452 5.767 ( 0.000 -0.000 1.452) 1.452 7.099 ( -0.000 0.000 -30.994) 30.994 ======================= Grid point 88 (59/145) ======================= q-point: ( 0.38 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.393 ( -4.314 4.314 9.813) 11.555 2.496 ( -12.613 12.613 10.161) 20.529 5.050 ( 2.150 -2.150 11.406) 11.805 5.749 ( 2.531 -2.531 2.647) 4.452 5.769 ( 1.179 -1.179 1.252) 2.085 6.806 ( -3.926 3.926 -27.282) 27.842 ======================= Grid point 89 (60/145) ======================= q-point: ( 0.44 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.631 ( -7.288 7.288 7.397) 12.686 2.999 ( -20.161 20.161 8.865) 29.858 5.027 ( 5.483 -5.483 0.675) 7.784 5.740 ( 2.395 -2.395 0.961) 3.521 5.757 ( 1.117 -1.117 6.472) 6.662 6.628 ( -5.118 5.118 -22.372) 23.514 ======================= Grid point 90 (61/145) ======================= q-point: ( 0.50 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.893 ( -8.728 8.728 5.075) 13.347 3.599 ( -21.834 21.834 8.400) 32.000 4.852 ( 6.259 -6.259 -5.484) 10.412 5.679 ( 3.635 -3.635 0.608) 5.176 5.829 ( 0.301 -0.301 7.835) 7.847 6.509 ( -4.833 4.833 -18.197) 19.438 ======================= Grid point 91 (62/145) ======================= q-point: ( 0.56 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.145 ( -8.603 8.603 2.979) 12.525 4.145 ( -12.953 12.953 9.737) 20.746 4.656 ( 0.550 -0.550 -9.866) 9.896 5.586 ( 4.742 -4.742 0.302) 6.714 5.914 ( 0.627 -0.627 8.493) 8.540 6.434 ( -5.622 5.622 -14.864) 16.857 ======================= Grid point 92 (63/145) ======================= q-point: (-0.38 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.352 ( -6.933 6.933 1.213) 9.879 4.195 ( 11.352 -11.352 2.773) 16.291 4.791 ( -15.001 15.001 -4.594) 21.706 5.469 ( 5.447 -5.447 0.110) 7.704 5.988 ( 1.614 -1.614 7.368) 7.713 6.459 ( -9.781 9.781 -9.402) 16.725 ======================= Grid point 93 (64/145) ======================= q-point: (-0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 620 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.487 ( -3.822 3.822 0.000) 5.405 3.898 ( 14.332 -14.332 0.000) 20.269 5.018 ( -5.947 5.947 0.000) 8.411 5.341 ( 5.403 -5.403 0.000) 7.641 5.984 ( 2.315 -2.315 0.000) 3.274 6.702 ( -15.243 15.243 0.000) 21.557 ======================= Grid point 104 (65/145) ======================= q-point: ( 0.38 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 349 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.448 ( 0.000 -0.000 8.142) 8.142 2.448 ( 0.000 -0.000 8.142) 8.142 5.176 ( 0.000 -0.000 4.193) 4.193 5.796 ( 0.000 -0.000 1.027) 1.027 5.796 ( 0.000 -0.000 1.027) 1.027 6.479 ( -0.000 0.000 -19.884) 19.884 ======================= Grid point 105 (66/145) ======================= q-point: ( 0.44 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.578 ( -3.919 3.919 6.007) 8.173 2.688 ( -12.884 12.884 6.267) 19.268 5.125 ( 3.916 -3.916 -4.232) 6.970 5.792 ( 1.209 -1.209 0.758) 1.871 5.816 ( -0.648 0.648 2.468) 2.633 6.330 ( -2.303 2.303 -12.281) 12.705 ======================= Grid point 106 (67/145) ======================= q-point: ( 0.50 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.762 ( -6.808 6.808 3.790) 10.347 3.157 ( -21.086 21.086 4.665) 30.183 4.940 ( 6.056 -6.056 -6.039) 10.480 5.757 ( 2.469 -2.469 0.448) 3.521 5.884 ( -1.906 1.906 3.731) 4.603 6.261 ( -1.928 1.928 -9.082) 9.483 ======================= Grid point 107 (68/145) ======================= q-point: ( 0.56 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.973 ( -8.322 8.322 1.734) 11.897 3.737 ( -23.962 23.962 3.188) 34.036 4.720 ( 7.778 -7.778 -4.062) 11.726 5.688 ( 3.709 -3.709 0.176) 5.248 5.990 ( -2.487 2.487 6.020) 6.972 6.197 ( -1.571 1.571 -8.982) 9.252 ======================= Grid point 108 (69/145) ======================= q-point: (-0.38 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 621 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.180 ( -8.345 8.345 0.000) 11.802 4.310 ( -22.419 22.419 0.000) 31.706 4.480 ( 10.424 -10.424 0.000) 14.742 5.590 ( 4.776 -4.776 0.000) 6.754 6.087 ( 2.039 -2.039 0.000) 2.883 6.190 ( -7.439 7.439 0.000) 10.521 ======================= Grid point 123 (70/145) ======================= q-point: ( 0.44 0.44 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.591 ( 0.000 -0.000 4.115) 4.115 2.591 ( 0.000 -0.000 4.115) 4.115 5.096 ( -0.000 0.000 -7.908) 7.908 5.814 ( 0.000 -0.000 0.525) 0.525 5.814 ( 0.000 -0.000 0.525) 0.525 6.225 ( -0.000 0.000 -3.591) 3.591 ======================= Grid point 124 (71/145) ======================= q-point: ( 0.50 0.44 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.671 ( -3.747 3.747 2.021) 5.672 2.786 ( -13.103 13.103 2.119) 18.651 4.978 ( 2.961 -2.961 -5.071) 6.576 5.804 ( 1.248 -1.248 0.252) 1.783 5.846 ( -2.228 2.228 0.346) 3.169 6.198 ( 0.391 -0.391 -1.371) 1.478 ======================= Grid point 125 (72/145) ======================= q-point: (-0.44 -0.56 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 620 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.806 ( -6.661 6.661 0.000) 9.420 3.212 ( -21.490 21.490 0.000) 30.392 4.852 ( 5.428 -5.428 0.000) 7.676 5.762 ( 2.507 -2.507 0.000) 3.546 5.924 ( -4.186 4.186 0.000) 5.920 6.167 ( 1.307 -1.307 0.000) 1.849 ======================= Grid point 143 (73/145) ======================= q-point: (-0.50 -0.50 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 205 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.639 ( 0.000 -0.000 0.000) 0.000 2.639 ( 0.000 -0.000 0.000) 0.000 4.979 ( 0.000 -0.000 0.000) 0.000 5.820 ( 0.000 -0.000 0.000) 0.000 5.820 ( 0.000 -0.000 0.000) 0.000 6.199 ( 0.000 -0.000 0.000) 0.000 ======================= Grid point 304 (74/145) ======================= q-point: ( 0.19 0.12 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.191 ( 0.000 12.288 19.279) 22.863 1.486 ( 0.000 29.829 15.545) 33.636 2.547 ( 0.000 11.307 43.963) 45.394 5.500 ( -0.000 -9.703 0.667) 9.726 5.662 ( 0.000 0.849 1.429) 1.662 9.047 ( -0.000 -1.154 -16.944) 16.984 ======================= Grid point 305 (75/145) ======================= q-point: ( 0.25 0.12 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.636 ( -11.167 14.091 16.974) 24.726 2.029 ( -4.556 27.171 15.050) 31.393 3.201 ( -5.291 13.870 37.368) 40.209 5.391 ( -0.172 -9.628 0.234) 9.632 5.661 ( 5.133 1.256 1.794) 5.580 8.815 ( -0.617 -3.286 -18.976) 19.268 ======================= Grid point 306 (76/145) ======================= q-point: ( 0.31 0.12 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.134 ( -13.280 14.938 14.065) 24.440 2.556 ( -6.304 20.941 15.847) 27.007 3.853 ( -10.092 15.732 29.430) 34.864 5.294 ( -0.607 -8.014 0.129) 8.038 5.621 ( 8.085 1.380 2.137) 8.476 8.545 ( -0.629 -5.383 -19.998) 20.719 ======================= Grid point 307 (77/145) ======================= q-point: ( 0.38 0.12 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.604 ( -11.772 15.193 11.092) 22.191 3.017 ( -5.744 13.237 16.317) 21.782 4.477 ( -14.096 16.549 20.627) 29.967 5.226 ( -0.901 -5.982 1.063) 6.142 5.557 ( 9.641 1.126 2.201) 9.953 8.242 ( -0.329 -7.088 -20.276) 21.481 ======================= Grid point 308 (78/145) ======================= q-point: ( 0.44 0.12 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.010 ( -8.758 14.343 8.039) 18.629 3.363 ( -3.465 5.356 14.685) 16.011 5.009 ( -13.490 13.529 7.489) 20.521 5.224 ( -3.320 -2.445 7.018) 8.140 5.478 ( 10.003 1.013 2.294) 10.313 7.923 ( -0.315 -7.920 -19.521) 21.069 ======================= Grid point 309 (79/145) ======================= q-point: ( 0.50 0.12 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.316 ( -4.714 11.344 4.602) 13.119 3.552 ( 0.072 -1.315 10.180) 10.265 5.167 ( -1.678 1.394 2.280) 3.156 5.340 ( 5.986 -0.795 0.931) 6.110 5.561 ( -9.540 8.220 6.896) 14.357 7.636 ( -2.233 -6.299 -16.275) 17.593 ======================= Grid point 310 (80/145) ======================= q-point: (-0.44 0.12 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1032 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.453 ( -0.042 0.000 -0.042) 0.059 3.589 ( 4.845 -0.000 4.845) 6.852 5.231 ( 2.859 -0.000 2.859) 4.043 5.260 ( 2.973 0.000 2.973) 4.205 5.730 ( -4.254 0.000 -4.254) 6.016 7.506 ( -8.812 0.000 -8.812) 12.462 ======================= Grid point 322 (81/145) ======================= q-point: ( 0.25 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.630 ( 0.000 10.864 18.074) 21.088 1.848 ( 0.000 25.777 15.325) 29.988 3.515 ( 0.000 5.563 38.393) 38.794 5.535 ( -0.000 -11.089 2.039) 11.275 5.700 ( 0.000 0.917 1.754) 1.979 8.554 ( -0.000 1.970 -25.038) 25.115 ======================= Grid point 323 (82/145) ======================= q-point: ( 0.31 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.020 ( -6.736 13.885 15.689) 22.007 2.367 ( -7.152 25.671 14.066) 30.133 4.013 ( -2.989 6.690 31.417) 32.260 5.425 ( 0.311 -11.201 2.375) 11.454 5.705 ( 4.718 1.315 1.940) 5.268 8.281 ( -1.671 0.884 -26.425) 26.492 ======================= Grid point 324 (83/145) ======================= q-point: ( 0.38 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.445 ( -7.921 15.512 12.454) 21.411 2.911 ( -10.550 20.047 14.692) 27.000 4.462 ( -6.073 7.835 21.673) 23.833 5.338 ( -0.020 -8.759 3.697) 9.507 5.671 ( 7.702 1.318 2.063) 8.082 7.997 ( -2.534 -1.014 -26.473) 26.613 ======================= Grid point 325 (84/145) ======================= q-point: ( 0.44 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.842 ( -7.079 15.490 9.275) 19.393 3.387 ( -8.762 10.799 15.387) 20.740 4.826 ( -8.960 8.743 8.228) 14.981 5.329 ( -1.373 -4.827 7.805) 9.279 5.608 ( 9.417 1.016 2.140) 9.710 7.703 ( -2.926 -2.748 -25.245) 25.562 ======================= Grid point 326 (85/145) ======================= q-point: ( 0.50 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.176 ( -5.227 13.423 6.286) 15.717 3.689 ( -2.566 0.275 12.947) 13.201 5.063 ( -8.795 7.831 -0.247) 11.779 5.411 ( -0.475 -2.494 5.225) 5.809 5.558 ( 5.506 2.535 5.317) 8.064 7.429 ( -4.033 -2.910 -21.996) 22.551 ======================= Grid point 327 (86/145) ======================= q-point: ( 0.56 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.407 ( -2.467 7.925 3.351) 8.951 3.753 ( 4.634 -5.834 6.681) 10.007 5.241 ( -6.540 5.591 1.875) 8.806 5.344 ( 6.365 -3.574 0.097) 7.301 5.647 ( -1.788 1.682 1.940) 3.129 7.258 ( -8.325 0.758 -15.465) 17.580 ======================= Grid point 328 (87/145) ======================= q-point: (-0.38 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.453 ( 1.859 -6.962 0.236) 7.210 3.654 ( 9.887 -0.742 0.715) 9.940 5.212 ( 2.625 -5.031 -3.405) 6.617 5.370 ( -0.656 2.203 0.498) 2.352 5.629 ( 3.642 -2.991 1.414) 4.920 7.309 ( -15.703 6.981 -8.172) 19.029 ======================= Grid point 329 (88/145) ======================= q-point: (-0.31 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.217 ( 8.845 -16.170 -3.581) 18.776 3.523 ( 11.210 -1.713 -2.277) 11.567 5.056 ( 1.075 -9.637 -8.762) 13.069 5.269 ( 9.728 0.376 -4.290) 10.638 5.643 ( 2.291 -1.881 10.245) 10.665 7.556 ( -20.534 9.826 -4.441) 23.193 ======================= Grid point 330 (89/145) ======================= q-point: (-0.25 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.833 ( 13.201 -17.197 -5.966) 22.486 3.277 ( 12.513 -10.955 -3.356) 16.967 4.715 ( 13.518 -13.348 -9.626) 21.297 5.187 ( 4.565 4.007 -0.774) 6.123 5.711 ( 2.043 -1.968 8.665) 9.117 7.879 ( -22.082 9.329 -2.731) 24.127 ======================= Grid point 331 (90/145) ======================= q-point: (-0.19 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.391 ( 16.039 -17.214 -5.430) 24.146 2.913 ( 12.992 -19.417 -2.555) 23.502 4.250 ( 23.703 -13.820 -4.992) 27.888 5.209 ( 1.596 6.557 -0.422) 6.762 5.744 ( 2.033 -1.953 4.940) 5.687 8.211 ( -21.953 7.471 -1.360) 23.229 ======================= Grid point 332 (91/145) ======================= q-point: (-0.12 0.19 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.959 ( 17.544 -16.504 -0.000) 24.087 2.480 ( 13.613 -25.156 0.000) 28.603 3.736 ( 31.519 -13.959 -0.000) 34.472 5.285 ( 0.598 8.759 -0.000) 8.779 5.730 ( 1.865 -1.640 -0.000) 2.483 8.528 ( -20.910 5.287 0.000) 21.568 ======================= Grid point 339 (92/145) ======================= q-point: ( 0.31 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.022 ( 0.000 10.557 15.449) 18.712 2.189 ( 0.000 22.519 13.797) 26.410 4.309 ( 0.000 2.013 29.360) 29.429 5.588 ( -0.000 -10.499 2.418) 10.774 5.741 ( 0.000 0.997 1.755) 2.019 7.903 ( -0.000 4.183 -30.292) 30.579 ======================= Grid point 340 (93/145) ======================= q-point: ( 0.38 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.356 ( -3.193 14.026 13.009) 19.395 2.679 ( -9.969 23.959 12.608) 28.851 4.629 ( -0.950 0.916 20.813) 20.854 5.498 ( 0.224 -9.893 4.021) 10.681 5.749 ( 4.431 1.398 1.761) 4.969 7.614 ( -2.646 3.757 -30.085) 30.434 ======================= Grid point 341 (94/145) ======================= q-point: ( 0.44 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.707 ( -4.160 15.457 9.970) 18.858 3.235 ( -14.411 19.769 12.978) 27.693 4.816 ( -1.584 -0.584 8.252) 8.423 5.474 ( -0.645 -5.947 7.876) 9.890 5.716 ( 7.382 1.398 1.880) 7.745 7.348 ( -4.306 2.026 -28.226) 28.625 ======================= Grid point 342 (95/145) ======================= q-point: ( 0.50 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.034 ( -4.178 14.750 7.138) 16.910 3.724 ( -10.930 9.472 13.335) 19.673 4.882 ( -4.820 2.729 -2.276) 5.989 5.524 ( 0.446 -3.994 6.673) 7.790 5.672 ( 7.112 1.455 4.282) 8.428 7.110 ( -5.491 0.744 -24.601) 25.217 ======================= Grid point 343 (96/145) ======================= q-point: ( 0.56 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.303 ( -3.134 11.873 4.614) 13.118 3.950 ( 0.531 -4.792 9.285) 10.462 5.003 ( -11.577 9.327 -5.393) 15.815 5.473 ( 6.241 -3.864 1.408) 7.474 5.707 ( 2.322 0.057 7.002) 7.377 6.944 ( -7.959 2.036 -18.908) 20.615 ======================= Grid point 344 (97/145) ======================= q-point: (-0.38 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.476 ( -0.856 5.803 2.625) 6.426 3.856 ( 9.173 -10.374 2.265) 14.032 5.144 ( -9.459 1.384 -8.291) 12.654 5.383 ( 6.366 0.159 1.680) 6.586 5.751 ( 2.490 -1.715 7.705) 8.277 6.934 ( -13.654 6.764 -12.275) 19.567 ======================= Grid point 345 (98/145) ======================= q-point: (-0.31 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.474 ( 5.434 -10.941 1.739) 12.339 3.626 ( 10.706 -2.904 -2.868) 11.457 5.060 ( -2.938 -8.325 -7.610) 11.655 5.327 ( 8.642 2.134 -1.221) 8.985 5.804 ( 1.958 -2.054 9.011) 9.447 7.129 ( -19.919 10.976 -6.785) 23.734 ======================= Grid point 346 (99/145) ======================= q-point: (-0.25 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.158 ( 14.028 -17.917 -1.337) 22.795 3.469 ( 8.181 -5.807 -1.683) 10.173 4.873 ( 4.334 -10.454 -4.720) 12.262 5.226 ( 8.888 1.893 -0.551) 9.104 5.842 ( 2.016 -2.325 5.244) 6.080 7.461 ( -22.804 11.450 -2.742) 25.664 ======================= Grid point 347 (100/145) ======================= q-point: (-0.19 0.25 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.758 ( 16.056 -18.092 0.000) 24.189 3.233 ( 10.391 -14.062 0.000) 17.485 4.600 ( 15.219 -11.907 -0.000) 19.323 5.178 ( 4.424 3.955 -0.000) 5.934 5.822 ( 2.066 -2.243 0.000) 3.050 7.839 ( -23.004 9.706 0.000) 24.968 ======================= Grid point 356 (101/145) ======================= q-point: ( 0.38 0.31 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.343 ( 0.000 10.679 12.079) 16.123 2.482 ( 0.000 20.724 11.201) 23.557 4.862 ( 0.000 -1.954 17.507) 17.616 5.650 ( -0.000 -8.319 3.142) 8.893 5.780 ( 0.000 1.078 1.490) 1.839 7.182 ( -0.000 5.653 -30.476) 30.996 ======================= Grid point 357 (102/145) ======================= q-point: ( 0.44 0.31 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.621 ( -0.631 14.193 9.664) 17.182 2.945 ( -12.517 23.323 10.061) 28.317 4.951 ( 1.291 -5.948 6.569) 8.955 5.623 ( -0.338 -5.621 6.638) 8.705 5.787 ( 4.062 1.492 1.549) 4.596 6.929 ( -3.391 5.090 -27.361) 28.036 ======================= Grid point 358 (103/145) ======================= q-point: ( 0.50 0.31 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.906 ( -1.866 15.194 6.985) 16.826 3.506 ( -17.661 20.654 10.076) 28.983 4.870 ( 1.971 -6.704 -2.589) 7.452 5.652 ( 1.401 -3.524 5.488) 6.671 5.775 ( 4.600 1.310 4.192) 6.360 6.737 ( -5.494 3.370 -23.075) 23.958 ======================= Grid point 359 (104/145) ======================= q-point: ( 0.56 0.31 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.171 ( -2.474 14.103 4.579) 15.033 4.004 ( -12.416 8.960 10.738) 18.701 4.748 ( -2.844 0.122 -8.713) 9.166 5.602 ( 6.066 -3.036 1.227) 6.893 5.822 ( 1.680 -0.284 7.626) 7.814 6.604 ( -7.134 3.030 -18.288) 19.863 ======================= Grid point 360 (105/145) ======================= q-point: (-0.38 0.31 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.388 ( -1.805 11.213 2.586) 11.647 4.119 ( 5.792 -10.705 5.198) 13.235 4.830 ( -14.960 9.667 -8.024) 19.535 5.498 ( 7.424 -2.358 0.769) 7.827 5.889 ( 1.418 -1.426 8.501) 8.735 6.567 ( -10.606 5.674 -13.382) 17.993 ======================= Grid point 362 (106/145) ======================= q-point: (-0.31 0.31 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.525 ( 0.381 6.583 1.292) 6.720 3.874 ( 11.753 -15.334 -0.127) 19.320 4.972 ( -10.092 -0.610 -4.614) 11.113 5.404 ( 8.344 0.671 0.319) 8.377 5.941 ( 1.720 -2.283 6.239) 6.862 6.703 ( -17.152 10.537 -6.505) 21.155 ======================= Grid point 363 (107/145) ======================= q-point: (-0.25 0.31 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.525 ( 13.709 -17.453 0.000) 22.194 3.556 ( 3.767 0.928 -0.000) 3.880 4.960 ( -2.898 -7.997 0.000) 8.506 5.316 ( 9.750 1.723 -0.000) 9.901 5.927 ( 2.027 -2.575 0.000) 3.277 7.038 ( -21.771 12.598 0.000) 25.153 ======================= Grid point 374 (108/145) ======================= q-point: ( 0.44 0.38 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.581 ( 0.000 10.924 8.265) 13.698 2.705 ( 0.000 20.002 7.881) 21.498 5.094 ( -0.000 -6.758 2.156) 7.094 5.741 ( 0.000 -3.594 4.444) 5.715 5.809 ( 0.000 1.145 1.056) 1.557 6.555 ( -0.000 5.190 -21.531) 22.148 ======================= Grid point 375 (109/145) ======================= q-point: ( 0.50 0.38 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.803 ( 0.960 14.327 5.950) 15.543 3.141 ( -14.442 23.554 6.553) 28.396 4.963 ( 2.621 -9.211 -4.200) 10.457 5.757 ( 2.329 -1.744 3.293) 4.394 5.836 ( 0.205 1.176 3.231) 3.445 6.407 ( -2.986 3.458 -16.302) 16.930 ======================= Grid point 376 (110/145) ======================= q-point: ( 0.56 0.38 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.030 ( -0.631 15.042 3.616) 15.483 3.695 ( -20.438 22.593 5.937) 31.039 4.751 ( 4.262 -9.381 -6.129) 11.988 5.709 ( 5.566 -1.742 0.686) 5.872 5.911 ( -1.174 0.508 6.167) 6.298 6.313 ( -4.435 2.173 -13.210) 14.102 ======================= Grid point 377 (111/145) ======================= q-point: ( 0.62 0.38 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.244 ( -1.630 13.891 1.599) 14.077 4.229 ( -16.361 12.738 8.358) 22.356 4.520 ( 1.485 -5.551 -9.319) 10.948 5.616 ( 6.854 -2.350 0.236) 7.249 6.008 ( -0.287 -0.720 8.208) 8.245 6.261 ( -6.981 3.712 -11.359) 13.840 ======================= Grid point 378 (112/145) ======================= q-point: (-0.31 0.38 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.419 ( -1.284 11.243 -0.000) 11.316 4.184 ( 10.106 -14.813 0.000) 17.932 4.718 ( -18.233 10.913 0.000) 21.250 5.507 ( 7.770 -1.674 -0.000) 7.948 6.032 ( 1.748 -2.612 0.000) 3.143 6.371 ( -12.637 8.116 0.000) 15.018 ======================= Grid point 391 (113/145) ======================= q-point: ( 0.50 0.44 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.726 ( 0.000 11.130 4.193) 11.894 2.845 ( 0.000 19.867 4.077) 20.281 5.012 ( -0.000 -6.801 -6.586) 9.467 5.819 ( 0.000 0.153 1.593) 1.601 5.828 ( 0.000 1.190 0.541) 1.307 6.237 ( -0.000 0.877 -6.201) 6.262 ======================= Grid point 392 (114/145) ======================= q-point: ( 0.56 0.44 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.896 ( 1.708 14.415 2.011) 14.654 3.245 ( -15.523 24.039 2.302) 28.708 4.844 ( 2.747 -8.344 -3.896) 9.609 5.787 ( 4.669 -0.377 0.248) 4.691 5.901 ( -4.728 2.314 1.386) 5.443 6.181 ( -0.129 -0.970 -3.903) 4.024 ======================= Grid point 393 (115/145) ======================= q-point: (-0.38 -0.56 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.072 ( -0.235 15.015 -0.000) 15.017 3.768 ( -21.917 24.100 -0.000) 32.576 4.664 ( 5.363 -10.595 0.000) 11.875 5.715 ( 5.855 -1.556 -0.000) 6.058 6.022 ( -7.021 4.138 0.000) 8.150 6.125 ( 1.135 -2.289 0.000) 2.555 ======================= Grid point 410 (116/145) ======================= q-point: (-0.44 -0.50 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 580 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.775 ( 0.000 11.209 0.000) 11.209 2.893 ( 0.000 19.890 0.000) 19.890 4.915 ( -0.000 -5.442 0.000) 5.442 5.830 ( 0.000 0.776 0.000) 0.776 5.834 ( 0.000 1.206 0.000) 1.206 6.183 ( -0.000 -1.380 0.000) 1.380 ======================= Grid point 608 (117/145) ======================= q-point: ( 0.38 0.25 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.348 ( 0.000 16.571 15.577) 22.743 2.774 ( 0.000 24.560 11.854) 27.271 4.394 ( 0.000 5.760 23.819) 24.505 5.336 ( -0.000 -9.528 3.703) 10.222 5.774 ( 0.000 1.792 1.876) 2.595 7.958 ( -0.000 -0.701 -27.429) 27.438 ======================= Grid point 609 (118/145) ======================= q-point: ( 0.44 0.25 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.736 ( -2.065 17.850 13.408) 22.420 3.196 ( -5.607 18.067 11.490) 22.133 4.687 ( -2.905 4.836 12.475) 13.691 5.309 ( -0.397 -5.718 7.217) 9.216 5.790 ( 4.844 2.076 1.871) 5.592 7.630 ( -2.304 -3.019 -27.121) 27.386 ======================= Grid point 610 (119/145) ======================= q-point: ( 0.50 0.25 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2052 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.100 ( -1.150 17.534 9.745) 20.093 3.563 ( -7.102 7.267 11.989) 15.715 4.872 ( -6.197 6.019 0.889) 8.685 5.377 ( -1.052 -2.623 9.978) 10.371 5.757 ( 8.396 1.978 1.960) 8.846 7.318 ( -4.770 -4.833 -25.064) 25.967 ======================= Grid point 611 (120/145) ======================= q-point: ( 0.56 0.25 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.379 ( 0.450 13.313 5.590) 14.446 3.779 ( -2.522 -2.236 10.120) 10.666 5.029 ( -9.922 8.072 -3.713) 13.319 5.445 ( 1.574 -3.226 5.252) 6.361 5.710 ( 7.369 1.758 3.611) 8.392 7.068 ( -7.949 -4.200 -20.293) 22.195 ======================= Grid point 612 (121/145) ======================= q-point: (-0.38 0.25 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1034 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.495 ( 2.236 0.000 2.236) 3.163 3.830 ( 4.708 -0.000 4.708) 6.658 5.145 ( -9.054 0.000 -9.054) 12.804 5.413 ( 4.745 0.000 4.745) 6.710 5.708 ( 4.173 0.000 4.173) 5.901 6.966 ( -13.380 0.000 -13.380) 18.922 ======================= Grid point 626 (122/145) ======================= q-point: ( 0.44 0.31 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.675 ( 0.000 16.998 12.438) 21.063 3.032 ( 0.000 23.443 10.212) 25.570 4.800 ( 0.000 -2.616 10.491) 10.813 5.466 ( 0.000 -6.307 7.641) 9.908 5.815 ( 0.000 1.947 1.603) 2.522 7.287 ( -0.000 2.215 -29.117) 29.201 ======================= Grid point 627 (123/145) ======================= q-point: ( 0.50 0.31 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.009 ( 1.284 18.255 9.923) 20.818 3.448 ( -9.312 17.114 10.053) 21.924 4.824 ( -0.863 -3.262 0.297) 3.387 5.522 ( -0.151 -3.203 9.717) 10.232 5.830 ( 4.593 2.185 1.610) 5.335 6.994 ( -3.408 -0.148 -26.359) 26.578 ======================= Grid point 628 (124/145) ======================= q-point: ( 0.56 0.31 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.290 ( 2.261 16.773 6.665) 18.190 3.828 ( -9.865 5.036 10.588) 15.322 4.802 ( -4.898 2.297 -6.531) 8.481 5.576 ( 2.250 -3.206 6.397) 7.501 5.811 ( 6.090 1.694 3.074) 7.029 6.764 ( -6.988 -1.150 -21.442) 22.581 ======================= Grid point 629 (125/145) ======================= q-point: ( 0.62 0.31 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.492 ( 2.665 11.319 4.206) 12.366 3.976 ( 0.336 -6.820 6.634) 9.520 4.864 ( -12.356 5.567 -9.961) 16.820 5.540 ( 6.400 -1.859 3.268) 7.423 5.818 ( 3.601 -0.112 5.805) 6.832 6.648 ( -10.905 0.850 -15.440) 18.922 ======================= Grid point 630 (126/145) ======================= q-point: (-0.31 0.31 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.569 ( 4.058 -0.980 4.173) 5.903 3.861 ( 8.307 -6.103 -1.961) 10.493 4.907 ( -10.400 -3.151 -8.990) 14.103 5.503 ( 6.837 2.563 2.309) 7.658 5.852 ( 2.186 -1.732 7.524) 8.025 6.691 ( -16.798 4.693 -9.954) 20.082 ======================= Grid point 631 (127/145) ======================= q-point: (-0.25 0.31 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.433 ( 9.494 -16.915 3.592) 19.727 3.662 ( 8.546 -2.024 -5.853) 10.554 4.858 ( -4.627 -5.676 -3.886) 8.290 5.457 ( 9.355 2.717 0.074) 9.742 5.885 ( 1.672 -2.454 5.301) 6.076 6.912 ( -23.100 6.717 -4.356) 24.448 ======================= Grid point 632 (128/145) ======================= q-point: (-0.19 0.31 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.099 ( 13.032 -18.750 0.000) 22.834 3.455 ( 8.843 -11.186 0.000) 14.259 4.800 ( 2.290 -3.761 0.000) 4.403 5.366 ( 10.481 2.742 -0.000) 10.833 5.867 ( 1.784 -2.463 0.000) 3.041 7.253 ( -26.329 5.441 0.000) 26.885 ======================= Grid point 643 (129/145) ======================= q-point: ( 0.50 0.38 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.922 ( 0.000 17.487 8.621) 19.497 3.243 ( 0.000 23.180 7.791) 24.454 4.897 ( -0.000 -8.923 -1.178) 9.000 5.664 ( 0.000 -2.928 8.208) 8.715 5.848 ( 0.000 2.081 1.142) 2.373 6.660 ( -0.000 2.938 -23.345) 23.529 ======================= Grid point 644 (130/145) ======================= q-point: ( 0.56 0.38 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.196 ( 3.730 18.468 6.137) 19.815 3.657 ( -12.575 17.902 7.690) 23.189 4.746 ( 0.709 -7.970 -6.112) 10.069 5.709 ( 2.401 -2.465 5.677) 6.639 5.872 ( 2.145 1.925 2.318) 3.699 6.461 ( -4.131 0.783 -18.388) 18.863 ======================= Grid point 645 (131/145) ======================= q-point: ( 0.62 0.38 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2052 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.411 ( 4.124 16.611 3.660) 17.502 4.059 ( -11.464 4.437 9.488) 15.528 4.585 ( -3.655 -1.771 -11.321) 12.027 5.670 ( 6.836 -1.575 2.048) 7.308 5.915 ( 0.830 0.041 5.913) 5.971 6.351 ( -7.715 1.051 -13.982) 16.003 ======================= Grid point 646 (132/145) ======================= q-point: (-0.31 0.38 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.569 ( 3.959 12.632 2.030) 13.393 4.086 ( 5.955 -14.648 2.686) 16.039 4.644 ( -16.601 5.586 -6.504) 18.684 5.593 ( 7.980 0.697 1.129) 8.089 5.973 ( 0.942 -2.001 6.390) 6.762 6.356 ( -12.383 4.129 -8.950) 15.827 ======================= Grid point 647 (133/145) ======================= q-point: (-0.25 0.38 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1098 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.651 ( 4.852 7.166 -0.000) 8.654 3.799 ( 9.655 -18.046 0.000) 20.466 4.796 ( -11.225 -2.490 0.000) 11.498 5.527 ( 8.289 2.789 -0.000) 8.746 5.968 ( 1.490 -2.830 0.000) 3.198 6.553 ( -18.721 6.969 0.000) 19.976 ======================= Grid point 660 (134/145) ======================= q-point: ( 0.56 0.44 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.074 ( 0.000 17.853 4.406) 18.389 3.388 ( -0.000 23.874 4.423) 24.281 4.803 ( -0.000 -10.552 -5.314) 11.814 5.809 ( 0.000 -1.123 3.863) 4.023 5.868 ( -0.000 2.171 0.588) 2.249 6.255 ( -0.000 0.427 -10.878) 10.887 ======================= Grid point 661 (135/145) ======================= q-point: ( 0.62 0.44 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.292 ( 4.894 18.663 2.103) 19.409 3.792 ( -14.974 20.681 3.422) 25.761 4.601 ( 2.315 -12.096 -4.505) 13.113 5.779 ( 6.786 -0.395 0.959) 6.865 5.947 ( -4.904 1.670 3.579) 6.297 6.156 ( -1.854 -1.149 -7.834) 8.132 ======================= Grid point 662 (136/145) ======================= q-point: (-0.31 0.44 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.455 ( 4.728 16.894 -0.000) 17.543 4.246 ( -20.235 15.467 0.000) 25.469 4.363 ( 5.214 -15.106 0.000) 15.981 5.692 ( 8.037 -0.275 -0.000) 8.042 6.064 ( 0.945 -2.850 0.000) 3.003 6.114 ( -9.762 3.712 0.000) 10.444 ======================= Grid point 679 (137/145) ======================= q-point: (-0.38 -0.50 0.12) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 589 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.125 ( 0.000 17.988 0.000) 17.988 3.441 ( 0.000 24.505 0.000) 24.505 4.727 ( -0.000 -10.644 0.000) 10.644 5.848 ( 0.000 0.571 0.000) 0.571 5.875 ( 0.000 2.203 0.000) 2.203 6.138 ( -0.000 -2.391 0.000) 2.391 ======================= Grid point 913 (138/145) ======================= q-point: ( 0.56 0.38 0.19) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.358 ( 0.000 19.362 9.127) 21.405 3.650 ( -0.000 9.239 8.016) 12.232 4.753 ( -0.000 -1.793 -5.656) 5.933 5.598 ( 0.000 -2.914 8.739) 9.212 5.904 ( 0.000 2.660 1.286) 2.955 6.668 ( -0.000 -2.431 -22.631) 22.761 ======================= Grid point 914 (139/145) ======================= q-point: ( 0.62 0.38 0.19) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.612 ( 7.355 14.738 6.238) 17.613 3.861 ( -8.653 0.806 8.499) 12.156 4.666 ( -3.909 0.900 -11.641) 12.313 5.643 ( 2.087 -3.004 6.755) 7.682 5.922 ( 4.217 2.272 1.606) 5.053 6.438 ( -6.092 -2.580 -16.851) 18.104 ======================= Grid point 915 (140/145) ======================= q-point: (-0.31 0.38 0.19) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1032 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.687 ( 6.153 0.000 6.153) 8.702 3.986 ( 1.964 -0.000 1.964) 2.777 4.619 ( -12.197 -0.000 -12.197) 17.249 5.648 ( 5.094 0.000 5.094) 7.205 5.916 ( 3.634 0.000 3.634) 5.140 6.363 ( -11.456 0.000 -11.456) 16.201 ======================= Grid point 931 (141/145) ======================= q-point: ( 0.62 0.44 0.19) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1096 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.520 ( 0.000 19.786 4.692) 20.335 3.828 ( -0.000 10.899 7.053) 12.982 4.565 ( -0.000 -8.247 -9.274) 12.411 5.766 ( 0.000 -2.420 5.814) 6.298 5.926 ( -0.000 2.798 0.668) 2.877 6.245 ( -0.000 -1.374 -13.184) 13.255 ======================= Grid point 932 (142/145) ======================= q-point: ( 0.69 0.44 0.19) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 2048 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.721 ( 8.520 16.783 2.718) 19.017 4.071 ( -8.673 -1.366 10.141) 13.414 4.356 ( -4.176 -4.535 -13.631) 14.960 5.768 ( 6.205 -0.426 3.367) 7.072 5.977 ( -0.942 0.753 3.848) 4.033 6.133 ( -5.980 -0.662 -9.708) 11.421 ======================= Grid point 933 (143/145) ======================= q-point: (-0.25 0.44 0.19) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 1088 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.803 ( 8.650 11.108 -0.000) 14.079 3.970 ( 4.980 -18.435 0.000) 19.096 4.458 ( -16.610 3.687 0.000) 17.014 5.711 ( 7.620 1.917 -0.000) 7.857 5.995 ( 0.786 -3.020 0.000) 3.121 6.193 ( -12.286 2.993 0.000) 12.645 ======================= Grid point 948 (144/145) ======================= q-point: (-0.31 -0.50 0.19) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 580 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.575 ( 0.000 19.938 0.000) 19.938 3.932 ( 0.000 16.010 0.000) 16.010 4.424 ( -0.000 -15.203 0.000) 15.203 5.852 ( -0.000 -0.182 0.000) 0.182 5.934 ( 0.000 2.849 0.000) 2.849 6.075 ( -0.000 -2.943 0.000) 2.943 ======================= Grid point 1216 (145/145) ======================= q-point: (-0.25 -0.50 0.25) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 1.64e-04 Number of triplets: 299 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 4.028 ( -0.000 -18.571 0.000) 18.571 4.028 ( 0.000 18.571 0.000) 18.571 4.127 ( 0.000 -0.000 -0.000) 0.000 5.848 ( 0.000 -0.000 0.000) 0.000 6.004 ( -0.000 -3.092 0.000) 3.092 6.004 ( 0.000 3.092 0.000) 3.092 =================== 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/24576 10.0 1009.929 1009.929 1009.929 -0.000 0.000 0.000 3/24576 20.0 360.309 360.309 360.309 -0.000 0.000 0.000 3/24576 30.0 166.364 166.364 166.364 -0.000 0.000 0.000 3/24576 40.0 95.673 95.673 95.673 -0.000 0.000 0.000 3/24576 50.0 65.305 65.305 65.305 -0.000 0.000 0.000 3/24576 60.0 49.361 49.361 49.361 -0.000 0.000 0.000 3/24576 70.0 39.706 39.706 39.706 -0.000 0.000 0.000 3/24576 80.0 33.272 33.272 33.272 -0.000 0.000 0.000 3/24576 90.0 28.690 28.690 28.690 -0.000 0.000 0.000 3/24576 100.0 25.260 25.260 25.260 -0.000 0.000 0.000 3/24576 110.0 22.594 22.594 22.594 -0.000 0.000 0.000 3/24576 120.0 20.460 20.460 20.460 -0.000 0.000 0.000 3/24576 130.0 18.710 18.710 18.710 -0.000 0.000 0.000 3/24576 140.0 17.248 17.248 17.248 -0.000 0.000 0.000 3/24576 150.0 16.006 16.006 16.006 -0.000 0.000 0.000 3/24576 160.0 14.937 14.937 14.937 -0.000 0.000 0.000 3/24576 170.0 14.007 14.007 14.007 -0.000 0.000 0.000 3/24576 180.0 13.189 13.189 13.189 -0.000 0.000 0.000 3/24576 190.0 12.464 12.464 12.464 -0.000 0.000 0.000 3/24576 200.0 11.817 11.817 11.817 -0.000 0.000 0.000 3/24576 210.0 11.235 11.235 11.235 -0.000 0.000 0.000 3/24576 220.0 10.710 10.710 10.710 -0.000 0.000 0.000 3/24576 230.0 10.232 10.232 10.232 -0.000 0.000 0.000 3/24576 240.0 9.796 9.796 9.796 -0.000 0.000 0.000 3/24576 250.0 9.396 9.396 9.396 -0.000 0.000 0.000 3/24576 260.0 9.028 9.028 9.028 -0.000 0.000 0.000 3/24576 270.0 8.688 8.688 8.688 -0.000 0.000 0.000 3/24576 280.0 8.373 8.373 8.373 -0.000 0.000 0.000 3/24576 290.0 8.081 8.081 8.081 -0.000 0.000 0.000 3/24576 300.0 7.808 7.808 7.808 -0.000 0.000 0.000 3/24576 310.0 7.554 7.554 7.554 -0.000 0.000 0.000 3/24576 320.0 7.316 7.316 7.316 -0.000 0.000 0.000 3/24576 330.0 7.092 7.092 7.092 -0.000 0.000 0.000 3/24576 340.0 6.882 6.882 6.882 -0.000 0.000 0.000 3/24576 350.0 6.684 6.684 6.684 -0.000 0.000 0.000 3/24576 360.0 6.498 6.498 6.498 -0.000 0.000 0.000 3/24576 370.0 6.321 6.321 6.321 -0.000 0.000 0.000 3/24576 380.0 6.154 6.154 6.154 -0.000 0.000 0.000 3/24576 390.0 5.996 5.996 5.996 -0.000 0.000 0.000 3/24576 400.0 5.846 5.846 5.846 -0.000 0.000 0.000 3/24576 410.0 5.703 5.703 5.703 -0.000 0.000 0.000 3/24576 420.0 5.566 5.566 5.566 -0.000 0.000 0.000 3/24576 430.0 5.437 5.437 5.437 -0.000 0.000 0.000 3/24576 440.0 5.313 5.313 5.313 -0.000 0.000 0.000 3/24576 450.0 5.195 5.195 5.195 -0.000 0.000 0.000 3/24576 460.0 5.082 5.082 5.082 -0.000 0.000 0.000 3/24576 470.0 4.974 4.974 4.974 -0.000 0.000 0.000 3/24576 480.0 4.870 4.870 4.870 -0.000 0.000 0.000 3/24576 490.0 4.770 4.770 4.770 -0.000 0.000 0.000 3/24576 500.0 4.675 4.675 4.675 -0.000 0.000 0.000 3/24576 510.0 4.583 4.583 4.583 -0.000 0.000 0.000 3/24576 520.0 4.495 4.495 4.495 -0.000 0.000 0.000 3/24576 530.0 4.411 4.411 4.411 -0.000 0.000 0.000 3/24576 540.0 4.329 4.329 4.329 -0.000 0.000 0.000 3/24576 550.0 4.250 4.250 4.250 -0.000 0.000 0.000 3/24576 560.0 4.174 4.174 4.174 -0.000 0.000 0.000 3/24576 570.0 4.101 4.101 4.101 -0.000 0.000 0.000 3/24576 580.0 4.031 4.031 4.031 -0.000 0.000 0.000 3/24576 590.0 3.962 3.962 3.962 -0.000 0.000 0.000 3/24576 600.0 3.896 3.896 3.896 -0.000 0.000 0.000 3/24576 610.0 3.833 3.833 3.833 -0.000 0.000 0.000 3/24576 620.0 3.771 3.771 3.771 -0.000 0.000 0.000 3/24576 630.0 3.711 3.711 3.711 -0.000 0.000 0.000 3/24576 640.0 3.653 3.653 3.653 -0.000 0.000 0.000 3/24576 650.0 3.597 3.597 3.597 -0.000 0.000 0.000 3/24576 660.0 3.543 3.543 3.543 -0.000 0.000 0.000 3/24576 670.0 3.490 3.490 3.490 -0.000 0.000 0.000 3/24576 680.0 3.439 3.439 3.439 -0.000 0.000 0.000 3/24576 690.0 3.389 3.389 3.389 -0.000 0.000 0.000 3/24576 700.0 3.341 3.341 3.341 -0.000 0.000 0.000 3/24576 710.0 3.294 3.294 3.294 -0.000 0.000 0.000 3/24576 720.0 3.248 3.248 3.248 -0.000 0.000 0.000 3/24576 730.0 3.204 3.204 3.204 -0.000 0.000 0.000 3/24576 740.0 3.161 3.161 3.161 -0.000 0.000 0.000 3/24576 750.0 3.118 3.118 3.118 -0.000 0.000 0.000 3/24576 760.0 3.078 3.078 3.078 -0.000 0.000 0.000 3/24576 770.0 3.038 3.038 3.038 -0.000 0.000 0.000 3/24576 780.0 2.999 2.999 2.999 -0.000 0.000 0.000 3/24576 790.0 2.961 2.961 2.961 -0.000 0.000 0.000 3/24576 800.0 2.924 2.924 2.924 -0.000 0.000 0.000 3/24576 810.0 2.888 2.888 2.888 -0.000 0.000 0.000 3/24576 820.0 2.853 2.853 2.853 -0.000 0.000 0.000 3/24576 830.0 2.819 2.819 2.819 -0.000 0.000 0.000 3/24576 840.0 2.785 2.785 2.785 -0.000 0.000 0.000 3/24576 850.0 2.752 2.752 2.752 -0.000 0.000 0.000 3/24576 860.0 2.720 2.720 2.720 -0.000 0.000 0.000 3/24576 870.0 2.689 2.689 2.689 -0.000 0.000 0.000 3/24576 880.0 2.659 2.659 2.659 -0.000 0.000 0.000 3/24576 890.0 2.629 2.629 2.629 -0.000 0.000 0.000 3/24576 900.0 2.600 2.600 2.600 -0.000 0.000 0.000 3/24576 910.0 2.571 2.571 2.571 -0.000 0.000 0.000 3/24576 920.0 2.543 2.543 2.543 -0.000 0.000 0.000 3/24576 930.0 2.516 2.516 2.516 -0.000 0.000 0.000 3/24576 940.0 2.490 2.490 2.490 -0.000 0.000 0.000 3/24576 950.0 2.463 2.463 2.463 -0.000 0.000 0.000 3/24576 960.0 2.438 2.438 2.438 -0.000 0.000 0.000 3/24576 970.0 2.413 2.413 2.413 -0.000 0.000 0.000 3/24576 980.0 2.388 2.388 2.388 -0.000 0.000 0.000 3/24576 990.0 2.364 2.364 2.364 -0.000 0.000 0.000 3/24576 1000.0 2.341 2.341 2.341 -0.000 0.000 0.000 3/24576 Thermal conductivity related properties were written into "kappa-m161616.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 18:31:58]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|