Spatial localization in heterogeneous systems

Bumps and oscillons in networks of spiking neurons

We study localized patterns in an exact mean-field description of a spatially extended network of quadratic integrate-and-fire neurons. We investigate conditions for the existence and stability of localized solutions, so-called bumps, and give an analytic estimate for the parameter range, where these solutions exist in parameter space, when one or more microscopic network parameters are varied. We develop Galerkin methods for the model equations, which enable numerical bifurcation analysis of stationary and time-periodic spatially extended solutions. We study the emergence of patterns composed of multiple bumps, which are arranged in a snake-and-ladder bifurcation structure if a homogeneous or heterogeneous synaptic kernel is suitably chosen. Furthermore, we examine time-periodic, spatially localized solutions (oscillons) in the presence of external forcing, and in autonomous, recurrently coupled excitatory and inhibitory networks. In both cases, we observe period-doubling cascades leading to chaotic oscillations.

Defectlike structures and localized patterns in the cubic-quintic-septic Swift-Hohenberg equation

We study numerically the cubic-quintic-septic Swift-Hohenberg (SH357) equation on bounded one-dimensional domains. Under appropriate conditions stripes with wave number k≈1 bifurcate supercritically from the zero state and form S-shaped branches resulting in bistability between small and large amplitude stripes. Within this bistability range we find stationary heteroclinic connections or fronts between small and large amplitude stripes, and demonstrate that the associated spatially localized defectlike structures either snake or fall on isolas. In other parameter regimes we also find heteroclinic connections to spatially homogeneous states and a multitude of dynamically stable steady states consisting of patches of small and large amplitude stripes with different wave numbers or of spatially homogeneous patches. The SH357 equation is thus extremely rich in the types of patterns it exhibits. Some of the features of the bifurcation diagrams obtained by numerical continuation can be understood using a conserved quantity, the spatial Hamiltonian of the system.

Connection between nonlinear energy optimization and instantons

How systems transit between different stable states under external perturbation is an important practical issue. We discuss here how a recently developed energy optimization method for identifying the minimal disturbance necessary to reach the basin boundary of a stable state is connected to the instanton trajectory from large deviation theory of noisy systems. In the context of the one-dimensional Swift-Hohenberg equation, which has multiple stable equilibria, we first show how the energy optimization method can be straightforwardly used to identify minimal disturbances-minimal seeds-for transition to specific attractors from the ground state. Then, after generalizing the technique to consider multiple, equally spaced-in-time perturbations, it is shown that the instanton trajectory is indeed the solution of the energy optimization method in the limit of infinitely many perturbations provided a specific norm is used to measure the set of discrete perturbations. Importantly, we find that the key features of the instanton can be captured by a low number of discrete perturbations (typically one perturbation per basin of attraction crossed). This suggests a promising new diagnostic for systems for which it may be impractical to calculate the instanton.

Non-monotonic resonance in a spatially forced Lengyel-Epstein model

We study resonant spatially periodic solutions of the Lengyel-Epstein model modified to describe the chlorine dioxide-iodine-malonic acid reaction under spatially periodic illumination. Using multiple-scale analysis and numerical simulations, we obtain the stability ranges of 2:1 resonant solutions, i.e., solutions with wavenumbers that are exactly half of the forcing wavenumber. We show that the width of resonant wavenumber response is a non-monotonic function of the forcing strength, and diminishes to zero at sufficiently strong forcing. We further show that strong forcing may result in a π/2 phase shift of the resonant solutions, and argue that the nonequilibrium Ising-Bloch front bifurcation can be reversed. We attribute these behaviors to an inherent property of forcing by periodic illumination, namely, the increase of the mean spatial illumination as the forcing amplitude is increased.

Spatial forcing of pattern-forming systems that lack inversion symmetry

The entrainment of periodic patterns to spatially periodic parametric forcing is studied. Using a weak nonlinear analysis of a simple pattern formation model we study the resonant responses of one-dimensional systems that lack inversion symmetry. Focusing on the first three n:1 resonances, in which the system adjusts its wavenumber to one nth of the forcing wavenumber, we delineate commonalities and differences among the resonances. Surprisingly, we find that all resonances show multiplicity of stable phase states, including the 1:1 resonance. The phase states in the 2:1 and 3:1 resonances, however, differ from those in the 1:1 resonance in remaining symmetric even when the inversion symmetry is broken. This is because of the existence of a discrete translation symmetry in the forced system. As a consequence, the 2:1 and 3:1 resonances show stationary phase fronts and patterns, whereas phase fronts within the 1:1 resonance are propagating and phase patterns are transients. In addition, we find substantial differences between the 2:1 resonance and the other two resonances. While the pattern forming instability in the 2:1 resonance is supercritical, in the 1:1 and 3:1 resonances it is subcritical, and while the inversion asymmetry extends the ranges of resonant solutions in the 1:1 and 3:1 resonances, it has no effect on the 2:1 resonance range. We conclude by discussing a few open questions.