In this article, we assess the ability of various density functionals to predict accurate values for some basic properties of the bond critical points of about SO small molecules, including the recently proposed reduced gradient variation rates and involving typical ionic and covalent bonds, agostic interactions, and van der Waals complexes. The relation between the computed deviations and the geometric variations are discussed, as well as the topology variations. The possible correlation of these descriptors to atomization energies is considered, and the relevance of an accurate QTAIM analysis for correct descriptions of potential energy surfaces is addressed. Finally, we provide typical margins of error for the evaluation of these quantities and discuss their consequences for computational applications.