# Magnetic Nanoparticles with Core Shell: Macroscopic Model and Coercive Field

Authors: O. Hudak; M. Hudak
Journal: Advances in Materials Science and Engineering
 (2010)

---

## Abstract

There exists nonhomogeneous structure in mechanosynthesized nanoparticles
consisting of an ordered core surrounded by the shell. The shell may be structurally and magnetically disordered, or it may be ordered. These nanoparticles
are found to be roughly spherical. We formulate the macroscopic model for the description of magnetic properties of nanoparticles with core-shell structure.
The case of spheroids oriented in the same direction of polar axes is considered. There exits two coercive fields:Hcγ and Hcψ. They are depending on quantities o
and x, where x=V2/(V1−2V2), and o=IF/IS , here V1 is the volume of the nanoparticle, and V2 is the volume of the core. IF is the magnetization of the core, and IS is the
magnetization of the shell. There exists a transition from one type of the coercive field to another, changing x and o. We have found qualitative explanation
of observed enhancement of coercive field and of its d-dependence in NiFe2O4 nanoparticles.


---

## Body

## 1. Introduction

Ferromagnetic and (ferri-) or antiferromagnetic nanoparticles, mechanochemically prepared, have a shell structuredifferent from the core structure. The shell has mostly disordered structural and magnetic properties, while the core is magnetically ordered. In this paper, we formulate a model describing interaction of ellipsoid with ferromagnetic/antiferro- (ferri-)magnetic/ordering with shell (ellipsoid part which is not in the volume of the core ellipsoid). We study here the description of coercive fields of magnetic nanoparticles with spheroid shape with polar axes in the same direction for both spheroids. Due to fact that nanoparticles are small (smaller than100A°), we consider a single domain magnetic system. While porosity of ceramics prepared from such nanoparticles is in the range up to nearly 33 percents, and such pores do not disappear during further sintering, they undergo coalescence together with grain growth. Only subsurface pores have the possibility to go out in the surface, see in [1]. For the case of a spheroid with a pore, it can be easily shown that the pore may enhance the coercive field, or the coercive field is smaller than the bulk coercive field.MechanosynthesizedNiFe2O4 nanocrystalline nanoparticles were studied in [2]. There exists nonhomogeneous structure in these nanoparticles consisting of an ordered core surrounded by a disordered grain boundary-surface region. The core has a fully inverse spinel structure with a collinear spin alignment. The surface layer (shell) is structurally and magnetically disordered due to nearly random distribution of cations and the canted spin structure. It exhibits a reduced nonsaturating magnetization, an enhanced coercivity, and a shifted hysteresis loop. The thickness of the shell of about 1 nm is estimated from size-dependent magnetization measurements. The shell thickness of crystallite radius ratio is nearly 23 percent. The nanoparticles are found to be roughly spherical with the core-shell structure. In [2], it was shown that the core part of the nanoparticle is structurally ordered and of nearly sphere shape, and the shell part is disordered structurally. The spins in the shell region are found to be canted (28° and 40° in magnetic sublattices), and the λ parameter (parameter which in spinels describes how cations are distributed in spinels in their sublattices) equals nearly to 0.67, in the core the λ parameter is equal to 1. It is observed that an enhanced magnetic hardness, that is, the coercive field of the material has the value HC≈0.35 T which is about 35 times larger than that for the bulk material HC≈0.01 T. The field-cooled hysteresis loop is not symmetrical and is shifted to the left, and the shift is δHC≈0.03 T. Such a shift and an asymmetric hysteresis loop were explained [3, 4] in the exchange model to be present due to the exchange coupling between the collinear spins in the core and canted spins in the shell. The saturation magnetization is dependent on the diameter. The growth of NiFe2O4 nanoparticle is accompanied by a decrease of coercivity field value.The nonuniform core-shell structure of nanoparticles has been reported inMgFe2O4 [5] and in ball-milled nanocrystalline LiNbO3 [6]. Mechanochemically prepared nanoparticles MgFe2O4 [5] have an unusual property: enhanced magnetization. They have core-shell structure. Inner core exhibits a partly inverse structure with a Neel-type collinear spin alignment. The shell has a nonequilibrium cation distribution and a canted spin arrangement. Although spin canting effects tend to reduce the magnetic moment, the magnetization enhancement is attributed by authors to the nearly random distribution of magnetic cations in the shell region. On heating above 623°K, the mechanochemical MgFe2O4 relaxes to the magnetic state which is similar to the bulk one [5]. The thickness t of the shell assuming the spherical shape is estimated to be t≈1 nm and is comparable to the lattice constant.NanocompositeCuFe2O4/CeO2 [7] has a ferrimagnetic behaviour. The nanoparticle size of CuFe2O4/CeO2 decreases with increasing concentration of CeO2. Added CeO2 is dissolved in the CuFe2O4, or it is coated on the surface of CuFe2O4 nanoparticles [8]. There is then a surface layer-shell of CeO2. This form of core (shell) structure is of different type, as in NiFe2O4. Nevertheless, also in this case there is a core and a shell, and they have magnetic properties [9]. It is possible to model such materials as materials in which there is a core-shell structure similar to that in nanoparticles prepared by mechanochemical preparation. Modelling of CeO2-x in the LDA+U formalism has been used to account for the strong on-site Coulomb repulsion among the localized Ce 4f electrons. Results show that by choosing an appropriate U, it is possible to consistently describe structural, thermodynamic, and electronic properties oCeO2-xf, which enables modelling of redox processes involving ceria-based materials. We may expect that similar nanocomposites may be prepared using CeAl2 as a shell. We have modelled CeAl2 nanoparticles [10] by 1d Hubbard model. Surface valence phase transition between magnetic and nonmagnetic states was studied by us [10] in these nanoparticles. Coated CeAl2 may form the core-shell structure. The surface properties of the shell may be described by our model [10] (surface quantum phase transition) and also by the model of magnetic nanoparticles with core-shell in which the coercive field is calculated here.There is different surface anisotropy for the cubic and spherical nanoparticles [11]. Low temperature properties: saturation magnetization, coercivity, and spin canting are similar, while the blocking temperature TB (temperature at which a maximum of susceptibility in nanoparticle is achieved) is different for γ-Fe2O3 for the cubic and spherical shape. Theoretical calculations and experimental magnetic studies for nanoparticles of different shapes are [12–17] scare, see in [11]. For cobalt ferrite CoFe2O4, it was found that blocking temperature TB and the saturation magnetization are independent of the shape, and the value of the coercivity field HC is larger for spherical shape than for the cubic one. In cobalt ferrite there is a larger magnetic anisotropy. For γ-Fe2O3, intrinsic magnetocrystalline anisotropy is small. Monte Carlo simulations are based on Heisenberg model for exchange interaction between spins in the core and in the shell experiencing different anisotropies [11]. The surface thickness is taken to be equal to one lattice constant. Experimentally it was found that the spherical nanoparticles are more ordered than cubic nanoparticles. Low temperature measurements show however that their microscopical magnetic structure properties are similar, indicating similar crystalline properties. The behaviour is ascribed to the effects of morphology on the effective surface anisotropy and the small value of intrinsic anisotropy of γ-Fe2O3. The magnetic state in the core is ferrimagnetic and in the surface layer is of the spin-glass type. We will take into account the morphology effects in our paper by considering the ellipsoidal shape of nanoparticles. This shape enables us to describe the spherical shape of nanoparticles, the disc-like shape, or rod-like shape for extremal values of lengths axes in the ellipsoid. The case of cubic nanoparticles can be easily considered in the same way.In the next section, we formulate the core-shell model for description of magnetic nanoparticles. In the third section, the Model Free Energy is found. In the next sections we study the free energy for spheroids and spheres oriented in the same direction. Experimental behaviour of coercive field dependence on diameter of nanoparticle is compared with theoretical predictions. We have found that there exist two coercive fieldsHcγ and Hcψ. They are depending on quantities o and x, defined by x=V2/(V1-2V2), and by o=IF/IS, here V1 is the volume of the larger ellipsoid, and V2 is the volume of the smaller ellipsoid. IF is the magnetization of the volume where smaller ellipsoid occurs, and IS is the magnetization of the shell. There exists a transition from one type of the coercive field to another changing x and o within range of their possible values. Critical values xc and oc for which values of both coercive fields are equal can be found. For spheroidal shapes of ellipsoids for x>xc, the field Hcγ is the coercive field which may be observed in experiments. For volumes V1 and V2 for which x<xc, the field Hcψ is the coercive field and should be observed. The ratio of the coercive field value for nanoparticle to the bulk coercive field value is different for x>xc and for x<xc. Enhancement of the coercive field may occur in nanoparticles changing the diameter d.


