Write now the hopping gate only accepts list, this is slightly awkward for a single global hopping parameter. Right now one has to write:
qc.hop_fermions([theta_j], list(range(0, 4)))
It would be more natural to be able to write:
qc.hop_fermions(theta_j, list(range(0, 4)))