Matrix commuting differential operators of rank 2 Vardan Oganesyan

arXiv:1610.02367v1 [math-ph] 7 Oct 2016

Abstract In this paper we propose a very effective method for constructing matrix commuting differential operators of rank 2 and vector rank (2,2). We find new matrix commuting differential operators L, M of orders 2 and 2g respectively. Introduction Let us consider two differential operators Ln =

n X

ui (x)∂xi ,

Lm =

i=0

m X

vi (x)∂xi ,

i=0

where coefficients ui (x) and vi (x) are scalar or matrix valued functions. The commutativity condition Ln Lm = Lm Ln is equivalent to a very complicated system of nonlinear differential equations. The theory of commuting ordinary differential operators was first developed in the beginning of the XX century in the works of Wallenberg [1], Schur [2]. If two differential operators with scalar or matrix valued coefficients commute, then there exists a nonzero polynomial R(z, w) such that R(Ln , Lm ) = 0 (see [3], [4]). The curve Γ defined by R(z, w) = 0 is called the spectral curve. If Ln ψ = zψ,

Lm ψ = wψ,

then (z, w) ∈ Γ. If coefficients are scalar functions, then for almost all (z, w) ∈ Γ, the dimension of the space of common eigenfunctions ψ is the same. The dimension of the space of common eigenfunctions of two commuting scalar differential operators is called the rank of this pair. The rank is a common divisor of m and n. The genus of the spectral curve of a pair of commuting operators is called the genus of this pair. This research was supported by the Russian Science Foundation under grant 16-11-10260 and was done at the Faculty of Mechanics and Mathematics, Department of Geometry and Topology of Lomonosov Moscow State University

1

If the rank of two commuting scalar differential operators equals 1, then there are explicit formulas for coefficients of commutative operators in terms of Riemann theta-functions (see [5]). The case when rank of scalar commuting operators is greater than 1 is much more difficult. The first examples of commuting scalar differential operators of the nontrivial rank 2 and the nontrivial genus g = 1 were constructed by Dixmier [6] for the nonsingular elliptic spectral curve. A general classification of commuting scalar differential operators was obtained by Krichever [7]. The general form of commuting scalar operators of rank 2 for an arbitrary elliptic spectral curve was found by Krichever and Novikov [8]. The general form of scalar commuting operators of rank 3 with arbitrary elliptic spectral curve was found by Mokhov [9],[10]. In [11] Mironov developed theory of self-adjoint scalar operators of rank 2 and found examples of commuting scalar operators of rank 2 and arbitrary genus. Using Mironov’s method many examples of scalar commuting operators of rank 2 and arbitrary genus were found (see [12], [13], [14], [15], [16]). Moreover, examples of commuting scalar differential operators of arbitrary genus and arbitrary rank with polynomial coefficients were constructed by Mokhov in [17], [18]. Theory of commuting differential operators helps to find solutions of nonlinear partial differential equations from mathematical physics (see [19], [20], [21], [22]). Also there are deep connections between theory of commuting scalar differential operators and Schottky problem(see [23], [24]). The theory of commuting differential operators with polynomial coefficients has connections with the Dixmier conjecture and Jacobian conjecture (see [25], [26]). A general classification of commuting matrix differential operators was obtained by Grinevich [4]. Grinevich considered two differential operators L=

m X

Ui ∂xi ,

M=

i=0

n X

Vi ∂xi ,

i=0

where Ui and Vi are smooth and complex-valued s × s matrices. Let us suppose the following conditions 1) det(Um ) 6= 0. 2) Eigenvalues λ1 (x), ..., λs (x) of Um are distinct. 3) Matrix Vn is diagonalizable. Let µ1 (x), ..., µs (x) be eigenvalues of matrix Vn . Supµm pose that functions in are distinct constants for all i = 1, ..., s. λi Easy to see that if L and M commute, then Um and Vn commute. If operators 2

L and M commute, then F LF −1 and F M F −1 commute, where F is matrix. We also can change variable. So, without loss of generality we can suppose that (Um )ij = δij λi ,

(Vn )ij = δij µi ,

tr(Um−1 ) = 0

Let Γ be the spectral curve of commuting matrix operators L, M . Spectral curve of matrix commuting operators can be reducible. Let Γi be an irreducible component of the spectral curve. The dimension of the space of common eigenfunctions Lψ = zψ,

M ψ = wψ,

(z, w) ∈ Γi

is called the rank of commuting pair on ΓS i . Grinevich discovered that the spectral curve Γ has s points at infinity. So, Γ = ki=1 Γi , where k 6 s. Let li be the rank of operators on Γi . Operators L, M are called commuting operators of vector rank (l1 , ..., lk ), where k 6 s. Numbers li are common divisors of m and n. For more details see [4]. Also see [27], [28], [29]. If the rank of commuting matrix differential operators equals 1, then there exists explicit formulas for coefficients in terms of Riemann theta-functions [30]. In this paper we propose a very effective method for constructing matrix commuting operators of rank 2 and vector rank (2,2). We find new commuting operators L, M of orders 2 and 2g respectively.

Acknowledgments The author wishes to express gratitude to Professor O. I. Mokhov for advices and help in writing this paper. Explicit examples of commuting matrix differential operators of rank 2 Let us consider the operator L = E(x)∂x2 + R(x)∂ + Q(x),

(1)

where   λ1 (x) λ3 (x) E= , 0 λ2 (x)

  r1 (x) r3 (x) R= , r2 (x) −r1 (x)

3

  q1 (x) q3 (x) Q= q2 (x) q4 (x)

We want to find an operator M of order 2g such that [L, M ] = 0. Let us consider the operator M = B0 (x)Lg + (A1 (x) ∂x + B1 (x)) Lg−1 + (A2 (x)∂x + B2 (x))Lg−2 + ...+ +(Ag−1 (x)∂x + Bg−1 (x))L + Ag (x)∂x + Bg (x),

(2)

where Ag−k

  g−k a1 (x) ag−k 3 (x) , = g−k ag−k 2 (x) a4 (x)

 Bg−k =

 g−k bg−k 1 (x) b3 (x) , g−k bg−k 2 (x) b4 (x)

 A0 =

 0 0 . 0 0

We note that in the formulas above the number g − k is index not degree. The index g − k says us that functions aig−k and big−k are elements of matrices Ag−k and Bg−k respectively. Let us note that if [L, N ] = 0, then [L, M N ] = LM N −M N L = LM N −M LN = [L, M ]N , where L, M, N are matrix differential operators. So, we see that g g X X k [L, M ] = [L, (Ag−k ∂x + Bg−k )L ] = [L, (Ag−k ∂x + Bg−k )Lk ] = k=0

k=0 g

=

X

[L, (Ag−k ∂x + Bg−k )]Lk .

k=0

