Ferromagnetic phase transitions in $SU(N)$

We study the thermodynamics of a non-abelian ferromagnet consisting of "atoms" each carrying a fundamental representation of $SU(N)$, coupled with long-range two-body quadratic interactions. We uncover a rich structure of phase transitions from non-magnetized, global $SU(N)$-invariant states to magnetized ones breaking global invariance to $SU(N-1) \times U(1)$. Phases can coexist, one being stable and the other metastable, and the transition between states involves latent heat exchange, unlike in usual $SU(2)$ ferromagnets. Coupling the system to an external non-abelian magnetic field further enriches the phase structure, leading to additional phases. The system manifests hysteresis phenomena both in the magnetic field, as in usual ferromagnets, and in the temperature, in analogy to supercooled water. Potential applications are in fundamental situations or as a phenomenological model.