Unraveling the origins of ferroelectricity in doped hafnia through carrier-mediated phase transitions

Introduction
After the discovery of ferroelectricity in doped HfO2 and (Hf,Zr)O2 (HZO) thin films1,2, they have drawn widespread attention for the next-generation nanoscale devices. This interest arises from their full compatibility with the complementary metal-oxide-semiconductor (CMOS) process, coupled with their excellent retention and scalability3,4,5. Unfortunately, the stable phase of bulk HfO2 and ZrO2 at room temperature and ambient pressure is monoclinic phase (m, P21/c)6,7. At high-temperature, the stable phase transforms into the tetragonal (t, P42/nmc) phase and cubic phase (c, Fm(bar{3}) m), while at high pressure, it changes into the antipolar orthorhombic phase (o, Pbca)6. Since the orthorhombic ferroelectric phase (f, Pca21) is absent in phase diagrams, it is crucial to determine the possible thermodynamic and kinetic driving factors behind the formation of this unexpected ferroelectric f-phase.
Doping with various cations or nitrogen is recognized as a pivotal method for inducing and controlling ferroelectricity in HfO2 and HZO thin films8,9,10,11,12,13, encompassing synergistic interplay among factors such as oxygen vacancies14,15,16, surface energy17,18, and stress and strain19,20. Theoretical predictions21,22 and experiments11,23,24 have verified that dopants with larger ionic radius and lower electronegativity are more effective in stabilizing ferroelectric f-phase at appropriate concentrations. Excessive doping reduces the f-phase fraction by weakening the transition from the t-phase during rapid thermal annealing, favoring the dominance of the t-phase25,26. Notably, a closer examination of the dopants that stabilize the f-phase of hafnia reveals a common characteristic: they are typically p-type dopants (Ca, Sr, Al, Y, La, Gd, etc), which inevitably introduce holes into the crystal. Interestingly, experiments unexpectedly found that ferroelectricity in n-type Ta-doped HfO2 films27 occurs only at high concentrations, with no current explanation for this observation. Unfortunately, the roles played by the carriers introduced by these dopants are often overlooked.
Ferroelectric (polar) metals have recently become a research hotspot due to their potential to elucidate the mechanism underlying the coexistence of two seemingly incompatible properties: ferroelectricity (polarity) and metallicity28,29,30,31. Experimental findings have identified a critical electron concentration in n-type BaTiO3-based ferroelectric-like metals32. Theoretical insights suggest that the persistence of polar distortion is linked to the incomplete screening of Coulomb forces within this concentration range33. Other interaction mechanisms also include the meta-screening effect and chemical bonding, which anomalously enhance the ferroelectricity of perovskites with carrier doping34,35,36. Recent studies have realized the crucial role played by holes introduced by dopants37 or interface engineering38 in stabilizing f-phase HfO2. Nevertheless, the true microscopic mechanisms through which carriers introduced by dopants stabilize the metastable f-phase remain elusive.
In this letter, we perform a first-principles analysis of the evolution of crystal structure, electronic structure, and chemical bonds in the various phases of hafnia under carrier electrostatic doping. Notably, our analysis intentionally disregards interference from intrinsic properties of the dopant, such as ionic radius and electronegativity. Our findings reveal that carrier-doped hafnia undergoes intricate phase transitions, posing challenges in unequivocally identifying the ferroelectric nature of the material. The presence of multiple coexisting phases, driven by small energy differences, further complicates this identification. Remarkably, the phase transitions contribute to the enhancement of f-phase stability at a specific concentration in the case of p-type doping, while n-type doping proves to be more challenging. To elucidate this phenomenon, we examine the interplay between carriers and ferroelectricity in terms of both screening effects and chemical bonds.
For p-type doping, polar triply coordinated oxygen (OIII) and nonpolar quadruply coordinated oxygen (OIV) assume distinctive roles, where the presence of a hole on the OIV atom significantly impacts the total energy of each phase, whereas a hole on OIII atoms stabilizes the polar distortion by reinforcing the Hf-OIII covalent bond. Consequently, hafnia undergoes a sequence of ground-state phase transitions (m → f → c) with an increasing hole concentration. This evolving pattern of phase transition closely parallels observations in La-doped and Y-doped HfO2 thin films12,23. Conversely, n-type doping follows an opposite mechanism due to the strong screening effect and the instability of the Hf-OIV covalent bond. This drives the metastable f-phase towards destabilization and ultimately leads to phase transitions to the stable antipolar Pbcn phase with increasing electron concentration. This mechanism can be indirectly confirmed by Ta-doped HfO2 films27, since the remanent polarization Pr is almost identical to theoretical predictions.
These unconventional concepts might appear perplexing at first, particularly considering the inherent challenge of measuring ferroelectric properties in metals. To address this, we draw comparisons between the electronic characteristics of SiO2/HfO2-based heterostructures with and without TiN electrodes. Our speculation is grounded in the idea that carriers play a preferential role in fostering the formation of the f-phase. The rationale behind this lies in the fact that, when capped with electrodes, carriers are directed toward the electrode material, ensuring ferroelectric insulation. This rationale sheds light on why a 1 nm ultrathin HZO layer without a metal capping layer does not exhibit ferroelectric signatures, but displays such signatures when paired with a metal capping layer39. Our results are consistent well with the existing experimental results, demonstrating that our work resolves the mystery of the origin of ferroelectricity in doped hafnia and suggests manipulating carriers as a universal strategy to tune ferroelectricity in hafnia.
Results and discussions
Rich phase transition under carrier conditions
In Fig. 1a, all OIV atoms in the high-symmetry c-phase and t-phase are equivalent, exhibiting a nonpolar nature. The OIV atoms in the non-centrosymmetric f-phase are equivalent, and polarization in the f-phase originates from the motion of OIII atoms upon polarization switching along the c-direction. The antipolar o-phase can be conceptualized as the OIII polarization reversal of neighboring unit cells in the f-phase. Consequently, the energy difference between the two phases is minimized, allowing for a facile transformation from the o-phase to the f-phase under the influence of an external electric field40. Due to this characteristic, the o-phase is often overlooked in theoretical calculations. Interestingly, although the lower symmetry m-phase features polar OIII atoms and boasts the lowest energy, its overall behavior remains nonpolar. This nonpolar nature is attributed to its inability to form a long-range polar order arrangement.

