Homoclinic snaking: structure and stability

Square waves and Bykov T-points in a delay algebraic model for the Kerr-Gires-Tournois interferometer

We study theoretically the mechanisms of square wave formation of a vertically emitting micro-cavity operated in the Gires-Tournois regime that contains a Kerr medium and that is subjected to strong time-delayed optical feedback and detuned optical injection. We show that in the limit of large delay, square wave solutions of the time-delayed system can be treated as relative homoclinic solutions of an equation with an advanced argument. Based on this, we use concepts of classical homoclinic bifurcation theory to study different types of square wave solutions. In particular, we unveil the mechanisms behind the collapsed snaking scenario of square waves and explain the formation of complex-shaped multistable square wave solutions through a Bykov T-point. Finally, we relate the position of the T-point to the position of the Maxwell point in the original time-delayed system.

VisualPDE: Rapid Interactive Simulations of Partial Differential Equations

Computing has revolutionised the study of complex nonlinear systems, both by allowing us to solve previously intractable models and through the ability to visualise solutions in different ways. Using ubiquitous computing infrastructure, we provide a means to go one step further in using computers to understand complex models through instantaneous and interactive exploration. This ubiquitous infrastructure has enormous potential in education, outreach and research. Here, we present VisualPDE, an online, interactive solver for a broad class of 1D and 2D partial differential equation (PDE) systems. Abstract dynamical systems concepts such as symmetry-breaking instabilities, subcritical bifurcations and the role of initial data in multistable nonlinear models become much more intuitive when you can play with these models yourself, and immediately answer questions about how the system responds to changes in parameters, initial conditions, boundary conditions or even spatiotemporal forcing. Importantly, VisualPDE is freely available, open source and highly customisable. We give several examples in teaching, research and knowledge exchange, providing high-level discussions of how it may be employed in different settings. This includes designing web-based course materials structured around interactive simulations, or easily crafting specific simulations that can be shared with students or collaborators via a simple URL. We envisage VisualPDE becoming an invaluable resource for teaching and research in mathematical biology and beyond. We also hope that it inspires other efforts to make mathematics more interactive and accessible.

Stationary broken parity states in active matter models

We demonstrate that several nonvariational continuum models commonly used to describe active matter as well as other active systems exhibit nongeneric behavior: each model supports asymmetric but stationary localized states even in the absence of pinning at heterogeneities. Moreover, such states only begin to drift following a drift-transcritical bifurcation as the activity increases. Asymmetric stationary states should only exist in variational systems, i.e., in models with gradient structure. In other words, such states are expected in passive systems, but not in active systems where the gradient structure of the model is broken by activity. We identify a "spurious" gradient dynamics structure of these models that is responsible for this nongeneric behavior, and determine the types of additional terms that render the models generic, i.e., with asymmetric states that appear via drift-pitchfork bifurcations and are generically moving. We provide detailed illustrations of our results using numerical continuation of resting and steadily drifting states in both generic and nongeneric cases.

Collisions of localized patterns in a nonvariational Swift-Hohenberg equation

The cubic-quintic Swift-Hohenberg equation (SH35) has been proposed as an order parameter description of several convective systems with reflection symmetry in the layer midplane, including binary fluid convection. We use numerical continuation, together with extensive direct numerical simulations (DNSs), to study SH35 with an additional nonvariational quadratic term to model the effects of breaking the midplane reflection symmetry. The nonvariational structure of the model leads to the propagation of asymmetric spatially localized structures (LSs). An asymptotic prediction for the drift velocity of such structures, derived in the limit of weak symmetry breaking, is validated numerically. Next, we present an extensive study of possible collision scenarios between identical and nonidentical traveling structures, varying a temperaturelike control parameter. These collisions are inelastic and result in stationary or traveling structures. Depending on system parameters and the types of structures colliding, the final state may be a simple bound state of the initial LSs, but it can also be longer or shorter than the sum of the two initial states as a result of nonlinear interactions. The Maxwell point of the variational system, where the free energy of the global pattern state equals that of the trivial state, is shown to have no bearing on which of these scenarios is realized. Instead, we argue that the stability properties of bound states are key. While individual LSs lie on a modified snakes-and-ladders structure in the nonvariational SH35, the multipulse bound states resulting from collisions lie on isolas in parameter space, disconnected from the trivial solution. In the gradient SH35, such isolas are always of figure-eight shape, but in the present nongradient case they are generically more complex, although the figure-eight shape is preserved in a small subset of cases. Some of these complex isolas are shown to terminate in T-point bifurcations. A reduced model is proposed to describe the interactions between the tails of the LSs. The model consists of two coupled ordinary differential equations (ODEs) capturing the oscillatory nature of SH35 profiles at the linear level. It contains three parameters: two interaction amplitudes and a phase, whose values are deduced from high-resolution DNSs using gradient descent optimization. For collisions leading to the formation of simple bound states, the reduced model reproduces the trajectories of LSs with high quantitative accuracy. When nonlinear interactions lead to the creation or deletion of wavelengths, the model performs less well. Finally, we propose an effective signature of a given interaction in terms of net attraction or repulsion relative to free propagation. It is found that interactions can be attractive or repulsive in the net, irrespective of whether the two closest interacting extrema are of the same or opposite signs. Our findings highlight the rich temporal dynamics described by this bistable nonvariational SH35, and show that the interactions in this system can be quantitatively captured, to a significant extent, by a highly reduced ODE model.

