Charge redistribution at metal-ZrO2 interfaces: A combined DFT and continuum electrostatic study

Nanoscale metallic inclusions (NMIs) are commonly observed within oxide scales formed during high-temperature oxidation, revealing the existence of chemical and electronic heterogeneity beyond conventional corrosion theories that assume homogeneous, fully oxidized films. Using tetragonal zirconia (tZrO2) facing a series of face-centered cubic (fcc) metals as the model system, this work investigates the short-range and long-range charge redistributions across metal-oxide interfaces by coupling density functional theory (DFT) calculations with continuum modeling. We show that metal-oxide contact induces a short-range charge redistribution confined to a few atomic layers and a long-range redistribution of space charge that can extend over macroscopic distances within weakly doped oxides. DFT calculations show that the short-range redistribution is dominated by metal induced gap states (MIGS) in tZrO2 facing noble metals like Au and Ag, and by chemical bonding in tZrO2 facing active metals like Al. DFT-informed continuum theoretical analysis shows that the range of space-charge redistribution is governed by the doping level of tZrO2, and that the Schottky barrier height (SBH) exhibits a stronger dependence on the metal work function than the doping level. Both the short-range and long-range charge redistributions can alter the transport of charge carriers via their associated electric fields, extending several nm to hundreds of nm from the interface, depending on the doping concentrations, suggesting possible heterogeneous oxide growth caused by NMIs.

💡 Research Summary

This paper investigates charge redistribution at metal–zirconia interfaces by coupling first‑principles density functional theory (DFT) with a continuum electrostatic model. The authors focus on tetragonal ZrO₂ (tZrO₂), a representative high‑temperature oxidation product of Zr‑based nuclear fuel cladding, interfaced with three face‑centered cubic (fcc) metals—Au, Ag, and Al—that span a wide range of work functions, chemical reactivity, and valence‑electron character.

Methodology
Metal (111) surfaces and the tZrO₂ (100) surface are combined into slab models with minimal lattice mismatch. DFT calculations (VASP, PAW‑PBE, DFT+U for Zr 4d) provide relaxed atomic structures, electronic band alignments, work functions, adhesion energies, and Bader charge analyses. To capture phenomena beyond the nanometer scale accessible to DFT, the authors embed DFT‑derived interfacial properties into a self‑consistent Poisson solver that treats the metal as an ideal electron reservoir and the oxide as a lightly doped ionic semiconductor. The continuum model yields space‑charge density, electrostatic potential profiles, Schottky barrier height (SBH), and depletion‑layer width (DLW) as functions of metal work function and oxide doping concentration.

Short‑Range (Atomic‑Scale) Findings

• Noble metals (Au, Ag): Weak metal–oxygen bonding leads to metal‑induced gap states (MIGS) that penetrate a few atomic layers into tZrO₂. MIGS dominate the interfacial dipole, producing charge transfers of 0.2–0.8 e nm⁻².
• Active metal (Al): Strong Al–O covalent bonding governs charge redistribution; MIGS play a minor role. The interfacial charge is still of comparable magnitude but originates mainly from chemical bond formation.
• The short‑range redistribution is confined to ≤ 3 Å on the oxide side and a few screening lengths in the metal.

Long‑Range (Continuum) Findings

• The continuum model shows that the spatial extent of space‑charge redistribution is set primarily by the bulk doping level of tZrO₂. For carrier concentrations of 10¹⁶–10¹⁸ cm⁻³, the Debye length ranges from ~10 nm to >200 nm, leading to depletion layers that can span tens to hundreds of nanometers.
• Metal work function strongly influences the SBH: high‑Φ Au (5.15 eV) and Ag (4.47 eV) raise the SBH to 1.2–1.5 eV, while low‑Φ Al (4.02 eV) lowers it to ~0.8 eV. The SBH dependence on doping is weaker but not negligible.
• Global charge neutrality forces the net interfacial charge (from DFT) to be balanced by the integrated space‑charge within the oxide, linking the atomic‑scale dipole to the macroscopic electric field.

Implications for NMIs and Oxide Growth
Nanoscale metallic inclusions (NMIs) embedded in oxide scales create internal metal–oxide interfaces that generate both a localized dipole and an extended space‑charge region. The resulting electric fields modify the drift and diffusion of charged species (electrons, holes, oxygen vacancies) over distances far exceeding the immediate interface. Consequently, transport rates become spatially heterogeneous, providing a mechanistic explanation for experimentally observed non‑uniform oxide thickening and accelerated local oxidation near NMIs.

Significance of the Integrated Approach
By unifying DFT‑level interfacial electronic structure with continuum electrostatics, the study bridges the gap between atomistic and macroscopic descriptions that have traditionally been treated separately. This framework enables quantitative prediction of how metal properties (work function, bonding character) and oxide doping jointly dictate SBH, depletion width, and electric‑field profiles. The insights are directly relevant to designing corrosion‑resistant alloys, optimizing protective oxide layers, and engineering nano‑composite materials where metal‑oxide interfaces dominate performance.

In summary, the paper demonstrates that metal‑tZrO₂ interfaces exhibit a dual‑scale charge redistribution: a short‑range, metal‑specific dipole driven by MIGS or chemical bonding, and a long‑range space‑charge region governed by oxide doping. Both contributions generate electric fields that can extend from a few nanometers to hundreds of nanometers, profoundly influencing charge‑carrier transport and heterogeneous oxide growth in the presence of nanoscale metallic inclusions.