Access the full text.
Sign up today, get DeepDyve free for 14 days.
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
We consider the nonlinear Helmholtz (NLH) equation describing the beam propagation in a planar waveguide with Kerr-like nonlinearity under non-paraxial approximation. By applying the Lie symmetry analysis, we determine the Lie point symmetries and the corresponding symmetry reductions in the form of ordinary differential equations (ODEs) with the help of the optimal systems of one-dimensional subalgebras. Our investigation reveals an important fact that in spite of the original NLH equation being non-integrable, its symmetry reductions are of Painleve´ inte- grable. We study the resulting sets of nonlinear ODEs analytically either by constructing the integrals of motion using the modiﬁed Prelle-Singer method or by obtaining explicit travelling wave-like solutions including solitary and sym- biotic solitary wave solutions. Also, we carry out a detailed numerical analysis of the reduced equations and obtain multi-peak nonlinear wave trains. As a special case of the NLH equation, we also make a comparison between the symmetries of the present NLH system and that of the standard nonlinear Schro¨dinger equation for which symmetries are long available in the literature. Keywords: Lie symmetry analysis, nonlinear Helmholtz equation, symmetry reduction, Painleve´ analysis, modiﬁed Prelle-Singer method, periodic and solitary waves 1. Introduction Investigation of physical systems described by nonlinear evolution equations (NLEEs) and exploring their underly- ing dynamics remain the central focus of research for the past few decades. Finding exact solutions of these nonlinear equations, which can be either ordinary differential equations (ODEs) or partial differential equations (PDEs), is one of the most important tasks and their further investigation plays a crucial role in the study of nonlinear physical phe- nomena. Lie symmetry analysis [1–3] has a time honoured history and has been proved to be a powerful tool for studying nonlinear problems arising in many scientiﬁc ﬁelds including mathematics, physics, and biology. A Lie (point) symmetry of a given system of NLEE(s) is nothing but an inﬁnitesimal transformation of all of its (their) independent and dependent variables, which leaves the corresponding nonlinear equations invariant under that trans- formation. From the Lie point symmetries, with the aid of the invariance conditions, a given system of equations is reduced to a set of ODEs/PDEs with lesser number of independent coordinates. The reduced equation becomes an (or a set of) ODE(s) if the independent variables involved in the original NLEE are two, while it becomes a (or a set of) PDE(s) when the NLEE contains more than two independent variables. Thus a similarity reduction of a differential equation is closely connected with the invariance of the respective equation. For a detailed background on the theory and application of the Lie symmetry methods one can refer to Refs. [1–9]. Corresponding author. Tel.: +91 431 2407093; Fax: +91 431 2407093 Email addresses: firstname.lastname@example.org (K. Sakkaravarthi), email@example.com (A.G. Johnpillai), firstname.lastname@example.org (T. Kanna), email@example.com; firstname.lastname@example.org (M. Lakshmanan) Preprint submitted to Applied Mathematics and Computation March 6, 2018 arXiv:1803.01622v1 [nlin.SI] 5 Mar 2018 This work deals with the Lie symmetry analysis and the construction of group invariant solutions of a physically important evolution equation, namely the following nonlinear Helmholtz (NLH) equation iq + k q + q + γ |q| q = 0, (1) t tt xx where q is a complex valued dependent variable representing the envelope ﬁeld, the independent variables t and x denote the longitudinal and transverse coordinates, respectively. Here k and γ are arbitrary real constants and represent the coefﬁcients of non-paraxial term and the cubic nonlinearity, respectively. Equation (1) is receiving much attention in different areas of physics and mathematics. In the context of nonlinear optics, Eq. (1) describes the ultra-broad beam propagation in a Kerr-like nonlinear medium under the non-paraxial approximation . Due to this reason, Eq. (1) can also be termed as non-paraxial nonlinear Schro¨dinger equation . Physically, Eq. (1) can be viewed as an optical system governing the propagation of beam in a Kerr type nonlinear medium experiencing diffraction in both the transverse and longitudinal directions. The NLH equation (1) seems to be non-integrable as it does not posseses a Lax pair. It also does not admit sufﬁcient number of conserved quantities. So far to the best of our knowledge, there is no claim for the integrability of the NLH equation (1) but it admits interesting special solutions. Depending upon the sign of the nonlinearity coefﬁcient γ, Eq. (1) reduces to a focusing-type NLH equation (for γ > 0) or a defocusing-type NLH equation (for γ < 0), for which respectively bright and dark solitary-wave solutions have been reported [12, 13]. In recent years the multicomponent NLH system has also received considerable interest. Special solitary waves and elliptic waves for a multicomponent NLH system exhibiting distinct dynamical behaviours have been reported in . On the other hand, under the paraxial approximation the Maxwell’s equations result in the following standard nonlinear Schro¨dinger (NLS) equation : iq + q + γ |q| q = 0. (2) t xx Mathematically, NLS equation (2) can be obtained from the NLH equation (1) as a special case, for k = 0. Physically, it can be derived to describe a beam propagating in a Kerr type nonlinear medium, where the beam radius is sufﬁciently large compared with the beam wavelength which supports only the paraxial approximation . Here the second order derivative of the ﬁeld variable with respect to propagation direction (q ) vanishes due to the exclusion of non- tt paraxial effect (k = 0). The complete integrability of NLS equation (2) has been estabilished via the inverse scattering transformation method, the Painleve´ analysis, admittance of inﬁnite number of conserved quantities, etc. . Also, the NLS equation admits multi-soliton solutions (bright solitons for γ > 0 and dark solitons for γ < 0) in addition to several other exact analytical solutions and display elastic interaction of solitons . The objectives of the present work are three-fold. First and the important objective is to obtain the symmetry reductions of the NLH system (1) using the inﬁnitesimal similarity transformations. These similarity transformations are constructed by utilizing the Lie point symmetry generators admitted by the complex PDE (1). In addition to the construction of symmetry generators, we obtain a set of associated nonlinear ODEs as the integrable reductions of the NLH system (1). Second, we investigate the integrable nature of the reduced nonlinear ODEs by applying the Painleve´ singularity structure analysis. Then, we derive the ﬁrst integrals/constants of motion of these nonlinear ODEs by adopting the modiﬁed Prelle-Singer method and study the reduced ODEs by means of direct numerical analysis. The third objective is to revisit the symmetry analysis of the NLS equation (2) [8, 9] and to compare the results with that of the present NLH equation (1). The outline of the paper is as follows. In Section 2, we brieﬂy discuss the general algorithm for the Lie point symmetry analysis. We present the symmetry reductions and group-invariant solutions of the NLH equation (1) in Sec. 3. Section 4 deals the symmetry reductions of the NLS equation (2) and their relevance to the NLH system. Finally, in Section 5, we summarize the main results of the present work. 2. General algorithm for Lie symmetries and invariants In this section, we brieﬂy explain the important steps involved in the determination of Lie symmetries for a given system of (complex or real) nonlinear PDEs [1–9]. Firstly, the given equation has to be rewritten as a function of 2 all dependent and independent variables by using a similarity transformation with an inﬁnitesimal parameter. For example, the nonlinear equation (having N dependent functions u and two independent functions t and x) of the form F (u , u , u , ..., u , x, t) = 0, j = 1, 2, 3, ..., N, (3) j 1 2 3 N can be transformed to a set of PDEs using the following inﬁnitesimal transformations: j 2 u → uˆ = u + ǫ η (x, t, u , u , u , ..., u ) + O(ǫ ), j = 1, 2, 3, · · · , N, (4a) j j j 1 2 3 N t → t = t + ǫ τ(x, t, u , u , u , ..., u ) + O(ǫ ), (4b) 1 2 3 N x → xˆ = x + ǫ ξ(x, t, u , u , u , ..., u ) + O(ǫ ), (4c) 1 2 3 N where ǫ is a small expansion parameter. According to the Lie symmetry algorithm, a vector ﬁeld X = τ(t, x, u , u , ..., u )∂ + ξ(t, x, u , u , ..., u )∂ + η (t, x, u , u , ..., u )∂ , (5) 1 2 N t 1 2 N x 1 2 N u j=1 is a generator of point symmetry of the equation (3) if  X F (u , u , u , ..., u , x, t) = 0. (6) j 1 2 3 N F =0  Here the operator X is the second prolongation of the operator X which is determined by the order of the given differential equation and is deﬁned as N N N X X X j j j  X = X + ζ ∂ + ζ ∂ + ζ ∂ , (7) u u u x xx t j,t j,x j,xx j=1 j=1 j=1 j j j j where the coefﬁcients ζ , ζ , and ζ correspond to the prolongation formulae, ζ = D (η ) − u D (τ) − u D (ξ), x xx t j,t t j,x t t t j j j ζ = D (η )− u D (τ)− u D (ξ), ζ = D (ζ )− u D (τ)− u D (ξ). Here the index ‘ j’ ( j = 1, 2, 3, ..., N) denotes x x j,t x j,x x xx x x j,xt x j,xx x the dependent variable and the subscripts represent the independent variables t and x, while D and D are the standard t x total derivative operators. One can obtain the explicit form of the inﬁnitesimal coefﬁcients τ, ξ and η by solving the determining equations (6) with required number of arbitrary constants through tedious but straightforward calculations. From these inﬁnites- imal coefﬁcients and Eq. (5), the Lie point symmetry generators (vector ﬁelds) admitted by the system of PDEs (3) corresponding to each arbitrary constants can be furnished in a direct way. The symmetry reductions of the given equations can be identiﬁed by solving the following characteristic equation dt dx du du du du 1 2 3 N = = = = = ... = . (8) 1 2 3 N τ ξ η η η η Here one of the important points to note is that we wish to minimize the search for invariant solutions by ﬁnding the nonequivalent branches of solutions resulting for different combinations (linear superposition) of symmetry genera- tors. For this purpose, an optimal system of one-dimensional subalgebras is constructed with the help of adjoint repre- sentations of the generators and this gives exact invariants as well as the associated symmetry reductions. Speciﬁcally, the vector ﬁelds (X ) are written as a superposition, with equal number of arbitrary coefﬁcients (c ), in the form X = j j c X +c X +c X +· · ·+c X ≡ c X , where j = 1, 2, 3, · · · , n, and n is the number of available vector ﬁelds. Here 1 1 2 2 3 3 n n j j j=1 one has to reduce this equation to the simplest forms for different choices of the coefﬁcients c by using the adjoint 2 3 ε ε characteristics of the vector ﬁelds deﬁned as Ad(εX , X ) = X − ε[X , X ]+ [X , [X , X ]] − [X , [X , [X , X ]]]+ . . . , i j j i j i i j i i i j 2! 3! where the square bracket [X , X ] represents the usual commutator. For more details regarding the one-dimensional i j subalgebras, see [3, 17] and references therein. Especially, such type of group classiﬁcation for certain systems of nonlinear reaction diffusion equations have been investigated in . Thus for each optimal subset of the generators we obtain appropriate invariants which reduce the original PDEs to a set of ODEs. Then we have to analyze these reductions and solve them to ﬁnd their explicit solutions which can reveal the real dynamics of the physical problem of interest. By following the above algorithm, in the following sections, we perform the Lie point symmetry analysis of the NLH equation and also compare the results with those of the NLS equation. 3 3. Lie symmetry analysis of the nonlinear Helmholtz equation (1) In this section, we systematically derive the Lie point symmetries of the NLH equation (1). For this purpose, we denote q(x, t) = u(x, t) + iv(x, t), and decompose equation (1) into real and imaginary parts to obtain the following system of partial differential equations (PDEs): 2 2 u + kv + v + γ v(u + v ) = 0 ≡ F , (9a) t tt xx 1 2 2 v − ku − u − γ u(u + v ) = 0 ≡ F . (9b) t tt xx 2 From the previous section, we write the symmetry transformations as u → u = u + ǫ η (x, t, u, v), v → u = 1 2 v + ǫ η (x, t, u, v), t → t = t + ǫ τ(x, t, u, v), and x → xˆ = x + ǫ ξ(x, t, u, v). Then the generalized vector ﬁeld associated with Eq. (9) can be written as 1 2 X = τ(t, x, u, v)∂ + ξ(t, x, u, v)∂ + η (t, x, u, v)∂ + η (t, x, u, v)∂ . (10) t x u v The Lie symmetry/ invariance conditions of the PDEs (9) read h i  2 2 X u + kv + v + γ v(u + v ) = 0, (11a) t tt xx F =0 h i  2 2 X v − ku − u − γ u(u + v ) = 0, (11b) t tt xx F =0  where X denotes the second prolongation operator and it takes the form (6) with N = 2. 1 2 We have obtained the coefﬁcient functions τ, ξ, η and η by solving the determining equations resulting from 1 2 (11). Particularly, we ﬁnd that τ, ξ, η and η are independent of the derivatives of u and v. So, the coefﬁcients of like derivatives of u and v in Eq. (11) are equated to yield an over determined system of linear PDEs. Solving these 1 2 equations recursively, we have obtained the form of the inﬁnitesimal coefﬁcients τ, ξ, η and η of the NLH equation (9) as τ = c + 2kc x, ξ = c − 2c t, (12a) 1 4 2 4 1 2 η = (c x − c )v, η = −(c x − c )u, (12b) 4 3 4 3 where c , c , c and c are arbitrary real constants. Thus, the explicit expression of the generalized vector ﬁeld (10) 1 2 3 4 can be written as X = (c + 2kc x)∂ + (c − 2c t)∂ + (c x − c )v∂ − (c x − c )u(t, x, u, v)∂ . (13) 1 4 t 2 4 x 4 3 u 4 3 v Also, from the explicit form of the inﬁnitesimal coefﬁcients (12), the Lie point symmetry generators for the four arbitrary constants admitted by the NLH system of PDEs (9) are obtained as below: X = ∂ , X = ∂ , X = −v ∂ + u∂ , (14a) 1 t 2 x 3 u v X = 2kx ∂ − 2t ∂ + xv ∂ − xu ∂ . (14b) 4 t x u v Note that the ﬁrst three Lie symmetry generators X , X , and X can be associated with translation in time, translation 1 2 3 in space, and phase transformations, respectively. Now, we utilize the aforementioned Lie point symmetry generators (14) of the NLH equation (1) to obtain its sym- metry reductions. Then we construct exact group-invariant solutions by reducing Eq. (9) into a set of ODEs. In this regard, we ﬁrst construct the optimal system of one-dimensional subalgebras from the generators (14). Here the com- mutator and adjoint property corresponding to each generator set plays a crucial role in determining the subalgebras. We have given such commutator and adjoint properties of the generators in Table 1 and Table 2, respectively. Table 1. Commutator table of the vector ﬁelds of Eq. (9) [X , X ] X X X X i j 1 2 3 4 X 0 0 0 −2X 1 2 X 0 0 0 2kX − X 2 1 3 X 0 0 0 0 X 2X X − 2kX 0 0 4 2 3 1 4 Table 2. Adjoint table of the vector ﬁelds of Eq. (9) Ad(εX , X ) X X X X i j 1 2 3 4 X X X X X + 2εX 1 1 2 3 4 2 X X X X X − ε(2kX − X ) 2 1 2 3 4 1 3 X X X X X 3 1 2 3 4 2 2 X X − 2εX + ε (X − 2kX ) X − ε(X − 2kX ) − 2ε kX X X 4 1 2 3 1 2 3 1 2 3 4 After a careful analysis, we ﬁnd that the present NLH system admits the following optimal set of one-dimensional subalgebras as a linear superposition of symmetry generators (14): (i) X + aX , (ii) X + aX , (iii) X + aX + bX 1 2 2 3 1 2 3 and (iv) X , where a and b are arbitrary real constants. By solving the resulting set of characteristic equations for the above subalgebras, we get the invariants and group-invariant solutions for Eq. (9) as given in the Table 3. Table 3. Subalgebra, group invariants, and group-invariant solutions of Eq. (9) Case Subalgebra Invariant ‘y’ Invariant solution (i) X + aX x − at u = A(y), v = B(y) 1 2 (ii) X + aX t u = −A(y) sin ax + B(y) cos ax 2 3 v = A(y) cos ax + B(y) sin ax (iii) X + aX + bX x − at u = −A(y) sin bt + B(y) cos bt 1 2 3 v = A(y) cos bt + B(y) sin bt t t 2 2 (iv) X kx + t u = A(y) cos( ) + B(y) sin( ) 2k 2k t t v = −A(y) sin( ) + B(y) cos( ) 2k 2k In what follows, we analyze the invariant reductions obtained in the form of coupled ODEs and the associated integrals for the above four cases. We obtain exact analytical solutions for certain equations in terms of travelling waves while the other equations are studied numerically to explore the nature of the solitons. Case (i): X + aX (Travelling wave solutions) 1 2 In this case, substitution of the group-invariant solutions (case (i) in Table 3) into the PDEs (9) results in the following system of nonlinear second-order ODEs: ′′ 2 ′ 2 2 (1 + ka )A + aB + γA(A + B ) = 0, (15a) ′′ 2 ′ 2 2 (1 + ka )B − aA + γB(A + B ) = 0. (15b) Here and in the following ‘prime’ appearing in A and B represents differentiation with respect to the invariant variable y. In order to know about the integrable nature of the reduced ODEs, we have performed the Painleve´ singularity structure analysis  and ﬁnd that the above coupled second-order ODE is Painleve´ integrable for arbitrary a, k and γ parameters. Details of the Painleve´ analysis to Eq. (15) are given in Appendix A. Next we investigate Eq. (15) by using the modiﬁed Prelle-Singer method [19–21] (a brief algorithm for this method is given in Appendix B), an effective method to obtain the ﬁrst integrals. The explicit forms of the ﬁrst integrals of Eq. (15) obtained by this method are given below: 2 2 2 γ(A + B ) ′2 ′2 I = 2(A + B ) + , (16a) 1 + ka 2 2 a(A + B ) ′ ′ I = (A B − AB ) + , (16b) 2(1 + ka ) where I and I are the two ﬁrst integrals of Eq. (15). As the equation (15) is a system of two coupled second 1 2 order autonomous ODEs, the existence of two functionally independent ‘time-independent’ integrals itself ensures the complete integrability of the system. One can also proceed to ﬁnd two additional time dependent integrals to deduce the explicit solutions as well. The above integrals given by Eq. (16) can also be viewed as conserved quantities associated with the symmetry generators X + aX . 1 2 5 Figure 1: Nonlinear periodic wave trains of Eq. (15) for different non-paraxial coefﬁcient k = 0, 0.5, 1 (from top to bottom) and γ = 2. The left and right panels correspond to a = 1 and a = −1, respectively. Here and in the ﬁgures follow, the solid-red line represents the solution of ‘A’ and the dashed-blue line denotes the solution of ‘B’. Our numerical analysis of Eq. (15) shows that it admits nonlinear periodic wave structures for A(y) and B(y). By changing the parameters k and a the period of oscillations and amplitudes can be altered as shown in Fig. 1. The top panels show the nonlinear periodic waves in the absence of non-paraxial effect (k = 0), which corresponds to the NLS case, while the middle (k = 0.5) and bottom (k = 1) panels display the nonlinear periodic waves in its presence. A noticeable observation is that due to the non-paraxial effect the width of the pulse increases. Also, the nonlinear wave structures approach a stable proﬁle with less ﬂuctuations in their amplitude, especially when a = −1. One can evidence the applications of such nonlinear periodic waves in ﬁber Bragg gratings. For a special case, k = −1/a , the above system of nonlinear second-order ODEs (15) becomes a set of ﬁrst-order ODEs ′ 2 2 aA − γB(A + B ) = 0, (17a) ′ 2 2 aB + γA(A + B ) = 0. (17b) Solving the system of nonlinear ODEs (17), we get the explicit solution A(y) = 2b sin((−2b γ/a)y + b ) and 1 1 2 B(y) = − 2b cos((−2b γ/a)y + b ), where b and b are arbitrary real constants. Thus a travelling wave solution to 1 1 2 1 2 the equations (9) is given by −2b γ u(x, t) = 2b sin (x − at) + b , (18a) 1 2 −2b γ v(x, t) = − 2b cos (x − at) + b . (18b) 1 2 These solutions are not the standard travelling waves due to the dependence of frequency (2b γ) and wave number (2b γ/a) on the amplitude ( 2b ). Also, we have shown the numerically obtained periodic wave solutions in the form 1 1 6 of sine and cosine functions for Eq. (17) in Fig. 2. From Fig. 2, we ﬁnd that the analytical and numerical solutions of Eq. (17) exactly coincide. This ultimately conﬁrms the validity of direct numerical results of other set of ODEs obtained in this work. Such non-standard periodic waves can also be viewed as the propagation of nonlinear waves in bimodal optical ﬁbers. Also in Fig. 3, we have shown the periodic solution for a different b (b ) parameters, which 1 2 inﬂuence the amplitude and frequency (phase-shift). It also contains the original travelling wave solution for Eq. (9). Figure 2: Periodic sine and cosine wave solutions of Eq. (17) for a = 1 (a,c) and a = −1 (b,d) with γ = 2. Here the top panel shows the analytical solution for b = 0.5 and b = 3, while the bottom panel represents its direct numerical solution. A: solid-red line and B: dashed-blue line. 1 2 Figure 3: Periodic wave solutions of Eq. (17) and the travelling waveform (y = x − at) of Eq. (9) for γ = 2, a = 1, b = 0.5 and b = 1.5. 1 2 Another subcase of case (i) can be obtained for a = 0 as a is an arbitrary constant. This choice reduces the subalgebra as X + aX → X . Here also we obtain the group-invariant solutions as u(x, t) = A(y) and v(x, t) = B(y) 1 4 1 with invariant y = x. Consequently the PDEs (9) result into the following ODEs: ′′ ′ 2 2 kA − B + γA(A + B ) = 0, (19a) ′′ ′ 2 2 kB + A + γB(A + B ) = 0. (19b) The above equation too falls into the category of Painleve´ integrable case and the forms of integrals obtained by using the modiﬁed Prelle-Singer method can be written as below, which follows from (16). 2 2 2 γ(A + B ) ′2 ′2 I = 2(A + B ) + , (20a) 2 2 (A + B ) ′ ′ I = (AB − A B) + . (20b) 2k 7 The numerical solutions of Eq. (19) are shown in Fig. 4 for k = 0.5 and γ = 2. It shows that system (19) admits triple peak nonlinear wave structures of varying amplitudes. Such type of multipeaked stationary solutions exist in laser-BEC interactions and it has been shown that they can also explain the formation of bound state. Also, these triple-peak nonlinear wave can be associated with a series of asymmetric optical multipeak solitons on a continuous wave background . Figure 4: Triple peak nonlinear wave trains of Eq. (19) for k = 0.5 and γ = 2. Case (ii): X + aX 2 3 The next subgroup X + aX in the optimal subalgebras reduces Eq. (9) to the following system of nonlinear 2 3 second-order ODEs ′′ ′ 2 2 2 kA + B − a A + γA(A + B ) = 0, (21a) ′′ ′ 2 2 2 kB − A − a B + γB(A + B ) = 0. (21b) The Painleve´ analysis again shows that Eq. (21) is integrable thereby admitting required number of arbitrary constants for all resonance values. Again by using the modiﬁed Prelle-Singer method we obtain the explicit forms of integrals as below: γ 2a ′2 ′2 2 2 2 2 2 I = (A + B ) + (A + B ) − (A + B ), (22a) k k 2 2 A + B ′ ′ I = (A B − AB ) + . (22b) 2k As in the above case, Eq. (21) is solved numerically and its solutions are shown in Fig. 5. Here we observe a special set of non-identical chirped pulse trains for A and B which repeats periodically. Such type of chirped pulses ﬁnd applications in nonlinear optical systems with higher order nonlinearities like quintic and septic effects. Here also one can observe the inﬂuence of non-paraxial coefﬁcient k, which is directly proportional to the pulse width, that is, the increase in k results in a increase of the pulse width of the nonlinear wave trains (see Fig. 5). Figure 5: Nonlinear wave trains of Eq. (21) for γ = 2, a = 1, (a) k = 0.5 and (b) k = 1.5. One can also obtain a subcase of case (ii) for the choice a = 0, which results into X + aX → X . Here we obtain 2 3 2 the group-invariant solutions as u(x, t) = A(y) and v(x, t) = B(y), where the invariant y = t, which reduces the PDEs 8 (9) into a set of ODEs ′′ 2 2 A + γA(A + B ) = 0, (23a) ′′ 2 2 B + γB(A + B ) = 0. (23b) We have performed the Painleve´ analysis and it shows that the above equation is P-integrable. A basic solution for the above equation can be obtained as A(t) = b sin( γb t + b ), (24a) 1 1 2 B(t) = b cos( γb t + b ), (24b) 1 1 2 where b and b are arbitrary constants. We have shown the above periodic wave solutions in Fig. 6 for two different 1 2 choices of b and b . Even though the amplitude of periodic solutions (24) are same in A and B, they are out of phase 1 2 to each other. In addition to the above periodic solutions, one can also obtain elliptic function solutions of Eq. (23). Figure 6: Periodic sine and cosine wave solutions (24) for (a) b = 1 and b = 0 and (b) b = 1.5 and b = 3, with γ = 2. 1 2 1 2 Case (iii): X + aX + bX 1 2 3 From the Table 3, we ﬁnd that the group-invariant solution resulting for the sub-algebra X + aX + bX gives rise 1 2 3 to the following system of nonlinear second-order ODEs to the NLH equations (9): ′′ 2 ′ 2 2 (1 + ka )A − a(1 + 2kb)B − b(1 + kb)A + γA(A + B ) = 0, (25a) ′′ 2 ′ 2 2 (1 + ka )B + a(1 + 2kb)A − b(1 + kb)B + γB(A + B ) = 0. (25b) Here also we ﬁnd that Eq. (25) passes the Painleve´ test and is integrable for arbitrary values of parameters k, γ, a, and b. We have derived the ﬁrst integrals of Eq. (25) by employing the modiﬁed Prelle-Singer method as 2 2 2 2(1 + kb)b γ(A + B ) 2 2 ′2 ′2 I = (A + B ) − − 2(A + B ), (26a) 2 2 (1 + ka ) (1 + ka ) a(1 + 2kb) ′ ′ 2 2 I = (AB − A B) + (A + B ). (26b) 2(1 + ka ) The numerical analysis of Eq. (25) displays several types of nonlinear wave trains for different choices of param- eters as shown in Fig. 7. The arbitrary system parameters (a, b, k and γ) can be suitably altered to manipulate the nature of the resulting periodic structures. Here the linear coupling coefﬁcient b changes the periodicity of wave trains in the absence of non-paraxial coefﬁcient (k = 0) which corresponds to the NLS case. If k = −1/a , then the system of nonlinear second-order ODEs (25) becomes ′ 2 2 αB + βA − γA(A + B ) = 0, (27a) ′ 2 2 αA − βB + γB(A + B ) = 0, (27b) 9 Figure 7: Numerically obtained nonlinear wave trains in Eq. (25) for (a) k = 0.5 and b = 1, (b) k = 0 and b = 1, and (c) k = 0 and b = −1, with the other parameters γ = 2 and a = 1. 2 2 where α = a − 2b/a and β = b − b /a . A solution to the system of nonlinear ﬁrst-order ODEs (27) is A(y) = √ √ b sin(((β − γb )/α)y + b ) and B(y) = b cos(((β − γb )/α)y + b ), where b and b are arbitrary constants. Thus 1 1 2 1 1 2 1 2 the group-invariant solution to the equations (9) takes the form ! ! β − γb u(x, t) = b cos (x − at) + bt + b , (28a) 1 2 ! ! β − γb v(x, t) = b sin (x − at) + bt + b . (28b) 1 2 We have plotted the above analytical solution of Eq. (27) and the corresponding travelling periodic wave structures of (9) in Fig. 8. 1.0 0.5 A,B 0.0 - 0.5 - 1.0 0 5 10 15 20 Figure 8: Nonlinear periodic wave solution (28) of (27) and its travelling waveform (y = x − at) of Eq. (9) for γ = 2, a = 1, b = 1, b = 1 and b = 1.5. −1 One can also get another special case of Eq. (25) for k = , where the resultant set of ODEs become 2b ′′ 2 2 2 (1 + ka )A − b(1 + kb)A + γA(A + B ) = 0, (29a) ′′ 2 2 2 (1 + ka )B − b(1 + kb)B + γB(A + B ) = 0. (29b) The above Eq. (29) admits a special bright solitary wave solution which can be written in the following form. A = a sech(k y + k ), (30a) 1 1 2 B = a sech(k y + k ), (30b) 2 1 2 b(1+kb) 2b(1+kb) 2 2 2 2 with the constraints k = and a = − a , where we take k > 0. In Fig. 9 we have shown the above 1 1+ka 2 γ 1 1 bright solitary wave solution for a focussing type nonlinearity (γ > 0). Note that the width as well as the central position of the solitary waves alter as we change the parameter ‘a’ without affecting their amplitudes. Interestingly, for particular a choice of parameters, namely 2b(1 + kb) = γa , the solitary wave appears only in the ﬁrst component (A) while it vanishes in the other (B) component, see Fig. 9c. We refer to these structures as symbiotic solitary waves. 10 Figure 9: Solitary wave solutions, including symbiotic one, of Eq. (29) for (a) γ = 2, and a = ±1, (b) γ = 2, and a = 0, and (c) γ = 1, and a = 1, along with b = 1, k = −7. Here we expect that Eq. (29) can also support dark solitary wave for defocussing type nonlinearity (γ < 0). It would be nice to have rigorous study on (29) in future dealing with various solutions. The above solutions (30) display interesting travelling solitary wave structures with y = x − at in the NLH system (9) which we have shown below. In fact we observe beating solitary wave propagation, where the amplitude of solitary wave oscillates periodically. Figure 10: Travelling solitary wave solutions of Eq. (9) for the choice in Fig. 9a with y = x − at. Figure 11: Stationary solitary wave solutions of Eq. (9) for the choice in Fig. 9b with y = x − at. Case (iv): X Finally, for the symmetry generator X and its group-invariant solution, the system of PDEs (9) can be reduced into a system of nonlinear second-order ODEs as ′′ 1 1 γ ′ 2 2 yA + A + A + A(A + B ) = 0, (31a) 2 16k 4k ′′ 1 1 γ ′ 2 2 yB + B + B + B(A + B ) = 0. (31b) 2 16k 4k 11 Our analysis shows that the above set of ODEs (31) is also Painleve´ integrable. By applying the modiﬁed Prelle- Singer method , we obtain the ﬁrst integrals of the above equations (31) as ′2 ′2 2 2 2 2 2 (A + B )y γ(A + B ) (A + B ) I = + + , (32a) 2 16k 32k 1/2 ′ ′ I = y (A B − AB ). (32b) The direct numerical analysis of Eq. (31) is carried out and the results are given in Fig. 12. It shows that due to the explicit appearance of ‘y’ in the model equation (31), the nonlinear wave train expands almost with constant amplitude as y increases. Figure 12: Numerical solutions of Eq. (31) for (a) k = 0.5 and (b) k = 1.0 with γ = 2. In spite of all the above similarity reduced ODEs being integrable, one can show that the original NLH system is non-integrable. More details will be presented elsewhere. 4. Lie symmetry analysis of the nonlinear Schro¨dinger equation (2) As mentioned in the introduction, in this section we brieﬂy outline the results of symmetry reductions of the NLS equation (2) which arises as a special case of the NLH equation (1), for comparison purpose. We wish to note that the symmetry analysis of the NLS system has already been carried out in detail in refs. [8, 9]. By rewriting Eq. (2) into real and imaginary parts, we get a set of PDEs as given in (9) with k = 0. 2 2 u + v + γ v(u + v ) = 0, (33a) t xx 2 2 v − u − γ u(u + v ) = 0. (33b) t xx By applying the inﬁnitesimal symmetry transformations (4c) and recursively solving the equations resulting from (9), τ = c + 2c t, ξ = c + 2c t + c x, (34a) 1 5 2 4 5 1 2 η = −(c + c x)v − c u, η = (c + c x)u − c v. (34b) 3 4 5 3 4 5 Then Eq. (34) leads to the following symmetry generators for the NLS equation (2) X = ∂ , X = ∂ , X = −v ∂ + u∂ , (35a) 1 t 2 x 3 u v X = 2t ∂ − xv ∂ + xu ∂ , (35b) 4 x u v X = 2t ∂ + x ∂ − u ∂ − v ∂ . (35c) 5 t x u v From the above generators we can note that the NLS equation (2) is invariant under translation in time (X ), trans- lation in space (X ), transformation of the phases (X ), Galilean boost (X ), and scaling transformation (X ). From 2 3 4 5 the optimal system of one-dimensional subalgebras, we ﬁnd that the following combinations of generators give group invariant solutions for the NLS equation (2): (i) X + aX , (ii) X + cX , (iii) X + cX , (iv) X + cX + bX (v) X + cX 5 3 1 2 2 3 1 2 3 4 1 and (vi) X , where a, b and c are arbitrary real constants. 12 −1/2 Case (i): The symmetry generator X +aX gives the group-invariant solution of the form u(x, t) = t [A(y) cos(a/2 ln t)+ 5 3 −1/2 −1/2 B(y) sin(a/2 ln t)] and v(x, t) = t [A(y) sin(a/2 ln t) − B(y) cos(a/2 ln t)] with y = xt . This reduces the system of PDEs (33) to the following system of nonlinear second-order ODEs ′′ ′ 2 2 2A − aA − yB − B + 2γA(A + B ) = 0, (36a) ′′ ′ 2 2 2B − aB + yA + A + 2γB(A + B ) = 0. (36b) The above system of ODEs is found to be Painleve´ integrable and sufﬁcient number of independent integrals can be found as in the case of Eq. (15) to show the integrability. The direct numerical analysis of the above coupled ODEs is shown below in Fig. 13. Here the explicit appearance of y in the ﬁrst derivative term leads to a compression of wave train as y increases Fig. 13(a)-(b). Meanwhile, the amplitude reaches a steady state. This resembles the behaviour of an attractor. In the absence of linear coupling terms, one can obtain a periodic waves with diminishing amplitude as y increases. We have shown such nonlinear wave trains in Fig. 13(c)-(d). Figure 13: Numerical nonlinear wave trains in Eq. (36) γ = 2. The upper and lower panels show the wave trains in the presence (a = 1) and absence (a = 0) of linear coupling terms. Case (ii): For the sub-algebra with generators X + cX , one can obtain the group-invariant solution u = A(y) and 1 2 v = B(y) with the invariant y = x − ct. From this, the system of PDEs (33) can be reduced to the system of nonlinear second-order ODEs, ′′ ′ 2 2 A + cB + γA(A + B ) = 0, (37a) ′′ ′ 2 2 B − cA + γB(A + B ) = 0. (37b) Note that the above equation (37) also passes the Painleve´ test and becomes integrable. We wish to mention that the above equation (37) is similar to Eq. (19) and we can obtain multi-peak nonlinear wave structures as shown in Fig. 4. If c = 0, the group-invariant solution corresponding to X of the system of PDEs (33) is u = A(y), v = B(y), where y = x, and A(y) and B(y) satisfy the following system of nonlinear second-order ODEs, ′′ 2 2 A + γA(A + B ) = 0, (38a) ′′ 2 2 B + γB(A + B ) = 0. (38b) Here we ﬁnd that the above system of ODEs (38) is the same as the integrable equation (23) and it admits periodic as well as elliptic function solutions which have already been addressed. 13 Case (iii): In this case, for X + cX , we obtain the group-invariant solutions u = −A(y) sin cx + B(y) cos cx and 2 3 v = A(y) cos cx + B(y) sin cx with the invariant y = t. From this and Eq. (33), we get the following set of coupled nonlinear ﬁrst-order ODEs, 2 2 2 A + c B − γB(A + B ) = 0, (39a) 2 2 2 B − c A + γA(A + B ) = 0. (39b) Equation (39) is equivalent to Eq. (27) and exhibits the standard periodic wave solutions as shown in Fig. 8. One can solve the above system of ODEs (39) and obtain the following form of group-invariant solutions for (33) u(x, t) = b cos((γb − c )t + cx + b ), (40a) 1 1 2 v(x, t) = b sin((γb − c )t + cx + b ), (40b) 1 1 2 where b and b are constants. 1 2 If c = 0, the group-invariant solution corresponding to X of the system of PDEs (33) is u = A(y), v = B(y), where y = t, and A(y), B(y) satisfy the following system of nonlinear ﬁrst-order ODEs, 2 2 A + γB(A + B ) = 0, (41a) 2 2 B − γA(A + B ) = 0. (41b) The above system of ODEs (41) is similar to Eq. (17). The solution of (41) gives rise to the following group-invariant solution for the system of PDEs (33), which is independent of x, u(x, t) = b cos(γb t + b ), (42a) 1 1 2 v(x, t) = b sin(γb t + b ), (42b) 1 1 2 where b and b are constants. 1 2 Case (iv): For X + bX + cX , we obtain the invariant y = x − ct and the corresponding invariant solution u = 1 3 2 −A(y) sin bt + B(y) cos bt and v = A(y) cos bt + B(y) sin bt. This invariant solution reduces the system of PDEs (33) to the system of coupled nonlinear second-order ODEs, ′′ ′ 2 2 A − bA − cB + γA(A + B ) = 0, (43a) ′′ ′ 2 2 B − bB + cA + γB(A + B ) = 0. (43b) Equation (43) is equivalent to Eq. (25) for k = 0 and is also Painleve´ integrable. Here one can obtain different kinds of nonlinear wave trains for various choices of the arbitrary constants b and c. If c = 0, then the symmetry generator X + bX reduces the system of PDEs (33) to the system of nonlinear 1 3 second-order ODEs ′′ 2 2 A − bA + γA(A + B ) = 0, (44a) ′′ 2 2 B − bB + γB(A + B ) = 0. (44b) Note that the above system of ODEs (44) can be viewed as a special case of Eq. (25) for k = a = 0, which is also a Painleve´ integrable. Interestingly, Eq. (44) admits the following solitary wave solutions: A = a sech(k y + k ), (45a) 1 1 2 B = a sech(k y + k ), (45b) 2 1 2 2 2 2b 2 with the constraints k = b and a = − a . This is shown in Fig. 14 for different choices of arbitrary constants. 1 2 γ 1 3 2 3 2 Case (v): We obtain the group-invariant solution u = A(y) sin(tx/c − 2t /3c ) and v = −A(y) cos(tx/c − 2t /3c ) from generator X + cX along with an invariant as y = cx − t . This symmetry invariant reduces Eq. (33) to a nonlinear 4 1 second-order ODE of the form ′′ 1 γ A − yA + A = 0. (46) 4 2 c c 14 Figure 14: (Top panel) Solitary wave solutions of Eq. (44) for (a) a = 1.5 and (b) a = 1.0 with other parameters as γ = 2, b = 1, and k = 5. 1 1 2 (Bottom panel) The stationary solitary waves with periodically varying amplitude of NLS equation (33) for the choice of (a) with y = x. Figure 15: Exponentially growing periodic wave trains of Eq. (46) for c = 1 and γ = 2. We ﬁnd that the above equation is also Painleve´ integrable. The direct numerical analysis shows that the explicit ap- pearance of ‘y’ in the model equation induces the modulation of the periodic nonlinear wave structure as it progresses and it can be evidenced from Fig. 15. Case (vi): The symmetry generator X does not give any group-invariant solutions. 5. Conclusions We have investigated the symmetry reductions of the nonlinear Helmholtz equation arising in the context of non- linear optics by using the Lie symmetry analysis. Especially, with an inﬁnitesimal transformation, we have obtained the symmetry generators (vector ﬁelds), identiﬁed a set of invariants from the optimal system of one-dimensional subalgebras and symmetry reductions in the form of coupled ODEs. We have studied the integrability property of the reduced ODEs by using the Painleve´ analysis and constructed explicit ﬁrst integrals by using the modiﬁed Prelle- Singer method. In the study of NLH system, we have also obtained nonlinear periodic and solitary wave solutions for some ODEs resulting for speciﬁc symmetry reductions. For the remaining ODEs, we have given a qualitative picture of the solutions by numerically solving the underlying ODEs. We have explored some interesting nonlinear wave structures resulting for different choices of arbitrary parameters. The obtained nonlinear periodic wave solutions can 15 ﬁnd applications in the context of ﬁber Bragg gratings, bimodal ﬁbers, etc. Further, we have compared the known Lie symmetry analysis of the nonlinear Schro¨dinger equation with the results of the present nonlinear Helmholtz equa- tion. Our analysis shows the possible symmetry reductions of a scalar nonlinear partial differential equation. Thus our study provides an analytical treatment of the interesting NLH system in view of Lie point symmetry, integrability and invariant solutions. As a future study, the present work will be extended to investigate the Lie symmetry and group invariant solutions of various coupled nonlinear evolution equations arising in different contexts. Acknowledgments KS is grateful to DST-SERB for the award of a National Post-Doctoral Fellowship (PDF/2016/000547). AGJ would like to express his sincere thanks and gratitude to the Centre for Nonlinear Dynamics, School of Physics, Bharathidasan University, Tiruchirappalli, India, for its warmest hospitality and support during his visit. TK is sup- ported by the DST-SERB research project (EMR/2015/001408). The research work of ML is supported by a NASI Senior Scientist Platinum Jubilee Fellowship (NAS 69/5/2016-17) and DST-SERB Distinguished Fellowship (DO No. SB/DF-09/2017). ML is also supported by the research projects of DAE-NBHM (2/48(5)/2015/NBHM(R.P.)/R&D II/14127), CSIR India (03/1331/15/EMR-II) and DST-SERB (EMR/2014/001076). Appendix A. Painleve´ analysis of the ODE (15) It is well known that the Painleve´ analysis gives evidence for the integrability nature of a given system of ODEs/PDEs through an algorithmic procedure, namely leading order analysis, identiﬁcation of resonances and ﬁnd- ing the existence of sufﬁcient number of arbitrary functions at each resonances called arbitrary analysis. For a detailed algorithm see  and references therein. In this Appendix, we carry out the Painleve´ analysis for the ODE (15) to study its integrability nature. We express the dependent variables of real equation (15) as generalized Laurent series ∞ ∞ X X j+α j+β expansion (A = a φ and B = b φ , where α and β are negative integers yet to be determined) in the j j j=0 j=0 neighbourhood of the non-characteristic singular manifold φ(y). By terminating the Laurent series upto the zeroth α β −3 order (A ≈ a φ and B ≈ b φ ), we get the leading order of the ODE (15) as φ which results for α = β = −1. Also, 0 0 2 2 2 −3 we obtain the leading order equation as 2(1 + ka ) + γ(a + b ) = 0 as the coefﬁcient of φ . 0 0 From the generalized Laurent series and by using the above leading order equation, we obtain the characteristic 2 2 equation as γ(1 + ka ) ( j + 1) j( j − 3)( j − 4) = 0. So, the resonances for Eq. (15) are obtained as j = −1, 0, 3, 4, where the resonance j = −1 corresponds to the arbitrariness of the non-characteristic manifold. The arbitrary analysis shows that there exists required number of arbitrary functions at each of these resonances. So, we conclude that the coupled second-order ODE (15) is Painleve´ integrable. Similar to the above analysis for Eq. (15), we have carried out the Painleve´ analysis for other system of ODEs given in this work to identify their integrability nature. The ﬁnal result of the analysis is given below the respective equations. Appendix B. Modiﬁed Prelle-Singer method: Solutions for coupled second order ODEs In this Appendix, we brieﬂy present the historical developments of the Prelle-Singer method and the important steps involved in the modiﬁed Prelle-Singer method. Prelle and Singer, in 1983, proposed a method for ﬁnding the general solutions of ﬁrst-order ODEs . This method guarantees the solution of a given ﬁrst-order ODE if the equation under consideration admits an elementary solution. Later, Duarte et al modiﬁed the technique developed by Prelle and Singer and applied it to second order ODEs . Their approach was based on the conjecture that if an elementary solution exists for the given second order ODE then there exists at least one elementary ﬁrst integral I(t, x, x˙) whose derivatives are all rational functions of t, x and x˙. This method was recently extended by Chandrasekar, Senthilvelan and Lakshmanan to second order ODEs . They have generalized the theory of Duarte et al and shown that for the second order ODEs one can isolate even two independent integrals of motion and obtain general solutions explicitly without any integration. This method is known as the modiﬁed Prelle-Singer method. The method has been extended by these authors to coupled ODEs and higher order ODEs as well. 16 Let us brieﬂy discuss the modiﬁed PS procedure applicable to second order ODEs  which we have adopted in this manuscript. We assume that the following general second order ODE, x¨ = φ(t, x, x˙), (B.1) where the over dot denotes differentiation with respect to time, admits a ﬁrst integral I(t, x, x˙) = constant. We know that the total derivative of I vanishes, that is, dI = I dt + I dx + I dx˙ = 0, (B.2) t x x˙ where the subscript denotes partial differentiation with respect to that variable. Let us now rewrite Eq. (B.2) in the form φdt − dx˙ = 0 and add a null term S (t, x, x˙)x˙dt - S (t, x, x˙)dx. Then we get (φ + S x˙)dt − S dx − dx˙ = 0. (B.3) Multiplying the above equation by a factor R(t, x, x˙), which acts as the integrating factor, Eq. (B.3) becomes an exact equation of the form R(φ + S x˙)dt − RS dx − Rdx˙ = 0. Now, comparing this equation with Eq. (B.2) we get I = R(φ + x˙S ), I = −RS, I = −R. (B.4) t x x˙ Therefore, the ﬁrst integral I of Eq. (B.1) can be obtained by integrating the above system of equations (B.4) as I = r − r − R + (r − r ) dx˙, (B.5) 1 2 1 2 dx˙ R R where r = R(φ + x˙S )dt, r = RS + r dx. Note that the ﬁrst integral I(t, x, x˙) is given in terms of R and 1 2 1 dx S which have to be determined. In order to ﬁnd the explicit forms of R and S for a given form of φ, we use the compatibility conditions (I = I , I = I , and I = I ) and the following determining equations for R and S : tx xt tx˙ x˙t xx˙ x˙x S + x˙S + φS = −φ + φ S + S , (B.6a) t x x˙ x x˙ R + x˙R + φR = −(φ + S )R, (B.6b) t s x˙ x˙ R − S R − RS = 0. (B.6c) x x˙ x˙ We would like to emphasis here that any two independent particular solutions (S, R) of the above system of PDEs (B.6) are enough to deduce two independent ﬁrst integrals of Eq. (B.1). In order to obtain a particular solution of (B.6), one can use suitable ansatz for S and R and ﬁnd a compatible solution. Chandrasekar, Senthilvelan and Lakshmanan have also developed a straightforward procedure to deduce the time independent conservative Hamiltonian structure from the time independent ﬁrst integral obtained by the modiﬁed PS procedure . Using this procedure the Liouville sense of integrability can be established for the systems for which one is unable to explicitly obtain the second integral of motion. References  Bluman GW, Kumei S. Symmetries and differential equations. Springer, New York, 1989.  Ibragimov NH. CRC Handbook of Lie group analysis of differential equations, Vol. 1-3. CRC Press, Florida, 1993-1996.  Olver PJ. Applications of Lie Groups to Differential Equations. Second Edition, Springer, New York, 1993.  Ovsiannikov LV. Group Analysis of Differential Equations. Academic Press, New York, 1982.  Kaliappan P. Ph.D. Thesis: Studies on the invariances and symmetries of certain nonlinear evolution equations. University of Madras 1981.  Lakshmanan M, Kaliappan P. J Math Phys 1983;24:795.  Olver PJ, Rosenau P. SIAM J App Math 1987;47:263. Group-Invariant Solutions of Differential Equations.  Gagnon L, Winternitz P. J Phys A: Math Gen 1988;21:1493, ibid 1989;22:469, ibid 1989;22:499, Fushchych W, Serov M. J Phys A: Math Gen 1987;20:L929, Fushchich WI, Moskaliuk SS. Lettere al Nuovo Cimento 1981;31:571.  Kumei S. J Math Phys 1977;18:256, Tajiri M. J Phys Soc Japan 1983;52:1908, Clarkson PA. Nonlinearity 1992;5:453, Hereman W. Math Comput Modelling 1997;25:115. 17  Crosignani B, Porto PD, Yariv A. Opt Lett 1997;22:778.  Chamorro-Posada P, McDonald GS, New GHC. J Modern Opt 1998;45:1111.  Blair S. Chaos. 2000;10:570, Crosignani B, Yariv A, Mookherjea S. Opt Lett 2004;29:1254.  Chamorro-Posada P, McDonald GS, New GHC. J Opt Soc Am B. 2002;19:1217, Ciattoni A, Crosignani B, Mookherjea S, Yariv A. Opt Lett 2005;30:516.  Tamilselvan K, Kanna T, Avinash Khare. Commun Nonlinear Sci Num Simul 2016;39:134.  Chiao RY, Garmire E, Townes CH. Phys Rev Lett 1964;13:479, Talanov VI. ZhETF Pis Red 1965;2;223 [JETP Lett 1965;2:138], Kelley PL. Phys Rev Lett 1965;15:1005, Wang CS. Phys Rev 1968;173:908, Zakharov VE. Zh Eksp Teor Fiz 1967;53:1735 [Sov Phys JETP 1968;26:994].  Zakharov VE, Shabat AB. Zh Eksp Teor Fiz 1971;61:118 [Sov Phys JETP 1972;34:62].  Cherniha R. Rep Math Phys 2000;46:63, Cherniha R, King JR. J Phys A: Math Gen 2000;33:267, ibid 2003;36:405, Nikitin AG, Wiltshire RJ. J Math Phys 2001;42:1666, Nikitin AG. J Math Anal Appl 2006;324:615.  Ablowitz MJ, Ramani A, Segur H. J Math Phys 1980;21:715, Ramani A, Grammaticos B, Bountis T. Phys Rep 1989;180:160, Lakshmanan M, Sahadevan R. Phys Rep 1993;224:1.  Prelle M, Singer M. Trans Amer Math Soc 1983;279:215.  Duarte LGS, Duarte SES, daMota LACP, Skea JEF. J Phys A: Math Gen 2001;34:3015.  Chandrasekar VK, Senthilvelan M, Lakshmanan M. Proc R Soc A 2005;461:2451, Chandrasekar VK, Senthilvelan M, Lakshmanan M. Proc R Soc A 2009;465:609.  Cattani F, Kim A, Hansson T, Anderson D, Lisak M. Euro Phys Lett 2011;94:53003, Liu C, Yang ZY, Zhao LC, Duan L, Yang G, and Yang WL. Phys Rev E 2016;94:042221.
Nonlinear Sciences – arXiv (Cornell University)
Published: Mar 5, 2018
Access the full text.
Sign up today, get DeepDyve free for 14 days.