Multiple_HS_fields_per_operator
There are cases, such as spin-Peierls type Hamiltonians where we need to link two fields per operator. To achieve this, we will have to update the fields class (I think). Here is an explicit example where such a feature would be needed.