Next Article in Journal
A Rational Approximation of the Two-Term Machin-like Formula for π
Previous Article in Journal
Long- and Medium-Term Financial Strategies on Equities Using Dynamic Bayesian Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

From Algebro Geometric Solutions of the Toda Equation to Sato Formulas

Faculté des Sciences Mirande, Université de Bourgogne Franche-Comté, 2100 Dijon, France
AppliedMath 2024, 4(3), 856-867; https://doi.org/10.3390/appliedmath4030046
Submission received: 26 April 2024 / Revised: 2 July 2024 / Accepted: 2 July 2024 / Published: 9 July 2024

Abstract

:
We know that the degeneracy of solutions to PDEs, given in terms of theta functions on Riemann surfaces, provides important results about particular solutions, as in the case of the NLS equation. Here, we degenerate the so called finite gap solutions of the Toda lattice equation from the general formulation in terms of abelian functions when the gaps tend to points. This degeneracy allows us to recover the Sato formulas without using inverse scattering theory or geometric or representation theoretic methods.

1. Introduction

Different formulations can be chosen to represent the Toda equation. In the original work of 1967 [1], the vibration of a uniform chain of particles with nonlinear interaction was studied, and Toda considered the following equation for the n-th particle in the chain:
m x 2 ( u n ) = r ( ϕ ) ( u n u n 1 ) + r ( ϕ ) ( u n + 1 u n ) ,
where m stands for the mass of the particles, and ϕ ( r ) is the interaction energy between adjacent particles. Some solutions were constructed for particular interactions ϕ .
The same equation was considered with m = 1 by Date and Tanaka in 1976 [2], and they constructed some general solutions but in terms of integrals that are difficult to use.
In the appendix of the paper of Dubrovin [3], published in 1981, Krichever considered the non-abelian version of this equation and provided solutions in terms of the theta function on Riemann surfaces.
In 1995, Matveev and Stahlhoffen [4] considered the following version of the Toda equation:
x 2 ( u n ) = exp ( u n 1 u n ) exp ( u n u n + 1 ) .
They used the Darboux transformation to contruct solutions to this equation in terms of Casoratis.
In 2014, Zhang and Zhou [5] used the following representation of the Toda equation
x   t   ( u n ) = t   ( u n ) + α ( t ) ( u n 1 2 u n + u n + 1 ) .
They used the generalization of the exp-function to contruct multiwave solutions to this equation.
More recently, Sun, Ma, and Yu [5] examined the following representation of the Toda equation:
Δ ( u n ) = 4 ( exp ( u n 1 u n ) exp ( u n u n + 1 ) ) .
They used a logarithmic transformation to obtain some particular solutions in terms of logarithms.
In 2020, Duarte [6] considered the following representation of the Toda equation:
α x 2 ( u ( x , y , n ) ) + β y 2 ( u ( x , y , n ) ) = exp ( u ( x , y , n 1 ) u ( x , y , n ) ) exp ( u ( x , y , n ) u ( x , y , n + 1 ) ) .
He used a particular ansatz combined with the properties of Laplace’s equation to contruct some solutions in terms of trigonometric functions.
Also in 2020, Schiebold and Nilson [7] studied another version of the Toda equation in the form
x 2 y ( log ( 1 + u n ) ) = u n + 1 2 u n + u n + 1 .
They constructed solutions in the frame of linear algebra by means of determinants.
Here, we consider the Toda equation in the following representation [8,9,10,11]:
r n r n + 1 c n ψ n + 1 + v n ψ n + r n r n 1 c n 1 ψ n 1 = λ ψ n t   ψ n = r n 2 r n + 1 c n ψ n + 1 + t   ( ln r n ) ψ n r n 2 r n 1 c n 1 ψ n 1 ,
with
c n = e 2 I 0 θ 2 ( t V + n D + Z ) θ ( t V + ( n + 1 ) D + Z ) θ ( t V + ( n 1 ) D + Z )
v n = R 0 + t ln θ ( t V + ( n 1 ) D + Z ) θ ( t V + n D + Z ) ,
where θ is the classical Riemann theta function.
We know that the degeneracy of the solutions to PDEs given in terms of Riemann theta functions provides some important particular solutions. In the case of the NLS equation, we have managed to construct a quasi-rational solution involving the determinant of order N for each positive integer N depending on 2 N 2 real parameters [12]. In the case of the KdV equation, we constructed solutions in terms of Fredholm determinants and Wronskians [13], from Baker–Akhiezer functions.
In this study, we recovered Sato formulas for the Toda equation using this method.
From finite gap solutions given in terms of Riemann theta functions, we constructed some quasi-rational solutions and recovered the Sato formulas, using degeneracy, as given in the frame of the NLS equation [14].
Precisely, we derived multisoliton solutions from finite gap solutions by a limit transition, i.e., by making gaps tend toward points in a certain Riemann surface. This was accomplished in the spirit of [14] or more recently [13]. One strength of this approach is that it does not rely on inverse scattering theory or geometric and representation theoretic methods, which offers a fresh perspective on the problem.
We consider the Riemann surface Γ represented by k = 1 g a k b k a k 1 b k 1 of the algebraic curve defined by [8,9]
ω 2 = j = 1 2 g + 2 ( z E j ) ,
with E j E k , j k .
Let us consider Ω 1 and Ω 0 [11] abelian integrals, verifying
a k d Ω 1 ( P ) = a k d Ω 0 ( P ) = 0 k = 1 , , g ,
with the following asymptotic behavior for P = ( ω , z ) Γ [11],
Ω 1 ( P ) = ± 1 2 ( z + R 0 + O ( z 1 ) ) , P P ±
Ω 0 ( P ) = ± ( ln z I 0 + O ( z 1 ) ) , P P ± ,
ω = ± ( z g + 1 + O ( z g ) ) , P P ± .
The vectors V , D, X, and Z are defined by
V j = 1 2 π i b j d Ω 1 ,
D = U ( P + ) U ( P ) ,
X = K + j = 1 g U ( P j ) , Z = U ( P + ) X ,
where K is the vector of Riemann constants.
U ( P ) is the classical abelian integral P 0 P d U .
The solutions to the system (6) can be written [10,11] as
ψ n ( P , t ) = r ˜ n ( t ) e n Ω 0 ( P ) + t Ω 1 ( P ) ( θ ( U ( P ) + t V + n D X ) θ ( U ( P ) X ) ) ,
with
r ˜ n ( t ) = r n ( t ) θ ( Z ) θ ( Z D ) θ ( t V + n D + Z ) θ ( t V + ( n 1 ) D + Z ) 1 2 .

