TeNPy Forum

Hi Johannes, Thanks a lot for your help! I try this code: self.add_coupling(-J, 0, 'Bd gJW', 0, 'B', -1, plus_hc=True) and it works well. Actually, I wanted to add \hat n_j in the exponent at first instead of \hat n_{j-1} . However, what I would like to mention is that \hat n_j and \hat n_{j-1} give...

Hi, I am trying to define the Anyon-Hubbard model(AHM): H = -J {\sum}^{L}_{j=2} ({\hat b}^{\dagger}_{j} {\hat b}_{j-1} e^{i \theta {\hat n}_{j}} + H.c.) + U \sum^{L}_{j=1} {\hat n}_{j} ({\hat n}_{j} - 1) import numpy as np import matplotlib.pyplot as plt import tenpy import tenpy.linalg.np_conserved...