Search

Theories of localised pattern formation are important to understand a broad range of natural patterns, but are less well-understood than more established mechanisms of domain-filling pattern formation. Here, we extend recent work on pattern localisation via slow spatial heterogeneity in reaction-diffusion systems to the Swift-Hohenberg equation. We use a WKB asymptotic approach to show that, in the limit of a large domain and slowly varying heterogeneity, conditions for Turing-type linear instability localise in a simple way, with the spatial variable playing the role of a parameter. For nonlinearities locally corresponding to supercritical bifurcations in the spatially homogeneous system, this analysis asymptotically predicts regions where patterned states are confined, which we confirm numerically. We resolve the inner region of this asymptotic approach, finding excellent agreement with the tails of these confined pattern regions. In the locally subcritical case, however, this theory is insufficient to fully predict such confined regions, and so we propose an approach based on numerical continuation of a local homogeneous analog system. Pattern localisation in the heterogeneous system can then be determined based on the Maxwell point of this system, with the spatial variable parameterizing this point. We compare this theory of localisation via spatial heterogeneity to localised patterns arising from homoclinic snaking, and suggest a way to distinguish between different localisation mechanisms in natural systems based on how these structures decay to the background state (i.e. how their tails decay). We also explore cases where both of these local theories of pattern formation fail to capture the interaction between spatial heterogeneity and underlying pattern-forming mechanisms, suggesting that more work needs to be done to fully disentangle exogenous and intrinsic heterogeneity., Comment: 36 pages, 17 figures

General conditions are established under which reaction-cross-diffusion systems can undergo spatiotemporal pattern-forming instabilities. Recent work has focused on designing systems theoretically and experimentally to exhibit patterns with specific features, but the case of non-diagonal diffusion matrices has yet to be analysed. Here, a framework is presented for the design of general $n$-component reaction-cross-diffusion systems that exhibit Turing and wave instabilities of a given wavelength. For a fixed set of reaction kinetics, it is shown how to choose diffusion matrices that produce each instability; conversely, for a given diffusion tensor, how to choose linearised kinetics. The theory is applied to several examples including a hyperbolic reaction-diffusion system, two different 3-component models, and a spatio-temporal version of the Ross-Macdonald model for the spread of malaria., Comment: 49 pages, 14 figures

Weakly nonlinear amplitude equations are derived for the onset of spatially extended patterns on a general class of $n$-component bulk-surface reaction-diffusion systems in a ball, under the assumption of linear kinetics in the bulk. Linear analysis shows conditions under which various pattern modes can become unstable to either generalised pitchfork or transcritical bifurcations depending on the parity of the spatial wavenumber. Weakly nonlinear analysis is used to derive general expressions for the multi-component amplitude equations of different patterned states. These reduced-order systems are found to agree with prior normal forms for pattern formation bifurcations with $O(3)$ symmetry and provide information on the stability of bifurcating patterns of different symmetry types. The analysis is complemented with numerical results using a dedicated finite-element method. The theory is illustrated in two examples; a bulk-surface version of the Brusselator, and a four-component cell-polarity model.

A detailed set of experiments are described that capture over a 1000 different instances of the bounce of a golf ball. Video analysis is used to capture velocity and spin immediately prior to and subsequent to each bounce for a wide variety of landing conditions. Data are presented from two different turfs; one artificial and one from a typical tee. Measurement errors and repeatability are analysed. The data are compared to predictions from models of rigid bounce with friction, including Penner's modification to account for elasto-plasticity. Coefficients of restitution and friction, and Penner's effective contact angle are fit from the data. A better fit to the data is found using a non-physical piecewise-affine landing to lift-off relationship, which distinguishes between cases that bounce in pure slip from those that undergo rolling. Nevertheless, even balls that undergo rolling are typically found to lift off slipping, having undergone spin reversal. The findings suggest that further effort needs to be spent on finding simple physics-based models of golf-ball bounce., Comment: 15 pages, 9 figures, 3 tables, 2 data files

Inspired by the turf-ball interaction in golf, this paper seeks to understand the bounce of a ball that can be modelled as a rigid sphere and the surface as supplying an elasto-plastic contact force in addition to Coulomb friction. A general formulation is proposed that models the finite time interval of bounce from touch-down to lift-off. Key to the analysis is understanding transitions between slip and roll during the bounce. Starting from the rigid-body limit with a an energetic or Poisson coefficient of restitution, it is shown that slip reversal during the contact phase cannot be captured in this case, which result generalises to the case of pure normal compliance. Yet, the introduction of linear tangential stiffness and damping, does enable slip reversal. This result is extended to general weakly nonlinear normal and tangential compliance. An analysis using Filippov theory of piecewise-smooth systems leads to an argument in a natural limit that lift-off while rolling is non-generic and that almost all trajectories that lift off, do so under slip conditions. Moreover, there is a codimension-one surface in the space of incoming velocity and spin which divides balls that lift off with backspin from those that lift off with topspin. The results are compared with recent experimental measurements on golf ball bounce and the theory is shown to capture the main features of the data., Comment: 22 pages, 13 figures

The benefits, both in terms of productivity and public health, are investigated for different levels of engagement with the test, trace and isolate procedures in the context of a pandemic in which there is little or no herd immunity. Simple mathematical modelling is used in the context of a single, relatively closed workplace such as a factory or back-office where, in normal operation, each worker has lengthy interactions with a fixed set of colleagues. A discrete-time SEIR model on a fixed interaction graph is simulated with parameters that are motivated by the recent COVID-19 pandemic in the UK during a post-peak phase, including a small risk of viral infection from outside the working environment. Two kinds of worker are assumed, transparents who regularly test, share their results with colleagues and isolate as soon as a contact tests positive for the disease, and opaques who do none of these. Moreover, the simulations are constructed as a ``playable model" in which the transparency level, disease parameters and mean interaction degree can be varied by the user. The model is analysed in the continuum limit. All simulations point to the double benefit of transparency in maximising productivity and minimising overall infection rates. Based on these findings, public policy implications are discussed on how to incentivise this mutually beneficial behaviour in different kinds of workplace, and simple recommendations are made.

An investigation is undertaken of coupled reaction-diffusion systems in one spatial dimension that are able to support, in different regions of their parameter space, either an isolated spike solution, or stable localized patterns with an arbitrary number of peaks. The distinction between the two cases is drawn through the behavior of the far field, where there is either an oscillatory or a monotonic decay. Several examples are studied, including the Lugiato-Lefever model, a generalized Schakenberg system that arises in cellular-level morphogensis and a continuum model of urban crime spread. In each, it is found that localized patterns connected via a so-called homoclinic snaking curve in parameter space transition into a single spike solution as a second parameter is varied, via a change in topology of the snake into a series of disconnected branches. In each case, the transition is caused by a so-called Belyakov-Devaney transition between complex and real spatial eigenvalues of the fair field of the primary pulse. A codimension-two problem is studied in detail where a non-transverse homoclinic orbit undergoes this transition. A Shilnikov-style analysis is undertaken which reveals the asymptotics of how the infinite family of folds of multi-pulse orbits are all destroyed at the same parameter value. The results are shown to be consistent with numerical experiments on two of the examples.

The emergence of localised vibrations in cyclic and symmetric rotating structures, such as bladed disks of aircraft engines, has challenged engineers in the past few decades. In the linear regime, localised states may arise due to a lack of symmetry, as for example induced by inhomogeneities. However, when structures deviate from the linear behaviour, e.g. due to material nonlinearities, geometric nonlinearities like large deformations, or other nonlinear elements like joints or friction interfaces, localised states may arise even in perfectly symmetric structures. In this paper, a system consisting of coupled Duffing oscillators with linear viscous damping is subjected to external travelling wave forcing. The system may be considered a minimal model for bladed disks in turbomachinery operating in the nonlinear regime, where such excitation may arise due to imbalance or aerodynamic excitation. We demonstrate that near the resonance, in this non-conservative regime, localised vibration states bifurcate from the travelling waves. Complex bifurcation diagrams result, comprising stable and unstable dissipative solitons. The localised solutions can also be continued numerically to a conservative limit, where solitons bifurcate from the backbone curves of the travelling waves at finite amplitudes.

