This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Consistency around a cube property of Hirota’s discrete KdV equation and the lattice sine-Gordon equation

Nobutaka Nakazono Institute of Engineering, Tokyo University of Agriculture and Technology, 2-24-16 Nakacho Koganei, Tokyo 184-8588, Japan. nakazono@go.tuat.ac.jp
Abstract.

It has been unknown whether Hirota’s discrete Korteweg-de Vries equation and the lattice sine-Gordon equation have the consistency around a cube (CAC) property. In this paper, we show that they have the CAC property. Moreover, we also show that they can be extended to systems on the 3-dimensional integer lattice.

Key words and phrases:
discrete integrable systems; partial difference equation; Lax pair; discrete KdV equation; lattice sine-Gordon equation
2020 Mathematics Subject Classification:
33E30, 35Q53, 37K10, 39A14, 39A36, 39A45

1. Introduction

In this paper, we focus on the following two integrable partial difference equations (PΔ\DeltaEs). One is Hirota’s discrete Korteweg-de Vries (dKdV) equation [7, 23, 12]:

ul+1,m+1ul,m=βm+1αlul,m+1βmαl+1ul+1,m,u_{l+1,m+1}-u_{l,m}=\dfrac{\beta_{m+1}-\alpha_{l}}{u_{l,m+1}}-\dfrac{\beta_{m}-\alpha_{l+1}}{u_{l+1,m}}, (1.1)

where (l,m)2(l,m)\in\mathbb{Z}^{2} is a lattice parameter, ul,mu_{l,m}\in\mathbb{C} is its function and αl,βm\alpha_{l},\beta_{m}\in\mathbb{C} are parameters that depend only on ll and mm, respectively, and the other is the lattice sine-Gordon (lsG) equation[13]:

xl+1,m+1xl,m=(pl+1qmxl+1,mqmpl+1xl+1,m)(qm+1plxl,m+1plqm+1xl,m+1),\dfrac{x_{l+1,m+1}}{x_{l,m}}=\left(\dfrac{p_{l+1}-q_{m}x_{l+1,m}}{q_{m}-p_{l+1}x_{l+1,m}}\right)\left(\dfrac{q_{m+1}-p_{l}x_{l,m+1}}{p_{l}-q_{m+1}x_{l,m+1}}\right), (1.2)

where (l,m)2(l,m)\in\mathbb{Z}^{2} is a lattice parameter, xl,mx_{l,m}\in\mathbb{C} is its function and pl,qmp_{l},q_{m}\in\mathbb{C} are parameters that depend only on ll and mm, respectively.

The PΔ\DeltaE (1.1) is an integrable non-autonomous discrete version of the Korteweg-de Vries (KdV) equation[14]:

wt+6wwx+wxxx=0,w_{t}+6ww_{x}+w_{xxx}=0, (1.3)

where (t,x)2(t,x)\in\mathbb{C}^{2} and w=w(t,x)w=w(t,x)\in\mathbb{C}, which is known as a mathematical model of waves on shallow water surfaces. On the other hand, the PΔ\DeltaE (1.2) is an integrable non-autonomous discrete version of the sine-Gordon equation:

ϕttϕxx+sinϕ=0,\phi_{tt}-\phi_{xx}+\sin{\phi}=0, (1.4)

where (t,x)2(t,x)\in\mathbb{C}^{2} and ϕ=ϕ(t,x)\phi=\phi(t,x)\in\mathbb{C}, which is known as a motion equation of a row of pendulums hanging from a rod and being coupled by torsion springs.

The consistency around a cube (CAC) property[19, 18, 22, 21, 20] is known as an integrability of quad-equations. (For the CAC property and the definition of quad-equation, see Appendix A.1.) As classifications of quad-equations by using the CAC property, the lists of equations by Adler-Bobenko-Suris (ABS)[1, 2], Boll[4, 5] and Hietarinta [8] are known. These classifications have been done under the following additional conditions:

  1. [1]:

    All face-equations on the same cube are similar equations differing only by the parameter values and are invariant under the group D4D_{4} of the square symmetries. Moreover, the CAC cube also has the tetrahedron property.

  2. [2, 4, 5]:

    The face-equations on the same cube are allowed to be different, and no assumption is made about symmetry. However, the tetrahedron property is still required.

  3. [8]:

    No assumptions are made for symmetry or the tetrahedron property, but all face-equations are assumed to be homogeneous quadratic. The face-equations on the same cube allow different equations on the three orthogonal planes but are the same equations on parallel planes.

Note here that a CAC-cube means a cube with the six quad-equations on its faces (face-equations) which has the CAC property. However, the PΔ\DeltaEs (1.1) and (1.2) are not included in these lists, and the structures of their consistency have not been reported.

Recently, it has been shown that the PΔ\DeltaEs (1.1) and (1.2) have the consistency around a broken cube (CABC) property[11, 15]. (For the CABC property, see Appendix A.2.) By considering the Bäcklund transformation (BT)

ul,mvl+1,m+1βmαl+1ul+1,m+γαlvl,m+1=0,\displaystyle u_{l,m}-v_{l+1,m+1}-\dfrac{\beta_{m}-\alpha_{l+1}}{u_{l+1,m}}+\dfrac{\gamma-\alpha_{l}}{v_{l,m+1}}=0, (1.5a)
(ul,mvl,m)(βmαlul,m+vl,m+1)βm+γ=0,\displaystyle(u_{l,m}-v_{l,m})\left(\dfrac{\beta_{m}-\alpha_{l}}{u_{l,m}}+v_{l,m+1}\right)-\beta_{m}+\gamma=0, (1.5b)
βmαlul,mγαlvl,mul+1,m+vl+1,m=0,\displaystyle\dfrac{\beta_{m}-\alpha_{l}}{u_{l,m}}-\dfrac{\gamma-\alpha_{l}}{v_{l,m}}-u_{l+1,m}+v_{l+1,m}=0, (1.5c)

from the PΔ\DeltaE (1.1) to the multi-quadratic equation

βm(vl,mvl+1,mvl,m+1vl+1,m+1)2+(vl,mvl+1,m+1)(vl+1,mvl,m+1)(γvl,mvl+1,m)(γvl,m+1vl+1,m+1)+(αl+1vl,mαlvl+1,m+1)(αlvl+1,mαl+1vl,m+1)(αl+αl+1)(γvl,mvl+1,m+γvl,m+1vl+1,m+12vl,mvl+1,mvl,m+1vl+1,m+1)+(αl+1vl,mvl,m+1+αlvl+1,mvl+1,m+1)(2γvl,mvl+1,mvl,m+1vl+1,m+1)=0,\begin{split}&\beta_{m}\left(v_{l,m}\,v_{l+1,m}-v_{l,m+1}\,v_{l+1,m+1}\right)^{2}\\ &\quad+\bigl{(}v_{l,m}-v_{l+1,m+1}\bigr{)}\bigl{(}v_{l+1,m}-v_{l,m+1}\bigr{)}\bigl{(}\gamma-v_{l,m}\,v_{l+1,m}\bigr{)}\bigl{(}\gamma-v_{l,m+1}\,v_{l+1,m+1}\bigr{)}\\ &\quad+\bigl{(}\alpha_{l+1}v_{l,m}-\alpha_{l}v_{l+1,m+1}\bigr{)}\bigl{(}\alpha_{l}v_{l+1,m}-\alpha_{l+1}v_{l,m+1}\bigr{)}\\ &\quad-(\alpha_{l}+\alpha_{l+1})\bigl{(}\gamma\,v_{l,m}\,v_{l+1,m}+\gamma\,v_{l,m+1}\,v_{l+1,m+1}-2\,v_{l,m}\,v_{l+1,m}\,v_{l,m+1}\,v_{l+1,m+1}\bigr{)}\\ &\quad+\bigl{(}\alpha_{l+1}\,v_{l,m}\,v_{l,m+1}+\alpha_{l}\,v_{l+1,m}\,v_{l+1,m+1}\bigr{)}\bigl{(}2\gamma-v_{l,m}v_{l+1,m}-v_{l,m+1}\,v_{l+1,m+1}\bigr{)}=0,\end{split} (1.6)

where γ\gamma\in\mathbb{C} is a constant parameter, the CABC property of the PΔ\DeltaE (1.1) was shown in [11]. Similarly, by considering the following BT from the PΔ\DeltaE (1.2) to Equation (1.6):

vl,m+1vl+1,m+1plpl+1xl,m(pl+1qmxl+1,m)pl+1xl+1,mqm=0,\displaystyle\dfrac{v_{l,m+1}v_{l+1,m+1}}{p_{l}p_{l+1}}-\dfrac{x_{l,m}(p_{l+1}-q_{m}x_{l+1,m})}{p_{l+1}x_{l+1,m}-q_{m}}=0, (1.7a)
vl,mvl,m+1pl2xl,m(plqmxl,m)plxl,mqm=0,\displaystyle\dfrac{v_{l,m}v_{l,m+1}}{{p_{l}}^{2}}-\dfrac{x_{l,m}(p_{l}-q_{m}x_{l,m})}{p_{l}x_{l,m}-q_{m}}=0, (1.7b)
vl,mvl+1,mplpl+1xl+1,m(plqmxl,m)plxl,mqm=0,\displaystyle\dfrac{v_{l,m}v_{l+1,m}}{p_{l}p_{l+1}}-\dfrac{x_{l+1,m}(p_{l}-q_{m}x_{l,m})}{p_{l}x_{l,m}-q_{m}}=0, (1.7c)

the CABC property of the PΔ\DeltaE (1.2) was shown in [15]. In this case, the relation between the parameters plp_{l}, qmq_{m}, αl\alpha_{l}, βm\beta_{m} and γ\gamma is given by

pl=αlγ,qm=βmγ.p_{l}=\sqrt{\alpha_{l}-\gamma},\quad q_{m}=\sqrt{\beta_{m}-\gamma}. (1.8)

Since having a CABC property does not imply not having a CAC property, the CABC property of the PΔ\DeltaEs (1.1) and (1.2) was shown, but whether they have the CAC property or not is an open problem. It is natural to consider that the CAC property can be shown using the CABC property because both are meant to show integrability and are the properties around a cube. This is the motivation for this study.

Remark 1.1.

The first condition in Equation (1.8) means that given the complex parameters αl\alpha_{l} and γ\gamma, we determine one value of the complex parameter plp_{l} that satisfies

pl2=αlγ.{p_{l}}^{2}=\alpha_{l}-\gamma. (1.9)

The same is true for the second condition in Equation (1.8). In addition, if the root symbol \sqrt{~{}\,~{}} is used in the following, it means the same thing.

In this paper, we introduce another CABC-type BT from the PΔ\DeltaE (1.1) to Equation (1.6), which is different from the CABC-type BT (1.5). Then, using a total of three CABC-type BTs, we show the following results:

  1. (i)

    derive two CAC-type BTs from the PΔ\DeltaE (1.1) to the PΔ\DeltaE (1.2);

  2. (ii)

    derive a CAC-type auto-Bäcklund transformation (auto-BT) of the PΔ\DeltaE (1.2);

  3. (iii)

    construct Lax pairs of the PΔ\DeltaEs (1.1) and (1.2);

  4. (iv)

    derive integrable systems of PΔ\DeltaEs on the 3-dimensional integer lattice containing the PΔ\DeltaEs (1.1) and (1.2).

The main contribution of this paper is not only to show that the systems (1.1) and (1.2) have the CAC property but also to find examples of quad-equations having the CAC property that do not satisfy the conditions for classifying the list of equations by ABS, Boll and Hietarinta above. Finding such examples leads to the study of making new lists of PΔ\DeltaEs that have the CAC property, i.e., it contributes to the development of the study of integrable systems.

