TeNPy Forum

Hi Johannes, thanks for this extensive answer! The example which I gave might not be the best one. I work in a sector of a fixed number of particles which equals half the number of sites. I initialise my state with initial_state = [0]*(L//2) + [1]*(L//2) wich effectively gives half a particle per si...

I use bose-hubbard model. I want to construct an MPS that will be in a state of a superposition of 0 and 1 boson at each site. I also want to use this state in a DMRG computation that will conserve number of particles. Here is an example of what I mean: from tenpy.models.bose_hubbard_chain import Bo...