On localized solutions of chains of oscillators with cubic nonlinearity

Save this PDF as:
Size: px
Start display at page:

Download "On localized solutions of chains of oscillators with cubic nonlinearity"

Transcription

1 On localized solutions of chains of oscillators with cubic nonlinearity Francesco Romeo, Giuseppe Rega Dipartimento di Ingegneria Strutturale e Geotecnica, SAPIENZA Università di Roma, Italia Keywords: localization, wave propagation, nonlinear maps. SUMMARY. Stationary localized solutions arising from the dynamics of chains of linearly coupled mechanical oscillators characterized by on site cubic nonlinearity are addressed aiming to explore their relationship with the underlying nonlinear wave propagation regions. 1 INTRODUCTION In this study the dynamics of chains of linearly coupled mechanical oscillators characterized by on site cubic nonlinearity is addressed. It is well known that such nonlinear periodic systems possess propagation and attenuation regions in the frequency domain. Nonlinear travelling wave solutions exist inside the propagation zones the boundaries of which can be analytically determined through a map approach when dealing with a discrete lattice structure [1]. Given the nature of the considered maps, namely reversible area preserving maps, periodic orbits can be also identified by resorting to their symmetry lines, i.e. fixed points of the involutions with which the map can be factorized. The analytical description of both, propagation regions boundaries and periodic orbits, enables to predict the nature of the periodic orbits and thereby to distinguish between hyperbolic and elliptic fixed points. Such a distinction plays a crucial role for the onset of stationary localized solutions in the nonlinear discrete lattice under investigation. A number of studies have indeed recently confirmed the ability of nonlinear lattices to support time-periodic and spatially localized solutions (discrete breathers) by linking them to the map homoclinic and heteroclinic orbits [2]. According to this geometric approach discrete breathers are identified with sequences of homoclinic or heteroclinic intersection points of the mapping so that the unstable fixed points govern their existence since stable and unstable manifolds emanate from them. The cases in which the trivial solution is hyperbolic or elliptic, namely outside or inside the linear passing band, respectively, were extensively investigated to show the existence of a variety of breather solutions. Simple, soliton-like solutions were identified with sequences of homoclinic or heteroclinic primary intersection points. Less insight is available for the other unstable fixed points belonging to the symmetry lines; in this respect, information drawn by both nonlinear propagation regions and internal bifurcation thresholds, merged with the higher order symmetry lines, are useful to locate and classify the relevant localized solutions, along with their features. In Section 1 the equations of motion of the mechanical system are introduced and the main aspects of its Hamiltonian dynamics described by a planar map are recalled. In particular, the emphasis is placed on the chain propagation regions showing how they affect the nature of the fixed point of the underlying map. In Section 2 the symmetry lines of the nonlinear map are obtained by means of recurrence formulas; these lines are used to identify the periodic orbits and to give more insight on the map phase portraits. In Section 3, the nonlinear propagation regions are used to locate and classify different types of localized solutions given by sequences of homoclinic or heteroclinic intersection points of the mapping. Some remarks and future developments are discussed in the concluding section. 1

2 2 THE OSCILLATORY CHAIN MODEL A simple mass and spring chain is considered to analyze the regions of existence of nonlinear localization phenomenon and their connection with the underlying nonlinear propagation regions. Discrete models similar to the studied one are often used to describe the nonlinear dynamics of micromechanical arrays [3, 4]. In this context the oscillation of a single cantilever is modeled following the Euler-Bernoulli beam theory, characterized by the cubic hardening nonlinearity. In essence through a Galerkin approach, in which the beam displacement is expanded in terms of the first linear mode, a sdof Duffing type equation is obtained, where the oscillator displacements describe the tip displacement amplitude of the n-th cantilever. Moreover, in the simplest model, the driven cantilever array accounts for nearest neighbor interaction only. In this study we consider a chain of linearly coupled nonlinear oscillators described by the following equations of motion m q n + k 1 q n + k 3 q 3 n + k(2q n q n 1 q n+1 )=0 (1) Time periodic solutions of equation (1) are sought for by assuming the harmonic solution q n = a n cos(ωt). Equating coefficients of cos(ωt) and setting α = mω2 k 1 k 2 and β = 3k3 4k the following second-order difference equation for the stationary amplitude is obtained (α + βa 2 n )a n + a n+1 + a n 1 =0 (2) In the frequency dependent parameter α in (2), k is the linear coupling stiffness while m and k 1 represent mass and linear stiffness of the single oscillator, respectively. In the nonlinear parameter β, k 3 represents the softening (β >0) or hardening (β <0) oscillator cubic stiffness. Equation (2), relating the amplitudes a in adjacent chain sites n 1, n and n +1, can be rewritten in matrix form as a n+1 = T(a n )a n, where a n+1 =(a n+1,a n ) and a n =(a n,a n 1 ); T(a n ) is the frequency dependent nonlinear transfer matrix belonging to the class of area preserving maps such that det(dt(a n )) = 1, where DT is the Jacobian or tangent map with reciprocal eigenvalues. Therefore, in order to study the stationary wave transmission properties of the one-dimensional nonlinear chain (1) of length N, it is convenient to rely on the eigenvalues of the linearized map equations in the neighborhood of an orbit ranging from (a 0,a 1 ) to (a N 1,a N ). Indeed, according to [5], by considering the introduced map as a dynamical system where the chain index n plays the role of discrete time n, the analysis of the transmission properties is equivalent to the stability analysis of the orbits. The interest lies in the linear stability of spatially periodic orbits a n+m = a n with cycle length m. The eigenvalues of DT are determined by its trace. If the eigenvalues lie on the unit circle, then stable elliptic periodic cycles or oscillating solutions occur; if the eigenvalues are real, then unstable hyperbolic periodic cycles or exponentially increasing solutions occur. By setting a n+1 = x n+1 and a n = y n+1, the map T can be written as T : x n+1 = y n x n (βx 2 n + α) y n+1 = x n (3) The stability of period-m orbits is governed by the condition tr(dt m ) =2; so, for m =1, the stability condition gives r := {(x, α, β) x 2 + α +2=0} s := {(x, α, β) x 2 + α 2=0} (4) 2

