pith. sign in

arxiv: 2604.05184 · v1 · submitted 2026-04-06 · ❄️ cond-mat.mtrl-sci · physics.app-ph

Zr Concentration-Dependent Sub-Lattice Phase-Field Model of Hf1-xZrxO2: Analysis of Phase Composition and Polarization Switching

Tae Ryong Kim , Sumeet K. Gupta This is my paper

Pith reviewed 2026-05-10 18:55 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.app-ph
keywords Hf1-xZrxO2phase-field modelferroelectricanti-ferroelectricpolarization switchingzirconium concentrationmixed phasecharge-voltage characteristics
0
0 comments X

The pith

Zr concentration creates mixed orthorhombic-tetragonal phases that cause gradual polarization reversal in Hf1-xZrxO2 at x=0.7-0.8

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors build a sub-lattice phase-field model of Hf1-xZrxO2 by extending the time-dependent Ginzburg-Landau equation with zirconium-fraction-dependent interaction parameters and gradient coefficients. Calibrated to measured charge-voltage curves, the model tracks how the orthorhombic phase favors ferroelectric switching at lower x while the tetragonal phase stabilizes anti-ferroelectric behavior at higher x. At intermediate concentrations the two phases have nearly equal energies, so local electric-field variations near domain walls produce mixed-phase regions and spatially staggered reversal. This yields the smoother voltage response seen experimentally in that window. The work supplies a spatially resolved picture that links bulk phase stability to observable device characteristics.

Core claim

Our sub-lattice phase-field model expands the time-dependent Ginzburg-Landau equation to the sub-lattice level and incorporates x-dependent interaction parameters and gradient coefficients. It reproduces the ferroelectric-to-anti-ferroelectric transition as x rises from 0.5 to 1.0. At low x the orthorhombic phase dominates and yields distinct ferroelectric loops; at high x the tetragonal phase is stabilized and anti-ferroelectric transitions appear. At intermediate x of 0.7-0.8 the comparable phase energies make the system sensitive to local electric-field variations from stray fields near domain walls, producing mixed-phase composition and spatially staggered polarization reversal that in t

What carries the argument

Sub-lattice phase-field model extending the time-dependent Ginzburg-Landau equation with x-dependent parameters, which resolves thermodynamic preferences, kinetic barriers, and spatially non-uniform polarization and electric-field profiles between orthorhombic and tetragonal phases.

Load-bearing premise

The x-dependent interaction parameters and gradient coefficients, once calibrated to measured charge-voltage curves, correctly capture the thermodynamic preference and kinetic barriers between orthorhombic and tetragonal phases without further hidden adjustments.

What would settle it

Spatially resolved diffraction or imaging at x=0.75 that shows either uniform single-phase orthorhombic or tetragonal regions together with sharp rather than gradual polarization reversal would contradict the predicted mixed-phase behavior driven by local field variations.

Figures

