Unlocking Strong Zero-Bias Magnetoelectric Effects in Bilayers of PZT and Self-Biased Hexagonal Ferrite–Spinel Ferrite Composites

1. Introduction

Multiferroic magneto-electric (ME) composites, which combine distinct ferromagnetic and ferroelectric phases, offer a promising route toward controlling magnetic states via electric fields, and vice versa, through strain-mediated interfacial coupling [1,2,3,4,5]. In these systems, the overall ME conversion efficiency depends heavily on the product of the piezoelectric coefficient d and the piezomagnetic coefficient q = dl/dH (where l is the magnetostriction) [6]. A long-standing practical bottleneck for the use of ME composites in useful technologies, however, is the reliance on bulky, power-consuming external DC bias fields H to maximize q [1,2]. Eliminating this external bias requirement is critical if these materials are to be integrated into compact magnetic sensors, energy harvesters, or high-resolution imaging arrays [2,7,8]. Consequently, engineering “self-biased” ME structures that can deliver strong zero-bias ME coupling strengths has become a central focus in the field [9].

Past attempts to engineer self-bias have typically relied on complex structural or processing interventions. These include compositionally grading the magnetic phase [10,11,12,13,14], using thermal annealing to induce an exchange bias [15], exploiting hysteretic remanence [9], or built-in stress engineering [9]. While these methods work, they often introduce processing complications, material fatigue, or strict geometry constraints. A more straightforward solution is to introduce a magnetically hard phase with a high magneto-crystalline anisotropy field H_A to act as a built-in magnetic bias [16]. M-type barium hexaferrite BaFe₁₂O₁₉ (BaM) is a prime candidate here due to its strong uniaxial anisotropy field H_A ~ 18 kOe and high remanence [17,18,19,20,21]. Yet, a fundamental material trade-off limits its utility: while BaM provides an excellent internal bias, its saturation magnetostriction l ~ 1–2 ppm is very low, making pure BaM composites highly inefficient for strain transfer [22].

On the other end of the spectrum, soft spinel ferrites like nickel-zinc ferrite Ni_0.8Zn_0.2O₄ (NZFO), exhibit large saturation magnetostriction l ~ 30 ppm and high permeability, but they lack any self-biasing capability because of their soft magnetic nature [23,24,25]. Interestingly, while such hard/soft ferrite composites have been widely investigated, the literature has focused almost exclusively on exchange-spring interactions in nanoscale systems aimed at optimizing permanent magnet performance [26,27,28,29,30]. The cooperative use of these two phases to deliberately bypass the anisotropy-magnetostriction trade-off in bulk ME structures has been largely overlooked. Specifically, the physics of decoupling these phases at a microstructural scale—retaining high bulk-like anisotropy while ensuring efficient strain transfer across a ferroelectric interface—remains unaddressed.

In this paper, we address this trade-off by developing a bulk ferromagnetic composite made of micron-sized BaM and NZFO grains, configured as a bilayer with lead zirconate titanate (PZT). By deliberately avoiding the nanometer-scale exchange-coupling regime, we preserve the large macroscopic anisotropy fields needed for self-biasing, while simultaneously capturing the strong magnetostrictive strain of the soft spinel phase [31]. This distinct phase separation allows the internal flux of the hard phase to bias the soft phase without damping its magneto-elastic compliance.

Polycrystalline, dual-phase ferromagnetic composites with (100 − x) wt.% BaM and x wt.% of NZFO (BNZx) (x = 10–90) were used in our study. Advances in thin film deposition of BaM enabled precise control over the microstructure and crystallographic orientation that led to films with H_A of 17.5 to 18.5 kOe [19]. Although achieving consistently high H_A values in polycrystalline BaM remained a challenge due to the influence of grain boundaries [20], there are past reports on the synthesis of polycrystalline BaM with H_A as high as 16.6 kOe [21]. The phase purity of the composites we prepared by ceramic techniques was confirmed by X-ray diffraction. The morphological features of the samples were examined by scanning electron microscopy (SEM). Magnetic characterization of the composites was done by magnetization M measurements, magnetostriction λ, and ferromagnetic resonance (FMR). The anisotropy field H_A was determined from M vs. H as well as FMR data. Platelets of BNZx were bonded to lead zirconate titanate (PZT) to form bilayers and the ME coupling strengths were measured by ME voltage coefficient (MEVC) as a function of frequency and external bias H. We report an exceptional zero-bias MEVC at low frequencies and at electro-mechanical resonance modes. Details on the synthesis of the ferromagnetic composites and characterization and results on ME measurements in bilayers with PZT are provided in the sections that follow.

