Home Ion-Acoustic Cnoidal Waves with the Density Effect of Spin-up and Spin-down Degenerate Electrons in a Dense Astrophysical Plasma
Article Publicly Available

Ion-Acoustic Cnoidal Waves with the Density Effect of Spin-up and Spin-down Degenerate Electrons in a Dense Astrophysical Plasma

  • Nimardeep Kaur , Rupinder Kaur and N. S. Saini EMAIL logo
Published/Copyright: December 7, 2019

Abstract

An investigation of nonlinear ion acoustic (IA) cnoidal waves in a magnetised quantum plasma is presented by using spin evolution quantum hydrodynamics model, in which inertial classical ions and degenerate inertialess electrons with both spin-up and spin-down states taken as separate species are considered. The Korteweg–de Vries equation is derived using the reductive perturbation method. Further, using the Sagdeev pseudopotential approach, the solution for IA cnoidal waves is derived with suitable boundary conditions. There is the formation of only positive potential cnoidal, and in the limiting case, positive solitary waves are observed. The effects of density polarisation and other plasma parameters on the characteristic features of cnoidal and solitary waves have been analysed numerically. It is seen that the spin density polarisation significantly affects the characteristics of cnoidal structures as we move from strongly spin-polarised (μ = 1) to a zero spin-polarisation case (μ = 0). The results obtained in the present investigation may be useful in comprehending various nonlinear excitations in dense astrophysical regions, such as white dwarfs, neutron stars, and so on.

1 Introduction

In recent years, dense quantum plasmas have emerged as an active field of research due to their great relevance in different areas of practical importance, e.g. nanoscale electromechanical systems [1], [2], laser interactions with atomic systems [3], and in dense astrophysical systems [4], such as neutron stars, white dwarfs, and so on. Due to the high number density and low particle temperature of particles, quantum plasmas are distinguished significantly from the classical plasmas, where the density of particles is relatively low and possesses high plasma temperature. To investigate the various astrophysical phenomena in interstellar compact objects, dense quantum plasmas would be helpful in establishing a suitable frame. The density of the interiors of the interstellar objects is significantly high such that the nonthermal pressure is provided by the degenerate fermion/electron pressure, as well as interaction of particles.

Mathematically, Chandrashekhar [5], [6], [7] deduced the equation of state in such compact interstellar objects for the degenerate electrons with Pene53 for the nonrelativistic case and Pene43 for the ultrarelativistic case, where Pe and ne are the pressure and number density of degenerate electrons. In highly compressed plasma species, the uncertainty in momenta is infinitely large, which implies that the degenerate plasma species must move very fast (despite that they are extremely cold), giving rise to a very high pressure, called as “degenerate pressure.” The degenerate pressure depends only on the number density of electrons and not on their temperature. In the quantum hydrodynamic (QHD) model, which is considered as the quantum counterpart of the classical fluid model [8], the inclusion of the Fermi pressure and the Bohm potential term modifies the momentum equation of the charged particles [9]. The quantum ion acoustic (IA) waves and the role of quantum diffraction effects have been studied by Haas et al. [10] using the QHD model. It has also been found that the system supports the travelling waves with periodic patterns in the fully nonlinear regime. Haas [11] devised an ideal quantum magnetohydrodynamic (QMHD) model with the incorporation of the quantum diffraction effects with relevance to the dense astrophysical objects such as interiors of white dwarfs. Using the QHD and QMHD models, various nonlinear electrostatic and electromagnetic waves have been studied in quantum plasmas.

The spin effects are considered as one of the most important properties of quantum plasmas due to great significance of highly magnetised quantum plasmas in the atmospheres of neutron stars [12], [13]. Marklund and Brodin [14], [15] extended the QMHD model and proposed the spin-12 QMHD model for hydrodynamic waves. They found that the spin effects can significantly alter the characteristics of low-frequency electromagnetic modes. The dynamics of fast and slow magnetosonic waves has been studied by Mushtaq and Vladimirov [16] in two-dimensional spin-12 quantum plasmas using QMHD model. The linear dispersion relation for both slow and fast quantum magnetosonic waves has been discussed in detail. Andreev [17] considered the propagation of waves with degenerate electrons in the magnetised plasmas and derived the QHD equations for spin-12 particles. The author examined the evolution of electrons and contribution of magnetic field in the Langmuir wave dispersion via the different occupation of spin-up and spin-down states. The existence of spin-electron-positron acoustic spin-electron acoustic, and positron acoustic solitons in degenerate electron-positron-ion (epi) plasmas were demonstrated by Iqbal and Andreev [18]. The influence of the spin polarisation on the properties of three kinds of solitons was also analysed. It has been reported that there exist different types of spin-up and spin-down electrons in metal in the absence of an external magnetic field, which yields magnetic moment of metal as a whole equal to zero [19]. But in the presence of an external uniform magnetic field, there are more electrons with spin-up aligned along the direction of magnetic field than that of spin-down electrons. Owing to the presence of spin-up and spin-down electrons in a degenerate quantum plasma, some investigations focusing on the study of solitary and shock waves have been reported. Ahmad et al. [20] investigated the low-frequency electrostatic waves in plasmas having inertialess degenerate electrons by employing the separated spin evolution QHD model. The pulse stability analysis was carried out, and it was also found that spin polarisation significantly affects the amplitude and width of the solitary waves. Hussain and Mahmood [21] investigated the propagation characteristics of IA shock waves in a dense magnetised plasma with relative density effects of spin-up and spin-down degenerate electrons. The Korteweg–de Vries Burgers equation was solved numerically to analyse the influence of the spin density polarisation ratio on the propagation characteristics of shock waves.

