Proving the existence of localized patterns, periodic solutions, and branches of periodic solutions in the 1D Thomas model
Pith reviewed 2026-05-10 17:03 UTC · model grok-4.3
The pith
A framework using Newton-Kantorovich theorems and approximate inverses proves the existence of localized and periodic solutions in the 1D Thomas model.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Given an approximate solution bar u, the framework constructs an approximate inverse of the linearization and establishes conditions under which the Newton-Kantorovich fixed point map contracts on a ball around bar u. This yields rigorous proof of existence and local uniqueness for stationary localized solutions, spatially periodic solutions, and branches of periodic solutions in the 1D Thomas model, with explicit upper bounds computed for the required quantities despite the non-polynomial nature of the nonlinearity.
What carries the argument
The Newton-Kantorovich approach with a constructed approximate inverse of the linearized operator around an approximate solution bar u, which allows verification of contraction conditions for the fixed-point map.
If this is right
- Existence and local uniqueness of stationary localized solutions near the approximate solution.
- Existence and local uniqueness of spatially periodic solutions.
- Existence of branches of spatially periodic solutions.
- Control over the linearization around the approximate solution.
- Applicability to models with non-polynomial nonlinearities via specialized bounding techniques.
Where Pith is reading between the lines
- This method could extend to proving stability or dynamics of these solutions in time-dependent settings.
- Similar computer-assisted techniques might apply to higher-dimensional versions of the Thomas model or other reaction-diffusion equations.
- By making the bounds explicit and computable, the approach enables systematic exploration of parameter ranges where patterns exist.
Load-bearing premise
That a sufficiently accurate numerical approximate solution can be computed such that the constructed approximate inverse yields a contraction constant small enough to satisfy the Newton-Kantorovich conditions.
What would settle it
For a specific approximate solution bar u, if the computed upper bounds on the operator norms or the contraction constant exceed the thresholds required by the theorem, then the existence proof for solutions near that bar u fails.
Figures
read the original abstract
In this paper, we present a general framework for constructively proving the existence and of stationary localized solutions, spatially periodic solutions, and branches of spatially periodic solutions in the 1D Thomas model. Specifically, we develop the necessary analysis to compute explicit upper bounds required in a Newton--Kantorovich approach. Given an approximate solution $\bar{\mathbf{u}}$, this approach relies on establishing that a well-chosen fixed point map is contracting on a neighborhood $\bar{\mathbf{u}}$. For this matter, we construct an approximate inverse of the linearization around $\bar{\mathbf{u}}$, and establish sufficient conditions under which the contraction is achieved. This provides a framework for which computer-assisted analysis can be applied to verify the existence and local uniqueness of solutions in a vicinity of $\bar{\mathbf{u}}$, and control the linearization around $\bar{\mathbf{u}}$. Furthermore, as the Thomas model has a non-polynomial nonlinearity, we will need to use different techniques to handle it during our analysis. The code to perform the rigorous proofs is available on Github.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a general framework for constructively proving the existence of stationary localized solutions, spatially periodic solutions, and branches of spatially periodic solutions in the 1D Thomas model. It develops the analysis needed to compute explicit upper bounds in a Newton-Kantorovich approach: given a numerical approximation bar u, an approximate inverse to the linearization is constructed and sufficient conditions are derived to ensure the fixed-point map is contracting on a neighborhood of bar u. Special techniques are introduced to bound the non-polynomial nonlinearity, and the implementation code using interval arithmetic is made available on GitHub.
Significance. If the case-by-case validations hold, the work supplies a reproducible, computer-assisted method for rigorously confirming the existence and local uniqueness of patterns in a reaction-diffusion equation with transcendental nonlinearity. The emphasis on explicit bounds, the handling of the non-polynomial term, and the public code repository constitute clear strengths that enhance verifiability and potential reuse in the validated-numerics literature.
major comments (2)
- The central claim rests on the contraction constant being strictly less than one and the validated radius being positive for each chosen bar u. The manuscript should therefore include, in the results section that presents the first localized solution, the explicit numerical values of these quantities together with the truncation and rounding-error bounds that were used to obtain them.
- In the analysis of the non-polynomial nonlinearity, the Taylor remainder estimate must be shown to be compatible with the interval-arithmetic enclosure; a concrete inequality or lemma number should be cited that guarantees the remainder term does not destroy the contraction property.
minor comments (3)
- Abstract: the phrase 'existence and of stationary localized solutions' contains a typographical error and should read 'existence of stationary localized solutions'.
- Abstract: the statement 'control the linearization around bar u' is imprecise; a brief clarification of which spectral or norm properties are being controlled would improve readability.
- The manuscript should state explicitly in the introduction or methods section that all truncation and rounding errors arising from the numerical approximation and interval arithmetic are rigorously bounded, rather than leaving this implicit.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We are pleased that the referee recognizes the strengths of our framework and the availability of the code. Below we address each major comment point by point.
read point-by-point responses
-
Referee: The central claim rests on the contraction constant being strictly less than one and the validated radius being positive for each chosen bar u. The manuscript should therefore include, in the results section that presents the first localized solution, the explicit numerical values of these quantities together with the truncation and rounding-error bounds that were used to obtain them.
Authors: We agree that providing these explicit values will improve the transparency of our results. In the revised version of the manuscript, we will add to the results section, immediately following the presentation of the first localized solution, the specific numerical values for the contraction constant (strictly less than one), the positive validated radius, the truncation level, and the rounding-error bounds. These quantities are computed using our interval arithmetic implementation and will be reported directly from the code. revision: yes
-
Referee: In the analysis of the non-polynomial nonlinearity, the Taylor remainder estimate must be shown to be compatible with the interval-arithmetic enclosure; a concrete inequality or lemma number should be cited that guarantees the remainder term does not destroy the contraction property.
Authors: We appreciate this observation. Our analysis in the section on bounding the non-polynomial nonlinearity already employs a Taylor remainder estimate that is enclosed using interval arithmetic. To make the compatibility explicit as requested, we will cite the relevant inequality (or add a reference to a lemma in the appendix if appropriate) that ensures the remainder term's enclosure preserves the contraction property. This clarification will be incorporated in the revised manuscript. revision: yes
Circularity Check
No significant circularity
full rationale
The paper constructs a Newton-Kantorovich fixed-point argument around a given numerical approximate solution bar u, deriving explicit bounds on the approximate inverse and contraction constant to certify existence and local uniqueness in a ball. These bounds are computed from the linearization at bar u and interval arithmetic on the non-polynomial nonlinearity; the target solution is not defined in terms of the bounds, nor is any prediction fitted to data and then re-labeled. No self-citation chain, imported uniqueness theorem, or ansatz smuggling appears in the derivation. The central claim remains conditional on case-by-case verification that the contraction constant is strictly less than one, which is an external numerical check rather than an internal reduction.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Banach fixed-point theorem in a complete metric space
- domain assumption Existence of a sufficiently accurate numerical approximate solution bar u
Reference graph
Works this paper leans on
-
[1]
F. Al Saadi, A. Champneys, and N. Verschueren. Localized patt erns and semi-strong inter- action, a unifying framework for reaction–diffusion systems. IMA Journal of Applied Mathe- matics, 86(5):1031–1065, 2021
work page 2021
-
[2]
F. Al Saadi, A. Worthy, H. Alriheieli, and M. Nelson. Localised spatia l structures in the Thomas model. Mathematics and Computers in Simulation , 194:141–158, 2022
work page 2022
-
[3]
G. Arioli and H. Koch. Computer-Assisted Methods for the Stud y of Stationary Solutions in Dissipative Systems, Applied to the Kuramoto-Sivashinski Equatio n. Archive for Rational Mechanics and Analysis , 197:1033–1051, 2010
work page 2010
- [4]
-
[5]
L. Benet and D. Sanders. Intervalarithmetic.jl, 2022. https://github.com/JuliaIntervals/IntervalArithmetic.jl
work page 2022
-
[6]
D. Blanco. D3D6solutionSH.jl, 2026. https://github.com/dominicblanco/D3D6solutionsSH.jl
work page 2026
- [7]
-
[8]
D. Blanco. ThomasProofs.jl, 2026. https://https://github.com/dominicblanco/ThomasProofs.jl
work page 2026
-
[9]
D. Blanco and M. Cadiot. Localizedpatternsshsymmetry.jl, 2025 . https://github.com/dominicblanco/LocalizedPatternsSHsymmetry.jl
work page 2025
-
[10]
D. Blanco and M. Cadiot. Proving symmetry of localized solutions a nd application to dihedral patterns in the planar Swift-Hohenberg PDE. arXiv:2509.10375, 2025
- [11]
-
[12]
D. Blanco, M. Cadiot, and D. Fassler. Proving the existence of lo calized patterns and saddle node bifurcations in 1D activator-inhibitor type models. arXiv:2509.17099, 2025
work page internal anchor Pith review arXiv 2025
-
[13]
G. M. Boissoniere, R. Choksi, and J.-P. Lessard. Microscopic pa tterns in the 2D phase-field crystal model. Nonlinearity, 35:1500–1520, 2022
work page 2022
-
[14]
M. Breden. Computer-assisted proofs for some nonlinear diffu sion problems. Communications in Nonlinear Science and Numerical Simulation , 109:106292, 2022
work page 2022
-
[15]
M. Breden. A Posteriori Validation of Generalized Polynomial Cha os Expansions. SIAM Journal on Applied Dynamical Systems , 22(2):765–801, 2023
work page 2023
-
[16]
M. Breden. Computer-assisted proofs for differential equat ions and dynamical systems, 2025. 39
work page 2025
-
[17]
M. Breden and H. Chu. Constructive proofs for some semilinear PDEs on H 2(e|x|2/ 4, Rd). Numerische Mathematik , pages 2097–2143, 2025
work page 2097
-
[18]
M. Breden and H. Chu. Solutions of differential equations in Freu d-weighted Sobolev spaces, 2025
work page 2025
- [19]
-
[20]
M. Breden, J.-P. Lessard, and M. Vanicat. Global bifurcation d iagrams of steady states of systems of PDEs via rigorous numerics: a 3-component reaction -diffusion system. Acta Applicandae Mathematicae, 128(1):113–152, 2013
work page 2013
-
[21]
M. Breden and M. Payan. Computer-assisted proofs for the m any steady states of a chemotaxis model with local sensing. Physica D: Nonlinear Phenomena , 466:134221, 2024
work page 2024
-
[22]
X. Cabr´ e, E. Fontich, and R. de la Llave. The parameterizationmethod for invariant manifolds. I. Manifolds associated to non-resonant subspaces. Indiana Univ. Math. J. , 52(2):283–328, 2003
work page 2003
-
[23]
X. Cabr´ e, E. Fontich, and R. de la Llave. The parameterizationmethod for invariant manifolds. II. Regularity with respect to parameters. Indiana Univ. Math. J. , 52(2):329–360, 2003
work page 2003
-
[24]
X. Cabr´ e, E. Fontich, and R. de la Llave. The parameterizationmethod for invariant manifolds. III. Overview and applications. J. Differential Equations , 218(2):444–515, 2005
work page 2005
-
[25]
M. Cadiot. Proofkawahara.jl, 2023. https://github.com/matthieucadiot/ProofKawahara.jl
work page 2023
-
[26]
M. Cadiot. Localizedpatternsh.jl, 2024. https://github.com/matthieucadiot/LocalizedPatternSH.jl
work page 2024
-
[27]
M. Cadiot. Constructive proofs of existence and stability of so litary waves in the Whitham and capillary–gravity Whitham equations. Nonlinearity, 38(3):035021, 2025
work page 2025
-
[28]
M. Cadiot and D. Blanco. Localizedpatterngs.jl, 2024. https://github.com/dominicblanco/LocalizedPatternGS.jl
work page 2024
-
[29]
M. Cadiot and D. Blanco. Localized stationary patterns in the 2d gray scott model: computer- assisted proofs of existence. Nonlinearity, 38(4):045016, 2025
work page 2025
-
[30]
M. Cadiot, J.-P. Lessard, and J.-C. Nave. Rigorous Computatio n of Solutions of Semilinear PDEs on Unbounded Domains via Spectral Methods. SIAM Journal on Applied Dynamical Systems, 23(3):1966–2017, 2024
work page 1966
-
[31]
M. Cadiot, J.-P. Lessard, and J.-C. Nave. Stationary non-rad ial localized patterns in the planar Swift-Hohenberg PDE: Constructive proofs of existence. Journal of Differential Equations , 414:555–608, 2025
work page 2025
-
[32]
R. Calleja, C. Garcia-Azpeitia, O. Henot, J.-P. Lessard, and J. Mireles-James. From the Lagrange Triangle to the Figure Eight Choreography: Proof of Marchal’s Conjecture. arXiv:2406.17564, 2024
-
[33]
A. Champneys, F. Al Saadi, V. F. Bre˜ na–Medina, V. A. Grieneis en, A. F. Mar´ ee, N. Ver- schueren, and B. Wuyts. Bistability, wave pinning and localisation in na tural reaction–diffusion systems. Physica D: Nonlinear Phenomena , 416:132735, 2021
work page 2021
-
[34]
K. E. Church, J.-Y. Dai, O. Henot, and P. Lappicy. Global Contin uation of Stable Periodic Orbits in Systems of Competing Predators. arXiv:2504.03058v2, 2026
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[35]
J. Cyranka and J.-P. Lessard. Validated forward integration s cheme for parabolic PDEs via Chebyshev series. Commun. Nonlinear Sci. Numer. Simul. , 109:Paper No. 106304, 32, 2022
work page 2022
-
[36]
J. Cyranka and P. B. Mucha. A construction of two different so lutions to an elliptic system. Journal of Mathematical Analysis and Applications , 465(1):500–530, 2018. 40
work page 2018
- [37]
-
[38]
J.-L. Figueras and R. de la Llave. Numerical computations and co mputer assisted proofs of periodic orbits of the Kuramoto-Sivashinsky equation. SIAM J. Appl. Dyn. Syst. , 16(2):834– 852, 2017
work page 2017
-
[39]
J.-L. Figueras, A. Haro, and A. Luque. Rigorous computer ass isted application of KAM theory: a modern approach. Foundations of Computational Mathematics , 17:1123–1193, 2017
work page 2017
-
[40]
M. Gameiro and J.-P. Lessard. Analytic estimates and rigorous c ontinuation for equilibria of higher-dimensional PDEs. Journal of Differential Equations , 249(9):2237–2268, 2010
work page 2010
-
[41]
M. Gameiro and J.-P. Lessard. Efficient rigorous numerics for hig her-dimensional PDEs via one-dimensional estimates. SIAM Journal on Numerical Analysis , 51(4):2063–2087, 2013
work page 2063
-
[42]
M. Gameiro and J.-P. Lessard. A posteriori verification of invar iant objects of evolution equations: Periodic orbits in the Kuramoto–Sivashinsky PDE. SIAM Journal on Applied Dynamical Systems , 16(1):687–728, 2017
work page 2017
-
[43]
P. Gray and S. K. Scott. Sustained oscillations and other exotic patterns of behavior in isothermal reactions. Journal of Physical Chemistry , 89(1):22–32, 1985
work page 1985
- [44]
-
[45]
O. Henot. On polynomial forms of nonlinear functional different ial equations. Journal of Computational Dynamics , 8:307–323, 2021
work page 2021
-
[46]
J. Hervagault, A. Friboulet, J. Kernevez, and D. Thomas. Spa tiotemporal behaviors in im- mobilized enzyme systems. Biochimie, 62:367–373, 1980
work page 1980
-
[47]
O. H´ enot. Radiipolynomial.jl, 2022.https://github.com/OlivierHnt/RadiiPolynomial.jl
work page 2022
-
[48]
J. Kernevez, E. Doedel, M. Duban, J. Hervagault, G. Joly, and D. Thomas. Spatio-Temporal Organization in Immobilized Enzyme Systems. Rhythms in Biology and Other Fields of Ap- plication, 49:50–75, 1983
work page 1983
-
[49]
J. Kernevez, G. Joly, M. Duban, B. Bunow, and D. Thomas. Hys teresis, oscillations, and pattern formation in realistic immobilized enzyme systems. Journal of Mathematical Biology , 7:41–56, 1979
work page 1979
-
[50]
J.-P. Lessard. Rigorous verification of saddle-node bifurcatio ns in ODEs. Indagationes Math- ematicae, 27:1013–1026, 2016
work page 2016
-
[51]
J.-P. Lessard. Continuation of solutions and studying delay diffe rential equations via rigorous numerics. Proceedings of Symposia in Applied Mathematics , pages 81–122, 2018
work page 2018
-
[52]
J.-P. Lessard, J. Mireles-James, and J. Ransford. Automatic differentiation for Fourier series and the radii polynomial approach. Physica D: Nonlinear Phenomena , 334:174–186, 2016
work page 2016
-
[53]
J.-P. Lessard and A. Pugliese. Cusp bifurcations: Numerical de tection via two-parameter continuation and computer-assisted proofs of existence. Discrete and Continuous Dynamical Systems - B , 30(6):2135–2158, 2025
work page 2025
-
[54]
J.-P. Lessard, E. Sander, and T. Wanner. Rigorous continuat ion of bifurcation points in the diblock copolymer equation. Journal of Computational Dynamics , 4(1&2):71–118, 2017
work page 2017
- [55]
-
[56]
J. Murray. On pattern formation mechanisms for lepidopteran wing patterns and mammalian coat markings. Philosophical Transactions B , 295:473–496, 1981. 41
work page 1981
-
[57]
J. Murray. A pre-pattern formation mechanism for animal coa t markings. Journal of Theo- retical Biology, 88:161–199, 1981
work page 1981
-
[58]
J. D. Murray. Mathematical Biology, volume 1. Springer Nature Link, 2003
work page 2003
-
[59]
J. D. Murray. Mathematical Biology II , volume 2. Springer Nature Link, 2003
work page 2003
-
[60]
M. T. Nakao, M. Plum, and Y. Watanabe. Numerical verification methods and computer- assisted proofs for partial differential equations . Springer, 2019
work page 2019
-
[61]
R. Payne and C. Grierson. A theoretical model for rop localisat ion by auxin in arabidopsis root hair cells. PLoS ONE , 2009
work page 2009
-
[62]
I. Prigogine and R. Lefever. Symmetry breaking instabilities in dis sipative systems. ii. The Journal of Chemical Physics , 1968
work page 1968
-
[63]
E. Sander and T. Wanner. Pattern formation in a nonlinear mode l for animal coats. Journal of Differential Equations , 191:143–174, 2003
work page 2003
-
[64]
E. Sander and T. Wanner. Equilibrium validation in models for patte rn formation based on Sobolev embeddings. Discrete and Continuous Dynamical Systems - B , 26:603–632, 2020
work page 2020
-
[65]
E. Sander and T. Wanner. Equilibrium validation in models for patte rn formation based on Sobolev embeddings. Discrete and Continuous Dynamical Systems - B , 26(1):603–632, 2021
work page 2021
-
[66]
J. Schnakenberg. Simple chemical reaction systems with limit cyc le behaviour. Journal of Theoretical Biology, 81(3):389–400, 1979
work page 1979
-
[67]
E. E. Sel’kov. Self-oscillations in glycolysis 1. a simple kinetic model. European Journal of Biochemistry, 4(1):79–86, 1968
work page 1968
- [68]
- [69]
-
[70]
J. B. van den Berg, M. Breden, J.-P. Lessard, J. Mireles-Jame s, and K. Mischaikow. Rigorous numerics for symmetric connecting orbits: even homoclinics of the G ray-Scott equation. SIAM Journal on Mathematical Analysis , 43(4):1557–1594, 2011
work page 2011
-
[71]
J. B. van den Berg, M. Breden, J.-P. Lessard, and M. Murray. Continuation of homoclinic orbits in the suspension bridge equation: A computer-assisted pro of. Journal of Differential Equations, 264(5):3086–3130, 2018
work page 2018
-
[72]
J. B. van den Berg, M. Breden, J.-P. Lessard, and M. Murray. Continuation of homoclinic orbits in the suspension bridge equation: A computer-assisted pro of. Journal of Differential Equations, 264:3086–3130, 2018
work page 2018
-
[73]
J. B. van den Berg, A. Deschenes, J.-P. Lessard, and J. Mirele s-James. Stationary coexistence of hexagons and rolls via rigorous computations. SIADS, 2:942–979, 2015
work page 2015
-
[74]
J. B. van den Berg, O. H´ enot, and J.-P. Lessard. Construct ive proofs for localized radial solutions of semilinear elliptic systems on Rd. Nonlinearity, 36(12):6476–6512, 2023
work page 2023
-
[75]
J.-B. van den Berg, J.-P. Lessard, and K. Mischaikow. Global sm ooth solution curves using rigorous branch following. Mathematics of Computation , 79(271):1565–1584, 2010
work page 2010
-
[76]
J.-B. van den Berg and J. Williams. Rigorously Computing Symmetric Stationary Sates of the Ohta-Kawasaki Problem in Three Dimensions. SIAM J. Math Anal. , 2019. 42
work page 2019
-
[77]
J.-B. van den Berg and J. Williams. Optimal periodic structures wit h general space group symmetries in the Ohta-Kawasaki problem. Physica D (Nonlinear Phenomena) , 2021
work page 2021
-
[78]
L. van der Aalst and M. Cadiot. Existence proofs of traveling wa ve solutions on an infinite strip for the suspension bridge equation and proof of orbital st. arXiv:2509.16693, 2026
work page internal anchor Pith review arXiv 2026
-
[79]
L. van der Aalst, J. B. van den Berg, and J.-P. Lessard. Period ic localised traveling waves in the two-dimensional suspension bridge equation. Nonlinearity, 38(7):075029, 2025
work page 2025
-
[80]
J. M. Wunderlich. Computer-assisted Existence Proofs for Navier-Stokes Equ ations on an Unbounded Strip with Obstacle . PhD thesis, Karlsruher Institut f¨ ur Technologie (KIT), 2022
work page 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.