Probing the stability landscape of cylindrical shells for buckling knockdown factors

The buckling response of axially compressed cylindrical shells is well known for its imperfection sensitivity. Mapping out a stability landscape by localized probing has recently been proposed as a rational means for establishing shell buckling knockdown factors. Probing using a point force directed radially inwards and perpendicular to the cylinder wall is based on the insight that a localized single dimple exists as an edge state in the basin boundary of the stable prebuckling equilibrium. Here, we extend the idea of probing to bi-directional inwards and outwards forces to trigger both single-dimple and double-dimple edge states. We identify key features of the ensuing probing stability landscape and generalize these to derive three design curves of varying conservatism that are a function of the non-dimensional Batdorf parameter only. Interestingly, the most conservative of the three knockdown curves bounds a large dataset of experimental buckling results from below, despite being derived from probing features of geometrically perfect cylinders. Overall, the three design curves permit a more nuanced design approach than legacy knockdown factors, as different levels of conservatism can be chosen based on expected manufacturing quality. For instance, the most and least conservative of the three design guidelines differ by a factor of 3 for the most slender cylinder geometries, and the associated reduction in safety factor has profound implications for efficient structural design. This article is part of the theme issue 'Probing and dynamics of shock sensitive shells'.

Transitions between dissipative localized structures in the simplified Gilad-Meron model for dryland plant ecology

Spatially extended patterns and multistability of possible different states are common in many ecosystems, and their combination has an important impact on their dynamical behaviors. One potential combination involves tristability between a patterned state and two different uniform states. Using a simplified version of the Gilad-Meron model for dryland ecosystems, we study the organization, in bifurcation terms, of the localized structures arising in tristable regimes. These states are generally related to the concept of wave front locking and appear in the form of spots and gaps of vegetation. We find that the coexistence of localized spots and gaps, within tristable configurations, yields the appearance of hybrid states. We also study the emergence of spatiotemporal localized states consisting of a portion of a periodic pattern embedded in a uniform Hopf-like oscillatory background in a subcritical Turing-Hopf dynamical regime.

Temporal dissipative solitons in the Morris-Lecar model with time-delayed feedback

We study the dynamics and bifurcations of temporal dissipative solitons in an excitable system under time-delayed feedback. As a prototypical model displaying different types of excitability, we use the Morris-Lecar model. In the limit of large delay, soliton like solutions of delay-differential equations can be treated as homoclinic solutions of an equation with an advanced argument. Based on this, we use concepts of classical homoclinic bifurcation theory to study different types of pulse solutions and to explain their dependence on the system parameters. In particular, we show how a homoclinic orbit flip of a single-pulse soliton leads to the destabilization of equidistant multi-pulse solutions and to the emergence of stable pulse packages. It turns out that this transition is induced by a heteroclinic orbit flip in the system without feedback, which is related to the excitability properties of the Morris-Lecar model.

Square-wave generation in vertical external-cavity Kerr-Gires-Tournois interferometers

We study theoretically the mechanisms of square-wave (SW) formation in vertical external-cavity Kerr-Gires-Tournois interferometers in the presence of anti-resonant injection. We provide simple analytical approximations for their plateau intensities and for the conditions of their emergence. We demonstrate that SWs may appear via a homoclinic snaking scenario, leading to the formation of complex-shaped multistable SW solutions. The resulting SWs can host localized structures and robust bound states.

Bifurcation structure of traveling pulses in type-I excitable media

