Emergence of rogue-like waves in a reaction-diffusion system: Stochastic output from deterministic dissipative dynamics.

Rogue waves are an intriguing nonlinear phenomenon arising across different scales, ranging from ocean waves through optics to Bose-Einstein condensates. We describe the emergence of rogue wave-like dynamics in a reaction-diffusion system that arise as a result of a subcritical Turing instability. This state is present in a regime where all time-independent states are unstable and consists of intermittent excitation of spatially localized spikes, followed by collapse to an unstable state and subsequent regrowth. We characterize the spatiotemporal organization of spikes and show that in sufficiently large domains the dynamics are consistent with a memoryless process.

Time-dependent localized patterns in a predator-prey model.

Numerical continuation is used to compute solution branches in a two-component reaction-diffusion model of Leslie-Gower type. Two regimes are studied in detail. In the first, the homogeneous state loses stability to supercritical spatially uniform oscillations, followed by a subcritical steady state bifurcation of Turing type. The latter leads to spatially localized states embedded in an oscillating background that bifurcate from snaking branches of localized steady states. Using two-parameter continuation, we uncover a novel mechanism whereby disconnected segments of oscillatory states zip up into a continuous snaking branch of time-periodic localized states, some of which are stable. In the second, the homogeneous state loses stability to supercritical Turing patterns, but steady spatially localized states embedded either in the homogeneous state or in a small amplitude Turing state are nevertheless present. We show that such behavior is possible when sideband Turing states are strongly subcritical and explain why this is so in the present model. In both cases, the observed behavior differs significantly from that expected on the basis of a supercritical primary bifurcation.

Suppression of Wall Modes in Rapidly Rotating Rayleigh-Bénard Convection by Narrow Horizontal Fins.

The heat transport by rapidly rotating Rayleigh-Bénard convection is of fundamental importance to many geophysical flows. Laboratory measurements are impeded by robust wall modes that develop along vertical walls, significantly perturbing the heat flux. We show that narrow horizontal fins along the vertical walls efficiently suppress wall modes ensuring that their contribution to the global heat flux is negligible compared with bulk convection in the geostrophic regime, thereby paving the way for new experimental studies of geophysically relevant regimes of rotating convection.

Front propagation and global bifurcations in a multivariable reaction-diffusion model.

We study the existence and stability of propagating fronts in Meinhardt's multivariable reaction-diffusion model of branching in one spatial dimension. We identify a saddle-node-infinite-period bifurcation of fronts that leads to episodic front propagation in the parameter region below propagation failure and show that this state is stable. Stable constant speed fronts exist only above this parameter value. We use numerical continuation to show that propagation failure is a consequence of the presence of a T-point corresponding to the formation of a heteroclinic cycle in a spatial dynamics description. Additional T-points are identified that are responsible for a large multiplicity of different unstable traveling front-peak states. The results indicate that multivariable models may support new types of behavior that are absent from typical two-variable models but may nevertheless be important in developmental processes such as branching and somitogenesis.

Universal Wrinkling of Supported Elastic Rings.

An exactly solvable family of models describing the wrinkling of substrate-supported inextensible elastic rings under compression is identified. The resulting wrinkle profiles are shown to be related to the buckled states of an unsupported ring and are therefore universal. Closed analytical expressions for the resulting universal shapes are provided, including the one-to-one relations between the pressure and tension at which these emerge. The analytical predictions agree with numerical continuation results to within numerical accuracy, for a large range of parameter values, up to the point of self-contact.

Instability mechanisms of repelling peak solutions in a multi-variable activator-inhibitor system.

We study the linear stability properties of spatially localized single- and multi-peak states generated in a subcritical Turing bifurcation in the Meinhardt model of branching. In one spatial dimension, these states are organized in a foliated snaking structure owing to peak-peak repulsion but are shown to be all linearly unstable, with the number of unstable modes increasing with the number of peaks present. Despite this, in two spatial dimensions, direct numerical simulations reveal the presence of stable single- and multi-spot states whose properties depend on the repulsion from nearby spots as well as the shape of the domain and the boundary conditions imposed thereon. Front propagation is shown to trigger the growth of new spots while destabilizing others. The results indicate that multi-variable models may support new types of behavior that are absent from typical two-variable models.

Spatially localized structures in the Gray-Scott model.

