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
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-
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 e − i 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 (
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,
Electrons are inertialess, and the momentum equations of electrons with spin-up (
The Poisson’s equation is given as
The Fermi pressure of degenerate electrons with spin-up
In component form, the normalised continuity and momentum equations for ions are written as
The component forms of the normalised momentum equations for spin-up and spin-down electrons are written as
The normalised Poisson’s equation is
where
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 e − i plasmas. Therefore,
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
where
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:
From the z-component of the ion momentum equation:
From the z-component of electron momentum equation having the spin-up state:
From the z-component of the electron momentum equation having the spin-down state:
From Poisson’s equation:
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
On integrating (16)–(19) and applying some algebraic manipulations, we get the following relations:
where
Collecting the next higher-order terms of ε in the ion continuity equation, we have
Collecting the next higher-order terms of the z-component of the momentum equation of ions, we have
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:
The next higher-order terms of the Poisson’s equation
On simplifying (26)–(29) by eliminating the second-order quantities and making use of the boundary conditions
where the nonlinear coefficient A, dispersion coefficient B, and the coefficient C are given as
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
where
ρ0 and E0 are the integration constants representing the charge density and the electric field, respectively.
Substituting (33) and (34) in (32) and after some mathematical calculations, we obtain
where ϕ1 and ϕ2 are the real roots of Sagdeev potential, which are given as
and
To find the periodic wave solution of IA waves, the inequalities
where cn is the Jacobian elliptic function, and the parameters m and G are defined as
For the case, m → 1 (
where
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
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

(Colour online) Variation of the phase velocity vph vs. μ for the different values of kz. Solid (black) curve:

(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,

(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,
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

(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:

(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:
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.

(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
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
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
©2020 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Frontmatter
- Atomic, Molecular & Chemical Physics
- Numerical Investigation of the Cooling Temperature of the InGaP/InGaAs/Ge Subcells Under the Concentrated Illumination
- Dynamical Systems & Nonlinear Phenomena
- Ion-Acoustic Cnoidal Waves with the Density Effect of Spin-up and Spin-down Degenerate Electrons in a Dense Astrophysical Plasma
- Hydrodynamics
- Landau Quantised Modification of Rayleigh–Taylor Instability in Dense Plasmas
- Interaction of a Singular Surface with a Characteristic Shock in a Relaxing Gas with Dust Particles
- Quantum Theory
- Path Integrals, Spontaneous Localisation, and the Classical Limit
- Proposal for a New Quantum Theory of Gravity III: Equations for Quantum Gravity, and the Origin of Spontaneous Localisation
- Quantum-Phase-Field: From de Broglie–Bohm Double-Solution Program to Doublon Networks
- Solid State Physics & Materials Science
- Photovoltaic Generator Based on Laser-Induced Reversible Aggregation of Magnetic Nanoparticles
- Thermodynamics & Statistical Physics
- Investigation of the Finite Size Properties of the Ising Model Under Various Boundary Conditions
Articles in the same Issue
- Frontmatter
- Atomic, Molecular & Chemical Physics
- Numerical Investigation of the Cooling Temperature of the InGaP/InGaAs/Ge Subcells Under the Concentrated Illumination
- Dynamical Systems & Nonlinear Phenomena
- Ion-Acoustic Cnoidal Waves with the Density Effect of Spin-up and Spin-down Degenerate Electrons in a Dense Astrophysical Plasma
- Hydrodynamics
- Landau Quantised Modification of Rayleigh–Taylor Instability in Dense Plasmas
- Interaction of a Singular Surface with a Characteristic Shock in a Relaxing Gas with Dust Particles
- Quantum Theory
- Path Integrals, Spontaneous Localisation, and the Classical Limit
- Proposal for a New Quantum Theory of Gravity III: Equations for Quantum Gravity, and the Origin of Spontaneous Localisation
- Quantum-Phase-Field: From de Broglie–Bohm Double-Solution Program to Doublon Networks
- Solid State Physics & Materials Science
- Photovoltaic Generator Based on Laser-Induced Reversible Aggregation of Magnetic Nanoparticles
- Thermodynamics & Statistical Physics
- Investigation of the Finite Size Properties of the Ising Model Under Various Boundary Conditions