Two Dimensional Heisenberg Exchange Interaction in the Magnetization Studies of Multiferroic ()
1. Introduction
Multiferroics are unique classes of materials whose fundamental state is both magnetic and ferroelectric [1]. In these materials magnetism and ferroelectric ordering coexist no matter how their course of existence varies [2]. The microscopic description of magnetism is based on the localized magnetic moments, mostly due to the partially filled d or f shell transition-metal or rare-earth elements electrons in which the exchange interactions between these localized moments exhibit magnetic orderings. Whereas the origin of ferroelectrics is wide range (can be due to: Charge ordering, lone pairs, etc.) and forbids the d shell electrons in contrast to magnetism [3,4]. Thus, multiferroics couples the two insurmountable phenomena. They remarkably show strong coupling between charges, spins, lattice and orbital degrees of freedom and the interplay between various competing orders gives rise to interesting features and varieties to these classes of materials. The overlap of strongly correlated magnetic and ferroelectric orders enables the spinnability of one domain over the other [5-7]. For instance, in ferromagnetic materials the magnetization M is easily switchable parameter on an application of external stimulating magnetic field, H and in ferroelectrics the polarization P is also a controllable parameter by an external electric field, E. However, on the other hand, in multiferroics these switchable order parameters, P and M are cross-linked, i.e. magnetization can be switched by external electric field and polarization can also be tunned by an applied magnetic field. The magnetoelectric spectra in GaFeO3 is found to be proportional to the photon energy which depends on the direction of magnetization [8-10].
2. Theoretical Model
We consider the spin Hamiltonian of the multiferroic thin film system that exhibits frustrated spin interactions [11,12]. The Hamiltonian of this anisotropic many body interaction system on the influence of external fields can be expressed by [13-15],
(1)
where Si(Sj) is spin operator at sites i(j). Jij is the exchange coupling that depends on the relative positions of neigbhoring spins and is defined by Jij = J(Ri − Rj). The constant g is the gyromagnetic ratio, and the second term of Equation (1) is the effective Zeeman field with an external deriving field while in this model chosen to be parallel to the +z axis to line up all the spins in the same direction. Whereas, the third term includes an electric field E that interacts with the spins in the system. The constant μ describes the coupling strength between the spins and the radiation field that is either intrinsic or extrinsic [16,17]. The terms μB and D are the Bohr magneton and magnetic anisotropic (or spin wave stiffness) constants, respectively [17,18]. The total spin operators in the Hamiltonian are:
The ground state of N identical atoms of spin S is defined as
The total spin operator at any site i has components Siz, Six and Siz that should be treated independently with an identity Si. Si = S(S + 1). Transferring the spin operator problems into many body interaction systems upon replacing with bosonic creation and annihilation operators expressed as
(2)
(3)
where ai† and ai are the creation and annihilation operators satisfying the commutation relation [ai, ai†] = δil It is also possible to transform this bosonic operators to magnon variable operators bk and bk† so that these operators can be related to the bosonic operators [19,20]
and
where the inverse transformation is
and
At low excitation i.e. where the interaction is predominated by bilinear spin variables, neglecting the interaction corresponding to higher order terms, the reduced Hamiltonian is of the form
. (4)
Here, z is defined to be the number of nearest neighbor spins. In this Hamiltonian we consider only the term bilinear to the magnon variables
(5)
where,.
This simplified interaction term is equivalently expressed as, if we write
(6)
The dispersion equation on the approximation that can be reduced to
, (7)
so that
(8)
The above expression is a dispersion relation with some additional term that tends to show how the dispersion term depends on the field components. Moreover, the constant.
3. Results and Discussion
The spin-wave dispersion function of the multiferroic thin film system is obtained as a function of the wave vector and the field components.
A significant change in the dispersion relation is observed as depicted by Figure 1 above at different γ values [21]. This can be used to modify the behavior of the multiferroics up on appropriate selection of the magnetic
Figure 1. (color online) The figure indicated shows the variation of dispersion function for different gamma values.
and electric field values in some orientations. The inter-dependence of the magnetic and electric fields (their couplings) in multiferroic materials is also analyzed based on the spin-wave occupation number,.
The change in magnetization is obtained to be
where and and the change in magnetization versus temperature at some specific γ is shown in Figure 2 below for low temperature limits.
Consequently, the magnetization is a function of the γ values at low temperature as indicated by Figure 3 shown below [22,23]. Thus, the electric field dependence of
Figure 2. Change of magnetization versus temperature at different gamma values.
Figure 3. Dependence of magnetization on the parameter γ, which can be tuned by explicitly on the applied fields and the nature of magnetic anisotropy.
magnetization can be easily deduced from the expression γ = ζ + μE where ζ = −2DS + gμB H for some specific values of the parameters, D, S, and H at a particular temperature.
4. Acknowledgements
This work was supported by the research programm of Aaddis Ababa University, Ethiopia and Mizan-Tepi University, Ethiopia.