Bulk Ferroelectric Heterostructures Overview

• Bulk ferroelectric heterostructures are systems where thick ferroelectric phases interact with engineered interfaces to control polarization and multifunctionality.
• They encompass architectures such as free‐standing multilayers, oxide superlattices, and internally modulated ceramics that use strain, charge, and chemical partitioning for tunable properties.
• These heterostructures enable active tuning of electronic, photovoltaic, spintronic, and thermoelectric responses, guiding next-generation high-performance device design.

Bulk ferroelectric heterostructures (BFH) are heterostructured systems in which ferroelectricity is carried by a bulk-like crystal, thick film, multilayer, superlattice, or monolithic ceramic, while electrostatic, elastic, chemical, and domain-coupled boundary conditions across internal interfaces determine the observable state. In recent work, the term spans free-standing ferroelectric–elastic multilayers that interpolate between epitaxial-film and bulk-crystal limits (Prokhorenko et al., 2010), 2D/3D stacks such as WSe$_2$/BiFeO$_3$ in which a 40 nm BiFeO$_3$ layer acts as a non-volatile internal gate (Salazar et al., 2022), oxide superlattices such as BaTiO$_3$/BaBiO$_3$ and HfO$_2$/XO$_2$ where coherency and polymorph selection stabilize new structural states (Baksi et al., 2024, Mukherjee et al., 2024), and monolithic BiFeO$_3$–BaTiO$_3$ ceramics with an approximately 8 nm periodic composition fluctuation produced by elemental partitioning (Li et al., 13 Jul 2025).

1. Scope and defining architectures

The current literature uses BFH for several closely related architectures in which the relevant ferroelectric degree of freedom is not confined to an interfacial monolayer. In one usage, a bulk-like ferroelectric oxide governs an adjacent functional layer through bound charge, as in WSe$_2$/BiFeO$_3$0, where the BiFeO$_3$1 is a continuous, laterally extended, 40 nm ferroelectric film and the top WSe$_3$2 is modulated by its polarization (Salazar et al., 2022). In another, the heterostructure itself is the periodic bulk phase, as in HfO$_3$3/XO$_3$4 superlattices treated with periodic boundary conditions in all three directions (Mukherjee et al., 2024). A third usage emphasizes monolithic ceramics in which nanoscale chemical partitioning produces internal heterointerfaces throughout a bulk body, as in BiFeO$_3$5–BaTiO$_3$6 BFH ceramics (Li et al., 13 Jul 2025) and domain-engineered BFH ceramics derived from 0.7BiFeO$_3$7–0.3BaTiO$_3$8 (Li et al., 7 Sep 2025).

BFH usage	Representative system	Defining feature
2D/3D electrostatic BFH	WSe$_3$9/BiFeO$_3$0	Bulk ferroelectric oxide as internal gate
Free-standing multilayer BFH	Ferroelectric–elastic multilayers	Tunable $_3$1, $_3$2, $_3$3
Atomically periodic BFH	HfO$_3$4/XO$_3$5 superlattices	Bulk phase with heterostructure unit cell
Monolithic ceramic BFH	BiFeO$_3$6–BaTiO$_3$7	Internal periodic composition fluctuation

Free-standing ferroelectric–elastic multilayers provide the clearest formal bridge between thin-film and bulk behavior. Their thermodynamic treatment uses the volume fractions

$_3$8

and an averaged free-energy density

$_3$9

so that $_3$0 reproduces epitaxial thin-film behavior and $_3$1 reproduces bulk-crystal behavior (Prokhorenko et al., 2010). This formulation is important because it makes BFH a continuum concept: internal elastic coupling, rather than an external substrate alone, becomes the key control variable.

The ceramic literature extends the term further by replacing layer-by-layer deposition with bottom-up heterostructuring. In BiFeO$_3$2–BaTiO$_3$3, spontaneous elemental partitioning creates coherently lattice-matched Bi-rich and Ba-rich regions with periodic modulation and strong polarization–strain–charge coupling across internal interfaces (Li et al., 13 Jul 2025). Domain-engineered BFH then uses thermal aging in a pre-engineered domain state so that the internal bias field imprints a chosen domain orientation distribution throughout the bulk ceramic (Li et al., 7 Sep 2025).

2. Thermodynamic, electrostatic, and elastic principles

