Anisotropy of Percolation Threshold of BaTiO3-Ni0.5Zn0.5Fe2O4 Composite Films

A systematic study on the magnetic and electrical percolation phenomena of BaTiO3 (BTO)-NiZnF2O4 (NZFO) composite films is presented in this work with the purpose of simplifying the preparation process of high-performance 1–3-type multiferroic composite films. Results show that the percolation threshold of the composite films depends on the macroscopic dimension of the material. The low-dimensional nature of the composite films results in different percolation thresholds with topological transition in vertical and horizontal directions. BTO-NZFO composite films with a grain size of 15 nm and a thickness of 100 nm exhibited a percolation threshold of 0.18 in the normal direction and a percolation threshold of 0.48 in the horizontal direction. In light of this intriguing feature, a novel multiferroic composite film with 1–3 structure and strong magnetoelectric coupling was easily prepared by a 0–3 process via controlling the NZFO content in the region between two percolation thresholds.

www.nature.com/scientificreports www.nature.com/scientificreports/ the ferroelectric-ferromagnetic composite films lies in the region between the two percolation thresholds of 2D and 3D systems.
Nevertheless, this is only a theoretical assumption that needs to be verified experimentally. Subject to limited observation area of transmission electron microscopes and thickness of samples, it is difficult to directly observe the percolation channel that normally connects the entire film. Therefore, a method to demonstrate the anisotropy of multiferroic film percolation thresholds in different directions is also of vital importance and practical significance.
Herein, a classical multiferroic film system, BaTiO 3 (BTO)-Ni 0.5 Zn 0.5 F 2 O 4 (NZFO), is utilized to verify the proposed strategy. The films were prepared by magnetron sputtering, which is a conventional process for 0-3 composites. The percolation behavior of the composite films was investigated by analyzing the electrical and magnetic properties under external fields with different directions. The anisotropy of the percolation threshold of the multiferroic composite thin films was systematically characterized for the first time. A novel way of preparing multiferroic composites with excellent performance by a simple process was verified. Figure 2(a,b) show the permittivity and dc conductivity of the BTO-NZFO composite films under out-of-plane and in-plane electric fields as a function of the volume fraction of NZFO f NZFO , respectively. As shown in Fig. 2(a), the permittivity of the composite film increases nonlinearly from 100 to 200 as the value of f NZFO increases, regardless of the direction of the electric field. However, for the composite films with the ferrite fraction lying between 0.2 and 0.8, the permittivity under an out-of-plane electric field is slightly higher than that under an in-plane field. It is known that the permittivity of the BTO-NZFO composite films conforms with the Kirkpatrick model 23 . In Fig. 2(a), it can be seen that the experimental values of the permittivity of the composite films under the out-of-plane electric field can be well fitted by the theoretical curve when the percolation threshold f c is set as 0.18, and the theoretical curve is consistent with the experimental results under the in-plane field when f c is set as 0.5. In addition, as shown in Fig. 2(b), the conductivities of the composite films under the both electric fields exhibit a leap at a certain point with increasing f NZFO . Below and above this point, the conductivity values remain almost the same in both cases. It is noteworthy that the abrupt-change point of the conductivity of the composite film under the out-of-plane electric field is f NZFO = 0.16, while that under the in-plane field appears f NZFO = 0.42.

