Abstract

We prove that an indecomposable principally polarized abelian variety
$X$
is the Jacobain of a curve if and only if there exist vectors
$U\ne 0,V$
such that the roots
${x}_{i}\left(y\right)$
of the theta-functional equation
$\theta (Ux+Vy+Z)=0$
satisfy the equations of motion of the formal infinite-dimensional Calogero-Moser system. 1 Introduction

The Riemann-Schottky problem on the characterization of the Jacobians of curves among abelian varieties is more that 120 years old. Quite a few geometrical characterizations of the Jacobians have been found. None of them provides an explicit system of equations for the image of the Jacobian locus in the projective space under the level two theta imbedding.

The first effective solution of the Riemann-Schottky problem was obtained by T.Shiota ([1] ), who proved the famous Novikov's conjecture:

An indecomposable principally polarized abelian variety
$(X,\theta )$
is the Jacobian of a curve of a genus g if and only if there exist
$g$
-dimensional vectors
$U\ne 0,V,W$
such that the function

is a solution of the Kadomtsev-Petviashvilii (KP) equation

Here
$\theta \left(Z\right)=\theta \left(Z\right|B)$
is the Riemann theta-function,

where
$B$
is the corresponding symmetric matrix with a positive definite imaginary part.

$$\begin{array}{c}u(x,y,t)=-2{\partial}_{x}^{2}ln\theta (Ux+Vy+Wt+Z)\end{array}$$ | (1.1) |

$$\begin{array}{c}3{u}_{yy}={\left(4{u}_{t}+6u{u}_{x}-{u}_{xxx}\right)}_{x}.\end{array}$$ | (1.2) |

$$\begin{array}{c}\theta \left(z\right)={\sum}_{m\in {\mathbb{Z}}^{g}}{e}^{2\pi i(z,m)+\pi i(Bm,m)},(z,m)={m}_{1}{z}_{1}+\cdots +{m}_{g}{z}_{g},\end{array}$$ | (1.3) |

It is easy to show ([2] ) that the KP equation with
$u$
of the form ( 1.1 ) is in fact equivalent to the following system of algebraic equations for the fourth order derivatives of the level two theta constants:

Here
$\Theta [\varepsilon ,0]=\Theta [\varepsilon ,0]\left(0\right)$
, where
$\Theta [\varepsilon ,0]\left(z\right)=\theta [\varepsilon ,0]\left(2z\right|2B)$
are level two theta-functions with half-integer characteristics
$\varepsilon \in \frac{1}{2}{\mathbb{Z}}_{2}^{g}$
.

$$\begin{array}{c}{\partial}_{U}^{4}\Theta [\varepsilon ,0]-{\partial}_{U}{\partial}_{W}\Theta [\varepsilon ,0]+{\partial}_{V}^{2}\Theta [\varepsilon ,0]+c\Theta [\varepsilon ,0]=0,c=const.\end{array}$$ | (1.4) |

The KP equation admits the so-called zero-curvature representation ([3, 4] ) which is the compatibility condition for the following over-determined system of linear equations:

The main goal of the present paper is to show that the KP equation contains an excessive information and that the Jacobians can be characterized in terms of only the first of its auxiliary linear equations.

$$\begin{array}{ccc}& \left({\partial}_{y}-{\partial}_{x}^{2}+u\right)\psi =0,& \end{array}$$ | (1.5) |

$$\begin{array}{ccc}& \left({\partial}_{t}-{\partial}_{x}^{3}+\frac{3}{2}{\partial}_{x}+w\right)\psi =0.& \end{array}$$ | (1.6) |

Theorem 1.1
An indecomposable principally polarized abelian variety
$(X,\theta )$
is the Jacobian of a curve of genus g if and only if there exist
$g$
-dimensional vectors
$U\ne 0,V,A$
such that equation ( 1.5 ) is satisfied for

and

where
$p,E$
are constants.

The “if ” part of this statement follows from the exact theta-functional expression for the Baker-Akhiezer function ([5, 6] ). $$\begin{array}{c}u=-2{\partial}_{x}^{2}ln\theta (Ux+Vy+Z)\end{array}$$ | (1.7) |

$$\begin{array}{c}\psi =\frac{\theta (A+Ux+Vy+Z)}{\theta (Ux+Vy+Z)}{e}^{px+Ey},\end{array}$$ | (1.8) |

The addition formula for the Riemann theta-function directly implies that equation ( 1.5 ) with
$u$
and
$\psi $
of the form ( 1.7 ) and ( 1.8 ) is equivalent to the system of equations.

Recently Theorem 1.1 was proved by E.Arbarello, G. Marini and the author ([7] ) under the additional assumption that the closure
$\langle A\rangle $
of the subgroup of
$X$
generated by
$A$
is irreducible.

$$\begin{array}{c}\left({\partial}_{V}-{\partial}_{U}^{2}-2p{\partial}_{U}+(E-{p}^{2})\right)\Theta [\varepsilon ,0](A/2)=0,\varepsilon \in \frac{1}{2}{\mathbb{Z}}_{2}^{g}.\end{array}$$ | (1.9) |

The geometric interpretation of Theorem 1.1 is equivalent to the characterization of the Jacobians via flexes of Kummer varieties (see details in [7] ), which is a particular case of the so-called trisecant conjecture, first formulated in [8] .

Theorem 1.1 is not the strongest form of our main result. What we really prove is that the Jacobian locus in the space of principally polarized abelian varieties is characterized by a system of equations which formally can be seen as the equations of motion of the infinite-dimensional Calogero-Moser system.

Let
$\tau (x,y)$
be an entire function of the complex variable
$x$
smoothly depending on a parameter
$y$
. Consider the equation

which means that the meromorphic function given by the left hand side of ( 1.10 ) has no residues in the
$x$
variable. If
${x}_{i}\left(y\right)$
is a simple zero of
$\tau $
, i.e.
$\tau \left({x}_{i}\right(y),y)=0,{\partial}_{x}\tau \left({x}_{i}\right(y),y)\ne 0$
, then ( 1.10 ) implies

where “dots” stands for the
$y$
-derivatives and
${w}_{i}$
is the third coefficient of the Laurent expansion of
$u(x,y)=-2{\partial}_{x}^{2}\tau (x,y)$
at
${x}_{i}$
, i.e.

Formally, if we represent
$\tau $
as an infinite product,

then equation ( 1.10 ) can be written as the infinite system of equations

Equations ( 1.14 ) are purely formal because, even if
$\tau $
has simple zeros at
$y=0$
, then in the general case there is no nontrivial interval in
$y$
where the zeros stay simple. At the moment, the only reason for representing ( 1.11 ) in the form ( 1.14 ) is to show that in the case when
$\tau $
is a rational, trigonometric or elliptic polynomial the system ( 1.11 ) coincides with the equations of motion for the rational, trigonometrical or elliptic Calogero-Moser systems, respectively.

$$\begin{array}{c}re{s}_{x}\left({\partial}_{y}^{2}ln\tau +2{\left({\partial}_{x}^{2}ln\tau \right)}^{2}\right)=0,\end{array}$$ | (1.10) |

$$\begin{array}{c}{\ddot{x}}_{i}=2{w}_{i},\end{array}$$ | (1.11) |

$$\begin{array}{c}u(x,y)=\frac{2}{(x-{x}_{i}(y){)}^{2}}+{v}_{i}\left(y\right)+{w}_{i}\left(y\right)(x-{x}_{i}(y\left)\right)+\cdots \end{array}$$ | (1.12) |

$$\begin{array}{c}\tau (x,y)=c\left(y\right){\prod}_{i}(x-{x}_{i}(y\left)\right),\end{array}$$ | (1.13) |

$$\begin{array}{c}{\ddot{x}}_{i}=-4{\sum}_{j\ne i}\frac{1}{({x}_{i}-{x}_{j}{)}^{3}}.\end{array}$$ | (1.14) |

Equations ( 1.11 ) for the zeros of the function
$\tau =\theta (Ux+Vy+Z)$
were derived in [7] as a direct corollary of the assumptions of Theorem 1.1. Simple expansion of
$\theta $
at the points of its divisor
$z\in \Theta :\theta \left(z\right)=0$
gives the equation

which is valid on
$\Theta $
. Here and below
$\Theta $
is the divisor on
$X$
defined by the equation
$\theta \left(Z\right)=0$
and
${\partial}_{1}$
and
${\partial}_{2}$
are constant vector field on
${\mathbb{C}}^{g}$
corresponding to the vectors
$U$
and
$V$
.

$$\begin{array}{c}\left[\right({\partial}_{2}\theta {)}^{2}-\left({\partial}_{1}^{2}\theta {)}^{2}\right]{\partial}_{1}^{2}\theta +2[{\partial}_{1}^{2}\theta {\partial}_{1}^{3}\theta -{\partial}_{2}\theta {\partial}_{1}{\partial}_{2}\theta ]{\partial}_{1}\theta +[{\partial}_{2}^{2}\theta -{\partial}_{1}^{4}\theta ]({\partial}_{1}\theta {)}^{2}=0(mod\theta )\end{array}$$ | (1.15) |

It would be very interesting to understand if any reasonable general theory of equations ( 1.10 ) exists. The following form of our main result shows that in any case such theory has to be interesting and non-trivial.

Let
${\Theta}_{1}$
be defined by the equations
${\Theta}_{1}=\{Z:\theta (Z)={\partial}_{1}\theta (Z)=0\}$
. The
${\partial}_{1}$
-invariant subset
$\Sigma $
of
${\Theta}_{1}$
will be called the singular locus.

Theorem 1.2
An indecomposable principally polarized abelian variety
$(X,\theta )$
is the Jacobian of a curve of genus g if and only if there exist
$g$
-dimensional vectors
$U\ne 0,V$
, such that for each
$Z\in {\mathbb{C}}^{g}\backslash \Sigma $
equation ( 1.10 ) for the function
$\tau (x,y)=\theta (Ux+Vy+Z)$
is satisfied, i.e.
The main idea of Shiota's proof of the Novikov's conjecture is to show that if
$u$
is as in ( 1.1 ) and satisfies the KP equation, then it can be extended to a
$\tau $
-function of the KP hierarchy, as a global holomorphic function of the infinite number of variables
$t=\left\{{t}_{i}\right\},{t}_{1}=x,{t}_{2}=y,{t}_{3}=t$
. Local existence of
$\tau $
directly follows from the KP equation. The global existence of the
$\tau $
-function is crucial. The rest is a corollary of the KP theory and the theory of commuting ordinary differential operators developed by Buchnall-Chaundy ([9, 10] ) and the author ([5, 6] ).

