## Add `positions` similar to `momenta` to `Model`

It would be useful for Hamiltonians with real space dependence (such as hopping terms in non-crystalline systems, or the continuum description of a quantum well) if we allowed definition of symbols that are `positions`

. The only difference compared to `momenta`

would be that they don't flip sign under antiunitary symmetry operations.

It could also be useful to allow multiple vectors that transform the same way, for example to treat interaction terms that have multiple momentum variables.

This is all possible in current `qsymm`

in a hacky way: extend the real space dimension to the total number of transforming symbols, use `momenta`

that is a single vector with all these vectors stacked, and make `R`

in `PointGroupElements`

block-diagonal with the appropriate transformation matrices in each block. The `positions`

type transformation can be faked by using `-R`

in antiunitary operators.

`<TLDR>`

We could consider supporting other similar quantities, like `angular_momenta`

which change sign under TR but do not change under spatial inversion (this really only makes sense in 3D, in lower dimensions it is ambiguous whether something is a rotation or an inversion). This would mostly be useful for classical dynamical variables, *not* for fixed external fields. For example one may ask, what is the symmetry of the electromagnetic term `E⋅B`

, a product of a `position`

and an `angular_momentum`

type variable (it is invariant under anything that doesn't change space-time orientation). However, when examining the symmetries of the quantum Hamiltonian `B⋅σ`

with fixed external field `B`

, `B`

should *not* be transformed. We should give this more thought, it is not clear to me how many types of such objects we need in general, whether the behaviour under PH and TR is independent.