The unifying mechanism in BFH is the renormalization of ferroelectric free energy by internal boundary conditions. In free-standing multilayers, mechanical equilibrium requires vanishing average in-plane stress across ferroelectric and passive layers, while strain compatibility enforces common interfacial deformation. The resulting effective Helmholtz free energy renormalizes the second- and fourth-order Landau coefficients as functions of misfit strain $_3$4, temperature, and passive-layer fraction $_3$5, producing phase diagrams that differ strongly from both bulk crystals and rigidly clamped films (Prokhorenko et al., 2010).

Electrostatic coupling is the most direct BFH control channel. In WSe$_3$6/BiFeO$_3$7, the ferroelectric surface bound charge is

$_3$8

so an UP-polarized BiFeO$_3$9 surface attracts electrons into WSe$_3$0, whereas a DOWN-polarized surface depletes them. The measured consequence is a rigid shift of the WSe$_3$1 band structure relative to $_3$2: for 3L WSe$_3$3, the valence-band maximum shifts from about $_3$4 on DOWN BFO to about $_3$5 on UP BFO, yielding $_3$6; for 2L WSe$_3$7, the shift is about $_3$8 (Salazar et al., 2022). This is a paradigmatic non-volatile ferroelectric gating effect.

A closely related principle appears in negative-capacitance heterostructures. In Pb(Zr$_3$9Ti$_2$0)O$_2$1/SrTiO$_2$2, the series relation

$_2$3

implies that if $_2$4, then $_2$5 can exceed the dielectric-layer capacitance alone. Experimentally, the composite capacitance exceeds that of the constituent SrTiO$_2$6 above about $_2$7, and Landau analysis shows that the dielectric layer stabilizes the otherwise unstable negative-curvature region of the ferroelectric free energy (Khan et al., 2011).

Structural adaptation can replace domain formation as the preferred route to electrostatic relief. In compressively strained BiFeO$_2$8 near its morphotropic phase boundary, the growth-temperature phase is an ideal tetragonal phase with no critical thickness for polarization, and capacitor-like metal|ferroelectric|metal heterostructures retain single-domain polarization during growth. The electrostatic energy is lowered not by domain formation but by a reduction in tetragonality, with density-functional calculations and STEM indicating stabilization of a metastable R-like monoclinic structure upon cooling (Nordlander et al., 2020). This shows that BFH design must track not only polarization magnitude but also the coupled structural manifold available to the ferroelectric.

3. Electronic, optical, spin, and thermal functionalities

The electronic manifestation of BFH is clearest in spectroscopic studies. In centimeter-scale WSe$_2$9/BiFeO$_2$0, ARPES directly resolves the WSe$_2$1 valence band, its K-valley spin splitting, and the polarization-controlled rigid energy shift. XPS shows that W 4f and Se 3d core levels shift rigidly by about $_2$2 between UP and DOWN BFO, with no evidence of WO$_2$3 or SeO$_2$4, so the dominant effect is electrostatic rather than chemical (Salazar et al., 2022). The combination of low strain, preserved epitaxy, and large $_2$5 establishes BFH as a route to direct band-offset engineering in 2D semiconductors.

Ferroelectric photovoltaic BFH expose a different competition: bulk symmetry breaking versus interface band bending. In Pt/BiFeO$_2$6/La$_2$7Sr$_2$8MnO$_2$9 vertical heterostructures, well above-bandgap excitation at 405 nm yields a sinusoidal photocurrent consistent with the bulk photovoltaic tensor,

$_3$0

whereas 520 nm band-edge excitation produces a response dominated by interface band bending and anisotropic absorption (Abdelsamie et al., 2021). The observed crossover resolves a longstanding ambiguity: in vertical BFH, switchable photocurrent does not by itself identify a bulk photovoltaic origin.

Spintronic BFH exploit bulk ferroelectric inversion breaking rather than only interface Rashba states. Fe/GeTe(111) is a ferromagnetic/ferroelectric-Rashba interface in which strong Fe–Te hybridization suppresses GeTe surface states, yet the bulk Rashba bands remain hardly altered beneath the interface (Slawinska et al., 2018). The Fe overlayer stabilizes both polarization orientations, and the “in” configuration is more stable than “out” by about $_3$1 eV/Å$_3$2 for 3ML and 6ML Fe. This suggests that spin-to-charge conversion at this interface is more likely to involve bulk Rashba states than surface states.

Nonlocal thermoelectric BFH extend the concept from static electrostatics to collective dipolar excitations. In a ferroelectric–metal heterostructure, the ferroelectric free energy is described by a Landau–Ginzburg–Devonshire functional,

$_3$3

