TeNPy Forum

I'm having issues with iDRMG convergence, and was hoping someone might have an idea of what the issue is. I'm modeling the following 1D Rydberg Chain Hamiltonian: H = \Sigma [ \Omega (b_j + b_j^+) + V_1 n_j n_{j+1} - \Delta n_j] with the following representations for operators in the Pauli basis: X_...

I'm new at TeNPy and am trying to create a model that is a CouplingMPOModel, with a new operator Nj. class myTFLIMModel(CouplingMPOModel): def init_sites(self, model_params): site = SpinHalfSite(conserve=None) Nj = np.array([[1,0],[0,0]]) site.add_op('Nj', Nj) return site def init_terms(self, model_...