On the soliton solutions of a family of Tzitzeica equations
Corina N. Babalic, Radu Constantinescu, Vladimir S. Gerdjikov

TL;DR
This paper investigates soliton solutions of a family of Tzitzeica equations using dressing and Hirota methods, classifying solutions by eigenvalues and discussing their singularities and transformations.
Contribution
It introduces two types of soliton solutions derived via dressing method, linking eigenvalues to soliton types, and outlines construction of multi-soliton solutions for Tzitzeica equations.
Findings
Two types of soliton solutions identified based on eigenvalues.
Construction method for general multi-soliton solutions provided.
Discussion of singularities and variable transformations included.
Abstract
We analyze several types of soliton solutions to a family of Tzitzeica equations. To this end we use two methods for deriving the soliton solutions: the dressing method and Hirota method. The dressing method allows us to derive two types of soliton solutions. The first type corresponds to a set of 6 symmetrically situated discrete eigenvalues of the Lax operator ; to each soliton of the second type one relates a set of 12 discrete eigenvalues of . We also outline how one can construct general soliton solution containing solitons of first type and solitons of second type, . The possible singularities of the solitons and the effects of change of variables that relate the different members of Tzitzeica family equations are briefly discussed. All equations allow quasi-regular as well as singular soliton solutions.
| T1 | T2 | T3 | T4 | |
|---|---|---|---|---|
| , | T1 | T2 | T3 | T4 |
| , | T2 | T1 | T4 | T3 |
| , | T3 | T4 | T1 | T2 |
| , | T4 | T3 | T2 | T1 |
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.
ON THE SOLITON SOLUTIONS OF A FAMILY OF TZITZEICA EQUATIONS
Corina N. Babalic1,2
Radu Constantinescu2 and Vladimir S. Gerdjikov3
Abstract
We analyze several types of soliton solutions to a family of Tzitzeica equations. To this end we use two methods for deriving the soliton solutions: the dressing method and Hirota method. The dressing method allows us to derive two types of soliton solutions. The first type corresponds to a set of 6 symmetrically situated discrete eigenvalues of the Lax operator ; to each soliton of the second type one relates a set of 12 discrete eigenvalues of . We also outline how one can construct general soliton solution containing solitons of first type and solitons of second type, . The possible singularities of the solitons and the effects of change of variables that relate the different members of Tzitzeica family equations are briefly discussed. All equations allow quasi-regular as well as singular soliton solutions.
MSC: 35Q51, 35Q53, 37K40
Key words: Tzitzeica equations, singular soliton solutions, Zakharov-Shabat dressing method, Hirota method
\comm
Communicated by Metin Gürses
Contents
1 Introduction
In the present paper we continue our investigations of the famous equation due to the Romanian mathematician Gheorghe Tzitzeica 111The name of the famous Romanian mathematician contains the Romanian letter Ţ, which may be spelled as Tz. The factor 2 in eq. (1) can be easily removed, but is kept for historical reasons. , which we call now as Tzitzeica 1 equation [27, 28] and a closely related equation which we call Tzitzeica 2; in what follows we will denote them by Ţ1 and Ţ2. It was initially proposed as an equation describing special surfaces in differential geometry for which the ratio is constant, where is the Gauss curvature of the surface and is the distance from the origin to the tangent plane at the given point. Later on it turned out that the equation has wider importance, being nowadays used as an important evolutionary equation in nonlinear dynamics.
The explicit form of Ţ1 and Ţ2 equations is
[TABLE]
i.e. Ţ1 and Ţ2 have different signs in the right hand sides. The transition between T1 and T2 can be performed by several simple changes of variables (see below), some of which substantially modify the singularity properties of their solutions.
Tzitzeica equations attracted a lot of attention at the end of the ’70-ies when for some time it was believed, that it is the only known equation, allowing a finite number of higher integrals of motion [8]. Soon however, it was proved that in fact, it possesses, like the other soliton equations, an infinite number of integrals of motion [30]. Next it was discovered that the equation has a hidden symmetry, which becomes evident in its Lax representation [22, 21]. This important discovery led Mikhailov to the notion of the reduction group and to the family of two-dimensional Toda field theories (TFT) related to the algebras [21]. Soon after it was established that: i) 2-dimensional TFT can be related to any of the simple Lie algebras [23, 24, 9, 19]; ii) other classes of integrable NLEE may also possess such symmetries [9, 12, 13, 7]; and iii) the expansions over the squared solutions and the theory of their recursion operators can be constructed [18, 29, 15].
In previous papers [4, 5] we presented in the derivation of the soliton solutions of Ţ1. Both versions of Tzitzeica equation allow Lax representation proposed by Mikhailov [22, 21]. This allows one to apply the dressing method of Zakharov-Shabat-Mikhailov [33, 21, 32] for calculating their soliton solutions. In fact all these equations are particular examples of 2-dimensional Toda field theories (TFT) [21, 9, 23, 24, 19]. They all can be solved exactly using the inverse scattering method [10, 31, 16].
In the present paper we start with the analysis of a more general class of equations, which we call Tzitzeica family equations. Their general form is
[TABLE]
where and and are some positive real constants. Obviously equation Ţ1 (resp. equation Ţ2) is obtained from (2) by putting , , (resp. , , ). We will call Ţ3 and Ţ4 the equations
[TABLE]
which follow from (2) with , and , respectively.
The paper is organized as follows. In Section 2 we study a class of changes of variables that interrelate different members of Tzitzeica family. We shall see that Ţ1 – Ţ4 equations allow Lax representations so they can be solved exactly by the inverse scattering method, [22, 6]. In Section 3 the Zakharov-Shabat dressing method [33], adapted to systems with deep reductions, [21, 22] is used to construct their soliton solutions. As a result we derive the soliton solutions of first and second types and analyze their singularities. Indeed, we find that even the simplest one-soliton solutions of first type may have an infinite number of singularities for finite values of . Such singularities are characteristic also for other soliton-type equations, e.g. for Liouville equation [1, 2, 25, 26], for sinh-Gordon equation and others, see e.g. [25, 11, 20] and the references therein. At the same time, using an appropriate change of variables we obtain a solution having singularities at only two points which we call ‘quasi-regular’. In Section 4 we outline how the dressing formalism can be extended to derive the -soliton solution of the considered model with solitons of first type and solitons of second type, . In Section 5 we demonstrate how Hirota method can be applied for deriving the soliton solutions of Tzitzeica eqs. and show that it results compatible with the ones of the dressing method. In Section 6 we briefly outline the spectral properties of the Lax operators . We demonstrate that the resolvent of has pole singularities that coincide with the poles of the dressing factor and its inverse. We end by a discussion and conclusions.
2 Lorentz (Anti-)Invariance in 2-Dimensions
Obviously each of the TFT mentioned above can be viewed as a member of a hierarchy of NLEE which can be solved by applying the ISM to the corresponding Lax operator. However the Lorentz invariance singles out the TFT models from all the other members of NLEE in the hierarchy. Indeed, the TFT models allow changes of variables which may drastically change, as we shall demonstrate below, the properties of the soliton solutions.
2.1 Changes of Variables and the Lorentz (Anti-)Invariance
Let us now consider how simple linear change of variables
[TABLE]
affect the solutions of Tzitzeica eqs. Obviously this transformations have to preserve, up to a sign, which means that
[TABLE]
which is equivalent to the relations:
[TABLE]
These relations are satisfied in two cases
[TABLE]
Here and can be, in general, arbitrary complex numbers. However, below we will consider two cases: i) and – real and ii) and – purely imaginary.
Second class of transformations involves shifts of the field
[TABLE]
where is a real constant and takes the values [math] and . If and then Ţ1 goes into
[TABLE]
and similar expression for the Ţ2 equation for , but with opposite signs for the terms in the right hand side.
If we now choose and then Ţ1 goes into
[TABLE]
which for coincides with Ţ3 equation. We have listed the results of several such transformations in Table 1.
2.2 The Lax Representation of Ţ2 Equation
Since different members of Tzitzeica family are related by changes of variables (see Table 1), then it will be enough to consider the Lax representation and soliton solutions of only one of them, say the second equation in (1) Ţ2. It admits the following Lax representation
[TABLE]
where
[TABLE]
The reductions of the Lax pair for Ţ2 equation are similar but not the same as for the well known Ţ1 eq. [4]:
-reduction
[TABLE]
which restricts , and by
[TABLE]
These conditions are satisfied identically. 2. 2.
First -reduction
[TABLE]
i.e.
[TABLE]
which means that . 3. 3.
Second -reduction
[TABLE]
i.e.
[TABLE]
3 The Dressing Method and Dressing Factors for Ţ2 Equation
Let us start with a Lax representation of the form
[TABLE]
The fundamental solution , known also as the ‘naked’ solution, has as potential the trivial solution of Ţ2 equation:
The ’dressed’ Lax pair, given by (11), admits the "dressed" fundamental solution , with the potential the nontrivial solution .
The fundamental solutions and are related by the dressing factor ,
[TABLE]
which means that must satisfy
[TABLE]
Since both Lax pairs (the dressed (11) and the naked one (19)) satisfy the three reductions, then also the dressing factor must satisfy them
[TABLE]
where is defined by eq. (17).
3.1 One Soliton Solution of First Type
A natural anzatz for the dressing factor with simple poles in is [21]
[TABLE]
where is a degenerate matrix of the form
[TABLE]
The first reduction (22a) on is automatically satisfied by the anzatz (23). The second reduction (22b) leads to
[TABLE]
Here and below is understood modulo and
[TABLE]
In addition we must have
[TABLE]
and
[TABLE]
where again all matrix indices are understood modulo 3. These relations can be rewritten as
[TABLE]
So we can consider with no limitations that and . More specifically we will assume that the vector is real, while the vector has purely imaginary components.
The third reduction (22c) on can be put in the form
[TABLE]
Let us now multiply (29) by , take the limit and take into account eq. (14). This gives
[TABLE]
where
[TABLE]
Thus, taking into account that , – real and , we obtain
[TABLE]
In order to obtain the vectors and in terms of and we first impose the limit in equation (21). We obtain that the residue must satisfy
[TABLE]
Since we find that (33) is satisfied if
[TABLE]
i.e.
[TABLE]
which means that is an eigenfunction for the "dressed" Lax pair , , while is an eigenfunction for the "naked" Lax pair , .
From (19), using direct calculation we obtain
[TABLE]
which means that
[TABLE]
Using the notations
[TABLE]
we obtain the following explicit forms for the components of vector
[TABLE]
where
[TABLE]
For and we can rewrite from (39) as the following real-valued functions
[TABLE]
The components of the vector in (32) become
[TABLE]
In order to obtain the solution of Ţ2 equation we impose the limit in (21) with the result:
[TABLE]
where
[TABLE]
which means that
[TABLE]
or
[TABLE]
After introducing from (41) we obtain the 1-soliton solution of the first type for
[TABLE]
where . We observe that this is not a traveling wave solution. Only if we take the limit we obtain a traveling wave solution of the form
[TABLE]
The solution is singular and it blows up for , . For ( and they are purely imaginary), the solution (48) becomes
[TABLE]
The above solution is also singular and it blows up for , .
Remark 1
It is easy to check, that the real parts of in eqs. (48) and (49) are in fact solutions to Ţ4 equation.
In order to get ‘quasi-regular’ solutions of Ţ2 equation, we can apply the changes of variables with or with . This gives the following solutions expressed in terms of hyperbolic functions
[TABLE]
and
[TABLE]
which are singular at the points for which
[TABLE]
or
[TABLE]
respectively. These solutions have also been found by Mikhailov in [21]. As compared with the previous solutions, that have an infinite number of singularities, these ones have singularities at only two points. That is why we took the liberty to call them ‘quasi-regular’.
3.2 The Singularity Properties of the Soliton Solutions
Here we will discuss the types of singularities of the one-soliton solutions and how they are influenced by the changes of variables. As we already mentioned, the singularities in the soliton solutions are not rare, see [11, 20].
Let us first see how the changes of variables affect the Lax representation (11) and, as a consequence, how they affect the fundamental solution. We will be particularly interested in the properties of the ‘naked’ Lax pair and its fundamental solution . This comes from the fact, that the soliton solution is constructed as a rational function of the elements of .
Let us start with the change of variables . Here the situation is simple; we readily get
[TABLE]
In other words this change of variables leaves invariant the compatibility of the Lax pair, so obviously it will map a solution of Ţ2 into a solution of Ţ2. However now we have to keep in mind, that the change of variables must be extended also to the spectral parameter and, of course to the discrete eigenvalues of : and therefore .
In particular, from eq. (40) we see, that
[TABLE]
i.e. and are invariant with respect to transformations provided
[TABLE]
Now it is a bit more interesting to analyze the changes .
[TABLE]
Let us apply a gauge transformation, i.e. change from to
[TABLE]
where and are defined in eqs. (16) and (17) respectively. This gives us
[TABLE]
So the change is equivalent to interchanging the Lax operators and , which again preserves their compatibility condition. Applied to and these transformations lead to
[TABLE]
Of course, analyzing the fundamental solutions we have to pay attention also whether the parameters and are real or purely imaginary. In addition we have to take into account, that could be purely imaginary as above, but for other cases it could also be real. It is precisely this choice of the parameters , and that may change the singularity properties of the solutions.
3.3 One Soliton Solutions of Second Type
Our anzatz for the dressing factor is
[TABLE]
which obviously satisfies the -reduction and the first -reduction.
In order to find how the components of the vector are expressed in terms of the vector we use the same procedure as in the -poles case. First we rewrite the dressing factor in the following form.
[TABLE]
where
[TABLE]
with
[TABLE]
We insert the dressing factor into The second -reduction, we multiply by , take the limit and obtain
[TABLE]
After direct calculation we obtain
[TABLE]
where
[TABLE]
We rewrite the result in a matrix form
[TABLE]
where
[TABLE]
The result is
[TABLE]
where
[TABLE]
[TABLE]
From the above equations we obtain in terms of
[TABLE]
In this case we choose a general form for the poles: ; without restrictions we can choose and determine the expressions of as
[TABLE]
where
[TABLE]
We determine the 1-soliton solution for the second kind of solitons using exactly the same technique
[TABLE]
4 The Generic -soliton Solution for Ţ2 Equation
Let us consider the dressing factor of the following form
[TABLE]
with complex poles from which are purely imaginary, satisfying the following condition: .
Then we write down the residues as degenerate matrices of the form
[TABLE]
From the second -reduction, , after taking the limit , we obtain algebraic equation for in terms of
[TABLE]
Below, for simplicity, we write down the matrix for
[TABLE]
[TABLE]
with
[TABLE]
In order to obtain the 2-soliton solution of the Tzitzeica equation we take the limit in the equations satisfied by the dressing factor and integrate to get
[TABLE]
The above formulae can be easily generalized for any and .
5 Hirota Method for Building 1-soliton Solution of Ţ2 Equation
There are many methods for deriving the soliton solutions; we have demonstrated two of the most used: the dressing method and the Hirota method [17, 3, 6]. Both methods give the same results both for the kinks and for the breathers.
We build the Hirota bilinear form of Ţ2 eq. using the substitution
[TABLE]
Introducing it into the second equation in (1) and decoupling in the bilinear dispersion relation and the soliton-phase constraint, we obtain the following system
[TABLE]
We impose that:
[TABLE]
where , - the wave number, - the angular frequency.
Using a software for analytical computation like MATHEMATICA, we obtain that
[TABLE]
where the dispersion relation is .
Using the above results our 1-soliton solution for Ţ2 acquires the following form
[TABLE]
This solution coincide with the one obtained by Mikhailov in [21] for . In this very direct manner, Hirota method gives immediately the 1-soliton solution of the first type, which we have obtained also in (50) through the dressing method, as a particular case of a more general form (47).
One can also use the standard Hirota technique to derive -soliton solution of first type each parametrized with real eigenvalue and a vector with . We believe, that using Hirota method one can derive also more complicated cases of one- and -soliton solutions. To this end, however one needs a more complicated ansatz for the functions and which would solve equation (82) but could not be reduced to functions of only.
Of course, the equation (82) can be solved in a more general case, but the only one solution we were able to obtain by now, using the well known ansatz (83), was (84), which corresponds to the soliton solutions of first type. To find and corresponding to the second type of soliton solutions is still an open problem for us and it will be tackled in a next paper. A possible approach could be to start from the second type solitons given by the dressing factor method and, on this basis, to guess the ansatz which should be imposed to obtain and verifying (82).
6 The Spectral Properties of the Dressed Lax Operator
Here we shall demonstrate that each dressing adds to the discrete spectrum of sets of discrete eigenvalues.
In our previous paper we showed that the Lax operator has a set of 6 fundamental analytic solutions. We will denote them by where denotes the number of the sector , i.e. those are the sectors closed by the rays . The dressing factor for solitons of first type (23) obviously has simple poles located at , . The inverse of this dressing factor has also simple poles located at , . Each dressing factor for soliton of second type (59) has 6 simple poles located at and , . The inverse of this dressing factor has also 6 simple poles located at and , .
The FAS can be used to construct the kernel of the resolvent of the Lax operator . In this section by we will denote
[TABLE]
where is a regular FAS and is a dressing factor of general form (74). Let us introduce
[TABLE]
[TABLE]
where is the step-function and , see the table 2.
Then the following theorem holds true [4]:
Theorem 2
Let be a Schwartz-type function and let be the simple zeroes of the dressing factor (74). Then
The functions are analytic for where
[TABLE]
having pole singularities at ; 2. 2.
* is a kernel of a bounded integral operator for ;* 3. 3.
* is uniformly bounded function for and provides a kernel of an unbounded integral operator;* 4. 4.
* satisfy the equation*
[TABLE]
Remark 3
The dressing factor has simple poles located at , and where , and . Its inverse has also poles located , and . In what follows for brevity we will denote them by , for .
It remains to show that the poles of coincide with the poles of the dressing factors and its inverse .
The proof follows immediately from the definition of and from eq. (86), taking into account that the limiting value commutes with the corresponding matrix .
Thus we have established that dressing by the factor , we in fact add to the discrete spectrum of the Lax operator discrete eigenvalues; for they are shown on Figure 1.
7 Conclusions
Shortly before finishing this paper we became aware of the fact, that appropriate combination of changes of variables, considered in Section 2 can take each member of Tzitzeica family (2) into one of its 4 versions that we introduced. Let us demonstrate how this can be done of the equation
[TABLE]
where and are real positive constants. Now we shall use somewhat more general change of variables. First we apply the transformation (8) with and . Then we change , where is also real positive constant taken to be . Easy calculation shows that as a result eq. (91) goes into Ţ2 for . Using in addition Table 1 we can transform each member of Tzitzeica family into Ţ2 and then solve it using the results above.
Let us consider the soliton solutions Tzitzeica eq. in a small neighborhood around the singularities, where tends to . Then the second term in the Ţ2 equation can be neglected and the asymptotically we get
[TABLE]
Similarly if in a small neighborhood around the singularity tends to we have
[TABLE]
In both cases we find equations, equivalent to the Liouville equation. Thus the asymptotical behavior of the solutions of Tzitzeica equation around the singularities must be the same as the singularities of Liouville equation [26].
Acknowledgements
One of us (VSG) is grateful to professor A. V. Mikhailov and professor A. K. Pogrebkov for useful discussions. This work has been supported in part by a joint project between the Bulgarian academy of sciences and the Romanian academy of sciences. One of the authors (CNB) acknowledges the support of the strategic grant POSDRU/159/1.5/S/133255, Project ID 133255 (2014), co-financed by the European Social Fund within the Sectorial Operational Program Human Resources Development 2007-2013, and also the support of the project IDEI, PN-II-ID-PCE-2011-3-0083 (MECTS).
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Arkad’ev V. A., Pogrehkov A. K., and Polivanov M. K., The Inverse Scattering Method Applied to Singular Solutions of Nonlinear Equations. II , Teor Mat. Fiz. 54 (1983) 23-37.
- 2[2] Arkad’ev V. A., Pogrebkov A. K. and Polivanov M. K., Singular Solutions of the Kd V Equation and the Inverse Scattering Method , Zap. Nauch. Sem. Leningr. Otd. Mat. Inst., 133 (1984) 17-37, (In Russian).
- 3[3] Babalic C. N., Carstea A. S., On some new forms of lattice integrable equations , Central European Journal of Physics 12 (2014) 341-347.
- 4[4] Babalic N. C., Constantinescu R., Gerdjikov V. S., On Tzitzeica Equation and Spectral Properties of Related Lax Operators , Balkan Journal of Geometry and Its Applications 19 (2014) 11-22.
- 5[5] Babalic N. C., Constantinescu R., Gerdjikov V. S., Two Soliton Solutions of Tzitzeica Equation , Physics AUC 23 (2014) 36-41.
- 6[6] Babalic N. C., Carstea A. S., Alternative Integrable Discretisation of Korteweg de Vries Equation , Physics AUC 21 (2011) 95-100.
- 7[7] Constantin A., Lenells J. and Ivanov R. I., Inverse Scattering Transform for the Degasperis-Procesi Equation , Nonlinearity, 23 (2010) 2559-2575.
- 8[8] Dodd R. K., Bullough R. K., Polynomial Conserved Densities for the Sine-Gordon Equations , Proc. Roy. Soc. London Ser. A, 352 (1977) 481-503.
