Hans Fangohr, Computational Science
  • Software home
  • Pages
  • Software
  • fidimag
  • Ubermag
  • nmag
  • nsim
  • ovf2mfm
  • ovf2vtk
  • OOMMF CubicAnis8
  • OOMMF UniAnis4
  • pycaml
  • vdsim
  • OCaml Tutorial
  • Home HF
  • |
  • Software
  • |
  • OOMMF UniAnis4

Oxs Extension Module: Southampton_UniaxialAnisotropy4

The previous version of this class was called CED_UniaxialAnisotropy.


Description:

This is a generic Oxs extension object, derived from the Oxs_Energy class. It enhances the Oxs_UniaxialAnisotropy class as part of the standard OOMMF package by a fourth order anisotropy energy term. The energy density is computed by

E = -K1(u·m)2 - K2(u·m)4,

where m is the reduced (unit) magnetization, u is the easy axis, and K1/K2 are the second/fourth order phenomenological anisotropy constants.

The corresponding anisotropy field H is given by the relation

µ0MsH = 2K1(u·m)u + 4K2(u·m)3u,

where Ms is the saturation magnetization.

For numerical reasons, the anisotropy energy density is calculated in the easy axis case (K1>0) using an equivalent form with the cross product:

E = (K1 +2 K2) |u×m|2 - K2|u×m|4.

Installation:

Download the header and source code files below, and follow the general Oxs extension installation instructions.

Usage:

MIF 2.x files written to use this class should include a Specify block of the form

Specify Southampton_UniaxialAnisotropy4:name {
K1 k1_value
K2 k2_value
axis anisotropy_axis
}

The values for the K1 and K2 parameters should be scalar field objects, and axis should be a vector field object. The only difference with respect to the stock Oxs_UniaxialAnisotropy class is the inclusion of the K2 term.

Details:

  • Authors: Jürgen Zimmermann, Richard Boardman, and Hans Fangohr
  • Affiliation: School of Engineering Sciences, University of Southampton
  • Oxs_Ext class: Southampton_UniaxialAnisotropy4
  • OOMMF releases: 1.2a3
  • External libraries: none
  • License: Public Domain
  • Release date: 10/04/2007
  • Version: 1.1 (Version 1.0 was named CED_UniaxialAnisotropy)

Download:

  • C++ header file: uniaxialanisotropy4.h
  • C++ main source file: uniaxialanisotropy4.cc
  • Example MIF input files:
    • K2 = 0.0 K1: quadrupole00pc.mif
    • K2 = 0.1 K1: quadrupole10pc.mif
    • K2 = 0.3 K1: quadrupole30pc.mif

Sample results:

Output from the three example MIF files, illustrating the effect of increasing K2 relative to K1:

M vs. B, for
             various K2


DISCLAIMER: This software is free to use. However, the authors do not assume responsibility whatsoever for its use, and make no guarantees, expressed or implied, about its quality, reliability, or any other characteristic.


Return to Top

Created with reStructuredText and Sphinx.