We study the scenario in which traveling pulses emerge in a prototypical type-I one-dimensional excitable medium, which exhibits two different routes to excitable behavior, mediated by a homoclinic (saddle-loop) and a saddle-node on the invariant cycle bifurcations. We characterize the region in parameter space in which traveling pulses are stable together with the different bifurcations behind either their destruction or loss of stability. In particular, some of the bifurcations delimiting the stability region have been connected, using singular limits, with the two different scenarios that mediated type-I local excitability. Finally, the existence of traveling pulses has been linked to a drift pitchfork instability of localized steady structures.

Ferrofluids and bio-ferrofluids: looking back and stepping forward

Ferrofluids investigated along for about five decades are ultrastable colloidal suspensions of magnetic nanoparticles, which manifest simultaneously fluid and magnetic properties. Their magnetically controllable and tunable feature proved to be from the beginning an extremely fertile ground for a wide range of engineering applications. More recently, biocompatible ferrofluids attracted huge interest and produced a considerable increase of the applicative potential in nanomedicine, biotechnology and environmental protection. This paper offers a brief overview of the most relevant early results and a comprehensive description of recent achievements in ferrofluid synthesis, advanced characterization, as well as the governing equations of ferrohydrodynamics, the most important interfacial phenomena and the flow properties. Finally, it provides an overview of recent advances in tunable and adaptive multifunctional materials derived from ferrofluids and a detailed presentation of the recent progress of applications in the field of sensors and actuators, ferrofluid-driven assembly and manipulation, droplet technology, including droplet generation and control, mechanical actuation, liquid computing and robotics.

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.

Two-dimensional localized states in an active phase-field-crystal model

The active phase-field-crystal (active PFC) model provides a simple microscopic mean field description of crystallization in active systems. It combines the PFC model (or conserved Swift-Hohenberg equation) of colloidal crystallization and aspects of the Toner-Tu theory for self-propelled particles. We employ the active PFC model to study the occurrence of localized and periodic active crystals in two spatial dimensions. Due to the activity, crystalline states can undergo a drift instability and start to travel while keeping their spatial structure. Based on linear stability analyses, time simulations, and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of resting and traveling states. We explore, for instance, how the slanted homoclinic snaking of steady localized states found for the passive PFC model is modified by activity. Morphological phase diagrams showing the regions of existence of various solution types are presented merging the results from all the analysis tools employed. We also study how activity influences the crystal structure with transitions from hexagons to rhombic and stripe patterns. This in-depth analysis of a simple PFC model for active crystals and swarm formation provides a clear general understanding of the observed multistability and associated hysteresis effects, and identifies thresholds for qualitative changes in behavior.

Dissipative Light Bullets in Kerr Cavities: Multistability, Clustering, and Rogue Waves

We report the existence of stable dissipative light bullets in Kerr cavities. These three-dimensional (3D) localized structures consist of either an isolated light bullet (LB), bound together, or could occur in clusters forming well-defined 3D patterns. They can be seen as stationary states in the reference frame moving with the group velocity of light within the cavity. The number of LBs and their distribution in 3D settings are determined by the initial conditions, while their maximum peak power remains constant for a fixed value of the system parameters. Their bifurcation diagram allows us to explain this phenomenon as a manifestation of homoclinic snaking for dissipative light bullets. However, when the strength of the injected beam is increased, LBs lose their stability and the cavity field exhibits giant, short-living 3D pulses. The statistical characterization of pulse amplitude reveals a long tail probability distribution, indicating the occurrence of extreme events, often called rogue waves.

Efficient calculation of phase coexistence and phase diagrams: application to a binary phase-field-crystal model

We show that one can employ well-established numerical continuation methods to efficiently calculate the phase diagram for thermodynamic systems described by a suitable free energy functional. In particular, this involves the determination of lines of phase coexistence related to first order phase transitions and the continuation of triple points. To illustrate the method we apply it to a binary phase-field-crystal model for the crystallisation of a mixture of two types of particles. The resulting phase diagram is determined for one- and two-dimensional domains. In the former case it is compared to the diagram obtained from a one-mode approximation. The various observed liquid and crystalline phases and their stable and metastable coexistence are discussed as well as the temperature-dependence of the phase diagrams. This includes the (dis)appearance of critical points and triple points. We also relate bifurcation diagrams for finite-size systems to the thermodynamics of phase transitions in the infinite-size limit.

Reversible Self-Replication of Spatiotemporal Kerr Cavity Patterns

We uncover a novel and robust phenomenon that causes the gradual self-replication of spatiotemporal Kerr cavity patterns in cylindrical microresonators. These patterns are inherently synchronized multifrequency combs. Under proper conditions, the axially localized nature of the patterns leads to a fundamental drift instability that induces transitions among patterns with a different number of rows. Self-replications, thus, result in the stepwise addition or removal of individual combs along the cylinder's axis. Transitions occur in a fully reversible and, consequently, deterministic way. The phenomenon puts forward a novel paradigm for Kerr frequency comb formation and reveals important insights into the physics of multidimensional nonlinear patterns.