The study of cnoidal waves has become one of the important areas of research because of their wide range of applications in nonlinear transport processes in plasmas [22], ionosphere plasmas [23], single-mode drift wave spectra [24], and so on. A variety of investigations have been reported by numerous researchers [25], [26], [27], [28] to study the characteristics of cnoidal waves in different plasma regimes. Kaladze et al. [29] deduced Korteweg–de Vries (KdV) equation by employing the reductive perturbation method and investigated electrostatic acoustic nonlinear periodic waves in unmagnetised pair-ion plasmas that constitute the same mass ion species with different temperatures. El-Shamy [25] investigated the propagation characteristics of IA cnoidal waves in a dense relativistic degenerate magnetoplasma consisting of relativistic degenerate electrons and nondegenerate cold ions. The various solutions of nonlinear cnoidal and solitary waves were presented numerically. Ur-Rehman et al. [26] derived the KdV equation in a magnetised ei plasma with cold ions and warm electrons. The impact of various plasma parameters on the characteristics of compressive magnetoacoustic cnoidal waves was studied. The propagation properties of dust acoustic cnoidal waves in an unmagnetised ion beam dusty plasma were investigated by Kaur et al. [27]. Using reductive perturbation technique, the KdV equation was derived, and the solution of nonlinear cnoidal waves was determined by applying the appropriate boundary conditions. The characteristic features of magnetosonic cnoidal and solitary waves were investigated in a magnetised electron-ion-dust (eid) plasma by Kaur et al. [28]. It was found that there exist only positive potential magnetosonic cnoidal and solitary waves in the limit of high β plasma.

Earlier investigations were focussed on the study of solitary and shock waves in a degenerate quantum plasma with the effects of the spin-up and spin-down of the electrons. To the best of our knowledge, the study of IA cnoidal waves with the effects of both spin-up and spin-down of degenerate electrons in a magnetised quantum plasma has not been reported yet. Our aim in the present investigation is to numerically analyse the effect of spin polarisation density and other physical parameters on the propagation characteristics of IA cnoidal waves in dense astrophysical plasma with implication to the region of white dwarfs. The layout of the manuscript is as follows: In Section 2, the basic fluid equations are introduced. In Section 3, the derivation of the nonlinear KdV equation by using reductive perturbation method is presented. The cnoidal wave solution of KdV equation is given in Section 4. In Section 5, the numerical analysis for IA cnoidal waves and in the limiting case solitary waves has been presented. Conclusions are highlighted in Section 6.

2 Basic Fluid Equations

We consider a magnetised quantum plasma in which a uniform magnetic field is applied along the positive z-direction, B=B0z^. The separated spin evolution QHD model has been utilised by considering both the spin-up (ne) and spin-down (ne) states for degenerate electrons and inertial nondegenerate ions. Here, we assume that wave propagates in three dimensions, i.e. =x,y,z. Nondegenerate ions are taken due to their heavy mass as compared to the mass of electrons. The separate evolution of spin-up and spin-down electrons helps to discover the new longitudinal waves, which depend on parallel and perpendicular propagation to an external magnetic field. Their existence is related to the different population of spin-up and spin-down electrons in equilibrium plasmas [17]. The continuity and momentum equations for ions are given as

(1)nit+(niui)=0,
(2)(t+ui)ui=emiE+eB0micui×z^.

Electrons are inertialess, and the momentum equations of electrons with spin-up (ne) and spin-down (ne) states are given as

(3)0=ene()(E+B0cue())Pe()+22mene()(2ne()ne()).

The Poisson’s equation is given as

(4)E=4πe(nine+ne).

The Fermi pressure of degenerate electrons with spin-up (ne) and spin-down (ne) is Pe=KBTFene535n023 and Pe=KBTFene535n023. TFe=6π2n02322KBme and TFe=6π2n02322KBme are the Fermi temperatures for spin-up (ne) and spin-down (ne) electrons, respectively; KB is the Boltzmann constant. In this model, ion temperature is ignored as compared to electron temperature.

In component form, the normalised continuity and momentum equations for ions are written as

(5)nit+x(niuix)+y(niuiy)+z(niuiz)=0,
(6)tuix+Duix=ϕx+Ωuiy,
(7)tuiy+Duiy=ϕyΩuix,
(8)tuiz+Duiz=ϕz.

