Total energies of crystal structures can be calculated to high precision using quantumbased density functional theory (DFT) methods, but the calculations can be time consuming and scale badly with system size. Boron carbide exhibits disorder in the distribution of boron and carbon atoms among the crystallographic sites. A cluster expansion of the DFT energy in a series of pairs, triplets, etc. ...