Pulses in a Gierer Meinhardt Equation with a Slow Nonlinearity
|
|
- Derrick Campbell
- 6 years ago
- Views:
Transcription
1 SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 1, No. 1, pp c 013 Society or Industrial and Applied Mathematics Pulses in a Gierer Meinhardt Equation with a Slow Nonlinearity Frits Veerman and Arjen Doelman Abstract. In this paper, we study in detail the existence and stability o localized pulses in a Gierer Meinhardt equation with an additional slow nonlinearity. This system is an explicit example o a general class o singularly perturbed, two component reaction-diusion equations that goes signiicantly beyond well-studied model systems such as Gray Scott and Gierer Meinhardt. We investigate the existence o these pulses using the methods o geometric singular perturbation theory. The additional nonlinearity has a proound impact on both the stability analysis o the pulse compared to Gray Scott/Gierer Meinhardt-type models a distinct extension o the Evans unction approach has to be developed and the stability properties o the pulse: several (de)stabilization mechanisms turn out to be possible. Moreover, it is shown by numerical simulations that, unlike the Gray Scott/Gierer Meinhardt-type models, the pulse solutions o the model exhibit a rich and complex behavior near the Hop biurcations. Key words. dynamical systems, pattern ormation, singular perturbations AMS subject classiications. 35K57, 35B5, 35B35, 34A6, 34C37 DOI / Introduction. The study o localized pulses in a two-component system o singularly perturbed reaction-diusion equations has been a very active ield o research since the nineties o the previous century. In its most general orm, a system that may exhibit such a pulse reads in one, unbounded, spatial dimension as { (1.1a) Ut = U xx + F (U, V ), (1.1b) V t = ε V xx + G(U, V ), with U, V : R R + R, and0<ε 1asymptotically small. The nonlinear reaction terms F, G : R R areassumedtosatisyf (Ū, V )=G(Ū, V ) = 0 or certain (Ū, V )such that the trivial background state (U, V ) (Ū, V ) is spectrally stable. However, research on pulses in equations o the type (1.1) has been restricted mostly to model equations. In particular two o these models have played a central role in the development o the theory: the (irreversible) Gray Scott (GS) equation or a class o autocatalytic reactions [17] that became the center o research attention by the intriguing observations in [5, 8] and the Gierer Meinhardt (GM) equation [16] or (biological) morphogeneses or which the existence problem had already been considered in the mathematical literature or a somewhat longer time [31]. For both the GS and the GM models, quite precise insight has been obtained in Received by the editors May 4, 01; accepted or publication (in revised orm) by K. Promislow September 5, 01; published electronically January 17, 013. This research was unded by NWO through the VIDI-grant Formation o singularities in natural systems, awarded to V. Rottschäer, Universiteit Leiden, The Netherlands. Universiteit Leiden, Mathematical Institute, 333 CA Leiden, The Netherlands (rits@math.leidenuniv.nl, doelman@math.leidenuniv.nl). 8
2 PULSES IN A SLOWLY NONLINEAR GM EQUATION 9 the existence, stability, and dynamics o localized (multi)pulses, also in more than one spatial dimension although one certainly cannot claim that the models are ully understood; see [4, 5, 6, 11, 9, 18, 1,, 3, 6, 7, 3, 34] and the reerences therein or the literature on one spatial dimension. The ast V -component o a localized (multi)pulse solution o a singularly perturbed model (1.1) is asymptotically localized: it decays exponentially to the V -component V o the background state on a spatial scale that is asymptotically shorter than the spatial scale associated with the slow U-component. As a consequence, the two-component (U, V )-low generated by (1.1) is governed by a scalar equation in the slow component U: (1.) U t = U xx + F (U, V ) except or the asymptotically small spatial regions in which the V -component is not exponentially close to V. Clearly, this is in general a nonlinear equation. However, or the GS and GM models, this slow reduced scalar equation is linear: (1.3) (GS) U t = U xx + A(1 U), A > 0 parameter, (Ū, V )=(1, 0), (GM) U t = U xx αu, α > 0 parameter, (Ū, V )=(0, 0). In act, as ar as we are aware, this the act that the counterpart o (1.) is linear is the case or all singularly perturbed two-component reaction-diusion equations with exponentially localized pulse solutions considered in the literature (including the Schnakenberg model [33, 30]). There are a number o papers in the literature in which more general classes o equations than the GS or GM models are considered; see [5, 10]. In these papers the background state (Ū, V )istranslatedto(0, 0) so that F (Ū, V )=F(0, 0) = 0 in (1.1). Moreover, the nonlinear part o F (U, V ) is assumed to be separable; i.e., F (U, V ) is written as αu + F 1 (U)F (V ). Thereore F ( V )=F (0) = 0, and these more general systems also reduce to linear slow scalar equations like (1.3) outside the asymptotically small regions where V is not close to V. In this paper, and in its companion paper [13], we consider the potential impact o the nonlinearity o F (U, V ) as a unction o U in comparison with the literature on slowly linear model systems such as GS and GM. Here, we consider a very explicit model problem, a GM equation with a slow nonlinearity (see (1.7) below), in ull analytical detail; in [13] we consider the existence and stability o pulses in a general setting, i.e., as solutions o (1.1). We reer the reader to Remark 1.1 or a more speciic motivation o our choice to study equations with slow nonlinearities. In the standard orm (1.1), the classical GM equation [16] isgivenby (1.4a) U t = U xx αu + σv, (1.4b) V t = ε V xx V + V U, in which α>0 is the main biurcation parameter and σ>0is most oten scaled to 1. The pulse-type solutions o (1.4) have an amplitude o O ( ) 1 ε [5, 18]. Thereore, we scale U and V and subsequently x and ε, (1.5) U U ε, V V ε, x εx, ε ε,
3 30 FRITS VEERMAN AND ARJEN DOELMAN to bring (1.4) in its normal orm [5] (1.6a) (1.6b) ε U t = U xx ε αu + σv, V t = ε V xx V + V U. In this paper, we study a slowly nonlinearized version o (1.6) that is obtained rom (1.6) by adding a very simple nonlinear term to its slow U-equation (1.6a): (1.7a) ( ε U t = U xx ε αu γu d) + σv, (1.7b) V t = ε V xx V + V U, with new parameters γ 0, d > 1. Moreover, we now allow σ R\{0}. Systems incorporating such a slow nonlinearity were already encountered in [4] (although no pulse-type solutions were considered in this paper). This equation indeed reduces to a nonlinear slow reduced scalar U-equation away rom the regions in which V is not exponentially close to V =0: (1.8) U t = U χχ αu + γu d, in which χ = εx is a superslow spatial coordinate; see section. Note that scaling back the additional slowly nonlinear term γu d through (1.5) introduces an O ( ε d 1), i.e., an asymptotically small, additional term to the GM equation in its classical orm (1.4). We will see in the upcoming analysis that this term has a signiicant impact on the dynamics generated by (1.4). Thus, in a way, our work can also be interpreted as a study o the vulnerability o the classical GM model (1.4) to asymptotically small slowly nonlinear changes to the model. In recent years, the analysis o localized pulses in one-dimensional singularly perturbed reaction-diusion equations has been ocused mostly on pulse dynamics and interactions; see [4, 10, 9, 3] and the reerences therein. However like the work on multipulse patterns [8, 18, 1,, 3, 33] this analysis is based on undamental insights on the existence and stability o stationary, solitary pulses [5, 6, 11, 31, 34]. On the unbounded domain, i.e., or x R, these pulses correspond to homoclinic solutions o the our-dimensional spatial dynamical system reduction o the partial dierential equation. Here, we restrict our analysis to the existence and stability o homoclinic stationary pulse solutions (U h (x),v h (x)) to (1.7) that are biasymptotic to the background state (0, 0), i.e., lim x ± (U h (x),v h (x)) = (0, 0). Especially the issue o stability requires a signiicant extension o the methods developed in the literature or slowly linear GS/GM-type models. The present results orm the oundation or a subsequent analysis o the multipulse patterns (see Remark 1.1) and pulse interactions. Moreover, already at the level o these most basic pulse solutions, we encounter novel phenomena in the dynamics generated by (1.7) that have not yet been observed in the literature on slowly linear models. The existence problem (see section ) can be studied directly along the lines developed in [5] or slowly linear normal orm models o GM type with a separable nonlinearity. Our main result on the existence o homoclinic pulses (U h (x),v h (x)), Theorem.1, can be established
4 PULSES IN A SLOWLY NONLINEAR GM EQUATION 31 U h V h (a) x Figure 1. The stationary homoclinic pulse (U h (x),v h (x)) o (1.7); (a) σ>0, (b) σ<0. by a direct application o the methods o geometric singular perturbation theory [14, 15]. In other words, at the existence level the slow nonlinearity in (1.7) doesnotrequirethe development o novel theory. However, it is established by Theorem.1 that (1.7) exhibits homoclinic pulse patterns that dier signiicantly rom those ound in slowly linear GS/GMtype models. Unlike linear slow reductions such as (1.3), the planar stationary problem associated with reduction (1.8) has orbits homoclinic to its saddle point (that corresponds to the background state o (1.7)). As a consequence, unlike the classical GM model (1.4), system (1.7) has homoclinic pulse solutions (U h (x),v h (x)) or σ<0. At leading order in ε, theslow U-component U h (x) ollows a large part o the homoclinic orbit o (1.8) sothatorσ<0the slow component o the solitary homoclinic 1-pulse solution has the leading order structure o two combined slow scalar pulses; see Figure 1(b). The spectral stability o (U h (x),v h (x)) is studied in section 3 by the Evans unction D(λ) associated with the linearized stability problem, ollowing the ideas developed in [5, 6]. As is to be expected rom the general theory [3], D(λ) can be decomposed into a slow and a ast component, and all nontrivial eigenvalues are determined by the slow component. In [5, 6], i.e., or the GM and GS models, the zeroes o this slow component are determined analytically by the NLEP method. The linearity o the slow scalar reduction (1.3) plays a central role in this approach as it does in all analytical studies o the spectral stability o pulses in GS/GMtype models (see [4, 18, 1,, 3, 3, 34] and the reerences therein). More explicitly, the act that the spectral stability problem is exponentially close to a constant coeicient eigenvalue problem outside the asymptotically small regions in which V is not close to V is a crucial ingredient o the stability analysis o GS/GM-type models. Due to the nonlinearity in the slow scalar reduction (1.8) this is not the case or (1.7): away rom the ast V -pulse, the linear operator associated with the stability problem still has coeicients that depend explicitly (and slowly) on x (on χ; see(1.8)). Its solution space is thereore not governed by simple, pure exponentials (as or GS/GM-type models). The key to the NLEP approach as developed in [5, 6] is constructing a set o basis unctions or the linear operator/system associated with the stability o the pulse or which the Evans unction D(λ) the determinant o this set can be evaluated, or, better, approximated, explicitly. In this paper and its companion paper [13], we show that the NLEP approach can be based on a set o basis unctions that is determined by the slowly varying problem outside the U h V h (b) x
5 3 FRITS VEERMAN AND ARJEN DOELMAN Im Λ Im Λ 15 Im Λ Re Λ (a) (b) Re Λ Re Λ Figure. The orbits through C o the critical eigenvalue λ associated with the spectral stability o (U h (x),v h (x)) as a unction o increasing α (γ =, σ =1) to leading order in ε. (a) d =< 3: The same scenario as in the GS and the GM models [5, 6]. Two real positive eigenvalues merge and become a pair o complex conjugate eigenvalues that travels through the imaginary axis: the pulse is stabilized by a Hop biurcation at a critical value o α. (b) d =5> 3: A signiicantly dierent scenario. The eigenvalues initially display the same behavior as in the case d<3: the pulse is again stabilized by a Hop biurcation. However, or α increasing urther, the orbits sharply turn around and ollow the imaginary axis closely in the negative direction; see (c), a zoom o (b). Eventually, the orbits branch o, head back to the imaginary axis, and again cross the imaginary axis at a second critical Hop value o α. Finally, the pair meets again at the positive real axis and splits up into two positive real valued eigenvalues. ast V -pulse region, in such a way that it is still possible to determine an analytical approximation or the zeroes o D(λ). Here, a central role is played by the χ-dependent Sturm Liouville problem associated with the linearization o (1.8) about its (stationary) homoclinic orbit, deined on a hal-line. This problem has a two-dimensional set o slowly varying solutions. We show that these solutions can take over the role o the slow exponentials coming rom the (slow) stability problem about the trivial state U = 0 o the linear constant coeicient GS/GM-type reductions (1.3). In the present paper, these solutions can be expressed in terms o Legendre unctions, due to the special/simple nature o the nonlinearity in (1.8). In the general setting o [13] the construction o the Evans unctions cannot be this explicit. The main novel analytical result o this paper is given by Theorem 3.10, in which indeed an explicit expression is given or the zeroes o D(λ), which is a generalization o the corresponding slowly linear results in [5, 6]. In section 4, we analyze and interpret the expression obtained in Theorem One o our irst, and quite straightorward, results is Corollary 4.3: the σ<0 double hump pulses o Figure 1(b) cannot be stable. The σ>0 pulses o Figure 1(a), however, can very well be stable. In Figure, a graphical description is given o our two main stability results, Theorems 4.6 and 4.7. The stability o the pulse (U h (x),v h (x)) depends strongly on the character o the slow nonlinearity in (1.7). As long as the exponent o the nonlinearity d is smaller than 3, the stability scenario is exactly like that o the slowly linear GS/GM-type models: (U h (x),v h (x)) stabilizes by a Hop biurcation or increasing α even the shape o the orbit o the critical eigenvalues λ(α) through C is very similar to its counterparts in [5, 6]. However, this orbit changes drastically when d becomes larger than 3; see Figure : or d > 3 there is a second Hop biurcation (as unction o α) that destabilizes (U h (x),v h (x)). Dierent rom the results (c) 1 1
6 PULSES IN A SLOWLY NONLINEAR GM EQUATION 33 on GS/GM-type models, or d>3, there is only a bounded α-region or which (U h (x),v h (x)) can be stable. In Theorem 4.8 this is established rigorously or d>3 large enough. Finally, in section 5, we present some simulations o (1.7). We have not attempted to perorm a systematic (numerical biurcation) analysis o the dynamics o (1.7). Apart rom checking (and conirming) the outcome o our asymptotic stability analysis, our goal has been to obtain an indication o whether or not the slow nonlinearity o (1.7) generates behavior that is not known rom the (vast) literature on GS/GM-type models. We are not aware o any examples in the literature on GS/GM-type models o stable nonmoving solitary pulses that are not completely stationary. A priori, one would expect that i a pulse is destabilized (or instance, by decreasing α in the GM model (1.4)), it may biurcate into a stable standing pulse with a periodically varying amplitude. However, this requires a supercritical Hop biurcation, and all Hop biurcations o stationary pulses in GS/GM-type models reported on in the literature seem to be subcritical: as α decreases through its critical Hop biurcation value, the standing pulse starts to oscillate up and down, but the amplitude o this oscillation grows, and ater a certain time the pulse is extinguished; see, or instance, Figure 11(a) in section 5. It should be noted that this statement is based on numerical observations; the nature o the Hop biurcation o solitary, standing pulses in GS/GM-type models has not been analyzed in the literature (or instance, by a center maniold reduction). Moreover, it should also be remarked that, or instance, the GS model does exhibit periodic and even chaotic pulse dynamics; see, or instance, [4, 7]. However, this richer type o behavior occurs only in the context o pulse interactions; it is governed by the interactions between traveling pulses, and/or between pulses and the boundary o the domain. We have not considered this type o dynamics here, as we have completely ocused on the behavior o standing, solitary spatially homoclinic pulses. Nevertheless, we have observed very rich dynamics, much richer than those exhibited by linear GS/GM-type models. In section 5 examples are given o periodically oscillating pulses, i.e., standing pulses with an amplitude that varies periodically in time; quasi-periodically oscillating pulses (the amplitude o the pulse oscillation is modulated periodically); and oscillating pulses o which the amplitude is modulated in an even more complex ashion. A simulation o such a chaotically oscillating pulse is shown in Figure 3. It has not been investigated whether the pulse dynamics o Figure 3 are chaotic or, or instance, are quasi-periodic with three or more independent requencies. In other words, we have not studied the details o the associated biurcation scenario and have not computed any measure by which the (possible) chaotic nature o the pulse dynamics can be quantiied. The analytic core o this paper, the analysis o the spectrum associated with the stability o (U h (x),v h (x)), can serve as an ideal starting point or a center maniold analysis o the nature o the Hop (and subsequent) biurcations or pulses and/or multipulses occurring in this model (and/or generalizations o (1.7)). This is expected to give analytic insight into (the possible route leading to) the complex/chaotic behavior observed in Figure 3. This will be the subject o uture research. Remark 1.1. Our research is strongly motivated by recent indings on the character o the destabilization o spatially periodic multipulse patterns with long wavelength L. In [8] itis established or GM-type models that these patterns can be destabilized only by two distinct types o Hop biurcations as L, one in which the linearly growing mode also has
7 34 FRITS VEERMAN AND ARJEN DOELMAN amplitude t (a) amplitude Figure 3. The dynamics o the maximum o the U-pulse as a unction o time in a simulation o (1.7) with γ =, σ =1, ε =0.00, d =5,andα =90.6 or x [ 5000, 5000] with homogeneous Neumann boundary conditions; (b) zooms in on a small part (in time) o (a). The position o the maximum does not vary in time. The value o α is close to the second Hop biurcation at which (U h (x),v h (x)) destabilizes; see Figure. wavelength L (the most commonly encountered destabilization in the literature) and another in which this mode has wavelength L. Moreover, these destabilizations alternate countably many times as L. This is called the Hop dance in [8]. This Hop dance also occurs in the GS model, as indicated by the AUTO-simulations in [8]. The GM analysis in [8]showsthat this dance is completely driven by the exponential expression E(L) =e L α+λ h associated with the slow reduced eigenvalue equation u xx αu = λ h u originating rom (1.3), in which λ h C is the (complex) eigenvalue o the homoclinic (L ) limit pattern. The rotation o E(L) C as L is the mechanism underpinning the Hop dance. From a generic point o view, it is not at all clear why this linear Hop dance should take place (this is even more obvious or the subsequent belly dance [8]). Hence, to really understand the subtleties involved in the destabilization o long wavelength spatially periodic patterns, one needs to go beyond slowly linear models or which the associated slow reduced eigenvalue problems are not governed by expressions such as E(L). In other words, one needs to study systems o the type (1.1) withf (U, V ) not linear as a unction o U.. Pulse construction. Our goal is to construct a stationary pulse solution which is homoclinic to the trivial background state (U, V )=(0, 0). To achieve this goal, we use the singularly perturbed nature o the system. The spatial dynamics o the stationary pulse are given by the our-dimensional system (.1a) (.1b) (.1c) (.1d) ( p x = σv + ε αu γu d), εv x = q, u x = p, εq x = v v u. (b) t 40
8 PULSES IN A SLOWLY NONLINEAR GM EQUATION 35 Along the lines o Fenichel theory, we can perorm a slow-ast decomposition in the spatial variable x: recognizing system (.1) astheslow system, we can deine the ast variable ξ = x ε to obtain the associated ast system (.a) u ξ = εp, ( (.b) p ξ = εσv + ε 3 αu γu d), (.c) v ξ = q, (.d) q ξ = v v u. The trivial background state is in these systems represented by the origin (u, p, v, q) = (0, 0, 0, 0). While the vector ield which generates the low o the system is not deined at the origin due to the singular v v u term in the v-equation, the ratio u will be well-deined or the constructed pulse..1. Geometric analysis. When ε 0, the slow and ast systems (.1) and(.) reduce to the reduced slow system (.3a) (.3b) and the reduced ast system (.4a) u xx = σv, q = v v u =0 u ξ = p ξ =0, (.4b) v ξξ = v v u. We see that in this limit, the slow and ast dynamics decouple completely. We deine M 0 = {(u, p, v, q) u > 0, v = q = 0} as the two-dimensional normally hyperbolic invariant maniold that consists o hyperbolic equilibria o the reduced ast system (.4); it has three-dimensional stable and unstable maniolds W s,u (M 0 ) which are the unions o the -parameter amilies o one-dimensional stable and unstable maniolds (ibers) at the saddle points (u 0,p 0, 0, 0) M 0. The reduced ast dynamics (.4) allow a -parameter amily o homoclinic solutions v 0,h : (.5) v h,0 (ξ; u 0,p 0 )= 3 u 0 sech ( 1 ξ The union over this amily as a bundle over M 0 orms the intersection W s (M 0 ) W u (M 0 ); see Figure 4(a). Fenichel persistence theory [14, 15, 19, 0] states that, or ε suiciently small, the ull system (.) has a locally invariant slow maniold M ε which is O(ε) closetom 0.SinceM 0 is also invariant under the nonreduced (ast) low o (.), we have already ound M ε = M 0.Moreover, Fenichel theory states the existence o three-dimensional stable and unstable maniolds W s,u (M ε )whichareo(ε) close to their unperturbed counterparts W s,u (M 0 ). The intersection W s (M ε ) W u (M ε ) exists and is transversal and thereore determines a two-dimensional ).
9 36 FRITS VEERMAN AND ARJEN DOELMAN (a) (b) Figure 4. Transversal intersection o the stable and unstable maniolds. (a) The amily o homoclinic orbits vh,0 (ξ; u0, p0 ), viewed as a bundle over M0. Both the u- and p-directions are along the vertical axis; M0 is indicated in blue. (b) For the perturbed system (ε > 0), W s (Mε ) and W u (Mε ) intersect transversally: γh (indicated in red) represents W s (Mε ) W u (Mε ), a one-parameter amily o orbits homoclinic to Mε recall that dim (Mε ) =, dim (W s,u (Mε )) = 3 so that dim (γh ) =. maniold. The existence and transversality are based on a Melnikov-type calculation in [5], which can be applied directly to system (.). Since the original model equations (1.7) are invariant under reflection in the spatial variable x x, this reflection is in the our-dimensional system (.1) equivalent to the momentum reflection (p, q) ( p, q). Because the coordinate reflection ξ ξ maps W s (Mε ) to W u (Mε ) and vice versa, it ollows that the intersection o these two maniolds is symmetric in the invariant subset o the momentum reflection, the twodimensional hyperplane {(u, p, v, q) p = q = 0}. The transversality o this hyperplane to Mε excludes the possibility that it has the intersection W s (Mε ) W u (Mε ) as a subset, rom which we can conclude that W s (Mε ) W u (Mε ) intersects the hyperplane {(u, p, v, q) p = q = 0} transversally. This determines a 1-parameter amily o orbits biasymptotic to Mε. Since both W s (Mε ) and W u (Mε ) are O(ε) close to W s,u (M0 ) where the -parameter amily o homoclinic orbits was parametrized by u0 and p0 (see (.5)), it is convenient to use u0 to parametrize the 1-parameter amily o orbits biasymptotic to Mε determined by W s (Mε ) W u (Mε ). For a sketch o the situation, see Figure 4(b). The next step is to use this structure to construct an orbit homoclinic to (0, 0, 0, 0) in the ull, perturbed system (.1) (.). For that purpose, it is necessary to consider the dynamics
10 PULSES IN A SLOWLY NONLINEAR GM EQUATION 37 on M ε. The low on M ε can be determined by substituting v = q =0in(.1) and yields (.6) u xx = ε ( αu γu d). Introducing a superslow coordinate χ = εx, this can be written as (.7) u χχ = αu γu d. This equation allows a solution (bi)asymptotic to the trivial background state: since γ>0, it is homoclinic to (0, 0, 0, 0) M ε and explicitly given by [ α(d +1) (.8) u h,0 (χ) = γ ( )] 1 1 sech (d 1) d 1 αχ. The superslow dynamics on M ε allow us to get a grip on picking exactly that orbit biasymptotic to M ε rom the intersection W s (M ε ) W u (M ε ) which is also homoclinic to (0, 0, 0, 0) M ε, that is, which is (mostly) asymptotically close to u h,0 M ε. This orbit will make a ast excursion through the V -ield, since this is where the ast dynamics take place(see(.), (.4)). Since our goal is to construct a symmetric pulse, we can choose an interval symmetric around the origin in which the ast jump occurs. The interval needs to be asymptotically small with respect to the slow variable x but asymptotically large with respect to the ast variable ξ: to be asymptotically close to M ε,thev-component o the pulse needs to be exponentially small. A standard choice [5] or this ast spatial region is { (.9) I = ξ R ξ < 1 }. ε Indeed, x 1andξ 1on I. For a sketch o the orbit, see Figure 5. Now, we deine the take-o and touchdown sets T o,d M ε to be the collection o basepoints o all Fenichel ibers in W u (M ε )(resp.,w s (M ε )) that have points in the transverse intersection W s (M ε ) W u (M ε ). Detailed inormation on T o,d can be obtained by studying the ast system (.) onm ε. First, we observe that p ξ = O(ε 3 )onm ε so that the p-coordinate on M ε remains constant to leading order during the ast excursion through the V -ield. Thereore, the change in the p-coordinate o the pulse is completely determined by its accumulated change during its excursion through the ast ield, and is given by Δ ξ p = p ξ dξ = εσv + O(ε 3 )dξ = εσv h,0 (ξ; u 0,p 0 ) dξ + O(ε ) I I (.10) = 6 εσu 0 + O(ε ), where we have used (.) and(.5). Moreover, since u ξ = εp and p = O(ε) oni,weseethat Δ ξ u = O(ε ). This means that during the jump through the ast ield, the u-coordinate o the pulse does not change to leading order. Since W s (M ε ) W u (M ε ) intersects the hyperplane {(u, p, v, q) p = q =0} transversally, we can deine the take-o and touchdown sets as curves (.11) T o = { (u, p, 0, 0) M ε p =3εσu } and T d = { (u, p, 0, 0) M ε p = 3 εσu }
11 38 FRITS VEERMAN AND ARJEN DOELMAN Figure 5. An asymptotic construction o the orbit γ h (ξ) o Theorem.1 drawn in three dimensions. The p- and q-directions are combined, since there is no direct interaction between them. The blue surace represents the persistent slow maniold M ε, while the ast dynamics take place on the red surace, which is spanned by the v- andq-directions. The slow homoclinic orbit u h,0 (χ) is drawn in blue, while the ast homoclinic orbit v h,0 (ξ; u, 0) is drawn in red. The jump through the ast ield projected on M ε is indicated by the purple line. p u h,0 T o T d (a) For σ>0 the jump is downwards; the resulting pulse is shown in Figure 1(a). u * u p u h,0 T d To (b) For σ<0 the jump is upwards; the resulting pulse is shown in Figure 1(b). Figure 6. The homoclinic orbit u h,0 (χ) isdrawninblueinthe(u, p)-plane. The take-o and touchdown curves T o = { (u, p) p =3εσu } and T d = { (u, p) p = 3 εσu } are drawn in green. The jump through the ast ield at u = u is indicated by the dashed purple line. at leading order. Note that i σ changes sign, the take-o and touchdown curves are interchanged: or σ>0the take-o curve has positive p-values, while or σ<0 the take-o curve has negative p-values. This also means the direction o the ast jump is reversed when σ changes sign; see (.10) andfigure6(a). u * u
12 PULSES IN A SLOWLY NONLINEAR GM EQUATION 39 An orbit o the system (.1) (.) is homoclinic to (0, 0, 0, 0) i its Fenichel iber basepoints in T o,d intersect the superslow homoclinic orbit u h,0 M ε ;seefigure6(a). This intersection can be determined by integrating (.6) once, (.1) 1 p = ε ( 1 αu γ d+1 ud+1), and substituting p = ±3 εσu rom (.11) toobtain γ (.13) d +1 ud 1 = α 9 σ u, which or α, γ, σ > 0andd>1 always has a unique real positive solution, denoted by u. Furthermore, we deine χ as the (unique) positive χ-value or which u h,0 (χ )=u,theucoordinate o the intersection. When σ<0, we obtain a slightly dierent pulse, since part o the slow homoclinic orbit u h,0 is covered twice; see Figure 6(b). This has its consequences or the ormulation o our main existence result. Theorem.1. Let ε>0 be suiciently small. Then, or all values o the parameters α>0, γ > 0, σ > 0, andd>1, there exists a unique orbit γ h (ξ) =(u h (ξ),p h (ξ),v h (ξ),q h (ξ)) as a solution o system (.) which is homoclinic to (0, 0, 0, 0) and lies in the intersection W s (M ε ) W u (M ε ). Moreover, (.14) v h (ξ) v h,0 (ξ; u, 0) = O(ε), q h (ξ) d dξ v h,0(ξ; u, 0) = O(ε) or all ξ R and (.15) u h (χ) u h,0 (χ sgn(σ) χ ) = O(ε), p h (χ) ε d dχ u h,0(χ sgn(σ) χ ) = O(ε) or all χ<0, while (.16) u h (χ) u h,0 (χ + sgn(σ) χ ) = O(ε), p h (χ) ε d dχ u h,0(χ + sgn(σ) χ ) = O(ε) or all χ>0. The orbit γ h corresponds to a homoclinic pulse solution (U h,v h ) o system (1.7). Proo. The missing details in the above geometric construction, especially in the precise estimates o (.14), (.15), and (.16), can be obtained in a manner identical to the corresponding result on slowly linear systems in [5]. 3. Pulse stability: Analysis. The linear stability o the stationary pulse solution (U h,v h ) o (1.7) ound in the previous section is determined by adding a perturbation o the orm (ū(x), v(x)) e λt and linearizing equation (1.7) around the stationary solution, obtaining in the ast variable ξ (3.1a) (3.1b) (3.1c) (3.1d) ū ξ = ε p, ( p ξ = εσv h (ξ) v + ε 3 α + λ γdu h (ξ) d 1) ū, v ξ = q, ( q ξ = 1+λ V ) h(ξ) v + V h(ξ) U h (ξ) U h (ξ) ū.
13 40 FRITS VEERMAN AND ARJEN DOELMAN We write the ast system (3.1) in vector orm d (3.) dξ φ = A(ξ; λ, ε)φ, where φ(ξ) =(ū(ξ), p(ξ), v(ξ), q(ξ)) T and 0 ε 0 0 ε 3 ( α + λ γdu h (ξ) d 1) 0 εσv h (ξ) 0 (3.3) A(ξ; λ, ε) = V h (ξ) 0 1+λ V h(ξ) U h (ξ) U h (ξ) 0 Since the V -component o the stationary pulse decays much aster than its U-component, the ratio V h U h is well-deined and converges to zero as ξ ±. This results in the constant coeicient matrix 0 ε 0 0 (3.4) A (λ, ε) = lim A(ξ; λ, ε) = ε 3 (α + λ) ξ λ 0 which has eigenvalues (3.5) ± Λ = ± 1+λ and ± ε Λ s = ±ε α + λ and associated eigenvectors (3.6) E,± = ( 0, 0, 1, ± T ( 1+λ) and E s,± = 1, ±ε T α + λ, 0, 0). The essential spectrum o the linear eigenvalue problem (3.1) thereore is (3.7) σ ess = {λ R λ max( α, 1)} ; see [9]. Since α>0, we can conclude that the stability o the pulse (U h,v h ) is determined by its discrete spectrum The Evans unction and its decomposition. The Evans unction, which is complex analytic outside the essential spectrum (see [9, 3] and the reerences therein) associated with system (3.1) can be deined by (3.8) D(λ, ε) =det[φ i (ξ; λ, ε)], where the unctions φ i, i =1,, 3, 4, satisy boundary conditions at ± (see below) and span the solution space o (3.1). The eigenvalues o (3.) outside σ ess coincide with the roots o D(λ, ε), including multiplicities. Deinition 3.1. A statement o the orm (x) cg(x) as x is true whenever the 1 limit lim x g(x) (x) =c exists and is well-deined.
14 PULSES IN A SLOWLY NONLINEAR GM EQUATION 41 Lemma 3.. For all λ C \ σ ess, there are solutions φ,l/r (ξ; λ, ε) and φ s,l/r (ξ; λ, ε) to (3.1) such that the set { φ,l/r (ξ; λ, ε),φ s,l/r (ξ; λ, ε) } spans the solution space o (3.1) and (3.9a) (3.9b) (3.9c) (3.9d) φ,l (ξ; λ, ε) E,+ e Λ ξ φ,r (ξ; λ, ε) E, e Λ ξ φ s,l (ξ; λ, ε) E s,+ e ε Λ sξ φ s,r (ξ; λ, ε) E s, e ε Λ sξ as as as as ξ, ξ, ξ, ξ. Moreover, there exist analytic transmission unctions t,+ (λ, ε) and t s,+ (λ, ε) such that (3.10a) (3.10b) φ,l (ξ; λ, ε) t,+ (λ, ε) E,+ e Λ ξ φ s,l (ξ; λ, ε) t s,+ (λ, ε) E s,+ e ε Λ sξ as as ξ, ξ, where t s,+ (λ, ε) is deined only i t,+ (λ, ε) 0. These choices, when possible, determine φ,l/r and φ s,l uniquely. Proo. Although the linearized system (3.) is not identical to its counterpart in [5], exactly thesameargumentsasin[5] can be applied here. Thereore, we reer the reader to [5] orthe details o the proo. The Evans unction can be determined by taking the limit ξ o the determinant o the unctions deined in Lemma 3., since the Evans unction itsel does not depend on ξ since the trace o A(ξ; λ, ε) vanishes (Abel s theorem). This yields (see [5]) (3.11) D(λ, ε) =4εt,+ (λ, ε) t s,+ (λ, ε) 1+λ α + λ. Corollary 3.3. The set o eigenvalues o (3.) is contained in the union o the sets o roots o t,+ (λ, ε) and t s,+ (λ, ε). Note that, due to the act that t s,+ (λ, ε) is deined only when t,+ (λ, ε) 0,theEvans unction D(λ, ε) does not necessarily vanish when t,+ (λ, ε) = 0. This is called the resolution to the NLEP paradox in [5, 6]. The roots o t,+ will be discussed later, in section The slow solution φ s,l outside I. To obtain more inormation about the roots o t s,+ (λ, ε), it is necessary to determine the leading order behavior o φ s,l (ξ; λ, ε) inthe dierent coordinate regimes. From Lemma 3. we know that φ s,l is slowly growing in ξ, since its leading order behavior or both ξ ± is determined by the exponential growth actor ε Λ s = O(ε ). However, the dynamics governing φ s,l dier signiicantly inside and outside the ast spatial region I. Based on our knowledge o the homoclinic solution stated in Theorem.1, we can iner the orm o the matrix A(ξ; λ, ε) both inside and outside I : 0 ε 0 0 ε 3 ( α + λ γdu d 1 ) 0 εσv h,0 (ξ; u, 0) 0 (3.1) A (ξ; λ, ε) = λ v h,0(ξ;u,0) u 0 v h,0 (ξ;u,0) u
15 4 FRITS VEERMAN AND ARJEN DOELMAN to leading order or ξ I and 0 ε 0 0 (3.13) A s (ξ; λ, ε) = ε 3 ( α + λ γdu h,0 ( ε ξ +sgn(σ) χ ) d 1) λ 0 or ξ/ I to leading order. Note that it is the act that this intermediate slow matrix exists, or, better, that it is not identical to A (λ, ε)(3.4), that distinguishes the slowly nonlinear GM problem rom slowly linear problems such as the classical GM or GS systems. Note also that a intermediate matrix as A s (ξ; λ, ε) was already encountered in [7], in the study o a system with nonexponential (algebraic) decay. Lemma 3.4. Consider the system d (3.14) dξ ψ = A s(ξ; λ, ε)ψ with A s (ξ; λ, ε) as given in (3.13). There exist solutions ψ,± (ξ; λ, ε) and ψ s,± (ξ; λ, ε) which span the solution space o (3.14) or ξ< 1 ε and (3.15a) (3.15b) ψ,+ (ξ; λ, ε) E,+ e Λ ξ, ψ s,+ (ξ; λ, ε) E s,+ e ε Λ sξ, ψ, (ξ; λ, ε) E, e Λ ξ, ψ s, (ξ; λ, ε) E s, e ε Λ sξ as ξ. Proo. The same arguments as in the proo o Lemma 3. apply, since lim ξ A s (ξ; λ, ε) = A (λ, ε). Since A(ξ; λ, ε) is to leading order equal to A s (ξ; λ, ε) orξ< 1 ε and both φ s,l and ψ s,+ E s,+ e ε Λ sξ as ξ, combining Lemmas 3. and 3.4 yields the ollowing corollary. Corollary 3.5. For ξ< 1 ε,wecanwrite φ s,l (ξ; λ, ε) =ψ s,+ (ξ; λ, ε) to leading order. The slow evolution o the ū-component o ψ s,± can be written, again using χ = ε ξ,as ( (3.16) ū χχ α + λ γdu h ( χ +sgn(σ) χ ) d 1) ū =0. We can introduce the coordinate transormation z = 1 d dχ u h,0(χ sgn(σ) χ ) α u h,0 (χ sgn(σ) χ ) = 1 d α dχ log 1 u h,0 (χ sgn(σ) χ ) ( ) 1 (3.17) =tanh (d 1) α (χ sgn(σ) χ )
16 PULSES IN A SLOWLY NONLINEAR GM EQUATION 43 (by (.8)) or the region χ<0toobtain (3.18) (1 z )ū zz z ū z + where ) (ν(ν +1) μ 1 z ū =0, (3.19a) ν = d +1 d 1, μ =+ (3.19b) 1+ λ d 1 α, where we have chosen the branch cut associated with σ ess such that Re μ>0; note that ν>1. Equation (3.18) is the Legendre dierential equation: its solutions are the associated Legendre unctions Pν μ (z) andq μ ν(z) (see[, 1]). Given the symmetry z z o the equation, we choose the basis o the solution space to be Pν μ (±z). The limit χ corresponds to the limit z 1. Taking into account the normalization o ψ s,+ rom Lemma 3.4, the correct expression or the ū-component o ψ s,+ is (3.0) ū(χ) =Γ(1+μ) e Λssgn(σ) χ Pν μ ( z(χ)) such that (3.1) lim χ 0 ū(χ) =Γ(1+μ) e Λssgn(σ) χ P μ ν (z ), where we deine (3.) z =sgn(σ) tanh ( 1 (d 1) αχ ). We can express z in terms o u using (.7): integrating once yields (3.3) u χ = αu γ d +1 ud+1, so, by (3.17), (3.4) z = 1 α u χ u =1 γ α(d +1) ud 1, and hence (3.5) α ( 1 z ) γ = d +1 ud 1 = α 9σ u by (.13); rom this, we conclude that (3.6) z = 3 σ α u. Note that z inherits the sign o σ, sinceχ is chosen to be positive; see section.1.
17 44 FRITS VEERMAN AND ARJEN DOELMAN Lemma 3.6. Let ū s (ξ; λ, ε) be the ū-component o φ s,l (ξ; λ, ε) as deined in Lemma 3.. Then (3.7) ū s (ξ) =Γ(1+μ) e Λssgn(σ) χ P μ ν (z )+O(ε ε) or ξ I. Moreover, there are two transmission unctions t s,+ (λ, ε) and t s, (λ, ε) such that (3.8) φ s,l (ξ; λ, ε) =t s,+ (λ, ε) ψ s, ( ξ; λ, ε)+t s, (λ, ε) ψ s,+ ( ξ; λ, ε) or ξ> 1 ε up to exponentially small terms in ξ, wheret s,+ was already introduced in Lemma 3.. Proo. The ū-component o φ s,l is constant on I,sinceboth d dξ ūs and d dξ p s,+ are asymptotically small on I. Thereore, we can determine its leading order value using Corollary 3.5 and (3.1). The matrix A s as deined in (3.13) is symmetric in ξ. For the region ξ> 1 ε we can thereore use the same ψ,± and ψ s,± rom Lemma 3.4 as a basis or the solution space in this region, under the relection ξ ξ. The role o ψ s,+ and ψ s, is reversed compared to the interval ξ< 1 ε :weseethatψ s, ( ξ) grows (slowly) exponentially as ξ,whereas ψ s,+ ( ξ) has an exponential (slow) decay under the same limit. The normalization o φ s,l or ξ, which by Lemma 3. introduces t s,+ (λ, ε) in(3.8), does not exclude the possibility that or ξ> 1 ε, φ s,l has components which decay (slowly) as ξ. Thereore, we write the leading order expression o φ s,l in this region as a linear combination o a slowly increasing and a slowly decreasing component and introduce t s, (λ, ε) to measure the decreasing component. A term containing the ast decreasing component is omitted, since or ξ> 1 ε this would give only an exponentially small correction to the result in (3.8). Based on the results o Lemma 3.6, wehave (3.9) lim ū s (χ) =Γ(1+μ) [ t s,+ (λ, ε) Pν μ ( z )+t s, (λ, ε) Pν μ (z ) ] χ 0 to leading order. Corollary 3.7. Combining (3.7) and (3.9) yields (3.30) t s,+ (λ, ε) P μ ν ( z )+t s, (λ, ε) Pν μ (z )= Pν μ (z )+O(ε ε) to leading order. This gives a (irst) relation between t s,+ (λ, ε) andt s, (λ, ε) The ast components o φ s,l inside I. Since ū s (ξ; λ, ε) is constant to leading order or ξ I (see Lemma 3.6), we can represent it by its value at 0 I. Moreover, the equation or the v-component in (3.1) decouples and yields an inhomogeneous Sturm Liouville problem, ( (3.31) v ξξ (1 + λ) ) v h,0 (ξ; u, 0) v = 1 u u v h,0 (ξ; u, 0) ū s (0), where we used that u h (ξ) =u and v h (ξ) =v h,0 (ξ; u, 0) or ξ I to leading order (see Theorem.1). Basedontheslowbehavioroφ s,l determined in Lemmas 3. and 3.6, we observe that the solution v o (3.31) must extinguish as ξ I, which implies that v must decay exponentially ast in ξ.
18 PULSES IN A SLOWLY NONLINEAR GM EQUATION 45 By the nature o the GM equation (1.4) and its slow nonlinearity the problem can be solved exactly along the same lines as done in section 3. or the slow problem. First, we introduce a coordinate transormation similar to (3.17), (3.3) ζ = d dξ v h,0(ξ; u, 0) v h,0 (ξ; u, 0) using (.5). In this coordinate, v h canbewrittenas (3.33) v h (ζ; u, 0) = 3 u and (3.31) is transormed to ( (3.34) (1 ζ ) v ζζ ζ v ζ + 1 = d ( ) dξ log 1 1 v h,0 (ξ; u, 0) =tanh ξ, ( 1 ζ ) ) 4(1 + λ) 1 ζ v =9ū s (0) (1 ζ ) and I = {ζ R ζ < 1} up to exponentially small terms; compare (.9). Its homogeneous reduction ( ) (3.35) (1 ζ 4(1 + λ) ) v ζζ ζ v ζ ζ v =0 can again be solved using associated Legendre unctions; it is a special case (α =1,d =)o the slow eigenvalue problem (3.18). It should be noted that there is a crucial dierence between (3.18) and(3.34). The slow equation (3.18) is deined only on part o the ull domain: z ( 1,z ) ( 1, 1). Thereore, the eigenvalues o (3.18) do not yield direct implications or the stability o the pulse (U h,v h ). This is very dierent rom the ast system (3.34). It has three eigenvalues; its corresponding eigenunctions are (3.36a) (3.36b) (3.36c) λ (0) = 5 4, v(0) (ζ) =( 1 ζ ) 3, λ (1) =0, v (1) (ζ) =ζ ( 1 ζ ) = dζ d 3u dξ dζ v h(ζ; u, 0), λ () = 3 4, v() (ζ (ζ) = 5) 1 1 ζ. Reerring to [5], we recall that the roots o t,+ (λ, ε) are to leading order given by the eigenvalues o (3.35), so we have the ollowing lemma. Lemma 3.8. There are unique λ (i) (ε) R such that lim ε 0 λ (i) (ε) = λ (i) t,+ (λ (i) (ε),ε)=0with multiplicity 1 or i =0, 1,. Proo. See [5]. Hence, the eigenvalues o (3.35) are to leading order zeroes o the ast component o the Evans unction D(λ, ε) given in (3.11) and thus, in principle, candidates or being zeroes o the ull Evans unction. and
19 46 FRITS VEERMAN AND ARJEN DOELMAN For all λ C \ σ ess, the solution space o (3.35) is spanned by the associated Legendre unctions (3.37) v ± (ζ; λ) =c ± (λ)p 1+λ 3 (±ζ), lim v ±(ζ; λ) =0, ζ ±1 where we normalize v ± (i.e., choose c ± ) such that their Wronskian is given by (3.38) W ( v, v + )(ζ; λ) = 1 1 ζ, which implies that (3.39) c + (λ) c (λ) = 1 (4+ Γ ) 1+λ ( Γ 3+ ) 1+λ. Indeed, the expression in (3.39) haspolesatλ = λ (i), i =0, 1,. that v ± (ζ; λ) cannot span the two-dimensional solution space or λ = λ (i) This is due to the act. Since we have normalized the Wronskian (3.38), this is now encoded in the values o c ± (λ). We know that the inhomogeneous equation (3.34) has a unique bounded solution v in (ξ; λ). It can be determined using the Green s unction v (s; λ) v + (ζ; λ) W ( v (3.40) G(ζ,s; λ) =, v + )(s; λ)(1 s ), s < ζ, v (ζ; λ) v + (s; λ) W ( v, v + )(s; λ)(1 s ), s > ζ, or all λ C \ σ ess and λ λ (0,1,) so that (3.41) v in (ζ; λ) = 1 1 =9ū s (0) 9ū s (0) (1 s ) G(ζ,s; λ)ds [ v + (ζ; λ) ζ 1 (1 s ) v (s; λ)ds + v (ζ; λ) 1 ζ ] (1 s ) v + (s; λ)ds. Note that the inhomogeneous term in (3.34) is orthogonal only to the eigenunction corresponding to λ (1) = 0; or the other two eigenvalues the solvability condition (3.4) 1 1 9ū s (0) (1 ζ ) v (i) (ζ)dζ =0, i =0, 1,, is not satisied, since both v (0,) (ζ) and9(1 ζ ) are even unctions in ζ. This means that v in as a unction o λ has a simple pole at λ (0) and λ () and is smooth at λ (1) =0. To summarize this section, the resulting expression o v in is restated in the ollowing lemma. Lemma 3.9. The unique solution v in (ζ; λ) to (3.34) is given by [ ζ 1 ] (3.43) v in (ζ; λ) =9ū s (0) v + (ζ; λ) (1 s ) v (s; λ)ds + v (ζ; λ) (1 s ) v + (s; λ)ds, with v ± (ζ; λ) as deined in (3.37) and subject to condition (3.39). 1 ζ
20 PULSES IN A SLOWLY NONLINEAR GM EQUATION The slow transmission unction t s,+ (λ, ε). In section 3. we studied φ s,l outside I, and in section 3.3 we considered its ast dynamics inside I. However, we did not yet combine these results. Using (3.1), we see that (3.44) ū ξξ = ε σv h (ξ) v + O(ε 4 ) = ε σv h,0 (ξ; u, 0) v in (ζ(ξ); λ) to leading order in I. Thus, the total change o ū ξ over I is given by Δ ξ ū ξ = u ξξ dξ I (3.45) = ε σ = ε σ = ε σ = 6 ε σu 1 v h,0 (ξ; u, 0) v in (ζ(ξ); λ) dξ v h,0 (ξ(ζ); u, 0) v in (ζ; λ) 3 u (1 ζ ) v in (ζ; λ) 1 v in (ζ; λ) dζ := Δ, dζ 1 ζ dζ 1 ζ all to leading order. Using the expression or v in (ζ; λ) rom Lemma 3.9 and the symmetry in ζ between v + and v, this can be rewritten as (3.46) Δ = 108 ε σu ū s (0) 1 ζ 1 1 v + (ζ; λ) v (s; λ)(1 s )ds dζ. The desired coupling between the slow and ast dynamics can now be obtained by realizing that this change in ū ξ should match with the slow behavior o φ s,l outside I. Using Corollary 3.5 and Lemma 3.6, ( ) ( 1 ε Δ =Δ ξ ū ξ =ū ξ ū ξ 1 ) ε = ε Γ(1 + μ) dz d [ ts,+ Pν μ ( z)+t s, Pν μ (z) ] dχ dz z=z ε Γ(1 + μ) dz d [ P μ ν ( z) ] dχ dz z= z = ε Γ(1 + μ) dz d [ (3.47) ts,+ Pν μ ( z)+(t s, +1)Pν μ (z) ] dχ dz z=z to leading order. Together, expressions (3.45) and(3.47) can be used to obtain a second relation between the two transmission unctions t s,± (λ, ε); see Corollary 3.7. Thus, we can eliminate t s, and obtain a leading order expression or t s,+ : [ ] t s,+ ε dz d Pν μ ( z) P ν μ ( z ) dχ dz (z ) P ν μ Δ dz d [ (z) ε P μ ν (z) ] Γ(1 + μ) dχ dz z=z P μ ν z=z =
21 48 FRITS VEERMAN AND ARJEN DOELMAN so that, using the Wronskian W ( Pν μ (z),pν μ ( z) ) (z ), (3.48) t s,+ = Pν μ (z ) Δ dz Γ(1+μ) ε dχ ε dz dχ W ( P μ ν d dz [ P μ ν ] (z) ) (z),pν μ ( z) z=z, (z ) which, using (3.45) and(3.6), leads to the ollowing theorem. Theorem Let ε>0 be suiciently small. The unction t s,+ (λ, ε) is meromorphic as a unction o λ outside σ ess. It has simple poles at λ (0) (ε) and λ () (ε) and is analytic elsewhere. The leading order behavior o t s,+ is given by (3.49) t s,+ (λ, 0) = Pν μ (z ) αz Γ(1+μ) 1 1 v in (ζ; λ) dζ + dz 1 dz dχ W ( P μ ν d dχ dz [ (z),pν μ ( z) ] (z) ) (z ) The nontrivial roots o the Evans unction D(λ, ε) coincide with the roots o t s,+ (λ, ε). These roots determine the stability o the pulse (U h (ξ),v h (ξ)). Note that it is clear rom (3.49) thatt s,+ inherits the poles o v in at λ = λ (0,). The roots o the Evans unction D(λ, ε) outside σ ess are given by the roots o the product t,+ (λ, ε) t s,+ (λ, ε). Based on orthogonality arguments, we have established that t s,+ (λ, 0) has simple poles at λ = λ (0,) ; see the solvability condition (3.4). These coincide with the (simple) roots o t,+ (λ, 0) (see Lemma 3.8), so the Evans unction will not necessarily be zero at these values o λ. Moreover, since the Evans unction is analytic, this statement continues to hold or ε>0. Note that λ = 0 is always a trivial eigenvalue or system (3.1), with eigenunction d dξ (U h(ξ),v h (ξ)); it does not appear as a zero o t s,+ (λ, 0). 4. Pulse stability: Results. The purpose o this section is to analyze the roots o t s,+ (λ, 0) as given in Theorem The Wronskian in the denominator is always inite or 1 <z < 1, because the underlying dierential equation (3.18) is singular only at z = 1, 1. We can thereore ocus at the numerator, which is zero whenever Pν μ (z )=0or 1 αz (4.1) v in (ζ; λ) dζ + dz d [ P μ ν (z) ] Γ(1 + μ) dχ dz z=z =0. Using (4.) 1 [ ] dz = 1 dχ z=z (d 1) α(1 z ) P μ ν and d [ (4.3) P μ ν (z) ] dz z=z = 1 ( ) 1 z (ν μ)p μ ν 1 (z ) z νpν μ (z ), (4.1) can be rewritten as (4.4) 18 z P μ ν (z ) 1 ζ 1 1 v + (ζ; λ) v (s; λ)(1 s )ds dζ + 1 ( ) (d 1) (ν μ)p μ ν 1 (z ) z νpν μ (z ) =0 z=z.
22 PULSES IN A SLOWLY NONLINEAR GM EQUATION 49 LHS 5 LHS RHS Λ Figure 7. Here, LHS(λ) is plotted in blue or λ ( 1, ); the red line is the graph o RHS( λ ; ν, z ). Inthe α let plot α =0.05, ν =,andz =0.75. In the right plot α =1.5, ν =,andz = 0.60; theletigureisa representation o the statement in Theorem 4.. using (3.46) and recalling that ū s (0) = Γ(1+μ) P μ ν this equation is relevant only i Pν μ u 0;see(3.6)) to obtain the ollowing. RHS 5 (z ) 0, we divide by z P μ ν Corollary 4.1. I Pν μ in (3.11) are given to leading order by the solutions o the equation (z ) to leading order by Lemma 3.6. Since (z ) (note that z 0,since (z ) 0, the nontrivial roots o the Evans unction D(λ, ε) as deined 1 ζ (4.5) 18 v + (ζ; λ) v (s; λ)(1 s )dsdζ = ν 1 ( ν (ν μ) P μ ν 1 (z ) ) z Pν μ, (z ) with μ, ν, z as given in (3.19) and (3.6). The let-hand side (LHS) o this equation is a unction o λ only; all parameters are contained in the right-hand side (RHS). Moreover, we have restricted our parameter space (α, γ, σ, d) R >0 R >0 R \{0} (1, ), a union o two orthants in R 4,to(α, ν, z ) R >0 (1, ) ( 1, 0) (0, 1), the union o two (semicompact) slabs in R 3. It is useul to deine the LHS and RHS o (4.5) separately: 1 ζ (4.6) LHS(λ) =18 v + (ζ; λ) v (s; λ)(1 s )dsdζ, 1 1 ( ) (4.7) RHS( λ α ; ν, z )= 1 ν 1 ν (ν μ( λ α ; ν)) P μ( λ α ;ν) ν 1 (z ) z P μ( λ α ;ν) ν (z ) In Figure 7, the graphs o LHS(λ) and RHS( λ α ; ν, z ) are plotted or real values o λ. It is worthwhile to note that LHS(λ) = 88 R(P = 1+λ;, ), as used in [5] Immediate results: σ < 0 and γ 0. In this subsection we present the irst immediate implications o the developed theory or the stability o the pulse (U h,v h ). Theorem 4.. Let ε>0 be suiciently small. For all σ<0, there is a real zero λ pos > λ (0) > 0 o the Evans unction associated with the stability problem (3.1). Proo. As λ,rom(3.19) weinerthatμ ν such that the ratio P μ ν 1 (z ) 1. μ 1 z 1. P μ ν (z ) Thereore, RHS( λ α ; ν, z ) ν 1 3σu λ as λ. Using an argument equivalent to that in [1, Lemma 4.1(ii)], one can show that LHS(λ) increases monotonically (rom ) to
23 50 FRITS VEERMAN AND ARJEN DOELMAN zero or λ>λ (0). Thereore, there is a λ> 5 4 or which the LHS and RHS intersect and which thereore solves (4.5) or all parameter values when σ<0; see Figure 7. Corollary 4.3. A pulse with a double hump in the U-component, as shown in Figure 1(b), is always unstable. A direct consequence o the above corollary is that in order to obtain any stability result, we have to conine ourselves to the interval 0 <z < 1, since sgn(z )=sgn(σ); see (3.6). It would be beneicial or a complete understanding o the linear stability o the constructed pulse i more would be known about the zeroes o Pν μ (z ). However, while some inormation can be obtained regarding the number o zeroes o Pν μ (z ) or real values o μ (see [1]; the general case will be treated in [13]), the authors are not aware o any general analytic expressions concerning zeroes o Pν μ (z )orcomplexμ. Notwithstanding, direct numerical evaluation o Pν μ (z ) or a broad parameter range has led to the ollowing conjecture. Conjecture 4.4. For all λ C or which Re λ>0, Pν μ (z ) 0or all 0 <z < 1. Moreover, or Im λ 0, Pν μ (z ) 0or all 0 < z < 1. Based on this observation, the study o linear stability o the pulse can be conined to the study o solutions o (4.5). Moreover, any additional eigenvalues originating rom zeroes o Pν μ (z ) would occur on the real line and be negative. Note that in the ollowing results, this conjecture is not needed. Equation (4.7) can be studied or dierent parameter values (and limits thereo) to obtain inormation about the pulse spectrum. Another direct result can be obtained by taking the limit γ 0 to remove the inluence o the slow nonlinearity in (1.7) and obtain the classical GM equations. α 3 σ (see (.13)), so z sgn(σ) using(3.6). Note that, while the limit As γ 0, u γ 0 reduces (1.7) to the classical GM equation where the slow evolution in U is linear, yielding a simple exponential instead o an associated Legendre unction the coordinate z is ill-deined or γ = 0; see(3.17) inrelationto(3.16). Thereore, some o the expressions in the ollowing will still depend on ν, while ν disappears rom (1.7) asγ 0. Since Pν μ (z ) 1 ( ) μ 1 z (4.8a) as z 1, Γ(1 + μ) ( ) Pν μ Γ(μ) 1+z μ (4.8b) (z ) Γ(μ ν)γ(1 + μ + ν) (see [, 1]), this means that (4.9) lim RHS( λ γ 0 α ; ν, z (α, γ, σ, d)) = as z 1 lim RHS( λ α ; ν, z )=sgn(σ) z sgn(σ) =sgn(σ) 1+ λ α. Moreover, Pν μ (z ) can be written as Pν μ (z )= ( 1 z ) μ expansion (see [, 1]). Near z =1,F (z) can be expanded as F (z) = ( μ ) k j (4.10) a k = ( ν) j(ν +1) j Γ(1 + j + μ)(k j)!j!. j=0 μ ν 1 F (z ), where F (z) has a regular ) k,with k=0 a k ( 1 z
FRITS VEERMAN, ARJEN DOELMAN
PULSES IN A GIERER-MEINHARDT EQUATION WITH A SLOW NONLINEARITY FRITS VEERMAN, ARJEN DOELMAN Abstract. In this paper, we study in detail the existence and stability o localized pulses in a Gierer- Meinhardt
More informationLarge Stable Pulse Solutions in Reaction-diffusion Equations
Large Stable Pulse Solutions in Reaction-diusion Equations ARJEN DOELMAN, ROBERT A. GARDNER & TASSO J. KAPER ABSTRACT. In this paper we study the existence and stability o asymptotically large stationary
More informationDestabilization mechanisms of periodic pulse patterns near a homoclinic limit
Destabilization mechanisms of periodic pulse patterns near a homoclinic limit Arjen Doelman, Jens Rademacher, Björn de Rijk, Frits Veerman Abstract It has been observed in the Gierer-Meinhardt equations
More information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationDestabilization mechanisms of periodic pulse patterns near a homoclinic limit
Destabilization mechanisms of periodic pulse patterns near a homoclinic limit Arjen Doelman, Jens Rademacher, Björn de Rijk, Frits Veerman Abstract It has been observed in the Gierer-Meinhardt equations
More information2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction
. ETA EVALUATIONS USING WEBER FUNCTIONS Introduction So ar we have seen some o the methods or providing eta evaluations that appear in the literature and we have seen some o the interesting properties
More informationsystem May 19, 2009 MS69: New Developments in Pulse Interactions SIAM Conference on Applications of Dynamical Systems Snowbird, Utah, USA
CWI, Amsterdam heijster@cwi.nl May 9, 29 MS69: New Developments in Pulse s SIAM Conference on Applications of Dynamical Systems Snowbird, Utah, USA Joint work: A. Doelman (CWI/UvA), T.J. Kaper (BU), K.
More informationThe Dynamics of Reaction-Diffusion Patterns
The Dynamics of Reaction-Diffusion Patterns Arjen Doelman (Leiden) (Rob Gardner, Tasso Kaper, Yasumasa Nishiura, Keith Promislow, Bjorn Sandstede) STRUCTURE OF THE TALK - Motivation - Topics that won t
More informationSEPARATED AND PROPER MORPHISMS
SEPARATED AND PROPER MORPHISMS BRIAN OSSERMAN The notions o separatedness and properness are the algebraic geometry analogues o the Hausdor condition and compactness in topology. For varieties over the
More informationThe Clifford algebra and the Chevalley map - a computational approach (detailed version 1 ) Darij Grinberg Version 0.6 (3 June 2016). Not proofread!
The Cliord algebra and the Chevalley map - a computational approach detailed version 1 Darij Grinberg Version 0.6 3 June 2016. Not prooread! 1. Introduction: the Cliord algebra The theory o the Cliord
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationSEPARATED AND PROPER MORPHISMS
SEPARATED AND PROPER MORPHISMS BRIAN OSSERMAN Last quarter, we introduced the closed diagonal condition or a prevariety to be a prevariety, and the universally closed condition or a variety to be complete.
More informationGroup Method. December 16, Oberwolfach workshop Dynamics of Patterns
CWI, Amsterdam heijster@cwi.nl December 6, 28 Oberwolfach workshop Dynamics of Patterns Joint work: A. Doelman (CWI/UvA), T.J. Kaper (BU), K. Promislow (MSU) Outline 2 3 4 Interactions of localized structures
More informationarxiv: v1 [math.ap] 5 Feb 2015
Spectra and stability of spatially periodic pulse patterns: Evans function factorization via Riccati transformation Björn de Rijk brijk@math.leidenuniv.nl Jens Rademacher rademach@math.uni-bremen.de Arjen
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More informationMS: Nonlinear Wave Propagation in Singular Perturbed Systems
MS: Nonlinear Wave Propagation in Singular Perturbed Systems P. van Heijster: Existence & stability of 2D localized structures in a 3-component model. Y. Nishiura: Rotational motion of traveling spots
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationThe Deutsch-Jozsa Problem: De-quantization and entanglement
The Deutsch-Jozsa Problem: De-quantization and entanglement Alastair A. Abbott Department o Computer Science University o Auckland, New Zealand May 31, 009 Abstract The Deustch-Jozsa problem is one o the
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationVALUATIVE CRITERIA FOR SEPARATED AND PROPER MORPHISMS
VALUATIVE CRITERIA FOR SEPARATED AND PROPER MORPHISMS BRIAN OSSERMAN Recall that or prevarieties, we had criteria or being a variety or or being complete in terms o existence and uniqueness o limits, where
More informationMath 248B. Base change morphisms
Math 248B. Base change morphisms 1. Motivation A basic operation with shea cohomology is pullback. For a continuous map o topological spaces : X X and an abelian shea F on X with (topological) pullback
More informationRoberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points
Roberto s Notes on Dierential Calculus Chapter 8: Graphical analysis Section 1 Extreme points What you need to know already: How to solve basic algebraic and trigonometric equations. All basic techniques
More informationHOPF DANCES NEAR THE TIPS OF BUSSE BALLOONS. Arjen Doelman. (Communicated by the associate editor name)
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX HOPF DANCES NEAR THE TIPS OF BUSSE BALLOONS Arjen Doelman Mathematisch Instituut, Universiteit
More informationLecture : Feedback Linearization
ecture : Feedbac inearization Niola Misovic, dipl ing and Pro Zoran Vuic June 29 Summary: This document ollows the lectures on eedbac linearization tought at the University o Zagreb, Faculty o Electrical
More informationFeedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.
Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 216 Ravi N Banavar banavar@iitbacin September 24, 216 These notes are based on my readings o the two books Nonlinear Control
More informationNONLINEAR ASYMPTOTIC STABILITY OF THE SEMI-STRONG PULSE DYNAMICS IN A REGULARIZED GIERER-MEINHARDT MODEL
NONLINEAR ASYMPTOTIC STABILITY OF THE SEMI-STRONG PULSE DYNAMICS IN A REGULARIZED GIERER-MEINHARDT MODEL ARJEN DOELMAN, TASSO J. KAPER, AND KEITH PROMISLOW Abstract. We use renormalization group RG techniques
More informationThe achievable limits of operational modal analysis. * Siu-Kui Au 1)
The achievable limits o operational modal analysis * Siu-Kui Au 1) 1) Center or Engineering Dynamics and Institute or Risk and Uncertainty, University o Liverpool, Liverpool L69 3GH, United Kingdom 1)
More informationVALUATIVE CRITERIA BRIAN OSSERMAN
VALUATIVE CRITERIA BRIAN OSSERMAN Intuitively, one can think o separatedness as (a relative version o) uniqueness o limits, and properness as (a relative version o) existence o (unique) limits. It is not
More informationFinite Dimensional Hilbert Spaces are Complete for Dagger Compact Closed Categories (Extended Abstract)
Electronic Notes in Theoretical Computer Science 270 (1) (2011) 113 119 www.elsevier.com/locate/entcs Finite Dimensional Hilbert Spaces are Complete or Dagger Compact Closed Categories (Extended bstract)
More information1 Introduction Self-replicating spots and pulses have been observed in excitable reaction-diusion systems [22, 17, 24, 23, 16, 9, 2, 3, 4, 25, 21, 19,
Slowly-modulated two pulse solutions and pulse splitting bifurcations Arjen Doelman Korteweg-deVries Instituut Universiteit van Amsterdam Plantage Muidergracht 24 1018TV Amsterdam, The Netherlands Wiktor
More informationNONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION
NONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION Steven Ball Science Applications International Corporation Columbia, MD email: sball@nmtedu Steve Schaer Department o Mathematics
More informationFront Speeds, Cut-Offs, and Desingularization: A Brief Case Study
Contemporary Mathematics Front Speeds, Cut-Offs, and Desingularization: A Brief Case Study Nikola Popović Abstract. The study of propagation phenomena in reaction-diffusion systems is a central topic in
More informationCHAPTER 5 Reactor Dynamics. Table of Contents
1 CHAPTER 5 Reactor Dynamics prepared by Eleodor Nichita, UOIT and Benjamin Rouben, 1 & 1 Consulting, Adjunct Proessor, McMaster & UOIT Summary: This chapter addresses the time-dependent behaviour o nuclear
More information(C) The rationals and the reals as linearly ordered sets. Contents. 1 The characterizing results
(C) The rationals and the reals as linearly ordered sets We know that both Q and R are something special. When we think about about either o these we usually view it as a ield, or at least some kind o
More information1. Introduction. In this paper, we continue our study of slowly modulated two-pulse solutions in the one-dimensional Gray Scott model: V t = D 2 V
SIAM J. APPL. MATH. Vol. 61, No. 6, pp. 2036 2062 c 2001 Society for Industrial and Applied Mathematics SLOWLY MODULATED TWO-PULSE SOLUTIONS IN THE GRAY SCOTT MODEL II: GEOMETRIC THEORY, BIFURCATIONS,
More informationBasic mathematics of economic models. 3. Maximization
John Riley 1 January 16 Basic mathematics o economic models 3 Maimization 31 Single variable maimization 1 3 Multi variable maimization 6 33 Concave unctions 9 34 Maimization with non-negativity constraints
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationCategories and Natural Transformations
Categories and Natural Transormations Ethan Jerzak 17 August 2007 1 Introduction The motivation or studying Category Theory is to ormalise the underlying similarities between a broad range o mathematical
More information( x) f = where P and Q are polynomials.
9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More informationSIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )
SIO B, Rudnick! XVIII.Wavelets The goal o a wavelet transorm is a description o a time series that is both requency and time selective. The wavelet transorm can be contrasted with the well-known and very
More informationLeast-Squares Spectral Analysis Theory Summary
Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,
More informationSupplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values
Supplementary material or Continuous-action planning or discounted ininite-horizon nonlinear optimal control with Lipschitz values List o main notations x, X, u, U state, state space, action, action space,
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationRoot Arrangements of Hyperbolic Polynomial-like Functions
Root Arrangements o Hyperbolic Polynomial-like Functions Vladimir Petrov KOSTOV Université de Nice Laboratoire de Mathématiques Parc Valrose 06108 Nice Cedex France kostov@mathunicer Received: March, 005
More informationTHE CRITICAL WAVE SPEED FOR THE FISHER-KOLMOGOROV-PETROWSKII-PISCOUNOV EQUATION WITH CUT-OFF
THE CRITICAL WAVE SPEED FOR THE FISHER-KOLMOGOROV-PETROWSKII-PISCOUNOV EQUATION WITH CUT-OFF FREDDY DUMORTIER, NIKOLA POPOVIĆ, AND TASSO J. KAPER Abstract. The Fisher-Kolmogorov-Petrowskii-Piscounov (FKPP)
More informationTraveling stripes in the Klausmeier model of vegetation pattern formation
Traveling stripes in the Klausmeier model of vegetation pattern formation Paul Carter Arjen Doelman Abstract The Klausmeier equation is a widely studied reaction-diffusion-advection model of vegetation
More information9.3 Graphing Functions by Plotting Points, The Domain and Range of Functions
9. Graphing Functions by Plotting Points, The Domain and Range o Functions Now that we have a basic idea o what unctions are and how to deal with them, we would like to start talking about the graph o
More informationScattering of Solitons of Modified KdV Equation with Self-consistent Sources
Commun. Theor. Phys. Beijing, China 49 8 pp. 89 84 c Chinese Physical Society Vol. 49, No. 4, April 5, 8 Scattering o Solitons o Modiied KdV Equation with Sel-consistent Sources ZHANG Da-Jun and WU Hua
More informationObjectives. By the time the student is finished with this section of the workbook, he/she should be able
FUNCTIONS Quadratic Functions......8 Absolute Value Functions.....48 Translations o Functions..57 Radical Functions...61 Eponential Functions...7 Logarithmic Functions......8 Cubic Functions......91 Piece-Wise
More informationSyllabus Objective: 2.9 The student will sketch the graph of a polynomial, radical, or rational function.
Precalculus Notes: Unit Polynomial Functions Syllabus Objective:.9 The student will sketch the graph o a polynomial, radical, or rational unction. Polynomial Function: a unction that can be written in
More informationLECTURE 6: FIBER BUNDLES
LECTURE 6: FIBER BUNDLES In this section we will introduce the interesting class o ibrations given by iber bundles. Fiber bundles lay an imortant role in many geometric contexts. For examle, the Grassmaniann
More informationFinding all Bessel type solutions for Linear Differential Equations with Rational Function Coefficients
Finding all essel type solutions or Linear Dierential Equations with Rational Function Coeicients Mark van Hoeij & Quan Yuan Florida State University, Tallahassee, FL 06-07, USA hoeij@math.su.edu & qyuan@math.su.edu
More informationMath 754 Chapter III: Fiber bundles. Classifying spaces. Applications
Math 754 Chapter III: Fiber bundles. Classiying spaces. Applications Laurențiu Maxim Department o Mathematics University o Wisconsin maxim@math.wisc.edu April 18, 2018 Contents 1 Fiber bundles 2 2 Principle
More informationSWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS
J. Appl. Math. & Computing Vol. 23(2007), No. 1-2, pp. 243-256 Website: http://jamc.net SWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS MIHAI POPESCU Abstract. This paper deals with determining
More informationOn High-Rate Cryptographic Compression Functions
On High-Rate Cryptographic Compression Functions Richard Ostertág and Martin Stanek Department o Computer Science Faculty o Mathematics, Physics and Inormatics Comenius University Mlynská dolina, 842 48
More informationCurve Sketching. The process of curve sketching can be performed in the following steps:
Curve Sketching So ar you have learned how to ind st and nd derivatives o unctions and use these derivatives to determine where a unction is:. Increasing/decreasing. Relative extrema 3. Concavity 4. Points
More information( 1) ( 2) ( 1) nan integer, since the potential is no longer simple harmonic.
. Anharmonic Oscillators Michael Fowler Landau (para 8) considers a simple harmonic oscillator with added small potential energy terms mα + mβ. We ll simpliy slightly by dropping the term, to give an equation
More informationsystem CWI, Amsterdam May 21, 2008 Dynamic Analysis Seminar Vrije Universiteit
CWI, Amsterdam heijster@cwi.nl May 21, 2008 Dynamic Analysis Seminar Vrije Universiteit Joint work: A. Doelman (CWI/UvA), T.J. Kaper (BU), K. Promislow (MSU) Outline 1 2 3 4 Outline 1 2 3 4 Paradigm U
More informationExtreme Values of Functions
Extreme Values o Functions When we are using mathematics to model the physical world in which we live, we oten express observed physical quantities in terms o variables. Then, unctions are used to describe
More informationUniversidad de Chile, Casilla 653, Santiago, Chile. Universidad de Santiago de Chile. Casilla 307 Correo 2, Santiago, Chile. May 10, 2017.
A Laplace transorm approach to linear equations with ininitely many derivatives and zeta-nonlocal ield equations arxiv:1705.01525v2 [math-ph] 9 May 2017 A. Chávez 1, H. Prado 2, and E.G. Reyes 3 1 Departamento
More informationNatural convection in a vertical strip immersed in a porous medium
European Journal o Mechanics B/Fluids 22 (2003) 545 553 Natural convection in a vertical strip immersed in a porous medium L. Martínez-Suástegui a,c.treviño b,,f.méndez a a Facultad de Ingeniería, UNAM,
More informationOnline Appendix: The Continuous-type Model of Competitive Nonlinear Taxation and Constitutional Choice by Massimo Morelli, Huanxing Yang, and Lixin Ye
Online Appendix: The Continuous-type Model o Competitive Nonlinear Taxation and Constitutional Choice by Massimo Morelli, Huanxing Yang, and Lixin Ye For robustness check, in this section we extend our
More informationAn Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator
An Alternative Poincaré Section or Steady-State Responses and Biurcations o a Duing-Van der Pol Oscillator Jang-Der Jeng, Yuan Kang *, Yeon-Pun Chang Department o Mechanical Engineering, National United
More informationCHAPTER 4 Reactor Statics. Table of Contents
CHAPTER 4 Reactor Statics Prepared by Dr. Benjamin Rouben, & Consulting, Adjunct Proessor, McMaster University & University o Ontario Institute o Technology (UOIT) and Dr. Eleodor Nichita, Associate Proessor,
More informationProblem Set. Problems on Unordered Summation. Math 5323, Fall Februray 15, 2001 ANSWERS
Problem Set Problems on Unordered Summation Math 5323, Fall 2001 Februray 15, 2001 ANSWERS i 1 Unordered Sums o Real Terms In calculus and real analysis, one deines the convergence o an ininite series
More information2 Coherent D-Modules. 2.1 Good filtrations
2 Coherent D-Modules As described in the introduction, any system o linear partial dierential equations can be considered as a coherent D-module. In this chapter we ocus our attention on coherent D-modules
More informationOn a Closed Formula for the Derivatives of e f(x) and Related Financial Applications
International Mathematical Forum, 4, 9, no. 9, 41-47 On a Closed Formula or the Derivatives o e x) and Related Financial Applications Konstantinos Draais 1 UCD CASL, University College Dublin, Ireland
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationCover Page. The handle holds various files of this Leiden University dissertation.
Cover Page The handle http://hdl.handle.net/1887/45233 holds various files of this Leiden University dissertation. Author: Rijk, B. de Title: Periodic pulse solutions to slowly nonlinear reaction-diffusion
More informationFigure 3.1 Effect on frequency spectrum of increasing period T 0. Consider the amplitude spectrum of a periodic waveform as shown in Figure 3.2.
3. Fourier ransorm From Fourier Series to Fourier ransorm [, 2] In communication systems, we oten deal with non-periodic signals. An extension o the time-requency relationship to a non-periodic signal
More informationBiplots in Practice MICHAEL GREENACRE. Professor of Statistics at the Pompeu Fabra University. Chapter 6 Offprint
Biplots in Practice MICHAEL GREENACRE Proessor o Statistics at the Pompeu Fabra University Chapter 6 Oprint Principal Component Analysis Biplots First published: September 010 ISBN: 978-84-93846-8-6 Supporting
More informationSimpler Functions for Decompositions
Simpler Functions or Decompositions Bernd Steinbach Freiberg University o Mining and Technology, Institute o Computer Science, D-09596 Freiberg, Germany Abstract. This paper deals with the synthesis o
More informationIn the index (pages ), reduce all page numbers by 2.
Errata or Nilpotence and periodicity in stable homotopy theory (Annals O Mathematics Study No. 28, Princeton University Press, 992) by Douglas C. Ravenel, July 2, 997, edition. Most o these were ound by
More informationBANDELET IMAGE APPROXIMATION AND COMPRESSION
BANDELET IMAGE APPOXIMATION AND COMPESSION E. LE PENNEC AND S. MALLAT Abstract. Finding eicient geometric representations o images is a central issue to improve image compression and noise removal algorithms.
More informationThu June 16 Lecture Notes: Lattice Exercises I
Thu June 6 ecture Notes: attice Exercises I T. Satogata: June USPAS Accelerator Physics Most o these notes ollow the treatment in the class text, Conte and MacKay, Chapter 6 on attice Exercises. The portions
More informationBreakdown of Pattern Formation in Activator-Inhibitor Systems and Unfolding of a Singular Equilibrium
Breakdown of Pattern Formation in Activator-Inhibitor Systems and Unfolding of a Singular Equilibrium Izumi Takagi (Mathematical Institute, Tohoku University) joint work with Kanako Suzuki (Institute for
More informationPattern formation in reaction-diffusion systems an explicit approach
Pattern formation in reaction-diffusion systems an explicit approach Arjen Doelman Mathematisch Instituut, Leiden University, P.O. Box 9512, 2300 RA Leiden, the Netherlands, doelman@math.leidenuniv.nl
More informationEXPANSIONS IN NON-INTEGER BASES: LOWER, MIDDLE AND TOP ORDERS
EXPANSIONS IN NON-INTEGER BASES: LOWER, MIDDLE AND TOP ORDERS NIKITA SIDOROV To the memory o Bill Parry ABSTRACT. Let q (1, 2); it is known that each x [0, 1/(q 1)] has an expansion o the orm x = n=1 a
More informationMath 216A. A gluing construction of Proj(S)
Math 216A. A gluing construction o Proj(S) 1. Some basic deinitions Let S = n 0 S n be an N-graded ring (we ollows French terminology here, even though outside o France it is commonly accepted that N does
More informationOn Picard value problem of some difference polynomials
Arab J Math 018 7:7 37 https://doiorg/101007/s40065-017-0189-x Arabian Journal o Mathematics Zinelâabidine Latreuch Benharrat Belaïdi On Picard value problem o some dierence polynomials Received: 4 April
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationThe Evans function and the stability of travelling waves
The Evans function and the stability of travelling waves Jitse Niesen (University of Leeds) Collaborators: Veerle Ledoux (Ghent) Simon Malham (Heriot Watt) Vera Thümmler (Bielefeld) PANDA meeting, University
More information8.4 Inverse Functions
Section 8. Inverse Functions 803 8. Inverse Functions As we saw in the last section, in order to solve application problems involving eponential unctions, we will need to be able to solve eponential equations
More informationUNDERSTANDING of unsteady combustion of solid
JOURNAL OF PROPULSION AND POWER Vol. 4, No. 5, September October 8 Analytical Solution or Pressure-Coupled Combustion Response Functions o Composite Solid Propellants Michael Shusser and Fred E. C. Culick
More informationTelescoping Decomposition Method for Solving First Order Nonlinear Differential Equations
Telescoping Decomposition Method or Solving First Order Nonlinear Dierential Equations 1 Mohammed Al-Reai 2 Maysem Abu-Dalu 3 Ahmed Al-Rawashdeh Abstract The Telescoping Decomposition Method TDM is a new
More informationA SIMPLE MODEL FOR BIOLOGICAL AGGREGATION WITH ASYMMETRIC SENSING
A SIMPLE MODEL FOR BIOLOGICAL AGGREGATION WITH ASYMMETRIC SENSING PAUL A. MILEWSKI, XU YANG Abstract. In this paper we introduce a simple continuum model or swarming o organisms in which there is a nonlocal
More informationA Brief Survey on Semi-supervised Learning with Graph Regularization
000 00 002 003 004 005 006 007 008 009 00 0 02 03 04 05 06 07 08 09 020 02 022 023 024 025 026 027 028 029 030 03 032 033 034 035 036 037 038 039 040 04 042 043 044 045 046 047 048 049 050 05 052 053 A
More informationProduct Matrix MSR Codes with Bandwidth Adaptive Exact Repair
1 Product Matrix MSR Codes with Bandwidth Adaptive Exact Repair Kaveh Mahdaviani, Soheil Mohajer, and Ashish Khisti ECE Dept, University o Toronto, Toronto, ON M5S3G4, Canada ECE Dept, University o Minnesota,
More information8. INTRODUCTION TO STATISTICAL THERMODYNAMICS
n * D n d Fluid z z z FIGURE 8-1. A SYSTEM IS IN EQUILIBRIUM EVEN IF THERE ARE VARIATIONS IN THE NUMBER OF MOLECULES IN A SMALL VOLUME, SO LONG AS THE PROPERTIES ARE UNIFORM ON A MACROSCOPIC SCALE 8. INTRODUCTION
More informationWave regimes and mass transfer in two-layer falling films
Wave regimes and mass transer in two-layer alling ilms by GÖKÇEN ÇEKİÇ A thesis submitted to The University o Birmingham or the degree o Doctor o Philosophy School o Mathematics The University o Birmingham
More informationIn many diverse fields physical data is collected or analysed as Fourier components.
1. Fourier Methods In many diverse ields physical data is collected or analysed as Fourier components. In this section we briely discuss the mathematics o Fourier series and Fourier transorms. 1. Fourier
More informationStrong Lyapunov Functions for Systems Satisfying the Conditions of La Salle
06 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 6, JUNE 004 Strong Lyapunov Functions or Systems Satisying the Conditions o La Salle Frédéric Mazenc and Dragan Ne sić Abstract We present a construction
More informationOptimal Control. with. Aerospace Applications. James M. Longuski. Jose J. Guzman. John E. Prussing
Optimal Control with Aerospace Applications by James M. Longuski Jose J. Guzman John E. Prussing Published jointly by Microcosm Press and Springer 2014 Copyright Springer Science+Business Media New York
More information2 Frequency-Domain Analysis
2 requency-domain Analysis Electrical engineers live in the two worlds, so to speak, o time and requency. requency-domain analysis is an extremely valuable tool to the communications engineer, more so
More informationStanding Waves If the same type of waves move through a common region and their frequencies, f, are the same then so are their wavelengths, λ.
Standing Waves I the same type o waves move through a common region and their requencies,, are the same then so are their wavelengths,. This ollows rom: v=. Since the waves move through a common region,
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationA PATHOLOGY OF ASYMPTOTIC MULTIPLICITY IN THE RELATIVE SETTING
A PATHOLOGY OF ASYMPTOTIC MULTIPLICITY IN THE RELATIVE SETTING JOHN LESIEUTRE Abstract. We point out an example o a projective amily π : X S, a π-pseudoeective divisor D on X, and a subvariety V X or which
More informationLecture 8 Optimization
4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional
More information