a Crystal structures of the bulk HfO2 with different phases, red and blue balls indicate quadruply coordinated (OIV) and triply coordinated (OIII) oxygen, the cyan arrow is the direction of polarization. The total energy with different phases as a function of (b) hole and (c) electron concentration, the inset in (b) shows that the f-phase with the lowest total energy. The lattice parameters (a, b, c, α, β, γ) of the bulk HfO2 as a function of the carrier concentration, (d) hole-doped c (t)-phase, (e) hole-doped f-phase, (f) electron-doped f-phase, (g) hole-doped m-phase. The vertical orange-yellow dashed line indicates that the structural phase transition occurs here. h A conceptual phase diagram of carriers doped of bulk HfO2, the stable crystal structure corresponding to different carrier concentrations is determined based on total energy and lattice parameters, the new Pbcm and Pbcn space groups originate from structural phase transitions of m- and f-phase, respectively.
In Fig. 1b, c, carriers profoundly influence the total energy of the various phases of hafnia. Observing the abrupt changes in the total energy of certain phases with increasing carrier concentration is more intuitively reflected in the sudden alterations of phase lattice parameters. These changes imply the occurrence of structural phase transitions. The multitude of phase transitions, coupled with the inherent complexity of the original multiple-phase structures, adds intricacy to the analysis of doped hafnia. In general, hole doping and electron doping contribute to increased stability and instability of the phase structure, respectively.
In Fig. 1b, it is evident that phases with higher symmetry exhibit initially steeper slopes in the decrease of total energy. Consequently, the c-phase, characterized by the highest symmetry, demonstrates the most pronounced decrease slope (Supplementary Fig. 1a). Throughout this progression, the t-phase gradually undergoes transformation into the c-phase. This transformation is visually apparent in Fig. 1d, where the lattice constant of the t-phase gradually diminishes until it aligns with that of the c-phase, occurring at a critical concentration (p = 0.16 h+/f.u.). Surprisingly, the m-phase undergoes a structural phase transition to the new high-symmetry phase (m′, space group: Pbcm) at p = 0.37 h+/f.u., reverting to the ground state with a steeper slope than the c-phase (Supplementary Fig. 1a). This transition is marked by a significant change in the lattice parameter lattice parameter β, shifting abruptly from a non-right angle to a right angle (Fig. 1g). The m′-phase is recognized as an intermediate transition structure in the polarization switching process in many theoretical calculations40,41. In addition, the low-symmetry f-phase and o-phase undergo a structural phase transition to high-symmetry c-phase (Fig. 1e) and m′-phase (Supplementary Fig. 2d) at higher hole concentrations. Upon comparing the total energy descending slopes of the c-phase and m’-phase, it becomes apparent that symmetry under hole doping may not be the primary determinant of the total energy change in each phase.
Comparing the total energies of these doped phase structures, HfO2 undergoes ground-state phase transitions of m → f (o) → c (t) → m′, respectively. This sequence of phase transitions aligns with observations in La-doped HfO2 thin films23. In experiments with increasing La doping concentration, the highest phase fractions in the GIXRD patterns exhibit m → f → c phase structures. The unexpected emergence of the m′-phase following the c-phase at high concentrations naturally would not receive much attention in experiments. Additionally, since our doping approach ignores the physical properties of the dopant atoms, experiments have demonstrated that smaller and larger dopants generally induce the t-phase and c-phase in HfO2, respectively42, so that the t-phase dominates at high Al doping concentrations.
Turning attention to ferroelectric f-phase, the energy difference between the f-phase and m-phase gradually decreases with the increase of hole concentration (Supplementary Fig. 3a). Notably, when the hole concentration reaches 0.16–0.26 h+/f.u. (the inset of Fig. 1b), the total energy of the f-phase becomes the lowest, signifying that suitable hole concentration can indeed stabilize the f-phase. These findings are consistent with theoretical predictions suggesting that p-type elements are the most suitable dopants for promoting ferroelectricity in hafnia, corroborated by the large remanent polarization measured in hafnia films doped with these elements21. It’s worth noting that the actual dopant may have varying effects on the crystal structure of each phase. Theoretical calculations have not identified a scenario where the dopant alone can stabilize the f-phase as the ground state21.
The reduction in lattice constants b and c is not significant in Fig. 1e. However, due to the rapid decrease in lattice constant a, thereby reducing the volume with increasing hole concentration23. These changes in the longest lattice constants are much stronger than that of the two shorter axes lattice constants, which can be found in Al- and Gd-doped HfO2 thin films11. Another crystallographic signature worth noting is the aspect ratio 2a/(b + c), akin to perovskite tetragonality (c/a)43. The aspect ratio exhibits a decreasing trend from 1.04 to 0.95 with increasing hole concentration (Supplementary Fig. 4), mirroring observations observed in La-doped HfO2 thin films, where the aspect ratio decreases from 1.03 to 0.98 with increasing La content23. The aspect ratio is smaller than the ideal value 1, a phenomenon previously attributed to the larger biaxial in-plane tensile stress23,44,45. However, our prediction attributes the change in aspect ratio to the evolution of lattice constants resulting from hole doping. At high hole concentrations, the lattice constants b = c ≠ a (Fig. 1e), indicating that the f-phase eventually transforms into a pseudo-cubic phase (Supplementary Fig. 4).
In contrast to hole doping, the energy difference between the phases with higher symmetry and the m-phase increases with growing electron doping concentration (Supplementary Fig. 3b). Judging from the increase of lattice constants under electron doping (Fig. 1f and Supplementary Fig. 2), the structures of HfO2 become progressively unstable. Particularly noteworthy is the behavior of the t-phase under high electron doping, the long Hf-OIV bonds break, causing a rapid increase in the lattice constant c (c > 10 Å), ultimately resulting in the formation of a ground-state 2D structure (Supplementary Fig. 2a). To distinguish it from the original t-phase, it is named t’-phase (Fig. 1h).
Luo et al.27 discovered n-type Ta-doped HfO2 films exhibited only the presence of the m-phase at low Ta concentrations (<11%, much larger than p-type dopants), without evidence of f-phase/t-phase coexistence. This observation will not appear in p-type doped HfO2 films, however, it supports our prediction that the energy difference between the f-phase/t-phase and the m-phase increases at low electron concentration (Supplementary Fig. 3b), reducing the likelihood of f-phase/t-phase formation. As another empirical demonstration, Xu et al.46 experimentally showed that pentavalent Nb doping (<3%) stabilizes the m-phase of HfO2 thin films more than dopants like N, Zr, Sc, and Y, increasing the m-phase fraction to 100%. At a high Ta doping concentration (16%), the m-phase is suppressed, and the f-phase eventually becomes dominant in the XRD patterns. At this point, the remanent polarization Pr can reach up to 53 μC/cm2 is nearly equal to our predicted 58 μC/cm2. This is also in line with our predictions, that is, the f-phase unexpectedly undergoes a structural phase transition, giving rise to a new ground-state phase with space group Pbcn (f ″).
The f ″-phase of HfO2 primarily emerges under large tensile strain (ε > 3%)47, aligning with the observed large lattice constant at high electron concentrations (Fig. 1f). Examining its crystal structure as shown in Fig. 1h, nonpolar OIV atoms become polar OIII atoms, but the polar direction is diametrically opposite to the original OIII, resulting in an antipolar configuration. Theoretical calculations by Zhou et al.48 reveal that polar-antipolar coupling is the origin of ferroelectricity in HfO2 with large tensile strain. Moreover, Choe et al.49 found that compared with the t-phase, the f ″-phase acts as an intermediate transition phase for polarization switching, with an ultralow energy barrier for the domain wall switching.
The t-phase is widely considered a precursor phase for crystalline annealing of HfO2 ferroelectric films, which induces a kinetically driven t-phase to f-phase transition during the annealing and cooling process. Through climbing image nudged elastic band (NEB) calculations50,51, we found that the energy barrier to be overcome for the transition of the f ”-phase to the f-phase is comparable to that for the transition of the t-phase to the f-phase (Supplementary Fig. 5), and significantly lower than that for the transition of the f ”-phase to the other (c, t, m) phases. This suggests that the antipolar f ”-phase can also transform into the f-phase through thermal annealing. High-temperature annealing under ambient conditions can stabilize the f-phase beyond the critical dopant concentration, as shown by Luo et al.27 Considering the significant challenges with n-type doping, further experimental validation of the polar-antipolar phase transition is undoubtedly necessary.
Carrier regulates ferroelectric properties
In Fig. 2a, we utilized a NEB to calculate the polarization switching barrier (ΔEB). A larger ΔEB indicates more challenging switching polarization, requiring a higher coercive field (Ec)52,53. Figure 2b illustrates that ΔEB first decreases slowly and then decreases rapidly with the increase of hole concentration, indicating that the polarization switching becomes easier, and Ec decreases accordingly. For instance, the Ec drops from 1.35 MV/cm at 5 cat% La to below 1.2 MV/cm at 13 cat% La in the La-doped HfO2 thin films23. Additionally, the Ec = 1.2 MV/cm of the HZO film without La doping decreases to Ec = 0.8 MV/cm of the HZO film with La doping24. In Fig. 2c, ΔEB shows an increasing trend followed by a decreasing trend, indicating that polarization switching becomes increasingly difficult at low electron concentrations, but because OIV tends to be unstable under high electron concentrations, the ΔEB decreases rapidly. Therefore, the coercive field of the 16% Ta-doped HfO2 thin film is relatively small, that is, Ec = 0.8 MV/cm27. Recently, Wang et al.54 discovered rhombohedral ferroelectric Hf(Zr)1+xO2 materials with the ultralow coercive field of Ec = 0.65 MV/cm. which aligns with our prediction our prediction because excess Hf(Zr) atom introduces more electrons than Ta atom. In addition, the energy difference between the f-phase and o-phase (ΔEf-o) decreases with an increasing hole concentration, indicating that the transition from the o-phase to the f-phase becomes easier. This is confirmed by XRD, which shows a phase transition from o-phase to f-phase with increasing Y doping concentration from 8% to 12% in bulk HfO2:Y single-crystalline12.