The core of the problem is that there is a homological obstruction for the global existence of
$\tau $
. It is controlled by the cohomology group
${H}^{1}({\mathbb{C}}^{g}\backslash \Sigma ,\mathcal{V})$
, where
$\mathcal{V}$
is the sheaf of
${\partial}_{1}$
-invariant meromorphic functions on
${\mathbb{C}}^{g}\backslash \Sigma $
with poles along
$\Theta $
(see details in [11] ). The hardest part of the Shiota's work (clarified in [11] ) is the proof that the locus
$\Sigma $
is empty.

That insures vanishing of
${H}^{1}({\mathbb{C}}^{g},\mathcal{V})$
. Analogous obstructions have occurred in all the other attempts to apply the theory of soliton equations to various characterization problems in the theory of abelian varieties. None of them had been completely successful. Only partial results were obtained. (Note that Theorem 1.1 in one of its equivalent forms was proved earlier in [12] under the additional assumption that
${\Theta}_{1}$
does not contain
${\partial}_{1}$
-invariant line.) Strictly speaking, the KP equation and the KP hierarchy are not used in the present paper. But our main construction of the formal wave solutions of ( 1.5 ) is reminiscent to the construction of the
$\tau $
-function. All its difficulties can be traced back to those in the Shiota's work. The wave solution of ( 1.5 ) is a solution of the form

At the beginning of the next section we show that the assumptions of Theorem 1.2 are necessary and sufficient conditions for the local existence of the wave solutions such that

where
${\tau}_{s}(Z,y)$
, as a function of
$Z$
, is holomorphic in some open domain in
${\mathbb{C}}^{g}$
. The functions
${\xi}_{s}$
are defined recurrently by the equation
$2{\partial}_{1}{\xi}_{s+1}={\partial}_{y}{\xi}_{s}-{\partial}_{1}^{2}{\xi}_{s}+u{\xi}_{s}$
. Therefore, the global existence of
${\xi}_{s}$
is controlled by the same cohomology group
${H}^{1}(\mathbb{C}\backslash \Sigma ,\mathcal{V})$
as above. At the local level the main problem is to find translational invariant normalization of
${\xi}_{s}$
which defines wave solutions uniquely up to a
${\partial}_{1}$
-invariant factor.

$$\begin{array}{c}\psi (x,y,k)={e}^{kx+({k}^{2}+b)y}\left(1+{\sum}_{s=1}^{\infty}{\xi}_{s}(x,y){k}^{-s}\right).\end{array}$$ | (1.16) |

$$\begin{array}{c}{\xi}_{s}=\frac{{\tau}_{s}(Ux+Vy+Z,y)}{\theta (Ux+Vy+Z)},Z/\in \Sigma ,\end{array}$$ | (1.17) |

In the case of periodic potentials
$u(x+T,y)=u\left(x\right)$
the normalization problem for the wave functions was solved D.Phong and the author in [13] . It was shown that the condition that
${\xi}_{s}$
is periodic completely determines the
$y$
-dependence of the integration constants and the corresponding wave solutions are related by a
$x$
-independent factor. In general, the potential
$u=-2{\partial}_{x}^{2}\theta (Ux+Vy+Z)$
is only quasi-periodic in
$x$
. In that case the solution of the normalization problem is technically more involved but mainly goes along the same lines as in the periodic case. The corresponding wave solutions are called
$\lambda $
-periodic.

In the last section we show that for each
$Z/\in \Sigma $
a local
$\lambda $
-periodic wave solution is the common eigenfunction of a commutative ring
${\mathcal{A}}^{Z}$
of ordinary differential operators. The coefficients of these operators are independent of ambiguities in the construction of
$\psi $
. The correspondence
$Z/\in \Sigma \mapsto {\mathcal{A}}^{Z}$
extends by linearity and defines a monomorphism
$j:X\to J(\Gamma )$
, where
$J(\Gamma )$
is the generalized Jacobian of the corresponding
$Z$
-independent spectral curve
$\Gamma $
of
${\mathcal{A}}^{Z}$
. The final step is the proof of the statement that the genus of the spectral curve is
$g$
. A'priory, from the construction it follows only that the arithmetic genus
$\hat{g}$
of
$\Gamma $
is not greater than
${2}^{g}$
. The proof of the equation
$g=\hat{g}$
does require the global existence of the wave solutions. In the framework of our approach this key problem is solved by taking the preimage
${j}^{*}\left({\psi}_{BA}\right)$
of the Baker-Akhiezer function on
$\Gamma $
which is known to be globally well-defined on
$J(\Gamma )$
.

2 $\lambda $ -periodic wave solutions

As it was mentioned above, formal Calogero-Moser equations ( 1.11 ) were derived in [7] as a necessary condition for the existence of a meromorphic solution to equation ( 1.5 ).

Let
$\tau (x,y)$
be a holomorphic function of the variable
$x$
in some open domain
$D\in \mathbb{C}$
smoothly depending on a parameter
$y$
. Suppose that for each
$y$
the zeros of
$\tau $
are simple,

(All coefficients in these expansions are smooth functions of the variable
$y$
). Substitution of ( 2.2 , 2.3 ) in ( 1.5 ) gives a system of equations. The first three of them are

Taking the
$y$
-derivative of the first equation and using two others we get ( 1.11 ).

$$\begin{array}{c}\tau \left({x}_{i}\right(y),y)=0,{\tau}_{x}\left({x}_{i}\right(y),y)\ne 0.\end{array}$$ | (2.1) |

Lemma 2.1
([7] ) If equation ( 1.5 ) with the potential
$u=-2{\partial}_{x}^{2}ln\tau (x,y)$
has a meromorphic in
$D$
solution
${\psi}_{0}(x,y)$
, then equations ( 1.11 ) hold.

Proof. Consider the Laurent expansions of
${\psi}_{0}$
and
$u$
in the neighborhood of one of the zeros
${x}_{i}$
of
$\tau $
:
$$\begin{array}{c}u=\frac{2}{(x-{x}_{i}{)}^{2}}+{v}_{i}+{w}_{i}(x-{x}_{i})+\dots \end{array}$$ | (2.2) |

$$\begin{array}{c}{\psi}_{0}=\frac{{\alpha}_{i}}{x-{x}_{i}}+{\beta}_{i}+{\gamma}_{i}(x-{x}_{i})+{\delta}_{i}(x-{x}_{i}{)}^{2}+\dots \end{array}$$ | (2.3) |

$$\begin{array}{c}{\alpha}_{i}{\dot{x}}_{i}+2{\beta}_{i}=0,\end{array}$$ | (2.4) |

$$\begin{array}{c}{\dot{\alpha}}_{i}+{\alpha}_{i}{v}_{i}+2{\gamma}_{i}=0,\end{array}$$ | (2.5) |

$$\begin{array}{c}{\dot{\beta}}_{i}+{v}_{i}{\beta}_{i}-{\gamma}_{i}{\dot{x}}_{i}+{\alpha}_{i}{w}_{i}=0.\end{array}$$ | (2.6) |

Let us show that equations ( 1.11 ) are sufficient for the existence of meromorphic wave solutions.

We are going to prove by induction that this system has meromorphic solutions with simple poles at all the zeros
${x}_{i}$
of
$\tau $
.

Lemma 2.2
Suppose that equations ( 1.11 ) for the zeros of
$\tau (x,y)$
hold. Then there exist meromorphic wave solutions of equation ( 1.5 ) that have simple poles at
${x}_{i}$
and are holomorphic everywhere else.

Proof. Substitution of ( 1.16 ) into ( 1.5 ) gives a recurrent system of equations
$$\begin{array}{c}2{\xi}_{s+1}^{\prime}={\partial}_{y}{\xi}_{s}+u{\xi}_{s}-{\xi}_{s}^{\prime \prime}\end{array}$$ | (2.7) |

Let us expand
${\xi}_{s}$
at
${x}_{i}$
:

where for brevity we omit the index
$i$
in the notations for the coefficients of this expansion.

$$\begin{array}{c}{\xi}_{s}=\frac{{r}_{s}}{x-{x}_{i}}+{r}_{s0}+{r}_{s1}(x-{x}_{i}),\end{array}$$ | (2.8) |

Suppose that
${\xi}_{s}$
are defined and equation ( 2.7 ) has a meromorphic solution. Then the right hand side of ( 2.7 ) has the zero residue at
$x={x}_{i}$
, i.e.,

We need to show that the residue of the next equation vanishes also. From ( 2.7 ) it follows that the coefficients of the Laurent expansion for
${\xi}_{s+1}$
are equal to

These equations imply

and the lemma is proved.

$$\begin{array}{c}re{s}_{{x}_{i}}\left({\partial}_{y}{\xi}_{s}+u{\xi}_{s}-{\xi}_{s}^{\prime \prime}\right)={\dot{r}}_{s}+{v}_{i}{r}_{s}+2{r}_{s1}=0\end{array}$$ | (2.9) |

$$\begin{array}{c}{r}_{s+1}=-{\dot{x}}_{i}{r}_{s}-2{r}_{s0},\end{array}$$ | (2.10) |

$$\begin{array}{c}2{r}_{s+1,1}={\dot{r}}_{s0}-{r}_{s1}+{w}_{i}{r}_{s}+{v}_{i}{r}_{s0}.\end{array}$$ | (2.11) |

$$\begin{array}{c}{\dot{r}}_{s+1}+{v}_{i}{r}_{s+1}+2{r}_{s+1,1}=-{r}_{s}({\ddot{x}}_{i}-2{w}_{i})-{\dot{x}}_{i}({\dot{r}}_{s}-{v}_{i}{r}_{s}s+2{r}_{s1})=0,\end{array}$$ | (2.12) |

Our next goal is to fix a translation-invariant normalization of
${\xi}_{s}$
which defines wave functions uniquely up to a
$x$
-independent factor. It is instructive to consider first the case of the periodic potentials
$u(x+1,y)=u(x,y)$
(see details in [13] ).

Equations ( 2.7 ) are solved recursively by the formulae

