arxiv: v2 [nlin.ps] 26 Sep 2018

Transcription

1 Dynamics of Dirac solitons in networks arxiv:7.577v [nlin.ps] 6 Sep 8 K.K. Sabirov a,b, D.B. Babajanov a, D.U. Matrasulov a and P.G. Kevrekidis c a Turin Polytechnic University in Tashkent, 7 Niyazov Str., 95, Tashkent, Uzbekistan b Tashkent University of Information Technologies, Amir Temur Str.,, Tashkent, Uzbekistan c Department of Mathematics and Statistics, University of Massachusetts, Amherst, Massachusetts USA d.b.babajanov@gmail.com August 8 Abstract. We study dynamics of Dirac solitons in prototypical networks modeling them by the nonlinear Dirac equation on metric graphs. Stationary soliton solutions of the nonlinear Dirac equation on simple metric graphs are obtained. It is shown that these solutions provide reflectionless verte transmission of the Dirac solitons under suitable conditions. The constraints for bond nonlinearity coefficients, conjectured to represent necessary conditions for allowing reflectionless transmission over a Y-junction are derived. The Y-junction considerations are also generalized to a tree network. The analytical results are confirmed by direct numerical simulations. PACS numbers:

2 Dynamics of Dirac solitons in networks. Introduction Nonlinear evolution equations on networks have attracted much attention recently [,, 3, 4, 5, 6, 7, 8, 9, ]. Such interest is caused by a broad variety of potential applications of the nonlinear wave dynamics and soliton transport in networks, such as Bose-Einstein condensates (BECs) in branched traps, Josephson junction networks, the DNA double heli, polymer chains, etc. Despite the rapidly growing interest in wave dynamics on networks, most of the studies are mainly focused on nonrelativistic wave equations such as the nonlinear Schrödinger (NLS) equation [,, 3, 4, 5, 6, 7, ]. Nevertheless, there is a number of works on the sine-gordon (sg) equation in branched systems [8, 9, ]. However, relativistic wave equations such as the nonlinear Klein-Gordon and Dirac equations are important in field theory and condensed matter physics and hence eploring them on metric graphs is of interest in its own right. These graphs consist of a system of bonds which are connected at one or more vertices (branching points). The connection rule is associated with the topology of a graph. When the bonds can be assigned a length, the graph is called a metric graph. In this paper we address the problem of the nonlinear Dirac equation on simple metric graphs by focusing on conservation laws and soliton transmission at the graph vertices. Our prototypical eample will be the Y-junction. Early studies of the nonlinear Dirac equation date back to Thirring [] and Gross-Neveu [3] models of field theory. Integrability, the nonrelativistic limit and eact solutions of the relevant models have been considered in, e.g., [4, 5, 6, 7, 8] among others. A potential application of the nonlinear Dirac equation to neutrino oscillations was discussed in [9]. Recently, the possibility of eperimental realization of Dirac solitons in Bose-Einstein condensates in honeycomb optical lattices was discussed in [] where the nonlinear Dirac equation was derived. Further analysis of this setting in Refs. [, 3, 4, 6] has ecited a rapidly growing interest in the nonlinear Dirac equation and its soliton solutions (see, e.g. [7]- [49]). In [7, 8, 3, 3, 33, 34], a detailed study of soliton solutions for different types of nonlinearity, their stability and discussions of conservation laws were presented. These developments have, in turn, had an impact also on the mathematical literature where the stability of solutions for different forms of the nonlinear Dirac equation was eplored in [9, 3, 36]. In [37]-[4], the stability of Gross-Neveu solitons in both d and 3d

3 Dynamics of Dirac solitons in networks 3 and for other nonlinear Dirac models an important energy-based stability criterion was developed. The numerical corroboration of the stability for solitary waves and vortices in D nonlinear Dirac models of the Gross-Neveu type was considered quite recently in [35]. Dirac solitons in branched systems can, in principle, be eperimentally realized in different systems of optics, but also envisioned elsewhere (e.g. in atomic physics etc.). A relevant such possibility consists of the discrete waveguide arrays of [44, 49] that can be formulated as a branched system to be described (in the appropriate long wavelength limit [44]) by a nonlinear Dirac equation on metric graphs. Such networks in optical systems have been proposed, e.g., earlier in [5]. More eotic possibilities can be envisioned in branched networks of honeycomb (i.e., armchair nanoribbon) optical lattices for atomic BECs [3], although we will not focus on the latter here. b b O b 3 Figure. Metric star graph In the present work, we focus on a prototypical eample of eact solutions and transport of Dirac solitons through the network vertices. In particular, we show that under certain constraints, eact soliton solutions of the nonlinear Dirac equation on simple metric graphs can be obtained. The identified soliton solutions provide reflectionless transmission of the Dirac solitons at the vertices. This renders possible the tuning of the transport properties of a network in such a way that it canprovide ballistic transport of the Dirac solitons. This paper is organized as follows: In the net section we give the formulation of the problem on a metric star graph which includes the derivation

