Search

Hafnia-based ferroelectrics hold promise for nonvolatile ferroelectric memory devices. However, the high coercive field required for polarization switching remains a prime obstacle to their practical applications. A notable reduction in coercive field has been achieved in ferroelectric Hf(Zr)$_{1+x}$O$_2$ films with interstitial Hf(Zr) dopants [Science 381, 558 (2023)], suggesting a less-explored strategy for coercive field optimization. Supported by density functional theory calculations, we demonstrate the $Pca2_1$ phase, with a moderate concentration of interstitial Hf dopants, serves as a minimal model to explain the experimental observations, rather than the originally assumed rhombohedral phase. Large-scale deep potential molecular dynamics simulations suggest that interstitial defects promote the polarization reversal by facilitating $Pbcn$-like mobile 180$^\circ$ domain walls. A simple pre-poling treatment could reduce the switching field to less than 1 MV/cm and enable switching on a subnanosecond timescale. High-throughput calculations reveal a negative correlation between the switching barrier and dopant size and identify a few promising interstitial dopants for coercive field reduction.

Since the first report of ferroelectricity in nanoscale HfO$_2$-based thin films in 2011, this silicon-compatible binary oxide has quickly garnered intense interest in academia and industry, and continues to do so. Despite its deceivingly simple chemical composition, the ferroelectric physics supported by HfO$_2$ is remarkably complex, arguably rivaling that of perovskite ferroelectrics. Computational investigations, especially those utilizing first-principles density functional theory (DFT), have significantly advanced our understanding of the nature of ferroelectricity in these thin films. In this review, we provide an in-depth discussion of the computational efforts to understand ferroelectric hafnia, comparing various metastable polar phases and examining the critical factors necessary for their stabilization. The intricate nature of HfO$_2$ is intimately related to the complex interplay among diverse structural polymorphs, dopants and their charge-compensating oxygen vacancies, and unconventional switching mechanisms of domains and domain walls, which can sometimes yield conflicting theoretical predictions and theoretical-experimental discrepancies. We also discuss opportunities enabled by machine-learning-assisted molecular dynamics and phase-field simulations to go beyond DFT modeling, probing the dynamical properties of ferroelectric HfO$_2$ and tackling pressing issues such as high coercive fields.

In the realm of ferroelectric memories, HfO2-based ferroelectrics stand out because of their exceptional CMOS compatibility and scalability. Nevertheless, their switchable polarization and switching speed are not on par with those of perovskite ferroelectrics. It is widely acknowledged that defects play a crucial role in stabilizing the metastable polar phase of HfO2. Simultaneously, defects also pin the domain walls and impede the switching process, ultimately rendering the sluggish switching of HfO2. Herein, we present an effective strategy involving acceptor-donor co-doping to effectively tackle this dilemma. Remarkably enhanced ferroelectricity and the fastest switching process ever reported among HfO2 polar devices are observed in La3+-Ta5+ co-doped HfO2 ultrathin films. Moreover, robust macro-electrical characteristics of co-doped films persist even at a thickness as low as 3 nm, expanding potential applications of HfO2 in ultrathin devices. Our systematic investigations further demonstrate that synergistic effects of uniform microstructure and smaller switching barrier introduced by co-doping ensure the enhanced ferroelectricity and shortened switching time. The co-doping strategy offers an effective avenue to control the defect state and improve the ferroelectric properties of HfO2 films.

Ferroelectrics and ionic conductors are important functional materials, each supporting a plethora of applications in information and energy technology. The underlying physics governing their functional properties is ionic motion, and yet studies of ferroelectrics and ionic conductors are often considered separate fields. Based on first-principles calculations and deep-learning-assisted large-scale molecular dynamics (MD) simulations, we report ferroelectric-switching-promoted oxygen ion transport in HfO$_2$, a wide-band-gap insulator with both ferroelectricity and ionic conductivity. Applying a unidirectional bias can activate multiple switching pathways in ferroelectric HfO$_2$, leading to polar-antipolar phase cycling that appears to contradict classical electrodynamics. This apparent conflict is resolved by the geometric-quantum-phase nature of electric polarization that carries no definite direction. Our MD simulations demonstrate bias-driven successive ferroelectric transitions facilitate ultrahigh oxygen ion mobility at moderate temperatures, highlighting the potential of combining ferroelectricity and ionic conductivity for the development of advanced materials and technologies.

Recent discovery of HfO2-based and nitride-based ferroelectrics that are compatible to the semiconductor manufacturing process have revitalized the field of ferroelectric-based nanoelectronics. Guided by a simple design principle of charge compensation and density functional theory calculations, we discover HfO2-like mixed-anion materials, TaON and NbON, can crystallize in the polar Pca21 phase with a strong thermodynamic driving force to adopt anion ordering spontaneously. Both oxynitrides possess large remnant polarization, low switching barriers, and unconventional negative piezoelectric effect, making them promising piezoelectrics and ferroelectrics. Distinct from HfO2 that has a wide band gap, both TaON and NbON can absorb visible light and have high charge carrier mobilities, suitable for ferroelectric photovoltaic and photocatalytic applications. This new class of multifunctional nonperovskite oxynitride containing economical and environmentally benign elements offer a platform to design and optimize high-performing ferroelectric semiconductors for integrated systems.