## 2. Core-Shell Model

As it was shown above, there are many nanoparticle systems in which there exists core-shell structure. They are found in nanoparticles prepared, for example, by mechanochemical method and also by other methods. For an overview of the mechanochemical methods see in [18].We will describe such nanoparticles by the following Core-Shell model. Due to small diameter of nanoparticles we will assume that they are magnetically monodomain nanoparticles. This is true for nanoparticles with diameter less than100A° [19]. The shape of nanoparticles described above is of the sphere type, mostly, there are deviations from this shape in the core and in the shell. To describe these observations, we assume that the nanoparticle has an ellipsoid shape in general. The ellipsoid shape is described by three parameters a, b, and c (lengths of axes). To describe existence of the core and of the shell, we consider another ellipsoidal “nanoparticle” which is located inside the first one. The centre of that another “nanoparticle” may have almost any position inside the first ellipsoid, assuming that the second ellipsoid is as a whole inside the first ones volume.We will assume that main axes in both ellipsoids coincide with the polar (light magnetic) crystallographic axis of the nanoparticle. In general, these two axes (light and main axes in both ellipsoids) are oriented in different directions. We will denote the angle between polar axis in both byθ. The case of zero angle, θ=0, corresponds to parallel polar axes. The case of angle θ=π/2 corresponds to perpendicular polar axes. Here in this paper we consider the case of θ=0 angle which corresponds to the shell and the oriented core in the same direction.To describe core-shell magnetic structure we will assume that the larger ellipsoid has a homogeneous magnetizationIS, and that the smaller ellipsoid has a homogeneous magnetization IF-IS. Thus, there is a total magnetization IF-IS+IS=IF in the volume where smaller ellipsoid is localized in the larger one, see Figure 1. The magnetic structure of the material in core and in shell may be ferromagnetic or ferri- (antiferro-) magnetic. We will assume that sublattice magnetizations are coupled strongly, and that then we may consider magnetizations IF and IS (resulting magnetization) only. We will neglect also surface magnetic phenomena. The magnetic field around an ellipsoid of the sphere shape with homogeneous magnetization IF-IS is given in our case by the dependence [19](1)(IF-IS)4π3(d/2)3r3(2ercos(θ)+eρsin(θ)),
where r is the distance from the centre of the sphere, and d is a diameter of the sphere. Here, er is a unit vector in the direction of radius directing from the centre, and eρ is a unit vector in the direction perpendicular to the vector er. For larger distances from the smaller ellipsoid surface, we may neglect this contribution. This approximation corresponds to smaller inner ellipsoid volume with respect to the larger ones volume.Figure 1
Core-shell magnetic structure. The larger ellipsoid has a homogeneous magnetizationIS, and the smaller ellipsoid has a homogeneous magnetization IF-IS. There is a total magnetization IF-IS+IS=IF in the volume where smaller ellipsoid (core) is present in the larger one. The angle of the magnetization of the shell volume with the direction of the main axis is ψ, and the angle of the magnetization of the core with the direction of the main axis is γ.Let us consider first the case in which both ellipsoids are spheres. Then we will consider the case of spheroids (ellipsoids with the same two main axes lengths). The angle of the magnetization of the first (larger) ellipsoid with the direction of the main axis isψ, and the angle of the magnetization of the second (smaller) ellipsoid with the direction of the main axis is γ, see Figure 1.


