Although the rutile structure of TiO$_2$ is stable at high temperatures, the
conventional quasiharmonic approximation predicts that several acoustic phonons
decrease anomalously to zero frequency with thermal expansion, incorrectly
predicting a structural collapse at temperatures well below 1000\,K. Inelastic
neutron scattering was used to measure the temperature dependence of the phonon
density of states (DOS) of rutile TiO$_2$ from 300 to 1373\,K. Surprisingly,
these anomalous acoustic phonons were found to increase in frequency with
temperature. First-principles calculations showed that with lattice expansion,
the potentials for the anomalous acoustic phonons transform from quadratic to
quartic, stabilizing the rutile phase at high temperatures. In these modes, the
vibrational displacements of adjacent Ti and O atoms cause variations in
hybridization of $3d$ electrons of Ti and $2p$ electrons of O atoms. With
thermal expansion, the energy variation in this "phonon-tracked hybridization"
flattens the bottom of the interatomic potential well between Ti and O atoms,
and induces a quarticity in the phonon potential.