2. Experiment

Conventional ceramic processing techniques were utilized to synthesize the ferrite composites. Initially, micrometer-sized polycrystalline powders of Ni_0.8Zn_0.2Fe₂O₄ and BaFe₁₂O₁₉ were prepared separately [17,32]. For NZFO synthesis, 99.9% pure NiO, ZnO, and Fe₂O₃ were weighed in stoichiometric proportions, ball-milled for 8 h, dried, and pre-sintered at 900 °C for 6 h. This was followed by a second ball-milling and sintering at 1200 °C. BaM powder was synthesized with 99.9% pure BaCO₃ and Fe₂O₃ using a similar procedure, including sintering the powder at 1200 C. Composites of (100 – x) wt.% BaM and x wt.% NZFO, denoted by BNZx, where x = 10–90, in steps of 10, were prepared by mixing the individual ferrite powders in a ball-mill for several hours, dried, and a 2 wt.% polyvinyl alcohol (PVA) solution was used as a binder. The mixtures were compacted in a die into disk-shaped pellets (diameter ≈ 18 mm, thickness ≈ 1.5 mm) under an uniaxial pressure of 250 MPa. The pellets were then sintered at 1250 °C for 6 h to ensure densification. Following sintering, the pellets were sliced into rectangular platelets and used for structural and magnetic characterization.

Structural characterization of the composites was performed using an X-ray diffractometer (Miniflex, Rigaku, Japan) and an SEM (JSM-6510/GS, JEOL, Tokyo, Japan). Magnetostriction λ as a function of in-plane static field H was measured on thin rectangular samples using a strain gage and a strain indicator (Model P3, Micro-Measurements, Wendell, NC, USA). A two-part epoxy (M-bond) was used to attach the strain gage. The magnetic field H was applied parallel to the length (direction-1) of the sample and the strain gage. The magnetostriction obtained in this configuration is designated as λ₁₁. Magnetization as a function of H up to 2 T was measured on thin rectangular samples with a Physical Property Measurement System (PPMS, Quantum Design Inc., San Diego, CA). Data were obtained for H parallel and perpendicular to the sample plane. Ferromagnetic resonance (FMR) measurements were carried out on thin rectangular platelets of the composites placed in an S-shaped coplanar waveguide. A 67 GHz vector network analyzer (E8361A, Agilent, Santa Clara, CA, USA) was employed to record the S₂₁ vs. frequency f profiles for a series of H applied parallel to the sample plane.

Measurements of ME coupling strengths were performed on bilayers of BNZx and vendor-supplied PZT (PZT-850, American Piezo Ceramics, PA). Rectangular composite platelets of lateral dimensions 18 mm × 3mm and 0.4mm to 1 mm in thickness were epoxy bonded to 0.3 mm-thick PZT plates of similar lateral dimensions. The ME voltage coefficient (MEVC) measurements were made with the samples placed in an aluminum sample holder that was an effective RF shield to reduce noise pickup. The measurement system consisted of an electromagnet to apply a static field H and a pair of Helmholtz coils (placed inside the sample holder) to apply an AC magnetic field h. The AC voltage induced across the PZT layer was measured with a lock-in amplifier. The MEVC was measured as a function of H, frequency f of the AC field, and for H and h parallel to each other and either parallel or perpendicular to the sample plane.

3. Results

3.1. Structural Characterization

