The potential energy surfaces of a series of clusters with the formula CAl₃E (E = P, As, Sb, Bi) are systematically explored using density functional theory and high level <italic>ab initio</italic> calculations.