Pith. sign in

REVIEW 3 major objections 6 minor 20 references

A General Model of Interfacial Chemical Equilibrium in Phase-Field Method

T0 review · 3 major / 6 minor · reviewed 2026-07-31 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Auxiliary non-conserved variables unify interfacial chemical equilibrium in phase-field models and treat phases with no shared composition range.

desk verdict Useful auxiliary-variable bridge between WBM and KKS that actually handles non-overlapping ordered phases with raw CALPHAD free energies; the KKS limit rests on a pragmatic, under-justified kinetic tweak. read the letter →

arxiv 2607.24620 v1 pith:O2MH6UON submitted 2026-07-27 cond-mat.mtrl-sci

Yanzhou Ji , Long-Qing Chen This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords Phase-FieldmodelInterfacialcompositionsequilibriumMicrostructureevolutionthermodynamicdatabasesorderedphasesdiffusionalphasetransformations
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Phase-field models of diffusional microstructure evolution must decide how to split a local composition among the phases that coexist inside a diffuse interface. The classic equal-composition and equal-chemical-potential assumptions either fail for ordered phases whose stability ranges do not overlap or require expensive nonlinear solves and free-energy approximations. This paper introduces auxiliary non-conserved fields that simply track the composition differences between phases and evolve them by the mismatch of diffusional chemical potentials. The construction recovers the two classic models as opposite kinetic limits, accepts unsimplified thermodynamic database expressions, and works for interdiffusion between non-stoichiometric ordered phases that share no common composition range. The result is a single, flexible framework for composition evolution and interfacial chemical equilibrium across a wide class of multiphase materials problems.

What carries the argument

The auxiliary variable φ = X_B^α − X_B^β (and its multi-phase, multi-component generalizations), whose evolution is driven by the difference of diffusional chemical potentials and whose kinetic coefficient continuously interpolates between equal-composition and equal-potential interfacial conditions.

What would settle it

Run the auxiliary-variable model on θ/η1 interdiffusion in Al–Cu (or an analogous pair of ordered phases with disjoint composition ranges) and check whether the final compositions and common-tangent chemical potentials disagree with independent thermodynamic calculations or with a carefully regularized equal-potential solution; any systematic mismatch falsifies the claim.

Watch

Extended reading notes

Core claim

Interfacial chemical equilibrium among solution phases is described by auxiliary non-conserved variables that measure composition differences between phases; these variables evolve by Allen-Cahn dynamics driven by differences in chemical diffusional potentials, recovering the Wheeler-Boettinger-McFadden equal-composition model and the Kim-Kim-Suzuki equal-potential model as limiting cases while remaining valid for phases whose composition ranges do not overlap and while using thermodynamic databases without further approximation.

Load-bearing premise

Dropping the interfacial weighting factor from the thermodynamically derived evolution law for the auxiliary variable still yields a faithful, artifact-free description of chemical equilibrium in both bulk and interface.

Editorial extensions

If this is right

  • Interdiffusion between non-stoichiometric ordered phases or oxides can be simulated directly from thermodynamic databases without free-energy fitting or extrapolation.
  • Multicomponent multiphase systems are obtained by defining one auxiliary variable per independent component per phase pair.
  • Three-dimensional simulations become cheaper than the equal-potential model because interface nonlinear solves are replaced by a simple relaxation step.
  • The continuous family of interfacial conditions parameterized by the auxiliary kinetic coefficient can be used to explore how real interfaces sit between the two classic extremes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same auxiliary-variable construction can be coupled to elastic or electrostatic fields to treat chemo-mechanical interfacial equilibria without forcing equal compositions.
  • The logarithmic intermediate variables introduced for the Al–Cu example supply a general numerical safeguard for any sublattice free energy that contains logarithms of site fractions.
  • Measuring how far experimental interface profiles lie from the equal-composition versus equal-potential limits would give a direct experimental calibration of the auxiliary kinetic coefficient.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript introduces an auxiliary non-conserved field φ equal to the composition difference between two phases, so that the two phase compositions are determined from the conserved composition, phase-field variable, and φ. Allen-Cahn/Cahn-Hilliard dynamics are supplemented by a relaxation equation for φ driven by the difference between phase diffusion potentials. With φ=0 and Lφ=0 the formulation gives the WBM equal-composition limit; sufficiently fast φ relaxation is presented as a KKS equal-potential limit. A variational equation (8) and a modified equation (9), obtained by deleting the interfacial factor h(1-h), are both considered. Ti-V growth is benchmarked against WBM and KKS, while an Al-Cu θ/η1 example uses log-transformed variables and database free energies to treat phases with disjoint composition domains. A supplementary 3-D test reports removal of most interfacial extra potential and a 3.3× speedup over KKS.

Significance. If the remaining dynamical and reproducibility issues are resolved, the model would be a useful addition to phase-field methodology. It gives a transparent interpolation between equal-composition and equal-diffusion-potential interface conditions, derives the variational φ equation explicitly, and avoids refitting CALPHAD free energies in the demonstrated examples. The log-variable construction for θ/η1 is a practical, domain-preserving treatment, and the recovery of common-tangent compositions in Al-Cu is a meaningful test. The potential 3-D efficiency advantage over Newton-based KKS implementations is also important, although the present evidence for that advantage is incomplete.

major comments (3)
  1. [§3.2 and Supplement S4, Eqs. (8), (12), (S3-6)-(S3-7)] §3.2/S4, Eqs. (8), (12a)-(12b): the stated derivation does not follow from Eq. (8). Inserting φ̇ from Eq. (8) into Eqs. (S3-6a,b) gives factors -h²(1-h)LφΔμ/vm in the Y1 equation and +h(1-h)²LφΔμ/vm in the Y2 equation. The printed Eqs. (12)/S3-7 instead contain -hLφΔμ/vm and +(1-h)LφΔμ/vm, which correspond to Eq. (9). Thus either Fig. 5 was not generated using the equation stated, or Eq. (12) is algebraically incorrect. Please identify the implemented equation, correct S4, and rerun if necessary.
  2. [§2, Eq. (9)] The deletion of h(ξ)[1-h(ξ)] is load-bearing for the claimed KKS limit and system-wide equalization of diffusion potentials, but Eq. (9) is no longer the constant-mobility Allen-Cahn equation for δF/δφ. Relative to Eq. (8) it is equivalent to a position-dependent mobility Lφ/[h(1-h)] that diverges toward the interface fringes, and it evolves virtual phase compositions even where the corresponding phase has zero weight. The manuscript should provide a free-energy-dissipation or stability analysis, explain treatment of CALPHAD domains for virtual compositions, and test timestep/grid sensitivity.
  3. [Figs. 2-4 and Supplement S3] Finite Lφ introduces an additional interfacial kinetic time scale. Figure 2(d) shows that composition profiles depend on Lφ, but no convergence of interface velocity or growth law is demonstrated, and the statement that model choice does not significantly affect kinetics is based on one thin-interface example. The S3 efficiency test uses Lφ=Lξ and Eq. (8), while the >99.99% extra-potential reduction is not precisely defined. Please report Lφ convergence at fixed Lξ, M_X, w, and κξ; Δt/Δx convergence; quantitative comparison with KKS; and guidance for choosing or calibrating Lφ.
minor comments (6)
  1. [§2 vs. §3.1, Figs. 2-3] §3.1 reverses the convention introduced in §2: there ξ=0 denotes α and ξ=1 denotes β, whereas the Ti-V simulation uses α-Ti at ξ=1 and β-Ti at ξ=0. With Eq. (2), the stated initial and equilibrium Ti-V compositions imply φ≈-0.0677 and -0.2151, but the text and Fig. 3 report positive values. Please explicitly state the variable relabeling and sign convention used in the simulation.
  2. [Figure 5 caption] The Fig. 5(d) caption says that φ profiles were obtained “using Eq. (7b),” which is the composition equation. This presumably should refer to Eq. (8) or Eq. (9); correcting it is especially important in view of the S4 inconsistency.
  3. [Supplement S3 and Fig. S3-1] For the S3 timing comparison, please specify hardware, timestep, Newton tolerance and maximum iteration count, whether CPU time is total or per time step, and whether the KKS and present-model implementations received comparable optimization. These details are needed to assess the reported 3.3× speedup.
  4. [§2, Eqs. (5)-(9) and Table 1] The notation μα and μβ is described as phase “chemical potentials,” although Eq. (5) uses them as molar Gibbs-energy densities before differentiation, while μ̃_B denotes a diffusion potential. Defining these quantities and their units explicitly would avoid ambiguity.
  5. [Reproducibility] Please deposit or otherwise make available the simulation inputs/scripts and identify database versions. The thermodynamic sources are cited, but the numerical examples are not currently independently reproducible from the manuscript alone.
  6. [Throughout] There are several typographical issues, including “has to be imposed,” “furthermore simplifications,” “wide ly,” and “yttrim-stabilized.” A careful copyedit is recommended.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: WBM/KKS reductions are algebraic limiting cases of a postulated free-energy model, not fitted or self-justifying predictions.

full rationale

This is a methods paper that postulates a free-energy functional and Allen–Cahn/Cahn–Hilliard kinetics in the standard phase-field style, then introduces an auxiliary non-conserved field φ for phase-composition differences. The claimed reductions to WBM (φ≡0, L_φ=0) and KKS (fast relaxation of φ to equal diffusion potentials) are direct algebraic/limiting identities of that construction, not empirical predictions forced by fitted inputs. Numerical demonstrations use external CALPHAD databases (Ghosh; Wang et al.) without tuning free parameters to manufacture the reported equilibria. Self-citations (authors’ prior SOFC/sublattice order-parameter work) appear only as background for possible multicomponent extensions and are not load-bearing for the binary two-phase derivation or the Al–Cu/Ti–V results. Concerns about the pragmatic deletion of the h(ξ)(1−h(ξ)) factor in Eq. (9) versus the variational Eq. (8) are justification/correctness issues, not circularity: nothing is redefined as its own input. The derivation chain is therefore self-contained against external benchmarks.

Assumptions & free parameters 3 free parameters · 5 assumptions · 1 invented entities

The claim rests on standard diffuse-interface thermodynamics plus one new modeling object (φ) and a pragmatic kinetic simplification (Eq. 9). Free parameters are the usual phase-field kinetic and interface coefficients plus the new L_φ; none are fitted to force the central equilibria. Invented entity is the auxiliary composition-difference field itself.

free parameters (3)
  • L_φ (relaxation coefficient for auxiliary composition difference) = case-dependent (0 … 100 L_ξ)
    Chosen by hand relative to L_ξ (e.g. 0 to 100 L_ξ in Ti–V; fixed value in Al–Cu) to interpolate between WBM and KKS; not derived from a microscopic mobility.
  • Interface parameters w, κ_ξ (and κ_X when used) = w=1.1e9 J/m³, κ_ξ=4.09e-9 J/m (examples)
    Set to target interfacial energy ~0.5 J/m² and chosen interface thickness (6 nm or 60 nm); standard phase-field calibration, not predicted.
  • L_ξ and M_X (interface and interdiffusion mobilities) = problem-specific constants listed in §3
    Assigned constant values for numerical examples; control growth kinetics but are not outputs of the interfacial-equilibrium theory.
assumptions (5)
  • domain assumption Non-conserved fields obey Allen–Cahn and conserved composition obeys Cahn–Hilliard dynamics derived from a free-energy functional.
    Standard phase-field postulate used for ξ, X_B, and φ (§2, Eqs. 7–9).
  • domain assumption At an interface the local composition is the interpolated mixture X_B = (1−h(ξ)) X_B^α + h(ξ) X_B^β with the usual quintic h.
    Shared with WBM/KKS; Eq. (1). Load-bearing for reconstructing phase compositions from (X_B, φ).
  • domain assumption Bulk free-energy density is the same interpolated mixture of phase chemical potentials μ_α(X_B^α), μ_β(X_B^β).
    Eq. (5); standard multiphase phase-field construction.
  • ad hoc to paper Removing the h(ξ)(1−h(ξ)) prefactor from the variational φ equation yields an acceptable evolution law (Eq. 9) that equalizes diffusion potentials in bulk as well as at interfaces.
    Authors note Eq. (8) is thermodynamically consistent but ineffective in bulk; Eq. (9) is proposed for practicality (§2). Central to matching KKS globally.
  • domain assumption Phase chemical potentials from CALPHAD databases may be used pointwise without parabola approximation or common-composition extrapolation.
    Stated goal and practice in §3.1–3.2; relies on databases remaining numerically well-behaved inside each phase’s allowed range.
