Hi,Matt This is the input parameter file with your ALPS Kondo lattice model, the DMRG parameters are chosen so that the calculation can finish in a short time, thus be easy to be debugged :
--------------------------------INPUT PARAMETER FILE (FILE NAME my_parm)-----------------------------------------
LATTICE_LIBRARY="my_lattices.xml" LATTICE="chain lattice" MODEL="Kondo lattice" mu=0 t=1 J=0.2 {L=32 SWEEPS=4 MAXSTATES=40 NUMBER_EIGENVALUES=2}
-----------------------------END OF INPUT PARAMETER FILE---------------------------------
-----------------------------my_lattices.xml--------------------------------------------------------------
<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="fermion" 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="open"/> </FINITELATTICE> <UNITCELL ref="fermion"/> </LATTICEGRAPH>
</LATTICES>
-------------------------END OF my_lattices.xml----------------------------------------
-------------------------COMMAND LINE INPUT-----------------------------------------
parameter2xml my_parm
dmrg --write-xml my_parm.in.xml
------------------------END OF COMMAND LINE INPUT-------------------------------
ERROR: site basis for type 2 not found. Hope you guys can find out what is going on.
By the way, if I use wild card and my custom fermion code to reprogram Model="Kondo lattice" in the fashion of your ALPS Model="spin", the DMRG and ED both works. But the problem is with mu=0 t=1 and J=1, by DMRG, my custom Kondo Lattice model shows the gap = 0.21, with 10^-10 truncation error, while K.Ueda in Phys.Rev.B. 46, 3175 FIG.1 shows a energy gap smaller than 0.2 under the same parameter condition by ED. (BOTH HALF FILLING) I recheck my custom Kondo lattice model and find no obvious mistakes. And I do not know why, can you help me here, too?
--------------------------INPUT FILE OF CUSTOM KONDO LATTICE MODEL------------------------------
LATTICE_LIBRARY="my_lattices.xml" LATTICE="chain lattice" MODEL_LIBRARY="my_models.xml" MODEL="Kondo lattice" CONSERVED_QUANTUMNUMBERS="N" N_total=128 MEASURE_LOCAL[local density]=n local_max0=2 local_min0=0 local_max1=1 local_min1=1 mu=0 t=1 J=1 {L=64 SWEEPS=10 MAXSTATES=500 NUMBER_EIGENVALUES=2}
------------------------CUSTOM KONDO LATTICE MODEL--------------------------
<MODELS> <SITEBASIS name="fermion"> <PARAMETER name="local_maximum" default="local_max"/> <PARAMETER name="local_max" default="2"/> <PARAMETER name="local_minimun" default="local_min"/> <PARAMETER name="local_min" default="0"/> <QUANTUMNUMBER name="N" min="local_minimum" max="local_maximum" type="fermionic"/> <QUANTUMNUMBER name="Sz" min="-N*(2-N)/2" max="N*(2-N)/2"/> <QUANTUMNUMBER name="Nup" min="(N+2*Sz)/2" max="(N+2*Sz)/2"/> <QUANTUMNUMBER name="Ndown" min="(N-2*Sz)/2" max="(N-2*Sz)/2"/> <OPERATOR name="Splus" matrixelement="1"> <CHANGE quantumnumber="Nup" change="1"/> <CHANGE quantumnumber="Ndown" change="-1"/> <CHANGE quantumnumber="Sz" change="1"/> </OPERATOR> <OPERATOR name="Sminus" matrixelement="1"> <CHANGE quantumnumber="Nup" change="-1"/> <CHANGE quantumnumber="Ndown" change="+1"/> <CHANGE quantumnumber="Sz" change="-1"/> </OPERATOR> <OPERATOR name="Sz" matrixelement="Sz"/> <OPERATOR name="c_down" matrixelement="1"> <CHANGE quantumnumber="Ndown" change="-1"/> <CHANGE quantumnumber="N" change="-1"/> <CHANGE quantumnumber="Sz" change="1/2"/> </OPERATOR> <OPERATOR name="cdag_down" matrixelement="1"> <CHANGE quantumnumber="Ndown" change="1"/> <CHANGE quantumnumber="N" change="1"/> <CHANGE quantumnumber="Sz" change="-1/2"/> </OPERATOR> <OPERATOR name="c_up" matrixelement="1"> <CHANGE quantumnumber="Nup" change="-1"/> <CHANGE quantumnumber="N" change="-1"/> <CHANGE quantumnumber="Sz" change="-1/2"/> </OPERATOR> <OPERATOR name="cdag_up" matrixelement="1"> <CHANGE quantumnumber="Nup" change="1"/> <CHANGE quantumnumber="N" change="1"/> <CHANGE quantumnumber="Sz" change="+1/2"/> </OPERATOR> <OPERATOR name="n" matrixelement="N"/> <OPERATOR name="n_up" matrixelement="(N+2*Sz)/2"/> <OPERATOR name="n_down" matrixelement="(N-2*Sz)/2"/> </SITEBASIS>
<BASIS name="Kondo lattice"> <SITEBASIS ref="fermion"> <PARAMETER name="local_maximum" value="local_max#"/> <PARAMETER name="local_max#" value="local_max"/> <PARAMETER name="local_max" value="2"/> <PARAMETER name="local_minimum" value="local_min#"/> <PARAMETER name="local_min#" value="local_min"/> <PARAMETER name="local_min" value="0"/> </SITEBASIS> <CONSTRAINT quantumnumber="N" value="N_total"/> <!-- <CONSTRAINT quantumnumber="Sz" value="Sz_total"/> --> </BASIS>
<BONDOPERATOR name="exchange" source="x" target="y"> Sz(x)*Sz(y)+1/2*(Splus(x)*Sminus(y)+Sminus(x)*Splus(y)) </BONDOPERATOR>
<BONDOPERATOR name="fermion_hop" source="x" target="y">
cdag_up(x)*c_up(y)+cdag_up(y)*c_up(x)+cdag_down(x)*c_down(y)+cdag_down(y)*c_down(x) </BONDOPERATOR>
<HAMILTONIAN name="Kondo lattice"> <PARAMETER name="mu" default="0"/> <PARAMETER name="t" default="1"/> <PARAMETER name="J" default="1"/> <BASIS ref="Kondo lattice"/> <SITETERM type="0" site="i"> -mu*n(i) </SITETERM> <SITETERM type="1"> 0 </SITETERM> <BONDTERM type="0" source="i" target="j"> -t*fermion_hop(i,j) </BONDTERM> <BONDTERM type="1" source="i" target="j"> J*exchange(i,j) </BONDTERM> </HAMILTONIAN>
</MODELS>
--------------------------CHAIN LATTICE IN my_lattices.xml --------------------------------------
<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="fermion" 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="open"/> </FINITELATTICE> <UNITCELL ref="fermion"/> </LATTICEGRAPH>
</LATTICES>
---------------------------------------END-------------------------------------------------------
Best, Bo-Nan