Pith. sign in

REVIEW 4 major objections 4 minor 36 references

The energy of fractional Allen--Cahn layers in dimension one

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

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

Pith's one-line read The energy of the one-dimensional fractional Allen-Cahn layer is continuous and strictly decreasing in the fractional exponent s, with an explicit pole at s=1/2 and an explicit linear expansion at s=1.

desk verdict A genuinely new quantitative description of E(s), with careful endpoint analysis and a conditional computer-assisted core; referee it, but make sure the code gets independently run. read the letter →

arxiv 2608.06226 v1 pith:VW3YZFDL submitted 2026-08-06 math.AP

Alvaro Carballeira , Damien Galant , Javier Gómez-Serrano This is my paper · ORCID
classification math.AP
keywords energymathcalexplicitallen--cahnendpointfracfractionalkappa
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

Imagine a smooth transition from -1 to +1 along a line, like a boundary between two materials. The fractional Allen-Cahn equation is a version where the boundary feels interactions over a distance, and the strength of that long-range effect is measured by a number s between 0 and 1. For s larger than 1/2, there is exactly one natural transition profile, called the layer, and this paper studies how much energy that layer costs as a function of s.

The answer is clean. The cost is infinite at the critical value s=1/2, and it blows up like 1 divided by (s-1/2). At the other end, s=1, the layer is the familiar tanh curve and the energy is exactly 2 times the square root of 2, divided by 3. As s approaches 1 from below, the energy approaches this value with a linear correction whose slope is given by an explicit integral. In between, the energy is continuous and always decreasing: making the interaction more local makes the interface cheaper.

Proving the always decreasing part is hard because the layer itself is not known in closed form between the endpoints. The authors prove the two endpoint behaviors by hand, with explicit error estimates. On the middle interval they use a computer-assisted argument: they build an accurate approximate layer, prove with interval arithmetic that the true layer is close to it, and then verify finitely many inequalities that force the derivative of the energy to stay negative. The code is attached and uses 200-bit interval arithmetic.

Extended reading notes

Core claim

The energy E(s) of the fractional Allen-Cahn layer is continuous and strictly decreasing on (1/2,1], with E(s)=2*sqrt(2)/3+kappa_1(1-s)+o(1-s) as s approaches 1 from below and E(s)=1/(pi(s-1/2))+O(1) as s approaches 1/2 from above. If correct, E extends to a strictly decreasing homeomorphism from [1/2,1] onto [2*sqrt(2)/3, +infinity), with E(1/2)=+infinity.

Load-bearing premise

The interior proof of strict decrease rests on the interval-arithmetic certificate that conditions (C1)-(C4) of Proposition 5.5 hold at the 56-point grid. In particular, the Schur inequality (85) must hold with theta<1 for the constructed weights, and the fixed-point radius tau(s_i) must be at most 1 with margin eta_i>0. These are finite computations carried out only in the attached code, not reproduced as tables in the text. If any enclosure in the certificate is wrong or the code has a bug, the strict decrease on [0.530,0.980] is not established, and Theorem 1.2 fails.

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

Summary. The paper studies the energy E(s)=E_s[Phi_s] of the one-dimensional fractional Allen-Cahn layer as a function of s in (1/2,1]. It claims that E is continuous and strictly decreasing, with an explicit first-order expansion at s=1 (coefficient kappa_1 given as a convergent integral), and a simple pole 1/(pi(s-1/2)) with bounded remainder at s=1/2. The proof combines analytic endpoint expansions (a corrector argument near s=1 and a Fourier-side pole computation near s=1/2), a semiconcavity argument on the interior interval, and a computer-assisted certificate that reduces strict decrease on [0.530,0.980] to four finite interval-arithmetic conditions (C1)-(C4) at a 56-point grid. The analytic parts are presented in detail with explicit constants; the interior part is the main delicate point.

