This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Linearization, separability and Lax pairs representation of a4(2)a_{4}^{\left(2\right)} Toda lattice

Bruce Lionnel LIETAP NDI University of Maroua
Faculty of Sciences, Department of Mathematics and Computer Sciences
P.O. Box 814 Maroua, Cameroon
[email protected]
Djagwa DEHAINSALA University of NDjamena
Faculty of Exact and Applied Sciences, Department of Mathematics
1 route de Farcha, PO Box 1027 NDjamena, Chad
[email protected]
 and  Joseph DONGHO University of Maroua
Faculty of Sciences, Department of Mathematics and Computer Sciences
P.O. Box 814 Maroua, Cameroon
[email protected]
Abstract.

The aim of this work is focused on linearizing and found the Lax Pairs of the algebraic complete integrability (a.c.i) Toda lattice associated with the twisted affine Lie algebra a4(2)a_{4}^{\left(2\right)}. Firstly, we recall that our case of a.c.i is a two-dimensional algebraic completely integrable systems for which the invariant (real) tori can be extended to complex algebraic tori (abelian surfaces). This implies that the geometry can be used to study this system. Secondly, we show that the lattice is related to the Mumford system and we construct an explicit morphism between these systems, leading to a new Poisson structure for the Mumford system. Finally, we give a new Lax equation for this Toda lattice and we construct an explicit linearization of the system.

Key words and phrases:
Toda lattice, integrable system, linearization, Lax representation.
1991 Mathematics Subject Classification:
34G20,34M55,37J35
This work was completed with the support of the Pr Joseph Dongho and Dr Djagwa Dehainsala.
Corresponding author: University of Maroua-Cameroon & University of NDjamena-Chad.

1. Introduction

Many integrable systems from classical mechanics admit a complexification, where phase space and time are complexified, and the geometry of the (complex) momentum map is the best possible complex analogue of the geometry that appears in the Liouville Theorem. Namely, in many relevant examples the generic complexified fiber is an affine part of an Abelian variety (a compact algebraic torus) and the integrable vector fields are translation invariant, when restricted to any of these tori. Such integrable systems are call them algebraic completely integrable systems, following the original definition of Adler and van Moerbeke.

Integrable systems have been integrated classically in terms of quadratures, usually through a sequence of very ingenious algebraic manipulations especially tailored to the problem. More recently, it was realized that whenever a system could be represented as a family of Lax pairs. the system could be linearized on the Jacobian of a spectral curve, defined by the characteristic polynomial of one of the matrices in the Lax pair.

To show that a Hamiltonian system linearizes on an Abelian variety, one may either construct a Lax representation of the differential equation depending on an extra-parameter and linearize on the Jacobian of the curve specified by its characteristic equation, or one may complete the complexified invariant manifolds by using the Laurent solutions of the differential equations. The latter method allows us in addition to identify the nature of the invariant manifolds and of the solutions of the system: in most examples the isospectral manifolds and the invariant manifolds are different.

In the previous work [3], we have prove that the a4(2)a_{4}^{\left(2\right)} is a two-dimensional integrable system. This system satisfies the linearization criterion [[1], theorem 6.41] and it is an algebraic completely integrable in the Adler-van Moerbeke sense. This system has a smooth hyperelliptic curve of genus two. According to Vanhaecke [7] and Mumford’s description of hyperelliptic Jacobians (see [[5], Section 3.1]), like Γ\Gamma is a hyperelliptic curve of genus two then the Riemann surface Γ¯\overline{\Gamma} is embedded in its jacobian such that Jac(Γ¯)ΓJac(\overline{\Gamma})\ \Gamma is isomorphic to the space of pairs of polynomials (u(λ);v(λ))(u(\lambda);v(\lambda)). u(λ)u(\lambda) is a monic of degree two and v(λ)v(\lambda) less than two. f(λ)v2(λ)f(\lambda)-v^{2}(\lambda) is divisible by u(λ)u(\lambda).

The aim of this paper is how we can linearize and find the Toda lattice a4(2)a_{4}^{(2)} Lax pair or Lax representation? To prove this, we construct an explicit map from the generic fiber 𝔽c\mathbb{F}_{c} into the Jacobian of the Riemann surface Γ¯c\overline{\Gamma}_{c}. After we find the kummer surface of Jac(𝒦c)Jac(\mathcal{K}_{c}), u(λ),v(λ)u(\lambda),v(\lambda) and f(λ)f(\lambda).

This paper is organized as follows. In section 2, preliminaries of this work, we give the basic notions of linearising, separating variables and Lax representation. In section 3, main part of the paper, we show that the a.c.i a4(2)a_{4}^{\left(2\right)} Toda lattice is related to the Mumford system and we construct an explicit morphism between these systems, leading to a new Poisson structure for the Mumford system. Finally, we give a new Lax equation with spectral parameter for this Toda lattice and we construct an explicit linearization of the system.

2. Preliminaries

Let n\mathbb{C}^{n} denote a complex vector space of dimension nn.

Definition 2.1.

[2] A lattice in n\mathbb{C}^{n} is a discrete subgroup of maximal rank in n\mathbb{C}^{n}. It is a free abelian group of rank 22.

A lattice Λ\Lambda in n\mathbb{C}^{n} acts in a natural way on the vector space n\mathbb{C}^{n} and the quotient 𝕋n=n/Λ\mathbb{T}^{n}=\mathbb{C}^{n}/\Lambda is called a complex torus.

In the theory of linear algebraic groups there is the notion of a torus. Such a torus is an affine group, whereas a complex torus is compact.

Definition 2.2.

[2] An abelian variety is a complex torus admitting a positive line bundle or equivalently a projective embedding.

Abelian varieties over the complex numbers are special complex tori, that is, quotients of finite-dimensional complex vector spaces modulo a lattice of maximal rank.

The Riemann Relations are necessary and sufficient conditions for a complex torus to be an abelian variety. They were introduced by Riemann in the special case of a Jacobian variety of a curve.

Let 𝕋n=n/Λ\mathbb{T}^{n}=\mathbb{C}^{n}/\Lambda be a complex torus.

Definition 2.3.

[2] A positive line bundle on 𝕋n\mathbb{T}^{n} is by definition a line bundle on 𝕋n\mathbb{T}^{n} whose first Chern class is a positive definite hermitian form on n\mathbb{C}^{n}.

A polarization on 𝕋n\mathbb{T}^{n} is by definition the first Chern class H=c1(L)H=c_{1}(L) of a positive line bundle LL on 𝕋n\mathbb{T}^{n}.

By abuse of notation we sometimes consider the line bundle LL on 𝕋n\mathbb{T}^{n} itself as a polarization. The type of LL is called the type of the polarization. A polarization is called principal if it is of type (1,,1)(1,\cdot\cdot\cdot,1).

Definition 2.4.

[2] An abelian variety is a complex torus 𝕋n\mathbb{T}^{n} admitting a polarization H=c1(L)H=c_{1}(L). The pair (𝕋n,H)(\mathbb{T}^{n},H) is called a polarized abelian variety.

According to [2], let Γ\Gamma be a smooth projective curve of genus gg over the field of complex numbers. the gg-dimensional \mathbb{C}-vector space H0(ωΓ)H^{0}(\omega_{\Gamma}) of holomorphic 11-forms on Γ\Gamma. The homology group H1(Γ,)H^{1}(\Gamma,\mathbb{Z}) is a free abelian group of rank 2g2g. For convenience we use the same letter for (topological) 11-cycles on Γ\Gamma and their corresponding classes in H1(Γ,)H^{1}(\Gamma,\mathbb{Z}). By Stoke’s theorem any element γH1(Γ,)\gamma\in H^{1}(\Gamma,\mathbb{Z}) yields in a canonical way a linear form on the vector space H0(ωΓ)H^{0}(\omega_{\Gamma}), which we also denote by:

γ:H0(ωΓ)ωγω\begin{array}[]{cccc}\gamma:&H^{0}(\omega_{\Gamma})&\longrightarrow&\mathbb{C}\\ &\omega&\longmapsto&\int_{\gamma}\omega\end{array}
Definition 2.5.

[2] the Jacobian variety or simply the Jacobian of Γ\Gamma, denote by Jac(Γ)Jac(\Gamma) is a complex torus of dimension gg such that

Jac(Γ):=H0(ωΓ)/H1(Γ,)Jac(\Gamma):=H^{0}(\omega_{\Gamma})^{\ast}/H^{1}(\Gamma,\mathbb{Z})
Definition 2.6.

[2] A theta divisor of the Jacobian Jac(Γ)Jac(\Gamma) is any divisor on Jac(Γ)Jac(\Gamma) such that the line bundle 𝒪Jac(Γ)(Θ)\mathcal{O}_{Jac(\Gamma)}(\Theta) defines the canonical polarization.

