Kinetic Monte Carlo simulation reveals defect charge accumulation favoring ferroelectric phase formation in Hf1−<i>x</i>Zr<i>x</i>O2 thin films

L Luis Azevedo Antunes (Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,) S Sebastian Obernberger (Institut National de la Recherche Scientifique (INRS), Center of Energy, Materials, and Telecommunications , Varennes, Quebec J3X 1P7,) P Paul Schwermer (Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,) J Joshua Hintz (Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,) R Richard Ganser (Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,) U Uwe Schroeder (NaMLab gGmbH 2 , Noethnitzer Strasse 64 a, 01187 Dresden,) T Thomas Mikolajick (1 Chair of Nanoelectronics, TU Dresden, 01187 Dresden, Germany) A Alfred Kersch (Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,)

Abstract

Charged oxygen vacancies are thought to be responsible for fatigue effects in HfO2- and ZrO2-based ferroelectrics, while also supporting the formation of the ferroelectric phase. We investigate the possible influence of intrinsic electric fields generated by predominantly doubly positively charged oxygen vacancies that accumulate near the electrode and violate charge neutrality. Our calculations are based on a meticulous kinetic Monte Carlo simulation, which simultaneously treats the movement of electrons, the diffusion of oxygen defects, and their interactions that change the charge state in a self-consistent manner. It is shown that with a realistic vacancy concentration, electric fields of several MVcm−1 can form within 2nm of the electrode, supporting the formation of the ferroelectric phase. This effect becomes significant for thin films below 10nm. The strength of the effect depends on the position of the defect energy level relative to the valence band offset. This dependency could be used to control the effect.

Article Details

Volume / Issue Vol. 137, Issue 22
Published June 14, 2025
ISSN 0021-8979
Publisher American Institute of Physics

Journal Info

Journal of Applied Physics

American Institute of Physics

ISSN: 0021-8979 Physical Sciences

Authors (8)

L

Luis Azevedo Antunes

Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,

S

Sebastian Obernberger

Institut National de la Recherche Scientifique (INRS), Center of Energy, Materials, and Telecommunications , Varennes, Quebec J3X 1P7,

P

Paul Schwermer

Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,

J

Joshua Hintz

Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,

R

Richard Ganser

Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,

U

Uwe Schroeder

NaMLab gGmbH 2 , Noethnitzer Strasse 64 a, 01187 Dresden,

T

Thomas Mikolajick

1 Chair of Nanoelectronics, TU Dresden, 01187 Dresden, Germany

A

Alfred Kersch

Department of Applied Sciences and Mechatronics, Munich University of Applied Sciences 1 , Lothstr. 34, 80335 Munich,