invented entities (1)
  • Auxiliary non-conserved composition-difference field φ = X_B^α − X_B^β (and its multicomponent generalizations)
    purpose: Closes the under-determined interface mixture rule and provides a continuous kinetic path between equal-composition and equal-diffusion-potential limits.
    Defined in Eqs. (2)–(3); evolved by Eqs. (8) or (9). Not a physical order parameter measured independently; a modeling device.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A General Model of Interfacial Chemical Equilibrium in Phase-Field Method." pith.science (2026). https://pith.science/paper/O2MH6UON

@misc{pith2026260724620,
  author       = {Pith},
  title        = {Pith review of: A General Model of Interfacial Chemical Equilibrium in Phase-Field Method},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/O2MH6UON}},
  note         = {Machine review of arXiv:2607.24620}
}
read the original abstract

We report a new approach to interfacial equilibria among multiple solution phases in the phase-field method. It employs auxiliary non-conserved variables to describe the composition differences among different phases, and their temporal evolution is driven by the differences among different chemical diffusional potentials. It is reduced to the Wheeler-Boettinger-McFadden (WBM) model of equal interfacial chemical compositions and to the Kim-Kim-Suzuki (KKS) model of equal interfacial chemical diffusion potentials as two limiting cases. It can directly incorporate thermodynamic databases without any further approximations and simplifications. It is generally applicable to a wide range of problems involving composition evolution and chemical equilibria, including processes such as interdiffusion between two ordered phases not sharing any common composition range, which would pose difficulty to the existing treatments using WBM and KKS models.

Figures

Figures reproduced from arXiv: 2607.24620 by the authors.

Figure 1
Figure 1. Illustration of (a) the WBM model, the KKS model and (b) the current phase-field model in terms of extra potential 𝜇𝜇𝑒𝑒𝑒𝑒 at the interface. The magenta and blue lines represent the phase chemical potentials μα and μβ as a function of composition, respectively; the black dashed lines represent the common tangent line. 3. Results and Discussions In this section, we present two application examples of our proposed mode… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references

  1. [2]

    J., Warren, J

    Boettinger, W. J., Warren, J. A., Beckermann, C. & Karma, A. Phase -field simulation of solidification. Annu. Rev. Mater. Sci. 32, 163–194 (2002)

    2002

  2. [3]

    Phase-Field Model for Microstructure Evolution at the Mesoscopic Scale

    Steinbach, I. Phase-Field Model for Microstructure Evolution at the Mesoscopic Scale. Annu. Rev. Mater. Res. 43, 89–107 (2013)

    2013

  3. [4]

    & LLorca, J

    Tourret, D., Liu, H. & LLorca, J. Phase -field modeling of microstructure evolution: Recent applications, perspectives and challenges. Prog. Mater. Sci. 123, 100810 (2022)

    2022

  4. [5]

    & Zhao, Y

    Chen, L.-Q. & Zhao, Y. From classical thermodynamics to phase-field method. Prog. Mater. Sci. 124, 100868 (2022)

    2022

  5. [6]

    I., Warren, J

    Gránásy, L., Tóth, G. I., Warren, J. A., Podmaniczky, F., Tegze, G., Rátkai, L. & Pusztai, T. Phase-field modeling of crystal nucleation in undercooled liquids – A review. Prog. Mater. Sci. 106, 100569 (2019). 21

    2019

  6. [7]

    Wang, Y. & Li, J. Phase field modeling of defects and deformation. Acta Mater. 58, 1212– 1235 (2010)

    2010

  7. [8]

    & Wollants, P

    Moelans, N., Blanpain, B. & Wollants, P. An introduction to phase -field modeling of microstructure evolution. Calphad Comput. Coupling Phase Diagrams Thermochem. 32, 268– 294 (2008)

    2008

  8. [9]

    A., Boettinger, W

    Wheeler, A. A., Boettinger, W. J. & McFadden, G. B. Phase-field model for isothermal phase transitions in binary alloys. Phys. Rev. A 45, 7424–7439 (1992)

    1992

