TeNPy Forum

Hi, This can not be done out of the box. TEBD will also be difficult, since your gates are not nearest neighbour (even if H_i are NN on their subsystem, they are NNN on the full system). You could force TEBD by grouping sites to a nearest neighbour model. I would go for MPO time evolution though. Yo...

Yes, you can have a look at the operator in the site In [3]: import tenpy as tp In [4]: site = tp.FermionSite(filling=0.5) In [5]: site.dN Out[5]: <npc.Array shape=(2, 2) labels=['p', 'p*']> In [6]: site.dN.to_ndarray() Out[6]: array([[-0.5, 0. ], [ 0. , 0.5]]) In[7]: # The basis is { |empty> , |ful...

Hi,
Yes it turns out those examples are not yet updated for v1.0.
They should run with v0.99.
Working on a fix in https://github.com/tenpy/tenpy/pull/400.
Thanks for posting!