## General rogue waves in the Boussinesq equation

Bo Yang and Jianke Yang

Department of Mathematics and Statistics, University of Vermont, Burlington, VT 05405, USA

We derive general rogue wave solutions of arbitrary orders in the Boussinesq equation by the bilinear Kadomtsev-Petviashvili (KP) reduction method. These rogue solutions are given as Gram determinants with 2 N -2 free irreducible real parameters, where N is the order of the rogue wave. Tuning these free parameters, rogue waves of various patterns are obtained, many of which have not been seen before. Compared to rogue waves in other integrable equations, a new feature of rogue waves in the Boussinesq equation is that the rogue wave of maximum amplitude at each order is generally asymmetric in space. On the technical aspect, our contribution to the bilinear KP-reduction method for rogue waves is a new judicious choice of differential operators in the procedure, which drastically simplifies the dimension reduction calculation as well as the analytical expressions of rogue wave solutions.

## I. INTRODUCTION

In 1871, Boussinesq introduced an equation which governs the propagation of long surface waves on water of constant depth [1, 2] (see also [3]). After variable normalizations, this equation can be written as

$$\eta _ { t t } - \eta _ { x x } - ( \eta ^ { 2 } ) _ { x x } - \frac { 1 } { 3 } \eta _ { x x x x } = 0.$$

Actually, the quantities η, x and t can be further scaled to produce any desired coefficients for the terms in this equation. The present choice of coefficients is taken following [4]. This Boussinesq equation also arises in many other physical contexts, such as continuum approximations of certain Fermi-Pasta-Ulam chains [5-7], and ion sound waves in a plasma [8, 9]. Remarkably, this equation is integrable. Indeed, its multi-soliton solutions (by the bilinear method) and its Lax pair were reported almost simultaneously by Hirota and Zakharov in 1973 [5, 10] (see also [11-13]).

Rogue waves are spontaneous large nonlinear waves which 'come from nowhere and disappear with no trace' [14]. In other words, they arise from a constant (uniform) background, run to high amplitudes, and then retreat back to the same background. These waves are associated with freak waves in the ocean and extreme events in optical fibers [15-18]. Thus, they have attracted a lot of attention in the physics and nonlinear waves communities in recent years. Analytical expressions of rogue waves have been derived in a large number of integrable nonlinear wave equations, such as the nonlinear Schr¨dinger (NLS) equation [19-28], the derivative NLS equation [29, 30], the three-wave interaction o equation [31], the Davey-Stewartson equations [32, 33], and many others [34-53]. Rogue waves have also been observed in water tanks [54, 55] and optical fibers [56-58].

In this paper, we consider rogue waves in the Boussinesq equation (1). Since these waves arise from and retreat back to a constant background, we let

$$\eta ( x, t ) \to \eta _ { 0 }, \ \ x, t \to \pm \infty,$$

where η 0 is a constant. For rogue waves to appear, the background solution η = η 0 must be unstable to long-wave perturbations [41]. Simple calculations show that this so-called baseband instability occurs when 1 + 2 η 0 ≤ 0. Under this condition, after a variable shift of η = η 0 + u and proper scalings of x, t and u , the Boussinesq equation (1) reduces to

$$u _ { t t } + u _ { x x } - ( u ^ { 2 } ) _ { x x } - \frac { 1 } { 3 } u _ { x x x x } = 0,$$

and boundary conditions (2) reduce to

$$u ( x, t ) \to 0, \ \ x, t \to \pm \infty.$$

Fundamental (first-order) rogue waves to the Boussinesq equation (3) were derived in [59, 60] by taking a longwave limit of the two-soliton solution, following a similar procedure introduced in [61] where several singular rational solutions for the Boussinesq equation (1) were presented. Higher-order rogue waves to the Boussinesq equation (3) were recently considered by Clarkson and Dowie [4]. By converting this equation into a bilinear one, assuming certain polynomial forms for the bilinear solution, equating powers of x, t in the bilinear equation and solving the resulting algebraic equations for the polynomial coefficients, the authors obtained second- to fifth-order rogue waves with two free real parameters at each order.

Despite this progress, many important questions on rogue waves of the Boussinesq equation are still open. One of the most significant questions is that analytical expressions for rogue waves of arbitrary orders are still unknown. A related question is how many free real parameters exist in general rogue waves of arbitrary orders. Recall that all higher-order rogue waves reported in [4] contained the same number of free real parameters (which is two). But from experience of rogue waves in other integrable equations (such as the NLS equation [21, 25, 26]), one anticipates that the number of free parameters should increase with the order of the rogue wave. Thus, general rogue waves of high order in the Boussinesq equation should contain more parameters than two, and those solutions with more free parameters are still unclear. A third open question is what is the maximum amplitude that can be attained in rogue waves of a given order, and what is the spatial-temporal profile of that rogue wave of maximum amplitude. Since the Boussinesq equation is a physically and mathematically interesting equation, these open questions on its rogue waves clearly merit thorough investigation.

In this article, we derive general rogue waves of arbitrary orders in the Boussinesq equation (3) under boundary conditions (4). The technique we will use is the bilinear Kadomtsev-Petviashvili (KP) reduction method [62]. Although this technique has been applied to derive rogue waves in several other integrable equations before [26, 32, 33, 42, 4951, 53], when it is applied to the Boussinesq equation, the previous choices of differential operators within this technique would cause considerable technical complications. Here, we propose a new judicious choice of those differential operators, which will drastically simplify the dimension reduction calculation as well as the analytical expressions of rogue wave solutions. Our rogue solutions are given as Gram determinants with 2 N -2 free irreducible real parameters, where N is the order of the rogue wave. Thus, rogue waves of higher orders do contain more free parameters, and the rogue waves reported in [4] are special cases. Tuning these free parameters, we obtain rogue waves of various interesting patterns, many of which have not been seen before. We also find that in the Boussinesq equation, the rogue wave of maximum amplitude at each order is generally asymmetric in space. This is a surprising feature which contrasts rogue waves in other integrable equations.

## II. GENERAL ROGUE-WAVE SOLUTIONS

Our general rogue waves in the Boussinesq equation (3) are given by the following two theorems.

Theorem 1. The Boussinesq equation (3) under boundary conditions (4) admits rational nonsingular rogue-wave solutions

$$u ( x, t ) = 2 \partial _ { x } ^ { 2 } \ln \sigma,$$

where

$$\sigma ( x, t ) = \det _ { 1 \leq i, j \leq N } \left ( m _ { 2 i - 1, 2 j - 1 } \right ),$$

the matrix elements in σ are defined by

$$m _ { i, j } = \sum _ { k = 0 } ^ { i } \sum _ { l = 0 } ^ { j } \frac { a _ { k } } { ( i - k )! } \frac { ( - 1 ) ^ { l } a _ { l } ^ { * } } { ( j - l )! } \left [ f ( p ) \partial _ { p } \right ] ^ { i - k } \left [ f ( q ) \partial _ { q } \right ] ^ { j - l } \left ( \frac { 1 } { p + q } e ^ { \frac { 1 } { k } ( p + q ) z - \frac { 1 } { k } ( p ^ { 2 } - q ^ { 2 } ) i t } \right ) \Big | _ { p = - 1, \ q = - 1 },$$

with

$$f ( p ) = \frac { \sqrt { p ^ { 2 } - 4 } } { 3 }, \ \ f ( q ) = \frac { \sqrt { q ^ { 2 } - 4 } } { 3 },$$

the asterisk '*' represents complex conjugation, and a k are free complex constants.

The matrix elements in the above theorem are expressed in terms of derivatives with respect to dummy parameters p and q . Purely algebraic expressions for these matrix elements can also be obtained through elementary Schur polynomials. For this purpose, we first introduce elementary Schur polynomials S k ( x ) which are defined via the generating function,

$$\sum _ { k = 0 } ^ { \infty } S _ { k } ( x ) \lambda ^ { k } = \exp \left ( \sum _ { j = 1 } ^ { \infty } x _ { j } \lambda ^ { j } \right ),$$

where x = ( x , x 1 2 , · · · ). Specifically, we have

$$S _ { 0 } ( x ) = 1, \ \ S _ { 1 } ( x ) = x _ { 1 }, \ \ S _ { 2 } ( x ) = \frac { 1 } { 2 } x _ { 1 } ^ { 2 } + x _ { 2 }, \ \cdots, \ \ S _ { k } ( x ) = \sum _ { l _ { 1 } + 2 l _ { 2 } + \cdots + m l _ { m } = k } \left ( \prod _ { j = 1 } ^ { m } \frac { x _ { j } ^ { l _ { j } } } { l _ { j } ^ { l _ { j } } } \right ). \quad \quad ( 1 0 )$$

In terms of these Schur polynomials, matrix elements m i,j in Eq. (7) can be replaced by formulae in the following theorem.

Theorem 2. Purely algebraic expressions of matrix elements m i,j in Eq. (7) can be written as

$$m _ { i, j } & = \sum _ { \nu = 0 } ^ { \min ( i, j ) } \Phi _ { i \nu } \Psi _ { j \nu },$$

where

$$\Phi _ { i \nu } & = \left ( \frac { i } { 2 \sqrt { 3 } } \right ) ^ { \nu } \sum _ { k = 0 } ^ { i - \nu } a _ { k } S _ { i - \nu - k } ( x ^ { + } + \nu s ),$$

$$\Psi _ { j \nu } & = \left ( \frac { i } { 2 \sqrt { 3 } } \right ) ^ { \nu } \sum _ { l = 0 } ^ { j - \nu } ( - 1 ) ^ { l } a _ { l } ^ { * } S _ { j - \nu - l } ( x ^ { - } + \nu s ),$$

and vectors x ± = ( x ± 1 , x ± 2 , · · · ) , r = ( r , r 1 2 , · · · ), s = ( s , s 1 2 , · · · ) are defined by

$$x _ { k } ^ { \pm } & = \frac { e ^ { 2 i \pi / 3 } + ( - 1 ) ^ { k } e ^ { - 2 i \pi / 3 } } { 2 \cdot 3 ^ { k } \cdot k! } \left [ x \pm ( - 2 ) ^ { k - 1 } i t \right ] + r _ { k }, \\ \sim & \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd.$$