3 Table 1: Hyperbolic (r) and reflection hyperbolic (s) thresholds obtained from tr(dt m ) =2and m-periodic orbits p. m r s p 1 ± 2 α x = y = ± β 2 2 α ±i α x = y = 2 α β 3 4 ± 2 α 1 α ± 1 α ±i α 2 α ± 2 α 2 α x = y = 1 α β x = ±y = α β and the curves r and s represent hyperbolic (λ 1 = λ 2 =1) and reflection hyperbolic (λ 1 = λ 2 = 1) thresholds, respectively [1]. For period-m orbits, with m =1, 2, 3, 4, the latter curves are reported in Table 1. 3 SYMMETRY LINES AND PERIODIC ORBITS In this section the special properties of the map (3) are exploited to define and compute its symmetry lines with a twofold goal: to identify the periodic orbits and to give more insight on the map phase portraits. We begin by considering that the planar map (3) is reversible since it satisfies the condition R T R = T 1 where denotes the composition of two maps and R is called a reversing operator of T [6]; reversible maps have a number of special properties and occur often in applications. In particular, T is an area preserving map (4) for which the reversor R is an involution I 0, namely I 2 0 = 1. The higher order involutions are defined as I m := T m I 0, m = ±1, ±2,... The fixed points of the involutions with which the map can be factorized provide the symmetry lines, defined as γ m = Fix(I m )={ξ I m (ξ) =(ξ), ξ R} (5) The set γ m of I m is called the m-th order symmetry line satisfying the recursion relation (l = 0, 1, 2,...) [6] The main symmetry lines of the map (3) are summarized in Table 2, where the change of notation = n+1, = n is introduced. By defining the spatially m-periodic orbits by p m := {ξ T m (ξ) =(ξ)} the following relation 3

4 T l (γ m )=γ 2l+m m =0 T l (γ 0 )=γ 2l even symmetry line m =1 T l (γ 1 )=γ 2l+1 odd symmetry line Table 2: Main symmetry lines of the map (3) for β>0 and β<0. β>0 I 0 : I 1 : x = y y = x x = x y = y x(α + βx 2 ) γ 0 : y = x γ 1 : y = 1 2 x(βx2 + α) β<0 I 0 I 2 : x = y y = x x = x(α + βx 2 )+y y = x x(α + βx 2 ) 2 y(α + βx 2 ) γ 0 : y = x γ 2 : y = x( 1+α + βx 2 ) between periodic orbits and symmetry lines holds γ m γ l p m l (6) The symmetry lines and the relation (6) can be used to find spatially periodic solution of any order; in the last column of Table 1 they are reported for m =1, 2, 3, 4. In Figure 1 and Figure 2 the first 10 symmetry lines, generated by the involutions reported in Table 2, are superimposed on the phase portraits of the map (3) for two values of α. The periodic solutions in the figures are discussed against the T propagation region in Fig. 3a. In Figure 1, where α = 0.5 and β =1.0, the elliptic fixed point at (0, 0) as well as the period-4 orbits, located at x = ±y = ± 0.5 (both inside the pass band in Fig. 3a), are clearly highlighted by the relevant symmetry lines intersections obeying to (6); the quasi-periodic orbits scenario surrounding these points is also clearly related to the symmetry lines (Figure 1a). In Figure 1b, the symmetry lines are shown on a wider x, y range showing a large number of intersections (i.e. periodic solutions) confined in the range 2.5 x 2.5 (light gray); these boundaries correspond to hyperbolic period-2 orbits. Figure 2 refers to α = 3.0 and β =1.0; in Figure 2a the phase portrait refers to orbits corresponding to positive initial conditions x 0,y 0. The period-1 orbit at x = y =1undergoes a quadrupling bifurcation and the arrangement of the associated quasi-periodic orbits is governed by the main symmetry lines. The global view in Figure 2b shows that the scenario is mirrored with respect to both the main diagonals and that the symmetry lines intersections (periodic solutions) occur only within two separated regions given by 5 x 1 and 1 x 5. These regions are bounded below by the elliptic period-1 orbits and above by the hyperbolic period-2 ones. The analytical description of both, propagation regions boundaries and periodic orbits, enables to predict the nature of the periodic orbits and thereby to distinguish between hyperbolic and elliptic 4

