Second-order IFCs, which define phonons, can be calculated through DFPT, finite- Our anharmonic phonon workflow 106 automatically calculates interatomic force constants up to 4th order from perturbed training supercells, and uses them to obtain lattice thermal conductivity, coefficient of thermal expansion, and vibrational free energy and entropy. The workflow starts with a primitive structure and adjustable parameters such as the force field