Multi-kink scattering in the double sine-Gordon model

Vakhid A. Gani; Aliakbar Moradi Marjaneh; Danial Saadatmand

arXiv:1901.07966·hep-th·September 27, 2019

Multi-kink scattering in the double sine-Gordon model

Vakhid A. Gani, Aliakbar Moradi Marjaneh, Danial Saadatmand

PDF

TL;DR

This paper investigates the complex interactions of multiple kinks in the double sine-Gordon model, revealing how initial conditions influence collision outcomes and energy distribution.

Contribution

It introduces a specific initial condition setup for multi-kink collisions and analyzes energy density dependencies and final states in these interactions.

Findings

01

Maximal energy densities depend on the model parameter.

02

Final collision states vary with initial conditions.

03

Multiple kink collisions can be synchronized at a single point.

Abstract

We study collisions of two, three, and four kinks of the double sine-Gordon model. The initial conditions are taken in a special form in order to provide collision of all kinks in one point. We obtain dependences of the maximal energy densities on the model parameter. We also analyze the final states observed in these collisions.

Figures40

Click any figure to enlarge with its caption.

Equations44

L = \frac{1}{2} (\frac{\partial ϕ}{\partial t})^{2} - \frac{1}{2} (\frac{\partial ϕ}{\partial x})^{2} - V (ϕ),

L = \frac{1}{2} (\frac{\partial ϕ}{\partial t})^{2} - \frac{1}{2} (\frac{\partial ϕ}{\partial x})^{2} - V (ϕ),

\frac{\partial ^{2} ϕ}{\partial t ^{2}} - \frac{\partial ^{2} ϕ}{\partial x ^{2}} + \frac{d V}{d ϕ} = 0.

\frac{\partial ^{2} ϕ}{\partial t ^{2}} - \frac{\partial ^{2} ϕ}{\partial x ^{2}} + \frac{d V}{d ϕ} = 0.

V (ϕ) = tanh^{2} R (1 - cos ϕ) + \frac{4}{cosh ^{2} R} (1 + cos \frac{ϕ}{2}) .

V (ϕ) = tanh^{2} R (1 - cos ϕ) + \frac{4}{cosh ^{2} R} (1 + cos \frac{ϕ}{2}) .

V (ϕ) = ⎩ ⎨ ⎧ 4 (1 + cos \frac{ϕ}{2}) \mbox f or R = 0, 4 - (1 - cos \frac{ϕ}{2})^{2} \mbox f or R = \mbox a r s inh 1, 1 - cos ϕ \mbox f or R \to + \infty.

V (ϕ) = ⎩ ⎨ ⎧ 4 (1 + cos \frac{ϕ}{2}) \mbox f or R = 0, 4 - (1 - cos \frac{ϕ}{2})^{2} \mbox f or R = \mbox a r s inh 1, 1 - cos ϕ \mbox f or R \to + \infty.

ϕ_{k (\overset{ˉ}{k})} (x) = 4 π n \pm 4 arctan \frac{sinh x}{cosh R} .

ϕ_{k (\overset{ˉ}{k})} (x) = 4 π n \pm 4 arctan \frac{sinh x}{cosh R} .

ϕ_{k (\overset{ˉ}{k})} (x) = 4 π n \pm [ϕ_{SGK} (x + R) - ϕ_{SGK} (R - x)],

ϕ_{k (\overset{ˉ}{k})} (x) = 4 π n \pm [ϕ_{SGK} (x + R) - ϕ_{SGK} (R - x)],

ϕ_{SGK} (x) = 4 arctan exp (x)

ϕ_{SGK} (x) = 4 arctan exp (x)

ϕ_{k (\overset{ˉ}{k})} (x, t) = ϕ_{k (\overset{ˉ}{k})} (γ (x - v t)),

ϕ_{k (\overset{ˉ}{k})} (x, t) = ϕ_{k (\overset{ˉ}{k})} (γ (x - v t)),

E [ϕ] = - \infty \int \infty [\frac{1}{2} (\frac{\partial ϕ}{\partial t})^{2} + \frac{1}{2} (\frac{\partial ϕ}{\partial x})^{2} + V (ϕ)] d x .

E [ϕ] = - \infty \int \infty [\frac{1}{2} (\frac{\partial ϕ}{\partial t})^{2} + \frac{1}{2} (\frac{\partial ϕ}{\partial x})^{2} + V (ϕ)] d x .

E = K + U + P .

E = K + U + P .

ε (x, t) = k (x, t) + u (x, t) + p (x, t),

ε (x, t) = k (x, t) + u (x, t) + p (x, t),

k (x, t) = \frac{1}{2} (\frac{\partial ϕ}{\partial t})^{2}, u (x, t) = \frac{1}{2} (\frac{\partial ϕ}{\partial x})^{2}, p (x, t) = V (ϕ)

k (x, t) = \frac{1}{2} (\frac{\partial ϕ}{\partial t})^{2}, u (x, t) = \frac{1}{2} (\frac{\partial ϕ}{\partial x})^{2}, p (x, t) = V (ϕ)

p (x, t) = \frac{8 cosh ^{2} R \cdot cosh ^{2} [ γ ( x - v t )]}{( cosh ^{2} R + sinh ^{2} [ γ ( x - v t )] ) ^{2}} = (1 - v^{2}) u (x, t) = \frac{1 - v ^{2}}{v ^{2}} k (x, t) = \frac{1 - v ^{2}}{2} ε (x, t) .

p (x, t) = \frac{8 cosh ^{2} R \cdot cosh ^{2} [ γ ( x - v t )]}{( cosh ^{2} R + sinh ^{2} [ γ ( x - v t )] ) ^{2}} = (1 - v^{2}) u (x, t) = \frac{1 - v ^{2}}{v ^{2}} k (x, t) = \frac{1 - v ^{2}}{2} ε (x, t) .

E_{K} = \int_{- \infty}^{+ \infty} ε (x, t) d x = \frac{M _{K}}{1 - v ^{2}},

E_{K} = \int_{- \infty}^{+ \infty} ε (x, t) d x = \frac{M _{K}}{1 - v ^{2}},

M_{K} = 16 (1 + \frac{2 R}{sinh 2 R})

M_{K} = 16 (1 + \frac{2 R}{sinh 2 R})

p_{max}^{(1)} = (1 - v^{2}) u_{max}^{(1)} = \frac{1 - v ^{2}}{v ^{2}} k_{max}^{(1)} = \frac{1 - v ^{2}}{2} ε_{max}^{(1)} = \frac{8}{cosh ^{2} R}

p_{max}^{(1)} = (1 - v^{2}) u_{max}^{(1)} = \frac{1 - v ^{2}}{v ^{2}} k_{max}^{(1)} = \frac{1 - v ^{2}}{2} ε_{max}^{(1)} = \frac{8}{cosh ^{2} R}

p_{max}^{(1)} = (1 - v^{2}) u_{max}^{(1)} = \frac{1 - v ^{2}}{v ^{2}} k_{max}^{(1)} = \frac{1 - v ^{2}}{2} ε_{max}^{(1)} = 2 coth^{2} R

p_{max}^{(1)} = (1 - v^{2}) u_{max}^{(1)} = \frac{1 - v ^{2}}{v ^{2}} k_{max}^{(1)} = \frac{1 - v ^{2}}{2} ε_{max}^{(1)} = 2 coth^{2} R

\frac{d ^{2} ϕ _{n}}{d t ^{2}}

\frac{d ^{2} ϕ _{n}}{d t ^{2}}

ϕ_{k \overset{ˉ}{k}} (x) = ϕ_{k} (\frac{x + x _{0} - v _{in} t}{1 - v _{in}^{2}}) + ϕ_{\overset{ˉ}{k}} (\frac{x - x _{0} + v _{in} t}{1 - v _{in}^{2}}) - 2 π,

ϕ_{k \overset{ˉ}{k}} (x) = ϕ_{k} (\frac{x + x _{0} - v _{in} t}{1 - v _{in}^{2}}) + ϕ_{\overset{ˉ}{k}} (\frac{x - x _{0} + v _{in} t}{1 - v _{in}^{2}}) - 2 π,

ϕ_{k \overset{ˉ}{k} k} (x) = ϕ_{k} (\frac{x + x _{0} - v _{in} t}{1 - v _{in}^{2}}) + ϕ_{\overset{ˉ}{k}} (x) + ϕ_{k} (\frac{x - x _{0} + v _{in} t}{1 - v _{in}^{2}}),

ϕ_{k \overset{ˉ}{k} k} (x) = ϕ_{k} (\frac{x + x _{0} - v _{in} t}{1 - v _{in}^{2}}) + ϕ_{\overset{ˉ}{k}} (x) + ϕ_{k} (\frac{x - x _{0} + v _{in} t}{1 - v _{in}^{2}}),

ϕ_{k \overset{ˉ}{k} k \overset{ˉ}{k}} (x) = ϕ_{k} (\frac{x + x _{02} - v _{in2} t}{1 - v _{in2}^{2}}) + ϕ_{\overset{ˉ}{k}} (\frac{x + x _{01} - v _{in1} t}{1 - v _{in1}^{2}}) +

ϕ_{k \overset{ˉ}{k} k \overset{ˉ}{k}} (x) = ϕ_{k} (\frac{x + x _{02} - v _{in2} t}{1 - v _{in2}^{2}}) + ϕ_{\overset{ˉ}{k}} (\frac{x + x _{01} - v _{in1} t}{1 - v _{in1}^{2}}) +

+ ϕ_{k} (\frac{x - x _{01} + v _{in1} t}{1 - v _{in1}^{2}}) + ϕ_{\overset{ˉ}{k}} (\frac{x - x _{02} + v _{in2} t}{1 - v _{in2}^{2}}) - 2 π,

+ ϕ_{k} (\frac{x - x _{01} + v _{in1} t}{1 - v _{in1}^{2}}) + ϕ_{\overset{ˉ}{k}} (\frac{x - x _{02} + v _{in2} t}{1 - v _{in2}^{2}}) - 2 π,

