Peculiar magnetism of BiFeO3 nanoparticles with size approaching the period of the spiral spin structure.

Size effect of multiferroics is important for its potential applications in new type miniaturized multifunctional devices and thus has been widely studied. However, is there special size effect in the materials with spiral modulated spin structure (such as BiFeO3)? It is still an issue to be investigated. In this report, structural, magnetic and magnetoelectric coupling properties are investigated for sol-gel prepared BiFeO3 nanoparticles with various sizes. It is found that a structural anomaly arises for the particles with size close to the 62 nm period of the spiral modulated spin structure, which induces an obviously increased ferromagnetism. In addition, large magnetoelectric coupling effect is observed in 62 nm BiFeO3 nanoparticles. Our result provides another insight into the size effect of BiFeO3, and also a clue to the magnetic structure at nanoscale.

Multiferroics are promising materials for the design and synthesis of multifunctional devices1234. They are noteworthy for their unique and strong coupling of electric, magnetic, and structural order parameters, giving rise to simultaneous (anti)ferroelectricity, (anti)ferromagnetism, and/or ferroelasticity56. As one of the most important single-phase multiferroic materials, BiFeO3 (BFO) has attracted much attention because it has both high ferroelectric Curie temperature of about 1103 K and antiferromagnetic Néel temperature of about 643 K7. The magnetic structure of BFO is of G-type antiferromagnetic, and the combined action of exchange and spin-orbit interactions produces spin canting away from perfect antiferromagnetic ordering, resulting in a spiral spin arrangement with a wavelength of about 62 nm, and thereby producing a helimagnetic order and a vanishing magnetization in the bulk89. Like many other materials, BFO shows a variation of properties with the decrease of size1011. Néel and Curie temperature have been shown to decrease with decreasing particle size1213. Enhanced magnetization has been reported in nanoparticles, nanowires and nanopowders, which is thought to originate from the uncompensated spins at surface14151617, an effect that is well known in classic antiferromagnets such as NiO18. Moreover, some other distinct phenomena in BFO are also associated with size effect. For instance, some researchers report the fading away of the magnetic anomaly at about 50 K with the decrease of size, while some others provide just the opposite observations151920. Correspondingly, the physical mechanism for this magnetic anomaly is still an controversial issue at present1519202122. In our earlier work, an interesting film thickness dependent magnetization oscillation with a period of about 62 nm has been observed23.

All these results suggest the exceptional and complex properties of BFO in nanoscale, and the special magnetic fluctuating behavior observed in our BFO films reminds us of the important correlations among crystal structure, macroscopic magnetic property and the unique spiral spin structure. However, in spite of intense studies about BFO, it is still unknown what will happen when the sample size is close to the spiral spin period of BFO. Thus, in this paper, we systematically investigate the magnetic properties of BFO particles with size from 18 nm to 248 nm. Similar to other reports, it is observed that the ferromagnetic property basically increases with decreasing particle size because of the increased surface to volume ratio. However, our most important result is that BFO particles with size close to the spiral spin period 62 nm possess exceptionally enhanced ferromagnetic property and magnetodielectric effect compared with others, which is considered to be closely related to the structural anomaly at this particle size, and is also possibly associated with its intrinsic spiral-modulated spin order.

Results

Figure 1(a) shows the X-ray diffraction (XRD) patterns of BFO nanoparticles annealed at different temperatures, where the peaks are indexed based on Rhombohedral structure (JCPDS File No. 71-2494). It is seen that all the nanoparticles present single-phase perovskite structure and no obvious impurity phase is observed within the resolution of XRD. Fig. 1(b) and supplementary Fig. S1 show that the full width at half maximum of (012) peak decreases as the annealing temperature increases from 425 to 575°C, indicating the increase of particle size. The average particle sizes (i.e. diameter of the particle, d) determined from scanning electron microscopy (SEM) (Supplementary Figs. S1–S3 and Table S1) are around 18 ± 5 nm, 31 ± 7 nm, 45 ± 7 nm, 62 ± 12 nm, 83 ± 20 nm, 170 ± 45 nm, 248 ± 100 nm for the BFO nanoparticles annealed at 425°C, 450°C, 475°C, 500°C, 525°C, 550°C, and 575°C, respectively. Typical TEM images and high resolution TEM images obtained from a portion of an individual 62 nm (Fig. 1 (c)–(e)) and 18 nm (Supplementary Fig. S4) particle further confirm the highly crystalline nature of our BFO nanoparticles. Size-dependent crystal structure of BFO nanoparticles is analyzed by Rietveld refinement based on hexagonal representation of space group R3c, and the results are summarized in Fig. 1 (f), where a and c are the lattice parameters in hexagonal representation. No impurity phase is found in the refinement process, and the small values of R and R indicate that all the BFO nanoparticles fit very well to the Rhombohedral structure.

