We report on the size and charge dependence of the C-O stretching frequency, ͑CO͒, in complexes of CO with gas phase anionic, neutral, and cationic cobalt clusters ͑Co n CO −/0/+ ͒, anionic, neutral, and cationic rhodium clusters ͑Rh n CO −/0/+ ͒, and cationic nickel clusters ͑Ni n CO + ͒ for n up to 37. We develop models, based on the established vibrational spectroscopy of organometallic carbonyl compounds, to understand how cluster size and charge relate to ͑CO͒ in these complexes. The dominating factor is the available electron density for backdonation from the metal to the CO * orbital. Electrostatic effects play a significant but minor role. For the charged clusters, the size trends are related to the dilution of the charge density at the binding site on the cluster as n increases. At large n, ͑CO͒ approaches asymptotes that are not the same as found for ͑CO͒ on the single crystal metal surfaces, reflecting differences between binding sites on medium sized clusters and the more highly coordinated metal surface sites.