CC BY
8
УДК 537.61
Continual Approach in the Mean Field Theory of Incommensurate Magnetic States in the Frustrated Heisenberg Ferromagnet
Sergey N. Martynov Vasiliy I. Tugarinov*
Kirensky Institute of Physics, SB RAS Akademgorodok, 50/38, Krasnoyarsk, 660036
Russia
Aleksandr S. Martynov
Faculty of Physics Moscow State University Moscow, 119991
Russia
Received 20.08.2016, received in revised form 10.10.2016, accepted 15.11.2016 The algorithm of approximate solution has been developed for the differential equation describing the anharmonical change of the spin orientation angle in the model of ferromagnet with the exchange competition between nearest and next nearest magnetic neighbors and the easy axis exchange anisotropy. The obvious dependence of the angle velocity from angle and initial condition has been derived by expanding the first integral of the equation in the Taylor series in vicinity of initial condition.
Keywords: incommensurate magnetic structures, soliton lattice, continual approach. DOI: 10.17516/1997-1397-2017-10-1-55-59.
The theoretical description of incommensurate magnetic structures (IMS) in antiferromag-netic dielectrics in the framework of Landau phenomenological theory of the phase transition was elaborated by Dzyaloshinskii [1]. For helimagnet with the Dzialoshinskii-Moria antisymmetrical exchange leading to a Lifshitz invariant in a free energy expansion (relativistic mechanism of forming IMS) the energy minimization for the solutions within the anisotropic plane is reduced to the solution of static sine-Gordon equation. The equation has the anharmonic solutions in the form of elliptical integrals describing the inhomogeneous helical structure with the step changing when moving along the helix vector - the so called soliton lattice [2]. For helimagnets with competing exchange interactions between the nearest and further magnetic neighbours (the exchange mechanism) it is necessary to take into consideration the second derivatives (and highter ones in the general case) of an order parameter and the energy minimization can not be reduced to an analiticaly integrable differential equation. The aim of the present work is to develop the algorithm of approximate solution for the equation describing the anharmonic change of the spin orientation angle in the easy axis ferromagnet with the exchange competition between the nearest and next nearest neighbors. We consider the magnetic structure with the spin orientation within the anisotropy plane (the flat anharmonic helix) and depends on the one coordinate.
* [email protected] © Siberian Federal University. All rights reserved
The Hamiltonian of the classical spins S=1 with the ferromagnetic and antiferromagnetic exchanges between nearest and next nearest neighbors accordingly and the easy axis exchange anisotropy (XXZ-model) has a form
H = J^(SjSj+i + SSZS+) + (SiSi+2 + SSZS+), Ji < 0, J2 > 0, S> 0. (1) i i
The same relative anisotropy in the both exchanges allows to reduce the number of model parameters and simply separate an effect of frustration and anisotropy on the IMS energy. The orientation of the spin Si in the plane with the easy axis z is determined by the total exchange field from neighbor spins Si±i and Si±2. At T=0 all spins have the equal length equal to saturation one. The local field on the site i normalized on the exchange Ji has the components
hz = (1 + 6)h°, 1
h° = 2 (cos 0i+i + cos 0i_i + R(cos Oi+2 + cos #¿-2)), (2)
hx = ^(sin #¿+1 +sin 6i_i + R(sin Oi+2 +sin 0_2)),
where h0z is the z-component of exchange field without anisotropy, 9 is a polar angle and R = J2/Ji < 0 is a frustration parameter. The transition to the continual description is carried out by the Taylor series expansion of the neighbor spins angles on the each site
di = 9, #i±i = 9 ± En + Si2, #¿±2 = # ± S2i + S22, (3)
where Eap are the sums of odd and even derivatives of the variable 9
~ #(2n_i) ~ #(2n)
Eii ^(2n - 1)!' Ei2
n=i n=i (4)
22n_i#(2n_i) 22n#(2n)
E2i = E (2n - 1)! ' E22 =
n=i v ^ n=i v ^
After substitution (3) the components (2) and longitudinal field on the spin (an energy density on the single interval in coordinate space) takes the forms
h0 = cos(9 + E i2) cos Eii + R(cos(9 + E22) cos E2i),
hx = sin(9 + Ei2) cos Eii + R(sin(9 + E22) cos E2i). (5)
h|| = e = hz cos 9 + hx sin 9 = cos 9 + e0,
where e° = cosEii cosEi2 + RcosE2i cosE22 is an energy in the isotropic case 6 = 0. An orientation of each spin is uniquely determined by the collinearity condition of spins and local fields from the neighbor spins [3]. The transverse field on the spin must be equal to zero
h± = hz sin 9 — hx cos 9 = 6h° sin 9 — A° = 0,
(6)
A° = cos Eii sin Ei2 + Rcos E2i sin E22.
Taking into account the collinearity constraint (6) the magnetic energy density takes a multiplicative form - the anisotropic and frustration components hold as product terms
1+6
• e°.
1 + 6 sin2 9
Rewriting the general equation (7) in the following form
cos En sin £12 + R cos £21 sin £22 6 sin 0 cos 0 cos £11 cos £12 + R cos £21 cos £22 1 + 6 sin2 0'
(7)
one can make a general conclusion that anharmonicity in the change of angle 0 (the derivatives of the second order and higher) is appeared at 6 > 0, takes a maximum value at 0 = (2n + 1)n/4 and vanish at 0 = nn/2.
Further solution of equation (6) will be carried out in the linear anharmonicity approximation (0" << 1), neglecting the derivatives which are higher than second order. In this approach the Eq. (6) takes an autonomous form
0' cos 0' + 4R cos 20' 26 sin 0 cos 0 ' cos 0' + R cos 20' = 1 + 6 sin2 0 '
and at the substitution z = 0'2/2 is integrated by quadratures
I(z,z0)= [ ^^dz = ln(1 + 6sin2 0), (9)
Jz 0 £0
where C(z) = cos%/2z + 4Rcos2%/2z, £0(z) = cos%/2z + Rcos2%/2z. The variable z range {z0, zmax} in the interval 0 G {0, n/2}. Expanding the integral in the Taylor series in vicinity of z0 and taking into account that a0 = I(z = z0) = 0 we obtain the series
I( ) ^ (z - zp)" (C(z) xCn-1)
I (z,z0)=y an----, an = --. (10)
^ n! V £0 ) zo
"=1
At linear anharmonicity keeping first two nonzero terms of series (10) (the quadratic approximation) we obtain the obvious dependence of the angle velocity from angle and initial condition z0
£0(z)
z = z0 - -
3R • K(z0)(1 + 2cos2 ,/2z0)
x (\JC2(z0) - 6R • K(z0)(1 + 2 cos2 ■ ln(1 + 6 sin2 0) - C(z0^, (11)
where K(z0) = sin^2z0)^2z0.
The energy of the quarter of period divided by the corresponding length is equal to the average energy of one spin
(1+6) r
0
/2 £0(z) d0
E(z 1 +6sin2 0 ^/Me^ E(z0) =-F/-2-d0-. (12)
0
V/2z(0, z0)
The ground state of the solutions is determined by an extreme of the function (12) (the maximum in our case, because the fields and energy density are normalized on the J1 < 0). The corresponding initial condition z0 = zextr parametrizes the ground state solution for each set of the parameters 6 and R and after substituting into Eq. (11) determines the functional dependence of z on the angle 0 in the ground state. The solutions with increasing helix pitch upon changing the angle between the spin and easy axis 0 from zero to n/2 (0" > 0) are exist
when the condition z0 > zmin = (1 + 4R)/(1 + 16R) > 0 is fulfilled. It is impossible for the IMS wave vector to tend to zero at the phase transition from incommensurate to commensurate state, which also follows from the phenomenological analysis based on the Landau theory [2]. It is means that the transition on the model parameters S and R is accompanied by step-like change of the magnetic structure vector and hence is a first order phase transition.
The range of the z variation from zextr to zmax is increased with increasing S and takes a maximal value on the phase boundary with collinear phases. The difference zmax — zextr takes the maximum value equal to 0.4 in the triple point R = 1/2, S « 1.34 where the energy of the soliton phase is equal to the energies of ferromagnetic and "up-up-down-down" phases.
To assess the application of the quadratic approximation at the expansion of the first integral (10) one makes a numerical integration in the triple point and compares the result with the quadratic and cubic decomposition (10) (Fig. 1).
I(z)
lii( l ôMiiij)
0,8
0,4
0,0
ln[l+ôsin2(7i/2)] 4/
R=—1/2, 6=1.34 ln[l+6sin2(7t/4)] / ' À 'A / / /y // /V ''3
A/ ■ Kr / ■ ■ ' Z max z
0,4
0,6
0,8
0,4
0,0
0,8
Fig. 1. The numerical integration of the first integral (9) (solid line 1) and the quadratic (dash line 2) and cubic (dash-dot line 3) approximations are shown with the right part of Eq. (9)
The coefficients a1 and a2 are positive. The third coefficient a3 is negative and together with the next terms of expansion forms an alternate series which provides a fast expansion convergence to the numerical integration result.
References
[1] I.E.Dzyaloshinskii, Theory of helical structures in antiferromagnets. 1. Nonmetals, Sov. Phys. JETP, 19(1964), 960.
[2] Y.A.Izyumov, Neitronnaya diffraktsiya v dlinno-periodicheskih structurah (Neutron Diffraction in Long-Periodic Structures), Energoatomizdat, Moscow, 1987 (in Russian).
[3] S.N.Martynov, Flat and conical incommensurate magnetic structures in the two-subsystem partially frustrated heisenberg ferrimagnet, JMMM, 398(2016), 121.
Континуальное приближение в теории среднего поля несоизмеримых магнитных состояний
во фрустрированном гейзенберговском ферромагнетике
Сергей Н. Мартынов Василий И. Тугаринов
Институт физики им. Л. В. Киренского СО РАН Академгородок, 50/38, Красноярск, 660036
Александр С. Мартынов
Физический факультет МГУ имени М.В.Ломоносова Ленинские горы, 1/2, Москва, ГСП-1, 119991
Россия
Разработан алгоритм приближенного решения дифференциального уравнения, описывающего ангармоническое изменение угла ориентации спина в модели легкоосного ферромагнетика с конкуренцией обменов между ближайшими и следующими за ближайшими магнитными соседями. Явная зависимость угловой скорости от угла и начального условия получена путем 'разложения первого интеграла уравнения в ряд Тейлора в окрестности начального условия.
Ключевые слова: несоизмеримые магнитные структуры, солитонная решетка, континуальное приближение.
