Pith. sign in

REVIEW 4 major objections 5 minor 22 references

Phase-field approximation of an epitaxial growth model with adatoms

T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read A smoothed phase-field energy is proved to converge to the sharp epitaxial-growth model with adatoms as the smoothing vanishes.

desk verdict A plausible Gamma-convergence framework for phase-field epitaxy with adatoms, but the mass-constrained version breaks because the recovery sequence never satisfies the stated constraint on u. read the letter →

arxiv 2411.18491 v3 pith:IK6S3DWS submitted 2024-11-27 math.AP

Gabriele Fissore This is my paper · ORCID
classification math.AP MSC 49J4574G6549Q20
keywords phase-fieldapproximationGamma-convergenceepitaxialgrowthadatomslinearisedelasticityfreediscontinuityproblemsmassconstraintsurfaceenergy
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

This paper proves a rigorous phase-field approximation, in the sense of $\Gamma$-convergence, for the sharp variational model of epitaxial growth with adatoms studied in the author's companion paper [9]. The smooth functionals $F_\varepsilon$ and their mass-constrained versions $F_\varepsilon^{m,M}$ $\Gamma$-converge, as $\varepsilon\to0$, to the sharp relaxed functionals $F$ and $F^{m,M}$: the phase-field interface collapses to the film boundary, the displacement converges in $L^2_{\rm loc}$, and the weighted adatom measure converges weak-* to the limit measure. The payoff is that a numerically friendly model whose variables are all functions retains the exact surface energetics of the sharp model, including vertical fractures, the convex sub-additive envelope of the adatom cost, and the singular part of the adatom measure. As a direct corollary, the minimal energies of the constrained phase-field problems converge to the minimum of the sharp constrained problem, and minimizing sequences pass to sharp minimizers.

What carries the argument

The carrying mechanism is the Modica–Mortola-type interfacial energy $\frac1\sigma\int_{Q_+}(\varepsilon|\nabla w|^2+\frac1\varepsilon P(w))\psi(u)\,dx$, where $P$ is a double-well potential vanishing only at $0$ and $1$ and $\sigma$ normalizes the one-dimensional profile. The double well forces $w$ to spend its energy in a transition layer of width $\varepsilon$ that collapses onto the film boundary, and multiplying by $\psi(u)$ turns that layer into the adatom measure $\mu_\varepsilon$; in the limit, the cost of a density $u$ on the regular interface is $\widetilde\psi(u)$, while vertical cuts are handled by splitting the layer left and right of the cut and combining the two one-sided costs through $\psi^c$. Recovery sequences are built from the almost optimal profile $\gamma_\varepsilon$ solving $\varepsilon^2|\gamma_\varepsilon'|^2=P(\gamma_\varepsilon)+\sqrt\varepsilon$ with $\gamma_\varepsilon(0)=1$, $\gamma_\varepsilon(1)=0$, composed with the signed distance to $\Omega$; the 'wriggling process' from the earlier papers lets one approximate densities above the linear threshold of $\widetilde\psi$ by repeated copies of the threshold value $s_0$.

What would settle it

