We note that the overall workflow requires the selection of a multitude of input parameters. Such parameters are typically adjusted manually, which is infeasible for high-throughput calculations. Here, we provide a concise summary of the most critical parameters of our workflow with benchmarking, with the goal of enabling automatic selection of optimal inputs and ensuring consistency across large-scale calculations. Functional: It is essential to obtain precise lattice parameters for accurate phonon calculations. In the workflow, all DFT calculations are performed using Perdew-Burke-Ernzerhof revised for solids (PBEsol) exchangecorrelation functional and projector-augmented wave pseudopotentials (PAW) . PBEsol is selected over PBE because the latter tends to overestimate lattice parameters. Furthermore, from the recent work by Ning et al., the PBE functional is shown to underestimate the phonon frequencies in their testing cases of GaAs, Fe, NiO. The r 2 SCAN functional is also a viable option, but it is not set as the default due to computational cost. According to the extensive benchmarking by Kingsbury et al. , meta-GGA (generalized gradient approximation) functionals such as r 2 SCAN are usually 2-4 times more expensive than GGA functionals such as PBE and PBEsol.