where
${c}_{s}\left(y\right)$
are arbitrary functions of the variable
$y$
. Let us show that the periodicity condition
${\xi}_{s}(x+1,y)={\xi}_{s}(x,y)$
defines the functions
${c}_{s}\left(y\right)$
uniquely up to an additive constant. Assume that
${\xi}_{s-1}$
is known and satisfies the condition that the corresponding function
${\xi}_{s}^{0}$
is periodic. The choice of the function
${c}_{s}\left(y\right)$
does not affect the periodicity property of
${\xi}_{s}$
, but it does affect the periodicity in
$x$
of the function
${\xi}_{s+1}^{0}(x,y)$
. In order to make
${\xi}_{s+1}^{0}(x,y)$
periodic, the function
${c}_{s}\left(y\right)$
should satisfy the linear differential equation

where
$B\left(y\right)={\int}_{{x}_{0}}^{{x}_{0}+1}udx$
. This defines
${c}_{s}$
uniquely up to a constant.

$$\begin{array}{c}{\xi}_{s+1}(x,y)={c}_{s+1}\left(y\right)+{\xi}_{s+1}^{0}(x,y),\end{array}$$ | (2.13) |

$$\begin{array}{c}{\xi}_{s+1}^{0}(x,y)=\frac{1}{2}{\int}_{{x}_{0}}^{x}({\partial}_{y}{\xi}_{s}-{\xi}_{s}^{\prime \prime}+u{\xi}_{s})dx,\end{array}$$ | (2.14) |

$$\begin{array}{c}{\partial}_{y}{c}_{s}\left(y\right)+B\left(y\right){c}_{s}\left(y\right)+{\int}_{{x}_{0}}^{{x}_{0}+1}\left({\partial}_{y}{\xi}_{s}^{0}(x,y)+u(x,y){\xi}_{s}^{0}(x,y)\right)dx,\end{array}$$ | (2.15) |

In the general case, when
$u$
is quasi-periodic, the normalization of the wave functions is defined along the same lines.

Let
${Y}_{U}=\langle Ux\rangle $
be the closure of the group
$Ux$
in
$X$
. Shifting
${Y}_{U}$
if needed, we may assume, without loss of generality, that
${Y}_{U}$
is not in the singular locus,
${Y}_{U}/\in \Sigma $
. Then, for a sufficiently small
$y$
, we have
${Y}_{U}+Vy/\in \Sigma $
as well. Consider the restriction of the theta-function onto the affine subspace
${\mathbb{C}}^{d}+Vy$
, where
${\mathbb{C}}^{d}={\pi}^{-1}\left({Y}_{U}\right)$
, and
$\pi :{\mathbb{C}}^{g}\to X={\mathbb{C}}^{g}/\Lambda $
is the universal cover of
$X$
:

The function
$u(z,y)=-2{\partial}_{1}^{2}ln\tau $
is periodic with respect to the lattice
${\Lambda}_{U}=\Lambda \cap {\mathbb{C}}^{d}$
and, for fixed
$y$
, has a double pole along the divisor
${\Theta}^{U}\left(y\right)=\left(\Theta -Vy\right)\cap {\mathbb{C}}^{d}$
.

A particular solution of the first equation
$2{\partial}_{1}{\xi}_{1}=u+b$
is given by the formula

where
$(l,z)$
is a linear form on
${\mathbb{C}}^{d}$
given by the scalar product of
$z$
with a vector
$l\in {\mathbb{C}}^{d}$
such that
$(l,U)=1$
, and
$(l,\lambda )=1.$
The periodicity condition for
${\xi}_{1}^{0}$
defines the constant
$b$

which depends only on a choice of the lattice vector
$\lambda $
. A change the potential by an additive constant does not affect the results of the previous lemma. Therefore, equations ( 1.11 ) are sufficient for the local solvability of ( 2.20 ) in any domain, where
$\tau (z+Ux,y)$
has simple zeros, i.e., outside of the set
${\Theta}_{1}^{U}\left(y\right)=\left({\Theta}_{1}-Vy\right)\cap {\mathbb{C}}^{d}$
. Recall that
${\Theta}_{1}=\Theta \cap {\partial}_{1}\Theta $
. This set does not contain a
${\partial}_{1}$
-invariant line because any such line is dense in
${Y}_{U}$
. Therefore, the sheaf
${\mathcal{V}}_{0}$
of
${\partial}_{1}$
-invariant meromorphic functions on
${\mathbb{C}}^{d}\backslash {\Theta}_{1}^{U}\left(y\right)$
with poles along the divisor
${\Theta}^{U}\left(y\right)$
coincides with the sheaf of holomorphic
${\partial}_{1}$
-invariant functions. That implies the vanishing of
${H}^{1}({C}^{d}\backslash {\Theta}_{1}^{U}(y),{\mathcal{V}}_{0})$
and the existence of global meromorphic solutions
${\xi}_{s}^{0}$
of ( 2.20 ) which have a simple pole at the divisor
${\Theta}^{U}\left(y\right)$
(see details in [1, 11] ). If
${\xi}_{s}^{0}$
are fixed, then the general global meromorphic solutions are given by the formula
${\xi}_{s}={\xi}_{s}^{0}+{c}_{s}$
, where the constant of integration
${c}_{s}(z,y)$
is a holomorphic
${\partial}_{1}$
-invariant function of the variable
$z$
.

$$\begin{array}{c}\tau (z,y)=\theta (z+Vy),z\in {\mathbb{C}}^{d}.\end{array}$$ | (2.16) |

Lemma 2.3
Let equation ( 1.10 ) for
$\tau (Ux+z,y)$
hold and let
$\lambda $
be a vector of the sublattice
${\Lambda}_{U}=\Lambda \cap {\mathbb{C}}^{d}\subset {\mathbb{C}}^{g}$
. Then: (i) equation ( 1.5 ) with the potential
$u(Ux+z,y)$
has a wave solution of the form
$\psi ={e}^{kx+{k}^{2}y}\phi (Ux+z,y,k)$
such that the coefficients
${\xi}_{s}(z,y)$
of the formal series

are
$\lambda $
-periodic meromorphic functions of the variable
$z\in {\mathbb{C}}^{d}$
with a simple pole at the divisor
${\Theta}^{U}\left(y\right)$
,

(ii)
$\phi (z,y,k)$
is unique up to a factor
$\rho (z,k)$
that is
${\partial}_{1}$
-invariant and holomorphic in
$z$
,

Proof. The functions
${\xi}_{s}\left(z\right)$
are defined recursively by the equations
$$\begin{array}{c}\phi (z,y,k)={e}^{by}\left(1+{\sum}_{s=1}^{\infty}{\xi}_{s}(z,y){k}^{-s}\right)\end{array}$$ | (2.17) |

$$\begin{array}{c}{\xi}_{s}(z+\lambda ,y)={\xi}_{s}(z,y)=\frac{{\tau}_{s}(z,y)}{\tau (z,y)};\end{array}$$ | (2.18) |

$$\begin{array}{c}{\phi}_{1}(z,y,k)=\phi (z,y,k)\rho (z,k),{\partial}_{1}\rho =0.\end{array}$$ | (2.19) |

$$\begin{array}{c}2{\partial}_{1}{\xi}_{s+1}={\partial}_{y}{\xi}_{s}+(u+b){\xi}_{s}-{\partial}_{1}^{2}{\xi}_{s}.\end{array}$$ | (2.20) |

$$\begin{array}{c}2{\xi}_{1}^{0}=-2{\partial}_{1}ln\tau +(l,z)b,\end{array}$$ | (2.21) |

$$\begin{array}{c}b=2{\partial}_{1}ln\tau (z+\lambda ,y)-2{\partial}_{1}ln\tau (z,y),\end{array}$$ | (2.22) |

Let us assume, as in the example above, that a
$\lambda $
-periodic solution
${\xi}_{s-1}$
is known and that it satisfies the condition that there exists a periodic solution
${\xi}_{s}^{0}$
of the next equation.

Let
${\xi}_{s+1}^{*}$
be a solution of ( 2.20 ) for fixed
${\xi}_{s}^{0}$
. Then it is easy to see that the function

is a solution of ( 2.20 ) for
${\xi}_{s}={\xi}_{s}^{0}+{c}_{s}$
. A choice of a
$\lambda $
-periodic
${\partial}_{1}$
-invariant function
${c}_{s}(z,y)$
does not affect the periodicity property of
${\xi}_{s}$
, but it does affect the periodicity of the function
${\xi}_{s+1}^{0}$
. In order to make
${\xi}_{s+1}^{0}$
periodic, the function
${c}_{s}(z,y)$
should satisfy the linear differential equation

This equation, together with an initial condition
${c}_{s}\left(z\right)={c}_{s}(z,0)$
uniquely defines
${c}_{s}(x,y)$
.

$$\begin{array}{c}{\xi}_{s+1}^{0}(z,y)={\xi}_{s+1}^{*}(z,y)+{c}_{s}(z,y){\xi}_{1}^{0}(z,y)+\frac{(l,z)}{2}{\partial}_{y}{c}_{s}(z,y),\end{array}$$ | (2.23) |

$$\begin{array}{c}{\partial}_{y}{c}_{s}(z,y)=2{\xi}_{s+1}^{*}(z+\lambda ,y)-2{\xi}_{s+1}^{*}(z,y).\end{array}$$ | (2.24) |

The induction step is then completed. We have shown that the ratio of two periodic formal series
${\phi}_{1}$
and
$\phi $
is
$y$
-independent. Therefore, equation ( 2.19 ), where
$\rho (z,k)$
is defined by the evaluation of the both sides at
$y=0$
, holds. The lemma is thus proven.

where the coefficients of
${\rho}_{\lambda}$
are
${\partial}_{1}$
-invariant holomorphic functions. Then the same arguments as in [1] show that there exists a
${\partial}_{1}$
-invariant series
$f(z,k)$
with holomorphic in
$z$
coefficients and formal series
${\mu}_{\lambda}^{0}\left(k\right)$
with constant coefficients such that the equation

holds. The ambiguity in the choice of
$f$
and
$\mu $
corresponds to the multiplication by the exponent of a linear form in
$z$
vanishing on
$U$
, i.e.

where
$b\left(k\right)={\sum}_{s}{b}_{s}{k}^{-s}$
is a formal series with vector-coefficients that are orthogonal to
$U$
.