Take $\psi\equiv1$ and put $u_\varepsilon=m/\varepsilon$ on the $\varepsilon$-transition layer of $w_\varepsilon$; the surface-energy term of $F_\varepsilon$ stays bounded while $\mu_\varepsilon$ has mass of order $1/\varepsilon$, so no subsequence can weak-* converge with support on the interface. This calculation makes concrete the failure flagged in Remark 6: without hypothesis (12), the compactness asserted in the equi-coercivity step behind Theorem 7 cannot hold.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that the phase-field functionals $$F_\varepsilon(w,v,u)=\int_{Q_+}(w+\eta_\varepsilon)W(E(v)-E_0(y))\,dx+\frac1\$\sigma$\int_{Q_+}\left(\varepsilon|\nabla w|^2+\frac1\varepsilon P(w)\right)\psi(u)\,dx,$$ with $\eta_\varepsilon=o(\varepsilon)$ and $\sigma=2\int_0^1\sqrt{P(t)}\,dt$, have $\Gamma$-limit equal to the relaxed sharp energy $F$ from the earlier work, under the convergence of Definition 21: $w_\varepsilon\to\chi_\Omega$ in $L^1_{\rm loc}$, $v_\varepsilon\to v$ in $L^2_{\rm loc}$, and the measures $\mu_\varepsilon=\frac{u_\varepsilon}{\sigma}(\varepsilon|\nabla w_\varepsilon|^2+\frac1\varepsilon P(w_\varepsilon))\,\mathcal L^2\llcorner Q_+$ converge weak-* to $\mu$. The limit separates the surface energy into three contributions: $\widetilde\psi(u)$ on the regular part $\widetilde\Gamma$, $\psi^c(u)$ on the cut, fracture-like part $\Gamma^c$, and $\theta\,\mu_s$ on the singular part, where $\widetilde\psi$ is the convex sub-additive envelope of $\psi$ and $\psi^c(s)=\min\{\widetilde\psi(r)+\widetilde\psi(t): r+t=s\}$. The same statement holds under the mass constraints $\int_{Q_+}w=M$ and $\int_{Q_+}u=m$, and then the minima of the constrained phase-field problems converge to the sharp constrained minimum.

Load-bearing premise

The paper assumes—and its own Remark 6 says the assumption is essential—that the adatom measures from bounded-energy phase-field sequences have weak-* limits supported on the limiting film interface; this support condition is what lets the mass constraint pass to the limit and is not a consequence of the energy bounds.

Editorial extensions

If this is right

  • For every admissible sharp configuration $(\Omega,v,\mu)$ with finite energy and fixed mass $m,M$, there is a recovery sequence of phase-field triples respecting the exact mass constraints whose energy converges to the sharp value, so the constrained phase-field problem is a faithful variational approximation.
  • Any sequence of minimizers of $F_\varepsilon^{m,M}$ with bounded elastic energy and with adatom measures converging to an interface-supported measure clusters at minimizers of $F^{m,M}$, and the minimal energies converge.
  • The phase-field limit sees fractures: the cut part of a BV graph contributes $\psi^c(u)$ per unit length, so vertical cuts inside the film are energetically visible even though the approximating sequence uses smooth functions.
  • The approximation works for any Borel cost $\psi\colon[0,\infty)\to(0,\infty)$ with positive infimum; the convexity and sub-additivity that the sharp limit needs are generated by the $\Gamma$-limit itself, not imposed on the approximating densities.
  • The role of $\eta_\varepsilon=o(\varepsilon)$ is to provide compactness for the displacement and phase-field sequences; removing it would jeopardize the convergence of minimizing sequences.

Reading between the lines

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

  • An unconditional compactness result would need a coupling between $u_\varepsilon$ and the interfacial layer, such as an $L^2$ or $L^1$ penalization of $u_\varepsilon$ away from the transition region; without it, bounded energy allows the adatom mass to concentrate off the interface, as the example $\psi\equiv1$, $u_\varepsilon=m/\varepsilon$ on the layer shows.
  • The same splitting-at-a-cut argument could be transplanted to higher-dimensional epitaxial models once one has a local $\Gamma$-convergence statement for pairs of measures and energy measures, suggesting the result is not tied to the one-dimensional graph setting.
  • Numerically, the theorem justifies using the smooth constrained functional as a computational proxy: choose $\varepsilon$ small compared with island features, enforce the mass constraints on $w$ and $u$, and minimize; the $\Gamma$-limit guarantees that sharp equilibria are recovered in the limit.
  • The infimal convolution $\psi^c$ on cuts could be read as a selection rule for fracture patterns: when adatom densities are high, the energetically cheapest way to open a crack is to split the density between the two sides, a tendency that may appear in phase-field simulations as symmetric crack openings.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The paper proposes a phase-field approximation for a sharp model of epitaxial growth with adatoms previously relaxed in [9]. The phase-field functional couples a double-well Modica-Mortola surface term to an adatom density u through a coefficient ψ(u), and the main claimed results are Γ-convergence of the unconstrained functionals F_ε to F (Theorem 4), Γ-convergence of the constrained functionals F_ε^{m,M} to F^{m,M} (Theorem 6), and convergence of minima (Theorem 7), all with respect to the topology of Definition 21. The proof strategy adapts the methods of [3] and [5] to the relaxed functional of [9].

