arxiv: v3 [nlin.ps] 16 Nov 2018

Transcription

1 arxiv:84.768v [nlin.ps] 6 Nov 8 Localized solutions of nonlinear network wave equations J. G. Caputo, I. Khames, A. Knippel and A. B. Aceves Laboratoire de Mathématiques, INSA de Rouen Normandie, 768 Saint-Etienne du Rouvray, France. Department of Mathematics, Southern Methodist University, Dallas, TX 7575, USA. Abstract We study localized solutions for the nonlinear graph wave equation on finite arbitrary networks. Assuming a large amplitude localized initial condition on one node of the graph, we approximate its evolution by the Duffing equation. The rest of the network satisfies a linear system forced by the excited node. This approximation is validated by reducing the nonlinear graph wave equation to the discrete nonlinear Schrödinger equation and by Fourier analysis. Finally, we examine numerically the condition for localization in the parameter plane, coupling versus amplitude and show that the localization amplitude depends on the maximal normal eigenfrequency. Introduction Intrinsic localized modes, also called discrete breathers are time periodic and (typically exponentially) spatially localized solutions that appear in nonlinear discrete systems arising in many physical, biological systems and networks. They were analytically studied first by Sievers and Takeno [], Page []. Later, MacKay and Aubry [] proved the existence of discrete breathers by considering a lattice model of coupled anharmonic oscillators in the limit of very weak interaction (anticontinuous limit). Localized solutions exist in nonlinear networks because of the interplay between nonlinearity and discreteness. In fact, the non-resonance of the breather frequency with the linear spectrum is a necessary condition for obtaining a time-periodic localized state [4]. caputo@insa-rouen.fr imene.khames@insa-rouen.fr arnaud.knippel@insa-rouen.fr aaceves@mail.smu.edu

2 Localized modes have been investigated theoretically and numerically for a variety of physical systems [5]. Experimentally observed reports for various systems include Josephson-junction arrays [6], optical waveguides [7, 8], photonic crystals [9], DNA double strand [], micromechanical oscillators [] and electronic circuits []. Recently, photonic lattices built with different configurations such as multi-core fibers [, 4] and waveguide arrays [8, ], allow the enhancement of cubic (Kerr), quadratic, photorefractive, local and non-local nonlinearities [5] under different discrete geometries [6]. Novel graphene arrays [7] suggest the existence of localized modes in the nano-meter scale. In all cases the underlying periodicity or discreteness leads to new families of optical discrete modes different from the ones of continuous systems. Localized modes have been thoroughly studied in the Fermi-Pasta-Ulam lattice [8] and in the discrete nonlinear Schrödinger equation [9]. They have been found in the dynamics of such complex objects as molecules []. Here we aim to describe the dynamics of arbitrary topology networks. The physical systems are electrical networks of nonlinear elements, for example diodes or Josephson junctions connected by inductances and systems of identical particles coupled by mechanical interactions (see the book by Scott [] for examples) and fluid flows on a network []. The general model for these systems is the graph wave equation, where the usual continuum Laplacian is replaced by the graph Laplacian []. It arises naturally from discrete conservation laws [4]. The graph Laplacian, being a symmetric positive matrix, has real eigenvalues and we can choose a basis of orthogonal eigenvectors corresponding to the normal modes [4] giving rise to periodic solutions. When nonlinearity is present, these normal modes generally couple together. Only special eigenvectors with coordinates in {,,} extend into nonlinear periodic orbits [5]. Other nonlinear solutions exist for lattices, in particular the ones that are localized. It is then natural to study them in general networks. In the present article, we study the nonlinear graph wave equation on a network of arbitrary topology and search for nonlinear localized solutions by introducing a large initial condition on only one of the nodes of network. We observe localized modes, i.e. large amplitude nonlinear excitations that do not decay significantly over thousands of periods. The Fourier spectrum reveals the presence of linear modes that are not eigenfrequencies of the Laplacian but of a reduced matrix that we characterize. Our main result is a full picture of a discrete breather on a general network. The excited node satisfies a Duffing equation while the rest of the network obeys a driven linear system which we present. We confirm this approximation by numerical simulations and by modulation theory for a natural frequency ω. Then the graph nonlinear wave equation reduces to a graph nonlinear Schrödinger equation. We illustrate the analysis by considering two graphs, the paw graph and cycle 6. The localization threshold was estimated in the parameter plane: coupling versus amplitude. As expected, waves delocalize as the coupling increases because we get closer to the resonance condition. However, different nodes behave differently depending on their degree and the graph configuration.

