# Semiclassical Modeling of Isotropic Non-Heisenberg Magnets for SpinS=1 and Linear Quadrupole Excitation Dynamics

Authors: Yousef Yousefi; Khikmat Kh. Muminov
Journal: Physics Research International
 (2013)

---

## Abstract

Equations describing one-dimensional non-Heisenberg model are studied by use of generalized coherent states in real parameterization, and then dissipative spin wave equation for dipole and quadrupole branches is obtained if there is a small linear excitation
from the ground state. Finally, it is shown that for such exchange-isotropy Hamiltonians, optical branch of spin wave is nondissipative.


---

## Body

## 1. Introduction

Many condensed matter systems are fully described by use of effective continuum field models. Topologically nontrivial field configurations have an important role in modeling systems with reduced spatial dimensionality [1]. Magnetic systems are usually modeled with the help of the Heisenberg exchange interaction [2–4].However, for spinS>1/2, the general isotropic exchange goes beyond the purely Heisenberg interaction bilinear in spin operators S→i and includes higher order terms of the type (S→iS→j)n with n up to 2S [5]. Due to the spin states, the 2S+1 complex parameters are necessary to describe each of them, and this corresponds with the 4S+2 degrees of freedom. Two degrees of freedom are omitted, one because of normalization condition and the other for arbitrary phase decrease, hence 4S parameters are required to completely modeled the remainder 4S degrees of freedom of spin states [6].Particularly, in caseS=1 with the isotropic nearest neighbor exchange on a lattice is derived by use of the Hamiltonian
(1)H^=-∑i‍J(S^→iS^→i+1)+K(S^→iS^→i+1)2.HereS^ix,S^iy,S^iz are the spin operators acting at a site i, J and K are, respectively, the bilinear (Heisenberg) and biquadratic exchange integrals. The model (1) has been discussed recently in connection with S=1 bosonic gases in optical lattices [7] and in the context of the deconfined quantum criticality [8, 9]. Hamiltonian (1) is a special form presented in [10] and because of importance of quadrupole excitation in ferromagnetic Materials, it is considered here. This paper does not consider the antiferromagnetic and nematic states.Considering the effects of both dipole and quadrupole branches gives a nonlinear approximation. If higher order multipole effects are considered, the approximation is more accurate but at the same time, deriving the equations is too complicated. In this paper, only the effect of quadrupole branch for Hamiltonians described by (1) is considered. Study of isotropic and anisotropic spin Hamiltonian with non-Heisenberg terms is complicated due to quadrupole excitation dynamics [5, 11, 12]. Antiferromagnetic property of this excitation in states near the ground proves the existence of it, and Dzyaloshinskii calculated the effect of this excitation [13]. Also, numerical calculations more accurately justify laboratory results if the effect of quadrupole excitation in nanoparticles Fe8 and Mn12 is considered [14, 15]. In addition, this method may be promising for description of the multispin configuration of the Fe in the different ligand coordinations [16].To calculate effect of quadrupole excitation, at first it is necessary to obtain classical equivalent of Hamiltonian (1), and then to find out the solution of spin wave, it is necessary to analyze resulted equations in case there is small linear excitation from the ground states. Therefore, the stages of process are the following.(1)
Obtaining coherent states for spins=1 which are coherent states of SU(3) group.(2)
Calculation of average values of spin operator.(3)
Classical spin Hamiltonian equation is obtained using previously calculated values.(4)
Calculating Lagrangian equation by use of Feynman path integral over coherent states and then computing classical equations of motion.(5)
It is necessary to substitute resulted Hamiltonian in classical equations of motion to obtain nonlinear equations of magnets. Solutions of these nonlinear equations result in soliton description of magnet that is not interested here.(6)
Now, ground states of magnets calculated and then nonlinear equations are linearized above the ground states for small linear excitation.(7)
Finally, spin wave equations and dispersion equations must be calculated.In this paper, we write coherent states in real parameters, because each parameter in this representation related to one degree of freedom. In complex parameterization, each parameter related to two or more degrees of freedom. Then in physical problems, the first representation is very helpful.In next sections, we develop mathematical descriptions of the pervious stages and analyze these descriptions.


## 2. Theory and Calculations

