----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 05:19:48]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [3 3 3] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fm-3m (225) Number of symmetry operations in supercell: 5184 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 2.811709859999999 2.811709859999999 b 2.811709859999999 0.000000000000000 2.811709859999999 c 2.811709859999999 2.811709859999999 0.000000000000000 Atomic positions (fractional): *1 Ar 0.00000000000000 0.00000000000000 0.00000000000000 39.948 -------------------------------- unit cell --------------------------------- Lattice vectors: a 5.623419719999999 0.000000000000000 0.000000000000000 b 0.000000000000000 5.623419719999999 0.000000000000000 c 0.000000000000000 0.000000000000000 5.623419719999999 Atomic positions (fractional): *1 Ar 0.00000000000000 0.00000000000000 0.00000000000000 39.948 > 1 2 Ar 0.00000000000000 0.50000000000000 0.50000000000000 39.948 > 1 3 Ar 0.50000000000000 0.00000000000000 0.50000000000000 39.948 > 1 4 Ar 0.50000000000000 0.50000000000000 0.00000000000000 39.948 > 1 -------------------------------- super cell -------------------------------- Lattice vectors: a 16.870259159999996 0.000000000000000 0.000000000000000 b 0.000000000000000 16.870259159999996 0.000000000000000 c 0.000000000000000 0.000000000000000 16.870259159999996 Atomic positions (fractional): *1 Ar 0.00000000000000 0.00000000000000 0.00000000000000 39.948 > 1 2 Ar 0.33333333333333 0.00000000000000 0.00000000000000 39.948 > 1 3 Ar 0.66666666666667 0.00000000000000 0.00000000000000 39.948 > 1 4 Ar 0.00000000000000 0.33333333333333 0.00000000000000 39.948 > 1 5 Ar 0.33333333333333 0.33333333333333 0.00000000000000 39.948 > 1 6 Ar 0.66666666666667 0.33333333333333 0.00000000000000 39.948 > 1 7 Ar 0.00000000000000 0.66666666666667 0.00000000000000 39.948 > 1 8 Ar 0.33333333333333 0.66666666666667 0.00000000000000 39.948 > 1 9 Ar 0.66666666666667 0.66666666666667 0.00000000000000 39.948 > 1 10 Ar 0.00000000000000 0.00000000000000 0.33333333333333 39.948 > 1 11 Ar 0.33333333333333 0.00000000000000 0.33333333333333 39.948 > 1 12 Ar 0.66666666666667 0.00000000000000 0.33333333333333 39.948 > 1 13 Ar 0.00000000000000 0.33333333333333 0.33333333333333 39.948 > 1 14 Ar 0.33333333333333 0.33333333333333 0.33333333333333 39.948 > 1 15 Ar 0.66666666666667 0.33333333333333 0.33333333333333 39.948 > 1 16 Ar 0.00000000000000 0.66666666666667 0.33333333333333 39.948 > 1 17 Ar 0.33333333333333 0.66666666666667 0.33333333333333 39.948 > 1 18 Ar 0.66666666666667 0.66666666666667 0.33333333333333 39.948 > 1 19 Ar 0.00000000000000 0.00000000000000 0.66666666666667 39.948 > 1 20 Ar 0.33333333333333 0.00000000000000 0.66666666666667 39.948 > 1 21 Ar 0.66666666666667 0.00000000000000 0.66666666666667 39.948 > 1 22 Ar 0.00000000000000 0.33333333333333 0.66666666666667 39.948 > 1 23 Ar 0.33333333333333 0.33333333333333 0.66666666666667 39.948 > 1 24 Ar 0.66666666666667 0.33333333333333 0.66666666666667 39.948 > 1 25 Ar 0.00000000000000 0.66666666666667 0.66666666666667 39.948 > 1 26 Ar 0.33333333333333 0.66666666666667 0.66666666666667 39.948 > 1 27 Ar 0.66666666666667 0.66666666666667 0.66666666666667 39.948 > 1 28 Ar 0.00000000000000 0.16666666666667 0.16666666666667 39.948 > 1 29 Ar 0.33333333333333 0.16666666666667 0.16666666666667 39.948 > 1 30 Ar 0.66666666666667 0.16666666666667 0.16666666666667 39.948 > 1 31 Ar 0.00000000000000 0.50000000000000 0.16666666666667 39.948 > 1 32 Ar 0.33333333333333 0.50000000000000 0.16666666666667 39.948 > 1 33 Ar 0.66666666666667 0.50000000000000 0.16666666666667 39.948 > 1 34 Ar 0.00000000000000 0.83333333333333 0.16666666666667 39.948 > 1 35 Ar 0.33333333333333 0.83333333333333 0.16666666666667 39.948 > 1 36 Ar 0.66666666666667 0.83333333333333 0.16666666666667 39.948 > 1 37 Ar 0.00000000000000 0.16666666666667 0.50000000000000 39.948 > 1 38 Ar 0.33333333333333 0.16666666666667 0.50000000000000 39.948 > 1 39 Ar 0.66666666666667 0.16666666666667 0.50000000000000 39.948 > 1 40 Ar 0.00000000000000 0.50000000000000 0.50000000000000 39.948 > 1 41 Ar 0.33333333333333 0.50000000000000 0.50000000000000 39.948 > 1 42 Ar 0.66666666666667 0.50000000000000 0.50000000000000 39.948 > 1 43 Ar 0.00000000000000 0.83333333333333 0.50000000000000 39.948 > 1 44 Ar 0.33333333333333 0.83333333333333 0.50000000000000 39.948 > 1 45 Ar 0.66666666666667 0.83333333333333 0.50000000000000 39.948 > 1 46 Ar 0.00000000000000 0.16666666666667 0.83333333333333 39.948 > 1 47 Ar 0.33333333333333 0.16666666666667 0.83333333333333 39.948 > 1 48 Ar 0.66666666666667 0.16666666666667 0.83333333333333 39.948 > 1 49 Ar 0.00000000000000 0.50000000000000 0.83333333333333 39.948 > 1 50 Ar 0.33333333333333 0.50000000000000 0.83333333333333 39.948 > 1 51 Ar 0.66666666666667 0.50000000000000 0.83333333333333 39.948 > 1 52 Ar 0.00000000000000 0.83333333333333 0.83333333333333 39.948 > 1 53 Ar 0.33333333333333 0.83333333333333 0.83333333333333 39.948 > 1 54 Ar 0.66666666666667 0.83333333333333 0.83333333333333 39.948 > 1 55 Ar 0.16666666666667 0.00000000000000 0.16666666666667 39.948 > 1 56 Ar 0.50000000000000 0.00000000000000 0.16666666666667 39.948 > 1 57 Ar 0.83333333333333 0.00000000000000 0.16666666666667 39.948 > 1 58 Ar 0.16666666666667 0.33333333333333 0.16666666666667 39.948 > 1 59 Ar 0.50000000000000 0.33333333333333 0.16666666666667 39.948 > 1 60 Ar 0.83333333333333 0.33333333333333 0.16666666666667 39.948 > 1 61 Ar 0.16666666666667 0.66666666666667 0.16666666666667 39.948 > 1 62 Ar 0.50000000000000 0.66666666666667 0.16666666666667 39.948 > 1 63 Ar 0.83333333333333 0.66666666666667 0.16666666666667 39.948 > 1 64 Ar 0.16666666666667 0.00000000000000 0.50000000000000 39.948 > 1 65 Ar 0.50000000000000 0.00000000000000 0.50000000000000 39.948 > 1 66 Ar 0.83333333333333 0.00000000000000 0.50000000000000 39.948 > 1 67 Ar 0.16666666666667 0.33333333333333 0.50000000000000 39.948 > 1 68 Ar 0.50000000000000 0.33333333333333 0.50000000000000 39.948 > 1 69 Ar 0.83333333333333 0.33333333333333 0.50000000000000 39.948 > 1 70 Ar 0.16666666666667 0.66666666666667 0.50000000000000 39.948 > 1 71 Ar 0.50000000000000 0.66666666666667 0.50000000000000 39.948 > 1 72 Ar 0.83333333333333 0.66666666666667 0.50000000000000 39.948 > 1 73 Ar 0.16666666666667 0.00000000000000 0.83333333333333 39.948 > 1 74 Ar 0.50000000000000 0.00000000000000 0.83333333333333 39.948 > 1 75 Ar 0.83333333333333 0.00000000000000 0.83333333333333 39.948 > 1 76 Ar 0.16666666666667 0.33333333333333 0.83333333333333 39.948 > 1 77 Ar 0.50000000000000 0.33333333333333 0.83333333333333 39.948 > 1 78 Ar 0.83333333333333 0.33333333333333 0.83333333333333 39.948 > 1 79 Ar 0.16666666666667 0.66666666666667 0.83333333333333 39.948 > 1 80 Ar 0.50000000000000 0.66666666666667 0.83333333333333 39.948 > 1 81 Ar 0.83333333333333 0.66666666666667 0.83333333333333 39.948 > 1 82 Ar 0.16666666666667 0.16666666666667 0.00000000000000 39.948 > 1 83 Ar 0.50000000000000 0.16666666666667 0.00000000000000 39.948 > 1 84 Ar 0.83333333333333 0.16666666666667 0.00000000000000 39.948 > 1 85 Ar 0.16666666666667 0.50000000000000 0.00000000000000 39.948 > 1 86 Ar 0.50000000000000 0.50000000000000 0.00000000000000 39.948 > 1 87 Ar 0.83333333333333 0.50000000000000 0.00000000000000 39.948 > 1 88 Ar 0.16666666666667 0.83333333333333 0.00000000000000 39.948 > 1 89 Ar 0.50000000000000 0.83333333333333 0.00000000000000 39.948 > 1 90 Ar 0.83333333333333 0.83333333333333 0.00000000000000 39.948 > 1 91 Ar 0.16666666666667 0.16666666666667 0.33333333333333 39.948 > 1 92 Ar 0.50000000000000 0.16666666666667 0.33333333333333 39.948 > 1 93 Ar 0.83333333333333 0.16666666666667 0.33333333333333 39.948 > 1 94 Ar 0.16666666666667 0.50000000000000 0.33333333333333 39.948 > 1 95 Ar 0.50000000000000 0.50000000000000 0.33333333333333 39.948 > 1 96 Ar 0.83333333333333 0.50000000000000 0.33333333333333 39.948 > 1 97 Ar 0.16666666666667 0.83333333333333 0.33333333333333 39.948 > 1 98 Ar 0.50000000000000 0.83333333333333 0.33333333333333 39.948 > 1 99 Ar 0.83333333333333 0.83333333333333 0.33333333333333 39.948 > 1 100 Ar 0.16666666666667 0.16666666666667 0.66666666666667 39.948 > 1 101 Ar 0.50000000000000 0.16666666666667 0.66666666666667 39.948 > 1 102 Ar 0.83333333333333 0.16666666666667 0.66666666666667 39.948 > 1 103 Ar 0.16666666666667 0.50000000000000 0.66666666666667 39.948 > 1 104 Ar 0.50000000000000 0.50000000000000 0.66666666666667 39.948 > 1 105 Ar 0.83333333333333 0.50000000000000 0.66666666666667 39.948 > 1 106 Ar 0.16666666666667 0.83333333333333 0.66666666666667 39.948 > 1 107 Ar 0.50000000000000 0.83333333333333 0.66666666666667 39.948 > 1 108 Ar 0.83333333333333 0.83333333333333 0.66666666666667 39.948 > 1 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 1.5997120 0.0000000 0.0000000 0.0000000 1.5997120 0.0000000 0.0000000 0.0000000 1.5997120 -------------------------- Born effective charges -------------------------- 1 Ar 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 ---------------------------------------------------------------------------- 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: 324/324 Permutation basis: 966/966 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 26 Number of blocks in projector: 26 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: 20 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 6 Use standard eigh solver. Tree of FC basis block matrices: - (26, 25), data: False |-- (6, 6), data: True |-- (20, 19), data: True ----- Solver_atoms: 1 -- 108 / 108 Time (Solver_compr_matrix_reshape): 0.002 Solver_block: 80 / 80 - Time: 0.007 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.010 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 324/324 Permutation basis: 966/966 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 26 Number of blocks in projector: 26 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: 20 Use standard eigh solver. --- Eigsh_solver_block: 2 / 2 --- Block_size: 6 Use standard eigh solver. Tree of FC basis block matrices: - (26, 25), data: False |-- (6, 6), data: True |-- (20, 19), 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:19:54]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 05:19:54]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [3 3 3] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fm-3m (225) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 2.811709859999999 2.811709859999999 b 2.811709859999999 0.000000000000000 2.811709859999999 c 2.811709859999999 2.811709859999999 0.000000000000000 Atomic positions (fractional): 1 Ar 0.00000000000000 0.00000000000000 0.00000000000000 39.948 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.870259159999996 0.000000000000000 0.000000000000000 b 0.000000000000000 16.870259159999996 0.000000000000000 c 0.000000000000000 0.000000000000000 16.870259159999996 Atomic positions (fractional): 1 Ar 0.00000000000000 0.00000000000000 0.00000000000000 39.948 > 1 2 Ar 0.33333333333333 0.00000000000000 0.00000000000000 39.948 > 1 3 Ar 0.66666666666667 0.00000000000000 0.00000000000000 39.948 > 1 4 Ar 0.00000000000000 0.33333333333333 0.00000000000000 39.948 > 1 5 Ar 0.33333333333333 0.33333333333333 0.00000000000000 39.948 > 1 6 Ar 0.66666666666667 0.33333333333333 0.00000000000000 39.948 > 1 7 Ar 0.00000000000000 0.66666666666667 0.00000000000000 39.948 > 1 8 Ar 0.33333333333333 0.66666666666667 0.00000000000000 39.948 > 1 9 Ar 0.66666666666667 0.66666666666667 0.00000000000000 39.948 > 1 10 Ar 0.00000000000000 0.00000000000000 0.33333333333333 39.948 > 1 11 Ar 0.33333333333333 0.00000000000000 0.33333333333333 39.948 > 1 12 Ar 0.66666666666667 0.00000000000000 0.33333333333333 39.948 > 1 13 Ar 0.00000000000000 0.33333333333333 0.33333333333333 39.948 > 1 14 Ar 0.33333333333333 0.33333333333333 0.33333333333333 39.948 > 1 15 Ar 0.66666666666667 0.33333333333333 0.33333333333333 39.948 > 1 16 Ar 0.00000000000000 0.66666666666667 0.33333333333333 39.948 > 1 17 Ar 0.33333333333333 0.66666666666667 0.33333333333333 39.948 > 1 18 Ar 0.66666666666667 0.66666666666667 0.33333333333333 39.948 > 1 19 Ar 0.00000000000000 0.00000000000000 0.66666666666667 39.948 > 1 20 Ar 0.33333333333333 0.00000000000000 0.66666666666667 39.948 > 1 21 Ar 0.66666666666667 0.00000000000000 0.66666666666667 39.948 > 1 22 Ar 0.00000000000000 0.33333333333333 0.66666666666667 39.948 > 1 23 Ar 0.33333333333333 0.33333333333333 0.66666666666667 39.948 > 1 24 Ar 0.66666666666667 0.33333333333333 0.66666666666667 39.948 > 1 25 Ar 0.00000000000000 0.66666666666667 0.66666666666667 39.948 > 1 26 Ar 0.33333333333333 0.66666666666667 0.66666666666667 39.948 > 1 27 Ar 0.66666666666667 0.66666666666667 0.66666666666667 39.948 > 1 28 Ar 0.00000000000000 0.16666666666667 0.16666666666667 39.948 > 1 29 Ar 0.33333333333333 0.16666666666667 0.16666666666667 39.948 > 1 30 Ar 0.66666666666667 0.16666666666667 0.16666666666667 39.948 > 1 31 Ar 0.00000000000000 0.50000000000000 0.16666666666667 39.948 > 1 32 Ar 0.33333333333333 0.50000000000000 0.16666666666667 39.948 > 1 33 Ar 0.66666666666667 0.50000000000000 0.16666666666667 39.948 > 1 34 Ar 0.00000000000000 0.83333333333333 0.16666666666667 39.948 > 1 35 Ar 0.33333333333333 0.83333333333333 0.16666666666667 39.948 > 1 36 Ar 0.66666666666667 0.83333333333333 0.16666666666667 39.948 > 1 37 Ar 0.00000000000000 0.16666666666667 0.50000000000000 39.948 > 1 38 Ar 0.33333333333333 0.16666666666667 0.50000000000000 39.948 > 1 39 Ar 0.66666666666667 0.16666666666667 0.50000000000000 39.948 > 1 40 Ar 0.00000000000000 0.50000000000000 0.50000000000000 39.948 > 1 41 Ar 0.33333333333333 0.50000000000000 0.50000000000000 39.948 > 1 42 Ar 0.66666666666667 0.50000000000000 0.50000000000000 39.948 > 1 43 Ar 0.00000000000000 0.83333333333333 0.50000000000000 39.948 > 1 44 Ar 0.33333333333333 0.83333333333333 0.50000000000000 39.948 > 1 45 Ar 0.66666666666667 0.83333333333333 0.50000000000000 39.948 > 1 46 Ar 0.00000000000000 0.16666666666667 0.83333333333333 39.948 > 1 47 Ar 0.33333333333333 0.16666666666667 0.83333333333333 39.948 > 1 48 Ar 0.66666666666667 0.16666666666667 0.83333333333333 39.948 > 1 49 Ar 0.00000000000000 0.50000000000000 0.83333333333333 39.948 > 1 50 Ar 0.33333333333333 0.50000000000000 0.83333333333333 39.948 > 1 51 Ar 0.66666666666667 0.50000000000000 0.83333333333333 39.948 > 1 52 Ar 0.00000000000000 0.83333333333333 0.83333333333333 39.948 > 1 53 Ar 0.33333333333333 0.83333333333333 0.83333333333333 39.948 > 1 54 Ar 0.66666666666667 0.83333333333333 0.83333333333333 39.948 > 1 55 Ar 0.16666666666667 0.00000000000000 0.16666666666667 39.948 > 1 56 Ar 0.50000000000000 0.00000000000000 0.16666666666667 39.948 > 1 57 Ar 0.83333333333333 0.00000000000000 0.16666666666667 39.948 > 1 58 Ar 0.16666666666667 0.33333333333333 0.16666666666667 39.948 > 1 59 Ar 0.50000000000000 0.33333333333333 0.16666666666667 39.948 > 1 60 Ar 0.83333333333333 0.33333333333333 0.16666666666667 39.948 > 1 61 Ar 0.16666666666667 0.66666666666667 0.16666666666667 39.948 > 1 62 Ar 0.50000000000000 0.66666666666667 0.16666666666667 39.948 > 1 63 Ar 0.83333333333333 0.66666666666667 0.16666666666667 39.948 > 1 64 Ar 0.16666666666667 0.00000000000000 0.50000000000000 39.948 > 1 65 Ar 0.50000000000000 0.00000000000000 0.50000000000000 39.948 > 1 66 Ar 0.83333333333333 0.00000000000000 0.50000000000000 39.948 > 1 67 Ar 0.16666666666667 0.33333333333333 0.50000000000000 39.948 > 1 68 Ar 0.50000000000000 0.33333333333333 0.50000000000000 39.948 > 1 69 Ar 0.83333333333333 0.33333333333333 0.50000000000000 39.948 > 1 70 Ar 0.16666666666667 0.66666666666667 0.50000000000000 39.948 > 1 71 Ar 0.50000000000000 0.66666666666667 0.50000000000000 39.948 > 1 72 Ar 0.83333333333333 0.66666666666667 0.50000000000000 39.948 > 1 73 Ar 0.16666666666667 0.00000000000000 0.83333333333333 39.948 > 1 74 Ar 0.50000000000000 0.00000000000000 0.83333333333333 39.948 > 1 75 Ar 0.83333333333333 0.00000000000000 0.83333333333333 39.948 > 1 76 Ar 0.16666666666667 0.33333333333333 0.83333333333333 39.948 > 1 77 Ar 0.50000000000000 0.33333333333333 0.83333333333333 39.948 > 1 78 Ar 0.83333333333333 0.33333333333333 0.83333333333333 39.948 > 1 79 Ar 0.16666666666667 0.66666666666667 0.83333333333333 39.948 > 1 80 Ar 0.50000000000000 0.66666666666667 0.83333333333333 39.948 > 1 81 Ar 0.83333333333333 0.66666666666667 0.83333333333333 39.948 > 1 82 Ar 0.16666666666667 0.16666666666667 0.00000000000000 39.948 > 1 83 Ar 0.50000000000000 0.16666666666667 0.00000000000000 39.948 > 1 84 Ar 0.83333333333333 0.16666666666667 0.00000000000000 39.948 > 1 85 Ar 0.16666666666667 0.50000000000000 0.00000000000000 39.948 > 1 86 Ar 0.50000000000000 0.50000000000000 0.00000000000000 39.948 > 1 87 Ar 0.83333333333333 0.50000000000000 0.00000000000000 39.948 > 1 88 Ar 0.16666666666667 0.83333333333333 0.00000000000000 39.948 > 1 89 Ar 0.50000000000000 0.83333333333333 0.00000000000000 39.948 > 1 90 Ar 0.83333333333333 0.83333333333333 0.00000000000000 39.948 > 1 91 Ar 0.16666666666667 0.16666666666667 0.33333333333333 39.948 > 1 92 Ar 0.50000000000000 0.16666666666667 0.33333333333333 39.948 > 1 93 Ar 0.83333333333333 0.16666666666667 0.33333333333333 39.948 > 1 94 Ar 0.16666666666667 0.50000000000000 0.33333333333333 39.948 > 1 95 Ar 0.50000000000000 0.50000000000000 0.33333333333333 39.948 > 1 96 Ar 0.83333333333333 0.50000000000000 0.33333333333333 39.948 > 1 97 Ar 0.16666666666667 0.83333333333333 0.33333333333333 39.948 > 1 98 Ar 0.50000000000000 0.83333333333333 0.33333333333333 39.948 > 1 99 Ar 0.83333333333333 0.83333333333333 0.33333333333333 39.948 > 1 100 Ar 0.16666666666667 0.16666666666667 0.66666666666667 39.948 > 1 101 Ar 0.50000000000000 0.16666666666667 0.66666666666667 39.948 > 1 102 Ar 0.83333333333333 0.16666666666667 0.66666666666667 39.948 > 1 103 Ar 0.16666666666667 0.50000000000000 0.66666666666667 39.948 > 1 104 Ar 0.50000000000000 0.50000000000000 0.66666666666667 39.948 > 1 105 Ar 0.83333333333333 0.50000000000000 0.66666666666667 39.948 > 1 106 Ar 0.16666666666667 0.83333333333333 0.66666666666667 39.948 > 1 107 Ar 0.50000000000000 0.83333333333333 0.66666666666667 39.948 > 1 108 Ar 0.83333333333333 0.83333333333333 0.66666666666667 39.948 > 1 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 1.5997120 0.0000000 0.0000000 0.0000000 1.5997120 0.0000000 0.0000000 0.0000000 1.5997120 -------------------------- Born effective charges -------------------------- 1 Ar 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 ---------------------------------------------------------------------------- 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] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000000 (yyz) 0.00000000 (yyz) 0.00000000 (yzy) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000000 (zz) 0.00000000 (zz) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 05:19:56]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 05:19:56]------------------------- Compiled with OpenMP support (max 128 threads). Running in phono3py.load mode. Python version 3.14.2 Spglib version 2.6.1 ----------------------------- General settings ----------------------------- Run mode: conductivity-RTA HDF5 data compression filter: gzip Crystal structure was read from "phono3py.yaml". Supercell (dim): [3 3 3] Primitive matrix: [0. 0.5 0.5] [0.5 0. 0.5] [0.5 0.5 0. ] Spacegroup: Fm-3m (225) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 0.000000000000000 2.811709859999999 2.811709859999999 b 2.811709859999999 0.000000000000000 2.811709859999999 c 2.811709859999999 2.811709859999999 0.000000000000000 Atomic positions (fractional): 1 Ar 0.00000000000000 0.00000000000000 0.00000000000000 39.948 -------------------------------- supercell --------------------------------- Lattice vectors: a 16.870259159999996 0.000000000000000 0.000000000000000 b 0.000000000000000 16.870259159999996 0.000000000000000 c 0.000000000000000 0.000000000000000 16.870259159999996 Atomic positions (fractional): 1 Ar 0.00000000000000 0.00000000000000 0.00000000000000 39.948 > 1 2 Ar 0.33333333333333 0.00000000000000 0.00000000000000 39.948 > 1 3 Ar 0.66666666666667 0.00000000000000 0.00000000000000 39.948 > 1 4 Ar 0.00000000000000 0.33333333333333 0.00000000000000 39.948 > 1 5 Ar 0.33333333333333 0.33333333333333 0.00000000000000 39.948 > 1 6 Ar 0.66666666666667 0.33333333333333 0.00000000000000 39.948 > 1 7 Ar 0.00000000000000 0.66666666666667 0.00000000000000 39.948 > 1 8 Ar 0.33333333333333 0.66666666666667 0.00000000000000 39.948 > 1 9 Ar 0.66666666666667 0.66666666666667 0.00000000000000 39.948 > 1 10 Ar 0.00000000000000 0.00000000000000 0.33333333333333 39.948 > 1 11 Ar 0.33333333333333 0.00000000000000 0.33333333333333 39.948 > 1 12 Ar 0.66666666666667 0.00000000000000 0.33333333333333 39.948 > 1 13 Ar 0.00000000000000 0.33333333333333 0.33333333333333 39.948 > 1 14 Ar 0.33333333333333 0.33333333333333 0.33333333333333 39.948 > 1 15 Ar 0.66666666666667 0.33333333333333 0.33333333333333 39.948 > 1 16 Ar 0.00000000000000 0.66666666666667 0.33333333333333 39.948 > 1 17 Ar 0.33333333333333 0.66666666666667 0.33333333333333 39.948 > 1 18 Ar 0.66666666666667 0.66666666666667 0.33333333333333 39.948 > 1 19 Ar 0.00000000000000 0.00000000000000 0.66666666666667 39.948 > 1 20 Ar 0.33333333333333 0.00000000000000 0.66666666666667 39.948 > 1 21 Ar 0.66666666666667 0.00000000000000 0.66666666666667 39.948 > 1 22 Ar 0.00000000000000 0.33333333333333 0.66666666666667 39.948 > 1 23 Ar 0.33333333333333 0.33333333333333 0.66666666666667 39.948 > 1 24 Ar 0.66666666666667 0.33333333333333 0.66666666666667 39.948 > 1 25 Ar 0.00000000000000 0.66666666666667 0.66666666666667 39.948 > 1 26 Ar 0.33333333333333 0.66666666666667 0.66666666666667 39.948 > 1 27 Ar 0.66666666666667 0.66666666666667 0.66666666666667 39.948 > 1 28 Ar 0.00000000000000 0.16666666666667 0.16666666666667 39.948 > 1 29 Ar 0.33333333333333 0.16666666666667 0.16666666666667 39.948 > 1 30 Ar 0.66666666666667 0.16666666666667 0.16666666666667 39.948 > 1 31 Ar 0.00000000000000 0.50000000000000 0.16666666666667 39.948 > 1 32 Ar 0.33333333333333 0.50000000000000 0.16666666666667 39.948 > 1 33 Ar 0.66666666666667 0.50000000000000 0.16666666666667 39.948 > 1 34 Ar 0.00000000000000 0.83333333333333 0.16666666666667 39.948 > 1 35 Ar 0.33333333333333 0.83333333333333 0.16666666666667 39.948 > 1 36 Ar 0.66666666666667 0.83333333333333 0.16666666666667 39.948 > 1 37 Ar 0.00000000000000 0.16666666666667 0.50000000000000 39.948 > 1 38 Ar 0.33333333333333 0.16666666666667 0.50000000000000 39.948 > 1 39 Ar 0.66666666666667 0.16666666666667 0.50000000000000 39.948 > 1 40 Ar 0.00000000000000 0.50000000000000 0.50000000000000 39.948 > 1 41 Ar 0.33333333333333 0.50000000000000 0.50000000000000 39.948 > 1 42 Ar 0.66666666666667 0.50000000000000 0.50000000000000 39.948 > 1 43 Ar 0.00000000000000 0.83333333333333 0.50000000000000 39.948 > 1 44 Ar 0.33333333333333 0.83333333333333 0.50000000000000 39.948 > 1 45 Ar 0.66666666666667 0.83333333333333 0.50000000000000 39.948 > 1 46 Ar 0.00000000000000 0.16666666666667 0.83333333333333 39.948 > 1 47 Ar 0.33333333333333 0.16666666666667 0.83333333333333 39.948 > 1 48 Ar 0.66666666666667 0.16666666666667 0.83333333333333 39.948 > 1 49 Ar 0.00000000000000 0.50000000000000 0.83333333333333 39.948 > 1 50 Ar 0.33333333333333 0.50000000000000 0.83333333333333 39.948 > 1 51 Ar 0.66666666666667 0.50000000000000 0.83333333333333 39.948 > 1 52 Ar 0.00000000000000 0.83333333333333 0.83333333333333 39.948 > 1 53 Ar 0.33333333333333 0.83333333333333 0.83333333333333 39.948 > 1 54 Ar 0.66666666666667 0.83333333333333 0.83333333333333 39.948 > 1 55 Ar 0.16666666666667 0.00000000000000 0.16666666666667 39.948 > 1 56 Ar 0.50000000000000 0.00000000000000 0.16666666666667 39.948 > 1 57 Ar 0.83333333333333 0.00000000000000 0.16666666666667 39.948 > 1 58 Ar 0.16666666666667 0.33333333333333 0.16666666666667 39.948 > 1 59 Ar 0.50000000000000 0.33333333333333 0.16666666666667 39.948 > 1 60 Ar 0.83333333333333 0.33333333333333 0.16666666666667 39.948 > 1 61 Ar 0.16666666666667 0.66666666666667 0.16666666666667 39.948 > 1 62 Ar 0.50000000000000 0.66666666666667 0.16666666666667 39.948 > 1 63 Ar 0.83333333333333 0.66666666666667 0.16666666666667 39.948 > 1 64 Ar 0.16666666666667 0.00000000000000 0.50000000000000 39.948 > 1 65 Ar 0.50000000000000 0.00000000000000 0.50000000000000 39.948 > 1 66 Ar 0.83333333333333 0.00000000000000 0.50000000000000 39.948 > 1 67 Ar 0.16666666666667 0.33333333333333 0.50000000000000 39.948 > 1 68 Ar 0.50000000000000 0.33333333333333 0.50000000000000 39.948 > 1 69 Ar 0.83333333333333 0.33333333333333 0.50000000000000 39.948 > 1 70 Ar 0.16666666666667 0.66666666666667 0.50000000000000 39.948 > 1 71 Ar 0.50000000000000 0.66666666666667 0.50000000000000 39.948 > 1 72 Ar 0.83333333333333 0.66666666666667 0.50000000000000 39.948 > 1 73 Ar 0.16666666666667 0.00000000000000 0.83333333333333 39.948 > 1 74 Ar 0.50000000000000 0.00000000000000 0.83333333333333 39.948 > 1 75 Ar 0.83333333333333 0.00000000000000 0.83333333333333 39.948 > 1 76 Ar 0.16666666666667 0.33333333333333 0.83333333333333 39.948 > 1 77 Ar 0.50000000000000 0.33333333333333 0.83333333333333 39.948 > 1 78 Ar 0.83333333333333 0.33333333333333 0.83333333333333 39.948 > 1 79 Ar 0.16666666666667 0.66666666666667 0.83333333333333 39.948 > 1 80 Ar 0.50000000000000 0.66666666666667 0.83333333333333 39.948 > 1 81 Ar 0.83333333333333 0.66666666666667 0.83333333333333 39.948 > 1 82 Ar 0.16666666666667 0.16666666666667 0.00000000000000 39.948 > 1 83 Ar 0.50000000000000 0.16666666666667 0.00000000000000 39.948 > 1 84 Ar 0.83333333333333 0.16666666666667 0.00000000000000 39.948 > 1 85 Ar 0.16666666666667 0.50000000000000 0.00000000000000 39.948 > 1 86 Ar 0.50000000000000 0.50000000000000 0.00000000000000 39.948 > 1 87 Ar 0.83333333333333 0.50000000000000 0.00000000000000 39.948 > 1 88 Ar 0.16666666666667 0.83333333333333 0.00000000000000 39.948 > 1 89 Ar 0.50000000000000 0.83333333333333 0.00000000000000 39.948 > 1 90 Ar 0.83333333333333 0.83333333333333 0.00000000000000 39.948 > 1 91 Ar 0.16666666666667 0.16666666666667 0.33333333333333 39.948 > 1 92 Ar 0.50000000000000 0.16666666666667 0.33333333333333 39.948 > 1 93 Ar 0.83333333333333 0.16666666666667 0.33333333333333 39.948 > 1 94 Ar 0.16666666666667 0.50000000000000 0.33333333333333 39.948 > 1 95 Ar 0.50000000000000 0.50000000000000 0.33333333333333 39.948 > 1 96 Ar 0.83333333333333 0.50000000000000 0.33333333333333 39.948 > 1 97 Ar 0.16666666666667 0.83333333333333 0.33333333333333 39.948 > 1 98 Ar 0.50000000000000 0.83333333333333 0.33333333333333 39.948 > 1 99 Ar 0.83333333333333 0.83333333333333 0.33333333333333 39.948 > 1 100 Ar 0.16666666666667 0.16666666666667 0.66666666666667 39.948 > 1 101 Ar 0.50000000000000 0.16666666666667 0.66666666666667 39.948 > 1 102 Ar 0.83333333333333 0.16666666666667 0.66666666666667 39.948 > 1 103 Ar 0.16666666666667 0.50000000000000 0.66666666666667 39.948 > 1 104 Ar 0.50000000000000 0.50000000000000 0.66666666666667 39.948 > 1 105 Ar 0.83333333333333 0.50000000000000 0.66666666666667 39.948 > 1 106 Ar 0.16666666666667 0.83333333333333 0.66666666666667 39.948 > 1 107 Ar 0.50000000000000 0.83333333333333 0.66666666666667 39.948 > 1 108 Ar 0.83333333333333 0.83333333333333 0.66666666666667 39.948 > 1 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 1.5997120 0.0000000 0.0000000 0.0000000 1.5997120 0.0000000 0.0000000 0.0000000 1.5997120 -------------------------- Born effective charges -------------------------- 1 Ar 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000000 (yyz) 0.00000000 (yyz) 0.00000000 (yzy) 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 ... [ 15 15 15 ] 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.17, Number of G-points: 307, Lambda: 0.15 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/120) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 120 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 ======================= Grid point 1 (2/120) ======================= q-point: ( 0.07 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.192 ( -5.040 5.040 5.040) 8.729 0.192 ( -5.040 5.040 5.040) 8.729 0.281 ( -7.608 7.608 7.608) 13.178 ======================= Grid point 2 (3/120) ======================= q-point: ( 0.13 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.350 ( -3.765 3.765 3.765) 6.520 0.350 ( -3.765 3.765 3.765) 6.520 0.533 ( -6.498 6.498 6.498) 11.255 ======================= Grid point 3 (4/120) ======================= q-point: ( 0.20 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.460 ( -2.478 2.478 2.478) 4.291 0.460 ( -2.478 2.478 2.478) 4.291 0.742 ( -5.256 5.256 5.256) 9.104 ======================= Grid point 4 (5/120) ======================= q-point: ( 0.27 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.532 ( -1.672 1.672 1.672) 2.896 0.532 ( -1.672 1.672 1.672) 2.896 0.909 ( -4.170 4.170 4.170) 7.222 ======================= Grid point 5 (6/120) ======================= q-point: ( 0.33 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.583 ( -1.233 1.233 1.233) 2.136 0.583 ( -1.233 1.233 1.233) 2.136 1.039 ( -3.134 3.134 3.134) 5.429 ======================= Grid point 6 (7/120) ======================= q-point: ( 0.40 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.619 ( -0.805 0.805 0.805) 1.395 0.619 ( -0.805 0.805 0.805) 1.395 1.131 ( -1.971 1.971 1.971) 3.414 ======================= Grid point 7 (8/120) ======================= q-point: ( 0.47 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 344 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.639 ( -0.280 0.280 0.280) 0.485 0.639 ( -0.280 0.280 0.280) 0.485 1.178 ( -0.672 0.672 0.672) 1.165 ======================= Grid point 16 (9/120) ======================= q-point: ( 0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 288 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.238 ( -0.000 0.000 9.464) 9.464 0.238 ( -0.000 0.000 9.464) 9.464 0.300 ( -0.000 0.000 12.101) 12.101 ======================= Grid point 17 (10/120) ======================= q-point: ( 0.13 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.351 ( -0.985 0.985 7.506) 7.634 0.364 ( -1.493 1.493 7.401) 7.696 0.497 ( -4.233 4.233 9.400) 11.145 ======================= Grid point 18 (11/120) ======================= q-point: ( 0.20 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.452 ( -1.101 1.101 5.085) 5.318 0.470 ( -1.178 1.178 5.123) 5.387 0.697 ( -4.387 4.387 6.840) 9.235 ======================= Grid point 19 (12/120) ======================= q-point: ( 0.27 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.525 ( -0.807 0.807 3.633) 3.808 0.543 ( -0.701 0.701 3.562) 3.697 0.865 ( -3.933 3.933 4.905) 7.416 ======================= Grid point 20 (13/120) ======================= q-point: ( 0.33 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.581 ( -0.568 0.568 3.052) 3.156 0.593 ( -0.450 0.450 2.646) 2.721 1.000 ( -3.307 3.307 3.388) 5.775 ======================= Grid point 21 (14/120) ======================= q-point: ( 0.40 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.625 ( -0.221 0.221 2.789) 2.806 0.628 ( -0.217 0.217 1.922) 1.946 1.101 ( -2.442 2.442 1.979) 3.980 ======================= Grid point 22 (15/120) ======================= q-point: ( 0.47 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.647 ( 0.136 -0.136 1.140) 1.156 0.655 ( 0.332 -0.332 2.511) 2.555 1.161 ( -1.330 1.330 0.560) 1.962 ======================= Grid point 23 (16/120) ======================= q-point: (-0.47 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.648 ( 0.555 -0.555 0.333) 0.852 0.668 ( 0.992 -0.992 2.196) 2.606 1.176 ( -0.068 0.068 -0.807) 0.812 ======================= Grid point 24 (17/120) ======================= q-point: (-0.40 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.629 ( 0.968 -0.968 -0.414) 1.430 0.660 ( 1.659 -1.659 1.874) 3.003 1.145 ( 1.243 -1.243 -2.024) 2.681 ======================= Grid point 25 (18/120) ======================= q-point: (-0.33 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.593 ( 1.305 -1.305 -0.991) 2.095 0.634 ( 2.237 -2.237 1.594) 3.543 1.071 ( 2.510 -2.510 -2.975) 4.632 ======================= Grid point 26 (19/120) ======================= q-point: (-0.27 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.545 ( 1.603 -1.603 -1.315) 2.621 0.593 ( 2.731 -2.731 1.358) 4.094 0.958 ( 3.699 -3.699 -3.545) 6.319 ======================= Grid point 27 (20/120) ======================= q-point: (-0.20 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.484 ( 2.216 -2.216 -1.492) 3.471 0.535 ( 3.434 -3.434 0.990) 4.956 0.812 ( 4.971 -4.971 -3.679) 7.935 ======================= Grid point 28 (21/120) ======================= q-point: (-0.13 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.400 ( 3.494 -3.494 -1.451) 5.150 0.449 ( 4.627 -4.627 0.415) 6.557 0.635 ( 6.564 -6.564 -3.072) 9.778 ======================= Grid point 29 (22/120) ======================= q-point: (-0.07 0.07 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 512 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.287 ( 5.076 -5.076 0.000) 7.178 0.324 ( 5.975 -5.975 0.000) 8.450 0.436 ( 8.257 -8.257 0.000) 11.677 ======================= Grid point 32 (23/120) ======================= q-point: ( 0.13 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 288 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.438 ( -0.000 0.000 7.252) 7.252 0.438 ( -0.000 0.000 7.252) 7.252 0.563 ( -0.000 0.000 9.943) 9.943 ======================= Grid point 33 (24/120) ======================= q-point: ( 0.20 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.513 ( 0.194 -0.194 5.876) 5.883 0.520 ( -0.260 0.260 5.600) 5.612 0.703 ( -2.338 2.338 7.816) 8.487 ======================= Grid point 34 (25/120) ======================= q-point: ( 0.27 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.571 ( 0.143 -0.143 4.747) 4.751 0.580 ( -0.076 0.076 4.079) 4.081 0.845 ( -2.855 2.855 5.563) 6.874 ======================= Grid point 35 (26/120) ======================= q-point: ( 0.33 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.621 ( 0.089 -0.089 4.317) 4.319 0.623 ( 0.043 -0.043 3.134) 3.134 0.968 ( -2.831 2.831 3.721) 5.466 ======================= Grid point 36 (27/120) ======================= q-point: ( 0.40 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.654 ( 0.168 -0.168 2.450) 2.461 0.669 ( 0.224 -0.224 4.234) 4.246 1.066 ( -2.438 2.438 2.084) 4.029 ======================= Grid point 37 (28/120) ======================= q-point: ( 0.47 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.672 ( 0.447 -0.447 1.718) 1.831 0.708 ( 0.666 -0.666 4.102) 4.209 1.131 ( -1.688 1.688 0.511) 2.442 ======================= Grid point 38 (29/120) ======================= q-point: ( 0.53 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.672 ( 0.833 -0.833 0.946) 1.511 0.733 ( 1.304 -1.304 3.850) 4.269 1.157 ( -0.685 0.685 -0.926) 1.341 ======================= Grid point 39 (30/120) ======================= q-point: (-0.40 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.655 ( 1.239 -1.239 0.233) 1.767 0.737 ( 1.996 -1.996 3.485) 4.485 1.141 ( 0.466 -0.466 -2.087) 2.189 ======================= Grid point 40 (31/120) ======================= q-point: (-0.33 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.620 ( 1.582 -1.582 -0.321) 2.260 0.721 ( 2.613 -2.613 2.989) 4.753 1.086 ( 1.676 -1.676 -2.813) 3.678 ======================= Grid point 41 (32/120) ======================= q-point: (-0.27 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.574 ( 1.781 -1.781 -0.578) 2.584 0.685 ( 3.035 -3.035 2.328) 4.883 0.998 ( 2.858 -2.858 -2.876) 4.960 ======================= Grid point 42 (33/120) ======================= q-point: (-0.20 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.523 ( 2.006 -2.006 -0.486) 2.878 0.632 ( 3.354 -3.354 1.344) 4.929 0.886 ( 4.061 -4.061 -2.018) 6.087 ======================= Grid point 43 (34/120) ======================= q-point: (-0.13 0.13 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 512 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.466 ( 2.633 -2.633 0.000) 3.724 0.554 ( 3.884 -3.884 0.000) 5.493 0.761 ( 5.409 -5.409 0.000) 7.650 ======================= Grid point 48 (35/120) ======================= q-point: ( 0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 288 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.580 ( -0.000 0.000 4.783) 4.783 0.580 ( -0.000 0.000 4.783) 4.783 0.769 ( -0.000 0.000 7.392) 7.392 ======================= Grid point 49 (36/120) ======================= q-point: ( 0.27 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.629 ( 0.134 -0.134 3.855) 3.860 0.629 ( 0.136 -0.136 3.908) 3.912 0.864 ( -1.254 1.254 5.813) 6.077 ======================= Grid point 50 (37/120) ======================= q-point: ( 0.33 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.664 ( 0.248 -0.248 3.054) 3.074 0.674 ( -0.020 0.020 3.982) 3.982 0.959 ( -1.736 1.736 4.078) 4.760 ======================= Grid point 51 (38/120) ======================= q-point: ( 0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.691 ( 0.357 -0.357 2.569) 2.618 0.725 ( -0.037 0.037 4.297) 4.297 1.040 ( -1.803 1.803 2.357) 3.472 ======================= Grid point 52 (39/120) ======================= q-point: ( 0.47 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.707 ( 0.604 -0.604 1.973) 2.150 0.774 ( 0.340 -0.340 4.525) 4.551 1.098 ( -1.510 1.510 0.603) 2.219 ======================= Grid point 53 (40/120) ======================= q-point: ( 0.53 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.707 ( 0.965 -0.965 1.277) 1.869 0.812 ( 1.034 -1.034 4.530) 4.760 1.125 ( -0.944 0.944 -1.015) 1.677 ======================= Grid point 54 (41/120) ======================= q-point: (-0.40 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.691 ( 1.355 -1.355 0.627) 2.017 0.830 ( 1.864 -1.864 4.239) 4.992 1.119 ( -0.184 0.184 -2.197) 2.213 ======================= Grid point 55 (42/120) ======================= q-point: (-0.33 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.659 ( 1.718 -1.718 0.121) 2.433 0.823 ( 2.654 -2.654 3.440) 5.091 1.083 ( 0.750 -0.750 -2.591) 2.799 ======================= Grid point 56 (43/120) ======================= q-point: (-0.27 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.615 ( 1.934 -1.934 -0.122) 2.737 0.785 ( 3.169 -3.169 1.935) 4.882 1.025 ( 1.822 -1.822 -1.826) 3.158 ======================= Grid point 57 (44/120) ======================= q-point: (-0.20 0.20 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 512 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.567 ( 1.952 -1.952 0.000) 2.760 0.720 ( 3.312 -3.312 0.000) 4.684 0.957 ( 2.957 -2.957 0.000) 4.182 ======================= Grid point 64 (45/120) ======================= q-point: ( 0.27 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 288 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.671 ( -0.000 0.000 3.116) 3.116 0.671 ( -0.000 0.000 3.116) 3.116 0.918 ( -0.000 0.000 5.310) 5.310 ======================= Grid point 65 (46/120) ======================= q-point: ( 0.33 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.704 ( 0.163 -0.163 2.705) 2.714 0.710 ( -0.352 0.352 2.971) 3.012 0.983 ( -0.571 0.571 4.240) 4.316 ======================= Grid point 66 (47/120) ======================= q-point: ( 0.40 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.729 ( 0.313 -0.313 2.407) 2.447 0.761 ( -0.747 0.747 3.329) 3.492 1.041 ( -0.729 0.729 2.835) 3.016 ======================= Grid point 67 (48/120) ======================= q-point: ( 0.47 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.744 ( 0.573 -0.573 1.931) 2.094 0.821 ( -0.695 0.695 3.741) 3.868 1.081 ( -0.622 0.622 1.102) 1.409 ======================= Grid point 68 (49/120) ======================= q-point: ( 0.53 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.746 ( 0.930 -0.930 1.289) 1.841 0.878 ( -0.088 0.088 4.103) 4.105 1.095 ( -0.416 0.416 -0.799) 0.992 ======================= Grid point 69 (50/120) ======================= q-point: (-0.40 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.731 ( 1.291 -1.291 0.678) 1.948 0.918 ( 1.032 -1.032 4.272) 4.515 1.084 ( -0.317 0.317 -2.400) 2.442 ======================= Grid point 70 (51/120) ======================= q-point: (-0.33 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.701 ( 1.632 -1.632 0.220) 2.318 0.922 ( 2.416 -2.416 3.030) 4.567 1.060 ( -0.225 0.225 -2.304) 2.326 ======================= Grid point 71 (52/120) ======================= q-point: (-0.27 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 512 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.660 ( 1.892 -1.892 0.000) 2.675 0.873 ( 3.088 -3.088 0.000) 4.367 1.043 ( 0.629 -0.629 0.000) 0.890 ======================= Grid point 80 (53/120) ======================= q-point: ( 0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 288 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.736 ( -0.000 0.000 2.438) 2.438 0.736 ( -0.000 0.000 2.438) 2.438 1.026 ( -0.000 0.000 3.879) 3.879 ======================= Grid point 81 (54/120) ======================= q-point: ( 0.40 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.762 ( 0.171 -0.171 2.195) 2.209 0.774 ( -0.814 0.814 2.397) 2.659 1.068 ( -0.043 0.043 3.020) 3.020 ======================= Grid point 82 (55/120) ======================= q-point: ( 0.47 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.778 ( 0.438 -0.438 1.733) 1.841 0.830 ( -1.447 1.447 2.412) 3.163 1.095 ( 0.251 -0.251 1.738) 1.774 ======================= Grid point 83 (56/120) ======================= q-point: ( 0.53 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.780 ( 0.791 -0.791 1.075) 1.551 0.896 ( -1.624 1.624 2.393) 3.316 1.095 ( 0.726 -0.726 0.201) 1.046 ======================= Grid point 84 (57/120) ======================= q-point: ( 0.60 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.766 ( 1.123 -1.123 0.450) 1.650 0.961 ( -1.346 1.346 2.527) 3.164 1.065 ( 1.196 -1.196 -1.389) 2.189 ======================= Grid point 85 (58/120) ======================= q-point: (-0.33 0.33 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 512 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.739 ( 1.402 -1.402 0.000) 1.982 1.004 ( 2.370 -2.370 0.000) 3.352 1.020 ( -1.454 1.454 0.000) 2.056 ======================= Grid point 96 (59/120) ======================= q-point: ( 0.40 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 288 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.788 ( -0.000 0.000 1.906) 1.906 0.788 ( -0.000 0.000 1.906) 1.906 1.103 ( -0.000 0.000 2.583) 2.583 ======================= Grid point 97 (60/120) ======================= q-point: ( 0.47 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.805 ( 0.253 -0.253 1.373) 1.419 0.821 ( -1.052 1.052 1.479) 2.097 1.125 ( 0.359 -0.359 1.664) 1.740 ======================= Grid point 98 (61/120) ======================= q-point: ( 0.53 0.40 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.807 ( 0.584 -0.584 0.656) 1.055 0.870 ( -1.853 1.853 0.908) 2.773 1.123 ( 0.928 -0.928 0.596) 1.442 ======================= Grid point 99 (62/120) ======================= q-point: (-0.40 -0.60 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 512 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.793 ( 0.892 -0.892 0.000) 1.262 0.926 ( -2.227 2.227 0.000) 3.149 1.096 ( 1.509 -1.509 0.000) 2.133 ======================= Grid point 113 (63/120) ======================= q-point: ( 0.47 0.47 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 288 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.821 ( -0.000 0.000 0.762) 0.762 0.821 ( -0.000 0.000 0.762) 0.762 1.145 ( -0.000 0.000 0.938) 0.938 ======================= Grid point 114 (64/120) ======================= q-point: (-0.47 -0.53 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 512 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.822 ( 0.306 -0.306 0.000) 0.432 0.839 ( -1.111 1.111 0.000) 1.571 1.145 ( 0.517 -0.517 0.000) 0.731 ======================= Grid point 262 (65/120) ======================= q-point: ( 0.20 0.13 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.430 ( -0.000 0.261 6.150) 6.155 0.473 ( -0.000 2.733 6.006) 6.599 0.635 ( -0.000 4.870 8.031) 9.392 ======================= Grid point 263 (66/120) ======================= q-point: ( 0.27 0.13 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.499 ( -0.708 0.010 4.461) 4.517 0.560 ( 1.015 2.536 4.531) 5.291 0.786 ( -2.194 4.156 5.968) 7.596 ======================= Grid point 264 (67/120) ======================= q-point: ( 0.33 0.13 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.554 ( -0.903 -0.114 3.476) 3.594 0.620 ( 1.746 2.347 3.652) 4.680 0.922 ( -2.971 3.382 4.194) 6.153 ======================= Grid point 265 (68/120) ======================= q-point: ( 0.40 0.13 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.602 ( -0.825 -0.084 2.960) 3.074 0.664 ( 2.173 2.122 3.215) 4.423 1.035 ( -2.981 2.598 2.679) 4.777 ======================= Grid point 266 (69/120) ======================= q-point: ( 0.47 0.13 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.640 ( -0.365 -0.088 2.365) 2.395 0.695 ( 2.477 1.562 2.921) 4.136 1.116 ( -2.416 1.654 1.222) 3.173 ======================= Grid point 267 (70/120) ======================= q-point: ( 0.53 0.13 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.662 ( 0.467 -0.056 1.482) 1.555 0.711 ( 2.662 0.587 2.747) 3.869 1.158 ( -1.434 0.567 -0.184) 1.553 ======================= Grid point 278 (71/120) ======================= q-point: ( 0.27 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.560 ( -0.000 -1.203 4.750) 4.900 0.598 ( -0.000 1.464 4.471) 4.705 0.806 ( -0.000 2.747 6.358) 6.926 ======================= Grid point 279 (72/120) ======================= q-point: ( 0.33 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.601 ( -0.460 -1.137 4.009) 4.192 0.658 ( 1.126 1.828 3.735) 4.308 0.912 ( -1.513 2.497 4.615) 5.461 ======================= Grid point 280 (73/120) ======================= q-point: ( 0.40 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.640 ( -0.554 -0.893 3.563) 3.715 0.706 ( 1.836 1.936 3.587) 4.471 1.008 ( -2.230 2.127 2.946) 4.264 ======================= Grid point 281 (74/120) ======================= q-point: ( 0.47 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.674 ( -0.102 -0.629 2.904) 2.973 0.746 ( 2.163 1.445 3.719) 4.539 1.082 ( -2.256 1.592 1.305) 3.054 ======================= Grid point 282 (75/120) ======================= q-point: ( 0.53 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.693 ( 0.776 -0.409 1.950) 2.138 0.776 ( 2.312 0.391 3.858) 4.515 1.127 ( -1.728 0.878 -0.256) 1.955 ======================= Grid point 283 (76/120) ======================= q-point: (-0.40 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.692 ( 1.690 -0.318 1.007) 1.993 0.790 ( 2.520 -0.806 3.807) 4.636 1.137 ( -0.802 0.063 -1.551) 1.747 ======================= Grid point 284 (77/120) ======================= q-point: (-0.33 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.670 ( 2.446 -0.366 0.255) 2.487 0.786 ( 2.862 -1.846 3.444) 4.843 1.111 ( 0.411 -0.803 -2.387) 2.552 ======================= Grid point 285 (78/120) ======================= q-point: (-0.27 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.633 ( 2.981 -0.399 -0.212) 3.015 0.760 ( 3.220 -2.547 2.699) 4.914 1.054 ( 1.788 -1.666 -2.535) 3.521 ======================= Grid point 286 (79/120) ======================= q-point: (-0.20 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.588 ( 3.296 -0.252 -0.294) 3.319 0.714 ( 3.479 -2.777 1.521) 4.704 0.973 ( 3.204 -2.489 -1.741) 4.415 ======================= Grid point 287 (80/120) ======================= q-point: (-0.13 0.20 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.543 ( 3.681 -0.043 0.000) 3.681 0.647 ( 3.824 -2.664 0.000) 4.660 0.881 ( 4.624 -3.337 0.000) 5.703 ======================= Grid point 294 (81/120) ======================= q-point: ( 0.33 0.27 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.657 ( -0.000 -1.056 3.494) 3.650 0.687 ( -0.000 1.314 3.192) 3.452 0.936 ( -0.000 1.428 4.684) 4.897 ======================= Grid point 295 (82/120) ======================= q-point: ( 0.40 0.27 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.688 ( -0.220 -0.857 3.325) 3.441 0.739 ( 0.626 1.836 3.112) 3.667 1.006 ( -0.898 1.239 3.285) 3.624 ======================= Grid point 296 (83/120) ======================= q-point: ( 0.47 0.27 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.717 ( 0.231 -0.564 2.835) 2.900 0.790 ( 0.849 1.653 3.502) 3.965 1.062 ( -1.349 0.925 1.637) 2.314 ======================= Grid point 297 (84/120) ======================= q-point: ( 0.53 0.27 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.733 ( 1.044 -0.380 2.009) 2.295 0.838 ( 1.099 0.743 3.950) 4.166 1.096 ( -1.360 0.567 -0.143) 1.480 ======================= Grid point 298 (85/120) ======================= q-point: ( 0.60 0.27 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.730 ( 1.808 -0.377 1.175) 2.189 0.872 ( 1.619 -0.525 4.150) 4.485 1.104 ( -1.016 0.260 -1.689) 1.988 ======================= Grid point 299 (86/120) ======================= q-point: (-0.33 0.27 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.710 ( 2.434 -0.496 0.511) 2.536 0.882 ( 2.310 -1.854 3.721) 4.756 1.088 ( -0.329 -0.031 -2.437) 2.459 ======================= Grid point 300 (87/120) ======================= q-point: (-0.27 0.27 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.674 ( 2.924 -0.626 0.085) 2.991 0.859 ( 2.916 -2.838 2.179) 4.616 1.056 ( 0.746 -0.557 -1.776) 2.005 ======================= Grid point 301 (88/120) ======================= q-point: (-0.20 0.27 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.630 ( 3.199 -0.568 0.000) 3.249 0.799 ( 3.267 -3.085 0.000) 4.493 1.017 ( 2.031 -1.467 0.000) 2.505 ======================= Grid point 310 (89/120) ======================= q-point: ( 0.40 0.33 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.731 ( -0.000 -0.414 2.799) 2.830 0.754 ( -0.000 1.483 2.511) 2.916 1.032 ( -0.000 0.447 3.391) 3.421 ======================= Grid point 311 (90/120) ======================= q-point: ( 0.47 0.33 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.757 ( 0.415 -0.214 2.520) 2.563 0.807 ( -0.433 1.901 2.564) 3.221 1.070 ( -0.267 0.035 2.146) 2.163 ======================= Grid point 312 (91/120) ======================= q-point: ( 0.53 0.33 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.772 ( 1.168 -0.134 1.794) 2.145 0.868 ( -0.701 1.665 2.875) 3.395 1.087 ( -0.249 -0.437 0.528) 0.729 ======================= Grid point 313 (92/120) ======================= q-point: ( 0.60 0.33 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.769 ( 1.805 -0.229 1.014) 2.083 0.927 ( -0.318 0.844 3.349) 3.468 1.077 ( -0.140 -0.666 -1.315) 1.480 ======================= Grid point 314 (93/120) ======================= q-point: (-0.33 0.33 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.749 ( 2.287 -0.392 0.394) 2.353 0.967 ( 1.101 -1.099 3.678) 3.993 1.049 ( -0.501 0.095 -2.817) 2.862 ======================= Grid point 315 (94/120) ======================= q-point: (-0.27 0.33 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.716 ( 2.670 -0.568 0.000) 2.730 0.941 ( 2.509 -2.901 -0.000) 3.836 1.045 ( -0.439 0.648 0.000) 0.783 ======================= Grid point 326 (95/120) ======================= q-point: ( 0.47 0.40 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.789 ( -0.000 0.062 2.071) 2.072 0.806 ( -0.000 1.496 1.859) 2.386 1.098 ( -0.000 -0.396 2.202) 2.237 ======================= Grid point 327 (96/120) ======================= q-point: ( 0.53 0.40 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.805 ( 1.091 0.141 1.408) 1.787 0.857 ( -1.440 1.923 1.555) 2.861 1.108 ( 0.290 -1.034 1.079) 1.522 ======================= Grid point 328 (97/120) ======================= q-point: ( 0.60 0.40 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.801 ( 1.680 -0.010 0.620) 1.791 0.918 ( -1.998 2.120 1.206) 3.153 1.091 ( 0.840 -1.694 -0.008) 1.891 ======================= Grid point 329 (98/120) ======================= q-point: (-0.33 0.40 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.781 ( 2.047 -0.208 0.000) 2.057 0.978 ( -2.056 2.090 0.000) 2.932 1.051 ( 1.540 -2.297 -0.000) 2.765 ======================= Grid point 342 (99/120) ======================= q-point: ( 0.53 0.47 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.824 ( -0.000 0.243 0.795) 0.832 0.838 ( -0.000 1.389 0.717) 1.563 1.134 ( -0.000 -0.940 0.786) 1.225 ======================= Grid point 343 (100/120) ======================= q-point: (-0.40 -0.53 0.07) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.822 ( 1.349 0.211 0.000) 1.366 0.876 ( -1.825 1.942 0.000) 2.665 1.121 ( 0.527 -1.490 -0.000) 1.580 ======================= Grid point 523 (101/120) ======================= q-point: ( 0.40 0.27 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.633 ( -0.000 -0.715 3.786) 3.853 0.732 ( -0.000 2.434 3.302) 4.102 0.977 ( -0.000 1.824 3.428) 3.883 ======================= Grid point 524 (102/120) ======================= q-point: ( 0.47 0.27 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.672 ( -0.062 -0.369 3.584) 3.604 0.793 ( 1.357 2.624 3.134) 4.307 1.035 ( -1.329 1.153 2.009) 2.671 ======================= Grid point 525 (103/120) ======================= q-point: ( 0.53 0.27 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.707 ( 0.372 -0.210 3.046) 3.076 0.835 ( 2.424 2.029 3.081) 4.414 1.081 ( -1.984 0.592 0.454) 2.119 ======================= Grid point 526 (104/120) ======================= q-point: ( 0.60 0.27 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.727 ( 1.196 -0.084 2.175) 2.484 0.856 ( 2.958 0.760 3.077) 4.335 1.105 ( -1.812 0.173 -0.927) 2.043 ======================= Grid point 539 (105/120) ======================= q-point: ( 0.47 0.33 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.719 ( -0.000 -0.487 3.395) 3.430 0.803 ( -0.000 2.542 2.627) 3.655 1.044 ( -0.000 0.497 2.227) 2.282 ======================= Grid point 540 (106/120) ======================= q-point: ( 0.53 0.33 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.751 ( 0.436 -0.275 3.009) 3.053 0.860 ( 0.597 2.511 2.503) 3.595 1.068 ( -0.744 -0.187 0.810) 1.115 ======================= Grid point 541 (107/120) ======================= q-point: ( 0.60 0.33 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.770 ( 1.350 -0.091 2.139) 2.531 0.906 ( 1.148 1.446 2.851) 3.396 1.076 ( -1.181 -0.447 -0.820) 1.506 ======================= Grid point 542 (108/120) ======================= q-point: (-0.33 0.33 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.767 ( 2.282 0.069 1.169) 2.564 0.931 ( 1.764 -0.530 3.139) 3.640 1.069 ( -1.184 0.013 -2.039) 2.358 ======================= Grid point 543 (109/120) ======================= q-point: (-0.27 0.33 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.746 ( 2.997 0.113 0.405) 3.026 0.920 ( 2.241 -2.448 2.013) 3.882 1.060 ( -0.423 0.515 -1.592) 1.726 ======================= Grid point 544 (110/120) ======================= q-point: (-0.20 0.33 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.711 ( 3.442 0.131 0.000) 3.445 0.870 ( 2.630 -3.000 -0.000) 3.990 1.051 ( 0.852 -0.080 0.000) 0.856 ======================= Grid point 555 (111/120) ======================= q-point: ( 0.53 0.40 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.790 ( -0.000 -0.036 2.469) 2.469 0.855 ( -0.000 2.519 1.754) 3.070 1.084 ( 0.000 -0.819 1.183) 1.439 ======================= Grid point 556 (112/120) ======================= q-point: ( 0.60 0.40 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.808 ( 1.241 0.124 1.708) 2.115 0.909 ( -0.796 2.428 1.572) 3.000 1.077 ( 0.030 -1.549 0.014) 1.549 ======================= Grid point 557 (113/120) ======================= q-point: ( 0.67 0.40 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.804 ( 2.206 0.231 0.740) 2.338 0.966 ( -0.957 1.688 1.990) 2.779 1.047 ( 0.271 -1.819 -1.322) 2.265 ======================= Grid point 558 (114/120) ======================= q-point: (-0.27 0.40 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.781 ( 2.735 0.195 0.000) 2.742 0.990 ( 1.585 -2.800 -0.000) 3.218 1.023 ( -1.315 1.730 0.000) 2.173 ======================= Grid point 571 (115/120) ======================= q-point: ( 0.60 0.47 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.831 ( 0.000 0.296 0.925) 0.971 0.884 ( -0.000 2.396 0.625) 2.476 1.102 ( -0.000 -1.773 0.354) 1.808 ======================= Grid point 572 (116/120) ======================= q-point: (-0.33 -0.53 0.13) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.829 ( 1.696 0.349 0.000) 1.731 0.929 ( -1.692 2.504 0.000) 3.022 1.076 ( 0.533 -2.282 -0.000) 2.343 ======================= Grid point 785 (117/120) ======================= q-point: ( 0.60 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.787 ( -0.000 -0.214 2.811) 2.819 0.921 ( -0.000 2.921 1.654) 3.356 1.059 ( 0.000 -1.276 -0.131) 1.283 ======================= Grid point 786 (118/120) ======================= q-point: ( 0.67 0.40 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 1688 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.809 ( 1.005 -0.110 2.102) 2.332 0.965 ( 1.019 1.955 1.559) 2.699 1.039 ( -1.050 -1.352 -1.214) 2.099 ======================= Grid point 801 (119/120) ======================= q-point: ( 0.67 0.47 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.836 ( 0.000 0.114 1.187) 1.192 0.947 ( 0.000 2.849 0.557) 2.903 1.052 ( 0.000 -2.383 -0.293) 2.401 ======================= Grid point 802 (120/120) ======================= q-point: (-0.27 0.47 0.20) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 3.51e-05 Number of triplets: 904 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.836 ( 1.520 0.174 0.000) 1.530 0.991 ( -1.048 2.654 0.000) 2.854 1.015 ( 0.535 -2.775 -0.000) 2.826 =================== 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/10125 10.0 4.598 4.598 4.598 0.000 -0.000 -0.000 3/10125 20.0 3.864 3.864 3.864 0.000 -0.000 -0.000 3/10125 30.0 2.969 2.969 2.969 0.000 -0.000 -0.000 3/10125 40.0 2.358 2.358 2.358 0.000 -0.000 -0.000 3/10125 50.0 1.942 1.942 1.942 0.000 -0.000 -0.000 3/10125 60.0 1.646 1.646 1.646 0.000 -0.000 -0.000 3/10125 70.0 1.427 1.427 1.427 0.000 -0.000 -0.000 3/10125 80.0 1.258 1.258 1.258 0.000 -0.000 -0.000 3/10125 90.0 1.124 1.124 1.124 0.000 -0.000 -0.000 3/10125 100.0 1.016 1.016 1.016 0.000 -0.000 -0.000 3/10125 110.0 0.927 0.927 0.927 0.000 -0.000 -0.000 3/10125 120.0 0.852 0.852 0.852 0.000 -0.000 -0.000 3/10125 130.0 0.788 0.788 0.788 0.000 -0.000 -0.000 3/10125 140.0 0.733 0.733 0.733 0.000 -0.000 -0.000 3/10125 150.0 0.685 0.685 0.685 0.000 -0.000 -0.000 3/10125 160.0 0.643 0.643 0.643 0.000 -0.000 -0.000 3/10125 170.0 0.606 0.606 0.606 0.000 -0.000 -0.000 3/10125 180.0 0.573 0.573 0.573 0.000 -0.000 -0.000 3/10125 190.0 0.543 0.543 0.543 0.000 -0.000 -0.000 3/10125 200.0 0.516 0.516 0.516 0.000 -0.000 -0.000 3/10125 210.0 0.492 0.492 0.492 0.000 -0.000 -0.000 3/10125 220.0 0.470 0.470 0.470 0.000 -0.000 -0.000 3/10125 230.0 0.450 0.450 0.450 0.000 -0.000 -0.000 3/10125 240.0 0.431 0.431 0.431 0.000 -0.000 -0.000 3/10125 250.0 0.414 0.414 0.414 0.000 -0.000 -0.000 3/10125 260.0 0.398 0.398 0.398 0.000 -0.000 -0.000 3/10125 270.0 0.384 0.384 0.384 0.000 -0.000 -0.000 3/10125 280.0 0.370 0.370 0.370 0.000 -0.000 -0.000 3/10125 290.0 0.358 0.358 0.358 0.000 -0.000 -0.000 3/10125 300.0 0.346 0.346 0.346 0.000 -0.000 -0.000 3/10125 310.0 0.335 0.335 0.335 0.000 -0.000 -0.000 3/10125 320.0 0.324 0.324 0.324 0.000 -0.000 -0.000 3/10125 330.0 0.314 0.314 0.314 0.000 -0.000 -0.000 3/10125 340.0 0.305 0.305 0.305 0.000 -0.000 -0.000 3/10125 350.0 0.297 0.297 0.297 0.000 -0.000 -0.000 3/10125 360.0 0.288 0.288 0.288 0.000 -0.000 -0.000 3/10125 370.0 0.281 0.281 0.281 0.000 -0.000 -0.000 3/10125 380.0 0.273 0.273 0.273 0.000 -0.000 -0.000 3/10125 390.0 0.266 0.266 0.266 0.000 -0.000 -0.000 3/10125 400.0 0.260 0.260 0.260 0.000 -0.000 -0.000 3/10125 410.0 0.253 0.253 0.253 0.000 -0.000 -0.000 3/10125 420.0 0.247 0.247 0.247 0.000 -0.000 -0.000 3/10125 430.0 0.242 0.242 0.242 0.000 -0.000 -0.000 3/10125 440.0 0.236 0.236 0.236 0.000 -0.000 -0.000 3/10125 450.0 0.231 0.231 0.231 0.000 -0.000 -0.000 3/10125 460.0 0.226 0.226 0.226 0.000 -0.000 -0.000 3/10125 470.0 0.221 0.221 0.221 0.000 -0.000 -0.000 3/10125 480.0 0.217 0.217 0.217 0.000 -0.000 -0.000 3/10125 490.0 0.212 0.212 0.212 0.000 -0.000 -0.000 3/10125 500.0 0.208 0.208 0.208 0.000 -0.000 -0.000 3/10125 510.0 0.204 0.204 0.204 0.000 -0.000 -0.000 3/10125 520.0 0.200 0.200 0.200 0.000 -0.000 -0.000 3/10125 530.0 0.196 0.196 0.196 0.000 -0.000 -0.000 3/10125 540.0 0.193 0.193 0.193 0.000 -0.000 -0.000 3/10125 550.0 0.189 0.189 0.189 0.000 -0.000 -0.000 3/10125 560.0 0.186 0.186 0.186 0.000 -0.000 -0.000 3/10125 570.0 0.183 0.183 0.183 0.000 -0.000 -0.000 3/10125 580.0 0.179 0.179 0.179 0.000 -0.000 -0.000 3/10125 590.0 0.176 0.176 0.176 0.000 -0.000 -0.000 3/10125 600.0 0.174 0.174 0.174 0.000 -0.000 -0.000 3/10125 610.0 0.171 0.171 0.171 0.000 -0.000 -0.000 3/10125 620.0 0.168 0.168 0.168 0.000 -0.000 -0.000 3/10125 630.0 0.165 0.165 0.165 0.000 -0.000 -0.000 3/10125 640.0 0.163 0.163 0.163 0.000 -0.000 -0.000 3/10125 650.0 0.160 0.160 0.160 0.000 -0.000 -0.000 3/10125 660.0 0.158 0.158 0.158 0.000 -0.000 -0.000 3/10125 670.0 0.155 0.155 0.155 0.000 -0.000 -0.000 3/10125 680.0 0.153 0.153 0.153 0.000 -0.000 -0.000 3/10125 690.0 0.151 0.151 0.151 0.000 -0.000 -0.000 3/10125 700.0 0.149 0.149 0.149 0.000 -0.000 -0.000 3/10125 710.0 0.147 0.147 0.147 0.000 -0.000 -0.000 3/10125 720.0 0.145 0.145 0.145 0.000 -0.000 -0.000 3/10125 730.0 0.143 0.143 0.143 0.000 -0.000 -0.000 3/10125 740.0 0.141 0.141 0.141 0.000 -0.000 -0.000 3/10125 750.0 0.139 0.139 0.139 0.000 -0.000 -0.000 3/10125 760.0 0.137 0.137 0.137 0.000 -0.000 -0.000 3/10125 770.0 0.135 0.135 0.135 0.000 -0.000 -0.000 3/10125 780.0 0.134 0.134 0.134 0.000 -0.000 -0.000 3/10125 790.0 0.132 0.132 0.132 0.000 -0.000 -0.000 3/10125 800.0 0.130 0.130 0.130 0.000 -0.000 -0.000 3/10125 810.0 0.129 0.129 0.129 0.000 -0.000 -0.000 3/10125 820.0 0.127 0.127 0.127 0.000 -0.000 -0.000 3/10125 830.0 0.126 0.126 0.126 0.000 -0.000 -0.000 3/10125 840.0 0.124 0.124 0.124 0.000 -0.000 -0.000 3/10125 850.0 0.123 0.123 0.123 0.000 -0.000 -0.000 3/10125 860.0 0.121 0.121 0.121 0.000 -0.000 -0.000 3/10125 870.0 0.120 0.120 0.120 0.000 -0.000 -0.000 3/10125 880.0 0.118 0.118 0.118 0.000 -0.000 -0.000 3/10125 890.0 0.117 0.117 0.117 0.000 -0.000 -0.000 3/10125 900.0 0.116 0.116 0.116 0.000 -0.000 -0.000 3/10125 910.0 0.115 0.115 0.115 0.000 -0.000 -0.000 3/10125 920.0 0.113 0.113 0.113 0.000 -0.000 -0.000 3/10125 930.0 0.112 0.112 0.112 0.000 -0.000 -0.000 3/10125 940.0 0.111 0.111 0.111 0.000 -0.000 -0.000 3/10125 950.0 0.110 0.110 0.110 0.000 -0.000 -0.000 3/10125 960.0 0.109 0.109 0.109 0.000 -0.000 -0.000 3/10125 970.0 0.107 0.107 0.107 0.000 -0.000 -0.000 3/10125 980.0 0.106 0.106 0.106 0.000 -0.000 -0.000 3/10125 990.0 0.105 0.105 0.105 0.000 -0.000 -0.000 3/10125 1000.0 0.104 0.104 0.104 0.000 -0.000 -0.000 3/10125 Thermal conductivity related properties were written into "kappa-m151515.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 05:20:35]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|