Peer Reviews

No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.

Videos

No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.

Full text

Multi-kink scattering in the double sine-Gordon model

Vakhid A. Gani

[email protected]

Department of Mathematics, National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), 115409 Moscow, Russia

Theory Department, National Research Center Kurchatov Institute, Institute for Theoretical and Experimental Physics, 117218 Moscow, Russia

Aliakbar Moradi Marjaneh

Corresponding author, [email protected]

Young Researchers and Elite Club, Quchan Branch, Islamic Azad university, Quchan, Iran

Danial Saadatmand

[email protected]

Department of Physics, University of Sistan and Baluchestan, Zahedan, Iran

Abstract

We study collisions of two, three, and four kinks of the double sine-Gordon model. The initial conditions are taken in a special form in order to provide collision of all kinks in one point. We obtain dependences of the maximal energy densities on the model parameter. We also analyze the final states observed in these collisions.

pacs:

11.10.Lm, 11.27.+d, 05.45.Yv, 03.50.-z

I Introduction

Interactions of kinks — solutions of the type of solitary waves of non-integrable field-theoretical models in (1+1)-dimensional space-time — are of growing interest Rajaraman.book.1982 ; Vilenkin.book.2000 ; Manton.book.2004 ; Vachaspati.book.2006 . Great amount of important results has been obtained recently for the kink-(anti)kink scattering Gani.PRD.2014 ; Gani.JHEP.2015 ; Gani.JPCS.2016.phi8 ; Bazeia.EPJC.2018.sinh ; Bazeia.JPCS.2017.sinh ; Gani.EPJC.2018.dsg ; Belendryasova.JPCS.2019.dsg ; Gomes.JHEP.2018 ; Bazeia.2018.hybrid ; Dorey.PLB.2018.quasinormal ; Campos.2019.quasinormal ; Zhong.2019.inner . In many cases rather non-trivial picture of interaction has been observed. In particular, it was found that the collision of kink and antikink crucially depends on the initial velocity $v_{\mathrm{in}}$ (in the numerical experiments kink and antikink are initially placed at large distance from each other and are moving towards each other with the initial velocities $v_{\mathrm{in}}$ in the laboratory frame of reference). There is a critical value of the initial velocity $v_{\mathrm{cr}}$ , which separates two different regimes of the collision. At $v_{\mathrm{in}}>v_{\mathrm{cr}}$ inelastic reflection of kinks (or passing through each other, if the model allows such process) is observed. At $v_{\mathrm{in}}<v_{\mathrm{cr}}$ after the collision the kink and antikink form a long-living bound state — a bion — which decays slowly radiating energy in the form of small-amplitude waves. Besides, in some models at $v_{\mathrm{in}}<v_{\mathrm{cr}}$ an interesting phenomenon was found — so-called escape windows. Escape windows are intervals of the initial velocity, at which kinks do not form a bion, but escape from each other after two or more collisions. The cause of this phenomenon is the resonance energy exchange between translational and vibrational modes of the kink (antikink), see, e.g., Refs. Campbell.PhysD.1986 ; Kudryavtsev.UFN.1997.eng ; Kudryavtsev.UFN.1997.rus . In some models the kinetic energy of the kinks can be stored in vibrational modes of the system “kink+antikink” Gani.PRE.1999 ; Dorey.PRL.2011 ; Belendryasova.CNSNS.2019.phi8 ; Belendryasova.JPCS.2017.phi8 .

Another very interesting subject that should be mentioned here is interactions of kinks having power-law tails, see, e.g., Refs. Belendryasova.CNSNS.2019.phi8 ; Belendryasova.JPCS.2017.phi8 ; Gomes.PRD.2012 ; Radomskiy.JPCS.2017.tails ; Christov.PRD.2019.long-range ; Christov.2018.long-tails ; Manton.2018.long-range ; Manton.JPA.2018.long-range ; Bazeia.JPC.2018.polynomial-tails . Such kinks arise in various models with both polynomial and non-polynomial potentials. The power-law asymptotics (“fat tails”) of kinks result in their long-range interaction, i.e. the kinks feel each other at very large distances. As a result, the traditionally used initial conditions in the form of superposition of kink and antikink seem to be not valid for simulation of the kink-antikink scattering Christov.PRD.2019.long-range .

The interactions of kinks are studied both numerically and analytically. On the one hand, the scattering processes can be simulated directly by using the numerical methods of solving equations of motion – partial differential equations of second order. On the other hand, there are analytical approximate methods which allow to estimate forces between kink and (anti)kink. These are (i) the collective coordinate approximation Gani.PRD.2014 ; Weigel.PRD.2016.cc ; Demirkaya.JHEP.2107.cc and (ii) Manton’s method (Manton.book.2004, , Ch. 5), Perring.NuclPhys.1962 ; Rajaraman.PRD.1977 ; Manton.NPB.1979 ; Kevrekidis.PRE.2004 .

