Magnetoelectric Properties of Aurivillius-Layered Perovskites

: In the present work, we have synthesized rare-earth ion modified Bi 4 − x RE x Ti 2 Fe0 .7 Co 0.3 O 12 − δ (RE = Dy, Sm, La) multiferroic compounds by the conventional solid-state route. Analysis of X-ray diffraction by Rietveld refinement confirmed the formation of a polycrystalline orthorhombic phase. The morphological features revealed a non-uniform, randomly oriented, plate-like grain structure. The peaks evident in the Raman spectra closely corresponded to those of orthorhombic Aurivillius phases. Dielectric studies and impedance measurements were carried out. Asymmetric complex impedance spectra suggested the relaxation of charge carriers belonging to the non-Debye type and controlled by a thermally activated process. Temperature-dependent AC conductivity data showed a change of slope in the vicinity of the phase transition temperature of both magnetic and electrical coupling natures. Based on the universal law and its exponent nature, one can suppose that the conduction process is governed by a small polaron hopping mechanism but significant distortion of TiO 6 octahedral. The doping of the A-sites with rare-earth element ions and changes in the concentrations of Fe and Co ions located on the B-sites manifested themselves in saturated magnetic hysteresis loops, indicating competitive interactions between ferroelectric and canted antiferromagnetic spins. The magnetic order in the samples is attributed to pair-wise interactions between adjacent Fe 3+ –O–Fe 3+ , Co 2+/3+ –O–Co 3+/2+ , and Co 2+/3+ –O–Fe 3+ ions or Dzyaloshinskii–Moriya interactions among magnetic ions in the adjacent sub-lattices. As a result, enhanced magnetoelectric coefficients of 42.4 mV/cm-Oe, 30.3 mV/cm-Oe, and 21.6 mV/cm-Oe for Bi 4 − x Dy x Ti 2 Fe0 .7 Co 0.3 O 12 − δ (DBTFC), Bi 4 − x La x Ti 2 Fe0 .7 Co 0.3 O 12 − δ (LBTFC), and Bi 4 − x Sm x Ti 2 Fe0 .7 Co 0.3 O 12 − δ (SBTFC), respectively, have been obtained at lower magnetic fields (<3 kOe). The strong coupling of the Aurivillius compounds observed in this study is beneficial to future multiferroic applications.