2. Degeneracy of Solutions

Let us suppose that E j are real, E m < E j if m < j , and try to evaluate the limits of all objects in Formula (16) when E 2 m , E 2 m 1 tends to α m , α m = κ m 2 , κ m > 0 , for 1 m g .
As in the previous section, θ is constructed from the matrix of the B-periods of the surface Γ , the coefficients c j k are related with abelian differentials d U j by
d U j = k = 1 g c j k z g k k = 1 2 g + 2 ( z E k ) d z ,
and the coefficients c j k can be obtained by solving the system of linear equations
a k d U j = δ j k , 1 j g , 1 k g .
In the remainder of this article, we use the following notations:
κ m j = ( α m + E 2 g + 1 ) ( α j + E 2 g + 2 ) , κ j m = ( α j + E 2 g + 1 ) ( α m + E 2 g + 2 ) , C k ( z ) = ( z E 2 g + 1 ) ( α k E 2 g + 2 ) , D k ( z ) = ( z E 2 g + 2 ) ( α k E 2 g + 1 ) , L k = α k + E 2 g + 1 , M k = α k + E 2 g + 2 , F ( z ) = z E 2 g + 1 , G ( z ) = z E 2 g + 2 .

2.1. Limit of P ( z ) = j = 1 2 g + 2 ( z E j )

Now, we study the limit of P ( z ) = j = 1 2 g + 2 ( z E j ) .
The limit of P ( z ) = j = 1 2 g + 2 ( z E j ) is obviously equal to P ˜ ( z ) = j = 1 g ( z + α j ) 2 ( z E 2 g + 1 ) ( z E 2 g + 2 ) or, with (19),
P ˜ ( z ) = j = 1 g ( ( z + α j ) F ( z ) G ( z ) ) 2 .

2.2. Limit of d U m = k = 1 g c m k z g k k = 1 2 g + 2 ( z E k ) d z

