+1 vote
asked by (370 points)
closed by

Hi there,

I wanna run tDMRG for Anderson-Holstein model.

Because the initial state is with zero boson, at the step of DMRG, I restrict "ConserveNf=true, ConserveSz=true, ConserveNb=true", and set initial boson number to 0, so that I can find the ground state of the decouped boson-fermion system easily.

The MPS I got from DMRG step are with all those QNs for fermion and Boson. But after the quench, the boson number is no longer conserved, and will grow, so I need to remove the conservation number associated with boson. But I still want to keep the conservation related to fermion, "ConserveNf=true, ConserveSz=true".

I know I should use removeQNs to remove QNs from MPS, but how should I remove QNs partially?

Thank you very much!

closed with the note: Done

1 Answer

0 votes
answered by (70.1k points)
selected by
Best answer

This is a great idea for a feature we should have, so I opened an issue for it on our issue tracker page on Github.

However, unfortunately we will not be able to add this feature soon, as once we thought about the technical steps they are somewhat involved and we are also working to complete the Julia version of ITensor.

As a short-term fix, would you be able to just not conserve boson number in your initial ground state calculation? Maybe the ground state of that Hamiltonian has zero boson number anyway? Or if not, maybe you can include a chemical potential term to tune to the zero-boson-number ground state case?

You can always measure the total expected boson number after DMRG of course, to check.

Best regards,

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.