The component forms of the normalised momentum equations for spin-up and spin-down electrons are written as

(9)0=ϕxΩuey()(2μ())233ne()13ne()x+H22x(2ne()ne()),
(10)0=ϕyΩuex()(2μ())233ne()13ne()y+H22y(2ne()ne()),
(11)0=ϕz(2μ())233ne()13ne()z+H22z(2ne()ne()).

The normalised Poisson’s equation is

(12)2ϕx2+2ϕy2+2ϕz2=μne+μneni

where μ=n0ni0 and μ=n0ni0. The charge neutrality condition yields ne0+ne0=ni0. The density polarisation ratio is given as

(13)μ=ne0ne0ne0+ne0,

It is seen that μ can be positive, zero, or negative. 0 < μ (μ < 0) corresponds to plasmas where the number density of spin-up electrons is higher (lower) than that of the spin-down electrons. When two electron species have equal number densities, then μ = 0. From (13), the negative value of μ has no physical significance because in the presence of an external uniform magnetic field there are more electrons with spin-up aligned along the direction of the magnetic field than that of spin-down electrons [19]. For the case μ = 1 (i.e. electrons have only one state of spin), one can simply recover the case of ei plasmas. Therefore, μ=1+μ2 and μ=1μ2. D=u is the convective derivative, Ω=ωciωpi, where Ω is the ion-cyclotron frequency, ωci=eB0mic is normalised by ion plasma frequency ωpi=4πni0e2mi, and H=ωpimimecs2 is the quantum parameter. The fluid velocities (uix,uiy,uiz) are normalised by the quantum IA speed cs=2kBTFmi; the number densities ni,ne,ne are normalised by their respective equilibrium number densities and electrostatic potential ϕ by 2kBTFe. The space coordinates are normalised by λi=csωpi and time by the inverse of the ion-plasma frequency (ωpi1).

3 Derivation of the KdV Equation

We have employed the reductive perturbation method to derive the KdV equation for the study of cnoidal waves in a magnetised dense astrophysical quantum plasma. The stretching coordinates are given as

(14)ξ=ε12(kxx+kyy+kzzvpht),τ=ε32t.

where kx,ky, and kz are the direction cosines such that kx2+ky2+kz2=1, and vph is the phase velocity of the waves. Here, ε is a small expansion parameter. The perturbed quantities ni,uix,uiy,uiz, and ϕ are expanded in terms of smallness parameter ε in the following form:

(15)ni,e=1+εni,e(1)+ε2ni,e(2)+u(i,e)x=ε32u(i,e)x(1)+ε2u(i,e)x(2)+u(i,e)y=ε32u(i,e)y(1)+ε2u(i,e)y(2)+u(i,e)z=εu(i,e)z(1)+ε2u(i,e)z(2)+ϕ=εϕ(1)+ε2ϕ(2)+

Using (14) and (15) in (5)–(12) and collecting the terms in the lowest order of ε, we get the following first-order evolution equations.

From ion continuity equation:

(16)vphni(1)ξ+kzuiz(1)ξ=0,

From the z-component of the ion momentum equation:

(17)kzϕ(1)ξ=vphuiz(1)ξ,

From the z-component of electron momentum equation having the spin-up state:

(18)ϕ(1)ξ=(2μ)233ne13ne(1)ξ,

From the z-component of the electron momentum equation having the spin-down state:

(19)ϕ(1)ξ=(2μ)233ne13ne(1)ξ.

From Poisson’s equation:

(20)μne(1)+μne(1)=ni(1),

On solving (16)–(19), the phase velocity for IA cnoidal waves having both spin-up and spin-down states of degenerate electrons is obtained as

(21)vph=kz2233(μ13+μ13)

On integrating (16)–(19) and applying some algebraic manipulations, we get the following relations:

(22)ni(1)=kz2vph2ϕ(1)+C1(τ),
(23)uiz(1)=kzvphϕ(1)+C2(τ),
(24)ne(1)=3(2μ)23ϕ(1)+C3(τ),
(25)ne(1)=3(2μ)23ϕ(1)+C4(τ).

where C1(τ), C2(τ), C3(τ), and C4(τ) are the integration constants, which may depend on τ but independent of ξ.

Collecting the next higher-order terms of ε in the ion continuity equation, we have

(26)vphni(2)ξ+ni(1)τ+kxuix(2)ξ+kyuiy(2)ξ+kzuiz(2)ξ+kzξ(ni(1)uiz(1))=0.

Collecting the next higher-order terms of the z-component of the momentum equation of ions, we have

(27)vphuiz(2)ξ+uiz(1)τ+kzuiz(1)uiz(1)ξ=kzϕ(2)ξ.

The next higher-order terms of the z-component of the momentum equations of electrons having spin-up and spin-down states give the following equations:

(28)kzϕ(2)ξ(2μ())233kzne()(2)ξ+H24kz33ne()(1)ξ3+(2μ())239kzne()(1)ne()(1)ξ=0.

The next higher-order terms of the Poisson’s equation