Figure 1 and Supplementary Figure S1 show representative X-ray diffraction (XRD) patterns for the BNZx samples. All major diffraction peaks have been indexed and are consistent with the characteristic reflections of the hexagonal BaM and cubic spinel NZFO phases [17,32,33]. For BNZ10, the diffraction pattern is dominated by the BaM reflections due to its high BaM content, while only weak NZFO peaks are present. As the NZFO concentration increases from BNZ40 through BNZ50 and BNZ90, the intensity of the characteristic peaks for NZFO progressively increases, whereas those of BaM gradually decrease in intensity. Although diffraction peaks due to any impurity phases did not appear for samples with x-values up to 30, a weak peak identified with antiferromagnetic BaFe₂O₄ (BFO) as seen for x ≥40 and did not contribute to the magnetization of BNZx [20]. In our system, a small amount of the impurity phase, BFO, was identified in the XRD data only for high weight percent of NZFO in the ferromagnetic composites. Despite the large ionic radius of Ba and Sr, previous reports in the literature reported on the incorporation of both elements in nickel ferrite [34]. Our efforts to completely eliminate the impurity were unsuccessful.

The estimated lattice constants for pure BaM and NZFO and for the hexagonal ferrite in the BNZx composites for selected x-values are listed in Table 1. The lattice constant a = 8.3844 A° for NZFO did not show any variation with x and is not listed in Table 1, whereas BaM showed an increase in both a and c with x-values, with the largest increase in a by 0.07% and in c by 0.13% for x = 10. The corresponding increase in the unit cell volume is 0.26%. With further increase in the amount of NZFO, however, both a and c show a gradual decrease to values slightly smaller than for pure BaM. The possible cause of the variations in the lattice constants for BaM in the composites are discussed in Section 4.

Representative SEM images for BNZx for x = 10, 20, 30, and 90 are shown in Figure 2. Images for other x-values are shown in the Supplementary Figure S2. Hexagonal crystallites of BaM are seen in the images for composites with low x-values. A gradual increase in hexagonal features with x is also evident for samples with 10 < x < 40. However, for 50 < x < 90 as the NZFO content increases (i.e.,), the hexagonal crystallites become small and less prominent, and larger grains likely belonging to NZFO appear in the SEM images. This behavior is consistent with our previous observations in composites of Nickel Ferrite and Barium Ferrite (BaM) [31].

3.2. Magnetic Characterization

Magnetization (M) measurements as a function of H were done on BNZx composites for in-plane and out-of-plane fields. Representative M vs. H data for x = 60, 70, and 80 are shown in Figure 3 (Figure S3 in the Supplement shows M vs. H data for samples with other x-values). For in-plane H, data in Figure 3 show a rapid increase in M to a saturation value M_s. For out-of-plane H, M attains saturation at a much higher value and is indicative of an in-plane crystalline anisotropy for all of the BNZx composites. Figure 4a shows M_s as a function of x, and one infers from Figure 4a that M_s increases with an increase in x for x ≥ 0.5. For comparison, we estimated the expected M_s values for BNZx from M_s for pure BaM and NZFO (from data in the Figure S3). For M_s = 64 emu/g for NZFO and 82 emu/g for BaM, the expected M_s vs. x are plotted in Figure 4b along with the measured values. The measured Ms are lower than the expected values and could be due to migration of Ni and nonmagnetic Zn from NZFO to BaM. Such diffusion of Zn and Ni could also account for the variation in the lattice constants of BaM in the composites (Table 1).

Figure 4b shows the anisotropy field H_A estimated from the intersection field (H_i) corresponding to the overlap of the in-plane and out-of-plane M vs. H hysteresis loops. Analysis of the M vs. H data indicates that for lower x compositions, the positive values of H_A suggest a predominant in-plane magnetic anisotropy. However, for higher x (BNZ80 and BNZ90), the observed negative H_A indicates a transition toward out-of-plane anisotropy. A general decrease in H_A with increasing x is evident from the data in Figure 4b. Since single-crystal M-type hexagonal ferrites have uniaxial anisotropy field along the c-axis, the positive H_A values imply that a significant fraction of crystallites of BaM in the polycrystalline composites with NZFO have in-plane (a,b) orientation, whereas the negative H_A at higher x suggests an increasing contribution of c-axis oriented grains (out-of-plane alignment) in the polycrystalline BNZ composites.

