We still refer to these as "phases" everywhere in the code.
Working with phase factors means that users do not have to
exponentiate anything (and potentially get an incorrect sign).
In addition this is closer to what we need for non-abelian
gauge fields.