We investigate the thermodynamic stability of carbon-rich icosahedral boron carbide at different compositions, ranging from B4C to B2C, using first-principles calculations. Apart from B4C, generally addressed in the literature, B2.5C, represented by B10C2p(C-C), where Cp and (C-C) denote a carbon atom occupying the polar site of the icosahedral cluster and a diatomic carbon chain, respectively, is predicted to be thermodynamically stable under high pressures with respect to B4C as well as pure boron and carbon phases. The thermodynamic stability of B2.5C is determined by the Gibbs free energy G as a function of pressure p and temperature T, in which the contributions from the lattice vibrations and the configurational disorder are obtained ...