Analysis of accurate experimental and theoretical structure factors of diamond and silicon reveals that the contraction of the core shell due to covalent bond formation causes significant perturbations of the total charge density that cannot be ignored in precise charge density studies. We outline that the nature and origin of core contraction/expansion and core polarization phenomena can be analyzed by experimental studies employing an extended Hansen-Coppens multipolar model. Omission or insufficient treatment of these subatomic charge density phenomena might yield erroneous thermal displacement parameters and high residual densities in multipolar refinements. Our detailed studies therefore suggest that the refinement of contraction/expansion and population parameters of all atomic shells is essential to the precise reconstruction of electron density distributions by a multipolar model. Furthermore, our results imply that also the polarization of the inner shells needs to be adopted, especially in cases where second row or even heavier elements are involved in covalent bonding. These theoretical studies are supported by direct multipolar refinements of X-ray powder diffraction data of diamond obtained from a third-generation synchrotron-radiation source (SPring-8, BL02B2).