Ferromagnetic-like behavior of Bi0.9La0.1FeO3-KBr nanocomposites.

We studied magnetostatic response of the Bi0.9La0.1FeO3- KBr composites (BLFO-KBr) consisting of nanosized (≈100 nm) ferrite Bi0.9La0.1FeO3 (BLFO) conjugated with fine grinded ionic conducting KBr. When the fraction of KBr is rather small (less than 15 wt%) the magnetic response of the composite is very weak and similar to that observed for the BLFO (pure KBr matrix without Bi1-xLaxFeO3 has no magnetic response as anticipated). However, when the fraction of KBr increases above 15%, the magnetic response of the composite changes substantially and the field dependence of magnetization reveals ferromagnetic-like hysteresis loop with a remanent magnetization about 0.14 emu/g and coercive field about 1.8 Tesla (at room temperature). Nothing similar to the ferromagnetic-like hysteresis loop can be observed in Bi1-zLazFeO3 ceramics with z ≤ 0.15, which magnetization quasi-linearly increases with magnetic field. Different physical mechanisms were considered to explain the unusual experimental results for BLFO-KBr nanocomposites, but only those among them, which are highly sensitive to the interaction of antiferromagnetic Bi0.9La0.1FeO3 with ionic conductor KBr, can be relevant.

Magnetic, structural and electrophysical properties of pristine BiFeO3

Multiferroics are ideal systems for fundamental studies of couplings among the order parameters of different nature, e.g. ferroelectric (FE) polarization, structural rotational antiferrodistortion (AFD), ferromagnetic (FM) and antiferromagnetic (AFM) long-range order parameters[1-6]. Magnetoelectric (ME) coupling is especially important for the most of multiferroic fundamental studies and applications[1-5].

Bismuth ferrite BiFeO3 (BFO) is the unique multiferroic with the large FE polarization (more than 60 μC/cm2 at RT) and AFM order coexisting up to room (RT) and elevated temperatures[7,8]. BFO exhibits unusual electrophysical properties, such as conduction and magnetotransport enhancement at domain walls[9-15]. Specifically, bulk BFO exhibits AFD long-range order at temperatures below 1200 K; it is FE below Curie temperature T ≈ 1145 K and is AFM below Neel temperature T ≈ 645 K[16]. It is known that pristine bulk BFO is characterized by a cycloidal modulation of the magnetization superimposed on the AFM G-type magnetic structure[3]. Although the linear magnetoelectric effect is symmetry forbidden on a macroscopic scale, the Dzyaloshinskii-Moriya mechanism can be used locally in BFO to achieve an electric switching of the spin cycloid[17-20].

Magnetic, structural and electrophysical properties of Bi1-zLazFeO3

Chemical doping of pristine BFO with lanthanum La ions having ionic radius (1.16 Å[21]) similar to that of bismuth Bi ions (1.17 Å[20]) causes a structural transition from the polar rhombohedral phase to the anti-polar orthorhombic phase which is accompanied with a slight decrease in the unit cell volume[22-24]. The structural transition is associated with significant changes in the unit cell parameters, namely, the modification of the Fe–O chemical bond lengths (a shortening of the bond length between the iron ions and apical oxygen ions and an elongation of the chemical bond between the iron ions and the oxygen ions located in the basal ab plane), an increase in the Fe–O–Fe angles, decrease in the oxygen octahedra rotations and tilts etc. The mentioned changes of the structural parameters lead to a reduction of the polar displacement of the ions and to significant change of the magnetic properties of the compounds near the rhombohedral-orthorhombic phase boundary. The diffraction data obtained at RT for the compounds Bi1-zLazFeO3 (BLFO) indicate a single phase rhombohedral state for the compounds with z < 0.16, two phase structural state assuming a coexistence of the polar rhombohedral and the anti-polar orthorhombic phases in the concentration range 0.16 < z < 0.19; further increase in the dopant content leads to a stabilization of single phase orthorhombic state (incommensurately modulated anti-polar orthorhombic phase, anti-polar or non-polar orthorhombic phase) depending on the La content as described elsewhere[22,25,26].

