The step-wise construction of solitary solutions to Riccati equations with di ff usive coupling

: A novel scheme based on the generalized di ff erential operator and computer algebra was used to construct solitary solutions to a system of Riccati di ff erential equations with di ff usive coupling. The presented approach yields necessary and su ffi cient existence conditions of solitary solutions with respect to the system parameters. The proposed stepwise approach enabled the derivation of the explicit analytic solution, which could not be derived using direct balancing techniques due to the complexity of algebraic relationships. Computational experiments were used to demonstrate the e ffi cacy of proposed scheme.


Introduction
Solitary solutions (or solitons) have been a mainstay in the field of nonlinear dynamics since their rediscovery in the Fermi-Pasta-Ulam-Tsingou experiment after almost a century of obscurity [1].While the roots of solitary solutions are in physics, the analysis of such phenomena has branched out into numerous other disciplines including (but not limited to) engineering [2], biology [3,4] and neuroscience [5].
Some examples of recent publications on solitary solutions in various research fields are given below.Symbolic computation is used in conjunction with the modified rational expansion method to construct several types of solitary solutions to the Jimo-Miwa equation [6].A soliton potential of Bose-Einstein condensates loaded into a one-dimensional four color lattice is constructed using analytical techniques in [7].Four types of solitary to the Camassa-Holm equation, which are used to model shallow water dynamics, are constructed in [8].The existence of polariton-type solitons in microcavity wires is investigated in [9].
Bifurcation and traveling wave solutions to the Manakov model (an extension of Schrödinger's equation) are derived in [10].A fractional time-coupled Schrödinger's model and its exact traveling wave solutions are considered in [11].Stationary optical solitons and exact solutions to a generalized Shcrödinger system important in the optical fiber theory are derived in [12], while such solutions to another optical system are considered in [13].Exact solutions to a stochastic fractional long-short wave interaction system are first obtained in [14].
A deep-learning based method for the construction of solitary solutions to higher-order Kortewegde-Vries (KdV) and Boussinesq equations is presented in [15].The wave breaking of solitary wave solutions to the Novikov equation, a higher dimension generalization of the KdV model, is considered in [16].Soliton solutions for the nonlinear sine-Gordon equation with Neumann boundary conditions are constructed using a meshless collocation method in [17].Bright and dark solitary solutions to a heat transport with relaxational effects at a nanoscale are shown to exist in [18].
A class of ansatz methods have been successfully applied to various types of Schrödinger's equations [19,20] and other models of quantum physics [21].More recently, ansatz techniques have been used to construct Gaussian solutions for a class of Schrödinger equations with logarithmic nonlinearity [22] and analytical traveling-wave solutions to a generalized Gross-Pitaevskii equation [23].
The main objective of this paper is to construct kink solitary solutions to the following system of Riccati equations with diffusive coupling via the generalized differential operator technique: where The construction of solitary solutions to Riccati-type systems similar to (1.1) have been the subject of intensive research during the last decades.A variety of approaches based on the generalized differential operator have been developed [24].Such techniques are superior to direct balancing (or ansatz) methods, as they can be used not only to construct solitary solutions, but also the necessary and sufficient existence conditions for the existence of such solutions.Bright and dark solitary solutions and their necessary and sufficient existence conditions in a Riccati-type with multiplicative coupling are constructed in [25].
Note that (1.1) has only diffusive coupling terms.Due to this absence of interaction between diffusive and multiplicative coupling, (1.1) possesses unique dynamics and properties with respect to its solitary solutions.Multiplicative coupling of the same system (terms γxy in both equations) has already been considered in [25].An important property of multiplicative coupling is that the coupling term can only become equal to zero if one of the solutions becomes zero.This is not the case for diffusive coupling -the coupling term becomes zero if the functions x and y coincide.Due to this effect, it can be expected that diffusive coupling would yield a larger nomenclature of solitary solutions compared to the multiplicative case.
Furthermore, the generalized differential operator scheme presented in this paper discusses a novel way of applying such operators in conjunction with computer algebra for the stepwise construction of kink solitary solutions.Moreover, the techniques outlined in detail in this paper can be applied to any system of differential equations and further extended to construct higher-order (such as bright/dark) solitary solutions.

Solitary solutions
Solitary solutions considered in this paper have the following form [26,24]: (2.1) where . ., n and n defines the order of the solitary solution.Note that n = 1 and n = 2 correspond to kink and dark/bright solitary solutions [26].Typical waveforms of kink and dark/bright solitary solutions are depicted in Figure 1.Note that kink solutions are monotonous, while higher-order solitary solutions can have local extrema: In general, an n-th order solitary solution (2.1) can have at most n − 1 minima or maxima [27].

The inverse balancing technique
The inverse balancing technique is applied before the construction of solitary solutions in order to determine the necessary existence conditions of such solutions in (1.1).Note that inverse balancing is distinct from direct balancing (also referred to as ansatz methods) in the sense that the goal of the procedure is not to construct solutions to differential equations, but to determine the necessary existence conditions for solutions of some fixed analytical form.
The first step of this procedure is to insert a solitary solution as an ansatz into the given system of differential equations, thus, obtaining a system of linear equations with respect to the parameters of the differential equations.If the resulting system is solvable with some conditions imposed on the parameters of the solitary solutions, then those conditions correspond to the existence criteria for the respective solitary solution.However, if the obtained linear system is degenerate, then the solitary solution of respective order cannot exist in the analyzed model.A broader discussion and examples of this technique are given in [25].

The generalized differential operator technique
Let us consider the following system of ordinary differential equations: where P and Q are analytic functions.
Let D α denote the partial differentiation operator with respect to the variable α.The generalized differential operator D cuv corresponding to system (2.8) is defined as follows [24]: (2.9) The solution to system (2.8) can be expressed using (2.9) in the following form [24]: Let: and where H (m) p denotes the m-th order Hankel determinant with respect to the sequence p j ; j = 0, 1, . . . .The sequence p j ; j = 0, 1, . . . is called an m-th order linear recurring sequence if there exists such m ∈ N that satisfies the following condition: The elements of the sequence p j ; j = 0, 1, . . .can then be expressed as follows [28]: where λ 1 , . . ., λ m are constant coefficients and ρ 1 , . . ., ρ m are roots of the following characteristic polynomial: then q j ; j = 0, 1, . . . is a linear recurring sequence of order m and can be expressed as: where µ, . . ., µ m are constant coefficients and ν 1 , . . ., ν m are roots of the characteristic polynomial analogous to (2.15).

Application of linear recurring sequences in (2.10)
It has been shown in [24] that if the following conditions hold true for all c, u, v: then series solutions (2.10) can be expressed in the closed-form: Note that if both sets of characteristic roots ρ k , ν k (k = 1, . . ., m) each have at least one zero, then expressions (2.20) correspond to the solitary solutions (2.6)-(2.7).In relation to this, the following property of linear recurrence sequences is presented.Let w j ; j = 0, 1, . . .be an m-th order linear recurring sequence with one zero characteristic root.Then, the truncated sequence w j ; j = 1, 2, . . . is (m − 1)-th order linear recurring sequence [25].In other words, when constructing solitary solutions (2.6)-(2.7),one may consider the sequences starting with the elements p 1 and q 1 , respectively.
The procedure used for the construction of solitary solutions is summarized in Figure 2. 3. The determination of the maximum possible order of solitary solutions to (1.1) In this section, the maximum order of the solitary solutions (2.4)-(2.5) to (1.1) is determined by applying the inverse balancing technique outlined in the Section 2.2.First, the system of differential equations (1.1) is transformed via the substitutions (2.3) as follows: Note that substitution (2.3) is necessary -since the sequence of coefficients (2.15) does not form a linear recurring sequence for the untransformed system (1.1).Let the order of the solitary solutions be n = 2 (dark/bright solitons): 3) into (3.1) and taking t = t 1 , t 2 , x 1 , x 2 , y 1 , y 2 , 0 yields the following system of linear equations with respect to the parameters The solutions to the obtained linear system (3.4) exist only if one or both of the following conditions hold true: The first case -corresponding to n = 1 -will be considered in the further sections, while the second condition is trivial since it leads to a linear system of differential equations (1.1).Therefore, the maximal order of solitary solutions to (1.1) is n = 1, and bright/dark or higher order solitary solutions do not exist in the analyzed system.Inverse balancing technique (see Section 2.2) is applied in this section in order to derive the necessary existence conditions of kink solitary solutions to (1.1).
Let the order of the solitary solutions be n = 1: Inserting (4.1)-(4.2) as an ansatz into (3.1) and taking t = t 1 , x 1 , y 1 , 0 yields the following system of linear equations with respect to the parameters γ 1 , γ 2 , a k , b k (k = 0, 1, 2): 3) The solution to the obtained linear system (4.3)reads: where γ 1 , γ 2 ∈ R are chosen arbitrarily.Note that the linear system (4.3) can be solved without any necessary existence conditions imposed on the parameters of the solitary solutions (4.1)-(4.2).

