A New 2 + 1-Dimensional Integrable Variable Coefficient Toda Equation

In this paper, a new integrable variable coefficient Toda equation is proposed by utilizing a generalized version of the dressing method. At the same time, we derive the Lax pair of the new integrable variable coefficient Toda equation. The compatibility condition is given, which insures that the new Toda equation is integrable. To further analyze the character of the Toda equation, we derive one soliton solution of the obtained Toda equation by using separation of variables.

• The Generalized Dressing Method
• Variable Coefficient Toda
• Separation of Variables

1. Introduction

Integrable variable coefficient equations describe the real world in many fields of physical and engineering sciences. Many researchers are devoted to discussing these equations by utilizing different methods ref. [1] - [6]. In ref. [7] [8], Dai and Jeffrey extended the dressing method to a generalized version for solving nonlinear evolution equations associated with matrix spectral problems and variable coefficient cases, in which a key is that variable coefficient dressing operators are transformed to different variable coefficient ones. By using the generalization, we have studied integrable variable coefficient coupled Hirota equation in ref. [9]. In ref. [10] [11], integrable variable coefficient Manakov model and cylindrical NLS equation are discussed in detailed, respectively. In ref. [12], we developed the generalized dressing method to the discrete system and an integrable variable coefficient Toda equation is researched. Recently, the dressing method is extended to a matrix Lax pair for Camassa-Holm equation in ref. [13], in which interactions between soliton and cuspon solutions of the system are studied. The dressing method as nonlinear superposition in Sigma models has been researched by Dimitrios Katsinis et al. in ref. [14]. Multi-lump solutions of KP equation with integrable boundary are discussed in ref. [15] by using the generalized dressing method. Nabelek etal. in ref. [16] studied Kaup-Broer system and derived its solutions.

In the present paper, we extend the generalized dressing method to discrete operators similar to ref. [12]. Through direct calculations, we derive a new integrable variable coefficient Toda equation

χ y y − χ n , t t − Δ n ( n − 1 ) ( e χ n − 1 − χ n − 1 ) = 0, (1.1)

where, the coefficient is related to n, Δ = E − E − 1 . Equation (1.1) is an extension of the well known two-dimensional Toda equation. We will construct one soliton solution of (1.1).

The present paper is organized as follows. In Section 2, we obtain a new integrable variable coefficient Toda equation based on the generalized dressing method. In Section 3, as an application, we derive one soliton solution of (1.1) by utilizing the separation of variables.

2. Integrable Variable Coefficient Toda Equation

In this section, we first summarize the variable coefficient version of the dressing method. We extend the generalized version of the dressing method to discrete systems and derive different integrable cylindrical Toda lattice equations by choosing different operators.

First, we consider three linear differential difference operators ref. [12]

F ( n , m , t , y ) ψ n = ∑ m = − ∞ ∞     F ( n , m , t , y ) ψ m , K + ( n , m , t , y ) ψ n = ∑ m = n ∞     K + ( n , m , t , y ) ψ m , K − ( n , m , t , y ) ψ n = ∑ m = − ∞ n     K − ( n , m , t , y ) ψ m . (2.1)

Similar to the generalized dressing method application to continuous system, we introduce the triangular factorization about the operator “ F ”

I + F = ( I + K + ) − 1 ( I + K − ) , (2.2)

where I is the identity operator, K + ( n , m , t , y ) = 0 for m < n and K − ( n , m , t , y ) = 0 for m > n . It is assumed that

sup ∑ m = n 0 ∞ | K ± ( n , m , t , y ) | ψ m < ∞ ,   sup ∑ m = n 0 ∞ | F ( n , m , t , y ) | ψ m < ∞ ,

for all n 0 > − ∞ . For convenience, we denote F ( n , m , t , y ) = F ( n , m ) , K ± ( n , m , t , y ) = K ± ( n , m ) . The discrete Gelfand-Levitan-Marchenko (GLM) equation can be obtained from (2.2), which reads in ref. [12]

F ( n , m ) + K + ( n , m ) + ∑ s = n ∞     K + ( n , s ) F ( s , m ) = 0. (2.3)

We introduce two differential-difference operators M 1 and M 2 defined by

M 1 = ∂ t + ∂ y − n E , (2.4)

M 2 = ∂ t − ∂ y + n E − 1 , (2.5)

where E is the shift operator of the discrete variable n, defined by E k f ( n ) = f ( n + k ) , k ∈ Z , t and y are continuous variables.

The dressing operators N 1 and N 2 can be derived from the relations

N 1 ( I + K + ( n , m ) ) − ( I + K + ( n , m ) ) M 1 = 0 , (2.6)

N 2 ( I + K + ( n , m ) ) − ( I + K + ( n , m ) ) M 2 = 0. (2.7)

Similar to a theorem ref. [7] for continuous systems, it can be proved that N 1 and N 2 are differential-difference operators. For sake of simplicity, we denote K + ( n , m ) = K ( n , m ) .

We write the dressing operators

N 1 = M 1 + D 1 , (2.8)

N 2 = M 2 + D 2 . (2.9)

Acting on function ϕ n on (2.6) and with aid of (2.8), which is reduced to

M 1 K ( n , m ) φ n + D 1 K ( n , m ) φ n + D 1 φ n − K ( n , m ) M 1 φ n = ∑ m = n ∞     K t ( n , m ) φ m + ∑ m = n ∞     K y ( n , m ) φ m − n ∑ m = n + 1 ∞     K ( n + 1 , m ) φ m         + D 1 ∑ m = n ∞     K ( n , m ) φ m + D 1 φ n + ∑ m = n + 1 ∞     K ( n , m − 1 ) φ m ,

from which, comparing coefficient of φ n , we have

K t ( n , n ) + K y ( n , n ) + D 1 K ( n , n ) + D 1 = 0. (2.10)

Letting D 2 = d 1 E − 1 , with aid of (2.7) and (2.9), we have

M 2 K ( n , m ) φ n + D 2 K ( n , m ) φ n + D 2 φ n − K ( n , m ) M 2 φ n = ∑ m = n ∞     K t ( n , m ) φ m − ∑ m = n ∞     K y ( n , m ) φ m + n ∑ m = n − 1 ∞     K ( n − 1 , m ) φ m         + d 1 ∑ m = n − 1 ∞     K ( n − 1 , m ) φ m + d 1 φ n − 1 − n ∑ m = n − 1 ∞ ( m + 1 ) K ( n , m + 1 ) φ m ,

from which, comparing coefficient of φ n − 1 , we have

n K ( n − 1, n − 1 ) − n K ( n , n ) + d 1 K ( n − 1, n − 1 ) + d 1 = 0, (2.11)

and we derive

d 1 = n K ( n , n ) − K ( n − 1, n − 1 ) 1 + K ( n − 1, n − 1 ) . (2.12)

The following theorem in ref. [7] is an extension of original dressing method, which can yield a wide range of integrable variable-coefficient nonlinear evolution equations.

Theorem: If the operators M 1 and M 2 satisfy a relation

[ M 1 , M 2 ] = ρ 1 M 1 + ρ 2 M 2 , (2.13)

where ρ 1 , ρ 2 are arbitrary functions of x , y , n , then their corresponding dressing operators will satisfy the relation

[ N 1 , N 2 ] = ρ 1 N 1 + ρ 2 N 2 . (2.14)

Proof: According to (2.6), (2.7) and (2.13), we can give simple proof as follows through simple calculation. In fact,

[ N 1 , N 2 ] ( I + K + ) = N 1 ( I + K + ) M 2 − N 2 ( I + K + ) M 1 = ( I + K + ) M 1 M 2 − ( I + K + ) M 2 M 1 = ( I + K + ) [ M 1 , M 2 ] = ( ρ 1 N 1 + ρ 2 N 2 ) ( I + K + ) .

Actually, variable-coefficient Toda equations are obtained from (2.14). From (2.14), we derived

d 1 t + d 1 y + ( n + d 1 ) ( 1 − E − 1 ) D 1 = 0, (2.15)

D 1 y − D 1 t − n E d 1 + ( n − 1 ) d 1 = 0. (2.16)

Letting

u n = 1 + K ( n , n ) 1 + K ( n − 1 , n − 1 ) ,   D 1 = v n , (2.17)

then the above Equations (2.15) and (2.16) are reduced to

v n , y − v n , t − Δ n ( n − 1 ) ( u n − 1 ) = 0, (2.18)

u n , y + u n , t + u n ( v n − v n − 1 ) = 0. (2.19)

According to (2.19), we assume that

u n = e χ n − 1 − χ n ,   v n = χ n , t + χ n , y , (2.20)

then (2.18) is reduced to a new integrable variable coefficient Toda equation

χ n , y y − χ n , t t − Δ n ( n − 1 ) ( e χ n − 1 − χ n − 1 ) = 0. (2.21)

Let ξ = y + t , η = y − t , then the above equation is reduced to a new 2 + 1 dimensional Toda lattice equation

4 χ n , ξ η − Δ n ( n − 1 ) ( e χ n − 1 − χ n − 1 ) = 0. (2.22)

The above equations are new and different to classical Toda lattice equation in ref. [17] [18] [19] [20] [21]. Because the coefficient of equation is related to n, this is an important physical meaning.

3. Explicit Solution of Integrable Variable Coefficient Toda Equation

In this section, we shall use the generalized dressing method to construct explicit solutions of the variable coefficient Toda Equation (2.21). Using the relation [ M 1 , F ] = 0, [ M 2 , F ] = 0 , we have

F t ( n , m ) + F y ( n , m ) − n F ( n + 1 , m ) + ( m − 1 ) F ( n , m − 1 ) = 0 , (3.1)

F t ( n , m ) − F y ( n , m ) + n F ( n − 1 , m ) − ( m + 1 ) F ( n , m + 1 ) = 0. (3.2)

Assume that (3.1) and (3.2) have N-soliton solutions in the form of separation of variables

F ( m , n ) = ∑ j = 1 N     f j ( n , t , y ) g j ( m , t , y ) , (3.3)

moveover, we suppose that

K ( m , n ) = ∑ j = 1 N     k j ( n , t , y ) g j ( m , t , y ) . (3.4)

Substituting (3.3) and (3.4) into the GLM (2.3) yields that

K ( n , n ) = − ( f 1 , f 2 , ⋯ , f N ) L − 1 ( g 1 , g 2 , ⋯ , g N ) T , (3.5)

where L is defined by

L j l = δ j l + ∑ s = n ∞     g j ( t , y , s ) f l ( t , y , s ) ,   1 ≤ j , l ≤ N ,

and δ j l is Kronecker’s delta.

In what follows, we will obtain one soliton solution of (2.21). First, we give separation of variables solutions for N = 1 in (3.3) and (3.4),

F ( m , n ) = f 1 g 1 = e p t + q y + n w + η 0 e m w ,   K ( m , n ) = k 1 g 1 = k 1 e m w . (3.6)

From (3.5), we derive

K ( n , n ) = e p t + q y + 2 n w + η 0 − e p t + q y + ( 2 n + 2 ) w + η 0 + e 2 p t + 2 q y + 4 n w + 2 η 0 1 − e 2 w , (3.7)

with p = c h w , q = − s h w , using (2.17), we have

u n = 1 − e 2 w + e p t + q y + 2 n w + η 0 − e p t + q y + ( 2 n + 2 ) w + η 0 + e 2 p t + 2 q y + 4 n w + 2 η 0 1 − e 2 w + e p t + q y + 2 ( n − 1 ) w + η 0 − e p t + q y + 2 n w + η 0 + e 2 p t + 2 q y + 4 ( n − 1 ) w + 2 η 0 . (3.8)

Under transformation u n = e χ n − 1 − χ n , we derive one soliton solution of (2.21)

χ n = χ 0 − ln ( u 1 u 2 ⋯ u n ) . (3.9)

Acknowledgements

The authors thank the authors of the references. The work described in this paper is supported by National Natural Science Foundation of China.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Funding

The work described in this paper is supported by National Natural Science Foundation of China (Grant No.11301149).

Authors’ Contributions

Yanan Huang and Ting Su do derivation and calculations. Junhong Yao mainly draw soliton solution picture.

Conflicts of Interest

There is no competition of interests among authors.

References

1. Tian, B. and Gao, Y.T. (1975) Variable-Coefficient Balancing-Act Method and Variable-Coefficient KdV Equation from Fluid Dynamics and Plasma Physics. The European Physical Journal B, 22, 351-360.

2. Hou, S. and Liu, X.D. (2005) A Numerical Method for Solving Variable Coefficient Elliptic Equation with Interfaces. Journal of Computational Physics, 202, 411-445. https://doi.org/10.1016/j.jcp.2004.07.016

3. Sun, Z.Y., Gao, Y.T., Liu, Y., et al. (2011) Soliton Management for a Variable-Coefficient Modified Korteweg—De Vries Equation. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 84, Article ID: 026606. https://doi.org/10.1103/PhysRevE.84.026606

4. Yan, Z. (1999) Exact Soliton Solutions of the Variable Coefficient KdV-mKdV Equation with Three Arbitrary Functions. Acta Physica Sinica, 48, 1957-1961. https://doi.org/10.7498/aps.48.1957

5. Gngr, F. (2002) Symmetries and Invariant Solutions of the Two-Dimensional Variable Coefficient Burgers Equation. Journal of Physics A General Physics, 34, 4313.

6. Manganaro, N. and Parker, D.F. (1993) Similarity Reductions for Variable-Coefficient Coupled Nonlinear Schr&#246;dinger Equations, Journal of Physics A: Mathematical and General, 26, 4093-4106.

7. Dai, H.H. and Jeffrey, A. (1989) The Inverse Scattering Transforms for Certain Types of Variable-Coefficient KdV Equations. Physics Letters A, 139, 369-372.

8. Jeffrey, A. and Dai, H.H. (1990) On the Application of a Generalizzed Version of the Dressing Method to the Integration of Variable-Coefficient KdV Equation. Rendiconti di Matematica e delle sue Applicazioni, 10, 439-455.

9. Su, T., Dai, H.H. and Geng, X.G. (2009) On the Application of a Generalized Dressing Method to the Integration of Variable-Coefficient Coupled Hirota Equations. Journal of Mathematical Physics, 50, 226.

10. Su, T., Dai, H.H. and Geng, X.G. (2013) A Variable-Coefficient Manakov Model and Its Explicit Solutions through the Generalized Dressing Method. Chinese Physics Letters, 30, Article ID: 060201.

11. Su, T., Ding, G.H. and Fang, J.Y. (2014) Integrable Variable-Coefficient Coupled Cylindrical NLS Equations and Their Explicit Solutions. Acta Mathematicae Applicatae Sinica, 30, 1017-1024. https://doi.org/10.1007/s10255-014-0439-z

12. Dai, H.H. and Su, T. (2010) The Generalized Dressing Method with Applications to the Integration of Variable-Coefficient Toda Equations. Proceedings of the Estonian Academy of Sciences, 59, 293-298.

13. Ivanov, R., Lyons, T. and Orr, N. (2019) Camassa-Holm Cuspons, Solitons and Their Interactions via the Dressing Method. Journal of Nonlinear Science, 30, 225-260.

14. Katsinis, D., Mitsoulas, I. and Pastras, G. (2021) The Dressing Method as Non Linear Superposition in Sigma Models. Journal of High Energy Physics, 2021, 1-26.

15. Dubrovsky, V.G. and Topovsky, A.V. (2020) Multi-Lump Solutions of KP Equation with Integrable Boundary via Dressing Method. Physica D Nonlinear Phenomena, 414, Article ID: 132740.

16. Nabelek, P.V. and Zakharov, V.E. (2020) Solutions to the Kaup? CBroer System and Its (2 + 1) Dimensional Integrable Generalization via the Dressing Method. Physica D Nonlinear Phenomena, 409, Article ID: 132478.

17. Date, E. and Tanaka, S. (1976) Analogue of Inverse Scattering Theory for the Discrete Hill’s Equation and Exact Solutions for the Periodic Toda Lattice. Progress of Theoretical and Experimental Physics, 56, 457-465.

18. Henon M. ,et al. (1974)Integrals of the Toda lattice Physics Review B 9, 1921-1922.

19. Berezansky, Y.M. and Gekhtman, M.I. (1990) Inverse Problem of the Spectral Analysis and Non-Abelian Chains of Nonlinear Equations. Ukrainian Mathematical Journal, 42, 645-658.

20. Ma, W.X. and You, Y.C. (2004) Rational Solutions of the Toda Lattice Equation. Chaos, Solitons, Fractals, 22, 395-406.

21. Ma, W.X. and Maruno, K.I. (2004) Complexiton Solutions of the Toda Lattice Equaiton in Casoration Form. Journal of Physics A, 343, 219-237.