a Energy landscape in the process of structure transition from polarization-down to polarization-up of the undoped f-phase HfO2. Points marked with green dashed circles indicate the energies of the initial, transition, and final states, and purple arrows indicate the polarization switching barrier (ΔEB), the intermediate transition phase adopts the Pcca phase. Barrier (ΔEB) and energy difference (ΔEf-o = Ef - Eo) between the f-phase and o-phase as a function of (b) hole and (c) electron concentration. d Relative f-, m-, and c-phases fraction as a function of doping concentration, the La-doped HfO2 thin films as a reference [ref. 23]. e The spontaneous polarization magnitude of the f-phase HfO2 as a function of the carrier concentration.
The small energy differences between phases commonly result in multiphase coexistence in doped HfO2 thin films11,23. In Fig. 2d, the three (m, f, c) phase fractions under the hole concentration were calculated using the Boltzmann distribution, compared with La-doped HfO2 thin film23. The f-phase fraction showed first increasing from zero to maximum and then decreasing to zero, which is highly consistent with the experimental results. In addition, the calculated dielectric constant k of multiphase HfO2 the well agrees with the films (Supplementary Fig. 6a). In the case of idealized structures, the suitable concentration range tends to extend beyond what is observed in experimental conditions.
In Fig. 2e, the estimated spontaneous polarization Ps slightly increases (decreases) with the growth of hole (electron) concentration before reaching the critical concentration. Although quantifying the polarization value in metallic materials is uncertain55,56, the positive correlation between the remanent polarization Pr and the f-phase fraction23 suggests that Ps persists more or less unaffected by the carrier concentration. Since the f-phase fraction is not 100%, the Pr normally is much smaller than theoretically predicted. The polarization magnitude of multiphase HfO2 follows the same trend as the f-phase fraction (Supplementary Fig. 6b), consistent with the behavior of Pr in the La-doped HfO2 thin films23. Moreover, due to the presence of holes, the Pr of La-doped HZO films is larger than that of HZO films24.
Origin mechanism of hafnia ferroelectricity
To understand the microscopic mechanism of how carriers stabilize f-phase, we plotted the density of states (DOS) for the c-, f-, and m-phases with p = 0.2 h+/f.u. in Fig. 3a–c. In the c-phase, the holes are uniformly distributed in the OIV atom (Fig. 3a), while in the f-phase, the holes preferentially occupy OIII atoms, and then are slightly distributed in the OIV atom (Fig. 3b). Similar hole distribution patterns are also reflected in the m-, t-, and o-phase (Fig. 3c and Supplementary Fig. 7). Compared with the DOS of the m-phase and f (o) phase, the holes of the m-phase are more concentrated in the OIII atom. Such a hole distribution pattern does not change unless a phase structure transition occurs (Supplementary Fig. 8).