Figure 1

(a) X-ray diffraction patterns of BFO nanoparticles annealed at various temperatures. (b) Annealing temperature dependent full width at half maximum (FWHM) of (012) peak. (c) The TEM image and (d–e) high resolution TEM images of the 62 nm BFO particles. (f) Parameters obtained from Rietveld refinement.

Figure 2 (a) shows the initial magnetization curves measured at room temperature. Nonlinear curve is only observed for the 62 nm BFO particles at the magnetic field less than 5 kOe, indicating its different magnetic structure compared with that of all the others, since the shape of initial magnetization curve is controlled by domain structure24. This is further confirmed by magnetization-field hysteresis measurements. Typical antiferromagnetic response similar to bulk BFO is obtained for the particles with diameters of about 170 nm and 248 nm (Supplementary Fig. S5), while ferromagnetic hysteresis loops are observed for the particles with diameter smaller than 83 nm (Fig. 2(b)). Moreover, as shown in the inset of Fig. 2 (b), except the abnormal point for the 62 nm particles, the magnetization measured at 7 T increases with the decreasing particle size and presents almost a linear dependence on 1/d (or the surface to volume ratio), indicating that these BFO nanoparticles can be modeled by a superposition of an antiferromagnetic core and a ferromagnetic surface, and their magnetization are mainly contributed by the uncompensated spin at grain surface1518. The large deviation from the expected behavior for 62 nm BFO particles reveals the presence of new contribution to its ferromagnetism in addition to the uncompensated spin at the grain surfaces. In order to get a clearer image of ferromagnetic contribution in the hysteresis curves, the linear contributions from antiferromagnetic and paramagnetic parts are deducted, and the results are shown in Fig. 2 (c). Compared with others, obviously increased saturation magnetization (M) while decreased coercive field (H) (inset of Fig. 2(c)) are observed for the 62 nm BFO particles, which further approves the distinct origin of its ferromagnetism. We would like to mention that, all the results shown in this paper can be fully repeated for different batches of our sample.

Figure 2

(a) Initial magnetization curves, (b) room temperature magnetic hysteresis loops, (c) the ferromagnetic parts in the magnetic hysteresis curves of BFO nanoparticles with various sizes. The inset in (b) gives the magnetization at 7 T as a function of d and 1/d, and in (c) shows size-dependent M and H.

Electron energy loss spectroscopy (EELS) inspection may provide the clue for the abnormal magnetism of 62 nm particles. As shown in Fig. 3 (a) and in supplementary Fig. S6, two typical peaks of oxygen K-edge located at about 531 eV and 540 eV are observed for all the particles. While interestingly, obvious furcation of the two peaks are observed only for the 62 nm BFO particles, indicating the changed bonding information of oxygen25 and implying that a structural anomaly of a so-far unknown type is present for the particle. Since it cannot be observed in XRD measurement and refinement, we argue that the structural anomaly is merely an enhanced distortion without symmetry change, such as the increased rotation distortion of FeO6 octahedron around the [111] direction, which can cause the strengthening of Dzyaloshinskii-Moriya (DM) interaction and the suppression of the spiral spin structure262728, and thus increase the ferromagnetism of the 62 nm BFO particles. Theoretical results further confirm this inference. First, dependence of spontaneous magnetization on the rotation angle of FeO6 octahedron around [111] direction is analyzed by first principles calculation29. We fix the positions of cations and change the coordinates of oxygen ions to obtain different rotation of FeO6 octahedron, and then calculate the spontaneous magnetization under different rotation angles. Fig. 3 (b) shows that the magnetization per unit cell increases with the increased angle, indicating that local ferromagnetism of BFO can be increased by the increased rotation distortion of FeO6 octahedron through the strengthened DM interaction. Second, based on the spherical particle model, the macroscopic magnetization are calculated as a function of particle size according to the method given in our former paper29. For the 62 nm BFO particles, the ferromagnetism of DM type is calculated by first principles calculation as mentioned above supposing the rotation angle increases by 2° from 12.8° – that reported for bulk BFO ceramics at room temperature30. As shown in Fig. 3 (c), obviously increased magnetization are also obtained for the 62 nm BFO particles, which is in reasonable agreement with the experimental behavior given in the inset of Fig. 2(c). The relatively larger calculated magnetization for all the particles compared with our experiments probably comes from the deviation of the rotation angle, the randomly oriented grains and the irregular particle shape. In addition, the enhanced DM interaction can also induce improved ferroelectric property31, which is always reflected by increased value, where anor and cnor are the lattice parameters of tetragonal structure calculated from , and 32. Figure 3(d) depicts the variation of with particle size. Similar to the result given by Selbach et al.12, with the decrease of particle size, the value first decreases slowly and then abruptly with 30 nm as an inflexion point. However, an abnormal point corresponding to an obviously increased is observed for the 62 nm BFO particles, indicating its ferroelectric property is improved probably by the enhanced DM interaction through magnetoelecric coupling. Taken together, a structural anomaly probably corresponding to the increased rotation distortion of FeO6 octahedron is present when the particle size approaches the period of the spiral spin structure. This structural anomaly could cause an increased DM interaction, and in turn a suppressed spiral spin structure, which finally result in an increased ferromagnetism of BFO particles with size around 62 nm. While BFO nanoparticles with other sizes can be modeled by a superposition of an antiferromagnetic core and a ferromagnetic surface, as a result, the ferromagnetism contributed by the uncompensated spin at grain surface increases as particle size decreases.