Sol-gel synthesis method used to prepare BLFO solid solution leads to a formation of oxygen stoichiometric compounds those magnetic properties are governed by the exchange interactions between Fe ions being in 3+ oxidation state. Increase in the concentration of the lanthanum La ions leads to a gradual reduction in the critical magnetic field associated with a disruption of the spatially modulated spin structure and the compound with lanthanum content z ~ 0.17 is characterized by weak ferromagnetic structure at RT assuming complete disruption of the spatially modulated structure. The structural transition from the rhombohedral to the orthorhombic phase is accompanied by a significant increase of the remanent magnetization and coercive field, which is observed for the compounds having dominant or single phase orthorhombic structure.

It should be noted that transport properties of the BLFO compounds are also strongly dependent on the La concentration, and the conductivity of these solid solutions gradually increases with the dopant content increase up to z = 0.15[27]. Impedance measurements have shown mainly p-type conductivity of the compounds having single phase rhombohedral structure which significantly increases with temperature[22].

This work studies magnetostatic response of the newly synthesized Bi0.9La0.1FeO3 - KBr composites (BLFO-KBr) consisting of nanosized (≈100 nm) ferrite Bi0.9La0.1FeO3 (BLFO) conjugated with fine grinded ionic conducting KBr. We revealed that when the fraction of KBr increases above 15%, it demonstrates ferromagnetic-like hysteresis loop with a remanent magnetization about 0.14 emu/g and coercive field about 1.8 Tesla (at room temperature). To the best of our knowledge nothing similar was reported previously in the literature. Different physical mechanisms were considered to explain the unusual experimental results for BLFO-KBr nanocomposites. The original part of the manuscript is structured as follows. Magnetostatic and electrophysical properties of BLFO–KBr nano-composites are discussed in section II. Theoretical explanation of the ferromagetism in BLFO–KBr nano-composites is given in section III. Section IV is a brief summary. Auxiliary experimental results and samples characterization are given in the Supplement.

Magnetic and Electrophysical Properties of BLFO–KBr Nano-Composite

Nanocomposite preparation and characterization

The samples of nanograined BLFO ceramics have been prepared by the aqueous sol-gel method using hydrates Bi(NO3)3·5H2O, Fe(NO3)3·9H2O and lanthanum nitrate as starting materials. The prepared sol-gel samples were heated at 800 °C for 1.5 h[28,29]. The composites have been prepared by thorough mechanical mixing of the single phase ferrite (Bi0.9La0.1FeO3) and grinded KBr taken in the next mass ratios – 90:10, 85:15, 80:20, 70:30, 50:50 respectively. The obtained powders have been uniaxially pressed at compacting pressures of 5 GPa to obtain tablets with 5 mm diameter and 1 mm thickness.