## 3. Model-Free Energy

Interaction of the core magnetization with the shell magnetization will be described by the effective-field method. We will assume that an effective magnetic field of the larger ellipsoid is acting on the smaller ellipsoid magnetization. On the other hand, the effective magnetic field of the smaller ellipsoid outside and inside this ellipsoid which is acting on the larger ellipsoid magnetization will be neglected.The magnetic anisotropy of the shell structure may be described by the unit vectoreS, and the contribution to the free energy is of the form KS(eS·IS)2, where KS is the anisotropy constant for the shell. The anisotropy of the core may be described by the unit vector eC, and the contribution to the free energy is of the form KC(eC·IF)2, where KC is the magnetic anisotropy constant for the core. When the structure is random then the anisotropy is random, and we can use mean anisotropy of the magnetic structure in the shell. Random anisotropy in the shell (for the constant λ unequal to one) may be described by random orientation of the vector e and by the different values of the constant K. We will neglect magnetic anisotropy in this paper as the magnetic anisotropy constants KS and KC are small in NiFe2O4.We will assume that the nanoparticle is monodomain, and the free energy of the magnetic system will contain contributions similar to those for the ellipsoid with an ellipsoidal hollow [19]. The inner ellipsoid is corresponding to the hollow in [19], however we consider this ellipsoid as a magnetic system with the magnetization IF-IS then the total magnetization in the volume V2 of smaller ellipsoid is IF-IS+IS=IF. The volume of the larger ellipsoid is V1.Let us denote susceptibilities for larger and smaller ellipsoids asχS0 and χF0 respectively, they are defined as(2)χS0=1Nata2+Nbtb2+Nctc2,
where Ni is a demagnetization factor for the direction i(i=a,b,c) for the larger ellipsoid, and where ti is a directional cosine of the magnetization vector IS for the direction i, and(3)χF0=1Nαtα2+Nβtβ2+Nγtγ2,
where Nj is a demagnetization factor for the direction j(j=α,β,ϵ) for the smaller ellipsoid, and where tj is a directional cosine of the magnetization vector IS for the direction j.The free energyFS of the larger ellipsoid is(4)FS=12V11χS0IS2.The free energyFF of the smaller ellipsoid is(5)FF=12V21χF0(IF-IS)2.The free energyFSF describing the interaction of the effective magnetic field –(1/χS0)IS of the larger ellipsoid magnetization with the magnetization IF-IS of the smaller ellipsoid within the volume V2 of the smaller ellipsoid is(6)FSF=12V21χF0IS·(IF-IS).The free energy describing the interaction of the magnetic fieldH with the magnetization IS of the larger ellipsoid within the volume V1 is(7)FH1=-V1IS·H.and the free energy describing the interaction of the magnetic fieldH with the magnetization IF-IS within the volume of “second nanoparticle” V2 is(8)FH2=-V2(IF-IS)·H.Let us introduce the spherical coordinates for magnetizationsIS and IF:(9)IS=IS(l,m,n),
where IS is the amplitude of the magnetization IS of the larger ellipsoid, and (l,m,n) are coordinates of the directional unit vector of magnetization IS, and(10)IF=IF(L,M,N),
where IF is the amplitude of the magnetization of the larger ellipsoid, and (L,M,N) are coordinates of the unit vector of the direction of magnetization IF. Note that l2+m2+n2=1 and that L2+M2+N2=1.


