SpinW / spinw

SpinW Matlab library for spin wave calculation
http://www.spinw.org
GNU General Public License v3.0
35 stars 15 forks source link

Add quartic single ion anisotropy #152

Open mducle opened 10 months ago

mducle commented 10 months ago

Currently SpinW supports a quadratic single-ion anisotropy term:

\sum_i\mathbf{S}_iA_i\mathbf{S}_i

This allows to define axial and planar anisotropy.

A quartic single-ion anisotropy term, like:

\sum_i\left(\mathbf{S}_iA_i\mathbf{S}_i\right)^2

analogous to the biquadratic interaction:

\sum_{ij}K_{ij}\left(\mathbf{S}_i\cdot\mathbf{S}_j\right)^2

would allow more complex anisotropy (such as fan or conical types) which could manifest in rare earth systems.

The implementation of the current (quadratic) anisotropy is to add the SIA term to the diagonal of the JJ interaction matrix. A similar approach might be possible for a quartic anisotropy (by adding to the bq biquadratic interaction matrix).