{'fraction_non_alphanumeric': 0.03814432989690722, 'fraction_numerical': 0.013092783505154639, 'mean_word_length': 4.2223418573351275, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 6, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Two-breather solutions for the class I infinitely extended nonlinear Schrodinger equation and their special cases

We derive the two-breather solution of the class I infinitely extended nonlinear Schrodinger equation (NLSE). We present a general form of this multi-parameter solution that includes infinitely many free parameters of the equation and free parameters of the two breather components. Particular cases of this solution include rogue wave triplets, and special cases of breather-to-soliton and rogue wave-to-soliton transformations. The presence of many parameters in the solution allows one to describe wave propagation problems with higher accuracy than with the use of the basic NLSE.

Adding higher-order terms generally results in the loss of integrability of the resulting equation. This means that exact solutions cannot be written in analytical form, making the treatment more complicated. However, a special choice of the higher-order operators in these extensions allows us to keep the integrability. The power of using these operators consists in the possibility of applying arbitrary real coefficients to each of these operators, thus significantly extending the range of physical problems that can be solved in exact form. It was found that the NLSE can be extended to arbitrarily high orders of these operators [36,37,38], and these operators have been explicitly presented up to eighth order [37]. Using their recurrence relations, they can be calculated to any order, although the explicit form quickly becomes cumbersome. Nevertheless, there are no conceptual difficulties in construction of these equations. Moreover, infinitely many terms can be considered when finding solutions of these equations.
Presently, there are two sets of these operators that can be used for infinite-order extensions of the NLSE. We call them the class I [36,37,38] and class II [39,40] infinite extensions of the NLSE [41]. The presence of two independent extensions enables the more accurate description of physical problems with greater flexibility. Here, we deal exclusively with the class I extension. The class II extension is more involved, and will be left beyond the scope of the present work.
In this paper, we find 2-breather solutions of the class I infinitely extended NLSE equation. These are multi-parameter solutions that involve both the free parameters of the equation, and free parameters of the solution, which together control the features of the two breather components, such as their localisation, propagation, and their relative position and frequencies. The presence of an infinite number of free parameters allows us to consider many particular cases, such as breatherto-soliton conversion, which is exclusive to higher-order extensions of the basic equation.
We also derive several limiting cases, the most important one of which is the general second-order rogue wave solution, a particular case of the 2-breather collision. However, only a limited number of special cases can be given in the frame of a single manuscript. We leave others for future work in this direction.

The class I infinitely extended NLSE
First, we give a brief exposition of the class I infinitely extended nonlinear Schrödinger equation. It is the integrable equation written in general form [36,37] where the operator F (ψ, ψ * ) is defined through with the operators K n defined recursively by the integrals of the nonlinear Schrödinger equation [36], and where each coefficient α n is an arbitrary real number; that is, where p n is the n-th integral of the basic nonlinear Schrödinger equation, and p n+1 can be defined recursively as The four lowest order operators K n (n = 2, 3, 4, 5) derived in this way are: A few others can be found in [38]. The operators K n involve linear terms with derivatives of order n, and nonlinear terms involving t-derivatives of the function ψ and its complex conjugate ψ * . As already mentioned, the numbers α n can take any values whatsoever, and do not need to be viewed as representing small perturbations for the equation (1) to be completely integrable. This allows us to find solutions for which any order of dispersion can be taken into account without the need for approximation or numerical techniques. This extension substantially widens the range of applicability of the NLSE for solving nonlinear wave evolution problems.
Again, the coefficients α n are finite and arbitrary. However, physical applications, in general, require dispersive effects to decrease rapidly in strength with increasing order n. Convergence will thus not be an issue in practice for series involving α n , and we will therefore be comfortable leaving the operator F for the whole equation (1), as well as any other associated parameters, in the form of an infinite series when necessary.
While the operators K n in (3) rapidly become more complicated and the resulting differential equation of order n becomes much harder to solve, exact solutions can be found explicitly by using already known solutions to the NLSE as a guide, and a large class of breather and soliton solutions are already known [37,38]. In previous works [37,38], we have seen that the effect of nonzero odd order operators is to transform t as t → t + vx with v being a function of all coefficients α 2n+1 . The effect of the nonzero even order operators is to transform x as x → Bx, with B being a function of the parameters α 2n .
In this work we extend this approach to a general family of second-order solutions, so we introduce parameters B 1 and B 2 , and v 1 and v 2 , to play an analogous role for the two distinct breather components. This enables us to generalise the 2-breather to the infinite extension of the NLSE, and we now proceed to the analysis of these solutions.

The 2-breather solution
Higher analogues of the Akhmediev breathers can be obtained through iterations of the Darboux transformation [45,46]. After transforming the plane wave solution e ix with a Darboux transformation, with an eigenvalue λ such that λ 2 = −1, and repeating this transformation twice, we get the 2-breather solution to the basic NLSE. This can then be generalised to the 2-breather solution of the extended equation. The general 2-breather solution is of the form where Here κ 1 and κ 2 are the modulation parameters, is the growth rate of the modulational instability for each breather component, and the shorthand notation t m indicates t m = t + v m x for m = 1, 2. Note that whenever t m appears, we have ignored a constant of integration, and we have also done the same whenever δ m B m x appears. The most general solution allows for the replacements where T m and X m are real constants which determine relative positions along the axes of t and x, respectively, which we might include to incorporate a time delay in one breather component, for instance. For the time being, we set these constants to be both zero without substantial loss, to address their significance later. The phase factor φ is independent of the modulation, since this part has no physical effect on the modulation when it is real, and here it takes the same real value as it does for the plane wave solution, i.e.
The values B m determine the modulation frequency of each component, and the parameters v m , although they cannot be considered velocities in the usual sense, introduce a tilt to |ψ| relative to the axes of x and t. They are given explicitly by with m = 1, 2, where F (a, b; c; z) is the Gaussian hypergeometric function. Note that there is a simple relationship between v m and B m : The first two terms of B m for the two-breather solution have been previously given in [47]. Our new solution extends these coefficients to arbitrary orders of dispersion and nonlinearity. Also notice that the parameters are the same functions of κ m , for both m = 1, 2. This is at least suggested by symmetry. If two successive Darboux transformations generate a 2-breather solution, then there must be two independent eigenvalues, corresponding to two independent modulation parameters. Physically, we could reason that there should be no way of knowing which breather component is which, so the order in which each component was generated by Darboux transformation should be equally irrelevant. If so, it should then follow that B 1 is the same function of κ 1 as B 2 is of κ 2 , and similar for v 1 and v 2 , and this would also imply that B m and v m are the same functions as for the single-breather solution, which are already known [38]. It is worth considering whether this property extends to the general n-breather solution: i.e. whether, in general, we can find B 1 , . . . , B n and v 1 , . . . , v n which are the same functions of their respective modulation parameters κ 1 , . . . , κ n , but we do not answer this question here.
The growth rate δ m in both components will be real when κ m is real, but the eigenvalues of the Darboux transformation are free to take any complex value at all, although the transformations are trivial when the eigenvalues are real, and thus so are the modulation parameters. Real-valued modulation parameters correspond to Akhmediev breathers, whereas imaginary-valued modulation parameters correspond to Kuznetsov-Ma solitons, the functional form of the breathers being otherwise equivalent. An example which shows the difference between real and imaginary modulation parameters is given in Fig. (2). In Fig. (3) we give an example of the effects of altering the ratio of the modulation of the two components.

Breather-to-soliton conversion
If we choose parameters α n such that B m = 0, the 2breather solution may then behave in a way which is unique to the extension of the nonlinear Schrödinger equation [47], in the sense that it is only when higher orders of dispersion and nonlinearity are accounted for that it is possible to take B m = 0 without obtaining a trivial or otherwise degenerate solution.
For example, if we choose α 2 such that B 2 = 0 for all κ 2 , then writing B 1 = B, it is easy to show that B must take the value where we define the function E as the difference of hypergeometric functions: We can then simplify the general 2-breather solution considerably. We obtain An example of this solution is given in Fig. 4. The difference of this solution from the one shown in Fig. 1 is that the wave profiles at x → ±∞ are not plane waves. The periodic set of tails from each breather maximum extends to infinity, reminiscent of periodically repeating solitons. This is the phenomenon that is known as breather-to-soliton conversion [47]. Clearly, these 'solitons' do not have a separate spectral parameter related to them. Fig. 4: A wave profile of a 'breather-to-soliton conversion'. We use the same set of parameters as in Fig. 3, except α 2 is now chosen such that B 2 = 0. This choice extends to infinity the tails of the breathers that would otherwise decay.

The 2-breather solution in the semirational limit
When one of the modulation parameters, say κ 2 , tends to zero, we obtain the semirational limit, i.e. a solution obtained as a combination of polynomials and circular or hyperbolic functions. Then, writing κ for κ 1 , and δ for δ 1 , the functions G, H, and D become and the parameters B 2 and v 2 are reduced to and This semirational 2-breather solution is a superposition of a Peregrine solution with the Akhmediev breather, since taking the limit κ 2 → 0 reduces the frequency of one of the breathers to zero, meaning that it is transformed to a Peregrine solution. A plot of this solution is shown in Fig. 5. Here, the central feature is roughly the second order rogue wave while the peaks away from the origin belong to the remaining first-order breather. Here the nonzero modulation parameter is κ = 1, with α n the same as in Fig. (1). It can be considered as a superposition of the Akhmediev breather with the Peregrine solution.
6 The degenerate two-breather limit If both eigenvalues of the Darboux transformation are taken to be equal, so that both modulation parameters κ m are also equal, we obtain the case of degenerate breathers. Direct calculations provide no solution. In this case, one modulation parameter should instead be taken as a small perturbation from the other, say |κ 1 − κ 2 | = ε. Then, we take the limit as the perturbation ε becomes arbitrarily small, so that the solution remains well-defined at all times. Namely, if we put κ 1 = κ, and κ 2 = κ + ε, we have Next, recalling the identity take the Maclaurin series of the G(x, t), H(x, t), and D(x, t) with respect to ε. In the limit as ε → 0, the ratio of these series will be a well-defined solution with equal eigenvalues; it is thus sufficient to consider only the lowest-order non-vanishing terms in the series expansion for D(x, t) in ε, which in this case happen to be the coefficients of ε 2 . By this method we obtain the degenerate 2-breather solution in the form (6) with where B = B 1 , and we use B and v to denote the partial derivatives of B 2 and v 2 with respect to ε evaluated at the point ε = 0, i.e. when κ 2 → κ. That is, etc. We drop the subscripts due to the fact that as ε → 0 both modulation parameters take equal values anyway. A plot of this solution is given in Fig.6. Fig. 6: The degenerate 2-breather solution. We take the set of α n the same as in Fig. 1, and the modulation parameters κ 1 = κ 2 = 1 2 . The two breathers collide with the high peak at the origin due to the synchronised phases.
The degenerate breather solution is a one-parameter family of solutions which represents the collision of two breathers with the same modulation parameter κ, or, equivalently, with equal frequencies. It can be considered a generalisation of the known 2-soliton solution for the class I extension of the nonlinear Schrödinger equation [50].

Second-order rogue wave solution
When the frequency of the degenerate breather tends to zero, the spacing between the successive peaks in Fig.  6 becomes infinitely large, pushing them out to infinity.
What remains at the origin is the second-order rogue wave. In order to derive this solution, we take the limit κ → 0 in the expressions (17). However, calculations show that this limit cannot be found directly. In order to find it, we repeatedly apply l'Hôpital's rule to the degenerate breather solution as κ → 0. The derivatives of G, H, and D with respect to κ at the point κ = 0 vanish up to O(ε 6 ). The resulting functions G, H, and D become polynomials: where, in the same limit as κ → 0, The first-order derivatives B and v vanish as κ → 0, but the second-order derivatives still remain, and in the limit as κ → 0 are This solution is shown in Fig. 7. It is, naturally, the second-order rogue wave, but slanted and rescaled in the (x, t)-plane relative to the second-order rogue wave of the NLSE [51,52]. Fig. 7: The second-order rogue wave, Eqs. (18), (19), (20) obtained from the degenerate 2-breather solution shown in Fig. 6 in the limit κ → 0, with some stretching due to higher-order effects.

Rogue wave triplets
It is well known that the general n-th order rogue wave has the remarkable property of being able to split into 1 2 n(n+1) first-order components [48]. The second-order rogue wave discussed above is only a particular case of a more general rogue wave structure, where all three first-order components are located at the origin, and have merged into one single peak. In order to obtain the more general solution where the three components are not merged together, known as the rogue wave triplet [49], we re-introduce the constants of integration into the general 2-breather solution, i.e.
where X m and T m are arbitrary, and the parameter ε is introduced to make sure that the Taylor series in the degenerate limit still vanishes up to O(ε 2 ). The values of X m and T m determine the location of the components of the breather components. They add additional free parameters to the solution which we have previously given for the restricted case in which X m = T m = 0. Notice also that we do not make the replacement x → x − X m ε directly, but, for simplicity's sake, instead define X m to account for the higher-order terms in B m . In order to further simplify parametrisation, we assume that X m and T m are functions of the modulation parameter κ, and are of the order O(κ). Then, defining free parameters ξ and η independent of κ, such that κξ = 48(X 1 − X 2 ), κη = 48(T 1 − T 2 ), we have in the limit as κ → 0 the rogue wave triplet solution in the form witĥ whereĜ,Ĥ andD now contain two new free parameters, ξ and η, which determine the separation of the fundamental rogue wave components in the triplet [49], and where G, H, and D are as given in Eqs. (18)-(20), for the particular case in which ξ = η = 0. An example of the formation of rogue wave triplets, corresponding to nonzero ξ and η, is shown in Fig. 8. When both ξ = 0 and η = 0, all three components merge at the origin, as in Fig. 7. Generally, the coefficient B in Eq. (21) determines the degree of localisation along the x-axis. Larger values of B will correspond to narrower peaks, whereas smaller values of B will correspond to broader peaks, and B = 0 to minimal localisation. A point of interest here is that it is again possible to choose a parametrisation for which B is any fixed constant. If we choose, for instance, we end up with B = c, where c is a free parameter. However, B is entirely independent of the choice of c, since the coefficient of α 2 in B is zero. As the simplest example, we consider the completely de-localised case, B = 0, with B remaining arbitrary. The rogue wave solution then reduces to (21) witĥ H(x, t) = H 0 (x, t) + 12ξ + 48ξ(t + vx) 2 , D(x, t) = D 0 (x, t) − (ξ 2 + η 2 ) − 48(ξB + ηv )x+ + 12η(t + vx) − 16η(t + vx) 3 .
where G 0 (x, t) = 12{−3 − 192v x(t + vx) + 24(t + vx) 2 + + 16(t + vx) 4 }, Here, G 0 , H 0 , and D 0 are as given for the case where the components are merged and B = 0, andĜ,Ĥ,D incorporate the shifting of the first-order components through ξ and η. When B = 0, the second-order rogue wave acquires soliton-like tails similar to those in Fig. 4. When, additionally, ξ = η = 0, rogue waves merge at the origin to form a second-order rogue wave with extended tails. This case is shown in Fig. 9. When ξ or η is not zero, the components split, resulting in the disappearance of the central peak. This case is shown in Fig. 10. Here, the central peak is absent but the long tails remain, consisting of de-localised first-order components.

Conclusions
We have derived the general 2-breather solution for the class I infinitely extended nonlinear Schrödinger equation, and given many limiting cases; namely, breatherto-soliton conversions, the semirational limit, the degenerate 2-breather, and, probably most importantly, the general second-order rogue wave solution. These solutions completely describe a large family of secondorder solutions to the class I extension of the NLSE, Fig. 9: The second-order rogue wave solution with 'soliton'-like tails when α 2 chosen such that B = 0, and ξ = η = 0. Other parameters are the same as in Fig. 8. and exhibit rich behaviour.
These results provide a more detailed analysis of the formation of nonlinear wave structures such as breathers and rogue waves when higher-order effects come into play, and leaves a large range of future related work wide open.