From a high-resolution X-ray data set (sin theta/lambda = 1.1 A(-1)) measured at 20 K the electron-density distribution of the nucleoside thymidine was derived by a classical multipole refinement and by application of the invariom formalism. Owing to the presence of the heteroaromatic thymine ring system two invariom models were compared which considered the nearest and next-nearest neighbors for the invariom assignments. Differences between the two invariom models were small for the bond topological and atomic properties - about five times smaller than differences with the classical multipole refinement. Even the latter differences are in the uncertainty ranges which are commonly observed in experimental charge-density work and were found in molecular regions involved in intermolecular contacts. The application of the constrained wavefunction-fitting approach allowed the electron localization function (ELF) to be obtained from the experimental X-ray data, which was graphically represented and topologically analyzed. ELF basin populations were derived from experiment for the first time. The electron populations in the disynaptic valence basins were related quantitatively to bond orders.