Figures reproduced from arXiv: 2604.05184 by Sumeet K. Gupta, Tae Ryong Kim.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Microscopic schematics of [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Calibration of model (line) for different Zr concentrations [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) The calibrated [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Evolution of [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Phase ( [PITH_FULL_IMAGE:figures/full_fig_p005_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Local [PITH_FULL_IMAGE:figures/full_fig_p006_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. (a) Local [PITH_FULL_IMAGE:figures/full_fig_p006_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Phase transition plots from [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
read the original abstract

We develop a sub-lattice phase-field model of Hf1-xZrxO2 incorporating zirconium (Zr) concentration (x)-dependence. Our framework expands the time-dependent Ginzburg-Landau (TDGL) equation to the sub-lattice level and incorporates x-dependent interaction parameters and gradient coefficients. Our experimentally calibrated model captures the evolution of charge-voltage (Q-V) characteristics for x ranging from 0.5 to 1.0. The sub-lattice formulation explains the thermodynamic preference and kinetic transition barriers of competing orthorhombic phase (o-phase) and tetragonal phase (t-phase), while the phase-field framework enables spatially resolved analysis of polarization (P) and electric-field (E-field) profiles, allowing multi-domain (MD) polarization and mixed-phase states to emerge naturally. Our model reproduces the experimentally observed ferroelectric (FE)-to-anti-ferroelectric (AFE) transition as x increases from 0.5 to 1.0. At low Zr concentration (x = 0.5-0.6), the o-phase dominates, yielding distinct FE behavior. At high concentration (x = 0.9-1.0), the t-phase is stabilized, leading to AFE transitions. A key finding of our work is the unique behavior at intermediate Zr concentrations (x = 0.7-0.8). Here, the o- and t-phase energies are comparable, making the system strongly influenced by local variations in the electric field (E-field), which arise from stray fields near the domain walls. This non-uniform field distribution results in a mixed-phase composition and spatially staggered polarization reversal, which manifests as a more gradual Q-V evolution compared to other values of x. By linking energy landscapes to spatial field effects, the model provides insights into the FE-to-AFE crossover in Hf1-xZrxO2.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript develops a sub-lattice phase-field model for Hf_{1-x}Zr_xO_2 by extending the time-dependent Ginzburg-Landau equation to incorporate Zr concentration (x) dependence through x-dependent interaction parameters and gradient coefficients. The experimentally calibrated model reproduces charge-voltage (Q-V) characteristics across x from 0.5 to 1.0, explains the ferroelectric-to-antiferroelectric transition, and attributes mixed-phase states and staggered polarization reversal at intermediate x (0.7-0.8) to comparable orthorhombic and tetragonal phase energies influenced by local electric-field variations near domain walls.

Significance. If the phase-energy comparability and resulting mixed-phase behavior can be shown to arise independently of the Q-V calibration, the sub-lattice TDGL framework with spatial resolution would offer useful mechanistic insight into the FE-AFE crossover in hafnia-based materials, which is relevant for ferroelectric device applications. The approach of allowing multi-domain and mixed-phase states to emerge naturally is a positive feature, but the current calibration procedure limits the strength of the central mechanistic claims.

major comments (2)
  1. [Abstract] Abstract and model formulation: The x-dependent interaction parameters and gradient coefficients are calibrated to experimental Q-V data to reproduce the FE-to-AFE transition. This makes the reported comparability of o- and t-phase energies specifically at x=0.7-0.8 a likely direct outcome of the fitting procedure rather than an independent prediction of the sub-lattice dynamics. The claim that local E-field variations then produce mixed-phase composition and staggered reversal therefore rests on an energy landscape that was adjusted to match the very observations being explained; independent constraints (e.g., composition-dependent formation energies from first-principles) are not described.
  2. [Model formulation] Model section (parameterization of x-dependence): The manuscript must specify the functional form and number of free parameters used for the x-dependent terms, and demonstrate that these terms are not solely determined by Q-V fitting. Without such details or additional validation against non-Q-V observables, the thermodynamic preference and kinetic barriers between o- and t-phases cannot be regarded as robustly captured by the sub-lattice TDGL model.
minor comments (1)
  1. [Abstract] The abstract would benefit from a brief statement clarifying which aspects of the phase energies and switching behavior are direct outputs of the calibrated model versus which are interpretive interpretations of the spatial field profiles.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [Abstract] Abstract and model formulation: The x-dependent interaction parameters and gradient coefficients are calibrated to experimental Q-V data to reproduce the FE-to-AFE transition. This makes the reported comparability of o- and t-phase energies specifically at x=0.7-0.8 a likely direct outcome of the fitting procedure rather than an independent prediction of the sub-lattice dynamics. The claim that local E-field variations then produce mixed-phase composition and staggered reversal therefore rests on an energy landscape that was adjusted to match the very observations being explained; independent constraints (e.g., composition-dependent formation energies from first-principles) are not described.

    Authors: We agree that the x-dependent parameters are determined through calibration to experimental Q-V data to reproduce the overall FE-to-AFE transition. The comparability of o- and t-phase energies at x=0.7-0.8 is indeed a consequence of the x-dependence chosen to enable the transition. However, the mixed-phase composition and staggered polarization reversal arise as emergent behavior from the spatially resolved sub-lattice phase-field dynamics, specifically through the influence of local E-field variations near domain walls; these spatial features are not directly imposed by the fitting but follow from the model structure once the energies are comparable. We will revise the abstract and relevant sections to clarify this distinction between the calibrated global transition and the predicted local mechanistic effects. First-principles constraints are not included in the present phenomenological model. revision: partial

  2. Referee: [Model formulation] Model section (parameterization of x-dependence): The manuscript must specify the functional form and number of free parameters used for the x-dependent terms, and demonstrate that these terms are not solely determined by Q-V fitting. Without such details or additional validation against non-Q-V observables, the thermodynamic preference and kinetic barriers between o- and t-phases cannot be regarded as robustly captured by the sub-lattice TDGL model.

    Authors: We accept that the current manuscript lacks explicit detail on the functional forms and parameter count. In the revised version we will expand the Model section to state the precise functional forms (e.g., the polynomial or linear interpolations) used for the x-dependent interaction parameters and gradient coefficients, together with the number of free parameters adjusted during calibration. The sub-lattice TDGL formulation encodes the thermodynamic preferences and kinetic barriers through the free-energy landscape and the evolution equations; the calibration to Q-V data across multiple x values constrains these quantities consistently with experiment. Additional validation against independent non-Q-V observables is not provided in the present study. revision: yes

Circularity Check

1 steps flagged

x-dependent parameters calibrated to Q-V data impose comparable o/t energies at intermediate concentrations

specific steps
  1. fitted input called prediction [Abstract]
    "Our experimentally calibrated model captures the evolution of charge-voltage (Q-V) characteristics for x ranging from 0.5 to 1.0. ... A key finding of our work is the unique behavior at intermediate Zr concentrations (x = 0.7-0.8). Here, the o- and t-phase energies are comparable, making the system strongly influenced by local variations in the electric field (E-field), which arise from stray fields near the domain walls."

    The x-dependent parameters are adjusted to match the experimental Q-V data that exhibits the FE-to-AFE crossover and gradual evolution at x=0.7-0.8. Declaring that the o- and t-phase energies are comparable at those concentrations is therefore a direct consequence of the calibration that was performed to reproduce the same Q-V behavior, rather than a prediction derived from independent inputs.

full rationale

The paper introduces x-dependent interaction parameters and gradient coefficients in the sub-lattice TDGL model, then calibrates them to reproduce measured Q-V curves across x=0.5-1.0. The central claim—that o- and t-phase energies become comparable specifically at x=0.7-0.8, enabling mixed-phase states and staggered reversal—is presented as an emergent finding from the calibrated energy landscape. Because the calibration targets the very Q-V evolution that the transition produces, the comparability of phase energies reduces to an outcome of the fitting procedure rather than an independent first-principles result. No external constraint (e.g., ab-initio formation energies) is invoked to fix the x-dependence independently of the target data. The spatial phase-field analysis of stray fields operates on this fitted landscape, so the mechanistic explanation inherits the circularity. This matches the 'fitted input called prediction' pattern at the load-bearing step.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on the standard TDGL evolution equation plus two new x-dependent functions (interaction parameters and gradient coefficients) whose values are obtained by fitting to measured charge-voltage curves rather than derived from first principles or independent measurements.

free parameters (2)
  • x-dependent interaction parameters
    Made functions of zirconium concentration and adjusted so that the model matches experimental Q-V loops for x from 0.5 to 1.0.
  • x-dependent gradient coefficients
    Introduced to control interface energy between phases and likewise calibrated to reproduce observed phase stability.
axioms (2)
  • standard math The time-dependent Ginzburg-Landau equation governs the evolution of the order parameter at the sub-lattice level.
    Invoked in the first paragraph as the starting point for the sub-lattice expansion.
  • domain assumption Orthorhombic and tetragonal phases are the only two competing structures whose relative free energies determine the observed ferroelectric or antiferroelectric response.
    Stated when the model is said to explain the thermodynamic preference of o-phase versus t-phase.

pith-pipeline@v0.9.0 · 5664 in / 1794 out tokens · 50875 ms · 2026-05-10T18:55:06.201875+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages

  1. [1]

    merlin.mbs aapmrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs aapmrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translat...

    work page 2010

  2. [2]

    merlin.mbs aipauth4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs aipauth4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translat...

    work page 2010

  3. [3]

    merlin.mbs aipnum4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs aipnum4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translati...

    work page 2010

  4. [4]

    merlin.mbs apsrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs apsrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translati...

    work page 2010

  5. [5]

    merlin.mbs apsrmp4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs apsrmp4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translati...

    work page 2010