Introduction
Multiferroic compounds are a fascinating type of material that have garnered remarkable attention in the realm of materials science.Multiferroics exhibit two or more orders of ferroelectric (FE), ferroelastic, and magnetic (ferromagnetism (FM)/anti-ferromagnetism (AFM)) orders in the single phase.The present Aurivillius phase materials have shown their potential usage in many novel applications, like memory storage devices, sensors, transducers, and actuators [1].Multiferroics are promising materials due to their potential applications that exploit the combination of ferroelectric and magnetic orders [1,2].The most fascinating prospect of research on these materials is tuning the electric dipoles by applying magnetic fields under superimposed small AC magnetic fields [3][4][5][6][7].The coupling between ferroelectric (FE) and ferromagnetic (FM) orders makes it easier and faster to address a data bit in memory devices [2,8].This could lead to the development of magnetoelectric materials owing to their speed in read/write operations.However, very few room-temperature magnetoelectric (ME) compounds have been reported, due to the mutual exclusion nature of FM and FE in single-phase compounds.Ferroelectricity is attributed to empty d-orbitals, while partially or half-filled d-orbitals are responsible for magnetic nature [9].In the path of searching for single-phase multiferroic compounds, Aurivillius phase multiferroics are a class of bismuth layer compound that have been extensively investigated for their flexibility and the tunable nature of their ferroelectric and magnetic properties [10][11][12][13][14][15][16][17].
The general formula of Aurivillius phase compounds can be expressed as (Bi 2 O 2 ) 2+ (A m−1 B m O 3n+1 ) 2− , where m represents the dimensionality of perovskites (A m−1 B m O 3n+1 ), which are sandwiched between fluorite-like (Bi 2 O 2 ) 2+ layers.Here, A is the cuboctahedral site, which can be substituted by mono-, di-, or tri-valent elements and B is the octahedral site (BO 6 ), which can be substituted by di-, tri-, tetra-, or pentavalent elements [18].The Aurivillius phase compounds were found to be familiar for their ferroelectric/magnetoelectric behavior, high transition temperatures, and strong anisotropic nature of spontaneous polarization [19].Rare earth-modified A-site Aurivillius phases were studied to explore the changes in structural and ferroelectric properties [20].These materials are well-suited for non-volatile memory devices due to their fatigue-free nature.Recently, many Aurivillius phase compounds have shown room temperature magnetoelectric output when Fe/Co/Ni were placed in the B-sites.However, the origin of the combined ferroelectric and magnetic and its coupling are still under discussion.Therefore, a strategic way is required in choosing magnetic ions, which triggers the enhanced magnetoelectric output in the materials.It is a known fact that the multiferroic properties can be altered by the concentration of magnetic ions in the B-sites.From previous literature, iron-doped Aurivillius phase compounds show strong multiferroic behavior at room temperature [12,21].These materials can be expressed as Bi m+1 Fe m−3 Ti 3 O 3m+3 (BFTO-m).More interestingly, the key factors of the Aurivillius phase materials are that: (i) the origin of ferroelectricity is a combination of oxygen octahedral rotation and polar distortion; (ii) different numbers of layers in perovskite plates show significant differences in physical and structural properties.
The three-layer structure of bismuth titanate (Bi 4 Ti 3 O 12 -BIT) has attracted researchers due to its high ferroelectric transition temperature (above 600 • C), low processing temperature compared to higher-layer structured compounds, and its anisotropic ferroelectric behavior [37][38][39][40][41].However, these compounds possess high conductivity due to the volatile nature of bismuth at higher temperatures.This results in leakage currents, which deteriorate the ferroelectric nature [42][43][44][45].The incorporation of rare-earth ions at the A-site has been demonstrated as an effective approach to mitigate this trade-off [13,46,47].A noteworthy observation by Lu et al. is that Fe-doped BIT ceramics exhibit weak room-temperature ferromagnetic behavior [48].Shigyo et al. reported that the La-, Nd-, and Sm-doped BIT samples show improved magnetic and ferroelectric properties due to controlled oxygen vacancies [49].According to Paul et al., Bi 4−x Sm x Ti 3−x Co x O 12−δ at x = 0.07 has shown that the coefficient of ME coupling is about 0.65 mVcm −1 [50].In recent studies, Yu et al. have reported an ME coefficient of 31.58 mVcm −1 Oe −1 for Bi 4 LaTi 3 Co 0.3 Fe 0.7 O 15 ceramics [51].Aurivillius phases possessing antiferromagnetic order or weak ferromagnetic order have limitations for their practical use [52][53][54].These defects can be controlled by RE ions for the A-sites of Aurivillius phases.It should be remembered that the stoichiometry and synthesis conditions play a vital role in optimizing the ME coupling.
The major endeavor of the present work was to diminish the leakage current by introducing RE ions in the A-sites of perovskites and enhancing the cross coupling.We prepared RE (Dy, Sm, La) ion-doped Bi 4 Ti 2 Fe 0.7 Co 0.3 O 12−δ ceramics.Electrical, magnetic, and ME measurements were carried out.The dispersion in AC-conductivity and relaxations in the low-frequency range above 200 • C were explained in terms of Maxwell-Wagner and multiple relaxations (non-Debye type).The presence of a small polaron conduction mechanism played a vital role.The dielectric relaxations observed near 200 • C can be attributed to the charge hopping of Ti 3+ to Ti 4+ /Fe 2+ to Fe 3+ ions.The low concentration doping of lanthanides enhanced the dielectric properties.Apart from the above aspects, understanding structural evolution is interesting in Aurivillius phase compounds.Systematic studies on X-ray, Raman, and impedance analysis confirm the inherited magnetoelectric properties.Remarkably, the Dy-doped Bi 4 Ti 2 Fe 0.7 Co 0.3 O 12−δ (DBTFC) sample showed an ME coefficient of ~42.4 mV/cm-Oe at room temperature.Such systematic studies have not been reported so far.were milled in a planetary mill for 24 h using zirconia balls followed by calcination at 850 • C for 4 h.The calcined powder was again grounded with an agate mortar and pestle.The powder mixed with 2 wt.% polyvinyl alcohol binder was pressed into circular pellets with a thickness of around 1 mm and a diameter of 10 mm by a uniaxial hydraulic press.The green pellets were sintered at 900 • C for 4 h to attain higher density.The RE (Dy, Sm, and La) ion-doped Bi 4 Ti 2 Fe 0.7 Co 0.3 O 12−δ compounds were named DBTFC, SBTFC, and LBTFC, respectively.

