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\"odinger 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.

Complex systems may morph between structures with different dimensionality and degrees of freedom. As a tool for their modelling, nonlinear embeddings are introduced that encompass objects with different dimensionality as a continuous parameter $\kappa \in \mathbb{R}$ is being varied, thus allowing the unification of vectors, matrices and tensors in single mathematical structures. This technique is applied to construct warped models in the passage from supergravity in 10 or 11-dimensional spacetimes to 4-dimensional ones. We also show how nonlinear embeddings can be used to connect cellular automata (CAs) to coupled map lattices (CMLs) and to nonlinear partial differential equations, deriving a class of nonlinear diffusion equations. Finally, by means of nonlinear embeddings we introduce CA connections, a class of CMLs that connect any two arbitrary CAs in the limits $\kappa \to 0$ and $\kappa \to \infty$ of the embedding.

We introduce "superheated integrability," which produces characteristic staircase transmission plots for barrier collisions of breathers of the nonlinear Schr\"odinger equation. The effect makes tangible the inverse scattering transform, which treats the velocities and norms of the constituent solitons as the real and imaginary parts of the eigenvalues of the Lax operator. If all the norms are much greater than the velocities, an integrability-breaking potential may nonperturbatively change the velocities while having no measurable effect on the norms. This could be used to improve atomic interferometers.

We study the power spectrum of a space-time dependent neural field which describes the average membrane potential of neurons in a single layer. This neural field is modelled by a dissipative integro-differential equation, the so-called Amari equation. By considering a small perturbation with respect to a stationary and uniform configuration of the neural field we derive a linearized equation which is solved for a generic external stimulus by using the Fourier transform into wavevector-freqency domain, finding an analytical formula for the power spectrum of the neural field. In addition, after proving that for large wavelengths the linearized Amari equation is equivalent to a diffusion equation which admits space-time dependent analytical solutions, we take into account the nonlinearity of the Amari equation. We find that for large wavelengths a weak nonlinearity in the Amari equation gives rise to a reaction-diffusion equation which can be formally derived from a neural action functional by introducing a dual neural field. For some initial conditions, we discuss analytical solutions of this reaction-diffusion equation.

This paper proposes a control theoretic framework to model and analyze the self-organized pattern formation of molecular concentrations in biomolecular communication networks, emerging applications in synthetic biology. In biomolecular communication networks, bionanomachines, or biological cells, communicate with each other using a cell-to-cell communication mechanism mediated by a diffusible signaling molecule, thereby the dynamics of molecular concentrations are approximately modeled as a reaction-diffusion system with a single diffuser. We first introduce a feedback model representation of the reaction-diffusion system and provide a systematic local stability/instability analysis tool using the root locus of the feedback system. The instability analysis then allows us to analytically derive the conditions for the self-organized spatial pattern formation, or Turing pattern formation, of the bionanomachines. We propose a novel synthetic biocircuit motif called activator-repressor-diffuser system and show that it is one of the minimum biomolecular circuits that admit self-organized patterns over cell population.

New exact solution class of Born -- Infeld type nonlinear scalar field model is obtained. The variational principle of this model has a specific form which is characteristic for extremal four-dimensional hypersurface or hyper-film in five-dimensional space-time. Obtained solutions are singular solitons propagating with speed of light and having energy, momentum, and angular momentum which can be calculated for explicit conditions. Such solitons will be called the lightlike ones. The soliton singularity has a form of moving two-dimensional surface or shell. The lightlike soliton can have a set of tubelike singular shells with the appropriate cavities. A twisted lightlike soliton is considered. It is notable that its energy is proportional to its angular momentum in high-frequency approximation. A case with one tubelike cavity is considered. In this case the soliton shell is diffeomorphic to a cylindrical surface with threads by multifilar helix. The shell transverse size of the appropriate finite energy soliton can be converging to zero at infinity. The ideal gas of such lightlike solitons with minimal twist parameter is considered in a finite volume. Explicit conditions provide that the angular momentum of each soliton in the volume equals Planck constant. The equilibrium energy spectral density for the solitons is obtained. It has the form of Planck distribution in some approximation. A beam of the twisted lightlike solitons is considered. The representation of arbitrary polarization for the beam with the twisted lightlike solitons is discussed. It is shown that the effect of mechanical angular momentum transfer to absorbent by the circularly polarized beam can be provided. This effect is well known for photon beam. Thus the soliton solution which have determinate likeness with photon is obtained in particular.

