----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-09 00:53:34]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [3 3 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4/mmm (139) Number of symmetry operations in supercell: 288 ------------------------------ primitive cell ------------------------------ Lattice vectors: a -2.642727035000000 2.642727035000000 8.727389710000001 b 2.642727035000000 -2.642727035000000 8.727389710000001 c 2.642727035000000 2.642727035000000 -8.727389710000001 Atomic positions (fractional): *1 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 *2 K 0.14793473746599 0.14793473746599 0.00000000000000 39.098 3 K 0.85206526253401 0.85206526253401 0.00000000000000 39.098 4 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 *5 I 0.35340590890074 0.35340590890074 0.00000000000000 126.904 6 I 0.64659409109926 0.64659409109926 0.00000000000000 126.904 *7 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 -------------------------------- unit cell --------------------------------- Lattice vectors: a 5.285454070000000 0.000000000000000 0.000000000000000 b 0.000000000000000 5.285454070000000 0.000000000000000 c 0.000000000000000 0.000000000000000 17.454779420000001 Atomic positions (fractional): *1 K 0.50000000000000 0.00000000000000 0.00000000000000 39.098 > 1 *2 K 0.00000000000000 0.00000000000000 0.14793473746599 39.098 > 2 3 K 0.50000000000000 0.50000000000000 0.35206526253401 39.098 > 3 4 K 0.00000000000000 0.50000000000000 0.00000000000000 39.098 > 4 5 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 > 1 6 K 0.50000000000000 0.50000000000000 0.64793473746599 39.098 > 2 7 K 0.00000000000000 0.00000000000000 0.85206526253401 39.098 > 3 8 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 > 4 *9 I 0.00000000000000 0.00000000000000 0.35340590890074 126.904 > 5 10 I 0.50000000000000 0.50000000000000 0.14659409109926 126.904 > 6 11 I 0.50000000000000 0.50000000000000 0.85340590890074 126.904 > 5 12 I 0.00000000000000 0.00000000000000 0.64659409109926 126.904 > 6 *13 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 > 7 14 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 7 -------------------------------- super cell -------------------------------- Lattice vectors: a 15.856362210000000 0.000000000000000 0.000000000000000 b 0.000000000000000 15.856362210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 17.454779420000001 Atomic positions (fractional): *1 K 0.16666666666667 0.00000000000000 0.00000000000000 39.098 > 1 2 K 0.50000000000000 0.00000000000000 0.00000000000000 39.098 > 1 3 K 0.83333333333333 0.00000000000000 0.00000000000000 39.098 > 1 4 K 0.16666666666667 0.33333333333333 0.00000000000000 39.098 > 1 5 K 0.50000000000000 0.33333333333333 0.00000000000000 39.098 > 1 6 K 0.83333333333333 0.33333333333333 0.00000000000000 39.098 > 1 7 K 0.16666666666667 0.66666666666667 0.00000000000000 39.098 > 1 8 K 0.50000000000000 0.66666666666667 0.00000000000000 39.098 > 1 9 K 0.83333333333333 0.66666666666667 0.00000000000000 39.098 > 1 *10 K 0.00000000000000 0.00000000000000 0.14793473746599 39.098 > 2 11 K 0.33333333333333 0.00000000000000 0.14793473746599 39.098 > 2 12 K 0.66666666666667 0.00000000000000 0.14793473746599 39.098 > 2 13 K 0.00000000000000 0.33333333333333 0.14793473746599 39.098 > 2 14 K 0.33333333333333 0.33333333333333 0.14793473746599 39.098 > 2 15 K 0.66666666666667 0.33333333333333 0.14793473746599 39.098 > 2 16 K 0.00000000000000 0.66666666666667 0.14793473746599 39.098 > 2 17 K 0.33333333333333 0.66666666666667 0.14793473746599 39.098 > 2 18 K 0.66666666666667 0.66666666666667 0.14793473746599 39.098 > 2 19 K 0.16666666666667 0.16666666666667 0.35206526253401 39.098 > 3 20 K 0.50000000000000 0.16666666666667 0.35206526253401 39.098 > 3 21 K 0.83333333333333 0.16666666666667 0.35206526253401 39.098 > 3 22 K 0.16666666666667 0.50000000000000 0.35206526253401 39.098 > 3 23 K 0.50000000000000 0.50000000000000 0.35206526253401 39.098 > 3 24 K 0.83333333333333 0.50000000000000 0.35206526253401 39.098 > 3 25 K 0.16666666666667 0.83333333333333 0.35206526253401 39.098 > 3 26 K 0.50000000000000 0.83333333333333 0.35206526253401 39.098 > 3 27 K 0.83333333333333 0.83333333333333 0.35206526253401 39.098 > 3 28 K 0.00000000000000 0.16666666666667 0.00000000000000 39.098 > 4 29 K 0.33333333333333 0.16666666666667 0.00000000000000 39.098 > 4 30 K 0.66666666666667 0.16666666666667 0.00000000000000 39.098 > 4 31 K 0.00000000000000 0.50000000000000 0.00000000000000 39.098 > 4 32 K 0.33333333333333 0.50000000000000 0.00000000000000 39.098 > 4 33 K 0.66666666666667 0.50000000000000 0.00000000000000 39.098 > 4 34 K 0.00000000000000 0.83333333333333 0.00000000000000 39.098 > 4 35 K 0.33333333333333 0.83333333333333 0.00000000000000 39.098 > 4 36 K 0.66666666666667 0.83333333333333 0.00000000000000 39.098 > 4 37 K 0.00000000000000 0.16666666666667 0.50000000000000 39.098 > 1 38 K 0.33333333333333 0.16666666666667 0.50000000000000 39.098 > 1 39 K 0.66666666666667 0.16666666666667 0.50000000000000 39.098 > 1 40 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 > 1 41 K 0.33333333333333 0.50000000000000 0.50000000000000 39.098 > 1 42 K 0.66666666666667 0.50000000000000 0.50000000000000 39.098 > 1 43 K 0.00000000000000 0.83333333333333 0.50000000000000 39.098 > 1 44 K 0.33333333333333 0.83333333333333 0.50000000000000 39.098 > 1 45 K 0.66666666666667 0.83333333333333 0.50000000000000 39.098 > 1 46 K 0.16666666666667 0.16666666666667 0.64793473746599 39.098 > 2 47 K 0.50000000000000 0.16666666666667 0.64793473746599 39.098 > 2 48 K 0.83333333333333 0.16666666666667 0.64793473746599 39.098 > 2 49 K 0.16666666666667 0.50000000000000 0.64793473746599 39.098 > 2 50 K 0.50000000000000 0.50000000000000 0.64793473746599 39.098 > 2 51 K 0.83333333333333 0.50000000000000 0.64793473746599 39.098 > 2 52 K 0.16666666666667 0.83333333333333 0.64793473746599 39.098 > 2 53 K 0.50000000000000 0.83333333333333 0.64793473746599 39.098 > 2 54 K 0.83333333333333 0.83333333333333 0.64793473746599 39.098 > 2 55 K 0.00000000000000 0.00000000000000 0.85206526253401 39.098 > 3 56 K 0.33333333333333 0.00000000000000 0.85206526253401 39.098 > 3 57 K 0.66666666666667 0.00000000000000 0.85206526253401 39.098 > 3 58 K 0.00000000000000 0.33333333333333 0.85206526253401 39.098 > 3 59 K 0.33333333333333 0.33333333333333 0.85206526253401 39.098 > 3 60 K 0.66666666666667 0.33333333333333 0.85206526253401 39.098 > 3 61 K 0.00000000000000 0.66666666666667 0.85206526253401 39.098 > 3 62 K 0.33333333333333 0.66666666666667 0.85206526253401 39.098 > 3 63 K 0.66666666666667 0.66666666666667 0.85206526253401 39.098 > 3 64 K 0.16666666666667 0.00000000000000 0.50000000000000 39.098 > 4 65 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 > 4 66 K 0.83333333333333 0.00000000000000 0.50000000000000 39.098 > 4 67 K 0.16666666666667 0.33333333333333 0.50000000000000 39.098 > 4 68 K 0.50000000000000 0.33333333333333 0.50000000000000 39.098 > 4 69 K 0.83333333333333 0.33333333333333 0.50000000000000 39.098 > 4 70 K 0.16666666666667 0.66666666666667 0.50000000000000 39.098 > 4 71 K 0.50000000000000 0.66666666666667 0.50000000000000 39.098 > 4 72 K 0.83333333333333 0.66666666666667 0.50000000000000 39.098 > 4 *73 I 0.00000000000000 0.00000000000000 0.35340590890074 126.904 > 5 74 I 0.33333333333333 0.00000000000000 0.35340590890074 126.904 > 5 75 I 0.66666666666667 0.00000000000000 0.35340590890074 126.904 > 5 76 I 0.00000000000000 0.33333333333333 0.35340590890074 126.904 > 5 77 I 0.33333333333333 0.33333333333333 0.35340590890074 126.904 > 5 78 I 0.66666666666667 0.33333333333333 0.35340590890074 126.904 > 5 79 I 0.00000000000000 0.66666666666667 0.35340590890074 126.904 > 5 80 I 0.33333333333333 0.66666666666667 0.35340590890074 126.904 > 5 81 I 0.66666666666667 0.66666666666667 0.35340590890074 126.904 > 5 82 I 0.16666666666667 0.16666666666667 0.14659409109926 126.904 > 6 83 I 0.50000000000000 0.16666666666667 0.14659409109926 126.904 > 6 84 I 0.83333333333333 0.16666666666667 0.14659409109926 126.904 > 6 85 I 0.16666666666667 0.50000000000000 0.14659409109926 126.904 > 6 86 I 0.50000000000000 0.50000000000000 0.14659409109926 126.904 > 6 87 I 0.83333333333333 0.50000000000000 0.14659409109926 126.904 > 6 88 I 0.16666666666667 0.83333333333333 0.14659409109926 126.904 > 6 89 I 0.50000000000000 0.83333333333333 0.14659409109926 126.904 > 6 90 I 0.83333333333333 0.83333333333333 0.14659409109926 126.904 > 6 91 I 0.16666666666667 0.16666666666667 0.85340590890074 126.904 > 5 92 I 0.50000000000000 0.16666666666667 0.85340590890074 126.904 > 5 93 I 0.83333333333333 0.16666666666667 0.85340590890074 126.904 > 5 94 I 0.16666666666667 0.50000000000000 0.85340590890074 126.904 > 5 95 I 0.50000000000000 0.50000000000000 0.85340590890074 126.904 > 5 96 I 0.83333333333333 0.50000000000000 0.85340590890074 126.904 > 5 97 I 0.16666666666667 0.83333333333333 0.85340590890074 126.904 > 5 98 I 0.50000000000000 0.83333333333333 0.85340590890074 126.904 > 5 99 I 0.83333333333333 0.83333333333333 0.85340590890074 126.904 > 5 100 I 0.00000000000000 0.00000000000000 0.64659409109926 126.904 > 6 101 I 0.33333333333333 0.00000000000000 0.64659409109926 126.904 > 6 102 I 0.66666666666667 0.00000000000000 0.64659409109926 126.904 > 6 103 I 0.00000000000000 0.33333333333333 0.64659409109926 126.904 > 6 104 I 0.33333333333333 0.33333333333333 0.64659409109926 126.904 > 6 105 I 0.66666666666667 0.33333333333333 0.64659409109926 126.904 > 6 106 I 0.00000000000000 0.66666666666667 0.64659409109926 126.904 > 6 107 I 0.33333333333333 0.66666666666667 0.64659409109926 126.904 > 6 108 I 0.66666666666667 0.66666666666667 0.64659409109926 126.904 > 6 *109 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 > 7 110 O 0.33333333333333 0.00000000000000 0.00000000000000 15.999 > 7 111 O 0.66666666666667 0.00000000000000 0.00000000000000 15.999 > 7 112 O 0.00000000000000 0.33333333333333 0.00000000000000 15.999 > 7 113 O 0.33333333333333 0.33333333333333 0.00000000000000 15.999 > 7 114 O 0.66666666666667 0.33333333333333 0.00000000000000 15.999 > 7 115 O 0.00000000000000 0.66666666666667 0.00000000000000 15.999 > 7 116 O 0.33333333333333 0.66666666666667 0.00000000000000 15.999 > 7 117 O 0.66666666666667 0.66666666666667 0.00000000000000 15.999 > 7 118 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 7 119 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 7 120 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 7 121 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 7 122 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 7 123 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 7 124 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 7 125 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 7 126 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 7 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.0759119 0.0000000 0.0000000 0.0000000 3.0759119 0.0000000 0.0000000 0.0000000 3.0891469 -------------------------- Born effective charges -------------------------- 1 K 1.1710065 0.0000000 0.0000000 0.0000000 1.1424984 0.0000000 0.0000000 0.0000000 0.9864279 2 K 1.1290615 0.0000000 0.0000000 0.0000000 1.1290615 0.0000000 0.0000000 0.0000000 1.2349842 3 K 1.1290615 0.0000000 0.0000000 0.0000000 1.1290615 0.0000000 0.0000000 0.0000000 1.2349842 4 K 1.1424984 0.0000000 0.0000000 0.0000000 1.1710065 0.0000000 0.0000000 0.0000000 0.9864279 5 I -1.1615912 0.0000000 0.0000000 0.0000000 -1.1615912 0.0000000 0.0000000 0.0000000 -1.1299427 6 I -1.1615912 0.0000000 0.0000000 0.0000000 -1.1615912 0.0000000 0.0000000 0.0000000 -1.1299427 7 O -2.2484455 0.0000000 0.0000000 0.0000000 -2.2484455 0.0000000 0.0000000 0.0000000 -2.1829389 ---------------------------------------------------------------------------- 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: 378/378 Permutation basis: 5061/5061 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 352 Number of blocks in projector: 358 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 6 --- Eigsh_solver_block: 1 / 6 --- Block_size: 144 Use standard eigh solver. --- Eigsh_solver_block: 2 / 6 --- Block_size: 92 Use standard eigh solver. --- Eigsh_solver_block: 3 / 6 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 4 / 6 --- Block_size: 34 Use standard eigh solver. --- Eigsh_solver_block: 5 / 6 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 6 / 6 --- Block_size: 1 Use standard eigh solver. Tree of FC basis block matrices: - (352, 343), data: False |-- (1, 1), data: True |-- (1, 1), data: True |-- (34, 34), data: True |-- (80, 80), data: True |-- (92, 88), data: True |-- (144, 139), data: True ----- Solver_atoms: 1 -- 126 / 126 Time (Solver_compr_matrix_reshape): 0.005 Solver_block: 80 / 80 - Time: 0.264 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.271 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 378/378 Permutation basis: 5061/5061 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 352 Number of blocks in projector: 358 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 6 --- Eigsh_solver_block: 1 / 6 --- Block_size: 144 Use standard eigh solver. --- Eigsh_solver_block: 2 / 6 --- Block_size: 92 Use standard eigh solver. --- Eigsh_solver_block: 3 / 6 --- Block_size: 80 Use standard eigh solver. --- Eigsh_solver_block: 4 / 6 --- Block_size: 34 Use standard eigh solver. --- Eigsh_solver_block: 5 / 6 --- Block_size: 1 Use standard eigh solver. --- Eigsh_solver_block: 6 / 6 --- Block_size: 1 Use standard eigh solver. Tree of FC basis block matrices: - (352, 343), data: False |-- (1, 1), data: True |-- (1, 1), data: True |-- (34, 34), data: True |-- (80, 80), data: True |-- (92, 88), data: True |-- (144, 139), data: True Max drift after symmetrization by symfc projector: -0.00000000 (zz) -0.00000000 (zz) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-09 00:53:38]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 00:53:39]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [3 3 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4/mmm (139) ------------------------------ primitive cell ------------------------------ Lattice vectors: a -2.642727035000000 2.642727035000000 8.727389710000001 b 2.642727035000000 -2.642727035000000 8.727389710000001 c 2.642727035000000 2.642727035000000 -8.727389710000001 Atomic positions (fractional): 1 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 2 K 0.14793473746599 0.14793473746599 0.00000000000000 39.098 3 K 0.85206526253401 0.85206526253401 0.00000000000000 39.098 4 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 5 I 0.35340590890074 0.35340590890074 0.00000000000000 126.904 6 I 0.64659409109926 0.64659409109926 0.00000000000000 126.904 7 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.856362210000000 0.000000000000000 0.000000000000000 b 0.000000000000000 15.856362210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 17.454779420000001 Atomic positions (fractional): 1 K 0.16666666666667 0.00000000000000 0.00000000000000 39.098 > 1 2 K 0.50000000000000 0.00000000000000 0.00000000000000 39.098 > 1 3 K 0.83333333333333 0.00000000000000 0.00000000000000 39.098 > 1 4 K 0.16666666666667 0.33333333333333 0.00000000000000 39.098 > 1 5 K 0.50000000000000 0.33333333333333 0.00000000000000 39.098 > 1 6 K 0.83333333333333 0.33333333333333 0.00000000000000 39.098 > 1 7 K 0.16666666666667 0.66666666666667 0.00000000000000 39.098 > 1 8 K 0.50000000000000 0.66666666666667 0.00000000000000 39.098 > 1 9 K 0.83333333333333 0.66666666666667 0.00000000000000 39.098 > 1 10 K 0.00000000000000 0.00000000000000 0.14793473746599 39.098 > 10 11 K 0.33333333333333 0.00000000000000 0.14793473746599 39.098 > 10 12 K 0.66666666666667 0.00000000000000 0.14793473746599 39.098 > 10 13 K 0.00000000000000 0.33333333333333 0.14793473746599 39.098 > 10 14 K 0.33333333333333 0.33333333333333 0.14793473746599 39.098 > 10 15 K 0.66666666666667 0.33333333333333 0.14793473746599 39.098 > 10 16 K 0.00000000000000 0.66666666666667 0.14793473746599 39.098 > 10 17 K 0.33333333333333 0.66666666666667 0.14793473746599 39.098 > 10 18 K 0.66666666666667 0.66666666666667 0.14793473746599 39.098 > 10 19 K 0.16666666666667 0.16666666666667 0.35206526253401 39.098 > 19 20 K 0.50000000000000 0.16666666666667 0.35206526253401 39.098 > 19 21 K 0.83333333333333 0.16666666666667 0.35206526253401 39.098 > 19 22 K 0.16666666666667 0.50000000000000 0.35206526253401 39.098 > 19 23 K 0.50000000000000 0.50000000000000 0.35206526253401 39.098 > 19 24 K 0.83333333333333 0.50000000000000 0.35206526253401 39.098 > 19 25 K 0.16666666666667 0.83333333333333 0.35206526253401 39.098 > 19 26 K 0.50000000000000 0.83333333333333 0.35206526253401 39.098 > 19 27 K 0.83333333333333 0.83333333333333 0.35206526253401 39.098 > 19 28 K 0.00000000000000 0.16666666666667 0.00000000000000 39.098 > 28 29 K 0.33333333333333 0.16666666666667 0.00000000000000 39.098 > 28 30 K 0.66666666666667 0.16666666666667 0.00000000000000 39.098 > 28 31 K 0.00000000000000 0.50000000000000 0.00000000000000 39.098 > 28 32 K 0.33333333333333 0.50000000000000 0.00000000000000 39.098 > 28 33 K 0.66666666666667 0.50000000000000 0.00000000000000 39.098 > 28 34 K 0.00000000000000 0.83333333333333 0.00000000000000 39.098 > 28 35 K 0.33333333333333 0.83333333333333 0.00000000000000 39.098 > 28 36 K 0.66666666666667 0.83333333333333 0.00000000000000 39.098 > 28 37 K 0.00000000000000 0.16666666666667 0.50000000000000 39.098 > 1 38 K 0.33333333333333 0.16666666666667 0.50000000000000 39.098 > 1 39 K 0.66666666666667 0.16666666666667 0.50000000000000 39.098 > 1 40 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 > 1 41 K 0.33333333333333 0.50000000000000 0.50000000000000 39.098 > 1 42 K 0.66666666666667 0.50000000000000 0.50000000000000 39.098 > 1 43 K 0.00000000000000 0.83333333333333 0.50000000000000 39.098 > 1 44 K 0.33333333333333 0.83333333333333 0.50000000000000 39.098 > 1 45 K 0.66666666666667 0.83333333333333 0.50000000000000 39.098 > 1 46 K 0.16666666666667 0.16666666666667 0.64793473746599 39.098 > 10 47 K 0.50000000000000 0.16666666666667 0.64793473746599 39.098 > 10 48 K 0.83333333333333 0.16666666666667 0.64793473746599 39.098 > 10 49 K 0.16666666666667 0.50000000000000 0.64793473746599 39.098 > 10 50 K 0.50000000000000 0.50000000000000 0.64793473746599 39.098 > 10 51 K 0.83333333333333 0.50000000000000 0.64793473746599 39.098 > 10 52 K 0.16666666666667 0.83333333333333 0.64793473746599 39.098 > 10 53 K 0.50000000000000 0.83333333333333 0.64793473746599 39.098 > 10 54 K 0.83333333333333 0.83333333333333 0.64793473746599 39.098 > 10 55 K 0.00000000000000 0.00000000000000 0.85206526253401 39.098 > 19 56 K 0.33333333333333 0.00000000000000 0.85206526253401 39.098 > 19 57 K 0.66666666666667 0.00000000000000 0.85206526253401 39.098 > 19 58 K 0.00000000000000 0.33333333333333 0.85206526253401 39.098 > 19 59 K 0.33333333333333 0.33333333333333 0.85206526253401 39.098 > 19 60 K 0.66666666666667 0.33333333333333 0.85206526253401 39.098 > 19 61 K 0.00000000000000 0.66666666666667 0.85206526253401 39.098 > 19 62 K 0.33333333333333 0.66666666666667 0.85206526253401 39.098 > 19 63 K 0.66666666666667 0.66666666666667 0.85206526253401 39.098 > 19 64 K 0.16666666666667 0.00000000000000 0.50000000000000 39.098 > 28 65 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 > 28 66 K 0.83333333333333 0.00000000000000 0.50000000000000 39.098 > 28 67 K 0.16666666666667 0.33333333333333 0.50000000000000 39.098 > 28 68 K 0.50000000000000 0.33333333333333 0.50000000000000 39.098 > 28 69 K 0.83333333333333 0.33333333333333 0.50000000000000 39.098 > 28 70 K 0.16666666666667 0.66666666666667 0.50000000000000 39.098 > 28 71 K 0.50000000000000 0.66666666666667 0.50000000000000 39.098 > 28 72 K 0.83333333333333 0.66666666666667 0.50000000000000 39.098 > 28 73 I 0.00000000000000 0.00000000000000 0.35340590890074 126.904 > 73 74 I 0.33333333333333 0.00000000000000 0.35340590890074 126.904 > 73 75 I 0.66666666666667 0.00000000000000 0.35340590890074 126.904 > 73 76 I 0.00000000000000 0.33333333333333 0.35340590890074 126.904 > 73 77 I 0.33333333333333 0.33333333333333 0.35340590890074 126.904 > 73 78 I 0.66666666666667 0.33333333333333 0.35340590890074 126.904 > 73 79 I 0.00000000000000 0.66666666666667 0.35340590890074 126.904 > 73 80 I 0.33333333333333 0.66666666666667 0.35340590890074 126.904 > 73 81 I 0.66666666666667 0.66666666666667 0.35340590890074 126.904 > 73 82 I 0.16666666666667 0.16666666666667 0.14659409109926 126.904 > 82 83 I 0.50000000000000 0.16666666666667 0.14659409109926 126.904 > 82 84 I 0.83333333333333 0.16666666666667 0.14659409109926 126.904 > 82 85 I 0.16666666666667 0.50000000000000 0.14659409109926 126.904 > 82 86 I 0.50000000000000 0.50000000000000 0.14659409109926 126.904 > 82 87 I 0.83333333333333 0.50000000000000 0.14659409109926 126.904 > 82 88 I 0.16666666666667 0.83333333333333 0.14659409109926 126.904 > 82 89 I 0.50000000000000 0.83333333333333 0.14659409109926 126.904 > 82 90 I 0.83333333333333 0.83333333333333 0.14659409109926 126.904 > 82 91 I 0.16666666666667 0.16666666666667 0.85340590890074 126.904 > 73 92 I 0.50000000000000 0.16666666666667 0.85340590890074 126.904 > 73 93 I 0.83333333333333 0.16666666666667 0.85340590890074 126.904 > 73 94 I 0.16666666666667 0.50000000000000 0.85340590890074 126.904 > 73 95 I 0.50000000000000 0.50000000000000 0.85340590890074 126.904 > 73 96 I 0.83333333333333 0.50000000000000 0.85340590890074 126.904 > 73 97 I 0.16666666666667 0.83333333333333 0.85340590890074 126.904 > 73 98 I 0.50000000000000 0.83333333333333 0.85340590890074 126.904 > 73 99 I 0.83333333333333 0.83333333333333 0.85340590890074 126.904 > 73 100 I 0.00000000000000 0.00000000000000 0.64659409109926 126.904 > 82 101 I 0.33333333333333 0.00000000000000 0.64659409109926 126.904 > 82 102 I 0.66666666666667 0.00000000000000 0.64659409109926 126.904 > 82 103 I 0.00000000000000 0.33333333333333 0.64659409109926 126.904 > 82 104 I 0.33333333333333 0.33333333333333 0.64659409109926 126.904 > 82 105 I 0.66666666666667 0.33333333333333 0.64659409109926 126.904 > 82 106 I 0.00000000000000 0.66666666666667 0.64659409109926 126.904 > 82 107 I 0.33333333333333 0.66666666666667 0.64659409109926 126.904 > 82 108 I 0.66666666666667 0.66666666666667 0.64659409109926 126.904 > 82 109 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 > 109 110 O 0.33333333333333 0.00000000000000 0.00000000000000 15.999 > 109 111 O 0.66666666666667 0.00000000000000 0.00000000000000 15.999 > 109 112 O 0.00000000000000 0.33333333333333 0.00000000000000 15.999 > 109 113 O 0.33333333333333 0.33333333333333 0.00000000000000 15.999 > 109 114 O 0.66666666666667 0.33333333333333 0.00000000000000 15.999 > 109 115 O 0.00000000000000 0.66666666666667 0.00000000000000 15.999 > 109 116 O 0.33333333333333 0.66666666666667 0.00000000000000 15.999 > 109 117 O 0.66666666666667 0.66666666666667 0.00000000000000 15.999 > 109 118 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 109 119 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 109 120 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 109 121 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 109 122 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 109 123 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 109 124 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 109 125 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 109 126 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.0759119 0.0000000 0.0000000 0.0000000 3.0759119 0.0000000 0.0000000 0.0000000 3.0891469 -------------------------- Born effective charges -------------------------- 1 K 1.1710065 0.0000000 0.0000000 0.0000000 1.1424984 0.0000000 0.0000000 0.0000000 0.9864279 2 K 1.1290615 0.0000000 0.0000000 0.0000000 1.1290615 0.0000000 0.0000000 0.0000000 1.2349842 3 K 1.1290615 0.0000000 0.0000000 0.0000000 1.1290615 0.0000000 0.0000000 0.0000000 1.2349842 4 K 1.1424984 0.0000000 0.0000000 0.0000000 1.1710065 0.0000000 0.0000000 0.0000000 0.9864279 5 I -1.1615912 0.0000000 0.0000000 0.0000000 -1.1615912 0.0000000 0.0000000 0.0000000 -1.1299427 6 I -1.1615912 0.0000000 0.0000000 0.0000000 -1.1615912 0.0000000 0.0000000 0.0000000 -1.1299427 7 O -2.2484455 0.0000000 0.0000000 0.0000000 -2.2484455 0.0000000 0.0000000 0.0000000 -2.1829389 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 10, 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[ 73, 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[ 109, 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.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-09 00:53:46]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 00:53:46]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [3 3 1] Primitive matrix: [-0.5 0.5 0.5] [ 0.5 -0.5 0.5] [ 0.5 0.5 -0.5] Spacegroup: I4/mmm (139) ------------------------------ primitive cell ------------------------------ Lattice vectors: a -2.642727035000000 2.642727035000000 8.727389710000001 b 2.642727035000000 -2.642727035000000 8.727389710000001 c 2.642727035000000 2.642727035000000 -8.727389710000001 Atomic positions (fractional): 1 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 2 K 0.14793473746599 0.14793473746599 0.00000000000000 39.098 3 K 0.85206526253401 0.85206526253401 0.00000000000000 39.098 4 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 5 I 0.35340590890074 0.35340590890074 0.00000000000000 126.904 6 I 0.64659409109926 0.64659409109926 0.00000000000000 126.904 7 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 -------------------------------- supercell --------------------------------- Lattice vectors: a 15.856362210000000 0.000000000000000 0.000000000000000 b 0.000000000000000 15.856362210000000 0.000000000000000 c 0.000000000000000 0.000000000000000 17.454779420000001 Atomic positions (fractional): 1 K 0.16666666666667 0.00000000000000 0.00000000000000 39.098 > 1 2 K 0.50000000000000 0.00000000000000 0.00000000000000 39.098 > 1 3 K 0.83333333333333 0.00000000000000 0.00000000000000 39.098 > 1 4 K 0.16666666666667 0.33333333333333 0.00000000000000 39.098 > 1 5 K 0.50000000000000 0.33333333333333 0.00000000000000 39.098 > 1 6 K 0.83333333333333 0.33333333333333 0.00000000000000 39.098 > 1 7 K 0.16666666666667 0.66666666666667 0.00000000000000 39.098 > 1 8 K 0.50000000000000 0.66666666666667 0.00000000000000 39.098 > 1 9 K 0.83333333333333 0.66666666666667 0.00000000000000 39.098 > 1 10 K 0.00000000000000 0.00000000000000 0.14793473746599 39.098 > 10 11 K 0.33333333333333 0.00000000000000 0.14793473746599 39.098 > 10 12 K 0.66666666666667 0.00000000000000 0.14793473746599 39.098 > 10 13 K 0.00000000000000 0.33333333333333 0.14793473746599 39.098 > 10 14 K 0.33333333333333 0.33333333333333 0.14793473746599 39.098 > 10 15 K 0.66666666666667 0.33333333333333 0.14793473746599 39.098 > 10 16 K 0.00000000000000 0.66666666666667 0.14793473746599 39.098 > 10 17 K 0.33333333333333 0.66666666666667 0.14793473746599 39.098 > 10 18 K 0.66666666666667 0.66666666666667 0.14793473746599 39.098 > 10 19 K 0.16666666666667 0.16666666666667 0.35206526253401 39.098 > 19 20 K 0.50000000000000 0.16666666666667 0.35206526253401 39.098 > 19 21 K 0.83333333333333 0.16666666666667 0.35206526253401 39.098 > 19 22 K 0.16666666666667 0.50000000000000 0.35206526253401 39.098 > 19 23 K 0.50000000000000 0.50000000000000 0.35206526253401 39.098 > 19 24 K 0.83333333333333 0.50000000000000 0.35206526253401 39.098 > 19 25 K 0.16666666666667 0.83333333333333 0.35206526253401 39.098 > 19 26 K 0.50000000000000 0.83333333333333 0.35206526253401 39.098 > 19 27 K 0.83333333333333 0.83333333333333 0.35206526253401 39.098 > 19 28 K 0.00000000000000 0.16666666666667 0.00000000000000 39.098 > 28 29 K 0.33333333333333 0.16666666666667 0.00000000000000 39.098 > 28 30 K 0.66666666666667 0.16666666666667 0.00000000000000 39.098 > 28 31 K 0.00000000000000 0.50000000000000 0.00000000000000 39.098 > 28 32 K 0.33333333333333 0.50000000000000 0.00000000000000 39.098 > 28 33 K 0.66666666666667 0.50000000000000 0.00000000000000 39.098 > 28 34 K 0.00000000000000 0.83333333333333 0.00000000000000 39.098 > 28 35 K 0.33333333333333 0.83333333333333 0.00000000000000 39.098 > 28 36 K 0.66666666666667 0.83333333333333 0.00000000000000 39.098 > 28 37 K 0.00000000000000 0.16666666666667 0.50000000000000 39.098 > 1 38 K 0.33333333333333 0.16666666666667 0.50000000000000 39.098 > 1 39 K 0.66666666666667 0.16666666666667 0.50000000000000 39.098 > 1 40 K 0.00000000000000 0.50000000000000 0.50000000000000 39.098 > 1 41 K 0.33333333333333 0.50000000000000 0.50000000000000 39.098 > 1 42 K 0.66666666666667 0.50000000000000 0.50000000000000 39.098 > 1 43 K 0.00000000000000 0.83333333333333 0.50000000000000 39.098 > 1 44 K 0.33333333333333 0.83333333333333 0.50000000000000 39.098 > 1 45 K 0.66666666666667 0.83333333333333 0.50000000000000 39.098 > 1 46 K 0.16666666666667 0.16666666666667 0.64793473746599 39.098 > 10 47 K 0.50000000000000 0.16666666666667 0.64793473746599 39.098 > 10 48 K 0.83333333333333 0.16666666666667 0.64793473746599 39.098 > 10 49 K 0.16666666666667 0.50000000000000 0.64793473746599 39.098 > 10 50 K 0.50000000000000 0.50000000000000 0.64793473746599 39.098 > 10 51 K 0.83333333333333 0.50000000000000 0.64793473746599 39.098 > 10 52 K 0.16666666666667 0.83333333333333 0.64793473746599 39.098 > 10 53 K 0.50000000000000 0.83333333333333 0.64793473746599 39.098 > 10 54 K 0.83333333333333 0.83333333333333 0.64793473746599 39.098 > 10 55 K 0.00000000000000 0.00000000000000 0.85206526253401 39.098 > 19 56 K 0.33333333333333 0.00000000000000 0.85206526253401 39.098 > 19 57 K 0.66666666666667 0.00000000000000 0.85206526253401 39.098 > 19 58 K 0.00000000000000 0.33333333333333 0.85206526253401 39.098 > 19 59 K 0.33333333333333 0.33333333333333 0.85206526253401 39.098 > 19 60 K 0.66666666666667 0.33333333333333 0.85206526253401 39.098 > 19 61 K 0.00000000000000 0.66666666666667 0.85206526253401 39.098 > 19 62 K 0.33333333333333 0.66666666666667 0.85206526253401 39.098 > 19 63 K 0.66666666666667 0.66666666666667 0.85206526253401 39.098 > 19 64 K 0.16666666666667 0.00000000000000 0.50000000000000 39.098 > 28 65 K 0.50000000000000 0.00000000000000 0.50000000000000 39.098 > 28 66 K 0.83333333333333 0.00000000000000 0.50000000000000 39.098 > 28 67 K 0.16666666666667 0.33333333333333 0.50000000000000 39.098 > 28 68 K 0.50000000000000 0.33333333333333 0.50000000000000 39.098 > 28 69 K 0.83333333333333 0.33333333333333 0.50000000000000 39.098 > 28 70 K 0.16666666666667 0.66666666666667 0.50000000000000 39.098 > 28 71 K 0.50000000000000 0.66666666666667 0.50000000000000 39.098 > 28 72 K 0.83333333333333 0.66666666666667 0.50000000000000 39.098 > 28 73 I 0.00000000000000 0.00000000000000 0.35340590890074 126.904 > 73 74 I 0.33333333333333 0.00000000000000 0.35340590890074 126.904 > 73 75 I 0.66666666666667 0.00000000000000 0.35340590890074 126.904 > 73 76 I 0.00000000000000 0.33333333333333 0.35340590890074 126.904 > 73 77 I 0.33333333333333 0.33333333333333 0.35340590890074 126.904 > 73 78 I 0.66666666666667 0.33333333333333 0.35340590890074 126.904 > 73 79 I 0.00000000000000 0.66666666666667 0.35340590890074 126.904 > 73 80 I 0.33333333333333 0.66666666666667 0.35340590890074 126.904 > 73 81 I 0.66666666666667 0.66666666666667 0.35340590890074 126.904 > 73 82 I 0.16666666666667 0.16666666666667 0.14659409109926 126.904 > 82 83 I 0.50000000000000 0.16666666666667 0.14659409109926 126.904 > 82 84 I 0.83333333333333 0.16666666666667 0.14659409109926 126.904 > 82 85 I 0.16666666666667 0.50000000000000 0.14659409109926 126.904 > 82 86 I 0.50000000000000 0.50000000000000 0.14659409109926 126.904 > 82 87 I 0.83333333333333 0.50000000000000 0.14659409109926 126.904 > 82 88 I 0.16666666666667 0.83333333333333 0.14659409109926 126.904 > 82 89 I 0.50000000000000 0.83333333333333 0.14659409109926 126.904 > 82 90 I 0.83333333333333 0.83333333333333 0.14659409109926 126.904 > 82 91 I 0.16666666666667 0.16666666666667 0.85340590890074 126.904 > 73 92 I 0.50000000000000 0.16666666666667 0.85340590890074 126.904 > 73 93 I 0.83333333333333 0.16666666666667 0.85340590890074 126.904 > 73 94 I 0.16666666666667 0.50000000000000 0.85340590890074 126.904 > 73 95 I 0.50000000000000 0.50000000000000 0.85340590890074 126.904 > 73 96 I 0.83333333333333 0.50000000000000 0.85340590890074 126.904 > 73 97 I 0.16666666666667 0.83333333333333 0.85340590890074 126.904 > 73 98 I 0.50000000000000 0.83333333333333 0.85340590890074 126.904 > 73 99 I 0.83333333333333 0.83333333333333 0.85340590890074 126.904 > 73 100 I 0.00000000000000 0.00000000000000 0.64659409109926 126.904 > 82 101 I 0.33333333333333 0.00000000000000 0.64659409109926 126.904 > 82 102 I 0.66666666666667 0.00000000000000 0.64659409109926 126.904 > 82 103 I 0.00000000000000 0.33333333333333 0.64659409109926 126.904 > 82 104 I 0.33333333333333 0.33333333333333 0.64659409109926 126.904 > 82 105 I 0.66666666666667 0.33333333333333 0.64659409109926 126.904 > 82 106 I 0.00000000000000 0.66666666666667 0.64659409109926 126.904 > 82 107 I 0.33333333333333 0.66666666666667 0.64659409109926 126.904 > 82 108 I 0.66666666666667 0.66666666666667 0.64659409109926 126.904 > 82 109 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 > 109 110 O 0.33333333333333 0.00000000000000 0.00000000000000 15.999 > 109 111 O 0.66666666666667 0.00000000000000 0.00000000000000 15.999 > 109 112 O 0.00000000000000 0.33333333333333 0.00000000000000 15.999 > 109 113 O 0.33333333333333 0.33333333333333 0.00000000000000 15.999 > 109 114 O 0.66666666666667 0.33333333333333 0.00000000000000 15.999 > 109 115 O 0.00000000000000 0.66666666666667 0.00000000000000 15.999 > 109 116 O 0.33333333333333 0.66666666666667 0.00000000000000 15.999 > 109 117 O 0.66666666666667 0.66666666666667 0.00000000000000 15.999 > 109 118 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 109 119 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 109 120 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 109 121 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 109 122 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 109 123 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 109 124 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 109 125 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 109 126 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 3.0759119 0.0000000 0.0000000 0.0000000 3.0759119 0.0000000 0.0000000 0.0000000 3.0891469 -------------------------- Born effective charges -------------------------- 1 K 1.1710065 0.0000000 0.0000000 0.0000000 1.1424984 0.0000000 0.0000000 0.0000000 0.9864279 2 K 1.1290615 0.0000000 0.0000000 0.0000000 1.1290615 0.0000000 0.0000000 0.0000000 1.2349842 3 K 1.1290615 0.0000000 0.0000000 0.0000000 1.1290615 0.0000000 0.0000000 0.0000000 1.2349842 4 K 1.1424984 0.0000000 0.0000000 0.0000000 1.1710065 0.0000000 0.0000000 0.0000000 0.9864279 5 I -1.1615912 0.0000000 0.0000000 0.0000000 -1.1615912 0.0000000 0.0000000 0.0000000 -1.1299427 6 I -1.1615912 0.0000000 0.0000000 0.0000000 -1.1615912 0.0000000 0.0000000 0.0000000 -1.1299427 7 O -2.2484455 0.0000000 0.0000000 0.0000000 -2.2484455 0.0000000 0.0000000 0.0000000 -2.1829389 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz) Max drift of fc2: -0.00000000 (zz) -0.00000000 (zz) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... Generalized regular grid: [ 3 9 18 ] Grid generation matrix: [ 0 9 9 ] [ 9 0 9 ] [ 3 3 0 ] 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.66, Number of G-points: 295, Lambda: 0.12 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/60) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 60 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.219 ( 0.000 0.000 0.000) 0.000 1.219 ( 0.000 0.000 0.000) 0.000 1.840 ( 0.000 0.000 0.000) 0.000 1.921 ( 0.000 0.000 0.000) 0.000 1.921 ( 0.000 0.000 0.000) 0.000 2.275 ( 0.000 0.000 0.000) 0.000 2.275 ( 0.000 0.000 0.000) 0.000 2.361 ( 0.000 0.000 0.000) 0.000 2.361 ( 0.000 0.000 0.000) 0.000 2.413 ( 0.000 0.000 0.000) 0.000 2.539 ( 0.000 0.000 0.000) 0.000 4.119 ( 0.000 0.000 0.000) 0.000 4.119 ( 0.000 0.000 0.000) 0.000 4.329 ( 0.000 0.000 0.000) 0.000 4.510 ( 0.000 0.000 0.000) 0.000 7.256 ( 0.000 0.000 0.000) 0.000 7.256 ( 0.000 0.000 0.000) 0.000 8.488 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/60) ======================= q-point: ( 0.00 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.931 ( 12.071 0.000 6.191) 13.567 0.992 ( 11.121 0.000 1.391) 11.207 1.275 ( 1.421 0.000 -0.477) 1.499 1.402 ( 2.924 0.000 -0.150) 2.928 1.462 ( 8.914 0.000 4.125) 9.822 1.640 ( -4.114 -0.000 -2.964) 5.071 1.781 ( -4.045 -0.000 0.745) 4.113 1.895 ( -4.063 -0.000 0.707) 4.124 2.295 ( 0.478 0.000 -2.995) 3.033 2.309 ( -0.840 -0.000 -0.647) 1.060 2.410 ( 0.725 0.000 0.991) 1.228 2.575 ( 6.651 0.000 0.376) 6.661 3.104 ( 11.386 0.000 -7.455) 13.610 3.360 ( 10.502 0.000 4.881) 11.581 3.950 ( -7.021 -0.000 -3.301) 7.758 3.951 ( -3.958 -0.000 0.039) 3.958 4.174 ( -5.643 -0.000 3.156) 6.466 5.376 ( 15.977 0.000 -1.402) 16.038 7.308 ( 1.024 0.000 0.031) 1.025 8.331 ( -2.724 -0.000 5.207) 5.877 9.013 ( 6.927 0.000 -2.878) 7.501 ======================= Grid point 3 (3/60) ======================= q-point: ( 0.00 0.00 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.354 ( 9.123 9.123 0.000) 12.902 0.569 ( 11.688 11.688 0.000) 16.529 0.849 ( 16.788 16.788 0.000) 23.742 1.256 ( 1.710 1.710 0.000) 2.418 1.281 ( 2.355 2.355 0.000) 3.330 1.783 ( -2.134 -2.134 0.000) 3.019 2.002 ( 3.184 3.184 0.000) 4.502 2.148 ( -6.159 -6.159 0.000) 8.710 2.248 ( -1.352 -1.352 0.000) 1.912 2.254 ( -3.452 -3.452 0.000) 4.882 2.302 ( 2.917 2.917 0.000) 4.125 2.533 ( 4.078 4.078 0.000) 5.768 2.680 ( 11.118 11.118 0.000) 15.724 2.765 ( 6.781 6.781 0.000) 9.589 4.223 ( 4.034 4.034 0.000) 5.704 4.254 ( -5.576 -5.576 0.000) 7.886 4.273 ( -2.700 -2.700 0.000) 3.818 4.826 ( 8.046 8.046 0.000) 11.378 7.353 ( 5.096 5.096 0.000) 7.207 8.416 ( -2.817 -2.817 0.000) 3.983 8.689 ( 1.397 1.397 0.000) 1.976 ======================= Grid point 4 (4/60) ======================= q-point: ( 0.00 0.33 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.173 ( 7.577 0.753 0.747) 7.651 1.211 ( 6.963 1.859 1.968) 7.471 1.322 ( 0.737 2.159 -1.009) 2.494 1.382 ( -4.095 0.654 1.129) 4.298 1.532 ( -1.567 -0.858 -0.012) 1.787 1.614 ( 2.003 -0.935 -1.869) 2.895 1.817 ( -1.002 1.870 0.591) 2.202 1.856 ( -1.070 6.951 1.311) 7.154 2.078 ( -1.597 -9.614 -1.185) 9.817 2.276 ( 0.965 -9.590 -0.023) 9.639 2.533 ( -0.141 8.729 -0.842) 8.771 2.685 ( 1.527 -0.125 0.455) 1.598 3.340 ( 3.726 6.833 -2.514) 8.179 3.531 ( 7.796 -1.438 5.282) 9.526 3.791 ( -6.336 -2.675 -6.269) 9.306 4.032 ( -2.258 3.700 2.319) 4.916 4.061 ( -2.379 6.235 0.627) 6.703 5.566 ( 5.342 -2.445 -0.337) 5.885 7.537 ( 0.993 18.350 -0.017) 18.377 8.307 ( 0.196 1.275 1.498) 1.977 9.069 ( 2.824 -4.883 -0.513) 5.664 ======================= Grid point 5 (5/60) ======================= q-point: ( 0.00 -0.33 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.741 ( -8.506 3.702 -12.501) 15.567 0.837 ( -10.711 7.943 -2.962) 13.660 1.208 ( -12.303 4.829 -1.034) 13.257 1.276 ( -2.702 2.704 0.977) 3.945 1.368 ( -3.986 1.227 -2.201) 4.716 1.702 ( 2.752 -1.182 2.488) 3.894 1.961 ( 4.994 2.255 1.071) 5.583 2.007 ( 5.080 4.656 0.099) 6.891 2.159 ( 1.548 -6.355 3.618) 7.475 2.223 ( -0.318 -4.590 -0.201) 4.605 2.392 ( -8.352 -2.968 3.713) 9.610 2.609 ( -1.132 10.435 -2.090) 10.703 2.899 ( -12.018 7.891 4.604) 15.096 3.063 ( -15.121 -0.266 -2.306) 15.298 4.077 ( 6.972 -3.010 1.664) 7.774 4.209 ( 3.380 11.061 0.233) 11.568 4.267 ( 5.620 -2.447 -4.201) 7.431 5.010 ( -16.875 2.785 5.129) 17.855 7.438 ( -3.266 12.819 -0.031) 13.229 8.308 ( -7.475 1.058 1.818) 7.765 8.862 ( 2.519 -6.635 -4.164) 8.228 ======================= Grid point 6 (6/60) ======================= q-point: ( 0.00 0.00 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.754 ( 9.548 9.548 0.000) 13.503 0.969 ( 7.508 7.508 0.000) 10.617 1.347 ( 2.327 2.327 0.000) 3.291 1.349 ( 7.076 7.076 0.000) 10.007 1.370 ( 1.659 1.659 0.000) 2.347 1.697 ( -1.738 -1.738 0.000) 2.458 1.915 ( -4.711 -4.711 0.000) 6.662 2.084 ( -1.344 -1.344 0.000) 1.900 2.084 ( -4.624 -4.624 0.000) 6.539 2.199 ( 0.537 0.537 0.000) 0.760 2.416 ( 1.973 1.973 0.000) 2.790 2.731 ( 4.991 4.991 0.000) 7.058 3.031 ( 5.532 5.532 0.000) 7.823 3.240 ( 13.246 13.246 0.000) 18.732 4.057 ( -3.723 -3.723 0.000) 5.265 4.090 ( -6.139 -6.139 0.000) 8.681 4.350 ( 1.003 1.003 0.000) 1.418 5.220 ( 11.071 11.071 0.000) 15.657 7.721 ( 12.234 12.234 0.000) 17.301 8.319 ( -1.107 -1.107 0.000) 1.566 8.748 ( 1.159 1.159 0.000) 1.639 ======================= Grid point 7 (7/60) ======================= q-point: ( 0.00 -0.67 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.243 ( -2.816 2.896 -1.938) 4.480 1.299 ( -4.976 2.728 -1.552) 5.883 1.330 ( 0.483 1.801 2.074) 2.789 1.421 ( 2.049 0.282 -1.698) 2.676 1.503 ( 1.479 -0.040 1.850) 2.369 1.592 ( -0.565 -1.391 -2.065) 2.553 1.779 ( -0.484 -3.243 0.931) 3.409 1.870 ( 0.505 -7.321 1.584) 7.507 2.053 ( 0.853 -10.056 0.078) 10.093 2.075 ( 1.536 9.089 -0.056) 9.218 2.549 ( -2.998 -4.576 -0.375) 5.483 2.806 ( 1.652 9.937 0.326) 10.079 3.433 ( -8.565 1.700 -1.519) 8.863 3.571 ( -6.569 0.533 1.304) 6.719 3.631 ( 7.944 -4.539 2.555) 9.499 4.010 ( 3.685 -2.334 -1.118) 4.503 4.259 ( 1.606 4.448 -0.062) 4.729 5.628 ( -5.509 10.598 0.441) 11.953 7.997 ( -1.859 23.360 0.023) 23.434 8.377 ( 2.601 1.444 -2.077) 3.628 8.933 ( -4.372 -5.757 0.912) 7.286 ======================= Grid point 9 (8/60) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.124 ( 7.633 7.633 0.000) 10.795 1.212 ( 4.178 4.178 0.000) 5.909 1.417 ( 0.594 0.594 0.000) 0.840 1.424 ( 1.093 1.093 0.000) 1.545 1.490 ( 0.561 0.561 0.000) 0.793 1.642 ( -1.056 -1.056 0.000) 1.494 1.763 ( -2.432 -2.432 0.000) 3.439 1.869 ( -5.442 -5.442 0.000) 7.696 1.922 ( -5.155 -5.155 0.000) 7.290 2.256 ( 1.403 1.403 0.000) 1.985 2.444 ( -0.610 -0.610 0.000) 0.862 2.928 ( 4.366 4.366 0.000) 6.174 3.210 ( 2.859 2.859 0.000) 4.044 3.650 ( 5.511 5.511 0.000) 7.794 3.768 ( -8.764 -8.764 0.000) 12.394 3.952 ( -0.996 -0.996 0.000) 1.408 4.300 ( -2.731 -2.731 0.000) 3.862 5.743 ( 12.811 12.811 0.000) 18.118 8.293 ( 13.387 13.387 0.000) 18.932 8.341 ( 1.803 1.803 0.000) 2.550 8.774 ( 0.065 0.065 0.000) 0.092 ======================= Grid point 10 (9/60) ======================= q-point: ( 0.00 -0.67 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.247 ( -5.039 5.039 -3.958) 8.152 1.262 ( -3.287 3.287 -2.443) 5.251 1.370 ( -2.059 2.059 2.558) 3.877 1.391 ( -0.370 0.370 0.593) 0.791 1.543 ( 0.699 -0.699 -2.204) 2.415 1.599 ( 0.522 -0.522 0.410) 0.844 1.688 ( 2.786 -2.786 2.058) 4.445 1.818 ( 3.390 -3.390 2.858) 5.581 1.931 ( 5.719 -5.719 -0.072) 8.088 2.316 ( -1.719 1.719 -1.786) 3.017 2.342 ( -0.851 0.851 2.654) 2.914 3.045 ( -2.705 2.705 -3.943) 5.494 3.187 ( -2.767 2.767 0.738) 3.983 3.551 ( -7.786 7.786 5.663) 12.382 3.669 ( 7.205 -7.205 1.602) 10.314 4.041 ( 3.852 -3.852 -2.619) 6.044 4.298 ( 2.735 -2.735 0.074) 3.868 5.713 ( -14.094 14.094 0.884) 19.951 8.294 ( -13.119 13.119 -0.039) 18.553 8.425 ( 0.771 -0.771 -3.459) 3.627 8.757 ( -1.448 1.448 1.311) 2.432 ======================= Grid point 12 (10/60) ======================= q-point: ( 0.00 0.00 0.44) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.334 ( 1.674 1.674 0.000) 2.368 1.345 ( 1.805 1.805 0.000) 2.552 1.441 ( -0.007 -0.007 0.000) 0.010 1.460 ( 2.751 2.751 0.000) 3.891 1.477 ( -0.509 -0.509 0.000) 0.720 1.601 ( -0.679 -0.679 0.000) 0.960 1.652 ( -4.497 -4.497 0.000) 6.360 1.704 ( -0.586 -0.586 0.000) 0.829 1.719 ( -3.453 -3.453 0.000) 4.883 2.304 ( 0.719 0.719 0.000) 1.017 2.397 ( -0.984 -0.984 0.000) 1.391 3.122 ( 5.274 5.274 0.000) 7.459 3.280 ( 0.683 0.683 0.000) 0.965 3.404 ( -7.924 -7.924 0.000) 11.207 3.736 ( 0.083 0.083 0.000) 0.118 3.967 ( 0.796 0.796 0.000) 1.125 4.189 ( -1.760 -1.760 0.000) 2.489 6.169 ( 5.933 5.933 0.000) 8.390 8.419 ( 1.283 1.283 0.000) 1.815 8.718 ( 5.674 5.674 0.000) 8.024 8.764 ( -0.323 -0.323 0.000) 0.456 ======================= Grid point 27 (11/60) ======================= q-point: ( 0.06 -0.06 -0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.249 ( -12.874 -0.000 0.000) 12.874 0.361 ( -16.756 -0.000 0.000) 16.756 0.649 ( -27.010 -0.000 0.000) 27.010 1.228 ( -0.835 -0.000 0.000) 0.835 1.263 ( -3.746 -0.000 0.000) 3.746 1.809 ( 2.499 0.000 -0.000) 2.499 1.903 ( 1.728 0.000 -0.000) 1.728 2.208 ( 9.000 0.000 -0.000) 9.000 2.249 ( -2.011 -0.000 0.000) 2.011 2.316 ( -3.782 -0.000 0.000) 3.782 2.355 ( 0.594 0.000 -0.000) 0.594 2.423 ( 0.565 0.000 -0.000) 0.565 2.625 ( -16.208 -0.000 0.000) 16.208 2.701 ( -11.309 -0.000 0.000) 11.309 4.097 ( 2.108 0.000 -0.000) 2.108 4.302 ( 2.466 0.000 -0.000) 2.466 4.325 ( 8.759 0.000 -0.000) 8.759 4.743 ( -12.358 -0.000 0.000) 12.358 7.264 ( -0.724 -0.000 0.000) 0.724 8.448 ( 3.498 0.000 -0.000) 3.498 8.718 ( -5.679 -0.000 0.000) 5.679 ======================= Grid point 28 (12/60) ======================= q-point: ( 0.06 0.28 -0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.703 ( 9.123 0.000 13.500) 16.293 0.723 ( 14.105 0.000 2.600) 14.343 1.156 ( 16.877 0.000 3.528) 17.242 1.241 ( 1.600 0.000 -0.976) 1.874 1.356 ( 4.320 0.000 1.226) 4.491 1.718 ( -3.107 -0.000 -2.627) 4.069 1.860 ( -3.256 -0.000 0.940) 3.389 2.013 ( -6.861 -0.000 -1.722) 7.073 2.257 ( 3.021 0.000 -2.873) 4.169 2.327 ( -0.797 -0.000 -1.003) 1.281 2.407 ( -0.089 -0.000 -0.041) 0.098 2.431 ( 6.502 0.000 0.353) 6.511 2.819 ( 14.856 0.000 -8.553) 17.142 3.090 ( 14.658 0.000 3.184) 15.000 4.034 ( -3.679 -0.000 0.020) 3.679 4.110 ( -7.741 -0.000 -1.403) 7.867 4.296 ( -5.879 -0.000 4.369) 7.325 4.996 ( 19.044 0.000 -4.427) 19.551 7.284 ( 1.140 0.000 0.029) 1.140 8.297 ( 10.163 0.000 -3.342) 10.699 8.950 ( -5.195 -0.000 4.884) 7.130 ======================= Grid point 29 (13/60) ======================= q-point: ( 0.06 -0.39 -0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.171 ( -5.170 -0.000 -0.445) 5.189 1.187 ( -11.582 -0.000 -1.505) 11.679 1.296 ( -0.471 -0.000 0.138) 0.491 1.375 ( 4.736 0.000 -1.356) 4.926 1.545 ( 2.304 0.000 0.952) 2.493 1.631 ( -3.891 -0.000 0.927) 4.000 1.708 ( 2.245 0.000 -0.298) 2.264 1.847 ( 0.737 0.000 -1.806) 1.951 2.288 ( 0.736 -0.000 2.244) 2.361 2.294 ( 0.490 0.000 0.229) 0.541 2.429 ( -0.691 -0.000 -1.048) 1.255 2.684 ( -3.073 -0.000 -0.170) 3.077 3.256 ( -1.961 -0.000 2.935) 3.530 3.549 ( -8.416 -0.000 -5.624) 10.122 3.819 ( 5.816 0.000 6.790) 8.940 3.885 ( 1.845 0.000 -0.015) 1.845 4.077 ( 2.730 0.000 -3.232) 4.231 5.618 ( -6.089 -0.000 0.242) 6.094 7.324 ( -0.419 -0.000 -0.014) 0.419 8.293 ( 0.014 0.000 -1.614) 1.614 9.124 ( -2.495 -0.000 0.701) 2.591 ======================= Grid point 31 (14/60) ======================= q-point: ( 0.06 0.28 0.06) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.952 ( 11.089 2.315 5.803) 12.728 1.051 ( 9.173 4.292 2.077) 10.338 1.299 ( 1.506 2.176 -1.015) 2.834 1.403 ( 2.537 0.070 -0.344) 2.561 1.473 ( 6.840 0.911 3.642) 7.803 1.634 ( -3.521 -0.344 -2.988) 4.631 1.863 ( -3.544 -0.396 -0.101) 3.567 1.912 ( -3.962 8.239 0.967) 9.194 2.118 ( -2.071 -8.257 -2.594) 8.899 2.247 ( 1.584 -8.042 0.445) 8.209 2.518 ( 2.766 3.743 -2.138) 5.122 2.631 ( 2.305 5.622 1.122) 6.179 3.159 ( 11.606 4.886 -4.495) 13.371 3.346 ( 10.997 0.269 3.219) 11.462 3.927 ( -6.928 -2.124 -3.449) 8.025 4.119 ( -4.822 9.176 1.172) 10.432 4.146 ( -5.430 0.402 1.998) 5.800 5.351 ( 14.323 -0.473 -1.728) 14.435 7.499 ( 2.481 16.591 -0.022) 16.775 8.334 ( -2.038 0.271 4.803) 5.225 8.950 ( 7.368 -5.348 -2.237) 9.375 ======================= Grid point 33 (15/60) ======================= q-point: ( 0.06 -0.06 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.597 ( 5.105 13.289 0.000) 14.236 0.807 ( 7.866 11.232 0.000) 13.712 1.194 ( 6.674 14.327 0.000) 15.805 1.284 ( 2.794 1.249 -0.000) 3.060 1.354 ( 0.372 3.638 0.000) 3.657 1.729 ( -1.312 -2.862 -0.000) 3.148 1.978 ( 9.758 -3.554 -0.000) 10.385 1.997 ( -2.775 -7.240 -0.000) 7.753 2.199 ( -6.598 -2.451 0.000) 7.038 2.220 ( -4.504 0.130 0.000) 4.505 2.453 ( -1.540 7.999 0.000) 8.146 2.571 ( 9.279 0.265 -0.000) 9.283 2.969 ( 1.194 11.530 0.000) 11.592 3.010 ( 5.662 16.794 0.000) 17.723 4.118 ( -1.581 -6.527 -0.000) 6.715 4.191 ( -3.357 -5.073 -0.000) 6.083 4.222 ( 10.527 -3.070 -0.000) 10.965 5.067 ( 3.114 14.102 0.000) 14.441 7.438 ( 12.749 3.316 -0.000) 13.173 8.345 ( -1.144 -3.499 -0.000) 3.681 8.799 ( -4.959 7.820 0.000) 9.260 ======================= Grid point 34 (16/60) ======================= q-point: ( 0.06 0.28 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.211 ( 6.436 2.891 1.235) 7.163 1.258 ( 4.977 2.470 2.077) 5.932 1.371 ( 0.390 2.200 -1.920) 2.946 1.400 ( -2.537 1.009 0.639) 2.804 1.514 ( -0.906 -0.784 0.839) 1.462 1.597 ( -0.276 -0.670 -2.570) 2.671 1.781 ( -0.083 -2.756 0.667) 2.836 1.898 ( -2.854 -7.455 -1.091) 8.057 2.046 ( -0.254 -7.306 0.783) 7.352 2.077 ( -1.613 6.542 -0.196) 6.741 2.562 ( 1.743 -4.498 -1.696) 5.113 2.796 ( -0.547 9.301 1.388) 9.420 3.414 ( 9.568 -2.387 0.558) 9.877 3.562 ( -0.133 5.322 1.011) 5.419 3.709 ( -6.355 -4.318 -4.004) 8.664 3.982 ( -1.473 -2.706 1.638) 3.490 4.258 ( -1.560 4.429 0.027) 4.696 5.637 ( 4.561 10.499 -0.432) 11.454 7.997 ( 1.827 23.361 -0.036) 23.432 8.343 ( 0.665 2.096 1.197) 2.503 8.944 ( 3.202 -6.303 -0.155) 7.071 ======================= Grid point 35 (17/60) ======================= q-point: ( 0.06 -0.39 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.857 ( -7.226 7.226 -9.432) 13.907 0.996 ( -7.154 7.154 -2.617) 10.450 1.325 ( -4.760 4.760 1.325) 6.861 1.327 ( -2.339 2.339 2.002) 3.866 1.418 ( -3.547 3.547 -3.869) 6.335 1.672 ( 1.639 -1.639 2.192) 3.190 1.890 ( 4.691 -4.691 2.086) 6.954 2.050 ( 4.092 -4.092 3.528) 6.777 2.107 ( 1.943 -1.943 -1.879) 3.329 2.205 ( -1.051 1.051 -0.652) 1.623 2.343 ( -2.740 2.740 5.612) 6.820 2.789 ( -5.075 5.075 -4.498) 8.470 3.029 ( -5.376 5.376 0.232) 7.606 3.211 ( -12.764 12.764 3.229) 18.338 3.986 ( 5.536 -5.536 2.767) 8.303 4.180 ( 5.106 -5.106 -4.682) 8.606 4.349 ( -1.020 1.020 0.084) 1.445 5.173 ( -12.386 12.386 4.288) 18.033 7.723 ( -12.236 12.236 -0.204) 17.306 8.338 ( -1.306 1.306 -2.412) 3.038 8.761 ( 1.171 -1.171 -0.553) 1.746 ======================= Grid point 36 (18/60) ======================= q-point: ( 0.06 -0.06 0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.978 ( 5.823 11.158 0.000) 12.586 1.121 ( 4.398 6.675 0.000) 7.994 1.373 ( 2.975 0.823 -0.000) 3.086 1.413 ( 0.334 1.174 0.000) 1.221 1.467 ( 1.440 3.835 0.000) 4.097 1.659 ( -0.597 -1.673 -0.000) 1.777 1.819 ( -2.618 -3.878 -0.000) 4.679 1.992 ( -6.408 -4.206 0.000) 7.665 2.049 ( -4.880 -1.196 0.000) 5.024 2.139 ( 6.796 -3.941 -0.000) 7.856 2.521 ( -5.212 5.154 0.000) 7.330 2.777 ( 9.311 0.212 -0.000) 9.313 3.199 ( -0.948 9.613 0.000) 9.659 3.477 ( 10.016 8.059 -0.000) 12.856 3.939 ( -7.280 -8.125 -0.000) 10.909 3.981 ( -1.373 -2.868 -0.000) 3.180 4.314 ( 3.169 -3.094 -0.000) 4.429 5.472 ( 10.569 11.262 0.000) 15.445 7.915 ( 20.812 6.302 -0.000) 21.745 8.308 ( 1.267 0.224 -0.000) 1.287 8.838 ( -5.475 6.106 0.000) 8.202 ======================= Grid point 37 (19/60) ======================= q-point: ( 0.06 -0.72 0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.307 ( -1.179 2.854 -1.388) 3.385 1.349 ( -2.300 1.633 -1.284) 3.100 1.376 ( -0.978 2.529 1.241) 2.982 1.424 ( -0.735 0.262 -0.388) 0.872 1.509 ( 2.237 0.706 0.831) 2.488 1.574 ( 0.399 -0.284 -1.256) 1.348 1.691 ( -2.212 -4.634 1.346) 5.308 1.756 ( 2.379 -3.809 0.821) 4.566 1.848 ( 2.093 -8.288 0.011) 8.548 2.253 ( 2.109 6.711 0.309) 7.041 2.437 ( -2.848 -5.319 -0.078) 6.034 3.008 ( 2.328 9.012 -0.768) 9.339 3.345 ( -8.729 -6.151 -1.446) 10.776 3.502 ( 7.669 -5.922 1.044) 9.746 3.689 ( -4.524 5.528 5.193) 8.832 3.967 ( 2.885 -1.820 -3.012) 4.551 4.249 ( 1.349 -3.278 0.011) 3.545 5.930 ( -5.667 15.058 0.247) 16.091 8.408 ( 0.798 1.522 -1.159) 2.073 8.448 ( -2.590 17.939 0.083) 18.125 8.826 ( -3.073 -4.077 0.222) 5.111 ======================= Grid point 39 (20/60) ======================= q-point: ( 0.06 -0.06 0.39) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.275 ( 4.131 5.954 0.000) 7.246 1.291 ( 2.432 3.170 0.000) 3.996 1.420 ( 1.171 -0.801 -0.000) 1.419 1.433 ( 1.225 0.611 -0.000) 1.369 1.484 ( -0.460 -1.324 -0.000) 1.401 1.617 ( -0.851 -1.170 -0.000) 1.447 1.726 ( -1.760 -1.249 0.000) 2.159 1.765 ( -5.928 -4.052 0.000) 7.181 1.846 ( -8.195 -1.919 0.000) 8.417 2.236 ( 5.937 -1.133 -0.000) 6.044 2.465 ( -5.308 0.764 0.000) 5.362 2.972 ( 8.531 1.062 -0.000) 8.596 3.324 ( -3.882 7.136 0.000) 8.124 3.577 ( -8.412 -9.302 -0.000) 12.541 3.691 ( 4.466 -1.217 -0.000) 4.629 3.956 ( 0.088 1.412 0.000) 1.415 4.250 ( -3.187 -1.426 0.000) 3.492 5.939 ( 14.757 4.833 -0.000) 15.528 8.381 ( 2.203 1.948 -0.000) 2.941 8.449 ( 18.275 2.540 -0.000) 18.451 8.825 ( -4.548 2.809 0.000) 5.345 ======================= Grid point 42 (21/60) ======================= q-point: ( 0.06 0.94 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.353 ( 0.539 -0.539 -0.000) 0.762 1.361 ( -0.751 0.751 0.000) 1.063 1.422 ( 1.752 -1.752 -0.000) 2.478 1.447 ( 2.040 -2.040 -0.000) 2.885 1.522 ( -0.886 0.886 0.000) 1.253 1.592 ( 0.078 -0.078 -0.000) 0.111 1.593 ( -0.730 0.730 0.000) 1.032 1.706 ( -0.858 0.858 0.000) 1.213 1.719 ( -3.451 3.451 0.000) 4.881 2.317 ( -0.442 0.442 0.000) 0.625 2.381 ( 0.512 -0.512 -0.000) 0.724 3.179 ( 2.714 -2.714 -0.000) 3.838 3.257 ( 0.951 -0.951 -0.000) 1.345 3.369 ( -5.584 5.584 0.000) 7.897 3.710 ( 2.501 -2.501 -0.000) 3.537 3.983 ( -0.730 0.730 0.000) 1.032 4.189 ( -1.682 1.682 0.000) 2.378 6.166 ( 6.257 -6.257 -0.000) 8.849 8.432 ( 0.009 -0.009 -0.000) 0.013 8.703 ( 5.373 -5.373 -0.000) 7.599 8.779 ( 0.007 -0.007 -0.000) 0.010 ======================= Grid point 45 (22/60) ======================= q-point: ( 0.06 0.61 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.196 ( 1.492 -4.017 1.424) 4.516 1.246 ( 2.149 -6.512 1.170) 6.957 1.302 ( 0.894 -0.125 -0.903) 1.276 1.411 ( 0.520 5.110 1.822) 5.450 1.506 ( -0.236 0.339 -1.967) 2.010 1.627 ( -1.707 -1.175 2.555) 3.289 1.820 ( 1.557 1.083 -0.436) 1.946 1.861 ( 6.966 1.426 -0.769) 7.152 2.056 ( -10.158 -0.419 -0.863) 10.203 2.282 ( -9.547 -0.435 0.568) 9.574 2.534 ( 8.771 0.170 0.950) 8.824 2.675 ( -0.284 -2.335 -1.428) 2.751 3.288 ( 8.066 -7.405 -2.533) 11.239 3.617 ( -2.406 -8.495 1.695) 8.990 3.722 ( -3.270 8.032 0.730) 8.703 4.048 ( 12.434 2.410 0.155) 12.667 4.061 ( -1.811 3.312 -1.170) 3.952 5.557 ( -2.558 -6.221 -0.461) 6.742 7.537 ( 18.315 -1.015 -0.012) 18.343 8.353 ( 0.826 4.224 2.940) 5.212 9.046 ( -4.397 -5.125 -1.679) 6.958 ======================= Grid point 48 (23/60) ======================= q-point: ( 0.06 0.61 -0.39) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.141 ( 4.539 -5.326 5.073) 8.643 1.191 ( 3.792 -5.720 2.876) 7.441 1.316 ( 2.444 -1.488 -2.428) 3.753 1.392 ( -0.043 -3.447 -0.639) 3.506 1.558 ( -0.783 -0.651 3.879) 4.010 1.589 ( 0.728 2.486 -1.039) 2.791 1.781 ( -4.815 1.464 -1.280) 5.193 1.903 ( -5.188 2.162 -4.279) 7.064 2.081 ( -8.289 1.954 0.560) 8.535 2.137 ( 10.370 4.019 -0.337) 11.126 2.436 ( -4.863 -7.643 -0.934) 9.107 2.862 ( 11.380 3.234 1.589) 11.936 3.220 ( 0.102 -10.984 1.630) 11.104 3.377 ( 6.031 -11.268 -5.541) 13.930 3.810 ( -5.874 8.168 -0.568) 10.076 4.117 ( -3.211 6.031 1.515) 6.999 4.312 ( 3.389 2.897 0.031) 4.458 5.415 ( 11.390 -14.089 -1.653) 18.193 7.916 ( 20.982 -6.013 0.061) 21.827 8.444 ( -1.029 2.941 6.096) 6.847 8.786 ( -3.536 -7.869 -3.420) 9.280 ======================= Grid point 53 (24/60) ======================= q-point: ( 0.06 0.28 -0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.603 ( 4.847 -4.847 18.020) 19.280 0.612 ( 11.095 -11.095 4.087) 16.214 0.898 ( 14.330 -14.330 5.456) 20.987 1.238 ( 1.659 -1.659 -1.806) 2.961 1.289 ( 2.967 -2.967 1.084) 4.333 1.750 ( -1.826 1.826 -3.085) 4.023 2.010 ( 2.760 -2.760 0.778) 3.980 2.136 ( -5.595 5.595 -1.147) 7.995 2.143 ( 2.603 -2.603 -8.197) 8.986 2.253 ( -1.319 1.319 0.470) 1.923 2.257 ( -1.677 1.677 -3.501) 4.229 2.501 ( 8.905 -8.905 -1.122) 12.644 2.758 ( 5.464 -5.464 2.019) 7.986 2.766 ( 6.843 -6.843 0.033) 9.678 4.200 ( -4.290 4.290 -1.715) 6.304 4.221 ( 4.024 -4.024 -0.136) 5.692 4.365 ( -3.156 3.156 4.967) 6.678 4.704 ( 10.015 -10.015 -8.778) 16.663 7.354 ( 5.133 -5.133 0.104) 7.260 8.031 ( 15.318 -15.318 -16.193) 27.046 9.072 ( -14.764 14.764 16.082) 26.355 ======================= Grid point 55 (25/60) ======================= q-point: ( 0.11 -0.11 -0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.544 ( -14.954 -0.000 0.000) 14.954 0.696 ( -14.868 -0.000 0.000) 14.868 1.117 ( -17.909 -0.000 0.000) 17.909 1.251 ( -1.338 -0.000 0.000) 1.338 1.350 ( -3.947 -0.000 0.000) 3.947 1.747 ( 3.299 0.000 -0.000) 3.299 1.850 ( 3.221 0.000 -0.000) 3.221 2.026 ( 8.100 0.000 -0.000) 8.100 2.303 ( -2.960 -0.000 0.000) 2.960 2.337 ( 0.985 0.000 -0.000) 0.985 2.401 ( 1.434 0.000 -0.000) 1.434 2.427 ( -6.483 -0.000 0.000) 6.483 2.958 ( -12.879 -0.000 0.000) 12.879 3.000 ( -16.713 -0.000 0.000) 16.713 4.034 ( 3.702 0.000 -0.000) 3.702 4.150 ( 7.570 0.000 -0.000) 7.570 4.227 ( 4.646 0.000 -0.000) 4.646 5.047 ( -16.294 -0.000 0.000) 16.294 7.284 ( -1.134 -0.000 0.000) 1.134 8.358 ( 4.494 0.000 -0.000) 4.494 8.872 ( -8.308 -0.000 0.000) 8.308 ======================= Grid point 56 (26/60) ======================= q-point: ( 0.11 0.22 -0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.415 ( 14.539 0.000 4.899) 15.342 0.543 ( 9.080 0.000 16.066) 18.455 0.746 ( 19.114 0.000 12.733) 22.966 1.213 ( 1.026 0.000 -1.465) 1.789 1.259 ( 4.412 0.000 -0.227) 4.417 1.773 ( -2.093 -0.000 -3.399) 3.992 1.913 ( -1.706 -0.000 0.946) 1.951 2.122 ( 5.612 0.000 -7.049) 9.010 2.205 ( -5.305 -0.000 -1.964) 5.657 2.314 ( 3.465 0.000 -0.165) 3.469 2.345 ( -0.130 0.000 -0.867) 0.876 2.349 ( 6.532 0.000 -1.584) 6.721 2.596 ( 2.553 0.000 -3.260) 4.141 2.800 ( 10.802 0.000 1.544) 10.911 4.096 ( -2.080 -0.000 -0.037) 2.081 4.246 ( -4.580 -0.000 -1.973) 4.987 4.395 ( 1.503 0.000 1.319) 1.999 4.623 ( 12.238 0.000 -3.939) 12.856 7.264 ( 0.732 0.000 0.014) 0.732 7.839 ( 32.588 0.000 -19.302) 37.875 9.289 ( -25.834 -0.000 16.932) 30.889 ======================= Grid point 57 (27/60) ======================= q-point: ( 0.11 0.56 -0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.180 ( 4.242 0.000 0.473) 4.268 1.220 ( 7.979 0.000 1.573) 8.133 1.293 ( 0.786 0.000 -0.164) 0.803 1.406 ( -6.674 -0.000 1.519) 6.845 1.510 ( 0.176 0.000 -1.752) 1.761 1.655 ( 1.722 0.000 2.798) 3.285 1.714 ( -2.814 -0.000 0.297) 2.830 1.851 ( -0.863 -0.000 -1.496) 1.727 2.276 ( 0.286 -0.000 1.037) 1.076 2.290 ( -0.036 0.000 -0.229) 0.232 2.429 ( 0.751 0.000 -0.974) 1.230 2.687 ( 2.761 0.000 0.168) 2.766 3.188 ( 7.551 0.000 -3.568) 8.352 3.644 ( 7.269 0.000 1.545) 7.431 3.757 ( -7.139 -0.000 2.327) 7.509 3.885 ( -1.876 -0.000 0.015) 1.876 4.086 ( -3.386 -0.000 -2.001) 3.933 5.611 ( 6.821 0.000 -0.422) 6.834 7.324 ( 0.393 0.000 0.014) 0.393 8.344 ( -4.917 -0.000 3.310) 5.927 9.096 ( 5.359 0.000 -2.041) 5.734 ======================= Grid point 62 (28/60) ======================= q-point: ( 0.11 0.22 0.11) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.033 ( 9.070 5.266 5.240) 11.724 1.140 ( 6.440 4.115 1.998) 7.899 1.356 ( 1.013 2.900 -1.610) 3.468 1.405 ( 2.060 0.146 -0.677) 2.173 1.494 ( 2.929 0.920 2.542) 3.986 1.629 ( -2.399 -0.101 -2.550) 3.503 1.806 ( -2.705 -3.603 -1.213) 4.666 1.970 ( -3.713 -5.952 -2.399) 7.414 2.064 ( -1.485 -6.357 1.171) 6.632 2.140 ( -4.084 8.481 0.030) 9.413 2.477 ( 6.190 -4.641 -3.175) 8.363 2.816 ( -0.900 9.433 2.949) 9.924 3.203 ( 10.000 -0.977 0.420) 10.056 3.457 ( 9.277 9.162 -2.439) 13.264 3.853 ( -6.871 -5.004 -4.173) 9.469 4.063 ( -5.439 -3.965 4.141) 7.902 4.313 ( -3.041 3.317 -0.087) 4.501 5.452 ( 12.290 10.842 -1.837) 16.491 7.915 ( 6.200 20.872 0.028) 21.774 8.344 ( -0.417 0.744 3.781) 3.876 8.830 ( 6.251 -5.060 -1.117) 8.119 ======================= Grid point 64 (29/60) ======================= q-point: ( 0.11 -0.11 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.890 ( 2.525 14.098 0.000) 14.322 1.030 ( 4.058 9.544 0.000) 10.371 1.310 ( 2.610 1.171 -0.000) 2.861 1.407 ( 0.232 0.811 0.000) 0.844 1.436 ( 1.335 8.505 0.000) 8.609 1.670 ( -0.457 -2.495 -0.000) 2.537 1.864 ( -1.530 -4.956 -0.000) 5.187 1.896 ( 9.998 -3.781 -0.000) 10.688 2.145 ( -8.195 -2.644 0.000) 8.611 2.244 ( -8.331 1.751 0.000) 8.513 2.563 ( 10.624 -0.757 -0.000) 10.650 2.599 ( -0.807 5.523 0.000) 5.582 3.222 ( -0.878 12.164 0.000) 12.196 3.292 ( 5.623 8.963 0.000) 10.581 4.001 ( -0.282 -4.083 -0.000) 4.093 4.063 ( -4.233 -7.120 -0.000) 8.283 4.135 ( 11.858 -4.521 -0.000) 12.690 5.369 ( -0.187 13.198 0.000) 13.199 7.499 ( 16.617 2.487 -0.000) 16.802 8.288 ( 0.586 -1.502 -0.000) 1.612 8.968 ( -5.845 7.277 0.000) 9.334 ======================= Grid point 65 (30/60) ======================= q-point: ( 0.11 0.22 0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.285 ( 3.486 3.633 0.912) 5.117 1.310 ( 4.376 2.399 1.942) 5.355 1.407 ( 0.347 1.212 -1.671) 2.094 1.422 ( -0.917 1.166 -0.137) 1.490 1.503 ( -0.699 -0.221 1.294) 1.487 1.588 ( -1.618 -0.175 -2.231) 2.761 1.717 ( 0.549 -3.002 -1.077) 3.236 1.764 ( -3.917 -5.376 -0.112) 6.653 1.847 ( -2.049 -8.191 0.046) 8.444 2.247 ( -1.494 6.364 0.819) 6.588 2.449 ( 1.559 -5.274 -1.225) 5.635 2.987 ( -0.068 8.471 1.284) 8.568 3.327 ( 7.604 -4.453 0.435) 8.822 3.536 ( -8.157 -6.470 -2.522) 10.713 3.737 ( 0.459 4.722 1.605) 5.008 3.935 ( -0.372 -1.687 -0.851) 1.926 4.250 ( -1.403 -3.215 -0.023) 3.508 5.936 ( 5.121 14.858 -0.270) 15.718 8.388 ( 1.098 2.103 0.768) 2.494 8.449 ( 2.546 18.143 -0.055) 18.321 8.827 ( 2.841 -4.479 0.110) 5.305 ======================= Grid point 67 (31/60) ======================= q-point: ( 0.11 -0.11 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.199 ( 2.702 9.111 0.000) 9.503 1.236 ( 2.625 3.929 0.000) 4.725 1.391 ( 2.215 0.722 -0.000) 2.330 1.396 ( 1.616 -2.362 -0.000) 2.861 1.501 ( -1.087 -0.623 0.000) 1.253 1.635 ( -0.802 -0.586 0.000) 0.993 1.767 ( -2.011 -1.036 0.000) 2.262 1.908 ( -7.714 -3.623 0.000) 8.522 2.036 ( -3.782 -0.051 0.000) 3.782 2.079 ( 2.983 -1.464 -0.000) 3.323 2.588 ( -5.599 0.976 0.000) 5.684 2.774 ( 10.191 0.081 -0.000) 10.192 3.409 ( -3.697 9.990 0.000) 10.652 3.551 ( 8.432 -0.976 -0.000) 8.488 3.752 ( -7.976 -9.567 -0.000) 12.456 3.962 ( -0.515 1.295 0.000) 1.394 4.258 ( 4.394 -1.581 -0.000) 4.670 5.642 ( 10.482 4.058 -0.000) 11.240 7.998 ( 23.366 1.798 -0.000) 23.435 8.333 ( 2.208 2.191 -0.000) 3.110 8.944 ( -6.420 2.752 0.000) 6.985 ======================= Grid point 68 (32/60) ======================= q-point: ( 0.11 -0.78 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.348 ( -0.900 0.900 -0.495) 1.365 1.359 ( 0.362 -0.362 -0.246) 0.568 1.427 ( -1.878 1.878 0.413) 2.688 1.445 ( -1.368 1.368 -0.186) 1.944 1.522 ( 0.250 -0.250 -0.198) 0.405 1.577 ( 0.013 -0.013 -0.943) 0.943 1.620 ( 0.581 -0.581 2.536) 2.666 1.694 ( 2.075 -2.075 -1.186) 3.165 1.720 ( 3.480 -3.480 0.038) 4.922 2.338 ( 0.498 -0.498 1.942) 2.065 2.363 ( -0.565 0.565 -1.733) 1.909 3.161 ( -3.348 3.348 -1.511) 4.970 3.263 ( -0.881 0.881 0.553) 1.363 3.373 ( 5.985 -5.985 0.472) 8.478 3.773 ( -2.442 2.442 5.461) 6.462 3.935 ( 1.110 -1.110 -4.383) 4.656 4.189 ( 1.701 -1.701 0.022) 2.406 6.166 ( -6.178 6.178 0.071) 8.737 8.428 ( -0.390 0.390 -0.345) 0.651 8.707 ( -5.448 5.448 0.360) 7.714 8.776 ( 0.131 -0.131 -0.331) 0.379 ======================= Grid point 70 (33/60) ======================= q-point: (-0.89 -0.11 0.44) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.332 ( 1.489 -0.318 -0.000) 1.523 1.341 ( 2.739 1.452 -0.000) 3.100 1.371 ( 2.879 -3.377 -0.000) 4.438 1.433 ( -0.882 -3.279 -0.000) 3.396 1.497 ( 1.585 4.808 0.000) 5.062 1.599 ( -0.800 -0.494 0.000) 0.940 1.675 ( -4.998 -2.580 0.000) 5.625 1.743 ( -2.853 1.016 0.000) 3.029 1.847 ( -8.349 2.055 0.000) 8.598 2.247 ( 6.508 2.264 -0.000) 6.890 2.441 ( -5.225 -3.284 0.000) 6.171 3.015 ( 9.553 3.509 -0.000) 10.177 3.366 ( -7.863 -9.449 -0.000) 12.293 3.493 ( -5.865 7.430 0.000) 9.466 3.622 ( 5.987 -5.080 -0.000) 7.852 4.001 ( -0.935 2.739 0.000) 2.894 4.249 ( -3.316 1.319 0.000) 3.569 5.928 ( 15.157 -5.927 -0.000) 16.275 8.422 ( 0.863 1.953 0.000) 2.136 8.446 ( 18.026 -2.704 -0.000) 18.227 8.822 ( -3.723 -3.308 0.000) 4.980 ======================= Grid point 77 (34/60) ======================= q-point: ( 0.11 0.56 -0.44) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 244 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.064 ( 2.882 -6.806 4.561) 8.685 1.103 ( 4.151 -7.525 3.050) 9.119 1.276 ( 1.248 -1.913 -1.178) 2.570 1.393 ( -0.061 -5.028 -0.699) 5.077 1.561 ( 0.536 -3.407 3.791) 5.125 1.592 ( -1.012 3.604 0.417) 3.767 1.854 ( -1.601 2.385 -1.181) 3.106 1.917 ( 8.809 3.568 -0.266) 9.508 2.060 ( -9.608 0.699 -3.020) 10.095 2.261 ( -7.817 -1.435 0.849) 7.993 2.506 ( 0.674 -4.612 0.845) 4.737 2.625 ( 8.657 -0.594 -1.676) 8.838 3.074 ( 8.994 -12.355 -3.999) 15.796 3.384 ( -2.733 -12.949 0.802) 13.258 3.906 ( -3.105 8.594 0.875) 9.179 4.117 ( 12.025 4.105 -0.611) 12.721 4.176 ( -1.518 6.791 0.601) 6.984 5.315 ( -0.907 -16.240 -1.656) 16.350 7.499 ( 16.562 -2.449 0.000) 16.742 8.466 ( -0.919 4.600 8.444) 9.659 8.872 ( -3.884 -9.251 -5.784) 11.581 ======================= Grid point 80 (35/60) ======================= q-point: ( 0.11 0.56 -0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.044 ( 4.018 -4.018 8.292) 10.052 1.059 ( 6.585 -6.585 3.646) 10.001 1.276 ( 2.253 -2.253 -3.052) 4.412 1.302 ( 4.558 -4.558 -1.132) 6.544 1.516 ( 2.448 -2.448 6.333) 7.217 1.633 ( -1.303 1.303 -1.510) 2.383 1.852 ( -4.888 4.888 -1.531) 7.080 1.954 ( -2.952 2.952 -6.095) 7.388 2.138 ( -3.040 3.040 1.150) 4.451 2.221 ( 2.223 -2.223 0.829) 3.252 2.251 ( 3.193 -3.193 -3.548) 5.743 2.875 ( 6.145 -6.145 4.620) 9.842 3.024 ( 5.054 -5.054 -0.243) 7.152 3.115 ( 11.734 -11.734 -6.359) 17.771 3.966 ( -6.443 6.443 0.202) 9.114 4.244 ( -5.521 5.521 2.006) 8.062 4.347 ( 1.054 -1.054 -0.083) 1.493 5.090 ( 14.563 -14.563 -3.627) 20.912 7.727 ( 12.242 -12.242 0.204) 17.314 8.416 ( 5.716 -5.716 5.228) 9.627 8.744 ( -5.285 5.285 -2.373) 7.841 ======================= Grid point 84 (36/60) ======================= q-point: ( 0.17 -0.17 -0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.865 ( -15.269 -0.000 0.000) 15.269 0.977 ( -11.635 -0.000 0.000) 11.635 1.279 ( -1.214 -0.000 0.000) 1.214 1.404 ( -0.625 -0.000 0.000) 0.625 1.420 ( -11.024 -0.000 0.000) 11.024 1.677 ( 3.039 0.000 -0.000) 3.039 1.773 ( 3.889 0.000 -0.000) 3.889 1.881 ( 5.308 0.000 -0.000) 5.308 2.316 ( 1.030 0.000 -0.000) 1.030 2.364 ( -2.654 -0.000 0.000) 2.654 2.365 ( 1.872 0.000 -0.000) 1.872 2.571 ( -6.693 -0.000 0.000) 6.693 3.231 ( -12.855 -0.000 0.000) 12.855 3.254 ( -6.764 -0.000 0.000) 6.764 3.950 ( 3.954 0.000 -0.000) 3.954 4.021 ( 4.256 0.000 -0.000) 4.256 4.108 ( 6.643 0.000 -0.000) 6.643 5.391 ( -15.018 -0.000 0.000) 15.018 7.308 ( -1.028 -0.000 0.000) 1.028 8.283 ( 2.244 0.000 -0.000) 2.244 9.038 ( -6.743 -0.000 0.000) 6.743 ======================= Grid point 85 (37/60) ======================= q-point: ( 0.17 0.17 -0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 45 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.205 ( -0.000 -0.000 9.979) 9.979 0.205 ( -0.000 -0.000 9.979) 9.979 0.602 ( -0.000 -0.000 30.644) 30.644 1.202 ( 0.000 0.000 -1.673) 1.673 1.202 ( 0.000 0.000 -1.673) 1.673 1.800 ( 0.000 0.000 -3.851) 3.851 1.931 ( -0.000 -0.000 0.914) 0.914 1.931 ( -0.000 -0.000 0.914) 0.914 2.275 ( 0.000 0.000 -0.032) 0.032 2.275 ( 0.000 0.000 -0.032) 0.032 2.349 ( 0.000 0.000 -1.154) 1.154 2.349 ( 0.000 0.000 -1.154) 1.154 2.535 ( 0.000 0.000 -0.340) 0.340 2.776 ( 0.000 0.000 -11.697) 11.697 4.119 ( 0.000 0.000 -0.068) 0.068 4.119 ( 0.000 0.000 -0.068) 0.068 4.294 ( 0.000 0.000 -2.298) 2.298 4.622 ( -0.000 -0.000 8.906) 8.906 7.256 ( -0.000 -0.000 0.005) 0.005 7.256 ( -0.000 -0.000 0.005) 0.005 9.751 ( 0.000 0.000 -2.540) 2.540 ======================= Grid point 86 (38/60) ======================= q-point: ( 0.17 0.50 -0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.022 ( 10.061 0.000 1.501) 10.172 1.057 ( 6.808 0.000 5.726) 8.896 1.263 ( 1.890 0.000 -0.602) 1.983 1.394 ( 5.695 0.000 -0.760) 5.746 1.553 ( 3.253 0.000 3.080) 4.480 1.608 ( -2.412 -0.000 1.680) 2.939 1.797 ( -4.389 -0.000 0.782) 4.459 1.883 ( -2.320 -0.000 -2.087) 3.121 2.274 ( 0.231 -0.000 0.895) 0.925 2.295 ( -0.435 -0.000 -0.677) 0.805 2.402 ( 1.504 0.000 -1.731) 2.293 2.583 ( 6.576 0.000 0.367) 6.586 2.966 ( 13.168 0.000 -6.330) 14.610 3.418 ( 12.884 0.000 1.330) 12.953 3.939 ( -8.804 -0.000 1.689) 8.965 3.952 ( -3.966 -0.000 0.039) 3.966 4.198 ( -6.807 -0.000 -0.438) 6.821 5.346 ( 17.611 0.000 -1.435) 17.669 7.309 ( 1.016 0.000 0.031) 1.017 8.477 ( -5.236 -0.000 9.434) 10.789 8.918 ( 9.061 0.000 -6.941) 11.414 ======================= Grid point 94 (39/60) ======================= q-point: ( 0.17 0.17 0.17) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.163 ( 6.740 6.740 3.873) 10.288 1.225 ( 3.938 3.938 1.365) 5.734 1.406 ( 0.511 0.511 -0.873) 1.133 1.410 ( 1.356 1.356 -1.469) 2.416 1.508 ( 0.296 0.296 1.566) 1.621 1.622 ( -0.953 -0.953 -1.719) 2.184 1.735 ( -2.532 -2.532 -2.479) 4.355 1.858 ( -4.867 -4.867 -1.264) 6.998 1.927 ( -5.385 -5.385 0.323) 7.623 2.277 ( 1.534 1.534 1.912) 2.893 2.405 ( 0.088 0.088 -3.383) 3.385 2.968 ( 3.696 3.696 3.605) 6.349 3.202 ( 2.830 2.830 -0.703) 4.063 3.628 ( 4.851 4.851 -2.193) 7.202 3.729 ( -5.337 -5.337 -3.781) 8.442 3.982 ( -3.518 -3.518 2.918) 5.768 4.299 ( -2.732 -2.732 -0.075) 3.865 5.733 ( 13.267 13.267 -0.983) 18.788 8.293 ( 13.298 13.298 0.038) 18.806 8.364 ( 1.147 1.147 2.429) 2.920 8.773 ( 0.325 0.325 -0.306) 0.552 ======================= Grid point 96 (40/60) ======================= q-point: ( 0.17 -0.17 0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.169 ( 0.283 11.678 0.000) 11.681 1.187 ( 1.847 4.821 0.000) 5.163 1.333 ( 2.771 1.256 -0.000) 3.042 1.371 ( 0.873 -4.086 -0.000) 4.179 1.526 ( -1.281 -0.175 0.000) 1.293 1.650 ( -0.594 0.865 0.000) 1.049 1.803 ( -1.246 -0.973 0.000) 1.581 1.839 ( 10.340 -1.184 -0.000) 10.407 2.091 ( -9.662 -2.381 0.000) 9.951 2.278 ( -9.975 1.261 0.000) 10.054 2.548 ( 11.166 -0.464 -0.000) 11.176 2.674 ( -2.065 1.401 0.000) 2.495 3.371 ( 7.483 -1.287 -0.000) 7.592 3.469 ( -1.909 11.051 0.000) 11.215 3.894 ( -5.048 -8.789 -0.000) 10.136 3.969 ( 0.441 1.317 0.000) 1.389 4.056 ( 11.626 -2.509 -0.000) 11.893 5.569 ( -2.347 4.885 0.000) 5.420 7.537 ( 18.370 0.978 -0.000) 18.396 8.294 ( 1.335 2.240 0.000) 2.608 9.071 ( -5.016 1.868 0.000) 5.353 ======================= Grid point 97 (41/60) ======================= q-point: ( 0.17 0.17 0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.338 ( 1.454 1.454 0.427) 2.100 1.351 ( 1.006 1.006 0.517) 1.514 1.440 ( 1.950 1.950 -1.016) 2.939 1.442 ( 0.386 0.386 0.089) 0.553 1.504 ( 0.366 0.366 1.576) 1.659 1.582 ( -0.583 -0.583 -1.324) 1.560 1.668 ( -3.788 -3.788 1.359) 5.526 1.674 ( -0.547 -0.547 -2.640) 2.751 1.720 ( -3.471 -3.471 0.018) 4.908 2.325 ( 0.668 0.668 1.922) 2.142 2.378 ( -0.768 -0.768 -1.858) 2.153 3.134 ( 4.598 4.598 1.128) 6.599 3.274 ( 0.747 0.747 -0.529) 1.181 3.389 ( -7.061 -7.061 -1.114) 10.048 3.795 ( 0.539 0.539 4.675) 4.737 3.919 ( 0.279 0.279 -4.096) 4.115 4.189 ( -1.741 -1.741 -0.023) 2.462 6.168 ( 6.016 6.016 -0.087) 8.508 8.422 ( 1.026 1.026 0.288) 1.480 8.714 ( 5.599 5.599 -0.360) 7.926 8.768 ( -0.295 -0.295 0.401) 0.579 ======================= Grid point 99 (42/60) ======================= q-point: (-0.83 -0.17 0.39) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.274 ( 3.215 -0.417 -0.000) 3.242 1.293 ( 2.119 -2.539 -0.000) 3.307 1.316 ( 1.938 -1.994 -0.000) 2.780 1.458 ( -1.235 -2.950 -0.000) 3.198 1.467 ( 1.388 5.939 0.000) 6.100 1.627 ( -1.755 -0.360 0.000) 1.792 1.760 ( -2.843 -0.153 0.000) 2.847 1.852 ( -7.676 -1.436 0.000) 7.809 2.051 ( -5.996 1.358 0.000) 6.147 2.076 ( 4.795 1.222 -0.000) 4.948 2.558 ( -5.360 -3.790 0.000) 6.565 2.798 ( 11.075 2.448 -0.000) 11.342 3.460 ( 9.140 -7.037 -0.000) 11.535 3.529 ( -6.326 -10.126 -0.000) 11.940 3.621 ( -6.702 8.570 0.000) 10.879 4.021 ( -0.980 4.030 0.000) 4.148 4.260 ( 4.435 1.671 -0.000) 4.740 5.623 ( 10.678 -5.952 -0.000) 12.225 7.996 ( 23.363 -1.864 -0.000) 23.437 8.402 ( 0.872 4.460 0.000) 4.544 8.921 ( -5.291 -5.213 0.000) 7.428 ======================= Grid point 103 (43/60) ======================= q-point: ( 0.17 0.83 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.286 ( 4.206 -4.206 -0.000) 5.948 1.309 ( 1.519 -1.519 -0.000) 2.149 1.322 ( 3.853 -3.853 -0.000) 5.449 1.385 ( 0.300 -0.300 -0.000) 0.425 1.576 ( -1.403 1.403 0.000) 1.985 1.599 ( -0.509 0.509 0.000) 0.719 1.667 ( -2.929 2.929 0.000) 4.142 1.777 ( -2.168 2.168 0.000) 3.065 1.931 ( -5.845 5.845 0.000) 8.266 2.316 ( 1.054 -1.054 -0.000) 1.491 2.334 ( 1.788 -1.788 -0.000) 2.529 3.091 ( 2.015 -2.015 -0.000) 2.850 3.179 ( 2.735 -2.735 -0.000) 3.868 3.479 ( 8.506 -8.506 -0.000) 12.029 3.657 ( -7.682 7.682 0.000) 10.864 4.068 ( -3.365 3.365 0.000) 4.759 4.297 ( -2.736 2.736 0.000) 3.869 5.704 ( 14.466 -14.466 -0.000) 20.458 8.294 ( 13.030 -13.030 -0.000) 18.427 8.466 ( -2.380 2.380 0.000) 3.367 8.739 ( 2.661 -2.661 -0.000) 3.763 ======================= Grid point 110 (44/60) ======================= q-point: ( 0.17 0.50 -0.39) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 243 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.899 ( 8.400 -9.177 3.330) 12.879 0.977 ( 2.588 -3.864 9.921) 10.957 1.212 ( 2.369 -8.223 -0.697) 8.585 1.242 ( 1.774 -4.646 -2.389) 5.517 1.467 ( 1.939 -4.642 7.956) 9.413 1.655 ( -0.908 1.979 -1.979) 2.943 1.933 ( -1.786 5.298 -2.538) 6.140 1.989 ( 7.873 2.893 -0.737) 8.420 2.084 ( -7.541 1.485 -3.895) 8.616 2.228 ( -3.812 -1.902 0.312) 4.272 2.345 ( -5.241 -9.076 -1.318) 10.563 2.638 ( 13.214 -0.458 1.288) 13.284 2.805 ( 11.759 -11.142 -4.984) 16.948 3.082 ( -2.457 -15.560 0.019) 15.753 4.073 ( -3.547 6.973 0.705) 7.855 4.201 ( 11.701 3.363 -0.517) 12.186 4.333 ( -2.806 7.800 2.433) 8.639 4.914 ( 1.833 -20.496 -3.988) 20.960 7.439 ( 12.955 -3.181 0.063) 13.340 8.312 ( 3.168 -19.332 1.797) 19.672 8.905 ( -7.725 12.880 0.336) 15.023 ======================= Grid point 114 (45/60) ======================= q-point: ( 0.22 -0.22 -0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.166 ( -5.602 -0.000 0.000) 5.602 1.171 ( -13.483 -0.000 0.000) 13.483 1.297 ( -0.349 -0.000 0.000) 0.349 1.362 ( 4.477 0.000 -0.000) 4.477 1.546 ( 0.200 0.000 -0.000) 0.200 1.659 ( -2.432 -0.000 0.000) 2.432 1.704 ( 1.968 0.000 -0.000) 1.968 1.816 ( 0.938 0.000 -0.000) 0.938 2.297 ( 0.711 0.000 -0.000) 0.711 2.326 ( 1.745 0.000 -0.000) 1.745 2.407 ( -1.230 -0.000 0.000) 1.230 2.682 ( -3.230 -0.000 0.000) 3.230 3.286 ( 3.150 0.000 -0.000) 3.150 3.489 ( -11.304 -0.000 0.000) 11.304 3.885 ( 1.829 0.000 -0.000) 1.829 3.949 ( 8.338 0.000 -0.000) 8.338 3.985 ( -0.814 -0.000 0.000) 0.814 5.620 ( -5.717 -0.000 0.000) 5.717 7.324 ( -0.431 -0.000 0.000) 0.431 8.280 ( -2.244 -0.000 0.000) 2.244 9.128 ( -1.285 -0.000 0.000) 1.285 ======================= Grid point 116 (46/60) ======================= q-point: ( 0.22 0.44 -0.22) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.779 ( 12.672 0.000 2.709) 12.959 0.948 ( 4.125 0.000 10.457) 11.242 1.194 ( 12.424 0.000 -0.813) 12.451 1.219 ( 2.173 0.000 -1.178) 2.472 1.445 ( 5.631 0.000 8.413) 10.124 1.667 ( -2.041 -0.000 -2.276) 3.057 1.881 ( -3.358 -0.000 1.026) 3.511 1.954 ( -4.593 -0.000 -4.119) 6.170 2.239 ( 3.719 0.000 1.026) 3.858 2.305 ( -0.497 -0.000 -1.046) 1.158 2.389 ( -0.315 0.000 -2.057) 2.081 2.437 ( 6.644 0.000 0.307) 6.651 2.657 ( 14.971 0.000 -7.664) 16.819 3.114 ( 15.775 0.000 0.014) 15.775 4.035 ( -3.633 -0.000 0.020) 3.633 4.111 ( -7.183 -0.000 0.909) 7.241 4.363 ( -8.542 -0.000 2.418) 8.878 4.916 ( 21.854 0.000 -3.283) 22.100 7.285 ( 1.152 0.000 0.029) 1.152 8.281 ( 23.990 0.000 1.101) 24.015 8.998 ( -16.689 -0.000 0.346) 16.693 ======================= Grid point 129 (47/60) ======================= q-point: (-0.78 -0.22 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.214 ( 2.381 -2.336 -0.000) 3.336 1.252 ( 1.521 -5.768 -0.000) 5.965 1.293 ( 0.413 -0.347 -0.000) 0.539 1.437 ( 1.263 8.089 0.000) 8.187 1.481 ( -0.887 -3.338 -0.000) 3.454 1.666 ( -1.746 0.081 0.000) 1.748 1.805 ( -1.502 0.747 0.000) 1.678 1.851 ( 9.898 2.038 -0.000) 10.106 2.048 ( -10.715 -1.628 0.000) 10.838 2.290 ( -9.935 -0.235 0.000) 9.938 2.552 ( 12.207 1.021 -0.000) 12.249 2.651 ( -3.294 -3.547 -0.000) 4.840 3.261 ( 8.944 -8.580 -0.000) 12.394 3.616 ( -2.389 -4.188 -0.000) 4.822 3.759 ( -5.396 3.077 0.000) 6.211 4.038 ( 1.407 3.223 0.000) 3.517 4.050 ( 10.728 3.505 -0.000) 11.286 5.552 ( -2.575 -6.633 -0.000) 7.115 7.537 ( 18.300 -1.022 -0.000) 18.329 8.388 ( 0.407 6.824 0.000) 6.836 9.025 ( -4.009 -6.707 -0.000) 7.814 ======================= Grid point 133 (48/60) ======================= q-point: (-0.78 -0.22 0.44) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.200 ( 3.712 -4.225 -0.000) 5.624 1.222 ( 5.273 -4.268 -0.000) 6.784 1.281 ( 1.108 -2.215 -0.000) 2.476 1.385 ( 0.039 -4.051 -0.000) 4.051 1.578 ( 0.406 3.614 0.000) 3.637 1.619 ( -1.355 -0.260 0.000) 1.380 1.765 ( -4.658 1.344 0.000) 4.848 1.846 ( -5.292 0.514 0.000) 5.317 2.086 ( -8.531 2.202 0.000) 8.811 2.132 ( 10.341 3.775 -0.000) 11.008 2.431 ( -5.049 -8.234 -0.000) 9.659 2.877 ( 12.705 4.762 -0.000) 13.568 3.269 ( 11.235 -10.942 -0.000) 15.683 3.282 ( -6.240 -12.728 -0.000) 14.176 3.809 ( -6.604 8.547 0.000) 10.801 4.129 ( -2.308 6.123 0.000) 6.543 4.313 ( 3.288 2.818 -0.000) 4.331 5.398 ( 11.661 -14.865 -0.000) 18.893 7.917 ( 21.031 -5.925 -0.000) 21.850 8.524 ( -3.082 7.235 0.000) 7.864 8.734 ( -1.616 -11.753 -0.000) 11.864 ======================= Grid point 144 (49/60) ======================= q-point: ( 0.22 0.44 -0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 140 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.699 ( 10.297 -10.297 4.282) 15.178 0.843 ( 10.842 -10.842 6.582) 16.686 1.048 ( 1.971 -1.971 6.556) 7.124 1.195 ( 1.455 -1.455 -2.276) 3.069 1.369 ( 4.489 -4.489 9.097) 11.092 1.686 ( -1.150 1.150 -3.196) 3.586 2.031 ( 1.958 -1.958 1.471) 3.135 2.032 ( 1.148 -1.148 -3.243) 3.626 2.107 ( -4.627 4.627 -1.781) 6.782 2.159 ( -1.562 1.562 -6.671) 7.027 2.262 ( -1.141 1.141 0.409) 1.665 2.490 ( 9.631 -9.631 -0.187) 13.622 2.723 ( 7.098 -7.098 -3.984) 10.799 2.766 ( 6.969 -6.969 0.033) 9.855 4.187 ( -3.764 3.764 -0.036) 5.323 4.219 ( 4.005 -4.005 -0.135) 5.665 4.413 ( 4.925 -4.925 -0.380) 6.976 4.611 ( 0.568 -0.568 -0.482) 0.936 7.356 ( 5.208 -5.208 0.104) 7.365 7.863 ( 18.981 -18.981 -3.569) 27.079 9.238 ( -15.746 15.746 3.445) 22.533 ======================= Grid point 145 (50/60) ======================= q-point: ( 0.28 0.72 -0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.185 ( 3.735 0.000 -0.000) 3.735 1.235 ( 6.363 0.000 -0.000) 6.363 1.291 ( 0.990 0.000 -0.000) 0.990 1.422 ( -9.548 -0.000 0.000) 9.548 1.493 ( 3.755 0.000 -0.000) 3.755 1.697 ( -0.053 -0.000 0.000) 0.053 1.717 ( -3.108 -0.000 0.000) 3.108 1.821 ( -0.931 -0.000 0.000) 0.931 2.287 ( 0.198 0.000 -0.000) 0.198 2.298 ( 0.761 0.000 -0.000) 0.761 2.411 ( 0.891 0.000 -0.000) 0.891 2.689 ( 2.606 0.000 -0.000) 2.606 3.148 ( 9.475 0.000 -0.000) 9.475 3.640 ( 2.508 0.000 -0.000) 2.508 3.819 ( -1.543 -0.000 0.000) 1.543 3.886 ( -1.891 -0.000 0.000) 1.891 4.047 ( -4.867 -0.000 0.000) 4.867 5.606 ( 7.169 0.000 -0.000) 7.169 7.324 ( 0.380 0.000 -0.000) 0.380 8.384 ( -7.839 -0.000 0.000) 7.839 9.070 ( 7.306 0.000 -0.000) 7.306 ======================= Grid point 148 (51/60) ======================= q-point: ( 0.28 0.39 -0.28) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 135 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.513 ( 11.660 0.000 4.484) 12.493 0.715 ( 20.341 0.000 5.295) 21.019 1.051 ( -1.488 -0.000 9.893) 10.004 1.179 ( 1.441 0.000 -1.731) 2.253 1.311 ( 7.249 0.000 9.301) 11.792 1.701 ( -1.292 0.000 -3.681) 3.901 1.934 ( -1.670 -0.000 1.062) 1.979 2.044 ( 2.373 0.000 -1.563) 2.841 2.152 ( -2.810 -0.000 -3.857) 4.772 2.308 ( 1.317 0.000 -0.518) 1.415 2.329 ( 2.320 0.000 -0.629) 2.403 2.337 ( 5.448 0.000 -0.106) 5.449 2.486 ( -0.378 0.000 -7.626) 7.636 2.774 ( 16.025 0.000 -2.022) 16.152 4.096 ( -2.026 -0.000 -0.037) 2.026 4.227 ( -3.727 -0.000 -0.346) 3.743 4.361 ( 16.358 0.000 -2.338) 16.524 4.662 ( -6.152 -0.000 4.241) 7.472 7.264 ( 0.747 0.000 0.014) 0.747 7.659 ( 31.005 0.000 -3.443) 31.195 9.432 ( -21.081 -0.000 1.953) 21.171 ======================= Grid point 165 (52/60) ======================= q-point: (-0.72 -0.28 0.39) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.107 ( 5.316 -7.155 -0.000) 8.914 1.138 ( 2.071 -4.329 -0.000) 4.799 1.263 ( 0.604 -2.211 -0.000) 2.292 1.384 ( 0.045 -6.362 -0.000) 6.362 1.568 ( 0.455 3.591 0.000) 3.620 1.650 ( -1.391 -1.280 0.000) 1.890 1.828 ( -1.731 1.576 0.000) 2.341 1.915 ( 9.155 3.766 -0.000) 9.900 2.026 ( -11.231 -0.435 0.000) 11.240 2.271 ( -7.995 -1.421 0.000) 8.120 2.529 ( -4.211 -7.889 -0.000) 8.942 2.593 ( 14.106 2.746 -0.000) 14.370 3.031 ( 10.457 -12.877 -0.000) 16.588 3.389 ( -3.615 -13.433 -0.000) 13.911 3.923 ( -4.000 8.476 0.000) 9.373 4.108 ( 10.736 4.216 -0.000) 11.534 4.179 ( 0.537 7.450 0.000) 7.469 5.298 ( -1.061 -17.045 -0.000) 17.077 7.499 ( 16.559 -2.426 -0.000) 16.735 8.587 ( -2.576 12.232 0.000) 12.500 8.780 ( -2.127 -16.642 -0.000) 16.778 ======================= Grid point 169 (53/60) ======================= q-point: ( 0.28 0.72 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.106 ( 6.726 -6.726 -0.000) 9.513 1.128 ( 2.682 -2.682 -0.000) 3.793 1.236 ( 1.798 -1.798 -0.000) 2.543 1.288 ( 4.987 -4.987 -0.000) 7.052 1.607 ( 0.364 -0.364 -0.000) 0.515 1.619 ( -0.283 0.283 0.000) 0.400 1.837 ( -5.017 5.017 0.000) 7.095 1.869 ( -2.336 2.336 0.000) 3.304 2.149 ( -3.521 3.521 0.000) 4.980 2.216 ( 3.419 -3.419 -0.000) 4.835 2.230 ( 2.779 -2.779 -0.000) 3.930 2.956 ( 4.850 -4.850 -0.000) 6.860 3.012 ( 13.121 -13.121 -0.000) 18.556 3.021 ( 4.888 -4.888 -0.000) 6.913 3.971 ( -6.712 6.712 0.000) 9.492 4.262 ( -5.800 5.800 0.000) 8.203 4.346 ( 1.071 -1.071 -0.000) 1.514 5.053 ( 15.456 -15.456 -0.000) 21.859 7.729 ( 12.245 -12.245 -0.000) 17.317 8.485 ( 10.281 -10.281 -0.000) 14.540 8.704 ( -9.649 9.649 0.000) 13.645 ======================= Grid point 177 (54/60) ======================= q-point: ( 0.33 0.67 -0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.038 ( 9.504 0.000 -0.000) 9.504 1.116 ( 4.392 0.000 -0.000) 4.392 1.257 ( 2.161 0.000 -0.000) 2.161 1.383 ( 7.298 0.000 -0.000) 7.298 1.562 ( -3.460 -0.000 0.000) 3.460 1.669 ( 2.050 0.000 -0.000) 2.050 1.805 ( -4.581 -0.000 0.000) 4.581 1.846 ( -1.465 -0.000 0.000) 1.465 2.288 ( -0.217 -0.000 0.000) 0.217 2.301 ( -1.225 -0.000 0.000) 1.225 2.368 ( 3.256 0.000 -0.000) 3.256 2.586 ( 6.543 0.000 -0.000) 6.543 2.898 ( 14.130 0.000 -0.000) 14.130 3.428 ( 13.521 0.000 -0.000) 13.521 3.952 ( -3.970 -0.000 0.000) 3.970 3.967 ( -8.248 -0.000 0.000) 8.248 4.185 ( -8.229 -0.000 0.000) 8.229 5.331 ( 18.297 0.000 -0.000) 18.297 7.309 ( 1.012 0.000 -0.000) 1.012 8.617 ( -14.481 -0.000 0.000) 14.481 8.805 ( 18.220 0.000 -0.000) 18.220 ======================= Grid point 180 (55/60) ======================= q-point: ( 0.33 0.33 -0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 60 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.365 ( -0.000 -0.000 6.351) 6.351 0.365 ( -0.000 -0.000 6.351) 6.351 1.143 ( -0.000 -0.000 24.779) 24.779 1.163 ( 0.000 0.000 -1.967) 1.967 1.163 ( 0.000 0.000 -1.967) 1.967 1.717 ( 0.000 0.000 -4.095) 4.095 1.951 ( -0.000 -0.000 1.031) 1.031 1.951 ( -0.000 -0.000 1.031) 1.031 2.276 ( -0.000 -0.000 0.201) 0.201 2.276 ( -0.000 -0.000 0.201) 0.201 2.322 ( 0.000 0.000 -1.502) 1.502 2.322 ( 0.000 0.000 -1.502) 1.502 2.495 ( 0.000 0.000 -15.924) 15.924 2.528 ( 0.000 0.000 -0.341) 0.341 4.117 ( 0.000 0.000 -0.067) 0.067 4.117 ( 0.000 0.000 -0.067) 0.067 4.266 ( 0.000 0.000 -0.738) 0.738 4.782 ( -0.000 -0.000 6.755) 6.755 7.256 ( -0.000 -0.000 0.005) 0.005 7.256 ( -0.000 -0.000 0.005) 0.005 9.697 ( 0.000 0.000 -2.594) 2.594 ======================= Grid point 202 (56/60) ======================= q-point: (-0.67 -0.33 0.44) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 164 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.935 ( 9.130 -8.748 -0.000) 12.645 1.084 ( 1.481 -0.757 -0.000) 1.663 1.194 ( 3.368 -11.891 -0.000) 12.359 1.211 ( 0.607 -2.483 -0.000) 2.556 1.585 ( 1.261 -1.546 -0.000) 1.995 1.630 ( -0.664 -0.315 0.000) 0.735 1.885 ( -1.044 4.315 0.000) 4.440 1.995 ( 9.037 3.607 -0.000) 9.731 2.031 ( -11.373 0.904 0.000) 11.409 2.232 ( -3.356 -3.036 0.000) 4.525 2.334 ( -5.987 -9.380 -0.000) 11.128 2.658 ( 16.073 3.099 -0.000) 16.368 2.743 ( 11.737 -13.715 -0.000) 18.052 3.080 ( -2.950 -15.851 -0.000) 16.123 4.083 ( -3.719 6.464 0.000) 7.457 4.195 ( 11.707 3.530 -0.000) 12.228 4.357 ( -2.840 9.029 0.000) 9.465 4.875 ( 1.327 -21.671 -0.000) 21.711 7.440 ( 13.021 -3.143 -0.000) 13.396 8.338 ( 3.946 -24.519 -0.000) 24.834 8.901 ( -8.012 17.243 0.000) 19.013 ======================= Grid point 210 (57/60) ======================= q-point: ( 0.39 0.61 -0.39) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.808 ( 11.984 0.000 -0.000) 11.984 1.068 ( 0.032 0.000 -0.000) 0.032 1.159 ( 14.498 0.000 -0.000) 14.498 1.206 ( 2.491 0.000 -0.000) 2.491 1.570 ( 2.363 0.000 -0.000) 2.363 1.639 ( 0.422 0.000 -0.000) 0.422 1.892 ( -3.429 -0.000 0.000) 3.429 1.896 ( -3.941 -0.000 0.000) 3.941 2.270 ( 6.351 0.000 -0.000) 6.351 2.294 ( -0.352 -0.000 0.000) 0.352 2.345 ( -2.585 -0.000 0.000) 2.585 2.441 ( 6.740 0.000 -0.000) 6.740 2.575 ( 15.335 0.000 -0.000) 15.335 3.112 ( 16.330 0.000 -0.000) 16.330 4.035 ( -3.610 -0.000 0.000) 3.610 4.123 ( -6.299 -0.000 0.000) 6.299 4.386 ( -10.166 -0.000 0.000) 10.166 4.883 ( 22.461 0.000 -0.000) 22.461 7.285 ( 1.157 0.000 -0.000) 1.157 8.300 ( 29.456 0.000 -0.000) 29.456 8.994 ( -20.978 -0.000 0.000) 20.978 ======================= Grid point 240 (58/60) ======================= q-point: ( 0.39 0.61 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 90 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.746 ( 10.113 -10.113 -0.000) 14.302 0.903 ( 13.337 -13.337 -0.000) 18.862 1.093 ( -0.937 0.937 0.000) 1.326 1.169 ( 1.240 -1.240 -0.000) 1.753 1.528 ( 3.086 -3.086 -0.000) 4.365 1.638 ( -0.827 0.827 0.000) 1.170 2.009 ( -4.834 4.834 0.000) 6.836 2.056 ( 3.785 -3.785 -0.000) 5.352 2.061 ( 3.633 -3.633 -0.000) 5.137 2.079 ( -6.399 6.399 0.000) 9.050 2.266 ( -1.026 1.026 0.000) 1.451 2.489 ( 9.874 -9.874 -0.000) 13.964 2.674 ( 8.332 -8.332 -0.000) 11.783 2.767 ( 7.032 -7.032 -0.000) 9.944 4.188 ( -3.452 3.452 0.000) 4.882 4.217 ( 3.995 -3.995 -0.000) 5.650 4.399 ( 8.996 -8.996 -0.000) 12.722 4.620 ( -5.050 5.050 0.000) 7.142 7.357 ( 5.245 -5.245 -0.000) 7.417 7.836 ( 19.599 -19.599 -0.000) 27.717 9.264 ( -15.221 15.221 0.000) 21.525 ======================= Grid point 244 (59/60) ======================= q-point: ( 0.44 0.56 -0.44) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 100 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.560 ( 10.620 0.000 -0.000) 10.620 0.770 ( 21.798 0.000 -0.000) 21.798 1.103 ( -2.687 -0.000 0.000) 2.687 1.160 ( 1.670 0.000 -0.000) 1.670 1.489 ( 4.407 0.000 -0.000) 4.407 1.650 ( -1.218 -0.000 0.000) 1.218 1.945 ( -1.656 -0.000 0.000) 1.656 2.040 ( -10.090 -0.000 0.000) 10.090 2.084 ( 10.006 0.000 -0.000) 10.006 2.301 ( -0.237 -0.000 0.000) 0.237 2.324 ( 4.028 0.000 -0.000) 4.028 2.340 ( 6.198 0.000 -0.000) 6.198 2.390 ( -1.061 -0.000 0.000) 1.061 2.754 ( 16.923 0.000 -0.000) 16.923 4.095 ( -1.998 -0.000 0.000) 1.998 4.225 ( -3.335 -0.000 0.000) 3.335 4.338 ( 17.271 0.000 -0.000) 17.271 4.709 ( -9.499 -0.000 0.000) 9.499 7.264 ( 0.754 0.000 -0.000) 0.754 7.633 ( 30.396 0.000 -0.000) 30.396 9.442 ( -19.183 -0.000 0.000) 19.183 ======================= Grid point 279 (60/60) ======================= q-point: (-0.50 -0.50 0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.64e-04 1.64e-04 1.64e-04 1.64e-04 0.00e+00 0.00e+00 3.36e-05 Number of triplets: 30 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.429 ( -0.000 -0.000 0.000) 0.000 0.429 ( -0.000 -0.000 0.000) 0.000 1.142 ( -0.000 -0.000 0.000) 0.000 1.142 ( -0.000 -0.000 0.000) 0.000 1.431 ( -0.000 -0.000 0.000) 0.000 1.669 ( 0.000 0.000 -0.000) 0.000 1.962 ( -0.000 -0.000 0.000) 0.000 1.962 ( -0.000 -0.000 0.000) 0.000 2.280 ( -0.000 -0.000 0.000) 0.000 2.280 ( -0.000 -0.000 0.000) 0.000 2.295 ( -0.000 -0.000 0.000) 0.000 2.303 ( -0.000 -0.000 0.000) 0.000 2.303 ( -0.000 -0.000 0.000) 0.000 2.524 ( -0.000 -0.000 0.000) 0.000 4.116 ( -0.000 -0.000 0.000) 0.000 4.116 ( -0.000 -0.000 0.000) 0.000 4.260 ( -0.000 -0.000 0.000) 0.000 4.850 ( -0.000 -0.000 0.000) 0.000 7.256 ( -0.000 -0.000 0.000) 0.000 7.256 ( -0.000 -0.000 0.000) 0.000 9.670 ( -0.000 -0.000 0.000) 0.000 =================== End of collection of collisions =================== ----------- Thermal conductivity (W/m-k) with tetrahedron method ----------- # T(K) xx yy zz yz xz xy #ipm 0.0 0.000 0.000 0.000 0.000 0.000 0.000 3/10206 10.0 86.617 86.617 54.051 0.000 -0.000 0.000 3/10206 20.0 13.944 13.944 8.943 0.000 -0.000 0.000 3/10206 30.0 6.312 6.312 3.954 0.000 -0.000 0.000 3/10206 40.0 4.109 4.109 2.465 0.000 -0.000 0.000 3/10206 50.0 3.127 3.127 1.791 0.000 -0.000 0.000 3/10206 60.0 2.570 2.570 1.412 0.000 -0.000 0.000 3/10206 70.0 2.204 2.204 1.169 0.000 -0.000 0.000 3/10206 80.0 1.941 1.941 0.999 0.000 -0.000 0.000 3/10206 90.0 1.739 1.739 0.874 0.000 -0.000 0.000 3/10206 100.0 1.578 1.578 0.777 0.000 -0.000 0.000 3/10206 110.0 1.446 1.446 0.700 0.000 -0.000 0.000 3/10206 120.0 1.335 1.335 0.638 0.000 -0.000 0.000 3/10206 130.0 1.240 1.240 0.586 0.000 -0.000 0.000 3/10206 140.0 1.158 1.158 0.542 0.000 -0.000 0.000 3/10206 150.0 1.087 1.087 0.504 0.000 -0.000 0.000 3/10206 160.0 1.023 1.023 0.471 0.000 -0.000 0.000 3/10206 170.0 0.967 0.967 0.442 0.000 -0.000 0.000 3/10206 180.0 0.917 0.917 0.417 0.000 -0.000 0.000 3/10206 190.0 0.871 0.871 0.395 0.000 -0.000 0.000 3/10206 200.0 0.830 0.830 0.374 0.000 -0.000 0.000 3/10206 210.0 0.792 0.792 0.356 0.000 -0.000 0.000 3/10206 220.0 0.758 0.758 0.340 0.000 -0.000 0.000 3/10206 230.0 0.726 0.726 0.325 0.000 -0.000 0.000 3/10206 240.0 0.697 0.697 0.311 0.000 -0.000 0.000 3/10206 250.0 0.671 0.671 0.298 0.000 -0.000 0.000 3/10206 260.0 0.646 0.646 0.287 0.000 -0.000 0.000 3/10206 270.0 0.623 0.623 0.276 0.000 -0.000 0.000 3/10206 280.0 0.601 0.601 0.266 0.000 -0.000 0.000 3/10206 290.0 0.581 0.581 0.257 0.000 -0.000 0.000 3/10206 300.0 0.562 0.562 0.248 0.000 -0.000 0.000 3/10206 310.0 0.545 0.545 0.240 0.000 -0.000 0.000 3/10206 320.0 0.528 0.528 0.233 0.000 -0.000 0.000 3/10206 330.0 0.513 0.513 0.225 0.000 -0.000 0.000 3/10206 340.0 0.498 0.498 0.219 0.000 -0.000 0.000 3/10206 350.0 0.484 0.484 0.212 0.000 -0.000 0.000 3/10206 360.0 0.471 0.471 0.207 0.000 -0.000 0.000 3/10206 370.0 0.458 0.458 0.201 0.000 -0.000 0.000 3/10206 380.0 0.447 0.447 0.196 0.000 -0.000 0.000 3/10206 390.0 0.435 0.435 0.191 0.000 -0.000 0.000 3/10206 400.0 0.425 0.425 0.186 0.000 -0.000 0.000 3/10206 410.0 0.414 0.414 0.181 0.000 -0.000 0.000 3/10206 420.0 0.405 0.405 0.177 0.000 -0.000 0.000 3/10206 430.0 0.396 0.396 0.173 0.000 -0.000 0.000 3/10206 440.0 0.387 0.387 0.169 0.000 -0.000 0.000 3/10206 450.0 0.378 0.378 0.165 0.000 -0.000 0.000 3/10206 460.0 0.370 0.370 0.161 0.000 -0.000 0.000 3/10206 470.0 0.362 0.362 0.158 0.000 -0.000 0.000 3/10206 480.0 0.355 0.355 0.155 0.000 -0.000 0.000 3/10206 490.0 0.348 0.348 0.152 0.000 -0.000 0.000 3/10206 500.0 0.341 0.341 0.149 0.000 -0.000 0.000 3/10206 510.0 0.334 0.334 0.146 0.000 -0.000 0.000 3/10206 520.0 0.328 0.328 0.143 0.000 -0.000 0.000 3/10206 530.0 0.322 0.322 0.140 0.000 -0.000 0.000 3/10206 540.0 0.316 0.316 0.137 0.000 -0.000 0.000 3/10206 550.0 0.310 0.310 0.135 0.000 -0.000 0.000 3/10206 560.0 0.305 0.305 0.133 0.000 -0.000 0.000 3/10206 570.0 0.299 0.299 0.130 0.000 -0.000 0.000 3/10206 580.0 0.294 0.294 0.128 0.000 -0.000 0.000 3/10206 590.0 0.289 0.289 0.126 0.000 -0.000 0.000 3/10206 600.0 0.285 0.285 0.124 0.000 -0.000 0.000 3/10206 610.0 0.280 0.280 0.122 0.000 -0.000 0.000 3/10206 620.0 0.276 0.276 0.120 0.000 -0.000 0.000 3/10206 630.0 0.271 0.271 0.118 0.000 -0.000 0.000 3/10206 640.0 0.267 0.267 0.116 0.000 -0.000 0.000 3/10206 650.0 0.263 0.263 0.114 0.000 -0.000 0.000 3/10206 660.0 0.259 0.259 0.112 0.000 -0.000 0.000 3/10206 670.0 0.255 0.255 0.111 0.000 -0.000 0.000 3/10206 680.0 0.251 0.251 0.109 0.000 -0.000 0.000 3/10206 690.0 0.248 0.248 0.108 0.000 -0.000 0.000 3/10206 700.0 0.244 0.244 0.106 0.000 -0.000 0.000 3/10206 710.0 0.241 0.241 0.105 0.000 -0.000 0.000 3/10206 720.0 0.238 0.238 0.103 0.000 -0.000 0.000 3/10206 730.0 0.234 0.234 0.102 0.000 -0.000 0.000 3/10206 740.0 0.231 0.231 0.100 0.000 -0.000 0.000 3/10206 750.0 0.228 0.228 0.099 0.000 -0.000 0.000 3/10206 760.0 0.225 0.225 0.098 0.000 -0.000 0.000 3/10206 770.0 0.222 0.222 0.096 0.000 -0.000 0.000 3/10206 780.0 0.219 0.219 0.095 0.000 -0.000 0.000 3/10206 790.0 0.217 0.217 0.094 0.000 -0.000 0.000 3/10206 800.0 0.214 0.214 0.093 0.000 -0.000 0.000 3/10206 810.0 0.211 0.211 0.092 0.000 -0.000 0.000 3/10206 820.0 0.209 0.209 0.090 0.000 -0.000 0.000 3/10206 830.0 0.206 0.206 0.089 0.000 -0.000 0.000 3/10206 840.0 0.204 0.204 0.088 0.000 -0.000 0.000 3/10206 850.0 0.201 0.201 0.087 0.000 -0.000 0.000 3/10206 860.0 0.199 0.199 0.086 0.000 -0.000 0.000 3/10206 870.0 0.197 0.197 0.085 0.000 -0.000 0.000 3/10206 880.0 0.195 0.195 0.084 0.000 -0.000 0.000 3/10206 890.0 0.192 0.192 0.083 0.000 -0.000 0.000 3/10206 900.0 0.190 0.190 0.082 0.000 -0.000 0.000 3/10206 910.0 0.188 0.188 0.082 0.000 -0.000 0.000 3/10206 920.0 0.186 0.186 0.081 0.000 -0.000 0.000 3/10206 930.0 0.184 0.184 0.080 0.000 -0.000 0.000 3/10206 940.0 0.182 0.182 0.079 0.000 -0.000 0.000 3/10206 950.0 0.180 0.180 0.078 0.000 -0.000 0.000 3/10206 960.0 0.178 0.178 0.077 0.000 -0.000 0.000 3/10206 970.0 0.177 0.177 0.076 0.000 -0.000 0.000 3/10206 980.0 0.175 0.175 0.076 0.000 -0.000 0.000 3/10206 990.0 0.173 0.173 0.075 0.000 -0.000 0.000 3/10206 1000.0 0.171 0.171 0.074 0.000 -0.000 0.000 3/10206 Thermal conductivity related properties were written into "kappa-m3918.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-09 00:53:59]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|