Bending and pinching of three-phase stripes: From secondary instabilities to morphological deformations in organic photovoltaics

Optimizing the properties of the mosaic nanoscale morphology of bulk heterojunction (BHJ) organic photovoltaics (OPV) is not only challenging technologically but also intriguing from the mechanistic point of view. Among the recent breakthroughs is the identification and utilization of a three-phase (donor-mixed-acceptor) BHJ, where the (intermediate) mixed phase can inhibit mesoscale morphological changes, such as phase separation. Using a mean-field approach, we reveal and distinguish between generic mechanisms that alter, through transverse instabilities, the evolution of stripes: the bending (zigzag mode) and the pinching (cross-roll mode) of the donor-acceptor domains. The results are summarized in a parameter plane spanned by the mixing energy and illumination, and show that donor-acceptor mixtures with higher mixing energy are more likely to develop pinching under charge-flux boundary conditions. The latter is notorious as it leads to the formation of disconnected domains and hence to loss of charge flux. We believe that these results provide a qualitative road map for BHJ optimization, using mixed-phase composition and, therefore, an essential step toward long-lasting OPV. More broadly, the results are also of relevance to study the coexistence of multiple-phase domains in material science, such as in ion-intercalated rechargeable batteries.

Phase-field-crystal description of active crystallites: Elastic and inelastic collisions

The active Phase-Field-Crystal (aPFC) model combines elements of the Toner-Tu theory for self-propelled particles and the classical Phase-Field-Crystal (PFC) model that describes the transition between liquid and crystalline phases. In the liquid-crystal coexistence region of the PFC model, crystalline clusters exist in the form of localized states that coexist with a homogeneous background. At sufficiently strong activity (related to self-propulsion strength), they start to travel. We employ numerical path continuation and direct time simulations to first investigate the existence regions of different types of localized states in one spatial dimension. The results are summarized in morphological phase diagrams in the parameter plane spanned by activity and mean density. Then we focus on the interaction of traveling localized states, studying their collision behavior. As a result, we distinguish "elastic" and "inelastic" collisions. In the former, localized states recover their properties after a collision, while in the latter, they may completely or partially annihilate, forming resting bound states or various traveling states.

Stripes on finite domains: Why the zigzag instability is only a partial story

Stationary periodic patterns are widespread in natural sciences, ranging from nano-scale electrochemical and amphiphilic systems to mesoscale fluid, chemical, and biological media and to macro-scale vegetation and cloud patterns. Their formation is usually due to a primary symmetry breaking of a uniform state to stripes, often followed by secondary instabilities to form zigzag and labyrinthine patterns. These secondary instabilities are well studied under idealized conditions of an infinite domain; however, on finite domains, the situation is more subtle since the unstable modes depend also on boundary conditions. Using two prototypical models, the Swift-Hohenberg equation and the forced complex Ginzburg-Landau equation, we consider finite size domains with no flux boundary conditions transversal to the stripes and reveal a distinct mixed-mode instability that lies in between the classical zigzag and the Eckhaus lines. This explains the stability of stripes in the mildly zigzag unstable regime and, after crossing the mixed-mode line, the evolution of zigzag stripes in the bulk of the domain and the formation of defects near the boundaries. The results are of particular importance for problems with large timescale separation, such as bulk-heterojunction deformations in organic photovoltaic and vegetation in semi-arid regions, where early temporal transients may play an important role.

Chaotic motion of localized structures

Mobility properties of spatially localized structures arising from chaotic but deterministic forcing of the bistable Swift-Hohenberg equation are studied and compared with the corresponding results when the chaotic forcing is replaced by white noise. Short structures are shown to possess greater mobility, resulting in larger root-mean-square speeds but shorter displacements than longer structures. Averaged over realizations, the displacement of the structure is ballistic at short times but diffusive at larger times. Similar results hold in two spatial dimensions. The effects of chaotic forcing on the stability of these structures is also quantified. Shorter structures are found to be more fragile than longer ones, and their stability region can be displaced outside the pinning region for constant forcing. Outside the stability region the deterministic fluctuations lead either to the destruction of the structure or to its gradual growth.

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.

Origin of localized snakes-and-ladders solutions of plane Couette flow