Corollary 2.1
Let
${\lambda}_{1},\dots ,{\lambda}_{d}$
be a set of linear independent vectors of the lattice
${\Lambda}_{U}$
and let
${z}_{0}$
be a point of
${\mathbb{C}}^{d}$
. Then, under the assumptions of the previous lemma, there is a unique wave solution of equation ( 1.5 ) such that the corresponding formal series
$\phi (z,y,k;{z}_{0})$
is quasi-periodic with respect to
${\Lambda}_{U}$
, i.e. for
$\lambda \in {\Lambda}_{U}$

and satisfies the normalization conditions

The proof is identical to that of the part (b) of the Lemma 12 in [1] . Let us briefly present its main steps. As shown above there exist wave solutions corresponding to
$\phi $
which are
${\lambda}_{1}$
-periodic. Moreover, from the statement (ii) above it follows that for any
${\lambda}^{\prime}\in {\Lambda}_{U}$
$$\begin{array}{c}\phi (z+\lambda ,y,k;{z}_{0})=\phi (z,y,k;{z}_{0}){\mu}_{\lambda}\left(k\right)\end{array}$$ | (2.25) |

$$\begin{array}{c}{\mu}_{{\lambda}_{i}}\left(k\right)=1,\phi ({z}_{0},0,k;{z}_{0})=1.\end{array}$$ | (2.26) |

$$\begin{array}{c}\phi (z+\lambda ,y,k)=\phi (z,y,k){\rho}_{\lambda}(z,k),\end{array}$$ | (2.27) |

$$\begin{array}{c}f(z+\lambda ,k){\rho}_{\lambda}(z,k)=f(z,k){\mu}_{\lambda}\left(k\right)\end{array}$$ | (2.28) |

$$\begin{array}{c}{f}^{\prime}(z,k)=f(z,k){e}^{\left(b\right(k),z)},{\mu}_{\lambda}^{\prime}\left(k\right)={\mu}_{\lambda}\left(k\right){e}^{\left(b\right(k),\lambda )},\left(b\right(k),U)=0,\end{array}$$ | (2.29) |

The vector
$U$
is in general position with respect to the lattice. Therefore, the ambiguity can be uniquely fixed by imposing
$(d-1)$
normalizing conditions
${\mu}_{{\lambda}_{i}}\left(k\right)=1,i>1$
(recall that
${\mu}_{{\lambda}_{1}}\left(k\right)=1$
by construction).

The formal series
$f\phi $
is quasi-periodic and its multiplicators satisfy ( 2.26 ). Then, by that properties it is defined uniquely up to a factor which is constant in
$z$
and
$y$
. Therefore, for the unique definition of
${\phi}_{0}$
it is enough to fix its evaluation at
${z}_{0}$
and
$y=0$
. The corollary is proved.

3 The spectral curve

In this section we show that
$\lambda $
-periodic wave solutions of equation ( 1.5 ), with
$u$
as in ( 1.7 ), are common eigenfunctions of rings of commuting operators and identify
$X$
with the Jacobian of the spectral curve of these rings.

Note that a simple shift
$z\to z+Z$
, where
$Z/\in \Sigma ,$
gives
$\lambda $
-periodic wave solutions with meromorphic coefficients along the affine subspaces
$Z+{\mathbb{C}}^{d}$
. Theses
$\lambda $
-periodic wave solutions are related to each other by
${\partial}_{1}$
-invariant factor. Therefore choosing, in the neighborhood of any
$Z/\in \Sigma ,$
a hyperplane orthogonal to the vector
$U$
and fixing initial data on this hyperplane at
$y=0,$
we define the corresponding series
$\phi (z+Z,y,k)$
as a local meromorphic function of
$Z$
and the global meromorphic function of
$z$
.

Lemma 3.1
Let the assumptions of Theorem 1.2 hold. Then there is a unique pseudo-differential operator

such that

where
$\psi ={e}^{kx+{k}^{2}y}\phi (Ux+Z,y,k)$
is a
$\lambda $
-periodic solution of ( 1.5 ). The coefficients
${w}_{s}\left(Z\right)$
of
$\mathcal{\mathcal{L}}$
are meromorphic functions on the abelian variety
$X$
with poles along the divisor
$\Theta $
.
Proof. The construction of
$\mathcal{\mathcal{L}}$
is standard for the KP theory. First we define
$\mathcal{\mathcal{L}}$
as a pseudo-differential operator with coefficients
${w}_{s}(Z,y)$
, which are functions of
$Z$
and
$y$
.

$$\begin{array}{c}\mathcal{\mathcal{L}}(Z,{\partial}_{x})={\partial}_{x}+{\sum}_{s=1}^{\infty}{w}_{s}\left(Z\right){\partial}_{x}^{-s}\end{array}$$ | (3.1) |

$$\begin{array}{c}\mathcal{\mathcal{L}}(Ux+Vy+Z,{\partial}_{x})\psi =k\psi ,\end{array}$$ | (3.2) |

Let
$\psi $
be a
$\lambda $
-periodic wave solution. The substitution of ( 2.17 ) in ( 3.2 ) gives a system of equations that recursively define
${w}_{s}(Z,y)$
as differential polynomials in
${\xi}_{s}(Z,y)$
. The coefficients of
$\psi $
are local meromorphic functions of
$Z$
, but the coefficients of
$\mathcal{\mathcal{L}}$
are well-defined global meromorphic functions of on
${\mathbb{C}}^{g}\backslash \Sigma $
, because different
$\lambda $
-periodic wave solutions are related to each other by
${\partial}_{1}$
-invariant factor, which does not affect
$\mathcal{\mathcal{L}}$
. The singular locus is of codimension
$\ge 2$
. Then Hartog's holomorphic extension theorem implies that
${w}_{s}(Z,y)$
can be extended to a global meromorpic function on
${\mathbb{C}}^{g}$
.

The translational invariance of
$u$
implies the translational invariance of the
$\lambda $
-periodic wave solutions. Indeed, for any constant
$s$
the series
$\phi (Vs+Z,y-s,k)$
and
$\phi (Z,y,k)$
correspond to
$\lambda $
-periodic solutions of the same equation. Therefore, they coincide up to a
${\partial}_{1}$
-invariant factor. This factor does not affect
$\mathcal{\mathcal{L}}$
. Hence,
${w}_{s}(Z,y)={w}_{s}(Vy+Z)$
.

The
$\lambda $
-periodic wave functions corresponding to
$Z$
and
$Z+{\lambda}^{\prime}$
for any
${\lambda}^{\prime}\in \Lambda $
are also related to each other by a
${\partial}_{1}$
-invariant factor:

Hence,
${w}_{s}$
are periodic with respect to
$\Lambda $
and therefore are meromorphic functions on the abelian variety
$X$
. The lemma is proved.

$$\begin{array}{c}{\partial}_{1}\left({\phi}_{1}(Z+{\lambda}^{\prime},y,k){\phi}^{-1}(Z,y,k)\right)=0.\end{array}$$ | (3.3) |

Consider now the differential parts of the pseudo-differential operators
${\mathcal{\mathcal{L}}}^{m}$
. Let
${\mathcal{\mathcal{L}}}_{+}^{m}$
be the differential operator such that
${\mathcal{\mathcal{L}}}_{-}^{m}={\mathcal{\mathcal{L}}}^{m}-{\mathcal{\mathcal{L}}}_{+}^{m}={F}_{m}{\partial}^{-1}+O\left({\partial}^{-2}\right)$
. The leading coefficient
${F}_{m}$
of
${\mathcal{\mathcal{L}}}_{-}^{m}$
is the residue of
${\mathcal{\mathcal{L}}}^{m}$
:

From the construction of
$\mathcal{\mathcal{L}}$
it follows that
$[{\partial}_{y}-{\partial}_{x}^{2}+u,{\mathcal{\mathcal{L}}}^{n}]=0$
. Hence,

The functions
${F}_{m}$
are differential polynomials in the coefficients
${w}_{s}$
of
$\mathcal{\mathcal{L}}$
. Hence,
${F}_{m}\left(Z\right)$
are meromorphic functions on
$X$
. Next statement is crucial for the proof of the existence of commuting differential operators associated with
$u$
.

The coefficients of
$\Phi $
are universal differential polynomials on
${\xi}_{s}$
. Therefore,
${\phi}_{s}(z+Z,y)$
is a global meromorphic function of
$z\in {C}^{d}$
and a local meromorphic function of
$Z/\in \Sigma $
. Note that
$\mathcal{\mathcal{L}}=\Phi \left({\partial}_{x}\right){\Phi}^{-1}$
.

$$\begin{array}{c}{F}_{m}=re{s}_{\partial}{\mathcal{\mathcal{L}}}^{m}.\end{array}$$ | (3.4) |

$$\begin{array}{c}[{\partial}_{y}-{\partial}_{x}^{2}+u,{\mathcal{\mathcal{L}}}_{+}^{m}]=-[{\partial}_{y}-{\partial}_{x}^{2}+u,{\mathcal{\mathcal{L}}}_{-}^{m}]=2{\partial}_{x}{F}_{m}.\end{array}$$ | (3.5) |

Lemma 3.2
The abelian functions
${F}_{m}$
have at most the second order pole on the divisor
$\Theta $
.

Proof. We need a few more standard constructions from the KP theory. If
$\psi $
is as in Lemma 3.1, then there exists a unique pseudo-differential operator
$\Phi $
such that
$$\begin{array}{c}\psi =\Phi {e}^{kx+{k}^{2}y},\Phi =1+{\sum}_{s=1}^{\infty}{\phi}_{s}(Ux+Z,y){\partial}_{x}^{-s}.\end{array}$$ | (3.6) |

Consider the dual wave function defined by the left action of the operator
${\Phi}^{-1}$
:
${\psi}^{+}=\left({e}^{-kx-{k}^{2}y}\right){\Phi}^{-1}$
. Recall that the left action of a pseudo-differential operator is the formal adjoint action under which the left action of
${\partial}_{x}$
on a function
$f$
is
$\left(f{\partial}_{x}\right)=-{\partial}_{x}f$
. If
$\psi $
is a formal wave solution of ( 3.5 ), then
${\psi}^{+}$
is a solution of the adjoint equation

The same arguments, as before, prove that if equations ( 1.11 ) for poles of
$u$
hold then
${\xi}_{s}^{+}$
have simple poles at the poles of
$u$
. Therefore, if
$\psi $
as in Lemma 2.3, then the dual wave solution is of the form
${\psi}^{+}={e}^{-kx-{k}^{2}y}{\phi}^{+}(Ux+Z,y,k)$
, where the coefficients
${\xi}_{s}^{+}(z+Z,y)$
of the formal series

