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

$$\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)

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

$$\begin{array}{c}{\ddot{x}}_i=2w_i,\end{array}$$

(1.11)

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.

$$\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)

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

$$\begin{array}{c}\tau (x,y)=c\left(y\right){\prod }_i(x-x_i(y\left)\right),\end{array}$$

(1.13)

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

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

(1.14)

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.

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

$$\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)

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$  .

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.

equation ( 1.15 ) is valid on $\Theta$  .

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

$$\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)

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

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

(1.17)

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.

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,

$$\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)

(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

$$\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)

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

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

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)

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

Let us expand ${\xi }_s$  at $x_i$  :

$$\begin{array}{c}{\xi }_s=\frac{r_s}{x-x_i}+r_{s0}+r_{s1}(x-x_i),\end{array}$$

(2.8)

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

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.,

$$\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+2r_{s1}=0\end{array}$$

(2.9)

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

$$\begin{array}{c}{r}_{s+1}=-{\dot{x}}_i{r}_s-2r_{s0},\end{array}$$

(2.10)

$$\begin{array}{c}2r_{s+1,1}={\dot{r}}_{s0}-r_{s1}+w_i{r}_s+v_i{r}_{s0}.\end{array}$$

(2.11)

These equations imply

$$\begin{array}{c}{\dot{r}}_{s+1}+v_i{r}_{s+1}+2r_{s+1,1}=-r_s({\ddot{x}}_i-2w_i)-{\dot{x}}_i({\dot{r}}_s-v_i{r}_{s}s+2r_{s1})=0,\end{array}$$

(2.12)

and the lemma is proved.

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

$$\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)

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

$$\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)

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

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$  :

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

(2.16)

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$  .

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

$$\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)

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}{\xi }_s(z+\lambda ,y)={\xi }_s(z,y)=\frac{{\tau }_s(z,y)}{\tau (z,y)};\end{array}$$

(2.18)

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

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

(2.19)

Proof. The functions ${\xi }_s\left(z\right)$  are defined recursively by the equations

$$\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)

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

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

(2.21)

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$ 

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

(2.22)

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$  .

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

$$\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)

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

$$\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)

This equation, together with an initial condition $c_s\left(z\right)=c_s(z,0)$  uniquely defines $c_s(x,y)$  .

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.

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$ 

$$\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)

and satisfies the normalization conditions

$$\begin{array}{c}{\mu }_{{\lambda }_i}\left(k\right)=1,\phi (z_0,0,k;z_0)=1.\end{array}$$

(2.26)

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)=\phi (z,y,k){\rho }_{\lambda }(z,k),\end{array}$$

(2.27)

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

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

(2.28)

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.

$$\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)

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

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

$$\begin{array}{c}\mathcal{ℒ}(Z,{\partial }_x)={\partial }_x+{\sum }_{s=1}^{\infty }{w}_s\left(Z\right){\partial }_x^{-s}\end{array}$$

(3.1)

such that

$$\begin{array}{c}\mathcal{ℒ}(Ux+Vy+Z,{\partial }_x)\psi =k\psi ,\end{array}$$

(3.2)

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{ℒ}$  are meromorphic functions on the abelian variety $X$  with poles along the divisor $\Theta$  .

Proof. The construction of $\mathcal{ℒ}$  is standard for the KP theory. First we define $\mathcal{ℒ}$  as a pseudo-differential operator with coefficients $w_s(Z,y)$  , which are functions of $Z$  and $y$  .

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{ℒ}$  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{ℒ}$  . 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{ℒ}$  . 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:

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

(3.3)

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

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

$$\begin{array}{c}{F}_m=re{s}_{\partial }{\mathcal{ℒ}}^m.\end{array}$$

(3.4)

From the construction of $\mathcal{ℒ}$  it follows that $[{\partial }_y-{\partial }_x^2+u,{\mathcal{ℒ}}^n]=0$  . Hence,

$$\begin{array}{c}[{\partial }_y-{\partial }_x^2+u,{\mathcal{ℒ}}_{+}^m]=-[{\partial }_y-{\partial }_x^2+u,{\mathcal{ℒ}}_{-}^m]=2{\partial }_x{F}_m.\end{array}$$

(3.5)

The functions $F_m$  are differential polynomials in the coefficients $w_s$  of $\mathcal{ℒ}$  . 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$  .

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)

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{ℒ}=\Phi \left({\partial }_x\right){\Phi }^{-1}$  .

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

$$\begin{array}{c}(-{\partial }_y-{\partial }_x^2+u){\psi }^{+}=0.\end{array}$$

(3.7)

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

$$\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)

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)$  .

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

$$\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)

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

$$\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)

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$  .

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

$$\begin{array}{c}re{s}_k\left({\psi }^{+}\left({\mathcal{ℒ}}^n\psi \right)\right)=re{s}_k\left({\psi }^{+}{k}^n\psi \right)=J_{n+1}.\end{array}$$

(3.11)

On the other hand, using the identity

$$\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)

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

$$\begin{array}{c}re{s}_k\left({\psi }^{+}{\mathcal{ℒ}}^n\psi \right)=re{s}_k\left(e^{-kx}{\Phi }^{-1}\right)\left({\mathcal{ℒ}}^n\Phi e^{kx}\right)=re{s}_{\partial }{\mathcal{ℒ}}^n=F_n.\end{array}$$

(3.13)

Therefore, $F_n=J_{n+1}$  and the lemma is proved.

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

$$\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)

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

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

$$\begin{array}{c}{L}_n={\mathcal{ℒ}}_{+}^n+{\sum }_{i=0}^{n-1}{c}_{i,n}{\mathcal{ℒ}}_{+}^{n-i}=0,n/\in I,\end{array}$$

(3.15)

the equations

$$\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)

where $a_{s,n}$  are constants, hold.

Proof. First note that from ( 3.5 ) it follows that

$$\begin{array}{c}[{\partial }_y-{\partial }_x^2+u,L_n]=0.\end{array}$$

(3.17)

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.

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,

$$\begin{array}{c}[L_n^Z,L_m^Z]=0,Z/\in \Sigma .\end{array}$$

(3.19)

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 ).

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

$$\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$  ,

$$\begin{array}{c}j:X\mapsto J(\Gamma ).\end{array}$$

(3.21)

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{ℱ}\}$  , 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{ℱ}$  is a torsion free rank 1 sheaf on $\Gamma$  such that

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

(3.22)

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] .

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{ℱ}$  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

$$\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)

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$  .

Then there is a polynomial $R$  of the form

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

(3.24)

such that $R(L_m,L_n)=0$  .

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)$ 

$$\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)

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.

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

$$\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)

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,

$$\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.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 )).

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\}$  ,

$$\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)

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.

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

$$\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)

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

$$\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 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 ).

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 $\tilde{\psi }={\psi }_1-{\psi }_0$  has at $P_0$  the form $\tilde{\psi }=e^{k(x-x_0)}O\left(k^{-1}\right)$  . Therefore, at the infinity the function ${\Delta }_{\tilde{\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 }_{\tilde{\psi }}=0$  and therefore $\tilde{\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

$$\begin{array}{c}{L}_a{\psi }_0=a{\psi }_0,a\in A(\Gamma ,P_0).\end{array}$$

(3.31)

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 ).

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{ℱ})$  . 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 $\tilde{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 $\tilde{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{ℒ}$  .

A differential operator commuting with $\mathcal{ℒ}$  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{ℒ}}_{+}^i$  , where the constants $c_{i,\alpha }\left(Z\right)$  are ${\partial }_1$  -invariant. It commutes with $\mathcal{ℒ}$  if and only if

$$\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)

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.

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{ℒ}$ 

$$\begin{array}{c}{\partial }_{t_n}\mathcal{ℒ}=[{\mathcal{ℒ}}_{+}^n,\mathcal{ℒ}].\end{array}$$

(3.33)

If the operator $\mathcal{ℒ}$  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

$$\begin{array}{c}{\partial }_{t_n}u={\partial }_x{F}_n.\end{array}$$

(3.34)

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] ).

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,

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

(3.35)

In terms of $u$  the deformation with respect to $z_i$  is given by the equation

$$\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)

Comparing with ( 3.34 ), we get that

$$\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)

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.

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

$$\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)

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$  .

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] )

$$\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)

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

$$\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)

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$  ).

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

1. T. Shiota, Characterization of Jacobian varieties in terms of soliton equations, Invent. Math., 83(2):333–382, 1986.
2. B.A. Dubrovin, Kadomtsev-Petvishvili equation and relations for matrix of periods on Riemann surfaces, Izvestiya AN SSSR, ser. mat., 45 (1981), 1015-1028.
3. 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.
4. V.S. Druma, On analytic solution of the two-dimensional Korteweg-de Vries equation, JETP Letters 19 (1974), . 12, 219-225.
5. I. M. Krichever, Integration of non-linear equations by methods of algebraic geometry, Funct. Anal. Appl., 11 (1977), n. 1, 12–26.
6. I. M. Krichever, Methods of algebraic geometry in the theory of non-linear equations, Russian Math. Surveys, 32 (1977), n. 6, 185–213.
7. E. Arbarello, I. Krichever, G. Marini, Characterizing Jacobians via flexes of the Kummer Variety, math.AG/0502138.
8. G. E. Welters, A criterion for Jacobi varieties, Ann. of Math., 120 (1984), n. 3, 497–504.
9. J.L. Burchnall, T.W. Chaundy, Commutative ordinary differential operators.I, Proc. London Math Soc. 21 (1922), 420–440.
10. J.L. Burchnall, T.W. Chaundy, Commutative ordinary differential operators.II, Proc. Royal Soc.London 118 (1928), 557–583.
11. 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.
12. G. Marini, A geometrical proof of Shiota's theorem on a conjecture of S.P. Novikov, Compositio Mathematica, 111:305–322, 1998.
13. 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;
14. L.A.Dickey, Soliton equations and Hamiltonian systems, Advanced Series in Mathematical Physics, Vol. 12 (1991) World Scientific, Singapore.
15. 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.
16. G.Sigal, G.Wilson, Loop groups and equations of KdV type, IHES Publ.Math. 61, 1985, 5-65.