It is worth to mention that kink-like solutions can also be obtained in more complicated models with two or more scalar fields GaKu.SuSy.2001.eng ; GaKu.SuSy.2001.rus ; Lensky.JETP.2001.eng ; Lensky.JETP.2001.rus ; Alonso-Izquierdo.PRD.2002 ; Kurochkin.CMMP.2004.eng ; Kurochkin.CMMP.2004.rus ; Alonso-Izquierdo.AHEP.2013.kink-like ; Alonso-Izquierdo.JHEP.2014.moduli ; Katsura.PRD.2014 ; Alonso-Izquierdo.Physica_D.2018.dynamics ; Alonso-Izquierdo.PRD.2018.collisions ; Alonso-Izquierdo.2018.MSTB ; Alonso-Izquierdo.2019.scattering ; Morris.AP.2018.mixing ; Alonso-Izquierdo.2019.non-topological ; Mohammadi.CNSNS.2019.two-field-sg . In particular, in GaKu.SuSy.2001.eng ; GaKu.SuSy.2001.rus ; Alonso-Izquierdo.PRD.2002 ; Alonso-Izquierdo.Physica_D.2018.dynamics ; Alonso-Izquierdo.PRD.2018.collisions ; Alonso-Izquierdo.2018.MSTB ; Alonso-Izquierdo.2019.scattering ; Morris.AP.2018.mixing ; Alonso-Izquierdo.2019.non-topological the models with two real scalar fields with polynomial potentials have been studied. In GaLiRa ; GaLiRaconf ; blyankinshtein ; GaKsKu01.eng ; GaKsKu01.rus ; GaKsKu02.eng ; GaKsKu02.rus configurations of the type of domain wall with additional fields localized on it were investigated. Many important results have also been obtained for non-topological field configurations — lumps, Q-balls, etc. Brihaye.PRD.2015 ; Loginov.YadFiz.2011.eng ; Loginov.YadFiz.2011.rus ; Bazeia.EPJC.2017 ; Bazeia.fermion.2017 ; Schweitzer.PRD.2012.1 ; Schweitzer.PRD.2012.2 ; Schweitzer.NPA.2016 ; Nugaev.PRD.2013 ; Bazeia.EPJC.2016.Q-balls ; Bazeia.PLB.2016.Q-balls ; Bazeia.PLB.2017.Q-balls ; Levkov.JHEP.2017.Q-balls ; Loiko.PRD.2018.Q-balls ; Kovtun.PRD.2018.Q-balls ; Loginov.PLB.2018.Q-balls ; Loginov.arXiv.2019.Q-balls .

Kinks, domain walls and other topological defects arise in a great amount of models, hence they are very important for various physical applications from condensed matter to high energy physics and cosmology poltis1 ; valle2016 ; Matsui.JCAP.2016 ; Wachter.PRL.2017 ; Ye.EPJC.2017 ; Ringeval.JCAP.2017 ; GaKiRu ; GaKiRu.conf , see also Vilenkin.book.2000 . Many physical phenomena could be effectively described by one-dimensional topological structures. For example, a three-dimensional domain wall in the direction perpendicular to it can be viewed as a one-dimensional topological field configuration — kink connecting two different vacua of the model. At the same time, the case of $(1+1)$ -dimensional space-time structure can be investigated more easily compared to $2+1$ or $3+1$ dimensions.

The processes of collision of several kinks in one point (with extracting maximal energy densities) have been actively investigated over the last few years saadatmand.prd.2015 ; Moradi.JHEP.2017.multikink ; Moradi.EPJB.2017.multikink ; Moradi.CNSNS.2017.multikink ; saad.WaveMotion.2018 . Such processes are of interest for numerous applications. Being topological defects, the kinks bear significant amount of energy. As a result, simultaneous collision of several kinks (antikinks) in one point can lead to formation of spatial domains with high energy density. The authors of Ref. saad.WaveMotion.2018 have shown that with knowing maximal energy densities in collision of $N$ kinks one can predict how many kink-antikink pairs can be produced in particle collisions.

This study deals with the collisions of several kinks of the double sine-Gordon (DSG) model in one point. This model has been actively studied due to its numerous applications, see, e.g., Gani.EPJC.2018.dsg ; Gani.PRE.1999 ; Campbell.PhysD.1986.dsg ; Malomed.PLA.1989.dsg ; Nazifkar.BJP.2010.dsg and references therein. Note that the subject of multi-kink collisions is very vast. The number of parameters that could be changed is very large. These are, in particular, the model parameter $R$ , the initial velocities and the initial positions of kinks. In this paper we study collisions of two three and four kinks in one point (in a small spatial region). We focus on the dependences of maximal energy densities on the model parameter $R$ . In collisions of four kinks (two kinks and two antikinks, to be more accurate) we investigate how the final configurations depend on $R$ . In our numerical simulations we use fixed initial positions and initial velocities of kinks.

Our paper is organized as follows. In Sect. II we recall the main features of the DSG model, write down its kinks, and discuss their main properties. Section III presents the results of the numerical simulations of multi-kink collisions. We conclude with summarizing in Sect. IV, where we also formulate some possible directions for further research.

II Double sine-Gordon model

Within the $(1+1)$ -dimensional double sine-Gordon model the dynamics of the real scalar field $\phi(x,t)$ is described by the Lagrangian density

