Exactly solvable model of a superconducting to rotational phase transition

We consider a many-fermion model which exhibits a transition from a superconducting to a rotational phase with variation of a parameter in its Hamiltonian. The model has analytical solutions in its two limits due to the presence of dynamical symmetries. However, the symmetries are basically incompatible with one another; no simple solution exists in intermediate situations. Exact (numerical) solutions are possible and enable one to study the behavior of competing but incompatible symmetries and the phase transitions that result in a semirealistic situation. The results are remarkably simple and shed light on the nature of phase transitions.