Next Article in Journal
The Most Refined Axiom for a Digital Covering Space and Its Utilities
Previous Article in Journal
Three-Dimensional Volume Integral Equation Method for Solving Isotropic/Anisotropic Inhomogeneity Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solutions of the Two-Wave Interactions in Quadratic Nonlinear Media

1
Department of Mathematics, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), P.O. Box 90950, Riyadh 11623, Saudi Arabia
2
Department of Physics, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), P.O. Box 90950, Riyadh 11623, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(11), 1867; https://doi.org/10.3390/math8111867
Submission received: 11 September 2020 / Revised: 17 October 2020 / Accepted: 20 October 2020 / Published: 26 October 2020

Abstract

:
In this paper, we propose a reliable treatment for studying the two-wave (symbiotic) solitons of interactions in nonlinear quadratic media. We investigate the Schauder’s fixed point theorem for proving the existence theorem. Additionally, the uniqueness solution for this system is proved. Also, a highly accurate approximate solution is presented via an iteration algorithm.

1. Introduction

In the simplest case of type-I second harmonic generation (SHG) without walk-off between harmonic waves soliton evolution is described by the normalized system [1,2,3,4,5,6]:
i φ z + r 2 φ x 2 φ + φ * ψ = 0 , i σ ψ z + s 2 ψ x 2 α ψ + 1 2 φ 2 = 0 ,
where α is the rescaled soliton parameter and satisfies α = σ ( 2 β + λ ) / β , the dimensionless parameter β is the normalized nonlinearity-induced shift to the propagation constant of the fundamental harmonic wave, σ and λ are the coefficient and phase mismatch parameter, respectively. This system represents the generic model of χ 2 solitons. There are other types of normalization also used in the literature see e.g., [7,8,9]. The solutions of this system have been discussed for one dimensional in [3,4,5,10,11,12,13] and multi-dimensional cases in [14]. Under suitable assumptions, the problem of the two-wave (symbiotic) solitons can be reduced to the solution of the following coupled system [1,2,3,4,6,13,15]
r d 2 φ d x 2 φ + φ ψ = 0 , s d 2 ψ d x 2 α ψ + 1 2 φ 2 = 0 ,
The properties of solitons described by system (2) are well known see [3,4,5], where the authors Buryak and Kivshar [3,4,5] looked for stationary (i.e., z-independent) localized solutions of the normalized system in the form of an asymptotic series in the parameter α 1 and found the real functions φ ( x ) and ψ ( x ) in the form of asymptotic series:
ϕ ( x ) = 2 α 1 2 sech ( x ) + 4 s α 1 2 tanh 2 ( x ) sech ( x ) + ,
ψ ( x ) = 2 sech 2 ( x ) + s α 1 16 sech 2 ( x ) 20 sech 4 ( x ) + ,
for bright solitons at r = + 1 , and
ϕ ( τ ) = 2 α 1 2 tanh ( τ ) 2 s α 1 2 τ sech 2 ( τ ) tanh ( τ ) sech 2 ( τ ) + ,
ψ ( x ) = tanh 2 ( τ ) + s α 1 2 2 τ tanh ( τ ) sech 2 ( τ ) 4 sech 2 ( τ ) + 5 sech 4 ( τ ) + ,
where τ = x 2 , for dark solitons at r = 1 .
Exact solutions of system (2) have been found at α = 1 for r = s = 1 , r = s = 1 and r = s = 1 in [10,13]. Another solution of (2) in the case r = 1 , s = 1 and α = 2 is provided in an explicit analytical form [11]. Also, different analytical approximation methods have been proposed to deal with the system (2). For example, an accurate approximate solution with the help of the variational method is obtained in [13] and a family of bright (dark) solitons for s = r = 1 , s = r = 1 and α > 0 was also discussed in [15] using the numerical shooting and relaxation techniques.
Recently, the authors [16] studied this coupled system subject to the following boundary conditions:
ϕ ( l 1 ) = ϕ ( l 2 ) = 0 and ψ ( l 1 ) = ψ ( l 2 ) = 0
and showed that the solutions of this system could be obtained numerically using the Green’s function method.
Since the exact analytical solutions of Equation (2) cannot be found for arbitrary values of α , then the purpose of this work is to present a result of the existence and uniqueness of solutions. Also, an exact implicit solution is derived using a useful procedure at α = 1 . Thus, the paper is organized as follows: in Section 2, we investigate the existence and uniqueness theorem of the two-wave solitons in quadratic media, where the problem is formulated in the context of two nonlinear coupled differential equations in one dimension. Then, in Section 3, we solve the coupled system by an appropriate technique with suitable boundary conditions. A systematic numerical procedure is proposed in Section 4. Finally, we conclude with some remarks in Section 5.

