“…Instead, other sampling methods such as fitting forces from randomly displaced supercells or molecular dynamics simulations can be used to obtain higher order force constants. Methods to study the phonon anharmonicity and finite-temperature thermodynamics include the temperature-dependent effective potential (TDEP) method, phonon gas model (PGM) with the quasiparticle frequency approach, , iterative self-consistent phonon calculations (SCPH), − and related methods such as self-consistent ab initio lattice dynamics (SCAILD) , and stochastic self-consistent harmonic approximation (SSCHA). , Titanate perovskites BaTiO 3 and SrTiO 3 have been investigated with SCPH and SCPH-based methods, yielding information on phase transitions and finite-temperature phonon dispersions. ,− The anharmonic lattice dynamics of niobate perovskites KNbO 3 and NaNbO 3 have also been investigated using SCPH methods …”