Spatially localized invariant solutions of plane Couette flow are organized in a snakes-and-ladders structure strikingly similar to that observed for simpler pattern-forming partial differential equations [Schneider, Gibson, and Burke, Phys. Rev. Lett. 104, 104501 (2010)PRLTAO0031-900710.1103/PhysRevLett.104.104501]. We demonstrate the mechanism by which these snaking solutions originate from well-known periodic states of the Taylor-Couette system. They are formed by a localized slug of wavy-vortex flow that emerges from a background of Taylor vortices via a modulational sideband instability. This mechanism suggests a close connection between pattern-formation theory and Navier-Stokes flow.

Tunable Kerr frequency combs and temporal localized states in time-delayed Gires-Tournois interferometers

In this Letter, we study theoretically a new setup allowing for the generation of temporal localized states (TLSs) and frequency combs. The setup is compact (a few centimeters) and can be implemented using established technologies, while offering tunable repetition rates and potentially high power operation. It consists of a vertically emitting micro-cavity, operated in the Gires-Tournois regime, containing a Kerr medium strong time-delayed optical feedback, and detuned optical injection. We disclose sets of multistable dark and bright TLSs coexisting on their respective bistable homogeneous backgrounds.

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.

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.

Localized Faraday patterns under heterogeneous parametric excitation

Faraday waves are a classic example of a system in which an extended pattern emerges under spatially uniform forcing. Motivated by systems in which uniform excitation is not plausible, we study both experimentally and theoretically the effect of heterogeneous forcing on Faraday waves. Our experiments show that vibrations restricted to finite regions lead to the formation of localized subharmonic wave patterns and change the onset of the instability. The prototype model used for the theoretical calculations is the parametrically driven and damped nonlinear Schrödinger equation, which is known to describe well Faraday-instability regimes. For an energy injection with a Gaussian spatial profile, we show that the evolution of the envelope of the wave pattern can be reduced to a Weber-equation eigenvalue problem. Our theoretical results provide very good predictions of our experimental observations provided that the decay length scale of the Gaussian profile is much larger than the pattern wavelength.

Chimera states in a Duffing oscillators chain coupled to nearest neighbors

Coupled nonlinear oscillators can present complex spatiotemporal behaviors. Here, we report the coexistence of coherent and incoherent domains, called chimera states, in an array of identical Duffing oscillators coupled to their nearest neighbors. The chimera states show a significant variation of amplitude in the desynchronized domain. These intriguing states are observed in the bistability region between a homogeneous state and a spatiotemporal chaotic one. These dynamical behaviors are characterized by their Lyapunov spectra and their global phase coherence order parameter. The local coupling between oscillators prevents one domain from invading the other one. Depending on initial conditions, a family of chimera states appear, organized in a snaking-like diagram.

Resting and traveling localized states in an active phase-field-crystal model

The conserved Swift-Hohenberg equation (or phase-field-crystal [PFC] model) provides a simple microscopic description of the thermodynamic transition between fluid and crystalline states. Combining it with elements of the Toner-Tu theory for self-propelled particles, Menzel and Löwen [Phys. Rev. Lett. 110, 055702 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.055702] obtained a model for crystallization (swarm formation) in active systems. Here, we study the occurrence of resting and traveling localized states, i.e., crystalline clusters, within the resulting active PFC model. Based on linear stability analyses and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of periodic and localized, resting and traveling states in a one-dimensional active PFC model. This allows us, for instance, to explore how the slanted homoclinic snaking of steady localized states found for the passive PFC model is amended by activity. A particular focus lies on the onset of motion, where we show that it occurs either through a drift-pitchfork or a drift-transcritical bifurcation. A corresponding general analytical criterion is derived.

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.

Tristability between stripes, up-hexagons, and down-hexagons and snaking bifurcation branches of spatial connections between up- and down-hexagons

Third-order amplitude equations on hexagonal lattices can be used for predicting the existence and stability of stripes, up-hexagons, and down-hexagons in pattern-forming systems. These amplitude equations predict the nonexistence of bistable ranges between up- and down-hexagons and tristable ranges between stripes, up-, and down-hexagons. In the present work we use fifth-order amplitude equations for finding such bistable and tristable ranges for a generalized Swift-Hohenberg equation and discuss stationary front connections between up- and down-hexagons.

Spontaneous motion of localized structures induced by parity symmetry breaking transition

We consider a paradigmatic nonvariational scalar Swift-Hohenberg equation that describes short wavenumber or large wavelength pattern forming systems. This work unveils evidence of the transition from stable stationary to moving localized structures in one spatial dimension as a result of a parity breaking instability. This behavior is attributed to the nonvariational character of the model. We show that the nature of this transition is supercritical. We characterize analytically and numerically this bifurcation scenario from which emerges asymmetric moving localized structures. A generalization for two-dimensional settings is discussed.

