Revisiting the holstein-primakov transformations

It is shown that in addition to the Holstein-Primakov transformations in the theory of magnetism, a number of other transformations can be proposed which also lead to very interesting and consistent results. In particular, with the help of the transformation proposed in the paper for spin operators, it turns out to be possible to strictly analytically calculate the temperature dependence of the ferromagnet magnetization in a wide temperature range from zero and up to the Curie temperature. It is shown that for these transformations in the Curie temperature area, the magnetization tends to zero exponentially. Thus, the aim of this study is to describe magnetization over a wide range of temperatures using a new transformation for electron spin operators


Introduction
The issue that we would like to focus on in this paper relates to the classical theory of magnetism.The essence of the problem which will be discussed below, is to clarify some types of transformations for operators of atomic spins in magnetic materials.The fact is that the well-known Holstein-Primakov transformations give a correct answer to the issue of the temperature dependence of magnetization in the low-temperature area which however does not mean their uniqueness at all.In fact, if we calculate the commutator of spin operators S + , S − for spin S = (S x , S y , S z ) S ± = S x ± iS y it is defi ned as [S − , S + ] = 2S z [1,2].Therefore, any set of transformations makes it possible to fi nd operators S z by the type of operators S + , S − .In this case, dependence S z can be any, but in the limiting case of large spins, it should lead to equality S z = S − a + a where a + (a) − is the operator of the birth (destruction) of the spin deviation (magnon).As a result of the averaged spin component, S z  we arrive at Bloch's law in the low-temperature area (the socalled «three-second law» for a ferromagnet [3,4].
In general, the form of Hermitian-conjugate operators S + , S − can be represented as where n = a + a and f(n) − is an arbitrary operator function.
According to the switching rule, the longitudinal spin operator therefore is defi ned as In the special case of the Holstein-Primakov transformations the function is , and after calculating the corresponding commutators we come to an exact expression for the longitudinal spin operator S z = S − a + a.
It is quite clear that such dependence f(n) is hardly the only and we can also consider some other quite possible cases one of which, in our opinion, should be discussed in more detail.

Non-root dependency of the operator f(n)
Citation: Then due to the switching equality [1] , Taking into account the defi nitions (1) after a simple calculation we have the following from here As a result for z − spin operator components we arrive at the following relation: From this, it can be seen that in the limiting case when the condition 2S >> 1 is met expression (4) is converted to the form According to the ratio (5) we also come to Bloch's «threehalves» law which does not contradict the Holstein-Primakov theory as well as a lot of experimental facts confi rming it.
However, as can be seen from ( 5) we can also calculate additional temperature corrections to the «three-halves» law.They are determined by the last summand standing on the right side of expression (5) but which is not the last in the literal sense since it is a consequence of the exponential expansion if we proceed from formula (4).
This fact can be considered one of the advantages of transformations (1).Another advantage in our opinion is that these transformations allow us to calculate the temperature behavior of magnetization not only in the low-temperature area but also in the high-temperature area in particular in the Curie temperature area.In order to prove this statement, we will focus on this issue in more detail.

Temperature behavior of magnetization
As is well known (see for example [3] as well as the monograph [5]) the average value of magnetization is obtained by averaging the observed magnitude according to the equilibrium density matrix which was fi rst introduced by Bloch and Landau.It is on the basis of this defi nition that the average value of magnetization can be easily calculated.
Indeed, we follow the general thermodynamic defi nition of the average magnetic moment of an atom then in accordance with formula (4) it should be represented as Where angle brackets mean averaging by the basic state of a ferromagnet with an equilibrium density matrix.
To introduce the magnon energy  (k) where k − is its wave vector we should write down a general expression for the energy of a ferromagnet placed in a constant and uniform magnetic fi eld H = (0,0, H) with an axis z selected along the ferromagnet anisotropy axis direction.
As an illustrative example consider a easy axial ferromagnet with cubic symmetry.The general expression for its magnetic energy can be represented then as the following Hamiltonian (see for example [6] as well as review [7]): Where  − the characteristic energy of magnetic anisotropy In a + a operator-quadratic approximation from (7) we have H 2 2 Or in a compact form Where and the energy operator 2 The transition in (10) to k − representation using the expansion of the exchange integral and the birth and destruction operators into a Fourier series allows us to rewrite operator (10) in a compact form as Where the magnon energy is Therefore taking into account (12) returning to formula (6) in which the averaging can now be concretized using a transition to k − space for averaged thermodynamic values we have the right to write that Where parameters are , 1 .
Here and everywhere below the Boltzmann constant k B is considered equal to one.
As a result, taking into account the ( 13) -( 15) formula ( 6) can be rewritten as follows .
Or after simplifi cation we have If S >> 1 then it follows from formula (16) in linear approximation by S that   Where .
To analyze the general formula ( 16) it is convenient to rewrite it in non-dimensional form taking into account explicit expressions for parameters  and  according to (15) as well as the magnon dispersion law (12).Since in the simplest case dependence J(k) can be approximated by a periodic function Where a − is the interatomic distance it follows from (12) that in the long-wave approximation, the magnon dispersion has a well-known quadratic form Therefore for the temperature area from the range It can be seen from the presented fi gures, in particular, that the proposed transformation for spin operators makes it possible to describe the dependence of the average spontaneous magnetization of a ferromagnet in a fairly large range of temperature changes almost to the Curie temperature.It can also be seen from these fi gures that for each spin value, a ferromagnetic material must have its own Curie temperature, which is quite natural.
To relate Curie's temperature T C to exchange energy J ex , which according to (18) is defi ned as Let us proceed as follows.
We will fi nd the asymptotic behavior of function y(x) according to the obtained dependence (19), and we will «link» them.
Indeed if x  1 then The Curie temperature is in the area T T 0 and its value can be estimated from Figure 1 which is shown by a dot with a designation T c along the abscissa axis.
If we use for example the classical monograph [1] then following it the Curie temperature is associated with the exchange integral by ratio Where Z − is the number of nearest neighbors?If for example a cubic crystal lattice is volume-centered then Z = 8, and if it is simply cubic then Z = 6.
From the comparison of formulas ( 23) and ( 24) it can be seen that ratio must be observed which is well performed for almost any value of spin S.
The temperature dependence for the average magnetization in the region of high temperatures as is well known (see [1]) is described by Curie's law, i.e.According to the Brillouin function [5].
As can be seen from the above fi gure, at high temperatures, according to dependence (22) the behavior of the average magnetization is fully correlated with the classical law (25).
The disadvantage of the proposed transformation is apparently the qualitative difference between the hyperbolic dependence (22) and the root dependence which is characteristic in close proximity to the Curie point.According to (19) in temperature area T > T c due to the asymptotics of (22), it is quite possible to assume that <S z > ≈ 0.

Conclusion
Concluding this message, the following is worth noting.
If we use the transformation of electron spin operators which differs from the Holstein-Primakov transformations, we are fi nding Bloch's «three halves» law in the low-temperature area.And in the fi eld of high temperatures, we come to a slightly refi ned Curie's law. Figure 1 shows the dependence of the average magnetization over the entire temperature region.

Figure 1 : 5 2S
Figure 1: Line 1 corresponds to the value of an infi nitely large spin, i.e. S 1. Curve 2 corresponds to the value of the spin 5 2 S  , and curve 1 corresponds to the value