Arxiv:2104.06123V2 [Nlin.SI] 19 Jul 2021

arXiv:2104.06123v2 [nlin.SI] 19 Jul 2021 ehdo the or Kadomtsev–Petvi sol method and to (KdV) 28] 4, Vries [3, Korteweg–de Santini the and including Ablowitz t Fokas, the by to proposed addition first in models, 2D as form scalar in written be can ralternatively, or 11 sotnrfre oa h w-iesoa oalatti Toda two-dimensional the as to referred often is (1.1) n h ztec equation Tzitzeica the and nyteeutosascae ihteKcMoyalgebras Kac–Moody the with associated their equations to the regard only with form multi-component similar take also qain hmevs o ntne h nre fteCart the of entries the instance, perspe For various themselves. from equations integrable proven were equations these nwhich in htadi a that uenme fnniereutosta oss nc”alge “nice” possess that equations wav nonlinear periodic m of and the number waves huge e.g. solitary as see such etc., phenomena method, sc nonlinear algebro-geometric inverse the the as instance, well for as models, math nonlinear many of this, solutions by Motivated solutions. explicit construct ..[2 8 5 2 3,ahg ls f2 oatp equatio Toda-type 2D of class huge a 33], 32, 25, 18, [12, e.g. (0 the and component 0th the namely matrix, n ihio 1]i 99 sa nerbeetnino h fam the of extension integrable an as 1979, in [17] Mikhailov where n hsi sotnrfre oa h w-iesoa 2)T (2D) two-dimensional the as to referred often (1 is it thus and Here . + eaeitne osuyte2 oatp qain ihnth within equations Toda-type 2D the study to intended are We ntemdr hoyo nerbesses h oino int of notion the systems, integrable of theory modern the In mn nerbesses eytpclmdli h (2 the is model typical very a systems, integrable Among h DTd qain(.)i lsl eae oterepresen the to related closely is (1.1) equation Toda 2D The e od n phrases. and words Key )dmninlsse opsdo nifiienme fcomp of number infinite an of composed system 1)-dimensional θ n ff a ar n h don a ar n iutnosyconst simultaneously and pairs Lax adjoint the and applicat pairs particular Lax As system. Toda two-dimensional the of fasto ieritga qain nacranfr.Sc a Such linea form. direct the certain induces a also in but equations structure, integral integrability linear of set a of i algebras Lie A ∂ ϕ erential 1 bstract n ϕ and = n − ϕ ∂ ¯ 1 ( h ietlnaiainfaeoki rsne o h tw the for presented is framework linearisation direct The . ∂ ehd nadto osliganniereuto,teD a DL the equation, nonlinear a solving to addition In method. x − / di − IERITGA QAIN N W-IESOA OASYSTEM TODA TWO-DIMENSIONAL AND EQUATIONS INTEGRAL LINEAR 1 1 , ϕ A ff x ∞ n tn o h ata eiaie ihrsetto respect with derivatives partial the for stand rneeuto seatyslal ne ninitial-bound an under solvable exactly is equation erence − n h oe eoetersetv oain ftecentra the of locations respective the denote boxes the and , , 1 ieritga qain nnt arx two-dimensiona matrix, infinite equation, integral linear B , ∞ n steptnilfnto ftocniuu ievariables time continuous two of function potential the is ) and C ∂ ∞ 1 ∂ swl steKcMoyalgebras Kac–Moody the as well as , − 1                            θ θ θ − ...... 0 1 1                            , = 0)-entry. − ∂                           . 1 ∂ . ∂ iigtp ouinptnilylaigt h nvra s universal the to leading potentially solution type rising . 1 2 ∂ ∂ ∂ − ϕ U I N E FU WEI AND YIN YUE 1 1 1 1 n ∂ ∂ ∂ ϕ − .I 1. . oso hsfaeokt h w-iesoa oalattices Toda two-dimensional the to framework this of ions − − − n = . nmti of matrix an 2 1 1 1 1 uttegnrlsdCuh arxsolutions. matrix Cauchy generalised the ruct . ceentol rvdsauie esetv ounderstand to perspective unified a provides only not scheme = ϕ ϕ ϕ teigtasom h abu rnfr,Hrt’ bilin Hirota’s transform, Darboux the transform, attering e(DL of (2DTL) ce e mtclmtoswr netdaddvlpdt erhfre for search to developed and invented were methods ematical tvs h hrceitco i ler lorflcsi th in reflects also algebra Lie of characteristic The ctives. epcieCra arcs ntewoecasfiaino t of classification whole the In matrices. Cartan respective ntroduction 0 0 0 + ϕ ri n emti properties. geometric and braic e A n − . ngah 1 1 6 6.Teetcnqe o nyexplain only not techniques These 26]. 16, 11, [1, onographs u n-iesoa 1)Td atc c.[31]) (cf. lattice Toda (1D) one-dimensional ous eiiilvlepolm fnnierprildi partial nonlinear of problems value initial ve = = = )dmninlequation 1)-dimensional ϕ − . swr bandbsdo h lsicto fLealgebras. Lie of classification the on based obtained were ns 1 (1) si auemteaial,btas oiae h discove the motivated also but mathematically, nature in es redmninl(D qain(1.1). equation (3D) hree-dimensional 2 1 n A svl K)eutos sakn fgnrlsto ftedr the of generalisation of kind a as equations, (KP) ashvili ϕ . -iesoa oaeutosascae ihteinfinite- the with associated equations Toda o-dimensional − 1 e d qain qain(.)cnatraieyb rte a written be alternatively can (1.1) Equation equation. oda n e e r (1) − ϕ 2 2 ϕ and 1 n ainter.I eiso aesfo aiu uhr,se authors, various from papers of series a In theory. tation ϕ ϕ gaiiyo olna qain fe eest h prope the to refers often equations nonlinear of egrability , 0 − − − − 0 0 1 A − oaltie a ar acymti solution. matrix Cauchy pair, Lax lattice, Toda l nnsa olw seeg [12]): e.g. (see follows as onents ocle ietlnaiain(L cee h Lwas DL The scheme. (DL) linearisation direct so-called e 1 2 1 − 2 (2) − − e r A e x A ϕ , e e e − 1 2 (2) n C ∞ − − + ϕ 2 − . 1 r (1) 2 ϕ ϕ n and . aeytesn–odnequation sinh–Gordon the namely , A − + ϕ per stecoe the as appears 1 0 0 . 1 ϕ 0 , , ∞ − and n ϕ , tp.Te2 oatp qain nohrclasses other in equations Toda-type 2D The -type. n − x . , . − D 1 . poc losmlaeul rvdsa insight an provides simultaneously also pproach 1 r odto,wihi aycssalw sto us allows cases many in which condition, ary r (2) + epciey hseuto a rpsdby proposed was equation This respectively. , 1 lmnso nifiievco n ninfinite an and vector infinite an of elements l . o rirr integers arbitrary for . . x                           1                            and e e e − − − ffi ...... θ θ θ x − 1 0 1 − insi 12,adhne equation hence, and (1.2), in cients 1                            n n iceesailvariable spatial discrete one and lto pc,frec class each for space, olution , r ∈ Z erdsoe the rediscover we , + rmteaspect the from , h underlying the S ff dimensional rnilequations erential a method, ear nonlinear e e2DTLs, he yo a of ry xplicit essing dthe ed (1.4b) (1.4a) (1.3) (1.1) (1.2) All a s rty e 2 YUE YIN AND WEI FU into understanding the underlying structures of integrable equations. This is reflected in many aspects including searching for integrable discretisation of differential equations (see e.g. the review papers [20,27]), constructing integrability characteristics of nonlinear equations such as Lax pairs [7,22,36], commuting symmetries [9], and master symmetries [21,23], etc. Recently, one of the authors and his collaborator further developed the DL and established the connection between linear integral equations and the Lie algebras [8–10], which extended the range of application of the DL. The main idea behind the DL approach, according to Fokas and Ablowitz [3], is to construct solutions of a nonlinear equation, by solving its relevant linear integral equation. This brings the benefit that no additional conditions are imposed at the beginning, and thus, a large class of solutions are obtainable as various degenerations from a so-called direct linearising solution, which takes the form of an integral defined on an arbitrary integration domain through an arbitrary integration measure in terms of the spectral variable(s). In this sense, the establishment of the DL scheme of the 2DTLs potentially provides a path towards their respective ‘universal’ solution spaces. In the present paper, we consider the 2D Toda-type equations associated with the infinite-dimensional algebras A∞, B∞, C∞, (1) (2) (1) (2) Z+ as well as the Kac–Moody Lie algebras Ar , A2r , Cr and Dr+1, for arbitrary integers r ∈ . Our starting point is a set of linear integral equations (namely a linear integral equation and its adjoint) taking the form of (2.5). By selecting suitable plane wave factors and kernel, i.e. (3.1) and (3.2), the A∞-type 2D Toda system is recovered as the nonlinearisation of the wave function of either of the linear integral equations, in the language of infinite matrix. From such a set of linear integral equations, we are able to construct its direct linearising type solution, and also the integrability characteristics including the Lax pair, the adjoint Lax pair as well as the tau function, for the 2DTL of A∞-type. By performing various constraints on the integration measure in terms (1) (2) (1) (2) of the spectral parameters, we obtain the respective DL schemes for the 2DTLs of B∞-, C∞-, Ar -, A2r -, Cr - and Dr+1-types and subsequently the corresponding direct linearising solutions, the linear problems, and also the tau functions. These form a unified framework for the integrability of the 2DTLs, which generalises the DL approach proposed by Fokas and Ablowitz to the 2D Toda systems of various types, and simultaneously provides an insight into the results obtained in the literature (cf. e.g. [6,25,32]) from a different perspective. In addition, we also discuss the so-called generalised Cauchy matrix solutions for the various 2D Toda-type equations. This is realised by taking a particular integration measure possessing an arbitrary fixed number of higher-order singularities, with respect to the spectral variable(s). As a consequence, the linear integral equation reduces to a linear algebraic equation, which allows us to establish the general formula of the Cauchy matrix type solution, for each 2DTL. This generalises the results obtained (1) in [19,34,35] from A1 -type to other types. The paper is organised as follows. The formal structure of our approach is introduced in section 2, illustrated in the language of infinite matrix. In section 3, we establish DL scheme of the 2DTL of A∞-type, from which we construct its direct linearising solution, tau function, and associated linear problems. The DL schemes and the relevant results of the 2DTLs of B∞-, C∞-, as (1) (2) (1) (2) well as Ar -, A2r -, Cr - and Dr+1-types are discussed in sections 4 and 5, respectively. In section 6, general formulae of the Cauchy matrix type solutions for the 2DTLs are presented.

2. Preliminaries 2.1. Inﬁnite matrices and vectors. Suppose that an inﬁnite matrix of size ∞×∞ is given by

......   ··· U−1,−1 U−1,0 U−1,1 ···   U = ··· U0,−1 U0,0 U0,1 ··· .   ··· U U U ···  1,−1 1,0 1,1   . . .   . . .      We define the operations for these objects by following the same rules in the finite-dimensional case, namely for arbitrary infinite matrices U = (Ui, j)∞×∞ and V = (Vi, j)∞×∞, the addition U +V , the scalar multiplication p U and the matrix multiplication UV give rise to infinite matrices of size ∞×∞, and their respective (i, j)-entries of the addition U + V , the scalar multiplication p U 1 and the matrix multiplication UV are defined as Ui, j + Vi, j, pUi, j and i′∈Z Ui,i′ Vi′, j . Similarly, an infinite column vector and its transpose are given by P

t t u = (··· , u−1, u0 , u1, ··· ) and u = (··· , u−1, u0 , u1, ··· ), respectively. Here we note that we adopt the notation t( · ) for the transpose of a matrix or a vector. For arbitrary infinite column vectors u = (ui)∞×1 and v = (vi)∞×1, u + v and p u are infinite column vectors whose ith components are defined as ui + vi t t and pui, respectively, u v is an ∞×∞ infinite matrix having its (i, j)-entry uiv j, and uv is a scalar defined by i∈Z uivi. The multiplication between an infinite matrix and an infinite vector is defined in the same way, i.e. Uv is an infinite columnP vector t with its ith component given by j∈Z Ui, jv j, and vU is an infinite row vector with its jth component defined as i∈Z viUi, j. P P

1The multiplication is a formal definition. In general, the infinite summation could raise the issue of divergence; however, we only deal with convergent infinite summation throughout the paper. LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 3

Next, we introduce some special inﬁnite matrices and vectors. The projection matrix is deﬁned as

...... ··· 0 0 0 ···   O ··· 0 1 0 ··· .   ··· 0 0 0 ···    . . .   . . .      It is easily veriﬁed that the projection matrix O satisﬁes the properties

(i, j) (i, j) 1, i = j, (OU) = δi,0U0, j and (UO) = Ui,0δ0, j, where δi, j = ∀i, j ∈ Z, ( 0, i , j, as well as O2 = O and tO = O. Here the notation ( · )(i, j) stands for the (i, j)-entry of an inﬁnite matrix. The index-raising matrix Λ and its transpose tΛ are deﬁned as

......  0 1  .     ..     0  Λ  0 1  and tΛ   ,    1 0   ..     0 .  1 0       ..   .. ..   .  . .         respectively. Multiplying U by Λ from the left and tΛ from the right, respectively, one obtains

(i, j) t (i, j) (Λ U) = Ui+1, j and (U Λ) = Ui, j+1, namely, the operation of the index raising matrices Λ and tΛ raise the row and column indices by 1, respectively. In our approach we also need the column and row inﬁnite vectors

c(k) t(··· , k−1, 1 , k, ··· ) and tc(k′) (··· , k′−1, 1 , k′, ··· ), (2.1) composed of monomials of k and k′. The operation of Λ and tΛ on these two inﬁnite vectors give rise to relations

Λc(k) = k c(k) and tc(k′)tΛ = k′tc(k′), (2.2) or equivalently, [Λc(k)](i) = ki+1 and [tc(k′)tΛ](i) = k′i+1, where the notation ( · )(i) denotes the ith component of an inﬁnite vector. In addition, it is also proven that

tc(k′)tΛ jOΛic(k) = kik′ j.

The trace of an arbitrary inﬁnite matrix U is deﬁned as

tr U = Ui,i, Xi∈Z and it possesses the property

tr(UV ) = tr(V U). (2.3)

The determinant of an inﬁnite matrix is deﬁned through

ln(det U) = tr(ln U).

In practice we only deal with the determinant of a very special inﬁnite matrix in the form of 1 + ∗ (where 1 denotes the identity inﬁnite matrix) in our paper, and thus, we are allowed to reformulate the determinant as

∞ (−1)i−1 det(1 + ∗) = exp{tr[ln(1 + ∗)]} = exp tr ∗i ,  i  Xi=1    which is meaningful as long as the inﬁnite matrix ∗ guarantees convergent expansion. The determinant satisﬁes the property

det(1 + UOV ) = 1 + (V U)(0,0), (2.4) namely the well-known Weinstein–Aronszajn formula for the rank 1 case. 4 YUE YIN AND WEI FU

2.2. Infinite matrix representation of the linear integral equations. We start with introducing the notations for the key objects in the linear integral equations. The infinite column vector u(k) and the infinite row vector tv(k′) are the wave functions, whose components are functions of the dynamical variables (which could be either continuous or discrete arguments), relying on the spectral parameters k and k′, respectively. For the plane wave factors, we adopt the notations ρ(k) and σ(k′), which are scalars dependent on the dynamical variables, as well as their respective spectral parameters k and k′. The Cauchy kernel Ω(k, k′)isa scalar function of only the spectral parameters k and k′. The linear integral equations we consider in our scheme are as follows:

u(k) + dζ(l, l′)ρ(k)Ω(k, l′)σ(l′)u(l) = ρ(k)c(k), (2.5a) "D tv(k′) + dζ(l, l′)ρ(l)Ω(l, k′)σ(k′)tv(l′) = tc(k′)σ(k′), (2.5b) "D where k and k′ (or l and l′) are the spectral parameters in the complex field C, dζ is the integration measure, and D is the integration domain. In the paper, we refer to equation (2.5b) as the ‘adjoint’ of (2.5a). This set of equations, i.e. (2.5), is the starting point for the DL approach, and at this stage, there is no restriction (such as restricted domain and measure, etc.) on (2.5), except that we require that the homogeneous equation corresponding to any of (2.5a) and (2.5b) has only zero solution. Next, we consider the infinite matrix representation of the linear integral equations. If we introduce the infinite matrices Ω and U, which are defined by Ω(k, k′) tc(k′)Ωc(k) (2.6) and

U dζ(k, k′)u(k)tc(k′)σ(k′), or equivalently U dζ(k, k′)ρ(k)c(k)tv(k′), (2.7) "D "D respectively, and replace the kernel Ω(k, k′) with the help of (2.6) and (2.7), the linear integral equations in (2.5) are reformulated as follows: u(k) = (1 − UΩ)c(k)ρ(k), (2.8a) tv(k′) = σ(k′)tc(k′)(1 − ΩU). (2.8b)

Notice the first definition of U given in (2.7). We immediately derive from (2.8a) the infinite matrix relation2 U = (1 − UΩ)C, or alternatively U = C(1 + ΩC)−1, (2.9) where the infinite matrix C is defined as

C dζ(k, k′)ρ(k)c(k)tc(k′)σ(k′). (2.10) "D In our framework, we also define the formal tau function as τ det(1 + ΩC), (2.11) or equivalently τ det(1 + CΩ). The infinite matrix U, the infinite vectors u(k) and tv(k′) and the tau function τ are the key quantities in the our framework. Once the kernel and the plane wave factors in the linear integral equations are fully determined, we are able to construct nonlinear equations and its integrability characteristics based on these quantities. More concretely, the entries of the infinite matrix U (or their combinations) form potentials of the resulting nonlinear integrable equations; the components of the infinite vectors u(k) and tv(k′) play roles of the eigenfunctions of the associated linear problems (i.e. Lax pairs and the adjoint ones); while the tau function is the potential for the corresponding bilinear equations. To put it another way, (2.7) subject to (2.5) provides us with a general solution to the resulting nonlinear integrable equations within our scheme, which contains a large class of exact solutions as special cases.

3. 2DTL associated with A∞

3.1. Closed-form nonlinear equations. We start with the 2DTL of A∞-type in this section from the formal structure of the DL. To establish the scheme, we select the plane wave factors and the Cauchy kernel given by

−1 ′ ′−1 kx1+k x−1 n ′ k x1+k x−1 ′ −n ρn(k) = e k and σn(k ) = e (−k ) (3.1) as well as 1 Ω(k, k′) = , (3.2) k + k′ respectively, where x1 and x−1 are the continuous temporal variables, n is the discrete spatial variable. Since we need to consider the discrete dynamical evolutions in our scheme, from now on we adopt the suﬃx n for various quantities, in order to explicitly illustrate how these quantities evolves with respect to n, when it is necessary.

2Similarly, one is able to derive U = C(1 − ΩU) from (2.8b) by making use of the second deﬁnition in (2.7). This leads to another relation U = (1 + CΩ)−1C, which is equivalent to (2.9). This shows why the two deﬁnitions in (2.7) are equivalent. LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 5

Notice the deﬁnition of Cn, namely (2.10), for given plane wave factors in (3.1). We can prove that Cn satisﬁes the following dynamical evolutions: t ∂1Cn = ΛCn + Cn Λ, (3.3a) −1 t −1 ∂−1Cn = Λ Cn + Cn Λ , (3.3b) t Cn+1(− Λ) = ΛCn. (3.3c)

The derivation of these equations is direct. For example, taking the derivative of Cn with respect to x−1 provides us with

′ ′ t ′ ∂−1Cn = dζ(k, k )c(k)∂−1[ρn(k)σn(k )] c(k ) "D ′ ′ t ′ ′ ′ t ′ ′ = dζ(k, k )kc(k)ρn(k)σn(k ) c(k ) + dζ(k, k )c(k)ρn(k)σn(k ) c(k )k . "D "D And then by making use of the property (2.2), we end up with (3.3b). Meanwhile, the given Cauchy kernel (3.2) implies that Ω ∞ tΛ −i−1 Λi = − i=0(− ) O , since ∞ ∞ P 1 tc(k′) − (−tΛ)−i−1OΛi c(k) = − (−k′)−i−1ki = ,   k + k′  Xi=0  Xi=0   according to (2.6). This in turn tells us that Ω satisﬁes the algebraic relation ΩΛ + tΛΩ = O. (3.4) From equations (3.3) and (3.4), we are able to prove that the inﬁnite matrix U obeys the dynamical evolutions as follows: t ∂1Un = ΛUn + Un Λ − UnOUn, (3.5a) −1 t −1 t −1 −1 ∂−1Un = Λ Un + Un Λ − Un Λ OΛ Un, (3.5b) t Un+1(− Λ) = ΛUn − Un+1OUn. (3.5c)

We take (3.5b) as an example to show the derivation. By diﬀerentiating equation (2.9) with respect to x−1, we have

∂−1Un = (1 − UnΩ)(∂−1Cn) − (∂−1Un)ΩCn, which is equivalent to −1 t −1 −1 t −1 −1 (∂−1Un)(1 + ΩCn) = (1 − UnΩ)(Λ Cn + Cn Λ ) = Λ Cn + Un Λ − UnΩΛ Cn. Notice that (3.4) is equivalent to ΩΛ−1 + tΛ−1Ω = tΛ−1OΛ−1. The above equation is reformulated as −1 t −1 t −1 −1 t −1 (∂−1Un)(1 + ΩCn) = Λ Cn + Un Λ − Un( Λ OΛ − Λ Ω)Cn −1 t −1 t −1 −1 = Λ Cn + Un Λ (1 + ΩCn) − Un Λ OΛ Cn. −1 Multiplying this equation by (1 + ΩCn) from the right, we end up with (3.5b) in virtue of (2.9). Equations (3.5a) and (3.5c) are proven similarly. Equations in (3.5), from our view point, form the infinite matrix representation of the 2DTL of A∞-type. Now we present how the scalar closed-form 2DTL (1.1) arises from (3.5). Taking the derivative of ln τn gives rise to −1 ∂1 ln τn = ∂1 ln[det(1 + ΩCn)] = ∂1 tr[ln(1 + ΩCn)] = tr[∂1 ln(1 + ΩCn)] = tr[(1 + ΩCn) Ω(∂1Cn)] −1 t −1 −1 t = tr[(1 + ΩCn) Ω(ΛCn + Cn Λ)] = tr[(1 + ΩCn) ΩΛCn + (1 + ΩCn) ΩCn Λ] −1 −1t = tr[Cn(1 + ΩCn) ΩΛ + Cn(1 + ΩCn) ΛΩ] = tr(UnO) = tr(OUn), where the properties (3.4) and (2.3) are used. Therefore, we have (0,0) (−1,−1) ∂1 ln τn = Un , and similarly, ∂−1 ln τn = Un . (3.6a) Simultaneously, by performing the shift operation on the tau function with respect to the discrete variable n, we obtain t −1 t t −1 t −1 τn+1 = det(1 + ΩCn+1) = det 1 + ΩΛCn(− Λ) = det 1 − (O − ΛΩ)Cn Λ = det 1 − Λ OCn + ΩCn Ω Ω −1tΛ−1 tΛ−1 (0,0) (0,−1) = det(1 + Cn)det 1 − (1 + Cn) OCn = τn 1 − (Un ) = τn 1 − Un , in which we have made use of (3.4) and (2.4); in other words, the tau function satisfies τn+1 (0,−1) τn−1 (−1,0) = 1 − Un , and similarly, = 1 − Un . (3.6b) τn τn Equations (3.6a) and (3.6b) transfer the dynamical evolutions of the tau function with respect to the independent arguments to the entries of the infinite matrix Un. Notice that (0,0) ∂1∂−1 ln τn = ∂−1(∂1 ln τn) = ∂−1Un , because of (3.6a). By taking the (0, 0)-entry of (3.5b), we obtain (0,−1) (−1,0) (0,−1) (−1,0) (0,−1) (−1,0) ∂1∂−1 ln τn = Un + Un − Un Un = 1 − 1 − Un 1 − Un . Therefore, a closed-form equation τn+1 τn−1 ∂1∂−1 ln τn = 1 − τn τn 6 YUE YIN AND WEI FU

(−1,0) (0,−1) is derived, if we replace Un and Un by the tau function according to (3.6b), which can alternatively be written as 1 D D τ · τ = τ2 − τ τ (3.7) 2 1 −1 n n n n+1 n−1 2 with the help of the bilinear identity 2∂1∂−1 ln τn = D1D−1τn · τn/τn, where the notations D1 and D−1 stand for Hirota’s bilinear 3 derivatives (see e.g. [11]) with respect to the corresponding arguments x1 and x−1, respectively. Equation (3.7) is the so-called (0,−1) bilinear 2DTL, see e.g. [11] and also [12,32]. If we introduce a new quantity ϕn ln(τn+1/τn) = ln 1 − Un , equation (3.7) turns out to be its nonlinear form, which is nothing but the nonlinear 2DTL of A∞-type, i.e. (1.1). We can further write down its Cartan matrix form based on the quantity θn ϕn−1 − ϕn, and subsequently obtain (1.2). 3.2. Lax pair and its adjoint. The linear equations in (2.5) not only play a key role in constructing exact solutions of the 2D Toda equation, but also help to derive Lax pairs of the 2DTL. We ﬁrst concentrate on equation (2.8a), i.e. the linear integral equation of un(k). Parallel to the equations in (3.5), the wave function un(k) satisﬁes the following dynamical relations:

∂1un(k) = Λun(k) − UnOun(k), (3.8a) −1 t −1 −1 ∂−1un(k) = Λ un(k) − Un Λ OΛ un(k), (3.8b)

un+1(k) = Λun(k) − Un+1Oun(k). (3.8c)

The derivation of these equations is straightforward. For instance, taking the derivative of (2.8a) with respect to x−1, we obtain

∂−1un(k) = − (∂−1Un)Ωc(k)ρ(k) + (1 − UnΩ)c(k)[∂−1ρn(k)] −1 t −1 t −1 −1 −1 = − (Λ Un + Un Λ − Un Λ OΛ Un)Ωc(k)ρ(k) + (1 − UnΩ)Λ c(k)ρn(k) −1 −1 t −1 t −1 −1 = Λ (1 − UnΩ)c(k)ρn(k) − Un(ΩΛ + Λ Ω)c(k)ρn(k) + Un Λ OΛ UnΩc(k)ρn(k). Notice that ΩΛ−1 + tΛ−1Ω = tΛ−1OΛ−1. This equation is rewritten as −1 t −1 −1 ∂−1un(k) = Λ (1 − UnΩ)c(k)ρn(k) − Un Λ OΛ (1 − UnΩ)c(k)ρn(k), which is nothing but (3.8b), according to equation (2.8a). Equations listed in (3.8) form the inﬁnite vector representation of a Lax pair of the 2DTL. To derive a closed-form scalar Lax pair, we need to eliminate Λun(k) in (3.8), and as a consequence, the following relations are obtained:

∂1un(k) = (Un+1 − Un)Oun(k) + un+1(k), −1 t −1 t −1 −1 ∂−1un(k) = un−1(k) − Λ UnOun−1(k) − Un Λ Oun−1(k) + Un Λ OΛ UnOun−1(k). The respective 0th components of the above two equations provide us with (0,0) (0,0) (0,−1) (−1,0) ∂1φn = Un+1 − Un φn + φn+1 and ∂−1φn = 1 − Un 1 − Un φn−1, (0) where φn [un(k)] . We therefore end up with τn+1τn−1 ∂ φ = ∂ τ − ∂ τ φ + φ = ∂ ϕ φ + φ ∂ φ = φ = ϕn−ϕn−1 φ , 1 n ( 1 ln n+1 1 ln n) n n+1 ( 1 n) n n+1 and −1 n 2 n−1 e n−1 τn t by making use of (3.6). By introducing the inﬁnite vector Φ = (··· ,φ−1, φ0 ,φ1 ··· ), we can alternatively express these two linear equations as a multi-component system composed of

∂1Φ= P Φ and ∂−1Φ= QΦ, (3.9) in which ......   −θ ∂1ϕ−1 1  e −1 0         −θ  P =  ∂1ϕ0 1  and Q =  e 0 0  .        −θ1   ∂1ϕ1 1   e 0       .. ..   .. ..   . .  . .         We refer to (3.9) as the Lax pair of the 2DTL of A∞-type (see [5,6,32]), since equation (1.1) arises as the compatibility condition of the equations in (3.9), i.e. ∂1∂−1Φn = ∂−1∂1Φn. Another Lax pair of (1.1) in our scheme is derived from the adjoint linear t ′ integral equation (2.5b). Performing a similar calculation on (2.8b), we can show that the adjoint wave function vn(k ) satisﬁes the following dynamical relations: t ′ t ′ t t ′ ∂1 vn(k ) = vn(k ) Λ − vn(k )OUn, (3.10a) t ′ t ′ t −1 t ′ t −1 −1 ∂−1 vn(k ) = vn(k ) Λ − vn(k ) Λ OΛ Un, (3.10b) t ′ t ′ t t ′ vn−1(k ) = − vn(k ) Λ + vn(k )OUn−1. (3.10c)

3 For two arbitrary diﬀerentiable functions f (x1, x−1) and g(x1, x−1), Hirota’s bilinear derivative D is deﬁned through i j i j ∂ ∂ ∂ ∂ ′ ′ D D f · g = − − f (x1, x−1)g(x , x− ) ′ ′ , 1 −1  ′   ′  1 1 x =x1,x =x−1 ∂x1 ∂x ∂x−1 ∂x 1 −1  1   −1  where i and j are positive integers.     LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 7

t ′ t −1 Following the same idea of deriving (3.9), we first eliminate vn(k ) Λ in (3.10). And then we can find that the quantity t ′ (0) ψn [ vn(k )] satisfies

τn+1τn−1 ∂ ψ = ∂ τ − ∂ τ ψ − ψ = − ∂ ϕ ψ − ψ ∂ ψ = − ψ = − ϕn−ϕn−1 ψ . 1 n ( 1 ln n−1 1 ln n) n n−1 ( 1 n−1) n n−1 and −1 n 2 n+1 e n+1 τn We can reformulate these two equations as

′ ′ ∂1Ψ= P Ψ and ∂−1Ψ= Q Ψ, (3.11)

t with Ψ= (··· , ψ−1, ψ0 , ψ−1 ··· ), where

......  − −   −θ−1   1 ∂1ϕ−2   0 −e  ′   ′      −θ0  P =  −1 −∂1ϕ−1  and Q =  0 −e  .      − −   0 −e−θ1   1 ∂1ϕ0       . .   .. ..   .. ..  . .           The linear problem (3.11) is referred to as the adjoint Lax pair of the 2DTL (cf. [25]), as the 2DTL equation (1.1) also arises as the compatibility condition ∂1∂−1Ψn = ∂−1∂1Ψn of (3.11). We conclude the main results of this section as theorems below.

Theorem 3.1. Consider a linear integral equation of the form

′ ′ ′ un(k) + dζ(l, l )ρn(k)Ω(k, l )σn(l )un(l) = ρn(k)c(k) (3.12) "D

′ ′ ′ with arbitrary measure dζ(k, k ) and integration domain D, where ρn(k), σn(k ) and Ω(k, k ) are deﬁned as (3.1) and (3.2), respectively. For any solution un(k) of (3.12), the quantity

(0,−1) ϕn ln 1 − Un (3.13) provides the direct linearising type solution to the nonlinear 2D Toda equation (1.1), where

′ t ′ ′ Un dζ(k, k )un(k) c(k )σn(k ). (3.14) "D

(0) Meanwhile, the wave function φn [un(k)] satisﬁes the linear problem (3.9).

Theorem 3.2. Consider the adjoint linear integral equation

t ′ ′ ′ ′ t ′ t ′ ′ vn(k ) + dζ(l, l )ρn(l)Ω(l, k )σn(k ) vn(l ) = c(k )σn(k ). (3.15) "D

t ′ For any solution vn(k ) of (3.15), the quantity ϕn deﬁned as (3.13) also provides a solution to the nonlinear equation (1.1), where Un is deﬁned by

′ t ′ Un dζ(k, k )ρn(k)c(k)vn(k ). "D

t ′ (0) Furthermore, the adjoint wave function ψn [ vn(k )] provides a solution to the adjoint linear problem (3.11).

Theorem 3.3. The formal tau function

τn det(1 + ΩCn) (3.16)

Ω − ∞ −tΛ −i−1OΛi C ′ c tc ′ ′ with i=0( ) and n D dζ(k, k )ρn(k) (k) (k )σn(k ) provides the direct linearising type solution to the bilinear 2D TodaP equation (3.7). ! We remark that the formal tau function is fully determined by the key ingredients of the linear integral equations (3.12) and (3.15), including the plane wave factors, the kernel as well as the integration measure, through Ω and Cn, which provides the most general solution within the DL scheme. The theorems illustrate the DL scheme for the 2D Toda equation of A∞-type, namely the most general case. In the coming sections 4 and 5, by performing various reductions on the integration measure dζ(k, k′) and the integration domain D, 2DTLs (1) (2) (1) (2) associated with B∞ and C∞, as well as Ar , A2r , Cr and Dr+1 are recovered as special reductions, from the view point of the linear integral equations. 8 YUE YIN AND WEI FU

4. 2DTLs associated with B∞ and C∞

4.1. B∞-type. To recover the B∞-type 2DTL, we identify that the measure in the linear integral equations (3.12) and (3.15) satisﬁes dζ(k, k′) = dζ′(k, k′)k, inwhich dζ′(k, k′) = −dζ′(k′, k), ∀(k, k′) ∈ D, (4.1) with symmetric D. Such a restriction results in the fact that n k ′ −1 ′−1 t = ′ ′ ′ t (k+k )x1+(k +k )x−1 Cn dζ (k, k )kc(k ) c(k) − ′ e "D k ! ′ −1−n k ′ −1 ′−1 ′ ′ ′ ′ t (k+k )x1+(k +k )x−1 = dζ (k , k)k c(k ) c(k) − e = C−1−n. "D k ! Notice that tΩ = Ω holds identically. By taking the transpose of (2.9) and (2.11), we can derive t Un = U−1−n and τn = τ−1−n, and as a result we have ϕn = −ϕ−n−2 and θn = θ−1−n. (4.2) The bilinear equation (3.7) in this case then reduces to 1 D D τ · τ = τ2 − τ τ , 2 1 −1 0 0 0 1 0 1 D D τ · τ = τ2 − τ τ , n = 1, 2, ··· , (4.3) 2 1 −1 n n n n+1 n−1 due to (4.2), and subsequently, it leads to a reduced nonlinear equation of (1.1) in the form of

ϕ0 ϕ1−ϕ0 ∂1∂−1ϕ0 = e − e ,

ϕn−ϕn−1 ϕn+1−ϕn ∂1∂−1ϕn = e − e , n = 1, 2, ··· , (4.4) which we refer to as the 2DTL of B∞-type. Equation (4.4) can also be written as its matrix form −(θ +ln 2) θ0 + ln 2 2 −1 e 0 −θ  θ1  −2 2 −1   e 1  ∂1∂−1  θ  = −  −1 2 −1   e−θ2  . (4.5)  2       .   . . .   .   .   ......   .              We can easily see that the coeﬃcient matrix is exactly the Cartan matrix associated  with B∞.  The Lax pair and the adjoint one of (4.4) still take the same form of (3.9) and (3.11), respectively, subject to the additional conditions for ϕn and θn given in (4.2). To be more precise, the linear problem takes the same form of (3.9), but with the Lax matrices given by ......    −θ1   −∂1ϕ0 1  e 0    −θ   0 1   e 0 0         −θ0  P =  ∂1ϕ0 1  and Q =  e 0  ,        e−θ1 0   ∂1ϕ1 1       −θ2   ∂1ϕ2 1   e 0       . .   .. ..   .. ..  . .         respectively. Similarly, the adjoint linear problem takes the form of (3.11), in which the Lax matrices P ′ and Q′ are given as follows: ......    −θ1   −1 ∂1ϕ1   0 −e     −θ0   −1 ∂1ϕ0   0 −e  ′   ′   P =  −1 0  , Q =  0 −e−θ0  .      −1 −∂ ϕ   0 −e−θ1   1 0     −1 −∂ ϕ   0 −e−θ2   1 1     . .   . .   .. ..  .. ..         4.2. C∞-type. For the 2D Toda equation of C∞-type, we consider a symmetric integration domain D and simultaneously impose the symmetry condition dζ(k, k′) = dζ(k′, k), ∀(k, k′) ∈ D. (4.6)

Taking the transpose of Cn, we observe that n ′ −n k ′ −1 ′−1 k ′ −1 ′−1 t = ′ ′ t (k+k )x1+(k +k )x−1 = ′ ′ t (k+k )x1+(k +k )x−1 = Cn dζ(k, k )c(k ) c(k) − ′ e dζ(k , k)c(k ) c(k) − e C−n, "D k ! "D k ! LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 9

t t t and subsequently Un = U−n. Again, observing that Ω = Ω and simultaneously making use of Cn = C−n, we obtain

τn = τ−n, and therefore, ϕn = −ϕ−n−1, and θn = θ−n. (4.7)

By performing the reduction conditions (4.7) on the 2DTL of A∞-type, the bilinear equation (3.7) and the nonlinear equation (1.1) reduce to multi-component systems 1 D D τ · τ = τ2 − τ2, 2 1 −1 0 0 0 1 1 D D τ · τ = τ2 − τ τ , n = 1, 2, ··· (4.8) 2 1 −1 n n n n+1 n−1 and

2ϕ0 ϕ1−ϕ0 ∂1∂−1ϕ0 = e − e ,

ϕn−ϕn−1 ϕn+1−ϕn ∂1∂−1ϕn = e − e , n = 1, 2, ··· , (4.9) respectively, in which the latter is referred to as the 2DTL of C∞-type. Equation (4.9) can alternatively be written in its the matrix form by the quantity θn as follows: −θ θ0 2 −2 e 0 −θ  θ1  −1 2 −1   e 1  ∂1∂−1  θ  = −  −1 2 −1   e−θ2  , (4.10)  2       .   . . .   .   .   ......   .              with the entries of the Cartan matrix of C∞ appearing as the coeﬃcients.    The Lax pair and its adjoint of (4.9) also follow from (3.9) and (3.11), respectively, subject to (4.7). In other words, the linear problem and the adjoint linear problem are still written in their respective forms (3.9) and (3.11); however, the Lax matrices in the C∞ class are given by ...... −θ  −∂1ϕ1 1   e 2 0         −θ1   −∂1ϕ0 1   e 0         −θ0  P =  ∂1ϕ0 1  , Q =  e 0  ,    −     e θ1 0   ∂1ϕ1 1       −θ1   ∂1ϕ2 1   e 0       . .   .. ..   .. ..  . .         and     ......   −θ −1 ∂1ϕ2  0 −e 2       −1 ∂ ϕ   −θ1   1 1   0 −e  ′   ′      −θ0  P =  −1 ∂1ϕ0  , Q =  0 −e  ,      − −   0 −e−θ1   1 ∂1ϕ0       −θ2   −1 −∂1ϕ1   0 −e       . .   .. ..   .. ..   . .         respectively.    The results in this section show that theorems 3.1, 3.2 and 3.3 are also applicable to the 2DTLs of B∞- and C∞-types, subject to the reductions that the arbitrary integration measure and domain in the theorems obey (4.1) and (4.6), respectively. To put it another way, we have successfully constructed the direct linearising type solutions to the B∞-type and C∞-type 2D Toda equations.

s associated with (1) (2) (1) and (2) 5. 2DTL Ar , A2r , Cr Dr+1

We have discussed the 2DTLs of A∞-, B∞- and C∞-types from the linear integral equations. In these cases, there are eﬀectively two spectral parameters, namely k and k′, playing roles in the solution structure, leading to (2+1)-dimensional integrable models. In this section, we consider the so-called periodic reduction on the integration measure, by restricting the spectral parameters k and k′ on the curve kN = (−k′)N for integer N, namely there is eﬀectively only one spectral parameter. As a result, the 2DTLs (1) (2) (1) (2) of Ar -, A2r -, Cr - and Dr+1-types are recovered, as (1+1)-dimensional integrable systems.

(1) r+1 ′ r+1 5.1. Ar -type. We perform the periodicity condition k = (−k ) on the A∞-type Toda equation. This is realised by taking a particular measure in the form of

ϕ(r+1) ′ 1 dλ j(k)dk dζ(k, k′) = (5.1) 2πi ′ ( j) Xj=1 k + ωr+1k 10 YUE YIN AND WEI FU and a particular domain D = Γ × Γ′ with Γ being an arbitrary integration domain for k and Γ′ being a closed contour in terms ′ ( j) ( j) of k that containing singularities −ωr+1k. Here ωN stand for the primitive Nth roots of unity and ϕ(·) denotes Euler’s totient (1) function. For example, when N = 4 we have ϕ(N) = 2; therefore, there are two primitive 4th roots of unity, given by ω4 = ω4 (2) 3 and ω4 = ω4, respectively, where ω4 = exp(2πi/4). By Cauchy’s residue theorem, the linear integral equations (3.12) and (3.15) under such a constraint hence turn out to be ϕ(r+1) ( j) ( j) un(k) + dλ j(l)ρn(k)Ω(k, −ωr+1l)σn(−ωr+1l)un(l) = ρn(k)c(k) (5.2) Xj=1 ZΓ j and ϕ(r+1) t ′ ′ ′ t ( j) t ′ ′ vn(k ) + dλ j(l)ρn(l)Ω(l, k )σn(k ) vn(−ωr+1l) = c(k )σn(k ), (5.3) Xj=1 ZΓ j namely linear integral equations only involve single integrations in terms of the spectral parameter k, in which Γ j and dλ j(l) are arbitrary integration domains and measures reduced from D and dζ(k, k′), respectively. We remark that at this stage we still do not fix the integration domains and measures, as this shall bring us the most general solution space for the reduced equations. Once explicit forms are selected, special classes of solutions can be constructed. For example, when the integration domains are chosen to be closed contours that contain a finite number of higher-order poles, the generalised Cauchy matrix solutions can be obtained, see section 6 for more detail. Simultaneously, the reduced measure (5.1) also leads to ϕ(r+1) t ( j) ( j) Un = dλ j(k)un(k) c(−ωr+1k)σn(−ωr+1k), (5.4) Xj=1 ZΓ j as well as ϕ(r+1) t ( j) ( j) Cn = dλ j(k)ρn(k)c(k)c(−ωr+1k)σn(−ωr+1k). (5.5) Xj=1 ZΓ j Notice that n+r+1 ( j) ( j) −1 ( j) −1 1 ( j) (k−ωr+1k)x1+(k −(ωr+1k) )x−1 ρn+r+1(k)σn+r+1(−ω + k) = e = ρn(k)σn(−ω + k). r 1 ω( j)  r 1  r+1    It is easily verified from (5.5) that Cn+r+1 = Cn. Subsequently,  

Un+r+1 = Un and τn+r+1 = τn, and futhermore we have ϕn+r+1 = ϕn and θn+r+1 = θn. (5.6)

The periodicity of τn then reduces (3.7) to a multi-component system of (1+1)-dimensional bilinear equations 1 D D τ · τ = τ2 − τ τ , 2 1 −1 0 0 0 1 r 1 D D τ · τ = τ2 − τ τ , n = 1, 2, ··· , r − 1, (5.7) 2 1 −1 n n n n+1 n−1 1 D D τ · τ = τ2 − τ τ , 2 1 −1 r r r 0 r−1 (1) (1) which is the bilinear form of the Ar -type 2D Toda system. The nonlinear 2DTL of Ar -type reads

ϕ0−ϕr ϕ1−ϕ0 ∂1∂−1ϕ0 = e − e ,

ϕn−ϕn−1 ϕn+1−ϕn ∂1∂−1ϕn = e − e , n = 1, 2, ··· , r − 1, (5.8)

ϕr−ϕr−1 ϕ0−ϕr ∂1∂−1ϕr = e − e , as it follows from (1.1) in virtue of the periodicity of ϕn given in (5.6). Equation (5.8) can alternatively be written as −θ θ0 2 −1 −1 e 0 −θ1  θ1  −1 2 −1   e   .   . . .   .  ∂1∂−1  .  = −  ......   .  , (5.9)        θ   −1 2 −1  e−θr−1   r−1           −θr   θr  −1 −1 2   e           (1)    (1) where the coeﬃcient matrix is exactly the same as the Cartan matrix for Ar . In general, the nonlinear equations for the Ar case (1) can only be expressed as multi-component systems like (5.8). However, the A1 case is an exception. Due to τn+2 = τn, we have ϕ1 = ln(τ2/τ1) = ln(τ0/τ1) = −ϕ0, and therefore, the ﬁrst equation in (5.8) is closed as a scalar equation, taking the form of (1.3), i.e. the sinh-Gordon equation. The Lax pair and the adjoint Lax pair of equation (5.8) are also derived from the measure reduction (5.1). Notice that r+1 ρn+r+1(k) = k ρn(k). Shifting n by r + 1 units in (2.8a), we obtain r+1 un+r+1(k) = (1 − Un+r+1Ω)c(k)ρn+r+1(k) = k (1 − UnΩ)c(k)ρn(k). LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 11

In other words, the wave function φn obeys the quasi-periodicity r+1 −(r+1) φn+r+1 = k φn, and similarly, ψn+r+1 = k ψn,

′ r+1 ′ r+1 where we have replaced k by k in the equation for ψn, since (5.6) implies that k = (−k ) . Hence, the Lax pair (3.9) reduces to

∂1Φ= P Φ, ∂−1Φ= QΦ, (5.10) t where Φ= (φ0,φ1, ··· ,φr), the matrices P and Q are given by −(r+1) −θ ∂1ϕ0 1 0 k e 0 −θ1  ∂1ϕ1 1  e 0   . .   . .  P =  .. ..  and Q =  .. ..  ,      ∂ ϕ 1   e−θr−1 0   1 r−1     kr+1 ∂ ϕ   −θr   1 r  e 0      respectively; the adjoint Lax pair (3.11) reduces to    ′ ′ ∂1Ψ= P Ψ, ∂−1Ψ= Q Ψ, (5.11) t ′ ′ where Ψ= (ψ0, ψ1, ··· , ψr), and P and Q are deﬁned as r+1 −θ −∂1ϕ0 −k 0 −e 0 −θ1  −1 −∂1ϕ1   0 −e  ′  . .  ′  . .  P =  .. ..  and Q =  .. ..  ,      −1 −∂ ϕ   0 −e−θr−1   1 r−1     −1 −∂ ϕ  −k−(r+1) −θr   1 r  e 0      respectively.    

(2) 5.2. A2r -type. We now take a special measure

ϕ(N) ′ ′ 1 dλ j(k)dk dkdλ j(k ) dζ(k, k′) = k − . (5.12) 2πi k′ + ω( j)k k + ω( j)k′  Xj=1  N N    This leads to the reduced linear integral equation  

ϕ(N) ( j) ( j) ( j) ( j) un(k) + dλ j(l) ρn(k)Ω(k, −ωN l)σn(−ωN l)lun(l) − ρn(k)Ω(k, l)σn(l)(−ωN l)un(−ωN l) = ρn(k)c(k), (5.13) ZΓ Xj=1 j h i and its adjoint equation

ϕ(N) t ′ ′ ′ t ( j) ( j) ( j) ′ ′ ( j) t t ′ ′ vn(k ) + dλ j(l) ρn(l)Ω(l, k )σn(k )l vn(−ωN l) − ρn(−ωN l)Ω(−ωN l, k )σn(k )(−ωN l) vn(l) = c(k )σn(k ). (5.14) ZΓ Xj=1 j h i

Likewise, the inﬁnite matrices Un and Cn under such a reduction become ϕ(N) t ( j) ( j) ( j) ( j) t Un = dλ j(k) kun(k) c(−ωN k)σn(−ωN k) − (−ωN k)un(−ωN k) c(k)σn(k) (5.15) ZΓ Xj=1 j h i and ϕ(N) t ( j) ( j) ( j) ( j) ( j) t Cn = dλ j(k) kρn(k)c(k)c(−ωN k)σn(−ωN k) − (−ωN k)ρn(−ωN k)c(−ωN k) c(k)σn(k) , (5.16) ZΓ Xj=1 j h i respectively. (2) t For the A2r case, we set N = 2r + 1. Straightforward computation shows that Cn = C−1−n, as well as Cn+2r+1 = Cn, which t imply that Un+2r+1 = Un and Un = U−1−n. Based on these relations, we obtain the reductions on τn, ϕn as well as θn as follows:

τn+2r+1 = τn, τn = τ−1−n, ϕn+2r+1 = ϕn, ϕn = −ϕ−2−n, θn+2r+1 = θn, θn = θ−1−n. (5.17)

Performing the reduction of τn on (3.7), we immediately obtain the bilinear system 1 D D τ · τ = τ2 − τ τ , 2 1 −1 0 0 0 1 0 1 D D τ · τ = τ2 − τ τ , n = 1, 2, ··· , r − 1, (5.18) 2 1 −1 n n n n+1 n−1 1 D D τ · τ = τ2 − τ2 . 2 1 −1 r r r r−1 12 YUE YIN AND WEI FU

(2) Similarly, the reduction on ϕn provides us with the 2DTL of A2r -type, i.e.

ϕ0 ϕ1−ϕ0 ∂1∂−1ϕ0 = e − e ,

ϕn−ϕn−1 ϕn+1−ϕn ∂1∂−1ϕn = e − e , n = 1, 2, ··· , r − 2, (5.19)

ϕr−1−ϕr−2 −2ϕr−1 ∂1∂−1ϕr−1 = e − e , and its Cartan matrix form is given by −(θ +ln 2) θ0 + ln 2 2 −1 e 0 −θ1  θ1  −2 2 −1   e   .   . . .   .  ∂ ∂−  .  = −  . . .   .  . (5.20) 1 1  .   . . .   .       −θr−1   θr−1   −1 2 −1  e       −θ   θr   −2 2   e r              When r = 1, equation (5.19) turns out to be a closed-form scalar equation due to ϕ1 = −ϕ0, which is the nothing but the Tzitzeica equation (1.4a). 2r+1 −(2r+1) Observing that the wave functions φn and ψn still satisfy the periodicity conditions φn+2r+1 = k φn and ψn+2r+1 = k ψn, we obtain the associated linear problems for equation (5.19), including the linear problem

∂1Φ= P Φ, ∂−1Φ= QΦ (5.21) t with Φ= (φ0,φ1, ··· ,φ2r) and Lax matrices

−(2r+1) −θ0 ∂1ϕ0 1 0 k e . . e−θ1 0  .. ..       .. ..   ∂1ϕr−1 1   . .     −θ  P =  −∂ ϕ − 1  and Q =  e r 0  ,  1 r 1     . .   e−θr−1 0   .. ..       . .   −∂ ϕ 1  .. ..   1 0    k2r+1 0  e−θ0 0          as well as the adjoint linear problem    ′ ′ ∂1Ψ= P Ψ, ∂−1Ψ= Q Ψ, (5.22)

t in which Ψ= (ψ0, ψ1, ··· , ψ2r), and 0 −k2r+1 0 −e−θ0 −1 −∂ ϕ . .  1 0   .. ..   . .     .. ..   0 −e−θr      ′ =  − −  ′ =  − −θr−1  P  1 ∂1ϕr−1  , Q  0 e  .  −1 ∂ ϕ   . .   1 r−1   .. ..   . .     .. ..   0 −e−θ1       −1 ∂ ϕ  −k−(2r+1)e−θ0 0   1 0        We comment that the measure (5.12) in a sense is a combination of (4.1) and (5.1), as it leads to a reduction which simulta- neously satisﬁes the anti-symmetry and periodicity conditions. However, taking such a measure is not the only reduction to the (2) 2DTL of A2r -type. In fact, we can alternatively consider the measure

ϕ(N) ′ ′ 1 dλ j(k)dk dk dλ j(k ) dζ(k, k′) = + , (5.23) 2πi k′ + ω( j)k k + ω( j)k′  Xj=1  N N    which should be thought of as imposing (4.6) and (5.1) simultaneously. This brings us the linear integral equation

ϕ(N) ( j) ( j) ( j) un(k) + dλ j(l) ρn(k)Ω(k, −ωN l)σn(−ωN l)un(l) + ρn(k)Ω(k, l)σn(l)un(−ωN l) = ρn(k)c(k), (5.24) ZΓ Xj=1 j h i and the adjoint linear integral equation

ϕ(N) t ′ ′ ′ t ( j) ( j) ( j) ′ ′ t ′ ′ vn(k ) + dλ j(l) ρn(l)Ω(l, k )σn(k ) vn(−ωN l) + ρn(−ωN l)Ω(−ωN l, k )σn(k ) vn(l) = c(k )σn(k ), (5.25) ZΓ Xj=1 j h i as well as the inﬁnite matrices ϕ(N) t ( j) ( j) ( j) t Un = dλ j(k) un(k) c(−ωN k)σn(−ωN k) + un(−ωN k) c(k)σn(k) (5.26) ZΓ Xj=1 j h i LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 13 and ϕ(N) t ( j) ( j) ( j) ( j) t Cn = dλ j(k) ρ(k)c(k)c(−ωN k)σ(−ωN k) + ρ(−ωN k)c(−ωN k) c(k)σ(k) . (5.27) ZΓ Xj=1 j h i Therefore, by setting N = 2r + 1 we can follow the same procedure and ﬁnd reduction conditions

τn+2r+1 = τn, τn = τ−n, ϕn+2r+1 = ϕn, ϕn = −ϕ−n−1, θn+2r+1 = θn, θn = θ−n as well as 2r+1 −(2r+1) φn+2r+1 = k φn and ψn+2r+1 = k ψn

(2) on the 2DTL of A∞-type, leading to the equivalent results in the A2r class, compared with (5.18), (5.19), (5.20), (5.21) as well as (5.22). For example, in this case, the simplest model expressed by the potential ϕn when r = 1 gives rise to another form of the Tzitzeica equation, namely (1.4b).

(1) (1) 5.3. Cr -type. The Cr -type 2D Toda equation follows from the 2r-periodic reduction of the C∞ case. For this reason, we set t N = 2r in (5.23), (5.24), (5.25), (5.26) and (5.27). We can directly verify that Cn+2r = Cn and Cn = C−n and subsequently t Un+2r = Un and Un = U−n, i.e. both quantities satisfy the symmetry and periodicity conditions simultaneously. Based on these conditions, we conclude that

τn = τn+2r, τn = τ−n, ϕn+2r = ϕn, ϕn = −ϕ−n−1, θn+2r = θn, and θn = θ−n. (5.28) Performing these reductions, we obtain the bilinear system 1 D D τ · τ = τ2 − τ2, 2 1 −1 0 0 0 1 1 D D τ · τ = τ2 − τ τ , n = 1, 2, ··· , r − 1, (5.29) 2 1 −1 n n n n+1 n−1 1 D D τ · τ = τ2 − τ2 , 2 1 −1 r r r r−1 and the nonlinear system

2ϕ0 ϕ1−ϕ0 ∂1∂−1ϕ0 = e − e ,

ϕn−ϕn−1 ϕn+1−ϕn ∂1∂−1ϕn = e − e , n = 1, 2, ··· , r − 2 (5.30)

ϕr−1−ϕr−2 −2ϕr−1 ∂1∂−1ϕr−1 = e − e ,

(1) from (3.7) and (1.1), respectively, in which the latter is referred to as the Cr -type 2DTL. Alternatively, we can reformulate (5.30) as −θ θ0 2 −2 e 0 −θ1  θ1  −1 2 −1   e   .   . . .   .  ∂1∂−1  .  = −  ......   .  ,        θ   −1 2 −1  e−θr−1   r−1           −θr   θr   −2 2   e          (1)     from which we can observe the Cartan matrix of Cr . We comment that the relation ϕ−1 = −ϕ0 holds identically for r = 1, and (1) in this case equation (5.30) turns out to be the sinh–Gordon equation (which has been discussed in the A1 class. Therefore, we normally start with r = 2 in this class. 2r −2r 2r ′ 2r Notice that the periodicity of Un provides us with φn+2r = k φn and ψn+2r = k ψn, where we have made use of k = (−k ) . We obtain from (3.9) and (3.11) the Lax pair

∂1Φ= P Φ, ∂−1Φ= QΦ (5.31) t with Φ= (φ0,φ1, ··· ,φ2r−1), where −2r −θ ∂1ϕ0 1 0 k e 0 . . e−θ1 0  .. ..       .. ..   ∂1ϕr−1 1   . .     −θ  P =  −∂ ϕ − 1  , Q =  e r 0  ,  1 r 1     . .   e−θr−1 0   .. ..       . .   −∂ ϕ 1   .. ..   1 1     k2r −∂ ϕ   e−θ1 0   1 0       as well as the adjoint Lax pair    ′ ′ ∂1Ψ= P Ψ, ∂−1Ψ= Q Ψ, (5.32) 14 YUE YIN AND WEI FU

t with Ψ= (ψ0, ψ1, ··· , ψ2r−1), in which

2r −θ ∂1ϕ0 −k 0 −e 0 −1 −∂ ϕ . .  1 0   .. ..   . .     .. ..   0 −e−θr      ′ =  − −  ′ =  − −θr−1  P  1 ∂1ϕr−1  , Q  0 e  ,  −1 ∂ ϕ   . .   1 r−1   .. ..   . .     . .   −θ2   . .   0 −e     −2r −θ1   −1 ∂1ϕ1 −k e 0          respectively, as the two linear problems for the 2DTL of C∞-type.

(2) (2) 5.4. Dr+1-type. Finally, we discuss the Dr+1-type 2DTL. This is done by taking N = 2r + 2 in (5.12), (5.13), (5.14), (5.15) and (2) (5.16), which results in the fact that in the Dr+1 class the reduction conditions are as follows:

τn = τn+2r+2, τn = τ−1−n, ϕn = ϕn+2r+2, ϕn = −ϕ−2−n, θn = θn+2r+2, θn = θ−1−n. (5.33) Making use of the constraints in (5.33), we obtain the reduced bilinear system 1 D D τ · τ = τ2 − τ τ , 2 1 −1 0 0 0 1 0 1 D D τ · τ = τ2 − τ τ , n = 1, 2, ··· , r − 1, (5.34) 2 1 −1 n n n n+1 n−1 1 D D τ · τ = τ2 − τ τ , 2 1 −1 r r r r r−1 the nonlinear system

ϕ0 ϕ1−ϕ0 ∂1∂−1ϕ0 = e − e ,

ϕn−ϕn−1 ϕn+1−ϕn ∂1∂−1ϕn = e − e , n = 1, 2, ··· , r − 2, (5.35)

ϕr−1−ϕr−2 −ϕr−1 ∂1∂−1ϕr−1 = e − e , as well as the Cartan matrix form as follows: −(θ +ln 2) θ0 + ln 2 2 −1 e 0 −θ1  θ1  −2 2 −1   e   .   . . .   .  ∂1∂−1  .  = −  ......   .  . (5.36)        θ   −1 2 −2  e−θr−1   r−1       θ +   −1 2   e−(θr+ln 2)   r ln 2             (1)      The r = 1 case in this class also leads to the A1 class. This is because when r = 1, we have ϕ1 = −ϕ1, i.e. ϕ1 = 0, according to (5.33). Thus, we obtain from (5.35) a closed-form scalar equation

ϕ0 −ϕ0 ∂1∂−1ϕ0 = e − e , which is exactly the same as the sinh–Gordon equation (1.3), up to a scaling transformation. In other words, we have to start (2) with D3 in order to obtain nontrivial equations. 2r+2 The linear problems for (5.35) are derived accordingly, with the help of the quasi-periodicity conditions φn+2r+2 = k φn and −(2r+2) ψn+2r+2 = k ψn. The Lax pair takes the form of

∂1Φ= P Φ, ∂−1Φ= QΦ, (5.37)

t for Φ= (φ0, ··· ,φ2r+1), where the Lax matrices P and Q are as follows: ∂ ϕ 1 1 0 0 k−(2r+2)e−θ0 . .  . .  −θ1 . . e 0     ∂1ϕr−1 1   . .     .. ..   0 1       −θr  P =   , Q =  e 0  .  −∂1ϕr−1 1       −θr     e 0   .. ..     . .   .. ..     . .   −∂1ϕ0 1      e−θ0 0  k2r+2 0           The adjoint Lax pair is composed of the linear equations 

′ ′ ∂1Ψ= P Ψ and ∂−1Ψ= Q Ψ, (5.38) LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 15

t ′ ′ where Ψ= (ψ0, ··· , ψ2r+1), and the inﬁnite matrices P and Q take their respective forms of 0 −k2r+2 0 −e−θ0 −1 −∂ ϕ  1 0  . .  . .   .. ..   .. ..       −θr     0 −e  ′  −1 −∂1ϕr−1  ′  −θ  P =   and Q =  0 −e r  .  −     1 0       .. ..   −1 ∂1ϕr−1   . .     −θ   . .   0 −e 1   .. ..      − −(2r+2) −θ0     k e 0   −1 ∂1ϕ0          (1) (2) (1) (2) We have established the DL scheme for the 2D Toda-type equations of Ar+1-, A2r -, Cr - and Dr+1-types for positive integers r, by identifying the corresponding linear integral equations for each class of nonlinear equations. The main results in the theorems in section 3 are inherited for these classes of 2DTLs, subject to their respective measure degenerations, leading to the direct linearising type solutions.

6. General formulae of Cauchy matrix type solutions We now discuss a special class of solutions called Cauchy matrix solutions for all the 2D Toda-type systems. These solutions are constructed by taking particular measures containing an arbitrary number of singularities of arbitrary order, as special cases of the direct linearising type solutions presented in the previous sections. In this section, our aim is to propose a general formula of the Cauchy matrix solution for each class of the 2DTL. Notice that the potential for the 2D Toda-type equations has a common form

(0,−1) ϕn ln 1 − Un , (6.1) ′ (¯i,¯j) with ρn(k) and σn(k ) deﬁned as (3.1). Hence, below we only present the general formula of Un for arbitrary integers ¯i and ¯j, class by class. For convenience, we omit the suﬃx n in our formulae, if not necessary.

6.1. A∞-type. We consider the measure N N′ ′ ′ 1 1 1 ′ ζ k, k = A ′ s − s ′ − k k , d ( ) j, j ( j 1)!( j 1)! 2 s s′ d d (6.2) (2πi) (k − k j) j ′ ′ j′ Xj=1 Xj′=1 (k − k j′ ) and simultaneously the integration domain D = Γ × Γ′ in which Γ and Γ′ are the two separate contours on the k- and k′-planes ′ ′ ′ ′ containing singularities k j of order s j and k j of order s j′ inside for j = 1, ··· , N and j = 1, ··· , N , respectively. The residue theorem then reduces (3.12) and (3.14) to N N′ ′ ′ − ρ(k)σ(k ′ ) s j′ 1 j s j−1 ′ u(k) + A j, j ∂k′ ′ ∂k u(k j) = ρ(k)c(k) (6.3)  j′ +  j ′ k k j′ Xj=1 Xj =1     and   N N′ s′ −1 s j−1 j′ t ′ ′ ′ U = A j, j ∂k u(k j) ∂ ′ c(k j′ )σ(k j′ ) , (6.4) j k j′ Xj=1 Xj′=1 respectively. For convenience below, we introduce ﬁnite matrices

A1,1 ··· A1,N′ M1,1 ··· M1,N ′ s′ −1 ρ(k j)σ(k ′ ) . . . . s j−1 j′ j A =  . .  and M =  . .  where M j′, j = ∂ ∂ ′ , . . . . k j k ′ ′     j k + k ′     j j A ··· A ′  M ′ ··· M ′   N,1 N,N   N ,1 N ,N     diagonal block matrices    C C C C ′ tC ′ tC ′ = diag( 1, ··· , N ) and = diag( 1 , ··· , N′ ), in which ′ s ′ −1 s j−1t j t ′ ∂ c(k j) ∂k′ c(k j′ ) C t t k j tC ′  j′ t ′  j =  c(k j), ··· ,  , j = ′ , ··· , c(k j′ ) , (s − 1)!  (s ′ − 1)!   j   j          diagonal matrices   ′ ′ ′ D = diag((s1 − 1)!, ··· , (sN − 1)!) and D = diag((s1 − 1)!, ··· , (sN′ − 1)!), as well as ﬁnite row and column vectors t t t t t t r = ( r1, ··· , rN ) and s = ( s1, ··· , sN′ ), 16 YUE YIN AND WEI FU where ′ s ′ −1 s j−1 j ′ ∂ ρ(k j) ∂k′ σ(k j′ ) t k j t ′ j′ r = , ··· , ρ(k ) and s ′ = σ(k ′ ), ··· ,  . j  − j  j  j ′ −   (s j 1)!   (s j′ 1)!          With these objects, equation (6.4) is expressed as   

t ′ ′ t s1−1 sN−1 U = U AD C s, where U = ∂ u(k1), ··· ,∂ u(kN) . (6.5) k1 kN Meanwhile, from equation (6.3) we can write down tU + tU AM = trC D. (6.6) This in turn implies that we are able to eliminate tU in (6.5) through (6.6) and obtain U = trC D(I + AM)−1AD′C ′s. From the ﬁrst glance, the sizes of the left and right sides of the above equation do not match. But in fact, the multiplication on the right hand side should be thought of as that of block matrices; in other words, the right hand side reads ′ tC ′ (s1 − 1)! 1 s1 (s − 1)!tr C , ··· , (s − 1)!tr C (1 + AM)−1A  .  1 1 1 N N N  .   ′ t ′   C ′  (sN′ − 1)! N′ sN    which is of size  

(block row vector with components of size ∞× 1)1×N

× (ﬁnite matrix with entries of size 1 × 1)N×N′

× (block column vector with components of size 1 ×∞)N′×1, because the jth component in the block row vector takes the form of an ∞× 1 column vector s −1−i s j−1 j i ∂k ρ(k j) ∂k c(k j) (s − 1)!tr C = (s − 1)! j j , j j j j (s − 1 − i)! i! Xi=0 j and the j′th component in the block column vector takes the form of a 1 ×∞ row vector ′ − − ′ ′ s j′ 1 i ′ s ′ −1 t ′ i ′ j ∂k′ c(k j′ ) ∂k′ σ(k j′ ) ′ t ′ ′ j′ j′ ′ C ′ ′ = ′ (s j − 1)! j s j (s j − 1)! ′ ′ ′ , (s ′ − 1 − i )! i ! Xi′=0 j according to the multiplication of block matrices. Thus, the right hand side of U is indeed of size ∞×∞. Taking the (¯i, ¯j)-entry of U, we ﬁnally reach to the general solution structure U (¯i,¯j) = trK(¯i)D(I + AM)−1AD′K′(¯j)s, (6.7) where K(¯i) and K′(¯j) are deﬁned as (¯i) (¯i) (¯i) ′(¯j) t ′(¯j) t ′(¯j) K = diag K1 , ··· , KN and K = diag K 1 , ··· , K N′ , respectively, with the blocks given as follows: ′ s ′ −1 ¯j s j−1 ¯i j ′ ∂ k ∂ ′ k ′ ¯ k j ¯ k ′ j ¯ K(i) = t k¯i , ··· , j , tK′( j) =  j , ··· , k′ j  . j  j −  j′  ′ − j′   (s j 1)! (s j′ 1)!              6.2. B∞-type. We respect (4.1) and take a degenerate measure N ′ 1 1 1 ′ dζ(k, k ) = A ′ (s − 1)!(s ′ − 1)! kdkdk , where A ′ = −A ′ , (6.8) j, j j j 2 s ′ s ′ j, j j , j (2πi) (k − k ) j (k − k ′ ) j jX, j′=1 j j ′ ′ and the integration domain D = Γ × Γ with Γ and Γ containing singularities k j of order s j and k j′ of order s j′ on their respective complex planes, for j = 1, ··· , N and j = 1, ··· , N. The residue theorem reduces the linear integral equation (3.12) and the potential (3.14) to N ′ s j′ −1 ρ(k)σ(k j ) s j−1 u + ′ ∂ ∂ u = ρ c (k) A j, j k ′ k k j (k j) (k) (k) (6.9) j k + k ′ ! j jX, j′=1 j and N s j−1 s j′ −1t U = A j, j′ ∂ k ju(k j) ∂ c(k j′ )σ(k j′ ) , (6.10) k j k j′ jX, j′=1 LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 17 respectively. For convenience, we introduce an antisymmetric matrix A and the Cauchy matrix M given by

A1,1 ··· A1,N M1,1 ··· M1,N . . . . A =  . .  and M =  . .  ,     AN,1 ··· AN,N MN,1 ··· MN,N          respectively, in which

′ s j−1 s j′ −1 k jρ(k j)σ(k j ) ′ = − ′ , ′ = ∂ ∂ , A j, j A j j M j , j k k ′ j j k j + k j′ as well as a diagonal matrix D and diagonal block matrices as follows: C C C C ′ tC ′ tC ′ D = diag((s1 − 1)!, ··· , (sN − 1)!), = diag( 1, ··· , N), = diag( 1 , ··· , N ), C tC ′ where the blocks j and j are deﬁned as follows:

s j−1 t s j′ −1t ∂ k j c(k j) ∂ c(k j′ ) C t t k j tC ′ k j′ t j = k j c(k j), ···  , j =  , ··· , c(k j′ ) . (s j − 1)!  (s j′ − 1)!              Simultaneously, we also adopt the notations tr and s for a row vector and a column vector, deﬁned by

t t t t t t r = ( r1, ··· , rN) and s = ( s1, ··· , sN ), respectively, where

s j−1 s j′ −1 ∂ ρ(k j) ∂ σ(k j′ ) t k j t k j′ r j =  , ··· , ρ(k j) and s j′ = σ(k j′ ), ··· ,  . (s − 1)! (s ′ − 1)!  j   j          With these new symbols, we can rewrite U in (6.10) as   

t ′ t s1−1 sN−1 U = U ADC s, where U = ∂ k1u(k1), ··· ,∂ kN u(kN) . (6.11) k1 kN At the same time, we obtain from (6.9) the matrix equation

tU + tU AM = trC D, (6.12) which together with (6.11) provides us with the Cauchy matrix expression of U as follows:

U = trC D(I + AM)−1ADC ′s.

Then the (¯i, ¯j)-entry of U takes the expression

U (¯i,¯j) = trK(¯i)D(I + AM)−1ADK′(¯j)s. (6.13)

Here K(¯i) and K′(¯j) are block diagonal matrices

(¯i) (¯i) (¯i) ′(¯j) t ′(¯j) t ′(¯j) K = diag K1 , ··· , KN and K = diag K 1 , ··· , K N , where the blocks are ﬁnite vectors as follows:

s j−1 ¯i+1 s j′ −1 ¯j ∂ k ∂ k ′ (¯i) t ¯i+1 k j j t ′(¯j) k j′ j ¯j K j = k j , ··· ,  , K j′ =  , ··· , k j′  .  (s j − 1)!  (s j′ − 1)!              6.3. C∞-type. We consider a symmetric measure in the C∞ class following (4.6) as follows:

N ′ 1 1 1 ′ dζ(k, k ) = A ′ (s − 1)!(s ′ − 1)! dkdk , where A ′ = A ′ , (6.14) j, j j j 2 s ′ s ′ j, j j , j (2πi) (k − k ) j (k − k ′ ) j jX, j′=1 j j and still let the domain D being the tensor of two contours containing the singularities k j and k j′ on their respective complex planes. By comparing (6.14) with (6.2), we observe that the only diﬀerence is that additional constraints

′ ′ ′ A j, j′ = A j′, j, k j′ = k j′ , s j′ = s j′ and N = N (6.15) are imposed in the C∞ case. Therefore, equation (6.7) subject to (6.15) provides the general formula for the Cauchy matrix solution of the 2DTL of C∞-type. 18 YUE YIN AND WEI FU

(1) (1) 6.4. Ar -type. Recall that the linear integral equations and the potential of Ar -type are given by (5.2) and (5.4). To construct the Cauchy matrix solution, we take a special measure in the form of N j 1 1 = ′ ′ − dλ j(k) A j, j (s j, j 1)! s ′ dk, (6.16) 2πi (k − k ′ ) j, j Xj′=1 j, j and let k j, j′ be the singularities of order s j, j′ surrounded by their respective contours Γ j. And then the residue theorem implies that the linear integral equation (5.2) reduces to

ϕ(r+1) N j ( j) s j, j′ −1 ′ d ρ(k)σ(−ωr+1k j, j ) u(k) + A j, j′ u(k j, j′ ) = ρ(k)c(k). (6.17) s j, j′ −1 ( j)  ′  ′ dk ′ k − ω k j, j Xj=1 Xj =1 j, j  r+1    Simultaneously, the potential U under (6.16) turns out to be 

ϕ(r+1) N j s j, j′ −1 d t ( j) ( j) U = A j, j′ u(k j, j′ ) c(−ω k j, j′ )σ(−ω k j, j′ ) . (6.18) s j, j′ −1 r+1 r+1 Xj=1 Xj′=1 dk j, j′ For convenience, we introduce block vectors

s j, j′ −1 ∂ u(k j, j′ ) t t t t t t k j, j′ U U U U U U ′ ′ = 1,1, ··· , 1,N1 ; ··· ; ϕ(r+1),1, ··· , ϕ(r+1),Nϕ(r+1) , in which j, j =  , ··· , u(k j, j ) , (s ′ − 1)!  j, j    and  

s j, j′ −1 ∂ ρ(k j, j′ ) t t t t t t k j, j′ ′ ′ r = r1,1, ··· , r1,N1 ; ··· ; rϕ(r+1),1, ··· , rϕ(r+1),Nϕ(r+1) , with r j, j =  , ··· , ρ(k j, j ) , (6.19) (s ′ − 1)!  j, j    as well as  

s j, j′ −1 ( j) ∂ σ(−ω k j, j′ ) t t t t t t ( j) k j, j′ r+1 = ··· ··· ··· ′ = − ′ ··· s s1,1, , s1,N1 ; ; sϕ(r+1),1, , sϕ(r+1),Nϕ(r+1) , with s j, j σ( ωr+1k j, j ), ,  , (6.20) (s ′ − 1)!  j, j    and also block diagonal matrices  

′ ′ A = diag A1,1, ··· , A1,N1 ; ··· ; Aϕ(r+1),1, ··· , Aϕ(r+1),Nϕ(r+1) , in which A j, j = A j, j Is j, j′ ×s j, j′ , and ′ ′ ′ ′ ′ C = diag C , ··· , C ; ··· ; C + , ··· , C + and C = diag C , ··· , C ; ··· ; C , ··· , C , 1,1 1,N1 ϕ(r 1),1 ϕ(r 1),Nϕ(r+1) 1,1 1,N1 ϕ(r+1),1 ϕ(r+1),Nϕ(r+1) where t ( j) c k ′ ′ ( j, j ) c(−ωr+1k j, j ) . . . . C  . ..  C ′  . ..  j, j′ =   and j, j′ =   ,  s j, j′ −1   s ′ −1   c ′   j, j tc ( j)   ∂k (k j, j )   ∂ (−ω k j, j′ )   j, j′   k j, j′ r+1 t ( j)   ··· c(k ′ )  ··· − ′   (s ′ −1)! j, j   ′ c( ω + k j, j )  j, j   (s j, j −1)! r 1      as well as    

′ ′ D = D1,1, ··· , D1,N1 ; ··· ; Dϕ(r+1),1, ··· , Dϕ(r+1),Nϕ(r+1) , with D j, j = (s j, j − 1)!Is j, j′ ×s j, j′ , and ﬁnally a Cauchy matrix

M1,1;1,1 ··· M1,1;1,N1 ··· M1,1;ϕ(r+1),1 ··· M1,1;ϕ(r+1),Nϕ(r+1) . . . .  . . . .     M ,N , ··· M ,N ,N ··· M ,N r+ , ··· M ,N r+ ,N   1 1;1 1 1 1;1 1 1 1;ϕ( 1) 1 1 1;ϕ( 1) ϕ(r+1)   . . . .  M =  . . . .  ,    M ··· M ··· M ··· M   ϕ(r+1),1;1,1 ϕ(r+1),1;1,N1 ϕ(r+1),1;ϕ(r+1),1 ϕ(r+1),1;ϕ(r+1),Nϕ(r+1)   . . . .   . . . .    M ··· M ··· M ··· M   ϕ(r+1),Nϕ(r+1);1,1 ϕ(r+1),Nϕ(r+1);1,N1 ϕ(r+1),Nϕ(r+1);ϕ(r+1),1 ϕ(r+1),Nϕ(r+1);ϕ(r+1),Nϕ(r+1)    whose blocks are deﬁned as  ( j) ( j) ′ ′ ′ ′ ′ ′ M j, j ;i,i = S(−ωr+1k j, j )G(ki,i , −ωr+1k j, j )R(ki,i ), in which

( j) si,i′ −1 σ(−ω k j, j′ ) ∂ ρ(ki,i′ ) r+1 ki,i′ ··· ρ k ′ . . (s ′ −1)! ( i,i ) ( j)  . .   i,i  ′ . . ′ S(−ωr+1k j, j ) =   , R(ki,i ) =  . .  ,  s j, j′ −1 ( j)   . .   ∂ σ(−ω k ′ )   . .   k ′ r+1 j, j     j, j ( j)     ··· σ(−ω k j, j′ )  ρ(ki,i′ )   (s j, j′ −1)! r+1            LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 19 and 0 0 Cs −1 C0 i,i′ ( j) ··· ( j) s ′ s ′ −1 k ′ −ω k ′ (k ′ −ω k ′ ) i,i (−1) i,i  i,i r+1 j, j i,i r+1 j, j  ( j)  . .  G(k ′ , −ω k ′ ) = (s ′ − 1)!  . .  . i,i r+1 j, j j, j  . .   ( j) s ′ −1 s j, j′ −1 ( j) s ′ −1 s j, j′ −1   (−ω ) j, j C (−ω ) j, j C   r+1 s ′ −1 r+1 s ′ +s ′ −2   j, j i,i j, j   ( j) s ′ s ′ −1 ··· ( j) s ′ +s ′ −1 s ′ +s ′ −2   ′ ′ j, j j, j ′ ′ i,i j, j i,i j, j   (ki,i −ωr+1k j, j ) (−1) (ki,i −ωr+1k j, j ) (−1)    With these symbols, we can reformulate equations (6.18) and (6.17) as  U = tU ADC ′s and tU (I + AM) = trC , respectively, which in turn implies U = trC (I + AM)−1ADC ′s; in other words, the (¯i, ¯j)-entry takes the expression U (¯i,¯j) = trK(¯i)(I + AM)−1ADK′(¯j)s, (6.21) where K(¯i) and K′(¯j) are diagonal block matrices as follows:

(¯i) (¯i) (¯i) (¯i) (¯i) K = diag K1,1, ··· , K1,N ; ··· ; Kϕ(r+1),1, ··· , Kϕ(r+1),N , 1 ϕ(r+1) ′(¯j) ′(¯j) ′(¯j) ′(¯j) ′(¯j) K = diag K1,1 , ··· , K1,N ; ··· ; Kϕ(r+1),1, ··· , Kϕ(r+1),N , 1 ϕ(r+1) in which their respective blocks are deﬁned by

¯i ( j) ¯j k ′ ′ j, j (−ωr+1k j, j ) . .  . .  (¯i)  . ..  ′(¯j) . . K ′ =  .  and K ′ =  . .  , j, j  s ′ −1  j, j    j, j ¯i   s j, j′ −1 ( j) ¯j   ∂ k ′   ∂ −ω k ′   k ′ j, j   k ′ ( r+1 j, j ) ¯j  j, j ¯i   j, j ( j)   ··· k ′   ··· − ′   (s ′ −1)! j, j   ′ ( ω k j, j )   j, j   (s j, j −1)! r+1      respectively.    

(2) (2) (2) (2) 6.5. A2r - and Dr+1-types. The linear integral equations for the A2r -type and Dr+1-type 2DTL can be written in a uniform way (2) (2) as (5.13), in which N = 2r + 1 and N = 2r + 2 for positive integers r corresponding to the A2r and Dr+1 classes, respectively. Similarly, the inﬁnite matrix U takes the form of (5.15). To derive the Cauchy matrix solution, we still take the measure (6.16) and replace this in (5.13) and (5.15) simultaneously. By following a similar derivation, we end up with the general formula for U (¯i,¯j) as follows: (¯i) (1,1) (1,2) −1 ′(¯j) (¯i,¯j) t t ′ K 0 A 0 M M A 0 D 0 K 0 s U = ( r, r ) ¯ I + ¯ . (6.22) 0 K′′(i)! " 0 −A′! M (2,1) M (2,2)!# 0 −A′! 0 D! 0 K′′′( j)! s′ ! In formula (6.22), the Ks are the diagonal block matrices given as follows:

(¯i) (¯i) (¯i) (¯i) (¯i) K = diag K1,1, ··· , K1,N ; ··· ; Kϕ(N),1, ··· , Kϕ(N),N , 1 ϕ(N) ′(¯i) ′(¯i) ′(¯i) ′(¯i) ′(¯i) K = diag K 1,1, ··· , K 1,N ; ··· ; K ϕ(N),1, ··· , K ϕ(N),N , 1 ϕ(N) ′′(¯i) ′′(¯i) ′′(¯i) ′′(¯i) ′′(¯i) K = diag K1,1 , ··· , K1,N ; ··· ; Kϕ(N),1, ··· , Kϕ(N),N , 1 ϕ(N) ′′′(¯i) ′′′(¯i) ′′′(¯i) ′′′(¯i) ′′′(¯i) K = diag K 1,1, ··· , K 1,N ; ··· ; K ϕ(N),1, ··· , K ϕ(N),N , 1 ϕ(N) in which the blocks lower triangle matrices are deﬁned as ¯i ( j) ¯i k j, j′ k j, j′ (−ωN k j, j′ ) . . . . (¯i)  . ..  ′(¯i)  . ..  K ′ =  .  , K ′ =  .  , j, j  s ′ −1  j, j  (s ′ −1)   j, j ¯i   j, j ( j) ¯i   ∂ k ′   ∂ k j, j′ (−ω k j, j′ )   k j, j′ j, j ¯i   k j, j′ N ( j) ¯i  ··· k ′   ··· k j, j′ (−ω k j, j′ )   (s j, j′ −1)! j, j   (s j, j′ −1)! N      and     ¯ ( j) i ( j) ¯i+1 (−ωN k j, j′ ) −ωN k j, j′  . .  . . ′′(¯i) . .. ′′′(¯i)  . ..  K ′ =  .  , K ′ =  .  , j, j  − ¯  j, j    s j, j′ 1 ( j) i   (s j, j′ −1) ( j) ¯   ′   i+1   ∂ (−ωN k j, j )   ∂ (−ω k ′ )   ( j) ¯  k j, j′ N j, j ( j) ¯   −ω k j, j′ ( j) i  i+1  N   ··· −ω k ′   ··· (−ω k j, j′ )   (s ′ −1)! N j, j   (s j, j′ −1)! N   j, j      respectively. The A, A′ take their respective forms of    A = diag A , ··· , A ; ··· ; A , ··· , A and A′ = diag A′ , ··· , A′ ; ··· ; A′ , ··· , A′ , 1,1 1,N1 ϕ(N),1 ϕ(N),Nϕ(N) 1,1 1,N1 ϕ(N),1 ϕ(N),Nϕ(N) (6.24) 20 YUE YIN AND WEI FU where

′ s ′ −1 ′ ′ ′ j, j ′ A j, j = A j, j Is j, j′ ×s j, j′ and A j, j′ = diag A j, j (−ωN ) , ··· , A j, j . The D is a diagonal matrix in the form of

′ ′ D = D1,1, ··· , D1,N1 ; ··· ; Dϕ(r+1),1, ··· , Dϕ(r+1),Nϕ(r+1) , with D j, j = (s j, j − 1)!Is j, j′ ×s j, j′ . (6.25) The tr and tr′ are ﬁnite row vectors deﬁned as tr = tr , ··· , tr ; ··· ; tr , ··· , tr and tr′ = tr′ , ··· , tr′ ; ··· ; tr′ , ··· , tr′ , (6.26) 1,1 1,N1 ϕ(N),1 ϕ(N),Nϕ(N) 1,1 1,N1 ϕ(N),1 ϕ(N),Nϕ(N) respectively, in which

′ ′ − s j, j −1 ( j) s j, j 1 ′ ′ ∂ ( j) ρ(−ω k j, j ) ∂ ′ ρ(k j, j ) −ω k ′ N t k j, j t ′ N j, j ( j) r j, j′ =  , ··· , ρ(k j, j′ ) , r j, j′ =  , ··· , ρ(−ωN k j, j′ ) .  (s j, j′ − 1)!   (s j, j′ − 1)!          The s and s′ are ﬁnite column vectors given by    t t t t t ′ t t ′ t ′ t ′ t ′ s = s , ··· , s ; ··· ; s N , ··· , s N and s = s , ··· , s ; ··· ; s , ··· , s , (6.27) 1,1 1,N1 ϕ( ),1 ϕ( ),Nϕ(N) 1,1 1,N1 ϕ(N),1 ϕ(N),Nϕ(N) respectively, where

s j, j′ −1 ( j) s j, j′ −1 ∂ σ(−ω k j, j′ ) ∂ σ(k j, j′ ) t ( j) k j, j′ N t ′ k j, j′ s j, j′ = σ(−ωN k j, j′ ), ··· ,  , s j, j′ = σ(k j, j′ ), ··· ,  .  (s j, j′ − 1)!   (s j, j′ − 1)!          The M is a block Cauchy matrix with its blocks M (α,β) deﬁned by   M(α,β) ··· M(α,β) ··· M(α,β) ··· M(α,β) 1,1;1,1 1,1;1,N1 1,1;ϕ(N),1 1,1;ϕ(N),Nϕ(N)  . . . .   . . . .   (α,β) (α,β) (α,β) (α,β)   M ··· M ··· M ··· M   1,N1;1,1 1,N1;1,N1 1,N1;ϕ(N),1 1,N1;ϕ(N),Nϕ(N)   . . . .  M (α,β) =  . . . .  , (6.28)  . . . .   (α,β) (α,β) (α,β) (α,β)   M ··· M ··· M ··· M   ϕ(N),1;1,1 ϕ(N),1;1,N1 ϕ(N),1;ϕ(N),1 ϕ(N),1;ϕ(N),Nϕ(N)   . . . .   . . . .    M(α,β) ··· M(α,β) ··· M(α,β) ··· M(α,β)   ϕ(N),Nϕ(N);1,1 ϕ(N),Nϕ(N);1,N1 ϕ(N),Nϕ(N);ϕ(N),1 ϕ(N),Nϕ(N);ϕ(N),Nϕ(N)     (1,1)  for α, β ∈ {1, 2}. Here the M j, j′;i,i′ takes the form of (1,1) ( j) ( j) ′ ′ ′ ′ M j, j′;i,i′ = S(−ωN k j, j )G(ki,i , −ωN k j, j )R(ki,i ), (6.29) in which

( j) si,i′ −1 σ(−ω k j, j′ ) ∂ ρ(ki,i′ ) N ki,i′ ··· ρ(k ′ ) . . (s ′ −1)! i,i ( j)  . ..   i,i  S(−ω k j, j′ ) =  .  , R(ki,i′ ) =  . .  , (6.30) N  (s ′ −1)   . .   j, j ( j)   . .   ∂ σ(−ω k j, j′ )     k j, j′ N ( j)     ··· − ′   ′   ′ σ( ω k j, j )  ρ(ki,i )   (s j, j −1)! N        and    

′ 0 ′ 0 k j, j Cs −1 k j, j C0 i,i′ ( j) ··· ( j) s ′ s ′ −1 k ′ −ω k ′ (k ′ −ω k ′ ) i,i (−1) i,i  i,i N j, j i,i N j, j  ( j)  . .  G(ki,i′ , −ωN k j, j′ ) = (s j, j′ − 1)!  . .   s −1 s −2   ( j) s j, j′ −1 j, j′ ( j) s j, j′ −2 j, j′   k ′ (−ω ) C (−ω ) C   j, j N s ′ −1 N s ′ −2   j, j j, j   ( j) s ′ s ′ −1 + ( j) s ′ −1 s ′ −2 · · · ♠   ′ ′ j, j j, j ′ ′ j, j j, j   (ki,i −ωN k j, j ) (−1) (ki,i −ωN k j, j ) (−1)    for   ′ ′ ( j) s j, j −1 s j, j′ −1 ( j) s j, j −2 s j, j′ −2 k ′ −ω C −ω C j, j ( N ) s ′ +s ′ −2 ( N ) s ′ +s ′ −3 ♠ = i,i j, j + i,i j, j ; ( j) si,i′ +s j, j′ −1 si,i′ +s j, j′ −2 ( j) si,i′ +s j, j′ −2 si,i′ +s j, j′ −3 (ki,i′ − ωN k j, j′ ) (−1) (ki,i′ − ωN k j, j′ ) (−1) (1,2) the M j, j′;i,i′ is deﬁned as (1,2) ( j) ( j) ( j) ( j) M j, j′;i,i′ = S(−ωN k j, j′ )G(−ωN ki,i′ , −ωN k j, j′ )R(−ωN ki,i′ ) (6.31) in which

si,i′ −1 ( j) ( j) − ′ − ′ ∂ ( j) ρ( ωN ki,i ) σ( ω k j, j ) −ω k N N i,i′ ( j) ··· − ′ . . ′ − ρ( ωN ki,i ) ( j)  . ..  ( j)  (si,i 1)!  S(−ωN k j, j′ ) =   , R(−ωN ki,i′ ) =  . .  , (6.32)  (s j, j′ −1) ( j)   . ..   ∂ σ(−ω k j, j′ )   .   k j, j′ N ( j)     ··· σ(−ω k j, j′ )  ( j)   (s j, j′ −1)! N   ρ(−ω ki,i′ )     N      LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 21 and

′ 0 ′ 0 k j, j Cs −1 k j, j C0 i,i′ ( j) ( j) ··· ( j) ( j) s ′ s ′ −1 −ω k ′ −ω k ′ (−ω k ′ −ω k ′ ) i,i (−1) i,i  N i,i N j, j N i,i N j, j  ( j) ( j)  . .  G(−ωN ki,i′ , −ωN k j, j′ ) = (s j, j′ − 1)!  . .   s −1 s −2   ( j) s j, j′ −1 j, j′ ( j) s j, j′ −2 j, j′   k ′ (−ω ) C (−ω ) C   j, j N s ′ −1 N s ′ −2   j, j j, j   ( j) ( j) s ′ s ′ −1 + ( j) ( j) s ′ −1 s ′ −2 ··· ♥   (−ω k ′ −ω k ′ ) j, j (−1) j, j (−ω k ′ −ω k ′ ) j, j (−1) j, j   N i,i N j, j N i,i N j, j    for  

′ ′ ( j) s j, j −1 s j, j′ −1 ( j) s j, j −2 s j, j′ −2 ′ − − k j, j ( ωN ) Cs ′ +s ′ −2 ( ωN ) Cs ′ +s ′ −3 ♥ = i,i j, j + i,i j, j ; ( j) ( j) si,i′ +s j, j′ −1 si,i′ +s j, j′ −2 ( j) ( j) si,i′ +s j, j′ −2 si,i′ +s j, j′ −3 (−ωN ki,i′ − ωN k j, j′ ) (−1) (−ωN ki,i′ − ωN k j, j′ ) (−1) (2,1) the M j, j′;i,i′ is expressed by (2,1) M j, j′;i,i′ = S(k j, j′ )G(ki,i′ , k j, j′ )R(ki,i′ ), (6.33) where

s ′ −1 ′ i,i σ(k j, j ) ∂ ρ(ki,i′ ) ki,i′ ··· ′ . . ′ ρ(ki,i )  . ..   (si,i −1)!  S(k j, j′ ) = . , R(ki,i′ ) = . , (6.34)  (s −1)   .   j, j′   . ..   ∂ σ(k j, j′ )   .   k j, j′     ··· σ(k ′ )  ′   (s ′ −1)! j, j   ρ(ki,i )   j, j        and    

( j) 0 ( j) 0 − ′ −ω k ′ C ωN k j, j Cs ′ −1 N j, j 0 ··· i,i k ′ +k ′ si,i′ si,i′ −1  i,i j, j (ki,i′ +k j, j′ ) (−1)   . .  G(k ′ , k ′ ) = (s ′ − 1)!  . .  i,i j, j j, j  . .   ( j) s j, j′ −1 ( j) s j, j′ −2   −ω k ′ C ω C   N j, j s ′ −1 N s ′ −2   j, j + j, j ··· ♣   s j, j′ s j, j′ −1 s j, j′ −1 s j, j′ −1   (ki,i′ +k j, j′ ) (−1) (ki,i′ +k j, j′ ) (−1)    for  

( j) s j, j′ −1 ( j) s j, j′ −2 − ′ ωN k j, j Cs ′ +s ′ −2 ωN Cs ′ +s ′ −3 ♣ = i,i j, j + i,i j, j ; s ′ +s ′ −1 si,i′ +s j, j′ −2 s ′ +s ′ −2 si,i′ +s j, j′ −3 (ki,i′ + k j, j′ ) i,i j, j (−1) (ki,i′ + k j, j′ ) i,i j, j (−1)

(2,2) and ﬁnally M j, j′;i,i′ is given as

(2,2) ( j) ( j) M j, j′;i,i′ = S(k j, j′ )G(−ωN ki,i′ , k j, j′ )R(−ωN ki,i′ ), (6.35) where

si,i′ −1 ( j) ′ σ k ′ ∂ ( j) ρ(−ωN ki,i ) ( j, j ) −ω k N i,i′ ( j) . . ··· ρ(−ω ki,i′ )  . ..  ( j)  (si,i′ −1)! N  S(k j, j′ ) = . , R(−ω ki,i′ ) =   , (6.36)  (s −1)  N  . .   j, j′   . .   ∂ σ(k j, j′ )   . .   k j, j′     ′   ( j)   ··· σ(k j, j )  ′   (s j, j′ −1)!   ρ(−ω ki,i )     N      and  

( j) ( j) ′ 0 ′ 0 −ωN k j, j C − −ωN k j, j C0 si,i′ 1 ( j) ··· ( j) s ′ s ′ −1 −ω k ′ +k ′ (−ω k ′ +k ′ ) i,i (−1) i,i  N i,i j, j N i,i j, j  ( j)  . .  G(−ω k ′ , k ′ ) = (s ′ − 1)!  . .  N i,i j, j j, j  . .   ( j) s j, j′ −1 ( j) s j, j′ −2   −ω k j, j′ C ω C   N s j, j′ −1 N s j, j′ −2   ( j) s ′ s ′ −1 + ( j) s ′ −1 s ′ −1 ··· ♦   (−ω k ′ +k ′ ) j, j (−1) j, j (−ω k ′ +k ′ ) j, j (−1) j, j   N i,i j, j N i,i j, j    for  

( j) s j, j′ −1 ( j) s j, j′ −2 − ′ ωN k j, j Cs ′ +s ′ −2 ωN Cs ′ +s ′ −3 ♦ = i,i j, j + i,i j, j . ( j) si,i′ +s j, j′ −1 si,i′ +s j, j′ −2 ( j) si,i′ +s j, j′ −2 si,i′ +s j, j′ −3 (−ωN ki,i′ + k j, j′ ) (−1) (−ωN ki,i′ + k j, j′ ) (−1)

(2) (1) (2) (1) 6.6. A2r - and Cr -types. The linear integral equations and the potentials in the A2r and Cr classes have the uniﬁed forms (5.24) and (5.26), respectively. Taking the same measure as (6.2) and performing this on (5.24) and (5.26), we are able to construct the Cauchy matrix formula of U (¯i,¯j), by following a similar derivation to that in the previous cases, which takes the form of

(¯i) (1,1) (1,2) −1 ′(¯j) ¯ ¯ K 0 A 0 M M A 0 D 0 K 0 s U (i, j) = (tr, tr′) I + . (6.37) 0 K′′(¯i)! " 0 A′! M (2,1) M (2,2)!# 0 A′! 0 D! 0 K(¯j)! s′ ! 22 YUE YIN AND WEI FU

Here A, A′, D, r, r′, s as well as s′ are exactly the same as those given in (6.24), (6.25), (6.26) and (6.27). The K(¯i), K′(¯i) and K′′(¯i) are given as follows: ¯ ¯ ¯ ¯ K(¯i) = diag(K(i) , ··· , K(i) ; ··· ; K(i) , ··· , K(i) ), 1,1 1,N1 ϕ(N),1 ϕ(N),Nϕ(N) ¯ ¯ ¯ ¯ K′(¯i) = diag(K′(i), ··· , K′(i) ; ··· ; K′(i) , ··· , K′(i) ), 1,1 1,N1 ϕ(N),1 ϕ(N),Nϕ(N) ¯ ¯ ¯ ¯ K′′(¯i) = diag(K′′(i), ··· , K′′(i) ; ··· ; K′′(i) , ··· , K′′(i) ), 1,1 1,N1 ϕ(N),1 ϕ(N),Nϕ(N) where

¯i ( j) ¯i k ′ j, j (−ωN k j, j′ ) . .  . .  (¯i)  . ..  ′(¯i) . . K ′ =  .  , K ′ =  . .  , j, j  s ′ −1  j, j    j, j ¯i   s j, j′ −1 ( j) ¯i   ∂ k   ′   k ′ j, j′   ∂ (−ωN k j, j ) ¯i  j, j ¯i   k j, j′ ( j)   ··· k ′   ··· (−ω k ′ )   (s j, j′ −1)! j, j   (s ′ −1)! N j, j     j, j      and    

( j) ¯i (−ωN k j, j′ )  . .  ′′(¯i)  . ..  K j, j′ =   .  s j, j′ −1 ( j) ¯i   ′   ∂ ( j) (−ωN k j, j )   −ω k ¯i  N j, j′ ( j)   ··· −ω k ′   (s ′ −1)! ( N j, j )   j, j    (α,β)  (α,β)  The Cauchy matrices M take their respect forms of (6.28), in which M j, j′;i,i′ are deﬁned in (6.29), (6.31), (6.33) and (6.35), respectively, with the same R and S matrices listed in (6.30), (6.32), (6.34) and (6.36), but diﬀerent G matrices as follows:

0 0 Cs −1 C0 i,i′ ( j) ··· ( j) s ′ s ′ −1 k ′ −ω k ′ (k ′ −ω k ′ ) i,i (−1) i,i  i,i N j, j i,i N j, j  ( j)  . .  G(k ′ , −ω k ′ ) = (s ′ − 1)!  . .  , i,i N j, j j, j  . .   s ′ −1 s , ′ −1 s ′ −1 s , ′ −1   (−ω( j)) j, j C j j (−ω( j)) j, j C j j   N s ′ −1 N s ′ +s ′ −2   j, j i,i j, j   ( j) s ′ s ′ −1 ··· ( j) s ′ +s ′ −1 s ′ +s ′ −2   (k ′ −ω k ′ ) j, j (−1) j, j (k ′ −ω k ′ ) i,i j, j (−1) i,i j, j   i,i N j, j i,i N j, j   0   0 Cs −1  C0 i,i′ ( j) ( j) ··· ( j) ( j) s ′ s ′ −1 −ω k ′ −ω k ′ (−ω k ′ −ω k ′ ) i,i (−1) i,i  N i,i N j, j N i,i N j, j  ( j) ( j)  . .  G(−ω k ′ , −ω k ′ ) = (s ′ − 1)!  . .  , N i,i N j, j j, j  . .   ( j) s ′ −1 s j, j′ −1 ( j) s ′ −1 s j, j′ −1   (−ω ) j, j C (−ω ) j, j C   N s ′ −1 N s ′ +s ′ −2   j, j i,i j, j   ( j) ( j) s ′ s ′ −1 ··· ( j) ( j) s ′ +s ′ −1 s ′ +s ′ −2   (−ω k ′ −ω k ′ ) j, j (−1) j, j (−ω k ′ −ω k ′ ) i,i j, j (−1) i,i j, j   N i,i N j, j N i,i N j, j   0   0 C   C s ′ −1  0 ··· i,i k ′ +k ′ si,i′ si,i′ −1  i,i j, j (ki,i′ +k j, j′ ) (−1)  ′  . .  ′ ′ = ′ −  . .  G (ki,i , k j, j ) (s j, j 1)!  . .  ,  s j, j′ −1 s j, j′ −1   C C   s ′ −1 s ′ +s ′ −2   j, j i,i j, j   s ′ s ′ −1 ··· s ′ +s ′ −1 s ′ +s ′ −2   (k ′ +k ′ ) j, j (−1) j, j (k ′ +k ′ ) i,i j, j (−1) i,i j, j   i,i j, j i,i j, j   0   0 Cs −1  C0 i,i′ ( j) ··· ( j) s ′ s ′ −1 −ω k ′ +k ′ (−ω k ′ +k ′ ) i,i (−1) i,i  N i,i j, j N i,i j, j  ( j)  . .  G(−ω k ′ , k ′ ) = (s ′ − 1)!  . .  . N i,i j, j j, j  . .   s , ′ −1 s , ′ −1   C j j C j j   s ′ −1 s ′ +s ′ −2   j, j i,i j, j   ( j) s ′ s ′ −1 ··· ( j) s ′ +s ′ −1 s ′ +s ′ −2   (−ω k ′ +k ′ ) j, j (−1) j, j (−ω k ′ +k ′ ) i,i j, j (−1) i,i j, j   N i,i j, j N i,i j, j    Then the expression ϕ = ln 1 − U (0,−1) provides the Cauchy matrix solutions to the corresponding 2D Toda-type equations. The (1) (2) N = 2r and N = 2r + 1 cases for positive integers r correspond to the 2D Toda systems of Cr - and A2r -types, respectively. 6.7. Examples. We have presented all the general formulae for Cauchy matrix solutions to the 2DTLs. Below we take the three scalar equations, namely the A∞-type 2D Toda equation, the sinh–Gordon equation and the Tzitzeica equation as examples, to show how soliton solution, multi-pole solution and their mixture are obtained as special cases of the Cauchy matrix solutions.

Example 6.1. The first example is the 2D Toda equation of A∞-type. We consider a special case of (6.2) as follows: 2 3 ′ 1 1 1 ′ ζ k, k = A ′ k k , d ( ) j, j 2 ′ ′ d d (2πi) k − k k − k ′ Xj=1 Xj=1 j j ′ ′ ′ namely the double integration has only first-order singularities k j for j = 1, 2 and k j′ for j = 1, 2, 3 on the k- and k - planes, ′ (0,−1) located in the contours Γ and Γ , respectively. This results in the expression of Un given by (0,−1) t −1 ′(−1) Un = r(I + AM) AK s, LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 23 where the matrices A, M and K′(−1) and the vectors r and s are defined as

′ ′ ρn(k1)σn(k1) ρn(k2)σn(k1) ′ ′ k1+k1 k2+k1 ρ (k )σ (k′ ) ρ (k )σ (k′ ) A1,1 A1,2 A1,3  n 1 n 2 n 2 n 2  ′(−1) ′−1 ′−1 ′−1 A = , M = k +k′ k +k′ , K = diag k 1 , k 2 , k 3 , A A A  1 2 2 2  2,1 2,2 2,3!  ρ (k )σ (k′ ) ρ (k )σ (k′ )   n 1 n 3 n 2 n 3   k +k′ k +k′   1 3 2 3    and   t t r = (ρn(k1), ρn(k2)) , as well as s = (σn(k1), σn(k2), σn(k3)) , ′ respectively, with ρn(k) and σn(k ) deﬁned by (3.1). We refer to the corresponding ϕn deﬁned in (6.1) as the (2, 3)-soliton solution for (1.1).

(1) (1) Example 6.2. In the A1 class, since there is only one primitive square root of unit ω2 = −1, we consider a special measure 1 1 = dλ1(k) A1,1 2 dk, 2πi (k − k1,1) namely, we let the contour Γ1 on the k-plane contain only a double-pole k1,1 inside. Then the general formula (6.21) provides (0,−1) t −1 ′(−1) U0 = r(I + AM) AK s, in which the matrices A, M and K(−1) are given as follows:

1 1 −1 − 2 k 0 A1,1 0 σ0(k1,1) 0 k1,1+k1,1 (k1,1+k1,1) ∂k1,1 ρ0(k1,1) ρ0(k1,1) ′(−1) 1,1 A = , M = 1 2 , K = −2 −1 , 0 A1,1 ∂k1,1 σ0(k1,1) σ0(k1,1) − 2 3  ρ0(k1,1) 0 −k k ! ! (k1,1+k1,1) (k1,1+k1,1) ! 1,1 1,1!   and the vectors r and s are deﬁned as  

t t r = ∂k1,1 ρ0(k1,1), ρ0(k1,1) and s = σ0(k1,1),∂k1,1 σ0(k1,1) , (0,−1) respectively. Thus, ϕ0 = ln 1 − U0 provides a double-pole solution to the sinh–Gordon equation (1.3). (2) Example 6.3. The third example is the Tzitzeica equation (1.4b), as the ﬁrst nontrivial example in the A2r class, namely the (1) (2) 2 r = 1 case. In this case, there are two primitive cube roots of unity ω3 = ω3 and ω3 = ω3, where ω3 = exp(2πi/3). We consider measures 1 1 1 1 λ k = A k and λ k = A k, d 1( ) 1,1 d d 2( ) 2,1 2 d 2πi k − k1,1 2πi (k − k2,1) as particular cases of (6.16), which implies that there are two singularities on the k-plane, including the ﬁrst-order singularity k1,1 located in the contour Γ1 and the second-order singularity k2,1 located in Γ2. Then equation (6.37) implies that the Tzitzeica (0,−1) (0,−1) equation has its Cauchy matrix solution ϕ0 = ln 1 − U0 with the expression of U0 given by − A 0 M (1,1) M (1,2) 1 A 0 K′(−1) 0 s U (0,−1) = (tr, tr′) I + . 0 " 0 A′! M (2,1) M (2,2)!# 0 A′! 0 K(−1)! s′ ! Here A and A′ are the 3 × 3 diagonal matrices

′ (2) A = Diag A1,1, A2,1, A2,1 and A = Diag A1,1, −ω3 A2,1, A2,1 . The 3-component vectors r, r′, s and s′are deﬁned by

t t ′ (1) (2) (2) (2) r = ρ0(k1,1),∂k2,1 ρ0(k2,1), ρ0(k2,1) , r = ρ0(−ω k1,1),∂ ρ0(−ω k2,1), ρ0(−ω k2,1) , 3 −ω3 k2,1 3 3 and t = − (1) − (2) − (2) t ′ = s σ0( ω3 k1,1), σ0( ω3 k2,1),∂k2,1 σ0( ω3 k2,1) , s σ0(k1,1), σ0(k2,1),∂k2,1 σ0(k2,1) , respectively. K(−1) and K′(−1) are 3 × 3 lower-triangle matrices −1 − (2) −1 k1,1 0 0 ω3 k1,1 0 0 K(−1) = 0 k−1 0 and K′(−1) =  0 −ω(1)k−1 0  ,  2,1   3 2,1   0 −k−2 k−1   0 ω(1)k−2 −ω(1)k−1   2,1 2,1  3 2,1 3 2,1     (−1)  ′(−1)    respectively. In either K or K , we observe that the upper-left corner is a scalar, corresponding to the singularity k1,1 of order 1; while the lower-right corner is a 2 × 2 Jordan matrix, corresponding to the singularity k2,1 of order 2. In other words, the two singularities results in the Cauchy matrix solution being the mixture of a simple-pole soliton and a double-pole soliton. M (α,β) for α, β ∈ {1, 2} in the Cauchy matrix is of the form M (α,β) M (α,β) M (α,β) = 1,1;1,1 1,1;2,1 , M (α,β) M (α,β)   2,1;1,1 2,1;2,1   24 YUE YIN AND WEI FU

(α,β) (α,β) (α,β) (α,β) in which M1,1;1,1 is a scalar, M1,1;2,1 and M2,1;1,1 are 2-component row and column vectors, respectively, and M2,1;2,1 is a 2 × 2 matrix. When α = 1 and β = 1, we have for M (1,1) the following blocks: (1) σ0(−ω k1,1)ρ0(k1,1) M (1,1) = 3 , 1,1;1,1 (1) k1,1 − ω3 k1,1 − (1,1) (1) 1 1 ∂k2,1 ρ0(k2,1) ρ0(k2,1) M = σ0(−ω k1,1) , , 1,1;2,1 3  (1) (1) 2  ρ (k ) 0 k2,1 − ω k1,1 (k2,1 − ω k1,1) 0 2,1 !  3 3   1  (2)   k −ω(2)k (1,1) σ0(−ω3 k2,1) 0 1,1 3 2,1 M = (2) (2)  ω(2)  ρ0(k1,1), 2,1;1,1 ∂ σ (−ω k ) σ (−ω k ) 3 k2,1 0 3 2,1 0 3 2,1 !  (2) 2   (k1,1−ω k2,1)   3   1  −1 (2) − (2) − (2) 2 (1,1) σ0(−ω k2,1) 0 k2,1 ω3 k2,1 (k2,1 ω3 k2,1) ∂k ρ0(k2,1) ρ0(k2,1) 3 (2) (2) 2,1 M2,1;2,1 = (2) (2)  ω −2ω  . ∂ σ (−ω k ) σ (−ω k ) 3 3 ρ0(k2,1) 0 k2,1 0 3 2,1 0 3 2,1 !  (2) 2 (2) 3  !  (k2,1−ω k2,1) (k2,1−ω k2,1)   3 3    When α = 1 and β = 2, the corresponding blocks in M (1,2) are as follows:  (1) (1) σ0(−ω k1,1)ρ0(−ω k1,1) M (1,2) = 3 3 , 1,1;1,1 (1) (1) −ω3 k1,1 − ω3 k1,1 (1) (1) 1 −1 ∂ (1) ρ (−ω k ) ρ (−ω k ) (1,2) (1) −ω k2,1 0 3 2,1 0 3 2,1 M = σ0(−ω k1,1) , 3 , 1,1;2,1 3  (1) (1) (1) (1) 2   (1)  −ω k2,1 − ω k1,1 (−ω k2,1 − ω k1,1) ρ0(−ω3 k2,1) 0  3 3 3 3     1    (2)     −ω(2)k −ω(2)k (1,2) σ0(−ω3 k2,1) 0 3 1,1 3 2,1 (2) M = (2) (2)  ω(2)  ρ0(−ω k1,1), 2,1;1,1 ∂ σ (−ω k ) σ (−ω k ) 3 3 k2,1 0 3 2,1 0 3 2,1 !  (2) (2) 2   (−ω k1,1−ω k2,1)   3 3   1  −1 (2)  (2) (2)  (2) (2) (2) (2) −ω k −ω k −ω k −ω k 2 (2) (1,2) σ0(−ω k2,1) 0 3 2,1 3 2,1 ( 3 2,1 3 2,1) ∂−ω k ρ0(−ω k2,1) ρ0(−ω k2,1) 3 (2) (2) 3 2,1 3 3 M2,1;2,1 = (2) (2)  ω −2ω  (2) . ∂k σ0(−ω k2,1) σ0(−ω k2,1)! 3 3  −  2,1 3 3  (2) (2) 2 (2) (2) 3  ρ0( ω3 k2,1) 0  (−ω k2,1−ω k2,1) (−ω k2,1−ω k2,1)     3 3 3 3        For α = 2 and β = 1, M (2,1) takes its corresponding blocks as follows: 

(2,1) σ0(k1,1)ρ0(k1,1) M1,1;1,1 = , k1,1 + k1,1 1 −1 ∂ ρ (k ) ρ (k ) M (2,1) = σ k , k2,1 0 2,1 0 2,1 , 1,1;2,1 0( 1,1) 2 k2,1 + k1,1 (k2,1 + k1,1) ! ρ0(k2,1) 0 ! 1 (2,1) σ0(k2,1) 0 k1,1+k2,1 M2,1;1,1 = −1 ρ0(k1,1), ∂k2,1 σ0(k2,1) σ0(k2,1)!  (k +k )2   1,1 2,1   1  −1   2 (2,1) σ0(k2,1) 0 k2,1+k2,1 (k2,1+k2,1) ∂k2,1 ρ0(k2,1) ρ0(k2,1) M2,1;2,1 = −1 2 . ∂k2,1 σ0(k2,1) σ0(k2,1)!  2 3  ρ0(k2,1) 0 !  (k2,1+k2,1) (k2,1+k2,1)    For α = 2 and β = 1, M (2,2) is composed of the following blocks:  (1) σ0(k1,1)ρ0(−ω k1,1) M (2,2) = 3 1,1;1,1 (1) , −ω3 k1,1 + k1,1 (1) (1) 1 −1 ∂ (1) ρ (−ω k ) ρ (−ω k ) (2,2) −ω k2,1 0 3 2,1 0 3 2,1 M = σ0(k1,1) , 3 , 1,1;2,1  (1) (1) 2   (1)  −ω k2,1 + k1,1 (−ω k2,1 + k1,1) ρ0(−ω k2,1) 0  3 3   3   1    σ (2) (2,2) 0(k2,1) 0 −ω3 k1,1+k2,1 (2) M2,1;1,1 = −1 ρ0(−ω3 k1,1), ∂k σ0(k2,1) σ0(k2,1)  (2)  2,1 ! (−ω k +k )2  3 1,1 2,1     1  −1 (2) (2) (2,2) σ0(k2,1) 0 −ω(2)k +k (−ω(2)k +k )2 ∂ (2) ρ0(−ω k2,1) ρ0(−ω k2,1) 3 2,1 2,1 3 2,1 2,1 −ω3 k2,1 3 3 M2,1;2,1 = −1 2 (2) . ∂k σ0(k2,1) σ0(k2,1)     2,1 ! (−ω(2)k +k )2 (−ω(2)k +k )3 ρ0(−ω3 k2,1) 0  3 2,1 2,1 3 2,1 2,1        (α,β) (α,β)    We comment that the blocks M1,1;1,1 and M2,1;2,1 describe the behaviour of the simple-pole soliton and the double-pole soliton (α,β) (α,β) themselves; while the blocks M1,1;2,1 and M2,1;1,1 explain the interaction between the two solitons.

7. Concluding remarks

We established the DL schemes for the 2D Toda-type equations associated with the infinite-dimensional algebras A∞, B∞ and (1) (2) (1) (2) C∞, as well as the Kac–Moody algebras Ar , A2r , Cr and Dr+1. This was realised by associating the 2DTLs with the linear integral equations taking the form of (2.5). The starting point was the linear integral equation (2.5a) and its adjoint (2.5b) with arbitrary integration measure and domain, whose nonlinearisation resulted in the A∞-type 2DTL; while the 2D Toda equations of other types were obtained by imposing various constraints on the measure and the domain in the A∞ class. The linear integral LINEAR INTEGRAL EQUATIONS AND TWO-DIMENSIONAL TODA SYSTEMS 25 equations also helped to rediscover the Lax scheme (including both Lax pair and its adjoint) and the bilinear formulation in the 2D Toda theory, given by Fordy and Gibbons [5,6], Nimmo and Willox [25], as well as the Kyoto school [12,32]. The key point in the construction of the 2D Toda-type equations is the introduction of the infinite matrix U. This allows us to transfer the discrete and continuous dynamical relations to purely algebraic ones. Then the problem of constructing integrable systems turns to be searching for identities of the entries of the infinite matrix. In other words, the DL approach is a method to construct nonlinear equations directly. There are also other methods to construct the 2D Toda systems such as the Lax scheme and the bilinear approach, but the motivation seems different compared with the DL, to the best of the authors’ knowledge. For example, the starting point in the fermionic approach (i.e. the bilinear method) is a bilinear identity, which can generate bilinear 2D Toda equationsas the main object, see [12,32]. As for the Lax scheme for the 2DTLs (see e.g. [6]), the key step is to construct linear spectral problems corresponding to various Kac–Moody algebras, while the nonlinear equations arise as the compatibility conditions of the corresponing linear problems. The linear integral equations theoretically provide a very large class of solutions for nonlinear integral equations. This is because the formal solution obtained from the scheme, namely (2.7), is an integral defined on a very general integration domain in the complex space in terms of the spectral parameter(s), compared with other methods. Our main achievement in this paper is that we successfully constructed direct linearising type solution for each class of the 2D Toda-type equation. As a degenerate case of the direct linearising type solution, we presented the general formula of the Cauchy matrix type solution for each class of 2DTL, by imposing the condition that the integration is defined on the domain containing an arbitrary (1) number of singularities of arbitrary order. The solutions of A∞- and A1 -types coincide with those obtained from the operator approach and the so-called generalised Cauchy matrix approach, see e.g. [24,29,30] and also [19,34,35]. Here we extended these results to nonlinear integrable systems associated with other Lie algebras. The procedure of deriving Cauchy matrix type solutions can actually reduce the formal tau functions in the form of (2.11) to the corresponding degenerate cases, resulting in determinant type solutions for the bilinear forms of the 2D Toda systems. These solutions are closely related to the so-called Grammian solutions obtained from the well-known Hirota method, see e.g. [11]. We also note that DL in some cases induces more general solutions, compared with Grammian solutions. For instance, the numbers ′ of the singularities on the k- and k -plane are not necessarily equal in our result for the A∞-type 2DTL. This simultaneously results in the existence of the finite matrix A in (6.7), see also example 6.1. In fact, the coefficient matrix A plays a crucial role in the classification of solitons, explaining a variety of soliton behaviours, see e.g. [2, 13–15]. The DL scheme of the 2DTL of D∞-type and its reduced systems, as well as the equations associated with the exceptional Lie algebras, are not covered by the scheme in this paper. Searching for other classes of exact solutions for the 2D Toda-type equations from the linear integral equations is also an interesting problem, but it then requires nontrivial spectral analysis for the linear integral equations. This remains for our future work.

Acknowledgments This project was supported by the National Natural Science Foundation of China (grant no. 11901198) and Shanghai Pujiang Program (grant no. 19PJ1403200). WF was also partially sponsored by the Science and Technology Commission of Shanghai Municipality (grant no. 18dz2271000) as well as the Fundamental Research Funds for the Central Universities.

References [1] M.J. Ablowitz and P.A. Clarkson. Solitons, Nonlinear Evolution Equations and Inverse Scattering. Cambridge University Press, Cambridge, 1991. [2] R. Bury, A.V. Mikhailov, and J.P. Wang. Wave fronts and cascades of soliton interactions in the periodic two dimensional Volterra system. Physica D, 347:21–41, 2017. [3] A.S. Fokas and M.J. Ablowitz. Linearization of the Korteweg–de Vries and PainlevéII equations. Phys. Rev. Lett., 47:1096–1110, 1981. [4] A.S. Fokas and M.J. Ablowitz. On the inverse scattering and direct linearizing transforms for the Kadomtsev–Petviashvili equation. Phys. Lett. A, 94:67–70, 1983. [5] A.P. Fordy and J. Gibbons. Integrable nonlinear Klein–Gordon equations and Toda lattices. Commun. Math. Phys., 77:21–30, 1980. [6] A.P. Fordy and J. Gibbons. Nonlinear Klein–Gordon equations and simple Lie algebras. Proc. R. Ir. Acad. A, 83:33–44, 1983. [7] W. Fu. Direct linearization approach to discrete integrable systems associated with ZN graded Lax pairs. Proc. R. Soc. A, 476:20200036, 2020. [8] W. Fu and F.W. Nijhoff. Direct linearizing transform for three-dimensional discrete integrable systems: the lattice AKP, BKP and CKP equations. Proc. R. Soc. A, 473:20160915, 2017. [9] W. Fu and F.W. Nijhoff. Linear integral equations, infinite matrices, and soliton hierarchies. J. Math. Phys., 59:071101, 2018. [10] W. Fu and F.W. Nijhoff. On nonautonomous differential-difference AKP, BKP and CKP equations. Proc. R. Soc. A, 477:20200717, 2021. [11] R. Hirota. The Direct Method in Soliton Theory. Cambridge University Press, Cambridge, 2004. [12] M. Jimbo and T. Miwa. Solitons and infinite dimensional Lie algebras. Publ. RIMS, 19:943–1001, 1983. [13] Y. Kodama and L.K. Williams. KP solitons, total positivity, and cluster algebras. Proc. Natl. Acad. Sci., 108:8984–8989, 2011. [14] Y. Kodama and L.K. Williams. The Deodhar decomposition of the Grassmannian and the regularity of KP solitons. Adv. Math., 244:979–1032, 2013. [15] Y. Kodama and L.K. Williams. KP solitons and total positivity for the Grassmannian. Invent. Math., 198:637–699, 2014. [16] V.B. Matveev and M.A. Salle. Darboux Transformations and solitons. Springer-Verlag, Berlin, 1991. [17] A.V. Mikhailov. Integrability of a two-dimensional generalization of the Toda chain. JETP Lett., 30:414–418, 1979. [18] A.V. Mikhailov, M.A. Olshanetsky, and A.M. Perelomov. Two-dimensional generalized Toda lattice. Commun. Math. Phys., 79:473–488, 1983. [19] F.W. Nijhoff, J. Atkinson, and J. Hietarinta. Soliton solutions for ABS lattice equations: I. Cauchy matrix approach. J. Phys. A: Math. Theor., 42:404005, 2009. [20] F.W. Nijhoff, H.W. Capel, and G.L. Wiersma. Integrable lattice systems in two and three dimensions. In R. Martini, editor, Geometric Aspects of the Einstein Equations and Integrable Systems, volume 239 of Lecture Notes in Physics, pages 263–302. Springer-Verlag, Berlin, 1985. [21] F.W. Nijhoff, A. Hone, and N. Joshi. On a Schwarzian PDE associated with the KdV hierarchy. Phys. Lett. A, 267:147–156, 2000. [22] F.W. Nijhoff, V.G. Papageorgiou, H.W. Capel, and G.R.W. Quispel. The lattice Gel’fand–Dikii hierarchy. Inverse Probl., 8:597–621, 1992. [23] F.W. Nijhoff, A. Ramani, B. Grammaticos, and Y. Ohta. On Discrete Painlevéequations associated with the lattice KdV systems and the PainlevéVI equation. Stud. Appl. Math., 106:261–314, 2001. 26 YUE YIN AND WEI FU

[24] T. Nilson and C. Schiebold. Solution formulas for the two-dimensional Toda lattice and particle-like solutions with unexpected asymptotic behaviour. J. Nonlinear Math. Phys., 27:57–94, 2020. [25] J.J.C. Nimmo and R. Willox. Darboux transformations for the two-dimensional Toda system. Proc. R. Soc. A, 453:2497–2525, 1997. [26] S.P. Novikov, S.V. Manakov, L.P. Pitaevskii, and V.E. Zakharov. Theory of Solitons: The Inverse Scattering Method. Consultants Bureau, New York, 1984. [27] G.R.W. Quispel, F.W. Nijhoff, H.W. Capel, and J. van der Linden. Linear integral equations and nonlinear difference-difference equations. Physica A, 125:344–380, 1984. [28] P.M. Santini, M.J. Ablowitz, and A.S. Fokas. The direct linearization of a class of nonlinear evolution equations. J. Math. Phys., 25:2614–2619, 1984. [29] C. Schiebold. From the non-abelian to the scalar two-dimensional Toda lattice. Glasgow Math. J., 47A:177–189, 2005. [30] C. Schiebold. Cauchy-type determinants and integrable systems. Linear Algebra Appl., 433:447–475, 2010. [31] M. Toda. Vibration of a chain with nonlinear interaction. J. Phys. Soc. Jpn., 22:431–436, 1967. [32] K. Ueno and K. Takasaki. Toda lattice hierarchy. In K. Okamoto, editor, Group Representations and Systems of Differential Equations, volume 4 of Advanced Studies in Pure Mathematics, pages 1–95. Kinokuniya, Tokyo, 1984. [33] G. Wilson. The modified Lax and two-dimensional Toda lattice equations associated with simple Lie algebras. Ergod. Theory Dyn. Syst., 1:361–380, 1981. [34] D.D. Xu, D.J. Zhang, and S.L. Zhao. The Sylvester equation and integrable equations: I. The Kortewegde–de Vries system and sine–Gordon equation. J. Nonlinear Math. Phys., 21:382–406, 2014. [35] D.J. Zhang and S.L. Zhao. Solutions to ABS lattice equations via generalized Cauchy matrix approach. Stud. Appl. Math., 131:72–103, 2013. [36] D.J. Zhang, S.L. Zhao, and F.W. Nijhoff. Direct linearisation of extended lattice BSQ systems. Stud. Appl. Math., 129:220–248, 2012.

(YY) School of Mathematical Sciences,East China Normal University, 500 Dongchuan Road,Shanghai 200241, People’s Republic of China

(WF) School of Mathematical Sciences and Shanghai Key Laboratory of Pure Mathematics and Mathematical Practice,East China Normal University, 500 Dongchuan Road,Shanghai 200241, People’s Republic of China