and the bosonic excitations of the polarization field, “ferrons,” mediate a nonlocal Peltier and reciprocal Seebeck effect through long-range charge–dipole coupling (Tang et al., 2022). For graphene/h-BN/monolayer-ferroelectric parameters at $_3$4 nm, $_3$5, and $_3$6 K, the predicted nonlocal Peltier coefficient is $_3$7, corresponding to a temperature gradient of about $_3$8 in the ferroelectric under modest current density (Tang et al., 2022). The formalism is explicitly described as readily extendable to $_3$9 and, more broadly, to quasi-3D or bulk multilayers.

4. Structural stabilization and phase selection in oxide and fluorite BFH

Oxide BFH often operate by suppressing bulk distortions while preserving ferroelectric order. In BaTiO$_3$0/BaBiO$_3$1, oxide MBE yields trilayers and superlattices in which 4 uc BBO layers relax toward bulk-like pseudocubic parameters, whereas 2 uc BBO layers can remain coherently compressed and develop a very large tetragonal distortion with

$_3$2

(Baksi et al., 2024). Raman spectroscopy shows that the $_3$3 breathing mode at about $_3$4, associated with the charge-density-wave distortion of bulk BBO, is strongly suppressed in the superlattice. The ferroelectric response of the BTO layers is confirmed by PFM, and the system is explicitly framed as a platform for ferroelectric-field-induced electrostatic doping of ultrathin BBO (Baksi et al., 2024).

Fluorite-based BFH pursue a different target: thermodynamic stabilization of polar hafnia. First-principles studies of HfO$_3$5/XO$_3$6 superlattices with $_3$7 Si, Ge, Ti, Sn, Zr, Pb, and Ce show that composition and stacking direction can stabilize fully polar, fully antipolar, or mixed states with improved thermodynamic stability relative to polar orthorhombic HfO$_3$8 in bulk form (Mukherjee et al., 2024). The central design rule is that combining HfO$_3$9 with an oxide lacking a monoclinic ground state generally drives the superlattice away from the non-polar monoclinic phase and toward structures that minimize elastic and electrostatic penalties. In Ge/Hf, Ti/Hf, and Ce/Hf, mixed polar configurations can become favorable; in Sn/Hf, fully antipolar states are favored (Mukherjee et al., 2024).

Bulk single-crystalline HfO$_2$0:Y shows that hafnia ferroelectricity is not intrinsically confined to thin films. Laser-diode-heated floating-zone growth yields transparent single crystals in which 12% Y stabilizes a pure polar orthorhombic $_2$1 phase, while 8–11% Y yields mixed orthorhombic antipolar $_2$2, monoclinic, and cubic phases (Xu et al., 2020). Neutron diffraction, atomic imaging, and PUND measurements reveal switchable polarization, abundant $_2$3 domains, and little wake-up effect. This removes the presumed upper size limit for ferroelectricity in HfO$_2$4 and supplies a bulk fluorite ferroelectric building block for future BFH (Xu et al., 2020).

5. Domain-engineered bulk ceramics and high-temperature electromechanics

The most expansive BFH reinterpretation occurs in BiFeO$_2$5–BaTiO$_2$6 ceramics, where heterostructuring is generated internally by cation diffusion rather than by external deposition. In 0.7BiFeO$_2$7–0.3BaTiO$_2$8, elemental partitioning produces periodic composition fluctuation with a wavelength of about 8 nm, forming interpenetrating Bi-rich and Ba-rich lamellar networks within an overall perovskite structure (Li et al., 13 Jul 2025). These BFH ceramics exhibit a record $_2$9 up to $_3$00 and a room-temperature piezoelectric coefficient $_3$01, together with less than 3% loss in $_3$02 after 1000 h at $_3$03 in poled BFH (Li et al., 13 Jul 2025). Distinct electromechanical boundary conditions during thermal aging generate distinct morphologies of aliovalent A-site segregated regions and distinct ferroelectric orders: a relaxor/antiferroelectric-like BFH-AS state with random local fields, and a biased hard-ferroelectric BFH-AQ state with an internal bias field of about $_3$04 (Li et al., 13 Jul 2025).

Domain-engineered BFH push this logic from passive stabilization to active texture programming. In 0.7BiFeO$_3$05–0.3BaTiO$_3$06, thermal aging in a mechanically or electrically pre-engineered state produces internal bias fields up to $_3$07, yielding macroscopic imprint of the chosen domain orientation distribution (Li et al., 7 Sep 2025). In mechanically domain-engineered BFH, the remanent compressive strain state is locked in and then reversibly recovered under field, giving axial electrostrain up to $_3$08 at $_3$09. In electrically domain-engineered BFH, transverse-field-driven non-$_3$10 reorientation yields shear strain $_3$11 unipolar and $_3$12 peak-to-peak at $_3$13, increasing to about $_3$14 at $_3$15, with effective $_3$16 of about $_3$17 at room temperature and about $_3$18 at $_3$19 (Li et al., 7 Sep 2025). The paper characterizes the shear result as a record high peak-to-peak shear strain up to 0.9% at intermediate field levels.