Results
The variation in the conductivity of the BTO-NZFO composite films versus NZFO content shown in Fig. 2(b) is a typical evidence of conductivity percolation that usually occurs in the composites composed of the constituents with different conductivities. The abrupt-change point of the conductivity represents the value of f c , which is the ferrite volume fraction at the time of the formation of conductive paths in the direction of the electric field. It is well known that under an electric field, the free charges inside the material will move to generate current.  www.nature.com/scientificreports www.nature.com/scientificreports/ The direction of the charge migration is from one electrode to the other following the minimum voltage drop principle; that is, the migration paths are always along the paths with smallest total resistance under a constant applied voltage. The results in Fig. 2(b) show that when f NZFO = 0.16, a path consisted of low resistivity phases forms out-of-plane, while it forms in-plane when f NZFO = 0.42 comparatively. In BTO-NZFO composites, the conductivity of the NZFO phase is much larger than that of the BTO phase. Therefore, the formation of the low-resistance paths means that the NZFO grains between the positive and negative electrodes are able to connect to form percolation channels. When the NZFO volume fraction lies between 0.2 and 0.4, the NZFO grains fully connect out-of-plane, despite they are not yet fully connected in-plane at this time so that the charges cannot find a path with less voltage drop for long-range migration.  Fig. 3(b), the values of μ i under vertical and parallel magnetic fields both increase nonlinearly from 2.9 to 25.4 with the increase in NZFO ferrite volume fraction. However, when the ferrite volume fraction lies between 0.2 and 0.8, the initial permeability μ i of the composite films under an in-plane magnetic field is higher than that under an out-of-plane field. Figure 4(a,b) show the field-cooled and zero-field-cooled (FC-ZFC) curves of the BTO-NZFO composite films with f NZFO = 0.329, 0.433, 0.534, 0.728 and 0.912 under vertical and parallel magnetic fields, respectively. The plots of the variation of the average blocking temperature, T blo , and its difference with the irreversible temperature T irr , i.e. T irr − T blo , versus f NZFO under different magnetic fields are shown in Fig. 4(c,d), respectively. As shown in Fig. 4(c), under an out-of-plane magnetic field, T blo of the composite films first increases from 58.4 K to 72.5 K as f NZFO increases from 0.3 to 0.5, and then begins to decrease. By contrast, the value of T blo first decreases from 60.6 K to 54.3 K with f NZFO increasing from 0.3 to 0.4, and then increases until a maximum value of 76.5 K is reached at f NZFO = 0.728 under an in-plane magnetic field. As shown in Fig. 4(d), as the ferrite content increases, the value of T irr − T blo under the vertical external magnetic field remains relatively stable first and then decreases from f NZFO = 0.433 to f NZFO = 0.728, and finally remains unchanged. The value of T irr − T blo under a parallel external magnetic field decreases linearly until f NZFO = 0.728, and then remains almost unchanged.
It has been documented that in BTO-NZFO composite films, the exchange between the ferrite nanocrystals and the magnitude of induced magnetic flux can change due to the topological changes of the composites, resulting in a percolation transition of the initial permeability and coercivity in such multiferroic systems 23 . When the ferrite content in the composite films is lower than the percolation threshold, the magnetic NZFO grains are dispersed in the non-magnetic BTO matrix. The transmission of the magnetic induction current between the magnetic particles is blocked by the nonmagnetic BTO particles, leading to a small value of the initial permeability of the composite due to the suppression of the equivalent demagnetizing field. Moreover, as there is no interaction occurring between the grains, the coercivity and the anisotropy field of the composite are the average of all magnetic grains. The variation in coercivity and anisotropy field with changing composition of the composite is actually controlled by grain size. Owing to the diffusion during the crystallization process, the grain size of NZFO in the composite films generally increases with increasing ferrite content. Increased grain size renders the average anisotropic energy (represented by T blo ) and the coercivity of the material to decrease with increasing f NZFO , while the difference between the highest anisotropic energy and the average anisotropic energy (represented by the value of T irr − T blo ) does not change.
When the ferrite content in the composite films reaches the percolation threshold, the magnetic NZFO grains are connected with each other, forming a magnetic flux channel and initiating the exchange effect. The www.nature.com/scientificreports www.nature.com/scientificreports/ formation of the magnetic flux channel causes the initial permeability of the composite to be less suppressed by the equivalent demagnetizing field, and such effect accounts for the nonlinear increase in initial permeability with increasing ferrite content. The exchange effect in a material also tends to impel the coercivity to exhibit the same value with that of the harder interconnected grains. Before NZFO is completely interconnected and a new matrix phase comes into being, as the ferrite content increases, the number of interconnected NZFO grains increases, and the magnetization of hard magnetic grains can be transferred to other grains. The coercivity of the average anisotropic energy of the composite is consequently increased, and the difference between the highest anisotropy energy and the average anisotropy energy is reduced. When the NZFO grains in the composite films are completely interconnected to form a matrix phase, the exchange between the grains has reached saturation, and T irr − T blo no longer increases with increasing ferrite content. The coercivity and T blo of the composite films is once again controlled separately by ferrite grain size, which decreases with increasing ferrite content.
In ligh of linear physical transport, the percolation threshold of the magnetic permeability in a composite can be given by fitting the experimental results to the Kirkpatrick model 24 . Under the out-of-plane and in-plane magnetic fields, the values of f c are 0.38 and 0.14, respectively, as shown in Fig. 2(b). The value of percolation threshold can also be directly obtained from the coercivity as depicted in Fig. 2(a); that is, f c ≈ 0.4 under vertical magnetic field and f c ≈ 0.2 under parallel magnetic field. These results are in good accordance with the phenomena demonstrated in Fig. 4(c,d) that the magnetic induction of the out-of-plane field starts to transfer between the magnetic particles when f NZFO ≈ 0.4, and that of the in-plane field can transfer when f NZFO < 0.3.
More interestingly, the values of f c for the electric and magnetic responses of the composite are different. Precisely speaking, their behaviors are opposite under external fields with out-of-plane and in-plane directions. In fact, such difference is caused by the same transition in the topological microstructure of the composite films as well as the nature of magnetic response. It can be speculated that both exchange and magnetic current are generated in a direction perpendicular to the external magnetic field. Hence, the percolation transition of the magnetic properties of the multiferroic composite under out-of-plane magnetic field is caused by a topological transition along in-plane direction, while the transition under a radial magnetic field is caused by a topological transformation in the normal direction of the films. In addition, since the magnetic induction can cross the non-magnetic barrier layer of thickness ~15 nm 19 , the magnetic percolation threshold in the composite films appears in a much lower ferrite content region than the threshold of electrical performance and geometrical percolation.