Significance. If the result is correct, it provides the first sharp quantitative description of the energy of one-dimensional fractional Allen-Cahn layers across the full range s in (1/2,1], including the exact pole residue and the linear correction at s=1. The endpoint arguments are valuable and are carried out with explicit, checkable constants; the pole residue is derived from a Pohozaev identity and the identity that the integral of the layer derivative equals 2, so no fitting parameter enters. The computer-assisted part is also a worthwhile methodological contribution, since it combines rigorous interval arithmetic with a fixed-point validation for a nonlocal semilinear equation on the whole line. The main caveat is that the central theorem depends on a finite interval-arithmetic certificate whose numerical output is not reproduced in the text; until that certificate data and code are independently auditable, the interior monotonicity claim remains conditional on the correctness of the attached computation.

major comments (4)
  1. [Section 5.3, Proposition 5.5; Appendix A] The strict decrease on I0=[0.530,0.980] is exactly the content of conditions (C1)-(C4), but the paper does not present the numerical certificate: no grid values s_i, no values of theta_i from (85), no tau(s_i), eta_i, alpha(s_i), or the cell widths are listed. A reader cannot verify or even audit these load-bearing inequalities without executing the attached code and reproducing the 2h8m run. Please include a table (or a machine-readable data file) of all certified values and margins for every grid point, together with the exact code version or checksum, and state explicitly that Theorem 1.2 rests on those computations. A single failed enclosure at one grid point would invalidate the monotonicity chain.
  2. [Appendix A.2, compact-part estimate] The displayed cell-mass bound 'integral over I_k of K_s(x,y) dy <= integral over I_k of R_s(x-y) dy - inf_{x in I_j} integral over I_k of R_s(x+y) dy' is not, as written, an upper bound uniform in x in I_j: the first term on the right must be the supremum over x in I_j of the integral of R_s(x-y), not an evaluation at an arbitrary or fixed x. If the implementation computes this term at a fixed representative (for example a cell midpoint) and not a rigorous enclosure of the supremum, the Schur certificate could fail. Please clarify the quantifier and confirm that the code encloses the maximum over x in I_j.
  3. [Proposition 5.4 and Appendix A.1] The fixed-point radius tau(s) in Proposition 5.4 depends on omega(s)=6 C*(2s)(||Phi_ap||_Linf + C*(2s)), so a rigorous upper bound on ||Phi_ap||_Linf is needed. The appendix states that the approximate layer is bounded and that the profile is used in the contraction certificate, but it does not record the certified L-infinity enclosure used for ||Phi_ap||_Linf. Since condition (C1) only implies the weak bound |Phi_ap| <= sqrt(5/3), the tighter value used in omega(s) is part of the missing certificate data and should be reported.
  4. [Appendix A.3] The evaluation of J(s,Phi_ap) and the weighted seminorm is a modal exact sum and is a strength of the paper. However, the final sign margin in (C3), eta_i = |J(s_i,Phi_ap)| - L_J(s_i) tau(s_i), is only declared positive; no numerical values of eta_i or of the transfer constant L_J(s_i) are printed. Because C4 uses eta_i in the denominator, the size of the margin controls the admissible cell widths. Please include these values in the certificate table required by the first major comment.
minor comments (4)
  1. [Title page and figures] The title page and figure captions contain obvious OCR artifacts ('LA YERS', 'JA VIER', '©ap', '∙1'); these should be cleaned before publication.
  2. [Equation (95) and grid] The 56-point grid is only specified in the accompanying code. It would improve readability to print the grid values and the chosen cells [a_i,b_i] in an appendix table.
  3. [Theorem 1.1(i)] The approximate value kappa_1 approx 1.073 is stated, but the proof only establishes the lower bound kappa_1 > 119/900. A rigorous enclosure of kappa_1 would make the expansion fully quantitative and is easy to compute from the explicit integral.
  4. [Section 3.3] In Lemma 3.9, the constants A1, A2, A3 are estimated by splitting at |xi|=1; the presentation would be clearer if the displayed values 67, 230, 125 were labeled as upper bounds rather than exact evaluations, since they are rounded. No mathematical issue is involved.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: endpoint asymptotics are derived from explicit integrals and identities, and the computer-assisted interior proof certifies the sign of E' without fitting it as an input.

full rationale

The derivation chain is self-contained. The identity E'(s)=J(s,Phi_s) (Lemma 2.6) follows from an external minimality theorem (Palatucci--Savin--Valdinoci) plus compactness, and the formula contains no fitted parameter. At s=1, the constant kappa_1 is defined as the explicit logarithmic integral in (44) and bounded from below by hand in Lemma 3.5; the expansion is obtained from a fixed-point argument around an explicit corrector with certified error bounds. At s=1/2, the residue 1/pi comes from the identity Integral Phi_s' = 2, the Pohozaev relation P_s=(2s-1)K_s of Lemma 4.2, and explicit Fourier remainders; no asymptotic coefficient is tuned to the target. On the interior interval, the approximate layer is obtained by minimizing the PDE residual, not the sign of J, and the certified inequalities (C1)--(C4) of Proposition 5.5 are used only as sufficient conditions for strict decrease. Thus no prediction is equivalent by construction to an input. The only caveat is an external verification gap: the 56-point interval-arithmetic certificate is implemented in attached code rather than tabulated, so the interior proof is conditional on correct execution of that code. That is a reproducibility concern, not circularity. The cited Cabre--Sire and Palatucci--Savin--Valdinoci results are external, and the self-citations in the survey and related-work discussion are not load-bearing.

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

The theorems carry no fitted constants; the only numbers chosen by hand are computational auxiliaries, approximate layer coefficients and Schur weight scales, that are certified by interval arithmetic before use. The prior theorems from [4,31,30] are standard inputs, and the main unproved computational assumption is the correct execution of the attached certificate code.

free parameters (2)
  • Approximate layer coefficients a_alpha(s_i) = not listed; produced by Gauss-Newton on the PDE residual
    For each of 56 grid points, 36 coefficients are fitted to minimize the L2 residual of the fractional Allen-Cahn equation. They are certified by interval arithmetic and are proof artifacts, not parameters of the theorem.
  • Schur weight scale a(s_i) = not listed; selected by floating-point optimization of the Schur ratio
    The weight w_a in Appendix A.2 is chosen to make the Schur bound theta<1. After selection, the inequality is verified rigorously; this is a computational tuning step, not a free parameter of the conclusion.
assumptions (5)
  • domain assumption Cabre-Sire theorem: existence, uniqueness, odd monotone layer Phi_s with tail decay |1-Phi_s^2| bounded by C|x|^{-2s}.
    Invoked as Theorem 1.5 to define E(s), to justify decay estimates, and to identify fixed-point solutions with the true layer. It is taken from [4] and not reproved.
  • domain assumption Palatucci-Savin-Valdinoci minimality of the layer among admissible profiles.
    Minimality (12) is used for Lipschitz continuity, the Pohozaev scaling identity, and competitor upper bounds. It is taken from [31].
  • domain assumption Spectral gap of the linearized Allen-Cahn operator J on odd functions: sigma(J|_odd)={3/2} union [2,+infinity).
    Used in Lemma 3.7 to exclude spectrum of L_s and in the corrector norm bound. Taken from [30, pp. 73-75].
  • domain assumption Correctness of interval arithmetic and special-function enclosures in Arb and python-flint at 200-bit precision.
    The computer-assisted certificates in Appendix A use rigorous integration, gamma, digamma, and interval evaluations; the proof assumes these are correctly implemented.
  • ad hoc to paper The finite Schur certificates (C1)-(C4) are exactly reproduced by the attached code and certified enclosures.
    Proposition 5.5 states the certificates hold, but the paper does not print the 56-row output table. Acceptance of the interior monotonicity depends on this computational claim.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The energy of fractional Allen--Cahn layers in dimension one." pith.science (2026). https://pith.science/paper/VW3YZFDL

@misc{pith2026260806226,
  author       = {Pith},
  title        = {Pith review of: The energy of fractional Allen--Cahn layers in dimension one},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VW3YZFDL}},
  note         = {Machine review of arXiv:2608.06226}
}
abstract

We study the energy $\mathcal{E}(s) := E_{s}[\Phi_s]$ of the one-dimensional fractional Allen--Cahn layer solution $\Phi_s$, defined as the unique odd, increasing solution of $(-\Delta)^{s} \Phi_s = \Phi_s - \Phi_s^{3}$ with $\Phi_s(\pm\infty)=\pm 1$ and $\Phi_s(0)=0$, for $s\in(1/2,1]$. Our main results are a sharp qualitative and quantitative description of the energy $\mathcal{E}$ on this interval. We show that the energy is continuous and strictly decreasing, and we obtain explicit asymptotic expansions at both endpoints. At the upper endpoint we prove $\mathcal{E}(s) = \frac{2\sqrt{2}}{3} + \kappa_1 (1-s) + o(1-s)$ with an explicit formula for $\kappa_1$. At the lower endpoint we prove that the energy goes to infinity as $\mathcal{E}(s)=\frac{1}{\pi(s-1/2)}+O(1)$. The strict decrease is proved with computer assistance. On an interior subinterval it is reduced to finitely many inequalities verified by interval-arithmetic computations. The proofs combine the Cabr\'e--Sire construction of the layer, the minimality theorem of Palatucci--Savin--Valdinoci, an identity for the $s$ derivative of the energy, and a computer-assisted coercivity estimate at explicit approximate layers.

Figures

Figures reproduced from arXiv: 2608.06226 by the authors.

Figure 1
Figure 1. shows the plot of the curve E(s). 1.2. Preliminaries. Throughout the paper we work in dimension n = 1 and with the param￾eter s ∈ (1/2, 1] unless otherwise stated. We write C for a positive constant that may change from line to line but is independent of any quantity declared variable in its context. We use the Fourier transform and inverse transform ub(ξ) = Z R u(x) e −iξx dx, u(x) = 1 2π Z R ub(ξ) e iξx dξ. 0.5 0.… view at source ↗
Figure 2
Figure 2. The layer tail 1 − Φs(x) on a log–log scale. For s < 1 the decay is algebraic with slope −2s, in agreement with the Cabr´e–Sire bounds (11), and it degenerates to the exponential decay of tanh(x/√ 2) only at s = 1. The curves are computed from the approximate layer Φap s of Appendix A.1. In the next lemma we control these remainders. We need explicit constants in the argument because we need to show that our estimat… view at source ↗
Figure 3
Figure 3. The approximate layer Φap s of (101), for several values of s ∈ [0.53, 0.98], against the exact profile Φ1 = tanh(x/√ 2). The dependence on s is carried mostly by the tails. cell. At the midpoint c, we compute the normalized Taylor jets of Φap s and (−∆)sΦ ap s . The required coefficients are explicit: for k ≥ 1, (102) g (k) (c) k! = 2 π (−1)k+1 k (1 + c 2 ) −k/2 sin kπ 2 − k arctan c  , whereas, for k ≥ 0 and any… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: The weight wa on the positive half-line, at the values a = a(s) selected by the minimization of ρ(a), for six values of s. The coefficient κ normalizes the weight to have a maximum of 1 + ε. Indeed, if b(z) := κz/(1 + z 2 ) 3 , then b ′ (z) = κ 1 − 5z 2 (1 + z 2) 4 , b…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 30 canonical work pages

  1. [1]

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

    work page 1979

  2. [2]

    Buckmaster, G

    T. Buckmaster, G. Cao-Labora, and J. G´ omez-Serrano. Smooth imploding solutions for 3D compressible fluids.Forum Math. Pi, 13:Paper No. e6, 2025

    work page 2025

  3. [3]

    Cabr´ e and Y

    X. Cabr´ e and Y. Sire. Nonlinear equations for fractional Laplacians, I: Regularity, maximum principles, and Hamiltonian estimates.Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire, 31(1):23–53, 2014

    work page 2014

  4. [4]

    Cabr´ e and Y

    X. Cabr´ e and Y. Sire. Nonlinear equations for fractional Laplacians, II: Existence, uniqueness, and qualitative properties of solutions.Trans. Amer. Math. Soc., 367(2):911–941, 2015

    work page 2015

  5. [5]

    Cabr´ e and J

    X. Cabr´ e and J. Sol` a-Morales. Layer solutions in a half-space for boundary reactions.Comm. Pure Appl. Math., 58(12):1678–1732, 2005

    work page 2005

  6. [6]

    M. Cadiot. Constructive proofs of existence and stability of solitary waves in the Whitham and capillary- gravity Whitham equations.Nonlinearity, 38(3):035021, 2025

    work page 2025

  7. [7]

    Global Bifurcation and the Constructive Existence of Overhanging Periodic Steady Water Waves

    M. Cadiot and S. V. Haziot. Global bifurcation and the constructive existence of overhanging periodic steady water waves.arXiv preprint arXiv:2607.13567, 2026

    work page Pith review arXiv 2026

  8. [8]

    Cadiot, J.-P

    M. Cadiot, J.-P. Lessard, and J.-C. Nave. Rigorous computation of solutions of semilinear PDEs on un- bounded domains via spectral methods.SIAM J. Appl. Dyn. Syst., 23(3):1966–2017, 2024

    work page 1966

Show all 36 references
  1. [9]

    L. A. Caffarelli and P. E. Souganidis. Convergence of nonlocal threshold dynamics approximations to front propagation.Arch. Ration. Mech. Anal., 195(1):1–23, 2010

    2010

  2. [10]

    Cao-Labora, M

    G. Cao-Labora, M. Colombo, M. Dolce, and P. Ventura. Instability of two-dimensional Taylor–Green vor- tices.arXiv preprint arXiv:2601.23040, 2026

    2026

  3. [11]

    Carballeira, D

    A. Carballeira, D. Galant, and J. G´ omez-Serrano. energy-fractional-allen-cahn.https://github.com/ AlvaroCarballeira/energy-fractional-allen-cahn, 2026

    2026

  4. [12]

    Castro, D

    A. Castro, D. C´ ordoba, and J. G´ omez-Serrano. Global smooth solutions for the inviscid SQG equation.Mem. Amer. Math. Soc., 266(1292):89 pp., 2020

    2020

  5. [13]

    Castro, J

    A. Castro, J. G´ omez-Serrano, and M. M. G. Pascual-Caballo. Linear instability of a Burgers–Hilbert traveling wave.arXiv preprint arXiv:2605.03920, 2026

    2026 arXiv

  6. [14]

    Cayama, C

    J. Cayama, C. M. Cuesta, and F. de la Hoz. Numerical approximation of the fractional Laplacian onRusing orthogonal families.Appl. Numer. Math., 158:164–193, 2020. ENERGY OF FRACTIONAL ALLEN–CAHN LAYERS IN DIMENSION ONE 49

    2020

  7. [15]

    Chen and T

    J. Chen and T. Y. Hou. Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data I: Analysis.arXiv preprint arXiv:2210.07191, 2022

    2022 arXiv

  8. [16]

    Chen and T

    J. Chen and T. Y. Hou. Stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth data II: Rigorous numerics.Multiscale Model. Simul., 23(1):25–130, 2025

    2025

  9. [17]

    J. Dahne. Highest cusped waves for the fractional KdV equations.J. Differential Equations, 401:550–670, 2024

    2024

  10. [18]

    Dahne and J.-L

    J. Dahne and J.-L. Figueras. Self-similar singular solutions to the nonlinear Schr¨ odinger and the complex Ginzburg–Landau equations.arXiv preprint arXiv:2410.05480, 2024

    2024 arXiv

  11. [19]

    Dahne and J

    J. Dahne and J. G´ omez-Serrano. Highest cusped waves for the Burgers–Hilbert equation.Arch. Ration. Mech. Anal., 247(5):74, 2023

    2023

  12. [20]

    Dahne, J

    J. Dahne, J. G´ omez-Serrano, and J. Pech-Alberich. Monotonicity of the first Dirichlet eigenvalue of regular polygons.arXiv preprint arXiv:2601.16285, 2026

    2026

  13. [21]

    Di Nezza, G

    E. Di Nezza, G. Palatucci, and E. Valdinoci. Hitchhiker’s guide to the fractional Sobolev spaces.Bulletin des sciences math´ ematiques, 136(5):521–573, 2012

    2012

  14. [22]

    Donninger and B

    R. Donninger and B. Sch¨ orkhuber. Self-similar blowup for the cubic Schr¨ odinger equation.arXiv preprint arXiv:2406.16597, 2024

    2024

  15. [23]

    G´ omez-Serrano

    J. G´ omez-Serrano. Computer-assisted proofs in PDE: a survey.SeMA J., 76(3):459–484, 2019

    2019

  16. [24]

    Johansson

    F. Johansson. Arb: efficient arbitrary-precision midpoint-radius interval arithmetic.IEEE Trans. Comput., 66(8):1281–1292, 2017

    2017

  17. [25]

    Kato.Perturbation Theory for Linear Operators, volume 132 ofGrundlehren der mathematischen Wis- senschaften

    T. Kato.Perturbation Theory for Linear Operators, volume 132 ofGrundlehren der mathematischen Wis- senschaften. Springer-Verlag, Berlin-Heidelberg-New York, 2 edition, 1976

    1976

  18. [26]

    P. M. Lushnikov, D. A. Silantyev, and M. Siegel. Collapse versus blow up and global existence in the generalized Constantin–Lax–Majda equation.J. Nonlinear Sci., 31(5):Art. 82, 2021

    2021

  19. [27]

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

    1987

  20. [28]

    Modica and S

    L. Modica and S. Mortola. Un esempio di Γ-convergenza.Boll. Un. Mat. Ital. B (5), 14:285–299, 1977

    1977

  21. [29]

    M. T. Nakao, M. Plum, and Y. Watanabe.Numerical Verification Methods and Computer-Assisted Proofs for Partial Differential Equations, volume 53 ofSpringer Series in Computational Mathematics. Springer, Singapore, 2019

    2019

  22. [30]

    A. F. Nikiforov and V. B. Uvarov.Special Functions of Mathematical Physics: A Unified Introduction with Applications. Birkh¨ auser, Basel, 1988

    1988

  23. [31]

    Palatucci, O

    G. Palatucci, O. Savin, and E. Valdinoci. Local and global minimizers for a variational energy involving a fractional norm.Ann. Mat. Pura Appl. (4), 192(4):673–718, 2013

    2013

  24. [32]

    Savin and E

    O. Savin and E. Valdinoci. Γ-convergence for nonlocal phase transitions.Ann. Inst. H. Poincar´ e Anal. Non Lin´ eaire, 29(4):479–500, 2012

    2012

  25. [33]

    Savin and E

    O. Savin and E. Valdinoci. Density estimates for a variational model driven by the Gagliardo norm.J. Math. Pures Appl. (9), 101(1):1–26, 2014

    2014

  26. [34]

    Sire and E

    Y. Sire and E. Valdinoci. Fractional Laplacian phase transitions and boundary reactions: a geometric in- equality and a symmetry result.J. Funct. Anal., 256(6):1842–1864, 2009

    2009

  27. [35]

    E. M. Stein.Singular Integrals and Differentiability Properties of Functions, volume 30 ofPrinceton Mathe- matical Series. Princeton University Press, Princeton, NJ, 1970

    1970

  28. [36]

    Tucker.Validated Numerics: A Short Introduction to Rigorous Computations

    W. Tucker.Validated Numerics: A Short Introduction to Rigorous Computations. Princeton University Press, Princeton, NJ, 2011. Department of Mathematics, Brown University, 151 Thayer Street, Providence, RI 02912, USA Email address:alvaro carballeira@brown.edu Department of Math...

    2011

Pith tools

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