Spatially localized structures in the one-dimensional Gray-Scott reaction-diffusion model are studied using a combination of numerical continuation techniques and weakly nonlinear theory, focusing on the regime in which the activator and substrate diffusivities are different but comparable. Localized states arise in three different ways: in a subcritical Turing instability present in this regime, and from folds in the branch of spatially periodic Turing states. They also arise from the fold of spatially uniform states. These three solution branches interconnect in complex ways. We use numerical continuation techniques to explore their global behaviour within a formulation of the model that has been used to describe dryland vegetation patterns on a flat terrain.This article is part of the theme issue 'Dissipative structures in matter out of equilibrium: from chemistry, photonics and biology (part 2)'.

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.

Origin and stability of dark pulse Kerr combs in normal dispersion resonators.

We analyze dark pulse Kerr frequency combs in optical resonators with normal group-velocity dispersion using the Lugiato-Lefever model. We show that in the time domain the combs correspond to interlocked switching waves between the upper and lower homogeneous states, and explain how this fact accounts for many of their experimentally observed properties. Modulational instability does not play any role in their existence. We provide a detailed map indicating for which parameters stable dark pulse Kerr combs can be found, and how they are destabilized for increasing values of frequency detuning.

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.

Upscale energy transfer in three-dimensional rapidly rotating turbulent convection.

Rotating Rayleigh-Bénard convection exhibits, in the limit of rapid rotation, a turbulent state known as geostrophic turbulence. This state is present for sufficiently large Rayleigh numbers representing the thermal forcing of the system, and is characterized by a leading order balance between the Coriolis force and pressure gradient. This turbulent state is itself unstable to the generation of depth-independent or barotropic vortex structures of ever larger scale through a process known as spectral condensation. This process involves an inverse cascade mechanism with a positive feedback loop whereby large-scale barotropic vortices organize small scale convective eddies. In turn, these eddies provide a dynamically evolving energy source for the large-scale barotropic component. Kinetic energy spectra for the barotropic dynamics are consistent with a k-3 downscale enstrophy cascade and an upscale cascade that steepens to k-3 as the box-scale condensate forms. At the same time the flow maintains a baroclinic convective component with an inertial range consistent with a k-5/3 spectrum. The condensation process resembles a similar process in two dimensions but is fully three-dimensional.

Wake dynamics in buoyancy-driven flows: Steady-state-Hopf-mode interaction with O(2) symmetry revisited.

We present a detailed mathematical study of a truncated normal form relevant to the bifurcations observed in wake flow past axisymmetric bodies, with and without thermal stratification. We employ abstract normal form analysis to identify possible bifurcations and the corresponding bifurcation diagrams in parameter space. The bifurcations and the bifurcation diagrams are interpreted in terms of symmetry considerations. Particular emphasis is placed on the presence of attracting robust heteroclinic cycles in certain parameter regimes. The normal form coefficients are computed for several examples of wake flows behind buoyant disks and spheres, and the resulting predictions compared with the results of direct numerical flow simulations. In general, satisfactory agreement is obtained.

Bifurcation delay and front propagation in the real Ginzburg-Landau equation on a time-dependent domain.

This work analyzes bifurcation delay and front propagation in the one-dimensional real Ginzburg-Landau equation with periodic boundary conditions on isotropically growing or shrinking domains. First, we obtain closed-form expressions for the delay of primary bifurcations on a growing domain and show that the additional domain growth before the appearance of a pattern is independent of the growth time scale. We also quantify primary bifurcation delay on a shrinking domain; in contrast with a growing domain, the time scale of domain compression is reflected in the additional compression before the pattern decays. For secondary bifurcations such as the Eckhaus instability, we obtain a lower bound on the delay of phase slips due to a time-dependent domain. We also construct a heuristic model to classify regimes with arrested phase slips, i.e., phase slips that fail to develop. Then, we study how propagating fronts are influenced by a time-dependent domain. We identify three types of pulled fronts: homogeneous, pattern spreading, and Eckhaus fronts. By following the linear dynamics, we derive expressions for the velocity and profile of homogeneous fronts on a time-dependent domain. We also derive the natural "asymptotic" velocity and front profile and show that these deviate from predictions based on the marginal stability criterion familiar from fixed domain theory. This difference arises because the time dependence of the domain lifts the degeneracy of the spatial eigenvalues associated with speed selection and represents a fundamental distinction from the fixed domain theory that we verify using direct numerical simulations. The effect of a growing domain on pattern spreading and Eckhaus front velocities is inspected qualitatively and found to be similar to that of homogeneous fronts. These more complex fronts can also experience delayed onset. Lastly, we show that dilution-an effect present when the order parameter is conserved-increases bifurcation delay and amplifies changes in the homogeneous front velocity on time-dependent domains. The study provides general insight into the effects of domain growth on pattern onset, pattern transitions, and front propagation in systems across different scientific fields.

Elastic fingering in a rotating Hele-Shaw cell.