Show all 20 references
  1. [10]

    G., Kim, W

    Kim, S. G., Kim, W. T. & Suzuki, T. Phase-field model for binary alloys. Phys. Rev. E - Stat. Physics, Plasmas, Fluids, Relat. Interdiscip. Top. 60, 7186–7197 (1999)

    1999

  2. [11]

    L., Larsson, H

    Wang, W., Chen, H. L., Larsson, H. & Mao, H. Thermodynamic constitution of the Al –Cu– Ni system modeled by CALPHAD and ab initio methodology for designing high entropy alloys. Calphad Comput. Coupling Phase Diagrams Thermochem. 65, 346–369 (2019)

    2019

  3. [12]

    Y., Murray, J., Weiland, H., Liu, Z

    Hu, S. Y., Murray, J., Weiland, H., Liu, Z. K. & Chen, L. Q. Thermodynamic description and growth kinetics of stoichiometric precipitates in the phase- field approach. Calphad 31, 303– 312 (2007)

    2007

  4. [13]

    D., Zhang, F., Simpson, T

    Ji, Y., Lou, Y., Qu, M., Rowatt, J. D., Zhang, F., Simpson, T. W. & Chen, L. Q. Predicting Coherency Loss of γ″ Precipitates in IN718 Superalloy. Metall. Mater. Trans. A Phys. Metall. Mater. Sci. 47, 3235–3247 (2016)

    2016

  5. [14]

    Allen, S. M. & Cahn, J. W. A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening. Acta Metall. 27, 1085–1095 (1979)

    1979

  6. [15]

    Cahn, J. W. & Hilliard, J. E. Free Energy of a Nonuniform System. I. Interfacial Free Energy. Chem. Phys. 258, (1958)

    1958

  7. [16]

    Thermodynamic and kinetic modeling of the Cr-Ti-V system

    Ghosh, G. Thermodynamic and kinetic modeling of the Cr-Ti-V system. J. Phase Equilibria 23, 310–328 (2002)

    2002

  8. [17]

    & Qiu, C

    Hillert, M. & Qiu, C. A reassessment of the Cr-Fe-Ni system. Metall. Trans. A 21, 1673–1680 (1990)

    1990

  9. [18]

    N., Hallstedt, B

    Grundy, A. N., Hallstedt, B. & Gauckler, L. J. Assessment of the La -Sr-Mn-O system. Calphad Comput. Coupling Phase Diagrams Thermochem. 28, 191–201 (2004)

    2004

  10. [19]

    Ji, Y., Abernathy, H. W. & Chen, L.- Q. Thermodynamic Models of Multicomponent Nonstoichiometric Solution Phases Using Internal Process Order Parameters. Acta Mater. 117462 (2021)

    2021

  11. [20]

    N., Poplawsky, J

    Picard, Y. N., Poplawsky, J. D., Lee, S. & Abernathy, H. W. Nanoscale Analysis of LSM/YSZ Interfaces within Composite Cathodes for Commercial Solid Oxide Fuel Cells. ECS Meet. Abstr. MA2021-03, 54–54 (2021)

    2021

  12. [21]

    Thermodynamic Investigation of La0.8Sr0.2MnO3±δ Cathode, including the Prediction of Defect Chemistry, Electrical Conductivity and Thermo- Mechanical Properties

    Darvish, S. Thermodynamic Investigation of La0.8Sr0.2MnO3±δ Cathode, including the Prediction of Defect Chemistry, Electrical Conductivity and Thermo- Mechanical Properties. (Florida International University, 2018). 22 Supplementary Materials S1. About the evolution equation o...

    2018

Pith tools

Reviewed July 31, 2026 · model on record in the stance chip above.