Nonlinear modulation of dispersive fast magnetosonic waves in an inhomogeneous rotating solar low- magnetoplasma
Abstract
We study the modulation of fast magnetosonic waves (MSWs) in rotating inhomogeneous low- magnetoplasmas with the effects of gravitation and the Coriolis force. By employing the standard multiple-scale reductive perturbation technique (RPT), we derive a nonlinear Schrödinger (NLS) equation that governs the evolution of slowly varying MSW envelopes. The fast MSW becomes dispersive by the effects of the Coriolis force in the fluid motion, and the magnetic field and density inhomogeneity effects favor the Jeans instability in self-gravitating plasmas in a larger domain of the wave number (, below the Jeans critical wave number, ) than homogeneous plasmas. The relative influence of the Jeans frequency (, associated with the gravitational force) and the angular frequency (, relating to the Coriolis force) on the Jeans carrier MSW mode and the modulational instability (MI) of the MSW envelope is studied. We show that the MSW envelope (corresponding to the unstable carrier Jeans mode with and ) is always unstable against the plane wave perturbation with no cut-offs for growth rates. In contrast, the stable Jeans mode with but manifests either modulational stability or MI having a finite growth rate before being cut off. We find an enhancement of the MI growth rate by the influence of magnetic field or density inhomogeneity. The case with constant gravity force (other than the self-gravity) perpendicular to the magnetic field is also briefly discussed to show that the fast magnetosonic carrier mode is always unstable, giving MI of slowly varying envelopes with no cut-offs for the growth rates. Possible applications of MI in solar plasmas, such as those in the X-ray corona, are also briefly discussed.
I Introduction
One-fluid magnetohydrodynamics deals with a compressible conducting fluid immersed in a magnetic field and it is often regarded as a reasonable description of the large-scale dynamics of a plasma. Magnetohydrodynamic (MHD) waves are generally described applying this theory and represent one of the macroscopic processes responsible for the transformation of energy and information in plasmas. The theory of MHD waves in an infinite conducting medium was first developed by Alfvén [1] with its application to sunspots, coronal heating, particle acceleration, and generation of cosmic radiation. Since then, the study of nonlinear MHD waves in plasma has been one of the most popular research topics among researchers, given its remarkable application and progress in laboratory experiments and techniques as well as various space and astrophysical plasma environments like pulsar magnetosphere, magnetars, solar corona, etc.
Magnetosonic waves (MSWs) are one of the fundamental MHD wave modes in plasmas, which propagate nearly perpendicular to the background magnetic field and are often observed in laboratory plasmas [2], earth magnetosphere [3, 4, 5, 6, 7], and solar wind plasmas [8], etc. These waves are of great interest because of their important roles in plasma heating [9] and charged particle acceleration [10]. Several authors have investigated the characteristics of MSW waves in different space and astrophysical plasma environments. To mention a few, Marklund et al. [11] studied the magnetosonic solitons in a quantum magnetoplasma, including the quantum Bohm potential and electron spin-1/2 effects using the Sagdeev potential approach. Haas and Mahmood [12] analyzed the propagation of linear and weakly nonlinear magnetosonic waves in a plasma with arbitrary degeneracy of electrons and with the inclusion of Bohm diffraction effects. Nonlinear properties of fast MSWs in dense dissipative plasmas with degenerate electrons were studied theoretically by Masood et al. [13] by deriving the Zabolotskaya-Khokhlov (ZK) equation for small but finite amplitude excitations. Hussain and Mahmood [14] investigated the nonlinear propagation of MSWs and showed that these waves may evolve into shock-like structures that may be responsible for heating the solar chromosphere and the solar corona.
The modulational instability (MI) of nonlinear waves in dispersive or diffractive media has been known to be one of the most important mechanisms of energy localization via the formation of different localized coherent structures like envelope solitons [15], envelope shocks [16], freak wave (or rouge waves) [17, 18], giant waves [19], etc., as well as the transfer of energy between waves and particles, leading to particle heating, e.g., heating of coronal loops [20]. In this context, Watanabe [21] first experimentally observed the MI of nonlinear wave envelopes in dispersive media in . Subsequently, the investigations of MI and associated nonlinear structures have gained significant attention among researchers, and a large number of theoretical and experimental investigations on MI of electrostatic and electromagnetic waves have been done to explain the effects of different physical parameters in various space and astrophysical plasmas [22, 23, 24, 25]. For example, Sahyouni et al. [26] studied the amplitude modulation of fast magnetosonic surface waves in solar flux tubes and showed that the fast wave-envelope admits dark envelope solitons solution and discussed the possibility of the existence of solitary waves in the solar atmosphere. Sakai [27] examined the MI of fast MSWs theoretically and discussed its applications to solar plasmas. Misra and Shukla [28] studied the MI of magnetosonic waves using a two-fluid quantum magnetohydrodynamic model that includes the effects of the electron-1/2 spin and the plasma resistivity. Panwar et al. have investigated the MI and associated rogue-wave structures of slow magnetosonic perturbations in a Hall-MHD plasma [29]. Wang et al. have studied the MI of MSWs in a Fermi-Dirac-Pauli plasma by the combined effects of the electron relativistic degeneracy, the quantum tunneling, electron spin via Pauli paramagnetism, and plasma resistivity [30].
To the best of our knowledge, the combined effects of the Coriolis force and the gravity force on the fast magnetosonic modes and the MI of slowly varying envelopes in inhomogeneous magnetoplasmas have not been studied before. Our aim is to study the propagation characteristics of fast magnetosonic modes and their nonlinear evolution as slowly varying envelopes through the modulational instability. To this end, we consider a self-consistent MHD model for electron and ion fluids that includes the combined effects of the gravity and Coriolis forces as well as the magnetic field and density inhomogeneities in absence of any viscosity or magnetic diffusivity. We show that while the fast magnetosonic modes exhibit instability in plasmas with constant gravity perpendicular to the magnetic field, the fast magnetosonic Jeans mode in self-gravitating plasmas can be stable or unstable depending on the competitive roles of the gravity and Coriolis forces and whether the Jeans wave number is below or above a critical value. The possibility of the emergence of MI of slowly varying magnetosonic envelopes in different cases of the stable or unstable fast carrier modes in presence of the self-gravitation and the constant gravity force with the effects of magnetic field or density inhomogeneities are studied and its applications to solar coronal plasmas are discussed.
The paper is organised in the following fashion: In Sec. II, we present the basic MHD models governing the dynamics of fast magnetosonic waves and the inhomogeneous equilibrium state of plasmas. Section III demonstrate the linear fast magnetosonic mode, the compatibility condition, and the derivation of the nonlinear Schrödinger (NLS) equation for the evolution of slowly varying magnetosonic envelopes. We study the MI in Sec. IV for three different cases of stable and unstable Jeans carrier modes as well as for the unstable carrier mode under the influence of gravity. The applications of our results in solar plasmas are discussed in Sec. V. Finally, Sec. VI is left to summarize and conclude our results.
II Basic equations
We consider the nonlinear propagation of fast magnetosonic waves in a rotating magnetized plasma with the effects of the Coriolis force and the gravitational force. The latter may be considered in two cases (i) when the gravitational acceleration is not a constant (self-gravity) and (ii) when is a constant. The plasma is supposed to be rotating with uniform angular velocity and immersed in an external static magnetic field along the -axis, i.e.. We assume the gravitational force to be acting vertically downwards and parpendicular to the magnetic field and the wave propagation along the -axis, for simplicity. Furthermore, the background magnetic field and the density are assumed to vary along the -axis (inhomogeneities). A schematic diagram for the model is shown in Fig. 1. We first consider the case of self-gravitating plasmas. The basic equations describing the dynamics of electron-ion fluids in self-gravitating magnetoplasmas in the center-of-mass frame are [31]
(1) |
(2) |
(3) |
(4) |
where is the gravitational force per unit mass of the fluid acting vertically downwards; , , and are, respectively, the fluid density, fluid velocity, and the thermal pressure. Also, , , , and , respectively, denote the magnetic field, uniform angular velocity of the rotating fluid, gravitational potential, and the universal gravitational constant. The pressure satisfies the equation of state: , where is the adiabatic sound speed with denoting the Boltzmann constant, the electron temperature, the adiabatic index, and the ion mass.
Next, we consider the one-dimensional propagation of fast MSWs along the -axis and normalize Eqs. (1)-(4) according to: , , , , , , , . Here, is the Alfvén speed, is the squared Jeans frequency, and is the ion-cyclotron frequency at . Furthermore, the space and time coordinates are normalized as and . Thus, from Eqs. (1)-(4), after separating the velocity components along the axes and noting that the frozen-in-field condition holds for Eqs. (1) and (3), we obtain the following reduced set of normalized equations
(5) |
(6) |
(7) |
(8) |
(9) |
II.1 Equilibrium state
At equilibrium, the background plasma state with density and magnetic field inhomogeneities is defined by the following pressure-balance and the Poisson equations
(10) |
(11) |
where the suffix ‘’ denotes the equilibrium value of the corresponding physical quantity. We define , and . By means of the normalization defined before, Eqs. (10) and (11) reduce to
(12) |
(13) |
Using the relation , and noting that (Frozen-in-field condition), we obtain from Eq. (12), the following relation (in normalized form)
(14) |
where . Typically, if is not too small, i.e., if it lies in , the function approaches from small negative values to zero as . In addition, if the magnetic force dominates over the pressure gradient force (in low- plasmas such as those in solar corona) or if the length scale of magnetic field inhomogeneity is much smaller than the density inhomogeneity, i.e., , where and , then the term proportional to in Eq. (14) can be neglected compared to the the term involving the factor (the second term on the left-hand side). Thus, from Eqs. (13) and (14), we obtain for the following approximate solutions for and .
(15) |
On the other hand, if the contribution from the pressure gradient force is much higher than the magnetic force (in high- plasmas such as those in the solar photosphere and solar wind acceleration region) or if the length scale of magnetic field inhomogeneity is much larger than the density inhomogeneity, i.e., , approximate solutions for and can then be obtained from Eqs. (13) and (14) as
(16) |
where . Furthermore, for plasmas where the plasma pressure is comparable to the magnetic pressure , such as those in the lower chromospheric region of the solar atmosphere, Eqs. (13) and (14) are to be solved numerically. However, we are not considering these two cases in the present investigation.
III Derivation of NLS equation
We study the modulation of weakly nonlinear slowly varying magnetosonic wave envelopes that are generated due to nonlinear self-interactions of carrier fast magnetosonic modes and higher harmonic modes in self-gravitating rotating magnetoplasmas. To this end, we employ the standard multiple-scale reductive perturbation technique (RPT) [32] to Eqs. (5)-(9) and derive the nonlinear Schrödinger equation for the evolution of slowly varying MSW envelopes. In the RPT, we define a new frame of reference in which the space and time variables are stretched as
(17) |
where is the group velocity of the wave envelope along the -axis and is a small () expansion parameter, which scales the weakness of amplitudes of perturbations. Due to the stretched coordinates (17), the space and time derivatives will be replaced according to the following transformations.
(18) | |||
The dynamical variables are divided into unperturbed (equilibrium) and perturbed parts. In the latter, the slow and fast scales (for space and time), respectively, enter the -th harmonic amplitudes and the phase . Thus, the variables can be expanded as
(19) |
where and are normalized by and respectively, and the slowly varying wave amplitudes , , , , , etc. satisfy the reality condition: , where the asterisk denotes the complex conjugate, and and are, respectively, the wave number and the wave frequency of fast carrier MSWs.
In what follows, we apply the transformations [Eqs. (17) and (III)] and substitute the expansions [Eq. (19)] into the normalized Eqs. (5)-(9), and then obtain equations for different harmonic modes corresponding to different powers of (For some details, see Appendix n-th order reduced equations). Here, we assume that the length scales of the density () and magnetic field () inhomogeneities are much larger than the length scale of magnetosonic perturbations, i.e., . With these assumptions, the imaginary contributions originating from the inhomogeneities can be ignored. The results are given in Secs. III.1-III.3.
III.1 First-order perturbations: Linear dispersion relation
For the first-order first harmonic perturbations, we obtain from the coefficients of , the following relations
(20) |
(21) |
where is the normalized inhomogeneous Alfvén velocity.
Next, eliminating the variables and looking for their nonzero solutions, we obtain from Eqs. (20)-(21), the following linear dispersion relation for the fast carrier magnetosonic modes in self-gravitating magnetoplasmas.
(22) |
where is the critical Jeans wave number modified by the Coriolis force, given by,
(23) |
From Eq. (22), we note that the fast magnetosonic wave becomes dispersive due to the presence of the term proportional to , associated with the Coriolis force, and the Jeans instability may occur due to the term proportional to by the influence of the self-gravitating force [33]. The instability occurs in the region , provided the self-gravity force dominates over the Coriolis force, i.e., . In the other region, i.e., , the magnetosonic wave can propagate as a real (stable) eigenmode with . On the other hand, in absence of the gravity effects or when the Coriolis force dominates over the self-gravity force with , the MSWs can also propagate as a real mode (without any instability) with the frequency being smaller or larger than the ion-cyclotron frequency and the dispersion relation in the form of high-frequency Langmuir waves in classical plasmas. Thus, it is reasonable to investigate the nonlinear modulation of slowly varying magnetosonic fields by means of a NLS equation. We also note that the effects of the magnetic field and the density inhomogeneities enter the coefficient of via the Alfvén velocity , implying that the wave dispersion is greatly modified by the effects of inhomogeneities, and hence the modifications of the phase velocity as well as the group velocity dispersion of magnetosonic envelopes.
Figure 2 displays the profiles of the real wave mode (, for ) and the growth rate of instability (, for ) by the effects of the density and magnetic field inhomogeneities [See Eq. (15)], the Coriolis force, and the self-gravity force such that . We find that while the real wave frequency increases with the wave number , the instability growth rate falls off from a nonzero value with having a cut-off at the critical Jeans wave number, i.e., . Such a critical value shifts towards a higher value of due to the effects of the inhomogeneities in which increases, but decreases with increasing values of . In this case, the instability domain for expands and the domain of the real wave mode reduces. Thus, it follows that the density and magnetic field inhomogenities in the background plasma favor the Jeans instability in a wide range of values of , not reported before. On the other hand, a small reduction of the angular frequency of the rotating fluid, associated with the Coriolis force, can significantly reduce both the growth rate of instability and the instability domain in . We do not consider the case of at which the fast magnetosonic wave becomes dispersionless. Such a strict condition may be applicable to low-frequency long-wavelength plasma oscillations.
III.2 Second-order perturbations: Compatibility condition
For the second order () reduced equations with , we obtain the following expressions for different harmonic modes
(24) |
together with the compatibility condition
(25) |
where is the phase velocity of the carrier magnetosonic waves. From Eq. (25), it is evident that the group velocity of the wave envelope has an inverse relationship with the phase velocity of the carrier wave, i.e., the magnitude of the group velocity increases with a reduction of the magnitude of the phase velocity. Because of the dependency of on the wave frequency , the group velocity can be real (imaginary) for when . It can also be real when for any real value of . The expression for is clearly modified by the effects of the background magnetic field or density inhomogeneity. Furthermore, having known the characteristics of (See Fig. 2), one can also investigate the features of the real and imaginary parts of . It is seen that the real part of the group velocity decreases with and it can be further reduced (increased) by the effects of the density or magnetic field inhomogeneity (Coriolis force). On the other hand, the imaginary part of is always negative in the domain and its magnitude increases with increasing values of . The detailed analysis of is not necessary at this stage, we will rather focus on the coefficients of the NLS equation to be derived in Sec. III.3.
III.3 Third-order perturbations: The NLS equation
The second-order harmonic modes for and appear due to the nonlinear self-interaction of the carrier waves, and they are found to be proportional to as
(26) |
The nonlinear self-interactions of the first-order carrier wave modes also result into the generation of the zeroth harmonic modes. Thus, for , we obtain
(27) |
Finally, substituting all the above derived expressions into the third-order harmonic modes ( and ), we obtain the following NLS equation
(28) |
where we have replaced by for simplicity. The group velocity dispersion coefficient and the nonlinear coefficient , appeared due to the carrier wave self-interactions, are given by
(29) |
(30) |
From the dispersion Eq. (22), we note that the wave frequency can be either real or purely imaginary depending on whether or . So, the coefficients and can also be either real or purely imaginary. While the coefficient appears due to the group velocity dispersion and the coupling of the thermal and magnetic pressures with the Coriolis and the Gravity forces, the coefficient appears due to the nonlinear self-interactions of higher harmonic modes and coupling between the zeroth- and second-order second harmonic carrier waves. Specifically, the first term of appears due to transverse velocity perturbations and the second term of is due to the coupling between the zeroth and second harmonic modes. Before proceeding to study the modulational instability and the evolution of magnetosonic wave envelopes, it is pertinent to mention there cases of interest:
-
•
Case I: , .
-
•
Case II: , .
-
•
Case III: , .
From the linear analysis in Sec. III.1, we have noted that while Cases I and II correspond to the Jeans stable mode, Case III corresponds to the Jeans instability. Thus, it is of interest to study the modulational instability conditions in these three cases. Specifically, we will examine whether the Jeans instability region of carrier waves can give rise to the modulational instability of magnetosonic envelopes.
IV MODULATIONAL INSTABILITY
We follow a similar technique as in [34] to study the modulation of magnetosonic wave envelopes against a plane wave perturbation. Though the analysis is standard, we reproduce it here for the sake of clarity to the readers. The NLS Eq. (28) admits a plane wave time-dependent solution of the form
(31) |
where is some constant, and and are real functions of and . Substituting Eq. (31) into Eq. (28), and separating the real and imaginary parts, we obtain
(32) |
(33) |
Next, we modulate the wave amplitude and phase by small plane-wave perturbations with the wave number and the wave frequency as
(34) |
where is a constant, and and for , are the amplitudes of perturbations.
Substituting the perturbation expansion (34) in Eqs. (32) and (33), and looking for nonzero solutions of and , we obtain the following dispersion relation for the perturbed wave of modulation in self-gravitating rotating magnetoplasmas
(35) |
In Subsections IV.1-IV.3, we will study the conditions for the modulational instability and the instability growth rate in the three different cases as mentioned above.
IV.1 Case I: ,
We consider the case when the contribution from the Coriolis force dominates over the self-gravitating force, i.e. when . In this case, the carrier Jeans wave frequency becomes real for any values of the wave number [See Eq. (22)], and so are and . Also, from Eq. (29), it is evident that that for . Furthermore, inspecting on the expression of [Eq. (30)], we find that when holds, the second term of becomes positive and larger than the first term, giving . Thus, in this case, , and Eq. (35) gives a real wave frequency for the perturbation of modulation, implying that the slowly varying magnetosonic wave envelope is always stable under the modulation. It follows that the modulated magnetosonic wave is also stable in absence of the self-gravity force.
IV.2 Case II: ,
We consider the case when the self-gravitating force dominates over the Coriolis force, i.e. and . In this case, as discussed before in Sec. III.1, since the carrier wave frequency is real, the coefficients and of the NLS equation are also real. However, in contrast to Case I, and can be both positive and negative for , so that one can have both stable and unstable regions. Figure 3 shows the plot of versus the carrier wave number . We find that due to the inhomogenity effects, the domain of in which (for which and the modulational instability occurs) expands and the domain shifts towards lower values of due to a reduction of the magnitude of . For example, for (i) , we have in (See the solid line of Fig. 3) (ii) , we have in (See the dashed line of Fig. 3), and (iii) , we have in (See the dotted line of Fig. 3). The other fixed parameter values are and .
In the case of for which the modulational instability occurs, the instability growth rate can be obtained from Eq. (35) for as
(36) |
where is the critical wave number of perturbation, and the maximum growth rate, is attained at .
Figure 4 displays the growth rate of instability [Eq. (36)] for different parameter values as in Fig. 2. We find that the instability growth rate can be enhanced and maximized by the effects of the magnetic field or density inhomogeneity. However, it can be reduced or minimized by reducing the angular frequency compared to the Jeans frequency.
IV.3 Case III: ,
Here, we consider the most interesting case when , but . From the analysis in Sec. III.1, it is evident that the carrier wave frequency is purely imaginary (Jeans instability) and so are the coefficients and of the NLS equation. Thus, assuming and , Eq. (35) gives
(37) |
where and are reals, given by,
(38) |
with
(39) |
From Eq. (37), we find that the instability condition depends on the sign of the product and/or on the value of smaller or larger than a critical value . When , Eq. (37) gives purely imaginary for all values of and so we have the modulational instability for all , but with . The corresponding growth rate of instability can be obtained by setting as
(40) |
On the other hand, when , the modulational instability occurs for , where is the critical wave number, given by,
(41) |
The corresponding instability growth rate is given by
(42) |
However, the modulated wave is stable for . The instability condition stated above is in contrast to the typical condition (, ) of modulational instability of wave envelopes in plasmas without self-gravity effects. Inspecting on the coefficients and , we find that is always positive for and . From Fig. 5, we find that can be both positive and negative in a finite domain of . The domain for in which the modulational instability occurs expands due to the effects of magnetic field or density inhomogeneity and shrinks due to reduction of the contribution from the Coriolis force compared to the self-gravitating force. For example, for (i) , we have in (See the solid line of Fig. 5) (ii) , we have in (See the dashed line of Fig. 5), and (iii) , we have in (See the dotted line of Fig. 5). The other fixed parameter values are and .
V Applications and discussions
The solar corona is typically a complex system in the self-gravitating field of the Sun, and has been known to be an active medium for the dynamics and stability of magnetosonic waves, since one of the most important problems of coronal heating is associated with the wave [35]. Recent Solar Optical Telescope (SOT) observations of small-scale oscillations also indicate the existence of fast magnetosonic waves in solar prominence threads and pillars [36]. Past observations in the X-ray corona also reported slowly moving perturbations that can propagate with the speed km/s during the first nine minutes of filament disruption [37]. Such a velocity seems to be close to the Alfvén speed of fast magnetosonic waves in a low- plasma. The speed was, however, reduced to km/s and km/s after twenty minutes and four hours respectively. To explain this observation, one can assume that the fast magnetosonic waves originating from disrupted filaments can propagate as slowly moving magnetosonic wave envelopes through the modulational instability [27].
Typical solar plasma parameters are [35] (i) the number density of electrons/ions, , (ii) the electron temperature, K, and (iii) the plasma at the magnetic field strength, G ( G). For illustration purpose, we calculate the maximum growth rate in the case of , (Case II in Sec. IV.2) at G as (for ) such that the typical growth time can be estimated as s, i.e., faster than predicted before [27] and may be reasonable to explain the transition time (twenty minutes) from fast magnetosonic mode to slowly varying envelopes. The typical Jeans wavelength and the Alfvén velocity can be estimated as m and m/s. We can also estimate the group velocity as m/s. These velocities are close to the observational velocity km/s of slowly moving perturbations. The instability growth rate can be significantly high by the effects of the magnetic field and density inhomogeneities where the Alfvén velocity gets reduced. Note that the above results are valid for low- plasmas where the length scale of magnetic field inhomogeneity is much smaller than that of the density inhomogeneity and the wavelength associated with the fast carrier magnetosonic mode is small compared to the wavelength of the slow modulation along the magnetic field. On the other hand, the typical growth rate can be high having no cut-off at any wavelength of modulation, and the corresponding time scale can be low in the case when and (Case III in Sec. IV.2). However, to be consistent with the small perturbations, the admissible growth rate should be at a wavelength (of slow modulation) above its critical value.
It has been established that the fast magnetosonic mode can act as a candidate for coronal loop heating through the exchange of energy and momentum [35]. The present theory of modulational instability should help predict the energy transfer rate of slowly varying magnetosonic envelopes in solar plasmas under the relative influences of the self-gravity force and the Coriolis force and the excitation of fast magnetosonic carrier modes at length scales below or above the Jeans critical length.
It is to be noted that in the case of plasmas with constant gravity, the gravity force, is to be replaced by and there will be no gravitational Poisson equation. Even if it is considered for the sake of clarity, it will result into the quasineutrality condition, already assumed for the derivation of the MHD equations. In this case, the term in the dispersion equation (22) for the fast magnetosonic carrier waves will be replaced by , which will result into the real wave frequency, given by,
(43) |
and the instability growth rate, . From the inhomogeneous equilibrium state one can obtain by the same assumption as for the fluid density in Sec. II.1 as . The coefficients of the group velocity dispersion and the nonlinearity in the NLS equation will be complex, given by,
(44) |
(45) |
where and . Consequently, instead of the purely growing instability as in Case III of Sec. IV.3, we have the frequency up-shift and the modulational instability growth rate , given by,
(46) |
(47) |
where the expressions for and are
(48) |
(49) |
For brevity, we have presented the characteristics of and as shown in Fig. 7. It is found that both assume higher values at a reduced angular frequency (the dotted lines) and by the effects of the density inhomogeneity (the dashed lines). Also, the growth rate is similar to Case III of Sec. IV.3 in self-gravitating plasmas, i.e., it grows with the wave number of modulation without any cut-offs.
Thus, the instabilities of fast magnetosonic modes and the MI of solwly varying wave envelopes under the influences of self-gravity force and the constant gravity force are quite distinctive. While the Jeans critical wavelength appears above which the Jeans instability of fast carrier modes occurs and hence the MI of magnetosonic envelopes in self-gravitating fields, the fast modes appears to be unstable by the constant gravity force without any restriction of the wavelength unless the length at which the collective behaviors disappear. Furthermore, in contrast to plasmas under constant gravity, the stable Jeans mode in self-gravitating fields may give rise to the MI of slowly varying envelopes with a finite growth rate and cut-off at a finite wave number of modulation (See Case II in Sec. IV.2).
VI Summary and conclusion
To summarize, we have studied the dispersion properties of fast magnetosonic waves and the modulation of slowly varying magnetosonic wave envelopes in inhomogeneous magnetoplasmas under the relative influence of the gravitational force and the Coriolis force due to rotating fluids. Using the multiple-scale reductive perturbation technique, we have derived the NLS equation for the evolution of slowly varying magnetosonic envelopes. It is found that the inhomogeneities in the background magnetic field and the fluid density favor the Jeans instability to occur in a larger domain of the carrier wave number below the Jeans critical wave number . The latter is, however, defined when the self-gravity force dominates over the Coriolis force or the Jean frequency becomes larger than the increased angular frequency . On the other hand, when the Coriolis force dominates over the self-gravity force or is dominated by the self-gravity force with , the fast magnetosonic carrier mode is always stable. It is interesting to note that the MSW envelope, corresponding to the unstable carrier Jeans mode with and , is always unstable under the plane wave modulation having no cut-offs in the growth rates. However, the stable Jeans mode with but can lead to modulational stability or instability with finite growth rates and cut-offs. Such an instability growth rate is enhanced by the effects of the magnetic field or density inhomogeneities. We have also discussed the possible applications of modulational instability in solar plasmas, such as those in the X-ray corona and solar prominence, and noted that the estimated Alfvén velocity of fast modes and the group velocity of slowly varying envelopes are close to the observational values. Also, the estimated growth time is found to be reasonable to explain the transition time from the fast magnetosonic modes to slowly varying envelopes at the nonlinear stage.
As an illustration, we have also examined the influence of the constant gravity force instead of the self-gravitation. It is found that the fast carrier magnetosonic mode is always unstable, leading to a frequency up-shift of the frequency of modulation (instead of a purely unstable mode) and the instability growth rate having no cut-offs, similar to Case III of self-gravitating plasmas.
To conclude, since the fast magnetosonic modes can be a candidate for coronal loop heating in which the momentum and energy transfer occur, the theoretical results should help predict the growth time of slowly varying magnetosonic envelopes that can be transformed from the fast magnetosonic modes emanating from disrupted filaments in the X-ray corona through nonlinear interactions. The higher the growth rate shorter is the time scale for the instability. It is to be noted that the proposed coronal heating mechanism has been the collisionless Landau damping. The fast magnetosonic modes propagating parallel or perpendicular to the magnetic field cannot suffer collisionless damping, but may be damped due to viscosity effects [27]. However, it has been shown that modulational instability is the dominant mechanism for energy transfer in the coronal loops than the viscous damping [27]. The detailed investigation in this direction is beyond the scope of the present study. Accordingly, multi-dimensional propagation of magnetosonic waves in dissipative plasmas could be a project of future study.
Acknowledgments
One of us, J. Turi, wishes to thank the Council of Scientific and Industrial Research (CSIR) for a Senior Research Fellowship (SRF) with reference number 09/202(0115)/2020-EMR-I. The authors also thank Sima Roy of the University of Engineering & Management, Kolkata, for rechecking the derivations of different expressions for the NLS equation.
Author declarations
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
J. Turi: Writing–draft, Methodology, Investigation, Formal analysis. A. P. Misra: Writing–draft, review & editing, Validation, Methodology, Investigation, Formal analysis, Conceptualization.
Data availability statement
All data that support the findings of this study are included within the article (and any supplementary files).
n-th order reduced equations
We apply the transformations [Eqs. (17) and (III)] and substitute the expansions [Eq. (19)] into Eqs. (5)-(9), to get the following set of reduced n-th order equations:
(50) |
(51) |
(52) |
(53) |
(54) |
References
- [1] H. Alfvén, Existence of electromagnetic-hydrodynamic waves, Nature 150 (3805) (1942) 405–406.
- [2] A. Hannan, T. Hellsten, T. Johnson, Fast wave current drive scenarios for demo, Nuclear Fusion 53 (4) (2013) 043005.
- [3] K. Stasiewicz, P. Bellan, C. Chaston, C. Kletzing, R. Lysak, J. Maggs, O. Pokhotelov, C. Seyler, P. Shukla, L. Stenflo, et al., Small scale alfvénic structure in the aurora, Space Science Reviews 92 (2000) 423–533.
- [4] K. Stasiewicz, P. Shukla, G. Gustafsson, S. Buchert, B. Lavraud, B. Thidé, Z. Klos, Slow magnetosonic solitons detected by the cluster spacecraft, Physical review letters 90 (8) (2003) 085002.
- [5] K. Stasiewicz, Theory and observations of slow-mode solitons in space plasmas, Physical review letters 93 (12) (2004) 125004.
- [6] P. K. Shukla, L. Stenflo, R. Bingham, B. Eliasson, Nonlinear effects associated with dispersive alfvén waves in plasmas, Plasma physics and controlled fusion 46 (12B) (2004) B349.
- [7] P. Shukla, B. Eliasson, L. Stenflo, Alfvénic shock waves in a collisional magnetoplasma, Physics Letters A 375 (24) (2011) 2371–2373.
- [8] K. Klein, G. Howes, J. TenBarge, S. Bale, C. Chen, C. Salem, Using synthetic spacecraft data to interpret compressible fluctuations in solar wind turbulence, The Astrophysical Journal 755 (2) (2012) 159.
- [9] J. Schmidt, L. Ofman, Slow magnetoacoustic wave oscillation of an expanding coronal loop, The Astrophysical Journal 739 (2) (2011) 75.
- [10] B. Rau, T. Tajima, Strongly nonlinear magnetosonic waves and ion acceleration, Physics of Plasmas 5 (10) (1998) 3575–3580.
- [11] M. Marklund, B. Eliasson, P. K. Shukla, Magnetosonic solitons in a fermionic quantum plasma, Physical Review E 76 (6) (2007) 067401.
- [12] F. Haas, S. Mahmood, Magnetosonic waves in a quantum plasma with arbitrary electron degeneracy, Physical Review E 97 (6) (2018) 063206.
- [13] W. Masood, R. Jahangir, B. Eliasson, M. Siddiq, A nonlinear model for magnetoacoustic waves in dense dissipative plasmas with degenerate electrons, Physics of Plasmas 21 (10) (2014).
- [14] S. Hussain, S. Mahmood, Korteweg-de vries burgers equation for magnetosonic wave in plasma, Physics of Plasmas 18 (5) (2011).
- [15] M. Salahuddin, H. Saleem, M. Saddiq, Ion-acoustic envelope solitons in electron-positron-ion plasmas, Physical Review E 66 (3) (2002) 036407.
- [16] S. Sultana, I. Kourakis, Electron-scale electrostatic solitary waves and shocks: the role of superthermal electrons, The European Physical Journal D 66 (2012) 1–12.
- [17] M. McKerr, I. Kourakis, F. Haas, Freak waves and electrostatic wavepacket modulation in a quantum electron–positron–ion plasma, Plasma Physics and Controlled Fusion 56 (3) (2014) 035007.
- [18] N. Chowdhury, A. Mannan, A. Mamun, Rogue waves in space dusty plasmas, Physics of Plasmas 24 (11) (2017).
- [19] M. S. Ruderman, Freak waves in laboratory and space plasmas: Freak waves in plasmas, The European Physical Journal Special Topics 185 (1) (2010) 57–66.
-
[20]
S. R. Habbal, E. Leer, T. E. Holzer,
Heating of coronal loops by fast
mode mhd waves, Solar Physics 64 (2) (1979) 287–301.
doi:10.1007/BF00151440.
URL https://doi.org/10.1007/BF00151440 - [21] S. Watanabe, Self-modulation of a nonlinear ion wave packet, Journal of Plasma Physics 17 (3) (1977) 487–501.
- [22] S. Sultana, I. Kourakis, Electrostatic solitary waves in the presence of excess superthermal electrons: modulational instability and envelope soliton modes, Plasma Physics and Controlled Fusion 53 (4) (2011) 045003.
- [23] S. Shalini, N. Saini, A. Misra, Modulation of ion-acoustic waves in a nonextensive plasma with two-temperature electrons, Physics of Plasmas 22 (9) (2015).
- [24] A. P. Misra, C. Bhowmik, Nonlinear wave modulation in a quantum magnetoplasma, Physics of plasmas 14 (1) (2007).
- [25] A. Bains, A. P. Misra, N. Saini, T. Gill, Modulational instability of ion-acoustic wave envelopes in magnetized quantum electron-positron-ion plasmas, Physics of Plasmas 17 (1) (2010).
- [26] W. Sahyouni, I. Zhelyazkov, P. Nenovski, Dark envelope solitons of fast magnetosonic surface waves in solar flux tubes, Solar physics 115 (1988) 17–32.
- [27] J.-I. Sakai, Modulational instability of fast magnetosonic waves in a solar plasma, Solar Physics 84 (1983) 109–118.
- [28] A. Misra, P. K. Shukla, Modulational instability of magnetosonic waves in a spin 1/ 2 quantum plasma, Physics of Plasmas 15 (5) (2008).
- [29] A. Panwar, C.-M. Ryu, Modulational instability and associated rogue structures of slow magnetosonic wave in hall magnetohydrodynamic plasmas, Physics of Plasmas 21 (6) (2014).
- [30] Y. Wang, X. Lü, B. Eliasson, Modulational instability of spin modified quantum magnetosonic waves in fermi-dirac-pauli plasmas, Physics of Plasmas 20 (11) (2013).
- [31] J. Turi, A. Misra, Magnetohydrodynamic instabilities in a self-gravitating rotating cosmic plasma, Physica Scripta 97 (12) (2022) 125603.
- [32] N. Asano, T. Taniuti, N. Yajima, Perturbation method for a nonlinear wave modulation. ii, Journal of Mathematical Physics 10 (11) (1969) 2020–2024.
- [33] J. H. Jeans, I. the stability of a spherical nebula, Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character 199 (312-320) (1902) 1–53.
- [34] Y. H. Ichikawa, C. nonlinear wave modulation of electrostatic waves, Progress of Theoretical Physics Supplement 55 (1974) 212–232.
-
[35]
D. Y. Kolotkov, D. I. Zavershinskii, V. M. Nakariakov,
The solar corona as an
active medium for magnetoacoustic waves, Plasma Physics and Controlled
Fusion 63 (12) (2021) 124008.
doi:10.1088/1361-6587/ac36a5.
URL https://dx.doi.org/10.1088/1361-6587/ac36a5 -
[36]
L. Ofman, T. A. Kucera, C. R. DeVore,
Nonlinear fast
magnetosonic waves in solar prominence pillars, The Astrophysical Journal
944 (2) (2023) 210.
doi:10.3847/1538-4357/acb13b.
URL https://dx.doi.org/10.3847/1538-4357/acb13b -
[37]
D. M. Rust, Z. Švestka, Slowly
moving disturbances in the x-ray corona, Solar Physics 63 (2) (1979)
279–295.
doi:10.1007/BF00174535.
URL https://doi.org/10.1007/BF00174535