$$\sum _ { k = 1 } ^ { \infty } r _ { k } \lambda ^ { k } = - \ln _ { \substack { \lambda \\ \rho } } \left [ \frac { 1 } { 2 } - \cosh \left ( \frac { \lambda } { 3 } + \frac { 2 i \pi } { 3 } \right ) \right ], \quad \ \quad \ \.$$

$$\sum _ { k = 1 } ^ { k = 1 } s _ { k } \lambda ^ { k } = \ln \left [ \frac { 2 i \sqrt { 3 } } { \lambda } \tanh \frac { \lambda } { 6 } \tanh \left ( \frac { \lambda } { 6 } + \frac { 2 i \pi } { 3 } \right ) \right ].$$

The first few ( r k , s k ) values are

$$r _ { 1 } = \frac { i } { 2 \sqrt { 3 } }, \quad r _ { 2 } = - \frac { 5 } { 7 2 }, \quad r _ { 3 } = - \frac { i } { 5 4 \sqrt { 3 } }, \quad s _ { 1 } = \frac { 2 i } { 3 \sqrt { 3 } }, \quad s _ { 2 } = - \frac { 5 } { 1 0 8 }, \quad s _ { 3 } = - \frac { 5 i } { 2 4 3 \sqrt { 3 } }.$$

Remark 1. The two expressions (7) and (11) for the matrix element m i,j in Theorems 1 and 2 differ by a constant factor multiplying an exponential of a linear function in x and t , but they give the same solution for u x, t ( ).

Remark 2. Even though the matrix elements in these theorems contain complex numbers, we will show that the determinant (6) and the resulting solution u x, t ( ) are real (see Sec. IVC).

Remark 3. The rogue waves in Theorem 1 contain 2 N free complex parameters a , a 0 1 , . . . , a 2 N -1 , but these free parameters are reducible. First, we can set a 0 = 1 without loss of generality. Second, through a determinant manipulation similar to that in [26], we can set

$$a _ { 2 } = a _ { 4 } = a _ { 6 } = \cdots = a _ { e v e n } = 0$$

without loss of generality. These a 0 and a even values will be taken throughout the rest of this article. Thirdly, with a shifting of x and t , we can also normalize a 1 = 0 (as was done in [26]). Thus, the N -th order rogue waves in the Boussinesq equation contain N -1 free irreducible complex parameters, or 2 N -2 free irreducible real parameters. This number is the same as that for rogue waves in the NLS equation [26]. Recalling that the second- and higherorder Boussinesq rogue waves reported in [4] contain only two free real parameters, this means that those solutions are special cases of ours in the above theorems. As an example, our third-order rogue waves for the Boussinesq equation contain 4 free irreducible real parameters, which doubles the number of free real parameters in third-order rogue waves of Ref. [4].

Remark 4. Using the algebraic expressions of matrix elements in Theorem 2, we can show that when a , a 1 3 , . . . , a odd are all purely imaginary (and a , a 0 2 , . . . a even are taken as in Remark 3), then the rogue solution (5) is symmetric in time, i.e.,

$$u ( x, t ) = u ( x, - t ).$$

The reason is that, if t →-t , then x + k → x -k . In addition, if all a odd are purely imaginary and all a even purely real, then a ∗ k = ( -1) k a k . In this case, we see from Eqs. (11)-(13) that m ij ( x, -t ) = m ji ( x, t ), and thus u x, t ( ) = u x, ( -t ).

## III. ANALYSIS OF GENERAL ROGUE WAVES

## A. First- and second-order rogue waves

Setting N = 1 and a 1 = 0, we get the first-order rogue wave

$$u _ { 1 } ( x, t ) = 2 \partial _ { x } ^ { 2 } \ln \left ( \hat { x } ^ { 2 } + t ^ { 2 } + 1 \right ) = \frac { 4 ( 1 - \hat { x } ^ { 2 } + t ^ { 2 } ) } { 1 + \hat { x } ^ { 2 } + t ^ { 2 } },$$

where ˆ = x x +1. This solution matches the one derived earlier in [4, 59, 60], and it does not contain any irreducible free parameters.

For second-order rogue waves, we set N = 2,

$$a _ { 1 } = 0, \ \ a _ { 3 } = \frac { s _ { 3 } + i r _ { 3 } } { 7 2 \sqrt { 3 } },$$

where s 3 and r 3 are free real constants. Then, the corresponding rogue waves in Eq. (5) become

$$u _ { 2 } ( x, t ) = 2 \partial _ { x } ^ { 2 } \ln \sigma _ { 2 },$$

where

$$\sigma _ { 2 } ( x, t ) & = x ^ { 6 } + t ^ { 6 } + 3 x ^ { 4 } t ^ { 2 } + 3 x ^ { 2 } t ^ { 4 } + 1 4 x ^ { 5 } + 1 4 x t ^ { 4 } + 2 8 x ^ { 3 } t ^ { 2 } + 9 0 x ^ { 4 } + 1 2 8 x ^ { 2 } t ^ { 2 } + 2 2 t ^ { 4 } + \\ & \quad 3 2 4 x ^ { 3 } + 3 1 6 x t ^ { 2 } + 6 4 8 x ^ { 2 } + 3 6 0 t ^ { 2 } + 6 4 8 x + 3 2 4 + \\ & \quad 2 r _ { 3 } \left ( x ^ { 3 } - 3 x t ^ { 2 } + 7 x ^ { 2 } - 7 t ^ { 2 } + 1 6 x + 8 \right ) + 2 s _ { 3 } t \left ( t ^ { 2 } - 3 x ^ { 2 } - 1 4 x - 1 8 \right ) + r _ { 3 } ^ { 2 } + s _ { 3 } ^ { 2 }.$$

This solution contains two irreducible free real parameters s 3 and r 3 . Recall that the second-order rogue waves derived in [4] also contain two free real parameters. It turns out that the second-order rogue waves in [4] and our solution (22) above are equivalent. Indeed, our function σ 2 in (23) can be rewritten as

$$\sigma _ { 2 } \left ( x, t \right ) = \sigma _ { 2 } ^ { ( s ) } ( \hat { x }, t ) + 2 \alpha t \left ( 3 \hat { x } ^ { 2 } - t ^ { 2 } + \frac { 5 } { 3 } \right ) + 2 \beta \hat { x } \left ( \hat { x } ^ { 2 } - 3 t ^ { 2 } - \frac { 1 } { 3 } \right ) + \alpha ^ { 2 } + \beta ^ { 2 },$$

where ˆ = x x + 7 3 , α = -s 3 , β = r 3 -106 27 , and

$$\sigma _ { 2 } ^ { ( s ) } ( x, t ) = x ^ { 6 } + \left ( 3 t ^ { 2 } + \frac { 2 5 } { 3 } \right ) x ^ { 4 } + \left ( 3 t ^ { 4 } + 3 0 t ^ { 2 } - \frac { 1 2 5 } { 9 } \right ) x ^ { 2 } + t ^ { 6 } + \frac { 1 7 } { 3 } t ^ { 4 } + \frac { 4 7 5 } { 9 } t ^ { 2 } + \frac { 6 2 5 } { 9 }$$

is a second-order ( x, t )-symmetric rogue wave in the Boussinesq equation [4]. The above form of σ 2 ( x, t ) matches the generalized second-order rogue waves derived in [4] (except for a space shift).

## B. The third-order rogue waves

To get general third-order rogue waves, we set N = 3 and a 1 = 0. The corresponding rogue waves (5) contain two irreducible free complex parameters a 3 and a 5 , or four irreducible free real parameters. This number of irreducible free parameters doubles that in the third-order rogue waves of Ref. [4]. Thus, our third-order rogue waves contain many new solutions. For example, two triangular patterns with different orientations as well as two mixed patterns, together with their corresponding a 3 and a 5 values, are displayed in Fig. 1. When compared to the third-order rogue patterns reported in [4], we see that the present wave patterns are quite different. It is noted that some of these patterns, such as the triangular ones in the upper row of Fig. 1, look quite similar to those in the local NLS equation [26]. But the mixed patterns in the lower row of Fig. 1 have not been seen in the NLS equation.

## C. Rogue waves of maximum amplitude

Next, we explore the maximum amplitude that can be reached by rogue waves of each order. For the NLS equation, this question has been addressed [20, 22, 26, 28]. But for the Boussinesq equation, this question is still open.

FIG. 1: Third-order rogue waves in the Boussinesq equation. Free parameters in these solutions are chosen as (a) a 3 = 1, a 5 = 0; (b) a 3 = i, a 5 = 0; (c) a 3 = i / 2, a 5 = i / 2; and (d) a 3 = -i / 2, a 5 = i / 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00109v1-figure-baf368e0-61d3-4d5d-bf18-698355c2b1f9.jpg)

For the first-order rogue wave (20), since it has no irreducible free parameters, its maximum amplitude can be easily seen as 4, which is attained at ˆ = x t = 0. Shifting the location of this maximum amplitude to the origin ( x, t ) = (0 , 0), which is equivalent to choosing a 1 = -i / ( 2 √ 3 ) in the rogue waves of Theorems 1 and 2, this first-order rogue wave of maximum amplitude is plotted in Fig. 2(a). Notice that this solution is symmetric in both x and t .

For second-order rogue waves, we find that their maximum amplitude is 5 5. If we require this maximum amplitude . to be located at the origin x = t = 0 (which means that we cannot normalize a 1 to be zero), then there are two such rogue waves u (1) ( x, t ) and u (2) ( x, t ) with this maximum amplitude, and their corresponding ( a , a 1 3 ) values are

$$\cdot \cdot \cdot & & \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \\ & a _ { 1 } = \mathfrak { i } \frac { 9 - 7 \sqrt { 3 } } { 1 8 }, \quad a _ { 3 } = \mathfrak { i } \frac { 7 2 - 4 7 \sqrt { 3 } } { 3 2 4 },$$

and

