搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

离散可积系统: 多维相容性

张大军

离散可积系统: 多维相容性

张大军
PDF
HTML
导出引用
导出核心图
  • 对比已有完善而系统理论的微分方程领域, 差分方程理论尚处于发展之中. 近年来离散可积理论的进展, 带来了差分方程理论的革命. 多维相容性是伴随离散可积系统研究出现的新的概念, 作为对离散可积性的一种理解, 提供了构造离散可积系统的Bäcklund变换、Lax对和精确解的工具. 本文旨在综述多维相容性的概念及其在离散可积系统研究中的应用.
      通信作者: 张大军, djzhang@staff.shu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11875040, 11631007)
    [1]

    Zabusky N J, Kruskal M D 1965 Phys. Rev. Lett. 15 240

    [2]

    Gardner C S, Greene J M, Kruskal M D, Miura R M 1967 Phys. Rev. Lett. 19 1095

    [3]

    Lax P D 1968 Commun. Pure Appl. Math. 21 467

    [4]

    Deift P 2019 arXiv: 1902.10267

    [5]

    Case K M, Kac M 1973 J. Math. Phys. 14 594

    [6]

    Ablowitz M J, Ladik J F 1975 J. Math. Phys. 16 598

    [7]

    Ablowitz M J, Ladik J F 1976 J. Math. Phys. 17 1011

    [8]

    Ablowitz M J, Ladik J F 1976 Stud. Appl. Math. 55 213

    [9]

    Hirota R 1977 J. Phys. Soc. Jpn. 43 1424

    [10]

    Hirota R 1977 J. Phys. Soc. Jpn. 43 2074

    [11]

    Hirota R 1977 J. Phys. Soc. Jpn. 43 2079

    [12]

    Date E, Jimbo M, Miwa T 1982 J. Phys. Soc. Jpn. 51 4116

    [13]

    Date E, Jimbo M, Miwa T 1982 J. Phys. Soc. Jpn. 51 4125

    [14]

    Date E, Jimbo M, Miwa T 1983 J. Phys. Soc. Jpn. 52 388

    [15]

    Date E, Jimbo M, Miwa T 1983 J. Phys. Soc. Jpn. 52 761

    [16]

    Date E, Jimbo M, Miwa T 1983 J. Phys. Soc. Jpn. 52 766

    [17]

    Ueno K, Takasaki K 1984 In: Okamoto K (ed) Group Representations and Systems of Differential Equations Advanced Studies in Pure Mathematics (Vol. 4) (Tokyo: Kinokuniya) pp1−95

    [18]

    Nijhoff F W, Quispel G R W, Capel H W 1983 Phys. Lett. A 97 125

    [19]

    Nijhoff F W, Quispel G R W, Capel H W 1983 Phys. Lett. A 98 83

    [20]

    Nijhoff F W, Capel H W, Wiersma G L, Quispel G R W 1984 Phys. Lett. A 105 267

    [21]

    Nijhoff F W 1985 Lett. Math. Phys. 9 235

    [22]

    Nijhoff F W, Capel H W, Wiersma G L 1985 In: Martini R (ed) Geometric Aspects of the Einstein Equations and Integrable Systems (Scheveningen 1984) Lecture Notes in Phys (Vol. 239) (Berlin: Springer) pp263−302

    [23]

    Quispel G R W, Nijhoff F W, Capel H W, van ver Linden J 1984 Physica A 125 344

    [24]

    Fokas A S, Ablowitz M 1981 Phys. Rev. Lett. 47 1096

    [25]

    Levi D, Benguria R 1980 Proc. Natl. Acad. Sci. U.S.A. 77 5025

    [26]

    Takahashi D, Satsuma J 1990 J. Phys. Soc. Jpn. 59 3514

    [27]

    Tokihiro T, Takahashi D, Matsukidaira J, Satsuma J 1996 Phys. Rev. Lett. 76 3247

    [28]

    Grammaticos B, Ramani A, Papageorgiou V G 1991 Phys. Rev. Lett. 67 1825

    [29]

    Ramani A, Grammaticos B, Hietarinta J 1991 Phys. Rev. Lett. 67 1829

    [30]

    Hietarinta J, Viallet C 1998 Phys. Rev. Lett. 81 325

    [31]

    Bellon M P, Viallet C 1999 Commun. Math. Phys. 204 425

    [32]

    Sakai H 2001 Commun. Math. Phys. 220 165

    [33]

    Bobenko A I, Suris Yu B 2002 Int. Math. Res. Not. 2002 573

    [34]

    Adler V E, Bobenko A I, Suris Yu B 2003 Commun. Math. Phys. 233 513

    [35]

    Nijhoff F W, Walker A J 2001 Glasg. Math. J. 43A 109

    [36]

    Nijhoff F W 2002 Phys. Lett. A 297 49

    [37]

    Nijhoff F W, Atkinson J, Hietarinta J 2009 J. Phys. A: Math. Theor. 42 404005

    [38]

    Hietarinta J, Zhang D J 2008 J. Phys. A: Math. Theor. 42 404006

    [39]

    Atkinson J, Nijhoff F W 2010 Commun. Math. Phys. 299 283

    [40]

    Nijhoff F W, Atkinson J 2010 Int. Math. Res. Not. 2010 3837

    [41]

    Butler S, Joshi N 2010 Inverse Prob. 26 115012

    [42]

    Butler S 2012 Nonlinearity 25 1613

    [43]

    Cao C W, Xu X X 2012 J. Phys. A: Math. Theor. 45 055213

    [44]

    Cao C W, Zhang G Y 2012 J. Phys. A: Math. Theor. 45 095203

    [45]

    Zhang D J, Zhao S L 2013 Stud. Appl. Math. 131 72

    [46]

    Bobenko A I, Its A 2016 Duke Math. J. 165 2607

    [47]

    Hietarinta J, Joshi N, Nijhoff F W 2016 (Cambridge: Camb. Univ. Press)

    [48]

    Zhang D J, Chen S T 2010 Stud. Appl. Math. 125 393

    [49]

    Zhang D J, Chen S T 2010 Stud. Appl. Math. 125 419

    [50]

    Fu W, Qiao Z J, Sun J W, Zhang D J 2013 arXiv: 1307.3671

    [51]

    Fu W, Qiao Z J, Sun J W, Zhang D J 2015 J. Nonlinear. Math. Phys. 22 321

    [52]

    Wahlquist H D, Estabrook F B 1973 Phys. Rev. Lett. 31 1386

    [53]

    Lamb JR G L 1971 Rev. Mod. Phys. 43 99

    [54]

    Chen H H 1974 Phys. Rev. Lett. 33 925

    [55]

    Orfanidis S J 1978 Phys. Rev. D 18 3828

    [56]

    Bianchi L 1892 Rend. Lincei 5 2

    [57]

    Bianchi L 1894 Lezioni di Geometria Differenziale (3rd Ed.) (Pisa: Enrico Spoerri)

    [58]

    Konopelchenko B G 1982 Phys. Lett. A 87 445

    [59]

    Levi D 1981 J. Phys. A: Math. Gen. 14 1083

    [60]

    Adler V E, Yamilov R I 1994 J. Phys. A: Math. Gen. 27 477

    [61]

    Merola I, Ragnisco O, Tu G Z 1994 Inverse Prob. 10 1315

    [62]

    Zhang H W, Tu G Z, Oevel W, Fuchssteiner B 1991 J. Math. Phys. 32 1908

    [63]

    Chen K, Deng X, Zhang D J 2017 J. Nonlinear. Math. Phys. 24(Suppl.1) 18

    [64]

    Cao C W, Zhang G Y 2012 Chin. Phys. Lett. 29 050202

    [65]

    Walker A J 2001 Ph.D. Thesis (Leeds: University of Leeds)

    [66]

    Adler V E 1998 Int. Math. Res. Not. 1998 1

    [67]

    Hietarinta J 2005 J. Nonlinear. Math. Phys. 12 223

    [68]

    Nijhoff F W, Papageorgiou V G, Capel H W, Quispel G R W 1992 Inverse Prob. 8 597

    [69]

    Nijhoff F W 1997 In: Fokas A S, Gel’fand I M (eds) Algebraic Aspects of Integrable Systems: In memory of Irene Dorfman (Boston: Birkhauser) pp237−260

    [70]

    Nijhoff F W 1999 In: Bobenko A I, Seiler R (eds) Discrete Integrable Geometry and Physics (Oxford: Clarendon Press) pp209−234

    [71]

    Hietarinta J 2011 J. Phys. A: Math. Theor. 44 165204

    [72]

    Hietarinta J, Zhang D J 2008 preprint

    [73]

    Atkinson J 2008 J. Phys. A: Math. Theor. 41 135202

    [74]

    Adler V E, Bobenko A I, Suris Yu B 2009 Funct. Anal. Appl. 43 3

    [75]

    Boll R 2011 J. Nonlinear. Math. Phys. 18 337

    [76]

    Boll R 2012 Ph.D Dissertation (Berlin: Technischen Universität Berlin)

    [77]

    Adler V E, Bobenko A I, Suris Yu B 2012 Int. Math. Res. Not. 2012 1822

    [78]

    Miwa T 1982 Proc. Jpn. Acad. 58A 9

    [79]

    Konopelchenko B G, Schief W K 2002 Stud. Appl. Math. 109 89

    [80]

    Hirota R 1981 J. Phys. Soc. Jpn. 50 3785

    [81]

    Bianchi L 1885 Ann. Matem. 13 177

    [82]

    Atkinson J, Nieszporski M 2014 Int. Math. Res. Not. 2014 4215

    [83]

    Zhang D D, Zhang D J 2018 J. Nonlinear. Math. Phys. 25 34

    [84]

    Zhang D J, Cheng J W, Sun Y Y 2013 J. Phys. A: Math. Theor. 46 265202

    [85]

    Xenitidis P 2011 J. Phys. A: Math. Theor. 44 435201

    [86]

    Bridgman T, Hereman W, Quispel G R W, van der Kamp P H 2013 Found. Comput. Math. 13 517

    [87]

    Hietarinta J, Zhang D J 2010 J. Math. Phys. 51 033505

    [88]

    Hietarinta J, Zhang D J 2011 SIGMA 7 061

    [89]

    Atkinson J, Hietarinta J, Nijhoff F W 2007 J. Phys. A: Math. Theor. 40 F1

    [90]

    Konopelchenko B G, Schief W K 2002 J. Phys. A: Math. Gen. 35 6125

  • 图 1  $[x_0, x]$上的数值离散

    Fig. 1.  Numerical discretisation on $[x_0, x]$

    图 2  Bäcklund变换解的交换性质

    Fig. 2.  Permutation property of Bäcklund transformation

    图 3  变换与平面网格

    Fig. 3.  Map and lattice

    图 4  相容立方体

    Fig. 4.  Consistent cube

    图 5  离散sine-Gordon和势mKdV方程的相容立方体

    Fig. 5.  The consistent cube for the discrete sine-Gordon equation and potential mKdV equation

    图 6  定义3维方程的6面体以及8面体

    Fig. 6.  Cube and octahedron for 3D equations

    图 7  围绕超立方体的4D相容性

    Fig. 7.  4D consistency around the hyper cube

    图 8  相容立方体

    Fig. 8.  The consistent cube

  • [1]

    Zabusky N J, Kruskal M D 1965 Phys. Rev. Lett. 15 240

    [2]

    Gardner C S, Greene J M, Kruskal M D, Miura R M 1967 Phys. Rev. Lett. 19 1095

    [3]

    Lax P D 1968 Commun. Pure Appl. Math. 21 467

    [4]

    Deift P 2019 arXiv: 1902.10267

    [5]

    Case K M, Kac M 1973 J. Math. Phys. 14 594

    [6]

    Ablowitz M J, Ladik J F 1975 J. Math. Phys. 16 598

    [7]

    Ablowitz M J, Ladik J F 1976 J. Math. Phys. 17 1011

    [8]

    Ablowitz M J, Ladik J F 1976 Stud. Appl. Math. 55 213

    [9]

    Hirota R 1977 J. Phys. Soc. Jpn. 43 1424

    [10]

    Hirota R 1977 J. Phys. Soc. Jpn. 43 2074

    [11]

    Hirota R 1977 J. Phys. Soc. Jpn. 43 2079

    [12]

    Date E, Jimbo M, Miwa T 1982 J. Phys. Soc. Jpn. 51 4116

    [13]

    Date E, Jimbo M, Miwa T 1982 J. Phys. Soc. Jpn. 51 4125

    [14]

    Date E, Jimbo M, Miwa T 1983 J. Phys. Soc. Jpn. 52 388

    [15]

    Date E, Jimbo M, Miwa T 1983 J. Phys. Soc. Jpn. 52 761

    [16]

    Date E, Jimbo M, Miwa T 1983 J. Phys. Soc. Jpn. 52 766

    [17]

    Ueno K, Takasaki K 1984 In: Okamoto K (ed) Group Representations and Systems of Differential Equations Advanced Studies in Pure Mathematics (Vol. 4) (Tokyo: Kinokuniya) pp1−95

    [18]

    Nijhoff F W, Quispel G R W, Capel H W 1983 Phys. Lett. A 97 125

    [19]

    Nijhoff F W, Quispel G R W, Capel H W 1983 Phys. Lett. A 98 83

    [20]

    Nijhoff F W, Capel H W, Wiersma G L, Quispel G R W 1984 Phys. Lett. A 105 267

    [21]

    Nijhoff F W 1985 Lett. Math. Phys. 9 235

    [22]

    Nijhoff F W, Capel H W, Wiersma G L 1985 In: Martini R (ed) Geometric Aspects of the Einstein Equations and Integrable Systems (Scheveningen 1984) Lecture Notes in Phys (Vol. 239) (Berlin: Springer) pp263−302

    [23]

    Quispel G R W, Nijhoff F W, Capel H W, van ver Linden J 1984 Physica A 125 344

    [24]

    Fokas A S, Ablowitz M 1981 Phys. Rev. Lett. 47 1096

    [25]

    Levi D, Benguria R 1980 Proc. Natl. Acad. Sci. U.S.A. 77 5025

    [26]

    Takahashi D, Satsuma J 1990 J. Phys. Soc. Jpn. 59 3514

    [27]

    Tokihiro T, Takahashi D, Matsukidaira J, Satsuma J 1996 Phys. Rev. Lett. 76 3247

    [28]

    Grammaticos B, Ramani A, Papageorgiou V G 1991 Phys. Rev. Lett. 67 1825

    [29]

    Ramani A, Grammaticos B, Hietarinta J 1991 Phys. Rev. Lett. 67 1829

    [30]

    Hietarinta J, Viallet C 1998 Phys. Rev. Lett. 81 325

    [31]

    Bellon M P, Viallet C 1999 Commun. Math. Phys. 204 425

    [32]

    Sakai H 2001 Commun. Math. Phys. 220 165

    [33]

    Bobenko A I, Suris Yu B 2002 Int. Math. Res. Not. 2002 573

    [34]

    Adler V E, Bobenko A I, Suris Yu B 2003 Commun. Math. Phys. 233 513

    [35]

    Nijhoff F W, Walker A J 2001 Glasg. Math. J. 43A 109

    [36]

    Nijhoff F W 2002 Phys. Lett. A 297 49

    [37]

    Nijhoff F W, Atkinson J, Hietarinta J 2009 J. Phys. A: Math. Theor. 42 404005

    [38]

    Hietarinta J, Zhang D J 2008 J. Phys. A: Math. Theor. 42 404006

    [39]

    Atkinson J, Nijhoff F W 2010 Commun. Math. Phys. 299 283

    [40]

    Nijhoff F W, Atkinson J 2010 Int. Math. Res. Not. 2010 3837

    [41]

    Butler S, Joshi N 2010 Inverse Prob. 26 115012

    [42]

    Butler S 2012 Nonlinearity 25 1613

    [43]

    Cao C W, Xu X X 2012 J. Phys. A: Math. Theor. 45 055213

    [44]

    Cao C W, Zhang G Y 2012 J. Phys. A: Math. Theor. 45 095203

    [45]

    Zhang D J, Zhao S L 2013 Stud. Appl. Math. 131 72

    [46]

    Bobenko A I, Its A 2016 Duke Math. J. 165 2607

    [47]

    Hietarinta J, Joshi N, Nijhoff F W 2016 (Cambridge: Camb. Univ. Press)

    [48]

    Zhang D J, Chen S T 2010 Stud. Appl. Math. 125 393

    [49]

    Zhang D J, Chen S T 2010 Stud. Appl. Math. 125 419

    [50]

    Fu W, Qiao Z J, Sun J W, Zhang D J 2013 arXiv: 1307.3671

    [51]

    Fu W, Qiao Z J, Sun J W, Zhang D J 2015 J. Nonlinear. Math. Phys. 22 321

    [52]

    Wahlquist H D, Estabrook F B 1973 Phys. Rev. Lett. 31 1386

    [53]

    Lamb JR G L 1971 Rev. Mod. Phys. 43 99

    [54]

    Chen H H 1974 Phys. Rev. Lett. 33 925

    [55]

    Orfanidis S J 1978 Phys. Rev. D 18 3828

    [56]

    Bianchi L 1892 Rend. Lincei 5 2

    [57]

    Bianchi L 1894 Lezioni di Geometria Differenziale (3rd Ed.) (Pisa: Enrico Spoerri)

    [58]

    Konopelchenko B G 1982 Phys. Lett. A 87 445

    [59]

    Levi D 1981 J. Phys. A: Math. Gen. 14 1083

    [60]

    Adler V E, Yamilov R I 1994 J. Phys. A: Math. Gen. 27 477

    [61]

    Merola I, Ragnisco O, Tu G Z 1994 Inverse Prob. 10 1315

    [62]

    Zhang H W, Tu G Z, Oevel W, Fuchssteiner B 1991 J. Math. Phys. 32 1908

    [63]

    Chen K, Deng X, Zhang D J 2017 J. Nonlinear. Math. Phys. 24(Suppl.1) 18

    [64]

    Cao C W, Zhang G Y 2012 Chin. Phys. Lett. 29 050202

    [65]

    Walker A J 2001 Ph.D. Thesis (Leeds: University of Leeds)

    [66]

    Adler V E 1998 Int. Math. Res. Not. 1998 1

    [67]

    Hietarinta J 2005 J. Nonlinear. Math. Phys. 12 223

    [68]

    Nijhoff F W, Papageorgiou V G, Capel H W, Quispel G R W 1992 Inverse Prob. 8 597

    [69]

    Nijhoff F W 1997 In: Fokas A S, Gel’fand I M (eds) Algebraic Aspects of Integrable Systems: In memory of Irene Dorfman (Boston: Birkhauser) pp237−260

    [70]

    Nijhoff F W 1999 In: Bobenko A I, Seiler R (eds) Discrete Integrable Geometry and Physics (Oxford: Clarendon Press) pp209−234

    [71]

    Hietarinta J 2011 J. Phys. A: Math. Theor. 44 165204

    [72]

    Hietarinta J, Zhang D J 2008 preprint

    [73]

    Atkinson J 2008 J. Phys. A: Math. Theor. 41 135202

    [74]

    Adler V E, Bobenko A I, Suris Yu B 2009 Funct. Anal. Appl. 43 3

    [75]

    Boll R 2011 J. Nonlinear. Math. Phys. 18 337

    [76]

    Boll R 2012 Ph.D Dissertation (Berlin: Technischen Universität Berlin)

    [77]

    Adler V E, Bobenko A I, Suris Yu B 2012 Int. Math. Res. Not. 2012 1822

    [78]

    Miwa T 1982 Proc. Jpn. Acad. 58A 9

    [79]

    Konopelchenko B G, Schief W K 2002 Stud. Appl. Math. 109 89

    [80]

    Hirota R 1981 J. Phys. Soc. Jpn. 50 3785

    [81]

    Bianchi L 1885 Ann. Matem. 13 177

    [82]

    Atkinson J, Nieszporski M 2014 Int. Math. Res. Not. 2014 4215

    [83]

    Zhang D D, Zhang D J 2018 J. Nonlinear. Math. Phys. 25 34

    [84]

    Zhang D J, Cheng J W, Sun Y Y 2013 J. Phys. A: Math. Theor. 46 265202

    [85]

    Xenitidis P 2011 J. Phys. A: Math. Theor. 44 435201

    [86]

    Bridgman T, Hereman W, Quispel G R W, van der Kamp P H 2013 Found. Comput. Math. 13 517

    [87]

    Hietarinta J, Zhang D J 2010 J. Math. Phys. 51 033505

    [88]

    Hietarinta J, Zhang D J 2011 SIGMA 7 061

    [89]

    Atkinson J, Hietarinta J, Nijhoff F W 2007 J. Phys. A: Math. Theor. 40 F1

    [90]

    Konopelchenko B G, Schief W K 2002 J. Phys. A: Math. Gen. 35 6125

  • [1] 潘军廷, 张宏. 极化电场对可激发介质中螺旋波的控制. 物理学报, 2020, (): . doi: 10.7498/aps.69.20191934
    [2] 刘婉馨, 陈瑞, 刘永杰, 王俊峰, 韩小涛, 杨明. 脉冲强磁场下的电极化测量系统. 物理学报, 2020, 69(5): 057502. doi: 10.7498/aps.69.20191520
    [3] 周瑜, 操礼阳, 马晓萍, 邓丽丽, 辛煜. 脉冲射频容性耦合氩等离子体的发射探针诊断. 物理学报, 2020, (): . doi: 10.7498/aps.69.20191864
    [4] 黄永峰, 曹怀信, 王文华. 共轭线性对称性及其对\begin{document}$ {\mathcal{P}}{\mathcal{T}} $\end{document}-对称量子理论的应用. 物理学报, 2020, 69(3): 030301. doi: 10.7498/aps.69.20191173
    [5] 左富昌, 梅志武, 邓楼楼, 石永强, 贺盈波, 李连升, 周昊, 谢军, 张海力, 孙艳. 多层嵌套掠入射光学系统研制及在轨性能评价. 物理学报, 2020, 69(3): 030702. doi: 10.7498/aps.69.20191446
    [6] 庄志本, 李军, 刘静漪, 陈世强. 基于新的五维多环多翼超混沌系统的图像加密算法. 物理学报, 2020, 69(4): 040502. doi: 10.7498/aps.69.20191342
  • 引用本文:
    Citation:
计量
  • 文章访问数:  106
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-28
  • 修回日期:  2019-11-14
  • 上网日期:  2019-12-05
  • 刊出日期:  2020-01-01

离散可积系统: 多维相容性

    基金项目: 国家自然科学基金(批准号: 11875040, 11631007)

摘要: 对比已有完善而系统理论的微分方程领域, 差分方程理论尚处于发展之中. 近年来离散可积理论的进展, 带来了差分方程理论的革命. 多维相容性是伴随离散可积系统研究出现的新的概念, 作为对离散可积性的一种理解, 提供了构造离散可积系统的Bäcklund变换、Lax对和精确解的工具. 本文旨在综述多维相容性的概念及其在离散可积系统研究中的应用.

English Abstract

    • 离散系统泛指含有离散自变量的常差分、微分差分、偏差分系统、以及变换和映射等. 由于缺少导数、积分等局部化的数学工具, 对于非线性离散系统的研究, 往往伴随着新的数学概念、理论和方法的出现.

      现代可积理论兴起于20世纪60年代中期孤立子的命名[1]、反散射变换方法的建立[2]和Lax对概念的提出[3], 参见文献[4]. 早期对离散可积系统的探索主要是可积离散化, 早在20世纪70年代: Case和Kac[5]对Schrödinger谱问题的离散以及Ablowitz和Ladik[68]对Ablowitz-Kaup-Newell-Segur(AKNS) 谱问题的离散 建立了基于差分的反散射变换过程; Hirota[9-11]对双线性方程的离散利用双线性Bäcklund变换与Lax对的联系获得了一系列离散可积系统. 进入80年代后, 对离散可积系统的研究逐渐向形成独立于连续系统的研究方法, 系统进展有: 京都学派Date等[12-16]和 Ueno等[17]对Sato理论的离散; 荷兰学者Nijhoff等[18-22]和Quispel等[23]基于Fokas和Abolwitz[24]的直接线性化格式以及Levi和Benguria[25]的变换与离散的同等性认识发展起来的构造和研究离散可积系统的系统方法. 进入90年代以后, 离散可积系统的显著进展包括: 超离散可积系统的提出及其连续极限的建立[26,27], 基于奇点囿禁的奇点理论与可积性的联系的发现[28,29], 奇点囿禁在可积性判别上的不充分性的发现[30], 以及基于代数熵理论对可积性的判别[31]. 这些进展相继推动了超离散可积系统、离散Painlevé方程、可积性检验等方面的发展.

      进入新世纪后, 离散可积系统继续迎来新的发展. Sakai[32]基于有理曲面理论和Blow-up分析对离散Painlevé方程的分类, 揭示了离散Painlevé方程丰富的代数几何结构, Bobenko和Suris[33]、Adler等[34]和Nijhoff等[35,36]学者对于“多维相容性”的理解以及对若干离散可积系统的分类, 各种精确求解方法 在离散可积系统中相继实现[3745], 等等, 一系列进展标志着对离散可积系统的研究进入到一个新的阶段. 2009年在著名的英国剑桥牛顿所(The Isaac Newton Institute for Mathematical Sciences)举办的为期半年的离散可积系统主题研讨活动, 是对当时离散可积系统蓬勃发展的一个反映.

      可积系统与数学和物理的众多分支都有联系, 已经渗透到数学物理的各个领域[4]. 在数学方面, 离散可积系统的发展伴随着离散几何和离散复分析的发展, 这与连续意义下可积系统丰富的几何背景以及分析的工具在其研究中的重要作用是分不开的. 对与离散可积系统相关的差分算子理论与复分析方面的研究, 如势函数的渐近性、散射理论、Riemann-Hilbert问题、解的长时间渐近分析、初边值问题等等, 建议参考文献[42,46]. 除了提及的上述进展, 还有如离散几何、对称理论、Galois理论、Lagrangian多形式理论等许多方面的重要进展. 值得指出的是, 最近由Hietarinta, Joshi和Nijhoff [47]完成的一部专著《Discrete Systems and Integrability》已经出版. 读者可以参考了解其更多内容. 国际上每两年一届的SIDE会议(Symmetries and Integrability of Difference Equations)是离散可积系统及其相关问题的首要国际会议, 1994年始于加拿大Montreal, 至今已经成功举办13届, 其中SIDE-10于2012年在我国宁波召开.

      本文将在第二、三节重点介绍多维相容性的概念及其应用. 希望通过具体的例子让更多读者了解离散可积系统及其中的概念和方法.

    • 离散可积系统从未独立于连续系统. 下面首先来介绍两种引入离散变量的途径.

    • 离散化是引入离散变量的方式之一. 从熟悉的AKNS谱问题开始:

      $ { {\varPhi}}_x ={ {M \varPhi}},\; \; { M} = \left(\!\!\!\begin{array}{cc}\eta & u \\ v & -\eta\end{array}\!\!\!\right),\; \; \; { {\varPhi}} = (\phi_1, \phi_2)^{\rm T}, $

      其中$ u = u(x, t), \; v = v(x, t) $, $ \eta $是谱参数, 为了方便, 记$ { M} = { M}(\eta, U) $, $ U = (u, v)^{\rm T} $. 定义

      $ { {\varPhi}}_{n+j} = { {\varPhi}}(x+j\epsilon), $

      然后利用差分$ ({ {\varPhi}}_{n+1}-{ {\varPhi}}_n)/\epsilon $替代导数$ { {\varPhi}}_x $, 得到

      $ \begin{split} & { {\varPhi}}_{n+1} = \left(\!\!\!\begin{array}{cc} 1+\epsilon \eta & \epsilon u \\\epsilon v & 1- \epsilon\eta \!\!\!\end{array} \right) { {\varPhi}}_n,\\ &{ {\varPhi}}_n = (\phi_{1,n}, \phi_{2,n})^{\rm T}. \end{split}$

      Ablowitz-Ladik (AL)谱问题指[7]

      $ \begin{split} &{ {\varPhi}}_{n+1} = \left(\!\!\!\begin{array}{cc} \lambda & Q_n \\ R_n & 1/\lambda \end{array}\!\!\!\right) { {\varPhi}}_n,\\&{ {\varPhi}}_n = (\phi_{1,n}, \phi_{2,n})^{\rm T}. \end{split}$

      在(2)式和

      $ (Q_n, R_n) = \epsilon (u,v),\; \; \; \lambda = {\rm e}^{\epsilon \eta} $

      之下, 对$ { {\varPhi}}_{n+1} = { {\varPhi}}(x+\epsilon) $$ \epsilon = 0 $展开, 并取$ \epsilon \to 0 $, 则AKNS谱问题(1)式可以作为领头项从AL谱问题(4)式中恢复出来, 同时谱参数从$ \lambda $-平面变到$ \eta $-平面.

      (2)式是我们所熟悉的数值(差分)离散: 如图1, 将区间$ [x_0, x] $等分成n份, 步长为$ \epsilon $. 在x点, $ {\varPhi}_n $$ {\varPhi}(x) $仍然表示相同的数值, 但是自变量的空间已经从实数域$ \mathbb{R} $变为整数域$ \mathbb{Z} $. 差分离散是计算连续极限的基础, 但是对于可积系统而言, 差分离散(也称为直接离散)不足以保持原有系统的可积特征. 关于由AL谱问题(4)式引出的半离散AKNS系统, 读者可以参考文献[4851].

      图  1  $[x_0, x]$上的数值离散

      Figure 1.  Numerical discretisation on $[x_0, x]$

      递推关系可以很自然地视为离散的方程. 例如Hermite正交多项式$ \{H_n(x)\} $满足递推关系

      $ H_{n+1}(x) = 2xH_n(x)-2nH_{n-1}(x). $

      第I型Bessel函数

      $ {J_\alpha }(x) = {\left( {\frac{x}{2}} \right)^\alpha }\sum\limits_{k = 0}^\infty {\frac{{{{( - {x^2}/4)}^k}}}{{k!{\mkern 1mu} \Gamma (\alpha + k + 1)}}} $

      满足方程常微分方程

      $ x^2 w''+xw'+(x^2-\alpha^2)w = 0, $

      其中x是自变量, $ \alpha $是参数. 由此方程可以建立Bessel函数的递推关系(参考文献[47]的第37页)

      $ x(w_{\alpha+1}+w_{\alpha-1})-2\alpha w_{\alpha} = 0, $

      此时x$ \alpha $互换了角色, x是参数, $ \alpha $成为自变量. 再如, Painlevé II方程

      $ f{''}(t) = 2f^3+tf-\alpha $

      的解满足递推关系

      $ f_{\alpha+1}(t) = -f_{\alpha}(t)-\dfrac{\alpha+\dfrac{1}{2}}{f^{'}_{\alpha}(t)-f^2_{\alpha}(t)-\dfrac{t}{2}}, $

      这里$ f'(t) $表示对t的导数.

      在上面Bessel函数和Painlevé II方程的解的例子中, 递推关系分别来自于两个微分方程的解之间的变换, 可以视为通过变换建立起来的解之间的叠加关系. 这是由连续系统到离散系统的常见过程. 尤其, KdV方程的Bäcklund变换的非线性叠加公式提供了对KdV方程的完美的离散化.

      对于著名的KdV方程,

      $ u_t+6uu_x+u_{xxx} = 0, $

      Wahlquist和Estabrook[52]在1973年发现形如

      $ (\widetilde w+w)_x = 2\lambda-\frac{1}{2}(\widetilde w -w)^2, \tag{7a}$

      $ (\widetilde w-w)_t = \frac{1}{2}[(\widetilde w-w)^3]_x-6\lambda (\widetilde w-w)_{x}-(\widetilde w-w)_{xxx}\tag{7b} $

      的非线性Bäcklund变换, 其中w满足势KdV方程

      $ w_t+3(w_x)^2+ w_{xxx} = 0, $

      $ u = w_x $满足KdV方程(6).

      利用Bäcklund变换(7a)式可以建立KdV方程的解的非线性叠加公式. 首先在(7a)式中, 从同样的种子解w出发, 分别记由$ \lambda = \lambda_1 $$ \lambda = \lambda_2 $引出的$ \widetilde w $$ w_1 $$ w_2 $, 即

      $ (w_1+w)_x = 2\lambda_1-\frac{1}{2}(w_1 -w)^2,\tag{9a}$

      $ (w_2+w)_x = 2\lambda_2-\frac{1}{2}(w_2 -w)^2.\tag{9b}$

      接下来, 在(7a)式中取$ w = w_1, \; \lambda = \lambda_2 $, 记$ \widetilde w = w_{12} $, 有

      $ (w_{12}+w_1)_x = 2\lambda_2-\frac{1}{2}(w_{12} -w_1)^2;\tag{10a}$

      $ w = w_2, \; \lambda = \lambda_1 $, 记$ \widetilde w = w_{21} $, 有

      $ (w_{21}+w_2)_x = 2\lambda_1-\frac{1}{2}(w_{21} -w_2)^2.\tag{10b}$

      上述过程可以描述为图2.

      图  2  Bäcklund变换解的交换性质

      Figure 2.  Permutation property of Bäcklund transformation

      可以证明$ w_{12} $$ w_{21} $能够相等(参考文献[49]). 进一步, 从(9)式和(10)式得到

      $4(\lambda_1-\lambda_2) = (w_1-w_2)(w_{12}-w),$

      称为(势)KdV方程解的非线性叠加公式, 也称为Bianchi等式, 还称为离散的势KdV方程[47]. 作为离散的方程时, (11)式通常写为

      $ (w_{n+1,m}-w_{n,m+1})(w_{n,m}-w_{n+1,m+1}) = p-q, $

      其中pq分别是对应于n–方向和m–方向的方向参数.

      modified KdV(mKdV)方程

      $ u_t = 6u^2u_x -u_{xxx} $

      和sine-Gordon方程

      $ \theta_{xx}-\theta_{tt} = \sin \theta $

      拥有形式相同的非线性叠加公式(文献[53,54]):

      $ \tan \frac{\theta_{12}-\theta}{4} = \frac{p+q}{p-q}\tan \frac{\theta_{1}-\theta_2}{4}, $

      $ p \sin \frac{\theta_{12}-\theta_1+\theta_2-\theta}{4} = q\sin \frac{\theta_{12}+\theta_{1}-\theta_2-\theta}{4}. $

      此方程经过变换$ u = {\rm e}^{\frac{{\rm i}\theta}{2}} $以后, 可以写为

      $ p(u_{12} u_2 - u_1 u) = q (u_{12} u_1 - u_2 u), $

      此方程也称为离散势mKdV方程. 有意思的是, 离散的sine-Gordon方程形如

      $ p \sin \frac{\varphi_{12}+\varphi-\varphi_1-\varphi_2}{4} = q\sin \frac{\varphi_{12}+\varphi+\varphi_{1}+\varphi_2}{4}, $

      由Hirota[11]和Orfanidis[55]先后得到, 与(16)式只有部分符号差别, 经过同样的变换$ u = {\rm e}^{\textstyle\frac{{\rm i}\varphi}{2}} $以后, 可以写为[18]

      $ p(u_{12} u - u_1 u_2) = q(u u_1u_2 u_{12} -1). $

      2阶AKNS方程组

      $ u_t = u_{xx}-2u^2v,\; \; v_t = -v_{xx}-2v^2u $

      的Bäcklund变换的非线性叠加公式形如

      $ (u_1-u_2)(uv_{12}+1)+(p-q)u = 0, \tag{21a}$

      $ (v_1-v_2)(uv_{12}+1)-(p-q)v_{12} = 0, \tag{21b}$

      由Konopelchenko[58]于1982年获得.

      Bäcklund变换提供了引入离散变量的一种方式. (9a)式中$ w_1 $可视为w沿一个方向上的平移, (9b)式中$ w_2 $可视为w沿另一个方向上的平移. 非线性叠加公式(11)式是这两个方向上平移相容性的结果. 同时, (9a)式作为一个独立方程, 也可视为一个微分-差分方程.

      Bäcklund变换应用于非线性模型, Darboux变换的目标对象则是在线性层面(Lax对). 对于给定的连续谱问题, 存在与之相应的连续等谱发展方程族; 它的Darboux变换可视为一个离散的谱问题, 与原有连续谱问题之间的相容性引出的微分差分方程 自然给出连续等谱发展方程族的解之间的Bäcklund变换. 这一漂亮的联系首先由Levi和Benguria[25,59]发现, 对于研究和理解离散可积系统具有重要意义. AKNS谱问题(1)式存在如下Darboux变换[60]:

      $ \widetilde { {\varPhi}} = { T}{ {\varPhi}},\; { T} = { T}(\gamma,{ U},\widetilde { U}) = \left(\!\!\!\begin{array}{cc} 2(\eta-\gamma) + u \widetilde v & u \\ \widetilde v & 1 \end{array}\!\!\!\right), $

      式中$ { U} = (u, v)^{\rm T} $. 作为Darboux变换, $ \widetilde { {\varPhi}} $满足AKNS谱问题(1)式, 即$ \widetilde { {\varPhi}} _x = \widetilde { M} { {\varPhi}} $, $ \widetilde { M} = { M}(\eta, \widetilde { U}) $, 这相当于要求(1)式与(22)式相容, 即$ (\widetilde { {\varPhi}} )_x = (\widetilde{{ {\varPhi}}}_x) $, 由此有

      $ { T}_x-\widetilde { M}{ T}+{ {TM}} = 0, $

      引出的方程

      $ u_{x} = \widetilde u- u- u^2 \widetilde v,\; \; v_{x} = v - \underset{\sim }{\mathop{v}}\, + \underset{\sim }{\mathop{u}}\, v^2, $

      可视为连续的等谱AKNS方程族的Bäcklund变换. 作为离散的谱问题, (22)式引出一个新的半离散AKNS族[61,62], 可视为微分-差分Kadomtsev-Petviashvili (KP) 方程族Lax三重组的对称约束的结果[63]. 若考虑两个不同方向上的Darboux变换:

      $ \widetilde { {\varPhi}} = { T}(p/{2},{ U},\widetilde { U}){ {\varPhi}},\; \; \widehat{{ {\varPhi}} } = { T}(q/{2},{ U},{\widehat{{ U} }}){ {\varPhi}}, $

      作为离散的Lax对, 其相容性得到方程(21), 且在适当的连续极限下, (21)式引出(20)式(参考文献[43]). 显然, 由Darboux变换作为Lax对引出的4点方程$ \, { Q}({ U}, \widetilde { U}, {\widehat{ U }}, {\widehat{\widetilde{ U }}}; p, q) = 0 $即为4个解之间的非线性叠加公式. 更多的关于Darboux变换作为离散谱问题的例子可见文献[64].

      当作为离散的方程时, 变换的符号意义可以理解为

      $\begin{split} & u = u_{n,m},\\ & \widetilde u = u_{n+1,m},\\ &{\widehat{u }} = u_{n,m+1}, \\ & {\widehat{\widetilde{u }}} = u_{n+1,m+1}, \end{split}$

      变换引入的参数pq分别表示n–方向和m–方向上的链参数. 如图3所示.

      图  3  变换与平面网格

      Figure 3.  Map and lattice

    • 对于离散可积系统的研究需要引入新的概念. 下面介绍离散可积系统的“多维相容性”的概念.

      回到(势)KdV方程解的非线性叠加公式(12)式, 这一叠加关系可以重复下去. 为了方便,我们采用(26)式中的符号, 将(12)式写为

      $ (u-{\widehat{\widetilde{u }}})(\widetilde u- {\widehat{u }}) = p-q. $

      这是一个定义在平面网格上的四方格方程. 我们引入第三个方向l, 该方向上的平移表示为$ \, {\overline{u}} = u_{n, m, l+1} $, 方向参数为r, 如图4所示.

      图  4  相容立方体

      Figure 4.  Consistent cube

      观察方程在$ (^{\widetilde{\; \; \;}}, p) $$ (^{\widehat{\; \; \;}}, q) $互换下保持不变, 要求$ (\, \overline{\; \; }, r) $与另外两个方向之间也存在对称性, 则有

      $ (u-{\widetilde{\overline{u}}})(\widetilde u- {\overline{u}}) = p-r,\tag{28a}$

      $ (u-{\widehat{\overline{u}}})({\overline{u}}- {\widehat{u }}) = r-q,\tag{28b}$

      视(27)式定义在图4中立方体的底面, (28a)式和(28b)式则分别定义在左侧和后侧. 如何体现这种维数扩充后的相容性? 首先, 由上述3个方程有

      $ {\widehat{\widetilde{u }}} = \frac{p-q}{\widehat{u}-\widetilde{u}} + u,\tag{29a}$

      $ {\widehat{\overline{u}}} = \frac{r-q}{{\widehat{u }}-{\overline{u}}}+u,\tag{29b}$

      $ {\widetilde{\overline{u}}} = \frac{p-r}{{\overline{u}}-\widetilde u}+u,\tag{29c}$

      引入方向l以后的相容性体现在$\overline{ (29{\rm{a}})} = \widetilde{ (29{\rm{b}})} = $$ \widehat{ (29{\rm{c}})} $. 将(29a)式和(29b)式代入到$ \widehat{ (29{\rm{c}})} $算得

      $ \widehat{\widetilde{\overline{u}}} = \frac{(q-p)\widetilde u{\widehat{u }}+(r-q){\widehat{u }}{\overline{u}}+(p-r){\overline{u}}\widetilde u}{(q-r)\widetilde u+(p-q){\overline{u}}+(r-p){\widehat{u }}}. $

      此式右端在任意交换$ (^{\widetilde{\; \;\; }}, p) $$ (^{\widehat{\; \;\; }}, q) $$ (\, \overline{\phantom{M}}, r) $下都保持不变, 说明从$ \overline{ (29{\rm{a}})} $$ \widetilde{ (29{\rm{b}})} $也将得到同样的$ \widehat{\widetilde{\overline{u}}} $.

      方程(27)的这种性质称~为该方程的3维相容性, 也称为立方体相容(consistent-around-the-cube (CAC)), 它体现了(27)式自身的性质. 方程(17)也具有同样的性质. 一般地, 对于一个定义在平面四方格上的方程

      ${ Q}(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = 0, $

      将其嵌入到图4中立方体的6个面中, 得到

      $ \begin{split} & { Q}(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = 0, \;\; { Q}({\overline{u}},\overline{\widetilde u},\overline{{\widehat{u }}},\overline{\widehat{\widetilde u}};p,q) = 0,\\ & { Q}(u,\widetilde u,{\overline{u}},\overline{\widetilde u};p,r) = 0, \;\; { Q}({\widehat{u }},\widehat{\widetilde u},\overline{{\widehat{u }}},\overline{\widehat{\widetilde u}};p,r) = 0, \\ & { Q}(u,{\widehat{u }},{\overline{u}},\overline{{\widehat{u }}};q,r) = 0, \;\;{ Q}(\widetilde u,\widehat{\widetilde u},\overline{\widetilde u},\overline{\widehat{\widetilde u}};q,r) = 0. \end{split} $

      如果从上式左侧的3个方程中分别解出$ {\widehat{\widetilde{u }}} $, $ {\widetilde{\overline{u}}} $$ {\widehat{\overline{u}}} $, 代入到右侧的3个方程后得到同样的$ \overline{\widehat{\widetilde u}} $, 则称方程(30)具有3维相容性.

      对于离散系统而言, 上述多维相容性将提供一系列的可积特征, 后文将进一步介绍. 对于一个d-维的离散方程, 如果嵌入到一个$ (d+1) $-维链后, 所有的d-维子链都是相容的, 则称原d-维离散方程具有$ (d+1) $-维相容性. 对于2维离散方程, 这种多维相容性即为3维相容性, 或CAC性质.

    • 本世纪初, 多维相容性逐渐被系统地认识并作为工具应用到离散系统的研究中[65,35,36,33]. 2003年, Adler, Bobenko和Suris (ABS)[34]发表了他们在多维相容性的基础上对四方格方程的分类, 他们得到的方程列表被统称为ABS方程.

      ABS在假设方程(30)具有3维相容性的基础上, 进一步要求(30)式满足:

      (i) 仿线性(affine linear),

      (ii) 在正方形对称群$ D_4 $下不变,

      (iii) 四面体性质(tetrahedron property).

      在前两条要求下方程(30)的一般形式可以设为

      $ \begin{split} \,& k_0 u\widetilde u{\widehat{u }}{\widehat{\widetilde{u }}} + k_1(u\widetilde u {\widehat{u }}+ u \widetilde u {\widehat{\widetilde{u }}} + u{\widehat{u }} {\widehat{\widetilde{u }}}+\widetilde u{\widehat{u }}{\widehat{\widetilde{u }}}) \\ & ~~ + k_2(\widetilde u{\widehat{u }}+u {\widehat{\widetilde{u }}}) + k_3(u\widetilde u+{\widehat{u }}{\widehat{\widetilde{u }}})+ k_4(u{\widehat{u }}+\widetilde u {\widehat{\widetilde{u }}})\\ & ~~+k_5(u+\widetilde u+{\widehat{u }}+{\widehat{\widetilde{u }}})+k_6 = 0,\\[-10pt]\end{split} $

      其中$ k_i $都是待定系数; 四面体性质指$ \widehat{\widetilde{{\overline{u}}}} $只与$ {\widehat{\widetilde{u }}} $$ {\widetilde{\overline{u}}} $$ {\widehat{\overline{u}}} $有关, 而与u无关, 即($ \widehat{\widetilde{{\overline{u}}}} $, $ {\widehat{\widetilde{u }}} $, $ {\widetilde{\overline{u}}} $, $ {\widehat{\overline{u}}} $)四个点在立方体(图4)中构成一个四面体. ABS证明了满足上面3条要求的3维相容方程(30)只有9个(允许存在Möbius变换), 并分别命名为H1, H2, H3(δ), A1(δ), A2, Q1(δ), Q2, Q3(δ), Q4:

      $ {\rm H}1: \ \ (u-\widehat{\widetilde{u}})(\widetilde{u}-\widehat{u})-p+q = 0,\tag{33a}$

      $ {\rm H}2: \ \ (u-\widehat{\widetilde{u}})(\widetilde{u}-\widehat{u})-(p-q)(u+\widetilde{u}+\widehat{u}+\widehat{\widetilde{u}}+p+q) = 0, \tag{33b}$

      $ {\rm H}3(\delta): \ \ p(u\widetilde{u}+\widehat{u}\widehat{\widetilde{u}})-q(u\widehat{u}+\widetilde{u}\widehat{\widetilde{u}})+\delta(p^2-q^2) = 0, \tag{33c}$

      $ {\rm A}1(\delta): \ \ p(u+\widehat{u})(\widetilde{u}+\widehat{\widetilde{u}})\!-\!q(u+\widetilde{u})(\widehat{u}+\widehat{\widetilde{u}})-\delta^2pq(p-q) \!=\! 0,\tag{33d}$

      $ \begin{split}{\rm A}2: \ \ p(1\, & -q^2)(u\widetilde{u}+\widehat{u}\widehat{\widetilde{u}})-q(1-p^2)(u\widehat{u}+\widetilde{u}\widehat{\widetilde{u}})\\ & -(p^2-q^2)(1+u\widetilde{u}\widehat{u}\widehat{\widetilde{u}}) = 0,\\[-10pt]\end{split}\tag{33e} $

      $ \begin{split}{\rm Q}1(\delta): \ \ p(u\, &-\widehat{u})(\widetilde{u}-\widehat{\widetilde{u}})-q(u-\widetilde{u})(\widehat{u}-\widehat{\widetilde{u}})\\ &+\delta^2pq(p-q) = 0,\end{split} \tag{33f}$

      $ \begin{split} {\rm Q}2: \ \ & p(u-\widehat{u})(\widetilde{u}-\widehat{\widetilde{u}})-q(u-\widetilde{u})(\widehat{u}-\widehat{\widetilde{u}}) \\ +\, &pq(p-q)(u+\widetilde{u}+\widehat{u}+\widehat{\widetilde{u}}-p^2+pq-q^2) = 0,\end{split} \tag{33g}$

      $ \begin{split}{\rm Q}3(\delta):\ \ \, & p(1-q^2)(u\widehat{u}+\widetilde{u}\widehat{\widetilde{u}})-q(1-p^2)(u\widetilde{u}+\widehat{u}\widehat{\widetilde{u}}) \\ -\, &(p^2\!-\!q^2)\left(\widetilde{u}\widehat{u}\!+\! u\widehat{\widetilde{u}}\!+\!\frac{\delta^2(1\! - \!p^2)(1\!-\!q^2)}{4pq}\right) \!=\! 0,\end{split}\tag{33h}$

      $ \begin{split}{\rm Q}4: \ \ &{\rm{sn}}(p)(u\widetilde{u}+\widehat{u}\widehat{\widetilde{u}})-{\rm{sn}}(q)(u\widehat{u}+\widetilde{u}\widehat{\widetilde{u}})\\-\,&{\rm{sn}}(p-q)(\widetilde{u}\widehat{u}+u\widehat{\widetilde{u}} ) \\ +\,&{\rm{sn}}(p){\rm{sn}}(q){\rm{sn}}(p-q)(1+k^2u\widetilde{u}\widehat{u}\widehat{\widetilde{u}}) = 0, \end{split} \tag{33i}$

      其中, H1即为离散势KdV方程(27), H3(δ = 0)为离散势mKdV方程(17) ($ u \to i^{n+m} u $), Q1(δ = 0)为离散的Schwarzian KdV方程, 也称为交比(cross-ratio)方程, A1(δ)在$ u\to (-1)^{n+m} u $下即为Q1(δ), A2在$ u\to u^{(-1)^{n+m}} $下为Q3(δ = 0), Q4是著名的Krichever-Novikov方程的Bäcklund变换的非线性叠加公式, 原型中系数用Weierstrass椭圆函数表示[34,66], (33i)中Jacobi椭圆函数的参数化形式来自于Hietarinta[67], 式中k为Jacobi椭圆函数的模(modulus).

    • 离散的势Boussinesq方程是平面9点方程, 定义在$ 3\times 3 $格子上:

      $ \begin{split} &\frac{P-Q}{x_{11}-x_{20}}-\frac{P-Q}{x_{02}-x_{11}}\\ =\, & (x_{00}-x_{21})(x_{10}-x_{01})\\ & -(x_{01}-x_{22})(x_{21}-x_{12}),\end{split}$

      式中$ P, Q $分别表示$ n, m $的方向参数, $ x_{ij} = x_{n+i, m+j} $. 这类方程早期由Nijhoff等[68-70]学者构造并研究. 借助于其他(因)变量, (34)式可以改写为(参考文献[65])

      $ {\widetilde{y}} - x{\widetilde{x}}+z = 0,\; \; {\widehat{y}} - x\widehat{x}+z = 0, \tag{35a}$

      $ y - x\widehat{\widetilde{x}}+\widehat{\widetilde{z}}-\frac{P-Q}{\widetilde{x}-\widehat{x}} = 0. \tag{35b}$

      取向量$ { u} = (x, y, z)^{\rm T} $, 方程可以视为向量意义下的四方格系统, 并具有3维相容性.

      2011年, Hietarinta[71]在(35)式形式的基础上, 寻找了可能的离散Boussinesq型的3维相容系统, 代表方程为(考虑到对称性和Möbius变换):

      $ [{\rm B}-2]:\; \; \widetilde{y} = x \widetilde{x}-z, \; \; \widehat{y} = x \widehat{x}-z, \tag{36a}$

      $ y = x \widehat{ \widetilde{x}}- \widehat{ \widetilde{z}} +b_0( \widehat{ \widetilde{x}}-x)+b_1+\frac{P-Q}{ \widetilde{x}- \widehat{x}},\tag{36b}$

      $ [{\rm A}-2]:\; \;\widetilde{y} = z \widetilde{x}-x, \; \; \widehat{y} = z \widehat{x}-x, \tag{37a}$

      $ y = x \widehat{ \widetilde{z}}-b_0 x+\frac{P \widetilde{x}-Q \widehat{x}}{ \widehat{z}- \widetilde{z}},\tag{37b}$

      $ [{\rm C}-3]:\; \; {\widetilde{y}} \,z = {\widetilde{x}} - x,\; \; {\widehat{y}} \,z = {\widehat{x}} - x, \tag{38a}$

      $ \widehat{\widetilde{z}}\,y = b_0\,x + b_1 + z \frac{P\,{\widetilde{y}}\,{\widehat{z}}-Q\,{\widehat{y}}\,{\widetilde{z}}}{{\widetilde{z}}-{\widehat{z}}},\tag{38b} $

      式中$ P, Q $仍然分别表示$ n, m $的方向参数, $ b_i $都是任意常数, 有的$ b_i $可以通过适当的变换去掉. 显然, [B-2]是(35)式的推广. 文献[71]中还得到[C-4]方程, 现已证明与[C-3]之间存在Möbius变换, 故不再列出.

    • 离散的sine-Gordon方程(19)不是3维相容的, 但是它与离散势mKdV方程(17)一起可以实现立方体相容[72]. 在符号(26)式下, 两个方程分别表示为

      ${\rm B}1(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = p(u {\widehat{\widetilde{u }}} -\widetilde{u} \widehat{u})-q( u\widetilde u{\widehat{u }}{\widehat{\widetilde{u }}} -1) = 0, $

      $ {\rm B}2(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = p({\widehat{u }} {\widehat{\widetilde{u }}} -u \widetilde{u})- q( \widetilde u \,{\widehat{\widetilde{u }}}- u{\widehat{u }}) = 0. $

      方程(39)在变换$ u\to (-1)^n u $下成为(40)式.

      图  5  离散sine-Gordon和势mKdV方程的相容立方体

      Figure 5.  The consistent cube for the discrete sine-Gordon equation and potential mKdV equation

      将方程(39)和方程(40)以如下方式放置于立方体的6个面:

      $ \begin{split} & {{\text{底}}: \;\;{\rm B}1(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = 0,}\\&{{\text{顶}}:\;\;{\rm B}1({\overline{u}},{{\widetilde{\overline{u}}}}, {{\widehat{\overline{u}}}}, {\widehat{{\widetilde{\overline{u}}}}};p,q) = 0,}\\& {{\text{左}}:\;\; {\rm B}1(u,\widetilde u,{\overline{u}}, {{\widetilde{\overline{u}}}};p,r) = 0,}\\&{{\text{右}}: \;\;{\rm B}1({\widehat{u }},\widehat{\widetilde u},{{\widehat{\overline{u}}}}, {\widehat{{\widetilde{\overline{u}}}}};p,r) = 0,}\\& {{\text{后}}: \;\;{\rm B}2(u,{\widehat{u }},{\overline{u}}, {{\widehat{\overline{u}}}};q,r) = 0,}\\&{{\text{前}}: \;\;{\rm B}2(\widetilde u,\widehat{\widetilde u}, {{\widetilde{\overline{u}}}}, {\widehat{{\widetilde{\overline{u}}}}};q,r) = 0.} \end{split}$

      可以检验, 在给定初值$ u, \widetilde u, {\widehat{u }}, {\overline{u}} $以后, $ {\widehat{{\widetilde{\overline{u}}}}} $是唯一的. 这说明两个方程按照如此方式放置于立方体的6个面时, 它们(在一起)具有3维相容性.

      ABS方程和离散Boussinesq型方程都是利用同一方程构成相容立方体, 而上面这个例子说明允许用 不同的方程构成相容立方体. 文献[73]给出了更多的例子. Adler等[74]和Boll[75]对这种情况进行了讨论. 详细结果列于Boll[76]的博士学位论文中.

    • 对于定义在6面体上的8点3维(标量)离散方程, 如图6(a), 目前已知的具有四维相容性(见图7)的方程有[34]

      图  6  定义3维方程的6面体以及8面体

      Figure 6.  Cube and octahedron for 3D equations

      图  7  围绕超立方体的4D相容性

      Figure 7.  4D consistency around the hyper cube

      $ u\widehat{\widetilde{{\overline{u}}}}-\widetilde u{\widehat{\overline{u}}}+ {\widehat{u }}{\widetilde{\overline{u}}} -{\overline{u}} {\widehat{\widetilde{u }}} = 0, \tag{42a}$

      $ \frac{(\widetilde u-{\overline{u}})({\widehat{u }}-\widehat{\widetilde{{\overline{u}}}})}{({\overline{u}}-{\widehat{u }})(\widehat{\widetilde{{\overline{u}}}}-\widetilde u)} = \frac{(u-{\widetilde{\overline{u}}})({\widehat{\widetilde{u }}}-{\widehat{\overline{u}}})}{({\widetilde{\overline{u}}}-{\widehat{\widetilde{u }}})({\widehat{\overline{u}}}-u)} . \tag{42b}$

      对于定义在8面体上的6点3维离散方程, 如图6(b), 要求这些方程具有4维相容性, Adler等[77]于2011年进行了分类. 满足条件的方程共有5个:

      $ {\rm i})\; \; \; \; \; {\widehat{\widetilde{u }}}{\overline{u}} - {\widehat{\overline{u}}} \widetilde u + {\widetilde{\overline{u}}} {\widehat{u }} = 0,\tag{43a}$

      $ {\rm ii})\; \; \; \; \frac{({\widehat{\widetilde{u }}}-{\widehat{u }})({\widehat{\overline{u}}}-{\overline{u}})({\widetilde{\overline{u}}}-\widetilde u)}{({\widehat{\widetilde{u }}} - \widetilde u)({\widehat{\overline{u}}}-{\widehat{u }})({\widetilde{\overline{u}}}-{\overline{u}})} = 1, \tag{43b}$

      $ {\rm iii})\; \; \; ({\widehat{\widetilde{u }}}-{\widehat{\overline{u}}}){\widehat{u }} + ({\widehat{\overline{u}}}- {\widetilde{\overline{u}}}){\overline{u}} + ({\widetilde{\overline{u}}}-{\widehat{\widetilde{u }}}) \widetilde u = 0, \tag{43c}$

      $ {\rm iv}) \; \; \; \frac{({\widehat{\widetilde{u }}}-{\widehat{\overline{u}}})}{{\widehat{u }}} + \frac{({\widehat{\overline{u}}}- {\widetilde{\overline{u}}})}{{\overline{u}}} + \frac{({\widetilde{\overline{u}}}-{\widehat{\widetilde{u }}})}{ \widetilde u} = 0, \tag{43d}$

      $ {\rm v}) \; \; \; \; \frac{({\widehat{\overline{u}}}-{\widetilde{\overline{u}}})}{{\overline{u}}} = {\widehat{\widetilde{u }}} \left( \frac{1}{\widetilde u} - \frac{1}{ {\widehat{u }}}\right). \tag{43e}$

      这些方程可以存在方向参数, 但是这些参数可以通过适当的变换移除. 这些方程均已经出现于早期的文献中. 如: (42a)式是离散的BKP方程, 又称为Miwa方程[78]; (42b)式被认为是离散Schwarzian BKP, 首先从几何中得到[79]; (43a)式由Hirota[80]在1981年给出, 又称为离散广义Toda方程; (43b)式, (43c)式和(43d)式分别称为离散Schwarzian KP、离散势KP和离散modified KP方程, 最早由Nijhoff等系统构造(见文献[20]中(4)式和(23)式和文献[22]中(4.16)式); (43e)式由Date等于1982年获得(见文献[13]中(N-1)方程).

    • “多维相容性”不仅作为一类离散系统的可积性的理解, 还提供了研究手段, 可用于构造离散可积系统的Bäcklund变换、Lax对、精确解等等.

    • Bäcklund变换始于对伪球面的构造. Bianchi[81]首先证明了Bäcklund给出的含任意参数的变换 可以保持sine-Gordon方程不变. 变换不仅是引入离散变量的一个渠道, 也是联系离散与连续的一座桥梁. 通常, auto-Bäcklund变换指同一方程的解之间的变换, nonauto-Bäcklund变换指不同方程的解之间的变换.

      “多维相容性”提供了非常直观的构造auto-Bäcklund变换的方法. 对于任意具有3维相容性的方程(30), 将其置于相容立方体的6个面中, 如图8所示.

      图  8  相容立方体

      Figure 8.  The consistent cube

      底面和顶面的方程分别为

      $ { Q}(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = 0, \tag{44a}$

      $ { Q}({\overline{u}},{{\widetilde{\overline{u}}}},{\widehat{\overline{u}}},\widehat{\widetilde{{\overline{u}}}};p,q) = 0.\tag{44b}$

      当将$ {\overline{u}} $也视为方程(44a)的解时, 相容立方体侧面的两个方程

      $ { Q}(u,\widetilde u,{\overline{u}},{\widetilde{\overline{u}}};p,r) = 0,\tag{45a}$

      $ { Q}(u,{\overline{u}}, {\widehat{u }},{\widehat{\overline{u}}};r,q) = 0.\tag{45b}$

      自然地构成方程(44a)的一个Bäcklund变换. 对于H1方程(33a)来说, 多维相容性为它提供的Bäcklund变换是(28)式.

      即使是对于两个方程合作构成相容立方体时, 上述思想依然有效. 对于离散sine-Gordon方程(39)和离散势mKdV方程(40)而言, 它们合作构成相容立方体(见(41)式). 此时, 左侧和后侧的两个方程, 即

      $ p(u {\widetilde{\overline{u}}} -\widetilde{u} \overline{u})-r( u\widetilde u{\overline{u}}{\widetilde{\overline{u}}} -1) = 0, \tag{46a}$

      $ r({\widehat{u }} {\widehat{\overline{u}}} -u \overline{u})- q( {\overline{u}} \,{\widehat{\overline{u}}}- u{\widehat{u }}) = 0, \tag{46b}$

      构成离散sine-Gordon方程(39)的Bäcklund变换. 不过这个变换没有对称性.

      多维相容性也可以用于构造不同方程间的nonauto-Bäcklund变换. 下面以Q2和Q1(δ)为例来解释如何利用相容立方体建立nonauto-Bäcklund变换. 首先, 将Q2方程(33g)置于相容立方体的6个面. 然后在顶层方程$ {\rm Q}2({\overline{u}}, \widetilde{{\overline{u}}}, \widehat{{\overline{u}}}, \widehat{\widetilde{{\overline{u}}}}; p, q) = $ 0中引入关系

      ${\overline{u}} = \frac{v}{\epsilon}+\frac{\delta^2}{4\epsilon^2}, $

      对方程乘以$ \epsilon^2 $以后取极限$ \epsilon \to 0 $, 顶层方程变为

      $ p(v-\widehat{v})(\widetilde{v}-\widehat{\widetilde{v}})-q(v-\widetilde{v})(\widehat{v}-\widehat{\widetilde{v}})+\delta^2 pq(p-q) = 0, $

      即Q1(δ). 此时, 立方体的左侧和后侧的两个方程, 经过替换(47)式以及取极限以后, 得到

      $ \delta(u-\widetilde u)(v-\widetilde v) = -p[\delta^2(u+\widetilde u)-2 v\widetilde v]+p^2\delta (v+\widetilde v+p\delta)],\tag{48a}$

      $ \delta(u-{\widehat{u }})(v-{\widehat{v}}) = -q[\delta^2(u+{\widehat{u }})-2 v{\widehat{v}}]+q^2\delta (v+{\widehat{v}}+q\delta)],\tag{48b}$

      这就是Q2与Q1(δ)之间的nonauto-Bäcklund变换.

      2008年, Atkinson构造了ABS方程之间一系列nonauto-Bäcklund变换, 详细结果请参考文献[73]中表3. 除了上述方法, Bäcklund变换可以利用Yang-Baxter映射以及ABS方程的分解性质来构造, 相关内容读者可以分别参考文献[82]和文献[83].

    • 众所周知, Bäcklund变换与Lax对存在密切联系, 对于多维相容系统来说, 这种联系更加直接. 下面以H1方程(33a)为例, 利用其3维相容性引出的Bäcklund变换构造它的Lax对.

      将H1的Bäcklund变换(28)式改写为

      $ {\widetilde{\overline{u}}} = \frac{p-r}{{\overline{u}}-\widetilde u}+u,\; \; {\widehat{\overline{u}}} = \frac{r-q}{{\widehat{u }}-{\overline{u}}}+u. $

      $ {\overline{u}} = g/f $, 得到

      $\begin{split} & \frac{\widetilde{g}}{{\widetilde{f}}} = \frac{ug+(p-r-u\widetilde u)f}{g-\widetilde u f},\\ & \frac{\widehat{g}}{{\widehat{f}}} = \frac{ug+(q-r-u {\widehat{u }}) f}{g-{\widehat{u }} f}. \end{split}$

      引入$ { {\varPhi}} = (g, f)^{\rm T} $, 上式改写为

      $ \widetilde{{ {\varPhi}}} = \gamma_1 { M}{ {\varPhi}},\; \; \; \widehat{{\rm {\varPhi}}} = \gamma_2 { {N\varPhi}}, \tag{51a}$

      $ { M} = \left(\!\!\!\begin{array}{cc} u & p-r-u\widetilde u\\ 1 & -\widetilde u \end{array}\!\!\!\right),\; \; \; { N} = \left(\!\!\!\begin{array}{cc} u & q-r-u {\widehat{u }}\\ 1 & -{\widehat{u }} \end{array}\!\!\!\right),\tag{51b} $

      式中$ \gamma_i $可视为 分式线性形式(50)式写成矩阵形式(51)式后保留的原分子、分母的公因子, 一般取为$ \gamma_1 = \dfrac{1}{\sqrt{|{ M}|}}, \; \gamma_2 = \dfrac{1}{\sqrt{|{ N}|}} $, 对于H1来说, $ \gamma_i $可以取为任意常数. 相容性$ {\widehat{\widetilde{{ {\varPhi}} }}} = \widetilde{\widehat{{ {\varPhi}} }} $引出

      $ {\widehat{ M}} { N} = {\widetilde{ N}} { M}, $

      计算后可得H1方程(33a).

      在H1的Lax对(51)式中, 取$ {\varPhi} = (\phi_1, \phi_2)^{\rm T} $, $ \gamma_1 = \gamma_2 = 1 $, 从$ \widetilde { {\varPhi}} = { M}{ {\varPhi}} $中消去$ \phi_1 $, 得到

      $ \widetilde{\widetilde{\phi}}_2 = (u-\widetilde{\widetilde{u}}) \widetilde{\phi}_{2}-(p-r)\phi_{2}, $

      此为H1方程谱问题的标量形式. 由此出发 可以构造H1方程的无穷守恒律[84].

      上面以H1为例给出了利用Bäcklund变换构造Lax对的方法, 事实上, 所有的ABS方程都具有相同结构的Lax对, 且存在一般形式. 具体地, 对于任一ABS方程$ { Q}(u, \widetilde u, {\widehat{u }}, {\widehat{\widetilde{u }}}, p, q) = 0 $, 其Lax对为(参考文献[85])

      $ \widetilde{{ {\varPhi}}} = \gamma_1 { M}{ {\varPhi}},\; \; \; \widehat{{ {\varPhi}}} = \gamma_2 { N}{ {\varPhi}}, \tag{54a}$

      其中

      $\begin{split} & { M} = \left(\begin{array}{cc} -\partial_{{\widehat{u }}}Q& -Q\\ \partial_{{\widehat{\widetilde{u }}}}\partial_{{\widehat{u }}}Q& \partial_{{\widehat{\widetilde{u }}}}Q \end{array} \right)_{{\widehat{u }} = {\widehat{\widetilde{u }}} = 0,\; q = r}, \\&{ N} = \left(\begin{array}{cc} -\partial_{\widetilde u}Q& -Q\\ \partial_{{\widehat{\widetilde{u }}}}\partial_{\widetilde u}Q& \partial_{{\widehat{\widetilde{u }}}}Q \end{array} \right)_{\widetilde u = {\widehat{\widetilde{u }}} = 0,\; p = r},\end{split}\tag{54b}$

      $ \gamma_1 = \sqrt{|{ M}|}, \; \gamma_2 = \sqrt{|{ N}|} $. 显然, 将矩阵谱问题改写成标量形式后, 可以发现所有的ABS方程都具有类似于(53)式的二阶离散谱问题.

      通常, 对于由一个方程构成相容立方体时, 由于方程的对称性, 其Bäcklund变换具有对称性, Lax对也具有对称性, 即在上式M中, 将$ (^{\widetilde{\; \; }}, p) $换为$ (^{\widehat{\; \; }}, q) $, 即得到N. 在文献[86]中搜集了更多的利用3维相容性构造Lax对的例子. 对于由两个方程构成相容立方体时, 例如2.5节中的离散sine-Gordon方程和势mKdV方程, 由它们的相容立方体构造的离散sine-Gordon方程的Bäcklund变换具有非对称性, 由此引出的离散sine-Gordon方程的Lax对也是非对称的:

      $ \widetilde { {\varPhi}} = \frac{-1}{\sqrt{u\widetilde u}}\left(\begin{array}{lr}-p \widetilde u & r\\ ru\widetilde u & -pu \end{array}\right){ {\varPhi}}, \tag{55a}$

      $ \widehat{{ {\varPhi}} } =\frac{1}{\sqrt{u{\widehat{u }}}}\left(\begin{array}{lr} -r u & q u {\widehat{u }}\\ q & -r {\widehat{u }} \end{array}\right){ {\varPhi}}. \tag{55b}$

    • 下面以离散势KdV (即H1)方程(33a)为例, 演示如何通过其Bäcklund变换构造单孤子解(1-soliton solution (1SS)). 首先需要一个种子解. 取

      $ a^2 = -p,\; \; \; b^2 = -q, $

      不难发现

      $ u_0 = a n +b m +\lambda, $

      是方程(33a)的一个解, 其中$ \gamma $是任意常数.

      方程(33a)的Bäcklund变换是(28)式, 即

      $ (u-\widetilde{\overline{u}})(\widetilde{u}-\overline{u}) = -a^2+k^2,\tag{58a}$

      $ (u-{\widehat{\overline{u}}})(\overline{u}-\widehat{u}) = -k^2+b^2. \tag{58b}$

      相应于参数化(56)式, 这里取$ r = -k^2 $. 取u为(57)作为种子解, 设新解$ \overline{u} $形如

      $ \overline{u} = \overline{u}_{0}+v, $

      其中$ \overline{u}_{} $是种子解(57)沿l方向的一次平移:

      $ \overline{u}_{0} = a n +b m +k+\lambda, $

      v是待定函数. 将(59)式和(57)式代入Bäcklund变换(58)式得到

      $ \widetilde{v} = \frac{Ev}{v+F},\; \; \; \widehat{v} = \frac{Gv}{v+H}, $

      其中

      $\begin{split} & E = -(a+k),\; \; F = -(a-k),\\ & G = -(b +k),\; \; H = -(b-k).\end{split}$

      引入$ v = f/g $$ { {\varPhi}} = (f, g)^{\rm T} $, 将(61)改写为

      $ \widetilde{ {\varPhi}} = \left(\begin{array}{cc} E & 0\\ 1 & F \end{array} \right){ {\varPhi}},\; \; \; \widehat{{ {\varPhi}}} = \left(\begin{array}{cc} G & 0\\ 1 & H \end{array} \right){ {\varPhi}}. $

      验算可见成立相容性$ \widehat{\widetilde{{ {\varPhi}}}} = \widetilde{\widehat{{ {\varPhi}}}} $, 否则需要在上式的矩阵前乘上平衡因子(类似于(51)式中的$ \gamma_i $). 进一步求解, 先后得到

      $ f_{n,m} = E^n f_{0,m},\; \; \; g_{n,m} = \frac{E^n-F^n}{-2k}f_{0,m}+F^n g_{0,m}, $

      以及

      $ f_{n,m} = G^m f_{n,0},\; \; \; g_{n,m} = \frac{G^m-H^m}{-2k}f_{n,0}+ H^m g_{n,0}. $

      从而有

      $ {\varPhi}_{n,m} = \left(\begin{array}{cc} E^nG^m & 0\\ \dfrac{E^nG^m-F^nH^m}{-2k} & F^nH^m \end{array} \right){ {\varPhi}}_{0,0}. $

      $ \rho_{n,m}\! =\! \biggl(\frac{E}{F})^n \! \biggl(\frac{G}{H})^m \rho_{0,0} \!=\! \biggl(\frac{a +k}{a -k})^n \!\biggl(\frac{b + k}{b-k})^m \rho_{0,0}, $

      其中$ \rho_{0, 0} $为常数, 从(63)式可得

      $ v_{n,m} = \frac{v_{0,0}\rho_{n,m}/\rho_{0,0} } {1+\dfrac{v_{0,0}}{2k} - \dfrac{v_{0,0}}{2k}\rho_{n,m}/\rho_{0,0}}. $

      重新定义常数$ \rho_{0, 0} $后, 有

      $ v_{n,m} = \frac{-2k \rho_{n,m}}{1+\rho_{n,m}}. $

      代入到(59)式, 离散势KdV方程(33a)的1-孤子解为

      $ u_{n,m}^{1SS} = \overline{u}_{0}+v_{n,m} = an + bm + \lambda +\frac{k(1- \rho_{n,m})}{1+\rho_{n,m}}. $

      依据上述过程, 若利用Bäcklund 变换获得2-孤子解显然不容易. 但是1-孤子解的结构往往可以“暗示”一些2-孤子解的结构信息, 有助于发现合适的变换公式将离散方程双线性化, 并进一步获得N-孤子解. 相关内容可以参考文献[38, 87, 88].

    • 对于(33)式中所列的ABS方程, $ u = 0 $一般都不是解, 而且也不易明显看出一些简单的解. “不动点方法”是在Bäcklund 变换的基础上求解0-孤子解的方法. 从上一节的求解过程来看, 对于方程$ { Q}(u, \widetilde u, {\widehat{u }}, {\widehat{\widetilde{u }}};p, q) = 0 $的Bäcklund变换(3.1), 其中的参数r是孤子参数, 将作为离散谱在$ {\overline{u}} $中引出一个孤子. 现在, 如果孤子参数r在变换中“失效”, 即$ {\overline{u}} = u $, 则有

      $ { Q}(u,\widetilde u, u,\widetilde u;p,r) = 0,\tag{68a}$

      $ { Q}(u, u, {\widehat{u }},{\widehat{u }};r,q) = 0. \tag{68b}$

      由此变换引出的原方程的解应当是最简单的, 可以作为0-孤子解.

      对于H1方程(33a), 有

      $ (\widetilde u-u)^2 = r-p,\quad ({\widehat{u }}-u)^2 = r-q. $

      由此很容易得到

      $ u = an+bm+\gamma,\; \; a^2 = r-p, \; \; b^2 = r-q, $

      即(57)式.

      如果方程$ { Q}(u, \widetilde u, {\widehat{u }}, {\widehat{\widetilde{u }}};p, q) = 0 $在变换$u\to T(u) $下不变, 可在变换(3.1)中取$ {\overline{u}} = T(u) $, 从而有

      $ { Q}(u,\widetilde u, T(u), T(\widetilde u);p,r) = 0, \tag{69a}$

      $ { Q}(u, T(u), {\widehat{u }},T({\widehat{u }});r,q) = 0. \tag{69b}$

      对于Q1(δ)方程(33f), 在$ T(u) = u+c $下不变, 对于的Bäcklund变换(69)式有

      $\begin{split} & r(u-\widetilde u)^2 = p(c^2+\delta^2r(p-r)),\\ & r(u-{\widehat{u }})^2 = q(c^2+\delta^2r(q-r)). \end{split}$

      经过参数化$ (p, q)\to(\alpha, \beta) $:

      $\begin{split} & p = \frac{c^2/r-\delta^2 r}{a^2-\delta^2},\quad q = \frac{c^2/r-\delta^2 r}{b^2-\delta^2},\\ & \alpha: = pa,\quad \beta: = qb, \end{split}$

      可得Q1(δ)的线性0-孤子解:

      $ u_0 = \alpha n+\beta m+\lambda, $

      其中$ \lambda $是任意常数. 此外, Q1(δ)方程(33f), 还在$ T(u) = -u $下不变. 此时的Bäcklund变换(69)式给出

      $\begin{split} & 4pu\widetilde u +r(u-\widetilde u)^2 = \delta^2pr(p-r),\\ & 4q u{\widehat{u }}+ r(u-{\widehat{u }})^2 = \delta^2r(q-r). \end{split}$

      经过参数化$ (p, q)\to(\alpha, \beta) $:

      $ p = \dfrac{1}{2} r(1- \cosh(\alpha')) = -\dfrac{1}{4}r(1-\alpha)^2/\alpha, $

      $ q = \dfrac{1}{2} r(1- \cosh(\beta')) = -\dfrac{1}{4}r(1-\beta)^2/\beta, $

      得到解

      $ u_0 = A \alpha^n \beta^m +B \alpha^{-n}\beta^{-m},\quad AB = \delta^2 r^2/16. $

      上述所描述的过程称为“不动点方法”. 首先用于构造Q4方程的种子解[89], 其后又系统地应用于其他ABS方程0-孤子解的构造[37,38]. 当然, 也可以利用其他的方法构造0-孤子解, 例如, 利用方程间的Bäcklund变换(如(48)式), 从一个方程的0-孤子解得到另一个方程的0-孤子解.

    • 通过简单的描述和例子介绍了多维相容性的概念及其在离散可积系统中的应用. 对于多维相容性, 存在一定的几何背景, 换言之, 经典的初等几何中的点线之间的关系蕴含着若干离散可积系统的多维相容性(如文献[79,90]). 此外, 多维相容性也可以从平面波因子的对称性来理解. Miwa映射提供的离散的平面波因子[78]

      $ \psi =\prod\limits_{i}{{{(1-{{p}_{i}}k)}^{{{n}_{i}}}}} $

      在各个方向上具有对称性; 考虑到不少离散可积系统可以从平面波因子出发通过Cauchy矩阵方法[37,45] 等途径来构造, 不难理解由此获得的离散可积系统具有多方向上的相容性.

      对比连续系统, 目前, 离散可积系统在代数结构和工具、几何背景、离散的分析工具等方面都仍待发展. 例如, 基于谱问题、零曲率方程和Kac-Moody代数的连续可积系统 理论并不适用于离散系统, 许多连续意义下的几何体系仍未实现离散化, 离散的复分析也尚未成熟.

      对于离散可积系统的研究也是一个在不断认识离散可积性、发展研究方法和工具、与新兴数学工具结合(如Tropical几何、Cluster代数等)的过程. 关于离散可积系统中值得关注的发展与研究方向, 建议读者关注SIDE会议的网站.

参考文献 (90)

目录

    /

    返回文章
    返回