Implementation of six-fold anisotropy.

[J]

Hi all,

I’m trying to implement six-fold in-plane anisotropy. =(λ/2)*(cos(6φ_a) + cos(6φ_B)) , this is the expression for MnTe bilayered material.

Could anyone advise:

1. How to define six-fold anisotropy (in .mat file or custom function) ?
2. ⁠Any example scripts or references?

Thanks for your help!

Best regards,
J

[gabo...@gmail.com]

Not sure how much it will help, but there is an article by J. B. Collings et al. titled "Generalized form of the magnetic anisotropy field in micromagnetic and atomistic spin models" [1] that may be of interest to you.  That is because it looks like the article implements different θ,ϕ orders of anisotropy in VAMPIRE using C++ code in .cpp files.  It looks like the associate code files, such as uniaxial_order_4.cpp, are in the magnetocrystalline-anisotropy branch of src/anisotropy at [2].

[1] https://doi.org/10.1103/PhysRevB.107.064413

[2] https://github.com/richard-evans/vampire/tree/magnetocrystalline-anisotropy/src/anisotropy

Kind Regards,

Gavin

VAMPIRE user