I know that the wave function should be invariant under the rotation/permutation and so I need to check that the ov returned by compute_K is close to one. However, I am not sure how to define the correct permutation.
e.g. for an initial 3x3 Neel product state unit cell in the 2D iDMRG, with a default ordering
Does this correspond to the identity?
Code: Select all
psi.compute_K(perm=[0, 1, 2, 3, 4, 5, 6, 7, 8])
Can you give a lattice like this?
Code: Select all
model_params = dict(cons_N='N', cons_Sz='Sz', t=t, U=U, mu=mu, V=V, lattice="Square",
bc_MPS='infinite', order='default', Lx=3, Ly=3, bc_y='cylinder', verbose=0)
M = FermionicHubbardModel(model_params)
psi.compute_K(perm=M.lat)
Code: Select all
psi.compute_K(perm=[8, 0, 1, 2, 3, 4, 5, 6, 7])
Code: Select all
psi.compute_K(perm=[2, 0, 1, 5, 3, 4, 8, 6, 7])