Hi,all I am trying to sparsediag the Kondo lattice in ALPS, but the command line return "Called get_quantumnumber with illegal index". I feel that there may be a misunderstanding of the constraint in my code, but I just can not find it out now. So, dear users and developers please help me!!!
The parameter file is,
LATTICE="chain lattice" MODEL="Kondo lattice" LATTICE_LIBRARY="mylattices.xml" MODEL_LIBRARY="mymodels.xml" CONSERVED_QUANTUMNUMBERS="Sz" Sz_total=0 t=1 J=2 {L=4}
the "Kondo lattice" code in mymodels.xml is almostly copied from the "Kondo lattice" code in ALPS models.xml, the "only" change I made (no changes of sitebasis, operators or hamiltonian) is adding a constraint of the total Sz (Sz_spin+Sz_fermion) to the BASIS as shown below:
<BASIS name="Kondo lattice"> <SITEBASIS type="0" ref="fermion"/> <SITEBASIS type="1" ref="spin"/> <CONSTRAINT quantumnumber="Sz" value="Sz_total"/> </BASIS>
The lattice in mylattices.xml is defined as:
<LATTICES>
<LATTICE name="chain lattice" dimension="1"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE> <UNITCELL name="hf" dimension="1"> <VERTEX type="0"/> <VERTEX type="1"/> <EDGE type="0"><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0"/><TARGET vertex="2" offset="0"/></EDGE> </UNITCELL> <LATTICEGRAPH name = "chain lattice"> <FINITELATTICE> <LATTICE ref="chain lattice"/> <PARAMETER name="L" default="4"/> <EXTENT size="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="hf"/> </LATTICEGRAPH>
</LATTICES>
where VERTEX type=0 stands for fermion and VERTEX type=1 stands for spin.
Any advice is welcome!
Bo-Nan