Sine-Gordon solitons in networks: Scattering and transmission at vertices

# Sine-Gordon solitons in networks: Scattering and transmission at vertices

Zarif Sobirov, Doniyor Babajanov, Davron Matrasulov, Katsuhiro Nakamura, Hannes Uecker sobirovzar@gmail.comhannes.uecker@uni-oldenburg.de Tashkent Financial Institute, 60A, Amir Temur Str., 100000, Tashkent, Uzbekistan
Turin Polytechnic University in Tashkent, 17 Niyazov Str., 100095, Tashkent, Uzbekistan
Faculty of Mathematics, National University of Uzbekistan, Vuzgorodok, Tashkent 100174,Uzbekistan
Faculty of Physics, National University of Uzbekistan, Vuzgorodok, Tashkent 100174,Uzbekistan
Department of Applied Physics, Osaka City University, Osaka 558-8585, Japan
Institut für Mathematik, Universität Oldenburg, D26111 Oldenburg, Germany
###### Abstract

We consider the sine-Gordon equation on metric graphs with simple topologies and derive vertex boundary conditions from the fundamental conservation laws together with successive space-derivatives of sine-Gordon equation. We analytically obtain traveling wave solutions in the form of standard sine-Gordon solitons such as kinks and antikinks for star and tree graphs. We show that for this case the sine-Gordon equation becomes completely integrable just as in case of a simple chain. This simple analysis provides a cornerstone for the numerical solution of the general case, including a quantification of the vertex scattering. Applications of the obtained results to Josephson junction networks and DNA double helix are discussed.

###### pacs:
05.45.Yv,02.30.Ik,42.65.Tg

Introduction. Nonlinear wave dynamics described by the sine-Gordon equation is of importance in a variety of topics in physics, such as as elastic and stress wave propagation in solids, liquids and tectonic plates (see, e.g., Ablowitz and Segur (2006); Rajaraman (1982); Drazin and Johnson (1989); Ablowitz and Clarkson (1991); Kivshar and Agrawal (2003); Dauxois and Peyrard (2006); Cuevas-Maraver et al. (2014)), transport in Josephson junctions Barone and Paternò (1982), and topological quantum fields Rajaraman (1982); Dauxois and Peyrard (2006). Continuous and discrete forms of the sine-Gordon equation have been used so far for the description of wave transport in different media. However, there are structures for which the wave dynamics cannot be described within the traditional continuous or discrete approaches. These are networks and branched structures where the transmission through a branching point (network vertex) should be described by vertex conditions. Early studies of nonlinear evolution equations in branched structures are Christodoulides and Eugenieva (2001); Cuevas and Kevrekidis (2004); Kevrekidis et al. (2003), and in recent few years one can observe rapidly growing interest in nonlinear waves and soliton transport in networks described by nonlinear Schrödinger equation Sobirov et al. (2010); Adami et al. (2011, 2012a, 2012b, 2013); Sabirov et al. (2013). Integrable boundary conditions following from the conservation laws were formulated, and soliton solutions yielding reflectionless transport across the graph vertex were derived in Sobirov et al. (2010), see also Nakamura et al. (2011) for the case of a discrete nonlinear Schrödinger equation. Burioni et al Bur (2006); Burioni et al. (2005) studied the discrete nonlinear Schrödinger equation and computed transport and reflection coefficients as a function of the wavenumber of a Gaussian wave packet and the length of a graph attached to a defect site.

In this paper we address the wave dynamics in networks described by the sine-Gordon equation on metric graphs, which can be used for modeling of soliton transport in DNA double helix, tectonic plates and Josephson junction networks. The latter has attracted much attention in condensed matter physics Giuliano and Sodano (2009, 2013). Another interesting application of sine-Gordon equations, or, more generally, nonlinear Klein-Gordon equations, on metric graphs can be networks of granular chains Kivshar et al. (2003); Kevrekidis et al. (2003). Recently, soliton dynamics in networks was studied by considering the 2D sine-Gordon equation on and junctions Caputo and Dutykh (2014), and the metric graph limit was also studied numerically. See also Uecker et al. (2015) for similar results for the 2D Nonlinear Schrödinger equation on “fat” graphs.

Here we focus on the problem of integrability of sine-Gordon equations on metric graphs and soliton transmission at the graph vertex. In particular, using an approach similar to that of Sobirov et al. (2010), we discuss the conditions under which the sine-Gordon equation is completely integrable and allows exact traveling wave solutions which provide reflectionless transmission of sine-Gordon solitons across vertices. Numerical solutions with scattering at a vertex when these conditions are violated are also presented.

