Asymptotics of large bound states of localized structures

Speed of wave packets and the nonlinear Schrödinger equation

The universal theory of weakly nonlinear wave packets given by the nonlinear Schrödinger equation is revisited. In the limit where the group and phase velocities are very close together, a multiple-scale analysis carried out beyond all orders reveals that a single soliton, bright or dark, can travel at a different speed than the group velocity. In an exponentially small but finite range of parameters, the envelope of the soliton is locked to the rapid oscillations of the carrier wave. Eventually, the dynamics is governed by an equation analogous to that of a pendulum, in which the center of mass of the soliton is subjected to a periodic potential. Consequently, the soliton speed is not constant and generally contains a periodic component. Furthermore, the interaction between two distant solitons can in principle be profoundly altered by the aforementioned effective periodic potential and we conjecture the existence of new bound states. These results are derived on a wide class of wave models and in such a general way that they are believed to be of universal validity.

Dissipative switching waves and solitons in the systems with spontaneously broken symmetry

The paper addresses the bistability caused by spontaneous symmetry breaking bifurcation in a one-dimensional periodically corrugated nonlinear waveguide pumped by coherent light at normal incidence. The formation and the stability of the switching waves connecting the states of different symmetries are studied numerically. It is shown that the switching waves can form stable resting and moving bound states (dissipative solitons). The protocols of the creation of the discussed nonlinear localized waves are suggested and verified by numerical simulations.

Formation of localized states in dryland vegetation: Bifurcation structure and stability

We study theoretically the emergence of localized states of vegetation close to the onset of desertification. These states are formed through the locking of vegetation fronts, connecting a uniform vegetation state with a bare soil state, which occurs nearby the Maxwell point of the system. To study these structures we consider a universal model of vegetation dynamics in drylands, which has been obtained as the normal form for different vegetation models. Close to the Maxwell point localized gaps and spots of vegetation exist and undergo collapsed snaking. The presence of gaps strongly suggest that the ecosystem may undergo a recovering process. In contrast, the presence of spots may indicate that the ecosystem is close to desertification.

Continuum Modeling of Discrete Plant Communities: Why Does It Work and Why Is It Advantageous?

Understanding ecosystem response to drier climates calls for modeling the dynamics of dryland plant populations, which are crucial determinants of ecosystem function, as they constitute the basal level of whole food webs. Two modeling approaches are widely used in population dynamics, individual (agent)-based models and continuum partial-differential-equation (PDE) models. The latter are advantageous in lending themselves to powerful methodologies of mathematical analysis, but the question of whether they are suitable to describe small discrete plant populations, as is often found in dryland ecosystems, has remained largely unaddressed. In this paper, we first draw attention to two aspects of plants that distinguish them from most other organisms—high phenotypic plasticity and dispersal of stress-tolerant seeds—and argue in favor of PDE modeling, where the state variables that describe population sizes are not discrete number densities, but rather continuous biomass densities. We then discuss a few examples that demonstrate the utility of PDE models in providing deep insights into landscape-scale behaviors, such as the onset of pattern forming instabilities, multiplicity of stable ecosystem states, regular and irregular, and the possible roles of front instabilities in reversing desertification. We briefly mention a few additional examples, and conclude by outlining the nature of the information we should and should not expect to gain from PDE model studies.

Localized structures in dispersive and doubly resonant optical parametric oscillators

We study temporally localized structures in doubly resonant degenerate optical parametric oscillators in the absence of temporal walk-off. We focus on states formed through the locking of domain walls between the zero and a nonzero continuous-wave solution. We show that these states undergo collapsed snaking and we characterize their dynamics in the parameter space.

Fluid-supported elastic sheet under compression: Multifold solutions