Implications of tristability in pattern-forming ecosystems

Many ecosystems show both self-organized spatial patterns and multistability of possible states. The combination of these two phenomena in different forms has a significant impact on the behavior of ecosystems in changing environments. One notable case is connected to tristability of two distinct uniform states together with patterned states, which has recently been found in model studies of dryland ecosystems. Using a simple model, we determine the extent of tristability in parameter space, explore its effects on the system dynamics, and consider its implications for state transitions or regime shifts. We analyze the bifurcation structure of model solutions that describe uniform states, periodic patterns, and hybrid states between the former two. We map out the parameter space where these states exist, and note how the different states interact with each other. We further focus on two special implications with ecological significance, breakdown of the snaking range and complex fronts. We find that the organization of the hybrid states within a homoclinic snaking structure breaks down as it meets a Maxwell point where simple fronts are stationary. We also discover a new series of complex fronts between the uniform states, each with its own velocity. We conclude with a brief discussion of the significance of these findings for the dynamics of regime shifts and their potential control.

Nonlinear mode interaction in equal-leg angle struts susceptible to cellular buckling

A variational model that describes the interactive buckling of a thin-walled equal-leg angle strut under pure axial compression is presented. A formulation combining the Rayleigh-Ritz method and continuous displacement functions is used to derive a system of differential and integral equilibrium equations for the structural component. Solving the equations using numerical continuation reveals progressive cellular buckling (or snaking) arising from the nonlinear interaction between the weak-axis flexural buckling mode and the strong-axis flexural-torsional buckling mode for the first time-the resulting behaviour being highly unstable. Physical experiments conducted on 10 cold-formed steel specimens are presented and the results show good agreement with the variational model.

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.

Effects of nonlinear gradient terms on the defect turbulence regime in weakly dissipative systems

We investigate the behavior of traveling waves in a defect turbulence regime with the periodic boundary conditions by using the lowest-order complex Ginzburg-Landau equation (CGLE), and we show the effect of the nonlinear gradient terms in the system. It is found that the nonlinear gradient terms which appear at the same order as the quintic term can change the behavior of the wave patterns. The presence of the nonlinear gradient terms can cause major changes in the behavior of the solution. They can be considered like the stabilizing terms. The system which was initially unstable or chaotic can become stable by including the nonlinear gradient terms.

Two-dimensional localized chaotic patterns in parametrically driven systems

We study two-dimensional localized patterns in weakly dissipative systems that are driven parametrically. As a generic model for many different physical situations we use a generalized nonlinear Schrödinger equation that contains parametric forcing, damping, and spatial coupling. The latter allows for the existence of localized pattern states, where a finite-amplitude uniform state coexists with an inhomogeneous one. In particular, we study numerically two-dimensional patterns. Increasing the driving forces, first the localized pattern dynamics is regular, becomes chaotic for stronger driving, and finally extends in area to cover almost the whole system. In parallel, the spatial structure of the localized states becomes more and more irregular, ending up as a full spatiotemporal chaotic structure.

Streamwise localization of traveling wave solutions in channel flow

Channel flow of an incompressible fluid at Reynolds numbers above 2400 possesses a number of different spatially localized solutions that approach laminar flow far upstream and downstream. We use one such relative time-periodic solution, which corresponds to a spatially localized version of a Tollmien-Schlichting wave, to illustrate how the upstream and downstream asymptotics can be computed analytically. In particular, we show that for these spanwise uniform states the asymptotics predict the exponential localization that has been observed for numerically computed solutions of several canonical shear flows but never properly understood theoretically.

Delay-induced depinning of localized structures in a spatially inhomogeneous Swift-Hohenberg model

We report on the dynamics of localized structures in an inhomogeneous Swift-Hohenberg model describing pattern formation in the transverse plane of an optical cavity. This real order parameter equation is valid close to the second-order critical point associated with bistability. The optical cavity is illuminated by an inhomogeneous spatial Gaussian pumping beam and subjected to time-delayed feedback. The Gaussian injection beam breaks the translational symmetry of the system by exerting an attracting force on the localized structure. We show that the localized structure can be pinned to the center of the inhomogeneity, suppressing the delay-induced drift bifurcation that has been reported in the particular case where the injection is homogeneous, assuming a continuous wave operation. Under an inhomogeneous spatial pumping beam, we perform the stability analysis of localized solutions to identify different instability regimes induced by time-delayed feedback. In particular, we predict the formation of two-arm spirals, as well as oscillating and depinning dynamics caused by the interplay of an attracting inhomogeneity and destabilizing time-delayed feedback. The transition from oscillating to depinning solutions is investigated by means of numerical continuation techniques. Analytically, we use an order parameter approach to derive a normal form of the delay-induced Hopf bifurcation leading to an oscillating solution. Additionally we model the interplay of an attracting inhomogeneity and destabilizing time delay by describing the localized solution as an overdamped particle in a potential well generated by the inhomogeneity. In this case, the time-delayed feedback acts as a driving force. Comparing results from the later approach with the full Swift-Hohenberg model, we show that the approach not only provides an instructive description of the depinning dynamics, but also is numerically accurate throughout most of the parameter regime.

Chimera-type states induced by local coupling

Coupled oscillators can exhibit complex self-organization behavior such as phase turbulence, spatiotemporal intermittency, and chimera states. The latter corresponds to a coexistence of coherent and incoherent states apparently promoted by nonlocal or global coupling. Here we investigate the existence, stability properties, and bifurcation diagram of chimera-type states in a system with local coupling without different time scales. Based on a model of a chain of nonlinear oscillators coupled to adjacent neighbors, we identify the required attributes to observe these states: local coupling and bistability between a stationary and an oscillatory state close to a homoclinic bifurcation. The local coupling prevents the incoherent state from invading the coherent one, allowing concurrently the existence of a family of chimera states, which are organized by a homoclinic snaking bifurcation diagram.

Localized auxin peaks in concentration-based transport models of the shoot apical meristem

We study the formation of auxin peaks in a generic class of concentration-based auxin transport models, posed on static plant tissues. Using standard asymptotic analysis, we prove that, on bounded domains, auxin peaks are not formed via a Turing instability in the active transport parameter, but via simple corrections to the homogeneous steady state. When the active transport is small, the geometry of the tissue encodes the peaks' amplitude and location: peaks arise where cells have fewer neighbours, that is, at the boundary of the domain. We test our theory and perform numerical bifurcation analysis on two models that are known to generate auxin patterns for biologically plausible parameter values. In the same parameter regimes, we find that realistic tissues are capable of generating a multitude of stationary patterns, with a variable number of auxin peaks, that can be selected by different initial conditions or by quasi-static changes in the active transport parameter. The competition between active transport and production rate determines whether peaks remain localized or cover the entire domain. In particular, changes in the auxin production that are fast with respect to the cellular life cycle affect the auxin peak distribution, switching from localized spots to fully patterned states. We relate the occurrence of localized patterns to a snaking bifurcation structure, which is known to arise in a wide variety of nonlinear media, but has not yet been reported in plant models.

Non-unique results of collisions of quasi-one-dimensional dissipative solitons

We investigate collisions of quasi-one-dimensional dissipative solitons (DSs) for a large class of initial conditions, which are not temporally asymptotic quasi-one-dimensional DSs. For the case of sufficiently small approach velocity and sufficiently large values of the dissipative cross-coupling between the counter-propagating DSs, we find non-unique results for the outcome of collisions. We demonstrate that these non-unique results are intrinsically related to a modulation instability along the crest of the quasi-one-dimensional objects. As a model, we use coupled cubic-quintic complex Ginzburg-Landau equations. Among the final results found are stationary and oscillatory compound states as well as more complex assemblies consisting of quasi-one-dimensional and localized states. We analyse to what extent the final results can be described by the solutions of one cubic-quintic complex Ginzburg-Landau equation with effective parameters.

Recent advances in symmetric and network dynamics

We summarize some of the main results discovered over the past three decades concerning symmetric dynamical systems and networks of dynamical systems, with a focus on pattern formation. In both of these contexts, extra constraints on the dynamical system are imposed, and the generic phenomena can change. The main areas discussed are time-periodic states, mode interactions, and non-compact symmetry groups such as the Euclidean group. We consider both dynamics and bifurcations. We summarize applications of these ideas to pattern formation in a variety of physical and biological systems, and explain how the methods were motivated by transferring to new contexts René Thom's general viewpoint, one version of which became known as "catastrophe theory." We emphasize the role of symmetry-breaking in the creation of patterns. Topics include equivariant Hopf bifurcation, which gives conditions for a periodic state to bifurcate from an equilibrium, and the H/K theorem, which classifies the pairs of setwise and pointwise symmetries of periodic states in equivariant dynamics. We discuss mode interactions, which organize multiple bifurcations into a single degenerate bifurcation, and systems with non-compact symmetry groups, where new technical issues arise. We transfer many of the ideas to the context of networks of coupled dynamical systems, and interpret synchrony and phase relations in network dynamics as a type of pattern, in which space is discretized into finitely many nodes, while time remains continuous. We also describe a variety of applications including animal locomotion, Couette-Taylor flow, flames, the Belousov-Zhabotinskii reaction, binocular rivalry, and a nonlinear filter based on anomalous growth rates for the amplitude of periodic oscillations in a feed-forward network.

Localization in a spanwise-extended model of plane Couette flow

We consider a nine-partial-differential-equation (1-space and 1-time) model of plane Couette flow in which the degrees of freedom are severely restricted in the streamwise and cross-stream directions to study spanwise localization in detail. Of the many steady Eckhaus (spanwise modulational) instabilities identified of global steady states, none lead to a localized state. Spatially localized, time-periodic solutions were found instead, which arise in saddle node bifurcations in the Reynolds number. These solutions appear global (domain filling) in narrow (small spanwise) domains yet can be smoothly continued out to fully spanwise-localized states in very wide domains. This smooth localization behavior, which has also been seen in fully resolved duct flow (S. Okino, Ph.D. thesis, Kyoto University, Kyoto, 2011), indicates that an apparently global flow structure does not have to suffer a modulational instability to localize in wide domains.

Localized states in periodically forced systems

The theory of stationary spatially localized patterns in dissipative systems driven by time-independent forcing is well developed. With time-periodic forcing, related but time-dependent structures may result. These may consist of breathing localized patterns, or states that grow for part of the cycle via nucleation of new wavelengths of the pattern followed by wavelength annihilation during the remainder of the cycle. These two competing processes lead to a complex phase diagram whose structure is a consequence of a series of resonances between the nucleation time and the forcing period. The resulting diagram is computed for the periodically forced quadratic-cubic Swift-Hohenberg equation, and its details are interpreted in terms of the properties of the depinning transition for the fronts bounding the localized state on either side. The results are expected to shed light on localized states in a large variety of periodically driven systems.

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.

Observations of highly localized oscillons with multiple crests and troughs

Two types of stable, highly localized Faraday resonant standing waves with multiple crests and troughs are observed in an ethanol-water solution partly filled in a Hele-Shaw cell vertically oscillated with a single frequency. Systematical experiments are performed to investigate the properties of these oscillons. It is found that the wave height of these oscillons is independent of fluid depth from 1 to 5 cm. In particular, some experiments are performed to indicate the high localization of the oscillons, which suggests that these oscillons may be regarded as a combination of the two elementary oscillons discovered by Rajchenbach et al. [Phys. Rev. Lett. 107, 024502 (2011)], for instance, (2,3)=(1,1)+(1,2), where (m,n) denotes an oscillon with m crests and n troughs. So, our experiments also reveal an elegant "arithmetic" of these oscillons. These experimental phenomena are helpful to deepen and enrich our understanding about Faraday waves.

Cellular buckling in stiffened plates

An analytical model based on variational principles for a thin-walled stiffened plate subjected to axial compression is presented. A system of nonlinear differential and integral equations is derived and solved using numerical continuation. The results show that the system is susceptible to highly unstable local–global mode interaction after an initial instability is triggered. Moreover, snap-backs in the response showing sequential destabilization and restabilization, known as cellular buckling or snaking, arise. The analytical model is compared with static finite element (FE) models for joint conditions between the stiffener and the main plate that have significant rotational restraint. However, it is known from previous studies that the behaviour, where the same joint is insignificantly restrained rotationally, is captured better by an analytical approach than by standard FE methods; the latter being unable to capture cellular buckling behaviour even though the phenomenon is clearly observed in laboratory experiments.

Class of compound dissipative solitons as a result of collisions in one and two spatial dimensions

We study the interaction of quasi-one-dimensional (quasi-1D) dissipative solitons (DSs). Starting with quasi-1D solutions of the cubic-quintic complex Ginzburg-Landau (CGL) equation in their temporally asymptotic state as the initial condition, we find, as a function of the approach velocity and the real part of the cubic interaction of the two counterpropagating envelopes: interpenetration, one compound state made of both envelopes or two compound states. For the latter class both envelopes show DSs superposed at two different locations. The stability of this class of compound states is traced back to the quasilinear growth rate associated with the coupled system. We show that this mechanism also works for 1D coupled cubic-quintic CGL equations. For quasi-1D states that are not in their asymptotic state before the collision, a breakup along the crest can be observed, leading to nonunique results after the collision of quasi-1D states.