This thesis is a theoretical study of configurationally disordered icosahedral boronrich solids, in particular boron carbides, using density functional theory and alloy theory. The goal is to resolve discrepancies, regarding the properties of boron carbides, between experiments and previous theoretical calculations which have been a controversial issue in the field of icosahedral boron-rich solids. For instance, B 13 C 2 is observed experimentally to be a semiconductor, meanwhile electronic band structure calculations reveal a metallic character of B 13 C 2 due to its electron deficiency. In B 4 C, on the other hand, the experimentally observed band gap is unexpectedly smaller, not the usual larger, than that of standard DFT calculations. Another example is given by the existence of a small structural distortion in B 4 C, as predicted in theoretical calculations, which reduces the crystal symmetry from the experimentally observed rhombohedral (R3m) to the based-centered monoclinic (Cm). Since boron carbide is stable as a single-phase over a broad composition range (∼8-20 at.% C), substitution of boron and carbon atoms for one another is conceivable. For this reason, the discrepancies have been speculated in the literature, without a proof, to originate from configurational disorder induced by substitutional defects. However, owing to its complex atomic structure, represented by 12-atom icosahedra and 3-atom intericosahedral chains, a practical alloy theory method for direct calculations of the properties of the relevant configurations of disordered boron carbides, as well as for a thermodynamic assessment of their stability has been missing.In this thesis, a new approach, the superatom-special quasirandom structure (SA-SQS), has been developed. The approach allows one to model configurational disorder in boron carbide, induced by high concentrations of low-energy B/C substitutional defects. B 13 C 2 and B 4 C are the two stoichiometries, mainly considered in this study, as they are of particular importance and have been in focus in the literature. The results demonstrate that, from thermodynamic considerations, both B 13 C 2 and B 4 C configurationally disorder at high temperature. In the case of B 13 C 2 , the configurational disorder splits off some valence states into the band gap that in turn compensates the electron deficiency in ordered B 13 C 2 , thus resulting in a semiconducting character. As for B 4 C, the configurational disorder eliminates the monoclinic distortion, thus resulting in the restoration of the higher rhombohedral symmetry. Configurational disorder can also account for an exceliii iv lent agreement on elastic moduli of boron carbide between theory and experiment. Thus, several of the previous discrepancies between theory and experiments are resolved.Inspired by attempts to enhance the mechanical properties of boron suboxide by fabricating boron suboxide-boron carbide composites, as recently suggested in the literature, the SA-SQS approach is used for modeling mixtures of boron sub...