The properties of a hinged floating elastic sheet of finite length under compression are considered. Numerical continuation is used to compute spatially localized buckled states with many spatially localized folds. Both symmetric and antisymmetric states are computed and the corresponding bifurcation diagrams determined. Weakly nonlinear analysis is used to analyze the transition from periodic wrinkles to singlefold and multifold states and to compute their energy. States with the same number of folds have energies that barely differ from each other and the energy gap decreases exponentially as localization increases. The stability of the different competing states is studied and the multifold solutions are all found to be unstable. However, the decay time into solutions with fewer folds can be so slow that multifolds may appear to be stable.

Snakes and ghosts in a parity-time-symmetric chain of dimers

We consider linearly coupled discrete nonlinear Schrödinger equations with gain and loss terms and with a cubic-quintic nonlinearity. The system models a parity-time (PT)-symmetric coupler composed by a chain of dimers. We study uniform states and site-centered and bond-centered spatially localized solutions and present that each solution has a symmetric and antisymmetric configuration between the arms. The symmetric solutions can become unstable due to bifurcations of asymmetric ones, that are called ghost states, because they exist only when an otherwise real propagation constant is taken to be complex valued. When a parameter is varied, the resulting bifurcation diagrams for the existence of standing localized solutions have a snaking behavior. The critical gain and loss coefficient above which the PT symmetry is broken corresponds to the condition when bifurcation diagrams of symmetric and antisymmetric states merge. Past the symmetry breaking, the system no longer has time-independent states. Nevertheless, equilibrium solutions can be analytically continued by defining a dual equation that leads to ghost states associated with growth or decay, that are also identified and examined here. We show that ghost localized states also exhibit snaking bifurcation diagrams. We analyze the width of the snaking region and provide asymptotic approximations in the limit of strong and weak coupling where good agreement is obtained.

Bifurcation structure of localized states in the Lugiato-Lefever equation with anomalous dispersion

The origin, stability, and bifurcation structure of different types of bright localized structures described by the Lugiato-Lefever equation are studied. This mean field model describes the nonlinear dynamics of light circulating in fiber cavities and microresonators. In the case of anomalous group velocity dispersion and low values of the intracavity phase detuning these bright states are organized in a homoclinic snaking bifurcation structure. We describe how this bifurcation structure is destroyed when the detuning is increased across a critical value, and determine how a bifurcation structure known as foliated snaking emerges.

Homoclinic snaking in the discrete Swift-Hohenberg equation

We consider the discrete Swift-Hohenberg equation with cubic and quintic nonlinearity, obtained from discretizing the spatial derivatives of the Swift-Hohenberg equation using central finite differences. We investigate the discretization effect on the bifurcation behavior, where we identify three regions of the coupling parameter, i.e., strong, weak, and intermediate coupling. Within the regions, the discrete Swift-Hohenberg equation behaves either similarly or differently from the continuum limit. In the intermediate coupling region, multiple Maxwell points can occur for the periodic solutions and may cause irregular snaking and isolas. Numerical continuation is used to obtain and analyze localized and periodic solutions for each case. Theoretical analysis for the snaking and stability of the corresponding solutions is provided in the weak coupling region.

Robust chimera states in SQUID metamaterials with local interactions

We report on the emergence of robust multiclustered chimera states in a dissipative-driven system of symmetrically and locally coupled identical superconducting quantum interference device (SQUID) oscillators. The "snakelike" resonance curve of the single SQUID is the key to the formation of the chimera states and is responsible for the extreme multistability exhibited by the coupled system that leads to attractor crowding at the geometrical resonance (inductive-capacitive) frequency. Until now, chimera states were mostly believed to exist for nonlocal coupling. Our findings provide theoretical evidence that nearest-neighbor interactions are indeed capable of supporting such states in a wide parameter range. SQUID metamaterials are the subject of intense experimental investigations, and we are highly confident that the complex dynamics demonstrated in this paper can be confirmed in the laboratory.

Motion of dissipative optical fronts under the action of an oscillating pump