$$a _ { 1 } = - \mathfrak { i } \frac { 9 + 7 \sqrt { 3 } } { 1 8 }, \quad a _ { 3 } = - \mathfrak { i } \frac { 7 2 + 4 7 \sqrt { 3 } } { 3 2 4 }.$$

The profile of the first rogue wave u (1) ( x, t ) is plotted in Fig. 2(b), and the second wave is related to the first one by u (2) ( x, t ) = u (1) ( -x, t ). Notice that these second-order rogue waves of maximum amplitude are symmetric in t , but asymmetric in x . This contrasts the first-order rogue wave, which is symmetric in both x and t [see Fig. 2(a)].

For third-order rogue waves, their maximum amplitude attained at the origin becomes approximately 6 784. . Two rogue waves feature this maximum amplitude, and they are related to each other by switching x to -x . The first wave has parameter values

$$a _ { 1 } \approx - 0. 1 1 6 3 i, \ \ a _ { 3 } \approx - 0. 0 2 3 8 i, \ \ a _ { 5 } \approx - 0. 0 0 0 2 i,$$

and its profile is plotted in Fig. 2(c). The second wave has parameter values

$$a _ { 1 } \approx - 2. 0 0 1 i, \ \ a _ { 3 } \approx - 2. 5 8 4 i, \ \ a _ { 5 } \approx - 3. 9 3 7 i.$$

For fourth-order rogue waves, their maximum amplitude attained at the origin is approximately 7 944. Again, two . rogue waves feature this maximum amplitude. The first one, with parameter values

$$a _ { 1 } \approx - 0. 0 7 9 9 i, \ \ a _ { 3 } \approx - 0. 0 2 1 3 i, \ \ a _ { 5 } \approx 0. 0 0 0 1 i, \ \ a _ { 7 } \approx - 0. 0 0 0 1 i,$$

is plotted in Fig. 2(d), and the second one has parameter values

$$a _ { 1 } \approx - 1. 0 8 0 i, \ \ a _ { 3 } \approx - 0. 3 4 8 i, \ \ a _ { 5 } \approx - 0. 1 3 3 i, \ \ a _ { 7 } \approx - 0. 0 5 1 i.$$

FIG. 2: Profiles of the first- to fourth-order rogue waves of maximum amplitude in the Boussinesq equation. The parameter values for these rogue waves are given in the text.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00109v1-figure-ac70a7a7-b298-464d-b204-f28f49d1abe7.jpg)

For fifth- and sixth-order rogue waves, their maximum amplitudes are approximately 9.017 and 10.025 respectively. Profiles and parameter values for these rogue waves of maximum amplitude are omitted.

All rogue waves of maximum amplitude from the second order up are asymmetric in x and symmetric in t . At each order, two such waves exist which are mirror images of each other around x = 0.

The above results are briefly summarized in the following table.

TABLE I: Maximum amplitudes of rogue waves at each order N and their x -symmetry

|   N | Maximum amplitude   | x -symmetry   |   Amplitude of symmetric wave |
|-----|---------------------|---------------|-------------------------------|
|   1 | 4                   | symmetric     |                         4     |
|   2 | 5 . 5               | asymmetric    |                         4.846 |
|   3 | 6 . 784             | asymmetric    |                         6.545 |
|   4 | 7 . 944             | asymmetric    |                         6.956 |
|   5 | 9 . 017             | asymmetric    |                         8.648 |
|   6 | 10 . 025            | asymmetric    |                         8.757 |

It should be noted that at each order, there does exist a rogue wave which is symmetric in both x and t . Such symmetric solutions up to order six have been derived and plotted in Ref. [4]. They are special solutions of our general rogue waves (5) as well. For instance, the first-order symmetric rogue wave is obtained from our formula (5) when we choose a 1 = -i / (2 √ 3); the second-order one is obtained when we choose a 1 = -7i / (6 √ 3) and a 3 = -31i / (216 √ 3); the third-order one is obtained when we choose a 1 = -11i / (6 √ 3), a 3 = -409i / (648 √ 3) and a 5 = -1619i / (6480 √ 3); and for the fourth-order, the parameters are a 1 = -5i / (2 √ 3), a 3 = -359i / (216 √ 3), a 5 = -5041i / (3888 √ 3), and a 7 = -29937263i / (29393280 √ 3); and so on.

However, except for order one, this symmetric rogue wave does not attain the maximum amplitude of that order. For instance, the peak amplitudes of symmetric rogue waves from the second to sixth orders are approximately 4.846, 6.545, 6.956, 8.648 and 8.757 respectively. These peak values of symmetric rogue waves are also listed in Table 1. Compared to the maximum amplitudes in that table, it is clear that these peak amplitudes of symmetric rogue waves are lower than the maximum amplitudes of x -asymmetric rogue waves (except for order one).

The fact that rogue waves of maximum amplitude in the Boussinesq equation are generally asymmetric rather than symmetric in space is very counterintuitive, considering that the Boussinesq equation itself is symmetric in space. This asymmetry also contrasts the NLS equation, where the maximum-amplitude rogue waves are always symmetric

in space [20, 22, 26, 28]. These asymmetric wave patterns, as shown in Fig. 2(b-d), are very novel and have not been seen in rogue waves of other integrable equations (to the authors' best knowledge).

## IV. DERIVATION OF ROGUE-WAVE SOLUTIONS

In this section, we derive the general rogue-wave solutions given in Theorems 1 and 2. This derivation uses the bilinear KP-reduction method in the soliton theory. This method is based on Hirota's bilinear form of an integrable equation [62], and the observation that this bilinear equation is often a member of the KP hierarchy [63] (possibly after certain reductions such as the dimension reduction). Thus, solutions of the KP hierarchy, under restrictions due to those reductions, will provide solutions to the original integrable system. This technique often gives elegant determinant-type solutions. In addition, it produces solutions to higher-dimensional integrable equations more easily than to lower-dimensional ones, which is remarkable. This bilinear KP-reduction technique has been applied to derive rogue waves in several other integrable equations before [26, 32, 33, 42, 49-51, 53]. Since rogue waves are rational solutions, the way to derive such rational terms in this bilinear KP method is to define matrix elements as certain differential operators (with respect to parameters) acting on an exponential term. When this technique is applied to the Boussinesq equation, however, the previous choices of such differential operators would cause considerable technical difficulties for the dimension reduction (see [49, 50] for examples). To overcome those difficulties, we will introduce a new and general way to choose these differential operators. As we will see, this new treatment will streamline the dimension reduction calculation and simplify the analytical expressions of rogue wave solutions.

The outline of our derivation is as follows.

First, we make the standard variable transformation

$$u ( x, t ) = 2 \partial _ { x } ^ { 2 } \ln \sigma ( x, t ),$$

where σ x, t ( ) is a real variable. Under this transformation, the Boussinesq equation (3) is converted into the bilinear form [4]

$$\left ( D _ { x } ^ { 4 } - 3 D _ { x } ^ { 2 } - 3 D _ { t } ^ { 2 } \right ) \sigma \cdot \sigma = 0,$$

where D is Hirota's bilinear differential operator defined by

$$& P \left ( D _ { x }, D _ { y }, D _ { t }, \cdots \right ) F ( x, y, t, \cdots ) \cdot G ( x, y, t, \cdots ) \\ & \equiv P \left ( \partial _ { x } - \partial _ { x ^ { \prime } }, \partial _ { y } - \partial _ { y ^ { \prime } }, \partial _ { t } - \partial _ { t ^ { \prime } }, \cdots \right ) F ( x, y, t, \cdots ) G ( x ^ { \prime }, y ^ { \prime }, t ^ { \prime }, \cdots ) | _ { x ^ { \prime } = x, y ^ { \prime } = y, t ^ { \prime } = t, \cdots },$$

and P is a polynomial of D x , D y , D , t · · · .

In order to derive solutions to the bilinear equation (33), we consider a higher-dimensional bilinear equation

$$\left ( D _ { x _ { 1 } } ^ { 4 } - 4 D _ { x _ { 1 } } D _ { x _ { 3 } } + 3 D _ { x _ { 2 } } ^ { 2 } \right ) \sigma \cdot \sigma = 0,$$

which is the bilinear form of the KP equation [62, 63]. We first construct a wide class of algebraic solutions for this higher-dimensional bilinear equation in the form of Gram determinants. Then, we restrict these solutions so that they satisfy the dimension-reduction condition

$$\left ( \partial _ { x _ { 3 } } - 3 \partial _ { x _ { 1 } } \right ) \sigma = C \sigma,$$

where C is some constant (the reason for our choice of the coefficient 3 in front of ∂ x 1 will be explained later). Under this condition, the higher-dimensional bilinear equation (34) reduces to

$$\left ( D _ { x _ { 1 } } ^ { 4 } - 1 2 D _ { x _ { 1 } } ^ { 2 } + 3 D _ { x _ { 2 } } ^ { 2 } \right ) \sigma \cdot \sigma = 0.$$

Finally, we define and impose the reality condition

$$\sigma ^ { * } = \sigma.$$

Then, the bilinear equation (36) becomes the bilinear equation (33) of the Boussinesq equation, and σ (with x 3 = 0) becomes its algebraic (rogue wave) solution.

Next, we follow the above outline to derive general rogue-wave solutions to the Boussinesq equation (3).

$$x _ { 1 } = \frac { 1 } { 2 } x, \ \ x _ { 2 } = - \frac { 1 } { 4 } i t,$$

## A. Gram solutions for the higher-dimensional bilinear system

First, we derive algebraic solutions to the higher-dimensional bilinear equation (34). From Lemma 1 of Ref. [43], we learn that if functions m ( n ) i,j , ϕ ( n ) i and ψ ( n ) j of variables ( x 1 , x 2 , x 3 ) satisfy the following differential and difference relations

$$\begin{array} { c } \partial _ { x _ { 1 } } m _ { i, j } ^ { ( n ) } = \varphi _ { i } ^ { ( n ) } \psi _ { j } ^ { ( n ) }, \\ \partial _ { x _ { 2 } } m _ { i, j } ^ { ( n ) } = \varphi _ { i } ^ { ( n + 1 ) } \psi _ { j } ^ { ( n ) } + \varphi _ { i } ^ { ( n ) } \psi _ { j } ^ { ( n - 1 ) }, \\ \partial _ { x _ { 3 } } m _ { i, j } ^ { ( n ) } = \varphi _ { i } ^ { ( n + 2 ) } \psi _ { j } ^ { ( n ) } + \varphi _ { i } ^ { ( n + 1 ) } \psi _ { j } ^ { ( n - 1 ) } + \varphi _ { i } ^ { ( n ) } \psi _ { j } ^ { ( n - 2 ) }, \\ m _ { i, j } ^ { ( n + 1 ) } = m _ { i, j } ^ { ( n ) } + \varphi _ { i } ^ { ( n ) } \psi _ { j } ^ { ( n + 1 ) }, \\ \partial _ { x _ { k } } \varphi _ { i } ^ { ( n ) } = \varphi _ { i } ^ { ( n + k ) }, \, \partial _ { x _ { k } } \psi _ { j } ^ { ( n ) } = - \psi _ { j } ^ { ( n - k ) }, \, ( k = 1, 2, 3 ), \end{array} \right ) \text{min} \end{array}$$