Now, we study the limit of d U m = k = 1 g c m k z g k k = 1 2 g + 2 ( z E k ) d z .
The limit of d U m is equal to d U ˜ m = φ m ( z ) j = 1 g ( z + α j ) ( z E 2 g + 1 ) ( z E 2 g + 2 ) d z , where φ m ( z ) = k = 1 g c ˜ m k z g k . The normalization condition takes the form in the limit
a k d U j 2 π i φ j ( α k ) m k   ( α m α k ) ( α k + E 2 g + 1 ) ( α k + E 2 g + 2 ) = δ k j ,
which proves that the numbers α m , m k are the zeros of the polynomials φ k ( z ) ; hence, φ k ( z ) can be written as φ k ( z ) = c ˜ k 1 m k ( z + α m ) .
By (21), we obtain in the limit
c k 1 ˜ = ( α k + E 2 g + 1 ) ( α k + E 2 g + 2 ) 2 π i .
So,
φ k ( z ) = ( α k + E 2 g + 1 ) ( α k + E 2 g + 2 ) 2 π i m k ( z + α m ) .
Moreover,
d U ˜ k ( z ) = φ k ( z ) d z ( z E 2 g + 1 ) ( z E 2 g + 2 ) m = 1 g ( z + α m ) ;
in other words,
d U ˜ k ( z ) = ( α k + E 2 g + 1 ) ( α k + E 2 g + 2 ) 2 π i ( z E 2 g + 1 ) ( z E 2 g + 2 ) ( z + α k ) d z ,
and with (19),
d U ˜ k ( z ) = L k M k 2 π i F ( z ) G ( z ) ( z + α k ) d z .

2.3. Limit of B m k

The subject of this subsection is the study of the limit of B m k .
We have
I = α m E 2 g + 2 d U k 1 2 B ˜ m k .
The integral I can be easily evaluated along the real axis on the upper sheet of surface Γ , and we obtain
B ˜ m k = i π ln ( α m + E 2 g + 1 ) ( α k + E 2 g + 2 ) + ( α k + E 2 g + 1 ) ( α m + E 2 g + 2 ) ( α m + E 2 g + 1 ) ( α k + E 2 g + 2 ) ( α k + E 2 g + 1 ) ( α m + E 2 g + 2 ) ,
or with the previous notations (19),
B ˜ m k = i π ln κ m k + κ k m κ m k κ k m .
So, i B k k tends to .
Moreover, we have
E 2 g + 2 P d U ˜ k = i 2 π ln ( z E 2 g + 1 ) ( α k + E 2 g + 2 ) + ( z E 2 g + 2 ) ( α k + E 2 g + 1 ) ( z E 2 g + 1 ) ( α k + E 2 g + 2 ) ( z E 2 g + 2 ) ( α k + E 2 g + 1 ) ,
or with the previous notations (19),
E 2 g + 2 P d U ˜ k = i 2 π ln C k ( z ) + D k ( z ) C k ( z ) D k ( z ) .

2.4. Limit of Ω 1

In this subsection, we study the limit of Ω 1 .
Ω 1 is an abelian integral of the second kind satisfying the conditions
a k d Ω 1 ( P ) = 0 , 1 k g ,
such that
Ω 1 ( P ) = ± 1 2 ( z + R 0 + O ( z 1 ) ) , P P ±
The limit of d Ω 1 ( P ) is equal to d Ω ˜ 1 ( P ) = ϕ   ( z ) j = 1 g ( z + α j ) ( z E 2 g + 1 ) ( z E 2 g + 2 ) d z , where ϕ   ( z ) = k = 0 g + 1 c ˜ k z g + 1 k .
For Ω 1 ( P ) , satisfying the condition Ω 1 ( P ) = ± 1 2 z + O ( 1 ) when z ± , we have c 0 ˜ = 1 2 . Moreover, the conditions
a k d Ω 1 ( P ) = 0 , 1 k g
prove that α 1 , , α g are the zeros of ϕ ; thus, ϕ ( z ) = 1 2 m = 1 g ( z + α m ) ( z β ) .
We have that P ( z ) = k = 1 2 g + 2 ( z E k ) tends to P ˜ ( z ) = ( z E 2 g + 1 ) ( z E 2 g + 2 ) m = 1 g   ( z + α m ) , and we obtain
d Ω ˜ 1 ( z ) = z 2 ( z E 2 g + 1 ) ( z E 2 g + 2 ) d z ;
Ω 1 ( P ) tends to
Ω ˜ 1 ( P ) = E 2 g + 2 z u 2 ( u E 2 g + 1 ) ( u E 2 g + 2 ) d u .
This can be evaluated, and it gives
Ω ˜ 1 ( P ) = E 2 g + 1 + E 2 g + 2 4 ln z E 2 g + 1 + z E 2 g + 2 z E 2 g + 1 z E 2 g + 2 + 1 2 ( z E 2 g + 1 ) ( z E 2 g + 2 ) .
With the notations defined in (19), it can be rewritten as
Ω ˜ 1 ( P ) = E 2 g + 1 + E 2 g + 2 4 ln F ( z ) + G ( z ) F ( z ) G ( z ) + 1 2 F ( z ) G ( z ) .