The dynamics of domain walls in optical bistable systems with pump and loss is considered. It is shown that an oscillating component of the pump affects the average drift velocity of the domain walls. The cases of harmonic and biharmonic pumps are considered. It is demonstrated that in the case of biharmonic pulse the velocity of the domain wall can be controlled by the mutual phase of the harmonics. The analogy between this phenomenon and the ratchet effect is drawn. Synchronization of the moving domain walls by the oscillating pump in discrete systems is studied and discussed.

Subcritical Turing bifurcation and the morphogenesis of localized patterns

Subcritical Turing bifurcations of reaction-diffusion systems in large domains lead to spontaneous onset of well-developed localized patterns via the homoclinic snaking mechanism. This phenomenon is shown to occur naturally when balancing source and loss effects are included in a typical reaction-diffusion system, leading to a super- to subcritical transition. Implications are discussed [corrected]for a range of physical problems, arguing that subcriticality leads to naturally robust phase transitions to localized patterns.

Analytical results for front pinning between an hexagonal pattern and a uniform state in pattern-formation systems

In pattern-forming systems, localized patterns are states of intermediate complexity between fully extended ordered patterns and completely irregular patterns. They are formed by stationary fronts enclosing an ordered pattern inside an homogeneous background. In two dimensions, the ordered pattern is most often hexagonal and the conditions for fronts to stabilize are still unknown. In this Letter, we show how the locking of these fronts depends on their orientation relative to the pattern. The theory rests on general asymptotic arguments valid when the spatial scale of the front is slow compared to that of the hexagonal pattern. Our analytical results are confirmed by numerical simulations with the Swift-Hohenberg equation, relevant to hydrodynamical and buckling instabilities, and a nonlinear optical cavity model.

Localized oscillatory states in magnetoconvection

Localized states are found in many pattern forming systems. The aim of this paper is to investigate the occurrence of oscillatory localized states in two-dimensional Boussinesq magnetoconvection. Initially considering an idealized model, in which the vertical structure of the system has been simplified by a projection onto a small number of Fourier modes, we find that these states are restricted to the low ζ regime (where ζ represents the ratio of the magnetic to thermal diffusivities). These states always exhibit bistability with another nontrivial solution branch; in other words, they show no evidence of subcritical behavior. This is due to the weak flux expulsion that is exhibited by these time-dependent solutions. Using the results of this parameter survey, we locate corresponding states in a fully resolved two-dimensional system, although the mode of oscillation is more complex in this case. This is the first time that a localized oscillatory state, of this kind, has been found in a fully resolved magnetoconvection simulation.

Homoclinic snaking of localized patterns in a spatially forced system

Dissipative localized structures exhibit intricate bifurcation diagrams. An adequate theory has been developed in one space dimension; however, discrepancies arise with the experiments. Based on an optical feedback with spatially modulated input beam, we set up a 1D forced configuration in a nematic liquid crystal layer. We characterize experimentally and theoretically the homoclinic snaking diagram of localized patterns, providing a reconciliation between theory and experiments.

Variational approximations to homoclinic snaking in continuous and discrete systems

Localized structures appear in a wide variety of systems, arising from a pinning mechanism due to the presence of a small-scale pattern or an imposed grid. When there is a separation of length scales, the width of the pinning region is exponentially small and beyond the reach of standard asymptotic methods. We show how this behavior can be obtained using a variational method, for two systems. In the case of the quadratic-cubic Swift-Hohenberg equation, this gives results that are in agreement with recent work using exponential asymptotics. In addition, the method is applied to a discrete system with cubic-quintic nonlinearity, giving results that agree well with numerical simulations.

Variational approximations to homoclinic snaking

We investigate the snaking of localized patterns, seen in numerous physical applications, using a variational approximation. This method naturally introduces the exponentially small terms responsible for the snaking structure, which are not accessible via standard multiple-scales asymptotic techniques. We obtain the symmetric snaking solutions and the asymmetric "ladder" states, and also predict the stability of the localized states. The resulting approximate formulas for the width of the snaking region show good agreement with numerical results.