The determination of the parameter η of kink solitary solutions (4.1)-(4.2) to (1.1)
In this section, parameter η of kink solitary solutions (4.1)-(4.2) is determined by applying the results presented in Sections 2.3 and 2.4.
The procedure outlined above is depicted in the schematic diagram in Figure 3.Note that the algorithm displayed in Figure 3 has to be initiated from the coefficient that corresponds to the highest power of u, thus solving the equation H (2,5)   p = 0. Initiating the computations via H (2,k) p = 0, k 5 would not yield the described result.
The Riccati system (1.1) and the generalized differential operator (4.5) are symmetrical with respect to the coefficients a k , b k (k = 0, 1, 2) and γ 1 , γ 2 .Thus, the computations presented above can be repeated for H (2)  q with results that are obtained by replacing a k with b k , γ 1 with γ 2 and vice versa in (4.20) and (4.22).
All explicit analytical expressions obtained during the application of the procedure presented above are given in the supplementary file one (accessed via the GitHub repository https://bit.ly/3BuIIo3).
The results listed above are summarized in the following lemma.Let the conditions (4.15) and (4.20)-(4.23)hold true, then the sequences p j ; j = 0, 1, . . .and q j ; j = 0, 1, . . .have a single zero characteristic root each (see the last paragraph of the Section 2.4): The remaining characteristic roots ρ 2 and ν 2 are computed from the following characteristic equations: Solving (4.25) yields: Relations (2.14) and (2.17) yield the following systems of linear equations with respect to λ 1 , λ 2 and µ 1 , µ 2 : (4.27) Note that 0 0 is taken to be equal to one [29].Solutions to (4.27)-(4.28)read: Thus, kink solitary solutions to (3.1) read: Rearranging (4.30), (4.31) yields: where the parameters read: Since x 1 = cx 1 , y 1 = cy 1 , t 1 = ct 1 , the parameters x 1 , y 1 , t 1 of kink solutions x, y read: where Then, kink solitary solutions to the original system of differential equations (1.1) can be written as: where k = 1, 2. Analogously, equalities in (4.4) are written in the form representing functional relations between parameters of the system and parameters of the solitary solutions: Note that parameters x 1 , y 1 , t 1 are related to x 1 , y 1 , t 1 through relationships Furthermore, solving (4.4) with respect to σ 1 , σ 2 , x 1 , y 1 yields: The parameters of (1.1) and the solitary solution can be grouped into two categories: • Parameters of the system of differential equations (1.1): γ 1 , γ 2 , a l , b l (l = 0, 1, 2).
According to Lemma 4.1, system (5.1)admits kink solitary solutions if, and only if, the following constraint holds true: where α and β are computed using the stepwise computational scheme outlined in section three.All explicit analytical expressions obtained during the application of this procedure are given in the supplementary file two (accessed via the GitHub repository https://bit.ly/3BuIIo3).
The analytical expressions of kink solitary solutions to (5.1) are then obtained using the results presented in Section 4.3 as follows: (5.4) (5.5) The system (5.1) is integrated using numerical techniques; results are displayed in Figure 4. Part (a) of Figure 4 illustrates the phase portrait of the system, where solitary and non-solitary solutions to (5.1) correspond to solid red and black lines respectively.Note, that the red line depicts the constraint (5.The validity of constraint (5.3) can be verified by the following computational experiment.Let x( jh, u, v) and y( jh, u, v) denote the approximate numerical solutions to (5.1) at c = 0, obtained using the constant step numerical integrator with the step-size h = 0.1 ( j = 0, . . ., 100).The difference between the approximate numerical solution and the kink solitary solution (5.4)-(5.5) is then defined as follows: x( jh, u, v) − x( jh, 0, u) + y( jh, u, v) − y( jh, 0, v) . (5.6) The distribution of ε (u, v) is displayed in Figure 5.Is is clear that the values of ε (u, v) are closest to zero on the line v = 1 2 u − 1 8 .

Conclusions
Kink solitary solutions to a system of Riccati differential equations with diffusive coupling were constructed via the generalized differential operator technique aided by computer algebra computations.The presented scheme was also used to derive necessary and sufficient existence conditions for kink solitary solutions with respect to the Riccati system parameters.This approach allowed us to determine analytical conditions that generate the constraints that the solitary solutions must satisfy in both the space of system parameters and the space of solution parameters.Furthermore, it was proven via the inverse balancing technique that (1.1) cannot admit bright/dark and higher-order solitary solutions under any conditions.
The presented stepwise approach (Section 4.2) is a powerful tool for the construction of analytical solutions to nonlinear differential equations.Note that this technique is not limited to system (1.1), but can potentially be applied to a variety of both Riccati-type equations and more general ordinary differential equations (ODEs) with polynomial nonlinearity.
The obtained results are counterintuitive.Even though diffusive coupling can be perceived as less limiting than multiplicative coupling, the multiplicative case yields both kink and bright-dark solitary solutions, while diffusive coupling results only in kink solitary solutions.
Solitary solutions have special properties that make them especially important in the analysis of propagating waves.Problems concerning propagating waves arise from partial differential equations (PDEs), which are transformed to ODEs via a linear wave-variable substitution.Analysis of these

.10) 2 . 4 .
The construction of closed-form solitary solutions 2.4.1.Linear recurring sequences for the construction of solitary solutions Properties of linear recurring sequences exploited for the determination of the existence of closedform solitary solutions are discussed in this section.

Figure 2 .
Figure 2. A schematic diagram of the procedure used for the construction of solitary solutions described in Sections 2.1-2.4.

4 .
The construction of solitary solutions to (1.1) 4.1.The derivation of necessary existence conditions of kink solitary solutions to(1.1)

Figure 3 .
Figure 3.A schematic diagram of the stepwise computational scheme for the derivation of solitary solutions.

Figure 4 .
Figure 4. Numerical integration of the system (5.1).Part (a) illustrates the phase portrait of the system.Solid red and black lines correspond to solitary and non-solitary solutions respectively.Blue circle denotes equilibrium point.Fixing v = −0.025yields solitary and non-solitary solutions to (5.1) depicted in the part (b).Analogously, Fixing u = 0.2 yields solitary and non-solitary solutions to (5.1) depicted in the part (c).

Figure 5 .
Figure 5.The distribution of ε (u, v) for the system (5.1) at c = 0.The linear relationship between u and v is v = 1 2 u − 1 8 coincides with the curve on which ε (u, v) = 0.
Note that (4.20) and (4.22) lead to the following functional relations: