The circular dichroism spectra of the tris-bidentate metal complexes Lambda-[M(phen)3]2+, with M = Fe, Ru, Os and phen = 1,10-tris-phenanthroline, are investigated computationally, employing time-dependent density functional theory. Good agreement with experimental spectra is obtained for Ru and Os. The Lambda-[Os(phen)3]2+ spectrum is analyzed in detail. It is shown how relativistic effects red shift CD bands where the Os 5d-orbital participates to a large extent in the excitations. Further, the participation of the metal in the ligand pi --> pi exciton CD is determined to be of the order of 10%. Though solvent effects can have a noticeable effect on individual transitions and rotatory strengths, they are demonstrated to have only a very small overall effect on the resulting simulated CD spectra. For Lambda-[Fe(phen)3]2+, the results are shown to be rather sensitive to the choice of the applied hybrid and nonhybrid density functionals, and the optimized geometries based thereupon. In particular, the sign pattern of the lower-energy part (up to 33 x 10(3) cm(-1)) of the Lambda-[Fe(phen)3]2+ CD spectrum is difficult to reproduce. Some combinations of functionals and geometries yield good agreement with experiment, but no "best" approach can be devised based on the available results. Possible sources of errors in the spectrum of Lambda-[Fe(phen)3]2+ due to deficiencies in the functionals and the exchange-correlation kernels are investigated.