The atomic density of states (DOS) and 3D charge density distribution of the bulk HfO2 with hole concentration of 0.2 h+/f.u., (a) c-phase, (b) f-phase, (c) m-phase, the orange-yellow vertical dashed line represents the Fermi level (EF). The green isosurface demonstration of charge density distributions of the doping holes is drawn with the same scale in the inserted pictures. In-plane and out-of-plane orbital occupation number and proportion of f-phase HfO2 as a function of carrier concentration, (d), (f) hole doping, (e), (g) electron doping, the structural snapshots in the inset (g) show that the transformation of OIV to OIII at high electron concentrations with the exact opposite polar to neighboring OIII. (h) Negative crystal orbital Hamiltonian population (COHP) for Hf-O bond of f-phase HfO2 with different carrier concentrations. The COHP is determined by weighting the electronic DOS with the Hamiltonian elements corresponding to orbital overlaps, which divides DOS by energy into bonding, anti-bonding, and non-bonding regions. The convention is to plot negative COHP (-COHP, which is dimensionless). Hence bonding interactions and anti-bonding interactions are positive and negative, respectively. The number of electrons transferred (x-axis) and shared (y-axis) between adjacent basins intuitively illustrates electronic interactions of bonds in materials; they serve as quantitative measures of ionic and covalent properties, respectively. i A 2D map of different bonds in different materials and (j) Hf-O bonds with different carrier concentrations.
We mentioned earlier that symmetry is not the main factor that determines the total energy decrease slope. Judging from the perspective of the hole distribution pattern, it seems to indicate that the holes on the OIV atoms are the main factor for determining the total energy decrease slope in each phase. The total energy of the c-phase with more holes distributed on the OIV atoms declines rapidly, while few holes occupied the OIV atoms in the m-phase, making its slope the slowest (Supplementary Tab. 2). The most direct evidence comes from the m′-phase, which has lower symmetry than the c-phase, but its hole concentration on the OIV atoms is twice as high as that of the c-phase (Supplementary Tab. 2), so m′-phase has the steepest slope in Supplementary Fig. 1. Therefore, a ground-state phase transition from m → f (o)→ c (t)→ m′ was observed. For n-type doping, since electrons occupy all the Hf atoms in each phase (Supplementary Fig. 9), and increasing electron doping only increases the number of electrons on the Hf atoms without changing the electron distribution pattern (Supplementary Fig. 10), even if the f-phase undergoes a phase transition.
Gauss’ law stipulates that, in electrostatic equilibrium with effective charge screening, the net electric field within a conductor is zero30. Puggioni et al.57 attributed the coexistence of polarity and metallicity to the weak coupling between carriers near the Fermi level and the (soft) phonon(s) responsible for removing inversion symmetry. The number of holes (electrons) occupying the O 2p (Hf 5 d) orbitals is calculated by integrating the DOS near the Fermi level58. In Fig. 3d, it is evident that the holes predominantly occupy the in-plane px/py orbitals. Although the number of holes occupying the pz orbitals increases slowly, the proportion decreases all the time (Fig. 3f). The difficulty of in-plane holes to eliminate the out-of-plane dipoles generated by the non-centrosymmetric structure results in a weak coupling between the holes and the polarization distortion, thereby exhibiting strong ferroelectricity. However, it is important to acknowledge that, at excess concentrations, a considerable number of holes still occupy the pz orbital. This phenomenon explains why the Ps has to decrease after a slow increase, reaching its maximum until it eventually drops to zero at a critical concentration (Fig. 2e).
In Fig. 3e, the number of electrons occupied by the out-of-plane ({{rm{d}}}_{{rm{xz}}}{/{rm{d}}}_{{rm{yz}}}{/{rm{d}}}_{{{rm{z}}}^{2}}) orbitals steadily rises with increasing electron concentration, the heightened electron screening effect propels the metastable f-phase to transform into the antipolar f ″-phase at a critical concentration. Due to the inversion symmetry, the proportion of in-plane and out-of-plane orbital electrons remains unchanged after the phase transition (Fig. 3g). Excessive out-of-plane electrons intensify the electron screening effect, making it increasingly difficult to reverse f ″-phase.
Alongside electrostatic screening of carriers, chemical bonding arguments can also serve as a framework for understanding the intricate interplay between polar distortion and free carriers34. Notably, strong Ti-O covalent bonding59 has been identified as a significant factor in the lattice instability of ATiO3 (A = Ba, Pb). Moreover, robust Pb-O covalent bonding60 is considered to be the key factor that leads to the stronger ferroelectricity of PbTiO3 compared to BaTiO3. In Fig. 3h, undoped HfO2 being an insulator, positions the Fermi level in the Hf-O non-bonding region. Hole doping shifts the Fermi level into the Hf-OIII bonding region, fostering stronger Hf-OIII covalent bonding. However, concerning the Hf-OIV bond, the Fermi level occupies a lower position on the -COHP curve, residing in the very weak Hf-OIV anti-bonding region. This indicates that hole doping has minimal impact on the properties of the Hf-OIV bond, even when holes predominantly distributed on the OIV atom of t-phase (Supplementary Fig. 11). With electron doping, the Fermi level falls within both the Hf-OIII and Hf-OIV anti-bonding regions, but the Hf-OIV anti-bonding region is more extensive. This renders the Hf-OIV bond more unstable compared to the Hf-OIII bond. Consequently, during the f-phase transition, it is OIV that moves rather than OIII (Fig. 3g), resulting in the formation of the f ″-phase instead of the t-phase.
The transition between covalent and ionic bonds is gradual and lacks clear delineated. Recently developed quantum mechanical descriptions of shared electrons and transferred electrons have effectively addressed this issue, culminating in a 2D diagram for analyzing chemical bonds in solid materials61,62. Figure 3i offers an intuitive depiction, identifying the fundamental properties of ionic, metallic, and covalent bonds in a compound. Taking the example of binary oxides MgO, the relatively higher electrons transferred (ET) and few electrons shared (ES) in Mg-O bonds, underscore strong ionic bond properties. Whereas for Si, the ES in Si-Si bonds are overwhelmingly dominant, with virtually no ET, which signifies very strong covalent bond properties.
In Fig. 3i, the presence of more ES in the Hf-OIII bond in undoped HfO2, suggests that Hf-OIII covalent bonds are considerably stronger than Hf-OIV covalent bonds, a characteristic that favors hafnia ferroelectricity. Furthermore, Zr-OIII covalent bonds are stronger than Hf-OIII covalent bonds, potentially contributing to the enhancement of ferroelectricity with Zr doping. In Fig. 3j, under p-type doping, there is minimal change in ET and ES of Hf-OIV bonds. While the ES of Hf-OIII bonds in the f-phase experiences a slight decrease, the ET undergoes a significant reduction, resulting in an increased ES/ET ratio (Supplementary Fig. 12). Thus, the f-phase can still maintain strong Hf-OIII covalent bonds, a favorable condition for ferroelectricity. At low electron concentrations, the ES of Hf-OIV bond decreases not only faster than that of Hf-OIII bond in the f-phase but also faster than that of Hf-OIV bond of other phases, which is consistent with the anti-bonding states in Fig. 3h.
For real element doping, our previous studies have shown that II- and III-valent elements, which introduce holes, can stabilize the ferroelectric f-phase of hafnia, while V-valent elements that introduce electrons have the opposite effect63. IV-valent dopants, which do not introduce carriers, generally have minimal impact on f-phase stability, in line with findings by Künneth et al.64 Zr is the most widely studied dopant because Zr is structurally similar to Hf, allowing the formation of solid solutions with a wide range of Hf:Zr atomic ratios. However, from the perspective of phase energy, numerous theoretical studies21,63,64,65 of bulk HfO2 consistently shown that Zr doping has little effect on the phase energy difference (ΔEf-m/ΔEt-m) between the f-phase/t-phase and the m-phase, regardless of the Zr doping position or concentration. The origin of ferroelectricity in Zr-doped HfO2 films is complex, involving various factors, including oxygen vacancies15, surface energy16, annealing temperature, and grain size65, etc. Another commonly used IV-valence dopant, Si, induces ferroelectricity in HfO2 due to its small ionic radius that favors the t-phase42. As the Si doping concentration increases, the ΔEt-m decreases significantly, stabilizing the t-phase. As mentioned earlier, the t-phase at the crystallization temperature, allowing a phase transition to the f-phase during the cooling process by kinetic effects within the appropriate doping concentration range25.
To determine the key characteristic parameters that affect the stability of the doped HfO2 ferroelectric f-phase, we investigated 21 dopants spanning II, III, IV, and V-valent elements across the periodic table, as shown in Supplementary Fig. 13a. Six properties (Supplementary Tab. 3)—carrier concentration, ionic radius, atomic radius, electronegativity, ionization energy, and electron affinity as input variable sets (descriptors), while the energy difference ΔEf-m as the observed value, grouped by valence state, and principal component analysis (PCA) was performed66. The scree plot (Supplementary Fig. 13b) shows that the first principal component (PC1) can explain 70.8% of the variance in the original descriptor, which means that PC1 can capture most of the changes in the data. The inclusion of the second principal component (PC2) can explain 87.5% of the variance in the descriptor data, indicating that they have captured the main trend of most variables.
In the loading plot (Supplementary Fig. 13c), PC1 and PC2 are primarily influenced by carrier concentration and ionization energy. In PC1, the most significant component, carrier concentration, is positively correlated, while ionization energy is negatively correlated. Since hole concentration is designated as positive in our carrier data, a higher hole concentration correlates with enhanced stability of the ferroelectric phase, while lower ionization energy also promotes this stability. The other four variables contribute less, with atomic and ionic radii positively correlated and electronegativity and electron affinity also showing a positive correlation with ionization energy.
Combining the load plot (Supplementary Fig. 13c) and the score plot (Supplementary Fig. 13d), we observe that the major axis of the IV-valent element ellipse aligns with the direction of ionization energy, indicating that IV-valent elements are strongly influenced by ionization energy. Dopants with lower ionization energy tend to stabilize the ferroelectric phase. For II and III-valent elements, their sample points cluster along the positive PC1 axis with some overlap, reflecting their similar behavior. The carrier load value in the positive PC1 axis is the highest, showing that the phase energy difference under II- and III-valent doping is primarily influenced by carriers. Moreover, the relatively scattered distribution of these scores suggests that ionic radius also contributes to the phase energy difference. Furthermore, because II-valent elements introduce more holes than III-valent elements, II-valent elements score higher than III-valent elements. Since electron concentration is marked as negative, V-valent element samples are mostly located in the third quadrant, with their ellipse’s major axis opposite to the carrier direction in the load plot, suggesting that V-valent elements are more sensitive to carrier concentration.
Interface carrier effect in hafnia-based films
Epitaxial atomic-scale interface with polar catastrophes creates unique phenomena, such as the formation of two-dimensional electron/hole gas (2DEG/2DHG)67,68,69. Shi et al.37 recently unveiled that hole doping in the HZO layer, resulting from the MnO2 interface termination in HZO/La0.67Sr0.33MnO3 thin film, stably reduces the clear ferroelectric features to a 1.5 nm thick HZO. Additionally, Cheema et al.39 demonstrate ferroelectricity in a 1 nm HZO ultrathin film grown by low-temperature atomic layer deposition (ALD) on Si/SiO2. Subsequently, Cheema et al.70,71 also found that the thickness of ferroelectric ZrO2 films and HfO2-ZrO2 multilayer could be maintained at 0.5 nm and 2 nm, indicating an inverse size effect for hafnia-based films. Its structural characterization shows that only metal-oxygen dangling bonds are present out-of-plane. But as we know, these metal-oxygen dangling bonds inevitably bring surface or interfacial holes. As we mentioned earlier, holes can enhance the ferroelectricity of bulk hafnia, so here we further explore the effect of holes on the ferroelectric properties of f-phase HfO2 thin films. In Supplementary Fig. 14, constructing several models of SiO2/f-HfO2 and TiN/f-HfO2 heterostructure (HS) with oxygen terminations, capping with or without TiN electrode.
In Supplementary Fig. 14, it is evident that only two O atoms are paired with Si at the SiO2/f-HfO2 interface, while the remaining two dangling O atoms will form a metallic state at the SiO2/f-HfO2 interface, as shown in Fig. 4. Undoubtedly, the most dominant holes in the heterostructures come from the dangling O atoms on the surface. Figure 4a clearly demonstrates that the SiO2/f-HfO2 HS, without being capped with TiN, can still maintain the polar structure at a thickness of 1.5 nm. Essentially, every O layer in f-HfO2 exhibits hole distribution, with a higher concentration on the polar OIII atoms, promoting the formation of the f-phase in the HfO2 thin film. Since the number of holes generated at the interface and the surface remains constant, thinner films have a higher hole concentration and greater polar distortion (Supplementary Fig. 15). In Fig. 5, we observe that the aspect ratio of HfO2 thin films shows nearly exponential growth with decreasing thickness, similar to the trend observed in HZO ultrathin films39. Additionally, with the same hole doping, a 3 nm BaTiO3 film exhibits antipolar structure, while a 3 nm HfO2 film maintains a strong polar structure (Supplementary Fig. 16). Shi et al. introduced holes into HfO2-based films through the work function difference at the metal/dielectric layer interface, whereas we utilized dangling oxygen at the interface, indicating that the existence of holes may be one of the reasons for the lack of critical thickness of HfO2-based films. Other methods of introducing holes at an interface also include using the mismatch of interface chemical valence to affect the electrostatic potential at the interface68.

