FOMP

The magnetocrystalline anisotropy energy of a ferromagnetic crystal can be expressed as a power series of direction cosines of the magnetic moment with respect to the crystal axes. The coefficient of those terms is the constant anisotropy. In general, the expansion is limited to a few terms. Normally the magnetization curve is continuous with respect to the applied field up to saturation but, in certain intervals of the anisotropy constant values, irreversible field-induced rotations of the magnetization are possible, implying first-order magnetization transition between equivalent magnetization minima, the so-called first-order magnetization process (FOMP).

Theory

The total energy of a uniaxial magnetic crystal in an applied magnetic field can be written as a summation of the anisotropy term up to six order, neglecting the sixfold planar contribution,

$$E_A=K_1{\sin }^2\theta +K_2{\sin }^4\theta +K_3{\sin }^6\theta$$

and the field dependent Zeeman energy term

$$E_H=-H{M}_s\cos (\theta -\gamma )$$

where:

• $K_1,K_2,K_3$ are the anisotropy constants up to six order
• $H$ is the applied magnetic field
• $M_s$ is the saturation magnetization
• $\theta$ is the angle between the magnetization and the easy c-axis
• $\gamma$ is the angle between the field and the easy c-axis

The total energy can be written:

$$E_T=E_A+E_H$$

Phase diagram of easy and hard directions

In order to determine the preferred directions the magnetization vector in the absence of the external magnetic field we analyze first the case of uniaxial crystal. The maxima and minima of energy with respect to $\theta$ must satisfy

$$\frac{\partial E_A}{\partial \theta }=0$$ while for the existence of the minima

$$\frac{{\partial }^2{E}_A}{\partial {\theta }^2}>0.$$

For symmetry reasons the c-axis $(\theta =0)$ and the basal plane $(\theta =\frac{\pi }{2})$ are always points of extrema and can be easy or hard directions depending on the anisotropy constant values. We can have two additional extrema along conical directions at angles given by:

$$\sin {\theta }_c={[-K_2\pm (K_2^2-3K1K_3{)}^{1/2}]/3K_3}^{1/2}$$ The Template:Math and Template:Math are the cones associated to the + and - sign. It can be verified that the only Template:Math is always a minimum and can be an easy direction, while Template:Math is always a hard direction.

A useful representation of the diagram of the easy directions and other extrema is the representation in term of reduced anisotropy constant Template:Math and Template:Math. The following figure shows the phase diagram for the two cases Template:Math and Template:Math. All the information concerning the easy directions and the other extrema are contained in a special symbol that marks every different region. It simulate a polar type of energy representation indicating existing extrema by concave (minimum) and convex tips (maximum). Vertical and horizontal stems refer to the symmetry axis and the basal plane respectively. The left-hand and the right-hand oblique stems indicate the Template:Math and Template:Math cones respectively. The absolute minimum (easy direction) is indicated by filling of the tip.

Magnetic phase diagram of the uniaxial ferromagnet. Left: For Template:Math. Right: For Template:Math. See text for symbol explanations.

Transformation of Anisotropy Constant into Conjugate Quantities

Before going into the details of the calculation of the various types of FOMP we call the readers attention to a convenient transformation ($K\rightleftharpoons R$) of the anisotropy constants Template:Math into conjugate quantities, denoted byTemplate:Math. This transformation can be found in such a way that all the results obtained for the case of Template:Mvar parallel to c-axis can be immediately transferred to the case of Template:Mvar perpendicular to c-axis and vice versa according to the following symmetrical dual correspondence:

The $K\rightleftharpoons R$ transformation of anisotropy constants.

basal plane	c-axis	$\Leftarrow$ DUAL $\Rightarrow$	c-axis	basal plane
$K_1$	$R_1=-K_1-2K_2-3K_3$	$\rightleftharpoons$	$R_1$	$K_1=-R_1-2R_2-3R_3$
$K_2$	$R_2=K_2+3K_3$	$\rightleftharpoons$	$R_2$	$K_2=R_2+3R_3$
$K_3$	$R_3=-K_3$	$\rightleftharpoons$	$R_3$	$K_3=-R_3$

The use of the table is very simple. If we have a magnetization curve obtained with Template:Mvar perpendicular to c-axis and with the anisotropy constant Template:Math, we can have exactly the same magnetization curve using Template:Math according to the table but applying the Template:Mvar parallel to c-axis and vice versa.

FOMP examples

The determination of the conditions for the existence of FOMP requires the analysis of the magnetization curve dependence on the anisotropy constant values, for different directions of the magnetic field. We limit the analysis to the cases for Template:Mvar parallel or perpendicular to the c-axis, hereafter indicated as A-case and P-case, where A denotes axial while P stands for planar. The analysis of the equilibrium conditions shows that two types of FOMP are possible, depending on the final state after the transition, in case of saturation we have (type-1 FOMP) otherwise (type-2 FOMP). In the case when an easy cone is present we add the suffix C to the description of the FOMP-type. So all possible cases of FOMP-type are: A1, A2, P1, P2, P1C, A1C. In the following figure some examples of FOMP-type are represented, i.e. P1, A1C and P2 for different anisotropy constants, reduced variable are given on the axes, in particular on the abscissa Template:Math and on the ordinate Template:Math.

