https://github.com/tenpy/tenpy/blob/mas ... py#L13-L80

It is stated that the Hamiltonian reads

However, looking at the implementation of on-site terms I see.. math ::

H = \sum_{ij} t_{ij} b_i^\dagger b_j + \sum_i \mu (n_{A, i} - n_{B, i})

+ V \sum_{\langle ij \rangle, i<j} n_{A, i} n_{B, j}

Code: Select all

```
for u in range(len(self.lat.unit_cell)):
self.add_onsite(mu, 0, 'N', category='mu N')
self.add_onsite(-mu, 1, 'N', category='mu N')
```