perm parameter in compute_K function
Posted: 05 Feb 2019, 11:17
What is the correct way to define perm in the compute_K function?
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?
The documentation states that "If a lattice is given, we use it to read out the lattice structure and shift each site by one lattice-vector in y-direction".
Can you give a lattice like this?
Does the documentation mean shifting the chain by one site, like this...
...or literally shifting the labels in the unit cell by one in the y-direction...
...or something else entirely? Occasionally the wave function is invariant under the permutations that I choose, but I can't see any pattern to it yet. Any clarifications would be greatly appreciated
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])