Observed bimodal tree cover distributions at particular environmental conditions and theoretical models indicate that some areas in the tropics can be in either of the alternative stable vegetation states forest or savanna. However, when including spatial interaction in nonspatial differential equation models of a bistable quantity, only the state with the lowest potential energy remains stable. Our recent reaction-diffusion model of Amazonian tree cover confirmed this and was able to reproduce the observed spatial distribution of forest versus savanna satisfactorily when forced by heterogeneous environmental and anthropogenic variables, even though bistability was underestimated. These conclusions were solely based on simulation results. Here, we perform an analytical and numerical analysis of the model. We derive the Maxwell point (MP) of the homogeneous reaction-diffusion equation without savanna trees as a function of rainfall and human impact and show that the front between forest and nonforest settles at this point as long as savanna tree cover near the front remains sufficiently low. For parameters resulting in higher savanna tree cover near the front, we also find irregular forest-savanna cycles and woodland-savanna bistability, which can both explain the remaining observed bimodality., Comment: 26 pages, 6 figures, 2 tables, supplementary info included

This paper analyses in detail the dynamics in a neighbourhood of a G\'enot-Brogliato point, colloquially termed the G-spot, which physically represents so-called dynamic jam in rigid body mechanics with unilateral contact and Coulomb friction. Such singular points arise in planar rigid body problems with slipping point contacts at the intersection between the conditions for onset of lift-off and for the Painlev\'e paradox. The G-spot can be approached in finite time by an open set of initial conditions in a general class of problems. The key question addressed is what happens next. In principle trajectories could, at least instantaneously, lift off, continue in slip, or undergo a so-called impact without collision. Such impacts are non-local in momentum space and depend on properties evaluated away from the G-spot. The results are illustrated on a particular physical example, namely the a frictional impact oscillator first studied by Leine et al. The answer is obtained via an analysis that involves a consistent contact regularisation with a stiffness proportional to $1/\varepsilon^2$. Taking a singular limit as $\varepsilon \to 0$, one finds an inner and an outer asymptotic zone in the neighbourhood of the G-spot. Two distinct cases are found according to whether the contact force becomes infinite or remains finite as the G-spot is approached. In the former case it is argued that there can be no such canards and so an impact without collision must occur. In the latter case, the canard trajectory acts as a dividing surface between trajectories that momentarily lift off and those that do not before taking the impact. The orientation of the initial condition set leading to each eventuality is shown to change each time a certain positive parameter $\beta$ passes through an integer.

A system of two Schnakenberg-like reaction-diffusion equations is investigated analytically and numerically. The system has previously been used as a minimal model for concentrations of GTPases involved in the process of cell polarisation. Source and loss terms are added, breaking the mass conservation, which was shown previously to be responsible for the generation of stable fronts via a so-called wave-pinning mechanism. The extended model gives rise to a unique homogeneous equilibrium in the parameter region of interest, which loses stability via a pattern formation, or Turing bifurcation. The bistable character of the reaction terms ensures that this bifurcation is subcrtical for sufficiently small values of the driving parameter multiplying the nonlinear kinetics. This subcriticality leads to the onset of a multitude of localised solutions, through the homoclinic snaking mechanism. As the driving parameter is further decreased, the multitude of solutions transforms into a single pulse through a Belyakov-Devaney transition in which there is the loss of a precursive pattern. An asymptotic analysis is used to probe the conservative limit in which the source and loss terms vanish. Matched asymptotic analysis shows that on an infinite domain the pulse solution transitions into a pair of fronts, with an additional weak quadratic core and exponential tails. On a finite domain, the core and tails disappear, leading to the mere wave-pinning front and its mirror image.

The 120-year old so-called Painleve paradox involves the loss of determinism in models of planar rigid bodies in point contact with a rigid surface, subject to Coulomb-like dry friction. The phenomenon occurs due to coupling between normal and rotational degrees-of-freedom such that the effective normal force becomes attractive rather than repulsive. Despite a rich literature, the forward evolution problem remains unsolved other than in certain restricted cases in 2D with single contact points. Various practical consequences of the theory are revisited, including models for robotic manipulators, and the strange behaviour of chalk when pushed rather than dragged across a blackboard. Reviewing recent theory, a general formulation is proposed, including a Poisson or energetic impact law. The general problem in 2D with a single point of contact is discussed and cases or inconsistency or indeterminacy enumerated. Strategies to resolve the paradox via contact regularisation are discussed from a dynamical systems point of view. By passing to the infinite stiffness limit and allowing impact without collision, inconsistent and indeterminate cases are shown to be resolvable for all open sets of conditions. However, two unavoidable ambiguities that can be reached in finite time are discussed in detail, so called dynamic jam and reverse chatter. A partial review is given of 2D cases with two points of contact showing how a greater complexity of inconsistency and indeterminacy can arise. Extension to fully three-dimensional analysis is briefly considered and shown to lead to further possible singularities. In conclusion, the ubiquity of the \pain paradox is highlighted and open problems are discussed.

