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