Hi,

Thank you so much for your mail. My input for kitaev-heisenberg is
LATTICE_LIBRARY="lattices.xml"
MODEL_LIBRARY="models.xml"

MODEL="kitaev heisenberg"
LATTICE="honeycomb lattice"

local_S = 0.5
W=3
J = 1
Jx = 0
Jy = 0
Jz = 0
{L=3;}

though i havnt changed the lattices.xml file. And for the spin

MODEL="spin"
LATTICE="honeycomb lattice"

local_S = 0.5
W=3
J = 1
{L=3;}

Regards

Sriluckshmy


On Mon, Jun 23, 2014 at 4:37 AM, Matthias Troyer <troyer@phys.ethz.ch> wrote:
Can you send the input files that you used?

On 22 Jun 2014, at 12:03, Sriluckshmy Viswanathan <sriluckus@gmail.com> wrote:

> Hi,
>
> I am writing a code for the computing gap in kitaev-heisenberg model on the honeycomb lattice using sparse diagonalization. When the coupling for the kitaev term is zero the model patch that I wrote should match the model that is inbuilt. But that doesnt happen. I may have made a mistake in the model patch which is as follows.
>
> <HAMILTONIAN name="kitaev heisenberg">
> <PARAMETER name="J" default="1"/>
> <PARAMETER name="Jz" default="0"/>
> <PARAMETER name="Jx" default="0"/>
> <PARAMETER name="Jy" default="0"/>
> <BASIS ref="spin"/>
> <BONDTERM source = "i" target="j">
>  J*Sz(i)*Sz(j) +J*exchange_xy(i,j)
> </BONDTERM>
> <BONDTERM type="0" source="i" target="j">
>     Jz*Sz(i)*Sz(j) + J*Sz(i)*Sz(j) +J*exchange_xy(i,j)
> </BONDTERM>
>
> <BONDTERM type="1" source="i" target="j">
>     Jx*Sx(i)*Sx(j) + J*Sz(i)*Sz(j) +J*exchange_xy(i,j)
> </BONDTERM>
>
> <BONDTERM type="2" source="i" target="j">
>     Jy*Sy(i)*Sy(j) + J*Sz(i)*Sz(j) +J*exchange_xy(i,j)
> </BONDTERM>
> </HAMILTONIAN>
>
> I have also attached the output that I got in dat format. Could you please help me with this problem?
>
> Sriluckshmy
>
> <kitaev_heisenberg.task1.out.xml.dat><spin.task1.out.xml.dat>




--
Regards
Sriluckshmy