Recognition: unknown
Singularities in phase separation models: a spectral element approach for the nonlocal Cahn-Hilliard equation
Pith reviewed 2026-05-10 01:57 UTC · model grok-4.3
The pith
The nonlocal Cahn-Hilliard equation with singular long-range interactions can be simulated at high resolution using a pseudospectral multishape method that requires only limited computational resources.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors introduce an efficient and error-controlled numerical framework for the nonlocal Cahn-Hilliard system with constant mobility, logarithmic potential, Newtonian interaction kernel, and no-flux boundary conditions. Their approach rests on a pseudospectral multishape method that accurately approximates the action of singular convolution operators, and they show that high-resolution numerical solutions can be achieved with limited computational resources.
What carries the argument
The pseudospectral multishape method, which divides the computational domain to approximate the dense singular convolution operator without special post-processing.
If this is right
- High-resolution numerical solutions of the nonlocal system become practical on standard hardware.
- The singular Newtonian kernel can be treated accurately within the chosen discretization without post-hoc adjustments.
- Constant-mobility logarithmic-potential cases with no-flux boundaries admit controlled error behavior at moderate cost.
- The framework supports systematic exploration of long-range interaction effects in phase separation dynamics.
Where Pith is reading between the lines
- The same discretization strategy could be tested on other singular kernels or boundary conditions to check broader applicability.
- Practical high-resolution runs might allow direct comparison against existence and regularity theorems for the continuous model.
- Reducing the computational barrier could enable parameter sweeps that reveal how interaction range affects coarsening rates.
- The multishape idea may transfer to related nonlocal evolution equations that appear in fluid or biological models.
Load-bearing premise
The multishape pseudospectral discretization controls approximation errors near the kernel singularity and stays stable for the Newtonian interaction under no-flux boundary conditions without hidden instabilities or extra fixes.
What would settle it
Running the scheme on a benchmark problem with a known reference solution and observing either growing pointwise errors near the singularity or unphysical oscillations in the phase field would show the accuracy claim does not hold.
Figures
read the original abstract
The nonlocal Cahn-Hilliard equation provides a natural extension of the classical model for phase separation by incorporating long-range interactions through a singular convolution kernel. While this formulation admits a rich existence and regularity theory, its numerical approximation remains challenging: discretisation of the nonlocal term leads to dense operators, and the singularity of the kernel requires special treatment in collocation-based schemes. In this work, we develop an efficient and error-controlled numerical framework for the nonlocal Cahn-Hilliard system with constant mobility, logarithmic potential, Newtonian interaction kernel, and no-flux boundary conditions. Our approach is based on a pseudospectral multishape method that accurately approximates the action of singular convolution operators. We present high-resolution numerical solutions for this nonlocal system of equations that can be achieved with limited computational resources.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a pseudospectral multishape discretization for the nonlocal Cahn-Hilliard equation with Newtonian interaction kernel, logarithmic potential, constant mobility, and no-flux boundary conditions. It claims this yields an efficient, error-controlled framework that produces high-resolution numerical solutions with limited computational resources.
Significance. If the error-control and accuracy claims for the singular convolution operator are substantiated, the work would provide a practical advance in the numerical treatment of nonlocal phase-field models, enabling reliable high-resolution simulations of long-range interactions that are otherwise hindered by dense operators and kernel singularities.
major comments (1)
- [Abstract] Abstract: The assertion of an 'error-controlled numerical framework' and 'accurate' approximation of singular convolution operators lacks any supporting a priori error bounds, convergence-rate analysis, benchmark comparisons, or quantitative tables (e.g., error versus polynomial degree or element count). This is load-bearing for the central claim, as the skeptic correctly identifies that the method's ability to control quadrature errors near the kernel origin and domain boundaries under no-flux conditions remains unverified.
Simulated Author's Rebuttal
We thank the referee for their thorough review and constructive criticism. We address the single major comment below and agree that revisions are needed to align the abstract claims more closely with the manuscript's content.
read point-by-point responses
-
Referee: [Abstract] Abstract: The assertion of an 'error-controlled numerical framework' and 'accurate' approximation of singular convolution operators lacks any supporting a priori error bounds, convergence-rate analysis, benchmark comparisons, or quantitative tables (e.g., error versus polynomial degree or element count). This is load-bearing for the central claim, as the skeptic correctly identifies that the method's ability to control quadrature errors near the kernel origin and domain boundaries under no-flux conditions remains unverified.
Authors: We acknowledge the validity of this observation. The manuscript presents numerical experiments in Section 4 demonstrating the practical performance of the pseudospectral multishape method through high-resolution simulations and comparisons with reference solutions, including observed convergence behavior with respect to polynomial degree. However, no a priori error bounds or formal convergence-rate analysis are derived, and the abstract's phrasing overstates the level of theoretical control. We will revise the abstract by replacing 'error-controlled numerical framework' with 'efficient numerical framework' and 'accurate' with 'effective', while adding explicit references to the numerical validation studies. We will also insert a brief paragraph in the conclusions noting that rigorous a priori analysis of quadrature errors near singularities and boundaries is left for future work. These changes ensure the claims are supported by the presented evidence without misrepresentation. revision: yes
Circularity Check
No significant circularity: numerical discretization method is self-contained
full rationale
The paper proposes and demonstrates a pseudospectral multishape discretization for approximating singular convolution operators in the nonlocal Cahn-Hilliard equation. No derivation chain exists in which a claimed prediction or first-principles result reduces by construction to fitted inputs, self-citations, or renamed ansatzes. The central claims concern the efficiency and error control of the numerical scheme for the given kernel, potential, and boundary conditions, presented via high-resolution solutions rather than any self-referential mathematical reduction. The approach is a standard numerical framework development without load-bearing circular steps of the enumerated kinds.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The nonlocal Cahn-Hilliard equation with Newtonian kernel admits a rich existence and regularity theory
- domain assumption Pseudospectral methods can be extended to multishape domains to approximate singular convolution operators
Reference graph
Works this paper leans on
-
[1]
2011 , PAGES =
Bahouri, Hajer and Chemin, Jean-Yves and Danchin, Rapha\"el , TITLE =. 2011 , PAGES =
2011
-
[3]
Lions and E
J.-L. Lions and E. Magenes , publisher =. Problèmes aux. 1968 , series =
1968
-
[4]
Journal of Differential Equations , VOLUME =
Davoli, Elisa and Scarpa, Luca and Trussardi, Lara , TITLE =. Journal of Differential Equations , VOLUME =. 2021 , PAGES =
2021
-
[5]
Archive for Rational Mechanics and Analysis , VOLUME =
Davoli, Elisa and Scarpa, Luca and Trussardi, Lara , TITLE =. Archive for Rational Mechanics and Analysis , VOLUME =. 2021 , NUMBER =
2021
-
[6]
Annales de l'Institut Henri Poincar\'e
Brezis, Ha\"im and Mironescu, Petru , TITLE =. Annales de l'Institut Henri Poincar\'e. 2018 , NUMBER =
2018
-
[7]
2026 , eprint =
Grasselli, Maurizio and Poiatti, Andrea , title =. 2026 , eprint =
2026
-
[8]
, TITLE =
Gajewski, Herbert and Griepentrog, Jens A. , TITLE =. Discrete and Continuous Dynamical Systems , VOLUME =. 2006 , NUMBER =
2006
-
[9]
Nonlinearity , VOLUME =
Della Porta, Francesco and Giorgini, Andrea and Grasselli, Maurizio , TITLE =. Nonlinearity , VOLUME =. 2018 , NUMBER =
2018
-
[10]
Analysis & PDE , VOLUME =
Poiatti, Andrea , TITLE =. Analysis & PDE , VOLUME =. 2025 , NUMBER =
2025
-
[11]
ESAIM: Control, Optimisation and Calculus of Variations , VOLUME =
Poiatti, Andrea and Signori, Andrea , TITLE =. ESAIM: Control, Optimisation and Calculus of Variations , VOLUME =. 2024 , PAGES =
2024
-
[12]
Gal, C. G. and Giorgini, A. and Grasselli, M. and Poiatti, A. , title =. Journal de Mathématiques Pures et Appliqu\'. 2023 , doi =
2023
-
[13]
and Poiatti, Andrea , TITLE =
Gal, Ciprian G. and Poiatti, Andrea , TITLE =. European Journal of Applied Mathematics , VOLUME =. 2025 , NUMBER =
2025
-
[17]
Elbar, Charles and Skrzeczkowski, Jakub , journal =. Degenerate. 2023 , issn =. doi:10.1016/j.jde.2023.03.057 , publisher =
-
[19]
Mazón and Julio D
Fuensanta Andreu--Vaillo and José M. Mazón and Julio D. Rossi and J. Julián Toledo--Melero , publisher =. Nonlocal. 2010 , isbn =
2010
-
[20]
Global and Trajectory Attractors for a nonlocal
Frigeri, Sergio and Grasselli, Maurizio , journal =. Global and Trajectory Attractors for a nonlocal. 2012 , issn =. doi:10.1007/s10884-012-9272-3 , publisher =
-
[21]
Frigeri, Sergio and Grasselli, Maurizio , journal =. Nonlocal. 2012 , issn =. doi:10.4310/dpde.2012.v9.n4.a1 , publisher =
-
[22]
and Grasselli, Maurizio , journal =
Frigeri, Sergio and Gal, Ciprian G. and Grasselli, Maurizio , journal =. On Nonlocal. 2016 , issn =. doi:10.1007/s00332-016-9292-y , publisher =
-
[23]
Nonlinear evolution inclusions arising from phase change models , year =
Colli, Pierluigi and Krejčí, Pavel and Rocca, Elisabetta and Sprekels, Jürgen , journal =. Nonlinear evolution inclusions arising from phase change models , year =. doi:10.1007/s10587-007-0114-0 , publisher =
-
[24]
and Warma, Mahamadi , journal =
Gal, Ciprian G. and Warma, Mahamadi , journal =. Reaction--diffusion equations with fractional diffusion on non--smooth domains with various boundary conditions , year =. doi:10.3934/dcds.2016.36.1279 , publisher =
-
[25]
Ciprian G. Gal , journal =. Non-local. 2017 , issn =. doi:10.1017/s0956792516000504 , publisher =
-
[26]
Gal, Ciprian G. , journal =. Doubly nonlocal. 2018 , issn =. doi:10.1016/j.anihpc.2017.05.001 , publisher =
-
[29]
Emmanuele. Real. 2016 , isbn =. doi:10.1007/978-1-4939-4005-9 , issn =
-
[30]
Adams and John J
Robert A. Adams and John J. F. Fournier , publisher =. Sobolev. 2003 , edition =
2003
-
[31]
Bertozzi , journal =
Jacob Bedrossian and Nancy Rodríguez and Andrea L. Bertozzi , journal =. Local and global well-posedness for aggregation equations and. 2011 , number =
2011
-
[33]
, publisher =
Evans, Lawrence C. , publisher =. Partial. 2010 , edition =
2010
-
[38]
Fast, High-Order Accurate Integral Equation Methods and Application to
Hai Zhu , school =. Fast, High-Order Accurate Integral Equation Methods and Application to. 2021 , type =
2021
-
[41]
Martinsson, Per--Gunnar , publisher =. Fast. 2019 , isbn =
2019
-
[55]
Classical
Loukas Grafakos , publisher =. Classical. 2014 , edition =
2014
-
[56]
Theory and
Li, Changpin and Cai, Min , publisher =. Theory and. 2019 , isbn =
2019
-
[59]
, publisher =
Pozrikidis, C. , publisher =. The. 2016 , isbn =
2016
-
[67]
Panja, Madan Mohan and Mandal, Birendra Nath , publisher =. Wavelet. 2020 , isbn =
2020
-
[71]
Functional
Haïm Brezis , publisher =. Functional. 2011 , doi =
2011
-
[73]
Properties of Convolution Operators on
Thomas Potts , school =. Properties of Convolution Operators on. 2012 , type =
2012
-
[75]
Boyd , publisher =
John P. Boyd , publisher =. Chebyshev and. 2001 , edition =
2001
-
[82]
Chan, Tony F. and Mathew, Tarek P. , journal =. Domain decomposition algorithms , year =. doi:10.1017/s0962492900002427 , publisher =
-
[83]
Deville, M. O. and Fischer, P. F. and Mund, E. H. , publisher =. High-. 2002 , isbn =
2002
-
[84]
Spectral/hp
Karniadakis, George and Sherwin, Spencer , publisher =. Spectral/hp. 2005 , edition =
2005
-
[89]
Bueno--Orovio, Alfonso and Pérez--García, Víctor M. and Fenton, Flavio H. , journal =. Spectral Methods for Partial Differential Equations in Irregular Domains: The Spectral Smoothed Boundary Method , year =. doi:10.1137/040607575 , publisher =
-
[90]
and Mills--Williams, Rory D
Roden, Jonna C. and Mills--Williams, Rory D. and Pearson, John W. and Goddard, Benjamin D. , journal =
-
[91]
Roden and Rory Mills--Williams and Benjamin D
Jonna C. Roden and Rory Mills--Williams and Benjamin D. Goddard and John W. Pearson , howpublished =. Multi. 2022 , url =
2022
-
[92]
2017 , doi =
Goddard, Benjamin and Nold, Andreas and Kalliadasis, Serafim , title =. 2017 , doi =
2017
-
[98]
Burden and Douglas J
Richard L. Burden and Douglas J. Faires and Annette M. Burden , publisher =. Numerical. 2016 , edition =
2016
-
[99]
, publisher =
Trefethen, Lloyd N. , publisher =. Spectral. 2000 , isbn =
2000
-
[100]
Algorithms from
Lange, Kenneth , publisher =. Algorithms from. 2020 , isbn =
2020
-
[101]
and Gould, Nicholas I
Conn, Andrew R. and Gould, Nicholas I. M. and Toint, Philippe L. , publisher =. Trust. 2000 , isbn =
2000
-
[102]
Introductory
Kreyszig, Erwin , publisher =. Introductory. 1989 , isbn =
1989
-
[103]
Elliptic
Grisvard, Pierre , publisher =. Elliptic. 2011 , isbn =
2011
-
[105]
2023 , type =
Ivan Krylov , school =. 2023 , type =
2023
-
[107]
2025 , type =
Andrés Miniguano--Trujillo , school =. 2025 , type =
2025
-
[109]
Table of
Daniel Zwillinger and Victor. Table of. 2014 , isbn =
2014
-
[110]
Andrés Miniguano--Trujillo , journal =. The. 2025 , copyright =
2025
-
[111]
Convergence to equilibrium of weak solutions to the
Grasselli, Maurizio and Poiatti, Andrea , journal =. Convergence to equilibrium of weak solutions to the. 2025 , copyright =
2025
-
[113]
L. af Klinteberg , Singularity swap quadrature for nearly singular line integrals on closed curves in two dimensions , BIT Numerical Mathematics, 64 (2024), p. Art. 11, https://doi.org/10.1007/s10543-024-01013-0
-
[114]
M. Ainsworth and Z. Mao , Analysis and approximation of a fractional C ahn-- H illiard equation , SIAM Journal on Numerical Analysis, 55 (2017), pp. 1689--1718, https://doi.org/10.1137/16m1075302
-
[115]
H. Alzahrani, G. Turkiyyah, O. Knio, and D. Keyes , Space-fractional diffusion with variable order and diffusivity: D iscretization and direct solution strategies , Communications on Applied Mathematics and Computation, 4 (2022), pp. 1416--1440, https://doi.org/10.1007/s42967-021-00184-9
-
[116]
T. G. Anderson, H. Zhu, and S. Veerapaneni , A fast, high-order scheme for evaluating volume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings , Journal of Computational Physics, 472 (2023), p. 111688, https://doi.org/10.1016/j.jcp.2022.111688
-
[117]
H. Antil and S. Bartels , Spectral approximation of fractional PDE s in image processing and phase field modeling , Computational Methods in Applied Mathematics, 17 (2017), pp. 661--678, https://doi.org/10.1515/cmam-2017-0039
-
[118]
T. Askham and A. J. Cerfon , An adaptive fast multipole accelerated P oisson solver for complex geometries , Journal of Computational Physics, 344 (2017), pp. 1--22, https://doi.org/10.1016/j.jcp.2017.04.063
-
[119]
P. Assari and M. Dehghan , Application of dual- C hebyshev wavelets for the numerical solution of boundary integral equations with logarithmic singular kernels , Engineering with Computers, 35 (2019), pp. 175--190, https://doi.org/10.1007/s00366-018-0591-9
-
[120]
H. Bahouri, J.-Y. Chemin, and R. Danchin , Fourier A nalysis and N onlinear P artial D ifferential E quations , vol. 343 of Grundlehren der mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences], Springer Berlin, Heidelberg, 2011, https://doi.org/10.1007/978-3-642-16830-7
-
[121]
J. Bedrossian, N. Rodríguez, and A. L. Bertozzi , Local and global well-posedness for aggregation equations and P atlak-- K eller-- S egel models with degenerate diffusion , Nonlinearity, 24 (2011), pp. 1683--1714, https://doi.org/10.1088/0951-7715/24/6/001
-
[122]
C. Bernardi and Y. Maday , Theory and N umerics of D ifferential E quations , Springer Berlin, Heidelberg, 2001, ch. Spectral, S pectral E lement and M ortar E lement M ethods, pp. 1--57, https://doi.org/10.1007/978-3-662-04354-7_1
-
[123]
J. P. Boyd , Chebyshev and F ourier S pectral M ethods , Dover Publications, 2nd (revised) ed., 2001
2001
-
[124]
H. Brezis , Functional A nalysis, S obolev S paces and P artial D ifferential E quations , Springer New York, 2011, https://doi.org/10.1007/978-0-387-70914-7
-
[125]
H. Brezis and P. Mironescu , Gagliardo-- N irenberg inequalities and non-inequalities: The full story , Annales de l'Institut Henri Poincar\'e C. Analyse Non Lin\'eaire, 35 (2018), pp. 1355--1376, https://doi.org/10.1016/j.anihpc.2017.11.007
-
[126]
J. W. Cahn and J. E. Hilliard , Free energy of a nonuniform system. I . I nterfacial free energy , The Journal of Chemical Physics, 28 (1958), pp. 258--267, https://doi.org/10.1063/1.1744102
-
[127]
C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang , Spectral M ethods: F undamentals in S ingle D omains , Scientific Computation, Springer Berlin, Heidelberg, 2006, https://doi.org/10.1007/978-3-540-30726-6
-
[128]
C. Canuto, A. Quarteroni, M. Y. Hussaini, and T. A. Zang , Spectral M ethods: E volution to C omplex G eometries and A pplications to F luid D ynamics , Scientific Computation, Springer Berlin, Heidelberg, 2007, https://doi.org/10.1007/978-3-540-30728-0
-
[129]
P. Cavaliere and A. Cianchi , Classical and approximate T aylor expansions of weakly differentiable functions , Annales Academiae Scientiarum Fennic Mathematica, 39 (2014), pp. 527--544, https://doi.org/10.5186/aasfm.2014.3933
-
[130]
C. Cavaterra, S. Frigeri, and M. Grasselli , A nonlocal C ahn-- H illiard-- D arcy system with singular potential, degenerate mobility, and sources , Applied Mathematics & Optimization, 91 (2025), p. Art. 39, https://doi.org/10.1007/s00245-025-10239-5
-
[131]
L. Q. Chen and J. Shen , Applications of semi-implicit F ourier-spectral method to phase field equations , Computer Physics Communications, 108 (1998), pp. 147--158, https://doi.org/10.1016/s0010-4655(97)00115-x
-
[132]
X. Chen, X. Zhang, L. Wei, and L. Huang , A fast explicit time-splitting spectral scheme for the viscous C ahn-- H illiard equation with nonlocal diffusion operator , Applied Mathematics Letters, 157 (2024), p. 109169, https://doi.org/10.1016/j.aml.2024.109169
-
[133]
E. Davoli, L. Scarpa, and L. Trussardi , Local asymptotics for nonlocal convective C ahn- H illiard equations with W^ 1,1 kernel and singular potential , Journal of Differential Equations, 289 (2021), pp. 35--58, https://doi.org/10.1016/j.jde.2021.04.016
-
[134]
E. Davoli, L. Scarpa, and L. Trussardi , Nonlocal-to-local convergence of C ahn-- H illiard equations: N eumann boundary conditions and viscosity terms , Archive for Rational Mechanics and Analysis, 239 (2021), pp. 117--149, https://doi.org/10.1007/s00205-020-01573-9
-
[135]
F. Della Porta, A. Giorgini, and M. Grasselli , The nonlocal C ahn-- H illiard-- H ele-- S haw system with logarithmic potential , Nonlinearity, 31 (2018), pp. 4851--4881, https://doi.org/10.1088/1361-6544/aad52a
-
[136]
M. O. Deville, P. F. Fischer, and E. H. Mund , High- O rder M ethods for I ncompressible F luid F low , Cambridge University Press, 2002, https://doi.org/10.1017/cbo9780511546792
-
[137]
Q. Du, J. Yang, and Z. Zhou , Time-fractional A llen-- C ahn equations: Analysis and numerical methods , Journal of Scientific Computing, 85 (2020), p. Art. 42, https://doi.org/10.1007/s10915-020-01351-5
-
[138]
M. D’Elia, Q. Du, C. Glusa, M. Gunzburger, X. Tian, and Z. Zhou , Numerical methods for nonlocal and fractional models , Acta Numerica, 29 (2020), pp. 1--124, https://doi.org/10.1017/s096249292000001x
-
[139]
A. Ern and J.-L. Guermond , Finite E lements I : A pproximation and I nterpolation , vol. 72 of Texts in Applied Mathematics, Springer Cham, 2021, https://doi.org/10.1007/978-3-030-56341-7
-
[140]
L. C. Evans , Partial D ifferential E quations , vol. 19 of Graduate Studies in Mathematics, American Mathematical Society, 2nd ed., 2010
2010
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.