Conservation laws and boundary conditions. For evolution equations on graphs, the connections of the bonds at the vertices are provided by vertex boundary conditions. In case of linear wave equations such conditions follow from self-adjointness of the problem Kostrykin and Schrader (1999); Exner and Kovařík (2015). For nonlinear evolution equations one should use such fundamental laws as energy, flux, momentum and (for sine-Gordon model) topological charge conservations Sobirov et al. (2010); Adami et al. (2011); Caputo and Dutykh (2014). Below we derive such conditions and show the existence of infinitely many conservation laws in our model, which yields the complete integrability of the system.

Most of the 1D sine-Gordon models follow from the Lagrangian density where and are positive constants.

We want to explore a sine-Gordon model on networks modeled by graphs, i.e., system of bonds which are connected at one or more vertices (branching points). The connection rule is called the topology of the graph. When the bonds can be assigned a length, the graph is called a metric graph. The sine-Gordon model on each bond , is given by Lagrangian density

 Lk=[12(u2kt−a2ku2kx)−βk(1−cosuk)],

where . In the following we consider a star graph with three semi-infinite bonds connected at the point called vertex of the graph, see Fig. 1. The coordinates are defined as and , where 0 corresponds to the vertex point. The equation of motion derived from the above Lagrangian density leads to the sine-Gordon equation on each bond given as

 uktt−a2kukxx+βksinuk=0. (1)

To formulate physically reasonable vertex boundary conditions (VBC) one can use the continuity of wave function

 u1(0,t)=u2(0,t)=u3(0,t) (2)

and fundamental conservation laws such as energy, charge and momentum conservations together with the asymptotic conditions at infinities: and as , and and as , for some integer , .

For the primary star graph in Fig 1, the energy is defined as

 E(t)=3∑k=11βk∫Bk[12(u2kt+a2ku2kx)+βk(1−cosuk)]dx, (3)

where . Then

 ˙E=a21β1u1xu1t∣∣ ∣∣x1=0−a22β2u2xu2t∣∣ ∣∣x2=0−a23β3u3xu3t∣∣ ∣∣x3=0,

and by (2) the energy conservation reduces to

 a21β1u1x∣∣ ∣∣x1=0=a22β2u2x∣∣ ∣∣x2=0+a23β3u3x∣∣ ∣∣x3=0. (4)

For the same star graph as in Fig 1, the charge is given by

 2πQ=a1√β10∫−∞u1xdx+3∑k=2ak√βk+∞∫0ukxdx. (5)

From and (2) we obtain the sum rule

 a1√β1=a2√β2+a3√β3. (6)

The initial boundary problem (IBVP) (1), (2) and (4) with appropriate initial conditions and asymptotic conditions at infinities is now well defined. However, to see the infinite number of constants of motion, we must search for other additional conditions for parameters.

Integrability and traveling wave solutions. We consider the momentum defined by

 P=3∑k=1akβk∫Bkukxuktdx, (7)

such that

 ˙P = a1β1[12(u21t+a21u21x)−β1(1−cosu1)]∣∣∣x1=0 (8) − 3∑k=2akβk[12(u2kt+a2ku2kx)−βk(1−cosuk)]∣∣∣xk=0.

For , we impose the condition,

 β1=β2=β3=β(>0), (9)

which simplifies the sum rule (6) to

 a1=a2+a3. (10)

Then (8) becomes

 2β˙P = a31u21x(0,t)−a32u22x(0,t)−a33u23x(0,t) = −a2a3a2+a3(a2u2x(0,t)−a3u3x(0,t))2,

where we used (4), (9) and (10) in obtaining the last expression. Thus, yields , which together with (4), (9) and (10) gives

 a1u1x(0,t)=a2u2x(0,t)=a3u3x(0,t). (11)

Conditions on higher-order space-derivatives may be available from higher-order conservations, where the analysis becomes more laborious. However, they can also be obtained directly from (1), (2) via Thus,

 a21u1xx(0,t)=a22u2xx(0,t)=a23u3xx(0,t), (12)

and similarly, taking successive derivatives of (1), we obtain the conditions on higher-order space-derivatives. It should be emphasized that the momentum conservation requires (9)-(11), from which (12) follows. Now we shall prove that equations (2), (11) and (12) give a scaling function which guarantees the infinite number of constants of motion.

Let us introduce two functions defined on the bonds from 1 to as for and for . Thanks to the vertex boundary conditions (VBC) in (11) and (12), together with the continuity condition in (2), both of with and satisfy , where is a solution of the dimensionless sine-Gordon equation

 vtt−vxx+sinv=0 (13)

defined on the real line. This fact is identical to the expression of in terms of the function as

 uk(x,t)=v(√βakx,√βt),x∈Bk(k=1,2,3). (14)

