Remote Epitaxy of BaTiO₃/Graphene

• The paper demonstrates that remote epitaxy of BaTiO₃/graphene harnesses substrate potentials transmitted through graphene to achieve epitaxial registry while enabling membrane exfoliation.
• Methodological advances include a two-step PLD process using inert and oxidizing atmospheres to mitigate both ballistic and oxidative damage to graphene.
• The study shows that optimizing graphene microstructure and transfer protocols yields high-quality, freestanding BaTiO₃ membranes with robust ferroelectric properties for device integration.

Remote epitaxy via two-dimensional (2D) interlayers such as graphene enables the synthesis of single-crystal oxide membranes that inherit the crystallographic registry of their substrates, yet can be mechanically exfoliated and transferred to diverse platforms. For the BaTiO$_3$/graphene system, fundamental advances center on the interplay between oxide nucleation, graphene integrity under pulsed laser deposition (PLD), and scalable transfer processes. This synthesis provides an avenue to combine functional perovskite oxides with silicon technologies or flexible electronics, contingent on controlling atomic-scale phenomena during thin-film growth (Haque et al., 2024, Haque et al., 14 Nov 2025).

1. Principles of Remote Epitaxy Through Graphene

Remote epitaxy exploits the partial transmission of long-range electrostatic or ionic substrate potentials through an atomically thin interlayer—such as monolayer graphene (thickness $t \approx 3.4$ Å). The substrate’s potential, typically characterized by $\Phi(z)\propto \Phi_0 e^{-kz}$ where $z$ is vertical distance and $k \approx 2\pi/a_\mathrm{sub}$ is determined by substrate lattice constant $a_\mathrm{sub}$, tailors adatom adsorption and diffusion atop the 2D sheet. In BaTiO$_3$/graphene/SrTiO$_3$ heterostructures, this potential tail exceeds typical thermal energies during film growth ($T \sim 700\,^\circ$C), establishing epitaxial registry through the graphene, even while van der Waals binding is weak. The net adatom–substrate interaction energy is modeled as $$E_\mathrm{int} = \int\rho_\mathrm{ad}(r)\Phi(r)\,d^3r$$, showing that the templating effect persists despite the nominal chemical inertness of graphene (Haque et al., 2024).

Misfit strain, parametrized by $t \approx 3.4$0 (with $t \approx 3.4$1 for BaTiO$t \approx 3.4$2/SrTiO$t \approx 3.4$3), is only partially relieved in the presence of graphene due to reduced film–substrate binding energy. This translates into unique strain states both before and after exfoliation, directly impacting functional properties.

2. Pulsed Laser Deposition Schemes and Damage Mitigation

PLD, the primary method for BaTiO$t \approx 3.4$4 remote epitaxy on graphene, involves ablation with a 248 nm KrF excimer laser. Key parameters include laser fluence ($t \approx 3.4$5 to $t \approx 3.4$6), repetition rate ($t \approx 3.4$7–$t \approx 3.4$8), and chamber background (O$t \approx 3.4$9 or Ar). The kinetic energy distribution of ionic plume species is bimodal: “fast” Ba-ion components ($\Phi(z)\propto \Phi_0 e^{-kz}$0) can directly induce ballistic defects, while “slow” scattered components ($\Phi(z)\propto \Phi_0 e^{-kz}$1 at $\Phi(z)\propto \Phi_0 e^{-kz}$2) are chemically reactive but less damaging (Haque et al., 14 Nov 2025).

Preserving graphene requires minimizing both thermal oxidation and ballistic impact. Autonomy-driven parameter mapping demonstrates that monolayer graphene survives at $\Phi(z)\propto \Phi_0 e^{-kz}$3C and $\Phi(z)\propto \Phi_0 e^{-kz}$4, but these conditions yield poorly crystallized BaTiO$\Phi(z)\propto \Phi_0 e^{-kz}$5. Consequently, optimized workflows employ two-step growth:

• Nucleation in inert Ar ($\Phi(z)\propto \Phi_0 e^{-kz}$6, $\Phi(z)\propto \Phi_0 e^{-kz}$7C, $\Phi(z)\propto \Phi_0 e^{-kz}$8), generating $\Phi(z)\propto \Phi_0 e^{-kz}$95–10 nm of c-axis BTO nuclei.
• Stoichiometric BTO growth in O$z$0 ($z$1), completing the remaining 40–50 nm for ferroelectricity, after much of the graphene is shielded.

Damage is further reduced by laser aperture control, restricting the ablation spot during nucleation, which decreases defect metrics quantified by Raman D:G intensity ratio ($z$2) and D-peak full width at half maximum ($z$3). Typical metrics with damage-mitigating slit: $z$4, $z$5; without slit, values are $z$6, $z$7 (Haque et al., 2024).

3. Microstructure of Graphene and Epitaxial Film Quality

Graphene microstructure, specifically grain size and density of grain boundaries ($z$8), directly governs local damage accumulation and resultant BaTiO$z$9 crystallinity. Wet PMMA transfer yields initially defect-free graphene (Raman $k \approx 2\pi/a_\mathrm{sub}$0; RMS roughness $k \approx 2\pi/a_\mathrm{sub}$1 post-anneal). Two regimes are observed:

• Small-grained graphene ($k \approx 2\pi/a_\mathrm{sub}$2, $k \approx 2\pi/a_\mathrm{sub}$3) accumulates more damage and leads to BaTiO$k \approx 2\pi/a_\mathrm{sub}$4 X-ray rocking curve (RC) FWHM $k \approx 2\pi/a_\mathrm{sub}$5.
• Large-grained graphene ($k \approx 2\pi/a_\mathrm{sub}$6, $k \approx 2\pi/a_\mathrm{sub}$7) yields $k \approx 2\pi/a_\mathrm{sub}$8 and RC FWHM $k \approx 2\pi/a_\mathrm{sub}$9, approaching substrate-direct epitaxy ($a_\mathrm{sub}$0). Grain boundaries act as “soft spots,” visible as defect peaks in Raman line scans every $a_\mathrm{sub}$1300 μm.

The correlation between grain-boundary density and film RC width can be modeled semi-empirically: $a_\mathrm{sub}$2; $a_\mathrm{sub}$3 represents a damage-sensitivity coefficient (Haque et al., 2024).

4. Exfoliation and Transfer of Freestanding BaTiO$a_\mathrm{sub}$4 Membranes

The exfoliation protocol involves sequential Ni stressor deposition (100 nm e-beam evaporation, 1 μm sputtering) atop BaTiO$a_\mathrm{sub}$5/graphene/SrTiO$a_\mathrm{sub}$6, followed by thermal release tape (TRT) application and heating ($a_\mathrm{sub}$790–110 °C). Exfoliation is driven by the difference in work of adhesion ($a_\mathrm{sub}$8) between interfaces:

• $a_\mathrm{sub}$9
• $_3$0

Resultant BTO/Ni/TRT stacks (up to $_3$1) can be transferred onto SiO$_3$2/Si substrates, retaining single-crystalline BTO as validated by cross-sectional HRTEM, STEM-EDS, and preservation of lattice fringes for Ba, Ti, O, and Ni (Haque et al., 2024).

Strain state evolves from initial compressive in-plane ($_3$3) in BTO/STO ($_3$4; $_3$5), toward near-relaxation post-exfoliation ($_3$6, $_3$7). Bilayer graphene further weakens substrate bonding, fostering membrane release and expanded unit cell volume ($_3$8 after exfoliation).

5. Functional Characterization and Integration Potential

Ferroelectricity is retained in freestanding BTO, as shown by piezoresponse force microscopy (PFM) domain writing: clear 180$_3$9 phase contrast between $_3$0 V and $_3$1 V poled regions, amplitude variation indicative of ferroelectric switching. SS-PFM measurements display hysteresis below Curie temperature ($_3$2C) with coercive bias $_3$3V; loops collapse above $_3$4, confirming a ferroelectric transition (Haque et al., 14 Nov 2025).

Obtaining large-area, single-crystal, freestanding oxide membranes via remote epitaxy and PLD (at $_3$5C) enables direct integration with CMOS or flexible platforms. Functional properties—high dielectric constant, piezo- and ferroelectric response—combine with mechanical flexibility and reusable substrates, supporting scalable manufacturing for oxide memories, sensors, and MEMS.

6. Mechanistic Insights, Process Optimization, and AI Workflow

Graphene’s damage during PLD arises from combined ballistic (plume-induced) and oxidative (chemical etching) mechanisms. The critical ballistic threshold is $_3$6eV per C-atom, exceeded by “fast” plume species but not “slow” ones. Raman spectroscopy tracks damage evolution ($_3$7 increases exponentially with $_3$8 or $_3$9), and MD modeling describes oxidative degradation as $T \sim 700\,^\circ$0 with $T \sim 700\,^\circ$1–$T \sim 700\,^\circ$2.

AI-human collaborative workflows (HAIC), integrating hypothesis generation, policy updates, and active-learning driven PLD experiments, accelerate mapping of growth conditions, mitigation strategies, and mechanistic understanding. For BaTiO$T \sim 700\,^\circ$3/graphene, HAIC facilitated identification of growth windows, refined Raman-based defect metrics, and generalized the two-step Ar/O$T \sim 700\,^\circ$4 process to other oxide/2D material systems (Haque et al., 14 Nov 2025). A plausible implication is that this paradigm applies broadly, provided in situ diagnostics and adaptive protocols are employed.

7. Generalization to Other Oxide/2D Material Systems

The interplay of high-energy plume impact and reactive oxidation applies to a wide spectrum of oxygen-rich oxide PLD systems (e.g., SrTiO$T \sim 700\,^\circ$5, VO$T \sim 700\,^\circ$6, HfO$T \sim 700\,^\circ$7). The two-step inert/oxidizing-atmosphere PLD protocol is directly transferable to systems where stoichiometric oxide growth and 2D interlayer preservation are mutually antagonistic. The mechanistic insights from BaTiO$T \sim 700\,^\circ$8/graphene remote epitaxy thus extend to a broad class of functional heterostructures, contingent on appropriately balancing process parameters and exploiting AI-driven optimization strategies (Haque et al., 14 Nov 2025).