[TABLE]

where the potential $V(\phi)$ is non-negative function of the field. From the Lagrangian (1) one can obtain the equation of motion for the field $\phi(x,t)$ :

[TABLE]

Below we simulate multi-kink scattering by solving this equation numerically.

In literature there are various modifications and parameterizations of the potential, which all can be called “double sine-Gordon potential”. We will use the so-called $R$ -parameterization, see, e.g., Refs. Gani.EPJC.2018.dsg ; Gani.PRE.1999 ; Campbell.PhysD.1986.dsg . In this parameterization the DSG potential is:

[TABLE]

In Fig. 1 we show (a) this potential for some values of the parameter and (b) the $R$ -dependence of the coefficients in front of $(1-\cos\phi)$ and $\displaystyle\left(1+\cos\frac{\phi}{2}\right)$ , which are responsible for the “mixing” of $\phi$ – and $\displaystyle\frac{\phi}{2}$ –sine-Gordon potentials.

Depending on $R$ , the potential (3) looks different. In particular,

[TABLE]

The double sine-Gordon model has static topological solution — kink (antikink):

[TABLE]

Notice that for better understanding of many processes in kinks collisions the following fact could be important: the DSG kink (antikink) (5) can be interpreted as a superposition of two sine-Gordon solitons:

[TABLE]

where

[TABLE]

is the sine-Gordon kink (soliton). Thus the DSG kink can be viewed as a superposition of two sine-Gordon kinks, which are separated by the distance $2R$ (if the DSG kink is centered at $x=0$ then the two sine-Gordon solitons are centered at $x=\pm R$ ), see Fig. 2.

Due to the Lorentz invariance, the kink (antikink) moving along the $x$ -axis with the velocity $v$ can be obtained by the Lorentz boost:

[TABLE]

where $\gamma=1/\sqrt{1-v^{2}}$ is the Lorentz factor.

The energy functional which corresponds to the Lagrangian (1) reads

[TABLE]

The total energy (9) can be split into three parts: (i) the kinetic energy $K$ , (ii) the gradient energy $U$ , and (iii) the potential energy $P$ :

[TABLE]

According to this, the integrand in (9), i.e. the total energy density $\varepsilon(x,t)$ , can be written as

[TABLE]

where

[TABLE]

are (i) the kinetic energy density, (ii) the gradient energy density, and (iii) the potential energy density, respectively.

In the case of one moving kink $\phi_{\mathrm{k}}(x,t)$ we have:

[TABLE]

Integrating the total energy density with respect to $x$ , we obtain the total energy of the moving kink:

[TABLE]

where

[TABLE]

is the mass of kink, i.e. the energy of the static DSG kink.

In the next section we present the results of the numerical simulation of collisions of two, three and four DSG kinks in the same point. We find the maximal (over the spatial coordinate $x$ and time $t$ ) values of the energy densities mentioned above: kinetic, gradient, potential and total. Note that for one moving kink the maximal values can be obtained analytically. Depending on the parameter $R$ , they are the following:

[TABLE]

at $R\leq\mbox{arsinh}\,1$ , and

[TABLE]

at $R\geq\mbox{arsinh}\,1$ .

III Multi-kink scattering

We studied the collisions of two kinks (a kink and an antikink), three kinks (two kinks and an antikink), and four kinks (two kinks and two antikinks) at one point. Our main goal was to obtain $R$ -dependences of the maximal values of the energy densities — kinetic, gradient, potential and total. To do this, we numerically solved the discretized version of the equation of motion (2),

[TABLE]

using Störmer method of integration. Here $h=0.025$ is the space step, $n=0,\pm 1,\pm 2,...$ and $\phi_{n}(t)=\phi(nh,t)$ . The second derivative of the field $\phi$ with respect to $x$ , $\phi_{xx}$ , is discretized with the accuracy $O(h^{4})$ . In our numerical experiments we used the range $-375\leq x\leq 375$ . Integration with respect to time was performed with the step $\tau=0.005$ . We used $0\leq t\leq 500$ for searching maximal energy densities over $x$ and $t$ .

The initial conditions for each case (collisions of two, three and four kinks) will be written below. In the case of collisions of three and four kinks, the initial conditions were fit in such a way that the kinks collided in one point.

From the point of view of physical applications, the energy distribution in collision of two or several kinks is of great importance. The point is that a kink is a topological defect, which describes, for example, a planar domain wall or another formation (“defect”) in a real physical system. The kinks collision thus represents mathematical model of collision of such defects, and the energy densities show space-time distribution of energy.

III.1 Kink-antikink scattering

First of all, note that collisions of the DSG kink and antikink have been studied in detail in recent paper Gani.EPJC.2018.dsg . An interesting phenomenon was found — oscillons in the final state. Besides, the dependence of the critical velocity $v_{\mathrm{cr}}$ on $R$ obtained in Gani.EPJC.2018.dsg was somewhat different from that presented in classical paper Campbell.PhysD.1986.dsg . We think that a series of local maxima of the dependence $v_{\mathrm{cr}}(R)$ observed in Gani.EPJC.2018.dsg is a consequence of more detailed investigation of appropriate range of the initial velocities of the colliding kinks.