Figure 3

(a) Electron energy loss spectroscopy of oxygen K edges for 31 nm, 62 nm and 170 nm BFO particles. (b) Dependence of magnetization per unit cell on the rotation angle of FeO6 octahedron. (c) Size-dependent spontaneous magnetization. (d) Size-dependent .

Besides the room temperature magnetic property, the particularity of 62 nm particles can also be gleaned from temperature dependent magnetization studies. Figure 4 shows the magnetization of our BFO nanoparticles as a function of temperature, measured at 200 Oe after zero-field cooling (ZFC) and also with 200 Oe field cooling (FC). There are three features worth mentioning: (1) The splitting in the ZFC and FC magnetization curves are observed for all the BFO nanoparticles. (2) With temperature decreasing from 300 K to 10 K, the ZFC magnetization value (M) monotonically decreases for the 170 nm and 83 nm BFO particles, while it first increases and then decreases, as a result, a broad peak appears in the ZFC curves for the BFO particles with size less than 62 nm. (3) Except for the absence in 62 nm particles, a sharp cusp (shown in Figure 4 by arrows) is observed in the ZFC and FC curves for other BFO nanoparticles. Here we try to explain the above features one by one.

Figure 4

Magnetization as a function of temperature measured under ZFC and 200 Oe FC conditions for BFO nanoparticles with various sizes.

1) The splitting of the ZFC and FC curves usually appears in the co-existent system of antiferromagnetic and ferromagnetic phases33, which coincides with the core-shell structure proposed above, that is, our BFO nanoparticles are composed of an antiferromagnetic core and a ferromagnetic surface. Though the spiral-modulated spin structure is probably suppressed for the 62 nm BFO particles by the increased rotation distortion of FeO6 octahedron, local antiferromagnetic interaction still exists in the particle. 2) The ferromagnetic behavior of the BFO particles with size less than 62 nm is responsible for the increased M with decreasing temperature, and the presence of the broad peak in ZFC curves can be attributed to the magnetic blocking mechanism induced by the competition between thermal energy and magnetic anisotropy energy. 3) The sharp cusp presented around 55 K (T) in ZFC and FC curves has been attributed to superparamagnetic relaxation process19 and the spin glass behavior20. If the cusp represents a blocking behavior associated with superparamagnetic relaxation processes of isolated magnetic nanoparticles, one would expect a stronger dependence of T on particle size34. However, as shown in Fig. 4, T is about 55 K for most of our BFO particles with different sizes (31 nm, 45 nm, 83 nm and 170 nm), and even for 248 nm particles (Supplementary Fig. S7) and bulk BFO ceramics15. On the other hand, Vijayanand et al. proved that the above-mentioned cusp is not associated with any spin glass or phase transition by a series of ZFC and FC magnetization measurement, and they believe that it corresponds to some changes associated with domain structure21. Here, we consider that the sharp cusp of our BFO nanoparticles may be related to the magnetic domain pinning effect, which is well known in magnetic materials35. However, why the cusp could not be observed in the 62 nm BFO particles? In order to answer this question, a series of ZFC and FC curves are measured under various magnetic fields. As shown in Fig. 5, the cusp at about 55 K can also be observed for the 62 nm BFO particles when the applied magnetic field is larger than 2 kOe, where most of the ferromagnetic domains are switched to the direction of the applied magnetic field. Therefore, the cusp at about 55 K, which can be observed in BFO nanoparticles, ceramics, and single crystals, is related to the antiferromagnetic domain pinning effect. The origin of the domain pining is usually ascribed to the non-uniformly distributed defects, which affects the movement of domain walls. As shown in Fig. 4, although the shape of the cusp is difficult to be distinguished from the ZFC curves because of the superimposition of the broad peak originated from magnetic blocking, it can be clearly seen from the FC curves that the cusp height is almost the same for the BFO nanoparticles with various sizes, indicating their similar pining behavior. As for the 18 nm BFO particles, lower T of about 35 K is obtained probably because of the possible interparticle interactions.

