Dear All, I am struggling with converging to the ground state of a BH model with unit filling with mps_optim with PBC. For low hopping it produces the same result as in the OBC case ignoring the boundary condition. I have the following input <pasted below>. The state to which the mps_optim converges: - has the same energy as the OBC hamiltonian for t<0.08 (the value of the constant varies it is larger for larger L) (see enclosed plot) and density artifacts at edges (similar to the OBC ground state) - If I observe how the energy of a minimized state varies as a function of sweep number, then there are “plateaus” (for several iterations the energy seems to be converged, then it jumps down).. It seems that for small t the GS computation of mps_optim is stuck at the first plateau. So for say t<0.07 after each sweep starting from 5 the energy is constant but for larger t, at say sweep number 10 there is a sharp drop. - none of the above happens for the OBC. - I include file containing the difference between the numerical result and the essentially exact result (strong coupling expansion up to 16th order - actually I copied the formula only until 10th order). Clearly for all t between 0 and 0.06 the results are wrong, then they are fine again (and start to depart for larger t again but this is because I was to lazy to copy higher orders too) - i have checked MAXSTATES=180 as well as MAXSTATES=320 (the results are exactly the same). For MAXSTATES=600 (and 120 sweeps and ietl_jcd_toll decreased to 1e-10) the results are almost the same again but now problems are for t<0.08 not for t<0.07 so in some sense situation is worse. *************** questions ****************** I have thus two questions: 1. Is there a way to resolve it? 2. It seems that propagation of the PBC across the sample length (L=40 sites) gets truncated somehow. Still increasing number of states to 600 did not help. Is it what is going on or is it a bug? Best, Mateusz Łącki ************* conf ****************** LATTICE="chain lattice" MODEL="boson Hubbard" L=40 MAXSTATES=180 measure_each=5 init_state="local_quantumnumbers" initial_local_N='1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1' nsweeps=60 t=0.06; U=1.0; N_total=40 Nmax=5; CONSERVED_QUANTUMNUMBERS='N' MEASURE[Entropy]=1 MEASURE_LOCAL[Local Density]='n' {}
