Abstract
The study of dense degenerate plasmas at extreme temperatures and densities motivates the investigation of relativistic and ultrarelativistic regimes. In this paper, we are presenting a comprehensive study of the linear and nonlinear propagation characteristics for kinetic Alfven acoustic waves (KAAW) at relativistic and ultra relativistic Fermi energies with small temperature corrections. The linear wave propagation of fast and slow modes for relativistic KAAW are studied and a comparative analysis of both the modes is plotted for relativistic, ultrarelativistic and nonrelativistic regimes. It has been investigated that for constant magnetic fields, the phase velocity of these waves decreases due to an increase in both the specific heat and thermal pressure at low plasma beta values and KAAW propagate more efficiently in outer regions with lower specific heat than in denser inner regions of neutron stars. Various theoretical approaches like Sagdeev pseudopotential approach, dynamical analysis, and finite amplitude expansion method are employed to analyse the nonlinear relativistic and ultrarelativistic KAAW. Phase space trajectories have been plotted to study the nonlinear structures corresponding to arbitrary amplitude perturbations. The existence regions of solitary kinetic Alfven acoustic waves (SKAAW) and the physical reasons for their dependence on different parameters have also been shown. Possible applications are envisaged in the intense laser-matter interactions, plasma photonics, dense astrophysical objects like the interior state of terrestrial planets, neutron stars, and white dwarfs where relativistic effects dominate the microscopic scale and may have macroscopic consequences.
Similar content being viewed by others
Introduction
Degenerate quantum plasma is a mature, yet dynamic, field of investigation. This captivating form of matter is ubiquitously studied on a microscale1,2. Interestingly, even while studying such exotic states3,4on an astrophysical scale5,6, its quantum effects are found to be significant and non-negligible. Quantum plasmas occur naturally in dense astrophysical environments like, active galactic nuclei, neutron stars7, white dwarfs8, and the inner region of black hole’s accretion disks. But it must be noted that due to their distance from the Earth, it becomes very difficult to collect their data, apart from a spectroscopic approach. Nevertheless, scientists have gathered substantial information by observations from space telescopes like Hubble9,10,11,12,13,14, James Webb15,16,17and X-ray observatories18,19,20, about black holes, neutron stars, and white dwarfs. However, there is still much more left to be discovered about these compact celestials as no in-situ data can be collected using the present-day technology. Therefore to study them in detail, they must be observed in laboratories21,22, where they are reported to exist in quantum dots, carbon nanotubes, micro-electronics, nanophononics23and intense laser-plasma interactions24. Such plasmas are characterized by extremely high densities and relatively low temperatures (\(\:T\ll\:{T}_{F}\))25. Due to the recent advancements in technology, dense plasma layer of micrometer order depth has been generated in the laboratory, with the average density of \(\:{10}^{29}{m}^{-3}\)and the electron temperatures at 2–4 keV23. This work has been verified at a virtual laser plasma laboratory26 by the usage of X-ray spectra and particle-in-cell (PIC) simulation, where all the formation phases were observed. This was done by the irradiation of femtosecond laser pulses of relativistic intensities on Au and Ni nanowires. Using PIC simulations, it was also predicted that this density can increase by one more order at \(\:{10}^{30}{m}^{-3}\), if the nanowires are made of Uranium, rather than Au or Ni. This tabletop laser plasma is eerily reminiscent to the dense astrophysical environments, and hence allows a closer peek at these celestial objects.
Degenerate plasma exhibits relativistic effects when the Fermi energy becomes comparable to the rest mass energy (\(\:{\epsilon\:}_{F}\:\sim\:{{m}_{o}c}^{2}\)), which happens due to an increase in mean energy of matter as a result of compression27. However, relativistic effects and their observations have always remained spellbinding for researchers28, in both astrophysical and laboratory environments. Nevertheless, intense plasma interactions are providing new avenues for research and development in the relativistic regimes29. This not only has direct applications in inertially confined fusion plasmas, but also for the study of ultra-high energy cosmic rays in astrophysics30. Needless to say, relativistic effects in plasmas are not limited to degeneracy only and can appear in other forms too. For instance, highly intense short laser pulses have been used to observe relativistic plasma waves31,32, aiding in the acceleration of the particles to relativistic energies. In astrophysical landscapes, relativistic and ultrarelativistic effects27plays a crucial role in shaping the dynamics of neutron stars, white dwarfs and swirling accretion disks of black holes, which comes about due to the extreme conditions of density, temperature19and magnetic fields33. It is imperative to note that thermonuclear fusion ceases to exist in a white dwarf and the star’s structure is held by a balance between the electron degeneracy pressure and the immense force due to gravity34. In these compressed states, relativistic degenerate electrons are surrounded by bare nuclei of carbon and other heavier elements due to which they can be modelled as a one component plasma35. Whereas neutron stars are held together by a balance between neutron degeneracy pressure and gravity36. Lenzen et al.37 have expressed the total number density of electrons from the nuclei density of condensed magnetic matter. The surface electron density at neutron star is proportional to \(\:{B}^{6/5}\), for example, \(\:{n}_{e}\approx\:{10}^{34}{m}^{-3}\) for \(\:B\approx\:{10}^{8}\:Tesla\)37.
Astrophysical environments are rich in various phenomena such as instabilities, anisotropies, high energy interaction, relativistic effects and X-ray bursts from stars etc. The existence of MHD waves in the surface layer of neutron stars are justified38,39. Hence it becomes crucial to deal with wave dispersions and nonlinearities to study plasma in such domains. In recent decades, substantial research has been conducted to investigate the aforementioned phenomena. Tsintsadze N. L. et al.40 showed that high frequency electrical fields in the transverse direction are not “steady” in a relativistic plasma. Using such instabilities, the existence of strong cosmic radiation can be explained. Flacco et al.41suggested strong cosmic magnetization is the consequence of ultra-high energy particles and studied magnetic fields driven by relativistic electrons in plasma. For an ultrarelativistic degenerate electron positron gas, the nonlinear screening effects have been investigated and it was observed that screening length was reduced due to strong spatial fluctuations in the potential field42. Due to the prevailing nonlinearities in the MHD equations, the presence of solitary structures becomes possible. Solitary structures are observed throughout nature43,44,45and are also present in Fermionic superfluids46. The particle nature of solitons has been verified experimentally46 and their existence is supported by theoretical models and experimental studies. Roy et al.47carried out the stability analysis of electromagnetic (EM) solitons and studied their dependence on degeneracy parameters for a dense relativistic laser plasma. Modulational instabilities of EM waves in partially degenerate relativistic plasma have also been carried out for nonlinear solitary structures48. El-Labany et al.49have shown simulation results of oblique interactions of solitons in a relativistic dense plasma. Trapped electron distribution50,51,52interestingly changes the order of nonlinearities in a plasma system, and this verified phenomenon53,54,55,56has also compelling implications in both the relativistic and ultrarelativistic regimes57. Following the work by Shah et al.57, Iqbal et al.58 investigated the effect of quantizing magnetic field on nonlinear ion acoustic waves and examined the solitary profiles by varying density, magnetic field, and Mach number.
Investigation of kinetic Alfven-acoustic waves (KAAW) in a partially degenerate relativistic plasma can find applications in laser-solid density plasma interactions23, superdense stars, and nano-scale photonics. Primarily, KAAW were theoretically studied for classical plasmas59,60,61. In the laboratory, KAAW have been investigated in spontaneous fluctuations of temperature filament62and in the experiments on magnetized helium plasmas63. In space, the Kelvin-Helmholtz instability64and 2D vortex structures in ionosphere observed by IC-B-1300 satellite65are considered as an effective source of KAAW. In dense astrophysical objects, research on Alfven waves highlights various phenomena. Studies explore spin effects on KAW66,67, SKAAW68in nonrelativistic dense plasma and Alfven waves driven winds in proto-neutron stars69. Observations on fast radio burst (FRBs) data from SGR 1935 + 2154 (SGR 1935 + 2154 is a galactic magnetar) consider these FRBs as a source of Alfven waves70. Simulations show their role in explaining quasi-periodic oscillations71in neutron star cores and for the generation of plasmoids72. These plasmoids are closed magnetic loops that accelerate and drive blast waves into magnetar winds. Observations connect Alfven wave activity with X-ray and radio emissions from magnetars and soft gamma repeaters73, providing insight into high-energy astrophysical events. Recently, Chen et al.74 proposed that instabilities in the magnetospheres of neutron stars could generate relativistic Alfven waves and boost the electric currents along the curved magnetic field lines.
In this paper, relativistic and ultra relativistic regimes of degenerate plasmas are considered using two-potential theory75. This approach is valid for low frequency (\(\:{\omega\:}^{2}<<{{\omega\:}_{ci}}^{2}\)) and low \(\:\beta\:\left(=\:2{{c}_{s}}^{2}/{{v}_{A}}^{2}\right)\) plasmas \(\:\left(1>\beta\:>\:{m}_{e}/{m}_{i}\right)\). At the outset in the following section, we furnish the set of MHD equations and present the linear dispersion relation of relativistic KAAW. We have also discussed the limiting cases for ultrarelativistic and nonrelativistic regimes. Next, in section III, we analyse the nonlinear dynamics of the relativistic KAAW by using Sagdeev pseudopotential approach, dynamical analysis, and finite amplitude expansion method in a stationary frame of reference and discussed different wave solutions. In section IV, the nonlinear analysis of ultrarelativistic KAAW is carried out to study the nonlinear structures. Finally, we consolidate our work and give the conclusions in section V.
Governing equations of magnetohydrodynamic (MHD) model
We study the propagation characteristics of KAAW for a dense plasma immersed in a uniform external magnetic field \(\:{\mathbf{B}}_{\mathbf{o}}={\text{B}}_{\text{o}}\widehat{\text{z}}\), where \(\:{\text{B}}_{\text{o}}\) is the magnetic field strength and \(\:\widehat{\text{z}}\) is the unit vector along the z-direction. Relativistic degeneracy implies that the Fermi energy of the electrons is comparable to the rest mass energy of electrons. Due to the finite ion Larmour radius effect, we consider electrons to be inertialess in comparison to heavy ions (\(\:{m}_{i}\gg\:{m}_{e}\)), which can be modelled classically. We choose the propagation vector \(\:\varvec{k}\) to lie in the x-z plane and apply the two-potential theory, such that, the gradients of the potentials \(\:\varnothing\:\) and \(\:\psi\:\) can be written as the components of electric field \(\:\varvec{E}\)75 and as pointed out in the preceding section, the two potential approximation is applicable in low \(\:\beta\:\) plasma. Most of the astrophysical plasmas are low \(\:\beta\:\) plasma. Thus, we have
The dynamics of low frequency KAAW76 are governed by the ion equation of motion,
The continuity equation and the low frequency approximation of Maxwell’s curl equations are given by
Where \(\:{n}_{\varvec{q}}\) and \(\:{\varvec{v}}_{\varvec{q}}(q=i,\:e)\) respectively be the density and velocity of the species. \(\:\varvec{j}(=\:e\left({n}_{i}{\varvec{v}}_{\varvec{i}}-\:{n}_{e}{\varvec{v}}_{\varvec{e}}\:\right))\)is the current density and the electron density for partially degenerate relativistic plasma is described by the following distribution57
We note that the expression \(\:\left(5\right)\) follows from the Fermi-Dirac distribution function and here \(\:{N}_{oR}\left(={\mu\:}^{3}/3{\pi\:}^{2}{\left(\hslash\:c\right)}^{3}\right)\) is the relativistic background density, \(\:{{\Psi\:}}_{R}\left(=e\:\psi\:/\mu\:\right)\) is the normalized potential, \(\:{\epsilon}_{o}\left(={m}_{e\:}{c}^{2}/\mu\:\right)\) is the total normalized energy, \(\:{T}_{R}\left(=T/\mu\:\right)\) is the normalized ambient temperature the presence of which indicates partial degeneracy, and \(\:\mu\:\left(={\epsilon\:}_{Fo}+{m}_{e\:}{c}^{2}\right)\) is the chemical potential for a relativistic quantum plasma. The subscript R denotes relativistic case. The Fermi energy is given by \(\:{\epsilon\:}_{Fo}=({\hslash\:}^{2}/2{m}_{e}){\left(3{\pi\:}^{2}{N}_{oR}\right)}^{2/3}\) and in the absence of perturbed potential the electron number density\(\:{\:n}_{eoR}\) is defined as,
At a certain point in a plasma system both the relativistic and nonrelativistic situations cannot coexist. However, as we move towards the core, the density increases and system transits from nonrelativistic to a relativistic regime.
In the given geometry, we linearize Eqs. (1–4) for the low frequency KAAW waves \(\:\left({\omega\:}^{2}<<{\varOmega\:}_{ci}^{2}\right)\) and obtain the components of momentum equation. The linearized ion and electron continuity equations come from quasi-neutrality condition \(\:({n}_{eoR}\approx\:\:{n}_{ioR})\) and Eq. (4) is solved for two scalar potentials \(\:\varnothing\:\) and \(\:\psi\:\). The current density \(\:{j}_{z}\) is obtained by using linearized electron continuity equation. Further, by eliminating the dimensionless potential \(\:{\varPhi\:}_{R}\left(=e\:\varnothing\:/\mu\:\right)\), we obtain the linear dispersion relation given by,
The linear dispersion relation of KAAW is crucial for understanding wave propagation, which is influenced by both the Alfven velocity and ion thermal velocity. We note that the two factors on the left-hand side of Eq. (7), individually set to zero when \(\:{\lambda\:}_{sFR}=0\) (i.e. when \(\:{k}_{x}\to\:0\)), give Alfven and ion acoustic wave dispersion equations for a relativistic quantum plasma. These two waves are coupled through the coupling parameter \(\:{\lambda\:}_{sFR}\left(={\:k}_{x}^{2}\:{C}_{sFR}^{2}/{\omega\:}_{ci}^{2}\right)\). Here, \(\:{C}_{sFR}\left(=\sqrt{\mu\:/{m}_{i}\:{A}_{R}}\right)\)is the relativistic Fermi ion sound velocity57, \(\:{\omega\:}_{ci}\left(=e{B}_{0}/{m}_{i}\right)\) is the ion cyclotron frequency and \(\:{v}_{AR}(={B}_{0}/{{\mu\:}_{o}\:n}_{eoR}\:{m}_{i})\) is the Alfven velocity for a relativistic quantum plasma. The factor \(\:{A}_{R}\) carries both the relativistic and small temperature effects, defined as
We obtain the frequencies of both the fast and slow modes of KAAW from the dispersion relation in Eq. (7) as,
The \(\:\pm\:\) sign corresponds to fast and slow modes respectively. Equation (9) shows that for perpendicular propagation i.e. at \(\:\theta\:={90}^{o}\), \(\:\omega\:/k\to\:\infty\:\) and there will be no wave momentum in the medium. Also, as \(\:\omega\:={\omega\:}_{ci}\) the group velocity goes to zero and wave dispersion will disappear.
The fast and slow features of relativistic KAAW are shown in Fig. 1. Figure 1(a) is plotted for the normalized wave frequency \(\:\left(\omega\:/{k\:v}_{AR}\right)\) against the normalized perpendicular wave number \(\:\left({k}_{x}/k\right)\) with \(\:k=\sqrt{{k}_{x}^{2}+{k}_{z}^{2}\:}\) and Fig. 1(b) is plotted for the normalized wave frequency \(\:\left(\omega\:/{\omega\:}_{ci}\right)\) against the normalized wave number \(\:\left(k/{k}_{o}\right)\), with \(\:{k}_{o}\) a reference wave number, for angle of propagation \(\:{85}^{o}\). It is observed that the wave frequencies are affected by the small temperature effects. At higher \(\:{T}_{R}\), the ion thermal velocity increases and modifies the phase velocity of KAAW.
Dispersion relation plots for slow and fast relativistic KAAW at \(\:{\text{T}}_{R}=0\) and \(\:{\text{T}}_{R}=0.5\). Other parameters are \(\:{n}_{0}=1\times\:{10}^{36}\:{m}^{-3}\), \(\:{B}_{0}=1\times\:{10}^{9}\text{T}\), and \(\:\:\theta\:={85}^{^\circ\:}\) (a) Normalized wave frequency \(\:\left(=\omega\:/{k\:v}_{AR}\right)\) against the normalized perpendicular wave number \(\:\left(={k}_{x}/k\right)\), (b) Normalized wave frequency \(\:\left(\omega\:/{\omega\:}_{ci}\right)\) against the normalized wave number \(\:\left(=k/{k}_{o}\right)\) at \(\:\theta\:={85}^{^\circ\:}\).
We now briefly consider the two limiting cases of ultrarelativistic and nonrelativistic regimes.
A. Ultrarelativistic case: In a superdense plasma, the Fermi energy becomes much greater than the rest mass energy of electrons, i.e. \(\:{\epsilon\:}_{Fo}\gg\:{m}_{e}{c}^{2}\), and degenerate plasma is characterized in ultrarelativistic regime. In this regime \(\:{\epsilon}_{o}\to\:0\) and Fermi energy is defined as \(\:{\epsilon\:}_{Fo}={cp}_{F}={\left(3{\pi\:}^{2}{n}_{0}\right)}^{1/3}\hslash\:c\). Thus, the distribution given by Eq. (5) becomes
Where \(\:{N}_{oU}=\frac{{{\epsilon\:}_{Fo}}^{3}}{3{\pi\:}^{2}{\left(hc\right)}^{3}}\) is the ultrarelativistic background density and in the absence of perturbed potential the electron number density is \(\:{n}_{eoU}={N}_{oU}\left[1+{\pi\:}^{2}{T}_{U}^{2}\right]\). Adopting the same procedure for \(\:{n}_{eU}\), the set of Eqs. (1–4) finally gives the linear dispersion relation for ultrarelativistic degenerate plasma, as
Where \(\:{\lambda\:}_{SFU}\left(={\:k}_{x}^{2}\:{C}_{sFU}^{2}/{\omega\:}_{ci}^{2}\right)\) is the coupling parameter, \(\:{c}_{sFU}\left(=\sqrt{{\epsilon\:}_{Fo}/{m}_{i}\:{A}_{U}}\right)\) is the ultrarelativistic Fermi ion sound velocity, \(\:{v}_{AU}(={B}_{0}/\sqrt{{{\mu\:}_{o}\:n}_{eoU}\:{m}_{i}})\) is the Alfven velocity for a relativistic quantum plasma, and \(\:{A}_{U}\) is the ultrarelativistic factor for a partially degenerate plasma, defined as
B. Nonrelativistic case: In the nonrelativistic case, the Fermi energy is much less than the rest mass energy of electrons, i.e. \(\:{\epsilon\:}_{Fo}\ll\:{m}_{e}{c}^{2}\) and we retrieve the results already presented by Sabeen et al. and Asam et al.68,77.
It is clearly seen that the relativistic, ultrarelativistic and nonrelativistic KAAW show dispersion for both the slow and fast modes in a degenerate plasma as shown in Fig. (2). These waves can traverse through ultra-high densities and transfer energy, parallel to the strong ambient magnetic fields along with the trail of charged particles. It also shows that, in constant magnetic fields, as plasma shifts from nonrelativistic to ultrarelativistic regimes, the phase velocity of KAAW decreases due to an increase in both the specific heat and thermal pressure at low plasma beta \(\:\left({\beta\:}_{FR}=\frac{2{{C}_{sFR}}^{2}}{{{v}_{AR}}^{2}}\right)\:\)values i.e. \(\:\frac{{m}_{e}}{{m}_{i}}<{\beta\:}_{FR}<1\:\). This intriguing result shows that KAAW propagate more efficiently in outer regions with lower specific heat than in denser inner regions of neutron stars.
Nonlinear analysis of relativistic KAAW
-
A.
Pseudopotential method.
We investigate the nonlinear spread of KAAW in a partially degenerate relativistic magnetoplasma by employing Sagdeev potential (pseudopotential) approach. In order to examine the existence conditions of nonlinear solitary profiles, we choose a normalized comoving frame of reference, such that \(\:\xi\:={\:k}_{x}\:x+{\:k}_{z}\:z-M\:t\)78. \(\:{\:k}_{x}\) and \(\:{\:k}_{z}\) are the directional cosines and \(\:M(=u/{C}_{sFR})\)is the Mach number. In this one-dimensional frame, the wave propagates obliquely to the ambient magnetic field \(\:{\varvec{B}}_{\varvec{o}}\) and will be stationary with respect to \(\:\xi\:\). In this stationary frame, the particle number densities are normalized by their respective equilibrium densities, velocities by relativistic Fermi ion sound velocity \(\:{C}_{sFR}\), scale length by ion gyroradius \(\:(={C}_{sFR}/{\omega\:}_{ci})\), and time scale by \(\:{\omega\:}_{ci}\). Thenceforth, the system of Eqs. (1–4) in the dimensionless frame \(\:\xi\:\) is as follows,
We integrate Eq. (15) by applying boundary conditions that as \(\:\xi\:\to\:\infty\:\), the perturbed quantities \(\:{v}_{x},{v}_{z},{{\Psi\:}}_{R}\to\:0\) and the number density \(\:{n}_{R}\to\:{\alpha\:}_{R}\), yields
Using Eqs. (5), (13) and (17), we get
\(Here \;\:{n}_{R}=\left(\raisebox{1ex}{${n}_{eR}$}\!\left/\:\!\raisebox{-1ex}{${N}_{oR}$}\right.\right),\:\:{\alpha\:}_{R}={\delta\:}_{2}^{-\frac{1}{2}}\:\left[\:{\delta\:}_{2}^{\:\:2}+{\delta\:}_{T}^{2}\:\left(2-{\epsilon\:}_{0}^{2}\right)\right]\:,\:\:\:\:{\delta\:}_{T}=\frac{\pi\:\:\:{T}_{R}}{\sqrt{2}}\:\:\:,\:\:{\delta\:}_{1}={\left\{\left(1+\left.{{\Psi\:}}_{R}\right)\right.\right.}^{2}-\left.{\epsilon}_{0}^{2}\right\}\)
.
Integrating Eq. (19) and applying the boundary conditions \(\:{{v}_{z},\:\psi\:}_{R}\to\:0\) as \(\:\xi\:\to\:\infty\:\), we obtain an expression for the parallel ion velocity \(\:{v}_{z}\)
The constant of integration
Solving Eqs. (14), (17) and (18) gives the following expression
We insert Eqs. (5) and (20) in the above expression and expand the term \(\:\raisebox{1ex}{${\alpha\:}_{R}$}\!\left/\:\!\raisebox{-1ex}{${n}_{R}$}\right.\) with respect to temperature. Here, we have ignored higher order terms and the terms containing \(\:{T}_{R}^{4}\) (for a quantum plasma \(\:{\text{T}}_{R}<1\)). Moreover, Eq. (22) is integrated to express through Sagdeev potential79 in the usual manner,
Therefore Eq. (22) yields,
\(\:{C}_{2}\) is the constant of integration and defined as,
The obtained expression of Sagdeev Pseudopotential \(\:V\left({{\Psi\:}}_{R}\right)\) in Eq. (24) does not give solitary KAAW except for fully degenerate relativistic plasma, i.e. \(\:{\text{T}}_{R}=0\) (see Figs. (3–5) in the next sub section). Therefore, for partially degenerate relativistic regime the nonlinear KAAW is further analysed by employing dynamical analysis.
-
B.
Autonomous Dynamical Analysis.
As seen in the previous subsection, the presence of fractional nonlinearities show that the Sagdeev potential approach does not yield solitary solution for KAAW and, therefore, it will be useful to perform a qualitative analysis of phase portraits in a dynamical system80,81. Such a system enables the understanding of wave trajectories in a phase space to provide a better analysis from the numerical solution. Considering \(\:Q={\partial\:}_{\xi\:}{{\Psi\:}}_{R}\), we get the following dynamical system
The Hamiltonian for the system carries effective potential82, which in our case is \(\:V\left({{\Psi\:}}_{R}\right)\). Here Hamiltonian function \(\:H\left({{\Psi\:}}_{R},W\right)\) is defined as,
The characteristics equation \(\:det\left(J-\lambda\:\:I\right)=0\) along with Jacobian matrix \(\:J=\left(\begin{array}{cc}0&\:1\\\:Q&\:0\end{array}\right)\) gives the required eigenvalues \(\:{\:\lambda\:}_{\:\text{1,2}}\) as,
Where “I” is the identity matrix and \(\:Q={\partial\:}_{{{\Psi\:}}_{R}}\left({\partial\:}_{\xi\:}W\right)\). Here, for different parametric values, the fixed points and eigenvalues are obtained numerically due to the transcendental nature of the dynamical system. For \(\:Q<0,\) the imaginary eigen values correspond to the centers and for \(\:Q>0\), real eigen values correspond to saddle points83.
The dynamical system in Eq. (26) has been analysed by phase space plots to study the wave solutions in Figs. (2–5) for relativistic KAAW in a dense plasma with the support of numerical computations. The physical parameters used for plotting phase portraits are taken from the precincts of white dwarfs and neutron stars. For studying relativistic and ultrarelativistic regimes, the number density ranges8 from \(\:1\times\:{10}^{32}\:{m}^{-3}\) to \(\:1\times\:{10}^{38}\:{m}^{-3}\). Here, we have taken the parameters \(\:{N}_{o}=1\times\:{10}^{36}\:{m}^{-3}\) and magnetic field \(\:{B}_{o}=1\times\:{10}^{9}\text{T}\:\) for the aforementioned dense stars. Mach number \(\:M=0.6\:\)is taken for the subsonic regime that means the phase velocity is less than the relativistic Fermi ion sound velocity is i.e. \(\:(\omega\:/k)<{C}_{sFR}\). The angle of propagation \(\:\theta\:\) is taken \(\:{75}^{^\circ\:}\) because these KAAW propagate at large angles with respect to ambient magnetic field\(\:{\:B}_{o}\).
Dynamical plots with \(\:{n}_{0}=1\times\:{10}^{36}\:{m}^{-3}\), \(\:{B}_{0}=1\times\:{10}^{9}\text{T}\) ,\(\:\:\theta\:={75}^{^\circ\:}\)and \(\:M=0.6\). (a) Phase portrait at \(\:{\text{T}}_{R}=0.1\) (b) The corresponding nonlinear solitary pulses (c) The corresponding nonlinear periodic KAAW (d) The Effective potential.
In Fig. 3, we have plotted the phase portraits, effective potential and the corresponding nonlinear structures for partially degenerate relativistic KAAW. The phase portrait of dynamical system in Fig. 3(a) gives nonlinear homoclinic orbit (NHO) for \(\:{\text{T}}_{R}=0.1\), which corresponds to the nonlinear solitary pulses in Fig. 3(b). The closed nonlinear periodic orbits (NPO) within the NHO corresponds to the nonlinear periodic kinetic Alfven acoustic wave (NPKAAW) in Fig. 3(c). Mostly conservative systems cause such trajectories. The outermost closed orbit, which encloses NHO, is the nonlinear super periodic orbit (NSPO). For these trajectories, we have numerically calculated the equilibrium points (fixed points) and the eigenvalues. Eigenvalues show that the fixed point \(\:{\text{s}}_{o}\)is the saddle point and the fixed points \(\:{\text{s}}_{1},\:{\text{s}}_{2}\) are the centers. Figure 3(d) is the effective potential of the aforementioned nonlinear structures.
In case of a fully degenerate relativistic KAAW, the phase portraits, effective potential, and the corresponding nonlinear structures have been plotted in Fig. 4. It is observed in Fig. 4(a) that for a NHO, there exists a saddle point \(\:{\text{s}}_{o}\left(\text{0,0}\right)\) at the origin at \(\:{\text{T}}_{R}=0\). In the phase portrait this corresponds to the solitary waves. Figure 4(b) shows homoclinic orbits corresponding solitary kinetic Alfven acoustic waves (SKAAW). We have numerically obtained the range of Mach number at \(\:{\text{T}}_{R}=0\) for which SKAAW are obtained, i.e. \(\:0.316\le\:M\le\:1.4\). The range of Mach number confirms the existence of SKAAW for both the sub and supersonic regimes of a fully degenerate relativistic plasma. NPO enclosing the fixed point \(\:{\text{s}}_{1\:}\)corresponds to the NPKAAW is shown in Fig. 4(c). The eigenvalues have shown that the equilibrium point \(\:{\text{s}}_{1}\) is the center and the effective potential of the nonlinear structures is shown in Fig. 4(d).
Dynamical plots with \(\:{n}_{0}=1\times\:{10}^{36}\:{m}^{-3}\), \(\:{B}_{0}=1\times\:{10}^{9}\text{T}\) ,\(\:\:\theta\:={75}^{^\circ\:}\)and \(\:M=0.6\). (a) Phase portrait at \(\:{\text{T}}_{R}=0\) (b) The corresponding solitary KAAW (c) The corresponding nonlinear periodic KAAW (d) The effective potential.
Sagdeev potential plots corresponding to solitary profiles at \(\:{\text{T}}_{R}=0\) (a) effect of number density for \(\:{B}_{o}=1\times\:{10}^{9}\text{T}\) ,\(\:\:\theta\:={75}^{^\circ\:}\)and \(\:M=0.6\) (b) effect of magnetic field for \(\:{n}_{o}=1\times\:{10}^{36}\:{m}^{-3}\) ,\(\:\:\theta\:={75}^{^\circ\:}\)and \(\:M=0.6\) (c) effect of Mach number for \(\:{n}_{o}=1\times\:{10}^{36}\:{m}^{-3}\), \(\:{B}_{o}=1\times\:{10}^{9}\text{T}\) and \(\:\theta\:={75}^{^\circ\:}\:\)(d) effect of angle of propagation for \(\:{n}_{o}=1\times\:{10}^{36}\:{m}^{-3}\), \(\:{B}_{o}=1\times\:{10}^{9}\text{T}\) and \(\:M=0.6\).
In Fig. 5(a), (c), and (d), the increase in \(\:{n}_{o},\:M\) and \(\:\theta\:\) decreases the width of the Sagdeev potential. This effective potential corresponds to the SKAAW, where both the amplitude and energy density decrease with increasing values of number density, Mach number, and angle of propagation. This occurs because these physical parameters enhance dispersive effects relative to the nonlinear effects. Furthermore, as the ambient magnetic field increases, the energy density and amplitude of SKAAW grow, as shown in Fig. 5(b). It is well known that \(\:{B}_{o}\) consistently contributes energy to KAAW, and a similar effect is observed for relativistic SKAAW.
It has been examined that above \(\:{\text{T}}_{R}=0\), the nonlinear solitary pulses are obtained corresponding to homoclinic orbits, and we do not get solitary structures, as shown in Fig. 3. A further increase in temperature gives rise to nonlinear periodic orbits, as shown in Fig. 6. For \(\:{\text{T}}_{R}=0.3\), the phase space trajectories show NPO, enclosing the fixed point \(\:{\text{s}}_{1}\),i.e. center, corresponding to NPKAAW (Fig. 6(a)). The corresponding travelling wave solution is shown in Fig. 6(b).
-
A.
Finite amplitude solitary KAAW.
In order to get an exact solution of our nonlinear KAAW in a relativistic degenerate plasma, we use small amplitude perturbation expansion method, i.e. \(\:{{\Psi\:}}_{R}<1\). So, Eq. (24) gets the form
Here, we can get a localized stationary wave solution that describes the formation and spread of a nonlinear structure. The analytical form of the well-recognized soliton solution is given as
Where \(\:\varDelta\:\) is the width and \(\:{{\Psi\:}}_{0R}\) is amplitude of solitary structure and is expressed as
We get analytic range of Mach number for the existance of finite amplitude SKAAW in a partially degenerate relativistic plasma as,
And in case of fully degenerate relativistic plasma the analytic range for SKAAW is given as,
At \(\:{n}_{0}=1\times\:{10}^{36}\:{m}^{-3}\), \(\:{B}_{0}=1\times\:{10}^{9}\text{T}\) (a) \(\:{{\Psi\:}}_{0R}\)vs. \(\:M\) gives compressive regime at \(\:{\text{T}}_{R}=0\) (b) \(\:{\text{T}}_{R}\)vs. \(\:M\) gives the existence region for SKAAW (c) 3D plot showing hump and dip for partial degeneracy (d) The rarefactive SKAAW for \(\:{\text{T}}_{R}=0.258\) .
Analytical ranges for nonlinear SKAAW, given in Eqs. (34–36), are in accordance with the linear dispersion relation in Eq. (9), that for \(\:\theta\:={90}^{o}\), \(\:{k}_{z}\to\:0\) and there will be no propagation of nonlinear structures in the medium. In Eq. (34), the upper limit depends on the angle of propagation, plasma beta, partial degeneracy and relativistic effects but the lower limit only depends on the angle of propagation. We have also shown these ranges of Mach number graphically to analyse the existence regions for SKAAW.
For the relativistic regime, we get both the compressive and the rarefactive finite amplitude SKAAW. In the case of fully degenerate plasma, Fig. 7(a) is plotted against the amplitude and Mach number, gives the.
compressive SKAAW only for the complete range of Mach number. The shaded area in Fig. 7(b) accounts for the existence region of SKAAW for the fully and partially degenerate plasmas. Here it is worth noticing that for every value of normalized ambient temperature \(\:{T}_{R}\) there is an upper limit of Mach number given by the boundary of shaded area. The increase in \(\:{T}_{R}\) decreases the upper limit of Mach number but the lower limit remains fixed for the specific values of density, magnetic field and angle of propagation (see Eq. (34)). The existence regions of compressive and rarefactive SKAAW are also distinguished in Fig. 7(b). The humps and dips are due to the partial degeneracy, as shown in Fig. 7(c). After \(\:{T}_{R}=0.258\), the range of Mach number changes (see Eq. (35)) and we only get rarefactive SKAAW, as shown in Fig. 7(d). In the partial degeneracy limits the rise in \(\:{T}_{R}\) decreases the depth and increases the amplitude of SKAAW, which means that normalized ambient temperature lowers the energy density and increases the dispersive effects in comparison to the nonlinear effects.
It is observed that an increase in the angle of propagation for SKAAW decreases the upper and the lower range of Mach number at fixed value of \(\:{T}_{R}\:\), and hence the existence region for SKAAW narrows down (see Fig. 8(a) and (b)). The increase in ambient magnetic field and density also decreases the upper range of Mach number (see Fig. 8(c) and (d)). The effect of varying density on the normalized ambient temperature \(\:{T}_{R}\left(=T/\mu\:\right)\) is obvious and clear in Fig. 8(d). Therefore, at constant angle of propagation the increase in density can give SKAAW for the greater values of \(\:{T}_{R}\) .
Dynamical plots with \(\:{n}_{0}=1\times\:{10}^{38}\:{m}^{-3}\), \(\:{B}_{0}=2\times\:{10}^{10}\text{T}\) ,\(\:\:\theta\:={85}^{^\circ\:}\)and \(\:M=0.6\). (a) Phase portrait at \(\:{\text{T}}_{R}=0\) (b) The corresponding nonlinear solitary pulses (c) The corresponding nonlinear periodic KAAW (d) The Effective potential.
Nonlinear analysis of Ultrarelativistic KAAW
In the ultrarelativistic regime (\(\:{\epsilon\:}_{Fo}\gg\:{m}_{e}{c}^{2}\)) of nonlinear KAAW, we follow the steps of Sagdeev pseudopotential approach of the previous sub section and derive an expression of Sagdeev potential by using Eq. (10),
Where \(\:{C}_{3}\), \(\:{C}_{4}\) and \(\:{\alpha\:}_{U}\) are the constants of integration and carry the ultrarelativistic effects with partial degeneracy, defined as
The Sagdeev Pseudopotential \(\:V\left({{\Psi\:}}_{U}\right)\) for the ultrarelativistic regime is different from the relativistic regime and does not contain fractional nonlinearities. Equation (37) does not give solitary KAAW for both fully and partially degenerate ultrarelativistic plasmas (see Fig. (8)). As discussed in the previous section, for relativistic SKAAW, an increase in number density enhances wave dispersion, leading to a decrease in energy density. In the ultrarelativistic regime, a further rise in number density intensifies dispersive effects, disrupting the balance between dispersion and nonlinearity. This imbalance causes the system to transition from SKAAW to nonlinear periodic KAAW.
To analyze the nonlinear KAAW for ultrarelativistic regime, we follow the steps of dynamical analysis of the previous sub section to understand the wave trajectories in a phase space. We have taken the parameters \(\:{N}_{o}=1\times\:{10}^{38}\:{m}^{-3}\), magnetic field \(\:{B}_{0}=2\times\:{10}^{10}\text{T}\:\), and angle of propagation is \(\:\theta\:={85}^{^\circ\:}\). Figure 9 comprises of the phase portraits, effective potential and the corresponding nonlinear structures. The phase portrait in Fig. 9(a) gives NHO for \(\:{\text{T}}_{R}=0\), corresponding to the nonlinear solitary pulse in Fig. 9(b). The NPO within the NHO corresponds to the nonlinear periodic kinetic Alfven acoustic waves in ultrarelativistic regime (see Fig. 9(c)). For these phase portraits, we have numerically calculated two fixed points, a center and a saddle point. Figure 9(d) is the effective potential of the nonlinear structures.
The exact solution of the aforementioned system can also be obtained by following the finite amplitude perturbation expansion method of the previous sub section to analyse the nonlinear solitary KAAW in ultrarelativistic regime.
Conclusions
We have investigated the linear and nonlinear analysis of relativistic and ultrarelativistic KAAW for the wave propagation in dense magnetized degenerate plasmas. The frequencies of fast and slow modes have shown variation with small temperature corrections. Wave dispersion for relativistic, ultrarelativistic and nonrelativistic KAAW have shown that these waves propagate more effectively in outer regions with lower specific heat than in denser inner regions of a neutron stars. Nonlinear analysis of arbitrary amplitude perturbations has shown SKAAW for the fully degenerate case and nonlinear solitary pulses or NPO in the relativistic regime for partial degeneracy. It has been observed that for the existence of nonlinear solitary waves, the range of Mach number is strongly coupled with the normalized ambient temperature \(\:{T}_{R}\). The range of Mach number has shown variation with magnetic field and density. The existence regions have been plotted for finite amplitude perturbations and it has been shown that the humps and dips in SKAAW are due to the small temperature corrections. Moreover, after a certain temperature, only rarefactive structures can propagate. It is due to the fact that the change in \(\:{\text{T}}_{R}\) shifts the balance between nonlinearity and dispersion, which changes the sign of amplitude \(\:\left(\:{{\Psi\:}}_{0R}\:\right)\)[see Eq. (33)]. It has been found that for the ultrarelativistic case rise in number density increases dispersive effects, disrupting the balance between dispersion and nonlinearity, and the system transit from SKAAW to nonlinear periodic KAAW. These linear and nonlinear wave propagation for the relativistic and ultrarelativistic KAAW have never been investigated before. Analytical and graphical results have shown that these waves can transfer energy by propagate through the regions of ultra-high densities and temperatures and may have potential applications in astrophysics and ultrahigh energy physics like relativistic plasma nanophononics23.
Data availability
The datasets used and/or analysed during the current study are available from the corresponding authoron reasonable request.
References
KILLIAN, T. C. Experiments in botany. Nature 441, 298–298 (2006).
De Marco, L. et al. A Fermi Degenerate Gas of Polar molecules. Science 363, 853–856 (2019).
Roth, J. R. Astrophysical implications of the continuity equation plasma oscillation. Nature 226, 626–628 (1970).
Flowers, E., Ruderman, M. & Sutherland, P. Neutrino pair emission from finite-temperature neutron superfluid and the cooling of young neutron stars. Astrophys. J. Vol 205 Apr 15 Pt 1 P 541–544 205, 541–544 (1976). (1976).
Lehnert, B. Electromagnetic Phenomena in Cosmical Physics (Cambridge University Press, 2016).
Blandford, R. D. The Phenomena of high energy astrophysics. in Symposium-International Astronomical Union vol. 214 3–20Cambridge University Press, (2003).
Gold, T. Rotating neutron stars as the origin of the pulsating radio sources. Nature 218, 731–732 (1968).
Koester, D. & Chanmugam, G. Physics of white dwarf stars. Rep. Prog Phys. 53, 837–915 (1990).
Chatterjee, S. & Cordes, J. M. Bow shocks from neutron stars: scaling laws and Hubble Space Telescope observations of the Guitar Nebula. Astrophys. J. 575, 407 (2002).
Guillot, S. et al. Hubble Space Telescope nondetection of PSR J2144–3933: the coldest known Neutron Star∗. Astrophys. J. 874, 175 (2019).
Gebhardt, K., Rich, R. M. & Ho, L. C. An intermediate-mass black hole in the globular cluster G1: improved significance from new keck and Hubble Space Telescope observations. Astrophys. J. 634, 1093 (2005).
Kormendy, J. et al. Hubble Space Telescope Spectroscopic evidence for a 1\times 109 M☉ black hole in NGC 4594. Astrophys. J. 473, L91 (1996).
Richer, H. B. et al. White dwarfs in globular clusters: Hubble Space Telescope observations of M4. Astrophys. J. 484, 741 (1997).
O’Brien, M. S., Bond, H. E. & Sion, E. M. Hubble Space Telescope spectroscopy of V471 Tauri: oversized K star, paradoxical white dwarf. Astrophys. J. 563, 971 (2001).
Kaltenegger, L. et al. The white dwarf opportunity: robust detections of molecules in earth-like exoplanet atmospheres with the James Webb Space Telescope. Astrophys. J. Lett. 901, L1 (2020).
Chatterjee, S. et al. Faint light of old neutron stars and detectability at the James Webb Space Telescope. Phys. Rev. D. 108, L021301 (2023).
Inayoshi, K., Onoue, M., Sugahara, Y., Inoue, A. K. & Ho, L. C. The age of Discovery with the James Webb Space Telescope: excavating the spectral signatures of the first massive black holes. Astrophys. J. Lett. 931, L25 (2022).
Woltjer, L., X-Rays, Type, I. S. & Remnants Astrophys. J. Vol 140 P 1309–1313 140, 1309–1313 (1964).
Patterson, J. Distinguishing between a white dwarf and a neutron star in an X-ray binary. Nature 292, 810–811 (1981).
Bochenek, C. D. et al. A fast radio burst associated with a Galactic magnetar. Nature 587, 59–62 (2020).
Chatterjee, G. et al. Magnetic turbulence in a table-top laser-plasma relevant to astrophysical scenarios. Nat. Commun. 8, 15970 (2017).
Nenstiel, C. et al. Electronic excitations stabilized by a degenerate electron gas in semiconductors. Commun. Phys. 1, 38 (2018).
Purvis, M. A. et al. Relativistic plasma nanophotonics for ultrahigh energy density physics. Nat. Photonics. 7, 796–800 (2013).
Ji, L. L., Snyder, J., Pukhov, A., Freeman, R. R. & Akli, K. U. towards manipulating relativistic laser pulses with micro-tube plasma lenses. Sci. Rep. 6, 23256 (2016).
Manfredi, G. How to model quantum plasmas. Preprint at (2005). http://arxiv.org/abs/quant-ph/0505004
Pukhov, A. Three-dimensional electromagnetic relativistic particle-in-cell code VLPL (virtual laser plasma lab). J. Plasma Phys. 61, 425–433 (1999).
Landau, L. D. & Lifshitz, E. M. Statistical Physics: Volume 5vol. 5 (Elsevier, 2013).
Balberg, S. & Shapiro, S. L. The Properties of Matter in White Dwarfs and Neutron Stars. Preprint at (2000). http://arxiv.org/abs/astro-ph/0004317
Seemann, O., Wan, Y., Tata, S. & Kroupp, E. Malka, V. Refractive plasma optics for relativistic laser beams. Nat. Commun. 14, 3296 (2023).
Mourou, G. A., Tajima, T. & Bulanov, S. V. Optics in the relativistic regime. Rev. Mod. Phys. 78, 309–371 (2006).
Modena, A. et al. Electron acceleration from the breaking of relativistic plasma waves. Nature 377, 606–608 (1995).
Joshi, C. et al. Ultrahigh gradient particle acceleration by intense laser-driven plasma density waves. Nature 311, 525–529 (1984).
Levy, E. H. & ROSE, W. K. Origin of neutron star magnetic fields. Nature 250, 40–41 (1974).
Joss, P. C. Stars at the quantum limit. Nature 360, 15–15 (1992).
Chabrier, G., Ashcroft, N. W. & DeWitt, H. E. White dwarfs as quantum crystals. Nature 360, 48–50 (1992).
Leung, Y. C. & Wang, C. G. Minimum Mass of a Neutron Star. Nat. Phys. Sci. 240, 132–133 (1972).
Lenzen, R. & Trümper, J. Reflection of X rays by neutron star surfaces. Nature 271, 216–220 (1978).
Pacini, F. Energy emission from a neutron star. Nature 216, 567–568 (1967).
Fleishman, G. D. & Toptygin, I. N. Cosmic Electrodynamics: Electrodynamics and Magnetic Hydrodynamics of Cosmic Plasmasvol. 388 (Springer New York, 2013).
Tsintsadze, N. L. & Tsikarishvili, E. G. Parametric instabilities in relativistic plasma. Astrophys. Space Sci. 39, 191–199 (1976).
Flacco, A. et al. Persistence of magnetic field driven by relativistic electrons in a plasma. Nat. Phys. 11, 409–413 (2015).
Tsintsadze, N. L., Rasheed, A., Shah, H. A. & Murtaza, G. Nonlinear screening effect in an ultrarelativistic degenerate electron-positron gas. Phys. Plasmas. 16, 112307 (2009).
Renninger, W. H. & Wise, F. W. Optical solitons in graded-index multimode fibres. Nat. Commun. 4, 1719 (2013).
Blaschke, F., Karpíšek, O. N. & Beneš, P. Solitons in the Peyrard–Bishop model of DNA and the Renormalization Group method. Prog. Theor. Exp. Phys. 063J02 (2020). (2020).
Georgiev, D. D. & Glazebrook, J. F. On the quantum dynamics of Davydov solitons in protein α-helices. Phys. Stat. Mech. Its Appl. 517, 257–269 (2019).
Yefsah, T. et al. Heavy solitons in a fermionic superfluid. Nature 499, 426–430 (2013).
Roy, S. & Misra, A. P. Stability and evolution of electromagnetic solitons in relativistic degenerate laser plasmas. J. Plasma Phys. 86, 905860611 (2020).
Roy, S., Misra, A. P. & Abdikian, A. Modulation of electromagnetic waves in a relativistic degenerate plasma at finite temperature. Phys. Fluids. 35, 066123 (2023).
El-Labany, S. K., El-Taibany, W. F., Behery, E. E. & Abd-Elbaki, R. Oblique collision of ion acoustic solitons in a relativistic degenerate plasma. Sci. Rep. 10, 16152 (2020).
Gurevich, A. V. Distribution of captured particles in a potential well in the absence of collisions. Sov Phys. JETP. 26, 575–580 (1968).
Pitaevskii, L. P. & Lifshitz, E. M. Physical Kinetics: Volume 10vol. 10 (Butterworth-Heinemann, 2012).
Erokhin, N. S., Zol’Nikova, N. N. & Mikhailovskaya, L. A. Asymptotic theory of the nonlinear interaction between a whistler and trapped electrons in a nonuniform magnetic field. Plasma Phys. Rep. 22, (1996).
Matsumoto, H. & Omura, Y. Computer simulation studies of VLF triggered emissions deformation of distribution function by trapping and detrapping. Geophys. Res. Lett. 10, 607–610 (1983).
Hernandez, J. P. Electron self-trapping in liquids and dense gases. Rev. Mod. Phys. 63, 675–697 (1991).
Romano, F. et al. Design and First Tests of the Trapped Electrons Experiment T-REX. Preprint at (2024). http://arxiv.org/abs/2406.19123
Everett, M. et al. Trapped electron acceleration by a laser-driven relativistic plasma wave. Nature 368, 527–529 (1994).
Shah, H. A., Masood, W., Qureshi, M. N. S. & Tsintsadze, N. L. Effects of trapping and finite temperature in a relativistic degenerate plasma. Phys. Plasmas. 18, 102306 (2011).
Shah, H. A., Iqbal, M. J., Tsintsadze, N., Masood, W. & Qureshi, M. N. S. Effect of trapping in a degenerate plasma in the presence of a quantizing magnetic field. Phys. Plasmas. 19, 092304 (2012).
Hasegawa, A. & Mima, K. Exact Solitary Alfvén Wave. Phys. Rev. Lett. 37, 690–693 (1976).
Yu, M. Y. & Shukla, P. K. Finite amplitude solitary Alfvén waves. Phys. Fluids. 21, 1457–1458 (1978).
Hollweg, J. V. Kinetic Alfvén wave revisited. J. Geophys. Res. Space Phys. 104, 14811–14819 (1999).
Burke, A. T., Maggs, J. E. & Morales, G. J. Spontaneous fluctuations of a temperature filament in a magnetized plasma. Phys. Rev. Lett. 84, 1451–1454 (2000).
Maggs, J. E. & Morales, G. J. Magnetic fluctuations associated with field-aligned striations. Geophys. Res. Lett. 23, 633–636 (1996).
Pu, Z. & Zhou, Y. The kinetic Alfvén wave instability driven by a sheared plasma flow and the associated anomalous transport. Sci. Sin Ser. Math. Phys. Tech. Sci. 29, 301–311 (1986).
Chmyrev, V. M. et al. Alfvén vortices and related phenomena in the ionosphere and the magnetosphere. Phys. Scr. 38, 841–854 (1988).
Sadiq, N. & Ahmad, M. Kinetic Alfven waves in dense quantum plasmas with effect of spin magnetization. Plasma Res. Express. 1, 025007 (2019).
Hussain, A., Iqbal, Z., Brodin, G. & Murtaza, G. On the kinetic Alfvén waves in nonrelativistic spin quantum plasmas. Phys. Lett. A. 377, 2131–2135 (2013).
Sabeen, A., Shah, H. A., Masood, W. & Qureshi, M. N. Finite amplitude solitary structures of coupled kinetic Alfven-acoustic waves in dense plasmas. Astrophys. Space Sci. 355, 225–232 (2015).
Suzuki, T. K. & Nagataki, S. Alfven Wave–driven proto–neutron star winds and r -Process nucleosynthesis. Astrophys. J. 628, 914–922 (2005).
Yuan, Y. et al. Magnetar bursts due to Alfvén Wave Nonlinear Breakout. Astrophys. J. 933, 174 (2022).
Van Hoven, M. & Levin, Y. Hydromagnetic waves in a superfluid neutron star with strong vortex pinning. Mon Not R Astron. Soc. 391, 283–289 (2008).
Yuan, Y., Beloborodov, A. M., Chen, A. Y. & Levin, Y. Plasmoid ejection by Alfvén Waves and the fast radio bursts from SGR 1935 + 2154. Astrophys. J. Lett. 900, L21 (2020).
Thompson, C. & Duncan, R. C. The soft gamma repeaters as very strongly magnetized neutron stars. II. Quiescent neutrino, X-ray, and Alfven wave emission. Astrophys. J. 473, 322 (1996).
Chen, A. Y., Yuan, Y., Beloborodov, A. M. & Li, X. Relativistic Alfvén waves entering charge-starvation in the magnetospheres of Neutron stars. Astrophys. J. 929, 31 (2022).
Cramer, N. F. The Physics of Alfvén Waves (Wiley, 2011).
Wu, D. J. & Chen, L. Kinetic Alfvén Waves in Laboratory, Space, and Astrophysical Plasmas (Springer Singapore, 2020). https://doi.org/10.1007/978-981-13-7989-5
Asam, M. T. et al. Effect of trapping in coupled kinetic Alfven-acoustic waves in a partially degenerate plasma with quantizing magnetic field. Phys. Scr. 98, 025612 (2023).
Wu, D. J., Wang, D. Y. & Fälthammar, C. G. An analytical solution of finite-amplitude solitary kinetic Alfvén waves. Phys. Plasmas. 2, 4476–4481 (1995).
Chatterjee, P., Roy, K., Muniandy, S. V., Yap, S. L. & Wong, C. S. Effect of ion temperature on arbitrary amplitude ion acoustic solitary waves in quantum electron-ion plasmas. Phys. Plasmas 16, (2009).
Chow, S. N. & Hale, J. K. Methods of Bifurcation Theoryvol. 251 (Springer Science & Business Media, 2012).
Guckenheimer, J. & Holmes, P. Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fieldsvol. 42 (Springer Science & Business Media, 2013).
Yadav, K., Prasad, A. & Shrimali, M. D. Control of coexisting attractors via temporal feedback. Phys. Lett. A. 382, 2127–2132 (2018).
Perko, L. Differential Equations and Dynamical Systemsvol. 7 (Springer Science & Business Media, 2013).
Acknowledgements
The Pakistan Academy of Sciences provided a stipend to one of us M.T.A. for the period when this work was carried out.
Author information
Authors and Affiliations
Contributions
“M.T.A. carried out the calculations and the numerical work along with preparation of the figures""H.A.S. suggested the problem and checked the calculations and the numerical work and the writing up of the manuscript""W.M. checked and helped with calculations, numerical work and contributed to the initial draft of the manuscript”.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Appendix: nonlinear analysis
Appendix: nonlinear analysis
We obtain Sagdeev potential (pseudopotential) by using the dimensionless frame \(\:\xi\:\:={\:k}_{x}\:x+{\:k}_{z}\:z-M\:t\) in the system of Eqs. (1–4) and obtain Eqs. (13–16). In deriving Eq. (19), we plug in Eq. (17) in Eq. (13) and get,
Using Eqs. (A.3) in (A.1), we get
In the limit \(\:{{\Psi\:}}_{R}<1\), the above expression takes the form,
The first term in (A.5) is simplified as follows,
And the rest of the terms have undergone the same procedure. Eq. (A.4) finally takes the form,
Integrating equation (A.8) and applying the boundary conditions \(\:{{v}_{z},\:\psi\:}_{R}\to\:0\) as \(\:\xi\:\to\:\infty\:\), we obtain the parallel ion velocity \(\:{v}_{z}\)
Here, the constant of integration \(\:{C}_{1}\) is,
Equation (A.9) corresponds to Eq. (20) in the manuscript where we have modified the terms with the powers of
Now solving Eqs. (14) and (17) in the manuscript, we get
Solving (A.11) and Eq. (18) of manuscript, we get
We insert Eqs. (A.3) and (A.10) in the above expression and expand the term \(\:\raisebox{1ex}{${\alpha\:}_{R}$}\!\left/\:\!\raisebox{-1ex}{${n}_{R}$}\right.\) with respect to temperature. Here, we have ignored higher order terms and the terms containing \(\:{T}_{R}^{4}\) (for a quantum plasma \(\:{\text{T}}_{R}<1\)). Moreover, equation (A.12) is integrated to express through Sagdeev potential in the usual manner,
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Asam, M.T., Shah, H.A. & Masood, W. Kinetic alfven-acoustic waves at relativistic and ultra-relativistic Fermi energies. Sci Rep 15, 6853 (2025). https://doi.org/10.1038/s41598-024-82416-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-82416-5