New and More Solitary Wave Solutions for the Klein-Gordon-Schrödinger Model Arising in Nucleon-Meson Interaction

This paper considers methods to extract exact, explicit, and new single soliton solutions related to the nonlinear Klein-Gordon-Schrödinger model that is utilized in the study of neutral scalar mesons associated with conserved scalar nucleons coupled through the Yukawa interaction. Three state of the art integration schemes, namely, the e−Φ(ξ)-expansion method, Kudryashov's method, and the tanh-coth expansion method are employed to extract bright soliton, dark soliton, periodic soliton, combo soliton, kink soliton, and singular soliton solutions. All the constructed solutions satisfy their existence criteria. It is shown that these methods are concise, straightforward, promising, and reliable mathematical tools to untangle the physical features of mathematical physics equations.


INTRODUCTION
Many of the problems arising in mathematically-oriented scientific fields such as physics and engineering are described by partial differential equations (PDEs). PDEs are used to depict an ample variety of phenomena such as dislocations in crystals, superconductivity, laser pulses in two-phase [1,2], waves in ferromagnetic materials, and many more [3][4][5][6]. Many theories such as electromagnetism, diffusion, fluid flow, etc. are presented to understand the dynamics of PDEs [7,8]. Therefore, exploring exact solutions for PDEs plays an important role in such fields. These solutions might be essential and important for exploring some physical phenomena. The majority of PDEs are not exactly solvable with existing mathematical techniques. Especially for higher order nonlinear PDEs, existing methods are not able to find exact solutions. However, due to the invention of algebraic system solvers such as Mathematica and Maple, many integrating schemes have been proposed, such as the e − (ξ ) -expansion method, Hirota's bilinear method, the homogeneous balance reduction of the PDE to a quadrature problem, the truncated Painlevé expansion, etc. [9][10][11][12][13].
Solitons are formed because of an interplay between nonlinear and dispersive effects. The importance of such waves lies in their roles in telecommunication systems as well as other physical sciences such as nonlinear optics, acoustics, convictive fluids, condensed matter, and solid-state and plasma physics [14][15][16][17][18][19].
The main successes of quantum mechanics in the quantitative description of non-relativistic systems are connected with the Schrödinger equation. Schrödinger used the Klein-Gordon equation (KGE) to model a system of motion of massive spinless particles in a quantum field hypothesis [20,21]. The quantum wave model is thought of as the non-relativistic limitation of the KGE. The KGE is a second order differential equation in both spatial and temporal coordinates, reduced to two coupled first order differential equations. To obtain explicit exact solutions, various techniques have been implemented such as the Jacobi elliptic expansion method, the mapping method, and the Fexpansion method but these methods produce very complex solution expressions [22][23][24].
The nonlinear Klein-Gordon-Schrödinger (KGS) framework depicts the association of neutral scalar mesons connecting with scalar nucleons. This model depicts complex processes and has attracted the consideration of researchers from different fields. Many authors examined the behavior of solutions by employing numerical and analytical techniques. Solitary wave solutions are studied in [20], the behavior of the equations is evaluated by a modified decomposition method in [25], Biswas and Triki [26] scrutinize the KGS model alongside power law nonlinearity to attain the solution in the form of solitons, the Chebyshev pseudo spectral multidomain strategy was taken into consideration for the mathematical solution of the given system in [27], Yumak et al. analyzed the exact periodic and solitary wave polynomial solutions of nonlinear KGS equations in [28], a high-order compact finite difference technique is examined for a governing model in [29]. In addition, various techniques are applied in [30,31].
The aim of this paper is to establish some novel and widely applicable traveling wave solutions of the non-linear KGS equation through effective methods such as the tanh-coth expansion method, e − (ξ ) -expansion method, and Kudryashov's method. This model depicts the scalar nucleons associating with neutral scalar mesons linked by the Yukawa potential in quantum field theory [32,33]. These methods are used to obtain a new explicit solution of the KGS system which can be useful to study the physical nature of many nonlinear phenomena in a number of fields such as modern physics, fluid dynamics, quantum mechanics, and plasma physics. This paper is arranged as follows: In section 2, we proposed the model. In section 3, solutions of the KGS system are formulated through integration schemes and graphical interpretation. Section 4 contains the results and discussion. In section 5, concluding points are expressed.

