Hello all, my name is Axel! I am currently working on a lattice model essentially consisting of a 1D chain of cavities with bosonic degrees of freedom, connected via two-level systems (qubits), like Q-C-Q-C-Q-C....

I need to perform a DMRG calculation on the system. As I understand it, I am supposed to define a SiteSet describing the combined degree of freedom of an adjacent cavity and qubit, then use AutoMPO to define the Hamiltonian. Presumably, I could somehow also construct the site set and the operators by hand.

My question is: what is the best way of using ITensor to do this? Am I missing something here? It would be really useful if there was some kind of Kronecker product which I could use to define the MPO:s for example, since there is already a SiteSet defined for spin 1/2. That would make the pre-defined SiteSets much more flexible, I guess.

I also want to take the opportunity to thank the ITensor developers: I think that you have developed a comprehensible and very promising tensor library, and look forward to using and contributing to it in the future.