Dual-mode spiral vortices

We show that spiral vortices in oscillatory systems can lose stability to secondary modes to form dual-mode spiral vortices. The secondary modes grow at the vortex core where the oscillation amplitude vanishes but are nonlinearly damped by the oscillatory mode away from the core. Gradients of the oscillation phase, induced by the hosted secondary mode, can lead to additional hosting events that culminate in periodic core oscillations or in a novel form of spatiotemporal chaos. The results of this study apply to physical, chemical, and biological systems that go through cusp-Hopf, fold-Hopf, and Hopf-Turing bifurcations.

Influence of boundaries on localized patterns

We analytically study the influence of boundaries on distant localized patterns generated by a Turing instability. To this end, we use the Swift-Hohenberg model with arbitrary boundary conditions. We find that the bifurcation diagram of these localized structures generally involves four homoclinic snaking branches, rather than two for infinite or periodic domains. Second, steady localized patterns only exist at discrete locations, and only at the center of the domain if their size exceeds a critical value. Third, reducing the domain size increases the pinning range.

Homoclinic snaking in bounded domains

Homoclinic snaking is a term used to describe the back and forth oscillation of a branch of time-independent spatially localized states in a bistable spatially reversible system as the localized structure grows in length by repeatedly adding rolls on either side. On the real line this process continues forever. In finite domains snaking terminates once the domain is filled but the details of how this occurs depend critically on the choice of boundary conditions. With periodic boundary conditions the snaking branches terminate on a branch of spatially periodic states. However, with non-Neumann boundary conditions they turn continuously into a large amplitude filling state that replaces the periodic state. This behavior, shown here in detail for the Swift-Hohenberg equation, explains the phenomenon of "snaking without bistability," recently observed in simulations of binary fluid convection by Mercader et al. Phys. Rev. E 80, 025201 (2009).

Localized pinning states in closed containers: Homoclinic snaking without bistability

Binary mixtures with a negative separation ratio are known to exhibit time-independent spatially localized convection when heated from below. Numerical continuation of such states in a closed two-dimensional container with experimental boundary conditions and parameter values reveals the presence of a pinning region in Rayleigh number with multiple stable localized states but no bistability between the conduction state and an independent container-filling state. An explanation for this unusual behavior is offered.

Localized beating between dynamically generated frequencies

We analyze the beating between intrinsic frequencies that are simultaneously generated by a modulation (Turing) instability in a nonlinear extended system. The model studied is that of a coherently driven photonic crystal fiber cavity. Beating in the form of a slow modulation of fast intensity oscillations is found to be stable for a wide range of parameters. We find that such beating can also be localized and contain only a finite number of slow modulations. These structures consist of dips in the amplitude of the fast intensity oscillations, which can either be isolated or regularly spaced. An asymptotic analysis close to the modulation instability threshold allows us to explain this phenomenon as a manifestation of homoclinic snaking for dissipative localized structures.

Homoclinic snaking in a semiconductor-based optical system

We report on experimental observations of homoclinic snaking in a vertical-cavity semiconductor optical amplifier. Our observations in a quasi-one-dimensional and two-dimensional configurations agree qualitatively well with what is expected from recent theoretical and numerical studies. In particular, we show the bifurcation sequence leading to a snaking bifurcation diagram linking single localized states to "localized patterns" or clusters of localized states and demonstrate a parameter region where cluster states are inhibited.

Eckhaus instability and homoclinic snaking

