----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 17:34: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: P3m1 (156) Number of symmetry operations in supercell: 96 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.525459381353349 0.000000000000000 0.000000000000000 b -2.262729690676674 3.919162788046609 0.000000000000000 c 0.000000000000000 0.000000000000000 20.990822990000002 Atomic positions (fractional): *1 Pb 0.00000000000000 0.00000000000000 0.66664923299199 207.200 *2 Pb 0.33333333333333 0.66666666666667 0.33296383836945 207.200 *3 Pb 0.33333333333333 0.66666666666667 0.00039363701072 207.200 *4 I 0.66666666666667 0.33333333333333 0.42191703386653 126.904 *5 I 0.66666666666667 0.33333333333333 0.08941835397938 126.904 *6 I 0.33333333333333 0.66666666666667 0.57765044013801 126.904 *7 I 0.00000000000000 0.00000000000000 0.91142516292628 126.904 *8 I 0.00000000000000 0.00000000000000 0.24394974604188 126.904 *9 I 0.66666666666667 0.33333333333333 0.75563253467577 126.904 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.525459381353349 0.000000000000000 0.000000000000000 b -2.262729690676674 3.919162788046609 0.000000000000000 c 0.000000000000000 0.000000000000000 20.990822990000002 Atomic positions (fractional): *1 Pb 0.00000000000000 0.00000000000000 0.66664923299199 207.200 > 1 *2 Pb 0.33333333333333 0.66666666666667 0.33296383836945 207.200 > 2 *3 Pb 0.33333333333333 0.66666666666667 0.00039363701072 207.200 > 3 *4 I 0.66666666666667 0.33333333333333 0.42191703386653 126.904 > 4 *5 I 0.66666666666667 0.33333333333333 0.08941835397938 126.904 > 5 *6 I 0.33333333333333 0.66666666666667 0.57765044013801 126.904 > 6 *7 I 0.00000000000000 0.00000000000000 0.91142516292628 126.904 > 7 *8 I 0.00000000000000 0.00000000000000 0.24394974604188 126.904 > 8 *9 I 0.66666666666667 0.33333333333333 0.75563253467577 126.904 > 9 -------------------------------- super cell -------------------------------- Lattice vectors: a 18.101837525413394 0.000000000000000 0.000000000000000 b -9.050918762706695 15.676651152186436 0.000000000000000 c 0.000000000000000 0.000000000000000 20.990822990000002 Atomic positions (fractional): *1 Pb 0.00000000000000 0.00000000000000 0.66664923299199 207.200 > 1 2 Pb 0.25000000000000 0.00000000000000 0.66664923299199 207.200 > 1 3 Pb 0.50000000000000 0.00000000000000 0.66664923299199 207.200 > 1 4 Pb 0.75000000000000 0.00000000000000 0.66664923299199 207.200 > 1 5 Pb 0.00000000000000 0.25000000000000 0.66664923299199 207.200 > 1 6 Pb 0.25000000000000 0.25000000000000 0.66664923299199 207.200 > 1 7 Pb 0.50000000000000 0.25000000000000 0.66664923299199 207.200 > 1 8 Pb 0.75000000000000 0.25000000000000 0.66664923299199 207.200 > 1 9 Pb 0.00000000000000 0.50000000000000 0.66664923299199 207.200 > 1 10 Pb 0.25000000000000 0.50000000000000 0.66664923299199 207.200 > 1 11 Pb 0.50000000000000 0.50000000000000 0.66664923299199 207.200 > 1 12 Pb 0.75000000000000 0.50000000000000 0.66664923299199 207.200 > 1 13 Pb 0.00000000000000 0.75000000000000 0.66664923299199 207.200 > 1 14 Pb 0.25000000000000 0.75000000000000 0.66664923299199 207.200 > 1 15 Pb 0.50000000000000 0.75000000000000 0.66664923299199 207.200 > 1 16 Pb 0.75000000000000 0.75000000000000 0.66664923299199 207.200 > 1 *17 Pb 0.08333333333333 0.16666666666667 0.33296383836945 207.200 > 2 18 Pb 0.33333333333333 0.16666666666667 0.33296383836945 207.200 > 2 19 Pb 0.58333333333333 0.16666666666667 0.33296383836945 207.200 > 2 20 Pb 0.83333333333333 0.16666666666667 0.33296383836945 207.200 > 2 21 Pb 0.08333333333333 0.41666666666667 0.33296383836945 207.200 > 2 22 Pb 0.33333333333333 0.41666666666667 0.33296383836945 207.200 > 2 23 Pb 0.58333333333333 0.41666666666667 0.33296383836945 207.200 > 2 24 Pb 0.83333333333333 0.41666666666667 0.33296383836945 207.200 > 2 25 Pb 0.08333333333333 0.66666666666667 0.33296383836945 207.200 > 2 26 Pb 0.33333333333333 0.66666666666667 0.33296383836945 207.200 > 2 27 Pb 0.58333333333333 0.66666666666667 0.33296383836945 207.200 > 2 28 Pb 0.83333333333333 0.66666666666667 0.33296383836945 207.200 > 2 29 Pb 0.08333333333333 0.91666666666667 0.33296383836945 207.200 > 2 30 Pb 0.33333333333333 0.91666666666667 0.33296383836945 207.200 > 2 31 Pb 0.58333333333333 0.91666666666667 0.33296383836945 207.200 > 2 32 Pb 0.83333333333333 0.91666666666667 0.33296383836945 207.200 > 2 *33 Pb 0.08333333333333 0.16666666666667 0.00039363701072 207.200 > 3 34 Pb 0.33333333333333 0.16666666666667 0.00039363701072 207.200 > 3 35 Pb 0.58333333333333 0.16666666666667 0.00039363701072 207.200 > 3 36 Pb 0.83333333333333 0.16666666666667 0.00039363701072 207.200 > 3 37 Pb 0.08333333333333 0.41666666666667 0.00039363701072 207.200 > 3 38 Pb 0.33333333333333 0.41666666666667 0.00039363701072 207.200 > 3 39 Pb 0.58333333333333 0.41666666666667 0.00039363701072 207.200 > 3 40 Pb 0.83333333333333 0.41666666666667 0.00039363701072 207.200 > 3 41 Pb 0.08333333333333 0.66666666666667 0.00039363701072 207.200 > 3 42 Pb 0.33333333333333 0.66666666666667 0.00039363701072 207.200 > 3 43 Pb 0.58333333333333 0.66666666666667 0.00039363701072 207.200 > 3 44 Pb 0.83333333333333 0.66666666666667 0.00039363701072 207.200 > 3 45 Pb 0.08333333333333 0.91666666666667 0.00039363701072 207.200 > 3 46 Pb 0.33333333333333 0.91666666666667 0.00039363701072 207.200 > 3 47 Pb 0.58333333333333 0.91666666666667 0.00039363701072 207.200 > 3 48 Pb 0.83333333333333 0.91666666666667 0.00039363701072 207.200 > 3 *49 I 0.16666666666667 0.08333333333333 0.42191703386653 126.904 > 4 50 I 0.41666666666667 0.08333333333333 0.42191703386653 126.904 > 4 51 I 0.66666666666667 0.08333333333333 0.42191703386653 126.904 > 4 52 I 0.91666666666667 0.08333333333333 0.42191703386653 126.904 > 4 53 I 0.16666666666667 0.33333333333333 0.42191703386653 126.904 > 4 54 I 0.41666666666667 0.33333333333333 0.42191703386653 126.904 > 4 55 I 0.66666666666667 0.33333333333333 0.42191703386653 126.904 > 4 56 I 0.91666666666667 0.33333333333333 0.42191703386653 126.904 > 4 57 I 0.16666666666667 0.58333333333333 0.42191703386653 126.904 > 4 58 I 0.41666666666667 0.58333333333333 0.42191703386653 126.904 > 4 59 I 0.66666666666667 0.58333333333333 0.42191703386653 126.904 > 4 60 I 0.91666666666667 0.58333333333333 0.42191703386653 126.904 > 4 61 I 0.16666666666667 0.83333333333333 0.42191703386653 126.904 > 4 62 I 0.41666666666667 0.83333333333333 0.42191703386653 126.904 > 4 63 I 0.66666666666667 0.83333333333333 0.42191703386653 126.904 > 4 64 I 0.91666666666667 0.83333333333333 0.42191703386653 126.904 > 4 *65 I 0.16666666666667 0.08333333333333 0.08941835397938 126.904 > 5 66 I 0.41666666666667 0.08333333333333 0.08941835397938 126.904 > 5 67 I 0.66666666666667 0.08333333333333 0.08941835397938 126.904 > 5 68 I 0.91666666666667 0.08333333333333 0.08941835397938 126.904 > 5 69 I 0.16666666666667 0.33333333333333 0.08941835397938 126.904 > 5 70 I 0.41666666666667 0.33333333333333 0.08941835397938 126.904 > 5 71 I 0.66666666666667 0.33333333333333 0.08941835397938 126.904 > 5 72 I 0.91666666666667 0.33333333333333 0.08941835397938 126.904 > 5 73 I 0.16666666666667 0.58333333333333 0.08941835397938 126.904 > 5 74 I 0.41666666666667 0.58333333333333 0.08941835397938 126.904 > 5 75 I 0.66666666666667 0.58333333333333 0.08941835397938 126.904 > 5 76 I 0.91666666666667 0.58333333333333 0.08941835397938 126.904 > 5 77 I 0.16666666666667 0.83333333333333 0.08941835397938 126.904 > 5 78 I 0.41666666666667 0.83333333333333 0.08941835397938 126.904 > 5 79 I 0.66666666666667 0.83333333333333 0.08941835397938 126.904 > 5 80 I 0.91666666666667 0.83333333333333 0.08941835397938 126.904 > 5 *81 I 0.08333333333333 0.16666666666667 0.57765044013801 126.904 > 6 82 I 0.33333333333333 0.16666666666667 0.57765044013801 126.904 > 6 83 I 0.58333333333333 0.16666666666667 0.57765044013801 126.904 > 6 84 I 0.83333333333333 0.16666666666667 0.57765044013801 126.904 > 6 85 I 0.08333333333333 0.41666666666667 0.57765044013801 126.904 > 6 86 I 0.33333333333333 0.41666666666667 0.57765044013801 126.904 > 6 87 I 0.58333333333333 0.41666666666667 0.57765044013801 126.904 > 6 88 I 0.83333333333333 0.41666666666667 0.57765044013801 126.904 > 6 89 I 0.08333333333333 0.66666666666667 0.57765044013801 126.904 > 6 90 I 0.33333333333333 0.66666666666667 0.57765044013801 126.904 > 6 91 I 0.58333333333333 0.66666666666667 0.57765044013801 126.904 > 6 92 I 0.83333333333333 0.66666666666667 0.57765044013801 126.904 > 6 93 I 0.08333333333333 0.91666666666667 0.57765044013801 126.904 > 6 94 I 0.33333333333333 0.91666666666667 0.57765044013801 126.904 > 6 95 I 0.58333333333333 0.91666666666667 0.57765044013801 126.904 > 6 96 I 0.83333333333333 0.91666666666667 0.57765044013801 126.904 > 6 *97 I 0.00000000000000 0.00000000000000 0.91142516292628 126.904 > 7 98 I 0.25000000000000 0.00000000000000 0.91142516292628 126.904 > 7 99 I 0.50000000000000 0.00000000000000 0.91142516292628 126.904 > 7 100 I 0.75000000000000 0.00000000000000 0.91142516292628 126.904 > 7 101 I 0.00000000000000 0.25000000000000 0.91142516292628 126.904 > 7 102 I 0.25000000000000 0.25000000000000 0.91142516292628 126.904 > 7 103 I 0.50000000000000 0.25000000000000 0.91142516292628 126.904 > 7 104 I 0.75000000000000 0.25000000000000 0.91142516292628 126.904 > 7 105 I 0.00000000000000 0.50000000000000 0.91142516292628 126.904 > 7 106 I 0.25000000000000 0.50000000000000 0.91142516292628 126.904 > 7 107 I 0.50000000000000 0.50000000000000 0.91142516292628 126.904 > 7 108 I 0.75000000000000 0.50000000000000 0.91142516292628 126.904 > 7 109 I 0.00000000000000 0.75000000000000 0.91142516292628 126.904 > 7 110 I 0.25000000000000 0.75000000000000 0.91142516292628 126.904 > 7 111 I 0.50000000000000 0.75000000000000 0.91142516292628 126.904 > 7 112 I 0.75000000000000 0.75000000000000 0.91142516292628 126.904 > 7 *113 I 0.00000000000000 0.00000000000000 0.24394974604188 126.904 > 8 114 I 0.25000000000000 0.00000000000000 0.24394974604188 126.904 > 8 115 I 0.50000000000000 0.00000000000000 0.24394974604188 126.904 > 8 116 I 0.75000000000000 0.00000000000000 0.24394974604188 126.904 > 8 117 I 0.00000000000000 0.25000000000000 0.24394974604188 126.904 > 8 118 I 0.25000000000000 0.25000000000000 0.24394974604188 126.904 > 8 119 I 0.50000000000000 0.25000000000000 0.24394974604188 126.904 > 8 120 I 0.75000000000000 0.25000000000000 0.24394974604188 126.904 > 8 121 I 0.00000000000000 0.50000000000000 0.24394974604188 126.904 > 8 122 I 0.25000000000000 0.50000000000000 0.24394974604188 126.904 > 8 123 I 0.50000000000000 0.50000000000000 0.24394974604188 126.904 > 8 124 I 0.75000000000000 0.50000000000000 0.24394974604188 126.904 > 8 125 I 0.00000000000000 0.75000000000000 0.24394974604188 126.904 > 8 126 I 0.25000000000000 0.75000000000000 0.24394974604188 126.904 > 8 127 I 0.50000000000000 0.75000000000000 0.24394974604188 126.904 > 8 128 I 0.75000000000000 0.75000000000000 0.24394974604188 126.904 > 8 *129 I 0.16666666666667 0.08333333333333 0.75563253467577 126.904 > 9 130 I 0.41666666666667 0.08333333333333 0.75563253467577 126.904 > 9 131 I 0.66666666666667 0.08333333333333 0.75563253467577 126.904 > 9 132 I 0.91666666666667 0.08333333333333 0.75563253467577 126.904 > 9 133 I 0.16666666666667 0.33333333333333 0.75563253467577 126.904 > 9 134 I 0.41666666666667 0.33333333333333 0.75563253467577 126.904 > 9 135 I 0.66666666666667 0.33333333333333 0.75563253467577 126.904 > 9 136 I 0.91666666666667 0.33333333333333 0.75563253467577 126.904 > 9 137 I 0.16666666666667 0.58333333333333 0.75563253467577 126.904 > 9 138 I 0.41666666666667 0.58333333333333 0.75563253467577 126.904 > 9 139 I 0.66666666666667 0.58333333333333 0.75563253467577 126.904 > 9 140 I 0.91666666666667 0.58333333333333 0.75563253467577 126.904 > 9 141 I 0.16666666666667 0.83333333333333 0.75563253467577 126.904 > 9 142 I 0.41666666666667 0.83333333333333 0.75563253467577 126.904 > 9 143 I 0.66666666666667 0.83333333333333 0.75563253467577 126.904 > 9 144 I 0.91666666666667 0.83333333333333 0.75563253467577 126.904 > 9 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.0485280 -0.0000000 0.0000000 -0.0000000 7.0485280 0.0000000 0.0000000 0.0000000 5.2609756 -------------------------- Born effective charges -------------------------- 1 Pb 4.0617570 -0.0000000 0.0000000 0.0000000 4.0617570 0.0000000 0.0000000 0.0000000 1.8421182 2 Pb 4.0693843 -0.0000000 0.0000000 0.0000000 4.0693843 0.0000000 0.0000000 0.0000000 1.8495626 3 Pb 4.0973248 -0.0000000 0.0000000 -0.0000000 4.0973248 0.0000000 0.0000000 0.0000000 1.8249840 4 I -2.0522570 0.0000000 0.0000000 -0.0000000 -2.0522570 0.0000000 0.0000000 0.0000000 -0.9382789 5 I -2.0408895 0.0000000 0.0000000 -0.0000000 -2.0408895 0.0000000 0.0000000 0.0000000 -0.9132740 6 I -2.0206285 0.0000000 0.0000000 0.0000000 -2.0206285 0.0000000 0.0000000 0.0000000 -0.9093271 7 I -2.0333309 0.0000000 0.0000000 -0.0000000 -2.0333309 0.0000000 0.0000000 0.0000000 -0.9081123 8 I -2.0314621 0.0000000 0.0000000 -0.0000000 -2.0314621 0.0000000 0.0000000 0.0000000 -0.9133378 9 I -2.0498979 0.0000000 0.0000000 0.0000000 -2.0498979 0.0000000 0.0000000 0.0000000 -0.9343347 ---------------------------------------------------------------------------- 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: 432/432 Permutation basis: 6426/6426 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 1089 Number of blocks in projector: 837 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: 468 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 405 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 216 Use standard eigh solver. Tree of FC basis block matrices: - (1089, 1071), data: False |-- (216, 207), data: True |-- (405, 405), data: True |-- (468, 459), data: True ----- Solver_atoms: 1 -- 144 / 144 Time (Solver_compr_matrix_reshape): 0.016 Solver_block: 100 / 240 - Time: 1.476 Solver_block: 200 / 240 - Time: 1.478 Solver_block: 240 / 240 - Time: 0.683 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 3.673 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 432/432 Permutation basis: 6426/6426 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 1089 Number of blocks in projector: 837 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: 468 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 405 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 216 Use standard eigh solver. Tree of FC basis block matrices: - (1089, 1071), data: False |-- (216, 207), data: True |-- (405, 405), data: True |-- (468, 459), 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-08 17:35:06]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 17:35:06]------------------------- 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: P3m1 (156) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.525459381353349 0.000000000000000 0.000000000000000 b -2.262729690676674 3.919162788046609 0.000000000000000 c 0.000000000000000 0.000000000000000 20.990822990000002 Atomic positions (fractional): 1 Pb 0.00000000000000 0.00000000000000 0.66664923299199 207.200 2 Pb 0.33333333333333 0.66666666666667 0.33296383836945 207.200 3 Pb 0.33333333333333 0.66666666666667 0.00039363701072 207.200 4 I 0.66666666666667 0.33333333333333 0.42191703386653 126.904 5 I 0.66666666666667 0.33333333333333 0.08941835397938 126.904 6 I 0.33333333333333 0.66666666666667 0.57765044013801 126.904 7 I 0.00000000000000 0.00000000000000 0.91142516292628 126.904 8 I 0.00000000000000 0.00000000000000 0.24394974604188 126.904 9 I 0.66666666666667 0.33333333333333 0.75563253467577 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 18.101837525413394 0.000000000000000 0.000000000000000 b -9.050918762706695 15.676651152186436 0.000000000000000 c 0.000000000000000 0.000000000000000 20.990822990000002 Atomic positions (fractional): 1 Pb 0.00000000000000 0.00000000000000 0.66664923299199 207.200 > 1 2 Pb 0.25000000000000 0.00000000000000 0.66664923299199 207.200 > 1 3 Pb 0.50000000000000 0.00000000000000 0.66664923299199 207.200 > 1 4 Pb 0.75000000000000 0.00000000000000 0.66664923299199 207.200 > 1 5 Pb 0.00000000000000 0.25000000000000 0.66664923299199 207.200 > 1 6 Pb 0.25000000000000 0.25000000000000 0.66664923299199 207.200 > 1 7 Pb 0.50000000000000 0.25000000000000 0.66664923299199 207.200 > 1 8 Pb 0.75000000000000 0.25000000000000 0.66664923299199 207.200 > 1 9 Pb 0.00000000000000 0.50000000000000 0.66664923299199 207.200 > 1 10 Pb 0.25000000000000 0.50000000000000 0.66664923299199 207.200 > 1 11 Pb 0.50000000000000 0.50000000000000 0.66664923299199 207.200 > 1 12 Pb 0.75000000000000 0.50000000000000 0.66664923299199 207.200 > 1 13 Pb 0.00000000000000 0.75000000000000 0.66664923299199 207.200 > 1 14 Pb 0.25000000000000 0.75000000000000 0.66664923299199 207.200 > 1 15 Pb 0.50000000000000 0.75000000000000 0.66664923299199 207.200 > 1 16 Pb 0.75000000000000 0.75000000000000 0.66664923299199 207.200 > 1 17 Pb 0.08333333333333 0.16666666666667 0.33296383836945 207.200 > 17 18 Pb 0.33333333333333 0.16666666666667 0.33296383836945 207.200 > 17 19 Pb 0.58333333333333 0.16666666666667 0.33296383836945 207.200 > 17 20 Pb 0.83333333333333 0.16666666666667 0.33296383836945 207.200 > 17 21 Pb 0.08333333333333 0.41666666666667 0.33296383836945 207.200 > 17 22 Pb 0.33333333333333 0.41666666666667 0.33296383836945 207.200 > 17 23 Pb 0.58333333333333 0.41666666666667 0.33296383836945 207.200 > 17 24 Pb 0.83333333333333 0.41666666666667 0.33296383836945 207.200 > 17 25 Pb 0.08333333333333 0.66666666666667 0.33296383836945 207.200 > 17 26 Pb 0.33333333333333 0.66666666666667 0.33296383836945 207.200 > 17 27 Pb 0.58333333333333 0.66666666666667 0.33296383836945 207.200 > 17 28 Pb 0.83333333333333 0.66666666666667 0.33296383836945 207.200 > 17 29 Pb 0.08333333333333 0.91666666666667 0.33296383836945 207.200 > 17 30 Pb 0.33333333333333 0.91666666666667 0.33296383836945 207.200 > 17 31 Pb 0.58333333333333 0.91666666666667 0.33296383836945 207.200 > 17 32 Pb 0.83333333333333 0.91666666666667 0.33296383836945 207.200 > 17 33 Pb 0.08333333333333 0.16666666666667 0.00039363701072 207.200 > 33 34 Pb 0.33333333333333 0.16666666666667 0.00039363701072 207.200 > 33 35 Pb 0.58333333333333 0.16666666666667 0.00039363701072 207.200 > 33 36 Pb 0.83333333333333 0.16666666666667 0.00039363701072 207.200 > 33 37 Pb 0.08333333333333 0.41666666666667 0.00039363701072 207.200 > 33 38 Pb 0.33333333333333 0.41666666666667 0.00039363701072 207.200 > 33 39 Pb 0.58333333333333 0.41666666666667 0.00039363701072 207.200 > 33 40 Pb 0.83333333333333 0.41666666666667 0.00039363701072 207.200 > 33 41 Pb 0.08333333333333 0.66666666666667 0.00039363701072 207.200 > 33 42 Pb 0.33333333333333 0.66666666666667 0.00039363701072 207.200 > 33 43 Pb 0.58333333333333 0.66666666666667 0.00039363701072 207.200 > 33 44 Pb 0.83333333333333 0.66666666666667 0.00039363701072 207.200 > 33 45 Pb 0.08333333333333 0.91666666666667 0.00039363701072 207.200 > 33 46 Pb 0.33333333333333 0.91666666666667 0.00039363701072 207.200 > 33 47 Pb 0.58333333333333 0.91666666666667 0.00039363701072 207.200 > 33 48 Pb 0.83333333333333 0.91666666666667 0.00039363701072 207.200 > 33 49 I 0.16666666666667 0.08333333333333 0.42191703386653 126.904 > 49 50 I 0.41666666666667 0.08333333333333 0.42191703386653 126.904 > 49 51 I 0.66666666666667 0.08333333333333 0.42191703386653 126.904 > 49 52 I 0.91666666666667 0.08333333333333 0.42191703386653 126.904 > 49 53 I 0.16666666666667 0.33333333333333 0.42191703386653 126.904 > 49 54 I 0.41666666666667 0.33333333333333 0.42191703386653 126.904 > 49 55 I 0.66666666666667 0.33333333333333 0.42191703386653 126.904 > 49 56 I 0.91666666666667 0.33333333333333 0.42191703386653 126.904 > 49 57 I 0.16666666666667 0.58333333333333 0.42191703386653 126.904 > 49 58 I 0.41666666666667 0.58333333333333 0.42191703386653 126.904 > 49 59 I 0.66666666666667 0.58333333333333 0.42191703386653 126.904 > 49 60 I 0.91666666666667 0.58333333333333 0.42191703386653 126.904 > 49 61 I 0.16666666666667 0.83333333333333 0.42191703386653 126.904 > 49 62 I 0.41666666666667 0.83333333333333 0.42191703386653 126.904 > 49 63 I 0.66666666666667 0.83333333333333 0.42191703386653 126.904 > 49 64 I 0.91666666666667 0.83333333333333 0.42191703386653 126.904 > 49 65 I 0.16666666666667 0.08333333333333 0.08941835397938 126.904 > 65 66 I 0.41666666666667 0.08333333333333 0.08941835397938 126.904 > 65 67 I 0.66666666666667 0.08333333333333 0.08941835397938 126.904 > 65 68 I 0.91666666666667 0.08333333333333 0.08941835397938 126.904 > 65 69 I 0.16666666666667 0.33333333333333 0.08941835397938 126.904 > 65 70 I 0.41666666666667 0.33333333333333 0.08941835397938 126.904 > 65 71 I 0.66666666666667 0.33333333333333 0.08941835397938 126.904 > 65 72 I 0.91666666666667 0.33333333333333 0.08941835397938 126.904 > 65 73 I 0.16666666666667 0.58333333333333 0.08941835397938 126.904 > 65 74 I 0.41666666666667 0.58333333333333 0.08941835397938 126.904 > 65 75 I 0.66666666666667 0.58333333333333 0.08941835397938 126.904 > 65 76 I 0.91666666666667 0.58333333333333 0.08941835397938 126.904 > 65 77 I 0.16666666666667 0.83333333333333 0.08941835397938 126.904 > 65 78 I 0.41666666666667 0.83333333333333 0.08941835397938 126.904 > 65 79 I 0.66666666666667 0.83333333333333 0.08941835397938 126.904 > 65 80 I 0.91666666666667 0.83333333333333 0.08941835397938 126.904 > 65 81 I 0.08333333333333 0.16666666666667 0.57765044013801 126.904 > 81 82 I 0.33333333333333 0.16666666666667 0.57765044013801 126.904 > 81 83 I 0.58333333333333 0.16666666666667 0.57765044013801 126.904 > 81 84 I 0.83333333333333 0.16666666666667 0.57765044013801 126.904 > 81 85 I 0.08333333333333 0.41666666666667 0.57765044013801 126.904 > 81 86 I 0.33333333333333 0.41666666666667 0.57765044013801 126.904 > 81 87 I 0.58333333333333 0.41666666666667 0.57765044013801 126.904 > 81 88 I 0.83333333333333 0.41666666666667 0.57765044013801 126.904 > 81 89 I 0.08333333333333 0.66666666666667 0.57765044013801 126.904 > 81 90 I 0.33333333333333 0.66666666666667 0.57765044013801 126.904 > 81 91 I 0.58333333333333 0.66666666666667 0.57765044013801 126.904 > 81 92 I 0.83333333333333 0.66666666666667 0.57765044013801 126.904 > 81 93 I 0.08333333333333 0.91666666666667 0.57765044013801 126.904 > 81 94 I 0.33333333333333 0.91666666666667 0.57765044013801 126.904 > 81 95 I 0.58333333333333 0.91666666666667 0.57765044013801 126.904 > 81 96 I 0.83333333333333 0.91666666666667 0.57765044013801 126.904 > 81 97 I 0.00000000000000 0.00000000000000 0.91142516292628 126.904 > 97 98 I 0.25000000000000 0.00000000000000 0.91142516292628 126.904 > 97 99 I 0.50000000000000 0.00000000000000 0.91142516292628 126.904 > 97 100 I 0.75000000000000 0.00000000000000 0.91142516292628 126.904 > 97 101 I 0.00000000000000 0.25000000000000 0.91142516292628 126.904 > 97 102 I 0.25000000000000 0.25000000000000 0.91142516292628 126.904 > 97 103 I 0.50000000000000 0.25000000000000 0.91142516292628 126.904 > 97 104 I 0.75000000000000 0.25000000000000 0.91142516292628 126.904 > 97 105 I 0.00000000000000 0.50000000000000 0.91142516292628 126.904 > 97 106 I 0.25000000000000 0.50000000000000 0.91142516292628 126.904 > 97 107 I 0.50000000000000 0.50000000000000 0.91142516292628 126.904 > 97 108 I 0.75000000000000 0.50000000000000 0.91142516292628 126.904 > 97 109 I 0.00000000000000 0.75000000000000 0.91142516292628 126.904 > 97 110 I 0.25000000000000 0.75000000000000 0.91142516292628 126.904 > 97 111 I 0.50000000000000 0.75000000000000 0.91142516292628 126.904 > 97 112 I 0.75000000000000 0.75000000000000 0.91142516292628 126.904 > 97 113 I 0.00000000000000 0.00000000000000 0.24394974604188 126.904 > 113 114 I 0.25000000000000 0.00000000000000 0.24394974604188 126.904 > 113 115 I 0.50000000000000 0.00000000000000 0.24394974604188 126.904 > 113 116 I 0.75000000000000 0.00000000000000 0.24394974604188 126.904 > 113 117 I 0.00000000000000 0.25000000000000 0.24394974604188 126.904 > 113 118 I 0.25000000000000 0.25000000000000 0.24394974604188 126.904 > 113 119 I 0.50000000000000 0.25000000000000 0.24394974604188 126.904 > 113 120 I 0.75000000000000 0.25000000000000 0.24394974604188 126.904 > 113 121 I 0.00000000000000 0.50000000000000 0.24394974604188 126.904 > 113 122 I 0.25000000000000 0.50000000000000 0.24394974604188 126.904 > 113 123 I 0.50000000000000 0.50000000000000 0.24394974604188 126.904 > 113 124 I 0.75000000000000 0.50000000000000 0.24394974604188 126.904 > 113 125 I 0.00000000000000 0.75000000000000 0.24394974604188 126.904 > 113 126 I 0.25000000000000 0.75000000000000 0.24394974604188 126.904 > 113 127 I 0.50000000000000 0.75000000000000 0.24394974604188 126.904 > 113 128 I 0.75000000000000 0.75000000000000 0.24394974604188 126.904 > 113 129 I 0.16666666666667 0.08333333333333 0.75563253467577 126.904 > 129 130 I 0.41666666666667 0.08333333333333 0.75563253467577 126.904 > 129 131 I 0.66666666666667 0.08333333333333 0.75563253467577 126.904 > 129 132 I 0.91666666666667 0.08333333333333 0.75563253467577 126.904 > 129 133 I 0.16666666666667 0.33333333333333 0.75563253467577 126.904 > 129 134 I 0.41666666666667 0.33333333333333 0.75563253467577 126.904 > 129 135 I 0.66666666666667 0.33333333333333 0.75563253467577 126.904 > 129 136 I 0.91666666666667 0.33333333333333 0.75563253467577 126.904 > 129 137 I 0.16666666666667 0.58333333333333 0.75563253467577 126.904 > 129 138 I 0.41666666666667 0.58333333333333 0.75563253467577 126.904 > 129 139 I 0.66666666666667 0.58333333333333 0.75563253467577 126.904 > 129 140 I 0.91666666666667 0.58333333333333 0.75563253467577 126.904 > 129 141 I 0.16666666666667 0.83333333333333 0.75563253467577 126.904 > 129 142 I 0.41666666666667 0.83333333333333 0.75563253467577 126.904 > 129 143 I 0.66666666666667 0.83333333333333 0.75563253467577 126.904 > 129 144 I 0.91666666666667 0.83333333333333 0.75563253467577 126.904 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 7.0485280 -0.0000000 0.0000000 -0.0000000 7.0485280 0.0000000 0.0000000 0.0000000 5.2609756 -------------------------- Born effective charges -------------------------- 1 Pb 4.0617570 -0.0000000 0.0000000 0.0000000 4.0617570 0.0000000 0.0000000 0.0000000 1.8421182 2 Pb 4.0693843 -0.0000000 0.0000000 0.0000000 4.0693843 0.0000000 0.0000000 0.0000000 1.8495626 3 Pb 4.0973248 -0.0000000 0.0000000 -0.0000000 4.0973248 0.0000000 0.0000000 0.0000000 1.8249840 4 I -2.0522570 0.0000000 0.0000000 -0.0000000 -2.0522570 0.0000000 0.0000000 0.0000000 -0.9382789 5 I -2.0408895 0.0000000 0.0000000 -0.0000000 -2.0408895 0.0000000 0.0000000 0.0000000 -0.9132740 6 I -2.0206285 0.0000000 0.0000000 0.0000000 -2.0206285 0.0000000 0.0000000 0.0000000 -0.9093271 7 I -2.0333309 0.0000000 0.0000000 -0.0000000 -2.0333309 0.0000000 0.0000000 0.0000000 -0.9081123 8 I -2.0314621 0.0000000 0.0000000 -0.0000000 -2.0314621 0.0000000 0.0000000 0.0000000 -0.9133378 9 I -2.0498979 0.0000000 0.0000000 0.0000000 -2.0498979 0.0000000 0.0000000 0.0000000 -0.9343347 ---------------------------------------------------------------------------- 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[ 17, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 33, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 49, 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[ 81, 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] Computing fc3[ 113, 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[ 129, 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.00000007 (yyy) 0.00000007 (yyy) 0.00000007 (yyy) 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-08 17:35:31]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 17:35:31]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [4 4 1] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P3m1 (156) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.525459381353349 0.000000000000000 0.000000000000000 b -2.262729690676674 3.919162788046609 0.000000000000000 c 0.000000000000000 0.000000000000000 20.990822990000002 Atomic positions (fractional): 1 Pb 0.00000000000000 0.00000000000000 0.66664923299199 207.200 2 Pb 0.33333333333333 0.66666666666667 0.33296383836945 207.200 3 Pb 0.33333333333333 0.66666666666667 0.00039363701072 207.200 4 I 0.66666666666667 0.33333333333333 0.42191703386653 126.904 5 I 0.66666666666667 0.33333333333333 0.08941835397938 126.904 6 I 0.33333333333333 0.66666666666667 0.57765044013801 126.904 7 I 0.00000000000000 0.00000000000000 0.91142516292628 126.904 8 I 0.00000000000000 0.00000000000000 0.24394974604188 126.904 9 I 0.66666666666667 0.33333333333333 0.75563253467577 126.904 -------------------------------- supercell --------------------------------- Lattice vectors: a 18.101837525413394 0.000000000000000 0.000000000000000 b -9.050918762706695 15.676651152186436 0.000000000000000 c 0.000000000000000 0.000000000000000 20.990822990000002 Atomic positions (fractional): 1 Pb 0.00000000000000 0.00000000000000 0.66664923299199 207.200 > 1 2 Pb 0.25000000000000 0.00000000000000 0.66664923299199 207.200 > 1 3 Pb 0.50000000000000 0.00000000000000 0.66664923299199 207.200 > 1 4 Pb 0.75000000000000 0.00000000000000 0.66664923299199 207.200 > 1 5 Pb 0.00000000000000 0.25000000000000 0.66664923299199 207.200 > 1 6 Pb 0.25000000000000 0.25000000000000 0.66664923299199 207.200 > 1 7 Pb 0.50000000000000 0.25000000000000 0.66664923299199 207.200 > 1 8 Pb 0.75000000000000 0.25000000000000 0.66664923299199 207.200 > 1 9 Pb 0.00000000000000 0.50000000000000 0.66664923299199 207.200 > 1 10 Pb 0.25000000000000 0.50000000000000 0.66664923299199 207.200 > 1 11 Pb 0.50000000000000 0.50000000000000 0.66664923299199 207.200 > 1 12 Pb 0.75000000000000 0.50000000000000 0.66664923299199 207.200 > 1 13 Pb 0.00000000000000 0.75000000000000 0.66664923299199 207.200 > 1 14 Pb 0.25000000000000 0.75000000000000 0.66664923299199 207.200 > 1 15 Pb 0.50000000000000 0.75000000000000 0.66664923299199 207.200 > 1 16 Pb 0.75000000000000 0.75000000000000 0.66664923299199 207.200 > 1 17 Pb 0.08333333333333 0.16666666666667 0.33296383836945 207.200 > 17 18 Pb 0.33333333333333 0.16666666666667 0.33296383836945 207.200 > 17 19 Pb 0.58333333333333 0.16666666666667 0.33296383836945 207.200 > 17 20 Pb 0.83333333333333 0.16666666666667 0.33296383836945 207.200 > 17 21 Pb 0.08333333333333 0.41666666666667 0.33296383836945 207.200 > 17 22 Pb 0.33333333333333 0.41666666666667 0.33296383836945 207.200 > 17 23 Pb 0.58333333333333 0.41666666666667 0.33296383836945 207.200 > 17 24 Pb 0.83333333333333 0.41666666666667 0.33296383836945 207.200 > 17 25 Pb 0.08333333333333 0.66666666666667 0.33296383836945 207.200 > 17 26 Pb 0.33333333333333 0.66666666666667 0.33296383836945 207.200 > 17 27 Pb 0.58333333333333 0.66666666666667 0.33296383836945 207.200 > 17 28 Pb 0.83333333333333 0.66666666666667 0.33296383836945 207.200 > 17 29 Pb 0.08333333333333 0.91666666666667 0.33296383836945 207.200 > 17 30 Pb 0.33333333333333 0.91666666666667 0.33296383836945 207.200 > 17 31 Pb 0.58333333333333 0.91666666666667 0.33296383836945 207.200 > 17 32 Pb 0.83333333333333 0.91666666666667 0.33296383836945 207.200 > 17 33 Pb 0.08333333333333 0.16666666666667 0.00039363701072 207.200 > 33 34 Pb 0.33333333333333 0.16666666666667 0.00039363701072 207.200 > 33 35 Pb 0.58333333333333 0.16666666666667 0.00039363701072 207.200 > 33 36 Pb 0.83333333333333 0.16666666666667 0.00039363701072 207.200 > 33 37 Pb 0.08333333333333 0.41666666666667 0.00039363701072 207.200 > 33 38 Pb 0.33333333333333 0.41666666666667 0.00039363701072 207.200 > 33 39 Pb 0.58333333333333 0.41666666666667 0.00039363701072 207.200 > 33 40 Pb 0.83333333333333 0.41666666666667 0.00039363701072 207.200 > 33 41 Pb 0.08333333333333 0.66666666666667 0.00039363701072 207.200 > 33 42 Pb 0.33333333333333 0.66666666666667 0.00039363701072 207.200 > 33 43 Pb 0.58333333333333 0.66666666666667 0.00039363701072 207.200 > 33 44 Pb 0.83333333333333 0.66666666666667 0.00039363701072 207.200 > 33 45 Pb 0.08333333333333 0.91666666666667 0.00039363701072 207.200 > 33 46 Pb 0.33333333333333 0.91666666666667 0.00039363701072 207.200 > 33 47 Pb 0.58333333333333 0.91666666666667 0.00039363701072 207.200 > 33 48 Pb 0.83333333333333 0.91666666666667 0.00039363701072 207.200 > 33 49 I 0.16666666666667 0.08333333333333 0.42191703386653 126.904 > 49 50 I 0.41666666666667 0.08333333333333 0.42191703386653 126.904 > 49 51 I 0.66666666666667 0.08333333333333 0.42191703386653 126.904 > 49 52 I 0.91666666666667 0.08333333333333 0.42191703386653 126.904 > 49 53 I 0.16666666666667 0.33333333333333 0.42191703386653 126.904 > 49 54 I 0.41666666666667 0.33333333333333 0.42191703386653 126.904 > 49 55 I 0.66666666666667 0.33333333333333 0.42191703386653 126.904 > 49 56 I 0.91666666666667 0.33333333333333 0.42191703386653 126.904 > 49 57 I 0.16666666666667 0.58333333333333 0.42191703386653 126.904 > 49 58 I 0.41666666666667 0.58333333333333 0.42191703386653 126.904 > 49 59 I 0.66666666666667 0.58333333333333 0.42191703386653 126.904 > 49 60 I 0.91666666666667 0.58333333333333 0.42191703386653 126.904 > 49 61 I 0.16666666666667 0.83333333333333 0.42191703386653 126.904 > 49 62 I 0.41666666666667 0.83333333333333 0.42191703386653 126.904 > 49 63 I 0.66666666666667 0.83333333333333 0.42191703386653 126.904 > 49 64 I 0.91666666666667 0.83333333333333 0.42191703386653 126.904 > 49 65 I 0.16666666666667 0.08333333333333 0.08941835397938 126.904 > 65 66 I 0.41666666666667 0.08333333333333 0.08941835397938 126.904 > 65 67 I 0.66666666666667 0.08333333333333 0.08941835397938 126.904 > 65 68 I 0.91666666666667 0.08333333333333 0.08941835397938 126.904 > 65 69 I 0.16666666666667 0.33333333333333 0.08941835397938 126.904 > 65 70 I 0.41666666666667 0.33333333333333 0.08941835397938 126.904 > 65 71 I 0.66666666666667 0.33333333333333 0.08941835397938 126.904 > 65 72 I 0.91666666666667 0.33333333333333 0.08941835397938 126.904 > 65 73 I 0.16666666666667 0.58333333333333 0.08941835397938 126.904 > 65 74 I 0.41666666666667 0.58333333333333 0.08941835397938 126.904 > 65 75 I 0.66666666666667 0.58333333333333 0.08941835397938 126.904 > 65 76 I 0.91666666666667 0.58333333333333 0.08941835397938 126.904 > 65 77 I 0.16666666666667 0.83333333333333 0.08941835397938 126.904 > 65 78 I 0.41666666666667 0.83333333333333 0.08941835397938 126.904 > 65 79 I 0.66666666666667 0.83333333333333 0.08941835397938 126.904 > 65 80 I 0.91666666666667 0.83333333333333 0.08941835397938 126.904 > 65 81 I 0.08333333333333 0.16666666666667 0.57765044013801 126.904 > 81 82 I 0.33333333333333 0.16666666666667 0.57765044013801 126.904 > 81 83 I 0.58333333333333 0.16666666666667 0.57765044013801 126.904 > 81 84 I 0.83333333333333 0.16666666666667 0.57765044013801 126.904 > 81 85 I 0.08333333333333 0.41666666666667 0.57765044013801 126.904 > 81 86 I 0.33333333333333 0.41666666666667 0.57765044013801 126.904 > 81 87 I 0.58333333333333 0.41666666666667 0.57765044013801 126.904 > 81 88 I 0.83333333333333 0.41666666666667 0.57765044013801 126.904 > 81 89 I 0.08333333333333 0.66666666666667 0.57765044013801 126.904 > 81 90 I 0.33333333333333 0.66666666666667 0.57765044013801 126.904 > 81 91 I 0.58333333333333 0.66666666666667 0.57765044013801 126.904 > 81 92 I 0.83333333333333 0.66666666666667 0.57765044013801 126.904 > 81 93 I 0.08333333333333 0.91666666666667 0.57765044013801 126.904 > 81 94 I 0.33333333333333 0.91666666666667 0.57765044013801 126.904 > 81 95 I 0.58333333333333 0.91666666666667 0.57765044013801 126.904 > 81 96 I 0.83333333333333 0.91666666666667 0.57765044013801 126.904 > 81 97 I 0.00000000000000 0.00000000000000 0.91142516292628 126.904 > 97 98 I 0.25000000000000 0.00000000000000 0.91142516292628 126.904 > 97 99 I 0.50000000000000 0.00000000000000 0.91142516292628 126.904 > 97 100 I 0.75000000000000 0.00000000000000 0.91142516292628 126.904 > 97 101 I 0.00000000000000 0.25000000000000 0.91142516292628 126.904 > 97 102 I 0.25000000000000 0.25000000000000 0.91142516292628 126.904 > 97 103 I 0.50000000000000 0.25000000000000 0.91142516292628 126.904 > 97 104 I 0.75000000000000 0.25000000000000 0.91142516292628 126.904 > 97 105 I 0.00000000000000 0.50000000000000 0.91142516292628 126.904 > 97 106 I 0.25000000000000 0.50000000000000 0.91142516292628 126.904 > 97 107 I 0.50000000000000 0.50000000000000 0.91142516292628 126.904 > 97 108 I 0.75000000000000 0.50000000000000 0.91142516292628 126.904 > 97 109 I 0.00000000000000 0.75000000000000 0.91142516292628 126.904 > 97 110 I 0.25000000000000 0.75000000000000 0.91142516292628 126.904 > 97 111 I 0.50000000000000 0.75000000000000 0.91142516292628 126.904 > 97 112 I 0.75000000000000 0.75000000000000 0.91142516292628 126.904 > 97 113 I 0.00000000000000 0.00000000000000 0.24394974604188 126.904 > 113 114 I 0.25000000000000 0.00000000000000 0.24394974604188 126.904 > 113 115 I 0.50000000000000 0.00000000000000 0.24394974604188 126.904 > 113 116 I 0.75000000000000 0.00000000000000 0.24394974604188 126.904 > 113 117 I 0.00000000000000 0.25000000000000 0.24394974604188 126.904 > 113 118 I 0.25000000000000 0.25000000000000 0.24394974604188 126.904 > 113 119 I 0.50000000000000 0.25000000000000 0.24394974604188 126.904 > 113 120 I 0.75000000000000 0.25000000000000 0.24394974604188 126.904 > 113 121 I 0.00000000000000 0.50000000000000 0.24394974604188 126.904 > 113 122 I 0.25000000000000 0.50000000000000 0.24394974604188 126.904 > 113 123 I 0.50000000000000 0.50000000000000 0.24394974604188 126.904 > 113 124 I 0.75000000000000 0.50000000000000 0.24394974604188 126.904 > 113 125 I 0.00000000000000 0.75000000000000 0.24394974604188 126.904 > 113 126 I 0.25000000000000 0.75000000000000 0.24394974604188 126.904 > 113 127 I 0.50000000000000 0.75000000000000 0.24394974604188 126.904 > 113 128 I 0.75000000000000 0.75000000000000 0.24394974604188 126.904 > 113 129 I 0.16666666666667 0.08333333333333 0.75563253467577 126.904 > 129 130 I 0.41666666666667 0.08333333333333 0.75563253467577 126.904 > 129 131 I 0.66666666666667 0.08333333333333 0.75563253467577 126.904 > 129 132 I 0.91666666666667 0.08333333333333 0.75563253467577 126.904 > 129 133 I 0.16666666666667 0.33333333333333 0.75563253467577 126.904 > 129 134 I 0.41666666666667 0.33333333333333 0.75563253467577 126.904 > 129 135 I 0.66666666666667 0.33333333333333 0.75563253467577 126.904 > 129 136 I 0.91666666666667 0.33333333333333 0.75563253467577 126.904 > 129 137 I 0.16666666666667 0.58333333333333 0.75563253467577 126.904 > 129 138 I 0.41666666666667 0.58333333333333 0.75563253467577 126.904 > 129 139 I 0.66666666666667 0.58333333333333 0.75563253467577 126.904 > 129 140 I 0.91666666666667 0.58333333333333 0.75563253467577 126.904 > 129 141 I 0.16666666666667 0.83333333333333 0.75563253467577 126.904 > 129 142 I 0.41666666666667 0.83333333333333 0.75563253467577 126.904 > 129 143 I 0.66666666666667 0.83333333333333 0.75563253467577 126.904 > 129 144 I 0.91666666666667 0.83333333333333 0.75563253467577 126.904 > 129 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 7.0485280 -0.0000000 0.0000000 -0.0000000 7.0485280 0.0000000 0.0000000 0.0000000 5.2609756 -------------------------- Born effective charges -------------------------- 1 Pb 4.0617570 -0.0000000 0.0000000 0.0000000 4.0617570 0.0000000 0.0000000 0.0000000 1.8421182 2 Pb 4.0693843 -0.0000000 0.0000000 0.0000000 4.0693843 0.0000000 0.0000000 0.0000000 1.8495626 3 Pb 4.0973248 -0.0000000 0.0000000 -0.0000000 4.0973248 0.0000000 0.0000000 0.0000000 1.8249840 4 I -2.0522570 0.0000000 0.0000000 -0.0000000 -2.0522570 0.0000000 0.0000000 0.0000000 -0.9382789 5 I -2.0408895 0.0000000 0.0000000 -0.0000000 -2.0408895 0.0000000 0.0000000 0.0000000 -0.9132740 6 I -2.0206285 0.0000000 0.0000000 0.0000000 -2.0206285 0.0000000 0.0000000 0.0000000 -0.9093271 7 I -2.0333309 0.0000000 0.0000000 -0.0000000 -2.0333309 0.0000000 0.0000000 0.0000000 -0.9081123 8 I -2.0314621 0.0000000 0.0000000 -0.0000000 -2.0314621 0.0000000 0.0000000 0.0000000 -0.9133378 9 I -2.0498979 0.0000000 0.0000000 0.0000000 -2.0498979 0.0000000 0.0000000 0.0000000 -0.9343347 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000007 (yyy) 0.00000007 (yyy) 0.00000007 (yyy) 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 ... [ 13 13 2 ] 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.58, Number of G-points: 307, Lambda: 0.15 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/42) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 42 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 0.328 ( 0.000 0.000 0.000) 0.000 0.328 ( 0.000 0.000 0.000) 0.000 0.347 ( 0.000 0.000 -0.000) 0.000 0.347 ( 0.000 0.000 -0.000) 0.000 0.595 ( 0.000 0.000 -0.000) 0.000 0.604 ( 0.000 0.000 0.000) 0.000 1.887 ( 0.000 0.000 0.000) 0.000 1.887 ( 0.000 0.000 0.000) 0.000 1.932 ( 0.000 0.000 -0.000) 0.000 1.932 ( 0.000 0.000 -0.000) 0.000 1.954 ( 0.000 0.000 -0.000) 0.000 1.954 ( 0.000 0.000 -0.000) 0.000 2.272 ( 0.000 0.000 0.000) 0.000 2.272 ( 0.000 0.000 0.000) 0.000 2.276 ( 0.000 0.000 0.000) 0.000 2.276 ( 0.000 0.000 0.000) 0.000 2.308 ( 0.000 0.000 -0.000) 0.000 2.308 ( 0.000 0.000 -0.000) 0.000 2.687 ( 0.000 0.000 0.000) 0.000 2.690 ( 0.000 0.000 0.000) 0.000 2.772 ( 0.000 0.000 0.000) 0.000 3.401 ( 0.000 0.000 0.000) 0.000 3.700 ( 0.000 0.000 0.000) 0.000 3.701 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/42) ======================= q-point: ( 0.08 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.164 ( 7.313 4.222 0.969) 8.500 0.225 ( 9.537 5.506 1.679) 11.140 0.393 ( 5.421 3.130 -3.890) 7.370 0.408 ( 16.946 9.784 -1.242) 19.607 0.434 ( 5.828 3.365 2.645) 7.231 0.457 ( 7.036 4.062 6.604) 10.470 0.485 ( 6.077 3.509 -6.803) 9.773 0.633 ( 4.683 2.704 3.708) 6.557 0.644 ( 5.749 3.319 -2.621) 7.137 1.894 ( 0.649 0.375 0.119) 0.759 1.941 ( 0.647 0.373 -0.618) 0.969 1.961 ( 0.584 0.337 0.525) 0.854 2.223 ( 18.810 10.860 -1.438) 21.768 2.239 ( 17.482 10.093 1.115) 20.217 2.253 ( -1.576 -0.910 -0.950) 2.053 2.259 ( -1.365 -0.788 1.026) 1.881 2.274 ( 1.001 0.578 0.605) 1.304 2.278 ( 0.907 0.524 -0.098) 1.052 2.290 ( -0.872 -0.503 -0.051) 1.008 2.292 ( -1.241 -0.717 -0.100) 1.437 2.672 ( -1.172 -0.677 -0.157) 1.362 2.674 ( -1.228 -0.709 0.154) 1.426 2.746 ( -2.182 -1.260 -0.007) 2.520 3.299 ( -0.781 -0.451 -6.083) 6.150 3.349 ( -4.375 -2.526 6.045) 7.878 3.620 ( -6.700 -3.868 -0.918) 7.791 3.623 ( -6.272 -3.621 0.774) 7.283 ======================= Grid point 2 (3/42) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.331 ( 7.294 4.211 0.662) 8.449 0.423 ( 7.712 4.453 1.219) 8.988 0.525 ( -0.220 -0.127 -5.911) 5.917 0.528 ( 5.763 3.327 -3.176) 7.374 0.530 ( 1.131 0.653 5.164) 5.327 0.564 ( 5.211 3.009 2.060) 6.360 0.768 ( 13.794 7.964 -0.774) 15.946 0.820 ( 9.509 5.490 1.205) 11.046 0.850 ( 9.437 5.449 -0.185) 10.899 1.907 ( 0.030 0.017 0.425) 0.427 1.947 ( -0.389 -0.224 -1.512) 1.577 1.976 ( 0.454 0.262 1.188) 1.298 2.212 ( -1.549 -0.895 -1.265) 2.191 2.229 ( -3.040 -1.755 -0.065) 3.511 2.230 ( -0.787 -0.454 1.483) 1.740 2.240 ( -2.416 -1.395 0.487) 2.832 2.250 ( -2.431 -1.404 -0.482) 2.849 2.263 ( -0.954 -0.551 -0.319) 1.147 2.587 ( 5.233 3.021 -2.919) 6.711 2.601 ( 4.244 2.450 2.716) 5.603 2.679 ( 4.345 2.509 -3.751) 6.265 2.683 ( -0.500 -0.289 -0.688) 0.898 2.691 ( 8.028 4.635 5.361) 10.709 3.198 ( -7.654 -4.419 -7.077) 11.323 3.267 ( -4.002 -2.310 7.191) 8.547 3.431 ( -9.003 -5.198 -2.122) 10.610 3.463 ( -6.278 -3.625 1.277) 7.361 ======================= Grid point 3 (4/42) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.462 ( 1.944 1.123 0.688) 2.348 0.515 ( 0.765 0.441 -0.685) 1.117 0.561 ( 2.948 1.702 -0.167) 3.408 0.571 ( 5.318 3.070 1.303) 6.277 0.634 ( 3.521 2.033 -2.182) 4.615 0.670 ( 4.260 2.459 0.954) 5.010 1.000 ( 6.778 3.913 -0.668) 7.855 1.005 ( 5.632 3.252 1.017) 6.582 1.047 ( 8.759 5.057 -0.080) 10.115 1.886 ( -1.555 -0.898 0.745) 1.943 1.920 ( -1.860 -1.074 -1.589) 2.672 1.963 ( -1.643 -0.949 0.889) 2.096 2.150 ( -3.760 -2.171 0.016) 4.342 2.169 ( -3.483 -2.011 0.281) 4.031 2.179 ( -3.441 -1.986 -0.458) 3.999 2.203 ( 0.722 0.417 -0.669) 1.069 2.233 ( 0.970 0.560 1.237) 1.669 2.259 ( 0.475 0.274 -0.614) 0.823 2.589 ( -0.983 -0.568 0.054) 1.137 2.603 ( -0.656 -0.379 0.141) 0.770 2.623 ( -1.436 -0.829 -0.090) 1.661 2.796 ( -1.186 -0.685 -2.301) 2.678 2.839 ( -0.673 -0.389 4.126) 4.198 3.009 ( -9.439 -5.450 -5.214) 12.082 3.145 ( -2.726 -1.574 6.583) 7.297 3.270 ( -3.805 -2.197 -4.884) 6.569 3.387 ( -1.220 -0.704 1.658) 2.176 ======================= Grid point 4 (5/42) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.453 ( 0.268 0.155 6.173) 6.181 0.553 ( 2.278 1.315 -6.543) 7.052 0.667 ( 3.345 1.932 1.849) 4.283 0.668 ( 4.821 2.784 1.191) 5.693 0.702 ( 2.801 1.617 -2.062) 3.836 0.758 ( 3.521 2.033 0.281) 4.076 1.100 ( 3.624 2.092 1.586) 4.475 1.127 ( 4.906 2.833 -2.228) 6.087 1.207 ( 5.317 3.070 0.705) 6.180 1.855 ( -0.970 -0.560 0.932) 1.457 1.877 ( -1.575 -0.909 -1.445) 2.323 1.925 ( -1.191 -0.688 0.493) 1.461 2.062 ( -4.084 -2.358 -0.030) 4.716 2.091 ( -3.505 -2.024 0.016) 4.048 2.103 ( -3.346 -1.932 -0.120) 3.866 2.223 ( 0.524 0.303 -0.165) 0.627 2.258 ( 0.852 0.492 0.911) 1.341 2.272 ( 0.359 0.207 -0.728) 0.838 2.561 ( -1.370 -0.791 -0.346) 1.619 2.591 ( -0.682 -0.394 1.219) 1.451 2.614 ( -0.304 -0.176 -0.731) 0.811 2.663 ( -7.441 -4.296 -0.950) 8.644 2.704 ( -7.942 -4.585 2.821) 9.594 2.766 ( -10.453 -6.035 -2.391) 12.305 3.168 ( 2.033 1.174 3.988) 4.628 3.248 ( 0.278 0.160 -4.923) 4.934 3.368 ( -0.905 -0.522 1.406) 1.752 ======================= Grid point 5 (6/42) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.521 ( 5.312 3.067 5.354) 8.142 0.618 ( 3.012 1.739 -6.690) 7.540 0.736 ( 2.969 1.714 1.334) 3.678 0.748 ( 1.976 1.141 1.764) 2.884 0.764 ( 2.588 1.494 -1.590) 3.385 0.827 ( 2.429 1.402 0.279) 2.819 1.190 ( 4.376 2.526 1.295) 5.216 1.231 ( 4.237 2.446 -2.242) 5.382 1.298 ( 2.853 1.647 0.981) 3.437 1.845 ( -0.004 -0.002 0.137) 0.137 1.856 ( -0.281 -0.162 -0.205) 0.384 1.918 ( 0.399 0.231 0.073) 0.467 1.964 ( -4.325 -2.497 -0.149) 4.996 2.001 ( -4.450 -2.569 0.355) 5.151 2.016 ( -4.430 -2.558 -0.265) 5.122 2.219 ( -0.730 -0.422 -0.094) 0.849 2.264 ( -0.341 -0.197 0.633) 0.746 2.270 ( -0.427 -0.247 -0.544) 0.734 2.484 ( -5.101 -2.945 -0.420) 5.905 2.501 ( -5.458 -3.151 1.072) 6.393 2.525 ( -4.941 -2.852 -0.462) 5.724 2.594 ( -1.144 -0.660 -0.174) 1.332 2.621 ( -0.508 -0.293 0.256) 0.640 2.636 ( -0.501 -0.289 -0.300) 0.651 3.198 ( 0.720 0.416 2.671) 2.798 3.246 ( -0.521 -0.301 -3.622) 3.671 3.340 ( -1.510 -0.872 1.017) 2.019 ======================= Grid point 6 (7/42) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.643 ( 3.982 2.299 0.715) 4.653 0.672 ( 1.771 1.023 -1.794) 2.721 0.755 ( -1.106 -0.639 0.976) 1.608 0.794 ( 1.749 1.010 0.707) 2.140 0.809 ( 1.060 0.612 -0.914) 1.528 0.863 ( 0.718 0.414 0.205) 0.854 1.282 ( 3.023 1.745 0.671) 3.554 1.306 ( 1.995 1.152 -1.243) 2.617 1.334 ( 0.178 0.103 0.617) 0.651 1.849 ( 0.226 0.130 -0.398) 0.476 1.857 ( 0.116 0.067 0.429) 0.449 1.887 ( -1.974 -1.140 -0.149) 2.284 1.912 ( -2.605 -1.504 0.155) 3.012 1.928 ( -2.442 -1.410 -0.037) 2.820 1.930 ( 0.439 0.254 -0.020) 0.508 2.201 ( -0.566 -0.327 -0.044) 0.655 2.251 ( -0.451 -0.260 0.136) 0.538 2.259 ( -0.400 -0.231 -0.101) 0.473 2.405 ( -1.705 -0.985 -0.218) 1.981 2.415 ( -1.879 -1.085 0.564) 2.241 2.450 ( -1.520 -0.877 -0.280) 1.777 2.588 ( 0.087 0.050 -0.356) 0.370 2.618 ( -0.083 -0.048 0.058) 0.111 2.634 ( 0.005 0.003 0.226) 0.226 3.211 ( 0.513 0.296 1.865) 1.957 3.230 ( -0.705 -0.407 -2.241) 2.385 3.308 ( -0.954 -0.551 0.387) 1.167 ======================= Grid point 14 (8/42) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.272 ( 3.898 6.752 0.000) 7.797 0.410 ( 5.502 9.529 0.000) 11.003 0.505 ( 2.267 3.927 0.000) 4.534 0.517 ( 0.562 0.974 0.000) 1.125 0.524 ( 4.171 7.224 0.000) 8.342 0.546 ( 3.304 5.723 0.000) 6.608 0.678 ( 8.076 13.988 0.000) 16.152 0.759 ( 4.841 8.385 0.000) 9.682 0.765 ( 5.008 8.675 0.000) 10.017 1.913 ( 0.760 1.316 0.000) 1.520 1.959 ( 0.669 1.158 0.000) 1.338 1.971 ( 0.443 0.766 0.000) 0.885 2.225 ( -1.058 -1.832 0.000) 2.115 2.228 ( -1.078 -1.868 0.000) 2.156 2.248 ( -1.421 -2.461 0.000) 2.842 2.254 ( -0.986 -1.707 0.000) 1.971 2.262 ( -1.070 -1.854 0.000) 2.140 2.270 ( -0.861 -1.491 0.000) 1.722 2.532 ( 7.966 13.798 0.000) 15.933 2.540 ( 7.693 13.325 0.000) 15.386 2.660 ( 0.767 1.329 0.000) 1.535 2.661 ( 0.742 1.285 0.000) 1.484 2.700 ( -1.870 -3.238 0.000) 3.739 3.250 ( -4.127 -7.148 0.000) 8.253 3.291 ( -0.515 -0.891 0.000) 1.029 3.487 ( -5.317 -9.209 0.000) 10.633 3.488 ( -5.295 -9.171 0.000) 10.590 ======================= Grid point 15 (9/42) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.411 ( 5.232 4.421 -0.200) 6.853 0.499 ( 0.171 -0.481 -1.893) 1.961 0.523 ( 1.072 -0.626 1.927) 2.292 0.597 ( 2.031 9.599 0.901) 9.853 0.660 ( 1.660 8.260 -2.301) 8.734 0.678 ( 2.656 7.667 1.318) 8.220 0.917 ( 6.671 4.859 -0.089) 8.253 0.925 ( 7.113 5.264 0.240) 8.852 0.948 ( 8.311 5.786 0.043) 10.126 1.931 ( -1.395 2.400 0.393) 2.804 1.964 ( -1.862 1.757 -1.334) 2.887 1.988 ( -0.162 0.906 1.020) 1.374 2.177 ( -2.221 -2.957 -0.001) 3.698 2.193 ( -1.064 -1.839 -0.353) 2.154 2.201 ( -2.314 -2.213 0.130) 3.205 2.206 ( -1.217 -1.876 -0.540) 2.301 2.220 ( -0.434 -1.174 0.863) 1.520 2.246 ( 0.155 -1.124 -0.252) 1.162 2.593 ( -0.660 -1.024 -0.347) 1.267 2.599 ( -0.198 -1.266 0.444) 1.356 2.631 ( -1.977 -3.220 -0.057) 3.779 2.809 ( 3.088 7.089 -4.053) 8.730 2.831 ( 5.261 6.945 4.890) 9.991 3.084 ( -5.903 -7.259 -2.306) 9.636 3.215 ( -8.623 -2.898 2.614) 9.465 3.292 ( -5.305 -8.443 -2.234) 10.219 3.343 ( 1.094 -7.933 1.074) 8.079 ======================= Grid point 16 (10/42) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.459 ( -1.189 -0.066 2.892) 3.128 0.520 ( 2.251 0.473 -3.592) 4.265 0.583 ( 5.884 1.179 0.735) 6.046 0.709 ( 0.130 9.218 1.451) 9.332 0.743 ( -0.645 7.875 -1.600) 8.062 0.780 ( 0.655 7.525 0.257) 7.558 1.049 ( 4.865 1.870 1.491) 5.421 1.060 ( 5.914 2.590 -1.601) 6.652 1.132 ( 6.845 3.585 0.263) 7.732 1.912 ( -3.465 3.079 0.655) 4.681 1.934 ( -3.533 2.023 -1.154) 4.231 1.968 ( -2.894 1.045 0.455) 3.110 2.104 ( -3.100 -2.088 0.088) 3.738 2.127 ( -3.006 -2.223 -0.008) 3.738 2.139 ( -2.808 -1.675 -0.069) 3.270 2.205 ( 1.418 -0.173 -0.639) 1.565 2.234 ( 1.873 -0.453 1.093) 2.215 2.253 ( 1.320 -0.847 -0.556) 1.664 2.564 ( -0.922 -1.964 -0.107) 2.173 2.582 ( 0.196 -1.552 0.635) 1.688 2.598 ( 0.717 -1.711 -0.431) 1.904 2.774 ( -8.785 -1.637 -1.133) 9.007 2.819 ( -8.453 -1.426 3.332) 9.197 2.906 ( -10.917 -4.767 -3.854) 12.521 3.110 ( 2.961 -2.449 5.366) 6.600 3.200 ( 2.132 -4.886 -5.007) 7.314 3.322 ( 2.610 -5.446 1.272) 6.171 ======================= Grid point 17 (11/42) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.473 ( 2.579 2.310 5.553) 6.544 0.575 ( 3.325 1.375 -6.610) 7.526 0.688 ( 4.792 0.421 1.542) 5.052 0.786 ( -1.139 8.778 1.407) 8.963 0.808 ( -0.859 8.013 -1.444) 8.187 0.852 ( -0.209 6.500 0.157) 6.505 1.145 ( 4.157 2.742 1.580) 5.225 1.178 ( 4.969 2.619 -2.319) 6.077 1.251 ( 4.398 1.518 0.799) 4.720 1.888 ( -2.658 2.947 0.549) 4.006 1.898 ( -2.765 2.276 -0.686) 3.647 1.932 ( -1.891 0.267 0.067) 1.911 2.033 ( -4.191 0.458 0.026) 4.216 2.052 ( -3.853 -1.531 0.056) 4.147 2.069 ( -4.024 -1.242 -0.085) 4.212 2.220 ( 0.316 -0.347 -0.167) 0.498 2.255 ( 0.959 -0.654 0.831) 1.428 2.262 ( 0.677 -1.038 -0.674) 1.411 2.533 ( -1.486 -1.745 -0.546) 2.356 2.561 ( -2.672 -2.215 1.402) 3.743 2.582 ( -3.131 -2.540 -0.691) 4.091 2.611 ( -3.918 -1.580 -0.388) 4.242 2.649 ( -6.048 -1.354 1.878) 6.476 2.682 ( -8.533 -2.647 -1.793) 9.112 3.149 ( 3.468 -2.935 3.138) 5.522 3.206 ( 2.643 -4.123 -4.051) 6.355 3.307 ( 1.702 -5.316 1.107) 5.690 ======================= Grid point 18 (12/42) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.579 ( 5.521 2.901 2.674) 6.786 0.643 ( 2.556 1.161 -4.154) 5.014 0.739 ( 1.519 -1.370 1.478) 2.524 0.849 ( -1.552 8.339 1.012) 8.543 0.864 ( -1.799 7.528 -1.146) 7.825 0.900 ( -1.260 5.548 0.145) 5.691 1.242 ( 3.832 2.568 0.894) 4.699 1.274 ( 3.341 1.788 -1.588) 4.108 1.315 ( 2.136 0.142 0.742) 2.266 1.880 ( -1.010 1.617 -0.020) 1.907 1.888 ( -1.075 1.737 -0.058) 2.044 1.900 ( -1.494 -0.766 0.005) 1.679 1.962 ( -4.948 0.149 -0.008) 4.950 1.976 ( -4.981 0.060 0.106) 4.983 1.990 ( -3.917 3.617 -0.029) 5.331 2.212 ( -0.711 -0.077 -0.086) 0.720 2.252 ( -0.162 -0.822 0.244) 0.872 2.256 ( 0.023 -1.037 -0.188) 1.055 2.448 ( -4.319 -0.552 -0.375) 4.370 2.462 ( -4.676 -0.702 0.990) 4.831 2.490 ( -3.856 -0.569 -0.479) 3.927 2.579 ( 0.148 -0.979 -0.448) 1.086 2.609 ( 0.314 -1.022 0.207) 1.089 2.623 ( 0.326 -1.005 0.087) 1.060 3.169 ( 2.682 -3.201 2.127) 4.687 3.201 ( 1.611 -4.029 -2.789) 5.158 3.278 ( 1.164 -5.261 0.695) 5.433 ======================= Grid point 19 (13/42) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.663 ( 0.065 -0.113 -0.000) 0.130 0.672 ( 0.499 -0.865 0.000) 0.998 0.727 ( 0.877 -1.519 0.000) 1.754 0.883 ( -3.927 6.801 -0.000) 7.853 0.886 ( -3.715 6.434 0.000) 7.429 0.915 ( -2.637 4.568 -0.000) 5.275 1.299 ( 0.146 -0.253 0.000) 0.293 1.317 ( 0.148 -0.257 -0.000) 0.297 1.325 ( 0.373 -0.647 0.000) 0.747 1.870 ( 0.307 -0.532 -0.000) 0.614 1.885 ( 0.620 -1.074 0.000) 1.240 1.903 ( 0.473 -0.820 -0.000) 0.946 1.911 ( -3.190 5.526 0.000) 6.381 1.917 ( -3.146 5.449 -0.000) 6.292 1.982 ( -2.642 4.576 0.000) 5.284 2.204 ( -0.352 0.610 -0.000) 0.705 2.245 ( 0.225 -0.390 -0.000) 0.450 2.250 ( 0.364 -0.630 0.000) 0.727 2.409 ( -0.848 1.469 0.000) 1.697 2.418 ( -0.871 1.509 -0.000) 1.743 2.454 ( -0.777 1.345 -0.000) 1.554 2.579 ( 0.520 -0.901 -0.000) 1.040 2.607 ( 0.590 -1.022 0.000) 1.180 2.622 ( 0.580 -1.005 -0.000) 1.160 3.180 ( 1.998 -3.460 -0.000) 3.996 3.190 ( 1.975 -3.421 -0.000) 3.950 3.258 ( 2.562 -4.438 -0.000) 5.124 ======================= Grid point 29 (14/42) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.486 ( -0.048 -0.084 0.000) 0.097 0.502 ( 1.352 2.343 0.000) 2.705 0.518 ( 1.400 2.425 0.000) 2.800 0.770 ( 4.257 7.372 0.000) 8.513 0.812 ( 3.695 6.400 0.000) 7.390 0.827 ( 3.906 6.765 0.000) 7.811 0.996 ( 2.652 4.594 0.000) 5.304 1.010 ( 2.573 4.457 0.000) 5.147 1.081 ( 3.762 6.517 0.000) 7.525 1.978 ( 0.755 1.308 0.000) 1.511 1.988 ( -0.174 -0.301 0.000) 0.348 1.999 ( 0.021 0.036 0.000) 0.042 2.120 ( -1.361 -2.358 0.000) 2.723 2.143 ( -1.624 -2.813 0.000) 3.248 2.156 ( -0.858 -1.487 0.000) 1.717 2.197 ( 0.036 0.063 0.000) 0.073 2.201 ( -0.119 -0.206 0.000) 0.238 2.229 ( -0.199 -0.344 0.000) 0.397 2.566 ( -1.675 -2.901 0.000) 3.349 2.571 ( -0.746 -1.292 0.000) 1.492 2.573 ( -0.636 -1.102 0.000) 1.272 2.906 ( 0.087 0.150 0.000) 0.174 2.910 ( -0.008 -0.014 0.000) 0.017 2.918 ( -4.736 -8.203 0.000) 9.472 3.122 ( -4.005 -6.937 0.000) 8.011 3.123 ( -3.931 -6.809 0.000) 7.863 3.224 ( -1.494 -2.587 0.000) 2.987 ======================= Grid point 30 (15/42) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.484 ( -0.379 2.549 3.162) 4.079 0.542 ( 2.410 2.080 -4.099) 5.190 0.609 ( 4.620 1.339 1.062) 4.926 0.891 ( 0.913 8.595 0.359) 8.651 0.909 ( -0.185 8.435 -0.691) 8.465 0.929 ( 0.241 7.077 0.269) 7.086 1.102 ( 4.285 3.780 1.996) 6.053 1.116 ( 5.275 3.359 -2.110) 6.600 1.193 ( 4.035 2.680 0.246) 4.850 1.950 ( -2.548 -1.337 0.106) 2.879 1.961 ( -2.521 -0.107 -0.187) 2.530 1.979 ( -2.671 -0.170 0.055) 2.677 2.088 ( -1.590 0.494 -0.095) 1.667 2.105 ( -1.935 1.104 -0.094) 2.230 2.126 ( -2.580 0.460 0.331) 2.642 2.203 ( 1.045 0.122 -0.759) 1.297 2.220 ( 1.921 -0.803 0.781) 2.224 2.233 ( 1.092 -0.873 -0.188) 1.410 2.523 ( -0.400 -1.745 -0.274) 1.811 2.550 ( -0.066 -1.597 1.210) 2.005 2.561 ( 0.619 -1.822 -0.893) 2.122 2.729 ( -8.926 -2.991 -0.358) 9.420 2.787 ( -9.376 -2.013 2.838) 10.001 2.820 ( -8.437 -4.251 -3.100) 9.944 3.050 ( 2.333 -3.283 4.055) 5.715 3.090 ( 3.964 -5.746 -4.045) 8.068 3.196 ( 2.715 -6.759 0.595) 7.309 ======================= Grid point 31 (16/42) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.538 ( 2.812 3.717 2.538) 5.307 0.615 ( 2.800 2.632 -3.847) 5.437 0.683 ( 3.435 -0.852 1.466) 3.831 0.972 ( -1.567 9.193 0.304) 9.331 0.980 ( -1.572 8.483 -0.538) 8.644 0.990 ( -1.027 7.533 0.158) 7.604 1.207 ( 3.699 3.279 1.126) 5.070 1.231 ( 4.164 2.559 -1.594) 5.141 1.270 ( 3.223 0.427 0.541) 3.296 1.895 ( -1.644 -2.043 -0.059) 2.623 1.916 ( -1.972 -1.291 0.104) 2.359 1.922 ( -2.169 -1.833 -0.053) 2.840 2.064 ( -3.569 3.369 -0.141) 4.910 2.077 ( -4.206 2.451 0.094) 4.869 2.094 ( -2.861 4.034 0.069) 4.946 2.217 ( 0.255 0.129 -0.153) 0.325 2.236 ( 0.975 -1.321 0.130) 1.647 2.239 ( 0.624 -1.056 -0.012) 1.227 2.510 ( -0.834 -0.351 -0.489) 1.029 2.526 ( -1.816 -1.114 1.135) 2.414 2.544 ( -2.050 -1.032 -0.508) 2.351 2.586 ( -3.971 -0.774 -0.346) 4.060 2.636 ( -5.811 0.067 1.769) 6.075 2.653 ( -6.624 -0.389 -1.655) 6.839 3.065 ( 3.936 -5.243 2.312) 6.952 3.100 ( 3.841 -6.308 -2.895) 7.933 3.172 ( 2.952 -8.110 0.693) 8.658 ======================= Grid point 32 (17/42) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.625 ( 3.403 1.404 0.526) 3.719 0.670 ( 1.022 1.470 -1.369) 2.254 0.700 ( 0.395 -1.954 0.865) 2.174 1.016 ( -2.594 7.045 0.206) 7.510 1.023 ( -3.126 8.063 -0.177) 8.650 1.032 ( -3.450 8.196 -0.063) 8.893 1.279 ( 2.054 0.864 0.140) 2.233 1.296 ( 1.788 0.127 -0.400) 1.836 1.310 ( 1.295 -0.236 0.267) 1.343 1.859 ( -0.262 -1.484 -0.075) 1.508 1.869 ( -0.752 -2.185 0.110) 2.314 1.885 ( -0.192 -1.902 -0.024) 1.912 2.038 ( -4.089 5.157 -0.326) 6.590 2.045 ( -3.971 5.659 0.327) 6.921 2.088 ( -3.230 5.331 0.020) 6.233 2.218 ( -0.343 0.590 -0.010) 0.682 2.233 ( 0.491 -1.228 -0.125) 1.328 2.236 ( 0.223 -0.859 0.103) 0.894 2.465 ( -3.391 2.231 -0.200) 4.064 2.477 ( -3.428 2.204 0.577) 4.116 2.505 ( -2.705 2.125 -0.273) 3.450 2.560 ( 0.303 -0.772 -0.358) 0.903 2.590 ( -0.338 -0.414 0.111) 0.546 2.607 ( -0.326 -0.180 0.138) 0.397 3.080 ( 4.100 -5.672 1.098) 7.084 3.098 ( 3.292 -6.211 -1.431) 7.173 3.142 ( 3.376 -8.275 0.345) 8.944 ======================= Grid point 43 (18/42) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.542 ( 1.731 2.999 0.000) 3.463 0.606 ( 2.310 4.001 0.000) 4.620 0.635 ( 0.876 1.518 0.000) 1.753 1.035 ( 3.359 5.817 0.000) 6.717 1.044 ( 2.792 4.837 0.000) 5.585 1.054 ( 3.428 5.937 0.000) 6.856 1.192 ( 2.536 4.392 0.000) 5.072 1.194 ( 2.754 4.771 0.000) 5.509 1.236 ( 1.012 1.752 0.000) 2.024 1.899 ( -1.628 -2.820 0.000) 3.256 1.927 ( -1.612 -2.792 0.000) 3.224 1.933 ( -2.071 -3.587 0.000) 4.141 2.117 ( 0.830 1.438 0.000) 1.661 2.135 ( 0.202 0.349 0.000) 0.403 2.144 ( 0.847 1.467 0.000) 1.694 2.210 ( -0.005 -0.009 0.000) 0.010 2.212 ( 0.369 0.639 0.000) 0.738 2.223 ( -0.169 -0.293 0.000) 0.338 2.507 ( 0.034 0.059 0.000) 0.068 2.525 ( -0.460 -0.797 0.000) 0.921 2.537 ( -0.269 -0.465 0.000) 0.537 2.660 ( -2.219 -3.844 0.000) 4.438 2.730 ( -2.334 -4.043 0.000) 4.668 2.737 ( -2.168 -3.756 0.000) 4.337 2.996 ( -1.641 -2.842 0.000) 3.281 2.998 ( -1.446 -2.505 0.000) 2.892 3.078 ( -2.877 -4.983 0.000) 5.754 ======================= Grid point 44 (19/42) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 170 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.598 ( 1.772 2.089 0.324) 2.759 0.655 ( 1.155 -0.064 -0.351) 1.209 0.685 ( 1.756 2.625 0.071) 3.159 1.106 ( -0.339 5.115 0.076) 5.127 1.130 ( -0.194 6.455 -0.175) 6.461 1.144 ( -0.209 6.836 0.087) 6.840 1.256 ( 1.305 0.608 -0.054) 1.441 1.266 ( 1.991 1.029 -0.068) 2.242 1.279 ( 2.550 1.174 0.109) 2.809 1.853 ( -0.959 -1.953 -0.019) 2.176 1.864 ( -1.564 -2.840 0.062) 3.243 1.875 ( -1.187 -2.531 -0.024) 2.795 2.127 ( -1.856 2.256 -0.361) 2.944 2.133 ( -1.286 2.708 0.303) 3.013 2.167 ( -0.843 3.030 0.081) 3.146 2.210 ( 0.592 -1.272 0.039) 1.403 2.218 ( 0.426 -0.746 -0.002) 0.859 2.227 ( 0.290 0.454 -0.072) 0.543 2.515 ( -0.279 0.409 -0.186) 0.529 2.519 ( 0.264 0.400 0.378) 0.611 2.538 ( 0.303 0.165 -0.162) 0.381 2.595 ( -3.832 2.030 -0.164) 4.339 2.658 ( -4.460 2.326 0.923) 5.114 2.672 ( -4.309 2.383 -0.825) 4.993 2.955 ( 2.258 -5.768 0.964) 6.269 2.973 ( 2.799 -6.335 -1.294) 7.046 3.002 ( 1.827 -8.893 0.368) 9.086 ======================= Grid point 45 (20/42) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.628 ( 0.506 -0.876 0.000) 1.012 0.662 ( 0.338 -0.586 0.000) 0.676 0.714 ( -0.938 1.625 -0.000) 1.876 1.129 ( -2.614 4.528 0.000) 5.228 1.158 ( -3.057 5.294 0.000) 6.113 1.177 ( -3.356 5.812 -0.000) 6.711 1.269 ( 0.845 -1.464 0.000) 1.691 1.289 ( 0.312 -0.541 -0.000) 0.625 1.306 ( 0.335 -0.581 0.000) 0.671 1.832 ( 0.651 -1.128 0.000) 1.302 1.837 ( 0.584 -1.011 -0.000) 1.167 1.851 ( 0.733 -1.270 -0.000) 1.466 2.121 ( -1.788 3.097 0.000) 3.577 2.135 ( -1.814 3.142 -0.000) 3.628 2.175 ( -1.722 2.982 0.000) 3.444 2.208 ( 0.783 -1.356 -0.000) 1.566 2.219 ( 0.393 -0.680 0.000) 0.786 2.231 ( -0.280 0.485 -0.000) 0.560 2.512 ( -0.756 1.309 0.000) 1.511 2.523 ( -0.704 1.219 -0.000) 1.407 2.544 ( -0.216 0.374 -0.000) 0.432 2.574 ( -2.544 4.406 0.000) 5.087 2.631 ( -2.929 5.074 0.000) 5.858 2.648 ( -2.798 4.846 0.000) 5.595 2.952 ( 4.207 -7.286 -0.000) 8.413 2.966 ( 4.300 -7.447 0.000) 8.599 2.967 ( 4.858 -8.414 -0.000) 9.715 ======================= Grid point 58 (21/42) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 98 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.616 ( 0.053 0.091 0.000) 0.105 0.658 ( 0.261 0.451 0.000) 0.521 0.732 ( 0.821 1.422 0.000) 1.642 1.181 ( 1.448 2.508 0.000) 2.896 1.219 ( 1.207 2.091 0.000) 2.414 1.241 ( 0.739 1.280 0.000) 1.477 1.246 ( 0.808 1.399 0.000) 1.616 1.279 ( -0.042 -0.073 0.000) 0.084 1.291 ( -0.121 -0.210 0.000) 0.243 1.822 ( -0.592 -1.026 0.000) 1.185 1.826 ( -0.566 -0.980 0.000) 1.132 1.838 ( -0.546 -0.945 0.000) 1.091 2.158 ( 0.566 0.980 0.000) 1.132 2.167 ( 0.464 0.804 0.000) 0.929 2.186 ( -0.535 -0.927 0.000) 1.071 2.199 ( -0.052 -0.089 0.000) 0.103 2.214 ( 0.418 0.724 0.000) 0.836 2.236 ( 0.158 0.273 0.000) 0.315 2.521 ( 0.124 0.214 0.000) 0.248 2.529 ( 0.204 0.353 0.000) 0.408 2.546 ( 0.258 0.447 0.000) 0.516 2.663 ( 2.522 4.369 0.000) 5.045 2.720 ( 1.955 3.387 0.000) 3.910 2.737 ( 2.097 3.633 0.000) 4.195 2.836 ( -3.845 -6.660 0.000) 7.690 2.851 ( -2.800 -4.850 0.000) 5.601 2.863 ( -2.637 -4.567 0.000) 5.274 ======================= Grid point 181 (22/42) ======================= q-point: ( 0.00 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 21 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.179 ( 0.000 -0.000 0.000) 0.000 0.179 ( 0.000 0.000 0.000) 0.000 0.192 ( 0.000 0.000 0.000) 0.000 0.192 ( 0.000 -0.000 0.000) 0.000 0.335 ( 0.000 0.000 0.000) 0.000 0.344 ( 0.000 -0.000 -0.000) 0.000 0.399 ( 0.000 -0.000 0.000) 0.000 0.399 ( 0.000 0.000 0.000) 0.000 0.700 ( 0.000 -0.000 0.000) 0.000 1.897 ( 0.000 -0.000 0.000) 0.000 1.897 ( 0.000 0.000 0.000) 0.000 1.908 ( 0.000 0.000 0.000) 0.000 1.908 ( 0.000 -0.000 0.000) 0.000 1.969 ( 0.000 -0.000 0.000) 0.000 1.969 ( 0.000 -0.000 0.000) 0.000 2.256 ( 0.000 -0.000 0.000) 0.000 2.256 ( 0.000 0.000 0.000) 0.000 2.295 ( -0.000 0.000 0.000) 0.000 2.295 ( 0.000 0.000 0.000) 0.000 2.303 ( 0.000 0.000 -0.000) 0.000 2.303 ( 0.000 -0.000 -0.000) 0.000 2.658 ( 0.000 -0.000 0.000) 0.000 2.745 ( 0.000 0.000 -0.000) 0.000 2.746 ( 0.000 0.000 0.000) 0.000 3.688 ( -0.000 0.000 0.000) 0.000 3.692 ( 0.000 0.000 0.000) 0.000 3.706 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 183 (23/42) ======================= q-point: ( 0.08 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.267 ( 7.456 4.304 -4.237) 9.595 0.283 ( 5.059 2.921 -0.515) 5.864 0.297 ( 4.491 2.593 -1.243) 5.332 0.322 ( 7.658 4.422 4.511) 9.927 0.467 ( 5.029 2.903 -0.685) 5.847 0.477 ( 11.774 6.798 6.473) 15.057 0.509 ( 13.057 7.538 -5.703) 16.119 0.574 ( 11.967 6.909 0.784) 13.840 0.681 ( -1.648 -0.952 -0.219) 1.916 1.906 ( 0.599 0.346 -0.603) 0.917 1.915 ( 0.688 0.397 0.761) 1.100 1.978 ( 0.625 0.361 -0.183) 0.745 2.184 ( 16.214 9.361 0.201) 18.724 2.236 ( -1.606 -0.927 0.141) 1.860 2.255 ( -0.315 -0.182 -0.090) 0.375 2.278 ( -1.347 -0.778 -0.408) 1.608 2.279 ( -0.553 -0.319 0.083) 0.644 2.288 ( -1.260 -0.727 0.292) 1.484 2.290 ( -0.523 -0.302 -0.026) 0.605 2.551 ( 24.807 14.322 -23.040) 36.761 2.559 ( 25.465 14.702 23.099) 37.392 2.644 ( -1.074 -0.620 0.011) 1.240 2.730 ( -0.157 -0.090 -0.416) 0.454 2.732 ( -0.041 -0.024 0.494) 0.497 3.579 ( -8.301 -4.793 -2.951) 10.029 3.586 ( -7.560 -4.365 2.753) 9.153 3.631 ( -6.145 -3.548 -0.064) 7.096 ======================= Grid point 185 (24/42) ======================= q-point: ( 0.15 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.391 ( 5.043 2.912 -4.745) 7.512 0.412 ( 5.448 3.146 4.774) 7.898 0.444 ( 7.318 4.225 -2.664) 8.860 0.488 ( 6.715 3.877 3.284) 8.420 0.583 ( 4.662 2.692 -0.722) 5.431 0.619 ( -3.825 -2.208 0.072) 4.417 0.770 ( 12.466 7.197 1.379) 14.460 0.801 ( 11.417 6.591 -1.707) 13.293 0.836 ( 10.170 5.871 0.069) 11.743 1.912 ( -0.301 -0.174 -0.850) 0.919 1.932 ( 0.405 0.234 1.418) 1.493 1.987 ( -0.096 -0.056 -0.674) 0.683 2.200 ( -1.013 -0.585 0.535) 1.286 2.226 ( -2.157 -1.245 -0.161) 2.496 2.241 ( -2.624 -1.515 0.160) 3.034 2.244 ( -1.367 -0.789 -1.056) 1.899 2.251 ( -2.729 -1.576 0.059) 3.152 2.259 ( -0.826 -0.477 0.627) 1.142 2.555 ( 11.433 6.601 0.476) 13.210 2.631 ( 2.334 1.347 0.090) 2.696 2.656 ( -2.404 -1.388 -2.125) 3.496 2.659 ( -2.218 -1.281 1.868) 3.170 2.905 ( 5.263 3.039 -13.142) 14.479 2.941 ( 7.219 4.168 14.037) 16.325 3.361 ( -10.354 -5.978 -3.935) 12.587 3.406 ( -7.589 -4.381 3.562) 9.459 3.466 ( -7.049 -4.069 -0.816) 8.180 ======================= Grid point 187 (25/42) ======================= q-point: ( 0.23 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.482 ( -1.968 -1.136 -1.383) 2.660 0.514 ( 2.287 1.321 1.534) 3.054 0.551 ( 4.072 2.351 0.016) 4.701 0.579 ( 4.493 2.594 -1.814) 5.497 0.625 ( 5.343 3.085 2.501) 6.658 0.671 ( 3.248 1.875 -0.762) 3.827 0.987 ( 5.851 3.378 0.034) 6.756 1.021 ( 8.671 5.006 -0.826) 10.047 1.032 ( 7.245 4.183 0.522) 8.382 1.887 ( -1.736 -1.002 -0.778) 2.150 1.920 ( -1.332 -0.769 1.657) 2.260 1.962 ( -2.061 -1.190 -0.924) 2.553 2.162 ( -3.416 -1.972 -0.268) 3.953 2.168 ( -3.626 -2.094 -0.094) 4.188 2.174 ( -3.503 -2.022 0.523) 4.078 2.206 ( 1.295 0.748 0.831) 1.711 2.231 ( 0.303 0.175 -1.394) 1.438 2.259 ( 0.639 0.369 0.609) 0.957 2.595 ( -1.140 -0.658 0.103) 1.321 2.605 ( -1.119 -0.646 -0.683) 1.461 2.613 ( -0.986 -0.570 0.478) 1.235 2.779 ( 0.332 0.192 0.848) 0.931 2.888 ( -5.816 -3.358 -4.727) 8.212 2.977 ( -3.815 -2.203 6.783) 8.088 3.143 ( -6.109 -3.527 -5.553) 8.978 3.293 ( -1.675 -0.967 5.618) 5.942 3.367 ( -1.779 -1.027 -2.939) 3.586 ======================= Grid point 189 (26/42) ======================= q-point: ( 0.31 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.396 ( -1.196 -0.691 -3.374) 3.646 0.595 ( 3.283 1.896 4.039) 5.540 0.657 ( 2.843 1.641 -1.114) 3.467 0.666 ( 4.633 2.675 -1.561) 5.573 0.723 ( 3.179 1.836 1.991) 4.177 0.746 ( 3.613 2.086 -0.946) 4.278 1.086 ( 3.742 2.161 -0.609) 4.364 1.158 ( 4.534 2.618 2.521) 5.811 1.187 ( 5.796 3.346 -1.975) 6.978 1.849 ( -1.299 -0.750 -0.528) 1.591 1.893 ( -0.716 -0.413 1.731) 1.918 1.914 ( -1.714 -0.989 -1.183) 2.305 2.076 ( -4.058 -2.343 -0.021) 4.685 2.082 ( -3.897 -2.250 0.065) 4.501 2.103 ( -3.084 -1.780 0.091) 3.562 2.234 ( 0.760 0.439 0.913) 1.267 2.243 ( 0.363 0.209 -1.293) 1.359 2.276 ( 0.608 0.351 0.361) 0.790 2.561 ( -1.555 -0.898 0.452) 1.852 2.584 ( -1.017 -0.587 -1.053) 1.577 2.618 ( 0.229 0.132 0.460) 0.530 2.660 ( -7.015 -4.050 0.767) 8.136 2.703 ( -8.678 -5.010 -2.309) 10.283 2.774 ( -10.432 -6.023 2.082) 12.224 3.147 ( 2.275 1.313 -2.500) 3.626 3.286 ( 0.069 0.040 4.931) 4.931 3.347 ( -0.668 -0.386 -2.921) 3.021 ======================= Grid point 191 (27/42) ======================= q-point: ( 0.38 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.473 ( 6.821 3.938 -2.781) 8.353 0.660 ( 2.480 1.432 4.769) 5.563 0.722 ( 2.986 1.724 -0.616) 3.502 0.742 ( 1.764 1.019 -2.437) 3.177 0.785 ( 2.491 1.438 1.124) 3.088 0.819 ( 2.516 1.452 -0.531) 2.953 1.180 ( 4.600 2.656 -0.735) 5.362 1.248 ( 3.462 1.999 2.131) 4.530 1.291 ( 3.401 1.964 -1.429) 4.179 1.832 ( -0.264 -0.152 -0.024) 0.306 1.890 ( 0.239 0.138 0.718) 0.769 1.896 ( 0.110 0.064 -0.698) 0.710 1.977 ( -4.520 -2.610 0.555) 5.249 1.987 ( -4.375 -2.526 -0.656) 5.094 2.019 ( -4.445 -2.566 0.160) 5.135 2.235 ( -0.524 -0.302 1.066) 1.225 2.240 ( -0.630 -0.364 -1.166) 1.374 2.279 ( -0.317 -0.183 0.104) 0.381 2.478 ( -5.094 -2.941 0.211) 5.886 2.502 ( -4.881 -2.818 -0.752) 5.686 2.528 ( -5.506 -3.179 0.357) 6.368 2.592 ( -0.490 -0.283 0.119) 0.578 2.620 ( -1.150 -0.664 -0.104) 1.332 2.641 ( -0.483 -0.279 0.197) 0.591 3.180 ( 0.739 0.426 -1.308) 1.561 3.281 ( -0.541 -0.312 3.544) 3.598 3.323 ( -1.476 -0.852 -2.304) 2.866 ======================= Grid point 193 (28/42) ======================= q-point: ( 0.46 0.00 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.630 ( 5.772 3.332 -0.175) 6.667 0.700 ( 0.778 0.449 2.128) 2.310 0.741 ( -1.689 -0.975 -1.850) 2.688 0.782 ( 1.771 1.022 -0.201) 2.055 0.831 ( 1.273 0.735 0.671) 1.616 0.853 ( 0.495 0.286 -0.468) 0.739 1.279 ( 3.526 2.036 -0.535) 4.107 1.308 ( 1.453 0.839 1.130) 2.023 1.334 ( 0.216 0.125 -0.640) 0.687 1.832 ( 0.108 0.062 0.058) 0.137 1.893 ( -2.277 -1.314 0.329) 2.649 1.897 ( 0.297 0.172 -0.393) 0.522 1.904 ( -2.323 -1.341 -0.321) 2.701 1.905 ( 0.348 0.201 0.325) 0.517 1.929 ( -2.530 -1.461 0.023) 2.922 2.221 ( -0.467 -0.270 0.694) 0.879 2.224 ( -0.496 -0.287 -0.700) 0.905 2.268 ( -0.426 -0.246 0.016) 0.492 2.399 ( -1.670 -0.964 0.034) 1.929 2.429 ( -1.381 -0.797 -0.900) 1.831 2.442 ( -1.956 -1.129 0.801) 2.397 2.596 ( 0.477 0.275 0.606) 0.818 2.613 ( 0.035 0.020 -0.340) 0.342 2.631 ( -0.519 -0.300 -0.194) 0.630 3.190 ( 0.289 0.167 -0.460) 0.568 3.270 ( -0.180 -0.104 2.141) 2.151 3.288 ( -1.240 -0.716 -1.691) 2.216 ======================= Grid point 209 (29/42) ======================= q-point: ( 0.08 0.08 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.341 ( 2.111 3.656 0.000) 4.222 0.352 ( 2.895 5.015 -0.000) 5.790 0.458 ( 5.095 8.825 0.000) 10.190 0.464 ( 5.670 9.820 0.000) 11.339 0.569 ( 4.032 6.983 0.000) 8.064 0.641 ( -1.729 -2.995 0.000) 3.458 0.698 ( 6.773 11.730 0.000) 13.545 0.714 ( 6.954 12.045 0.000) 13.908 0.744 ( 5.429 9.403 -0.000) 10.858 1.926 ( 0.721 1.248 0.000) 1.441 1.931 ( 0.707 1.224 -0.000) 1.414 1.989 ( 0.473 0.819 0.000) 0.945 2.205 ( -1.075 -1.862 0.000) 2.150 2.236 ( -0.995 -1.723 0.000) 1.990 2.251 ( -1.027 -1.780 0.000) 2.055 2.255 ( -1.211 -2.097 0.000) 2.421 2.264 ( -0.859 -1.487 0.000) 1.717 2.267 ( -1.111 -1.924 0.000) 2.222 2.470 ( 9.177 15.895 0.000) 18.354 2.628 ( -0.225 -0.390 0.000) 0.450 2.669 ( -1.271 -2.202 -0.000) 2.542 2.671 ( -1.280 -2.217 0.000) 2.559 2.875 ( 5.882 10.187 -0.000) 11.763 2.878 ( 5.973 10.346 0.000) 11.947 3.440 ( -5.062 -8.768 -0.000) 10.125 3.442 ( -5.100 -8.834 0.000) 10.200 3.499 ( -5.336 -9.242 0.000) 10.672 ======================= Grid point 211 (30/42) ======================= q-point: ( 0.15 0.08 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 169 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.434 ( 4.730 2.315 0.137) 5.268 0.459 ( 4.326 3.286 0.597) 5.466 0.552 ( -3.355 -4.100 -0.517) 5.323 0.606 ( 1.803 9.471 -1.511) 9.759 0.640 ( 3.162 9.085 1.780) 9.782 0.695 ( 1.291 7.702 -0.234) 7.813 0.905 ( 6.732 4.417 0.012) 8.052 0.923 ( 8.167 6.261 -0.179) 10.292 0.943 ( 8.032 5.871 -0.033) 9.949 1.936 ( -1.846 2.399 -0.756) 3.120 1.953 ( -0.661 1.894 1.299) 2.390 1.996 ( -1.129 0.814 -0.627) 1.527 2.181 ( 0.302 -1.319 0.254) 1.377 2.190 ( -2.036 -2.325 0.016) 3.090 2.197 ( -2.380 -3.062 -0.595) 3.924 2.201 ( -2.283 -2.422 0.802) 3.424 2.231 ( -0.664 -0.810 -0.844) 1.346 2.243 ( 0.374 -1.098 0.523) 1.272 2.592 ( 0.103 -0.600 0.167) 0.631 2.614 ( -1.800 -2.587 -1.137) 3.350 2.618 ( -1.418 -2.635 0.931) 3.134 2.768 ( 5.776 8.935 0.469) 10.650 2.952 ( -2.903 0.926 -6.915) 7.557 2.992 ( 0.898 0.073 7.425) 7.479 3.233 ( -9.985 -6.458 -2.581) 12.168 3.290 ( -3.149 -6.774 2.947) 8.030 3.336 ( -0.578 -8.707 -1.328) 8.827 ======================= Grid point 213 (31/42) ======================= q-point: ( 0.23 0.08 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 169 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.431 ( -3.740 -1.986 -1.677) 4.554 0.534 ( 4.315 1.152 2.351) 5.047 0.586 ( 5.726 1.703 -0.714) 6.016 0.706 ( -0.510 8.930 -1.054) 9.007 0.757 ( 0.633 8.127 1.715) 8.330 0.774 ( -0.214 7.513 -0.777) 7.556 1.031 ( 5.170 1.669 -0.319) 5.442 1.096 ( 6.004 3.109 1.473) 6.920 1.110 ( 6.724 3.627 -1.306) 7.750 1.911 ( -3.773 3.001 -0.545) 4.852 1.939 ( -2.914 2.266 1.278) 3.906 1.964 ( -3.266 1.072 -0.693) 3.506 2.114 ( -3.104 -2.731 -0.222) 4.141 2.122 ( -3.089 -2.084 0.152) 3.730 2.137 ( -2.757 -1.438 0.070) 3.111 2.208 ( 2.296 -0.404 0.793) 2.463 2.230 ( 0.860 -0.316 -1.153) 1.473 2.254 ( 1.549 -0.776 0.454) 1.791 2.565 ( -0.881 -2.087 0.169) 2.271 2.578 ( -0.212 -2.060 -0.613) 2.159 2.599 ( 1.065 -1.071 0.346) 1.549 2.770 ( -8.526 -1.021 0.822) 8.626 2.823 ( -9.587 -3.283 -2.822) 10.519 2.919 ( -10.928 -3.308 4.062) 12.119 3.080 ( 3.522 -3.394 -4.111) 6.389 3.242 ( 2.294 -4.116 5.025) 6.889 3.299 ( 2.780 -5.504 -2.952) 6.836 ======================= Grid point 215 (32/42) ======================= q-point: ( 0.31 0.08 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 169 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.425 ( 2.779 3.844 -3.164) 5.702 0.609 ( 3.732 0.448 4.481) 5.849 0.688 ( 4.527 0.458 -1.832) 4.905 0.775 ( -1.425 9.018 -0.703) 9.157 0.826 ( -0.830 7.375 0.969) 7.485 0.846 ( -0.091 6.907 -0.397) 6.919 1.131 ( 4.461 2.801 -0.730) 5.318 1.204 ( 4.075 2.167 2.322) 5.167 1.238 ( 5.013 1.957 -1.651) 5.629 1.881 ( -3.023 3.053 -0.170) 4.300 1.916 ( -1.908 1.579 0.680) 2.568 1.925 ( -2.353 1.289 -0.442) 2.719 2.037 ( -4.008 -0.820 -0.032) 4.091 2.047 ( -4.267 -0.727 -0.037) 4.329 2.070 ( -3.888 -1.319 0.073) 4.106 2.232 ( 0.736 -0.532 0.897) 1.277 2.238 ( 0.274 -0.528 -1.061) 1.216 2.268 ( 0.941 -0.976 0.173) 1.367 2.530 ( -1.954 -1.844 0.560) 2.744 2.554 ( -1.918 -1.992 -1.077) 2.967 2.587 ( -3.637 -2.495 0.389) 4.427 2.617 ( -3.923 -1.494 0.685) 4.254 2.636 ( -5.593 -1.887 -1.545) 6.101 2.693 ( -8.911 -2.349 1.130) 9.284 3.128 ( 3.850 -2.927 -1.669) 5.116 3.243 ( 2.206 -4.272 3.928) 6.209 3.290 ( 1.922 -5.152 -2.457) 6.023 ======================= Grid point 217 (33/42) ======================= q-point: ( 0.38 0.08 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 169 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.558 ( 6.203 4.405 -1.488) 7.752 0.668 ( 2.738 -0.148 3.556) 4.491 0.732 ( 0.912 -1.414 -2.074) 2.671 0.837 ( -1.732 8.697 -0.361) 8.875 0.878 ( -1.422 6.863 0.618) 7.036 0.898 ( -1.413 5.952 -0.262) 6.123 1.236 ( 4.131 2.923 -0.596) 5.096 1.280 ( 2.778 1.185 1.384) 3.322 1.314 ( 2.384 0.424 -0.836) 2.562 1.873 ( -1.403 3.084 0.015) 3.389 1.896 ( -1.091 -1.151 0.098) 1.589 1.905 ( -1.001 -0.573 -0.041) 1.154 1.967 ( -4.699 0.735 0.351) 4.769 1.972 ( -4.429 3.207 -0.607) 5.502 1.983 ( -4.936 1.100 0.188) 5.061 2.227 ( -0.377 -0.513 0.860) 1.071 2.230 ( -0.522 -0.352 -0.848) 1.056 2.265 ( 0.057 -1.136 0.016) 1.137 2.443 ( -4.372 -0.495 0.139) 4.402 2.468 ( -3.913 -0.490 -0.936) 4.053 2.487 ( -4.466 -0.812 0.664) 4.587 2.582 ( 0.708 -0.772 0.526) 1.172 2.604 ( 0.096 -1.057 -0.298) 1.103 2.626 ( -0.044 -1.149 -0.078) 1.152 3.152 ( 2.627 -3.055 -0.853) 4.118 3.233 ( 1.867 -4.357 2.666) 5.438 3.262 ( 0.995 -5.075 -1.846) 5.491 ======================= Grid point 219 (34/42) ======================= q-point: ( 0.46 0.08 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.668 ( 0.227 -0.393 0.000) 0.454 0.690 ( 0.419 -0.725 0.000) 0.838 0.703 ( 0.907 -1.571 0.000) 1.814 0.870 ( -4.045 7.006 0.000) 8.089 0.903 ( -3.251 5.632 0.000) 6.503 0.913 ( -3.010 5.213 0.000) 6.020 1.302 ( 0.111 -0.193 -0.000) 0.223 1.309 ( 0.304 -0.527 0.000) 0.608 1.328 ( 0.216 -0.375 0.000) 0.433 1.872 ( 0.406 -0.703 -0.000) 0.812 1.881 ( 0.464 -0.803 0.000) 0.927 1.897 ( -3.334 5.774 0.000) 6.668 1.905 ( 0.515 -0.893 0.000) 1.031 1.952 ( -2.843 4.925 0.000) 5.686 1.960 ( -2.876 4.981 0.000) 5.752 2.219 ( 0.041 -0.071 -0.000) 0.082 2.222 ( -0.161 0.278 0.000) 0.321 2.258 ( 0.433 -0.751 0.000) 0.867 2.404 ( -0.860 1.489 0.000) 1.719 2.435 ( -0.765 1.325 0.000) 1.530 2.444 ( -0.842 1.458 0.000) 1.683 2.589 ( 0.560 -0.970 0.000) 1.120 2.602 ( 0.621 -1.076 0.000) 1.242 2.618 ( 0.490 -0.849 0.000) 0.980 3.160 ( 1.819 -3.150 0.000) 3.637 3.229 ( 2.368 -4.102 -0.000) 4.737 3.239 ( 2.350 -4.069 0.000) 4.699 ======================= Grid point 239 (35/42) ======================= q-point: ( 0.15 0.15 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.460 ( -1.633 -2.829 0.000) 3.267 0.501 ( 1.865 3.230 0.000) 3.729 0.534 ( 2.094 3.627 0.000) 4.188 0.775 ( 4.210 7.292 0.000) 8.420 0.801 ( 3.993 6.916 0.000) 7.986 0.841 ( 3.559 6.165 0.000) 7.119 0.971 ( 2.258 3.912 0.000) 4.517 1.045 ( 3.562 6.170 0.000) 7.124 1.063 ( 3.424 5.930 0.000) 6.847 1.979 ( 0.593 1.028 0.000) 1.187 1.989 ( 0.321 0.556 0.000) 0.642 2.001 ( -0.226 -0.391 0.000) 0.451 2.126 ( -1.954 -3.385 -0.000) 3.908 2.139 ( -1.395 -2.416 0.000) 2.789 2.164 ( -0.707 -1.224 0.000) 1.413 2.170 ( 0.091 0.158 -0.000) 0.183 2.221 ( -0.155 -0.268 0.000) 0.309 2.225 ( -0.132 -0.229 -0.000) 0.265 2.562 ( -1.418 -2.455 0.000) 2.835 2.565 ( -1.293 -2.239 0.000) 2.586 2.579 ( -0.356 -0.616 0.000) 0.712 2.910 ( -3.457 -5.988 -0.000) 6.914 2.915 ( -2.252 -3.900 0.000) 4.503 2.917 ( 1.396 2.418 0.000) 2.792 3.089 ( -6.029 -10.442 -0.000) 12.057 3.185 ( -1.877 -3.251 0.000) 3.754 3.189 ( -1.858 -3.218 -0.000) 3.716 ======================= Grid point 241 (36/42) ======================= q-point: ( 0.23 0.15 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 169 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.444 ( 0.305 3.407 -1.388) 3.692 0.565 ( 2.316 1.895 2.263) 3.752 0.614 ( 4.226 1.235 -1.008) 4.517 0.891 ( 0.214 9.138 -0.322) 9.146 0.911 ( 0.965 7.245 0.662) 7.339 0.929 ( -0.476 7.769 -0.274) 7.788 1.077 ( 5.291 3.441 -0.431) 6.326 1.159 ( 3.960 3.327 1.594) 5.412 1.173 ( 4.518 2.927 -1.298) 5.537 1.952 ( -2.627 -0.735 -0.126) 2.731 1.966 ( -2.528 -0.087 0.299) 2.547 1.976 ( -2.669 -0.977 -0.146) 2.846 2.087 ( -1.389 1.020 0.102) 1.726 2.102 ( -2.042 0.380 0.021) 2.078 2.130 ( -2.969 0.876 -0.273) 3.108 2.198 ( 2.287 -0.424 0.544) 2.388 2.224 ( 0.713 -0.232 -0.548) 0.929 2.234 ( 1.409 -1.052 0.179) 1.767 2.525 ( -0.488 -1.660 0.511) 1.804 2.537 ( 0.202 -1.776 -0.799) 1.958 2.570 ( 0.474 -1.648 0.245) 1.732 2.738 ( -9.215 -2.316 1.155) 9.572 2.756 ( -8.304 -3.585 -1.934) 9.249 2.855 ( -11.090 -3.407 1.619) 11.714 3.006 ( 5.383 -3.757 -1.780) 6.802 3.146 ( 1.942 -5.221 3.430) 6.542 3.172 ( 3.531 -6.783 -2.475) 8.038 ======================= Grid point 243 (37/42) ======================= q-point: ( 0.31 0.15 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 169 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.527 ( 3.250 5.689 -1.993) 6.848 0.620 ( 2.682 0.542 3.188) 4.201 0.686 ( 3.172 -0.382 -1.354) 3.470 0.969 ( -1.828 9.497 -0.143) 9.672 0.981 ( -1.222 7.922 0.306) 8.021 0.993 ( -1.192 7.847 -0.086) 7.938 1.198 ( 4.322 3.743 -0.599) 5.749 1.242 ( 3.068 1.569 1.363) 3.705 1.266 ( 3.802 0.892 -0.836) 3.994 1.903 ( -1.691 -1.714 0.225) 2.418 1.909 ( -2.155 -2.242 -0.243) 3.119 1.922 ( -1.961 -1.421 0.026) 2.422 2.067 ( -3.998 2.809 0.223) 4.891 2.079 ( -2.894 4.087 -0.221) 5.012 2.087 ( -3.771 3.194 -0.022) 4.942 2.219 ( 0.702 -0.595 0.212) 0.944 2.230 ( 0.218 -0.213 -0.156) 0.343 2.243 ( 0.956 -1.531 -0.021) 1.806 2.506 ( -1.235 -0.417 0.308) 1.339 2.526 ( -1.223 -0.624 -0.824) 1.601 2.547 ( -2.038 -1.401 0.383) 2.502 2.595 ( -5.129 -0.398 0.874) 5.218 2.611 ( -4.253 -0.659 -1.181) 4.463 2.671 ( -7.395 -0.025 0.537) 7.415 3.046 ( 4.592 -5.039 -0.998) 6.890 3.130 ( 3.074 -6.892 2.598) 7.981 3.159 ( 3.244 -7.749 -1.713) 8.574 ======================= Grid point 245 (38/42) ======================= q-point: ( 0.38 0.15 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 169 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.632 ( 3.436 2.502 -0.910) 4.347 0.659 ( 1.733 -0.725 1.526) 2.420 0.703 ( -0.424 -0.860 -0.639) 1.152 1.014 ( -2.976 7.850 -0.122) 8.396 1.025 ( -2.874 7.553 0.084) 8.082 1.033 ( -3.282 7.835 0.072) 8.495 1.280 ( 2.142 0.825 -0.160) 2.301 1.292 ( 1.633 0.250 0.325) 1.683 1.313 ( 1.406 -0.269 -0.172) 1.442 1.861 ( -0.617 -1.777 0.126) 1.885 1.866 ( -0.370 -1.923 -0.154) 1.964 1.886 ( -0.200 -1.795 0.016) 1.806 2.030 ( -4.376 5.289 0.084) 6.865 2.065 ( -3.526 5.200 -0.138) 6.284 2.076 ( -3.446 5.683 0.034) 6.647 2.218 ( -0.001 -0.350 0.015) 0.350 2.229 ( -0.195 0.157 0.029) 0.252 2.239 ( 0.527 -1.401 -0.012) 1.497 2.462 ( -3.471 2.388 0.104) 4.214 2.484 ( -2.891 2.093 -0.794) 3.656 2.499 ( -3.125 1.970 0.598) 3.742 2.567 ( 0.621 -0.513 0.557) 0.979 2.582 ( 0.073 -0.824 -0.353) 0.900 2.610 ( -1.061 0.078 -0.105) 1.069 3.069 ( 3.966 -5.158 -0.405) 6.519 3.117 ( 3.756 -7.286 1.268) 8.295 3.134 ( 3.083 -7.709 -0.874) 8.349 ======================= Grid point 267 (39/42) ======================= q-point: ( 0.23 0.23 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.537 ( 2.983 5.166 0.000) 5.965 0.599 ( 0.829 1.436 0.000) 1.658 0.645 ( 1.277 2.213 0.000) 2.555 1.033 ( 2.985 5.170 0.000) 5.969 1.043 ( 3.350 5.803 0.000) 6.700 1.060 ( 3.222 5.580 0.000) 6.443 1.171 ( 3.288 5.695 0.000) 6.576 1.222 ( 1.600 2.771 -0.000) 3.200 1.225 ( 1.411 2.445 0.000) 2.823 1.909 ( -1.621 -2.807 0.000) 3.241 1.916 ( -1.981 -3.432 0.000) 3.963 1.934 ( -1.784 -3.091 0.000) 3.569 2.124 ( 0.567 0.982 0.000) 1.134 2.128 ( 1.176 2.037 0.000) 2.352 2.148 ( 0.187 0.323 0.000) 0.373 2.199 ( 0.324 0.560 0.000) 0.647 2.218 ( -0.244 -0.422 0.000) 0.487 2.224 ( 0.120 0.208 0.000) 0.241 2.507 ( -0.098 -0.169 0.000) 0.196 2.520 ( 0.038 0.066 0.000) 0.076 2.542 ( -0.617 -1.068 0.000) 1.234 2.677 ( -2.205 -3.819 0.000) 4.410 2.682 ( -2.132 -3.693 0.000) 4.264 2.776 ( -2.620 -4.538 0.000) 5.239 2.953 ( -0.792 -1.372 0.000) 1.584 3.054 ( -2.562 -4.438 0.000) 5.124 3.056 ( -2.375 -4.114 -0.000) 4.750 ======================= Grid point 269 (40/42) ======================= q-point: ( 0.31 0.23 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 169 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.615 ( 1.303 1.824 -1.035) 2.469 0.631 ( 1.684 1.352 1.012) 2.385 0.693 ( 1.579 1.456 -0.020) 2.147 1.109 ( -0.186 5.537 -0.145) 5.542 1.122 ( -0.227 5.753 0.202) 5.761 1.149 ( -0.558 7.144 -0.044) 7.166 1.255 ( 1.307 0.517 0.056) 1.407 1.265 ( 2.105 1.162 0.064) 2.405 1.278 ( 2.716 1.107 -0.106) 2.934 1.857 ( -1.227 -2.298 0.327) 2.626 1.858 ( -1.202 -2.424 -0.355) 2.729 1.877 ( -1.185 -2.518 0.009) 2.783 2.123 ( -1.967 2.061 0.133) 2.852 2.148 ( -1.474 2.704 -0.050) 3.080 2.159 ( -0.790 3.300 -0.105) 3.395 2.205 ( 0.755 -1.272 -0.013) 1.479 2.218 ( 0.449 -0.454 -0.018) 0.638 2.229 ( 0.236 -0.011 0.065) 0.245 2.512 ( -0.011 0.705 0.018) 0.706 2.526 ( -0.354 0.162 -0.423) 0.575 2.533 ( 0.514 0.139 0.375) 0.652 2.612 ( -3.920 2.276 0.836) 4.609 2.623 ( -3.565 2.153 -0.946) 4.271 2.694 ( -5.165 2.284 0.177) 5.650 2.945 ( 3.188 -4.986 -0.425) 5.933 2.983 ( 1.647 -7.828 0.963) 8.057 3.000 ( 2.252 -8.156 -0.578) 8.481 ======================= Grid point 271 (41/42) ======================= q-point: ( 0.38 0.23 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.632 ( 0.632 -1.095 -0.000) 1.264 0.657 ( -0.037 0.064 0.000) 0.074 0.714 ( -0.719 1.246 0.000) 1.438 1.135 ( -2.753 4.769 0.000) 5.507 1.150 ( -2.872 4.974 0.000) 5.743 1.178 ( -3.370 5.837 0.000) 6.740 1.267 ( 0.799 -1.383 0.000) 1.597 1.288 ( 0.530 -0.918 -0.000) 1.060 1.311 ( 0.154 -0.267 0.000) 0.308 1.834 ( 0.528 -0.915 0.000) 1.056 1.834 ( 0.671 -1.163 -0.000) 1.342 1.853 ( 0.748 -1.296 0.000) 1.497 2.114 ( -1.692 2.930 0.000) 3.383 2.148 ( -1.689 2.925 0.000) 3.378 2.169 ( -1.972 3.415 0.000) 3.944 2.206 ( 0.655 -1.134 0.000) 1.309 2.219 ( 0.246 -0.426 0.000) 0.492 2.231 ( 0.042 -0.073 0.000) 0.084 2.513 ( -0.864 1.496 0.000) 1.727 2.523 ( -0.528 0.914 0.000) 1.056 2.540 ( -0.482 0.835 0.000) 0.964 2.592 ( -2.658 4.604 -0.000) 5.316 2.606 ( -2.448 4.240 0.000) 4.896 2.659 ( -2.976 5.154 0.000) 5.952 2.954 ( 3.808 -6.595 -0.000) 7.615 2.958 ( 4.853 -8.405 0.000) 9.705 2.973 ( 4.705 -8.150 0.000) 9.410 ======================= Grid point 297 (42/42) ======================= q-point: ( 0.31 0.31 -0.50) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.94e-05 1.94e-05 1.94e-05 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Number of triplets: 91 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.620 ( -0.013 -0.022 0.000) 0.025 0.657 ( 0.367 0.636 0.000) 0.735 0.728 ( 0.755 1.307 0.000) 1.510 1.191 ( 1.594 2.760 -0.000) 3.187 1.205 ( 1.569 2.717 0.000) 3.137 1.243 ( -0.545 -0.945 0.000) 1.091 1.246 ( 0.811 1.405 -0.000) 1.622 1.268 ( 0.404 0.699 0.000) 0.808 1.303 ( 0.222 0.384 0.000) 0.444 1.823 ( -0.570 -0.987 0.000) 1.139 1.825 ( -0.509 -0.882 0.000) 1.019 1.840 ( -0.586 -1.015 0.000) 1.172 2.146 ( 0.315 0.545 0.000) 0.629 2.175 ( -0.262 -0.454 0.000) 0.524 2.191 ( 0.377 0.654 0.000) 0.755 2.205 ( 0.765 1.326 0.000) 1.531 2.214 ( -0.093 -0.161 -0.000) 0.185 2.227 ( -0.138 -0.240 0.000) 0.277 2.522 ( 0.132 0.229 0.000) 0.264 2.529 ( 0.172 0.298 0.000) 0.344 2.545 ( 0.307 0.532 -0.000) 0.614 2.680 ( 2.423 4.196 -0.000) 4.845 2.692 ( 2.453 4.249 0.000) 4.907 2.750 ( 1.643 2.846 0.000) 3.287 2.836 ( -3.690 -6.392 0.000) 7.381 2.850 ( -3.192 -5.528 0.000) 6.383 2.863 ( -2.355 -4.078 0.000) 4.709 =================== 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/9126 10.0 15.988 15.988 3.638 -0.000 0.000 -0.000 3/9126 20.0 9.022 9.022 1.549 -0.000 0.000 -0.000 3/9126 30.0 6.028 6.028 1.025 -0.000 0.000 -0.000 3/9126 40.0 4.553 4.553 0.808 -0.000 0.000 -0.000 3/9126 50.0 3.682 3.682 0.680 -0.000 0.000 -0.000 3/9126 60.0 3.099 3.099 0.589 -0.000 0.000 -0.000 3/9126 70.0 2.678 2.678 0.519 -0.000 0.000 -0.000 3/9126 80.0 2.358 2.358 0.464 -0.000 0.000 -0.000 3/9126 90.0 2.106 2.106 0.419 -0.000 0.000 -0.000 3/9126 100.0 1.902 1.902 0.381 -0.000 0.000 -0.000 3/9126 110.0 1.735 1.735 0.350 -0.000 0.000 -0.000 3/9126 120.0 1.594 1.594 0.323 -0.000 0.000 -0.000 3/9126 130.0 1.474 1.474 0.300 -0.000 0.000 -0.000 3/9126 140.0 1.371 1.371 0.279 -0.000 0.000 -0.000 3/9126 150.0 1.281 1.281 0.262 -0.000 0.000 -0.000 3/9126 160.0 1.202 1.202 0.246 -0.000 0.000 -0.000 3/9126 170.0 1.133 1.133 0.232 -0.000 0.000 -0.000 3/9126 180.0 1.070 1.070 0.220 -0.000 0.000 -0.000 3/9126 190.0 1.015 1.015 0.209 -0.000 0.000 -0.000 3/9126 200.0 0.965 0.965 0.199 -0.000 0.000 -0.000 3/9126 210.0 0.919 0.919 0.189 -0.000 0.000 -0.000 3/9126 220.0 0.878 0.878 0.181 -0.000 0.000 -0.000 3/9126 230.0 0.840 0.840 0.173 -0.000 0.000 -0.000 3/9126 240.0 0.805 0.805 0.166 -0.000 0.000 -0.000 3/9126 250.0 0.773 0.773 0.160 -0.000 0.000 -0.000 3/9126 260.0 0.744 0.744 0.154 -0.000 0.000 -0.000 3/9126 270.0 0.716 0.716 0.148 -0.000 0.000 -0.000 3/9126 280.0 0.691 0.691 0.143 -0.000 0.000 -0.000 3/9126 290.0 0.667 0.667 0.138 -0.000 0.000 -0.000 3/9126 300.0 0.645 0.645 0.133 -0.000 0.000 -0.000 3/9126 310.0 0.624 0.624 0.129 -0.000 0.000 -0.000 3/9126 320.0 0.605 0.605 0.125 -0.000 0.000 -0.000 3/9126 330.0 0.587 0.587 0.121 -0.000 0.000 -0.000 3/9126 340.0 0.569 0.569 0.118 -0.000 0.000 -0.000 3/9126 350.0 0.553 0.553 0.115 -0.000 0.000 -0.000 3/9126 360.0 0.538 0.538 0.111 -0.000 0.000 -0.000 3/9126 370.0 0.523 0.523 0.108 -0.000 0.000 -0.000 3/9126 380.0 0.510 0.510 0.106 -0.000 0.000 -0.000 3/9126 390.0 0.497 0.497 0.103 -0.000 0.000 -0.000 3/9126 400.0 0.484 0.484 0.100 -0.000 0.000 -0.000 3/9126 410.0 0.472 0.472 0.098 -0.000 0.000 -0.000 3/9126 420.0 0.461 0.461 0.096 -0.000 0.000 -0.000 3/9126 430.0 0.451 0.451 0.093 -0.000 0.000 -0.000 3/9126 440.0 0.440 0.440 0.091 -0.000 0.000 -0.000 3/9126 450.0 0.431 0.431 0.089 -0.000 0.000 -0.000 3/9126 460.0 0.421 0.421 0.087 -0.000 0.000 -0.000 3/9126 470.0 0.412 0.412 0.086 -0.000 0.000 -0.000 3/9126 480.0 0.404 0.404 0.084 -0.000 0.000 -0.000 3/9126 490.0 0.395 0.395 0.082 -0.000 0.000 -0.000 3/9126 500.0 0.388 0.388 0.080 -0.000 0.000 -0.000 3/9126 510.0 0.380 0.380 0.079 -0.000 0.000 -0.000 3/9126 520.0 0.373 0.373 0.077 -0.000 0.000 -0.000 3/9126 530.0 0.366 0.366 0.076 -0.000 0.000 -0.000 3/9126 540.0 0.359 0.359 0.075 -0.000 0.000 -0.000 3/9126 550.0 0.352 0.352 0.073 -0.000 0.000 -0.000 3/9126 560.0 0.346 0.346 0.072 -0.000 0.000 -0.000 3/9126 570.0 0.340 0.340 0.071 -0.000 0.000 -0.000 3/9126 580.0 0.334 0.334 0.069 -0.000 0.000 -0.000 3/9126 590.0 0.329 0.329 0.068 -0.000 0.000 -0.000 3/9126 600.0 0.323 0.323 0.067 -0.000 0.000 -0.000 3/9126 610.0 0.318 0.318 0.066 -0.000 0.000 -0.000 3/9126 620.0 0.313 0.313 0.065 -0.000 0.000 -0.000 3/9126 630.0 0.308 0.308 0.064 -0.000 0.000 -0.000 3/9126 640.0 0.303 0.303 0.063 -0.000 0.000 -0.000 3/9126 650.0 0.298 0.298 0.062 -0.000 0.000 -0.000 3/9126 660.0 0.294 0.294 0.061 -0.000 0.000 -0.000 3/9126 670.0 0.289 0.289 0.060 -0.000 0.000 -0.000 3/9126 680.0 0.285 0.285 0.059 -0.000 0.000 -0.000 3/9126 690.0 0.281 0.281 0.058 -0.000 0.000 -0.000 3/9126 700.0 0.277 0.277 0.058 -0.000 0.000 -0.000 3/9126 710.0 0.273 0.273 0.057 -0.000 0.000 -0.000 3/9126 720.0 0.269 0.269 0.056 -0.000 0.000 -0.000 3/9126 730.0 0.266 0.266 0.055 -0.000 0.000 -0.000 3/9126 740.0 0.262 0.262 0.054 -0.000 0.000 -0.000 3/9126 750.0 0.259 0.259 0.054 -0.000 0.000 -0.000 3/9126 760.0 0.255 0.255 0.053 -0.000 0.000 -0.000 3/9126 770.0 0.252 0.252 0.052 -0.000 0.000 -0.000 3/9126 780.0 0.249 0.249 0.052 -0.000 0.000 -0.000 3/9126 790.0 0.245 0.245 0.051 -0.000 0.000 -0.000 3/9126 800.0 0.242 0.242 0.050 -0.000 0.000 -0.000 3/9126 810.0 0.239 0.239 0.050 -0.000 0.000 -0.000 3/9126 820.0 0.236 0.236 0.049 -0.000 0.000 -0.000 3/9126 830.0 0.234 0.234 0.049 -0.000 0.000 -0.000 3/9126 840.0 0.231 0.231 0.048 -0.000 0.000 -0.000 3/9126 850.0 0.228 0.228 0.047 -0.000 0.000 -0.000 3/9126 860.0 0.225 0.225 0.047 -0.000 0.000 -0.000 3/9126 870.0 0.223 0.223 0.046 -0.000 0.000 -0.000 3/9126 880.0 0.220 0.220 0.046 -0.000 0.000 -0.000 3/9126 890.0 0.218 0.218 0.045 -0.000 0.000 -0.000 3/9126 900.0 0.215 0.215 0.045 -0.000 0.000 -0.000 3/9126 910.0 0.213 0.213 0.044 -0.000 0.000 -0.000 3/9126 920.0 0.211 0.211 0.044 -0.000 0.000 -0.000 3/9126 930.0 0.209 0.209 0.043 -0.000 0.000 -0.000 3/9126 940.0 0.206 0.206 0.043 -0.000 0.000 -0.000 3/9126 950.0 0.204 0.204 0.042 -0.000 0.000 -0.000 3/9126 960.0 0.202 0.202 0.042 -0.000 0.000 -0.000 3/9126 970.0 0.200 0.200 0.042 -0.000 0.000 -0.000 3/9126 980.0 0.198 0.198 0.041 -0.000 0.000 -0.000 3/9126 990.0 0.196 0.196 0.041 -0.000 0.000 -0.000 3/9126 1000.0 0.194 0.194 0.040 -0.000 0.000 -0.000 3/9126 Thermal conductivity related properties were written into "kappa-m13132.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 17:35:46]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|