5 Figure 1: Phase portrait of the map for α = 0.5, β =1.0 and the first 10 symmetry lines generated by the involutions I 0, I 1 in Table 2. solutions. As it will be shown in the next section, such a distinction plays a crucial role for the onset of stationary localized solutions in the nonlinear discrete lattice under investigation. 4 LOCALIZED SOLUTIONS AGAINST PROPAGATION REGIONS According to [7] the discrete breathers, i.e. time-periodic spatially localized solutions, are solutions of (1) with period T b =2πΩ b that can be expanded in a Fourier series as follows q n (t) = k A kn e ikω bt with A k, n 0 (7) where the assumption of Fourier coefficients localized in space is introduced. Following the approach proposed in [8] and [9] for hardening and softening type of nonlinearity, respectively, we look for localized states by relating them to the intersections of the invariant manifolds emanating from unstable fixed points of the map (3). As known, hyperbolic fixed points P possess stable W s (P ) and unstable W u (P ) manifolds such that a point belonging to W s (P )(W u (P )) approaches P under map iteration T n for n (n ). A homoclinic point Q is defined as Q W s W u ; thus, any Q belongs to both the stable and unstable manifold of P so, when iterated either forward or backward, it will converge to P and it corresponds to a breather solution. Due to the entangled pattern of the manifolds, additional intersection points are generated implying an infinite number of different breathers and multi-breathers, where more than one site exhibits large amplitude oscillations [7]. In Figure 3 the regions of existence of localized solutions are shown on the (α, x) plane for β =1.0 (Figure 3a) and β = 1.0 (Figure 3b); these regions are classified according to the nature of the map T fixed points, namely the trivial solution (x = y =0) and the period-1 solutions (x = y = β ). For the sake of completeness, the period-2 orbits are also reported as they play a twofold role. As shown in Figures 1b and 2b, they represent the upper boundary of all the periodic solutions; moreover, their location on the x = y symmetry line affects the topology of the invariant manifolds. Regions I +,I correspond to discrete breathers associated to the homoclinic connections emanating 5

6 Figure 2: Phase portrait of the map for α = 3.0, β =1.0 and the first 10 symmetry lines generated by the involutions I 0, I 1 in Table 2. from the trivial solutions for β>0 and β<0, respectively. Region II corresponds to nested homoclinic connections arising from the period doubling bifurcations series of period-1 orbits beginning at the intersection with the propagation region boundary s; for β =1, the latter bifurcation series occurs in the range α 4.0. As shown in [10], each period doubling bifurcation can be identified by the intersections of the newborn elliptic solutions and the relevant s-type higher order internal thresholds. Region III refers to the heteroclinic connections between the hyperbolic period-1 orbits. To validate the analytical predictions, numerical simulations pertaining to localized solutions belonging to regions I +, II and III are reported in the following. Figure 3: Projection on the (α, x) plane of the regions of existence of localized solutions against the T propagation regions. The dashed and solid lines correspond to hyperbolic and elliptic points, respectively; a) β = 1.0 (softening nonlinearity); b) β = 1.0 (hardening nonlinearity). In Figure 4 the localization scenario corresponding to the parameters α = 2.8, β = 1 be- 6

7 longing to region I + is illustrated. In Figure 4a stable (blue) and unstable (red) manifolds of (0, 0) surrounding the elliptic period-1 points (green) located at x = y = ±0.894 are shown; in Figure 4b a closer view of the manifolds is displayed with two families of homoclinic orbits; the ensuing discrete breathers are also shown in Figures 4c,d. Figure 4: Stable (blue) and unstable (red) manifolds of (0, 0) at α = 2.8, β =1(region I + ) and two homoclinic orbits; a) global view; b) closer view with intersections forming two families of homoclinic orbits; c) discrete breather corresponding to black dots in b); discrete breather corresponding to white dots in b). Figure 5 refers to region II localization, characterized by nested homoclinic orbits at α = 4.05). As shown by the global view of the invariant manifolds (Figure 5a), besides the global homoclinic connection emanating from the trivial solution, two more homoclinic tangles arise. They stem from new born period-1 hyperbolic branch following the period doubling bifurcations occurring at (± 2.0, ± 2.0) for α = 4.0. The closer view in Figure 5b shows one of the two homoclinic connections in which the newborn elliptic points are located along the γ 1 (see Table 2) symmetry line (dashed line). The associated homoclinic orbit (black dots) gives rise to the multibreather displayed in Figure 5d, whereas Figure 5c refers to a breather corresponding to a global homoclinic orbit. The numerical results presented in Figure 6 are related to the localized solutions associated with the heteroclinic connections occurring in region III. In particular, the invariant manifolds of ( 2, 2) and ( 2, 2) shown in Figure 6a are obtained for the parameters α =0.0, β = 1. The black 7

8 Figure 5: Stable (blue) and unstable (red) manifolds of (0, 0) and (± 2.05, ± 2.05) at α = 4.05, β =1and homoclinic orbits; a) global view; b) closer view of the manifolds of ( 2.05, 2.05); c) discrete breather corresponding to a homoclinic orbit of (0, 0); discrete multibreather corresponding to a homoclinic orbit of ( 2.05, 2.05). dots correspond to the discrete breather depicted in Figure 6b. 5 CONCLUSIONS The nonlinear propagation regions of chains of oscillators with cubic nonlinearity are used to identify the regions of existence of discrete breathers, i.e. time periodic and spatially localized solutions, and to guide their analysis. By describing the map symmetry lines through recurrence formulas, the phase portrait topological features and the periodic orbits arrangement have been discussed. Moreover, by merging the analytical descriptions of the map symmetry lines and propagation bands, the regions where the spatially periodic orbits are hyperbolic (saddle) or elliptic are determined. Since discrete breathers correspond to sequences of homoclinic or heteroclinic intersection points, the unstable fixed points govern their existence. Besides the known regions corresponding to localized solutions stemming from the trivial fixed points, located either outside (homoclinic connections) or inside (heteroclinic connections) the linear passing band, a further region has been identified. It corresponds to nested homoclinic orbits, emanating from period-1 unstable solutions, giving rise to multibreathers. In agreement with recent literature findings, the larger region of existence of 8