A simple nonlinear transmission-line model of the cochlea with longitudinal coupling is introduced that can reproduce Basilar membrane response and neural tuning in the chinchilla. It is found that the middle ear has little effect on cochlear resonances, and hence conclude that the theory of coherent reflections is not applicable to the model. The model also provides an explanation of the emergence of spontaneous otoacoustic emissions (SOAEs). It is argued that SOAEs arise from Hopf bifurcations of the transmission-line model and not from localized instabilities. The paper shows that emissions can become chaotic, intermittent and fragile to perturbations., Comment: 9 pages, 4 figures

An analytic approximation for the critical clearing time (CCT) metric is derived from direct methods for power system stability. The formula has been designed to incorporate as many features of transient stability analysis as possible such as different fault locations and different post-fault network states. The purpose of this metric is to analyse trends in stability (in terms of CCT) of power systems under the variation of a system parameter. We demonstrate the performance of this metric to measure stability trends on an aggregated power network, the so-called two machine infinite bus network, by varying load parameters in the full bus admittance matrix using numerical continuation. Our metric is compared to two other expressions for the CCT which incorporate additional non-linearities present in the model., Comment: 11 pages

The well-studied Hindmarsh-Rose model of neural action potential is revisited from the point of view of global bifurcation analysis. This slow-fast system of three paremeterised differential equations is arguably the simplest reduction of Hodgkin-Huxley models capable of exhibiting all qualitatively important distinct kinds of spiking and bursting behaviour. First, keeping the singular perturbation parameter fixed, a comprehensive two-parameter bifurcation diagram is computed by brute force. Of particular concern is the parameter regime where lobe-shaped regions of irregular bursting undergo a transition to stripe-shaped regions of periodic bursting. The boundary of each stripe represents a fold bifurcation that causes a smooth spike-adding transition where the number of spikes in each burst is increased by one. Next, numerical continuation studies reveal that the global structure is organised by various curves of homoclinic bifurcations. In particular the lobe to stripe transition is organised by a sequence of codimension-two orbit- and inclination-flip points that occur along {\em each} homoclinic branch. Each branch undergoes a sharp turning point and hence approximately has a double-cover of the same curve in parameter space. The sharp turn is explained in terms of the interaction between a two-dimensional unstable manifold and a one-dimensional slow manifold in the singular limit. Finally, a new local analysis is undertaken using approximate Poincar\'{e} maps to show that the turning point on each homoclinic branch in turn induces an inclination flip that gives birth to the fold curve that organises the spike-adding transition. Implications of this mechanism for explaining spike-adding behaviour in other excitable systems are discussed., Comment: 32 pages, 18 figures, submitted to SIAM Journal on Applied Dynamical Systems

We address the existence and properties of discrete embedded solitons (ESs), i.e., localized waves existing inside the phonon band in a nonlinear dynamical-lattice model. The model describes a one-dimensional array of optical waveguides with both the quadratic (second-harmonic generation) and cubic nonlinearities. A rich family of ESs was previously known in the continuum limit of the model. First, a simple motivating problem is considered, in which the cubic nonlinearity acts in a single waveguide. An explicit solution is constructed asymptotically in the large-wavenumber limit. The general problem is then shown to be equivalent to the existence of a homoclinic orbit in a four-dimensional reversible map. From properties of such maps, it is shown that (unlike ordinary gap solitons), discrete ESs have the same codimension as their continuum counterparts. A specific numerical method is developed to compute homoclinic solutions of the map, that are symmetric under a specific reversing transformation. Existence is then studied in the full parameter space of the problem. Numerical results agree with the asymptotic results in the appropriate limit and suggest that the discrete ESs may be semi-stable as in the continuous case., Comment: A revtex4 text file and 51 eps figure files. To appear in Nonlinearity

