Journal of Nonlinear Mathematical Physics

Volume 26, Issue 1, December 2018, Pages 107 - 132

Rational solutions to Q3δ in the Adler-Bobenko-Suris list and degenerations

Song-lin Zhao*
Department of Applied Mathematics, Zhejiang University of Technology Hangzhou, 310023, Zhejiang, P.R. China,
Da-jun Zhang*
Department of Mathematics, Shanghai University Shanghai, 200444, P.R. China,
*Corresponding author
Corresponding Authors
Song-lin Zhao, Da-jun Zhang
Received 15 October 2017, Accepted 25 August 2018, Available Online 6 January 2021.
10.1080/14029251.2019.1544793How to use a DOI?
NQC equation; ABS list; Casoratian; rational solutions

We derive rational solutions in Casoratian form for the Nijhoff-Quispel-Capel (NQC) equation by using the lattice potential Korteweg-de Vries (lpKdV) equation and two Miura transformations between the lpKdV and the lattice potential modified KdV (lpmKdV) and the NQC equation. This allows us to present rational solutions for the whole Adler-Bobenko-Suris (ABS) list except Q4. The known Miura transformation for soliton solutions between the NQC equation and Q3δ and the known degenerations for solitons from Q3δ to Q2, Q1δ, H3δ, H2 and H1 in the ABS list are used. We show that the Miura transformation and degenerations are valid as well for rational solutions which are usually considered as “long-wave-limit” of solitons. All the rational solutions can be expressed in terms of {zj} which are linear functions of (n, m).

© 2019 The Authors. Published by Atlantis and Taylor & Francis
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (

1. Introduction

In recent years multidimensional consistency [22] has become increasingly popular as one of interpretations of integrability for lattice equations. With this property and two mild additional requirements on lattice equations: symmetry and the so-called ‘tetrahedron property’, Adler, Bobenko and Suris (ABS) classified integrable affine linear models defined on an elementary quadrilateral [3]. Their results are known as the ABS list, which consists of nine lattice equations: Q4, Q3δ, Q2, Q1δ, A2, A1δ, H3δ, H2, H1. Some of these equations have been known before, for example, H1 is the lattice potential Korteweg-de Vries (lpKdV) equation [21], H3δ=0 is the lattice potential modified KdV (lpmKdV) equation [21], Q1δ=0 is the lattice Schwarzian KdV (lSKdV) equation [20] and Q4 is known as the Adler’s equation [2] which is the nonlinear superposition formula of the Krichever-Novikov equation. After introducing some new parameters [19], the lattice equations given originally by ABS [3] can be written as

where in (1.1a) (p, P) = 𝔭 and (q, Q) = 𝔮 are the points on the elliptic curve
and in (1.1d)

Q4 will be considered elsewhere, which is not listed here. We omit A1δ and A2 from the above list because of the equivalence between A1δ and Q1δ by u → (−1)n+mu, and between A2 and Q3δ=0 by uu(−1)n+m. In Eqs. (1.1), δ is a constant; u = un,m := u(n, m) denotes dependent variable of lattice points labeled by (n, m) ∈ 𝕑 2; p and q are continuous lattice parameters associated with the grid size in the directions of the lattice given by the independent variables n and m, respectively; notations with elementary lattice shifts are denoted by


There are many ways of degenerations among the lattice equations in the list (1.1) [3, 4, 19].

Various approaches have been shown to be significant in deriving soliton solutions for the ABS list as evidenced by a series of papers. Atkinson, Hietarinta and Nijhoff constructed N-soliton solutions to Q3δ in terms of the τ-function of the Hirota-Miwa equation. The corresponding solutions were expressed by the usual Hirota’s polynomial of exponentials [5]. By developing Hirota’s direct method, Hietarinta and Zhang derived N-soliton solutions to H-series of equations and Q1δ [14]. Their method is algorithmic and based on multidimensional consistency, progressing in each case from background solution to 1-soliton solution and to N-soliton solutions, where many Casoratian shift formulae were established. Meanwhile, Nijhoff and collaborators proposed Cauchy matrix approach [19] to catch the N-soliton solutions for lattice equations in the ABS list except Q4. The authors of the present paper extended Cauchy matrix approach to a generalised case [30], which can be used to construct more kinds of exact solutions beyond soliton solutions for integrable systems (see also Ref. [25]), e.g., multiple-pole solutions. Inverse Scattering Transform was also established to solve some ABS equations [8, 9]. As the ‘master’ and the most complicate equation in this list, Q4 was solved by using Bäcklund transformation [6, 7].

Different from soliton solutions, rational solutions are usually expressed by fraction of polynomials of independent variables. Generally speaking, such type of solutions can be derived from soliton solutions through a special limit procedure (see Refs. [1, 26, 31] as examples). Compared with the case in continuous integrable systems, it is more difficult to get rational solutions of lattice equations. In spite of this, until now much progress has been got. Algebraic solutions and lump-like solutions for the Hirota-Miwa equation were, respectively, given in Refs. [17] and [12]. With the help of bilinear method [14], rational solutions for H3δ and Q1δ as well as the lattice Boussinesq equation were obtained in recent papers [23, 24]. Besides, by imposing reduction conditions on rational solutions of the Hirota-Miwa equation, rational solutions for the lpKdV equation and two semi-discrete lpKdV equations were obtained [10].

Recently, in [28] a transformation approach was employed to construct rational solutions for the ABS list (1.1) except Q3δ. Those transformations used in [28] are nonauto-Bäcklund transformations (Miura-type transformation) in which spectral parameters are absent. This is just the case of rational solutions that requires spectral parameters varnished. More examples can be found in [28]. However, the transformation connected to Q3δ is too complicated to be used for generating rational solutions.

If we forget about Q4, then Q3δ can act as a top equation in the ABS list in the sense that other “lower” equations can be obtained as degenerations of Q3δ [19]. Note that solutions in terms of Cauchy matrix given in [19, 30] are not available to generate rational solutions by taking “long-wave-limit” as done in Casoratian form (cf. [23]).

In this paper we aim to construct rational solutions to Q3δ. It is known that the Nijhoff-Quispel-Capel (NQC) equation (cf. [21]) is in some sense Q30 and by transformation solutions of Q3δ can be expressed in terms of solutions of the NQC equation [19]. Our strategy is the following. First we construct Casoratian solutions for the NQC equation so that rational solutions can be included. Then we examine the degeneration procedure given in [19] and present rational solutions for Q3δ and “lower equations” in the ABS list.

The paper is organised as follows. In Sec. 2, we explain our plan of solving the NQC equation and describe relations between the NQC equation and Q3δ. In Sec. 3, we solve the NQC equation by means of bilinear method and derive its rational solutions in Casoratian form. In Sec. 4, degenerations of Q3δ are analyzed and rational solutions for “lower equations” Q2, Q1δ, H3δ, H2 and H1 in (1.1) are obtained. Sec. 5 is for conclusions. In addition, two appendices are given as complements to the paper.

2. NQC and Q3δ

2.1. Plan of solving the NQC equation

The celebrated NQC equation [21] takes the form

where S(a, b) = S(b, a) are functions of (n, m) with (a, b) being branch point parameters, p and q are spacing parameters of n and m, respectively.

The known solutions of the NQC equation are obtained by means of Cauchy matrix [19, 30], which do not allow to take “long-wave-limit”. To derive Casoratian solutions of this equation, we introduce the following system

together with assuming symmetric property
where V(a) is a function of (n, m) with a as a parameter, V(b) = V(a)|ab, w is a function of (n, m) but independent of (a, b). Before we proceed, let us give some remarks on the above system. Firstly, this system is consistent and solvable. In fact, all the equations in (2.2) appeared in a same Cauchy matrix scheme (e.g. Eqs. (9.47, 9.53, 9.54) in [13] or (2.21, 2.37, 2.41) in [19]), which means they have solutions in terms of Cauchy matrix. Secondly, (2.2) contains Miura (nonauto-Bäcklund) transformations of several integrable equations. In fact, (2.2g) is known as the lpKdV equation [20]; the right hand sides of (2.2c,2.2d) yield a parameter-extended lpmKdV equation
which comes from the right hand sides of (2.2e,2.2f) as well as from eliminating w from (2.2c)(2.2f); with assumption (2.3), eliminating V from (2.2a, 2.2b) one gets the NQC equation (2.1). Thus, (2.2a)(2.2f) are understood as Miura transformations connecting the lpKdV (2.2g), extended lpmKdV (2.4) and NQC equation (2.1). Finally, we can conclude that under assumption (2.3) once the system (2.2) is solved, so is the NQC (2.1).

2.2. From NQC to Q3δ

In Ref. [19], soliton solutions for Q3δ (1.1a) were written as a linear combination of four terms each of which contains as an essential ingredient the soliton solution of the NQC equation (2.1). Since solutions of the NQC equation can be provided by system (2.2) with assumption (2.3), in the following we describe relations between Q3δ and the system (2.2) together with (2.3).

Theorem 2.1.

The solution of Q3δ (1.1a) is formulated by

in which S(a, b) satisfies the system (2.2) and symmetry (2.3), function Ϝ(a, b) is defined as
and P, Q are defined by (1.2); A, B, C and D are constants subject to the constraint

The proof is similar to the one given in [19]. We skip it here and leave it in Appendix A.

According to this Theorem, any solution S(a, b) solved from the system (2.2) with (2.3), including rational solutions, will generates a solution to Q3δ via formula (2.5).

3. Rational solutions to the NQC equation (2.1)

In this section, we construct rational solutions for the NQC equation (2.1) by solving system (2.2). Bilinear method will be employed and solutions will be presented in terms of Casoratians. Some Casoratian techniques developed in the literatures [14, 16] will be adopted.

3.1. Preliminary

Casoratian can be viewed as a discrete version of Wronskian. Let us consider functions ϕj with 5 independent variables n, m, α, β, l ∈ 𝕑:

where ρj±, kj, p, q, a, b, c ∈ ℂ and we denote it by ϕj(l) without confusion. Let

Define a N × N Casoratian

where M^ is a compact form standing for consequtive columns (0, 1, 2,...,M) (cf. [11]). With this notation it is easy to understand |N2^,N|=|0,1,,N2,N|, etc. In addition to the above notations, we need the following Laplace expansion identity for Casoratian verification [11].

Lemma 3.1.

Suppose that G is a N × (N − 2) matrix, and a,b,c,d are Nth-order column vectors, then


Note that ϕj(l) defined in (3.1) can be regarded as a 5-dimensional function isotropically defined on 5 directions (n, m, α, β, l) together their spacing parameters (p, q, a, b, c). For convenience, we label these 5 directions by

together with their spacing parameters by

We also introduce shift operators T±ni,


It is then easy to find shift relations for ϕ(l)a


Besides ϕj(l) in (3.1), we also introduce auxiliary functions and vecors


These vectors are necessary in Casoratian verifications (cf. [14]). They obey shift relations slightly different from each other. Suppose σ(l) is one of the above vectors and ε = (ε1, ε2, ε3, ε4). Then shift relations of σ(l), including (3.6), can be expressed through a universal formula

where ε varies with σ(l):

Note also that these vectors are related to ϕ(l) by


In addition, it is remarkable that ϕ(l) composed by (3.1) is not the only vector that obeys the relation (3.6). One can easily find that taking ϕ1(l) from (3.1) and defining new φj(l)=k1j1φ1(l) for j = 2, 3,...,N, ϕ(l) constructed in such a way satisfies (3.6) as well. More choices for ϕ(l) can be found in Sec. 3.3.

3.2. Bilinearization of (2.2) and Casoratian solutions

Under dependent variable transformations

system (2.2) is transformed into bilinear forms
where (3.13f) and (3.13g) compose the bilinear lpKdV equation [14] and in the derivation of (3.13c)(3.13e), we have made use of (3.13f) and (3.13g).

Casoratian solutions to bilinear system (3.13) can be summarized in the following Theorem.

Theorem 3.1.

The Casoratians

composed by ϕ(l), solve the bilinear system (3.13), where x and y are arbitrary nonzero constants and we require the basic column vector ϕ(l) satisfies shift relation (3.6) and assume there are invertible matrices A[ni] to define auxiliary vectors σ[ni](l) via (3.9) that obey shift relations (3.8). To meet the symmetric relation (3.12) and the assumption that V (a) is related to parameter a but independent of b, V(b) = V(a)|ab and w is independent of (a, b), one needs to impose (3.14)|a=b=0.


Due to the shift relation (3.6), the Casoratian |N1^| can be defined in terms of shifts of any variable ni and they are same (refer to Eqs.(2.22) and (2.24) in [14]). For the Eqs. (3.13a)(3.13e), with the notations in (3.14) and (3.5), they can be expressed as the following:

where i, j, k ∈ {1, 2, 3, 4}, ijk. For a certain triplet {ni, nj, nk}, this is the well known Hirota-Miwa equation [18]. It is also one of bilinear equations of H3 (see (5.20a,b) in [14]). Thus the proof for (3.15) under the condition of Theorem 3.1 has been given in [24].

The bilinear lpKdV equations (3.13f) and (3.13g) can be proved following the procedure given in [27] where the case c = 0 was handled. Here we only prove (3.13f), and the other can be treated similarly. The down-tilde-hat version of 31 is


For (pc)f˜+g˜, f̂, (qc)f̂+ĝ, f˜ and (pq)f̂˜, we use (B.1b) and (B.1a) with i = 2, respectively, i = 1 and (B.1c) with i = 1, j = 2. Then we have

where we have utilized Lemma 3.1, in which G=N3^, (a, b, c, d) = (ϕ(N − 1), Tnϕ(N − 2), ϕ(N − 2), Tmϕ(N − 2)).

3.3. Solutions for system (3.6) and (3.8)

Note that system (3.6) and (3.8) with general invertible matrices A[ni] are difference equations for unknown functions ϕ(l) and σ[ni](l). Solutions for these systems can be classified according to the canonical forms of A[ni]. Let us list them out case by case.

3.3.1. Soliton solutions

When A[ni] are diagonal matrices defined as (3.10), we take ϕj(l) to be (3.1) and φj(l), ψj(l), ϖj(l), χj(l) to be (3.7). Parameters x and y in (3.14) are taken as


This is the case to generate solitons.

There is a singularity for S(a, b) defined in (3.11) when b = −a. In the following, we define S(a, −a). For ϕj(l) given in (3.1) and x, y defined above, it is easy to check


Thus, by means of the L’Hopital rule, we can define


3.3.2. Jordan block solutions

To present elements of the basic Casoratian column vector of this case, we first introduce lower triangular Toeplitz (LTT) matrices which are defined as


Note that all the LTT matrices of same order compose a commutative set in terms of matrix product. Canonical form of such a matrix is a Jordan matrix. LTT matrices play an important role in generating multiple-pole (or limit) solutions (cf. [25, 26, 30, 31]).

When A[ni] takes a LTT form

the vector ϕ(l) that is used to construct Casoratians can be taken as
where ρ± ∈ ℂ, 𝒜± are two arbitrary LTT matrices of N-th order. Accordingly, auxiliary vectors σ[ni](l) are defined through (3.9).

It is known that Jordan block solutions can be understood as limit solutions of solitons (cf. [26]).b Employing a same limit procedure on the soliton case of Sec. 3.3.1, we can find in this case that

and S(a, −a) is given by (3.18) but now with column vector (3.20).

3.3.3. Rational solutions

Rational solutions are formally generated from Jordan block solutions by taking k = 0 in (3.20a) and (3.19). However, to avoid trivial solutions, in practice we do the following. Consider a generating function

and expand ϕ(l, k) in terms of k at k = 0, where in practice we usually replace (c ± k)l with (c ± k)l+l0 and suppose that l0 is either an integer large enough or a non-integer so that the derivative ks(c+k)l+l0|k=00.

In special cases when ρ+ = −ρ the generating function ϕ(l, k) is an odd function of k and the remained in its expansion are only odd order terms of k. We denote coefficients of these terms by


Then Casoratian column ϕ(l) for rational solutions can be taken as


A[ni] = (γs,j)N×N is taken as

and those auxiliary vectors σ[ni](l) are defined through (3.9) accordingly. When ρ+ = ρ one can also take
and in this case matrix A[ni] is still defined through (3.25) and σ[ni](l) are defined through (3.9). Since both cases share a same A[ni], a more general choice for ϕ(l) can be
where ϕ±(l) are defined in (3.24) and (3.26), and 𝒜± are two arbitrary LTT matrices of N-th order.

In this case, parameters x and y are

and S(a,a) is still given in the form (3.18) but with column vector (3.27).

3.3.4. Rational solutions: revisit

It will be interesting to have a close look at the explicit formulae of ηj defined through (3.23). For more convenience we consider vector

where ϕj(l) is defined in (3.1) but here we take c = 1 in this mini section. Denote Casoratians |N1^| and |N2^,N| composed by the above ϕ′(l) by f′ and g′ + Nf′, respectively. Then, compared with those Casoratians f and g that are composed by (3.1) with c = 1, we have

Due to gauge property of discrete Hirota bilinear equations (cf. [14, 15]), both (f, g) and (f′, g′) solve bilinear equations (3.13). Consequently, system (2.2) admits an alternative expressions for soliton solutions

where f′ and g′ are defined through (3.31),
with finally taking α = β = 0. For S(a, −a),

With regard to rational solutions, we start from (3.30) with a general k, i.e.

where we specially take
with arbitrary constants γj, and we expand ϕ′±(l, k) as

Similar to the treatment in [28], we rewrite ϕ′±(l, k) as


Then by comparison of (3.37) and (3.38), we can find all {ηj±} can be expressed in terms of {xj}. The formula is given by (cf. [28])


Explicit forms of some ηj+ are


With these results in hand, let us summarize rational solutions of the system (2.2) by the following Theorem.

Theorem 3.2.


where {ηj+} are defined in (3.40), and 𝒜± are two arbitrary LTT matrices of N-th order. Then, rational solutions of the system (2.2) are expressed by (3.32), where
are Casoratians composed by (3.41c),
and finally we need to take α = β = 0 in (3.42). S(a, −a) is defined by

We remark that in special case f′(ϕ′odd(l)) leads to a rational solution one order higher than f′(ϕ′even(l) does (see Lemma 5.4 in Ref. [28]). Here we write out explicit forms of some f′ and g′ composed by ϕ′odd(l) without any restriction on α, β:


A second remark is presented through the following Proposition.

Proposition 1.

For the Casoratian f′(ϕ′odd(l)), the following relations hold:



For ϕ′+(l, k) defined in (3.35), it is easy to see


It then follows from (3.37) and (3.38) that

and consequently we have
which leads to
and moreover

Making use of these relations, (3.45) can be verified directly.

We also remark that some properties of f′(ϕ′odd(l)) can be found in Theorem C.2 in Appendix C of Ref. [28]. One of the properties is:

Proposition 2.

f′(ϕ′odd(l)) is a polynomial of {x1, x3,...,x2N−1} and coefficients are independent of {p, q, a, b}. (So are g′(ϕ′odd(l)) and ζ′(ϕ′odd(l)) due to Proposition 1.)

Thanks to such a property, S(a, −a) is well defined by the formula


In fact, write f′ = f′[x] = f′[x1, x3,...,x2N−1]. Then we have


By Taylor expanding TaTb f′[x] at (ε1,...,ε2N−1) = (0,...,0), we have


Noticing that each ε2j−1 has a factor a + b, i.e.

after substituting (3.54) into (3.52), it is easy to see that the term a + b in the denominator of (3.52) can be eliminated and as a definition (3.52) is valid to S(a, −a).

The first two rational solutions of the NQC equation (2.1) are


Rational solutions of Q3δ (1.1a) are given by (2.5), where Ϝ(a, b) is defined in (2.6) and A, B, C, D obey the constraint (2.7). The first two rational solutions of Q3δ are


4. Degeneration of rational solutions

We now consider the problem of degeneration of rational solutions of Q3δ down to those of the “lower” equations Q2, Q1δ, H3δ, H2 and H1 in the ABS list (1.1). To do so we follow the degenerations given in Ref. [19] which are limits on the parameters a and b and the dependent variable u, where a small parameter ε is introduced, and all degenerations are obtained in the limit ε → 0. The degeneration relations between Q3δ and the “lower equations” are depicted as Fig.1 [19].

Fig. 1.

Degeneration relation

4.1. Q3δ −→ Q2

The degeneration from Q3δ to Q2 is implemented through taking

and making the following replacements of constants in (2.5),

Then, rational solutions for Q2 are given as

in which
and ξ0, A and D are arbitrary constants; ρ(a) is defined by (A.10). In (4.4), when N = 1 and N = 2, we have
where zj is defined by (3.57). However, it is hard to give an explicit formula for Z(a, −a) in terms of Casoratians.

Pure rational solutions of Q2 are obtained by taking A = D = 0. The first two of them are


4.2. Q2 −→ Q1δ

To achieve rational solutions for Q1δ, consider degeneration of (4.3) by taking


Meanwhile, we replace the constants appearing in solution (4.3) by


Then the rational solutions for Q1δ can be described as

where ρ(a) is defined by (A.10) and constants A, B, D are chosen to satisfy the constraint

Some explicit S(a, b) can be seen from (3.56). The first two solutions given by (4.8) are


4.3. Q3δ −→ H3δ

By setting

rational solutions to H3δ can be obtained from (2.5) and take the form
in which V1a) come from Taylor expressions:
where P, Q are defined by (1.3) and constants A, B, C and D are subject to the constraint

Let us have a close look at V1(a). Considering S(a, b) defined by (3.52), which is valid for b = −a as well, we find


Since f′[x] is a polynomial of {x1, x3,...,x2N−1}, (see Proposition 2), substituting b=1ε2 and after taking ε → 0, all the terms of 1bj vanish, i.e. all the shifts with respect to β do not make sense any longer. Thus, we immediately reach


When N = 1, V1(a) reads

and when N = 2, V1(a) is

The first two solutions of H3δ are


4.4. Q2 −→ H2

The degeneration from Q2 to H2 can be obtained by setting


Substituting (4.19) into (4.3), combined with

and the following choice for the constants
with unconstrained constants γ′, γ−1, rational solutions for H2 can be obtained, which are given by

In the following we derive explicit Casoratian forms for S(0) and S(1). For convenience we take f′ = f′(ϕ′odd(l)). Recalling the analysis for f′ in Sec. 3.3.4, noticing Taylor expansion (3.54) and relation (3.55), we have

where f′[z] = f′[z1, z3,...,z2N−1]. When taking a=1ε we have

Thus, compared with (4.20) and making use of Proposition 1, yield


The case of S(a, a) is little bit complicated. Let εj=α1ajβ1bj. Then TaTbf′[x] can be written as


Taking b=a=1ε and then β = α = 0 yields


Then we have

where we have taken a=1ε. After a comparison with (4.20) and making use of Proposition 1, we reach

The first two S(0) are


The first two S(1) are


And the first two rational solutions of H2 are


4.5. Q1δ −→ H1

The rational solution to H1 can be obtained from (4.8) through degeneration. Substituting

into (4.8) and using
with constants γ′, γ−1, we have
which presents rational solutions of H1, where γ′, A and B satisfy A2B2 = −1. Here S(0) is same as the one in solution (4.22). The first two rational solutions given by (4.31) are

5. Conclusions

In this paper we have derived solutions in terms of Casoratians for the NQC equation. It turns out that the τ-function f with auxiliary directions (α, β) satisfies the Hirota-Miwa equations that are defined in any three directs of (n, m, α, β), but finally (α, β) are restricted to be (0,0). This indicates that, with regard to generating soliton solutions for the ABS equations (except Q4), the Hirota-Miwa equation acts as a master equation to govern the τ-function f.

Rational solutions of the NQC equation are obtained by considering the case that auxiliary matrices A[ni] take the form (3.25) with k = 0. Such rational solutions can be equivalently derived by taking “long-wave-limit”. In this paper, we examined the significance of our rational solutions in the Miura transformation (2.5) from the NQC equation to Q3δ and in the degenerations from Q3δ to the “lower” equations in the ABS list. As a result, we obtained rational solutions in terms of Casoratians for the NQC equation, Q3δ, Q1δ, H3δ, H2 and H1. However, for Q2 we did not find an explicit form of Z(a, −a) in terms of Casoratians. This will be considered in the future. Besides, it is worthy to mention that {xj} defined through (3.38) play an important role in the analysis of rational solutions.

Compared with [28], one can find that some rational solutions of the ABS equations were obtained by means of Bäcklund transformations in [28] can not be derived in this paper. For example, for the rational solutions of Q1δ, (5.28c) in [28] can not be obtained from (4.8) of the present paper, and vice versa.

A further remark is about solving the NQC equation (2.1). In this paper we solve the equation through solving the system (2.2) that contains the lpKdV equation and two Miura transformations, and is based on the Cauchy matrix scheme that provides many clues for bilinearisation. This system is also used in deriving the Miura transformation (2.5) between the NQC equation and Q3δ (see Appendix A). In some sense, (3.13a3.13e) can be viewed as a bilinear form of the NQC equation. If we only solve the NQC equation, (2.2a, 2.2b) together with extensions in α-and β-direction are enough (cf. Sec. 9.5.3 in [13]). Note that (2.2a, 2.2b) can be considered as a generalisation (with parameters (a, b)) of the class of Bäcklund transformations investigated in [29] (see (3.1) with (3.5) in [29]). This motivates us to reconsider [29] with parameters, which will be done in the future.

Appendix A. Proof of Theorem 2.1


The proof is similar to the one given in Ref. [19]. We need three steps. Let us proceed them step by step.

Step # 1. We first introduce a new associated dependent variable U that is given by:


Lemma A.1.

For u defined in (2.5) and the associated variable U defined in (A.1) the following hold:

where w satisfies system (2.2). The relations (A.2) hold for arbitrary coefficients A, B, C, D.


We denote u and U, respectively, defined by (2.5) and (A.1) as

where uA(a, b) = AϜ(a, b)[1 − (a + b)S(a, b)] and UA(a, b) = (a + b)AϜ(a, b)V(a)V(b). By direct computation, we get
where relations (2.2a), (2.2b) and (2.2d) have been used. In terms of symmetric property, we also havex

By adding the four equations in (A.4) and (A.5), we arrive at (A.2a). Similar analysis can be done to (A.2b)(A.2d).

Step # 2.

Lemma A.2.

The following identities hold

in which the 2 × 2 matrix 𝒜 is given by
provided that Sa, ±b), va) and vb) satisfy the system (2.2) together with symmetry S(a, b) = S(b, a) with parameters ±a, ±b.