1.1. Notation and Terminology

  • For simplicity, we use the following shorthand notations:

    ¯:ll+1,~:mm+1,^:nn+1.\overline{\rule{0.0pt}{5.0pt}\hskip 6.00006pt}:\,l\to l+1,\quad\widetilde{}\,:\,m\to m+1,\quad\widehat{}\,:\,n\to n+1. (1.10)

    For example,

    u¯=ul+1,m,u~=ul,m+1,u¯~=ul+1,m+1,u¯¯=ul+2,m,u~~=ul,m+2\overline{u}=u_{l+1,m},\quad\widetilde{u}=u_{l,m+1},\quad\widetilde{\overline{u}}=u_{l+1,m+1},\quad\overline{\overline{u}}=u_{l+2,m},\quad\widetilde{\widetilde{u}}=u_{l,m+2} (1.11)

    for u=ul,mu=u_{l,m} and

    u¯=ul+1,m,n,u~=ul,m+1,n,u^=ul,m,n+1,u¯~^=ul+1,m+1,n+1\overline{u}=u_{l+1,m,n},\quad\widetilde{u}=u_{l,m+1,n},\quad\widehat{u}=u_{l,m,n+1},\quad\widehat{\widetilde{\overline{u}}}=u_{l+1,m+1,n+1} (1.12)

    for u=ul,m,nu=u_{l,m,n}, and the like.

  • If there is a substitution notation (e.g., ll+1l\to l+1) in the subscript of an equation number, it means the equation with the substitution applied to all symbols in the corresponding equation. For example, (1.1)l→l+1 implies the following equation:

    ul+2,m+1ul+1,m=βm+1αl+1ul+1,m+1βmαl+2ul+2,m.u_{l+2,m+1}-u_{l+1,m}=\dfrac{\beta_{m+1}-\alpha_{l+1}}{u_{l+1,m+1}}-\dfrac{\beta_{m}-\alpha_{l+2}}{u_{l+2,m}}.
  • The following terminologies are used in this paper.

    1. (1)(\mathrm{1})

      A multivariate polynomial with complex coefficients that is affine linear with respect to each variable is called a multilinear polynomial. For example, the general form of a multilinear polynomial in four variables is given by

      A1xyzw+A2xyz+A3xyw+A4xzw+A5yzw+A6xy+A7xz+A8xw\displaystyle A_{1}xyzw+A_{2}xyz+A_{3}xyw+A_{4}xzw+A_{5}yzw+A_{6}xy+A_{7}xz+A_{8}xw
      +A9yz+A10yw+A11zw+A12x+A13y+A14z+A15w+A16,\displaystyle\quad+A_{9}yz+A_{10}yw+A_{11}zw+A_{12}x+A_{13}y+A_{14}z+A_{15}w+A_{16}, (1.13)

      where AiA_{i}\in\mathbb{C} are parameters.

    2. (2)(\mathrm{2})

      Let Q(x,y,z,w)Q(x,y,z,w) be a multilinear polynomial and x=f(y,z,w)x=f(y,z,w) be the solution of Q(x,y,z,w)=0Q(x,y,z,w)=0. We call Q(x,y,z,w)Q(x,y,z,w) an irreducible multilinear polynomial in four variables if the following holds:

      yf(y,z,w)0,zf(y,z,w)0,wf(y,z,w)0.\dfrac{\partial}{\partial y}f(y,z,w)\neq 0,\quad\dfrac{\partial}{\partial z}f(y,z,w)\neq 0,\quad\dfrac{\partial}{\partial w}f(y,z,w)\neq 0. (1.14)
    3. (3)(\mathrm{3})

      Consider a PΔ\DeltaE of the form

      Q(u,u¯,u~,u¯~)=0,Q(u,\overline{u},\widetilde{u},\widetilde{\overline{u}}\,)=0, (1.15)

      where u=ul,mu=u_{l,m} and QQ is an irreducible multilinear polynomial. By putting the variable ul,mu_{l,m} on the grid point (l,m)(l,m) of the integer lattice 2\mathbb{Z}^{2}, the PΔ\DeltaE (1.15) can be regarded as a relation on the quadrilateral given by

      (l,m),(l+1,m),(l,m+1),(l+1,m+1)(l,m),\quad(l+1,m),\quad(l,m+1),\quad(l+1,m+1)

      as shown in Figure 1.1. For this reason, a PΔ\DeltaE of the form (1.15) is called a quad-equation. Note that quad-equations are not necessarily autonomous but may contain parameters that evolve with ll or mm.

Refer to caption
Figure 1.1. A quadrilateral for the quad-equation (1.15).

1.2. Plan of the paper

This paper is organized as follows. In §2, using the CABC property of the dKdV equation (1.1) and the lsG equation (1.2), we show that they have the CAC property. In §3, we show how the PΔ\DeltaEs (1.1) and (1.2) can be extended to systems of PΔ\DeltaEs on the 3-dimensional integer lattice. Some concluding remarks are given in §4. In Appendix A, we recall the definitions of the CAC and CABC properties. In Appendix B, we list the Lax pairs of the PΔ\DeltaEs (1.1) and (1.2) constructed in this study. In Appendix C, we give the proof of Theorem 3.3.

2. Consistencies of the dKdV equation and the lsG equation

In this section, we give the new consistencies of the dKdV equation (1.1) and the lsG equation (1.2). See Appendix A for descriptions of the CABC and CAC properties we deal with here.

2.1. A CABC-type BT from the PΔ\DeltaE (1.1) to the multi-quadratic equation (1.6)

The following system of equations is a new CABC-type BT from the PΔ\DeltaE (1.1) to Equation (1.6):

uv~βmαl+1u¯+γαl+1v¯~=0,\displaystyle u-\widetilde{v}-\dfrac{\beta_{m}-\alpha_{l+1}}{\overline{u}}+\dfrac{\gamma-\alpha_{l+1}}{\widetilde{\overline{v}}}=0, (2.1a)
(uv~)(βmαlu+v)βm+γ=0,\displaystyle(u-\widetilde{v})\left(\dfrac{\beta_{m}-\alpha_{l}}{u}+v\right)-\beta_{m}+\gamma=0, (2.1b)
βmαluγαl+1v¯u¯+v=0,\displaystyle\dfrac{\beta_{m}-\alpha_{l}}{u}-\dfrac{\gamma-\alpha_{l+1}}{\overline{v}}-\overline{u}+v=0, (2.1c)

where u=ul,mu=u_{l,m} and v=vl,mv=v_{l,m}. By direct calculation, we can verify that the system of equations (1.1) and (2.1) has the CABC and tetrahedron properties. Indeed, the following relations can be derived from the CABC-system (1.1) and (2.1):

αl+1γv¯~=αl+1βmu¯(βmγ)uαlβmuv,v~=u+(βmγ)v¯αl+1γu¯v¯.\dfrac{\alpha_{l+1}-\gamma}{\widetilde{\overline{v}}}=\dfrac{\alpha_{l+1}-\beta_{m}}{\overline{u}}-\dfrac{(\beta_{m}-\gamma)u}{\alpha_{l}-\beta_{m}-uv},\quad\widetilde{v}=u+\dfrac{(\beta_{m}-\gamma)\overline{v}}{\alpha_{l+1}-\gamma-\overline{u}\,\overline{v}}. (2.2)

As shown in [11, 15], the Lax pair (B.1) with (B.2a) of the PΔ\DeltaE (1.1) can be constructed from the equations (2.1b) and (2.1c) by setting

vl,m=Fl,mGl,m,Φl,m=(Fl,mGl,m).v_{l,m}=\dfrac{F_{l,m}}{G_{l,m}},\quad\Phi_{l,m}=\begin{pmatrix}F_{l,m}\\ G_{l,m}\end{pmatrix}. (2.3)

2.2. A CAC-type BT from the PΔ\DeltaE (1.1) to the PΔ\DeltaE (1.2) (I)

In this subsection, using the CABC-type BTs (1.5) and (1.7), we obtain a CAC-type BT from the PΔ\DeltaE (1.1) to the PΔ\DeltaE (1.2).

Eliminating v~\widetilde{v} from the equations (1.5b) and (1.7b), we obtain

((qmplx)v+plux)((pl2qm2)vplu(plqmx))=0.\Big{(}(q_{m}-p_{l}x)v+p_{l}ux\Big{)}\Big{(}({p_{l}}^{2}-{q_{m}}^{2})v-p_{l}u(p_{l}-q_{m}x)\Big{)}=0. (2.4)

Let us assume that the following holds:

v=pluxplxqm,v=\dfrac{p_{l}ux}{p_{l}x-q_{m}}, (2.5)

which satisfies Equation (2.4). Eliminating vv from Equation (1.5b) by using Equation (2.5) and v~\widetilde{v} from Equation (1.5c) by using Equation (2.5)m→m+1, we obtain

u~u+(plqmx)(qm+1plx~)x~=0,\displaystyle\widetilde{u}u+\dfrac{(p_{l}-q_{m}x)(q_{m+1}-p_{l}\widetilde{x})}{\widetilde{x}}=0, (2.6a)
u¯u+(plqmx)(qmpl+1x¯)x=0,\displaystyle\overline{u}u+\dfrac{(p_{l}-q_{m}x)(q_{m}-p_{l+1}\overline{x})}{x}=0, (2.6b)

respectively. The equations above, together with the PΔ\DeltaEs (1.1) and (1.2), lead to the following theorem.

Theorem 2.1.

The following system has the CAC property but does not have the tetrahedron property:

u¯~u=βm+1αlu~βmαl+1u¯,\displaystyle\widetilde{\overline{u}}-u=\dfrac{\beta_{m+1}-\alpha_{l}}{\widetilde{u}}-\dfrac{\beta_{m}-\alpha_{l+1}}{\overline{u}}, (2.7a)
u~u+(plqmx)(qm+1plx~)x~=0,\displaystyle\widetilde{u}u+\dfrac{(p_{l}-q_{m}x)(q_{m+1}-p_{l}\widetilde{x})}{\widetilde{x}}=0, (2.7b)
u¯u+(plqmx)(qmpl+1x¯)x=0,\displaystyle\overline{u}u+\dfrac{(p_{l}-q_{m}x)(q_{m}-p_{l+1}\overline{x})}{x}=0, (2.7c)
x¯~x=(pl+1qmx¯qmpl+1x¯)(qm+1plx~plqm+1x~),\displaystyle\dfrac{\,\widetilde{\overline{x}}\,}{x}=\left(\dfrac{p_{l+1}-q_{m}\overline{x}}{q_{m}-p_{l+1}\overline{x}}\right)\left(\dfrac{q_{m+1}-p_{l}\widetilde{x}}{p_{l}-q_{m+1}\widetilde{x}}\right), (2.7d)

where plp_{l} and qmq_{m} are given by Equation (1.8). Since the equations (2.7a) and (2.7d) are equal to the PΔ\DeltaEs (1.1) and (1.2), respectively, we can also say that the pair of equations (2.7b) and (2.7c) is a CAC-type BT from the PΔ\DeltaE (1.1) to the PΔ\DeltaE (1.2).

Proof.

Consider the cube in Figure 2.1. In the three different ways (see Appendix A.1), x¯~\widetilde{\overline{x}} can be uniquely represented by the initial values {u,u¯,u~,x}\{u,\overline{u},\widetilde{u},x\} as

x¯~=qm+1u~((αl+1βm)(plqmx)qmuu¯x)pl+1u¯((αlβm+1)(plqmx)pluu~).\widetilde{\overline{x}}=\dfrac{q_{m+1}\widetilde{u}\,\Big{(}(\alpha_{l+1}-\beta_{m})(p_{l}-q_{m}x)-q_{m}u\overline{u}x\Big{)}}{p_{l+1}\overline{u}\,\Big{(}(\alpha_{l}-\beta_{m+1})(p_{l}-q_{m}x)-p_{l}u\widetilde{u}\Big{)}}. (2.8)

Since x¯~\widetilde{\overline{x}} depends on uu, the tetrahedron property does not hold. Therefore, we have completed the proof. ∎

Theorem 2.1 gives the following corollary.

Corollary 2.2.

Both of the PΔ\DeltaEs (1.1) and (1.2) have the CAC property.

Refer to caption
Figure 2.1. A cube for the CAC-systems (2.7) and (2.12). The uu- and xx-variables are assigned at the bottom and top vertices, respectively. The corresponding equations of a CAC-system are assigned to the faces, e.g. the equations (2.7a) and (2.7d) are assigned to the bottom and top faces, respectively. Also, as a convention, the parameters in a CAC-system are placed on the line segments. Thus, the parameters {αl,αl+1}\{\alpha_{l},\alpha_{l+1}\} are assigned to the line segments connecting the variables that represent a shift in the ll-direction (e.g., uu and u¯\overline{u}) since {αl,αl+1}\{\alpha_{l},\alpha_{l+1}\} depend on only ll. Similarly, the parameters {βm,βm+1}\{\beta_{m},\beta_{m+1}\} are assigned to the line segments connecting the variables that represent a shift in the mm-direction. On the other hand, the parameter γ\gamma is assigned to the line segments connecting the uu- and xx-varibales (e.g., uu and xx) since γ\gamma is a parameter associated with the BT from the quad-equation given by uu-variable to that given by xx-variable. For simplicity, only one parameter is shown for each in the figure, but they can also be placed on three other parallel line segments. For example, γ\gamma can also be placed on the line segment connecting u¯\overline{u} and x¯\overline{x}.
Remark 2.3.

If we choose

v=plu(plqmx)pl2qm2v=\dfrac{p_{l}u(p_{l}-q_{m}x)}{{p_{l}}^{2}-{q_{m}}^{2}} (2.9)

as the solution of Equation (2.4) instead of Equation (2.5) and construct another system of PΔ\DeltaEs, then we can verify that the resulting system does not have CAC property.

As shown in [3, 9, 16, 24], the Lax pair (B.1) with (B.2b) of the PΔ\DeltaE (1.1) can be constructed from the equations (2.7b) and (2.7c) by setting

xl,m=Fl,mGl,m,Φl,m=(Fl,mGl,m).x_{l,m}=\dfrac{F_{l,m}}{G_{l,m}},\quad\Phi_{l,m}=\begin{pmatrix}F_{l,m}\\ G_{l,m}\end{pmatrix}. (2.10)
Remark 2.4.

Similar to the Lax pair of the PΔ\DeltaE (1.1), we can also construct a Lax pair of the PΔ\DeltaE (1.2) from the equations (2.7b) and (2.7c) by setting

ul,m=Fl,mGl,m,Φl,m=(Fl,mGl,m).u_{l,m}=\dfrac{F_{l,m}}{G_{l,m}},\quad\Phi_{l,m}=\begin{pmatrix}F_{l,m}\\ G_{l,m}\end{pmatrix}. (2.11)

However, we still have not figured out how to introduce a spectral parameter in this case.

2.3. A CAC-type BT from the PΔ\DeltaE (1.1) to the PΔ\DeltaE (1.2) (II)

In this subsection, using the CABC-type BTs (1.7) and (2.1), we obtain a CAC-type BT from the PΔ\DeltaE (1.1) to the PΔ\DeltaE (1.2). Since the process for demonstrating the result is the same as that for the CABC-type BTs (1.5) and (1.7) discussed in §2.2, we omit detailed arguments.

Theorem 2.5.

The following system has the CAC property but does not have the tetrahedron property:

u¯~u=βm+1αlu~βmαl+1u¯,\displaystyle\widetilde{\overline{u}}-u=\dfrac{\beta_{m+1}-\alpha_{l}}{\widetilde{u}}-\dfrac{\beta_{m}-\alpha_{l+1}}{\overline{u}}, (2.12a)
u~u+(qmplx)(plqm+1x~)x=0,\displaystyle\widetilde{u}u+\dfrac{(q_{m}-p_{l}x)(p_{l}-q_{m+1}\widetilde{x})}{x}=0, (2.12b)
u¯u+(qmplx)(pl+1qmx¯)x¯=0,\displaystyle\overline{u}u+\dfrac{(q_{m}-p_{l}x)(p_{l+1}-q_{m}\overline{x})}{\overline{x}}=0, (2.12c)
x¯~x=(pl+1qmx¯qmpl+1x¯)(qm+1plx~plqm+1x~),\displaystyle\dfrac{\,\widetilde{\overline{x}}\,}{x}=\left(\dfrac{p_{l+1}-q_{m}\overline{x}}{q_{m}-p_{l+1}\overline{x}}\right)\left(\dfrac{q_{m+1}-p_{l}\widetilde{x}}{p_{l}-q_{m+1}\widetilde{x}}\right), (2.12d)

where plp_{l} and qmq_{m} are given by Equation (1.8). Since the equations (2.12a) and (2.12d) are equal to the PΔ\DeltaEs (1.1) and (1.2), respectively, we can also say that the pair of equations (2.12b) and (2.12c) is a CAC-type BT from the PΔ\DeltaE (1.1) to the PΔ\DeltaE (1.2).

Proof.

Consider the cube in Figure 2.1. In the three different ways (see Appendix A.1), x¯~\widetilde{\overline{x}} can be uniquely represented by the initial values {u,u¯,u~,x}\{u,\overline{u},\widetilde{u},x\} as

x¯~=pl+1u¯((αlβm+1)(qmplx)+pluu~x)qm+1u~((αl+1βm)(qmplx)+qmuu¯).\widetilde{\overline{x}}=\dfrac{p_{l+1}\overline{u}\,\Big{(}(\alpha_{l}-\beta_{m+1})(q_{m}-p_{l}x)+p_{l}u\widetilde{u}x\Big{)}}{q_{m+1}\widetilde{u}\,\Big{(}(\alpha_{l+1}-\beta_{m})(q_{m}-p_{l}x)+q_{m}u\overline{u}\Big{)}}. (2.13)

Since x¯~\widetilde{\overline{x}} depends on uu, the tetrahedron property does not hold. Therefore, we have completed the proof. ∎

Setting

xl,m=Fl,mGl,m,Φl,m=(Fl,mGl,m),x_{l,m}=\dfrac{F_{l,m}}{G_{l,m}},\quad\Phi_{l,m}=\begin{pmatrix}F_{l,m}\\ G_{l,m}\end{pmatrix}, (2.14)

from the equations (2.12b) and (2.12c), we obtain the Lax pair (B.1) with (B.2c) of the PΔ\DeltaE (1.1).

Remark 2.6.

Similar to the Lax pair of the PΔ\DeltaE (1.1), we can also construct a Lax pair of the PΔ\DeltaE (1.2) from the equations (2.12b) and (2.12c) by setting

ul,m=Fl,mGl,m,Φl,m=(Fl,mGl,m).u_{l,m}=\dfrac{F_{l,m}}{G_{l,m}},\quad\Phi_{l,m}=\begin{pmatrix}F_{l,m}\\ G_{l,m}\end{pmatrix}. (2.15)

However, in this case, how to introduce a spectral parameter is still unsolved at this point.

2.4. A CAC-type auto-BT of the PΔ\DeltaE (1.2)

In this subsection, using two CAC-systems (2.7) and (2.12), we obtain a CAC-system, which implies a CAC-type auto-BT of the PΔ\DeltaE (1.2).

Here, let us use the following CAC-system:

u¯~u=βm+1αlu~βmαl+1u¯,\displaystyle\widetilde{\overline{u}}-u=\dfrac{\beta_{m+1}-\alpha_{l}}{\widetilde{u}}-\dfrac{\beta_{m}-\alpha_{l+1}}{\overline{u}}, (2.16a)
u~u+(PlQmX)(Qm+1PlX~)X~=0,\displaystyle\widetilde{u}u+\dfrac{(P_{l}-Q_{m}X)(Q_{m+1}-P_{l}\widetilde{X})}{\widetilde{X}}=0, (2.16b)
u¯u+(PlQmX)(QmPl+1X¯)X=0,\displaystyle\overline{u}u+\dfrac{(P_{l}-Q_{m}X)(Q_{m}-P_{l+1}\overline{X})}{X}=0, (2.16c)
X¯~X=(Pl+1QmX¯QmPl+1X¯)(Qm+1PlX~PlQm+1X~),\displaystyle\dfrac{\,\widetilde{\overline{X}}\,}{X}=\left(\dfrac{P_{l+1}-Q_{m}\overline{X}}{Q_{m}-P_{l+1}\overline{X}}\right)\left(\dfrac{Q_{m+1}-P_{l}\widetilde{X}}{P_{l}-Q_{m+1}\widetilde{X}}\right), (2.16d)

where X=Xl,mX=X_{l,m} and

Pl=αlγλ,Qm=βmγλ,P_{l}=\sqrt{\alpha_{l}-\gamma-\lambda},\quad Q_{m}=\sqrt{\beta_{m}-\gamma-\lambda}\,, (2.17)

instead of the CAC-system (2.7) with the variable xl,mx_{l,m} and the parameters plp_{l}, qmq_{m} in the system (2.7) replaced by the variable Xl,mX_{l,m} and the parameters PlP_{l}, QmQ_{m}, respectively. Eliminating u~u\widetilde{u}u from the equations (2.12b) and (2.16b) and u¯u\overline{u}u from the equations (2.12c) and (2.16c), we obtain

(qmplx)(plqm+1x~)x=(PlQmX)(Qm+1PlX~)X~,\displaystyle\dfrac{(q_{m}-p_{l}x)(p_{l}-q_{m+1}\widetilde{x})}{x}=\dfrac{(P_{l}-Q_{m}X)(Q_{m+1}-P_{l}\widetilde{X})}{\widetilde{X}}, (2.18a)
(qmplx)(pl+1qmx¯)x¯=(PlQmX)(QmPl+1X¯)X,\displaystyle\dfrac{(q_{m}-p_{l}x)(p_{l+1}-q_{m}\overline{x})}{\overline{x}}=\dfrac{(P_{l}-Q_{m}X)(Q_{m}-P_{l+1}\overline{X})}{X}, (2.18b)

respectively. The equations above, together with the equations (2.12d) and (2.16d), lead to the following theorem.

Theorem 2.7.

The following system has the CAC property but does not have the tetrahedron property:

x¯~x=(pl+1qmx¯qmpl+1x¯)(qm+1plx~plqm+1x~),\displaystyle\dfrac{\,\widetilde{\overline{x}}\,}{x}=\left(\dfrac{p_{l+1}-q_{m}\overline{x}}{q_{m}-p_{l+1}\overline{x}}\right)\left(\dfrac{q_{m+1}-p_{l}\widetilde{x}}{p_{l}-q_{m+1}\widetilde{x}}\right), (2.19a)
(qmplx)(plqm+1x~)x=(PlQmX)(Qm+1PlX~)X~,\displaystyle\dfrac{(q_{m}-p_{l}x)(p_{l}-q_{m+1}\widetilde{x})}{x}=\dfrac{(P_{l}-Q_{m}X)(Q_{m+1}-P_{l}\widetilde{X})}{\widetilde{X}}, (2.19b)
(qmplx)(pl+1qmx¯)x¯=(PlQmX)(QmPl+1X¯)X,\displaystyle\dfrac{(q_{m}-p_{l}x)(p_{l+1}-q_{m}\overline{x})}{\overline{x}}=\dfrac{(P_{l}-Q_{m}X)(Q_{m}-P_{l+1}\overline{X})}{X}, (2.19c)
X¯~X=(Pl+1QmX¯QmPl+1X¯)(Qm+1PlX~PlQm+1X~),\displaystyle\dfrac{\,\widetilde{\overline{X}}\,}{X}=\left(\dfrac{P_{l+1}-Q_{m}\overline{X}}{Q_{m}-P_{l+1}\overline{X}}\right)\left(\dfrac{Q_{m+1}-P_{l}\widetilde{X}}{P_{l}-Q_{m+1}\widetilde{X}}\right), (2.19d)

where

Pl=pl2λ,Qm=qm2λ.P_{l}=\sqrt{{p_{l}}^{2}-\lambda},\quad Q_{m}=\sqrt{{q_{m}}^{2}-\lambda}\,. (2.20)

Since both of the equations (2.19a) and (2.19d) are equal to the PΔ\DeltaE (1.2), we can also say that the pair of equations (2.19b) and (2.19c) is a CAC-type auto-BT of the PΔ\DeltaE (1.2).

Proof.

Consider the cube in Figure 2.2. In the three different ways (see Appendix A.1), X¯~\widetilde{\overline{X}} can be uniquely represented by the initial values {x,x¯,x~,X}\{x,\overline{x},\widetilde{x},X\} as

X¯~=Qm+1(plqm+1x~)((Pl+12Qm2)(PlQmX)x¯+Qm(pl+1qmx¯)(qmplx)X)Pl+1(pl+1qmx¯)((Pl2Qm+12)(PlQmX)x+Pl(qmplx)(plqm+1x~)).\widetilde{\overline{X}}=\dfrac{Q_{m+1}(p_{l}-q_{m+1}\widetilde{x})\Big{(}({P_{l+1}}^{2}-{Q_{m}}^{2})(P_{l}-Q_{m}X)\overline{x}+Q_{m}(p_{l+1}-q_{m}\overline{x})(q_{m}-p_{l}x)X\Big{)}}{P_{l+1}(p_{l+1}-q_{m}\overline{x})\Big{(}({P_{l}}^{2}-{Q_{m+1}}^{2})(P_{l}-Q_{m}X)x+P_{l}(q_{m}-p_{l}x)(p_{l}-q_{m+1}\widetilde{x})\Big{)}}. (2.21)

Since X¯~\widetilde{\overline{X}} depends on xx, the tetrahedron property does not hold. Moreover, the condition (2.20) can be obtained from the conditions (1.8) and (2.17). Therefore, we have completed the proof. ∎

Refer to caption
Figure 2.2. A cube for the CAC-system (2.19). See the caption of Figure 2.1 for explaining the placement of variables and parameters in the figure.

Setting

Xl,m=Fl,mGl,m,Ψl,m=(Fl,mGl,m),X_{l,m}=\dfrac{F_{l,m}}{G_{l,m}},\quad\Psi_{l,m}=\begin{pmatrix}F_{l,m}\\ G_{l,m}\end{pmatrix}, (2.22)

from the equations (2.19b) and (2.19c), we obtain the Lax pair (B.4) of the PΔ\DeltaE (1.2).

Remark 2.8.

Since Equation (2.19d) is also the lsG equation, setting

xl,m=Fl,mGl,m,Ψl,m=(Fl,mGl,m),x_{l,m}=\dfrac{F_{l,m}}{G_{l,m}},\quad\Psi_{l,m}=\begin{pmatrix}F_{l,m}\\ G_{l,m}\end{pmatrix}, (2.23)

we can construct another Lax pair of the lsG equation. However, since the system (2.19) and the condition (2.20) are invariant under the transformation

xl,m1Xl,m,plPl,qmQm,λλ,x_{l,m}\leftrightarrow\dfrac{1}{X_{l,m}},\quad p_{l}\leftrightarrow P_{l},\quad q_{m}\leftrightarrow Q_{m},\quad\lambda\leftrightarrow-\lambda, (2.24)

the resulting Lax pair is essentially the same as the Lax pair (B.4).

3. Integrable systems on the 3-dimensional integer lattice

In general, if a quad-equation has a CAC-type auto-BT, it can be easily extended to a system of PΔ\DeltaEs on the 3-dimensional integer lattice 3\mathbb{Z}^{3}. (See [1, 9, 20].) However, in the case of non-auto-BTs, a little ingenuity is required. (See [4].) In this section, we show how quad-equations that have CABC-type BTs, which are non-auto-BTs, can be extended to a system of PΔ\DeltaEs on the lattice 3\mathbb{Z}^{3}. Moreover, we also show how the resulting system can be reduced to systems on the lattice 3\mathbb{Z}^{3} associated with the CAC-type BTs.

3.1. A system of PΔ\DeltaEs on the lattice 3\mathbb{Z}^{3} associated with the CABC property

In this subsection, using the three CABC-systems (1.5), (1.7) and (2.1), we construct an integrable system of PΔ\DeltaEs on the lattice 3\mathbb{Z}^{3}. We also show the τ\tau-function of the resulting system.

Remark 3.1.

Since Equation (1.5a) can be obtained from the equations (1.5b) and (1.5c), we here omit Equation (1.5a) when using the CABC-system (1.5) for simplicity. For the same reason, the equations (1.7a) and (2.1a) are also omitted when using the CABC-systems (1.7) and (2.1), respectively.

Let us denote the cubes for the three CABC-systems (1.5), (1.7) and (2.1) by the symbols Cl,m(1)C^{(1)}_{l,m}, Cl,m(2)C^{(2)}_{l,m} and Cl,m(3)C^{(3)}_{l,m}, respectively. (See §A.2 for the correspondence between a CABC-system and a cube.) Overlap the cubes as in Figure 3.1 (Left), in the order Cl,m(1)C^{(1)}_{l,m}, Cl,m(2)C^{(2)}_{l,m}, Cl,m(2)C^{(2)}_{l,m}, Cl,m(3)C^{(3)}_{l,m}, so that the facing surfaces have the same equations. We here also denote the resulting rectangular constructed by the set of four stacked cubes by the symbol Cl,m(1223)C^{(1223)}_{l,m}. Because both of the face-equations at the bottom and top of Cl,m(1223)C^{(1223)}_{l,m} are the dKdV equation, it can be placed repeatedly up and down as shown in Figure 3.1 (Right). Note that since the parameter γ\gamma is not included in the dKdV equation, a different γ\gamma can be chosen for each stack of Cl,m(1223)C^{(1223)}_{l,m}. This repeated placement yields the following system of PΔ\DeltaEs on the lattice 3\mathbb{Z}^{3}:

u¯~u=βm+1αlu~βmαl+1u¯,\displaystyle\widetilde{\overline{u}}-u=\dfrac{\beta_{m+1}-\alpha_{l}}{\widetilde{u}}-\dfrac{\beta_{m}-\alpha_{l+1}}{\overline{u}}, (3.1a)
x¯~x=(pl+1,nqm,nx¯qm,npl+1,nx¯)(qm+1,npl,nx~pl,nqm+1,nx~),\displaystyle\dfrac{\,\widetilde{\overline{x}}\,}{x}=\left(\dfrac{p_{l+1,n}-q_{m,n}\overline{x}}{q_{m,n}-p_{l+1,n}\overline{x}}\right)\left(\dfrac{q_{m+1,n}-p_{l,n}\widetilde{x}}{p_{l,n}-q_{m+1,n}\widetilde{x}}\right), (3.1b)
(uv)(βmαlu+v~)βm+γn=0,\displaystyle(u-v)\left(\dfrac{\beta_{m}-\alpha_{l}}{u}+\widetilde{v}\right)-\beta_{m}+\gamma_{n}=0, (3.1c)
βmαluγnαlvu¯+v¯=0,\displaystyle\dfrac{\beta_{m}-\alpha_{l}}{u}-\dfrac{\gamma_{n}-\alpha_{l}}{v}-\overline{u}+\overline{v}=0, (3.1d)
vv~pl,n2x(pl,nqm,nx)pl,nxqm,n=0,\displaystyle\dfrac{v\widetilde{v}}{{p_{l,n}}^{2}}-\dfrac{x(p_{l,n}-q_{m,n}x)}{p_{l,n}x-q_{m,n}}=0, (3.1e)
vv¯pl,npl+1,nx¯(pl,nqm,nx)pl,nxqm,n=0,\displaystyle\dfrac{v\overline{v}}{p_{l,n}p_{l+1,n}}-\dfrac{\overline{x}(p_{l,n}-q_{m,n}x)}{p_{l,n}x-q_{m,n}}=0, (3.1f)
(u^v~)(βmαlu^+v)βm+γn=0,\displaystyle(\widehat{u}-\widetilde{v})\left(\dfrac{\beta_{m}-\alpha_{l}}{\widehat{u}}+v\right)-\beta_{m}+\gamma_{n}=0, (3.1g)
βmαlu^γnαl+1v¯u¯^+v=0,\displaystyle\dfrac{\beta_{m}-\alpha_{l}}{\widehat{u}}-\dfrac{\gamma_{n}-\alpha_{l+1}}{\overline{v}}-\widehat{\overline{u}}+v=0, (3.1h)

where u=ul,m,nu=u_{l,m,n}, v=vl,m,nv=v_{l,m,n}, x=xl,m,nx=x_{l,m,n} and

pl,n=αlγn,qm,n=βmγn.p_{l,n}=\sqrt{\alpha_{l}-\gamma_{n}},\quad q_{m,n}=\sqrt{\beta_{m}-\gamma_{n}}. (3.2)
Remark 3.2.

The correspondence between the cubes Cl,m(1)C^{(1)}_{l,m}, Cl,m(2)C^{(2)}_{l,m}, Cl,m(3)C^{(3)}_{l,m} and the equations in the system (3.1) is as follows.

Cl,m(1)C^{(1)}_{l,m} : (3.1a), (3.1c), (3.1d)
Cl,m(2)C^{(2)}_{l,m} : (3.1b), (3.1e), (3.1f)
Cl,m(3)C^{(3)}_{l,m} : (3.1a)n→n+1, (3.1g), (3.1h)
Refer to caption
Refer to caption
Figure 3.1. (Left): The rectangular Cl,m(1223)C^{(1223)}_{l,m} constructed by the four cubes. Since the uu-variable in the bottom cube Cl,m(1)C^{(1)}_{l,m} and that in the top cube Cl,m(3)C^{(3)}_{l,m} can be treated as different, we used the symbol UU for the uu-variable of Cl,m(3)C^{(3)}_{l,m} to differentiate them.
(Right): Two rectangulars Cl,m(1223)C^{(1223)}_{l,m}. Since we choose a different parameter γ\gamma for each Cl,m(1223)C^{(1223)}_{l,m}, we place the parameters γn\gamma_{n} on the vertical line segments of the rectangles.

The τ\tau-function is one of the most important objects in the theory of integrable systems. We show the τ\tau-function of the system (3.1) in the following theorem.

Theorem 3.3.

Let τ=τl,m,n\tau=\tau_{l,m,n} be the solution of the following bilinear equations:

pl,nτ~^τ¯qm,nτ~τ¯^+ττ¯~^=0,\displaystyle p_{l,n}\widehat{\widetilde{\tau}}\,\overline{\tau}-q_{m,n}\widetilde{\tau}\,\widehat{\overline{\tau}}+\tau\,\widehat{\widetilde{\overline{\tau}}}=0, (3.3a)
qm,nτ~^τ¯pl,nτ~τ¯^τ^τ¯~=0.\displaystyle q_{m,n}\widehat{\widetilde{\tau}}\,\overline{\tau}-p_{l,n}\widetilde{\tau}\,\widehat{\overline{\tau}}-\widehat{\tau}\,\widetilde{\overline{\tau}}=0. (3.3b)

Then, the τ\tau-function of the system (3.1) is given as

u=ττ¯~τ¯τ~,v=pl,nττ¯^τ¯τ^,x=τ¯^τ~τ¯τ~^.u=-\dfrac{\,\tau\,\widetilde{\overline{\tau}}\,}{\overline{\tau}\,\widetilde{\tau}},\quad v=p_{l,n}\,\dfrac{\,\tau\,\widehat{\overline{\tau}}\,}{\overline{\tau}\,\widehat{\tau}},\quad x=\dfrac{\,\widehat{\overline{\tau}}\,\widetilde{\tau}\,}{\overline{\tau}\,\widehat{\widetilde{\tau}}}. (3.4)
Proof.

The proof is given in Appendix C. ∎

3.2. Systems of PΔ\DeltaEs on the lattice 3\mathbb{Z}^{3} associated with the CAC property

In this subsection, we derive integrable systems on the lattice 3\mathbb{Z}^{3} associated with the CAC property from the system (3.1).

As in §2.2 and §2.3, eliminating the vv-variable from the system (3.1) we obtain the following system of PΔ\DeltaEs depending only on the uu- and xx-variables:

u¯~u=βm+1αlu~βmαl+1u¯,\displaystyle\widetilde{\overline{u}}-u=\dfrac{\beta_{m+1}-\alpha_{l}}{\widetilde{u}}-\dfrac{\beta_{m}-\alpha_{l+1}}{\overline{u}}, (3.5a)
x¯~x=(pl+1,nqm,nx¯qm,npl+1,nx¯)(qm+1,npl,nx~pl,nqm+1,nx~),\displaystyle\dfrac{\,\widetilde{\overline{x}}\,}{x}=\left(\dfrac{p_{l+1,n}-q_{m,n}\overline{x}}{q_{m,n}-p_{l+1,n}\overline{x}}\right)\left(\dfrac{q_{m+1,n}-p_{l,n}\widetilde{x}}{p_{l,n}-q_{m+1,n}\widetilde{x}}\right), (3.5b)
u~u+(pl,nqm,nx)(qm+1,npl,nx~)x~=0,\displaystyle\widetilde{u}\,u+\dfrac{(p_{l,n}-q_{m,n}x)(q_{m+1,n}-p_{l,n}\widetilde{x})}{\widetilde{x}}=0, (3.5c)
u¯u+(pl,nqm,nx)(qm,npl+1,nx¯)x=0,\displaystyle\overline{u}\,u+\dfrac{(p_{l,n}-q_{m,n}x)(q_{m,n}-p_{l+1,n}\overline{x})}{x}=0, (3.5d)
u~^u^+(qm,npl,nx)(pl,nqm+1,nx~)x=0,\displaystyle\widehat{\widetilde{u}}~{}\widehat{u}+\dfrac{(q_{m,n}-p_{l,n}x)(p_{l,n}-q_{m+1,n}\widetilde{x})}{x}=0, (3.5e)
u¯^u^+(qm,npl,nx)(pl+1,nqm,nx¯)x¯=0.\displaystyle\widehat{\overline{u}}~{}\widehat{u}+\dfrac{(q_{m,n}-p_{l,n}x)(p_{l+1,n}-q_{m,n}\overline{x})}{\overline{x}}=0. (3.5f)

Indeed, eliminating v~\widetilde{v} from the equations (3.1c) and (3.1e), we obtain

v=pl,nuxpl,nxqm,n.v=\dfrac{p_{l,n}ux}{p_{l,n}x-q_{m,n}}. (3.6)

Then, eliminating vv-variable from the equations (3.1c) and (3.1d) by using the equation above, we obtain the equations (3.5c) and (3.5d), respectively. Moreover, eliminating v~\widetilde{v} from the equations (3.1e) and (3.1g), we obtain

v=pl,n(pl,nqm,nx)u^,v=\dfrac{p_{l,n}(p_{l,n}-q_{m,n}x)}{\widehat{u}}, (3.7)

and then eliminating vv-variable from the equations (3.1g) and (3.1h) by using the equation above, we obtain the equations (3.5e) and (3.5f), respectively.

Remark 3.4.

As shown in Figure 3.2 (Left), the system (3.5) can be thought of as a superposition of the CAC-systems (2.7) and (2.12) alternating in the nn-direction.

Refer to caption
Refer to caption
Figure 3.2. (Left): Cubes for the CAC-system (3.5). The system of equations (3.5a), (3.5b), (3.5c) and (3.5d) correspond to the lower cube, and the system of equations (3.5a)n→n+1, (3.5b), (3.5e) and (3.5f) correspond to the upper cube.
(Right): A cube for the CAC-system (3.8).

As in §2.4, we eliminate u~^u^\widehat{\widetilde{u}}\,\widehat{u} from the equations (3.5c)n→n+1 and (3.5e). Also, we eliminate u¯^u^\widehat{\overline{u}}\,\widehat{u} from the equations (3.5d)n→n+1 and (3.5f). The resulting equations, together with Equation (3.5b), lead to the following system of PΔ\DeltaEs depending only on the xx-variable:

x¯~x=(pl+1,nqm,nx¯qm,npl+1,nx¯)(qm+1,npl,nx~pl,nqm+1,nx~),\displaystyle\dfrac{\,\widetilde{\overline{x}}\,}{x}=\left(\dfrac{p_{l+1,n}-q_{m,n}\overline{x}}{q_{m,n}-p_{l+1,n}\overline{x}}\right)\left(\dfrac{q_{m+1,n}-p_{l,n}\widetilde{x}}{p_{l,n}-q_{m+1,n}\widetilde{x}}\right), (3.8a)
(qm,npl,nx)(pl,nqm+1,nx~)x=(pl,n+1qm,n+1x^)(qm+1,n+1pl,n+1x~^)x~^,\displaystyle\dfrac{(q_{m,n}-p_{l,n}x)(p_{l,n}-q_{m+1,n}\widetilde{x})}{x}=\dfrac{(p_{l,n+1}-q_{m,n+1}\widehat{x})(q_{m+1,n+1}-p_{l,n+1}\widehat{\widetilde{x}})}{\widehat{\widetilde{x}}}, (3.8b)
(qm,npl,nx)(pl+1,nqm,nx¯)x¯=(pl,n+1qm,n+1x^)(qm,n+1pl+1,n+1x¯^)x^.\displaystyle\dfrac{(q_{m,n}-p_{l,n}x)(p_{l+1,n}-q_{m,n}\overline{x})}{\overline{x}}=\dfrac{(p_{l,n+1}-q_{m,n+1}\widehat{x})(q_{m,n+1}-p_{l+1,n+1}\widehat{\overline{x}})}{\widehat{x}}. (3.8c)
Remark 3.5.
  1. (1)(\mathrm{1})

    As shown in Figure 3.2 (Right), the system (3.8) can be thought of as a superposition of the CAC-system (2.19) in the nn-direction.

  2. (2)(\mathrm{2})

    Both equations (3.8b) and (3.8c) are also the lsG equation. Indeed, by setting

    ym,n=qm,npl,nxl,m,nαlβmxl,m,n,am=αlβm,bn=pl,n,y_{m,n}=\dfrac{q_{m,n}-p_{l,n}x_{l,m,n}}{\sqrt{\alpha_{l}-\beta_{m}}\,x_{l,m,n}},\quad a_{m}=-\sqrt{\alpha_{l}-\beta_{m}},\quad b_{n}=p_{l,n}, (3.9)

    Equation (3.8b) can be rewritten as the following standard form of the lsG equation in (m,n)(m,n)-direction:

    ym+1,n+1ym,n=(bn+1amym,n+1ambn+1ym,n+1)(am+1bnym+1,nbnam+1ym+1,n),\dfrac{y_{m+1,n+1}}{y_{m,n}}=\left(\dfrac{b_{n+1}-a_{m}y_{m,n+1}}{a_{m}-b_{n+1}y_{m,n+1}}\right)\left(\dfrac{a_{m+1}-b_{n}y_{m+1,n}}{b_{n}-a_{m+1}y_{m+1,n}}\right), (3.10)

    while by setting

    zl,n=pl,nxl,m,nqm,nβmαl,cl=βmαl,dn=qm,n,z_{l,n}=\dfrac{p_{l,n}x_{l,m,n}-q_{m,n}}{\sqrt{\beta_{m}-\alpha_{l}}},\quad c_{l}=-\sqrt{\beta_{m}-\alpha_{l}},\quad d_{n}=q_{m,n}, (3.11)

    Equation (3.8c) can also be rewritten as the following standard form of the lsG equation in (l,n)(l,n)-direction:

    zl+1,n+1zl,n=(dn+1clzl,n+1cldn+1zl,n+1)(cl+1dnzl+1,ndncl+1zl+1,n).\dfrac{z_{l+1,n+1}}{z_{l,n}}=\left(\dfrac{d_{n+1}-c_{l}z_{l,n+1}}{c_{l}-d_{n+1}z_{l,n+1}}\right)\left(\dfrac{c_{l+1}-d_{n}z_{l+1,n}}{d_{n}-c_{l+1}z_{l+1,n}}\right). (3.12)

4. Concluding remarks

In this paper, using the CABC properties of the dKdV equation (1.1) and the lsG equation (1.2), we have shown that they have the CAC property. Using the CABC and CAC properties, we also obtained their new Lax pairs. Moreover, we have constructed the integrable system (3.1) on the lattice 3\mathbb{Z}^{3} by stacking the cubes associated with the CABC property in the third direction and have clarified the structure of its τ\tau-function.

In addition to the dKdV and lsG equations, there exist other integrable PΔ\DeltaEs for which we still do not know whether they have the CAC property. We have shown that it is useful to consider the CABC property to clarify the CAC property of PΔ\DeltaEs. Many studies have been done to find PΔ\DeltaEs which have the CAC property (see §1). However, the dKdV and lsG equations have the CAC property but are not included in these previous studies. For the reasons above, classifying PΔ\DeltaEs with the CABC property may be important for clarifying the CAC property of known integrable PΔ\DeltaEs. Moreover, it may lead to finding new integrable PΔ\DeltaEs. This should be a subject of future research.

Acknowledgment

The author would like to thank Dr. Pavlos Kassotakis and Dr. Yang Shi for fruitful discussions. This research was supported by a JSPS KAKENHI Grant Number JP19K14559. The author used Wolfram Mathematica for checking computations.

Appendix A The CAC and CABC properties

In this appendix, we explain the consistency around a cube (CAC) property of a system of PΔ\DeltaEs in a particular form. In addition, we also explain the consistency around a broken cube (CABC) property of a system of PΔ\DeltaEs in a particular form. For more general definitions, see [3, 9, 16, 24] for the CAC property and see [11, 15] for the CABC property.

Remark A.1.

For simplicity, we omit the parameters of the equations in this appendix. Moreover, the equations in this appendix are not limited to autonomous type. For example, the relation between Q(u,u¯,u~,u¯~)=0Q(u,\overline{u},\widetilde{u},\widetilde{\overline{u}})=0 and Q(u¯,u¯¯,u¯~,u¯¯~)=0Q(\overline{u},\overline{\overline{u}},\widetilde{\overline{u}},\widetilde{\overline{\overline{u}}})=0 is not simple replacement of the corresponding uu-variables; applying ll+1l\to l+1 to Q(u,u¯,u~,u¯~)=0Q(u,\overline{u},\widetilde{u},\widetilde{\overline{u}})=0 gives Q(u¯,u¯¯,u¯~,u¯¯~)=0Q(\overline{u},\overline{\overline{u}},\widetilde{\overline{u}},\widetilde{\overline{\overline{u}}})=0. See §2 for examples involving parameters.

A.1. The CAC property

Let us consider the following system of PΔ\DeltaEs:

A(u,u¯,u~,u¯~)=0,\displaystyle A(u,\overline{u},\widetilde{u},\widetilde{\overline{u}})=0, (A.1a)
B(u,u~,v,v~)=0,\displaystyle B(u,\widetilde{u},v,\widetilde{v})=0, (A.1b)
C(u,u¯,v,v¯)=0,\displaystyle C(u,\overline{u},v,\overline{v})=0, (A.1c)
A(v,v¯,v~,v¯~)=0,\displaystyle A^{\prime}(v,\overline{v},\widetilde{v},\widetilde{\overline{v}})=0, (A.1d)

where u=ul,mu=u_{l,m} and v=vl,mv=v_{l,m}. Here, the functions AA, BB, CC and AA^{\prime} are irreducible multilinear polynomials in four variables. Consider the following sublattice of the lattice 3\mathbb{Z}^{3}:

{(l,m,0)3|l,m}{(l,m,1)3|l,m}\left\{\left.(l,m,0)\in\mathbb{Z}^{3}~{}\right|~{}l,m\in\mathbb{Z}\right\}\cup\left\{\left.(l,m,1)\in\mathbb{Z}^{3}~{}\right|~{}l,m\in\mathbb{Z}\right\} (A.2)

and assign the uu- and vv-variables on the vertices of the sublattice by the following correspondences:

(l,m,0)ul,m,(l,m,1)vl,m.(l,m,0)\,\leftrightarrow\,u_{l,m},\qquad(l,m,1)\,\leftrightarrow\,v_{l,m}. (A.3)

Then, focusing on the cube given by the 8 points

(l,m,0),(l+1,m,0),(l,m+1,0),(l+1,m+1,0),\displaystyle(l,m,0),\quad(l+1,m,0),\quad(l,m+1,0),\quad(l+1,m+1,0),
(l,m,1),(l+1,m,1),(l,m+1,1),(l+1,m+1,1),\displaystyle(l,m,1),\quad(l+1,m,1),\quad(l,m+1,1),\quad(l+1,m+1,1), (A.4)

from the system (A.1), we obtain the following face-equations of the cube (see Figure A.1):

𝒜=A(u,u¯,u~,u¯~)=0,\displaystyle{\mathcal{A}}=A(u,\overline{u},\widetilde{u},\widetilde{\overline{u}})=0, (A.5a)
𝒜=A(v,v¯,v~,v¯~)=0,\displaystyle{\mathcal{A}}^{\prime}=A^{\prime}(v,\overline{v},\widetilde{v},\widetilde{\overline{v}})=0, (A.5b)
=B(u,u~,v,v~)=0,\displaystyle{\mathcal{B}}=B(u,\widetilde{u},v,\widetilde{v})=0, (A.5c)
¯=B(u¯,u¯~,v¯,v¯~)=0,\displaystyle\overline{{\mathcal{B}}}=B(\overline{u},\widetilde{\overline{u}},\overline{v},\widetilde{\overline{v}})=0, (A.5d)
𝒞=C(u,u¯,v,v¯)=0,\displaystyle{\mathcal{C}}=C(u,\overline{u},v,\overline{v})=0, (A.5e)
𝒞~=C(u~,u¯~,v~,v¯~)=0.\displaystyle\widetilde{{\mathcal{C}}}=C(\widetilde{u},\widetilde{\overline{u}},\widetilde{v},\widetilde{\overline{v}})=0. (A.5f)

The CAC and tetrahedron properties are defined as follows.

  1. (1)(\mathrm{1})

    There are following three ways to calculate v¯~\widetilde{\overline{v}} by using all the equations in the system (A.5) with {u,u¯,u~,v}\{u,\overline{u},\widetilde{u},v\} as initial values.

    1. (a)

      Express v¯~\widetilde{\overline{v}} as a rational function in terms of {v,v¯,v~}\{v,\overline{v},\widetilde{v}\} by using Equation (A.5b). Then, eliminate v~\widetilde{v} by using Equation (A.5c) and v¯\overline{v} by using Equation (A.5e) from it.

    2. (b)

      Express v¯~\widetilde{\overline{v}} as a rational function in terms of {u¯,u¯~,v¯}\{\overline{u},\widetilde{\overline{u}},\overline{v}\} by using Equation (A.5d). Then, eliminate u¯~\widetilde{\overline{u}} by using Equation (A.5a) and v¯\overline{v} by using Equation (A.5e) from it.

    3. (c)

      Express v¯~\widetilde{\overline{v}} as a rational function in terms of {u~,u¯~,v~}\{\widetilde{u},\widetilde{\overline{u}},\widetilde{v}\} by using Equation (A.5f). Then, eliminate u¯~\widetilde{\overline{u}} by using Equation (A.5a) and v~\widetilde{v} by using Equation (A.5c) from it.

    When v¯~\widetilde{\overline{v}} is uniquely determined as a rational function with the initial values {u,u¯,u~,v}\{u,\overline{u},\widetilde{u},v\}, then the system (A.1) is said to have the CAC property or said to be a CAC-system.

  2. (2)(\mathrm{2})

    Let the system (A.1) be a CAC-system. When v¯~\widetilde{\overline{v}} is represented as a rational function with the initial values {u,u¯,u~,v}\{u,\overline{u},\widetilde{u},v\}, the system (A.1) is said to have the tetrahedron property if the rational function does not depend on uu.

  3. (3)(\mathrm{3})

    We say that the PΔ\DeltaE (A.1a) (or the PΔ\DeltaE (A.1d)) has the CAC property if the system (A.1) has the CAC property. The pair of equations (A.1b) and (A.1c) is then referred to as a CAC-type BT from the PΔ\DeltaE (A.1a) to the PΔ\DeltaE (A.1d) (or from the PΔ\DeltaE (A.1d) to the PΔ\DeltaE (A.1a)). When the equations (A.1a) and (A.1d) are the same PΔ\DeltaE, the CAC-type BT (A.1b) and (A.1c) is specifically called a CAC-type auto-BT of the PΔ\DeltaE (A.1a).

Remark A.2.

Even if Equation (A.1a) has the CAC property, it is not necessarily integrable[8]. Even if a Lax pair of Equation (A.1a) is constructed from the CAC-type BT (A.1b) and (A.1c) by the method in [3, 9, 16, 24], it may be fake [6] or weak [10]. From this perspective, in [8], the CAC-type BT (A.1b) and (A.1c) are characterized as follows. Calculate v¯~\widetilde{\overline{v}} in the following two ways.

  1. (a)

    Express v¯~\widetilde{\overline{v}} as a rational function in terms of {u¯,u¯~,v¯}\{\overline{u},\widetilde{\overline{u}},\overline{v}\} by using Equation (A.5d). Then, eliminate v¯\overline{v} from it by using Equation (A.5e). Finally, let v¯~=f1(u,u¯,u¯~,v)\widetilde{\overline{v}}=f_{1}(u,\overline{u},\widetilde{\overline{u}},v) be the resulting relation.

  2. (b)

    Express v¯~\widetilde{\overline{v}} as a rational function in terms of {u~,u¯~,v~}\{\widetilde{u},\widetilde{\overline{u}},\widetilde{v}\} by using Equation (A.5f). Then, eliminate v~\widetilde{v} from it by using Equation (A.5c). Finally, let v¯~=f2(u,u~,u¯~,v)\widetilde{\overline{v}}=f_{2}(u,\widetilde{u},\widetilde{\overline{u}},v) be the resulting relation.

If f1(u,u¯,u¯~,v)=f2(u,u~,u¯~,v)f_{1}(u,\overline{u},\widetilde{\overline{u}},v)=f_{2}(u,\widetilde{u},\widetilde{\overline{u}},v) is automatically equal, then the CAC-type BT (A.1b) and (A.1c) is said to be trivial. If only Equation (A.5a) is obtained from f1(u,u¯,u¯~,v)=f2(u,u~,u¯~,v)f_{1}(u,\overline{u},\widetilde{\overline{u}},v)=f_{2}(u,\widetilde{u},\widetilde{\overline{u}},v), then the CAC-type BT (A.1b) and (A.1c) is said to be strong. If f1(u,u¯,u¯~,v)=f2(u,u~,u¯~,v)f_{1}(u,\overline{u},\widetilde{\overline{u}},v)=f_{2}(u,\widetilde{u},\widetilde{\overline{u}},v) yields Equation (A.5a) plus other PΔ\DeltaEs for {u,u¯,u~,u¯~}\{u,\overline{u},\widetilde{u},\widetilde{\overline{u}}\}, then the CAC-type BT (A.1b) and (A.1c) is said to be weak. A trivial BT gives a fake Lax pair, and a weak BT gives a weak Lax pair. It is claimed in [8] that if the CAC-type BT is strong, then the integrability of (A.1a) should be guaranteed.

Note that the CAC-type BTs given in Theorems 2.1, 2.5 and 2.7 are all strong.

Refer to caption
Figure A.1. A cube for the system (A.1). The uu- and vv-variables are assigned at the bottom and top vertices of the cube, respectively. Also, the equations (A.5) are assigned to the faces of the cube.

A.2. The CABC property

Let us consider the following system of PΔ\DeltaEs:

A(u,u¯,u~,u¯~)=0,\displaystyle A(u,\overline{u},\widetilde{u},\widetilde{\overline{u}})=0, (A.6a)
S(u,u¯,v~,v¯~)=0,\displaystyle S(u,\overline{u},\widetilde{v},\widetilde{\overline{v}})=0, (A.6b)
B(u,v,v~)=0,\displaystyle B(u,v,\widetilde{v})=0, (A.6c)
C(u,u¯,v,v¯)=0,\displaystyle C(u,\overline{u},v,\overline{v})=0, (A.6d)

where u=ul,mu=u_{l,m} and v=vl,mv=v_{l,m}. Here, the functions AA, SS and CC are irreducible multilinear polynomials in four variables. The function B=B(x,y,z)B=B(x,y,z) is a polynomial in three variables satisfying the following:

  1. 1)

    degxB1\deg_{x}B\geq 1, degyB=degzB=1\deg_{y}B=\deg_{z}B=1;

  2. 2)

    Let y=f(x,z)y=f(x,z) be the solution of B=0B=0. Then, f(x,z)f(x,z) is a rational function that depends on xx and zz, that is, the following hold:

    xf(x,z)0,zf(x,z)0.\dfrac{\partial}{\partial x}f(x,z)\neq 0,\quad\dfrac{\partial}{\partial z}f(x,z)\neq 0. (A.7)

Assign the uu- and vv-variables on the sublattice (A.2) by the correspondence (A.3). Then, by considering the cube consisting of the eight vertices (A.4) (see Figure A.2), the system of equations around the cube obtained from the system (A.6) is the following:

𝒜=A(u,u¯,u~,u¯~)=0,\displaystyle{\mathcal{A}}=A(u,\overline{u},\widetilde{u},\widetilde{\overline{u}})=0, (A.8a)
𝒮=S(u,u¯,v~,v¯~)=0,\displaystyle{\mathcal{S}}=S(u,\overline{u},\widetilde{v},\widetilde{\overline{v}})=0, (A.8b)
=B(u,v,v~)=0,\displaystyle{\mathcal{B}}=B(u,v,\widetilde{v})=0, (A.8c)
¯=B(u¯,v¯,v¯~)=0,\displaystyle\overline{{\mathcal{B}}}=B(\overline{u},\overline{v},\widetilde{\overline{v}})=0, (A.8d)
𝒞=C(u,u¯,v,v¯)=0,\displaystyle{\mathcal{C}}=C(u,\overline{u},v,\overline{v})=0, (A.8e)
𝒞~=C(u~,u¯~,v~,v¯~)=0.\displaystyle\widetilde{{\mathcal{C}}}=C(\widetilde{u},\widetilde{\overline{u}},\widetilde{v},\widetilde{\overline{v}})=0. (A.8f)
Refer to caption
Figure A.2. A cube for the system (A.6). The uu- and vv-variables are assigned at the bottom and top vertices of the cube, respectively. Also, each of the equations (A.8) is assigned to a face of the cube. Note that the face for Equation (A.8b) corresponds to the cutting plane bisecting the cube diagonally, and the faces for the equations (A.8c) and (A.8d) correspond to the upper triangles.

The CABC and tetrahedron properties are defined as follows.

  1. (1)(\mathrm{1})

    There are following three ways to calculate v¯~\widetilde{\overline{v}} by using all the equations in the system (A.8) with {u,u¯,u~,v}\{u,\overline{u},\widetilde{u},v\} as initial values.

    1. (a)

      Express v¯~\widetilde{\overline{v}} as a rational function in terms of {u,u¯,v~}\{u,\overline{u},\widetilde{v}\} by using Equation (A.8b). Then, eliminate v~\widetilde{v} from it by using Equation (A.8c).

    2. (b)

      Express v¯~\widetilde{\overline{v}} as a rational function in terms of {u¯,v¯}\{\overline{u},\overline{v}\} by using Equation (A.8d). Then, eliminate v¯\overline{v} from it by using Equation (A.8e).

    3. (c)

      Express v¯~\widetilde{\overline{v}} as a rational function in terms of {u~,u¯~,v~}\{\widetilde{u},\widetilde{\overline{u}},\widetilde{v}\} by using Equation (A.8f). Then, eliminate u¯~\widetilde{\overline{u}} by using Equation (A.8a) and v~\widetilde{v} by using Equation (A.8c) from it.

    When v¯~\widetilde{\overline{v}} is uniquely determined as a rational function with the initial values {u,u¯,u~,v}\{u,\overline{u},\widetilde{u},v\}, then the system (A.6) is said to have the CABC property or said to be a CABC-system; the PΔ\DeltaE (A.6a) is said to have the CABC property.

  2. (2)(\mathrm{2})

    Let the system (A.6) be a CABC-system. Then, the system (A.6) is said to have the tetrahedron property, if the equations

    K1(u,u¯,v,v¯~)=0,K2(u,u¯,v¯,v~)=0,K_{1}(u,\overline{u},v,\widetilde{\overline{v}})=0,\quad K_{2}(u,\overline{u},\overline{v},\widetilde{v})=0, (A.9)

    where K1K_{1} and K2K_{2} are irreducible multilinear polynomials in four variables, can be obtained from the system (A.6).

  3. (3)(\mathrm{3})

    If there exists a PΔ\DeltaE given only by {v,v¯,v~,v¯~}\{v,\overline{v},\widetilde{v},\widetilde{\overline{v}}\}, then the tuple of equations (A.6b), (A.6c) and (A.6d) is referred to as a CABC-type BT from the PΔ\DeltaE (A.6a) to the PΔ\DeltaE given by {v,v¯,v~,v¯~}\{v,\overline{v},\widetilde{v},\widetilde{\overline{v}}\}.

Appendix B Lax pairs of the dKdV equation and the lsG equation

In this appendix, we show Lax pairs of the dKdV equation (1.1) and the lsG equation (1.2) obtained in this study.

B.1. Lax pairs of the PΔ\DeltaE (1.1)

A Lax pair of the PΔ\DeltaE (1.1) is given by

Φl+1,m=Ll,mΦl,m,Φl,m+1=Ml,mΦl,m,\Phi_{l+1,m}=L_{l,m}\Phi_{l,m},\quad\Phi_{l,m+1}=M_{l,m}\Phi_{l,m}, (B.1)

where the matrices Ll,mL_{l,m} and Ml,mM_{l,m} are given by one of the following pairs:

Ll,m=(0γαl+11βmαluu¯),Ml,m=(uγαl1βmαlu),\displaystyle L_{l,m}=\begin{pmatrix}0&\gamma-\alpha_{l+1}\\ 1&\dfrac{\beta_{m}-\alpha_{l}}{u}-\overline{u}\end{pmatrix},\quad M_{l,m}=\begin{pmatrix}u&\gamma-\alpha_{l}\\ 1&\dfrac{\beta_{m}-\alpha_{l}}{u}\end{pmatrix}, (B.2a)
Ll,m=(qm2uu¯plqmupl+1qmupl+1plu),Ml,m=(qm+1qmuplqm+1uplqmuu~pl2u),\displaystyle L_{l,m}=\begin{pmatrix}\dfrac{{q_{m}}^{2}}{u}-\overline{u}&-\dfrac{p_{l}q_{m}}{u}\\[5.0pt] \dfrac{p_{l+1}q_{m}}{u}&-\dfrac{p_{l+1}p_{l}}{u}\end{pmatrix},\quad M_{l,m}=\begin{pmatrix}\dfrac{q_{m+1}q_{m}}{u}&-\dfrac{p_{l}q_{m+1}}{u}\\[5.0pt] \dfrac{p_{l}q_{m}}{u}&\widetilde{u}-\dfrac{{p_{l}}^{2}}{u}\end{pmatrix}, (B.2b)
Ll,m=(pl+1plupl+1qmuplqmuu¯qm2u),Ml,m=(pl2uu~plqmuplqm+1uqm+1qmu).\displaystyle L_{l,m}=\begin{pmatrix}\dfrac{p_{l+1}p_{l}}{u}&-\dfrac{p_{l+1}q_{m}}{u}\\[5.0pt] \dfrac{p_{l}q_{m}}{u}&\overline{u}-\dfrac{{q_{m}}^{2}}{u}\end{pmatrix},\quad M_{l,m}=\begin{pmatrix}\dfrac{{p_{l}}^{2}}{u}-\widetilde{u}&-\dfrac{p_{l}q_{m}}{u}\\[5.0pt] \dfrac{p_{l}q_{m+1}}{u}&-\dfrac{q_{m+1}q_{m}}{u}\end{pmatrix}. (B.2c)

Here, plp_{l} and qmq_{m} are given by Equation (1.8) and γ\gamma\in\mathbb{C} is a spectral parameter. Indeed, for each pair in (B.2), the compatibility condition

Ll,m+1Ml,m=Ml+1,mLl,mL_{l,m+1}M_{l,m}=M_{l+1,m}L_{l,m} (B.3)

gives the PΔ\DeltaE (1.1).

B.2. A Lax pair of the PΔ\DeltaE (1.2)

The linear system

Ψl+1,m=l,mΨl,m,Ψl,m+1=l,mΨl,m\Psi_{l+1,m}={\mathcal{L}}_{l,m}\Psi_{l,m},\quad\Psi_{l,m+1}={\mathcal{M}}_{l,m}\Psi_{l,m} (B.4a)
where
l,m=(pl+1qmx¯+Qm2x¯qmplxPlQmx¯plxqmPl+1Qmx¯qmplxPlPl+1x¯plxqm),\displaystyle{\mathcal{L}}_{l,m}=\begin{pmatrix}p_{l+1}-q_{m}\overline{x}+\dfrac{{Q_{m}}^{2}\overline{x}}{q_{m}-p_{l}x}&\dfrac{P_{l}Q_{m}\overline{x}}{p_{l}x-q_{m}}\\[10.00002pt] \dfrac{P_{l+1}Q_{m}\overline{x}}{q_{m}-p_{l}x}&\dfrac{P_{l}P_{l+1}\overline{x}}{p_{l}x-q_{m}}\end{pmatrix}, (B.4b)
l,m=(QmQm+1qmplxPlQm+1plxqmPlQmqmplxPl2plxqm+qm+1x~plx),\displaystyle{\mathcal{M}}_{l,m}=\begin{pmatrix}\dfrac{Q_{m}Q_{m+1}}{q_{m}-p_{l}x}&\dfrac{P_{l}Q_{m+1}}{p_{l}x-q_{m}}\\[10.00002pt] \dfrac{P_{l}Q_{m}}{q_{m}-p_{l}x}&\dfrac{{P_{l}}^{2}}{p_{l}x-q_{m}}+\dfrac{q_{m+1}\widetilde{x}-p_{l}}{x}\end{pmatrix}, (B.4c)

is a Lax pair of the PΔ\DeltaE (1.2). Here, PlP_{l} and QmQ_{m} are given by Equation (2.20) and λ\lambda\in\mathbb{C} is a spectral parameter. Indeed, the compatibility condition

l,m+1l,m=l+1,ml,m{\mathcal{L}}_{l,m+1}{\mathcal{M}}_{l,m}={\mathcal{M}}_{l+1,m}{\mathcal{L}}_{l,m} (B.5)

gives the PΔ\DeltaE (1.2).

Appendix C Proof of Theorem 3.3

Define ω=ωl,m,n\omega=\omega_{l,m,n} by

ω=τ^τ,\omega=\dfrac{\,\widehat{\tau}\,}{\tau}, (C.1)

where τ=τl,m,n\tau=\tau_{l,m,n} is a function satisfying the bilinear equations (3.3). Then, using Equation (3.4), we can express the variables u=ul,m,nu=u_{l,m,n}, v=vl,m,nv=v_{l,m,n} and x=xl,m,nx=x_{l,m,n} in the ω\omega-variable as

u=ττ¯~τ¯τ~=pl,nτ¯^τ¯qm,nτ~^τ~τ^τ=pl,nω¯qm,nω~ω,\displaystyle u=-\dfrac{\,\tau\,\widetilde{\overline{\tau}}\,}{\overline{\tau}\,\widetilde{\tau}}=\cfrac{p_{l,n}\dfrac{\,\widehat{\overline{\tau}}\,}{\overline{\tau}}-q_{m,n}\dfrac{\,\widehat{\widetilde{\tau}}\,}{\widetilde{\tau}}}{\dfrac{\,\widehat{\tau}\,}{\tau}}=\dfrac{p_{l,n}\overline{\omega}-q_{m,n}\widetilde{\omega}}{\omega}, (C.2a)
v=pl,nττ¯^τ¯τ^=pl,nω¯ω,\displaystyle v=p_{l,n}\,\dfrac{\,\tau\,\widehat{\overline{\tau}}\,}{\overline{\tau}\,\widehat{\tau}}=p_{l,n}\dfrac{\,\overline{\omega}\,}{\omega}, (C.2b)
x=τ¯^τ~τ¯τ~^=ω¯ω~.\displaystyle x=\dfrac{\,\widehat{\overline{\tau}}\,\widetilde{\tau}\,}{\overline{\tau}\,\widehat{\widetilde{\tau}}}=\dfrac{\,\overline{\omega}\,}{\widetilde{\omega}}. (C.2c)