Significance. If the results were valid, the paper would provide a useful phase-field approximation of a nontrivial free-boundary model with adatoms, including a mass constraint, and would extend the approaches of Bonnetier-Chambolle and Caroccia-Cristoferi to the relaxed functional with cuts. The paper has no free parameters and relies on a previously published relaxation theorem rather than circular reasoning. However, the constrained version of the result is currently not established: the adatom mass constraint in Definition 19 is imposed on the unweighted L1 integral of u, while the object that passes to the limit is the energy-weighted measure μ_ε of Definition 21(iii). This mismatch breaks the compactness argument, the claimed equi-coercivity, and the admissibility of the recovery sequence.

major comments (4)
  1. [§5, Eq. (13)] The equality μ_ε(Q+)=m is not a consequence of the constraint ∫_{Q+}u_ε dx=m. By Definition 19 the constraint is on the unweighted L1 integral of u_ε, while μ_ε is the energy-weighted measure defined in Definition 21(iii). Since ψ is only assumed bounded below, the energy bound (10) does not control μ_ε(Q+). Taking ψ≡1 and u_ε=m/|A_ε| on the transition layer A_ε with |A_ε|≈εH^1(Γ) gives ∫_{Q+}u_ε dx=m and sup F_ε<∞, but μ_ε(Q+) tends to infinity; thus (13) is false even when the weak-* convergence in the hypotheses is assumed, because the limit of μ_ε(Q+) need not equal m. Consequently Theorem 5 does not establish that the limiting measure has mass m.
  2. [Before Theorem 7, Definition 2] The remark that the family (F_ε^{m,M})_ε is equi-coercive is false under the topology of Definition 21. The functional does not control the total mass of μ_ε: the example in the previous comment gives sup F_ε<∞ and ∫_{Q+}u_ε dx=m but no weak-* convergent subsequence of μ_ε. The same example shows that Remark 5 is incorrect, since strict positivity of ψ is not sufficient for coerciveness. A growth condition on ψ or a constraint imposed directly on μ_ε(R^2) would be needed before Theorem 7 can be derived from Theorems 2, 5 and 6.
  3. [§7, Theorem 15, Step 3 (Eqs. (51)–(54))] The recovery sequence is not checked against the admissibility condition in Definition 19. The computation verifies μ_ε(R^2)=m, but A_p(m,M) requires ∫_{Q+}u_ε dx=m. With the definition u_ε|_{R_j}=u_j H^1(Γ∩R_j)/p_j^ε on each grid rectangle R_j, one has ∫_{Q+}u_ε dx≈Σ_j u_j |R_j∩Q_+| (since p_j^ε→H^1(Γ∩R_j)), whereas m=Σ_j u_j H^1(Γ∩R_j). For a flat interface segment contained in a tall rectangle these two quantities differ in general. Hence the constructed (w_ε,v_ε,u_ε) is generally not in A_p(m,M), and the limsup inequality for Theorem 6 is not proven.
  4. [§6, Theorem 10, Eq. (15)] The lower bound for the bulk term is not justified by the stated hypotheses. The sequence is only assumed to satisfy v_ε→v in L2_loc and w_ε→χ_Ω in L1_loc; no bound on E(v_ε) is available. The displayed argument 'by convexity and liminf w_ε=1' is not sufficient for the weak lower semicontinuity of the quadratic bulk energy. If the liminf is finite, one can extract a weighted L2 control on E(v_ε) on compact subsets of Ω and then use a truncation or Egorov argument, but this reasoning is absent from the proof. This gap affects the liminf inequality for Theorem 4 as well.
