0 votes
asked by (190 points)

Dear Sir / Ma'am,

I was calculating the ground state enrgy for extended hubbard model using the code provided in sample folder. I used the range for U1 to be -8 to 8, and the range for V1 (nearest neighbour spin-spin interaction) to be in the range -6 to 6. in the above range, for some region of negative V1, the energy was not converging. It will be very helpful if someone can tell me what is going wrong here.

Some other points about the parameters I used:

  1. I used a lattice site of 10, with half-filling

  2. I used periodic boundary condition

  3. The sweeps were following:

nsweeps = 60
sweeps
{
maxdim mindim cutoff niter noise
50 10 1E-12 2 1E-7
100 20 1E-12 2 1E-8
200 20 1E-12 2 1E-10
400 20 1E-12 2 0
800 20 1E-12 2 1E-11
800 20 1E-12 2 0
}

Best,
Sreetama

commented by (46.8k points)
Hi Sreetama,
Thanks for the question. I think before we can provide a good answer to you, we'll need to ask some questions of you.

1. what is the reason for choosing periodic boundary conditions? You may know that DMRG does not work nearly as well for periodic compared to open boundary conditions. It can make a very large difference in the ability of DMRG to give accurate results. Would you be able to use open boundary conditions in your study or is there a special reason you need periodic?

2. when you say that the results were not converging, what did you see in the output? Was the energy still changing by a large amount in the last few sweeps? Or were you not reaching the truncation error you desired? Or what was the criterion for convergence you were using? If you could show a printout of the last few sweeps that would be useful to see.

Best regards,
Miles
commented ago by (190 points)
Hi Miles,

thanks for your reply. I will try to answer your questions below.

1. Yes I know that DMRG does not work well with periodic boundary conditions. I thought for small lattice size, it should not be a problem. But after the ground state energy did not converge for some values of V, I removed the periodic boundary condition. But the non-convergence still persisted. The PBC is not necessary for my work, so I can work without it.

2. The energy was changing in the 6th or 7th decimal place. When I calculated the entanglement, it was showing an absurd change with respect to nearby regions.


3. The energy is not converging even in some regions where U and V both are positive.
Sincerely,
Sreetama

1 Answer

0 votes
answered ago by (920 points)

Hi Streetama,

When I was computing E.H.M. ground states for V > U > 0, which nearly approaches a product state as V >>> U, the entanglement goes down.

This sounds to us as a good thing, but in fact with DMRG it is a bad thing because it takes an extremely large number of fast sweeps in order to converge. In essence, the small number of nonzero entries in the SVD leads to poor communication of information across the lattice with each sweep.

For instance, I believe I needed to run hundreds of sweeps to get the energy to settle down, even for small system sizes, and I observed the entanglement entropy changed very slowly with each sweep.

In order to test this, I'd recommend changing your sweeps parameters so that you stay at a single maxdim for several hundred sweeps and wait until the energy settles down before increasing maxdim.

On the other hand, if your boundary conditions are incorrect (for instance, anti-periodic instead of periodic) you can also see wild changes in sweep to sweep energy; I experienced this error one time due to a mistake in my definition of the Hamiltonian.

If you figure it out please post back here as I'm interested to know what you find, as I am continuing to use these methods.

Good luck!

commented ago by (190 points)
Dear Jon,

So I was trying to get the ground state for U=0.4 and V=-6, and used about 600 sweeps with the final ones having maxdim 2400. The energy converged only upto 5th decimal place, but of course for higher sweeps it will change in the 5th decimal place as well. So, I can say that the energy was not at all converged. As a result, when I am calculating two-site reduced density matrix, it is providing a wrong state (because the probabilities are becoming negative).

I think the region you mentioned i.e. V>U>0, will have the same convergence problem.

Kindly let me know if you have any other insight about this. Were you able to get an well converged energy for V>U>0? Was the state a valid one?

Regards,
Sreetama
Welcome to ITensor Support Q&A, where you can ask questions and receive answers from other members of the community.

Formatting Tips:
  • To format code, indent by four spaces
  • To format inline LaTeX, surround it by @@ on both sides
  • To format LaTeX on its own line, surround it by $$ above and below
  • For LaTeX, it may be necessary to backslash-escape underscore characters to obtain proper formatting. So for example writing \sum\_i to represent a sum over i.
If you cannot register due to firewall issues (e.g. you cannot see the capcha box) please email Miles Stoudenmire to ask for an account.

To report ITensor bugs, please use the issue tracker.

Categories

...