Note that we use Equation (3.3b) for the deformation in Equation (C.2a). Moreover, the following lemma holds.

Lemma C.1.

The following relations hold:

ω¯~ω=qm,nω¯pl,nω~qm,nω~pl,nω¯,\displaystyle\dfrac{\,\widetilde{\overline{\omega}}\,}{\omega}=\dfrac{q_{m,n}\overline{\omega}-p_{l,n}\widetilde{\omega}}{q_{m,n}\widetilde{\omega}-p_{l,n}\overline{\omega}}, (C.3a)
ωω^ω¯ω~=pl,n+1ω¯^qm,n+1ω~^pl,nω~qm,nω¯.\displaystyle\dfrac{\omega\,\widehat{\omega}}{\overline{\omega}\,\widetilde{\omega}}=\dfrac{p_{l,n+1}\widehat{\overline{\omega}}-q_{m,n+1}\widehat{\widetilde{\omega}}}{p_{l,n}\widetilde{\omega}-q_{m,n}\overline{\omega}}. (C.3b)
Proof.

Using Equation (3.3a), we obtain

ω¯~=τ¯~^τ¯~=qm,nτ¯^τ~pl,mτ¯τ~^τ¯~τ,\widetilde{\overline{\omega}}=\dfrac{\,\widehat{\widetilde{\overline{\tau}}}\,}{\widetilde{\overline{\tau}}}=\dfrac{q_{m,n}\widehat{\overline{\tau}}~{}\widetilde{\tau}-p_{l,m}\overline{\tau}~{}\widehat{\widetilde{\tau}}}{\widetilde{\overline{\tau}}~{}\tau}, (C.4)

and using Equation (3.3b), we obtain

ω(qm,nω¯pl,mω~)qm,nω~pl,mω¯=τ^(qm,nτ¯^τ~pl,mτ¯τ~^)τ(qm,nτ¯τ~^pl,mτ~τ¯^)=qm,nτ¯^τ~pl,mτ¯τ~^τ¯~τ.\dfrac{\omega(q_{m,n}\overline{\omega}-p_{l,m}\widetilde{\omega})}{q_{m,n}\widetilde{\omega}-p_{l,m}\overline{\omega}}=\dfrac{\widehat{\tau}(q_{m,n}\widehat{\overline{\tau}}~{}\widetilde{\tau}-p_{l,m}\overline{\tau}~{}\widehat{\widetilde{\tau}})}{\tau(q_{m,n}\overline{\tau}~{}\widehat{\widetilde{\tau}}-p_{l,m}\widetilde{\tau}~{}\widehat{\overline{\tau}})}=\dfrac{q_{m,n}\widehat{\overline{\tau}}~{}\widetilde{\tau}-p_{l,m}\overline{\tau}~{}\widehat{\widetilde{\tau}}}{\widetilde{\overline{\tau}}~{}\tau}. (C.5)

Equation (C.3a) follows from these two equations. Moreover, using the equations (3.3b)n→n+1 and (3.3a) in turn, we obtain

pl,n+1ω¯^qm,n+1ω~^ω^(pl,nω~qm,nω¯)\displaystyle\dfrac{p_{l,n+1}\widehat{\overline{\omega}}-q_{m,n+1}\widehat{\widetilde{\omega}}}{\widehat{\omega}(p_{l,n}\widetilde{\omega}-q_{m,n}\overline{\omega})} =τ¯τ~τ^(pl,n+1τ~^τ¯^^qm,n+1τ¯^τ~^^)τ¯^τ~^τ^^(pl,nτ¯τ~^qm,nτ~τ¯^)=τ¯τ~τ^τ¯~^τ¯^τ~^(pl,nτ¯τ~^qm,nτ~τ¯^)\displaystyle=\dfrac{\overline{\tau}~{}\widetilde{\tau}~{}\widehat{\tau}(p_{l,n+1}\widehat{\widetilde{\tau}}~{}\widehat{\widehat{\overline{\tau}}}-q_{m,n+1}\widehat{\overline{\tau}}~{}\widehat{\widehat{\widetilde{\tau}}})}{\widehat{\overline{\tau}}~{}\widehat{\widetilde{\tau}}~{}\widehat{\widehat{\tau}}(p_{l,n}\overline{\tau}~{}\widehat{\widetilde{\tau}}-q_{m,n}\widetilde{\tau}~{}\widehat{\overline{\tau}})}=-\dfrac{\overline{\tau}~{}\widetilde{\tau}~{}\widehat{\tau}~{}\widehat{\widetilde{\overline{\tau}}}}{\widehat{\overline{\tau}}~{}\widehat{\widetilde{\tau}}(p_{l,n}\overline{\tau}~{}\widehat{\widetilde{\tau}}-q_{m,n}\widetilde{\tau}~{}\widehat{\overline{\tau}})}
=τ¯τ~τ^ττ¯^τ~^=ωω¯ω~,\displaystyle=\dfrac{~{}\overline{\tau}~{}\widetilde{\tau}~{}\widehat{\tau}~{}}{\tau~{}\widehat{\overline{\tau}}~{}\widehat{\widetilde{\tau}}}=\dfrac{\omega}{\,\overline{\omega}\,\widetilde{\omega}\,}, (C.6)

which is equivalent to Equation (C.3b). ∎

Remark C.2.

Equation (C.3a) is known as the lattice modified KdV equation [17, 19, 1].

Consider the system (3.1) expressed in the ω\omega-variable by using the equations (C.2). It is sufficient to show that the system (3.1) expressed in the ω\omega-variable can be solved by the equations (C.3). For simplicity, we omit the description of the specific equations written in the ω\omega-variable and itemize the method below by using only words.

  • It is obvious that the equations (3.1c), (3.1d), (3.1e) and (3.1f) written in terms of the ω\omega-variable are equivalent to Equation (C.3a).

  • Firstly eliminating ω¯¯~\widetilde{\overline{\overline{\omega}}} by using Equation (C.3a)l→l+1, secondly eliminating ω¯~~\widetilde{\widetilde{\overline{\omega}}} by using Equation (C.3a)m→m+1, and finally eliminating ω¯~\widetilde{\overline{\omega}} by using Equation (C.3a), we can show that the equations (3.1a) and (3.1b) written in terms of the ω\omega-variable hold.

  • Eliminating ω¯~\widetilde{\overline{\omega}} by using Equation (C.3a) and then eliminating ω¯^\widehat{\overline{\omega}} by using (C.3b), we can show that Equation (3.1g) written in terms of the ω\omega-variable holds.

  • Firstly eliminating ω¯¯^\widehat{\overline{\overline{\omega}}} by using Equation (C.3b)l→l+1, secondly eliminating ω¯~^\widehat{\widetilde{\overline{\omega}}} by using Equation (C.3b)m→m+1, thirdly eliminating ω¯^\widehat{\overline{\omega}} by using Equation (C.3b), and finally eliminating ω¯~\widetilde{\overline{\omega}} by using Equation (C.3a), we can show that Equation (3.1h) written in terms of the ω\omega-variable holds.

Therefore, we have completed the proof of Theorem 3.3.

References

  • [1] V. E. Adler, A. I. Bobenko, and Y. B. Suris. Classification of integrable equations on quad-graphs. The consistency approach. Comm. Math. Phys., 233(3):513–543, 2003.
  • [2] V. E. Adler, A. I. Bobenko, and Y. B. Suris. Discrete nonlinear hyperbolic equations: classification of integrable cases. Funktsional. Anal. i Prilozhen., 43(1):3–21, 2009.
  • [3] A. I. Bobenko and Y. B. Suris. Integrable systems on quad-graphs. Int. Math. Res. Not. IMRN, (11):573–611, 2002.
  • [4] R. Boll. Classification of 3D consistent quad-equations. J. Nonlinear Math. Phys., 18(3):337–365, 2011.
  • [5] R. Boll. Corrigendum: Classification of 3D consistent quad-equations. J. Nonlinear Math. Phys., 19(4):1292001, 3, 2012.
  • [6] S. Butler and M. Hay. Two definitions of fake Lax pairs. In AIP Conference Proceedings, volume 1648, page 180006. AIP Publishing LLC, 2015.
  • [7] H. W. Capel, F. W. Nijhoff, and V. G. Papageorgiou. Complete integrability of Lagrangian mappings and lattices of KdV type. Phys. Lett. A, 155(6-7):377–387, 1991.
  • [8] J. Hietarinta. Search for CAC-integrable homogeneous quadratic triplets of quad equations and their classification by BT and Lax. J. Nonlinear Math. Phys., 26(3):358–389, 2019.
  • [9] J. Hietarinta, N. Joshi, and F. W. Nijhoff. Discrete systems and integrability. Cambridge Texts in Applied Mathematics. Cambridge University Press, Cambridge, 2016.
  • [10] J. Hietarinta and C. Viallet. Weak Lax pairs for lattice equations. Nonlinearity, 25(7):1955–1966, 2012.
  • [11] N. Joshi and N. Nakazono. On the three-dimensional consistency of Hirota’s discrete Korteweg-de Vries equation. Stud. Appl. Math., 147(4):1409–1424, 2021.
  • [12] K. Kajiwara and Y. Ohta. Bilinearization and Casorati determinant solution to the non-autonomous discrete KdV equation. Journal of the Physical Society of Japan, 77(5):054004, 2008.
  • [13] P. Kassotakis and M. Nieszporski. Difference systems in bond and face variables and non-potential versions of discrete integrable systems. J. Phys. A, Math. Theor., 51(38):21, 2018. Id/No 385203.
  • [14] D. J. Korteweg and G. De Vries. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Phil. Mag. (5), 39:422–443, 1895.
  • [15] N. Nakazono. Properties of the non-autonomous lattice sine-Gordon equation: consistency around a broken cube property. SIGMA. Symmetry, Integrability and Geometry: Methods and Applications, 18:032, 2022.
  • [16] F. W. Nijhoff. Lax pair for the Adler (lattice Krichever-Novikov) system. Phys. Lett. A, 297(1-2):49–58, 2002.
  • [17] F. W. Nijhoff and H. W. Capel. The discrete Korteweg-de Vries equation. Acta Appl. Math., 39(1-3):133–158, 1995. KdV ’95 (Amsterdam, 1995).
  • [18] F. W. Nijhoff, H. W. Capel, G. L. Wiersma, and G. R. W. Quispel. Bäcklund transformations and three-dimensional lattice equations. Phys. Lett. A, 105(6):267–272, 1984.
  • [19] F. W. Nijhoff, G. R. W. Quispel, and H. W. Capel. Direct linearization of nonlinear difference-difference equations. Phys. Lett. A, 97(4):125–128, 1983.
  • [20] F. W. Nijhoff and A. J. Walker. The discrete and continuous Painlevé VI hierarchy and the Garnier systems. Glasg. Math. J., 43A:109–123, 2001. Integrable systems: linear and nonlinear dynamics (Islay, 1999).
  • [21] J. J. C. Nimmo and W. K. Schief. An integrable discretization of a (2+1)(2+1)-dimensional sine-Gordon equation. Stud. Appl. Math., 100(3):295–309, 1998.
  • [22] G. R. W. Quispel, F. W. Nijhoff, H. W. Capel, and J. van der Linden. Linear integral equations and nonlinear difference-difference equations. Phys. A, 125(2-3):344–380, 1984.
  • [23] S. Tremblay, B. Grammaticos, and A. Ramani. Integrable lattice equations and their growth properties. Phys. Lett., A, 278(6):319–324, 2001.
  • [24] A. Walker. Similarity reductions and integrable lattice equations. Ph.D. Thesis, University of Leeds, 2001.