Definition 2.7.

[4] A system of ordinary differential equations over \mathbb{R} is called algebraic complete integrable (a.c.i.) when it is completely integrable and the complexified invariant manifolds complete into algebraic tori (Abelian varieties), whose (complexified) commuting flows extend holomorphically.

According to [4], Let 𝕋n=n/Λ\mathbb{T}^{n}=\mathbb{C}^{n}/\Lambda be a complex algebraic torus, (Abelian variety) with an origin 0 chosen. Let ii be the inverse morphism which coincides with the (1)(-1)-reflection about 0.

Definition 2.8.

[4] The Kumrner variety of 𝕋n\mathbb{T}^{n}, denoted by 𝒦c\mathcal{K}_{c}, is the quotient of 𝕋n\mathbb{T}^{n} by the action of the group (1,i)(1,i).

The Kummer variety bears the moduli information and has the advantage of possessing a lower degree of embedding in projective space. According to [4], let 𝒟\mathcal{D} be a divisor on 𝕋n\mathbb{T}^{n}. Denote by (𝒟)\mathcal{L}(\mathcal{D}) the invertible sheaf associated to 𝒟\mathcal{D}.

(𝒟)\displaystyle\mathcal{L}(\mathcal{D}) =\displaystyle= { the vector space of functions ff such that
(f)=divisor of zeroes-divisor of poles𝒟}\displaystyle(f)=\mbox{divisor of zeroes-divisor of poles}\geq-\mathcal{D}\}

According to [7] Let Γ\Gamma be a smooth curve of genus gg. We define two divisor 𝒟\mathcal{D} and 𝒟\mathcal{D}^{\prime} in Div(Γ)Div(\Gamma), the divisor group of Γ\Gamma, to be linearly equivalent, 𝒟l𝒟\mathcal{D}\sim_{l}\mathcal{D}^{\prime}, if and only if there exists a meromorphic function ff on Γ\Gamma.

According to [4], let 𝒟\mathcal{D} be an ample divisor on 𝕋n\mathbb{T}^{n}. We denote by 𝒞(𝒟)\mathcal{C(D)} the set of all divisors 𝒟\mathcal{D}^{\prime} on 𝕋n\mathbb{T}^{n} such that there are two positive numbers n,nn,n^{\prime} and n𝒟n\mathcal{D} is algebraically equivalent to n𝒟n^{\prime}\mathcal{D}^{\prime}.

Definition 2.9.

[1] A compact Riemann surface for which the Kodaira map is not an embedding is called a hyperelliptic Riemann surface (a compact Riemann surfaces of genus 11 being called an elliptic Riemann surface), while any curve whose (compact) Riemann surface is hyperelliptic is called a hyperelliptic curve (one speaks of an elliptic curve in the genus 11 case).

3. Separability and linearization of two-dimensional Toda lattice a4(2)a_{4}^{(2)}

3.1. Linearization procedure

According to [1], since Jac(Γ)Jac(\Gamma) is a principally polarized Abelian variety of dimension gg, the Lefschetz Theorem implies that it can be embedded in 3g1\mathbb{P}^{3^{g}-1} , by using the sections of [3Θ]\left[3\Theta\right]. However, the sections of [2Θ]\left[2\Theta\right] never embed Jac(Γ)Jac(\Gamma) in projective space, but rather they embed its Kummer variety Kc(Γ)K_{c}(\Gamma) in projective space. An important particular case is that of the Kummer surface Kc(Γ)K_{c}(\Gamma), where Γ\Gamma is a hyperelliptic Riemann surface of genus 22. The line bundle [2Θ]\left[2\Theta\right] that corresponds to twice the principal polarization on Jac(Γ)Jac(\Gamma) has in this case 44 independent sections and the associated Kodaira map, which maps Jac(Γ)Jac(\Gamma) into 3\mathbb{P}^{3} , factors through Kc(Γ)K_{c}(\Gamma), realizing the Kummer surface as a surface in 3\mathbb{P}^{3}.

Being two-dimensional the image is given by a single equation; to compute the degree of this equation, we use the fact that this degree is given by Kc(Γ)ω\displaystyle\int_{K_{c}(\Gamma)}\omega, where ω\omega is associated (1,1)(1,1)-form ofthe standard Kahler structure on 3\mathbb{P}^{3} . Clearly this is twice the volume of Kc(Γ)K_{c}(\Gamma), which itself is half the volume of the Jacobi surface (with the polarization of type (1,1)(1,1)).

In the two-dimensional case, the invariant manifolds complete into Abelian surfaces by adding one (or several) curves to the affine surfaces. In this case, Vanhaecke proposed in [7] a method which leads to an explicit linearization of the vector field of the a.c.i. system. The computation of the first few terms of the Laurent solutions to the differential equations enables us to construct an embedding of the invariant manifolds in the projective space N\mathbb{P}^{N}. From this embedding, one deduces the structure of the divisors 𝒟c\mathcal{D}_{c} to be adjoined to the generic affine in order to complete them into Abelian surfaces 𝕋c\mathbb{T}_{c}. Thus, the system is a.c.i.. The different steps of the algorithm of Vanhaecke are given by:
case 1

  • a)

    If one of the components of 𝒟c\mathcal{D}_{c} is a smooth curve Γc\Gamma_{c} of genus two, compute the image of the rational map φ[2Γc]:𝕋c23\varphi_{[2\Gamma_{c}]}:\mathbb{T}_{c}^{2}\rightarrow\mathbb{P}^{3} which is a singular surface in 3\mathbb{P}^{3}, the Kummer surface 𝒦c\mathcal{K}_{c} of jacobian Jac(Γc)Jac(\Gamma_{c}) of the curve Γc\Gamma_{c}.

  • b)

    Otherwise, if one of the components of 𝒟c\mathcal{D}_{c} is a d:1d:1 unramified cover 𝒞c\mathcal{C}_{c} of a smooth curve Γc\Gamma_{c} of genus two, the map p:𝒞cΓcp:\mathcal{C}_{c}\rightarrow\Gamma_{c} extends to the map p~:𝕋c2Jac(Γc)\widetilde{p}:\mathbb{T}_{c}^{2}\rightarrow Jac(\Gamma_{c}). In this case, let 𝒞c\mathcal{C}_{c} denote the (non complete) linear system p~[2Γc][2𝒞c]\widetilde{p}[2\Gamma_{c}]\subset[2\mathcal{C}_{c}] which corresponds to the complete linear system [2𝒞c][2\mathcal{C}_{c}] and compute now the Kummer surface 𝒞c\mathcal{C}_{c} of Jac(Γc)Jac(\Gamma_{c}) as image of φεc:𝕋c23\varphi_{\varepsilon_{c}}:\mathbb{T}_{c}^{2}\rightarrow\mathbb{P}^{3}.

  • c)

    Otherwise, change the divisor at infinity so as to arrive in case (a) or (b). This can always be done for any irreducible Abelian surface.

case 2. Choose a Weierstrass point WW on the curve Γc\Gamma_{c} and coordinates (z0:z1:z2:z3)(z_{0}:z_{1}:z_{2}:z_{3}) for 3\mathbb{P}^{3} such φ[2Γc](W)=(0:0:0:1)\varphi_{[2\Gamma_{c}]}(W)=(0:0:0:1) in case 1.(a) and φεc(W)=(0:0:0:1)\varphi_{\varepsilon_{c}}(W)=(0:0:0:1) in case 1.(b). Then this point will be a singular point (node) for the Kummer surface 𝒦c\mathcal{K}_{c} whose equation is p2(zo;z1;z2)z32+p3(zo;z1;z2)z3+p4(zo;z1;z2)=0\small{p_{2}(z_{o};z_{1};z_{2})z^{2}_{3}+p_{3}(z_{o};z_{1};z_{2})z_{3}+p_{4}(z_{o};z_{1};z_{2})=0}

where the pip_{i} are polynomials of degree ii. After a projective transformation which fixes (0:0:0:1)(0:0:0:1), we may assume that p2(zo;z1;z2)=z124z0z2p_{2}(z_{o};z_{1};z_{2})=z^{2}_{1}-4z_{0}z_{2}.

case 3. Finally, let s1s_{1} and s2s_{2} be the roots of the quadractic equation z0s2+z1s+z2=0z_{0}s^{2}+z_{1}s+z_{2}=0, whose discriminant is p2(zo;z1;z2)p^{2}(z_{o};z_{1};z_{2}), with the ziz_{i} expressed in terms of the original variables. Then the differential equations describing the vector field of the system are rewritten by direct computation in the classical Weierstrass form

