It should work by specifying the parameter
D="0 0 1"
The three components are Dx, Dy and Dz. Please be aware however that in that case only local updates can be used.
Matthias
On Mar 27, 2007, at 1:30 AM, David Miller wrote:
Dear all,
I am doing classical MC simulations using the Heisenberg model with Alps-1.3b2. I would like to include an easy axis/plane anisotropy in the Heisenberg Hamiltonian of the form (-/+)D*sum_i((Sz_i)^2), how would I acheive this?
Thanks,
David Miller
Solve the Conspiracy and win fantastic prizes. http:// www.theconspiracygame.co.uk/