In this study, we investigate the typical systems modeled by the (3 + 1)-dimensional as well as (1 + 1)-dimensional Schrödinger equations incorporating third-order dispersion effects, higher-order scattering effects, and cubic–fifth–seventh degree nonlinear interactions. We use the F-expansion method and the self-similar method to solve the higher-order Schrödinger equation for one-dimensional and three-dimensional settings, respectively, identifying typical bright soliton solutions under appropriate system settings. The bright soliton features are demonstrated analytically in regions around the soliton peak region. Pictorial bright soliton features are demonstrated for the three-dimensional setting as well as one-dimensional setting. Our work shows the applicability of the theoretical treatment utilized in studying bright soliton dynamics for systems with third-order dispersion and seventh degree nonlinearity.
Due to the balance of the dispersion effect and scattering effect,1–3 solitons exhibit robust stability in the process of propagation and interaction.4 It is an extremely important phenomenon in modern physics, which has attracted extensive attention in theoretical and experimental research due to its unique properties. The nonlinear Schrödinger equation (NLSE) is an appropriate choice for theoretical study of soliton behavior and plays an important role in nonlinear optics,5,6 condensed matter physics,7,8 and plasma physics.9,10 Specifically, when the dispersion effects of the NLSE are considered up to the third order, besides the cubic–fifth degree nonlinearity (in the atomic system) caused by the two-body and three-body effects,11 the nonlinearity up to the seventh degree12–15 by incorporating the multi-body effect should be considered. In addition, the multi-body effect of the quantum system should also be considered, so the higher-order scattering effect should be incorporated.
In this work, we study the (3 + 1)-dimensional as well as (1 + 1)-dimensional higher-order Schrödinger equations incorporating the third-order dispersion effect and the cubic–fifth–seventh degree nonlinear interaction,16 and we use the F-expansion17,18 method and the self-similar method19–21 to solve the three-dimensional as well as one-dimensional higher-order NLSE under appropriate parametric settings. We first identify the bright soliton solutions by deriving the one-dimensional solution of the bright soliton type, then the three-dimensional bright soliton solution is identified, and we then derive the key characteristics of bright solitons around the location of the soliton peak, demonstrating typical bright soliton features. The obtained analytical results can be used to guide the experimental detection of bright solitons in (3 + 1)-dimensional as well as (1 + 1)-dimensional systems with high-order dispersion and seventh degree nonlinear effects.
This work is organized as follows: in Sec. II, the third-order NLSE incorporating cubic–fifth–seventh degree nonlinearity and high-order dispersion terms is analyzed, followed by the introduction of the F-expansion method. In Sec. III, we use the F-expansion method to solve such a NLSE in the one-dimensional setting, with typical bright soliton analytical features demonstrated and pictorially displayed. Section IV presents the derivation of the bright soliton solution of this category of NLSE in the three-dimensional setting, with typical bright soliton features pictorially demonstrated. Section V presents typical application scenario demonstration of our work and stability analysis of our derived bright soliton solutions. Section VI gives the conclusion remarks.
II. NONLINEAR SCHRÖDINGER EQUATION MODEL AND F-EXPANSION METHOD
A. Nonlinear Schrödinger equation involving cubic–fifth–seventh degree nonlinearity and higher-order dispersion terms
Incorporating higher-order dispersion and cubic–fifth–seventh degree nonlinear interaction effects, the (1 + 1)-dimensional Schrödinger equation is expressed as follows:
The function ψ(x, t) denotes the complex valued function of x and t, where x and t denote the space and time coordinates, respectively, in cold atomic systems for example. Additionally, t denotes the spatial coordinate, and x denotes time in the optical system. a1 and a2 denote second and third-order dispersion coefficients, respectively. Furthermore, b1, b2, and b3 denote the parameters of the cubic, fifth degree, and seventh degree nonlinear interactions that are attributed to two-body and three-body effects, respectively, in the atomic system. The parameters λ1, λ2, and λ3 denote the coefficients of the leading and higher-order scattering terms. To proceed, we give a general introduction of the F-expansion method to be used.
B. F-expansion method
In order to solve Eq. (1), we use the F-expansion method, which solves the nonlinear partial differential equation taking the following form:
Equation (2) is expressed as the polynomial of an unknown function u(x, t) and its partial derivatives of various orders. The unknown function u(x, t) is to be expressed by the polynomial of base-function F(ξ), with F(ξ) defined as a function of ξ = px + qt as
where n, n − 1, and …, a3, a2, a1, and a0 are constants. Furthermore, if we differentiate both sides of Eq. (3) with respect to ξ, we obtain the transformed equation for F(ξ) as
u(x, t) is expressed as
where usually hi(t) is just a parametric constant (also constant in this work) hi, which is to be determined by the ensuing equation solving steps. By substituting Eq. (5) into Eq. (2) and using Eq. (3) or Eq. (4), we reach the polynomial in terms of F(ξ) and F(ξ) times dF(ξ)/dξ. m is set by balancing between the highest-order nonlinear term and highest-order differential term. For example, if in Eq. (2), the highest order for the nonlinear term is k (k ≥ 2, nonlinear term uk, for example), the highest-order for the differential term is r (, for example), we need at least two highest order terms Fs of order s in Eq. (2) by substituting Eq. (5) into Eq. (2) and utilizing formula (4) or (5), so that hm can be nonzero. The highest order (of ) from the nonlinear term is s1 = km, the highest order (of ) from the differential term is s2 = r(n/2 − 1) + m, by letting s1 = s2 = km = r(n/2 − 1) + m, we can solve for m. For example, for a traditional cubic nonlinear Schrödinger equation and n = 4 in definition (3), the solution is m = 1 [so that ansatz for u(x, t) is u(x, t) = h1F(ξ) + h0]. The resultant analytical form for Eq. (2) is a polynomial of F on the left-hand side (LHS) and zero on the right-hand side, and the coefficients of the LHS polynomial are expressions of hi(t) and ai.
In order to solve Eq. (2), the coefficients of all terms are set to zero. This generates a series of ordinary differential equations [ODEs, for functions hi(t)] or algebraic equations (for parameters hi), whose solutions determine the formulation of Eq. (5), so that Eq. (1) is solved accordingly. According to the definitions (4) and (5) of the base function F(ξ), it is of elliptical function feature in most scenarios with its precise analytical form determined by the coefficient values solved by the ODEs mentioned above.
III. SOLITON SOLUTION FOR (1 + 1)-DIMENSIONAL NLSE WITH SEVENTH DEGREE NONLINEARITY AND HIGHER-ORDER DISPERSION
We consider the traveling wave solution format of Eq. (1) as follows:
where ξ = px + qt, A, B, p, and q are constants to be determined. φ is the modulus of wave function Ψ, that is φ = |Ψ|. Substituting Eq. (6) into Eq. (1), the expressions of each differential term take the following form:
The imaginary part of Eq. (1) takes the following form:
Utilizing and with the following definition:
After integration (we require that the integration constant is zero, since ξ → 0, φ → 0, and ), Eq. (11) is transformed to the following form:
where , and it can be seen that w is symmetrical about axis ξ = 0 in Eq. (14), when ξ varies from 0 to +∞, w(ξ) monotonically increases from the minimum value c1 + c2 to infinity, according to Eq. (13), φ monotonically decreases from the maximum value at ξ = 0 to the zero value when ξ → ∞, which are typical features of bright solitons. Such a bright soliton feature is visually shown (in Fig. 1) via numerical evaluation of Eq. (14). According to the traveling wave ansatz [analytical solution (6)] of φ, the solution of φ is just a bright soliton solution (with a positive peak of a constant shape and tends to zero as ξ → 0, which is the feature of bright soliton).
When ξ → 0, w → c1 + c2, near ξ = 0, ignoring higher-order terms of ξ, the approximate analytical form of Eq. (14) is
It is just the classical form of solution of the bright soliton category, which is a traveling wave with stable shape that comes into being by the balance between nonlinear and dispersion effects with the center peak possessing the highest positive amplitude. Figure 1 shows the bright soliton solution plot for |ψ| by numerically evaluating Eq. (14) with a typical setting of ci(i = 1, 2, 3). With such a typical bright soliton solution for one-dimensional Eq. (1), we can build its three-dimensional solution based on the self-similar method, which will be studied in Sec. IV.
We use the method22 to verify the correctness of our analytical solution. Let ,
For bright solitons near ξ = 0, φ′ = 0, the above formula expands Taylor near ξ = 0, φ = cv, and we can get
It is not difficult to get
IV. BRIGHT SOLITON BEHAVIOR OF THREE-DIMENSIONAL SEVENTH DEGREE NONLINEAR SYSTEM WITH HIGHER-ORDER DISPERSION
Incorporating the third-order dispersion and seventh degree nonlinear effects, the three-dimensional nonlinear Schrödinger equation, which is the three-dimensional analog of Eq. (1), takes the following form:
where is the velocity direction vector of system flow, (x1, x2, x3) = (x, y, z). To solve Eq. (20), We consider its self-similar ansatz (analytical solution, …) of the form as follows:
The self-similar formulation requires that Eq. (22) has the same form as Eq. (1), so the coefficient of the first derivative of ψ1D with respect to ν in Eq. (22) must be zero, this requires that A′(x, y, z) and ν(x, y, z) are linear functions of x, y, z, and τ(t) is a linear function of t as follows:
The degree of freedom for is 2, and we can choose the direction of as (1, 0, 0); Eq. (27c) becomes
where q, d, and B are constants determined by the initial conditions of the system. ψ1D is the form of bright soliton solution. When pν + qt → 0, , so . The precise waveforms are shown in Fig. 2 by precise numerical evaluation of self-similar ansatz ψ1D. From Eq. (29), we can see that the flow velocity is along the direction of , the peak moving velocity is along the direction , and and are generally non-collinear, so solution (29) is truly a three-dimensional solution.
V. PRACTICAL APPLICATION OF SOLITON DYNAMICS IN SEVENTH DEGREE NONLINEAR SYSTEM AND STABILITY ANALYSIS
The seventh degree nonlinearity related analysis has important applications in the non-Kerr medium, where the propagation of optical pulse with ultrashort femtosecond features has the following form:16
This (1 + 1)-dimensional equation [Eq. (30)] (with seventh degree nonlinearity but with only second-order dispersion) is very similar to the (1 + 1)-dimensional equation model presented in this work. In addition, the bright soliton profile of the complex value electric field E in Eq. (30) as shown in Fig. 3 in the non-Kerr medium work16 is very similar to that shown in Fig. 1 of this work.
With regard to the applicability of our theoretical analysis to experimental observation, it is important to perform some stability analysis of the base bright soliton solution |ψ1D| = φ≃ in Eq. (16), with φ’s amplitude , soliton peak distribution width . According to the analytical expressions for c1, c2, and c3 following Eq. (14), δ only depends on the cubic nonlinearity constant λ1 via c3 through β1, so it is stable relative to the variation of fifth degree and seventh degree nonlinear strength constants λ2 and λ3. This means that the soliton is stable (against wave shape spreading) relative to variations of quintic and seventh degree nonlinearity strength constants and combining Eq. (12), , and . Since β2 ≪ β1 (λ2 ≪ λ1) and β3 ≪ β1 (λ3 ≪ λ1), the soliton peak amplitude is more sensitive to the variation of fifth degree and seventh degree nonlinear strength constants. However, this will not affect the stability (determined by the sensitivity of δ, which is independent of λ2 and λ3) of the bright soliton solution ψ1D derived in this work.
In this study, we investigated the (3 + 1)-dimensional as well as (1 + 1)-dimensional higher-order cubic–fifth–seventh degree nonlinear Schrödinger equation (NLSE) with higher-order dispersion effects. We focused on identifying the bright soliton behavior of the system modeled by the higher-order NLSE. We use the F-expansion method and self-similar approach to determine the bright soliton solution. Through numerical evaluation of the inexplicit analytical form of the bright soliton solution, the typical bright soliton characteristics of the three-dimensional solution as well as the one-dimensional solution are identified, and then, the analytical expression of bright soliton solution around the peak region is derived analytically, with the exact characteristics of bright soliton solutions and stability feature illustrated. The theoretical results obtained from the NLSE model with seventh degree nonlinear interaction and the high-order dispersion effect can be used as an experimental guide for detecting the behavior of bright solitons in optical or ultracold atomic systems with seventh degree nonlinearity and higher-order dispersion.
Y.G. and Q.C. contributed equally to this work.
This work was supported by the National Natural Science Foundation (NSF) of China under Grant No. 11547024 and the Postgraduate Research and Practice Innovation Program of Jiangsu Province under Grant No. KYCX20_3113.
Data sharing is not applicable to this article as no new data were created or analyzed in this study. Our study is theoretical and analytical work of the corresponding theoretical models. We do not use previously published data in our work.