----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-07 23:10:49]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [4 4 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.214008208050170 0.000000000000000 0.000000000000000 b -2.107004104025085 3.649438159927587 0.000000000000000 c 0.000000000000000 0.000000000000000 15.904157939999999 Atomic positions (fractional): *1 S 0.33333333333333 0.66666666666667 0.90485241533818 32.065 2 S 0.66666666666667 0.33333333333333 0.40485241533818 32.065 3 S 0.33333333333333 0.66666666666667 0.59514758466182 32.065 4 S 0.66666666666667 0.33333333333333 0.09514758466182 32.065 *5 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 6 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 *7 Pr 0.00000000000000 0.00000000000000 0.00000000000000 140.908 8 Pr 0.00000000000000 0.00000000000000 0.50000000000000 140.908 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.214008208050170 0.000000000000000 0.000000000000000 b -2.107004104025085 3.649438159927587 0.000000000000000 c 0.000000000000000 0.000000000000000 15.904157939999999 Atomic positions (fractional): *1 S 0.33333333333333 0.66666666666667 0.90485241533818 32.065 > 1 2 S 0.66666666666667 0.33333333333333 0.40485241533818 32.065 > 2 3 S 0.33333333333333 0.66666666666667 0.59514758466182 32.065 > 3 4 S 0.66666666666667 0.33333333333333 0.09514758466182 32.065 > 4 *5 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 > 5 6 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 > 6 *7 Pr 0.00000000000000 0.00000000000000 0.00000000000000 140.908 > 7 8 Pr 0.00000000000000 0.00000000000000 0.50000000000000 140.908 > 8 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.856032832200679 0.000000000000000 0.000000000000000 b -8.428016416100339 14.597752639710349 0.000000000000000 c 0.000000000000000 0.000000000000000 15.904157939999999 Atomic positions (fractional): *1 S 0.08333333333333 0.16666666666667 0.90485241533818 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.90485241533818 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.90485241533818 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.90485241533818 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.90485241533818 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.90485241533818 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.90485241533818 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.90485241533818 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.90485241533818 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.90485241533818 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.90485241533818 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.90485241533818 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.90485241533818 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.90485241533818 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.90485241533818 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.90485241533818 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.40485241533818 32.065 > 2 18 S 0.41666666666667 0.08333333333333 0.40485241533818 32.065 > 2 19 S 0.66666666666667 0.08333333333333 0.40485241533818 32.065 > 2 20 S 0.91666666666667 0.08333333333333 0.40485241533818 32.065 > 2 21 S 0.16666666666667 0.33333333333333 0.40485241533818 32.065 > 2 22 S 0.41666666666667 0.33333333333333 0.40485241533818 32.065 > 2 23 S 0.66666666666667 0.33333333333333 0.40485241533818 32.065 > 2 24 S 0.91666666666667 0.33333333333333 0.40485241533818 32.065 > 2 25 S 0.16666666666667 0.58333333333333 0.40485241533818 32.065 > 2 26 S 0.41666666666667 0.58333333333333 0.40485241533818 32.065 > 2 27 S 0.66666666666667 0.58333333333333 0.40485241533818 32.065 > 2 28 S 0.91666666666667 0.58333333333333 0.40485241533818 32.065 > 2 29 S 0.16666666666667 0.83333333333333 0.40485241533818 32.065 > 2 30 S 0.41666666666667 0.83333333333333 0.40485241533818 32.065 > 2 31 S 0.66666666666667 0.83333333333333 0.40485241533818 32.065 > 2 32 S 0.91666666666667 0.83333333333333 0.40485241533818 32.065 > 2 33 S 0.08333333333333 0.16666666666667 0.59514758466182 32.065 > 3 34 S 0.33333333333333 0.16666666666667 0.59514758466182 32.065 > 3 35 S 0.58333333333333 0.16666666666667 0.59514758466182 32.065 > 3 36 S 0.83333333333333 0.16666666666667 0.59514758466182 32.065 > 3 37 S 0.08333333333333 0.41666666666667 0.59514758466182 32.065 > 3 38 S 0.33333333333333 0.41666666666667 0.59514758466182 32.065 > 3 39 S 0.58333333333333 0.41666666666667 0.59514758466182 32.065 > 3 40 S 0.83333333333333 0.41666666666667 0.59514758466182 32.065 > 3 41 S 0.08333333333333 0.66666666666667 0.59514758466182 32.065 > 3 42 S 0.33333333333333 0.66666666666667 0.59514758466182 32.065 > 3 43 S 0.58333333333333 0.66666666666667 0.59514758466182 32.065 > 3 44 S 0.83333333333333 0.66666666666667 0.59514758466182 32.065 > 3 45 S 0.08333333333333 0.91666666666667 0.59514758466182 32.065 > 3 46 S 0.33333333333333 0.91666666666667 0.59514758466182 32.065 > 3 47 S 0.58333333333333 0.91666666666667 0.59514758466182 32.065 > 3 48 S 0.83333333333333 0.91666666666667 0.59514758466182 32.065 > 3 49 S 0.16666666666667 0.08333333333333 0.09514758466182 32.065 > 4 50 S 0.41666666666667 0.08333333333333 0.09514758466182 32.065 > 4 51 S 0.66666666666667 0.08333333333333 0.09514758466182 32.065 > 4 52 S 0.91666666666667 0.08333333333333 0.09514758466182 32.065 > 4 53 S 0.16666666666667 0.33333333333333 0.09514758466182 32.065 > 4 54 S 0.41666666666667 0.33333333333333 0.09514758466182 32.065 > 4 55 S 0.66666666666667 0.33333333333333 0.09514758466182 32.065 > 4 56 S 0.91666666666667 0.33333333333333 0.09514758466182 32.065 > 4 57 S 0.16666666666667 0.58333333333333 0.09514758466182 32.065 > 4 58 S 0.41666666666667 0.58333333333333 0.09514758466182 32.065 > 4 59 S 0.66666666666667 0.58333333333333 0.09514758466182 32.065 > 4 60 S 0.91666666666667 0.58333333333333 0.09514758466182 32.065 > 4 61 S 0.16666666666667 0.83333333333333 0.09514758466182 32.065 > 4 62 S 0.41666666666667 0.83333333333333 0.09514758466182 32.065 > 4 63 S 0.66666666666667 0.83333333333333 0.09514758466182 32.065 > 4 64 S 0.91666666666667 0.83333333333333 0.09514758466182 32.065 > 4 *65 Cs 0.16666666666667 0.08333333333333 0.75000000000000 132.905 > 5 66 Cs 0.41666666666667 0.08333333333333 0.75000000000000 132.905 > 5 67 Cs 0.66666666666667 0.08333333333333 0.75000000000000 132.905 > 5 68 Cs 0.91666666666667 0.08333333333333 0.75000000000000 132.905 > 5 69 Cs 0.16666666666667 0.33333333333333 0.75000000000000 132.905 > 5 70 Cs 0.41666666666667 0.33333333333333 0.75000000000000 132.905 > 5 71 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 > 5 72 Cs 0.91666666666667 0.33333333333333 0.75000000000000 132.905 > 5 73 Cs 0.16666666666667 0.58333333333333 0.75000000000000 132.905 > 5 74 Cs 0.41666666666667 0.58333333333333 0.75000000000000 132.905 > 5 75 Cs 0.66666666666667 0.58333333333333 0.75000000000000 132.905 > 5 76 Cs 0.91666666666667 0.58333333333333 0.75000000000000 132.905 > 5 77 Cs 0.16666666666667 0.83333333333333 0.75000000000000 132.905 > 5 78 Cs 0.41666666666667 0.83333333333333 0.75000000000000 132.905 > 5 79 Cs 0.66666666666667 0.83333333333333 0.75000000000000 132.905 > 5 80 Cs 0.91666666666667 0.83333333333333 0.75000000000000 132.905 > 5 81 Cs 0.08333333333333 0.16666666666667 0.25000000000000 132.905 > 6 82 Cs 0.33333333333333 0.16666666666667 0.25000000000000 132.905 > 6 83 Cs 0.58333333333333 0.16666666666667 0.25000000000000 132.905 > 6 84 Cs 0.83333333333333 0.16666666666667 0.25000000000000 132.905 > 6 85 Cs 0.08333333333333 0.41666666666667 0.25000000000000 132.905 > 6 86 Cs 0.33333333333333 0.41666666666667 0.25000000000000 132.905 > 6 87 Cs 0.58333333333333 0.41666666666667 0.25000000000000 132.905 > 6 88 Cs 0.83333333333333 0.41666666666667 0.25000000000000 132.905 > 6 89 Cs 0.08333333333333 0.66666666666667 0.25000000000000 132.905 > 6 90 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 > 6 91 Cs 0.58333333333333 0.66666666666667 0.25000000000000 132.905 > 6 92 Cs 0.83333333333333 0.66666666666667 0.25000000000000 132.905 > 6 93 Cs 0.08333333333333 0.91666666666667 0.25000000000000 132.905 > 6 94 Cs 0.33333333333333 0.91666666666667 0.25000000000000 132.905 > 6 95 Cs 0.58333333333333 0.91666666666667 0.25000000000000 132.905 > 6 96 Cs 0.83333333333333 0.91666666666667 0.25000000000000 132.905 > 6 *97 Pr 0.00000000000000 0.00000000000000 0.00000000000000 140.908 > 7 98 Pr 0.25000000000000 0.00000000000000 0.00000000000000 140.908 > 7 99 Pr 0.50000000000000 0.00000000000000 0.00000000000000 140.908 > 7 100 Pr 0.75000000000000 0.00000000000000 0.00000000000000 140.908 > 7 101 Pr 0.00000000000000 0.25000000000000 0.00000000000000 140.908 > 7 102 Pr 0.25000000000000 0.25000000000000 0.00000000000000 140.908 > 7 103 Pr 0.50000000000000 0.25000000000000 0.00000000000000 140.908 > 7 104 Pr 0.75000000000000 0.25000000000000 0.00000000000000 140.908 > 7 105 Pr 0.00000000000000 0.50000000000000 0.00000000000000 140.908 > 7 106 Pr 0.25000000000000 0.50000000000000 0.00000000000000 140.908 > 7 107 Pr 0.50000000000000 0.50000000000000 0.00000000000000 140.908 > 7 108 Pr 0.75000000000000 0.50000000000000 0.00000000000000 140.908 > 7 109 Pr 0.00000000000000 0.75000000000000 0.00000000000000 140.908 > 7 110 Pr 0.25000000000000 0.75000000000000 0.00000000000000 140.908 > 7 111 Pr 0.50000000000000 0.75000000000000 0.00000000000000 140.908 > 7 112 Pr 0.75000000000000 0.75000000000000 0.00000000000000 140.908 > 7 113 Pr 0.00000000000000 0.00000000000000 0.50000000000000 140.908 > 8 114 Pr 0.25000000000000 0.00000000000000 0.50000000000000 140.908 > 8 115 Pr 0.50000000000000 0.00000000000000 0.50000000000000 140.908 > 8 116 Pr 0.75000000000000 0.00000000000000 0.50000000000000 140.908 > 8 117 Pr 0.00000000000000 0.25000000000000 0.50000000000000 140.908 > 8 118 Pr 0.25000000000000 0.25000000000000 0.50000000000000 140.908 > 8 119 Pr 0.50000000000000 0.25000000000000 0.50000000000000 140.908 > 8 120 Pr 0.75000000000000 0.25000000000000 0.50000000000000 140.908 > 8 121 Pr 0.00000000000000 0.50000000000000 0.50000000000000 140.908 > 8 122 Pr 0.25000000000000 0.50000000000000 0.50000000000000 140.908 > 8 123 Pr 0.50000000000000 0.50000000000000 0.50000000000000 140.908 > 8 124 Pr 0.75000000000000 0.50000000000000 0.50000000000000 140.908 > 8 125 Pr 0.00000000000000 0.75000000000000 0.50000000000000 140.908 > 8 126 Pr 0.25000000000000 0.75000000000000 0.50000000000000 140.908 > 8 127 Pr 0.50000000000000 0.75000000000000 0.50000000000000 140.908 > 8 128 Pr 0.75000000000000 0.75000000000000 0.50000000000000 140.908 > 8 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.7978134 -0.0000000 0.0000000 0.0000000 5.7978134 0.0000000 0.0000000 0.0000000 4.8403995 -------------------------- Born effective charges -------------------------- 1 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 2 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 3 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 4 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 5 Cs 1.0916780 -0.0000000 0.0000000 0.0000000 1.0916780 0.0000000 0.0000000 0.0000000 2.1481525 6 Cs 1.0916780 -0.0000000 0.0000000 0.0000000 1.0916780 0.0000000 0.0000000 0.0000000 2.1481525 7 Pr 4.7109525 -0.0000000 0.0000000 -0.0000000 4.7109525 0.0000000 0.0000000 0.0000000 2.2247901 8 Pr 4.7109525 -0.0000000 0.0000000 -0.0000000 4.7109525 0.0000000 0.0000000 0.0000000 2.2247901 ---------------------------------------------------------------------------- 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: 182 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.008 Solver_block: 80 / 80 - Time: 0.192 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.201 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) 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: 182 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 (xx) 0.00000000 (xx) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-07 23:10:53]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 23:10:54]------------------------- 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.214008208050170 0.000000000000000 0.000000000000000 b -2.107004104025085 3.649438159927587 0.000000000000000 c 0.000000000000000 0.000000000000000 15.904157939999999 Atomic positions (fractional): 1 S 0.33333333333333 0.66666666666667 0.90485241533818 32.065 2 S 0.66666666666667 0.33333333333333 0.40485241533818 32.065 3 S 0.33333333333333 0.66666666666667 0.59514758466182 32.065 4 S 0.66666666666667 0.33333333333333 0.09514758466182 32.065 5 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 6 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 7 Pr 0.00000000000000 0.00000000000000 0.00000000000000 140.908 8 Pr 0.00000000000000 0.00000000000000 0.50000000000000 140.908 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.856032832200679 0.000000000000000 0.000000000000000 b -8.428016416100339 14.597752639710349 0.000000000000000 c 0.000000000000000 0.000000000000000 15.904157939999999 Atomic positions (fractional): 1 S 0.08333333333333 0.16666666666667 0.90485241533818 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.90485241533818 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.90485241533818 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.90485241533818 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.90485241533818 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.90485241533818 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.90485241533818 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.90485241533818 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.90485241533818 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.90485241533818 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.90485241533818 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.90485241533818 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.90485241533818 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.90485241533818 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.90485241533818 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.90485241533818 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.40485241533818 32.065 > 17 18 S 0.41666666666667 0.08333333333333 0.40485241533818 32.065 > 17 19 S 0.66666666666667 0.08333333333333 0.40485241533818 32.065 > 17 20 S 0.91666666666667 0.08333333333333 0.40485241533818 32.065 > 17 21 S 0.16666666666667 0.33333333333333 0.40485241533818 32.065 > 17 22 S 0.41666666666667 0.33333333333333 0.40485241533818 32.065 > 17 23 S 0.66666666666667 0.33333333333333 0.40485241533818 32.065 > 17 24 S 0.91666666666667 0.33333333333333 0.40485241533818 32.065 > 17 25 S 0.16666666666667 0.58333333333333 0.40485241533818 32.065 > 17 26 S 0.41666666666667 0.58333333333333 0.40485241533818 32.065 > 17 27 S 0.66666666666667 0.58333333333333 0.40485241533818 32.065 > 17 28 S 0.91666666666667 0.58333333333333 0.40485241533818 32.065 > 17 29 S 0.16666666666667 0.83333333333333 0.40485241533818 32.065 > 17 30 S 0.41666666666667 0.83333333333333 0.40485241533818 32.065 > 17 31 S 0.66666666666667 0.83333333333333 0.40485241533818 32.065 > 17 32 S 0.91666666666667 0.83333333333333 0.40485241533818 32.065 > 17 33 S 0.08333333333333 0.16666666666667 0.59514758466182 32.065 > 33 34 S 0.33333333333333 0.16666666666667 0.59514758466182 32.065 > 33 35 S 0.58333333333333 0.16666666666667 0.59514758466182 32.065 > 33 36 S 0.83333333333333 0.16666666666667 0.59514758466182 32.065 > 33 37 S 0.08333333333333 0.41666666666667 0.59514758466182 32.065 > 33 38 S 0.33333333333333 0.41666666666667 0.59514758466182 32.065 > 33 39 S 0.58333333333333 0.41666666666667 0.59514758466182 32.065 > 33 40 S 0.83333333333333 0.41666666666667 0.59514758466182 32.065 > 33 41 S 0.08333333333333 0.66666666666667 0.59514758466182 32.065 > 33 42 S 0.33333333333333 0.66666666666667 0.59514758466182 32.065 > 33 43 S 0.58333333333333 0.66666666666667 0.59514758466182 32.065 > 33 44 S 0.83333333333333 0.66666666666667 0.59514758466182 32.065 > 33 45 S 0.08333333333333 0.91666666666667 0.59514758466182 32.065 > 33 46 S 0.33333333333333 0.91666666666667 0.59514758466182 32.065 > 33 47 S 0.58333333333333 0.91666666666667 0.59514758466182 32.065 > 33 48 S 0.83333333333333 0.91666666666667 0.59514758466182 32.065 > 33 49 S 0.16666666666667 0.08333333333333 0.09514758466182 32.065 > 49 50 S 0.41666666666667 0.08333333333333 0.09514758466182 32.065 > 49 51 S 0.66666666666667 0.08333333333333 0.09514758466182 32.065 > 49 52 S 0.91666666666667 0.08333333333333 0.09514758466182 32.065 > 49 53 S 0.16666666666667 0.33333333333333 0.09514758466182 32.065 > 49 54 S 0.41666666666667 0.33333333333333 0.09514758466182 32.065 > 49 55 S 0.66666666666667 0.33333333333333 0.09514758466182 32.065 > 49 56 S 0.91666666666667 0.33333333333333 0.09514758466182 32.065 > 49 57 S 0.16666666666667 0.58333333333333 0.09514758466182 32.065 > 49 58 S 0.41666666666667 0.58333333333333 0.09514758466182 32.065 > 49 59 S 0.66666666666667 0.58333333333333 0.09514758466182 32.065 > 49 60 S 0.91666666666667 0.58333333333333 0.09514758466182 32.065 > 49 61 S 0.16666666666667 0.83333333333333 0.09514758466182 32.065 > 49 62 S 0.41666666666667 0.83333333333333 0.09514758466182 32.065 > 49 63 S 0.66666666666667 0.83333333333333 0.09514758466182 32.065 > 49 64 S 0.91666666666667 0.83333333333333 0.09514758466182 32.065 > 49 65 Cs 0.16666666666667 0.08333333333333 0.75000000000000 132.905 > 65 66 Cs 0.41666666666667 0.08333333333333 0.75000000000000 132.905 > 65 67 Cs 0.66666666666667 0.08333333333333 0.75000000000000 132.905 > 65 68 Cs 0.91666666666667 0.08333333333333 0.75000000000000 132.905 > 65 69 Cs 0.16666666666667 0.33333333333333 0.75000000000000 132.905 > 65 70 Cs 0.41666666666667 0.33333333333333 0.75000000000000 132.905 > 65 71 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 > 65 72 Cs 0.91666666666667 0.33333333333333 0.75000000000000 132.905 > 65 73 Cs 0.16666666666667 0.58333333333333 0.75000000000000 132.905 > 65 74 Cs 0.41666666666667 0.58333333333333 0.75000000000000 132.905 > 65 75 Cs 0.66666666666667 0.58333333333333 0.75000000000000 132.905 > 65 76 Cs 0.91666666666667 0.58333333333333 0.75000000000000 132.905 > 65 77 Cs 0.16666666666667 0.83333333333333 0.75000000000000 132.905 > 65 78 Cs 0.41666666666667 0.83333333333333 0.75000000000000 132.905 > 65 79 Cs 0.66666666666667 0.83333333333333 0.75000000000000 132.905 > 65 80 Cs 0.91666666666667 0.83333333333333 0.75000000000000 132.905 > 65 81 Cs 0.08333333333333 0.16666666666667 0.25000000000000 132.905 > 81 82 Cs 0.33333333333333 0.16666666666667 0.25000000000000 132.905 > 81 83 Cs 0.58333333333333 0.16666666666667 0.25000000000000 132.905 > 81 84 Cs 0.83333333333333 0.16666666666667 0.25000000000000 132.905 > 81 85 Cs 0.08333333333333 0.41666666666667 0.25000000000000 132.905 > 81 86 Cs 0.33333333333333 0.41666666666667 0.25000000000000 132.905 > 81 87 Cs 0.58333333333333 0.41666666666667 0.25000000000000 132.905 > 81 88 Cs 0.83333333333333 0.41666666666667 0.25000000000000 132.905 > 81 89 Cs 0.08333333333333 0.66666666666667 0.25000000000000 132.905 > 81 90 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 > 81 91 Cs 0.58333333333333 0.66666666666667 0.25000000000000 132.905 > 81 92 Cs 0.83333333333333 0.66666666666667 0.25000000000000 132.905 > 81 93 Cs 0.08333333333333 0.91666666666667 0.25000000000000 132.905 > 81 94 Cs 0.33333333333333 0.91666666666667 0.25000000000000 132.905 > 81 95 Cs 0.58333333333333 0.91666666666667 0.25000000000000 132.905 > 81 96 Cs 0.83333333333333 0.91666666666667 0.25000000000000 132.905 > 81 97 Pr 0.00000000000000 0.00000000000000 0.00000000000000 140.908 > 97 98 Pr 0.25000000000000 0.00000000000000 0.00000000000000 140.908 > 97 99 Pr 0.50000000000000 0.00000000000000 0.00000000000000 140.908 > 97 100 Pr 0.75000000000000 0.00000000000000 0.00000000000000 140.908 > 97 101 Pr 0.00000000000000 0.25000000000000 0.00000000000000 140.908 > 97 102 Pr 0.25000000000000 0.25000000000000 0.00000000000000 140.908 > 97 103 Pr 0.50000000000000 0.25000000000000 0.00000000000000 140.908 > 97 104 Pr 0.75000000000000 0.25000000000000 0.00000000000000 140.908 > 97 105 Pr 0.00000000000000 0.50000000000000 0.00000000000000 140.908 > 97 106 Pr 0.25000000000000 0.50000000000000 0.00000000000000 140.908 > 97 107 Pr 0.50000000000000 0.50000000000000 0.00000000000000 140.908 > 97 108 Pr 0.75000000000000 0.50000000000000 0.00000000000000 140.908 > 97 109 Pr 0.00000000000000 0.75000000000000 0.00000000000000 140.908 > 97 110 Pr 0.25000000000000 0.75000000000000 0.00000000000000 140.908 > 97 111 Pr 0.50000000000000 0.75000000000000 0.00000000000000 140.908 > 97 112 Pr 0.75000000000000 0.75000000000000 0.00000000000000 140.908 > 97 113 Pr 0.00000000000000 0.00000000000000 0.50000000000000 140.908 > 113 114 Pr 0.25000000000000 0.00000000000000 0.50000000000000 140.908 > 113 115 Pr 0.50000000000000 0.00000000000000 0.50000000000000 140.908 > 113 116 Pr 0.75000000000000 0.00000000000000 0.50000000000000 140.908 > 113 117 Pr 0.00000000000000 0.25000000000000 0.50000000000000 140.908 > 113 118 Pr 0.25000000000000 0.25000000000000 0.50000000000000 140.908 > 113 119 Pr 0.50000000000000 0.25000000000000 0.50000000000000 140.908 > 113 120 Pr 0.75000000000000 0.25000000000000 0.50000000000000 140.908 > 113 121 Pr 0.00000000000000 0.50000000000000 0.50000000000000 140.908 > 113 122 Pr 0.25000000000000 0.50000000000000 0.50000000000000 140.908 > 113 123 Pr 0.50000000000000 0.50000000000000 0.50000000000000 140.908 > 113 124 Pr 0.75000000000000 0.50000000000000 0.50000000000000 140.908 > 113 125 Pr 0.00000000000000 0.75000000000000 0.50000000000000 140.908 > 113 126 Pr 0.25000000000000 0.75000000000000 0.50000000000000 140.908 > 113 127 Pr 0.50000000000000 0.75000000000000 0.50000000000000 140.908 > 113 128 Pr 0.75000000000000 0.75000000000000 0.50000000000000 140.908 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.7978134 -0.0000000 0.0000000 0.0000000 5.7978134 0.0000000 0.0000000 0.0000000 4.8403995 -------------------------- Born effective charges -------------------------- 1 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 2 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 3 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 4 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 5 Cs 1.0916780 -0.0000000 0.0000000 0.0000000 1.0916780 0.0000000 0.0000000 0.0000000 2.1481525 6 Cs 1.0916780 -0.0000000 0.0000000 0.0000000 1.0916780 0.0000000 0.0000000 0.0000000 2.1481525 7 Pr 4.7109525 -0.0000000 0.0000000 -0.0000000 4.7109525 0.0000000 0.0000000 0.0000000 2.2247901 8 Pr 4.7109525 -0.0000000 0.0000000 -0.0000000 4.7109525 0.0000000 0.0000000 0.0000000 2.2247901 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 65, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 97, 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.00000002 (xyx) 0.00000002 (xyx) 0.00000002 (xxy) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-07 23:11:00]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-07 23:11:01]------------------------- 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.214008208050170 0.000000000000000 0.000000000000000 b -2.107004104025085 3.649438159927587 0.000000000000000 c 0.000000000000000 0.000000000000000 15.904157939999999 Atomic positions (fractional): 1 S 0.33333333333333 0.66666666666667 0.90485241533818 32.065 2 S 0.66666666666667 0.33333333333333 0.40485241533818 32.065 3 S 0.33333333333333 0.66666666666667 0.59514758466182 32.065 4 S 0.66666666666667 0.33333333333333 0.09514758466182 32.065 5 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 6 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 7 Pr 0.00000000000000 0.00000000000000 0.00000000000000 140.908 8 Pr 0.00000000000000 0.00000000000000 0.50000000000000 140.908 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.856032832200679 0.000000000000000 0.000000000000000 b -8.428016416100339 14.597752639710349 0.000000000000000 c 0.000000000000000 0.000000000000000 15.904157939999999 Atomic positions (fractional): 1 S 0.08333333333333 0.16666666666667 0.90485241533818 32.065 > 1 2 S 0.33333333333333 0.16666666666667 0.90485241533818 32.065 > 1 3 S 0.58333333333333 0.16666666666667 0.90485241533818 32.065 > 1 4 S 0.83333333333333 0.16666666666667 0.90485241533818 32.065 > 1 5 S 0.08333333333333 0.41666666666667 0.90485241533818 32.065 > 1 6 S 0.33333333333333 0.41666666666667 0.90485241533818 32.065 > 1 7 S 0.58333333333333 0.41666666666667 0.90485241533818 32.065 > 1 8 S 0.83333333333333 0.41666666666667 0.90485241533818 32.065 > 1 9 S 0.08333333333333 0.66666666666667 0.90485241533818 32.065 > 1 10 S 0.33333333333333 0.66666666666667 0.90485241533818 32.065 > 1 11 S 0.58333333333333 0.66666666666667 0.90485241533818 32.065 > 1 12 S 0.83333333333333 0.66666666666667 0.90485241533818 32.065 > 1 13 S 0.08333333333333 0.91666666666667 0.90485241533818 32.065 > 1 14 S 0.33333333333333 0.91666666666667 0.90485241533818 32.065 > 1 15 S 0.58333333333333 0.91666666666667 0.90485241533818 32.065 > 1 16 S 0.83333333333333 0.91666666666667 0.90485241533818 32.065 > 1 17 S 0.16666666666667 0.08333333333333 0.40485241533818 32.065 > 17 18 S 0.41666666666667 0.08333333333333 0.40485241533818 32.065 > 17 19 S 0.66666666666667 0.08333333333333 0.40485241533818 32.065 > 17 20 S 0.91666666666667 0.08333333333333 0.40485241533818 32.065 > 17 21 S 0.16666666666667 0.33333333333333 0.40485241533818 32.065 > 17 22 S 0.41666666666667 0.33333333333333 0.40485241533818 32.065 > 17 23 S 0.66666666666667 0.33333333333333 0.40485241533818 32.065 > 17 24 S 0.91666666666667 0.33333333333333 0.40485241533818 32.065 > 17 25 S 0.16666666666667 0.58333333333333 0.40485241533818 32.065 > 17 26 S 0.41666666666667 0.58333333333333 0.40485241533818 32.065 > 17 27 S 0.66666666666667 0.58333333333333 0.40485241533818 32.065 > 17 28 S 0.91666666666667 0.58333333333333 0.40485241533818 32.065 > 17 29 S 0.16666666666667 0.83333333333333 0.40485241533818 32.065 > 17 30 S 0.41666666666667 0.83333333333333 0.40485241533818 32.065 > 17 31 S 0.66666666666667 0.83333333333333 0.40485241533818 32.065 > 17 32 S 0.91666666666667 0.83333333333333 0.40485241533818 32.065 > 17 33 S 0.08333333333333 0.16666666666667 0.59514758466182 32.065 > 33 34 S 0.33333333333333 0.16666666666667 0.59514758466182 32.065 > 33 35 S 0.58333333333333 0.16666666666667 0.59514758466182 32.065 > 33 36 S 0.83333333333333 0.16666666666667 0.59514758466182 32.065 > 33 37 S 0.08333333333333 0.41666666666667 0.59514758466182 32.065 > 33 38 S 0.33333333333333 0.41666666666667 0.59514758466182 32.065 > 33 39 S 0.58333333333333 0.41666666666667 0.59514758466182 32.065 > 33 40 S 0.83333333333333 0.41666666666667 0.59514758466182 32.065 > 33 41 S 0.08333333333333 0.66666666666667 0.59514758466182 32.065 > 33 42 S 0.33333333333333 0.66666666666667 0.59514758466182 32.065 > 33 43 S 0.58333333333333 0.66666666666667 0.59514758466182 32.065 > 33 44 S 0.83333333333333 0.66666666666667 0.59514758466182 32.065 > 33 45 S 0.08333333333333 0.91666666666667 0.59514758466182 32.065 > 33 46 S 0.33333333333333 0.91666666666667 0.59514758466182 32.065 > 33 47 S 0.58333333333333 0.91666666666667 0.59514758466182 32.065 > 33 48 S 0.83333333333333 0.91666666666667 0.59514758466182 32.065 > 33 49 S 0.16666666666667 0.08333333333333 0.09514758466182 32.065 > 49 50 S 0.41666666666667 0.08333333333333 0.09514758466182 32.065 > 49 51 S 0.66666666666667 0.08333333333333 0.09514758466182 32.065 > 49 52 S 0.91666666666667 0.08333333333333 0.09514758466182 32.065 > 49 53 S 0.16666666666667 0.33333333333333 0.09514758466182 32.065 > 49 54 S 0.41666666666667 0.33333333333333 0.09514758466182 32.065 > 49 55 S 0.66666666666667 0.33333333333333 0.09514758466182 32.065 > 49 56 S 0.91666666666667 0.33333333333333 0.09514758466182 32.065 > 49 57 S 0.16666666666667 0.58333333333333 0.09514758466182 32.065 > 49 58 S 0.41666666666667 0.58333333333333 0.09514758466182 32.065 > 49 59 S 0.66666666666667 0.58333333333333 0.09514758466182 32.065 > 49 60 S 0.91666666666667 0.58333333333333 0.09514758466182 32.065 > 49 61 S 0.16666666666667 0.83333333333333 0.09514758466182 32.065 > 49 62 S 0.41666666666667 0.83333333333333 0.09514758466182 32.065 > 49 63 S 0.66666666666667 0.83333333333333 0.09514758466182 32.065 > 49 64 S 0.91666666666667 0.83333333333333 0.09514758466182 32.065 > 49 65 Cs 0.16666666666667 0.08333333333333 0.75000000000000 132.905 > 65 66 Cs 0.41666666666667 0.08333333333333 0.75000000000000 132.905 > 65 67 Cs 0.66666666666667 0.08333333333333 0.75000000000000 132.905 > 65 68 Cs 0.91666666666667 0.08333333333333 0.75000000000000 132.905 > 65 69 Cs 0.16666666666667 0.33333333333333 0.75000000000000 132.905 > 65 70 Cs 0.41666666666667 0.33333333333333 0.75000000000000 132.905 > 65 71 Cs 0.66666666666667 0.33333333333333 0.75000000000000 132.905 > 65 72 Cs 0.91666666666667 0.33333333333333 0.75000000000000 132.905 > 65 73 Cs 0.16666666666667 0.58333333333333 0.75000000000000 132.905 > 65 74 Cs 0.41666666666667 0.58333333333333 0.75000000000000 132.905 > 65 75 Cs 0.66666666666667 0.58333333333333 0.75000000000000 132.905 > 65 76 Cs 0.91666666666667 0.58333333333333 0.75000000000000 132.905 > 65 77 Cs 0.16666666666667 0.83333333333333 0.75000000000000 132.905 > 65 78 Cs 0.41666666666667 0.83333333333333 0.75000000000000 132.905 > 65 79 Cs 0.66666666666667 0.83333333333333 0.75000000000000 132.905 > 65 80 Cs 0.91666666666667 0.83333333333333 0.75000000000000 132.905 > 65 81 Cs 0.08333333333333 0.16666666666667 0.25000000000000 132.905 > 81 82 Cs 0.33333333333333 0.16666666666667 0.25000000000000 132.905 > 81 83 Cs 0.58333333333333 0.16666666666667 0.25000000000000 132.905 > 81 84 Cs 0.83333333333333 0.16666666666667 0.25000000000000 132.905 > 81 85 Cs 0.08333333333333 0.41666666666667 0.25000000000000 132.905 > 81 86 Cs 0.33333333333333 0.41666666666667 0.25000000000000 132.905 > 81 87 Cs 0.58333333333333 0.41666666666667 0.25000000000000 132.905 > 81 88 Cs 0.83333333333333 0.41666666666667 0.25000000000000 132.905 > 81 89 Cs 0.08333333333333 0.66666666666667 0.25000000000000 132.905 > 81 90 Cs 0.33333333333333 0.66666666666667 0.25000000000000 132.905 > 81 91 Cs 0.58333333333333 0.66666666666667 0.25000000000000 132.905 > 81 92 Cs 0.83333333333333 0.66666666666667 0.25000000000000 132.905 > 81 93 Cs 0.08333333333333 0.91666666666667 0.25000000000000 132.905 > 81 94 Cs 0.33333333333333 0.91666666666667 0.25000000000000 132.905 > 81 95 Cs 0.58333333333333 0.91666666666667 0.25000000000000 132.905 > 81 96 Cs 0.83333333333333 0.91666666666667 0.25000000000000 132.905 > 81 97 Pr 0.00000000000000 0.00000000000000 0.00000000000000 140.908 > 97 98 Pr 0.25000000000000 0.00000000000000 0.00000000000000 140.908 > 97 99 Pr 0.50000000000000 0.00000000000000 0.00000000000000 140.908 > 97 100 Pr 0.75000000000000 0.00000000000000 0.00000000000000 140.908 > 97 101 Pr 0.00000000000000 0.25000000000000 0.00000000000000 140.908 > 97 102 Pr 0.25000000000000 0.25000000000000 0.00000000000000 140.908 > 97 103 Pr 0.50000000000000 0.25000000000000 0.00000000000000 140.908 > 97 104 Pr 0.75000000000000 0.25000000000000 0.00000000000000 140.908 > 97 105 Pr 0.00000000000000 0.50000000000000 0.00000000000000 140.908 > 97 106 Pr 0.25000000000000 0.50000000000000 0.00000000000000 140.908 > 97 107 Pr 0.50000000000000 0.50000000000000 0.00000000000000 140.908 > 97 108 Pr 0.75000000000000 0.50000000000000 0.00000000000000 140.908 > 97 109 Pr 0.00000000000000 0.75000000000000 0.00000000000000 140.908 > 97 110 Pr 0.25000000000000 0.75000000000000 0.00000000000000 140.908 > 97 111 Pr 0.50000000000000 0.75000000000000 0.00000000000000 140.908 > 97 112 Pr 0.75000000000000 0.75000000000000 0.00000000000000 140.908 > 97 113 Pr 0.00000000000000 0.00000000000000 0.50000000000000 140.908 > 113 114 Pr 0.25000000000000 0.00000000000000 0.50000000000000 140.908 > 113 115 Pr 0.50000000000000 0.00000000000000 0.50000000000000 140.908 > 113 116 Pr 0.75000000000000 0.00000000000000 0.50000000000000 140.908 > 113 117 Pr 0.00000000000000 0.25000000000000 0.50000000000000 140.908 > 113 118 Pr 0.25000000000000 0.25000000000000 0.50000000000000 140.908 > 113 119 Pr 0.50000000000000 0.25000000000000 0.50000000000000 140.908 > 113 120 Pr 0.75000000000000 0.25000000000000 0.50000000000000 140.908 > 113 121 Pr 0.00000000000000 0.50000000000000 0.50000000000000 140.908 > 113 122 Pr 0.25000000000000 0.50000000000000 0.50000000000000 140.908 > 113 123 Pr 0.50000000000000 0.50000000000000 0.50000000000000 140.908 > 113 124 Pr 0.75000000000000 0.50000000000000 0.50000000000000 140.908 > 113 125 Pr 0.00000000000000 0.75000000000000 0.50000000000000 140.908 > 113 126 Pr 0.25000000000000 0.75000000000000 0.50000000000000 140.908 > 113 127 Pr 0.50000000000000 0.75000000000000 0.50000000000000 140.908 > 113 128 Pr 0.75000000000000 0.75000000000000 0.50000000000000 140.908 > 113 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 5.7978134 -0.0000000 0.0000000 0.0000000 5.7978134 0.0000000 0.0000000 0.0000000 4.8403995 -------------------------- Born effective charges -------------------------- 1 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 2 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 3 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 4 S -2.9013152 0.0000000 0.0000000 0.0000000 -2.9013152 0.0000000 0.0000000 0.0000000 -2.1864713 5 Cs 1.0916780 -0.0000000 0.0000000 0.0000000 1.0916780 0.0000000 0.0000000 0.0000000 2.1481525 6 Cs 1.0916780 -0.0000000 0.0000000 0.0000000 1.0916780 0.0000000 0.0000000 0.0000000 2.1481525 7 Pr 4.7109525 -0.0000000 0.0000000 -0.0000000 4.7109525 0.0000000 0.0000000 0.0000000 2.2247901 8 Pr 4.7109525 -0.0000000 0.0000000 -0.0000000 4.7109525 0.0000000 0.0000000 0.0000000 2.2247901 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000002 (xyx) 0.00000002 (xyx) 0.00000002 (xxy) Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) --------------------------- 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.66, Number of G-points: 291, Lambda: 0.16 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: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 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.129 ( 0.000 0.000 0.000) 0.000 1.129 ( 0.000 0.000 0.000) 0.000 1.374 ( 0.000 0.000 0.000) 0.000 1.374 ( 0.000 0.000 0.000) 0.000 1.778 ( 0.000 0.000 0.000) 0.000 1.778 ( 0.000 0.000 0.000) 0.000 2.202 ( 0.000 0.000 0.000) 0.000 2.558 ( 0.000 0.000 0.000) 0.000 3.001 ( 0.000 0.000 0.000) 0.000 4.699 ( 0.000 0.000 0.000) 0.000 4.699 ( 0.000 0.000 0.000) 0.000 4.716 ( 0.000 0.000 0.000) 0.000 4.716 ( 0.000 0.000 0.000) 0.000 6.041 ( 0.000 0.000 0.000) 0.000 6.041 ( 0.000 0.000 0.000) 0.000 6.063 ( 0.000 0.000 0.000) 0.000 6.063 ( 0.000 0.000 0.000) 0.000 6.397 ( 0.000 0.000 0.000) 0.000 7.535 ( 0.000 0.000 0.000) 0.000 7.656 ( 0.000 0.000 0.000) 0.000 8.037 ( 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: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.423 ( 18.467 10.662 0.000) 21.324 0.453 ( 18.782 10.844 0.000) 21.688 0.780 ( 32.016 18.484 0.000) 36.969 1.245 ( 9.371 5.410 0.000) 10.821 1.352 ( 16.549 9.555 0.000) 19.110 1.373 ( -0.633 -0.366 0.000) 0.731 1.381 ( 0.608 0.351 0.000) 0.703 1.800 ( 1.916 1.106 0.000) 2.212 1.868 ( 8.326 4.807 0.000) 9.614 2.214 ( 2.005 1.158 0.000) 2.316 2.511 ( -1.323 -0.764 0.000) 1.528 2.986 ( -1.396 -0.806 0.000) 1.612 4.746 ( 4.023 2.323 0.000) 4.646 4.760 ( 3.758 2.170 0.000) 4.339 5.153 ( 32.527 18.780 0.000) 37.559 5.997 ( -3.787 -2.187 0.000) 4.373 6.015 ( -2.159 -1.246 0.000) 2.493 6.019 ( -3.719 -2.147 0.000) 4.294 6.037 ( -2.251 -1.300 0.000) 2.599 6.385 ( -0.961 -0.555 0.000) 1.109 7.463 ( -6.131 -3.540 0.000) 7.080 7.633 ( -1.606 -0.927 0.000) 1.855 7.895 ( -12.001 -6.929 0.000) 13.857 8.090 ( 0.236 0.136 0.000) 0.272 ======================= Grid point 2 (3/48) ======================= q-point: ( 0.14 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.828 ( 14.070 8.123 0.000) 16.247 0.830 ( 17.395 10.043 0.000) 20.087 1.319 ( 1.813 1.047 0.000) 2.093 1.396 ( 16.026 9.253 0.000) 18.505 1.401 ( 1.106 0.638 0.000) 1.277 1.502 ( 12.422 7.172 0.000) 14.343 1.713 ( 12.616 7.284 0.000) 14.568 1.867 ( 4.023 2.322 0.000) 4.645 2.185 ( 19.955 11.521 0.000) 23.042 2.322 ( 7.343 4.240 0.000) 8.479 2.600 ( 10.286 5.938 0.000) 11.877 2.937 ( -2.966 -1.713 0.000) 3.425 4.867 ( 6.205 3.583 0.000) 7.165 4.874 ( 5.847 3.376 0.000) 6.752 5.826 ( 20.882 12.056 0.000) 24.113 5.879 ( -4.829 -2.788 0.000) 5.576 5.910 ( -5.511 -3.182 0.000) 6.364 5.952 ( -3.128 -1.806 0.000) 3.612 5.970 ( -3.350 -1.934 0.000) 3.868 6.362 ( -0.839 -0.484 0.000) 0.969 7.276 ( -10.042 -5.798 0.000) 11.596 7.539 ( -18.430 -10.641 0.000) 21.281 7.618 ( 1.029 0.594 0.000) 1.188 8.085 ( -0.849 -0.490 0.000) 0.981 ======================= Grid point 3 (4/48) ======================= q-point: ( 0.21 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.088 ( 9.200 5.311 0.000) 10.623 1.201 ( 15.225 8.790 0.000) 17.580 1.314 ( 1.324 0.764 0.000) 1.529 1.429 ( 1.277 0.737 0.000) 1.475 1.706 ( 9.807 5.662 0.000) 11.324 1.778 ( 11.627 6.713 0.000) 13.426 1.889 ( 4.027 2.325 0.000) 4.649 1.980 ( 5.840 3.372 0.000) 6.743 2.507 ( 8.039 4.641 0.000) 9.282 2.719 ( 25.407 14.669 0.000) 29.338 2.837 ( -3.651 -2.108 0.000) 4.215 2.962 ( 17.630 10.179 0.000) 20.357 4.999 ( 4.683 2.704 0.000) 5.408 4.999 ( 4.954 2.860 0.000) 5.720 5.716 ( -6.972 -4.025 0.000) 8.051 5.785 ( -5.236 -3.023 0.000) 6.046 5.891 ( -1.857 -1.072 0.000) 2.144 5.903 ( -2.254 -1.302 0.000) 2.603 6.118 ( 4.401 2.541 0.000) 5.082 6.354 ( -0.135 -0.078 0.000) 0.156 7.019 ( -12.612 -7.282 0.000) 14.563 7.111 ( -18.284 -10.556 0.000) 21.112 7.684 ( 4.004 2.312 0.000) 4.623 8.048 ( -2.470 -1.426 0.000) 2.852 ======================= Grid point 4 (5/48) ======================= q-point: ( 0.29 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.257 ( 5.993 3.460 0.000) 6.920 1.399 ( 6.007 3.468 0.000) 6.936 1.456 ( 1.113 0.643 0.000) 1.286 1.511 ( 12.170 7.026 0.000) 14.053 1.872 ( 5.297 3.058 0.000) 6.116 1.942 ( 1.222 0.705 0.000) 1.411 2.020 ( 9.687 5.593 0.000) 11.186 2.121 ( 6.445 3.721 0.000) 7.442 2.648 ( 3.934 2.271 0.000) 4.542 2.710 ( -6.683 -3.859 0.000) 7.717 3.276 ( 23.127 13.352 0.000) 26.704 3.389 ( 18.331 10.583 0.000) 21.167 5.073 ( 1.901 1.097 0.000) 2.195 5.078 ( 2.104 1.215 0.000) 2.430 5.586 ( -4.013 -2.317 0.000) 4.634 5.693 ( -2.487 -1.436 0.000) 2.872 5.874 ( 0.116 0.067 0.000) 0.134 5.875 ( -0.373 -0.215 0.000) 0.431 6.135 ( -1.370 -0.791 0.000) 1.581 6.338 ( -1.268 -0.732 0.000) 1.465 6.697 ( -14.637 -8.451 0.000) 16.902 6.766 ( -12.736 -7.353 0.000) 14.706 7.749 ( 0.933 0.539 0.000) 1.077 7.971 ( -4.361 -2.518 0.000) 5.035 ======================= Grid point 5 (6/48) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.367 ( 3.796 2.192 0.000) 4.383 1.479 ( 0.970 0.560 0.000) 1.120 1.561 ( 7.541 4.354 0.000) 8.708 1.750 ( 9.001 5.197 0.000) 10.393 1.959 ( 0.444 0.256 0.000) 0.512 1.962 ( 3.028 1.748 0.000) 3.496 2.212 ( 7.227 4.173 0.000) 8.345 2.260 ( 5.615 3.242 0.000) 6.484 2.545 ( -7.652 -4.418 0.000) 8.836 2.674 ( -1.539 -0.889 0.000) 1.777 3.740 ( 17.415 10.054 0.000) 20.109 3.759 ( 13.978 8.070 0.000) 16.140 5.100 ( 0.873 0.504 0.000) 1.008 5.110 ( 1.134 0.655 0.000) 1.310 5.553 ( 1.237 0.714 0.000) 1.428 5.690 ( 2.535 1.464 0.000) 2.927 5.871 ( -0.310 -0.179 0.000) 0.358 5.880 ( 0.076 0.044 0.000) 0.088 6.089 ( -2.446 -1.412 0.000) 2.825 6.309 ( -1.073 -0.620 0.000) 1.239 6.411 ( -10.711 -6.184 0.000) 12.368 6.507 ( -9.889 -5.709 0.000) 11.419 7.717 ( -3.411 -1.969 0.000) 3.939 7.854 ( -5.813 -3.356 0.000) 6.712 ======================= Grid point 6 (7/48) ======================= q-point: ( 0.43 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.428 ( 1.596 0.921 0.000) 1.843 1.501 ( 0.897 0.518 0.000) 1.036 1.704 ( 4.452 2.571 0.000) 5.141 1.911 ( 4.418 2.551 0.000) 5.102 1.976 ( 1.947 1.124 0.000) 2.248 2.026 ( 2.853 1.647 0.000) 3.294 2.341 ( 3.989 2.303 0.000) 4.606 2.364 ( 3.373 1.947 0.000) 3.895 2.378 ( -6.650 -3.840 0.000) 7.679 2.599 ( -4.218 -2.436 0.000) 4.871 4.004 ( 7.439 4.295 0.000) 8.590 4.043 ( 9.098 5.253 0.000) 10.506 5.120 ( 0.820 0.473 0.000) 0.947 5.136 ( 1.052 0.607 0.000) 1.215 5.630 ( 4.874 2.814 0.000) 5.628 5.806 ( 7.293 4.211 0.000) 8.422 5.857 ( -0.770 -0.445 0.000) 0.890 5.873 ( -0.577 -0.333 0.000) 0.667 6.025 ( -3.185 -1.839 0.000) 3.678 6.196 ( -8.757 -5.056 0.000) 10.112 6.292 ( -0.417 -0.241 0.000) 0.482 6.335 ( -5.082 -2.934 0.000) 5.868 7.623 ( -4.052 -2.340 0.000) 4.679 7.728 ( -4.749 -2.742 0.000) 5.484 ======================= Grid point 7 (8/48) ======================= q-point: (-0.50 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 64 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.444 ( 0.000 0.000 0.000) 0.000 1.513 ( 0.000 0.000 0.000) 0.000 1.752 ( 0.000 0.000 0.000) 0.000 1.943 ( 0.000 0.000 0.000) 0.000 2.031 ( -0.000 -0.000 0.000) 0.000 2.069 ( 0.000 0.000 0.000) 0.000 2.286 ( 0.000 0.000 0.000) 0.000 2.387 ( 0.000 0.000 0.000) 0.000 2.404 ( 0.000 0.000 0.000) 0.000 2.538 ( 0.000 0.000 0.000) 0.000 4.089 ( 0.000 0.000 0.000) 0.000 4.145 ( 0.000 0.000 0.000) 0.000 5.130 ( 0.000 0.000 0.000) 0.000 5.150 ( 0.000 0.000 0.000) 0.000 5.707 ( 0.000 0.000 0.000) 0.000 5.846 ( 0.000 0.000 0.000) 0.000 5.864 ( 0.000 0.000 0.000) 0.000 5.974 ( 0.000 0.000 0.000) 0.000 5.976 ( 0.000 0.000 0.000) 0.000 6.017 ( -0.000 -0.000 0.000) 0.000 6.279 ( 0.000 0.000 0.000) 0.000 6.288 ( 0.000 0.000 0.000) 0.000 7.571 ( 0.000 0.000 0.000) 0.000 7.669 ( 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: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.721 ( 10.177 17.627 0.000) 20.354 0.768 ( 10.110 17.512 0.000) 20.221 1.219 ( 11.569 20.038 0.000) 23.137 1.336 ( -1.652 -2.862 0.000) 3.305 1.423 ( 2.678 4.639 0.000) 5.357 1.445 ( 7.814 13.534 0.000) 15.627 1.619 ( 8.470 14.670 0.000) 16.939 1.855 ( 2.530 4.382 0.000) 5.059 2.063 ( 9.150 15.849 0.000) 18.301 2.277 ( 3.325 5.759 0.000) 6.650 2.546 ( 3.614 6.260 0.000) 7.229 2.953 ( -1.471 -2.548 0.000) 2.942 4.839 ( 3.961 6.860 0.000) 7.921 4.846 ( 3.650 6.322 0.000) 7.300 5.674 ( 15.984 27.685 0.000) 31.968 5.913 ( -3.543 -6.136 0.000) 7.085 5.942 ( -3.086 -5.345 0.000) 6.172 5.972 ( -1.703 -2.949 0.000) 3.406 5.990 ( -1.896 -3.284 0.000) 3.792 6.366 ( -0.740 -1.282 0.000) 1.481 7.334 ( -5.297 -9.174 0.000) 10.594 7.613 ( -0.243 -0.422 0.000) 0.487 7.657 ( -9.574 -16.582 0.000) 19.147 8.077 ( -0.843 -1.459 0.000) 1.685 ======================= Grid point 17 (10/48) ======================= q-point: ( 0.14 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.065 ( 6.140 13.714 0.000) 15.025 1.071 ( 12.332 13.863 0.000) 18.554 1.289 ( 0.031 -1.810 0.000) 1.810 1.411 ( 3.276 2.002 0.000) 3.839 1.625 ( 7.826 10.087 0.000) 12.767 1.724 ( 8.023 14.290 0.000) 16.388 1.839 ( 5.606 5.012 0.000) 7.520 1.964 ( 3.461 7.464 0.000) 8.227 2.429 ( 6.623 6.401 0.000) 9.211 2.491 ( 19.436 19.271 0.000) 27.370 2.789 ( 12.432 13.060 0.000) 18.031 2.887 ( -2.430 -3.243 0.000) 4.052 4.989 ( 4.055 8.120 0.000) 9.077 4.992 ( 4.157 8.855 0.000) 9.782 5.771 ( -5.399 -6.499 0.000) 8.449 5.824 ( -4.225 -5.368 0.000) 6.831 5.910 ( -2.106 -2.325 0.000) 3.138 5.921 ( -2.228 -2.921 0.000) 3.673 6.070 ( 7.560 9.478 0.000) 12.124 6.344 ( -0.020 -1.360 0.000) 1.360 7.118 ( -8.797 -10.033 0.000) 13.344 7.285 ( -15.094 -14.548 0.000) 20.964 7.638 ( 2.649 1.370 0.000) 2.982 8.034 ( -0.754 -4.656 0.000) 4.717 ======================= Grid point 18 (11/48) ======================= q-point: ( 0.21 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.240 ( 3.700 6.482 0.000) 7.464 1.330 ( 3.013 4.135 0.000) 5.116 1.395 ( 11.341 10.685 0.000) 15.582 1.478 ( 1.742 3.959 0.000) 4.326 1.807 ( 5.647 4.269 0.000) 7.079 1.915 ( 2.384 2.062 0.000) 3.152 1.997 ( 6.013 13.152 0.000) 14.461 2.113 ( 3.646 10.081 0.000) 10.720 2.589 ( 5.638 3.377 0.000) 6.572 2.773 ( -5.199 -4.030 0.000) 6.578 3.032 ( 22.006 16.875 0.000) 27.731 3.196 ( 17.256 13.111 0.000) 21.672 5.116 ( 0.517 9.109 0.000) 9.124 5.125 ( 0.324 9.621 0.000) 9.626 5.625 ( -4.544 -4.887 0.000) 6.674 5.710 ( -3.072 -4.492 0.000) 5.442 5.869 ( -0.426 -1.278 0.000) 1.347 5.874 ( -0.472 -1.478 0.000) 1.552 6.155 ( -0.960 1.049 0.000) 1.422 6.335 ( 0.199 -1.708 0.000) 1.719 6.846 ( -11.797 -9.967 0.000) 15.444 6.903 ( -14.282 -10.506 0.000) 17.730 7.701 ( 3.359 -1.058 0.000) 3.522 7.962 ( -0.911 -7.172 0.000) 7.229 ======================= Grid point 19 (12/48) ======================= q-point: ( 0.29 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.361 ( 4.179 4.698 0.000) 6.288 1.424 ( 2.693 3.627 0.000) 4.517 1.559 ( 1.873 6.839 0.000) 7.091 1.659 ( 9.036 7.799 0.000) 11.936 1.911 ( 3.886 1.086 0.000) 4.035 1.949 ( 0.875 0.402 0.000) 0.963 2.212 ( 3.808 12.901 0.000) 13.451 2.274 ( 2.487 11.576 0.000) 11.840 2.630 ( -6.881 -3.813 0.000) 7.867 2.666 ( 1.468 -0.380 0.000) 1.516 3.534 ( 19.387 12.784 0.000) 23.222 3.596 ( 15.535 10.462 0.000) 18.730 5.183 ( -3.178 9.806 0.000) 10.308 5.191 ( -3.271 9.721 0.000) 10.256 5.543 ( -0.779 -1.834 0.000) 1.992 5.651 ( 0.854 -2.815 0.000) 2.942 5.864 ( 0.832 -1.204 0.000) 1.464 5.871 ( 0.514 0.086 0.000) 0.521 6.120 ( -2.560 -0.669 0.000) 2.646 6.314 ( -0.567 -1.465 0.000) 1.571 6.545 ( -12.320 -7.184 0.000) 14.262 6.607 ( -10.348 -8.457 0.000) 13.364 7.707 ( 0.497 -4.668 0.000) 4.695 7.863 ( -1.894 -8.097 0.000) 8.315 ======================= Grid point 20 (13/48) ======================= q-point: ( 0.36 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.455 ( 0.443 5.747 0.000) 5.764 1.485 ( 1.743 1.624 0.000) 2.382 1.674 ( 3.554 5.956 0.000) 6.936 1.854 ( 6.450 5.228 0.000) 8.303 1.962 ( 0.501 0.247 0.000) 0.558 1.973 ( 3.397 -0.325 0.000) 3.413 2.359 ( 0.352 6.881 0.000) 6.890 2.412 ( 0.257 11.484 0.000) 11.487 2.487 ( -7.067 2.140 0.000) 7.383 2.635 ( -3.000 -2.370 0.000) 3.823 3.905 ( 10.475 6.662 0.000) 12.414 3.915 ( 13.123 7.677 0.000) 15.203 5.216 ( -4.500 10.610 0.000) 11.524 5.221 ( -4.158 9.765 0.000) 10.614 5.568 ( 4.051 0.761 0.000) 4.122 5.693 ( 6.326 -1.170 0.000) 6.433 5.861 ( 0.513 -1.996 0.000) 2.061 5.880 ( -0.351 1.284 0.000) 1.331 6.063 ( -3.158 -1.109 0.000) 3.347 6.294 ( -0.359 -0.847 0.000) 0.920 6.299 ( -9.624 -5.756 0.000) 11.214 6.390 ( -6.651 -5.731 0.000) 8.780 7.638 ( -1.984 -5.614 0.000) 5.954 7.746 ( -2.667 -7.010 0.000) 7.500 ======================= Grid point 21 (14/48) ======================= q-point: ( 0.43 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.491 ( -2.045 5.014 0.000) 5.415 1.519 ( 0.226 1.525 0.000) 1.542 1.760 ( 0.856 2.970 0.000) 3.091 1.948 ( 0.265 1.283 0.000) 1.310 2.002 ( 3.336 1.180 0.000) 3.538 2.026 ( 3.211 -1.360 0.000) 3.487 2.310 ( -3.444 -2.095 0.000) 4.032 2.494 ( -4.885 12.232 0.000) 13.172 2.498 ( -1.391 7.998 0.000) 8.118 2.551 ( -4.757 0.951 0.000) 4.851 4.070 ( 3.632 2.205 0.000) 4.249 4.115 ( 4.713 2.062 0.000) 5.144 5.238 ( -5.382 10.728 0.000) 12.002 5.241 ( -4.613 9.290 0.000) 10.373 5.661 ( 4.697 -0.113 0.000) 4.698 5.796 ( 5.113 -2.657 0.000) 5.762 5.849 ( 0.661 -1.917 0.000) 2.028 5.918 ( 4.083 4.165 0.000) 5.833 5.996 ( -2.910 -0.761 0.000) 3.008 6.093 ( -8.259 -5.659 0.000) 10.012 6.279 ( -1.751 -2.460 0.000) 3.020 6.286 ( 0.050 -0.377 0.000) 0.380 7.561 ( -0.769 -3.573 0.000) 3.654 7.657 ( -0.736 -4.106 0.000) 4.171 ======================= Grid point 31 (15/48) ======================= q-point: ( 0.14 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.228 ( 1.464 2.536 0.000) 2.928 1.349 ( 8.087 14.007 0.000) 16.173 1.363 ( 5.422 9.392 0.000) 10.845 1.464 ( 1.977 3.424 0.000) 3.954 1.781 ( 3.504 6.069 0.000) 7.008 1.907 ( 1.423 2.465 0.000) 2.847 2.018 ( 8.488 14.702 0.000) 16.976 2.143 ( 6.021 10.429 0.000) 12.042 2.552 ( 3.266 5.657 0.000) 6.532 2.800 ( -2.904 -5.030 0.000) 5.808 2.915 ( 13.379 23.173 0.000) 26.758 3.104 ( 10.345 17.917 0.000) 20.689 5.168 ( 5.431 9.406 0.000) 10.861 5.183 ( 5.583 9.669 0.000) 11.165 5.640 ( -3.653 -6.328 0.000) 7.307 5.713 ( -3.254 -5.636 0.000) 6.508 5.866 ( -1.131 -1.959 0.000) 2.262 5.869 ( -1.148 -1.989 0.000) 2.296 6.176 ( 0.924 1.600 0.000) 1.847 6.316 ( -0.713 -1.235 0.000) 1.426 6.894 ( -7.273 -12.598 0.000) 14.547 6.991 ( -8.707 -15.082 0.000) 17.415 7.659 ( 0.183 0.318 0.000) 0.367 7.919 ( -3.788 -6.561 0.000) 7.577 ======================= Grid point 32 (16/48) ======================= q-point: ( 0.21 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.301 ( 4.460 3.272 0.000) 5.531 1.495 ( -0.644 8.639 0.000) 8.663 1.567 ( 4.851 5.048 0.000) 7.001 1.604 ( 7.663 10.331 0.000) 12.862 1.871 ( 2.512 2.314 0.000) 3.416 1.942 ( 0.833 1.100 0.000) 1.379 2.280 ( 5.193 13.998 0.000) 14.930 2.353 ( 4.230 13.429 0.000) 14.079 2.637 ( 2.246 1.947 0.000) 2.972 2.695 ( -4.253 -3.087 0.000) 5.256 3.368 ( 16.014 17.171 0.000) 23.480 3.465 ( 12.960 14.024 0.000) 19.095 5.325 ( 0.711 11.025 0.000) 11.048 5.335 ( 0.261 10.236 0.000) 10.240 5.545 ( -1.881 -2.576 0.000) 3.190 5.619 ( -1.520 -4.392 0.000) 4.648 5.838 ( 0.038 -1.896 0.000) 1.896 5.856 ( 0.921 0.053 0.000) 0.922 6.163 ( -2.309 0.053 0.000) 2.309 6.298 ( 0.212 -1.943 0.000) 1.955 6.635 ( -8.870 -11.685 0.000) 14.671 6.711 ( -10.254 -8.811 0.000) 13.519 7.643 ( 0.340 -4.493 0.000) 4.506 7.792 ( -1.703 -9.305 0.000) 9.460 ======================= Grid point 33 (17/48) ======================= q-point: ( 0.29 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 200 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.427 ( 5.212 4.139 0.000) 6.655 1.552 ( -1.325 6.979 0.000) 7.103 1.687 ( 2.824 6.003 0.000) 6.634 1.805 ( 6.091 6.919 0.000) 9.218 1.921 ( 2.289 -0.006 0.000) 2.289 1.961 ( 0.347 0.742 0.000) 0.819 2.416 ( -1.008 2.477 0.000) 2.674 2.510 ( 0.734 6.842 0.000) 6.881 2.649 ( -3.235 10.800 0.000) 11.274 2.684 ( -0.588 7.598 0.000) 7.621 3.766 ( 14.113 10.549 0.000) 17.620 3.790 ( 11.443 8.890 0.000) 14.491 5.400 ( -2.857 9.049 0.000) 9.490 5.409 ( -3.220 11.368 0.000) 11.815 5.532 ( 0.436 1.860 0.000) 1.911 5.593 ( 2.719 -2.396 0.000) 3.624 5.823 ( 0.759 -2.944 0.000) 3.040 5.888 ( 1.087 1.805 0.000) 2.107 6.117 ( -3.474 0.719 0.000) 3.548 6.283 ( 0.293 -1.673 0.000) 1.699 6.381 ( -8.231 -9.842 0.000) 12.830 6.477 ( -9.537 -4.414 0.000) 10.508 7.588 ( -0.034 -6.639 0.000) 6.639 7.684 ( -0.256 -9.272 0.000) 9.276 ======================= Grid point 34 (18/48) ======================= q-point: ( 0.36 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.530 ( 1.709 4.076 0.000) 4.420 1.588 ( -2.155 6.259 0.000) 6.620 1.778 ( 1.519 4.432 0.000) 4.685 1.944 ( 3.683 3.697 0.000) 5.219 1.959 ( 2.939 -1.001 0.000) 3.104 1.971 ( 0.467 0.641 0.000) 0.794 2.349 ( -3.699 -2.735 0.000) 4.601 2.503 ( -1.004 -4.101 0.000) 4.222 2.745 ( -5.241 18.692 0.000) 19.413 2.776 ( -5.251 18.185 0.000) 18.928 4.007 ( 6.817 3.189 0.000) 7.527 4.030 ( 8.389 3.733 0.000) 9.183 5.433 ( -4.236 9.960 0.000) 10.824 5.454 ( -4.576 11.392 0.000) 12.276 5.582 ( 3.350 0.795 0.000) 3.444 5.659 ( 6.615 -1.867 0.000) 6.874 5.807 ( 1.131 -3.213 0.000) 3.406 5.918 ( 0.465 2.265 0.000) 2.313 6.058 ( -4.363 1.036 0.000) 4.484 6.169 ( -7.421 -6.489 0.000) 9.858 6.277 ( 0.405 -0.895 0.000) 0.982 6.312 ( -5.633 -2.075 0.000) 6.003 7.525 ( -0.117 -5.244 0.000) 5.245 7.606 ( 0.464 -6.799 0.000) 6.815 ======================= Grid point 35 (19/48) ======================= q-point: ( 0.43 0.14 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.564 ( -1.743 3.019 0.000) 3.486 1.601 ( -3.281 5.684 0.000) 6.563 1.813 ( -1.310 2.270 0.000) 2.621 1.969 ( -0.539 0.933 0.000) 1.078 1.986 ( 1.442 -2.498 0.000) 2.884 2.009 ( 0.273 -0.472 0.000) 0.545 2.295 ( -0.235 0.408 0.000) 0.471 2.463 ( 2.044 -3.541 0.000) 4.088 2.796 ( -10.154 17.588 0.000) 20.309 2.824 ( -10.412 18.034 0.000) 20.824 4.080 ( 0.853 -1.477 0.000) 1.705 4.119 ( 1.028 -1.781 0.000) 2.056 5.447 ( -5.824 10.088 0.000) 11.648 5.472 ( -6.331 10.966 0.000) 12.662 5.623 ( 1.844 -3.193 0.000) 3.687 5.721 ( 3.062 -5.304 0.000) 6.124 5.803 ( 1.594 -2.760 0.000) 3.187 5.968 ( 2.164 -3.749 0.000) 4.329 6.016 ( -1.898 3.287 0.000) 3.795 6.026 ( -2.844 4.926 0.000) 5.688 6.255 ( -0.128 0.223 0.000) 0.257 6.277 ( 0.354 -0.614 0.000) 0.709 7.496 ( 1.699 -2.943 0.000) 3.398 7.575 ( 2.448 -4.240 0.000) 4.896 ======================= Grid point 47 (20/48) ======================= q-point: ( 0.21 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.400 ( 3.797 6.577 0.000) 7.594 1.612 ( 2.184 3.782 0.000) 4.368 1.684 ( 3.610 6.253 0.000) 7.220 1.787 ( 4.713 8.163 0.000) 9.426 1.901 ( 0.546 0.946 0.000) 1.092 1.962 ( 0.528 0.914 0.000) 1.056 2.436 ( -0.039 -0.068 0.000) 0.079 2.528 ( 0.769 1.332 0.000) 1.538 2.743 ( 5.646 9.779 0.000) 11.292 2.747 ( 6.796 11.771 0.000) 13.592 3.687 ( 8.605 14.905 0.000) 17.211 3.727 ( 7.073 12.250 0.000) 14.145 5.446 ( 0.807 1.398 0.000) 1.615 5.487 ( 1.980 3.429 0.000) 3.960 5.569 ( 0.656 1.135 0.000) 1.311 5.572 ( 2.661 4.609 0.000) 5.322 5.796 ( -1.344 -2.327 0.000) 2.687 5.885 ( 1.600 2.771 0.000) 3.200 6.167 ( 0.187 0.324 0.000) 0.375 6.262 ( -0.907 -1.572 0.000) 1.815 6.400 ( -6.793 -11.766 0.000) 13.587 6.555 ( -4.268 -7.393 0.000) 8.536 7.547 ( -2.668 -4.622 0.000) 5.337 7.627 ( -4.041 -6.999 0.000) 8.082 ======================= Grid point 48 (21/48) ======================= q-point: ( 0.29 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.537 ( 3.756 6.760 0.000) 7.733 1.667 ( -0.305 4.538 0.000) 4.548 1.787 ( 2.415 4.196 0.000) 4.841 1.916 ( 1.097 -0.292 0.000) 1.135 1.920 ( 3.501 4.570 0.000) 5.757 1.977 ( 0.183 0.830 0.000) 0.850 2.375 ( -2.770 -3.360 0.000) 4.355 2.481 ( -2.071 -5.072 0.000) 5.478 2.972 ( 2.305 18.625 0.000) 18.767 3.002 ( 2.735 19.543 0.000) 19.733 3.922 ( 6.030 3.934 0.000) 7.200 3.929 ( 7.533 5.578 0.000) 9.373 5.480 ( 2.106 1.767 0.000) 2.749 5.512 ( 2.281 -0.087 0.000) 2.283 5.625 ( 0.130 4.306 0.000) 4.308 5.631 ( -1.900 5.020 0.000) 5.368 5.757 ( 0.471 -3.628 0.000) 3.658 5.932 ( 0.062 1.027 0.000) 1.029 6.156 ( -3.815 3.067 0.000) 4.895 6.207 ( -4.294 -5.762 0.000) 7.186 6.249 ( 1.437 -1.795 0.000) 2.300 6.415 ( -6.597 -2.301 0.000) 6.987 7.477 ( -0.417 -4.128 0.000) 4.149 7.528 ( -0.099 -6.152 0.000) 6.153 ======================= Grid point 49 (22/48) ======================= q-point: ( 0.36 0.21 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.629 ( -0.594 5.699 0.000) 5.730 1.695 ( -1.520 4.523 0.000) 4.772 1.844 ( 0.379 2.244 0.000) 2.276 1.936 ( 1.853 -1.233 0.000) 2.226 1.984 ( -0.207 0.709 0.000) 0.739 1.993 ( 1.639 1.150 0.000) 2.002 2.312 ( -1.658 -0.914 0.000) 1.893 2.401 ( -0.168 -5.386 0.000) 5.389 3.137 ( -6.384 20.479 0.000) 21.451 3.177 ( -6.628 21.401 0.000) 22.404 4.008 ( 4.268 -3.596 0.000) 5.581 4.046 ( 4.698 -2.615 0.000) 5.376 5.547 ( 3.060 1.975 0.000) 3.642 5.579 ( 4.305 0.956 0.000) 4.409 5.612 ( -3.870 2.030 0.000) 4.371 5.638 ( -1.479 0.246 0.000) 1.499 5.744 ( 1.732 -3.215 0.000) 3.652 5.891 ( 0.109 -5.918 0.000) 5.919 6.118 ( -4.364 4.714 0.000) 6.424 6.168 ( -2.768 7.147 0.000) 7.664 6.260 ( 0.957 -0.866 0.000) 1.291 6.297 ( -5.265 0.186 0.000) 5.268 7.449 ( 0.657 -2.350 0.000) 2.440 7.493 ( 1.785 -4.545 0.000) 4.883 ======================= Grid point 62 (23/48) ======================= q-point: ( 0.29 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 114 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.671 ( 3.509 6.078 0.000) 7.019 1.730 ( 1.185 2.053 0.000) 2.371 1.848 ( 1.110 1.922 0.000) 2.220 1.913 ( 0.010 0.017 0.000) 0.020 1.982 ( 1.100 1.905 0.000) 2.199 1.990 ( 0.289 0.500 0.000) 0.578 2.323 ( -1.018 -1.763 0.000) 2.036 2.374 ( -2.935 -5.083 0.000) 5.870 3.321 ( 9.178 15.896 0.000) 18.355 3.369 ( 9.653 16.719 0.000) 19.305 3.928 ( -1.803 -3.123 0.000) 3.606 3.968 ( -0.960 -1.664 0.000) 1.921 5.525 ( 1.693 2.933 0.000) 3.387 5.525 ( 1.047 1.814 0.000) 2.095 5.665 ( 0.183 0.317 0.000) 0.366 5.676 ( -0.312 -0.540 0.000) 0.624 5.692 ( -1.290 -2.234 0.000) 2.580 5.858 ( -3.864 -6.693 0.000) 7.729 6.216 ( 2.605 4.511 0.000) 5.209 6.216 ( 0.811 1.405 0.000) 1.622 6.218 ( -0.111 -0.191 0.000) 0.221 6.372 ( -1.427 -2.471 0.000) 2.853 7.431 ( -0.500 -0.866 0.000) 1.000 7.449 ( -1.224 -2.121 0.000) 2.449 ======================= Grid point 63 (24/48) ======================= q-point: ( 0.36 0.29 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 112 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.732 ( -2.305 3.993 0.000) 4.610 1.757 ( -0.826 1.431 0.000) 1.653 1.864 ( 0.034 -0.059 0.000) 0.068 1.915 ( 0.403 -0.699 0.000) 0.807 1.995 ( -0.232 0.403 0.000) 0.465 2.000 ( 0.056 -0.098 0.000) 0.113 2.306 ( -0.017 0.030 0.000) 0.035 2.317 ( 1.501 -2.599 0.000) 3.001 3.513 ( -9.530 16.506 0.000) 19.060 3.569 ( -9.867 17.091 0.000) 19.735 3.862 ( 5.995 -10.384 0.000) 11.990 3.918 ( 5.652 -9.790 0.000) 11.305 5.556 ( 0.846 -1.465 0.000) 1.692 5.571 ( -0.554 0.959 0.000) 1.107 5.640 ( -0.700 1.212 0.000) 1.399 5.667 ( -1.450 2.511 0.000) 2.899 5.682 ( 1.700 -2.944 0.000) 3.400 5.770 ( 2.997 -5.191 0.000) 5.994 6.203 ( -1.850 3.204 0.000) 3.700 6.242 ( 0.594 -1.029 0.000) 1.188 6.303 ( -0.927 1.606 0.000) 1.855 6.312 ( -2.012 3.485 0.000) 4.024 7.426 ( 0.205 -0.356 0.000) 0.411 7.432 ( 0.960 -1.663 0.000) 1.921 ======================= Grid point 211 (25/48) ======================= q-point: ( 0.00 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 48 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.451 ( -0.000 -0.000 20.943) 20.943 0.451 ( 0.000 0.000 20.943) 20.943 0.840 ( -0.000 0.000 39.397) 39.397 0.864 ( 0.000 0.000 -18.074) 18.074 0.864 ( 0.000 0.000 -18.074) 18.074 1.555 ( 0.000 0.000 10.040) 10.040 1.555 ( 0.000 0.000 10.040) 10.040 1.634 ( 0.000 0.000 -35.720) 35.720 1.720 ( -0.000 -0.000 -5.470) 5.470 1.720 ( 0.000 0.000 -5.470) 5.470 2.993 ( 0.000 0.000 22.975) 22.975 3.378 ( 0.000 0.000 -13.023) 13.023 4.703 ( 0.000 -0.000 0.379) 0.379 4.703 ( 0.000 0.000 0.379) 0.379 4.712 ( -0.000 0.000 -0.391) 0.391 4.712 ( 0.000 -0.000 -0.391) 0.391 6.046 ( 0.000 -0.000 0.499) 0.499 6.046 ( 0.000 0.000 0.499) 0.499 6.058 ( 0.000 0.000 -0.490) 0.490 6.058 ( -0.000 0.000 -0.490) 0.490 7.565 ( 0.000 -0.000 2.569) 2.569 7.625 ( 0.000 0.000 -2.661) 2.661 8.046 ( -0.000 0.000 0.692) 0.692 8.060 ( -0.000 0.000 -0.524) 0.524 ======================= Grid point 212 (26/48) ======================= q-point: ( 0.07 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.638 ( 13.217 7.631 14.709) 21.196 0.675 ( 11.669 6.737 16.493) 21.297 0.976 ( 8.589 4.959 -16.317) 19.094 1.027 ( 9.436 5.448 -16.779) 20.006 1.093 ( 20.901 12.067 26.715) 36.002 1.582 ( 2.519 1.455 9.979) 10.394 1.633 ( 7.612 4.395 10.765) 13.898 1.723 ( 8.425 4.864 -30.506) 32.020 1.744 ( 2.101 1.213 -5.294) 5.823 1.807 ( 7.986 4.611 -5.797) 10.892 2.844 ( -7.039 -4.064 15.456) 17.463 3.022 ( -9.285 -5.361 -1.184) 10.787 4.749 ( 3.958 2.285 0.311) 4.581 4.757 ( 3.826 2.209 -0.321) 4.429 5.275 ( 36.411 21.022 12.488) 43.859 5.752 ( 37.117 21.429 -34.742) 55.171 6.004 ( -3.633 -2.097 0.458) 4.220 6.014 ( -3.673 -2.121 -0.432) 4.263 6.021 ( -2.183 -1.260 0.476) 2.565 6.032 ( -2.229 -1.287 -0.469) 2.616 7.523 ( -4.882 -2.819 3.733) 6.762 7.597 ( -2.691 -1.554 -3.050) 4.354 7.920 ( -9.713 -5.608 2.158) 11.421 7.987 ( -3.470 -2.004 -4.685) 6.165 ======================= Grid point 213 (27/48) ======================= q-point: ( 0.14 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.931 ( 11.724 6.769 8.970) 16.240 0.936 ( 12.013 6.936 9.524) 16.826 1.164 ( 3.787 2.186 -11.904) 12.681 1.180 ( 8.205 4.737 -12.709) 15.852 1.549 ( 15.576 8.993 14.252) 22.948 1.678 ( 6.116 3.531 8.614) 11.139 1.796 ( 5.858 3.382 -0.648) 6.795 1.818 ( 4.492 2.593 -4.601) 6.933 2.084 ( 20.927 12.082 -8.952) 25.769 2.122 ( 20.067 11.586 -5.830) 23.893 2.785 ( 2.066 1.193 9.425) 9.723 2.913 ( -2.227 -1.286 -2.954) 3.916 4.869 ( 6.118 3.532 0.152) 7.066 4.872 ( 5.939 3.429 -0.158) 6.859 5.884 ( -5.245 -3.028 1.297) 6.193 5.902 ( -4.663 -2.692 -0.318) 5.394 5.942 ( 17.439 10.068 9.193) 22.136 5.957 ( -3.185 -1.839 0.400) 3.699 5.966 ( -3.296 -1.903 -0.396) 3.826 6.200 ( 8.453 4.880 -12.925) 16.196 7.331 ( -11.344 -6.550 4.528) 13.860 7.441 ( -12.829 -7.407 -5.636) 15.849 7.748 ( -2.368 -1.367 8.759) 9.176 7.954 ( -0.038 -0.022 -10.099) 10.099 ======================= Grid point 214 (28/48) ======================= q-point: ( 0.21 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.161 ( 7.892 4.556 6.439) 11.158 1.203 ( 11.987 6.921 1.120) 13.887 1.261 ( 5.542 3.200 -3.851) 7.469 1.323 ( 4.600 2.656 -7.979) 9.585 1.794 ( 7.441 4.296 6.889) 11.013 1.851 ( 8.740 5.046 5.124) 11.318 1.899 ( 3.308 1.910 -1.394) 4.066 1.944 ( 6.501 3.754 -3.291) 8.197 2.493 ( 12.328 7.118 0.813) 14.258 2.590 ( 15.361 8.869 -7.487) 19.253 2.925 ( 7.858 4.537 4.101) 9.957 2.969 ( 13.751 7.939 -0.217) 15.880 4.999 ( 4.752 2.744 0.003) 5.487 4.999 ( 4.887 2.822 -0.005) 5.643 5.736 ( -6.572 -3.795 1.648) 7.766 5.771 ( -5.697 -3.289 -1.338) 6.713 5.894 ( -1.957 -1.130 0.248) 2.273 5.900 ( -2.156 -1.245 -0.246) 2.502 6.172 ( 3.195 1.845 4.765) 6.026 6.289 ( 0.913 0.527 -5.416) 5.517 7.044 ( -13.585 -7.843 2.020) 15.816 7.089 ( -16.333 -9.430 -1.904) 18.956 7.771 ( 2.485 1.435 7.642) 8.163 7.953 ( -0.570 -0.329 -8.102) 8.129 ======================= Grid point 215 (29/48) ======================= q-point: ( 0.29 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.303 ( 4.980 2.875 4.097) 7.061 1.402 ( 2.634 1.521 -4.497) 5.429 1.412 ( 7.826 4.518 1.503) 9.160 1.465 ( 10.673 6.162 -3.257) 12.747 1.918 ( 3.881 2.241 3.048) 5.420 1.951 ( 1.490 0.860 0.099) 1.723 2.050 ( 8.534 4.927 2.435) 10.151 2.100 ( 7.003 4.043 -1.917) 8.311 2.645 ( 2.581 1.490 0.235) 2.990 2.675 ( -2.615 -1.510 -2.396) 3.855 3.313 ( 20.931 12.084 2.872) 24.339 3.368 ( 18.799 10.854 -1.989) 21.798 5.074 ( 1.951 1.127 0.113) 2.256 5.077 ( 2.053 1.185 -0.114) 2.374 5.613 ( -3.750 -2.165 2.352) 4.927 5.667 ( -3.007 -1.736 -2.278) 4.153 5.874 ( -0.006 -0.003 0.023) 0.024 5.875 ( -0.250 -0.145 -0.022) 0.290 6.177 ( -1.428 -0.825 3.832) 4.172 6.276 ( -1.506 -0.869 -4.857) 5.159 6.725 ( -13.907 -8.029 1.897) 16.170 6.756 ( -13.044 -7.531 -0.966) 15.093 7.803 ( -0.276 -0.159 4.779) 4.790 7.915 ( -2.935 -1.695 -4.867) 5.931 ======================= Grid point 216 (30/48) ======================= q-point: ( 0.36 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.394 ( 3.137 1.811 2.389) 4.339 1.450 ( 1.734 1.001 -2.470) 3.180 1.594 ( 7.863 4.540 3.141) 9.608 1.683 ( 8.497 4.906 -4.822) 10.933 1.973 ( 0.630 0.364 1.094) 1.313 1.982 ( 1.936 1.118 0.664) 2.332 2.225 ( 6.750 3.897 1.096) 7.871 2.249 ( 5.954 3.437 -0.973) 6.943 2.577 ( -5.466 -3.156 2.813) 6.910 2.642 ( -2.488 -1.437 -2.798) 4.011 3.745 ( 16.489 9.520 0.485) 19.046 3.755 ( 14.773 8.529 -0.368) 17.062 5.103 ( 0.937 0.541 0.224) 1.105 5.108 ( 1.067 0.616 -0.227) 1.253 5.586 ( 1.566 0.904 2.897) 3.415 5.654 ( 2.187 1.263 -3.041) 3.953 5.873 ( -0.212 -0.122 0.199) 0.315 5.878 ( -0.019 -0.011 -0.196) 0.197 6.134 ( -2.301 -1.328 3.978) 4.783 6.236 ( -2.129 -1.229 -5.150) 5.706 6.452 ( -9.893 -5.711 2.433) 11.679 6.493 ( -9.975 -5.759 -1.377) 11.600 7.753 ( -3.908 -2.256 3.046) 5.444 7.821 ( -5.121 -2.957 -2.872) 6.574 ======================= Grid point 217 (31/48) ======================= q-point: ( 0.43 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.446 ( 1.431 0.826 1.581) 2.287 1.482 ( 1.083 0.625 -1.594) 2.026 1.745 ( 5.033 2.906 3.815) 6.952 1.847 ( 5.815 3.357 -5.172) 8.476 1.984 ( 0.478 0.276 1.265) 1.380 2.015 ( 0.925 0.534 -1.257) 1.649 2.347 ( 3.822 2.207 0.529) 4.445 2.359 ( 3.515 2.029 -0.494) 4.089 2.450 ( -5.079 -2.932 5.496) 8.037 2.557 ( -4.267 -2.464 -3.924) 6.299 4.014 ( 7.840 4.526 0.853) 9.093 4.033 ( 8.671 5.006 -0.810) 10.045 5.124 ( 0.876 0.506 0.350) 1.070 5.132 ( 0.992 0.573 -0.358) 1.200 5.676 ( 5.704 3.293 3.875) 7.642 5.763 ( 6.939 4.007 -3.686) 8.820 5.861 ( -0.720 -0.416 0.345) 0.900 5.869 ( -0.623 -0.360 -0.337) 0.795 6.064 ( -4.119 -2.378 3.359) 5.823 6.147 ( -6.351 -3.667 -3.984) 8.346 6.303 ( -2.933 -1.694 0.890) 3.502 6.324 ( -4.643 -2.681 -0.927) 5.441 7.653 ( -4.151 -2.396 2.454) 5.385 7.705 ( -4.511 -2.604 -2.057) 5.600 ======================= Grid point 218 (32/48) ======================= q-point: (-0.50 0.00 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 89 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.461 ( 0.000 0.000 1.481) 1.481 1.496 ( 0.000 0.000 -1.482) 1.482 1.804 ( 0.000 0.000 4.688) 4.688 1.918 ( 0.000 0.000 -4.305) 4.305 2.003 ( 0.000 0.000 -1.453) 1.453 2.016 ( 0.000 0.000 -1.189) 1.189 2.385 ( 0.000 0.000 6.424) 6.424 2.391 ( 0.000 0.000 0.382) 0.382 2.400 ( 0.000 0.000 -0.359) 0.359 2.497 ( 0.000 0.000 -3.853) 3.853 4.103 ( 0.000 0.000 1.237) 1.237 4.132 ( 0.000 0.000 -1.211) 1.211 5.135 ( 0.000 0.000 0.412) 0.412 5.145 ( 0.000 0.000 -0.423) 0.423 5.772 ( 0.000 0.000 5.649) 5.649 5.851 ( 0.000 0.000 0.393) 0.393 5.860 ( 0.000 0.000 -0.381) 0.381 5.903 ( 0.000 0.000 -5.706) 5.706 5.988 ( 0.000 0.000 1.104) 1.104 6.013 ( -0.000 -0.000 -0.925) 0.925 6.273 ( -0.000 -0.000 -0.092) 0.092 6.281 ( 0.000 0.000 -0.594) 0.594 7.600 ( 0.000 0.000 2.350) 2.350 7.649 ( 0.000 0.000 -1.856) 1.856 ======================= Grid point 227 (33/48) ======================= q-point: ( 0.07 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.852 ( 6.273 10.866 10.964) 16.662 0.885 ( 8.447 14.631 10.044) 19.655 1.130 ( 2.283 3.955 -13.681) 14.423 1.147 ( 6.013 10.415 -13.179) 17.841 1.430 ( 10.025 17.363 17.054) 26.321 1.664 ( 4.058 7.028 9.291) 12.336 1.743 ( 3.036 5.259 0.263) 6.079 1.808 ( 3.274 5.670 -3.146) 7.264 1.958 ( 10.911 18.898 -14.546) 26.225 1.997 ( 9.223 15.975 -5.783) 19.331 2.778 ( -0.444 -0.768 11.294) 11.329 2.929 ( -1.925 -3.334 -3.311) 5.077 4.841 ( 3.885 6.729 0.154) 7.771 4.844 ( 3.729 6.459 -0.162) 7.460 5.800 ( 14.521 25.151 11.813) 31.352 5.922 ( -3.348 -5.799 0.680) 6.730 5.937 ( -3.116 -5.397 -0.458) 6.248 5.977 ( -1.755 -3.039 0.389) 3.531 5.986 ( -1.839 -3.186 -0.381) 3.698 6.135 ( 7.304 12.650 -17.637) 22.900 7.398 ( -5.729 -9.923 4.996) 12.500 7.513 ( -5.090 -8.816 -5.644) 11.639 7.771 ( -3.676 -6.368 7.125) 10.239 7.947 ( -0.794 -1.375 -9.311) 9.445 ======================= Grid point 228 (34/48) ======================= q-point: ( 0.14 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.084 ( 8.621 9.341 2.582) 12.971 1.134 ( 6.115 8.770 3.551) 11.266 1.214 ( 3.847 5.761 -3.259) 7.656 1.326 ( 3.785 8.160 -7.704) 11.843 1.693 ( 8.226 6.625 8.065) 13.289 1.825 ( 3.797 7.673 2.986) 9.067 1.853 ( 5.318 5.664 0.485) 7.784 1.937 ( 3.583 8.210 -2.410) 9.276 2.361 ( 13.514 14.093 -2.125) 19.641 2.417 ( 17.226 17.377 -5.445) 25.067 2.854 ( 6.932 6.743 3.788) 10.386 2.895 ( 0.958 0.226 -0.126) 0.992 4.989 ( 4.082 8.307 0.070) 9.256 4.991 ( 4.132 8.676 -0.069) 9.610 5.787 ( -5.101 -6.160 1.344) 8.110 5.814 ( -4.502 -5.590 -0.970) 7.242 5.913 ( -2.158 -2.482 0.228) 3.297 5.918 ( -2.217 -2.770 -0.225) 3.555 6.133 ( 5.891 7.120 5.610) 10.811 6.270 ( 2.126 1.748 -6.238) 6.819 7.157 ( -9.967 -10.807 3.342) 15.077 7.237 ( -12.744 -12.674 -3.785) 18.368 7.733 ( 1.616 -0.329 8.114) 8.280 7.927 ( 0.420 -2.832 -8.854) 9.306 ======================= Grid point 229 (35/48) ======================= q-point: ( 0.21 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.263 ( 5.295 4.365 1.826) 7.101 1.313 ( 6.497 4.883 -2.515) 8.507 1.379 ( 5.907 9.908 1.735) 11.665 1.450 ( 2.377 8.072 -3.344) 9.055 1.854 ( 5.490 2.645 3.854) 7.210 1.915 ( 2.890 1.041 -1.059) 3.249 2.038 ( 4.792 12.135 3.046) 13.398 2.093 ( 3.768 10.707 -1.846) 11.500 2.600 ( 5.100 3.604 1.858) 6.515 2.686 ( 0.999 0.920 -5.777) 5.934 3.101 ( 17.576 13.203 4.472) 22.433 3.174 ( 16.769 12.575 -2.238) 21.079 5.118 ( 0.469 9.236 0.187) 9.250 5.123 ( 0.371 9.492 -0.185) 9.501 5.648 ( -4.290 -4.805 1.984) 6.740 5.691 ( -3.554 -4.611 -1.695) 6.063 5.870 ( -0.445 -1.339 0.090) 1.414 5.872 ( -0.471 -1.412 -0.112) 1.492 6.194 ( -0.844 0.325 3.516) 3.630 6.283 ( -0.340 -1.059 -4.255) 4.398 6.866 ( -12.056 -9.970 1.495) 15.715 6.893 ( -13.422 -10.244 -0.922) 16.909 7.763 ( 2.425 -2.526 5.516) 6.534 7.894 ( 0.299 -5.578 -5.780) 8.038 ======================= Grid point 230 (36/48) ======================= q-point: ( 0.29 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.385 ( 3.751 4.285 2.062) 6.057 1.423 ( 2.533 3.497 -0.795) 4.391 1.554 ( 4.251 7.736 0.736) 8.858 1.608 ( 6.850 7.700 -3.617) 10.922 1.941 ( 3.297 0.500 1.697) 3.742 1.956 ( 1.525 0.104 0.216) 1.544 2.231 ( 3.274 12.514 1.517) 13.024 2.261 ( 2.653 11.883 -1.130) 12.227 2.632 ( -4.214 -2.464 0.428) 4.901 2.651 ( -0.077 -0.782 -1.164) 1.405 3.552 ( 18.250 12.030 1.453) 21.906 3.583 ( 16.347 10.880 -1.237) 19.676 5.185 ( -3.203 9.789 0.170) 10.301 5.189 ( -3.249 9.747 -0.167) 10.275 5.570 ( -0.459 -2.046 2.352) 3.151 5.624 ( 0.339 -2.536 -2.322) 3.455 5.866 ( 0.722 -0.887 0.127) 1.151 5.869 ( 0.572 -0.275 -0.175) 0.659 6.160 ( -2.157 -0.894 3.617) 4.305 6.254 ( -1.390 -1.342 -4.654) 5.039 6.572 ( -11.533 -7.462 1.763) 13.849 6.600 ( -10.720 -8.135 -0.784) 13.479 7.746 ( 0.016 -5.484 3.385) 6.444 7.824 ( -1.192 -7.205 -3.364) 8.041 ======================= Grid point 231 (37/48) ======================= q-point: ( 0.36 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.465 ( 0.459 4.816 0.880) 4.917 1.482 ( 0.987 2.741 -0.480) 2.952 1.700 ( 4.757 5.668 2.681) 7.870 1.784 ( 6.120 5.153 -4.813) 9.336 1.974 ( 0.718 -0.062 0.961) 1.201 1.986 ( 2.161 -0.286 0.130) 2.184 2.378 ( 0.649 9.111 1.311) 9.227 2.402 ( 0.294 10.918 -0.900) 10.959 2.526 ( -5.727 0.263 3.395) 6.663 2.601 ( -3.697 -1.755 -3.035) 5.095 3.908 ( 11.135 6.885 0.236) 13.093 3.913 ( 12.456 7.391 -0.200) 14.485 5.217 ( -4.428 10.378 0.104) 11.284 5.219 ( -4.256 9.962 -0.119) 10.834 5.601 ( 4.804 0.609 2.852) 5.620 5.664 ( 5.976 -0.307 -2.589) 6.520 5.863 ( 0.195 -1.380 0.260) 1.418 5.872 ( -0.266 0.144 -0.525) 0.606 6.105 ( -3.239 -1.607 3.702) 5.174 6.202 ( -3.633 -2.531 -5.042) 6.710 6.354 ( -6.503 -4.211 1.878) 7.972 6.381 ( -6.777 -5.212 -0.827) 8.590 7.667 ( -2.050 -5.955 2.475) 6.767 7.721 ( -2.404 -6.655 -2.224) 7.417 ======================= Grid point 232 (38/48) ======================= q-point: ( 0.43 0.07 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.499 ( -1.545 4.195 0.674) 4.521 1.513 ( -0.427 2.467 -0.552) 2.564 1.801 ( 1.754 2.645 3.833) 4.976 1.904 ( 2.912 2.164 -4.960) 6.145 1.985 ( 1.000 -0.043 0.557) 1.145 2.009 ( 0.764 -0.686 -1.189) 1.572 2.382 ( -1.683 -1.934 5.006) 5.624 2.466 ( -1.334 3.333 -2.536) 4.395 2.511 ( -4.559 11.273 0.728) 12.182 2.533 ( -5.376 6.182 -1.347) 8.302 4.081 ( 3.905 2.155 0.996) 4.570 4.104 ( 4.446 2.085 -0.982) 5.007 5.238 ( -5.169 10.292 0.040) 11.517 5.240 ( -4.792 9.592 -0.107) 10.723 5.713 ( 6.003 0.076 4.275) 7.370 5.785 ( 6.581 -1.664 -1.583) 6.971 5.851 ( 0.250 -1.450 0.339) 1.510 5.881 ( 1.714 2.085 -2.652) 3.784 6.017 ( -4.036 -1.787 1.813) 4.772 6.064 ( -6.648 -4.322 -2.331) 8.265 6.277 ( -1.466 -1.463 0.112) 2.074 6.283 ( -0.903 -0.992 -0.348) 1.386 7.589 ( -0.688 -3.752 2.262) 4.435 7.636 ( -0.682 -4.012 -1.865) 4.477 ======================= Grid point 242 (39/48) ======================= q-point: ( 0.14 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.232 ( 2.663 4.612 0.723) 5.374 1.278 ( 5.626 9.745 -3.811) 11.880 1.400 ( 5.519 9.559 3.319) 11.526 1.470 ( 3.948 6.838 -2.019) 8.150 1.816 ( 3.128 5.418 3.810) 7.325 1.893 ( 1.451 2.513 -2.186) 3.633 2.067 ( 7.192 12.457 3.474) 14.797 2.125 ( 6.206 10.749 -1.779) 12.539 2.569 ( 3.826 6.628 2.455) 8.037 2.675 ( 2.760 4.780 -7.163) 9.043 3.012 ( 8.765 15.182 4.974) 18.222 3.085 ( 9.455 16.376 -1.993) 19.014 5.172 ( 5.466 9.467 0.315) 10.936 5.179 ( 5.541 9.597 -0.316) 11.086 5.661 ( -3.585 -6.209 1.783) 7.388 5.698 ( -3.385 -5.863 -1.391) 6.912 5.867 ( -1.115 -1.931 0.064) 2.230 5.868 ( -1.124 -1.946 -0.030) 2.247 6.206 ( 0.429 0.742 2.744) 2.875 6.276 ( -0.394 -0.682 -3.305) 3.398 6.921 ( -7.477 -12.950 2.254) 15.123 6.969 ( -8.188 -14.183 -1.918) 16.489 7.720 ( -0.748 -1.296 5.444) 5.646 7.851 ( -2.713 -4.699 -5.833) 7.966 ======================= Grid point 243 (40/48) ======================= q-point: ( 0.21 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.341 ( 4.627 4.668 3.640) 7.513 1.439 ( 3.317 6.722 -5.011) 9.016 1.563 ( 3.944 8.342 1.765) 9.395 1.607 ( 3.797 7.849 -1.404) 8.832 1.899 ( 2.581 1.854 2.135) 3.828 1.934 ( 1.395 1.015 -0.950) 1.969 2.303 ( 4.816 13.600 1.837) 14.544 2.339 ( 4.358 13.388 -1.297) 14.139 2.642 ( 1.089 1.437 0.732) 1.946 2.671 ( -2.127 -1.032 -1.769) 2.953 3.396 ( 14.923 15.963 2.260) 21.968 3.443 ( 13.467 14.467 -1.897) 19.855 5.328 ( 0.612 10.873 0.218) 10.892 5.333 ( 0.390 10.487 -0.203) 10.496 5.564 ( -1.867 -3.157 1.689) 4.039 5.602 ( -1.688 -4.073 -1.536) 4.669 5.844 ( 0.205 -1.266 0.456) 1.361 5.853 ( 0.654 -0.311 -0.320) 0.792 6.189 ( -1.764 -0.648 2.406) 3.053 6.254 ( -0.592 -1.737 -3.376) 3.842 6.662 ( -9.017 -10.762 2.013) 14.184 6.698 ( -9.803 -9.414 -1.211) 13.644 7.679 ( -0.084 -5.627 3.178) 6.463 7.754 ( -1.110 -8.043 -3.308) 8.767 ======================= Grid point 244 (41/48) ======================= q-point: ( 0.29 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 296 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.463 ( 3.851 4.629 3.070) 6.759 1.529 ( 0.390 5.924 -2.424) 6.412 1.697 ( 3.242 6.544 1.360) 7.429 1.751 ( 4.469 6.549 -3.515) 8.673 1.948 ( 2.152 0.323 1.338) 2.555 1.961 ( 0.943 0.482 -0.098) 1.064 2.443 ( -0.399 3.474 2.255) 4.161 2.490 ( 0.417 5.635 -1.824) 5.938 2.657 ( -2.508 10.310 0.726) 10.636 2.675 ( -1.207 8.651 -0.811) 8.772 3.772 ( 13.434 10.097 0.519) 16.814 3.784 ( 12.100 9.265 -0.514) 15.249 5.403 ( -3.054 9.777 0.256) 10.246 5.408 ( -3.223 10.914 -0.157) 11.382 5.547 ( 1.225 0.716 1.366) 1.970 5.578 ( 2.357 -1.409 -1.296) 3.036 5.838 ( 0.593 -1.910 1.333) 2.403 5.871 ( 0.712 0.478 -1.452) 1.686 6.147 ( -2.696 -0.160 2.710) 3.826 6.222 ( -1.386 -1.983 -4.084) 4.746 6.424 ( -7.781 -7.495 2.500) 11.090 6.464 ( -8.958 -5.381 -1.251) 10.524 7.613 ( -0.008 -7.290 2.118) 7.591 7.661 ( -0.125 -8.608 -2.045) 8.848 ======================= Grid point 245 (42/48) ======================= q-point: ( 0.36 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.548 ( 0.652 4.494 1.409) 4.755 1.576 ( -1.349 5.592 -1.082) 5.853 1.801 ( 2.335 4.314 2.283) 5.410 1.871 ( 3.677 3.375 -4.160) 6.497 1.977 ( 0.433 0.389 0.552) 0.802 1.983 ( 1.737 0.118 0.216) 1.755 2.395 ( -2.652 -3.183 3.645) 5.518 2.470 ( -1.465 -3.806 -2.945) 5.031 2.755 ( -5.169 18.509 0.729) 19.231 2.769 ( -5.193 18.288 -0.577) 19.020 4.012 ( 7.216 3.321 0.505) 7.960 4.024 ( 8.003 3.593 -0.513) 8.787 5.437 ( -4.437 10.239 0.370) 11.165 5.447 ( -4.622 10.871 -0.498) 11.823 5.611 ( 4.865 0.505 2.190) 5.359 5.648 ( 6.476 -1.006 -1.122) 6.649 5.824 ( 0.486 -2.601 1.706) 3.148 5.878 ( 0.152 0.031 -3.021) 3.025 6.085 ( -4.787 0.391 2.203) 5.284 6.137 ( -5.998 -3.048 -2.469) 7.167 6.291 ( -1.958 -2.072 0.943) 3.003 6.307 ( -4.595 -2.256 -0.499) 5.144 7.547 ( 0.105 -5.694 1.846) 5.987 7.587 ( 0.387 -6.465 -1.656) 6.685 ======================= Grid point 246 (43/48) ======================= q-point: ( 0.43 0.14 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.574 ( -2.103 3.642 0.877) 4.296 1.593 ( -2.884 4.995 -0.716) 5.812 1.844 ( -0.920 1.594 2.891) 3.427 1.924 ( 0.100 -0.174 -4.149) 4.154 1.984 ( -0.036 0.062 0.715) 0.718 2.001 ( 0.139 -0.241 -0.722) 0.774 2.349 ( 0.654 -1.133 4.073) 4.278 2.429 ( 1.638 -2.836 -3.069) 4.488 2.805 ( -10.149 17.578 0.695) 20.310 2.818 ( -10.297 17.835 -0.512) 20.600 4.090 ( 0.896 -1.552 0.851) 1.983 4.109 ( 0.984 -1.704 -0.853) 2.145 5.450 ( -5.859 10.148 0.363) 11.723 5.462 ( -6.015 10.418 -0.697) 12.050 5.671 ( 2.261 -3.915 3.457) 5.692 5.717 ( 3.022 -5.235 -0.712) 6.086 5.814 ( 1.289 -2.233 1.497) 2.982 5.888 ( 1.173 -2.032 -5.099) 5.612 6.032 ( -2.486 4.305 0.481) 4.995 6.037 ( -1.657 2.870 0.385) 3.336 6.261 ( -0.059 0.101 0.522) 0.535 6.272 ( 0.198 -0.343 -0.430) 0.585 7.518 ( 1.946 -3.371 1.823) 4.298 7.557 ( 2.313 -4.007 -1.596) 4.894 ======================= Grid point 258 (44/48) ======================= q-point: ( 0.21 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.451 ( 3.622 6.274 4.420) 8.486 1.556 ( 2.930 5.075 -4.741) 7.538 1.703 ( 3.507 6.075 1.719) 7.222 1.750 ( 3.740 6.477 -2.534) 7.897 1.929 ( 0.807 1.398 1.680) 2.329 1.955 ( 0.602 1.043 -0.727) 1.407 2.462 ( 0.202 0.350 2.134) 2.172 2.508 ( 0.595 1.030 -1.838) 2.190 2.744 ( 6.001 10.395 0.076) 12.003 2.746 ( 6.578 11.394 -0.098) 13.157 3.697 ( 8.199 14.201 0.866) 16.421 3.717 ( 7.434 12.875 -0.867) 14.893 5.457 ( 0.949 1.644 0.954) 2.125 5.478 ( 1.516 2.625 -0.835) 3.144 5.569 ( 1.311 2.271 0.012) 2.622 5.571 ( 2.295 3.975 -0.112) 4.591 5.820 ( -0.758 -1.313 2.032) 2.535 5.865 ( 0.724 1.253 -1.846) 2.346 6.178 ( -0.197 -0.342 1.110) 1.178 6.218 ( -1.078 -1.867 -2.667) 3.429 6.457 ( -5.590 -9.682 3.875) 11.833 6.528 ( -4.652 -8.057 -2.502) 9.634 7.567 ( -2.984 -5.169 1.725) 6.213 7.607 ( -3.672 -6.361 -1.723) 7.544 ======================= Grid point 259 (45/48) ======================= q-point: ( 0.29 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.571 ( 2.824 6.081 2.901) 7.305 1.636 ( 0.777 4.884 -2.730) 5.649 1.804 ( 2.381 4.295 1.592) 5.163 1.853 ( 2.522 3.688 -3.016) 5.390 1.959 ( 1.456 1.022 1.206) 2.149 1.974 ( 0.572 0.768 -0.317) 1.009 2.405 ( -2.568 -3.866 2.473) 5.259 2.458 ( -2.234 -4.705 -2.096) 5.614 2.980 ( 2.431 18.826 0.661) 18.994 2.995 ( 2.644 19.287 -0.619) 19.477 3.923 ( 6.411 4.360 0.156) 7.755 3.927 ( 7.163 5.182 -0.167) 8.842 5.488 ( 2.192 1.301 0.727) 2.651 5.504 ( 2.262 0.382 -0.671) 2.390 5.627 ( -0.275 4.300 0.216) 4.314 5.631 ( -0.975 4.764 -0.117) 4.865 5.790 ( -0.130 -2.899 3.154) 4.286 5.875 ( -0.413 -0.927 -4.330) 4.447 6.167 ( -2.216 2.283 0.841) 3.291 6.187 ( -1.568 -0.733 -1.039) 2.018 6.308 ( -3.246 -3.928 3.997) 6.476 6.384 ( -5.787 -2.900 -2.826) 7.063 7.491 ( -0.323 -4.651 1.134) 4.799 7.516 ( -0.164 -5.663 -1.078) 5.767 ======================= Grid point 260 (46/48) ======================= q-point: ( 0.36 0.21 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 294 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.646 ( -0.816 5.312 1.446) 5.565 1.679 ( -1.291 4.714 -1.404) 5.085 1.862 ( 0.722 1.906 1.704) 2.657 1.907 ( 1.488 0.419 -2.254) 2.733 1.987 ( 0.181 0.603 0.297) 0.696 1.993 ( 0.996 0.651 -0.127) 1.197 2.337 ( -1.146 -2.310 2.065) 3.303 2.381 ( -0.443 -4.475 -1.761) 4.830 3.147 ( -6.422 20.669 0.885) 21.662 3.167 ( -6.545 21.132 -0.843) 22.138 4.017 ( 4.373 -3.331 0.830) 5.559 4.036 ( 4.588 -2.841 -0.815) 5.457 5.553 ( 2.882 1.779 0.554) 3.431 5.567 ( 2.923 1.536 -0.749) 3.386 5.628 ( -1.608 0.522 0.904) 1.917 5.639 ( -1.304 0.231 -0.130) 1.330 5.768 ( 0.972 -2.782 2.482) 3.853 5.841 ( 0.262 -3.844 -3.836) 5.437 6.137 ( -4.435 3.959 1.367) 6.100 6.161 ( -3.629 5.632 -0.760) 6.743 6.270 ( -0.260 -0.069 0.881) 0.921 6.289 ( -3.434 0.391 -0.756) 3.537 7.461 ( 0.963 -2.953 0.974) 3.255 7.483 ( 1.526 -4.048 -0.928) 4.424 ======================= Grid point 273 (47/48) ======================= q-point: ( 0.29 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 162 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.685 ( 2.914 5.048 1.275) 5.967 1.715 ( 1.748 3.028 -1.278) 3.723 1.863 ( 0.995 1.723 1.293) 2.373 1.895 ( 0.516 0.894 -1.498) 1.819 1.985 ( 0.747 1.294 0.209) 1.509 1.989 ( 0.399 0.691 -0.117) 0.806 2.337 ( -1.567 -2.715 1.161) 3.343 2.362 ( -2.516 -4.357 -1.053) 5.140 3.333 ( 9.287 16.085 1.046) 18.603 3.357 ( 9.523 16.495 -1.054) 19.076 3.938 ( -1.580 -2.736 0.872) 3.277 3.958 ( -1.160 -2.008 -0.836) 2.465 5.524 ( 1.167 2.022 -0.039) 2.335 5.524 ( 1.486 2.575 -0.035) 2.973 5.657 ( -0.519 -0.899 0.081) 1.042 5.668 ( -0.480 -0.831 -0.681) 1.176 5.740 ( -0.954 -1.652 3.485) 3.973 5.817 ( -2.357 -4.082 -3.401) 5.812 6.201 ( 0.593 1.026 0.132) 1.192 6.210 ( 0.634 1.098 -0.521) 1.371 6.277 ( 0.604 1.046 3.486) 3.689 6.344 ( -0.824 -1.428 -2.522) 3.013 7.436 ( -0.690 -1.195 0.400) 1.437 7.445 ( -1.052 -1.822 -0.383) 2.138 ======================= Grid point 274 (48/48) ======================= q-point: ( 0.36 0.29 0.33) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.65e-04 1.65e-04 1.65e-04 1.65e-04 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 161 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.738 ( -1.928 3.339 0.505) 3.889 1.750 ( -1.190 2.060 -0.570) 2.446 1.876 ( 0.110 -0.191 1.108) 1.129 1.902 ( 0.293 -0.507 -1.124) 1.268 1.996 ( -0.160 0.277 0.143) 0.351 1.999 ( -0.017 0.030 -0.097) 0.103 2.309 ( 0.373 -0.646 0.230) 0.781 2.314 ( 1.130 -1.957 -0.217) 2.270 3.527 ( -9.589 16.608 1.203) 19.215 3.555 ( -9.749 16.887 -1.243) 19.538 3.877 ( 5.880 -10.184 1.246) 11.825 3.905 ( 5.716 -9.901 -1.175) 11.493 5.560 ( 0.540 -0.935 0.336) 1.131 5.567 ( -0.209 0.362 -0.300) 0.514 5.632 ( -0.160 0.278 -0.374) 0.493 5.638 ( 0.193 -0.334 -1.027) 1.097 5.729 ( 0.418 -0.724 2.488) 2.624 5.766 ( 1.751 -3.034 -0.785) 3.590 6.203 ( -1.151 1.993 0.308) 2.322 6.223 ( -0.058 0.101 -1.382) 1.387 6.312 ( -1.135 1.966 0.586) 2.344 6.317 ( -1.542 2.671 0.152) 3.088 7.427 ( 0.392 -0.679 0.116) 0.792 7.430 ( 0.769 -1.333 -0.118) 1.543 =================== 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 507.424 507.424 336.712 -0.000 -0.000 -0.000 3/14112 20.0 103.227 103.227 66.407 -0.000 -0.000 -0.000 3/14112 30.0 53.841 53.841 31.518 -0.000 -0.000 -0.000 3/14112 40.0 35.414 35.414 19.614 -0.000 -0.000 -0.000 3/14112 50.0 25.485 25.485 13.778 -0.000 -0.000 -0.000 3/14112 60.0 19.483 19.483 10.424 -0.000 -0.000 -0.000 3/14112 70.0 15.607 15.607 8.307 -0.000 -0.000 -0.000 3/14112 80.0 12.962 12.962 6.877 -0.000 -0.000 -0.000 3/14112 90.0 11.070 11.070 5.860 -0.000 -0.000 -0.000 3/14112 100.0 9.660 9.660 5.103 -0.000 -0.000 -0.000 3/14112 110.0 8.573 8.573 4.522 -0.000 -0.000 -0.000 3/14112 120.0 7.710 7.710 4.061 -0.000 -0.000 -0.000 3/14112 130.0 7.009 7.009 3.687 -0.000 -0.000 -0.000 3/14112 140.0 6.429 6.429 3.378 -0.000 -0.000 -0.000 3/14112 150.0 5.940 5.940 3.119 -0.000 -0.000 -0.000 3/14112 160.0 5.523 5.523 2.897 -0.000 -0.000 -0.000 3/14112 170.0 5.162 5.162 2.706 -0.000 -0.000 -0.000 3/14112 180.0 4.846 4.846 2.539 -0.000 -0.000 -0.000 3/14112 190.0 4.568 4.568 2.392 -0.000 -0.000 -0.000 3/14112 200.0 4.321 4.321 2.262 -0.000 -0.000 -0.000 3/14112 210.0 4.100 4.100 2.145 -0.000 -0.000 -0.000 3/14112 220.0 3.901 3.901 2.040 -0.000 -0.000 -0.000 3/14112 230.0 3.721 3.721 1.945 -0.000 -0.000 -0.000 3/14112 240.0 3.558 3.558 1.859 -0.000 -0.000 -0.000 3/14112 250.0 3.408 3.408 1.781 -0.000 -0.000 -0.000 3/14112 260.0 3.271 3.271 1.708 -0.000 -0.000 -0.000 3/14112 270.0 3.144 3.144 1.642 -0.000 -0.000 -0.000 3/14112 280.0 3.027 3.027 1.581 -0.000 -0.000 -0.000 3/14112 290.0 2.919 2.919 1.524 -0.000 -0.000 -0.000 3/14112 300.0 2.818 2.818 1.471 -0.000 -0.000 -0.000 3/14112 310.0 2.724 2.724 1.422 -0.000 -0.000 -0.000 3/14112 320.0 2.636 2.636 1.376 -0.000 -0.000 -0.000 3/14112 330.0 2.554 2.554 1.332 -0.000 -0.000 -0.000 3/14112 340.0 2.477 2.477 1.292 -0.000 -0.000 -0.000 3/14112 350.0 2.404 2.404 1.254 -0.000 -0.000 -0.000 3/14112 360.0 2.336 2.336 1.218 -0.000 -0.000 -0.000 3/14112 370.0 2.271 2.271 1.184 -0.000 -0.000 -0.000 3/14112 380.0 2.210 2.210 1.153 -0.000 -0.000 -0.000 3/14112 390.0 2.152 2.152 1.122 -0.000 -0.000 -0.000 3/14112 400.0 2.097 2.097 1.094 -0.000 -0.000 -0.000 3/14112 410.0 2.045 2.045 1.066 -0.000 -0.000 -0.000 3/14112 420.0 1.996 1.996 1.040 -0.000 -0.000 -0.000 3/14112 430.0 1.949 1.949 1.016 -0.000 -0.000 -0.000 3/14112 440.0 1.904 1.904 0.992 -0.000 -0.000 -0.000 3/14112 450.0 1.861 1.861 0.970 -0.000 -0.000 -0.000 3/14112 460.0 1.820 1.820 0.948 -0.000 -0.000 -0.000 3/14112 470.0 1.780 1.780 0.928 -0.000 -0.000 -0.000 3/14112 480.0 1.743 1.743 0.908 -0.000 -0.000 -0.000 3/14112 490.0 1.707 1.707 0.889 -0.000 -0.000 -0.000 3/14112 500.0 1.672 1.672 0.871 -0.000 -0.000 -0.000 3/14112 510.0 1.639 1.639 0.854 -0.000 -0.000 -0.000 3/14112 520.0 1.607 1.607 0.837 -0.000 -0.000 -0.000 3/14112 530.0 1.576 1.576 0.821 -0.000 -0.000 -0.000 3/14112 540.0 1.547 1.547 0.806 -0.000 -0.000 -0.000 3/14112 550.0 1.519 1.519 0.791 -0.000 -0.000 -0.000 3/14112 560.0 1.491 1.491 0.777 -0.000 -0.000 -0.000 3/14112 570.0 1.465 1.465 0.763 -0.000 -0.000 -0.000 3/14112 580.0 1.439 1.439 0.750 -0.000 -0.000 -0.000 3/14112 590.0 1.415 1.415 0.737 -0.000 -0.000 -0.000 3/14112 600.0 1.391 1.391 0.725 -0.000 -0.000 -0.000 3/14112 610.0 1.368 1.368 0.713 -0.000 -0.000 -0.000 3/14112 620.0 1.346 1.346 0.701 -0.000 -0.000 -0.000 3/14112 630.0 1.324 1.324 0.690 -0.000 -0.000 -0.000 3/14112 640.0 1.303 1.303 0.679 -0.000 -0.000 -0.000 3/14112 650.0 1.283 1.283 0.668 -0.000 -0.000 -0.000 3/14112 660.0 1.264 1.264 0.658 -0.000 -0.000 -0.000 3/14112 670.0 1.245 1.245 0.648 -0.000 -0.000 -0.000 3/14112 680.0 1.226 1.226 0.639 -0.000 -0.000 -0.000 3/14112 690.0 1.208 1.208 0.629 -0.000 -0.000 -0.000 3/14112 700.0 1.191 1.191 0.620 -0.000 -0.000 -0.000 3/14112 710.0 1.174 1.174 0.612 -0.000 -0.000 -0.000 3/14112 720.0 1.158 1.158 0.603 -0.000 -0.000 -0.000 3/14112 730.0 1.142 1.142 0.595 -0.000 -0.000 -0.000 3/14112 740.0 1.126 1.126 0.587 -0.000 -0.000 -0.000 3/14112 750.0 1.111 1.111 0.579 -0.000 -0.000 -0.000 3/14112 760.0 1.096 1.096 0.571 -0.000 -0.000 -0.000 3/14112 770.0 1.082 1.082 0.564 -0.000 -0.000 -0.000 3/14112 780.0 1.068 1.068 0.556 -0.000 -0.000 -0.000 3/14112 790.0 1.055 1.055 0.549 -0.000 -0.000 -0.000 3/14112 800.0 1.041 1.041 0.542 -0.000 -0.000 -0.000 3/14112 810.0 1.028 1.028 0.536 -0.000 -0.000 -0.000 3/14112 820.0 1.016 1.016 0.529 -0.000 -0.000 -0.000 3/14112 830.0 1.004 1.004 0.523 -0.000 -0.000 -0.000 3/14112 840.0 0.992 0.992 0.516 -0.000 -0.000 -0.000 3/14112 850.0 0.980 0.980 0.510 -0.000 -0.000 -0.000 3/14112 860.0 0.968 0.968 0.504 -0.000 -0.000 -0.000 3/14112 870.0 0.957 0.957 0.499 -0.000 -0.000 -0.000 3/14112 880.0 0.946 0.946 0.493 -0.000 -0.000 -0.000 3/14112 890.0 0.936 0.936 0.487 -0.000 -0.000 -0.000 3/14112 900.0 0.925 0.925 0.482 -0.000 -0.000 -0.000 3/14112 910.0 0.915 0.915 0.477 -0.000 -0.000 -0.000 3/14112 920.0 0.905 0.905 0.471 -0.000 -0.000 -0.000 3/14112 930.0 0.895 0.895 0.466 -0.000 -0.000 -0.000 3/14112 940.0 0.886 0.886 0.461 -0.000 -0.000 -0.000 3/14112 950.0 0.876 0.876 0.456 -0.000 -0.000 -0.000 3/14112 960.0 0.867 0.867 0.452 -0.000 -0.000 -0.000 3/14112 970.0 0.858 0.858 0.447 -0.000 -0.000 -0.000 3/14112 980.0 0.849 0.849 0.442 -0.000 -0.000 -0.000 3/14112 990.0 0.841 0.841 0.438 -0.000 -0.000 -0.000 3/14112 1000.0 0.832 0.832 0.434 -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-07 23:11:13]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|