0 votes
asked by (180 points)
edited by

Hi everyone,

I am trying to implement a system with mixed Sites ( two bosonic sites at positions i, j, the rest fermionic spin 1/2 sites). The bosonic sites should not conserve QNs, the fermionic sites do.
I pass arguments like {"ConserveQNs",true} to the fermionic subsystem and {"ConserveQNs", false, "MaxOcc", MaxOcc} to the bosonic subsystem.

Creating the SiteSet, creating an initial state InitState(siteset) and setting the different sites also works without an problem.
However the created state does not seem to work with the function MPS(state) to create an MPS from it. It gives a segmentation fault.
Is this because of the different QN conservation on the subsystems? And is there an easy fix/ alternative way to create an MPS from the SiteSet?

Thank you very much in advance for any help!

commented by (14.1k points)
We would need to investigate this in more detail, but just a quick comment that making MPS with mixed types of sites is generally much easier and smoother in the Julia version. Do you have a particular reason for using the C++ version in this case?
commented by (180 points)
Thanks for the comment! I work with C++ as it is the program commonly used in my group and it is easier to compare with work done previously by others.

1 Answer

0 votes
answered by (180 points)

For anyone interested, I solved the issue by now by adding "empty" quantum numbers in the file ITensor/itensor/mps/sites/boson.h
Replacing one line
s = Index(1+maxOcc,tags);
s = Index(QN(),1+maxOcc,tags);
did the trick.

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.