2. An Existence and Uniqueness Theorem

Rewrite (2) in the following system
d 2 φ d x 2 = f 1 ( φ , ψ ) , d 2 ψ d x 2 = f 2 ( φ , ψ ) ,
where f i : [ l 1 , l 2 ] × R × R R , i = 1 , 2 are defined by f 1 ( φ , ψ ) = 1 r ( φ φ ψ ) and f 2 ( φ , ψ ) = 1 s ( α ψ 1 2 φ 2 ) .
Existence-In this section, we shall deal with the existence of solutions of the BVP (2) with (7). First off all, we shall prove the following lemmas, which are useful tools in the proof of the existence and uniqueness theorem.
Lemma 1.
If we assume that φ , ψ C [ l 1 , l 2 ] , l 2 > l 1 . Then, the functions f i ( φ , ψ ) , i = 1 , 2 are Lipschitz continuous functions of φ and ψ .
Proof. 
From the definition of f i ( φ , ψ ) , we have
f 1 ( φ 1 , ψ 1 ) f 1 ( φ 2 , ψ 2 ) = 1 r ( φ 1 ϕ 2 ) + ψ 2 ( φ 2 φ 1 ) + φ 1 ( ψ 2 ψ 1 ) .
Since φ i , ψ i C [ l 1 , l 2 ] , then there exist M i , M i * > 0 , i = 1 , 2 such that φ i ( x ) M i and ψ i ( x ) M i * , M i , M i * > 0 , i = 1 , 2 for all x [ l 1 , l 2 ] .
f 1 ( φ 1 , ψ 1 ) f 1 ( φ 2 , ψ 2 ) 1 r φ 2 φ 1 + M 2 * φ 2 φ 1 + M 1 ψ 2 ψ 1 .
Hence,
f 1 ( φ 1 , ψ 1 ) f 1 ( φ 2 , ψ 2 ) L 1 , 1 φ 2 φ 1 + L 1 , 2 ψ 2 ψ 1 .
Similarly, we obtain
f 2 ( φ 1 , ψ 1 ) f 2 ( φ 2 , ψ 2 ) L 2 , 1 φ 2 φ 1 + L 2 , 2 ψ 2 ψ 1 ,
where L 1 , 1 = 1 r ( 1 + M 2 * ) , L 1 , 2 = M 1 r , L 2 , 1 = M 1 + M 2 2 s and L 2 , 2 = α s . □
Lemma 2.
(See pp. 70–71 [17]). Let g : [ a , b ] R be a continuous function. The unique solution u of the following boundary value problem
u = g ( x )
subject to the Dirichlet boundary conditions u ( a ) = u ( b ) = 0 is given by
u ( x ) = a b G ( x , y ) g ( y ) d y ,
where G ( x , y ) is the Green function given by
G ( x , y ) = 1 b a ( x a ) ( b y ) , a x y b , 1 b a ( y a ) ( b x ) , a y x b
and a b G ( x , y ) d y ( b a ) 2 8 .
Replacing g ( x ) , a and b by f i ( φ , ψ ) ,   l 1 and l 2 in Lemma 2, respectively, we obtain an equivalent integral system
φ ( x ) = l 1 l 2 G ( x , s ) f 1 ( φ ( s ) , ψ ( s ) ) d s , x [ l 1 , l 2 ] , ψ ( x ) = l 1 l 2 G ( x , s ) f 2 ( φ ( s ) , ψ ( s ) ) d s , x [ l 1 , l 2 ] .
Define the Banach space X = C [ l 1 , l 2 ] with norm ( φ , ψ ) = φ + ψ , where u = max l 1 x l 2 u ( x ) and the operator T : X X by T ( φ , ψ ) = T 1 ( φ , ψ ) , T 2 ( φ , ψ ) , where
T 1 ( φ , ψ ) = l 1 l 2 G ( x , s ) f 1 ( φ ( s ) , ψ ( s ) ) d s
and
T 2 ( φ , ψ ) = l 1 l 2 G ( x , s ) f 2 ( φ ( s ) , ψ ( s ) ) d s .
Since φ ( x ) M and ψ ( x ) M * , x [ l 1 , l 2 ] , consider the closed and convex set
S = ( φ , ψ ) X : ( φ , ψ ) M + M * .
Furthermore, assume that l 2 l 1 8 r s ( M + M * ) s ( M + M M * ) + r ( α M * + M 2 2 ) 3 .
In the theory of differential equations, there are a lot of methods to establish the existence of solutions. Theorems concerning the existence and properties of fixed points are known as fixed-point theorems. Such theorems are the most important tools for proving the existence and uniqueness of the solution. The fundamental theorem used in this theory is Schauder’ s theorem. In order to make use of this theorem, it is sufficient to prove the following lemma.
Lemma 3.
For any ( φ , ψ ) S , T ( φ , ψ ) is contained in S .
Proof. 
It follows by the definition of T ( φ , ψ ) that
T 1 ( φ ( x ) , ψ ( x ) ) l 1 l 2 G ( x , s ) f 1 ( φ ( s ) , ψ ( s ) ) d s ( l 2 l 1 ) 3 8 r ( M + M M * ) ,
T 2 ( φ ( x ) , ψ ( x ) ) l 1 l 2 G ( x , s ) f 2 ( φ ( s ) , ψ ( s ) ) d s ( l 2 l 1 ) 3 8 s ( α M * + M 2 2 ) .
Thus
T ( φ ( x ) , ψ ( x ) ) ( l 2 l 1 ) 3 8 M + M M * r + α M * + M 2 2 s .
Since l 2 l 1 is defined by the above condition, thus T ( φ ( x ) , ψ ( x ) ) M + M * . On account of the continuity of f i ( φ ( x ) , ψ ( x ) ) , φ and ψ , it follows that T ( φ ( x ) , ψ ( x ) ) is continuous. This shows that T ( φ ( x ) , ψ ( x ) ) is also contained in S .  □
In order to prove that T ( φ ( x ) , ψ ( x ) ) is equicontinuous, it is easy to see from its definition that
T ( φ ( x ) , ψ ( x ) ) T ( φ ( x ) , ψ ( x ) ) K x x , for any x , x [ l 1 , l 2 ] .
where K = ( l 2 l 1 ) 2 8 M + M M * r + α M * + M 2 2 s .
Therefore T is compact by the classical Ascoli lemma, and Schauder’s fixed point theorem yield the fixed point of T . Thus, we have proved:
Theorem 1.
There exists a continuous solution ( φ , ψ ) which satisfies system Equations (2) and (7) with the condition on l 1 and l 2 .
Uniqueness-A uniqueness theorem can also be obtained from the Lipschitz continuous in φ and ψ .
Theorem 2.
If max max ( L 1 , 1 , L 2 , 1 ) , max ( L 1 , 2 , L 2 , 2 ) < 8 ( l 2 l 1 ) 2 , then the system Equation (2) with (7) has a unique solution ( φ ( x ) , ψ ( x ) ) .
Proof. 
Let ( φ 1 , ψ 1 ) and ( φ 2 , ψ 2 ) be two solutions of (2)–(7). Then, for x [ l 1 , l 2 ] ,
φ 2 φ 1 ( l 2 l 1 ) 2 8 L 1 , 1 max l 1 x l 2 φ 2 ( x ) φ 1 ( x ) + L 1 , 2 max l 1 x l 2 ψ 2 ( x ) ψ 1 ( x ) ,
ψ 2 ψ 1 ( l 2 l 1 ) 2 8 L 2 , 1 max l 1 x l 2 φ 2 ( x ) φ 1 ( x ) + L 2 , 2 max l 1 x l 2 ψ 2 ( x ) ψ 1 ( x ) .
Consequently,
φ 2 φ 1 + ψ 2 ψ 1 A φ 2 φ 1 + ψ 2 ψ 1 ,
where A = ( l 2 l 1 ) 2 8 max max ( L 1 , 1 , L 2 , 1 ) , max ( L 1 , 2 , L 2 , 2 ) . We now apply the condition A < 1 to this inequality, we get φ 2 φ 1 = 0 and ψ 2 ψ 1 = 0 . Therefore ( φ 1 , ψ 1 ) = ( φ 2 , ψ 2 ) .  □
The proof is complete.

3. On the Decoupling of the System (2)

In this section, first of all, we are concerned with the norms estimate for the functions φ and ψ when r = s = 1 .
Lemma 4.
Let φ and ψ be two functions in L 2 ( I ) , where I = [ l 1 , l 2 ] . Then
φ 1 < 2 ψ 1 , i f α < 1 , φ 1 = 2 ψ 1 , i f α = 1 , φ 1 > 2 ψ 1 , i f α > 1 ,
where . 1 is the norm defined in the Sobolev space H 0 1 ( I ) by
u 1 2 = l 1 l 2 u 2 ( x ) + u 2 ( x ) d x , w i t h u ( l 1 ) = u ( l 2 ) = 0 .
Proof. 
Multiplying both sides of the first equation of system (2) by φ and integrating from l 1 to l 2 , we obtain
l 1 l 2 φ ( x ) φ ( x ) d x l 1 l 2 φ 2 ( x ) d x + l 1 l 2 ψ ( x ) φ 2 ( x ) = 0 .
From the second equation of system (2), we have
φ 2 = 2 d 2 ψ d x 2 + 2 α ψ .
By substitution into the last term of (29), we obtain
l 1 l 2 φ ( x ) φ ( x ) d x l 1 l 2 φ 2 ( x ) d x 2 l 1 l 2 ψ ( x ) ψ ( x ) d x + 2 α l 1 l 2 ψ 2 ( x ) d x = 0 .
Integrating by parts and taking into account the given boundary conditions, we obtain
l 1 l 2 φ 2 ( x ) d x + l 1 l 2 φ 2 ( x ) d x = 2 α l 1 l 2 ψ 2 ( x ) d x + 2 l 1 l 2 ψ 2 ( x ) d x .
This gives (27). □
Let us now consider the case α = 1 [6,10,13] and in view of Lemma 6 if φ ( x ) = 2 ψ ( x ) , then it may be shown that the two equations of system (2) can be separated into the following nonlinear equation
ψ ( x ) ψ ( x ) + ψ 2 ( x ) = 0 .
The exact solution to Equation (33) follows by simply multiplying both sides of Equation (33) by ψ .
ψ ψ ψ ψ + ψ 2 = 0 ,
which can be written as follows
1 2 d d x ψ 2 1 2 d d x ψ 2 + 1 3 d d x ψ 3 = 0 .
and integrating with respect to x , we obtain
ψ 2 ψ 2 + 2 3 ψ 3 = c 1 ,
where c 1 is an arbitrary constant of integration. Thus
ψ = ± ψ 2 2 3 ψ 3 + c 1 .
In view of ψ ( x ) = d ψ d x , we have
± d ψ ψ 2 2 3 ψ 3 + c 1 = d x .
Consequently,
± d ψ ψ 2 2 3 ψ 3 + c 1 = x + c 2 ,
where c 2 is also a constant of integration.
The LHS of Equation (39) can be evaluated direct from the integrals of irrational functions. Indeed, if we choose c 1 = 0 , then
d ψ ψ 2 2 3 ψ 3 = 2 tanh 1 1 2 3 ψ = x + c 2 .
Since tanh 2 ( x ) = tanh 2 ( x ) . Hence, a simple computation leads to the implicit solution
( φ , ψ ) = 3 2 1 tanh 2 ( 1 2 ( x + c 2 ) ) , 3 2 1 tanh 2 ( 1 2 ( x + c 2 ) ) .
Thus, we have
Lemma 5.
The system (2) can be decoupled without increase the order of the system into the nonlinear Equation (33) when α = 1 . Furthermore, the solution ( φ , ψ ) is given by (41).
In Figure 1, we display the variation of the exact solutions Equation (41) in terms of the independent variables x for different values of the constant of integration c 2 (Equation (39)). it can be seen that this constant of integration shifted left or right the distribution away from the origin with negative or positive values of c 2 , respectively. Also, it does not affect the behavior of the solutions, and the maximum value of the solution remains unchanged. Thus it can be chosen c 2 = 0 .

4. Numerical Analysis

An integral system equivalent to Equation (2) with Equation (7) can be derived. Indeed, integrating (2) twice from l 1 to x and taking into account the boundary conditions φ ( l 1 ) = ψ ( l 1 ) = 0 , we obtain
φ = β ( x l 1 ) + l 1 x l 1 y f 1 ( φ ( s ) , ψ ( s ) ) d s d y , ψ = γ ( x l 1 ) + l 1 x l 1 y f 2 ( φ ( s ) , ψ ( s ) ) d s d y ,
where β = φ ( l 1 ) and γ = ψ ( l 1 ) are unknown constants to be determined from the second boundary conditions φ ( l 2 ) = ψ ( l 2 ) = 0 .
We now construct a sequence of approximation of the solution that converges to the solution. The components ( φ n , ψ n ) can be elegantly determined by setting the recursion scheme
φ 0 = β ( x l 1 ) , ψ 0 = γ ( x l 1 ) , φ n = φ 0 + l 1 x l 1 y f 1 ( φ n 1 ( s ) , ψ n 1 ( s ) ) d s d y , n 1 , ψ n = ψ 0 + l 1 x l 1 y f 2 ( φ n 1 ( s ) , ψ n 1 ( s ) ) d s d y . n 1 .
Theorem 3.
The sequence ( φ n , ψ n ) defined by (43) converges uniformly on I to the unique solution ( φ , ψ ) .
Proof. 
We shall construct an upper bound for ( φ n + 1 , ψ n + 1 ) ( φ n , ψ n ) by induction.
φ 1 φ 0 = l 1 x l 1 y f 1 ( φ 0 ( s ) , ψ 0 ( s ) ) d s d y K 1 ( x l 1 ) 2 2 ! , ψ 1 ψ 0 = l 1 x l 1 y f 1 ( φ 0 ( s ) , ψ 0 ( s ) ) d s d y K 2 ( x l 1 ) 2 2 ! ,
where K 1 = M + M M * r and K 2 = α M * + M 2 2 s . Proceeding in the same manner, we obtain by induction
φ n + 1 φ n K 1 n + 1 ( x l 1 ) n + 2 ( n + 2 ) ! , ψ n + 1 ψ n K 2 n + 1 ( x l 1 ) n + 2 ( n + 2 ) ! .
The two series n = 0 K 1 n + 1 ( x l 1 ) n + 2 ( n + 2 ) ! and n = 0 K 2 n + 1 ( x l 1 ) n + 2 ( n + 2 ) ! are absolutely convergent series. Moreover, these series dominate the two series n = 0 [ φ n + 1 φ n ] and n = 0 [ ψ n + 1 ψ n ] . Hence, by the Weierstrass test, the last two infinite series converge absolutely and uniformly on I . If we consider the m t h partial sum of these series, we see that n = 0 m [ φ n + 1 φ n ] = φ m + 1 φ 0 and n = 0 m [ ψ n + 1 ψ n ] = ψ m + 1 ψ 0 , that is, ( φ n , ψ n ) converges absolutely and uniformly on I . If we now define ( φ , ψ ) = lim n ( φ n , ψ n ) , then taking the limit as n , we obtain
φ = β ( x l 1 ) + l 1 x l 1 y f 1 ( φ ( s ) , ψ ( s ) ) d s d y , ψ = γ ( x l 1 ) + l 1 x l 1 y f 2 ( φ ( s ) , ψ ( s ) ) d s d y .
It follows that upon differentiation of this system that ( φ , ψ ) is the solution of Equation (2). Furthermore, it is clear that φ ( l i ) = ψ ( l i ) = 0 , i = 1 , 2 .  □
In view of (43), the numerical solutions are then given by
φ = β ( x l 1 ) + 1 r β ( x l 1 ) 3 3 ! β γ ( x l 1 ) 4 4 ! + , ψ = γ ( x l 1 ) + 1 s α γ ( x l 1 ) 3 3 ! β 2 2 ( x l 1 ) 4 4 ! + .
If we match ( φ 1 , ψ 1 ) at x = l 2 , then we need to solve
β ( l 2 l 1 ) + 1 r β ( l 2 l 1 ) 3 3 ! β γ ( l 2 l 1 ) 4 4 ! = 0 , γ ( l 2 l 1 ) + 1 s α γ ( l 2 l 1 ) 3 3 ! β 2 2 ( l 2 l 1 ) 4 4 ! = 0 ,
we obtain β = 2 ( 4 ! ) 2 r s [ 1 + ( l 2 l 1 ) 2 3 ! s ] [ 1 + α ( l 2 l 1 ) 2 3 ! s ] ( l 2 l 1 ) 6 and γ = 4 ! r [ 1 + ( l 2 l 1 ) 2 3 ! r ] ( l 2 l 1 ) 3 .
In Figure 2, we present the variation of the numerical solutions Equation (47) against of the independent variables x for different values of the rescaled soliton parameter α .
Now if we insert the solutions of Equation (47) in the second member of Equation (43) and performing the integrals, then using the boundary condition for ( φ 2 , ψ 2 ) at x = l 2 , we can obtain the second solutions φ 2 and ψ 2 . As the mathematical expressions are more cumbersome, we plot the numerical solutions in Figure 3. It is clear from these plots that when we increase the order of the recurrence, the numerical solutions converge rapidly to the exact solutions.

5. Conclusions

This paper is concerned with the treatment of the interaction of two-wave solitons in nonlinear quadratic media. These kinds of problems appear in different applications of nonlinear physics and play a crucial role in the stability problems of solitary waves. The problem is presented within the framework of two coupled nonlinear differential equations, which can be solved numerically with specific boundary conditions. But the generalization to any type of boundary condition constitutes a great challenge.
With concordance to real physical problems, the boundary conditions can be chosen properly. Thus, within this framework, we have proved a theorem of existence and uniqueness for the two-wave solitons in nonlinear quadratic media. Furthermore, we have suggested a useful technique of separation of the coupled system, and we have revealed that the formalism leads to analytic solutions.
Moreover, we have explored an interesting numerical technique, and we have used it to obtain the numerical solutions of the coupled system with suitable boundary conditions. The obtained results are in good agreement with those analytically achieved.
These crucial results open a novel class of investigations, which involve solitary waves with more coupled differential equations and more coupling terms. Some other examples of two or more-coupled solitary waves can be treated with the proposed techniques, and the results will be reported elsewhere.

Author Contributions

Investigation, L.B. and S.B.; Methodology, L.B.; Writing—review and editing, L.B. and S.B.; writing—original draft preparation, L.B. and S.B.; Supervision, S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Buryak, A.V.; Akhmediev, N.N. Internal friction between solitons in near-integrable systems. Phys. Rev. E 1994, 50, 31263133. [Google Scholar]
  2. Buryak, A.V.; Akhmediev, N.N. Influence of radiation on soliton dynamics in nonlinear fibre couplers. Opt. Commun. 1994, 110, 287–292. [Google Scholar]
  3. Buryak, A.V.; Kivshar, Y.S. Spatial optical solitons governed by quadratic nonlinearity. Opt. Lett. 1994, 19, 1612–1614, Erratum in 1995, 20, 1080–1080. [Google Scholar]
  4. Buryak, A.V.; Kivshar, Y.S. Solitons due to second harmonic generation. Phys. Lett. A 1995, 197, 407–412. [Google Scholar]
  5. Buryak, A.V.; Kivshar, Y.S. Twin-hole dark solitons. Phys. Rev. A 1995, 51, R41–R44. [Google Scholar]
  6. Buryak, A.V.; Trapani, P.D.; Skryabin, D.V.; Trillo, S. Optical solitons due to quadratic nonlinearities: From basic physics to futuristic applications. Phys. Rep. 2002, 370, 63–235. [Google Scholar]
  7. He, H.; Werner, M.J.; Drummond, P.D. Simultaneous solitary-wave solutions in a nonlinear parametric waveguide. Phys. Rev. E 1996, 54, 896–911. [Google Scholar]
  8. Skryabin, D.V.; Firth, W.J. Generation and stability of optical bullets in quadratic nonlinear media. Opt. Commun. 1998, 148, 79–84. [Google Scholar]
  9. Peschel, T.; Peschel, U.; Lederer, F.; Malomed, B.A. Solitary waves in Bragg gratings with a quadratic nonlinearity. Phys. Rev. E 1997, 55, 4730–4739. [Google Scholar]
  10. Hayata, K.; Koshiba, M. Multidimensional solitons in quadratic nonlinear media. Phys. Rev. Lett. 1993, 71, 3275. [Google Scholar]
  11. Werner, M.J.; Drummond, P.D. Stongly coupled nonlinear parametric solitary waves. Opt. Lett. 1994, 19, 613–615. [Google Scholar]
  12. Torner, L.; Wright, E.M. Soliton excitation and mutual locking of light beams in bulk quadratic nonlinear crystals. J. Opt. Soc. Am. B 1996, 13, 864–875. [Google Scholar]
  13. Sukhorukov, A.A. Approximate solutions and scaling transformations for quadratic solitons. Phys. Rev. E 2000, 61, 4530–4539. [Google Scholar]
  14. Firth, W.J.; Skryabin, D.V. Optical solitons carrying orbital angular momentum. Phys. Rev. Lett. 1997, 79, 2450. [Google Scholar]
  15. Nikolov, N.I.; Neshev, D.; Bang, O.; Królikowski, W.Z. Quadratic solitons as nonlocal solitons. Phys. Rev. E 2003, 68, 036614. [Google Scholar]
  16. Chen, J.; Ge, J.; Lu, D.; Hu, W. A simple approach to study the boundary-induced trajectory evolution of spatial nonlocal quadratic solitons: Basrd on the Green’s function method. Appl. Math. Lett. 2020, 102, 106–108. [Google Scholar]
  17. Hunter, J.K.; Nachtergaele, B. Applied Analysis; World Scientific Publisching: Singapore, 2001. [Google Scholar]
Figure 1. (Color online) The exact solutions (Equation (41)) for different values of the integration constant c 2 . (a) for c 2 = 0 . (b) for c 2 = 2 and (c) for c 2 = 2 . For other parameters (see text) r = s = 1 and α = 1 . Solid red line φ ( x ) , Dashed black line ψ ( x ) .
Figure 1. (Color online) The exact solutions (Equation (41)) for different values of the integration constant c 2 . (a) for c 2 = 0 . (b) for c 2 = 2 and (c) for c 2 = 2 . For other parameters (see text) r = s = 1 and α = 1 . Solid red line φ ( x ) , Dashed black line ψ ( x ) .
Mathematics 08 01867 g001
Figure 2. (Color online) The numerical solutions of Equation (47) for different values of the rescaled soliton parameter α . (a) for α = 1 . and (b) for α = 0.1 . The other parameters are r = 1 and s = 1 . The solid red line represents the solution ϕ while the black dashed is for ψ .
Figure 2. (Color online) The numerical solutions of Equation (47) for different values of the rescaled soliton parameter α . (a) for α = 1 . and (b) for α = 0.1 . The other parameters are r = 1 and s = 1 . The solid red line represents the solution ϕ while the black dashed is for ψ .
Mathematics 08 01867 g002
Figure 3. (Color online) The numerical solutions of Equation (47) for different values of the integration constant α . (a) for α = 1 . and (b) for α = 0.1 . The other parameters are r = 1 and s = 1 . The solid red line represents the solution φ ( x ) while the black dashed is for ψ ( x ) .
Figure 3. (Color online) The numerical solutions of Equation (47) for different values of the integration constant α . (a) for α = 1 . and (b) for α = 0.1 . The other parameters are r = 1 and s = 1 . The solid red line represents the solution φ ( x ) while the black dashed is for ψ ( x ) .
Mathematics 08 01867 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bougoffa, L.; Bougouffa, S. Solutions of the Two-Wave Interactions in Quadratic Nonlinear Media. Mathematics 2020, 8, 1867. https://doi.org/10.3390/math8111867

AMA Style

Bougoffa L, Bougouffa S. Solutions of the Two-Wave Interactions in Quadratic Nonlinear Media. Mathematics. 2020; 8(11):1867. https://doi.org/10.3390/math8111867

Chicago/Turabian Style

Bougoffa, Lazhar, and Smail Bougouffa. 2020. "Solutions of the Two-Wave Interactions in Quadratic Nonlinear Media" Mathematics 8, no. 11: 1867. https://doi.org/10.3390/math8111867

APA Style

Bougoffa, L., & Bougouffa, S. (2020). Solutions of the Two-Wave Interactions in Quadratic Nonlinear Media. Mathematics, 8(11), 1867. https://doi.org/10.3390/math8111867

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop