“…The simulation box comprised 12 × 12 × 12 unit cells of Zn 4 Sb 3 containing 14 904 atoms, with periodic boundary conditions applied in all three orthogonal directions. The temperature ranged from 300 to 700 K in intervals of 25 K. The Nosé–Hoover thermostat was used to control the temperature, enabling the system to reach thermal equilibrium at different temperatures before performing any statistical analysis. , The atomic interaction was described by a three-body potential that was used to study the thermoelectric, structural, and mechanical properties of the Zn 4 Sb 3 crystals. ,, Energy minimization was performed using the conjugate gradient method to relax the initial structure, which was built to correspond to the measurements of crystal structure for the samples used in our experimental study. The velocity Verlet algorithm was used for time integration, with a step of 1.0 fs in an isothermal–isobaric (NPT) ensemble for 10 ns, to allow the system to reach thermal equilibrium.…”