## 4. Free Energy: Spheroids Oriented in the Same Direction

Let us consider the case of spheroids with polar axes in the same direction. This case corresponds to the cases of mechanochemically prepared nanoparticles in which the shell has the direction of polar axis the same as the direction of polar axis of the core spheroid andθ=0. The free energy for this case has the form F=FS+FF+FSF+FH1+FH2.The total free energyF for the case of spheroids with the polar axis in the same direction has the form(11)F=F0+12IS2(Na-Nb)(V1-2V2)cos2(ψ)+12IF2(Na-Nb)V2cos2(γ)-V2ISIFcos(γ+ψ)(Na-Nb)(cos2(γ)-cos2(ψ))-(V1-V2)ISHcos(ψ)-V2IFHcos(γ),
where F0 is the free energy which does not depend on angles γ and ψ. We assume that the deviation of magnetization vectors IS and IF from polar axes in core and in shell is small, thus angles ψ and γ are small angles with values near 0.The free energy constantF0 has the form(12)F0=12V1NbIS2+12V2NbIS2(IF2+IS2)-V2NbIS2.


## 5. The Case of Spheres

Let us now discuss simple case of two spheres. Note that centres of these spheres are not necessary in the same point. In this case we haveNa-Nb=0 and Nα-Nβ=0. Thus in the case of the spheres there are only last two terms present in (11). The total free energy F for the case of spheroids with the polar axis in the same direction is(13)F=F0-(V1-V2)ISHcos(ψ)-V2IFHcos(γ),
where F0 is the free energy which does not depend on the angles γ and ψ (see above).Note that forH>0, the lowest free energy (13) corresponds to angles ψ=0 and γ=0. For H<0, the lowest free energy (13) corresponds to angles ψ=π and γ=π. Thus spins in both spheres are oriented in the direction of the magnetic field. The magnetization M, its component in the direction of main axis, has in this simple case the form:(14)M=(V1-V2)IScos(ψ)+V2IFcos(γ).Thus the difference of magnetization forH>0 and for H<0 has the form for the angles ψ and γ which minimizes free energy:(15)ΔM=2(V1-V2)IS+2V2IF.The coercive field has the zero valueHc,±=0. The magnetic susceptibility χ, has in this case value(16)χ=δMδH=2ΔMδ(H).In fact the delta function above will be a peak like (a gaussian peak) with a finite and with an amplitude given by2(V1-V2)IS+2V2IF.