are
$\lambda $
-periodic meromorphic functions of the variable
$z\in {\mathbb{C}}^{d}$
with a simple pole at the divisor
${\Theta}^{U}\left(y\right)$
.

$$\begin{array}{c}(-{\partial}_{y}-{\partial}_{x}^{2}+u){\psi}^{+}=0.\end{array}$$ | (3.7) |

$$\begin{array}{c}{\phi}^{+}(z+Z,y,k)={e}^{-by}\left(1+{\sum}_{s=1}^{\infty}{\xi}_{s}^{+}(z+Z,y){k}^{-s}\right)\end{array}$$ | (3.8) |

The ambiguity in the definition of
$\psi $
does not affect the product

Therefore, although each factor is only a local meromorphic function on
${\mathbb{C}}^{g}\backslash \Sigma $
, the coefficients
${J}_{s}$
of the product

are global meromorphic functions of
$Z$
. Moreover, the translational invariance of
$u$
implies that they have the form
${J}_{s}(Z,y)={J}_{s}(Z+Vy)$
. Each of the factors in the left hand side of ( 3.10 ) has a simple pole on
$\Theta -Vy$
. Hence,
${J}_{s}\left(Z\right)$
is a meromorphic function on
$X$
with a second order pole at
$\Theta $
.

$$\begin{array}{c}{\psi}^{+}\psi =\left({e}^{-kx-{k}^{2}y}{\Phi}^{-1}\right)\left(\Phi {e}^{kx+{k}^{2}y}\right).\end{array}$$ | (3.9) |

$$\begin{array}{c}{\psi}^{+}\psi ={\phi}^{+}(Z,y,k)\phi (Z,y,k)=1+{\sum}_{s=2}^{\infty}{J}_{s}(Z,y){k}^{-s}.\end{array}$$ | (3.10) |

From the definition of
$\mathcal{\mathcal{L}}$
it follows that

On the other hand, using the identity

which holds for any two pseudo-differential operators ([14] ), we get

Therefore,
${F}_{n}={J}_{n+1}$
and the lemma is proved.

$$\begin{array}{c}re{s}_{k}\left({\psi}^{+}\left({\mathcal{\mathcal{L}}}^{n}\psi \right)\right)=re{s}_{k}\left({\psi}^{+}{k}^{n}\psi \right)={J}_{n+1}.\end{array}$$ | (3.11) |

$$\begin{array}{c}re{s}_{k}\left({e}^{-kx}{\mathcal{D}}_{1}\right)\left({\mathcal{D}}_{2}{e}^{kx}\right)=re{s}_{\partial}\left({\mathcal{D}}_{2}{\mathcal{D}}_{1}\right),\end{array}$$ | (3.12) |

$$\begin{array}{c}re{s}_{k}\left({\psi}^{+}{\mathcal{\mathcal{L}}}^{n}\psi \right)=re{s}_{k}\left({e}^{-kx}{\Phi}^{-1}\right)\left({\mathcal{\mathcal{L}}}^{n}\Phi {e}^{kx}\right)=re{s}_{\partial}{\mathcal{\mathcal{L}}}^{n}={F}_{n}.\end{array}$$ | (3.13) |

Let
$\hat{\mathbf{F}}$
be a linear space generated by
$\{{F}_{m},m=0,1,\dots \}$
, where we set
${F}_{0}=1$
. It is a subspace of the
${2}^{g}$
-dimensional space of the abelian functions that have at most second order pole at
$\Theta $
. Therefore, for all but
$\hat{g}=dim\hat{\mathbf{F}}$
positive integers
$n$
, there exist constants
${c}_{i,n}$
such that

Let
$I$
denote the subset of integers
$n$
for which there are no such constants. We call this subset the gap sequence.

Hence, if
$\psi $
is a
$\lambda $
-periodic wave solution of ( 1.5 ) corresponding to
$Z/\in \Sigma $
, then
${L}_{n}\psi $
is also a formal solution of the same equation. That implies the equation
${L}_{n}\psi ={a}_{n}(Z,k)\psi $
, where
$a$
is
${\partial}_{1}$
-invariant. The ambiguity in the definition of
$\psi $
does not affect
${a}_{n}$
. Therefore, the coefficients of
${a}_{n}$
are well-defined global meromorphic functions on
${\mathbb{C}}^{g}\backslash \Sigma $
. The
${\partial}_{1}$
invariance of
${a}_{n}$
implies that
${a}_{n}$
, as a function of
$Z$
, is holomorphic outside of the locus. Hence it has an extension to a holomorphic function on
${\mathbb{C}}^{g}$
. Equations ( 3.3 ) imply that
${a}_{n}$
is periodic with respect to the lattice
$\Lambda $
. Hence
${a}_{n}$
is
$Z$
-independent. Note that
${a}_{s,n}={c}_{s,n},s\le n$
. The lemma is proved.

$$\begin{array}{c}{F}_{n}\left(Z\right)+{\sum}_{i=0}^{n-1}{c}_{i,n}{F}_{i}\left(Z\right)=0.\end{array}$$ | (3.14) |

Lemma 3.3
Let
$\mathcal{\mathcal{L}}$
be the pseudo-differential operator corresponding to a
$\lambda $
-periodic wave function
$\psi $
constructed above. Then, for the differential operators

the equations

where
${a}_{s,n}$
are constants, hold.

Proof. First note that from ( 3.5 ) it follows that
$$\begin{array}{c}{L}_{n}={\mathcal{\mathcal{L}}}_{+}^{n}+{\sum}_{i=0}^{n-1}{c}_{i,n}{\mathcal{\mathcal{L}}}_{+}^{n-i}=0,n/\in I,\end{array}$$ | (3.15) |

$$\begin{array}{c}{L}_{n}\psi ={a}_{n}\left(k\right)\psi ,{a}_{n}\left(k\right)={k}^{n}+{\sum}_{s=1}^{\infty}{a}_{s,n}{k}^{n-s}\end{array}$$ | (3.16) |

$$\begin{array}{c}[{\partial}_{y}-{\partial}_{x}^{2}+u,{L}_{n}]=0.\end{array}$$ | (3.17) |

The operator
${L}_{m}$
can be regarded as a
$Z/\in \Sigma $
-parametric family of ordinary differential operators
${L}_{m}^{Z}$
whose coefficients have the form

$$\begin{array}{c}{L}_{m}^{Z}={\partial}_{x}^{n}+{\sum}_{i=1}^{m}{u}_{i,m}(Ux+Z){\partial}_{x}^{m-i},m/\in I.\end{array}$$ | (3.18) |

Corollary 3.1
The operators
${L}_{m}^{Z}$
commute with each other,

From ( 3.16 ) it follows that
$[{L}_{n}^{Z},{L}_{m}^{Z}]\psi =0$
. The commutator is an ordinary differential operator. Hence, the last equation implies ( 3.19 ). $$\begin{array}{c}[{L}_{n}^{Z},{L}_{m}^{Z}]=0,Z/\in \Sigma .\end{array}$$ | (3.19) |

Remark. A change of the local coordinates
${k}^{\prime}=k+O\left({k}^{-1}\right)$
acts on the space of wave solutions in an obvious way. In order to simplify notations below it is convenient to normalized the local coordinate as follows. Let us choose an integer
$n/\in I$
and choose
$k$
such that the corresponding eigenvalue
${a}_{n}\left(k\right)$
in ( 3.16 ) is equal to

The correspondence becomes one-to-one if the rings
$\mathcal{A}$
are considered modulo conjugation
${\mathcal{A}}^{\prime}=g\left(x\right)\mathcal{A}{g}^{-1}\left(x\right)$
. Note that in [9, 10, 5, 6] the main attention was paid to the generic case of the commutative rings corresponding to smooth algebraic curves. The invariant formulation of the correspondence given above is due to Mumford [15] .

$$\begin{array}{c}{a}_{n}\left(k\right)={k}^{n}.\end{array}$$ | (3.20) |

Lemma 3.4
Let
${\mathcal{A}}^{Z},Z/\in \Sigma ,$
be a commutative ring of ordinary differential operators spanned by the operators
${L}_{n}^{Z}$
. Then there is an irreducible algebraic curve
$\Gamma $
with a fixed smooth point
${P}_{0}$
such that
${\mathcal{A}}^{Z}$
is isomorphic to the ring
$A(\Gamma ,{P}_{0})$
of the meromorphic functions on
$\Gamma $
with the only pole at the puncture
${P}_{0}$
. The arithmetic genus of
$\Gamma $
is equal to
$\hat{g}=dim\hat{\mathbf{F}}$
. The correspondence
$Z\to {\mathcal{A}}^{Z}$
defines an imbedding of
$X$
into the generalized Jacobian of
$\Gamma $
,

It is the fundamental fact of the theory of commuting differential operators [9, 10, 5, 6] that there is a natural correspondence between the set of commutative rings
$\mathcal{A}$
of ordinary linear differential operators containing a pair of monic operators of co-prime orders and the set of algebro-geometric data
$\{\Gamma ,{P}_{0},[{k}^{-1}{]}_{1},\mathcal{\mathcal{F}}\}$
, where
$\Gamma $
is an algebraic curve with fixed first jet
$[{k}^{-1}{]}_{1}$
of a local coordinate
${k}^{-1}$
in the neighborhood of a smooth point
${P}_{0}\in \Gamma $
and
$\mathcal{\mathcal{F}}$
is a torsion free rank 1 sheaf on
$\Gamma $
such that
$$\begin{array}{c}j:X\mapsto J(\Gamma ).\end{array}$$ | (3.21) |

$$\begin{array}{c}{H}^{0}(\Gamma ,\mathcal{\mathcal{F}})=1,{H}^{1}(\Gamma ,\mathcal{\mathcal{F}})=0.\end{array}$$ | (3.22) |

The algebraic curve corresponding to a commutative ring
$\mathcal{A}$
is called the spectral curve of
$\mathcal{A}$
. The ring
$\mathcal{A}$
is isomorphic to the ring
$A(\Gamma ,{P}_{0})$
of meromorphic functions on
$\Gamma $
with the only pole at the puncture
${P}_{0}$
. If the ring
$\mathcal{A}$
is maximal, i.e. is not contained in any bigger commutative ring, then in the general case the sheaf
$\mathcal{\mathcal{F}}$
is a line bundle and is defined by a non-special effective divisor
$D$
of degree
$\hat{g}$
which is equal to the arithmetic genus of
$\Gamma $
.

