A Neumann System of the Third Order Differential Operator Associated with the Boussinesq Equation

Jianli Cao^*, Jingfang Han

College of Science, Henan University of Technology, Zhengzhou, China.

DOI: 10.4236/jamp.2020.812211                

Abstract

Finite-dimensional integrable Hamiltonian system, obtained through the nonlinearization of the 3 × 3 spectral problem associated with the Boussinesq equation, is investigated. A generating function method starting from the Lax-Moser matrix is used to give an effective way to prove the involutivity of integrals. Finite-parameter solution of the Boussinesq equation is calculated based on the commutative system of ordinary differential equations with these integrals as Hamiltonians. The problem of the third order differential operator associated with the Boussinesq Neumann system put forward by H. Flaschka in 1983 is solved.

Keywords

Boussinesq Equation, Conserved Integrals, Lax-Moser Matrix

1. Introduction

Boussinesq equation

$v_{tt}=-\frac{1}{3}{\left(v_{xx}-4v^2\right)}_{xx}$ (1)

is one of the difficult soliton equations, which has been paid common attention in physical and mathematical fields [1] - [6]. In 1983, H. Flaschka put forward a problem of the third order differential operator associated with the Boussinesq Neumann system [7]. Some works focus on the decomposition and the structures of the Modified Boussinesq equation [8] [9] [10] [11]. The decomposition of the Boussinesq Neumann system has not been done thoroughly for a long time. A Neumann system of the Boussinesq equation associated with the third order differential operator is obtained in this paper, which is the extension of the famous KdV Neumann system associated with the second order differential operator. There are many methods to deal with the integrability and involutivity [12] [13] [14]. A generating function method starting from the Lax-Moser matrix [15] - [20] is used to give an effective way to prove the involutivity of integrals.

2. The Generating Function of Integrals

Let $R^{6N}$ be the phase space. The canonical coordinates are

$q^i=\left(q_1^i,\cdots ,q_N^i\right),p^i=\left(p_1^i,\cdots ,p_N^i\right),i=1,2,3.$

Write

$q=\left(q^1,q^2,q^3\right),p=\left(p^1,p^2,p^3\right);q_k=\left(q_k^1,q_k^2,q_k^3\right),p_k=\left(p_k^1,p_k^2,p_k^3\right),$

for short. Take $A=diag\left({\alpha }_1,\cdots ,{\alpha }_N\right)$, where ${\alpha }_1\mathrm{,}\cdots \mathrm{,}{\alpha }_N$ are distinct spectral constants. We denote:

$\langle X,Y\rangle \triangleq {\displaystyle \underset{k=1}{\overset{N}{\sum }}}{X}_k{Y}_k,r^{ij}=\langle q^i,p^j\rangle ,A^{ij}=\langle A{q}^i,p^j\rangle ,$

where X and Y are two N dimensional vectors.

For any matrix $M=\left(M^{ij}\right)$, the element of its adjoint $M^{\mathrm{<mo>\; *\; </mo>}}$ is $M_{ij}^{\mathrm{<mo>\; *\; </mo>}}={\left(-1\right)}^{i+j}{M}_{ji}$, where $M_{ji}$ is the cofactor of the element $M^{ij}$ in M.

In order to proof the Liouville integrability of the Hamiltonian system, consider the Lax-Moser matrix defined as:

$V_{\lambda }=V_{\lambda }\left(p,q\right)={\left(Q_{\lambda }^{ij}\right)}_{3\times 3}+l_{\lambda },$ (2)

where

$Q_{\lambda }^{ij}=Q_{\lambda }\left(q^i,p^j\right)\triangleq {\displaystyle \underset{k=1}{\overset{N}{\sum }}}\frac{q_k^i{p}_k^j}{\lambda -{\alpha }_k},l_{\lambda }=l_{\lambda }\left(p,q\right)=\left(\begin{array}{ccc}0& 0& r^{32}\\ r^{32}& 0& 2r^{12}-r^{31}\\ 0& 0& 0\end{array}\right).$

Let $L_{\lambda \xi }=\xi I_3+V_{\lambda }$. Then we have the generating function of integrals:

$F_{\lambda \xi }=\det L_{\lambda \xi }={\xi }^3+F_1\left(\lambda \right){\xi }^2+F_2\left(\lambda \right)\xi +F_3\left(\lambda \right),$ (3)

where

$\begin{array}{l}{F}_1\left(\lambda \right)=Q_{\lambda }^{11}+Q_{\lambda }^{22}+Q_{\lambda }^{33},\\ F_2\left(\lambda \right)=\left|\begin{array}{cc}{Q}_{\lambda }^{11}& Q_{\lambda }^{12}\\ Q_{\lambda }^{21}+r^{32}& Q_{\lambda }^{22}\end{array}\right|+\left|\begin{array}{cc}{Q}_{\lambda }^{22}& Q_{\lambda }^{23}+2r^{12}-r^{31}\\ Q_{\lambda }^{32}& Q_{\lambda }^{33}\end{array}\right|+\left|\begin{array}{cc}{Q}_{\lambda }^{11}& Q_{\lambda }^{13}+r^{32}\\ Q_{\lambda }^{31}& Q_{\lambda }^{33}\end{array}\right|,\\ F_3\left(\lambda \right)=\left|\begin{array}{ccc}{Q}_{\lambda }^{11}& Q_{\lambda }^{12}& Q_{\lambda }^{13}+r^{32}\\ Q_{\lambda }^{21}+r^{32}& Q_{\lambda }^{22}& Q_{\lambda }^{23}+2r^{12}-r^{31}\\ Q_{\lambda }^{31}& Q_{\lambda }^{32}& Q_{\lambda }^{33}\end{array}\right|.\end{array}$

A series of polynomials $F_{jm}=F_{jm}\left(q\mathrm{,}p\right)$, are defined as the coefficients of the power series expansions as $\left|\lambda \right|>\max \left(\left|{\alpha }_1\right|,\cdots ,\left|{\alpha }_N\right|\right)$ :

$F_j\left(\lambda \right)={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{\lambda }^{-m-1}{F}_{jm},j=1,2,3.$

The first few are expressed as $\left|\lambda \right|>\max \left(\left|{\alpha }_1\right|,\cdots ,\left|{\alpha }_N\right|\right)$ :

$\begin{array}{l}{F}_{10}=r^{11}+r^{22}+r^{33},F_{20}=-3r^{12}{r}^{32},F_{30}={\left(r^{32}\right)}^3,\\ F_{21}=-r^{32}\left(A^{12}+A^{31}\right)-\left(2r^{12}-r^{31}\right)A^{32}+{\displaystyle \underset{1\le i<j\le 3}{\sum }}\left(r^{ii}{r}^{jj}-r^{ij}{r}^{ji}\right),\\ F_{31}={\left(r^{32}\right)}^2{A}^{32}+r^{32}\left(r^{21}{r}^{32}-r^{22}{r}^{31}+r^{13}{r}^{32}-r^{12}{r}^{33}\right)\\ -\left(2r^{12}-r^{31}\right)\left(r^{11}{r}^{32}-r^{12}{r}^{31}\right).\end{array}$

By comparing the coefficients of the same power of $\lambda$, general explicit formulas are obtained:

$F_{1m}=\langle A^m{q}^1,p^1\rangle +\langle A^m{q}^2,p^2\rangle +\langle A^m{q}^3,p^3\rangle ,m\ge 1;$

$\begin{array}{l}{F}_{2m}=-r^{32}\left(\langle A^m{q}^1,p^2\rangle +\langle A^m{q}^3,p^1\rangle \right)-\left(2r^{12}-r^{31}\right)\langle A^m{q}^3,p^2\rangle \\ +{\displaystyle \underset{\begin{array}{c}k+l=m-1\\ k,l\ge 0\end{array}}{\sum }}{\displaystyle \underset{1\le i<j\le 3}{\sum }}\left|\begin{array}{cc}\langle A^k{q}^i,p^i\rangle & \langle A^l{q}^i,p^j\rangle \\ \langle A^k{q}^j,p^i\rangle & \langle A^l{q}^j,p^j\rangle \end{array}\right|,m\ge 2;\end{array}$

$\begin{array}{l}{F}_{3m}={\left(r^{32}\right)}^2\langle A^m{q}^3,p^2\rangle \\ +r^{32}{\displaystyle \underset{\begin{array}{c}k+l=m-1\\ k,l\ge 0\end{array}}{\sum }}\left(\left|\begin{array}{cc}\langle A^k{q}^2,p^1\rangle & \langle A^l{q}^2,p^2\rangle \\ \langle A^k{q}^3,p^1\rangle & \langle A^l{q}^3,p^2\rangle \end{array}\right|-\left|\begin{array}{cc}\langle A^k{q}^1,p^2\rangle & \langle A^l{q}^1,p^3\rangle \\ \langle A^k{q}^3,p^2\rangle & \langle A^l{q}^3,p^3\rangle \end{array}\right|\right)\\ -\left(2r^{12}-r^{31}\right){\displaystyle \underset{\begin{array}{c}k+l=m-1\\ k,l\ge 0\end{array}}{\sum }}\left|\begin{array}{cc}\langle A^k{q}^1,p^1\rangle & \langle A^l{q}^1,p^2\rangle \\ \langle A^k{q}^3,p^1\rangle & \langle A^l{q}^3,p^2\rangle \end{array}\right|\\ +{\displaystyle \underset{\begin{array}{c}j+k+l=m-2\\ j,k,l\ge 0\end{array}}{\sum }}\left|\begin{array}{ccc}\langle A^j{q}^1,p^1\rangle & \langle A^k{q}^1,p^2\rangle & \langle A^l{q}^1,p^3\rangle \\ \langle A^j{q}^2,p^1\rangle & \langle A^k{q}^2,p^2\rangle & \langle A^l{q}^2,p^3\rangle \\ \langle A^j{q}^3,p^1\rangle & \langle A^k{q}^3,p^2\rangle & \langle A^l{q}^3,p^3\rangle \end{array}\right|,m\ge 2.\end{array}$

Expand $F_j\left(\lambda \right)$ in non-negative power of $\lambda$ :

$F_j\left(\lambda \right)={\displaystyle \underset{m=1}{\overset{\infty }{\sum }}}{\lambda }^{m-1}{F}_{j,-m},j=1,2,3$

as $\left|\lambda \right|<\max \left(\left|{\alpha }_1\right|,\cdots ,\left|{\alpha }_N\right|\right)$. By comparing the coefficients of the same power of $\lambda$, we have:

$F_{1,-m}=-\left(\langle A^{-m}{q}^1,p^1\rangle +\langle A^{-m}{q}^2,p^2\rangle +\langle A^{-m}{q}^3,p^3\rangle \right),m\ge 1;$

$\begin{array}{l}{F}_{2,-m}=r^{32}\left(\langle A^{-m}{q}^1,p^2\rangle +\langle A^{-m}{q}^3,p^1\rangle \right)\\ +\left(2r^{12}-r^{31}\right)\langle A^{-m}{q}^3,p^2\rangle \\ +{\displaystyle \underset{\begin{array}{c}k+l=m+1\\ k,l\ge 1\end{array}}{\sum }}{\displaystyle \underset{1\le i<j\le 3}{\sum }}\left|\begin{array}{cc}\langle A^{-k}{q}^i,p^i\rangle & \langle A^{-l}{q}^i,p^j\rangle \\ \langle A^{-k}{q}^j,p^i\rangle & \langle A^{-l}{q}^j,p^j\rangle \end{array}\right|,m\ge 1;\end{array}$

$\begin{array}{l}{F}_{3,-m}=-{\left(r^{32}\right)}^2\langle A^m{q}^3,p^2\rangle \\ +r^{32}{\displaystyle \underset{\begin{array}{c}k+l=m+1\\ k,l\ge 1\end{array}}{\sum }}\left(\left|\begin{array}{cc}\langle A^{-k}{q}^2,p^1\rangle & \langle A^{-l}{q}^2,p^2\rangle \\ \langle A^{-k}{q}^3,p^1\rangle & \langle A^{-l}{q}^3,p^2\rangle \end{array}\right|-\left|\begin{array}{cc}\langle A^{-k}{q}^1,p^2\rangle & \langle A^{-l}{q}^1,p^3\rangle \\ \langle A^{-k}{q}^3,p^2\rangle & \langle A^{-l}{q}^3,p^3\rangle \end{array}\right|\right)\\ -\left(2r^{12}-r^{31}\right){\displaystyle \underset{\begin{array}{c}k+l=m+1\\ k,l\ge 1\end{array}}{\sum }}\left|\begin{array}{cc}\langle A^{-k}{q}^1,p^1\rangle & \langle A^{-l}{q}^1,p^2\rangle \\ \langle A^{-k}{q}^3,p^1\rangle & \langle A^{-l}{q}^3,p^2\rangle \end{array}\right|\\ -{\displaystyle \underset{\begin{array}{c}j+k+l=m+2\\ j,k,l\ge 1\end{array}}{\sum }}\left|\begin{array}{ccc}\langle A^{-j}{q}^1,p^1\rangle & \langle A^{-k}{q}^1,p^2\rangle & \langle A^{-l}{q}^1,p^3\rangle \\ \langle A^{-j}{q}^2,p^1\rangle & \langle A^{-k}{q}^2,p^2\rangle & \langle A^{-l}{q}^2,p^3\rangle \\ \langle A^{-j}{q}^3,p^1\rangle & \langle A^{-k}{q}^3,p^2\rangle & \langle A^{-l}{q}^3,p^3\rangle \end{array}\right|,m\ge 1.\end{array}$

3. The Involutivity of Integrals

The involutivity is critical to the integrability of the Hamiltonian system, which is defined as the Poisson bracket of two conserved integrals being zero. A direct calculation gives:

Proposition 1. The Hamiltonian system for the $F_{\lambda \xi }$ -flow is expressed as:

$\frac{\text{d}{q}_k}{\text{d}{t}_{\lambda \xi }}=W_{\lambda \xi \mu }{q}_k,\frac{\text{d}{p}_k}{\text{d}{t}_{\lambda \xi }}=-W_{\lambda \xi \mu }^{\text{T}}{p}_k,$ (4)

where

$W_{\lambda \xi \mu }=\frac{1}{\lambda -{\alpha }_k}{L}_{\lambda \xi }^{*}+L_{0\lambda \xi },L_{0\lambda \xi }=\left(\begin{array}{ccc}0& 0& L_{23}\\ -2L_{23}& 0& L_{13}-L_{21}\\ 0& 0& 0\end{array}\right).$

Proof. From the Equation (3) and the property of Hamiltonian system, we calculate the partial derivatives of $q_k$ and $p_k$ with respect to $t_{\lambda \xi }$, then the results are obtained. $\square$

Lemma 2. ${\epsilon }_k=q_k{p}_k^{\text{T}}$ satisfies the Lax equation along the $t_{\lambda \xi }$ -flow:

$\frac{\text{d}{\epsilon }_k}{\text{d}{t}_{\lambda \xi }}=\left[W_{\lambda \xi \mu },{\epsilon }_k\right]=\frac{1}{\lambda -{\alpha }_k}\left[L_{\lambda \xi }^{*},{\epsilon }_k\right].$

Proof.

${\stackrel{\dot{}}{\epsilon }}_k={\stackrel{\dot{}}{q}}_k{p}_k^{\text{T}}+q_k{\stackrel{\dot{}}{p}}_k^{\text{T}}=\left(W{q}_k\right)p_k^{\text{T}}+q_k{\left(-W^{\text{T}}{q}_k\right)}^{\text{T}}=W{\epsilon }_k-{\epsilon }_{k}W=\left[W,{\epsilon }_k\right].$ $\square$

Proposition 3. $L_{\mu \eta }$ satisfies the Lax equation along the $t_{\lambda \xi }$ -flow:

$\frac{\text{d}{L}_{\mu \eta }}{\text{d}{t}_{\lambda \xi }}=\left[W_{\lambda \xi \mu }\mathrm{,}{L}_{\mu \eta }\right]\mathrm{.}$ (5)

Proof. $Q_{\mu }$ can be expressed by ${\epsilon }_k$ as:

$Q_{\mu }={\displaystyle \underset{k=1}{\overset{N}{\sum }}}\frac{{\epsilon }_k}{\mu -{\alpha }_k}.$

Resorting to Lemma 2, we have:

$\frac{\text{d}{Q}_{\mu }}{\text{d}{t}_{\lambda \xi }}=\left[W_{\lambda \xi \mu },L_{\mu \eta }\right]-\left[L_{0\lambda \xi },l_{\mu }\right],$

$\frac{\text{d}{L}_{\mu \eta }}{\text{d}{t}_{\lambda \xi }}-\left[W_{\lambda \xi \mu },L_{\mu \eta }\right]=\frac{\text{d}{l}_{\mu }}{\text{d}{t}_{\lambda \xi }}-\left[L_{0\lambda \xi },l_{\mu }\right]=0.$ $\square$

Lemma 4. The two determinants are true as $a,b,c,d=1,2,3$ :

$\left|\begin{array}{cc}{Q}_{\lambda }\left(q^a,p^b\right)& Q_{\lambda }\left(q^a,p^d\right)\\ Q_{\lambda }\left(q^c,p^b\right)& Q_{\lambda }\left(q^c,p^d\right)\end{array}\right|=\frac{1}{2}{\displaystyle \underset{k=1}{\overset{N}{\sum }}}{\displaystyle \underset{l=1}{\overset{N}{\sum }}}\frac{1}{\left(\lambda -{\alpha }_k\right)\left(\lambda -{\alpha }_l\right)}\left|\begin{array}{cc}{q}_k^a& q_l^a\\ q_k^c& q_l^c\end{array}\right|\left|\begin{array}{cc}{p}_k^b& p_k^d\\ p_l^b& p_l^d\end{array}\right|.$ (6)

$\begin{array}{l}det{\left(Q_{\lambda }\left(q^i,p^j\right)\right)}_{3\times 3}\\ =\frac{1}{6}{\displaystyle \underset{j=1}{\overset{N}{\sum }}}{\displaystyle \underset{k=1}{\overset{N}{\sum }}}{\displaystyle \underset{l=1}{\overset{N}{\sum }}}\frac{1}{\left(\lambda -{\alpha }_j\right)\left(\lambda -{\alpha }_k\right)\left(\lambda -{\alpha }_l\right)}\left|\begin{array}{ccc}{q}_j^1& q_k^1& q_l^1\\ q_j^2& q_k^2& q_l^2\\ q_j^3& q_k^3& q_l^3\end{array}\right|\left|\begin{array}{ccc}{p}_j^1& p_k^2& p_l^3\\ p_j^1& p_k^2& p_l^3\\ p_j^1& p_k^2& p_l^3\end{array}\right|.\end{array}$ (7)

From (6) and (7), we have:

$F_1\left(\lambda \right)=F_{10}\frac{{\beta }_1\left(\lambda \right)}{\alpha \left(\lambda \right)}=\left(r^{11}+r^{22}+r^{33}\right)\frac{{\beta }_1\left(\lambda \right)}{\alpha \left(\lambda \right)}={\displaystyle \underset{k=1}{\overset{N}{\sum }}}\frac{E_{1k}}{\lambda -{\alpha }_k},$ (8)

$F_2\left(\lambda \right)=F_{20}\frac{{\beta }_2\left(\lambda \right)}{\alpha \left(\lambda \right)}=-3r^{12}{r}^{32}\frac{{\beta }_2\left(\lambda \right)}{\alpha \left(\lambda \right)}={\displaystyle \underset{k=1}{\overset{N}{\sum }}}\frac{E_{2k}}{\lambda -{\alpha }_k},$ (9)

$F_3\left(\lambda \right)=F_{30}\frac{{\beta }_3\left(\lambda \right)}{\alpha \left(\lambda \right)}={\left(r^{32}\right)}^3\frac{{\beta }_3\left(\lambda \right)}{\alpha \left(\lambda \right)}={\displaystyle \underset{k=1}{\overset{N}{\sum }}}\frac{E_{3k}}{\lambda -{\alpha }_k},$ (10)

where $\alpha \left(\lambda \right)={\displaystyle {\prod }_{k=1}^N}\left(\lambda -{\alpha }_k\right)$, and ${\beta }_1\left(\lambda \right)\mathrm{,}{\beta }_2\left(\lambda \right)\mathrm{,}{\beta }_3\left(\lambda \right)$ are polynomials of degree $N-1$. From (8)-(10), another group of conserved integrals $E_{1k}\mathrm{,}{E}_{2k}\mathrm{,}{E}_{3k}$ is obtained:

$E_{1k}={\displaystyle \underset{i=1}{\overset{3}{\sum }}}{q}_k^i{p}_k^i,1\le k\le N;$

$\begin{array}{l}{E}_{2k}=-r^{32}\left(q_k^1{p}_k^2+q_k^3{p}_k^1\right)-\left(2r^{12}-r^{31}\right)q_k^3{p}_k^2\\ +{\displaystyle \underset{\begin{array}{c}1\le l\le N\\ l\ne k\end{array}}{\sum }}\frac{1}{{\alpha }_k-{\alpha }_l}{\displaystyle \underset{1\le i\mathrm{<mo>\; <\; </mo>}j\le 3}{\sum }}\left|\begin{array}{cc}{q}_k^i& q_l^i\\ q_k^j& q_l^j\end{array}\right|\left|\begin{array}{cc}{p}_k^i& p_k^j\\ p_l^i& p_l^j\end{array}\right|\mathrm{,}1\le k\le N\mathrm{<mo>\; ;\; </mo>}\end{array}$

$\begin{array}{l}{E}_{3k}={\left(r^{32}\right)}^2{q}_k^3{p}_k^2+r^{32}{\displaystyle \underset{\begin{array}{c}1\le l\le N\\ l\ne k\end{array}}{\sum }}\frac{1}{{\alpha }_k-{\alpha }_l}\left(\left|\begin{array}{cc}{q}_k^2& q_l^2\\ q_k^3& q_l^3\end{array}\right|\left|\begin{array}{cc}{p}_k^1& p_k^2\\ p_l^1& p_l^2\end{array}\right|-\left|\begin{array}{cc}{q}_k^1& q_l^1\\ q_k^3& q_l^3\end{array}\right|\left|\begin{array}{cc}{p}_k^2& p_k^3\\ p_l^2& p_l^3\end{array}\right|\right)\\ -\left(2r^{12}-r^{31}\right){\displaystyle \underset{\begin{array}{c}1\le l\le N\\ l\ne k\end{array}}{\sum }}\frac{1}{{\alpha }_k-{\alpha }_l}\left|\begin{array}{cc}{q}_k^1& q_l^1\\ q_k^3& q_l^3\end{array}\right|\left|\begin{array}{cc}{p}_k^1& p_k^2\\ p_l^1& p_l^2\end{array}\right|\\ +{\displaystyle \underset{\begin{array}{c}1\le j\le N\\ j\ne k\end{array}}{\sum }}{\displaystyle \underset{\begin{array}{c}1\le l\le N\\ l\ne j,k\end{array}}{\sum }}\frac{1}{2\left({\alpha }_j-{\alpha }_k\right)\left({\alpha }_l-{\alpha }_k\right)}\left|\begin{array}{ccc}{q}_j^1& q_k^1& q_l^1\\ q_j^2& q_k^2& q_l^2\\ q_j^3& q_k^3& q_l^3\end{array}\right|\left|\begin{array}{ccc}{p}_j^1& p_k^2& p_l^3\\ p_j^1& p_k^2& p_l^3\\ p_j^1& p_k^2& p_l^3\end{array}\right|,1\le k\le N.\end{array}$

From the definition of involutivity of two conserved integrals and direct calcultion, we have:

Proposition 5. The integrals

$\left\{F_{1m}\mathrm{,}{F}_{2m}\mathrm{,}{F}_{3m}\mathrm{,}m\in Z\mathrm{<mo>\; ;\; </mo>}{E}_{1k}\mathrm{,}{E}_{2k}\mathrm{,}{E}_{3k}\mathrm{,1}\le k\le N\right\}$

are involutive in pairs:

$\left(F_{im},F_{jn}\right)=0,\forall i,j=1,2,3;\forall m,n=0,1,2,\cdots .$ (11)

$\left(F_{i,-m},F_{j,-n}\right)=0,\forall i,j=1,2,3;\forall m,n=,1,2,\cdots .$ (12)

$\left(F_{im},F_{j,-n}\right)=0,\forall i,j=1,2,3;\forall m=0,1,2,\cdots ;n=1,2,\cdots .$ (13)

$\left(E_{ik},E_{jl}\right)=0,\forall i,j=1,2,3;\forall 1\le k,l\le N.$ (14)

$\left(F_{im},E_{jl}\right)=0,\forall i,j=1,2,3;\forall m\in Z;1\le l\le N.$ (15)

4. Hamiltonian Systems

By direct calculations, the canonical equations of the Hamiltonian systems $\left(F_{31}\right)$ and $\left(F_{21}\right)$ can be expressed as:

$\begin{array}{l}{q}_x=\frac{\partial F_{31}}{\partial p}=Uq,p_x=-\frac{\partial F_{31}}{\partial q}=-U^{\text{T}}p;\\ q_t=\frac{\partial F_{21}}{\partial p}=Vq,p_t=-\frac{\partial F_{21}}{\partial q}=-V^{\text{T}}p;\end{array}$ (16)

where $U\mathrm{,}V$ are $3N\times 3N$ matrices:

$U=\left(\begin{array}{ccc}{r}^{31}& 1& r^{11}-r^{22}\\ -\left[2r^{11}+r^{33}+{\left(r^{31}\right)}^2\right]& -r^{31}& A+2A^{32}+2r^{13}+2r^{21}+\left(r^{11}-r^{22}\right)r^{31}\\ 1& 0& 0\end{array}\right);$

$V=\left(\begin{array}{ccc}{r}^{22}+r^{33}& 0& -A+A^{32}-r^{13}\\ -\left(A+2A^{32}+r^{21}\right)& r^{11}+r^{33}& r^{31}A-A^{12}-A^{31}-r^{23}\\ -r^{31}& -1& r^{11}+r^{22}\end{array}\right)\mathrm{.}$

Proposition 6. Let $q\left(x\mathrm{,}t\right)\mathrm{,}p\left(x\mathrm{,}t\right)$ be compatible solution of the $\left(F_{31}\right)$ and $\left(F_{21}\right)$ flow on the level set $\left\{\left(p,q\right)\in {ℝ}^{6N}:r^{32}=1,r^{12}=0\right\}$. Then the Boussinesq Equation (1) has a finite-parameter solution $\left(u\mathrm{,}v\right)=B\left(q\mathrm{,}p\right)$ given as:

$u=-\frac{3}{2}{A}^{32}+\frac{1}{2}{F}_{31},v=\frac{3}{2}{r}^{11}-\frac{1}{2}{F}_{10}.$ (17)

Proof. A direct calculation gives:

$\begin{array}{l}{r}_x^{11}=-r^{13}+r^{21}+\left(r^{11}-r^{22}\right)r^{31},\\ r_x^{22}=3A^{32}+2r^{13}+r^{21}+\left(r^{11}-r^{22}\right)r^{31},\\ r_x^{33}=-3A^{32}-r^{13}-2r^{21}-2\left(r^{11}-r^{22}\right)r^{31},\\ r_x^{13}=-A^{12}+r^{23}+r^{13}{r}^{31}+\left(r^{11}-r^{22}\right)\left(r^{33}-r^{11}\right),\\ r_x^{21}=A^{31}+2r^{31}{A}^{32}-r^{23}+2r^{13}{r}^{31}-\left(r^{11}-r^{22}\right)\left(2r^{11}+r^{33}\right),\end{array}$

$\begin{array}{l}{r}_x^{23}=A^{33}-A^{22}+2\left(r^{33}-r^{22}\right)A^{32}-r^{23}{r}^{31}+\left(-2r^{11}-2r^{22}+r^{33}\right)r^{13},\\ +\left(-r^{11}-r^{22}+2r^{33}\right)r^{21}-r^{13}{\left(r^{31}\right)}^2+\left(r^{11}-r^{22}\right)\left(r^{33}-r^{22}\right)r^{31},\end{array}$

$\begin{array}{l}{r}_x^{31}=3r^{11},A_x^{12}=-A^{11}+A^{22}+2r^{31}{A}^{21}+\left(r^{11}-r^{22}\right)A^{32},\\ A_x^{31}=A^{11}-A^{33}-r^{31}{A}^{31}+\left[2r^{11}+r^{33}+{\left(r^{31}\right)}^2\right]A^{32},\\ A_x^{32}=A^{12}-A^{31}+r^{31}{A}^{32}.\end{array}$

So from the above calculations and (16), we have:

$v_t=-u_x,u_t=\frac{1}{3}{\left(v_{xx}-4v^2\right)}_x,$

and

$v_{tt}={\left(-u_x\right)}_t=-{\left(u_t\right)}_x=-\frac{1}{3}{\left(v_{xx}-4v^2\right)}_{xx}$

which is exactly the “good” Boussinesq Equation (1). $\square$

5. Conclusion

The third order differential operator associated with the Boussinesq Neumann system is a critical point in researching problems of integrable system. In this paper, we obtain a Neumann system of the Boussinesq equation associated with the third order differential operator, which is the extension of the famous KdV Neumann system associated with the second order differential operator. By means of generating function from the Lax-Moser matrix, we prove the involutivity of integrals successfully. Meanwhile, a finite-parameter solution to the Boussinesq equation is obtained naturally.

Acknowledgements

This work is supported by the Foundation (Grant No. 11601123 and 202010463050).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Tu, G.Z. (1981) The Backlund Transformation and Conserved Law for the Boussinesq Equation. Acta Mathematicae Applicatae Sinica, 1, 63-68.
[2]	Deift, P., Tomei, C. and Trubowitz, E. (1982) Inverse Scattering and the Boussinesq Equation. Communications on Pure and Applied Mathematics, 35, 567-628. https://doi.org/10.1002/cpa.3160350502
[3]	Ortega, T. and Sanz-Serna, J.M. (1990) Nonlinear Stability and Convergence of Finite-Difference Methods for the “Good” Boussinesq Equation. Numerische Mathematik, 58, 215-229. https://doi.org/10.1007/BF01385620
[4]	Schaffer, H.A., Madsen, P.A. and Deigaard, R. (1993) A Boussinesq Model for Waves Breaking in Shallow Water. Coastal Engineering, 20, 185-202. https://doi.org/10.1016/0378-3839(93)90001-O
[5]	Hu, X.B., Li, Y. and Liu, Q.M. (1993) Nonlinear Superposition Formula of the Boussinesq Hierarchy. Acta Mathematicae Applicatae Sinica, 1, 17-27.
[6]	Ma, W.X. (1993) On the Complete Integrability of Nonlinearized Lax Systems for the Classical Boussinesq Hierarchy. Acta Mathematicae Applicatae Sinica, 1, 92-96.
[7]	Flaschka, H. (1983) Relations between Infinite-Dimensional and Finite-Dimensional Isospectral Equations. Proceedings of RIMS Symposium on Nonlinear Integrable Systems, Kyoto, Japan, World Sci. Pub., Singapore, 219-239.
[8]	Zeng, Y.B. (1995) The Variable Separation of the Modified Boussinesq. Annual Mathematics, 1, 50-59.
[9]	Conte, R., Musette, M. and Pickering, A. (1995) The Two-Singular Manifold Method: II. Classical Boussinesq System. Journal of Physics A: Mathematical and General, 28, 179-187. https://doi.org/10.1088/0305-4470/28/1/020
[10]	Dai, H.H. and Geng, X.G. (2003) Finite-Dimensional Integrable Systems through the Decomposition of a Modified Boussinesq Equation. Physics Letters A, 317, 389-400. https://doi.org/10.1016/j.physleta.2003.08.049
[11]	Dougalis, V.A., Durn, A., Lpez-Marcos, M.A. and Mitsotakis, D.E. (2007) A Numerical Study of the Stability of Solitary Waves of the Bona-Smith Family of Boussinesq Systems. Journal of Nonlinear Science, 17, 569-607. https://doi.org/10.1007/s00332-007-9004-8
[12]	Lax, P.D. (1968) Integrals of Nonlinear Equations of Evolution and Solitary Waves. Communications on Pure and Applied Mathematics, 21, 467-490. https://doi.org/10.1002/cpa.3160210503
[13]	Ma, W.X. (1990) The Confocal Involutive Systems and the Integrability of the Nonlinearized Lax Systems of AKNS Hierarchy. In: Gu, C., Li, Y., Tu, G. and Zeng, Y., Eds., Nonlinear Physics, Research Reports in Physics, Springer-Verlag, Berlin, 79-84. https://doi.org/10.1007/978-3-642-84148-4_10
[14]	Wang, X., Cao, J.L. and Chen, Y. (2015) Higher-Order Rogue Wave Solutions of the Three-Wave Resonant Interaction Equation via the Generalized Darboux Transformation. Physica Scripta, 90, 10. https://doi.org/10.1088/0031-8949/90/10/105201
[15]	Cao, C.W. (1990) Nonlinearization of the Lax System for AKNS Hierarchy. Science China A, 33, 528-536.
[16]	Cao, C.W. and Cao, J.L. (2006) On an Integrable (2+1) Dimensional Korteweg-de Vries Equation with a Discrete Variable. Il Nuovo Cimento B, 121, 675-687.
[17]	Cao, C.W. and Cao, J.L. (2007) An Integrable (2+1)-Dimensional Toda Equation with Two Discrete Variables. Physics Letters A, 365, 301-308. https://doi.org/10.1016/j.physleta.2007.01.028
[18]	Cao, J.L. (2007) Finite Parameter Solutions to the Manakov and the Derivative Manakov Equations. Chinese Quarterly Journal of Mathematics, 22, 236-244.
[19]	Cao, J.L., Zhang, H. and Jiao, W.T. (2008) Coupled KdV Equations and Their Explicit Solutions Through Two-Dimensional Hamiltonian System with a Quartic Potential. Communications in Theoretical Physics, 49, 1379-1382. https://doi.org/10.1088/0253-6102/49/6/04
[20]	Cao, J.L. and Shi, W. (2013) A New Integrable Generalized Neumann System Related to Soliton Hierarchy. Chaos, Solitons and Fractals, 47, 13-17. https://doi.org/10.1016/j.chaos.2012.11.005