Localized Film–Substrate Reactions

• Localized film–substrate reaction is a process at the interface where a thin film undergoes chemical or physical transformations due to substrate interactions, resulting in spatially confined microstructural alterations.
• Mechanisms include thermally induced reactions, diffusively driven interfacial changes, and defect-mediated buckling that enable precise control of film morphology and compositional gradients.
• Modeling these reactions involves coupled equations for species transport, heat transfer, reaction kinetics, and mechanical response to predict features like halo formation and phase transformation.

A localized film–substrate reaction refers to a chemically or physically mediated process occurring at an interface zone where a thin film is in close proximity to, or in direct contact with, a structurally or compositionally distinct substrate. These reactions generate microstructural, compositional, or morphological alterations that are spatially confined—often determined by externally patterned regions, local energy input, or intrinsic material inhomogeneity. Such phenomena underpin the functionality and reliability of a wide spectrum of device architectures and provide critical pathways for interface engineering, phase selection, catalytic activation, and direct-write patterning in various material systems.

1. Types of Localized Film–Substrate Reactions

Localized film–substrate reactions are triggered by a variety of mechanisms, broadly classifiable by the driver of localization:

• Thermally Induced Chemical Reaction: Focused energy input (e.g., laser irradiation) on coated substrates can locally exceed critical temperatures for decomposition or reaction, selectively transforming the film or forming new phases with the substrate (Afrin et al., 2016).
• Diffusively Driven Interfacial Reaction: Exposure of distinct substrate regions (e.g., via FIB-milled apertures) initiates interfacial reactions—such as metal–silicide formation—only at selected sites. The resulting reaction zones (“halos”) are governed by the interplay of local diffusion and reaction kinetics (Peddiraju et al., 2023, Peddiraju et al., 10 Jan 2026).
• Mechanically Localized Buckling and Patterning: Arrays of defects in substrates (e.g., patterned holes) promote spatially localized delamination, wrinkling, or folding of films through the local reduction of constraint, reversing typical bifurcation sequences (Liao et al., 2017).
• Interfacial Charge Transfer and Bonding Modulation: Thermal annealing or chemical treatment modulates bonding strength, strain, and charge transfer at the monolayer–substrate interface, producing heterogeneous property landscapes (Su et al., 2016).
• Substrate-Mediated Catalytic Activation: Ultrathin films supported on reactive substrates can promote unique dissociative adsorption or redox chemistry confined to the interface, which is absent in thicker films or on inert supports (Li et al., 2016).

2. Governing Equations and Modeling Approaches

Localized film–substrate reactions are described by coupled partial differential equations governing:

• Species transport: $\partial c_i/\partial t = D_i \nabla^2 c_i$ for each relevant species (metal atom, oxidant, etc.) in the film or reaction zone. Spatial gradients are especially steep and localized adjacent to active interfaces (Peddiraju et al., 10 Jan 2026).
• Heat transfer: $$\rho C_p \partial T/\partial t = \nabla\cdot[k(T)\nabla T]$$ in the solid domain, directly coupled to exothermic or endothermic reactions, and to substrate thermal properties (Afrin et al., 2016).
• Chemical reaction kinetics: Often zero- or first-order Arrhenius expressions are used, e.g., $k(T) = (R T / 2\pi M_{O_2})^{1/2} E \exp(-E_a/RT)$ for thermally triggered decomposition (Afrin et al., 2016); interfacial phase transformations (e.g., silicide growth) are described by Stefan-type moving boundary conditions incorporating local mass flux (Peddiraju et al., 10 Jan 2026).
• Mechanical response: For patterned substrates, elastic energy functionals include bending (film), substrate deformation, and work of imposed stress; defect sites modulate local stiffness and nucleate buckling/folding at lower critical strains (Liao et al., 2017).
• Electronic structure contributions: For catalytic interfaces, first-principles methods (DFT) quantify charge transfer, preferred adsorption geometries, and activation barriers at specific local configurations (Li et al., 2016).

Numerical implementation varies: finite-volume methods with moving meshes for time-evolving geometries (Afrin et al., 2016), axisymmetric analytical models for reaction-diffusion coupling (Peddiraju et al., 10 Jan 2026), and large-scale atomistic or first-principles calculations for ultrathin film interfaces (Li et al., 2016).

3. Microstructural and Property Outcomes

Distinct spatial features result from localized film–substrate reactions, as demonstrated in controlled experiments and simulations:

Localization Type	Major Microstructural Trait	Typical Characteristic Length Scale
Thermally-induced film decomposition	Paint crater, pigment parcel ejection	$\sim$10–50 µm (laser spot)
Silicide/halo formation (FIB patterning)	Core–shell (silicide + Ag/halo) rings	$\sim$0.5–2 µm (halo width)
Substrate-defect-induced folding	Periodic ridges/folds above holes	$\sim d$ (hole diameter)
Monolayer 2D interface reorganization	Inhomogeneous strain/doping domains	$\sim$1–10 µm (PL map resolution)
Catalytically enhanced O₂ dissociation	Localized interfacial O atoms	Sub-nm (interface region)

Key observations include:

• Sharp chemical gradients: Near FIB-milled apertures, composition transitions from nearly pure Ag in the “halo” to equilibrium Ag–Cu spinodal domains in the far field (Peddiraju et al., 2023, Peddiraju et al., 10 Jan 2026).
• Pattern reversibility and tunability: In elastic bilayers, the presence, geometry, and spacing of substrate holes reverses the typical wrinkle → fold sequence, enabling the design of local vs. global surface reliefs (Liao et al., 2017).
• Non-classical kinetics: Lateral extent of reaction/diffusion zones (e.g., halo width $W_H$) scale as $t^{0.46}$ (halo) and $t^{0.30}$ (core silicide) with time, rather than $t^{1/2}$, due to nontrivial coupling between mass balance and V-shaped growth geometry (Peddiraju et al., 10 Jan 2026).
• Defect-driven electronic and optoelectronic inhomogeneity: Annealing-induced variations in strain, doping, and adhesion in 2D films create mesoscale spatial heterogeneity in PL and carrier transport (Su et al., 2016).

4. Kinetic Theory and Modeling of Halo Formation

The evolution of locally reacted zones, particularly in thin film–substrate systems displaying phase transformations (e.g., silicide formation upon annealing), requires incorporating both mass transport through the film and phase boundary movement.

The key kinetic model (Peddiraju et al., 10 Jan 2026) includes:

• Species conservation: Initial Cu inventory in a local region partitions into Cu₃Si product and a Cu-depleted/Ag-rich “halo”, requiring explicit accounting for geometric volumes ($V_F$, $V_P$, $V_H$) and concentration differences.
• Diffusive transport: The Cu gradient across the halo drives the net flux $J_{\rm Cu} = -D_{\rm Cu}(C_F - C_I)/W_H$.
• Modified Stefan condition: The reaction front’s propagation rate is determined by matching the diffusive supply of Cu to the growing silicide core, yielding:

$$\frac{dW_P}{dt} = \frac{D_{\rm Cu}(C_F-C_I)/W_H}{(C_F-C_I)2h_F + C_F W_P \tan\theta} W_P$$

where $W_P$ is the product width, $h_F$ the film thickness, and $\theta$ the interface angle.

Predicted exponents and experimental results (e.g., $W_P \propto t^{0.29}$, $W_H \propto t^{0.46}$) deviate from classical Fickian kinetics due to the strong constraint imposed by finite film thickness and geometric coupling between product and diffusive zone (Peddiraju et al., 10 Jan 2026). This modeling approach is broadly transferable to any finite-thickness system with an advancing reaction front and restricted axial diffusion.

5. Substrate and Interface Selection for Reaction Control

The spatial confinement and nature of localized reactions are fundamentally determined by substrate selection, interfacial engineering, and external patterning:

• Reactive substrate exposure (e.g., Si exposed through Si₃N₄ by FIB) leads to robust local phase transformations (Cu₃Si, Ag-rich halo), while inert substrates (e.g., Si₃N₄ windows) suppress such reactions—only thin Ag enrichment is observed (Peddiraju et al., 2023).
• Patterned surface energy/modulus modification: Substrate holes introduce strain-relief sites, localizing buckling in stiff-film/soft-substrate systems (Liao et al., 2017).
• Catalytic function tuning: In ultrathin oxide–metal systems, catalytic activation (e.g., O₂ dissociation) and the formation of reactive interfacial species (e.g., O_int, OH) occur only for films below a critical thickness and on supports that accommodate charge transfer (Li et al., 2016).

A plausible implication is that by judicious substrate and interface design, it is possible to spatially modulate device properties (plasmonic, catalytic, magnetic, electronic) at the sub-micron to nanometer scale, introducing spatially varying functional “hot spots” or catalytically active domains.

6. Applications and Engineering Implications

Localized film–substrate reactions present a versatile strategy for controlled material patterning and interface property manipulation, enabling:

• Direct-write patterning and selective phase formation: FIB milling followed by annealing enables phase and composition modulation with sub-micron precision (Peddiraju et al., 2023, Peddiraju et al., 10 Jan 2026).
• Surface catalysis enhancement: Ultrathin oxide films on metals can harness interface-induced reactivity for environmental catalysis or green synthesis applications (Li et al., 2016).
• Strain/charge patterning in 2D materials: Thermal annealing protocols can be engineered to produce controlled landscapes of strain and carrier doping, optimizing optoelectronic device function (Su et al., 2016).
• Controllable wrinkling/folding: Surface morphologies can be engineered by defect patterning in elastomeric substrates to achieve customizable wave, ridge, or fold states over device-relevant areas (Liao et al., 2017).
• Microstructural tailoring for functional device modules: The ability to generate halos, decomposed zones, or reaction-driven ejecta supports integration of localized functional regions in sensors, interconnects, Schottky-junction arrays, or catalytic spots (Afrin et al., 2016, Peddiraju et al., 2023).

The extensive tunability of the localization—via anneal parameters, substrate geometry, film thickness, or interfacial chemistry—provides a robust platform for advancing nanofabrication, interface science, and micro/nanoscale device engineering.