Device-area selection of memristive transport regimes in epitaxial Hf_(0.5)Zr_(0.5)O₂-based ferroelectric devices
Pith reviewed 2026-05-10 10:18 UTC · model grok-4.3
The pith
Small epitaxial hafnia ferroelectric devices exhibit area-distributed tunneling transport while larger devices form localized conductive channels, with a crossover at roughly 1000 square micrometers that also marks the onset of wake-up.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In epitaxial Hf0.5Zr0.5O2/La0.67Sr0.33MnO3 heterostructures with Pt electrodes, devices smaller than a few thousand square micrometers display low-resistance states whose resistance scales inversely with area, indicating area-distributed tunneling, whereas larger devices exhibit area-independent resistance from localized conductive channels. A statistical nucleation model reproduces the observed resistance-area dependence and identifies a crossover area A* ≈ 10^3 μm² that also correlates with the onset of ferroelectric wake-up, providing a unified picture in which channel nucleation and oxygen-vacancy redistribution compete with polarization-driven changes.
What carries the argument
The statistical nucleation model for localized conductive channels, which treats the probability of forming one or more area-independent paths as a function of device area and thereby predicts the transition from inverse scaling to constant resistance.
If this is right
- Lateral device area can be used as a design parameter to select between distributed tunneling and filamentary transport regimes in hafnia memristors.
- The same nucleation process that sets the resistance crossover also controls when ferroelectric wake-up appears.
- Oxygen-vacancy redistribution and conductive-channel formation occur together once the device area allows stable nucleation.
- Neuromorphic circuits can exploit this size dependence to achieve either uniform or stochastic switching by choosing appropriate lateral dimensions.
Where Pith is reading between the lines
- For scaled-down neuromorphic arrays, devices kept below the crossover area should give more reproducible, area-scalable conductance states.
- Variability in large devices will be dominated by the random locations of the few nucleated channels rather than by film thickness fluctuations.
- Changing electrode materials or oxygen-vacancy engineering could shift the crossover area in a predictable way, allowing the same film to be used for either regime.
- The model implies that cycling-induced wake-up is not purely ferroelectric but partly the result of the same vacancy motion that creates the channels.
Load-bearing premise
That the observed inverse scaling of resistance with area is caused solely by uniform tunneling and that constant resistance with area is caused solely by localized channels, without significant size-dependent changes in interface quality or defect density.
What would settle it
Fabricating a new set of devices with areas both well below and well above 10^3 μm² and finding either that the resistance-area product stops scaling inversely below the crossover or that the wake-up onset does not shift with the same area threshold.
Figures
read the original abstract
Ferroelectric memristive devices based on hafnia are promising systems for neuromorphic electronics, yet the interplay between polarization-modulated resistive changes and defect-mediated transport often leads to complex and debated switching mechanisms. Here, we investigate this competition in epitaxial Hf$_{0.5}$Zr$_{0.5}$O$_2$/La$_{0.67}$Sr$_{0.33}$MnO$_3$ heterostructures with Pt top electrodes by combining structural, ferroelectric, and memristive characterization with a statistical analysis across a broad range of device areas spanning three orders of magnitude. We identify two distinct memristive regimes with opposite resistance--voltage chiralities. Small devices exhibit a low-resistance state that scales inversely with area, consistent with area-distributed tunneling transport, while larger devices display an area-independent resistance indicative of localized conductive channels. A statistical nucleation model quantitatively captures this behavior and yields a crossover characteristic area $A^* \approx 10^3~\mu\mathrm{m}^2$. This crossover also correlates with the onset of ferroelectric wake-up in larger devices, linking conductive-channel nucleation and oxygen-vacancy redistribution within a unified physical picture. These results establish lateral device size as a key parameter controlling the dominant transport mechanism in epitaxial hafnia-based devices.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates memristive transport regimes in epitaxial Hf_{0.5}Zr_{0.5}O_2/La_{0.67}Sr_{0.33}MnO_3 heterostructures with Pt top electrodes. By performing structural, ferroelectric, and electrical characterization over device areas spanning three orders of magnitude, the authors identify two regimes with opposite resistance-voltage chiralities: small devices exhibit low-resistance states scaling inversely with area (attributed to distributed tunneling), while larger devices show area-independent resistance (attributed to localized conductive channels). A statistical nucleation model is reported to quantitatively capture the crossover at a characteristic area A^* ≈ 10^3 μm², which also correlates with the onset of ferroelectric wake-up in larger devices.
Significance. If the supporting data and model validation hold, the work offers a useful framework for understanding the competition between polarization-modulated and defect-mediated transport in hafnia-based ferroelectric memristors. The broad device-area range and the proposed link between nucleation, oxygen-vacancy redistribution, and wake-up effects could inform geometry-dependent design rules for neuromorphic applications. The multi-technique approach is a clear strength.
major comments (3)
- [Abstract] Abstract: The claim that the statistical nucleation model 'quantitatively captures' the crossover behavior requires explicit presentation of the model equations, the number of devices measured per area bin, resistance-area data with error bars, fit quality metrics (e.g., R² or χ²), and device exclusion criteria. Without these, it is not possible to assess whether the model is predictive or partly fitted to the same dataset from which A^* is extracted.
- [Abstract] Abstract / Results: The interpretation that inverse resistance-area scaling unambiguously signals area-distributed tunneling transport (while area-independent resistance signals localized channels) does not address possible confounding size-dependent effects such as variations in electrode interface quality, lithography-induced edge defects, or processing-dependent oxygen-vacancy distributions. These alternatives must be ruled out with additional controls or simulations for the central claim to be load-bearing.
- [Abstract] Abstract: The reported correlation between the crossover area A^* and the onset of ferroelectric wake-up is presented as evidence for a unified physical picture, but the manuscript should clarify whether the wake-up data are independent of the resistance-area dataset used to determine A^* and provide a quantitative statistical measure of the correlation strength.
minor comments (1)
- [Abstract] Abstract: Notation for the crossover area should be standardized (A^* vs. A*) and the units (μm²) consistently applied throughout.
Simulated Author's Rebuttal
We thank the referee for the thorough and constructive review of our manuscript. We address each of the major comments point by point below, indicating the revisions we will make to strengthen the paper.
read point-by-point responses
-
Referee: [Abstract] Abstract: The claim that the statistical nucleation model 'quantitatively captures' the crossover behavior requires explicit presentation of the model equations, the number of devices measured per area bin, resistance-area data with error bars, fit quality metrics (e.g., R² or χ²), and device exclusion criteria. Without these, it is not possible to assess whether the model is predictive or partly fitted to the same dataset from which A^* is extracted.
Authors: We agree with the referee that these details are essential for a rigorous assessment of the model. In the revised manuscript, we will add the explicit equations of the statistical nucleation model in the methods or results section. We will also include a table summarizing the number of devices measured per area bin, plot the resistance-area data with error bars, report fit quality metrics such as R² and χ², and clearly state the device exclusion criteria. These additions will demonstrate that the model is not merely fitted to the data but provides a predictive framework. revision: yes
-
Referee: [Abstract] Abstract / Results: The interpretation that inverse resistance-area scaling unambiguously signals area-distributed tunneling transport (while area-independent resistance signals localized channels) does not address possible confounding size-dependent effects such as variations in electrode interface quality, lithography-induced edge defects, or processing-dependent oxygen-vacancy distributions. These alternatives must be ruled out with additional controls or simulations for the central claim to be load-bearing.
Authors: This is a valid concern, as size-dependent fabrication effects could influence transport. The manuscript relies on the clear scaling behaviors observed over three orders of magnitude in area, combined with consistent ferroelectric characterization, to support the distributed tunneling vs. localized filament interpretation. To strengthen this, we will expand the discussion to explicitly consider and argue against the listed confounding factors, using perimeter-to-area ratio analysis to estimate edge defect contributions and referencing literature on interface quality in LSMO/HZO systems. While we cannot perform new simulations or experiments in this revision, the existing data across multiple techniques provide supporting evidence. We believe this addresses the comment without overclaiming. revision: partial
-
Referee: [Abstract] Abstract: The reported correlation between the crossover area A^* and the onset of ferroelectric wake-up is presented as evidence for a unified physical picture, but the manuscript should clarify whether the wake-up data are independent of the resistance-area dataset used to determine A^* and provide a quantitative statistical measure of the correlation strength.
Authors: The wake-up measurements were performed on an independent set of devices not included in the resistance-area dataset used to extract A^*. In the revised version, we will explicitly note this independence and include a quantitative statistical analysis of the correlation, such as the Spearman rank correlation coefficient between the crossover area and the wake-up onset, along with the associated p-value. This will provide a more rigorous basis for the unified physical picture linking nucleation, oxygen-vacancy redistribution, and wake-up effects. revision: yes
Circularity Check
No significant circularity in the derivation chain
full rationale
The paper reports direct experimental measurements of resistance versus device area over three orders of magnitude, identifies the change from inverse-area to area-independent scaling, and introduces a statistical nucleation model whose role is to describe the observed crossover at A* ≈ 10^3 μm². This is standard phenomenological modeling fitted to the same dataset rather than an independent first-principles prediction that reduces to the inputs by construction. No self-citations, uniqueness theorems, or ansatzes are invoked to justify the core claims; the wake-up correlation is presented as an additional empirical observation. The derivation chain therefore remains self-contained as empirical characterization plus descriptive modeling.
Axiom & Free-Parameter Ledger
free parameters (1)
- crossover area A* =
~10^3 μm²
axioms (2)
- domain assumption Low-resistance state scaling inversely with area indicates area-distributed tunneling transport
- domain assumption Area-independent resistance indicates localized conductive channels
Reference graph
Works this paper leans on
-
[1]
A. Mehonic and A. Kenyon, Brain-inspired computing needs a master plan, Nature604, 255 (2022)
work page 2022
-
[2]
A. Mehonic, D. Ielmini, K. Roy, O. Mutlu, S. Kvatinsky, T. Serrano-Gotarredona, B. Linares- Barranco, S. Spiga, S. Savel’ev, A. G. Balanov, N. Chawla, G. Desoli, G. Malavena, C. Monzio Compagnoni, Z. Wang, J. J. Yang, S. G. Sarwat, A. Sebastian, T. Mikolajick, B. Noheda, S. Slesazeck, B. Dieny, T.-H. Hou, A. Varri, F. Bruckerhoff-Pl¨ uckelmann, W. Per- ni...
work page 2024
-
[3]
M. Zolfagharinejad, U. Alegre-Ibarra, T. Chen, S. Kinge, and W. G. van der Wiel, Brain- inspired computing systems: a systematic literature review, The European Physical Journal B97, 70 (2024)
work page 2024
-
[4]
F. L. Aguirre, A. Sebastian, M. Le Gallo, W. Song, T. Wang, J. J. Yang, W. Lu, M.-F. Chen, D. Ielmini, Y. Yang, A. Mehonic, A. J. Kenyon, I. Valov, R. Waser, M. Lanza,et al., Hardware implementation of memristor-based artificial neural networks, Nature Communications15, 1974 (2024)
work page 1974
-
[5]
S. Schmidgall, R. Ziaei, J. Achterberg, L. Kirsch, S. P. Hajiseyedrazi, and J. Eshraghian, Brain-inspired learning in artificial neural networks: A review, APL Machine Learning2, 021501 (2024)
work page 2024
-
[6]
Sawa, Resistive switching in transition metal oxides, Materials Today11, 28 (2008)
A. Sawa, Resistive switching in transition metal oxides, Materials Today11, 28 (2008)
work page 2008
-
[7]
S. Yu,Neuro-inspiring computing using resistive synaptic devices(Springer International Pub- lishing, 2017). 13
work page 2017
-
[8]
M. J. Rozenberg, M. J. S´ anchez, R. Weht, C. Acha, F. Gomez-Marlasca, and P. Levy, Mecha- nism for bipolar resistive switching in transition-metal oxides, Phys. Rev. B81, 115101 (2010)
work page 2010
-
[9]
A. Herpers, C. Lenser, C. Park, F. Offi, F. Borgatti, G. Panaccione, S. Menzel, R. Waser, and R. Dittmann, Spectroscopic proof of the correlation between redox-state and charge-carrier transport at the interface of resistively switching Ti/PCMO devices, Advanced Materials26, 2730 (2014)
work page 2014
-
[10]
R. Degraeve, A. Fantini, N. Raghavan, L. Goux, S. Clima, B. Govoreanu, A. Belmonte, D. Linten, and M. Jurczak, Causes and consequences of the stochastic aspect of filamentary RRAM, Microelectronic Engineering147, 171 (2015)
work page 2015
-
[11]
J. L. Rieck, F. V. E. Hensling, and R. Dittmann, Trade-off between variability and retention of memristive epitaxial SrTiO 3 devices, APL Materials9, 021110 (2021)
work page 2021
-
[12]
E. Y. Tsymbal and H. Kohlstedt, Tunneling across a ferroelectric, Science313, 181 (2006)
work page 2006
- [13]
-
[14]
M. Y. Zhuravlev, Y. Wang, S. Maekawa, and E. Y. Tsymbal, Tunneling electroresistance in ferroelectric tunnel junctions with a composite barrier, Appl. Phys. Lett.95, 052902 (2009)
work page 2009
-
[15]
A. Chanthbouala, V. Garcia, R. O. Cherifi, K. Bouzehouane, S. Fusil, X. Moya, S. Xavier, H. Yamada, C. Deranlot, N. D. Mathur,et al., A ferroelectric memristor, Nature Mater.11, 860 (2012)
work page 2012
-
[16]
P. W. M. Blom, R. M. Wolf, J. F. M. Cillessen, and M. P. C. M. Krijn, Ferroelectric Schottky diode, Phys. Rev. Lett.73, 2107 (1994)
work page 1994
-
[17]
R. Meyer and R. Waser, Hysteretic resistance concepts in ferroelectric thin films, J. Appl. Phys.100, 051611 (2006)
work page 2006
-
[18]
L. Pintilie, V. Stancu, L. Trupina, and I. Pintilie, Ferroelectric Schottky diode behavior from a SrRuO3-Pb(Zr0.2Ti0.8)O3-Ta structure, Phys. Rev. B82, 085319 (2010)
work page 2010
-
[19]
C. Ferreyra, M. Rengifo, M. S´ anchez, A. Everhardt, B. Noheda, and D. Rubi, Key role of oxygen-vacancy electromigration in the memristive response of ferroelectric devices, Phys. Rev. Applied14, 044045 (2020)
work page 2020
-
[20]
M. Rengifo, M. H. Aguirre, M. Sirena, U. L¨ uders, and D. Rubi, Epitaxial ferroelectric mem- ristors integrated with silicon, Frontiers in Nanotechnology4(2022). 14
work page 2022
-
[21]
T. S. B¨ oscke, J. M¨ uller, D. Br¨ auhaus, U. Schr¨ oder, and U. B¨ ottger, Ferroelectricity in hafnium oxide thin films, Applied Physics Letters99, 102903 (2011)
work page 2011
-
[22]
U. Schroeder, C. S. Hwang, and H. Funakubo, eds.,Ferroelectricity in Doped Hafnium Oxide: Materials, Properties and Devices, Woodhead Publishing Series in Electronic and Optical Materials (Woodhead Publishing, Cambridge, UK, 2019)
work page 2019
- [23]
- [24]
-
[25]
P. Nukala, M. Ahmadi, Y. Wei, S. de Graaf, E. Stylianidis, T. Chakrabortty, S. Matzen, H. W. Zandbergen, A. Bj¨ orling, D. Mannix, D. Carbone, B. Kooi, and B. Noheda, Reversible oxygen migration and phase transitions in hafnia-based ferroelectric devices, Science372, 630 (2021)
work page 2021
-
[26]
M. C. Sulzbach, S. Estand´ ıa, X. Long, J. Lyu, N. Dix, J. G` azquez, M. F. Chisholm, F. S´ anchez, I. Fina, and J. Fontcuberta, Unraveling ferroelectric polarization and ionic contributions to electroresistance in epitaxial Hf 0.5Zr0.5O2 tunnel junctions, Adv. Electr. Mater. , 1900852 (2019)
work page 2019
-
[27]
J. Lee, K. Yang, J. Y. Kwon, J. E. Kim, D. I. Han, D. H. Lee, J. H. Yoon, and M. H. Park, Role of oxygen vacancies in ferroelectric or resistive switching hafnium oxide, Nano Convergence 10, 55 (2023)
work page 2023
-
[28]
K. Liu, K. Liu, X. Zhang, F. Jin, J. Fang, E. Hua, H. Ye, J. Zhang, Z. Liang, Q. Lv, W. Wu, C. Ma, and L. Wang, Optimizing the ferroelectric performance of Hf 0.5Zr0.5O2 epi- taxial film by La 0.67Sr0.33MnO3 capping layer, Advanced Electronic Materials10, 2400136 (2024), https://advanced.onlinelibrary.wiley.com/doi/pdf/10.1002/aelm.202400136
-
[29]
A. Jan, T. Rembert, S. Taper, J. Symonowicz, N. Strkalj, T. Moon, Y. S. Lee, H. Bae, H. J. Lee, D.-H. Choe, J. Heo, J. MacManus-Driscoll, B. Monserrat, and G. Di Martino, In operando optical tracking of oxygen vacancy migration and phase change in few-nanometer ferroelectric HZO memories, Advanced Functional Materials33, 2214970 (2023)
work page 2023
-
[30]
Y. Wei, P. Nukala, M. Salverda, S. Matzen, H. J. Zhao, J. Momand, A. S. Everhardt, G. Agnus, G. R. Blake, P. Lecoeur, B. J. Kooi, J. ´I˜ niguez, B. Dkhil, and B. Noheda, A rhombohedral ferroelectric phase in epitaxially strained Hf 0.5Zr0.5O2 thin films, Nature Mater.17, 1095 (2018). 15
work page 2018
-
[31]
J. Knabe, F. Berg, K. Thorben Goss, U. Boettger, and R. Dittmann, “Dual-mode operation of epitaxial Hf 0.5Zr0.5O2: Ferroelectric and filamentary-type resistive switching,”physica status solidi (a)230, 2300409 (2023)
work page 2023
-
[32]
X. Long, H. Tan, F. S´ anchez, I. Fina, and J. Fontcuberta, Ferroelectric electroresistance after a breakdown in epitaxial Hf 0.5Zr0.5O2 tunnel junctions, ACS Applied Electronic Materials5, 740 (2023)
work page 2023
-
[33]
A. Petraru, O. Gronenberg, U. Sch¨ urmann, L. Kienle, R. Droopad, and H. Kohlst- edt, Distinguishing the rhombohedral phase from orthorhombic phases in epitaxial doped HfO 2 ferroelectric films, ACS Applied Materials & Interfaces16, 42534 (2024), https://doi.org/10.1021/acsami.4c10423
-
[34]
A. El Boutaybi, T. Maroutian, L. Largeau, S. Matzen, and P. Lecoeur, Stabilization of the epi- taxial rhombohedral ferroelectric phase in ZrO 2 by surface energy, Physical Review Materials 6, 074406 (2022)
work page 2022
-
[35]
X. Lyu, J. Zou, F. Ali, I. Fina, and F. S´ anchez, Intrinsic ferroelectric polarization in doped HfO2 and ZrO2, Small21, e2506334 (2025)
work page 2025
-
[36]
D. Rubi, F. Tesler, I. Alposta, A. Kalstein, N. Ghenzi, F. Gomez-Marlasca, M. Rozenberg, and P. Levy, Two resistive switching regimes in thin film manganite memory devices on silicon, Appl. Phys. Lett.103, 163506 (2013)
work page 2013
-
[37]
C. Ferreyra, M. S´ anchez, M. Aguirre, C. Acha, S. Bengi´ o, J. Lecourt, U. L¨ uders, and D. Rubi, Selective activation of memristive interfaces in TaO x-based devices by controlling oxygen va- cancies dynamics at the nanoscale, Nanotechnology31, 155204 (2020)
work page 2020
-
[38]
W. Rom´ an Acevedo, C. Ferreyra, M. J. S´ anchez, C. Acha, R. Gay, and D. Rubi, Concurrent ionic migration and electronic effects at the memristive TiO x/La1/3Ca2/3MnO3−x interface, J. Phys. D: Appl. Phys.51, 125304 (2018)
work page 2018
-
[39]
Y. Wei, S. Matzen, C. P. Quinteros, T. Maroutian, G. Agnus, P. Lecoeur, and B. Noheda, Magneto-ionic control of spin polarization in multiferroic tunnel junctions, npj Quantum Mater.4, 62 (2019)
work page 2019
-
[40]
S. Monaghan, P. K. Hurley, K. Cherkaoui, M. A. Negara, and A. Schenk, Determination of electron effective mass and electron affinity in HfO 2 using MOS and MOSFET structures, Solid-State Electronics53, 438 (2009)
work page 2009
-
[41]
S. M. Ross,Introduction to Probability Models, 11th ed. (Academic Press, 2014). 16
work page 2014
- [42]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.