Equilibrium bond lengths, harmonic frequencies, electron affinities, and dissociation energies (as appropriate) were determined for I, I-,I-2, I-2(-), and I-3(-) by CCSD(T) and four DFT methods (BLYP, BPW91, B3LYP, and B3PW91) using a basis set consisting of a relativistic effective core potential and a triple-zeta plus 2df polarization functions (ECP-TZ(2df)) for the valence electrons. Comparison of the DFT results with the CCSD(T) results and available experimental information indicates that the B3PW91 approach does particularly well describing the bonding in these species. B3PW91/ECP-TZ(2df) calculations of I-5(-) in linear (D-infinity h) and bent (C-2v) geometries indicate that the linear structure is a low-energy transition state lying only about 0.1-0.2 eV above the C-2v global minimum-energy structure. Harmonic frequencies and infrared and Raman intensities are calculated for both structures and compared to available experimental information.