fbuessen / SpinParser

Pseudofermion functional renormalization group solver for (frustrated) quantum magnets in two and three spatial dimensions.
MIT License
24 stars 6 forks source link

Error while implementing Pyrochlore with J3 bonds #10

Closed Sourin-chatterjee closed 2 years ago

Sourin-chatterjee commented 2 years ago

Hi,

I was trying to implement pyrochlore with third nearest neighbour interactions. ( J3 interactions are of two types, I have labelled them as j3a and j3b). The code is working fine for j3a=j3b, but if I take j3a =/ j3b, it's showing an error message "Lattice symmetry calculation has failed. (Internal error. Could not establish identity as a valid symmetry operation)".

My j3 interactions are::

[model name="pyrochlore-heisenberg-j3"]

    <interaction parameter="j3a" from="0,0,0,0" to="0,0,1,0" type="heisenberg" />
<interaction parameter="j3a" from="0,0,0,0" to="0,1,0,0" type="heisenberg" />
<interaction parameter="j3a" from="0,0,0,0" to="1,0,0,0" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,0" to="0,-1,1,0" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,0" to="-1,1,0,0" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,0" to="-1,0,1,0" type="heisenberg" />

<interaction parameter="j3a" from="0,0,0,1" to="0,0,1,1" type="heisenberg" />
<interaction parameter="j3a" from="0,0,0,1" to="0,1,0,1" type="heisenberg" />
<interaction parameter="j3a" from="0,0,0,1" to="1,0,0,1" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,1" to="0,-1,1,1" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,1" to="-1,1,0,1" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,1" to="-1,0,1,1" type="heisenberg" />

<interaction parameter="j3a" from="0,0,0,2" to="0,0,1,2" type="heisenberg" />
<interaction parameter="j3a" from="0,0,0,2" to="0,1,0,2" type="heisenberg" />
<interaction parameter="j3a" from="0,0,0,2" to="1,0,0,2" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,2" to="0,-1,1,2" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,2" to="-1,1,0,2" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,2" to="-1,0,1,2" type="heisenberg" />

<interaction parameter="j3a" from="0,0,0,3" to="0,0,1,3" type="heisenberg" />
<interaction parameter="j3a" from="0,0,0,3" to="0,1,0,3" type="heisenberg" />
<interaction parameter="j3a" from="0,0,0,3" to="1,0,0,3" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,3" to="0,-1,1,3" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,3" to="-1,1,0,3" type="heisenberg" />
<interaction parameter="j3b" from="0,0,0,3" to="-1,0,1,3" type="heisenberg" />   

[/model]

Sourin-chatterjee commented 2 years ago

Hi Finn, The issue is resolved.

Thanks