## 6. The Case of Spheroids

Let us now discuss the case of spheroids with the same shape, in whichNa-Nb is different from 0 and Nα-Nβ is different from 0, Na-Nb=Nα-Nβ. Assuming small angles ψ and γ we can find the following form of the free energy expansion in ψ and γ:(17)F=F1+IS2(Na-Nb)(V1-2V2)×[ψ2(-12-x·o+H1+x2IS(Na-Nb))+γ2(-x1+o22+Hxo12IS(Na-Nb))],
where F1 is a constant not dependent on angles γ and ψ. As we can see from (17) using the method from [19], there are two coercive fields, which we denote as Hcγ and Hcψ, at which coefficients of ψ2 and γ2 become zero:(18)Hcγ=IS(Na-Nb)(1+o2)o(forx≠0 ) and(19)Hcψ=IS(Na-Nb)1+2xo1+x.Here the quantitieso and x are defined by x=V2/(V1-2V2) and o=IF/IS. The coercive field which will be measured in experiments is that field from coercive fields Hcγ and Hcψ, which is smaller and in this case both magnetizations IS and IF are reoriented. When this field is approached decreasing the magnetic field from above, there will be a change in the sign of the coefficient corresponding free energy terms. Note that for small values of x for which the core is small, this field becomes Hcψ=IS(Na-Nb)(1+x(2o-1)) and thus is increasing with increasing x for 2o>1, is larger than the bulk coercive field. For 2o<1 it is smaller than the bulk coercive fieldFrom (18) and (19) we can find conditions for critical values of x and o quantities for which values of both coercive fields Hcγ and Hcψ are equal:(20)1+2xo1+x=(1+o2)o.The coercive fieldsHcψ and Hcγ change their role as the observed coercive field at those values of o and x for which (20) holds.The critical valuexc at which this change takes place is given by:(21)xc=1-o+o2o2-1.We see that ato=0 where there is zero magnetization in the core, we have the critical value of x such that V1=V2. At o=∞ we have the critical value of x such that V2=V1/3. There is zero magnetization IS=0 in the shell now. For o=1+, we find xc=+∞ and V2=V1/2, for o=1- we obtain xc=-∞ and again V2=V1/2. There is the same magnetic moment in both core and shell for o=1+ and o=1-.For the same spheroidal shapes of ellipsoids forx>xc the field Hcγ is the coercive field observed in experiments. For volumes V1 and V2 for which x<xc, the field Hcψ is the coercive field observed, now the smaller spheroid is very small. The bulk value of the coercive field as can be easily shown is given by(22)Hcbψ=IS(Na-Nb).


