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–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–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–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 F-expansion method but these methods produce very complex solution expressions [22–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.

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.

In this section, three state of the art integration schemes, Raza et al. [34], Asokan and Vinodh [35], and Ullah et al. [36] are employed to extract bright soliton, dark soliton, dipole and combo soliton, kink soliton, and singular soliton solutions.

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

By substituting Equation (9) into Equation (7) a system of equations for a₀ and a₁ is retrieved by comparing the coefficients of e^−Φ(ξ) equal to zero. By finding unknowns a₀ and a₁ 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

SET 2

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

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² − 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² − 4m.

According to the the tanh-coth method [35], the estimated solution of Equation (7) has the form

The homogenous balance method, gives N = 1. For N = 1, the above equation takes the form

Substituting Equation (26) into Equation (7), a system of nonlinear equations for a₀, a₁, and b₁ is obtained by the comparison of the coefficients of tanh(mξ) to zero. Upon solving the obtained system for a₀, a₁, and b₁ 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 2

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₀ and a₁ is obtained by equating every coefficient of different powers of 11+d(cosh(ξ)+sinh(ξ)) to zero. The obtained system is solved for a₀ and a₁, 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.

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 dual-power 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, dark-singular combo soliton, and singular soliton solutions. The third method employed here obtains bright soliton and kink soliton solutions.

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.

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

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

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.