PROPOSED MODEL
The proposed model [20] has the form This coupled system describes the interplay of a meson field with a nucleon field and is significant in modern physics. Here N = N(x, t) is a meson field, W = W(x, t) is a complex scalar nucleon field, and c is a real constant.
The following wave transformation is applied, where ξ = x − αt, to obtain traveling wave solutions for the proposed model given by Equations (1) and (2). In the above transformation the wave number is ω, k is the frequency, and α is the velocity of the soliton. Plugging Equation (3) into Equations (1) and (2), then equating the real parts gives The imaginary part of Equation (1) gives the velocity of soliton, i.e., α = −2k. Solving Equation (4) for n(ξ ), we get After plugging the value of n(ξ ) in Equation (5), we obtain the following ODE as In accordance with the e − (ξ ) -expansion scheme [34], the solution of Equation (7) has the form The homogenous balance method gives N = 1. For N = 1, Equation (8) becomes here (ξ ) is the solution of the following ODE (11) and (12) with suitable choose of values of parameters as l = 3, s = 1, m = 2, k = 1, ξ = x + 2t and C = 0.
By substituting Equation (9) into Equation (7) a system of equations for a 0 and a 1 is retrieved by comparing the coefficients of e − (ξ ) equal to zero. By finding unknowns a 0 and a 1 from the obtained system and inserting them into Equation (9), solutions of the coupled KGS Equations (1) and (2) are obtained.
The obtained results are summarized in the following sets. SET 1 Soliton solutions for Set 1 are calculated. When s > 0 and m = 0, then and The Figure 1, depicts the graphical representation of the absolute values of the obtained soliton solutions given in Equation (11) and Equation (12). In Figure 1A represents the kink soliton and Figure 1B represents the bright soliton. When s < 0 and m = 0, then and The Figure 2, depicts the graphical representation of the absolute values of the obtained soliton solutions given in Equation (13) and Equation (14). In Figure 2A represents the periodic soliton and Figure 2B represents the singular soliton. When s > 0 and m = 0 and l = 0 then and  (13) and (14) with values of parameters as l = 1, s = − 3, m = 1, k = 1, ξ = x + 2t and C = 0.
When s = 0 and m = 0 and l = 0, then and When s = 0 and m = 0 and l = 0, then and where C is the constant of integration and s = l 2 − 4m. Soliton solutions for Set 2 are calculated. When s > 0 and m = 0, then provided that m < 0.
provided that m < 0. When s < 0 and m = 0, then provided that m > 0.
provided that m > 0, where C is the constant of integration and s = l 2 − 4m.
Substituting Equation (26) into Equation (7), a system of nonlinear equations for a 0 , a 1 , and b 1 is obtained by the comparison of the coefficients of tanh(mξ ) to zero. Upon solving the obtained system for a 0 , a 1 , and b 1 and plugging them in Equation (26), the solutions of the coupled KGS Equations (1) and (2) are obtained. The obtained results are summarized in the following sets. SET 1

SET 3
Singular soliton solutions relative to SET 1 are obtained as and A dark-singular combo soliton solution relative to SET 2 is attained as and a singular soliton is calculated as The Figure 3, depicts the graphical representation of the absolute values of the obtained soliton solutions given in Equation (32) and Equation (33). In Figure 3A represents the dark-singular combo soliton and Figure 3B represents the singular soliton.
A dark soliton solution relative to SET 3 is attained as and bright soliton is calculated as According to Kudryashov's method [36], the predicted solution of Equation (7) has the following form The homogenous balance method gives N = 1. For N = 1, Equation (36) becomes Inserting Equation (37) into ODE Equation (7), an algebraic system of equations for a 0 and a 1 is obtained by equating every coefficient of different powers of 1 1+d(cosh(ξ )+sinh(ξ )) to zero. The obtained system is solved for a 0 and a 1 , and replacing these values in Equation (37) gives solutions of the coupled KGS Equations (1) and (2). The following solution set arises SET 1 A kink soliton solution is given as and a bright soliton solution is obtained as The Figure 4, depicts the graphical representation of the absolute values of the obtained soliton solutions given in Equation (38) and Equation (39). In Figure 4A represents the kink soliton and Figure 4B represents the bright soliton.

Novelty of the Results
It is worth mentioning here that the proposed model has been solved for the first time by the e − (ξ ) -expansion method, tanhcoth expansion technique, and Kudryashov's method to extract solitonic structures. The results presented in this piece of research could be very useful in discussing the physical properties of the different nonlinear evolution equations emerging in quantum mechanics, fluid dynamics, and plasma physics. The solitonic structures obtained in this study could attract the attention of researchers working in the field of optical fiber communication systems. The comparison of our results, with the outcomes of [20,21], show that bright and dark solitons as well as dipole soliton, singular soliton, and kink soliton solutions have been found for the first time in this article.

RESULTS AND DISCUSSION
It is important to clarify that the analytical methods utilized in this article are truly state of the art techniques for extracting the soliton solution of the non-linear Klein-Gordon-Schrödinger model. It is important to note here that each integration method has its own benefits and disadvantages compared to other accessible strategies. For example, the inverse scattering method is not useful for log-law, power law, and dualpower law nonlinearities. Only bright solitons are recovered by the semi-inverse variational algorithm. Likewise, here, the e − (ξ ) -expansion technique gives bright soliton, kink soliton, periodic soliton, and singular soliton solutions. The second method applied in this research extracts dark soliton, darksingular combo soliton, and singular soliton solutions. The third method employed here obtains bright soliton and kink soliton solutions.

CONCLUSION
In this article, new soliton solutions have been obtained by utilizing three well-known integration architectures namely, the tanh-coth expansion strategy, Kudryashov's strategy, and the e − (ξ ) -expansion strategy. To the best of our knowledge, these fresh examples of soliton solutions have been obtained for the first time for the KGS model. Since the invention of symbolic computation tools, the solution procedures have been simplified, and therefore the described methods are becoming more efficient in solving many physical problems. The outcomes of this paper consist of dispersive solitons incorporating CQS and cubic nonlinearities. Kudryashov's method along with the generalized tanh method extract singular, bright, singular periodic, and a combo type of solitons for the given model. The advantage of these techniques is quite evident as they have no limitations in finding such wave profiles.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
All authors contributed equally to the writing of this paper and read and approved the final version of the manuscript.