Coupled circularly polarized electromagnetic soliton states in magnetized plasmas

The interaction between two co-propagating electromagnetic pulses in a magnetized plasma is considered, from first principles, relying on a fluid-Maxwell model. Two circularly polarized wavepackets by same group velocities are considered, characterized by opposite circular polarization, to be identified as left-hand- or right hand circularly polarized (i.e. LCP or RCP, respectively). A multiscale perturbative technique is adopted, leading to a pair of coupled nonlinear Schrodinger-type (NLS) equations for the modulated amplitudes of the respective vector potentials associated with the two pulses. Systematic analysis reveals the existence, in certain frequency bands, of three different types of vector soliton modes: an LCP-bright/RCP-bright coupled soliton pair state, an LCP bright/RCP-dark soliton pair, and an LCP-dark/RCP-bright soliton pair. The value of the magnetic field plays a critical role since it determines the type of vector solitons that may occur in certain frequency bands and, on the other hand, it affects the width of those frequency bands that are characterized by a specific type of vector soliton (type). The magnetic field (strength) thus arises as an order parameter, affecting the existence conditions of each type of solution (in the form of an envelope soliton pair). An exhaustive parametric investigation is presented in terms of frequency bands and in a wide range of magnetic field (strength) values, leading to results that may be applicable in beam-plasma interaction scenarios as well as in space plasmas and in the ionosphere.

The interaction between two co-propagating electromagnetic pulses in a magnetized plasma is considered, from first principles, relying on a fluid-Maxwell model.Two circularly polarized wavepackets by same group velocities are considered, characterized by opposite circular polarization, to be identified as left-hand-or right hand circularly polarized (i.e.LCP or RCP, respectively).A multiscale perturbative technique is adopted, leading to a pair of coupled nonlinear Schr odinger-type (NLS) equations for the modulated amplitudes of the respective vector potentials associated with the two pulses.Systematic analysis reveals the existence, in certain frequency bands, of three different types of vector soliton modes: an LCP-bright/RCP-bright coupled soliton pair state, an LCP-bright/RCPdark soliton pair, and an LCP-dark/RCP-bright soliton pair.The value of the magnetic field plays a critical role since it determines the type of vector solitons that may occur in certain frequency bands and, on the other hand, it affects the width of those frequency bands that are characterized by a specific type of vector soliton (type).The magnetic field (strength) thus arises as an order parameter, affecting the existence conditions of each type of solution (in the form of an envelope soliton pair).An exhaustive parametric investigation is presented in terms of frequency bands and in a wide range of magnetic field (strength) values, leading to results that may be applicable in beam-plasma interaction scenarios as well as in space plasmas and in the ionosphere.

I. INTRODUCTION
Electromagnetic (EM) pulse propagation in a plasma has been in the focus of researchers for decades.From a basic theoretical point of view, it is associated with fundamental physical processes involving various nonlinear mechanisms and instabilities, which have been attracting increasing interest since the early days of plasma research.Thanks to a series of seminal early studies [1][2][3][4] which set the foundations of contemporary laserplasma interaction science [5,6], the study of electromagnetic pulse propagation has been shown to be relevant with phenomena like parametric instabilities [7,8], harmonic generation [9], self-focusing [10], intense electric and magnetic field [11,12], and wakefield [13,14] generation, among others.As laser power and technological sophistication increased, in the last decades, new regimes were attained where ultrashort ultrastrong electromagnetic (laser) pulses could be realized in the laboratory, thus entering the realm of relativistic dynamics (where the electron quiver velocity is comparable to the speed of light) [15] and, more recently, even quantumelectrodynamics [16][17][18][19].EM beam propagation is nowadays a cutting edge topic in science, due to its relevance with applications including laser-assisted (inertial confinement fusion, ICF) schemes for energy production, to ion acceleration [20] and to sophisticated experimental diagnostic techniques [21] and thus one of the hottest topics in contemporary research.
The generation of relativistic electromagnetic solitons is governed by various physical mechanisms, including dispersion due to the finite electron inertia and nonlinearity due to a relativistic electron mass variation as well as to ponderomotive forces.Various analytical studies have focused on elucidating the properties of relativistic EM solitons, most of them focusing on a one-dimensional (1D) relativistic fluid model comprising of coupled nonlinear equations for the vector and scalar potentials, coupled to the plasma.Circularly polarized (CP) EM solitons in a cold plasma were investigated early on by Kozlov et al [2], who relied on a quasineutrality hypothesis to establish the occurrence of small-amplitude localized solutions (solitary waves).Later, computer simulations have shed some light on large-amplitude (relativistic) pulses, wherein charge separation in the plasma may be significant [2].Drifting envelope solitons of CP light were later modeled by Kaw and coworkers [3].A stationary (zero group velocity) relativistic EM soliton solution with was shown to exist within a 1D cold plasma model in Ref. 4. Further remarkable works on 1D solitary waves in cold plasmas based on a perturbation technique include the studies by Kuehl and Zhang [43] and by Fa-rina et al [42,44], who incorporated the effect of ion motion in a fluid-Maxwell description.Later, Poornakala et al [45] discussed the existence of bright or dark type envelope solitons and identified associated regions in parameter space.Thermal effects were taken into account in Refs.47 and 48.Single peak and multi-peak structures and their stability and mutual interaction properties were studied by Saxena et al, considering a nonuniform plasma [46], and also later in Refs.49 and 50.Note, for rigor, that the aforementioned studies have focused on non-magnetized plasma.
Magnetization during laser plasma interactions has attracted significant attention, as it impacts both the laser pulse dynamics and the substrate i.e. the background plasma properties.Inertial confinement fusion (ICF) has triggered interest in the study of magnetized plasma-laser interaction [51,52].Experimental observations have confirmed the generation of intense magnetic fields (up to hundreds of MG) [12,[53][54][55] during LPI.
As regards EM pulse propagation in magnetized plasma, Shukla and Stenflo [56] studied stationary solitary wave electromagnetic field structures adopting a slowly varying envelope approximation, followed by Nagesha et al [57], who adopted a slowly varying envelop approximation and neglected perturbations in the longitudinal component of the electron momentum, to find stationary solutions for CP EM waves in cold magnetized plasma.Weakly super-acoustic (supersonic) EM envelope waves were studied by Rao [58].The occurrence of standing 1D relativistic solitons in cold magnetized plasmas and the magnetic field's role on soliton stability was discussed by Farina et al [59], who found that the frequency interval for stable solitons depends on both the magnitude and the orientation of the magnetic field in a significant manner.Furthermore, the maximum field amplitude of the solitonic pulses essentially depends on the ambient magnetic field (strength).
In the article at hand, we have undertaken an analytical and numerical study of the interaction between left-and right-hand circularly polarized solitons in magnetized plasmas.Based on a fluid-Maxwell model, as starting point, we have established a system of scalar equations governing the propagation of circularly polarized EM waves in magnetized plasma.First, analysis of the linear regime has shown that, in certain frequency bands, left-hand circularly polarized (LCP-) and right-hand circularly polarized (RCP-) modes can propagate at the same group velocity.Next, we adopt a large-amplitude (nonlinear) approach, treating the model equations within an asymptotic multiscale expansion method.Assuming the existence of coupled localized (soliton) states for the magnetic vector potential, the analysis leads to a system of two coupled nonlinear Schrödinger (NLS) type equations for the respective EM field envelopes, i.e. for the right-hand circularly polarized (RCP-) and the left-hand circularly polarized (LCP-) mode respectively.The analysis of the coupled NLS equations shows that there are three types coupled (vec-tor) solitons for certain frequency bands, namely, (a) a RCP-propagating bright soliton coupled with a LCP-propagating bright soliton, (b) a RCP-propagating bright soliton coupled with a LCP-propagating dark soliton, and (c) a RCP-propagating dark soliton coupled with a LCPpropagating bright soliton.
This paper is organized as follows.In Section II, we introduce the fluid-Maxwell model for magnetized plasma and discuss the underlying physics.In Section III, the system is linearized and the dispersive modes obtained are identified.In Section IV, we assume a coupled-mode solution, and apply a multiscale perturbative technique to derive a system of (two) coupled NLS equations for the respective wavepacket amplitudes.In Section V, we present analytical results for each type of vector soliton.Finally, in the concluding Section VI, we summarize our results.

II. FLUID-MAXWELL MODEL FOR EM WAVES
We consider an electron-ion plasma, permeated by a uniform magnetic field along x.The instantaneous state of the plasma is given by a number of dynamical state variables, namely: the electron number density n, fluid speed u = u x and momentum (vector) p = (p y , p z ).These are subject to the influence of the electrostatic (scalar) potential φ and the magnetic (vector) potential A = (A y , A z ).All of these variables are functions of space x and time t, where one-dimensional (1D) propagation was considered, in the direction parallel to the magnetic field, i.e. along x.Given the (high) frequency of interest, the ions will be assumed to be stationary.
Our analysis will be based on the following closed system of scalar equations which govern the propagation of the circularly polarized EM waves in magnetized plasma: All physical quantities have been normalized by appropriate scales, namely: the scalar and vector potentials are normalized by mc 2 /e, the electric field E by mcω pe /e, the magnetic field by B by mω pe /e, the momentum by mc, the density by the n 0 , the electron velocity by the light velocity c.Thus, in this framework the length is measured in units of the skin length c/ω p0 , and time is measured in units of the plasma period (inverse plasma frequency) (ω −1 pe ), where ω pe = 4πn 0 e 2 /m e (here n 0 denotes the equilibrium electron density).Also, the electron momentum can be expressed as , for a circularly polarized EM pulse, where γ is the relativistic factor and we have used the notation α = +1 (α = −1) for left-(right-) hand circularly polarized electromagnetic waves.

III. MULTIPLE SCALE ANALYSIS FOR A SINGLE CIRCULARLY POLARIZED MONOCHROMATIC WAVEPACKET
Let the state vector S = (n, u, φ; A y , A z ; p y , p z ) describe the state of the system at a given position x and time t.Assuming a moderate (but finite) deviation from the equilibrium state S (0) = (1, 0, 0; 0, 0; 0, 0), one may adopt a multiple scales perturbation technique by expanding the system's state as a series in terms of powers of a small parameter ǫ ≪ 1 and subsequently considering multiharmonic analysis in each order.The method was described exhaustively in the past, e.g. in Ref. 60 in general, and also, for electrostatic and electromagnetic waves in plasmas, in Ref. 61 and 62 respectively, hence lengthy details will be omitted here.
The system's state S is expanded as where S (0) = (1, 0, 0; 0, 0; 0, 0) demonstrates the equilibrium state of the system and ǫ ≪ 1 is a dimensionless, real and small parameter.We now introduce new independent spatial and temporal variables, and accordingly expand the space and time derivative operators ∂ x = ∂ x0 + ǫ∂ x1 + . . .and ∂ t = ∂ t0 + ǫ∂ t1 + . ... Now, we use the assumption that the wavenumber k and the frequency ω (fast variables) affect the perturbed state only via the phase kx−ωt of the pulse, which are normalized by ω pe and k pe = ω/c, respectively.On the other hand, the amplitude of various frequency harmonics is assumed to be a slowly varying function in space and time, and will thus assumed to vary only on the slower scales (x 2 , x 3 , ...; t 2 , t 3 , ...).The solutions to be sought of Eqs.(1-7) are therefore expressed in the form: The methodology adopted in this work relies on a Newell type multiscale technique, well known from nonlinear optics (where, interestingly space and time are reversed, hence dispersion occurs in the frequency domain, rather than the wavenumber space).The fundamental aim of the method it to allow for a slow variation of the wavepacket's amplitude (envelope) in space and time.An excitation is considered around the equilibrium state (represented by the zeroth order state in a smallness parameter (ǫ 0 ), while the linear (harmonic) carrier wave regime is obtained in first order (∼ ǫ 1 ) and the group velocity frame is defined in the second order (∼ ǫ 2 ), upon imposing certain compatibility conditions (namely, dictating the annihilation of secular terms occurring in that order).In each order in ǫ, a secondary expansion is considered in harmonics, from the zeroth to the n-th harmonic.Within this (tedious but straightforward) scheme, the first harmonic (carrier) amplitudes arise as solutions in the first order, while the second harmonic amplitudes will arise via the second-order second-harmonic amplitudes.Zeroth-order amplitudes will originate from the 2nd order 2nd harmonic expressions, in combination with their 3rd order counterparts.Finally, an evolution equation is obtained for the envelope, by eliminating the secular terms occurring in 3rd order (1st harmonics).Linear analysis is therefore embedded in this analysis, thus generalized to account for the generation of higher harmonics.Substituting Eqs. ( 9)-( 10) into Eqs.( 1)-(-7)) and collecting the terms arising in each order in ǫ, we obtain the amplitude evolution equations at successive orders.
In order O(ǫ 2 ), we have the equations Same as in the Ref. [62], in the first harmonic ((ℓ = 1), the condition for annihilation of secular terms leads to y/z or A y/z ), implying that the envelope moves at the group velocity v g = ω ′ (k), given by: Following the same procedure in order O(ǫ 3 ), and using the variables X = x 1 −v g t 1 ≡ ǫ(x−v g t) and T = t 2 ≡ ǫ 2 t, we derive from the non-secularity condition at O(ǫ 3 ) the NLS Eq. ( 37) where ψ denotes the amplitude A (11) y , the (slow) time and space variables are τ = t 2 and ξ = x 1 − v g t 1 , and the dispersion coefficient P and the nonlinear self-phase modulation (SPM) coefficient Q are respectively given by the (real) expressions:

IV. MULTIPLE SCALE ANALYSIS FOR COUPLED MODES
Now, we consider for a solution in the form where subscripts j = 1 and j = 2 correspond to the LCP and RCP waves, S 11 j is an unknown complex function, θ j = k j x 0 − ω j t 0 , while the wavenumbers k j and frequencies ω j satisfy the dispersion relation provided in Eq. ( 26).In Fig. 1, we show the dispersion relation for the Ω = 0.2, Ω = 0.8, and Ω = 2.It is clear, from Fig. 1, that for LCP EM wave (top panel) there exist two frequency bands where the propagation is possible: the high-frequency band and the low-frequency band.For Ω = 0.2 (depicted by solid (blue line) the propagation is possible for ω > 1.105 (high-frequency band) and for ω < 0.2 (low-frequency band) namely a gap appears between 0.2 < ω < 1.105 where EM wave propagation is not possible.For Ω = 0.8 (depicted by dashed (red) line) the propagation is possible for ω > 1.477 (highfrequency band) and for ω < 0.8 (low-frequency band) where the gap appears between 0.8 < ω < 1.477.Also, for Ω = 2 (depicted by thick (green) line) the propagation is possible for ω > 2.414 (high-frequency band) and for ω < 2 (low-frequency band).In this case the gap appears between 2 < ω < 2.414.Moreover, for the RCP EM wave (bottom panel) propagation is possible for frequencies ω > 0.9047 when Ω = 0.2 (depicted by solid (blue line).In case where Ω = 0.8(depicted by dashed (red) line) and Ω = 2 (depicted by thick (green) line) the propagation is possible for frequencies ω > 0.6767 and ω > 0.414,respectively.Furthermore, the gap which appears for LCP EM waves increases (decreases) as the Ω decreases (increases).It is important to notice that there is a frequency region, which is independent from polarization of EM waves, where the propagation cannot take place for values of Ω < 0.707.
Below we will demonstrate that in the nonlinear setting, coupling between LCP and RCP EM waves with equal group velocities is possible, and we consider the interaction between LCP and RCP waves.Observing Fig. 1 and knowing that the group velocity at a given wavenumber represents a tangent (line) to the dispersion curve, it is possible to identify domains, with equal group velocities.Thus, using Eq. ( 26) we obtain the group velocity v g = ∂ω/∂k: i.e.Eq. ( 36) above.This result is presented in Fig. 2, where the group velocity v g is plotted as a function of the normalized frequency ω.Notice that the dashed (red) solid (blue) and solid (blue) lines show the group velocity υ g for RCPand LCP-EM waves, respectively.
As seen in Fig. 2 using a horizontal cut we can to observe that there are a RCP and a LCP electromagnetic wave which have the same group velocity (and are coupled in the nonlinear regime).
The point A,B, C and D, in Fig. 2, have the same group velocity, namely, v g = 0.4.Also, we see that there is a maximum possible common v ℓ gmax in the low frequency band for LCP wave.Thus, one can divide each of the group-velocity curves into three sub-regions, depending on the sign of the group-velocity dispersion (GVD), ∂v g /∂ω, where the interaction with equal group velocities may appear.As depicted in Fig. 2 these subregions are: (a) the sub-bands I and II for the low frequency band of the LCP EM wave, characterized by positive and negative GVD respectively, and the sub-bands III for the high frequency band of the LCP EM wave, characterized by positive GVD, (b) the band of the RCP EM wave, again characterized by positive GVD respectively.Thus, nonlinear LCP and RCP modes of equal v g can feature the following three different possible interactions: • Case 1:band of the RCP-mode and the LCP-mode in band I, both featuring positive GVD for v g ≤ v ℓ gmax .
• Case 2: band of the RCP-mode and the LCP-mode in band II; here, the RCP (LCP) mode features positive (negative) GVD for v g ≤ v ℓ gmax .
• Case 3: band of the RCP-mode and the LCP-mode in band III, both featuring positive GVD for 0 < v g ≤ 1.
As shown in Fig. 3, upon inspection of the groupvelocity curves for Ω = 0.8 (middle panel) we see that Notice that, v g = 1 is the group velocity for unmagnetized plasma.It is obvious that the above considerations is the result of the existence of the frequency gap as shows in Fig. 3 for all values of Ω, i.e, for Ω = 0.2, Ω = 0.8and Ω = 2. Furthermore, we can observe that as Ω increases the v ℓ g is raised and the gap is smaller.When the Ω ≫ 1 then v g → 1 and the LCP and RCP wave are degenerated.

A. Nonlinear analysis: the coupled NLS equations
In our case, since we study the interaction between a LCP and a RCP nonlinear mode with equal group velocities, we seek for a solution of Eqs.(1-7) in the form: where "c.c." denotes complex conjugate.In Eq. ( 41), subscripts j = 1, 2 correspond to the LCP and RCP mode, S j (X, T ) are unknown (continuous) slowlyvarying envelope functions depending on the slow scales X = ǫ(x − v g t) (where v g is the common group velocity) and T = ǫ 2 t, while exp(iθ j ), with θ j = k j x − ω j t, are the carriers of frequencies ω j and wavenumbers k j .The ǫ is a formal small parameter related to the soliton amplitude (see below).According to Eq. ( 41), the field S (n) j is the leading-order form of a more general ansatz employing multiple time and space scales.In this context, use of a formal multi-scale expansion method leads to a hierarchy of equations at various powers of ǫ, which are solved up to the third-order.Indeed, at orders O(ǫ) (linear limit) and O(ǫ 2 ), we derive the dispersion relation, Eq. ( 26), and the group velocity, Eq. ( 36), respectively.At the order, O(ǫ 3 ) we obtain the nonlinear coupled NLS equations as where the amplitude Ψ j , the normalized GVD coefficients D j , the self-phase modulation (SPM) coefficients g jj , and the cross-phase modulation (CPM) coefficients g j,3−j (with j = 1, 2) are respectively given by: Note that, using the Eqs.( 26) and ( 36) the above relations be cast in the form: As shown from Eqs. ( 46)-( 47) the coefficients g jj , and g j,3−j are always positive.Next, using scale transformations, we measure normalized time T and densities |Ψ j | 2 in units of |D 1 | −1 and |D 1 /g jj | respectively, and cast Eqs. ( 42)-( 43) in the form: 51) where where the coefficients λ 1,2 and d are positive.As arises from Eqs. ( 50) when the coupling coefficients (λ j = 0) are zero then the evolution of either the LCP wave Ψ 1 or the RCP wave Ψ 2 is described by a single NLS equation; the latter, supports soliton solutions of the dark or the bright type, depending on the relative signs of dispersion and nonlinearity coefficients (see, e.g., Ref. [70]).Particularly, the mode Ψ 1 supports dark solitons for s < 0 or bright solitons for s > 0. The mode Ψ 2 supports only bright solitons because the coefficient d is always positive.However, because the above conditions are modified when λ j = 0 various types of vector (coupled) solitons can be found as seen from the Eqs.( 50) and ( 51).
According to the Ref. [71], four types of vector solitons are possible: • bright-bright (BB) solitons, in the form: • bright-dark (BD) solitons, in the form: • dark-bright (DB) solitons, in the form: • dark-dark (DD) solitons, in the form: In the above equations, η j (j=1,2) and Ψ j,0 declare the frequencies and amplitudes of each soliton, while b is the (common) inverse width of the solitons.Now, each of the above ansatz is substituted in Eqs. ( 50) and ( 51), leading to a set of equations connecting the soliton parameters.Particularly, the equations connecting the amplitudes Ψ j,0 and the inverse width b are of the form: for the BB, BD, DB and DD solitons respectively, where parameters α j (j = 1, 2) are given by: It is worth noting that a similar set of coupled NLS equation has been derived in the past for modulational wavepacket interaction in electron-ion plasmas, considering either upper-hybrid [72] or, more recently, electrostatic [73][74][75] waves.Although the mathematical setting is analogous, and may therefore interest the reader of this paper, the physical background in this paper is distinct from those works; in those earlier studies, for instance, an unmagnetized plasma was considered, in the electrostatic approximation (i.e.neglecting magnetization, thus precluding electromagnetic excitations), with the focus being on the nonthermal (kappa-distributed) background electron statistics (a common occurrence in Space plasmas).Contrary to that picture, our focus here is on electromagnetic solitary waves propagating in magnetized plasma and, in particular, on the role of the ambient magnetic field on their structural and propagation characteristics.

V. SOLITON INTERACTIONS IN DIFFERENT FREQUENCY BANDS
It is important to understand that the sign of the parameters α j (j = 1, 2) defined above determines the type of vector solitons that may occur in the plasma, as dictated by Eqs. ( 61)- (64).Also, the parameters α j depend on the frequency and the magnetic field according to Eq. ( 65), through coefficients λ 1,2 , s and d.Consequently, we expect that the precise solitons solutions will arise from the investigation of the sign of the parameters α j and the value of magnetic field, for all cases (Case 1, 2 and 3) as have been defined in Section V. On the other hand, the Eqs.( 50)-( 51) are no longer of the Manakov type and, thus, generally, they are not completely integrable.Nevertheless, standing wave soliton solutions whose exact analytical form depend on by sign of the parameters α j can still be found.That is valid for all cases of interactions (Case 1, 2, and 3).Initially, we study the coupling between a propagating soliton, with a frequency that is in band I, and a propagating soliton, with a frequency that is in RCP band.In Fig. 4 we depict the dependence of the parameters λ 1 , λ 2 and d, as a function of the normalized frequency ω (for Ω = 0.2, Ω = 0.8 and Ω = 2).Now, we have s = +1 (cf.Fig. 3) and λ 1,2 > 0 and d > 0, as can be seen in the Fig. 4. Also, the parameters α j (j = 1, 2) are given by: In this case, the Eqs.( 50)-( 51) are no longer of the Manakov type and, thus, generally, they are not completely integrable.Nevertheless, standing wave solitons solutions can still be found whose exact analytical form depend on by sign of the parameters α j .1.When α 1 > 0 and α 2 > 0 then we have a case of coupled solitons, a dark soliton and a bright soliton, whose exact analytical form are given by Eqs. ( 57)-( 58), where the soliton amplitude parameters η 1,2 and the inverse width b are connected via the following equations: It is thus clear that the above solutions have one free parameter.Employing the solutions ( 57)-( 58), we can approximate the vector potential A y (x, t), in terms of the original coordinates x and t, as follows: The amplitude A 0 and the frequencies ω 0j (j = 1, 2) of the soliton are given by: In this case, In Fig. 5 we show the profiles (at t=0) of the dark (LCP mode) and bright (RCP mode) solitons in the absence of coupling (top panel) as well as the dark-bright soliton (bottom panel), for Ω = 0.8.Now, the group velocity of the DB soliton (common for both components) is v g = 0.2965, which occurs when the angular frequencies for the modes Ψ 1 and Ψ 2 take, respectively, the values ω 1 = 0.0335 (k 1 = 0.2119) and ω 2 = 0.7007 (k 2 = 0.1551).2. When α 1 < 0 and α 2 < 0 then we have a case of coupled bright solitons, whose exact analytical form are given by Eqs. ( 53)- (54).Now, the parameters η 1,2 and b are connected via the following equations: The vector potential A y (x, t), in terms of the original coordinates x and t and soliton parameters, are given by the Eq. ( 71) and Eqs. ( 72)-( 73) respectively, where, now In Fig. 6 we show the profiles (at t=0) of the bright (LCP mode) and bright (RCP mode) solitons in the absence of coupling (top panel) as well as the bright-bright soliton (bottom panel), for Ω = 0.8.Now, the group velocity of BB soliton (common for both components) is v g = 0.4572, which occurs when the angular frequencies for the modes Ψ 1 and Ψ 2 take, respectively, the values ω 1 = 0.1454 (k 1 = 0.4932) and ω 2 = 0.77414 (k 2 = 0.2621).
3. When α 1 > 0 and α 2 < 0 then we have a case of coupled solitons, a bright soliton and a dark soliton, whose exact analytical form are given by Eqs. ( 55)- (56), where the soliton amplitude parameters η 1,2 and the inverse width b are connected via the following equations: Now, we can approximate the vector potential A y (x, t), in terms of the original coordinates x and t, as in the Eq.(71) where, in this case, In this case, the solution amplitude A 0 and the frequencies ω 0j (j = 1, 2) are given by Eqs. ( 72)- (73).
In Fig. 7 we show the profiles (at t=0) of the bright (LCP mode) and dark (RCP mode) solitons in the absence of coupling (top panel) as well as the bright-dark soliton (bottom panel), for Ω = 0.8.Now, the group velocity of the BD soliton (common for both components) is v g = 0.4647, which occurs when the angular frequencies for the modes Ψ 1 and Ψ 2 take, respectively, the values ω 1 = 0.1743 (k 1 = 0.5558) and ω 2 = 0.744 (k 2 = 0.2677).
The above table shows the existence of vector solitons according to the sign of the coefficients α 1 and α 2 , for the case 1.Moreover, as are depicted in Fig. 4, for Ω = 0.8, the above solutions are defined only for a narrow frequency band (0−0.129),(0.136−0.166) and (0.17−0.198) for dark-bright (DB), bright-bright (BB) and bright-dark (BD) interactions , respectively.Also, as is observed in the Fig. 4, as Ω is decreased (increased) lead to an increase (decrease) of the width of the DB band where DB solitons can be formed, while the width of BB and BD band is decreased (increased).As seen in this figure, for Ω = 0.2 the DB, BB and BD band possess, approximately, the 80%, 10% and 6% of I-band, respectively.Notice that for Ω = 0.8 (Ω = 2) the DB, BB and BD band possess, the 65% (47%), 15% (21%) and 14% (30%) of I-band, respectively.
B. Case 2: Bright-dark solitons in bands II and RCP.
In this case s = −1, we study the coupling between a propagating soliton, with a frequency that is in band II, and a propagating soliton, with a frequency that is in RCP band.In Fig. 8 we depict the dependence of the parameters λ 1 , λ 2 and d, as a function of the normalized frequency ω (for Ω = 0.2, Ω = 0.8 and Ω = 2).Now, the coefficients α 1 and α 2 take the following form: Following our previous considerations, here we have a case of coupled bright-dark solitons where the exact analytical form is given by Eqs. ( 55)- (56).The soliton amplitude parameters η 1,2 and the inverse width b are connected via the following equations: The analytical form for vector potential A y (x, t) and the solitons parameters are same with the previous case, namely, in the interaction in bands I and RCP [c.f Eqs.(84)-( 85)].
C. Case 3: Solitons in bands III and RCP.
Finally, we study the coupling between a propagating soliton, with a frequency that is in band III, and a propagating soliton, with a frequency that is in RCP band.Now, s = +1 (cf.Fig. 3) while the other dispersion and nonlinearity coefficients are shown in Fig. 9 as functions of the normalized frequency ω for Ω = 0.2, Ω = 0.8 and Ω = 2. Also, as is observed in the Fig. 9, as Ω is increased lead to an increase of the width of the BD band where BD solitons can be formed, while the width of BB is decreased.

D. Discussions
The analysis of previous section end in some important results which confirm our expectations.First of all in the Case 1 (Solitons in bands I and RCP) we obtain three types of vector solitons, namely bright-bright, brightdark and dark-bright.Also, we observe that for low frequencies and Ω values the width of the dark-bright band where dark-bright solitons can be formed possesses the larger percentage of the band.However as Ω is increased the occurrence of the dark-bright type is decreased while the width of the bright-bright and mainly also that of the bright-dark category are increased.In Case 2 (Case 2: Bright-dark solitons in bands II and RCP) for all the frequencies and Ω values, we have only the existence of bright-dark solitons.Finally, in Case 3 (Solitons in bands III and RCP) we have only bright-bright solitons in all the frequencies for low Ω values.As Ω is increased bright-dark solitons appeared which dominate in all frequency band for the large values of Ω.
According to the table in Fig. 10, we found that brightdark solitons exist in all frequency bands for all values of Ω.Furthermore, we observe bright-bright solitons in bands (I-RCP) and (III-RCP) for the small values of Ω.Also, we found that the dark-bright solitons exist only in I-RCP frequency band for small values of ω and Ω.
Generally our detailed analysis can support various coupled NLS systems in magnetized plasmas with the equal group velocities between the propagating waves.Our analytical investigation reveals the existence of various type vector solitons taking into account the certain frequency bands and the strength of the magnetic field.It is very important to note that from our suggested method the capability control of a solitonic mode e.g.LCP (RCP) from the change the amplitude or width of other mode RCP (LCP) emerges.Although our present work has an analogous mathematical settings with the works of Refs.[62], [71], [73], and [75] it has serious differences as: a) The physical background in this paper is different compared to Refs.[71], [73], and [75] .In Refs.[73] and [75] the study was carried out in the framework of an unmagnetized plasma while in the Ref. [71] the physical background was completely different where a nonlinear composite right left handed transmission line was studied, b) With regards to Ref. [62], where the physical background is same, namely the magnetized plasma, in this paper we study the interaction between left hand circularly polarized (LCP) and right hand circularly polarized (RCP) solitons while in Ref. [62], we investigated only the beam-plasma interaction, c) As previously mentioned, in this paper, we investigate the interaction between of a LCP and RCP electromagnetic wave which have the same group velocity in different frequency bands.In Refs [73] and [75] the electrostatic wave packets are studied in a regime with different group velocities, d) In this work, we study the vector soliton solutions in different frequency bands and for different magnetic fields.The analysis in Ref. [71] was carried out only for different frequency bands, in Ref. [73], the instability growth rate was investigated while in the Ref. [75], the analysis was focused on the values of the κ-index of k-distribution function.

VI. CONCLUSION
In conclusion, we have performed analytical techniques to study the existence of coupled left-and right-circularly polarized solitons in magnetized plasmas.Our analysis started with the derivation of a closed system of scalar equations which the govern of the propagation of LCP and RCP EM waves in magnetized plasmas.Upon the linear limit, the dispersion relation has been derived and we have shown that there are LCP-and RCP modes, with the same group velocity, which can be propagate for specific frequency bands.Next, in the nonlinear regime we used the multiple scale perturbation method in order to investigated the coupling between LCP-and RCP modes.In this framework, we obtain, a system of two coupled nonlinear Schrödinger (NLS) equations for the unknown vector potential envelope functions.The above system of equations was used in order to predict the existence of coupled left-and right-circularly polarized solitons in magnetized plasmas, of the bright-bright, brightdark and dark-bright type.Thus, we found vector solitons of the same type (namely, bright-bright) as well as ones of the mixed type (namely, bright-dark and dark-bright) in specific frequency bands.Also, the values of Ω change the width of frequency band where the vector solitons can be formed.Thus, we found vector solitons of the type dark-bright, bright-bright and bright-dark in the bands I and RCP where the increase of the Ω lead to a decrease of the width of the band of dark-bright solitons with simulta-neous increase of the width of the band of the brightbright and bright-dark solitons.In the bands II and RCP only bright-dark solitons were found.Finally, we found bright-dark and bright-bright solitons in the bands III and RCP where the increase of the Ω lead to a decrease of the width of the band of bright-bright solitons with simultaneous increase of the width of the band of the bright-dark solitons.
Our results should be of interest to the plasma science community, where beam-plasma interactions are a hot topic of research nowadays, both theoretically and in the laboratory (laser plasma interaction experiments), and also in nonlinear optics, where coupled beam propagation is an ubiquitous point of focus.Beyond these physical settings, coupled NLS equations such as our Eqs.( 24)-( 25) arise in various contexts, ranging from hydrodynamics to supraconductivity (Bose-Einstein condensates) and left-hand materials, to mention a few.The mathematical setting and the methodology of our study is thus expected to apply in a wider context.

FIG. 2 .
FIG.2.The group velocity υg as a function of the normalized frequency ω.The dashed (red) and solid (blue) dashed (red) lines show the group velocity υg for RCP-and LCP-EM waves, respectively.

A. Case 1 :
Solitons in bands I and RCP.

FIG. 10 .
FIG. 10.Different types of solitons as a function of the angular (normalized) frequency ω are depicted in the table, for the different values of Ω.