We consider the steady-state fingering instability of an elastic membrane separating two fluids of different density under external pressure in a rotating Hele-Shaw cell. Both inextensible and highly extensible membranes are considered, and the role of membrane tension is detailed in each case. Both systems exhibit a centrifugally driven Rayleigh-Taylor-like instability when the density of the inner fluid exceeds that of the outer one, and this instability competes with the restoring forces arising from curvature and tension, thereby setting the finger scale. Numerical continuation is used to compute not only strongly nonlinear primary finger states up to the point of self-contact, but also secondary branches of mixed modes and circumferentially localized folds as a function of the rotation rate and the externally imposed pressure. Both reflection-symmetric and symmetry-broken chiral states are computed. The results are presented in the form of bifurcation diagrams. The ratio of system scale to the natural length scale is found to determine the ordering of the primary bifurcations from the unperturbed circle state as well as the solution profiles and onset of secondary bifurcations.

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.

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.

Heat transport in low-Rossby-number Rayleigh-Bénard convection.

We demonstrate, via simulations of asymptotically reduced equations describing rotationally constrained Rayleigh-Bénard convection, that the efficiency of turbulent motion in the fluid bulk limits overall heat transport and determines the scaling of the nondimensional Nusselt number Nu with the Rayleigh number Ra, the Ekman number E, and the Prandtl number σ. For E << 1 inviscid scaling theory predicts and simulations confirm the large Ra scaling law Nu-1 ≈ C(1)σ(-1/2)Ra(3/2)E(2), where C(1) is a constant, estimated as C(1) ≈ 0.04 ± 0.0025. In contrast, the corresponding result for nonrotating convection, Nu-1 ≈ C(2)Ra(α), is determined by the efficiency of the thermal boundary layers (laminar: 0.28 ≤ α ≤ 0.31, turbulent: α ~ 0.38). The 3/2 scaling law breaks down at Rayleigh numbers at which the thermal boundary layer loses rotational constraint, i.e., when the local Rossby number ≈ 1. The breakdown takes place while the bulk Rossby number is still small and results in a gradual transition to the nonrotating scaling law. For low Ekman numbers the location of this transition is independent of the mechanical boundary conditions.

Electrolyte stability in a nanochannel with charge regulation.

The stability of an electrolyte confined in one dimension between two solid surfaces is analyzed theoretically in the case where overlapping double layers produce nontrivial interactions. Within the Poisson-Boltzmann-Nernst-Planck description of the electrostatic interaction and transport of electrical charges, the presence of Stern layers can enrich the set of possible solutions. Our analytical and numerical study of the stability properties of the trivial state of this system identified an instability to a new antisymmetric state. This state is stable for a range of gap widths that depends on the Debye and Stern lengths, but for smaller gap widths, where the Stern layers overlap, a second transition takes place and the stable nontrivial solution diverges. The origin of this divergence is explained and its properties analyzed using asymptotic techniques which are in good agreement with numerical results. The relevance of our results to confined electrolytes at nanometer scales is discussed in the context of energy storage in nanometric systems.

Localized and extended patterns in the cubic-quintic Swift-Hohenberg equation on a disk.

Axisymmetric and nonaxisymmetric patterns in the cubic-quintic Swift-Hohenberg equation posed on a disk with Neumann boundary conditions are studied via numerical continuation and bifurcation analysis. Axisymmetric localized solutions in the form of spots and rings known from earlier studies persist and snake in the usual fashion until they begin to interact with the boundary. Depending on parameters, including the disk radius, these states may or may not connect to the branch of domain-filling target states. Secondary instabilities of localized axisymmetric states may create multiarm localized structures that grow and interact with the boundary before broadening into domain-filling states. High azimuthal wave number wall states referred to as daisy states are also found. Secondary bifurcations from these states include localized daisies, i.e., wall states localized in both radius and angle. Depending on parameters, these states may snake much as in the one-dimensional Swift-Hohenberg equation, or invade the interior of the domain, yielding states referred to as worms, or domain-filling stripes.

Origin of jumping oscillons in an excitable reaction-diffusion system.

Oscillons, i.e., immobile spatially localized but temporally oscillating structures, are the subject of intense study since their discovery in Faraday wave experiments. However, oscillons can also disappear and reappear at a shifted spatial location, becoming jumping oscillons (JOs). We explain here the origin of this behavior in a three-variable reaction-diffusion system via numerical continuation and bifurcation theory, and show that JOs are created via a modulational instability of excitable traveling pulses (TPs). We also reveal the presence of bound states of JOs and TPs and patches of such states (including jumping periodic patterns) and determine their stability. This rich multiplicity of spatiotemporal states lends itself to information and storage handling.