Ligand chemical shifts are calculated and analyzed for three paramagnetic transition metal tris-acetylacetonato (acac) complexes, namely high-spin Fe(III) and Cr(III), and low-spin Ru(III), using scalar relativistic density functional theory (DFT). The signs and magnitudes of the paramagnetic NMR ligand chemical shifts are directly related to the extent of covalent acac oxygen-to-metal sigma donation involving unoccupied metal valence d(sigma) acceptor orbitals. The role of delocalization of metal-centered spin density over the ligand atoms plays a minor secondary role. Of particular interest is the origin of the sign and magnitude of the methyl carbon chemical shift in the acac ligands, and the role played by the DFT delocalization error when calculating such shifts. It is found that the alpha versus beta spin balance of oxygen sigma donation to metal valence d acceptor orbitals is responsible for the sign and the magnitude of the ligand methyl carbon chemical shift. A problematic case is the methyl carbon shift of Fe(acac)(3). Most functionals produce shifts in excess of 1400 ppm, whereas the experimental shift is approximately 279 ppm. Range-separated hybrid functionals that are optimally tuned for Fe(acac)(3) based on DFT energetic criteria predict a lower limit of about 2000 ppm for the methyl carbon shift of the high-spin electronic configuration. Since the experimental value is based on a very strongly broadened signal it is possibly unreliable.