Our next goal is to show that the rings
${\mathcal{A}}^{Z}$
are maximal and that the corresponding spectral curve is
$Z$
-independent, i.e.
$\Gamma ={\Gamma}^{Z}$
. To begin with, let us briefly recall the key steps of the construction of the algebro-geometric spectral transform.

Let
${L}_{n}$
and
${L}_{m}$
be a pair of monic commuting differential operators of co-prime orders,
$(n,m)=1$
. The commutativity equation
$[{L}_{n},{L}_{m}]=0$
is invariant under gauge transformations of the form
${L}_{n}^{\prime}=g{L}_{n}{g}^{-1},{L}_{m}^{\prime}=g{L}_{m}{g}^{-1}$
, where
$g=g\left(x\right)$
. Therefore, without loss of generality we may assume that
${L}_{n}$
and
${L}_{m}$
have the form

The commutativity equations are translational invariant. Therefore, without loss of generality, we may assume also that
${u}_{i},{v}_{j}$
are smooth functions of
$x$
in the neighborhood of
${x}_{0}=0$
.

$$\begin{array}{c}{L}_{n}={\partial}_{x}^{n}+{\sum}_{i=2}^{n}{u}_{i}\left(x\right){\partial}_{x}^{n-i},{L}_{m}={\partial}_{x}^{m}+{\sum}_{j=1}^{m}{v}_{j}\left(x\right){\partial}_{x}^{n-i}.\end{array}$$ | (3.23) |

Then there is a polynomial
$R$
of the form

such that
$R({L}_{m},{L}_{n})=0$
.

$$\begin{array}{c}R(\mu ,E)={\mu}^{n}+{E}^{m}+{\sum}_{mi+nj<nm}{r}_{ij}{\mu}^{i}{E}^{j}\end{array}$$ | (3.24) |

The polynomial
$R$
is constructed as follows. Consider the space
$V\left(E\right)$
of solutions of the ordinary differential equation
${L}_{n}c=Ec$
, where
$E$
is a complex number. It is a linear space of dimension
$n$
. The space
$V\left(E\right)$
is invariant for
${L}_{m}$
. The restriction of
${L}_{m}$
onto
$V\left(E\right)$
can be seen as a finite-dimensional linear operator
${L}_{m}\left(E\right):V\left(E\right)\mapsto V\left(E\right)$
. Then
$R$
is the characteristic polynomial of
${L}_{m}\left(E\right)$
, i.e.
$R(\mu ,E)=det\left(\mu \cdot Id-{L}_{m}\left(E\right)\right)$
. By definition it is a polynomial of degree
$n$
in the variable
$\mu $
. In order to prove that
$R$
is polynomial in
$E$
it is enough to consider the canonical basis
${c}_{i}(x,E;{x}_{0})\in V\left(E\right)$
defined by the initial conditions
${\partial}_{x}^{i}{c}_{j}(x,E;{x}_{0}){|}_{x={x}_{0}}={\delta}_{j}^{i},i,j=0,\dots ,n-1$
. In the basis
${c}_{j}$
the operator
${L}_{m}\left(E\right)$
is represented by a matrix
${L}_{m}^{ij}(E;{x}_{0})$
whose entries are polynomials in
$E$
(see details in [5, 6] ). Hence,
$R$
is polynomial in both the variables. It is independent of a choice of the normalization point
${x}_{0}$
.

Consider now the formal wave solution of the equation
${L}_{n}{\psi}_{0}(x,k;{x}_{0})={k}^{n}{\psi}_{0}(x,{k}_{0};{x}_{0})$

normalized by the condition
${\psi}_{0}({x}_{0},k;{x}_{0})=1$
. If
${L}_{n}$
is as in ( 3.23 ), then
${\psi}_{0}$
exists and unique.

$$\begin{array}{c}{\psi}_{0}={e}^{k(x-{x}_{0})}\left(1+{\sum}_{s=1}^{\infty}{\xi}_{s}(x;{x}_{0}){k}^{-s}\right),\end{array}$$ | (3.25) |

The coefficients
${\xi}_{s}$
are defined recursively by a first order differential equation and the initial condition
${\xi}_{s}({x}_{0};{x}_{0})=0$
. The uniqueness of
${\psi}_{0}$
implies that if
${L}_{m}$
commutes with
${L}_{n}$
, then there is a formal series
${a}_{m}\left(k\right)$
with constant coefficients such that

Let
${V}_{\infty}\left(E\right)$
be the linear space spanned by the formal solutions
${\psi}_{0}(x,{\varepsilon}_{j}k;{x}_{0})$
over the field of Laurent series in
${k}^{-1},{\varepsilon}_{j}^{n}=1,{k}^{n}=E$
. In this space one can introduce the basis satisfying the same initial conditions as the canonical basis
${c}_{j}$
above. In this basis
${L}_{n}$
is represented by the same matrix
${L}_{n}^{ij}(E;{x}_{0})$
. Hence,

From ( 3.27 ) it follows that
$R$
has the form ( 3.24 ). Equation ( 3.24 ) implies that if
$n$
and
$m$
are co-prime, then the affine curve defined in
${C}^{2}$
by the equation
$R(\mu ,E)=0$
is compactified at the infinity by one smooth point
${P}_{0}$
. A branch
${k}^{-1}$
of
${E}^{-1/n}$
is a local coordinate in the neighborhood of
${P}_{0}$
(compare with ( 3.20 )).

$$\begin{array}{c}{L}_{m}{\psi}_{0}={a}_{m}\left(k\right){\psi}_{0},{a}_{m}={k}^{m}+{\sum}_{s=1}^{\infty}{a}_{s,m}{k}^{-s}.\end{array}$$ | (3.26) |

$$\begin{array}{c}R(\mu ,E){=}^{n-1}{\prod}_{j=0}(\mu -{a}_{m}({\varepsilon}_{j}k\left)\right).\end{array}$$ | (3.27) |

From ( 3.24 ) it follows that if
$E$
is big enough, then the operator
${L}_{m}\left(E\right)$
has
$n$
distinct eigenvalues. Hence, it has simple eigenvalues for almost all
$E$
. Therefore, for almost all points
$P=(\mu ,E)$
of the curve
${\Gamma}_{0}$
there is a unique normalized eigenvector
$f(P;{x}_{0})=\left\{{f}_{i}\right(P;{x}_{0}\left)\right\}$
,

The coordinates
${f}_{i},1\le i\le n-1,$
are rational functions of the variables
$E$
and
$\mu $
. Therefore, if
${\Gamma}_{0}$
is smooth, then
$f(P,{x}_{0})$
extends to all points of
${\Gamma}_{0}$
as a meromorphic vector-function.

$$\begin{array}{c}{L}_{m}(E;{x}_{0})f(P;{x}_{0})=\mu f(P;{x}_{0}),{f}_{0}(P,{x}_{0})=1.\end{array}$$ | (3.28) |

In the general case, when the spectral curve is singular, its extensions on branches of
${\Gamma}_{0}$
passing through a singular point might have different values at the singular point.

The function
${\psi}_{0}(x,P;{x}_{0})={\sum}_{j=0}^{n-1}{f}_{i}\left(P\right){c}_{j}(x,E;{x}_{0})$
is the solution of the equations

normalized by the condition
${\psi}_{0}({x}_{0},P;{x}_{0})=1$
. From the definition of the canonical basis it follows that

The first coordinate of
$f$
in the neighborhood of the puncture
${P}_{0}$
has a simple pole
${f}_{1}=k+O\left({k}^{-1}\right)$
. Hence,
${\psi}_{0}$
in the neighborhood of
${P}_{0}$
has the form ( 3.25 ).

$$\begin{array}{c}{L}_{n}{\psi}_{0}(x,P;{x}_{0})=E\left(P\right){\psi}_{0}(x,P;{x}_{0}),{L}_{m}{\psi}_{0}(x,P)=\mu \left(P\right){\psi}_{0}(x,P;{x}_{0}),\end{array}$$ | (3.29) |

$$\begin{array}{c}{\psi}_{0}(x,P;{x}_{0})=exp\left({\int}_{{x}_{0}}^{x}{f}_{1}(P;x)dx\right).\end{array}$$ | (3.30) |

The basis functions
${c}_{j}$
are entire functions of the variable
$E$
. Therefore, if
${\Gamma}_{0}$
is smooth, then
${\psi}_{0}$
is meromorphic on
$\Gamma \backslash {P}_{0}$
and the divisor
$D$
of its poles coincides with the pole divisor of
$f$
. Hence, it is
$x$
-independent. Let us show that the degree of
$D$
is equal to the genus
${g}_{0}$
of
${\Gamma}_{0}$
.

The curve
${\Gamma}_{0}$
is realized as an
$n$
-sheeted cover of the
$E$
-plane. For a generic
$E$
the equation
$R(\mu ,E)=0$
has
$n$
different roots. Let
${P}_{j}\left(E\right)=\left({\mu}_{j}\right(E),E)$
be the corresponding points of
${\Gamma}_{0}$
. Consider the square of the Wronskian of the functions
$\left\{{\psi}_{0}\right(x,{P}_{j}\left(E\right);{x}_{0}\left)\right\}$
, i.e.,
${\Delta}_{{\psi}_{0}}(E;{x}_{0})={det}^{2}\left({\partial}_{x}^{i}{\psi}_{0}(x;{P}_{j}(E);{x}_{0})\right)$
. It does not depend on
$x$
and an ordering of the points
${P}_{j}\left(E\right)$
. Therefore,
${\Delta}_{{\psi}_{0}}(E;{x}_{0})$
is a meromoprhic function of the variable
$E$
. It has double poles at the projections of the poles of
${\psi}_{0}$
and a pole of order
$n-1$
at the infinity
$E=\infty $
. Its zeros are the projections of the branch points of the affine part
${\Gamma}_{0}$
over
$E$
.

The multiplicities
${\nu}_{s}$
of the zeros are equal to the multiplicities of the branch points. Hence
$2degD+n-1={\sum}_{s}{\nu}_{s}$
. The infinity is the branch point of the multiplicity
$n-1$
. Then the Riemann-Hurwitz formula implies
$degD={g}_{0}$
.

