2.1. From Ill-Posedness of the HDE to Its Symmetries
Equation (
1) can be considered naturally as an evolution equation, where, say,
and
play the role of space variables, and
is the time one, i.e., one can consider the Cauchy problem:
where
is some given function. But it is easy to check directly that this problem has two trivial solutions:
and
, i.e., solutions that obey
or
. Thus the initial problem (
23) for Equation (
1) is ill-posed. In References [
12,
13], we resolved this ill-definiteness by assuming the linear growth of the solution at infinity. More exactly, we imposed on
v asymptotic behavior in Equations (
6)–(
8), which excludes the cases
. In Reference [
12], we proved by means of the inverse scattering transform that this condition is also sufficient for solvability of the above mentioned Cauchy problem in the class of rapidly decaying
, as in Equation (
7).
Formulation of the inverse problem in Equations (
13) and (
14) demonstrates a specific property of the Hirota difference equation: for an arbitrary set of pairwise different real constants
,
, one can introduce a corresponding number of independent variables
,
, in a way that with respect to any three different variables function
obeys the same Equation (
1) (see also [
11]). Here, we prove this directly by means of the Lax representation (Equations (
3)–(
5)). In order to explicitly specify the set of independent variables involved in the Lax operators and the HDE, we use notation
for the equation as:
Compatibility of the pair of these equations,
and
, with one common index is equivalent to the HDE (Equation (
1)):
with respect to the independent variables
,
, and
, so we denote this equation by
. Summing up the left-hand (and right-hand) sides of equalities (Equation (
24)), we have identity:
that is valid independently of the validity of these equalities themselves. Thus, every HDE
is condition of compatibility of any two equations of the three:
,
, and
, and the third one is compatible, as well.
For a given
evolution with respect to any other
, where
, is nothing but a discrete symmetry. Action of this symmetry on the dependent variable
is given by means of the corresponding equations
,
, and
(see Equation (
24)). These equations must be mutually compatible and compatible with equations
,
, and
, that generates
itself. As we mentioned above, compatibility of equations
and
gives equation
, and the same is valid for all other cases with one common index in operators
L. Thus, we get that also equations
, and
are satisfied. Next, we have to consider compatibility of the pair, say,
and
, where all indexes are different. This results in equality:
that looks to be an equation with respect to four independent variables. It is not the case because, thanks to Equation (
26), we can write
, so that compatibility of the left hand side (l.h.s.) with
follows from compatibility of
with two equations in the right hand side that is already established. In fact, it is easy to check directly that Equation (
27) =
.
Setting that function
depends on an arbitrary number of independent variables, we have to substitute Equation (
7) by means of the equality:
where
at
m-infinity, and where
are real constants that parametrize linear growth of
at infinity. All these parameters must obey the condition in Equation (
6) for any
. As well,
in Equation (
11) sounds now as:
Formulation of the inverse problem (Equations (
12)–(
14)), where
is given above, demonstrates that in the limit
, the spectral data
and then the Jost solution and function
depend on the sum
only. Thanks to Equation (
28) the same is valid for the function
(see Reference [
13] for more detail). Thus, in this limit,
and corresponding equation
becomes identity, as we discussed above. In References [
14,
15], we mentioned that this specific degeneracy of the HDE enables introduction of the continuous symmetries, i.e., symmetries parametrized by the continuous variables. Let us consider limit
as an example. We introduce function
by means of the limit procedure:
In a generic case, we introduce “time”
by equality:
and the same is definition of
-dependence of
,
, and
. In other words, in all these cases, we write, say,
etc., and consider the first order terms with respect to
. Thus, thanks to Equation (
28), we have that:
and thanks to Equation (
10):
Taking that the
t-dependence appears as the limit procedure for discrete symmetries into account, it is natural to expect that it gives continuous symmetries of the HDE. To show this explicitly, we have to consider corresponding limits of Equations (
24) and (
25). It is clear that
if
,
, as well as
for any
. Taking that
and
become identities in this limit, we define:
so that by Equations (
24) and (
28) we get:
Next, in the first order of
equation
(see Equation (
25)) reduces to:
which is exactly the condition of compatibility of the equation
with equations
,
.
Similarly, let us start with the set of six discrete variables, say,
under the condition that all
,
, are valid. Then we can consider limits
,
. Denoting corresponding “times” (continuous parameters) as
,
, and
, we derive three equations
,
, (see Equation (
35)) that guarantee that dependence on these parameters gives symmetries of the
. On the other side, these symmetries must be mutually compatible, i.e., equations
and
must be compatible for all
if
is fulfilled. Thanks to Equation (
31) this gives:
It is necessary to mention that Equation (
27) in the first order of
under substitution (Equation (
31)) reduces to:
that like Equation (
27) can be written as
by means of notation (Equation (
34)) (cf. Equation (
36)). Next, setting in analogy to Equation (
30)
, we again derive Equation (
37) (up to change
) in the first order of Equation (
38) with respect to
. Notice that Equation (
37) involves four independent variables:
,
,
, and
similarly to Equations (
27) and (
38), but in contrast, it is unclear if Equation (
37) can be written as combination of two three-dimensional equations.
2.2. Symmetries and Integrable Equations
Existence of the symmetries introduced above and their mutual compatibility demonstrates that we can consider them equally with the independent variables of the HDE. To be more exact, let now u and depend on six independent variables: Three discrete variables and three continuous ones , .
We know that choosing discrete variables
,
, and
as independent ones, we get the Hirota difference equation
, i.e., Equation (
1). In this case variables
,
are parameters of the continuous symmetries. Now we choose two discrete and one continuous variables to be independent ones, say,
,
, and
. Zero curvature condition of these three variables is given by compatibility of the equations
and
, so it is given by equation
in Equation (
36). In this case, variables
,
, and
play the role of symmetries of this equation, correspondingly the discrete and continuous ones. Next, following Reference [
15], we choose one discrete and two continuous variables:
,
, and
. This choice of variables determines the choice of
in Equation (
35) as the first operator of the Lax pair. But, in this case, we can use neither
nor
as the second operator of the Lax pair in contrast to the above: We have no shift with respect to the second discrete variable. This shift of the Jost solution can be excluded from
by means of the derivative of the equation
(see Equation (
24)), with respect to
. Thus, we get the second equation of the Lax pair in the form:
where we denoted
. In terms of this evolution,
denotes the shift
in the argument of the function
in correspondence to Equation (
2a), while the upper index of
now denotes only a function, different from
v. Compatibility condition now sounds as:
where the first equation can be considered as the evolution one on the function
, while the second stands as definition of an auxiliary function
. Complimentary variables
,
, and
are parameters of the symmetries of this system.
Finally, we consider the case where all three continuous variables are chosen to be independent variables of some differential equation. In this case, neither equations
nor
can be used in the Lax pair: Shifts of the function
with respect to the discrete variables are not allowed. So, we have to use the same trick as above: To exclude these shifts, we differentiate
, as in Equation (
35), with respect to
,
, that gives
. Then the compatibility, i.e., the zero curvature condition reads as:
Substituting difference in the l.h.s. by means of Equation (
24), we derive exactly three equations, Equations (
20)–(
22), for different choices of
. In other words, we get the Lax representation for the Darboux system. Now discrete variables are parameters of the discrete symmetries, i.e., Darboux transformations of the Darboux system, while from the point of view of the latter equations functions
,
are just different functions that are not obliged to be related by any transformation. In the next section, we briefly consider relation of this approach with the inverse scattering transform.
2.3. Symmetries and the Inverse Scattering Transform
The above consideration proves that the equations HDE in Equation (
1),
in Equation (
36), system Equations (
40) and (
41), and the Darboux system (Equation (
19)) have in common solution
depending on all six variables
,
,
and evolutions with respect to any pair of these variables are compatible. But this does not mean that any solution of one of these equations with respect to the corresponding three variables admits compatible introduction of other three variables, in the way that
obeys other equations of the list. To clarify this point, we briefly consider the relation of these symmetries with the corresponding scattering problems. To proceed, notice that because of Equation (
33), we have to modify relations (Equations (
10) and (
11)) as:
and to impose condition that asymptotically:
so that Equations (
24) and (
35) take the form:
Asymptotically, see Equations (
8) and (
32), where we have that:
The inverse problem for the HDE, given by Equations (
12)–(
14) admits switching on dependence on continuous variables
by replacing
in Equation (
13) with:
where
is defined in Equation (
43). The same is valid for the equation
in Equation (
36), because its linear problem is given by the same equation
as the linear problem of HDE, i.e., the inverse problem is also given by Equations (
12) and (
13) with the above substitution. But situation changes, if we consider the system equation (Equations (
40) and (
41)), where the linear problem (i.e., the Lax operator) is given by equation
in Equation (
35), i.e., Equation (
46) with
in terms of
. Taking the asymptotic behavior in Equation (
47) into account, we consider function
in Equation (
32) as perturbation. So the “bare” equation, i.e., equation on
that corresponds to
, sounds as:
Let function
be the Green’s function of this equation:
Then solution of Equation (
46) with the normalization of Equation (
44) is given by integral equation:
It is easy to see that the Green’s function is given by means of:
that has the only departure from analyticity in the
-plane, given by continuous (with exception to the point
)
-derivative:
where
is a function, explicit form of which is not relevant here. This property of the Green’s function assumes that solution of the integral Equation (
51) also obeys the continuous (with the same exception)
-derivative. So, its inverse problem is also given by Equations (
12) and (
13) and admits introduction of other discrete and continuous variables by means of Equation (
48).
Consideration of the Darboux problem is more involved. Choosing Equation (
3) as the first equation of the Lax pair, we rewrite it, thanks to Equation (
42) as equation on
:
Because of asymptotic behavior in Equation (
47), the constant (corresponding to
) part of the operator equals:
Let us introduce the Green’s function:
where in order to obey Equation (
44), we have to construct the Green’s function that decays when
. This condition is satisfied if:
It is easy to see that besides the nonzero
-derivative with respect to
in the complex domain, continuous with exception to the points
and
, this function has discontinuity on the circle:
Correspondingly to the above discussion, it is natural to expect that the inverse problem for the function
in this case will be different from the one given by Equations (
12) and (
13): The departure from analyticity of the Green’s function (Equation (
57)) means that the inverse problem is a combination of the
-problem (Equation (
13)) and the nonlocal Riemann–Hilbert problem on the circle (Equation (
58)), cf. Reference [
24], where the analogous linear problem was considered in detail. To control both these defects of the Jost solution, we need (besides the scattering data
in Equation (
48)) some function
, where
p is a real parameter and support of this function on the
-plane belongs to the circle (Equation (
58)).
In summary, let
be a solution of the Darboux system parametrized via the inverse problem by means of the two kinds of the scattering data described above. If we want to switch on dependence on discrete variables that are compatible with the original ones, shifts of function
give
with respect to Equation (
2a), and this function obeys the HDE (Equation (
1)) for all
, we have to impose the condition that the scattering data
that control discontinuity of the Jost solution on the circle (Equation (
58)) vanish.