Homoclinic snaking is a term used to describe the back and forth oscillation of a branch of time-independent spatially localized states in a bistable, spatially reversible system as the localized structure grows in length by repeatedly adding rolls on either side. This behavior is simplest to understand within the subcritical Swift-Hohenberg equation, but is also present in the subcritical regime of doubly diffusive convection driven by horizontal gradients. In systems that are unbounded in one spatial direction homoclinic snaking continues indefinitely as the localized structure grows to resemble a spatially periodic state of infinite extent. In finite domains or in periodic domains with finite spatial period the process must terminate. In this paper we show that the snaking branches in general turn over once the length of the localized state becomes comparable to the domain, and examine the factors that determine the location of the termination point or points, and their relation to the Eckhaus instability of the spatially periodic state.

Secondary Turing-type instabilities due to strong spatial resonance

We investigate the dynamics of pattern-forming systems in large domains near a codimension-two point corresponding to a ‘strong spatial resonance’ where competing instabilities with wavenumbers in the ratio 1 : 2 or 1 : 3 occur. We supplement the standard amplitude equations for such a mode interaction with Ginzburg–Landau-type modulational terms, appropriate to pattern formation in a large domain. In cases where the coefficients of these new diffusive terms differ substantially from each other, we show that spatially periodic solutions found near onset may be unstable to two long-wavelength modulational instabilities. Moreover, these instabilities generically occur near the codimension-two point only in the 1 : 2 and 1 : 3 cases, and not when weaker spatial resonances arise. The first instability is ‘amplitude-driven’ and is the analogue of the well-known Turing instability of reaction–diffusion systems. The second is a phase instability for which the subsequent nonlinear development is described, at leading order, by the Cahn–Hilliard equation.

The normal forms for strong spatial resonances are also well known to permit uniformly travelling wave solutions. We also show that these travelling waves may be similarly unstable.

Proposed resolution of theory-experiment discrepancy in homoclinic snaking

In spatially extended Turing-unstable systems, parameter variation should, in theory, produce only fully developed patterns. In experiment, however, localized patterns or solitons sitting on a smooth background often appear. Addition of a nonlocal nonlinearity can resolve this discrepancy by tilting the "snaking" bifurcation diagram characteristic of such problems.

Homoclinic snaking: structure and stability

The bistable Swift-Hohenberg equation exhibits multiple stable and unstable spatially localized states of arbitrary length in the vicinity of the Maxwell point between spatially homogeneous and periodic states. These states are organized in a characteristic snakes-and-ladders structure. The origin of this structure in one spatial dimension is reviewed, and the stability properties of the resulting states with respect to perturbations in both one and two dimensions are described. The relevance of the results to several different physical systems is discussed.

Nonvariational real Swift-Hohenberg equation for biological, chemical, and optical systems

We derive asymptotically an order parameter equation in the limit where nascent bistability and long-wavelength modulation instabilities coalesce. This equation is a variant of the Swift-Hohenberg equation that generally contains nonvariational terms of the form psinabla(2)psi and /nablapsi/(2). We briefly review some of the properties already derived for this equation and derive it on three examples taken from chemical, biological, and optical contexts. Finally, we derive the equation on a general class of partial differential systems.

On homoclinic snaking in optical systems

The existence of localized structures, including so-called cavity solitons, in driven optical systems is discussed. In theory, they should exist only below the threshold of a subcritical modulational instability, but in experiment they often appear spontaneously on parameter variation. The addition of a nonlocal nonlinearity may resolve this discrepancy by tilting the "snaking" bifurcation diagram characteristic of such problems.

Reciprocal oscillons and nonmonotonic fronts in forced nonequilibrium systems

The formation of oscillons in a synchronously oscillating background is studied in the context of both damped and self-exciting oscillatory media. Using the forced complex Ginzburg-Landau equation we show that such states bifurcate from finite amplitude homogenous states near the 2:1 resonance boundary. In each case we identify a region in parameter space containing a finite multiplicity of coexisting stable oscillons with different structure. Stable time-periodic monotonic and nonmonotonic frontlike states are present in an overlapping region. Both types of structure are related to the presence of a Maxwell point between the zero and finite amplitude homogeneous states.