The analytical properties of
${\psi}_{0}$
described above are the defining properties of Baker-Akhiezer functions. The same simple arguments as the ones used for the counting of the degree of
$D$
prove the uniqueness of the Baker-Akhiezer functions. Suppose that there exists a function
${\psi}_{1}$
which has the form ( 3.25 ) in the neighborhood of
${P}_{0}$
and is meromorphic on
${\Gamma}_{0}\backslash {P}_{0}$
, where it has the same divisor of poles
$D$
as
${\psi}_{0}$
. The function
$\stackrel{~}{\psi}={\psi}_{1}-{\psi}_{0}$
has at
${P}_{0}$
the form
$\stackrel{~}{\psi}={e}^{k(x-{x}_{0})}O\left({k}^{-1}\right)$
. Therefore, at the infinity the function
${\Delta}_{\stackrel{~}{\psi}}$
has a pole of order less or equal than
$n-2$
. Hence, the total number of its zeros is greater than the total number of poles. That implies
${\Delta}_{\stackrel{~}{\psi}}=0$
and therefore
$\stackrel{~}{\psi}=0$
.

The uniqueness of the Baker-Akhiezer function implies that for each meromorphic function
$a\left(P\right)\in A(\Gamma ,{P}_{0})$
there is a unique operator
${L}_{a}$
such that

Indeed, for each
$a$
there exists a unique operator
${L}_{a}$
such that in the neighborhood of
${P}_{0}$
the following congruence
$({L}_{a}-a){\psi}_{0}={e}^{kx}O\left({k}^{-1}\right)$
holds. The function
$({L}_{a}-a){\psi}_{0}$
is meromorphic outside of
${P}_{0}$
and has the same divisor of poles
$D$
. As it was shown above that implies ( 3.31 ).

$$\begin{array}{c}{L}_{a}{\psi}_{0}=a{\psi}_{0},a\in A(\Gamma ,{P}_{0}).\end{array}$$ | (3.31) |

The degree of the operator
${L}_{a}$
is equal to the order of the pole of
$a\left(P\right)$
at
${P}_{0}$
. Therefore, by the Weierstrass gap theorem the number of integers
$m$
such that the ring of the operators
${L}_{a}$
does not contain an operator of order
$m$
is equal to the genus
${g}_{0}$
of
${\Gamma}_{0}$
. The set of such integers will be called the gap sequence of the commutative ring.

Consider now the case when
${\Gamma}_{0}$
is singular. There are almost no changes in the previous consideration if
${\psi}_{0}$
extends to a single-valued meromorphic function outside of
${P}_{0}$
, i.e.
${\psi}_{0}\in {H}^{0}({\Gamma}_{0}\backslash {P}_{0},\mathcal{O}(D\left)\right)$
. Note that
${\psi}_{0}$
can not be equal to zero identically at any point, because
$f\ne 0$
. Therefore, changing if needed the normalization point
${x}_{0}$
, we can always assume without loss of generality that all the points of
$D$
are smooth. Then the counting of the degree of
$D$
gives the equation
$degD+n-1={\sum}_{s}{\nu}_{s}+{\sum}_{j}{\mu}_{j}$
, where
${\nu}_{s}$
are multiplicities of the branch points and
${\mu}_{j}$
are multiplicities of the singular points of
${\Gamma}_{0}$
. Therefore,
$degD={g}_{0}$
, where
${g}_{0}$
is the arithmetic genus of
${\Gamma}_{0}$
. As in the smooth case the corresponding operators
${L}_{a}$
span a maximal commutative ring. The length of its gap sequence is equal to the arithmetic genus
${g}_{0}$
of
${\Gamma}_{0}$
.

If
${\psi}_{0}$
is not a single-valued meromorphic function outside of
${P}_{0}$
, then it can be regarded as a section of a rank 1 sheaf,
${\psi}_{0}\in {H}^{0}({\Gamma}_{0}\backslash {P}_{0},\mathcal{\mathcal{F}})$
. If that case the multiplicity of the zero of
${\Delta}_{{\psi}_{0}}$
at a singular point might be less than the multiplicity
${\mu}_{j}$
of the point. That implies
$degD<{g}_{0}$
. This inequality corresponds to the case, when the ring generated by
${L}_{n}$
and
${L}_{m}$
is not maximal.

Let
${L}_{n},{L}_{m},{L}_{i},i=1,\dots ,r,$
be generators of the maximal commutative ring. The spectral curve
$\Gamma $
of this ring is defined in
${\mathbb{C}}^{r+2}$
by equations
$R(\mu ,E)=0$
and
${R}_{i}({\mu}_{i},E)=0$
, where
${R}_{i}$
are characteristic polynomials for the operators
${L}_{i}\left(E\right)$
. The projection of
${\mathbb{C}}^{2+r}$
onto
${\mathbb{C}}^{2}$
defines a regular map
$\Gamma \to {\Gamma}_{0}$
, which is one-to-one away from the singular points. The common eigenfunction
${\psi}_{0}$
of the commuting operators becomes a single-valued meromorphic function on
$\Gamma \backslash {P}_{0}$
. The degree of its pole divisor is equal to the arithmetic genus
$\stackrel{~}{g}$
of
$\Gamma $
.

Now we are in the position to present a proof of the lemma. Let
${\hat{\mathcal{A}}}^{Z}$
be the maximum commutative ring containing
${\mathcal{A}}^{Z}$
. The arithmetic genus of the corresponding spectral curve
${\Gamma}^{Z}$
is equal to the degree of the pole divisor of the
$\lambda $
-periodic wave solution
$\psi $
. It is continuous in
$Z$
. Hence, all spectral curves
${\Gamma}^{Z}$
are of the same arithmetic genus
$\stackrel{~}{g}$
. That implies that the lengths of the gap sequences
${I}^{Z}$
of the rings
${\hat{\mathcal{A}}}^{Z}$
are all equal. Let us show that the gap sequences by themselves are
$Z$
-independent. For that it is enough to show that
${I}^{Z}$
is locally
$Z$
-independent.

Suppose that the gap sequences are not constant in the neighborhood of a point
${Z}_{0}\in X\backslash \Sigma $
. Then there is a sequence of points
${Z}_{i}\in X\backslash \Sigma $
, which converges to
${Z}_{0}$
, such that
${I}^{{Z}_{i}}={I}^{{Z}_{j}}\ne {I}^{{Z}_{0}}$
. The lengths of
${I}^{Z}$
are equal to each other. Hence, there exists
$m\in {I}^{{Z}_{0}}$
suchthat
$m/\in {I}^{{Z}_{i}}$
. The last condition means that there are operators
${L}_{m}^{\left(i\right)}$
of order
$m$
such that
$[{L}_{n}^{{Z}_{i}},{L}_{m}^{\left(i\right)}]=0$
. Note that
${L}_{m}^{\left(i\right)}$
belongs to
${\hat{\mathcal{A}}}^{{Z}_{i}}$
but not necessary to
${\mathcal{A}}^{{Z}_{i}}$
. The coefficients of the commuting operators are bounded holomorphic functions at the neighborhood of some point
${x}_{0}$
. Therefore, we may choose a subsequence such that the operators
${L}_{m}^{\left(i\right)}$
converge to an operator
${L}_{m}^{0}$
. By continuity, this operator commutes with
${L}_{n}^{{Z}_{0}}$
, and this contradicts to the assumption that
$m\in {I}^{{Z}_{0}}$
.

Let us show that
${\mathcal{A}}^{Z}$
is a maximal commutative ring. Suppose that it is not. Then there exits
$\alpha /\in {I}_{0}$
, where
${I}_{0}$
is the gap sequence of
${\mathcal{A}}^{Z}$
, such that for each
$Z/\in \Sigma $
there exist an operator
${L}_{\alpha}^{Z}$
of order
$\alpha $
which commutes with
${L}_{n}^{Z},n/\in {I}_{0}$
. Therefore, it commutes with
$\mathcal{\mathcal{L}}$
.

A differential operator commuting with
$\mathcal{\mathcal{L}}$
up to the order
$O\left(1\right)$
can be represented in the form
${L}_{\alpha}={\sum}_{m<\alpha}{c}_{i,\alpha}\left(Z\right){\mathcal{\mathcal{L}}}_{+}^{i}$
, where the constants
${c}_{i,\alpha}\left(Z\right)$
are
${\partial}_{1}$
-invariant. It commutes with
$\mathcal{\mathcal{L}}$
if and only if

Note the difference between ( 3.14 ) and ( 3.32 ). In the first equation the coefficients
${c}_{i,n}$
are constants. The
$\lambda $
-periodic wave solution of equation ( 1.5 ) is a common eigenfunction of all commuting operators, i.e.
${L}_{\alpha}\psi ={a}_{\alpha}(Z,k)\psi $
, where
${a}_{\alpha}={k}^{\alpha}+{\sum}_{s=1}^{\infty}{a}_{s,\alpha}\left(Z\right){k}^{\alpha -s}$
is
${\partial}_{1}$
-invariant. The same arguments as the ones used in the proof of equation ( 3.16 ) show that the eigenvalue
${a}_{\alpha}$
is
$Z$
-independent. We have
${a}_{s,\alpha}={c}_{s,\alpha},s\le \alpha $
. Therefore, the coefficients in ( 3.32 ) are
$Z$
-independent. That contradicts to the assumption that
${\mathcal{A}}^{Z}$
is not maximal.

$$\begin{array}{c}{F}_{\alpha}\left(Z\right)+{\sum}_{i=0}^{n-1}{c}_{i,\alpha}\left(Z\right){F}_{i}\left(Z\right)=0,{\partial}_{1}{c}_{i,\alpha}=0.\end{array}$$ | (3.32) |

Let
$\Gamma $
be a spectral curve of the rings
${\mathcal{A}}^{Z}$
. It is
$Z$
-independent because the eigenvalues of all operators
${L}_{n}^{Z}\in {\mathcal{A}}^{Z}$
are
$Z$
-independent. The set of all maximal commuting rings corresponding to
$\Gamma $
is parameterized by points of the generalized Jacobian
$J(\Gamma )$
, the open set of which is just the set of non-special divisors of degree
$\hat{g}=dim\hat{\mathbf{F}}.$
In terms of the KP theory, the set of maximal commuting operators corresponding to a fixed algebraic curve with a puncture is identified with the orbit of the KP hierarchy. It is finite-dimensional and a choice of the basic flows defines an affine structure on the universal cover
${\mathbb{C}}^{\hat{g}}\mapsto J(\Gamma )$
.

