“…The first-order self-consistent phonon (SC1) theory is one of the most successful methods, which determines the renormalized phonon frequencies by the variational principle applied to the first-order cumulant expansion of the Helmholtz free energy [16,17]. Since the SC1 theory can, to a large extent, remedy the negative frequency problems of the harmonic approximation, it has been actively employed in first-principles calculations of phononrelated physics of anharmonic materials, including thermal transport [13,[18][19][20][21], phonon-limited mobility [22], bandgap renormalization [23,24], thermal expansion [25,26], and conventional superconductivity [7,12,27]. While these studies clearly demonstrate the advantage of the SC1 theory over the quasiharmonic theory and purely perturbative approaches, the quantitative accuracy of SC1 is still inadequate for strongly anharmonic materials.…”