Aurivillius multiferroic ceramics Bi
Phase conformation of all the samples was performed by Powder X-ray diffraction using a PAN analytic X'Pert diffractometer (Almelo, The Netherlands).The sintered discs were crushed into powders and annealed at 400 • C to relieve mechanical strain.Rietveld refinement was also performed for all samples using Full-Prof Suit software (Version 3).The phase of all samples was well matched with standard three-layered compounds of XRD data (BIT, ICDD#00-008-0258).To study surface morphology and elemental analysis (Energy Dispersive X-ray spectroscopy, EDAX) the samples were analyzed using a ZEISS-EVO185 SEM instrument (Oberkochen, Germany).Image J 1.54h software was used to calculate the average grain size of the sintered pellets, employing histograms with Gaussian distribution.Raman spectra were measured using a Jobin Yvon spectrometer (Horiba, Japan) with a 532 nm excitation laser.The magnetization-magnetic field measurements were carried out using a vibrating sample magnetometer (VSM, Cryogenic, London, UK).ME studies were performed using a micro-measurement Group Strain Indicator, Model 3800, and a series WK strain gauge.The samples were electrically and magnetically poled before ME was performed.Ferroelectric characterization was carried out by measuring polarizationelectric field hysteresis loops under various drive voltages using a custom-made loop tracer.Before electrical measurements, the samples were poled under an electric field of 20 kV/cm by placing them in a silicon oil bath at 120 • C and cooling them in the presence of the field.Magnetic poling was performed by placing the samples in the magnetic field at 5 kOe at room temperature (RT).

Powder X-ray Diffraction and Morphological Studies
The powder X-ray diffraction patterns of DBTFC, SBTFC, and LBTFC were obtained in the range of 20-80 • at room temperature.The diffraction peaks of the compounds were in good agreement with the standard XRD pattern of Bi 4 Ti 3 O 12 (ICDD#00-008-0258) with space group Fmmm as shown in Figure 1a-c [55,56].The XRD data were fitted into Rietveld refinement with Full-Prof software (version 3).The refining analysis revealed that all samples consisted mainly of a single-phase structure with an admixture of another phase of Bi 12 Fe 0.5 Ti 0.5 O 20 belonging to the I23 space group.From the XRD studies, it could be concluded that the pyrochlore (Bi 12 FeO 20 ) formation was mainly due to the substitution of heterovalent elements (Ti 4+ , Fe 3+ , and Co 2+ ) in the B-site of the ABO 3 perovskites of the Aurivillius phases [55,57].The reliability parameters (R p , R wp , and R exp ) and goodness of fit (χ 2 ) were within the specified range and well-matched between theoretical and experimental diffraction patterns.The experimental and theoretical values of reliability (R) factors revealed the quality of the phase formation of compounds and the observed error was found to be less than 15% for higher present structured phases [58].Here, the pseudo-Voigt function was adopted to refine the XRD patterns (the original curve is a combination of both Gaussian and Lorentzian).The tolerance factor (τ) explains the perovskite nature of the samples.It is reported that τ < 1 indicates perovskite nature.From these values, one can understand the tilt of TiO 6 octahedral, which is mainly responsible for dielectric relaxations in the layered perovskite compounds.The tolerance factor for the present compounds was calculated using the following equation [59,60]: The terms R A , R B , and R O represent the ionic radii of A-site, B-site, and Oxygen ions respectively.The variation of the tolerance factor with the ratio of R Bi /R RE is shown in Figure 1d.The size of the RE ionic radii is smaller than the Bi 3+ (1.35 Å) and follows Hume-Rothery conditions [59], and therefore it is reasonable to infer that RE ions can readily be substituted with Bi 3+ sites [60].The changes observed in the tolerance factor are due to A-O and B-O bond lengths, which would be reflected in the octahedral tilt.More aspects of this are seen in the Raman spectra, near 800 cm −1 which corresponds to the stretching of TiO 6 .The lattice parameters (a, b, c), volumes are given in Table 1 and structural parameters are given in Table 2.The most intense peak was observed around 30 • and indexed as (117).The intensity of the peak associated with the (117) plane of Bi 4 Ti 3 O 12 was highest, indicating that the prepared ceramic composition conformed to the three-layered structure (m = 3).The result is consistent with common reporting of the strongest diffraction reflection, corresponding to (1, 1, 2m + 1) reflection in the Aurivillius phase three-layered compound [56].

Morphological Studies
The morphological images of the prepared ceramics were carried out by FESEM (field emission scanning electron microscope) as shown in Figure 2. Randomly oriented plate-like grains were observed in all ceramics.This kind of surface morphology is a characteristic of Aurivillius phase ceramics [60].It was also observed that all grains were non-uniform and closely packed.The experimental densities of pellets were found to be roughly 98%.Based on the mentioned values, the samples were found to be denser.The histograms of grain sizes were drawn and well-fitted by the Gaussian distribution function as shown in the inset of Figure 2. The mean grain sizes of the pellets were observed to be ~0.492,0.490, and 0.496 µm for DBTFC, SBTFC, and LBTFC, respectively.
Furthermore, grains associated with impurity phases were found in all samples, circled in the SEM pictures in Figure 2. The presence of an additional magnetic phase was obvious due to heterovalent element substitution in the B-sites of perovskites [61].EDAX (energy dispersive X-ray analysis) studies were performed on all samples to calculate the elemental composition.The experimental (measured) and theoretical (estimated) percentage of elements are presented in the inserted table of Figure 3a-c.The measured atomic percentages of the elements were consistent with the expected formula.Based on the SEM photographs, one can speculate the formation of pyrochlore or secondary phase or mixture of Fe 2+ /Co 2+ (marked as red circle in Figure 2).Similar results were more pronounced in the intergrowth Aurivillius phase compounds [62,63].It was also observed that grain size and uniformity were found to be greater in LBTFC.The appearance of lamellar morphology along with the micro-pores could not be executed, which indicated a low density (>97%) of the sample.As a result, the element composition in all samples was found to be closer to the expected values, and it can be concluded that Fe 3+ and Co 2+ ions were segregated into the grain boundary region.

Raman Spectroscopic Studies
Raman scattering spectra of the prepared ceramics are shown in Figure 4a-c to investigate lattice vibrational modes, lattice distortion, and occupancy of doped ions at the specified A-and B-sites.All observed modes of the prepared samples showed similar characteristic modes of the orthorhombic structure of three-layered Aurivillius compounds [57,62,63].The overlapped peaks were de-convoluted using Lorentzian peak fittings, for

Raman Spectroscopic Studies
Raman scattering spectra of the prepared ceramics are shown in Figure 4a-c to investigate lattice vibrational modes, lattice distortion, and occupancy of doped ions at the specified A-and B-sites.All observed modes of the prepared samples showed similar characteristic modes of the orthorhombic structure of three-layered Aurivillius compounds [57,62,63].The overlapped peaks were de-convoluted using Lorentzian peak fittings, for better identification of peak positions.The Raman modes ν 1 , ν 2 , and ν 3 below 200 cm −1 were attributed to the A-site of ABO 3 perovskite blocks' cation vibrations.Slight changes observed for SBTFC, below 200 cm −1 , were associated with the vibration of Bi 3+ ions of (Bi 2 O 2 ) 2− layer slabs.Changes observed below 200 cm −1 for SBTFC were due to strains caused in the layer-perovskite slabs.The phonon modes at high frequencies in the range of 200-800 cm −1 were attributed to the vibrational modes of BO 6 octahedral [64].The Raman modes ν 4 , ν 5 , and ν 6 (around 225 cm −1 , 260 cm −1 , and 312 cm −1 ) demonstrated the torsional bending of BO 6 octahedral.The phonon modes ν 7 , ν 8 , and ν 13 were ascribed to the stretching of the O-Ti/O-Fe/O-Co bond.The broad peaks, namely ν 7 , ν 8 , and ν 13, were used to reveal the site occupancies of doped ions in the proposed sites.The broadening of the peaks by merging ν 8 , and ν 9 for the SBTFC sample was due to the inhomogeneous distribution in the compositions of the A-and B-sites of perovskite blocks.Finally, this in turn led to lattice distortion.From this, one can say that the RE and Co/Fe ions readily occupy the A-and B-sites of perovskites [50,65].The Raman modes ν 4 , ν 5 , and ν 6 were found to vary slightly depending on the RE ion.The same phenomenon was observed in vibration modes above 200 cm −1 for SBTFC.This was once again attributed to the lattice distortion caused by Fe and Co ion doping in Ti-sites of BO 6 octahedral.The phonon modes ν 10 and ν 11 , around 730 cm −1 in all samples, were attributed to the FeO 6 and CoO 6 octahedral.Based on this, one can say that the Fe and Co ions are partially occupied by the Ti-sites of the perovskite blocks.

Dielectric Studies
The dielectric measurements of all compounds were carried out on the three samples at different frequencies and shown in Figure 5a-c.The variation of dielectric loss (tan δ) with the temperature in the range of 30 to 500 • C at 10 kHz, 20 kHz, 30 kHz, 50 kHz, and 100 kHz frequencies are shown in the inset of Figure 5a-c.The higher the temperature was, the greater the dielectric constant.A notable dielectric dispersion observed with frequency at high temperatures was due to the presence of thermally activated polarization of charge species such as space charges, charge defects, and oxygen vacancies.A slow increase of dielectric nature especially in the DBTFC sample indicated the presence of oxygen vacancies, as pointed out by many researchers.Temperature independent of the dielectric constant up to 250 • C was attributed to the hopping conduction mechanism in between the neighboring sites of Ti 3+ and Ti 4+ or Fe 2+ and Fe 3+ .Typically, the formation of dipoles with doubly ionized oxygen vacancies ( V ′′ O is less likely than singly ionized oxygen vacancy ( V ′ O .The doubly ionized oxygen vacancies (V ′′ O ) and electrons together create defective or complex dipoles, such as Fe 3+ − V ′′ O , Co 3+ − V ′′ O , and Ti 4+ − V ′′ O .The doubly ionized oxygen vacancies predominate in the conduction process over mono ionized oxygen vacancies [66][67][68].Generally, the dielectric transition temperature of ceramics depends on complex or defective ions.Since the measurements were carried out in the range of 30 to 500 • C, the transition temperatures (T C ) of all samples could not be observed.It was observed that Fe-/Co-doped Bi 4 Ti 3 O 12 compounds showed transition temperatures [67][68][69][70].
In the present analysis, SBTFC compounds showed a higher dielectric constant than DBTFC and LBTFC.This may be due to an increase in the interfacial or space charge polarization.This increase in the interfacial polarization was observed in the FESEM images (Figure 3).This reveals a decrease in grain size in SBTFC.The decrement in grain size led to an increase in the interfacial charges due to an increment in the volume fraction of grain boundaries.The same phenomenon was also applicable to the DBTFC and LBTFC samples.A similar conclusion can be inferred from the variation of dielectric loss as a function of temperature, as shown in the inset of Figure 5a-c.

Electrical Studies
AC-conductivity analysis was performed for a better understanding of the transport mechanism in the samples.Figure 6a-c depicts the AC-conductivity (σ ac ) as a function of frequency plots (f ) at various temperatures.The AC-conductivity data were obtained by considering the following expression [69].
where σ ac is AC-conductivity, and Z ′ and Z ′′ are real and imaginary values of impedance.
The terms d and S denote the thickness and area of the sample respectively.The ACconductivity data were found to obey the following Jonscher's power law [69] where σ ac is ac conductivity, σ dc is DC-conductivity, A(T) is temperature dependent coefficient, and the exponent n represents slope.The term A(T)ω n characterizes the dispersion phenomena observed in the sample, and ω is the angular frequency.The exponent n is a dimensionless quantity that generally lies between 0 and 1.This quantity depicts the degree of interaction between mobile ions and the surrounding lattice.Attempts were made to fit the AC-conductivity data with Equation ( 2), but no data were found to fit successfully throughout the frequency region.From the AC-conductivity plots, we can observe two regions, namely high-frequency and intermediate-frequency regions, as shown in Figure 6a-c.The AC-conductivity in the low frequency (region-I) is merely independent of frequency.To account for better fitting, the experimental data have been fitted with the following relation [69]: where n 1 (0 < n 1 < 1) and n 2 (1< n 2 < 2) are frequency exponents, which give slopes of corresponding regions.The exponent n 1 depicts the intermediate-frequency dispersion attributing to the ion hopping mechanism.The exponent n 2 describes high-frequency dispersion, attributing to the localized relaxation process [70][71][72].The experimental data were in good agreement with fitting data using Equation (4) as shown in Figure 6a-c.The AC-conductivity of all samples was strongly independent of frequencies above 400 • C. The AC-conductivity ≤ 10 4 Hz was independent of frequency (ω) at all temperatures and provided DC-conductivity.At a temperature below 400 • C, the AC-conductivity became strongly dependent on frequency for all samples.Strong coalitions or merging into a single curve at higher frequencies for all samples indicated the migration of oxygen vacancies.The exponent values n 1 and n 2 were plotted against temperature (Figure 6d).The frequency-independent conductivity at lower frequencies for all samples can be illustrated using the ion-jump relaxation model [73,74].The defective ions, which were oxygen vacancies, hopped between their neighborhood vacant sites of Ti 3+ and Ti 4+ .This hoping of ions may take a long time and results in a long-range transitional motion of ions.This process is called the DC-conductivity of the samples.The correlated barrier hopping model can generally explain the frequency-dependent conductivity.It was deduced from the fitting, shown in Figure 6d, that the frequency exponent n 1 decreased with increasing temperature.A cation in the lattice was hopping either a forward or backward vacant site.Thermally active oxygen vacancies yield to more dispersive nature at lower frequencies [75].Frequency exponent n illustrates the motion or localization of charge carriers in the samples.In region-I (R1), where 0 < n 1 < 1, the charge carriers can have translation motion and a sudden long-hop takes place, whereas 1 < n 2 < 2 illustrates a localized relaxation or hopping mechanism in the vicinity of lattice sites.Here, exponent n 1 values are less than 1, so the charge carriers have translational ion hopping within the ceramics [71,76,77].On the other side, the values of n 2 are greater than 1 at a higher frequency range, which illustrates the existence of a localized relaxation process.From this analysis, one can speculate that the conductivity is mainly due to short-range hopping or a kind of competitive interaction via doubly ionized oxygen vacancies.The DC-conductivity of the samples was calculated from AC-conductivity plots by extrapolating the frequency-independent term to 1 Hz. Figure 7 shows DC-conductivity as a function of temperature.The increase in conductivity with increasing temperature indicates the negative temperature coefficient of resistance behavior.The variation of DC-conductivity with temperature can explain the overall conduction (bulk) of the samples.The activation energies were calculated using the following Arrhenius relation: where σ dc is DC-conductivity, E a is activation energy, T is temperature, and K is Boltzmann constant.The variation of DC-conductivity plots reveals the conduction of a thermally activated rotation of dipoles.From the plots, activation energies were observed to be 0.931, 0.828, and 0.807 eV for DBTFC, SBTFC, and LBTFC, respectively.From the previous reports, the activation energies of Aurivillius compounds were observed to lie in the order of 0.87-1.4eV.The present activation energy values were in good agreement with previous reports [78].Aurivillius phase compounds have the volatile nature of bismuth during high sintering temperatures.To conserve charge neutrality, certain oxygen loss occurs, as per the following Kroger-Vink notations: where, V

′′
Bi is doubly ionized bismuth vacancy, e ′ is electron released, and V ′ O and V

′′
O are singly and doubly ionized oxygen vacancies, respectively.The overall conductivity can be attributed to oxygen vacancies and ion charge fluctuations.The ion charge fluctuations in the present samples can be explained as follows: From the above equations, Fe 2+ and Ti 4+ sites form dipoles with vacant defects (oxygen vacancy).These dipoles try to orient themselves using an electron-hopping mechanism.The charge carriers get trapped near localized sites and may form large polarons.The same phenomenon can apply to the Fe 2+ ions, where the conduction mechanism is attributed to small polarons.
Another aspect of the above defect mechanism can be illustrated by complex impedance curves (Cole-Cole) as shown in Figure 8a-c.The first big and second small semi-circles represent the grain (g) and grain boundary (gb) contribution of the samples.Large grain resistance (Rg) indicates defects such as oxygen vacancies and complex dipoles at grain interfaces.This interim affects the single domain or ferroelectric nature of the samples.All fitting parameters are depicted in Table 3.The inset of Figure 8a-c exhibits the room temperature hysteresis (polarization vs. electric field) loops of the DBTFC, SBTFC, and LBTFC samples under the applied strengths of 500 V/cm, 750 V/cm, and 1000 V/cm with a constant frequency of 50 Hz.The unsaturated loops revealed very high coercive and saturation fields.Generally, the Aurivillius phase compounds showed low remanent polarization and high coercive and saturation fields.In addition, the polarization of Aurivillius phase compounds depended on factors such as crystal symmetry, strain effects, phase transitions, and the competition of dopants.The Aurivillius phase compounds caused lattice distortion in the BO 6 octahedron, which affected the electric dipole moments [57].The asymmetric and unsaturated ferroelectric hysteresis (P-E) loops illustrated the leaky nature or accounted for the presence of a greater number of oxygen vacancies or defective charge carriers.It is a known fact that ferroelectric (hysteresis) nature is considered as a collective nature rather than a single ionic migration.Since the overall conductivity was due to the migration of oxygen vacancies or complex defect dipoles at lower temperatures, a saturated hysteresis loop could not be obtained.This result was an enhancement of the magnetic nature of the samples.From this Bi 3.25 La 0.75 Ti 3 O 12 impedance data observation, it is evident that complex-defect dipoles, which were accumulated at the grain interfaces, became relaxed easily during the conduction mechanism.Therefore, similar compound Sm 0.75 Bi 3.25 Ti 2.9625 V 0.03 O 12 (SBVT), rare-earth doped Aurivillius multiferroics showed wellsaturated P-E loops [59,79].This suggests that when Ti 4+ ions were replaced by Fe 3+ or Co 3+ , these defective charge complexes were trapped by the oxygen vacancies, so higher electric fields were required to activate these trapped charges to acquire saturated P-E hysteresis loops.Within the applied electric fields, all samples showed less loss nature.

Magnetic Studies
Figure 9a-c shows the variation of magnetization with the applied magnetic field (M-H loops) for all DBTFC, SBTFC, and LBTFC samples obtained at room temperature.Small hysteresis loop areas indicate a canted antiferromagnetic (AFM) nature.This nature correlates to the tilt-canted magnetic dipole moments.It should be noted that the canted nature of two sub-lattices yields a small hysteresis loop area.These materials consist of spontaneous magnetized domains.These materials possess low retentively and coercivity values even at high magnetic fields.The introduction of magnetic ions Ni/Fe/Co in the B-site of Aurivillius phases raised the spontaneous magnetization [80].The magnetization in the Co-/Fe-doped Aurivillius phase compounds main contributions, viz (i) exchange interactions between the neighboring ions like Fe 3+ −O-Co 2+ or Fe 3+ −O-Co 3+ in iron-rich regions, and (ii) tilt-canted spins or antisymmetric DM (Dzyaloshinskii-Moriya) interactions.Threelayered single-phase Aurivillius compounds with orthorhombic Fmmm structure exhibited AFM or weak ferromagnetism, which emerged from the localized magnetic (Fe/Co) rich regions.This suggests a super-exchange interaction between neighboring Fe and Co ions via oxygen vacancies [81,82].This can also be due to the long-range magnetic order in the ceramics.The structure evolution and Fe/Co-O-Co/Fe bond angles in ABO 3 perovskites in the layered compounds also affect the magnetization of the ceramics.To understand more about FM nature, the following law of approach is used [79,80]: The term A H represents the inhomogeneity of magnetization; χ P H K is the forced field-induced magnetization; and B H 2 parameters explain magneto-anisotropic nature.By considering the importance of the above LAS, Equation ( 8) is fitted in the low magnetic field regions and shown in the inset of Figure 9a-c.By using Equation ( 8), the output of fitting curves (R 2 ) was found to be ~0.99 and thus the results of fitting were highly reliable.
In the present investigation, the DBTFC sample showed an unsaturated hysteresis nature, indicating a combination of both FM and AFM [83].The formation of Bi-O-Dy bonds enhances magnetization.The Sm 3+ -doped REBTFC sample showed a higher value of magnetization.The results are consistent with our recent work on rare-earth modified Aurivillius intergrowths whose magnetization follows the similar rare-earth dependence [57].The magnetic moment per formula unit (µ f .u. ) in the Bohr Magnetron is calculated using the following formula: where M (emu/g) is measured magnetization from the fitting of the law of saturation (shown in the 4th quadrant of Figure 9); mol.weight (g) is the molecular weight of the formula; µ B. is the Bohr magnetron; 9.274 × 10 −24 J/T, N is Avogadro's number = 6.023 × 10 23 ; and the conversion factor is 10 −3 J/T.The magnetic moment per formula units (in µ B ) for DBTFC, SBTFC, and LBTFC were found to be 0.30, 0.49, and 0.21 µ B , respectively.The SBTFC sample exhibited remanent magnetization (2M r ) of about 550 × 10 −3 emu/g at room temperature, which was bigger than the earlier reported value [84,85].The inset of Figure 9a-c (left) illustrates the plots of the derivation of magnetization (dM/dH) as a function of the applied magnetic field.It can be noticed that a single sharp peak is observed in all curves.This can be attributed to the behavior of soft kind magnetism and uniform magnetic grain nature.A single sharp peak is observed in the left side inset of Figure 9a-c.From this, one can anticipate that all prepared samples were of a soft magnet.

Magnetoelectric Studies
The variation of the magnetoelectric (ME) coefficient with the applied magnetic field of prepared samples is shown in Figure 10.The values of the magnetoelectric coefficient were found to be 42.4 mV/cm-Oe, 30.3 mV/cm-Oe, and 21.6 mV/cm-Oe for DBTFC, LBTFC, and SBTFC, respectively.[64,86].The obtained ME values of the prepared samples were found to be higher compared to the other Aurivillius compounds.Furthermore, the high ME coefficient was obtained at lower magnetic fields (3 kOe), as shown in Figure 10.The ME coupling mainly arises from two aspects, namely, (i) spinexchange or spin-orbit interactions, driven by inverse DM interaction, and (ii) spin-lattice interaction.The contribution of spin-orbit or spin-exchange interaction via inverse DM interactions is smaller than the spin-lattice interactions [87].In particular, the magnetic ions in BO 6 octahedral sites are slightly shifted from their regular sites under the application of a magnetic field.This induces a strain in the lattice structure and develops the voltage via ferroelectric accumulated changes.However, further studies are needed to establish the plausible reason for obtaining a high ME coefficient.As it is reported, the impurity phases of Co/Fe compositions influence the ME nature in Aurivillius compounds [57].The presence of circled regions of SEM micrograph dots (nanoregions) represents a Co/Fe-rich magnetic phase.Based on these observations, such Co/Fe-rich phases could contribute partly to the enhanced ME properties of the prepared samples.
In particular, when an external magnetic field is applied, the magnetic ions are slightly displaced from their normal octahedral position, which induces a voltage in the lattice, resulting in ferroelectric and ferroelastic changes, and finally producing an output voltage.The parameters of orthorhombic distortion, tetragonal stain, and orthorhombicity play a significant role in ME properties [88].In our study, the DBTFC sample showed higher ME coefficients.This is in good agreement with the lattice structure of the prepared samples.Higher orthorhombicity and lower orthorhombic distortion values were observed for the Dy-doped sample (DBTFC), which are favorable to enhancing the ME coefficient under magnetic fields.In this study, the RE (Dy, La, and Sm) element substitution in the A-site of Aurivillius compounds induced lattice distortion, and it led to accompanying electric polarization under magnetic fields.It can be concluded that the RE-doped three-layered Aurivillius phase compounds show strong ME coupling.These multiferroic materials are useful for understanding the ME phenomenon.

Conclusions
In summary, we investigated the structural, morphological, electrical, magnetic, and magnetoelectric properties of the RE-doped (Dy, Sm, and La) Bi 3 RETi 2 Fe 0.7 Co 0.3 O 12−δ Aurivillius multiferroic compounds.The XRD and Raman spectroscopic studies revealed that the prepared samples were formed in single phase and orthorhombic structure with space group Fmmm.This revealed that the ions were properly substituted into lattice cells.The plate-like, anisotropic, randomly orientated grains were observed in FESEM studies, which is a characteristic feature of Aurivillius phase ceramics.The substitution of RE significantly enhanced the dielectric and ferroelectric properties with improved dielectric constants in the frequency range, and reduced dielectric loss values were found at lower frequency ranges.Improved remanent polarization was observed in all samples at room temperature.This is attributed to the reduction in the oxygen vacancies by doping RE in the Bi-sites of the samples.Ferroelectric and magnetic studies revealed good multiferroic behavior in the samples.The M-H loops of the samples suggested an exchange interaction between adjacent ions or DM interactions among magnetic ions.Strong ferroelectricmagnetic coupling accompanied by a high coefficient of magnetoelectric coupling of 42.4 mV/cm-Oe was exhibited by the DBTFC sample.This strong ME coefficient in a single phase may be applicable in multiferroic devices and can be interpreted from the perspective of previous studies and of the working hypotheses.

Figure 10 .
Figure 10.ME coefficient vs. applied magnetic field plots of DBTFC, SBTFC, and LBTFC samples.It is observed from the previous reports that the Bi 4 Ti 3 Fe 0.7 Co 0.3 O 15 (BFTO) sample showed an ME coefficient of 16.45 mV cm −1 Oe −1[64,86].The obtained ME values of the prepared samples were found to be higher compared to the other Aurivillius compounds.Furthermore, the high ME coefficient was obtained at lower magnetic fields (3 kOe), as shown in Figure10.The ME coupling mainly arises from two aspects, namely, (i) spinexchange or spin-orbit interactions, driven by inverse DM interaction, and (ii) spin-lattice interaction.The contribution of spin-orbit or spin-exchange interaction via inverse DM interactions is smaller than the spin-lattice interactions[87].In particular, the magnetic ions 4-x RE x Ti 2 Fe0 .7 Co 0.3 O 12−δ (RE = Dy, Sm, La) were synthesized by conventional solid-state reaction method.For this purpose, stoichiometric amounts of reagents (purity > 99.9%) Bi 2 O 3 , Sm 2 O 3 , Dy 2 O 3 , La 2 O 3 , TiO 2 , Fe 2 O 3 , and Co 3 O 4