Thanks for your support so far! I have a question about creating a mixed site set involving spin and bosons。I have recently started learning physics and ITensor, so I have limited understanding of how to represent mixed-type MPOs (Matrix Product Operators). Could you provide me with some guidance or suggestions to help me better understand and represent these mixed-type MPOs?
make a site array that contains site indices of different site types (site types are special Index tags such as “S=1/2” or “Boson” that our system looks for in order to make operators)
use our OpSum system to make your Hamiltonian, and input appropriate strings of operators that you want to make or include in the sum
where I’ve made some slightly different choices than in your math above, the main one being that I’ve made the spins live on different sites than the bosons, instead of in your math where you treat them as living on the same site. (Both are mathematically ok, but it’s practically more convenient and efficient in DMRG / MPS methods to split them across different sites.)
Finally, please note that if you use “Boson” sites in ITensorMPS, that they come with a default maximum occupancy. If you need a larger occupancy for your system you should increase it.
You can also make site arrays in a more “manual” way by making the site indices one-by-one in a loop and push! -ing them into an array, which lets you control more easily the details of each site index.