
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 Faradayinstability 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 Weberequation 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 11dimensional spacetimes to 4dimensional 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
integrabilitybreaking 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 spacetime dependent neural field which
describes the average membrane potential of neurons in a single layer. This
neural field is modelled by a dissipative integrodifferential equation, the
socalled 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 wavevectorfreqency 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 spacetime 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
reactiondiffusion 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 reactiondiffusion equation.

This paper proposes a control theoretic framework to model and analyze the
selforganized 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 celltocell communication mechanism
mediated by a diffusible signaling molecule, thereby the dynamics of molecular
concentrations are approximately modeled as a reactiondiffusion system with a
single diffuser. We first introduce a feedback model representation of the
reactiondiffusion 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
selforganized spatial pattern formation, or Turing pattern formation, of the
bionanomachines. We propose a novel synthetic biocircuit motif called
activatorrepressordiffuser system and show that it is one of the minimum
biomolecular circuits that admit selforganized 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 fourdimensional hypersurface or hyperfilm in
fivedimensional spacetime. 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
twodimensional 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 highfrequency 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 Hertzlike
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 quasiparticle collisions. Next we consider wave
propagation along Hertzlike granular chains of beads in the presence of
disorder and periodicity in the interaction exponents present in the Hertzlike
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 GrossPitaevskii equation in one dimension is considered with
a complex periodic potential satisfying the conditions of the PT (paritytime
reversal) symmetry. Under rather general assumptions on the potentials we prove
bifurcations of PTsymmetric nonlinear bound states from the end points of a
real interval in the spectrum of the nonselfadjoint linear Schrodinger
operator with a complex PTsymmetric 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 Kochtype
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 noisedriven 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
stepwise adjustable IFS to introduce some variance (such as
nonselfsimilarity) 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 (Cantorlike
sets, Sierpinskilike Triangles, etc) with the calculations of their
corresponding dimensions are studied.

By studying the coarsening dynamics of a onedimensional spin1 BoseHubbard
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 powerlaw 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 flowforming 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 flowforming instability follows GinzburgLandau (GL)
dynamics. Here, the SSD framework is used to study the equilibration of this
flowforming instability for small supercriticality. First, we compare the
predictions of the weakly nonlinear GL 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 GL
branch. This new branch at weak supercriticalities involves jets with larger
amplitude compared to the ones of the GL branch. Furthermore, this new branch
continues even for subcritical values with respect to the linear flowforming
instability. Thus, a new nonlinear flowforming instability out of homogeneous
turbulence is revealed. Second, we investigate how both the linear flowforming
instability and the novel nonlinear flowforming 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 GL dynamics that is able to
capture the asymmetric evolution for weak jets at scales other than the
marginal scale (sideband instabilities) for the linear flowforming
instability.

We study the waves at the interface between two thin horizontal layers of
immiscible liquids subject to highfrequency tangential vibrations. Nonlinear
governing equations are derived for the cases of two and threedimensional
flows and arbitrary ratio of layer thicknesses. The derivation is performed
within the framework of the longwavelength approximation, which is relevant as
the linear instability of a thinlayers system is longwavelength. 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 volumefraction asymptotic properties of a nonlocal
isoperimetric functional with a confinement term, derived as the sharp
interface limit of a variational model for selfassembly 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 twostage 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 nonexistence 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 springmassdamper
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 indegree $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 GRNbased 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 dissipationfree
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 KuramotoSivashinsky equation,
where the frozen spatial disorder of parameters has been reported to lead to
excitation of localized patterns. Imposed advection in the modified
KuramotoSivashinsky 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 ("washingout") 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 Yjunction are
derived. The analytical results are confirmed by direct numerical simulations.

Traveling localized spots represent an important class of selforganized
twodimensional patterns in reactiondiffusion systems. We study openloop
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
speedingup substantially the otherwise timeconsuming 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 GrossPitaevskii 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 singularlyperturbed, period2 Toda
lattice and use exponential asymptotics to examine `nanoptera', which are
nonlocal solitary waves with constantamplitude, exponentially small wave
trains. With this approach, we isolate the exponentially small,
constantamplitude 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
solitarywave solutions. In the limit of small mass ratio, we derive a simple
antiresonance condition for the manifestation these wavefree solutions.

We report the results of a new method of measuring the capillary length for
the dendritic crystal growth of nonfaceted 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 CahnHilliard/NavierStokes model. The topology of
the resulting intricate bicontinuous microstructure is analyzed through the
principal curvatures to prove selfsimilar morphological evolution. We find
that the self similarity exists for both systems: isoviscous 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
SchottkyKlein prime functions, we extended the existing integrable theory from
a single interface to multiple moving interfaces. After applying this extended
theory to the twodimensional Laplacian growth, we derive a new rich class of
exact (nonsingular) solutions for the unsteady dynamics of an arbitrary
assembly of air bubbles within a layer of a viscous fluid in a HeleShaw 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 powerlaw asymptotics. The slow
powerlaw falloff leads to a longrange 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.