Direct calculations show that   L, Ag−k ∂x + Bg−k ] = [E∂x2 + R∂ + Q, Ag−k ∂x + Bg−k = +(EAg−k − Ag−k E)∂x3 + (2EA0g−k + EBg−k + RAg−k − Ag−k E 0 − Ag−k R − Bg−k E)∂x2 + 0 +(EA00g−k + 2EBg−k + Rg−k A0 + RBg−k + QAg−k − Ag−k R0 − Ag−k Q − Bg−k R)∂x + 00 0 +(EBg−k + RBg−k + QBg−k − Ag−k Q0 − Bg−k Q) =

= Kg−k ∂x3 + Pg−k ∂x2 + Tg−k ∂x + Fg−k , where

Tg−k =

Kg−k = EAg−k − Ag−k E, 0 Pg−k = 2EAg−k + EBg−k + RAg−k − Ag−k E 0 − Ag−k R − Bg−k E, 0 EA00g−k + 2EBg−k + Rg−k A0 + RBg−k + QAg−k − Ag−k R0 − Ag−k Q − 00 0 Fg−k = EBg−k + RBg−k + QBg−k − Ag−k Q0 − Bg−k Q. 4

Bg−k R,

Using the fact that ∂x L = E∂x3 + (E 0 + R)∂x2 + (R0 + Q)∂x + Q0 we get Kg−k ∂x3 + Pg−k ∂x2 + Tg−k ∂x + Fg−k =
Kg−k E −1 ∂x L + Pg−k − Kg−k E −1 (E 0 + R) ∂x2 +
Tg−k − Kg−k E −1 (R0 + Q) ∂x + (Fg−k − Kg−k E −1 Q0 ) =
Kg−k E −1 ∂x L + Pg−k − Kg−k E −1 (E 0 + R) E −1 L+

Tg−k − Kg−k E −1 (R0 + Q) − Pg−k − Kg−k E −1 (E 0 + R) E −1 R ∂x + +

Fg−l − Kg−k E −1 Q0 − Pg−k − Kg−k E −1 (E 0 + R) E −1 Q = e g−k ∂x L + Peg−k L + Teg−k ∂x + Feg−k , K where

e g−k = Kg−k E −1 K
Peg−k = Pg−k − Kg−k E −1 (E 0 + R) E −1 Teg−k = Tg−k − Kg−k E −1 (R0 + Q) − Kg−k E −1 (E 0 + R)E −1 R
Feg−k = Fg−k − Kg−k E −1 Q0 − Kg−k E −1 (E 0 + R)E −1 Q .




(3)

Finally we obtain [L, M ] =

g X

[L, (Ag−k ∂x + Bg−k )]Lk =

k=0

   g+1 e e e e e e + (T0 + K1 )∂x + (F0 + P1 ) Lg−1 + = K0 ∂x + P0 L   e 2 )∂x + (Fe1 + Pe2 ) Lg−1 + ...+ + (Te1 + K   e g )∂x + (Feg−1 + Peg ) L + Tg ∂x + Fg . + (Teg−1 + K 

So, if e 0 = 0, Pe0 = 0, K e 1 = −Te0 , Pe1 = −Fe0 , ..., K e m = −Tem−1 , Pem = −Fem−1 , ..., K e g = −Teg−1 , Peg = −Feg−1 , K Teg = 0, Feg = 0, 5

(4)

then [L, M ] = 0. e g−k , Peg−k , Teg−k , Feg−k from (3). This formulas are too hard for Let us calculate K analyzing but in the sequel only special cases are considered. We are going to show that in some cases formulas (4) give a very effective methods for finding commuting operators. Let us describe the main idea. We know that A0 = 0 because M is e 0 = 0 and Pe0 = 0. Then we can operator of order 2g. Let us take B0 such that K e e e e calculate T0 , F0 . Using (4) we can find K1 , P1 , then a12 , a13 , a14 , b11 , b12 , b13 , b14 . And using (4) we can find Te1 and Fe1 . So, we get recurrence relations ( m m m m m m m i = 1, 2, 3, 4 am+1 = gi (am 1 , a2 , a3 , a4 , b1 , b2 , b3 , b4 , r1 , r2 , r3 , q1 , q2 , q3 , q4 ), i m+1 m m m m m m m m bi = hi (a1 , a2 , a3 , a4 , b1 , b2 , b3 , b4 , r1 , r2 , r3 , q1 , q2 , q3 , q4 ) i = 1, 2, 3, 4. We will see that if there exists g such that ( ag+1 =0 i i = 1, ..., 4 g+1 bi = 0 then the operator L commutes with operator M . Let us suppose that λ1 = 1, λ2 = −1, λ3 = 0, r2 (x) = r3 (x) = 0, q3 = q2 and q4 = −q1 . Then we have     0 −2ag−k 2(a1g−k )0 −2bg−k − 2(a3g−k )0 3 3 e e Kg−k = , Pg−k = , 0 −2ag−k 0 −2b2g−k − 2(a2g−k )0 2(ag−k 2 4 ) ! ! 1 1 1 1 e e e e Fg−k,2 F Tg−k,2 T , , Feg−k = eg−k,1 Teg−k = eg−k,1 2 2 2 2 e Fg−k,1 Feg−k,2 Tg−k,1 Tg−k,2 g−k 0 g−k 0 g−k 0 g−k 00 1 Teg−k,1 = a2g−k q2 + ag−k 3 q2 − r1 (a1 ) + 2(b1 ) − a1 r1 + (a1 ) , g−k g−k 0 g−k 0 g−k 0 g−k 00 Te1 = −ag−k 1 q2 + a4 q2 − r1 (a3 ) + 2(b3 ) − a3 r + (a3 ) , g−k,2

2 Teg−k,1 = 2 Teg−k,2 =

1

a1g−k q2 ag−k 2 q2

− +

ag−k 4 q2 g−k a3 q 2

+ +

r1 (a2g−k )0 r1 (a4g−k )0

− −

0 2(bg−k 2 ) 0 2(bg−k 4 )

+ +

a2g−k r10 a4g−k r10

− (a2g−k )00 , − (a4g−k )00 ,

g−k g−k 0 g−k 0 g−k 0 00 0 g−k 1 + q20 ag−k + (bg−k Feg−k,1 = bg−k 3 1 ) , 2 q2 + b3 q2 − 2q1 (a1 ) + 2q2 (a3 ) + r1 (b1 ) − q1 a1 g−k g−k 0 g−k 0 g−k 0 1 00 0 g−k − q20 a1g−k + (bg−k Feg−k,2 = −bg−k 3 ) , 1 q2 + b4 q2 − 2q2 (a1 ) − 2q1 (a3 ) + r1 (b3 ) − q1 a3 g−k g−k 0 g−k 0 g−k 0 2 0 g−k 00 Feg−k,1 = bg−k − q20 a4g−k − (bg−k 1 q2 − b4 q2 + 2q1 (a2 ) − 2q2 (a4 ) − r1 (b2 ) + q1 a2 2 ) , g−k g−k 0 g−k 0 g−k 0 00 0 g−k 2 + q20 a2g−k − (bg−k Feg−k,2 = bg−k 4 ) 2 q2 + b3 q2 + 2q2 (a2 ) + 2q1 (a4 ) − r1 (b4 ) + q1 a4

6

From (4) we get e g−k+1 = −Teg−k , Peg−k+1 = −Feg−k . K So, we obtain recurrence relations a0i (x) ≡ 0, i = 1, ...4, 0 0 b2 (x) = b3 ≡ 0

(5)

1 =0⇔ Teg−k,1

bg−k 1

1 =− 2

Z 

a2g−k q2

+

ag−k 3 q2

−

0 (ag−k 1 ) r1

−

0 ag−k 1 r1

+

00 (ag−k 1 )



dx + C1g−k ,

(6)

2 Teg−k,2 =0⇔

bg−k = 4

Z 1  2

 g−k 0 g−k 0 g−k 00 g−k ag−k q + a q + (a ) r + a r − (a ) dx + C2g−k , 2 2 1 2 3 4 4 4 1

(7)

1 −2ag−k+1 = −Teg−k,2 ⇔ 3

ag−k+1 3

 1  g−k g−k g−k 0 g−k 0 g−k 0 g−k 00 −a1 q2 + a4 q2 − (a3 ) r1 + 2(b3 ) − a3 r1 + (a3 ) , = 2

(8)

2 −2ag−k+1 = −Teg−k,1 ⇔ 2

a2g−k+1

 1  g−k g−k g−k 0 g−k 0 g−k 0 g−k 00 a q2 − a4 q2 + (a2 ) r1 − 2(b2 ) + a2 r1 − (a2 ) , = 2 1

(9)

1 2(a1g−k+1 )0 = −Feg−k,1 ⇔

ag−k+1 =− 1

Z 1 

ag−k+1 =− 4

Z 1 

g−k g−k 0 g−k 0 g−k 0 bg−k 2 q2 + b3 q2 − 2(a1 ) q1 + 2q2 (a3 ) + (b1 ) r1 − 2  g−k 0 g−k 0 g−k 00 −a1 q1 + a3 q2 + (b1 ) dx + C3g−k+1 ,

(10)

2 ⇔ 2(a4g−k+1 )0 = −Feg−k,2

2

g−k 0 g−k 0 g−k g−k 0 bg−k 2 q2 + b3 q2 + 2(a2 ) q2 + 2(a4 ) q1 − (b4 ) r1 +  g−k 0 g−k 00 0 +ag−k q + a q − (b ) dx + C4g−k+1 , 4 2 4 1 2

7

(11)

1 −2b3g−k+1 − (a3g−k+1 )0 = −Feg−k,2 ⇔

bg−k+1 = 3

1  g−k g−k 0 g−k 0 g−k 0 −b1 q2 + bg−k 4 q2 − 2(a1 ) q2 − 2(a3 ) q1 + (b3 ) r1 − 2  g−k 00 g−k 0 0 )0 , ) − (ag−k+1 q + (b q − a −ag−k 3 3 1 3 2 1

(12)

2 ⇔ −2bg−k+1 − (ag−k+1 )0 = −Feg−k,1 2 2

bg−k+1 = 2

1  g−k g−k 0 g−k 0 0 b q2 − b4g−k q2 + 2(ag−k 2 ) q1 − 2(a4 ) q2 − (b2 ) r1 + 2 1  g−k 0 g−k 0 g−k 00 +a2 q1 − a4 q2 − (b2 ) − (a2g−k+1 )0 .

We see that

( Teg = 0 Feg = 0

⇔

(13)

( e g+1 = 0 K Peg+1 = 0

where  e g+1 = K

 0 −2ag+1 3 , 2ag+1 0 2

 Peg+1 =

0 2(ag+1 1 ) g+1 0 −2b2 − 2(ag+1 2 )

0 −2bg+1 − 2(ag+1 3 3 ) g+1 0 2(a4 ) ,



where Cij are arbitrary constants. Let us note that if a0i (x) = 0, b02 (x) = 0, b03 (x) = 0 for all i = 1, ..., 4, then from (6) and (7) we get that b01 (x) = const and b04 (x) = const. We obtain the following theorem Theorem 1. If there exists number g and constants of integration Cjm such that ag+1 = 0, bg+1 = 0, bg+1 = 0 for all i = 1, ..., 4, then the operator 2 3 i       1 0 r1 (x) 0 q1 (x) q2 (x) 2 L= ∂ + ∂ + 0 −1 x 0 −r1 (x) x q2 (x) −q1 (x) commutes with operator M = B0 Lg + (A1 ∂x + B1 )Lg−1 + ... + A0 ∂x + B0 ,   C31 r1 (x) (µ1 − µ2 )q2   1 − µ1 0 C1 +  2 21 B0 = , B1 =  , (µ1 − µ2 )q2 C r (x) 0 µ2 1 C14 + 4 2 2 8

(14)

where µ1 , µ2 are arbitrary constants and Cij are some constants. We see that if µ1 6= µ2 and q2 6= const, then B1 is not constant matrix hence M is not polynomial in L. Theorem 2. The operator      2  1 0 α 2 x2 + α 0 0 βx + α2 x γx 2 L= ∂ + ∂ + , 0 −1 x 0 −α2 x2 − α0 x γx −βx2 − α2 x where γ 2 = −n2 α22 ,

n∈N

and α2 , α0 , β are arbitrary constants, commutes with differential operator (14) of order 4n, where g = 2n. The order of operator M equals 4n. Remark. Calculations show that if n 6 3, then the spectral curve of operators L, M from Theorem 2 is nonsingular for almost all α0 , α2 , β and is hyperelliptic. Hence L and M are operators of rank 2. In some cases spectral curve is reducible and we get commuting operators of rank (2, 2). Note that operators from Theorem 2 can’t be operators of rank 1. Also note that from Theorem 1 we see that the matrix operator M from Theorem 2 is operator with polynomial coefficients. Example 1. If n = 1 and µ1 = 1, µ2 = −1, then the operator      2  1 0 α 2 x2 + α 0 0 βx + α2 x iα2 x 2 L= ∂ + ∂ + 0 −1 x 0 −α2 x2 − α0 x iα2 x −βx2 − α2 x commutes with operator M = B0 L2 + A1 ∂x L + B1 L + A2 ∂x + B2 . Calculations show that     1 0 2(α2 x2 + α0 ) 0 4 M= ∂ + ∂ 3+ 0 −1 x 0 −2(α2 x2 + α0 ) x   2 4 iα2 x α2 x + 2(α0 α2 + β)x2 + 6α2 x + α02 ∂ 2+ + iα2 x −α22 x4 − 2(α0 α2 + β)x2 − 6α2 x − α02 x       m1 m2 h1 h2 C0 0 + ∂ + + C1 L + , m2 −m1 x h2 −h1 + 2β − α0 α2 0 C0 m1 = 2α2 βx4 + 4α22 x3 + 2α0 βx2 + 4(α0 α2 + β)x + 4α2 , m2 = iα22 x3 + iα0 α2 x + iα2 , 3α2 h1 = β22 x4 + 4α2 βx3 + 2 x2 + 2α0 βx + 4β, 2 9

3 iα2 α0 h2 = iα2 βx3 + iα22 x2 + , 2 2 where C1 and C0 are arbitrary constants. The spectral curve of operators L, M has the form 2  α2 α0 − 2β ) = z 4 − (α0 α2 − 2β)z 2 − α2 α0 β + β 2 . w − C1 z − (C0 − 2 α2 α0 − 2β If we take C0 = , C1 = 0, then we get 2 w2 = z 4 − (α0 α2 − 2β)z 2 − α2 α0 β + β 2 . This spectral curve is nonsingular if α2 α0 β(α2 α0 − β) 6= 0. So in nonsingular case we get that operators L, M are operators of rank 2. If α0 = 0, then the spectral curve has the form w2 = (z 2 + β)2 ⇔ (w − z 2 − β)(w + z 2 + β) = 0 (15) We see that if α0 = 0, then the spectral curve is reducible. Note that M 6= L2 + β and M 6= −L2 − β but (M − L2 − β)(M + L2 + β) = 0 and we have operators of vector rank (2,2). Example 2. If n = 1 and µ1 = 1, µ2 = 2, then the operator      2  1 0 α 2 x2 + α 0 0 βx + α2 x iα2 x 2 L= ∂ + ∂ + 0 −1 x 0 −α2 x2 − α0 x iα2 x −βx2 − α2 x commutes with operator M = B0 L2 + A1 ∂x L + B1 L + A2 ∂x + B2 . Direct calculations show that     1 0 4 2(α2 x2 + α0 ) 0 M= ∂ + ∂ 3+ 0 2 x 0 4(α2 x2 + α0 ) x   iα2 x 2 4 2 2 −  2 α2 x + 2(α0 α2 + β)x + 6α2 x + α0 2 +  ∂x + iα2 x 2α22 x4 + 4(α0 α2 + β)x2 + 12α2 x + 2α02 − 2      m1 m3 h1 h3 C0 0 + ∂ + + C1 L + , m2 m4 x h2 h4 0 C0 m1 = 2α2 βx4 + 4α22 x3 + 2α0 βx2 + 4(α0 α2 + β)x + 4α2 , iα2 x3 iα0 α2 x 7iα2 m2 = − 2 − − , 2 2 2 iα2 x3 iα0 α2 x 5iα2 m3 = − 2 − + , 2 2 2 10

m4 = 4α2 βx4 + 8α22 x3 + 4α0 βx2 + 8(α0 α2 + β)x + 8α2 , 3α2 α0 3α2 , h1 = β 2 x4 + 4α2 βx3 + 2 x2 + 2α0 βx + β + 4 2 iα2 βx3 9iα22 x2 7iα2 α0 h2 = − − − , 2 4 4 iα2 βx3 3iα22 x2 5iα2 α0 + + , h3 = − 2 4 4 9α2 h4 = 2β 2 x4 + 8α2 βx3 + 2 x2 + 4α0 βx + 4β + 2α2 α0 , 4 where C1 and C0 are arbitrary constants. If we take C1 = 0 and C0 = 0, then the spectral curve of operators L, M has the form 16w2 − 8w(α2 α0 − 2β + 6z 2 ) + 32z 4 + 16(α2 α0 − 2β)z 2 + α22 α02 = 0 We see that the spectral curve is nonsingular for almost all α2 , α0 .β and L, M are operators of rank 2. If α0 = 0, then the spectral curve has the form (w − 2z 2 )(w − z 2 + β) = 0 and L, M are operators of vector rank (2, 2). Theorem 3. Let ℘(x) be the Weierstrass elliptic function satisfying the equation (℘0 (x))2 = 4℘3 (x) + g2 ℘(x). The operator     1 0 0 α℘(x) 2 L= ∂ + , 0 −1 x α℘(x) 0 α2 = 64n4 − 4n2 ,

n∈N

commutes with a differential operator (14), of order 4n, where g = 2n. The order of operator M equals 4n. Proofs of Theorem 2 and Theorem 3 C2k

We prove Theorem 2 and Theorem 3 using Theorem 1. Let us suppose that = C3k = C4k = 0 and C12k+1 = 0 for all k. We know that   µ1 0 A0 ≡ 0, B0 = . 0 µ2

11

Direct calculations using (6) − (13) show that a11 = a12 = a13 = a14 = 0, µ1 − µ2 µ1 − µ2 b12 = q2 , b13 = − q2 = −b12 . 2 2 Then

(16)

µ1 − µ2 0 q2 = −(b12 )0 , 2 2 2 0 r a − (a 1 2 2) b22 = b23 = 2 3 3 3 a1 = a2 = a3 = a34 = 0, b32 = −b33

a21 = a24 = 0,

a22 = a23 = −

(17)

Lemma 1. If k = 2m + 1, then ak1 = ak2 = ak3 = ak4 = 0,

bk2 = −bk3

(18)

If k = 2m, then ak1 = ak4 = 0,

ak2 = ak3 = −(b2k−1 )0 ,

bk3 = bk2 =

r1 ak2 − (ak2 )0 2

(19)

Proof We see that relations (18) and (19) is true when k = 1 and k = 2. Let us suppose that (18) and (19) is true for some k. If k = 2m + 1, then using (6) − (13) we get a2m+2 = a2m+2 = 0, 1 4

a2m+2 = a32m+2 = −(b2m+1 )0 , 2 2

b32m+2 = b2m+2 = 2

r1 a2m+2 − (a2m+2 )0 2 2 . 2

If k = 2m, then again using (6) − (13) we have a2m+1 = a2m+1 = a2m+1 = a2m+1 =0 1 2 3 4

b32m+1 = −b2m+1 . 2

The Lemma is proved. Proof of Theorem 2 From (16) and (17) we get a11 = a12 = a13 = a14 = 0, µ1 − µ2 µ1 − µ2 µ1 − µ2 µ1 − µ2 b12 = q2 = γx, b13 = − q2 = − γx. 2 2 2 2 12

Then µ1 − µ2 γ = −(b12 )0 , 2 µ1 − µ2 µ1 − µ2 b22 = b23 = − α0 γ − α2 γx2 , 4 4 a31 = a32 = a33 = a34 = 0, 2C12 + (µ1 − µ2 )(α0 α2 − 2β) µ1 − µ2 b32 = −b33 = γx + γ(γ 2 + α22 )x3 4 4 a21 = a24 = 0,

a22 = a23 = −

(20)

We want to prove that L commutes with differential operator (14), where g = 2n. From Theorem 1 and Lemma 1 we know that we must prove that there exists constants C12k such that b2n+1 ≡ 0. Let us note that recurrence relations (6) − (13) are 2 k+1 k+1 linear in ai and bi . Assume that b2m−1 = x2m−1 . Then we have i α0 (2m − 1)x2m−2 α2 (2m − 1)x2m − = −b2m 3 , 2 2 = −(b2m−1 )0 = −(2m − 1)x2m−2 = a2m 2 3 .

2m−3 b2m − 2 = (2m − 1)(m − 1)x

a2m 2

Again using (6) − (13) we obtain a2m+1 = a22m+1 = a32m+1 = a2m+1 = 0, 1 4 2 2 2 (2m − 1)(α2 m + γ ) 2m+1 b2m+1 = x + 2 2m (α2 α0 − 2β)(2m − 1)2 2m−1 α02 (2m − 2)(2m − 1) 2m−3 + x + x − 2 4 (2m − 1)(2m − 2)(2m − 3)(2m − 4) 2m−5 C12m β − x + x. 4 2

(21)

From (20) we see that b32 = K13 x+K33 x3 , where K13 = K13 (C12 ) is constant and depends µ1 − µ2 on C12 , K33 = γ(γ 2 + α22 ). Let us suppose that for some m 4 2m−1 2m−1 b2m−1 = K12m−1 x + K32m−1 x3 + ... + K2m−1 x , 2

where (µ1 − µ2 ) 2m−1 K2m−1 =−

m−1 Q

((2j − 1)(α2 j 2 + γ 2 ))

j=1

2m (m − 1)!

,

2m−1 2m−1 is constant and depends on C12 , K2m−5 is constant and depends on C12 , C14 , K2m−3 2j 2m−1 depends on C12 , ..., C1 and K12m−1 depends on C12 , C14 , ..., C12m−2 . We see K2m−2j−1

13

that it is true when m = 2. Using (21) we get 2m+1 2m−1 2m+1 2m+1 b2m+1 = K12m+1 x + K32m+1 x3 + ... + K2m−1 x + K2m+1 x , 2 2m 2 3α α2 α0 − 2β 2m−1 C β K12m+1 = 1 − 30K52m−1 + 0 K32m−1 + K1 2 2 2 α2 + γ 2 2m−1 9(α2 α0 − 2β) 2m−1 K32m+1 = 2 + + 5α02 K52m−1 − 210K72m−1 K1 K3 2 2 ...

(2m − 3)(α22 (m − 1)2 + γ 2 ) 2m−1 (α2 α0 − 2β)(2m − 1)2 2m−1 K2m−3 + K2m−1 2m − 2 2 (2m − 1)(α22 m2 + γ 2 ) 2m−1 2m+1 = K2m+1 K2m−1 . 2m 2m+1 2m−1 Easy to see that K2m−1 depends on constant of integration C12 because K2m−3 de2m+1 2m+1 2 2 4 pends on C1 . And K2m−3 depends on C1 , C3 . The last coefficient K1 depends on constants of integrations C12 , ..., C12m . Now let us consider 2m+1 K2m−1 =

2n+1 2n−1 2n+1 2n+1 b2n+1 = K12n+1 x + K32n+1 x3 + ... + K2n−1 x + K2n+1 x 2 2n+1 We know that γ 2 + n2 α2 = 0 and hence K2n+1 = 0. To prove Theorem 2 we must 2n+1 2n+1 2n+1 2n 2 = K3 = ... = K2n−1 = 0. It is always find constants C1 , ..., C1 such that K1 2n+1 2n+1 possible because K2n−1 depends n C12 and K2n−3 depends on C12 , C32 etc. The last coefficient K12n+1 depends on constants of integration C12 , ..., C12n .

Theorem 2 is proved. Proof of Theorem 3 The proof of Theorem 3 coincides with the proof of Theorem 2. Let us prove that principle parts of functions a2n+1 and bi2n+1 equals zero. We see from (6) − (13) that i aji and bji are elliptic functions for any i, j. Hence if principle parts of functions ai2n+1 and bi2n+1 equals zero, then a2n+1 and b2n+1 havn’t poles and hence are constants. In i i our case these constants are zeroes. From (16) and (17) we get a11 = a12 = a13 = a14 = 0, µ1 − µ2 µ1 − µ2 µ1 − µ2 q2 = α + O(x2 ), b13 = − α + O(x2 ). b12 = 2 2 2x 2x2 14

Then µ1 − µ2 α + O(x) = −(b12 )0 , x3 3(µ1 − µ2 ) (µ1 − µ2 )g2 α b22 = b23 = + O(x4 ), α+ 4 2x 40 a31 = a32 = a33 = a34 = 0, (µ1 − µ2 )α(α2 − 60) α(80C12 + 6(µ1 − µ2 )g2 α2 ) b32 = −b33 = + + O(x2 ) 4x6 160x2 a21 = a24 = 0,

a22 = a23 =

(22)

We mentioned before that recurrence relations (6) − (13) are linear in ak+1 and bk+1 . i i 1 2m−1 2m−1 2m−1 2m−1 2m−1 2m−1 = a2 = a3 = a4 = 0. Assume that b3 = −b2 = 4m−2 and a1 x Then we have 8m2 − 6m + 1 2m 2m , b2 = b3 = x4m 4m − 2 2m−1 0 a2m ) = 4m−1 . 2 = −(b2 x Again using (6) − (13) we obtain a2m+1 = a22m+1 = a32m+1 = a2m+1 = 0, 1 4 2 2 4 (2m − 1)(α + 4m − 64m ) const const C12m α b2m+1 = + + + ... + + O(x2 ) 2 2mx4m+2 x4m−2 x4m−6 2x2 (23) 3 3 K K From (22) we see that b32 = 66 + 22 + O(x2 ), where K23 = K23 (C12 ) is constant and x x 2 (µ − µ )α(α − 60) 1 2 depends on C12 , K63 = . Let us suppose that for some m 4 b2m−1 = 2

2m−1 2m−1 K4m−6 K4m−2 K22m−1 + + ... + + O(x2 ), x4m−2 x4m−6 x2

where (µ1 − µ2 ) 2m−1 K4m−2 =−

m−1 Q

((2j − 1)(α2 + 4j 2 − 64j 4 )))

j=1

2m (m − 1)!

,

2m−1 2m−1 K4m−6 is constant and depends on C12 , K4m−10 is constant and depends on C12 , C14 , 2m−1 K4m−4j−2 depends on C12 , ..., C12j and K22m−1 depends on C12 , C14 , ..., C12m−2 . We see that it is true when m = 2. Using (23) we get 2m+1 2m+1 2m+1 2m−1 x + K2m+1 x . b2m+1 = K12m+1 x + K32m+1 x3 + ... + K2m−1 2

15

2m+1 2m+1 Easy to see that K4m−2 depends on constant of integration C12 . And K4m−6 de2m+1 2 2 pends on C1 , C3 . The last coefficient K2 depends on constants of integrations C12 , ..., C12m . Now let us consider 2n+1 2n−1 2n+1 2n+1 b2n+1 = K12n+1 x + K32n+1 x3 + ... + K2n−1 x + K2n+1 x 2 4n+2 We know that α2 + 4n2 − 64n4 = 0 hence K2n+1 = 0. To prove Theorem 3 we must 2n+1 2n+1 2n+1 2n 2 = K6 = ... = K4n−2 = 0. It is always find constants C1 , ..., C1 such that K2 2n+1 2n+1 2 2 possible because K4n−2 depends on C1 and K4n−2 depends on C1 , C32 etc. The last coefficient K22n+1 depends on constants of integration C12 , ..., C12n .

Theorem 3 is proved.

References [1] Wallenberg, G.: Uber die Vertauschbarkeit homogener linearer Differentialausdrucke. Arch. Math. Phys. 4 (1903), 252–268. [2] Schur, J.: Uber vertauschbare lineare Differentialausdrucke. Sitzungsber. der Berliner Math. Gesell. 4 (1905), 2–8. [3] Burchnall J.-L., Chaundy T.W. Commutative ordinary differential operators. Proc. London Math. Soc. 21 (1923), 420-440; Proc. Royal Soc. London (A) 118 (1928), 557-583. [4] Grinevich, P. G.: Vector rank of commuting matrix differential operators. Proof of S. P. Novikov’s criterion, Math USSR IZV, 1987, 28 (3), 445–465. [5] Krichever, I. M.: Integration of nonlinear equations by the methods of algebraic geometry, Functional Analysis and Its Applications, 11: 1 (1977), 12–26. [6] Dixmier J.: Sur les algebres de Weyl. Bulletin de la Societe Mathematique de France 96, 209–242 (1968) [7] Krichever, I.M.: Commutative rings of ordinary linear differential operators, Functional Functional Analysis and Its Applications, 12:3 (1978), 175–185. [8] Krichever, I. M., Novikov, S.P.: Holomorphic bundles over algebraic curves and non-linear equations, Russian Mathematical Surveys, 1980, 35:6, 53–79.

16

[9] Mokhov, O. I.: Commuting ordinary differential operators of rank 3 corresponding to an elliptic curve, Russian Math. Surveys, 37:4 (1982), 129–130. [10] Mokhov, O. I.: Commuting differential operators of rank 3, and nonlinear differential equations, Math. USSR, Izvestiya,35:3 (1990), 629–655. [11] Mironov, A. E.: Self-adjoint commuting differential operators and commutative subalgebras of the Weyl algebra, Invent. math. (2014) 197:417-431. [12] Mironov, A. E.: Periodic and rapid decay rank two self-adjoint commuting differential operators, Amer. Math. Soc. Transl. Ser. 2, V. 234, 2014, P. 309– 322. [13] Oganesyan, V.: Commuting differential operators of rank 2 with polynomial coefficients, Functional Analysis and Its Applications, 2016, 50:1, 54–61. [14] Oganesyan, V.: Commuting differential operators of rank 2 and arbitrary genus g with polynomial coefficients, International Mathematics Research Notices (2016), doi:10.1093/imrn/rnw085. [15] Davletshina V.N.: Commuting differential operators of rank 2 with trigonometric coeffcients, Siberian Mathematical Journal, 56, 405-410 (2015) [16] Mironov, A. E., Zheglov, A. B.: Commuting Ordinary Differential Operators with Polynomial Coefficients and Automorphisms of the First Weyl Algebra, International Mathematics Research Notices 2015 : rnv218. [17] Mokhov, O. I.: Commuting ordinary differential operators of arbitrary genus and arbitrary rank with polynomial coefficients, American Mathematical Society Translations, Volume 234 (2014), 323-336. [18] Mokhov, O. I.: On Commutative Subalgebras of the Weyl Algebra Related to Commuting Operators of Arbitrary Rank and Genus, Mathematical Notes, vol. 94:2 (2013), 298-300. [19] Smirnov, A. O.: Finite-gap elliptic solutions of the KdV equation. Acta Appl. Math., 36 (1994), 125-166. [20] Smirnov, A. O.: Real finite-gap regular solutions of the Kaup-Boussinesq equation, Theoret. and Math. Phys., 66 (1986), 19-31.

17

[21] Dubrovin, B. A, Natanzon, S. M.: Real theta-function solutions of the Kadomtsev–Petviashvili equation, Mathematics of the USSR-Izvestiya, 1989, 32:2, 269–288. [22] Dubrovin, B. A, Natanzon, S. M.: Real two-zone solutions of the sine-Gordon equation, Functional Analysis and Its Applications, 1982, 16:1, 21–33. [23] Dubrovin, B. A.: The kadomcev-petviasvili equation and the relations between the periods of holomorphic differentials on riemann surfaces, Mathematics of the USSR-Izvestiya, 1982, 19:2, 285–296. [24] Shiota, T., Characterization of Jacobian varieties in terms of soliton equations. Invent. Math. 83 (1986), no. 2, 333–382. [25] A. Belov-Kanel, M. Kontsevich, The Jacobian conjecture is stably equivalent to the Dixmier conjecture, Moscow Mathematical Journal 7 (2): 209–218. [26] Y. Tsuchimoto, Endomorphisms of Weyl algebra and p-curvatures, Osaka J. Math. 42: 435–452. [27] Gesztesy, F., Weikard, R.: A characterization of all elliptic algebro-geometric solutions of the AKNS hierarchy, Acta Mathematica, September 1998, Volume 181, Issue 1, pp 63-108. [28] Weikard, R.: On commuting matrix differential operators, New York J. Math., vol. 8, 2002, pp. 9-30. [29] Dubrovin, B. A.: Completely integrable Hamiltonian systems associated with matrix operators and Abelian varieties, Functional Analysis and Its Applications, 1977, 11:4, 265–277. [30] Krichever, I.M.: Algebraic curves and commuting matricial differential operators, Functional Analysis and Its Applications, 1976, 10:2, 144–146. Department of Geometry and Topology, Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, Moscow, 119991 Russia. E-mail address: [email protected]

18