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"/>

</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

--
Stay foolish,Stay hungry.