Nematic order in a degenerate Hubbard model with spin–orbit coupling