The compound Bi0.9La0.1FeO3, which has been mainly used to prepare the composites, is characterized by the grains with average radius of about 100 nm, and the small amount of the grains have the average radius of about 35–40 nm (the image obtained by scanning electron microscopy (SEM) and the histogram of the grain sizes distribution are shown in the Fig. 1(a). SEM measurements performed for the composites have allowed to clarify their morphology and to specify a distribution of the constituent phases [Fig. 1(b)]. The crystallines having micrometer size are attributed to the KBr phase as determined by EDS analysis carried out for the selected areas of the SEM images [see Supplement, Fig. S1]. Elemental composition of the constituent components of the composites is confirmed by EDS data presented in the Supplement.

Figure 1

(a) SEM image of the compound Bi0.9La0.1FeO3 obtained at room temperature. The inset shows a distribution of crystalline size over the SEM image of the compound. (b) SEM image of the composite 50% Bi0.9La0.1FeO3 −50% KBr, the spots marked by different colors denote the areas attributed to the ferrite phase, KBr phase or their mixture as determined by EDS measurements (detailed results of EDS experiments are presented in the Supplement).

Characteristic size of the grains attributed to KBr phase is in the range of 0.100–10 μm. The crystallines ascribed to the ferrite phase are characterized by rectangular-shaped form, ferrite grains are mainly agglomerated into clusters having size of about 1–2 microns, while separate grains are characterized by nanometer size as confirmed by the SEM images obtained for the composites and the single phase ferrite [compare Fig. 1(a) and 1(b)]. The SEM data have confirmed homogeneous distribution of the constituent phases in the composites which facilitates active chemical processes in the boundary region of both phases.

The modified sol-gel synthesis method and preparation conditions[28,29] used to synthesize nanoscale powder of Bi1-zLazFeO3 (z = 0, 0.05, 0.1) lead to a formation of single phase stoichiometric ferrites utilized to form the composites. The X-ray diffraction measurements performed for the ferrite compounds have confirmed their single phase rhombohedral structure, the XRD patterns recorded at room temperature have been successfully refined assuming non-centrosymmetric rhombohedral space group R3c.

Under standard conditions, potassium bromide (KBr) is a white nonmagnetic crystalline powder which exhibits high ionic conductivity. Phase purity of the commercial KBr powder has been confirmed by the XRD measurements performed by the authors, the obtained results are in accordance with the structural data provided by the producer (see Supplement, Fig. S2). The XRD patterns obtained for the composites have confirmed a coexistence of two constituents having the phase ratios in accordance with the mentioned chemical formulas. The diffraction peaks attributed to the constituents do not show any modification of the reflections positions, their width, asymmetry etc. as compared to the reflections observed on the diffraction patterns obtained separately for BLFO and KBr (see Supplement, Fig. S2). Phase purity of the constituents assumes that the magnetic properties of the composites are determined only by the ferrite component.

Under standard conditions, potassium bromide (KBr) is a white nonmagnetic crystalline powder which exhibits high ionic conductivity. Phase purity of the commercial KBr powder has been confirmed by the XRD measurements performed by the authors, the obtained results are in accordance with the structural data provided by the producer (see Supplement, Fig. S2). The XRD patterns obtained for the composites have confirmed a coexistence of two constituents having the phase ratios in accordance with the mentioned chemical formulas. The diffraction peaks attributed to the constituents do not show any modification of the reflections positions, their width, asymmetry etc. as compared to the reflections observed on the diffraction patterns obtained separately for BLFO and KBr (see Supplement, Fig. S2). Phase purity of the constituents assumes that the magnetic properties of the composites are determined only by the ferrite component.

For comparison with the composite, the temperature dependences of magnetization for the nanograined ferrites Bi1-zLazFeO3 (z = 0, 0.05 and 0.1) have been measured in the temperature range 300–1000 K using MPMS SQUID VSM magnetometer (Supplement, Fig. S3). The isothermal magnetization measurements have been done at room temperature and at T = 5 K (Supplement, Fig. S4). The obtained results demonstrate a quasi-linear behavior of the magnetization dependencies (see Supplement, Fig. S4). An increase of the magnetization observed for the compounds with the dopant content up to ~13% in strong magnetic fields (5–8 Tesla) is associated with a partial disruption of the spatially modulated magnetic structure, the transition is reversible and after exposure to magnetic field the magnetization of the compounds returns to nearly nullified values[30,31]. The very small values of the remanent magnetization (0.015–0.02 emu/g) observed for the lightly-doped compounds mainly depends on the dopant content and is caused by uncompensated spin magnetic moments formed in the surface layer of the crystallites due to disruption of the modulated magnetic structure. The compounds Bi1-zLazFeO3 with the dopant content z ≥ 0.15 are characterized by the anti-polar orthorhombic structure and have remanent magnetization of ~0.15–0.2 emu/g[31]. The compound Bi0.9La0.1FeO3 having rhombohedral crystal structure and modulated magnetic structure with almost zero remanent magnetization has been selected as the most appropriate compound for formation of the composites to estimate an influence of the constituents on the physical properties of nanocomposites.

Magnetic and electrophysical properties of (Bi0.9La0.1FeO3)x–(KBr)1-x nano-composites

The isothermal dependencies of magnetization (M) on magnetic field (H) of the nanocomposites (Bi0.9La0.1FeO3)x - (KBr)1-x (shortly (BLFO)-(KBr)) are shown in Fig. 2 for the fractions x = 100% (red symbols), 90% (blue), 85% (green), 80% (violet), 70% (brown) and 50% (black). The M-H-dependencies have been obtained in the fields up to 14 T using a magnetometer from Cryogenic Ltd. We have also confirmed experimentally that KBr matrix without ferrite compound has no magnetic response. Figure 2, left inset shows KBr-fraction (1-x)-dependencies of the remanent magnetization (M), maximal magnetization (M) at the maximal field Hm = 14 T and the coercive field (H) for (BLFO)x-(KBr)1-x composites.

Figure 2

Dependencies of magnetic moment versus applied magnetic field for the composites (BLFO)x – (KBr)1-x at RT (the data for composites with the fractions of x = 100% are denoted by red symbols, 90% - blue, 85% - green, 80% - violet, 70% - brown and 50% - by black symbols. Left inset shows KBr fraction (1-x)-dependencies of Mr, Mm and Hc for the composites; right inset shows enlarged part of the M(H) dependences.

The M(H) dependencies of the composites with the fraction of KBr less than 15% are similar to those observed for the compounds without KBr content (see red, blue and green symbols in Fig. 2). Further increase in the KBr fraction significantly changes the magnetic properties of the composites (BLFO)x–(KBr)1-x, namely leads to the appearance of ferromagnetic-like hysteresis of M(H) loops (see violet, brown and black symbols in Fig. 2). Note that the increase of the KBr fraction (1-x) above 15% leads to the step-like decrease of Mm, and increase of Mr and Hc values of the composites (see Fig. 2 insets). The composite containing 20% of KBr is characterized by Mr ≈ 0.15 emu/g and Hc ≈ 1.8 T at RT (violet symbols in Fig. 2 and inset). Further increasing in the KBr fraction leads to the gradual decrease of Mr, which is mainly associated with a dilution of the magnetically active ferrite component, wherein Hc attributed to the composites remains nearly unchanged (see violet, brown and black symbols in Fig. 2 and inset). The composite with equal fractions of the constituents (x = 50%) is characterized by the opened magnetization loop specific for ferromagnetic-like materials. Magnetization measurements performed at temperature ~5 K do not reveal any notable changes of the composites’ magnetic properties, slight increase in magnetization observed at temperature of ~5 K can be justified by a reinforcement of the magnetic interactions occurred at low temperatures.

Different physical mechanisms can be responsible for the appearance of ferromagnetic-like response of (BLFO)x–(KBr)1-x nanocomposites, where BLFO with La content less than 16% is purely antiferromagnetic in the bulk at temperatures lower than 650 K and KBr is magnetically inactive material. Anyway, the changes observed in the magnetization behavior of the composites can be completely attributed to the changes occurred in the ferrite component. One of the most plausible models, describing the changes of the magnetic properties observed for the composites (BLFO)x–(KBr)1-x depending on the KBr content assumes intensive chemical processes in the vicinity of ferrite particle - salt interface occurred during the composites formation. As a result of the redox reactions the Fe ions near the complex oxide – alkali halide interface can partially change their effective oxidation state from 3+ down to 2+. An uncompensated magnetic moment formed due to a difference between the magnetic moments attributed to the Fe3+ and Fe2+ ions within the antiferromagnetic matrix should lead to a notable increase of remanent magnetization of the composites and to a decrease of coercivity. Remanent magnetization estimated for the composites does not significantly change with an increase of KBr and it remains nearly stable for the composites with KBr fraction of about (20–50)%, which can be explained assuming that the chemical processes occur only in the thin sub-surface layer of the crystallites of the ferrite component. An average size of the ferrite crystallites is about 100 nm as confirmed by scanning electron microscopy measurements (see Fig. 1), that causes the presence of a concentration threshold (~20%) affecting the chemical processes and the composites having smaller KBr fraction does not show significant difference in magnetic and transport properties as compared to the ferrite compounds (see Supplement, Fig. S4). The manifestation of this effect might be a threshold of about ~20% of KBr content. This idea is verified by the isothermal magnetization curves obtained for the composites with KBr phase ratios 15% and 10% which have revealed almost absence of the remanent magnetization, while distinct opened hysteresis loop is observed for the composite having 20% of KBr content, so pointing at certain threshold. Notably, that the single-crystalline BiFeO3 powder synthesized by using molten KCl-KBr salt at 750 °C showed a very weak ferrimagnetic nature at low magnetic field[32].

In order to get more information concerning the effect of KBr matrix we have measured the electroresistivity of the nanocomposites (BLFO)x-(KBr)1-x and nanograined ceramics BLFO produced by the sol-gel method. Current-voltage (I-V) characteristics of all these samples are quasi-linear (see Supplement, Fig. S5), and the values of the electro-resistivity determined from I-V characteristics are listed in Table 1. The I (V) curves and resistivity values presented in the Table 1 demonstrate the pronounced decrease of electroresistivity with an increase of KBr content. Since the I(V) curves represent only electronic conductance it is not possible to itemize a contribution of the ferrite and KBr phases into the conductivity of the composites; dc conductivity of KBr phase is caused by a number of defects formed within the KBr grains, while an increase in the conductivity of the ferrite phase can be induced by a formation of Fe2+ ions in the surface layer of the ferrite grains due to redox reactions occurred because of the KBr constituent. It should be noted that a formation of a significant amount of the Fe2+ ions should lead to a rapid decrease of coercivity of the composites[33] while the isothermal magnetization data (Fig. 2) do not show any notable decrease of coercivity.

Table 1

Electro-resistivity the nanocomposite x(BLFO)–(1-x)KBr.

	Composite xBi0.9La0.1FeO3–(1-x)KBr
Resistivity (in 109 Ω/cm)	x = 1	x = 0.9	x = 0.7	x = 0.5	x = 0
	>10	5	2	0.1	<0.01

Theoretical Explanation

Discussion of the possible mechanisms

Different physical mechanisms (analyzed in the next subsection) can be responsible for the appearance of room-temperature weak ferromagnetism of Bi1-zLazFeO3 nanosized inclusions. Recently we constructed a comprehensive Landau-Ginzburg-Devonshire (LGD) thermodynamic potential and using it modelled the phase diagram of pristine BFO[34]. The role of the AFD, rotomagnetic (RM), and rotoelectric (RE) couplings was established in ref.[35]. This in complex allows reconstructing the phase diagram of BFO and long-range order parameter distributions including the temperature stability of the AFM, FE, and AFD phases, as well as prediction of novel intermediate structural phases.

Actually, Bi1-zLazFeO3 can exhibit weak ferromagnetic properties near the surface of micro- and especially nanoparticles via structural distortions unlocing the cycloid modulated spin structure and resulting in Fe3+ spins canting due to non-vanishing Dzyaloshinskii-Moriya (DM) interaction[17,18,36-40]. Different physical mechanisms can be responsible for the appearance of room-temperature ferromagnetism of ferroelectric Bi1-zLazFeO3 inclusions, while Bi1-zLazFeO3 is FE and purely AFM in the bulk[3]. Such mechanisms, in particular, are surface piezomagnetic effect existing near the surface of any antiferromagnetic due to the absence of spatial inversion center at the surface[41], flexo-magnetoelectric[42] and linear antiferrodistortive-antiferromagnetic[43] couplings, Vegard effect (so called “chemical pressure”)[44,45] as well as omnipresent magnetic defects and/or impurities accumulation at the surface due to the strong lowering of their formation energy at the surface[46,47].

According to this, different physical mechanisms should be considered to explain the unexpected experimental result for 0.5(BLFO)–0.5(KBr) and 0.7(BLFO)–0.3(KBr) nanocomposites (the model is shown in Fig. 3), but only those among them, which are highly sensitive to the interaction of antiferromagnetic BLFO with ionic-conductor KBr, can be relevant.

Figure 3

(a) Schematics of the spherical grain with radius R covered by the shell of thickness R, where the defects are accumulated. Adapted from ref.[33].

Below we will try to prove that the appropriate mechanism turned out to be ferro-magneto-ionic coupling, that is a natural extension of ferro-ionic and anti-ferro-ionic couplings[48-51] revealed earlier in ultra-thin ferroelectric films exposed to ionic exchange with ambient media (see original papers[52,53], review[54] and Refs therein). Following Stephenson and Highland (SH) model[52,53] the screening by ions is electrically coupled to the electrochemical processes at the ferroelectric surface, and thus, the stabilization of ferroelectric state in ultrathin PbTiO3 films occurs due to the chemical switching[55-57]. The analysis[49,51] leads to the elucidation of the ferro-ionic, anti-ferro-ionic and electret-like non-ferroelectric states, which are the result of nonlinear electrostatic interaction between the ferroelectric polarization and absorbed ions.

In accordance with experiments, the critical size of ferroelectricity and antiferromagnetism disappearance in BiFeO3 nanoparticles is very small – about 10 nm[58]. So that the “core” of 100-nm Bi0.9La0.1FeO3 particle for sure should be ferroelectric (FE) and anti-ferromagnetic (AFM). The ferromagnetism can be induced by the ferro-ionic exchange only in the vicinity of surface[59,60], namely in a thin “shell” of thickness 2–5 nm (see Fig. 3).

In accordance with experimental results presented above (see Fig. 2), one can assume that in (BLFO)x–(KBr)1-x nanocomposites the concentration of the “frustrated” spins in the shell of Bi0.9La0.1FeO3 nanoparticle first increases and then saturates with the increasing fraction of KBr. Thus, the idea of ferro-magneto-ionic coupling is that the surface of ferroelectric Bi0.9La0.1FeO3 is in the dynamic direct and indirect electron exchange with ionic conductor KBr. The role of the exchanges is twofold. The first one concerns the electronic states near the ferrite surface, which can lead to a rather weak ferromagnetic exchange due to the linear ME effect at the surface (as well as to DM-interaction, and/or geometric frustration). The mechanism will not be considered below.

The second is related to the fact that the interaction with nonmagnetic KBr matrix, via the rotostrictive and magnetostrictive coupling[35], can induce the magnetic dipoles in the vicinity of surface of the Bi0.9La0.1FeO3 inclusions. At 25% concentration of ferromagnetic spins at the ferrite surface and under the surface, magnetic interaction between them in the particle “shell” cannot be excluded. In this case, the appearance of a loop with a relatively large coercive field occurs due to the mechanism of blocking the elementary spins.

Theoretical modeling

To quantify the above idea, we use the thermodynamic potential of LGD-type that describes AFM, FE, and AFD properties of BFO nanoceramics, including the RM, RE and ME biquadratic couplings, and the AFD, FE, AFM contributions, as well as elastic energy in the form[33].

Gibbs free energy of the nanoparticle (NP) has the following form:where G and G are the magnetic and magneto-elastic (magnetostriction and rotosriction) contributions to the NP energy, respectively.

The magnetic part of the free energy G describing the AFM-FM order-disorder type transition in the NP has the form:

The NP average volume is . The dimensionless order parameter is a degree of the magnetic dipoles ordering in the NP. Subscript i = 1, 2, 3 indicates Cartesian coordinates x, y and z. The value is statistically averaged over the orientations of elementary magnetic dipoles l in the NP. The concentration N of magnetic dipoles is related with the concentration of ions at the NP surface. J is the exchange constant that is positive () for the considered case of the AFM core, and related with Neel temperature T as in the mean field approximation (Boltzmann constant is ). J is the nonlinear exchange constant. The macroscopic magnetization components, , are coupled with the applied magnetic field H.

The magneto-elastic contribution to the free energy (1) is:

In Equation (3) σ is the elastic stress tensor, Z is the magnetostriction stress tensor, and s is the elastic compliances tensor of the AFM material. The summation is performed over all repeated indices. The bracket means the statistical averaging (summation) that is regarded equivalent to the intergation over the NP volume, in the ergodic case. The last term in Eq. (3) is the Vegard-type energy densitywhere we denote the elastic dipole tensor (Vegard expansion) of a surface defect as W[61]. Tensor is the magnetostriction strain tensor, c is the elastic stiffness. The value is the random concentration of elastic defects in the NP shell (including complexes with the surface ions and related with them). The approximate equality in Eq. (4) is valid if the main magnetization-dependent part of the stress is due to the magnetostriction mechanism. The function designates the small high order terms.

Allowing for the presence of the term , the energy (3) changes the coefficient in the term in Eq. (2). Thus, the substitution of elastic fields (4) into the Eq. (3) and then to the Gibbs potential Eq. (1) leads to the renormalization of the coefficient in Eq. (2), namely:

The statistical averaging over defect distribution in expression (5) gives:

The summation in Eq. (6) is performed over defect sites and the averaging of the function in the equation leads to the integration over the shell region, where the defects are accumulated. For the clarity we assume that the distribution function of defects N(r) depends on the distances from the NP surface (as the strongest inhomogeneity), and has exponential decay far from the surface [see Fig. 3]:

here is the decay length of defect concentration under the surface. The amplitude depends on the KBr fraction (1-x) of NP surrounding, and preparation conditions of the composite. Indeed the ionic surrounding affects the defect formation energies, in accordance with e.g. Stephenson-Highland ionic adsorption model[50,53]. From Eq. (7), the average concentration of defects is

The approximate equality is valid at .

The temperature of the possible FM transition induced by surface ions can be defined from the expansion of the free energy (2) in the mean-field approximation, namely from the expansion up to quadratic terms of the expression

The coefficients α and β have the following form:

The critical temperature satisfies the equation and acquires the form:

Since , the product should be negative in order to make positive under the condition . Hereinafter we assume that the Vegard tensor is isotropic and diagonal, i.e. , and magnetostriction tensor symmetry is cubic, i.e. . So that . We should underline that positive can originate from magnetostivie coupling.

For a semi-quantitative description of the magnetization curves observed in our experiments [see Fig. 2] we model the dynamics of magnetization dependence on the quasi-static magnetic field from the relaxation time-dependent equation

Hereinafter is the relaxation coefficient. In the static case the right-hand side of Eq. (11) can be expanded in l-series as . The values of the remanent magnetization and “intrinsic” thermodynamic coercive field H can be estimated from the expansion as:

The dependence of magnetization of (BLFO3)x–(KBr)1-x nanocomposite on quasi-static magnetic field was calculated from Eq. (11) at RT. Results are presented in Fig. 4(a). Different curves correspond to the gradual decrease of the BLFO fraction x = 1, 0.9, 0.8, 0.7, 0.5, 0.25, 0.1 and 0.01 (see labels near the curves). Two insets Fig. 4(b,c) illustrate the dependence of the remanent magnetization (Mr), maximal in-field magnetization (Mm) and coercive field (Hc) on the fraction 1-x of KBr.

Figure 4

(a) Magnetization dependence on quasi-static magnetic field calculated for (Bi0.9La0.1FeO3)x – (KBr)1-x nanocomposite at RT. Different curves (1–7) correspond to the gradual decrease of the Bi0.9La0.1FeO3 fraction x = 1, 0.9, 0.8, 0.7, 0.5, 0.25, 0.1 and 0.01 (see labels near the curves). Two insets (b,c) show the dependence of the remanent magnetization (Mr), maximal in-field magnetization (Mm) and coercive field (Hc) on the fraction x of KBr. The average radius R of Bi0.9La0.1FeO3 nanoparticles was 50 nm, parameters J/k = −650 K, z = 5 × 10−3 magnetic units,  = (1.4 − 1.6) × 1023 m−3 for the curves (1–7), respectively, W = −10 Å3, and R = 2 nm.

There are evident qualitative similarities (such as the loop order and shape changes with x increase) between Figs 2 and 4. However they differ quantitatively. In particular, the calculated concentration dependence of the remanent magnetization and the saturation law at high magnetic fields significantly differs from the experiment. It appears that the experimental data shown in Fig. 2 can be rigorously fitted (with not more than several % error) by the function

Here M is the magnetization, H is the applied magnetic field. The fitting parameters M0, H, Δ, and M depend on BLFO concentration x and are listed in Table S-1 in the Supplement. They are different for the “up” and “down” parts of the loops, and the small values of M can be regarded as the error. Notably that Eq. (13) contains a small linear term not included in the theoretical model (11), which magnitude depends on x% of BLFO. The presence of the linear term may be associated with the net response of the Bi0.9La0.1FeO3 ceramics, which magnetization quasi-linearly increases with H. The presence of the nonlinear term corroborates the validity of the Ising-type order-disorder Eq. (11) used for the qualitative understanding of the experimental results. Actually, the approximate analytical solution of the static equation (11) is that is rigorous at negligibly small contributions of the linear and nonlinear terms and , respectively. The nonlinear contribution to the magnetostatic response is attributed with ferro-magneto-ionic coupling between BLFO and KBr, mediated by magnetostriction.

Conclusion

We studied magnetostatic response of the composites consisting of nanosized ferrite Bi0.9La0.1FeO3 conjugated with fine grinded ionic conducting KBr powder. When the fraction of KBr is rather small (less than 15%) the magnetic response of the composite is very weak and similar to that observed for the compound Bi0.9La0.1FeO3; pure KBr matrix without Bi1-xLaxFeO3 has no magnetic response at all as anticipated. However, when the fraction of KBr increases above 15%, the magnetic response of the composite changes substantially and the field dependence of magnetization discloses ferromagnetic-like hysteresis loops with the remanent magnetization about 0.14 emu/g and the coercive field about 1.8 Tesla (at room temperature). Nothing similar to the ferromagnetic-like hysteresis loop can be observed in Bi0.9La0.1FeO3 ceramics, which magnetization quasi-linearly increases with magnetic field.

Different physical mechanisms were proposed to explain the unusual experimental results for nanocomposites Bi0.9La0.1FeO3-KBr, but only those among them, which are highly sensitive to the interaction of antiferromagnetic Bi0.9La0.1FeO3 with ionic conductor KBr, can be relevant. We have shown that the ferromagnetic behavior appears as a synergetic effect driven by the ferro-magneto-ionic and magnetostrictive coupling.

Supplement