An explanation is given for previous numerical results which suggest a certain bifurcation of `vector solitons' from scalar (single-component) solitary waves in coupled nonlinear Schroedinger (NLS) systems. The bifurcation in question is nonlocal in the sense that the vector soliton does not have a small-amplitude component, but instead approaches a solitary wave of one component with two infinitely far-separated waves in the other component. Yet, it is argued that this highly nonlocal event can be predicted from a purely local analysis of the central solitary wave alone. Specifically the linearisation around the central wave should contain asymptotics which grow at precisely the speed of the other-component solitary waves on the two wings. This approximate argument is supported by both a detailed analysis based on matched asymptotic expansions, and numerical experiments on two example systems. The first is the usual coupled NLS system involving an arbitrary ratio between the self-phase and cross-phase modulation terms, and the second is a coupled NLS system with saturable nonlinearity that has recently been demonstrated to support stable multi-peaked solitary waves. The asymptotic analysis further reveals that when the curves which define the proposed criterion for scalar nonlocal bifurcations intersect with boundaries of certain local bifurcations, the nonlocal bifurcation could turn from scalar to non-scalar at the intersection. This phenomenon is observed in the first example. Lastly, we have also selectively tested the linear stability of several solitary waves just born out of scalar nonlocal bifurcations. We found that they are linearly unstable. However, they can lead to stable solitary waves through parameter continuation., Comment: To appear in Nonlinearity

We demonstrate that weak parametric interaction of a fundamental beam with its third harmonic field in Kerr media gives rise to a rich variety of families of non-fundamental (multi-humped) solitary waves. Making a comprehensive comparison between bifurcation phenomena for these families in bulk media and planar waveguides, we discover two novel types of soliton bifurcations and other interesting findings. The later includes (i) multi-humped solitary waves without even or odd symmetry and (ii) multi-humped solitary waves with large separation between their humps which, however, may not be viewed as bound states of several distinct one-humped solitons., Comment: 9 pages, 17 figures, submitted to Phys. Rev. E

The first theoretical results are reported predicting {\em moving} solitons residing inside ({\it embedded} into) the continuous spectrum of radiation modes. The model taken is a Bragg-grating medium with Kerr nonlinearity and additional second-derivative (wave) terms. The moving embedded solitons (ESs) are doubly isolated (of codimension 2), but, nevertheless, structurally stable. Like quiescent ESs, moving ESs are argued to be stable to linear approximation, and {\it semi}-stable nonlinearly. Estimates show that moving ESs may be experimentally observed as $\sim$10 fs pulses with velocity $\leq 1/10$th that of light., Comment: 9 pages 2 figures

We show that soliton excitation in intracavity downconversion naturally selects a strictly defined frequency difference between the signal and idler fields. In particular, this phenomenon implies that if the signal has smaller losses than the idler then its frequency is pulled away from the cavity resonance and the idler frequency is pulled towards the resonance and {\em vice versa}. The frequency selection is shown to be closely linked with the relative energy balance between the idler and signal fields., Comment: 5 pages, 3 figures. To appear in Phys Rev Lett

We develop {\em the first analytical theory of multikinks} for strongly {\em dispersive nonlinear systems}, considering the examples of the weakly discrete sine-Gordon model and the generalized Frenkel-Kontorova model with a piecewise parabolic potential. We reveal that there are no $2\pi$-kinks for this model, but there exist {\em discrete sets} of $2\pi N$-kinks for all N>1. We also show their bifurcation structure in driven damped systems., Comment: 4 pages 5 figures. To appear in Phys Rev E

We report a rich spectrum of isolated solitons residing inside ({\it embedded } into) the continuous radiation spectrum in a simple model of three-wave spatial interaction in a second-harmonic-generating planar optical waveguide equipped with a quasi-one-dimensional Bragg grating. An infinite sequence of fundamental embedded solitons are found, which differ by the number of internal oscillations. Branches of these zero-walkoff spatial solitons give rise, through bifurcations, to several secondary branches of walking solitons. The structure of the bifurcating branches suggests a multistable configuration of spatial optical solitons, which may find straightforward applications for all-optical switching., Comment: 5 pages 5 figures. To appear in Phys Rev E

The effect of dispersion or diffraction on zero-velocity solitons is studied for the generalized massive Thirring model describing a nonlinear optical fiber with grating or parallel-coupled planar waveguides with misaligned axes. The Thirring solitons existing at zero dispersion/diffraction are shown numerically to be separated by a finite gap from three isolated soliton branches. Inside the gap, there is an infinity of multi-soliton branches. Thus, the Thirring solitons are structurally unstable. In another parameter region (far from the Thirring limit), solitons exist everywhere., Comment: 12 pages, Latex. To appear in Phys. Rev. Lett