Dear Alps Users,
                          I am trying to simulate a binary ferromagnetic system i.e FeCo. I am running the classical monte carlo simulations by using the classical Heisenberg model . I am using the spinmc binary for my task.
            
                   My confusion is regarding the double counting aspect of the Heisenberg model. since i  have FeCo as the binary system. I will have two set of edges , which are as follows:
Set 1 :  { Fe-Fe, Fe-Co, Co-Co}

Set 2 : {Fe-Fe, Fe-Co, Co-Fe, Co-Co}.

My question is, In order to simulate the FeCo system correctly should i consider the edge Fe-Co ( J= 25 meV) only once OR should i consider both F e-Co and Co-Fe (J= 50 mev).

I have tried to locate the exact equation of the Heisenberg Hamiltonian which is being used for the classical monte carlo simulation, but unfortunately i could not find this information. is the double counting aspect already taken care in the ALPS simulations. if so, then i think i have to define both Fe-Co and Co-Fe edges for the simulation. I am kindly requesting you to provide me insight into this aspect.
 
I am caught up with this confusion, which is stopping me from proceeding further. Therefore, I am kindly requesting you to provide me your suggestions and help regarding the above mentioned issue.

I thank you in advance for your kind consideration and time.

with best regards
Sunil