3.3. Ferromagnetic Resonance

We also carried out FMR measurements on the rectangular platelets of the BNZx composites. The samples were placed in an S-shaped coplanar waveguide and excited with microwave power from a VNA. A static magnetic field H was applied parallel to the sample plane. Profiles of the scattering matrix S₂₁ as a function of frequency f were recorded for a series of H-values. FMR originating from NZFO was observed in the frequency range 5–30 GHz for composites for x ≥ 30. Resonances associated with magneto-dielectric modes were observed in the 30–60 GHz frequency range and are not considered here for further analysis [35]. Representative profiles of S₂₁ vs. f for BNZ60 and BNZ90 for in-plane H are plotted in Figure 5, and the dependence of the resonance frequency (f_r) on H is shown in Figure 5c. FMR profiles for other x values are shown in the Supplementary Figure S4.

One can fit the f_r vs. H data to the Kittle equation given by

$$f_r=\gamma {\left(\left(H+\left(N_z-N_x\right)4\pi M_{eff}\right)\left(H+\left(N_y-N_x\right)4\pi M_{eff}\right)\right)}^{{\displaystyle \raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}}$$

(1)

where g is the gyromagnetic ratio, H is the magnetic field along the x-direction, and N_x, N_y, and N_z are the demagnetization factors along the platelet’s length, width, and thickness, respectively. These demagnetization factor values for each composite are listed in Table S1 in the Supplement. The effective magnetization 4πM_eff and the gyromagnetic ratio (g) were directly obtained from the fitting results. The anisotropy field H_A was then calculated from 4πM_eff and 4πM values at 5 kOe. The calculated H_A values from FMR and M vs. H are plotted in Figure 4b (also listed in Table S2 in the Supplement). FMR values of H_A are smaller than the measured values from M vs. H data. The discrepancy in the two HA values could be due error in the demagnetization factors we used assuming a near-ellipsoidal shape for our rectangular samples [36]. However, there is overall agreement in the trend in the variation of H_A with x obtained from M vs. H and FMR. The most important inference from H_A values is the high crystalline in-plane anisotropy field in the composites that is a key ingredient for realizing ME coupling under zero external magnetic bias in a composite with a ferroelectric. The values of g range from 2.5 MHz/Oe for BNZ30 to 2.7 MHz/Oe for BNZ90. These values are smaller than 3 MHz/Oe reported for NZFO [23,37].

Magnetostriction of the ferromagnetic phase is a key parameter that determines the strength of ME coupling in the BNZx-PZT bilayers. We measured the variation in l with H in the platelets of BNZx using a strain gage bonded to the sample surface. A strain indicator was used for the measurements of magnetostriction λ₁₁ parallel to the sample length and to the applied magnetic field H. Figure 6 shows the λ₁₁ vs. H data for BNZx samples and pure NZFO. All the samples show negative values as expected [38]. With increasing NZFO content, λ₁₁ gradually increases. For BNZ10 highest value of λ₁₁ is −2 ppm, comparable to −1ppm for pure BaM [22,31]. For BNZ90, λ₁₁ increases with H to −24ppm, a value comparable to the pure NZFO.

3.4. Magneto-Electric Effects in the BNZx-PZT Bilayers

The BNZx platelets were bonded to vendor-supplied PZT with d₃₃ = −400 pm/V and their ME response was evaluated by measurements of the MEVC that depends on the piezomagnetic coupling coefficient q = dλ/dH, and the piezoelectric coefficient d of PZT [1,39]. The H-dependence of MEVC was measured for two different field orientations. (i) α₃₁ for H and ac magnetic field h parallel (direction1) to the sample length (direction 1), and V measured across PZT thickness (direction 3). (ii) α₃₃ for both H and h perpendicular to the bilayer plane, with V measured across the PZT thickness. The MEVC is given by MEVC =V/(h × t), where t is the thickness of the PZT.

