Global bifurcation of localised 2D patterns emerging from spatial heterogeneity
Pith reviewed 2026-05-18 11:11 UTC · model grok-4.3
The pith
Compact radial heterogeneity in the 2D Swift-Hohenberg equation produces branches of fully localised patterns that alternate between spots and dipoles.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors prove the existence of local bifurcation branches of fully localised patterns in the two-dimensional Swift-Hohenberg equation whose linear bifurcation parameter is perturbed by a radially-symmetric potential function with compact support. The trivial state is unstable in a compact neighbourhood of the origin and linearly stable outside. They characterise stability and bifurcation structure along these branches and then continue the solutions rigorously to large amplitude via analytic global bifurcation theory. The primary bifurcating branch alternates between an axisymmetric spot and a non-axisymmetric dipole pattern depending on the width of the spatial heterogeneity.
What carries the argument
Analytic global bifurcation theory applied within function spaces adapted to the compact radially symmetric heterogeneity, which allows local bifurcation branches of localised solutions to be continued globally in amplitude.
If this is right
- Fully localised multi-dimensional patterns exist near a Turing instability when a compact spatial heterogeneity is present.
- The primary bifurcation branch switches symmetry type with changes in the width of the heterogeneity.
- Solutions on the branches can be continued to large amplitudes while remaining localised.
- Stability properties of the localised patterns can be tracked along the bifurcating branches.
Where Pith is reading between the lines
- The same function-space setup and global bifurcation arguments could apply to other reaction-diffusion systems that possess compact inhomogeneities in two or more dimensions.
- The observed alternation between spot and dipole patterns may control which localised structures are selected in experiments where the size of an inhomogeneity can be tuned.
- Combining the analytic continuation with numerical path-following methods would allow quantitative checks of the large-amplitude behaviour.
Load-bearing premise
The perturbation to the linear bifurcation parameter must be radially symmetric with compact support so the trivial state is linearly unstable only inside a bounded region and stable outside it.
What would settle it
A numerical continuation that finds no global branch of localised solutions or that shows no switch between axisymmetric spots and non-axisymmetric dipoles when the heterogeneity radius is varied would falsify the global continuation and primary-branch claims.
Figures
read the original abstract
We present a general approach to prove the existence, both locally and globally in amplitude, of fully localised multi-dimensional patterns in partial differential equations containing a compact spatial heterogeneity. While one-dimensional localised patterns induced by spatial heterogeneities have been well-studied, proving the existence of fully localised patterns emerging from a Turing instability in higher dimensions remains a key open problem in pattern formation. In order to demonstrate the approach, we consider the two-dimensional Swift--Hohenberg equation, whose linear bifurcation parameter is perturbed by a radially-symmetric potential function. In this case, the trivial state is unstable in a compact neighbourhood of the origin and linearly stable outside. We prove the existence of local bifurcation branches of fully localised patterns, characterise their stability and bifurcation structure, and then rigorously continue solutions to large amplitude via analytic global bifurcation theory. Notably, the primary bifurcating branch in the Swift--Hohenberg equation alternates between an axisymmetric spot and a non-axisymmetric `dipole' pattern, depending on the width of the spatial heterogeneity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a general approach to prove existence of fully localised multi-dimensional patterns in PDEs with compact spatial heterogeneities. Applied to the 2D Swift-Hohenberg equation with a radially symmetric compact perturbation to the linear bifurcation parameter (rendering the trivial state unstable only inside a ball and stable outside), the authors establish local bifurcations of localised patterns, characterise their stability and bifurcation structure, and apply analytic global bifurcation theory to continue branches to large amplitude. A key observation is that the primary branch alternates between axisymmetric spots and non-axisymmetric dipole patterns depending on the heterogeneity width.
Significance. If the proofs hold, this addresses a longstanding open problem in pattern formation by providing the first rigorous local-to-global existence results for fully localised 2D Turing patterns induced by compact heterogeneity. The technical advance lies in adapting analytic global bifurcation theorems to weighted spaces for a 2D operator whose essential spectrum touches the imaginary axis, with the compact radial perturbation used to isolate discrete eigenvalues. The alternation between spot and dipole patterns on the primary branch is a concrete, falsifiable prediction that strengthens the work. The approach could template similar analyses in other heterogeneous systems.
major comments (2)
- §2.3 (Function space setup) and §5 (Global continuation): The central claim relies on the linearised operator being Fredholm with the essential spectrum strictly to the left of the imaginary axis uniformly along the branch in the chosen weighted spaces. The compact radial potential is invoked to localise the instability, but the manuscript must supply an explicit estimate (or lemma) showing that the 2D dispersion relation of the unperturbed Swift-Hohenberg operator is shifted sufficiently far left for all amplitudes on the global branch; without this, the analytic global bifurcation theorem cannot be applied directly as branches may accumulate on the essential spectrum rather than continuing to large amplitude while remaining localised.
- Theorem 4.2 (Local bifurcation) and the primary branch characterisation: The claim that the primary branch alternates between axisymmetric and dipole patterns depending on heterogeneity width is load-bearing for the novelty of the result. The proof should explicitly verify the crossing direction and the dimension of the kernel for each width regime, including confirmation that the transversality condition holds after the radial perturbation; the current sketch leaves open whether higher-order terms or the specific form of the potential could alter the alternation for certain parameter values.
minor comments (3)
- [Abstract and §4] The abstract states that stability is characterised, yet the main text would benefit from a dedicated subsection or corollary summarising the stability conclusions for the local branches (e.g., which segments are stable).
- [§2] Notation for the weighted Sobolev spaces (e.g., the precise form of the weight function) should be introduced earlier and used consistently; occasional switches between H^2_w and other spaces obscure the Fredholm property arguments.
- [Figures] Figure captions for the dipole patterns should include the specific heterogeneity width value used and a brief note on how the pattern was obtained (analytic vs. numerical continuation).
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for highlighting these important technical points. We address each major comment below and will incorporate clarifications and additions in the revised version.
read point-by-point responses
-
Referee: §2.3 (Function space setup) and §5 (Global continuation): The central claim relies on the linearised operator being Fredholm with the essential spectrum strictly to the left of the imaginary axis uniformly along the branch in the chosen weighted spaces. The compact radial potential is invoked to localise the instability, but the manuscript must supply an explicit estimate (or lemma) showing that the 2D dispersion relation of the unperturbed Swift-Hohenberg operator is shifted sufficiently far left for all amplitudes on the global branch; without this, the analytic global bifurcation theorem cannot be applied directly as branches may accumulate on the essential spectrum rather than continuing to large amplitude while remaining localised.
Authors: We agree that a more explicit estimate strengthens the application of the analytic global bifurcation theorem. The manuscript already selects weighted spaces in which the essential spectrum of the unperturbed operator lies strictly to the left of the imaginary axis, with the compact radial heterogeneity ensuring the perturbation remains relatively compact. In the revision we will add an explicit lemma in §5 deriving a uniform lower bound on the distance of the essential spectrum from the imaginary axis for all solutions along the branch with amplitude bounded by any fixed large constant, using the exponential decay guaranteed by the weighted norm. Continuation to arbitrarily large amplitudes then follows from the analyticity of the branch and a priori localisation estimates derived from the equation. revision: yes
-
Referee: Theorem 4.2 (Local bifurcation) and the primary branch characterisation: The claim that the primary branch alternates between axisymmetric and dipole patterns depending on heterogeneity width is load-bearing for the novelty of the result. The proof should explicitly verify the crossing direction and the dimension of the kernel for each width regime, including confirmation that the transversality condition holds after the radial perturbation; the current sketch leaves open whether higher-order terms or the specific form of the potential could alter the alternation for certain parameter values.
Authors: The alternation result in Theorem 4.2 follows from the radial symmetry of the heterogeneity and the ordering of the unperturbed eigenvalues under the 2D dispersion relation. We will expand the proof to include explicit verification of the crossing directions and kernel dimensions for the two width regimes: for narrow heterogeneity the axisymmetric mode crosses first with simple kernel, while for wider heterogeneity the first dipole pair crosses first. The transversality condition is preserved because the radial perturbation is a compact, symmetry-preserving operator that shifts the eigenvalues continuously without destroying simplicity at the critical values; this is confirmed by the Fredholm index and the cubic nonlinearity determining the leading-order bifurcation direction. Higher-order terms do not reverse the ordering for the parameter ranges considered. revision: yes
Circularity Check
No circularity: standard global bifurcation theorems applied to explicitly defined perturbed operator
full rationale
The derivation begins from the perturbed Swift-Hohenberg equation with a given compactly supported radial potential that defines the linear operator's spectrum (unstable inside a ball, stable outside). Local bifurcation branches are obtained from the standard Crandall-Rabinowitz theorem applied to this operator at the critical eigenvalue; global continuation follows from analytic versions of Rabinowitz/Dancer theorems in the weighted spaces constructed from the same compact support assumption. Neither step reduces to a fitted parameter renamed as prediction, nor to a self-citation chain, nor to an ansatz smuggled from prior work by the same authors. The radial symmetry is an explicit modeling choice used to equip the function spaces, not a result derived from the patterns themselves. The paper is therefore self-contained against external mathematical benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Standard analytic global bifurcation theorems apply to the nonlinear operator obtained after the radial perturbation of the Swift-Hohenberg equation.
- domain assumption The heterogeneity is radially symmetric and compactly supported, making the trivial state unstable only inside a bounded region.
Reference graph
Works this paper leans on
- [1]
-
[2]
M. Avery and A. Scheel. Sharp decay rates for localized perturbations to the critical front in the Ginzburg-Landau equation.J. Dynam. Differ. Equations, 36:S287–S322, 2024
work page 2024
-
[3]
D. Avitabile and H. Schmidt. Snakes and ladders in an inhomogeneous neural field model.Phys. D, 294:24–36, 2015
work page 2015
-
[4]
J. Baltaev and M. Kuˇ cera. Global bifurcation for quasivariational inequalities of reaction–diffusion type. J. Math. Anal. Appl., 345(2):917–928, 2008
work page 2008
-
[5]
R. Bastiaansen, M. Chirilus-Bruckner, and A. Doelman. Pulse solutions for an extended Klausmeier model with spatially varying coefficients.SIAM J. Appl. Dyn. Syst., 19(1):1–57, 2020
work page 2020
-
[6]
H. Berestycki and P.-L. Lions. Nonlinear Scalar Field Equations, I Existence of a Ground State.Arch. Rat. Mech. Anal., 82(4):313–345, Dec 1983
work page 1983
-
[7]
J. J. Bramburger, D. J. Hill, and D. J. B. Lloyd. Localized patterns.SIAM Review, 2025
work page 2025
-
[8]
P. C. Bressloff and Z. P. Kilpatrick. Two-dimensional bumps in piecewise smooth neural fields with synaptic depression.SIAM J. Appl. Math., 71(2):379–408, 2011
work page 2011
- [9]
- [10]
-
[11]
B. Buffoni and J. Toland.Analytic Theory of Global Bifurcation. Princeton University Press, 2003
work page 2003
- [12]
-
[13]
S. B¨ ohmer, B. Hilder, and J. Jansen. Pattern formation and film rupture in a two-dimensional thermo- capillary thin-film model of the B´ enard-Marangoni problem, 2025. arXiv preprint
work page 2025
-
[14]
M. Cadiot, J.-P. Lessard, and J.-C. Nave. Stationary non-radial localized patterns in the planar Swift- Hohenberg PDE: constructive proofs of existence.J. Differ. Equations, 414:555–608, 2025
work page 2025
-
[15]
D. Casta˜ no, M. Navarro, and H. Herrero. Routes to chaos from axisymmetric vertical vortices in a rotating cylinder.Appl. Math. Model., 54:1–20, 2018
work page 2018
-
[16]
K. Chen, Z. Deiman, R. Goh, S. Jankovic, and A. Scheel. Strain and defects in oblique stripe growth. Multiscale Model. Simul., 19(3):1236–1260, 2021
work page 2021
-
[17]
R. M. Chen, S. Walsh, and M. H. Wheeler. Global bifurcation for monotone fronts of elliptic equations. J. Eur. Math. Soc., 2024
work page 2024
-
[18]
A. Constantin, W. Strauss, and E. V˘ arv˘ aruc˘ a. Global bifurcation of steady gravity water waves with critical layers.Acta Math., 217(2):195–262, 2016
work page 2016
-
[19]
M. G. Crandall and P. H. Rabinowitz. Bifurcation from simple eigenvalues.J. Funct. Anal., 8(2):321– 340, 1971. 23
work page 1971
-
[20]
E. N. Dancer. Global structure of the solutions of non-linear real analytic eigenvalue problems.Proc. London Math. Soc., s3-27(4):747–765, 1973
work page 1973
- [21]
-
[22]
A. Doelman, P. van Heijster, and J. Shen. Pulse dynamics in reaction-diffusion equations with strong spatially localized impurities.Philos. Trans. Roy. Soc. A, 376(2117):20170183, 20, 2018
work page 2018
-
[23]
A. Doelman, P. van Heijster, and F. Xie. A geometric approach to stationary defect solutions in one space dimension.SIAM J. Appl. Dyn. Syst., 15(2):655–712, 2016
work page 2016
-
[24]
H. Fujii and Y. Nishiura. Global bifurcation diagram in nonlinear diffusion systems. In H. Fujita, P. D. Lax, and G. Strang, editors,Nonlinear Partial Differential Equations in Applied Science; Proceedings of The U.S.-Japan Seminar, Tokyo, 1982, volume 81 ofNorth-Holland Mathematics Studies, pages 17–35. North-Holland, 1983
work page 1982
-
[25]
C. Garcia-Azpeitia, Z. Ghanem, and W. Krawcewicz. Global bifurcation of spiral wave solutions to the complex Ginzburg–Landau equation, 2025. arXiv preprint
work page 2025
-
[26]
R. Goh, T. J. Kaper, and A. Scheel. Pitchfork bifurcation along a slow parameter ramp: coherent structures in the critical scaling.Stud. Appl. Math., 153(2):Paper No. e12702, 21, 2024
work page 2024
- [27]
- [28]
-
[29]
M. D. Groves and D. J. Hill. On function spaces for radial functions, 2024. arXiv preprint
work page 2024
-
[30]
M. D. Groves and J. Horn. Small-amplitude static periodic patterns at a fluid–ferrofluid interface.Proc. R. Soc. A, 474(2216):20180038, 2018
work page 2018
-
[31]
S. V. Haziot, V. M. Hur, W. A. Strauss, J. F. Toland, E. Wahl´ en, S. Walsh, and M. H. Wheeler. Traveling water waves—the ebb and flow of two centuries.Quart. Appl. Math., 80(2):317–401, 2022
work page 2022
-
[32]
S. V. Haziot and M. H. Wheeler. Large-amplitude steady solitary water waves with constant vorticity. Arch. Rat. Mech. Anal., 247(2):27, Mar 2023
work page 2023
-
[33]
D. J. Hill. Predicting the emergence of localised dihedral patterns in models for dryland vegetation.J. Nonlinear Sci., 34(4):67, May 2024
work page 2024
-
[34]
D. J. Hill, J. J. Bramburger, and D. J. B. Lloyd. Approximate localised dihedral patterns near a Turing instability.Nonlinearity, 36(5):2567–2630, 2023
work page 2023
-
[35]
D. J. Hill, J. J. Bramburger, and D. J. B. Lloyd. Dihedral rings of patterns emerging from a Turing bifurcation.Nonlinearity, 37(3):035015, 39, 2024
work page 2024
-
[36]
D. J. Hill and D. J. B. Lloyd. Radial amplitude equations for fully localized planar patterns.SIAM J. Appl. Math., 84(6):2590–2611, 2024
work page 2024
-
[37]
D. J. Hill, D. J. B. Lloyd, and M. R. Turner. Localised radial patterns on the surface of a ferrofluid.J. Nonlinear Sci., 31(79), 2021
work page 2021
-
[38]
D. J. Hill, D. J. B. Lloyd, and M. R. Turner.Spatially_Varying_SH, 2025. Github repository, https://github.com/Dan-Hill95/Spatially_Varying_SH. 24
work page 2025
-
[39]
O. Ja¨ ıbi, A. Doelman, M. Chirilus-Bruckner, and E. Meron. The existence of localized vegetation patterns in a systematically reduced model for dryland vegetation.Phys. D, 412:132637, 2020
work page 2020
- [40]
-
[41]
G. Jaramillo and A. Scheel. Deformation of striped patterns by inhomogeneities.Math. Methods Appl. Sci., 38(1):51–65, 2015
work page 2015
-
[42]
G. Jaramillo, A. Scheel, and Q. Wu. The effect of impurities on striped phases.Proc. Roy. Soc. Edinburgh Sect. A, 149(1):131–168, 2019
work page 2019
-
[43]
G. Jaramillo and S. C. Venkataramani. Target patterns in a 2D array of oscillators with nonlocal coupling.Nonlinearity, 31(9):4162–4201, 2018
work page 2018
-
[44]
J. Kamphuis and M. Chirilus-Bruckner. Pattern formation in a Swift-Hohenberg equation with spatially periodic coefficients, 2025. arXiv preprint
work page 2025
-
[45]
H.-C. Kao, C. Beaume, and E. Knobloch. Spatial localization in heterogeneous systems.Phys. Rev. E, 89:012903, Jan 2014
work page 2014
-
[46]
T. Kapitula and K. Promislow.Spectral and dynamical stability of nonlinear waves, volume 457. Springer, 2013
work page 2013
-
[47]
Kato.Perturbation Theory for Linear Operators
T. Kato.Perturbation Theory for Linear Operators. Classics in Mathematics. Springer Berlin Heidelberg, Berlin, Heidelberg, 1995
work page 1995
-
[48]
C. J. K. Knight, G. Derks, A. Doelman, and H. Susanto. Stability of stationary fronts in a non-linear wave equation with spatial inhomogeneity.J. Differ. Equations, 254(2):408–468, 2013
work page 2013
-
[49]
V. Kozlov. On non-uniqueness of solitary waves on two-dimensional rotational flow, 2025
work page 2025
-
[50]
D. J. B. Lloyd, C. Gollwitzer, I. Rehberg, and R. Richter. Homoclinic snaking near the surface instability of a polarisable fluid.J. Fluid Mech., 783:283–305, 2015
work page 2015
-
[51]
G. Martine-La Boissoni` ere, R. Choksi, and J.-P. Lessard. Microscopic patterns in the 2D phase-field- crystal model.Nonlinearity, 35(3):1500–1520, 2022
work page 2022
-
[52]
Y. Mau, L. Haim, A. Hagberg, and E. Meron. Competing resonances in spatially forced pattern-forming systems.Phys. Rev. E, 88:032917, Sep 2013
work page 2013
-
[53]
M. Navarro, D. Casta˜ no, and H. Herrero. Thermoconvective instabilities to explain the main charac- teristics of a dust devil-like vortex.Phys. D, 308:109–115, 2015
work page 2015
- [54]
-
[55]
M. R. Owen, C. R. Laing, and S. Coombes. Bumps and rings in a two-dimensional neural field: splitting and rotational instabilities.New J. Phys., 9(10):378, oct 2007
work page 2007
-
[56]
B. C. Ponedel and E. Knobloch. Forced snaking: Localized structures in the real Ginzburg-Landau equation with spatially periodic parametric forcing.Eur. Phys. J.: Spec. Top., 225(13):2549–2561, 2016
work page 2016
-
[57]
P. H. Rabinowitz. Some global results for nonlinear eigenvalue problems.J. Funct. Anal., 7(3):487–513, 1971. 25
work page 1971
- [58]
-
[59]
R. Richter and I. V. Barashenkov. Two-dimensional solitons on the surface of magnetic fluids.Phys. Rev. Lett., 94:184503, May 2005
work page 2005
- [60]
-
[61]
L. van der Aalst, J. B. van den Berg, and J.-P. Lessard. Periodic localised traveling waves in the two-dimensional suspension bridge equation.Nonlinearity, 38(7):Paper No. 075029, 2025
work page 2025
-
[62]
P. van Heijster and B. Sandstede. Planar radial spots in a three-component Fitzhugh–Nagumo system. J. Nonlinear Sci., 21(5):705–745, 2011
work page 2011
-
[63]
Verhulst.Nonlinear Differential Equations and Dynamical Systems
F. Verhulst.Nonlinear Differential Equations and Dynamical Systems. Universitext. Springer, 1996
work page 1996
-
[64]
M. H. Wheeler. Solitary water waves of large amplitude generated by surface pressure.Arch. Rat. Mech. Anal., 218(2):1131–1187, Nov 2015. 26
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.