Whereas low-temperature ferroelectrics have a well understood ordered spatial dipole arrangement, the fate of these dipoles in paraelectric phases remains poorly understood. This is studied here as an energy minimization problem using both static and molecular dynamic (MD) density functional theory (DFT). We find that considering the non-thermal internal energy already reveals the formation of a distribution of static local displacements that (i) mimic the symmetries of the low temperature phases, while (ii) being the precursors of what high temperature DFT MD finds as thermal motifs.