The KP hierarchy in the Sato form is a system of commuting differential equation for a pseudo-differential operator
$\mathcal{\mathcal{L}}$

If the operator
$\mathcal{\mathcal{L}}$
is as above. i.e., if it is defined by
$\lambda $
-periodic wave solutions equation ( 1.5 ), then equations ( 3.33 ) are equivalent to the equations

The first two times of the hierarchy are identified with the variables
${t}_{1}=x,{t}_{2}=y$
. Equations ( 3.14 ) imply that the corresponding linear combination of the flows is trivial. The flows defined by linearly independent vector-fields
${\partial}_{\alpha}={\partial}_{{t}_{\alpha}},\alpha \in I,$
acts transitively on
$J(\Gamma )$
(see details in [5, 6, 1, 16] ).

$$\begin{array}{c}{\partial}_{{t}_{n}}\mathcal{\mathcal{L}}=[{\mathcal{\mathcal{L}}}_{+}^{n},\mathcal{\mathcal{L}}].\end{array}$$ | (3.33) |

$$\begin{array}{c}{\partial}_{{t}_{n}}u={\partial}_{x}{F}_{n}.\end{array}$$ | (3.34) |

The correspondence
$Z\mapsto {A}^{Z}$
defines an imbedding of
$j:X\backslash \Sigma \mapsto J(\Gamma )$
. Under this imbedding the vector-fields
${\partial}_{{z}_{i}}=\partial /\partial {z}_{i}$
, where
$Z=({z}_{1},\dots ,{z}_{g})$
, map to tangent vector-fields on
$J(\Gamma )$
.

Therefore,

In terms of
$u$
the deformation with respect to
${z}_{i}$
is given by the equation

Comparing with ( 3.34 ), we get that

Note that a'priori the coefficients
${h}_{\beta ,i}$
are
${\partial}_{1}$
-invariant functions of the variable
$Z$
. They are constants because
${H}_{i}$
, as well as
${F}_{\beta}$
, are meromorphic functions on
$X$
with a double pole along
$\Theta $
. Equation ( 3.35 ) implies that the imbedding
$j:X\backslash \Sigma \mapsto J(\Gamma )$
extends by linearity onto
$X$
. The lemma is proved.

$$\begin{array}{c}{j}_{*}\left({\partial}_{{z}_{i}}\right)={\sum}_{\beta \in I}{h}_{\beta ,i}{\partial}_{\beta}.\end{array}$$ | (3.35) |

$$\begin{array}{c}{\partial}_{{z}_{i}}u={\partial}_{x}{H}_{i},{H}_{i}={\partial}_{x}{\partial}_{{z}_{i}}ln\theta \left(Z\right).\end{array}$$ | (3.36) |

$$\begin{array}{c}{H}_{i}\left(Z\right)={\sum}_{\beta \in I}{h}_{\beta ,i}{F}_{\beta}\left(Z\right)\mapsto {H}_{i}\in \mathbf{F}.\end{array}$$ | (3.37) |

Proof of Theorem 1.2. Lemma 3.4 implies that if
$\hat{g}=dim\mathbf{F}=g$
, then the linear imbedding ( 3.21 ) is an isomorphism. For the proof of the last equation it is enough to show that all functions
${F}_{n}$
have the form
${F}_{n}={\partial}_{x}{Q}_{n}$
, where
${Q}_{n}$
is a global meromorphic function with poles along
$\Theta $
. The space of such functions
${Q}_{n}$
that have a simple pole along
$\Theta $
is
$g$
-dimensional and is spanned by the functions
${\partial}_{{z}_{i}}ln\theta .$
It is an easy matter to show that outside of the singular locus
$\Sigma $
each function
${F}_{n}$
is locally a
${\partial}_{1}$
-derivative of a meromorphic function with a pole along
$\Theta $
. Indeed, as shown above
${F}_{n}={J}_{n+1}$
, where
${J}_{s}$
are the coefficients of expansion ( 3.10 ) for
${\psi}^{+}\psi $
. The left and right action of pseudo-differential operators are formally adjoint, i.e., for any two operators the equality
$\left({e}^{-kx}{\mathcal{D}}_{1}\right)\left({\mathcal{D}}_{2}{e}^{kx}\right)={e}^{-kx}\left({\mathcal{D}}_{1}{\mathcal{D}}_{2}{e}^{kx}\right)+{\partial}_{x}\left({e}^{-kx}\left({\mathcal{D}}_{3}{e}^{kx}\right)\right)$
holds. Here
${\mathcal{D}}_{3}$
is a pseudo-differential operator whose coefficients are differential polynomials in the coefficients of
${\mathcal{D}}_{1}$
and
${\mathcal{D}}_{2}$
. Therefore, from ( 3.9 ) it follows that

The coefficients of the series
$Q$
are differential polynomials in the coefficients
${\xi}_{l},{\xi}_{l}^{+},l\le s,$
of
$\psi $
and
${\psi}^{+}$
. Recall that the
$\lambda $
-periodic wave functions are defined only locally on
$X\backslash \Sigma $
.

$$\begin{array}{c}{\psi}^{+}\psi =1+{\sum}_{s=2}^{\infty}{F}_{s-1}{k}^{-s}=1+{\partial}_{x}\left({\sum}_{s=2}^{\infty}{Q}_{s}{k}^{-s}\right).\end{array}$$ | (3.38) |

To complete the proof of the theorem, it remains to show that there exists a
$\lambda $
-periodic wave solution which is defined globally. The existence of the imbedding ( 3.21 ) allows us to make that crucial step.

It is instructive to consider first the case when the spectral curve
$\Gamma $
of the rings
${\mathcal{A}}^{Z}$
is smooth. The corresponding Baker-Akhiezer function is of the form ([5, 6] )

Here
$\hat{\theta}\left(\hat{Z}\right)$
is the Riemann theta-function constructed with the help of the matrix of
$b$
-periods of normalized holomorphic differentials on
$\Gamma $
;
$\hat{A}:\Gamma \to J(\Gamma )$
is the Abel map;
$\Omega $
is the abelian integral corresponding to the second kind meromorphic differential
$d\Omega $
with the only pole of the form
$dk$
at the puncture
${P}_{0}$
. The last factors in the numerator and the denominator are
$x$
-independent. Therefore, the function

is also a common eigenfunction of the commuting operators. Its image
$\psi ={j}^{*}\left(\hat{\psi}\right)$
is the desired wave solution of ( 1.5 ) (at
$y=0$
).

$$\begin{array}{c}{\psi}_{0}=\frac{\hat{\theta}\left(\hat{A}\right(P)+\hat{U}x+\hat{Z})\hat{\theta}\left(\hat{Z}\right)}{\hat{\theta}(\hat{U}x+\hat{Z})\hat{\theta}\left(\hat{A}\right(P)+\hat{Z})}{e}^{x\Omega \left(P\right)}.\end{array}$$ | (3.39) |

$$\begin{array}{c}\hat{\psi}=\frac{\hat{\theta}\left(\hat{A}\right(P)+\hat{U}x+\hat{Z})}{\hat{\theta}(\hat{U}x+\hat{Z})}{e}^{x\Omega \left(P\right)}\end{array}$$ | (3.40) |

In the general case when
$\Gamma $
might be singular, the Riemann theta-function has to be replaced by the corresponding
$\tau $
-function, which is known to be holomorphic, see [16] . The rest of the arguments is the same as in the smooth case. The theorem is thus proved.

Acknowledgments. The author would like to thank Enrico Arbarello for very useful conversations and help during the preparation of this paper.

References

- T. Shiota, Characterization of Jacobian varieties in terms of soliton equations, Invent. Math., 83(2):333–382, 1986.
- B.A. Dubrovin, Kadomtsev-Petvishvili equation and relations for matrix of periods on Riemann surfaces, Izvestiya AN SSSR, ser. mat., 45 (1981), 1015-1028.
- Zakharov V.E., Shabat A.B. Integration method of nonlinear equations of mathematical physics with the help of inverse scattering problem, Funk. anal i pril. 8 (1974), . 3, 43-53.
- V.S. Druma, On analytic solution of the two-dimensional Korteweg-de Vries equation, JETP Letters 19 (1974), . 12, 219-225.
- I. M. Krichever, Integration of non-linear equations by methods of algebraic geometry, Funct. Anal. Appl., 11 (1977), n. 1, 12–26.
- I. M. Krichever, Methods of algebraic geometry in the theory of non-linear equations, Russian Math. Surveys, 32 (1977), n. 6, 185–213.
- E. Arbarello, I. Krichever, G. Marini, Characterizing Jacobians via flexes of the Kummer Variety, math.AG/0502138.
- G. E. Welters, A criterion for Jacobi varieties, Ann. of Math., 120 (1984), n. 3, 497–504.
- J.L. Burchnall, T.W. Chaundy, Commutative ordinary differential operators.I, Proc. London Math Soc. 21 (1922), 420–440.
- J.L. Burchnall, T.W. Chaundy, Commutative ordinary differential operators.II, Proc. Royal Soc.London 118 (1928), 557–583.
- E.Arbarello and C. De Concini, Another proof of a conjecture of S.P. Novikov on periods of abelian integrals on Riemann surfaces, Duke Math. Journal, 54 (1987), 163–178.
- G. Marini, A geometrical proof of Shiota's theorem on a conjecture of S.P. Novikov, Compositio Mathematica, 111:305–322, 1998.
- I. Krichever and D.H. Phong, Symplectic forms in the theory of solitons, Surveys in Differential Geometry IV (1998), edited by C.L. Terng and K. Uhlenbeck, 239-313, International Press, hep-th/9708170;
- L.A.Dickey, Soliton equations and Hamiltonian systems, Advanced Series in Mathematical Physics, Vol. 12 (1991) World Scientific, Singapore.
- D. Mumford, An algebro-geometric construction of commuting operators and of solutions to the Toda lattice equation, Korteweg-de Vries equation and related non-linear equations – Proceedings Int.Symp. Algebraic Geometry, Kyoto, 1977, 115–153, Kinokuniya Book Store, Kyoto, 1978.
- G.Sigal, G.Wilson, Loop groups and equations of KdV type, IHES Publ.Math. 61, 1985, 5-65.