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>