Figure 5

ZFC and FC magnetization curves measured with various magnetic fields for the 62 nm BFO particles.

Raman spectra are performed and shown in Fig. 6(a) since it can give some indications of spin-phonon coupling36. Similar to the literatures1537, peaks at about 136 cm−1 (peak A), 168 cm−1 (peak B) and 218 cm−1 (peak C) are observed, which can be attributed to the manifestation of A1 modes. Since the suppressed intensity of peak A implies the enhanced coupling of magnetic, ferroelectric and/or structural order parameters15, we fit the Raman spectra using WiRE2.0 software and estimate the suppressed degree of peak A by IA/IB, where IA and IB represent the integral intensities of peak A and peak B, respectively. It is seen from Fig. 6(a) that obviously decreased IA/IB values are observed for the BFO particles smaller than 62 nm, which coincides well with their enhanced ferromagnetic properties. Furthermore, the smallest IA/IB value of the 62 nm BFO particles probably indicates the increased magnetoelectric coupling in this special particle.

Figure 6

(a) Raman spectra of BFO nanoparticles with various sizes, in which the black lines and color lines correspond to the measured and fitted results, respectively. Temperature dependent dielectric constant under zero and 9.5 kOe magnetic field for the (b) 62 nm, (c) 31 nm and 18 nm particles. The inset in (b) gives the temperature dependent dielectric constant of the 62 nm BFO particles measured at different frequencies. (d) Temperature dependent MD effect for the BFO nanoparticles with various sizes.

To gain insight into the magnetoelectric coupling of our BFO nanoparticles, magnetodielectric (MD) effect, i.e the variation in dielectric constant ε under an applied bias magnetic field (, where εr(H) and εr(0) are the dielectric constant measured under and without applied magnetic field, respectively), is investigated. It is said that the MD effect can be attributed by not only the intrinsic coupling of spontaneous polarization and magnetization, but also the extrinsic coupling contributed by a combination of magnetoresistance and Maxwell-Wagner effect3839, and the MD effect measured below 700 kHz cannot truly reflect the intrinsic magnetoelectric coupling383940. Thus, the MD effect present in Fig. 6 (b–d) and Fig. S8 (Supplementary information) of our BFO nanoparticles are measured at 1 MHz for the temperature range from 4 K to 310 K. In fact, as shown in the inset of Fig. 6 (b), almost frequency independent dielectric constant is observed for the 62 nm BFO particles when the frequency range from 104.5 to 106 Hz, indicating that the MD effect of our BFO nanoparticles are mainly come from the intrinsic coupling between spontaneous polarization and magnetization. The MD effect over the whole temperature range is clearly demonstrated for the 62 nm, 31 nm and 18 nm BFO particles (Fig. 6 (b) and (c)), which is much obvious than that observed in other particles (Supplementary Fig. S8). More clearly, as shown in Fig. 6 (d), larger MD values are obtained in the temperature range 4–310 K for the 62 nm BFO particles, which is much larger than that of the nanoparticles with other size and that reported for some doped BFO materials4041. The greater MD effect of 62 nm BFO particles corresponds to its improved ferromagnetic, ferroelectric and magnetoelectric coupling properties probably related with the enhanced rotation distortion of FeO6 octahedron. In addition, the dielectric anomaly observed at around 250 K for our BFO nanoparticles (shown in Fig. 6 (b) and (c) by arrows, and also in supplementary Fig. S9) could be ascribed to surface effect according to the recently study reported by Jarrier et al.42 since no magnetic anomaly is observed in neutron data for BFO single crystal below room temperature4344.

Discussion