Figure 7 and Supplementary Figure S5 show the H dependence of α₃₁ for the ac field at 100 Hz. The ME response rises with increasing magnetic field, reaches a maximum, and then gradually decreases and its magnitude essentially tracks variation in q with H. When λ reaches near saturation MEVC drops to very small values [1,39]. For low NZFO content (x ≤ 40), the α₃₁ vs. H shows a broad peak, and for (x ≥ 50) the peak value of MEVC increases as the Ni-Zn ferrite content increases. The maximum value of α₃₁ increases from 3 mV/cm Oe for BNZ10-PZT to 161mV/cm Oe for BNZ90-PZT bilayer. This trend agrees well with the magnetostriction data, where the magnitude of λ₁₁ increases with increasing NZFO content [39].

Figure 8 and Figure S6 in the Supplement show α₃₁ vs. H around H = 0 on an expanded scale. A zero-bias ME response is observed in all bilayers, as expected (Supplementary Table S2), and α₃₁ increases with NZFO content, from 1 mV/cm Oe for BNZ10-PZT to 114 mV/cm Oe for BNZ90-PZT. Hysteresis is observed in the H dependence of α₃₁ for all bilayers. Upon reversal of the H direction, the magnitude of α₃₁ initially decreases and then gradually increases with H to a maximum, after which it falls off rapidly at higher H [1,39,40]. A weak secondary maximum in α₃₁ near H = 0 is observed in Figure 8. This feature could be due to the dual-phase nature of the ferrite composite. The hard BaM and soft-magnetic NZFO phases exhibit different coercive fields, anisotropies, and magnetostrictive responses, leading to sequential magnetization reversal during field sweeping.

Figure 9 (and Supplementary Figure S7) shows the H dependence of α₃₃ for the ac field at 100 Hz and α₃₃ are much lower compared α₃₁ due to a large demagnetization factor for out-of-plane DC and AC magnetic fields [41]. The magnitude of α₃₃ increases with increasing NZFO content in the composites, from ~2 mV/cm Oe for BNZ10-PZT to 12 mV/cm Oe for BNZ60-PZT bilayer. The small and nearly constant magnetostriction of BNZ10 makes ME response weak.

The zero-bias region of α₃₃ vs. H is shown in Figure 10. A self-biased signal is still present, but it is weaker than in the longitudinal mode. The Supplementary Figure S8 shows the same trend for the remaining compositions. Some intermediate compositions show sharper transitions near zero field, which may be associated with stronger anisotropy-driven magnetic switching behavior. However, these do not lead to a large ME voltage, indicating that anisotropy alone is not enough. A large anisotropy field helps create internal bias, but it does not contribute to ME response. Effective ME coupling requires a balance between anisotropy-induced bias and magnetostrictive strain transfer [17,39].

The frequency dependence of MEVC was measured for BNZx-PZT bilayers. Results of the measurement of MEVC vs. f for BNZx-PZT (x = 20, 60, 90) are shown in Figure 11 and in the Supplementary Figure S9 for other x values. Prior to these measurements, we obtained the resonance frequencies for the bending and longitudinal acoustic modes by measuring the frequency dependence of the impedance. The data were obtained for both at zero-bias and for H-values corresponding to the maximum in low-frequency MEVC.

The frequency-dependent ME response exhibits two distinct resonance peaks. The first (lower-frequency) peak corresponds to the bending (flexural) resonance mode, whereas the second (higher-frequency) peak corresponds to the longitudinal acoustic resonance mode. The resonance modes were identified from impedance measurements and are consistent with the resonance behavior reported for magnetoelectric composites [1,4]. For BNZx—PZT bilayers with x = 20, 40, 50, 80, 90, the first resonance (bending mode) produces the larger MEVC, which can be attributed to enhanced flexural deformation leading to amplified strain in the piezoelectric layer, particularly away from the neutral axis, as discussed in Ref. [4]. In contrast, for BNZx (x = 10, 30, 60, 70) samples, the second resonance corresponding to the longitudinal acoustic mode exhibits the higher MEVC, indicating that axial strain transfer is more efficient in these compositions. This variation in the dominant resonance mode indicates that the interplay between bending and longitudinal vibrations is governed by changes in magnetostriction, elastic stiffness, and magnetic anisotropy with BNZx content, which determine the efficiency of strain transfer across the bilayer.

