----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 18:32:49]------------------------- 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: [4 4 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) Number of symmetry operations in supercell: 384 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.175399543208685 0.000000000000000 0.000000000000000 b -2.087699771604342 3.616002075368662 0.000000000000000 c 0.000000000000000 0.000000000000000 6.885862750000000 Atomic positions (fractional): *1 Cu 0.33333333333333 0.66666666666667 0.62650592587408 63.546 2 Cu 0.66666666666667 0.33333333333333 0.12650592587408 63.546 *3 I 0.33333333333333 0.66666666666667 0.99849407412592 126.904 4 I 0.66666666666667 0.33333333333333 0.49849407412592 126.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.175399543208685 0.000000000000000 0.000000000000000 b -2.087699771604342 3.616002075368662 0.000000000000000 c 0.000000000000000 0.000000000000000 6.885862750000000 Atomic positions (fractional): *1 Cu 0.33333333333333 0.66666666666667 0.62650592587408 63.546 > 1 2 Cu 0.66666666666667 0.33333333333333 0.12650592587408 63.546 > 2 *3 I 0.33333333333333 0.66666666666667 0.99849407412592 126.904 > 3 4 I 0.66666666666667 0.33333333333333 0.49849407412592 126.904 > 4 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.701598172834739 0.000000000000000 0.000000000000000 b -8.350799086417368 14.464008301474648 0.000000000000000 c 0.000000000000000 0.000000000000000 13.771725500000001 Atomic positions (fractional): *1 Cu 0.08333333333333 0.16666666666667 0.31325296293704 63.546 > 1 2 Cu 0.33333333333333 0.16666666666667 0.31325296293704 63.546 > 1 3 Cu 0.58333333333333 0.16666666666667 0.31325296293704 63.546 > 1 4 Cu 0.83333333333333 0.16666666666667 0.31325296293704 63.546 > 1 5 Cu 0.08333333333333 0.41666666666667 0.31325296293704 63.546 > 1 6 Cu 0.33333333333333 0.41666666666667 0.31325296293704 63.546 > 1 7 Cu 0.58333333333333 0.41666666666667 0.31325296293704 63.546 > 1 8 Cu 0.83333333333333 0.41666666666667 0.31325296293704 63.546 > 1 9 Cu 0.08333333333333 0.66666666666667 0.31325296293704 63.546 > 1 10 Cu 0.33333333333333 0.66666666666667 0.31325296293704 63.546 > 1 11 Cu 0.58333333333333 0.66666666666667 0.31325296293704 63.546 > 1 12 Cu 0.83333333333333 0.66666666666667 0.31325296293704 63.546 > 1 13 Cu 0.08333333333333 0.91666666666667 0.31325296293704 63.546 > 1 14 Cu 0.33333333333333 0.91666666666667 0.31325296293704 63.546 > 1 15 Cu 0.58333333333333 0.91666666666667 0.31325296293704 63.546 > 1 16 Cu 0.83333333333333 0.91666666666667 0.31325296293704 63.546 > 1 17 Cu 0.08333333333333 0.16666666666667 0.81325296293704 63.546 > 1 18 Cu 0.33333333333333 0.16666666666667 0.81325296293704 63.546 > 1 19 Cu 0.58333333333333 0.16666666666667 0.81325296293704 63.546 > 1 20 Cu 0.83333333333333 0.16666666666667 0.81325296293704 63.546 > 1 21 Cu 0.08333333333333 0.41666666666667 0.81325296293704 63.546 > 1 22 Cu 0.33333333333333 0.41666666666667 0.81325296293704 63.546 > 1 23 Cu 0.58333333333333 0.41666666666667 0.81325296293704 63.546 > 1 24 Cu 0.83333333333333 0.41666666666667 0.81325296293704 63.546 > 1 25 Cu 0.08333333333333 0.66666666666667 0.81325296293704 63.546 > 1 26 Cu 0.33333333333333 0.66666666666667 0.81325296293704 63.546 > 1 27 Cu 0.58333333333333 0.66666666666667 0.81325296293704 63.546 > 1 28 Cu 0.83333333333333 0.66666666666667 0.81325296293704 63.546 > 1 29 Cu 0.08333333333333 0.91666666666667 0.81325296293704 63.546 > 1 30 Cu 0.33333333333333 0.91666666666667 0.81325296293704 63.546 > 1 31 Cu 0.58333333333333 0.91666666666667 0.81325296293704 63.546 > 1 32 Cu 0.83333333333333 0.91666666666667 0.81325296293704 63.546 > 1 33 Cu 0.16666666666667 0.08333333333333 0.06325296293704 63.546 > 2 34 Cu 0.41666666666667 0.08333333333333 0.06325296293704 63.546 > 2 35 Cu 0.66666666666667 0.08333333333333 0.06325296293704 63.546 > 2 36 Cu 0.91666666666667 0.08333333333333 0.06325296293704 63.546 > 2 37 Cu 0.16666666666667 0.33333333333333 0.06325296293704 63.546 > 2 38 Cu 0.41666666666667 0.33333333333333 0.06325296293704 63.546 > 2 39 Cu 0.66666666666667 0.33333333333333 0.06325296293704 63.546 > 2 40 Cu 0.91666666666667 0.33333333333333 0.06325296293704 63.546 > 2 41 Cu 0.16666666666667 0.58333333333333 0.06325296293704 63.546 > 2 42 Cu 0.41666666666667 0.58333333333333 0.06325296293704 63.546 > 2 43 Cu 0.66666666666667 0.58333333333333 0.06325296293704 63.546 > 2 44 Cu 0.91666666666667 0.58333333333333 0.06325296293704 63.546 > 2 45 Cu 0.16666666666667 0.83333333333333 0.06325296293704 63.546 > 2 46 Cu 0.41666666666667 0.83333333333333 0.06325296293704 63.546 > 2 47 Cu 0.66666666666667 0.83333333333333 0.06325296293704 63.546 > 2 48 Cu 0.91666666666667 0.83333333333333 0.06325296293704 63.546 > 2 49 Cu 0.16666666666667 0.08333333333333 0.56325296293704 63.546 > 2 50 Cu 0.41666666666667 0.08333333333333 0.56325296293704 63.546 > 2 51 Cu 0.66666666666667 0.08333333333333 0.56325296293704 63.546 > 2 52 Cu 0.91666666666667 0.08333333333333 0.56325296293704 63.546 > 2 53 Cu 0.16666666666667 0.33333333333333 0.56325296293704 63.546 > 2 54 Cu 0.41666666666667 0.33333333333333 0.56325296293704 63.546 > 2 55 Cu 0.66666666666667 0.33333333333333 0.56325296293704 63.546 > 2 56 Cu 0.91666666666667 0.33333333333333 0.56325296293704 63.546 > 2 57 Cu 0.16666666666667 0.58333333333333 0.56325296293704 63.546 > 2 58 Cu 0.41666666666667 0.58333333333333 0.56325296293704 63.546 > 2 59 Cu 0.66666666666667 0.58333333333333 0.56325296293704 63.546 > 2 60 Cu 0.91666666666667 0.58333333333333 0.56325296293704 63.546 > 2 61 Cu 0.16666666666667 0.83333333333333 0.56325296293704 63.546 > 2 62 Cu 0.41666666666667 0.83333333333333 0.56325296293704 63.546 > 2 63 Cu 0.66666666666667 0.83333333333333 0.56325296293704 63.546 > 2 64 Cu 0.91666666666667 0.83333333333333 0.56325296293704 63.546 > 2 *65 I 0.08333333333333 0.16666666666667 0.49924703706296 126.904 > 3 66 I 0.33333333333333 0.16666666666667 0.49924703706296 126.904 > 3 67 I 0.58333333333333 0.16666666666667 0.49924703706296 126.904 > 3 68 I 0.83333333333333 0.16666666666667 0.49924703706296 126.904 > 3 69 I 0.08333333333333 0.41666666666667 0.49924703706296 126.904 > 3 70 I 0.33333333333333 0.41666666666667 0.49924703706296 126.904 > 3 71 I 0.58333333333333 0.41666666666667 0.49924703706296 126.904 > 3 72 I 0.83333333333333 0.41666666666667 0.49924703706296 126.904 > 3 73 I 0.08333333333333 0.66666666666667 0.49924703706296 126.904 > 3 74 I 0.33333333333333 0.66666666666667 0.49924703706296 126.904 > 3 75 I 0.58333333333333 0.66666666666667 0.49924703706296 126.904 > 3 76 I 0.83333333333333 0.66666666666667 0.49924703706296 126.904 > 3 77 I 0.08333333333333 0.91666666666667 0.49924703706296 126.904 > 3 78 I 0.33333333333333 0.91666666666667 0.49924703706296 126.904 > 3 79 I 0.58333333333333 0.91666666666667 0.49924703706296 126.904 > 3 80 I 0.83333333333333 0.91666666666667 0.49924703706296 126.904 > 3 81 I 0.08333333333333 0.16666666666667 0.99924703706296 126.904 > 3 82 I 0.33333333333333 0.16666666666667 0.99924703706296 126.904 > 3 83 I 0.58333333333333 0.16666666666667 0.99924703706296 126.904 > 3 84 I 0.83333333333333 0.16666666666667 0.99924703706296 126.904 > 3 85 I 0.08333333333333 0.41666666666667 0.99924703706296 126.904 > 3 86 I 0.33333333333333 0.41666666666667 0.99924703706296 126.904 > 3 87 I 0.58333333333333 0.41666666666667 0.99924703706296 126.904 > 3 88 I 0.83333333333333 0.41666666666667 0.99924703706296 126.904 > 3 89 I 0.08333333333333 0.66666666666667 0.99924703706296 126.904 > 3 90 I 0.33333333333333 0.66666666666667 0.99924703706296 126.904 > 3 91 I 0.58333333333333 0.66666666666667 0.99924703706296 126.904 > 3 92 I 0.83333333333333 0.66666666666667 0.99924703706296 126.904 > 3 93 I 0.08333333333333 0.91666666666667 0.99924703706296 126.904 > 3 94 I 0.33333333333333 0.91666666666667 0.99924703706296 126.904 > 3 95 I 0.58333333333333 0.91666666666667 0.99924703706296 126.904 > 3 96 I 0.83333333333333 0.91666666666667 0.99924703706296 126.904 > 3 97 I 0.16666666666667 0.08333333333333 0.24924703706296 126.904 > 4 98 I 0.41666666666667 0.08333333333333 0.24924703706296 126.904 > 4 99 I 0.66666666666667 0.08333333333333 0.24924703706296 126.904 > 4 100 I 0.91666666666667 0.08333333333333 0.24924703706296 126.904 > 4 101 I 0.16666666666667 0.33333333333333 0.24924703706296 126.904 > 4 102 I 0.41666666666667 0.33333333333333 0.24924703706296 126.904 > 4 103 I 0.66666666666667 0.33333333333333 0.24924703706296 126.904 > 4 104 I 0.91666666666667 0.33333333333333 0.24924703706296 126.904 > 4 105 I 0.16666666666667 0.58333333333333 0.24924703706296 126.904 > 4 106 I 0.41666666666667 0.58333333333333 0.24924703706296 126.904 > 4 107 I 0.66666666666667 0.58333333333333 0.24924703706296 126.904 > 4 108 I 0.91666666666667 0.58333333333333 0.24924703706296 126.904 > 4 109 I 0.16666666666667 0.83333333333333 0.24924703706296 126.904 > 4 110 I 0.41666666666667 0.83333333333333 0.24924703706296 126.904 > 4 111 I 0.66666666666667 0.83333333333333 0.24924703706296 126.904 > 4 112 I 0.91666666666667 0.83333333333333 0.24924703706296 126.904 > 4 113 I 0.16666666666667 0.08333333333333 0.74924703706296 126.904 > 4 114 I 0.41666666666667 0.08333333333333 0.74924703706296 126.904 > 4 115 I 0.66666666666667 0.08333333333333 0.74924703706296 126.904 > 4 116 I 0.91666666666667 0.08333333333333 0.74924703706296 126.904 > 4 117 I 0.16666666666667 0.33333333333333 0.74924703706296 126.904 > 4 118 I 0.41666666666667 0.33333333333333 0.74924703706296 126.904 > 4 119 I 0.66666666666667 0.33333333333333 0.74924703706296 126.904 > 4 120 I 0.91666666666667 0.33333333333333 0.74924703706296 126.904 > 4 121 I 0.16666666666667 0.58333333333333 0.74924703706296 126.904 > 4 122 I 0.41666666666667 0.58333333333333 0.74924703706296 126.904 > 4 123 I 0.66666666666667 0.58333333333333 0.74924703706296 126.904 > 4 124 I 0.91666666666667 0.58333333333333 0.74924703706296 126.904 > 4 125 I 0.16666666666667 0.83333333333333 0.74924703706296 126.904 > 4 126 I 0.41666666666667 0.83333333333333 0.74924703706296 126.904 > 4 127 I 0.66666666666667 0.83333333333333 0.74924703706296 126.904 > 4 128 I 0.91666666666667 0.83333333333333 0.74924703706296 126.904 > 4 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.1532900 -0.0000000 0.0000000 -0.0000000 7.1532900 0.0000000 0.0000000 0.0000000 7.1937858 -------------------------- Born effective charges -------------------------- 1 Cu 1.0026013 -0.0000000 0.0000000 -0.0000000 1.0026013 0.0000000 0.0000000 0.0000000 1.0429109 2 Cu 1.0026013 -0.0000000 0.0000000 -0.0000000 1.0026013 0.0000000 0.0000000 0.0000000 1.0429109 3 I -1.0026013 0.0000000 0.0000000 0.0000000 -1.0026013 0.0000000 0.0000000 0.0000000 -1.0429109 4 I -1.0026013 0.0000000 0.0000000 0.0000000 -1.0026013 0.0000000 0.0000000 0.0000000 -1.0429109 ---------------------------------------------------------------------------- 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: 384/384 Permutation basis: 2856/2856 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 224 Number of blocks in projector: 188 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: 98 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 46 Use standard eigh solver. Tree of FC basis block matrices: - (224, 220), data: False |-- (46, 44), data: True |-- (80, 80), data: True |-- (98, 96), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.010 Solver_block: 80 / 80 - Time: 0.158 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.169 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (yy) 0.00000000 (yy) Permutation basis: 384/384 Permutation basis: 2856/2856 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 224 Number of blocks in projector: 188 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: 98 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 46 Use standard eigh solver. Tree of FC basis block matrices: - (224, 220), data: False |-- (46, 44), data: True |-- (80, 80), data: True |-- (98, 96), data: True Max drift after symmetrization by symfc projector: 0.00000000 (yy) 0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-08 18:32:53]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 18:32: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: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.175399543208685 0.000000000000000 0.000000000000000 b -2.087699771604342 3.616002075368662 0.000000000000000 c 0.000000000000000 0.000000000000000 6.885862750000000 Atomic positions (fractional): 1 Cu 0.33333333333333 0.66666666666667 0.62650592587408 63.546 2 Cu 0.66666666666667 0.33333333333333 0.12650592587408 63.546 3 I 0.33333333333333 0.66666666666667 0.99849407412592 126.904 4 I 0.66666666666667 0.33333333333333 0.49849407412592 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.701598172834739 0.000000000000000 0.000000000000000 b -8.350799086417368 14.464008301474648 0.000000000000000 c 0.000000000000000 0.000000000000000 13.771725500000001 Atomic positions (fractional): 1 Cu 0.08333333333333 0.16666666666667 0.31325296293704 63.546 > 1 2 Cu 0.33333333333333 0.16666666666667 0.31325296293704 63.546 > 1 3 Cu 0.58333333333333 0.16666666666667 0.31325296293704 63.546 > 1 4 Cu 0.83333333333333 0.16666666666667 0.31325296293704 63.546 > 1 5 Cu 0.08333333333333 0.41666666666667 0.31325296293704 63.546 > 1 6 Cu 0.33333333333333 0.41666666666667 0.31325296293704 63.546 > 1 7 Cu 0.58333333333333 0.41666666666667 0.31325296293704 63.546 > 1 8 Cu 0.83333333333333 0.41666666666667 0.31325296293704 63.546 > 1 9 Cu 0.08333333333333 0.66666666666667 0.31325296293704 63.546 > 1 10 Cu 0.33333333333333 0.66666666666667 0.31325296293704 63.546 > 1 11 Cu 0.58333333333333 0.66666666666667 0.31325296293704 63.546 > 1 12 Cu 0.83333333333333 0.66666666666667 0.31325296293704 63.546 > 1 13 Cu 0.08333333333333 0.91666666666667 0.31325296293704 63.546 > 1 14 Cu 0.33333333333333 0.91666666666667 0.31325296293704 63.546 > 1 15 Cu 0.58333333333333 0.91666666666667 0.31325296293704 63.546 > 1 16 Cu 0.83333333333333 0.91666666666667 0.31325296293704 63.546 > 1 17 Cu 0.08333333333333 0.16666666666667 0.81325296293704 63.546 > 1 18 Cu 0.33333333333333 0.16666666666667 0.81325296293704 63.546 > 1 19 Cu 0.58333333333333 0.16666666666667 0.81325296293704 63.546 > 1 20 Cu 0.83333333333333 0.16666666666667 0.81325296293704 63.546 > 1 21 Cu 0.08333333333333 0.41666666666667 0.81325296293704 63.546 > 1 22 Cu 0.33333333333333 0.41666666666667 0.81325296293704 63.546 > 1 23 Cu 0.58333333333333 0.41666666666667 0.81325296293704 63.546 > 1 24 Cu 0.83333333333333 0.41666666666667 0.81325296293704 63.546 > 1 25 Cu 0.08333333333333 0.66666666666667 0.81325296293704 63.546 > 1 26 Cu 0.33333333333333 0.66666666666667 0.81325296293704 63.546 > 1 27 Cu 0.58333333333333 0.66666666666667 0.81325296293704 63.546 > 1 28 Cu 0.83333333333333 0.66666666666667 0.81325296293704 63.546 > 1 29 Cu 0.08333333333333 0.91666666666667 0.81325296293704 63.546 > 1 30 Cu 0.33333333333333 0.91666666666667 0.81325296293704 63.546 > 1 31 Cu 0.58333333333333 0.91666666666667 0.81325296293704 63.546 > 1 32 Cu 0.83333333333333 0.91666666666667 0.81325296293704 63.546 > 1 33 Cu 0.16666666666667 0.08333333333333 0.06325296293704 63.546 > 33 34 Cu 0.41666666666667 0.08333333333333 0.06325296293704 63.546 > 33 35 Cu 0.66666666666667 0.08333333333333 0.06325296293704 63.546 > 33 36 Cu 0.91666666666667 0.08333333333333 0.06325296293704 63.546 > 33 37 Cu 0.16666666666667 0.33333333333333 0.06325296293704 63.546 > 33 38 Cu 0.41666666666667 0.33333333333333 0.06325296293704 63.546 > 33 39 Cu 0.66666666666667 0.33333333333333 0.06325296293704 63.546 > 33 40 Cu 0.91666666666667 0.33333333333333 0.06325296293704 63.546 > 33 41 Cu 0.16666666666667 0.58333333333333 0.06325296293704 63.546 > 33 42 Cu 0.41666666666667 0.58333333333333 0.06325296293704 63.546 > 33 43 Cu 0.66666666666667 0.58333333333333 0.06325296293704 63.546 > 33 44 Cu 0.91666666666667 0.58333333333333 0.06325296293704 63.546 > 33 45 Cu 0.16666666666667 0.83333333333333 0.06325296293704 63.546 > 33 46 Cu 0.41666666666667 0.83333333333333 0.06325296293704 63.546 > 33 47 Cu 0.66666666666667 0.83333333333333 0.06325296293704 63.546 > 33 48 Cu 0.91666666666667 0.83333333333333 0.06325296293704 63.546 > 33 49 Cu 0.16666666666667 0.08333333333333 0.56325296293704 63.546 > 33 50 Cu 0.41666666666667 0.08333333333333 0.56325296293704 63.546 > 33 51 Cu 0.66666666666667 0.08333333333333 0.56325296293704 63.546 > 33 52 Cu 0.91666666666667 0.08333333333333 0.56325296293704 63.546 > 33 53 Cu 0.16666666666667 0.33333333333333 0.56325296293704 63.546 > 33 54 Cu 0.41666666666667 0.33333333333333 0.56325296293704 63.546 > 33 55 Cu 0.66666666666667 0.33333333333333 0.56325296293704 63.546 > 33 56 Cu 0.91666666666667 0.33333333333333 0.56325296293704 63.546 > 33 57 Cu 0.16666666666667 0.58333333333333 0.56325296293704 63.546 > 33 58 Cu 0.41666666666667 0.58333333333333 0.56325296293704 63.546 > 33 59 Cu 0.66666666666667 0.58333333333333 0.56325296293704 63.546 > 33 60 Cu 0.91666666666667 0.58333333333333 0.56325296293704 63.546 > 33 61 Cu 0.16666666666667 0.83333333333333 0.56325296293704 63.546 > 33 62 Cu 0.41666666666667 0.83333333333333 0.56325296293704 63.546 > 33 63 Cu 0.66666666666667 0.83333333333333 0.56325296293704 63.546 > 33 64 Cu 0.91666666666667 0.83333333333333 0.56325296293704 63.546 > 33 65 I 0.08333333333333 0.16666666666667 0.49924703706296 126.904 > 65 66 I 0.33333333333333 0.16666666666667 0.49924703706296 126.904 > 65 67 I 0.58333333333333 0.16666666666667 0.49924703706296 126.904 > 65 68 I 0.83333333333333 0.16666666666667 0.49924703706296 126.904 > 65 69 I 0.08333333333333 0.41666666666667 0.49924703706296 126.904 > 65 70 I 0.33333333333333 0.41666666666667 0.49924703706296 126.904 > 65 71 I 0.58333333333333 0.41666666666667 0.49924703706296 126.904 > 65 72 I 0.83333333333333 0.41666666666667 0.49924703706296 126.904 > 65 73 I 0.08333333333333 0.66666666666667 0.49924703706296 126.904 > 65 74 I 0.33333333333333 0.66666666666667 0.49924703706296 126.904 > 65 75 I 0.58333333333333 0.66666666666667 0.49924703706296 126.904 > 65 76 I 0.83333333333333 0.66666666666667 0.49924703706296 126.904 > 65 77 I 0.08333333333333 0.91666666666667 0.49924703706296 126.904 > 65 78 I 0.33333333333333 0.91666666666667 0.49924703706296 126.904 > 65 79 I 0.58333333333333 0.91666666666667 0.49924703706296 126.904 > 65 80 I 0.83333333333333 0.91666666666667 0.49924703706296 126.904 > 65 81 I 0.08333333333333 0.16666666666667 0.99924703706296 126.904 > 65 82 I 0.33333333333333 0.16666666666667 0.99924703706296 126.904 > 65 83 I 0.58333333333333 0.16666666666667 0.99924703706296 126.904 > 65 84 I 0.83333333333333 0.16666666666667 0.99924703706296 126.904 > 65 85 I 0.08333333333333 0.41666666666667 0.99924703706296 126.904 > 65 86 I 0.33333333333333 0.41666666666667 0.99924703706296 126.904 > 65 87 I 0.58333333333333 0.41666666666667 0.99924703706296 126.904 > 65 88 I 0.83333333333333 0.41666666666667 0.99924703706296 126.904 > 65 89 I 0.08333333333333 0.66666666666667 0.99924703706296 126.904 > 65 90 I 0.33333333333333 0.66666666666667 0.99924703706296 126.904 > 65 91 I 0.58333333333333 0.66666666666667 0.99924703706296 126.904 > 65 92 I 0.83333333333333 0.66666666666667 0.99924703706296 126.904 > 65 93 I 0.08333333333333 0.91666666666667 0.99924703706296 126.904 > 65 94 I 0.33333333333333 0.91666666666667 0.99924703706296 126.904 > 65 95 I 0.58333333333333 0.91666666666667 0.99924703706296 126.904 > 65 96 I 0.83333333333333 0.91666666666667 0.99924703706296 126.904 > 65 97 I 0.16666666666667 0.08333333333333 0.24924703706296 126.904 > 97 98 I 0.41666666666667 0.08333333333333 0.24924703706296 126.904 > 97 99 I 0.66666666666667 0.08333333333333 0.24924703706296 126.904 > 97 100 I 0.91666666666667 0.08333333333333 0.24924703706296 126.904 > 97 101 I 0.16666666666667 0.33333333333333 0.24924703706296 126.904 > 97 102 I 0.41666666666667 0.33333333333333 0.24924703706296 126.904 > 97 103 I 0.66666666666667 0.33333333333333 0.24924703706296 126.904 > 97 104 I 0.91666666666667 0.33333333333333 0.24924703706296 126.904 > 97 105 I 0.16666666666667 0.58333333333333 0.24924703706296 126.904 > 97 106 I 0.41666666666667 0.58333333333333 0.24924703706296 126.904 > 97 107 I 0.66666666666667 0.58333333333333 0.24924703706296 126.904 > 97 108 I 0.91666666666667 0.58333333333333 0.24924703706296 126.904 > 97 109 I 0.16666666666667 0.83333333333333 0.24924703706296 126.904 > 97 110 I 0.41666666666667 0.83333333333333 0.24924703706296 126.904 > 97 111 I 0.66666666666667 0.83333333333333 0.24924703706296 126.904 > 97 112 I 0.91666666666667 0.83333333333333 0.24924703706296 126.904 > 97 113 I 0.16666666666667 0.08333333333333 0.74924703706296 126.904 > 97 114 I 0.41666666666667 0.08333333333333 0.74924703706296 126.904 > 97 115 I 0.66666666666667 0.08333333333333 0.74924703706296 126.904 > 97 116 I 0.91666666666667 0.08333333333333 0.74924703706296 126.904 > 97 117 I 0.16666666666667 0.33333333333333 0.74924703706296 126.904 > 97 118 I 0.41666666666667 0.33333333333333 0.74924703706296 126.904 > 97 119 I 0.66666666666667 0.33333333333333 0.74924703706296 126.904 > 97 120 I 0.91666666666667 0.33333333333333 0.74924703706296 126.904 > 97 121 I 0.16666666666667 0.58333333333333 0.74924703706296 126.904 > 97 122 I 0.41666666666667 0.58333333333333 0.74924703706296 126.904 > 97 123 I 0.66666666666667 0.58333333333333 0.74924703706296 126.904 > 97 124 I 0.91666666666667 0.58333333333333 0.74924703706296 126.904 > 97 125 I 0.16666666666667 0.83333333333333 0.74924703706296 126.904 > 97 126 I 0.41666666666667 0.83333333333333 0.74924703706296 126.904 > 97 127 I 0.66666666666667 0.83333333333333 0.74924703706296 126.904 > 97 128 I 0.91666666666667 0.83333333333333 0.74924703706296 126.904 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.1532900 -0.0000000 0.0000000 -0.0000000 7.1532900 0.0000000 0.0000000 0.0000000 7.1937858 -------------------------- Born effective charges -------------------------- 1 Cu 1.0026013 -0.0000000 0.0000000 -0.0000000 1.0026013 0.0000000 0.0000000 0.0000000 1.0429109 2 Cu 1.0026013 -0.0000000 0.0000000 -0.0000000 1.0026013 0.0000000 0.0000000 0.0000000 1.0429109 3 I -1.0026013 0.0000000 0.0000000 0.0000000 -1.0026013 0.0000000 0.0000000 0.0000000 -1.0429109 4 I -1.0026013 0.0000000 0.0000000 0.0000000 -1.0026013 0.0000000 0.0000000 0.0000000 -1.0429109 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000029 (xzx) -0.00000029 (xzx) -0.00000029 (xxz) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 18:32:58]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 18:32:59]------------------------- 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): [4 4 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3mc (186) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.175399543208685 0.000000000000000 0.000000000000000 b -2.087699771604342 3.616002075368662 0.000000000000000 c 0.000000000000000 0.000000000000000 6.885862750000000 Atomic positions (fractional): 1 Cu 0.33333333333333 0.66666666666667 0.62650592587408 63.546 2 Cu 0.66666666666667 0.33333333333333 0.12650592587408 63.546 3 I 0.33333333333333 0.66666666666667 0.99849407412592 126.904 4 I 0.66666666666667 0.33333333333333 0.49849407412592 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.701598172834739 0.000000000000000 0.000000000000000 b -8.350799086417368 14.464008301474648 0.000000000000000 c 0.000000000000000 0.000000000000000 13.771725500000001 Atomic positions (fractional): 1 Cu 0.08333333333333 0.16666666666667 0.31325296293704 63.546 > 1 2 Cu 0.33333333333333 0.16666666666667 0.31325296293704 63.546 > 1 3 Cu 0.58333333333333 0.16666666666667 0.31325296293704 63.546 > 1 4 Cu 0.83333333333333 0.16666666666667 0.31325296293704 63.546 > 1 5 Cu 0.08333333333333 0.41666666666667 0.31325296293704 63.546 > 1 6 Cu 0.33333333333333 0.41666666666667 0.31325296293704 63.546 > 1 7 Cu 0.58333333333333 0.41666666666667 0.31325296293704 63.546 > 1 8 Cu 0.83333333333333 0.41666666666667 0.31325296293704 63.546 > 1 9 Cu 0.08333333333333 0.66666666666667 0.31325296293704 63.546 > 1 10 Cu 0.33333333333333 0.66666666666667 0.31325296293704 63.546 > 1 11 Cu 0.58333333333333 0.66666666666667 0.31325296293704 63.546 > 1 12 Cu 0.83333333333333 0.66666666666667 0.31325296293704 63.546 > 1 13 Cu 0.08333333333333 0.91666666666667 0.31325296293704 63.546 > 1 14 Cu 0.33333333333333 0.91666666666667 0.31325296293704 63.546 > 1 15 Cu 0.58333333333333 0.91666666666667 0.31325296293704 63.546 > 1 16 Cu 0.83333333333333 0.91666666666667 0.31325296293704 63.546 > 1 17 Cu 0.08333333333333 0.16666666666667 0.81325296293704 63.546 > 1 18 Cu 0.33333333333333 0.16666666666667 0.81325296293704 63.546 > 1 19 Cu 0.58333333333333 0.16666666666667 0.81325296293704 63.546 > 1 20 Cu 0.83333333333333 0.16666666666667 0.81325296293704 63.546 > 1 21 Cu 0.08333333333333 0.41666666666667 0.81325296293704 63.546 > 1 22 Cu 0.33333333333333 0.41666666666667 0.81325296293704 63.546 > 1 23 Cu 0.58333333333333 0.41666666666667 0.81325296293704 63.546 > 1 24 Cu 0.83333333333333 0.41666666666667 0.81325296293704 63.546 > 1 25 Cu 0.08333333333333 0.66666666666667 0.81325296293704 63.546 > 1 26 Cu 0.33333333333333 0.66666666666667 0.81325296293704 63.546 > 1 27 Cu 0.58333333333333 0.66666666666667 0.81325296293704 63.546 > 1 28 Cu 0.83333333333333 0.66666666666667 0.81325296293704 63.546 > 1 29 Cu 0.08333333333333 0.91666666666667 0.81325296293704 63.546 > 1 30 Cu 0.33333333333333 0.91666666666667 0.81325296293704 63.546 > 1 31 Cu 0.58333333333333 0.91666666666667 0.81325296293704 63.546 > 1 32 Cu 0.83333333333333 0.91666666666667 0.81325296293704 63.546 > 1 33 Cu 0.16666666666667 0.08333333333333 0.06325296293704 63.546 > 33 34 Cu 0.41666666666667 0.08333333333333 0.06325296293704 63.546 > 33 35 Cu 0.66666666666667 0.08333333333333 0.06325296293704 63.546 > 33 36 Cu 0.91666666666667 0.08333333333333 0.06325296293704 63.546 > 33 37 Cu 0.16666666666667 0.33333333333333 0.06325296293704 63.546 > 33 38 Cu 0.41666666666667 0.33333333333333 0.06325296293704 63.546 > 33 39 Cu 0.66666666666667 0.33333333333333 0.06325296293704 63.546 > 33 40 Cu 0.91666666666667 0.33333333333333 0.06325296293704 63.546 > 33 41 Cu 0.16666666666667 0.58333333333333 0.06325296293704 63.546 > 33 42 Cu 0.41666666666667 0.58333333333333 0.06325296293704 63.546 > 33 43 Cu 0.66666666666667 0.58333333333333 0.06325296293704 63.546 > 33 44 Cu 0.91666666666667 0.58333333333333 0.06325296293704 63.546 > 33 45 Cu 0.16666666666667 0.83333333333333 0.06325296293704 63.546 > 33 46 Cu 0.41666666666667 0.83333333333333 0.06325296293704 63.546 > 33 47 Cu 0.66666666666667 0.83333333333333 0.06325296293704 63.546 > 33 48 Cu 0.91666666666667 0.83333333333333 0.06325296293704 63.546 > 33 49 Cu 0.16666666666667 0.08333333333333 0.56325296293704 63.546 > 33 50 Cu 0.41666666666667 0.08333333333333 0.56325296293704 63.546 > 33 51 Cu 0.66666666666667 0.08333333333333 0.56325296293704 63.546 > 33 52 Cu 0.91666666666667 0.08333333333333 0.56325296293704 63.546 > 33 53 Cu 0.16666666666667 0.33333333333333 0.56325296293704 63.546 > 33 54 Cu 0.41666666666667 0.33333333333333 0.56325296293704 63.546 > 33 55 Cu 0.66666666666667 0.33333333333333 0.56325296293704 63.546 > 33 56 Cu 0.91666666666667 0.33333333333333 0.56325296293704 63.546 > 33 57 Cu 0.16666666666667 0.58333333333333 0.56325296293704 63.546 > 33 58 Cu 0.41666666666667 0.58333333333333 0.56325296293704 63.546 > 33 59 Cu 0.66666666666667 0.58333333333333 0.56325296293704 63.546 > 33 60 Cu 0.91666666666667 0.58333333333333 0.56325296293704 63.546 > 33 61 Cu 0.16666666666667 0.83333333333333 0.56325296293704 63.546 > 33 62 Cu 0.41666666666667 0.83333333333333 0.56325296293704 63.546 > 33 63 Cu 0.66666666666667 0.83333333333333 0.56325296293704 63.546 > 33 64 Cu 0.91666666666667 0.83333333333333 0.56325296293704 63.546 > 33 65 I 0.08333333333333 0.16666666666667 0.49924703706296 126.904 > 65 66 I 0.33333333333333 0.16666666666667 0.49924703706296 126.904 > 65 67 I 0.58333333333333 0.16666666666667 0.49924703706296 126.904 > 65 68 I 0.83333333333333 0.16666666666667 0.49924703706296 126.904 > 65 69 I 0.08333333333333 0.41666666666667 0.49924703706296 126.904 > 65 70 I 0.33333333333333 0.41666666666667 0.49924703706296 126.904 > 65 71 I 0.58333333333333 0.41666666666667 0.49924703706296 126.904 > 65 72 I 0.83333333333333 0.41666666666667 0.49924703706296 126.904 > 65 73 I 0.08333333333333 0.66666666666667 0.49924703706296 126.904 > 65 74 I 0.33333333333333 0.66666666666667 0.49924703706296 126.904 > 65 75 I 0.58333333333333 0.66666666666667 0.49924703706296 126.904 > 65 76 I 0.83333333333333 0.66666666666667 0.49924703706296 126.904 > 65 77 I 0.08333333333333 0.91666666666667 0.49924703706296 126.904 > 65 78 I 0.33333333333333 0.91666666666667 0.49924703706296 126.904 > 65 79 I 0.58333333333333 0.91666666666667 0.49924703706296 126.904 > 65 80 I 0.83333333333333 0.91666666666667 0.49924703706296 126.904 > 65 81 I 0.08333333333333 0.16666666666667 0.99924703706296 126.904 > 65 82 I 0.33333333333333 0.16666666666667 0.99924703706296 126.904 > 65 83 I 0.58333333333333 0.16666666666667 0.99924703706296 126.904 > 65 84 I 0.83333333333333 0.16666666666667 0.99924703706296 126.904 > 65 85 I 0.08333333333333 0.41666666666667 0.99924703706296 126.904 > 65 86 I 0.33333333333333 0.41666666666667 0.99924703706296 126.904 > 65 87 I 0.58333333333333 0.41666666666667 0.99924703706296 126.904 > 65 88 I 0.83333333333333 0.41666666666667 0.99924703706296 126.904 > 65 89 I 0.08333333333333 0.66666666666667 0.99924703706296 126.904 > 65 90 I 0.33333333333333 0.66666666666667 0.99924703706296 126.904 > 65 91 I 0.58333333333333 0.66666666666667 0.99924703706296 126.904 > 65 92 I 0.83333333333333 0.66666666666667 0.99924703706296 126.904 > 65 93 I 0.08333333333333 0.91666666666667 0.99924703706296 126.904 > 65 94 I 0.33333333333333 0.91666666666667 0.99924703706296 126.904 > 65 95 I 0.58333333333333 0.91666666666667 0.99924703706296 126.904 > 65 96 I 0.83333333333333 0.91666666666667 0.99924703706296 126.904 > 65 97 I 0.16666666666667 0.08333333333333 0.24924703706296 126.904 > 97 98 I 0.41666666666667 0.08333333333333 0.24924703706296 126.904 > 97 99 I 0.66666666666667 0.08333333333333 0.24924703706296 126.904 > 97 100 I 0.91666666666667 0.08333333333333 0.24924703706296 126.904 > 97 101 I 0.16666666666667 0.33333333333333 0.24924703706296 126.904 > 97 102 I 0.41666666666667 0.33333333333333 0.24924703706296 126.904 > 97 103 I 0.66666666666667 0.33333333333333 0.24924703706296 126.904 > 97 104 I 0.91666666666667 0.33333333333333 0.24924703706296 126.904 > 97 105 I 0.16666666666667 0.58333333333333 0.24924703706296 126.904 > 97 106 I 0.41666666666667 0.58333333333333 0.24924703706296 126.904 > 97 107 I 0.66666666666667 0.58333333333333 0.24924703706296 126.904 > 97 108 I 0.91666666666667 0.58333333333333 0.24924703706296 126.904 > 97 109 I 0.16666666666667 0.83333333333333 0.24924703706296 126.904 > 97 110 I 0.41666666666667 0.83333333333333 0.24924703706296 126.904 > 97 111 I 0.66666666666667 0.83333333333333 0.24924703706296 126.904 > 97 112 I 0.91666666666667 0.83333333333333 0.24924703706296 126.904 > 97 113 I 0.16666666666667 0.08333333333333 0.74924703706296 126.904 > 97 114 I 0.41666666666667 0.08333333333333 0.74924703706296 126.904 > 97 115 I 0.66666666666667 0.08333333333333 0.74924703706296 126.904 > 97 116 I 0.91666666666667 0.08333333333333 0.74924703706296 126.904 > 97 117 I 0.16666666666667 0.33333333333333 0.74924703706296 126.904 > 97 118 I 0.41666666666667 0.33333333333333 0.74924703706296 126.904 > 97 119 I 0.66666666666667 0.33333333333333 0.74924703706296 126.904 > 97 120 I 0.91666666666667 0.33333333333333 0.74924703706296 126.904 > 97 121 I 0.16666666666667 0.58333333333333 0.74924703706296 126.904 > 97 122 I 0.41666666666667 0.58333333333333 0.74924703706296 126.904 > 97 123 I 0.66666666666667 0.58333333333333 0.74924703706296 126.904 > 97 124 I 0.91666666666667 0.58333333333333 0.74924703706296 126.904 > 97 125 I 0.16666666666667 0.83333333333333 0.74924703706296 126.904 > 97 126 I 0.41666666666667 0.83333333333333 0.74924703706296 126.904 > 97 127 I 0.66666666666667 0.83333333333333 0.74924703706296 126.904 > 97 128 I 0.91666666666667 0.83333333333333 0.74924703706296 126.904 > 97 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 7.1532900 -0.0000000 0.0000000 -0.0000000 7.1532900 0.0000000 0.0000000 0.0000000 7.1937858 -------------------------- Born effective charges -------------------------- 1 Cu 1.0026013 -0.0000000 0.0000000 -0.0000000 1.0026013 0.0000000 0.0000000 0.0000000 1.0429109 2 Cu 1.0026013 -0.0000000 0.0000000 -0.0000000 1.0026013 0.0000000 0.0000000 0.0000000 1.0429109 3 I -1.0026013 0.0000000 0.0000000 0.0000000 -1.0026013 0.0000000 0.0000000 0.0000000 -1.0429109 4 I -1.0026013 0.0000000 0.0000000 0.0000000 -1.0026013 0.0000000 0.0000000 0.0000000 -1.0429109 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000029 (xzx) -0.00000029 (xzx) -0.00000029 (xxz) 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 ... [ 14 14 7 ] fc3-r2q-transformation over three atoms: True --------------------------- Phonon calculations ---------------------------- Use NAC by Gonze et al. (no real space sum in current implementation) PRB 50, 13035(R) (1994), PRB 55, 10355 (1997) G-cutoff distance: 0.88, Number of G-points: 313, Lambda: 0.25 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/96) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 96 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 1.084 ( 0.000 0.000 0.000) 0.000 1.084 ( -0.000 -0.000 -0.000) 0.000 3.473 ( 0.000 -0.000 -0.000) 0.000 4.486 ( -0.000 -0.000 -0.000) 0.000 4.486 ( -0.000 0.000 -0.000) 0.000 4.486 ( 0.000 0.000 0.000) 0.000 4.486 ( 0.000 0.000 -0.000) 0.000 4.544 ( -0.000 -0.000 -0.000) 0.000 4.722 ( 0.000 0.000 -0.000) 0.000 ======================= Grid point 1 (2/96) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.288 ( 12.559 7.251 0.000) 14.502 0.293 ( 12.435 7.179 0.000) 14.359 0.689 ( 29.560 17.067 0.000) 34.133 1.125 ( 3.399 1.963 0.000) 3.925 1.256 ( 13.393 7.732 0.000) 15.465 3.456 ( -1.296 -0.748 0.000) 1.496 4.480 ( -0.697 -0.403 0.000) 0.805 4.486 ( 0.339 0.196 0.000) 0.391 4.500 ( 0.813 0.469 0.000) 0.939 4.548 ( 0.376 0.217 0.000) 0.434 4.740 ( 1.549 0.894 0.000) 1.788 4.788 ( -0.218 -0.126 0.000) 0.252 ======================= Grid point 2 (3/96) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.558 ( 10.688 6.171 0.000) 12.341 0.571 ( 12.180 7.032 0.000) 14.064 1.225 ( 5.069 2.927 0.000) 5.854 1.335 ( 26.725 15.429 0.000) 30.859 1.602 ( 15.263 8.812 0.000) 17.624 3.430 ( -0.417 -0.241 0.000) 0.481 4.448 ( -2.221 -1.282 0.000) 2.564 4.502 ( 1.006 0.581 0.000) 1.162 4.522 ( 1.113 0.642 0.000) 1.285 4.561 ( 0.780 0.450 0.000) 0.900 4.775 ( 0.056 0.032 0.000) 0.065 4.797 ( 2.352 1.358 0.000) 2.716 ======================= Grid point 3 (4/96) ======================= q-point: ( 0.21 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.773 ( 8.043 4.643 0.000) 9.287 0.842 ( 11.527 6.655 0.000) 13.311 1.341 ( 4.781 2.760 0.000) 5.521 1.851 ( 13.373 7.721 0.000) 15.442 1.944 ( 19.087 11.020 0.000) 22.040 3.460 ( 3.361 1.941 0.000) 3.881 4.378 ( -3.834 -2.214 0.000) 4.428 4.528 ( 1.237 0.714 0.000) 1.428 4.552 ( 1.541 0.890 0.000) 1.779 4.583 ( 1.076 0.621 0.000) 1.242 4.767 ( -0.552 -0.319 0.000) 0.637 4.861 ( 2.788 1.610 0.000) 3.219 ======================= Grid point 4 (5/96) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.923 ( 5.097 2.943 0.000) 5.885 1.094 ( 10.465 6.042 0.000) 12.084 1.434 ( 3.270 1.888 0.000) 3.775 2.004 ( 2.694 1.555 0.000) 3.111 2.403 ( 18.820 10.866 0.000) 21.732 3.585 ( 7.353 4.245 0.000) 8.490 4.274 ( -5.242 -3.026 0.000) 6.053 4.557 ( 1.256 0.725 0.000) 1.450 4.591 ( 1.755 1.013 0.000) 2.026 4.605 ( 0.664 0.383 0.000) 0.767 4.753 ( -0.569 -0.328 0.000) 0.657 4.915 ( 1.848 1.067 0.000) 2.133 ======================= Grid point 5 (6/96) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.009 ( 2.592 1.497 0.000) 2.993 1.314 ( 8.713 5.030 0.000) 10.061 1.488 ( 1.529 0.883 0.000) 1.766 2.015 ( -1.089 -0.629 0.000) 1.258 2.785 ( 14.588 8.422 0.000) 16.844 3.773 ( 8.521 4.920 0.000) 9.839 4.143 ( -6.145 -3.548 0.000) 7.095 4.585 ( 1.172 0.677 0.000) 1.353 4.599 ( -1.416 -0.818 0.000) 1.635 4.628 ( 1.473 0.851 0.000) 1.701 4.749 ( 0.343 0.198 0.000) 0.396 4.945 ( 0.794 0.458 0.000) 0.917 ======================= Grid point 6 (7/96) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.048 ( 0.959 0.554 0.000) 1.108 1.481 ( 5.605 3.236 0.000) 6.472 1.509 ( 0.397 0.229 0.000) 0.458 1.976 ( -1.904 -1.099 0.000) 2.199 3.061 ( 9.159 5.288 0.000) 10.576 3.950 ( 6.447 3.722 0.000) 7.445 4.007 ( -5.225 -3.017 0.000) 6.034 4.546 ( -2.763 -1.595 0.000) 3.190 4.607 ( 0.747 0.431 0.000) 0.862 4.655 ( 0.829 0.479 0.000) 0.958 4.764 ( 0.764 0.441 0.000) 0.882 4.954 ( 0.162 0.094 0.000) 0.187 ======================= Grid point 7 (8/96) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 128 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.058 ( 0.000 0.000 0.000) 0.000 1.512 ( 0.000 0.000 0.000) 0.000 1.550 ( 0.000 0.000 0.000) 0.000 1.950 ( 0.000 0.000 0.000) 0.000 3.173 ( 0.000 0.000 0.000) 0.000 3.939 ( 0.000 0.000 0.000) 0.000 4.032 ( 0.000 0.000 0.000) 0.000 4.506 ( 0.000 0.000 0.000) 0.000 4.616 ( 0.000 0.000 0.000) 0.000 4.665 ( 0.000 0.000 0.000) 0.000 4.774 ( 0.000 0.000 0.000) 0.000 4.956 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 16 (9/96) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.495 ( 7.062 12.232 0.000) 14.125 0.528 ( 7.908 13.698 0.000) 15.817 1.127 ( 11.937 20.676 0.000) 23.874 1.246 ( 7.136 12.359 0.000) 14.271 1.486 ( 8.391 14.534 0.000) 16.783 3.434 ( -0.584 -1.012 0.000) 1.169 4.457 ( -1.092 -1.891 0.000) 2.183 4.506 ( 0.589 1.020 0.000) 1.178 4.512 ( 0.791 1.370 0.000) 1.582 4.557 ( 0.383 0.663 0.000) 0.766 4.775 ( 1.439 2.492 0.000) 2.878 4.780 ( -0.342 -0.593 0.000) 0.684 ======================= Grid point 17 (10/96) ======================= q-point: ( 0.14 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.745 ( 8.933 10.059 0.000) 13.453 0.780 ( 4.922 14.220 0.000) 15.047 1.334 ( 3.619 7.753 0.000) 8.556 1.692 ( 16.802 18.960 0.000) 25.334 1.777 ( 11.859 8.669 0.000) 14.690 3.436 ( 1.294 1.119 0.000) 1.711 4.401 ( -2.415 -3.453 0.000) 4.214 4.531 ( 0.501 2.156 0.000) 2.213 4.547 ( 1.171 1.910 0.000) 2.240 4.574 ( 0.755 0.812 0.000) 1.108 4.766 ( -0.316 -0.898 0.000) 0.952 4.832 ( 2.362 2.085 0.000) 3.151 ======================= Grid point 18 (11/96) ======================= q-point: ( 0.21 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.969 ( 1.086 13.656 0.000) 13.699 0.993 ( 9.226 8.400 0.000) 12.477 1.455 ( 1.201 8.224 0.000) 8.311 1.950 ( 5.635 2.848 0.000) 6.314 2.195 ( 18.609 13.374 0.000) 22.916 3.515 ( 5.222 3.653 0.000) 6.373 4.308 ( -3.594 -4.813 0.000) 6.007 4.561 ( 0.251 2.375 0.000) 2.389 4.587 ( 0.978 2.430 0.000) 2.619 4.596 ( 0.890 0.646 0.000) 1.099 4.754 ( -0.133 -0.701 0.000) 0.713 4.889 ( 2.230 1.123 0.000) 2.496 ======================= Grid point 19 (12/96) ======================= q-point: ( 0.29 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.088 ( -2.190 12.418 0.000) 12.609 1.223 ( 8.590 6.892 0.000) 11.013 1.535 ( -1.122 7.755 0.000) 7.835 2.012 ( 0.935 -0.531 0.000) 1.075 2.609 ( 16.151 9.875 0.000) 18.931 3.675 ( 8.087 4.683 0.000) 9.345 4.185 ( -4.529 -5.898 0.000) 7.436 4.577 ( -0.103 -0.093 0.000) 0.139 4.611 ( 0.024 1.439 0.000) 1.439 4.630 ( 0.917 2.544 0.000) 2.704 4.751 ( -0.041 0.499 0.000) 0.501 4.927 ( 1.499 0.183 0.000) 1.510 ======================= Grid point 20 (13/96) ======================= q-point: ( 0.36 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.146 ( -4.356 11.216 0.000) 12.032 1.416 ( 6.690 5.283 0.000) 8.525 1.573 ( -2.862 7.093 0.000) 7.649 1.993 ( -1.406 -1.473 0.000) 2.036 2.936 ( 12.163 6.911 0.000) 13.989 3.861 ( 8.045 3.942 0.000) 8.959 4.044 ( -4.575 -6.232 0.000) 7.731 4.547 ( -1.172 -3.440 0.000) 3.634 4.628 ( -0.537 2.630 0.000) 2.685 4.667 ( 0.215 2.809 0.000) 2.818 4.762 ( 0.249 1.106 0.000) 1.134 4.945 ( 0.761 -0.390 0.000) 0.855 ======================= Grid point 21 (14/96) ======================= q-point: ( 0.43 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.167 ( -5.508 10.446 0.000) 11.809 1.541 ( 2.574 2.995 0.000) 3.949 1.583 ( -3.704 6.556 0.000) 7.530 1.954 ( -1.055 -1.065 0.000) 1.499 3.142 ( 5.193 2.946 0.000) 5.970 3.924 ( -1.196 -5.158 0.000) 5.295 4.004 ( 3.799 1.800 0.000) 4.204 4.490 ( -0.521 -3.448 0.000) 3.487 4.641 ( -0.967 2.372 0.000) 2.562 4.690 ( -0.930 2.805 0.000) 2.955 4.775 ( 0.153 0.571 0.000) 0.591 4.949 ( 0.403 -0.595 0.000) 0.719 ======================= Grid point 31 (15/96) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.954 ( 6.251 10.826 0.000) 12.501 1.052 ( 7.090 12.280 0.000) 14.179 1.498 ( 4.672 8.092 0.000) 9.344 1.907 ( 3.172 5.494 0.000) 6.344 2.086 ( 11.171 19.349 0.000) 22.342 3.487 ( 2.463 4.266 0.000) 4.926 4.313 ( -3.069 -5.315 0.000) 6.137 4.579 ( 1.478 2.560 0.000) 2.956 4.588 ( 0.582 1.007 0.000) 1.163 4.591 ( 1.054 1.826 0.000) 2.108 4.750 ( -0.242 -0.419 0.000) 0.484 4.873 ( 1.084 1.878 0.000) 2.168 ======================= Grid point 32 (16/96) ======================= q-point: ( 0.21 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.166 ( 6.967 8.982 0.000) 11.368 1.247 ( 1.002 13.164 0.000) 13.202 1.634 ( 0.988 9.067 0.000) 9.121 1.980 ( 2.052 0.400 0.000) 2.091 2.460 ( 13.318 13.556 0.000) 19.004 3.605 ( 5.486 5.479 0.000) 7.753 4.193 ( -3.982 -6.701 0.000) 7.795 4.581 ( -0.924 -1.263 0.000) 1.565 4.626 ( 0.718 2.812 0.000) 2.902 4.637 ( 1.450 2.392 0.000) 2.797 4.758 ( 0.435 1.417 0.000) 1.482 4.904 ( 1.235 0.395 0.000) 1.297 ======================= Grid point 33 (17/96) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 400 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.352 ( -3.248 12.816 0.000) 13.222 1.364 ( 6.186 7.558 0.000) 9.767 1.710 ( -2.484 9.182 0.000) 9.512 1.990 ( 0.068 -1.588 0.000) 1.590 2.799 ( 12.576 9.631 0.000) 15.840 3.769 ( 7.203 4.837 0.000) 8.677 4.050 ( -4.304 -7.458 0.000) 8.610 4.535 ( -1.548 -3.510 0.000) 3.836 4.661 ( 0.048 2.825 0.000) 2.825 4.685 ( 1.139 2.772 0.000) 2.997 4.781 ( -0.267 2.430 0.000) 2.445 4.923 ( 1.117 -0.662 0.000) 1.298 ======================= Grid point 34 (18/96) ======================= q-point: ( 0.36 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.395 ( -5.984 12.591 0.000) 13.941 1.523 ( 4.271 5.525 0.000) 6.983 1.737 ( -4.412 8.676 0.000) 9.733 1.961 ( -1.143 -1.748 0.000) 2.088 3.065 ( 8.604 6.425 0.000) 10.738 3.877 ( -0.110 -6.237 0.000) 6.238 3.960 ( 3.142 2.053 0.000) 3.754 4.466 ( -1.430 -4.604 0.000) 4.821 4.683 ( -0.759 2.718 0.000) 2.822 4.729 ( 0.500 3.154 0.000) 3.193 4.790 ( -1.045 1.738 0.000) 2.027 4.930 ( 0.841 -1.055 0.000) 1.349 ======================= Grid point 35 (19/96) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.404 ( -7.036 12.187 0.000) 14.073 1.594 ( -1.353 2.344 0.000) 2.706 1.739 ( -4.794 8.304 0.000) 9.588 1.937 ( 0.318 -0.550 0.000) 0.635 3.182 ( -0.735 1.273 0.000) 1.470 3.815 ( 3.220 -5.577 0.000) 6.440 4.027 ( -0.354 0.613 0.000) 0.708 4.422 ( 1.867 -3.233 0.000) 3.733 4.689 ( -1.405 2.434 0.000) 2.810 4.754 ( -1.914 3.315 0.000) 3.827 4.785 ( -0.332 0.575 0.000) 0.663 4.932 ( 0.642 -1.112 0.000) 1.284 ======================= Grid point 47 (20/96) ======================= q-point: ( 0.21 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.343 ( 5.091 8.818 0.000) 10.182 1.459 ( 4.644 8.044 0.000) 9.289 1.789 ( 3.615 6.262 0.000) 7.231 1.974 ( -0.502 -0.870 0.000) 1.005 2.732 ( 7.986 13.832 0.000) 15.972 3.728 ( 3.952 6.845 0.000) 7.904 4.049 ( -4.405 -7.630 0.000) 8.811 4.531 ( -2.013 -3.487 0.000) 4.026 4.677 ( 1.260 2.183 0.000) 2.520 4.687 ( 1.481 2.565 0.000) 2.961 4.801 ( 1.386 2.401 0.000) 2.772 4.906 ( -0.098 -0.170 0.000) 0.196 ======================= Grid point 48 (21/96) ======================= q-point: ( 0.29 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.508 ( 4.637 7.255 0.000) 8.610 1.569 ( -0.919 8.472 0.000) 8.522 1.871 ( -1.656 6.611 0.000) 6.815 1.951 ( 0.058 -2.207 0.000) 2.207 2.996 ( 8.484 10.410 0.000) 13.429 3.826 ( -0.011 -3.176 0.000) 3.176 3.939 ( 0.886 0.514 0.000) 1.025 4.451 ( -2.337 -4.860 0.000) 5.393 4.714 ( 0.246 2.588 0.000) 2.600 4.735 ( 1.519 2.162 0.000) 2.642 4.834 ( -0.766 2.664 0.000) 2.772 4.902 ( 0.805 -1.401 0.000) 1.615 ======================= Grid point 49 (22/96) ======================= q-point: ( 0.36 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 392 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.606 ( -3.608 8.225 0.000) 8.981 1.633 ( 1.222 5.578 0.000) 5.710 1.865 ( -2.704 1.866 0.000) 3.286 1.945 ( -1.017 2.066 0.000) 2.302 3.199 ( 3.500 7.180 0.000) 7.987 3.718 ( -0.571 -8.613 0.000) 8.632 4.035 ( 2.286 4.512 0.000) 5.058 4.368 ( -0.598 -5.217 0.000) 5.251 4.734 ( -0.881 2.451 0.000) 2.604 4.782 ( 1.414 2.006 0.000) 2.455 4.824 ( -2.432 1.387 0.000) 2.800 4.907 ( 0.989 -1.184 0.000) 1.542 ======================= Grid point 62 (23/96) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 228 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.644 ( 3.653 6.327 0.000) 7.306 1.686 ( 2.050 3.550 0.000) 4.099 1.886 ( -1.943 -3.366 0.000) 3.887 1.978 ( 1.821 3.153 0.000) 3.641 3.204 ( 5.917 10.249 0.000) 11.834 3.678 ( -5.060 -8.764 0.000) 10.119 4.038 ( 3.697 6.403 0.000) 7.394 4.346 ( -3.211 -5.562 0.000) 6.423 4.765 ( 1.269 2.198 0.000) 2.538 4.768 ( 0.854 1.479 0.000) 1.707 4.866 ( -0.977 -1.692 0.000) 1.953 4.879 ( 0.704 1.219 0.000) 1.407 ======================= Grid point 63 (24/96) ======================= q-point: ( 0.36 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 224 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.718 ( -1.829 3.167 0.000) 3.657 1.730 ( -2.195 3.802 0.000) 4.390 1.832 ( 1.717 -2.973 0.000) 3.433 2.009 ( -1.240 2.149 0.000) 2.481 3.345 ( -4.098 7.099 0.000) 8.197 3.552 ( 4.371 -7.571 0.000) 8.742 4.129 ( -2.683 4.647 0.000) 5.366 4.265 ( 2.796 -4.844 0.000) 5.593 4.778 ( -1.085 1.879 0.000) 2.170 4.806 ( -0.193 0.334 0.000) 0.385 4.832 ( 0.254 -0.440 0.000) 0.508 4.893 ( 0.066 -0.114 0.000) 0.132 ======================= Grid point 211 (25/96) ======================= q-point: ( 0.00 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.309 ( 0.000 0.000 14.352) 14.352 0.309 ( 0.000 -0.000 14.352) 14.352 0.782 ( 0.000 -0.000 36.988) 36.988 1.074 ( 0.000 0.000 -1.002) 1.002 1.074 ( -0.000 0.000 -1.002) 1.002 3.386 ( 0.000 0.000 -8.332) 8.332 4.484 ( -0.000 0.000 -0.156) 0.156 4.484 ( -0.000 0.000 -0.156) 0.156 4.485 ( 0.000 -0.000 -0.158) 0.158 4.485 ( -0.000 -0.000 -0.158) 0.158 4.731 ( 0.000 0.000 0.762) 0.762 4.859 ( 0.000 -0.000 -0.706) 0.706 ======================= Grid point 212 (26/96) ======================= q-point: ( 0.07 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.427 ( 8.611 4.971 10.484) 14.449 0.549 ( 9.565 5.522 14.716) 18.400 0.962 ( 18.201 10.509 24.873) 32.564 1.116 ( 3.446 1.989 -0.945) 4.089 1.243 ( 13.198 7.620 -1.189) 15.286 3.369 ( -1.273 -0.735 -8.290) 8.419 4.484 ( 0.180 0.104 0.066) 0.219 4.484 ( 0.038 0.022 0.214) 0.218 4.498 ( 0.783 0.452 -0.227) 0.932 4.530 ( 2.402 1.387 -0.798) 2.886 4.743 ( 1.156 0.667 0.312) 1.370 4.812 ( -2.445 -1.412 0.674) 2.903 ======================= Grid point 213 (27/96) ======================= q-point: ( 0.14 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.640 ( 9.399 5.427 7.127) 12.984 0.727 ( 7.612 4.395 13.436) 16.055 1.218 ( 5.121 2.956 -0.827) 5.971 1.474 ( 23.381 13.499 12.070) 29.573 1.587 ( 15.590 9.001 -0.928) 18.026 3.345 ( -0.239 -0.138 -8.009) 8.014 4.467 ( -1.562 -0.902 1.665) 2.454 4.500 ( 1.000 0.577 -0.108) 1.159 4.519 ( 1.092 0.631 -0.278) 1.292 4.572 ( 1.404 0.811 0.800) 1.808 4.769 ( -0.617 -0.356 -0.597) 0.929 4.789 ( 1.971 1.138 -0.637) 2.363 ======================= Grid point 214 (28/96) ======================= q-point: ( 0.21 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.836 ( 7.500 4.330 5.588) 10.307 0.910 ( 8.297 4.790 6.889) 11.800 1.334 ( 4.819 2.782 -0.725) 5.612 1.841 ( 9.007 5.200 -2.253) 10.642 2.028 ( 21.741 12.552 8.708) 26.571 3.381 ( 3.759 2.170 -7.202) 8.409 4.412 ( -3.164 -1.826 2.957) 4.700 4.527 ( 1.244 0.718 -0.111) 1.441 4.549 ( 1.511 0.873 -0.327) 1.775 4.600 ( 1.164 0.672 1.532) 2.038 4.752 ( -0.708 -0.409 -1.334) 1.564 4.847 ( 2.551 1.473 -1.284) 3.213 ======================= Grid point 215 (29/96) ======================= q-point: ( 0.29 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.977 ( 4.852 2.802 4.872) 7.425 1.100 ( 8.207 4.738 1.347) 9.572 1.428 ( 3.304 1.908 -0.645) 3.869 1.958 ( 1.817 1.049 -4.419) 4.892 2.490 ( 18.431 10.641 7.751) 22.650 3.515 ( 7.564 4.367 -6.176) 10.697 4.326 ( -4.332 -2.501 4.391) 6.656 4.556 ( 1.272 0.735 -0.081) 1.471 4.586 ( 1.723 0.995 -0.395) 2.029 4.623 ( 0.618 0.357 1.640) 1.789 4.736 ( -0.659 -0.380 -1.631) 1.799 4.896 ( 1.681 0.970 -1.718) 2.592 ======================= Grid point 216 (30/96) ======================= q-point: ( 0.36 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.060 ( 2.489 1.437 4.547) 5.379 1.278 ( 7.246 4.184 -2.659) 8.779 1.483 ( 1.569 0.906 -0.564) 1.897 1.953 ( -1.717 -0.991 -5.916) 6.239 2.864 ( 14.311 8.262 7.016) 17.953 3.698 ( 7.775 4.489 -5.872) 10.728 4.222 ( -4.514 -2.606 6.254) 8.141 4.584 ( 1.180 0.681 -0.052) 1.364 4.612 ( -1.771 -1.022 1.151) 2.346 4.624 ( 1.454 0.839 -0.450) 1.738 4.733 ( 0.546 0.316 -1.451) 1.582 4.923 ( 0.688 0.397 -2.016) 2.167 ======================= Grid point 217 (31/96) ======================= q-point: ( 0.43 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.097 ( 0.924 0.534 4.417) 4.544 1.420 ( 4.881 2.818 -5.091) 7.595 1.504 ( 0.428 0.247 -0.484) 0.692 1.902 ( -2.256 -1.302 -6.877) 7.354 3.134 ( 8.961 5.174 6.504) 12.222 3.833 ( 3.259 1.882 -6.009) 7.090 4.145 ( -1.519 -0.877 7.874) 8.066 4.552 ( -2.944 -1.700 0.509) 3.437 4.607 ( 0.747 0.431 -0.043) 0.864 4.650 ( 0.821 0.474 -0.479) 1.063 4.752 ( 0.860 0.496 -1.061) 1.453 4.931 ( 0.119 0.069 -2.176) 2.181 ======================= Grid point 218 (32/96) ======================= q-point: (-0.50 0.00 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.107 ( 0.000 0.000 4.383) 4.383 1.480 ( 0.000 0.000 -5.914) 5.914 1.508 ( 0.000 0.000 -0.447) 0.447 1.872 ( -0.000 -0.000 -7.212) 7.212 3.244 ( 0.000 0.000 6.267) 6.267 3.857 ( -0.000 -0.000 -4.654) 4.654 4.142 ( -0.000 -0.000 7.077) 7.077 4.511 ( 0.000 0.000 0.349) 0.349 4.616 ( 0.000 0.000 -0.043) 0.043 4.660 ( 0.000 0.000 -0.487) 0.487 4.763 ( 0.000 0.000 -0.984) 0.984 4.932 ( -0.000 -0.000 -2.221) 2.221 ======================= Grid point 227 (33/96) ======================= q-point: ( 0.07 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.614 ( 6.851 11.866 7.391) 15.568 0.674 ( 4.125 7.145 14.502) 16.685 1.193 ( 4.061 7.033 0.944) 8.176 1.337 ( 13.024 22.558 13.193) 29.199 1.471 ( 8.222 14.241 -1.414) 16.505 3.348 ( -0.518 -0.898 -8.130) 8.196 4.472 ( -0.705 -1.221 1.266) 1.895 4.505 ( 0.582 1.007 -0.087) 1.166 4.509 ( 0.771 1.335 -0.267) 1.564 4.562 ( 0.889 1.539 0.347) 1.811 4.770 ( 1.195 2.070 -0.355) 2.416 4.779 ( -0.878 -1.521 -0.215) 1.770 ======================= Grid point 228 (34/96) ======================= q-point: ( 0.14 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.822 ( 4.920 9.043 7.014) 12.457 0.853 ( 6.191 10.591 7.138) 14.193 1.331 ( 3.238 7.941 -0.508) 8.591 1.738 ( 9.278 9.027 -1.261) 13.006 1.822 ( 18.018 17.141 9.334) 26.563 3.354 ( 1.547 1.377 -7.580) 7.858 4.429 ( -1.920 -2.907 2.457) 4.264 4.529 ( 0.485 2.166 -0.114) 2.223 4.544 ( 1.201 1.850 -0.299) 2.226 4.589 ( 0.939 0.948 1.295) 1.860 4.754 ( -0.552 -1.054 -1.078) 1.606 4.821 ( 2.135 1.858 -1.059) 3.022 ======================= Grid point 229 (35/96) ======================= q-point: ( 0.21 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.992 ( 4.440 8.503 2.710) 9.968 1.045 ( 4.017 10.377 4.839) 12.134 1.450 ( 1.095 8.398 -0.543) 8.487 1.913 ( 4.841 2.020 -3.668) 6.401 2.287 ( 18.052 13.292 8.378) 23.932 3.441 ( 5.550 3.884 -6.603) 9.460 4.352 ( -2.933 -4.204 3.744) 6.348 4.559 ( 0.302 2.375 -0.109) 2.397 4.583 ( 0.945 2.381 -0.354) 2.586 4.614 ( 0.901 0.622 1.589) 1.930 4.738 ( -0.263 -0.736 -1.465) 1.661 4.872 ( 2.060 0.991 -1.545) 2.759 ======================= Grid point 230 (36/96) ======================= q-point: ( 0.29 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.122 ( -0.460 10.328 2.548) 10.647 1.210 ( 5.678 6.768 0.146) 8.835 1.530 ( -1.176 7.831 -0.544) 7.937 1.956 ( 0.275 -1.037 -5.295) 5.403 2.692 ( 15.744 9.780 7.427) 19.967 3.605 ( 7.938 4.510 -5.899) 10.870 4.250 ( -3.448 -5.011 5.368) 8.112 4.579 ( -0.011 0.357 0.076) 0.365 4.616 ( 0.140 1.542 0.086) 1.551 4.633 ( 0.553 1.687 0.724) 1.917 4.734 ( -0.032 0.660 -1.535) 1.671 4.907 ( 1.372 0.116 -1.881) 2.331 ======================= Grid point 231 (37/96) ======================= q-point: ( 0.36 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.185 ( -3.648 10.229 3.338) 11.361 1.365 ( 5.453 4.762 -3.827) 8.189 1.567 ( -2.869 7.038 -0.581) 7.622 1.924 ( -1.871 -1.767 -6.467) 6.960 3.013 ( 11.866 6.832 6.790) 15.283 3.771 ( 6.249 2.603 -6.052) 9.081 4.149 ( -2.126 -4.355 7.364) 8.816 4.551 ( -1.151 -3.609 0.304) 3.801 4.630 ( -0.628 2.578 0.170) 2.659 4.664 ( 0.066 2.777 -0.271) 2.791 4.748 ( 0.360 1.229 -1.193) 1.750 4.923 ( 0.680 -0.394 -2.084) 2.227 ======================= Grid point 232 (38/96) ======================= q-point: ( 0.43 0.07 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.208 ( -5.043 9.694 3.556) 11.492 1.474 ( 2.230 2.699 -5.629) 6.629 1.577 ( -3.643 6.450 -0.641) 7.435 1.878 ( -1.251 -1.104 -7.057) 7.251 3.213 ( 5.044 2.904 6.354) 8.617 3.837 ( 1.337 -1.502 -4.774) 5.180 4.119 ( 1.474 -1.553 7.159) 7.472 4.493 ( -0.489 -3.584 0.205) 3.623 4.641 ( -1.015 2.422 0.090) 2.627 4.685 ( -0.962 2.780 -0.469) 2.978 4.763 ( 0.185 0.602 -1.006) 1.187 4.926 ( 0.359 -0.554 -2.165) 2.263 ======================= Grid point 242 (39/96) ======================= q-point: ( 0.14 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.977 ( 4.590 7.951 3.118) 9.696 1.099 ( 6.804 11.785 4.276) 14.264 1.497 ( 4.712 8.162 -0.242) 9.428 1.870 ( 2.716 4.704 -3.545) 6.486 2.186 ( 10.831 18.760 8.961) 23.443 3.412 ( 2.660 4.606 -6.818) 8.648 4.355 ( -2.637 -4.567 3.532) 6.347 4.579 ( 1.490 2.581 -0.039) 2.980 4.583 ( 1.055 1.827 -0.429) 2.153 4.608 ( 0.556 0.963 1.483) 1.853 4.735 ( -0.302 -0.523 -1.343) 1.472 4.857 ( 0.974 1.688 -1.461) 2.435 ======================= Grid point 243 (40/96) ======================= q-point: ( 0.21 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.138 ( 5.326 7.228 -1.472) 9.098 1.291 ( 1.528 12.368 3.816) 13.033 1.633 ( 0.778 9.234 -0.241) 9.269 1.927 ( 1.622 -0.377 -4.990) 5.261 2.550 ( 12.899 13.368 7.980) 20.218 3.534 ( 5.533 5.482 -6.074) 9.878 4.252 ( -3.258 -5.733 4.934) 8.235 4.586 ( -1.004 -1.081 0.246) 1.496 4.630 ( 0.882 2.559 0.239) 2.717 4.637 ( 1.244 2.108 0.265) 2.462 4.743 ( 0.445 1.574 -1.363) 2.129 4.885 ( 1.138 0.298 -1.768) 2.124 ======================= Grid point 244 (41/96) ======================= q-point: ( 0.29 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.294 ( 4.351 6.913 -4.495) 9.324 1.402 ( -1.955 11.603 3.380) 12.243 1.705 ( -2.718 9.069 -0.612) 9.487 1.925 ( -0.273 -1.996 -6.033) 6.360 2.881 ( 12.199 9.535 7.226) 17.087 3.690 ( 6.359 4.023 -5.893) 9.558 4.139 ( -2.774 -5.843 6.748) 9.347 4.538 ( -1.609 -3.651 0.200) 3.994 4.662 ( 0.086 2.615 0.141) 2.621 4.682 ( 1.042 2.715 -0.189) 2.914 4.768 ( -0.255 2.608 -1.146) 2.860 4.902 ( 1.046 -0.660 -1.949) 2.309 ======================= Grid point 245 (42/96) ======================= q-point: ( 0.36 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.399 ( -2.110 9.346 -2.154) 9.820 1.478 ( 0.385 7.490 -0.992) 7.565 1.726 ( -4.383 7.764 -1.471) 9.037 1.890 ( -1.474 -1.270 -6.212) 6.509 3.140 ( 8.318 6.312 6.603) 12.355 3.789 ( 2.399 -1.194 -4.872) 5.560 4.075 ( 1.008 -2.364 7.219) 7.662 4.467 ( -1.431 -4.693 0.078) 4.907 4.683 ( -0.703 2.589 0.056) 2.683 4.724 ( 0.435 3.033 -0.456) 3.098 4.779 ( -1.047 1.862 -0.986) 2.353 4.909 ( 0.777 -0.956 -2.002) 2.351 ======================= Grid point 246 (43/96) ======================= q-point: ( 0.43 0.14 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.422 ( -5.888 10.199 0.261) 11.780 1.530 ( -2.010 3.482 -4.067) 5.719 1.722 ( -3.951 6.844 -2.276) 8.224 1.868 ( -0.376 0.651 -5.804) 5.853 3.253 ( -0.727 1.260 6.284) 6.451 3.784 ( 2.013 -3.487 -2.288) 4.631 4.091 ( 0.659 -1.141 5.237) 5.400 4.424 ( 1.892 -3.276 0.063) 3.784 4.690 ( -1.346 2.331 0.040) 2.692 4.748 ( -1.852 3.209 -0.601) 3.753 4.775 ( -0.380 0.659 -0.910) 1.186 4.911 ( 0.563 -0.975 -2.000) 2.295 ======================= Grid point 258 (44/96) ======================= q-point: ( 0.21 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.283 ( 4.276 7.407 -4.651) 9.735 1.494 ( 4.512 7.815 3.097) 9.541 1.790 ( 3.600 6.235 -0.058) 7.200 1.906 ( -0.911 -1.578 -6.349) 6.605 2.817 ( 7.818 13.542 7.510) 17.346 3.652 ( 3.551 6.151 -5.849) 9.201 4.133 ( -3.437 -5.952 6.520) 9.474 4.534 ( -2.102 -3.641 0.182) 4.208 4.674 ( 1.135 1.966 -0.175) 2.277 4.687 ( 1.495 2.589 0.047) 2.990 4.789 ( 1.451 2.513 -1.072) 3.094 4.885 ( -0.118 -0.204 -1.894) 1.908 ======================= Grid point 259 (45/96) ======================= q-point: ( 0.29 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.425 ( 3.991 6.551 -6.842) 10.279 1.602 ( -0.589 8.002 2.835) 8.510 1.844 ( -3.080 1.770 -4.988) 6.124 1.898 ( 1.069 2.033 -2.582) 3.456 3.075 ( 8.190 10.142 6.885) 14.742 3.746 ( 1.880 0.632 -4.795) 5.189 4.046 ( -0.515 -2.453 7.129) 7.557 4.452 ( -2.363 -4.925 0.031) 5.463 4.710 ( 0.492 2.247 -0.446) 2.343 4.733 ( 1.176 2.263 -0.096) 2.553 4.824 ( -0.788 2.696 -0.962) 2.969 4.882 ( 0.862 -1.284 -1.837) 2.401 ======================= Grid point 260 (46/96) ======================= q-point: ( 0.36 0.21 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.533 ( 0.803 5.442 -7.694) 9.458 1.646 ( -2.979 7.620 2.585) 8.580 1.792 ( -1.799 -1.498 -7.421) 7.781 1.928 ( -2.054 4.973 -0.966) 5.466 3.269 ( 3.352 6.871 6.227) 9.860 3.717 ( -0.021 -5.649 -0.163) 5.652 4.075 ( 1.820 2.220 3.656) 4.649 4.369 ( -0.587 -5.186 0.044) 5.219 4.730 ( -0.708 2.238 -0.359) 2.375 4.775 ( 1.191 1.974 -0.563) 2.374 4.814 ( -2.272 1.211 -0.984) 2.756 4.889 ( 0.864 -0.813 -1.583) 1.978 ======================= Grid point 273 (47/96) ======================= q-point: ( 0.29 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.551 ( 3.434 5.948 -8.129) 10.642 1.713 ( 1.941 3.362 2.192) 4.459 1.797 ( -2.119 -3.670 -7.890) 8.956 1.976 ( 1.771 3.067 -0.274) 3.553 3.275 ( 5.622 9.738 6.184) 12.833 3.687 ( -3.378 -5.850 0.651) 6.786 4.071 ( 2.442 4.230 3.089) 5.779 4.348 ( -3.169 -5.488 0.059) 6.338 4.753 ( 1.134 1.964 -1.083) 2.513 4.769 ( 0.858 1.486 0.074) 1.718 4.848 ( -0.919 -1.592 -1.706) 2.508 4.870 ( 0.742 1.285 -0.790) 1.681 ======================= Grid point 274 (48/96) ======================= q-point: ( 0.36 0.29 0.14) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.630 ( -2.431 4.210 -8.544) 9.830 1.736 ( 1.862 -3.226 -8.623) 9.393 1.749 ( -1.476 2.556 2.541) 3.895 2.005 ( -1.255 2.174 -0.362) 2.537 3.407 ( -3.770 6.530 5.417) 9.285 3.590 ( 3.659 -6.338 3.268) 8.015 4.143 ( -2.285 3.958 1.307) 4.753 4.268 ( 2.691 -4.661 0.264) 5.389 4.772 ( -1.150 1.992 -0.779) 2.428 4.798 ( -0.061 0.106 -0.395) 0.413 4.815 ( 0.365 -0.632 -1.600) 1.759 4.885 ( -0.098 0.170 -0.795) 0.819 ======================= Grid point 422 (49/96) ======================= q-point: ( 0.00 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.585 ( 0.000 0.000 12.078) 12.078 0.585 ( 0.000 0.000 12.078) 12.078 1.037 ( 0.000 0.000 -2.764) 2.764 1.037 ( 0.000 0.000 -2.764) 2.764 1.521 ( 0.000 0.000 33.926) 33.926 3.131 ( 0.000 -0.000 -16.112) 16.112 4.481 ( 0.000 0.000 -0.196) 0.196 4.481 ( -0.000 -0.000 -0.196) 0.196 4.481 ( 0.000 0.000 -0.199) 0.199 4.481 ( 0.000 0.000 -0.199) 0.199 4.752 ( -0.000 -0.000 1.233) 1.233 4.839 ( 0.000 0.000 -1.220) 1.220 ======================= Grid point 423 (50/96) ======================= q-point: ( 0.07 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.657 ( 5.721 3.303 10.853) 12.705 0.796 ( 13.281 7.668 10.180) 18.407 1.080 ( 3.596 2.076 -2.635) 4.917 1.202 ( 12.705 7.335 -2.796) 14.935 1.576 ( 6.068 3.503 31.046) 31.826 3.117 ( -0.954 -0.551 -15.815) 15.853 4.482 ( 0.388 0.224 -0.090) 0.457 4.492 ( 0.700 0.404 -0.332) 0.874 4.493 ( 0.940 0.543 0.388) 1.152 4.518 ( 2.732 1.577 -0.442) 3.185 4.754 ( 0.263 0.152 0.707) 0.769 4.809 ( -2.165 -1.250 -0.651) 2.583 ======================= Grid point 424 (51/96) ======================= q-point: ( 0.14 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.815 ( 7.534 4.350 8.904) 12.448 1.024 ( 6.133 3.541 13.401) 15.157 1.186 ( 5.297 3.058 -2.363) 6.557 1.520 ( 13.423 7.750 -3.553) 15.902 1.844 ( 16.868 9.739 20.742) 28.454 3.108 ( 0.799 0.461 -14.593) 14.623 4.498 ( 0.985 0.569 -0.071) 1.140 4.506 ( -0.102 -0.059 1.926) 1.930 4.512 ( 1.046 0.604 -0.415) 1.277 4.584 ( 2.667 1.540 0.120) 3.082 4.755 ( -1.623 -0.937 -0.532) 1.948 4.775 ( 0.994 0.574 -0.689) 1.339 ======================= Grid point 425 (52/96) ======================= q-point: ( 0.21 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.979 ( 6.518 3.763 7.555) 10.664 1.117 ( 2.977 1.719 12.217) 12.692 1.306 ( 4.954 2.860 -2.109) 6.097 1.756 ( 6.560 3.787 -5.865) 9.580 2.272 ( 19.345 11.169 13.956) 26.339 3.179 ( 5.638 3.255 -11.945) 13.604 4.484 ( -1.798 -1.038 3.587) 4.144 4.525 ( 1.266 0.731 -0.073) 1.463 4.540 ( 1.440 0.831 -0.479) 1.730 4.633 ( 1.728 0.998 1.054) 2.257 4.721 ( -1.283 -0.741 -1.167) 1.886 4.813 ( 1.860 1.074 -1.895) 2.864 ======================= Grid point 426 (53/96) ======================= q-point: ( 0.29 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.104 ( 4.366 2.521 6.797) 8.462 1.182 ( 2.882 1.664 6.917) 7.675 1.403 ( 3.414 1.971 -1.919) 4.384 1.822 ( -0.117 -0.068 -8.549) 8.550 2.693 ( 17.146 9.899 11.118) 22.707 3.355 ( 9.149 5.282 -8.911) 13.821 4.429 ( -2.845 -1.642 5.130) 6.091 4.554 ( 1.323 0.764 -0.011) 1.528 4.576 ( 1.639 0.946 -0.588) 1.981 4.666 ( 0.916 0.529 1.660) 1.968 4.691 ( -1.258 -0.726 -1.939) 2.423 4.849 ( 1.210 0.698 -2.682) 3.024 ======================= Grid point 427 (54/96) ======================= q-point: ( 0.36 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.179 ( 2.269 1.310 6.413) 6.928 1.252 ( 3.165 1.827 1.151) 3.831 1.461 ( 1.679 0.969 -1.759) 2.618 1.781 ( -2.933 -1.693 -10.245) 10.790 3.044 ( 13.514 7.802 9.648) 18.347 3.564 ( 8.576 4.951 -6.801) 12.013 4.360 ( -3.086 -1.782 6.532) 7.441 4.584 ( 1.208 0.697 0.053) 1.396 4.611 ( 1.400 0.809 -0.685) 1.756 4.631 ( -2.911 -1.681 -0.321) 3.376 4.704 ( 1.330 0.768 -0.388) 1.584 4.867 ( 0.393 0.227 -3.201) 3.233 ======================= Grid point 428 (55/96) ======================= q-point: ( 0.43 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.213 ( 0.841 0.486 6.239) 6.314 1.320 ( 2.618 1.512 -3.213) 4.412 1.484 ( 0.512 0.296 -1.621) 1.726 1.711 ( -2.793 -1.612 -10.946) 11.412 3.299 ( 8.415 4.859 8.677) 13.028 3.720 ( 4.604 2.658 -5.132) 7.389 4.302 ( -1.638 -0.946 7.043) 7.292 4.557 ( -3.186 -1.840 -0.446) 3.706 4.607 ( 0.750 0.433 0.075) 0.870 4.637 ( 0.801 0.463 -0.732) 1.180 4.731 ( 0.937 0.541 -0.569) 1.222 4.870 ( -0.000 -0.000 -3.464) 3.464 ======================= Grid point 429 (56/96) ======================= q-point: (-0.50 0.00 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.222 ( 0.000 0.000 6.187) 6.187 1.354 ( 0.000 0.000 -4.904) 4.904 1.489 ( 0.000 0.000 -1.561) 1.561 1.674 ( -0.000 -0.000 -11.053) 11.053 3.402 ( 0.000 0.000 8.229) 8.229 3.769 ( -0.000 -0.000 -3.943) 3.943 4.287 ( 0.000 0.000 6.614) 6.614 4.513 ( 0.000 0.000 -0.359) 0.359 4.616 ( 0.000 0.000 0.076) 0.076 4.646 ( 0.000 0.000 -0.744) 0.744 4.742 ( 0.000 0.000 -0.695) 0.695 4.870 ( -0.000 -0.000 -3.530) 3.530 ======================= Grid point 438 (57/96) ======================= q-point: ( 0.07 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.793 ( 5.375 9.310 9.079) 14.071 0.966 ( 4.227 7.322 12.135) 14.790 1.173 ( 3.736 6.470 -2.103) 7.762 1.424 ( 7.759 13.439 -3.261) 15.856 1.744 ( 8.111 14.049 23.842) 28.837 3.106 ( -0.069 -0.120 -15.059) 15.060 4.495 ( 0.318 0.550 0.385) 0.743 4.503 ( 0.715 1.239 -0.274) 1.457 4.509 ( 0.410 0.710 0.919) 1.232 4.564 ( 1.553 2.690 -0.212) 3.114 4.763 ( 0.549 0.951 -0.252) 1.127 4.771 ( -1.291 -2.235 -0.442) 2.619 ======================= Grid point 439 (58/96) ======================= q-point: ( 0.14 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.978 ( 3.492 10.983 7.657) 13.837 1.075 ( 4.065 2.105 12.592) 13.399 1.307 ( 2.968 8.391 -1.790) 9.079 1.669 ( 7.953 6.453 -4.891) 11.350 2.099 ( 14.841 15.243 16.065) 26.658 3.136 ( 2.792 2.712 -13.089) 13.656 4.486 ( -0.780 -1.695 2.658) 3.248 4.528 ( 0.426 2.110 0.017) 2.152 4.537 ( 1.099 1.821 -0.180) 2.135 4.614 ( 1.716 1.447 0.603) 2.324 4.730 ( -1.186 -1.462 -0.803) 2.047 4.793 ( 1.497 1.174 -1.500) 2.422 ======================= Grid point 440 (59/96) ======================= q-point: ( 0.21 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.101 ( 1.909 5.004 7.446) 9.172 1.179 ( 2.403 7.418 7.902) 11.101 1.427 ( 0.921 8.748 -1.772) 8.973 1.797 ( 3.068 0.403 -7.466) 8.082 2.509 ( 16.296 12.361 12.237) 23.835 3.263 ( 7.173 5.207 -10.140) 13.468 4.440 ( -1.687 -3.122 4.340) 5.605 4.557 ( 0.388 2.174 -0.127) 2.212 4.575 ( 0.775 2.468 -0.360) 2.612 4.652 ( 1.372 0.900 1.374) 2.140 4.701 ( -0.889 -1.091 -1.474) 2.038 4.830 ( 1.590 0.599 -2.389) 2.931 ======================= Grid point 441 (60/96) ======================= q-point: ( 0.29 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.173 ( 2.450 2.698 3.620) 5.136 1.278 ( -0.673 9.137 5.563) 10.719 1.507 ( -1.256 7.835 -1.845) 8.147 1.800 ( -1.130 -1.907 -9.449) 9.706 2.884 ( 14.678 9.319 10.318) 20.218 3.460 ( 8.993 5.287 -7.628) 12.924 4.372 ( -2.054 -3.938 5.861) 7.354 4.577 ( 0.202 0.132 -0.402) 0.469 4.610 ( 0.313 2.526 -0.401) 2.577 4.663 ( -0.905 -0.074 1.179) 1.488 4.695 ( 0.662 1.341 -1.427) 2.067 4.855 ( 1.025 -0.083 -2.977) 3.149 ======================= Grid point 442 (61/96) ======================= q-point: ( 0.36 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.238 ( 1.471 3.181 0.730) 3.580 1.340 ( -1.690 7.934 3.139) 8.698 1.542 ( -2.716 6.346 -2.191) 7.242 1.744 ( -2.816 -1.624 -10.163) 10.671 3.185 ( 11.123 6.489 9.125) 15.783 3.647 ( 7.048 3.365 -5.835) 9.749 4.301 ( -1.321 -3.836 6.907) 8.011 4.550 ( -0.996 -3.910 -0.648) 4.086 4.633 ( -0.595 2.713 0.047) 2.778 4.659 ( -0.813 2.629 -0.208) 2.760 4.723 ( 0.790 1.237 -0.740) 1.644 4.865 ( 0.458 -0.402 -3.313) 3.368 ======================= Grid point 443 (62/96) ======================= q-point: ( 0.43 0.07 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.277 ( -1.274 4.206 1.039) 4.516 1.381 ( -1.677 5.716 -0.332) 5.966 1.546 ( -2.922 4.717 -2.968) 6.292 1.691 ( -2.085 0.339 -9.679) 9.906 3.373 ( 4.694 2.696 8.326) 9.931 3.748 ( 2.512 0.204 -4.033) 4.756 4.265 ( 0.785 -2.585 6.622) 7.152 4.492 ( -0.346 -3.803 -0.521) 3.854 4.645 ( -1.113 2.562 0.213) 2.802 4.672 ( -1.088 2.709 -0.687) 2.999 4.742 ( 0.272 0.565 -0.760) 0.986 4.866 ( 0.233 -0.431 -3.431) 3.466 ======================= Grid point 453 (63/96) ======================= q-point: ( 0.14 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.107 ( 1.261 2.184 9.225) 9.563 1.211 ( 6.189 10.720 6.031) 13.769 1.481 ( 4.859 8.416 -1.413) 9.821 1.759 ( 1.498 2.594 -7.136) 7.739 2.421 ( 9.779 16.938 12.994) 23.481 3.225 ( 3.685 6.382 -10.707) 12.998 4.438 ( -1.810 -3.135 4.092) 5.464 4.571 ( 0.930 1.610 -0.679) 1.980 4.580 ( 1.510 2.616 0.179) 3.026 4.642 ( 0.855 1.480 1.169) 2.071 4.703 ( -0.646 -1.118 -1.221) 1.777 4.817 ( 0.659 1.141 -2.259) 2.616 ======================= Grid point 454 (64/96) ======================= q-point: ( 0.21 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.156 ( 2.003 2.330 3.950) 5.004 1.390 ( 1.705 11.118 5.268) 12.421 1.616 ( 0.134 9.487 -1.555) 9.615 1.779 ( 0.910 -1.956 -8.956) 9.212 2.755 ( 11.857 12.601 11.040) 20.524 3.381 ( 6.559 6.678 -8.279) 12.497 4.366 ( -2.222 -4.319 5.503) 7.340 4.579 ( -0.716 -0.863 -1.064) 1.545 4.631 ( 0.935 3.093 0.115) 3.233 4.663 ( 0.522 0.277 1.927) 2.016 4.706 ( 0.549 2.159 -1.953) 2.963 4.836 ( 0.875 0.022 -2.804) 2.937 ======================= Grid point 455 (65/96) ======================= q-point: ( 0.29 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.222 ( 2.710 2.662 -1.174) 3.976 1.490 ( -2.370 10.819 4.612) 11.997 1.663 ( -3.361 5.487 -4.815) 8.036 1.771 ( -0.698 0.822 -6.929) 7.012 3.064 ( 11.319 9.077 9.666) 17.434 3.562 ( 7.121 4.961 -6.314) 10.733 4.283 ( -1.850 -4.802 6.691) 8.440 4.534 ( -1.545 -3.781 -0.922) 4.187 4.667 ( 0.008 2.510 0.246) 2.522 4.684 ( 0.629 1.963 0.496) 2.120 4.739 ( -0.095 2.975 -1.408) 3.293 4.847 ( 0.862 -0.629 -3.078) 3.258 ======================= Grid point 456 (66/96) ======================= q-point: ( 0.36 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.293 ( 2.117 2.799 -4.826) 5.967 1.534 ( -4.222 10.089 3.778) 11.571 1.629 ( -3.155 0.323 -8.463) 9.038 1.772 ( -3.376 5.979 -3.554) 7.731 3.306 ( 7.673 5.880 8.596) 12.936 3.696 ( 4.141 1.619 -4.221) 6.131 4.223 ( 0.133 -3.738 6.739) 7.708 4.463 ( -1.331 -4.778 -0.737) 5.014 4.685 ( -0.522 2.124 0.140) 2.192 4.714 ( -0.003 2.781 -0.428) 2.813 4.756 ( -0.855 2.042 -1.044) 2.447 4.854 ( 0.586 -0.610 -3.104) 3.217 ======================= Grid point 457 (67/96) ======================= q-point: ( 0.43 0.14 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.330 ( -0.941 1.630 -6.086) 6.370 1.551 ( -5.709 9.889 4.139) 12.146 1.593 ( 0.057 -0.099 -10.020) 10.021 1.772 ( -4.163 7.210 -2.786) 8.779 3.409 ( -0.639 1.107 8.045) 8.146 3.734 ( 0.753 -1.304 -2.530) 2.945 4.211 ( 1.489 -2.579 5.858) 6.572 4.420 ( 1.924 -3.332 -0.580) 3.891 4.691 ( -1.099 1.904 0.086) 2.200 4.731 ( -1.635 2.832 -0.989) 3.416 4.756 ( -0.568 0.984 -0.688) 1.329 4.856 ( 0.297 -0.514 -3.068) 3.125 ======================= Grid point 469 (68/96) ======================= q-point: ( 0.21 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.211 ( 1.893 3.279 -1.016) 3.921 1.572 ( 3.975 6.885 3.749) 8.790 1.726 ( -1.494 -2.589 -10.100) 10.533 1.781 ( 3.615 6.261 -0.967) 7.294 3.006 ( 7.329 12.694 9.979) 17.732 3.522 ( 4.223 7.315 -6.524) 10.672 4.275 ( -2.699 -4.674 6.627) 8.547 4.528 ( -2.127 -3.683 -1.027) 4.375 4.676 ( 0.620 1.073 0.597) 1.376 4.690 ( 1.513 2.620 0.256) 3.036 4.761 ( 1.607 2.783 -1.473) 3.535 4.833 ( -0.150 -0.260 -2.982) 2.997 ======================= Grid point 470 (69/96) ======================= q-point: ( 0.29 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.284 ( 2.479 3.625 -5.296) 6.880 1.632 ( -1.709 -0.532 -8.074) 8.270 1.698 ( -1.329 4.691 0.677) 4.923 1.872 ( -0.068 6.107 -1.145) 6.214 3.246 ( 7.514 9.388 8.886) 14.952 3.652 ( 3.787 3.903 -4.363) 6.972 4.194 ( -1.299 -3.800 6.853) 7.943 4.446 ( -2.319 -4.922 -0.852) 5.508 4.699 ( 0.579 1.325 -0.373) 1.493 4.733 ( 0.678 2.353 0.074) 2.450 4.798 ( -0.855 2.490 -1.461) 3.011 4.833 ( 1.096 -0.622 -2.603) 2.892 ======================= Grid point 471 (70/96) ======================= q-point: ( 0.36 0.21 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.354 ( 0.966 3.341 -7.995) 8.719 1.579 ( -0.848 -2.981 -11.622) 12.028 1.722 ( -3.396 7.123 3.980) 8.838 1.907 ( -2.689 6.149 -1.244) 6.826 3.423 ( 3.024 6.042 7.823) 10.337 3.707 ( 0.976 -0.796 -1.042) 1.635 4.170 ( 1.015 -1.196 5.055) 5.293 4.366 ( -0.508 -4.980 -0.558) 5.037 4.717 ( -0.166 1.356 -0.792) 1.580 4.762 ( 0.414 1.855 -0.817) 2.069 4.788 ( -1.598 0.801 -1.063) 2.080 4.849 ( 0.378 0.403 -2.068) 2.141 ======================= Grid point 484 (71/96) ======================= q-point: ( 0.29 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.362 ( 2.403 4.162 -8.521) 9.783 1.571 ( -2.242 -3.883 -11.914) 12.730 1.788 ( 1.760 3.048 3.811) 5.187 1.963 ( 1.699 2.943 -1.047) 3.556 3.427 ( 4.907 8.499 7.719) 12.486 3.693 ( -0.199 -0.345 -0.372) 0.545 4.155 ( 0.166 0.287 4.677) 4.688 4.344 ( -2.968 -5.141 -0.523) 5.960 4.725 ( 0.753 1.304 -1.307) 1.994 4.767 ( 0.546 0.946 -0.705) 1.300 4.807 ( -0.421 -0.729 -1.557) 1.770 4.849 ( 0.811 1.405 -1.215) 2.027 ======================= Grid point 485 (72/96) ======================= q-point: ( 0.36 0.29 0.29) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.422 ( -1.879 3.255 -10.122) 10.797 1.514 ( 1.844 -3.194 -11.612) 12.184 1.817 ( -1.400 2.425 3.802) 4.721 1.991 ( -1.261 2.184 -1.110) 2.755 3.539 ( -2.961 5.129 6.481) 8.779 3.662 ( 1.904 -3.298 3.025) 4.864 4.181 ( -1.159 2.007 2.275) 3.247 4.273 ( 2.326 -4.028 0.068) 4.652 4.743 ( -0.685 1.187 -1.746) 2.219 4.771 ( 0.438 -0.759 -2.433) 2.586 4.800 ( -0.314 0.544 0.451) 0.774 4.863 ( -0.337 0.584 -1.205) 1.381 ======================= Grid point 633 (73/96) ======================= q-point: ( 0.00 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 96 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.804 ( 0.000 0.000 8.878) 8.878 0.804 ( 0.000 0.000 8.878) 8.878 0.953 ( 0.000 0.000 -5.520) 5.520 0.953 ( 0.000 0.000 -5.520) 5.520 2.178 ( 0.000 0.000 29.172) 29.172 2.722 ( 0.000 -0.000 -23.122) 23.122 4.477 ( -0.000 0.000 -0.089) 0.089 4.477 ( -0.000 0.000 -0.089) 0.089 4.477 ( 0.000 0.000 -0.089) 0.089 4.477 ( 0.000 0.000 -0.089) 0.089 4.780 ( 0.000 -0.000 1.445) 1.445 4.810 ( 0.000 0.000 -1.455) 1.455 ======================= Grid point 634 (74/96) ======================= q-point: ( 0.07 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.858 ( 4.481 2.587 8.367) 9.837 0.984 ( 12.966 7.486 7.862) 16.911 1.000 ( 3.895 2.249 -5.272) 6.930 1.120 ( 12.622 7.287 -5.202) 15.475 2.194 ( 1.921 1.109 27.917) 28.005 2.719 ( -0.001 -0.001 -22.447) 22.447 4.482 ( 0.468 0.270 0.070) 0.545 4.485 ( 0.583 0.337 -0.260) 0.721 4.501 ( 1.842 1.063 0.373) 2.159 4.509 ( 2.471 1.426 -0.432) 2.885 4.771 ( -0.647 -0.373 0.958) 1.215 4.792 ( -1.474 -0.851 -0.972) 1.960 ======================= Grid point 635 (75/96) ======================= q-point: ( 0.14 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.987 ( 6.336 3.658 7.370) 10.385 1.113 ( 5.661 3.269 -4.748) 8.080 1.260 ( 9.189 5.305 9.387) 14.167 1.420 ( 11.717 6.765 -6.180) 14.874 2.305 ( 8.658 4.999 22.609) 24.721 2.750 ( 3.493 2.017 -19.705) 20.114 4.498 ( 0.980 0.566 0.123) 1.139 4.503 ( 1.001 0.578 -0.341) 1.205 4.543 ( 1.451 0.838 1.638) 2.343 4.573 ( 2.687 1.551 -1.100) 3.292 4.752 ( -1.406 -0.812 0.249) 1.642 4.762 ( -0.238 -0.138 -0.564) 0.628 ======================= Grid point 636 (76/96) ======================= q-point: ( 0.21 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.129 ( 5.738 3.313 6.523) 9.298 1.241 ( 5.239 3.024 -4.260) 7.399 1.377 ( 1.694 0.978 12.046) 12.204 1.599 ( 3.762 2.172 -9.200) 10.175 2.583 ( 14.687 8.479 15.512) 22.983 2.900 ( 9.521 5.497 -14.674) 18.336 4.525 ( 1.303 0.752 0.138) 1.511 4.531 ( 1.362 0.786 -0.386) 1.619 4.555 ( -0.446 -0.257 3.225) 3.265 4.616 ( 0.874 0.504 -2.449) 2.649 4.730 ( -0.367 -0.212 1.775) 1.825 4.772 ( 0.872 0.504 -2.033) 2.268 ======================= Grid point 637 (77/96) ======================= q-point: ( 0.29 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.240 ( 3.929 2.269 5.977) 7.504 1.343 ( 3.616 2.088 -3.916) 5.725 1.375 ( -1.133 -0.654 11.082) 11.159 1.613 ( -1.734 -1.001 -11.251) 11.428 2.929 ( 14.935 8.623 11.321) 20.630 3.154 ( 11.994 6.924 -10.308) 17.264 4.529 ( -1.638 -0.946 4.457) 4.842 4.556 ( 1.411 0.814 0.231) 1.645 4.564 ( 1.524 0.880 -0.500) 1.829 4.615 ( -0.814 -0.470 -3.754) 3.870 4.730 ( 0.150 0.086 2.897) 2.902 4.790 ( 0.608 0.351 -2.895) 2.979 ======================= Grid point 638 (78/96) ======================= q-point: ( 0.36 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.307 ( 2.042 1.179 5.660) 6.132 1.345 ( -1.233 -0.712 7.767) 7.897 1.405 ( 1.842 1.063 -3.680) 4.250 1.551 ( -3.262 -1.883 -11.258) 11.872 3.242 ( 12.198 7.042 9.093) 16.765 3.416 ( 10.407 6.008 -7.646) 14.243 4.483 ( -2.323 -1.341 5.330) 5.967 4.580 ( -2.225 -1.285 -3.909) 4.678 4.587 ( 1.259 0.727 0.327) 1.490 4.597 ( 1.329 0.767 -0.611) 1.651 4.735 ( 0.282 0.163 2.654) 2.674 4.798 ( 0.085 0.049 -3.302) 3.304 ======================= Grid point 639 (79/96) ======================= q-point: ( 0.43 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.326 ( -0.448 -0.258 4.138) 4.170 1.338 ( 0.737 0.426 5.483) 5.548 1.432 ( 0.623 0.360 -3.512) 3.585 1.482 ( -2.460 -1.420 -10.179) 10.568 3.471 ( 7.503 4.332 7.612) 11.533 3.609 ( 6.179 3.567 -5.779) 9.181 4.431 ( -1.952 -1.127 5.344) 5.800 4.521 ( -2.612 -1.508 -3.142) 4.356 4.611 ( 0.760 0.439 0.366) 0.951 4.622 ( 0.778 0.449 -0.661) 1.115 4.743 ( 0.346 0.200 1.761) 1.806 4.797 ( -0.056 -0.033 -3.314) 3.315 ======================= Grid point 640 (80/96) ======================= q-point: (-0.50 0.00 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 203 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.322 ( 0.000 0.000 2.403) 2.403 1.345 ( 0.000 0.000 5.420) 5.420 1.438 ( 0.000 0.000 -3.444) 3.444 1.451 ( -0.000 -0.000 -9.423) 9.423 3.562 ( 0.000 0.000 6.884) 6.884 3.682 ( 0.000 0.000 -4.768) 4.768 4.407 ( 0.000 0.000 4.869) 4.869 4.485 ( -0.000 -0.000 -2.540) 2.540 4.620 ( 0.000 0.000 0.371) 0.371 4.631 ( 0.000 0.000 -0.671) 0.671 4.747 ( 0.000 0.000 1.349) 1.349 4.797 ( -0.000 -0.000 -3.233) 3.233 ======================= Grid point 649 (81/96) ======================= q-point: ( 0.07 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.968 ( 4.432 7.676 7.403) 11.549 1.090 ( 3.857 6.680 -3.886) 8.637 1.190 ( 6.107 10.577 7.682) 14.429 1.332 ( 7.081 12.264 -5.632) 15.240 2.258 ( 3.598 6.232 24.540) 25.573 2.732 ( 1.152 1.995 -20.732) 20.860 4.497 ( 0.587 1.016 0.011) 1.174 4.498 ( 0.653 1.130 -0.153) 1.314 4.533 ( 1.006 1.742 1.108) 2.296 4.553 ( 1.552 2.687 -0.833) 3.213 4.761 ( -0.261 -0.451 0.027) 0.522 4.764 ( -0.989 -1.713 -0.236) 1.992 ======================= Grid point 650 (82/96) ======================= q-point: ( 0.14 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.128 ( 3.417 9.758 6.414) 12.167 1.233 ( 3.629 7.999 -3.100) 9.315 1.342 ( 3.605 2.774 10.157) 11.129 1.537 ( 5.963 4.217 -7.937) 10.786 2.459 ( 10.019 10.415 18.017) 23.097 2.824 ( 5.621 5.764 -16.689) 18.529 4.518 ( 0.589 0.956 0.472) 1.218 4.525 ( 0.736 1.613 -0.281) 1.795 4.558 ( 0.469 0.785 2.139) 2.327 4.601 ( 1.489 1.179 -1.695) 2.545 4.734 ( -0.711 -0.934 1.028) 1.560 4.762 ( 0.535 0.171 -1.480) 1.583 ======================= Grid point 651 (83/96) ======================= q-point: ( 0.21 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.256 ( 0.762 7.877 6.519) 10.253 1.349 ( 1.537 5.083 1.422) 5.498 1.397 ( 0.357 4.257 5.536) 6.993 1.609 ( 1.200 -0.939 -10.384) 10.495 2.773 ( 13.526 10.381 12.853) 21.353 3.032 ( 10.153 7.682 -11.962) 17.469 4.518 ( -0.285 -1.798 2.847) 3.379 4.551 ( 0.548 0.984 -0.646) 1.299 4.578 ( 0.165 2.590 1.037) 2.795 4.623 ( 0.185 0.758 -2.776) 2.884 4.724 ( 0.230 -0.431 2.503) 2.550 4.778 ( 0.962 0.060 -2.582) 2.756 ======================= Grid point 652 (84/96) ======================= q-point: ( 0.29 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.297 ( -0.073 0.366 7.707) 7.716 1.399 ( -1.190 6.163 4.356) 7.640 1.450 ( -1.627 6.780 -2.279) 7.336 1.588 ( -2.062 -0.453 -10.147) 10.364 3.098 ( 13.076 8.321 9.996) 18.443 3.292 ( 11.024 6.833 -8.688) 15.610 4.481 ( -0.930 -3.259 4.532) 5.659 4.553 ( 0.091 -2.044 -2.208) 3.010 4.609 ( -0.130 3.258 0.447) 3.291 4.633 ( -0.769 2.563 -1.965) 3.320 4.728 ( 0.510 -0.184 2.812) 2.863 4.790 ( 0.624 -0.307 -3.136) 3.212 ======================= Grid point 653 (85/96) ======================= q-point: ( 0.36 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.290 ( 0.408 -1.501 5.423) 5.641 1.425 ( -2.078 4.222 -0.450) 4.727 1.462 ( -3.168 4.474 -1.850) 5.786 1.556 ( -3.786 4.768 -6.547) 8.940 3.369 ( 10.029 5.768 8.227) 14.196 3.521 ( 8.587 4.657 -6.531) 11.750 4.428 ( -0.994 -3.795 5.210) 6.522 4.514 ( -0.871 -4.015 -2.984) 5.078 4.635 ( -0.560 3.119 0.289) 3.182 4.647 ( -0.954 3.096 -0.832) 3.345 4.737 ( 0.491 0.205 2.022) 2.091 4.794 ( 0.256 -0.320 -3.264) 3.290 ======================= Grid point 654 (86/96) ======================= q-point: ( 0.43 0.07 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.287 ( 0.791 -1.343 2.832) 3.233 1.413 ( -0.180 -1.413 -7.909) 8.036 1.456 ( -4.719 8.092 4.265) 10.293 1.548 ( -4.322 7.219 -4.464) 9.525 3.536 ( 4.175 2.217 6.994) 8.442 3.658 ( 3.435 1.398 -4.883) 6.131 4.385 ( 0.380 -3.155 4.851) 5.799 4.462 ( -0.022 -3.649 -2.535) 4.443 4.649 ( -1.123 2.697 0.270) 2.934 4.658 ( -1.175 2.677 -0.591) 2.983 4.746 ( 0.189 0.194 1.339) 1.366 4.795 ( 0.109 -0.201 -3.154) 3.162 ======================= Grid point 664 (87/96) ======================= q-point: ( 0.14 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.326 ( 0.138 0.239 10.139) 10.143 1.331 ( 5.503 9.531 5.023) 12.097 1.440 ( 4.591 7.952 -1.759) 9.350 1.578 ( 0.169 0.292 -10.049) 10.054 2.702 ( 7.941 13.755 13.724) 20.991 2.980 ( 5.653 9.791 -12.787) 17.069 4.513 ( -0.868 -1.503 2.882) 3.364 4.553 ( 0.358 0.620 -1.183) 1.383 4.590 ( 1.380 2.390 0.994) 2.933 4.624 ( 0.893 1.547 -2.072) 2.736 4.719 ( -0.295 -0.510 2.143) 2.223 4.768 ( 0.215 0.372 -2.406) 2.444 ======================= Grid point 665 (88/96) ======================= q-point: ( 0.21 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.303 ( -1.019 -1.398 9.592) 9.747 1.478 ( 0.531 6.343 -0.294) 6.372 1.547 ( -1.204 5.517 -3.806) 6.810 1.604 ( 2.622 4.528 -5.391) 7.512 2.985 ( 10.331 11.071 10.804) 18.601 3.197 ( 8.417 8.888 -9.524) 15.510 4.469 ( -1.416 -3.138 4.284) 5.496 4.541 ( -0.847 -1.947 -2.641) 3.389 4.639 ( 0.741 3.298 0.718) 3.455 4.663 ( 0.577 2.830 -1.552) 3.279 4.716 ( 0.456 -0.196 2.619) 2.666 4.775 ( 0.576 -0.292 -2.963) 3.032 ======================= Grid point 666 (89/96) ======================= q-point: ( 0.29 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.276 ( -0.310 -1.436 6.417) 6.583 1.461 ( -2.157 -1.924 -10.013) 10.422 1.604 ( -2.472 9.522 3.772) 10.536 1.681 ( -1.887 8.486 -3.354) 9.318 3.259 ( 10.125 8.103 8.816) 15.681 3.424 ( 8.645 6.629 -7.146) 13.029 4.407 ( -1.419 -4.192 5.102) 6.754 4.492 ( -1.317 -3.850 -3.091) 5.109 4.675 ( -0.275 3.007 0.586) 3.076 4.689 ( -0.372 2.443 -0.544) 2.530 4.732 ( 0.666 1.143 1.415) 1.937 4.782 ( 0.699 -0.358 -2.995) 3.097 ======================= Grid point 667 (90/96) ======================= q-point: ( 0.36 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.266 ( 0.450 -0.867 2.868) 3.030 1.406 ( -1.274 -2.568 -9.713) 10.127 1.636 ( -4.753 9.797 4.182) 11.665 1.711 ( -4.357 8.991 -3.044) 10.444 3.474 ( 6.788 5.018 7.275) 11.144 3.602 ( 5.605 3.630 -5.157) 8.437 4.346 ( -0.411 -4.247 5.046) 6.608 4.428 ( -0.936 -4.576 -2.776) 5.434 4.691 ( -0.577 1.908 0.468) 2.047 4.704 ( -0.762 2.415 -0.669) 2.619 4.751 ( -0.120 1.452 0.891) 1.707 4.791 ( 0.298 0.188 -2.708) 2.731 ======================= Grid point 668 (91/96) ======================= q-point: ( 0.43 0.14 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.268 ( 0.401 -0.695 1.078) 1.344 1.379 ( 0.922 -1.598 -8.994) 9.182 1.642 ( -5.590 9.682 4.212) 11.947 1.717 ( -5.189 8.988 -2.973) 10.796 3.563 ( -0.412 0.714 6.425) 6.478 3.670 ( 0.007 -0.013 -3.930) 3.930 4.321 ( 1.798 -3.114 4.480) 5.744 4.392 ( 1.898 -3.288 -2.328) 4.454 4.695 ( -0.866 1.500 0.391) 1.776 4.710 ( -1.206 2.088 -0.961) 2.596 4.757 ( -0.674 1.167 0.937) 1.641 4.794 ( -0.210 0.363 -2.572) 2.606 ======================= Grid point 680 (92/96) ======================= q-point: ( 0.21 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.272 ( -0.800 -1.386 6.862) 7.046 1.469 ( -1.855 -3.213 -11.152) 11.753 1.680 ( 3.747 6.490 3.630) 8.327 1.746 ( 3.704 6.416 -2.461) 7.806 3.208 ( 6.522 11.297 9.114) 15.912 3.379 ( 5.419 9.385 -7.431) 13.140 4.398 ( -2.256 -3.907 5.115) 6.820 4.485 ( -2.048 -3.546 -3.190) 5.191 4.698 ( 0.258 0.447 1.055) 1.174 4.698 ( 1.390 2.408 0.522) 2.829 4.738 ( 1.407 2.437 -0.250) 2.825 4.770 ( 0.009 0.016 -2.772) 2.772 ======================= Grid point 681 (93/96) ======================= q-point: ( 0.29 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.255 ( 0.027 -0.575 3.095) 3.148 1.404 ( -1.806 -3.207 -10.418) 11.049 1.772 ( -0.577 6.833 3.737) 7.809 1.836 ( -0.406 6.505 -2.410) 6.949 3.420 ( 6.592 8.169 7.546) 12.928 3.554 ( 5.366 6.386 -5.402) 9.937 4.321 ( -1.685 -4.265 5.253) 6.973 4.407 ( -2.028 -4.601 -2.984) 5.847 4.703 ( 0.087 0.410 0.872) 0.968 4.728 ( 0.000 1.150 -0.937) 1.483 4.773 ( -0.012 1.946 -0.362) 1.979 4.788 ( 1.044 1.441 -1.541) 2.354 ======================= Grid point 682 (94/96) ======================= q-point: ( 0.36 0.21 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 686 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.258 ( 0.394 0.095 -0.226) 0.464 1.354 ( -0.378 -2.567 -8.809) 9.183 1.805 ( -3.138 6.813 3.724) 8.374 1.870 ( -2.933 6.474 -2.448) 7.517 3.570 ( 2.554 4.721 6.010) 8.057 3.665 ( 1.868 2.630 -3.148) 4.507 4.270 ( 0.322 -3.235 4.260) 5.359 4.338 ( -0.220 -4.364 -2.260) 4.920 4.709 ( 0.002 0.395 0.246) 0.465 4.731 ( -0.173 0.242 -1.725) 1.750 4.787 ( -0.616 1.838 0.720) 2.068 4.811 ( -0.258 1.641 -1.517) 2.250 ======================= Grid point 695 (95/96) ======================= q-point: ( 0.29 0.29 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 372 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.256 ( 0.365 0.631 -0.731) 1.032 1.344 ( -1.572 -2.723 -8.672) 9.224 1.867 ( 1.705 2.954 3.538) 4.914 1.929 ( 1.670 2.893 -2.299) 4.055 3.572 ( 3.863 6.691 5.849) 9.690 3.662 ( 2.283 3.955 -2.826) 5.370 4.250 ( -1.456 -2.522 4.161) 5.078 4.318 ( -2.431 -4.211 -2.226) 5.347 4.710 ( 0.205 0.355 0.074) 0.417 4.732 ( -0.180 -0.312 -1.873) 1.907 4.803 ( 0.654 1.134 0.671) 1.471 4.823 ( 0.822 1.423 -1.181) 2.024 ======================= Grid point 696 (96/96) ======================= q-point: ( 0.36 0.29 0.43) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 2.11e-04 2.11e-04 0.00e+00 0.00e+00 Number of triplets: 371 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.270 ( -0.541 0.938 -3.330) 3.502 1.309 ( 1.039 -1.800 -6.842) 7.151 1.895 ( -1.321 2.288 3.490) 4.377 1.956 ( -1.274 2.207 -2.309) 3.438 3.652 ( -1.781 3.085 4.029) 5.378 3.696 ( -0.125 0.217 -0.123) 0.279 4.230 ( 0.309 -0.535 2.219) 2.303 4.265 ( 1.557 -2.697 -0.974) 3.263 4.715 ( -0.195 0.338 -0.569) 0.690 4.726 ( 0.240 -0.415 -1.574) 1.645 4.814 ( -0.435 0.753 0.940) 1.281 4.837 ( -0.443 0.767 -1.231) 1.516 =================== 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/16464 10.0 647.256 647.256 698.730 -0.000 -0.000 -0.000 3/16464 20.0 351.687 351.687 362.470 -0.000 -0.000 -0.000 3/16464 30.0 202.561 202.561 213.149 -0.000 -0.000 -0.000 3/16464 40.0 120.255 120.255 129.636 -0.000 -0.000 -0.000 3/16464 50.0 79.911 79.911 87.604 -0.000 -0.000 -0.000 3/16464 60.0 58.548 58.548 64.872 -0.000 -0.000 -0.000 3/16464 70.0 45.948 45.948 51.268 -0.000 -0.000 -0.000 3/16464 80.0 37.804 37.804 42.387 -0.000 -0.000 -0.000 3/16464 90.0 32.154 32.154 36.180 0.000 -0.000 -0.000 3/16464 100.0 28.016 28.016 31.608 0.000 -0.000 -0.000 3/16464 110.0 24.856 24.856 28.099 0.000 -0.000 -0.000 3/16464 120.0 22.362 22.362 25.320 0.000 -0.000 -0.000 3/16464 130.0 20.342 20.342 23.062 0.000 -0.000 -0.000 3/16464 140.0 18.670 18.670 21.189 0.000 -0.000 -0.000 3/16464 150.0 17.262 17.262 19.608 0.000 -0.000 -0.000 3/16464 160.0 16.059 16.059 18.255 0.000 -0.000 -0.000 3/16464 170.0 15.019 15.019 17.082 0.000 -0.000 -0.000 3/16464 180.0 14.109 14.109 16.056 0.000 -0.000 -0.000 3/16464 190.0 13.307 13.307 15.150 0.000 -0.000 -0.000 3/16464 200.0 12.594 12.594 14.344 0.000 -0.000 -0.000 3/16464 210.0 11.955 11.955 13.621 0.000 -0.000 -0.000 3/16464 220.0 11.380 11.380 12.969 0.000 -0.000 -0.000 3/16464 230.0 10.858 10.858 12.378 0.000 -0.000 -0.000 3/16464 240.0 10.384 10.384 11.840 0.000 -0.000 -0.000 3/16464 250.0 9.950 9.950 11.348 0.000 -0.000 -0.000 3/16464 260.0 9.552 9.552 10.896 0.000 -0.000 -0.000 3/16464 270.0 9.185 9.185 10.479 0.000 -0.000 -0.000 3/16464 280.0 8.845 8.845 10.093 0.000 -0.000 -0.000 3/16464 290.0 8.531 8.531 9.735 0.000 -0.000 -0.000 3/16464 300.0 8.238 8.238 9.402 0.000 -0.000 -0.000 3/16464 310.0 7.965 7.965 9.091 0.000 -0.000 -0.000 3/16464 320.0 7.710 7.710 8.801 0.000 -0.000 -0.000 3/16464 330.0 7.470 7.470 8.529 0.000 -0.000 -0.000 3/16464 340.0 7.246 7.246 8.273 0.000 -0.000 -0.000 3/16464 350.0 7.034 7.034 8.032 0.000 -0.000 -0.000 3/16464 360.0 6.835 6.835 7.805 0.000 -0.000 -0.000 3/16464 370.0 6.647 6.647 7.591 0.000 -0.000 -0.000 3/16464 380.0 6.469 6.469 7.388 0.000 -0.000 -0.000 3/16464 390.0 6.301 6.301 7.196 0.000 -0.000 -0.000 3/16464 400.0 6.141 6.141 7.014 0.000 -0.000 -0.000 3/16464 410.0 5.989 5.989 6.841 0.000 -0.000 -0.000 3/16464 420.0 5.844 5.844 6.676 0.000 -0.000 -0.000 3/16464 430.0 5.707 5.707 6.519 0.000 -0.000 -0.000 3/16464 440.0 5.575 5.575 6.369 0.000 -0.000 -0.000 3/16464 450.0 5.450 5.450 6.226 0.000 -0.000 -0.000 3/16464 460.0 5.330 5.330 6.090 0.000 -0.000 -0.000 3/16464 470.0 5.216 5.216 5.959 0.000 -0.000 -0.000 3/16464 480.0 5.106 5.106 5.834 0.000 -0.000 -0.000 3/16464 490.0 5.001 5.001 5.714 0.000 -0.000 -0.000 3/16464 500.0 4.900 4.900 5.599 0.000 -0.000 -0.000 3/16464 510.0 4.803 4.803 5.488 0.000 -0.000 -0.000 3/16464 520.0 4.710 4.710 5.382 0.000 -0.000 -0.000 3/16464 530.0 4.621 4.621 5.280 0.000 -0.000 -0.000 3/16464 540.0 4.534 4.534 5.182 0.000 -0.000 -0.000 3/16464 550.0 4.451 4.451 5.087 0.000 -0.000 -0.000 3/16464 560.0 4.371 4.371 4.995 0.000 -0.000 -0.000 3/16464 570.0 4.294 4.294 4.907 0.000 -0.000 -0.000 3/16464 580.0 4.220 4.220 4.822 0.000 -0.000 -0.000 3/16464 590.0 4.148 4.148 4.740 0.000 -0.000 -0.000 3/16464 600.0 4.078 4.078 4.661 0.000 -0.000 -0.000 3/16464 610.0 4.011 4.011 4.584 0.000 -0.000 -0.000 3/16464 620.0 3.946 3.946 4.510 0.000 -0.000 -0.000 3/16464 630.0 3.883 3.883 4.438 0.000 -0.000 -0.000 3/16464 640.0 3.822 3.822 4.369 0.000 -0.000 -0.000 3/16464 650.0 3.763 3.763 4.301 0.000 -0.000 -0.000 3/16464 660.0 3.706 3.706 4.236 0.000 -0.000 -0.000 3/16464 670.0 3.650 3.650 4.172 0.000 -0.000 -0.000 3/16464 680.0 3.596 3.596 4.111 0.000 -0.000 -0.000 3/16464 690.0 3.544 3.544 4.051 0.000 -0.000 -0.000 3/16464 700.0 3.493 3.493 3.993 0.000 -0.000 -0.000 3/16464 710.0 3.444 3.444 3.937 0.000 -0.000 -0.000 3/16464 720.0 3.396 3.396 3.882 0.000 -0.000 -0.000 3/16464 730.0 3.349 3.349 3.828 0.000 -0.000 -0.000 3/16464 740.0 3.304 3.304 3.777 0.000 -0.000 -0.000 3/16464 750.0 3.260 3.260 3.726 0.000 -0.000 -0.000 3/16464 760.0 3.217 3.217 3.677 0.000 -0.000 -0.000 3/16464 770.0 3.175 3.175 3.629 0.000 -0.000 -0.000 3/16464 780.0 3.134 3.134 3.583 0.000 -0.000 -0.000 3/16464 790.0 3.094 3.094 3.537 0.000 -0.000 -0.000 3/16464 800.0 3.056 3.056 3.493 0.000 -0.000 -0.000 3/16464 810.0 3.018 3.018 3.450 0.000 -0.000 -0.000 3/16464 820.0 2.981 2.981 3.408 0.000 -0.000 -0.000 3/16464 830.0 2.945 2.945 3.366 0.000 -0.000 -0.000 3/16464 840.0 2.910 2.910 3.326 0.000 -0.000 -0.000 3/16464 850.0 2.875 2.875 3.287 0.000 -0.000 -0.000 3/16464 860.0 2.842 2.842 3.249 0.000 -0.000 -0.000 3/16464 870.0 2.809 2.809 3.211 0.000 -0.000 -0.000 3/16464 880.0 2.777 2.777 3.175 0.000 -0.000 -0.000 3/16464 890.0 2.746 2.746 3.139 0.000 -0.000 -0.000 3/16464 900.0 2.715 2.715 3.104 0.000 -0.000 -0.000 3/16464 910.0 2.686 2.686 3.070 0.000 -0.000 -0.000 3/16464 920.0 2.656 2.656 3.037 0.000 -0.000 -0.000 3/16464 930.0 2.628 2.628 3.004 0.000 -0.000 -0.000 3/16464 940.0 2.600 2.600 2.972 0.000 -0.000 -0.000 3/16464 950.0 2.572 2.572 2.941 0.000 -0.000 -0.000 3/16464 960.0 2.546 2.546 2.910 0.000 -0.000 -0.000 3/16464 970.0 2.519 2.519 2.880 0.000 -0.000 -0.000 3/16464 980.0 2.494 2.494 2.851 0.000 -0.000 -0.000 3/16464 990.0 2.468 2.468 2.822 0.000 -0.000 -0.000 3/16464 1000.0 2.444 2.444 2.794 0.000 -0.000 -0.000 3/16464 Thermal conductivity related properties were written into "kappa-m14147.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 18:33:06]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|