In quantum mechanics, coherent states are special kinds of quantum states that their dynamics are very close to their corresponding classical system. Type of coherent state is used in problem and depends on the operators symmetry in Hamiltonian. Due to the operators symmetry in Hamiltonian (1), for more detailed description and considering all multipolar excitation of coherent states we used coherent states in SU(3) group. Vacuum state of this group is (1,0,0)T and coherent state for a single site in this group is introduced as [11]:
(2)|ψ〉=D1(θ,ϕ)e-iγS^ze2igQ^xy|0〉,
where D1(θ,ϕ) is wigner function and Qxy is quadrupole moment. Two angles, θ and ϕ, the Euler angles, determine the direction of classical spin vector in spherical coordinate system. The angle γ determines the direction of quadrupole moment around the spin vector and parameter g shows change of magnitude of spin vector and average value of quadrupole moment. Two angles ϕ and γ change between 0 to π and angle θ changes between -π to π.In order to derive Lagrangian from path integral, we drive the path integral from the following transition amplitude:(3)P(ψ1,t1;ψ,t)=〈ψ1|exp(-ih-(H^(t1-t)))|ψ〉.
Using completeness relation and doing some mathematical work, the following equation for transition amplitude obtained:
(4)P(ψ1,t1;ψ,t)=limn→∞∑j‍∫‍Dμj(ψ)exp(-ih-∫tt1‍Lj(θ,ϕ,g,γ)dτ).In the pervious equationL is Lagrangian and has the following form [11]:
(5)L=h-cos2g(cosθϕt+γt)-H(θ,ϕ,g,γ).x
t
=
(
∂
x
/
∂
t
)
(
x
=
θ
,
ϕ
) and H(θ,ϕ,g,γ) is classical Hamiltonian. By use of (4) and (5) and the action stationary principle, the classical equations of motion are obtained.When Lagrangian is obtained from path integral, another two terms appear, one is kinetic term and has “Berry phase” properties that are important in spin tunneling phenomena and the other depends on boundary condition values. In this paper, we do not consider these two terms.Now classical equivalent of spin vector and their products must be computed so that the classical equivalent of Hamiltonian (1) is obtained. So consider the following:
(6)S→=〈ψ|S^→|ψ〉
as classical spin vector, and also consider the following:
(7)Q^ij=12(S^iS^j+S^jS^i-43δijI)
components of quadrupole moment. Because we can write any coherent state as product of single site coherent states, namely,
(8)|ψ〉=∏i‍|ψ〉i.Then spin operators in ground state of nonsingle ions Hamiltonian can be commuted in different lattices [11], so
(9)〈ψ|S^niS^n+1j|ψ〉=〈ψ|S^ni|ψ〉〈ψ|S^n+1i|ψ〉,
where |ψ〉=|ψ〉n|ψ〉n+1.The average spin values in SU(3) group for coherent states (2) are defined as [17]
(10)S+=eiϕcos(2g)sinθ,S-=e-iϕcos(2g)sinθ,Sz=cos(2g)cosθ.
And also
(11)S2=cos2(2g),q2=sin2(2g),S2+q2=1.In the pervious relation,S2 is related to dipole moment and q2 is related to quadrupole moment. Classical Hamiltonian can be obtained from the average calculation of Hamiltonian (1) over coherent states. The classical continuous limit of Hamiltonian in SU(3) group is
(12)H=-∫‍dxa0(Jcos2(2g)+Kcos4(2g)-a022(4gx2sin2(2g)(J+2Kcos2(2g))+cos2(2g)(J+2Kcos2(2g))×(θx2+ϕx2sin2θ))),
where a0 is length of crystal sites. The pervious classical Hamiltonian is substituted in equation of motion that obtained from the Lagrangian, and the result is classical equations of motion:
(13)1ω0θt=-a02cos(2g)(J+K+Kcos(4g))ϕxxsinθ,1ω0ϕt=a02cos(2g)(J+K+Kcos(4g))(ϕx2cosθ+θxxcscθ),1ω0gt=0,1ω0γt=4cos(2g)(J+K+Kcos(4g))+({12}Kcos3(2g)(16gx2-8θx2-5ϕx2+3ϕx2cos(2θ)-2θxx)+cos(2g)×({12}8gx2(J-K)-2Jθx2+8gx2Kcos(4g)+12Jϕx2(-3+cos(2θ))-Jθxxcotθ)+4gxxJsin(2g)+2gxxK(sin(2g)+sin(6g)){12})a02.In these equations,ω0 is h-a0. These equations describe nonlinear dynamics of non-Heisenberg ferromagnetic chain completely. If we omitted quadrupole excitation (g=0) in the pervious equations, these equations reduced to Landau-Lifshitz equations. Then in comparison with Landau-Lifshitz equations, these equations are more complete and contain more degrees of freedom. Note that solutions of these equations are different forms of magnetic solitons.In this paper, only the linearized form of (13) for small excitation above the ground states is considered. To this end,at first, classical ground states must be calculated, so in above Hamiltonian only nonderivative part is considered:
(14)H0=-∫‍dxa0(Jcos22g+Kcos42g).To find the smallest value of theH0 we vary it respect to all the parameters, with the ground state is obtained at;
(15)g=0,g=π2.In this paper, only dispersion of spin wave in neighborhood of the ground states is studied. For this purpose, small linear excitations from the ground states, as shown in (15), are defined:
(16)2g→π+g.In this situation, the linearized classical equations motion are(17)1ω0θt=-a02(J+2K)ϕxx,1ω0ϕt=a02(J+2K)θxx,1ω0gt=0,1ω0γt=-4(J+2K).Consider functionsθ and ϕ as plane waves to obtain dispersion equation:
(18)ϕ=ϕ0ei(ωt-kx)+ϕ-0e-i(ωt-kx),θ=θ0ei(ωt-kx)+θ-0e-i(ωt-kx).Substitution of these equations in (17) results in dispersion equation of spin wave near the ground states.
(19)ω12=(J+2K)2a02k4ω02,ω2=-4(J+2K)ω0.It is evident from (19) that in addition to the dispersion acoustic branch, there exist nondispersion optical branches which is related to the dipole and quadrupole excitations.


## 3. Conclusion

In this paper, describing equations of one-dimensional isotropic non-Heisenberg Hamiltonians are obtained using real-parameter coherent states. It is shown that both dipole and quadrupole excitations have different dispersion if there is small linear excitation from the ground state.In addition, it is shown that for isotropic ferromagnets, the magnitude of average quadrupole moment is constant (gt=0) and its dynamic, is rotational dynamics around the classical spin vector (γt≠0).



---

*Source: 634073-2013-03-27.xml*