Abstract
The Painlevé integrability of the higher-order Boussinesq equation is proved by using the standard Weiss-Tabor-Carnevale (WTC) method. The multisoliton solutions of the higher-order Boussinesq equation are obtained by introducing dependent variable transformation. The soliton molecule and asymmetric soliton of the higher-order Boussinesq equation can be constructed by the velocity resonance mechanism. Lump solution can be derived by solving the bilinear form of the higher-order Boussinesq equation. By some detailed calculations, the lump wave of the higher-order Boussinesq equation is just the bright form. These types of the localized excitations are exhibited by selecting suitable parameters.
1. Introduction
The soliton molecule as a boundary state is comprised by a balance of repulsive and attractive forces between solitons caused by the nonlinear and dispersive effects [1]. The characteristics of the soliton molecule in both experiment and simulation have attracted considerable attention [2–4]. The soliton molecules are first discovered by theoretical analysis of the nonlinear Schröinger equation [5] and the complex Ginzburg-Landau equation [6]. The theoretical frameworks to address the soliton molecule have been proposed [7, 8]. Recently, Lou proposed the velocity resonance mechanism to form the soliton molecule [9]. The high-order dispersive terms play a key role in the velocity resonance mechanism [10]. The velocity resonance mechanism is developed to some integrable systems, the (2+1)-dimensional fifth-order Korteweg-de Vries (KdV) equation [11], the complex modified KdV equation [12], the (3+1)-dimensional Boiti-Leon-Manna-Pempinelli equation [13], and so on [14–16]. Combining the Darboux transformation and the variable separation approach, some interactions between soliton molecules and breather solutions and between soliton molecules and dromions are explored [11–15, 17]. In addition to the soliton molecule, lump solutions are a kind of rational function solutions which have become a hot field in nonlinear systems [18–22]. Lump solutions will decay polynomially in all directions of space [23]. The Hirota bilinear method is a useful method to find lump solutions. Lump waves of the nonlinear systems have been validated by the Hirota bilinear method [24–29]. In this paper, our objective is to explore a higher-order Boussinesq equation which is considered as the combination between the fourth-order and sixth-order Boussinesq forms. We get lots of interesting results for the higher-order Boussinesq equation, such as the multisoliton, the soliton molecule, and lump solution. Particularly, the soliton molecule of the higher-order Boussinesq equation is not valid for just the fourth-order Boussinesq equation or the sixth-order Boussinesq equation.
The higher-order Boussinesq equation reads where , , and are arbitrary constants. The higher-order Boussinesq equation (1) will become the fourth-order Boussinesq equation and the sixth-order Boussinesq equation with and , respectively. The higher-order Boussinesq equation (1) can describe the propagation of long waves in shallow water. The integrability of the fourth-order Boussinesq equation was handled by the inverse scattering transform method [30]. The group-invariant solutions of the fourth-order Boussinesq equation were obtained by similarity reductions [31]. The high-dimensional Boussinesq equations can be constructed by means of the fourth-order Boussinesq equation [32]. The -soliton solutions for the fourth-order [33] and the sixth-order Boussinesq equations [34] were obtained by the Hirota bilinear method.
This paper is organized as follows. In Section 2, the Painlevé integrability of the higher-order Boussinesq equation is proved by the standard Weiss-Tabor-Carnevale (WTC) approach. In Section 3, the higher-order Boussinesq equation can be transformed to the bilinear form by the dependent variable transformation. The Hirota’s bilinear method is employed to derive the multisoliton by handling the bilinear form of the higher-order Boussinesq equation. The soliton molecule of the higher-order Boussinesq equation is constructed by a new resonance condition. In Section 4, lump solution of the higher-order Boussinesq equation is obtained by solving the corresponding Hirota bilinear form. Finally, the conclusions of this paper follow in Section 5.
2. Painlevé Analysis for the Higher-Order Boussinesq Equation
According to the WTC approach [35], the solution of the higher-order Boussinesq equation can be written as with being a positive integer. The solution of the model is single valued about the arbitrary movable singularity manifold . From the leading order analysis, we get
By inserting the expression into (1) and vanishing the coefficient term of , the polynomial equation in is derived as
The corresponding resonances occur at
The resonance at corresponds to the fact that the singularity manifold is an arbitrary function. By vanishing subsequent coefficient terms of , we infer that the number of arbitrary functions (,) is the same as the number of resonances (). The higher-order Boussinesq equation is thus the Painlevé integrability in the sense of the WTC test.
3. Multisoliton and Soliton Molecule for the Higher-Order Boussinesq Equation
To determine the multisoliton solutions of the higher-order Boussinesq equation (1), the dependent variable transformation reads
By substituting (7) into (1), the bilinear form of the higher-order Boussinesq equation reads where is bilinear derivative operator [36].
The -soliton solutions of the higher-order Boussinesq equation can be calculated as where is the summation with possible combinations of and . By substituting (7) and (10) into (1), the phases shift and the dispersion relation read as
We shall study the soliton molecule with a new resonance condition. The new resonance condition of velocity resonance reads
By solving the above condition (13), the velocity resonant condition becomes
A soliton molecule and an asymmetric soliton can be constructed by selecting appropriate parameters in (13) or (14). We take two-soliton () in (10) and the new resonance condition to describe these phenomena. For the Figure 1, the parameters are selected as

(a)

(b)
For the Figure 2, the parameters are selected as

(a)

(b)
The soliton molecule and the asymmetric soliton are described in Figures 1 and 2, respectively. Two solitons in molecule are different amplitude, while two solitons in molecule possess the same velocity simultaneously.
4. Lump Solution of the Higher-Order Boussinesq Equation
To obtain lump solution of the higher-order Boussinesq equation (1), we take a quadratic function as
By substituting (17) into the Hirota bilinear form (8) and balancing the different powers of and , the parameters are constrained as
The function of can be localized in -plane with the parameters satisfying
By substituting (17) into (7), a lump wave of (1) is generated with . To describe the lump wave of the higher-order Boussinesq equation (1), the parameters are selected as
The spatiotemporal structure and the corresponding density of a lump wave are described in Figures 3(a) and 3(b), respectively. The critical point of the lump wave can be calculated by solving

(a)

(b)
By solving the above condition (22), we find that the function reaches the maximum value at the point and the minimum values at two points . By substituting the above three points values into (20), the maximum and minimum values of the function are and , respectively. We can only get the bright lump form of the higher-order Boussinesq equation by the above detail analysis.
5. Conclusion
In summary, the Painlevé property, the soliton molecule, and the lump solution of the higher-order Boussinesq equation (1) are studied by the standard WTC and the Hirota bilinear methods. The multisoliton solutions of (1) are obtained by introducing dependent variable transformation. The soliton molecule and the asymmetric soliton are constructed by the velocity resonance mechanism. The lump solution can be derived by using a positive quadratic function. By detail calculations of the maximum and minimum values of the function , lump wave of the higher-order Boussinesq equation is just the bright form.
In this paper, the higher-order Boussinesq equation which possesses Painlevé integrability is constructed by introducing the Hirota bilinear operator based on the fourth-order Boussinesq equation. Similar introducing high-order Hirota bilinear operator procedure, we propose one equation
with being arbitrary constants. The integrabilityproperties and nonlinear excitations of (23) are worthy to study further.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work is supported by the National Natural Science Foundation of China No. 11775146.