(3.1) s1˙f(s1)+s2˙f(s2)=α1dts1s1˙f(s1)+s2s2˙f(s2)=α2dt\begin{array}[]{l}\frac{\dot{s_{1}}}{\sqrt{f(s_{1})}}+\frac{\dot{s_{2}}}{\sqrt{f(s_{2})}}=\alpha_{1}dt\\ \frac{s_{1}\dot{s_{1}}}{\sqrt{f(s_{1})}}+\frac{s_{2}\dot{s_{2}}}{\sqrt{f(s_{2})}}=\alpha_{2}dt\end{array}

where α1\alpha_{1} and α2\alpha_{2} depend on the torus. From it, the symmetric functions s1+s2:=z1z0s_{1}+s_{2}:=-\frac{z_{1}}{z_{0}} , s1s2:=z2z0s_{1}s_{2}:=\frac{z_{2}}{z_{0}} and the original variables can be written in terms of the Riemann theta function associated to the curve y2=f(x)y^{2}=f(x).

3.2. A.C.I of a4(2)a_{4}^{\left(2\right)} Toda lattice

In this section, we recall, according to [3], some results relating the two-dimensional a4(2)a_{4}^{\left(2\right)} Toda lattice. It is well known that this system is a.c.i.

The Toda lattice, introduced by Morikazu Toda in 19671967 [6], is a simple model for a one-dimensional crystal in solid-state physics. It is famous because it is one of the first examples of a completely integrable nonlinear system. It is described by a chain of particles with nearest-neighbor interaction, and its dynamics are governed by the Hamiltonian

H(p,q)=n(p2(n,t)2+V(q(n+1,t)q(n,t))),H\left(p,q\right)=\displaystyle\sum_{n\in\mathbb{Z}}\left(\frac{p^{2}\left(n,t\right)}{2}+V\left(q\left(n+1,t\right)-q\left(n,t\right)\right)\right),

and the equations of motion