minor comments (5)
  1. [Abstract/Introduction] There is a typo in the introduction: 'we condiderv∈H1(Ω;R2)' should read 'we consider v∈H1(Ω;R2)'.
  2. [§7, Eq. (33)] The initial conditions γ_ε(0)=1, γ_ε(1)=0 in the almost optimal profile problem are stated without comment, while the earlier equation (9) uses the opposite convention. The sign convention is consistent with the subsequent extension γ_ε=1 for x<0 and γ_ε=0 for x>1, but the inconsistency should be resolved explicitly.
  3. [§6, Step 2.1] In the proof of Theorem 10, the rectangle R_r is defined with the same x-interval as R_l; it should presumably be (x_c, x_c+δ) instead of (x_c-δ, x_c).
  4. [§6, Step 2] The sentence 'C_{ξγ} is monotonically converging to the empty set, as ε→0' should refer to γ→0, not ε→0.
  5. [§7, Theorem 14] The statement 'lim_{ε→∞} H(Ω_ε,v_ε,μ_ε)=...' should be lim_{ε→0}, and the reference 'where H is defined (14)' appears to point to the wrong equation, since (14) is the Γ-limit E from [5].

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the phase-field Gamma-convergence is argued from external published results ([3], [5], [20]) and the author's prior relaxation theorem [9]; no equation reduces to its own input.

full rationale

The central Gamma-convergence claim is not equivalent to its inputs. The sharp functional G is taken from the externally published relaxation result [9]; the phase-field surface-energy Gamma-convergence is quoted from [5, Theorem 3.10] and [20, Proposition 2], and the elasticity coupling follows the strategy of [3]. These are independent, stated-assumption results that do not contain the target theorem. The author's use of [9] is self-citation, but it is load-bearing only as a source of the sharp model and of approximation lemmas whose proofs are given there; it is not used to assert the phase-field Gamma-limit. No parameter is fitted to data and no limiting object is defined in terms of the functional being approximated. The main correctness risk, highlighted in the paper itself by Remark 6 and Figure 1, is that bounded F_eps energy does not control the mass of the weighted measure mu_eps, so the family is not equi-coercive and Theorem 7 is not justified by Theorem 5; additionally, the recovery sequence in Theorem 15 verifies mu_eps(R^2)=m but not the stated admissibility condition int_{Q+} u_eps dx = m (Definition 19 versus Definition 21(iii)). These are mathematical gaps, not circular reductions: the claimed results are not assumed in the hypotheses, and fixing them would require additional arguments rather than exposing a tautology. Hence the circularity score is low.

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

No free parameters or invented entities. The central claim rests on external theorems from [5] and [9], including the author's own prior relaxation result, and on an unstated compactness condition that is not implied by the stated hypotheses.

assumptions (5)
  • domain assumption Relaxation theorem from [9]: the sharp functional F is the relaxation of H under the topology of Definition 17
    The target Gamma-limit G is imported from the author's prior paper; the present paper does not reprove this relaxation result.
  • domain assumption Phase-field Gamma-convergence result from [5] (Theorem 9 here) for surface energies with adatom density
    Used as a black box in the liminf inequality for the regular part of the interface.
  • domain assumption Approximation results from [9] (Propositions 7.6, 7.7, 7.8) reducing recovery sequences to smooth graphs with grid-constant densities
    The limsup proof depends on these prior approximation theorems.
  • ad hoc to paper Unstated growth or compactness control for the adatom measure
    Theorem 7 equi-coercivity requires either a growth condition on psi, such as psi(u) bounded below by c u for large u, or a precompactness hypothesis on mu_eps; with only inf psi > 0 the energy does not bound mu_eps mass.
  • standard math Standard BV, Korn, and Gamma-convergence compactness facts
    Used in Theorem 8 for compactness of w_eps and v_eps; cited to [1], [4], [18], and [20].