We study solitary wave propagation in 1D granular crystals with Hertz-like interaction potentials. We consider interfaces between media with different exponents in the interaction potential. For an interface with increasing interaction potential exponent along the propagation direction we obtain mainly transmission with delayed secondary transmitted and reflected pulses. For interfaces with decreasing interaction potential exponent we observe both significant reflection and transmission of the solitary wave, where the transmitted part of the wave forms a multipulse structure. We also investigate impurities consisting of beads with different interaction exponents compared to the media they are embedded in, and we find that the impurities cause both reflection and transmission, including the formation of multipulse structures, independent of whether the exponent in the impurities is smaller than in the surrounding media. We explain wave propagation effects at interfaces and impurities in terms of quasi-particle collisions. Next we consider wave propagation along Hertz-like granular chains of beads in the presence of disorder and periodicity in the interaction exponents present in the Hertz-like potential, modelling, for instance, inhomogeneity in the contact geometry between beads in the granular chain. We find that solitary waves in media with randomised interaction exponents (which models disorder in the contact geometry) experience exponential decay, where the dependence of the decay rate is similar to the case of randomised bead masses. In the periodic case of chains with interaction exponents alternating between two fixed values, we find qualitatively different propagation properties depending on the choice of the two exponents. In particular, we find regimes with either exponential decay or stable solitary wave propagation with pairwise collective behaviour.

The stationary Gross-Pitaevskii equation in one dimension is considered with a complex periodic potential satisfying the conditions of the PT (parity-time reversal) symmetry. Under rather general assumptions on the potentials we prove bifurcations of PT-symmetric nonlinear bound states from the end points of a real interval in the spectrum of the non-selfadjoint linear Schrodinger operator with a complex PT-symmetric periodic potential. The nonlinear bound states are approximated by the effective amplitude equation, which bears the form of the cubic nonlinear Schrodinger equation. In addition we provide sufficient conditions for the appearance of complex spectral bands when the complex $\PT$-symmetric potential has an asymptotically small imaginary part.

Beginning with addition and multiplication which are intrinsic to a Koch-type curve, I formulate and solve a wave equation that describes wave propagation along a fractal coastline. As opposed to the examples known from the literature I do not replace the fractal by the continuum in which it is embedded. This seems to be the first example of a truly intrinsic description of wave propagation along a fractal curve.

We report measurements for the initial stages of sidebranching during the dendritic growth of ammonium chloride from supersaturated aqueous solution. The earliest sidebranches are approximately periodic; they are first evident about $36 \rho$ behind the tip, where $\rho$ is the tip radius, and have an average initial spacing of about $5 \rho$, though both values show considerable variation. The initial sidebranch amplitude grows approximately exponentially, but quickly saturates as sidebranches compete and coarsening sets in. This initial sidebranch growth is reasonably consistent with what would be expected for noise-driven sidebranches, though there are some quantitative differences.

In this article we provide a systematic way of creating generalized Moran sets using an analogous iterated function system (IFS) procedure. We use a step-wise adjustable IFS to introduce some variance (such as non-self-similarity) in the fractal limit sets. The process retains the computational simplicity of a standard IFS procedure. In our construction of the generalized Moran sets, we also weaken the fourth Moran Structure Condition that requires the same pattern of diameter ratios be used across a generation. Moreover, we provide upper and lower bounds for the Hausdorff dimension of the fractals created from this generalized process. Specific examples (Cantor-like sets, Sierpinski-like Triangles, etc) with the calculations of their corresponding dimensions are studied.

By studying the coarsening dynamics of a one-dimensional spin-1 Bose-Hubbard model in a superfluid regime, we analytically find an unconventional universal dynamical scaling for the growth of the spin correlation length, which is characterized by the exponential integral unlike the conventional power-law or simple logarithmic behavior, and numerically confirmed with the truncated Wigner approximation.

