In this Letter we propose that for Lax integrable nonlinear partial differential equations the natural concept of weak solutions is implied by the compatibility condition for the respective distributional Lax pairs. We illustrate our proposal by comparing two concepts of weak solutions of the modified Camassa-Holm equation pointing out that in the <i>peakon</i> sector (a family of non-smooth solitons) only one of them, namely the one obtained from the distributional compatibility condition, supports the time invariance of the Sobolev <i>H</i>¹ norm.

Motivated by the paper of Beals, Sattinger and Szmigielski (2000) [3], we propose an extension of the CamassaHolm equation, which also admits the multipeakon solutions. The novel aspect is that our approach is mainly based on classic determinant technique. Furthermore, the proposed equation is shown to possess a nonisospectral Lax pair.

By introducing some special bi-orthogonal polynomials, we derive the so-called discrete hungry quotient-difference (dhQD) algorithm and a system related to the QD-type discrete hungry LotkaVolterra (QD-type dhLV) system, together with their Lax pairs. These two known equations can be regarded as extensions of the QD algorithm. When this idea is applied to a higher analogue of the discrete-time Toda (HADT) equation and the quotientquotient-difference (QQD) scheme proposed by Spicer, Nijhoff and van der Kamp, two extended systems are constructed. We call these systems the hungry forms of the higher analogue discrete-time Toda (hHADT) equation and the quotient-quotient-difference (hQQD) scheme, respectively. In addition, the corresponding Lax pairs are provided.

In two recent papers by Barry (2010) [29] and (2011) [30], it is conjectured that Somos-4 admits a solution expressed in terms of Hankel determinant with its elements satisfying a convolution recursion relation. In this paper, Barrys conjecture on Somos-4 is firstly confirmed. Actually, we present a more generalized result. The proof is mainly based on new findings on properties for so-called BlockHankel determinants. The method can also be used to prove another conjecture proposed by Michael Somos, which has been solved by Guoce Xin.

Skew-orthogonal polynomials (SOPs) arise in the study of the <i>n</i>-point distribution function for orthogonal and symplectic random matrix ensembles. Motivated by the average of characteristic polynomials of the Bures random matrix ensemble studied in Forrester and Kieburg (Commun Math Phys 342(1):151187, 2016), we propose the concept of <i>partial-skew-orthogonal polynomials</i> (PSOPs) as a modification of the SOPs, and then the PSOPs with a variety of special skew-symmetric kernels and weight functions are addressed. By considering appropriate deformations of the weight functions, we derive nine integrable lattices in different dimensions. As a consequence, the tau-functions for these systems are shown to be expressed in terms of Pfaffians and the wave vectors PSOPs. In fact, the tau-functions also admit the multiple integral representations. Among these integrable lattices, some of them are