how reviews work

0 comments
Cite this review

Pith. "Pith review of Phase-field approximation of an epitaxial growth model with adatoms." pith.science (2026). https://pith.science/paper/IK6S3DWS

@misc{pith2026241118491,
  author       = {Pith},
  title        = {Pith review of: Phase-field approximation of an epitaxial growth model with adatoms},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IK6S3DWS}},
  note         = {Machine review of arXiv:2411.18491}
}
read the original abstract

We perform a mass constrained phase-field approximation for a model that describes the epitaxial growth of a two-dimensional thin film on a substrate in the context of linearised elasticity. The approximated model encodes a variable on the free surface of the film, that physically is interpreted as an adatom density.

Figures

Figures reproduced from arXiv: 2411.18491 by the authors.

Figure 1
Figure 1. The phase-field variable converges to χΩ , whereas the uε , being independent from the other variables, can concentrate in an area away from Γ. Proof. Step 1. The compactness of the phase-field variable wε follows a standard argument, which makes use of the first bound in (10), and can be found in [20]. Therefore, we have the existence of w ∈ BV(Q+; {0, 1}) such that, up to a subsequence wε → χΩ in L 1 (Q+). Moreove… view at source ↗
Figure 2
Figure 2. For each x, the function weε(x, y) = wε(x, y − ℓy) allows us to se the potential term W , where ∇weε = 0. = Z Q∩Kε |∇v| 2 dx ⩽ ||∇v||L2(Q) . We estimate the second one and then conclude for the entire right-hand side of (41). We have that Z Q∩K |vε(x, y) ⊗ ∇weε(x, y)| 2 dx ⩽ Z Q∩K |vε(x, y)| 2 |∇weε(x, y)| 2 dx ⩽ Z Q∩Kε |v(x, y)| 2 [PITH_FULL_IMAGE:figures/full_fig_p021_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

  1. [9]

    Cristoferi and G

    R. Cristoferi and G. Fissore , Two-dimensional graph model for epitaxial crystal growth with adatoms , Adv. Calc. Var., 18 (2025), pp. 511--554

    work page 2025

  2. [3]

    Bonnetier and A

    E. Bonnetier and A. Chambolle , Computing the equilibrium configuration of epitaxially strained crystalline films , SIAM J. Appl. Math., 62 (2002), pp. 1093--1121

    work page 2002

  3. [5]

    Caroccia and R

    M. Caroccia and R. Cristoferi , On the gamma convergence of functionals defined over pairs of measures and energy-measures , J. Nonlinear Sci., 30 (2020), pp. 1723--1769

    work page 2020

  4. [1]

    Ambrosio, N

    L. Ambrosio, N. Fusco, and D. Pallara , Functions of bounded variation and free discontinuity problems , Oxford Mathematical Monographs, The Clarendon Press Oxford University Press, New York, 2000

    work page 2000

  5. [2]

    Ambrosio and V

    L. Ambrosio and V. M. Tortorelli , Approximation of functionals depending on jumps by elliptic functionals via -convergence , Comm. Pure Appl. Math., 43 (1990), pp. 999--1036

    work page 1990

  6. [4]

    Brezis , Functional analysis, S obolev spaces and partial differential equations , Universitext, Springer, New York, 2011

    H. Brezis , Functional analysis, S obolev spaces and partial differential equations , Universitext, Springer, New York, 2011

    work page 2011

  7. [6]

    Caroccia, R

    M. Caroccia, R. Cristoferi, and L. Dietrich , Equilibria configurations for epitaxial crystal growth with adatoms , Arch. Ration. Mech. Anal., 230 (2018), pp. 785--838

    work page 2018

  8. [7]

    Conti, M

    S. Conti, M. Focardi, and F. Iurlano , Phase field approximation of cohesive fracture models , Ann. Inst. H. Poincar\'e C Anal. Non Lin\'eaire, 33 (2016), pp. 1033--1067

    work page 2016

Show all 22 references
  1. [8]

    Crismale, M

    V. Crismale, M. Friedrich, and F. Solombrino , Integral representation for energies in linear elasticity with surface discontinuities , Adv. Calc. Var., 15 (2022), pp. 705--733

    2022

  2. [10]

    Dal Maso , An introduction to -convergence , Progress in Nonlinear Differential Equations and their Applications, 8, Birkh\"auser Boston, Inc., Boston, MA, 1993

    G. Dal Maso , An introduction to -convergence , Progress in Nonlinear Differential Equations and their Applications, 8, Birkh\"auser Boston, Inc., Boston, MA, 1993

    1993

  3. [11]

    height 2pt depth -1.6pt width 23pt, Generalised functions of bounded deformation , J. Eur. Math. Soc. (JEMS), 15 (2013), pp. 1943--1997

    2013

  4. [12]

    De Lellis , Rectifiable sets, densities and tangent measures , Zurich Lectures in Advanced Mathematics, European Mathematical Society (EMS), Z\"urich, 2008

    C. De Lellis , Rectifiable sets, densities and tangent measures , Zurich Lectures in Advanced Mathematics, European Mathematical Society (EMS), Z\"urich, 2008

    2008

  5. [13]

    Focardi , On the variational approximation of free-discontinuity problems in the vectorial case , Math

    M. Focardi , On the variational approximation of free-discontinuity problems in the vectorial case , Math. Models Methods Appl. Sci., 11 (2001), pp. 663--684

    2001

  6. [14]

    Fonseca, N

    I. Fonseca, N. Fusco, G. Leoni, and M. Morini , Equilibrium configurations of epitaxially strained crystalline films: existence and regularity results , Arch. Ration. Mech. Anal., 186 (2007), pp. 477--537

    2007

  7. [15]

    Fonseca and G

    I. Fonseca and G. Leoni , Modern methods in the calculus of variations: L^p spaces , Springer Monographs in Mathematics, Springer, New York, 2007

    2007

  8. [16]

    Fried and M

    E. Fried and M. E. Gurtin , A unified treatment of evolving interfaces accounting for small deformations and atomic transport with emphasis on grain-boundaries and epitaxy , Advances in applied mechanics, 40 (2004), pp. 1--177

    2004

  9. [17]

    M. A. Grinfeld , Stress driven instability in elastic crystals: mathematical models and physical manifestations , J. Nonlinear Sci., 3 (1993), pp. 35--83

    1993

  10. [18]

    Leoni , A first course in S obolev spaces , vol

    G. Leoni , A first course in S obolev spaces , vol. 181 of Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, second ed., 2017

    2017

  11. [19]

    Maggi , Sets of finite perimeter and geometric variational problems , vol

    F. Maggi , Sets of finite perimeter and geometric variational problems , vol. 135 of Cambridge Studies in Advanced Mathematics, Cambridge University Press, Cambridge, 2012. An introduction to geometric measure theory

    2012

  12. [20]

    Modica , The gradient theory of phase transitions and the minimal interface criterion , Arch

    L. Modica , The gradient theory of phase transitions and the minimal interface criterion , Arch. Rational Mech. Anal., 98 (1987), pp. 123--142

    1987

  13. [21]

    Mumford and J

    D. Mumford and J. Shah , Optimal approximations by piecewise smooth functions and associated variational problems , Comm. Pure Appl. Math., 42 (1989), pp. 577--685

    1989

  14. [22]

    Spencer and J

    B. Spencer and J. Tersoff , Equilibrium shapes and properties of epitaxially strained islands , Phys. Rev. Lett., 79 (1997), pp. 4858--4861

    1997

Pith tools

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