Examples of FOMP-types. Left: For Template:Math,Template:Math, Template:Math. Center: for Template:Math,Template:Math, Template:Math. Right: For Template:Math,Template:Math, Template:Math for Template:Math with Template:Math and Template:Math. See text for label explanations.

FOMP diagram

Tedious calculations allow now to determine completely the regions of existence of type 1 or type 2 FOMP. As in the case of the diagram of the easy directions and other extrema is convenient the representation in term of reduced anisotropy constant Template:Math and Template:Math. In the following figure we summarize all the FOMP-types distinguished by the labels A1, A2, P1, P2, P1C, A1C which specifies the magnetic field directions (A axial; P planar) and the type of FOMP (1 and 2) and the easy cone regions with type 1 FOMP (A1C, P1C).

FOMP phase diagram of the uniaxial ferromagnet. Left: For Template:Math. Right: For Template:Math. See text for label explanations.

Polycrystalline system

Computer simulation of the magnetization of polycrystal sample with Template:Math,Template:Math, Template:Math (Template:Math, red curve). The green is the first derivative Template:Math. While the black ones are the magnetization Template:Math of the corresponding single crystal with Template:Mvar applied along c-axis Template:Math, Template:Math and Template:Math.

Since the FOMP transition represents a singular point in the magnetization curve of a single crystal, we analyze how this singularity is transformed when we magnetize a polycrystalline sample. The result of the mathematical analysis shows the possibility of carrying out the measurement of critical field (Template:Math) where the FOMP transition takes place in the case of polycrystalline samples.

For determining the characteristics of FOMP when the magnetic field is applied at a variable angle Template:Mvar with respect to the c axis, we have to examine the evolution of the total energy of the crystal with increasing field for different values of Template:Mvar between Template:Mvar and Template:Math. The calculations are complicated and we report only the conclusions. The sharp FOMP transition, evident in single crystal, in the case of polycrystalline samples moves at higher fields for Template:Mvar different from hard direction and then becomes smeared out. For higher value of Template:Mvar the magnetization curve becomes smooth, as is evident from computer magnetization curves obtained by summation of all curves corresponding to all angles Template:Mvar between Template:Mvar and Template:Math.

Origin of high order anisotropy constants

The origin of high anisotropy constant can be found in the interaction of two sublattices (Template:Mvar and Template:Mvar) each of them having a competing high anisotropy energy, i.e. having different easy directions. In particular we can no longer consider the system as a rigid collinear magnetic structure, but we must allow for substantial deviations from the equilibrium configuration present at zero field. Limiting up to fourth order, neglecting the in plane contribution, the anisotropy energy becomes:

$$\begin{array}{cc}{E}_t& =-J\overrightarrow{M_A}\cdot \overrightarrow{M_B}+K_{1A}{\sin }^2{\theta }_A+K_{1B}{\sin }^2{\theta }_B\\ & +K_{2A}{\sin }^4{\theta }_A+K_{2B}{\sin }^4{\theta }_B-(\overrightarrow{M_A}+\overrightarrow{M_B})\cdot \overrightarrow{H}\end{array}$$

where:

• $J$ is the exchange integral (Template:Math) in case of ferromagnetism
• $K_{1A},K_{2A}$ are the anisotropy constants of the Template:Mvar sublattice
• $K_{1B},K_{2B}$ are the anisotropy constants of the Template:Mvar sublattice respectively
• $H$ is the applied field
• $M_A,M_B$ are the saturation magnetizations of Template:Mvar and Template:Mvar sublattices
• ${\theta }_A,{\theta }_B$ are the angles between the magnetization of Template:Mvar and Template:Mvar sublattices with the easy c-axis

The equilibrium equation of the anisotropy energy has not a complete analytical solution so computer analysis is useful. The interesting aspect regards the simulation of the resulted magnetization curves, analytical or discontinuous with FOMP. By computer it is possible to fit the obtained results by an equivalent anisotropy energy expression:

$$E_A=K_1^{*}{\sin }^2\theta +K_2^{*}{\sin }^4\theta +K_3^{*}{\sin }^6\theta$$

where:

• $K_1^{*},K_2^{*},K_3^{*}$ are the equivalent anisotropy constants up to six order
• $\theta$ is the angle between the magnetization and the easy c-axis

So starting from a forth order anisotropy energy expression we obtain an equivalent expression of sixth order, that is higher anisotropy constant can be derived from the competing anisotropy of different sublattices.

FOMP in other symmetries

The problem for cubic crystal system has been approached by Bozorth, and partial results have been obtained by different authors, but exact complete phase diagrams with anisotropy contributions up to sixth and eighth order have only been determined more recently.

The FOMP in the trigonal crystal system has been analyzed for case of the anisotropy energy expression up to forth order: $$E_A=K_1{\sin }^2\theta +K_2{\sin }^4\theta +K_t{\sin }^3\theta \cos \theta \sin 3\varphi$$ where $\theta$ and $\varphi$ are the polar angles of the magnetization vector with respect to c-axis. The study of the energy derivatives allows the determination of the magnetic phase and the FOMP-phase as in the hexagonal case, see the reference for the diagrams.