----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-09 04:49:49]------------------------- Compiled with OpenMP support (max 128 threads). Running in phonopy.load mode. Python version 3.14.2 Spglib version 2.6.1 Crystal structure was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". Unit of length: angstrom Settings: Supercell: [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 4.549011500000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.549011500000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.549011500000000 Atomic positions (fractional): *1 Na 0.00000000000000 0.50000000000000 0.50000000000000 22.990 2 Na 0.50000000000000 0.00000000000000 0.50000000000000 22.990 3 Na 0.50000000000000 0.50000000000000 0.00000000000000 22.990 *4 Br 0.00000000000000 0.00000000000000 0.00000000000000 79.904 *5 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.549011500000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.549011500000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.549011500000000 Atomic positions (fractional): *1 Na 0.00000000000000 0.50000000000000 0.50000000000000 22.990 > 1 2 Na 0.50000000000000 0.00000000000000 0.50000000000000 22.990 > 2 3 Na 0.50000000000000 0.50000000000000 0.00000000000000 22.990 > 3 *4 Br 0.00000000000000 0.00000000000000 0.00000000000000 79.904 > 4 *5 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 5 -------------------------------- super cell -------------------------------- Lattice vectors: a 13.647034500000000 0.000000000000000 0.000000000000000 b 0.000000000000000 13.647034500000000 0.000000000000000 c 0.000000000000000 0.000000000000000 13.647034500000000 Atomic positions (fractional): *1 Na 0.00000000000000 0.16666666666667 0.16666666666667 22.990 > 1 2 Na 0.33333333333333 0.16666666666667 0.16666666666667 22.990 > 1 3 Na 0.66666666666667 0.16666666666667 0.16666666666667 22.990 > 1 4 Na 0.00000000000000 0.50000000000000 0.16666666666667 22.990 > 1 5 Na 0.33333333333333 0.50000000000000 0.16666666666667 22.990 > 1 6 Na 0.66666666666667 0.50000000000000 0.16666666666667 22.990 > 1 7 Na 0.00000000000000 0.83333333333333 0.16666666666667 22.990 > 1 8 Na 0.33333333333333 0.83333333333333 0.16666666666667 22.990 > 1 9 Na 0.66666666666667 0.83333333333333 0.16666666666667 22.990 > 1 10 Na 0.00000000000000 0.16666666666667 0.50000000000000 22.990 > 1 11 Na 0.33333333333333 0.16666666666667 0.50000000000000 22.990 > 1 12 Na 0.66666666666667 0.16666666666667 0.50000000000000 22.990 > 1 13 Na 0.00000000000000 0.50000000000000 0.50000000000000 22.990 > 1 14 Na 0.33333333333333 0.50000000000000 0.50000000000000 22.990 > 1 15 Na 0.66666666666667 0.50000000000000 0.50000000000000 22.990 > 1 16 Na 0.00000000000000 0.83333333333333 0.50000000000000 22.990 > 1 17 Na 0.33333333333333 0.83333333333333 0.50000000000000 22.990 > 1 18 Na 0.66666666666667 0.83333333333333 0.50000000000000 22.990 > 1 19 Na 0.00000000000000 0.16666666666667 0.83333333333333 22.990 > 1 20 Na 0.33333333333333 0.16666666666667 0.83333333333333 22.990 > 1 21 Na 0.66666666666667 0.16666666666667 0.83333333333333 22.990 > 1 22 Na 0.00000000000000 0.50000000000000 0.83333333333333 22.990 > 1 23 Na 0.33333333333333 0.50000000000000 0.83333333333333 22.990 > 1 24 Na 0.66666666666667 0.50000000000000 0.83333333333333 22.990 > 1 25 Na 0.00000000000000 0.83333333333333 0.83333333333333 22.990 > 1 26 Na 0.33333333333333 0.83333333333333 0.83333333333333 22.990 > 1 27 Na 0.66666666666667 0.83333333333333 0.83333333333333 22.990 > 1 28 Na 0.16666666666667 0.00000000000000 0.16666666666667 22.990 > 2 29 Na 0.50000000000000 0.00000000000000 0.16666666666667 22.990 > 2 30 Na 0.83333333333333 0.00000000000000 0.16666666666667 22.990 > 2 31 Na 0.16666666666667 0.33333333333333 0.16666666666667 22.990 > 2 32 Na 0.50000000000000 0.33333333333333 0.16666666666667 22.990 > 2 33 Na 0.83333333333333 0.33333333333333 0.16666666666667 22.990 > 2 34 Na 0.16666666666667 0.66666666666667 0.16666666666667 22.990 > 2 35 Na 0.50000000000000 0.66666666666667 0.16666666666667 22.990 > 2 36 Na 0.83333333333333 0.66666666666667 0.16666666666667 22.990 > 2 37 Na 0.16666666666667 0.00000000000000 0.50000000000000 22.990 > 2 38 Na 0.50000000000000 0.00000000000000 0.50000000000000 22.990 > 2 39 Na 0.83333333333333 0.00000000000000 0.50000000000000 22.990 > 2 40 Na 0.16666666666667 0.33333333333333 0.50000000000000 22.990 > 2 41 Na 0.50000000000000 0.33333333333333 0.50000000000000 22.990 > 2 42 Na 0.83333333333333 0.33333333333333 0.50000000000000 22.990 > 2 43 Na 0.16666666666667 0.66666666666667 0.50000000000000 22.990 > 2 44 Na 0.50000000000000 0.66666666666667 0.50000000000000 22.990 > 2 45 Na 0.83333333333333 0.66666666666667 0.50000000000000 22.990 > 2 46 Na 0.16666666666667 0.00000000000000 0.83333333333333 22.990 > 2 47 Na 0.50000000000000 0.00000000000000 0.83333333333333 22.990 > 2 48 Na 0.83333333333333 0.00000000000000 0.83333333333333 22.990 > 2 49 Na 0.16666666666667 0.33333333333333 0.83333333333333 22.990 > 2 50 Na 0.50000000000000 0.33333333333333 0.83333333333333 22.990 > 2 51 Na 0.83333333333333 0.33333333333333 0.83333333333333 22.990 > 2 52 Na 0.16666666666667 0.66666666666667 0.83333333333333 22.990 > 2 53 Na 0.50000000000000 0.66666666666667 0.83333333333333 22.990 > 2 54 Na 0.83333333333333 0.66666666666667 0.83333333333333 22.990 > 2 55 Na 0.16666666666667 0.16666666666667 0.00000000000000 22.990 > 3 56 Na 0.50000000000000 0.16666666666667 0.00000000000000 22.990 > 3 57 Na 0.83333333333333 0.16666666666667 0.00000000000000 22.990 > 3 58 Na 0.16666666666667 0.50000000000000 0.00000000000000 22.990 > 3 59 Na 0.50000000000000 0.50000000000000 0.00000000000000 22.990 > 3 60 Na 0.83333333333333 0.50000000000000 0.00000000000000 22.990 > 3 61 Na 0.16666666666667 0.83333333333333 0.00000000000000 22.990 > 3 62 Na 0.50000000000000 0.83333333333333 0.00000000000000 22.990 > 3 63 Na 0.83333333333333 0.83333333333333 0.00000000000000 22.990 > 3 64 Na 0.16666666666667 0.16666666666667 0.33333333333333 22.990 > 3 65 Na 0.50000000000000 0.16666666666667 0.33333333333333 22.990 > 3 66 Na 0.83333333333333 0.16666666666667 0.33333333333333 22.990 > 3 67 Na 0.16666666666667 0.50000000000000 0.33333333333333 22.990 > 3 68 Na 0.50000000000000 0.50000000000000 0.33333333333333 22.990 > 3 69 Na 0.83333333333333 0.50000000000000 0.33333333333333 22.990 > 3 70 Na 0.16666666666667 0.83333333333333 0.33333333333333 22.990 > 3 71 Na 0.50000000000000 0.83333333333333 0.33333333333333 22.990 > 3 72 Na 0.83333333333333 0.83333333333333 0.33333333333333 22.990 > 3 73 Na 0.16666666666667 0.16666666666667 0.66666666666667 22.990 > 3 74 Na 0.50000000000000 0.16666666666667 0.66666666666667 22.990 > 3 75 Na 0.83333333333333 0.16666666666667 0.66666666666667 22.990 > 3 76 Na 0.16666666666667 0.50000000000000 0.66666666666667 22.990 > 3 77 Na 0.50000000000000 0.50000000000000 0.66666666666667 22.990 > 3 78 Na 0.83333333333333 0.50000000000000 0.66666666666667 22.990 > 3 79 Na 0.16666666666667 0.83333333333333 0.66666666666667 22.990 > 3 80 Na 0.50000000000000 0.83333333333333 0.66666666666667 22.990 > 3 81 Na 0.83333333333333 0.83333333333333 0.66666666666667 22.990 > 3 *82 Br 0.00000000000000 0.00000000000000 0.00000000000000 79.904 > 4 83 Br 0.33333333333333 0.00000000000000 0.00000000000000 79.904 > 4 84 Br 0.66666666666667 0.00000000000000 0.00000000000000 79.904 > 4 85 Br 0.00000000000000 0.33333333333333 0.00000000000000 79.904 > 4 86 Br 0.33333333333333 0.33333333333333 0.00000000000000 79.904 > 4 87 Br 0.66666666666667 0.33333333333333 0.00000000000000 79.904 > 4 88 Br 0.00000000000000 0.66666666666667 0.00000000000000 79.904 > 4 89 Br 0.33333333333333 0.66666666666667 0.00000000000000 79.904 > 4 90 Br 0.66666666666667 0.66666666666667 0.00000000000000 79.904 > 4 91 Br 0.00000000000000 0.00000000000000 0.33333333333333 79.904 > 4 92 Br 0.33333333333333 0.00000000000000 0.33333333333333 79.904 > 4 93 Br 0.66666666666667 0.00000000000000 0.33333333333333 79.904 > 4 94 Br 0.00000000000000 0.33333333333333 0.33333333333333 79.904 > 4 95 Br 0.33333333333333 0.33333333333333 0.33333333333333 79.904 > 4 96 Br 0.66666666666667 0.33333333333333 0.33333333333333 79.904 > 4 97 Br 0.00000000000000 0.66666666666667 0.33333333333333 79.904 > 4 98 Br 0.33333333333333 0.66666666666667 0.33333333333333 79.904 > 4 99 Br 0.66666666666667 0.66666666666667 0.33333333333333 79.904 > 4 100 Br 0.00000000000000 0.00000000000000 0.66666666666667 79.904 > 4 101 Br 0.33333333333333 0.00000000000000 0.66666666666667 79.904 > 4 102 Br 0.66666666666667 0.00000000000000 0.66666666666667 79.904 > 4 103 Br 0.00000000000000 0.33333333333333 0.66666666666667 79.904 > 4 104 Br 0.33333333333333 0.33333333333333 0.66666666666667 79.904 > 4 105 Br 0.66666666666667 0.33333333333333 0.66666666666667 79.904 > 4 106 Br 0.00000000000000 0.66666666666667 0.66666666666667 79.904 > 4 107 Br 0.33333333333333 0.66666666666667 0.66666666666667 79.904 > 4 108 Br 0.66666666666667 0.66666666666667 0.66666666666667 79.904 > 4 *109 O 0.16666666666667 0.16666666666667 0.16666666666667 15.999 > 5 110 O 0.50000000000000 0.16666666666667 0.16666666666667 15.999 > 5 111 O 0.83333333333333 0.16666666666667 0.16666666666667 15.999 > 5 112 O 0.16666666666667 0.50000000000000 0.16666666666667 15.999 > 5 113 O 0.50000000000000 0.50000000000000 0.16666666666667 15.999 > 5 114 O 0.83333333333333 0.50000000000000 0.16666666666667 15.999 > 5 115 O 0.16666666666667 0.83333333333333 0.16666666666667 15.999 > 5 116 O 0.50000000000000 0.83333333333333 0.16666666666667 15.999 > 5 117 O 0.83333333333333 0.83333333333333 0.16666666666667 15.999 > 5 118 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 5 119 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 5 120 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 5 121 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 5 122 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 5 123 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 5 124 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 5 125 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 5 126 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 5 127 O 0.16666666666667 0.16666666666667 0.83333333333333 15.999 > 5 128 O 0.50000000000000 0.16666666666667 0.83333333333333 15.999 > 5 129 O 0.83333333333333 0.16666666666667 0.83333333333333 15.999 > 5 130 O 0.16666666666667 0.50000000000000 0.83333333333333 15.999 > 5 131 O 0.50000000000000 0.50000000000000 0.83333333333333 15.999 > 5 132 O 0.83333333333333 0.50000000000000 0.83333333333333 15.999 > 5 133 O 0.16666666666667 0.83333333333333 0.83333333333333 15.999 > 5 134 O 0.50000000000000 0.83333333333333 0.83333333333333 15.999 > 5 135 O 0.83333333333333 0.83333333333333 0.83333333333333 15.999 > 5 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.3407171 0.0000000 0.0000000 0.0000000 3.3407171 0.0000000 0.0000000 0.0000000 3.3407171 -------------------------- Born effective charges -------------------------- 1 Na 0.8715923 0.0000000 0.0000000 0.0000000 1.1178688 0.0000000 0.0000000 0.0000000 1.1178688 2 Na 1.1178688 0.0000000 0.0000000 0.0000000 0.8715923 0.0000000 0.0000000 0.0000000 1.1178688 3 Na 1.1178688 0.0000000 0.0000000 0.0000000 1.1178688 0.0000000 0.0000000 0.0000000 0.8715923 4 Br -1.2258001 0.0000000 0.0000000 0.0000000 -1.2258001 0.0000000 0.0000000 0.0000000 -1.2258001 5 O -1.8815299 0.0000000 0.0000000 0.0000000 -1.8815299 0.0000000 0.0000000 0.0000000 -1.8815299 ---------------------------------------------------------------------------- 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.004 Solver_block: 80 / 80 - Time: 0.060 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.064 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (yy) -0.00000000 (yy) 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 (yy) -0.00000000 (yy) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-09 04:49:53]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:49: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: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.549011500000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.549011500000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.549011500000000 Atomic positions (fractional): 1 Na 0.00000000000000 0.50000000000000 0.50000000000000 22.990 2 Na 0.50000000000000 0.00000000000000 0.50000000000000 22.990 3 Na 0.50000000000000 0.50000000000000 0.00000000000000 22.990 4 Br 0.00000000000000 0.00000000000000 0.00000000000000 79.904 5 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 -------------------------------- supercell --------------------------------- Lattice vectors: a 13.647034500000000 0.000000000000000 0.000000000000000 b 0.000000000000000 13.647034500000000 0.000000000000000 c 0.000000000000000 0.000000000000000 13.647034500000000 Atomic positions (fractional): 1 Na 0.00000000000000 0.16666666666667 0.16666666666667 22.990 > 1 2 Na 0.33333333333333 0.16666666666667 0.16666666666667 22.990 > 1 3 Na 0.66666666666667 0.16666666666667 0.16666666666667 22.990 > 1 4 Na 0.00000000000000 0.50000000000000 0.16666666666667 22.990 > 1 5 Na 0.33333333333333 0.50000000000000 0.16666666666667 22.990 > 1 6 Na 0.66666666666667 0.50000000000000 0.16666666666667 22.990 > 1 7 Na 0.00000000000000 0.83333333333333 0.16666666666667 22.990 > 1 8 Na 0.33333333333333 0.83333333333333 0.16666666666667 22.990 > 1 9 Na 0.66666666666667 0.83333333333333 0.16666666666667 22.990 > 1 10 Na 0.00000000000000 0.16666666666667 0.50000000000000 22.990 > 1 11 Na 0.33333333333333 0.16666666666667 0.50000000000000 22.990 > 1 12 Na 0.66666666666667 0.16666666666667 0.50000000000000 22.990 > 1 13 Na 0.00000000000000 0.50000000000000 0.50000000000000 22.990 > 1 14 Na 0.33333333333333 0.50000000000000 0.50000000000000 22.990 > 1 15 Na 0.66666666666667 0.50000000000000 0.50000000000000 22.990 > 1 16 Na 0.00000000000000 0.83333333333333 0.50000000000000 22.990 > 1 17 Na 0.33333333333333 0.83333333333333 0.50000000000000 22.990 > 1 18 Na 0.66666666666667 0.83333333333333 0.50000000000000 22.990 > 1 19 Na 0.00000000000000 0.16666666666667 0.83333333333333 22.990 > 1 20 Na 0.33333333333333 0.16666666666667 0.83333333333333 22.990 > 1 21 Na 0.66666666666667 0.16666666666667 0.83333333333333 22.990 > 1 22 Na 0.00000000000000 0.50000000000000 0.83333333333333 22.990 > 1 23 Na 0.33333333333333 0.50000000000000 0.83333333333333 22.990 > 1 24 Na 0.66666666666667 0.50000000000000 0.83333333333333 22.990 > 1 25 Na 0.00000000000000 0.83333333333333 0.83333333333333 22.990 > 1 26 Na 0.33333333333333 0.83333333333333 0.83333333333333 22.990 > 1 27 Na 0.66666666666667 0.83333333333333 0.83333333333333 22.990 > 1 28 Na 0.16666666666667 0.00000000000000 0.16666666666667 22.990 > 28 29 Na 0.50000000000000 0.00000000000000 0.16666666666667 22.990 > 28 30 Na 0.83333333333333 0.00000000000000 0.16666666666667 22.990 > 28 31 Na 0.16666666666667 0.33333333333333 0.16666666666667 22.990 > 28 32 Na 0.50000000000000 0.33333333333333 0.16666666666667 22.990 > 28 33 Na 0.83333333333333 0.33333333333333 0.16666666666667 22.990 > 28 34 Na 0.16666666666667 0.66666666666667 0.16666666666667 22.990 > 28 35 Na 0.50000000000000 0.66666666666667 0.16666666666667 22.990 > 28 36 Na 0.83333333333333 0.66666666666667 0.16666666666667 22.990 > 28 37 Na 0.16666666666667 0.00000000000000 0.50000000000000 22.990 > 28 38 Na 0.50000000000000 0.00000000000000 0.50000000000000 22.990 > 28 39 Na 0.83333333333333 0.00000000000000 0.50000000000000 22.990 > 28 40 Na 0.16666666666667 0.33333333333333 0.50000000000000 22.990 > 28 41 Na 0.50000000000000 0.33333333333333 0.50000000000000 22.990 > 28 42 Na 0.83333333333333 0.33333333333333 0.50000000000000 22.990 > 28 43 Na 0.16666666666667 0.66666666666667 0.50000000000000 22.990 > 28 44 Na 0.50000000000000 0.66666666666667 0.50000000000000 22.990 > 28 45 Na 0.83333333333333 0.66666666666667 0.50000000000000 22.990 > 28 46 Na 0.16666666666667 0.00000000000000 0.83333333333333 22.990 > 28 47 Na 0.50000000000000 0.00000000000000 0.83333333333333 22.990 > 28 48 Na 0.83333333333333 0.00000000000000 0.83333333333333 22.990 > 28 49 Na 0.16666666666667 0.33333333333333 0.83333333333333 22.990 > 28 50 Na 0.50000000000000 0.33333333333333 0.83333333333333 22.990 > 28 51 Na 0.83333333333333 0.33333333333333 0.83333333333333 22.990 > 28 52 Na 0.16666666666667 0.66666666666667 0.83333333333333 22.990 > 28 53 Na 0.50000000000000 0.66666666666667 0.83333333333333 22.990 > 28 54 Na 0.83333333333333 0.66666666666667 0.83333333333333 22.990 > 28 55 Na 0.16666666666667 0.16666666666667 0.00000000000000 22.990 > 55 56 Na 0.50000000000000 0.16666666666667 0.00000000000000 22.990 > 55 57 Na 0.83333333333333 0.16666666666667 0.00000000000000 22.990 > 55 58 Na 0.16666666666667 0.50000000000000 0.00000000000000 22.990 > 55 59 Na 0.50000000000000 0.50000000000000 0.00000000000000 22.990 > 55 60 Na 0.83333333333333 0.50000000000000 0.00000000000000 22.990 > 55 61 Na 0.16666666666667 0.83333333333333 0.00000000000000 22.990 > 55 62 Na 0.50000000000000 0.83333333333333 0.00000000000000 22.990 > 55 63 Na 0.83333333333333 0.83333333333333 0.00000000000000 22.990 > 55 64 Na 0.16666666666667 0.16666666666667 0.33333333333333 22.990 > 55 65 Na 0.50000000000000 0.16666666666667 0.33333333333333 22.990 > 55 66 Na 0.83333333333333 0.16666666666667 0.33333333333333 22.990 > 55 67 Na 0.16666666666667 0.50000000000000 0.33333333333333 22.990 > 55 68 Na 0.50000000000000 0.50000000000000 0.33333333333333 22.990 > 55 69 Na 0.83333333333333 0.50000000000000 0.33333333333333 22.990 > 55 70 Na 0.16666666666667 0.83333333333333 0.33333333333333 22.990 > 55 71 Na 0.50000000000000 0.83333333333333 0.33333333333333 22.990 > 55 72 Na 0.83333333333333 0.83333333333333 0.33333333333333 22.990 > 55 73 Na 0.16666666666667 0.16666666666667 0.66666666666667 22.990 > 55 74 Na 0.50000000000000 0.16666666666667 0.66666666666667 22.990 > 55 75 Na 0.83333333333333 0.16666666666667 0.66666666666667 22.990 > 55 76 Na 0.16666666666667 0.50000000000000 0.66666666666667 22.990 > 55 77 Na 0.50000000000000 0.50000000000000 0.66666666666667 22.990 > 55 78 Na 0.83333333333333 0.50000000000000 0.66666666666667 22.990 > 55 79 Na 0.16666666666667 0.83333333333333 0.66666666666667 22.990 > 55 80 Na 0.50000000000000 0.83333333333333 0.66666666666667 22.990 > 55 81 Na 0.83333333333333 0.83333333333333 0.66666666666667 22.990 > 55 82 Br 0.00000000000000 0.00000000000000 0.00000000000000 79.904 > 82 83 Br 0.33333333333333 0.00000000000000 0.00000000000000 79.904 > 82 84 Br 0.66666666666667 0.00000000000000 0.00000000000000 79.904 > 82 85 Br 0.00000000000000 0.33333333333333 0.00000000000000 79.904 > 82 86 Br 0.33333333333333 0.33333333333333 0.00000000000000 79.904 > 82 87 Br 0.66666666666667 0.33333333333333 0.00000000000000 79.904 > 82 88 Br 0.00000000000000 0.66666666666667 0.00000000000000 79.904 > 82 89 Br 0.33333333333333 0.66666666666667 0.00000000000000 79.904 > 82 90 Br 0.66666666666667 0.66666666666667 0.00000000000000 79.904 > 82 91 Br 0.00000000000000 0.00000000000000 0.33333333333333 79.904 > 82 92 Br 0.33333333333333 0.00000000000000 0.33333333333333 79.904 > 82 93 Br 0.66666666666667 0.00000000000000 0.33333333333333 79.904 > 82 94 Br 0.00000000000000 0.33333333333333 0.33333333333333 79.904 > 82 95 Br 0.33333333333333 0.33333333333333 0.33333333333333 79.904 > 82 96 Br 0.66666666666667 0.33333333333333 0.33333333333333 79.904 > 82 97 Br 0.00000000000000 0.66666666666667 0.33333333333333 79.904 > 82 98 Br 0.33333333333333 0.66666666666667 0.33333333333333 79.904 > 82 99 Br 0.66666666666667 0.66666666666667 0.33333333333333 79.904 > 82 100 Br 0.00000000000000 0.00000000000000 0.66666666666667 79.904 > 82 101 Br 0.33333333333333 0.00000000000000 0.66666666666667 79.904 > 82 102 Br 0.66666666666667 0.00000000000000 0.66666666666667 79.904 > 82 103 Br 0.00000000000000 0.33333333333333 0.66666666666667 79.904 > 82 104 Br 0.33333333333333 0.33333333333333 0.66666666666667 79.904 > 82 105 Br 0.66666666666667 0.33333333333333 0.66666666666667 79.904 > 82 106 Br 0.00000000000000 0.66666666666667 0.66666666666667 79.904 > 82 107 Br 0.33333333333333 0.66666666666667 0.66666666666667 79.904 > 82 108 Br 0.66666666666667 0.66666666666667 0.66666666666667 79.904 > 82 109 O 0.16666666666667 0.16666666666667 0.16666666666667 15.999 > 109 110 O 0.50000000000000 0.16666666666667 0.16666666666667 15.999 > 109 111 O 0.83333333333333 0.16666666666667 0.16666666666667 15.999 > 109 112 O 0.16666666666667 0.50000000000000 0.16666666666667 15.999 > 109 113 O 0.50000000000000 0.50000000000000 0.16666666666667 15.999 > 109 114 O 0.83333333333333 0.50000000000000 0.16666666666667 15.999 > 109 115 O 0.16666666666667 0.83333333333333 0.16666666666667 15.999 > 109 116 O 0.50000000000000 0.83333333333333 0.16666666666667 15.999 > 109 117 O 0.83333333333333 0.83333333333333 0.16666666666667 15.999 > 109 118 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 109 119 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 109 120 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 109 121 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 109 122 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 109 123 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 109 124 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 109 125 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 109 126 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 109 127 O 0.16666666666667 0.16666666666667 0.83333333333333 15.999 > 109 128 O 0.50000000000000 0.16666666666667 0.83333333333333 15.999 > 109 129 O 0.83333333333333 0.16666666666667 0.83333333333333 15.999 > 109 130 O 0.16666666666667 0.50000000000000 0.83333333333333 15.999 > 109 131 O 0.50000000000000 0.50000000000000 0.83333333333333 15.999 > 109 132 O 0.83333333333333 0.50000000000000 0.83333333333333 15.999 > 109 133 O 0.16666666666667 0.83333333333333 0.83333333333333 15.999 > 109 134 O 0.50000000000000 0.83333333333333 0.83333333333333 15.999 > 109 135 O 0.83333333333333 0.83333333333333 0.83333333333333 15.999 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 3.3407171 0.0000000 0.0000000 0.0000000 3.3407171 0.0000000 0.0000000 0.0000000 3.3407171 -------------------------- Born effective charges -------------------------- 1 Na 0.8715923 0.0000000 0.0000000 0.0000000 1.1178688 0.0000000 0.0000000 0.0000000 1.1178688 2 Na 1.1178688 0.0000000 0.0000000 0.0000000 0.8715923 0.0000000 0.0000000 0.0000000 1.1178688 3 Na 1.1178688 0.0000000 0.0000000 0.0000000 1.1178688 0.0000000 0.0000000 0.0000000 0.8715923 4 Br -1.2258001 0.0000000 0.0000000 0.0000000 -1.2258001 0.0000000 0.0000000 0.0000000 -1.2258001 5 O -1.8815299 0.0000000 0.0000000 0.0000000 -1.8815299 0.0000000 0.0000000 0.0000000 -1.8815299 ---------------------------------------------------------------------------- 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.0000 0.0100 0.0000] [ 0.0000 -0.0100 0.0000] [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 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-09 04:49:58]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:49:58]------------------------- 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: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: Pm-3m (221) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.549011500000000 0.000000000000000 0.000000000000000 b 0.000000000000000 4.549011500000000 0.000000000000000 c 0.000000000000000 0.000000000000000 4.549011500000000 Atomic positions (fractional): 1 Na 0.00000000000000 0.50000000000000 0.50000000000000 22.990 2 Na 0.50000000000000 0.00000000000000 0.50000000000000 22.990 3 Na 0.50000000000000 0.50000000000000 0.00000000000000 22.990 4 Br 0.00000000000000 0.00000000000000 0.00000000000000 79.904 5 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 -------------------------------- supercell --------------------------------- Lattice vectors: a 13.647034500000000 0.000000000000000 0.000000000000000 b 0.000000000000000 13.647034500000000 0.000000000000000 c 0.000000000000000 0.000000000000000 13.647034500000000 Atomic positions (fractional): 1 Na 0.00000000000000 0.16666666666667 0.16666666666667 22.990 > 1 2 Na 0.33333333333333 0.16666666666667 0.16666666666667 22.990 > 1 3 Na 0.66666666666667 0.16666666666667 0.16666666666667 22.990 > 1 4 Na 0.00000000000000 0.50000000000000 0.16666666666667 22.990 > 1 5 Na 0.33333333333333 0.50000000000000 0.16666666666667 22.990 > 1 6 Na 0.66666666666667 0.50000000000000 0.16666666666667 22.990 > 1 7 Na 0.00000000000000 0.83333333333333 0.16666666666667 22.990 > 1 8 Na 0.33333333333333 0.83333333333333 0.16666666666667 22.990 > 1 9 Na 0.66666666666667 0.83333333333333 0.16666666666667 22.990 > 1 10 Na 0.00000000000000 0.16666666666667 0.50000000000000 22.990 > 1 11 Na 0.33333333333333 0.16666666666667 0.50000000000000 22.990 > 1 12 Na 0.66666666666667 0.16666666666667 0.50000000000000 22.990 > 1 13 Na 0.00000000000000 0.50000000000000 0.50000000000000 22.990 > 1 14 Na 0.33333333333333 0.50000000000000 0.50000000000000 22.990 > 1 15 Na 0.66666666666667 0.50000000000000 0.50000000000000 22.990 > 1 16 Na 0.00000000000000 0.83333333333333 0.50000000000000 22.990 > 1 17 Na 0.33333333333333 0.83333333333333 0.50000000000000 22.990 > 1 18 Na 0.66666666666667 0.83333333333333 0.50000000000000 22.990 > 1 19 Na 0.00000000000000 0.16666666666667 0.83333333333333 22.990 > 1 20 Na 0.33333333333333 0.16666666666667 0.83333333333333 22.990 > 1 21 Na 0.66666666666667 0.16666666666667 0.83333333333333 22.990 > 1 22 Na 0.00000000000000 0.50000000000000 0.83333333333333 22.990 > 1 23 Na 0.33333333333333 0.50000000000000 0.83333333333333 22.990 > 1 24 Na 0.66666666666667 0.50000000000000 0.83333333333333 22.990 > 1 25 Na 0.00000000000000 0.83333333333333 0.83333333333333 22.990 > 1 26 Na 0.33333333333333 0.83333333333333 0.83333333333333 22.990 > 1 27 Na 0.66666666666667 0.83333333333333 0.83333333333333 22.990 > 1 28 Na 0.16666666666667 0.00000000000000 0.16666666666667 22.990 > 28 29 Na 0.50000000000000 0.00000000000000 0.16666666666667 22.990 > 28 30 Na 0.83333333333333 0.00000000000000 0.16666666666667 22.990 > 28 31 Na 0.16666666666667 0.33333333333333 0.16666666666667 22.990 > 28 32 Na 0.50000000000000 0.33333333333333 0.16666666666667 22.990 > 28 33 Na 0.83333333333333 0.33333333333333 0.16666666666667 22.990 > 28 34 Na 0.16666666666667 0.66666666666667 0.16666666666667 22.990 > 28 35 Na 0.50000000000000 0.66666666666667 0.16666666666667 22.990 > 28 36 Na 0.83333333333333 0.66666666666667 0.16666666666667 22.990 > 28 37 Na 0.16666666666667 0.00000000000000 0.50000000000000 22.990 > 28 38 Na 0.50000000000000 0.00000000000000 0.50000000000000 22.990 > 28 39 Na 0.83333333333333 0.00000000000000 0.50000000000000 22.990 > 28 40 Na 0.16666666666667 0.33333333333333 0.50000000000000 22.990 > 28 41 Na 0.50000000000000 0.33333333333333 0.50000000000000 22.990 > 28 42 Na 0.83333333333333 0.33333333333333 0.50000000000000 22.990 > 28 43 Na 0.16666666666667 0.66666666666667 0.50000000000000 22.990 > 28 44 Na 0.50000000000000 0.66666666666667 0.50000000000000 22.990 > 28 45 Na 0.83333333333333 0.66666666666667 0.50000000000000 22.990 > 28 46 Na 0.16666666666667 0.00000000000000 0.83333333333333 22.990 > 28 47 Na 0.50000000000000 0.00000000000000 0.83333333333333 22.990 > 28 48 Na 0.83333333333333 0.00000000000000 0.83333333333333 22.990 > 28 49 Na 0.16666666666667 0.33333333333333 0.83333333333333 22.990 > 28 50 Na 0.50000000000000 0.33333333333333 0.83333333333333 22.990 > 28 51 Na 0.83333333333333 0.33333333333333 0.83333333333333 22.990 > 28 52 Na 0.16666666666667 0.66666666666667 0.83333333333333 22.990 > 28 53 Na 0.50000000000000 0.66666666666667 0.83333333333333 22.990 > 28 54 Na 0.83333333333333 0.66666666666667 0.83333333333333 22.990 > 28 55 Na 0.16666666666667 0.16666666666667 0.00000000000000 22.990 > 55 56 Na 0.50000000000000 0.16666666666667 0.00000000000000 22.990 > 55 57 Na 0.83333333333333 0.16666666666667 0.00000000000000 22.990 > 55 58 Na 0.16666666666667 0.50000000000000 0.00000000000000 22.990 > 55 59 Na 0.50000000000000 0.50000000000000 0.00000000000000 22.990 > 55 60 Na 0.83333333333333 0.50000000000000 0.00000000000000 22.990 > 55 61 Na 0.16666666666667 0.83333333333333 0.00000000000000 22.990 > 55 62 Na 0.50000000000000 0.83333333333333 0.00000000000000 22.990 > 55 63 Na 0.83333333333333 0.83333333333333 0.00000000000000 22.990 > 55 64 Na 0.16666666666667 0.16666666666667 0.33333333333333 22.990 > 55 65 Na 0.50000000000000 0.16666666666667 0.33333333333333 22.990 > 55 66 Na 0.83333333333333 0.16666666666667 0.33333333333333 22.990 > 55 67 Na 0.16666666666667 0.50000000000000 0.33333333333333 22.990 > 55 68 Na 0.50000000000000 0.50000000000000 0.33333333333333 22.990 > 55 69 Na 0.83333333333333 0.50000000000000 0.33333333333333 22.990 > 55 70 Na 0.16666666666667 0.83333333333333 0.33333333333333 22.990 > 55 71 Na 0.50000000000000 0.83333333333333 0.33333333333333 22.990 > 55 72 Na 0.83333333333333 0.83333333333333 0.33333333333333 22.990 > 55 73 Na 0.16666666666667 0.16666666666667 0.66666666666667 22.990 > 55 74 Na 0.50000000000000 0.16666666666667 0.66666666666667 22.990 > 55 75 Na 0.83333333333333 0.16666666666667 0.66666666666667 22.990 > 55 76 Na 0.16666666666667 0.50000000000000 0.66666666666667 22.990 > 55 77 Na 0.50000000000000 0.50000000000000 0.66666666666667 22.990 > 55 78 Na 0.83333333333333 0.50000000000000 0.66666666666667 22.990 > 55 79 Na 0.16666666666667 0.83333333333333 0.66666666666667 22.990 > 55 80 Na 0.50000000000000 0.83333333333333 0.66666666666667 22.990 > 55 81 Na 0.83333333333333 0.83333333333333 0.66666666666667 22.990 > 55 82 Br 0.00000000000000 0.00000000000000 0.00000000000000 79.904 > 82 83 Br 0.33333333333333 0.00000000000000 0.00000000000000 79.904 > 82 84 Br 0.66666666666667 0.00000000000000 0.00000000000000 79.904 > 82 85 Br 0.00000000000000 0.33333333333333 0.00000000000000 79.904 > 82 86 Br 0.33333333333333 0.33333333333333 0.00000000000000 79.904 > 82 87 Br 0.66666666666667 0.33333333333333 0.00000000000000 79.904 > 82 88 Br 0.00000000000000 0.66666666666667 0.00000000000000 79.904 > 82 89 Br 0.33333333333333 0.66666666666667 0.00000000000000 79.904 > 82 90 Br 0.66666666666667 0.66666666666667 0.00000000000000 79.904 > 82 91 Br 0.00000000000000 0.00000000000000 0.33333333333333 79.904 > 82 92 Br 0.33333333333333 0.00000000000000 0.33333333333333 79.904 > 82 93 Br 0.66666666666667 0.00000000000000 0.33333333333333 79.904 > 82 94 Br 0.00000000000000 0.33333333333333 0.33333333333333 79.904 > 82 95 Br 0.33333333333333 0.33333333333333 0.33333333333333 79.904 > 82 96 Br 0.66666666666667 0.33333333333333 0.33333333333333 79.904 > 82 97 Br 0.00000000000000 0.66666666666667 0.33333333333333 79.904 > 82 98 Br 0.33333333333333 0.66666666666667 0.33333333333333 79.904 > 82 99 Br 0.66666666666667 0.66666666666667 0.33333333333333 79.904 > 82 100 Br 0.00000000000000 0.00000000000000 0.66666666666667 79.904 > 82 101 Br 0.33333333333333 0.00000000000000 0.66666666666667 79.904 > 82 102 Br 0.66666666666667 0.00000000000000 0.66666666666667 79.904 > 82 103 Br 0.00000000000000 0.33333333333333 0.66666666666667 79.904 > 82 104 Br 0.33333333333333 0.33333333333333 0.66666666666667 79.904 > 82 105 Br 0.66666666666667 0.33333333333333 0.66666666666667 79.904 > 82 106 Br 0.00000000000000 0.66666666666667 0.66666666666667 79.904 > 82 107 Br 0.33333333333333 0.66666666666667 0.66666666666667 79.904 > 82 108 Br 0.66666666666667 0.66666666666667 0.66666666666667 79.904 > 82 109 O 0.16666666666667 0.16666666666667 0.16666666666667 15.999 > 109 110 O 0.50000000000000 0.16666666666667 0.16666666666667 15.999 > 109 111 O 0.83333333333333 0.16666666666667 0.16666666666667 15.999 > 109 112 O 0.16666666666667 0.50000000000000 0.16666666666667 15.999 > 109 113 O 0.50000000000000 0.50000000000000 0.16666666666667 15.999 > 109 114 O 0.83333333333333 0.50000000000000 0.16666666666667 15.999 > 109 115 O 0.16666666666667 0.83333333333333 0.16666666666667 15.999 > 109 116 O 0.50000000000000 0.83333333333333 0.16666666666667 15.999 > 109 117 O 0.83333333333333 0.83333333333333 0.16666666666667 15.999 > 109 118 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 109 119 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 109 120 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 109 121 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 109 122 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 109 123 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 109 124 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 109 125 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 109 126 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 109 127 O 0.16666666666667 0.16666666666667 0.83333333333333 15.999 > 109 128 O 0.50000000000000 0.16666666666667 0.83333333333333 15.999 > 109 129 O 0.83333333333333 0.16666666666667 0.83333333333333 15.999 > 109 130 O 0.16666666666667 0.50000000000000 0.83333333333333 15.999 > 109 131 O 0.50000000000000 0.50000000000000 0.83333333333333 15.999 > 109 132 O 0.83333333333333 0.50000000000000 0.83333333333333 15.999 > 109 133 O 0.16666666666667 0.83333333333333 0.83333333333333 15.999 > 109 134 O 0.50000000000000 0.83333333333333 0.83333333333333 15.999 > 109 135 O 0.83333333333333 0.83333333333333 0.83333333333333 15.999 > 109 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 3.3407171 0.0000000 0.0000000 0.0000000 3.3407171 0.0000000 0.0000000 0.0000000 3.3407171 -------------------------- Born effective charges -------------------------- 1 Na 0.8715923 0.0000000 0.0000000 0.0000000 1.1178688 0.0000000 0.0000000 0.0000000 1.1178688 2 Na 1.1178688 0.0000000 0.0000000 0.0000000 0.8715923 0.0000000 0.0000000 0.0000000 1.1178688 3 Na 1.1178688 0.0000000 0.0000000 0.0000000 1.1178688 0.0000000 0.0000000 0.0000000 0.8715923 4 Br -1.2258001 0.0000000 0.0000000 0.0000000 -1.2258001 0.0000000 0.0000000 0.0000000 -1.2258001 5 O -1.8815299 0.0000000 0.0000000 0.0000000 -1.8815299 0.0000000 0.0000000 0.0000000 -1.8815299 ---------------------------------------------------------------------------- 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 (yy) -0.00000000 (yy) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 11 11 11 ] fc3-r2q-transformation over three atoms: True --------------------------- Phonon calculations ---------------------------- Use NAC by Gonze et al. (no real space sum in current implementation) PRB 50, 13035(R) (1994), PRB 55, 10355 (1997) G-cutoff distance: 0.91, Number of G-points: 305, Lambda: 0.17 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) -------------------- ======================= Grid point 0 (1/56) ======================= q-point: ( 0.00 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 56 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 2.612 ( 0.000 0.000 0.000) 0.000 2.612 ( 0.000 0.000 0.000) 0.000 2.612 ( 0.000 0.000 0.000) 0.000 3.346 ( 0.000 0.000 0.000) 0.000 3.346 ( 0.000 0.000 0.000) 0.000 3.346 ( 0.000 0.000 0.000) 0.000 6.280 ( 0.000 0.000 0.000) 0.000 6.280 ( 0.000 0.000 0.000) 0.000 6.280 ( 0.000 0.000 0.000) 0.000 11.210 ( 0.000 0.000 0.000) 0.000 11.210 ( 0.000 0.000 0.000) 0.000 11.210 ( 0.000 0.000 0.000) 0.000 ======================= Grid point 1 (2/56) ======================= q-point: ( 0.09 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.536 ( 26.144 0.000 0.000) 26.144 0.536 ( 26.144 0.000 0.000) 26.144 0.941 ( 43.514 0.000 0.000) 43.514 2.611 ( -0.089 0.000 0.000) 0.089 2.611 ( -0.089 0.000 0.000) 0.089 3.368 ( 2.228 0.000 0.000) 2.228 3.368 ( 2.228 0.000 0.000) 2.228 3.373 ( 2.671 0.000 0.000) 2.671 3.566 ( 12.443 0.000 0.000) 12.443 6.256 ( -2.394 0.000 0.000) 2.394 6.256 ( -2.394 0.000 0.000) 2.394 7.148 ( 5.481 0.000 0.000) 5.481 11.211 ( 0.078 0.000 0.000) 0.078 11.211 ( 0.078 0.000 0.000) 0.078 12.726 ( -5.928 0.000 0.000) 5.928 ======================= Grid point 2 (3/56) ======================= q-point: ( 0.18 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.032 ( 23.186 0.000 0.000) 23.186 1.032 ( 23.186 0.000 0.000) 23.186 1.687 ( 30.515 0.000 0.000) 30.515 2.605 ( -0.673 0.000 0.000) 0.673 2.605 ( -0.673 0.000 0.000) 0.673 3.435 ( 4.448 0.000 0.000) 4.448 3.435 ( 4.448 0.000 0.000) 4.448 3.446 ( 4.390 0.000 0.000) 4.390 3.920 ( 21.932 0.000 0.000) 21.932 6.187 ( -4.335 0.000 0.000) 4.335 6.187 ( -4.335 0.000 0.000) 4.335 7.320 ( 11.957 0.000 0.000) 11.957 11.213 ( 0.134 0.000 0.000) 0.134 11.213 ( 0.134 0.000 0.000) 0.134 12.552 ( -11.353 0.000 0.000) 11.353 ======================= Grid point 3 (4/56) ======================= q-point: ( 0.27 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.452 ( 18.654 0.000 0.000) 18.654 1.452 ( 18.654 0.000 0.000) 18.654 2.168 ( 18.289 0.000 0.000) 18.289 2.581 ( -1.752 0.000 0.000) 1.752 2.581 ( -1.752 0.000 0.000) 1.752 3.539 ( 4.640 0.000 0.000) 4.640 3.541 ( 5.964 0.000 0.000) 5.964 3.541 ( 5.964 0.000 0.000) 5.964 4.376 ( 21.855 0.000 0.000) 21.855 6.089 ( -5.269 0.000 0.000) 5.269 6.089 ( -5.269 0.000 0.000) 5.269 7.627 ( 18.431 0.000 0.000) 18.431 11.216 ( 0.149 0.000 0.000) 0.149 11.216 ( 0.149 0.000 0.000) 0.149 12.282 ( -15.161 0.000 0.000) 15.161 ======================= Grid point 4 (5/56) ======================= q-point: ( 0.36 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.770 ( 12.910 0.000 0.000) 12.910 1.770 ( 12.910 0.000 0.000) 12.910 2.444 ( 9.867 0.000 0.000) 9.867 2.536 ( -2.677 0.000 0.000) 2.677 2.536 ( -2.677 0.000 0.000) 2.677 3.622 ( 3.453 0.000 0.000) 3.453 3.660 ( 5.476 0.000 0.000) 5.476 3.660 ( 5.476 0.000 0.000) 5.476 4.740 ( 13.824 0.000 0.000) 13.824 5.988 ( -4.545 0.000 0.000) 4.545 5.988 ( -4.545 0.000 0.000) 4.545 8.021 ( 19.493 0.000 0.000) 19.493 11.219 ( 0.116 0.000 0.000) 0.116 11.219 ( 0.116 0.000 0.000) 0.116 11.973 ( -14.836 0.000 0.000) 14.836 ======================= Grid point 5 (6/56) ======================= q-point: ( 0.45 0.00 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 126 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.955 ( 5.111 0.000 0.000) 5.111 1.955 ( 5.111 0.000 0.000) 5.111 2.486 ( -1.720 0.000 0.000) 1.720 2.486 ( -1.720 0.000 0.000) 1.720 2.573 ( 3.176 0.000 0.000) 3.176 3.670 ( 1.269 0.000 0.000) 1.269 3.742 ( 2.280 0.000 0.000) 2.280 3.742 ( 2.280 0.000 0.000) 2.280 4.921 ( 4.454 0.000 0.000) 4.454 5.921 ( -1.845 0.000 0.000) 1.845 5.921 ( -1.845 0.000 0.000) 1.845 8.326 ( 8.959 0.000 0.000) 8.959 11.221 ( 0.044 0.000 0.000) 0.044 11.221 ( 0.044 0.000 0.000) 0.044 11.742 ( -6.746 0.000 0.000) 6.746 ======================= Grid point 12 (7/56) ======================= q-point: ( 0.09 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.753 ( 18.071 18.071 0.000) 25.557 0.808 ( 18.573 18.573 0.000) 26.267 1.259 ( 27.448 27.448 0.000) 38.817 2.614 ( 0.253 0.253 0.000) 0.358 2.675 ( 2.924 2.924 0.000) 4.134 3.308 ( -1.817 -1.817 0.000) 2.569 3.388 ( 2.060 2.060 0.000) 2.913 3.440 ( 3.292 3.292 0.000) 4.656 3.607 ( 10.078 10.078 0.000) 14.253 6.230 ( -2.547 -2.547 0.000) 3.602 6.323 ( 2.251 2.251 0.000) 3.183 7.196 ( 4.976 4.976 0.000) 7.037 11.148 ( -3.078 -3.078 0.000) 4.352 11.212 ( 0.056 0.056 0.000) 0.079 12.698 ( -4.486 -4.486 0.000) 6.345 ======================= Grid point 13 (8/56) ======================= q-point: ( 0.18 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.148 ( 20.030 10.672 0.000) 22.696 1.184 ( 18.719 12.079 0.000) 22.278 1.807 ( 24.193 11.130 0.000) 26.631 2.618 ( 0.044 1.337 0.000) 1.338 2.695 ( -0.217 6.902 0.000) 6.905 3.276 ( -1.365 -5.491 0.000) 5.658 3.451 ( 4.179 1.605 0.000) 4.477 3.574 ( 8.757 1.660 0.000) 8.913 3.948 ( 22.691 7.885 0.000) 24.022 6.157 ( -4.622 -2.989 0.000) 5.504 6.313 ( -2.603 9.529 0.000) 9.878 7.355 ( 11.201 3.669 0.000) 11.787 11.114 ( -0.741 -7.882 0.000) 7.917 11.213 ( 0.096 -0.003 0.000) 0.096 12.537 ( -11.342 -2.693 0.000) 11.658 ======================= Grid point 14 (9/56) ======================= q-point: ( 0.27 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.520 ( 16.694 6.287 0.000) 17.839 1.538 ( 16.212 6.835 0.000) 17.594 2.192 ( 14.700 2.354 0.000) 14.887 2.613 ( -0.663 3.214 0.000) 3.281 2.686 ( -0.526 8.745 0.000) 8.761 3.255 ( -0.801 -12.252 0.000) 12.278 3.552 ( 5.750 1.079 0.000) 5.850 3.752 ( 8.567 5.172 0.000) 10.007 4.436 ( 24.042 8.636 0.000) 25.546 6.052 ( -5.643 -3.660 0.000) 6.726 6.236 ( -4.816 12.784 0.000) 13.661 7.646 ( 17.567 2.052 0.000) 17.687 11.111 ( 0.484 -9.419 0.000) 9.431 11.215 ( 0.108 -0.081 0.000) 0.135 12.259 ( -15.975 -2.935 0.000) 16.243 ======================= Grid point 15 (10/56) ======================= q-point: ( 0.36 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.801 ( 11.212 2.838 0.000) 11.566 1.815 ( 11.153 3.496 0.000) 11.688 2.413 ( 7.784 -2.659 0.000) 8.225 2.592 ( -1.298 5.743 0.000) 5.888 2.677 ( -0.390 12.330 0.000) 12.336 3.242 ( -0.505 -19.744 0.000) 19.750 3.668 ( 5.418 0.786 0.000) 5.474 3.902 ( 6.176 6.865 0.000) 9.235 4.850 ( 16.595 11.842 0.000) 20.387 5.943 ( -4.898 -4.414 0.000) 6.593 6.133 ( -5.180 13.384 0.000) 14.351 8.023 ( 18.656 0.391 0.000) 18.660 11.136 ( 2.038 -8.143 0.000) 8.394 11.217 ( 0.084 -0.155 0.000) 0.176 11.923 ( -16.847 -5.025 0.000) 17.580 ======================= Grid point 16 (11/56) ======================= q-point: ( 0.45 0.09 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.957 ( 4.107 0.163 0.000) 4.110 1.970 ( 4.040 1.229 0.000) 4.223 2.513 ( 2.365 -3.758 0.000) 4.440 2.569 ( -0.800 8.105 0.000) 8.145 2.672 ( -0.067 15.071 0.000) 15.072 3.235 ( -0.217 -25.097 0.000) 25.098 3.749 ( 2.292 0.758 0.000) 2.414 3.988 ( 2.248 7.176 0.000) 7.519 5.080 ( 6.044 15.274 0.000) 16.427 5.871 ( -1.999 -4.957 0.000) 5.345 6.049 ( -2.546 12.395 0.000) 12.654 8.314 ( 8.536 -0.896 0.000) 8.583 11.190 ( 2.748 -3.472 0.000) 4.428 11.219 ( 0.032 -0.199 0.000) 0.201 11.640 ( -9.359 -9.746 0.000) 13.512 ======================= Grid point 24 (12/56) ======================= q-point: ( 0.18 0.18 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.407 ( 14.071 14.071 0.000) 19.900 1.431 ( 12.160 12.160 0.000) 17.197 2.057 ( 12.016 12.016 0.000) 16.994 2.657 ( 2.432 2.432 0.000) 3.439 2.766 ( -0.562 -0.562 0.000) 0.795 3.254 ( 0.839 0.839 0.000) 1.187 3.500 ( 3.368 3.368 0.000) 4.763 3.590 ( 4.104 4.104 0.000) 5.804 4.271 ( 22.151 22.151 0.000) 31.326 6.071 ( -5.463 -5.463 0.000) 7.726 6.470 ( 5.152 5.152 0.000) 7.286 7.478 ( 8.985 8.985 0.000) 12.706 10.969 ( -5.785 -5.785 0.000) 8.181 11.213 ( -0.004 -0.004 0.000) 0.005 12.410 ( -10.154 -10.154 0.000) 14.360 ======================= Grid point 25 (13/56) ======================= q-point: ( 0.27 0.18 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.668 ( 11.466 5.851 0.000) 12.873 1.673 ( 11.818 8.152 0.000) 14.356 2.249 ( 7.040 2.858 0.000) 7.598 2.712 ( 2.908 6.650 0.000) 7.258 2.734 ( -1.444 -7.093 0.000) 7.239 3.218 ( -3.874 9.996 0.000) 10.720 3.585 ( 4.995 2.230 0.000) 5.471 3.755 ( 10.091 -2.557 0.000) 10.409 4.757 ( 24.704 21.793 0.000) 32.942 5.947 ( -6.760 -6.770 0.000) 9.567 6.514 ( -0.540 13.320 0.000) 13.331 7.722 ( 15.249 6.079 0.000) 16.416 10.894 ( -1.646 -11.189 0.000) 11.309 11.213 ( -0.003 -0.136 0.000) 0.136 12.139 ( -16.560 -9.480 0.000) 19.082 ======================= Grid point 26 (14/56) ======================= q-point: ( 0.36 0.18 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.864 ( 7.176 3.008 0.000) 7.781 1.874 ( 8.619 2.053 0.000) 8.860 2.351 ( 3.609 -3.327 0.000) 4.908 2.721 ( 0.002 -14.823 0.000) 14.823 2.766 ( 2.356 11.478 0.000) 11.718 3.130 ( -4.358 14.258 0.000) 14.910 3.690 ( 5.156 1.424 0.000) 5.349 3.932 ( 6.981 -2.445 0.000) 7.397 5.202 ( 18.808 21.876 0.000) 28.850 5.815 ( -5.984 -8.297 0.000) 10.230 6.463 ( -3.994 18.196 0.000) 18.629 8.055 ( 16.598 3.477 0.000) 16.959 10.912 ( 3.747 -13.583 0.000) 14.091 11.213 ( -0.001 -0.258 0.000) 0.258 11.769 ( -19.731 -10.408 0.000) 22.308 ======================= Grid point 27 (15/56) ======================= q-point: ( 0.45 0.18 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.959 ( 2.324 -0.078 0.000) 2.325 1.994 ( 3.097 0.957 0.000) 3.242 2.404 ( 1.647 -6.368 0.000) 6.577 2.722 ( -0.007 -22.784 0.000) 22.784 2.799 ( 0.891 14.585 0.000) 14.612 3.062 ( -2.005 19.126 0.000) 19.231 3.770 ( 2.318 1.275 0.000) 2.645 4.026 ( 2.365 -2.130 0.000) 3.182 5.475 ( 7.603 22.559 0.000) 23.806 5.726 ( -2.487 -9.451 0.000) 9.773 6.384 ( -2.830 20.036 0.000) 20.235 8.314 ( 7.531 1.646 0.000) 7.709 11.053 ( 10.061 -11.038 0.000) 14.935 11.213 ( -0.000 -0.331 0.000) 0.331 11.389 ( -16.787 -14.246 0.000) 22.018 ======================= Grid point 36 (16/56) ======================= q-point: ( 0.27 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.780 ( 5.660 5.660 0.000) 8.005 1.819 ( 6.126 6.126 0.000) 8.663 2.290 ( 1.008 1.008 0.000) 1.425 2.550 ( -8.856 -8.856 0.000) 12.525 2.870 ( 8.717 8.717 0.000) 12.328 3.425 ( 6.116 6.116 0.000) 8.649 3.642 ( 3.540 3.540 0.000) 5.006 3.739 ( 2.898 2.898 0.000) 4.098 5.230 ( 23.927 23.927 0.000) 33.839 5.790 ( -8.578 -8.578 0.000) 12.132 6.726 ( 7.407 7.407 0.000) 10.475 7.905 ( 12.339 12.339 0.000) 17.450 10.701 ( -7.373 -7.373 0.000) 10.427 11.210 ( -0.147 -0.147 0.000) 0.207 11.879 ( -16.309 -16.309 0.000) 23.065 ======================= Grid point 37 (17/56) ======================= q-point: ( 0.36 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.897 ( 5.818 0.409 0.000) 5.833 1.913 ( 3.287 1.811 0.000) 3.753 2.287 ( -1.068 -3.022 0.000) 3.205 2.408 ( -5.420 -13.847 0.000) 14.870 3.034 ( 7.012 14.642 0.000) 16.234 3.447 ( -1.533 14.133 0.000) 14.216 3.724 ( 4.386 1.925 0.000) 4.789 3.855 ( 6.325 -4.372 0.000) 7.689 5.620 ( -7.889 -10.886 0.000) 13.444 5.663 ( 18.417 22.933 0.000) 29.413 6.806 ( 0.939 14.882 0.000) 14.911 8.189 ( 14.601 10.450 0.000) 17.955 10.628 ( 0.797 -13.758 0.000) 13.781 11.208 ( -0.111 -0.278 0.000) 0.300 11.504 ( -20.602 -16.112 0.000) 26.154 ======================= Grid point 38 (18/56) ======================= q-point: ( 0.45 0.27 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.955 ( 0.998 -0.289 0.000) 1.039 1.993 ( 3.045 -1.796 0.000) 3.536 2.273 ( 0.333 -7.223 0.000) 7.230 2.321 ( -3.387 -14.948 0.000) 15.327 3.131 ( 2.501 18.016 0.000) 18.189 3.418 ( -0.931 15.091 0.000) 15.120 3.798 ( 2.318 1.392 0.000) 2.704 3.944 ( 2.312 -5.319 0.000) 5.800 5.501 ( -3.408 -12.824 0.000) 13.270 5.930 ( 7.415 22.061 0.000) 23.274 6.794 ( -1.072 19.362 0.000) 19.391 8.418 ( 6.617 9.658 0.000) 11.708 10.752 ( 11.499 -17.819 0.000) 21.207 11.092 ( -19.369 -15.917 0.000) 25.070 11.206 ( -0.041 -0.356 0.000) 0.359 ======================= Grid point 48 (19/56) ======================= q-point: ( 0.36 0.36 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.896 ( -0.426 -0.426 0.000) 0.603 1.938 ( 0.829 0.829 0.000) 1.173 2.196 ( -6.707 -6.707 0.000) 9.485 2.231 ( -2.989 -2.989 0.000) 4.227 3.316 ( 12.513 12.513 0.000) 17.696 3.642 ( 4.202 4.202 0.000) 5.943 3.769 ( 2.624 2.624 0.000) 3.712 3.797 ( 0.015 0.015 0.000) 0.021 5.398 ( -10.682 -10.682 0.000) 15.107 6.075 ( 17.389 17.389 0.000) 24.592 7.026 ( 6.985 6.985 0.000) 9.878 8.461 ( 15.362 15.362 0.000) 21.725 10.410 ( -6.695 -6.695 0.000) 9.468 11.135 ( -20.233 -20.233 0.000) 28.614 11.203 ( -0.211 -0.211 0.000) 0.298 ======================= Grid point 49 (20/56) ======================= q-point: ( 0.45 0.36 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.884 ( -0.539 -8.034 0.000) 8.052 1.948 ( 0.230 -0.289 0.000) 0.370 2.125 ( -0.923 -3.159 0.000) 3.291 2.180 ( -1.840 -4.168 0.000) 4.556 3.492 ( 4.495 17.235 0.000) 17.811 3.650 ( -0.814 8.204 0.000) 8.244 3.823 ( 2.166 1.105 0.000) 2.431 3.837 ( 1.810 -4.873 0.000) 5.198 5.230 ( -5.003 -13.583 0.000) 14.475 6.324 ( 6.757 16.557 0.000) 17.883 7.103 ( 1.456 10.872 0.000) 10.969 8.713 ( 7.572 19.135 0.000) 20.578 10.406 ( 6.348 -15.187 0.000) 16.460 10.737 ( -17.556 -19.372 0.000) 26.143 11.200 ( -0.078 -0.269 0.000) 0.280 ======================= Grid point 60 (21/56) ======================= q-point: ( 0.45 0.45 0.00) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 216 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.760 ( -3.366 -3.366 0.000) 4.760 1.944 ( -0.118 -0.118 0.000) 0.166 2.102 ( -0.043 -0.043 0.000) 0.061 2.121 ( -1.806 -1.806 0.000) 2.554 3.750 ( 1.270 1.270 0.000) 1.796 3.765 ( 8.020 8.020 0.000) 11.342 3.770 ( -0.746 -0.746 0.000) 1.054 3.842 ( 0.962 0.962 0.000) 1.361 5.004 ( -7.264 -7.264 0.000) 10.273 6.560 ( 6.366 6.366 0.000) 9.004 7.236 ( 2.962 2.962 0.000) 4.189 9.083 ( 13.531 13.531 0.000) 19.136 10.209 ( -2.839 -2.839 0.000) 4.014 10.362 ( -15.736 -15.736 0.000) 22.254 11.196 ( -0.100 -0.100 0.000) 0.141 ======================= Grid point 133 (22/56) ======================= q-point: ( 0.09 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 0.948 ( 14.048 14.048 14.048) 24.332 0.948 ( 14.048 14.048 14.048) 24.332 1.496 ( 20.466 20.466 20.466) 35.449 2.666 ( 1.397 1.397 1.397) 2.420 2.666 ( 1.397 1.397 1.397) 2.420 3.241 ( -3.504 -3.504 -3.504) 6.069 3.512 ( 5.692 5.692 5.692) 9.859 3.512 ( 5.692 5.692 5.692) 9.859 3.574 ( 5.507 5.507 5.507) 9.539 6.299 ( 0.753 0.753 0.753) 1.305 6.299 ( 0.753 0.753 0.753) 1.305 7.238 ( 4.438 4.438 4.438) 7.686 11.148 ( -2.099 -2.099 -2.099) 3.636 11.148 ( -2.099 -2.099 -2.099) 3.636 12.670 ( -4.057 -4.057 -4.057) 7.026 ======================= Grid point 134 (23/56) ======================= q-point: ( 0.18 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.251 ( 16.537 9.083 9.083) 20.939 1.280 ( 16.934 8.624 8.624) 20.869 1.915 ( 18.763 9.483 9.483) 23.064 2.646 ( -1.146 1.747 1.747) 2.723 2.728 ( 2.782 2.736 2.736) 4.766 3.187 ( -2.060 -6.184 -6.184) 8.985 3.486 ( -1.183 6.166 6.166) 8.800 3.715 ( 11.054 7.170 7.170) 15.001 3.931 ( 24.756 3.129 3.129) 25.148 6.233 ( -4.233 2.403 2.403) 5.428 6.324 ( -1.329 3.906 3.906) 5.681 7.384 ( 10.439 3.063 3.063) 11.302 11.090 ( -1.626 -4.170 -4.170) 6.117 11.147 ( -0.068 -3.270 -3.270) 4.625 12.517 ( -11.118 -2.913 -2.913) 11.856 ======================= Grid point 135 (24/56) ======================= q-point: ( 0.27 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.573 ( 15.074 5.182 5.182) 16.761 1.599 ( 14.395 4.972 4.972) 16.021 2.212 ( 11.245 1.677 1.677) 11.493 2.623 ( -1.158 1.506 1.506) 2.425 2.788 ( 3.144 7.272 7.272) 10.754 3.154 ( -1.278 -9.064 -9.064) 12.882 3.473 ( -0.281 0.463 0.463) 0.712 3.940 ( 11.053 9.855 9.855) 17.787 4.458 ( 25.861 4.963 4.963) 26.797 6.136 ( -5.219 2.529 2.529) 6.327 6.259 ( -4.804 6.562 6.562) 10.450 7.659 ( 16.798 1.510 1.510) 16.933 11.079 ( 0.577 -5.766 -5.766) 8.175 11.145 ( -0.079 -3.459 -3.459) 4.892 12.238 ( -16.338 -2.717 -2.717) 16.784 ======================= Grid point 136 (25/56) ======================= q-point: ( 0.36 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.835 ( 10.679 2.784 2.784) 11.382 1.837 ( 9.108 2.018 2.018) 9.544 2.384 ( 6.389 -2.563 -2.563) 7.346 2.600 ( -1.112 2.330 2.330) 3.478 2.847 ( 2.604 12.990 12.990) 18.554 3.134 ( -0.818 -12.047 -12.047) 17.056 3.470 ( -0.119 -6.056 -6.056) 8.565 4.140 ( 8.523 11.803 11.803) 18.742 4.909 ( 18.454 7.413 7.413) 21.225 6.035 ( -4.613 2.588 2.588) 5.889 6.145 ( -6.154 6.820 6.820) 11.441 8.021 ( 17.970 0.092 0.092) 17.971 11.117 ( 3.421 -5.233 -5.233) 8.154 11.144 ( -0.064 -3.638 -3.638) 5.145 11.886 ( -18.167 -3.808 -3.808) 18.948 ======================= Grid point 137 (26/56) ======================= q-point: ( 0.45 0.09 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.957 ( 2.968 -0.040 -0.040) 2.969 1.985 ( 3.988 1.371 1.371) 4.434 2.472 ( 2.305 -3.898 -3.898) 5.975 2.581 ( -0.575 3.645 3.645) 5.186 2.885 ( 1.004 17.303 17.303) 24.490 3.122 ( -0.306 -14.185 -14.185) 20.063 3.467 ( -0.106 -10.982 -10.982) 15.531 4.262 ( 3.291 12.866 12.866) 18.490 5.172 ( 7.242 9.964 9.964) 15.843 5.966 ( -1.918 2.562 2.562) 4.100 6.040 ( -3.461 5.586 5.586) 8.625 8.302 ( 8.202 -0.938 -0.938) 8.308 11.143 ( -0.025 -3.747 -3.747) 5.299 11.221 ( 6.356 -1.342 -1.342) 6.633 11.557 ( -12.841 -7.655 -7.655) 16.796 ======================= Grid point 145 (27/56) ======================= q-point: ( 0.18 0.18 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.472 ( 12.102 12.102 5.857) 18.089 1.482 ( 10.999 10.999 4.971) 16.330 2.112 ( 8.766 8.766 4.576) 13.215 2.680 ( 1.352 1.352 1.450) 2.400 2.752 ( -0.046 -0.046 -1.448) 1.449 3.098 ( -2.549 -2.549 -9.869) 10.507 3.706 ( 5.673 5.673 14.773) 16.811 3.715 ( 5.291 5.291 13.370) 15.322 4.227 ( 22.602 22.602 -3.550) 32.161 6.188 ( -4.526 -4.526 10.019) 11.889 6.449 ( 5.278 5.278 -2.005) 7.729 7.492 ( 8.114 8.114 1.515) 11.575 10.967 ( -5.809 -5.809 -0.218) 8.218 11.125 ( -0.309 -0.309 -7.672) 7.684 12.391 ( -9.898 -9.898 -2.570) 14.232 ======================= Grid point 146 (28/56) ======================= q-point: ( 0.27 0.18 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.700 ( 10.716 5.713 3.168) 12.550 1.710 ( 10.746 6.907 3.246) 13.180 2.246 ( 4.896 1.461 -0.218) 5.114 2.670 ( -1.452 1.678 -2.061) 3.028 2.782 ( 2.278 -6.849 1.703) 7.416 3.078 ( 0.748 2.137 -12.788) 12.987 3.624 ( -4.651 13.610 10.691) 17.921 4.011 ( 14.173 -1.792 17.693) 22.740 4.726 ( 25.459 20.589 -2.115) 32.810 6.080 ( -6.070 -5.575 11.780) 14.377 6.497 ( -0.509 13.136 -1.467) 13.227 7.720 ( 14.636 5.116 -0.039) 15.504 10.888 ( -1.887 -10.761 -0.703) 10.947 11.121 ( -0.129 -0.487 -8.226) 8.242 12.126 ( -16.212 -8.822 -1.833) 18.548 ======================= Grid point 147 (29/56) ======================= q-point: ( 0.36 0.18 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.877 ( 6.351 1.949 1.195) 6.750 1.899 ( 8.144 2.984 2.341) 8.983 2.323 ( 3.172 -3.217 -2.390) 5.111 2.656 ( 0.290 -2.839 -3.468) 4.491 2.794 ( -1.369 -8.024 -0.195) 8.143 3.138 ( 5.072 7.726 -8.470) 12.536 3.523 ( -5.445 11.090 4.062) 13.005 4.256 ( 10.033 0.405 19.580) 22.005 5.186 ( 19.600 19.282 -0.575) 27.501 5.958 ( -5.815 -6.871 12.926) 15.751 6.442 ( -4.411 18.197 -1.850) 18.815 8.044 ( 16.295 2.797 -0.899) 16.558 10.900 ( 3.389 -13.461 -1.261) 13.938 11.121 ( 0.130 -0.991 -8.440) 8.499 11.764 ( -19.393 -8.554 -1.095) 21.223 ======================= Grid point 148 (30/56) ======================= q-point: ( 0.45 0.18 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.957 ( 1.870 0.024 0.021) 1.870 2.014 ( 3.033 1.253 1.736) 3.713 2.373 ( 1.525 -5.603 -2.856) 6.471 2.682 ( 2.223 -8.863 -3.531) 9.796 2.743 ( -3.171 -4.382 -2.982) 6.177 3.258 ( 6.251 13.705 3.614) 15.490 3.407 ( -5.923 6.047 -8.522) 12.011 4.398 ( 3.767 1.457 20.805) 21.193 5.478 ( 8.542 19.227 0.972) 21.061 5.864 ( -3.042 -8.178 13.063) 15.709 6.355 ( -3.143 20.424 -2.783) 20.851 8.299 ( 7.428 1.374 -1.260) 7.658 11.017 ( 6.513 -10.439 -3.198) 12.713 11.141 ( 3.492 -6.181 -7.458) 10.297 11.389 ( -16.653 -7.590 -0.001) 18.301 ======================= Grid point 157 (31/56) ======================= q-point: ( 0.27 0.27 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.806 ( 5.354 5.354 2.506) 7.976 1.846 ( 5.953 5.953 2.458) 8.771 2.266 ( 0.475 0.475 -2.099) 2.203 2.544 ( -8.949 -8.949 -0.551) 12.668 2.810 ( 5.263 5.263 -5.947) 9.527 3.128 ( 2.797 2.797 -18.147) 18.574 3.922 ( 4.669 4.669 17.535) 18.737 3.947 ( 5.985 5.985 18.990) 20.791 5.183 ( 23.219 23.219 -3.747) 33.050 5.946 ( -7.379 -7.379 14.073) 17.520 6.713 ( 7.587 7.587 -1.369) 10.817 7.885 ( 11.637 11.637 -1.847) 16.560 10.699 ( -7.296 -7.296 -0.125) 10.319 11.111 ( -0.608 -0.608 -8.944) 8.985 11.884 ( -15.114 -15.114 -0.215) 21.375 ======================= Grid point 158 (32/56) ======================= q-point: ( 0.36 0.27 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.908 ( 4.258 1.010 1.353) 4.580 1.943 ( 4.124 1.387 2.556) 5.046 2.264 ( -0.485 -2.607 -2.053) 3.354 2.397 ( -5.683 -13.599 -1.009) 14.773 2.897 ( 2.926 9.337 -10.762) 14.545 3.209 ( 5.427 1.188 -20.651) 21.385 3.825 ( -5.557 16.096 17.650) 24.525 4.210 ( 11.285 -3.923 21.468) 24.569 5.589 ( 15.171 16.864 -2.880) 22.866 5.811 ( -4.442 -4.886 14.047) 15.522 6.792 ( 0.753 15.369 -1.397) 15.451 8.162 ( 14.637 9.676 -2.451) 17.716 10.623 ( 0.286 -13.268 -0.479) 13.280 11.090 ( -1.997 -2.506 -9.801) 10.311 11.546 ( -17.614 -12.747 2.354) 21.869 ======================= Grid point 159 (33/56) ======================= q-point: ( 0.45 0.27 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.957 ( 0.953 -0.214 0.436) 1.070 2.019 ( 2.618 -1.578 2.198) 3.765 2.260 ( 0.437 -5.975 -1.242) 6.118 2.310 ( -3.134 -15.049 -1.050) 15.408 2.929 ( 0.595 11.182 -13.194) 17.305 3.316 ( 3.915 -0.619 -25.199) 25.509 3.728 ( -3.361 17.783 20.072) 27.026 4.371 ( 4.323 -3.126 23.052) 23.662 5.654 ( -2.266 -10.781 14.515) 18.222 5.895 ( 6.222 19.581 -3.051) 20.771 6.777 ( -1.179 20.023 -1.754) 20.134 8.394 ( 6.799 9.153 -2.079) 11.590 10.724 ( 9.414 -17.098 -2.735) 19.709 10.979 ( -10.417 -10.917 -8.007) 17.083 11.273 ( -6.747 -4.076 3.497) 8.624 ======================= Grid point 169 (34/56) ======================= q-point: ( 0.36 0.36 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.908 ( -1.070 -1.070 1.150) 1.901 1.963 ( 0.803 0.803 2.321) 2.584 2.193 ( -6.301 -6.301 -0.337) 8.917 2.219 ( -2.289 -2.289 -1.101) 3.419 3.073 ( 6.823 6.823 -18.060) 20.476 3.233 ( 2.056 2.056 -23.280) 23.461 4.082 ( 3.665 3.665 21.964) 22.568 4.153 ( 4.040 4.040 20.838) 21.608 5.563 ( -8.815 -8.815 15.488) 19.882 6.027 ( 15.559 15.559 -4.402) 22.440 7.018 ( 7.102 7.102 -0.748) 10.072 8.425 ( 15.257 15.257 -3.369) 21.838 10.414 ( -6.546 -6.546 0.349) 9.263 10.977 ( -10.474 -10.474 -9.864) 17.796 11.309 ( -8.737 -8.737 4.634) 13.196 ======================= Grid point 170 (35/56) ======================= q-point: ( 0.45 0.36 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.881 ( -1.042 -8.372 -0.257) 8.440 1.973 ( 0.166 -0.373 2.230) 2.267 2.131 ( -0.602 -2.567 0.594) 2.703 2.179 ( -1.478 -3.492 -0.018) 3.791 3.140 ( 0.862 9.190 -21.103) 23.033 3.298 ( 2.883 -1.028 -26.069) 26.248 4.041 ( -1.992 13.317 22.223) 25.984 4.302 ( 4.871 -3.291 22.515) 23.270 5.418 ( -4.454 -11.828 18.041) 22.027 6.254 ( 6.171 15.125 -6.895) 17.731 7.095 ( 1.426 11.175 -0.739) 11.290 8.681 ( 7.844 18.702 -3.018) 20.503 10.400 ( 5.047 -14.017 -0.635) 14.912 10.694 ( -14.436 -17.095 -3.823) 22.699 11.219 ( -1.495 -1.589 1.329) 2.555 ======================= Grid point 181 (36/56) ======================= q-point: ( 0.45 0.45 0.09) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.751 ( -3.561 -3.561 -0.973) 5.128 1.972 ( 0.035 0.035 2.661) 2.661 2.112 ( 0.024 0.024 1.035) 1.036 2.130 ( -1.559 -1.559 0.880) 2.374 3.268 ( 2.419 2.419 -25.592) 25.820 3.285 ( 0.616 0.616 -25.924) 25.938 4.236 ( 3.398 3.398 22.386) 22.896 4.261 ( 1.346 1.346 21.430) 21.514 5.223 ( -6.139 -6.139 21.010) 22.733 6.469 ( 5.751 5.751 -8.992) 12.125 7.231 ( 2.988 2.988 -0.475) 4.253 9.041 ( 13.145 13.145 -3.796) 18.973 10.218 ( -2.756 -2.756 0.856) 3.990 10.354 ( -14.280 -14.280 -0.661) 20.206 11.201 ( -0.414 -0.414 0.113) 0.597 ======================= Grid point 266 (37/56) ======================= q-point: ( 0.18 0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.617 ( 8.227 8.227 8.227) 14.249 1.617 ( 8.227 8.227 8.227) 14.249 2.198 ( 3.288 3.288 3.288) 5.695 2.706 ( -0.644 -0.644 -0.644) 1.115 2.706 ( -0.644 -0.644 -0.644) 1.115 2.928 ( -6.859 -6.859 -6.859) 11.880 4.017 ( 10.493 10.493 10.493) 18.174 4.017 ( 10.493 10.493 10.493) 18.174 4.254 ( 16.404 16.404 16.404) 28.412 6.392 ( 2.537 2.537 2.537) 4.394 6.392 ( 2.537 2.537 2.537) 4.394 7.556 ( 5.539 5.539 5.539) 9.594 10.961 ( -4.038 -4.038 -4.038) 6.994 10.961 ( -4.038 -4.038 -4.038) 6.994 12.281 ( -8.767 -8.767 -8.767) 15.184 ======================= Grid point 267 (38/56) ======================= q-point: ( 0.27 0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.784 ( 7.863 4.603 4.603) 10.208 1.793 ( 8.900 5.221 5.221) 11.565 2.238 ( 1.222 -0.671 -0.671) 1.547 2.636 ( -3.499 -1.512 -1.512) 4.101 2.772 ( 4.518 -5.702 -5.702) 9.243 2.819 ( -4.361 -7.900 -7.900) 11.993 3.958 ( -5.631 19.742 19.742) 28.482 4.321 ( 14.241 9.058 9.058) 19.154 4.761 ( 27.039 9.601 9.601) 30.257 6.308 ( -4.568 6.167 6.167) 9.845 6.496 ( 1.284 4.756 4.756) 6.847 7.737 ( 12.763 2.369 2.369) 13.196 10.844 ( -3.609 -5.241 -5.241) 8.244 10.953 ( -0.409 -5.947 -5.947) 8.420 12.041 ( -14.962 -7.107 -7.107) 18.025 ======================= Grid point 268 (39/56) ======================= q-point: ( 0.36 0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.910 ( 4.461 2.058 2.058) 5.327 1.955 ( 6.900 2.951 2.951) 8.064 2.267 ( 1.883 -2.983 -2.983) 4.619 2.582 ( -1.587 -5.715 -5.715) 8.236 2.741 ( -3.647 -6.060 -6.060) 9.314 2.869 ( 4.795 -9.414 -9.414) 14.150 3.841 ( -5.628 20.685 20.685) 29.789 4.571 ( 10.319 9.554 9.554) 17.002 5.245 ( 20.656 8.497 8.497) 23.897 6.218 ( -4.128 6.741 6.741) 10.388 6.450 ( -5.196 7.868 7.868) 12.281 8.036 ( 15.615 0.783 0.783) 15.654 10.828 ( 2.166 -8.236 -8.236) 11.847 10.945 ( -0.317 -6.010 -6.010) 8.505 11.704 ( -18.017 -5.260 -5.260) 19.492 ======================= Grid point 269 (40/56) ======================= q-point: ( 0.45 0.18 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.963 ( 1.102 0.582 0.582) 1.376 2.053 ( 2.583 1.862 1.862) 3.689 2.304 ( 1.238 -3.797 -3.797) 5.510 2.574 ( 0.425 -10.807 -10.807) 15.289 2.677 ( -2.271 -2.321 -2.321) 3.991 2.944 ( 2.233 -12.540 -12.540) 17.875 3.753 ( -2.604 22.032 22.032) 31.267 4.716 ( 3.884 9.654 9.654) 14.194 5.556 ( 9.179 7.916 7.916) 14.477 6.156 ( -1.756 7.143 7.143) 10.254 6.335 ( -4.575 8.925 8.925) 13.425 8.283 ( 7.259 0.442 0.442) 7.286 10.923 ( 6.082 -11.982 -11.982) 18.003 10.941 ( -0.119 -6.044 -6.044) 8.548 11.373 ( -12.642 -2.069 -2.069) 12.976 ======================= Grid point 278 (41/56) ======================= q-point: ( 0.27 0.27 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.876 ( 4.535 4.535 4.284) 7.713 1.905 ( 5.498 5.498 3.224) 8.418 2.220 ( -0.874 -0.874 -2.293) 2.605 2.526 ( -9.378 -9.378 -1.396) 13.335 2.665 ( -0.075 -0.075 -7.326) 7.327 2.784 ( -0.269 -0.269 -15.793) 15.798 4.318 ( 6.835 6.835 17.724) 20.188 4.353 ( 6.150 6.150 22.237) 23.877 5.117 ( 21.694 21.694 -1.610) 30.722 6.286 ( -4.613 -4.613 18.415) 19.537 6.673 ( 8.098 8.098 -2.489) 11.719 7.849 ( 9.621 9.621 -0.985) 13.641 10.696 ( -7.081 -7.081 -0.192) 10.016 10.900 ( -1.254 -1.254 -11.237) 11.377 11.839 ( -12.816 -12.816 -4.777) 18.743 ======================= Grid point 279 (42/56) ======================= q-point: ( 0.36 0.27 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.951 ( 2.522 1.563 2.786) 4.070 2.002 ( 4.211 1.441 3.092) 5.419 2.215 ( 0.449 -1.915 -2.623) 3.279 2.365 ( -6.304 -12.584 -2.472) 14.290 2.674 ( 0.339 0.147 -10.649) 10.656 2.794 ( 1.216 -1.428 -19.460) 19.550 4.238 ( -4.876 18.005 21.123) 28.180 4.623 ( 11.459 -2.430 18.734) 22.095 5.520 ( 17.476 16.587 -3.224) 24.309 6.180 ( -5.540 -5.820 21.450) 22.905 6.753 ( 0.275 16.344 -2.211) 16.496 8.110 ( 14.888 7.560 -2.028) 16.820 10.607 ( -1.387 -11.372 -1.324) 11.533 10.873 ( -1.706 -3.194 -10.950) 11.533 11.552 ( -14.979 -9.619 -2.258) 17.944 ======================= Grid point 280 (43/56) ======================= q-point: ( 0.45 0.27 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.976 ( 0.385 0.293 1.505) 1.581 2.069 ( 2.025 -1.456 2.459) 3.502 2.230 ( 0.821 -3.642 -1.547) 4.041 2.276 ( -2.648 -13.532 -2.586) 14.029 2.676 ( -0.036 1.520 -11.723) 11.821 2.819 ( 0.907 -2.746 -23.107) 23.287 4.163 ( -2.122 18.320 21.580) 28.387 4.783 ( 4.240 -1.547 17.441) 18.015 5.785 ( 8.102 13.397 -4.450) 16.277 6.081 ( -3.698 -7.342 23.987) 25.356 6.728 ( -1.495 21.586 -2.824) 21.821 8.354 ( 7.362 7.726 -1.180) 10.737 10.641 ( 3.812 -14.381 -5.775) 15.959 10.819 ( -3.615 -7.803 -7.432) 11.366 11.305 ( -7.652 -4.387 -0.662) 8.845 ======================= Grid point 290 (44/56) ======================= q-point: ( 0.36 0.36 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.939 ( -3.101 -3.101 1.823) 4.749 2.022 ( 0.833 0.833 3.262) 3.468 2.182 ( -4.813 -4.813 -0.759) 6.848 2.194 ( -0.557 -0.557 -1.233) 1.463 2.699 ( 1.587 1.587 -18.302) 18.439 2.778 ( -0.227 -0.227 -21.745) 21.747 4.541 ( 4.175 4.175 17.740) 18.697 4.544 ( 3.556 3.556 22.416) 22.973 5.819 ( 11.524 11.524 -3.161) 16.601 6.060 ( -4.867 -4.867 20.053) 21.202 6.997 ( 7.440 7.440 -1.348) 10.608 8.348 ( 15.131 15.131 -3.456) 21.676 10.424 ( -6.138 -6.138 0.637) 8.704 10.780 ( -6.786 -6.786 -9.591) 13.569 11.350 ( -9.332 -9.332 -0.718) 13.218 ======================= Grid point 291 (45/56) ======================= q-point: ( 0.45 0.36 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.874 ( -2.162 -10.032 -0.503) 10.275 2.027 ( -0.106 0.502 2.967) 3.011 2.150 ( 0.350 -1.229 1.319) 1.836 2.183 ( -0.668 -1.835 0.424) 1.999 2.719 ( 0.397 2.413 -20.556) 20.701 2.783 ( 0.458 -1.232 -24.897) 24.932 4.472 ( -2.012 12.443 19.997) 23.638 4.715 ( 4.615 -4.149 18.566) 19.576 5.900 ( -1.269 -6.651 27.188) 28.018 6.064 ( 3.189 9.905 -10.251) 14.607 7.075 ( 1.345 12.038 -1.306) 12.183 8.615 ( 8.549 17.661 -2.740) 19.812 10.382 ( 1.680 -10.734 -1.071) 10.918 10.609 ( -8.374 -13.130 -4.183) 16.125 11.222 ( -3.155 -3.376 -1.761) 4.945 ======================= Grid point 302 (46/56) ======================= q-point: ( 0.45 0.45 0.18) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.723 ( -4.053 -4.053 -1.687) 5.975 2.042 ( 0.502 0.502 4.098) 4.159 2.141 ( 0.203 0.203 1.793) 1.816 2.156 ( -1.047 -1.047 1.713) 2.264 2.756 ( 0.887 0.887 -25.165) 25.196 2.768 ( -0.196 -0.196 -25.476) 25.478 4.648 ( 1.542 1.542 18.428) 18.557 4.650 ( 1.352 1.352 17.397) 17.501 5.782 ( -3.623 -3.623 32.903) 33.300 6.215 ( 4.147 4.147 -15.792) 16.845 7.218 ( 3.072 3.072 -0.814) 4.420 8.957 ( 12.384 12.384 -3.668) 17.893 10.242 ( -2.533 -2.533 1.499) 3.883 10.342 ( -11.002 -11.002 -0.398) 15.564 11.178 ( -1.040 -1.040 -3.015) 3.355 ======================= Grid point 399 (47/56) ======================= q-point: ( 0.27 0.27 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.969 ( 3.860 3.860 3.860) 6.685 1.969 ( 3.860 3.860 3.860) 6.685 2.178 ( -1.697 -1.697 -1.697) 2.940 2.449 ( -8.688 -8.688 -8.688) 15.049 2.547 ( -4.728 -4.728 -4.728) 8.189 2.547 ( -4.728 -4.728 -4.728) 8.189 4.644 ( 9.403 9.403 9.403) 16.286 4.644 ( 9.403 9.403 9.403) 16.286 5.214 ( 12.742 12.742 12.742) 22.071 6.617 ( 4.880 4.880 4.880) 8.453 6.617 ( 4.880 4.880 4.880) 8.453 7.891 ( 6.566 6.566 6.566) 11.373 10.692 ( -4.581 -4.581 -4.581) 7.934 10.692 ( -4.581 -4.581 -4.581) 7.934 11.688 ( -10.207 -10.207 -10.207) 17.678 ======================= Grid point 400 (48/56) ======================= q-point: ( 0.36 0.27 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.011 ( 0.933 2.688 2.688) 3.915 2.053 ( 3.228 1.602 1.602) 3.943 2.169 ( 1.065 -1.681 -1.681) 2.605 2.271 ( -7.005 -7.689 -7.689) 12.935 2.486 ( -2.160 -8.838 -8.838) 12.684 2.497 ( -1.605 -7.759 -7.759) 11.089 4.616 ( -3.921 16.167 16.167) 23.197 4.890 ( 9.837 6.328 6.328) 13.298 5.514 ( 14.252 4.028 4.028) 15.349 6.557 ( -2.694 9.683 9.683) 13.956 6.753 ( 1.310 6.872 6.872) 9.806 8.125 ( 15.198 4.827 4.827) 16.660 10.554 ( -4.767 -4.669 -4.669) 8.144 10.689 ( -0.121 -6.484 -6.484) 9.170 11.452 ( -12.454 -7.472 -7.472) 16.333 ======================= Grid point 401 (49/56) ======================= q-point: ( 0.45 0.27 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 2.013 ( -0.190 1.746 1.746) 2.476 2.079 ( -0.327 -4.163 -4.163) 5.897 2.198 ( -0.740 -3.541 -3.541) 5.062 2.203 ( 1.488 -1.123 -1.123) 2.176 2.455 ( -0.953 -11.471 -11.471) 16.250 2.473 ( -0.716 -8.527 -8.527) 12.080 4.549 ( -2.009 16.566 16.566) 23.514 5.026 ( 3.531 5.855 5.855) 9.001 5.724 ( 5.814 -0.163 -0.163) 5.819 6.517 ( -1.115 10.282 10.282) 14.584 6.732 ( -1.634 10.277 10.277) 14.625 8.388 ( 8.267 5.733 5.733) 11.579 10.499 ( -1.017 -7.774 -7.774) 11.041 10.688 ( -0.041 -6.299 -6.299) 8.909 11.246 ( -6.483 -4.946 -4.946) 9.537 ======================= Grid point 411 (50/56) ======================= q-point: ( 0.36 0.36 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.974 ( -6.748 -6.748 1.441) 9.652 2.080 ( 0.977 0.977 2.093) 2.507 2.151 ( -1.428 -1.428 -3.276) 3.848 2.168 ( 1.473 1.473 -1.680) 2.676 2.377 ( -3.556 -3.556 -12.725) 13.683 2.393 ( -3.350 -3.350 -15.308) 16.025 4.855 ( 3.772 3.772 13.557) 14.569 4.921 ( 3.025 3.025 14.198) 14.829 5.658 ( 8.266 8.266 -8.739) 14.595 6.562 ( -3.059 -3.059 24.605) 24.983 6.966 ( 7.945 7.945 -1.620) 11.352 8.333 ( 15.091 15.091 3.319) 21.599 10.438 ( -5.591 -5.591 0.760) 7.944 10.602 ( -3.512 -3.512 -8.136) 9.532 11.276 ( -9.159 -9.159 -6.531) 14.506 ======================= Grid point 412 (51/56) ======================= q-point: ( 0.45 0.36 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 666 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.858 ( -3.374 -11.741 -1.177) 12.273 2.074 ( -0.535 2.693 0.609) 2.812 2.177 ( 2.073 -0.650 0.688) 2.279 2.196 ( 0.583 0.425 0.742) 1.036 2.327 ( -1.285 -3.464 -17.896) 18.273 2.345 ( -1.534 -4.243 -17.062) 17.649 4.818 ( -2.014 10.243 14.386) 17.775 5.031 ( 3.431 -3.530 12.632) 13.557 5.784 ( 3.534 4.871 -13.598) 14.870 6.509 ( -1.650 -3.596 28.050) 28.328 7.047 ( 1.315 12.922 -1.224) 13.046 8.614 ( 9.456 16.453 3.996) 19.392 10.357 ( -2.188 -5.906 -1.387) 6.449 10.535 ( -2.515 -9.454 -3.297) 10.323 11.134 ( -4.064 -5.318 -7.077) 9.741 ======================= Grid point 423 (52/56) ======================= q-point: ( 0.45 0.45 0.27) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.686 ( -4.642 -4.642 -1.910) 6.837 2.120 ( 1.346 1.346 2.602) 3.223 2.170 ( 0.522 0.522 -1.197) 1.406 2.194 ( -0.964 -0.964 1.751) 2.219 2.284 ( -0.682 -0.682 -20.454) 20.477 2.294 ( -1.378 -1.378 -19.314) 19.412 4.952 ( 1.162 1.162 12.671) 12.777 4.964 ( 0.005 0.005 13.127) 13.127 5.865 ( 2.369 2.369 -18.284) 18.588 6.450 ( -1.741 -1.741 31.385) 31.482 7.200 ( 3.216 3.216 -0.909) 4.637 8.941 ( 11.884 11.884 3.315) 17.130 10.275 ( -2.242 -2.242 1.707) 3.601 10.340 ( -7.555 -7.555 0.146) 10.685 11.060 ( -1.865 -1.865 -9.130) 9.504 ======================= Grid point 532 (53/56) ======================= q-point: ( 0.36 0.36 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.947 ( -7.500 -7.500 -7.500) 12.990 2.039 ( -4.332 -4.332 -4.332) 7.503 2.039 ( -4.332 -4.332 -4.332) 7.503 2.150 ( 1.063 1.063 1.063) 1.842 2.240 ( -1.928 -1.928 -1.928) 3.339 2.240 ( -1.928 -1.928 -1.928) 3.339 5.077 ( 4.956 4.956 4.956) 8.584 5.077 ( 4.956 4.956 4.956) 8.584 5.547 ( -1.390 -1.390 -1.390) 2.407 6.935 ( 5.204 5.204 5.204) 9.014 6.935 ( 5.204 5.204 5.204) 9.014 8.519 ( 14.722 14.722 14.722) 25.499 10.452 ( -3.177 -3.177 -3.177) 5.502 10.452 ( -3.177 -3.177 -3.177) 5.502 11.109 ( -9.303 -9.303 -9.303) 16.113 ======================= Grid point 533 (54/56) ======================= q-point: ( 0.45 0.36 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.802 ( -4.328 -6.369 -6.369) 9.993 1.938 ( -2.811 -8.236 -8.236) 11.982 2.035 ( 0.166 -11.426 -11.426) 16.159 2.174 ( 0.917 1.221 1.221) 1.955 2.222 ( -0.440 0.690 0.690) 1.071 2.228 ( -0.072 0.887 0.887) 1.257 5.047 ( -1.972 8.547 8.547) 12.248 5.193 ( 2.947 2.685 2.685) 4.806 5.563 ( 1.277 -6.626 -6.626) 9.458 6.915 ( -0.547 8.867 8.867) 12.552 7.061 ( 2.265 5.310 5.310) 7.844 8.810 ( 10.322 15.298 15.298) 23.971 10.331 ( -4.259 -1.209 -1.209) 4.589 10.461 ( 0.243 -4.767 -4.767) 6.747 10.957 ( -4.585 -9.561 -9.561) 14.278 ======================= Grid point 544 (55/56) ======================= q-point: ( 0.45 0.45 0.36) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 366 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.651 ( -5.049 -5.049 -1.571) 7.311 1.871 ( -2.238 -2.238 -16.336) 16.640 1.877 ( -2.735 -2.735 -15.956) 16.418 2.205 ( 1.191 1.191 2.699) 3.181 2.235 ( 0.260 0.260 1.591) 1.633 2.238 ( 0.143 0.143 1.697) 1.709 5.155 ( 0.895 0.895 7.660) 7.764 5.172 ( -0.494 -0.494 7.627) 7.659 5.524 ( 0.424 0.424 -14.688) 14.700 6.955 ( -0.147 -0.147 17.986) 17.987 7.183 ( 3.386 3.386 -0.713) 4.842 9.128 ( 12.046 12.046 15.367) 22.943 10.306 ( -1.977 -1.977 1.352) 3.106 10.344 ( -4.484 -4.484 0.068) 6.342 10.812 ( -3.797 -3.797 -15.191) 16.112 ======================= Grid point 665 (56/56) ======================= q-point: ( 0.45 0.45 0.45) Boundary mean free path (millimeter): 1000.000 Mass variance parameters: 0.00e+00 0.00e+00 0.00e+00 1.56e-04 3.36e-05 Number of triplets: 146 Calculating Gamma of ph-isotope with tetrahedron method Frequency group velocity (x, y, z) |gv| 1.617 ( -3.034 -3.034 -3.034) 5.255 1.649 ( -4.974 -4.974 -4.974) 8.616 1.649 ( -4.974 -4.974 -4.974) 8.616 2.246 ( 1.318 1.318 1.318) 2.284 2.259 ( 0.519 0.519 0.519) 0.899 2.259 ( 0.519 0.519 0.519) 0.899 5.257 ( 1.323 1.323 1.323) 2.291 5.257 ( 1.323 1.323 1.323) 2.291 5.328 ( -3.260 -3.260 -3.260) 5.647 7.173 ( 2.249 2.249 2.249) 3.895 7.173 ( 2.249 2.249 2.249) 3.895 9.473 ( 14.525 14.525 14.525) 25.159 10.325 ( -1.043 -1.043 -1.043) 1.806 10.325 ( -1.043 -1.043 -1.043) 1.806 10.520 ( -10.171 -10.171 -10.171) 17.616 =================== 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/19965 10.0 456.968 456.968 456.968 0.000 0.000 0.000 3/19965 20.0 69.914 69.914 69.914 0.000 0.000 0.000 3/19965 30.0 29.599 29.599 29.599 0.000 -0.000 0.000 3/19965 40.0 18.590 18.590 18.590 0.000 -0.000 0.000 3/19965 50.0 13.782 13.782 13.782 0.000 -0.000 0.000 3/19965 60.0 11.137 11.137 11.137 0.000 -0.000 0.000 3/19965 70.0 9.460 9.460 9.460 0.000 -0.000 0.000 3/19965 80.0 8.288 8.288 8.288 0.000 -0.000 0.000 3/19965 90.0 7.409 7.409 7.409 0.000 -0.000 0.000 3/19965 100.0 6.719 6.719 6.719 0.000 -0.000 0.000 3/19965 110.0 6.158 6.158 6.158 0.000 -0.000 0.000 3/19965 120.0 5.690 5.690 5.690 0.000 -0.000 0.000 3/19965 130.0 5.293 5.293 5.293 0.000 -0.000 0.000 3/19965 140.0 4.950 4.950 4.950 0.000 -0.000 0.000 3/19965 150.0 4.650 4.650 4.650 0.000 -0.000 0.000 3/19965 160.0 4.386 4.386 4.386 0.000 -0.000 0.000 3/19965 170.0 4.151 4.151 4.151 0.000 -0.000 0.000 3/19965 180.0 3.941 3.941 3.941 0.000 -0.000 0.000 3/19965 190.0 3.750 3.750 3.750 0.000 -0.000 0.000 3/19965 200.0 3.578 3.578 3.578 0.000 -0.000 0.000 3/19965 210.0 3.421 3.421 3.421 0.000 -0.000 0.000 3/19965 220.0 3.276 3.276 3.276 0.000 -0.000 0.000 3/19965 230.0 3.144 3.144 3.144 0.000 -0.000 0.000 3/19965 240.0 3.022 3.022 3.022 0.000 -0.000 0.000 3/19965 250.0 2.909 2.909 2.909 0.000 -0.000 0.000 3/19965 260.0 2.803 2.803 2.803 0.000 -0.000 0.000 3/19965 270.0 2.706 2.706 2.706 0.000 -0.000 0.000 3/19965 280.0 2.614 2.614 2.614 0.000 -0.000 0.000 3/19965 290.0 2.529 2.529 2.529 0.000 -0.000 0.000 3/19965 300.0 2.449 2.449 2.449 0.000 -0.000 0.000 3/19965 310.0 2.374 2.374 2.374 0.000 -0.000 0.000 3/19965 320.0 2.303 2.303 2.303 0.000 -0.000 0.000 3/19965 330.0 2.236 2.236 2.236 0.000 -0.000 0.000 3/19965 340.0 2.173 2.173 2.173 0.000 -0.000 0.000 3/19965 350.0 2.113 2.113 2.113 0.000 -0.000 0.000 3/19965 360.0 2.057 2.057 2.057 0.000 -0.000 0.000 3/19965 370.0 2.003 2.003 2.003 0.000 -0.000 0.000 3/19965 380.0 1.952 1.952 1.952 0.000 -0.000 0.000 3/19965 390.0 1.904 1.904 1.904 0.000 -0.000 0.000 3/19965 400.0 1.858 1.858 1.858 0.000 -0.000 0.000 3/19965 410.0 1.814 1.814 1.814 0.000 -0.000 0.000 3/19965 420.0 1.772 1.772 1.772 0.000 -0.000 0.000 3/19965 430.0 1.732 1.732 1.732 0.000 -0.000 0.000 3/19965 440.0 1.694 1.694 1.694 0.000 -0.000 0.000 3/19965 450.0 1.657 1.657 1.657 0.000 -0.000 0.000 3/19965 460.0 1.622 1.622 1.622 0.000 -0.000 0.000 3/19965 470.0 1.589 1.589 1.589 0.000 -0.000 0.000 3/19965 480.0 1.556 1.556 1.556 0.000 -0.000 0.000 3/19965 490.0 1.525 1.525 1.525 0.000 -0.000 0.000 3/19965 500.0 1.495 1.495 1.495 0.000 -0.000 0.000 3/19965 510.0 1.467 1.467 1.467 0.000 -0.000 0.000 3/19965 520.0 1.439 1.439 1.439 0.000 -0.000 0.000 3/19965 530.0 1.412 1.412 1.412 0.000 -0.000 0.000 3/19965 540.0 1.387 1.387 1.387 0.000 -0.000 0.000 3/19965 550.0 1.362 1.362 1.362 0.000 -0.000 0.000 3/19965 560.0 1.338 1.338 1.338 0.000 -0.000 0.000 3/19965 570.0 1.315 1.315 1.315 0.000 -0.000 0.000 3/19965 580.0 1.293 1.293 1.293 0.000 -0.000 0.000 3/19965 590.0 1.271 1.271 1.271 0.000 -0.000 0.000 3/19965 600.0 1.250 1.250 1.250 0.000 -0.000 0.000 3/19965 610.0 1.230 1.230 1.230 0.000 -0.000 0.000 3/19965 620.0 1.211 1.211 1.211 0.000 -0.000 0.000 3/19965 630.0 1.192 1.192 1.192 0.000 -0.000 0.000 3/19965 640.0 1.173 1.173 1.173 0.000 -0.000 0.000 3/19965 650.0 1.156 1.156 1.156 0.000 -0.000 0.000 3/19965 660.0 1.138 1.138 1.138 0.000 -0.000 0.000 3/19965 670.0 1.122 1.122 1.122 0.000 -0.000 0.000 3/19965 680.0 1.105 1.105 1.105 0.000 -0.000 0.000 3/19965 690.0 1.089 1.089 1.089 0.000 -0.000 0.000 3/19965 700.0 1.074 1.074 1.074 0.000 -0.000 0.000 3/19965 710.0 1.059 1.059 1.059 0.000 -0.000 0.000 3/19965 720.0 1.045 1.045 1.045 0.000 -0.000 0.000 3/19965 730.0 1.030 1.030 1.030 0.000 -0.000 0.000 3/19965 740.0 1.017 1.017 1.017 0.000 -0.000 0.000 3/19965 750.0 1.003 1.003 1.003 0.000 -0.000 0.000 3/19965 760.0 0.990 0.990 0.990 0.000 -0.000 0.000 3/19965 770.0 0.978 0.978 0.978 0.000 -0.000 0.000 3/19965 780.0 0.965 0.965 0.965 0.000 -0.000 0.000 3/19965 790.0 0.953 0.953 0.953 0.000 -0.000 0.000 3/19965 800.0 0.941 0.941 0.941 0.000 -0.000 0.000 3/19965 810.0 0.930 0.930 0.930 0.000 -0.000 0.000 3/19965 820.0 0.918 0.918 0.918 0.000 -0.000 0.000 3/19965 830.0 0.908 0.908 0.908 0.000 -0.000 0.000 3/19965 840.0 0.897 0.897 0.897 0.000 -0.000 0.000 3/19965 850.0 0.886 0.886 0.886 0.000 -0.000 0.000 3/19965 860.0 0.876 0.876 0.876 0.000 -0.000 0.000 3/19965 870.0 0.866 0.866 0.866 0.000 -0.000 0.000 3/19965 880.0 0.856 0.856 0.856 0.000 -0.000 0.000 3/19965 890.0 0.847 0.847 0.847 0.000 -0.000 0.000 3/19965 900.0 0.837 0.837 0.837 0.000 -0.000 0.000 3/19965 910.0 0.828 0.828 0.828 0.000 -0.000 0.000 3/19965 920.0 0.819 0.819 0.819 0.000 -0.000 0.000 3/19965 930.0 0.811 0.811 0.811 0.000 -0.000 0.000 3/19965 940.0 0.802 0.802 0.802 0.000 -0.000 0.000 3/19965 950.0 0.794 0.794 0.794 0.000 -0.000 0.000 3/19965 960.0 0.786 0.786 0.786 0.000 -0.000 0.000 3/19965 970.0 0.777 0.777 0.777 0.000 -0.000 0.000 3/19965 980.0 0.770 0.770 0.770 0.000 -0.000 0.000 3/19965 990.0 0.762 0.762 0.762 0.000 -0.000 0.000 3/19965 1000.0 0.754 0.754 0.754 0.000 -0.000 0.000 3/19965 Thermal conductivity related properties were written into "kappa-m111111.hdf5". Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-09 04:50:04]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_|