(29)2ϕ(1)ξ2=μni(2)+μni(2)ni(2).

On simplifying (26)–(29) by eliminating the second-order quantities and making use of the boundary conditions C1(τ)τ=C2(τ)τ=C3(τ)τ=C4(τ)τ=0, we obtain the following form of the KdV equation:

(30)ϕ(1)τ+Aϕ(1)ϕ(1)ξ+B3ϕ(1)ξ3+Cϕ(1)ξ=0

where the nonlinear coefficient A, dispersion coefficient B, and the coefficient C are given as

(31)A=3kz2(11243(μ13+μ13)),B=vph32kz2(1+(1kz2)Ω29H2243(μ13+μ13)),C=vph32kz2(kz4vph4C2μ13223C3μ13223C4).

In the limiting case, when coefficient C vanishes in (38) of [20] and in (29) of [21], one can get the KdV equation, which agrees with (30) (for C = 0) in the present investigation.

4 Cnoidal Wave Solution of the KdV Equation

The stationary solution of (30) is obtained by using the transformation η=ξU1τ; here, U1 is velocity of nonlinear cnoidal waves. By changing (30) into η coordinate and applying some algebraic simplifications, the energy balance equation is obtained as

(32)12(ϕη)2+V(ϕ)=0,

where V(ϕ) is the Sagdeev potential, which is given as

(33)V(ϕ)=A6Bϕ3U2Bϕ2+ρ0ϕ12E02

ρ0 and E0 are the integration constants representing the charge density and the electric field, respectively. U=U1C and E02/2 are the total energy of oscillations. By using the initial boundary conditions ϕ(0)=ϕ0 and ϕ(0)η=0, the expression for electric field E0 is obtained as

(34)E02=A3Bϕ03UBϕ02+2ρ0ϕ0

Substituting (33) and (34) in (32) and after some mathematical calculations, we obtain

(35)(ϕη)2=A3B(ϕ0ϕ)(ϕϕ1)(ϕϕ2)

where ϕ1 and ϕ2 are the real roots of Sagdeev potential, which are given as

(36)ϕ1,2=32[UAϕ03±12(β1ϕ0)(ϕ0β2)],

and

(37)β1,2=UA±2U2A22ρ0BA,

To find the periodic wave solution of IA waves, the inequalities β2ϕ0β1 or β1ϕ0β2 must be satisfied. The periodic wave solution (31) of [30] is given as

(38)ϕ(η)=ϕ1+ϕcncn2(Gη,m)

where cn is the Jacobian elliptic function, and the parameters m and G are defined as m2=(ϕ0ϕ1)(ϕ0ϕ2) and G=A12B(ϕ0ϕ2), where m2 is varied as 0m1. The amplitude and the wavelength of the cnoidal wave are defined as ϕcn=(ϕ0ϕ1) and λ=43BA(ϕ0ϕ1)K(m), where K(m) is the complete elliptic integral of the first kind.

For the case, m → 1 (E0=ρ0=0) at ϕ1=ϕ2=0, ϕcn=ϕ0=3UA=ϕm, G=(Aϕ012B)1/2=(U4B)1/2=1w, and cnδ=sechδ, the cnoidal wave solution reduces to the following solitary wave solution [31]

(39)ϕ(η)=ϕmsech2(ηw),

where ϕm=3UA and w=4BU are the peak amplitude and the width of the IA solitary waves, respectively. It is evident from (21), (31), and (39) that the propagation characteristics of IA cnoidal and solitary waves show the dependence on various plasma parameters such as density polarisation ratio (via μ), direction cosine (via kz), quantum parameter (via H), and ion-cyclotron frequency Ω.

5 Numerical Analysis

In the previous section, we have described the dynamics of IA cnoidal waves analytically in a magnetised quantum plasma in the presence of degenerate electrons with separated spin-up and spin-down population. In the present investigation, numerically it is confirmed that the nonlinear coefficient (A) and dispersion coefficient (B) of (31) are always positive. Hence, there is the formation of only compressive (positive potential) IA cnoidal and solitary waves. However, in the previous study, Ahmad et al. [20] derived the Zarkharov-Kuznetsov equation and obtained only negative potential solitary structures. In other study, Hussain and Mahmood [21] studied the oscillatory and monotonic shock structures by varying the spin-polarisation density. We have focused our analysis to study numerically the effects of various plasma parameters and spin polarisation density on the propagation characteristics of IA cnoidal waves and also discussed the existence of solitary waves in the limiting case in the given magnetised quantum degenerate plasma. In our present work, we have mainly analysed the effects of density of spin-up and spin-down electrons through the density polarisation ratio (i.e. μ) on the maximum amplitude of cnoidal waves. For numerical analysis, the plasma parameters are chosen as n0=10261029cm3, B0=103104G, and the temperature of the system is 105K in the dense astrophysical plasma regions of white dwarfs [32], [33]. We have chosen H = 0.2 for various plots with magnetic field intensity B0=104G, Fermi temperature TF=105K, and plasma density n0=1027cm3.

