The distribution of oxygen atoms around their equilibrium positions from AIMD at 2500 K is shown in Fig. . Here, no signs of ion o↵-centering is observed. We do observe strong local correlations corresponding to the displacement of ions along the soft-mode, which exist along all three cartesian directions, corresponding to the three possible symmetry lowering pathways from cubic to tetragonal. However, this is exactly the expected behaviour from a low energy phonon mode. Thus, there are no signs of the cubic phase being an average over multiple locally tetragonal phases. While there is no evidence to support their formation, we cannot rule out that such domains do emerge over longer length-and/or time scales.