Discussion
The electrical and magnetic properties of the composites (Figs 2-4) confirm the mathematical model simulation of the lattice network (Fig. 1), indicating that the BTO-NZFO composite films possess different percolation thresholds in the normal and radial directions, i.e. 0.16-0.18 and 0.42-0.5, respectively. In low-dimensional BTO-NZFO composite films, the out-of-plane percolation threshold will be lower than the in-plane percolation www.nature.com/scientificreports www.nature.com/scientificreports/ threshold. This fact implies that the nanowire-like channels connecting the NZFO nanocrystals come into existence in the normal direction instead of the radial direction in the composite. If the composition of the composite films is tailored to lie between the percolation thresholds in two directions (x = 0.2-0.4), some NZFO paths that are arranged in the normal direction and not in contact with each other in the radial direction can be possibly formed in the composite films. Each NZFO path can be approximated as a one-dimensional nanowire in the ferromagnetic phase. Conclusively, it is entirely possible to form a 1-3 structure in which the ferromagnetic phase is 1D nanowire and the ferroelectric phase is a matrix. In other words, the 0.6BTO-0.4NZFO composite thin film prepared by a simple 0-3 process may have a large area of anisotropic contact between the ferroelectric and the ferromagnetic phases as shown in Fig. 5(a), thus giving rise to a magnetoelectric coupling as strong as 1-3 multiferroic composites.
As can be seen from Fig. 5(c), at the Curie point (~375 °C) of BTO where the first-order phase transition from tetragonal to cubic structure occurs, the magnetization of the 0.6BTO-0.4NZFO composite film prepared by our  www.nature.com/scientificreports www.nature.com/scientificreports/ 0-3 process indeed changes abruptly as expected, exhibiting a magnetoelectric coupling strength as large as 9.23% comparable to conventional 1-3 multiferroic composite films 5,14 . Compared with the magnitudes of the change in the magnetization of 0.8BTO-0.2NZFO and 0.4BTO-0.6NZFO composite thin films, it is evidenced that the magnetization of 0.6BTO-0.4NZFO composite is much higher. It indicates that only when the NZFO content is controlled in the region of x = 0.2-0.4, which locates between the percolation thresholds in two directions with the formation of 1-3 structure, can the magnetoelectric coupling strength of the multiferroic composites achieve a satisfactorily high level. More importantly, it has been verified that the 1-3 multiferroic composite films with strong magnetoelectric coupling can be easily prepared via a traditional 0-3 process by tailoring the anisotropy of the percolation threshold of composite films.
It should be noted that the grain sizes as well as the defects in composite films are changeable depending on the volume fraction of the components, according to our previous researches [23][24][25] . Nevertheless, the increase of NZFO grain size can be considered as a result of increasing ferrite volume. In particular, this work mainly focuses on the different percolative behaviors of the same BTO-NZFO composite thin films under electric and magnetic fields with different directions. In this scenario, the grain sizes and defects as well as their effect are supposed to be the same. Therefore, the influence of grain size and defects in lattices can be excluded in the course of interpreting the different percolative behaviors of the thin films under in-plane and out-of-plane external fields.

Methods
The BTO-NZFO composite films were prepared by magnetron sputtering as described in previous work. When the thicknesses of the BTO-NZFO composite films were controlled as 100 nm, the experimental compositions of the composite films were very closed to the nominal one (Table S1 in Supporting Information), and the grain sizes of BTO and NZFO are distributed in the range of 23.9-15.8 nm and 14.3-20.0 nm, respectively, as shown in Fig. 6 [19][20][21] .
The electrical and magnetic properties of the composite films were tested with HP4292A-LRC Agilent 4292A Precision Impedance Analyzer and MPMS-XL-5 Magnetic Measurement System manufactured by Quantum Design Instruments, USA.