## 7. Discussion

The coercive field depends on the diameterd in nanoparticles prepared by mechanochemical methods. The nanoparticle growth of NiFe2O4 is accompanied by the decrease of coercivity [2]. The surface to volume fraction in this material decreases with increasing annealing temperature. In [20], it was found that increasing the diameter d of the nanoparticle the width of the shell decreases for BaFe12O19 and for Co0.56Zn1.44W. While this fact was found to be true for nanoparticles with the diameter d from the interval 100A° to 400A° we will assume that the same qualitative dependence of thickness t on the diameter d is present for diameter d<100A° in NiFe2O4 nanoparticles. From two coercive fields, which we have found above, the field Hcψ depends on the volume of nanoparticles: Hcψ=Hcbψ((1+2xo)/(1+x)). Here Hcbψ is a bulk value of the coercive field, Hcbψ=IS(Na-Nb). This coercive field changes with the diameter d as(23)dddHcψ=Hcb2o-1(1+x)2dxdd.As we can see the change of the coercive fieldHcψ on the diameter d is given by dx/dd. This change is found to be(24)dxdd=V1V2′-V2V1′(V1-2V2)2.HereV1′=dV1/dd and V2′=dV2/dd. We can find that the coercive field Hcψ is decreasing for o>1/2 with increasing diameter d when the condition t(d)>d(dt(d)/dd), where t(d) is d-dependent width of the shell. d-dependent width of the shell is decreasing with increasing diameter d as is observed in some experiments [20], in [2], it is assumed that the thickness of the shell in NiFe2O4 does not depend on d. In both cases the condition t(d)>d(dt(d)/dd) is fulfilled. While the spin IS of the shell is canted from the spin IF of the core by angles 20° and 40° depending on the sublattice, the value of the spin in the shell may be expected to be smaller than the value of the spin in the core due to disorder in the shell and the condition o>1/2 may be assumed to be satisfied. Thus our theory qualitatively explains observed decrease of the coercive field with increasing diameter d of nanoparticles for nanoparticles with the width of the shell independent of diameter or with decreasing width of the shell with increasing diameter d. In the shell of nanocrystalline complex oxides a nonequilibrium cation distribution was found accompanied by a deformation of polyhedron geometries [21]. Note for example that in some nanoparticles like BPFTO ((Bi0.9Pb0.1)(Fe0.9Ti0.1)O3), the width of the shell with increasing diameter may increase. For small increase of the shell width which authors are describing as almost constant, see [22].The coercive field observed in experiment inNiFe2O4 is enhanced by the factor 35 with respect to the coercive field of the bulk [2]. The ratio of the coercive field Hcψ to the bulk coercive field Hcbψ is given as (1+2xo)/(1+x). We then obtain that the value of o corresponding to the enhancement factor in NiFe2O4 is o=35+17/x and we see that for small volumes V2 of the core the ratio o of spin values is too large. The spins in the shell do not have large value. For larger values of the volume V2 and larger values of shell spins the enhancement factor may be expected to have the value 35. In our approximation of a small x the enhancement of the coercive field may exist, however this enhancement does not reach the value 35; one can find (d/dx)((1+2xo)/(1+x))=2(2o-1)/(1+x)2>0 for o>1/2 and for x≈0. We have (1+2xo)/(1+x)≈1+x(2o-1)>1, however this ratio does not reach the value 35 in NiFe2O4.



---

*Source: 909810-2010-07-04.xml*