CC BY
11
УДК 53.537.6
S. B. Leble, A. Yu. Chychkalo
HYSTERESIS LOOPS FOR A BULK FERROMAGNETIC BY HEISENBERG MODEL
In a framework of Heisenberg theory, that link quantum-statistical description taking Gauss distribution into account, explicit form is obtained by permutations group theory, the division paramagnetic/ferromagnetic is studied. The magnetization curves are built for a given set of parameters, such as exchange integral, number of closest neighbours and temperature. In a special range of the parameters a transition from unique solution of the resulting Heisenberg equation to multi-valued is observed. For this case exemplary hysteresis loops are built. The expression for Curie temperature allows to evaluate the exchange integral and proceed into temperature range above the critical temperature.
В рамках теории Гейзенберга, связывающей квантово-статистиче-ское описание с учетом распределения Гаусса, явная форма получается с помощью теории групп перестановок, изучается парамагнитное/ферромагнитное разделение. Кривые намагниченности построены для заданного набора параметров, таких как обменный интеграл, число ближайших соседей и температура. В специальной области параметров наблюдается переход от единственного решения полученного уравнения Гейзенберга к многозначному. Для этого случая построены примерные петли гистерезиса. Выражение для температуры Кюри позволяет учесть обменный интеграл и перейти в температурный диапазон выше критической температуры.
Keywords: hysteresis loop, ferromagnetism, paramagnetism, Curie temperature, Curie point, magnetisation curve, Curie — Weiss law, transcedental equation.
Ключевые слова: петля гистерезиса, ферромагнетизм, парамагнетизм, температура Кюри, точка Кюри, кривая намагничивания, закон Кюри — Вейсса, трансцедентное уравнение.
35
Introduction
The problem of theoretical description of magnetization takes the significant place in modern studies, beginning from seminal results of Weiss [2]. In the famous paper of Heisenberg [3] on ferromagnetism it was established that the Weiss electric forces are originated from exchange effect of quantum mechanics primarily introduced in textbooks' Heitler — London results. This paper [3] contains also very deep results of general significance. Multielectron terms theory was built using the very common symmetry: in respect to group of electrons permutations. It has its extension based on joint symmetry group of permutations and space symmetry group [9], where the exchange integral notion is «lifted» up the Hartree — Fock equations level. The permutation group theory allows to express energy via its irreducible
© Leble S. B., Chychkalo A. Yu., 2020
Вестник Балтийского федерального университета им. И. Канта. Сер.: Физико-математические и технические науки. 2020. № 3. С. 35 — 42.
representations characters [6] and, in same context, statistical distribution function is constructed. Its derivative gives the internal parameter (magnetisation M value conjugate to magnetic field H ,that results in equation of state M(H)).
The microscopic Heisenberg model goes to macroscopic Landau — Lif-shits — Gilbert (LLG) equations as in [7]. For example, in the article of Rivas [1] the ferromagnetic alloy Co66 Fe4 Mo2 Si16 B12is studied. They explore such LLG-based model for hysteresis loop building with the help of LLG differential equations of the second order. Experimentally, the alloy was isothermally annealed: the following material was heated, then slowly 36 cooled, in order to change the structural organisation of molecules in it. When the lattice (or other amorphous analog of it) changes, that changes the distance between units (atom or molecule), the properties of a material change too. The result is obtained at 530 degrees Celsium heating — there were some agglomerates of only 1 — 6 particles surrounded by the amorphous material. The group proceeded to the tests of its magnetic properties, changing the angle of applied magnetic field and other characteristics. Moreover, the bias of HL was detected.
In our paper we experienced two alike ways of graphical solutions of the resulting transcendent Heisenberg equations for the paramagnetism and ferromagnetism cases. It admits unique solution in the paramagnetism range, that gives the magnetisation curve. In the ferromagnetism domain few branches of solutions take place, that yields the loops.
Practically, notching the intersection points (IP) of functions into the so called matrix of values, which depends on the changes of external magnetic field strength and characteristics of the material. Connecting the extreme values of the matrix in case of multiple intersection between the same graphs, we managed either to build a hysteresis loop for a bulk ferromagnetic, or to identify some special cases, such as double loops with the central symmetry [10; 11].
In the last section we study a transition between para- and ferro-magne-tics, parametrized by Curie temperature, that generalizes the Curie — Weiss relation.
Para-magnetic materials. On theoretical base
The Heisenberg equations, derived in [3—5] result in a couple of explicit functions, depending on M [9]:
y1 = M, (1)
32 M „2 M3
a+pM-p1 — + p ~2
where
У = tanh(--2z-), (2)
ей H, (3)
mkBT
— = J, kBT'
(4)
z — a number of closest neighbours; e — an electron charge; h — Plank constant; m — an electron mass; kB — Boltzmann constant; T — temperature; J — an exchange integral.
The values of the parameters e, m, h are quite definite, while the z, J, T ones choice need some discussion. So, Heisenberg [3—5] estimates the parameter J by the note, that the exchange energy should be about kBT, giving the value J «10 13 erg . For such the value the parameter — is of order of 1 for the room temperature.
The paramagnetic case domain in zfi plane is fixed by the condition:
3(1 --) < 2.
z
(5)
37
The magnetisation curve construction
To build the magnetisation curve the values of M for each a are required. We take them by the means of solving the system of equations (1) and (2) graphically (Fig. 1). There is one intersection point between the line and the curve for para-magnetic materials, for the choice of z, —, such that (5) holds.
Fig. 1. The choice of parameters: — = 1; z = 10
The intersection points (IP) of a diagonal line with curves are the solution of the system of equations (1) and (2). The curves correspond to a, that is proportional to the magnetic field H, substitution from 0 value (the lowest curve) to 3 (the highest curve). The value of a between the closest curves differs in 0.1.
Next, the magnetization curve by the mentioned algorithm is built in Figure 2.
38
Fig. 2. Magnetisation curve for para-magnets looks typical [8]
The dependence of M on a in a range of a = 0..2, shown for the same parameter values p = 1; z = 10 as in Figure 1. With the increase of a value the curve changes like a linear. The thickening of intersection points is observed.
Ferromagnetic materials
In the case of ferromagnetic state the domain by (5) is above the curve
z(—-—), whence the solution of the Eq. 2 is not unique. In order to find out
kBT
the values of joint points of (1) and (2) graphs, we should follow the bends of hysteresis loop curves more accurately, using the alternative variant of graphical solution of the following equivalent equations:
y1 = arctan h( M ),
(6)
У 2 =
07l. „2 M 02 M3
a + PM -p — + p -
_z_
2 '
(7)
As it is observed in Figure 3 — the resultant loops, built by the means of (6) and (7) system, are rather convex and curved, comparing with (1) and (2) variant. The results of two different cases are presented in Figures 3—5.
Grey curves show the change of a value in 0.25 one by one. The black curve (6), unlike its possible alternative (1), gives the resultant thickening of
intersection points in fracture areas of the function. There are cases of 5 intersections of a single grey curve with the black — this is a critically needed condition to obtain the loop, illustrated at Figure 4.
Fig. 3. The intersection points of (6) and (7) are carried out for a = -0.75..0.75; — = 10; z = 8
Fig. 4. The hysteresis loop for common ferromagnetic materials. The values for its compilation are obtained from Figure 3. The presence of 5 intersection points (IP) described at Figure 3, gives the expansion of a loop at the approaching area
to abscissa axis
The further away from each other are the extreme points from the set of 5 IP (for a single grey tangensoid curve), the wider is the entire loop. The medium IP are not pictured and are not taken into account here. Though, medium IP are all inside the loop area.
The Heisenberg theory do not take into account domain walls presence in explicit form. It uses the Gauss distribution of states per energy level with
parameters expressed in terms of characters of the universal symmetry group, the group of permutations. Few interesting properties however are exhibited in some range of the parameters z and T, see again [10; 11].
The represented graph is the Hysteresis loop in case of maximum 3 intersections of (1) with the single curve (2) (a = -2.50..2.50;^ = 11; z = 8). There is no tangible difference if we take either ((6) and (7)) or ((1) and (2)) system for compilation. The vivid peculiarity is in 3 IP at once for a single curve ((2) and (7)) with the auxiliary functions from appropriate for each system of equations. The critical condition of double loop compilation: all three IP must be found in the same quarter area. The curves have a central symmetry for positive and negative a values, so the sub-loop at the third quarter (-a) always repeats the sub-loop at the first (+a) if a value is the same, having the difference only in its sign.
To add, the Ferromagnetism existence is demonstrated via the set of features, including the Curie law [8]. Its analog is derived at Heisenberg's paper [3] as well. According to it, there exists the critical temperature for switching from paramagnetic into ferromagnetic properties appearance and vice versa, by the boundary of the domain (5). It is marked as 0 at original text:
As well, the author emphasizes the similarity of the Curie Law with modification to the Weiss theory. Thus, the critical point is taken from the following conclusion, where M, T, z are mathematically connected.
40
Fig. 5. The double loop with the center symmetry of sub-loops
Towards the Curie law
(8)
As a result of following the Heisenberg's algorithm of transformations, we have obtained the resembling formula (9). We use it for graphical analysis and further substitutions of contingent constants, proceeding the imitations of practical situations with real materials.
M =-
4aT2
z(1 -,/1 - z )(T ~0)(T(1 + J1 - z )-0(1 -J 1 - z ))
(9)
Iron (Fe) was chosen as the primary sample of the corresponding section of this article out of its lattice structure. The approximation formula got vastly simplified (10) throughout the graphing procedure, because of the number of closest neighbours, which is 8 in our graphical imitation.
M =-
aT2
2(T -0)2
(10)
Fig. 6. The IP-printing for building (A) from Figure 7 with the usage of (1) and (2) system of equations. The horizontal axis shows the designation of magnetisation recession (z = 8;a = 0.1; 6 = 1043[K]; T = 1143..1903[K])
41
Fig. 7. The dependence of magnetisation on temperature (for Iron), where T >6. The comparison of two methods is graphed (z = 8;a = 0.1;6 = 1043[K ]; T = 1143..1703[K ]) (A) — the solutions of the equation system (1) and (2); (B) — the result plotted by (10)
42
Conclusions
We have built the magnetisation curves for paramagnetic and ferromagnetic ranges for temperatures below and above Curie temрerature, appropriate for chosen material. We would note that microscopic Heisenberg model has its continuous analog, named Landau — Lifshits — Gilbert equations. It could be used in making description of domain walls dynamics [12] and hysteresis phenomenon as well [1].
References
1. Martinez-Garzia J. C., Rivas M, Garcia J. A. Induced ferro-ferromagnetic exchange bias in nanocrystalline systems // Journal of magnetism and magnetic materials. 2015. Vol. 377. P. 424—229.
2. Weiss P. Molecular field and ferromagnetic property // Journ. de phys. 1907. Vol. (4) 6. P. 661.
3. Heisenberg W. Zur Theorie des Ferromagnetismus // Zs. Phys. 1928. Vol. 49. P. 619—636.
4. Heisenberg W. Zur Quantentheorie des Ferromagnetismus // Probleme der Modernen Physik, A. Sommerfeld Festschrift. Leipzig, 1928. P. 114—122.
5. Heisenberg W. Zur Theorie des Magnetpstriktin und der Magnetisierungkurve // Zs. Phys. 1931. Vol. 69. P. 287—297.
6. Heitler W. Staorungsenergie und Austausch beim Mehrkaorperproblem // Zs. Phys. 1927. Vol. 46. P. 47.
7. Lakshmanan M. The fascinating world of the Landau — Lifshitz — Gilbert equation: an overview // Phil. Trans. R. Soc. A. 2011. Vol. 369. P. 1280 — 1300.
8. Kittel C. Introduction to Solid State Physics. Wiley, 2004.
9. Leble S. Heisenberg chain equations in terms of Fockian covariance with electric field account and multiferroics in nanoscale // Nanosystems: physics, chemistry, mathematics. 2019. Vol. 10 (1). P. 18 — 30.
10. Hirsch A. Double hysteresis loops in ferromagnetic crystals // J. Phys. Radium. 1959. Vol. 20 (2—3). P. 262— 263.
11. Brandao J., Dugato D.A., Seeger R.L. et al. Observation of magnetic skyrmions in unpatterned symmetric multilayers at room temperature and zero magnetic field // Sci Rep. 2009. Vol. 9. 4144. doi:10.1038/s41598-019-40705-4.
12. Leble S. Waveguide Propagation of Nonlinear Waves. Impact of Inhomoge-neity and Accompanying Effects. Springer, 2019.
The authors
Prof. Sergey B. Leble, Immanuel Kant Baltic Federal University, Russia.
E-mail: lebleu@mail.ru
Anastasiia Yu. Chychkalo, Master's Student, Immanuel Kant Baltic Federal University, Russia.
E-mail: nanacappella23@gmail.com
Об авторах
Сергей Борисович Лебле — д-р физ.-мат. наук, проф., Балтийский федеральный университет им. И. Канта, Россия.
E-mail: lebleu@mail.ru
Анастасия Юрьевна Чичкало — магистрант, Балтийский федеральный университет им. И. Канта, Россия.
E-mail: nanacappella23@gmail.com