3 The article is organized as follows: We introduce the graph nonlinear wave equation and the localized modes in section. In section, we reduce the graph nonlinear wave equation for ω to a discrete nonlinear Schrödinger equation and determine nonlinear localized solutions. Section 4 confirms this analysis by studying the dynamics in real and Fourier space of two main networks; in particular we examine the localization vs delocalization regimes in the parameter plane coupling vs amplitude. Section 5 summarizes the paper. The graph nonlinear wave equation : Localized modes We study the nonlinear wave equation on a connected graph G with N nodes d u dt = ( ω I+ǫ ) u u, () whereu = (u (t),u (t),...,u N (t)) T is the field amplitude, u = (u,u,...,u N )T, I is the identity matrix and is the graph Laplacian []. This N N matrix is = A D, where A is the adjacency matrix such that A ij = if nodes i and j are connected (i j) and A ij = otherwise, and D is the diagonal matrix where the entry d i = N j= A ij is the degree of vertex i. For weighted graphs the in the adjacency matrix will be replaced by the weight α ij of the edge e i,j. The other parameters are the natural frequency ω and the linear coupling coefficient ǫ which is assumed small in this article. Equation () is well-posed because the cubic on-site nonlinearity guarantees the existence of the solution for all times. Note that we use bold-face capitals for matrices and bold-face lower-case letters for vectors. Equation () is an extension to a graph of the Φ 4 well-known model in condensed matter physics []. In the literature, for example [6] and references therein, the discrete Φ 4 model was studied only in the particular case of lattices where the graph Laplacian ( u) i = u i+ u i + u i (for a one-dimensional lattice) is a finite difference discretization of the continuous Laplacian. This formulation is natural since the linear graph wave equation arises from discrete conservation laws [4]. We formulate the discrete Φ 4 model using the graph Laplacian to describe general networks of arbitrary topology, like for example an electrical network. The model can describe networks of nonlinear oscillators, such as Josephson junctions or diodes [6]. In mechanical engineering, several aerospace structures e.g. turbine rotors or space antennas, are composed of weakly coupled sectors assembled in a cyclic and symmetric configuration. Such a complex system can be reduced to () [7]. The graph Laplacian is a real symmetric negative-semi definite matrix. It has real non positive eigenvalues ν i where ν = < ν ν N.

4 In our previous work [5], we constructed nonlinear periodic orbits which are extension of some linear normal modes (associated to the eigenvectors) of the graph Laplacian. Here instead, we take a different approach, we assume a large amplitude localized initial condition and search for nonlinear localized solutions. An important remark is that this work can be generalized to any odd power of the nonlinearity.. Natural frequency ω = We consider a large amplitude initial condition localized at node j and examine its evolution. First, we consider the anti-continuum limit, ǫ =, the evolution of u j satisfies d u j dt = u j, () where u j () = ρ. The other nodes u k, verify u k () = and therefore u k (t) =. The solution of () can be written in terms of the Jacobi cosine elliptic function [8] ( u j (t) = ρ cn ρt, ), () where the modulus of cn is κ = (Appendix A) and where we assumed du j dt () =. The period of oscillations (Appendix B) is T = Γ( ) 4 ρ π, (4) where Γ(.) is the gamma function and Γ ( 4).656. The frequency of oscillations is Ω = π = π π T Γ ( )ρ. (5) Now examine the weak coupling limit ǫ. The nearest neighbors k of j solve the forced system d u N k dt = ǫ kp u p u k = ǫd ku k +ǫu j +ǫ p= 4 p k, p j u p u k, (6) where d k is the degree of the node k, the notation p k indicates the adjacency of vertices and the sum is taken over the other neighbors p of k. We assume that u k is small and will find a condition on ρ for this to hold. If u k is small, it is natural to neglect the cubic term u k. The part of the solution for u k due to the forcing is d u f ( k dt = ǫu j (t) = ǫρ cn ρt, ), (7) 4

5 where the forcing u j is periodic of frequency Ω and amplitude ρ. Then, the response u k to this periodic forcing will be of amplitude ( ) u f ǫ k = O. (8) ρ This amplitude is small if ρ. Similarly, the next nearest neighbors l of node j exhibit a forced oscillation given by d u f l dt u f l = O = ǫu k (t) and this gives ( ) ǫ. (9) For simplicity and without loss of generality, we assume an initial excitation of node j =. The evolution of the nodes {,...,N} is described by the forced system of linear ordinary differential equations ρ d v dt = ǫ v+f, () where v = (u,u,...,u N ) T, is the matrix obtained by removing the first line and the first column from the graph Laplacian and wheref = (f,f,...,f N ) T is the forcing term such that f k = ǫu if k adjacent to (k ) and otherwise. The matrix is a reduction of the graph Laplacian. It is therefore real symmetric and negative, then has real eigenvalues > ω ω ω N and a basis of orthonormal eigenvectors z,z,...z N. These verify z m = ω m z m, for m {,...,N }. We expand v using a basis of the eigenvectors z m as v = N m= a m z m. () Substituting () into () and projecting on each eigenvector z m, we get d N a m dt = ǫωma m + f p zp m, where we have used the orthonormality of the eigenvectors of. The sum can be written as N f p zp m = ǫu zk m. p= We then get a set of (N ) second order inhomogeneous ordinary differential equations d a m dt = ǫωma m +ǫu zk, m () 5 k p= k