Up to now, most studies on the size effect of ferroelectric and magnetic materials have revealed that their properties usually vary monotonically with decreasing size. For example, the ferromagnetism contributed by the uncompensated spins at grain surface is enhanced with decreasing sample size in antiferromagnetic materials, which is also observed in our BFO nanoparticles. However, we find some unexpected changes of the properties when the particle size is around 62 nm, such as the obviously increased ferromagnetism and magnetoelectric coupling, which are related with a structural anomaly probably being an enhanced rotation distortion of FeO6 octahedron. This report reveals for the first time that, special distortion of lattice and magnetic structure can be induced when the sample size is close to the spiral period, which can further lead to extraordinary macroscopic properties. On the one hand, the results provide novel approaches to improve the ferromagnetism and magnetoelectric coupling of BFO, which is important for its applications since the antiferromagnetic nature and poor coupling between electric and magnetic properties are the primary weakness of BFO as room temperature multiferroics. On the other hand, our work not only highlights the important regulatory role of the spiral spin period on the properties of BFO but also calls forth more attention on the size effect of materials with periodically modulated spin structure. Further studies both theoretical and experimental are urgently demanded in order to understand what really happens when the sample size approaches the period of the spin-structure.

In summary, BFO nanoparticles with various sizes are prepared by a sol-gel method. For the particles with size close to the spiral spin period of BFO, a structural anomaly arises, which improves the multiferroic and magnetoelectric coupling properties of the material. The present work highlights the role of the spiral spin structure on the modulation of microscopic properties, and opens a new window for the application of BFO.

Methods

Sample preparation

BFO nanoparticles with size from 18 to 248 nm were prepared by a sol-gel method. The precursor solution was prepared by dissolving 10 mmol Bi(NO3)3·5H2O and 10 mmol Fe(NO3)3·9H2O in distilled water with the addition of dilute nitric acid, wherein tartaric acid as organic complex was added with 1:1 molar ratio to the cations. The mixture was stirred for several hours, and then evaporated and dried at 140°C for 15 h to obtain the xerogel powder. Finally, the xerogel were calcined at temperatures varied from 425 to 575°C for 2 h to obtain BFO nanoparticles with various sizes.

Structure characterization

X-ray diffraction with Cu kα radiation was used for analyzing the structure of the BFO nanoparticles. Scanning electron microscopy, transmission electron microscopy (FEI Tecnai F20) were used for analyzing the size and microstructure of the particles. Raman spectra were performed on micro-Raman spectrometer (Horiba Jobin Yvon HR800). The electron energy loss spectroscopy (EELS) were conducted on a FEI TITAN 80–300 field emission TEM equipped with a Gatan imaging filter system (model: Tridiem 863). The microscope was operated at an accelerating voltage of 300 kV, and the spectra were acquired in diffraction mode.

Magnetic and dielectric measurements

The magnetic properties were measured with a commercial magnetic property measurement system. For the measurement of dielectric properties, the nanoparticles were pressed into a tetragonal shape (3 × 3 × 0.3 mm3) at 10 MPa, then, silver paste was painted on both sides of the sample as a capacitance. The dielectric properties as well as the MD effect were investigated using HP4194A impedance analyzer with home-built sample holder coupled with a helium closed cycle refrigerator (temperature range from 4 K to 310 K) and applied magnetic field.

First-principles calculation

Our first-principles calculations were performed using the projector-augmented wave (PAW) method within the local spin-density approximation plus the on-site repulsion (LSDA + U), as implemented in the Vienna ab initio simulation package (VASP). We used U = 4 eV on Fe 3d states and a plane-wave cutoff energy of 500 eV. The calculations were performed with a 5 × 5 × 5 k-point sampling. We switched on the spin-orbit coupling (SOC) and employed the convergence threshold of 10−7 eV due to the small value of SOC. For the calculations of spontaneous magnetization, we set the initial spins collinear, and allow the spins relax freely within the self-consistency cycle to induce canted magnetic moment.

Author Contributions

F.Z.H. and X.M.L. designed the experiments and wrote the paper. F.Z.H. and Z.J.W. prepared the sample and performed the structural measurements. F.Z.H. and K.L.M. contributed to the measurements of dielectric and MD effect. R.X.T., T.T.X., J.H. and C.Y. performed the magnetic measurements. J.T.Z. conducted the calculations. F.Z.H., X.M.L., W.W.L., J.T.Z. and J.S.Z. analyzed and discussed the data. All authors contributed through scientific discussions.