Here we present an automated workflow designed to (i) provide a bridge between atomic quantum simulation at zero temperature and macroscopic thermal properties at finite temperature, (ii) help realize large-scale automated anharmonic phonon calculations, and (iii) serve as a user-friendly pipeline with benchmarked input parameters. This comprehensive lattice dynamic (LD) workflow calculates beyond harmonic phonon properties of interest, including lattice thermal conductivity (LTC), coefficient of thermal expansion (CTE), finite-temperature phonons, and vibrational free energy at finite temperatures. It uses a suite of software packages and is optimized with suitable input parameters (supercell size, cutoff radius for fitting IFCs, etc) that balance accuracy and computational costs. Some immediate potential applications of the results are for (i) high-throughput calculations of thermoelectric figure of merit (zT ), using combined anharmonic phonons and electron scattering properties obtained from the recent development of AMSET , (ii) phase transitions of metastable materials at finite temperatures, and (iii) searching for negative thermal expansion (NTE) materials.