The scaling function in (14) together with the sum rule (10) guarantees the infinite number of constants of motion and thereby the complete integrability of the sine-Gordon equation on the network.

In fact, from (14) and the sum rule (10) we find that all the conservation laws Sanuki and Konno (1974); Geicke (1981)

 +∞∫−∞g(v,∂xv,∂tv,∂2xv,...,∂nx∂ltv)dx=const (15)

of the sine-Gordon equation (13) on the real line also hold on the star graph, because

 3∑k=1 ∫Bkg(uk,akβ−12∂xuk,β−12∂tuk,...,ankβ−n+l2∂nx∂ltuk)dx = a1+∞∫−∞g(v,∂xv,∂tv,…,∂nx∂ltv)dx +(a2+a3−a1)+∞∫0g(v,∂xv,∂tv,…,∂nx∂ltv)dx =const. (16)

The conservation of energy , charge and momentum are just special cases.

From now on, we shall prescribe without loss of generality. Eq.(13) has a number of explicit soliton solutions, for instance: kink “+” and anti-kink “-” solutions which can be written as Drazin and Johnson (1989); Ablowitz and Clarkson (1991)

 (17)

where is the velocity of the kink. Other soliton solutions include breathers, kink-kink collisions and kink-antikink collisions Ablowitz and Clarkson (1991), to name just a few, see also Saadatmand et al. (2015) for further multi-soliton type solutions. If the sum rule in (10) holds, then all these solutions, or, more generally, all solutions of (13), transfer via (14) to solutions on the metric graph. For instance, the kinks then provide reflectionless transmission of energy through the graph vertex, where the speed and energy of a kink traveling in positive direction () split according to the ratios and , respectively. On the other hand, launching two suitably fine tuned kinks on bonds 2 and 3 in negative direction, their joint energy is transmitted to bond 1.

Before entering into the numerical analysis of kink dynamics, we comment on other VBCs originating in the local scattering properties at each vertex. The VBC (2) of continuity and (4) of local flux conservation with are also called VBC. They naturally appear (Caputo and Dutykh (2014), see also Uecker et al. (2015) for a similar construction for the case of the NLS, and (Exner and Kovařík, 2015, Chapter 8) for an overview of related linear results) by considering the 2D sine–Gordon equation on a “fat” graph, i.e., a 2D branched domain with Neumann boundary conditions, where and are the relative widths of the (fat) bonds.

Similarly, the so–called VBCs (Exner and Kovařík, 2015, Chapter 8) consist of (11) and

 a1u1(0,t)−a2u2(0,t)−a3u3(0,t)=0, (18)

which conserve charge and energy for all values of the . A simple calculation shows both and VBCs can be derived from (10) and (14), but the inverse derivation is not possible. We note that (11) and (18) conserve and , but if conservation of is enforced, then Eq.(18) reduces to (2). Most importantly, (10) and (14) give the existence of the infinite number of constants of motion (as shown in (15),(16)), which is equivalent to the complete integrability of the sine–Gordon equation on the graph.

Vertex transmission. An important issue for wave dynamics in networks is the scattering at vertices. The sum rule in (10) allows the tuning of the vertex scattering to achieve reflectionless transmission. We now give numerical solutions of (1) with , using 2nd order in space and time finite differences, where we first focus on (2) and (4) as VBC, i.e., the case. Figure 2 shows the reflectionless propagation of a kink in the special case that the sum rule (10) holds. Figure 2: (Color online). Numerical solution of (1), (2), (4), with βk=1, k=1,2,3, and a1=1, a2=0.7, a3=0.3 fulfilling (10). The initial conditions belongs to the kink (17) on bond 1 with x0=−5 and ν=0.9, while u2,3=∂tu2,3≡0. Panels (a),(b),(d) arranged corresponding to Fig. 1, while (c) shows the t–dependence of the energies.

In Fig. 3 we numerically treat the transmission of solitons through the graph vertex when the sum rule (10) is violated. In (a)-(c) we consider the “natural” case , , and the same kink initial condition as in Fig.2. The total energy is still conserved (by (4)), but in contrast to the reflectionless case from Fig.2, there now is significant reflection at the vertex. To demonstrate and quantify the dependence of the vertex transmission on the in some more detail, in Fig. 3(d) we essentially return to the situation of Fig. 2. That is, we set , , but let vary and plot the reflection coefficient , defined as the ratio of the energies in the incoming bond at initial time and at . At , corresponding to Fig. 2, we have , i.e. zero reflection.

