A systematic investigation of rhodium-doped silicon clusters, Rh 2 Si n q with n = 1-10 and q = 0, ±1, in the neutral, anionic, and cationic states is performed using density functional theory approach at B3LYP/GENECP level. According to the optimum Rh 2 Si n q clusters, mostly equilibrium geometries prefer the three-dimensional structures for n = 2-10. When n = 10, one Rh atom in Rh 2 Si 10 0,±1 clusters completely falls into the center of Si frame, and cage-like Rh 2 Si 10 0,±1 geometries are formed. The Rh 2 Si 1,6-9 ? and Rh 2 Si 5,7,9 -clusters significantly deform their corresponding neutral geometries, which are in line with the calculated ionization potential and electron affinity values. The relative stabilities of Rh 2 Si n q clusters for the lowest-energy structures are analyzed on the basis of binding energy, fragmentation energy, second-order energy difference, and HOMO-LUMO gaps. The theoretical results confirm that the Rh 2 Si 6 -, Rh 2 Si 6 , and Rh 2 Si 6 ? clusters are more stable than their neighboring ones. The natural population analysis reveals that the charges in Rh 2 Si n q clusters transfer from the Si atoms to the Rh atoms except Rh 2 Si ? . In addition, the relationship between static polarizability and HOMO-LUMO gaps is discussed.