Calculated layer-resolved DOS along with its 3D charge density projected on the bands forming the metallic states near the Fermi level for (a) SiO2/f-HfO2, (b) SiO2/f-HfO2/TiN, (c) TiN/f-HfO2 HS models. Since the coordination number of some O atoms cannot be determined after relaxation, the different types of O atoms are no longer distinguished by color in the HfO2-based HS. The holes and electrons are depicted by green and yellow isosurfaces, respectively. The blue arrows indicate the direction of polarization, which is directed along the line from the negative charge toward the positive charge. The cyan line represents the electrostatic potential. The purple box in (b) shows a thickness of 1 nm.

Aspect ratio 2a/(b + c) as a function of HfO2 film thickness, with experimental HZO film as comparison39. Dashed line represents the aspect ratio of bulk f-phase HfO2. The inset shows the 3D hole distribution in the 1 nm HfO2 ultra-film.
Because it is not feasible to measure ferroelectricity in metals, the annealed HZO region without metal capping does not exhibit ferroelectric signatures. However, clear PFM phase contrast images emerge on the HZO region with capped metal during phase annealing, and 180° phase contrast and butterfly-shaped amplitude are observed only for metal-capped HZO39. This phenomenon is elucidated by the layer-resolved partial DOS in SiO2/f-HfO2/TiN HS (Fig. 4b). Besides forming a metallic state on SiO2/f-HfO2 (p-type) and f-HfO2/TiN (n-type) interface layer, the interior of the HfO2 thin film remains insulating. Due to the asymmetry interface, the robust electrostatic potential (cyan line in Fig. 4b) propels the polarization direction switching of the initially constructed structure. Moreover, Cheema et al.39 showed PFM phase contrast and topographical imaging for 1 nm HZO on TiN-buffered Si without metal capping. Figure 4c illustrates a small number of holes on the surface of the HfO2 film in the TiN/f-HfO2 HS, while the film’s interior retains insulation, affirming that the 1 nm HZO on TiN exhibits ferroelectric signatures. Capping with a TiN electrode also leads to n-type conductivity at the top interface, with minimal alterations in the f-phase structure morphology (Supplementary Fig. 17). Moreover, La doping introduces more holes into SiO2/f-HfO2@La HS (Supplementary Fig. 18a), the high hole concentration resulting from the superposition of doped holes and interface holes, may explain the formation of a very thin interfacial c/t-phase layer in the La-doped23 and Gd-doped72 thin films. Once capped with TiN, almost all holes are absorbed (Supplementary Fig. 18b), and even electrons brought by Ta doping are absorbed (Supplementary Fig. 18c), ensuring insulation.
In general, as the TiN top electrode is typically sputter-deposited on already-existing HfO2-based thin films, the top electrode cannot affect the growth behavior of dielectric layers. We speculate that these carriers play a pivotal role in initiating the formation of the f-phase in the uncapped HfO2-based films. Once this task is accomplished, in the presence of capped electrodes, these carriers are assimilated by the capping electrodes to maintain charge neutrality within the films. This could explain why, in most instances, capped HfO2 layers possess better ferroelectricity than those without capping5. In addition, this phenomenon is supported by Pt/Ta-doped HfO2/Pt capacitors27. Even with significant Ta doping concentrations (Ta atomic ratio of 11%–21%), strong ferroelectricity can still be measured. Finally and most importantly, due to the limited screening length of the electrodes on charge carriers, hafnia films with strong ferroelectricity can be found to have thicknesses of essentially ~10 nm or less, much smaller than that of perovskite ferroelectric materials. Of course, apart from the compensation with surface charges, the more or less existing oxygen vacancies in the doped HfO2-based thin film play this role because they can locally neutralize excess charges by generating oxygen vacancy-dopant pairs8. In addition, co-doping of n-p type elements can not only achieve charge neutrality but also enhance the ferroelectric properties of hafnia-based films57,73.
In summary, our first-principles calculations unveil distinct mechanisms underlying the induction of ferroelectricity in hafnia when doped with p-type and n-type dopants. For p-type doping, the holes concentrated on the polar OIII atoms in the f-phase significantly enhance robust ferroelectricity, driven by a weak hole screening effect and augmented polar Hf-OIII covalent bonds. However, the presence of holes on nonpolar OIV atoms crucially influences the phase total energy change, leading to a ground state m → f → c phase transition as hole concentration increases. As a result, the f-phase fraction and polarization magnitude exhibit an initial increase followed by a decrease. In contrast, n-type doping involves strong electron screening and a weaker Hf-OIV covalent bond, inducing a transition from the polar metastable f-phase to the antipolar stable f ″-phase. This distinction might explain the challenges in identifying ferroelectricity in n-type doped hafnia. Lastly, we emphasize the pivotal role of electrode capping in ensuring the ferroelectric insulation of hafnia-based films.
While the internal energy calculated by the standard density-functional theory (DFT) approach has certain limitations in fully capturing the complexity of ferroelectric f-phase stability in hafnia—particularly under the influence of temperature, grain size, and domain structure—our findings are broadly consistent with experimental measurements of doped hafnia films. This alignment underscores the potential of precisely modulating carrier concentration as a promising strategy to significantly enhance the ferroelectric properties of hafnia-based devices.
To further improve the precision and depth of such studies, overcoming current theoretical and computational limitations will be essential. Free energy differences, being more reliable than internal energies alone in predicting ferroelectric phase stability, highlight the potential of advanced computational approaches—such as ab initio molecular dynamics (MD), stochastic self-consistent harmonic approximations (SSCHA), and temperature-dependent effective potential (TDEP) methods. These methods have shown promise in addressing critical factors such as temperature effects and anharmonicity but remain computationally prohibitive for systematically exploring multiple dopants and phases.
Looking ahead, the development of machine-learning-based force-field models presents a compelling opportunity to overcome these challenges. By integrating machine learning with traditional quantum mechanical frameworks, it may become possible to efficiently explore phase stability in doped systems and uncover the intricate mechanisms by which dopants modulate ferroelectric properties. Such advancements could pave the way for more precise theoretical insights and optimized designs of ferroelectric materials and devices.
Methods
First-principles calculations of bulk structures
We consider five competing phases, including ground nonpolar monoclinic m-phase, cubic c-phases, tetragonal t-phases, polar orthorhombic f-phase, and antipolar orthorhombic o-phase, in this work. Our work is based on DFT as implemented in the Vienna ab initio simulation package (VASP)74,75, using a plane-wave basis set with a kinetic energy of 500 eV, employing the standard projector-augmented wave (PAW) methodology and the generalized gradient approximation (GGA)76 parametrized by the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. For structure relaxations of each phase hafnia, a 7 × 7 × 7 Gamma-centered k-point meshes are used, all lattice parameters and atomic positions were allowed to relax until the energy differences were within 10−6 eV and the Hellman-Feynman forces were less than 0.01 eV/Å. For the DOS calculations, denser k-point meshes with an 11 × 11 × 11 Gamma-centered grid are used. Supplementary Tab. 1 lists the lattice parameters of the five phases of hafnia we used in this work together with previously reported theoretical and experimental results, as predicted by our DFT calculations, which are in good agreement with the experimental data.
We use the smallest cells describing the equilibrium structures of the five phases of hafnia, simulate the effect of doping by varying the number of electrons in the cell, and add a neutralizing homogeneous charge background. Therefore, the doping holes were introduced to the system with a negative background to achieve charge neutrality, while doping electrons did the opposite. This approach does not take into account the fundamental properties of dopant atoms, such as ionic radius, electronegativity, etc., which greatly simplifies calculations and provides insight into the origin of ferroelectricity in p-type and n-type doped hafnia. In addition, crystal orbital Hamilton population (COHP) analyses are carried out with the LOBSTER code77. Moreover, the energy barriers for polarization switching are obtained with the aid of the climbing image nudged elastic band (NEB) method. Crystal structures and the isosurfaces of the three-dimensional (3D) charge density are visualized using the VESTA package78.
The polarization magnitude of f-phase hafnia was calculated by Born effective charges (BEC, Z*)79 multiplied by atomic displacements (δd) with respect to the referenced centrosymmetric t-phase. Formally, Ps simply took the form of:
Where i represents the ith atom, e and Ω are the charge of the electron and the cell volume considered, respectively. The BEC along the c-axis derived from density-functional perturbation theory calculations80. Our spontaneous polarization Ps = 58.41 μC/cm2 of undoped f-phase is close to the Ps = 53.31 μC/cm2 predicted by Zhou et al.47
Chemical bond calculations
For analyzing electronic interactions and bonding in materials, the initial DFT wavefunctions were post-processed in Quantum Espresso81, then followed by a transformation into maximally localized Wannier functions (MLWFs)82 and integration within the Critic2 code83.
The integral of the nonclassical exchange-(correlation) part of the electron-pair density function satisfies the summation rule84:
Where the localization index LI(Ω), measuring the number of electrons that are fully localized in the Bader’s basin, depending on the atom type and its specific chemical environment. while the delocalization index δ(Ω′,Ω), measures the number of electron pairs being shared between the two basins, therefore, their double corresponds to the amount of shared electrons (ES). N(Ω) is the electron population of the Ω atom, and the difference between it and the total number of electrons of Ω atom is expressed as the transferred electron (ET) between two atomic basins.
Heterostructure modeling and calculation
In addition, we constructed SiO2/HfO2 HS and SiO2/HfO2/TiN HS to explore the influence of interface carriers on the f-phase structure and ferroelectricity. Among them, SiO2 adopts the most stable crystal structure with space group I(bar{4}) 2d (a = b = 4.95 Å, c = 7.31 Å, α = β = γ = 90°) found in the Materials Project web85, which has a relatively small in-plane lattice mismatch with the f-phase structure. In 2020, Cheema et al.39 prepared ferroelectric Si/SiO2/HZO/TiN thin films and used out-of-plane one-dimensional GI-XRD patterns and two-dimensional pole figure patterns to reflect the stacking of the (311) and (111) densely packed planes in the fluorite structure. The results confirmed that all the cation sites lie in-planes, which minimizes surface energy effects because only metal-oxygen dangling bonds are present out-of-plane. So here, both the surface and the interface of the f-phase film in our heterojunction are O-terminated. The TiN adopts the crystal structure of space group Fm(bar{3}) m (a = b = c = 4.22 Å), to reduce the in-plane lattice mismatch between the TiN electrode and the f-phase hafnia (a = 5.27 Å, b = 5.05 Å) film, we constructed the TiN electrode along the (111) direction (a = 5.97 Å, b = 5.17 Å) compared to the (001) direction (a = b = 4.22 Å). In addition, TiN in the (111) direction is alternately stacked with Ti terminals and N terminals, which makes it easier to adapt to the O terminals of the f-phase at the interface. We construct SiO2/f-HfO2, SiO2/f-HfO2/TiN(111), and TiN(111)/f-phase HfO2 HS along the (001) direction with a 20 Å vacuum layer, and fix the in-plane lattice constants of all constructed HS to those of the f-phase. The 6×6×1 k-mesh and Hellman-Feynman forces were less than 0.05 eV/Å are used in VASP calculation.
Responses