pith. sign in

arxiv: 2412.06762 · v2 · pith:MR2LE5S4new · submitted 2024-12-09 · 🧮 math.AP

Interface dynamics in a degenerate Cahn-Hilliard model for viscoelastic phase separation

Katharina Hopf , John King , Andreas M\"unch , Barbara Wagner This is my paper

Pith reviewed 2026-05-23 07:34 UTC · model grok-4.3

classification 🧮 math.AP
keywords sharp-interface asymptoticsCahn-Hilliard modelviscoelastic phase separationsurface diffusion flowgradient flow structuredouble-obstacle potentialmatched asymptoticsnon-local evolution laws
0
0 comments X

The pith

Degenerate Cahn-Hilliard asymptotics produce non-local third-order interface flows

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

This paper performs formal matched asymptotics on a degenerate Cahn-Hilliard model with cross-diffusive coupling to stress for viscoelastic phase separation. It demonstrates that the sharp-interface limit yields non-local versions of surface diffusion flow, reducing to the intermediate surface diffusion flow when coupling is constant. For monotonic non-constant coupling the motion is governed by a new family of third-order laws derived from a constrained elliptic problem whose solution encodes the gradient structure. The result matters for understanding how viscoelastic effects alter phase separation dynamics at the interface scale.

Core claim

In the deep quench regime the formal sharp-interface asymptotics of the model with rank-deficient mobility lead to non-local lower-order counterparts of classical surface diffusion flow. Constant coupling gives the intermediate surface diffusion flow. Non-constant monotonic coupling produces a new family of third-order evolution laws in which the normal velocity is the Lagrange multiplier of a constrained elliptic equation on the interface; this equation is solved rigorously by a variational argument and encodes the gradient structure of the effective law.

What carries the argument

Constrained elliptic equation on the interface whose Lagrange multiplier is the normal velocity and which encodes the gradient structure of the third-order flow

If this is right

  • The effective interface evolution inherits the Onsager gradient-flow structure of the diffuse model.
  • The leading-order operator for non-constant coupling behaves like the square root of the negative Laplace-Beltrami operator.
  • The analysis is carried out for the double-obstacle potential in the deep-quench limit.
  • Constant coupling recovers the known intermediate surface diffusion flow.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach may extend to other degenerate phase-field models with similar mobility structures.
  • The variational elliptic constraint offers a route to stable numerical approximations of the interface motion.
  • These non-local flows could appear in related models of phase separation in soft matter with internal variables.

Load-bearing premise

The coupling functions are monotonic and connect the two phases and the formal matched asymptotics hold without higher-order corrections that alter the leading non-local operator.

What would settle it

A simulation of the original model in which the interface velocity does not match the normal velocity obtained from the constrained elliptic problem would falsify the asymptotic claim.

Figures

Figures reproduced from arXiv: 2412.06762 by Andreas M\"unch, Barbara Wagner, John King, Katharina Hopf.

Figure 1
Figure 1. Figure 1: Free transition layer ΩI ε (t) and sharp interface Γ(t). 11 [PITH_FULL_IMAGE:figures/full_fig_p011_1.png] view at source ↗
read the original abstract

The formal sharp-interface asymptotics in a degenerate Cahn-Hilliard model for viscoelastic phase separation with cross-diffusive coupling to a bulk stress variable are shown to lead to non-local lower-order counterparts of the classical surface diffusion flow. The diffuse-interface model is a variant of the Zhou-Zhang-E model and has an Onsager gradient-flow structure with a rank-deficient mobility matrix reflecting the ODE character of stress relaxation. In the case of constant coupling, we find that the evolution of the zero level set of the order parameter approximates the so-called intermediate surface diffusion flow. For non-constant coupling functions monotonically connecting the two phases, our asymptotic analysis leads to a new family of third-order evolution laws with associated propagation operators behaving, at leading order, like the square root of the minus Laplace-Beltrami operator. In this case, the normal velocity of the moving sharp interface arises as the Lagrange multiplier in a constrained elliptic equation, which is at the core of our derivation. The constrained elliptic problem can be solved rigorously by a variational argument, and is shown to encode the gradient structure of the effective geometric evolution law. The asymptotics are presented for deep quench, an intermediate free boundary problem based on the double-obstacle potential.

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 / 0 minor

Summary. The paper performs formal sharp-interface asymptotics on a degenerate Cahn-Hilliard model (variant of Zhou-Zhang-E) with rank-deficient mobility and cross-diffusive coupling to a bulk stress variable. For constant coupling the zero level set evolves by the intermediate surface diffusion flow; for non-constant monotonic coupling functions the analysis yields a new family of third-order non-local geometric laws whose normal velocity is the Lagrange multiplier of a constrained elliptic problem solved variationally, with the propagation operator behaving like the square root of the minus Laplace-Beltrami operator. The asymptotics are carried out in the deep-quench regime with the double-obstacle potential.

Significance. If the formal derivation is accurate, the manuscript supplies a new class of non-local interface evolution equations that inherit a gradient-flow structure from the underlying PDE and extend classical surface diffusion. The rigorous variational treatment of the auxiliary constrained elliptic problem is a clear technical strength.

major comments (2)
  1. [Abstract and matched asymptotics] Abstract (non-constant coupling paragraph) and the matched-asymptotics derivation: the claim that the leading-order inner problem is exactly the stated constrained elliptic equation whose multiplier is the normal velocity rests on the assumption that O(ε) curvature and outer-variation contributions from the cross-diffusive term do not enter the solvability condition at the same order. No explicit expansion or estimate is supplied showing that these corrections remain higher order for the rank-deficient mobility and singular inner ODEs at the interface endpoints.
  2. [Abstract] Abstract (final sentence on variational argument): while the constrained elliptic problem is solved rigorously, the manuscript provides neither error estimates for the formal asymptotics nor any numerical verification that the derived non-local operator is recovered in the sharp-interface limit; the central claim therefore depends on unverified formal steps.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for highlighting the potential significance of the derived non-local interface laws. We address the two major comments below. The derivation remains formal, consistent with the manuscript's scope.

read point-by-point responses

  1. Referee: [Abstract and matched asymptotics] Abstract (non-constant coupling paragraph) and the matched-asymptotics derivation: the claim that the leading-order inner problem is exactly the stated constrained elliptic equation whose multiplier is the normal velocity rests on the assumption that O(ε) curvature and outer-variation contributions from the cross-diffusive term do not enter the solvability condition at the same order. No explicit expansion or estimate is supplied showing that these corrections remain higher order for the rank-deficient mobility and singular inner ODEs at the interface endpoints.

    Authors: We agree that the manuscript would benefit from a more explicit order-by-order expansion to confirm that O(ε) curvature and outer cross-diffusive contributions remain higher order. The rank-deficient mobility and the singular behavior of the inner ODEs at the endpoints are central to why these terms do not affect the leading solvability condition, but this justification is only sketched. In the revision we will insert a dedicated subsection detailing the full inner expansion up to O(ε) and verifying the higher-order status of the indicated corrections. revision: yes

  2. Referee: [Abstract] Abstract (final sentence on variational argument): while the constrained elliptic problem is solved rigorously, the manuscript provides neither error estimates for the formal asymptotics nor any numerical verification that the derived non-local operator is recovered in the sharp-interface limit; the central claim therefore depends on unverified formal steps.

    Authors: The asymptotics are explicitly formal, as stated in the title and throughout the text; the manuscript does not claim or attempt a rigorous convergence proof. The variational treatment of the auxiliary constrained elliptic problem is rigorous and establishes the gradient-flow structure of the effective law, but error estimates between the diffuse-interface model and the sharp-interface limit lie outside the present scope. Numerical verification of the non-local operator is likewise not included. We therefore do not plan to add either error estimates or simulations, as both would require substantial additional work beyond the formal derivation. revision: no

Circularity Check

0 steps flagged

Formal matched asymptotics derivation is self-contained from the PDE model with no reduction to inputs by construction

full rationale

The paper derives the sharp-interface limit via formal matched asymptotics applied to the given degenerate Cahn-Hilliard model with rank-deficient mobility and double-obstacle potential. The key step identifies the normal velocity as the Lagrange multiplier in a constrained elliptic problem for non-constant monotonic coupling, solved variationally to encode the gradient structure. This follows directly from the inner/outer expansions and solvability conditions without fitting parameters, renaming known results, or load-bearing self-citations. The derivation chain begins from the Onsager structure of the original PDE and produces the non-local operators as output; no equation reduces to an input by definition or statistical forcing. This is the standard case of an independent asymptotic analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the model possessing an Onsager gradient-flow structure with rank-deficient mobility, monotonic coupling functions, and the validity of formal matched asymptotics in the deep-quench double-obstacle setting; no free parameters or invented entities are mentioned.

axioms (2)
  • domain assumption The system admits an Onsager gradient-flow structure with rank-deficient mobility matrix reflecting ODE character of stress relaxation
    Stated directly in the abstract as the structural property of the Zhou-Zhang-E variant.
  • domain assumption Coupling functions are monotonic and connect the two phases
    Required for the non-constant case leading to the new third-order laws.

pith-pipeline@v0.9.0 · 5750 in / 1475 out tokens · 22162 ms · 2026-05-23T07:34:10.006793+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