{ddtp(n,t)=H(p,q)q(n,t)=e(q(n,t)q(n1,t))e(q(n+1,t)q(n,t))ddtq(n,t)=H(p,q)p(n,t)=p(n,t)\left\{\begin{array}[]{ll}\frac{d}{dt}p\left(n,t\right)=-\frac{\partial H\left(p,q\right)}{\partial q\left(n,t\right)}=e^{-\left(q\left(n,t\right)-q\left(n-1,t\right)\right)}-e^{-\left(q\left(n+1,t\right)-q\left(n,t\right)\right)}\\ \frac{d}{dt}q\left(n,t\right)=\frac{\partial H\left(p,q\right)}{\partial p\left(n,t\right)}=p\left(n,t\right)\end{array}\right.

where q(n,t)q\left(n,t\right) is the displacement of the nn-th particle from its equilibrium position, and p(n,t)p\left(n,t\right) is its momentum (with mass m=1m=1), and the Toda potential is given by V(r)=er+r1V\left(r\right)=e^{-r}+r-1. The classical Toda lattice is a system of particles with unit mass, connected by exponential springs. Its equations of motion derived from the Hamiltonian.

(3.2) H=12j=1npj2+j=1n1eqjqj+1.H=\frac{1}{2}\sum_{j=1}^{n}p_{j}^{2}+\sum_{j=1}^{n-1}e^{q_{j}-q_{j+1}}.

where qjq_{j} is the position of the j-th particle and pjp_{j} is its amount of movement. This type of Hamiltonian was considered first by Morikazu Toda [6]. The equation (3.2) is known as the finite classic no periodic Toda lattice to distinguish other versions of various forms of the system. The periodic version of (3.2) is given by

H=12j=1npj2+j=1neqjqj+1,qn+1=q1.H=\frac{1}{2}\sum_{j=1}^{n}p_{j}^{2}+\sum_{j=1}^{n}e^{q_{j}-q_{j+1}},q_{n+1}=q_{1}.

where the equations of motion are given by

p˙j=Hqj=e(qj1qj)e(qjqj+1) and q˙j=Hpj=pj,1jn.\dot{p}_{j}=-\frac{\partial H}{\partial q_{j}}=e^{\left(q_{j-1}-q_{j}\right)}-e^{\left(q_{j}-q_{j+1}\right)}\mbox{ and }\dot{q}_{j}=\frac{\partial H}{\partial p_{j}}=p_{j},1\leq j\leq n.

The differential equations of the periodic Toda lattice a4(2)a_{4}^{\left(2\right)} are given on the five dimensions hyperplane ={(x0,x1,x2,y0,y1,y2)6|y0+2y1+2y2=0} of 6\mathcal{H}=\{\left(x_{0},x_{1},x_{2},y_{0},y_{1},y_{2}\right)\in\mathbb{C}^{6}|y_{0}+2y_{1}+2y_{2}=0\}\mbox{ of }\mathbb{C}^{6} by

{x˙=x.yy˙=Ax\left\{\begin{array}[]{l}\dot{x}=x.y\\ \dot{y}=Ax\end{array}\right.

where x=(x0,x1,x2)x=\left(x_{0},x_{1},x_{2}\right)^{\top} , y=(y0,y1,y2)y=\left(y_{0},y_{1},y_{2}\right)^{\top} and AA is the Cartan matrix of the twisted affine Lie algebra a4(2)a_{4}^{\left(2\right)} given in [1] by

(220122012)\left(\begin{array}[]{ccc}2&-2&0\\ -1&2&-2\\ 0&-1&2\\ \end{array}\right)

and ε=(1,2,2)\varepsilon=\left(1,2,2\right)^{\top} is the normalized null vector of AA^{\top}. The equations of motion of the Toda lattice a4(2)a_{4}^{\left(2\right)} are given in [1] by :

(3.3) x0˙=x0y0y0˙=2x02x1x1˙=x1y1y1˙=x0+2x12x2x2˙=x2y2y2˙=x1+2x2\begin{array}[]{llllll}\dot{x_{0}}=x_{0}y_{0}&&&&&\dot{y_{0}}=2x_{0}-2x_{1}\\ \dot{x_{1}}=x_{1}y_{1}&&&&&\dot{y_{1}}=-x_{0}+2x_{1}-2x_{2}\\ \dot{x_{2}}=x_{2}y_{2}&&&&&\dot{y_{2}}=-x_{1}+2x_{2}\end{array}

We denote by 𝒱1\mathcal{V}_{1} the vector field defined by the above differential equations (3.3). Then 𝒱1\mathcal{V}_{1} is the Hamiltonian vector field, with Hamiltonian function F2=y02+4y224x08x116x2F_{2}=y_{0}^{2}+4y_{2}^{2}-4x_{0}-8x_{1}-16x_{2}
with respect to the Poisson structure {,}\{\cdot,\cdot\} defined by the following skew-symmetric matrix

(3.4) J=18(0004x02x000002x12x1x10000x2x24x02x100002x02x1x20000x1x2000)J=\frac{1}{8}\left(\begin{array}[]{cccccc}0&0&0&4x_{0}&-2x_{0}&0\\ 0&0&0&-2x_{1}&2x_{1}&-x_{1}\\ 0&0&0&0&-x_{2}&x_{2}\\ -4x_{0}&2x_{1}&0&0&0&0\\ 2x_{0}&-2x_{1}&x_{2}&0&0&0\\ 0&x_{1}&-x_{2}&0&0&0\end{array}\right)

This Poisson structure is given on 6\mathbb{C}^{6}; the function F0=y0+2y1+2y2F_{0}=y_{0}+2y_{1}+2y_{2} is a Casimir, so that the hyperplane \mathcal{H} is a Poisson subvariety. The rank of this Poisson structure {,}\{\cdot,\cdot\} is 0 on the three-dimensional subspace {x0=x1=x2=0}\{x_{0}=x_{1}=x_{2}=0\}; the rank is 22 on the three four-dimensional subspaces: {x0=x1=0}\{x_{0}=x_{1}=0\}, {x0=x2=0}\{x_{0}=x_{2}=0\} and {x1=x2=0}\{x_{1}=x_{2}=0\}. Thus, for all points of \mathcal{H} except the four subspaces above the rank is 44. The vector field 𝒱1\mathcal{V}_{1} admits also the following two constants of motion:

(3.5) F1=x0x12x22F2=y02+4y224x08x116x2F3=(y024x0)(y224x2)4x1(y0y24x2x1)\begin{array}[]{l}F_{1}=x_{0}x_{1}^{2}x_{2}^{2}\\ F_{2}=y_{0}^{2}+4y_{2}^{2}-4x_{0}-8x_{1}-16x_{2}\\ F_{3}=\left(y_{0}^{2}-4x_{0}\right)\left(y_{2}^{2}-4x_{2}\right)-4x_{1}\left(y_{0}y_{2}-4x_{2}-x_{1}\right)\end{array}

F1F_{1} is a Casimir for {,}\{\cdot,\cdot\}, and the function F3F_{3} generates a second Hamiltonian vector field 𝒱2\mathcal{V}_{2}, which commutes with 𝒱1\mathcal{V}_{1}, given by the differential equations

(3.6) x0=x0y2(y0y22x1)4x0x2y0x1=x1y1y2(y1+y2)x12y1+x1(x0y2+2x2y0)x2=x2(y1+y2)((y1+y2)y2+x1)+x0x2y0y0=2(2x1x2+x0y22)+x1(2x1y0y2)8x0x2y1=x0y22+2x2(3x0x1)+y0y2(x1+x2)2x12+x2y0y1y2=x1y2(y1+y2)+x12x2(y1+y2)2x2x0\begin{array}[]{l}x_{0}^{{}^{\prime}}=x_{0}y_{2}\left(y_{0}y_{2}-2x_{1}\right)-4x_{0}x_{2}y_{0}\\ x_{1}^{{}^{\prime}}=-x_{1}y_{1}y_{2}\left(y_{1}+y_{2}\right)-x_{1}^{2}y_{1}+x_{1}\left(x_{0}y_{2}+2x_{2}y_{0}\right)\\ x_{2}^{{}^{\prime}}=x_{2}\left(y_{1}+y_{2}\right)\left(\left(y_{1}+y_{2}\right)y_{2}+x_{1}\right)+x_{0}x_{2}y_{0}\\ y_{0}^{{}^{\prime}}=2\left(2x_{1}x_{2}+x_{0}y_{2}^{2}\right)+x_{1}\left(2x_{1}-y_{0}y_{2}\right)-8x_{0}x_{2}\\ y_{1}^{{}^{\prime}}=-x_{0}y_{2}^{2}+2x_{2}\left(3x_{0}-x_{1}\right)+y_{0}y_{2}\left(x_{1}+x_{2}\right)-2x_{1}^{2}+x_{2}y_{0}y_{1}\\ y_{2}^{{}^{\prime}}=x_{1}y_{2}\left(y_{1}+y_{2}\right)+x_{1}^{2}-x_{2}\left(y_{1}+y_{2}\right)-2x_{2}x_{0}\end{array}

Hence the system (3.3) is completely integrable in the Livouille sense. It can be written as a Hamiltonian vector fields

z˙=JHz,z=(z1,,z6)=(x0,x1,x2,y0,y1,y2)\dot{z}=J\frac{\partial H}{\partial z},z=\left(z_{1},\cdots,z_{6}\right)^{\top}=\left(x_{0},x_{1},x_{2},y_{0},y_{1},y_{2}\right)^{\top}

where H=F2.H=F_{2}. the Hamiltonian structure is defined by the following Poisson bracket

{F,H}=Fz,JHz=i,k=16JikFziHzk\{F,H\}=\left\langle\frac{\partial F}{\partial z},J\frac{\partial H}{\partial z}\right\rangle=\displaystyle\sum_{i,k=1}^{6}J_{ik}\frac{\partial F}{\partial z_{i}}\frac{\partial H}{\partial z_{k}}

where Hz=(Hx0,Hx1,Hx2,Hy0,Hy1,Hy2)\frac{\partial H}{\partial z}=\left(\frac{\partial H}{\partial x_{0}},\frac{\partial H}{\partial x_{1}},\frac{\partial H}{\partial x_{2}},\frac{\partial H}{\partial y_{0}},\frac{\partial H}{\partial y_{1}},\frac{\partial H}{\partial y_{2}}\right)^{\top} and JJ is an antisymmetric matrix.

The vector field 𝒱2\mathcal{V}_{2} admits the same constants of motion (3.5) and is in involution with 𝒱1\mathcal{V}_{1} therefore {F2,F3}=0\{F_{2},F_{3}\}=0. The involution σ\sigma defined on 6\mathbb{C}^{6} by

σ(x0,x1,x2,y0,y1,y2)=(x0,x1,x2,y0,y1,y2)\sigma\left(x_{0},x_{1},x_{2},y_{0},y_{1},y_{2}\right)=\left(x_{0},x_{1},x_{2},-y_{0},-y_{1},-y_{2}\right)

preserves the constants of motion F1,F2F_{1},F_{2} and F3F_{3}, hence leave the fibers of the momentum map FF invariant. This involution can be restricts to the hyperplane \mathcal{H}.

Lemma 3.1.

[3] The system of differential equation (3.3) of the vector field 𝒱1\mathcal{V}_{1} has three distincts families of homogeneous Laurent solutions with weights depending on four (dim1)\left(dim\mathcal{H}-1\right) free parameters.

The set of regular values of the momentum map 𝐅\mathbf{F} is the Zariski open subset Ω\Omega defined by

Ω={c=(c1,c2,c3)3c10 and 256(3200000c12+2000c32c2c1225c3c23c1+c35)+1728c25c132c34c22+c33c240}.\begin{array}[]{cl}\Omega=&\left\{c=\left(c_{1},c_{2},c_{3}\right)\in\mathbb{C}^{3}\mid c_{1}\neq 0\mbox{ and }\right.\\ &\left.256\left(3200000c_{1}^{2}+2000c_{3}^{2}c_{2}c_{1}-225c_{3}c_{2}^{3}c_{1}+c_{3}^{5}\right)+1728c_{2}^{5}c_{1}-32c_{3}^{4}c_{2}^{2}+c_{3}^{3}c_{2}^{4}\neq 0\right\}.\end{array}

At a generic point c=(c1,c2,c3)3c=\left(c_{1},c_{2},c_{3}\right)\in\mathbb{C}^{3}, the fiber on cΩc\in\Omega of 𝐅\mathbf{F} is therefore:

𝔽c:=𝔽1(c)=i=13{m:Fi(m)=ci}\mathbb{F}_{c}:=\mathbb{F}^{-1}\left(c\right)=\bigcap_{i=1}^{3}\{m\in\mathcal{H}:F_{i}(m)=c_{i}\}

Hence we have the following result which prove that Toda lattice a4(2)a_{4}^{\left(2\right)} is a completely integrable system in the Liouville sense.

Proposition 3.2.

[3] For cΩc\in\Omega, the fiber 𝔽c\mathbb{F}_{c} over cc of the momentum FF is a smooth affine variety of dimension 22 and the rank of the Poisson structure (3.4) is maximal and equal to 44 at each point of 𝔽c\mathbb{F}_{c} ; moreover the vector fields 𝒱1\mathcal{V}_{1} and 𝒱2\mathcal{V}_{2} are independent at each point of the fiber 𝔽c\mathbb{F}_{c}.

Proposition 3.3.

[3] (,{,},𝐅)\left(\mathcal{H},\{\cdot,\cdot\},\mathbf{F}\right) is a completely integrable system describing the Toda lattice a4(2)a_{4}^{\left(2\right)} where 𝐅=(F1,F2,F3)\mathbf{F}=\left(F_{1},F_{2},F_{3}\right) and {,}\{\cdot,\cdot\} are given respectively by (3.5) and (3.4) with commuting vector fields (3.3) and (3.6).

The algebraic complete integrability of the a4(2)a_{4}^{\left(2\right)} Toda lattice was established in [3] by the following theorem

Theorem 3.4.

[3] Let (,{,},𝐅)\left(\mathcal{H},\{\cdot,\cdot\},\mathbf{F}\right) be an integrable system describing the Toda lattice a4(2)a_{4}^{\left(2\right)} where 𝐅=(F1,F2,F3)\mathbf{F}=\left(F_{1},F_{2},F_{3}\right) and {,}\{\cdot,\cdot\} are given respectively by (3.5) and (3.4) with commuting vector fields (3.3).

  • i)

    (,{,},𝐅)\left(\mathcal{H},\{\cdot,\cdot\},\mathbf{F}\right) is a weight homogeneous algebraical completely integrable system.

  • ii)

    For cΩc\in\Omega, the fiber 𝔽c\mathbb{F}_{c} of its momentum map is completed in an abelian surface 𝕋c2\mathbb{T}^{2}_{c} (the Jacobian of the hyperelliptic curve (of genus two) Γc¯(2)\overline{\Gamma_{c}}^{\left(2\right)}) by the addition of a singular divisor 𝒟c\mathcal{D}_{c} composed of three irreducible components: 𝒟c(0)\mathcal{D}_{c}^{\left(0\right)} defined by:

    Γc(0):16d2a8(256d3+8d2c2)a6+(1536d2+96dc2+8c3+c22)d2a4((8(8c3+48dc2+c22+512d2)d+2c2c3)d2+64c1)a2+(8d(c2c3+16dc3+64d2c2+512d3+2dc22)+c32)d2=0\begin{array}[]{ll}\Gamma_{c}^{\left(0\right)}:&16d^{2}a^{8}-\left(256d^{3}+8d^{2}c_{2}\right)a^{6}+\left(1536d^{2}+96dc_{2}+8c_{3}+c_{2}^{2}\right)d^{2}a^{4}-\left(\left(8\left(8c_{3}+48dc_{2}+c_{2}^{2}\right.\right.\right.\\ &\left.\left.\left.+512d^{2}\right)d+2c_{2}c_{3}\right)d^{2}+64c_{1}\right)a^{2}+\left(8d\left(c_{2}c_{3}+16dc_{3}+64d^{2}c_{2}+512d^{3}+2dc_{2}^{2}\right)+c_{3}^{2}\right)d^{2}=0\end{array}

    and 𝒟c(1)\mathcal{D}_{c}^{\left(1\right)} defined by:

    Γc(1):256ad3((4a2c2)216c3)d2+64c1=0,\Gamma_{c}^{\left(1\right)}:256ad^{3}-\left(\left(4a^{2}-c_{2}\right)^{2}-16c_{3}\right)d^{2}-+64c_{1}=0,

    two singular curves of respective genus 33 and 44 and one smooth curve and 𝒟c(2)\mathcal{D}_{c}^{\left(2\right)} defined by

    Γc(2):e4a4(8c1+c2e2)a2e264e5+4e2c1c2+4c3e4+16c12=0.\Gamma_{c}^{\left(2\right)}:e^{4}a^{4}-\left(8c_{1}+c_{2}e^{2}\right)a^{2}e^{2}-64e^{5}+4e^{2}c_{1}c_{2}+4c_{3}e^{4}+16c_{1}^{2}=0.

    of genus 22 and isomorphic to Γc(2)¯\overline{\Gamma_{c}^{\left(2\right)}}. The curves intercept each other as indicated in figure:

    [Uncaptioned image]

    Figure: Curves completing the invariant surfaces 𝔽c\mathbb{F}_{c} of the Toda lattice a4(2)a_{4}^{\left(2\right)} in abelian surfaces where 𝒟i\mathcal{D}_{i} is the curve 𝒟c(i)\mathcal{D}_{c}^{\left(i\right)} .

4. Linearization and Lax pairs of the a4(2)a_{4}^{\left(2\right)} Toda lattice

The involution (1)(-1) on the abelian surface give a singular surface, his Kummer surface. Here we give an equation of Kummer surface lie with the Jacobi surface 𝒯c2=Jac(Γ¯c)\mathcal{T}_{c}^{2}=Jac(\overline{\Gamma}_{c}) where Γ¯c\overline{\Gamma}_{c} is a hyperelliptic Riemann surface of genus 22 define above. The surface 𝒯c2\mathcal{T}_{c}^{2} is an abelian principal polarisation and the section of the line bundle [2𝒟c(2)]\left[2\mathcal{D}_{c}^{(2)}\right] embed his Kummer surface in the projectif space 𝒫6\mathcal{P}^{6}

Consider the functions which have a double pole on one of component of divisor 𝒟c\mathcal{D}_{c}, namely 𝒟c(2)\mathcal{D}_{c}^{(2)} and no pole on the other.

Now, we find a basis function on \mathcal{H} which has a double pole in tt when we substitute the principal balance x(t;m2)x(t;m_{2}) and no poles when the other principal balances are substituted. Using x(t,m0)x(t,m_{0}), x(t,m1)x(t,m_{1}) and x(t,m2)x(t,m_{2}) give in [3], we obtain a basis of these functions constituate by the functions θi\theta_{i} give in the following table:

k dimk\dim\mathcal{F}^{k} dimk\dim\mathcal{H}^{k} dim𝒵ρk\dim\mathcal{Z}^{k}_{\rho} dep\sharp dep ζk\zeta^{k} indep. functions
0 1 1 1 0 1 θ0\theta_{0}
1 2 0 0 0 0 -
2 6 1 2 1 1 θ1\theta_{1}
3 10 0 0 0 0 -
4 20 2 4 3 1 θ2\theta_{2}
5 30 0 0 0 0 -
6 50 2 6 5 1 θ3\theta_{3}
7 70 0 0 0 0 -
8 105 3 8 8 0 -
(4.1) θ0=1θ1=x2θ2=x1x2+4x22y22x2θ3=x1x22\begin{array}[]{cl}\theta_{0}=&1\\ \theta_{1}=&x_{2}\\ \theta_{2}=&x_{1}x_{2}+4x_{2}^{2}-y_{2}^{2}x_{2}\\ \theta_{3}=&x_{1}x_{2}^{2}\end{array}

The four functions θi\theta_{i} are the line bundle section [2𝒟c(2)]\left[2\mathcal{D}_{c}^{(2)}\right].
Hence we can formulate the following result:

Proposition 4.1.

The Koidara map which correspond to these functions:

ψc:Jac(Γ¯c)𝒫3m=(x0,x1,x2,y0,y2)(θ0(m):θ1(m):θ2(m):θ3(m)),\begin{array}[]{cccc}\psi_{c}:&Jac(\overline{\Gamma}_{c})&\longrightarrow&\mathcal{P}^{3}\\ &m=(x_{0},x_{1},x_{2},y_{0},y_{2})&\longmapsto&(\theta_{0}(m):\theta_{1}(m):\theta_{2}(m):\theta_{3}(m)),\end{array}

applied the Jacobi surface 𝒯c2=Jac(Γ¯c)\mathcal{T}_{c}^{2}=Jac(\overline{\Gamma}_{c}) on his Kummer surface, which is a singular quartic in the projective space 𝒫3\mathcal{P}^{3}. The basis (θ0:θ1:θ2:θ3)(\theta_{0}:\theta_{1}:\theta_{2}:\theta_{3}) is taking convenably.

Proof.

By substitute the balance x(t,m2)x(t,m_{2}) in the θi\theta_{i}, i=0,,3i=0,...,3 functions and taking the coefficients of t2t^{-2} of Laurent series θi(t,m2)\theta_{i}(t,m_{2}), the map ψc\psi_{c} induce on Γc\Gamma_{c} a map

ψc(2):(a,e)(0:1:14e2(a2e2c2e24c1):e).\begin{array}[]{cccc}\psi_{c}^{(2)}:&(a,e)&\longmapsto&\left(0:1:\frac{1}{4e^{2}}\left(a^{2}e^{2}-c_{2}e^{2}-4c_{1}\right):e\right)\end{array}.

Consider a Weierstrass point on Γ¯c\overline{\Gamma}_{c} :a=ς1  , e=164(ς4c2ς2+4c3+O(ς6))\infty:a=\varsigma^{-1}\mbox{ }\mbox{ , }e=\frac{1}{64}\left(\varsigma^{-4}-c_{2}\varsigma^{-2}+4c_{3}+O\left(\varsigma^{6}\right)\right). we obtain

(4.4) ψc(2)()=limς0(0:64ς4:16ς216c2ς4+O(ς6):1c2ς2+O(ς4))=(0:0:0:1)\displaystyle\begin{array}[]{ccl}\psi_{c}^{(2)}(\infty)&=&\lim\limits_{\varsigma\rightarrow 0}\left(0:64\varsigma^{4}:16\varsigma^{2}-16c_{2}\varsigma^{4}+O(\varsigma^{6}):1-c_{2}\varsigma^{2}+O(\varsigma^{4})\right)\\ &=&\left(0:0:0:1\right)\end{array}

hence a basis (θ0:θ1:θ2:θ3)(\theta_{0}:\theta_{1}:\theta_{2}:\theta_{3}) is take convenably. ∎

Consider the constants of motion

(4.5) F1=x0x12x22=c1F2=y02+4y224x08x116x2=c2F3=(y024x0)(y224x2)4x1(y0y24x2x1)=c3\begin{array}[]{cl}F_{1}=&x_{0}x_{1}^{2}x_{2}^{2}=c_{1}\\ F_{2}=&y_{0}^{2}+4y_{2}^{2}-4x_{0}-8x_{1}-16x_{2}=c_{2}\\ F_{3}=&(y_{0}^{2}-4x_{0})(y_{2}^{2}-4x_{2})-4x_{1}(y_{0}y_{2}-4x_{2}-x_{1})=c_{3}\end{array}

and eliminating the variables (x0,x1,x2,y0,y2)(x_{0},x_{1},x_{2},y_{0},y_{2}) in the principals balances x(t,m0)x(t,m_{0}), x(t,m1)x(t,m_{1}) and x(t,m2)x(t,m_{2}) in [3] we obtain:

(4.6) x0=c1θ12θ32  ,  x1=θ3θ12  ,  x2=θ1x_{0}=\frac{c_{1}\theta_{1}^{2}}{\theta_{3}^{2}}\mbox{ }\mbox{ , }\mbox{ }x_{1}=\frac{\theta_{3}}{\theta_{1}^{2}}\mbox{ }\mbox{ , }\mbox{ }x_{2}=\theta_{1}

Using the second equations of (4.1) and (4.5), we obtain

(4.7) y02=1θ12θ32(4c1θ14+4θ33+θ1θ32(c2θ1+4θ2))    ,   y22=1θ12(4θ13θ1θ2+θ3)y_{0}^{2}=\frac{1}{\theta_{1}^{2}\theta_{3}^{2}}\left(4c_{1}\theta_{1}^{4}+4\theta_{3}^{3}+\theta_{1}\theta_{3}^{2}\left(c_{2}\theta_{1}+4\theta_{2}\right)\right)\mbox{ }\mbox{ }\mbox{ }\mbox{ , }\mbox{ }\mbox{ }y_{2}^{2}=\frac{1}{\theta_{1}^{2}}\left(4\theta_{1}^{3}-\theta_{1}\theta_{2}+\theta_{3}\right)

Rewriting the last equation of (4.5) on the follow form

4x1y0y2=((y024x0)(y224x2)c3)+4x1(4x2+x1)4x_{1}y_{0}y_{2}=\left((y_{0}^{2}-4x_{0})(y_{2}^{2}-4x_{2})-c_{3}\right)+4x_{1}\left(4x_{2}+x_{1}\right)

we obtain a Kummer surface of Jac(Γ¯c\overline{\Gamma}_{c}). It can be put in the follow form

(4.8) ((c2+16θ1)216(16θ2+4c2θ1+c3))θ32+2θ3f3(θ1,θ2)+f4(θ1,θ2)=0\left(\left(c_{2}+16\theta_{1}\right)^{2}-16\left(16\theta_{2}+4c_{2}\theta_{1}+c_{3}\right)\right)\theta_{3}^{2}+2\theta_{3}f_{3}\left(\theta_{1},\theta_{2}\right)+f_{4}\left(\theta_{1},\theta_{2}\right)=0

where f3f_{3} is a polynomial of degre 33, f4f_{4} of degre 44 in θ1\theta_{1} and θ2\theta_{2} given by

f3(θ1,θ2)\displaystyle f_{3}\left(\theta_{1},\theta_{2}\right) =\displaystyle= (c2+16θ1)(θ2(θ1c2+4θ2)+c3θ12)64c1\displaystyle-\left(c_{2}+16\theta_{1}\right)\left(\theta_{2}\left(\theta_{1}c_{2}+4\theta_{2}\right)+c_{3}\theta_{1}^{2}\right)-64c_{1}
f4(θ1,θ2)\displaystyle f_{4}\left(\theta_{1},\theta_{2}\right) =\displaystyle= (c3θ12+4θ22)2θ1(2θ12c2θ2c3+256c1θ12θ1c22θ2264c1θ28θ23c2)\displaystyle\left(c_{3}\theta_{1}^{2}+4\theta_{2}^{2}\right)^{2}-\theta_{1}\left(-2\theta_{1}^{2}c_{2}\theta_{2}c_{3}+256c_{1}\theta_{1}^{2}-\theta_{1}c_{2}^{2}\theta_{2}^{2}-64c_{1}\theta_{2}-8\theta_{2}^{3}c_{2}\right)

Hence we have the following results:

Proposition 4.2.

A quartic equation of the Kummer surface of Jac(Γ¯c\overline{\Gamma}_{c}), in terms of θi\theta_{i} is given by

((c2+16θ1)216(16θ2+4c2θ1+c3))θ32+2θ3f3(θ1,θ2)+f4(θ1,θ2)=0\left(\left(c_{2}+16\theta_{1}\right)^{2}-16\left(16\theta_{2}+4c_{2}\theta_{1}+c_{3}\right)\right)\theta_{3}^{2}+2\theta_{3}f_{3}\left(\theta_{1},\theta_{2}\right)+f_{4}\left(\theta_{1},\theta_{2}\right)=0

where f3f_{3} is a polynomial of degre 33, f4f_{4} of degre 44 in θ1\theta_{1} and θ2\theta_{2} given by

f3(θ1,θ2)\displaystyle f_{3}\left(\theta_{1},\theta_{2}\right) =\displaystyle= (c2+16θ1)(θ2(θ1c2+4θ2)+c3θ12)64c1\displaystyle-\left(c_{2}+16\theta_{1}\right)\left(\theta_{2}\left(\theta_{1}c_{2}+4\theta_{2}\right)+c_{3}\theta_{1}^{2}\right)-64c_{1}
f4(θ1,θ2)\displaystyle f_{4}\left(\theta_{1},\theta_{2}\right) =\displaystyle= (c3θ12+4θ22)2θ1(2θ12c2θ2c3+256c1θ12θ1c22θ2264c1θ28θ23c2)\displaystyle\left(c_{3}\theta_{1}^{2}+4\theta_{2}^{2}\right)^{2}-\theta_{1}\left(-2\theta_{1}^{2}c_{2}\theta_{2}c_{3}+256c_{1}\theta_{1}^{2}-\theta_{1}c_{2}^{2}\theta_{2}^{2}-64c_{1}\theta_{2}-8\theta_{2}^{3}c_{2}\right)
Theorem 4.3.

The vector field 𝒱1\mathcal{V}_{1} 3.3 extends to a linear vector field on the abelian surface 𝕋c2\mathbb{T}^{2}_{c} and the Jacobi form for the differentials equation can be written as

{λ1˙f(λ1)+λ2˙f(λ2)=0λ1λ1˙f(λ1)+λ2λ2˙f(λ2)=12idt\left\{\begin{array}[]{l}\frac{\dot{\lambda_{1}}}{\sqrt{f(\lambda_{1})}}+\frac{\dot{\lambda_{2}}}{\sqrt{f(\lambda_{2})}}=0\\ \frac{\lambda_{1}\dot{\lambda_{1}}}{\sqrt{f(\lambda_{1})}}+\frac{\lambda_{2}\dot{\lambda_{2}}}{\sqrt{f(\lambda_{2})}}=\frac{1}{2i}dt\end{array}\right.

with f(λ)=λi5+2c2λi4+(8c3+c22)λi3+8c2c3λi2+16c32λi16384c1f(\lambda)=\lambda_{i}^{5}+2c_{2}\lambda_{i}^{4}+\left(8c_{3}+c_{2}^{2}\right)\lambda_{i}^{3}+8c_{2}c_{3}\lambda_{i}^{2}+16c_{3}^{2}\lambda_{i}-16384c_{1} and v2=f(λ)v^{2}=f(\lambda) is birational equivalent to the hyperelliptic curve of genus two 𝒦c\mathcal{K}_{c}

Proof.

Consider coefficient of θ32\theta_{3}^{2} in equation (4.8) with the variables xix_{i} and yiy_{i}

Δ=(c2+16x2)24(4x2(16y22+64x2+16x1+4c2)+4c3)\Delta=\left(c_{2}+16x_{2}\right)^{2}-4\left(4x_{2}\left(-16y_{2}^{2}+64x_{2}+16x_{1}+4c_{2}\right)+4c_{3}\right)

Let u(λ)u(\lambda) an unitary polynomial in λ\lambda such that the discriminant is Δ\Delta, hence we have:

u(λ)\displaystyle u(\lambda) =\displaystyle= λ2+(c2+16x2)λ+4x2(16y22+64x2+16x1+4c2)+4c3\displaystyle\lambda^{2}+\left(c_{2}+16x_{2}\right)\lambda+4x_{2}\left(-16y_{2}^{2}+64x_{2}+16x_{1}+4c_{2}\right)+4c_{3}
=\displaystyle= λ2+(y02+4y224x08x1)λ+(4x12y0y2)216x0y22\displaystyle\lambda^{2}+\left(y_{0}^{2}+4y_{2}^{2}-4x_{0}-8x_{1}\right)\lambda+\left(4x_{1}-2y_{0}y_{2}\right)^{2}-16x_{0}y_{2}^{2}

.

Lets λ1\lambda_{1} and λ2\lambda_{2} roots of polynomial f(λ)f(\lambda), we have:

(4.9) λ1+λ2=16x2c2     ,    λ1λ2=4x2(16y22+64x2+16x1+4c2)+4c3\lambda_{1}+\lambda_{2}=-16x_{2}-c_{2}\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ , }\mbox{ }\mbox{ }\mbox{ }\lambda_{1}\lambda_{2}=4x_{2}\left(-16y_{2}^{2}+64x_{2}+16x_{1}+4c_{2}\right)+4c_{3}

that imply,with respect with 𝒱1\mathcal{V}_{1}

(4.10) λ˙1+λ˙2=16x2y2     ,    λ˙1λ2+λ1λ˙2=16x2(y2(y024x0)+2x1y0)\dot{\lambda}_{1}+\dot{\lambda}_{2}=-16x_{2}y_{2}\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ , }\mbox{ }\mbox{ }\mbox{ }\dot{\lambda}_{1}\lambda_{2}+\lambda_{1}\dot{\lambda}_{2}=-16x_{2}\left(-y_{2}\left(y_{0}^{2}-4x_{0}\right)+2x_{1}y_{0}\right)

Let v(λ)v\left(\lambda\right) a polynom define,up to a multiplicative constante, by :

(4.13) v(λ)=32i[x2y2λ+x2(y2(y024x0)2x1y0)]=2i(λ˙1+λ˙2)λ+2i(λ˙1λ2+λ1λ˙2).\displaystyle\begin{array}[]{ccc}v\left(\lambda\right)&=&32i\left[x_{2}y_{2}\lambda+x_{2}\left(y_{2}\left(y_{0}^{2}-4x_{0}\right)-2x_{1}y_{0}\right)\right]\\ &=&-2i\left(\dot{\lambda}_{1}+\dot{\lambda}_{2}\right)\lambda+2i\left(\dot{\lambda}_{1}\lambda_{2}+\lambda_{1}\dot{\lambda}_{2}\right).\end{array}

by substituation (4.9) and (4.10) in (4.5), and by eliminating variables x0,x1,x2,y0x_{0},x_{1},x_{2},y_{0} and y2y_{2}, we obtain two quadrics polynoms in λ˙i2\dot{\lambda}_{i}^{2} given by

λ˙i2=λi5+2c2λi4+(8c3+c22)λi3+8c2c3λi2+16c32λi16384c14(λ1λ2)2 , i=1,2\dot{\lambda}_{i}^{2}=\frac{\lambda_{i}^{5}+2c_{2}\lambda_{i}^{4}+\left(8c_{3}+c_{2}^{2}\right)\lambda_{i}^{3}+8c_{2}c_{3}\lambda_{i}^{2}+16c_{3}^{2}\lambda_{i}-16384c_{1}}{4\left(\lambda_{1}-\lambda_{2}\right)^{2}}\mbox{ , }i=1,2

verify

(4.14) {λ1˙f(λ1)+λ2˙f(λ2)=0λ1λ1˙f(λ1)+λ2λ2˙f(λ2)=12idt\left\{\begin{array}[]{l}\frac{\dot{\lambda_{1}}}{\sqrt{f(\lambda_{1})}}+\frac{\dot{\lambda_{2}}}{\sqrt{f(\lambda_{2})}}=0\\ \frac{\lambda_{1}\dot{\lambda_{1}}}{\sqrt{f(\lambda_{1})}}+\frac{\lambda_{2}\dot{\lambda_{2}}}{\sqrt{f(\lambda_{2})}}=\frac{1}{2i}dt\end{array}\right.

with

f(λ)=λi5+2c2λi4+(8c3+c22)λi3+8c2c3λi2+16c32λi16384c1f(\lambda)=\lambda_{i}^{5}+2c_{2}\lambda_{i}^{4}+\left(8c_{3}+c_{2}^{2}\right)\lambda_{i}^{3}+8c_{2}c_{3}\lambda_{i}^{2}+16c_{3}^{2}\lambda_{i}-16384c_{1}

and like v2=f(λ)v^{2}=f(\lambda) then:

f(λl)\displaystyle\sqrt{f\left(\lambda_{l}\right)} =\displaystyle= v(λl)\displaystyle v\left(\lambda_{l}\right)
=\displaystyle= 2i[16x2y2λl+16ix2(y2(y024x0)2x1y0)]\displaystyle 2i\left[16x_{2}y_{2}\lambda_{l}+16ix_{2}\left(y_{2}\left(y_{0}^{2}-4x_{0}\right)-2x_{1}y_{0}\right)\right]
=\displaystyle= 2i(λ˙1+λ˙2)λl+2i(λ˙1λ2+λ1λ˙2)\displaystyle-2i\left(\dot{\lambda}_{1}+\dot{\lambda}_{2}\right)\lambda_{l}+2i\left(\dot{\lambda}_{1}\lambda_{2}+\lambda_{1}\dot{\lambda}_{2}\right)

hence

{f(λ1)=2i(λ˙1+λ˙2)λ1+2i(λ˙1λ2+λ1λ˙2)f(λ2)=2i(λ˙1+λ˙2)λ2+2i(λ˙1λ2+λ1λ˙2){f(λ1)=2i(λ1λ2)λ˙1f(λ2)=2i(λ1λ2)λ˙2\left\{\begin{array}[]{l}\sqrt{f\left(\lambda_{1}\right)}=-2i\left(\dot{\lambda}_{1}+\dot{\lambda}_{2}\right)\lambda_{1}+2i\left(\dot{\lambda}_{1}\lambda_{2}+\lambda_{1}\dot{\lambda}_{2}\right)\\ \sqrt{f\left(\lambda_{2}\right)}=-2i\left(\dot{\lambda}_{1}+\dot{\lambda}_{2}\right)\lambda_{2}+2i\left(\dot{\lambda}_{1}\lambda_{2}+\lambda_{1}\dot{\lambda}_{2}\right)\end{array}\right.\Longrightarrow\left\{\begin{array}[]{l}\sqrt{f\left(\lambda_{1}\right)}=-2i\left(\lambda_{1}-\lambda_{2}\right)\dot{\lambda}_{1}\\ \sqrt{f\left(\lambda_{2}\right)}=2i\left(\lambda_{1}-\lambda_{2}\right)\dot{\lambda}_{2}\end{array}\right.

This show that the Toda is linearising on the Jacobian variety of the curve 𝒦c\mathcal{K}_{c}. It is able to see how 𝒦c\mathcal{K}_{c} and v2=f(s)v^{2}=f(s) are related.

Like

𝒦c:z2=h(t)=t52c2t4+(8c3+c22)t38c2c3t2+16c32t+16384c1\mathcal{K}_{c}:z^{2}=h(t)=t^{5}-2c_{2}t^{4}+\left(8c_{3}+c_{2}^{2}\right)t^{3}-8c_{2}c_{3}t^{2}+16c_{3}^{2}t+16384c_{1}

and

v2=f(λ)=λ5+2c2λ4+(8c3+c22)λ3+8c2c3λ2+16c32λ16384c1v^{2}=f(\lambda)=\lambda^{5}+2c_{2}\lambda^{4}+\left(8c_{3}+c_{2}^{2}\right)\lambda^{3}+8c_{2}c_{3}\lambda^{2}+16c_{3}^{2}\lambda-16384c_{1}

then we easy verify by taking λ=t\lambda=-t that z=ivz=iv .

one verifies, by a direct computation, that the expression f(λ)v2(λ)f(\lambda)-v^{2}(\lambda) is divisible by u(λ)u(\lambda) with

f(λ)=λ5+2c2λ4+(8c3+c22)λ3+8c2c3λ2+16c32λ16384c1f(\lambda)=\lambda^{5}+2c_{2}\lambda^{4}+\left(8c_{3}+c_{2}^{2}\right)\lambda^{3}+8c_{2}c_{3}\lambda^{2}+16c_{3}^{2}\lambda-16384c_{1}

Hence y2=f(λ)y^{2}=f(\lambda) is birational to the affine curve Γc\Gamma_{c} by adding the Weierstrass points at infinity a=±t52c2t4+(8c3+c22)t38c2c3t2+16c32t+16384c1(t2c2t+4c3)2,   e=164(t2c2t+4c3)a=\pm\sqrt{\frac{t^{5}-2c_{2}t^{4}+\left(8c_{3}+c_{2}^{2}\right)t^{3}-8c_{2}c_{3}t^{2}+16c_{3}^{2}t+16384c_{1}}{\left(t^{2}-c_{2}t+4c_{3}\right)^{2}}},\mbox{ }\mbox{ }\mbox{ }e=\frac{1}{64}\left(t^{2}-c_{2}t+4c_{3}\right). ∎

The form 4.14 is ewuivalent to

ddt(k=120kQkω)=(02i)\frac{d}{dt}\left(\sum^{2}_{k=1}\int_{0_{k}}^{Q_{k}}\overrightarrow{\omega}\right)=\left(\begin{array}[]{c}0\\ 2i\\ \end{array}\right)

where ω=(dxf(x),xdxf(x))\overrightarrow{\omega}=\left(\frac{dx}{\sqrt{f(x)}},\frac{xdx}{\sqrt{f(x)}}\right)^{\top} is a basis for holomorphic differentials on Γ¯c\overline{\Gamma}_{c}, Q1:=(λ1,f(λ1))Q_{1}:=\left(\lambda_{1},\sqrt{f(\lambda_{1})}\right) and Q2:=(λ2,f(λ2))Q_{2}:=\left(\lambda_{2},\sqrt{f(\lambda_{2})}\right) two points of Γc\Gamma_{c} and Q1+Q2=(λ1,f(λ1))+(λ2,f(λ2))Q_{1}+Q_{2}=\left(\lambda_{1},\sqrt{f(\lambda_{1})}\right)+\left(\lambda_{2},\sqrt{f(\lambda_{2})}\right) viewed as a divisor on the genus 22 hyperelliptic curve Γc\Gamma_{c}. Thus, by integrating 4.14, we see that the flow of 𝒱¯1\overline{\mathcal{V}}_{1} is linear on the Jacobian of the curve Γc\Gamma_{c}. By using [[5], Theorem 5.3], one shows that the symmetric functions λ1\lambda_{1} and λ2\lambda_{2}, and hence the original phase variables can be written in terms of theta functions.

Now we also establish a link between the a4(2)a_{4}^{\left(2\right)} Toda lattice and the Mumford system [5]. By using a method due to Vanhaecke [8], we construct an explicit morphism between these two systems. Thus, we obtain a new Poisson structure for the Mumford system and then derive a new Lax equation for the a4(2)a_{4}^{\left(2\right)} Toda lattice.

According the fact that the expression f(λ)v2(λ)f(\lambda)-v^{2}(\lambda) is divisible by u(λ)u(\lambda) such that the above formulas define a point of Jac(Γ¯c2)\Γc2Jac(\overline{\Gamma}^{2}_{c})\backslash\Gamma^{2}_{c}, there exist a polynomial ww in λ\lambda of degree 3=3=deg u+1u+1. By direct calculation, we obtain:

w(λ)\displaystyle w(\lambda) =\displaystyle= f(λ)v2(λ)u(λ)\displaystyle\frac{f(\lambda)-v^{2}(\lambda)}{u(\lambda)}
=\displaystyle= λ3+w2λ2+w1λ+w0,\displaystyle\lambda^{3}+w_{2}\lambda^{2}+w_{1}\lambda+w_{0},

where

w0\displaystyle w_{0} =\displaystyle= 256y02x221024x0x22\displaystyle 256y_{0}^{2}x_{2}^{2}-1024x_{0}x_{2}^{2}
w1\displaystyle w_{1} =\displaystyle= 16x12+4y02y2232y02x216x0y22+128x0x2+256x22+128x2x116x1y0y2\displaystyle 16x_{1}^{2}+4y_{0}^{2}y_{2}^{2}-32y_{0}^{2}x_{2}-16x_{0}y_{2}^{2}+128x_{0}x_{2}+256x_{2}^{2}+128x_{2}x_{1}-16x_{1}y_{0}y_{2}
w2\displaystyle w_{2} =\displaystyle= y028x132x24x0+4y22.\displaystyle y_{0}^{2}-8x_{1}-32x_{2}-4x_{0}+4y_{2}^{2}.

The linearizing variables 4.9 and 4.10 suggest a morphism φ\varphi from the a4(2)a^{(2)}_{4} Toda lattice to genus 22 odd Mumford system:

{(v(λ)u(λ)w(λ)v(λ))M2([λ]) such that deg(u)=2=deg(w)1deg(v)<2;u,w are monic }7,\left\{\left(\begin{array}[]{cc}v(\lambda)&u(\lambda)\\ w(\lambda)&-v(\lambda)\\ \end{array}\right)\in M_{2}(\mathbb{C}[\lambda])\mbox{ such that }\begin{array}[]{c}deg(u)=2=deg(w)-1\\ deg(v)<2;\mbox{u,w are monic }\end{array}\right\}\cong\mathbb{C}^{7},

where 7\mathbb{C}^{7} is a phase space of Mumford system. The morphism φ\varphi is given by:

(4.15) (x0,x1,x2,y0,y2){u(λ)=λ2+u1λ+u0v(λ)=v1λ+v0w(λ)=λ3+w2λ2+w1λ+w0\begin{array}[]{ccc}(x_{0},x_{1},x_{2},y_{0},y_{2})&\longmapsto&\left\{\begin{array}[]{ll}u(\lambda)=&\lambda^{2}+u_{1}\lambda+u_{0}\\ v(\lambda)=&v_{1}\lambda+v_{0}\\ w(\lambda)=&\lambda^{3}+w_{2}\lambda^{2}+w_{1}\lambda+w_{0}\end{array}\right.\end{array}

with

u0=(4x12y0y2)216x0y22v0=16x2(y2(y024x0)2x1y0)u1=(y02+4y224x08x1)v1=16x2y2\begin{array}[]{llllll}u_{0}=&\left(4x_{1}-2y_{0}y_{2}\right)^{2}-16x_{0}y_{2}^{2}&&&v_{0}=&16x_{2}\left(y_{2}\left(y_{0}^{2}-4x_{0}\right)-2x_{1}y_{0}\right)\\ u_{1}=&-\left(y_{0}^{2}+4y_{2}^{2}-4x_{0}-8x_{1}\right)&&&v_{1}=&16x_{2}y_{2}\end{array}
w0=256y02x221024x0x22w1=16x12+4y02y2232y02x216x0y22+128x0x2+256x22+128x2x116x1y0y2w2=y028x132x24x0+4y22\begin{array}[]{ll}w_{0}=&256y_{0}^{2}x_{2}^{2}-1024x_{0}x_{2}^{2}\\ w_{1}=&16x_{1}^{2}+4y_{0}^{2}y_{2}^{2}-32y_{0}^{2}x_{2}-16x_{0}y_{2}^{2}+128x_{0}x_{2}+256x_{2}^{2}+128x_{2}x_{1}-16x_{1}y_{0}y_{2}\\ w_{2}=&y_{0}^{2}-8x_{1}-32x_{2}-4x_{0}+4y_{2}^{2}\end{array}
Theorem 4.4.

A Lax representation of the vector field 𝒱1=𝒳F1\mathcal{V}_{1}=\mathcal{X}_{F_{1}} is given by:

X˙=[X(λ),Y(λ)]\dot{X}=\left[X(\lambda),Y(\lambda)\right]

by taking

X(λ)=(v(λ)u(λ)w(λ)v(λ)) and Y(λ)=(01b(λ)0)X(\lambda)=\left(\begin{array}[]{cc}v(\lambda)&u(\lambda)\\ w(\lambda)&-v(\lambda)\\ \end{array}\right)\mbox{ and }Y(\lambda)=\left(\begin{array}[]{cc}0&1\\ b(\lambda)&0\\ \end{array}\right)

where u(λ),v(λ)u(\lambda),v(\lambda) and w(λ)w(\lambda) are the polynomials defined above. The coefficient b(λ)=λ32x2b(\lambda)=\lambda-32x_{2} of the matrix Y(λ)Y(\lambda) is the polynomial part of the rational function w(λ)/u(λ)w(\lambda)/u(\lambda).

Acknowledgements

We would like to extend our sincere gratitude to Professor Pol Vanhaecke at University of Poitiers for his particular contributions in providing clarifications and guidance on our research theme, for the enriching exchanges and thoughtful advice he generously offered us throughout this project.

References

  • [1] Adler, M., Moerbeke,P.V and Vanhaecke, P., Intégrabilité algébrique, géométrie Painlevé et algèbres de Lie, Ergeb. Math. Grenzgeb. 47 (3) (2004) Berlin- Heidelberg : Springer.
  • [2] Herbert,L. and Birkenhake, C.,Abelian Varieties over the Complex Numbers,A Graduate Course, Grundlehren Text Editions, Springer,https://doi.org/10.1007/978-3-031-25570-0.
  • [3] Lietap, N.,Dehainsala, D. and Dongho, J., Algebraic complete integrability of the a4(2)a_{4}^{(2)} Toda lattice, University of maroua,SIGMA 20 (2024), 087, 26 pages, https://doi.org/10.3842/SIGMA.2024.087.
  • [4] Piovan, L., Algebraically completely integrable systems and Kummer varieties, Brandeis University, Math. Ann. 290, pp.349-403 (1991)
  • [5] Mumford, D.: Tata Lectures on Theta: 2. Jacobian Theta Functions and Differential Equations, Progr. Math., vol. 43, Boston, MA: Birkhäuser, 1984.
  • [6] Toda, M., One-Dimensional Dual Transformation, J. Phys. Soc. Japan, 1965, vol. 20, no. 11, pp. 2095A (see also Progr. Theoret. Phys. Suppl., 1966, no. 36, pp. 113-119).
  • [7] Vanhaecke, P., Linearising two-dimensional integrable systems and the construction of action-angle variables. Math. Z., 1992, vol. 211, no. 2, pp. 265-313.
  • [8] Vanhaecke, P.,Integrable Systems in the Realm of Algebraic Geometry, 2nd ed., Lecture Notes in Math., vol.1638, Berlin: Springer, 2001.