Abstract
We study certain deformations of the integrable sine-Gordon model (DSG). It is found analytically and numerically several towers of infinite number of anomalous charges for soliton solutions possessing a special space–time symmetry. Moreover, it is uncovered exact conserved charges associated to two-solitons with a definite parity under space-reflection symmetry, i.e. kink-kink (odd parity) and kink-antikink (even parity) scatterings with equal and opposite velocities. Moreover, we provide a linear formulation of the modified SG model and a related tower of infinite number of exact non-local conservation laws. We back up our results with extensive numerical simulations for kink-kink, kink-antikink and breather configurations of the Bazeia et al. potential V q w = 64 q 2 tan 2 w 2 1 − sin w 2 q 2 , q ∈ R , which contains the usual SG potential V 2 w = 2 1 − cos 2 w .
Keywords
- quasi-integrability
- solitons
- deformed sine-Gordon
- anomalous charges
- non-local charges
1. Introduction
Solitons can be regarded as isolated waves that travel without loss of energy. The solitons emerge with their velocities and shapes completely unchanged after collision to each other, the only outcome being their phase shifts. The soliton solution is the main feature of the integrable models [1, 2, 3]. However, certain non-linear models in physics, with solitary wave solutions, are not integrable. Recently, certain deformations of integrable models such as the sine-Gordon (SG), nonlinear Schrödinger (NLS), Korteweg-de Vries (KdV) and Toda models have been introduced, such that they exhibit soliton-type solutions with some properties resembling to their counterparts of the truly integrable ones. In this context the so-called quasi-integrability concept has been put forward [4]. These properties have been examined in the frameworks of the anomalous zero-curvature [4, 5, 6, 7] and the Riccati-type pseudo-potential approaches [8, 9, 10], respectively.
The main developments have been focused on the construction of infinite number of quasi-conservation laws which give rise to asymptotically conserved charges, i.e. conserved charges, such that their values vary during the scattering of the solitons only. The main observation in the both approaches to quasi-integrability is that, in general, the conserved charges of the standard integrable systems turn out to be the so-called asymptotically conserved charges in the deformed models. In fact, the exact conservation laws of the usual integrable systems become quasi-conservation laws of the deformed integrable models. The non-homogeneous terms of the quasi-conservation laws are dubbed as anomalies such that they vanish when integrated on the space–time plane, provided that the fields satisfy a special space–time symmetry.
The properties of the soliton-like configurations in the quasi-integrable models are, so far, largely unknown. We summarize the main results. First, the one-soliton sectors exhibit infinite conserved charges. Second, the space–time integration of the anomalies vanish when one-soliton like solutions are located far away from each other. The anomalies are significant around the space–time regions of their interaction. Third, a sufficient condition for the vanishing of the space–time integrated anomalies is that the
In the context of the Riccati-type method there have been shown that the deformed SG, KdV and NLS models [8, 9, 10], respectively, possess linear system formulations and that they exhibit infinite towers of exact non-local conservation laws. The NLS-type, KdV-type and SG-type models share the same importance due to their potential applications, since they are ubiquitous in all areas of nonlinear physics, such as Bose-Einsten condensation and superconductivity [12, 13, 14], soliton gas and soliton turbulence in fluid dynamics [15, 16, 17, 18, 19, 20], the Alice-Bob physics [21, 22] and the understanding of a kind of triality among the gauge theories, integrable models and gravity theories [23].
Here, we discuss the previous results in the field by utilizing a deformed sine-Gordon model. We will introduce the relationship between the space–time parity and asymptotically conserved charges. Next, we clarified on the space-reflection parity related to the linear combination of the dual sets of anomalous quantities. In addition, it is focused on the space-reflection symmetry of some two-soliton solutions of deformed sine-Gordon models. Then one proceeds to construct a tower of exactly conserved charges for each solution possessing a definite space-reflection parity. Lastly, by considering linear combinations of the anomalous conserved charges it is showed, through analytical and numerical methods, that there is a subset of exactly conserved charges.
A modified SG model and the space–time symmetries are presented in the next section. In Section 3, the towers of quasi-conservation laws are presented. In Section 4 our numerical simulations are described. The linear formulation and the non-local conservation laws are discussed in the Riccati-type pseudo-potential approach in Section 5. Finally, in Section 6 we present some conclusions.
2. A deformation of the sine-Gordon model
Let us consider the relativistic field theories in
where
where
So, we introduce the deformation parameter
The model (1) possesses several towers of anomalous charges associated to quasi-conservation laws [4, 8, 11]. In [11] it has been introduced a subset of exactly conserved charges associated to space-reflection eigenstates as kink-antikink, kink-kink and breather configurations, respectively. New types of two sets of dual towers of asymptotically conserved charges have been uncovered [8]. Remarkably, even the usual sine-Gordon models possesses anomalous charges. So far, it is attributed to the space–time symmetry properties of the solitons. Those charges can be relevant in the study of soliton gases and formation of certain structures, such as soliton turbulence, soliton gas dynamics and rogue waves [16].
The quasi-integrability has been introduced for deformed sine-Gordon models such that the field
under the special space–time reflection
defined around a given point
and assume that the scalar field is an eigenstate of the operator
In addition, consider an even potential
Several towers of quasi-conservation laws, with anomaly terms possessing odd parities under (6)–(8), have been found [8, 11]. Next, we consider those quasi-conservation laws and examine their anomalies in view of the symmetries (4)–(5) and (6)–(8), respectively.
3. Quasi-conservation laws of the deformed SG model
We will discuss some of the infinite towers of quasi-conservation laws of the deformed SG model (1).
3.1 First type of towers: The SG-type quasi-conservation laws
The usual SG charges turn out to be the anomalous charges of the DSG. So, one has the infinite set of quasi-conservation laws [4, 11].
where the quantities
In
where the quantities
These towers of quasi-conservation laws reproduce the same polynomial form as in the usual sine-Gordon charge densities. In fact, the anomalies
The importance and the relevance of such a dual construction will become clear below when the linear combinations of the charges in (9) and (10) give rise to infinite towers of exactly conserved charges, provided that the space-integral of the linear combination of the anomaly densities
3.1.1 Space-reflection parity and conserved charges
The above dual sets of quasi-conservation laws are used to construct a sequence of conserved charges and vanishing anomalies. The space-reflection symmetry of some soliton solutions of the deformed SG model will imply the existence of an infinite tower of conserved charges. So, let us examine a linear combination, at each order
with the charges
in which the quantities
Since the theory (1) is invariant under space–time translations one has that the energy momentum tensor is conserved. In fact, one has
where
The first non-trivial anomalies become [11].
Notice that for the SG potential (3) the factor
The properties of the quantities
Let us write the anomalies in terms of the
where we have defined the anomaly density
Following analogous procedure as above one has
where we have defined the anomaly density
The anomalies
By direct construction it has been found new towers of anomalous charges in [8]. In the next subsections we will discuss those charges and anomalies in relation to the symmetry (4) and (5).
3.2 Second type of towers
The quasi-conservation laws [8].
define the asymptotically conserved charges
The dual quasi-conservation laws become
where we have introduced the dual asymptotically conserved charges
The densities of the anomalies
3.3 Third type of towers
Let us define the quasi-conservation laws [8].
where we have introduced the asymptotically conserved charges
The interchange
where we have defined the dual asymptotically conserved charges
Similarly, the densities of the anomalies
The relevant anomalies of the lowest order quasi-conservation laws of the above towers will be simulated below for 2-soliton interactions.
Remarkably, the above charges turn out to be anomalous even for the standard sine-Gordon model. In fact, the relevant 2-soliton solutions have been constructed analytically [4, 11] which possess a definite parity under (4)–(5), such that the odd anomaly densities vanish upon space–time integration. The usual explanation for the appearance of novel anomalous charges in the standard sine-Gordon model is the symmetry argument. The anomalous charges also appear in the standard KdV and its deformations [9].
These charges have been computed for soliton collisions in the treatment of soliton gases and formation of some structures in integrable systems, such as integrable turbulence and rogue waves. In the context of the usual KdV model it has been analyzed the behavior of the statistical moments defined by (see e.g. [16, 17])
4. Numerical simulations
Here we will check numerically the lowest order expressions of the various towers of quasi-conservation laws presented above. For this purpose we will numerically solve the Eq. (1) with the particular deformed potential (2). In the Figures 1 and 2 we plot the kink-kink and kink-antikink collisions, respectively. Moreover, we show the first conserved charges, i.e. the energy and momentum for these field configurations.

Figure 1.
Kink-kink with velocities v 2 = − v 1 = 0.15 and q = 2.01 in (2), for initial (green), collision (blue) and final (red) times. Bottom, the energy ( E ) and momentum ( P ) charges of the kink-kink.

Figure 2.
Kink-antikink with velocities v 2 = − v 1 = 0.15 and q = 2.01 in (2), for initial (green), collision (blue) and final (red) times. Bottom, the energy ( E ) and momentum ( P ) charges of the kink-antikink.
4.1 First non-trivial anomalies of the SG-type quasi-conservation laws
We have checked our results by numerical simulation of the anomalies
So, let us write (11) in the form
where
The simulations of the kink-antikink, kink-kink and breather systems of the deformed SG model will consider, as the initial condition, two analytic solitary wave solutions presented in Eq. (1.2) of [4], located some distance apart and stitched together at the middle point.
4.1.1 Kink-antikink
In the Figures 3 and 4 we show the results for kink-antikink system with velocities

Figure 3.
f + 3 , α + 3 and ∫ dt α + 3 in (18) and (19) for kink-antikink with velocities v 2 = − v 1 = 0.5 and ε = 0.06 . The density figure shows initial t i (green), collision t c (blue) and final t f (red) times of the kink-antikink.

Figure 4.
f − 3 , α − 3 and ∫ dt α − 3 in (20)–(21) for kink-antikink with velocities v 2 = − v 1 = 0.5 and ε = 0.06 . The density figure shows initial t i (green), collision t c (blue) and final t f (red) times of the kink-antikink.
4.1.2 kink-kink
In the Figures 5 and 6 we show the results for kink-kink system with velocities

Figure 5.
f + 3 , α + 3 and ∫ dt α + 3 in (18) and (19) for kink-kink with velocities v 2 = − v 1 = 0.5 and ε = 0.06 . The density figure shows initial t i (green), collision t c (blue) and final t f (red) times of the kink-kink.

Figure 6.
f − 3 , α − 3 and ∫ dt α − 3 in (20) and (21) for kink-kink with velocities v 2 = − v 1 = 0.5 and ε = 0.06 . The density figure shows initial t i (green), collision t c (blue) and final t f (red) times of the kink-kink.
So, one can conclude that for kink-antikink (kink-kink) solution the definite parity related to the space-reflection symmetry is a necessary condition in order to achieve a conserved
The both kink-antikink and kink-kink solitons of the SG model with opposite and different velocities do not possess the required parity symmetry. However, it has been shown that in the center-of-mass reference frame (
4.1.3 Breather: kink-antikink bound state
Figures 7 and 8 show the results for breather (kink-antikink bound state) with

Figure 7.
f + 3 , α + 3 and ∫ dt α + 3 in (18) and (19) for breather with ε = 0.06 . The density is shown for three times t ∈ t f − T 0 t f , T 0 = 7.025 . The long-lived breather for t f ≈ 10 5 .

Figure 8.
f − 3 , α − 3 and ∫ dt α − 3 in (20) and (21) for breather with ε = 0.06 . The density is shown for three times t ∈ t f − T 0 t f , T 0 = 7.025 . The long-lived breather for t f ≈ 10 5 .
4.2 Lowest order anomalies of the second and third types of towers
We will compute the linear combinations of the lowest order anomalies of the second and third types of towers in (22)–(27) and (28)–(33), respectively,
4.2.1 Second and third types of towers and lowest order anomalies
The two anomalies in (35) can be written as
Similarly, the two anomalies in (36) can be written as
Notice that under the space–time reflection transformation (4) and (5), the densities of the above anomalies
Under the space-reflection symmetry (6) and (8), some of the densities of the above anomalies will present odd parities; therefore, they must vanish upon space integration. So, in such cases one can have exact conserved charges. These results will be verified for certain solutions as we will see below in the numerical simulations for the kink-kink and kink-antikink solutions.
Figures 9–12 show the anomalies

Figure 9.
Top: The anomaly densities (37) and (38), respectively, plotted in x − coordinate for three times t i (green), t c (blue) and t f (red). Bottom: The anomalies a ± vs t , for kink-antikink collision shown in Figure 2.

Figure 10.
Top: Anomaly densities (39) and (40), respectively, plotted in x − coordinate for three times t i (green), t c (blue) and t f (red). Bottom: Anomalies γ ± vs t , for kink-antikink shown in Figure 2.

Figure 11.
Top: Anomaly densities of (37) and (38), respectively, plotted in x − coordinate for three successive times t i (green), t c (blue) and t f (red). Bottom figures show the relevant anomalies a ± vs t , for kink-kink shown in Figure 1.

Figure 12.
Top: Anomaly densities of (39) and (40), respectively, plotted in x − coordinate for three successive times t i (green), t c (blue) and t f (red). Bottom: Anomalies γ ± vs t , for kink-kink shown in Figure 1.
These results suggest that the quasi-integrable models set forward in the literature [4, 6, 7], and in particular the model (1), would possess more specific integrability structures, such as an infinite set of exactly conserved charges, and some type of linear formulations for certain deformed potentials. So, in the next section we will tackle the problem of extending the Riccati-type pseudo-potential formalism to the deformed sine-Gordon model (1).
5. Riccati-type pseudo-potentials and non-local conservation laws
The Lax equations and Backlund transformations, as well as the conservation laws for the well-known non-linear evolution equations can be generated from the pseudo-potentials and the properties of the Riccati Equation [25, 26, 27, 28, 29].
So, in the next steps we consider a convenient deformation of the usual pseudo-potential approach to integrable field theories. Let us consider the system of Riccati-type equations
and the next linear first order equation for
The compatibility condition
From the system (41) and (42) one can get a quasi-conservation law
This equation has been used to construct a tower of infinite number of quasi-conservation laws [8]. For the standard SG one has
5.1 Pseudo-potentials and a linear system associated to DSG
In this section we search for a linear system formulation of the DSG model. It is achieved by taking into account the Riccati Eq. (41) and the conservation law (44), as well as the Eq. (43). So, the following system of equations has been proposed as a linear formulation of the deformed SG model [8].
where the auxiliary non-local field
In fact, taking into account the expression for the auxiliary field
with
In (49) the coefficient of the linear term in
5.2 Non-local conservation laws
For non-linear equations, not necessarily integrable, which can be derived from a compatibility condition of an associated linear system with spectral parameter, explicit expressions of local and non-local currents can be obtained (see e.g. [30, 31]). In the non-linear
So, one can construct an infinite set of non-local conserved currents through an inductive procedure. Let us define the currents
where
where
Then one can show by an inductive procedure that the (non-local) currents
The first current conservation law
The next conservation law
The construction of analogous linear systems have been performed for deformations of the KdV and NLS models [9, 10]. The construction of the classical Yangian as a Poisson-Hopf type algebra [34] for those non-local currents is worth to pursue in a future work.
6. Conclusions
Our work presents an in-depth demonstration of the quasi-integrability property of the modified sine-Gordon models and the presence of several towers of infinite number of asymptotically conserved charges for soliton configurations satisfying the space–time symmetry (4) and (5). In addition, it is observed that there exist a subset of towers of infinite number of exactly conserved charges, provided that some two-soliton configurations are eigenstates (even or odd) of the space-reflection symmetry (6)–(8).
Moreover, we have uncovered a linear system formulation (45) of the modified SG model, and an infinite set of exact non-local conservation laws (58) associated to that linear formulation.
The space–time and internal symmetries related to quasi-integrability deserve further investigations, due to their applications in several areas of non-linear science, but we hope that the results reported here have opened new lines of research in the context of the quasi-integrability phenomena.
Acknowledgments
JPRC acknowledges brazilian CAPES for financial support. HB thanks FC-UNI (Lima-Perú) and FC-UNASAM (Huaraz-Perú) for partial support and kind hospitality. BMC and CR thank UNASAM for partial financial support.
Notes
- In the x and t laboratory coordinates: η = t + x 2 , ξ = t − x 2 , ∂ η = ∂ t + ∂ x , ∂ ξ = ∂ t − ∂ x , ∂ η ∂ ξ = ∂ t 2 − ∂ x 2