CC BY
3Russian Journal of Nonlinear Dynamics, 2019, vol. 15, no. 1, pp. 21-34. Full-texts are available at http://nd.ics.org.ru DOI: 10.20537/nd190103
NONLINEAR PHYSICS AND MECHANICS
MSC 2010: 35C08, 35Q51, 65M06
Excitation of Large-Amplitude Localized Nonlinear Waves by the Interaction of Kinks of the Sine-Gordon Equation with Attracting Impurity
A. M. Gumerov, E. G. Ekomasov, R. V. Kudryavtsev, M. I. Fakhretdinov
The generation and evolution of localized waves on an impurity in the scattering of a kink of the sine-Gordon equation are studied. It is shown that the problem can be considered as excitation of oscillations of a harmonic oscillator by a short external impulse. The external impulse is modeled by the scattering of a kink on an impurity. The influence of the modes of motion of a kink on the excitation energy of localized waves is numerically and analytically studied. The method of collective coordinate for the analytical study is used. The value of this energy is determined by the ratio of the impurity parameters and the initial kink velocity. It is shown that the dependence of the energy (and amplitude) of the generated localized waves
Received January 27, 2019 Accepted March 21, 2019
This work was supported by Act 211 of the Government of the Russian Federation, contract No. 02.A03.21.0011. For A.M. Gumerov and R. V. Kudryavtsev the work was supported by the RFBR grant, project 18-31-00122.
Azamat M. Gumerov article@solitonlab.com Marat I. Fakhretdinov article@solitonlab.com
Bashkir State University
ul. Zaki Validi 32, Ufa, 450076 Russia
Evgenii G. Ekomasov article@solitonlab.com
Bashkir State University
ul. Zaki Validi 32, Ufa, 450076 Russia
South Ural State University (National Research University)
pr. Lenina 76, Chelyabinsk, 454080 Russia
Roman V. Kudryavtsev article@solitonlab.com
Bashkir State University ul. Zaki Validi 32, Ufa, 450076 Russia Institute of Molecule and Crystal Physics pr. Oktyabrya 151, Ufa, 450075 Russia
on the initial kink velocity has only one maximum. This behavior is observed for the cases of point and extended impurities. Analytical expression for the amplitude of the localized wave in the case of point impurity is obtained. This allows controlling the excitation energy of localized waves using the initial kink velocity and impurity parameters. The study of the evolution of localized impurities under the action of an external force and damping has shown a good agreement with the nondissipative case. It is shown that small values of the external force have no significant effect on the oscillations of localized waves. An analytical expression for the logarithmic decrement of damping is obtained. This study may help to control the parameters of the excited waves in real physical systems.
Keywords: sine-Gordon equation, impurity, kink, wave generation
1. Introduction
Soliton solutions of nonlinear differential equations attract increased attention of researchers because of an increasing use in physical applications [1-18]. For example, solitons of the sine-Gordon equation(the SGE) are frequently used. Only in solid state physics they describe domain boundaries in magnets, dislocations in crystals, fluxons in Josephson junctions. The sine-Gordon equation appears in the modeling of wave processes in rocks, molecular biology, models of field theory and particle physics [2, 3]. Accounting for the influence of perturbations on the SGE solutions leads to a change in the structure of solitons and to the excitation of the internal degrees of freedom of solitons. It is important to study the SGE solutions considering the perturbations arising in real physical applications. The effect of small perturbations on the SGE solutions can be studied using a well-developed perturbation theory for solitons [5, 19, 20]. The effect of perturbations can generally be investigated only with the help of numerical methods [20-23]. For example, the influence of a coordinate- and time-dependent external force has been studied [24-26].
Numerous works have studied spatially dependent perturbations in the sine-Gordon model [5, 6, 18-43]. Many papers study the effect of spatial modulation of the periodic potential (or the presence of impurities in the system) on the dynamics of the SGE solitons [1, 2, 5, 6, 19, 28]. The SGE model with impurities describes, for example, a multilayer ferro-magnet [40-43]. The importance of impurity modes in the kink dynamics was shown in [5, 28, 36]. The structure and properties of localized nonlinear waves excited on an impurity were analyzed numerically in [21, 23, 32]. The case of several point impurities that are of interest in some physical applications was considered in [44] and even the case of spatially modulated harmonic potential was considered [31]. The possibility of exciting localized waves as a result of the kink-impurity interaction has been also investigated in some cases [5, 6]. At first, a simple classical particle model was used for the theoretical description of the kink-impurity interaction when the excitation of an impurity mode was neglected, i.e., a localized oscillatory state on the impurity [5]. Accounting for the excitation of impurity modes in the kink-impurity interaction leads to an interesting effect, i.e., the reflection of a kink by an attracting impurity due to the resonance energy exchange between the translational mode of the kink and the impurity mode. The work [22] has shown that this effect can be observed under certain conditions even in the presence of damping that always occurs in real physical systems. The value of energy consumed for the excitation of localized waves by the kink will determine the value of a new contribution to the effective damping of a moving kink [45]. It may change depending on its initial velocity. It is known that, in this case, the frequency of excited localized nonlinear waves is practically
independent of the initial kink velocity. But this velocity must influence the amplitude of the excited wave. In this paper, we study the problem of the possibility of controlling the structure, dynamics, and energy of localized nonlinear waves excited by a moving kink on an impurity.
2. Basic equations and results
Consider a system determined by the Lagrangian
l= f f H^O2" 2"[1"(l"cosu)+4/?cos ^)dx' (2-l)
The Rayleigh dissipation function has the form
'' 1 f dU
R= J 2aU> ^
— oo
The equation of motion for the scalar field u (x, t) taking into account damping in the system takes the following form:
d2u d2u u du
—r - ttt + [1 - sin u + 2/?.sin - + a— = 0. (2.3)
dt2 dx2 2 dt
where 5 (x) is the Dirac delta function, e is the constant, h is a value of the external force, a is a parameter of damping. Equation (2.3) is the modified sine-Gordon equation (MSGE).
In the case e = a = h = 0 Eq. (2.3) has a solution in the form of a topological soliton (or kink):
u(x,t) = 4arctan (exp [A(u0)_1(x — v01)]) , (2.4)
where A(u) = (1 — u2)1/2, u0 is a continuous parameter (0 < u0 < 1) determining the velocity of a kink. There is also a solution (2.3) in the form of a spatially localized soliton, namely, a quiescent breather:
. . , \/l — w2 sinwi , Ubreather(x,t) = 4 arctan---r -T , (2.5)
w
cosh
a/1 — U)'2(x — Xq)
where w is the breather frequency and x0 is the coordinate of its center.
In the research [5] the case of a "point impurity" was studied in detail (e5(x)) and it was shown that for the corresponding sign of the constant e it acts as an attracting potential on a kink, therefore, soliton can be localized. In the case of a spatially extended impurity the kink-impurity interaction was also studied in both undeformable and deformable kink models [41, 46, 47].
First, we consider the simplest case of point impurity. In this case, for e, h, a ^ 1, using the method of collective coordinate for the coordinate of the center X = X (t) and the amplitude of localized wave a = a (t) the following system of equations has been obtained (see, e.g., [22]):
8X (t) + U' (X) — ea (t) F' (X) = — 8X (t) a + 8h,
a (t) + Q2a (t) - f F (X) = -a (t) a,
where
F(X) = 2 SiDh/Yffl , U{X) = 8- 2t"
'cosh2 X (t)' cosh2 X (t)'
and the frequency of "impurity mode" is determined in the form
F2
Q2 = 1-T (2.7)
The system of equations (2.6) is a system of second-order ordinary differential equations. It is much easier to study than the original equation (2.3). The first equation of (2.6) describes the dynamics of the kink taking into account the presence of the impurity and localized wave on it. The second equation of (2.6) describes the dynamics of a localized wave on an impurity taking into account the interaction of the kink.
2.1. Nondissipative case
Consider the problem of the generation of the impurity mode of maximum possible amplitude for the given parameters of the system. It is known that, when a kink is scattered by an impurity, part of its energy is spent on the excitation of a nonlinear wave localized on an impurity or a "breather" (if the kink passes through the impurity) [48, 49]. The value of this energy can vary depending on the initial kink velocity u0. Since the frequency of the excited wave is practically independent of u0 [50], u0 should influence its amplitude. Consider the simplest case, i.e., the case of a point impurity in the framework of the model (2.6) for h = a = 0. In Fig. 1 the result of the integration is given (2.6) with different modeling parameters, whence it is seen that the breather oscillation amplitude can vary significantly. Suppose that in this model the motion of a kink occurs according to the known law in the form
Xi(t) = uot - 10. (2.8)
If we consider the influence of impurities and neglect the only excitation of localized wave by analogy with solution (2.14), the law of motion of a kink can be represented as
X2 (t) = arcsinh
^+e/2.sinh(t*i-10)
Uo
(2.9)
0.8 n
0.4e 0.0 -0.4 H -0.8
0
(a)
200
400
(b)
(c)
Fig. 1. Time dependence of a(t) calculated by numerical integration (2.6) with initial conditions: X(0) = -10, X(0) = 0.7, a(0) = 0, a(0) = 0 and parameters (a) e = 0.3, (b) e = 0.5, (c) e = 0.7.
(a)
(b)
-10
-10
25
50 t
75 100
25
50 t
75 100
(c)
(d)
100
100
Fig. 2. X(t) = Xi(t) (a) and X(t) = X2(t) (b) dependencies from expressions (2.8) and (2.9), respectively, and dependencies F(t) at X(t) = X1(t) (c) and X(t) = X2(t) (d) constructed from the expression (2.11). The parameters are e = 0.7: 1) v0 = 0.2, 2) v0 = 0.5, 3) v0 = 0.7.
Figures 2a, 2b presents the graph of dependence X(t) corresponding to (2.8)-(2.9) which shows that these dependences are slightly different for large values of u0 (see curves 3 in Figs. 2a, 2b). However, in the region of small values of u0 the dependence (2.9) shows a better qualitative correspondence with the results of numerical simulation (see curves 1 in Figs. 2a, 2b).
Since the collective coordinate X(t) is described by the known law, the system (2.6) reduces to one equation for a(t):
a + Q2a = F (t), (2.10)
where
F (t) = e:
2 sinh (X(t)) cosh2 (X(t))
(2.11)
Then this task can be considered as the excitation of a harmonic oscillator by an external force F(t). Figures 2c, 2d present the dependence F(t) for some particular cases. In this case, the initial conditions are assumed as follows:
X (0) = —10, X(0) = U0, a(0) = 0, a(0) = 0.
We use the Runge-Kutta method of numerical integration to study (2.8)-(2.11). The obtained time evolution of a(t) is presented in Fig. 3, which shows that the oscillation amplitude a(t) differs
0.81 0.4 e 0.0 -0.4 -0.8 H
0
(a)
25
50
t
75 100
(c)
0.81 0.4 e 0.0 -0.4 -0.8 H
(b)
0
25
50 t
(d)
75 100
Fig. 3. Time dependence of a(t) calculated by numerical integration (2.10), for the law of motion of a kink (2.8) (a), (c), (e) and (2.9) — (b), (d), (f) and parameters e = 0.7: (a), (b) u0 = 0.2, (c), (d) uo = 0.5, (e), (f) uo = 0.7.
significantly in the case of small u0 (Figs. 3a, 3b). We also calculate the maximum breather oscillation amplitude Amax for cases (2.8) and (2.9) and compare it with the values obtained by numerical integration of the original system (2.6). Figure 4 shows that for u0 — 1 when the energy of the excited breather is considerably less than the total energy of a kink, both laws of motion of a kink (2.8) and (2.9) represent well the calculations of the original system (2.6). However, for small u0 when the laws of motion differ significantly (see Figs. 2a, 2b, curve 1), the law of motion X(t) = X2(t) (2.9) gives considerably more accurate results. Therefore, it is used in further calculations.
It is possible to find an analytical expression for the breather oscillation amplitude in the case of a point defect using the law of conservation of energy. Then the localized oscillatory mode is given as in [6]:
u(x, t) = B(t) exp(-e\x\/2), (2.12)
where B(t) = a0 cos(Qt + d0), Q is a frequency of the impurity mode (2.7), d0 being an initial phase. Then the energy stored by the breather can be calculated [5, 51]:
E*n = \+j{ + [1 dx = (2.13)
(a)
0.3 0.6 0.9
vo
(b) (c)
0.6 3 0.6 2 I
0.4 S 0.4-B
0.2 0.2 1
0.0 J y 0.0 J , , ,
0.3 0.6 0.9
vq
0.3 0.6 0.9
vo
Fig. 4. Maximum breather oscillation amplitude Amax depending on the initial kink velocity v0 at e = 0.3 (a), e = 0.5 (b) and e = 0.7 (c), calculated by numerical integration: 1) Eqs. (2.10) for the law of motion of a kink (2.8), 2) Eqs. (2.10) for the law of motion of a kink (2.9), 3) systems of equations (2.6).
Suppose that the variation X(t) is described by the solution [6]:
X (t) = arcsinh(7 sinh(f0t)),
(2.14)
where 7 = y/l + ej(2uq). Then the force acting on the breather (oscillator) is given by the expression (taking into account (2.11)):
f (t) =
e2 Asinh (vot) 1 + A2sinh (v0t)
(2.15)
If the complex variable £(t) = à + ifit is used, then the equation of oscillations (2.10) can be reduced to the form
£ - ifi£(t) = f (t), (2.16)
which has the solution [51]:
£(t) = e
int
f(t)e-inTdT
(2.17)
with initial conditions a (-x>) = d(-x>) = 0, i.e., the oscillator (breather) is absent before the interaction. The total energy transferred from the particle (kink) to the oscillator (breather) can be found as in [51]:
Eim = 2n2 e2 sinh2
Q
■ arccos
2"o ~ g
2 Un +1
cosh
2
Ott 2Üö
(2.18)
Then we obtain the following expression from (2.13) and (2.18) for the breather amplitude:
a0
2tr2e3
sinh2
Q
- arccos , 0
2uo V 2u0 + e
2v2- e
cosh-2
/fin
\2uo
(2.19)
The expression (2.19) describes well the results obtained by numerical integration of the system (2.6) (Fig. 5). However, the dependences obtained are very far from the breather amplitude of the model (2.3). It was calculated by the numerical method described in [21]. Nevertheless,
t
(a)
(b)
0.6
|0.4
0.2 0.0
0.6-
¡0.4^
0.2 H 0.0
^ 3
0.6-1 |0.4 0.2 0.0
(c) o
3
□□□□□□□□
0.3 0.6 0.9
0.3 0.6 0.9
vq
0.3 0.6 0.9
Fig. 5. Maximum breather oscillation amplitude Amax depending on the initial kink velocity v0 at e = = 0.3 (a), e = 0.5 (b) and e = 0.7 (c), calculated: 1) by analytical expressions (2.19), 2) by numerical integration of the system of equations (2.6), 3) numerical simulation of Eq. (2.3).
(a)
(b)
0.6-
0.5-
T?
C3 0.4-
X « 0.3-
H
0.2-
o.i-
0.0 J
0.0 0.2
0.4
0.6 0.8
0.60.5-S 0.4-¡0.3-^ 0.20.10.0
0.0 0.2
0.4
vq
0.6 0.8
Fig. 6. Maximum breather amplitude Amax in the center of the extended impurity depending on the initial kink velocity u0 at W =1 (a) and W =1.5 (b) in the framework of the model (2.20). Curves 1-6 correspond to the cases AK = 0.5, 0.75, 1, 1.25, 1.5, 1.75, respectively.
the nature of the dependences is qualitatively the same: one maximum is present in all the curves in Fig. 5, which is presumably related to the ratio of the shape of the external force F(t) and the shape of the potential U(X). Moreover, the maximum point depends on the parameters of the impurity. Similar behavior is observed for the case of extended impurities of rectangular form with a width W and a depth AK, (AK = 0 at x < 0, x > W):
d2u d2u r . _ u du ,
~ q-2 +11 - AR]+2/?sm 2 + am = (2"20)
Figure 6 shows the maximum amplitude of the generated breather obtained numerically from Eq. (2.20). It is also shown that all curves have one maximum. The change in amplitude in this case is more significant than in the case of a point impurity.
2.2. Interaction of a kink with an impurity mode taking into account dissipation and external force
Consider the solution of the system of equations (2.6) in the presence of damping and external force. Since parameter h is absent in the second equation of the system (2.6) (due to the chosen approximation for the effective Lagrangian), the action of the external force does
0.8 0.40
-0.4-1 -0.8
(a)
0
100 200
t
(b)
100 200
t
(c)
100 200
t
Fig. 7. Time evolution of the collective coordinate a(t) (curve 1) calculated from (2.6) and 0(x*) (curve 2) calculated from (2.3), for the case where the kink crosses by inertia the impurity with center at the point x* = 0. The simulation parameters are v0 = 0.7, a = 0.02: (a) e = 0.3, (b) e = 0.5, (c) e = 0.7.
not affect the oscillations of the impurity mode within this analytical model. Therefore, in this section, we restrict ourselves to the assumption that the kink moving by inertia passes through the impurity region.
Figure 7 presents the result of the integration of the analytical model (2.6) (curves 1) and the initial equation (2.3) (curves 2) for three different values of e. It shows that the oscillation amplitude decreases under the effect of damping and the curves 1 and 2 are qualitatively the same. However, a quantitative comparison of their amplitude varies several times with increasing parameter e.
Analysis of the analytical model (in the absence of dissipation) has shown that the excitation of an impurity mode as a result of the scattering of a kink on an impurity can be considered as the excitation of an effective oscillator by a short-term impulse F(t). Figures 2c, 2d present an example of an approximate view of the dependence F(t), where the graphs of this dependence are given for particular cases of the laws of motion of a kink X(t). Thus, after the impurity mode is excited F (t) ^ 0 further its oscillations occur under the influence of damping in the absence of external influence. Then the second equation of the system (2.6) can be represented (in our case) in the form
a(t) + Q2a (t) = -a(t) a, (2.21)
which is a well-known equation of oscillations of a harmonic oscillator with damping [52], its solution can be represented in the form
a(t) = B(t) cos(fii + <po), B(t) = Boexp^-^y (2.22)
where is an initial phase and B(t) is an amplitude of the oscillations. Then the logarithmic
decrement of damping A of the dependence a(t) taking into account the expression Q2 = 1 — can be found as follows:
an an .
A = — = . 2.23
ft \/l —t2/4
The formulas (2.22)-(2.23) describe sufficiently accurately the dependence a(t) only from a particular moment of time. Figure 8 shows that the oscillation amplitude of the dependence a(t) (solid curves) coincides quantitatively with the analytical expression for B(t) (dashed curves), starting from the third oscillation period, since the effect of a nonzero value of the
(a)
200
(b)
200
Fig. 8. Curve 1 is a time evolution of collective coordinate of a(t) from (2.6) for the case where the kink crosses by inertia the impurity at u0 = 0.7, a = 0.02: a) e = 0.5, b) e = 0.7. Curve 2 is a dependence of B(t) calculated by the formula (2.22) at a = 0.02: (a) B0 = 0.38169, (b) B0 = 0.79306.
impulse F(t) is still manifested in the first two periods. Therefore, the initial part of the dependence should be excluded from the analysis when calculating the decrement of damping directly (for example, for the model (2.3)).
Next calculate the logarithmic decrement of damping A of the oscillations of the impurity mode within the framework of the initial equation (2.3). Figure 9 presents the dependence A(e), calculated by integration of Eq. (2.3) (solid curve 1) and the analytical dependence, constructed using the formula (2.23) (dashed curve 2). This figure shows that the curves obtained have a good quantitative agreement. The error of the curve is due to the fact that the amplitude of the excited impurity mode is extremely low at small values of e (see Fig. 7). Therefore, the relative error in calculating the amplitude of the impurity mode is significant.
Thus, the effect of damping on the impurity mode is almost identical to the model (2.3) despite admitted approximation used in the construction of analytical models of (2.6). This suggests that the dissipation of the system, which takes into account the waste of energy on the excitation of the impurity mode, is one of the main damping channels (at least for the point impurities under consideration). However, it is possible that the radiation emitted by a localized wave can also play a significant role for large oscillation amplitudes of the impurity modes. Similar assumptions are discussed, for example, in [50, 53].
1.6 1.2
X
11 / >
1 t 1 \ X',V>—A,
Fig. 9. Logarithmic decrement of damping A depending on e when a = 0.02, vq = 0.7. Curve 1 is calculated in the integration of Eq. (2.3). Curve 2 is built according to the formula (2.23).
3 0.!
<o
0.
0.0
0.1 0.3 0.5 0.7 0.9
£
Fig. 10. Dependence of the relative deviation of the oscillation frequency of the impurity mode for a point impurity on the parameter e at a = 0.02. The deviation is calculated from the exact value, which is determined by the formula (2.7). The impurity mode is excited as the kink moves by inertia (curve 1) under the action of an external force (curve 2) h = 0.0196039 (uo = 0.7).
The significant difference in the amplitude of the excited impurity modes in both models (Fig. 7) is due to the approximations allowed in the derivation of (2.3).
The effect of an external force can have a certain effect on the oscillations of the impurity mode as it was shown in [54]. For example, the action of the field can lead to a change in the oscillation frequency, however, this change may be insignificant in the case of a small external force. Two graphs of the dependence of the relative error in the oscillation frequency of the impurity mode for the cases of motion of the kink under the action of an external force (curve 1) and by inertia (curve 2) are shown for comparison in Fig. 10. The occurrence of damping did not affect considerably the deviation of the frequency of the impurity mode from the nondissipative case (in both cases < 1%). Also, Fig. 10 shows that the relative deviation of the frequency of the impurity mode from the nondissipative case does not exceed 2% for the case of motion of the kink under the action of an external force h = 0.0196039. Thus, the presence of a small external force has no significant effect on the oscillations of the impurity mode. The approximation in the analytical model, when h is absent in the second equation of (2.6), is justified.
3. Conclusion
It was shown that the scattering of a kink on an impurity can be considered as an effective way to excite localized nonlinear waves. The dependence of the breather amplitude on the initial kink velocity when there is only one maximum was calculated for both point and extended impurities. This allows controlling the excitation energy of localized waves using the initial kink velocity and impurity parameters. This behavior remains qualitatively the same in the presence of a small external force and damping. The analytical expression for the decrement of damping was calculated. It is shown that the considered small values of the external force do not have a significant effect on the oscillations of localized waves. This opens up the possibility of using this method of generation of localized impurities in real physical systems. A small external force can be used to set the initial kink velocity.
References
[1] The Sine-Gordon Model and Its Applications: From Pendula and Josephson Junctions to Gravity and High-Energy Physics, J. Cuevas-Maraver, P. Kevrekidis, F.Williams (Eds.), Cham: Springer, 2014.
[2] Encyclopedia of Nonlinear Science, A.Scott (Ed.), New York: Routledge, 2005.
[3] Yakushevich, L.V., Savin, A. V., and Manevitch, L.I., On the Internal Dynamics of Topological Solitons in DNA, Phys. Rev. E, 2002, vol.66, no. 1, 016614, 14pp.
[4] Shamsutdinov, M. A., Khabibullin, I. T., Kharisov, A. T., and Tankeyev, A. P., Dynamics of Magnetic Kinks in Exchange-Coupled Ferromagnetic Layers, Phys. Metals Metallogr., 2009, vol. 108, no. 4, Art. 327; see also: Fiz. Met. i Metalloved, 2009, vol.108, no. 4, pp. 345-358.
[5] Braun, O.M. and Kivshar, Yu. S., The Frenkel-Kontorova Model: Concepts, Methods, and Applications, Berlin: Springer, 2004.
[6] Dauxois, Th. and Peyrard, M., Physics of Solitons, Cambridge: Cambridge Univ. Press, 2010.
[7] Gani, V. A., Lensky, V., and Lizunova, M. A., Kink Excitation Spectra in the (1 + 1)-Dimensional y8 Model, J. High Energ. Phys., 2015, vol. 2015, no. 8, Art. 147, 21 pp.
[8] Gani, V. A., Lizunova, M. A., and Radomskiy, R. V., Scalar Triplet on a Domain Wall, J. Phys. Conf. Ser, 2016, vol.675, 012020, 4pp.
[9] Gani, V. A., Kudryavtsev, A. E., and Lizunova, M. A., Kink Interactions in the (1 + 1)-Dimensional y6 Model, Phys. Rev. D, 2014, vol.89, no. 12, 125009, 12pp.
[10] Kavitha, L., Parasuraman, E., Gopi, D., Prabhu, A., and Vicencio, R. A., Nonlinear Nano-Scale Localized Breather Modes in a Discrete Weak Ferromagnetic Spin Lattice, J. Magn. Magn. Mater.,
2016, vol.401, pp. 394-405.
[11] Kavitha, L., Mohamadou, A., Parasuraman, E., Gopi, D., Akila, N., and Prabhu, A., Modulational Instability and Nano-Scale Energy Localization in Ferromagnetic Spin Chain with Higher Order Dispersive Interactions, J. Magn. Magn. Mater., 2016, vol.404, pp. 91-118.
[12] Golovchan, A. V., Kruglyak, V. V., Tkachenko, V. S., and Kuchko, A. N., Magnonic Band Spectrum of Spin Waves in an Elliptical Helix, R. Soc. Open Sci., 2018, vol.5, no. 1, 172285, 9pp.
[13] Askari, A., Saadatmand, D., and Javidan, K., Collective Coordinate System in (2 + 1) Dimensions: CP1 Lumps-Potential Interaction, Waves Random Complex Media, 2019, vol. 29, no. 2, pp. 368-381.
[14] Moradi Marjaneh, A., Askari, A., Saadatmand, D., and Dmitriev, S. V., Extreme Values of Elastic Strain and Energy in Sine-Gordon Multi-Kink Collisions, Eur. Phys. J. B, 2018, vol. 91, no. 1, Paper No. 22, 8 pp.
[15] Popov, S.P., Compactons and Riemann Waves of an Extended Modified Korteweg-de Vries Equation with Nonlinear Dispersion, Comput. Math. Math. Phys., 2018, vol. 58, no. 3, pp. 437-448.
[16] Popov, S. P., Nonautonomous Soliton Solutions of the Modified Korteweg-de Vries-Sine-Gordon Equation, Comput. Math. Math. Phys., 2016, vol.56, no. 11, pp. 1929-1937.
[17] Dubovik, M. N., Korzunin, L. G., and Filippov, B. N., Asymmetrical Pinning of Vortex Domain Walls in Ferromagnetic Films in Areas with Increased Saturation Magnetization, Phys. Metals Metallogr., 2015, vol. 116, no. 7, pp. 656-662; see also: Fuz. Met. i Metalloved, 2015, vol. 116, no. 7, pp. 694-700.
[18] Ekomasov, E.G., Gumerov, A.M., Kudryavtsev, R. V., Dmitriev, S.V., and Nazarov, V. N., Multi-soliton Dynamics in the Sine-Gordon Model with Two Point Impurities, Braz. J. Phys., 2018, vol. 48, no. 6, pp. 576-584.
[19] Nonlinear Science at the Dawn of the 21st Century, P.L.Christiansen, M. P. S0rensen, A.C.Scott (Eds.), Lecture Notes in Phys., vol.542, Berlin: Springer, 2000.
[20] Currie, J. P., Trullinger, S. E., Bishop, A. R., and Krumhandl, J. A., Numerical Simulation of Sine-Gordon Soliton Dynamics in the Presence of Perturbations, Phys. Rev. B, 1977, vol. 15, no. 12, pp. 5567-5580.
[21] Ekomasov, E. G., Gumerov, A. M., and Murtazin, R. R., Interaction of Sine-Gordon Solitons in the Model with Attracting Impurities, Math. Models Methods Appl. Sci., 2016, vol.40, no. 17, pp.6178-6186.
[22] Ekomasov, E. G., Gumerov, A. M., and Kudryavtsev, R. V., On the Possibility of the Observation of the Resonance Interaction between Kinks of the Sine-Gordon Equation and Localized Waves in Real Physical Systems, JETP Lett., 2015, vol.101, no. 12, pp. 835-839; see also: Pis'ma v Zh. Eksper. Teoret. Fiz, 2015, vol. 101, no. 12, pp. 935-939.
[23] Ekomasov, E. G., Gumerov, A. M., and Kudryavtsev, R. V., Resonance Dynamics of Kinks in the Sine-Gordon Model with Impurity, External Force and Damping, J. Comput. Appl. Math.,
2017, vol.312, pp. 198-208.
[24] Kivshar, Yu. S., Pelinovsky, D. E., Cretegny, T., and Peyrard, M., Internal Modes of Solitary Waves, Phys. Rev. Lett., 1998, vol. 80, no. 23, pp. 5032-5035.
[25] Kivshar, Yu. S., Malomed, B.A., Zhang, F., and Vazquez, L., Creation of Sine-Gordon Solitons by a Pulse Force, Phys. Rev. B, 1991, vol.43, no. 1, pp. 1098-1109.
[26] Gonzalez, J. A., Bellorin, A., and Guerrero, L. E., Internal Modes of Sine-Gordon Solitons in the Presence of Spatiotemporal Perturbations, Phys. Rev. E (3), 2002, vol.65, no. 6, 065601, 4pp.
[27] Kivshar, Yu. S. and Malomed, B. A., Addendum: Dynamics of Solitons in Nearly Integrable Systems, Rev. Mod.. Phys., 1991, vol.63, no. 1, pp. 211-212.
[28] Javidan, K., Analytical Formulation for Soliton-Potential Dynamics, Phys. Rev. E, 2008, vol. 78, no. 4, 046607, 8 pp.
[29] Chacon, R., Bellorin, A., Guerrero, L.E., and Gonzalez, J. A., Spatiotemporal Chaos in Sine-Gordon Systems Subjected to Wave Fields: Onset and Suppression, Phys. Rev. E, 2008, vol. 77, no. 4, 046212, 4 pp.
[30] Gonzalez, J. A., Bellorin, A., Reyes, L.I., Vasquez, C., and Guerrero, L.E., Geometrical Resonance in Spatiotemporal Systems, Europhys. Lett., 2003, vol. 64, no. 6, pp. 743-749.
[31] Gonzalez, J. A., Cuenda, S., and Sanchez, A., Kink Dynamics in Spatially Inhomogeneous Media: The Role of Internal Modes, Phys. Rev. E (3), 2007, vol.75, no. 3, 036611, 7pp.
[32] Gumerov, A. M., Ekomasov, E. G., Murtazin, R. R., and Nazarov, V. N., Transformation of sine-Gordon solitons in models with variable coefficients and damping, Comput. Math. Math. Phys., 2015, vol.55, no. 4, pp. 628-637.
[33] Gonzalez, J. A., Bellorin, A., Garcia-Nustes, M.A., Guerrero, L.E., Jimenez, S., and Vazquez, L., Arbitrarily Large Numbers of Kink Internal Modes in Inhomogeneous Sine-Gordon Equations, Phys. Lett. A, 2017, vol. 381, no. 24, pp. 1995-1998.
[34] Gonzalez, J. A. Jimenez, S., Bellorin, A., Guerrero, L.E., and Vazquez, L., Internal Degrees of Freedom, Long-Range Interactions and Nonlocal Effects in Perturbed Klein-Gordon Equations, Phys. A, 2012, vol.391, no. 3, pp. 515-527.
[35] Saadatmand, D., Dmitriev, S. V., Borisov, D. I., and Kevrekidis, P. G., Interaction of Sine-Gordon Kinks and Breathers with a Parity-Time-Symmetric Defect, Phys. Rev. E, 2014, vol. 90, no. 5, 052902, 10 pp.
[36] Belova, T. I. and Kudryavtsev, A. E., Solitons and Their Interactions in Classical Field Theory, Physics-Uspekhi, 1997, vol.40, no.4, pp. 359-386; see also: Uspekhi Fiz. Nauk, 1997, vol.167, no.4, pp. 377-406.
[37] Popov, S.P., Influence of Dislocations on Kink Solutions of the Double Sine-Gordon Equation, Comput. Math. Math. Phys., 2013, vol. 53, no. 12, pp. 1891-1899.
[38] Malomed, B. A., Dynamics of Quasi-One-Dimensional Kinks in the Two-Dimensional Sine-Gordon Model, Phys. D, 1991, vol.52, nos. 2-3, pp. 157-170.
[39] Saadatmand, D. and Javidan, K., Collective-Coordinate Analysis of Inhomogeneous Nonlinear Klein-Gordon Field Theory, Braz. J. Phys., 2013, vol.43, nos. 1-2, pp.48-56.
[40] Ekomasov, E. G., Gumerov, A. M., Murtazin, R. R., Kudryavtsev, R. V., Ekomasov, A. E., and Abakumova, N. N., Resonant Dynamics of the Domain Walls in Multilayer Ferromagnetic Structure, Solid State Phenom, 2015, vol. 233-234, pp. 51-54.
[41] Ekomasov, E. G. and Shabalin, M. A., Simulation the Nonlinear Dynamics of Domain Walls in Weak Ferromagnets, Phys. Metals Metallogr, 2006, vol.101, Suppl. 1, pp. S48-S50.
[42] Ekomasov, E.G., Murtazin, R. R., Bogomazova, O.B., and Gumerov, A.M., One-Dimensional Dynamics of Domain Walls in Two-Layer Ferromagnet Structure with Different Parameters of Magnetic Anisotropy and Exchange, J. Magn. Magn. Mater., 2013, vol. 339, pp. 133-137.
[43] Ekomasov, E. G., Murtazin, R. R., Bogomazova, O. B., and Nazarov, V. N., Excitation and Dynamics of Domain Walls in Three-Layer Ferromagnetic Structure with Different Parameters of the Magnetic Anisotropy and Exchange, Mater. Sci. Forum, 2016, vol. 845, pp. 195-198.
[44] Gulevich, D.R. and Kusmartsev, F.V., Perturbation Theory for Localized Solutions of the Sine-Gordon Equation: Decay of a Breather and Pinning by a Microresistor, Phys. Rev. B, 2006, vol. 74, no. 21, 214303, 5 pp.
[45] Gumerov, A. M., Ekomasov, E. G., Kudryavtsev, R. V., and Fakhretdinov, M. I., Localized Magnetic Inhomogeneities Generation on Defects As a New Channel of Damping for a Moving Domain Wall, Letters on Materials, 2018, vol.8, no. 3, pp. 299-304 (Russian).
[46] Paul, D. I., Soliton Theory and the Dynamics of a Ferromagnetic Domain Wall, J. Phys. C, 1979, vol. 12, no. 3, pp. 585-593.
[47] Piette, B. and Zakrzewski, W. J., Scattering of Sine-Gordon Kinks on Potential Wells, J. Phys. A, 2007, vol.40, no. 22, pp. 5995-6010.
[48] Gumerov, A. M. and Ekomasov, E. G., Study of the Effect of Point Defects on the Nonlinear Dynamics of Magnetic Nonuniformities, Letters on Materials, 2013, vol. 3, no. 2, pp. 103-105 (Russian).
[49] Ekomasov, E.G., Gumerov, A.M., and Kudryavtsev, R. V., Dynamics of Localized Magnetic In-homogeneities in the Five-Layer Ferromagnetic Structure, Letters on Materials, 2016, vol. 6, no. 2, pp. 138-140 (Russian).
[50] Ekomasov, E. G., Azamatov, S. A., and Murtazin, R. R., Studying the Nucleation and Evolution of Magnetic Inhomogeneities of the Soliton and Breather Type in Magnetic Materials with Local Inhomogeneities of Anisotropy, Phys. Metals Metallogr., 2008, vol. 105, no. 4, pp. 313-321; see also: Fiz. Met. i Metalloved, 2008, vol.105, no. 4, pp. 341-349.
[51] Zhang, F., Kivshar, Yu. S., and Vazquez, L., Resonant Kink-Impurity Interactions in the Sine-Gordon Model, Phys. Rev. A, 1992, vol.45, no. 8, pp. 6019-6030.
[52] Landa, P. S., Nonlinear Oscillations and Waves in Dynamical ¡Systems, Math. Appl., vol.360, Dordrecht: Springer, 2013.
[53] Goodman, R. H., Holmes, P. J., and Weinstein, M. I., Interaction of Sine-Gordon Kinks with Defects: Phase Space Transport in a Two-Mode Model, Phys. D, 2002, vol. 161, no. 1, pp. 21-44.
[54] Ekomasov, E.G., Murtazin, R. R., and Nazarov, V. N., Excitation of Magnetic Inhomogeneities in Three-Layer Ferromagnetic Structure with Different Parameters of the Magnetic Anisotropy and Exchange, J. Magn. Magn. Mater., 2015, vol.385, pp. 217-221.