Zonal jets in a barotropic setup emerge out of homogeneous turbulence through a flow-forming instability of the homogeneous turbulent state (`zonostrophic instability') which occurs as the turbulence intensity increases. This has been demonstrated using the statistical state dynamics (SSD) framework with a closure at second order. Furthermore, it was shown that for small supercriticality the flow-forming instability follows Ginzburg-Landau (G-L) dynamics. Here, the SSD framework is used to study the equilibration of this flow-forming instability for small supercriticality. First, we compare the predictions of the weakly nonlinear G-L dynamics to the fully nonlinear SSD dynamics closed at second order for a wide ranges of parameters. A new branch of jet equilibria is revealed that is not contiguously connected with the G-L branch. This new branch at weak supercriticalities involves jets with larger amplitude compared to the ones of the G-L branch. Furthermore, this new branch continues even for subcritical values with respect to the linear flow-forming instability. Thus, a new nonlinear flow-forming instability out of homogeneous turbulence is revealed. Second, we investigate how both the linear flow-forming instability and the novel nonlinear flow-forming instability are equilibrated. We identify the physical processes underlying the jet equilibration as well as the types of eddies that contribute in each process. Third, we propose a modification of the diffusion coefficient of the G-L dynamics that is able to capture the asymmetric evolution for weak jets at scales other than the marginal scale (side-band instabilities) for the linear flow-forming instability.

We study the waves at the interface between two thin horizontal layers of immiscible liquids subject to high-frequency tangential vibrations. Nonlinear governing equations are derived for the cases of two- and three-dimensional flows and arbitrary ratio of layer thicknesses. The derivation is performed within the framework of the long-wavelength approximation, which is relevant as the linear instability of a thin-layers system is long-wavelength. The dynamics of equations is integrable and the equations themselves can be compared to the Boussinesq equation for the gravity waves in shallow water, which allows one to compare the action of the vibrational field to the action of the gravity and its possible effective inversion.

We address small volume-fraction asymptotic properties of a nonlocal isoperimetric functional with a confinement term, derived as the sharp interface limit of a variational model for self-assembly of diblock copolymers under confinement by nanoparticle inclusion. We introduce a small parameter $\eta$ to represent the size of the domains of the minority phase, and study the resulting droplet regime as $\eta\to 0$. By considering confinement densities which are spatially variable and attain a nondegenerate maximum, we present a two-stage asymptotic analysis wherein a separation of length scales is captured due to competition between the nonlocal repulsive and confining attractive effects in the energy. A key role is played by a parameter $M$ which gives the total volume of the droplets at order $\eta^3$ and its relation to existence and non-existence of Gamow's Liquid Drop model on $\mathbb{R}^3$. For large values of $M$, the minority phase splits into several droplets at an intermediate scale $\eta^{1/3}$, while for small $M$ minimizers form a single droplet converging to the maximum of the confinement density.

Gene regulatory network (GRN)-based morphogenetic models have recently gained an increasing attention. However, the relationship between microscopic properties of intracellular GRNs and macroscopic properties of morphogenetic systems has not been fully understood yet. Here we propose a theoretical morphogenetic model representing an aggregation of cells, and reveal the relationship between criticality of GRNs and morphogenetic pattern formation. In our model, the positions of the cells are determined by spring-mass-damper kinetics. Each cell has an identical Kauffman's $NK$ random Boolean network (RBN) as its GRN. We varied the properties of GRNs from ordered, through critical, to chaotic by adjusting node in-degree $K$. We randomly assigned four cell fates to the attractors of RBNs for cellular behaviors. By comparing diverse morphologies generated in our morphogenetic systems, we investigated what the role of the criticality of GRNs is in forming morphologies. We found that nontrivial spatial patterns were generated most frequently when GRNs were at criticality. Our finding indicates that the criticality of GRNs facilitates the formation of nontrivial morphologies in GRN-based morphogenetic systems.

The effect of spatial localization of states in distributed parameter systems under frozen parametric disorder is well known as the Anderson localization and thoroughly studied for the Schr\"odinger equation and linear dissipation-free wave equations. Some similar (or mimicking) phenomena can occur in dissipative systems such as the thermal convection ones. Specifically, many of these dissipative systems are governed by a modified Kuramoto-Sivashinsky equation, where the frozen spatial disorder of parameters has been reported to lead to excitation of localized patterns. Imposed advection in the modified Kuramoto-Sivashinsky equation can affect the localized patterns in a nontrivial way; it changes the localization properties and suppresses the pattern. The latter effect is considered in this paper by means of both numerical simulation and model reduction, which turns out to be useful for a comprehensive understanding of the bifurcation scenarios in the system. Two possible bifurcation scenarios of advective suppression ("washing-out") of localized patterns are revealed and characterised.

We study dynamics of Dirac solitons in prototypical networks modeling them by the nonlinear Dirac equation on metric graphs. Soliton solutions of the nonlinear Dirac equation on simple metric graphs are obtained. It is shown that these solutions provide reflectionless vertex transmission of the Dirac solitons under suitable conditions. The constraints for bond nonlinearity coefficients, allowing reflectionless transmission over a Y-junction are derived. The analytical results are confirmed by direct numerical simulations.

Traveling localized spots represent an important class of self-organized two-dimensional patterns in reaction-diffusion systems. We study open-loop control intended to guide a stable spot along a desired trajectory with desired velocity. Simultaneously, the spot's concentration profile does not change under control. For a given protocol of motion, we first express the control signal analytically in terms of the Goldstone modes and the propagation velocity of the uncontrolled spot. Thus, detailed information about the underlying nonlinear reaction kinetics is unnecessary. Then, we confirm the optimality of this solution by demonstrating numerically its equivalence to the solution of a regularized, optimal control problem. To solve the latter, the analytical expressions for the control are excellent initial guesses speeding-up substantially the otherwise time-consuming calculations.

The dynamics of single carrier wavepackets in nonlinear wave problems over periodic structures can be often formally approximated by the constant coefficient nonlinear Schr\"odinger equation (NLS) as an effective model for the wavepacket envelope. We provide a detailed proof of this approximation result for the Gross-Pitaevskii equation (GP) and a semilinear wave equation, both with periodic coefficients in $\mathbb{N}\ni d$ spatial dimensions and with cubic nonlinearities. The proof is carried out in Bloch expansion variables with estimates in an $L^1$-type norm, which translates to an estimate of the supremum norm of the error. The regularity required from the periodic coefficients in order to ensure a small residual and a small error is discussed. We also present a numerical example in two spatial dimensions confirming the approximation result and presenting an approximate traveling solitary wave in the GP with periodic coefficients.

We study asymptotic solutions to a singularly-perturbed, period-2 Toda lattice and use exponential asymptotics to examine `nanoptera', which are nonlocal solitary waves with constant-amplitude, exponentially small wave trains. With this approach, we isolate the exponentially small, constant-amplitude waves, and we elucidate the dynamics of these waves in terms of the Stokes phenomenon. We find a simple asymptotic expression for the waves, and we study configurations in which these waves vanish, producing localized solitary-wave solutions. In the limit of small mass ratio, we derive a simple anti-resonance condition for the manifestation these wave-free solutions.

We report the results of a new method of measuring the capillary length for the dendritic crystal growth of non-faceted materials. This method uses a nearly spherical crystal held near unstable equilibrium in an oscillating temperature field. For the growth of ammonium chloride crystals from aqueous solution, previous published estimates of the capillary length varied by over a factor of 20. With this new method, we find the product of the diffusion constant and capillary length $D d_0$ to be $0.78 \pm 0.07 \,\mathrm{\mu m^3/s}$, similar to that obtained for ammonium bromide crystals.

We present a computational study of the hydrodynamic coarsening in 3D of a critical mixture using the Cahn-Hilliard/Navier-Stokes model. The topology of the resulting intricate bicontinuous microstructure is analyzed through the principal curvatures to prove self-similar morphological evolution. We find that the self similarity exists for both systems: iso-viscous and with variable viscosity. However the two system have distinct topological character. Our simulations confirm that the predicted viscous growth regime exists in both cases. Moreover the coarsening rate is inversely proportional to an \textit{effective viscosity} that is the geometrical average of the viscosities of the two phases.

We predict a novel selection phenomenon in nonlinear interface dynamics out of equilibrium. Using a recently developed formalism based on the Schottky-Klein prime functions, we extended the existing integrable theory from a single interface to multiple moving interfaces. After applying this extended theory to the two-dimensional Laplacian growth, we derive a new rich class of exact (non-singular) solutions for the unsteady dynamics of an arbitrary assembly of air bubbles within a layer of a viscous fluid in a Hele-Shaw cell. These solutions demonstrate that all bubbles reach an asymptotic velocity, $U$, which is {\it precisely twice} greater than the velocity, $V$, of the uniform background flow, i.e., $U=2V$. The result does not depend on the number of bubbles. It is worth to mention that contrary to common belief, the predicted velocity selection does not require surface tension.

We study the scattering of the $\varphi^8$ kinks off each other, namely, we consider those $\varphi^8$ kinks that have power-law asymptotics. The slow power-law fall-off leads to a long-range interaction between the kink and the antikink. We investigate how the scattering scenarios depend on the initial velocities of the colliding kinks. In particular, we observe the `escape windows' -- the escape of the kinks after two or more collisions, explained by the resonant energy exchange between the translational and vibrational modes. In order to elucidate this phenomenon, we also analyze the excitation spectra of a solitary kink and of a composite kink+antikink configuration.