Showing papers in "arXiv: Pattern Formation and Solitons in 2010"

Slime mould logical gates: exploring ballistic approach

Abstract: University of the West of England, Bristol, United Kingdomandrew.adamatzky@uwe.ac.ukAbstractPlasmodium of Physarum polycephalum is a single cell visible by unaided eye. On anon-nutrient substrate the plasmodium propagates as a traveling localization, as acompact wave-fragment of protoplasm. The plasmodium-localization travels in itsoriginally predetermined direction for a substantial period of time even when nogradient of chemo-attractants is present. We utilize this property of Physarum lo-calizations to design a two-input two-output Boolean logic gates hx;yi!hxy;x+yiand hx;yi!hx;xyi. We verify the designs in laboratory experiments and computersimulations. We cascade the logical gates into one-bit half-adder and simulate itsfunctionality.Keywords: Physarum polycephalum, logical gate, adder, unconventional computer,chemical computer, biological computer

KAM tori in 1D random discrete nonlinear Schr\"odinger model?

Abstract: We suggest that KAM theory could be extended for certain infinite-dimensional systems with purely discrete linear spectrum. We provide empirical arguments for the existence of square summable infinite-dimensional invariant tori in the random discrete nonlinear Schr\"odinger equation, appearing with a finite probability for a given initial condition with sufficiently small norm. Numerical support for the existence of a fat Cantor set of initial conditions generating almost-periodic oscillations is obtained by analyzing (i) sets of recurrent trajectories over successively larger time scales, and (ii) finite-time Lyapunov exponents. The norm region where such KAM-like tori may exist shrinks to zero when the disorder strength goes to zero and the localization length diverges.

Can a traveling wave connect two unstable states? The case of the nonlocal Fisher equation

Abstract: This note investigates the properties of the traveling waves solutions of the nonlocal Fisher equation. The existence of such solutions has been proved recently in \cite{BNPR} but their asymptotic behavior was still unclear. We use here a new numerical approximation of these traveling waves which shows that some traveling waves connect the two homogeneous steady states $0$ and $1$, which is a striking fact since $0$ is dynamically unstable and $1$ is unstable in the sense of Turing.

On electrical correlates of Physarum polycephalum spatial activity: Can we see Physarum Machine in the dark?

Abstract: Plasmodium of Physarum polycephalum is a single cell visible by unaided eye, which spans sources of nutrients with its protoplasmic network. In a very simple experimental setup we recorded electric potential of the propagating plasmodium. We discovered a complex interplay of short range oscillatory behaviour combined with long range, low frequency oscillations which serve to communicate information between different parts of the plasmodium. The plasmodium's response to changing environmental conditions forms basis patterns of electric activity, which are unique indicators of the following events: plasmodium occupies a site, plasmodium functions normally, plasmodium becomes `agitated' due to drying substrate, plasmodium departs a site, and plasmodium forms sclerotium. Using a collective particle approximation of Physarum polycephalum we found matching correlates of electrical potential in computational simulations by measuring local population flux at the node positions, generating trains of high and low frequency oscillatory behaviour. Motifs present in these measurements matched the response `grammar' of the plasmodium when encountering new nodes, simulated consumption of nutrients, exposure to simulated hazardous illumination and sclerotium formation. The distributed computation of the particle collective was able to calculate beneficial network structures and sclerotium position by shifting the active growth zone of the simulated plasmodium. The results show future promise for the non-invasive study of the complex dynamical behaviour within --- and health status of --- living systems.

Coherently coupled bright optical solitons and their collisions

Abstract: We obtain explicit bright one- and two-soliton solutions of the integrable case of the coherently coupled nonlinear Schr{o}dinger equations by applying a non-standard form of the Hirota's direct method. We find that the system admits both degenerate and non-degenerate solitons in which the latter can take single hump, double hump, and flat-top profiles. Our study on the collision dynamics of solitons in the integrable case shows that the collision among degenerate solitons and also the collision of non-degenerate solitons are always standard elastic collisions. But the collision of a degenerate soliton with a non-degenerate soliton induces switching in the latter leaving the former unaffected after collision, thereby showing a different mechanism from that of the Manakov system.

Traveling plateaus for a hyperbolic Keller-Segel system with attraction and repulsion: existence and branching instabilities

Abstract: How can repulsive and attractive forces, acting on a conservative system, create stable traveling patterns or branching instabilities? We have proposed to study this question in the framework of the hyperbolic Keller-Segel system with logistic sensitivity. This is a model system motivated by experiments on cell communities auto-organization, a field which is also called socio-biology. We continue earlier modeling work, where we have shown numerically that branching patterns arise for this system and we have analyzed this instability by formal asymptotics for small diffusivity of the chemo-repellent. Here we are interested in the more general situation, where the diffusivities of both the chemo-attractant and the chemo-repellent are positive. To do so, we develop an appropriate functional analysis framework. We apply our method to two cases. Firstly we analyze steady states. Secondly we analyze traveling waves when neglecting the degradation coefficient of the chemo-repellent; the unique wave speed appears through a singularity cancelation which is the main theoretical difficulty. This shows that in different situations the cell density takes the shape of a plateau. The existence of steady states and traveling plateaus are a symptom of how rich the system is and why branching instabilities can occur. Numerical tests show that large plateaus may split into smaller ones, which remain stable.

On modelling transitional turbulent flows using under-resolved direct numerical simulations: The case of plane Couette flow

Abstract: Direct numerical simulations have proven of inestimable help to our understanding of the transition to turbulence in wall-bounded flows. While the dynamics of the transition from laminar flow to turbulence via localised spots can be investigated with reasonable computing resources in domains of limited extent, the study of the decay of turbulence in conditions approaching those in the laboratory requires consideration of domains so wide as to exclude the recourse to fully resolved simulations. Using Gibson's C++ code ChannelFlow, we scrutinize the effects of a controlled lowering of the numerical resolution on the decay of turbulence in plane Couette flow at a quantitative level. We show that the number of Chebyshev polynomials describing the cross-stream dependence can be drastically decreased while preserving all the qualitative features of the solution. In particular, the oblique turbulent band regime experimentally observed in the upper part of the transitional range is extremely robust. In terms of Reynolds numbers, the resolution lowering is seen to yield a regular downward shift of the upper and lower thresholds Rt and Rg where the bands appear and break down. The study is illustrated with the results of two preliminary experiments.

Global existence for semilinear reaction-diffusion systems on evolving domains

Abstract: We present global existence results for solutions of reaction-diffusion systems on evolving domains. Global existence results for a class of reaction-diffusion systems on fixed domains are extended to the same systems posed on spatially linear isotropically evolving domains. The results hold without any assumptions on the sign of the growth rate. The analysis is valid for many systems that commonly arise in the theory of pattern formation. We present numerical results illustrating our theoretical findings.

On polymorphic logical gates in sub-excitable chemical medium

Abstract: In a sub-excitable light-sensitive Belousov-Zhabotinsky chemical medium an asymmetric disturbance causes the formation of localized traveling wave-fragments. Under the right conditions these wave-fragment can conserve their shape and velocity vectors for extended time periods. The size and life span of a fragment depend on the illumination level of the medium. When two or more wave-fragments collide they annihilate or merge into a new wave-fragment. In computer simulations based on the Oregonator model we demonstrate that the outcomes of inter-fragment collisions can be controlled by varying the illumination level applied to the medium. We interpret these wave-fragments as values of Boolean variables and design collision-based polymorphic logical gates. The gate implements operation XNOR for low illumination, and it acts as NOR gate for high illumination. As a NOR gate is a universal gate then we are able to demonstrate that a simulated light sensitive BZ medium exhibits computational universality.

Flow-induced channelization in a porous medium

Abstract: We propose a theory for erosional channelization induced by fluid flow in a saturated granular porous medium. When the local fluid flow-induced stress is larger than a critical threshold, grains are dislodged and carried away so that the porosity of the medium is altered by erosion. This in turn affects the local hydraulic conductivity and pressure in the medium and results in the growth and development of channels that preferentially conduct the flow. Our multiphase model involves a dynamical porosity field that evolves along with the volume fraction of the mobile and immobile grains in response to fluid flow that couples the spatiotemporal dynamics of the three phases. Numerical solutions of the resulting initial boundary value problem show how channels form in porous media and highlights how heterogeneity in the erosion threshold dictates the form of the patterns and thus the ability to control them.

Domain structure of bulk ferromagnetic crystals in applied fields near saturation

Abstract: We investigate the ground state of a uniaxial ferromagnetic plate with perpendicular easy axis and subject to an applied magnetic field normal to the plate. Our interest is the asymptotic behavior of the energy in macroscopically large samples near the saturation field. We establish the scaling of the critical value of the applied field strength below saturation at which the ground state changes from the uniform to a branched domain magnetization pattern and the leading order scaling behavior of the minimal energy. Furthermore, we derive a reduced sharp-interface energy giving the precise asymptotic behavior of the minimal energy in macroscopically large plates under a physically reasonable assumption of small deviations of the magnetization from the easy axis away from domain walls. On the basis of the reduced energy, and by a formal asymptotic analysis near the transition, we derive the precise asymptotic values of the critical field strength at which non-trivial minimizers (either local or global) emerge. The non-trivial minimal energy scaling is achieved by magnetization patterns consisting of long slender needle-like domains of magnetization opposing the applied field

Localized standing waves in inhomogeneous Schrodinger equations

Abstract: A nonlinear Schrodinger equation arising from light propagation down an inhomogeneous medium is considered. The inhomogeneity is reflected through a non-uniform coefficient of the non-linear term in the equation. In particular, a combination of self-focusing and self-defocusing nonlinearity, with the self-defocusing region localized in a finite interval, is investigated. Using numerical computations, the extension of linear eigenmodes of the corresponding linearized system into nonlinear states is established, particularly nonlinear continuations of the fundamental state and the first excited state. The (in)stability of the states is also numerically calculated, from which it is obtained that symmetric nonlinear solutions become unstable beyond a critical threshold norm. Instability of the symmetric states is then investigated analytically through the application of a topological argument. Determination of instability of positive symmetric states is reduced to simple geometric properties of the composite phase plane orbit of the standing wave. Further the topological argument is applied to higher excited states and instability is again reduced to straightforward geometric calculations. For a relatively high norm, it is observed that asymmetric states bifurcate from the symmetric ones. The stability and instability of asymmetric states is also considered.

Localized bioconvection of Euglena caused by phototaxis in the lateral direction

Abstract: Euglena, a swimming micro-organism, exhibited a characteristic bioconvection that was localized at the center of a sealed chamber under bright illumination to induce negative phototaxis. This localized pattern consisted of high-density spots, in which convection was found. These observations were reproduced by a mathematical model that was based on the phototaxis of individual cells in both the vertical and lateral directions. Our results indicate that this convection is maintained by upward swimming, as with general bioconvection, and the localization originates from lateral phototaxis.

Nonlinear waves in Newton's cradle and the discrete p-Schroedinger equation

Abstract: We study nonlinear waves in Newton's cradle, a classical mechanical system consisting of a chain of beads attached to linear pendula and interacting nonlinearly via Hertz's contact forces. We formally derive a spatially discrete modulation equation, for small amplitude nonlinear waves consisting of slow modulations of time-periodic linear oscillations. The fully-nonlinear and unilateral interactions between beads yield a nonstandard modulation equation that we call the discrete p-Schroedinger (DpS) equation. It consists of a spatial discretization of a generalized Schroedinger equation with p-Laplacian, with fractional p>2 depending on the exponent of Hertz's contact force. We show that the DpS equation admits explicit periodic travelling wave solutions, and numerically find a plethora of standing wave solutions given by the orbits of a discrete map, in particular spatially localized breather solutions. Using a modified Lyapunov-Schmidt technique, we prove the existence of exact periodic travelling waves in the chain of beads, close to the small amplitude modulated waves given by the DpS equation. Using numerical simulations, we show that the DpS equation captures several other important features of the dynamics in the weakly nonlinear regime, namely modulational instabilities, the existence of static and travelling breathers, and repulsive or attractive interactions of these localized structures.

Multibreather and vortex breather stability in Klein--Gordon lattices: Equivalence between two different approaches

Abstract: In this work, we revisit the question of stability of multibreather configurations, i.e., discrete breathers with multiple excited sites at the anti-continuum limit of uncoupled oscillators. We present two methods that yield quantitative predictions about the Floquet multipliers of the linear stability analysis around such exponentially localized in space, time-periodic orbits, based on the Aubry band method and the MacKay effective Hamiltonian method and prove that their conclusions are equivalent. Subsequently, we showcase the usefulness of the methods by a series of case examples including one-dimensional multi-breathers, and two-dimensional vortex breathers in the case of a lattice of linearly coupled oscillators with the Morse potential and in that of the discrete $\phi^4$ model.

Computational Modelling of Nonlinear Calcium Waves

Abstract: The calcium transport in biological systems is modelled as a reaction-diffusion process. Nonlinear calcium waves are then simulated using a stochastic cellular automaton whose rules are derived from the corresponding coupled partial differential equations. Numerical simulations show self-organized criticality in the complex calcium waves and patterns. Both the stochastic cellular automaton approach and the equation-based simulations can predict the characteristics of calcium waves and complex pattern formation. The implication of locality of calcium distribution with positional information in biological systems is also discussed.

Creating very slow optical gap solitons with inter-fiber coupling

Abstract: We show that gap-acoustic solitons, i.e., optical gap solitons with electrostrictive coupling to sound modes, can be produced with velocities down to less than 2.5% of the speed of light using a fiber Bragg grating that is linearly coupled to a non-Bragg fiber over a finite domain. Forward- and backward-moving light pulses in the non-Bragg fiber that reach the coupling region simultaneously couple into the Bragg fiber and form a moving soliton, which then propagates beyond the coupling region.

The large core limit of spiral waves in excitable media: A numerical approach

Abstract: We modify the freezing method introduced by Beyn & Thuemmler, 2004, for analyzing rigidly rotating spiral waves in excitable media. The proposed method is designed to stably determine the rotation frequency and the core radius of rotating spirals, as well as the approximate shape of spiral waves in unbounded domains. In particular, we introduce spiral wave boundary conditions based on geometric approximations of spiral wave solutions by Archimedean spirals and by involutes of circles. We further propose a simple implementation of boundary conditions for the case when the inhibitor is non-diffusive, a case which had previously caused spurious oscillations. We then utilize the method to numerically analyze the large core limit. The proposed method allows us to investigate the case close to criticality where spiral waves acquire infinite core radius and zero rotation frequency, before they begin to develop into retracting fingers. We confirm the linear scaling regime of a drift bifurcation for the rotation frequency and the core radius of spiral wave solutions close to criticality. This regime is unattainable with conventional numerical methods.

Collapse in the nonlocal nonlinear Schr\"odinger equation

Abstract: We discuss spatial dynamics and collapse scenarios of localized waves governed by the nonlinear Schr\"{o}dinger equation with nonlocal nonlinearity. Firstly, we prove that for arbitrary nonsingular attractive nonlocal nonlinear interaction in arbitrary dimension collapse does not occur. Then we study in detail the effect of singular nonlocal kernels in arbitrary dimension using both, Lyapunoff's method and virial identities. We find that for for a one-dimensional case, i.e. for $n=1$, collapse cannot happen for nonlocal nonlinearity. On the other hand, for spatial dimension $n\geq2$ and singular kernel $\sim 1/r^\alpha$, no collapse takes place if $\alpha<2$, whereas collapse is possible if $\alpha\ge2$. Self-similar solutions allow us to find an expression for the critical distance (or time) at which collapse should occur in the particular case of $\sim 1/r^2$ kernels. Moreover, different evolution scenarios for the three dimensional physically relevant case of Bose Einstein condensate are studied numerically for both, the ground state and a higher order toroidal state with and without an additional local repulsive nonlinear interaction. In particular, we show that presence of an additional local repulsive term can prevent collapse in those cases.

A primer of swarm equilibria

Abstract: We study equilibrium configurations of swarming biological organisms subject to exogenous and pairwise endogenous forces. Beginning with a discrete dynamical model, we derive a variational description of the corresponding continuum population density. Equilibrium solutions are extrema of an energy functional, and satisfy a Fredholm integral equation. We find conditions for the extrema to be local minimizers, global minimizers, and minimizers with respect to infinitesimal Lagrangian displacements of mass. In one spatial dimension, for a variety of exogenous forces, endogenous forces, and domain configurations, we find exact analytical expressions for the equilibria. These agree closely with numerical simulations of the underlying discrete model.The exact solutions provide a sampling of the wide variety of equilibrium configurations possible within our general swarm modeling framework. The equilibria typically are compactly supported and may contain $\delta$-concentrations or jump discontinuities at the edge of the support. We apply our methods to a model of locust swarms, which are observed in nature to consist of a concentrated population on the ground separated from an airborne group. Our model can reproduce this configuration; quasi-two-dimensionality of the model plays a critical role.

Approximating Mexican highways with slime mould

Abstract: Plasmodium of Physarum polycephalum is a single cell visible by unaided eye. During its foraging behavior the cell spans spatially distributed sources of nutrients with a protoplasmic network. Geometrical structure of the protoplasmic networks allows the plasmodium to optimize transport of nutrients between remote parts of its body. Assuming major Mexican cities are sources of nutrients how much structure of Physarum protoplasmic network correspond to structure of Mexican Federal highway network? To find an answer undertook a series of laboratory experiments with living Physarum polycephalum. We represent geographical locations of major cities by oat flakes, place a piece of plasmodium in Mexico city area, record the plasmodium's foraging behavior and extract topology of nutrient transport networks. Results of our experiments show that the protoplasmic network formed by Physarum is isomorphic, subject to limitations imposed, to a network of principle highways. Ideas and results of the paper may contribute towards future developments in bio-inspired road planning.

Time-delayed-feedback control of pattern formation in symmetrical bistable media

Abstract: Effects of time-delayed-feedback on pattern formation are studied in symmetrical bistable media. The results show that the time delay alters the behavior of the front bifurcation remarkably. The critical point of the Nonequilibrium Ising-Bloch bifurcation is shifted, which results from the competition between the feedbacks of the two variables. The stability of the front to transverse perturbation is also changed by the time delay, leading to a transition between spiral and labyrinthine patterns. One can control the pattern formation in bistable media by using time-delayed-feedback.

Breathing modes of long Josephson junctions with phase-shifts

Abstract: We consider a spatially inhomogeneous sine-Gordon equation with a time-periodic drive, modeling a microwave driven long Josephson junction with phase-shifts. Under appropriate conditions, Josephson junctions with phase-shifts can have a spatially nonuniform ground state. In recent reports, it is experimentally shown that a microwave drive can be used to measure the eigenfrequency of a junction's ground state. Such a microwave spectroscopy is based on the observation that when the frequency of the applied microwave is in the vicinity of the natural frequency of the ground state, the junction can switch to a resistive state, characterized by a non-zero junction voltage. It was conjectured that the process is analogous to the resonant phenomenon in a simple pendulum motion driven by a time periodic external force. In the case of long junctions with phase-shifts, it would be a resonance between the internal breathing mode of the ground state and the microwave field. Nonetheless, it was also reported that the microwave power needed to switch the junction into a resistive state depends on the magnitude of the eigenfrequency to be measured. Using multiple scale expansions, we show here that an infinitely long Josephson junction with phase-shifts cannot be switched to a resistive state by microwave field with frequency close to the system's eigenfrequency, provided that the applied microwave amplitude is small enough, which confirms the experimental observations. It is because higher harmonics with frequencies in the continuous spectrum are excited, in the form of continuous wave radiation. The presence of applied microwaves balances the nonlinear damping, creating a stable breather mode oscillation. We confirm our analytical results numerically.

Stability of non-time-reversible phonobreathers

Abstract: Non-time reversible phonobreathers are non-linear waves that can transport energy in coupled oscillator chains by means of a phase-torsion mechanism. In this paper, the stability properties of these structures have been considered. It has been performed an analytical study for low-coupling solutions based upon the so called {\em multibreather stability theorem} previously developed by some of the authors [Physica D {\bf 180} 235]. A numerical analysis confirms the analytical predictions and gives a detailed picture of the existence and stability properties for arbitrary frequency and coupling.

Emergence of steady and oscillatory localized structures in a phytoplankton-nutrient model

Abstract: Co-limitation of marine phytoplankton growth by light and nutrient, both of which are essential for phytoplankton, leads to complex dynamic behavior and a wide array of coherent patterns. The building blocks of this array can be considered to be deep chlorophyll maxima, or DCMs, which are structures localized in a finite depth interior to the water column. From an ecological point of view, DCMs are evocative of a balance between the inflow of light from the water surface and of nutrients from the sediment. From a (linear) bifurcational point of view, they appear through a transcritical bifurcation in which the trivial, no-plankton steady state is destabilized. This article is devoted to the analytic investigation of the weakly nonlinear dynamics of these DCM patterns, and it has two overarching themes. The first of these concerns the fate of the destabilizing stationary DCM mode beyond the center manifold regime. Exploiting the natural singularly perturbed nature of the model, we derive an explicit reduced model of asymptotically high dimension which fully captures these dynamics. Our subsequent and fully detailed study of this model - which involves a subtle asymptotic analysis necessarily transgressing the boundaries of a local center manifold reduction - establishes that a stable DCM pattern indeed appears from a transcritical bifurcation. However, we also deduce that asymptotically close to the original destabilization, the DCM looses its stability in a secondary bifurcation of Hopf type. This is in agreement with indications from numerical simulations available in the literature. Employing the same methods, we also identify a much larger DCM pattern. The development of the method underpinning this work - which, we expect, shall prove useful for a larger class of models - forms the second theme of this article.

Nonlinear Waves in Lattices: Past, Present, Future

Abstract: In the present work, we attempt a brief summary of various areas where nonlinear waves have been emerging in the phenomenology of lattice dynamical systems. These areas include nonlinear optics, atomic physics, mechanical systems, electrical lattices, nonlinear metamaterials, plasma dynamics and granular crystals. We give some of the recent developments in each one of these areas and speculate on some of the potentially interesting directions for future study.

Impulse Stability of Large Flocks: an Example

Abstract: Consider a string of N + 1 damped oscillators moving in R of which the motion of the first (called the “leader”) is independent of the others Each of the followers ‘observes’ the relative velocity and position of only its nearest neighbors Inasmuch as these are different from 0, this information is then used to determine its own acceleration Fix all parameters except the number N in such a way that the system is asymptotically stable Now as N tends tends we consider the following problem At t = 0 the leader gets kicked and starts moving with unit velocity away from the flock Due to asymptotic stability the followers will eventually fall in behind the leader and travel each at its own predetermined distance from the leader In this note we conjecture that before equilibrium ensues, the perturbations to the orbit of the last oscillator grow exponentially in N except when there is a symmetry in the interactions and the growth is then linear in N There are two cases We prove the conjecture in one case, and give a strong heuristic argument in the other

From Phase Space Representation to Amplitude Equations in a Pattern Forming Experiment

Abstract: We describe and demonstrate a method to reconstruct an amplitude equation from the nonlinear relaxation dynamics in the succession of the Rosensweig instability. A flat layer of a ferrofluid is cooled such that the liquid has a relatively high viscosity. Consequently, the dynamics of the formation of the Rosensweig pattern becomes very slow. By sudden switching of the magnetic induction, the system is pushed to an arbitrary point in the phase space spanned by the pattern amplitude and the magnetic induction. Afterwards, it is allowed to relax to its equilibrium point. From the dynamics of this relaxation, we reconstruct the underlying fully nonlinear equation of motion of the pattern amplitude. The measured nonlinear dynamics serves to select the best weakly nonlinear expansion which describes this hysteretic transition.

3D Dune Skeleton Model as a Coupled Dynamical System of 2D Cross-Sections

Abstract: To analyze theoretically the stability of the shape and the migration process of transverse dunes and barchans, we propose a {\it skeleton model} of 3D dunes described with coupled dynamics of 2D cross-sections. First, 2D cross-sections of a 3D dune parallel to the wind direction are extracted as elements of a skeleton of the 3D dune, hence, the dynamics of each and interaction between them is considered. This model simply describes the essential dynamics of 3D dunes as a system of coupled ordinary differential equations. Using the model we study the stability of the shape of 3D transversal dunes and their deformation to barchans depending on the amount of available sand in the dune field, sand flow in parallel and perpendicular to wind direction.

First and second order semi-strong interaction in reaction-diffusion systems

Abstract: Spatial scale separation often leads to sharp interfaces that can be fully localized pulses or transition layer fronts connecting different states. This paper concerns the asymptotic interaction laws of pulses and fronts in the so-called semi-strong regime of strongly differing diffusion lengths for reaction-diffusion systems in one space dimension. An asymptotic expansion and matching approach is applied in a model independent common framework. First order semi-strong interaction is introduced as a general interface interaction type. It is distinct from the semi-strong interaction studied over the past decade, which is referred to as `second order' here. Laws of motion are derived for pulses as well as fronts in abstract systems with attention to the effect of symmetries. First order interaction for pulses is shown to be gradient-like under conditions that are numerically checked for a class of equations including the Gray-Scott and Schnakenberg models.