2.5. Limit of Ω 0

We consider Ω 0 , and we study its limit.
Ω 0 is an abelian integral of the third kind satisfying the conditions
a k d Ω 0 ( P ) = 0 , 1 k g ,
such that
Ω 0 ( P ) = ± ( ln z I 0 + O ( z 1 ) ) , P P ± .
The limit of d Ω 0 ( P ) is equal to
d Ω ˜ 0 ( P ) = k = 0 g c ˜ k z g k j = 1 g ( z + α j ) ( z E 2 g + 1 ) ( z E 2 g + 2 ) d z .
Ω 0 ( P ) satisfies the condition Ω 0 ( P ) = ± ln ( z ) + O ( 1 ) , and when z ± , we have c 0 ˜ = 1 .
Moreover, the conditions
a k d Ω 0 ( P ) = 0 , 1 k g
prove that α 1 , , α g are the zeros of ϕ ; thus, ϕ ( z ) , defined by ϕ ( z ) = k = 0 g c ˜ k z g k , can be written as ϕ ( z ) = m = 1 g ( z + α m ) .
As P ( z ) = k = 1 2 g + 2 ( z E k ) tends to P ˜ ( z ) = ( z E 2 g + 1 ) ( z E 2 g + 2 ) m = 1 g ( z + α m ) , we obtain
d Ω ˜ 0 ( z ) = 1 ( z E 2 g + 1 ) ( z E 2 g + 2 ) d z ,
and Ω 0 ( P ) tends to Ω ˜ 1 ( P ) = E 2 g + 2 z 1 ( u E 2 g + 1 ) ( u E 2 g + 2 ) d u .
This can be calculated, and it gives
Ω ˜ 0 ( P ) = ln z E 2 g + 1 + z E 2 g + 2 z E 2 g + 1 z E 2 g + 2 .
It can be written with (19) as
Ω ˜ 0 ( P ) = ln F ( z ) + G ( z ) F ( z ) G ( z ) .

2.6. Limit of X

We deduce the limit of X.
From the previous sections,
B ˜ m k = i π ln κ m k + κ k m κ m k κ k m ,
U ˜ k ( P ) = i 2 π ln C k ( z ) + D k ( z ) ) C k ( z ) D k ( z ) .
As K is defined by K m = j = 1 g B j m m 2 , and X = K + j = 1 g U ( P j ) , we can write
X m X m ˜ = i π j = 1 , j m g ln κ m j + κ j m κ m j κ j m m 2 + i 2 π j = 1 g ln C m ( z j ) + D m ( z j ) C m ( z j ) D m ( z j ) + i .

2.7. Limit of U ( P + ) , U ( P ) , D and Z

In this subsection, we study the limit of U ( P + ) , U ( P ) , D, and Z.
From the previous sections, it is easy to obtain the limits of U ˜ k as
lim P P + U ˜ k ( P ) = i 2 π ln α k + E 2 g + 2 + α k + E 2 g + 1 α k + E 2 g + 2 α k + E 2 g + 1 ,
or
lim P P + U ˜ k ( P ) = i 2 π ln M k + L k M k + L k .
D k = U k ( P + ) U k ( P ) D k ˜ = i 2 π ln M k + i L k M k i L k .
Thus,
lim P P U ˜ k ( P ) = i 2 π ln α k + E 2 g + 2 + α k + E 2 g + 1 α k + E 2 g + 2 α k + E 2 g + 1 D k ,
or
lim P P U ˜ k ( P ) = i 2 π ln M k + L k M k + L k + i 2 π ln M k + i L k M k + i L k .
Therefore, the limit of Z = U ( P + ) X is given by
Z k Z k ˜ = i π j = 1 , j k g ln κ k j + κ j k κ k j κ j k + k 2 i 2 π j = 1 , g ln C k ( z j ) + D k ( z j ) C k ( z j ) D k ( z j ) i i 2 π ln M k + L k M k L k