30 extracted references · 30 canonical work pages

  1. [1]

    N. D. Alikakos, P. W. Bates, and X. Chen. Convergence of the C ahn- H illiard equation to the H ele- S haw model. Arch. Rational Mech. Anal. , 128(2):165--205, 1994

    work page 1994

  2. [2]

    Abels, H

    H. Abels, H. Garcke, and G. Gr\" u n. Thermodynamically consistent, frame indifferent diffuse interface models for incompressible two-phase flows with different densities. Math. Models Methods Appl. Sci. , 22(3):1150013, 40, 2012

    work page 2012

  3. [3]

    A. D. Aleksandrov. Uniqueness theorems for surfaces in the large. Vestnik Leningrad. Univ. , 11(19):5--17, 1956

    work page 1956

  4. [4]

    Ambrosio

    L. Ambrosio. Geometric evolution problems, distance function and viscosity solutions. In Calculus of variations and partial differential equations ( P isa, 1996) , pages 5--93. Springer, Berlin, 2000

    work page 1996

  5. [5]

    J. F. Blowey and C. M. Elliott. The C ahn- H illiard gradient theory for phase separation with nonsmooth free energy. I . M athematical analysis. European J. Appl. Math. , 2(3):233--280, 1991

    work page 1991

  6. [6]

    Brunk and M

    A. Brunk and M. Luk\'a c ov\'a-Medvid'ov \'a . Global existence of weak solutions to viscoelastic phase separation: part II . D egenerate case. Nonlinearity , 35(7):3459--3486, 2022

    work page 2022

  7. [7]

    Bretin, S

    E. Bretin, S. Masnou, A. Sengers, and G. Terii. Approximation of surface diffusion flow: a second-order variational C ahn- H illiard model with degenerate mobilities. Math. Models Methods Appl. Sci. , 32(4):793--829, 2022

    work page 2022

  8. [8]

    J. W. Cahn, C. M. Elliott, and A. Novick-Cohen. The C ahn- H illiard equation with a concentration dependent mobility: motion by minus the L aplacian of the mean curvature. European J. Appl. Math. , 7(3):287--301, 1996

    work page 1996

  9. [9]

    J. W. Cahn and J. E. Taylor. Surface motion by surface diffusion. Acta Met. Mat. , 42(4):1045--1063, 1994

    work page 1994

  10. [10]

    Doi and A

    M. Doi and A. Onuki. Dynamic coupling between stress and composition in polymer solutions and blends. Journal de Physique II , 2(8):1631--1656, 1992

    work page 1992

  11. [11]

    C. M. Elliott and H. Garcke. On the C ahn- H illiard equation with degenerate mobility. SIAM J. Math. Anal. , 27(2):404--423, 1996

    work page 1996

  12. [12]

    C. M. Elliott and H. Garcke. Existence results for diffusive surface motion laws. Adv. Math. Sci. Appl. , 7(1):467--490, 1997

    work page 1997

  13. [13]

    Escher, Y

    J. Escher, Y. Giga, and K. Ito. On a limiting motion and self-intersections of curves moved by the intermediate surface diffusion flow. In Proceedings of the T hird W orld C ongress of N onlinear A nalysts, P art 6 ( C atania, 2000) , volume 47, pages 3717--3728, 2001

    work page 2000

  14. [14]

    Escher, Y

    J. Escher, Y. Giga, and K. Ito. On a limiting motion and self-intersections for the intermediate surface diffusion flow. J. Evol. Equ. , 2(3):349--364, 2002

    work page 2002

  15. [15]

    Escher and K

    J. Escher and K. Ito. On the intermediate surface diffusion flow. In Free boundary problems ( T rento, 2002) , volume 147 of Internat. Ser. Numer. Math. , pages 131--138. Birkh\" a user, Basel, 2004

    work page 2002

  16. [16]

    Escher, U

    J. Escher, U. F. Mayer, and G. Simonett. The surface diffusion flow for immersed hypersurfaces. SIAM J. Math. Anal. , 29(6):1419--1433, 1998

    work page 1998

  17. [17]

    Escher and G

    J. Escher and G. Simonett. Moving surfaces and abstract parabolic evolution equations , pages 183--212. Birkh \"a user, Basel, 1999

    work page 1999

  18. [18]

    Gugenberger, R

    C. Gugenberger, R. Spatschek, and K. Kassner. Comparison of phase-field models for surface diffusion. Phys. Rev. E (3) , 78(1):016703, 17, 2008

    work page 2008

  19. [19]

    Gilbarg and N

    D. Gilbarg and N. S. Trudinger. Elliptic partial differential equations of second order . Classics in Mathematics. Springer-Verlag, Berlin, 2001

    work page 2001

  20. [20]

    Luckhaus and L

    S. Luckhaus and L. Modica. The G ibbs- T hompson relation within the gradient theory of phase transitions. Arch. Rational Mech. Anal. , 107(1):71--83, 1989

    work page 1989

  21. [21]

    unch, and E. S\

    A. A. Lee, A. M\"unch, and E. S\"uli. Sharp-interface limits of the C ahn- H illiard equation with degenerate mobility. SIAM J. Appl. Math. , 76(2):433--456, 2016

    work page 2016

  22. [22]

    R. L. Pego. Front migration in the nonlinear C ahn- H illiard equation. Proc. Roy. Soc. London Ser. A , 422(1863):261--278, 1989

    work page 1989

  23. [23]

    u ss and G. Simonett. Moving interfaces and quasilinear parabolic evolution equations , volume 105 of Monographs in Mathematics . Birkh\

    J. Pr\" u ss and G. Simonett. Moving interfaces and quasilinear parabolic evolution equations , volume 105 of Monographs in Mathematics . Birkh\" a user/Springer, 2016

    work page 2016

  24. [24]

    P. J. Strasser, G. Tierra, B. D\"unweg, and M. Luk\'a c ov\'a-Medvid'ov \'a . Energy-stable linear schemes for polymer-solvent phase field models. Comput. Math. Appl. , 77(1):125--143, 2019

    work page 2019

  25. [25]

    H. Tanaka. Viscoelastic phase separation. Journal of Physics: Condensed Matter , 12(15):R207, 2000

    work page 2000

  26. [26]

    H. Tanaka. Viscoelastic phase separation in biological cells. Communications Physics , 5(1):167, 2022

    work page 2022

  27. [27]

    J. E. Taylor and J. W. Cahn. Linking anisotropic sharp and diffuse surface motion laws via gradient flows. J. Statist. Phys. , 77(1-2):183--197, 1994

    work page 1994

  28. [28]

    Taniguchi and A

    T. Taniguchi and A. Onuki. Network domain structure in viscoelastic phase separation. Phys. Rev. Lett. , 77:4910--4913, Dec 1996

    work page 1996

  29. [29]

    E. Zeidler. Nonlinear functional analysis and its applications. III . Springer-Verlag, New York, 1985. Variational methods and optimization

    work page 1985

  30. [30]

    D. Zhou, P. Zhang, and W. E. Modified models of polymer phase separation. Phys. Rev. E , 73(6):061801, 2006

    work page 2006