then the determinant

$$\tau _ { n } = \det _ { 1 \leq i, j \leq N } \left ( m _ { i, j } ^ { ( n ) } \right )$$

satisfies the bilinear equation (34), i.e.,

$$\left ( D _ { x _ { 1 } } ^ { 4 } - 4 D _ { x _ { 1 } } D _ { x _ { 3 } } + 3 D _ { x _ { 2 } } ^ { 2 } \right ) \tau _ { n } \cdot \tau _ { n } = 0.$$

Next, we introduce functions m ( n ) , ϕ ( n ) and ψ ( n ) as

$$m ^ { ( n ) } = \frac { 1 } { p + q } \left ( - \frac { p } { q } \right ) ^ { n } e ^ { \xi + \eta }, \quad \varphi ^ { ( n ) } = p ^ { n } e ^ { \xi }, \quad \psi ^ { ( n ) } = ( - q ) ^ { - n } e ^ { \eta },$$

where

$$\xi = p x _ { 1 } + p ^ { 2 } x _ { 2 } + p ^ { 3 } x _ { 3 }, \quad \eta = q x _ { 1 } - q ^ { 2 } x _ { 2 } + q ^ { 3 } x _ { 3 }.$$

It is easy to check that these functions satisfy the differential and difference relations

$$t \text{ these functions satisfy the differential and difference relations} \\ \partial _ { x _ { 1 } } m ^ { ( n ) } & = \varphi ^ { ( n ) } \psi ^ { ( n ) }, \\ \partial _ { x _ { 2 } } m ^ { ( n ) } & = \varphi ^ { ( n + 1 ) } \psi ^ { ( n ) } + \varphi ^ { ( n ) } \psi ^ { ( n - 1 ) }, \\ \partial _ { x _ { 3 } } m ^ { ( n ) } & = \varphi ^ { ( n + 2 ) } \psi ^ { ( n ) } + \varphi ^ { ( n + 1 ) } \psi ^ { ( n - 1 ) } + \varphi ^ { ( n ) } \psi ^ { ( n - 2 ) }, \\ m ^ { ( n + 1 ) } & = m ^ { ( n ) } + \varphi ^ { ( n ) } \psi ^ { ( n + 1 ) }, \\ \partial _ { x _ { k } } \varphi ^ { ( n ) } & = \varphi ^ { ( n + k ) }, \quad \partial _ { x _ { k } } \psi ^ { ( n ) } = - \psi ^ { ( n - k ) }, \quad ( k = 1, 2, 3 ).$$

Therefore, by defining

$$m _ { i j } ^ { ( n ) } = \mathcal { A } _ { i } \mathcal { B } _ { j } m ^ { ( n ) }, \quad \varphi _ { i } ^ { ( n ) } = \mathcal { A } _ { i } \varphi ^ { ( n ) }, \quad \psi _ { j } ^ { ( n ) } = \mathcal { B } _ { j } \psi ^ { ( n ) },$$

where A i and B j are differential operators with respect to p and q respectively as

$$\mathcal { A } _ { i } & = \sum _ { k = 0 } ^ { i } \frac { a _ { k } } { ( i - k )! } \left [ f _ { 1 } ( p ) \partial _ { p } \right ] ^ { i - k } \\ \mathcal { B } _ { j } & = \sum _ { l = 0 } ^ { j } \frac { b _ { l } } { ( j - l )! } \left [ f _ { 2 } ( q ) \partial _ { q } \right ] ^ { j - l } \ \right ],$$

f 1 ( p ) and f 2 ( q ) are arbitrary functions of p and q , and a k , b k are arbitrary complex constants, we would see that these m ( n ) ij , ϕ ( n ) i and ψ ( n ) j obey the differential and difference relations (39) since the operators A i and B j commute with differentials ∂ x k . Then, Lemma 1 of Ref. [43] tells us that for an arbitrary sequence of indices ( i 1 , i 2 , · · · , i N ; j 1 , j 2 , · · · , j N ), the determinant

$$\tau _ { n } = \det _ { 1 \leq \nu, \mu \leq N } \left ( m _ { i _ { \nu }, j _ { \mu } } ^ { ( n ) } \right )$$

satisfies the higher-dimensional bilinear equation (41).

The above solutions (46) are a very broad class of algebraic solutions to the bilinear equation (41) which contain a huge amount of freedom. For instance, parameters p and q are totally arbitrary, so are the functions f 1 ( p ) and f 2 ( q ), complex constants a k and b k , as well as sequences of indices ( i 1 , i 2 , · · · , i N ; j 1 , j 2 , · · · , j N ). Only a small portion of such solutions can satisfy the dimension reduction condition (35) and the reality condition (38), which then make them algebraic solutions to the bilinear equation (33) of the the Boussinesq equation (3) under variable connections (37).

Next, we will restrict solutions (46) so that they meet the dimension reduction condition (35) and the reality condition (38).

## B. Dimensional reduction through the W -p treatment

We first condider the dimensional reduction (35) for the bilinear equation (41), which will introduce restrictions on parameters p and q , functions f 1 ( p ) and f 2 ( q ), as well as sequences of indices ( i 1 , i 2 , · · · , i N ; j 1 , j 2 , · · · , j N ). We will also explain the choice of the coefficient 3 in front of ∂ x 1 in that reduction (35).

Dimension reduction is a crucial step in the bilinear KP-reduction procedure. In the past, the functions f 1 ( p ) and f 2 ( q ) in this procedure were always chosen to be linear functions [26, 32, 33, 42, 49-51, 53]. In many cases, such choices were appropriate. However, for a coupled NLS-Boussinesq system and the Yajima-Oikawa system studied in [49, 50], dimension reduction under such choices became cumbersome and complicated. Indeed, the authors were forced to introduce two additional indices in the coefficients of the differential operators and matrix elements, and those coefficients had to be obtained from nontrivial recurrence relations. For the present Boussinesq equation (3), choices of linear functions for f 1 ( p ) and f 2 ( q ) would encounter a similar difficulty as in [49, 50].

In this article, we will choose functions f 1 ( p ) and f 2 ( q ) differently. We will show that under our judicious choices of f 1 ( p ) and f 2 ( q ), dimension reduction will simplify dramatically. As a result, recurrence relations for coefficients in the solution will be eliminated, and the solution expressions of rogue waves will become more clean and concise.

We start from a general dimension reduction condition

$$\left ( \partial _ { x _ { 3 } } + \alpha \partial _ { x _ { 1 } } \right ) \tau _ { n } = C \tau _ { n },$$

where α and C are constants to be determined. In order to calculate the left side of the above equation, we notice from the definitions (42) and (44) of m ( n ) and m ( n ) i,j that

$$\left ( \partial _ { x _ { 3 } } + \alpha \partial _ { x _ { 1 } } \right ) m _ { i, j } ^ { ( n ) } = \mathcal { A } _ { i } \mathcal { B } _ { j } \left [ Q _ { 1 } ( p ) + Q _ { 2 } ( q ) \right ] m ^ { ( n ) },$$

where

$$Q _ { 1 } ( p ) = p ^ { 3 } + \alpha p, \ \ Q _ { 2 } ( q ) = q ^ { 3 } + \alpha q.$$

To proceed, we introduce new variables W 1 and W 2 through

$$Q _ { 1 } ( p ) = \mathcal { W } _ { 1 } ( p ) + \frac { 1 } { \mathcal { W } _ { 1 } ( p ) }, \ \ Q _ { 2 } ( q ) = \mathcal { W } _ { 2 } ( q ) + \frac { 1 } { \mathcal { W } _ { 2 } ( q ) }.$$

In terms of these new variables, our new choices of functions f 1 ( p ) and f 2 ( q ) in the differential operators A i and B j are

$$f _ { 1 } ( p ) = \frac { \mathcal { W } _ { 1 } ( p ) } { \mathcal { W } _ { 1 } ^ { \prime } ( p ) }, \ \ f _ { 2 } ( q ) = \frac { \mathcal { W } _ { 2 } ( q ) } { \mathcal { W } _ { 2 } ^ { \prime } ( q ) }.$$

More explicit expressions for f 1 ( p ) and f 2 ( q ) will be provided later after the value of α is ascertained [see Eq. (69)]. The motivation behind this choice of f 1 ( p ) is that, under this choice,

$$f _ { 1 } ( p ) \partial _ { p } = \frac { \mathcal { W } _ { 1 } ( p ) } { \mathcal { W } _ { 1 } ^ { \prime } ( p ) } \partial _ { p } = \mathcal { W } _ { 1 } ( p ) \partial _ { \mathcal { W } _ { 1 } ( p ) }.$$

Thus,

$$\mathcal { A } _ { i } Q _ { 1 } ( p ) m ^ { ( n ) } = \left [ \sum _ { k = 0 } ^ { i } \frac { a _ { k } } { ( i - k )! } \left ( \mathcal { W } _ { 1 } \partial _ { \mathcal { W } _ { 1 } } \right ) ^ { i - k } \right ] \left ( \mathcal { W } _ { 1 } + \frac { 1 } { \mathcal { W } _ { 1 } } \right ) m ^ { ( n ) }.$$

One can recognize that the operators on the right side of this equation are the same as those in [26], except for a notation of W 1 instead of p . Using results of [26], we immediately get

$$\mathcal { A } _ { i } Q _ { 1 } ( p ) m ^ { ( n ) } = \sum _ { k = 0 } ^ { i } \frac { 1 } { k! } \left [ \mathcal { W } _ { 1 } ( p ) + ( - 1 ) ^ { k } \frac { 1 } { \mathcal { W } _ { 1 } ( p ) } \right ] \mathcal { A } _ { i - k } m ^ { ( n ) }.$$

For exactly the same reasons, we also have

$$\mathcal { B } _ { j } Q _ { 2 } ( q ) m ^ { ( n ) } = \sum _ { l = 0 } ^ { j } \frac { 1 } { l! } \left [ \mathcal { W } _ { 2 } ( q ) + ( - 1 ) ^ { l } \frac { 1 } { \mathcal { W } _ { 2 } ( q ) } \right ] \mathcal { B } _ { j - l } m ^ { ( n ) }.$$

Substituting these two equations into (48), we then get

$$( \partial _ { x _ { 3 } } + \alpha \partial _ { x _ { 1 } } ) \, m _ { i, j } ^ { ( n ) } = \sum _ { k = 0 } ^ { i } \frac { 1 } { k! } \left [ \mathcal { W } _ { 1 } ( p ) + ( - 1 ) ^ { k } \frac { 1 } { \mathcal { W } _ { 1 } ( p ) } \right ] m _ { i - k, j } ^ { ( n ) } + \sum _ { l = 0 } ^ { j } \frac { 1 } { l! } \left [ \mathcal { W } _ { 2 } ( q ) + ( - 1 ) ^ { l } \frac { 1 } { \mathcal { W } _ { 2 } ( q ) } \right ] m _ { i, j - l } ^ { ( n ) }.$$

Now, it is time to select values of p , q and α so that the above equation can be further simplified (the selected values of p and q will be denoted as p 0 and q 0 ). Since treatments for p 0 and q 0 values are the same, we will consider p 0 only. Motivated by [26], we require W p 1 ( 0 ) = 1, so that the oddk terms in the above summation drop out. This W p 1 ( 0 ) = 1 condition leads to Q 1 ( p 0 ) = 2 in view of Eq. (50). Differentiating the first equation in (50) with respect to p , we get

$$\mathcal { W } _ { 1 } ^ { \prime } ( p ) = \frac { Q _ { 1 } ^ { \prime } ( p ) } { 1 - \mathcal { W } _ { 1 } ^ { - 2 } ( p ) }.$$

At the selected p 0 value where W p 1 ( 0 ) = 1, W ′ 1 ( p 0 ) needs to be well defined in view of the definition of f 1 ( p ) in Eq. (51). Then, the above equation requires Q ′ 1 ( p 0 ) = 0. Thus, the two conditions for the p 0 and α values are

$$Q _ { 1 } ( p _ { 0 } ) = 2, \ \ Q _ { 1 } ^ { \prime } ( p _ { 0 } ) = 0.$$

Inserting the Q 1 ( p ) function from (49) into these two constraints, we get

$$p _ { 0 } ^ { 3 } + \alpha p _ { 0 } = 2, \ \ 3 p _ { 0 } ^ { 2 } + \alpha = 0,$$

whose solutions are

$$p _ { 0 } = - 1, \quad \alpha = - 3.$$

This explains the coefficient 3 we have chosen in the dimension reduction condition (35). The same consideration for q 0 leads to q 0 = -1.

Under the above choices of p, q and α values, Eq. (56) further simplifies to

$$( \partial _ { x _ { 3 } } - 3 \partial _ { x _ { 1 } } ) \, m _ { i, j } ^ { ( n ) } \Big | _ { p = - 1, \ q = - 1 } = 2 \sum _ { \substack { k = 0, \\ k \colon e v e n } } ^ { i } \frac { 1 } { k! } \, m _ { i - k, j } ^ { ( n ) } \Big | _ { p = - 1, \ q = - 1 } + 2 \sum _ { \substack { l = 0, \\ l \colon e v e n } } ^ { j } \frac { 1 } { l! } \, m _ { i, j - l } ^ { ( n ) } \Big | _ { p = - 1, \ q = - 1 }.$$

This is an important relation which shows that, at the selected ( p, q ) values, ( ∂ x 3 -3 ∂ x 1 ) m ( n ) i,j is a linear combination of m ( n ) i,j and other m ( n ) ˆ ˆ i,j terms of lower row and column indices with jumps of 2. Due to this relation, if we choose indices ( i 1 , i 2 , · · · , i N ; j 1 , j 2 , · · · , j N ) in the determinant (46) as

$$\tau _ { n } = \det _ { 1 \leq i, j \leq N } \left ( m _ { 2 i - 1, 2 j - 1 } ^ { ( n ) } \right ),$$

then the same calculation as in [26] would show that this τ n function satisfies the dimension reduction condition

$$\left ( \partial _ { x _ { 3 } } - 3 \partial _ { x _ { 1 } } \right ) \tau _ { n } = 4 N \tau _ { n }.$$

We note by passing that another index choice of

$$\tau _ { n } = \det _ { 1 \leq i, j \leq N } \left ( m _ { 2 i - 2, 2 j - 2 } ^ { ( n ) } \right )$$

would also satisfy the dimension reduction condition (63). But as we have shown in a different but similar context [53], this other index choice would lead to solutions which are equivalent to those from (62).

When the dimension reduction condition (63) is substituted into the higher-dimensional bilinear equation (41) and setting n = 0 and x 3 = 0, we get

$$\left ( D _ { x _ { 1 } } ^ { 4 } - 1 2 D _ { x _ { 1 } } ^ { 2 } + 3 D _ { x _ { 2 } } ^ { 2 } \right ) \tau _ { 0 } \cdot \tau _ { 0 } = 0.$$

Thus, the third dimension x 3 has been eliminated, and dimension reduction has been completed.

During the above dimension reduction, the α value in (47) has been ascertained [see (60)]. Thus, we can now derive more explicit expressions for functions f 1 ( p ) and f 2 ( q ) as defined in Eq. (51). From Eq. (50), we have

$$\left [ \mathcal { W } _ { 1 } ( p ) - \frac { 1 } { \mathcal { W } _ { 1 } ( p ) } \right ] ^ { 2 } = Q _ { 1 } ^ { 2 } ( p ) - 4.$$

Differentiating the first equation in (50) with respect to p , we get

$$\frac { \mathcal { W } _ { 1 } ^ { \prime } ( p ) } { \mathcal { W } _ { 1 } ( p ) } \left [ \mathcal { W } _ { 1 } ( p ) - \frac { 1 } { \mathcal { W } _ { 1 } ( p ) } \right ] = Q _ { 1 } ^ { \prime } ( p ).$$

Utilizing these two equations and the definition of f 1 ( p ) in (51), we get

$$f _ { 1 } ( p ) = \frac { \sqrt { Q _ { 1 } ^ { 2 } ( p ) - 4 } } { Q _ { 1 } ^ { \prime } ( p ) }.$$

In view of the definition of Q 1 ( p ) in (49) and the α value in (47), the above f 1 ( p ) formula can be further simplified. Following the same calculation, the f 2 ( q ) function can also be derived. The final results are

$$f _ { 1 } ( p ) = \frac { \sqrt { p ^ { 2 } - 4 } } { 3 }, \ \ f _ { 2 } ( q ) = \frac { \sqrt { q ^ { 2 } - 4 } } { 3 }.$$

## C. The reality condition

In the bilinear equation (65), when the ( x , x 1 2 ) variables are linked to ( x, t ) through (37), i.e., x 1 = x/ 2 and x 2 = -i t/ 4, then this bilinear equation becomes (33) of the Boussinesq equation. The only remaining condition is the reality condition (38), i.e.,

$$\tau _ { 0 } ^ { * } & = \tau _ { 0 }.$$

Notice that when p = q = -1, f 1 ( p ) and f 2 ( q ) are purely imaginary in view of Eq. (69). Hence,

$$\left ( \left [ f _ { 1 } ( p ) \partial _ { p } \right ] ^ { k } \right ) ^ { * } = ( - 1 ) ^ { k } \left [ f _ { 1 } ( p ) \partial _ { p } \right ] ^ { k }, \quad \left ( \left [ f _ { 2 } ( q ) \partial _ { q } \right ] ^ { l } \right ) ^ { * } = ( - 1 ) ^ { l } \left [ f _ { 2 } ( q ) \partial _ { q } \right ] ^ { l }.$$

Thus, if we constrain parameters a k and b k by

$$b _ { k } = ( - 1 ) ^ { k } a _ { k } ^ { * },$$

then since x 1 = x/ 2 is real and x 2 = -i t/ 4 is imaginary, we can easily show that

$$\left [ m _ { i, j } ^ { ( 0 ) } \right ] ^ { * } = ( - 1 ) ^ { i + j } m _ { j, i } ^ { ( 0 ) },$$

and therefore the reality condition (70) holds.

Regarding the regularity of solutions given in Theorem 1, notice from Eqs. (6) and (73) that σ x, t ( ) is the determinant of a Hermitian matrix M = mat 1 ≤ i,j ≤ N ( m 2 i -1 2 , j -1 ). Then, using techniques similar to that in [26], we can show that the matrix M is positive definite, so that σ x, t ( ) &gt; 0, which proves that the solution u x, t ( ) is nonsingular. Summarizing the above results, Theorem 1 is then proved, except for the boundary conditions (4), which will be discussed in the end of the next subsection.

## D. Algebraic representation of rogue waves

Finally, we derive purely algebraic expressions of rogue waves and prove Theorem 2. These algebraic expressions are useful for multiple purposes. For instance, they can produce explicit formulae of rogue waves much more quickly than the expressions in Theorem 1 where repeated differentiations have to be performed. For another instance, these algebraic expressions allow us to derive highest-power polynomial terms of rogue waves [26, 52, 53], which can be used to analytically prove the boundary conditions (4) of these solutions.

The basic idea of this derivation is the same as that in [26], except that we will work with variables W 1 ( p ) and W 2 ( q ) instead of p and q . Introducing the generator G of the differential operators ( f 1 ( p ∂ ) p ) k ( f 2 ( q ∂ ) q ) l as

$$\mathcal { G } = \sum _ { k = 0 } ^ { \infty } \sum _ { l = 0 } ^ { \infty } \frac { \kappa ^ { k } } { k! } \frac { \lambda ^ { l } } { l! } \left [ f _ { 1 } ( p ) \partial _ { p } \right ] ^ { k } \left [ f _ { 2 } ( q ) \partial _ { q } \right ] ^ { l },$$

and utilizing Eq. (52), we get

$$\mathcal { G } = \sum _ { k = 0 } ^ { \infty } \sum _ { l = 0 } ^ { \infty } \frac { \kappa ^ { k } } { k! } \frac { \lambda ^ { l } } { l! } \left [ \partial _ { \ln } \wta _ { 1 } \right ] ^ { k } \left [ \partial _ { \ln } \wta _ { 2 } \right ] ^ { l } = \exp \left ( \kappa \partial _ { \ln } \wta _ { 1 } + \lambda \partial _ { \ln } \wta _ { 2 } \right ).$$

Thus, for any function F ( W W 1 , 2 ), we have [26]

$$\mathcal { G } F ( \mathcal { W } _ { 1 }, \mathcal { W } _ { 2 } ) = F ( e ^ { \kappa } \mathcal { W } _ { 1 }, e ^ { \lambda } \mathcal { W } _ { 2 } ).$$

Next, we will apply this generator on the function m (0) . After dimension reduction ( x 3 = 0) and variable relations (37), this m (0) reduces from (42) to

$$m ^ { ( 0 ) } = \frac { 1 } { p + q } e ^ { \frac { 1 } { 2 } ( p + q ) x - \frac { 1 } { 4 } ( p ^ { 2 } - q ^ { 2 } ) i t }.$$

To utilize Eq. (76), we need to express p and q in this m (0) as functions of W 1 and W 2 . Equations (49)-(50) tell us that

$$p ^ { 3 } - 3 p = \mathcal { W } _ { 1 } + \frac { 1 } { \mathcal { W } _ { 1 } }, \quad q ^ { 3 } - 3 q = \mathcal { W } _ { 2 } + \frac { 1 } { \mathcal { W } _ { 2 } }.$$

These equations for p and q can be solved and there are three roots. Due to our earlier conditions that W 1 = W 2 = 1 when p = q = -1, the suitable roots for p and q are

$$p ( \mathcal { W } _ { 1 } ) = c _ { 1 } \mathcal { W } _ { 1 } ^ { 1 / 3 } + c _ { 2 } \mathcal { W } _ { 1 } ^ { - 1 / 3 }, \ \ q ( \mathcal { W } _ { 2 } ) = c _ { 1 } \mathcal { W } _ { 2 } ^ { 1 / 3 } + c _ { 2 } \mathcal { W } _ { 2 } ^ { - 1 / 3 },$$

∗

Now, we apply Eq. (76) on m and get where c 1 = exp(2i π/ 3) and c 2 = c 1 . (0)

$$\mathcal { G } m ^ { ( 0 ) } = \frac { 1 } { p ( e ^ { \kappa } \mathcal { W } _ { 1 } ) + q ( e ^ { \lambda } \mathcal { W } _ { 2 } ) } \exp \left ( \frac { 1 } { 2 } \left [ p ( e ^ { \kappa } \mathcal { W } _ { 1 } ) + q ( e ^ { \lambda } \mathcal { W } _ { 2 } ) \right ] x - \frac { 1 } { 4 } \left [ p ^ { 2 } ( e ^ { \kappa } \mathcal { W } _ { 1 } ) - q ^ { 2 } ( e ^ { \lambda } \mathcal { W } _ { 2 } ) \right ] i t \right ).$$

At p = q = -1, W 1 = W 2 = 1. Thus,

$$\frac { 1 } { m ^ { ( 0 ) } } \, \mathcal { G } m ^ { ( 0 ) } \Big | _ { p = q = - 1 } = \frac { ( - 2 ) } { p ( e ^ { \kappa } ) + q ( e ^ { \lambda } ) } \exp \left ( \frac { 1 } { 2 } \left [ p ( e ^ { \kappa } ) + q ( e ^ { \lambda } ) + 2 \right ] x - \frac { 1 } { 4 } \left [ p ^ { 2 } ( e ^ { \kappa } ) - q ^ { 2 } ( e ^ { \lambda } ) \right ] i t \right ).$$

We need to expand the right side of this equation into double Taylor series in κ and λ . To expand the fraction in front of the exponential term, we notice that for any functions f ( κ ) and g λ ( ),

$$& \frac { - 2 } { f ( \kappa ) + g ( \lambda ) } = \frac { - 2 \left [ f ( 0 ) + g ( 0 ) \right ] } { \left [ f ( \kappa ) + g ( 0 ) \right ] \left [ g ( \lambda ) + f ( 0 ) \right ] } \frac { \cdot \cdot \cdot } { 1 - \frac { f ( \kappa ) - f ( 0 ) } { f ( \kappa ) + g ( 0 ) } \, \frac { g ( \lambda ) - g ( 0 ) } { g ( \lambda ) + f ( 0 ) } } \\ & = \exp \left ( - \ln \frac { \left [ f ( \kappa ) + g ( 0 ) \right ] \left [ g ( \lambda ) + f ( 0 ) \right ] } { - 2 \left [ f ( 0 ) + g ( 0 ) \right ] } \right ) \sum _ { \nu = 0 } ^ { \infty } \left [ \frac { f ( \kappa ) - f ( 0 ) } { f ( \kappa ) + g ( 0 ) } \, \frac { g ( \lambda ) - g ( 0 ) } { g ( \lambda ) + f ( 0 ) } \right ] ^ { \nu }.$$

$$= \exp \left ( - \ln \frac { \left [ f ( \kappa ) + g ( 0 ) \right ] \left [ g ( \lambda ) + f ( 0 ) \right ] } { - 2 \left [ f ( 0 ) + g ( 0 ) \right ] } \right ) \sum _ { \nu = 0 } ^ { \infty } \left [ \frac { f ( \kappa ) - f ( 0 ) } { f ( \kappa ) + g ( 0 ) } \, \frac { g ( \lambda ) - g ( 0 ) } { g ( \lambda ) + f ( 0 ) } \right ] ^ { \nu }.$$

Thus, substituting

$$f ( \kappa ) = p ( e ^ { \kappa } ) = c _ { 1 } e ^ { \kappa / 3 } + c _ { 2 } e ^ { - \kappa / 3 }, \ \ g ( \lambda ) = q ( e ^ { \lambda } ) = c _ { 1 } e ^ { \lambda / 3 } + c _ { 2 } e ^ { - \lambda / 3 }$$

into the above equation, we get

$$\L o v e \equation, \, & \text{ve} \, \text{get} \\ & \exp \left ( - \ln \frac { [ f ( \kappa ) + g ( 0 ) ] \, [ g ( \lambda ) + f ( 0 ) ] } { - 2 \, [ f ( 0 ) + g ( 0 ) ] } \right ) \\ & = \exp \left ( - \ln \left [ \frac { 1 } { 2 } - \cosh \left ( \frac { \kappa } { 3 } + \frac { 2 i \pi } { 3 } \right ) \right ] - \ln \left [ \frac { 1 } { 2 } - \cosh \left ( \frac { \lambda } { 3 } + \frac { 2 i \pi } { 3 } \right ) \right ] \right ) \\ & = \exp \left ( \sum _ { k = 1 } ^ { \infty } r _ { k } \left ( \kappa ^ { k } + \lambda ^ { k } \right ) \right ),$$

where r k are coefficients of the Taylor expansion given by (15) in Theorem 2. In addition, we get

$$& \text{where} \, r _ { k } \, \text{are coefficients of the $1$aylor expansion given by (15) in $1$Theorem $2$.  In addition, we get} \\ & \quad \sum _ { \nu = 0 } ^ { \infty } \left [ \frac { f ( \kappa ) - f ( 0 ) } { f ( \kappa ) + g ( 0 ) } \frac { g ( \lambda ) - g ( 0 ) } { g ( \lambda ) + f ( 0 ) } \right ] ^ { \nu } \\ & \quad = \sum _ { \nu = 0 } ^ { \infty } \left ( - \frac { \kappa \lambda } { 1 2 } \right ) ^ { \nu } \exp \left ( \nu \ln \left [ \frac { 2 i \sqrt { 3 } } { \kappa } \tanh \frac { \kappa } { 6 } \tanh \left ( \frac { \kappa } { 6 } + \frac { 2 i \pi } { 3 } \right ) \right ] + \nu \ln \left [ \frac { 2 i \sqrt { 3 } } { \lambda } \tanh \frac { \lambda } { 6 } \tanh \left ( \frac { \lambda } { 6 } + \frac { 2 i \pi } { 3 } \right ) \right ] \right ) \\ & \quad = \sum _ { \nu = 0 } ^ { \infty } \left ( - \frac { \kappa \lambda } { 1 2 } \right ) ^ { \nu } \exp \left ( \nu \sum _ { k = 1 } ^ { \infty } s _ { k } \left ( \kappa ^ { k } + \lambda ^ { k } \right ) \right ), \\ & \text{where} \, s _ { k } \, \text{are defined by} \, ( 16 ) \, \text{in Theorem $2$.  Regarding the exponential term in $Eq.$  (81).  when functions $p(-} ) \, \text{and} \, q(-} \right )$$

where s k are defined by (16) in Theorem 2. Regarding the exponential term in Eq. (81), when functions p ( ) and · q ( ) · in Eq. (79) are inserted, we get

$$\left ( \frac { 1 } { 2 } \left [ p ( e ^ { \kappa } ) + q ( e ^ { \kappa } ) + 2 \right ] x - \frac { 1 } { 4 } \left [ p ^ { 2 } ( e ^ { \kappa } ) - q ^ { 2 } ( e ^ { \kappa } ) \right ] i \right ) \\ = & \exp \left [ \frac { 1 } { 2 } \left ( c _ { 1 } e ^ { \kappa / 3 } + c _ { 2 } e ^ { - \kappa / 3 } + c _ { 1 } e ^ { \lambda / 3 } + c _ { 2 } e ^ { - \lambda / 3 } + 2 \right ) x - \frac { 1 } { 4 } \left ( c _ { 2 } e ^ { 2 \kappa / 3 } + c _ { 1 } e ^ { - 2 \kappa / 3 } - c _ { 2 } e ^ { 2 \lambda / 3 } - c _ { 1 } e ^ { - 2 \lambda / 3 } \right ) i \right ] \\ = & \exp \left ( \sum _ { k = 1 } ^ { \infty } \frac { \kappa ^ { k } } { k! } \frac { c _ { 1 } + ( - 1 ) ^ { k } c _ { 2 } } { 2 \cdot 3 ^ { k } } \left [ x + ( - 2 ) ^ { k - 1 } i t \right ] + \sum _ { k = 1 } ^ { \infty } \frac { \lambda ^ { k } } { k! } \frac { c _ { 1 } + ( - 1 ) ^ { k } c _ { 2 } } { 2 \cdot 3 ^ { k } } \left [ x - ( - 2 ) ^ { k - 1 } i t \right ] \right ). \\ \tilde { \tau } \text{minin} \, \text{all theoo non-title} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$E_{n}$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$} \, \text{$I$}}$$

Combining all these results, Eq. (81) reduces to

$$\frac { 1 } { m ^ { ( 0 ) } } \, \mathcal { G } m ^ { ( 0 ) } \Big | _ { p = q = - 1 } = \sum _ { \nu = 0 } ^ { \infty } \left ( - \frac { \kappa \lambda } { 1 2 } \right ) ^ { \nu } \exp \left ( \sum _ { k = 1 } ^ { \infty } \left ( x _ { k } ^ { + } + \nu s _ { k } \right ) \kappa ^ { k } + \sum _ { k = 1 } ^ { \infty } \left ( x _ { k } ^ { - } + \nu s _ { k } \right ) \lambda ^ { k } \right ),$$

where x ± k are given by (14) in Theorem 2. Then, taking the coefficients of κ λ k l on both sides, we find

$$\frac { 1 } { m ^ { ( 0 ) } } \frac { 1 } { k! l! } \left [ f _ { 1 } ( p ) \partial _ { p } \right ] ^ { k } \left [ f _ { 2 } ( q ) \partial _ { q } ^ { l } \ m ^ { ( 0 ) } \right ] _ { p = - 1, q = - 1 } = \sum _ { \nu = 0 } ^ { \min ( k, l ) } \left ( - \frac { \kappa \lambda } { 1 2 } \right ) ^ { \nu } S _ { k = - \nu } \left ( x ^ { + } + \nu s \right ) S _ { l = \nu } \left ( x ^ { - } + \nu s \right ).$$

From this result, we get

$$\frac { 1 } { m ^ { ( 0 ) } } \mathcal { A } _ { i } \mathcal { B } _ { j } \ m ^ { ( 0 ) } \Big | _ { p = - 1, q = - 1 } = \sum _ { \nu = 0 } ^ { \min ( i, j ) } \Phi _ { i \nu } \Psi _ { j \nu },$$

where Φ iν and Ψ jν are defined in Eqs. (12)-(13) of Theorem 2. Since the Boussinesq solution (5) is invariant when the function σ is multiplied by an exponential factor of a linear function in x and t , the above expression (with the exponential factor 1 /m (0) ) as the matrix element of σ , rather than the original A B i j m (0) in Eq. (7), would give the same solution. This finishes the proof of Theorem 2.

In the end, we show the boundary conditions (4) of these rogue waves. For this purpose, we first rewrite the determinant (6) of σ x, t ( ) with Schur-polynomial matrix elements in Theorem 2 into a 3 N × 3 N determinant, as was done for the NLS equation in [26]. Using this larger determinant, one can directly obtain the leading-order terms in the polynomial function σ x, t ( ) as

$$\sigma ( x, t ) = c _ { N } \left ( x ^ { 2 } + t ^ { 2 } \right ) ^ { N ( N + 1 ) / 2 } + \text{lower-degree terms},$$

where c N is a N -dependent constant. From this result, one can see quickly that the solution u x, t ( ) in Eq. (5) satisfies the boundary conditions (4).

## V. CONCLUSION AND DISCUSSION

In this article, we have derived rogue wave solutions in the Boussineq equation (3) through the bilinear KPreduction method, and these solutions are given explicitly as Gram determinants with matrix elements in terms of Schur polynomials. Our solutions contain more free parameters than those reported before [4], and they exhibit new

interesting wave patterns. We have also shown that the rogue wave of maximum amplitude at each order is generally asymmetric in space, which is quite unusual in integrable equations.

Technically, our main contribution to the bilinear KP-reduction method for rogue waves is a new judicious choice of differential operators for matrix elements (the W -p treatment). Compared to the previous choices, this W -p treatment drastically simplifies the dimension reduction calculation as well as the analytical expressions of rogue wave solutions. Can this W -p treatment be applied to other integrable equations, such as the NLS-Boussinesq equation and the Yajima-Oikawa system considered in [49, 50]? The answer is a definite yes, maybe with minor modifications possibly. For instance, for the Yajima-Oikawa system, we can slightly modify our definitions of the W 1 ( p ) and W 2 ( q ) functions as Q k = β k [ W k + W -1 k ], where Q 1 ( p ) and Q 2 ( q ) are certain functions arising from the dimension reduction of the Yajima-Oikawa system, and β 1 , β 2 are certain constants. For these W 1 ( p ) and W 2 ( q ) functions, the differential operators (45) with f 1 ( p ) and f 2 ( q ) defined in (51) will again simplify the dimension reduction calculation and produce rogue wave expressions similar to those in this article. Comparatively, the rogue wave expressions derived in [49, 50] with old choices of differential operators had to involve additional indices in the coefficients of the differential operators and matrix elements, and those coefficients had to be obtained by complicated recurrence relations. This W -p treatment can also be applied to the NLS equation. In this special case, we would get f 1 ( p ) = p and f 2 ( q ) = q , which reproduces the old differential operators used in Ref. [26]. Thus, this W -p treatment is a general and useful technique to streamline the bilinear derivation of rogue waves in integrable systems when the dimension reduction is needed.

## Acknowledgement

This material is based upon work supported by the Air Force Office of Scientific Research under award number FA9550-18-1-0098 and the National Science Foundation under award number DMS-1616122.

- [1] J. Boussinesq, 'Th´ eorie de l'intumescence liquide, appel´e onde solitaire ou de translation se propagente dans un canal e rectangulaire', Comptes Rendus de l'Acad´mie des Sciences 72, 755-759 (1871). e
- [2] J. Boussinesq, 'Th´orie des ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal, en communie quant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond', Journal de Math´matiques e Pures et Appliqu´es 17, 55-108 (1872). e
- [3] F. Ursell, 'The long-wave paradox in the theory of gravity waves', Proc. Camb. Phil. Soc. 49, 685-694 (1953).
- [4] P. A. Clarkson, E. Dowie, 'Rational solutions of the Boussinesq equation and applications to rogue waves,' Transactions of Mathematics and Its Applications 1, 1-26 (2017).
- [5] V.E. Zakharov, 'On stocastization of one-dimensional chains of nonlinear oscillations', Zh. Eksp. Teor. Fiz. 65, 219-225 (1973) [Sov. Phys. JETP 38, 108-110 (1974)].
- [6] M. Toda, 'Studies of a nonlinear lattice', Phys. Rep. 8, 1-125 (1975).
- [7] N.J. Zabusky, 'A synergetic approach to problems of nonlinear dispersive wave propagation and interaction', in Nonlinear Partial Differential Equations [W. F. Ames, Editor] (Academic Press, New York, 1967, pp. 233-258).
- [8] E. Infeld and G. Rowlands, Nonlinear Waves, Solitons and Chaos (Cambridge University Press, Cambridge, 1990).
- [9] A.C. Scott, 'The application of B¨cklund transforms to physical problems', in a B¨ acklund Transformations [R.M. Miura, Editor] Lect. Notes Math. Vol. 515 (Springer-Verlag, Berlin, 1975, pp. 80-105).
- [10] R. Hirota, 'Exact N-soliton solution of the wave equation of long waves in shallow-water and in nonlinear lattices', J. Math. Phys. 14, 810-815 (1973).
- [11] P. Deift, C. Tomei and E. Trubowitz, 'Inverse scattering and the Boussinesq equation', Commun. Pure Appl. Math. 35, 567-628 (1982).
- [12] M.J. Ablowitz and R. Haberman, 'Resonantly coupled nonlinear evolution equations', J. Math. Phys. 16, 2301-2305 (1975).
- [13] M.J. Ablowitz and P.A. Clarkson, Solitons, Nonlinear Evolution Equations and Inverse Scattering (Cambridge University Press, Cambridge, 1991).
- [14] N. Akhmediev, A. Ankiewicz and M. Taki, 'Waves that appear from nowhere and disappear without a trace', Phys. Lett. A 373, 675-678 (2009).
- [15] K. Dysthe, H. E. Krogstad and P. M¨ller, 'Oceanic Rogue Waves', Annu. Rev. Fluid Mech. 40, 287-310 (2008). u
- [16] C. Kharif, E. Pelinovsky, and A. Slunyaev, Rogue Waves in the Ocean (Springer, Berlin, 2009).
- [17] D.R. Solli, C. Ropers, P. Koonath and B. Jalali, 'Optical rogue waves', Nature 450, 1054-1057 (2007).
- [18] S. Wabnitz, Nonlinear Guided Wave Optics: A Testbed for Extreme Waves (IOP Publishing, Bristol, 2017).

- [19] D.H. Peregrine, 'Water waves, nonlinear Schrodinger equations and their solutions', J. Aust. Math. Soc. B 25, 16-43 (1983).
- [20] N. Akhmediev, A. Ankiewicz and J. M. Soto-Crespo, 'Rogue waves and rational solutions of the nonlinear Schr¨dinger o equation,' Phys. Rev. E 80, 026601 (2009).
- [21] P. Dubard, P. Gaillard, C. Klein. and V. B. Matveev, 'On multi-rogue wave solutions of the NLS equation and positon solutions of the KdV equation', Eur. Phys. J. Spec. Top. 185, 247-258 (2010).
- [22] A. Ankiewicz, P. A. Clarkson and N. Akhmediev, 'Rogue waves, rational solutions, the patterns of their zeros and integral relations', J. Phys. A 43, 122002 (2010).
- [23] P. Dubard and V. B. Matveev, 'Multi-rogue waves solutions to the focusing NLS equation and the KP-I equation', Nat. Hazards Earth Syst. Sci. 11, 667-672 (2011).
- [24] D. J. Kedziora , A. Ankiewicz and N. Akhmediev, 'Circular rogue wave clusters', Phys. Rev. E 84, 056611 (2010).
- [25] B. L. Guo, L. M. Ling and Q. P. Liu, 'Nonlinear Schr¨dinger equation: o generalized Darboux transformation and rogue wave solutions,' Phys. Rev. E 85, 026607 (2012).
- [26] Y. Ohta and J. Yang, 'General high-order rogue waves and their dynamics in the nonlinear Schr¨dinger equation', Proc. o R. Soc. Lond. A 468, 1716-1740 (2012).
- [27] P. Dubard and V. B. Matveev, 'Multi-rogue waves solutions: from the NLS to the KP-I equation,' Nonlinearity 26, R93-R125 (2103).
- [28] L. H. Wang, C. Yang, J. Wang, and J. S. He, 'The height of an n th-order fundamental rogue wave for the nonlinear Schr¨dinger equation', Phys. Lett. A 381, 1714-1718 (2017). o
- [29] S. W. Xu, J. S. He and L. H. Wang, 'The Darboux transformation of the derivative nonlinear Schr¨dinger equation,' J. o Phys. A 44, 305203 (2011).
- [30] B. L. Guo, L. M. Ling and Q. P. Liu, 'High-order solutions and generalized Darboux transformations of derivative nonlinear Schr¨dinger equations', Stud. Appl. Math. 130, 317-344 (2013). o
- [31] F. Baronio, M. Conforti, A. Degasperis and S. Lombardo, 'Rogue waves emerging from the resonant interaction of three waves,' Phys. Rev. Lett. 111, 114101 (2103).
- [32] Y. Ohta and J. Yang, 'Rogue waves in the Davey-Stewartson I equation', Phys. Rev. E 86, 036604 (2012).
- [33] Y. Ohta and J. Yang, 'Dynamics of rogue waves in the Davey-Stewartson II equation', J. Phys. A 46, 105202 (2013).
- [34] A. Ankiewicz, N. Akhmediev and J. M. Soto-Crespo, 'Discrete rogue waves of the Ablowitz-Ladik and Hirota equations,' Phys. Rev. E 82, 026602 (2010).
- [35] A. Ankiewicz, J. M. Soto-Crespo and N. Akhmediev, 'Rogue waves and rational solutions of the Hirota equation,' Phys. Rev. E 81, 046602 (2010).
- [36] Y. S. Tao and J. S. He, 'Multisolitons, breathers, and rogue waves for the Hirota equation generated by the Darboux transformation', Phys. Rev. E 85, 026601 (2012).
- [37] F. Baronio, A. Degasperis, M. Conforti and S. Wabnitz, 'Solutions of the vector nonlinear Schr¨dinger equations: evidence o for deterministic rogue waves', Phys. Rev. Lett. 109, 044102 (2012).
- [38] K. W. Chow , H. N. Chan , D. J. Kedziora and R. H. J. Grimshaw, 'Rogue wave modes for the long wave-short wave resonance model,' J. Phys. Soc. Jpn. 82, 074001 (2013).
- [39] A. Degasperis and S. Lombardo, 'Rational solitons of wave resonant-interaction models,' Phys. Rev. E 88, 052914 (2013).
- [40] N. V. Priya, M. Senthilvelan and M. Lakshmanan, 'Akhmediev breathers, Ma solitons, and general breathers from rogue waves: a case study in the Manakov system,' Phys. Rev. E 88, 022918 (2013).
- [41] F. Baronio, M. Conforti, A. Degasperis, S. Lombardo, M. Onorato and S. Wabnitz, 'Vector rogue waves and baseband modulation instability in the defocusing regime', Phys. Rev. Lett. 113, 034101 (2014).
- [42] Y. Ohta and J. Yang, 'General rogue waves in the focusing and defocusing Ablowitz-Ladik equations,' J. Phys. A 47, 255201 (2014).
- [43] G. Mu and Z. Qin, 'Two spatial dimensional N -rogue waves and their dynamics in Melnikov equation', Nonlinear Analysis: Real World Applications. 18, 1-13 (2014).
- [44] G. Mu, Z. Qin, and R. Grimshaw, 'Dynamics of rogue waves on a multi-soliton background in a vector nonlinear Schr¨dinger o equation', SIAM J. Appl. Math. 75, 1-20 (2015).
- [45] G. Mu and Z. Qin, 'Dynamic patterns of high-order rogue waves for Sasa-Satsuma equation', Nonlinear Anal. Real World Appl. 31, 179-209 (2016).
- [46] L. M. Ling, 'The algebraic representation for high order solution of Sasa-Satsuma equation', Discrete Continuous Dyn Syst Ser B 9, 1975-2010 (2016).
- [47] L. M. Ling, B. F. Feng and Z. N. Zhu, 'Multi-soliton, multi-breather and higher order rogue wave solutions to the complex short pulse equation,' Physica D 327, 13-29 (2016).
- [48] A. Degasperis. and S. Lombardo, 'Integrability in Action: Solitons, Instability and Rogue Waves', In: M. Onorato, S. Resitori and F. Baronio, Rogue and Shock Waves in Nonlinear Dispersive Media , Lecture Notes in Physics, Vol. 926, pp. 23-53, (Springer, Cham, 2016).
- [49] X. Zhang and Y. Chen 'General high-order rogue waves to nonlinear Schr¨dinger-Boussinesq equation with the dynamical o analysis,' Nonlinear Dynamics 93, 2169-2184 (2018).
- [50] J. Chen, Y. Chen, B. F. Feng, K. Maruno and Y. Ohta, 'General high-order rogue waves of the (1+1)-dimensional Yajima-Oikawa system,' J. Phys. Soc. Jpn. 87, 094007 (2018).
- [51] B. Sun and Z. Lian, 'Rogue waves in the multicomponent Mel'nikov system and multicomponent Schr¨dinger-Boussinesq o system,' Pramana 90, 23 (2018).
- [52] B. Yang and J. Yang, 'Rogue waves in the nonlocal PT -symmetric nonlinear Schr¨dinger equation,' Lett. Math. Phys. o

109, 945-973 (2019).

- [53] B. Yang and J. Yang, 'On general rogue waves in the parity-time-symmetric nonlinear Schr¨dinger o equation', arXiv:1903.06203 (2019).
- [54] A. Chabchoub, N. P. Hoffmann and N. Akhmediev, 'Rogue wave observation in a water wave tank,' Phys. Rev. Lett. 106, 204502 (2011).
- [55] A. Chabchoub, N. Hoffmann, M. Onorato, A. Slunyaev, A. Sergeeva, E. Pelinovsky and N. Akhmediev, 'Observation of a hierarchy of up to fifth-order rogue waves in a water tank,' Phys. Rev. E 86, 056601 (2012).
- [56] B. Kibler, J. Fatome, C. Finot, G. Millot, F. Dias, G. Genty, N. Akhmediev, J. M. Dudley, 'The Peregrine soliton in nonlinear fibre optics,' Nat. Phys. 6, 790-795 (2010).
- [57] B. Frisquet, B. Kibler, P. Morin, F. Baronio, M. Conforti, G. Millot and S. Wabnitz, 'Optical dark rogue wave', Sci. Rep. 6, 20785 (2016).
- [58] F. Baronio, B. Frisquet, S. Chen, G. Millot, S. Wabnitz and B. Kibler, 'Observation of a group of dark rogue waves in a telecommunication optical fiber', Phys. Rev. A 97, 013852 (2018).
- [59] M. Tajiri and Y. Murakami, 'Rational growing mode: exact solutions to the Boussinesq equation', J. Phys. Soc. Japan 60, 2791-2792 (1991).
- [60] J. Rao, Y. Liu, C. Qian and J. He, 'Rogue waves and hybrid solutions of the Boussinesq equation,' Z. Naturforsch. A 72, 307-314 (2017).
- [61] M. J. Ablowitz, and J. Satsuma, 'Solitons and rational solutions of nonlinear evolution equations,' J. Math. Phys. 19, 2180-2186 (1978).
- [62] R. Hirota, The direct method in soliton theory (Cambridge University Press, Cambridge, 2004).
- [63] M. Jimbo and T. Miwa, 'Solitons and Infinite Dimensional Lie Algebras', Publ RIMS, Kyoto Univ. 19, 943-1001 (1983).
Bo Yang Jianke Yang
General rogue waves in the Boussinesq equation
We derive general rogue wave solutions of arbitrary orders in the Boussinesq
equation by the bilinear Kadomtsev-Petviashvili (KP) reduction method. These
rogue solutions are given as Gram determinants with $2N-2$ free irreducible
real parameters, where $N$ is the order of the rogue wave. Tuning these free
parameters, rogue waves of various patterns are obtained, many of which have
not been seen before. Compared to rogue waves in other integrable equations, a
new feature of rogue waves in the Boussinesq equation is that the rogue wave of
maximum amplitude at each order is generally asymmetric in space. On the
technical aspect, our contribution to the bilinear KP-reduction method for
rogue waves is a new judicious choice of differential operators in the
procedure, which drastically simplifies the dimension reduction calculation as
well as the analytical expressions of rogue wave solutions.