----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-09 04:49:19]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) Number of symmetry operations in supercell: 384 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.031293939063714 0.000000000000000 0.000000000000000 b -2.015646969531857 3.491202961351413 0.000000000000000 c 0.000000000000000 0.000000000000000 15.830686910000001 Atomic positions (fractional): *1 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 2 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 *3 Dy 0.00000000000000 0.00000000000000 0.50000000000000 162.500 4 Dy 0.00000000000000 0.00000000000000 0.00000000000000 162.500 *5 S 0.33333333333333 0.66666666666667 0.90806073924251 32.065 6 S 0.66666666666667 0.33333333333333 0.40806073924251 32.065 7 S 0.66666666666667 0.33333333333333 0.09193926075749 32.065 8 S 0.33333333333333 0.66666666666667 0.59193926075749 32.065 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.031293939063714 0.000000000000000 0.000000000000000 b -2.015646969531857 3.491202961351413 0.000000000000000 c 0.000000000000000 0.000000000000000 15.830686910000001 Atomic positions (fractional): *1 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 > 1 2 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 > 2 *3 Dy 0.00000000000000 0.00000000000000 0.50000000000000 162.500 > 3 4 Dy 0.00000000000000 0.00000000000000 0.00000000000000 162.500 > 4 *5 S 0.33333333333333 0.66666666666667 0.90806073924251 32.065 > 5 6 S 0.66666666666667 0.33333333333333 0.40806073924251 32.065 > 6 7 S 0.66666666666667 0.33333333333333 0.09193926075749 32.065 > 7 8 S 0.33333333333333 0.66666666666667 0.59193926075749 32.065 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.125175756254855 0.000000000000000 0.000000000000000 b -8.062587878127427 13.964811845405652 0.000000000000000 c 0.000000000000000 0.000000000000000 15.830686910000001 Atomic positions (fractional): *1 Cs 0.16666666666667 0.08333333333333 0.75000000000000 132.905 > 1 2 Cs 0.41666666666667 0.08333333333333 0.75000000000000 132.905 > 1 3 Cs 0.66666666666667 0.08333333333333 0.75000000000000 132.905 > 1 4 Cs 0.91666666666667 0.08333333333333 0.75000000000000 132.905 > 1 5 Cs 0.16666666666667 0.33333333333333 0.75000000000000 132.905 > 1 6 Cs 0.41666666666667 0.33333333333333 0.75000000000000 132.905 > 1 7 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 > 1 8 Cs 0.91666666666667 0.33333333333333 0.75000000000000 132.905 > 1 9 Cs 0.16666666666667 0.58333333333333 0.75000000000000 132.905 > 1 10 Cs 0.41666666666667 0.58333333333333 0.75000000000000 132.905 > 1 11 Cs 0.66666666666667 0.58333333333333 0.75000000000000 132.905 > 1 12 Cs 0.91666666666667 0.58333333333333 0.75000000000000 132.905 > 1 13 Cs 0.16666666666667 0.83333333333333 0.75000000000000 132.905 > 1 14 Cs 0.41666666666667 0.83333333333333 0.75000000000000 132.905 > 1 15 Cs 0.66666666666667 0.83333333333333 0.75000000000000 132.905 > 1 16 Cs 0.91666666666667 0.83333333333333 0.75000000000000 132.905 > 1 17 Cs 0.08333333333333 0.16666666666667 0.25000000000000 132.905 > 2 18 Cs 0.33333333333333 0.16666666666667 0.25000000000000 132.905 > 2 19 Cs 0.58333333333333 0.16666666666667 0.25000000000000 132.905 > 2 20 Cs 0.83333333333333 0.16666666666667 0.25000000000000 132.905 > 2 21 Cs 0.08333333333333 0.41666666666667 0.25000000000000 132.905 > 2 22 Cs 0.33333333333333 0.41666666666667 0.25000000000000 132.905 > 2 23 Cs 0.58333333333333 0.41666666666667 0.25000000000000 132.905 > 2 24 Cs 0.83333333333333 0.41666666666667 0.25000000000000 132.905 > 2 25 Cs 0.08333333333333 0.66666666666667 0.25000000000000 132.905 > 2 26 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 > 2 27 Cs 0.58333333333333 0.66666666666667 0.25000000000000 132.905 > 2 28 Cs 0.83333333333333 0.66666666666667 0.25000000000000 132.905 > 2 29 Cs 0.08333333333333 0.91666666666667 0.25000000000000 132.905 > 2 30 Cs 0.33333333333333 0.91666666666667 0.25000000000000 132.905 > 2 31 Cs 0.58333333333333 0.91666666666667 0.25000000000000 132.905 > 2 32 Cs 0.83333333333333 0.91666666666667 0.25000000000000 132.905 > 2 *33 Dy 0.00000000000000 0.00000000000000 0.50000000000000 162.500 > 3 34 Dy 0.25000000000000 0.00000000000000 0.50000000000000 162.500 > 3 35 Dy 0.50000000000000 0.00000000000000 0.50000000000000 162.500 > 3 36 Dy 0.75000000000000 0.00000000000000 0.50000000000000 162.500 > 3 37 Dy 0.00000000000000 0.25000000000000 0.50000000000000 162.500 > 3 38 Dy 0.25000000000000 0.25000000000000 0.50000000000000 162.500 > 3 39 Dy 0.50000000000000 0.25000000000000 0.50000000000000 162.500 > 3 40 Dy 0.75000000000000 0.25000000000000 0.50000000000000 162.500 > 3 41 Dy 0.00000000000000 0.50000000000000 0.50000000000000 162.500 > 3 42 Dy 0.25000000000000 0.50000000000000 0.50000000000000 162.500 > 3 43 Dy 0.50000000000000 0.50000000000000 0.50000000000000 162.500 > 3 44 Dy 0.75000000000000 0.50000000000000 0.50000000000000 162.500 > 3 45 Dy 0.00000000000000 0.75000000000000 0.50000000000000 162.500 > 3 46 Dy 0.25000000000000 0.75000000000000 0.50000000000000 162.500 > 3 47 Dy 0.50000000000000 0.75000000000000 0.50000000000000 162.500 > 3 48 Dy 0.75000000000000 0.75000000000000 0.50000000000000 162.500 > 3 49 Dy 0.00000000000000 0.00000000000000 0.00000000000000 162.500 > 4 50 Dy 0.25000000000000 0.00000000000000 0.00000000000000 162.500 > 4 51 Dy 0.50000000000000 0.00000000000000 0.00000000000000 162.500 > 4 52 Dy 0.75000000000000 0.00000000000000 0.00000000000000 162.500 > 4 53 Dy 0.00000000000000 0.25000000000000 0.00000000000000 162.500 > 4 54 Dy 0.25000000000000 0.25000000000000 0.00000000000000 162.500 > 4 55 Dy 0.50000000000000 0.25000000000000 0.00000000000000 162.500 > 4 56 Dy 0.75000000000000 0.25000000000000 0.00000000000000 162.500 > 4 57 Dy 0.00000000000000 0.50000000000000 0.00000000000000 162.500 > 4 58 Dy 0.25000000000000 0.50000000000000 0.00000000000000 162.500 > 4 59 Dy 0.50000000000000 0.50000000000000 0.00000000000000 162.500 > 4 60 Dy 0.75000000000000 0.50000000000000 0.00000000000000 162.500 > 4 61 Dy 0.00000000000000 0.75000000000000 0.00000000000000 162.500 > 4 62 Dy 0.25000000000000 0.75000000000000 0.00000000000000 162.500 > 4 63 Dy 0.50000000000000 0.75000000000000 0.00000000000000 162.500 > 4 64 Dy 0.75000000000000 0.75000000000000 0.00000000000000 162.500 > 4 *65 S 0.08333333333333 0.16666666666667 0.90806073924251 32.065 > 5 66 S 0.33333333333333 0.16666666666667 0.90806073924251 32.065 > 5 67 S 0.58333333333333 0.16666666666667 0.90806073924251 32.065 > 5 68 S 0.83333333333333 0.16666666666667 0.90806073924251 32.065 > 5 69 S 0.08333333333333 0.41666666666667 0.90806073924251 32.065 > 5 70 S 0.33333333333333 0.41666666666667 0.90806073924251 32.065 > 5 71 S 0.58333333333333 0.41666666666667 0.90806073924251 32.065 > 5 72 S 0.83333333333333 0.41666666666667 0.90806073924251 32.065 > 5 73 S 0.08333333333333 0.66666666666667 0.90806073924251 32.065 > 5 74 S 0.33333333333333 0.66666666666667 0.90806073924251 32.065 > 5 75 S 0.58333333333333 0.66666666666667 0.90806073924251 32.065 > 5 76 S 0.83333333333333 0.66666666666667 0.90806073924251 32.065 > 5 77 S 0.08333333333333 0.91666666666667 0.90806073924251 32.065 > 5 78 S 0.33333333333333 0.91666666666667 0.90806073924251 32.065 > 5 79 S 0.58333333333333 0.91666666666667 0.90806073924251 32.065 > 5 80 S 0.83333333333333 0.91666666666667 0.90806073924251 32.065 > 5 81 S 0.16666666666667 0.08333333333333 0.40806073924251 32.065 > 6 82 S 0.41666666666667 0.08333333333333 0.40806073924251 32.065 > 6 83 S 0.66666666666667 0.08333333333333 0.40806073924251 32.065 > 6 84 S 0.91666666666667 0.08333333333333 0.40806073924251 32.065 > 6 85 S 0.16666666666667 0.33333333333333 0.40806073924251 32.065 > 6 86 S 0.41666666666667 0.33333333333333 0.40806073924251 32.065 > 6 87 S 0.66666666666667 0.33333333333333 0.40806073924251 32.065 > 6 88 S 0.91666666666667 0.33333333333333 0.40806073924251 32.065 > 6 89 S 0.16666666666667 0.58333333333333 0.40806073924251 32.065 > 6 90 S 0.41666666666667 0.58333333333333 0.40806073924251 32.065 > 6 91 S 0.66666666666667 0.58333333333333 0.40806073924251 32.065 > 6 92 S 0.91666666666667 0.58333333333333 0.40806073924251 32.065 > 6 93 S 0.16666666666667 0.83333333333333 0.40806073924251 32.065 > 6 94 S 0.41666666666667 0.83333333333333 0.40806073924251 32.065 > 6 95 S 0.66666666666667 0.83333333333333 0.40806073924251 32.065 > 6 96 S 0.91666666666667 0.83333333333333 0.40806073924251 32.065 > 6 97 S 0.16666666666667 0.08333333333333 0.09193926075749 32.065 > 7 98 S 0.41666666666667 0.08333333333333 0.09193926075749 32.065 > 7 99 S 0.66666666666667 0.08333333333333 0.09193926075749 32.065 > 7 100 S 0.91666666666667 0.08333333333333 0.09193926075749 32.065 > 7 101 S 0.16666666666667 0.33333333333333 0.09193926075749 32.065 > 7 102 S 0.41666666666667 0.33333333333333 0.09193926075749 32.065 > 7 103 S 0.66666666666667 0.33333333333333 0.09193926075749 32.065 > 7 104 S 0.91666666666667 0.33333333333333 0.09193926075749 32.065 > 7 105 S 0.16666666666667 0.58333333333333 0.09193926075749 32.065 > 7 106 S 0.41666666666667 0.58333333333333 0.09193926075749 32.065 > 7 107 S 0.66666666666667 0.58333333333333 0.09193926075749 32.065 > 7 108 S 0.91666666666667 0.58333333333333 0.09193926075749 32.065 > 7 109 S 0.16666666666667 0.83333333333333 0.09193926075749 32.065 > 7 110 S 0.41666666666667 0.83333333333333 0.09193926075749 32.065 > 7 111 S 0.66666666666667 0.83333333333333 0.09193926075749 32.065 > 7 112 S 0.91666666666667 0.83333333333333 0.09193926075749 32.065 > 7 113 S 0.08333333333333 0.16666666666667 0.59193926075749 32.065 > 8 114 S 0.33333333333333 0.16666666666667 0.59193926075749 32.065 > 8 115 S 0.58333333333333 0.16666666666667 0.59193926075749 32.065 > 8 116 S 0.83333333333333 0.16666666666667 0.59193926075749 32.065 > 8 117 S 0.08333333333333 0.41666666666667 0.59193926075749 32.065 > 8 118 S 0.33333333333333 0.41666666666667 0.59193926075749 32.065 > 8 119 S 0.58333333333333 0.41666666666667 0.59193926075749 32.065 > 8 120 S 0.83333333333333 0.41666666666667 0.59193926075749 32.065 > 8 121 S 0.08333333333333 0.66666666666667 0.59193926075749 32.065 > 8 122 S 0.33333333333333 0.66666666666667 0.59193926075749 32.065 > 8 123 S 0.58333333333333 0.66666666666667 0.59193926075749 32.065 > 8 124 S 0.83333333333333 0.66666666666667 0.59193926075749 32.065 > 8 125 S 0.08333333333333 0.91666666666667 0.59193926075749 32.065 > 8 126 S 0.33333333333333 0.91666666666667 0.59193926075749 32.065 > 8 127 S 0.58333333333333 0.91666666666667 0.59193926075749 32.065 > 8 128 S 0.83333333333333 0.91666666666667 0.59193926075749 32.065 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.6604814 -0.0000000 0.0000000 0.0000000 5.6604814 0.0000000 0.0000000 0.0000000 5.0469838 -------------------------- Born effective charges -------------------------- 1 Cs 1.0703510 0.0000000 0.0000000 -0.0000000 1.0703510 0.0000000 0.0000000 0.0000000 2.2487040 2 Cs 1.0703510 0.0000000 0.0000000 -0.0000000 1.0703510 0.0000000 0.0000000 0.0000000 2.2487040 3 Dy 4.4466869 -0.0000000 0.0000000 0.0000000 4.4466869 0.0000000 0.0000000 0.0000000 2.1178773 4 Dy 4.4466869 -0.0000000 0.0000000 0.0000000 4.4466869 0.0000000 0.0000000 0.0000000 2.1178773 5 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 6 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 7 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 8 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 384/384 Permutation basis: 5136/5136 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 259 Number of blocks in projector: 213 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 123 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 76 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 60 Use standard eigh solver. Tree of FC basis block matrices: - (259, 253), data: False |-- (60, 57), data: True |-- (76, 76), data: True |-- (123, 120), data: True ----- Solver_atoms: 1 -- 128 / 128 Time (Solver_compr_matrix_reshape): 0.009 Solver_block: 80 / 80 - Time: 0.184 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.194 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (yy) 0.00000000 (yy) Permutation basis: 384/384 Permutation basis: 5136/5136 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 259 Number of blocks in projector: 213 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 123 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 76 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 60 Use standard eigh solver. Tree of FC basis block matrices: - (259, 253), data: False |-- (60, 57), data: True |-- (76, 76), data: True |-- (123, 120), data: True Max drift after symmetrization by symfc projector: 0.00000000 (yy) 0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-09 04:49:23]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:49:24]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.031293939063714 0.000000000000000 0.000000000000000 b -2.015646969531857 3.491202961351413 0.000000000000000 c 0.000000000000000 0.000000000000000 15.830686910000001 Atomic positions (fractional): 1 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 2 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 3 Dy 0.00000000000000 0.00000000000000 0.50000000000000 162.500 4 Dy 0.00000000000000 0.00000000000000 0.00000000000000 162.500 5 S 0.33333333333333 0.66666666666667 0.90806073924251 32.065 6 S 0.66666666666667 0.33333333333333 0.40806073924251 32.065 7 S 0.66666666666667 0.33333333333333 0.09193926075749 32.065 8 S 0.33333333333333 0.66666666666667 0.59193926075749 32.065 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.125175756254855 0.000000000000000 0.000000000000000 b -8.062587878127427 13.964811845405652 0.000000000000000 c 0.000000000000000 0.000000000000000 15.830686910000001 Atomic positions (fractional): 1 Cs 0.16666666666667 0.08333333333333 0.75000000000000 132.905 > 1 2 Cs 0.41666666666667 0.08333333333333 0.75000000000000 132.905 > 1 3 Cs 0.66666666666667 0.08333333333333 0.75000000000000 132.905 > 1 4 Cs 0.91666666666667 0.08333333333333 0.75000000000000 132.905 > 1 5 Cs 0.16666666666667 0.33333333333333 0.75000000000000 132.905 > 1 6 Cs 0.41666666666667 0.33333333333333 0.75000000000000 132.905 > 1 7 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 > 1 8 Cs 0.91666666666667 0.33333333333333 0.75000000000000 132.905 > 1 9 Cs 0.16666666666667 0.58333333333333 0.75000000000000 132.905 > 1 10 Cs 0.41666666666667 0.58333333333333 0.75000000000000 132.905 > 1 11 Cs 0.66666666666667 0.58333333333333 0.75000000000000 132.905 > 1 12 Cs 0.91666666666667 0.58333333333333 0.75000000000000 132.905 > 1 13 Cs 0.16666666666667 0.83333333333333 0.75000000000000 132.905 > 1 14 Cs 0.41666666666667 0.83333333333333 0.75000000000000 132.905 > 1 15 Cs 0.66666666666667 0.83333333333333 0.75000000000000 132.905 > 1 16 Cs 0.91666666666667 0.83333333333333 0.75000000000000 132.905 > 1 17 Cs 0.08333333333333 0.16666666666667 0.25000000000000 132.905 > 17 18 Cs 0.33333333333333 0.16666666666667 0.25000000000000 132.905 > 17 19 Cs 0.58333333333333 0.16666666666667 0.25000000000000 132.905 > 17 20 Cs 0.83333333333333 0.16666666666667 0.25000000000000 132.905 > 17 21 Cs 0.08333333333333 0.41666666666667 0.25000000000000 132.905 > 17 22 Cs 0.33333333333333 0.41666666666667 0.25000000000000 132.905 > 17 23 Cs 0.58333333333333 0.41666666666667 0.25000000000000 132.905 > 17 24 Cs 0.83333333333333 0.41666666666667 0.25000000000000 132.905 > 17 25 Cs 0.08333333333333 0.66666666666667 0.25000000000000 132.905 > 17 26 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 > 17 27 Cs 0.58333333333333 0.66666666666667 0.25000000000000 132.905 > 17 28 Cs 0.83333333333333 0.66666666666667 0.25000000000000 132.905 > 17 29 Cs 0.08333333333333 0.91666666666667 0.25000000000000 132.905 > 17 30 Cs 0.33333333333333 0.91666666666667 0.25000000000000 132.905 > 17 31 Cs 0.58333333333333 0.91666666666667 0.25000000000000 132.905 > 17 32 Cs 0.83333333333333 0.91666666666667 0.25000000000000 132.905 > 17 33 Dy 0.00000000000000 0.00000000000000 0.50000000000000 162.500 > 33 34 Dy 0.25000000000000 0.00000000000000 0.50000000000000 162.500 > 33 35 Dy 0.50000000000000 0.00000000000000 0.50000000000000 162.500 > 33 36 Dy 0.75000000000000 0.00000000000000 0.50000000000000 162.500 > 33 37 Dy 0.00000000000000 0.25000000000000 0.50000000000000 162.500 > 33 38 Dy 0.25000000000000 0.25000000000000 0.50000000000000 162.500 > 33 39 Dy 0.50000000000000 0.25000000000000 0.50000000000000 162.500 > 33 40 Dy 0.75000000000000 0.25000000000000 0.50000000000000 162.500 > 33 41 Dy 0.00000000000000 0.50000000000000 0.50000000000000 162.500 > 33 42 Dy 0.25000000000000 0.50000000000000 0.50000000000000 162.500 > 33 43 Dy 0.50000000000000 0.50000000000000 0.50000000000000 162.500 > 33 44 Dy 0.75000000000000 0.50000000000000 0.50000000000000 162.500 > 33 45 Dy 0.00000000000000 0.75000000000000 0.50000000000000 162.500 > 33 46 Dy 0.25000000000000 0.75000000000000 0.50000000000000 162.500 > 33 47 Dy 0.50000000000000 0.75000000000000 0.50000000000000 162.500 > 33 48 Dy 0.75000000000000 0.75000000000000 0.50000000000000 162.500 > 33 49 Dy 0.00000000000000 0.00000000000000 0.00000000000000 162.500 > 49 50 Dy 0.25000000000000 0.00000000000000 0.00000000000000 162.500 > 49 51 Dy 0.50000000000000 0.00000000000000 0.00000000000000 162.500 > 49 52 Dy 0.75000000000000 0.00000000000000 0.00000000000000 162.500 > 49 53 Dy 0.00000000000000 0.25000000000000 0.00000000000000 162.500 > 49 54 Dy 0.25000000000000 0.25000000000000 0.00000000000000 162.500 > 49 55 Dy 0.50000000000000 0.25000000000000 0.00000000000000 162.500 > 49 56 Dy 0.75000000000000 0.25000000000000 0.00000000000000 162.500 > 49 57 Dy 0.00000000000000 0.50000000000000 0.00000000000000 162.500 > 49 58 Dy 0.25000000000000 0.50000000000000 0.00000000000000 162.500 > 49 59 Dy 0.50000000000000 0.50000000000000 0.00000000000000 162.500 > 49 60 Dy 0.75000000000000 0.50000000000000 0.00000000000000 162.500 > 49 61 Dy 0.00000000000000 0.75000000000000 0.00000000000000 162.500 > 49 62 Dy 0.25000000000000 0.75000000000000 0.00000000000000 162.500 > 49 63 Dy 0.50000000000000 0.75000000000000 0.00000000000000 162.500 > 49 64 Dy 0.75000000000000 0.75000000000000 0.00000000000000 162.500 > 49 65 S 0.08333333333333 0.16666666666667 0.90806073924251 32.065 > 65 66 S 0.33333333333333 0.16666666666667 0.90806073924251 32.065 > 65 67 S 0.58333333333333 0.16666666666667 0.90806073924251 32.065 > 65 68 S 0.83333333333333 0.16666666666667 0.90806073924251 32.065 > 65 69 S 0.08333333333333 0.41666666666667 0.90806073924251 32.065 > 65 70 S 0.33333333333333 0.41666666666667 0.90806073924251 32.065 > 65 71 S 0.58333333333333 0.41666666666667 0.90806073924251 32.065 > 65 72 S 0.83333333333333 0.41666666666667 0.90806073924251 32.065 > 65 73 S 0.08333333333333 0.66666666666667 0.90806073924251 32.065 > 65 74 S 0.33333333333333 0.66666666666667 0.90806073924251 32.065 > 65 75 S 0.58333333333333 0.66666666666667 0.90806073924251 32.065 > 65 76 S 0.83333333333333 0.66666666666667 0.90806073924251 32.065 > 65 77 S 0.08333333333333 0.91666666666667 0.90806073924251 32.065 > 65 78 S 0.33333333333333 0.91666666666667 0.90806073924251 32.065 > 65 79 S 0.58333333333333 0.91666666666667 0.90806073924251 32.065 > 65 80 S 0.83333333333333 0.91666666666667 0.90806073924251 32.065 > 65 81 S 0.16666666666667 0.08333333333333 0.40806073924251 32.065 > 81 82 S 0.41666666666667 0.08333333333333 0.40806073924251 32.065 > 81 83 S 0.66666666666667 0.08333333333333 0.40806073924251 32.065 > 81 84 S 0.91666666666667 0.08333333333333 0.40806073924251 32.065 > 81 85 S 0.16666666666667 0.33333333333333 0.40806073924251 32.065 > 81 86 S 0.41666666666667 0.33333333333333 0.40806073924251 32.065 > 81 87 S 0.66666666666667 0.33333333333333 0.40806073924251 32.065 > 81 88 S 0.91666666666667 0.33333333333333 0.40806073924251 32.065 > 81 89 S 0.16666666666667 0.58333333333333 0.40806073924251 32.065 > 81 90 S 0.41666666666667 0.58333333333333 0.40806073924251 32.065 > 81 91 S 0.66666666666667 0.58333333333333 0.40806073924251 32.065 > 81 92 S 0.91666666666667 0.58333333333333 0.40806073924251 32.065 > 81 93 S 0.16666666666667 0.83333333333333 0.40806073924251 32.065 > 81 94 S 0.41666666666667 0.83333333333333 0.40806073924251 32.065 > 81 95 S 0.66666666666667 0.83333333333333 0.40806073924251 32.065 > 81 96 S 0.91666666666667 0.83333333333333 0.40806073924251 32.065 > 81 97 S 0.16666666666667 0.08333333333333 0.09193926075749 32.065 > 97 98 S 0.41666666666667 0.08333333333333 0.09193926075749 32.065 > 97 99 S 0.66666666666667 0.08333333333333 0.09193926075749 32.065 > 97 100 S 0.91666666666667 0.08333333333333 0.09193926075749 32.065 > 97 101 S 0.16666666666667 0.33333333333333 0.09193926075749 32.065 > 97 102 S 0.41666666666667 0.33333333333333 0.09193926075749 32.065 > 97 103 S 0.66666666666667 0.33333333333333 0.09193926075749 32.065 > 97 104 S 0.91666666666667 0.33333333333333 0.09193926075749 32.065 > 97 105 S 0.16666666666667 0.58333333333333 0.09193926075749 32.065 > 97 106 S 0.41666666666667 0.58333333333333 0.09193926075749 32.065 > 97 107 S 0.66666666666667 0.58333333333333 0.09193926075749 32.065 > 97 108 S 0.91666666666667 0.58333333333333 0.09193926075749 32.065 > 97 109 S 0.16666666666667 0.83333333333333 0.09193926075749 32.065 > 97 110 S 0.41666666666667 0.83333333333333 0.09193926075749 32.065 > 97 111 S 0.66666666666667 0.83333333333333 0.09193926075749 32.065 > 97 112 S 0.91666666666667 0.83333333333333 0.09193926075749 32.065 > 97 113 S 0.08333333333333 0.16666666666667 0.59193926075749 32.065 > 113 114 S 0.33333333333333 0.16666666666667 0.59193926075749 32.065 > 113 115 S 0.58333333333333 0.16666666666667 0.59193926075749 32.065 > 113 116 S 0.83333333333333 0.16666666666667 0.59193926075749 32.065 > 113 117 S 0.08333333333333 0.41666666666667 0.59193926075749 32.065 > 113 118 S 0.33333333333333 0.41666666666667 0.59193926075749 32.065 > 113 119 S 0.58333333333333 0.41666666666667 0.59193926075749 32.065 > 113 120 S 0.83333333333333 0.41666666666667 0.59193926075749 32.065 > 113 121 S 0.08333333333333 0.66666666666667 0.59193926075749 32.065 > 113 122 S 0.33333333333333 0.66666666666667 0.59193926075749 32.065 > 113 123 S 0.58333333333333 0.66666666666667 0.59193926075749 32.065 > 113 124 S 0.83333333333333 0.66666666666667 0.59193926075749 32.065 > 113 125 S 0.08333333333333 0.91666666666667 0.59193926075749 32.065 > 113 126 S 0.33333333333333 0.91666666666667 0.59193926075749 32.065 > 113 127 S 0.58333333333333 0.91666666666667 0.59193926075749 32.065 > 113 128 S 0.83333333333333 0.91666666666667 0.59193926075749 32.065 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.6604814 -0.0000000 0.0000000 0.0000000 5.6604814 0.0000000 0.0000000 0.0000000 5.0469838 -------------------------- Born effective charges -------------------------- 1 Cs 1.0703510 0.0000000 0.0000000 -0.0000000 1.0703510 0.0000000 0.0000000 0.0000000 2.2487040 2 Cs 1.0703510 0.0000000 0.0000000 -0.0000000 1.0703510 0.0000000 0.0000000 0.0000000 2.2487040 3 Dy 4.4466869 -0.0000000 0.0000000 0.0000000 4.4466869 0.0000000 0.0000000 0.0000000 2.1178773 4 Dy 4.4466869 -0.0000000 0.0000000 0.0000000 4.4466869 0.0000000 0.0000000 0.0000000 2.1178773 5 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 6 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 7 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 8 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000009 (xxy) 0.00000009 (xxy) 0.00000009 (xyx) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-09 04:49:31]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:49:31]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P6_3/mmc (194) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.031293939063714 0.000000000000000 0.000000000000000 b -2.015646969531857 3.491202961351413 0.000000000000000 c 0.000000000000000 0.000000000000000 15.830686910000001 Atomic positions (fractional): 1 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 2 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 3 Dy 0.00000000000000 0.00000000000000 0.50000000000000 162.500 4 Dy 0.00000000000000 0.00000000000000 0.00000000000000 162.500 5 S 0.33333333333333 0.66666666666667 0.90806073924251 32.065 6 S 0.66666666666667 0.33333333333333 0.40806073924251 32.065 7 S 0.66666666666667 0.33333333333333 0.09193926075749 32.065 8 S 0.33333333333333 0.66666666666667 0.59193926075749 32.065 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.125175756254855 0.000000000000000 0.000000000000000 b -8.062587878127427 13.964811845405652 0.000000000000000 c 0.000000000000000 0.000000000000000 15.830686910000001 Atomic positions (fractional): 1 Cs 0.16666666666667 0.08333333333333 0.75000000000000 132.905 > 1 2 Cs 0.41666666666667 0.08333333333333 0.75000000000000 132.905 > 1 3 Cs 0.66666666666667 0.08333333333333 0.75000000000000 132.905 > 1 4 Cs 0.91666666666667 0.08333333333333 0.75000000000000 132.905 > 1 5 Cs 0.16666666666667 0.33333333333333 0.75000000000000 132.905 > 1 6 Cs 0.41666666666667 0.33333333333333 0.75000000000000 132.905 > 1 7 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 > 1 8 Cs 0.91666666666667 0.33333333333333 0.75000000000000 132.905 > 1 9 Cs 0.16666666666667 0.58333333333333 0.75000000000000 132.905 > 1 10 Cs 0.41666666666667 0.58333333333333 0.75000000000000 132.905 > 1 11 Cs 0.66666666666667 0.58333333333333 0.75000000000000 132.905 > 1 12 Cs 0.91666666666667 0.58333333333333 0.75000000000000 132.905 > 1 13 Cs 0.16666666666667 0.83333333333333 0.75000000000000 132.905 > 1 14 Cs 0.41666666666667 0.83333333333333 0.75000000000000 132.905 > 1 15 Cs 0.66666666666667 0.83333333333333 0.75000000000000 132.905 > 1 16 Cs 0.91666666666667 0.83333333333333 0.75000000000000 132.905 > 1 17 Cs 0.08333333333333 0.16666666666667 0.25000000000000 132.905 > 17 18 Cs 0.33333333333333 0.16666666666667 0.25000000000000 132.905 > 17 19 Cs 0.58333333333333 0.16666666666667 0.25000000000000 132.905 > 17 20 Cs 0.83333333333333 0.16666666666667 0.25000000000000 132.905 > 17 21 Cs 0.08333333333333 0.41666666666667 0.25000000000000 132.905 > 17 22 Cs 0.33333333333333 0.41666666666667 0.25000000000000 132.905 > 17 23 Cs 0.58333333333333 0.41666666666667 0.25000000000000 132.905 > 17 24 Cs 0.83333333333333 0.41666666666667 0.25000000000000 132.905 > 17 25 Cs 0.08333333333333 0.66666666666667 0.25000000000000 132.905 > 17 26 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 > 17 27 Cs 0.58333333333333 0.66666666666667 0.25000000000000 132.905 > 17 28 Cs 0.83333333333333 0.66666666666667 0.25000000000000 132.905 > 17 29 Cs 0.08333333333333 0.91666666666667 0.25000000000000 132.905 > 17 30 Cs 0.33333333333333 0.91666666666667 0.25000000000000 132.905 > 17 31 Cs 0.58333333333333 0.91666666666667 0.25000000000000 132.905 > 17 32 Cs 0.83333333333333 0.91666666666667 0.25000000000000 132.905 > 17 33 Dy 0.00000000000000 0.00000000000000 0.50000000000000 162.500 > 33 34 Dy 0.25000000000000 0.00000000000000 0.50000000000000 162.500 > 33 35 Dy 0.50000000000000 0.00000000000000 0.50000000000000 162.500 > 33 36 Dy 0.75000000000000 0.00000000000000 0.50000000000000 162.500 > 33 37 Dy 0.00000000000000 0.25000000000000 0.50000000000000 162.500 > 33 38 Dy 0.25000000000000 0.25000000000000 0.50000000000000 162.500 > 33 39 Dy 0.50000000000000 0.25000000000000 0.50000000000000 162.500 > 33 40 Dy 0.75000000000000 0.25000000000000 0.50000000000000 162.500 > 33 41 Dy 0.00000000000000 0.50000000000000 0.50000000000000 162.500 > 33 42 Dy 0.25000000000000 0.50000000000000 0.50000000000000 162.500 > 33 43 Dy 0.50000000000000 0.50000000000000 0.50000000000000 162.500 > 33 44 Dy 0.75000000000000 0.50000000000000 0.50000000000000 162.500 > 33 45 Dy 0.00000000000000 0.75000000000000 0.50000000000000 162.500 > 33 46 Dy 0.25000000000000 0.75000000000000 0.50000000000000 162.500 > 33 47 Dy 0.50000000000000 0.75000000000000 0.50000000000000 162.500 > 33 48 Dy 0.75000000000000 0.75000000000000 0.50000000000000 162.500 > 33 49 Dy 0.00000000000000 0.00000000000000 0.00000000000000 162.500 > 49 50 Dy 0.25000000000000 0.00000000000000 0.00000000000000 162.500 > 49 51 Dy 0.50000000000000 0.00000000000000 0.00000000000000 162.500 > 49 52 Dy 0.75000000000000 0.00000000000000 0.00000000000000 162.500 > 49 53 Dy 0.00000000000000 0.25000000000000 0.00000000000000 162.500 > 49 54 Dy 0.25000000000000 0.25000000000000 0.00000000000000 162.500 > 49 55 Dy 0.50000000000000 0.25000000000000 0.00000000000000 162.500 > 49 56 Dy 0.75000000000000 0.25000000000000 0.00000000000000 162.500 > 49 57 Dy 0.00000000000000 0.50000000000000 0.00000000000000 162.500 > 49 58 Dy 0.25000000000000 0.50000000000000 0.00000000000000 162.500 > 49 59 Dy 0.50000000000000 0.50000000000000 0.00000000000000 162.500 > 49 60 Dy 0.75000000000000 0.50000000000000 0.00000000000000 162.500 > 49 61 Dy 0.00000000000000 0.75000000000000 0.00000000000000 162.500 > 49 62 Dy 0.25000000000000 0.75000000000000 0.00000000000000 162.500 > 49 63 Dy 0.50000000000000 0.75000000000000 0.00000000000000 162.500 > 49 64 Dy 0.75000000000000 0.75000000000000 0.00000000000000 162.500 > 49 65 S 0.08333333333333 0.16666666666667 0.90806073924251 32.065 > 65 66 S 0.33333333333333 0.16666666666667 0.90806073924251 32.065 > 65 67 S 0.58333333333333 0.16666666666667 0.90806073924251 32.065 > 65 68 S 0.83333333333333 0.16666666666667 0.90806073924251 32.065 > 65 69 S 0.08333333333333 0.41666666666667 0.90806073924251 32.065 > 65 70 S 0.33333333333333 0.41666666666667 0.90806073924251 32.065 > 65 71 S 0.58333333333333 0.41666666666667 0.90806073924251 32.065 > 65 72 S 0.83333333333333 0.41666666666667 0.90806073924251 32.065 > 65 73 S 0.08333333333333 0.66666666666667 0.90806073924251 32.065 > 65 74 S 0.33333333333333 0.66666666666667 0.90806073924251 32.065 > 65 75 S 0.58333333333333 0.66666666666667 0.90806073924251 32.065 > 65 76 S 0.83333333333333 0.66666666666667 0.90806073924251 32.065 > 65 77 S 0.08333333333333 0.91666666666667 0.90806073924251 32.065 > 65 78 S 0.33333333333333 0.91666666666667 0.90806073924251 32.065 > 65 79 S 0.58333333333333 0.91666666666667 0.90806073924251 32.065 > 65 80 S 0.83333333333333 0.91666666666667 0.90806073924251 32.065 > 65 81 S 0.16666666666667 0.08333333333333 0.40806073924251 32.065 > 81 82 S 0.41666666666667 0.08333333333333 0.40806073924251 32.065 > 81 83 S 0.66666666666667 0.08333333333333 0.40806073924251 32.065 > 81 84 S 0.91666666666667 0.08333333333333 0.40806073924251 32.065 > 81 85 S 0.16666666666667 0.33333333333333 0.40806073924251 32.065 > 81 86 S 0.41666666666667 0.33333333333333 0.40806073924251 32.065 > 81 87 S 0.66666666666667 0.33333333333333 0.40806073924251 32.065 > 81 88 S 0.91666666666667 0.33333333333333 0.40806073924251 32.065 > 81 89 S 0.16666666666667 0.58333333333333 0.40806073924251 32.065 > 81 90 S 0.41666666666667 0.58333333333333 0.40806073924251 32.065 > 81 91 S 0.66666666666667 0.58333333333333 0.40806073924251 32.065 > 81 92 S 0.91666666666667 0.58333333333333 0.40806073924251 32.065 > 81 93 S 0.16666666666667 0.83333333333333 0.40806073924251 32.065 > 81 94 S 0.41666666666667 0.83333333333333 0.40806073924251 32.065 > 81 95 S 0.66666666666667 0.83333333333333 0.40806073924251 32.065 > 81 96 S 0.91666666666667 0.83333333333333 0.40806073924251 32.065 > 81 97 S 0.16666666666667 0.08333333333333 0.09193926075749 32.065 > 97 98 S 0.41666666666667 0.08333333333333 0.09193926075749 32.065 > 97 99 S 0.66666666666667 0.08333333333333 0.09193926075749 32.065 > 97 100 S 0.91666666666667 0.08333333333333 0.09193926075749 32.065 > 97 101 S 0.16666666666667 0.33333333333333 0.09193926075749 32.065 > 97 102 S 0.41666666666667 0.33333333333333 0.09193926075749 32.065 > 97 103 S 0.66666666666667 0.33333333333333 0.09193926075749 32.065 > 97 104 S 0.91666666666667 0.33333333333333 0.09193926075749 32.065 > 97 105 S 0.16666666666667 0.58333333333333 0.09193926075749 32.065 > 97 106 S 0.41666666666667 0.58333333333333 0.09193926075749 32.065 > 97 107 S 0.66666666666667 0.58333333333333 0.09193926075749 32.065 > 97 108 S 0.91666666666667 0.58333333333333 0.09193926075749 32.065 > 97 109 S 0.16666666666667 0.83333333333333 0.09193926075749 32.065 > 97 110 S 0.41666666666667 0.83333333333333 0.09193926075749 32.065 > 97 111 S 0.66666666666667 0.83333333333333 0.09193926075749 32.065 > 97 112 S 0.91666666666667 0.83333333333333 0.09193926075749 32.065 > 97 113 S 0.08333333333333 0.16666666666667 0.59193926075749 32.065 > 113 114 S 0.33333333333333 0.16666666666667 0.59193926075749 32.065 > 113 115 S 0.58333333333333 0.16666666666667 0.59193926075749 32.065 > 113 116 S 0.83333333333333 0.16666666666667 0.59193926075749 32.065 > 113 117 S 0.08333333333333 0.41666666666667 0.59193926075749 32.065 > 113 118 S 0.33333333333333 0.41666666666667 0.59193926075749 32.065 > 113 119 S 0.58333333333333 0.41666666666667 0.59193926075749 32.065 > 113 120 S 0.83333333333333 0.41666666666667 0.59193926075749 32.065 > 113 121 S 0.08333333333333 0.66666666666667 0.59193926075749 32.065 > 113 122 S 0.33333333333333 0.66666666666667 0.59193926075749 32.065 > 113 123 S 0.58333333333333 0.66666666666667 0.59193926075749 32.065 > 113 124 S 0.83333333333333 0.66666666666667 0.59193926075749 32.065 > 113 125 S 0.08333333333333 0.91666666666667 0.59193926075749 32.065 > 113 126 S 0.33333333333333 0.91666666666667 0.59193926075749 32.065 > 113 127 S 0.58333333333333 0.91666666666667 0.59193926075749 32.065 > 113 128 S 0.83333333333333 0.91666666666667 0.59193926075749 32.065 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 5.6604814 -0.0000000 0.0000000 0.0000000 5.6604814 0.0000000 0.0000000 0.0000000 5.0469838 -------------------------- Born effective charges -------------------------- 1 Cs 1.0703510 0.0000000 0.0000000 -0.0000000 1.0703510 0.0000000 0.0000000 0.0000000 2.2487040 2 Cs 1.0703510 0.0000000 0.0000000 -0.0000000 1.0703510 0.0000000 0.0000000 0.0000000 2.2487040 3 Dy 4.4466869 -0.0000000 0.0000000 0.0000000 4.4466869 0.0000000 0.0000000 0.0000000 2.1178773 4 Dy 4.4466869 -0.0000000 0.0000000 0.0000000 4.4466869 0.0000000 0.0000000 0.0000000 2.1178773 5 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 6 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 7 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 8 S -2.7585189 0.0000000 0.0000000 -0.0000000 -2.7585189 0.0000000 0.0000000 0.0000000 -2.1832907 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000009 (xxy) 0.00000009 (xxy) 0.00000009 (xyx) Max drift of fc2: 0.00000000 (yy) 0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 14 14 3 ] 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.69, Number of G-points: 291, Lambda: 0.17 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/48) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 48 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.101 ( 0.000 0.000 0.000) 0.000 1.101 ( 0.000 0.000 0.000) 0.000 1.417 ( 0.000 0.000 0.000) 0.000 1.417 ( 0.000 0.000 0.000) 0.000 1.784 ( 0.000 0.000 0.000) 0.000 1.784 ( 0.000 0.000 0.000) 0.000 2.200 ( 0.000 0.000 0.000) 0.000 2.691 ( 0.000 0.000 0.000) 0.000 3.086 ( 0.000 0.000 0.000) 0.000 4.930 ( 0.000 0.000 0.000) 0.000 4.930 ( 0.000 0.000 0.000) 0.000 4.953 ( 0.000 0.000 0.000) 0.000 4.953 ( 0.000 0.000 0.000) 0.000 6.336 ( 0.000 0.000 0.000) 0.000 6.336 ( 0.000 0.000 0.000) 0.000 6.354 ( 0.000 0.000 0.000) 0.000 6.354 ( 0.000 0.000 0.000) 0.000 6.783 ( 0.000 0.000 0.000) 0.000 8.112 ( 0.000 0.000 0.000) 0.000 8.234 ( 0.000 0.000 0.000) 0.000 8.409 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/48) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.438 ( 18.379 10.611 0.000) 21.223 0.466 ( 18.525 10.696 0.000) 21.391 0.800 ( 32.073 18.518 0.000) 37.035 1.218 ( 9.033 5.215 0.000) 10.430 1.319 ( 15.503 8.950 0.000) 17.901 1.427 ( 0.909 0.525 0.000) 1.049 1.433 ( 0.664 0.383 0.000) 0.767 1.804 ( 1.754 1.013 0.000) 2.026 1.869 ( 7.278 4.202 0.000) 8.404 2.219 ( 2.593 1.497 0.000) 2.995 2.625 ( -2.776 -1.603 0.000) 3.205 3.064 ( -1.849 -1.068 0.000) 2.136 4.983 ( 4.360 2.517 0.000) 5.035 5.002 ( 4.084 2.358 0.000) 4.716 5.411 ( 33.336 19.247 0.000) 38.493 6.279 ( -4.569 -2.638 0.000) 5.276 6.307 ( -3.728 -2.153 0.000) 4.305 6.317 ( -1.552 -0.896 0.000) 1.792 6.335 ( -1.604 -0.926 0.000) 1.852 6.796 ( 1.005 0.580 0.000) 1.161 7.984 ( -7.851 -4.533 0.000) 9.065 8.193 ( -3.162 -1.826 0.000) 3.651 8.239 ( -2.133 -1.231 0.000) 2.463 8.277 ( -6.471 -3.736 0.000) 7.472 ======================= Grid point 2 (3/48) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.855 ( 14.178 8.186 0.000) 16.371 0.864 ( 17.488 10.096 0.000) 20.193 1.408 ( 4.572 2.640 0.000) 5.279 1.458 ( 1.624 0.938 0.000) 1.875 1.475 ( 11.745 6.781 0.000) 13.562 1.480 ( 17.310 9.994 0.000) 19.988 1.681 ( 12.709 7.337 0.000) 14.675 1.867 ( 3.543 2.046 0.000) 4.091 2.142 ( 16.362 9.447 0.000) 18.893 2.354 ( 8.478 4.895 0.000) 9.789 2.658 ( 6.443 3.720 0.000) 7.440 2.999 ( -3.494 -2.017 0.000) 4.034 5.123 ( 7.039 4.064 0.000) 8.128 5.134 ( 6.642 3.835 0.000) 7.670 6.126 ( 0.119 0.069 0.000) 0.138 6.179 ( 19.112 11.034 0.000) 22.068 6.197 ( -5.231 -3.020 0.000) 6.040 6.268 ( -2.376 -1.371 0.000) 2.743 6.284 ( -2.504 -1.446 0.000) 2.892 6.827 ( 1.725 0.996 0.000) 1.992 7.762 ( -11.025 -6.365 0.000) 12.730 7.920 ( -16.189 -9.347 0.000) 18.694 8.115 ( -2.621 -1.513 0.000) 3.026 8.299 ( 1.017 0.587 0.000) 1.174 ======================= Grid point 3 (4/48) ======================= q-point: ( 0.21 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.134 ( 9.654 5.573 0.000) 11.147 1.252 ( 15.088 8.711 0.000) 17.422 1.438 ( 2.731 1.577 0.000) 3.154 1.500 ( 1.865 1.077 0.000) 2.153 1.744 ( 10.740 6.201 0.000) 12.402 1.854 ( 11.305 6.527 0.000) 13.054 1.876 ( 4.418 2.551 0.000) 5.102 1.969 ( 5.001 2.888 0.000) 5.775 2.569 ( 9.050 5.225 0.000) 10.450 2.625 ( 22.474 12.975 0.000) 25.951 2.843 ( -0.366 -0.211 0.000) 0.422 2.976 ( 9.868 5.697 0.000) 11.394 5.292 ( 6.756 3.901 0.000) 7.801 5.294 ( 6.478 3.740 0.000) 7.480 5.978 ( -6.004 -3.466 0.000) 6.933 6.079 ( -4.314 -2.490 0.000) 4.981 6.216 ( -1.863 -1.075 0.000) 2.151 6.227 ( -2.154 -1.244 0.000) 2.487 6.554 ( 7.168 4.139 0.000) 8.277 6.881 ( 2.533 1.463 0.000) 2.925 7.450 ( -15.070 -8.700 0.000) 17.401 7.563 ( -13.426 -7.752 0.000) 15.503 8.095 ( 0.889 0.513 0.000) 1.027 8.304 ( -0.519 -0.300 0.000) 0.599 ======================= Grid point 4 (5/48) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.323 ( 6.568 3.792 0.000) 7.584 1.542 ( 1.610 0.930 0.000) 1.859 1.556 ( 6.913 3.991 0.000) 7.982 1.564 ( 11.168 6.448 0.000) 12.896 1.935 ( 1.331 0.768 0.000) 1.537 1.975 ( 8.756 5.055 0.000) 10.111 2.050 ( 5.867 3.387 0.000) 6.774 2.095 ( 5.433 3.137 0.000) 6.273 2.721 ( -6.711 -3.875 0.000) 7.750 2.725 ( 3.851 2.223 0.000) 4.446 3.156 ( 21.436 12.376 0.000) 24.752 3.325 ( 15.620 9.018 0.000) 18.037 5.425 ( 4.487 2.591 0.000) 5.181 5.427 ( 4.589 2.649 0.000) 5.299 5.876 ( -2.230 -1.288 0.000) 2.575 6.011 ( -1.149 -0.663 0.000) 1.327 6.185 ( -0.799 -0.461 0.000) 0.922 6.187 ( -1.239 -0.715 0.000) 1.431 6.603 ( -1.268 -0.732 0.000) 1.464 6.915 ( 0.236 0.136 0.000) 0.272 7.108 ( -13.196 -7.619 0.000) 15.238 7.267 ( -11.689 -6.749 0.000) 13.498 8.125 ( 0.760 0.439 0.000) 0.878 8.273 ( -2.106 -1.216 0.000) 2.432 ======================= Grid point 5 (6/48) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.452 ( 4.372 2.524 0.000) 5.048 1.575 ( 1.220 0.705 0.000) 1.409 1.742 ( 8.213 4.742 0.000) 9.484 1.775 ( 6.644 3.836 0.000) 7.672 1.957 ( 0.751 0.434 0.000) 0.867 2.151 ( 3.023 1.745 0.000) 3.490 2.156 ( 6.429 3.712 0.000) 7.424 2.216 ( 4.658 2.689 0.000) 5.378 2.559 ( -6.657 -3.843 0.000) 7.687 2.743 ( -2.189 -1.264 0.000) 2.527 3.612 ( 16.789 9.693 0.000) 19.387 3.666 ( 12.733 7.351 0.000) 14.703 5.511 ( 2.966 1.712 0.000) 3.424 5.515 ( 3.107 1.794 0.000) 3.587 5.880 ( 2.449 1.414 0.000) 2.828 6.037 ( 3.427 1.979 0.000) 3.957 6.161 ( -1.180 -0.681 0.000) 1.362 6.168 ( -0.842 -0.486 0.000) 0.972 6.540 ( -3.616 -2.088 0.000) 4.176 6.827 ( -10.923 -6.306 0.000) 12.613 6.907 ( -0.628 -0.363 0.000) 0.725 7.027 ( -8.218 -4.744 0.000) 9.489 8.110 ( -1.979 -1.143 0.000) 2.286 8.208 ( -3.282 -1.895 0.000) 3.790 ======================= Grid point 6 (7/48) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.528 ( 2.059 1.189 0.000) 2.378 1.600 ( 0.878 0.507 0.000) 1.014 1.877 ( 2.128 1.228 0.000) 2.457 1.911 ( 5.351 3.090 0.000) 6.179 1.980 ( 1.248 0.721 0.000) 1.442 2.214 ( 2.767 1.597 0.000) 3.195 2.274 ( 3.513 2.028 0.000) 4.057 2.306 ( 2.759 1.593 0.000) 3.185 2.425 ( -4.274 -2.468 0.000) 4.936 2.646 ( -5.153 -2.975 0.000) 5.951 3.907 ( 7.334 4.234 0.000) 8.468 3.928 ( 9.447 5.454 0.000) 10.908 5.569 ( 1.922 1.110 0.000) 2.219 5.578 ( 2.118 1.223 0.000) 2.445 5.971 ( 4.497 2.597 0.000) 5.193 6.130 ( -1.259 -0.727 0.000) 1.454 6.143 ( -1.148 -0.663 0.000) 1.326 6.166 ( 7.196 4.154 0.000) 8.309 6.448 ( -3.772 -2.178 0.000) 4.356 6.582 ( -9.877 -5.703 0.000) 11.405 6.886 ( -3.770 -2.177 0.000) 4.354 6.893 ( -0.529 -0.306 0.000) 0.611 8.049 ( -2.554 -1.475 0.000) 2.949 8.132 ( -2.713 -1.566 0.000) 3.133 ======================= Grid point 7 (8/48) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 64 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.551 ( 0.000 0.000 0.000) 0.000 1.612 ( 0.000 0.000 0.000) 0.000 1.897 ( 0.000 0.000 0.000) 0.000 1.974 ( 0.000 0.000 0.000) 0.000 2.000 ( 0.000 0.000 0.000) 0.000 2.264 ( 0.000 0.000 0.000) 0.000 2.317 ( 0.000 0.000 0.000) 0.000 2.340 ( 0.000 0.000 0.000) 0.000 2.372 ( 0.000 0.000 0.000) 0.000 2.564 ( 0.000 0.000 0.000) 0.000 3.996 ( 0.000 0.000 0.000) 0.000 4.042 ( 0.000 0.000 0.000) 0.000 5.594 ( 0.000 0.000 0.000) 0.000 5.605 ( 0.000 0.000 0.000) 0.000 6.039 ( 0.000 0.000 0.000) 0.000 6.113 ( 0.000 0.000 0.000) 0.000 6.126 ( 0.000 0.000 0.000) 0.000 6.355 ( 0.000 0.000 0.000) 0.000 6.364 ( -0.000 -0.000 0.000) 0.000 6.393 ( 0.000 0.000 0.000) 0.000 6.842 ( 0.000 0.000 0.000) 0.000 6.886 ( 0.000 0.000 0.000) 0.000 8.015 ( 0.000 0.000 0.000) 0.000 8.097 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 16 (9/48) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.750 ( 10.220 17.702 0.000) 20.441 0.792 ( 10.119 17.527 0.000) 20.238 1.270 ( 12.104 20.964 0.000) 24.208 1.415 ( -0.983 -1.703 0.000) 1.966 1.416 ( 7.350 12.731 0.000) 14.701 1.486 ( 3.679 6.372 0.000) 7.358 1.584 ( 8.244 14.279 0.000) 16.487 1.858 ( 2.335 4.044 0.000) 4.669 2.038 ( 7.426 12.862 0.000) 14.851 2.299 ( 3.926 6.799 0.000) 7.851 2.624 ( 1.920 3.326 0.000) 3.840 3.020 ( -1.830 -3.170 0.000) 3.661 5.089 ( 4.295 7.439 0.000) 8.590 5.100 ( 3.976 6.886 0.000) 7.952 5.984 ( 17.471 30.260 0.000) 34.942 6.178 ( -3.865 -6.695 0.000) 7.731 6.227 ( -3.060 -5.299 0.000) 6.119 6.287 ( -1.093 -1.893 0.000) 2.186 6.301 ( -1.268 -2.196 0.000) 2.536 6.815 ( 0.668 1.158 0.000) 1.337 7.831 ( -5.672 -9.824 0.000) 11.344 8.031 ( -8.872 -15.367 0.000) 17.744 8.132 ( -1.992 -3.449 0.000) 3.983 8.276 ( 0.349 0.605 0.000) 0.699 ======================= Grid point 17 (10/48) ======================= q-point: ( 0.14 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.107 ( 6.386 13.998 0.000) 15.386 1.118 ( 12.344 13.914 0.000) 18.601 1.393 ( 0.998 -0.857 0.000) 1.315 1.481 ( 4.093 2.625 0.000) 4.862 1.688 ( 7.547 12.930 0.000) 14.971 1.754 ( 9.292 11.818 0.000) 15.033 1.818 ( 5.936 5.273 0.000) 7.940 1.959 ( 2.854 6.595 0.000) 7.186 2.413 ( 16.897 16.896 0.000) 23.895 2.482 ( 7.486 7.184 0.000) 10.376 2.795 ( 8.587 9.763 0.000) 13.002 2.940 ( -2.338 -3.682 0.000) 4.362 5.266 ( 5.033 8.784 0.000) 10.124 5.266 ( 5.189 9.562 0.000) 10.880 6.027 ( -5.028 -6.275 0.000) 8.041 6.109 ( -3.693 -5.170 0.000) 6.354 6.239 ( -1.700 -1.240 0.000) 2.105 6.252 ( -1.852 -1.474 0.000) 2.367 6.476 ( 10.360 12.732 0.000) 16.415 6.845 ( 1.972 0.566 0.000) 2.052 7.590 ( -10.334 -10.087 0.000) 14.440 7.686 ( -11.972 -12.978 0.000) 17.656 8.080 ( -0.320 -1.650 0.000) 1.681 8.274 ( 0.717 -2.875 0.000) 2.963 ======================= Grid point 18 (11/48) ======================= q-point: ( 0.21 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.306 ( 4.023 8.363 0.000) 9.280 1.448 ( 3.380 2.812 0.000) 4.396 1.451 ( 10.847 10.360 0.000) 14.999 1.581 ( 3.108 5.438 0.000) 6.263 1.897 ( 2.988 3.216 0.000) 4.389 1.954 ( 5.246 9.865 0.000) 11.173 1.971 ( 6.169 4.746 0.000) 7.783 2.091 ( 3.172 8.442 0.000) 9.018 2.664 ( 5.877 3.852 0.000) 7.027 2.796 ( -5.830 -3.393 0.000) 6.745 2.922 ( 20.065 15.517 0.000) 25.365 3.150 ( 14.609 10.173 0.000) 17.802 5.435 ( 2.859 9.755 0.000) 10.166 5.443 ( 2.639 10.444 0.000) 10.772 5.896 ( -3.076 -4.092 0.000) 5.120 6.008 ( -1.850 -4.220 0.000) 4.608 6.203 ( -1.405 -0.176 0.000) 1.416 6.215 ( -1.529 0.208 0.000) 1.543 6.624 ( -0.248 2.433 0.000) 2.445 6.890 ( 2.299 -0.555 0.000) 2.365 7.267 ( -13.579 -8.888 0.000) 16.229 7.378 ( -10.118 -9.912 0.000) 14.163 8.084 ( 2.056 -1.879 0.000) 2.785 8.245 ( 0.766 -5.306 0.000) 5.361 ======================= Grid point 19 (12/48) ======================= q-point: ( 0.29 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.448 ( 3.576 6.337 0.000) 7.276 1.534 ( 2.671 2.678 0.000) 3.783 1.703 ( 7.454 6.851 0.000) 10.124 1.711 ( 3.757 7.978 0.000) 8.818 1.937 ( 1.244 -0.007 0.000) 1.244 2.089 ( 4.117 1.228 0.000) 4.296 2.150 ( 3.763 10.781 0.000) 11.419 2.232 ( 2.219 9.239 0.000) 9.502 2.646 ( -6.966 -2.970 0.000) 7.573 2.745 ( 0.884 -0.186 0.000) 0.903 3.407 ( 18.203 11.824 0.000) 21.706 3.506 ( 13.983 8.505 0.000) 16.367 5.559 ( -0.225 10.140 0.000) 10.142 5.567 ( -0.571 10.451 0.000) 10.467 5.851 ( 1.045 -0.920 0.000) 1.392 5.972 ( 2.202 -2.794 0.000) 3.557 6.183 ( -1.071 0.069 0.000) 1.074 6.201 ( -1.515 2.086 0.000) 2.579 6.589 ( -3.627 -0.540 0.000) 3.666 6.903 ( 0.416 -1.087 0.000) 1.164 6.959 ( -11.823 -7.353 0.000) 13.923 7.121 ( -8.795 -6.865 0.000) 11.157 8.087 ( 1.215 -4.025 0.000) 4.205 8.193 ( 0.129 -6.171 0.000) 6.172 ======================= Grid point 20 (13/48) ======================= q-point: ( 0.36 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.553 ( 0.821 6.267 0.000) 6.320 1.593 ( 1.149 2.320 0.000) 2.589 1.853 ( 3.413 3.661 0.000) 5.005 1.866 ( 4.731 6.413 0.000) 7.969 1.956 ( 1.270 -0.257 0.000) 1.296 2.154 ( 3.155 -0.533 0.000) 3.200 2.298 ( 1.525 8.623 0.000) 8.757 2.353 ( 0.481 9.517 0.000) 9.530 2.505 ( -6.636 -0.177 0.000) 6.638 2.697 ( -4.008 -2.687 0.000) 4.825 3.787 ( 13.224 7.088 0.000) 15.004 3.797 ( 10.349 5.250 0.000) 11.604 5.641 ( -2.201 10.357 0.000) 10.588 5.645 ( -2.218 10.087 0.000) 10.328 5.908 ( 4.559 1.392 0.000) 4.767 6.027 ( 6.039 -1.932 0.000) 6.341 6.159 ( -1.365 -0.571 0.000) 1.479 6.201 ( -1.531 3.766 0.000) 4.065 6.506 ( -4.571 -0.862 0.000) 4.652 6.692 ( -10.062 -6.920 0.000) 12.212 6.893 ( -0.193 -0.992 0.000) 1.011 6.940 ( -5.519 -3.432 0.000) 6.499 8.046 ( -0.430 -4.675 0.000) 4.695 8.128 ( -0.541 -5.352 0.000) 5.379 ======================= Grid point 21 (14/48) ======================= q-point: ( 0.43 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.602 ( -1.986 5.509 0.000) 5.856 1.625 ( -0.141 2.087 0.000) 2.092 1.904 ( 0.071 1.335 0.000) 1.337 1.976 ( 1.504 3.078 0.000) 3.426 1.980 ( 1.435 -0.487 0.000) 1.516 2.201 ( 3.122 -1.625 0.000) 3.520 2.369 ( -0.005 1.957 0.000) 1.957 2.429 ( -5.751 7.567 0.000) 9.505 2.443 ( -1.714 10.096 0.000) 10.240 2.581 ( -4.457 -2.890 0.000) 5.311 3.963 ( 4.341 0.974 0.000) 4.449 4.000 ( 5.409 1.466 0.000) 5.604 5.691 ( -4.094 9.883 0.000) 10.697 5.693 ( -3.713 9.125 0.000) 9.851 6.004 ( 3.734 0.503 0.000) 3.768 6.091 ( 2.271 -1.818 0.000) 2.909 6.131 ( -0.756 -0.492 0.000) 0.902 6.269 ( 5.352 4.420 0.000) 6.941 6.425 ( -3.456 0.402 0.000) 3.479 6.462 ( -8.736 -5.294 0.000) 10.215 6.853 ( -1.817 -1.098 0.000) 2.123 6.882 ( 0.043 -0.736 0.000) 0.737 7.997 ( 0.171 -3.287 0.000) 3.292 8.078 ( 0.338 -3.459 0.000) 3.475 ======================= Grid point 31 (15/48) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.304 ( 2.653 4.596 0.000) 5.306 1.406 ( 7.906 13.694 0.000) 15.812 1.481 ( 5.422 9.391 0.000) 10.844 1.549 ( 2.444 4.232 0.000) 4.887 1.891 ( 1.362 2.359 0.000) 2.724 1.943 ( 4.005 6.936 0.000) 8.009 1.961 ( 7.384 12.790 0.000) 14.768 2.116 ( 4.849 8.398 0.000) 9.697 2.626 ( 3.700 6.408 0.000) 7.400 2.812 ( 12.184 21.103 0.000) 24.367 2.840 ( -2.972 -5.148 0.000) 5.945 3.058 ( 8.410 14.566 0.000) 16.819 5.469 ( 5.980 10.357 0.000) 11.960 5.484 ( 6.283 10.882 0.000) 12.565 5.902 ( -3.072 -5.322 0.000) 6.145 6.001 ( -2.870 -4.972 0.000) 5.741 6.218 ( -0.502 -0.869 0.000) 1.004 6.231 ( -0.256 -0.444 0.000) 0.513 6.648 ( 2.212 3.832 0.000) 4.424 6.853 ( 0.325 0.562 0.000) 0.649 7.357 ( -7.443 -12.891 0.000) 14.885 7.416 ( -7.199 -12.468 0.000) 14.397 8.050 ( -0.818 -1.416 0.000) 1.635 8.195 ( -2.523 -4.370 0.000) 5.046 ======================= Grid point 32 (16/48) ======================= q-point: ( 0.21 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.403 ( 4.623 3.835 0.000) 6.006 1.605 ( -0.849 8.738 0.000) 8.779 1.658 ( 6.811 9.452 0.000) 11.650 1.705 ( 6.342 6.909 0.000) 9.378 1.923 ( 0.783 0.736 0.000) 1.074 2.049 ( 2.625 2.805 0.000) 3.842 2.197 ( 4.824 11.518 0.000) 12.487 2.291 ( 3.595 10.529 0.000) 11.126 2.725 ( 2.086 2.250 0.000) 3.068 2.729 ( -4.662 -2.473 0.000) 5.277 3.244 ( 14.672 15.547 0.000) 21.377 3.372 ( 11.335 11.176 0.000) 15.918 5.662 ( 2.736 11.374 0.000) 11.699 5.681 ( 2.361 11.602 0.000) 11.840 5.829 ( -0.449 -1.953 0.000) 2.004 5.922 ( -0.589 -3.732 0.000) 3.778 6.202 ( -0.844 -0.189 0.000) 0.865 6.239 ( 0.041 2.193 0.000) 2.193 6.656 ( -2.868 0.935 0.000) 3.016 6.870 ( 1.518 -1.127 0.000) 1.891 7.072 ( -10.197 -10.829 0.000) 14.875 7.191 ( -7.025 -7.427 0.000) 10.223 8.018 ( 0.435 -4.301 0.000) 4.323 8.110 ( -0.183 -7.161 0.000) 7.164 ======================= Grid point 33 (17/48) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 200 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.537 ( 4.988 4.453 0.000) 6.687 1.669 ( -1.356 7.743 0.000) 7.861 1.832 ( 4.136 5.505 0.000) 6.885 1.866 ( 3.855 6.984 0.000) 7.977 1.940 ( 0.578 0.254 0.000) 0.631 2.104 ( 2.106 0.245 0.000) 2.121 2.350 ( 1.615 6.221 0.000) 6.427 2.442 ( 1.409 9.487 0.000) 9.591 2.649 ( -5.264 5.513 0.000) 7.623 2.743 ( -1.939 1.771 0.000) 2.626 3.625 ( 13.513 9.131 0.000) 16.308 3.663 ( 10.757 6.448 0.000) 12.542 5.793 ( -0.734 10.844 0.000) 10.869 5.795 ( -0.749 8.594 0.000) 8.627 5.859 ( 2.406 3.642 0.000) 4.366 5.919 ( 3.529 -1.694 0.000) 3.915 6.177 ( -1.407 -0.921 0.000) 1.681 6.266 ( -0.964 3.943 0.000) 4.059 6.597 ( -4.937 1.918 0.000) 5.296 6.794 ( -9.659 -8.278 0.000) 12.721 6.876 ( 0.718 -1.537 0.000) 1.696 7.013 ( -6.657 -3.452 0.000) 7.499 7.978 ( 1.105 -5.949 0.000) 6.051 8.047 ( 1.514 -7.471 0.000) 7.623 ======================= Grid point 34 (18/48) ======================= q-point: ( 0.36 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.646 ( 1.665 4.627 0.000) 4.917 1.712 ( -2.310 6.922 0.000) 7.297 1.917 ( 0.816 2.378 0.000) 2.514 1.954 ( 1.063 0.069 0.000) 1.065 1.984 ( 2.267 4.819 0.000) 5.326 2.139 ( 2.383 -0.875 0.000) 2.539 2.381 ( -0.540 0.066 0.000) 0.544 2.512 ( 0.024 1.078 0.000) 1.079 2.661 ( -6.073 14.426 0.000) 15.652 2.722 ( -6.262 9.936 0.000) 11.745 3.866 ( 7.413 1.113 0.000) 7.496 3.888 ( 9.018 2.427 0.000) 9.339 5.860 ( -1.800 8.447 0.000) 8.636 5.868 ( -2.213 9.377 0.000) 9.635 5.947 ( 2.828 3.076 0.000) 4.179 5.996 ( 5.314 -0.462 0.000) 5.334 6.139 ( -1.115 -1.347 0.000) 1.748 6.277 ( -1.164 2.402 0.000) 2.669 6.528 ( -5.805 3.662 0.000) 6.863 6.575 ( -8.549 -2.235 0.000) 8.836 6.869 ( 0.282 -1.297 0.000) 1.328 6.889 ( -4.185 -1.578 0.000) 4.473 7.941 ( 1.226 -5.061 0.000) 5.207 8.009 ( 1.954 -5.888 0.000) 6.204 ======================= Grid point 35 (19/48) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.685 ( -2.084 3.610 0.000) 4.168 1.728 ( -3.535 6.123 0.000) 7.071 1.930 ( -0.666 1.154 0.000) 1.333 1.968 ( 0.355 -0.614 0.000) 0.709 2.030 ( -1.169 2.025 0.000) 2.338 2.162 ( 1.228 -2.127 0.000) 2.456 2.373 ( -0.074 0.128 0.000) 0.148 2.503 ( 1.600 -2.771 0.000) 3.200 2.686 ( -8.680 15.035 0.000) 17.361 2.723 ( -8.597 14.890 0.000) 17.193 3.939 ( 1.980 -3.429 0.000) 3.960 3.982 ( 1.891 -3.275 0.000) 3.781 5.890 ( -4.896 8.481 0.000) 9.793 5.900 ( -4.825 8.358 0.000) 9.651 5.992 ( 0.788 -1.364 0.000) 1.575 6.049 ( 1.574 -2.726 0.000) 3.148 6.123 ( 0.220 -0.381 0.000) 0.440 6.284 ( 1.710 -2.962 0.000) 3.420 6.475 ( -4.675 8.097 0.000) 9.349 6.489 ( -3.581 6.203 0.000) 7.163 6.843 ( -0.111 0.192 0.000) 0.222 6.864 ( 0.546 -0.946 0.000) 1.093 7.925 ( 2.113 -3.660 0.000) 4.227 7.996 ( 2.604 -4.509 0.000) 5.207 ======================= Grid point 47 (20/48) ======================= q-point: ( 0.21 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.511 ( 3.840 6.651 0.000) 7.679 1.735 ( 2.578 4.465 0.000) 5.155 1.823 ( 3.715 6.435 0.000) 7.431 1.860 ( 4.286 7.424 0.000) 8.572 1.937 ( 0.354 0.614 0.000) 0.708 2.087 ( 0.649 1.125 0.000) 1.299 2.367 ( 2.360 4.088 0.000) 4.720 2.478 ( 3.510 6.079 0.000) 7.019 2.738 ( 2.501 4.331 0.000) 5.001 2.780 ( 2.495 4.322 0.000) 4.990 3.541 ( 7.561 13.097 0.000) 15.123 3.590 ( 5.600 9.700 0.000) 11.200 5.811 ( 0.887 1.536 0.000) 1.773 5.848 ( 2.064 3.575 0.000) 4.128 5.885 ( 1.865 3.230 0.000) 3.729 5.894 ( 4.299 7.446 0.000) 8.598 6.186 ( -0.896 -1.552 0.000) 1.792 6.298 ( 1.734 3.003 0.000) 3.468 6.673 ( 0.426 0.738 0.000) 0.852 6.840 ( -6.203 -10.744 0.000) 12.406 6.850 ( -0.415 -0.719 0.000) 0.831 7.070 ( -2.911 -5.042 0.000) 5.822 7.928 ( -2.271 -3.933 0.000) 4.541 7.980 ( -2.916 -5.051 0.000) 5.833 ======================= Grid point 48 (21/48) ======================= q-point: ( 0.29 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.656 ( 3.845 6.907 0.000) 7.905 1.803 ( -0.270 5.089 0.000) 5.096 1.919 ( 1.798 2.820 0.000) 3.344 1.947 ( 0.196 0.474 0.000) 0.513 1.988 ( 2.891 4.748 0.000) 5.558 2.104 ( 0.849 -0.149 0.000) 0.862 2.389 ( -0.411 -0.319 0.000) 0.520 2.509 ( -0.967 -2.324 0.000) 2.518 2.877 ( 0.409 14.338 0.000) 14.344 2.915 ( 0.705 13.955 0.000) 13.973 3.752 ( 5.708 1.811 0.000) 5.988 3.763 ( 7.270 3.879 0.000) 8.240 5.867 ( 2.842 2.793 0.000) 3.984 5.882 ( 2.802 0.732 0.000) 2.896 5.985 ( 1.195 5.141 0.000) 5.278 6.004 ( 0.432 5.701 0.000) 5.717 6.142 ( -1.364 -2.304 0.000) 2.677 6.317 ( -1.893 -1.066 0.000) 2.173 6.676 ( -3.726 5.202 0.000) 6.399 6.683 ( -3.736 -0.453 0.000) 3.763 6.843 ( 0.864 -1.645 0.000) 1.858 6.968 ( -5.041 -1.446 0.000) 5.244 7.871 ( 0.447 -4.035 0.000) 4.060 7.912 ( 1.124 -5.257 0.000) 5.376 ======================= Grid point 49 (22/48) ======================= q-point: ( 0.36 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.760 ( -0.304 6.165 0.000) 6.172 1.834 ( -1.799 4.723 0.000) 5.054 1.949 ( -0.268 0.860 0.000) 0.901 1.957 ( 0.585 0.203 0.000) 0.619 2.055 ( 0.681 2.002 0.000) 2.115 2.120 ( 1.447 -0.914 0.000) 1.711 2.378 ( -0.358 0.051 0.000) 0.362 2.453 ( -0.220 -4.130 0.000) 4.136 2.996 ( -6.142 17.222 0.000) 18.284 3.036 ( -6.291 17.674 0.000) 18.761 3.828 ( 4.925 -5.103 0.000) 7.092 3.873 ( 5.428 -4.145 0.000) 6.829 5.952 ( 2.892 2.573 0.000) 3.871 5.962 ( 4.168 1.406 0.000) 4.399 6.020 ( -1.379 1.459 0.000) 2.008 6.022 ( -2.878 2.270 0.000) 3.665 6.111 ( 0.232 -1.339 0.000) 1.359 6.237 ( -1.257 -5.001 0.000) 5.156 6.656 ( -5.423 7.755 0.000) 9.463 6.683 ( -4.091 10.426 0.000) 11.200 6.843 ( 0.704 -1.133 0.000) 1.334 6.871 ( -3.912 -0.324 0.000) 3.925 7.854 ( 1.562 -3.169 0.000) 3.534 7.898 ( 2.647 -4.828 0.000) 5.506 ======================= Grid point 62 (23/48) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.799 ( 3.721 6.444 0.000) 7.441 1.876 ( 1.240 2.147 0.000) 2.479 1.952 ( 0.445 0.771 0.000) 0.890 1.956 ( 0.170 0.294 0.000) 0.339 2.055 ( 0.969 1.678 0.000) 1.938 2.102 ( 0.036 0.063 0.000) 0.072 2.382 ( -0.063 -0.109 0.000) 0.126 2.432 ( -2.253 -3.902 0.000) 4.506 3.168 ( 7.483 12.960 0.000) 14.965 3.214 ( 7.852 13.601 0.000) 15.705 3.730 ( -2.091 -3.622 0.000) 4.182 3.776 ( -1.386 -2.400 0.000) 2.772 5.923 ( 1.804 3.124 0.000) 3.608 5.937 ( 2.196 3.804 0.000) 4.392 6.043 ( 0.558 0.967 0.000) 1.116 6.065 ( 0.345 0.598 0.000) 0.690 6.094 ( -1.176 -2.037 0.000) 2.352 6.222 ( -3.557 -6.161 0.000) 7.114 6.761 ( 3.341 5.786 0.000) 6.681 6.770 ( 1.709 2.961 0.000) 3.419 6.816 ( -0.287 -0.496 0.000) 0.573 6.936 ( -1.209 -2.094 0.000) 2.417 7.821 ( -0.605 -1.048 0.000) 1.210 7.839 ( -1.122 -1.943 0.000) 2.243 ======================= Grid point 63 (24/48) ======================= q-point: ( 0.36 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.875 ( -2.408 4.170 0.000) 4.815 1.899 ( -0.881 1.525 0.000) 1.761 1.959 ( -0.116 0.201 0.000) 0.232 1.960 ( -0.018 0.031 0.000) 0.036 2.064 ( 0.420 -0.727 0.000) 0.840 2.104 ( 0.254 -0.440 0.000) 0.508 2.382 ( -0.102 0.177 0.000) 0.204 2.386 ( 1.063 -1.842 0.000) 2.127 3.332 ( -8.352 14.466 0.000) 16.704 3.386 ( -8.677 15.029 0.000) 17.354 3.663 ( 5.847 -10.127 0.000) 11.694 3.721 ( 5.720 -9.907 0.000) 11.439 5.974 ( 0.094 -0.163 0.000) 0.188 5.996 ( -0.867 1.502 0.000) 1.734 6.045 ( -0.538 0.931 0.000) 1.075 6.049 ( -1.231 2.132 0.000) 2.462 6.079 ( 1.070 -1.854 0.000) 2.141 6.138 ( 2.320 -4.018 0.000) 4.640 6.787 ( -2.316 4.012 0.000) 4.632 6.824 ( 0.447 -0.774 0.000) 0.893 6.869 ( -2.624 4.545 0.000) 5.248 6.871 ( -0.687 1.190 0.000) 1.374 7.813 ( 0.545 -0.944 0.000) 1.090 7.822 ( 1.267 -2.195 0.000) 2.534 ======================= Grid point 211 (25/48) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.450 ( -0.000 0.000 20.769) 20.769 0.450 ( 0.000 -0.000 20.769) 20.769 0.859 ( 0.000 0.000 -17.424) 17.424 0.859 ( 0.000 0.000 -17.424) 17.424 0.861 ( -0.000 0.000 40.058) 40.058 1.573 ( 0.000 0.000 9.255) 9.255 1.573 ( -0.000 0.000 9.255) 9.255 1.665 ( 0.000 0.000 -35.549) 35.549 1.729 ( -0.000 -0.000 -5.171) 5.171 1.729 ( -0.000 -0.000 -5.171) 5.171 3.086 ( -0.000 0.000 22.259) 22.259 3.465 ( 0.000 -0.000 -12.808) 12.808 4.936 ( 0.000 0.000 0.478) 0.478 4.936 ( -0.000 0.000 0.478) 0.478 4.947 ( 0.000 0.000 -0.487) 0.487 4.947 ( 0.000 -0.000 -0.487) 0.487 6.341 ( -0.000 0.000 0.405) 0.405 6.341 ( -0.000 -0.000 0.405) 0.405 6.350 ( -0.000 0.000 -0.398) 0.398 6.350 ( 0.000 0.000 -0.398) 0.398 8.140 ( 0.000 -0.000 2.493) 2.493 8.201 ( -0.000 0.000 -2.748) 2.748 8.418 ( 0.000 0.000 0.636) 0.636 8.429 ( 0.000 0.000 -0.323) 0.323 ======================= Grid point 212 (26/48) ======================= q-point: ( 0.07 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.647 ( 13.254 7.652 14.421) 21.029 0.692 ( 12.149 7.014 16.441) 21.612 0.978 ( 8.866 5.119 -15.696) 18.740 1.041 ( 10.772 6.219 -16.300) 20.504 1.117 ( 20.873 12.051 26.934) 36.144 1.598 ( 2.189 1.264 9.279) 9.618 1.639 ( 6.766 3.907 8.369) 11.450 1.751 ( 1.907 1.101 -5.021) 5.482 1.758 ( 7.788 4.496 -28.624) 30.004 1.811 ( 6.855 3.958 -5.449) 9.610 2.929 ( -7.462 -4.308 15.158) 17.435 3.106 ( -9.637 -5.564 -0.998) 11.173 4.988 ( 4.293 2.478 0.405) 4.973 4.998 ( 4.155 2.399 -0.412) 4.815 5.535 ( 37.403 21.595 12.789) 45.043 6.039 ( 39.488 22.798 -37.298) 58.908 6.293 ( -3.189 -1.841 -0.118) 3.684 6.302 ( -3.774 -2.179 -0.501) 4.387 6.321 ( -1.566 -0.904 0.391) 1.850 6.330 ( -1.592 -0.919 -0.385) 1.878 8.074 ( -7.344 -4.240 4.024) 9.387 8.150 ( -5.135 -2.965 -3.196) 6.736 8.280 ( -8.058 -4.652 1.569) 9.436 8.291 ( -6.575 -3.796 0.499) 7.609 ======================= Grid point 213 (27/48) ======================= q-point: ( 0.14 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.962 ( 12.310 7.107 9.385) 17.033 0.972 ( 12.285 7.093 9.528) 17.088 1.205 ( 9.069 5.236 -12.678) 16.443 1.222 ( 5.348 3.088 -13.066) 14.452 1.630 ( 18.033 10.411 14.427) 25.332 1.682 ( 5.069 2.926 8.534) 10.348 1.820 ( 3.916 2.261 -4.454) 6.347 1.832 ( 6.920 3.995 5.874) 9.917 2.064 ( 17.400 10.046 -9.376) 22.172 2.088 ( 17.917 10.344 -11.626) 23.731 2.844 ( 0.382 0.220 9.635) 9.645 2.976 ( -3.115 -1.798 -2.981) 4.672 5.126 ( 6.941 4.008 0.227) 8.018 5.131 ( 6.743 3.893 -0.230) 7.789 6.150 ( -5.836 -3.370 1.469) 6.897 6.177 ( -4.050 -2.338 -0.933) 4.769 6.272 ( -2.409 -1.391 0.345) 2.803 6.280 ( -2.473 -1.428 -0.343) 2.876 6.301 ( 21.440 12.379 11.070) 27.120 6.614 ( 12.439 7.182 -16.319) 21.740 7.815 ( -12.827 -7.406 3.830) 15.298 7.889 ( -14.868 -8.584 -2.754) 17.388 8.153 ( -2.270 -1.311 3.331) 4.239 8.240 ( -0.105 -0.061 -4.441) 4.443 ======================= Grid point 214 (28/48) ======================= q-point: ( 0.21 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.208 ( 8.498 4.906 6.550) 11.798 1.270 ( 12.598 7.273 2.337) 14.733 1.359 ( 6.902 3.985 -5.538) 9.705 1.378 ( 5.577 3.220 -8.608) 10.750 1.834 ( 7.451 4.302 5.708) 10.325 1.913 ( 7.077 4.086 4.911) 9.534 1.933 ( 5.542 3.200 -3.371) 7.233 1.950 ( 3.594 2.075 2.826) 5.021 2.507 ( 14.564 8.408 -3.305) 17.138 2.531 ( 17.836 10.298 -5.137) 21.226 2.924 ( 4.149 2.395 4.454) 6.541 2.978 ( 7.576 4.374 -0.667) 8.774 5.293 ( 6.687 3.861 0.047) 7.722 5.294 ( 6.548 3.781 -0.048) 7.561 6.004 ( -5.639 -3.256 2.237) 6.885 6.055 ( -4.788 -2.764 -2.120) 5.921 6.218 ( -1.936 -1.118 0.242) 2.248 6.224 ( -2.082 -1.202 -0.241) 2.416 6.627 ( 6.027 3.480 6.403) 9.457 6.787 ( 3.614 2.086 -7.575) 8.649 7.490 ( -14.346 -8.283 2.936) 16.823 7.543 ( -13.644 -7.877 -1.800) 15.857 8.141 ( 0.860 0.496 4.164) 4.280 8.246 ( 0.196 0.113 -4.853) 4.858 ======================= Grid point 215 (29/48) ======================= q-point: ( 0.29 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.372 ( 5.619 3.244 4.341) 7.807 1.480 ( 3.311 1.912 -5.029) 6.317 1.541 ( 9.392 5.422 -0.200) 10.847 1.555 ( 9.936 5.736 -1.133) 11.529 1.984 ( 0.999 0.577 3.727) 3.902 2.014 ( 7.426 4.287 2.938) 9.064 2.044 ( 3.171 1.831 -1.302) 3.886 2.071 ( 5.931 3.424 -2.140) 7.175 2.699 ( -1.402 -0.810 -0.995) 1.901 2.706 ( 3.274 1.890 -1.330) 4.007 3.211 ( 18.292 10.561 4.210) 21.537 3.293 ( 15.962 9.215 -2.993) 18.672 5.425 ( 4.512 2.605 0.044) 5.211 5.426 ( 4.563 2.635 -0.044) 5.270 5.909 ( -2.020 -1.166 2.844) 3.678 5.976 ( -1.496 -0.864 -2.939) 3.409 6.186 ( -0.909 -0.525 0.048) 1.050 6.187 ( -1.129 -0.652 -0.048) 1.304 6.666 ( -1.340 -0.774 5.562) 5.774 6.811 ( -1.319 -0.761 -7.352) 7.508 7.175 ( -11.663 -6.733 4.119) 14.083 7.244 ( -11.617 -6.707 -2.210) 13.595 8.161 ( 0.181 0.104 3.172) 3.179 8.236 ( -1.263 -0.729 -3.236) 3.549 ======================= Grid point 216 (30/48) ======================= q-point: ( 0.36 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.481 ( 3.668 2.118 2.571) 4.955 1.543 ( 2.115 1.221 -2.751) 3.679 1.746 ( 7.918 4.572 0.610) 9.163 1.766 ( 7.259 4.191 -0.955) 8.436 2.000 ( 0.625 0.361 3.781) 3.849 2.094 ( 1.457 0.841 -4.475) 4.781 2.173 ( 5.863 3.385 1.412) 6.916 2.203 ( 5.000 2.887 -1.196) 5.896 2.616 ( -4.615 -2.665 4.457) 6.947 2.705 ( -2.696 -1.557 -3.381) 4.596 3.626 ( 15.703 9.066 1.165) 18.169 3.653 ( 13.675 7.895 -1.126) 15.830 5.512 ( 3.000 1.732 0.097) 3.466 5.514 ( 3.071 1.773 -0.097) 3.547 5.919 ( 2.783 1.607 3.320) 4.621 5.996 ( 3.277 1.892 -3.400) 5.088 6.163 ( -1.094 -0.632 0.161) 1.274 6.166 ( -0.926 -0.534 -0.160) 1.081 6.596 ( -4.338 -2.504 4.938) 7.034 6.724 ( -6.236 -3.600 -6.621) 9.782 6.967 ( -5.548 -3.203 2.965) 7.060 7.012 ( -7.568 -4.369 -1.369) 8.846 8.136 ( -2.241 -1.294 2.194) 3.393 8.185 ( -2.900 -1.674 -2.022) 3.911 ======================= Grid point 217 (31/48) ======================= q-point: ( 0.43 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.546 ( 1.783 1.029 1.533) 2.567 1.582 ( 1.195 0.690 -1.582) 2.099 1.882 ( 3.076 1.776 0.617) 3.605 1.901 ( 4.152 2.397 -0.930) 4.884 2.025 ( 1.664 0.961 4.014) 4.450 2.132 ( 1.960 1.132 -5.604) 6.044 2.283 ( 3.300 1.905 0.711) 3.876 2.299 ( 2.925 1.689 -0.649) 3.440 2.507 ( -4.004 -2.312 5.540) 7.216 2.608 ( -4.745 -2.740 -3.522) 6.514 3.912 ( 7.871 4.544 0.440) 9.099 3.923 ( 8.928 5.155 -0.435) 10.319 5.571 ( 1.969 1.137 0.188) 2.282 5.576 ( 2.067 1.194 -0.192) 2.395 6.023 ( 5.346 3.086 4.400) 7.581 6.121 ( 6.820 3.937 -3.998) 8.832 6.133 ( -1.230 -0.710 0.276) 1.446 6.140 ( -1.174 -0.678 -0.271) 1.383 6.475 ( -5.417 -3.128 2.424) 6.708 6.540 ( -8.426 -4.865 -3.283) 10.269 6.889 ( -2.948 -1.702 0.239) 3.413 6.893 ( -1.703 -0.983 -0.079) 1.968 8.073 ( -2.557 -1.476 1.931) 3.528 8.114 ( -2.642 -1.525 -1.634) 3.461 ======================= Grid point 218 (32/48) ======================= q-point: (-0.50 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 89 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.566 ( 0.000 0.000 1.300) 1.300 1.597 ( 0.000 0.000 -1.326) 1.326 1.913 ( 0.000 0.000 1.391) 1.391 1.947 ( 0.000 0.000 -1.652) 1.652 2.055 ( 0.000 0.000 4.253) 4.253 2.166 ( 0.000 0.000 -5.989) 5.989 2.323 ( 0.000 0.000 0.513) 0.513 2.334 ( 0.000 0.000 -0.476) 0.476 2.456 ( 0.000 0.000 4.792) 4.792 2.536 ( 0.000 0.000 -2.644) 2.644 4.008 ( 0.000 0.000 0.988) 0.988 4.031 ( 0.000 0.000 -0.966) 0.966 5.597 ( 0.000 0.000 0.247) 0.247 5.602 ( 0.000 0.000 -0.253) 0.253 6.109 ( 0.000 0.000 6.161) 6.161 6.116 ( 0.000 0.000 0.299) 0.299 6.123 ( 0.000 0.000 -0.292) 0.292 6.260 ( 0.000 0.000 -7.030) 7.030 6.376 ( -0.000 -0.000 1.050) 1.050 6.398 ( -0.000 -0.000 -0.610) 0.610 6.854 ( 0.000 0.000 0.974) 0.974 6.876 ( 0.000 0.000 -0.895) 0.895 8.039 ( 0.000 0.000 1.940) 1.940 8.080 ( 0.000 0.000 -1.595) 1.595 ======================= Grid point 227 (33/48) ======================= q-point: ( 0.07 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.884 ( 6.493 11.247 11.184) 17.139 0.907 ( 8.535 14.783 9.751) 19.658 1.167 ( 6.517 11.287 -12.768) 18.246 1.175 ( 3.200 5.542 -14.326) 15.690 1.484 ( 10.786 18.681 17.188) 27.582 1.677 ( 4.100 7.101 9.015) 12.186 1.763 ( 3.346 5.796 5.176) 8.461 1.814 ( 3.405 5.897 -1.738) 7.028 1.964 ( 9.298 16.104 -20.138) 27.411 1.970 ( 7.241 12.542 -6.062) 15.700 2.847 ( -1.103 -1.911 11.258) 11.472 2.998 ( -2.374 -4.112 -3.199) 5.726 5.091 ( 4.219 7.307 0.242) 8.441 5.097 ( 4.059 7.030 -0.242) 8.121 6.116 ( 15.070 26.102 11.710) 32.335 6.189 ( -3.433 -5.945 0.778) 6.909 6.225 ( -1.881 -3.259 0.602) 3.811 6.290 ( -1.126 -1.950 0.319) 2.275 6.298 ( -1.229 -2.129 -0.308) 2.477 6.519 ( 9.620 16.662 -21.326) 28.723 7.897 ( -6.920 -11.986 4.550) 14.569 7.985 ( -7.765 -13.449 -3.604) 15.942 8.168 ( -2.319 -4.017 2.740) 5.388 8.233 ( -0.827 -1.432 -3.175) 3.580 ======================= Grid point 228 (34/48) ======================= q-point: ( 0.14 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.134 ( 8.581 9.747 3.087) 13.348 1.178 ( 7.306 9.475 3.793) 12.552 1.293 ( 4.453 6.775 -4.853) 9.448 1.381 ( 4.897 9.512 -7.936) 13.320 1.751 ( 7.588 6.404 6.825) 12.049 1.853 ( 3.490 5.402 2.144) 6.779 1.890 ( 6.004 6.909 3.680) 9.865 1.947 ( 3.081 6.779 0.469) 7.461 2.347 ( 16.000 16.804 -4.618) 23.658 2.358 ( 14.258 14.609 -6.313) 21.368 2.881 ( 3.926 3.932 4.809) 7.348 2.940 ( -0.336 -1.341 -0.783) 1.589 5.266 ( 5.070 8.987 0.019) 10.319 5.266 ( 5.148 9.377 -0.007) 10.697 6.049 ( -4.688 -5.998 1.853) 7.835 6.090 ( -4.002 -5.413 -1.695) 6.942 6.244 ( -1.830 -1.273 0.320) 2.252 6.250 ( -1.870 -1.378 -0.230) 2.334 6.558 ( 8.782 10.388 7.267) 15.422 6.740 ( 4.613 4.466 -8.543) 10.687 7.621 ( -10.604 -11.036 2.400) 15.492 7.668 ( -11.432 -12.375 -1.703) 16.934 8.120 ( 0.065 -1.753 3.659) 4.058 8.216 ( 0.691 -2.269 -4.646) 5.216 ======================= Grid point 229 (35/48) ======================= q-point: ( 0.21 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.332 ( 5.414 6.018 2.156) 8.377 1.391 ( 7.480 4.473 -3.337) 9.332 1.475 ( 5.918 10.522 2.167) 12.266 1.546 ( 3.526 9.224 -3.450) 10.460 1.932 ( 4.535 3.263 1.939) 5.914 1.942 ( 4.721 1.913 0.952) 5.182 2.041 ( 2.984 7.414 2.894) 8.499 2.076 ( 2.957 8.724 -1.138) 9.281 2.649 ( 6.443 5.025 -0.437) 8.182 2.693 ( 3.834 3.935 -4.947) 7.393 3.030 ( 12.682 9.411 5.872) 16.849 3.122 ( 13.405 9.405 -2.768) 16.607 5.437 ( 2.798 9.928 0.172) 10.316 5.441 ( 2.687 10.273 -0.164) 10.620 5.924 ( -2.832 -4.159 2.426) 5.586 5.980 ( -2.228 -4.225 -2.385) 5.339 6.207 ( -1.490 -0.021 0.269) 1.514 6.212 ( -1.550 0.173 -0.218) 1.575 6.681 ( 0.025 1.648 5.085) 5.346 6.810 ( 0.953 0.138 -6.250) 6.324 7.309 ( -11.964 -9.121 2.933) 15.327 7.360 ( -10.568 -9.648 -1.696) 14.410 8.122 ( 1.924 -2.609 3.310) 4.633 8.202 ( 1.278 -4.328 -3.605) 5.776 ======================= Grid point 230 (36/48) ======================= q-point: ( 0.29 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.471 ( 3.534 5.383 2.019) 6.748 1.517 ( 2.900 3.619 -1.777) 4.966 1.688 ( 6.308 7.990 0.239) 10.183 1.710 ( 4.988 8.014 -1.051) 9.498 1.974 ( 1.591 -0.552 3.079) 3.510 2.044 ( 3.370 0.146 -3.160) 4.622 2.182 ( 2.553 10.122 1.916) 10.613 2.217 ( 2.137 9.587 -1.306) 9.909 2.669 ( -3.882 -1.631 2.071) 4.693 2.719 ( -0.171 -0.364 -2.195) 2.232 3.434 ( 16.853 10.773 2.255) 20.129 3.483 ( 14.786 9.140 -2.022) 17.500 5.561 ( -0.311 10.214 0.159) 10.220 5.565 ( -0.485 10.369 -0.162) 10.382 5.881 ( 1.374 -1.307 2.617) 3.232 5.942 ( 1.959 -2.233 -2.621) 3.961 6.187 ( -1.278 0.591 0.369) 1.456 6.196 ( -1.505 1.584 -0.405) 2.223 6.645 ( -3.543 -1.221 5.021) 6.265 6.780 ( -3.419 -2.338 -7.152) 8.265 7.044 ( -8.290 -5.950 3.874) 10.915 7.103 ( -8.644 -6.553 -1.743) 10.986 8.114 ( 1.039 -4.534 2.316) 5.196 8.167 ( 0.488 -5.610 -2.266) 6.070 ======================= Grid point 231 (37/48) ======================= q-point: ( 0.36 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.564 ( 0.778 5.306 0.910) 5.440 1.584 ( 0.880 3.351 -0.800) 3.556 1.850 ( 4.002 4.661 0.218) 6.147 1.862 ( 4.183 5.071 -0.641) 6.605 1.996 ( 1.707 -0.089 3.561) 3.950 2.089 ( 2.508 -0.442 -4.711) 5.355 2.316 ( 0.938 9.292 1.177) 9.413 2.340 ( 0.405 9.503 -1.003) 9.565 2.568 ( -5.224 -1.228 4.746) 7.164 2.660 ( -4.242 -2.312 -3.381) 5.897 3.790 ( 12.507 6.622 0.202) 14.153 3.795 ( 11.068 5.702 -0.212) 12.452 5.642 ( -2.211 10.276 0.081) 10.512 5.644 ( -2.219 10.146 -0.093) 10.386 5.943 ( 5.366 0.890 2.907) 6.168 6.003 ( 6.224 -0.818 -2.241) 6.665 6.165 ( -1.737 0.262 0.627) 1.865 6.185 ( -1.903 2.498 -1.144) 3.342 6.545 ( -5.521 -2.133 3.437) 6.844 6.634 ( -7.858 -4.912 -4.454) 10.282 6.913 ( -2.552 -2.167 1.304) 3.593 6.934 ( -4.778 -3.136 -0.623) 5.749 8.068 ( -0.395 -4.850 1.879) 5.217 8.109 ( -0.457 -5.189 -1.673) 5.471 ======================= Grid point 232 (38/48) ======================= q-point: ( 0.43 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.608 ( -1.530 4.624 0.476) 4.894 1.619 ( -0.613 2.923 -0.510) 3.029 1.916 ( 0.528 1.497 1.107) 1.935 1.945 ( 1.139 1.518 -1.463) 2.396 2.037 ( 1.730 0.613 3.763) 4.187 2.132 ( 2.377 -0.540 -4.854) 5.432 2.397 ( -0.858 6.294 1.350) 6.494 2.417 ( -2.697 8.674 -0.696) 9.110 2.486 ( -3.791 2.749 3.417) 5.797 2.554 ( -4.143 -1.784 -2.494) 5.154 3.972 ( 4.618 1.095 0.815) 4.815 3.991 ( 5.152 1.341 -0.807) 5.385 5.691 ( -3.986 9.639 0.022) 10.431 5.692 ( -3.801 9.277 -0.068) 10.026 6.053 ( 4.880 -0.088 3.557) 6.040 6.091 ( 3.006 -1.665 -0.273) 3.447 6.140 ( -1.012 0.316 1.416) 1.769 6.217 ( 3.286 3.065 -4.524) 6.376 6.428 ( -5.029 -0.849 0.571) 5.132 6.449 ( -6.797 -2.912 -1.079) 7.473 6.861 ( -1.358 -0.994 0.675) 1.814 6.875 ( -0.444 -0.829 -0.581) 1.106 8.020 ( 0.254 -3.362 1.876) 3.858 8.060 ( 0.332 -3.444 -1.591) 3.808 ======================= Grid point 242 (39/48) ======================= q-point: ( 0.14 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.307 ( 3.297 5.710 0.521) 6.614 1.346 ( 5.872 10.171 -3.295) 12.198 1.503 ( 6.212 10.760 2.461) 12.665 1.559 ( 4.579 7.931 -1.398) 9.264 1.900 ( 3.967 6.872 0.484) 7.950 1.915 ( 1.187 2.055 0.671) 2.466 2.056 ( 4.965 8.599 3.958) 10.689 2.103 ( 4.819 8.347 -1.034) 9.693 2.611 ( 5.112 8.855 -0.228) 10.227 2.668 ( 5.823 10.086 -5.887) 13.050 2.970 ( 4.651 8.055 5.246) 10.679 3.042 ( 6.944 12.028 -1.779) 14.002 5.473 ( 6.051 10.481 0.339) 12.107 5.481 ( 6.203 10.744 -0.321) 12.410 5.928 ( -3.061 -5.301 2.173) 6.496 5.977 ( -2.963 -5.132 -2.101) 6.288 6.223 ( -0.401 -0.694 0.388) 0.890 6.229 ( -0.281 -0.488 -0.175) 0.589 6.694 ( 1.677 2.904 4.066) 5.271 6.795 ( 0.712 1.233 -4.747) 4.956 7.378 ( -7.290 -12.627 1.586) 14.666 7.407 ( -7.193 -12.459 -0.953) 14.418 8.083 ( -1.133 -1.962 2.920) 3.696 8.155 ( -1.981 -3.431 -3.346) 5.186 ======================= Grid point 243 (40/48) ======================= q-point: ( 0.21 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.437 ( 4.780 5.306 3.170) 7.813 1.527 ( 3.849 7.482 -4.939) 9.756 1.692 ( 3.697 9.768 2.076) 10.649 1.720 ( 4.435 7.653 -0.035) 8.846 1.947 ( 1.179 0.386 2.009) 2.361 2.001 ( 2.485 1.985 -3.065) 4.417 2.239 ( 3.994 10.584 2.565) 11.600 2.279 ( 3.626 10.529 -1.181) 11.199 2.716 ( 0.028 1.264 -0.545) 1.377 2.717 ( -0.739 0.687 -0.654) 1.203 3.282 ( 13.239 13.767 3.041) 19.341 3.345 ( 11.732 11.758 -2.441) 16.788 5.667 ( 2.641 11.421 0.402) 11.729 5.676 ( 2.453 11.532 -0.417) 11.797 5.852 ( -0.475 -2.434 1.966) 3.164 5.898 ( -0.549 -3.330 -2.017) 3.932 6.214 ( -0.716 0.575 0.957) 1.327 6.232 ( -0.260 1.756 -0.620) 1.880 6.695 ( -2.346 -0.103 3.640) 4.331 6.796 ( -0.805 -1.554 -5.311) 5.592 7.121 ( -8.169 -9.196 3.175) 12.703 7.174 ( -7.234 -7.910 -1.670) 10.848 8.040 ( 0.360 -4.952 1.935) 5.329 8.086 ( 0.049 -6.388 -2.026) 6.702 ======================= Grid point 244 (41/48) ======================= q-point: ( 0.29 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 296 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.570 ( 3.714 5.187 2.866) 6.994 1.637 ( 0.662 6.723 -2.861) 7.336 1.840 ( 3.342 6.197 0.837) 7.090 1.860 ( 3.038 6.205 -0.734) 6.947 1.971 ( 1.361 0.492 2.783) 3.137 2.047 ( 1.946 0.391 -4.081) 4.538 2.387 ( 1.523 7.311 2.431) 7.854 2.428 ( 1.375 8.943 -1.350) 9.148 2.674 ( -4.230 4.321 2.127) 6.410 2.721 ( -2.581 2.476 -1.955) 4.076 3.635 ( 12.790 8.431 0.813) 15.340 3.654 ( 11.413 7.089 -0.813) 13.460 5.794 ( -0.853 10.594 0.120) 10.629 5.796 ( -0.912 9.618 -0.002) 9.661 5.873 ( 3.067 1.865 1.292) 3.815 5.904 ( 3.570 -0.665 -1.296) 3.856 6.199 ( -1.620 0.179 1.923) 2.521 6.244 ( -1.394 2.653 -1.905) 3.550 6.629 ( -4.853 0.312 2.944) 5.685 6.714 ( -5.612 -3.616 -4.814) 8.230 6.939 ( -3.953 -4.057 3.443) 6.629 6.995 ( -6.011 -3.825 -1.765) 7.340 7.996 ( 1.265 -6.333 1.513) 6.633 8.030 ( 1.466 -7.094 -1.440) 7.385 ======================= Grid point 245 (42/48) ======================= q-point: ( 0.36 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.663 ( 0.615 5.109 1.412) 5.337 1.696 ( -1.392 6.234 -1.439) 6.548 1.923 ( 1.036 2.438 0.605) 2.717 1.939 ( 1.283 1.817 -0.853) 2.382 2.019 ( 2.048 2.558 2.774) 4.293 2.090 ( 2.256 0.539 -3.596) 4.279 2.427 ( -0.442 -0.211 3.253) 3.290 2.488 ( -0.233 0.563 -2.223) 2.305 2.677 ( -6.005 13.302 1.364) 14.658 2.708 ( -6.145 11.066 -1.258) 12.720 3.872 ( 7.823 1.450 0.471) 7.970 3.883 ( 8.625 2.107 -0.481) 8.892 5.862 ( -2.099 8.748 0.157) 8.997 5.865 ( -2.309 9.177 -0.173) 9.465 5.968 ( 4.325 2.113 1.491) 5.039 5.991 ( 5.365 0.279 -0.594) 5.405 6.163 ( -1.738 -0.780 2.386) 3.053 6.233 ( -1.624 1.043 -3.580) 4.067 6.540 ( -6.504 2.789 1.028) 7.151 6.564 ( -7.779 -0.064 -0.993) 7.843 6.875 ( -0.968 -1.541 0.511) 1.890 6.885 ( -3.168 -1.639 -0.374) 3.586 7.959 ( 1.457 -5.311 1.549) 5.721 7.993 ( 1.817 -5.721 -1.401) 6.164 ======================= Grid point 246 (43/48) ======================= q-point: ( 0.43 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.694 ( -2.420 4.191 0.851) 4.914 1.715 ( -3.139 5.438 -0.995) 6.357 1.940 ( -0.504 0.872 0.861) 1.325 1.959 ( -0.044 0.076 -0.762) 0.767 2.057 ( -0.610 1.057 2.378) 2.673 2.119 ( 0.421 -0.729 -3.189) 3.298 2.413 ( 0.521 -0.903 3.059) 3.232 2.475 ( 1.190 -2.062 -2.412) 3.389 2.698 ( -8.499 14.720 0.909) 17.021 2.716 ( -8.513 14.746 -0.675) 17.040 3.950 ( 1.956 -3.387 0.939) 4.022 3.972 ( 1.911 -3.310 -0.931) 3.934 5.891 ( -4.813 8.336 0.062) 9.626 5.894 ( -4.697 8.136 -0.281) 9.399 6.028 ( 1.232 -2.134 2.222) 3.319 6.050 ( 1.565 -2.710 -0.126) 3.132 6.140 ( 0.090 -0.157 2.230) 2.238 6.224 ( 0.907 -1.571 -4.706) 5.044 6.481 ( -4.318 7.478 0.481) 8.649 6.489 ( -3.702 6.412 -0.165) 7.406 6.849 ( 0.041 -0.070 0.519) 0.525 6.860 ( 0.372 -0.644 -0.401) 0.845 7.944 ( 2.274 -3.939 1.625) 4.830 7.980 ( 2.515 -4.357 -1.455) 5.237 ======================= Grid point 258 (44/48) ======================= q-point: ( 0.21 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.561 ( 3.801 6.584 4.321) 8.745 1.667 ( 3.432 5.945 -5.041) 8.516 1.849 ( 3.080 5.335 1.325) 6.301 1.863 ( 3.606 6.245 -0.027) 7.211 1.960 ( 0.671 1.163 2.361) 2.716 2.032 ( 0.815 1.412 -4.017) 4.336 2.409 ( 2.614 4.527 2.912) 5.984 2.460 ( 3.215 5.569 -1.700) 6.652 2.748 ( 2.572 4.454 0.877) 5.218 2.769 ( 2.559 4.432 -0.934) 5.202 3.554 ( 7.033 12.182 1.064) 14.107 3.578 ( 6.055 10.487 -1.037) 12.154 5.822 ( 0.779 1.349 0.899) 1.799 5.841 ( 1.162 2.013 -0.739) 2.439 5.884 ( 3.098 5.366 0.070) 6.196 5.890 ( 4.105 7.110 -0.341) 8.217 6.217 ( -0.420 -0.728 2.594) 2.727 6.273 ( 0.912 1.579 -2.248) 2.894 6.687 ( -0.235 -0.407 1.527) 1.598 6.743 ( -2.007 -3.476 -3.868) 5.574 6.958 ( -3.448 -5.972 4.868) 8.441 7.040 ( -3.141 -5.440 -2.766) 6.864 7.941 ( -2.413 -4.180 1.111) 4.953 7.967 ( -2.737 -4.741 -1.094) 5.583 ======================= Grid point 259 (45/48) ======================= q-point: ( 0.29 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.691 ( 2.903 6.446 3.107) 7.722 1.765 ( 0.915 5.550 -3.231) 6.487 1.926 ( 1.343 2.362 0.587) 2.780 1.939 ( 0.770 1.507 -0.581) 1.789 2.013 ( 2.214 3.360 2.067) 4.524 2.066 ( 1.392 1.382 -2.768) 3.393 2.428 ( -0.723 -1.273 2.965) 3.307 2.486 ( -0.953 -2.120 -2.125) 3.149 2.887 ( 0.490 14.231 0.839) 14.264 2.906 ( 0.637 14.040 -0.796) 14.077 3.755 ( 6.102 2.344 0.231) 6.541 3.760 ( 6.883 3.378 -0.235) 7.671 5.872 ( 2.942 2.250 0.350) 3.720 5.879 ( 2.909 1.227 -0.278) 3.169 5.989 ( 1.037 5.170 0.375) 5.286 5.999 ( 0.706 5.478 -0.427) 5.539 6.180 ( -1.718 -2.051 3.396) 4.323 6.266 ( -1.974 -1.538 -4.088) 4.793 6.678 ( -2.952 3.768 0.257) 4.794 6.683 ( -3.065 2.093 -0.147) 3.714 6.879 ( -1.424 -2.000 2.937) 3.828 6.940 ( -4.101 -1.726 -2.424) 5.067 7.882 ( 0.624 -4.353 0.917) 4.492 7.902 ( 0.962 -4.963 -0.864) 5.129 ======================= Grid point 260 (46/48) ======================= q-point: ( 0.36 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.777 ( -0.725 5.803 1.507) 6.039 1.814 ( -1.486 5.079 -1.670) 5.549 1.952 ( 0.031 0.683 0.246) 0.726 1.956 ( 0.440 0.350 -0.099) 0.571 2.069 ( 0.833 1.553 1.242) 2.156 2.101 ( 1.211 0.195 -1.524) 1.957 2.400 ( -0.323 -1.331 1.723) 2.201 2.436 ( -0.262 -3.327 -1.455) 3.641 3.006 ( -6.160 17.286 0.887) 18.372 3.026 ( -6.236 17.513 -0.864) 18.611 3.839 ( 5.048 -4.839 0.996) 7.064 3.862 ( 5.300 -4.362 -0.968) 6.932 5.954 ( 3.070 2.090 0.278) 3.724 5.961 ( 3.759 1.759 -0.204) 4.155 6.020 ( -1.353 1.261 0.258) 1.867 6.026 ( -1.804 1.623 -0.072) 2.428 6.137 ( -0.554 -1.409 2.488) 2.912 6.201 ( -1.142 -3.229 -2.983) 4.542 6.662 ( -5.226 7.655 0.574) 9.286 6.676 ( -4.763 9.060 -0.623) 10.255 6.852 ( -0.336 -0.514 0.737) 0.959 6.866 ( -2.627 -0.123 -0.490) 2.675 7.866 ( 1.851 -3.626 0.956) 4.182 7.887 ( 2.393 -4.453 -0.915) 5.138 ======================= Grid point 273 (47/48) ======================= q-point: ( 0.29 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.818 ( 3.121 5.405 1.625) 6.449 1.856 ( 1.879 3.255 -1.669) 4.113 1.953 ( 0.390 0.675 0.079) 0.783 1.955 ( 0.250 0.434 -0.087) 0.508 2.066 ( 0.802 1.390 0.983) 1.882 2.090 ( 0.358 0.620 -1.062) 1.281 2.396 ( -0.735 -1.273 1.142) 1.862 2.421 ( -1.808 -3.131 -1.006) 3.753 3.179 ( 7.552 13.080 0.977) 15.135 3.202 ( 7.734 13.396 -1.010) 15.501 3.742 ( -1.893 -3.278 1.001) 3.915 3.765 ( -1.543 -2.672 -0.944) 3.227 5.926 ( 1.850 3.204 0.283) 3.710 5.933 ( 2.044 3.540 -0.355) 4.102 6.041 ( 0.131 0.228 0.251) 0.363 6.054 ( 0.165 0.285 -0.815) 0.879 6.133 ( -1.140 -1.974 2.948) 3.726 6.194 ( -2.424 -4.199 -2.444) 5.429 6.751 ( 2.003 3.470 -0.090) 4.008 6.761 ( 1.605 2.779 -0.713) 3.287 6.859 ( 0.172 0.298 2.917) 2.937 6.914 ( -0.662 -1.147 -2.004) 2.402 7.826 ( -0.743 -1.287 0.413) 1.542 7.835 ( -1.001 -1.733 -0.393) 2.040 ======================= Grid point 274 (48/48) ======================= q-point: ( 0.36 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 5.20e-05 5.20e-05 1.65e-04 1.65e-04 1.65e-04 1.65e-04 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.880 ( -2.027 3.511 0.465) 4.080 1.892 ( -1.271 2.202 -0.569) 2.605 1.960 ( -0.084 0.145 0.074) 0.184 1.960 ( -0.044 0.075 0.042) 0.097 2.074 ( 0.377 -0.653 0.858) 1.142 2.094 ( 0.293 -0.508 -0.870) 1.049 2.383 ( 0.190 -0.330 0.088) 0.391 2.385 ( 0.773 -1.338 -0.092) 1.548 3.345 ( -8.377 14.510 1.089) 16.790 3.371 ( -8.520 14.757 -1.210) 17.083 3.679 ( 5.744 -9.948 1.322) 11.563 3.708 ( 5.700 -9.872 -1.185) 11.461 5.980 ( -0.104 0.180 0.457) 0.502 5.989 ( -0.640 1.108 -0.459) 1.359 6.031 ( -0.092 0.158 -0.703) 0.726 6.031 ( -0.058 0.102 -0.694) 0.704 6.115 ( 0.315 -0.546 2.048) 2.143 6.141 ( 1.273 -2.205 -0.274) 2.561 6.779 ( -1.594 2.762 -0.131) 3.192 6.795 ( -0.629 1.089 -1.623) 2.054 6.884 ( -1.128 1.955 0.729) 2.372 6.886 ( -1.672 2.896 0.653) 3.407 7.815 ( 0.723 -1.253 0.193) 1.459 7.820 ( 1.084 -1.878 -0.194) 2.178 =================== 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/14112 10.0 557.660 557.660 364.344 -0.000 -0.000 0.000 3/14112 20.0 117.575 117.575 74.231 -0.000 -0.000 0.000 3/14112 30.0 59.928 59.928 34.725 -0.000 -0.000 0.000 3/14112 40.0 40.064 40.064 21.906 -0.000 -0.000 0.000 3/14112 50.0 29.432 29.432 15.675 -0.000 0.000 0.000 3/14112 60.0 22.791 22.791 12.036 -0.000 0.000 0.000 3/14112 70.0 18.369 18.369 9.693 -0.000 0.000 0.000 3/14112 80.0 15.289 15.289 8.082 -0.000 0.000 0.000 3/14112 90.0 13.059 13.059 6.920 -0.000 0.000 0.000 3/14112 100.0 11.386 11.386 6.047 -0.000 0.000 0.000 3/14112 110.0 10.092 10.092 5.371 -0.000 0.000 0.000 3/14112 120.0 9.065 9.065 4.832 -0.000 0.000 0.000 3/14112 130.0 8.231 8.231 4.394 -0.000 0.000 0.000 3/14112 140.0 7.541 7.541 4.030 -0.000 0.000 0.000 3/14112 150.0 6.961 6.961 3.723 -0.000 0.000 0.000 3/14112 160.0 6.466 6.466 3.461 -0.000 0.000 0.000 3/14112 170.0 6.039 6.039 3.234 -0.000 0.000 0.000 3/14112 180.0 5.666 5.666 3.036 -0.000 0.000 0.000 3/14112 190.0 5.338 5.338 2.861 -0.000 0.000 0.000 3/14112 200.0 5.047 5.047 2.706 -0.000 0.000 0.000 3/14112 210.0 4.787 4.787 2.568 -0.000 0.000 0.000 3/14112 220.0 4.553 4.553 2.443 -0.000 0.000 0.000 3/14112 230.0 4.342 4.342 2.330 -0.000 0.000 0.000 3/14112 240.0 4.149 4.149 2.227 -0.000 0.000 0.000 3/14112 250.0 3.974 3.974 2.133 -0.000 0.000 0.000 3/14112 260.0 3.813 3.813 2.047 -0.000 0.000 0.000 3/14112 270.0 3.665 3.665 1.967 -0.000 0.000 0.000 3/14112 280.0 3.528 3.528 1.894 -0.000 0.000 0.000 3/14112 290.0 3.401 3.401 1.826 -0.000 0.000 0.000 3/14112 300.0 3.283 3.283 1.763 -0.000 0.000 0.000 3/14112 310.0 3.173 3.173 1.704 -0.000 0.000 0.000 3/14112 320.0 3.071 3.071 1.649 -0.000 0.000 0.000 3/14112 330.0 2.975 2.975 1.598 -0.000 0.000 0.000 3/14112 340.0 2.884 2.884 1.549 -0.000 0.000 0.000 3/14112 350.0 2.800 2.800 1.504 -0.000 0.000 0.000 3/14112 360.0 2.720 2.720 1.461 -0.000 0.000 0.000 3/14112 370.0 2.645 2.645 1.420 -0.000 0.000 0.000 3/14112 380.0 2.573 2.573 1.382 -0.000 0.000 0.000 3/14112 390.0 2.506 2.506 1.346 -0.000 0.000 0.000 3/14112 400.0 2.442 2.442 1.312 -0.000 0.000 0.000 3/14112 410.0 2.381 2.381 1.279 -0.000 0.000 0.000 3/14112 420.0 2.323 2.323 1.248 -0.000 0.000 0.000 3/14112 430.0 2.268 2.268 1.218 -0.000 0.000 0.000 3/14112 440.0 2.216 2.216 1.190 -0.000 0.000 0.000 3/14112 450.0 2.166 2.166 1.163 -0.000 0.000 0.000 3/14112 460.0 2.118 2.118 1.138 -0.000 0.000 0.000 3/14112 470.0 2.072 2.072 1.113 -0.000 0.000 0.000 3/14112 480.0 2.029 2.029 1.090 -0.000 0.000 0.000 3/14112 490.0 1.987 1.987 1.067 -0.000 0.000 0.000 3/14112 500.0 1.946 1.946 1.046 -0.000 0.000 0.000 3/14112 510.0 1.908 1.908 1.025 -0.000 0.000 0.000 3/14112 520.0 1.871 1.871 1.005 -0.000 0.000 0.000 3/14112 530.0 1.835 1.835 0.986 -0.000 0.000 0.000 3/14112 540.0 1.801 1.801 0.967 -0.000 0.000 0.000 3/14112 550.0 1.767 1.767 0.949 -0.000 0.000 0.000 3/14112 560.0 1.736 1.736 0.932 -0.000 0.000 0.000 3/14112 570.0 1.705 1.705 0.916 -0.000 0.000 0.000 3/14112 580.0 1.675 1.675 0.900 -0.000 0.000 0.000 3/14112 590.0 1.647 1.647 0.885 -0.000 0.000 0.000 3/14112 600.0 1.619 1.619 0.870 -0.000 0.000 0.000 3/14112 610.0 1.592 1.592 0.855 -0.000 0.000 0.000 3/14112 620.0 1.566 1.566 0.841 -0.000 0.000 0.000 3/14112 630.0 1.541 1.541 0.828 -0.000 0.000 0.000 3/14112 640.0 1.517 1.517 0.815 -0.000 0.000 0.000 3/14112 650.0 1.493 1.493 0.802 -0.000 0.000 0.000 3/14112 660.0 1.471 1.471 0.790 -0.000 0.000 0.000 3/14112 670.0 1.449 1.449 0.778 -0.000 0.000 0.000 3/14112 680.0 1.427 1.427 0.767 -0.000 0.000 0.000 3/14112 690.0 1.406 1.406 0.755 -0.000 0.000 0.000 3/14112 700.0 1.386 1.386 0.745 -0.000 0.000 0.000 3/14112 710.0 1.367 1.367 0.734 -0.000 0.000 0.000 3/14112 720.0 1.347 1.347 0.724 -0.000 0.000 0.000 3/14112 730.0 1.329 1.329 0.714 -0.000 0.000 0.000 3/14112 740.0 1.311 1.311 0.704 -0.000 0.000 0.000 3/14112 750.0 1.293 1.293 0.695 -0.000 0.000 0.000 3/14112 760.0 1.276 1.276 0.686 -0.000 0.000 0.000 3/14112 770.0 1.260 1.260 0.677 -0.000 0.000 0.000 3/14112 780.0 1.243 1.243 0.668 -0.000 0.000 0.000 3/14112 790.0 1.228 1.228 0.659 -0.000 0.000 0.000 3/14112 800.0 1.212 1.212 0.651 -0.000 0.000 0.000 3/14112 810.0 1.197 1.197 0.643 -0.000 0.000 0.000 3/14112 820.0 1.182 1.182 0.635 -0.000 0.000 0.000 3/14112 830.0 1.168 1.168 0.627 -0.000 0.000 0.000 3/14112 840.0 1.154 1.154 0.620 -0.000 0.000 0.000 3/14112 850.0 1.141 1.141 0.613 -0.000 0.000 0.000 3/14112 860.0 1.127 1.127 0.605 -0.000 0.000 0.000 3/14112 870.0 1.114 1.114 0.598 -0.000 0.000 0.000 3/14112 880.0 1.102 1.102 0.592 -0.000 0.000 0.000 3/14112 890.0 1.089 1.089 0.585 -0.000 0.000 0.000 3/14112 900.0 1.077 1.077 0.578 -0.000 0.000 0.000 3/14112 910.0 1.065 1.065 0.572 -0.000 0.000 0.000 3/14112 920.0 1.054 1.054 0.566 -0.000 0.000 0.000 3/14112 930.0 1.042 1.042 0.560 -0.000 0.000 0.000 3/14112 940.0 1.031 1.031 0.554 -0.000 0.000 0.000 3/14112 950.0 1.020 1.020 0.548 -0.000 0.000 0.000 3/14112 960.0 1.010 1.010 0.542 -0.000 0.000 0.000 3/14112 970.0 0.999 0.999 0.537 -0.000 0.000 0.000 3/14112 980.0 0.989 0.989 0.531 -0.000 0.000 0.000 3/14112 990.0 0.979 0.979 0.526 -0.000 0.000 0.000 3/14112 1000.0 0.969 0.969 0.520 -0.000 0.000 0.000 3/14112 Thermal conductivity related properties were written into "kappa-m14143.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-09 04:49:44]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|