Figure 1 presents the variation of phase velocity vph with μ for different values of direction cosine (via kz). It has been found that with increase in the value of μ, the phase velocity (vph) decreases, whereas as the value of kz rises, there is enhancement in the value of phase velocity. It can be well explained from the expression of phase velocity given by (21). The phase velocity is directly proportional to the direction cosine kz and inversely proportional to the density polarisation ratio, which lead to the aforementioned variation in the characteristics of phase velocity. The variation in the maximum amplitude of cnoidal waves and the depth of corresponding Sagdeev potential is analysed using (33) for different values of ion-cyclotron frequency (Ω) with μ = 0 and is depicted in Figure 2a. Clearly, the solid (black) curve represents the Sagdeev potential corresponding to IA solitary waves with ρ0 = 0 and E0 = 0, whereas for finite value of ρ0(=0.3) and E0(=0.2), the other curves represent the Sagdeev potential corresponding to the IA cnoidal waves in a quantum degenerate plasma. The roots ϕ0, ϕ1, and ϕ2 given in (38) depict the motion of the pseudoparticle. The domain between the roots ϕ1<ϕ<ϕ2 confines the back-and-forth periodic oscillations of the pseudoparticle, and it cannot reach point ϕ0 because of the potential barrier. In Figure 2a, the rise in the value of Ω enhances the maximum amplitude and depth of the Sagdeev potential V(ϕ). This occurs due to the fact that V(ϕ) is dependent on the dispersion coefficient (B), which contains the quantity ion-cyclotron frequency (Ω). The maximum amplitude of IA cnoidal wave increases as Ω (ion-cyclotron frequency) increases, and the depth of the Sagdeev potential increases along the negative axis of V(ϕ). As V(ϕ) is dependent on coefficient B (which is sensitive to any variation in Ω), any change in Ω leads to change in B, which subsequently changes V(ϕ). In this case, when Ω increases, B decreases [see (33)]. With the small value of B in the denominator of (33) makes V(ϕ) larger. Also, the Sagdeev potential V(ϕ) corresponding to the IA cnoidal waves does not vanish at ϕ = 0, whereas for the IA solitary wave V(ϕ)|ϕ=0=0 (for ρ0 = 0 and E0 = 0). The phase plot, which shows the similar variation in the characteristics of IA cnoidal waves with the varying value of ion-cyclotron frequency (via Ω), is shown in Figure 2b. The black solid curve, called the separatrix, represents the characteristics of IA solitary waves, whereas all other curves lying inside the separatrix are characterised as the IA cnoidal waves. The pseudoparticle for IA cnoidal wave oscillates between the two real zeros of the Sagdeev potential, ϕ1 and ϕ2, as clearly depicted in Figure 2a. This means that potential structure is reappearing, and the distance between recurrences of wave shapes corresponds to one wavelength. Also, in Figure 3a, Sagdeev potential V(ϕ) is plotted against ϕ for μ = 1 for different values of ion-cyclotron frequency Ω. From Figure 3, we have depicted that with the increase in value of Ω, there is also enhancement in both maximum amplitude of cnoidal waves and depth of the Sagdeev potential. On comparing both Figures 2 and 3, it is observed that the maximum amplitude of cnoidal waves for zero spin polarisation case μ = 0 is smaller than the maximum amplitude of cnoidal waves for strongly spin-polarised case μ = 1.

Figure 1: (Colour online) Variation of the phase velocity vph vs. μ for the different values of kz. Solid (black) curve: kz=0.75${k_{z}}=0.75$; dotted (red) curve: kz=0.85${k_{z}}=0.85$; dashed (blue) curve: kz=0.95${k_{z}}=0.95$.
Figure 1:

(Colour online) Variation of the phase velocity vph vs. μ for the different values of kz. Solid (black) curve: kz=0.75; dotted (red) curve: kz=0.85; dashed (blue) curve: kz=0.95.

Figure 2: (Colour online) Variation of (a) Sagdeev potential of IA cnoidal waves V(ϕ) vs. ϕ, (b) phase plot for ion IA cnoidal waves for different values of Ω with μ = 0, kz=0.30${k_{z}}=0.30$, H = 0.2, U = 0.3, E0 = 0.2, and ρ = 0.3. Dashed (red) curve: Ω = 0.35; dotted (blue) curve: Ω = 0.37; dot-dashed (green) curve: Ω = 0.39; solid (black) curve: soliton with Ω = 0.35, ρ0 = 0 and E0 = 0.
Figure 2:

(Colour online) Variation of (a) Sagdeev potential of IA cnoidal waves V(ϕ) vs. ϕ, (b) phase plot for ion IA cnoidal waves for different values of Ω with μ = 0, kz=0.30, H = 0.2, U = 0.3, E0 = 0.2, and ρ = 0.3. Dashed (red) curve: Ω = 0.35; dotted (blue) curve: Ω = 0.37; dot-dashed (green) curve: Ω = 0.39; solid (black) curve: soliton with Ω = 0.35, ρ0 = 0 and E0 = 0.

Figure 3: (Colour online) Variation of (a) Sagdeev potential of IA cnoidal waves V(ϕ) vs. ϕ, (b) phase plot for ion IA cnoidal waves for different values of Ω with μ = 1, kz=0.30${k_{z}}=0.30$, H = 0.2, U = 0.3, and ρ = 0.2. Dashed (red) curve: Ω = 0.35; dotted (blue) curve: Ω = 0.37; dot-dashed (green) curve: Ω = 0.39; solid (black) curve: soliton with Ω = 0.35, ρ0 = 0, and E0 = 0.
Figure 3:

(Colour online) Variation of (a) Sagdeev potential of IA cnoidal waves V(ϕ) vs. ϕ, (b) phase plot for ion IA cnoidal waves for different values of Ω with μ = 1, kz=0.30, H = 0.2, U = 0.3, and ρ = 0.2. Dashed (red) curve: Ω = 0.35; dotted (blue) curve: Ω = 0.37; dot-dashed (green) curve: Ω = 0.39; solid (black) curve: soliton with Ω = 0.35, ρ0 = 0, and E0 = 0.

The influence of the varying value of the parameter kz (i.e. direction cosine) on the characteristics of Sagdeev potential and phase plane plot is depicted in Figure 4. With the increase in the value of kz, both the maximum amplitude and depth of the Sagdeev potential are reduced as shown in Figure 4a. The phase plane plot for the different values of kz is shown in Figure 4b. It is clearly visible from both figures that for ρ00 and E00, the phase curves show the repetitive behaviour on the same path. The pseudoparticle gets reflected back due to the potential force and starts oscillating between the two points whenever dV(ϕ)/dη=0 (i.e. when the pseudoparticle velocity becomes zero) as dV(ϕ)/dϕ does not get vanished. In physical space, one wavelength is equivalent to one cycle of the phase plot completed by the pseudoparticle. Hence, as the value of direction cosine kz rises, there is fall in the wavelength of IA cnoidal waves as depicted in Figure 4b. Figure 5a depicts the variation of Sagdeev potential for the different values of kz for μ = 1. Here, with the increase in value of kz, the maximum amplitude and depth of the Sagdeev potential of IA cnoidal wave are decreased. But again, on comparing both Figures 4 and 5, there is significant variation in the maximum amplitude of Sagdeev potential of IA cnoidal wave. The maximum amplitude of Sagdeev potential of IA cnoidal wave for zero spin polarisation ratio μ = 0 is smaller than the case for spin polarisation ratio μ = 1.

Figure 4: (Colour online) Variation of (a) Sagdeev potential of IA cnoidal waves V(ϕ) vs. ϕ, (b) phase plot for ion IA cnoidal waves for different values of direction cosine kz with μ = 0, Ω = 0.35, U = 0.3, H = 0.2, E0 = 0.2, and ρ = 0.3. Dashed (red) curve: kz=0.30${k_{z}}=0.30$; dotted (blue) curve: kz=0.32${k_{z}}=0.32$; dot-dashed (green) curve: kz=0.34${k_{z}}=0.34$; solid (black) curve: soliton with kz=0.3${k_{z}}=0.3$, ρ0 = 0, and E0 = 0.
Figure 4:

(Colour online) Variation of (a) Sagdeev potential of IA cnoidal waves V(ϕ) vs. ϕ, (b) phase plot for ion IA cnoidal waves for different values of direction cosine kz with μ = 0, Ω = 0.35, U = 0.3, H = 0.2, E0 = 0.2, and ρ = 0.3. Dashed (red) curve: kz=0.30; dotted (blue) curve: kz=0.32; dot-dashed (green) curve: kz=0.34; solid (black) curve: soliton with kz=0.3, ρ0 = 0, and E0 = 0.

Figure 5: (Colour online) Variation of (a) Sagdeev potential of IA cnoidal wave V(ϕ) vs. ϕ (b), phase plot for ion IA cnoidal waves for different values of direction cosine kz with μ = 1, Ω = 0.35, U = 0.3, H = 0.2, E0 = 0.2, and ρ = 0.2. Dashed (red) curve: kz=0.30${k_{z}}=0.30$; dotted (blue) curve: kz=0.32${k_{z}}=0.32$; dot-dashed (green) curve: kz=0.34${k_{z}}=0.34$; solid (black) curve: soliton with kz=0.3${k_{z}}=0.3$, ρ0 = 0, and E0 = 0.
Figure 5:

(Colour online) Variation of (a) Sagdeev potential of IA cnoidal wave V(ϕ) vs. ϕ (b), phase plot for ion IA cnoidal waves for different values of direction cosine kz with μ = 1, Ω = 0.35, U = 0.3, H = 0.2, E0 = 0.2, and ρ = 0.2. Dashed (red) curve: kz=0.30; dotted (blue) curve: kz=0.32; dot-dashed (green) curve: kz=0.34; solid (black) curve: soliton with kz=0.3, ρ0 = 0, and E0 = 0.

Figure 6a and b illustrate the variation of pulse profile of positive potential IA cnoidal waves for different values of ion-cyclotron frequency (via Ω) for μ = 0 and μ = 1, respectively. The maximum amplitude of IA cnoidal waves increases with the increasing value of Ω. It can be seen that the amplitude for μ = 1 is larger as compared to the amplitude of pulse profile for μ = 0. In the limiting case (for ρ0 = 0 and E0 = 0), we have also depicted the pulse profile of IA solitary waves in both figures. This can be verified from the fact that (33) is dependent on dispersion coefficient B, which is sensitive to ion-cyclotron frequency Ω, and any change in Ω causes a variation in dispersion coefficient B, which leads to the change in the maximum amplitude of IA cnoidal wave profile. The solid (black) curve represents the wave profile of IA solitary waves, which does not repeat itself like IA cnoidal waves and hence represents the characteristics of IA solitary waves. It is remarked that density spin polarisation has a significant influence on the propagation characteristics of cnoidal as well as solitary waves.

Figure 6: (Colour online) Variation of pulse profile of IA cnoidal wave (a) μ = 0, (b) μ = 1 for different values of Ω with U = 0.3, H = 0.2, E0 = 0.2, and ρ0=0.3${\rho_{0}}=0.3$. Dashed (red) curve: Ω = 0.35; dotted (blue) curve: Ω = 0.37; dot-dashed (green) curve: Ω = 0.39; solid (black) curve: soliton with Ω = 0.35, ρ0 = 0, and E0 = 0.
Figure 6:

(Colour online) Variation of pulse profile of IA cnoidal wave (a) μ = 0, (b) μ = 1 for different values of Ω with U = 0.3, H = 0.2, E0 = 0.2, and ρ0=0.3. Dashed (red) curve: Ω = 0.35; dotted (blue) curve: Ω = 0.37; dot-dashed (green) curve: Ω = 0.39; solid (black) curve: soliton with Ω = 0.35, ρ0 = 0, and E0 = 0.

6 Conclusions

The propagation characteristics of IA cnoidal waves are studied in the presence of magnetic field considering degenerate electrons with spin-up and spin-down states in quantum plasmas. We have considered that the Fermi step of the spin-up electrons is shorter than the Fermi step of the spin-down electrons, and due to this fact, there is no contribution of the outer-species collisions. Also, we are considering spin-up and spin-down electrons as different species, and their direction is related to a preferable direction in space for uniform external magnetic field, which lies along the z-direction in the present case. The equations of state for the spin-up and spin-down electrons are different due to the presence of external magnetic field, which changes an equilibrium concentration of each species. For nonlinear analysis, we have derived KdV equation by using the reductive perturbation method. Only positive potential (compressive) IA cnoidal waves and in the limiting case solitary waves are observed in the given plasma system. The influence of density spin polarisation via μ [for strongly spin-polarised (μ=1), zero spin-polarisation density (μ=0)] and other plasma parameters such as ion-cyclotron frequency (via Ω), quantum parameter (via H), and direction cosine (via kz) has been analysed numerically on the dynamics of compressive cnoidal waves. The phase velocity for IA cnoidal waves rises with the increment in the value of kz. The maximum amplitude of IA cnoidal waves and depth of the Sagdeev potential are enhanced (reduced) with the rising value of Ω (kz). On comparing the variation of Sagdeev potential with the different values of Ω and kz for both cases of spin polarisation (i.e. for μ = 0 and μ = 1), we have concluded that there is growth in the maximum amplitude of IA cnoidal waves as the value of μ changes from μ = 0 to μ = 1. It is emphasised that spin polarisation plays a very vital role for the variation in the characteristics of IA cnoidal waves. The findings of the present investigation may have paramount importance in different environments of dense astrophysical plasma regions, specifically like neutron stars, white dwarfs, and so on.

Acknowledgement

The authors are thankful to the reviewers for their critical and positive comments. R.K. acknowledges DST, Government of India, under DST-Purse scheme for financial support. This work is also supported by DRS-II (SAP) No. F/530/17/DRS-II/2015 (SAP-1) University Grants Commission, New Delhi, India.

References

[1] F. Calvayrac, P. G. Reinhand, E. Suraud, and C. Ullrich, Phys. Rep. 337, 493 (2000).10.1016/S0370-1573(00)00043-0Search in Google Scholar

[2] L. Stenflo, P. K. Shukla, and M. Marklund, Europhys. Lett. 74, 844 (2006).10.1209/epl/i2006-10032-xSearch in Google Scholar

[3] Y. I. Salamin, S. X. Hu, Z. Hatsagortsyam, and C. H. Keital, Phys. Rep. 427, 41 (2006).10.1016/j.physrep.2006.01.002Search in Google Scholar

[4] A. K. Harding and D. Lai, Rep. Prog. Phys. 69, 2631 (2006).10.1088/0034-4885/69/9/R03Search in Google Scholar

[5] S. Chandrashekhar, Astrophys. J. 74, 81 (1931).10.1086/143324Search in Google Scholar

[6] S. Chandrashekhar, Philos. Mag. 11, 592 (1931).10.1080/14786443109461710Search in Google Scholar

[7] S. Chandrashekhar, Mon. Not. R. Astron. Soc. 170, 405 (1935).Search in Google Scholar

[8] P. K. Shukla and B. Eliasson, Rev. Mod. Phys. 83, 885 (2011).10.1103/RevModPhys.83.885Search in Google Scholar

[9] D. Bohm, Phys. Rev. 85, 166 (1952).10.1103/PhysRev.85.166Search in Google Scholar

[10] F. Haas, L. G. Garcia, J. Goedert, and G. Manfredi, Phys. Plasmas 10, 3858 (2003).10.1063/1.1609446Search in Google Scholar

[11] F. Haas, Phys. Plasmas 12, 062117 (2005).10.1063/1.1939947Search in Google Scholar

[12] P. K. Shukla, B. Eliasson, and L. Stenflo, Phys. Plasmas 19, 072302 (2002).10.1063/1.4731733Search in Google Scholar

[13] G. Charbrier, F. Douchin, and A. Y. Potekhin, J. Phys. Condens. Matter 14, 9133 (2002).10.1088/0953-8984/14/40/307Search in Google Scholar

[14] M. Marklund and G. Brodin, Phys. Rev. Lett. 98, 025001 (2007).10.1103/PhysRevLett.98.025001Search in Google Scholar PubMed

[15] G. Brodin and M. Marklund, New J. Phys. 9, 277 (2007).10.1088/1367-2630/9/8/277Search in Google Scholar

[16] A. Mushtaq and S. V. Vladimirov, Phys. Plasmas 17, 102310 (2010).10.1063/1.3493632Search in Google Scholar

[17] P. A. Andreev, Phys. Rev. E 91, 033111 (2015).10.1103/PhysRevE.91.033111Search in Google Scholar PubMed

[18] Z. Iqbal and P. A. Andreev, Phys. Plasmas 23, 062320 (2016).10.1063/1.4954908Search in Google Scholar

[19] B. M. Askerov and S. R. Figarova, Thermodynamics, Gibbs Method and Statistical Physics of Electron Gases, Springer-Verlag, Berlin 2010.10.1007/978-3-642-03171-7Search in Google Scholar

[20] R. Ahmad, N. Gul, M. Adnan, and F. Y. Khattak, Phys. Plasmas 23, 112112 (2016).10.1063/1.4967768Search in Google Scholar

[21] S. Hussain and S. Mahmood, Phys. Plasmas 24, 102106 (2017).10.1063/1.4987002Search in Google Scholar

[22] K. Konno, T. Mitsuhashi, and Y. H. Ichikawa, J. Phys. Soc. Jpn. 46, 1907 (1979).10.1143/JPSJ.46.1907Search in Google Scholar

[23] A. V. Gurevich and L. Stenflo, Phys. Scripta 38, 855 (1988).10.1088/0031-8949/38/6/017Search in Google Scholar

[24] U. Kauschke and H. Schluter, Plasma Phys. Contr. F. 33, 1309 (1991).10.1088/0741-3335/33/11/006Search in Google Scholar

[25] E. F. El-Shamy, Phys. Rev. E 91, 033105 (2015).10.1103/PhysRevE.91.033105Search in Google Scholar PubMed

[26] H. Ur-Rehman, S. Mahmood, and S. Hussain, Phys. Plasmas 24, 012106 (2017).10.1063/1.4973836Search in Google Scholar

[27] N. Kaur, M. Singh, R. Kohli, and N. S. Saini, IEEE T. Plasma Sci. 46, 768 (2018).10.1109/TPS.2017.2760099Search in Google Scholar

[28] N. Kaur, M. Singh, and N. S. Saini, Phys. Plasmas 25, 043704 (2018).10.1063/1.5011948Search in Google Scholar

[29] T. Kaladze, S. Mahmood and H. Ur-Rehman, Phys. Scripta 86, 035506 (2012).10.1088/0031-8949/86/03/035506Search in Google Scholar

[30] T. Kaladze and S. Mahmood, Phys. Plasmas 21, 032306 (2016).10.1063/1.4868228Search in Google Scholar

[31] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Function, Dover, New York 1972.Search in Google Scholar

[32] D. Koester and G. Chanmugam, Rep. Prog. Phys. 53, 837 (1990).10.1088/0034-4885/53/7/001Search in Google Scholar

[33] G. Chabrier, D. Saumon, and A. Y. Potekhin, J. Phys. A Math. Gen. 39, 4411 (2006).10.1088/0305-4470/39/17/S16Search in Google Scholar

Received: 2019-04-26
Accepted: 2019-11-04
Published Online: 2019-12-07
Published in Print: 2020-02-25

©2020 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 25.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/zna-2019-0140/html
Scroll to top button