Ferroelectric dynamic-field-driven nucleation and growth model for predictive materials-to-circuit co-design
Pith reviewed 2026-05-22 11:31 UTC · model grok-4.3
The pith
A compact model fits ferroelectric switching data under arbitrary time-varying voltages to extract domain wall velocity and growth parameters for direct use in circuit simulations.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that a compact dynamic-field-driven nucleation and growth model enables quantitative fits to switching transients across multiple ferroelectric materials to extract time-varying domain wall velocity and growth dimensionality even under arbitrary voltage waveforms, and that coupling this model to application-related waveforms and a circuit-level simulation platform creates a predictive materials-circuit co-design framework linking the nucleation and growth parameters to memory window, disturb error, speed, and energy dissipation.
What carries the argument
The dynamic-field-driven nucleation and growth (DFNG) model, which relaxes the constant-field assumption of earlier frameworks so that time-dependent electric fields can be incorporated when fitting switching data and extracting velocity and dimensionality parameters.
If this is right
- The model supports quantitative device modeling under complex signals that span disparate time and frequency scales.
- Nucleation and growth parameters extracted from transients can be linked directly to memory window and disturb error in circuit simulations.
- The framework enables predictive co-design by connecting material parameters to speed and energy dissipation for next-generation ferroelectric technologies.
- Fits remain possible even when the driving voltage waveform is arbitrary rather than a simple pulse.
Where Pith is reading between the lines
- The extracted parameters could be used to simulate and optimize voltage waveforms that minimize energy dissipation in larger memory arrays.
- The same fitting procedure might be applied to switching data from additional ferroelectric compositions to test how broadly the velocity and dimensionality trends generalize.
- Circuit-level simulations built on this model could reveal trade-offs between disturb immunity and write speed that are not apparent from material characterization alone.
Load-bearing premise
Parameters such as time-varying domain wall velocity extracted from fits to switching transients can be transferred directly into circuit simulations to predict quantities like memory window and disturb error without additional calibration or independent experimental checks.
What would settle it
Fabricate a ferroelectric capacitor or memory cell, apply a complex non-constant voltage waveform, measure its actual memory window and disturb error rates, then compare those measurements to the model's predictions that use only the parameters fitted from simple switching transients; systematic mismatch without refitting would show the transfer does not hold.
read the original abstract
Real ferroelectric devices operate under mixed and distorted time-varying voltages, yet the standard nucleation-growth frameworks used to interpret ferroelectric switching - most notably the Kolmogorov-Avrami-Ishibashi (KAI) and nucleation-limited switching models (NLS) - are derived under the critically limiting assumption of a constant electric field. Thus, the prevailing interpretation of ferroelectric switching dynamics fails under real operating conditions. Here we introduce a compact dynamic-field-driven nucleation and growth (DFNG) model that enables quantitative fits to switching transients across multiple ferroelectric materials to extract time-varying domain wall velocity and growth dimensionality, even under arbitrary voltage waveform. This capability then motivates its use in device modeling under complex signals spanning disparate time and frequency scales. Coupling the compact model to application-related waveforms and circuit-level simulation platform facilitates a predictive materials-circuit co-design framework by linking nucleation and growth parameters to memory window, disturb error, speed, and energy dissipation for next-generation ferroelectric technologies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a compact dynamic-field-driven nucleation and growth (DFNG) model for ferroelectric switching under arbitrary time-varying electric fields. It claims this enables quantitative fits to polarization switching transients across multiple materials to extract time-varying domain wall velocity v(t) and growth dimensionality n(t), overcoming limitations of constant-field models like KAI and NLS. The extracted parameters are then coupled to application waveforms and circuit-level simulations to enable predictive materials-to-circuit co-design, linking nucleation/growth parameters to device metrics including memory window, disturb error, speed, and energy dissipation.
Significance. If the transferability and predictive accuracy claims are substantiated, this work would be significant for ferroelectric device engineering. It directly addresses the mismatch between standard models (derived for constant E) and real-device operation under mixed/distorted signals spanning disparate timescales. The compact formulation and explicit linkage to SPICE-style circuit platforms represent a practical advance for co-design workflows in next-generation ferroelectric memories and related technologies. The multi-material applicability and focus on waveform-independent extraction are strengths if supported by clear validation.
major comments (2)
- [§4] §4 (Circuit co-design and simulation results): The central predictive claim—that parameters v(t) and n(t) extracted from fits to switching transients can be directly transferred to forecast memory window and disturb error under complex application waveforms—lacks independent experimental validation. The presented circuit predictions appear generated from the same fitted quantities used for the transient fits, creating a circularity risk for the co-design framework as the model may embed waveform-specific effects (e.g., nucleation statistics or inhomogeneity) not explicitly separated.
- [Model formulation] Model formulation and fitting procedure (around the DFNG equations): While the model relaxes the constant-field assumption, the manuscript does not demonstrate that the extracted time-varying v(t) and n(t) remain invariant when the input waveform changes in frequency content or shape. A concrete cross-validation test—fitting on one set of transients and predicting an independent set under a dissimilar waveform with quantitative error metrics—is required to support the transferability needed for the co-design application.
minor comments (3)
- [Figures] Figure captions and axis labels should explicitly state the voltage waveforms (e.g., pulse shape, frequency range) used for each transient fit and simulation to improve reproducibility.
- [Abstract/Introduction] The abstract and introduction would benefit from a brief quantitative example (e.g., one material, one error metric) rather than only qualitative assertions about 'quantitative fits'.
- [Methods] Notation for time-dependent quantities (v(t), n(t)) should be consistently distinguished from their constant-field counterparts in the methods or appendix to avoid reader confusion.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments, which help clarify the scope and limitations of our claims. We address each major point below, indicating where revisions will be made to improve the manuscript.
read point-by-point responses
-
Referee: [§4] §4 (Circuit co-design and simulation results): The central predictive claim—that parameters v(t) and n(t) extracted from fits to switching transients can be directly transferred to forecast memory window and disturb error under complex application waveforms—lacks independent experimental validation. The presented circuit predictions appear generated from the same fitted quantities used for the transient fits, creating a circularity risk for the co-design framework as the model may embed waveform-specific effects (e.g., nucleation statistics or inhomogeneity) not explicitly separated.
Authors: We agree that the current presentation risks appearing circular because the circuit-level forecasts in §4 rely on parameters fitted to the same class of transients. The DFNG formulation is intended to be waveform-agnostic by extracting instantaneous v(t) and n(t) from the local field history, but this transferability is not yet demonstrated against independent device-level measurements under complex signals. In the revised manuscript we will add a dedicated paragraph in §4 explicitly stating that the reported memory-window and disturb-error results are model-based projections rather than experimentally validated predictions, and we will include a forward-looking statement on the experimental data sets required for full validation. revision: yes
-
Referee: [Model formulation] Model formulation and fitting procedure (around the DFNG equations): While the model relaxes the constant-field assumption, the manuscript does not demonstrate that the extracted time-varying v(t) and n(t) remain invariant when the input waveform changes in frequency content or shape. A concrete cross-validation test—fitting on one set of transients and predicting an independent set under a dissimilar waveform with quantitative error metrics—is required to support the transferability needed for the co-design application.
Authors: The referee correctly notes that invariance of the extracted v(t) and n(t) across waveform shapes has not been shown. Although the model equations are written to accept arbitrary E(t), the manuscript presents fits only to the waveforms used for parameter extraction. We will therefore add a new subsection (or appendix) that performs the requested cross-validation: the DFNG parameters will be fitted to transients recorded under one waveform family and then used to predict the polarization response under a second, dissimilar waveform (different frequency content and shape). Quantitative metrics (RMSE and maximum deviation) between predicted and measured transients will be reported to quantify the degree of transferability. revision: yes
Circularity Check
Fitted DFNG parameters directly linked to circuit predictions without independent validation
specific steps
-
fitted input called prediction
[Abstract]
"This capability then motivates its use in device modeling under complex signals spanning disparate time and frequency scales. Coupling the compact model to application-related waveforms and circuit-level simulation platform facilitates a predictive materials-circuit co-design framework by linking nucleation and growth parameters to memory window, disturb error, speed, and energy dissipation for next-generation ferroelectric technologies."
Nucleation and growth parameters (time-varying domain wall velocity and dimensionality) are obtained via quantitative fits to switching transients; these same parameters are then directly linked to circuit quantities in the co-design framework. The circuit predictions are therefore generated from the fitted inputs by construction rather than from independent measurements on the target waveforms.
full rationale
The paper extracts time-varying domain wall velocity and growth dimensionality by fitting the DFNG model to polarization switching transients under given waveforms. These fitted quantities are then used without additional calibration or separate experimental checks to predict circuit-level metrics such as memory window and disturb error under complex application waveforms. This matches the fitted-input-called-prediction pattern: the co-design outputs are generated from the same fitted parameters, reducing the predictive claim to model application rather than new empirical support. The derivation chain is therefore partially circular at the materials-to-circuit transfer step, though the underlying model equations themselves are not self-referential.
Axiom & Free-Parameter Ledger
free parameters (2)
- time-varying domain wall velocity
- growth dimensionality
axioms (1)
- domain assumption The nucleation and growth process in ferroelectrics remains describable by the same physical framework when the electric field is time-varying rather than constant.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We build our derivation from Cahn’s time-cone model... ⟨N(t)⟩ = ... ∫ J(τ) [∫ v(t') dt']^d dτ ... v(t') = v0 exp(−Va/V(t')) ... Merz law
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The DFNG model ... links nucleation and growth parameters to memory window, disturb error, speed, and energy dissipation
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Determinism under CMOS conditions: 99% switching fraction should be completed within 10 ns (delay) with a pulse amplitude (Vmin) ≤ 1.5 V4
-
[2]
Memory stability against disturb: characterized by the polarization loss fraction Δ𝑓𝑓 under a write disturb pulse in the V DD/2 scheme 30. Nominally Δ𝑓𝑓≤ 5% is acceptable for memory device such as FeRAM (ferroelectric random-access memory), and Δ𝑓𝑓≤ 1% for CiM devices as they encounter more frequent write operations
-
[3]
Energy and delay for CiM under write operation
-
[4]
𝜖𝜖𝐹𝐹𝐹𝐹 is the dielectric constant of the ferroelectric (subtracting the linear part)
Memory window for FeCaps 31: 𝐶𝐶𝑛𝑛𝑛𝑛𝑛𝑛𝑛𝑛= � 𝑑𝑑𝑑𝑑 𝑑𝑑𝑉𝑉� 𝑑𝑑=0.5 = 1 𝑃𝑃𝑆𝑆𝑡𝑡𝐹𝐹𝐹𝐹 𝜖𝜖𝐹𝐹𝐹𝐹. 𝜖𝜖𝐹𝐹𝐹𝐹 is the dielectric constant of the ferroelectric (subtracting the linear part)
-
[5]
We locate the loop of the current HZO sample according to its fit parameters (Fig
Read margin / Memory window for FeFETs32: VC = Vf=0.5. We locate the loop of the current HZO sample according to its fit parameters (Fig. 4a, purple star). This immediately suggests that a modest increase in nucleation density and growth dimension would adapt the material for a CiM application (orange star). Physically, such tuning could be achieved by in...
-
[6]
Martin, L. W. & Rappe, A. M. Thin-film ferroelectric materials and their applications. Nat. Rev. Mater. 2, 16087 (2016)
work page 2016
- [7]
-
[8]
Wang, J. et al. Ferroelectric domain-wall logic units. Nat. Commun. 13, 3255 (2022)
work page 2022
-
[9]
Khan, A. I., Keshavarzi, A. & Datta, S. The future of ferroelectric field-effect transistor technology. Nat. Electron. 3, 588–597 (2020)
work page 2020
-
[10]
Mikolajick, T., Schroeder, U. & Slesazeck, S. The Past, the Present, and the Future of Ferroelectric Memories. IEEE Trans. Electron Devices 67, 1434–1443 (2020)
work page 2020
-
[11]
Beyond CMOS and Emerging Materials Integration
IEEE International Roadmap for Devices and Systems. Beyond CMOS and Emerging Materials Integration. (2023) doi:10.60627/0P45-ZJ55
-
[12]
Kolmogorov, A. N. On the Statistical Theory of Crystallization of Metals [in Russian]. Izv Akad Nauk SSSR Ser Mat 3, 355–359 (1937)
work page 1937
-
[13]
Avrami, M. Kinetics of Phase Change. I General Theory. J. Chem. Phys. 7, 1103–1112 (1939)
work page 1939
-
[14]
Ishibashi, Y . & Takagi, Y . Note on Ferroelectric Domain Switching. J. Phys. Soc. Jpn. 31, 506–510 (1971)
work page 1971
-
[15]
Tagantsev, A. K. Non-Kolmogorov-Avrami switching kinetics in ferroelectric thin films. Phys. Rev. B 66, (2002)
work page 2002
-
[16]
L. D. Landau & I. M. Khalatnikov. On the anomalous absorption of sound near a second- order phase transition point. Dokl. Akad. Nauk SSSR 469–472 (1954)
work page 1954
-
[17]
Devonshire, A. F. XCVI. Theory of barium titanate: Part I. Lond. Edinb. Dublin Philos. Mag. J. Sci. 40, 1040–1063 (1949)
work page 1949
-
[18]
Devonshire, A. F. CIX. Theory of barium titanate—Part II. Lond. Edinb. Dublin Philos. Mag. J. Sci. 42, 1065–1079 (1951)
work page 1951
-
[19]
Merz, W. J. Domain Formation and Domain Wall Motions in Ferroelectric BaTiO3 Single Crystals. Phys. Rev. 95, 690–698 (1954)
work page 1954
-
[20]
Cahn, J. W. The Time Cone method for Nucleation and Growth Kinetics on a Finite Domain. MRS Online Proc. Libr. 398, 425–437 (1995)
work page 1995
-
[21]
On the theory of the coercive field of single-domain crystals of BaTiO3
Janovec, V . On the theory of the coercive field of single-domain crystals of BaTiO3. Cechoslov- Fiz. Zurnal 8, 3–15 (1958)
work page 1958
-
[22]
Kay, H. F. & Dunn, J. W. Thickness dependence of the nucleation field of triglycine sulphate. Philos. Mag. J. Theor. Exp. Appl. Phys. 7, 2027–2034 (1962)
work page 2027
-
[23]
Chiang, T. et al. Material-Limited Switching in Nanoscale Ferroelectrics. Preprint at https://doi.org/10.48550/arXiv.2507.12353 (2025)
-
[24]
Gao, P. et al. Revealing the role of defects in ferroelectric switching with atomic resolution. Nat. Commun. 2, 591 (2011)
work page 2011
-
[25]
Zhang, L. et al. Mapping electric fields and observation of ferroelectric domain switching in hafnia-zirconia devices by electron holography. Nat. Commun. 16, 11233 (2025)
work page 2025
- [26]
-
[27]
Nelson, C. T. et al. Domain Dynamics During Ferroelectric Switching. Science 334, 968– 971 (2011)
work page 2011
-
[28]
Fields, S. S. et al. Compositional and phase dependence of elastic modulus of crystalline and amorphous Hf1-xZrxO2 thin films. Appl. Phys. Lett. 118, 102901 (2021)
work page 2021
-
[29]
Beyer, S. et al. FeFET: A versatile CMOS compatible device with game-changing potential. in 2020 IEEE International Memory Workshop (IMW) 1–4 (2020). doi:10.1109/IMW48823.2020.9108150
-
[30]
Mukherjee, S. et al. Pulse-Based Capacitive Memory Window with High Non- Destructive Read Endurance in Fully BEOL Compatible Ferroelectric Capacitors. in 2023 International Electron Devices Meeting (IEDM) 1–4 (2023). doi:10.1109/IEDM45741.2023.10413879
-
[31]
Luo, Z. et al. High-precision and linear weight updates by subnanosecond pulses in ferroelectric tunnel junction for neuro-inspired computing. Nat. Commun. 13, 699 (2022)
work page 2022
- [32]
-
[33]
Yu, E., K, G. K., Saxena, U. & Roy, K. Ferroelectric capacitors and field-effect transistors as in-memory computing elements for machine learning workloads. Sci. Rep. 14, 9426 (2024)
work page 2024
-
[34]
Luo, Y .-C. et al. Experimental Demonstration of Non-volatile Capacitive Crossbar Array for In-memory Computing. in 2021 IEEE International Electron Devices Meeting (IEDM) 1–4 (2021). doi:10.1109/IEDM19574.2021.9720508
-
[35]
Mueller, S. et al. From MFM Capacitors Toward Ferroelectric Transistors: Endurance and Disturb Characteristics of HfO 2-Based FeFET Devices. IEEE Trans. Electron Devices 60, 4199– 4205 (2013)
work page 2013
-
[36]
Zhang, X. et al. Single-crystalline BaTiO3-based ferroelectric capacitive memory via membrane transfer. Sci. Adv. 11, eadz2553 (2025)
work page 2025
-
[37]
Toprasertpong, K., Takenaka, M. & Takagi, S. Memory Window in Ferroelectric Field- Effect Transistors: Analytical Approach. IEEE Trans. Electron Devices 69, 7113–7119 (2022)
work page 2022
-
[38]
Pramanick, A., Prewitt, A. D., Forrester, J. S. & Jones, J. L. Domains, Domain Walls and Defects in Perovskite Ferroelectric Oxides: A Review of Present Understanding and Recent Contributions. Crit. Rev. Solid State Mater. Sci. 37, 243–275 (2012)
work page 2012
-
[39]
Shin, Y .-H., Grinberg, I., Chen, I.-W. & Rappe, A. M. Nucleation and growth mechanism of ferroelectric domain-wall motion. Nature 449, 881–884 (2007)
work page 2007
-
[40]
Liu, S., Grinberg, I. & Rappe, A. M. Intrinsic ferroelectric switching from first principles. Nature 534, 360–363 (2016)
work page 2016
-
[41]
Theoretical Lower Limit of Coercive Field in Ferroelectric Hafnia
Yang, J. Theoretical Lower Limit of Coercive Field in Ferroelectric Hafnia. Phys. Rev. X 15, (2025)
work page 2025
-
[42]
Lenox, M. K., Jaszewski, S. T., Fields, S. S., Shukla, N. & Ihlefeld, J. F. Impact of Electric Field Pulse Duration on Ferroelectric Hafnium Zirconium Oxide Thin Film Capacitor Endurance. Phys. Status Solidi A 221, 2300566 (2024). Acknowledgements Y .L., T. C., and J. T. H. acknowledge partial support from the Intel FEINMAN 2.0 program. Y .L. and J. T. H....
work page 2024
-
[43]
The conventional form of the KAI model only applies to the constant electric field case
Deficiency of KAI model in capturing switching behavior under time-dependent voltage KAI model describes the transformed fraction f of polarization at time t by: 𝑓𝑓(𝑡𝑡) = 1 − exp �− � 𝑡𝑡 𝑡𝑡0 � 𝑛𝑛 � (𝑆𝑆1) where t0 is the characteristic time related to voltage-independent nucleation rates and domain wall velocities, and n is the Avrami exponent related to d...
-
[44]
Detailed derivation of DFNG model with voltage-dependent homogeneous nucleation rate The voltage-dependent homogeneous nucleation rate is derived from the Janovec-Kay-Dunn model, which gives the probability of nucleation per site governed by thermodynamic barrier by exp � − Δ𝐺𝐺∗ 𝑘𝑘𝐵𝐵𝑇𝑇� and Δ𝐺𝐺∗ is the critical Gibbs free energy to create a nucleus. Assum...
work page 1918
-
[45]
The extraction of material parameters for BTO Following the fit procedure in Methods, the experimentally measured switching transients of a 3 μm BTO capacitor is fit to the model, as shown in Fig. S2. The extracted parameters are 𝑑𝑑= 1.144 ± 0.004, 𝑉𝑉𝑎𝑎= 0.536 ± 0.003 𝑉𝑉, 𝜎𝜎′2 = 1.47 ± 0.04 𝑉𝑉, 𝐴𝐴= 6.13 ± 0.07 ns−1, and 𝐵𝐵= 713 ± 44 ns−1 d−1. Specifically...
-
[46]
The truncated time cone also accounts for the deviation between the model and the experimental data when approaching the end of the switching. Besides a higher growth dimension compared to HZO, BTO also exhibits a much lower Merz barrier (V a) and heterogeneous nucleation density (A). The domain wall energy 𝜎𝜎 of BTO is derived to be ~4.4 mJ/m 2, accordin...
-
[47]
The method to generate the polarization switching during a voltage ramp is shown in Fig
Parameter space exploration and application metrics extraction 4.1 Simulation of hysteresis loops Leakage, linear dielectric constant, and depolarization in the ferroelectric capacitor are ignored and a preset -PS state is assumed. The method to generate the polarization switching during a voltage ramp is shown in Fig. S3a and S3b, which is essentially ha...
-
[48]
Parameter sets filtering for FeFETs and FeCaps We use a nominal criterion of VC ≥ 1.4 V to select favorable FeFETs gate materials. It is inferred from Fig. 4b that only the delay has a broad distribution and needs consideration under this selection. Table S1 is a table of the satisfying parameter sets, in the order of increasing delay. The parameter set t...
-
[49]
On the theory of the coercive field of single-domain crystals of BaTiO3
Janovec, V . On the theory of the coercive field of single-domain crystals of BaTiO3. Czech J Phys 8, 3–15 (1958)
work page 1958
-
[50]
Kay, H. F. & Dunn, J. W. Thickness dependence of the nucleation field of triglycine sulphate. The Philosophical Magazine: A Journal of Theoretical Experimental and Applied Physics 7, 2027–2034 (1962)
work page 2027
-
[51]
Deck, L.-T. & Mazzotti, M. Conceptual Validation of Stochastic and Deterministic Methods To Estimate Crystal Nucleation Rates. Crystal Growth & Design 23, 899–914 (2023)
work page 2023
-
[52]
in Kinetics of Materials 533–542 (John Wiley & Sons, Ltd, 2005)
Concurrent Nucleation and Growth. in Kinetics of Materials 533–542 (John Wiley & Sons, Ltd, 2005). doi:10.1002/0471749311.ch21
-
[53]
Ab initio study of ferroelectric domain walls in PbTiO3
Meyer, B. Ab initio study of ferroelectric domain walls in PbTiO3. Phys. Rev. B 65, (2002)
work page 2002
-
[54]
Hlinka, J. et al. Soft-mode spectroscopy of epitaxial BaTiO3/SrTiO3 superlattices. Phys. Rev. B 82, 224102 (2010)
work page 2010
-
[55]
Kashida, S., Hatta, I., Ikushima, A. & Yamada, Y . Ultrasonic Velocities in BaTiO3. J. Phys. Soc. Jpn. 34, 997–1001 (1973)
work page 1973
- [56]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.