Abstract
Macroscopic dynamics of soliton gases can be analytically described by the thermodynamic limit of the Whitham equations, yielding an integro-differential kinetic equation for the density of states. Under a delta-functional ansatz, the kinetic equation for soliton gas reduces to a non-diagonalisable system of hydrodynamic type whose matrix consists of several \(2\times 2\) Jordan blocks. Here we demonstrate the integrability of this system by showing that it possesses a hierarchy of commuting hydrodynamic flows and can be solved by an extension of the generalised hodograph method. Our approach is a generalisation of Tsarev’s theory of diagonalisable systems of hydrodynamic type to quasilinear systems with non-trivial Jordan block structure.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Quasilinear systems of the form
have been thoroughly investigated in the literature. Here \(u=(u^1, \dots , u^n)^T\) is a column vector of the dependent variables and V is a \(n\times n\) matrix. The main emphasis has always been on the strictly hyperbolic case where the matrix V has real distinct eigenvalues. Under the additional condition that the Haantjes tensor of matrix V vanishes, any such system can be reduced to a diagonal form,
\(i=1, \dots , n\), in specially adapted coordinates \(r^1, \dots , r^n\) known as Riemann invariants. We recall that the Haantjes tensor \(H^i_{jk}\) of matrix \(V=\{v^i_j(u)\}\) is defined by the formula
where \(N^i_{jk}\) is the Nijenhuis tensor,
Systems of type (2) govern a wide range of problems in pure and applied mathematics, see, e.g., (Tsarev 1985, 1991; Dubrovin and Novikov 1989; Serre 1999, 2000; Sévennec 1994). It was shown by Tsarev (1985, 1991) that under the so-called semi-Hamiltonian constraint,
system (2) possesses infinitely many commuting flows and conservation laws and can be solved by the generalised hodograph method (here \(i\ne j\ne k\) and \(v^i_{r^j}=\partial v^i/\partial r^j\), etc). Recall that commuting flows of diagonal systems (2) are also diagonal,
where the requirement of commutativity, \(r^i_{ty}=r^i_{yt}\), implies the relations (\(i\ne j\)):
Introducing the notation \(a_{ij}=\frac{v^i_{r^j}}{v^j-v^i}\) and rewriting the equations for commuting flows in the form
the requirement of their compatibility, \((w^i_{r^j})_{r^k}=(w^i_{r^k})_{r^j}\), implies the integrability conditions
which are equivalent to the semi-Hamiltonian property (3). Under conditions (5), system (4) for commuting flows possesses infinitely many solutions parametrised by n arbitrary functions of one variable (Tsarev 1985).
In this paper, we study integrability aspects of quasilinear systems (1) whose matrix V consists of n Jordan blocks of size \(2\times 2\):
\(i=1, \dots , n\), where the coefficients \(v^i(r, \eta )\) and \(p^i(r,\eta )\) are functions of the 2n dependent variables \(r=(r^1, \dots , r^n)\) and \(\eta =(\eta ^1, \dots , \eta ^n)\). Systems of this form have appeared recently as delta-functional reductions in the kinetic equation for soliton gas (Pavlov et al. 2012), see below, as well as hydrodynamic reductions in linearly degenerate dispersionless integrable systems in multidimensions, see Sect. 5.2. In Sect. 2, we derive Jordan block analogues of equations for commuting flows (4) and integrability conditions (5). Although obtained in exactly the same way as for diagonalisable systems, these conditions are considerably more complicated.
Block-diagonal systems of type (6) share the following general properties:
-
The Haantjes tensor of system (6) vanishes identically, which makes these systems natural parabolic analogues of diagonal systems (2). Note that systems (2) can be obtained from (6) by setting \(\eta ^i=const\).
-
The class of systems (6) is invariant under changes of variables of the form
$$\begin{aligned} r^i\rightarrow R^i(r^i, \eta ^i), \qquad \eta ^i\rightarrow {{\mathcal {N}}}^i(\eta ^i), \end{aligned}$$(7)where \(R^i\) and \({{\mathcal {N}}}^i\) are arbitrary functions of the indicated arguments. Note that the group preserving the class of diagonal systems (2) is more narrow, generated by transformations of the form \(r^i\rightarrow R^i(r^i)\), functions of one variable only.
-
Matrices V corresponding to systems (6) form a commutative family. In particular, hydrodynamic commuting flows of systems (6) are also of type (6). This observation allows one to develop the integrability theory of such systems in full analogy with Tsarev’s theory of diagonalisable systems (2) (Tsarev 1985, 1991), by requiring the existence of a hierarchy of commuting flows.
Remark
Integrable quasilinear systems of a single Jordan block type (of arbitrary size) have appeared in the literature as degenerations of hydrodynamic systems associated with multi-dimensional hypergeometric functions (Kodama and Konopelchenko 2016), in the context of parabolic regularisation of the Riemann equation (Konopelchenko and Ortenzi 2018) and as reductions in hydrodynamic chains and linearly degenerate dispersionless PDEs in 3D (Pavlov 2018). A connection of such systems with the modified KP hierarchy was established in Xue and Ferapontov (2020).
Our interest in systems (6) stems from the study of El’s integro-differential kinetic equation for dense soliton gas (El 2003; El and Kamchatnov 2005; El et al. 2011):
where \(f(\eta )=f(\eta , x, t)\) is the distribution function and \(s(\eta )=s(\eta , x, t)\) is the associated transport velocity. Here the variable \(\eta \) is a spectral parameter in the Lax pair associated with the dispersive hydrodynamics; the function \(S(\eta )\) (free soliton velocity) and the kernel \(G(\mu , \eta )\) (phase shift due to pairwise soliton collisions) are independent of x and t. The kernel \(G(\mu , \eta )\) is assumed to be symmetric: \(G(\mu , \eta )=G(\eta , \mu )\). Equation (8) describes the evolution of a dense soliton gas and represents a broad generalisation of Zakharov’s kinetic equation for rarefied soliton gas (Zakharov 1971). It has appeared independently in the context of generalised hydrodynamics of multi-body quantum integrable systems (Doyon and Yoshimura 2018). Note that both \(S(\eta )\) and \(G(\mu ,\eta )\) are system specific. In the special case
system (8) was derived in El (2003) as thermodynamic limit of the KdV Whitham equations and generalised in El and Kamchatnov (2005), El and Tovbis (2020) to the NLS case. It was demonstrated in Pavlov et al. (2012) that under a delta-functional ansatz,
system (8) reduces to a \(2n\times 2n\) quasilinear system for \(u^i(x, t)\) and \(\eta ^i(x, t)\),
where \(v^i\equiv -s(\eta ^i, x, t)\) can be recovered from the linear system
To make the paper more self-contained, we sketch the derivation of Eqs. (10), (11) in Appendix. The special choice \(\eta ^i(x, t)=const\) was discussed previously in El et al. (2011). In this case, the last n equations (10) are satisfied identically, while the first n equations constitute an integrable diagonalisable linearly degenerate system. The main aim of this paper is to demonstrate integrability of the full system (10). This task is not so straightforward since system (10) is not diagonalisable; thus, Tsarev’s theory of diagonalisable systems of hydrodynamic type (Tsarev 1985, 1991, see also Dubrovin and Novikov 1989) does not apply, and a suitable generalisation is required. Following (Pavlov et al. 2012), let us introduce the new variables \(r^i\) by the formula
In the dependent variables \(r^i, \eta ^i\), system (10) reduces to block-diagonal form (6) where the coefficients \(v^i\) and \(p^i\) can be expressed in terms of \((r, \eta )-\)variables as follows. Let us introduce the \(n\times n\) matrix \({\hat{\epsilon }}\) with diagonal entries \(r^1, \dots , r^n\) (so that \(\epsilon ^{ii}=r^i\)) and off-diagonal entries \(\epsilon ^{ik}={G(\eta ^i, \eta ^k)}, \ k\ne i\). Note that this matrix is symmetric due to the symmetry of the kernel G. Define another symmetric matrix \({\hat{\beta }}=-\hat{\epsilon }^{-1}\). Explicitly, for \(n=2\) we have
Denote \(\beta _{ik}\) the matrix elements of \({\hat{\beta }}\) (indices i and k are allowed to coincide). Introducing the notation \(\xi ^k(\eta ^k)=-S(\eta ^k)\), we have the following formulae for \(u^i, v^i\) and \(p^i\) (Pavlov et al. 2012):
see Sect. 5.1 for explicit form of \(v^i\) and \(p^i\) in the simplest cases \(n=1, 2\). Here and in what follows we use the notation \(\epsilon ^{ki}_{,\eta ^i}\) to indicate partial derivative with respect to \(\eta ^i\). Although the explicit form of \(v^i\) and \(p^i\) gets increasingly more complicated as n grows, there exists a remarkably compact formula for the general solution of system (6) corresponding to reduction (10) that works for arbitrary n:
here \(\mu ^i(\eta ^i)\) are arbitrary functions of their arguments and the functions \(\varphi ^i(\eta ^1, \dots , \eta ^n)\) satisfy the relations \(\varphi ^i_{,\eta ^j}=\epsilon ^{ji}(\eta ^i, \eta ^j)\, \mu ^j(\eta ^j)\), \(i\ne j\), no summation. The last n equations (13) define \(\eta ^i(x, t)\) as implicit functions of x and t; then, the first n equations define \(r^i(x, t)\) explicitly. Formula (13) results from the generalised hodograph approach outlined in Sect. 4.
We demonstrate integrability of system (10) by explicitly presenting its commuting flows (Sect. 2) and conservation laws (Sect. 3) parametrised by 2n arbitrary functions of one variable. We emphasise that integrability of system (10) holds for arbitrary functions \(S(\eta )\) and \(G(\mu , \eta )\) in the kinetic equation (8), even without the assumption of symmetry of the kernel G. Finally, we point out that formula (13) does not possess any obvious limit to the case \(\eta ^i=const\) discussed in El et al. (2011). The main reason for this is the structure of commuting flows of systems (6), see Sect. 2.2.
2 Commuting Flows
Given system (6), let us look for commuting flows in the form
The requirement of commutativity of (6) and (14), \(r^i_{ty}=r^i_{yt}\) and \(\eta ^i_{ty}=\eta ^i_{yt}\), leads to the two groups of conditions. First of all, for every \(i\in \{1, \dots , n\}\) one has the relations
Secondly, for every \(i\ne j\in \{1, \dots , n\}\) one has the relations
Here \(a_i, b_i, a_{ij}, b_{ij}, c_{ij}, d_{ij}\) in (15) and (16) is just the notation for the coefficients shared by commuting flows of the hierarchy. Commutativity conditions (15), (16) are presented in symmetric form where the left-hand sides (respectively, right-hand sides) are expressed in terms of system (6) (respectively, its commuting flow (14)).
2.1 Integrability Conditions
Integrability conditions of system (6) can be obtained by calculating consistency conditions of the linear system (15), (16) governing commuting flows. Thus, we regard \(v^i, p^i\) as given and \(w^i, q^i\) as the unknowns. First of all, for every \(i\in \{1, \dots , n\}\) one has the relations
Secondly, for every \(i\ne j\in \{1, \dots , n\}\) one has the relations
Note that the relations \((w^i_{r^i})_{\eta ^i}=(w^i_{\eta ^i})_{r^i}\) give second-order derivatives \(q^i_{r^i r^i}\):
Expressed in terms of the coefficients \(a_i, b_i, a_{ij}, b_{ij}, c_{ij}, d_{ij}\), the consistency conditions of relations (17), (18) and (19) are the required integrability conditions for system (6). There are two types of integrability conditions, involving two and three distinct indices, respectively.
Two-index conditions There are several groups thereof, each involving two distinct indices \(i\ne j\).
The consistency condition \((w^i_{r^i})_{r^j}=(w^i_{r^j})_{r^i}\) implies
The consistency condition \((w^i_{r^i})_{\eta ^j}=(w^i_{\eta ^j})_{r^i}\) implies
The consistency condition \((w^i_{\eta ^i})_{r^j}=(w^i_{r^j})_{\eta ^i}\) implies
The consistency condition \((w^i_{\eta ^i})_{\eta ^j}=(w^i_{\eta ^j})_{\eta ^i}\) implies
The consistency condition \((w^i_{r^j})_{\eta ^j}=(w^i_{\eta ^j})_{r^j}\) implies
The consistency condition \((q^i_{r^j})_{\eta ^j}=(q^i_{\eta ^j})_{r^j}\) implies
Finally, direct calculation shows that the remaining consistency conditions between (19) and the last two equations (18), namely, \((q^i_{r^j})_{r^ir^i}=(q^i_{r^ir^i})_{r^j}\) and \((q^i_{\eta ^j})_{r^ir^i}=(q^i_{r^ir^i})_{\eta ^j}\), are satisfied identically. Relations (20)–(25) form a complete set of integrability conditions for \(4\times 4\) systems (6) (case of two Jordan blocks). For more than two Jordan blocks, additional three-index conditions are required.
Three-index conditions There are several groups thereof, each involving three distinct indices \(i\ne j\ne k\).
The consistency condition \((w^i_{r^j})_{r^k}=(w^i_{r^k})_{r^j}\) implies
The consistency condition \((w^i_{r^j})_{\eta ^k}=(w^i_{\eta ^k})_{r^j}\) implies
The consistency condition \((w^i_{\eta ^j})_{\eta ^k}=(w^i_{\eta ^k})_{\eta ^j}\) implies
The consistency condition \((q^i_{r^j})_{r^k}=(q^i_{r^k})_{r^j}\) implies
The consistency condition \((q^i_{r^j})_{\eta ^k}=(q^i_{\eta ^k})_{r^j}\) implies
The consistency condition \((q^i_{\eta ^j})_{\eta ^k}=(q^i_{\eta ^k})_{\eta ^j}\) implies
Relations (20)–(31) form a complete set of integrability conditions for systems (6). They can be viewed as Jordan block analogues of the semi-Hamiltonian property (5). For integrable systems (6), the general solution of equations (17) – (19) for commuting flows depends on 2n arbitrary functions of one variable.
For system (10) governing reductions in the soliton gas equation, all integrability conditions can be verified by direct calculation (which is quite cumbersome, see Sect. 5.1). In fact, in this particular case commuting flows can be found explicitly: the general commuting flow (14) of system (6) corresponding to reduction (10) has the coefficients
where \(\mu ^i(\eta ^i)\) are n arbitrary functions of one variable and the functions \(\varphi ^i(\eta ^1,\dots ,\eta ^n)\) satisfy the relations \(\partial _{\eta ^{j}}\varphi ^{i}=\epsilon ^{ji}\mu ^j, \quad j\ne i\), no summation (same functions as in (13)). The general commuting flow depends on 2n arbitrary functions of one variable: n functions \(\mu ^i(\eta ^i)\), plus extra n functions coming from \(\varphi ^i\). This demonstrates integrability of the system in question.
2.2 Genuinely Nonlinear and Linearly Degenerate Systems
The analysis of consistency conditions derived in Sect. 2.1 naturally splits into two cases, depending on whether the coefficients \(a_i\) are zero or not. In view of the first relation (15), this is equivalent to \(v^i_{r^i}\ne 0\) (genuine nonlinearity) or \(v^i_{r^i}=0\) (linear degeneracy).
Genuinely nonlinear case \(a_i\ne 0\). Keeping in mind the consistency conditions \(a_{i, r^j}=a_{i, \eta ^j}=0\), one can set \(a_i=1\) via a suitable equivalence transformation (7). In this case the first equation (17) implies \(q^i=w^i_{r^i}\), so that the second equation (17) and the first two equations (18) take the form
while the last two equations (18) will be satisfied identically due to the consistency conditions (20), (21). Thus, in the genuinely nonlinear case commuting flows can be parametrised in the form
where \(w^i\) solve the linear system (33). Note that the first equation (33) is the Lax operator of the modified KP (mKP) hierarchy. Thus, genuinely nonlinear integrable systems of type (14) are governed by the mKP hierarchy. The compatibility conditions of system (33) can be written in the form
where \( {b_{ij}}=b_{j}a_{ij}-{a_{ij}^{2}}-a_{ij, r^j}. \) This provides a compact formulation of integrability conditions (20)–(31).
Linearly degenerate case \(a_i= 0\). Then, the first relation (15) implies \(w^i_{r^i}=0\); thus, commuting flows of linearly degenerate systems (6) are also linearly degenerate. Note that this is the case for system (10) governing reductions in the soliton gas equation. This explains the observation that formula (13) does not possess any obvious limit to the case \(\eta ^i=const\); indeed, in the diagonal case, generic commuting flows of linearly degenerate semi-Hamiltonian systems are not linearly degenerate.
For linearly degenerate systems, integrability conditions (20)–(31) constitute a Darboux-integrable system which can be ‘solved explicitly’, we hope to report the details elsewhere.
3 Conservation Laws
Conservation laws provide an alternative way to derive integrability conditions. Recall that conservation laws of system (6) are relations of the form
here \(h(r, \eta )\) is the conserved density and \(g(r, \eta )\) is the corresponding flux. The requirement that relation (34) holds identically modulo (6) implies the equations
Eliminating the flux g, we obtain a system of second-order linear PDEs for the density h:
here \(i\ne j\) and the coefficients \(a_i, b_i, a_{ij}, b_{ij}, c_{ij}, d_{ij}\) are the same as in Sect. 2, which confirms the well-known fact that commuting flows share the same conserved densities. One can show that the compatibility conditions of equations (35) coincide with integrability conditions (20)–(31) obtained in Sect. 2. For integrable systems (6), linear system (35) is in involution and its general solution h depends on 2n arbitrary functions of one variable.
Remark
In the genuinely nonlinear case, \(a_i=1\), the first two equations (35) reduce to
while the last two equations (35) will be satisfied identically due to integrability conditions (20)–(31). Thus, conserved densities h are governed by the adjoint Lax pair (36) of the mKP hierarchy, compared with (33).
System (10) possesses infinitely many conservation laws of the form
where \(f^i(\eta ^i)\) are arbitrary functions of their arguments. Although these conservation laws involve only n (rather than the required 2n) arbitrary functions of one variable, their existence already implies the integrability of system (10), as well as of the corresponding system (6). This is a consequence of the following general proposition.
Proposition 1
Suppose that system (6) possesses 2n conservation laws with functionally independent densities. Then, the system is integrable (and automatically possesses infinitely many conservation laws parametrised by 2n arbitrary functions of one variable).
Proof
This statement is a straightforward generalisation of the analogous fact known for diagonalisable systems: if system (2) possesses n conservation laws with functionally independent densities, then it is semi-Hamiltonian (and, therefore, possesses infinitely many conservation laws parametrised by n arbitrary functions of one variable), see (Sévennec 1994). In the present context, a proof can be summarised as follows. Note that all compatibility conditions of system (35) for conserved densities h have the form
where \(\varphi _i\) and \(\psi _i\) depend on the coefficients \(a_i, b_i, a_{ij}, b_{ij}, c_{ij}, d_{ij}\) and partial derivatives thereof. We emphasise that the only ‘free’ second-order partial derivative \(h_{r^i\eta ^i}\) will not enter any of these conditions. Requiring that all compatibility conditions are satisfied identically, that is, \(\varphi _i=\psi _i=0\), we recover the full set of integrability conditions (20)–(31). It remains to note that it is sufficient to require the existence of 2n functionally independent densities h to reach the same conclusion, indeed, any non-trivial compatibility condition imposes a linear homogeneous relation on the gradient of h, thus contradicting functional independence.
Setting in (37) \(f^i=1\) and \(f^i=\eta ^i\), we obtain 2n functionally independent densities \(u^i\) and \(u^i\eta ^i\). Thus, by the above proposition, system (10), as well as the corresponding system (6), are integrable: all integrability conditions (20)–(31) are satisfied identically. This proof of integrability can be viewed as alternative to the one from Sect. 2.1.
In fact, all conservation laws of system (10) can be found explicitly: the general conservation law has the form
here \(\sigma ^i(\eta ^i)\) are arbitrary functions of one variable, the functions \(\tau ^i(\eta ^i)\) can be recovered from the equations \((\tau ^i)'=(\sigma ^i)'\xi ^i\) and the functions \(\psi ^i(\eta ^1, \dots , \eta ^n)\) satisfy the equations \(\psi ^i_{,\eta ^j}=(\sigma ^j)'\epsilon ^{ij}, \ j\ne i\). The general conservation law (38) depends on 2n arbitrary functions of one variable: n functions \(\sigma ^i(\eta ^i)\), plus extra n functions coming from \(\psi ^i\). Setting in (38) \(\sigma ^i=0\) we recover conservation laws (37). \(\square \)
4 Generalised Hodograph Formula
Consider system (1) together with its commuting flow,
where V(u) and W(u) are \(n\times n\) matrices (the commutativity conditions \(u_{ty}=u_{yt}\) impose differential constraints on V and W). Then the matrix relation
where I is the \(n\times n\) identity matrix, defines an implicit solution u(x, t) of system (1) (Tsarev 1985, 1991). Note that, due to the commutativity conditions, only n out of \(n^2\) relations (39) will be functionally independent. Since any integrable system (1) possesses infinitely many commuting flows parametrised by n arbitrary functions of one variable, formula (39) provides a general (multivalued) solution. For commuting systems (6) and (14), the hodograph formula (39) specialises to
which is a system of 2n implicit relations for the 2n dependent variables \(r, \eta \).
Proposition 2
Substituting into formula (40) \(v^i, p^i\) from (12) and \(w^i, q^i\) from (32), one obtains formula (13) for the general solution of system (6) corresponding to reduction (10) of the soliton gas equation:
Proof
Using
the relations \(w^i=x+v^i\, t\) take the form
equivalently,
which gives the second half of the required equations due to non-degeneracy of matrix \({\hat{\beta }}\).
Similarly, using
the relations \(q^i=p^i\, t\) take the form
Using \(w^{i}=x+tv^{i}\), both sums in the above relation cancel out, leaving
the first half of the required equations. \(\square \)
5 Examples
In Sect. 5.1, we give some more details on integrable system (6) coming from reductions in the kinetic equation for soliton gas. In Sect. 5.2, we discuss hydrodynamic reductions of type (6) of the 3D Mikhalev system.
5.1 Reductions in the Kinetic Equation for Soliton Gas
Let us begin with the simplest cases \(n=1\) and \(n=2\).
Case \(n=1\). In this case, system (6) is a single \(2\times 2\) Jordan block,
recall that \(\xi =-S(\eta )\). The general commuting flow of this system has the form
where \(\mu (\eta )\) and \(\varphi (\eta )\) are two arbitrary functions. The corresponding generalised hodograph formula (40) simplifies to
here the second relation specifies \(\eta \) as an implicit function of x and t; the first relation gives an explicit formula for r.
Case \(n=2\). The corresponding \(4\times 4\) system (6) has two Jordan blocks; setting \(\epsilon ^{12}=\epsilon \) in (12) we have
and
where \(\xi ^1=\xi ^1(\eta ^1), \ \xi ^2=\xi ^2(\eta ^2)\) and \(\epsilon (\eta ^1, \eta ^2)={G(\eta ^1, \eta ^2)}\). Direct computation of the invariants from Sect. 2 gives
note that the ‘soliton velocities’ \(\xi ^1\) and \(\xi ^2\) do not enter these expressions. In fact, modulo equivalence transformations (7), one can normalise \(\xi ^i(\eta ^i)\equiv \eta ^i\). Direct calculation shows that the corresponding \(4\times 4\) system is integrable: it satisfies all two-index integrability conditions (20)–(25). Note that the same conclusion is true in the more general case \(G(\eta ^1, \eta ^2)\ne G(\eta ^2, \eta ^1)\), where the matrix \(\epsilon \) is not symmetric.
For \(n=2\), equations (17), (18) for commuting flows are straightforward to solve. The general commuting flow is given by the formulae
here \(\varphi ^1(\eta ^1, \eta ^2)\) and \(\varphi ^2(\eta ^1, \eta ^2)\) are two functions such that \(\varphi ^1_{,\eta ^2}=\epsilon \, \mu ^2\) and \(\varphi ^2_{,\eta ^1}=\epsilon \, \mu ^1\), where \(\mu ^1(\eta ^1)\) and \(\mu ^2(\eta ^2)\) are arbitrary functions of the indicated arguments. Thus, the general commuting flow depends on four arbitrary functions of a single argument: \(\mu ^1(\eta ^1)\), \(\mu ^2(\eta ^2)\), plus two more functions coming from the equations for \(\varphi ^1\) and \(\varphi ^2\). The corresponding generalised hodograph formula (40) simplifies to
here the last two relations specify \(\eta ^1\) and \(\eta ^2\) as implicit functions of x and t; the first two relations give explicit formulae for \(r^1\) and \(r^2\). This is a special case of the general formula (13).
For \(n=2\), Eq. (35) for conserved densities is also straightforward to solve: a general conserved density has the form
where \(\sigma ^1(\eta ^1)\) and \(\sigma ^2(\eta ^2)\) are arbitrary functions of the indicated arguments and \(\psi ^1(\eta ^1, \eta ^2)\) and \(\psi ^2(\eta ^1, \eta ^2)\) are two functions such that \(\psi ^1_{,\eta ^2}=(\sigma ^2)'\epsilon \) and \(\psi ^2_{,\eta ^1}= (\sigma ^1)'\epsilon \). Thus, conserved densities depend on four arbitrary functions of one variable. Note that the above conserved density h and the corresponding flux g can be written in a simple form
where \((\tau ^1)'=(\sigma ^1)'\xi ^1, \ (\tau ^2)'=(\sigma ^2)'\xi ^2\). Formulae (42) and (44) generalise to the case of arbitrary n in a straightforward way, leading to the general commuting flows (32) and the general conservation laws (38); in fact, this is the way formulae (32) and (38) were originally found.
Case \(n>2\). Using computer algebra one can show that for \(n=3\) and \(n=4\), which correspond to \(6\times 6\) and \(8\times 8\) systems (6), all integrability conditions (20)–(31) are also satisfied identically (we do not present the explicit forms of \(v^i\) and \(p^i\) due to their complexity). A general proof (alternative to the ones given in Sects. 2 and 3) that works for arbitrary n is based on the following explicit formulae for the coefficients \(a_i, b_i, a_{ij}, b_{ij}, c_{ij}, d_{ij}\):
Since all these coefficients are expressed in terms of \(u^i\) and \(\beta _{ij}\), the verification of integrability conditions (20)–(31), which is a direct calculation, requires the following differentiation rules:
Let us give a direct proof that formula (13) provides a general solution of reduction (10) of the soliton gas equation (which would not appeal to the generalised hodograph method).
Proposition 3
The general solution of system (6) corresponding to reduction (10) of the soliton gas equation is given by formula (13):
here \(\mu ^i(\eta ^i)\) are arbitrary functions of their arguments and the functions \(\varphi ^i(\eta ^1, \dots , \eta ^n)\) satisfy the relations \(\varphi ^i_{,\eta ^j}=\epsilon ^{ji}(\eta ^i, \eta ^j)\, \mu ^j(\eta ^j)\), \(i\ne j\), no summation.
Proof
Differentiating the last set of relations (13) by x and t, we obtain
respectively. Using \(\varphi _{,\eta ^{i}}^{i}-t(\xi ^{i})^{\prime }=r^{i}\mu ^{i}\), we can rewrite these relations in the form
Taking into account that \(\epsilon ^{ii}=r^{i}\), we come to
Contracting these formulae (in the index i) with the inverse matrix \(\beta _{ik}\) and using \(\epsilon ^{ji}\beta _{ik}=-\delta ^j_k\), we obtain
which readily implies \(\eta _{t}^{k}=v^{k}\eta _{x}^{k}\). \(\square \)
Similarly, differentiating \(r^{i}\mu ^{i}=\varphi _{,\eta ^{i}}^{i}-t(\xi ^{i})^{\prime }\) by x and t we obtain (using \(\varphi _{,\eta ^{i}\eta ^{j}}^{i}=\varphi _{,\eta ^{j}\eta ^{i}}^{i}=\epsilon _{,\eta ^{i}}^{ji}\mu ^{j}\)):
respectively. Equivalently,
Replacing \(\eta _{t}^{i}\) by \(v^{i}\eta _{x}^{i}\) and subtracting the first of the above equations (multiplied by \(v^{i}\)) from the second, we obtain:
To verify the equations \(r_{t}^{i}=v^{i}r_{x}^{i}+p^{i}\eta _{x}^{i}\), it remains to show that the following identity holds:
This directly follows from the formula \(\mu ^{i}\eta _{x}^{i}=-u^{i}\), see (45); indeed, under this substitution the above formula reduces to
which identically coincides with (12).
Remark
Reductions in the kinetic equation discussed in this paper satisfy the conditions \(a_i=0\), which are equivalent to \(v^i_{r^i}=0\). These can be seen as Jordan block analogues of the conditions of linear degeneracy which are known to prevent breakdown of classical solutions. Thus, Eq. (8) for soliton gas should be viewed as a linearly degenerate kinetic equation and is expected to have good ‘regularity properties’.
5.2 Reductions in the Mikhalev System
Integrable systems of type (6), as well as integrable quasilinear systems (1) with a more general Jordan block structure, typically arise as hydrodynamic reductions of multi-dimensional linearly degenerate integrable PDEs. As a simple illustration, here we consider the linearly degenerate system
which apparently first appeared in Mikhalëv (1992) in the context of Hamiltonian formalism of KdV type hierarchies and was subsequently investigated in the framework of multi-dimensional dispersionless integrability, see e.g. (Pavlov 2018; Ferapontov and Moss 2015). Hydrodynamic reductions are sought in the form
where the variables \(r=(r^1, \dots , r^n)\) and \(\eta =(\eta ^1, \dots , \eta ^n)\) satisfy a pair of commuting systems (6), (14). Remarkably, all such reductions can be found explicitly: they are parametrised by n arbitrary functions of one variable. The substitution of (47) into (46) gives
as well as the dispersion relations,
Substituting (49) into commutativity conditions (15) and using (48), one obtains the relations
Substituting (49) into commutativity conditions (16) and using (48), one obtains the relations
\(i\ne j\). Furthermore, calculating consistency conditions of equations (48) we obtain
\(i\ne j\) (the last remaining consistency condition, \(n_{r^i\eta ^i}=n_{\eta ^i r^i}\), will be taken into account later). Equations (50), (51) and (52) are straightforward to solve: they imply
Without any loss of generality, modulo equivalence transformations (7), one can set \(m=\sum _k r^k\), \(v^i=m+\eta ^i\). Then the last remaining consistency condition, \(n_{r^i\eta ^i}=n_{\eta ^i r^i}\), implies \(p^i_{r^i}=1\), which we will solve in the form \(p^i=r^i+({s^i}(\eta ^i))'\). Ultimately, we arrive at the following formulae:
here \(s^k(\eta ^k)\) are arbitrary functions of their arguments. Integrability of these reductions can be verified by direct calculation.
References
Bogoyavlenskij, O.I.: Block-diagonalizability problem for hydrodynamic type systems. J. Math. Phys. 47(6), 063502 (2006)
Doyon, B., Yoshimura, T.: Soliton gases and generalized hydrodynamics, J.-S. Caux. Phys. Rev. Lett. 120, 045301 (2018)
Dubrovin, B.A., Novikov, S.P.: Hydrodynamics of weakly deformed soliton lattices: differential geometry and Hamiltonian theory. Russ. Math. Surv. 44, 35–124 (1989)
El, G.A., Kamchatnov, A.M., Pavlov, M.V., Zykov, S.A.: Kinetic equation for a soliton gas and its hydrodynamic reductions. J. Nonlinear Sci. 21(2), 151–191 (2011)
El, G.A.: The thermodynamic limit of the Whitham equations. Phys. Lett. A 311(4–5), 374–383 (2003)
El, G.A., Kamchatnov, A.M.: Kinetic equation for a dense soliton gas. Phys. Rev. Lett. 95, 204101 (2005)
El, G.A., Tovbis, A.: Spectral theory of soliton and breather gases for the focusing nonlinear Schrödinger equation. Phys. Rev. E 101, 052207 (2020)
Ferapontov, E.V., Moss, J.: Linearly degenerate PDEs and quadratic line complexes. Commun. Anal. Geom. 23(1), 91–127 (2015)
Kodama, Yu., Konopelchenko, B.G.: Confluence of hypergeometric functions and integrable hydrodynamic-type systems. Theor. Math. Phys. 188(3), 1334–1357 (2016)
Konopelchenko, B.G., Ortenzi, G.: Parabolic regularization of the gradient catastrophes for the Burgers–Hopf equation and Jordan chain. J. Phys. A 51(27), 275201 (2018)
Mikhalëv, V.G.: Hamiltonian formalism of Korteweg-de Vries-type hierarchies. Funct. Anal. Appl. 26(2), 140–142 (1992)
Pavlov, M.V.: Integrability of exceptional hydrodynamic-type systems. Proc. Steklov Inst. Math. 302(1), 325–335 (2018)
Pavlov, M.V., Taranov, V.B., El, G.A.: Generalized hydrodynamic reductions of the kinetic equation for a soliton gas. Theor. Math. Phys. 171(2), 675–682 (2012)
Serre, D.: Systems of Conservation Laws. 1. Hyperbolicity, Entropies. Shock Waves. Cambridge University Press, Cambridge (1999)
Serre, D.: Systems of Conservation Laws. 2. Geometric Structures, Oscillations, and Initial-Boundary Value Problems. Cambridge University Press, Cambridge (2000)
Sévennec, B.: Géométrie des systèmes hyperboliques de lois de conservation, Mémoire (nouvelle série) N56, Supplément au. Bull. Soc. Math. France 122, 1–125 (1994)
Tsarev, S.P.: Poisson brackets and one-dimensional Hamiltonian systems of hydrodynamic type. Sov. Math. Dokl. 31, 488–491 (1985)
Tsarev, S.P.: The geometry of Hamiltonian systems of hydrodynamic type. The generalized hodograph method. Math. USSR Izvestiya 37, 397–419 (1991)
Xue, L., Ferapontov, E.V.: Quasilinear systems of Jordan block type and the mKP hierarchy. J. Phys. A Math. Theor. 53, 205202 (2020)
Zakharov, V.E.: Kinetic equation for solitons. Sov. Phys. JETP 33, 538–541 (1971)
Acknowledgements
We thank A. Bolsinov, G. El, and D. Tunitsky for clarifying discussions. We also thank the referees for useful comments. The research of EVF was supported by a grant from the Russian Science Foundation No. 21-11-00006, https://rscf.ru/project/21-11-00006/. MVP was supported by the Ministry of Science and Higher Education of the Russian Federation (Agreement No. 075-02-2021-1748).
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Peter Miller.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendix: Derivation of System (10)
Although delta-functional reduction in the kinetic equation for soliton gas was first obtained in Pavlov et al. (2012), to make the paper more self-contained, below we present a short derivation of system (10).
Substituting ansatz (9),
into the first equation (8),
we obtain
Multiplying this equality by a test function \(\varphi (\eta )\) and integrating over \(\eta \) (using integration by parts to eliminate \(\delta '\)) gives
Introducing the notation \(v^i=-s(\eta ^i, x, t)\) and equating to zero the coefficient at \(\varphi (\eta ^i)\), we obtain \(u^i_t=(u^iv^i)_x\). Similarly, equating to zero the coefficient at \(\varphi '(\eta ^i)\) gives \(\eta ^i_t=v^i\eta ^i_x\). Thus, we have recovered equations (10).
Substituting ansatz (9) into the second equation (8),
gives
Setting \(\eta =\eta ^k\) and using \(v^i=-s(\eta ^i, x, t)\), we obtain n linear relations for \(v^i\),
which are identical to (11) under the interchange of indices \(i\leftrightarrow k\).
Concluding Remarks
The approach to hydrodynamic integrability developed in this paper can be generalised to the whole class of non-semisimple quasilinear systems (1) whose matrix V has block-diagonal structure such that:
-
each diagonal block has upper-triangular Toeplitz form as in Xue and Ferapontov (2020), with the unique eigenvalue; here is a \(3\times 3\) Toeplitz matrix with eigenvalue v:
$$\begin{aligned} \left( \begin{array}{ccc} v &{} \quad p &{} \quad q\\ 0&{}\quad v&{}\quad p\\ 0&{}\quad 0&{}\quad v \end{array}\right) . \end{aligned}$$ -
different blocks have different eigenvalues.
Matrices of this type (known as \(\mathrm gl\)-regular) form a commutative family, and the corresponding systems (1) automatically have zero Haantjes tensor. Commuting flows of such systems necessarily have Toeplitz blocks of the same sizes. An \(n\times n\) system of this type should be called integrable if it possesses a hierarchy of hydrodynamic commuting flows/conservation laws parametrised by n arbitrary functions of one variable. In this connection, let us mention that a criterion for block-diagonalisability of quasilinear systems was obtained in Bogoyavlenskij (2006) (theorems 1, 2).
Finally, it would be interesting to consider continuum versions of ansatz (9), system (10), and its general solution (13), with an ultimate goal of establishing integrability of the full kinetic equation (8).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ferapontov, E., Pavlov, M.V. Kinetic Equation for Soliton Gas: Integrable Reductions. J Nonlinear Sci 32, 26 (2022). https://doi.org/10.1007/s00332-022-09782-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00332-022-09782-0
Keywords
- Soliton gas
- Kinetic equation
- Hydrodynamic reduction
- Commuting flow
- Conservation law
- Generalised hodograph formula