Pith. sign in

REVIEW 3 major objections 5 minor 1 cited by

Crystalline elastic flow of polygonal curves: long time behaviour and convergence to stationary solutions

T0 review · 3 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash

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

Pith's one-line read The crystalline elastic flow of polygonal curves always converges to a stationary curve.

desk verdict Serious paper with a real gap in the irregular-case convergence proof (Prop. 7.5); the regular case and the classification are solid, but Theorem 7.3 needs repair. read the letter →

arxiv 2506.15869 v1 pith:K3NBWZEX submitted 2025-06-18 math.AP math.CA

Giovanni Bellettini , Shokhrukh Yu. Kholmatov , Matteo Novaga This is my paper · ORCID
classification math.APmath.CA MSC 53E1053E4046N20
keywords crystallinecurvatureelasticflowpolygonalcurvesLojasiewicz-Simoninequalitystationarysolutionstranslatinganisotropicenergysquareanisotropy
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

The paper studies the gradient flow of a crystalline elastic energy on polygonal curves, an energy that adds a squared crystalline-curvature term to the anisotropic length. It establishes that, under the standard assumption that segments translate in parallel, a unique flow exists, that it can be restarted when zero-curvature segments vanish, and that for closed curves it converges as time goes to infinity to a stationary, or generalized stationary, polygonal curve. This matters because it extends long-time convergence results known for smooth elastic flows to the crystalline, nonsmooth setting, and because it gives a complete classification of stationary solutions for the square anisotropy.

What carries the argument

The carrier object is the signed-height ODE system (3.2): each segment $S_i$ of an admissible polygon translates along its normal with speed set by the first variation of $F_\alpha$, so the fourth-order geometric evolution reduces to $n$ ordinary differential equations in the signed heights $h_i$. The load-bearing identities are (3.1), which turns the length and curvature contributions into the transition number $c_i$ times the facet length of the Wulff shape, and the crystalline Lojasiewicz-Simon inequalities (Propositions 6.6 and 7.5), which give quantitative energy decay near stationary curves and force full convergence rather than mere subconvergence. The restart mechanism at the maximal time removes only zero-curvature segments, preserves admissibility, and keeps the index of the curve unchanged.

What would settle it

Numerically integrate (3.2) for a closed admissible polygon with square Wulff shape, choosing initial data for which a zero-curvature segment vanishes at a finite restart time; if the restarted heights do not converge to values satisfying (7.11) for all nondegenerate segments, for instance if they oscillate periodically instead of settling, the long-time convergence claim fails. Alternatively, an explicit initial polygonal curve that develops a new segment direction under the crystalline elastic evolution would falsify the admissibility restriction itself.

Watch

Extended reading notes

Core claim

For a closed admissible polygonal curve, the crystalline elastic flow is globally defined through finitely many restarts and converges in the Kuratowski sense to a closed admissible polygonal curve that is a generalized stationary solution: its nondegenerate segments satisfy the stationarity system (7.11), while any degenerate segment has zero crystalline curvature. If all segment lengths stay bounded away from zero, the limit is a true stationary curve parallel to the initial one, obtained through a crystalline Lojasiewicz-Simon inequality. For the square Wulff shape the paper classifies all stationary solutions as staircases, right-angle chains, double-right-angle chains, and the square Wulff shape of sidelength $\sqrt{4\alpha}$, and it partially classifies translating solutions.

Load-bearing premise

The entire development assumes segments only translate in parallel and that no new edges appear, no segment bends or breaks, and no facet-breaking occurs, so the fourth-order flow is exactly a finite ODE system; if that geometric restriction fails, the flow as defined does not apply.

Editorial extensions

If this is right

  • Every closed admissible polygonal curve has a unique global crystalline elastic evolution, and the index of the curve is preserved through all restarts.
  • If the initial curve is convex, the flow exists globally without restarts and converges to a Wulff shape of radius $\sqrt{\alpha}$, possibly covered multiple times according to the index.
  • In the square-anisotropy case, the only stationary curve with nonzero index is the square Wulff shape of sidelength $\sqrt{4\alpha}$; all other stationary curves have index zero.
  • Bounded polygonal curves cannot be translating solutions, while the one-segment and one-rectangle unbounded examples give explicit grim-reaper-type translating curves.
  • The flow preserves convexity, in contrast with the Euclidean elastic flow, under the admissibility and no-facet-breaking assumptions.

Reading between the lines

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

  • If facet-breaking or spontaneous creation of new segments occurred, the ODE formulation would miss it; the convergence statements therefore indicate stability inside the admissible class, not for arbitrary polygonal evolutions.
  • The same Lojasiewicz-Simon machinery may extend to other crystalline fourth-order flows, such as crystalline surface diffusion, but only under a comparable admissible-class restriction, since facet-breaking is known to occur there.
  • A testable extension suggested by the paper's own conjectures is that unbounded curves with nonparallel half-lines diverge to infinity, while parallel co-directed half-lines should converge to a translating solution after rescaling.
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 / 5 minor

Summary. The paper studies the gradient flow of the crystalline elastic energy F_alpha(Gamma)=∫_Gamma(1+alpha(kappa^phi)^2)phi^o(nu)dH^1, restricted to phi-admissible polygonal curves whose segments translate in the normal direction and for which no facet-breaking, edge creation, or segment bending occurs. The flow is reduced to the finite ODE system (3.2) for signed heights. The authors prove short-time existence and uniqueness (Theorem 4.1), show that at a finite maximal time only zero-curvature segments can vanish, and construct a unique globally defined flow with finitely many restarts (Theorem 4.2). Under a uniform lower bound on segment lengths they prove Kuratowski convergence to a stationary curve via a crystalline Lojasiewicz-Simon inequality (Theorem 7.1); dropping that bound, they prove convergence to a possibly degenerate generalized stationary curve (Theorem 7.3). The final section gives a complete classification of stationary curves and a partial classification of translating curves for the square anisotropy.

Significance. If the proofs are completed, this is a solid contribution to the crystalline analogue of the elastic flow. The paper is largely self-contained: it proves a crystalline Lojasiewicz-Simon-type inequality from scratch, carefully derives the first-variation formula and the energy-dissipation identity, and gives explicit examples including self-similar Wulff-shape evolution and translating solutions. The classification results for the square anisotropy are concrete and checkable. The main theorems are honestly stated under the explicit parallel-translation and no-facet-breaking assumptions, so the restriction to polygonal curves with fixed combinatorial type is a limitation of the model rather than a hidden circularity. The central reservation is the proof of the second Lojasiewicz-Simon inequality, Proposition 7.5, which is the key ingredient for the irregular-case convergence in Theorem 7.3.

major comments (3)
  1. [§7.2, Proposition 7.5, proof of openness of U] The proof asserts that for any ĥ in U, an open neighbourhood of ĥ is obtained by applying Lemma 2.5 to the associated curve Γ̂. Lemma 2.5 requires every segment of Γ̂ to have positive length, because the bound in (2.4) uses min_i H1(S_i). In the setting of Proposition 7.5 the reference curve Γ is a generalized stationary curve and may contain degenerate segments with H1(S_i)=0; Theorem 7.3 explicitly allows exactly this situation. Moreover, for an index i with c_i=0 and |c_{i-1}|=|c_{i+1}|=1, the length formula (2.3) shows that H1(S_i) is independent of h_i and is determined by h_{i-1} and h_{i+1}; the requirement H1(S_i)≥0 is then a one-sided constraint, so the admissible heights form a set with boundary rather than an open subset of R^n. Consequently U is not proved to be open, and the subsequent application of the Lojasiewicz inequality to O=g(U) is not justified as written. The proof needs either a degenerate-segment version of Lemma 2.5 or a direct subanalytic Lojasiewicz argument on a possibly lower-dimensional or boundary-contained set U. This is load-bearing: Proposition 7.5 is the key input in Step 4 of the proof of Theorem 7.3.
  2. [§7.2, Theorem 7.3 and Definitions 2.3–2.5] The paper uses the notions of 'parallel', 'admissible', and 'signed height' for polygonal curves with zero-length segments, but Definition 2.3 explicitly requires all segments to be nondegenerate, and the admissibility and index definitions in Section 2 are stated for ordinary polygonal curves. Theorem 7.3 and Proposition 7.5 need a rigorous extension of these notions to generalized polygonal curves with H1(S_i)=0, including a definition of the associated straight line for a degenerate segment and a precise statement of what it means for a nondegenerate curve to be parallel to such a generalized curve. The sentence in Step 1 that Γ∞ is 'not necessarily parallel to Γ∞' also appears to be a typo, presumably for 'parallel to Γ0'. Without these definitions, the transitivity of parallelness used in Step 4 and the identities (7.14) and (7.17) for degenerate configurations are not fully verifiable.
  3. [§7.2, Step 4 of Theorem 7.3] In the proof of Theorem 7.3, the sets I_k are defined using only the indices i in J, while the indices outside J are controlled later through (7.25). The argument for i outside J requires that the neighbouring segments with nonzero transition numbers have lengths bounded below uniformly in k and in the interval [t_k,r]; this follows from (3.10), but the constants in (7.25) and in the subsequent estimate (7.26) should be shown to be independent of k and of the exit time r. As written the reader must assemble these uniformity claims from the preceding estimates; a short explicit statement would make the proof complete.
minor comments (5)
  1. [§2.5] There is a typo: 'vector vields' should be 'vector fields'.
  2. [§4.1, proof of Theorem 4.1(c), Case 2.1] The displayed condition 'c_{i-1}=c_{i+1}+0' should read 'c_{i-1}=c_{i+1}=0'.
  3. [§7.2, Step 2] The sentence 'The proof runs along the same lines of Proposition 7.4' should refer to Proposition 6.3, not to Proposition 7.4, which is the proposition being proved.
  4. [§4.1, proof of Theorem 4.1] The phrase 'Theferfore it has a unique fixed point' contains a typo; it should be 'Therefore'.
  5. [§2.8, Lemma 2.5] It would help to state explicitly that Lemma 2.5 is applied in later sections only to curves all of whose segments have positive length, since the bound (2.4) is vacuous when min_i H1(S_i)=0. This would clarify the gap identified in the first major comment.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the derivation is self-contained and the long-time convergence conclusions are not forced by fitted parameters, definitional identities, or load-bearing self-citations.

full rationale

Walking the derivation chain shows no circular step of the kind that would make a 'prediction' equal to an input by construction. The crystalline elastic flow is defined by the ODE system (3.2), which is obtained from the first variation of F_alpha; there is no parameter fitted to a subset of data and then relabelled as a prediction. Short-time existence, uniqueness, and the restart mechanism are proved in Theorems 4.1 and 4.2 by contraction arguments inside the paper, not imported from the authors' prior work. The Lojasiewicz-Simon inequalities (Propositions 6.6 and 7.5) are proved internally by applying the standard external Lojasiewicz inequality from [26] to subanalytic sets derived from the length map g; the stationary limit in Theorem 7.3 is obtained by passing to the limit in the evolution equation at a sequence where h'_i(t_k)->0, so stationarity of the limit is a consequence of the dynamics rather than an assumed conclusion. Self-citations such as [4], [5], and [19] appear as background or as pointers to elementary facts that the paper re-derives, e.g. identity (3.1) is verified in the text after mentioning [32]. The only genuine weakness found is a correctness gap rather than circularity: in Proposition 7.5 the proof claims U is open using Lemma 2.5, whose small-height bound (2.4) requires all segment lengths to be positive, while generalized stationary curves may contain degenerate segments and, for i not in J, the heights are determined by admissibility and closure, so U may not be open in R^n. This affects the support of the Lojasiewicz argument in the irregular case of Theorem 7.3, but it does not make the theorem's conclusion equivalent to an input by definition or by self-citation.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new physical entities and fits no free parameters. Its central results rest on standard tools (subanalytic geometry, ODE theory) plus the explicit modeling assumption of non-breaking facets.

assumptions (4)
  • standard math Lojasiewicz inequality for subanalytic functions (see [26])
    Used in Propositions 6.6 and 7.5 to obtain a growth estimate relating the energy gap to the stationarity defects.
  • standard math Existence and uniqueness theory for systems of ODEs with smooth coefficients
    Used in Theorem 4.1 to construct the flow via a contraction mapping on signed heights.
  • standard math Kuratowski convergence and compactness of connected compact sets (see [16])
    Used in Theorems 4.1, 7.1, and 7.3 to pass to subsequential limits of the evolving curves.
  • ad hoc to paper No facet-breaking, no new edge creation, and no segment bending during the evolution
    This modeling restriction, stated in the introduction and Section 3.1, reduces the flow to a finite-dimensional ODE system but is not derived from the energy or from more general principles.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Crystalline elastic flow of polygonal curves: long time behaviour and convergence to stationary solutions." pith.science (2026). https://pith.science/paper/K3NBWZEX

@misc{pith2026250615869,
  author       = {Pith},
  title        = {Pith review of: Crystalline elastic flow of polygonal curves: long time behaviour and convergence to stationary solutions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/K3NBWZEX}},
  note         = {Machine review of arXiv:2506.15869}
}
read the original abstract

Given a planar crystalline anisotropy, we study the crystalline elastic flow of immersed polygonal curves, possibly also unbounded. Assuming that the segments evolve by parallel translation (as it happens in the standard crystalline curvature flow), we prove that a unique regular flow exists until a maximal time when some segments having zero crystalline curvature disappear. Furthermore, for closed polygonal curves, we analyze the behaviour at the maximal time, and show that it is possible to restart the flow finitely many times, yielding a globally in time evolution, that preserves the index of the curve. Next, we investigate the long-time properties of the flow using a Lojasiewicz-Simon-type inequality, and show that, as time tends to infinity, the flow fully converges to a stationary curve. We also provide a complete classification of the stationary solutions and a partial classification of the translating solutions in the case of the square anisotropy.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. One thousand and one higher-order geometric flows of networks

    math.AP 2026-07 conditional novelty 7.0 of 10

    For every order m≥2 and every admissible set of junction/endpoint conditions, the L²-gradient flow of ∫|∂s^m γ|²/2 + length/2 exists uniquely up to reparametrization for a short time, and either runs forever or degene...

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages · cited by 1 Pith paper

  1. [1]

    Andrews: Singularities in crystalline curvature flows

    B. Andrews: Singularities in crystalline curvature flows. Asian J. Math. 6 (2002), 101–121

    work page 2002

  2. [2]

    Bellettini: On facet breaking for crystalline mean curvature flow in 3D

    G. Bellettini: On facet breaking for crystalline mean curvature flow in 3D. Period. Math. Hungar. 48 (2004), 185–206

    work page 2004

  3. [3]

    Bellettini, L

    G. Bellettini, L. Mugnai: Anisotropic geometric functionals and gradient flows. In: Nonlocal and abstract parabolic equations and their applications. Banach Center Publ. 86, 21–43, Polish Acad. Sci. Inst. Math., Warsaw, 2009

    work page 2009

  4. [4]

    Bellettini, M

    G. Bellettini, M. Novaga, M. Paolini: Facet-breaking for three-dimensional crystals evolving by mean curvature. Interfaces Free Bound. 1 (1999), 39–55

    work page 1999

  5. [5]

    Bellettini, M

    G. Bellettini, M. Paolini: Anisotropic motion by mean curvature in the context of Finsler geometry. Hokkaido Math. J. 25 (1996), 537–566

    work page 1996

  6. [6]

    Benois, T

    O. Benois, T. Bodineau, P. Butt´a, E. Presutti: On the validity of the van der Waals theory of surface tension. Markov Process. Related Fields 3 (1997), 175–198

    work page 1997

  7. [7]

    Bierstone, P.D

    E. Bierstone, P.D. Milman: Semianalytic and subanalytic sets. Inst. Hautes ´Etudes Sci. Publ. Math. 67 (1988), 5–42

    work page 1988

  8. [8]

    Cahn: Stability, microstructural evolution, grain growth, and coarsening in a two-dimensional two-phase microstructure

    J. Cahn: Stability, microstructural evolution, grain growth, and coarsening in a two-dimensional two-phase microstructure. Acta Metall. Mater. 39 (1991), 2189–2199

    work page 1991

Show all 32 references
  1. [9]

    Cahn, C.A

    J.W. Cahn, C.A. Handwerker, J.E. Taylor: Overiew no. 98. I – Geometric models of crystal growth. Acta Metall. Mater. 40 (1992), 1443–1474

    1992

  2. [10]

    Cahn, J.E

    J.W. Cahn, J.E. Taylor: Overview no. 113 – Surface motion by surface diffusion. Acta Metall. Mater. 42 (1994), 1045–1063

    1994

  3. [11]

    Carter, A.R

    W.C. Carter, A.R. Roosen, J.W. Cahn, J.E. Taylor: Shape evolution by surface diffusion and surface attachment limited kinetics on completely faceted surfaces. Acta Metall. Mater. 43 (1995), 4309–4323

    1995

  4. [12]

    Dziuk, E

    G. Dziuk, E. Kuwert, R. Sch ¨atzle: Evolution of elastic curves in Rn: existence and computation. SIAM J. Math. Anal. 33 (2002), 1228–1245

    2002

  5. [13]

    Epshteyn, Ch

    Y . Epshteyn, Ch. Liu, M. Mizuno: Motion of grain boundaries with dynamic lattice misorientations and with triple junctions drag. SIAM J. Math. Anal. 53 (2021), 3072–3097

    2021

  6. [14]

    Fonseca, N

    I. Fonseca, N. Fusco, G. Leoni, V . Millot: Material voids in elastic solids with anisotropic surface energies. J. Math. Pures Appl. 96 (2011), 591–639

    2011

  7. [15]

    Garcke, B

    H. Garcke, B. Nestler: A mathematical model for grain growth in thin metallic films. Math. Models Methods Appl. Sci. 10 (2000), 895–921. CRYSTALLINE ELASTIC FLOW 39

    2000

  8. [16]

    Falconer: The Geometry of Fractal Sets

    K. Falconer: The Geometry of Fractal Sets. Cambridge University Press, Cambridge, 1985

    1985

  9. [17]

    M.-H. Giga, Y . Giga: Crystalline surface diffusion flow for graph-like curves. Discrete Contin. Dyn. Syst.43 (2023), 1436–1468

    2023

  10. [18]

    Y . Giga, M. Gurtin, J. Matias: On the dynamics of crystalline motion. Japan J. Indust. Appl. Math. 15 (1998), 7–50

    1998

  11. [19]

    Y . Giga, N. Poˇz´ar: Motion by crystalline-like mean curvature: a survey. Bull. Math. Sc. 12 (2022), 2230004

    2022

  12. [20]

    Y . Giga, P. Rybka:Facet bending in the driven crystalline curvature flow in the plane. J. Geom. Anal.18 (2008), 109–147

    2008

  13. [21]

    Herring: Some theorems on the free energies of crystal surfaces

    C. Herring: Some theorems on the free energies of crystal surfaces. Phys. Rev. 82 (1951), 87–93

    1951

  14. [22]

    Herring: Surface tension as a motivation for sintering

    C. Herring: Surface tension as a motivation for sintering. Fundamental Contributions to the Continuum Theory of Evolving Phase Interfaces in Solids. Springer, Berlin, Heidelberg, 1999

    1999

  15. [23]

    Hironaka: Subanalytic sets

    H. Hironaka: Subanalytic sets. In: Number theory, algebraic geometry and commutative algebra, in honor of Yasuo Akizuki,453–493 (1973), Kinokuniya Book Store, Tokyo

    1973

  16. [24]

    Ishii, H.M

    K. Ishii, H.M. Soner: Regularity and convergence of crystalline motion. SIAM J. Math. Anal. 30 (1998), 19–37

    1998

  17. [25]

    Kinderlehrer, C

    D. Kinderlehrer, C. Liu: Evolution of grain boundaries. Math. Models Methods Appl. Sci. 11 (2001), 713–729

    2001

  18. [26]

    Lojasiewicz: On semi-analytic and subanalytic geometry

    S. Lojasiewicz: On semi-analytic and subanalytic geometry. In: Panoramas of mathematics (Warsaw, 1992/1994). Banach Center Publ. 34 (1995), 89–104, Polish Acad. Sci. Inst. Math., Warsaw

    1995

  19. [27]

    Mantegazza: Smooth geometric evolutions of hypersurfaces

    C. Mantegazza: Smooth geometric evolutions of hypersurfaces. Geom. Funct. Anal. 12 (2002), 138–182

    2002

  20. [28]

    Mantegazza, A

    C. Mantegazza, A. Pluda, M. Pozzetta: A survey of the elastic flow of curves and networks. Milan J. Math. 89 (2021), 59–121

    2021

  21. [29]

    Mantegazza, M

    C. Mantegazza, M. Pozzetta: The Lojasiewicz-Simon inequality for the elastic flow. Calc. Var. Partial Differential Equations 60 (2021)

    2021

  22. [30]

    Mullins: Two-dimensional motion of idealized grain boundaries

    W. Mullins: Two-dimensional motion of idealized grain boundaries. J. Appl. Phys. 27 (1956), 900–904

    1956

  23. [31]

    Taylor: Crystalline variational problems

    J.E. Taylor: Crystalline variational problems. Bull. Amer. Math. Soc. 84 (1978), 568–588

    1978

  24. [32]

    Taylor: Overview no

    J. Taylor: Overview no. 98. II – mean curvature and weighted mean curvature. Acta Metall. Mater. 40 (1992), 1475–1485. (G.Bellettini) U NIVERSITY OF SIENA , V IA ROMA 56, 53100 S IENA , I TALY & I NTERNATIONAL CENTRE FOR THEORETICAL PHYSICS (ICTP), S TRADA COSTIERA 11, 34151 T...

    1992

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.