Unlike the previous works Gani.EPJC.2018.dsg ; Gani.PRE.1999 ; Campbell.PhysD.1986.dsg , we focus on study of maximal energy densities in collisions of the kink and antikink at the initial velocity $v_{\mathrm{in}}=0.1$ and various $R$ . To do this we used the initial condition in the form of the kink and antikink, which are initially placed at $x=-x_{0}$ and $x=+x_{0}$ and moving towards each other with the velocities $v=+v_{\mathrm{in}}$ and $v=-v_{\mathrm{in}}$ , respectively,

[TABLE]

where $\phi_{\mathrm{k}(\mathrm{\bar{k}})}(x)$ is defined by Eq. (5) with $n=0$ . This initial configuration is shown in Fig. 3(a).

In our simulations we used $x_{0}=25$ and $v_{\mathrm{in}}=0.1$ . Space-time picture of the collision for some selected values of $R$ is shown in Fig. 4.

We performed numerical simulation of the kink-antikink collisions within a wide range of the parameter $R$ . We have obtained the dependences of the maximal energy densities on the parameter $R$ , see Fig. 5(d).

These dependences are rather complicated. At $R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}0.5$ and $R\stackrel{{\scriptstyle\scriptstyle>}}{{\phantom{}{}_{\sim}}}2.5$ we see smooth and monotonously decreasing curves. At the same time, within the range $0.5\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}2.5$ the dependences look quite stochastic. The observed difference in behavior of the maximal energy densities at various $R$ could be a consequence of different character of kink-antikink interaction at given initial conditions. In particular, depending on the parameter $R$ , the initial velocity $v_{\mathrm{in}}=0.1$ can be either less or more than the critical value $v_{\mathrm{cr}}$ . The point is that the critical velocity depends on $R$ in a complicated way Gani.EPJC.2018.dsg .

Analyzing the dependence $v_{\mathrm{cr}}(R)$ which is presented in Fig. 6 of Ref. Gani.EPJC.2018.dsg , we can see that the initial velocity $v_{\mathrm{in}}=0.1$ is less than the critical in the range $0.5\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}2.5$ . Thus complicated and almost stochastic behavior of the maximal energy densities is observed for those $R$ ’s at which the kinks form a bound state, i.e., a bion. On the other hand, at those $R$ ’s at which the kinks escape from each other after a collision, the smooth and monotonous dependences are observed. We can assume that the complicated dependences of the maximal energy densities on $R$ can be a consequence of complex resonance processes in bion. As a result, energy is redistributed that leads to abrupt changes of the maximal energy densities.

In order to confirm our hypothesis, we performed numerical simulations of the kink-antikink collisions at $v_{\mathrm{in}}=0.17$ and $v_{\mathrm{in}}=0.2$ . The maximal energy densities are shown in Figs. 6 and 7. From these figures one can see that regions of the stochastic behavior of the maximal energy densities become narrower. As it was for $v_{\mathrm{in}}=0.1$ , at these new initial velocities the stochastic behavior of the energy densities appears at those $R$ ’s at which $v_{\mathrm{in}}<v_{\mathrm{cr}}(R)$ and the kink and antikink capture each other and form a bion in the final state.

Another interesting fact is the presence of smooth and monotonous segments of the maximal energy densities in the range $1.5\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}1.6$ , see Fig. 5(d). We can assume that it is a consequence of the fact that at these $R$ ’s the bound state of kinks evolves differently. The space-time picture of the kink-antikink collision at some selected values of the parameter $R$ from the range $1.49\leq R\leq 1.63$ is shown in Fig. 8. It it clearly seen that at $1.5\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}1.6$ two escaping oscillons in the final state are formed after two oscillations of the kink-antikink bound state. At other values of $R$ in Fig. 8 the oscillons in the final state are either formed after more complicated evolution of the kink-antikink bound state or not formed at all.

III.2 Collision of two kinks and an antikink

We performed numerical simulation of the collision of three DSG kinks in one point. We used the following initial configuration:

[TABLE]

where $\phi_{\mathrm{k}(\mathrm{\bar{k}})}(x)$ is defined by Eq. (5) with $n=0$ . This initial condition corresponds to two kinks and an antikink. The antikink is at rest at the origin $x=0$ , and two kinks are moving towards the antikink with the initial velocities $v=+v_{\mathrm{in}}$ and $v=-v_{\mathrm{in}}$ , starting from $x=-x_{0}$ and $x=+x_{0}$ , respectively. This initial configuration is shown in Fig. 3(b). In our simulations we used $x_{0}=25$ and $v_{\mathrm{in}}=0.1$ . Space-time picture of the collision for some selected values of $R$ is shown in Fig. 9.

The interaction of three kinks looks rather complicated. In particular, in a pair “kink+antikink” resonance interaction can occur. This leads to appearance of various intermediate and final states, see Fig. 9. (Note that the term “final state” stands here for configuration observed after a long time after the collision. “Long time”, in turn, means a period of time which is much longer than life-times of various bound states (of oscillons, of kink and oscillon, and so on), if such states do not have infinite life-times.)

We have obtained the dependencies of the maximal energy densities on the parameter $R$ , see Fig. 10.

The maximal potential energy density is monotonously decreasing, while the maximal kinetic, gradient and total energy densities behave in a complicated way, thereby reflecting complex interaction of the three colliding kinks.

Note that in the range $1.7\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}2.1$ we see a smooth and monotonous segments of the maximal energy densities. We can assume that in this range of $R$ the field evolution is somewhat different. Actually, from Fig. 11 we see that at these $R$ ’s two escaping oscillons in the final state are formed right after the kinks collision. At the same time, at $R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}1.7$ and $R\stackrel{{\scriptstyle\scriptstyle>}}{{\phantom{}{}_{\sim}}}2.1$ the escape of oscillons occur after some time during which the oscillons remain bound.

III.3 Collision of two kinks and two antikinks

Finally, we studied the collision of four DSG kinks in one point. To do this, we used the initial configuration of the form of

[TABLE]

where $\phi_{\mathrm{k}(\mathrm{\bar{k}})}(x)$ is defined by Eq. (5) with $n=0$ . The initial configuration (21) corresponds to two kinks and two antikinks moving to the collision point, see Fig. 3(c). We used the following initial values of positions and velocities of the kinks and antikinks: $x_{01}=10$ , $x_{02}=25$ , $v_{\mathrm{in1}}=0.05$ , $v_{\mathrm{in2}}=0.1$ . This choice of the parameters guarantees almost simultaneous collision of all four waves in the origin $x=0$ .

The space-time picture of the collision for some selected values of the parameter $R$ is presented in Fig. 12.

We can see a vast variety of final states. At $R=1.5$ , $R=2.0$ and $R=2.5$ we observe four oscillons in the final state. We could propose the following explanation of the formation of oscillons. Each DSG kink/antikink consists of two subkinks, see Eq. (6). In the kinks collision the constituent subkinks interact with each other in a complicated way, thus forming bound states — oscillons. This means that we consider an oscillon as a bound state of subkinks (at the same time, remind that a bion is a bound state of kink and antikink). At $R=0.5$ we observe formation of two bions. At $R=1.0$ the final state has the form of escaping kink and antikink together with two oscillons also escaping from the collision point. At $R=3.0$ in the final state we see a bion localized near the origin $x=0$ and escaping DSG kink and antikink. Maximal values of the energy densities in the case of four kinks collisions are highly dependent on the initial positions of the kinks. That is why we did not study the energy densities here but briefly summarized the observed final states.

IV Discussion and conclusion

We have studied collisions of several kinks of the double sine-Gordon model in one point (in a small spatial region). We have performed numerical simulations of

(a)

collisions of a kink and an antikink,

(b)

collisions of two kinks and an antikink in one point,

(c)

collisions of two kinks and two antikinks in one point.

Collisions of two kinks (i.e. kink-antikink scattering) are thoroughly investigated in literature Gani.EPJC.2018.dsg ; Gani.PRE.1999 ; Campbell.PhysD.1986.dsg . In this paper we focused on energy distributions in the kink-antikink collisions at various values of the parameter $R$ . We have obtained $R$ -dependences of the maximal (over $x$ and $t$ ) kinetic, gradient, potential and total energy densities, Fig. 5(d). The maximal kinetic, gradient and total energy densities behave rather stochastic at $0.5\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}2.5$ , while at $R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}0.5$ and $R\stackrel{{\scriptstyle\scriptstyle>}}{{\phantom{}{}_{\sim}}}2.5$ smooth and monotonous segments are observed. We suppose that such behavior could be a consequence of the fact that the initial velocity of the kinks $v_{\mathrm{in}}=0.1$ is less than the critical velocity $v_{\mathrm{cr}}$ , see Gani.EPJC.2018.dsg . This, in turn, leads to kinks’ capture and formation of a bion. Energy redistribution in a bion is complicated and can be viewed as chaotic. At the same time, at $R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}0.5$ and $R\stackrel{{\scriptstyle\scriptstyle>}}{{\phantom{}{}_{\sim}}}2.5$ we have $v_{\mathrm{in}}>v_{\mathrm{cr}}$ , and colliding kinks bounce off each other after the first impact. In this case the complex energy redistribution does not take place, and $R$ -dependences of the maximal energy densities are smooth and monotonous.

Besides, at $1.5\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}1.6$ we also observed small smooth and monotonous segments of the maximal energy densities in Fig. 5(d). In this range of the parameter, unlike other close values of $R$ , in the final state we have two escaping oscillons formed after two oscillations of the kink-antikink bound state.

In three kinks collisions (in collisions of two kinks and an antikink) we observed various final states, in particular, similar to (i) excited kink ( $R=0.5$ , Fig. 9(a)), (ii) excited kink and escaping oscillons ( $R=1.0,1.5,2.0,2.5$ , Figs. 9(b)–9(e)), and (iii) two escaping kinks and an antikink at rest ( $R=3.0$ , Fig. 9(f)). The maximal kinetic, gradient and total energy densities behave stochastic. At the same time, within the range $1.7\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}R\stackrel{{\scriptstyle\scriptstyle<}}{{\phantom{}{}_{\sim}}}2.1$ we observed smooth and monotonous segments of the curves. We found that it could be a consequence of the formation of two escaping oscillons in the final state right after the kinks collision. Notice that the maximal total energy density in three-kink collisions is considerably larger than in two-kink collisions, see Figs. 5(d) and 10. At the same time, the maximal potential energy density is nearly the same for the two- and three-kink collisions. It means that the collision process does not affect on the maximal potential energy density. It is interesting that similar situation was observed in the multi-kink collisions in the sine-Gordon model Moradi.EPJB.2017.multikink (up to seven kinks collision) and in the collisions of odd number of kinks in the $\varphi^{4}$ (up to five kinks collision) Moradi.CNSNS.2017.multikink , while in the collisions of even number of kinks in the $\varphi^{4}$ (up to four kinks collision) and in the $\varphi^{6}$ multi-kink collisions the maximal potential energy depends on the number of colliding kinks Moradi.JHEP.2017.multikink ; Moradi.CNSNS.2017.multikink .

Finally, we have studied collisions of four kinks (two kinks and two antikinks). The initial positions and velocities of the kinks were fit in such a way in order to provide almost simultaneous collision of all four waves in one point. After collision we observed a vast variety of the final states. Besides, the maximal energy densities strongly depend on the initial positions of the colliding kinks. That is why the analysis of the maximal energy densities requires detailed investigation of their dependences on parameters of the initial condition, which could be a subject of a separate publication. In this paper we have limited ourselves to the classification of the most typical final configurations. In particular, at $R=0.5$ in the final state we observed two oscillating structures which we classified as bound states of the DSG kink and antikink. At $R=1.0$ we observed escaping kink and antikink along with two oscillons — bound states of subkinks, which are components of the DSG kink (antikink). At $R=1.5$ and $R=2.0$ — only oscillons, at $R=3.0$ — escaping antikink and kink, while the other kink-antikink pair forms a bion localized near the collision point.

So we see a vast variety of final states in collisions of more than two kinks. Moreover, the behavior of the maximal energy densities strongly depends on oscillating structures being formed in the final state.

In conclusion, we would like to mention the following issues which are beyond the scope of our paper, nevertheless, in our opinion, could become subject to future work.

•

First, the multi-kink scattering within the field-theoretic models with polynomial potentials of 8th or higher degrees. As we have already mentioned in the Introduction, such models can have kinks with long-range interaction between them due to the kinks’ power-law tails. This long-range interaction could lead to new collective phenomena in multi-kink systems.

•

Second, the multi-kink scattering within models with non-polynomial potentials, e.g., sinh-deformed $\varphi^{4}$ or $\varphi^{6}$ , could be of interest, and will be studied in a separate publication.

Acknowledgments

The work of MEPhI group was supported by the MEPhI Academic Excellence Project (Contract No. 02.a03.21.0005, 27.08.2013). V.A. Gani also acknowledges the support of the Russian Foundation for Basic Research under Grant No. 19-02-00971.

Bibliography93

The reference list from the paper itself. Each links out to its DOI / PubMed record.

1(1) R. Rajaraman, Solitons and Instantons: An Introduction to Solitons and Instantons in Quantum Field Theory (North-Holland, Amsterdam, 1982)
2(2) A. Vilenkin, E.P.S. Shellard, Cosmic Strings and Other Topological Defects (Cambridge University Press, Cambridge, 2000)
3(3) N. Manton, P. Sutcliffe, Topological Solitons (Cambridge University Press, Cambridge, 2004)
4(4) T. Vachaspati, Kinks and Domain Walls: An Introduction to Classical and Quantum Solitons (Cambridge University Press, Cambrifge, 2006)
5(5) V.A. Gani, A.E. Kudryavtsev, M.A. Lizunova, Kink interactions in the (1+1)-dimensional φ 6 superscript 𝜑 6 \varphi^{6} model. Phys. Rev. D 89 , 125009 (2014) . ar Xiv:1402.5903 · doi ↗
6(6) V.A. Gani, V. Lensky, M.A. Lizunova, Kink excitation spectra in the (1+1)-dimensional φ 8 superscript 𝜑 8 \varphi^{8} model. JHEP 08 , 147 (2015) . ar Xiv:1506.02313 · doi ↗
7(7) V.A. Gani, V. Lensky, M.A. Lizunova, E.V. Mrozovskaya, Excitation spectra of solitary waves in scalar field models with polynomial self-interaction. J. Phys.: Conf. Ser. 675 , 012019 (2016) . ar Xiv:1602.02636 · doi ↗
8(8) D. Bazeia, E. Belendryasova, V.A. Gani, Scattering of kinks of the sinh-deformed φ 4 superscript 𝜑 4 \varphi^{4} model. Eur. Phys. J. C 78 , 340 (2018) . ar Xiv:1710.04993 · doi ↗