2.8. Limit of V

In this paragraph, the limit of V is studied.
V k = 1 2 π i b k d Ω 1 is defined in (13).
We have I = α k E 2 g + 2 d Ω 1 1 2 b k d Ω ˜ 1 .
As in the determination of the limit of Ω 1 , we have
α k E 2 g + 2 d Ω ˜ 1 ( P ) = E 2 g + 1 + E 2 g + 2 4 ln α k + E 2 g + 1 + α k + E 2 g + 2 α k + E 2 g + 1 α k + E 2 g + 2
1 2 ( α k + E 2 g + 1 ) ( α k + E 2 g + 2 ) .
So,
V k = 1 2 π i b k d Ω 1 V ˜ = E 2 g + 1 + E 2 g + 2 4 π i ln α k + E 2 g + 1 + α k + E 2 g + 2 α k + E 2 g + 1 α k + E 2 g + 2
1 2 π i ( α k + E 2 g + 1 ) ( α k + E 2 g + 2 ) ,
V ˜ k = E 2 g + 1 + E 2 g + 2 4 π i ln L k M k L k + M k 1 2 π i L k M k .

2.9. Limit of θ ( p )

We determine the limit of θ ( p ) .
Let us denote A as the argument of θ ( p ) = k Z g exp { π i ( B k | k ) + 2 π i ( k | p ) } .
A can be rewritten in the form
A = π i j = 1 g B j j k j ( k j 1 ) + 2 π i j > m B m j k m k j + j = 1 g π i ( 2 p j + B j j ) k j .
Using the inequality k j ( k j 1 ) 0 for all k Z g and the fact that i B k k tends to , we can reduce the limit θ ˜ of θ ( p ) to a finite sum taken over vectors k Z g , such that each k j must be equal to 0 or 1.
In this section, we compute the limit of all the terms in the expression of the solution ψ ( P , t ) in θ . We denote p i as the arguments of these different expressions.
We first study the terms of r ˜ n corresponding to the arguments p 1 = Z = U ( P + ) X , p 2 = Z D , p 3 = t V + n D Z , p 4 = t V + ( n 1 ) D Z .
Then, we study the remaining terms of ψ n corresponding to the arguments p 5 = U ( P ) X , p 6 = U ( P ) + t V + n D X .
With these notations, the solution ψ ( n , , t , z ) can be written as
ψ n ( P , z ) = r n ( t ) θ ( p 1 ) θ ( p 2 ) θ ( p 3 ) θ ( p 4 ) 1 2 θ ( p 6 ) θ ( p 5 ) e n Ω 0 ( P ) + t Ω 1 ( P ) .

2.10. Limit of θ ( p 1 )

We denote p 1 as the term Z = U ( P + ) X .
We study its limit.
Then, θ ( p 1 ) θ ˜ 1 :
θ ˜ 1 = k Z g , k j = 0 o r 1 exp j > l 2 ln κ l j κ j l κ l j + κ j l k j k l + j = 1 g k j ln M j + L j M j L j
+ l j 2 ln κ j l + κ l j κ j l κ l j + j π i + l = 1 g ln C j ( z l ) + D j ( z l ) C j ( z l ) D j ( z l ) .

2.11. Limit of θ ( p 2 )

We consider p 2 = Z D . We study its limit.
Then, θ ( p 2 ) θ ˜ 2 . As D j D ˜ j , it is easy to see that
θ ˜ 2 = k Z g , k j = 0 o r 1 exp j > l 2 ln κ l j κ j l κ l j + κ j l k j k l + j = 1 g k j ln M j + L j M j L j
+ l j 2 ln κ j l + κ l j κ j l κ l j + j π i + l = 1 g ln C j ( z l ) + D j ( z l ) C j ( z l ) D j ( z l ) + j = 1 g ln M j + i L j M j i L j .

