3. 2D & 3D Coherent Rotations
Author: Prof. Alexandru STANCU, Al. I. Cuza University, Iasi, Romania


3. 2D & 3D coherent magnetization rotation model

We shall add to the previous presentation a vector generalization of the Stoner-Wohlfarth model, developed by Thiaville [8].

To start, one presents the 2D generalization. The same principles will apply to the 3D generalization. The total free energy density will be written in the general form:

where is the anisotropy energy density, which depends in the 2D case on the angle between the easy axis and the polarization vector. This expression could be written as:
if one makes the notations:
where is a constant (in the uniaxial case the anisotropy constant), and is a unit vector on polarization vector direction.
The equilibrium and stability conditions are obtained deriving (2); one obtains:

(e) (8)
(s) (9)

and is a unit vector, perpendicular to .

If one consider the inverse problem (for a given find the locus of the applied field vector tip), (e) and (s) represent the equations of two straight lines:
(e) is a line perpendicular to the vector (parallel to the vector)
(s) is a line perpendicular to the vector (parallel to the vector).
So, the (e) and (s) lines are reciprocally perpendicular, (e) being parallel to the equilibrium orientation of the polarization vector. The system defining the critical curve is given by:
and has the solution:
. (12)
Example: For the uniaxial case:

. (13)
which represents the equation of an astroid.
If one makes the derivative of the in order to find the direction of the tangent to the critical curve, one finds:
that is, the tangent is parallel to the equilibrium orientation of the polarization vector, as the (e) line. We can introduce an orientation of the critical curve: the positive sense of the critical curve may be given by.

Prove that the orientation is reversing in the cusp points of the critical curve.

One observes that all the properties we derived for the uniaxial case in the previous sections applies to the general 2D case.

Constant energy curves

The locus of the applied field vector tip for one given value of the energy, w, in the equilibrium state is the solution of the system:
that is,
. (16)


# MAPLE program
> restart: with(plots):
> wa := sin(theta)^2:
> wa_d1 := diff(wa, theta):
> wa_d2 := diff(wa_d1, theta):
> h0c_x := -(wa_d1*sin(theta)+wa_d2*cos(theta))/2:
> h0c_y := +(wa_d1*cos(theta)-wa_d2*sin(theta))/2:
> p1 := plot([h0c_x, h0c_y, theta=0..2*Pi]):
> h0w_x := ((wa-w)*cos(theta)-wa_d1*sin(theta))/2:
> h0w_y := ((wa-w)*sin(theta)+wa_d1*cos(theta))/2:
> w := 0.01:
> p2 := plot([h0w_x, h0w_y, theta=0..2*Pi], color=blue):
> w := 'w':
> display([p1,p2]);

Fig. 15 Exterior (astroid) – critical curve, Interior – constant energy curve.

The tangent to the constant energy curve is given by:
that is, it is perpendicular to the equilibrium orientation, (see Figure 16).
Defining by the arc-length on the constant energy curve we have:

Fig. 16

where is the curvature radius of the constant energy curve.
In order to find the position of the curvature centre, noted with, we can use the relation (see Figure 16):
. (19)
Introducing in (19) the expressions (16) and (18) one obtains:
that is, the curvature centre is on the critical curve. From the mathematical point of view one may say that the critical curve is the common evolute of all constant energy curves.

Geometric construction for energy calculations

Thiaville [9] also proved that it is possible to use a simple geometric construction to estimate the free energy density when a field is applied to the ferromagnetic particle. The normalized total free energy density may be expressed as:
and if one uses (11), one obtains:
Taking into account (14), the arc-length on the critical curve is given by:
which is equivalent to:
. (24)
In conformity with (24) the difference between the values of in two equilibrium states on the critical curve characterized by q0 and q is twice the arc-length measured on the critical curve. This property could be used in (22) for expressing the sum between the first and the second term. The final term in (22) is minus twice the length of the segment on the equilibrium line (tangent to the critical curve) linking the point on the critical curve with the tip of the applied field vector (see Figure). In fact, in this way one can calculate the difference of energy for two states characterized by two values of the applied field. For example, when the field is applied as presented in the Figure, the equilibrium (stable) orientation of the polarization vector is the orientation of the vector . Using the consequences of the previous geometrical discussion we can say that the difference of energy between the “C” state and the “A” state (see Figure, the “A” state is characterized by an applied field vector whose tip is on the point “A”) is twice the sum of the lengths of the arc AB and the segment BC, that is:

Fig. 17

The points characterized by (A, see Figure 18) and (B) have the same energy; the points A and B could be connected by two quarters of the astroid one positive (AC) and one negative (CB).

Fig. 18


When the tip of the applied field vector is inside the astroid there are three tangents to the critical curve passing through it: two for stable equilibrium states and one for an unstable equilibrium state (see Figure 19).

3D coherent magnetization rotation model

The method applied in the previous section could be used for developing a similar set of geometrical rules for the 3D coherent magnetization model [10].
The most general form of the free energy density (normalized to 2K this time in order to obtain a simpler form) is:
is the unit vector of the polarization vector.
To obtain the stable equilibrium conditions we need the following expressions:
In the calculus we used:
where is a set of three unit vectors which are an orthogonal basis.
From the equilibrium conditions:
one can find a general expression for the vector. From (38) and (39), which represent the equations of two planes, we can calculate the components on and directions:
The intersection of these planes is the straight line oriented on the direction. The general solution of the system (38)-(39) is:
If (,) are the spherical angles defining the orientation of the polarization vector at the equilibrium the free energy density could be developed in the Taylor series near this equilibrium orientation as follows:
, (44)
. (45)
Taking into account the equilibrium conditions (38) and (39) the stability condition,
may be written as:
or, in a simpler form;
In order to have (48) for any value of x, the determinant should be negative and the coefficient of should be positive:
Introducing the equilibrium condition given by (42) in the first equation in (50) we are founding two solutions:
From the other equation in (50) we can say that the stable equilibrium part of the line defined by (42) is given by the condition:
. (52)
Introducing in (42) one finds the equations of two surfaces, .


# Maple program: Surfaces: Case 1: biaxial anisotropy
> restart: with(plots):
> px := sin(theta)*cos(phi):
> py := sin(theta)*sin(phi):
> pz := cos(theta):
> wa := px^2+0.5*py^2:
> wa_t := diff(wa, theta):
> wa_p := diff(wa, phi):
> wa_tt := diff(wa_t, theta):
> wa_pp := diff(wa_p, phi):
> wa_tp := diff(wa_t, phi):
> wa_pt := diff(wa_p, theta):
> c1 := wa_tt+wa_t*cos(theta)/sin(theta)+wa_pp/sin(theta)^2:
> c2 := wa_tt-wa_t*cos(theta)/sin(theta)-wa_pp/sin(theta)^2:
> c3 := wa_tp/sin(theta)-wa_p*cos(theta)/sin(theta)^2:
> l_plus := (1/2)*(-c1+sqrt(c2^2+4*c3^2)):
> px := sin(theta)*cos(phi):
> py := sin(theta)*sin(phi):
> pz := cos(theta):
> ux := cos(theta)*cos(phi):
> uy := cos(theta)*sin(phi):
> uz := -sin(theta):
> vx := -sin(phi):
> vy := cos(phi):
> vz := 0.0:
> hx := l_plus*px + wa_t*ux + wa_p/sin(theta)*vx:
> hy := l_plus*py + wa_t*uy + wa_p/sin(theta)*vy:
> hz := l_plus*pz + wa_t*uz + wa_p/sin(theta)*vz:
> plot3d([hx, hy, hz], theta=0..Pi, phi=0..2*Pi, grid = [50,50],axes=BOXED, orientation = [40, 80], color = BLUE);

> l_minus := (1/2)*(-c1-sqrt(c2^2+4*c3^2)):
> hx := l_minus*px + wa_t*ux + wa_p/sin(theta)*vx:
> hy := l_minus*py + wa_t*uy + wa_p/sin(theta)*vy:
> hz := l_minus*pz + wa_t*uz + wa_p/sin(theta)*vz:
> plot3d([hx, hy, hz], theta=0..Pi, phi=0..2*Pi, grid = [50,50], axes=BOXED, orientation = [60, 80], color = BLUE);

# Case 2: Cubic anisotropy
> wa := (px*py)^2 + (py*pz)^2 + (pz*px)^2: