Investigation of Electromagnetic Wave Structures for a Coupled Model in Anti-ferromagnetic Spin Ladder Medium

The article studies the extraction of electromagnetic wave structures in a spin ladder anti-ferromagnetic medium with a coupled generalized non-linear Schrodinger model. The direct algebraic technique is used to extract the wave solutions. The solutions are obtained in the form of dark, singular, kink, and dark-singular under different constraint conditions. Moreover, the dynamic behavior of the structures have depicted in 3D graphs and their corresponding counterplots. The results are helpful for the understanding of wave propagation study and are also vital for numerical and experimental verifications in the field of electromagnetic wave theory.

The magnetic moments of molecules and atoms, normally linked to the spins of electrons, adjusted in an ordinary pattern with neighboring spins spell in reverse directions. This is like ferrimagnetism and ferromagnetism, a manifestation of ordered magnetism. Theoretical and mathematical theories argue that the spin ladder system is an excellent medium through which the interaction between different spins can be mapped to an approximate Heisinberg-type coupling with a coupling parameter that is inversely proportional to the distance between two separated spins. The dynamics of electromagnetic solitons with the coupled model in an anti-ferromagnetic spin ladder medium is of great interest among researchers due to its variety of applications. The spin ladder systems are a great source with which to develop significant interest in both experimental and theoretical points of view. Anti-ferromagnets test different ideas that involve a strong correlated system [6][7][8][9]. Spin ladders have many applications in different fields of quark physics, superconductors, and ultra-cold atoms, etc. The study of anti-ferromagnetic is still in its early stages. In anti-ferromagnets, the staggered magnetization variable M contributes the first derivative in this manner, i.e., (∇M). This is because, in the case of anti-ferromagnets, we have taken limits in specified intervals for two different sublattices individually. The parent cuprate insulators are the best example of anti-ferromagnets with isotropic and predominantly nearest-neighbor coupling. They satisfied the theory that gives an ordered ground state for (s = 1 2 AF) by showing simple long range anti-ferromagnetic order at low temperature [19][20][21][22][23][24][25][26].
In this article, a coupled generalized derivative nonlinear Schrödinger that describes the dynamical behavior of electromagnetic waves in a spin ladder antiferromagnetic medium system is considered and under investigation. The Heisenberg model was studied in Chen et al. [15] and Xu et al. [16] and considered with an anisotropic spin ladder and two ferromagnetic lattices. This lattices consist of N spins and are directed in the same direction? For more details see also Kavitha et al. [26]. The system is read as where q j for j = 1, 2, and n = 3 − j, are the wave profiles, and α k for k = 1, 2, · · · , 5, represent the real coefficients and are defined by Where j a is ferromagnetic spin exchange interaction, j b is antiferromagnetic coupling, j c is the exchange coupling, and A represents the single-ion uniaxial anisotropy. The last equation is reduced to generalized coupled derivative non-linear Schrodinger system by considering α m = 0, for m = 2, · · · , 5. In the following section, the considered model is analyzed.

THE MODIFIED DIRECT ALGEBRAIC METHOD
The section studies, MDAM [27] to investigate the wave structures of NLPDs. Thus, we consider NLPDs in following form: where q is a profile of wave structure and H is called a polynomial of q and its partial derivatives along with non-linear terms.
To extract wave structures, the method is followed by using the steps as discussed under.
Step 1: First, the NPDEs is converted into non-linear ODEs using the following transformation.
q(x, t) = U(ξ ), and ξ = B(x − wt), where B and w are arbitrary parameters. It allows us to reduce the above equation in an ODE of U and have the form Step 2: It is supposed that the solution of above equation satisfies the following ansatze: where γ is a parameter and its value is determined, Step 3: The homogeneous balance technique is followed where the highest order derivative is balanced with non-linear terms, to find the value of m, and where m ∈ Z + .
Step 4: The use of the above equation and collecting the terms of the same order of ϕ j together. Equate each term of ϕ j to zero, which produce the system of algebraic equations.
Step 5: The solution of the system of algebraic equations along with the following wave structures are general solutions.
In the following section, the exact wave structures of Equation (1.1) can be obtained.

ANALYTICAL ANALYSIS
We consider the complex transformation q j (x, t) = U j (ξ ) × e iφ , where ξ = B(x − νt) and φ = −kx + wt + θ . It reduces the partial differential equation to an ordinary differential equation. After some mathematical work, the following real part of Equation The imaginary part equation gives the constraint condition for α 2 (α 3 + α 5 ) > 0. To find the solution of Equation (1.1), let us consider the solution (see also [18]), where Z ′ = γ + Z 2 , a ′ i s, b ′ i s, and γ are the real parameters. The parameters are to be determined later. It is also noted that Z = Z(ξ ), and so does Z ′ = dZ dξ . To investigate the electromagnetic waves of the system, we find the solution of Equation (3.2) by finding the homogenous balance m = 1 between the non-linear term and highest order derivatives present in this equation. We have the following value of U after substituting the homogenous balance: where a 0 , a 1 and b 1 are real parameters. To calculate real parameters, we put U and required derivatives in Equation (3.2). After simplification and by equating the coefficients of same power of Z, the system of equations is obtained. To get the values of parameters and the solutions against these different parameters, we solve this system by using Maple. Thus, different cases along with solutions are discussed below.

For case 1:
The values of parameters are and the corresponding dark and singular wave structures can be obtained for different values of γ . The constraint condition, for the existence of these solutions, is given by For γ < 0, one may have the following wave structures of Equation (1.1) The following two cases are obtained from the above solution and considered as the diagonal components of the spin ladder.
For γ > 0, one may have the following periodic solutions. and For γ = 0, one may have the following periodic solutions where, The graphical representations and contour plots of the solutions for q 1 to q 4 are shown in Figure 1 for different values of parameters α 1 = 1, r = 1.25, B = 5, k = 0.1, ξ = 0.01, θ = 0.2, and b = 0.5.

For case 2:
The values of the parameters are and the corresponding combined dark-singular wave structures are constructed. For γ < 0, the following type of exact solutions of Equation (1.1) is written: The following two cases are obtained from the above solution and considered as the diagonal components of the spin ladder.
and we also have For γ > 0, one may have the following periodic solutions and For γ = 0, one may have the following periodic solutions where, The pattern of the solutions for q 6 to q 9 are shown in Figure 2 for the values of parameters α 1 = 0.007, α 2 = 0.98, α 3 = 0.01, r = 0.76, B = 0.98, k = 0.98, ξ = 0.01, θ = 0.2, and b = 1.5.

For case 3: The values of parameters are
and the corresponding dark and singular wave structures can be obtained. For γ < 0, the following forms of the exact solutions to Equation (1.1) are obtained. and The following two cases are obtained from the above solution and considered as the diagonal components of the spin ladder.
For γ > 0, one may have the following periodic solutions and For γ = 0, one may have the following periodic solutions where, The pattern of the solutions for q 11 to q 14 are shown in Figure 3 for the values of parameters α 1 = 2, r = 1.5, B = 3.9, k = 0.98, ξ = 0.01, θ = 0.2, and b = 0.25.

For case 4: The values of parameters are
and the corresponding dark and singular wave structures can be obtained. For γ < 0, the following exact solutions to Equation (1.1) are obtained.
The following two cases are obtained from the above solution and are considered the diagonal components of the spin ladder.
For γ > 0, one may have the following periodic solutions Frontiers in Physics | www.frontiersin.org and For γ = 0, one may have the following periodic solutions The pattern of the solutions for q 16 to q 19 are shown in Figure 4 for the values of parameters α 1 = 0.002, r = 0.5, B = 0.9, k = 0.98, ξ = 0.01, θ = 0.2, and b = 0.25.

For case 5:
The values of parameters are , and the corresponding dark and singular wave structures are obtained. For γ < 0, one can obtain the following exact solutions to Equation (1.1).
The following two cases are obtained from the above solution and considered as the diagonal components of the spin ladder.
For γ > 0, one may have the following periodic solutions For γ = 0, one may have the following periodic solutions These are the new solitons and periodic wave structures.

CONCLUSIONS
The article gives single and combined electromagnetic wave structures for the coupled non-linear Schrödinger equations along with the coefficients of ferromagnetic spin exchange interaction, antiferromagnetic coupling, exchange coupling, and single-ion uniaxial anisotropy. The model under investigation describes the dynamic behavior of electromagnetic waves in a spin ladder antiferromagnetic medium. First the complex transformation is used and then modified extended direct algebraic method is utilized to find dark, singular, and dark-singular wave structures. Some other solutions (singular periodic) are also fall out during the analytical analysis. The constraint conditions for the existence of wave structures for different parameters are also observed. Moreover, the 3D plots and corresponding contour plots of the real part of solutions are drawn by choosing suitable parameters. It is also observed that the method used is effective, powerful, reliable, and much more practical in obtaining the exact wave structures for non-linear phenomena that arise in fields like telecommunication engineering, mathematical

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