----------------------------- ------- calculate fc2 ------- ----------------------------- _ _ __ | |__ ___ _ __ ___ _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | | | |_) | | | | (_) | | | | (_) || |_) | |_| | | .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, | |_| |_| |___/ 2.47.1 -------------------------[time 2026-01-08 08:42:34]------------------------- 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 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P4/nmm (129) Number of symmetry operations in supercell: 288 ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.209387229999999 0.000000000000000 0.000000000000000 b 0.000000000000000 4.209387229999999 0.000000000000000 c 0.000000000000000 0.000000000000000 7.033745940000000 Atomic positions (fractional): *1 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 2 O 0.50000000000000 0.50000000000000 0.00000000000000 15.999 *3 Cl 0.50000000000000 0.00000000000000 0.63129275280323 35.453 4 Cl 0.00000000000000 0.50000000000000 0.36870724719677 35.453 *5 Ac 0.00000000000000 0.50000000000000 0.82042027633109 227.000 6 Ac 0.50000000000000 0.00000000000000 0.17957972366891 227.000 -------------------------------- unit cell --------------------------------- Lattice vectors: a 4.209387229999999 0.000000000000000 0.000000000000000 b 0.000000000000000 4.209387229999999 0.000000000000000 c 0.000000000000000 0.000000000000000 7.033745940000000 Atomic positions (fractional): *1 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 > 1 2 O 0.50000000000000 0.50000000000000 0.00000000000000 15.999 > 2 *3 Cl 0.50000000000000 0.00000000000000 0.63129275280323 35.453 > 3 4 Cl 0.00000000000000 0.50000000000000 0.36870724719677 35.453 > 4 *5 Ac 0.00000000000000 0.50000000000000 0.82042027633109 227.000 > 5 6 Ac 0.50000000000000 0.00000000000000 0.17957972366891 227.000 > 6 -------------------------------- super cell -------------------------------- Lattice vectors: a 12.628161689999997 0.000000000000000 0.000000000000000 b 0.000000000000000 12.628161689999997 0.000000000000000 c 0.000000000000000 0.000000000000000 14.067491880000000 Atomic positions (fractional): *1 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 > 1 2 O 0.33333333333333 0.00000000000000 0.00000000000000 15.999 > 1 3 O 0.66666666666667 0.00000000000000 0.00000000000000 15.999 > 1 4 O 0.00000000000000 0.33333333333333 0.00000000000000 15.999 > 1 5 O 0.33333333333333 0.33333333333333 0.00000000000000 15.999 > 1 6 O 0.66666666666667 0.33333333333333 0.00000000000000 15.999 > 1 7 O 0.00000000000000 0.66666666666667 0.00000000000000 15.999 > 1 8 O 0.33333333333333 0.66666666666667 0.00000000000000 15.999 > 1 9 O 0.66666666666667 0.66666666666667 0.00000000000000 15.999 > 1 10 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 1 11 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 1 12 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 1 13 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 1 14 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 1 15 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 1 16 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 1 17 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 1 18 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 1 19 O 0.16666666666667 0.16666666666667 0.00000000000000 15.999 > 2 20 O 0.50000000000000 0.16666666666667 0.00000000000000 15.999 > 2 21 O 0.83333333333333 0.16666666666667 0.00000000000000 15.999 > 2 22 O 0.16666666666667 0.50000000000000 0.00000000000000 15.999 > 2 23 O 0.50000000000000 0.50000000000000 0.00000000000000 15.999 > 2 24 O 0.83333333333333 0.50000000000000 0.00000000000000 15.999 > 2 25 O 0.16666666666667 0.83333333333333 0.00000000000000 15.999 > 2 26 O 0.50000000000000 0.83333333333333 0.00000000000000 15.999 > 2 27 O 0.83333333333333 0.83333333333333 0.00000000000000 15.999 > 2 28 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 2 29 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 2 30 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 2 31 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 2 32 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 2 33 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 2 34 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 2 35 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 2 36 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 2 *37 Cl 0.16666666666667 0.00000000000000 0.31564637640161 35.453 > 3 38 Cl 0.50000000000000 0.00000000000000 0.31564637640161 35.453 > 3 39 Cl 0.83333333333333 0.00000000000000 0.31564637640161 35.453 > 3 40 Cl 0.16666666666667 0.33333333333333 0.31564637640161 35.453 > 3 41 Cl 0.50000000000000 0.33333333333333 0.31564637640161 35.453 > 3 42 Cl 0.83333333333333 0.33333333333333 0.31564637640161 35.453 > 3 43 Cl 0.16666666666667 0.66666666666667 0.31564637640161 35.453 > 3 44 Cl 0.50000000000000 0.66666666666667 0.31564637640161 35.453 > 3 45 Cl 0.83333333333333 0.66666666666667 0.31564637640161 35.453 > 3 46 Cl 0.16666666666667 0.00000000000000 0.81564637640161 35.453 > 3 47 Cl 0.50000000000000 0.00000000000000 0.81564637640161 35.453 > 3 48 Cl 0.83333333333333 0.00000000000000 0.81564637640161 35.453 > 3 49 Cl 0.16666666666667 0.33333333333333 0.81564637640161 35.453 > 3 50 Cl 0.50000000000000 0.33333333333333 0.81564637640161 35.453 > 3 51 Cl 0.83333333333333 0.33333333333333 0.81564637640161 35.453 > 3 52 Cl 0.16666666666667 0.66666666666667 0.81564637640161 35.453 > 3 53 Cl 0.50000000000000 0.66666666666667 0.81564637640161 35.453 > 3 54 Cl 0.83333333333333 0.66666666666667 0.81564637640161 35.453 > 3 55 Cl 0.00000000000000 0.16666666666667 0.18435362359839 35.453 > 4 56 Cl 0.33333333333333 0.16666666666667 0.18435362359839 35.453 > 4 57 Cl 0.66666666666667 0.16666666666667 0.18435362359839 35.453 > 4 58 Cl 0.00000000000000 0.50000000000000 0.18435362359839 35.453 > 4 59 Cl 0.33333333333333 0.50000000000000 0.18435362359839 35.453 > 4 60 Cl 0.66666666666667 0.50000000000000 0.18435362359839 35.453 > 4 61 Cl 0.00000000000000 0.83333333333333 0.18435362359839 35.453 > 4 62 Cl 0.33333333333333 0.83333333333333 0.18435362359839 35.453 > 4 63 Cl 0.66666666666667 0.83333333333333 0.18435362359839 35.453 > 4 64 Cl 0.00000000000000 0.16666666666667 0.68435362359839 35.453 > 4 65 Cl 0.33333333333333 0.16666666666667 0.68435362359839 35.453 > 4 66 Cl 0.66666666666667 0.16666666666667 0.68435362359839 35.453 > 4 67 Cl 0.00000000000000 0.50000000000000 0.68435362359839 35.453 > 4 68 Cl 0.33333333333333 0.50000000000000 0.68435362359839 35.453 > 4 69 Cl 0.66666666666667 0.50000000000000 0.68435362359839 35.453 > 4 70 Cl 0.00000000000000 0.83333333333333 0.68435362359839 35.453 > 4 71 Cl 0.33333333333333 0.83333333333333 0.68435362359839 35.453 > 4 72 Cl 0.66666666666667 0.83333333333333 0.68435362359839 35.453 > 4 *73 Ac 0.00000000000000 0.16666666666667 0.41021013816555 227.000 > 5 74 Ac 0.33333333333333 0.16666666666667 0.41021013816555 227.000 > 5 75 Ac 0.66666666666667 0.16666666666667 0.41021013816555 227.000 > 5 76 Ac 0.00000000000000 0.50000000000000 0.41021013816555 227.000 > 5 77 Ac 0.33333333333333 0.50000000000000 0.41021013816555 227.000 > 5 78 Ac 0.66666666666667 0.50000000000000 0.41021013816555 227.000 > 5 79 Ac 0.00000000000000 0.83333333333333 0.41021013816555 227.000 > 5 80 Ac 0.33333333333333 0.83333333333333 0.41021013816555 227.000 > 5 81 Ac 0.66666666666667 0.83333333333333 0.41021013816555 227.000 > 5 82 Ac 0.00000000000000 0.16666666666667 0.91021013816555 227.000 > 5 83 Ac 0.33333333333333 0.16666666666667 0.91021013816555 227.000 > 5 84 Ac 0.66666666666667 0.16666666666667 0.91021013816555 227.000 > 5 85 Ac 0.00000000000000 0.50000000000000 0.91021013816555 227.000 > 5 86 Ac 0.33333333333333 0.50000000000000 0.91021013816555 227.000 > 5 87 Ac 0.66666666666667 0.50000000000000 0.91021013816555 227.000 > 5 88 Ac 0.00000000000000 0.83333333333333 0.91021013816555 227.000 > 5 89 Ac 0.33333333333333 0.83333333333333 0.91021013816555 227.000 > 5 90 Ac 0.66666666666667 0.83333333333333 0.91021013816555 227.000 > 5 91 Ac 0.16666666666667 0.00000000000000 0.08978986183445 227.000 > 6 92 Ac 0.50000000000000 0.00000000000000 0.08978986183445 227.000 > 6 93 Ac 0.83333333333333 0.00000000000000 0.08978986183445 227.000 > 6 94 Ac 0.16666666666667 0.33333333333333 0.08978986183445 227.000 > 6 95 Ac 0.50000000000000 0.33333333333333 0.08978986183445 227.000 > 6 96 Ac 0.83333333333333 0.33333333333333 0.08978986183445 227.000 > 6 97 Ac 0.16666666666667 0.66666666666667 0.08978986183445 227.000 > 6 98 Ac 0.50000000000000 0.66666666666667 0.08978986183445 227.000 > 6 99 Ac 0.83333333333333 0.66666666666667 0.08978986183445 227.000 > 6 100 Ac 0.16666666666667 0.00000000000000 0.58978986183445 227.000 > 6 101 Ac 0.50000000000000 0.00000000000000 0.58978986183445 227.000 > 6 102 Ac 0.83333333333333 0.00000000000000 0.58978986183445 227.000 > 6 103 Ac 0.16666666666667 0.33333333333333 0.58978986183445 227.000 > 6 104 Ac 0.50000000000000 0.33333333333333 0.58978986183445 227.000 > 6 105 Ac 0.83333333333333 0.33333333333333 0.58978986183445 227.000 > 6 106 Ac 0.16666666666667 0.66666666666667 0.58978986183445 227.000 > 6 107 Ac 0.50000000000000 0.66666666666667 0.58978986183445 227.000 > 6 108 Ac 0.83333333333333 0.66666666666667 0.58978986183445 227.000 > 6 ---------------------------------------------------------------------------- NAC parameters were read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". --------------------------- Dielectric constant ---------------------------- 4.2012710 0.0000000 0.0000000 0.0000000 4.2012710 0.0000000 0.0000000 0.0000000 4.0787131 -------------------------- Born effective charges -------------------------- 1 O -2.8552025 0.0000000 0.0000000 0.0000000 -2.8552025 0.0000000 0.0000000 0.0000000 -2.2085931 2 O -2.8552025 0.0000000 0.0000000 0.0000000 -2.8552025 0.0000000 0.0000000 0.0000000 -2.2085931 3 Cl -1.3775560 0.0000000 0.0000000 0.0000000 -1.3775560 0.0000000 0.0000000 0.0000000 -1.4905265 4 Cl -1.3775560 0.0000000 0.0000000 0.0000000 -1.3775560 0.0000000 0.0000000 0.0000000 -1.4905265 5 Ac 4.2327586 0.0000000 0.0000000 0.0000000 4.2327586 0.0000000 0.0000000 0.0000000 3.6991195 6 Ac 4.2327586 0.0000000 0.0000000 0.0000000 4.2327586 0.0000000 0.0000000 0.0000000 3.6991195 ---------------------------------------------------------------------------- Displacement-force dataset was read from "phonopy_mlp_eval_fc2_dataset.yaml.xz". -------------------------------- Symfc start ------------------------------- Symfc version 1.5.4 (https://github.com/symfc/symfc) Citation: A. Seko and A. Togo, Phys. Rev. B, 110, 214302 (2024) Computing [2] order force constants. Permutation basis: 324/324 Permutation basis: 3366/3366 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 256 Number of blocks in projector: 262 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 4 --- Eigsh_solver_block: 1 / 4 --- Block_size: 96 Use standard eigh solver. --- Eigsh_solver_block: 2 / 4 --- Block_size: 72 Use standard eigh solver. --- Eigsh_solver_block: 3 / 4 --- Block_size: 64 Use standard eigh solver. --- Eigsh_solver_block: 4 / 4 --- Block_size: 24 Use standard eigh solver. Tree of FC basis block matrices: - (256, 250), data: False |-- (24, 24), data: True |-- (64, 61), data: True |-- (72, 72), data: True |-- (96, 93), data: True ----- Solver_atoms: 1 -- 108 / 108 Time (Solver_compr_matrix_reshape): 0.008 Solver_block: 80 / 80 - Time: 0.202 Solver: Calculate X.T @ X and X.T @ y (disp @ compr @ eigvecs).T @ (disp @ compr @ eigvecs): 0.223 --------------------------------- Symfc end -------------------------------- Max drift of force constants: 0.00000000 (xx) 0.00000000 (xx) Permutation basis: 324/324 Permutation basis: 3366/3366 Construct permutation basis matrix. Finding block diagonal structure in projector. Using scipy connected_components. Rank of projector: 256 Number of blocks in projector: 262 Finding block diagonal structure in projector. Using scipy connected_components. Number of blocks in projector (Sum rule): 4 --- Eigsh_solver_block: 1 / 4 --- Block_size: 96 Use standard eigh solver. --- Eigsh_solver_block: 2 / 4 --- Block_size: 72 Use standard eigh solver. --- Eigsh_solver_block: 3 / 4 --- Block_size: 64 Use standard eigh solver. --- Eigsh_solver_block: 4 / 4 --- Block_size: 24 Use standard eigh solver. Tree of FC basis block matrices: - (256, 250), data: False |-- (24, 24), data: True |-- (64, 61), data: True |-- (72, 72), data: True |-- (96, 93), 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-08 08:42:38]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate fc3 ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 08:42:38]------------------------- 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 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P4/nmm (129) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.209387229999999 0.000000000000000 0.000000000000000 b 0.000000000000000 4.209387229999999 0.000000000000000 c 0.000000000000000 0.000000000000000 7.033745940000000 Atomic positions (fractional): 1 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 2 O 0.50000000000000 0.50000000000000 0.00000000000000 15.999 3 Cl 0.50000000000000 0.00000000000000 0.63129275280323 35.453 4 Cl 0.00000000000000 0.50000000000000 0.36870724719677 35.453 5 Ac 0.00000000000000 0.50000000000000 0.82042027633109 227.000 6 Ac 0.50000000000000 0.00000000000000 0.17957972366891 227.000 -------------------------------- supercell --------------------------------- Lattice vectors: a 12.628161689999997 0.000000000000000 0.000000000000000 b 0.000000000000000 12.628161689999997 0.000000000000000 c 0.000000000000000 0.000000000000000 14.067491880000000 Atomic positions (fractional): 1 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 > 1 2 O 0.33333333333333 0.00000000000000 0.00000000000000 15.999 > 1 3 O 0.66666666666667 0.00000000000000 0.00000000000000 15.999 > 1 4 O 0.00000000000000 0.33333333333333 0.00000000000000 15.999 > 1 5 O 0.33333333333333 0.33333333333333 0.00000000000000 15.999 > 1 6 O 0.66666666666667 0.33333333333333 0.00000000000000 15.999 > 1 7 O 0.00000000000000 0.66666666666667 0.00000000000000 15.999 > 1 8 O 0.33333333333333 0.66666666666667 0.00000000000000 15.999 > 1 9 O 0.66666666666667 0.66666666666667 0.00000000000000 15.999 > 1 10 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 1 11 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 1 12 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 1 13 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 1 14 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 1 15 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 1 16 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 1 17 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 1 18 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 1 19 O 0.16666666666667 0.16666666666667 0.00000000000000 15.999 > 19 20 O 0.50000000000000 0.16666666666667 0.00000000000000 15.999 > 19 21 O 0.83333333333333 0.16666666666667 0.00000000000000 15.999 > 19 22 O 0.16666666666667 0.50000000000000 0.00000000000000 15.999 > 19 23 O 0.50000000000000 0.50000000000000 0.00000000000000 15.999 > 19 24 O 0.83333333333333 0.50000000000000 0.00000000000000 15.999 > 19 25 O 0.16666666666667 0.83333333333333 0.00000000000000 15.999 > 19 26 O 0.50000000000000 0.83333333333333 0.00000000000000 15.999 > 19 27 O 0.83333333333333 0.83333333333333 0.00000000000000 15.999 > 19 28 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 19 29 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 19 30 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 19 31 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 19 32 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 19 33 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 19 34 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 19 35 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 19 36 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 19 37 Cl 0.16666666666667 0.00000000000000 0.31564637640161 35.453 > 37 38 Cl 0.50000000000000 0.00000000000000 0.31564637640161 35.453 > 37 39 Cl 0.83333333333333 0.00000000000000 0.31564637640161 35.453 > 37 40 Cl 0.16666666666667 0.33333333333333 0.31564637640161 35.453 > 37 41 Cl 0.50000000000000 0.33333333333333 0.31564637640161 35.453 > 37 42 Cl 0.83333333333333 0.33333333333333 0.31564637640161 35.453 > 37 43 Cl 0.16666666666667 0.66666666666667 0.31564637640161 35.453 > 37 44 Cl 0.50000000000000 0.66666666666667 0.31564637640161 35.453 > 37 45 Cl 0.83333333333333 0.66666666666667 0.31564637640161 35.453 > 37 46 Cl 0.16666666666667 0.00000000000000 0.81564637640161 35.453 > 37 47 Cl 0.50000000000000 0.00000000000000 0.81564637640161 35.453 > 37 48 Cl 0.83333333333333 0.00000000000000 0.81564637640161 35.453 > 37 49 Cl 0.16666666666667 0.33333333333333 0.81564637640161 35.453 > 37 50 Cl 0.50000000000000 0.33333333333333 0.81564637640161 35.453 > 37 51 Cl 0.83333333333333 0.33333333333333 0.81564637640161 35.453 > 37 52 Cl 0.16666666666667 0.66666666666667 0.81564637640161 35.453 > 37 53 Cl 0.50000000000000 0.66666666666667 0.81564637640161 35.453 > 37 54 Cl 0.83333333333333 0.66666666666667 0.81564637640161 35.453 > 37 55 Cl 0.00000000000000 0.16666666666667 0.18435362359839 35.453 > 55 56 Cl 0.33333333333333 0.16666666666667 0.18435362359839 35.453 > 55 57 Cl 0.66666666666667 0.16666666666667 0.18435362359839 35.453 > 55 58 Cl 0.00000000000000 0.50000000000000 0.18435362359839 35.453 > 55 59 Cl 0.33333333333333 0.50000000000000 0.18435362359839 35.453 > 55 60 Cl 0.66666666666667 0.50000000000000 0.18435362359839 35.453 > 55 61 Cl 0.00000000000000 0.83333333333333 0.18435362359839 35.453 > 55 62 Cl 0.33333333333333 0.83333333333333 0.18435362359839 35.453 > 55 63 Cl 0.66666666666667 0.83333333333333 0.18435362359839 35.453 > 55 64 Cl 0.00000000000000 0.16666666666667 0.68435362359839 35.453 > 55 65 Cl 0.33333333333333 0.16666666666667 0.68435362359839 35.453 > 55 66 Cl 0.66666666666667 0.16666666666667 0.68435362359839 35.453 > 55 67 Cl 0.00000000000000 0.50000000000000 0.68435362359839 35.453 > 55 68 Cl 0.33333333333333 0.50000000000000 0.68435362359839 35.453 > 55 69 Cl 0.66666666666667 0.50000000000000 0.68435362359839 35.453 > 55 70 Cl 0.00000000000000 0.83333333333333 0.68435362359839 35.453 > 55 71 Cl 0.33333333333333 0.83333333333333 0.68435362359839 35.453 > 55 72 Cl 0.66666666666667 0.83333333333333 0.68435362359839 35.453 > 55 73 Ac 0.00000000000000 0.16666666666667 0.41021013816555 227.000 > 73 74 Ac 0.33333333333333 0.16666666666667 0.41021013816555 227.000 > 73 75 Ac 0.66666666666667 0.16666666666667 0.41021013816555 227.000 > 73 76 Ac 0.00000000000000 0.50000000000000 0.41021013816555 227.000 > 73 77 Ac 0.33333333333333 0.50000000000000 0.41021013816555 227.000 > 73 78 Ac 0.66666666666667 0.50000000000000 0.41021013816555 227.000 > 73 79 Ac 0.00000000000000 0.83333333333333 0.41021013816555 227.000 > 73 80 Ac 0.33333333333333 0.83333333333333 0.41021013816555 227.000 > 73 81 Ac 0.66666666666667 0.83333333333333 0.41021013816555 227.000 > 73 82 Ac 0.00000000000000 0.16666666666667 0.91021013816555 227.000 > 73 83 Ac 0.33333333333333 0.16666666666667 0.91021013816555 227.000 > 73 84 Ac 0.66666666666667 0.16666666666667 0.91021013816555 227.000 > 73 85 Ac 0.00000000000000 0.50000000000000 0.91021013816555 227.000 > 73 86 Ac 0.33333333333333 0.50000000000000 0.91021013816555 227.000 > 73 87 Ac 0.66666666666667 0.50000000000000 0.91021013816555 227.000 > 73 88 Ac 0.00000000000000 0.83333333333333 0.91021013816555 227.000 > 73 89 Ac 0.33333333333333 0.83333333333333 0.91021013816555 227.000 > 73 90 Ac 0.66666666666667 0.83333333333333 0.91021013816555 227.000 > 73 91 Ac 0.16666666666667 0.00000000000000 0.08978986183445 227.000 > 91 92 Ac 0.50000000000000 0.00000000000000 0.08978986183445 227.000 > 91 93 Ac 0.83333333333333 0.00000000000000 0.08978986183445 227.000 > 91 94 Ac 0.16666666666667 0.33333333333333 0.08978986183445 227.000 > 91 95 Ac 0.50000000000000 0.33333333333333 0.08978986183445 227.000 > 91 96 Ac 0.83333333333333 0.33333333333333 0.08978986183445 227.000 > 91 97 Ac 0.16666666666667 0.66666666666667 0.08978986183445 227.000 > 91 98 Ac 0.50000000000000 0.66666666666667 0.08978986183445 227.000 > 91 99 Ac 0.83333333333333 0.66666666666667 0.08978986183445 227.000 > 91 100 Ac 0.16666666666667 0.00000000000000 0.58978986183445 227.000 > 91 101 Ac 0.50000000000000 0.00000000000000 0.58978986183445 227.000 > 91 102 Ac 0.83333333333333 0.00000000000000 0.58978986183445 227.000 > 91 103 Ac 0.16666666666667 0.33333333333333 0.58978986183445 227.000 > 91 104 Ac 0.50000000000000 0.33333333333333 0.58978986183445 227.000 > 91 105 Ac 0.83333333333333 0.33333333333333 0.58978986183445 227.000 > 91 106 Ac 0.16666666666667 0.66666666666667 0.58978986183445 227.000 > 91 107 Ac 0.50000000000000 0.66666666666667 0.58978986183445 227.000 > 91 108 Ac 0.83333333333333 0.66666666666667 0.58978986183445 227.000 > 91 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py_mlp_eval_fc3_disp.yaml.xz". --------------------------- Dielectric constant ---------------------------- 4.2012710 0.0000000 0.0000000 0.0000000 4.2012710 0.0000000 0.0000000 0.0000000 4.0787131 -------------------------- Born effective charges -------------------------- 1 O -2.8552025 0.0000000 0.0000000 0.0000000 -2.8552025 0.0000000 0.0000000 0.0000000 -2.2085931 2 O -2.8552025 0.0000000 0.0000000 0.0000000 -2.8552025 0.0000000 0.0000000 0.0000000 -2.2085931 3 Cl -1.3775560 0.0000000 0.0000000 0.0000000 -1.3775560 0.0000000 0.0000000 0.0000000 -1.4905265 4 Cl -1.3775560 0.0000000 0.0000000 0.0000000 -1.3775560 0.0000000 0.0000000 0.0000000 -1.4905265 5 Ac 4.2327586 0.0000000 0.0000000 0.0000000 4.2327586 0.0000000 0.0000000 0.0000000 3.6991195 6 Ac 4.2327586 0.0000000 0.0000000 0.0000000 4.2327586 0.0000000 0.0000000 0.0000000 3.6991195 ---------------------------------------------------------------------------- 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[ 37, 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[ 73, 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 (zzz) 0.00000000 (zzz) 0.00000000 (zzz) fc3 was written into "fc3.hdf5". Max drift of fc2: -0.00000000 (xx) -0.00000000 (xx) fc2 was written into "fc2.hdf5". --------------------------- Calculation settings --------------------------- Non-analytical term correction (NAC): True NAC unit conversion factor: 14.39965 BZ integration: Tetrahedron-method Temperatures: 0.0 300.0 Cutoff frequency: 0.01 Frequency conversion factor to THz: 15.63330 ----------- None of ph-ph interaction calculation was performed. ----------- Summary of calculation was written in "phono3py.yaml". -------------------------[time 2026-01-08 08:42:42]------------------------- _ ___ _ __ __| | / _ \ '_ \ / _` | | __/ | | | (_| | \___|_| |_|\__,_| ---------------------------- ------ calculate LTC ------- ---------------------------- _ _____ _ __ | |__ ___ _ __ ___|___ / _ __ _ _ | '_ \| '_ \ / _ \| '_ \ / _ \ |_ \| '_ \| | | | | |_) | | | | (_) | | | | (_) |__) | |_) | |_| | | .__/|_| |_|\___/|_| |_|\___/____/| .__/ \__, | |_| |_| |___/ 3.23.0 -------------------------[time 2026-01-08 08:42:43]------------------------- 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 2] Primitive matrix: [1. 0. 0.] [0. 1. 0.] [0. 0. 1.] Spacegroup: P4/nmm (129) ------------------------------ primitive cell ------------------------------ Lattice vectors: a 4.209387229999999 0.000000000000000 0.000000000000000 b 0.000000000000000 4.209387229999999 0.000000000000000 c 0.000000000000000 0.000000000000000 7.033745940000000 Atomic positions (fractional): 1 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 2 O 0.50000000000000 0.50000000000000 0.00000000000000 15.999 3 Cl 0.50000000000000 0.00000000000000 0.63129275280323 35.453 4 Cl 0.00000000000000 0.50000000000000 0.36870724719677 35.453 5 Ac 0.00000000000000 0.50000000000000 0.82042027633109 227.000 6 Ac 0.50000000000000 0.00000000000000 0.17957972366891 227.000 -------------------------------- supercell --------------------------------- Lattice vectors: a 12.628161689999997 0.000000000000000 0.000000000000000 b 0.000000000000000 12.628161689999997 0.000000000000000 c 0.000000000000000 0.000000000000000 14.067491880000000 Atomic positions (fractional): 1 O 0.00000000000000 0.00000000000000 0.00000000000000 15.999 > 1 2 O 0.33333333333333 0.00000000000000 0.00000000000000 15.999 > 1 3 O 0.66666666666667 0.00000000000000 0.00000000000000 15.999 > 1 4 O 0.00000000000000 0.33333333333333 0.00000000000000 15.999 > 1 5 O 0.33333333333333 0.33333333333333 0.00000000000000 15.999 > 1 6 O 0.66666666666667 0.33333333333333 0.00000000000000 15.999 > 1 7 O 0.00000000000000 0.66666666666667 0.00000000000000 15.999 > 1 8 O 0.33333333333333 0.66666666666667 0.00000000000000 15.999 > 1 9 O 0.66666666666667 0.66666666666667 0.00000000000000 15.999 > 1 10 O 0.00000000000000 0.00000000000000 0.50000000000000 15.999 > 1 11 O 0.33333333333333 0.00000000000000 0.50000000000000 15.999 > 1 12 O 0.66666666666667 0.00000000000000 0.50000000000000 15.999 > 1 13 O 0.00000000000000 0.33333333333333 0.50000000000000 15.999 > 1 14 O 0.33333333333333 0.33333333333333 0.50000000000000 15.999 > 1 15 O 0.66666666666667 0.33333333333333 0.50000000000000 15.999 > 1 16 O 0.00000000000000 0.66666666666667 0.50000000000000 15.999 > 1 17 O 0.33333333333333 0.66666666666667 0.50000000000000 15.999 > 1 18 O 0.66666666666667 0.66666666666667 0.50000000000000 15.999 > 1 19 O 0.16666666666667 0.16666666666667 0.00000000000000 15.999 > 19 20 O 0.50000000000000 0.16666666666667 0.00000000000000 15.999 > 19 21 O 0.83333333333333 0.16666666666667 0.00000000000000 15.999 > 19 22 O 0.16666666666667 0.50000000000000 0.00000000000000 15.999 > 19 23 O 0.50000000000000 0.50000000000000 0.00000000000000 15.999 > 19 24 O 0.83333333333333 0.50000000000000 0.00000000000000 15.999 > 19 25 O 0.16666666666667 0.83333333333333 0.00000000000000 15.999 > 19 26 O 0.50000000000000 0.83333333333333 0.00000000000000 15.999 > 19 27 O 0.83333333333333 0.83333333333333 0.00000000000000 15.999 > 19 28 O 0.16666666666667 0.16666666666667 0.50000000000000 15.999 > 19 29 O 0.50000000000000 0.16666666666667 0.50000000000000 15.999 > 19 30 O 0.83333333333333 0.16666666666667 0.50000000000000 15.999 > 19 31 O 0.16666666666667 0.50000000000000 0.50000000000000 15.999 > 19 32 O 0.50000000000000 0.50000000000000 0.50000000000000 15.999 > 19 33 O 0.83333333333333 0.50000000000000 0.50000000000000 15.999 > 19 34 O 0.16666666666667 0.83333333333333 0.50000000000000 15.999 > 19 35 O 0.50000000000000 0.83333333333333 0.50000000000000 15.999 > 19 36 O 0.83333333333333 0.83333333333333 0.50000000000000 15.999 > 19 37 Cl 0.16666666666667 0.00000000000000 0.31564637640161 35.453 > 37 38 Cl 0.50000000000000 0.00000000000000 0.31564637640161 35.453 > 37 39 Cl 0.83333333333333 0.00000000000000 0.31564637640161 35.453 > 37 40 Cl 0.16666666666667 0.33333333333333 0.31564637640161 35.453 > 37 41 Cl 0.50000000000000 0.33333333333333 0.31564637640161 35.453 > 37 42 Cl 0.83333333333333 0.33333333333333 0.31564637640161 35.453 > 37 43 Cl 0.16666666666667 0.66666666666667 0.31564637640161 35.453 > 37 44 Cl 0.50000000000000 0.66666666666667 0.31564637640161 35.453 > 37 45 Cl 0.83333333333333 0.66666666666667 0.31564637640161 35.453 > 37 46 Cl 0.16666666666667 0.00000000000000 0.81564637640161 35.453 > 37 47 Cl 0.50000000000000 0.00000000000000 0.81564637640161 35.453 > 37 48 Cl 0.83333333333333 0.00000000000000 0.81564637640161 35.453 > 37 49 Cl 0.16666666666667 0.33333333333333 0.81564637640161 35.453 > 37 50 Cl 0.50000000000000 0.33333333333333 0.81564637640161 35.453 > 37 51 Cl 0.83333333333333 0.33333333333333 0.81564637640161 35.453 > 37 52 Cl 0.16666666666667 0.66666666666667 0.81564637640161 35.453 > 37 53 Cl 0.50000000000000 0.66666666666667 0.81564637640161 35.453 > 37 54 Cl 0.83333333333333 0.66666666666667 0.81564637640161 35.453 > 37 55 Cl 0.00000000000000 0.16666666666667 0.18435362359839 35.453 > 55 56 Cl 0.33333333333333 0.16666666666667 0.18435362359839 35.453 > 55 57 Cl 0.66666666666667 0.16666666666667 0.18435362359839 35.453 > 55 58 Cl 0.00000000000000 0.50000000000000 0.18435362359839 35.453 > 55 59 Cl 0.33333333333333 0.50000000000000 0.18435362359839 35.453 > 55 60 Cl 0.66666666666667 0.50000000000000 0.18435362359839 35.453 > 55 61 Cl 0.00000000000000 0.83333333333333 0.18435362359839 35.453 > 55 62 Cl 0.33333333333333 0.83333333333333 0.18435362359839 35.453 > 55 63 Cl 0.66666666666667 0.83333333333333 0.18435362359839 35.453 > 55 64 Cl 0.00000000000000 0.16666666666667 0.68435362359839 35.453 > 55 65 Cl 0.33333333333333 0.16666666666667 0.68435362359839 35.453 > 55 66 Cl 0.66666666666667 0.16666666666667 0.68435362359839 35.453 > 55 67 Cl 0.00000000000000 0.50000000000000 0.68435362359839 35.453 > 55 68 Cl 0.33333333333333 0.50000000000000 0.68435362359839 35.453 > 55 69 Cl 0.66666666666667 0.50000000000000 0.68435362359839 35.453 > 55 70 Cl 0.00000000000000 0.83333333333333 0.68435362359839 35.453 > 55 71 Cl 0.33333333333333 0.83333333333333 0.68435362359839 35.453 > 55 72 Cl 0.66666666666667 0.83333333333333 0.68435362359839 35.453 > 55 73 Ac 0.00000000000000 0.16666666666667 0.41021013816555 227.000 > 73 74 Ac 0.33333333333333 0.16666666666667 0.41021013816555 227.000 > 73 75 Ac 0.66666666666667 0.16666666666667 0.41021013816555 227.000 > 73 76 Ac 0.00000000000000 0.50000000000000 0.41021013816555 227.000 > 73 77 Ac 0.33333333333333 0.50000000000000 0.41021013816555 227.000 > 73 78 Ac 0.66666666666667 0.50000000000000 0.41021013816555 227.000 > 73 79 Ac 0.00000000000000 0.83333333333333 0.41021013816555 227.000 > 73 80 Ac 0.33333333333333 0.83333333333333 0.41021013816555 227.000 > 73 81 Ac 0.66666666666667 0.83333333333333 0.41021013816555 227.000 > 73 82 Ac 0.00000000000000 0.16666666666667 0.91021013816555 227.000 > 73 83 Ac 0.33333333333333 0.16666666666667 0.91021013816555 227.000 > 73 84 Ac 0.66666666666667 0.16666666666667 0.91021013816555 227.000 > 73 85 Ac 0.00000000000000 0.50000000000000 0.91021013816555 227.000 > 73 86 Ac 0.33333333333333 0.50000000000000 0.91021013816555 227.000 > 73 87 Ac 0.66666666666667 0.50000000000000 0.91021013816555 227.000 > 73 88 Ac 0.00000000000000 0.83333333333333 0.91021013816555 227.000 > 73 89 Ac 0.33333333333333 0.83333333333333 0.91021013816555 227.000 > 73 90 Ac 0.66666666666667 0.83333333333333 0.91021013816555 227.000 > 73 91 Ac 0.16666666666667 0.00000000000000 0.08978986183445 227.000 > 91 92 Ac 0.50000000000000 0.00000000000000 0.08978986183445 227.000 > 91 93 Ac 0.83333333333333 0.00000000000000 0.08978986183445 227.000 > 91 94 Ac 0.16666666666667 0.33333333333333 0.08978986183445 227.000 > 91 95 Ac 0.50000000000000 0.33333333333333 0.08978986183445 227.000 > 91 96 Ac 0.83333333333333 0.33333333333333 0.08978986183445 227.000 > 91 97 Ac 0.16666666666667 0.66666666666667 0.08978986183445 227.000 > 91 98 Ac 0.50000000000000 0.66666666666667 0.08978986183445 227.000 > 91 99 Ac 0.83333333333333 0.66666666666667 0.08978986183445 227.000 > 91 100 Ac 0.16666666666667 0.00000000000000 0.58978986183445 227.000 > 91 101 Ac 0.50000000000000 0.00000000000000 0.58978986183445 227.000 > 91 102 Ac 0.83333333333333 0.00000000000000 0.58978986183445 227.000 > 91 103 Ac 0.16666666666667 0.33333333333333 0.58978986183445 227.000 > 91 104 Ac 0.50000000000000 0.33333333333333 0.58978986183445 227.000 > 91 105 Ac 0.83333333333333 0.33333333333333 0.58978986183445 227.000 > 91 106 Ac 0.16666666666667 0.66666666666667 0.58978986183445 227.000 > 91 107 Ac 0.50000000000000 0.66666666666667 0.58978986183445 227.000 > 91 108 Ac 0.83333333333333 0.66666666666667 0.58978986183445 227.000 > 91 ---------------------------------------------------------------------------- NAC parameters were read from "phono3py.yaml". --------------------------- Dielectric constant ---------------------------- 4.2012710 0.0000000 0.0000000 0.0000000 4.2012710 0.0000000 0.0000000 0.0000000 4.0787131 -------------------------- Born effective charges -------------------------- 1 O -2.8552025 0.0000000 0.0000000 0.0000000 -2.8552025 0.0000000 0.0000000 0.0000000 -2.2085931 2 O -2.8552025 0.0000000 0.0000000 0.0000000 -2.8552025 0.0000000 0.0000000 0.0000000 -2.2085931 3 Cl -1.3775560 0.0000000 0.0000000 0.0000000 -1.3775560 0.0000000 0.0000000 0.0000000 -1.4905265 4 Cl -1.3775560 0.0000000 0.0000000 0.0000000 -1.3775560 0.0000000 0.0000000 0.0000000 -1.4905265 5 Ac 4.2327586 0.0000000 0.0000000 0.0000000 4.2327586 0.0000000 0.0000000 0.0000000 3.6991195 6 Ac 4.2327586 0.0000000 0.0000000 0.0000000 4.2327586 0.0000000 0.0000000 0.0000000 3.6991195 ---------------------------------------------------------------------------- fc3 was read from "fc3.hdf5". fc2 was read from "fc2.hdf5". ----------------------------- Force constants ------------------------------ Max drift of fc3: 0.00000000 (zzz) 0.00000000 (zzz) 0.00000000 (zzz) 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 ... [ 12 12 7 ] 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.83, Number of G-points: 303, Lambda: 0.18 Running harmonic phonon calculations... -------------------- Lattice thermal conductivity (RTA) --------------------