{"id":"f5f5548c-3409-49d4-b8a2-86934d61c2c9","arxiv_id":"2411.16430","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A single variational framework reproduces both Cahn-Hilliard phase separation and CALPHAD-style sharp-interface diffusion as limiting cases, with a stabilized mixed finite element method.","lead":"The paper introduces a unified simulation framework that treats two standard ways of modeling phase changes, Cahn-Hilliard and CALPHAD-based reactive diffusion, as two ends of one variational model controlled by interface energy. It adds a stabilized finite element scheme that can handle both smooth and sharp interfaces, and demonstrates the approach on binary alloys and a vacancy-driven ternary alloy.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Continuous P2 mixed elements in Eqs. (49)-(50) do not inherit the stability of [BDM85]; the central claim of a stable continuous FEM for sharp interfaces is unsupported and likely false.","rationale":"I read the paper as making two intertwined claims: (1) reactive diffusion and Cahn-Hilliard models are unified by a single Lagrangian with flux as an independent variable, and (2) this formulation yields a stable continuous finite element method for both smooth and sharp interfaces. The first claim is plausible and largely follows from standard Onsager variational arguments, modulo the sign inconsistency between Eqs. (26) and (46) and the dropped time-integration constant in Eq. (31). The second claim is the load-bearing pillar: without a stable discretization, the numerical demonstrations in Section 3 do not support the methodology. The paper explicitly relies on [BDM85] for stability, yet it replaces the H(div)-conforming BDM spaces with continuous P2 spaces. That substitution is known to destroy the inf-sup condition for mixed Laplace-type problems; equal-order continuous P2/P2 for flux and Lagrange multiplier is not a stable pair. The transient term and the nonlinear relation µ̄ = -∂f/∂x do not repair the inf-sup deficiency because the affinity space remains too rich relative to the divergence of the flux space. Thus the central claim 'stable continuous FEM' is not established and is likely false as written. This is not a disagreement with consensus; it is an internal mismatch between the cited theory and the chosen finite element spaces. The reader's weakest assumption, Eq. (39), is an acknowledged approximation that affects quantitative fidelity but does not threaten the mathematical core as directly as the stability failure. I therefore recommend REJECT: the main numerical contribution needs either a different, inf-sup-stable element choice or a genuine stability analysis, not just a citation to BDM mixed elements.","tokens_in":15537,"tokens_out":13435,"duration_ms":137583,"concrete_test":"Assemble the linearized steady mixed system from Eqs. (48)-(50) with a quadratic free energy (so µ̄ = -a x) on a unit-square mesh and refine h = 1/8, 1/16, 1/32, 1/64. Compute the discrete inf-sup constant β_h = inf_{q_h ∈ P2} sup_{v_h ∈ [P2]^2} (∫ q_h ∇·v_h) / (‖v_h‖_{H(div)} ‖q_h‖_{L2}). If β_h decays with h instead of remaining bounded below, equal-order continuous P2/P2 is not inf-sup stable and the BDM-based stability claim fails. A complementary check is to solve the steady mixed problem with a smooth manufactured solution and monitor L2 errors; if checkerboard oscillations appear or convergence rates degrade, the instability is confirmed.","verdict_should_be":"REJECT","load_bearing_attack":"The paper's core numerical contribution is the claim that the mixed variational forms (49) and (50) stabilize the solution for convex free energies, citing [BDM85]. But [BDM85] establishes stability for H(div)-conforming mixed finite elements with appropriately chosen pressure spaces, not for equal-order continuous P2 spaces for both the flux j and the affinity µ̄ as used in Section 2.4. For the resulting mixed Poisson saddle-point problem, the discrete inf-sup condition for the bilinear form b(j, µ̄) = ∫ µ̄ ∇·j is not satisfied by continuous P2/P2: the divergence of P2 vector fields lies in P1, while the affinity space is P2, allowing spurious checkerboard modes. The nonlinear coupling in Eq. (48), which ties µ̄ to the composition x, does not restore a uniform inf-sup constant; it merely slaves the unstable pressure modes to x, so the discrete Laplacian B M_j^{-1} B^T retains a near-kernel. Thus the method's advertised stability for sharp interfaces is not demonstrated, and the oscillations visible in Figure 3 may be the expected numerical artifact rather than physical behavior. Separately, the discrete Lagrangian (46) changes the sign of the constraint term relative to the continuous Lagrangian (26), so Eqs. (48)-(50) do not exactly discretize the claimed variational principle, but the inf-sup issue is the more fundamental threat to the central claim.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a unified variational framework for diffusional phase transformations, claiming that Cahn-Hilliard models and reactive diffusion (CALPHAD-based) models are the same physical model distinguished only by the choice of molar free energy fm. The authors derive a Lagrangian formulation with fluxes and chemical affinities as independent variables, prove thermodynamic consistency via nonnegative entropy production and a Gibbs-Duhem check, and introduce a relation between Onsager coefficients and interdiffusion coefficients. They then propose a continuous Galerkin mixed finite element discretization of degree two for the resulting variational forms and demonstrate the method on binary Cahn-Hilliard, sharp-interface, and ternary vacancy-diffusion examples. The central claims are the conceptual unification of the two modeling paradigms and the existence of a stable continuous finite element method for both smooth and sharp interface problems.","tokens_in":107,"tokens_out":10656,"duration_ms":360775,"significance":"If the claims hold, the paper would make a useful conceptual contribution by showing that Cahn-Hilliard and reactive diffusion models are not fundamentally different, and by providing a single variational principle from which both can be recovered. The algebraic derivation of the Gibbs-Duhem consistency and the explicit link between Onsager and interdiffusion coefficients are valuable and clearly presented. The numerical stabilization claim, however, is the main advertised novelty and is currently not supported by the analysis or the experiments; the manuscript needs substantial revision before the stability and consistency of the discrete formulation can be accepted.","major_comments":[{"comment":"The stabilization claim is not established. The paper cites [BDM85] as the origin of forms (49)-(50), but [BDM85] proves stability for H(div)-conforming mixed finite elements with appropriately selected pressure spaces, not for continuous P2/P2 equal-order elements as used here. For continuous P2 vector fields and P2 scalar affinities, the discrete inf-sup condition for b(j, \\bar{\\mu}) = \\int \\bar{\\mu} \\nabla\\cdot j is not satisfied: the divergence of a P2 vector field lies in P1, so any P2 affinity component orthogonal to P1 belongs to the discrete kernel, giving an inf-sup constant of zero and permitting checkerboard modes. The oscillations visible in Figure 3 may therefore be the expected numerical artifact of this unstable discretization rather than physical interfacial fluctuations. This directly undermines the advertised stable continuous FEM for sharp interfaces.","section":"Section 2.4, Eqs. (49)-(50)"},{"comment":"There is a sign inconsistency between the continuous and discrete Lagrangians. The continuous Lagrangian (26)/(40) contains the constraint term with a minus sign, while the scaled discrete Lagrangian (45)/(46) uses a plus sign. As written, Eq. (48) enforces \\bar{\\mu} = -\\delta F/\\delta x rather than Eq. (31)'s \\bar{\\mu} = \\delta f_m/\\delta x, and Eq. (49), after integration by parts, corresponds to j = +(D/k)\\nabla\\bar{\\mu} rather than Eq. (42)'s j = -(D/(k\\Omega))\\nabla\\bar{\\mu}. The two sign changes happen to cancel if \\bar{\\mu} is reinterpreted as the negative affinity, but the manuscript does not make this reinterpretation, and the variational forms are therefore not a discretization of the stated variational principle. The authors must correct the signs or redefine \\bar{\\mu} consistently, then rerun the simulations.","section":"Section 2.4, Eqs. (45)-(49) and Eq. (42)"},{"comment":"The relation \\tilde{L} = \\tilde{D} T/(k\\Omega) is derived under the assumption that the molar free energy is quadratic in the bulk phases, and the paper explicitly states that the relation fails at interfaces. Nevertheless, the scaled Lagrangian (41) and all numerical experiments apply this relation globally. Since the paper's new capability is specifically the treatment of sharp interfaces, the interface region—where this relation is invalid—may contribute appreciably to the kinetics. The manuscript needs a quantitative justification, such as a convergence study in interface width or an estimate of the interface contribution, before the kinetic predictions built on Eq. (39) can be considered reliable.","section":"Section 2.2, Eq. (39), used in Section 2.3 and all simulations"}],"minor_comments":[{"comment":"The step from \\partial_t(\\delta f_m/\\delta x - \\bar{\\mu})=0 to \\bar{\\mu}=\\delta f_m/\\delta x silently drops a time-independent integration constant; the constant should be fixed by initial or boundary conditions, or the argument should state why it is irrelevant for the subsequent flux law.","section":"Eqs. (30)-(31)"},{"comment":"The caption states that the solution \"just behaves well when the phases are entirely separated,\" which directly contradicts the claimed stabilization of the sharp-interface case; the authors should clarify whether the observed fluctuations are numerical artifacts or physical features.","section":"Figure 3 caption"},{"comment":"There are numerous typographical errors, including \"adressed\" (p. 3), \"Therby\" (p. 9), \"subsititutional\" (p. 11), \"entirely seperated\" (Fig. 7 caption), and \"Therefore\" capitalized mid-sentence in Section 2.3.","section":"Throughout"},{"comment":"The keyword \"Theorem of Minimum Entropy Production\" is listed, but the paper does not state or prove such a theorem; either add a relevant discussion or remove the keyword.","section":"Keywords"},{"comment":"The use of \\tilde{D}(x_t) in Eq. (51) is presented as a way to avoid variation of the dissipative term, but the paper does not discuss the consistency of this explicit-in-time treatment with the variational principle; a brief comment on the resulting time-discretization error would be helpful.","section":"Section 2.5"}],"recommendation":"major_revision","confidential_remarks":"The conceptual unification and thermodynamic derivations are interesting, but the numerical stability claim and the sign consistency of the discrete variational principle are load-bearing and currently unresolved. The issues appear fixable within the manuscript's scope, but the authors should be asked to provide a rigorous discrete stability analysis (or use genuinely stable elements) and to correct the sign conventions before the paper can be considered further."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here is my take on 2411.16430. The paper is genuinely useful as a conceptual bridge: it shows Cahn-Hilliard and CALPHAD-based reactive diffusion can be derived from one variational Lagrangian, with the choice of molar free energy as the knob that moves between smooth and sharp interfaces. The thermodynamic consistency argument is standard but correctly assembled, and the authors are honest about the bulk-phase-only status of the Onsager-Fick relation (Eq. 39). The vacancy example is a good stress test for the framework. I agree with the reader that the central derivation is sound in spirit.\n\nThe problems are concentrated in the discrete formulation. First, the discrete Lagrangian (46) has a plus sign on the constraint term while the continuous one (26) has a minus sign, so the discrete variational forms (48)-(50) do not actually discretize the claimed principle with the same definition of μ̄. In particular, (48) forces μ̄ = -δfm/δx, not +δfm/δx. This may be fixable by redefining μ̄, but as written it is an inconsistency that has to be resolved.\n\nSecond, and more serious: the paper claims stabilization by citing [BDM85], but then uses continuous P2 elements for both flux and affinity. That is not the BDM mixed space, and the inf-sup condition for the mixed Poisson structure is not satisfied for P2/P2. The oscillations visible in Figure 3 are exactly what a near-kernel of the discrete Laplacian produces. The authors even write that the method is 'at the limits of numerical stability'—that is not a stable method. The central numerical claim of the paper is therefore unsupported and likely false. Switching to a stable element pair (e.g., H(div) flux with a discontinuous pressure-like multiplier, or a P2/P1-like choice) or adding stabilized terms is necessary.\n\nThere are also smaller issues: Eq. (31) drops a time-integration constant without comment; no code or data are provided; and the numerical validation is qualitative only. The Onsager-Fick relation is applied globally, but since the authors flag that it fails at interfaces, I treat that as a limitation rather than a fatal flaw.\n\nWho is this for? Researchers who want to see the two modeling communities connected, and who can separate the conceptual message from the numerics. The conceptual part deserves referee time; the numerics need major revision before the paper can be accepted. I would send it to peer review rather than desk-reject, because the underlying idea is worthwhile and the flaws are identifiable and fixable.","headline":"A genuinely useful conceptual unification of Cahn-Hilliard and reactive diffusion, but the advertised stable continuous FEM is likely unstable with P2/P2 and the discrete variational forms contain a sign inconsistency.","tokens_in":16383,"tokens_out":8735,"would_cite":false,"duration_ms":80120,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M60","35K55","80A22"],"pacs":[],"model":"deepseek-v4-flash","headline":"Reactive diffusion and the Cahn-Hilliard model are the same variational problem, differing only in the molar free energy.","keywords":["Cahn-Hilliard equation","reactive diffusion","CALPHAD","interface energy","Onsager coefficients","mixed finite element method","phase-field modeling","vacancy diffusion"],"falsifier":"Run the sharp-interface model on a domain with many closely spaced interfaces, so that interface regions occupy a measurable fraction of the volume, and compare the resulting phase-growth kinetics against the same initial state modeled with Onsager coefficients evaluated separately inside each bulk phase; any significant discrepancy in interface velocity or concentration profiles would falsify the global use of Equation (39).","tokens_in":15299,"feed_emoji":"⚛️","tokens_out":9458,"duration_ms":85089,"temperature":0.7,"pith_summary":"This paper argues that two modeling traditions that materials scientists usually treat as separate, Cahn-Hilliard phase-field models and CALPHAD-based reactive diffusion models, are in fact one and the same mathematical object once both are written as a variational problem with the diffusional flux as an independent variable. The two differ only in what molar free energy $f_m$ is inserted: a double-well potential with a gradient term for Cahn-Hilliard, or its convex hull (the CALPHAD common-tangent energy) for reactive diffusion. On that basis the authors build a continuous finite element method, stabilized by solving explicitly for the flux and the chemical affinities, that remains stable for both smooth interfaces and sharp, regularization-free interfaces. If the unification holds, a single simulation framework can dial interface energy from zero to arbitrary strength and can treat phase transformations, ternary alloys, and vacancy-mediated diffusion without switching models.","feed_headline":"One variational principle unifies Cahn-Hilliard and CALPHAD models","feed_subtitle":"A single Lagrangian treats smooth and sharp interfaces, putting interface energy under the modeler's control.","key_machinery":"The load-bearing object is the Lagrangian functional with the flux as an independent variable, together with the mixed variational form that solves for flux explicitly. Because the affinity never appears under a differential operator in these forms, the solution no longer requires differentiating a free energy that is only piecewise smooth or $C^1$ (the convex hull), which is what previously forced either regularization or discontinuous Galerkin discretizations. The relation $\\tilde{L}=\\tilde{D}T/(k\\Omega)$, obtained from in-phase quadratic free energies, is what converts measured diffusion coefficients into Onsager coefficients and produces the scaled Lagrangian in which the flux has units of velocity.","core_discovery":"The central claim is that reactive diffusion and the Cahn-Hilliard phase-field model are inherently identical, distinguishable only by their definitions of the molar free energy $f_m$. The proof runs through a Lagrangian in which the fluxes are independent variables and the chemical affinities enter as Lagrange multipliers that enforce mass conservation; varying the Lagrangian with respect to flux, mole fraction, and affinity reproduces the flux law, the identity that the affinity is the variational derivative of the molar free energy, and the mass-conservation equations. Combined for a binary system this gives $\\partial c/\\partial t = \\nabla\\cdot(-\\tilde{L}/T\\,\\nabla(\\delta f_m/\\delta x))$, which is the Cahn-Hilliard equation when $f_m$ contains a double-well potential and the CALPHAD/reactive-diffusion equation when $f_m$ is its convex hull. The paper further proposes the relation $\\tilde{L} = \\tilde{D}T/(k\\Omega)$ between Onsager coefficients, interdiffusion coefficients, the curvature $k$ of the quadratic phase free energies, temperature, and molar volume, and shows how the resulting scaled Lagrangian supports a continuous Galerkin finite element method that stays stable for both smooth and sharp interfaces.","pith_inferences":["An untested corollary of the claimed equivalence is that every Cahn-Hilliard simulation with vanishing gradient coefficient has a sharp-interface reactive-diffusion twin with the same macroscopic kinetics; running a known spinodal-decomposition benchmark in both limits would test the unification directly.","Because Equation (39) is exact only inside bulk phases, the error of kinetic predictions should grow with the volume fraction of interfaces; for lamellar or nanoprecipitate microstructures with dense interfaces, the global relation may need to be replaced by phase-local Onsager coefficients.","The Lagrangian's structure of free energy plus dissipation plus constraint invites direct coupling to mechanics, heat conduction, and other rate processes by adding terms, so the method could become a template for multiphysics phase-transformation simulation.","The convex-hull/double-well duality is not alloy-specific; the same variational identification might apply to any system where a free energy and its convex envelope describe different interface physics."],"forward_implications":["A single finite-element code can reproduce classical Cahn-Hilliard spinodal decomposition, regularized smooth-interface transformations, and sharp-interface CALPHAD-type reactive diffusion by switching the molar free energy between double-well, convex-hull with $\\kappa>0$, and convex-hull with $\\kappa=0$.","Interface energy strength can be varied continuously in the same framework, so morphology changes, such as the observed transition from spherical coarsening to non-area-minimizing growth, can be studied as a function of $\\Delta F_{\\mathrm{int}}$ rather than as a change of model.","The thermodynamic-consistency proof and the explicit flux solve guarantee nonnegative entropy production and mass conservation, so results are not hostage to the stability of a particular discretization of the chemical potential.","The Onsager-to-diffusion-coefficient relation turns experimentally measured interdiffusion coefficients and phase-curvature data into model inputs, enabling ternary and multi-component simulations with non-diagonal, coupled fluxes.","The vacancy model with sources and sinks, including non-conserved vacancy generation and annihilation and the Kirkendall effect, is realized within the same variational formalism, showing that the method reaches beyond binary diffusion."],"supporting_citations":[{"why":"Defines the Cahn-Hilliard free-energy functional with the double-well potential and gradient-energy term that the paper compares against.","marker":"[CH58]"},{"why":"Supplies the linear nonequilibrium thermodynamics and reciprocal relations underlying the Onsager-coefficient flux law.","marker":"[Ons31]"},{"why":"Gives the CALPHAD-based reactive-diffusion treatment and its finite-difference implementation, the sharp-interface counterpart being unified with Cahn-Hilliard.","marker":"[SF13]"},{"why":"Provides the vacancy diffusion model with sources and sinks used for the ternary example and the Kirkendall-effect demonstration.","marker":"[SFF06]"},{"why":"Introduces the mixed finite-element spaces for second-order elliptic problems that the stabilized variational form draws on.","marker":"[BDM85]"},{"why":"Provides the chemical-potential formula in terms of molar free energy and the common-tangent (CALPHAD) construction of the convex hull.","marker":"[Lup83]"},{"why":"Supplies the standard form of the Onsager flux equations, the entropy production rate, and the chemical-affinity definitions used in the Lagrangian.","marker":"[KP14]"},{"why":"Previous numerical treatment of reactive diffusion by the discontinuous Galerkin method whose analysis of affinity smoothness motivates the new continuous scheme.","marker":"[Fla+24]"}],"fun_headline_variants":["Single Lagrangian merges phase-field and reactive diffusion","Interface energy now tunable in both phase-field and CALPHAD","Unified model bridges sharp and smooth interface methods","Cahn-Hilliard and CALPHAD under one variational roof"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the relation $\\tilde{L}=\\tilde{D}T/(k\\Omega)$, derived inside bulk phases where the molar free energy is quadratic, can be applied globally in the simulations, even though the paper itself notes it fails at interfaces where the convex free energy is linear.","fun_headline_variants_meta":{"raw":{"variants":["Single Lagrangian merges phase-field and reactive diffusion","Interface energy now tunable in both phase-field and CALPHAD","Unified model bridges sharp and smooth interface methods","Cahn-Hilliard and CALPHAD under one variational roof"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000166,"raw_usage":{"total_tokens":1275,"prompt_tokens":986,"completion_tokens":289,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":222}},"tokens_in":602,"tokens_out":289,"duration_ms":11958,"temperature":1.0,"reasoning_tokens":222,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:08:24.763051+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the sharp-interface model on a domain with many closely spaced interfaces, so that interface regions occupy a measurable fraction of the volume, and compare the resulting phase-growth kinetics against the same initial state modeled with Onsager coefficients evaluated separately inside each bulk phase; any significant discrepancy in interface velocity or concentration profiles would falsify the global use of Equation (39).","supporting_citations":[],"review_version":1}