Dear ITensor community,

I want to ask for help in defining the siteset for a two-component Bose-Hubbard Hamiltonian. Hamiltonian of the interest is as follows:

$$

\begin{equation}

H = -J\sum_i (a^\dagger_{i+1}a_i + b^{\dagger}_{i+1}b_i + h.c.) + U_{aa}/2\sum_i n^{a}_i(n^{a}_i-1) + U_{bb}/2\sum_i n^{b}_i (n^{b}_i-1) + U_{ab}\sum_{i}n^{a}_i n^{b}_i

\end{equation}

$$

where operators @@a_i@@ fulfill bosonic commutation relations, as well as @@b@@.

Can I ask for some steps of modification "bosons.h" file to construct mentioned Hamiltonian?

The additional constraints are:

1. Total number of particles is conserved N*a + N*b = const

2. I would like to consider the case that, in general, the magnetization S*z = N*a - N_b is not conserved.

I really appreciate any help you can provide.