We consider the following 2 × 2 matrices:


Evaluating the entries in these matrices we obtain

in which again the dots in each line on the right hand sides stand for similar terms with (a, b) replaced by (a, −b), (−a, b), (−a, −b) and A replaced by B, C, D respectively. On the right-hand side of this expression we have introduced vectors
in which the plane-wave factor ρ(a) is given byx

A similar computation as above yields

whilst U and Ũ can be written as

Thus, we find that the matrix L in (A.8) can be written as

and similarly for the matrix M. Using now the general determinantal identity
for any collection of r pairs of r-component column vectors xi, yi (the superindex T denoting transposition), we obtain the following result:

It remains to compute the determinant of matrix (r(a),r˜(a)) whose columns are the 2-component vectors r(a) and r˜(a). This is done as follows

where we have used the fact that S(a, b) = S(b, a). Thus, putting everything together we obtain the result:

On the other hand, a direct computation of the determinant gives:


Comparing the two expressions for det(L) from (A.13) and (A.14) we obtain the first equation in Lemma A.2.

Step # 3. The last step is by combining the relations (A.2) and (A.6) as well as the lpKdV equation (2.2g) to assert that u solves the Q3δ equation. In fact, multiplying for instance (A.2b) by (A.2d) and using (A.6b), where we identify

according to (2.7), we obtain from the lpKdV equation:
where we have used the fact that 𝔭, 𝔮 are on the elliptic curve (1.2) and the identification of U˜U˜^ via the result (A.6) of Lemma A.2. From the last step, after some cancellation of terms, we obtain Q3δ, in the form (1.1a), for the function u, which completes the proof of the theorem.

Appendix B. Casoratian shift formulae

We list some shift formulae for the Casoratians (3.14), where the basic column vector ϕ(l) satisfies shift relation (3.6) and auxiliary relations (3.9). For convenience, we introduce notations T¯niφ(l)=A[ni]Tni(A[ni]1φ(l)) with i = 1, 2, 3, 4. These Casoratian shift formulae are

where i = 1, 2, 3, 4.


The authors are grateful to the referee for the invaluable comments. This project is supported by the NSF of China (Nos. 11401529, 11301483, 11371241, 11631007) and the NSF of Zhejiang Province (Nos. LY17A010024, LY18A010033).



We keep the variable l in the notation ϕ(l) since Casoratians are defined in terms of shifts in l.


One can consider solitons where we take ϕ1 to be (3.1) and take ϕj (where j ≥ 2) by just replacing k1 with kj in ϕ1. Then we expand kj = k1 + εj at εj = 0 and take limit εj → 0 in the rational forms (3.11) successively for j = 2, 3,...,N. As a result, we get Jordan block solutions. For more details one can refer to [26].


[10]W. Feng, S.L. Zhao, and Y. Shi, Rational solutions for lattice potential KdV equation and two semi-discrete lattice potential KdV equations, Z. Naturforsch, Vol. 71, No. 2)a, 2016, pp. 121-128.
[24]Y. Shi and D.J. Zhang, Rational solutions of the H3 and Q1 models in the ABS lattice list, SIGMA, Vol. 7, 2011, pp. 11. 046
[26]D.J. Zhang, Notes on solutions in Wronskian form to soliton equations: KdV-type, arXiv:nlin.SI/0603008, 2006, pp. 45.
[27]D.J. Zhang and J. Hietarinta, Generalized solutions for the H1 model in ABS list of lattice equations, Nonl. Mod. Math. Phys: Proceedings of the First International Workshop, AIP Conference Proceedings, Vol. 1212, 2010, pp. 154-161.
[28]D.D. Zhang and D.J. Zhang, Rational solutions to the ABS list: Transformation approach, SIGMA, Vol. 13, 2017, pp. 24. 078
Journal of Nonlinear Mathematical Physics
26 - 1
107 - 132
Publication Date
ISSN (Online)
ISSN (Print)
10.1080/14029251.2019.1544793How to use a DOI?
© 2019 The Authors. Published by Atlantis and Taylor & Francis
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (

Cite this article

AU  - Song-lin Zhao
AU  - Da-jun Zhang
PY  - 2021
DA  - 2021/01/06
TI  - Rational solutions to Q3δ in the Adler-Bobenko-Suris list and degenerations
JO  - Journal of Nonlinear Mathematical Physics
SP  - 107
EP  - 132
VL  - 26
IS  - 1
SN  - 1776-0852
UR  -
DO  - 10.1080/14029251.2019.1544793
ID  - Zhao2021
ER  -