Additionally, the red line in Fig. 3(d) shows the analogous simulation for the case of vertex conditions (11) and (18). Again we have zero reflection at , while violating the sum rule gives qualitatively similar but slightly stronger reflections than the case. Two points should be noted: the simulations in Figs. 2 and 3 have also confirmed the conservation up to numerical discretization effects) of and so long as the sum rule (10) holds (see Fig.4); the simulations do not use the soliton properties of the kinks, but only the fact that they are traveling wave solutions for which we have formulas for the initial conditions. Thus, these numerical results can be transfered to general nonlinear Klein-Gordon equations that admit travelling wave solutions. Figure 3: (Color online). (a)-(c) Numerical solution of (1), (2), (4), with βk=1 for k=1,2,3. (a)-(c) Reflection of incoming kink at the vertex in case that a1=a2=a3=1, violating (10), initial conditions as in Fig. 2. u3 is identical to u2, and hence u3 and E3 in (c) are omitted. (d) (blue line) Dependence of the vertex reflection coefficient R=E1|t=15/E1|t=0 on a3, where a1=1,a2=0.7, hence a3=0.3 corresponding to Fig. 2. The red line is similarly obtained from solving (1) with δ′ VBC (11) and (18).

Other graph topologies. Our results can be extended to other simple topologies such as general star graphs, tree graphs, loop graphs and their combinations. Figure 4: (Color online). Time evolution of momenta. Left panel: a1=1, a2=0.7, a3=0.3; right panel: a1=a2=a3=1.

Exact traveling wave solutions of sine-Gordon models on such graphs with one incoming semi-infinite bond can be obtained similarly to the above case of a star graph with three bonds, leading to generalizations of the sum rule. We illustrate this for the tree graph from Fig. 5, consisting of three “layers” , where run over the given bonds.

On each bond we have a sine-Gordon equation given by (1). Setting for all , the and have to be determined from the sum rule like (10) at each vertex. For instance, at the three nodes in Fig. 5 we need

 end of b1:a0=a11+a12, end of b11:a11=a111+a112+a113, end of b12:a12=a121+a122, (19)

and this continues through the layers. By (13) and (14) this is based on scalings such as

 u1(x,t)=v(x/a1,t) and u1i(x,t)=v(x/a1i,t), (20)

where at subsequent bonds we also need to take into account the finite propagation length in the previous bonds, for instance

 u111(x,t)=v((x+x0)/a111,t), x0/a111=L1/a11, (21)

i.e. . Necessarily, the speeds and energies of, e.g., an incoming kink, also split according to rules like (19), such that on each final bond we only have slow and small energy kinks. A similar construction has been done and formalized for the propagation of Nonlinear Schrödinger solitons on tree graphs in Sobirov et al. (2010).

Another graph for which soliton solutions of sine-Gordon models can be obtained is a graph with a loop (see Fig.6), which consists of two semi-infinite bonds connected by bonds having finite lengths . Requiring the conditions

 a0=n∑k=1ak=an+1 (22)

for the coefficients, and () with a constant , we can write soliton solutions in a similar way as in (20) and (21).

Finally, it can be shown that the above approach can be applied to obtain exact traveling wave solutions of sine-Gordon models on other (than above) graphs consisting of at least two semi-infinite bonds and any subgraph between them. In this case one needs to impose the pertinent vertex conditions like (20) or (22) at the vertices connecting the semi-infinite bonds with the subgraph.

Conclusions. In this work we studied sine-Gordon equations on simple metric graphs, and derived vertex boundary conditions for charge, energy and momentum conservation, and additionally conditions on parameters, for which the problem has explicit analytical soliton solutions. We find the sum rule (10) for bond-dependent coefficients at each vertex of the graph, which makes the sine-Gordon equation on the graph completely integrable. It is shown that the obtained solutions provide the reflectionless transmission of solitons at the graph vertex. This is also illustrated numerically by quantifying the reflections for a case where these conditions are violated, and we discussed how to generalize the results to other graph topologies. The results can be directly applied to several important problems such as Josephson junction network and DNA double helix. In such approach our model corresponds to continuous version of the system considered in Giuliano and Sodano (2009, 2013). Finally, a very important application can be DNA double helix models where the energy transport is described in terms of sine-Gordon equations Yomosa (1983); Yakushevich (2004). Base pairs of the DNA double helix can be considered as a branched system and modeled by a star graph Yakushevich (2004). Then the -bond energy between two base pairs in such system can be characterized by the parameter, .

Acknowledgement. We thank Panayotis Kevrekidis for his valuable comments on this paper. This work is supported by a grant of the Volkswagen Foundation. The work of DM is partially supported by the grant of the Committee for the Coordination Science and Technology Development (Ref.Nr. F3-003).

## References

You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters   