Two Dimensional Heisenberg Exchange Interaction in the Magnetization Studies of Multiferroic

Multiferroics are novel classes of materials that exhibit cross-coupling of mutually excluding phenomena, i.e. magnetism and ferroelectricity. In recent years, the coexistence of ferroelectricity and magnetic orderings has become a hot issue and drawn considerable attentions due to the promising applications to these days technology and the fundamental science involved in these classes of materials. The microscopic origins of magnetism and ferroelectricity differ fundamentally, while the real mechanism of ferroelectricity is still under debate. In the present work, we have started from a simple method Heisengerg hamiltonian and an interaction term resulting from electric field coupling with the magnetic spins with anisotropic limit, demonstrated that magnetization can be manipulated by electric field and anisotropic field in agreement with results experimentally observed. In the multiferroic thin film system the magnetic field tends to play a role in stabilizing the spins in preferred orientations and induces a coupling of magnetism and ferroelectricity that opens a route to switch magnetization with electric polarazation and vice-versa.


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][6][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 ex-ternal 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 GaFeO 3 is found to be proportional to the photon energy which depends on the direction of magnetization [8][9][10].

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][14][15], )   where S i (S j ) is spin operator at sites i(j).J ij is the exchange coupling that depends on the relative positions of neigbhoring spins and is defined by . 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 Two Dimensional Heisenberg Exchange Interaction in the Magnetization Studies of Multiferroic 216 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 S iz , S ix and S iz that should be treated independently with an identity S i .S i = S(S + 1).Transferring the spin operator problems into many body interaction systems upon replacing with bosonic creation and annihilation operators expressed as

2S
1 2S where a i † and a i are the creation and annihilation operators satisfying the commutation relation [a i , a i † ] = δ il It is also possible to transform this bosonic operators to magnon variable operators b k and b k † so that these operators can be related to the bosonic operators [19,20] where the inverse transformation is 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 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 This simplified interaction term is equivalently expressed as The dispersion equation on the approximation that 1 The above expression is a dispersion relation with some additional term 2 B γ = DS+gμ H + μE  that tends to show how the dispersion term depends on the field components.Moreover, the constant .

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 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, 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 1 .
Figure 1.(color online) The figure indicated shows the variation of dispersion function for different gamma values.
in magnetization versus temperature at some specific γ is shown in Figure2below for low temperature limits.

Figure 2 .
Figure 2. Change of magnetization versus temperature at different gamma values.

Figure 3 .
Figure 3. Dependence of magnetization on the parameter γ, which can be tuned by explicitly on the applied fields and the nature of magnetic anisotropy.