Alternating Heisenberg coupling at infinite chain

[Yu-Jie Liu]

Hi,

I wanted to implement infinite DMRG to solve an alternating Heisenberg chain with nearest neighbour couplings J_ on the even and J on the odd sites. I created a CouplingModel class and used

Code: Select all

self.add_coupling([J_, J],0,"Sigmaz",0,"Sigmaz",[1])
self.add_coupling([J_,J], 0,"Sigmax",0,"Sigmax",[1])
self.add_coupling([J_,J], 0,"Sigmay",0,"Sigmay",[1])

Is this correct? I only specify the coupling strengths for the first two sites but my understanding is TeNPy will lift this to fill the entire chain.

Thanks,
Yu-Jie

[Johannes]

Hi,
yes, that is correct. It works if the number of couplings is even (ie. odd number of sites for finite systems, even number of sites for infinite MPS) , otherwise it complains about incommensurate arrays for the expansion.

[Inderpreet Kaur]

Hi, I am very new to tenpy. I am sorry if answer to this question is already there but I could not find any particular solution anywhere. I need to incorporate alternate strengths in one term of the Hamiltonian and compute it's ground state using the finite matrix-product state (MPS) representation. If I consider even coupling and even number of sites, then for e.g., the following statement works fine:

self.add_coupling([-t,0],0,'Bd',1,'B',0,plus_hc=True)#For even lengths, say L=40

It is also according to the following link,
"https://tenpy.readthedocs.io/en/latest/ ... d_coupling"

However, I am getting incommensurate error when I consider odd number of sites and finite MPS, opposite to what you have suggested in the previous comment. Can you please suggest, if there is any solution when one considers odd number of sites and finite MPS (with the motivation that the strengths of coupling is symmetric throughout the whole chain if one considers open boundaries)?