Among the BNZx-PZT bilayers, the highest MEVC of 6V/cm Oe is measured at the bending resonance frequency of 10.6 kHz for BNZ90-PZT at H= 41 Oe, and the corresponding zero-bias MEVC = 3.6V/cm.Oe. For lower x-values, x = 10–50, the MEVC values exhibit relatively low values at resonance: less than 1 V/cm Oe, suggesting weaker magnetoelectric coupling. This could be attributed to lower piezoelectric activity at these concentrations. The MEVC values increase significantly with increasing Ni-Zn composition for x = 60–90 due to strong magnetostrictive response in the ferromagnetic phase. Overall, the results in Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 and in Supplementary Figures S5–S9 show a clear enhancement in the strengths of ME coupling as the amount of NZFO is increased.

4. Discussion

The experimental results demonstrate that phase-pure composites combining spinel and M-type hexagonal ferrites can be synthesized without detrimental secondary phases, ultimately yielding a strong zero-bias ME response when coupled with PZT. The signature XRD peaks for both constituent phases—NZFO and BaM—confirm that the structural integrity of each phase is preserved during sintering. Retaining the phase purity of the highly anisotropic BaM ensures the development of a strong internal magnetocrystalline anisotropy field, while the unreacted NZFO phase provides the necessary magneto-elastic compliance. Although minor antiferromagnetic impurities were detected in the BNZ40, BNZ70, and BNZ80 samples, they did not discernibly alter the macroscopic magnetic characteristics of the bulk composites. This structural distribution is supported by the SEM surface morphology, which reveals distinct hexagonal BaM grains that naturally become more prominent as the NZFO weight fraction decreases.

As expected, the NZFO-rich composites (x > 60) exhibit high magnetostriction levels approaching those of pure zinc-substituted nickel ferrite, whereas the BaM-rich sample (BNZ10) displays poor magnetostrictive strain. To evaluate how these mixed phases influence the overall magnetic behavior, the magnetic anisotropy was evaluated using both FMR data and macroscopic M vs. H loops. The internal anisotropy field H_A derived from FMR yielded consistently higher values for the NZFO-rich samples (x = 60, 70, 80, 90) compared to the values extracted via static M vs. H (Table S1: Supplement). These fields operate in-plane, suggesting a preferential alignment of the BaM grains within the (a,b) crystallographic planes. Concurrently, H_A from the intersection of the in-plane and out-of-plane magnetization curves (H_A= H_i − 4pM_s) steadily decreases with increasing NZFO content. For the BaM-rich compositions, the positive H_A values confirm an in-plane easy axis alignment. However, as the soft NZFO phase begins to dominate, this internal field drops, ultimately flipping to small, negative values for BNZ80 and BNZ90. This negative sign signifies a clear transition toward an out-of-plane easy axis orientation. While highly textured BaM films typically display a strong out-of-plane anisotropy depending on grain orientation [42], the bulk BNZ composites studied here show a predominant in-plane anisotropy for x < 70 (peaking near BNZ30) before shifting to an out-of-plane character at high spinel concentrations (x = 80, 90). The discrepancy between H_A values estimated from FMR and M vs. H might be due to sample response under high-frequency and static magnetic measurements.

While ferrite-PZT bilayers are widely documented, this configuration differs from conventional designs by utilizing a compositionally engineered magnetic layer that systematically balances a hard hexaferrite with a soft, highly magnetostrictive spinel ferrite. Rather than substituting one ferrite phase entirely, varying the BaM/NZFO weight ratio serves as a mechanism to balance internal anisotropy against macro-magnetostriction, optimizing both the internal magnetic bias and the resulting strain-mediated ME coupling. This approach bypasses the need for an external DC bias field. Furthermore, the selection of a zinc-substituted spinel ferrite provides a noticeable performance advantage over previously studied BaM-NiFe₂O₄ configurations [31].