The multicomponent oxide solid solution is a versatile platform to tune the delicate balance between competing spin, charge, orbital, and lattice degrees of freedom for materials design and discovery. The development of compositionally complex oxides with superior functional properties has been largely empirical and serendipitous, in part due to the exceedingly complex chemistry and structure of solid solutions that span a range of length scales. The classical molecular dynamics (MD), as a powerful statistical method to investigate materials properties over large spatial and temporal scales, often plays a secondary role in computer-aided materials discovery because of the limited availability and accuracy of classical force fields. Here, we introduce the strategy of ``modular developing deep potential" (ModDP) that enables a systematic development and improvement of deep neural network-based model potential, termed as deep potential, for complex solid solutions with minimum human intervention. The converged training database associated with an end-member material is treated as an independent module and is reused to train the deep potential of solid solutions via a concurrent learning procedure. We apply ModDP to obtain classical force fields of two technologically important solid solutions, Pb$_x$Sr$_{1-x}$TiO$_3$ and Hf$_x$Zr$_{1-x}$O$_2$. For both materials systems, a single model potential is capable of predicting various properties of solid solutions including temperature-driven and composition-driven phase transitions over a wide range of compositions. In particular, the deep potential of Pb$_x$Sr$_{1-x}$TiO$_3$ reproduces a few known topological textures such as polar vortex lattice and electric dipole waves in PbTiO$_3$/SrTiO$_3$ superlattices, paving the way for MD investigations on the dynamics of topological structures in response to external stimuli., Comment: 32 pages, 9 figures

Defects such as oxygen vacancy are widely considered to be critical for the performance of HfO2-based devices, and yet atomistic mechanisms underlying various exotic effects such as wake-up and fluid imprint remain elusive. Here, guided by a lattice-mode-matching criterion, we systematically study the phase transitions between different polymorphs of hafnia under the influences of neutral and positively charged oxygen vacancies by mapping out the minimum energy pathways using a first-principles-based variable-cell nudged elastic band technique. We find that the positively charged oxygen vacancy can substantially promote the transition of various nonpolar phases to the polar phase kinetically, enabled by a transient high-energy tetragonal phase and extreme charge-carrier-inert ferroelectricity of the polar $Pca2_1$ phase. The intricate coupling between structural polymorphism kinetics and the charge state of the oxygen vacancy has important implications for the origin of ferroelectricity in HfO$_2$-based thin films as well as wake-up, fluid imprint, and inertial switching.

Ferroelectric domains are mesoscale structures that mediate between synchronized atomic-scale ion displacements and switchable macroscopic polarization. Here, we evaluated the randomness of the domain distribution at the onset of ferroelectricity. First-principle calculations combined with atomic-scale imaging demonstrate that oxygen vacancies that serve as pinning sites for the ferroic domain walls remain immobile above the Curie temperature. Thus, upon cooling to a ferroelectric state, these oxygen vacancies dictate reproducible domain-wall patterning. Domain-scale imaging with variable-temperature piezoresponse force microscopy confirmed the memory effect, questioning the spontaneity of domain distribution under thermotropic transitions.

The discovery of two-dimensional (2D) ferroelectrics with switchable out-of-plane polarization such as monolayer $\alpha$-In$_2$Se$_3$ offers a new avenue for ultrathin high-density ferroelectric-based nanoelectronics such as ferroelectric field effect transistors and memristors. The functionality of ferroelectrics depends critically on the dynamics of polarization switching in response to an external electric/stress field. Unlike the switching dynamics in bulk ferroelectrics that have been extensively studied, the mechanisms and dynamics of polarization switching in 2D remain largely unexplored. Molecular dynamics (MD) using classical force fields is a reliable and efficient method for large-scale simulations of dynamical processes with atomic resolution. Here we developed a deep neural network-based force field of monolayer In$_2$Se$_3$ using a concurrent learning procedure that efficiently updates the first-principles-based training database. The model potential has accuracy comparable with density functional theory (DFT), capable of predicting a range of thermodynamic properties of In$_2$Se$_3$ polymorphs and lattice dynamics of ferroelectric In$_2$Se$_3$. Pertinent to the switching dynamics, the model potential also reproduces the DFT kinetic pathways of polarization reversal and 180$^\circ$ domain wall motions. Moreover, isobaric-isothermal ensemble MD simulations predict a temperature-driven $\alpha \rightarrow \beta$ phase transition at the single-layer limit, as revealed by both local atomic displacement and Steinhardt's bond orientational order parameter $Q_4$. Our work paves the way for further research on the dynamics of ferroelectric $\alpha$-In$_2$Se$_3$ and related systems., Comment: 27pages, 10figures