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,

\displaystyle E_A=K_1\sin^2\theta+K_2\sin^4\theta+K_3\sin^6\theta

and the field dependent Zeeman energy term

\displaystyle E_H=-HM_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:

\displaystyle 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:

\displaystyle \sin\theta_c={[-K_2\pm (K_2^2-3K1K_3)^{1/2}]/3K_3}^{1/2} The C + and C − are the cones associated to the + and - sign. It can be verified that the only C + is always a minimum and can be an easy direction, while C − 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 K2 / K1 and K3 / K1. The following figure shows the phase diagram for the two cases K10 and {{math| K1 − and C + cones respectively. The absolute minimum (easy direction) is indicated by filling of the tip.

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 K1 , K2 , K3 into conjugate quantities, denoted byR1 , R2 , R3. This transformation can be found in such a way that all the results obtained for the case of H parallel to c-axis can be immediately transferred to the case of H perpendicular to c-axis and vice versa according to the following symmetrical dual correspondence:

The use of the table is very simple. If we have a magnetization curve obtained with H perpendicular to c-axis and with the anisotropy constant K1, K2, K3, we can have exactly the same magnetization curve using R1, R2, R3 according to the table but applying the H 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 H 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 and on the ordinate .

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 K2 / K1 and K3 / K1. 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).

Polycrystalline system

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 (Hcr) 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 γ 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 γ between 0 and π/2. 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 H different from hard direction and then becomes smeared out. For higher value of γ the magnetization curve becomes smooth, as is evident from computer magnetization curves obtained by summation of all curves corresponding to all angles γ between 0 and π/2.

Origin of high order anisotropy constants

The origin of high anisotropy constant can be found in the interaction of two sublattices (A and B) 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:

\displaystyle \begin{align}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{align}

where:

• ,J is the exchange integral (J0) in case of ferromagnetism
• ,K_{1A},\ K_{2A} are the anisotropy constants of the A sublattice
• ,K_{1B},\ K_{2B} are the anisotropy constants of the B sublattice respectively
• ,H is the applied field
• ,M_A,\ M_B are the saturation magnetizations of A and B sublattices
• ,\theta_A,\ \theta_B are the angles between the magnetization of A and B 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:

\displaystyle 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 fourth 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 fourth order: \displaystyle E_A=K_1\sin^2\theta+K_2\sin^4\theta+K_t\sin^3\theta\cos\theta\sin3\phi where \theta and \phi 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.

References

References

1. (2008). "Determining anisotropy constants from a first-order magnetization process in Tb2Fe17". Physical Review B.
2. (1990). "Ferromagnetic materials: a handbook on the properties of magnetically ordered substances.". North-Holland.
3. (1990). "Competing anisotropies and first-order magnetization processes". Journal of Applied Physics.
4. (1 December 1936). "Determination of Ferromagnetic Anisotropy in Single Crystals and in Polycrystalline Sheets". Physical Review.
5. (1964). "Über die magnetische Anisotropieenergie kubischer Kristalle". Physica Status Solidi B.
6. (1966). "Physical Principles of Magnetism". D. Van Nostrand Company.
7. (1971). "Transition du premier ordre dans les cristaux cubiques induite par un champ magnétique orienté suivant une direction de difficile aimantation". Journal de Physique Colloques.
8. (January 1977). "The magnetization process in cubic ferromagnetic single crystals". Physica B+C.