----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 05:39:03]------------------------- Compiled with OpenMP support (max 48 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) Number of symmetry operations in supercell: 1296 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.958652480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.958652480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.958652480000000 Atomic positions (fractional): *1 H 0.00000000000000 0.50000000000000 0.00000000000000 1.008 2 H 0.00000000000000 0.00000000000000 0.50000000000000 1.008 3 H 0.50000000000000 0.00000000000000 0.00000000000000 1.008 *4 Li 0.00000000000000 0.00000000000000 0.00000000000000 6.941 *5 Ba 0.50000000000000 0.50000000000000 0.50000000000000 137.327 -------------------------------- unit cell --------------------------------- Lattice vectors: a 3.958652480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.958652480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.958652480000000 Atomic positions (fractional): *1 H 0.00000000000000 0.50000000000000 0.00000000000000 1.008 > 1 2 H 0.00000000000000 0.00000000000000 0.50000000000000 1.008 > 2 3 H 0.50000000000000 0.00000000000000 0.00000000000000 1.008 > 3 *4 Li 0.00000000000000 0.00000000000000 0.00000000000000 6.941 > 4 *5 Ba 0.50000000000000 0.50000000000000 0.50000000000000 137.327 > 5 -------------------------------- super cell -------------------------------- Lattice vectors: a 11.875957440000001 0.000000000000000 0.000000000000000 b 0.000000000000000 11.875957440000001 0.000000000000000 c 0.000000000000000 0.000000000000000 11.875957440000001 Atomic positions (fractional): *1 H 0.00000000000000 0.16666666666667 0.00000000000000 1.008 > 1 2 H 0.33333333333333 0.16666666666667 0.00000000000000 1.008 > 1 3 H 0.66666666666667 0.16666666666667 0.00000000000000 1.008 > 1 4 H 0.00000000000000 0.50000000000000 0.00000000000000 1.008 > 1 5 H 0.33333333333333 0.50000000000000 0.00000000000000 1.008 > 1 6 H 0.66666666666667 0.50000000000000 0.00000000000000 1.008 > 1 7 H 0.00000000000000 0.83333333333333 0.00000000000000 1.008 > 1 8 H 0.33333333333333 0.83333333333333 0.00000000000000 1.008 > 1 9 H 0.66666666666667 0.83333333333333 0.00000000000000 1.008 > 1 10 H 0.00000000000000 0.16666666666667 0.33333333333333 1.008 > 1 11 H 0.33333333333333 0.16666666666667 0.33333333333333 1.008 > 1 12 H 0.66666666666667 0.16666666666667 0.33333333333333 1.008 > 1 13 H 0.00000000000000 0.50000000000000 0.33333333333333 1.008 > 1 14 H 0.33333333333333 0.50000000000000 0.33333333333333 1.008 > 1 15 H 0.66666666666667 0.50000000000000 0.33333333333333 1.008 > 1 16 H 0.00000000000000 0.83333333333333 0.33333333333333 1.008 > 1 17 H 0.33333333333333 0.83333333333333 0.33333333333333 1.008 > 1 18 H 0.66666666666667 0.83333333333333 0.33333333333333 1.008 > 1 19 H 0.00000000000000 0.16666666666667 0.66666666666667 1.008 > 1 20 H 0.33333333333333 0.16666666666667 0.66666666666667 1.008 > 1 21 H 0.66666666666667 0.16666666666667 0.66666666666667 1.008 > 1 22 H 0.00000000000000 0.50000000000000 0.66666666666667 1.008 > 1 23 H 0.33333333333333 0.50000000000000 0.66666666666667 1.008 > 1 24 H 0.66666666666667 0.50000000000000 0.66666666666667 1.008 > 1 25 H 0.00000000000000 0.83333333333333 0.66666666666667 1.008 > 1 26 H 0.33333333333333 0.83333333333333 0.66666666666667 1.008 > 1 27 H 0.66666666666667 0.83333333333333 0.66666666666667 1.008 > 1 28 H 0.00000000000000 0.00000000000000 0.16666666666667 1.008 > 2 29 H 0.33333333333333 0.00000000000000 0.16666666666667 1.008 > 2 30 H 0.66666666666667 0.00000000000000 0.16666666666667 1.008 > 2 31 H 0.00000000000000 0.33333333333333 0.16666666666667 1.008 > 2 32 H 0.33333333333333 0.33333333333333 0.16666666666667 1.008 > 2 33 H 0.66666666666667 0.33333333333333 0.16666666666667 1.008 > 2 34 H 0.00000000000000 0.66666666666667 0.16666666666667 1.008 > 2 35 H 0.33333333333333 0.66666666666667 0.16666666666667 1.008 > 2 36 H 0.66666666666667 0.66666666666667 0.16666666666667 1.008 > 2 37 H 0.00000000000000 0.00000000000000 0.50000000000000 1.008 > 2 38 H 0.33333333333333 0.00000000000000 0.50000000000000 1.008 > 2 39 H 0.66666666666667 0.00000000000000 0.50000000000000 1.008 > 2 40 H 0.00000000000000 0.33333333333333 0.50000000000000 1.008 > 2 41 H 0.33333333333333 0.33333333333333 0.50000000000000 1.008 > 2 42 H 0.66666666666667 0.33333333333333 0.50000000000000 1.008 > 2 43 H 0.00000000000000 0.66666666666667 0.50000000000000 1.008 > 2 44 H 0.33333333333333 0.66666666666667 0.50000000000000 1.008 > 2 45 H 0.66666666666667 0.66666666666667 0.50000000000000 1.008 > 2 46 H 0.00000000000000 0.00000000000000 0.83333333333333 1.008 > 2 47 H 0.33333333333333 0.00000000000000 0.83333333333333 1.008 > 2 48 H 0.66666666666667 0.00000000000000 0.83333333333333 1.008 > 2 49 H 0.00000000000000 0.33333333333333 0.83333333333333 1.008 > 2 50 H 0.33333333333333 0.33333333333333 0.83333333333333 1.008 > 2 51 H 0.66666666666667 0.33333333333333 0.83333333333333 1.008 > 2 52 H 0.00000000000000 0.66666666666667 0.83333333333333 1.008 > 2 53 H 0.33333333333333 0.66666666666667 0.83333333333333 1.008 > 2 54 H 0.66666666666667 0.66666666666667 0.83333333333333 1.008 > 2 55 H 0.16666666666667 0.00000000000000 0.00000000000000 1.008 > 3 56 H 0.50000000000000 0.00000000000000 0.00000000000000 1.008 > 3 57 H 0.83333333333333 0.00000000000000 0.00000000000000 1.008 > 3 58 H 0.16666666666667 0.33333333333333 0.00000000000000 1.008 > 3 59 H 0.50000000000000 0.33333333333333 0.00000000000000 1.008 > 3 60 H 0.83333333333333 0.33333333333333 0.00000000000000 1.008 > 3 61 H 0.16666666666667 0.66666666666667 0.00000000000000 1.008 > 3 62 H 0.50000000000000 0.66666666666667 0.00000000000000 1.008 > 3 63 H 0.83333333333333 0.66666666666667 0.00000000000000 1.008 > 3 64 H 0.16666666666667 0.00000000000000 0.33333333333333 1.008 > 3 65 H 0.50000000000000 0.00000000000000 0.33333333333333 1.008 > 3 66 H 0.83333333333333 0.00000000000000 0.33333333333333 1.008 > 3 67 H 0.16666666666667 0.33333333333333 0.33333333333333 1.008 > 3 68 H 0.50000000000000 0.33333333333333 0.33333333333333 1.008 > 3 69 H 0.83333333333333 0.33333333333333 0.33333333333333 1.008 > 3 70 H 0.16666666666667 0.66666666666667 0.33333333333333 1.008 > 3 71 H 0.50000000000000 0.66666666666667 0.33333333333333 1.008 > 3 72 H 0.83333333333333 0.66666666666667 0.33333333333333 1.008 > 3 73 H 0.16666666666667 0.00000000000000 0.66666666666667 1.008 > 3 74 H 0.50000000000000 0.00000000000000 0.66666666666667 1.008 > 3 75 H 0.83333333333333 0.00000000000000 0.66666666666667 1.008 > 3 76 H 0.16666666666667 0.33333333333333 0.66666666666667 1.008 > 3 77 H 0.50000000000000 0.33333333333333 0.66666666666667 1.008 > 3 78 H 0.83333333333333 0.33333333333333 0.66666666666667 1.008 > 3 79 H 0.16666666666667 0.66666666666667 0.66666666666667 1.008 > 3 80 H 0.50000000000000 0.66666666666667 0.66666666666667 1.008 > 3 81 H 0.83333333333333 0.66666666666667 0.66666666666667 1.008 > 3 *82 Li 0.00000000000000 0.00000000000000 0.00000000000000 6.941 > 4 83 Li 0.33333333333333 0.00000000000000 0.00000000000000 6.941 > 4 84 Li 0.66666666666667 0.00000000000000 0.00000000000000 6.941 > 4 85 Li 0.00000000000000 0.33333333333333 0.00000000000000 6.941 > 4 86 Li 0.33333333333333 0.33333333333333 0.00000000000000 6.941 > 4 87 Li 0.66666666666667 0.33333333333333 0.00000000000000 6.941 > 4 88 Li 0.00000000000000 0.66666666666667 0.00000000000000 6.941 > 4 89 Li 0.33333333333333 0.66666666666667 0.00000000000000 6.941 > 4 90 Li 0.66666666666667 0.66666666666667 0.00000000000000 6.941 > 4 91 Li 0.00000000000000 0.00000000000000 0.33333333333333 6.941 > 4 92 Li 0.33333333333333 0.00000000000000 0.33333333333333 6.941 > 4 93 Li 0.66666666666667 0.00000000000000 0.33333333333333 6.941 > 4 94 Li 0.00000000000000 0.33333333333333 0.33333333333333 6.941 > 4 95 Li 0.33333333333333 0.33333333333333 0.33333333333333 6.941 > 4 96 Li 0.66666666666667 0.33333333333333 0.33333333333333 6.941 > 4 97 Li 0.00000000000000 0.66666666666667 0.33333333333333 6.941 > 4 98 Li 0.33333333333333 0.66666666666667 0.33333333333333 6.941 > 4 99 Li 0.66666666666667 0.66666666666667 0.33333333333333 6.941 > 4 100 Li 0.00000000000000 0.00000000000000 0.66666666666667 6.941 > 4 101 Li 0.33333333333333 0.00000000000000 0.66666666666667 6.941 > 4 102 Li 0.66666666666667 0.00000000000000 0.66666666666667 6.941 > 4 103 Li 0.00000000000000 0.33333333333333 0.66666666666667 6.941 > 4 104 Li 0.33333333333333 0.33333333333333 0.66666666666667 6.941 > 4 105 Li 0.66666666666667 0.33333333333333 0.66666666666667 6.941 > 4 106 Li 0.00000000000000 0.66666666666667 0.66666666666667 6.941 > 4 107 Li 0.33333333333333 0.66666666666667 0.66666666666667 6.941 > 4 108 Li 0.66666666666667 0.66666666666667 0.66666666666667 6.941 > 4 *109 Ba 0.16666666666667 0.16666666666667 0.16666666666667 137.327 > 5 110 Ba 0.50000000000000 0.16666666666667 0.16666666666667 137.327 > 5 111 Ba 0.83333333333333 0.16666666666667 0.16666666666667 137.327 > 5 112 Ba 0.16666666666667 0.50000000000000 0.16666666666667 137.327 > 5 113 Ba 0.50000000000000 0.50000000000000 0.16666666666667 137.327 > 5 114 Ba 0.83333333333333 0.50000000000000 0.16666666666667 137.327 > 5 115 Ba 0.16666666666667 0.83333333333333 0.16666666666667 137.327 > 5 116 Ba 0.50000000000000 0.83333333333333 0.16666666666667 137.327 > 5 117 Ba 0.83333333333333 0.83333333333333 0.16666666666667 137.327 > 5 118 Ba 0.16666666666667 0.16666666666667 0.50000000000000 137.327 > 5 119 Ba 0.50000000000000 0.16666666666667 0.50000000000000 137.327 > 5 120 Ba 0.83333333333333 0.16666666666667 0.50000000000000 137.327 > 5 121 Ba 0.16666666666667 0.50000000000000 0.50000000000000 137.327 > 5 122 Ba 0.50000000000000 0.50000000000000 0.50000000000000 137.327 > 5 123 Ba 0.83333333333333 0.50000000000000 0.50000000000000 137.327 > 5 124 Ba 0.16666666666667 0.83333333333333 0.50000000000000 137.327 > 5 125 Ba 0.50000000000000 0.83333333333333 0.50000000000000 137.327 > 5 126 Ba 0.83333333333333 0.83333333333333 0.50000000000000 137.327 > 5 127 Ba 0.16666666666667 0.16666666666667 0.83333333333333 137.327 > 5 128 Ba 0.50000000000000 0.16666666666667 0.83333333333333 137.327 > 5 129 Ba 0.83333333333333 0.16666666666667 0.83333333333333 137.327 > 5 130 Ba 0.16666666666667 0.50000000000000 0.83333333333333 137.327 > 5 131 Ba 0.50000000000000 0.50000000000000 0.83333333333333 137.327 > 5 132 Ba 0.83333333333333 0.50000000000000 0.83333333333333 137.327 > 5 133 Ba 0.16666666666667 0.83333333333333 0.83333333333333 137.327 > 5 134 Ba 0.50000000000000 0.83333333333333 0.83333333333333 137.327 > 5 135 Ba 0.83333333333333 0.83333333333333 0.83333333333333 137.327 > 5 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.3490772 0.0000000 0.0000000 0.0000000 5.3490772 0.0000000 0.0000000 0.0000000 5.3490772 -------------------------- Born effective charges -------------------------- 1 H -1.1841729 0.0000000 0.0000000 0.0000000 -0.8761288 0.0000000 0.0000000 0.0000000 -1.1841729 2 H -1.1841729 0.0000000 0.0000000 0.0000000 -1.1841729 0.0000000 0.0000000 0.0000000 -0.8761288 3 H -0.8761288 0.0000000 0.0000000 0.0000000 -1.1841729 0.0000000 0.0000000 0.0000000 -1.1841729 4 Li 0.8715575 0.0000000 0.0000000 0.0000000 0.8715575 0.0000000 0.0000000 0.0000000 0.8715575 5 Ba 2.3729170 0.0000000 0.0000000 0.0000000 2.3729170 0.0000000 0.0000000 0.0000000 2.3729170 ---------------------------------------------------------------------------- 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: 405/405 Permutation basis: 3615/3615 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 104 Number of blocks in projector: 104 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 74 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 30 Use standard eigh solver. Tree of FC basis block matrices: - (104, 100), data: False |-- (30, 30), data: True |-- (74, 70), data: True ----- Solver_atoms: 1 -- 135 / 135 Time (Solver_compr_matrix_reshape): 0.003 Solver_block: 80 / 80 - Time: 0.055 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.059 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 405/405 Permutation basis: 3615/3615 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 104 Number of blocks in projector: 104 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 2 --- Eigsh_solver_block: 1 / 2 --- Block_size: 74 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 30 Use standard eigh solver. Tree of FC basis block matrices: - (104, 100), data: False |-- (30, 30), data: True |-- (74, 70), 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 05:39:07]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 05:39:08]------------------------- Compiled with OpenMP support (max 48 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.958652480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.958652480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.958652480000000 Atomic positions (fractional): 1 H 0.00000000000000 0.50000000000000 0.00000000000000 1.008 2 H 0.00000000000000 0.00000000000000 0.50000000000000 1.008 3 H 0.50000000000000 0.00000000000000 0.00000000000000 1.008 4 Li 0.00000000000000 0.00000000000000 0.00000000000000 6.941 5 Ba 0.50000000000000 0.50000000000000 0.50000000000000 137.327 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.875957440000001 0.000000000000000 0.000000000000000 b 0.000000000000000 11.875957440000001 0.000000000000000 c 0.000000000000000 0.000000000000000 11.875957440000001 Atomic positions (fractional): 1 H 0.00000000000000 0.16666666666667 0.00000000000000 1.008 > 1 2 H 0.33333333333333 0.16666666666667 0.00000000000000 1.008 > 1 3 H 0.66666666666667 0.16666666666667 0.00000000000000 1.008 > 1 4 H 0.00000000000000 0.50000000000000 0.00000000000000 1.008 > 1 5 H 0.33333333333333 0.50000000000000 0.00000000000000 1.008 > 1 6 H 0.66666666666667 0.50000000000000 0.00000000000000 1.008 > 1 7 H 0.00000000000000 0.83333333333333 0.00000000000000 1.008 > 1 8 H 0.33333333333333 0.83333333333333 0.00000000000000 1.008 > 1 9 H 0.66666666666667 0.83333333333333 0.00000000000000 1.008 > 1 10 H 0.00000000000000 0.16666666666667 0.33333333333333 1.008 > 1 11 H 0.33333333333333 0.16666666666667 0.33333333333333 1.008 > 1 12 H 0.66666666666667 0.16666666666667 0.33333333333333 1.008 > 1 13 H 0.00000000000000 0.50000000000000 0.33333333333333 1.008 > 1 14 H 0.33333333333333 0.50000000000000 0.33333333333333 1.008 > 1 15 H 0.66666666666667 0.50000000000000 0.33333333333333 1.008 > 1 16 H 0.00000000000000 0.83333333333333 0.33333333333333 1.008 > 1 17 H 0.33333333333333 0.83333333333333 0.33333333333333 1.008 > 1 18 H 0.66666666666667 0.83333333333333 0.33333333333333 1.008 > 1 19 H 0.00000000000000 0.16666666666667 0.66666666666667 1.008 > 1 20 H 0.33333333333333 0.16666666666667 0.66666666666667 1.008 > 1 21 H 0.66666666666667 0.16666666666667 0.66666666666667 1.008 > 1 22 H 0.00000000000000 0.50000000000000 0.66666666666667 1.008 > 1 23 H 0.33333333333333 0.50000000000000 0.66666666666667 1.008 > 1 24 H 0.66666666666667 0.50000000000000 0.66666666666667 1.008 > 1 25 H 0.00000000000000 0.83333333333333 0.66666666666667 1.008 > 1 26 H 0.33333333333333 0.83333333333333 0.66666666666667 1.008 > 1 27 H 0.66666666666667 0.83333333333333 0.66666666666667 1.008 > 1 28 H 0.00000000000000 0.00000000000000 0.16666666666667 1.008 > 28 29 H 0.33333333333333 0.00000000000000 0.16666666666667 1.008 > 28 30 H 0.66666666666667 0.00000000000000 0.16666666666667 1.008 > 28 31 H 0.00000000000000 0.33333333333333 0.16666666666667 1.008 > 28 32 H 0.33333333333333 0.33333333333333 0.16666666666667 1.008 > 28 33 H 0.66666666666667 0.33333333333333 0.16666666666667 1.008 > 28 34 H 0.00000000000000 0.66666666666667 0.16666666666667 1.008 > 28 35 H 0.33333333333333 0.66666666666667 0.16666666666667 1.008 > 28 36 H 0.66666666666667 0.66666666666667 0.16666666666667 1.008 > 28 37 H 0.00000000000000 0.00000000000000 0.50000000000000 1.008 > 28 38 H 0.33333333333333 0.00000000000000 0.50000000000000 1.008 > 28 39 H 0.66666666666667 0.00000000000000 0.50000000000000 1.008 > 28 40 H 0.00000000000000 0.33333333333333 0.50000000000000 1.008 > 28 41 H 0.33333333333333 0.33333333333333 0.50000000000000 1.008 > 28 42 H 0.66666666666667 0.33333333333333 0.50000000000000 1.008 > 28 43 H 0.00000000000000 0.66666666666667 0.50000000000000 1.008 > 28 44 H 0.33333333333333 0.66666666666667 0.50000000000000 1.008 > 28 45 H 0.66666666666667 0.66666666666667 0.50000000000000 1.008 > 28 46 H 0.00000000000000 0.00000000000000 0.83333333333333 1.008 > 28 47 H 0.33333333333333 0.00000000000000 0.83333333333333 1.008 > 28 48 H 0.66666666666667 0.00000000000000 0.83333333333333 1.008 > 28 49 H 0.00000000000000 0.33333333333333 0.83333333333333 1.008 > 28 50 H 0.33333333333333 0.33333333333333 0.83333333333333 1.008 > 28 51 H 0.66666666666667 0.33333333333333 0.83333333333333 1.008 > 28 52 H 0.00000000000000 0.66666666666667 0.83333333333333 1.008 > 28 53 H 0.33333333333333 0.66666666666667 0.83333333333333 1.008 > 28 54 H 0.66666666666667 0.66666666666667 0.83333333333333 1.008 > 28 55 H 0.16666666666667 0.00000000000000 0.00000000000000 1.008 > 55 56 H 0.50000000000000 0.00000000000000 0.00000000000000 1.008 > 55 57 H 0.83333333333333 0.00000000000000 0.00000000000000 1.008 > 55 58 H 0.16666666666667 0.33333333333333 0.00000000000000 1.008 > 55 59 H 0.50000000000000 0.33333333333333 0.00000000000000 1.008 > 55 60 H 0.83333333333333 0.33333333333333 0.00000000000000 1.008 > 55 61 H 0.16666666666667 0.66666666666667 0.00000000000000 1.008 > 55 62 H 0.50000000000000 0.66666666666667 0.00000000000000 1.008 > 55 63 H 0.83333333333333 0.66666666666667 0.00000000000000 1.008 > 55 64 H 0.16666666666667 0.00000000000000 0.33333333333333 1.008 > 55 65 H 0.50000000000000 0.00000000000000 0.33333333333333 1.008 > 55 66 H 0.83333333333333 0.00000000000000 0.33333333333333 1.008 > 55 67 H 0.16666666666667 0.33333333333333 0.33333333333333 1.008 > 55 68 H 0.50000000000000 0.33333333333333 0.33333333333333 1.008 > 55 69 H 0.83333333333333 0.33333333333333 0.33333333333333 1.008 > 55 70 H 0.16666666666667 0.66666666666667 0.33333333333333 1.008 > 55 71 H 0.50000000000000 0.66666666666667 0.33333333333333 1.008 > 55 72 H 0.83333333333333 0.66666666666667 0.33333333333333 1.008 > 55 73 H 0.16666666666667 0.00000000000000 0.66666666666667 1.008 > 55 74 H 0.50000000000000 0.00000000000000 0.66666666666667 1.008 > 55 75 H 0.83333333333333 0.00000000000000 0.66666666666667 1.008 > 55 76 H 0.16666666666667 0.33333333333333 0.66666666666667 1.008 > 55 77 H 0.50000000000000 0.33333333333333 0.66666666666667 1.008 > 55 78 H 0.83333333333333 0.33333333333333 0.66666666666667 1.008 > 55 79 H 0.16666666666667 0.66666666666667 0.66666666666667 1.008 > 55 80 H 0.50000000000000 0.66666666666667 0.66666666666667 1.008 > 55 81 H 0.83333333333333 0.66666666666667 0.66666666666667 1.008 > 55 82 Li 0.00000000000000 0.00000000000000 0.00000000000000 6.941 > 82 83 Li 0.33333333333333 0.00000000000000 0.00000000000000 6.941 > 82 84 Li 0.66666666666667 0.00000000000000 0.00000000000000 6.941 > 82 85 Li 0.00000000000000 0.33333333333333 0.00000000000000 6.941 > 82 86 Li 0.33333333333333 0.33333333333333 0.00000000000000 6.941 > 82 87 Li 0.66666666666667 0.33333333333333 0.00000000000000 6.941 > 82 88 Li 0.00000000000000 0.66666666666667 0.00000000000000 6.941 > 82 89 Li 0.33333333333333 0.66666666666667 0.00000000000000 6.941 > 82 90 Li 0.66666666666667 0.66666666666667 0.00000000000000 6.941 > 82 91 Li 0.00000000000000 0.00000000000000 0.33333333333333 6.941 > 82 92 Li 0.33333333333333 0.00000000000000 0.33333333333333 6.941 > 82 93 Li 0.66666666666667 0.00000000000000 0.33333333333333 6.941 > 82 94 Li 0.00000000000000 0.33333333333333 0.33333333333333 6.941 > 82 95 Li 0.33333333333333 0.33333333333333 0.33333333333333 6.941 > 82 96 Li 0.66666666666667 0.33333333333333 0.33333333333333 6.941 > 82 97 Li 0.00000000000000 0.66666666666667 0.33333333333333 6.941 > 82 98 Li 0.33333333333333 0.66666666666667 0.33333333333333 6.941 > 82 99 Li 0.66666666666667 0.66666666666667 0.33333333333333 6.941 > 82 100 Li 0.00000000000000 0.00000000000000 0.66666666666667 6.941 > 82 101 Li 0.33333333333333 0.00000000000000 0.66666666666667 6.941 > 82 102 Li 0.66666666666667 0.00000000000000 0.66666666666667 6.941 > 82 103 Li 0.00000000000000 0.33333333333333 0.66666666666667 6.941 > 82 104 Li 0.33333333333333 0.33333333333333 0.66666666666667 6.941 > 82 105 Li 0.66666666666667 0.33333333333333 0.66666666666667 6.941 > 82 106 Li 0.00000000000000 0.66666666666667 0.66666666666667 6.941 > 82 107 Li 0.33333333333333 0.66666666666667 0.66666666666667 6.941 > 82 108 Li 0.66666666666667 0.66666666666667 0.66666666666667 6.941 > 82 109 Ba 0.16666666666667 0.16666666666667 0.16666666666667 137.327 > 109 110 Ba 0.50000000000000 0.16666666666667 0.16666666666667 137.327 > 109 111 Ba 0.83333333333333 0.16666666666667 0.16666666666667 137.327 > 109 112 Ba 0.16666666666667 0.50000000000000 0.16666666666667 137.327 > 109 113 Ba 0.50000000000000 0.50000000000000 0.16666666666667 137.327 > 109 114 Ba 0.83333333333333 0.50000000000000 0.16666666666667 137.327 > 109 115 Ba 0.16666666666667 0.83333333333333 0.16666666666667 137.327 > 109 116 Ba 0.50000000000000 0.83333333333333 0.16666666666667 137.327 > 109 117 Ba 0.83333333333333 0.83333333333333 0.16666666666667 137.327 > 109 118 Ba 0.16666666666667 0.16666666666667 0.50000000000000 137.327 > 109 119 Ba 0.50000000000000 0.16666666666667 0.50000000000000 137.327 > 109 120 Ba 0.83333333333333 0.16666666666667 0.50000000000000 137.327 > 109 121 Ba 0.16666666666667 0.50000000000000 0.50000000000000 137.327 > 109 122 Ba 0.50000000000000 0.50000000000000 0.50000000000000 137.327 > 109 123 Ba 0.83333333333333 0.50000000000000 0.50000000000000 137.327 > 109 124 Ba 0.16666666666667 0.83333333333333 0.50000000000000 137.327 > 109 125 Ba 0.50000000000000 0.83333333333333 0.50000000000000 137.327 > 109 126 Ba 0.83333333333333 0.83333333333333 0.50000000000000 137.327 > 109 127 Ba 0.16666666666667 0.16666666666667 0.83333333333333 137.327 > 109 128 Ba 0.50000000000000 0.16666666666667 0.83333333333333 137.327 > 109 129 Ba 0.83333333333333 0.16666666666667 0.83333333333333 137.327 > 109 130 Ba 0.16666666666667 0.50000000000000 0.83333333333333 137.327 > 109 131 Ba 0.50000000000000 0.50000000000000 0.83333333333333 137.327 > 109 132 Ba 0.83333333333333 0.50000000000000 0.83333333333333 137.327 > 109 133 Ba 0.16666666666667 0.83333333333333 0.83333333333333 137.327 > 109 134 Ba 0.50000000000000 0.83333333333333 0.83333333333333 137.327 > 109 135 Ba 0.83333333333333 0.83333333333333 0.83333333333333 137.327 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 5.3490772 0.0000000 0.0000000 0.0000000 5.3490772 0.0000000 0.0000000 0.0000000 5.3490772 -------------------------- Born effective charges -------------------------- 1 H -1.1841729 0.0000000 0.0000000 0.0000000 -0.8761288 0.0000000 0.0000000 0.0000000 -1.1841729 2 H -1.1841729 0.0000000 0.0000000 0.0000000 -1.1841729 0.0000000 0.0000000 0.0000000 -0.8761288 3 H -0.8761288 0.0000000 0.0000000 0.0000000 -1.1841729 0.0000000 0.0000000 0.0000000 -1.1841729 4 Li 0.8715575 0.0000000 0.0000000 0.0000000 0.8715575 0.0000000 0.0000000 0.0000000 0.8715575 5 Ba 2.3729170 0.0000000 0.0000000 0.0000000 2.3729170 0.0000000 0.0000000 0.0000000 2.3729170 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] Computing fc3[ 82, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Computing fc3[ 109, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: -0.00000000 (xxx) -0.00000000 (xxx) -0.00000000 (xxx) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (yy) -0.00000000 (yy) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 05:39:11]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 05:39:11]------------------------- Compiled with OpenMP support (max 48 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [3 3 3] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.958652480000000 0.000000000000000 0.000000000000000 b 0.000000000000000 3.958652480000000 0.000000000000000 c 0.000000000000000 0.000000000000000 3.958652480000000 Atomic positions (fractional): 1 H 0.00000000000000 0.50000000000000 0.00000000000000 1.008 2 H 0.00000000000000 0.00000000000000 0.50000000000000 1.008 3 H 0.50000000000000 0.00000000000000 0.00000000000000 1.008 4 Li 0.00000000000000 0.00000000000000 0.00000000000000 6.941 5 Ba 0.50000000000000 0.50000000000000 0.50000000000000 137.327 -------------------------------- supercell --------------------------------- Lattice vectors: a 11.875957440000001 0.000000000000000 0.000000000000000 b 0.000000000000000 11.875957440000001 0.000000000000000 c 0.000000000000000 0.000000000000000 11.875957440000001 Atomic positions (fractional): 1 H 0.00000000000000 0.16666666666667 0.00000000000000 1.008 > 1 2 H 0.33333333333333 0.16666666666667 0.00000000000000 1.008 > 1 3 H 0.66666666666667 0.16666666666667 0.00000000000000 1.008 > 1 4 H 0.00000000000000 0.50000000000000 0.00000000000000 1.008 > 1 5 H 0.33333333333333 0.50000000000000 0.00000000000000 1.008 > 1 6 H 0.66666666666667 0.50000000000000 0.00000000000000 1.008 > 1 7 H 0.00000000000000 0.83333333333333 0.00000000000000 1.008 > 1 8 H 0.33333333333333 0.83333333333333 0.00000000000000 1.008 > 1 9 H 0.66666666666667 0.83333333333333 0.00000000000000 1.008 > 1 10 H 0.00000000000000 0.16666666666667 0.33333333333333 1.008 > 1 11 H 0.33333333333333 0.16666666666667 0.33333333333333 1.008 > 1 12 H 0.66666666666667 0.16666666666667 0.33333333333333 1.008 > 1 13 H 0.00000000000000 0.50000000000000 0.33333333333333 1.008 > 1 14 H 0.33333333333333 0.50000000000000 0.33333333333333 1.008 > 1 15 H 0.66666666666667 0.50000000000000 0.33333333333333 1.008 > 1 16 H 0.00000000000000 0.83333333333333 0.33333333333333 1.008 > 1 17 H 0.33333333333333 0.83333333333333 0.33333333333333 1.008 > 1 18 H 0.66666666666667 0.83333333333333 0.33333333333333 1.008 > 1 19 H 0.00000000000000 0.16666666666667 0.66666666666667 1.008 > 1 20 H 0.33333333333333 0.16666666666667 0.66666666666667 1.008 > 1 21 H 0.66666666666667 0.16666666666667 0.66666666666667 1.008 > 1 22 H 0.00000000000000 0.50000000000000 0.66666666666667 1.008 > 1 23 H 0.33333333333333 0.50000000000000 0.66666666666667 1.008 > 1 24 H 0.66666666666667 0.50000000000000 0.66666666666667 1.008 > 1 25 H 0.00000000000000 0.83333333333333 0.66666666666667 1.008 > 1 26 H 0.33333333333333 0.83333333333333 0.66666666666667 1.008 > 1 27 H 0.66666666666667 0.83333333333333 0.66666666666667 1.008 > 1 28 H 0.00000000000000 0.00000000000000 0.16666666666667 1.008 > 28 29 H 0.33333333333333 0.00000000000000 0.16666666666667 1.008 > 28 30 H 0.66666666666667 0.00000000000000 0.16666666666667 1.008 > 28 31 H 0.00000000000000 0.33333333333333 0.16666666666667 1.008 > 28 32 H 0.33333333333333 0.33333333333333 0.16666666666667 1.008 > 28 33 H 0.66666666666667 0.33333333333333 0.16666666666667 1.008 > 28 34 H 0.00000000000000 0.66666666666667 0.16666666666667 1.008 > 28 35 H 0.33333333333333 0.66666666666667 0.16666666666667 1.008 > 28 36 H 0.66666666666667 0.66666666666667 0.16666666666667 1.008 > 28 37 H 0.00000000000000 0.00000000000000 0.50000000000000 1.008 > 28 38 H 0.33333333333333 0.00000000000000 0.50000000000000 1.008 > 28 39 H 0.66666666666667 0.00000000000000 0.50000000000000 1.008 > 28 40 H 0.00000000000000 0.33333333333333 0.50000000000000 1.008 > 28 41 H 0.33333333333333 0.33333333333333 0.50000000000000 1.008 > 28 42 H 0.66666666666667 0.33333333333333 0.50000000000000 1.008 > 28 43 H 0.00000000000000 0.66666666666667 0.50000000000000 1.008 > 28 44 H 0.33333333333333 0.66666666666667 0.50000000000000 1.008 > 28 45 H 0.66666666666667 0.66666666666667 0.50000000000000 1.008 > 28 46 H 0.00000000000000 0.00000000000000 0.83333333333333 1.008 > 28 47 H 0.33333333333333 0.00000000000000 0.83333333333333 1.008 > 28 48 H 0.66666666666667 0.00000000000000 0.83333333333333 1.008 > 28 49 H 0.00000000000000 0.33333333333333 0.83333333333333 1.008 > 28 50 H 0.33333333333333 0.33333333333333 0.83333333333333 1.008 > 28 51 H 0.66666666666667 0.33333333333333 0.83333333333333 1.008 > 28 52 H 0.00000000000000 0.66666666666667 0.83333333333333 1.008 > 28 53 H 0.33333333333333 0.66666666666667 0.83333333333333 1.008 > 28 54 H 0.66666666666667 0.66666666666667 0.83333333333333 1.008 > 28 55 H 0.16666666666667 0.00000000000000 0.00000000000000 1.008 > 55 56 H 0.50000000000000 0.00000000000000 0.00000000000000 1.008 > 55 57 H 0.83333333333333 0.00000000000000 0.00000000000000 1.008 > 55 58 H 0.16666666666667 0.33333333333333 0.00000000000000 1.008 > 55 59 H 0.50000000000000 0.33333333333333 0.00000000000000 1.008 > 55 60 H 0.83333333333333 0.33333333333333 0.00000000000000 1.008 > 55 61 H 0.16666666666667 0.66666666666667 0.00000000000000 1.008 > 55 62 H 0.50000000000000 0.66666666666667 0.00000000000000 1.008 > 55 63 H 0.83333333333333 0.66666666666667 0.00000000000000 1.008 > 55 64 H 0.16666666666667 0.00000000000000 0.33333333333333 1.008 > 55 65 H 0.50000000000000 0.00000000000000 0.33333333333333 1.008 > 55 66 H 0.83333333333333 0.00000000000000 0.33333333333333 1.008 > 55 67 H 0.16666666666667 0.33333333333333 0.33333333333333 1.008 > 55 68 H 0.50000000000000 0.33333333333333 0.33333333333333 1.008 > 55 69 H 0.83333333333333 0.33333333333333 0.33333333333333 1.008 > 55 70 H 0.16666666666667 0.66666666666667 0.33333333333333 1.008 > 55 71 H 0.50000000000000 0.66666666666667 0.33333333333333 1.008 > 55 72 H 0.83333333333333 0.66666666666667 0.33333333333333 1.008 > 55 73 H 0.16666666666667 0.00000000000000 0.66666666666667 1.008 > 55 74 H 0.50000000000000 0.00000000000000 0.66666666666667 1.008 > 55 75 H 0.83333333333333 0.00000000000000 0.66666666666667 1.008 > 55 76 H 0.16666666666667 0.33333333333333 0.66666666666667 1.008 > 55 77 H 0.50000000000000 0.33333333333333 0.66666666666667 1.008 > 55 78 H 0.83333333333333 0.33333333333333 0.66666666666667 1.008 > 55 79 H 0.16666666666667 0.66666666666667 0.66666666666667 1.008 > 55 80 H 0.50000000000000 0.66666666666667 0.66666666666667 1.008 > 55 81 H 0.83333333333333 0.66666666666667 0.66666666666667 1.008 > 55 82 Li 0.00000000000000 0.00000000000000 0.00000000000000 6.941 > 82 83 Li 0.33333333333333 0.00000000000000 0.00000000000000 6.941 > 82 84 Li 0.66666666666667 0.00000000000000 0.00000000000000 6.941 > 82 85 Li 0.00000000000000 0.33333333333333 0.00000000000000 6.941 > 82 86 Li 0.33333333333333 0.33333333333333 0.00000000000000 6.941 > 82 87 Li 0.66666666666667 0.33333333333333 0.00000000000000 6.941 > 82 88 Li 0.00000000000000 0.66666666666667 0.00000000000000 6.941 > 82 89 Li 0.33333333333333 0.66666666666667 0.00000000000000 6.941 > 82 90 Li 0.66666666666667 0.66666666666667 0.00000000000000 6.941 > 82 91 Li 0.00000000000000 0.00000000000000 0.33333333333333 6.941 > 82 92 Li 0.33333333333333 0.00000000000000 0.33333333333333 6.941 > 82 93 Li 0.66666666666667 0.00000000000000 0.33333333333333 6.941 > 82 94 Li 0.00000000000000 0.33333333333333 0.33333333333333 6.941 > 82 95 Li 0.33333333333333 0.33333333333333 0.33333333333333 6.941 > 82 96 Li 0.66666666666667 0.33333333333333 0.33333333333333 6.941 > 82 97 Li 0.00000000000000 0.66666666666667 0.33333333333333 6.941 > 82 98 Li 0.33333333333333 0.66666666666667 0.33333333333333 6.941 > 82 99 Li 0.66666666666667 0.66666666666667 0.33333333333333 6.941 > 82 100 Li 0.00000000000000 0.00000000000000 0.66666666666667 6.941 > 82 101 Li 0.33333333333333 0.00000000000000 0.66666666666667 6.941 > 82 102 Li 0.66666666666667 0.00000000000000 0.66666666666667 6.941 > 82 103 Li 0.00000000000000 0.33333333333333 0.66666666666667 6.941 > 82 104 Li 0.33333333333333 0.33333333333333 0.66666666666667 6.941 > 82 105 Li 0.66666666666667 0.33333333333333 0.66666666666667 6.941 > 82 106 Li 0.00000000000000 0.66666666666667 0.66666666666667 6.941 > 82 107 Li 0.33333333333333 0.66666666666667 0.66666666666667 6.941 > 82 108 Li 0.66666666666667 0.66666666666667 0.66666666666667 6.941 > 82 109 Ba 0.16666666666667 0.16666666666667 0.16666666666667 137.327 > 109 110 Ba 0.50000000000000 0.16666666666667 0.16666666666667 137.327 > 109 111 Ba 0.83333333333333 0.16666666666667 0.16666666666667 137.327 > 109 112 Ba 0.16666666666667 0.50000000000000 0.16666666666667 137.327 > 109 113 Ba 0.50000000000000 0.50000000000000 0.16666666666667 137.327 > 109 114 Ba 0.83333333333333 0.50000000000000 0.16666666666667 137.327 > 109 115 Ba 0.16666666666667 0.83333333333333 0.16666666666667 137.327 > 109 116 Ba 0.50000000000000 0.83333333333333 0.16666666666667 137.327 > 109 117 Ba 0.83333333333333 0.83333333333333 0.16666666666667 137.327 > 109 118 Ba 0.16666666666667 0.16666666666667 0.50000000000000 137.327 > 109 119 Ba 0.50000000000000 0.16666666666667 0.50000000000000 137.327 > 109 120 Ba 0.83333333333333 0.16666666666667 0.50000000000000 137.327 > 109 121 Ba 0.16666666666667 0.50000000000000 0.50000000000000 137.327 > 109 122 Ba 0.50000000000000 0.50000000000000 0.50000000000000 137.327 > 109 123 Ba 0.83333333333333 0.50000000000000 0.50000000000000 137.327 > 109 124 Ba 0.16666666666667 0.83333333333333 0.50000000000000 137.327 > 109 125 Ba 0.50000000000000 0.83333333333333 0.50000000000000 137.327 > 109 126 Ba 0.83333333333333 0.83333333333333 0.50000000000000 137.327 > 109 127 Ba 0.16666666666667 0.16666666666667 0.83333333333333 137.327 > 109 128 Ba 0.50000000000000 0.16666666666667 0.83333333333333 137.327 > 109 129 Ba 0.83333333333333 0.16666666666667 0.83333333333333 137.327 > 109 130 Ba 0.16666666666667 0.50000000000000 0.83333333333333 137.327 > 109 131 Ba 0.50000000000000 0.50000000000000 0.83333333333333 137.327 > 109 132 Ba 0.83333333333333 0.50000000000000 0.83333333333333 137.327 > 109 133 Ba 0.16666666666667 0.83333333333333 0.83333333333333 137.327 > 109 134 Ba 0.50000000000000 0.83333333333333 0.83333333333333 137.327 > 109 135 Ba 0.83333333333333 0.83333333333333 0.83333333333333 137.327 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 5.3490772 0.0000000 0.0000000 0.0000000 5.3490772 0.0000000 0.0000000 0.0000000 5.3490772 -------------------------- Born effective charges -------------------------- 1 H -1.1841729 0.0000000 0.0000000 0.0000000 -0.8761288 0.0000000 0.0000000 0.0000000 -1.1841729 2 H -1.1841729 0.0000000 0.0000000 0.0000000 -1.1841729 0.0000000 0.0000000 0.0000000 -0.8761288 3 H -0.8761288 0.0000000 0.0000000 0.0000000 -1.1841729 0.0000000 0.0000000 0.0000000 -1.1841729 4 Li 0.8715575 0.0000000 0.0000000 0.0000000 0.8715575 0.0000000 0.0000000 0.0000000 0.8715575 5 Ba 2.3729170 0.0000000 0.0000000 0.0000000 2.3729170 0.0000000 0.0000000 0.0000000 2.3729170 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: -0.00000000 (xxx) -0.00000000 (xxx) -0.00000000 (xxx) 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 13 ] 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: 1.05, Number of G-points: 305, Lambda: 0.25 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/84) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 84 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 8.713 ( 0.000 0.000 0.000) 0.000 8.713 ( 0.000 0.000 0.000) 0.000 8.713 ( 0.000 0.000 0.000) 0.000 18.250 ( 0.000 0.000 0.000) 0.000 18.250 ( 0.000 0.000 0.000) 0.000 18.250 ( 0.000 0.000 0.000) 0.000 19.749 ( 0.000 0.000 0.000) 0.000 19.749 ( 0.000 0.000 0.000) 0.000 19.749 ( 0.000 0.000 0.000) 0.000 27.302 ( 0.000 0.000 0.000) 0.000 27.302 ( 0.000 0.000 0.000) 0.000 27.302 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/84) ======================= q-point: ( 0.08 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.639 ( 31.724 0.000 0.000) 31.724 0.639 ( 31.724 0.000 0.000) 31.724 0.935 ( 47.133 0.000 0.000) 47.133 8.750 ( 3.793 0.000 0.000) 3.793 8.750 ( 3.793 0.000 0.000) 3.793 9.119 ( 7.979 0.000 0.000) 7.979 18.345 ( 9.437 0.000 0.000) 9.437 18.345 ( 9.437 0.000 0.000) 9.437 19.829 ( 8.234 0.000 0.000) 8.234 19.829 ( 8.234 0.000 0.000) 8.234 19.951 ( 20.245 0.000 0.000) 20.245 26.796 ( 7.211 0.000 0.000) 7.211 27.290 ( -1.173 0.000 0.000) 1.173 27.290 ( -1.173 0.000 0.000) 1.173 29.148 ( -2.075 0.000 0.000) 2.075 ======================= Grid point 2 (3/84) ======================= q-point: ( 0.15 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.214 ( 26.977 0.000 0.000) 26.977 1.214 ( 26.977 0.000 0.000) 26.977 1.813 ( 42.760 0.000 0.000) 42.760 8.854 ( 6.673 0.000 0.000) 6.673 8.854 ( 6.673 0.000 0.000) 6.673 9.338 ( 14.125 0.000 0.000) 14.125 18.588 ( 14.611 0.000 0.000) 14.611 18.588 ( 14.611 0.000 0.000) 14.611 20.066 ( 15.979 0.000 0.000) 15.979 20.066 ( 15.979 0.000 0.000) 15.979 20.499 ( 34.872 0.000 0.000) 34.872 26.993 ( 12.612 0.000 0.000) 12.612 27.260 ( -1.835 0.000 0.000) 1.835 27.260 ( -1.835 0.000 0.000) 1.835 29.091 ( -3.711 0.000 0.000) 3.711 ======================= Grid point 3 (4/84) ======================= q-point: ( 0.23 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.678 ( 20.660 0.000 0.000) 20.660 1.678 ( 20.660 0.000 0.000) 20.660 2.581 ( 35.860 0.000 0.000) 35.860 9.000 ( 7.993 0.000 0.000) 7.993 9.000 ( 7.993 0.000 0.000) 7.993 9.647 ( 17.042 0.000 0.000) 17.042 18.881 ( 14.931 0.000 0.000) 14.931 18.881 ( 14.931 0.000 0.000) 14.931 20.430 ( 20.671 0.000 0.000) 20.671 20.430 ( 20.671 0.000 0.000) 20.671 21.247 ( 40.549 0.000 0.000) 40.549 27.224 ( -1.673 0.000 0.000) 1.673 27.224 ( -1.673 0.000 0.000) 1.673 27.267 ( 15.044 0.000 0.000) 15.044 29.009 ( -4.632 0.000 0.000) 4.632 ======================= Grid point 4 (5/84) ======================= q-point: ( 0.31 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.015 ( 14.106 0.000 0.000) 14.106 2.015 ( 14.106 0.000 0.000) 14.106 3.194 ( 26.968 0.000 0.000) 26.968 9.153 ( 7.487 0.000 0.000) 7.487 9.153 ( 7.487 0.000 0.000) 7.487 9.975 ( 16.068 0.000 0.000) 16.068 19.151 ( 12.524 0.000 0.000) 12.524 19.151 ( 12.524 0.000 0.000) 12.524 20.831 ( 19.553 0.000 0.000) 19.553 20.831 ( 19.553 0.000 0.000) 19.553 22.013 ( 36.830 0.000 0.000) 36.830 27.199 ( -0.839 0.000 0.000) 0.839 27.199 ( -0.839 0.000 0.000) 0.839 27.556 ( 14.245 0.000 0.000) 14.245 28.916 ( -4.751 0.000 0.000) 4.751 ======================= Grid point 5 (6/84) ======================= q-point: ( 0.38 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.229 ( 8.055 0.000 0.000) 8.055 2.229 ( 8.055 0.000 0.000) 8.055 3.620 ( 16.700 0.000 0.000) 16.700 9.279 ( 5.289 0.000 0.000) 5.289 9.279 ( 5.289 0.000 0.000) 5.289 10.247 ( 11.411 0.000 0.000) 11.411 19.358 ( 8.495 0.000 0.000) 8.495 19.358 ( 8.495 0.000 0.000) 8.495 21.156 ( 13.284 0.000 0.000) 13.284 21.156 ( 13.284 0.000 0.000) 13.284 22.628 ( 25.392 0.000 0.000) 25.392 27.192 ( -0.025 0.000 0.000) 0.025 27.192 ( -0.025 0.000 0.000) 0.025 27.801 ( 10.492 0.000 0.000) 10.492 28.830 ( -3.937 0.000 0.000) 3.937 ======================= Grid point 6 (7/84) ======================= q-point: ( 0.46 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 196 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.332 ( 2.602 0.000 0.000) 2.602 2.332 ( 2.602 0.000 0.000) 2.602 3.838 ( 5.650 0.000 0.000) 5.650 9.351 ( 1.905 0.000 0.000) 1.905 9.351 ( 1.905 0.000 0.000) 1.905 10.401 ( 4.122 0.000 0.000) 4.122 19.472 ( 3.045 0.000 0.000) 3.045 19.472 ( 3.045 0.000 0.000) 3.045 21.332 ( 4.597 0.000 0.000) 4.597 21.332 ( 4.597 0.000 0.000) 4.597 22.968 ( 9.024 0.000 0.000) 9.024 27.194 ( 0.130 0.000 0.000) 0.130 27.194 ( 0.130 0.000 0.000) 0.130 27.946 ( 4.005 0.000 0.000) 4.005 28.773 ( -1.675 0.000 0.000) 1.675 ======================= Grid point 14 (8/84) ======================= q-point: ( 0.08 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.830 ( 20.686 20.686 0.000) 29.255 0.888 ( 21.173 21.173 0.000) 29.943 1.340 ( 32.072 32.072 0.000) 45.356 8.789 ( 3.859 3.859 0.000) 5.457 8.861 ( 7.349 7.349 0.000) 10.392 9.111 ( 3.705 3.705 0.000) 5.240 18.387 ( 6.767 6.767 0.000) 9.570 18.449 ( 10.184 10.184 0.000) 14.403 19.803 ( 2.691 2.691 0.000) 3.806 19.905 ( 7.841 7.841 0.000) 11.089 20.232 ( 23.133 23.133 0.000) 32.715 26.941 ( 10.392 10.392 0.000) 14.697 27.192 ( -5.650 -5.650 0.000) 7.990 27.277 ( -1.301 -1.301 0.000) 1.839 29.168 ( -0.118 -0.118 0.000) 0.167 ======================= Grid point 15 (9/84) ======================= q-point: ( 0.15 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.301 ( 24.953 9.816 0.000) 26.815 1.331 ( 22.702 11.118 0.000) 25.279 2.008 ( 34.997 17.741 0.000) 39.237 8.894 ( 6.753 4.003 0.000) 7.851 8.975 ( 5.533 11.609 0.000) 12.860 9.312 ( 15.152 -1.593 0.000) 15.236 18.605 ( 14.458 2.946 0.000) 14.755 18.711 ( 15.823 12.199 0.000) 19.980 19.958 ( 13.030 -5.362 0.000) 14.090 20.132 ( 15.406 6.821 0.000) 16.848 20.790 ( 33.284 24.282 0.000) 41.200 27.038 ( -2.153 -2.321 0.000) 3.166 27.223 ( 10.192 3.087 0.000) 10.649 27.243 ( -2.077 -1.672 0.000) 2.667 29.132 ( -3.131 2.886 0.000) 4.259 ======================= Grid point 16 (10/84) ======================= q-point: ( 0.23 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.734 ( 18.361 5.485 0.000) 19.162 1.744 ( 20.047 6.832 0.000) 21.179 2.659 ( 31.269 7.517 0.000) 32.159 9.041 ( 8.032 4.118 0.000) 9.026 9.100 ( 7.174 9.974 0.000) 12.286 9.643 ( 18.048 -0.084 0.000) 18.048 18.905 ( 15.434 2.604 0.000) 15.652 19.028 ( 16.064 14.641 0.000) 21.735 20.291 ( 20.526 -10.572 0.000) 23.089 20.487 ( 20.332 5.813 0.000) 21.146 21.478 ( 36.241 20.388 0.000) 41.582 26.983 ( -2.530 -15.064 0.000) 15.275 27.202 ( -1.996 -2.242 0.000) 3.002 27.446 ( 11.471 8.666 0.000) 14.376 29.062 ( -3.779 4.917 0.000) 6.202 ======================= Grid point 17 (11/84) ======================= q-point: ( 0.31 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.037 ( 12.760 2.138 0.000) 12.938 2.072 ( 13.680 5.712 0.000) 14.825 3.202 ( 24.228 0.803 0.000) 24.242 9.194 ( 7.468 4.131 0.000) 8.534 9.243 ( 7.208 9.035 0.000) 11.558 9.987 ( 16.656 1.246 0.000) 16.702 19.182 ( 12.735 3.094 0.000) 13.106 19.316 ( 13.289 16.564 0.000) 21.236 20.714 ( 21.725 -10.685 0.000) 24.210 20.885 ( 19.560 5.458 0.000) 20.307 22.150 ( 31.785 12.828 0.000) 34.276 26.950 ( -0.852 -19.706 0.000) 19.724 27.170 ( -1.175 -2.911 0.000) 3.139 27.642 ( 8.233 2.897 0.000) 8.728 28.997 ( -2.642 8.280 0.000) 8.692 ======================= Grid point 18 (12/84) ======================= q-point: ( 0.38 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.232 ( 7.334 0.194 0.000) 7.336 2.279 ( 7.745 4.974 0.000) 9.205 3.588 ( 15.222 -3.181 0.000) 15.551 9.320 ( 5.246 4.062 0.000) 6.635 9.367 ( 5.255 8.776 0.000) 10.229 10.267 ( 11.712 2.027 0.000) 11.886 19.390 ( 8.457 3.133 0.000) 9.018 19.534 ( 8.918 17.739 0.000) 19.854 21.089 ( 15.829 -6.809 0.000) 17.231 21.212 ( 13.420 5.621 0.000) 14.550 22.676 ( 21.528 4.437 0.000) 21.980 26.949 ( 0.558 -21.243 0.000) 21.251 27.157 ( -0.280 -3.516 0.000) 3.527 27.758 ( 3.782 -7.018 0.000) 7.973 28.963 ( -0.868 13.136 0.000) 13.165 ======================= Grid point 19 (13/84) ======================= q-point: ( 0.46 0.08 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.325 ( 2.374 -0.706 0.000) 2.477 2.377 ( 2.486 4.553 0.000) 5.188 3.787 ( 5.181 -5.033 0.000) 7.223 9.391 ( 1.884 3.994 0.000) 4.415 9.438 ( 1.915 8.762 0.000) 8.969 10.425 ( 4.215 2.387 0.000) 4.844 19.503 ( 2.998 2.965 0.000) 4.216 19.654 ( 3.181 18.289 0.000) 18.563 21.302 ( 5.672 -3.153 0.000) 6.490 21.390 ( 4.665 5.839 0.000) 7.474 22.963 ( 7.608 -1.085 0.000) 7.685 26.962 ( 0.508 -21.235 0.000) 21.241 27.156 ( 0.035 -3.876 0.000) 3.876 27.799 ( 0.829 -14.577 0.000) 14.601 28.957 ( -0.037 16.716 0.000) 16.716 ======================= Grid point 28 (14/84) ======================= q-point: ( 0.15 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.583 ( 17.790 17.790 0.000) 25.159 1.594 ( 14.766 14.766 0.000) 20.882 2.413 ( 22.393 22.393 0.000) 31.669 9.003 ( 6.923 6.923 0.000) 9.791 9.246 ( 11.821 11.821 0.000) 16.718 9.335 ( 7.907 7.907 0.000) 11.182 18.737 ( 10.649 10.649 0.000) 15.060 19.028 ( 19.335 19.335 0.000) 27.344 19.944 ( 4.308 4.308 0.000) 6.093 20.332 ( 13.682 13.682 0.000) 19.349 21.339 ( 30.861 30.861 0.000) 43.645 26.868 ( -10.887 -10.887 0.000) 15.396 27.198 ( -2.784 -2.784 0.000) 3.937 27.402 ( 10.652 10.652 0.000) 15.064 29.168 ( 0.794 0.794 0.000) 1.123 ======================= Grid point 29 (15/84) ======================= q-point: ( 0.23 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.870 ( 13.078 7.945 0.000) 15.302 1.935 ( 17.289 12.342 0.000) 21.242 2.842 ( 21.209 10.580 0.000) 23.702 9.152 ( 8.099 7.038 0.000) 10.730 9.366 ( 5.675 16.817 0.000) 17.748 9.659 ( 19.423 2.041 0.000) 19.530 18.993 ( 14.675 6.792 0.000) 16.171 19.416 ( 19.564 24.129 0.000) 31.064 20.125 ( 14.179 -5.212 0.000) 15.106 20.655 ( 19.046 11.456 0.000) 22.226 21.940 ( 29.750 25.505 0.000) 39.186 26.672 ( -8.846 -16.477 0.000) 18.702 27.141 ( -2.942 -3.876 0.000) 4.866 27.562 ( 5.025 2.501 0.000) 5.613 29.181 ( 1.084 7.144 0.000) 7.226 ======================= Grid point 30 (16/84) ======================= q-point: ( 0.31 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.090 ( 9.418 3.118 0.000) 9.921 2.227 ( 12.369 9.825 0.000) 15.797 3.221 ( 17.358 0.982 0.000) 17.386 9.305 ( 7.405 7.013 0.000) 10.199 9.489 ( 6.568 15.770 0.000) 17.083 10.027 ( 17.715 2.942 0.000) 17.958 19.266 ( 12.647 5.492 0.000) 13.788 19.762 ( 15.630 28.156 0.000) 32.204 20.480 ( 21.433 -11.972 0.000) 24.550 21.039 ( 19.375 10.244 0.000) 21.916 22.450 ( 21.764 16.375 0.000) 27.236 26.535 ( -5.054 -22.032 0.000) 22.605 27.090 ( -2.167 -5.176 0.000) 5.611 27.590 ( -1.810 -7.859 0.000) 8.065 29.221 ( 3.018 13.935 0.000) 14.258 ======================= Grid point 31 (17/84) ======================= q-point: ( 0.38 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.235 ( 5.484 0.106 0.000) 5.485 2.413 ( 6.931 8.533 0.000) 10.994 3.503 ( 11.387 -5.323 0.000) 12.570 9.429 ( 5.132 6.888 0.000) 8.590 9.606 ( 5.100 15.313 0.000) 16.140 10.324 ( 12.361 3.725 0.000) 12.910 19.468 ( 7.978 4.658 0.000) 9.239 20.014 ( 10.163 30.531 0.000) 32.178 20.905 ( 20.611 -11.535 0.000) 23.620 21.369 ( 13.746 10.282 0.000) 17.166 22.765 ( 10.642 3.366 0.000) 11.161 26.474 ( -1.416 -26.094 0.000) 26.133 27.059 ( -1.042 -6.369 0.000) 6.454 27.519 ( -4.585 -15.785 0.000) 16.437 29.286 ( 3.201 18.443 0.000) 18.719 ======================= Grid point 32 (18/84) ======================= q-point: ( 0.46 0.15 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.305 ( 1.782 -1.332 0.000) 2.225 2.501 ( 2.189 7.859 0.000) 8.158 3.654 ( 3.961 -8.317 0.000) 9.212 9.498 ( 1.829 6.778 0.000) 7.021 9.676 ( 1.903 15.178 0.000) 15.297 10.490 ( 4.435 4.181 0.000) 6.095 19.573 ( 2.720 3.879 0.000) 4.738 20.150 ( 3.574 31.586 0.000) 31.787 21.209 ( 8.870 -6.426 0.000) 10.953 21.552 ( 4.846 10.621 0.000) 11.674 22.883 ( 2.398 -8.657 0.000) 8.982 26.465 ( 0.037 -28.221 0.000) 28.221 27.047 ( -0.250 -7.087 0.000) 7.092 27.446 ( -2.278 -19.261 0.000) 19.395 29.332 ( 1.305 20.280 0.000) 20.322 ======================= Grid point 42 (19/84) ======================= q-point: ( 0.23 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.023 ( 7.432 7.432 0.000) 10.511 2.196 ( 13.617 13.617 0.000) 19.257 3.049 ( 10.352 10.352 0.000) 14.640 9.302 ( 8.101 8.101 0.000) 11.457 9.719 ( 11.543 11.543 0.000) 16.325 9.733 ( 12.756 12.756 0.000) 18.039 19.169 ( 10.957 10.957 0.000) 15.495 19.911 ( 25.333 25.333 0.000) 35.827 20.111 ( 4.076 4.076 0.000) 5.764 20.929 ( 16.464 16.464 0.000) 23.284 22.411 ( 21.494 21.494 0.000) 30.397 26.364 ( -14.728 -14.728 0.000) 20.828 27.058 ( -4.426 -4.426 0.000) 6.259 27.531 ( -5.400 -5.400 0.000) 7.637 29.336 ( 8.608 8.608 0.000) 12.174 ======================= Grid point 43 (20/84) ======================= q-point: ( 0.31 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.150 ( 5.467 2.858 0.000) 6.169 2.436 ( 10.582 11.192 0.000) 15.403 3.238 ( 8.989 0.789 0.000) 9.023 9.454 ( 7.283 7.995 0.000) 10.815 9.832 ( 5.862 18.856 0.000) 19.747 10.100 ( 18.150 4.556 0.000) 18.713 19.390 ( 10.766 7.058 0.000) 12.873 20.294 ( 14.870 -6.182 0.000) 16.104 20.356 ( 19.791 31.464 0.000) 37.171 21.275 ( 18.344 13.859 0.000) 22.991 22.719 ( 9.012 9.530 0.000) 13.116 26.117 ( -10.314 -20.528 0.000) 22.973 26.977 ( -3.759 -6.235 0.000) 7.280 27.384 ( -8.692 -11.824 0.000) 14.675 29.502 ( 8.044 14.117 0.000) 16.248 ======================= Grid point 44 (21/84) ======================= q-point: ( 0.38 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.234 ( 3.191 -0.211 0.000) 3.198 2.597 ( 5.933 9.943 0.000) 11.578 3.391 ( 6.496 -5.964 0.000) 8.819 9.575 ( 4.980 7.831 0.000) 9.280 9.938 ( 4.720 18.232 0.000) 18.833 10.407 ( 12.839 4.685 0.000) 13.667 19.559 ( 6.420 4.454 0.000) 7.814 20.668 ( 12.230 35.281 0.000) 37.341 20.669 ( 22.234 -11.820 0.000) 25.181 21.600 ( 13.988 13.117 0.000) 19.176 22.744 ( -5.871 -7.101 0.000) 9.214 25.968 ( -5.110 -25.400 0.000) 25.909 26.918 ( -2.294 -7.948 0.000) 8.272 27.231 ( -6.303 -11.823 0.000) 13.399 29.635 ( 5.402 16.732 0.000) 17.583 ======================= Grid point 45 (22/84) ======================= q-point: ( 0.46 0.23 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.275 ( 1.034 -1.730 0.000) 2.015 2.671 ( 1.819 9.269 0.000) 9.446 3.480 ( 2.409 -9.215 0.000) 9.525 9.642 ( 1.762 7.710 0.000) 7.908 10.004 ( 1.794 17.883 0.000) 17.973 10.580 ( 4.637 4.897 0.000) 6.744 19.640 ( 2.029 2.848 0.000) 3.497 20.828 ( 4.188 36.835 0.000) 37.072 21.055 ( 13.508 -9.209 0.000) 16.348 21.789 ( 5.083 13.370 0.000) 14.303 22.573 ( -7.939 -23.867 0.000) 25.153 25.907 ( -1.416 -28.352 0.000) 28.388 26.889 ( -0.726 -9.008 0.000) 9.037 27.151 ( -2.000 -9.271 0.000) 9.485 29.706 ( 1.860 17.629 0.000) 17.727 ======================= Grid point 56 (23/84) ======================= q-point: ( 0.31 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.198 ( 1.988 1.988 0.000) 2.812 2.640 ( 9.349 9.349 0.000) 13.222 3.254 ( 1.039 1.039 0.000) 1.469 9.604 ( 7.113 7.113 0.000) 10.060 10.188 ( 11.852 11.852 0.000) 16.761 10.204 ( 11.093 11.093 0.000) 15.688 19.529 ( 6.788 6.788 0.000) 9.599 20.272 ( 4.593 4.593 0.000) 6.496 20.926 ( 25.870 25.870 0.000) 36.586 21.567 ( 15.803 15.803 0.000) 22.349 22.756 ( -6.870 -6.870 0.000) 9.715 25.758 ( -15.911 -15.911 0.000) 22.502 26.858 ( -5.773 -5.773 0.000) 8.164 27.186 ( -7.496 -7.496 0.000) 10.601 29.748 ( 10.927 10.927 0.000) 15.454 ======================= Grid point 57 (24/84) ======================= q-point: ( 0.38 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.227 ( 1.075 -0.479 0.000) 1.177 2.784 ( 5.264 8.959 0.000) 10.391 3.281 ( 1.643 -5.084 0.000) 5.343 9.722 ( 4.822 6.941 0.000) 8.452 10.288 ( 4.078 17.095 0.000) 17.575 10.499 ( 12.773 4.590 0.000) 13.573 19.634 ( 3.816 3.135 0.000) 4.938 20.487 ( 17.157 -5.837 0.000) 18.122 21.329 ( 15.340 31.337 0.000) 34.890 21.863 ( 13.564 13.533 0.000) 19.161 22.454 ( -22.816 -22.634 0.000) 32.139 25.512 ( -9.306 -20.900 0.000) 22.878 26.762 ( -3.945 -7.776 0.000) 8.719 27.103 ( -1.452 -1.367 0.000) 1.994 29.919 ( 6.576 12.305 0.000) 13.952 ======================= Grid point 58 (25/84) ======================= q-point: ( 0.46 0.31 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.240 ( 0.329 -1.738 0.000) 1.769 2.848 ( 1.484 8.626 0.000) 8.753 3.309 ( 0.924 -8.038 0.000) 8.091 9.786 ( 1.698 6.831 0.000) 7.039 10.345 ( 1.555 16.553 0.000) 16.626 10.672 ( 4.653 4.353 0.000) 6.372 19.680 ( 1.090 1.296 0.000) 1.693 20.864 ( 17.063 -9.942 0.000) 19.749 21.525 ( 4.996 33.521 0.000) 33.891 21.985 ( -19.945 -34.652 0.000) 39.982 22.053 ( 5.286 13.294 0.000) 14.306 25.393 ( -2.988 -23.826 0.000) 24.013 26.709 ( -1.372 -9.071 0.000) 9.175 27.100 ( 0.371 3.147 0.000) 3.169 30.004 ( 2.156 12.832 0.000) 13.012 ======================= Grid point 70 (26/84) ======================= q-point: ( 0.38 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.217 ( -0.479 -0.479 0.000) 0.678 2.929 ( 5.551 5.551 0.000) 7.850 3.203 ( -2.727 -2.727 0.000) 3.857 9.836 ( 4.688 4.688 0.000) 6.630 10.573 ( 7.632 7.632 0.000) 10.793 10.584 ( 7.909 7.909 0.000) 11.184 19.678 ( 1.430 1.430 0.000) 2.023 20.479 ( 5.562 5.562 0.000) 7.866 21.837 ( 19.902 19.902 0.000) 28.145 21.909 ( -31.543 -31.543 0.000) 44.609 22.110 ( 11.547 11.547 0.000) 16.330 25.177 ( -13.180 -13.180 0.000) 18.640 26.628 ( -5.687 -5.687 0.000) 8.042 27.143 ( 4.080 4.080 0.000) 5.769 30.109 ( 7.238 7.238 0.000) 10.237 ======================= Grid point 71 (27/84) ======================= q-point: ( 0.46 0.38 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.210 ( -0.199 -1.291 0.000) 1.306 2.993 ( 1.312 6.072 0.000) 6.212 3.177 ( -0.278 -5.378 0.000) 5.385 9.899 ( 1.647 4.610 0.000) 4.895 10.623 ( 1.244 11.622 0.000) 11.688 10.743 ( 4.416 2.844 0.000) 5.253 19.694 ( 0.349 0.289 0.000) 0.453 20.695 ( 14.740 -6.528 0.000) 16.121 21.323 ( -25.128 -31.171 0.000) 40.039 22.085 ( 5.966 23.089 0.000) 23.847 22.285 ( 5.310 9.988 0.000) 11.312 25.006 ( -4.425 -15.581 0.000) 16.198 26.551 ( -2.084 -6.881 0.000) 7.190 27.214 ( 2.295 6.919 0.000) 7.289 30.202 ( 2.354 7.557 0.000) 7.915 ======================= Grid point 84 (28/84) ======================= q-point: ( 0.46 0.46 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 343 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.193 ( -0.477 -0.477 0.000) 0.675 3.073 ( 1.871 1.871 0.000) 2.646 3.107 ( -1.569 -1.569 0.000) 2.219 9.960 ( 1.618 1.618 0.000) 2.289 10.778 ( 2.750 2.750 0.000) 3.890 10.782 ( 2.385 2.385 0.000) 3.373 19.697 ( 0.012 0.012 0.000) 0.017 20.647 ( 2.500 2.500 0.000) 3.536 20.843 ( -17.129 -17.129 0.000) 24.224 22.391 ( 7.681 7.681 0.000) 10.863 22.425 ( 4.258 4.258 0.000) 6.022 24.800 ( -5.368 -5.368 0.000) 7.591 26.456 ( -2.595 -2.595 0.000) 3.670 27.320 ( 3.162 3.162 0.000) 4.471 30.299 ( 2.467 2.467 0.000) 3.489 ======================= Grid point 183 (29/84) ======================= q-point: ( 0.08 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.031 ( 16.570 16.570 16.570) 28.700 1.031 ( 16.570 16.570 16.570) 28.700 1.618 ( 24.355 24.355 24.355) 42.184 8.894 ( 5.843 5.843 5.843) 10.120 8.894 ( 5.843 5.843 5.843) 10.120 9.112 ( 2.896 2.896 2.896) 5.016 18.482 ( 7.569 7.569 7.569) 13.110 18.482 ( 7.569 7.569 7.569) 13.110 19.639 ( -4.188 -4.188 -4.188) 7.254 20.305 ( 17.698 17.698 17.698) 30.654 20.305 ( 17.698 17.698 17.698) 30.654 27.076 ( 10.762 10.762 10.762) 18.640 27.180 ( -4.184 -4.184 -4.184) 7.247 27.180 ( -4.184 -4.184 -4.184) 7.247 29.192 ( 0.813 0.813 0.813) 1.409 ======================= Grid point 184 (30/84) ======================= q-point: ( 0.15 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.407 ( 21.268 10.184 10.184) 25.685 1.428 ( 21.116 9.966 9.966) 25.388 2.158 ( 29.438 14.057 14.057) 35.522 8.984 ( 5.840 6.497 6.497) 10.887 9.028 ( 5.993 7.606 7.606) 12.313 9.304 ( 15.545 -0.177 -0.177) 15.547 18.703 ( 16.209 5.677 5.677) 18.088 18.716 ( 13.499 6.353 6.353) 16.216 19.688 ( 9.546 -12.841 -12.841) 20.516 20.504 ( 12.982 21.721 21.721) 33.349 20.927 ( 34.811 15.283 15.283) 40.975 27.019 ( -7.161 -4.087 -4.087) 9.202 27.148 ( -1.963 -5.746 -5.746) 8.359 27.326 ( 12.423 8.561 8.561) 17.347 29.175 ( -2.298 3.222 3.222) 5.103 ======================= Grid point 185 (31/84) ======================= q-point: ( 0.23 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.795 ( 17.800 6.187 6.187) 19.834 1.811 ( 17.759 6.496 6.496) 19.994 2.719 ( 27.532 5.841 5.841) 28.745 9.112 ( 7.047 5.662 5.662) 10.667 9.161 ( 7.743 7.699 7.699) 13.361 9.646 ( 18.680 0.547 0.547) 18.696 18.961 ( 11.043 3.597 3.597) 12.158 19.069 ( 20.391 7.056 7.056) 22.702 20.004 ( 21.961 -16.714 -16.714) 32.264 20.796 ( 16.537 18.927 18.927) 31.463 21.615 ( 34.870 13.910 13.910) 40.036 26.899 ( -5.099 -9.342 -9.342) 14.161 27.108 ( -1.974 -5.943 -5.943) 8.633 27.544 ( 9.398 6.828 6.828) 13.475 29.121 ( -2.664 5.569 5.569) 8.314 ======================= Grid point 186 (32/84) ======================= q-point: ( 0.31 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.085 ( 12.008 3.713 3.713) 13.106 2.107 ( 12.615 4.586 4.586) 14.185 3.205 ( 21.939 0.218 0.218) 21.941 9.248 ( 6.632 4.811 4.811) 9.501 9.317 ( 7.873 8.021 8.021) 13.808 10.000 ( 17.100 1.452 1.452) 17.223 19.144 ( 7.995 -1.141 -1.141) 8.156 19.459 ( 18.886 11.938 11.938) 25.332 20.478 ( 24.980 -15.079 -15.079) 32.844 21.121 ( 16.119 15.222 15.222) 26.893 22.244 ( 29.027 9.382 9.382) 31.916 26.826 ( -2.217 -13.225 -13.225) 18.833 27.075 ( -1.375 -6.336 -6.336) 9.066 27.673 ( 3.589 0.452 0.452) 3.646 29.087 ( -0.629 8.962 8.962) 12.689 ======================= Grid point 187 (33/84) ======================= q-point: ( 0.38 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.265 ( 6.649 1.945 1.945) 7.196 2.301 ( 7.350 3.594 3.594) 8.936 3.557 ( 13.995 -3.117 -3.117) 14.672 9.360 ( 4.696 4.123 4.123) 7.486 9.452 ( 5.746 8.489 8.489) 13.309 10.288 ( 11.967 2.074 2.074) 12.321 19.274 ( 5.271 -4.887 -4.887) 8.693 19.777 ( 13.298 17.606 17.606) 28.227 20.912 ( 18.461 -11.466 -11.466) 24.571 21.393 ( 11.333 11.964 11.964) 20.364 22.717 ( 19.107 3.883 3.883) 19.880 26.812 ( 0.542 -15.167 -15.167) 21.457 27.055 ( -0.650 -6.894 -6.894) 9.771 27.689 ( -1.348 -8.119 -8.119) 11.561 29.095 ( 1.154 12.707 12.707) 18.008 ======================= Grid point 188 (34/84) ======================= q-point: ( 0.46 0.08 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.349 ( 2.103 0.983 0.983) 2.521 2.395 ( 2.396 3.175 3.175) 5.090 3.740 ( 4.793 -4.655 -4.655) 8.143 9.423 ( 1.693 3.739 3.739) 5.551 9.530 ( 2.091 8.795 8.795) 12.613 10.449 ( 4.298 2.387 2.387) 5.465 19.344 ( 1.890 -6.829 -6.829) 9.841 19.957 ( 4.825 21.310 21.310) 30.521 21.161 ( 6.657 -8.674 -8.674) 13.957 21.544 ( 3.980 10.107 10.107) 14.838 22.970 ( 6.685 0.153 0.153) 6.689 26.831 ( 0.893 -15.258 -15.258) 21.596 27.048 ( -0.165 -7.344 -7.344) 10.388 27.653 ( -1.510 -14.379 -14.379) 20.391 29.117 ( 0.750 14.863 14.863) 21.032 ======================= Grid point 197 (35/84) ======================= q-point: ( 0.15 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.673 ( 16.135 16.135 8.709) 24.424 1.677 ( 14.348 14.348 8.204) 21.887 2.492 ( 19.001 19.001 7.615) 27.929 9.096 ( 6.333 6.333 9.257) 12.881 9.268 ( 11.597 11.597 2.315) 16.563 9.352 ( 8.716 8.716 2.131) 12.509 18.816 ( 9.605 9.605 8.116) 15.823 18.995 ( 17.527 17.527 -1.662) 24.843 19.552 ( -0.064 -0.064 -22.796) 22.796 20.850 ( 15.032 15.032 35.276) 41.186 21.394 ( 30.091 30.091 5.504) 42.910 26.855 ( -10.912 -10.912 -1.321) 15.488 27.068 ( -3.088 -3.088 -10.998) 11.833 27.516 ( 9.250 9.250 8.903) 15.824 29.223 ( 1.671 1.671 4.798) 5.348 ======================= Grid point 198 (36/84) ======================= q-point: ( 0.23 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.944 ( 12.612 8.292 7.264) 16.751 1.993 ( 15.672 11.861 5.639) 20.447 2.865 ( 18.904 8.639 2.255) 20.907 9.231 ( 7.302 6.980 7.870) 12.805 9.405 ( 6.664 16.134 4.042) 17.918 9.676 ( 19.409 2.754 1.843) 19.690 19.028 ( 9.852 4.670 2.634) 11.217 19.338 ( 18.043 19.124 -6.188) 27.010 19.744 ( 19.814 -8.947 -22.991) 31.642 21.136 ( 14.731 16.065 33.288) 39.789 21.989 ( 29.334 23.015 5.140) 37.638 26.647 ( -9.870 -15.134 -2.678) 18.265 27.010 ( -2.901 -4.651 -11.476) 12.718 27.636 ( 2.341 1.523 5.428) 6.104 29.256 ( 2.232 7.945 7.097) 10.885 ======================= Grid point 199 (37/84) ======================= q-point: ( 0.31 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.153 ( 8.781 3.586 6.193) 11.327 2.257 ( 11.256 9.838 3.054) 15.258 3.209 ( 16.081 0.233 -1.132) 16.122 9.367 ( 6.430 7.319 6.031) 11.458 9.547 ( 7.586 15.056 6.065) 17.917 10.045 ( 17.869 3.229 1.940) 18.262 19.164 ( 4.936 3.011 -10.513) 11.998 19.737 ( 22.296 15.445 -0.312) 27.125 20.218 ( 26.292 -9.576 -17.617) 33.066 21.413 ( 13.287 14.652 27.578) 33.938 22.488 ( 21.230 14.488 3.978) 26.009 26.488 ( -6.148 -19.828 -5.036) 21.361 26.959 ( -2.295 -6.122 -11.579) 13.298 27.600 ( -5.671 -7.994 -0.547) 9.816 29.320 ( 4.206 14.063 9.547) 17.510 ======================= Grid point 200 (38/84) ======================= q-point: ( 0.38 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.285 ( 4.847 0.625 4.822) 6.866 2.428 ( 6.421 8.562 1.637) 10.827 3.474 ( 10.773 -5.172 -2.908) 12.299 9.474 ( 4.376 7.294 4.490) 9.619 9.681 ( 5.832 14.599 7.593) 17.458 10.346 ( 12.502 3.805 2.206) 13.253 19.241 ( 3.112 1.223 -17.898) 18.207 20.147 ( 18.380 19.160 9.499) 28.199 20.694 ( 21.261 -9.876 -15.786) 28.263 21.633 ( 8.944 12.550 20.254) 25.450 22.800 ( 10.885 3.527 3.274) 11.902 26.414 ( -1.623 -23.670 -6.581) 24.622 26.924 ( -1.221 -7.236 -12.037) 14.098 27.448 ( -8.824 -15.257 -7.733) 19.247 29.404 ( 4.009 17.667 11.366) 21.387 ======================= Grid point 201 (39/84) ======================= q-point: ( 0.46 0.15 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.346 ( 1.512 -0.794 3.822) 4.187 2.510 ( 2.066 7.834 1.158) 8.185 3.617 ( 3.777 -7.700 -3.663) 9.326 9.532 ( 1.555 7.172 3.678) 8.209 9.761 ( 2.158 14.492 8.348) 16.863 10.514 ( 4.489 4.174 2.386) 6.578 19.283 ( 1.147 0.205 -20.889) 20.922 20.404 ( 7.114 23.810 17.603) 30.453 20.993 ( 8.377 -8.571 -15.969) 19.966 21.751 ( 3.050 11.060 14.239) 18.285 22.927 ( 2.824 -6.066 3.729) 7.660 26.408 ( 0.378 -25.910 -6.288) 26.665 26.911 ( -0.187 -7.461 -12.505) 14.563 27.306 ( -4.576 -18.927 -13.732) 23.827 29.461 ( 1.583 19.009 12.251) 22.670 ======================= Grid point 211 (40/84) ======================= q-point: ( 0.23 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.099 ( 7.322 7.322 7.439) 12.750 2.241 ( 12.858 12.858 4.536) 18.741 3.040 ( 9.136 9.136 -0.883) 12.950 9.378 ( 7.793 7.793 7.622) 13.399 9.751 ( 12.700 12.700 1.750) 18.046 9.761 ( 11.983 11.983 4.386) 17.505 19.145 ( 6.431 6.431 -5.463) 10.609 19.665 ( 6.119 6.119 -26.956) 28.311 19.800 ( 23.092 23.092 -9.976) 34.147 21.430 ( 14.304 14.304 35.599) 40.945 22.434 ( 20.884 20.884 2.237) 29.620 26.350 ( -14.710 -14.710 -1.377) 20.848 26.919 ( -4.738 -4.738 -12.793) 14.442 27.582 ( -6.992 -6.992 3.263) 10.412 29.424 ( 9.194 9.194 8.533) 15.552 ======================= Grid point 212 (41/84) ======================= q-point: ( 0.31 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.220 ( 4.994 3.068 6.894) 9.049 2.466 ( 9.746 11.037 2.945) 15.016 3.215 ( 8.658 0.381 -2.334) 8.975 9.521 ( 6.680 8.205 6.673) 12.509 9.877 ( 6.646 18.235 4.577) 19.941 10.125 ( 18.040 4.907 2.611) 18.877 19.244 ( 3.668 4.585 -18.588) 19.493 19.903 ( 18.835 0.761 -18.308) 26.278 20.288 ( 24.704 17.944 -12.629) 33.042 21.687 ( 11.859 13.512 31.102) 35.924 22.742 ( 9.799 9.839 2.079) 14.041 26.097 ( -10.796 -19.684 -1.999) 22.539 26.831 ( -4.063 -6.870 -13.287) 15.500 27.391 ( -11.708 -12.402 -0.969) 17.083 29.600 ( 8.489 14.067 9.547) 19.002 ======================= Grid point 213 (42/84) ======================= q-point: ( 0.38 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.293 ( 2.635 0.209 5.839) 6.410 2.613 ( 5.415 9.968 1.618) 11.459 3.366 ( 6.548 -5.686 -2.498) 9.025 9.630 ( 4.384 8.444 5.444) 10.961 9.998 ( 5.358 17.323 6.050) 19.115 10.431 ( 12.813 4.803 2.423) 13.897 19.299 ( 2.222 4.104 -25.559) 25.981 20.340 ( 24.294 -3.711 -13.867) 28.218 20.708 ( 17.366 15.792 -9.904) 25.477 21.877 ( 7.400 12.440 23.025) 27.196 22.802 ( -3.047 -4.528 5.049) 7.435 25.942 ( -5.226 -24.125 -2.795) 24.843 26.769 ( -2.131 -8.448 -13.688) 16.226 27.166 ( -10.584 -12.253 -7.553) 17.866 29.741 ( 5.740 16.274 10.268) 20.080 ======================= Grid point 214 (43/84) ======================= q-point: ( 0.46 0.23 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.325 ( 0.788 -1.205 4.995) 5.199 2.680 ( 1.675 9.293 1.036) 9.500 3.456 ( 2.453 -8.535 -2.376) 9.193 9.688 ( 1.524 8.506 4.646) 9.811 10.072 ( 2.030 16.868 6.892) 18.335 10.604 ( 4.635 4.873 2.376) 7.132 19.330 ( 0.836 3.897 -27.774) 28.059 20.762 ( 16.545 -0.842 -12.699) 20.874 20.912 ( 2.233 12.720 -6.697) 14.548 21.972 ( 2.427 11.608 14.571) 18.787 22.695 ( -5.061 -17.559 10.638) 21.145 25.883 ( -1.231 -27.129 -2.610) 27.282 26.751 ( -0.025 -8.871 -13.292) 15.980 27.009 ( -4.724 -9.975 -14.110) 17.914 29.817 ( 1.994 16.988 10.681) 20.166 ======================= Grid point 225 (44/84) ======================= q-point: ( 0.31 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.268 ( 1.854 1.854 6.973) 7.450 2.665 ( 8.947 8.947 2.424) 12.883 3.229 ( 1.218 1.218 -2.572) 3.096 9.673 ( 7.093 7.093 6.900) 12.175 10.221 ( 11.130 11.130 1.779) 15.840 10.238 ( 11.786 11.786 5.084) 17.426 19.324 ( 3.461 3.461 -25.309) 25.778 19.973 ( 8.513 8.513 -16.415) 20.356 20.702 ( 21.694 21.694 -20.908) 37.127 21.931 ( 11.243 11.243 27.632) 31.880 22.810 ( -3.721 -3.721 4.417) 6.870 25.747 ( -15.750 -15.750 -1.063) 22.299 26.703 ( -5.956 -5.956 -14.125) 16.446 27.162 ( -10.339 -10.339 -4.190) 15.210 29.846 ( 10.979 10.979 9.561) 18.235 ======================= Grid point 226 (45/84) ======================= q-point: ( 0.38 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.293 ( 0.771 -0.234 6.506) 6.555 2.800 ( 4.832 8.866 1.555) 10.216 3.263 ( 2.110 -4.669 -1.826) 5.439 9.788 ( 4.627 7.510 6.638) 11.039 10.330 ( 4.522 16.181 4.206) 17.319 10.526 ( 12.467 4.814 2.821) 13.659 19.380 ( 2.369 3.859 -27.399) 27.770 20.239 ( 18.111 -3.291 -16.303) 24.589 21.057 ( 13.107 15.664 -25.710) 32.835 22.110 ( 6.812 11.177 20.482) 24.307 22.611 ( -14.872 -13.938 11.868) 23.586 25.503 ( -9.228 -20.333 -0.947) 22.349 26.613 ( -3.136 -7.262 -14.128) 16.192 27.004 ( -5.884 -4.388 -10.779) 13.041 30.019 ( 6.700 12.135 9.737) 16.940 ======================= Grid point 227 (46/84) ======================= q-point: ( 0.46 0.31 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.301 ( 0.174 -1.251 6.015) 6.146 2.857 ( 1.336 8.578 0.984) 8.737 3.298 ( 1.125 -7.441 -1.112) 7.607 9.849 ( 1.588 7.738 6.285) 10.095 10.392 ( 1.738 15.533 4.955) 16.397 10.695 ( 4.551 4.313 2.343) 6.694 19.413 ( 0.942 4.302 -26.701) 27.062 20.608 ( 16.037 -9.149 -21.504) 28.343 21.155 ( -2.462 2.635 -25.459) 25.713 22.194 ( 2.197 10.949 12.619) 16.851 22.353 ( -8.461 -12.093 22.253) 26.702 25.387 ( -2.883 -23.196 -0.746) 23.387 26.578 ( -0.724 -8.671 -12.915) 15.573 26.929 ( -1.888 0.927 -17.272) 17.399 30.106 ( 2.221 12.549 9.915) 16.147 ======================= Grid point 239 (47/84) ======================= q-point: ( 0.38 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.286 ( -0.413 -0.413 6.793) 6.818 2.941 ( 5.330 5.330 1.260) 7.642 3.194 ( -2.239 -2.239 -0.929) 3.300 9.912 ( 5.014 5.014 7.543) 10.352 10.591 ( 7.656 7.656 1.958) 11.003 10.622 ( 7.475 7.475 4.000) 11.303 19.445 ( 2.870 2.870 -25.193) 25.518 20.259 ( 5.729 5.729 -16.729) 18.588 21.192 ( -4.346 -4.346 -38.521) 39.008 22.295 ( 7.447 7.447 14.864) 18.217 22.370 ( -7.433 -7.433 20.027) 22.618 25.176 ( -12.938 -12.938 -0.195) 18.298 26.493 ( -4.974 -4.974 -13.234) 14.987 26.975 ( 0.594 0.594 -17.156) 17.176 30.207 ( 7.222 7.222 9.591) 14.011 ======================= Grid point 240 (48/84) ======================= q-point: ( 0.46 0.38 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.279 ( -0.200 -0.942 6.835) 6.902 3.002 ( 1.170 5.984 0.842) 6.155 3.175 ( -0.028 -4.969 -0.116) 4.971 9.978 ( 1.736 5.312 7.935) 9.706 10.654 ( 1.387 10.869 3.166) 11.406 10.765 ( 4.151 2.813 2.304) 5.518 19.488 ( 1.286 3.252 -21.572) 21.853 20.453 ( 12.504 -6.085 -21.790) 25.849 20.935 ( -16.671 -20.521 -31.102) 40.821 22.336 ( 1.370 7.507 17.147) 18.768 22.389 ( 2.076 8.195 8.566) 12.035 25.007 ( -4.325 -15.320 0.137) 15.919 26.426 ( -1.812 -6.621 -12.303) 14.089 26.995 ( 0.829 4.522 -22.158) 22.630 30.300 ( 2.367 7.456 9.599) 12.383 ======================= Grid point 253 (49/84) ======================= q-point: ( 0.46 0.46 0.08) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.266 ( -0.350 -0.350 7.296) 7.313 3.080 ( 1.795 1.795 0.694) 2.632 3.111 ( -1.392 -1.392 0.449) 2.020 10.049 ( 1.879 1.879 8.947) 9.333 10.797 ( 2.753 2.753 2.032) 4.392 10.804 ( 2.094 2.094 2.348) 3.779 19.535 ( 1.437 1.437 -16.611) 16.735 20.402 ( 1.695 1.695 -23.137) 23.261 20.565 ( -14.423 -14.423 -26.387) 33.351 22.479 ( 4.854 4.854 7.101) 9.877 22.501 ( 2.868 2.868 6.210) 7.417 24.805 ( -5.303 -5.303 0.460) 7.513 26.334 ( -2.508 -2.508 -11.959) 12.474 27.066 ( 2.128 2.128 -25.754) 25.929 30.396 ( 2.445 2.445 9.501) 10.111 ======================= Grid point 366 (50/84) ======================= q-point: ( 0.15 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.889 ( 12.696 12.696 12.696) 21.990 1.889 ( 12.696 12.696 12.696) 21.990 2.684 ( 11.522 11.522 11.522) 19.957 9.332 ( 8.679 8.679 8.679) 15.032 9.332 ( 8.679 8.679 8.679) 15.032 9.445 ( 8.900 8.900 8.900) 15.416 19.043 ( 10.595 10.595 10.595) 18.351 19.043 ( 10.595 10.595 10.595) 18.351 19.200 ( -11.179 -11.179 -11.179) 19.363 21.541 ( 21.933 21.933 21.933) 37.989 21.541 ( 21.933 21.933 21.933) 37.989 26.819 ( -8.110 -8.110 -8.110) 14.047 26.819 ( -8.110 -8.110 -8.110) 14.047 27.674 ( 5.803 5.803 5.803) 10.052 29.322 ( 5.124 5.124 5.124) 8.874 ======================= Grid point 367 (51/84) ======================= q-point: ( 0.23 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.131 ( 11.197 10.198 10.198) 18.259 2.136 ( 12.165 9.542 9.542) 18.168 2.925 ( 13.107 3.847 3.847) 14.191 9.422 ( 4.987 10.018 10.018) 15.020 9.540 ( 9.571 10.793 10.793) 18.016 9.738 ( 18.725 4.853 4.853) 19.943 18.982 ( -6.115 -9.143 -9.143) 14.303 19.318 ( 19.998 6.027 6.027) 21.739 19.451 ( 22.680 -4.617 -4.617) 23.602 21.737 ( 10.352 25.772 25.772) 37.888 22.156 ( 27.958 13.053 13.053) 33.502 26.550 ( -13.210 -8.304 -8.304) 17.675 26.770 ( -2.699 -11.303 -11.303) 16.212 27.697 ( -4.035 -0.574 -0.574) 4.116 29.429 ( 6.112 10.180 10.180) 15.640 ======================= Grid point 368 (52/84) ======================= q-point: ( 0.31 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.303 ( 6.486 7.292 7.292) 12.183 2.347 ( 9.320 7.322 7.322) 13.931 3.180 ( 12.690 -1.685 -1.685) 12.911 9.518 ( 4.698 8.978 8.978) 13.537 9.726 ( 9.343 12.214 12.214) 19.638 10.104 ( 17.966 4.126 4.126) 18.889 18.913 ( -1.410 -13.252 -13.252) 18.794 19.787 ( 27.028 3.793 3.793) 27.555 19.927 ( 25.245 -9.934 -9.934) 28.891 21.920 ( 7.793 22.470 22.470) 32.718 22.609 ( 18.371 8.448 8.448) 21.914 26.326 ( -9.361 -11.958 -11.958) 19.329 26.718 ( -2.570 -11.814 -11.814) 16.904 27.521 ( -13.312 -8.036 -8.036) 17.503 29.559 ( 6.973 14.183 14.183) 21.235 ======================= Grid point 369 (53/84) ======================= q-point: ( 0.38 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.390 ( 2.774 4.565 4.565) 7.026 2.493 ( 5.670 6.005 6.005) 10.212 3.397 ( 9.111 -4.739 -4.739) 11.311 9.598 ( 3.324 8.085 8.085) 11.907 9.887 ( 6.785 13.028 13.028) 19.634 10.408 ( 12.725 4.109 4.109) 13.989 18.909 ( 0.539 -14.602 -14.602) 20.658 20.313 ( 25.450 6.317 6.317) 26.972 20.384 ( 20.549 -14.241 -14.241) 28.773 22.019 ( 2.174 17.899 17.899) 25.407 22.872 ( 9.043 3.463 3.463) 10.284 26.203 ( -3.083 -15.129 -15.129) 21.617 26.673 ( -1.927 -12.447 -12.447) 17.708 27.223 ( -16.010 -14.659 -14.659) 26.193 29.683 ( 5.373 16.194 16.194) 23.524 ======================= Grid point 370 (54/84) ======================= q-point: ( 0.46 0.15 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.422 ( 0.716 2.973 2.973) 4.264 2.567 ( 1.886 5.408 5.408) 7.877 3.520 ( 3.274 -6.059 -6.059) 9.173 9.643 ( 1.197 7.570 7.570) 10.772 9.979 ( 2.462 13.436 13.436) 19.161 10.579 ( 4.588 4.212 4.212) 7.519 18.921 ( 0.427 -14.921 -14.921) 21.106 20.674 ( 8.096 -16.157 -16.157) 24.241 20.696 ( 11.581 11.255 11.255) 19.684 22.019 ( -1.092 12.785 12.785) 18.114 22.978 ( 2.405 0.082 0.082) 2.408 26.192 ( 0.874 -16.470 -16.470) 23.309 26.646 ( -0.737 -12.899 -12.899) 18.256 26.966 ( -8.269 -19.178 -19.178) 28.354 29.756 ( 1.973 16.913 16.913) 24.000 ======================= Grid point 380 (55/84) ======================= q-point: ( 0.23 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.290 ( 7.267 7.267 11.460) 15.394 2.359 ( 11.039 11.039 7.200) 17.192 3.017 ( 5.818 5.818 -1.267) 8.324 9.584 ( 6.945 6.945 13.050) 16.333 9.798 ( 12.547 12.547 3.045) 18.003 9.893 ( 12.898 12.898 9.233) 20.444 18.793 ( -7.899 -7.899 -25.534) 27.871 19.458 ( 10.226 10.226 2.201) 14.628 19.607 ( 16.390 16.390 -7.666) 24.414 22.086 ( 11.777 11.777 30.210) 34.497 22.494 ( 18.843 18.843 3.827) 26.922 26.311 ( -14.697 -14.697 -2.539) 20.939 26.613 ( -5.664 -5.664 -17.779) 19.500 27.585 ( -10.934 -10.934 -4.200) 16.024 29.637 ( 10.929 10.929 12.611) 19.948 ======================= Grid point 381 (56/84) ======================= q-point: ( 0.31 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.400 ( 3.875 3.695 11.103) 12.326 2.542 ( 7.613 10.781 4.654) 13.994 3.153 ( 7.842 -0.859 -3.796) 8.755 9.699 ( 4.930 9.159 11.093) 15.207 10.001 ( 9.471 15.173 7.636) 19.448 10.209 ( 16.945 6.848 6.334) 19.343 18.702 ( -1.692 -7.834 -31.281) 32.291 19.746 ( 18.398 -1.788 -1.711) 18.563 19.988 ( 22.149 10.248 -15.467) 28.893 22.254 ( 5.851 12.340 25.099) 28.574 22.789 ( 10.504 8.615 2.646) 13.840 26.037 ( -12.507 -16.729 -4.298) 21.324 26.519 ( -4.202 -8.800 -17.719) 20.225 27.293 ( -18.198 -14.851 -9.797) 25.450 29.840 ( 9.572 14.068 14.230) 22.182 ======================= Grid point 382 (57/84) ======================= q-point: ( 0.38 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.445 ( 1.034 1.628 9.177) 9.378 2.657 ( 4.352 9.735 2.979) 11.072 3.300 ( 6.686 -4.958 -4.022) 9.245 9.776 ( 3.065 9.901 9.325) 13.942 10.165 ( 7.008 14.934 10.681) 19.653 10.502 ( 12.443 5.447 4.890) 14.436 18.705 ( 1.363 -6.177 -31.622) 32.249 20.115 ( 17.779 -10.900 -11.985) 24.053 20.441 ( 23.355 4.632 -14.024) 27.633 22.304 ( -0.964 11.451 18.522) 21.797 22.900 ( 1.526 -1.432 4.222) 4.712 25.854 ( -6.019 -19.757 -6.506) 21.654 26.451 ( -2.638 -10.554 -17.645) 20.729 26.925 ( -18.354 -14.842 -16.924) 29.045 29.997 ( 6.361 15.419 15.091) 22.493 ======================= Grid point 383 (58/84) ======================= q-point: ( 0.46 0.23 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.453 ( 0.084 0.667 7.540) 7.569 2.713 ( 1.435 8.843 2.535) 9.310 3.393 ( 2.544 -6.702 -3.888) 8.155 9.817 ( 1.055 9.998 8.436) 13.124 10.260 ( 2.565 14.911 11.977) 19.296 10.671 ( 4.544 4.958 4.391) 8.031 18.732 ( 0.909 -4.393 -30.457) 30.786 20.373 ( 7.383 -14.509 -23.181) 28.326 20.818 ( 12.415 0.430 -6.275) 13.918 22.240 ( -3.696 10.035 11.896) 15.996 22.893 ( -0.918 -8.572 7.702) 11.560 25.797 ( -0.662 -22.721 -6.698) 23.697 26.426 ( 0.044 -10.245 -18.576) 21.214 26.638 ( -9.372 -13.270 -22.778) 27.978 30.081 ( 2.204 15.834 15.458) 22.238 ======================= Grid point 394 (59/84) ======================= q-point: ( 0.31 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.451 ( 1.550 1.550 11.387) 11.596 2.730 ( 7.939 7.939 4.063) 11.940 3.159 ( 1.727 1.727 -4.431) 5.059 9.859 ( 7.009 7.009 11.843) 15.443 10.271 ( 11.331 11.331 3.270) 16.354 10.384 ( 11.574 11.574 9.830) 19.093 18.620 ( -0.591 -0.591 -40.585) 40.593 19.802 ( 6.884 6.884 -4.686) 10.805 20.221 ( 13.254 13.254 -26.008) 32.058 22.426 ( 5.669 5.669 21.289) 22.749 22.883 ( 0.171 0.171 2.463) 2.475 25.717 ( -15.320 -15.320 -2.079) 21.765 26.368 ( -6.710 -6.710 -19.316) 21.521 26.977 ( -16.755 -16.755 -15.305) 28.208 30.088 ( 11.130 11.130 14.356) 21.303 ======================= Grid point 395 (60/84) ======================= q-point: ( 0.38 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.464 ( -0.043 0.433 10.623) 10.632 2.841 ( 3.723 8.679 2.634) 9.805 3.215 ( 3.466 -3.594 -2.953) 5.801 9.965 ( 3.982 9.096 11.186) 14.957 10.443 ( 6.859 13.060 6.972) 16.316 10.616 ( 10.677 6.237 6.720) 14.073 18.663 ( 4.155 1.497 -41.651) 41.884 19.972 ( 9.386 -3.725 -13.281) 16.684 20.483 ( 12.964 -0.733 -30.307) 32.972 22.475 ( -0.829 6.393 15.019) 16.344 22.807 ( -6.496 -7.558 7.268) 12.335 25.473 ( -9.236 -18.468 -2.167) 20.762 26.268 ( -3.313 -7.974 -20.137) 21.910 26.672 ( -14.039 -11.073 -23.260) 29.338 30.264 ( 6.917 11.809 14.549) 19.975 ======================= Grid point 396 (61/84) ======================= q-point: ( 0.46 0.31 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.459 ( -0.235 0.089 9.778) 9.781 2.885 ( 1.012 8.379 1.857) 8.642 3.269 ( 1.652 -5.814 -1.715) 6.282 10.016 ( 1.263 10.064 10.471) 14.578 10.536 ( 2.525 12.986 9.624) 16.359 10.764 ( 4.040 4.428 4.726) 7.633 18.739 ( 2.628 4.615 -39.098) 39.457 20.118 ( 4.502 -11.407 -28.505) 31.031 20.681 ( 6.081 -14.742 -23.055) 28.033 22.411 ( -3.705 7.575 8.852) 12.226 22.707 ( -2.676 -8.445 13.006) 15.736 25.360 ( -2.653 -21.116 -2.102) 21.386 26.245 ( 0.591 -8.254 -20.314) 21.935 26.457 ( -7.061 -5.969 -30.211) 31.594 30.354 ( 2.323 12.020 14.694) 19.126 ======================= Grid point 408 (62/84) ======================= q-point: ( 0.38 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.465 ( -0.251 -0.251 11.104) 11.109 2.976 ( 4.775 4.775 2.280) 7.127 3.169 ( -0.878 -0.878 -1.466) 1.921 10.114 ( 5.911 5.911 12.766) 15.260 10.650 ( 7.838 7.838 4.040) 11.798 10.741 ( 6.343 6.343 8.192) 12.148 18.748 ( 6.761 6.761 -42.762) 43.818 19.949 ( 0.267 0.267 -15.833) 15.838 20.406 ( -6.592 -6.592 -41.375) 42.412 22.556 ( 1.747 1.747 10.366) 10.656 22.655 ( -6.387 -6.387 9.568) 13.158 25.168 ( -12.263 -12.263 -0.575) 17.352 26.151 ( -4.016 -4.016 -21.109) 21.860 26.490 ( -7.679 -7.679 -31.929) 33.725 30.449 ( 7.175 7.175 14.384) 17.603 ======================= Grid point 409 (63/84) ======================= q-point: ( 0.46 0.38 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.459 ( -0.208 -0.062 11.153) 11.155 3.025 ( 0.811 5.788 1.593) 6.058 3.174 ( 0.661 -3.852 0.015) 3.908 10.190 ( 1.945 7.471 13.342) 15.414 10.748 ( 2.495 8.557 6.529) 11.049 10.837 ( 2.812 2.978 5.208) 6.625 18.882 ( 5.151 9.304 -38.557) 39.997 19.938 ( -1.489 -6.961 -31.020) 31.826 20.290 ( -3.530 -23.079 -33.112) 40.515 22.534 ( -1.841 4.935 5.652) 7.726 22.609 ( -0.077 -1.894 9.628) 9.813 25.009 ( -4.095 -14.416 -0.057) 14.986 26.106 ( -0.857 -5.874 -19.682) 20.558 26.381 ( -3.197 -2.436 -39.903) 40.105 30.542 ( 2.381 7.264 14.372) 16.278 ======================= Grid point 422 (64/84) ======================= q-point: ( 0.46 0.46 0.15) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.458 ( -0.041 -0.041 11.853) 11.853 3.101 ( 1.607 1.607 1.418) 2.679 3.125 ( -0.895 -0.895 0.997) 1.611 10.292 ( 2.750 2.750 15.524) 16.003 10.860 ( 2.807 2.807 4.511) 6.009 10.877 ( 1.185 1.185 5.251) 5.512 19.060 ( 7.062 7.062 -31.603) 33.144 19.842 ( -4.206 -4.206 -32.765) 33.301 19.910 ( -11.386 -11.386 -39.184) 42.363 22.597 ( 1.233 1.233 4.074) 4.431 22.604 ( 0.789 0.789 3.434) 3.611 24.817 ( -5.055 -5.055 0.743) 7.188 26.026 ( -2.202 -2.202 -18.786) 19.042 26.350 ( -0.766 -0.766 -46.670) 46.683 30.636 ( 2.404 2.404 14.273) 14.672 ======================= Grid point 549 (65/84) ======================= q-point: ( 0.23 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.506 ( 8.528 8.528 8.528) 14.770 2.506 ( 8.528 8.528 8.528) 14.770 3.007 ( 0.927 0.927 0.927) 1.606 9.864 ( 9.381 9.381 9.381) 16.249 9.864 ( 9.381 9.381 9.381) 16.249 10.113 ( 13.066 13.066 13.066) 22.630 18.369 ( -16.335 -16.335 -16.335) 28.293 19.566 ( 5.883 5.883 5.883) 10.190 19.566 ( 5.883 5.883 5.883) 10.190 22.572 ( 11.660 11.660 11.660) 20.196 22.572 ( 11.660 11.660 11.660) 20.196 26.254 ( -10.944 -10.944 -10.944) 18.955 26.254 ( -10.944 -10.944 -10.944) 18.955 27.393 ( -15.548 -15.548 -15.548) 26.930 29.890 ( 12.905 12.905 12.905) 22.353 ======================= Grid point 550 (66/84) ======================= q-point: ( 0.31 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.618 ( 1.859 8.336 8.336) 11.934 2.641 ( 6.171 7.397 7.397) 12.146 3.081 ( 6.530 -3.097 -3.097) 7.863 9.930 ( 3.148 11.951 11.951) 17.192 10.153 ( 13.958 8.082 8.082) 18.041 10.380 ( 13.759 11.010 11.010) 20.779 18.174 ( -4.138 -21.338 -21.338) 30.459 19.716 ( 7.617 -3.286 -3.286) 8.923 19.786 ( 17.757 -3.249 -3.249) 18.342 22.620 ( 0.461 11.988 11.988) 16.960 22.845 ( 8.927 3.203 3.203) 10.011 25.915 ( -15.541 -9.142 -9.142) 20.215 26.187 ( -3.448 -14.927 -14.927) 21.390 27.004 ( -23.430 -19.514 -19.514) 36.201 30.122 ( 10.644 14.188 14.188) 22.713 ======================= Grid point 551 (67/84) ======================= q-point: ( 0.38 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.612 ( -1.429 6.587 6.587) 9.425 2.740 ( 3.928 6.340 6.340) 9.789 3.220 ( 6.655 -3.995 -3.995) 8.730 9.982 ( 2.101 11.541 11.541) 16.457 10.386 ( 9.546 11.229 11.229) 18.529 10.623 ( 10.578 7.649 7.649) 15.130 18.174 ( 2.911 -21.505 -21.505) 30.552 19.848 ( 5.590 -13.822 -13.822) 20.330 20.212 ( 24.250 -10.755 -10.755) 28.625 22.568 ( -6.019 8.858 8.858) 13.898 22.942 ( 1.825 -0.291 -0.291) 1.871 25.679 ( -8.101 -11.580 -11.580) 18.270 26.126 ( -2.709 -15.094 -15.094) 21.518 26.537 ( -23.076 -21.945 -21.945) 38.674 30.293 ( 6.825 14.770 14.770) 21.975 ======================= Grid point 552 (68/84) ======================= q-point: ( 0.46 0.23 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.587 ( -0.793 5.469 5.469) 7.775 2.792 ( 1.337 5.837 5.837) 8.362 3.312 ( 2.545 -4.220 -4.220) 6.488 10.010 ( 0.729 11.249 11.249) 15.925 10.513 ( 3.346 13.305 13.305) 19.112 10.769 ( 3.991 5.549 5.549) 8.804 18.235 ( 2.248 -19.159 -19.159) 27.187 19.923 ( 2.047 -22.236 -22.236) 31.514 20.618 ( 13.623 -15.064 -15.064) 25.287 22.425 ( -6.190 7.236 7.236) 11.959 22.950 ( -0.247 -1.934 -1.934) 2.746 25.604 ( -0.544 -13.395 -13.395) 18.951 26.088 ( -1.040 -15.183 -15.183) 21.498 26.181 ( -11.235 -22.725 -22.725) 34.045 30.383 ( 2.329 14.964 14.964) 21.290 ======================= Grid point 563 (69/84) ======================= q-point: ( 0.31 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.688 ( 1.387 1.387 12.392) 12.546 2.815 ( 6.751 6.751 4.552) 10.576 3.065 ( 2.333 2.333 -5.062) 6.042 10.119 ( 6.801 6.801 14.703) 17.570 10.343 ( 11.862 11.862 4.005) 17.246 10.605 ( 11.263 11.263 12.379) 20.173 17.904 ( -6.376 -6.376 -30.899) 32.188 19.706 ( 1.183 1.183 -6.193) 6.415 19.773 ( 2.550 2.550 -18.556) 18.904 22.732 ( 0.653 0.653 10.292) 10.333 22.897 ( 1.155 1.155 -1.011) 1.920 25.668 ( -14.745 -14.745 -2.865) 21.049 25.981 ( -7.109 -7.109 -20.011) 22.395 26.568 ( -24.178 -24.178 -26.264) 43.116 30.373 ( 11.289 11.289 14.331) 21.454 ======================= Grid point 564 (70/84) ======================= q-point: ( 0.38 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.684 ( -1.304 1.614 11.392) 11.579 2.897 ( 2.544 8.425 2.981) 9.292 3.156 ( 5.474 -2.333 -2.908) 6.624 10.208 ( 2.895 11.143 13.490) 17.735 10.580 ( 11.369 8.901 6.940) 16.020 10.790 ( 7.281 8.879 10.808) 15.769 17.902 ( 4.878 -6.256 -34.207) 35.115 19.644 ( -5.865 -6.221 -20.624) 22.326 19.967 ( 14.379 -15.299 -22.347) 30.663 22.672 ( -6.444 2.548 5.605) 8.913 22.883 ( -1.938 -5.558 0.617) 5.919 25.415 ( -10.005 -14.698 -3.939) 18.211 25.871 ( -4.707 -11.590 -19.770) 23.395 26.121 ( -20.183 -19.445 -32.480) 42.900 30.552 ( 7.053 11.548 14.438) 19.788 ======================= Grid point 565 (71/84) ======================= q-point: ( 0.46 0.31 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.659 ( -0.803 2.345 10.219) 10.516 2.930 ( 0.845 7.593 2.772) 8.127 3.235 ( 2.250 -3.576 -1.685) 4.548 10.243 ( 0.848 12.418 12.736) 17.808 10.748 ( 5.036 10.251 11.379) 16.122 10.879 ( 1.946 5.639 7.221) 9.366 18.017 ( 4.642 -3.133 -32.733) 33.209 19.532 ( -4.100 -17.556 -30.634) 35.545 20.221 ( 8.698 -25.514 -25.623) 37.191 22.532 ( -5.537 4.128 3.816) 7.890 22.849 ( -1.074 -7.179 2.018) 7.535 25.298 ( -2.368 -16.746 -4.549) 17.514 25.764 ( -6.106 -17.828 -33.796) 38.694 25.871 ( -4.277 -9.615 -23.003) 25.296 30.644 ( 2.376 11.612 14.492) 18.722 ======================= Grid point 577 (72/84) ======================= q-point: ( 0.38 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.696 ( -0.067 -0.067 12.153) 12.153 3.027 ( 4.120 4.120 2.806) 6.467 3.140 ( 1.052 1.052 -1.476) 2.096 10.389 ( 6.903 6.903 15.102) 17.983 10.747 ( 8.477 8.477 5.891) 13.357 10.934 ( 5.210 5.210 11.330) 13.514 17.917 ( 6.930 6.930 -40.062) 41.244 19.598 ( -7.272 -7.272 -21.091) 23.465 19.624 ( -11.148 -11.148 -38.068) 41.204 22.680 ( -2.227 -2.227 2.615) 4.093 22.757 ( -5.643 -5.643 1.405) 8.103 25.152 ( -11.288 -11.288 -1.073) 16.000 25.646 ( -11.568 -11.568 -35.383) 38.982 25.834 ( -9.544 -9.544 -28.969) 31.959 30.735 ( 7.120 7.120 14.322) 17.507 ======================= Grid point 578 (73/84) ======================= q-point: ( 0.46 0.38 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.691 ( -0.233 1.061 12.169) 12.217 3.062 ( 0.443 5.583 2.072) 5.972 3.177 ( 1.564 -2.308 0.355) 2.811 10.473 ( 1.995 10.582 15.226) 18.649 10.897 ( 6.246 5.341 8.344) 11.712 10.977 ( -0.804 3.926 9.458) 10.272 18.097 ( 8.214 10.456 -39.592) 41.765 19.256 ( -12.943 -10.678 -37.703) 41.268 19.666 ( 2.429 -29.966 -32.233) 44.078 22.593 ( -3.720 2.265 0.616) 4.398 22.712 ( -0.562 -6.297 1.352) 6.465 25.003 ( -3.842 -12.645 -0.568) 13.228 25.471 ( -5.141 -12.016 -49.830) 51.515 25.723 ( -2.618 -5.848 -21.173) 22.121 30.827 ( 2.382 7.122 14.301) 16.153 ======================= Grid point 591 (74/84) ======================= q-point: ( 0.46 0.46 0.23) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.704 ( 0.316 0.316 12.906) 12.914 3.134 ( 1.391 1.391 1.980) 2.791 3.150 ( -0.179 -0.179 1.491) 1.513 10.626 ( 4.394 4.394 18.061) 19.100 10.975 ( 3.049 3.049 7.381) 8.549 11.013 ( -0.362 -0.362 8.749) 8.764 18.362 ( 13.558 13.558 -38.279) 42.813 19.097 ( -10.614 -10.614 -43.227) 45.759 19.163 ( -13.252 -13.252 -36.772) 41.273 22.621 ( -0.622 -0.622 -1.525) 1.761 22.625 ( -1.085 -1.085 -1.103) 1.890 24.833 ( -4.558 -4.558 0.801) 6.496 25.309 ( -4.405 -4.405 -58.432) 58.763 25.644 ( -2.086 -2.086 -19.882) 20.100 30.919 ( 2.376 2.376 14.251) 14.642 ======================= Grid point 732 (75/84) ======================= q-point: ( 0.31 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.899 ( 5.149 5.149 5.149) 8.918 2.899 ( 5.149 5.149 5.149) 8.918 2.986 ( -0.208 -0.208 -0.208) 0.360 10.419 ( 9.656 9.656 9.656) 16.724 10.419 ( 9.656 9.656 9.656) 16.724 10.838 ( 10.873 10.873 10.873) 18.832 17.473 ( -12.747 -12.747 -12.747) 22.078 19.536 ( -8.559 -8.559 -8.559) 14.825 19.536 ( -8.559 -8.559 -8.559) 14.825 22.843 ( -1.344 -1.344 -1.344) 2.329 22.843 ( -1.344 -1.344 -1.344) 2.329 25.609 ( -10.470 -10.470 -10.470) 18.134 25.609 ( -10.470 -10.470 -10.470) 18.134 25.998 ( -30.754 -30.754 -30.754) 53.268 30.626 ( 11.378 11.378 11.378) 19.708 ======================= Grid point 733 (76/84) ======================= q-point: ( 0.38 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.873 ( -3.069 6.364 6.364) 9.508 2.965 ( 2.653 5.051 5.051) 7.620 3.112 ( 7.035 -1.369 -1.369) 7.297 10.476 ( 2.223 13.487 13.487) 19.203 10.713 ( 13.964 6.895 6.895) 17.032 11.003 ( 5.622 10.195 10.195) 15.475 17.400 ( 3.988 -16.436 -16.436) 23.583 19.258 ( -16.358 -17.052 -17.052) 29.140 19.544 ( 5.100 -22.828 -22.828) 32.684 22.718 ( -8.037 -0.200 -0.200) 8.042 22.840 ( -1.479 -4.698 -4.698) 6.807 25.296 ( -14.500 -10.494 -10.494) 20.748 25.455 ( -22.570 -30.980 -30.980) 49.284 25.554 ( -2.444 -13.679 -13.679) 19.499 30.807 ( 7.092 11.402 11.402) 17.616 ======================= Grid point 734 (77/84) ======================= q-point: ( 0.46 0.31 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.830 ( -1.162 6.653 6.653) 9.480 3.000 ( 0.925 4.716 4.716) 6.734 3.206 ( 2.527 -1.253 -1.253) 3.087 10.504 ( 0.710 13.902 13.902) 19.673 10.940 ( 8.867 7.886 7.886) 14.248 11.044 ( -1.650 9.052 9.052) 12.907 17.537 ( 6.610 -15.542 -15.542) 22.952 18.970 ( -10.184 -26.294 -26.294) 38.555 19.661 ( 4.590 -31.942 -31.942) 45.405 22.572 ( -5.151 0.645 0.645) 5.231 22.812 ( -0.940 -5.088 -5.088) 7.256 25.055 ( -7.856 -32.186 -32.186) 46.191 25.191 ( -5.714 -11.008 -11.008) 16.583 25.520 ( -0.937 -13.407 -13.407) 18.984 30.899 ( 2.386 11.402 11.402) 16.300 ======================= Grid point 746 (78/84) ======================= q-point: ( 0.38 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.918 ( 0.131 0.131 10.284) 10.286 3.081 ( 3.565 3.565 2.659) 5.701 3.115 ( 3.053 3.053 -0.942) 4.419 10.690 ( 7.321 7.321 15.653) 18.767 10.870 ( 9.887 9.887 6.476) 15.409 11.156 ( 4.985 4.985 10.760) 12.864 17.273 ( 2.789 2.789 -24.641) 24.955 18.964 ( -18.912 -18.912 -29.079) 39.508 19.116 ( -14.817 -14.817 -28.445) 35.330 22.679 ( -4.407 -4.407 -2.019) 6.551 22.729 ( -4.843 -4.843 -3.760) 7.813 24.829 ( -25.638 -25.638 -43.697) 56.781 25.128 ( -10.183 -10.183 -1.366) 14.465 25.399 ( -4.351 -4.351 -18.436) 19.436 30.987 ( 7.081 7.081 11.347) 15.134 ======================= Grid point 747 (79/84) ======================= q-point: ( 0.46 0.38 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 1099 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.913 ( -0.280 2.426 10.277) 10.563 3.104 ( 0.397 5.158 2.195) 5.620 3.187 ( 2.249 -0.700 0.637) 2.440 10.769 ( 1.689 12.818 15.101) 19.879 11.063 ( 8.316 4.977 8.528) 12.910 11.185 ( -1.933 4.647 10.908) 12.013 17.430 ( 9.176 4.169 -27.329) 29.128 18.539 ( -16.971 -17.619 -34.716) 42.470 19.006 ( -1.264 -34.168 -35.669) 49.410 22.571 ( -4.098 -0.489 -2.474) 4.812 22.683 ( -0.953 -7.336 -3.718) 8.279 24.480 ( -9.293 -24.172 -49.673) 56.018 24.987 ( -3.840 -9.913 -1.226) 10.701 25.347 ( -1.263 -5.320 -17.131) 17.982 31.079 ( 2.378 7.065 11.336) 13.567 ======================= Grid point 760 (80/84) ======================= q-point: ( 0.46 0.46 0.31) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.941 ( 0.625 0.625 11.090) 11.125 3.174 ( 1.213 1.213 2.050) 2.673 3.181 ( 0.586 0.586 1.611) 1.811 10.967 ( 6.196 6.196 16.563) 18.738 11.143 ( 3.714 3.714 9.560) 10.908 11.211 ( -1.750 -1.750 11.200) 11.470 17.648 ( 14.853 14.853 -33.449) 39.498 18.276 ( -12.056 -12.056 -40.143) 43.614 18.429 ( -18.919 -18.919 -38.082) 46.541 22.558 ( -1.547 -1.547 -4.418) 4.930 22.568 ( -2.488 -2.488 -4.216) 5.492 24.158 ( -8.473 -8.473 -57.619) 58.852 24.847 ( -3.880 -3.880 0.685) 5.529 25.284 ( -1.527 -1.527 -16.750) 16.888 31.171 ( 2.370 2.370 11.334) 11.819 ======================= Grid point 915 (81/84) ======================= q-point: ( 0.38 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.062 ( 2.651 2.651 2.651) 4.591 3.126 ( 2.749 2.749 2.749) 4.762 3.126 ( 2.749 2.749 2.749) 4.762 10.984 ( 9.556 9.556 9.556) 16.552 10.984 ( 9.556 9.556 9.556) 16.552 11.322 ( 5.644 5.644 5.644) 9.775 16.997 ( -3.614 -3.614 -3.614) 6.260 18.515 ( -24.544 -24.544 -24.544) 42.512 18.515 ( -24.544 -24.544 -24.544) 42.512 22.627 ( -4.559 -4.559 -4.559) 7.897 22.627 ( -4.559 -4.559 -4.559) 7.897 24.064 ( -32.644 -32.644 -32.644) 56.540 25.103 ( -6.507 -6.507 -6.507) 11.271 25.103 ( -6.507 -6.507 -6.507) 11.271 31.167 ( 7.068 7.068 7.068) 12.242 ======================= Grid point 916 (82/84) ======================= q-point: ( 0.46 0.38 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.071 ( -0.104 5.312 5.312) 7.513 3.151 ( 0.667 3.000 3.000) 4.295 3.200 ( 2.305 0.626 0.626) 2.469 11.051 ( 1.760 13.416 13.416) 19.055 11.217 ( 8.782 6.995 6.995) 13.228 11.362 ( -1.354 6.437 6.437) 9.204 17.084 ( 8.810 -7.571 -7.571) 13.865 17.952 ( -21.374 -24.689 -24.689) 40.938 18.304 ( -5.188 -35.396 -35.396) 50.325 22.515 ( -3.931 -2.872 -2.872) 5.653 22.586 ( -1.173 -5.891 -5.891) 8.414 23.629 ( -11.326 -35.667 -35.667) 51.696 24.946 ( -4.104 -3.732 -3.732) 6.686 25.086 ( -0.458 -8.568 -8.568) 12.126 31.259 ( 2.379 7.069 7.069) 10.276 ======================= Grid point 929 (83/84) ======================= q-point: ( 0.46 0.46 0.38) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 595 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.122 ( 0.889 0.889 7.328) 7.435 3.209 ( 1.165 1.165 1.234) 2.058 3.210 ( 1.095 1.095 1.503) 2.158 11.257 ( 6.218 6.218 13.059) 15.744 11.328 ( 4.949 4.949 8.888) 11.313 11.415 ( -0.967 -0.967 8.793) 8.899 17.125 ( 9.683 9.683 -19.007) 23.427 17.587 ( -15.024 -15.024 -29.329) 36.216 17.721 ( -18.777 -18.777 -33.461) 42.718 22.469 ( -2.223 -2.223 -4.259) 5.293 22.480 ( -3.052 -3.052 -4.375) 6.146 23.163 ( -11.725 -11.725 -42.492) 45.613 24.858 ( -3.229 -3.229 0.457) 4.590 25.010 ( -1.090 -1.090 -11.025) 11.133 31.351 ( 2.381 2.381 7.090) 7.848 ======================= Grid point 1098 (84/84) ======================= q-point: ( 0.46 0.46 0.46) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 1.15e-04 1.15e-04 1.15e-04 1.46e-03 6.24e-05 Number of triplets: 231 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 3.215 ( 1.802 1.802 1.802) 3.121 3.230 ( 0.871 0.871 0.871) 1.508 3.230 ( 0.871 0.871 0.871) 1.508 11.458 ( 5.390 5.390 5.390) 9.335 11.458 ( 5.390 5.390 5.390) 9.335 11.520 ( 1.476 1.476 1.476) 2.556 16.943 ( 0.414 0.414 0.414) 0.716 17.186 ( -15.417 -15.417 -15.417) 26.704 17.186 ( -15.417 -15.417 -15.417) 26.704 22.409 ( -2.319 -2.319 -2.319) 4.017 22.409 ( -2.319 -2.319 -2.319) 4.017 22.591 ( -14.609 -14.609 -14.609) 25.304 24.865 ( -1.836 -1.836 -1.836) 3.180 24.865 ( -1.836 -1.836 -1.836) 3.180 31.444 ( 2.393 2.393 2.393) 4.144 =================== 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/32955 10.0 1915.671 1915.671 1915.671 0.000 0.000 0.000 3/32955 20.0 448.217 448.217 448.217 0.000 0.000 0.000 3/32955 30.0 259.446 259.446 259.446 0.000 0.000 0.000 3/32955 40.0 215.887 215.887 215.887 0.000 0.000 0.000 3/32955 50.0 195.847 195.847 195.847 0.000 0.000 0.000 3/32955 60.0 182.080 182.080 182.080 0.000 0.000 0.000 3/32955 70.0 170.518 170.518 170.518 0.000 0.000 0.000 3/32955 80.0 159.074 159.074 159.074 0.000 -0.000 0.000 3/32955 90.0 145.905 145.905 145.905 0.000 -0.000 0.000 3/32955 100.0 129.946 129.946 129.946 0.000 -0.000 0.000 3/32955 110.0 111.750 111.750 111.750 0.000 -0.000 0.000 3/32955 120.0 93.162 93.162 93.162 0.000 -0.000 0.000 3/32955 130.0 76.074 76.074 76.074 0.000 -0.000 0.000 3/32955 140.0 61.564 61.564 61.564 0.000 -0.000 0.000 3/32955 150.0 49.852 49.852 49.852 0.000 -0.000 0.000 3/32955 160.0 40.660 40.660 40.660 0.000 -0.000 0.000 3/32955 170.0 33.527 33.527 33.527 0.000 -0.000 0.000 3/32955 180.0 27.999 27.999 27.999 0.000 -0.000 0.000 3/32955 190.0 23.692 23.692 23.692 0.000 -0.000 0.000 3/32955 200.0 20.309 20.309 20.309 0.000 -0.000 0.000 3/32955 210.0 17.622 17.622 17.622 0.000 -0.000 0.000 3/32955 220.0 15.465 15.465 15.465 0.000 -0.000 0.000 3/32955 230.0 13.714 13.714 13.714 0.000 -0.000 0.000 3/32955 240.0 12.276 12.276 12.276 0.000 -0.000 0.000 3/32955 250.0 11.084 11.084 11.084 0.000 -0.000 0.000 3/32955 260.0 10.085 10.085 10.085 0.000 -0.000 0.000 3/32955 270.0 9.240 9.240 9.240 0.000 -0.000 0.000 3/32955 280.0 8.518 8.518 8.518 0.000 -0.000 0.000 3/32955 290.0 7.898 7.898 7.898 0.000 -0.000 0.000 3/32955 300.0 7.360 7.360 7.360 0.000 -0.000 0.000 3/32955 310.0 6.890 6.890 6.890 0.000 -0.000 0.000 3/32955 320.0 6.477 6.477 6.477 0.000 -0.000 0.000 3/32955 330.0 6.112 6.112 6.112 0.000 -0.000 0.000 3/32955 340.0 5.787 5.787 5.787 0.000 -0.000 0.000 3/32955 350.0 5.496 5.496 5.496 0.000 -0.000 0.000 3/32955 360.0 5.235 5.235 5.235 0.000 -0.000 0.000 3/32955 370.0 4.999 4.999 4.999 0.000 -0.000 0.000 3/32955 380.0 4.784 4.784 4.784 0.000 -0.000 0.000 3/32955 390.0 4.589 4.589 4.589 0.000 -0.000 0.000 3/32955 400.0 4.410 4.410 4.410 0.000 -0.000 0.000 3/32955 410.0 4.246 4.246 4.246 0.000 -0.000 0.000 3/32955 420.0 4.095 4.095 4.095 0.000 -0.000 0.000 3/32955 430.0 3.955 3.955 3.955 0.000 -0.000 0.000 3/32955 440.0 3.826 3.826 3.826 0.000 -0.000 0.000 3/32955 450.0 3.705 3.705 3.705 0.000 -0.000 0.000 3/32955 460.0 3.593 3.593 3.593 0.000 -0.000 0.000 3/32955 470.0 3.488 3.488 3.488 0.000 -0.000 0.000 3/32955 480.0 3.390 3.390 3.390 0.000 -0.000 0.000 3/32955 490.0 3.297 3.297 3.297 0.000 -0.000 0.000 3/32955 500.0 3.211 3.211 3.211 0.000 -0.000 0.000 3/32955 510.0 3.129 3.129 3.129 0.000 -0.000 0.000 3/32955 520.0 3.052 3.052 3.052 0.000 -0.000 0.000 3/32955 530.0 2.978 2.978 2.978 0.000 -0.000 0.000 3/32955 540.0 2.909 2.909 2.909 0.000 -0.000 0.000 3/32955 550.0 2.843 2.843 2.843 0.000 -0.000 0.000 3/32955 560.0 2.781 2.781 2.781 0.000 -0.000 0.000 3/32955 570.0 2.721 2.721 2.721 0.000 -0.000 0.000 3/32955 580.0 2.665 2.665 2.665 0.000 -0.000 0.000 3/32955 590.0 2.611 2.611 2.611 0.000 -0.000 0.000 3/32955 600.0 2.559 2.559 2.559 0.000 -0.000 0.000 3/32955 610.0 2.509 2.509 2.509 0.000 -0.000 0.000 3/32955 620.0 2.462 2.462 2.462 0.000 -0.000 0.000 3/32955 630.0 2.416 2.416 2.416 0.000 -0.000 0.000 3/32955 640.0 2.373 2.373 2.373 0.000 -0.000 0.000 3/32955 650.0 2.331 2.331 2.331 0.000 -0.000 0.000 3/32955 660.0 2.291 2.291 2.291 0.000 -0.000 0.000 3/32955 670.0 2.252 2.252 2.252 0.000 -0.000 0.000 3/32955 680.0 2.214 2.214 2.214 0.000 -0.000 0.000 3/32955 690.0 2.178 2.178 2.178 0.000 -0.000 0.000 3/32955 700.0 2.144 2.144 2.144 0.000 -0.000 0.000 3/32955 710.0 2.110 2.110 2.110 0.000 -0.000 0.000 3/32955 720.0 2.078 2.078 2.078 0.000 -0.000 0.000 3/32955 730.0 2.046 2.046 2.046 0.000 -0.000 0.000 3/32955 740.0 2.016 2.016 2.016 0.000 -0.000 0.000 3/32955 750.0 1.987 1.987 1.987 0.000 -0.000 0.000 3/32955 760.0 1.958 1.958 1.958 0.000 -0.000 0.000 3/32955 770.0 1.931 1.931 1.931 0.000 -0.000 0.000 3/32955 780.0 1.904 1.904 1.904 0.000 -0.000 0.000 3/32955 790.0 1.878 1.878 1.878 0.000 -0.000 0.000 3/32955 800.0 1.853 1.853 1.853 0.000 -0.000 0.000 3/32955 810.0 1.828 1.828 1.828 0.000 -0.000 0.000 3/32955 820.0 1.804 1.804 1.804 0.000 -0.000 0.000 3/32955 830.0 1.781 1.781 1.781 0.000 -0.000 0.000 3/32955 840.0 1.759 1.759 1.759 0.000 -0.000 0.000 3/32955 850.0 1.737 1.737 1.737 0.000 -0.000 0.000 3/32955 860.0 1.716 1.716 1.716 0.000 -0.000 0.000 3/32955 870.0 1.695 1.695 1.695 0.000 -0.000 0.000 3/32955 880.0 1.675 1.675 1.675 0.000 -0.000 0.000 3/32955 890.0 1.655 1.655 1.655 0.000 -0.000 0.000 3/32955 900.0 1.636 1.636 1.636 0.000 -0.000 0.000 3/32955 910.0 1.617 1.617 1.617 0.000 -0.000 0.000 3/32955 920.0 1.599 1.599 1.599 0.000 -0.000 0.000 3/32955 930.0 1.581 1.581 1.581 0.000 -0.000 0.000 3/32955 940.0 1.563 1.563 1.563 0.000 -0.000 0.000 3/32955 950.0 1.546 1.546 1.546 0.000 -0.000 0.000 3/32955 960.0 1.530 1.530 1.530 0.000 -0.000 0.000 3/32955 970.0 1.514 1.514 1.514 0.000 -0.000 0.000 3/32955 980.0 1.498 1.498 1.498 0.000 -0.000 0.000 3/32955 990.0 1.482 1.482 1.482 0.000 -0.000 0.000 3/32955 1000.0 1.467 1.467 1.467 0.000 -0.000 0.000 3/32955 Thermal conductivity related properties were written into "kappa-m131313.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 05:39:22]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|