Propagation of shock structures in a high density plasma
Abstract
A theoretical investigation has been made to study the cylindrical
and spherical electronacoustic shock waves (EASWs) in an
unmagnetized, collisionless degenerate quantum plasma system
containing two distinct groups of electrons (one inertial
nonrelativistic cold electrons and other inertialess
ultrarelativistic hot electrons) and positively charged static
ions. By employing well known reductive perturbation method the
modified Burgers (mB) equation has been derived. It is seen that
only rarefactive shock waves can propagate in such a quantum
plasma system. The effects of degenerate plasma pressure and
number density of hot and cold electron fluids, nonplanar
geometry, and positively charged static ions are responsible to
modify the fundamental properties of EASWs. It is also observed
that the properties of planar mB shocks are quite different from
those of nonplanar mB shocks. The findings of the present
investigation should be useful in understanding the nonlinear
phenomena associated with nonplanar EAWs in both space and laboratory plasmas.
Keywords: Electronacoustic waves, modified Burgers equation, Shock waves, Degenerate pressure, Relativity, Compact objects
I Introduction
The apparent interest for quantum plasmas has been retained their interest due to their existence both in laboratory and space plasma environments Haas2003 ; Shukla2008 ; Eliasson2006 ; Khan2007 ; Misra2007 ; Sah2009b ; Hossain2013 ; a1 ; a2 ; a3 ; a4 . Generally, electronacoustic (EA) waves (EAWs) occur in a plasma environments (white dwarfs, neutron star, etc.) consisting of two distinct temperature electrons (as hot and cold electrons) Watanabe1977 ; Yu1983 ; Tokar1984 ; Mace1990 ; Kourakis2004 ; Mamun2004 . EAWs are nothing but highfrequency electrostatic mode, for which inertia is provided by the cold electron motion, while the restoring force comes from the hot electron thermal pressure. The positively charged ions may be safely assumed to be stationary, simply maintaining the quasineutrality condition of the plasma system. In such acoustic mode, the frequency lies in the range between the plasma frequency of the cold and hot electron fluids. Watanabe and Taniuti Watanabe1977 have first shown the existence of the electronacoustic (EA) mode in a plasma of twotemperature (cold and hot) electrons. Some past decades, EA waves has received a great deal of renewed interest not only because of the two distinct group of electron plasma is very common in laboratory experiments Derfler1969 ; Henry1972 ; Kadomtsev1971 ; Armstrong1979 ; Sheridan1991 and in space Dubouloz1993 ; Pottelette1999 ; Singh2001 ; Ang2007 ; Barnes2003 ; Fedele2002a ; Fedele2002b but also because of the potential importance of the EA waves in interpreting electrostatic component of the broadband electrostatic noise (BEN) observed in the cusp of the terrestrial magnetosphere Tokar1984 , in the geomagnetic tail Schriver1989 , in white dwarfs and neutron stars Shapiro1983 , etc.
Now a days, researchers of plasma community gives great attention to study the nonlinear behavior of astrophysical compact objects e.g. white dwarfs, neutron stars, etc. The plasma particle number density for such compact objects is so high (in white dwarfs it can be of the order of , even more) Hossen2014 ; Hossen2014c ; Hossen2014d ; Hossen2014e that the electron Fermi energy is comparable to the electron mass energy and the electron speed is comparable to the speed of light in a vacuum. Chandrasekhar Chandrasekhar1931a ; Chandrasekhar1931b presented a general expression for the relativistic ion and electron pressures in his classical papers. The pressure for electron fluid can be given by the following equation
(1) 
where is the electron number density and
(2) 
for the nonrelativistic limit (where , and is the Planck constant divided by ). And
(3) 
where
(4) 
for the ultrarelativistic limit Chandrasekhar1931a ; Chandrasekhar1931b ; Chandrasekhar1935 ; Chandrasekhar1939 ; Mamun2010a ; Mamun2010b ; Hossen2014b ; Hossen2014m ; Hossen2014n ; Hossen2014o .
A large number number of works on relativistic degenerate quantum plasma have been accomplished considering different acoustic waves in the recent years Shah2015a ; Shah2015b ; Shah2015c ; Shah2015d ; Hossen2014i ; Hossen2015c ; Hossen2014k ; Ema2015 ; H2014a ; H2014b ; H2014c ; H2014d ; H2014e ; H2015a ; H2016a ; H2016b . Han et al. Han2013 investigate the existence of electronacoustic shock waves and their interactions in a nonMaxwellian plasma with qnonextensive distributed electrons . Later on, Han et al. Han2014 theoretically investigated the nonlinear electronacoustic solitary and shock waves in a dissipative, nonplanar space plasma with superthermal hot electrons. Sahu and Tribeche Sahu2013 considered electron acoustic shock waves (EASWs) in an unmagnetized plasma whose constituents are cold electrons, immobile ions and Boltzmann distributed hot electrons and studied the effects of several parameters and ion kinematic viscosity on the basic features of EA shock waves. By considering quantum plasma ElLabany et al. ELLabany2013 investigated the effects of Bohm potential on the head on collision between two quantum electronacoustic solitary waves using the extended Poincar LighthillKuo method. Mahmood and Masood Mahmood2008 illustrated that an increase in quantum diffraction parameter broadens the nonlinear structure. Recently, Sah Sah2009 demonstrated that the width, the amplitude, and the velocity of electronacoustic double layers, in three component dense quantum plasmas consisting of stationary background ions and two electron populations: one cold and the other hot , are significantly affected by the ratio of unperturbed cold to hot electron densities. Again, the effect of static ions is very common in plasma physics literature Pakzad2011 ; Tribeche2010 ; Amour2012 . To the best of our knowledge, none of the authors did consider the combine effects of nonplanar geometry, effects of relativistic limits (i.e., both nonrelativistic and ultrarelativistic) and degenerate plasma pressure which can significantly modify the propagation of solitary and shock waves.
Ii Governing Equations
We consider a cylindrical and spherical EA waves in an unmagnetized, collisionless plasma, which is composed of nonrelativistic inertial cold electrons, both nonrelativistic and ultrarelativistic degenerate hot electron fluids, and static positive ions. Thus at equilibrium, we have , where is the equilibrium number density of the species ( for cold electrons, hot electrons, positive ions, respectively). The nonlinear dynamics of the electrostatic waves propagating in such a degenerate quantum plasma system is governed by the following normalized equations
(5)  
(6)  
(7)  
(8)  
(9) 
where for one dimensional planar geometry, (2) for nonplanar cylindrical (spherical) geometry, (s=c, h, i) is the the plasma species number density normalized by its equilibrium value , is the plasma fluid speed normalized by with () being the hot electron (cold electron) rest mass, is the speed of light in vacuum, is the electrostatic wave potential normalized by . Here is the ratio of iontocold electron number density. The time variable () is normalized by , and the space variable () is normalized by . The coefficient of viscosity is a normalized quantity given by . We have defined and .
Iii Derivation of modified Burgers Equation
We derive a dynamical modified Burgers (mB) equation for the nonlinear propagation of the EA waves by using equations (5)(9). To do so, we employ a reductive perturbation technique to examine electrostatic perturbations propagating in the relativistic degenerate dense plasma system due to the effect of dissipation, we first introduce the stretched coordinates Maxon1974 :
(10) 
where is the wave phase speed ( with being the angular frequency and being the wave number), and is a smallness parameter measuring the weakness of the dissipation (). We expand the parameters , , , , and in power series of as:
(11)  
(12)  
(13)  
(14)  
(15) 
Now, expressing equations (5)(9) (using equation (10), in terms of and , and substituting equations (11)(15), one can easily develop different sets of equations in various powers of . To the lowest order in , we have: , , , , where and . The relation represents the dispersion relation as well as the phase speed for the EA type electrostatic waves in the degenerate quantum plasma under consideration.
To the next higher order in , we obtain a sets of equations
(16)  
(17)  
(18)  
(19) 
Iv Discussion and Results
In this section, our first intention to numerically analyze the Burgers equation. However, for clear understanding, we first briefly discuss about the stationary shock wave solution for equation (20) with , though the solution is similar for both IA and EA waves (excluding the values of A and B). We should note that for a large value of , the term is negligible. So, in our numerical analysis, we start with a large value of (viz. ), and at this large (negative) value of , we choose the stationary shock wave solution of equation (23) [without the term ] as our initial pulse. The stationary shock wave solution of this standard Burgers equation is obtained by considering a frame (moving with speed which is the ion fluid speed at equilibrium) and the solution is H2016m ; H2016n ; H2016o
(23) 
where and .
Now, we investigate the dynamical properties of EA shock waves in terms of the intrinsic parameters of our model, namely the iontocold electron density ratio and the (cold electron) kinematic viscosity . It is important to not that the dissipation term only depends on the electron kinematic viscosity . The shock profile is nothing but sudden increase or decrease of the permanent profile of the waves H2016p ; H2016q ; H2016r ; H2016s . It is observed that EA waves are significantly modified when cold electron being nonrelativistic degenerate () and hot electron being ultrarelativistic degenerate () than both cold and hot electron being nonrelativistic degenerate (). It is important to mention that modified Burgers equation derived here is valid only for the limits , and . It is also important to note that for () rarefactive shock waves are found but no compressive shock waves exist at . We have considered for our numerical analysis of EASWs for the plasma system under investigation here.
Figure 1 shows the variation of phase speed () with ion to cold electron number density ratio . It is found that the phase speed decreases with the increasing values of . It is expected as the phase speed (derived from this considered plasma) is higher for lower values of (see the expression of ). The variation of the rarefactive amplitude of shock structures for both planar and nonplanar geometry is shown in Fig. 2 where cold electron and hot electron being nonrelativistic degenerate. Figure 3 shows the variation of the rarefactive amplitude of shock structures for both planar and nonplanar geometry where cold electron being nonrelativistic degenerate and hot electron being ultrarelativistic degenerate. The cylindrical and spherical variation of the amplitude of EASWs for both nonrelativistic and ultrarelativistic limits is shown in Figs. 47. Finally, the results that we have found in this investigation can be summarized as follows:

The cylindrical and spherical plasma system under consideration supports only rarefactive shock waves with negative potential, but no compressive shock waves exist.

The fundamental properties of EASWs are found to be significantly modified by the relativistic parameters, nonplanar geometry and plasma particle number densities.

It is observed that the phase speed (Vp) of these EA shocks inversely proportional to the square root of ion to cold electron number densities ratio .

It is also found that the phase speed () of EA waves decreases with the increasing values of (see Fig. 1).

It is observed that the amplitude of the shock is maximum for the spherical geometry, intermediate for cylindrical geometry, while it is minimum for the planar geometry (see Figs 23).

The amplitude of shocks proportional to the fluid speed but inversely proportional to the constant A.
In conclusion, our simplified theoretical model represents a small yet steady step towards the rigorous understanding of the behavior of cylindrical and spherical EA shocks in degenerate plasma environments, which appear to be of fundamental importance in a wide range of astrophysical Dubouloz1993 ; Pottelette1999 ; Singh2001 ; Ang2007 ; Barnes2003 ; Fedele2002a ; Fedele2002b and laboratory scenarios Derfler1969 ; Henry1972 ; Kadomtsev1971 ; Armstrong1979 ; Sheridan1991 .
References
 (1) F. Haas, L. G. Garcia, J. Goedert and G. Manfredi, Phys. Plasmas 10, 3858 (2003).
 (2) P. K. Shukla, J. Plasma Phys. 74, 107 (2008).
 (3) P. K. Shukla and B. Eliasson, Phys. Rev. Lett. 96, 245001 (2006).
 (4) S. A. Khan and A. Mushtaq, Phys. Plasmas 14, 083703 (2007).
 (5) A. P. Misra, P. K. Shukla and C. Bhowmik, Phys. Plasmas 14, 082309 (2007).
 (6) O. P. Sah and J. Manta, Phys. Plasmas 16, 032304 (2009).
 (7) M. Hasan, M. M. Hossain and A. A. Mamun, Astrophys. Space Sci. 345, 113 (2013).
 (8) B. Hosen, M. G. Shah, M. R. Hossen and A. A. Mamun, arXiv preprint arXiv:1707.05569 (2017).
 (9) M. Sarker, B. Hosen, M. R. Hossen and A. A. Mamun, arXiv preprint arXiv:1707.05571 (2017).
 (10) B. Hosen, M. G. Shah, M. R. Hossen and A. A. Mamun, arXiv preprint arXiv:1707.05567 (2017).
 (11) M. G. Shah, M. M. Rahman, M. R. Hossen and A. A. Mamun, arXiv preprint arXiv:1707.05568 (2017).
 (12) K. Watanabe and T. Taniuti, J. Phys. Soc. Jpn. 43, 1819 (1977).
 (13) M. Y. Yu and P. K. Shukla, J. Plasma Phys. 29, 409 (1983).
 (14) R. L. Tokar and S. P. Gary, Geophys. Res. Lett. 11, 1180 (1984).
 (15) R. L. Mace, M. A. Hellberg, J. Plasma Phys. 43, 239 (1990).
 (16) I. Kourakis, P. K. Shukla, Phys. Rev. E 69, 036411 (2004).
 (17) A. A. Mamun, P. K. Shukla and L. Stenflo, Phys. Plasmas 9, 1474 (2004).
 (18) H. Derfler and T. C. Simonen, Phys. Fluids 12, 269 (1969).
 (19) D. Henry and J. P. Treguier, J. Plasma Phys. 8, 311 (1972).
 (20) B. B. Kadomtsev and O. P. Pogutse, Nucl. Fusion 11, 67 (1971).
 (21) R. J. Armstrong, W. J. Weber and J. Trulsen, Phys. Lett. A 74, 319 (1979).
 (22) T. E. Sheridan, M. J. Goeckner, J. Goree and J. Vac. Sci. Technol. 9, 688 (1991).
 (23) N. Dubouloz et al., J. Geophys. Res. 98, 17415 (1993).
 (24) R. Pottelette et al., Geophys. Res. Lett. 26, 2629 (1999).
 (25) S. V. Singh and G. S. Lakhina, Planet. Space Sci. 49, 107 (2001).
 (26) L. K. Ang, P. Zhang, Phys. Rev. Lett. 98, 164802 (2007).
 (27) W. L. Barnes, A. Dereux and T. W. Ebbesen, Nature 424, 824 (2003).
 (28) R. Fedele and H. Schamel, Eur. Phys. J. B 27, 313 (2002).
 (29) R. Fedele, Phys. Scr. 65, 502 (2002).
 (30) D. Schriver and M. AshourAbdalla, Geophys. Res. Lett. 16, 899 (1989).
 (31) L. Shapiro and S. A. Teukolsky, Black Holes, White Dwarfs and Neutron Stars: The Physics of Compact Objects. John Wiley and Sons, New York, 1983.
 (32) M. R. Hossen, L. Nahar, S. Sultana and A. A. Mamun, Astrophys. Space Sci. 353, 123 (2014).
 (33) M. R. Hossen, L. Nahar, S. Sultana and A. A. Mamun, High Energ. Density Phys. 13, 13 (2014).
 (34) M. R. Hossen and A. A. Mamun, Braz. J. Phys. 44, 673 (2014).
 (35) M. R. Hossen, L. Nahar and A. A. Mamun, Phys. Scr. 89, 105603 (2014).
 (36) S. Chandrasekhar, Phi. Mag. 11, 592 (1931).
 (37) S. Chandrasekhar, Astrophys. J. 74, 81 (1931).
 (38) S. Chandrasekhar, Mon. Not. R. Astron. Soc. 170, 405 (1935).
 (39) S. Chandrasekhar, An Introduction to the Study of Stellar Structure. Dover Publication, New York, 1939.
 (40) P. K. Shukla and A. A. Mamun, Phys. Lett. A 324, 42380 (2010).
 (41) P. K. Shukla and A. A. Mamun, Phys. Plasmas 17, 1045040 (2010).
 (42) M. R. Hossen, L. Nahar and A. A. Mamun, Braz. J. Phys. 44, 638 (2014).
 (43) M. R. Hossen, L. Nahar and A. A. Mamun, J. Astrophys. 2014, 653065 (2014).
 (44) M. R. Hossen, Nonlinear excitations in degenerate quantum plasmas. Laplambert Academic Publishing Company, Germany, 2014 ISBN10: 3659597392.
 (45) M. R. Hossen, S. A. Ema and A. A. Mamun, Commu. Theor. Phys. 62, 888 (2014).
 (46) M. G. Shah, M. R. Hossen, S. Sultana and A. A. Mamun, Chin. Phys. Lett. 32, 085203 (2015).
 (47) M. G. Shah, M. R. Hossen and A. A. Mamun, Braz. J. Phys. 45, 219 (2015).
 (48) M. G. Shah, M. R. Hossen and A. A. Mamun, J. Plasma Phys. 81, 905810517 (2015).
 (49) M. G. Shah, A. A. Mamun and M. R. Hossen, J. Korean Phys. Soc. 66, 1239 (2015).
 (50) M. R. Hossen, L. Nahar and A. A. Mamun, J. Korean Phys. Soc. 65, 1863 (2014).
 (51) M. R. Hossen and A. A. Mamun, Plasma Sci. Technol. 17, 177 (2015).
 (52) M. R. Hossen and A. A. Mamun, J. Korean Phys. Soc. 65, 2045 (2014).
 (53) S. A. Ema, M. R. Hossen and A. A. Mamun, Contrib. Plasma Phys. 55 551 (2015).
 (54) M. R. Hossen, L. Nahar, S. Sultana and A. A. Mamun, Astrophys. Space Sci. 353 123 (2014).
 (55) M. R. Hossen, L. Nahar, S. Sultana and A. A. Mamun, High Energy Density Phys. 13, 13 (2014).
 (56) M. R. Hossen, L. Nahar and A. A. Mamun, Phys. Scr. 89, 105603 (2014).
 (57) M. R. Hossen and A. A. Mamun, Braz. J. Phys. 44, 673 (2014).
 (58) M. R. Hossen and A. A. Mamun, SOP Trans. Theor. Phys. 1, 57 (2014).
 (59) M. R. Hossen and A. A. Mamun, Braz. J. Phys. 45, 200 (2015).
 (60) B. Hosen, M. G. Shah, M. R. Hossen and A. A. Mamun, Eur. Phys. J. Plus 131, 81 (2016).
 (61) B. Hosen, M. Amina, M. R. Hossen and A. A. Mamun, J. Korean Phys. Soc. 69, 1762 (2016).
 (62) JiuNing Han, JunXiu Li, YongLin He, ZhenHai Han, GuangXing Dong and YaGong Nan, Phys. Plasmas 20, 072109 (2013).
 (63) JiuNing Han, WenShan Duan, JunXiu Li, YongLin He, JunHua Luo, YaGong Nan, ZhenHai Han and GuangXing Dong, Phys. Plasmas 21, 012102 (2014).
 (64) B. Sahu and M. Tribeche, Adv. Space Res. 51, 2353 (2013).
 (65) S. K. ELLabany, E. F. ELShamy and M. G. ElMahgoub, Astrophys. Space Sci. 339 195 (2012).
 (66) S. Mahmood and W. Masood, Phys. Plasmas 15, 122302 (2008).
 (67) O. P. Sah, Phys. Plasmas 16, 012105 (2009).
 (68) H. R. Pakzad and M. Tribeche, Astrophys. Space Sci. 334, 45 (2011).
 (69) M. Tribeche and L. Djebarni, Phys. Plasmas 17, 124502 (2010).
 (70) R. Amour, M. Tribeche and P. K. Shukla, Astrophys. Space Sci. 338, 287 (2012).
 (71) S. Maxon and J. Viecelli, Phys. Rev. Lett. 32, 4 (1974).
 (72) M. R. Hossen, S. A. Ema and A. A. Mamun, Chin. Phys. Lett. 33, 065203 (2016).
 (73) S. A. Ema, M. R. Hossen and A. A. Mamun, Contrib. Plasma Phys. 55, 596 (2015).
 (74) M. G. Shah, M. M. Rahman, M.R. Hossen and A.A. Mamun, Plasma Phys. Rep. 42, 168 (2016).
 (75) S. A. Ema, M. R. Hossen and A. A. Mamun, Phys. Plasmas 22, 092108 (2015).
 (76) M. R. Hossen, S. A. Ema and A. A. Mamun, IEEE Trans. Plasma Phys. 44, 492 (2016).
 (77) S. A. Ema, M. R. Hossen and A. A. Mamun, Plasma Phys. Rep. 42, 362 (2016).
 (78) M. G. Shah, M. M. Rahman, M. R. Hossen and A. A. Mamun, Commun Theor. Phys. 64, 208 (2015).