6 where m {,...,N }. At this level, we just rewrote equation () in the basis z m. Initially, a m () = dam dt () = so that we only observe the forced response of the system. In particular, the modes a m such that k zm k = will remain zero. This reveals that the harmonic frequencies of the solutions u,...,u N that will be observed are ǫ ωm, () for m {,...,N } such that k zm k. This is the sum of the eigenvector components on the neighboring nodes of the excited node. In general, for initial excitation of node j, the harmonic frequencies () of the matrix j (obtained by removing the j line and the j column from the graph Laplacian ) will be observed if zk m + zk m, (4) k j, k<j k j, k>j where z m are the eigenvectors of j. This condition on the eigenvectors of j depends on the topology. For example, it could be not satisfied when there are symmetries in the graph. Adding weights will usually break the symmetries, as shown below. Then the condition would be satisfied for all eigenvectors and all eigenfrequencies would be observed. Following the interlacing theorem [9], the eigenvalues ω, ω,..., ω N of j, the submatrix of the graph Laplacian, interlace the eigenvalues of, ν, ν,..., ν N as ν ω ν ω ν ν N ω N ν N. (5) For a general graph, the spectrum of the graph Laplacian needs to be computed numerically and the spectrum of j as well. For some special cases however, like cycles, chains and grids, the eigenvalues and eigenvectors of the Laplacian have an explicit formula. In [], we determined a part of the spectrum of for special configurations, through the graph topology.. Natural frequency ω Now, we consider the equation () with a natural frequency ω. In the anti-continuum limit ǫ =, the evolution at the excited node j satisfies d u j dt = ω u j u j, (6) whereu j () = ρ. The solution can be written in terms of cosine elliptic functions [8] ) u j (t) = ρ cn( ω +ρ t, κ, (7) 6

7 where the modulus κ = ρ (ω +ρ ) duj and we assumed dt () =. As above, we examine the weak coupling limit ǫ and assume an excitation at node j =. Following the same procedure, the evolution at nodes {,...,N} is described by the forced system of linear ordinary differential equations d v dt = ( ǫ ω I ) v+f. (8) Substituting () into (8) and projecting on each eigenvector z m of, we get d a m dt = ( ǫωm +ω) a m +ǫu zk m, (9) where m {,...,N }. The harmonic frequencies are ǫω m +ω, () for m {,...,N } such that k zm k. k Modulation theory ω When the natural frequency is not zero, following [], we reduce the discrete Φ 4 equation () to the discrete nonlinear Schrödinger equation. We write u(t) = ǫψ(t)e iωt + ǫψ (T)e iωt, () where T = ǫt, ψ = (ψ (t),ψ (t),...,ψ N (t)) T is the field vector and ψ is the complex conjugate of ψ. Plugging () into () and collecting terms in order of ǫ,ǫ,..., we obtain for the order ǫ the graph nonlinear Schrödinger equation (see Appendix C) iω dψ ǫ dt = ψ ψ ψ. () This model describes the coupling between waveguides in an optical array. In [], we examined how linear normal modes couple due to the cubic nonlinearity in (). Here instead, we assume a large amplitude localized initial condition. This is a natural and relevant consideration that parallels classical studies of discrete solitons in the nonlinear Schrödinger equation and light localization in nonlinear photonic structures. We assume that ψ j = r = ρ constant at a given node where ρ = ǫ u j (), and ψ k =, k j. The evolution of the excited node j is given by The solution of () is iω ǫ dψ j dt = r ψ j. () ψ j (t) = r e iǫr ω t = ρ ǫ eiρ 8ω t. 7

8 Thus, the solution u j can be approximated using () by where the nonlinear frequency is (( ) ) ρ u j (t) ρ cos 8ω +ω t, Ω ρ 8ω +ω. (4) This regime is valid when the correction to the frequency of oscillation due to the nonlinearity is smaller than the natural frequency 8ω ρ ω. (5) This means ω large enough. Discrete breathers for () in chains were studied by Panayotaros in [9]. He used a continuation argument in ǫ starting from the anticontinuous limit, to show the existence of discrete breathers. The nearest neighbors k of j solve a forced system iω ǫ dψ k dt N = kp ψ p = d k ψ k +ψ j + ψ p, (6) p= p k, p j where we neglected the cubic terms. A particular solution of the forced part iω dψ f k ǫ dt = ψ j is ψ f k (t) = ( e iǫr t) ω = ( ) ǫ e iρ 8ω t. (7) r ρ Similarly, the evolution of the next nearest neighbors l is given by iω ǫ dψ N l dt = lp ψ p = d l ψ l +ψ k + p= p l, p k We have the proportional scalings for ψ k and ψ l ( ) ( ) ǫ ǫ ǫ ψ k = O, ψ l = O ρ ρ, corresponding to the scalings for u k and u l ( ) ( ) ǫ ǫ u k = O, u l = O ρ ρ, similarly to the scalings for ω = (8,9). ψ p. (8) 8

9 As above, for simplicity and without loss of generality, we assume an excitation of node j =. The evolution of the nodes {,...,N} is described by the forced system of linear ordinary differential equations dϕ dt = iǫ ω ( ϕ+f ), (9) where ϕ = (ψ,ψ,...,ψ N ) T, and where f = (f,f,...,f N ) T is the forcing term such that f k = ψ if k adjacent to (k ) and otherwise. We expand ϕ using a basis of the eigenvectors z m of ϕ = N m= β m z m. () Substituting () into (9) and projecting on each eigenvector z m, we get dβ m dt = i ǫ ω ω m β m i ǫ ω ψ zk m. () The harmonic frequencies of ϕ are ǫ ω ω m for m {,...,N } such that k zm k. The harmonic frequencies of v = (u,u,...,u N ) T using () are ǫ ω ω m +ω, () for m {,...,N } such that k zm k. Notice that () and () are almost equal for large ω and small ǫ, and these are the conditions of validity of the approximation by modulation theory. k 4 Numerical results We illustrate our findings on two graphs: a cycle joined to a single isolated node known as the paw graph (Fig.), and the cycle 6 (Fig.7). The first graph is not regular and has one symmetry, the permutation of nodes and 4. The second is a cycle invariant under cyclic permutations. We will see how symmetries affect the observed modes and how localized solutions destabilize. The system of ordinary differential equations () is solved in double precision, using a Runge-Kutta 4-5 method with a time step and a relative error of 8. To check the validity of the solutions, we calculated the Fourier transform û k of each u k, k {,...,N}. This revealed the frequencies of the motion and allowed a detailed comparison with the analysis of sections and. In practise, we used the fast Fourier transform (FFT) of Matlab on a time-series of n = points on a time t f = to approximate the continuum Fourier transform. The data was multiplied by a Hamming window u k (m) (.54.46cos 9 ( π (m ) n )),

10 for m {,...,n}, k {,...,N}. 4. Paw graph We consider the paw graph studied in [], shown in Fig.. 4 Figure : Paw graph. This graph is invariant by the permutation symmetry of nodesand4. Then the components u (t) and u 4 (t) of the solution u(t) are equal. 4.. Natural frequency ω = Exciting node We solve equation () with ω =, ǫ =. for the paw graph. The left panel of Fig. shows the time evolution of the solutions u k, k {,...,4} when exciting the system at node j = with initial amplitude u () = ρ =. We see clearly a localized solution; it can be observed over more than a thousand periods with no significant decay. The logarithm with base of the modulus of the discrete Fourier transform of the solutions log ( û k ), k {,...,4} are shown on the right panel of Fig.. The Fourier components of the neighbor u corresponding to the linear mode and the nonlinear mode are about equal. As we go to the next nearest neighbor the Fourier component due to the nonlinear excitation of node is a times smaller than the linear response of the network. This is a general feature that we see on all the systems we have analyzed. It confirms the exponential localization of the nonlinear mode. From the Fourier spectrum, we determine that u oscillates at the nonlinear frequency (5) Ω = 6π π Γ ( 4 ).54 and at the odd harmonics of Ω (Ω and weakly at 5Ω ) due to the Fourier expansion of the solution () (formula (54)

11 Appendix.A) where u (t) = ρ cn ( ) ρt, 4 Ω [b cos(ω t)+ + b cos(ω t)+b 5 cos(5ω t)+...], π 5π +e 5π. b = e π +e π, b = e +e π, b 5 = e ( ǫ The solutions shown in Fig. are such thatu = O(ρ) and u = O ρ ( ǫ u 4 = O ρ ). To describe the evolution of u,u and u 4, it is then natural to reduce the system () to the linear system forced by u d dt u u = ǫ u 4 ), u = d u dt = u, () u u + ǫu (4) u 4 The Fourier representation shows that the nearest neighbor u and similarly the next nearest neighbors u and u 4 oscillate at the nonlinear frequency Ω and at the eigenfrequencies () of the matrix ǫ ǫω =. ǫω =. ( )., ( + ).86, where =. (5) z = The eigenvectors of are 6, z =, z = 6+ + (6). The absence of the eigenfrequency ǫω =..77 is due to zk = z =. k

12 The equations () are then d a ( dt =. ) a +.( ) u, (7) 6 d a dt =. ( + ) a +.( +) 6+ u, (8) To conclude, when exciting node with a large amplitude ρ, the evolution of u is given by () and the evolution of u, u and u 4 by u (t) u (t) = a (t)z +a (t)z, (9) u 4 (t) wherea anda are solutions of equations (7,8). Note that solving the reduced system (,4) yields the same results as the ones shown in Fig. u j - 4 log FFT(u j ) t -6 frequency Figure : Solution of equation () for the paw graph and an initial condition u () =, u () = u () = u 4 () =. Left panel: time evolution of u (red online), u (blue online), u = u 4 (black online). Right panel: Fourier transform of the solutions û (red online), û (blue online), û = û 4 (black online). The parameters are ω =, ǫ =.. Exciting node To observe the different response of the system, we now excite node j = with initial amplitude u () = ρ = and ω =, ǫ =.. The left panel of Fig. shows the time evolution of the solutions u k, k {,...,4}. The logarithm with base of the modulus of the discrete Fourier transform of the solutions log ( û k ), k {,...,4} are shown on the right panel of Fig..

13 u j - 4 log FFT(u j ) t -4 4 frequency Figure : Solution of equation () for the paw graph and an initial condition u () =, u () = u () = u 4 () =. Left panel: time evolution of u (blue online), u = u = u 4 (black online). Right panel: Fourier transform of the solutions û (blue online), û = û = û 4 (black online). The parameters are ω =, ǫ =.. The Fourier spectrum shows that u oscillates at the nonlinear frequency (5) Ω = 6π π Γ ( 4 ).54 and at the odd harmonics of Ω (Ω and weakly at 5Ω ). The solutions shown in Fig. are such that u = O(ρ) and u k ρ for k =, and 4. To describe the evolution of u,u and u 4, it is then natural to reduce the system () to the linear system forced by u d dt u u = ǫ u u +ǫu. (4) u 4 u 4 The Fourier representation shows that the nearest neighbors u, u and u 4 oscillate at the nonlinear frequency Ω and at the eigenfrequencies () of the matrix ǫ ǫω = ǫω =..447, (4) where =. (4) The eigenvectors of are z =, z =, z =.

14 Here, only one linear frequency exists. The absence of the eigenfrequency ǫω =. is due to k,k< z k + k,k> z k = z +z +z =. Then, the evolution of the nearest neighbors u, u and u 4 is given by u (t) a (t) u (t) = a (t)z +a (t)z = a (t), (4) u 4 (t) a (t) where a and a are solutions of equations d a dt = ǫω a + ǫu, (44) d a dt = ǫω a +ǫu. (45) This explains why the solution u is equal to u and u 4. Exciting node We now excite node with initial amplitude u () = ρ = and ω =, ǫ =.. The left panel of Fig.4 shows the time evolution of the solutions u k, k {,...,4}. The logarithm with base of the modulus of the discrete Fourier transform of the solutions log ( û k ), k {,...,4} are shown on the right panel of Fig.4. u j t log FFT(u j ) frequency Figure 4: Solution of equation () for the paw graph and an initial condition u () =, u () = u () = u 4 () =. Left panel: time evolution of u (red online), u (blue online), u (green online) and u 4 (black online). Right panel: Fourier transform of the solutions û (red online), û (blue online), û (green online) and û 4 (black online). The parameters are ω =, ǫ =.. 4

15 From the Fourier spectrum, we can determine that u oscillates at the nonlinear frequency (5) Ω = 6π π Γ ( 4 ).54 and at the odd harmonics of Ω (Ω and weakly at 5Ω ). The solutions shown in Fig.4 are such that u = O(ρ) and u k ρ for k =, and 4. To describe the evolution of u,u and u 4, it is then natural to reduce the system () to the linear system forced by u d dt u u = ǫ u u +ǫu (46) u 4 The Fourier representation shows that the nearest neighbors u and u 4 (and similarly the next nearest neighbor u ) oscillate at the nonlinear frequency Ω and at the eigenfrequencies () of the matrix ǫ where u 4 ( ( )) π ǫω =. cos.59, 9 ( ( )) 4π ǫω =. cos.5749, 9 ( ( π ǫω =. +cos )) =. (47) The eigenvectors of are z = , z =.449.9, z = Natural frequency ω We now analyze a non zero natural frequency, chooseω = andǫ =. and solve the graph nonlinear wave equation () for the paw graph. The left of Fig.5 shows the time evolution of the localized solutions at j = with amplitude u () =. Again this has been observed for over a thousand periods with no significant decay. The right of Fig.5 shows the logarithm with base of the modulus of the discrete Fourier transform of the solutions log ( û k ), k {,...,4}. Note that u oscillates at frequencies Ω, Ω and weakly at 5Ω where Ω 4. The nearest neighbor u (and similarly the next nearest neighbors u and u 4 ) oscillate at the nonlinear frequency Ω and at eigenfrequencies () which are 5

16 almost equal to those in () ǫ ǫω +ω ω ω +ω, ǫ ǫω +ω ω ω +ω.. We proceed to validate the modulation theory by solving the graph nonlinear Schrödinger equation () with ω = and ǫ =.. From this solution, we calculate u using the change of variables (). The comparison of the left panels of Figs. 5 and 6 confirm the approximation by the modulation theory. u j t log FFT(u j ) frequency Figure 5: Solution of equation () for the paw graph and an initial condition u () =, u () = u () = u 4 () =. Left panel: time evolution of u (red online), u (blue online), u = u 4 (black online). Right panel: Fourier transform of the solutions û (red online), û (blue online), û = û 4 (black online). The parameters are ω =, ǫ =.. 6

17 u j - 4 log FFT(u j ) t frequency Figure 6: Solution of equation () for the paw graph plotted in the u variables using () for an initial condition u () =, u () = u () = u 4 () =. Left panel: time evolution ofu (red online), u (blue online), u = u 4 (black online). Right panel: Fourier transform of the solutions û (red online), û (blue online), û = û 4 (black online). Same parameters as in Fig.5. When exciting initially node, we note that u oscillates at frequencies Ω, Ω and weakly at 5Ω where Ω 4. The nearest neighbors u, u and u 4 oscillate at the nonlinear frequency Ω and at eigenfrequencies () ǫω +ω = ǫω +ω.. Similarly, exciting initially node, we observe that u oscillates at frequencies Ω, Ω and weakly at 5Ω where Ω 4. The nearest neighbors u and u 4 (and similarly the next nearest neighbor u ) oscillate at the nonlinear frequency Ω and at eigenfrequencies () ǫω +ω.56, ǫω +ω.546, ǫω +ω.66. 7

18 4. Cycle Figure 7: Cycle Natural frequency ω = We now consider the cycle 6 where all the nodes are invariant under cyclic permutations, so that they behave the same way. We solve equation () with ω = and ǫ =. when exciting the system at site j = with initial amplitude u () =. The left panel of Fig.8 shows the time evolution of the solutions u k, k {,...,6}. Notice how u is large while u, u and u 4 are small, indicating a localized oscillation. The right panel of Fig.8 shows the logarithm with base of the modulus log ( û k ), k {,...,6} of the discrete Fourier transform of the solutions. The permutation symmetry of nodes 6 and 5 is reflected in the solutions u = u 6 and u = u 5. This network should show five linear modes, nevertheless due to the symmetry only three linear modes are present. The dynamics at nodes {,...,6} is described by the linear system u u ǫu d u dt u 4 u 5 = ǫ u u 4 u 5 + u 6 u 6 ǫu (48) Using the result of the Fourier spectrum, we can determine that u oscillates at frequencies Ω.54, Ω and weakly at 5Ω. The nearest neighbors (u and u 6 ) and next nearest neighbors (u and u 5 ) oscillate at frequencies Ω and at the eigenfrequencies () of the matrix ǫ 8

19 ǫω =. ǫω =..6, ǫω5 =. ( )., ( + ).86, where The eigenvectors of are z =, z = z 4 = =. (49), z5 =, z =. The absence of the frequencies ǫω and ǫω 4 is due to zl = z +z 5 =, l zl 4 = z 4 +z5 4 = l, 9

20 u j t log FFT(u j ) frequency Figure 8: Left panel: time evolution of u (red online), u (blue online), u (green online) and u 4 (black online) solutions of () with ω =, ǫ =. in cycle 6 for initial amplitudes u () =, u () = u () = u 4 () = u 5 () = u 6 () =. Right panel: Fourier transform of the solutions û (red online), û (blue online), û (green online) and û 4 (black online). 4.. Natural frequency ω We consider the equation () with ω = and ǫ =. for cycle 6 when exciting the system at site j = with initial amplitude u () =. The left panel of Fig.9 shows the time evolution of the localized solutions at j = with amplitude u () =. Again this has been observed for over a thousand periods with no significant decay. The right panel of Fig.9 shows the logarithm with base of the modulus of the discrete Fourier transform of the solutions log ( û k ), k {,...,6}. Note that u oscillates at frequencies Ω, Ω and weakly at 5Ω where Ω 4. The nearest neighbors and next nearest neighbors oscillate at frequencies Ω and at the eigenfrequencies () of the matrix ǫ ω I ǫω +ω, ǫω +ω.6, ǫω 5 +ω..

21 u j t log FFT(u j ) frequency 4 Figure 9: Left panel: time evolution of u (red online), u (blue online), u (green online) and u 4 (black online) solutions of () with ω =, ǫ =. in cycle 6 for initial amplitudes u () =, u () = u () = u 4 () = u 5 () = u 6 () =. Right panel: Fourier transform of the solutions û (red online), û (blue online), û (green online) and û 4 (black online). Solving the graph nonlinear Schrödinger equation () with ω = and ǫ =. for cycle 6 and calculating u using the change of variables (), yields the same results as the ones shown in Fig Weighted graphs We now assume that the graph has weights on its edges. Consider the paw graph studied in section 4. with weights α, α, α 4, α 4 on all the edges e,, e,, e,4, e,4 respectively. This is done to show the effect of a symmetry break in this graph. We will see that our formalism extends to this case. α α α 4 α 4 4 Figure : A weighted paw graph.

22 The Laplacian of the weighted paw graph α α = α (α +α +α 4 ) α α 4 α (α +α 4 ) α 4 (5) α 4 α 4 (α 4 +α 4 ) The reduced matrix is = (α +α +α 4 ) α α 4 α (α +α 4 ) α 4 (5) α 4 α 4 (α 4 +α 4 ) We solve equation () with ω =, ǫ =. for the paw graph Fig. with coupling coefficients α, =, α, =, α,4 = 4, α,4 =. The left panel of Fig. shows the time evolution of the solutionsu k, k {,...,4} when exciting the system at node j = with initial amplitude u () = ρ =. The right panel of Fig. shows the logarithm with base of the modulus log ( û k ), k {,..., 4} of the discrete Fourier transform of the solutions. u j - 4 log FFT(u j ) t frequency Figure : Solution of equation () for the paw graph Fig. with α, =, α, =, α,4 = 4, α,4 = and an initial condition u () =, u () = u () = u 4 () =. Left panel: time evolution of u (red online), u (blue online), u (green online) and u 4 (black online). Right panel: Fourier transform of the solutions û (red online), û (blue online), û (green online) and û 4 (black online). The parameters are ω =, ǫ =.. Fig. shows the exponential decay of the Fourier amplitude of the nonlinear component as we go from the excitation node to its nearest neighbor and next nearest neighbor, even though the coupling between and is larger than the one between and. The weight between nodes and is larger than

23 the one between nodes and 4 so that the symmetry between and 4 is broken. Comparing to section 4., here we observe all the eigenfrequencies of ǫω.69, ǫω.6, ǫω.89. (5) The eigenvectors of are z = , z = , z = (5) Localization vs delocalization Up to now, we choose a large amplitude ρ and a small couplingǫ. This leads to a localized solution. For a fixed amplitude, as we increase the coupling, the linear spectrum of the matrix j will collide with the nonlinear frequency Ω. For lattices, this causes the disappearance of the localized nonlinear solutions, see [] and [4]. Aubry and Mackay [] point out that localized solutions also exist for general networks, as long as there is no resonance with the linear spectrum. This resonance argument explains why the threshold of the localization depends on the maximal normal eigenfrequency. Numerically, we observe that the localized solution disappears and there is a strong coupling with the neighboring nodes. Also, the spectrum does not show well defined frequencies. To illustrate the delocalization regime, we choose ω =, ǫ =.5 and solve equation () for the paw graph with initial amplitude u () =. Fig. shows the time evolution of the solutions, there is a strong exchange of energy between nodes. 4 u j Figure : Time evolution of u (red online), u (blue online), u and u 4 (black online) solutions of () with ω =, ǫ =.5 in the paw graph for initial amplitudes u () =, u () = u () = u 4 () =. t Using the localized character of the solution and the Fourier spectrum as

24 indicators, we examined the parameter plane (ǫ, ρ) and plotted the regions of localization versus delocalization. First we consider the paw graph and plot these regions for initial excitations of nodes, and. This is shown in Fig.. 8 localization amplitude ρ 6 4 node node node cycle6 delocalization Figure : Regions of delocalization in the (ǫ,ρ) plane with ω =, exciting initially nodes, and for the paw graph and any node for cycle 6 with amplitude ρ. The lines separates the regions of localization (large ρ) and the regions of delocalization (small ρ). ǫ The separation curves behave as O( ǫ) since the eigenfrequencies of ǫ j scale like ǫ. As one can see, the amplitudes for localization decrease as one excites node then node and finally node. The maximal linear frequencies are ǫω N where ω N =.9,.96 and, respectively for nodes, and. Then the linear spectrum is closer to the nonlinear frequency Ω for nodes and and farther for node. We observe that the curve for node is above the one for node in Fig. despite the fact that the maximal linear frequencies are very close. This could be due to their different degrees, d = and d =. Similarly, the curve for node (d = ) of the paw graph is above the one for any node in cycle 6 (d = ), see Fig., even though the maximal eigenfrequencies are equal. For equal degrees, exciting node ( d = ) and any node of cycle 6 gives very close curves in Fig.. Taking ω leads to very close results as the linear spectrum and the nonlinear frequency are both approximately shifted by ω. When choosing a quintic nonlinearity in equation (), the nonlinear frequency is much larger while the linear spectrum is unchanged. Then, we need to increase ǫ considerably to observe coupling to the linear modes. 4

25 5 Conclusion We studied localized solutions for a nonlinear graph wave equation. These are approximated by a nonlinear ordinary differential equation for the excited node and a forced linear system for the neighboring nodes. We validate this approximation by calculating the Fourier spectrum of the numerical solution. This shows the nonlinear frequency of the excited node together with the normal eigenfrequencies of the linear system describing the neighboring nodes. The existence of these localized solutions is confirmed using numerical simulations and modulation theory. We also examined the localization / delocalization regions in the parameter plane (ǫ, ρ). Localization holds when the nonlinear frequency Ω is well above the linear spectrum which is bounded by O( ǫ ω N ), where ω N is the maximal normal eigenfrequency of j. This condition explains the observed localization amplitude ρ. Acknowledgment This work is part of the XTerM project, co-financed by the European Union with the European Regional Development Fund (ERDF) and by the Normandie Regional Council. A.B. Aceves thanks the Laboratoire de Mathématiques de l INSA de Rouen Normandie for its hospitality during a visit in 6. We acknowledge the support of Agence Nationale de la Recherche through the project Fractal Grid. Appendices A Basic properties of cosine Jacobi elliptic functions Consider the upper limit φ of the integral x = as a function of x. The function φ dy κ sin y cn(x,κ) = cos(φ), 5

26 is called the cosine Jacobi elliptic functions [8], with elliptic modulus κ [, ]. cn(x,κ) is periodic function with period 4K, where K K(κ) = π dθ κ sin(θ), ( ) which for κ = gives K = Γ ( 4) 4 π.854. The Fourier series of the cosine elliptic function is given by cn(x,κ) = π κk m= q ( m+ +q m+ cos where q = e πk K so that q = e π for κ =. Γ ( 4 ) ρ π (m+) πx K The Fourier series of the periodic solution () with period T = 4K is given by ( ) u j (t) = ρ cn ρt, = 8π T where Ω = π T and ). ( ) ρ = b m+ cos((m+)ω t), (54) m= b m+ = e π(m+ ) +e π(m+). B Period of the Duffing oscillator The period of u j solution of () d u j dt = u j, dropping the j index for clarity and integrating the equation, we get ( ) du = ( ρ 4 u 4), dt where we assumed u() = ρ and du dt () =. This yields dt = ± du ρ4 u 4. The period T of u j (t) is then given by the elliptic integral T = 4 ρ du ρ4 u 4 = 4 ρ 6 dz z 4.

27 We know that dz = Γ ( 4 ) z 4 4 π, where Γ(.) is the gamma function and Γ( 4 ).656. Then, the period of oscillations is T = Γ ( 4 ) ρ π. C Derivation of the Graph nonlinear Schrödinger equation We consider the continuum equation of () To eliminate the term ω u, we write d u dt +ω u = ǫ u u. (55) u = ǫψ(t)e iωt + ǫψ (T)e iωt, where T = ǫt and ψ is the complex conjugate of ψ. We have d ( u d dt = ψ ǫ5 dt eiωt + d ψ ) ( ) dt e iωt +iωǫ dψ dt eiωt dψ dt e iωt ω ǫ ( ψe iωt +ψ e iωt). The left hand side of the equation (55) is d ( u d ψ dt +ω u = ǫ 5 dt eiωt + d ψ ) dt e iωt +iωǫ ( ) dψ dt eiωt dψ dt e iωt. The right hand side of the equation (55) gives ǫ u u ( = ǫ ψe iωt + ψ e iωt ψ e iωt (ψ ) e iωt We obtain for the order O(ǫ ) ψ ψe iωt ψ ψ e iωt). iω dψ dt eiωt iω dψ dt e iωt = ψe iωt + ψ e iωt ψ e iωt (ψ ) e iωt ψ ψe iωt ψ ψ e iωt. 7

28 Multiplying by e iωt, we get iω dψ dt iωdψ dt e iωt = ψ + ψ e iωt ψ e iωt (ψ ) e 4iωt ψ ψ ψ ψ e iωt. The terms with a non zero phase are rotating fast and average to zero on the slow time scale. Only the terms that have phase contribute. This is the rotating wave approximation [6]. We obtain the nonlinear Schrödinger equation References iω dψ dt = ψ ψ ψ. [] A. J. Sievers and S. Takeno, Intrinsic localized modes in anharmonic crystals, Phys. Rev. Lett. 6, (988). [] J. B. Page, Asymptotic solutions for localized vibrational modes in strongly anharmonic periodic systems, Phys. Rev. B 4, (99). [] R. S. MacKay and S. Aubry, Proof of existence of breathers for timereversible or Hamiltonian networks of weakly coupled oscillators, Nonlinearity 7, 6-64 (994). [4] S. Flach, Conditions on the existence of localized excitations in nonlinear discrete systems, Phys. Rev. E 5, 4 (994). [5] S. Flach and C. R Willis, Discrete breathers, Physics reports 95, 8-64 (998); S. Flach and A. V. Gorbach, Discrete breathers - Advances in theory and applications, Physics Reports 467, -6 (8). [6] P. Binder, D. Abraimov, A. V. Ustinov, S. Flach and Y. Zolotaryuk, Observation of breathers in Josephson ladders, Phys. Rev. Lett. 84, 745 (). [7] A. B. Aceves, C. De Angelis, T. Peschel, R. Muschall, F. Lederer, S. Trillo, and S. Wabnitz, Discrete self-trapping, soliton interactions, and beam steering in nonlinear waveguide arrays, Phys. Rev. E 5, 7 (996). [8] H. S. Eisenberg, Y. Silberberg, R. Morandotti, A. R. Boyd, and J. S. Aitchison, Discrete spatial optical solitons in waveguide arrays, Phys. Rev. Lett. 8, 8-86 (998). [9] J. W. Fleischer, M. Segev, N. K. Efremidis, and D. N. Christodoulides, Observation of two-dimensional discrete solitons in optically induced nonlinear photonic lattices, Nature (London) 4, 47 (). 8

29 [] M. Peyrard, Nonlinear dynamics and statistical physics of DNA, Nonlinearity 7 R (4). [] M. Sato, B. E. Hubbard, A. J. Sievers, B. Ilic, D. A. Czaplewski, and H. G. Craighead, Observation of locked intrinsic localized vibrational modes in a micromechanical oscillator array, Phys. Rev. Lett. 9, 44 (). [] M. Sato, S. Yasui, M. Kimura, T. Hikihara and A. J. Sievers, Management of localized energy in discrete nonlinear transmission lines, Europhysics Letters 8, (7). [] D. Christodoulides, Intrinsic Localized Modes in Optical Photonic Lattices and Arrays, APS Meeting Abstracts (6). [4] Y. Shen, P. G. Kevrekidis, G. Srinivasan and A. B. Aceves, Existence, stability and dynamics of discrete solitary waves in a binary waveguide array, J. Phys. A: Math. Theor. 49, 955 (6). [5] M. Jenkinson and M.I. Weinstein, Discrete Solitary Waves in Systems with Nonlocal Interactions and the Peierls-Nabarro Barrier, Communications in Mathematical Physics Volume 5, Issue, (7). [6] R.A. Vicencio and M. Johansson, Discrete flat-band solitons in the kagome lattice, Physical Review A 87, 68 (). [7] Y. Kou and J. Förstner, Discrete plasmonic solitons in graphene-coated nanowire arrays, Optics Express 4, Issue 5, (6); Y. Fan, B. Wang, K. Wang, H. Long and P. Lu. Plasmonic Zener tunneling in binary graphene sheet arrays, Optics Letters 4, Issue, (6). [8] S. Flach and A. Gorbach, Discrete breathers in Fermi-Pasta-Ulam lattices, Chaos 5, 5 (5). [9] P. Panayotaros, Continuation of normal modes in finite NLS lattices, Physics Letters A, 74, 9, (). [] A. A. Ovchinnikov, N.S. Erikhman, K.A. Pronin, Vibrational-Rotational Excitations in Nonlinear Molecular Systems, Kluwer Academic/Plenum Publishers, (). [] A. C. Scott, Nonlinear Science: Emergence and Dynamics of Coherent Structures, Oxford Texts in Applied and Engineering Mathematics, nd edn, Oxford-New York: Oxford University Press (). [] C. Maas, Transportation in graphs and the admittance spectrum, Discrete Applied Mathematics, 6, -49, (987). [] D. Cvetkovic, P. Rowlinson and S. Simic, An Introduction to the Theory of Graph Spectra, London Mathematical Society Student Texts 75, Cambridge: Cambridge University Press (). 9

30 [4] J-G. Caputo, A. Knippel and E. Simo, Oscillations of networks: the role of soft nodes, J. Phys. A: Math. Theor. 46, 5 (). [5] J-G. Caputo, I. Khames, A. Knippel and P. Panayotaros, Periodic orbits in nonlinear wave equations on networks, J. Phys. A: Math. Theor. 5, 75 (7). [6] A. C. Scott, Encyclopedia of nonlinear science, London: Routledge, Taylor and Francis Group (5). [7] A. Grolet, N. Hoffmann, F. Thouverez C. Schwingshackl, Travelling and standing envelope solitons in discrete non-linear cyclic structures. Mechanical Systems and Signal Processing, Vol 8, (6). [8] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions. New York: Dover (965). [9] Steve Fisk, A very short proof of Cauchy s interlace theorem for eigenvalues of Hermitian matrices. (5). [] J-G. Caputo, I. Khames and A. Knippel. On graph Laplacians eigenvectors with components in {,,}. "Discrete and Applied Mathematics, in revision, 8. [] Y.S. Kivshar and M. Peyrard, Modulational instabilities in discrete lattices, Phys. Rev. A (99); Y.S. Kivshar, Localized modes in a chain with nonlinear on-site potential, Physics Letters A, volume 7, Issue, 7-78 (99). [] A. B. Aceves and J-G. Caputo, Mode dynamics in nonuniform waveguide arrays: A graph Laplacian approach, Journal of Optics, Volume 6, Issue, article id. 5 (4).