4 Dynamics of Dirac solitons in networks 4 of the verte boundary conditions. Section III presents soliton solutions of the nonlinear Dirac equation on a metric star graph for fied and moving frames. Also, an analysis of the verte transmission of the Dirac solitons on the basis of the numerical solution of the nonlinear Dirac equation is presented. In section IV we etend the treatment to a metric tree graph and discuss the etension to other simple topologies. Finally, section V presents some concluding remarks and a number of important directions for future work. These include a conjecture about the sufficient, yet not necessary nature of our verte conditions towards ensuring reflectionless transmission.. Conservation laws and verte boundary conditions The nonlinear Dirac system that we are going to eplore, i.e., specifically the Gross- Neveu model, follows from the Lagrangian (in the units h = c = ) [7, 8, 3] L = Ψγ i µ µ Ψ m ΨΨ+ g ), ( ΨΨ µ =,, () and describes interacting Dirac fields with and corresponding to time and coordinate variables, respectively. Here g is the nonlinearity coefficient characterizing the strength of the nonlinear interaction. The field equations for this Lagrangian lead to the nonlinear Dirac model of the form where and γ = (iγ µ µ m)ψ+g ( ΨΨ ) Ψ =, () Ψ(,t) = φ(,t) ; Ψ = Ψ γ = (φ (,t), χ (,t)), (3) χ(, t), γ =. For a metric star graph consisting of three semi-infinite bonds (see, e.g., Fig. ), the Lagrangian density of the nonlinear Dirac field for each bond can be written as L j where j parametrizes the bond and L for each is of the form of Eq. (). The field equation following from this Lagrangian density can be written as Eq. () before, where the spatial coordinates are defined as (,] and,3 [, ), while coincides with the graph verte. The formulation of an evolution set of equations on metric graphs requires imposing verte boundary conditions which provide gluing of the graph bonds at the graph

5 Dynamics of Dirac solitons in networks 5 vertices. For the linear Dirac equation on a metric graph studied earlier by Bolte and Harrison in [5], the general verte boundary conditions have been derived from the self-adjointness of the Dirac operator on a graph. In that case such conditions led to Kirchhoff rules and continuity of the wave function at the graph verte [5]. However, for the nonlinear problem in order to derive verte BCs, it is arguably more natural to use suitable conservation laws of the nonlinear flow, which give rise to appropriate BCs. For the nonlinear Dirac system fundamental conservation laws can be presented in terms of the momentum-energy tensor given by [7, 8] T j µν = i [ Ψ j γ µ ν Ψ j ν Ψj γ µ Ψ j ] g µν L j, where g µν =. The energy on each bond of a star graph can be written as E j = Td j = [ i b j b j [ Ψ j γ Ψ j Ψj γ Ψ j ] m Ψ j Ψ j + g j ( Ψ j Ψ j ) ]d,(4) where integration is performed along the bond b j. Here, to derive the verte boundary conditions (VBC) we use conservations of charge and energy. The total charge for the star graph is defined as Q = Q +Q +Q 3, where the charge for each bond is given by ( Q j = φj + χ j ) d. (5) b j The conservation of charge(kirchhoff rule for charge current), Q = yields the following verte boundary condition: Re[φ χ ] = = Re[φ χ ] = + Re[φ 3χ 3 ] =. (6) Here we used the asymptotic conditions Ψ at and Ψ,3 at. (7) We note that the boundary conditions (6) can be derived also from the conservation of the current density given by j = ϕχ +χϕ. The energy conservation, Ė = leads to Im[φ t χ +χ t φ ] = = Im[φ t χ +χ t φ ] = + Im[φ 3 t χ 3 +χ 3 t φ 3 ] =. (8)

6 Dynamics of Dirac solitons in networks 6 As per the above analysis, the VBC for NLDE on metric graphs have been obtained from the energy and charge conservation. However, the VBC given by Eqs. (6) can be fulfilled, if the following linear relations at the vertices are imposed (see Appendi for details): α φ = = α φ = +α 3 φ 3 =, α χ = = α χ = = α 3 χ 3 =, (9) where α,α,α 3 are the real constants which will be determined below. In the following we will use Eqs. (9) as the verte boundary conditions for Eq. () on a metric star graph. When applied to the eact soliton solutions of the NLDE, the VBC will lead to algebraic conditions connecting α j and g j, as will be derived in the net section. We note that in the linear limit (g j ) the verte boundary conditions given by Eq. (9) preserve the self-adjointness of the (linear) Dirac equation on metric graph, i.e. belong to the class of general boundary conditions derived by Bolte and Harrison in [5]. 3. Soliton dynamics and verte transmission In [7, 8, 3] the soliton solutions of the nonlinear Dirac equation on the line have been obtained in the form of standing wave solutions asymptoting to zero as ±. Here, we use the same definition for soliton solutions of the nonlinear Dirac equation on the metric star graph presented in Fig.. We note in passing that for soliton solutions of the NLDE derived in the Refs.[7, 8, 3] all elements of the energy-momentum tensor, ecept T become zero, i.e., T = T = T =. We look for the soliton solution of NLDE on a metric star graph in the form Ψ j (,t) = ψ j(,t) = e iωt A j(). () ψ j (,t) ib j () Then, from Eqs. () and () we have da j d +(m+ω)b j gj(a j Bj)B j =, () db j d +(m ω)a j gj (A j B j )A j =. () The verte boundary conditions for the functions A j and B j can be written as α A = = α A = +α 3 A 3 =, α B = = α B = = α 3 B 3 =. (3)

7 Dynamics of Dirac solitons in networks 7 Figure. (Color online) Charge density, Ψ j = Ψ Ψ = A j () + B j (), plotted using the standing wave soliton solutions given by Eqs. (4) and (5)( = ) for g =, g = g 3 =. A prototypical static solution of the system (), () vanishing at ± can be written as [8] A j () = (m+ω)cosh (β( )) β m+ωcosh(β( )) gj(m+ωcosh(β( ))), (4) B j () = (m ω)sinh (β( )) β m+ωcosh(β( )) gj(m+ωcosh(β( ))), (5) where is the position of the soliton s center of mass and β = m ω. In order for these solutions of Eqs. () and () to solve the problem on the metric graph, they need to also satisfy the verte boundary conditions (3). This can be achieved if the constants α j and coupling parameters g j fulfill the following conditions: α g = α g + α 3 g 3. (6) which stems from the first of Eqs. (9) and α,3 α = g g,3, (7) in accordance with the second one of conditions (9). The combination of the two leads to the sum rule : = +. (8) g g g3 It is important to note that this sum rule is derived by assuming that incoming wave comes from the first bond, b, while if it comes from the bond b (or b 3 ) one should replace g in Eq. (8) with g (or g 3 ). In Fig. plots of the soliton solutions of NLDE on

8 Dynamics of Dirac solitons in networks 8 metric star graph corresponding to Eqs. (4) and (5), satisfying the verte conditions are given. We note that soliton solutions given by Eqs.(4) and (5) describe the standing waves in metric graphs. The traveling wave (soliton) solutions of NLDE can be obtained by considering the case of the moving frame and taking into account the Lorentz invariance [8], in the case of the homogeneous domain (i.e., without a Y-junction). Assuming implicitly that the solitary wave is centered around a position well to the left of the origin, the Lorentz transformations between the frames moving with relative velocity, v can be written as [8] = γ( vt); t = γ(t v), (9) where γ = v = coshη, sinhη = v v. () Using these transformations, the traveling wave (soliton) solutions of the nonlinear Dirac equation in the moving frame determined by the constraints (6) can be written as ψ j (,t) = [cosh(η/)a j ( )+isinh(η/)b j ( )]e iωt, Ψ j (,t) = [sinh(η/)a j ( )+icosh(η/)b j ( )]e iωt () where, t and η are given by Eqs. (9) and (). It is important to note that this transformation does not affect the scaling dependence on the coupling parameters g j. For that reason, we epect the form of Eqs. (6) to still be valid in the case of traveling/moving wave solutions. To analyze the dynamics of the Dirac solitons on a metric graph, we solve numerically the time-dependent NLDE given by Eq. () for the verte boundary conditions given by Eq. (9) by considering the an initially traveling waveform localized in bond, far from the verte. We do this for the case when the constraint given by Eq.(8)isfulfilled, aswellasforthecaseofarbitraryg j whichdonotrespecttherelevant constraints. Fig. 3(a) presents the plots of the solution of Eqs. () and (9) obtained numerically for the values of g j obeying the constraint (8). For verte boundary conditions given by Eq.(3) parameters α j are chosen as α = / g, α 3 = / g 3, α = g (α / g +α 3 / g 3 ). One can observe the absence of the verte reflection in this case. Fig. 3(b) presents soliton solutions of nonlinear Dirac equation for the boundary

9 Dynamics of Dirac solitons in networks 9 (a) bond t bond t bond t (b) 5 bond bond t bond 3 t t. Figure 3. (Color online) Time and coordinate dependence of the charge density obtained from the numerical solution of Eq. () ( = ) for boundary conditions (9) for the case (a) when constraints (8) are fulfilled (g =, g = g 3 = ), (b) when these constraints are not fulfilled (g =, g = ;g 3 = 3.5). conditions (9) obtained numerically for those values of g j which do not fulfill Eq. (8). Here, reflection can be clearly discerned. Fig. 4(a) shows the conservation of the energy and reflectionless transmission through the Y-junction during the propagation of Dirac soliton in graph. The emergence of the verte reflection when the constraints (8) are not fulfilled, can also be seen from the Fig. 4(b). Here the verte reflection coefficient, which is determined according to the definition R = E (t = 7)/E(t = 7) (with E(t = 7) = E (t = 7) + E (t = 7) + E 3 (t = 7)) is plotted as a function of g 3 for fied values of g and g. This systematic analysis clearly illustrates the necessity of the symmetric scenario of g 3 =, consonant with our verte sum rule, for reflection to be absent. A similar phenomenon was observed in the case of other nonlinear PDEs on metric graphs, such as the NLS and sg equations, considered earlier in the Refs.[] and [], respectively. For the initial condition corresponding to a traveling wave of

10 Dynamics of Dirac solitons in networks Figure 4. (Color online) Numerical solution of Eq.() for boundary conditions (9). (a) Total energy and time-dependence of the energy for each bond with g =, g = g 3 =. (b) Reflection coefficient as a function of g 3 for g =, g =. the NLDE in the form of () centered at, with initial velocity v > for t +, our numerical computations provide strong evidence to the following fact. The solution can be asymptotically represented as a superposition of solitary waves in the branches associated with the second and third bonds, only if the condition (9) holds. In other words, under the necessary constraints (8), Eqs. (9) asymptotically lead to refectionless transmission through the verte. Figure 5. Tree graph. (,], i [,L i ], and ij [,+ ) with i =,;j =,,3.

11 Dynamics of Dirac solitons in networks 4. Etension for other graphs The above treatment of the NLDE on the metric star graph can be etended to the tree graph presented in Fig. 5. It consists of three subgraphs b,(b i ),(b ij ), where i,j run over the given bonds of a subgraph. On each bond b,b i,b ij we posit that the nonlinear Dirac equation is satisfied as given by Eq. (). The verte boundary conditions can be written similarly to those in Eq. (3). Soliton solutions have a similar form to those in Eqs. (4)-(5) where is replaced by +s b, with s = s i = l, s ij = l+l i, with l being position of the center of mass of soliton. The sum rules in this case generalize according to: g = i= gi, g i 3 = j= gij, i =,. Figure 6. Triangle graph. (,], i [,L i ] with i =,3,4, and j [,+ ) with j = 5,6. Similarly to that for star graph, one can obtain soliton solutions of the nonlinear Dirac equation on the tree graph for the moving frame. Another graph for which a soliton solution of the nonlinear Dirac equation can be obtained, is presented in Fig. 6. It has the form of a triangle whose vertices are connected to outgoing semi-infinite leads. The verte boundary conditions for Eqs. () and () which follow from the conservation of current and energy can be written as α φ = = α φ = +α 3 φ 3 =, α χ = = α χ = = α 3 χ 3 =, α φ =L = α 4 φ 4 = +α 5 φ 5 =, α χ =L = α 4 χ 4 = = α 5 χ 5 =, α 3 φ 3 =L3 +α 4 φ 4 =L4 = α 6 φ 6 =, α 3 χ 3 =L3 = α 4 χ 4 =L4 = α 6 χ 6 =. In short, the conditions matching the of charge and energy conservations must be applied to each node connecting the multiple vertices of the graph. The soliton solution

12 Dynamics of Dirac solitons in networks (b) (a) bond bond bond.4. 3 bond bond 4. 3 bond bond t.6 5. bond bond bond 4. 3 bond bond t Figure 7. (Color online) Time and coordinate dependence of the charge density obtained from the numerical solution of Eq.() ( = )for the case (a) when constraints (4) (6) are fulfilled (g =, g = g 3 =, g 4 = g 5 =, g 6 = / 3), (b) when these constraints are not fulfilled (g = g 3 = g 5 =, g = g 4 = g 6 = ). obeying these boundary conditions can be written as A b () = g b cosh(β(+s m+ω b)) m+ωcosh(β(+s b )) β m+ωcosh(β(+s b )), () B b () = g b sinh(β(+s m ω b)) m+ωcosh(β(+s b )) β m+ωcosh(β(+s b )), (3) where b is considered to be an inde running from to 6, s = s = s 3 = l, s 4 =

13 Dynamics of Dirac solitons in networks 3 l+l, s 5 = l+l, s 6 = l+l 3, L 3 = L +L 4 (this condition is necessary for the waves traveling along b 3 and those split along b and b 4 to arrive at the verte between b 3 and b 4 concurrently) and coefficients g j fulfill the following constraints: = +, (4) g g g3 = +, (5) g g4 g5 + =. (6) g3 g4 g6 Again, the traveling wave (soliton) solutions can be obtained in case of the moving frame using the Lorentz transformations. This approach for designing soliton solutions of the nonlinear Dirac equation can be applied to other graph topologies with multiple junctions, provided a graph consists of finite parts connected to outgoing semi-infinite bonds. In Fig. 7 the solution of NLDE on each bond of this triangle graph is plotted for the caseswhenthesumrulegivenbyeqs.(4)-(6)isfulfilled(a)andisbroken(b). Absence of the verte reflection in Fig. 7(a) is clearly seen, while in Fig. 7(b) the transmission of Dirac solitons through the graph vertices is clearly impeded by the mismatch at the VBC. Thus, reflection events are clearly discernible in bonds such as, 3 and Conclusions & Future Challenges In this paper we studied dynamics of Dirac solitons in networks by considering the case of metric graphs. We obtained soliton solutions of the nonlinear Dirac equation on some of the simplest metric graphs such as the star, tree and triangle graphs. Constraints enabling such eact solutions to eist are derived in the form of sum rules for bond nonlinearity coefficients. The boundary conditions at the branching points (vertices) are derived from the fundamental conservation laws. It is shown via direct numerical simulations that within the corresponding dynamics the obtained constraints provide for reflectionless transmission of solitary waves at the graph verte. In the case where the relevant conditions are violated, nontrivial reflections ensue. Our computations clearly suggest that the considered verte boundary conditions are necessary (yet not necessarily sufficient) conditions for reflectionless transmission through the junctions of interest. This is also in line with the recent homogenization calculation of [54]. It is thus natural to formulate the following conjecture. Assume

14 Dynamics of Dirac solitons in networks 4 that the initial data correspond to a traveling wave of the NLDE in the form of Eq. () centered at, with v >. Then for t +, the solution can consist asymptotically of a superposition of solitary waves in branches and 3, only if the condition of Eq. (8) holds. This conjecture poses a challenging mathematical question for further rigorous study in the problem at hand, both at the level of the NLDE, but also by direct analogy in the contet of the NLS and other similar models. It should be noted that an important, additional set of challenges emerges as regards the study of stationary solutions involving the verte. In the language of the very recent work of [54], the states considered herein are the so-called half-soliton states. Yet in the NLS framework of the latter work, additional so-called shifted (non-monotone in the different branches) states also eist, yet they may be spectrally unstable, depending on whether they are monotonic or non-monotonic in the outgoing edges of the metric graph. The half-soliton considered here is especially interesting as at the NLS level it turns out to be spectrally stable, yet nonlinearly unstable. This issue of the (more involved at the NLDE level [55]) spectral stability of these stationary states is also an especially interesting direction for future study. Naturally, the above study can be etended for other simple topologies, such as a graph with one or multiple loops (e.g. the dumbbell graph), or other combinations of star, tree and loop graphs. It will be etremely interesting if these ideas induce eperimental interest in systems such as suitably tailored discrete optical waveguides [5, 53] or in atomic settings, in the same way as they have for instance for the propagation of traveling waves through networks of granular crystals [56]. Acknowledgement This work is partially supported by the grant of the Ministry for Innovation Developments of Uzbekistan (Ref. No.BF--). 6. Appendi A Here we will show that the linear verte boundary conditions given by Eqs. (9) lead to the ones given by (6)-(8). Consider the following (linear) relations given at the verte

15 Dynamics of Dirac solitons in networks 5 of the a metric star graph α φ = = α φ = +α 3 φ 3 =, (7) χ = = χ = = χ 3 =, α α α 3 (8) where α, α, α 3 are real constants. Multiplying both sides of Eq. (7) by α χ =, from Eq. (8) we have From Eq. (9) we get φ χ = = φ χ = +φ 3 χ 3 =. (9) Re[φ χ ] = = Re[φ χ ] = +Re[φ 3 χ 3] =, (3) which is nothing but the verte boundary condition given by Eq. (6). Taking the time derivative from Eq. (8) and multiplying both sides of (7) by α t χ =, we obtain φ t χ = = φ t χ = +φ 3 t χ 3 =. (3) Taking the time-derivative from Eq. (7) and multiplying the comple conjugate of the obtained epression by α χ = we get χ t φ = = χ t φ = +χ 3 t φ 3 =. (3) Adding Eqs. (3) and (3) leads to [φ t χ +χ t φ ] = = [φ t χ +χ t φ ] = + [φ 3 t χ 3 +χ 3 t φ 3] =. (33) The last equation yields the verte boundary condition given by Eq. (8). [] Z. Sobirov, D. Matrasulov, K. Sabirov, S. Sawada, and K. Nakamura, Phys. Rev. E 8, 666 (). [] K. Nakamura, Z. A. Sobirov, D. U. Matrasulov, and S. Sawada, Phys. Rev. E 84, 669 (). [3] R. Adami, C. Cacciapuoti, D. Finco, and D. Noja, Rev. Math. Phys. 3, 49 (). [4] R. Adami, C. Cacciapuoti, D. Finco, and D. Noja, Europhys. Lett., 3 (). [5] R. Adami, D. Noja, and C. Ortoleva, J. Math. Phys. 54, 35 (3). [6] K.K. Sabirov, Z.A. Sobirov, D. Babajanov, and D.U. Matrasulov, Phys. Lett. A, 377, 86 (3). [7] D. Noja, Philos. Trans. R. Soc. A 37, 3 (4). [8] G. H. Susanto and S.A. van Gils, Phys. Lett. A 338, 39 (5). [9] J.-G Caputo, D. Dutykh, Phys. Rev. E 9, 9 (4). [] H. Uecker, D. Grieser, Z. Sobirov, D. Babajanov, and D. Matrasulov, Phys. Rev. E 9, 39 (5). [] Z. Sobirov, D. Babajanov, D. Matrasulov, K. Nakamura, and H. Uecker, EPL 5, 5 (6). [] W. Thirring, Ann. Phys. 3, 9 (958).

16 Dynamics of Dirac solitons in networks 6 [3] D. J. Gross, A. Neveu, Phys. Rev. D, 335 (974). [4] W.H. Steeb, W. Oevel, and W. Strampp, J. Math. Phys. 5, 33 (984). [5] I. V. Barashenkov and B. S. Getmanov, Commun. Math. Phys., 43 (987). [6] W.I. Fushchych and R.Z. Zhdanov, J. Math. Phys. 3, 3488 (99). [7] I. V. Barashenkov and B. S. Getmanov, J. Math. Phys. 34, 354 (99). [8] F.M. Toyama, Y. Hosono, B. Ilyas and Y. Nogami, J. Phys. A: Math. Gen. 7, 339 (994). [9] W.K. Ng, Int. J. Mod. Phys. A, 4, 3476 (9). [] W.K. Ng and R.R. Parwani, SIGMA 5, 3 (9). [] L. H. Haddad and L. D. Carr, Physica D 38, 43 (9). [] L. H. Haddad and L. D. Carr, Europhys. Lett (). [3] L.H. Haddad, C.M. Weaver and L.D Carr, New J. Phys. 7, 6333 (5). [4] L.H. Haddad, L.D Carr, New J. Phys. 7, 6334 (5). [5] L.H. Haddad, L.D Carr, New J. Phys. 7, 3 (5). [6] L.H. Haddad, K. M. O Hara, L.D Carr, Phys. Rev. A 9, 4369 (5). [7] F.Cooper, A. Khare, B. Mihaila and A. Saena, Phys. Rev. E 8, 3664 (). [8] F. G. Mertens, N. R. Quintero, F. Cooper, A. Khare, and A. Saena, Phys. Rev. E 86, 466 (). [9] D.E. Pelinovsky, A.Stefanov, J. Math.Phys.,53, 7375 (). [3] D.E. Pelinovsky, Y. Shimabukuro, Lett. Mat. Phys.,4, (4). [3] S. Shao, F. G. Mertens, N. R. Quintero, F. Cooper, A. Khare, and A. Saena, Phys. Rev. E 9, 395 (4). [3] J. Cuevas-Maraver, P.G Kevrekidis and A. Saena, J. Phys. A: Math. Gen. 48, 554 (5). [33] F. G. Mertens, F. Cooper, N.R. Quintero, S. Shao, A. Khare, and A. Saena, J. Phys. A: Math. Gen. 49, 654 (6). [34] J. Cuevas-Maraver, P.G. Kevrekidis and A. Saena, F. Cooper, A. Khare, A. Comech, C.M.Bender, IEEE J. Selected Topics Quant. Electr.,, 59 (6). [35] J. Cuevas-Maraver, P.G Kevrekidis and A. Saena, A.Comech, R.Lan, Phys. Rev. Lett. 6, 4 (6). [36] A. Contreras, D.E. Pelinovsky, Y. Shimabukuro, Commun. Part. Diff. Eq.,4, 7 (6). [37] A. Komech, A. Komech, SIAM J. Math. Anal., 4, 944 (). [38] G. Berkolaiko, A. Comech, Math. Mod. of Natural Phenomena, 7, 3 (). [39] A. Comech, M. Guan and S. Gustafson, Ann. Henri Poincare(Analyse non lineaire), 3 639(4). [4] G. Berkolaiko, A. Comech, A. Sukhtayev, Nonlinearity, (5). [4] J. Cuevas-Maraver, P.G. Kevrekidis, A. Saena, A. Comech, R. Lan, Phys. Rev. Lett. 6 4 (6). [4] Y Zhang, Nonlinear Analysis, 8 8 (5). [43] Ch-W. Lee, P. Kurzynski, H. Nha, Phys. Rev. A (5). [44] T.X. Tran, X.N. Nguyen, F. Biancalana, Phys. Rev. A 9, 384 (8). [45] A.N. Andriotis, M. Menon, Appl.Phys.Lett., 9 45 (8). [46] Y.P. Chen, Y.E. Xie, L. Z. Sun and J. Zhong, Appl. Phys. Lett., (8).

17 Dynamics of Dirac solitons in networks 7 [47] J. G. Xu, L. Wang and M. Q. Weng, J. Appl. Phys., (3). [48] A.G. Kvashnin, D.G. Kvashnin, O.P. Kvashnina and L.A. Chernozatonskii, Nanotechnology, (5). [49] T.X. Tran, S. Longhi, F. Biancalana, Ann. Phys (4). [5] D.N. Christodoulides and E.D. Eugenieva Phys. Rev. Lett. 87, 339 () [5] J. Bolte and J. Harrison, J. Phys. A: Math. Gen (3). [5] M. Conforti, C. De Angelis, T. R. Akylas, and A.B. Aceves Phys. Rev. A (). [53] A. Marini, S. Longhi, and F. Biancalana Phys. Rev. Lett (4). [54] A. Kairzhan, D.E. Pelinovsky, J. Phys. A: Math. Theor. 5, 953 (8). [55] J. Cuevas-Maraver, N. Boussaid, A. Comech, R.Lan, P.G. Kevrekidis, A. Saena, arxiv: [56] A. Leonard, L. Ponson, C. Daraio, J. Mech. Phys. Solids 73, 3 (4).