Using a self-biased magnetic phase containing both hexagonal and spinel structures clearly provide a viable pathway to strong zero-bias ME behavior. Previous studies on Zn-substituted spinel ferrite-PZT multilayers reported low-frequency ME voltage coefficients of 280 mV/cm Oe for cobalt zinc ferrite-PZT and 600 mV/cm Oe for NZFO-PZT under an applied magnetic field [44]. Although these systems exhibit larger MEVC than BNZx-PZT, their MEVC under zero-bias was essentially zero. In contrast, incorporating BaM with high magnetocrystalline anisotropy into the NZFO-PZT bilayers introduces an internal magnetic bias that enables substantial ME coupling without an external dc magnetic field. In our previous study on BaM/NiFe₂O₄—PZT bilayers, the maximum values of a₃₁ under zero-bias were 22 mV/cm Oe and 0.9 V/cm Oe for low-frequency and at EMR, respectively [31]. In comparison, in the present system, the corresponding zero-bias values are a₃₁ = 113 mV/cm Oe and 3.6 V/cm Oe. Thus, the enhancement in l due to Zn-substitution in nickel ferrite resulted in a substantial strengthening of ME coupling across both low-frequency and resonance regimes, outperforming previous values for hexaferrite-based ME composites [16,31]. This confirms that Zn substitution acts as a primary driver for increasing magnetostrictive strain while allowing the composite to retain enough collective internal anisotropy for self-biased operation, a result of the enhanced magneto-elastic compliance of the modified spinel phase [5,44]. Ultimately, these results underscore the role of compositional tuning in the BNZ layer for tailoring the magnetic and ME properties of the bilayer system. Interlocking the high l of the Ni-Zn ferrite with the permanent internal flux of the BaM phase yields a composite system capable of high sensitivity under zero-bias conditions. This combination makes these BNZx-PZT bilayers highly practical options for compact magnetic field sensors and stress-mediated energy-harvesting systems [45].

Although magnetoelectric bilayers based on ferrite-PZT have been extensively investigated, the present work differs from previous studies by employing a compositionally engineered magnetic layer consisting of both a hard M-type hexaferrite (BaM) and a highly magnetostrictive Ni-Zn spinel ferrite (NZFO). Rather than simply replacing one ferrite with another, the BaM/NZFO ratio was systematically varied to tailor the balance between magnetocrystalline anisotropy and magnetostriction, thereby optimizing the internal magnetic bias and magnetoelectric coupling. This compositional engineering enables strong self-biased magnetoelectric behavior without the use of externally applied DC magnetic bias and demonstrates that Zn substitution in the spinel ferrite is an effective route for significantly improving the magnetoelectric performance compared with previously reported BaM-NiFe₂O₄ composites.

5. Conclusions

In summary, dual-phase ferromagnetic composites comprising M-type barium hexaferrite, BaM, and nickel-zinc spinel ferrite (NZFO) were successfully synthesized via conventional solid-state reaction and integrated with PZT to form bilayer structures. The primary goal was to achieve strong ME coupling without the need for an external DC bias field. Ferromagnetic resonance and a law-of-approach analysis of the in-plane and out-of-plane M vs. H hysteresis loops confirmed that the composites maintain strong macroscopic magneto-crystalline anisotropy, as evidenced by the high anisotropy field resolved in the BNZ30 sample.

ME voltage coefficient measurements revealed a strong compositional dependence. The BNZ90-PZT bilayer (rich in NZFO) exhibited the highest overall performance, yielding a low-frequency MEVC of 161.5 mV/cm Oe and a peak value of 6 V/cm Oe at the electromechanical resonance frequency of 10.6 kHz under a residual bias field of ~40 Oe. Crucially, this same composition demonstrated a strong zero-bias values of 114 mV/cm Oe at 100 HZ and 3.9 V/cm Oe at resonance. In contrast, the lowest performance was observed for the BNZ10-PZT bilayer, which registered just 3.3 mV/cm Oe at 100 Hz. These findings show that compositional tuning of the magnetic phase is a viable path for optimizing self-biased ME composites. Balancing the internal field of hard ferrites with the magneto-elastic compliance of soft ferrites allows for strong zero-bias cross-coupling without external biasing elements. Consequently, these BNZx-PZT bilayers, particularly those with higher spinel content, hold distinct potential for compact magnetic sensors, tunable microwave devices, and energy-harvesting technologies.