A different bulk route to multifunctionality is the self-grown ferroelectric–ferromagnetic composite in the BiFeO$_3$20–BaTiO$_3$21 alloy system. Detailed magnetic separation and XRD show that the observed ferromagnetism originates from spontaneous precipitation of a minor BaFe$_3$22O$_3$23 phase at about $_3$24, not from homogeneous intrinsic ferromagnetism of the perovskite matrix (Kumar et al., 2017). Despite this small fraction, the composite exhibits enhanced ferroelectric switching under magnetic field and a dc magnetoelectric coupling of about $_3$25, comparable to layered laminates and bilayer thin-film heterostructures (Kumar et al., 2017). This is a BFH in the literal bulk-composite sense: a dense monolith in which a continuous ferroelectric matrix is internally coupled to sparse ferrimagnetic inclusions.

6. Experimental probes, interpretive debates, and open directions

BFH research is methodologically heterogeneous because the relevant observables range from reciprocal-space order parameters to local polarization vectors and spectroscopically resolved band dispersions. PFM, AFM, XRD, reciprocal-space mapping, synchrotron CTR analysis, Raman spectroscopy, ARPES, XPS, ISHG, HAADF-STEM, STEM-EDS, STEM-EELS, PUND, in-situ synchrotron XRD, and digital image correlation each address a different part of the coupled electrostatic–structural–domain problem (Salazar et al., 2022, Baksi et al., 2024, Abdelsamie et al., 2021, Li et al., 7 Sep 2025). The field increasingly relies on joint structural and functional datasets because single measurements are often insufficient to distinguish bulk, interface, and defect-mediated mechanisms.

Several recurrent controversies are now better constrained. First, switchable response is not a sufficient identifier of a bulk ferroelectric mechanism: the BFO photovoltaic work shows that bulk photovoltaic and interface-Schottky responses can coexist in the same vertical heterostructure and exchange dominance with photon energy (Abdelsamie et al., 2021). Second, apparent interfacial capacitance enhancement requires discrimination from leakage and Maxwell–Wagner artifacts; the PZT/STO negative-capacitance study explicitly rules out a Maxwell–Wagner origin by frequency and resistance analysis up to 1 MHz (Khan et al., 2011). Third, magnetic functionality in nominally single-phase BiFeO$_3$26–BaTiO$_3$27 alloys can arise from an undetected minor hexaferrite precipitate, so phase-purity claims require magnetic separation or equally sensitive methods (Kumar et al., 2017). A broader implication is that BFH should not be reduced either to interfacial electrostatics alone or to nominal chemical composition alone; the operative state is usually a coupled product of phase purity, coherency, internal bias, and accessible domain texture.

The immediate research trajectory is toward operando control and architecture generalization. Static UP/DOWN spectroscopy in WSe$_3$28/BiFeO$_3$29 motivates operando or time-resolved ARPES during ferroelectric switching (Salazar et al., 2022). BTO/BBO indicates that controlled intermixing, coherent strain, and BBO thickness remain decisive for realizing the predicted interfacial electronic phases (Baksi et al., 2024). HfO$_3$30-based studies suggest a route from kinetic stabilization in bulk single crystals to thermodynamic stabilization in superlattices (Xu et al., 2020, Mukherjee et al., 2024). The ceramic BFH program points toward scalable high-temperature actuators, sensors, energy harvesters, multiple-state memory devices, and domain-wall switches, while also implying that analogous concepts may be extended across the wider ferroic family (Li et al., 13 Jul 2025, Li et al., 7 Sep 2025).

A plausible synthesis of the literature is that BFH are best understood not as one specific geometry but as a design regime. In that regime, a ferroelectric order parameter with bulk character is embedded in a heterostructured environment that is sufficiently extended to renormalize its free-energy landscape, yet sufficiently coherent that electrostatic, elastic, and domain-mediated coupling remain deterministic. Under that definition, free-standing multilayers, oxide superlattices, 2D/3D electronic stacks, and internally modulated ceramics are not competing meanings of BFH but experimentally distinct realizations of the same materials principle (Prokhorenko et al., 2010, Li et al., 13 Jul 2025).