----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-09 04:23: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: [2 2 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) Number of symmetry operations in supercell: 144 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.530182459846612 2.038151793480932 6.179755803333327 b -3.530182459846604 2.038151793480932 6.179755803333327 c 0.000000000000003 -4.076303586961870 6.179755803333327 Atomic positions (fractional): *1 Pm 0.64637016383452 0.64637016383452 0.64637016383452 145.000 2 Pm 0.85362983616548 0.85362983616548 0.85362983616548 145.000 3 Pm 0.14637016383452 0.14637016383452 0.14637016383452 145.000 4 Pm 0.35362983616548 0.35362983616548 0.35362983616548 145.000 *5 S 0.25000000000000 0.95456286963452 0.54543713036548 32.065 6 S 0.95456286963452 0.54543713036548 0.25000000000000 32.065 7 S 0.54543713036548 0.25000000000000 0.95456286963452 32.065 8 S 0.45456286963452 0.75000000000000 0.04543713036548 32.065 9 S 0.04543713036548 0.45456286963452 0.75000000000000 32.065 10 S 0.75000000000000 0.04543713036548 0.45456286963452 32.065 -------------------------------- unit cell --------------------------------- Lattice vectors: a 7.060364919693217 0.000000000000000 0.000000000000000 b -3.530182459846608 6.114455380442803 0.000000000000000 c 0.000000000000000 0.000000000000000 18.539267410000001 Atomic positions (fractional): *1 Pm 0.33333333333333 0.66666666666667 0.31303683050119 145.000 > 1 2 Pm 0.66666666666667 0.33333333333333 0.18696316949882 145.000 > 2 3 Pm 0.00000000000000 0.00000000000000 0.14637016383452 145.000 > 3 4 Pm 0.33333333333333 0.66666666666667 0.02029650283215 145.000 > 4 5 Pm 0.00000000000000 0.00000000000000 0.64637016383452 145.000 > 1 6 Pm 0.33333333333333 0.66666666666667 0.52029650283215 145.000 > 2 7 Pm 0.66666666666667 0.33333333333333 0.47970349716785 145.000 > 3 8 Pm 0.00000000000000 0.00000000000000 0.35362983616548 145.000 > 4 9 Pm 0.66666666666667 0.33333333333333 0.97970349716785 145.000 > 1 10 Pm 0.00000000000000 0.00000000000000 0.85362983616548 145.000 > 2 11 Pm 0.33333333333333 0.66666666666667 0.81303683050119 145.000 > 3 12 Pm 0.66666666666667 0.33333333333333 0.68696316949881 145.000 > 4 *13 S 0.00000000000000 0.70456286963452 0.25000000000000 32.065 > 5 14 S 0.70456286963452 0.00000000000000 0.25000000000000 32.065 > 6 15 S 0.29543713036548 0.29543713036548 0.25000000000000 32.065 > 7 16 S 0.37122953630119 0.03789620296785 0.08333333333333 32.065 > 8 17 S 0.96210379703215 0.33333333333333 0.08333333333333 32.065 > 9 18 S 0.66666666666667 0.62877046369881 0.08333333333333 32.065 > 10 19 S 0.66666666666667 0.03789620296785 0.58333333333333 32.065 > 5 20 S 0.37122953630119 0.33333333333333 0.58333333333333 32.065 > 6 21 S 0.96210379703215 0.62877046369881 0.58333333333333 32.065 > 7 22 S 0.03789620296785 0.37122953630119 0.41666666666667 32.065 > 8 23 S 0.62877046369881 0.66666666666667 0.41666666666667 32.065 > 9 24 S 0.33333333333333 0.96210379703215 0.41666666666667 32.065 > 10 25 S 0.33333333333333 0.37122953630119 0.91666666666667 32.065 > 5 26 S 0.03789620296785 0.66666666666667 0.91666666666667 32.065 > 6 27 S 0.62877046369881 0.96210379703215 0.91666666666667 32.065 > 7 28 S 0.70456286963452 0.70456286963452 0.75000000000000 32.065 > 8 29 S 0.29543713036548 0.00000000000000 0.75000000000000 32.065 > 9 30 S 0.00000000000000 0.29543713036548 0.75000000000000 32.065 > 10 -------------------------------- super cell -------------------------------- Lattice vectors: a 14.120729839386433 0.000000000000000 0.000000000000000 b -7.060364919693216 12.228910760885606 0.000000000000000 c 0.000000000000000 0.000000000000000 18.539267410000001 Atomic positions (fractional): *1 Pm 0.16666666666667 0.33333333333333 0.31303683050119 145.000 > 1 2 Pm 0.66666666666667 0.33333333333333 0.31303683050119 145.000 > 1 3 Pm 0.16666666666667 0.83333333333333 0.31303683050119 145.000 > 1 4 Pm 0.66666666666667 0.83333333333333 0.31303683050119 145.000 > 1 5 Pm 0.33333333333333 0.16666666666667 0.18696316949882 145.000 > 2 6 Pm 0.83333333333333 0.16666666666667 0.18696316949882 145.000 > 2 7 Pm 0.33333333333333 0.66666666666667 0.18696316949882 145.000 > 2 8 Pm 0.83333333333333 0.66666666666667 0.18696316949882 145.000 > 2 9 Pm 0.00000000000000 0.00000000000000 0.14637016383452 145.000 > 3 10 Pm 0.50000000000000 0.00000000000000 0.14637016383452 145.000 > 3 11 Pm 0.00000000000000 0.50000000000000 0.14637016383452 145.000 > 3 12 Pm 0.50000000000000 0.50000000000000 0.14637016383452 145.000 > 3 13 Pm 0.16666666666667 0.33333333333333 0.02029650283215 145.000 > 4 14 Pm 0.66666666666667 0.33333333333333 0.02029650283215 145.000 > 4 15 Pm 0.16666666666667 0.83333333333333 0.02029650283215 145.000 > 4 16 Pm 0.66666666666667 0.83333333333333 0.02029650283215 145.000 > 4 17 Pm 0.00000000000000 0.00000000000000 0.64637016383452 145.000 > 1 18 Pm 0.50000000000000 0.00000000000000 0.64637016383452 145.000 > 1 19 Pm 0.00000000000000 0.50000000000000 0.64637016383452 145.000 > 1 20 Pm 0.50000000000000 0.50000000000000 0.64637016383452 145.000 > 1 21 Pm 0.16666666666667 0.33333333333333 0.52029650283215 145.000 > 2 22 Pm 0.66666666666667 0.33333333333333 0.52029650283215 145.000 > 2 23 Pm 0.16666666666667 0.83333333333333 0.52029650283215 145.000 > 2 24 Pm 0.66666666666667 0.83333333333333 0.52029650283215 145.000 > 2 25 Pm 0.33333333333333 0.16666666666667 0.47970349716785 145.000 > 3 26 Pm 0.83333333333333 0.16666666666667 0.47970349716785 145.000 > 3 27 Pm 0.33333333333333 0.66666666666667 0.47970349716785 145.000 > 3 28 Pm 0.83333333333333 0.66666666666667 0.47970349716785 145.000 > 3 29 Pm 0.00000000000000 0.00000000000000 0.35362983616548 145.000 > 4 30 Pm 0.50000000000000 0.00000000000000 0.35362983616548 145.000 > 4 31 Pm 0.00000000000000 0.50000000000000 0.35362983616548 145.000 > 4 32 Pm 0.50000000000000 0.50000000000000 0.35362983616548 145.000 > 4 33 Pm 0.33333333333333 0.16666666666667 0.97970349716785 145.000 > 1 34 Pm 0.83333333333333 0.16666666666667 0.97970349716785 145.000 > 1 35 Pm 0.33333333333333 0.66666666666667 0.97970349716785 145.000 > 1 36 Pm 0.83333333333333 0.66666666666667 0.97970349716785 145.000 > 1 37 Pm 0.00000000000000 0.00000000000000 0.85362983616548 145.000 > 2 38 Pm 0.50000000000000 0.00000000000000 0.85362983616548 145.000 > 2 39 Pm 0.00000000000000 0.50000000000000 0.85362983616548 145.000 > 2 40 Pm 0.50000000000000 0.50000000000000 0.85362983616548 145.000 > 2 41 Pm 0.16666666666667 0.33333333333333 0.81303683050119 145.000 > 3 42 Pm 0.66666666666667 0.33333333333333 0.81303683050119 145.000 > 3 43 Pm 0.16666666666667 0.83333333333333 0.81303683050119 145.000 > 3 44 Pm 0.66666666666667 0.83333333333333 0.81303683050119 145.000 > 3 45 Pm 0.33333333333333 0.16666666666667 0.68696316949881 145.000 > 4 46 Pm 0.83333333333333 0.16666666666667 0.68696316949881 145.000 > 4 47 Pm 0.33333333333333 0.66666666666667 0.68696316949881 145.000 > 4 48 Pm 0.83333333333333 0.66666666666667 0.68696316949881 145.000 > 4 *49 S 0.00000000000000 0.35228143481726 0.25000000000000 32.065 > 5 50 S 0.50000000000000 0.35228143481726 0.25000000000000 32.065 > 5 51 S 0.00000000000000 0.85228143481726 0.25000000000000 32.065 > 5 52 S 0.50000000000000 0.85228143481726 0.25000000000000 32.065 > 5 53 S 0.35228143481726 0.00000000000000 0.25000000000000 32.065 > 6 54 S 0.85228143481726 0.00000000000000 0.25000000000000 32.065 > 6 55 S 0.35228143481726 0.50000000000000 0.25000000000000 32.065 > 6 56 S 0.85228143481726 0.50000000000000 0.25000000000000 32.065 > 6 57 S 0.14771856518274 0.14771856518274 0.25000000000000 32.065 > 7 58 S 0.64771856518274 0.14771856518274 0.25000000000000 32.065 > 7 59 S 0.14771856518274 0.64771856518274 0.25000000000000 32.065 > 7 60 S 0.64771856518274 0.64771856518274 0.25000000000000 32.065 > 7 61 S 0.18561476815059 0.01894810148393 0.08333333333333 32.065 > 8 62 S 0.68561476815059 0.01894810148393 0.08333333333333 32.065 > 8 63 S 0.18561476815059 0.51894810148393 0.08333333333333 32.065 > 8 64 S 0.68561476815059 0.51894810148393 0.08333333333333 32.065 > 8 65 S 0.48105189851607 0.16666666666667 0.08333333333333 32.065 > 9 66 S 0.98105189851607 0.16666666666667 0.08333333333333 32.065 > 9 67 S 0.48105189851607 0.66666666666667 0.08333333333333 32.065 > 9 68 S 0.98105189851607 0.66666666666667 0.08333333333333 32.065 > 9 69 S 0.33333333333333 0.31438523184941 0.08333333333333 32.065 > 10 70 S 0.83333333333333 0.31438523184941 0.08333333333333 32.065 > 10 71 S 0.33333333333333 0.81438523184941 0.08333333333333 32.065 > 10 72 S 0.83333333333333 0.81438523184941 0.08333333333333 32.065 > 10 73 S 0.33333333333333 0.01894810148393 0.58333333333333 32.065 > 5 74 S 0.83333333333333 0.01894810148393 0.58333333333333 32.065 > 5 75 S 0.33333333333333 0.51894810148393 0.58333333333333 32.065 > 5 76 S 0.83333333333333 0.51894810148393 0.58333333333333 32.065 > 5 77 S 0.18561476815059 0.16666666666667 0.58333333333333 32.065 > 6 78 S 0.68561476815059 0.16666666666667 0.58333333333333 32.065 > 6 79 S 0.18561476815059 0.66666666666667 0.58333333333333 32.065 > 6 80 S 0.68561476815059 0.66666666666667 0.58333333333333 32.065 > 6 81 S 0.48105189851607 0.31438523184941 0.58333333333333 32.065 > 7 82 S 0.98105189851607 0.31438523184941 0.58333333333333 32.065 > 7 83 S 0.48105189851607 0.81438523184941 0.58333333333333 32.065 > 7 84 S 0.98105189851607 0.81438523184941 0.58333333333333 32.065 > 7 85 S 0.01894810148393 0.18561476815059 0.41666666666667 32.065 > 8 86 S 0.51894810148393 0.18561476815059 0.41666666666667 32.065 > 8 87 S 0.01894810148393 0.68561476815059 0.41666666666667 32.065 > 8 88 S 0.51894810148393 0.68561476815059 0.41666666666667 32.065 > 8 89 S 0.31438523184941 0.33333333333333 0.41666666666667 32.065 > 9 90 S 0.81438523184941 0.33333333333333 0.41666666666667 32.065 > 9 91 S 0.31438523184941 0.83333333333333 0.41666666666667 32.065 > 9 92 S 0.81438523184941 0.83333333333333 0.41666666666667 32.065 > 9 93 S 0.16666666666667 0.48105189851607 0.41666666666667 32.065 > 10 94 S 0.66666666666667 0.48105189851607 0.41666666666667 32.065 > 10 95 S 0.16666666666667 0.98105189851607 0.41666666666667 32.065 > 10 96 S 0.66666666666667 0.98105189851607 0.41666666666667 32.065 > 10 97 S 0.16666666666667 0.18561476815059 0.91666666666667 32.065 > 5 98 S 0.66666666666667 0.18561476815059 0.91666666666667 32.065 > 5 99 S 0.16666666666667 0.68561476815059 0.91666666666667 32.065 > 5 100 S 0.66666666666667 0.68561476815059 0.91666666666667 32.065 > 5 101 S 0.01894810148393 0.33333333333333 0.91666666666667 32.065 > 6 102 S 0.51894810148393 0.33333333333333 0.91666666666667 32.065 > 6 103 S 0.01894810148393 0.83333333333333 0.91666666666667 32.065 > 6 104 S 0.51894810148393 0.83333333333333 0.91666666666667 32.065 > 6 105 S 0.31438523184941 0.48105189851607 0.91666666666667 32.065 > 7 106 S 0.81438523184941 0.48105189851607 0.91666666666667 32.065 > 7 107 S 0.31438523184941 0.98105189851607 0.91666666666667 32.065 > 7 108 S 0.81438523184941 0.98105189851607 0.91666666666667 32.065 > 7 109 S 0.35228143481726 0.35228143481726 0.75000000000000 32.065 > 8 110 S 0.85228143481726 0.35228143481726 0.75000000000000 32.065 > 8 111 S 0.35228143481726 0.85228143481726 0.75000000000000 32.065 > 8 112 S 0.85228143481726 0.85228143481726 0.75000000000000 32.065 > 8 113 S 0.14771856518274 0.00000000000000 0.75000000000000 32.065 > 9 114 S 0.64771856518274 0.00000000000000 0.75000000000000 32.065 > 9 115 S 0.14771856518274 0.50000000000000 0.75000000000000 32.065 > 9 116 S 0.64771856518274 0.50000000000000 0.75000000000000 32.065 > 9 117 S 0.00000000000000 0.14771856518274 0.75000000000000 32.065 > 10 118 S 0.50000000000000 0.14771856518274 0.75000000000000 32.065 > 10 119 S 0.00000000000000 0.64771856518274 0.75000000000000 32.065 > 10 120 S 0.50000000000000 0.64771856518274 0.75000000000000 32.065 > 10 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 6.6541000 -0.0000000 0.0000000 -0.0000000 6.6541000 0.0000000 0.0000000 0.0000000 6.3895731 -------------------------- Born effective charges -------------------------- 1 Pm 4.3360558 0.2362889 0.0000000 -0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 2 Pm 4.3360558 -0.2362889 0.0000000 0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 3 Pm 4.3360558 -0.2362889 0.0000000 0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 4 Pm 4.3360558 0.2362889 0.0000000 -0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 5 S -2.7493294 0.2448679 -0.7336526 0.2448679 -3.0320784 -0.4235745 -0.9211400 -0.5318204 -2.6377016 6 S -3.1734529 0.0000000 0.0000000 -0.0000000 -2.6079548 0.8471491 0.0000000 1.0636409 -2.6377016 7 S -2.7493294 -0.2448679 0.7336526 -0.2448679 -3.0320784 -0.4235745 0.9211400 -0.5318204 -2.6377016 8 S -2.7493294 -0.2448679 0.7336526 -0.2448679 -3.0320784 -0.4235745 0.9211400 -0.5318204 -2.6377016 9 S -3.1734529 0.0000000 0.0000000 -0.0000000 -2.6079548 0.8471491 0.0000000 1.0636409 -2.6377016 10 S -2.7493294 0.2448679 -0.7336526 0.2448679 -3.0320784 -0.4235745 -0.9211400 -0.5318204 -2.6377016 ---------------------------------------------------------------------------- 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: 360/360 Permutation basis: 7392/7392 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 501 Number of blocks in projector: 268 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 236 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 192 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 73 Use standard eigh solver. Tree of FC basis block matrices: - (501, 493), data: False |-- (73, 71), data: True |-- (192, 190), data: True |-- (236, 232), data: True ----- Solver_atoms: 1 -- 120 / 120 Time (Solver_compr_matrix_reshape): 0.012 Solver_block: 100 / 160 - Time: 0.438 Solver_block: 160 / 160 - Time: 0.302 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.756 --------------------------------- Symfc end -------------------------------- Max drift of force constants: -0.00000000 (zz) -0.00000000 (zz) Permutation basis: 360/360 Permutation basis: 7392/7392 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 501 Number of blocks in projector: 268 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 3 --- Eigsh_solver_block: 1 / 3 --- Block_size: 236 Use standard eigh solver. --- Eigsh_solver_block: 2 / 3 --- Block_size: 192 Use standard eigh solver. --- Eigsh_solver_block: 3 / 3 --- Block_size: 73 Use standard eigh solver. Tree of FC basis block matrices: - (501, 493), data: False |-- (73, 71), data: True |-- (192, 190), data: True |-- (236, 232), data: True Max drift after symmetrization by symfc projector: 0.00000000 (xx) 0.00000000 (xx) Force constants are written into "force_constants.hdf5". ---------------------------------------------------------------------------- One of the following run modes may be specified for phonon calculations. - Mesh sampling (MESH, --mesh) - Q-points (QPOINTS, --qpoints) - Band structure (BAND, --band) - Animation (ANIME, --anime) - Modulation (MODULATION, --modulation) - Characters of Irreps (IRREPS, --irreps) - Create displacements (CREATE_DISPLACEMENTS, -d) ---------------------------------------------------------------------------- Summary of calculation was written in "phonopy.yaml". -------------------------[time 2026-01-09 04:23:57]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:23: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: force constants HDF5 data compression filter: gzip Crystal structure was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". Supercell (dim): [2 2 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.530182459846612 2.038151793480932 6.179755803333327 b -3.530182459846604 2.038151793480932 6.179755803333327 c 0.000000000000003 -4.076303586961870 6.179755803333327 Atomic positions (fractional): 1 Pm 0.64637016383452 0.64637016383452 0.64637016383452 145.000 2 Pm 0.85362983616548 0.85362983616548 0.85362983616548 145.000 3 Pm 0.14637016383452 0.14637016383452 0.14637016383452 145.000 4 Pm 0.35362983616548 0.35362983616548 0.35362983616548 145.000 5 S 0.25000000000000 0.95456286963452 0.54543713036548 32.065 6 S 0.95456286963452 0.54543713036548 0.25000000000000 32.065 7 S 0.54543713036548 0.25000000000000 0.95456286963452 32.065 8 S 0.45456286963452 0.75000000000000 0.04543713036548 32.065 9 S 0.04543713036548 0.45456286963452 0.75000000000000 32.065 10 S 0.75000000000000 0.04543713036548 0.45456286963452 32.065 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.120729839386433 0.000000000000000 0.000000000000000 b -7.060364919693216 12.228910760885606 0.000000000000000 c 0.000000000000000 0.000000000000000 18.539267410000001 Atomic positions (fractional): 1 Pm 0.16666666666667 0.33333333333333 0.31303683050119 145.000 > 1 2 Pm 0.66666666666667 0.33333333333333 0.31303683050119 145.000 > 1 3 Pm 0.16666666666667 0.83333333333333 0.31303683050119 145.000 > 1 4 Pm 0.66666666666667 0.83333333333333 0.31303683050119 145.000 > 1 5 Pm 0.33333333333333 0.16666666666667 0.18696316949882 145.000 > 5 6 Pm 0.83333333333333 0.16666666666667 0.18696316949882 145.000 > 5 7 Pm 0.33333333333333 0.66666666666667 0.18696316949882 145.000 > 5 8 Pm 0.83333333333333 0.66666666666667 0.18696316949882 145.000 > 5 9 Pm 0.00000000000000 0.00000000000000 0.14637016383452 145.000 > 9 10 Pm 0.50000000000000 0.00000000000000 0.14637016383452 145.000 > 9 11 Pm 0.00000000000000 0.50000000000000 0.14637016383452 145.000 > 9 12 Pm 0.50000000000000 0.50000000000000 0.14637016383452 145.000 > 9 13 Pm 0.16666666666667 0.33333333333333 0.02029650283215 145.000 > 13 14 Pm 0.66666666666667 0.33333333333333 0.02029650283215 145.000 > 13 15 Pm 0.16666666666667 0.83333333333333 0.02029650283215 145.000 > 13 16 Pm 0.66666666666667 0.83333333333333 0.02029650283215 145.000 > 13 17 Pm 0.00000000000000 0.00000000000000 0.64637016383452 145.000 > 1 18 Pm 0.50000000000000 0.00000000000000 0.64637016383452 145.000 > 1 19 Pm 0.00000000000000 0.50000000000000 0.64637016383452 145.000 > 1 20 Pm 0.50000000000000 0.50000000000000 0.64637016383452 145.000 > 1 21 Pm 0.16666666666667 0.33333333333333 0.52029650283215 145.000 > 5 22 Pm 0.66666666666667 0.33333333333333 0.52029650283215 145.000 > 5 23 Pm 0.16666666666667 0.83333333333333 0.52029650283215 145.000 > 5 24 Pm 0.66666666666667 0.83333333333333 0.52029650283215 145.000 > 5 25 Pm 0.33333333333333 0.16666666666667 0.47970349716785 145.000 > 9 26 Pm 0.83333333333333 0.16666666666667 0.47970349716785 145.000 > 9 27 Pm 0.33333333333333 0.66666666666667 0.47970349716785 145.000 > 9 28 Pm 0.83333333333333 0.66666666666667 0.47970349716785 145.000 > 9 29 Pm 0.00000000000000 0.00000000000000 0.35362983616548 145.000 > 13 30 Pm 0.50000000000000 0.00000000000000 0.35362983616548 145.000 > 13 31 Pm 0.00000000000000 0.50000000000000 0.35362983616548 145.000 > 13 32 Pm 0.50000000000000 0.50000000000000 0.35362983616548 145.000 > 13 33 Pm 0.33333333333333 0.16666666666667 0.97970349716785 145.000 > 1 34 Pm 0.83333333333333 0.16666666666667 0.97970349716785 145.000 > 1 35 Pm 0.33333333333333 0.66666666666667 0.97970349716785 145.000 > 1 36 Pm 0.83333333333333 0.66666666666667 0.97970349716785 145.000 > 1 37 Pm 0.00000000000000 0.00000000000000 0.85362983616548 145.000 > 5 38 Pm 0.50000000000000 0.00000000000000 0.85362983616548 145.000 > 5 39 Pm 0.00000000000000 0.50000000000000 0.85362983616548 145.000 > 5 40 Pm 0.50000000000000 0.50000000000000 0.85362983616548 145.000 > 5 41 Pm 0.16666666666667 0.33333333333333 0.81303683050119 145.000 > 9 42 Pm 0.66666666666667 0.33333333333333 0.81303683050119 145.000 > 9 43 Pm 0.16666666666667 0.83333333333333 0.81303683050119 145.000 > 9 44 Pm 0.66666666666667 0.83333333333333 0.81303683050119 145.000 > 9 45 Pm 0.33333333333333 0.16666666666667 0.68696316949881 145.000 > 13 46 Pm 0.83333333333333 0.16666666666667 0.68696316949881 145.000 > 13 47 Pm 0.33333333333333 0.66666666666667 0.68696316949881 145.000 > 13 48 Pm 0.83333333333333 0.66666666666667 0.68696316949881 145.000 > 13 49 S 0.00000000000000 0.35228143481726 0.25000000000000 32.065 > 49 50 S 0.50000000000000 0.35228143481726 0.25000000000000 32.065 > 49 51 S 0.00000000000000 0.85228143481726 0.25000000000000 32.065 > 49 52 S 0.50000000000000 0.85228143481726 0.25000000000000 32.065 > 49 53 S 0.35228143481726 0.00000000000000 0.25000000000000 32.065 > 53 54 S 0.85228143481726 0.00000000000000 0.25000000000000 32.065 > 53 55 S 0.35228143481726 0.50000000000000 0.25000000000000 32.065 > 53 56 S 0.85228143481726 0.50000000000000 0.25000000000000 32.065 > 53 57 S 0.14771856518274 0.14771856518274 0.25000000000000 32.065 > 57 58 S 0.64771856518274 0.14771856518274 0.25000000000000 32.065 > 57 59 S 0.14771856518274 0.64771856518274 0.25000000000000 32.065 > 57 60 S 0.64771856518274 0.64771856518274 0.25000000000000 32.065 > 57 61 S 0.18561476815059 0.01894810148393 0.08333333333333 32.065 > 61 62 S 0.68561476815059 0.01894810148393 0.08333333333333 32.065 > 61 63 S 0.18561476815059 0.51894810148393 0.08333333333333 32.065 > 61 64 S 0.68561476815059 0.51894810148393 0.08333333333333 32.065 > 61 65 S 0.48105189851607 0.16666666666667 0.08333333333333 32.065 > 65 66 S 0.98105189851607 0.16666666666667 0.08333333333333 32.065 > 65 67 S 0.48105189851607 0.66666666666667 0.08333333333333 32.065 > 65 68 S 0.98105189851607 0.66666666666667 0.08333333333333 32.065 > 65 69 S 0.33333333333333 0.31438523184941 0.08333333333333 32.065 > 69 70 S 0.83333333333333 0.31438523184941 0.08333333333333 32.065 > 69 71 S 0.33333333333333 0.81438523184941 0.08333333333333 32.065 > 69 72 S 0.83333333333333 0.81438523184941 0.08333333333333 32.065 > 69 73 S 0.33333333333333 0.01894810148393 0.58333333333333 32.065 > 49 74 S 0.83333333333333 0.01894810148393 0.58333333333333 32.065 > 49 75 S 0.33333333333333 0.51894810148393 0.58333333333333 32.065 > 49 76 S 0.83333333333333 0.51894810148393 0.58333333333333 32.065 > 49 77 S 0.18561476815059 0.16666666666667 0.58333333333333 32.065 > 53 78 S 0.68561476815059 0.16666666666667 0.58333333333333 32.065 > 53 79 S 0.18561476815059 0.66666666666667 0.58333333333333 32.065 > 53 80 S 0.68561476815059 0.66666666666667 0.58333333333333 32.065 > 53 81 S 0.48105189851607 0.31438523184941 0.58333333333333 32.065 > 57 82 S 0.98105189851607 0.31438523184941 0.58333333333333 32.065 > 57 83 S 0.48105189851607 0.81438523184941 0.58333333333333 32.065 > 57 84 S 0.98105189851607 0.81438523184941 0.58333333333333 32.065 > 57 85 S 0.01894810148393 0.18561476815059 0.41666666666667 32.065 > 61 86 S 0.51894810148393 0.18561476815059 0.41666666666667 32.065 > 61 87 S 0.01894810148393 0.68561476815059 0.41666666666667 32.065 > 61 88 S 0.51894810148393 0.68561476815059 0.41666666666667 32.065 > 61 89 S 0.31438523184941 0.33333333333333 0.41666666666667 32.065 > 65 90 S 0.81438523184941 0.33333333333333 0.41666666666667 32.065 > 65 91 S 0.31438523184941 0.83333333333333 0.41666666666667 32.065 > 65 92 S 0.81438523184941 0.83333333333333 0.41666666666667 32.065 > 65 93 S 0.16666666666667 0.48105189851607 0.41666666666667 32.065 > 69 94 S 0.66666666666667 0.48105189851607 0.41666666666667 32.065 > 69 95 S 0.16666666666667 0.98105189851607 0.41666666666667 32.065 > 69 96 S 0.66666666666667 0.98105189851607 0.41666666666667 32.065 > 69 97 S 0.16666666666667 0.18561476815059 0.91666666666667 32.065 > 49 98 S 0.66666666666667 0.18561476815059 0.91666666666667 32.065 > 49 99 S 0.16666666666667 0.68561476815059 0.91666666666667 32.065 > 49 100 S 0.66666666666667 0.68561476815059 0.91666666666667 32.065 > 49 101 S 0.01894810148393 0.33333333333333 0.91666666666667 32.065 > 53 102 S 0.51894810148393 0.33333333333333 0.91666666666667 32.065 > 53 103 S 0.01894810148393 0.83333333333333 0.91666666666667 32.065 > 53 104 S 0.51894810148393 0.83333333333333 0.91666666666667 32.065 > 53 105 S 0.31438523184941 0.48105189851607 0.91666666666667 32.065 > 57 106 S 0.81438523184941 0.48105189851607 0.91666666666667 32.065 > 57 107 S 0.31438523184941 0.98105189851607 0.91666666666667 32.065 > 57 108 S 0.81438523184941 0.98105189851607 0.91666666666667 32.065 > 57 109 S 0.35228143481726 0.35228143481726 0.75000000000000 32.065 > 61 110 S 0.85228143481726 0.35228143481726 0.75000000000000 32.065 > 61 111 S 0.35228143481726 0.85228143481726 0.75000000000000 32.065 > 61 112 S 0.85228143481726 0.85228143481726 0.75000000000000 32.065 > 61 113 S 0.14771856518274 0.00000000000000 0.75000000000000 32.065 > 65 114 S 0.64771856518274 0.00000000000000 0.75000000000000 32.065 > 65 115 S 0.14771856518274 0.50000000000000 0.75000000000000 32.065 > 65 116 S 0.64771856518274 0.50000000000000 0.75000000000000 32.065 > 65 117 S 0.00000000000000 0.14771856518274 0.75000000000000 32.065 > 69 118 S 0.50000000000000 0.14771856518274 0.75000000000000 32.065 > 69 119 S 0.00000000000000 0.64771856518274 0.75000000000000 32.065 > 69 120 S 0.50000000000000 0.64771856518274 0.75000000000000 32.065 > 69 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 6.6541000 -0.0000000 0.0000000 -0.0000000 6.6541000 0.0000000 0.0000000 0.0000000 6.3895731 -------------------------- Born effective charges -------------------------- 1 Pm 4.3360558 0.2362889 0.0000000 -0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 2 Pm 4.3360558 -0.2362889 0.0000000 0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 3 Pm 4.3360558 -0.2362889 0.0000000 0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 4 Pm 4.3360558 0.2362889 0.0000000 -0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 5 S -2.7493294 0.2448679 -0.7336526 0.2448679 -3.0320784 -0.4235745 -0.9211400 -0.5318204 -2.6377016 6 S -3.1734529 0.0000000 0.0000000 -0.0000000 -2.6079548 0.8471491 0.0000000 1.0636409 -2.6377016 7 S -2.7493294 -0.2448679 0.7336526 -0.2448679 -3.0320784 -0.4235745 0.9211400 -0.5318204 -2.6377016 8 S -2.7493294 -0.2448679 0.7336526 -0.2448679 -3.0320784 -0.4235745 0.9211400 -0.5318204 -2.6377016 9 S -3.1734529 0.0000000 0.0000000 -0.0000000 -2.6079548 0.8471491 0.0000000 1.0636409 -2.6377016 10 S -2.7493294 0.2448679 -0.7336526 0.2448679 -3.0320784 -0.4235745 -0.9211400 -0.5318204 -2.6377016 ---------------------------------------------------------------------------- Sets of supercell forces were read from "FORCES_FC3.xz". Displacement dataset for fc3 was read from "phono3py_mlp_eval_fc3_disp.yaml.xz". ----------------------------- Force constants ------------------------------ Computing fc3[ 1, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Computing fc3[ 49, x, x ] using numpy.linalg.pinv. Displacements (in Angstrom): [ 0.0100 0.0000 0.0000] [-0.0100 0.0000 0.0000] [ 0.0000 0.0000 0.0100] [ 0.0000 0.0000 -0.0100] Expanding fc3. Symmetrizing fc3 by traditional approach (N=3). Symmetrizing fc2 by traditional approach (N=3). Max drift of fc3: 0.00000000 (xxy) 0.00000000 (xxy) 0.00000000 (xyx) fc3 was written into "fc3.hdf5". Max drift of fc2: 0.00000001 (xy) 0.00000001 (xy) 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:24:05]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-09 04:24:05]------------------------- 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): [2 2 1] Primitive matrix: [ 0.66666667 -0.33333333 -0.33333333] [ 0.33333333 0.33333333 -0.66666667] [0.33333333 0.33333333 0.33333333] Spacegroup: R-3c (167) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 3.530182459846612 2.038151793480932 6.179755803333327 b -3.530182459846604 2.038151793480932 6.179755803333327 c 0.000000000000003 -4.076303586961870 6.179755803333327 Atomic positions (fractional): 1 Pm 0.64637016383452 0.64637016383452 0.64637016383452 145.000 2 Pm 0.85362983616548 0.85362983616548 0.85362983616548 145.000 3 Pm 0.14637016383452 0.14637016383452 0.14637016383452 145.000 4 Pm 0.35362983616548 0.35362983616548 0.35362983616548 145.000 5 S 0.25000000000000 0.95456286963452 0.54543713036548 32.065 6 S 0.95456286963452 0.54543713036548 0.25000000000000 32.065 7 S 0.54543713036548 0.25000000000000 0.95456286963452 32.065 8 S 0.45456286963452 0.75000000000000 0.04543713036548 32.065 9 S 0.04543713036548 0.45456286963452 0.75000000000000 32.065 10 S 0.75000000000000 0.04543713036548 0.45456286963452 32.065 -------------------------------- supercell --------------------------------- Lattice vectors: a 14.120729839386433 0.000000000000000 0.000000000000000 b -7.060364919693216 12.228910760885606 0.000000000000000 c 0.000000000000000 0.000000000000000 18.539267410000001 Atomic positions (fractional): 1 Pm 0.16666666666667 0.33333333333333 0.31303683050119 145.000 > 1 2 Pm 0.66666666666667 0.33333333333333 0.31303683050119 145.000 > 1 3 Pm 0.16666666666667 0.83333333333333 0.31303683050119 145.000 > 1 4 Pm 0.66666666666667 0.83333333333333 0.31303683050119 145.000 > 1 5 Pm 0.33333333333333 0.16666666666667 0.18696316949882 145.000 > 5 6 Pm 0.83333333333333 0.16666666666667 0.18696316949882 145.000 > 5 7 Pm 0.33333333333333 0.66666666666667 0.18696316949882 145.000 > 5 8 Pm 0.83333333333333 0.66666666666667 0.18696316949882 145.000 > 5 9 Pm 0.00000000000000 0.00000000000000 0.14637016383452 145.000 > 9 10 Pm 0.50000000000000 0.00000000000000 0.14637016383452 145.000 > 9 11 Pm 0.00000000000000 0.50000000000000 0.14637016383452 145.000 > 9 12 Pm 0.50000000000000 0.50000000000000 0.14637016383452 145.000 > 9 13 Pm 0.16666666666667 0.33333333333333 0.02029650283215 145.000 > 13 14 Pm 0.66666666666667 0.33333333333333 0.02029650283215 145.000 > 13 15 Pm 0.16666666666667 0.83333333333333 0.02029650283215 145.000 > 13 16 Pm 0.66666666666667 0.83333333333333 0.02029650283215 145.000 > 13 17 Pm 0.00000000000000 0.00000000000000 0.64637016383452 145.000 > 1 18 Pm 0.50000000000000 0.00000000000000 0.64637016383452 145.000 > 1 19 Pm 0.00000000000000 0.50000000000000 0.64637016383452 145.000 > 1 20 Pm 0.50000000000000 0.50000000000000 0.64637016383452 145.000 > 1 21 Pm 0.16666666666667 0.33333333333333 0.52029650283215 145.000 > 5 22 Pm 0.66666666666667 0.33333333333333 0.52029650283215 145.000 > 5 23 Pm 0.16666666666667 0.83333333333333 0.52029650283215 145.000 > 5 24 Pm 0.66666666666667 0.83333333333333 0.52029650283215 145.000 > 5 25 Pm 0.33333333333333 0.16666666666667 0.47970349716785 145.000 > 9 26 Pm 0.83333333333333 0.16666666666667 0.47970349716785 145.000 > 9 27 Pm 0.33333333333333 0.66666666666667 0.47970349716785 145.000 > 9 28 Pm 0.83333333333333 0.66666666666667 0.47970349716785 145.000 > 9 29 Pm 0.00000000000000 0.00000000000000 0.35362983616548 145.000 > 13 30 Pm 0.50000000000000 0.00000000000000 0.35362983616548 145.000 > 13 31 Pm 0.00000000000000 0.50000000000000 0.35362983616548 145.000 > 13 32 Pm 0.50000000000000 0.50000000000000 0.35362983616548 145.000 > 13 33 Pm 0.33333333333333 0.16666666666667 0.97970349716785 145.000 > 1 34 Pm 0.83333333333333 0.16666666666667 0.97970349716785 145.000 > 1 35 Pm 0.33333333333333 0.66666666666667 0.97970349716785 145.000 > 1 36 Pm 0.83333333333333 0.66666666666667 0.97970349716785 145.000 > 1 37 Pm 0.00000000000000 0.00000000000000 0.85362983616548 145.000 > 5 38 Pm 0.50000000000000 0.00000000000000 0.85362983616548 145.000 > 5 39 Pm 0.00000000000000 0.50000000000000 0.85362983616548 145.000 > 5 40 Pm 0.50000000000000 0.50000000000000 0.85362983616548 145.000 > 5 41 Pm 0.16666666666667 0.33333333333333 0.81303683050119 145.000 > 9 42 Pm 0.66666666666667 0.33333333333333 0.81303683050119 145.000 > 9 43 Pm 0.16666666666667 0.83333333333333 0.81303683050119 145.000 > 9 44 Pm 0.66666666666667 0.83333333333333 0.81303683050119 145.000 > 9 45 Pm 0.33333333333333 0.16666666666667 0.68696316949881 145.000 > 13 46 Pm 0.83333333333333 0.16666666666667 0.68696316949881 145.000 > 13 47 Pm 0.33333333333333 0.66666666666667 0.68696316949881 145.000 > 13 48 Pm 0.83333333333333 0.66666666666667 0.68696316949881 145.000 > 13 49 S 0.00000000000000 0.35228143481726 0.25000000000000 32.065 > 49 50 S 0.50000000000000 0.35228143481726 0.25000000000000 32.065 > 49 51 S 0.00000000000000 0.85228143481726 0.25000000000000 32.065 > 49 52 S 0.50000000000000 0.85228143481726 0.25000000000000 32.065 > 49 53 S 0.35228143481726 0.00000000000000 0.25000000000000 32.065 > 53 54 S 0.85228143481726 0.00000000000000 0.25000000000000 32.065 > 53 55 S 0.35228143481726 0.50000000000000 0.25000000000000 32.065 > 53 56 S 0.85228143481726 0.50000000000000 0.25000000000000 32.065 > 53 57 S 0.14771856518274 0.14771856518274 0.25000000000000 32.065 > 57 58 S 0.64771856518274 0.14771856518274 0.25000000000000 32.065 > 57 59 S 0.14771856518274 0.64771856518274 0.25000000000000 32.065 > 57 60 S 0.64771856518274 0.64771856518274 0.25000000000000 32.065 > 57 61 S 0.18561476815059 0.01894810148393 0.08333333333333 32.065 > 61 62 S 0.68561476815059 0.01894810148393 0.08333333333333 32.065 > 61 63 S 0.18561476815059 0.51894810148393 0.08333333333333 32.065 > 61 64 S 0.68561476815059 0.51894810148393 0.08333333333333 32.065 > 61 65 S 0.48105189851607 0.16666666666667 0.08333333333333 32.065 > 65 66 S 0.98105189851607 0.16666666666667 0.08333333333333 32.065 > 65 67 S 0.48105189851607 0.66666666666667 0.08333333333333 32.065 > 65 68 S 0.98105189851607 0.66666666666667 0.08333333333333 32.065 > 65 69 S 0.33333333333333 0.31438523184941 0.08333333333333 32.065 > 69 70 S 0.83333333333333 0.31438523184941 0.08333333333333 32.065 > 69 71 S 0.33333333333333 0.81438523184941 0.08333333333333 32.065 > 69 72 S 0.83333333333333 0.81438523184941 0.08333333333333 32.065 > 69 73 S 0.33333333333333 0.01894810148393 0.58333333333333 32.065 > 49 74 S 0.83333333333333 0.01894810148393 0.58333333333333 32.065 > 49 75 S 0.33333333333333 0.51894810148393 0.58333333333333 32.065 > 49 76 S 0.83333333333333 0.51894810148393 0.58333333333333 32.065 > 49 77 S 0.18561476815059 0.16666666666667 0.58333333333333 32.065 > 53 78 S 0.68561476815059 0.16666666666667 0.58333333333333 32.065 > 53 79 S 0.18561476815059 0.66666666666667 0.58333333333333 32.065 > 53 80 S 0.68561476815059 0.66666666666667 0.58333333333333 32.065 > 53 81 S 0.48105189851607 0.31438523184941 0.58333333333333 32.065 > 57 82 S 0.98105189851607 0.31438523184941 0.58333333333333 32.065 > 57 83 S 0.48105189851607 0.81438523184941 0.58333333333333 32.065 > 57 84 S 0.98105189851607 0.81438523184941 0.58333333333333 32.065 > 57 85 S 0.01894810148393 0.18561476815059 0.41666666666667 32.065 > 61 86 S 0.51894810148393 0.18561476815059 0.41666666666667 32.065 > 61 87 S 0.01894810148393 0.68561476815059 0.41666666666667 32.065 > 61 88 S 0.51894810148393 0.68561476815059 0.41666666666667 32.065 > 61 89 S 0.31438523184941 0.33333333333333 0.41666666666667 32.065 > 65 90 S 0.81438523184941 0.33333333333333 0.41666666666667 32.065 > 65 91 S 0.31438523184941 0.83333333333333 0.41666666666667 32.065 > 65 92 S 0.81438523184941 0.83333333333333 0.41666666666667 32.065 > 65 93 S 0.16666666666667 0.48105189851607 0.41666666666667 32.065 > 69 94 S 0.66666666666667 0.48105189851607 0.41666666666667 32.065 > 69 95 S 0.16666666666667 0.98105189851607 0.41666666666667 32.065 > 69 96 S 0.66666666666667 0.98105189851607 0.41666666666667 32.065 > 69 97 S 0.16666666666667 0.18561476815059 0.91666666666667 32.065 > 49 98 S 0.66666666666667 0.18561476815059 0.91666666666667 32.065 > 49 99 S 0.16666666666667 0.68561476815059 0.91666666666667 32.065 > 49 100 S 0.66666666666667 0.68561476815059 0.91666666666667 32.065 > 49 101 S 0.01894810148393 0.33333333333333 0.91666666666667 32.065 > 53 102 S 0.51894810148393 0.33333333333333 0.91666666666667 32.065 > 53 103 S 0.01894810148393 0.83333333333333 0.91666666666667 32.065 > 53 104 S 0.51894810148393 0.83333333333333 0.91666666666667 32.065 > 53 105 S 0.31438523184941 0.48105189851607 0.91666666666667 32.065 > 57 106 S 0.81438523184941 0.48105189851607 0.91666666666667 32.065 > 57 107 S 0.31438523184941 0.98105189851607 0.91666666666667 32.065 > 57 108 S 0.81438523184941 0.98105189851607 0.91666666666667 32.065 > 57 109 S 0.35228143481726 0.35228143481726 0.75000000000000 32.065 > 61 110 S 0.85228143481726 0.35228143481726 0.75000000000000 32.065 > 61 111 S 0.35228143481726 0.85228143481726 0.75000000000000 32.065 > 61 112 S 0.85228143481726 0.85228143481726 0.75000000000000 32.065 > 61 113 S 0.14771856518274 0.00000000000000 0.75000000000000 32.065 > 65 114 S 0.64771856518274 0.00000000000000 0.75000000000000 32.065 > 65 115 S 0.14771856518274 0.50000000000000 0.75000000000000 32.065 > 65 116 S 0.64771856518274 0.50000000000000 0.75000000000000 32.065 > 65 117 S 0.00000000000000 0.14771856518274 0.75000000000000 32.065 > 69 118 S 0.50000000000000 0.14771856518274 0.75000000000000 32.065 > 69 119 S 0.00000000000000 0.64771856518274 0.75000000000000 32.065 > 69 120 S 0.50000000000000 0.64771856518274 0.75000000000000 32.065 > 69 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 6.6541000 -0.0000000 0.0000000 -0.0000000 6.6541000 0.0000000 0.0000000 0.0000000 6.3895731 -------------------------- Born effective charges -------------------------- 1 Pm 4.3360558 0.2362889 0.0000000 -0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 2 Pm 4.3360558 -0.2362889 0.0000000 0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 3 Pm 4.3360558 -0.2362889 0.0000000 0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 4 Pm 4.3360558 0.2362889 0.0000000 -0.2362889 4.3360558 0.0000000 0.0000000 0.0000000 3.9565523 5 S -2.7493294 0.2448679 -0.7336526 0.2448679 -3.0320784 -0.4235745 -0.9211400 -0.5318204 -2.6377016 6 S -3.1734529 0.0000000 0.0000000 -0.0000000 -2.6079548 0.8471491 0.0000000 1.0636409 -2.6377016 7 S -2.7493294 -0.2448679 0.7336526 -0.2448679 -3.0320784 -0.4235745 0.9211400 -0.5318204 -2.6377016 8 S -2.7493294 -0.2448679 0.7336526 -0.2448679 -3.0320784 -0.4235745 0.9211400 -0.5318204 -2.6377016 9 S -3.1734529 0.0000000 0.0000000 -0.0000000 -2.6079548 0.8471491 0.0000000 1.0636409 -2.6377016 10 S -2.7493294 0.2448679 -0.7336526 0.2448679 -3.0320784 -0.4235745 -0.9211400 -0.5318204 -2.6377016 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000000 (xxy) 0.00000000 (xxy) 0.00000000 (xyx) Max drift of fc2: 0.00000000 (xx) 0.00000000 (xx) --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperature: 0.0 10.0 20.0 30.0 40.0 ... 1000.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 Length for sampling mesh generation: 50.00 Generating grid system ... [ 9 9 9 ] 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.65, Number of G-points: 295, Lambda: 0.17 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) --------------------