9 Figure 6: Stable (blue) and unstable (red) manifolds of ( 2, 2) and ( 2, 2) at α =0.0, β = 1; a) heteroclinic orbits; b) discrete breather corresponding to the heteroclinic orbit highlighted in a). breathers for the hardening cubic nonlinearity with respect to the softening one has been confirmed. Further studies will address the stability analysis in time of the different families of localized solutions. References [1] Romeo, F., Rega, G., Wave propagation properties of chains of oscillators with cubic nonlinearities via nonlinear map approach, Chaos, Solitons and Fractals, 27, (2006). [2] Hennig, D., Rasmussen, K.., Gabriel, H. and Blow, A., Soliton-like solutions of the generalized discrete nonlinear Schrödinger equation, Phys. Rev. E, 54, (1996). [3] Sato, M., Hubbard, B.E., Sievers, A.J., Nonlinear energy localization and its manipulation in micromechanical oscillator arrays, Rev. Modern Phys., 78, , (2006). [4] Sato, M., Sievers, A.J., Driven localized excitations in the acoustic spectrum of small nonlinear macroscopic and microscopic lattices, Phys. Rev. Lett., 98, , (2007). [5] Hennig, D., Tsironis, G.P., Wave transmission in nonlinear lattices, Physics Reports, 307, , (1999). [6] Shi, Y.G., Chen, L., Reversible maps and their symmetry lines, Commun Nonlinear Sci Numer Simulat., 16, (2011). [7] Flach, S., Gorbach, A.V., Discrete breathers Advances in theory and applications, Phys. Reports, 467, (2008). [8] Bountis, T., Capel, H.W., Kollmann, M., Ross, J.C., Bergamin, J.M., van der Weele J.P., Multibreathers and homoclinic orbits in 1-dimensional nonlinear lattices, Phys. Lett. A, 268, 50-60, (2000). 9

10 [9] Panagopoulos, P., Bountis, T., Skokos, C., Existence and stability of localized oscillations in 1-dimensional lattices with soft-spring and hard-spring potentials, J. of Vib. and Acoustics, 126, , (2004). [10] Romeo, F., Rega, G., Propagation properties of bi-coupled nonlinear oscillatory chains: analytical prediction and numerical validation, Int. J. of Bif. and Chaos, 18, , (2008). 10

arxiv:nlin/ v1 [nlin.cd] 17 Jan 2002

arxiv:nlin/ v1 [nlin.cd] 17 Jan 2002 Symbolic Dynamics of Homoclinic Orbits in a Symmetric Map arxiv:nlin/0201031v1 [nlin.cd] 17 Jan 2002 Zai-Qiao Bai and Wei-Mou Zheng Institute of Theoretical Physics, Academia Sinica, Beijing 100080, China

More information

Hamiltonian Chaos and the standard map

Hamiltonian Chaos and the standard map Hamiltonian Chaos and the standard map Outline: What happens for small perturbation? Questions of long time stability? Poincare section and twist maps. Area preserving mappings. Standard map as time sections

More information

NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM

NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM CYBERNETICS AND PHYSICS, VOL. 5, NO. 3, 016, 91 95 NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM Leonid I. Manevitch Polymer and Composite Materials Department N. N. Semenov Institute

More information

