Nonlinear Sciences

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 defect-like 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.

We present previously unknown solutions to the 3D Gross--Pitaevskii equation describing atomic Bose-Einstein condensates. This model supports elaborate patterns, including excited states bearing vorticity. The discovered coherent structures exhibit striking topological features, involving combinations of vortex rings and multiple, possibly bent vortex lines. Although unstable, many of them persist for long times in dynamical simulations. These solutions were identified by a state-of-the-art numerical technique called deflation, which is expected to be applicable to many problems from other areas of physics.

The time-dependent Ginzburg-Landau equation and the Swift-Hohenberg equation, both added with a stochastic term, are proposed to describe cloud pattern formation and cloud regime phase transitions of shallow convective clouds organized in mesoscale systems. The starting point is the Hottovy-Stechmann linear spatio-temporal stochastic model for tropical precipitation, used to describe the dynamics of water vapor and tropical convection. By taking into account that shallow stratiform clouds are close to a self-organized criticallity and that water vapor content is the order parameter, it is observed that sources must have non-linear terms in the equation to include the dynamical feedback due to precipitation and evaporation. The inclusion of this non-linearity leads to a kind of time-dependent Ginzburg-Landau stochastic equation, originally used to describe superconductivity phases. By performing numerical simulations, pattern formation is observed, specially for cellular convective phases. These patterns are much better compared with real satellite observations than the pure linear model. This is done by comparing the spatial Fourier transform of real and numerical cloud fields. Finally, by considering fluctuation terms for the turbulent eddy diffusion we arrive to a Hohenberg-Swift equation. The obtained patterns are much more organized that the patterns obtained from the Ginzburg-Landau equation in the case of closed cellular and roll convection.

Recent experimental results indicate that mixing is enhanced by a reciprocal flow induced inside a levitated droplet with an oscillatory deformation [T. Watanabe et al. Sci. Rep. 8, 10221 (2018)]. Generally, reciprocal flow cannot convect the solutes in time average, and agitation cannot take place. In the present paper, we focus on the diffusion process coupled with the reciprocal flow. We theoretically derive that the diffusion process can be enhanced by the reciprocal flow, and the results are confirmed via numerical calculation of the over-damped Langevin equation with a reciprocal flow.

Sedentism was a decisive moment in the history of humankind. In a review article Kay and Kaplan quantified land use for early human settlements and found that sedentism and the emergence of farming go hand in hand. For these settlements two primary land use categories, farming and living, can be identified, whereas for hunter gatherer societies no distinct differences can be made. It is natural to search for this in the behavior of two different groups, settlers and farmers. The development of two distinct zones and the two groups lead us to the hypothesis that the emergence of settlements is the result of diffusion-driven Turing instability. In this short communication we further specify this and show that this results in a regular settlement arrangement as can still be seen today in agricultural regions.

Spatial organization of proteins in cells is important for many biological functions. In general, the nonlinear, spatially coupled models for protein-pattern formation are only accessible to numerical simulations, which has limited insight into the general underlying principles. To overcome this limitation, we adopt the setting of two diffusively coupled, well-mixed compartments that represents the elementary feature of any pattern -- an interface. For intracellular systems, the total numbers of proteins are conserved on the relevant timescale of pattern formation. Thus, the essential dynamics is the redistribution of the globally conserved mass densities between the two compartments. We present a phase-portrait analysis in the phase-space of the redistributed masses that provides insights on the physical mechanisms underlying pattern formation. We demonstrate this approach for several paradigmatic model systems. In particular, we show that the pole-to-pole Min oscillations in Escherichia coli are relaxation oscillations of the MinD polarity orientation. This reveals a close relation between cell polarity oscillatory patterns in cells. Critically, our findings suggest that the design principles of intracellular pattern formation are found in characteristic features in these phase portraits (nullclines and fixed points). These features are not uniquely determined by the topology of the protein-interaction network but depend on parameters (kinetic rates, diffusion constants) and distinct networks can give rise to equivalent phase portrait features.

Scientists have observed and studied diffusive waves in contexts as disparate as population genetics and cell signaling. Often, these waves are propagated by discrete entities or agents, such as individual cells in the case of cell signaling. For a broad class of diffusive waves, we characterize the transition between the collective propagation of diffusive waves -- in which the wave speed is well-described by continuum theory -- and the propagation of diffusive waves by individual agents. We show that this transition depends heavily on the dimensionality of the system in which the wave propagates and that disordered systems yield dynamics largely consistent with lattice systems. In some system dimensionalities, the intuition that closely packed sources more accurately mimic a continuum can be grossly violated.

Whitham modulation theory describes the zero dispersion limit of nonlinear waves by a system of conservation laws for the parameters of modulated periodic traveling waves. Here, admissible, discontinuous, weak solutions of the Whitham modulation equations--termed Whitham shocks--are identified with zero dispersion limits of traveling wave solutions to higher order dispersive partial differential equations (PDEs). The far-field behavior of the traveling wave solutions satisfies the Rankine-Hugoniot jump conditions. Generally, the numerically computed traveling waves represent heteroclinic connections between two periodic orbits of an ordinary differential equation. The focus here is on the fifth order Korteweg-de Vries equation. The three admissible one-parameter families of Whitham shocks are used as solution components for the generalized Riemann problem of the Whitham modulation equations. Admissible KdV5-Whitham shocks are generally undercompressive, i.e., all characteristic families pass through the shock. The heteroclinic traveling waves that limit to admissible Whitham shocks are found to be ubiquitous in numerical simulations of smoothed step initial conditions for other higher order dispersive equations including the Kawahara equation (with third and fifth order dispersion) and a nonlocal model of weakly nonlinear gravity-capillary waves with full dispersion. Whitham shocks are linked to recent studies of nonlinear higher order dispersive waves in optics and ultracold atomic gases. The approach presented here provides a novel method for constructing new traveling wave solutions to dispersive nonlinear wave equations and a framework to identify physically relevant, admissible shock solutions of the Whitham modulation equations.

In this paper, we analyze the formation and dynamical properties of discrete light bullets (dLBs) in an array of passively mode-locked lasers coupled via evanescent fields in a ring geometry. Using a generic model based upon a system of nearest-neighbor coupled Haus master equations we show numerically the existence of dLBs for different coupling strengths. In order to reduce the complexity of the analysis, we approximate the full problem by a reduced set of discrete equations governing the dynamics of the transverse profile of the dLBs. This effective theory allows us to perform a detailed bifurcation analysis via path-continuation methods. In particular, we show the existence of multistable branches of discrete localized states (dLSs), corresponding to different number of active elements in the array. These branches are either independent of each other or are organized into a snaking bifurcation diagram where the width of the dLS grows via a process of successive increase and decrease of the gain. Mechanisms are revealed by which the snaking branches can be created and destroyed as a second parameter, e.g., the linewidth enhancement factor or the coupling strength are varied. For increasing couplings, the existence of moving bright and dark dLSs is also demonstrated.

Numerical simulation has indicated that vortex structures can exist for a long time in the form of quantized filaments on arrays of coupled weakly dissipative nonlinear oscillators in a finite three-dimensional domain under a resonant external force applied at the boundary of the domain. Ranges of the parameters of the system and an external signal favorable for the formation of modulationally stable quasi-uniform energy background, which is a decisive factor for the occurrence of this phenomenon, have been qualitatively revealed.