2.12. Limit of θ ( p 3 )

Let p 3 be t V + n D + Z . We study its limit.
Then, θ ( p 3 ) θ ˜ 3 :
θ ˜ 3 = k Z g , k j = 0 o r 1 exp j > l 2 ln κ l j κ j l κ l j + κ j l k j k l + j = 1 g k j ln L j + M j L j M j
+ t E g + 1 + E 2 g + 2 2 ln L j M j L j + M j L j M j + n ln M j + i L j M j i L j
+ l j 2 ln κ j l + κ l j κ j l κ l j + j π i + l = 1 g ln C j ( z l ) + D j ( z l ) C j ( z l ) D j ( z l ) .

2.13. Limit of θ ( p 4 )

We consider p 4 = t V + ( n 1 ) D + Z ; θ ( p 4 ) θ ˜ 4 . We determine its limit.
As D j D ˜ j , it is easy to verify that
θ ˜ 4 = k Z g , k j = 0 o r 1 exp j > l 2 ln κ l j κ j l κ l j + κ j l k j k l + j = 1 g k j ln L j + M j L j M j
+ t E g + 1 + E 2 g + 2 2 ln L j M j L j + M j + L j M j + ( n 1 ) ln M j + i L j M j i L j
+ l j 2 ln κ j l + κ l j κ j l κ l j + j π i + l = 1 g ln C j ( z l ) + D j ( z l ) C j ( z l ) D j ( z l ) .

2.14. Limit of θ ( p 5 )