Title Coupled Cantilever Arrays( Disserta. https://doi.org/ /doctor.k14

Title Coupled Cantilever Arrays( Disserta. https://doi.org/ /doctor.k14 Title Studies on the Manipulation of Intr Coupled Cantilever Arrays( Disserta Author(s) Kimura, Masayuki Citation Kyoto University ( 京都大学 ) Issue Date 29-3-23 URL https://doi.org/1.14989/doctor.k14 Right

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A

More information

Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics

Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics Chapter 23 Predicting Chaos We have discussed methods for diagnosing chaos, but what about predicting the existence of chaos in a dynamical system. This is a much harder problem, and it seems that the

More information

Bifurcations of Traveling Wave Solutions for a Generalized Camassa-Holm Equation

Bifurcations of Traveling Wave Solutions for a Generalized Camassa-Holm Equation Computational and Applied Mathematics Journal 2017; 3(6): 52-59 http://www.aascit.org/journal/camj ISSN: 2381-1218 (Print); ISSN: 2381-1226 (Online) Bifurcations of Traveling Wave Solutions for a Generalized

More information

Invariant manifolds of the Bonhoeffer-van der Pol oscillator

Invariant manifolds of the Bonhoeffer-van der Pol oscillator Invariant manifolds of the Bonhoeffer-van der Pol oscillator R. Benítez 1, V. J. Bolós 2 1 Dpto. Matemáticas, Centro Universitario de Plasencia, Universidad de Extremadura. Avda. Virgen del Puerto 2. 10600,

More information

Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber

Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber J.C. Ji, N. Zhang Faculty of Engineering, University of Technology, Sydney PO Box, Broadway,

More information

Symbolic dynamics and chaos in plane Couette flow

Symbolic dynamics and chaos in plane Couette flow Dynamics of PDE, Vol.14, No.1, 79-85, 2017 Symbolic dynamics and chaos in plane Couette flow Y. Charles Li Communicated by Y. Charles Li, received December 25, 2016. Abstract. According to a recent theory

More information

Barcelona, Spain. RTBP, collinear points, periodic orbits, homoclinic orbits. Resumen

Barcelona, Spain.   RTBP, collinear points, periodic orbits, homoclinic orbits. Resumen XX Congreso de Ecuaciones Diferenciales y Aplicaciones X Congreso de Matemática Aplicada Sevilla, 24-28 septiembre 27 (pp. 1 8) The dynamics around the collinear point L 3 of the RTBP E. Barrabés 1, J.M.

More information

Additive resonances of a controlled van der Pol-Duffing oscillator

Additive resonances of a controlled van der Pol-Duffing oscillator Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University

More information

Chaotic Synchronization of Symbolic Information. in the Discrete Nonlinear Schrödinger Equation. Abstract

Chaotic Synchronization of Symbolic Information. in the Discrete Nonlinear Schrödinger Equation. Abstract Chaotic Synchronization of Symbolic Information in the Discrete Nonlinear Schrödinger Equation. C. L. Pando L. IFUAP, Universidad Autónoma de Puebla, arxiv:nlin/3511v1 [nlin.cd] 8 May 23 Apdo. Postal J-48.

More information

Example of a Blue Sky Catastrophe

Example of a Blue Sky Catastrophe PUB:[SXG.TEMP]TRANS2913EL.PS 16-OCT-2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of

More information

Approximation of Solitons in the Discrete NLS Equation

Approximation of Solitons in the Discrete NLS Equation University of Massachusetts Amherst ScholarWorks@UMass Amherst Mathematics and Statistics Department Faculty Publication Series Mathematics and Statistics 8 Approximation of Solitons in the Discrete NLS

More information

Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model

Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model Letter Forma, 15, 281 289, 2000 Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model Yasumasa NISHIURA 1 * and Daishin UEYAMA 2 1 Laboratory of Nonlinear Studies and Computations,

More information

arxiv: v1 [nlin.cd] 20 Jul 2010

arxiv: v1 [nlin.cd] 20 Jul 2010 Invariant manifolds of the Bonhoeffer-van der Pol oscillator arxiv:1007.3375v1 [nlin.cd] 20 Jul 2010 R. Benítez 1, V. J. Bolós 2 1 Departamento de Matemáticas, Centro Universitario de Plasencia, Universidad

More information

One Dimensional Dynamical Systems

One 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 information

Chapter 3. Gumowski-Mira Map. 3.1 Introduction

Chapter 3. Gumowski-Mira Map. 3.1 Introduction Chapter 3 Gumowski-Mira Map 3.1 Introduction Non linear recurrence relations model many real world systems and help in analysing their possible asymptotic behaviour as the parameters are varied [17]. Here

More information

Section 9.3 Phase Plane Portraits (for Planar Systems)

Section 9.3 Phase Plane Portraits (for Planar Systems) Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Half-line solutions Unstable solution Stable

More information

A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term

A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science

More information

Yulin, A. V., & Champneys, A. R. (2009). Discrete snaking: multiple cavity solitons in saturable media.

Yulin, A. V., & Champneys, A. R. (2009). Discrete snaking: multiple cavity solitons in saturable media. Yulin, A. V., & Champneys, A. R. (2009). Discrete snaking: multiple cavity solitons in saturable media. Early version, also known as pre-print Link to publication record in Explore Bristol Research PDF-document

More information

Effect of various periodic forces on Duffing oscillator

Effect of various periodic forces on Duffing oscillator PRAMANA c Indian Academy of Sciences Vol. 67, No. 2 journal of August 2006 physics pp. 351 356 Effect of various periodic forces on Duffing oscillator V RAVICHANDRAN 1, V CHINNATHAMBI 1, and S RAJASEKAR

More information

TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR

TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR Samo Lasič, Gorazd Planinšič,, Faculty of Mathematics and Physics University of Ljubljana, Slovenija Giacomo Torzo, Department of Physics, University

More information

Hopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators

Hopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators The Open Acoustics Journal 8 9-3 9 Open Access Hopf ifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators Jianping Cai *a and Jianhe Shen b a Department of

More information

8.1 Bifurcations of Equilibria

8.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 information

7 Planar systems of linear ODE

7 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 information

Phonons and lattice dynamics

Phonons and lattice dynamics Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed

More information

Difference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay

Difference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay Difference Resonances in a controlled van der Pol-Duffing oscillator involving time delay This paper was published in the journal Chaos, Solitions & Fractals, vol.4, no., pp.975-98, Oct 9 J.C. Ji, N. Zhang,

More information

INVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS

INVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS International Journal of Bifurcation and Chaos, Vol. 13, No. 7 (2003) 1841 1852 c World Scientific Publishing Company INVARIANT CURVES AND FOCAL POINTS IN A LYNESS ITERATIVE PROCESS L. GARDINI and G. I.

More information

COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS

COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 9 No. III (September, 2015), pp. 197-210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL

More information

Strange dynamics of bilinear oscillator close to grazing

Strange dynamics of bilinear oscillator close to grazing Strange dynamics of bilinear oscillator close to grazing Ekaterina Pavlovskaia, James Ing, Soumitro Banerjee and Marian Wiercigroch Centre for Applied Dynamics Research, School of Engineering, King s College,

More information

Daba Meshesha Gusu and O.Chandra Sekhara Reddy 1

Daba Meshesha Gusu and O.Chandra Sekhara Reddy 1 International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations

More information

Homoclinic bifurcations in Chua s circuit

Homoclinic bifurcations in Chua s circuit Physica A 262 (1999) 144 152 Homoclinic bifurcations in Chua s circuit Sandra Kahan, Anibal C. Sicardi-Schino Instituto de Fsica, Universidad de la Republica, C.C. 30, C.P. 11 000, Montevideo, Uruguay

More information

Growth of Sobolev norms for the cubic defocusing NLS Lecture 1

Growth of Sobolev norms for the cubic defocusing NLS Lecture 1 Growth of Sobolev norms for the cubic defocusing NLS Lecture 1 Marcel Guardia February 8, 2013 M. Guardia (UMD) Growth of Sobolev norms February 8, 2013 1 / 40 The equations We consider two equations:

More information

Introduction LECTURE 1

Introduction LECTURE 1 LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in

More information

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.

More information

Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.

Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e. Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may

More information

Kronig-Penney model. 2m dx. Solutions can be found in region I and region II Match boundary conditions

Kronig-Penney model. 2m dx. Solutions can be found in region I and region II Match boundary conditions Kronig-Penney model I II d V( x) m dx E Solutions can be found in region I and region II Match boundary conditions Linear differential equations with periodic coefficients Have exponentially decaying solutions,

More information

Pseudo-Chaotic Orbits of Kicked Oscillators

Pseudo-Chaotic Orbits of Kicked Oscillators Dynamical Chaos and Non-Equilibrium Statistical Mechanics: From Rigorous Results to Applications in Nano-Systems August, 006 Pseudo-Chaotic Orbits of Kicked Oscillators J. H. Lowenstein, New York University

More information

Simulation of two-level quantum system using classical coupled oscillators. Abstract

Simulation of two-level quantum system using classical coupled oscillators. Abstract Simulation of two-level quantum system using classical coupled oscillators Namit Anand and Ashok Mohapatra School of Physical Sciences, National Institute of Science Education and Research, Bhubaneswar

More information

Symmetry Properties of Confined Convective States

Symmetry Properties of Confined Convective States Symmetry Properties of Confined Convective States John Guckenheimer Cornell University 1 Introduction This paper is a commentary on the experimental observation observations of Bensimon et al. [1] of convection

More information

Introduction to Applied Nonlinear Dynamical Systems and Chaos

Introduction to Applied Nonlinear Dynamical Systems and Chaos Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium

More information

Torus Maps from Weak Coupling of Strong Resonances

Torus Maps from Weak Coupling of Strong Resonances Torus Maps from Weak Coupling of Strong Resonances John Guckenheimer Alexander I. Khibnik October 5, 1999 Abstract This paper investigates a family of diffeomorphisms of the two dimensional torus derived

More information

AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS

AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS POINCARÉ-MELNIKOV-ARNOLD METHOD FOR TWIST MAPS AMADEU DELSHAMS AND RAFAEL RAMíREZ-ROS 1. Introduction A general theory for perturbations of an integrable planar map with a separatrix to a hyperbolic fixed

More information

TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1

TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1 COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2

More information

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip

More information

2.034: Nonlinear Dynamics and Waves. Term Project: Nonlinear dynamics of piece-wise linear oscillators Mostafa Momen

2.034: Nonlinear Dynamics and Waves. Term Project: Nonlinear dynamics of piece-wise linear oscillators Mostafa Momen 2.034: Nonlinear Dynamics and Waves Term Project: Nonlinear dynamics of piece-wise linear oscillators Mostafa Momen May 2015 Massachusetts Institute of Technology 1 Nonlinear dynamics of piece-wise linear

More information

Models Involving Interactions between Predator and Prey Populations

Models Involving Interactions between Predator and Prey Populations Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate

More information

Stabilization of Hyperbolic Chaos by the Pyragas Method

Stabilization of Hyperbolic Chaos by the Pyragas Method Journal of Mathematics and System Science 4 (014) 755-76 D DAVID PUBLISHING Stabilization of Hyperbolic Chaos by the Pyragas Method Sergey Belyakin, Arsen Dzanoev, Sergey Kuznetsov Physics Faculty, Moscow

More information

= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :

= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F : 1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change

More information

Igor A. Khovanov,Vadim S. Anishchenko, Astrakhanskaya str. 83, Saratov, Russia

Igor A. Khovanov,Vadim S. Anishchenko, Astrakhanskaya str. 83, Saratov, Russia Noise Induced Escape from Dierent Types of Chaotic Attractor Igor A. Khovanov,Vadim S. Anishchenko, Dmitri G. Luchinsky y,andpeter V.E. McClintock y Department of Physics, Saratov State University, Astrakhanskaya

More information

Dynamics of four coupled phase-only oscillators

Dynamics of four coupled phase-only oscillators Communications in Nonlinear Science and Numerical Simulation 13 (2008) 501 507 www.elsevier.com/locate/cnsns Short communication Dynamics of four coupled phase-only oscillators Richard Rand *, Jeffrey

More information

MAPPING CLASS ACTIONS ON MODULI SPACES. Int. J. Pure Appl. Math 9 (2003),

MAPPING CLASS ACTIONS ON MODULI SPACES. Int. J. Pure Appl. Math 9 (2003), MAPPING CLASS ACTIONS ON MODULI SPACES RICHARD J. BROWN Abstract. It is known that the mapping class group of a compact surface S, MCG(S), acts ergodically with respect to symplectic measure on each symplectic

More information

Hamiltonian Dynamics

Hamiltonian Dynamics Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;

More information

Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy

Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy John H. CANTRELL 1, Sean A. CANTRELL 2 1 NASA Langley Research Center, Hampton, Virginia 23681, USA 2 NLS Analytics, LLC, Glencoe,

More information

NBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011

NBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets

More information

APPLIED SYMBOLIC DYNAMICS AND CHAOS

APPLIED SYMBOLIC DYNAMICS AND CHAOS DIRECTIONS IN CHAOS VOL. 7 APPLIED SYMBOLIC DYNAMICS AND CHAOS Bai-Lin Hao Wei-Mou Zheng The Institute of Theoretical Physics Academia Sinica, China Vfö World Scientific wl Singapore Sinaaoore NewJersev

More information

7 Two-dimensional bifurcations

7 Two-dimensional bifurcations 7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed

More information

An Exactly Solvable 3 Body Problem

An Exactly Solvable 3 Body Problem An Exactly Solvable 3 Body Problem The most famous n-body problem is one where particles interact by an inverse square-law force. However, there is a class of exactly solvable n-body problems in which

More information

Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting

Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd

More information

Exploiting pattern transformation to tune phononic band gaps in a two-dimensional granular crystal

Exploiting pattern transformation to tune phononic band gaps in a two-dimensional granular crystal Exploiting pattern transformation to tune phononic band gaps in a two-dimensional granular crystal The Harvard community has made this article openly available. Please share how this access benefits you.

More information

16 Period doubling route to chaos

16 Period doubling route to chaos 16 Period doubling route to chaos We now study the routes or scenarios towards chaos. We ask: How does the transition from periodic to strange attractor occur? The question is analogous to the study of

More information

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x

More information

Nonlinear dynamics & chaos BECS

Nonlinear dynamics & chaos BECS Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes

More information

Solitons. Nonlinear pulses and beams

Solitons. Nonlinear pulses and beams Solitons Nonlinear pulses and beams Nail N. Akhmediev and Adrian Ankiewicz Optical Sciences Centre The Australian National University Canberra Australia m CHAPMAN & HALL London Weinheim New York Tokyo

More information

Sufficient conditions for a period incrementing big bang bifurcation in one-dimensional maps.

Sufficient conditions for a period incrementing big bang bifurcation in one-dimensional maps. Sufficient conditions for a period incrementing big bang bifurcation in one-dimensional maps. V. Avrutin, A. Granados and M. Schanz Abstract Typically, big bang bifurcation occur for one (or higher)-dimensional

More information

STABILITY OF ORBITS NEAR LARGE MASS RATIO BINARY SYSTEMS

STABILITY OF ORBITS NEAR LARGE MASS RATIO BINARY SYSTEMS IAA-AAS-DyCoSS2-05-08 STABILITY OF ORBITS NEAR LARGE MASS RATIO BINARY SYSTEMS Natasha Bosanac, Kathleen C. Howell and Ephraim Fischbach INTRODUCTION With recent scientific interest into the composition,

More information

A classification of explosions in dimension one

A classification of explosions in dimension one Ergod. Th. & Dynam. Sys. (29), 29, 715 731 doi:1.117/s143385788486 c 28 Cambridge University Press Printed in the United Kingdom A classification of explosions in dimension one E. SANDER and J. A. YORKE

More information

Extremely slow edge waves in mechanical graphene with rotating grains

Extremely slow edge waves in mechanical graphene with rotating grains Phononic Crystals and Acoustic Metamaterials: Paper ICA2016-104 Extremely slow edge waves in mechanical graphene with rotating grains Li-Yang Zheng (a), Vincent Tournat (b), Georgios Theocharis (c),vitalyi

More information

On Riddled Sets and Bifurcations of Chaotic Attractors

On Riddled Sets and Bifurcations of Chaotic Attractors Applied Mathematical Sciences, Vol. 1, 2007, no. 13, 603-614 On Riddled Sets and Bifurcations of Chaotic Attractors I. Djellit Department of Mathematics University of Annaba B.P. 12, 23000 Annaba, Algeria

More information

2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers 2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

More information

The SD oscillator and its attractors

The SD oscillator and its attractors The SD oscillator and its attractors Qingjie Cao, Marian Wiercigroch, Ekaterina Pavlovskaia Celso Grebogi and J Michael T Thompson Centre for Applied Dynamics Research, Department of Engineering, University

More information

Calculus and Differential Equations II

Calculus and Differential Equations II MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,

More information

Geometry of Crystal Lattice

Geometry of Crystal Lattice 0 Geometry of Crystal Lattice 0.1 Translational Symmetry The crystalline state of substances is different from other states (gaseous, liquid, amorphous) in that the atoms are in an ordered and symmetrical

More information

ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems

ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions

More information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

More information

Complex Behavior in Coupled Nonlinear Waveguides. Roy Goodman, New Jersey Institute of Technology

Complex Behavior in Coupled Nonlinear Waveguides. Roy Goodman, New Jersey Institute of Technology Complex Behavior in Coupled Nonlinear Waveguides Roy Goodman, New Jersey Institute of Technology Nonlinear Schrödinger/Gross-Pitaevskii Equation i t = r + V (r) ± Two contexts for today: Propagation of

More information

Physics 211B : Problem Set #0

Physics 211B : Problem Set #0 Physics 211B : Problem Set #0 These problems provide a cross section of the sort of exercises I would have assigned had I taught 211A. Please take a look at all the problems, and turn in problems 1, 4,

More information

Hardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field

Hardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field Int. J. Adv. Appl. Math. and Mech. 3(2) (215) 16 32 (ISSN: 2347-2529) IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Hardening nonlinearity

More information

Influence of moving breathers on vacancies migration

Influence of moving breathers on vacancies migration Physics Letters A 315 (2003) 364 371 www.elsevier.com/locate/pla Influence of moving breathers on vacancies migration J. Cuevas a,, C. Katerji a, J.F.R. Archilla a,j.c.eilbeck b,f.m.russell b a Grupo de

More information

TWO VAN DER POL-DUFFING OSCILLATORS WITH HUYGENS COUPLING

TWO VAN DER POL-DUFFING OSCILLATORS WITH HUYGENS COUPLING ENOC-8, Saint Petersburg, Russia, June, 3 July, 8 TWO VAN DER POL-DUFFING OSCILLATORS WITH HUYGENS COUPLING V.N. Belykh, E.V. Pankratova Mathematics Department Volga State Academy Nizhny Novgorod, 3, Russia.

More information

Two models for the parametric forcing of a nonlinear oscillator

Two models for the parametric forcing of a nonlinear oscillator Nonlinear Dyn (007) 50:147 160 DOI 10.1007/s11071-006-9148-3 ORIGINAL ARTICLE Two models for the parametric forcing of a nonlinear oscillator Nazha Abouhazim Mohamed Belhaq Richard H. Rand Received: 3

More information

UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM

UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume 00, Number 0, Xxxx XXXX pp. 000 000 UNIFORM SUBHARMONIC ORBITS FOR SITNIKOV PROBLEM CLARK ROBINSON Abstract. We highlight the

More information

On low speed travelling waves of the Kuramoto-Sivashinsky equation.

On low speed travelling waves of the Kuramoto-Sivashinsky equation. On low speed travelling waves of the Kuramoto-Sivashinsky equation. Jeroen S.W. Lamb Joint with Jürgen Knobloch (Ilmenau, Germany) Marco-Antonio Teixeira (Campinas, Brazil) Kevin Webster (Imperial College

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,

More information

Computational Neuroscience. Session 4-2

Computational Neuroscience. Session 4-2 Computational Neuroscience. Session 4-2 Dr. Marco A Roque Sol 06/21/2018 Two-Dimensional Two-Dimensional System In this section we will introduce methods of phase plane analysis of two-dimensional systems.

More information

Lecture 2. Contents. 1 Fermi s Method 2. 2 Lattice Oscillators 3. 3 The Sine-Gordon Equation 8. Wednesday, August 28

Lecture 2. Contents. 1 Fermi s Method 2. 2 Lattice Oscillators 3. 3 The Sine-Gordon Equation 8. Wednesday, August 28 Lecture 2 Wednesday, August 28 Contents 1 Fermi s Method 2 2 Lattice Oscillators 3 3 The Sine-Gordon Equation 8 1 1 Fermi s Method Feynman s Quantum Electrodynamics refers on the first page of the first

More information

dynamical zeta functions: what, why and what are the good for?

dynamical zeta functions: what, why and what are the good for? dynamical zeta functions: what, why and what are the good for? Predrag Cvitanović Georgia Institute of Technology November 2 2011 life is intractable in physics, no problem is tractable I accept chaos

More information

Delay Coordinate Embedding

Delay Coordinate Embedding Chapter 7 Delay Coordinate Embedding Up to this point, we have known our state space explicitly. But what if we do not know it? How can we then study the dynamics is phase space? A typical case is when

More information

Bifurcation Trees of Periodic Motions to Chaos in a Parametric, Quadratic Nonlinear Oscillator

Bifurcation Trees of Periodic Motions to Chaos in a Parametric, Quadratic Nonlinear Oscillator International Journal of Bifurcation and Chaos, Vol. 24, No. 5 (2014) 1450075 (28 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500758 Bifurcation Trees of Periodic Motions to Chaos

More information

Hamiltonian dynamics of breathers with third-order dispersion

Hamiltonian dynamics of breathers with third-order dispersion 1150 J. Opt. Soc. Am. B/ Vol. 18, No. 8/ August 2001 S. Mookherjea and A. Yariv Hamiltonian dynamics of breathers with third-order dispersion Shayan Mookherjea* and Amnon Yariv Department of Applied Physics

More information

How Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos?

How Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos? Progress of Theoretical Physics Supplement No. 150, 2003 439 How Does the Smaller Alignment Index (SALI) Distinguish Order from Chaos? Charalampos Skokos, 1,2, ) Chris Antonopoulos, 1, ) Tassos C. Bountis

More information

Breather Modes Induced by Localized RF Radiation: Analytical and Numerical Approaches

Breather Modes Induced by Localized RF Radiation: Analytical and Numerical Approaches Proceedings of the 5th International Conference on Nonlinear Dynamics ND-KhPI2016 September 27-30, 2016, Kharkov, Ukraine Breather Modes Induced by Localized RF Radiation: Analytical and Numerical Approaches

More information

On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I

On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I IOP PUBLISHING Nonlinearity 2 (28) 923 972 NONLINEARITY doi:.88/95-775/2/5/3 On dynamical properties of multidimensional diffeomorphisms from Newhouse regions: I S V Gonchenko, L P Shilnikov and D V Turaev

More information

EXACT DARK SOLITON, PERIODIC SOLUTIONS AND CHAOTIC DYNAMICS IN A PERTURBED GENERALIZED NONLINEAR SCHRODINGER EQUATION

EXACT DARK SOLITON, PERIODIC SOLUTIONS AND CHAOTIC DYNAMICS IN A PERTURBED GENERALIZED NONLINEAR SCHRODINGER EQUATION CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 1, Spring 9 EXACT DARK SOLITON, PERIODIC SOLUTIONS AND CHAOTIC DYNAMICS IN A PERTURBED GENERALIZED NONLINEAR SCHRODINGER EQUATION JIBIN LI ABSTRACT.

More information

DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS

DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University

More information

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

More information

Part 1: Fano resonances Part 2: Airy beams Part 3: Parity-time symmetric systems

Part 1: Fano resonances Part 2: Airy beams Part 3: Parity-time symmetric systems Lecture 3 Part 1: Fano resonances Part 2: Airy beams Part 3: Parity-time symmetric systems Yuri S. Kivshar Nonlinear Physics Centre, Australian National University, Canberra, Australia http://wwwrsphysse.anu.edu.au/nonlinear/

More information

Simplest Chaotic Flows with Involutional Symmetries

Simplest Chaotic Flows with Involutional Symmetries International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries

More information