The term p 5 = U ( P ) X tends to p ˜ 5 = U ˜ ( P ) X ˜ .
We determine its limit.
Then, θ ( p 5 ) θ ˜ 5 :
θ ˜ 5 = k Z g , k j = 0 o r 1 exp j > l 2 ln κ l j κ j l κ l j + κ j l k j k l + j = 1 g k j ln C j ( z ) D j ( z ) C j ( z ) + D j ( z )
+ l j 2 ln κ j l + κ l j κ j l κ l j + j π i + l = 1 g ln C j ( z l + D j ( z l C j ( z l D j ( z l .

2.15. Limit of θ ( p 6 )

The term p 6 = U ( P ) + t V + n D X tends to p ˜ 6 = U ˜ ( P ) + t V ˜ X ˜ + n D ˜ .
Then, we determine its limit. θ ( p 6 ) θ ˜ 6 :
θ ˜ 6 = k Z g , k j = 0 o r 1 exp j > l 2 ln κ l j κ j l κ l j + κ j l k j k l + j = 1 g k j ln L j + M j L j M j
+ t E g + 1 + E 2 g + 2 2 ln L j M j L j + M j L j M j + n ln M j + i L j M j i L j
+ l j 2 ln κ j l + κ l j κ j l κ l j + j π i l = 1 g ln C j ( z l ) + D j ( z l ) C j ( z l ) D j ( z l ) + ln C j ( z ) D j ( z ) C j ( z ) + D j ( z ) .

2.16. Limit of Solutions Ψ

From the previous section, we can give the limit ψ ˜ of ψ . It takes the following form:
ψ ˜ ( z , n , t ) ) = r n ( t ) θ 1 ˜ θ 2 ˜ 1 2 θ 5 ˜ × e n ln F ( z ) + G ( z ) F ( z ) G ( z ) + t E 2 g + 1 + E 2 g + 2 4 ln F ( z ) + G ( z ) F ( z ) G ( z ) + 1 2 F ( z ) G ( z ) θ 6 ˜ θ 3 ˜ θ 4 ˜ 1 2 .
In the previous expression, θ ˜ 1 , θ ˜ 2 , and θ ˜ 5 are independent of n and t.
Only θ ˜ 3 , θ ˜ 4 , and θ ˜ 6 depend on n and t.
We choose the particular case in which t = 0 .
We replace n by x.
We denote c as the coefficient defined by
c = r n ( 0 ) θ 1 ˜ θ 2 ˜ 1 2 θ 5 ˜ .
We denote H as the function defined by
H ( z ) = F ( z ) G ( z ) F ( z ) + G ( z ) .
In the different sums involving θ , only terms with k j = 1 remain. So, the sums can be reduced only on the subsets J of [ 1 ; g ] of integer N . We denote a j as the term a j = 1 2 M j + i L j M j i L j . Then, the term θ ˜ 6 can be written as
θ ˜ 6 = J { 1 , , g } j , l J , j < l ln κ l j κ j l κ l j + κ j l 2 j J L j + M j L j M j
l J l j ln κ j l + κ l j κ j l κ l j 2 j J C j ( z l ) + D j ( z l ) C j ( z l ) + D j ( z l ) C j ( z ) D j ( z ) C j ( z ) + D j ( z ) e j J j π i e 2 x j J a j .
We denote τ as the function defined by
τ ( x ) = θ ˜ 3 ( x ) = J { 1 , , g } j , l J , j < l ln κ l j κ j l κ l j + κ j l 2 j J L j + M j L j M j
l J l j ln κ j l + κ l j κ j l κ l j 2 j J C j ( z l ) + D j ( z l ) C j ( z l ) D j ( z l ) e j J j π i e 2 x j J a j .
Thus, the solution Ψ ˜ can be written as
Ψ ˜ ( x , z ) = c exp ( x H ( z ) ) χ ( x , z ) τ ( x ) τ ( x 1 ) .
Here, we return to the Sato formulation of the solution, as given in [15].
This expression is very similar to the solution expressed on page 5830, and the expression of the solution is given by
Ψ ˜ ( x , z ) = exp ( x z ) χ ( x , z ) τ ( x ) τ ( x 1 ) .
The difference between these two statements comes from the fact that these two equations are treated differently:
Van Diejen considers the equation
a ( x ) f ( x + 1 ) + b ( x ) f ( x ) + a ( x 1 ) f ( x 1 ) = λ f ( x ) ;
here, we consider (1)
r n r n + 1 c n ψ n + 1 + v n ψ n + r n r n 1 c n 1 ψ n 1 = λ ψ n .

2.17. Limit of the Associated Potentials

In this subsection, we determine the limit of the associated potentials.
The potentials u n = c n and v n are described in (6) and (7) as
c n = e 2 I 0 θ 2 ( t V + n D + Z ) θ ( t V + ( n + 1 ) D + Z ) θ ( t V + ( n 1 ) D + Z )
v n = R 0 + t ln θ ( t V + ( n 1 ) D + Z ) θ ( t V + n D + Z ) .
So, the limit of u n is equal to
u n ˜ = k 1 τ ( x ) τ ( x + 1 ) τ ( x 1 ) ,
which is the same as in [15] (p. 5830), up to the constant k 1 .
Denoting σ ( x ) = t θ ( t V ˜ + x D ˜ + Z ˜ ) t = 0 , the limit of v n is equal to
v n ˜ = k 2 + σ ( x 1 ) τ ( x 1 ) σ ( x ) τ ( x ) ,
which is similar to the one written in [15] (p. 5830), up to the constant k 2 .

3. Conclusions

We have used the degeneracy of solutions of some PDEs given in terms of Riemann theta functions to obtain some important solutions. In particular, in the case of the NLS equation, we have managed to construct a quasi-rational solution involving the determinant of order N for each positive integer N depending on 2 N 2 real parameters [12]. In the case of the KdV equation, from abelian functions, we constructed solutions in terms of the Fredholm determinants and Wronskians [13].
In this study, we have managed to recover Sato formulas using this method for the Toda equation. From solutions given in terms of the Baker–Akiezer functions, we succeeded to construct by degeneracy, as given in the frame of the NLS equation [14], some quasi-rational solutions, and we recovered the Sato formulas for the Toda equation.
Precisely, we derived multisoliton solutions from finite gap solutions by a limit transition, i.e., by making gaps tend toward points in a certain Riemann surface. One strength of this approach is that it does not rely on inverse scattering theory or geometric and representation theoretic methods, which offers another perspective on the problem.
The degeneracy of the solutions to the NLS equation has allowed building new quasi-rational solutions of order N depending on 2 N 2 real parameters and their construction up to order 23, which were previously unknown. In the case of the KdV equation, the degeneracy of the solutions made it possible to find the solutions given by the Darboux method, which constituted a bridge between the geometric algebra approach and the Darboux transformations framework. In this article, this is somewhat the same situation as that of the KdV equation, where we linked the algebro–geometric approach to the framework of the inverse scattering method.

Funding

This research received no external funding.

Data Availability Statement

The data that support the findings of this study are available within the article.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Toda, M. Vibration of a chain with nonlinear interaction. J. Phys. Soc. Jap. 1967, 6, 22–431. [Google Scholar] [CrossRef]
  2. Date, E.; Tanaka, S. Periodic Multi-Soliton Solutions of Korteweg-de Vries Equation and Toda Lattice. Suppl. Prog. Theor. Phys. 1976, 59, 107–125. [Google Scholar] [CrossRef]
  3. Dubrovin, B.A. Theta functions and nonlinear equations. Russ. Math. Surv. 1981, 36, 11. [Google Scholar] [CrossRef]
  4. Stahlhofen, A.A.; Matveev, V.B. Positons for Toda lattice and related spectral problems. Phys. A Math. Gen. 1995, 28, 1957–1965. [Google Scholar] [CrossRef]
  5. Zhang, S.; Zhou, Y.Y. Multiwave solutions for the Toda lattice equation by generalizing exp function method. Int. J. Appl. Math. 2014, 44, 1–6. [Google Scholar]
  6. Duarte, V.N. A class of solutions of the two dimensional Toda lattice equation. Phys. Lett. 2021, 385, 126979. [Google Scholar] [CrossRef]
  7. Nilson, T.; Schiebold, C. Solution formulas for the two dimensional Toda lattice and particle like solitons with unexpected asymptotic behaviour. J. Nonliner. Math. Phys. 2020, 27, 57–94. [Google Scholar] [CrossRef]
  8. Bellokolos, E.D.; Bobenko, A.I.; Enol’skii, V.Z.; Its, A.R.; Matveev, V.B. Algebro-Geometric Approach to Nonlinear Integrable Equations; Series in Nonlinear Dynamics; Springer-Verlag: Berlin/Heidelberg, Germany, 1996. [Google Scholar]
  9. Dubrovin, B.A.; Matveev, V.B.; Novikov, S.P. Integrable Systems: Non-Linear Equations of Korteweg-De Vries Type, Finite-Zone Linear Operators, and Abelian Varieties; Lecture Notes Series; London Mathematical Society: London, UK, 1976; Volume 60, pp. 53–140. [Google Scholar]
  10. Krichever, I. The periodic non abelian Toda lattice and two dimensional generalization. Uspekhi Math Nauka 1981, 36, 72–77. [Google Scholar]
  11. Smirnov, A.O. Elliptic Solitons of the Integrable Nonlinear Equations; St.-Petersburg State University: Saint Petersburg, Russia, 2002; pp. 1–170. [Google Scholar]
  12. Gaillard, P. Families of quasi-rational solutions of the NLS equation and multi-rogue waves. J. Phys. A Meth. Theor. 2011, 44, 435204. [Google Scholar] [CrossRef]
  13. Gaillard, P. Degenerate Riemann theta functions, Fredholm and wronskian representations of the solutions to the KdV equation and the degenerate rational case. J. Geom. Phys. 2021, 161, 104059. [Google Scholar] [CrossRef]
  14. Its, A.R.; Rybin, A.V.; Salle, M.A. Exact integration of NLS equation. Theor. Math. Phys. 1988, 74, 29–32. [Google Scholar] [CrossRef]
  15. Diejen, J.F.V. A Sato formula for reflectionless finite difference operators. J. Math. Phys. 1999, 40, 5822–5834. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gaillard, P. From Algebro Geometric Solutions of the Toda Equation to Sato Formulas. AppliedMath 2024, 4, 856-867. https://doi.org/10.3390/appliedmath4030046

AMA Style

Gaillard P. From Algebro Geometric Solutions of the Toda Equation to Sato Formulas. AppliedMath. 2024; 4(3):856-867. https://doi.org/10.3390/appliedmath4030046

Chicago/Turabian Style

Gaillard, Pierre. 2024. "From Algebro Geometric Solutions of the Toda Equation to Sato Formulas" AppliedMath 4, no. 3: 856-867. https://doi.org/10.3390/appliedmath4030046

APA Style

Gaillard, P. (2024). From Algebro Geometric Solutions of the Toda Equation to Sato Formulas. AppliedMath, 4(3), 856-867. https://doi.org/10.3390/appliedmath4030046

Article Metrics

Back to TopTop