Nonlinear modulation of dispersive fast magnetosonic waves in an inhomogeneous rotating solar low-β𝛽\betaitalic_β magnetoplasma

Jyoti Turi jyotituri.maths@gmail.com Department of Mathematics, Siksha Bhavana, Visva-Bharati University, Santiniketan-731 235, West Bengal, India    A. P. Misra apmisra@visva-bharati.ac.in Department of Mathematics, Siksha Bhavana, Visva-Bharati University, Santiniketan-731 235, West Bengal, India
(June 19, 2024)
Abstract

We study the modulation of fast magnetosonic waves (MSWs) in rotating inhomogeneous low-β𝛽\betaitalic_β 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 (k𝑘kitalic_k, below the Jeans critical wave number, kJsubscript𝑘𝐽k_{J}italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT) than homogeneous plasmas. The relative influence of the Jeans frequency (ωJsubscript𝜔𝐽\omega_{J}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT, associated with the gravitational force) and the angular frequency (Ω0subscriptΩ0\Omega_{0}roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, 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 ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT) is always unstable against the plane wave perturbation with no cut-offs for growth rates. In contrast, the stable Jeans mode with ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT but k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT 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 1977197719771977. 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 g𝑔gitalic_g is not a constant (self-gravity) and (ii) when g𝑔gitalic_g is a constant. The plasma is supposed to be rotating with uniform angular velocity 𝛀=(0,Ω0cosλ,Ω0sinλ)𝛀0subscriptΩ0𝜆subscriptΩ0𝜆\mathbf{\Omega}=(0,\Omega_{0}\cos\lambda,\Omega_{0}\sin\lambda)bold_Ω = ( 0 , roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_λ , roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_λ ) and immersed in an external static magnetic field along the z𝑧zitalic_z-axis, i.e.𝐁𝟎=B0z^subscript𝐁0subscript𝐵0^𝑧\mathbf{B_{0}}=B_{0}\hat{z}bold_B start_POSTSUBSCRIPT bold_0 end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT over^ start_ARG italic_z end_ARG. We assume the gravitational force to be acting vertically downwards and parpendicular to the magnetic field and the wave propagation along the x𝑥xitalic_x-axis, for simplicity. Furthermore, the background magnetic field and the density are assumed to vary along the x𝑥xitalic_x-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]

ρt+(ρ𝐯)=0,𝜌𝑡𝜌𝐯0\frac{\partial\rho}{\partial t}+\mathbf{\nabla}\cdot\left(\rho\mathbf{v}\right% )=0,divide start_ARG ∂ italic_ρ end_ARG start_ARG ∂ italic_t end_ARG + ∇ ⋅ ( italic_ρ bold_v ) = 0 , (1)
𝐯t+(𝐯)𝐯=1ρ(P+𝐁22μ0)+1ρμ0(𝐁)𝐁2𝛀×𝐯+𝐠,𝐯𝑡𝐯𝐯1𝜌𝑃superscript𝐁22subscript𝜇01𝜌subscript𝜇0𝐁𝐁2𝛀𝐯𝐠\frac{\partial\mathbf{v}}{\partial t}+\left(\mathbf{v}\cdot\mathbf{\nabla}% \right)\mathbf{v}=-\frac{1}{\rho}\mathbf{\nabla}\left(P+\frac{\mathbf{B}^{2}}{% 2\mu_{0}}\right)\\ +\frac{1}{\rho\mu_{0}}\left(\mathbf{B}\cdot\mathbf{\nabla}\right)\mathbf{B}-2% \mathbf{\Omega}\times\mathbf{v}+\mathbf{g},start_ROW start_CELL divide start_ARG ∂ bold_v end_ARG start_ARG ∂ italic_t end_ARG + ( bold_v ⋅ ∇ ) bold_v = - divide start_ARG 1 end_ARG start_ARG italic_ρ end_ARG ∇ ( italic_P + divide start_ARG bold_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) end_CELL end_ROW start_ROW start_CELL + divide start_ARG 1 end_ARG start_ARG italic_ρ italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ( bold_B ⋅ ∇ ) bold_B - 2 bold_Ω × bold_v + bold_g , end_CELL end_ROW (2)
𝐁t+(𝐯)𝐁=(B)𝐯(𝐯)𝐁,𝐁𝑡𝐯𝐁B𝐯𝐯𝐁\frac{\partial\mathbf{B}}{\partial t}+\left(\mathbf{v}\cdot\mathbf{\nabla}% \right)\mathbf{B}=\left(\mathbf{\textbf{B}}\cdot\mathbf{\nabla}\right)\mathbf{% v}-\left(\mathbf{\nabla}\cdot\mathbf{v}\right)\mathbf{B},divide start_ARG ∂ bold_B end_ARG start_ARG ∂ italic_t end_ARG + ( bold_v ⋅ ∇ ) bold_B = ( B ⋅ ∇ ) bold_v - ( ∇ ⋅ bold_v ) bold_B , (3)
𝐠=4πGρ,𝐠4𝜋𝐺𝜌\mathbf{\nabla}\cdot\mathbf{g}=-4\pi G\rho,∇ ⋅ bold_g = - 4 italic_π italic_G italic_ρ , (4)

where 𝐠=ψ𝐠𝜓\mathbf{g}=-\mathbf{\nabla}\psibold_g = - ∇ italic_ψ is the gravitational force per unit mass of the fluid acting vertically downwards; ρ𝜌\rhoitalic_ρ, 𝐯𝐯\mathbf{v}bold_v, and P𝑃Pitalic_P are, respectively, the fluid density, fluid velocity, and the thermal pressure. Also, 𝐁𝐁\mathbf{B}bold_B, 𝛀𝛀\mathbf{\Omega}bold_Ω, ψ𝜓\psiitalic_ψ, and G𝐺Gitalic_G, respectively, denote the magnetic field, uniform angular velocity of the rotating fluid, gravitational potential, and the universal gravitational constant. The pressure P𝑃Pitalic_P satisfies the equation of state: P=cs2ρ𝑃superscriptsubscript𝑐𝑠2𝜌\mathbf{\nabla}P=c_{s}^{2}\mathbf{\nabla}\rho∇ italic_P = italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∇ italic_ρ, where cs=γkBTe/misubscript𝑐𝑠𝛾subscript𝑘𝐵subscript𝑇𝑒subscript𝑚𝑖c_{s}=\sqrt{\gamma k_{B}T_{e}/m_{i}}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = square-root start_ARG italic_γ italic_k start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT / italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG is the adiabatic sound speed with kBsubscript𝑘𝐵k_{B}italic_k start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT denoting the Boltzmann constant, Tesubscript𝑇𝑒T_{e}italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT the electron temperature, γ𝛾\gammaitalic_γ the adiabatic index, and misubscript𝑚𝑖m_{i}italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT the ion mass.

Refer to caption
Figure 1: A schematic diagram showing the directions of wave propagation 𝐤𝐤\mathbf{k}bold_k, the external magnetic field 𝐁𝟎subscript𝐁0\mathbf{B_{0}}bold_B start_POSTSUBSCRIPT bold_0 end_POSTSUBSCRIPT, and the gravitational force 𝐠𝐠\mathbf{g}bold_g. The background plasma state is determined by Eqs. (10), (11). The magnetic field and density inhomogeneities are along the x𝑥xitalic_x-axis.

Next, we consider the one-dimensional propagation of fast MSWs along the x𝑥xitalic_x-axis and normalize Eqs. (1)-(4) according to: ρρ/ρ0𝜌𝜌subscript𝜌0\rho\rightarrow\rho/\rho_{0}italic_ρ → italic_ρ / italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, BB/B0𝐵𝐵subscript𝐵0B\rightarrow B/B_{0}italic_B → italic_B / italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, ωω/ωci𝜔𝜔subscript𝜔𝑐𝑖\omega\rightarrow\omega/\omega_{ci}italic_ω → italic_ω / italic_ω start_POSTSUBSCRIPT italic_c italic_i end_POSTSUBSCRIPT, (vx,vy,vz)(vx,vy,vz)/VAsubscript𝑣𝑥subscript𝑣𝑦subscript𝑣𝑧subscript𝑣𝑥subscript𝑣𝑦subscript𝑣𝑧subscript𝑉𝐴\left(v_{x},v_{y},v_{z}\right)\rightarrow\left(v_{x},v_{y},v_{z}\right)/V_{A}( italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) → ( italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT , italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) / italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT, cscs/VAsubscript𝑐𝑠subscript𝑐𝑠subscript𝑉𝐴{c}_{s}\rightarrow c_{s}/V_{A}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT → italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT / italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT, Ω0Ω0/ωcisubscriptΩ0subscriptΩ0subscript𝜔𝑐𝑖{\Omega}_{0}\rightarrow\Omega_{0}/\omega_{ci}roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT → roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / italic_ω start_POSTSUBSCRIPT italic_c italic_i end_POSTSUBSCRIPT, ωJωJ/ωcisubscript𝜔𝐽subscript𝜔𝐽subscript𝜔𝑐𝑖{\omega}_{J}\rightarrow\omega_{J}/\omega_{ci}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT → italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT / italic_ω start_POSTSUBSCRIPT italic_c italic_i end_POSTSUBSCRIPT, ψψ/VA2𝜓𝜓superscriptsubscript𝑉𝐴2\psi\rightarrow\psi/V_{A}^{2}italic_ψ → italic_ψ / italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Here, VA=B0/μ0ρ0subscript𝑉𝐴subscript𝐵0subscript𝜇0subscript𝜌0V_{A}={B_{0}}/{\sqrt{\mu_{0}\rho_{0}}}italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT = italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / square-root start_ARG italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG is the Alfvén speed, ωJ2=4πGρ0superscriptsubscript𝜔𝐽24𝜋𝐺subscript𝜌0\omega_{J}^{2}=4\pi G\rho_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 4 italic_π italic_G italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the squared Jeans frequency, and ωci=eB0/misubscript𝜔𝑐𝑖𝑒subscript𝐵0subscript𝑚𝑖\omega_{ci}=\sqrt{eB_{0}/m_{i}}italic_ω start_POSTSUBSCRIPT italic_c italic_i end_POSTSUBSCRIPT = square-root start_ARG italic_e italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG is the ion-cyclotron frequency at x=0𝑥0x=0italic_x = 0. Furthermore, the space and time coordinates are normalized as xxωci/VA𝑥𝑥subscript𝜔𝑐𝑖subscript𝑉𝐴x\rightarrow x\omega_{ci}/V_{A}italic_x → italic_x italic_ω start_POSTSUBSCRIPT italic_c italic_i end_POSTSUBSCRIPT / italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT and ttωci𝑡𝑡subscript𝜔𝑐𝑖t\rightarrow t\omega_{ci}italic_t → italic_t italic_ω start_POSTSUBSCRIPT italic_c italic_i end_POSTSUBSCRIPT. Thus, from Eqs. (1)-(4), after separating the velocity components along the axes and noting that the frozen-in-field condition ρ/ρ0=B/B0𝜌subscript𝜌0𝐵subscript𝐵0\rho/\rho_{0}=B/B_{0}italic_ρ / italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_B / italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT holds for Eqs. (1) and (3), we obtain the following reduced set of normalized equations

Bt+x(vxB)=0,𝐵𝑡𝑥subscript𝑣𝑥𝐵0\displaystyle\frac{\partial B}{\partial t}+\frac{\partial}{\partial x}\left(v_% {x}B\right)=0,divide start_ARG ∂ italic_B end_ARG start_ARG ∂ italic_t end_ARG + divide start_ARG ∂ end_ARG start_ARG ∂ italic_x end_ARG ( italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT italic_B ) = 0 , (5)
vxt+vxvxx=cs21BBx121BB2xsubscript𝑣𝑥𝑡subscript𝑣𝑥subscript𝑣𝑥𝑥superscriptsubscript𝑐𝑠21𝐵𝐵𝑥121𝐵superscript𝐵2𝑥\displaystyle\frac{\partial v_{x}}{\partial t}+v_{x}\frac{\partial v_{x}}{% \partial x}=-{c}_{s}^{2}\frac{1}{B}\frac{\partial B}{\partial x}-\frac{1}{2}% \frac{1}{B}\frac{\partial B^{2}}{\partial x}divide start_ARG ∂ italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_t end_ARG + italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT divide start_ARG ∂ italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_x end_ARG = - italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG 1 end_ARG start_ARG italic_B end_ARG divide start_ARG ∂ italic_B end_ARG start_ARG ∂ italic_x end_ARG - divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG 1 end_ARG start_ARG italic_B end_ARG divide start_ARG ∂ italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_x end_ARG
2(Ω0vzcosλΩ0vysinλ)+ψx,2subscriptΩ0subscript𝑣𝑧𝜆subscriptΩ0subscript𝑣𝑦𝜆𝜓𝑥\displaystyle-2\left({\Omega}_{0}v_{z}\cos\lambda-{\Omega}_{0}v_{y}\sin\lambda% \right)+\frac{\partial\psi}{\partial x},- 2 ( roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT roman_cos italic_λ - roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT roman_sin italic_λ ) + divide start_ARG ∂ italic_ψ end_ARG start_ARG ∂ italic_x end_ARG , (6)
vyt+vxvyx=2Ω0vxsinλ,subscript𝑣𝑦𝑡subscript𝑣𝑥subscript𝑣𝑦𝑥2subscriptΩ0subscript𝑣𝑥𝜆\displaystyle\frac{\partial v_{y}}{\partial t}+v_{x}\frac{\partial v_{y}}{% \partial x}=-2{\Omega}_{0}v_{x}\sin\lambda,divide start_ARG ∂ italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_t end_ARG + italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT divide start_ARG ∂ italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_x end_ARG = - 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT roman_sin italic_λ , (7)
vzt+vxvzx=2Ω0vxcosλ,subscript𝑣𝑧𝑡subscript𝑣𝑥subscript𝑣𝑧𝑥2subscriptΩ0subscript𝑣𝑥𝜆\displaystyle\frac{\partial v_{z}}{\partial t}+v_{x}\frac{\partial v_{z}}{% \partial x}=2{\Omega}_{0}v_{x}\cos\lambda,divide start_ARG ∂ italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_t end_ARG + italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT divide start_ARG ∂ italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG ∂ italic_x end_ARG = 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT roman_cos italic_λ , (8)
2ψx2+ωJ2B=0.superscript2𝜓superscript𝑥2superscriptsubscript𝜔𝐽2𝐵0\displaystyle\frac{\partial^{2}\psi}{\partial x^{2}}+{\omega}_{J}^{2}B=0.divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ψ end_ARG start_ARG ∂ italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B = 0 . (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

ddx(P0(x)+B02(x)2μ0)=ρ0(x)dψ0(x)dx,𝑑𝑑𝑥subscript𝑃0𝑥superscriptsubscript𝐵02𝑥2subscript𝜇0subscript𝜌0𝑥𝑑subscript𝜓0𝑥𝑑𝑥\frac{d}{dx}\left(P_{0}(x)+\frac{B_{0}^{2}(x)}{2\mu_{0}}\right)=-\rho_{0}(x)% \frac{d\psi_{0}(x)}{dx},divide start_ARG italic_d end_ARG start_ARG italic_d italic_x end_ARG ( italic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) + divide start_ARG italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) end_ARG start_ARG 2 italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ) = - italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG italic_d italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG italic_d italic_x end_ARG , (10)
d2ψ0(x)dx2=4πGρ0(x),superscript𝑑2subscript𝜓0𝑥𝑑superscript𝑥24𝜋𝐺subscript𝜌0𝑥\frac{d^{2}\psi_{0}(x)}{dx^{2}}=4\pi G\rho_{0}(x),divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG italic_d italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = 4 italic_π italic_G italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) , (11)

where the suffix ‘00’ denotes the equilibrium value of the corresponding physical quantity. We define B0(0)=B0subscript𝐵00subscript𝐵0B_{0}(0)=B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 0 ) = italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, ρ0(0)=ρ0subscript𝜌00subscript𝜌0\rho_{0}(0)=\rho_{0}italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 0 ) = italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and ψ0(0)=ψ0subscript𝜓00subscript𝜓0\psi_{0}(0)=\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 0 ) = italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. By means of the normalization defined before, Eqs. (10) and (11) reduce to

ddx(cs2ρ0(x)+12B02(x))=ρ0(x)dψ0(x)dx,𝑑𝑑𝑥superscriptsubscript𝑐𝑠2subscript𝜌0𝑥12superscriptsubscript𝐵02𝑥subscript𝜌0𝑥𝑑subscript𝜓0𝑥𝑑𝑥\frac{d}{dx}\left({c}_{s}^{2}\rho_{0}(x)+\frac{1}{2}B_{0}^{2}(x)\right)=-\rho_% {0}(x)\frac{d\psi_{0}(x)}{dx},divide start_ARG italic_d end_ARG start_ARG italic_d italic_x end_ARG ( italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ) = - italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG italic_d italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG italic_d italic_x end_ARG , (12)
d2ψ0(x)dx2=ωJ2ρ0(x).superscript𝑑2subscript𝜓0𝑥𝑑superscript𝑥2superscriptsubscript𝜔𝐽2subscript𝜌0𝑥\frac{d^{2}\psi_{0}(x)}{dx^{2}}={\omega}_{J}^{2}\rho_{0}(x).divide start_ARG italic_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG italic_d italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) . (13)

Using the relation P=cs2ρ𝑃superscriptsubscript𝑐𝑠2𝜌\nabla P=c_{s}^{2}\nabla\rho∇ italic_P = italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∇ italic_ρ, and noting that B0(x)/B0=ρ0(x)/ρ0subscript𝐵0𝑥subscript𝐵0subscript𝜌0𝑥subscript𝜌0B_{0}(x)/B_{0}=\rho_{0}(x)/\rho_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) / italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) / italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT (Frozen-in-field condition), we obtain from Eq. (12), the following relation (in normalized form)

cs2lnρ0(x)+2ρ0(x)=ψ~0ψ0(x),superscriptsubscript𝑐𝑠2subscript𝜌0𝑥2subscript𝜌0𝑥subscript~𝜓0subscript𝜓0𝑥c_{s}^{2}\ln\rho_{0}(x)+2\rho_{0}(x)=\tilde{\psi}_{0}-\psi_{0}(x),italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_ln italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) + 2 italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) = over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) , (14)

where ψ~0=2+ψ0subscript~𝜓02subscript𝜓0\tilde{\psi}_{0}=2+\psi_{0}over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 2 + italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Typically, if ρ0(x)subscript𝜌0𝑥\rho_{0}(x)italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) is not too small, i.e., if it lies in 0.5ρ0(x)1less-than-or-similar-to0.5subscript𝜌0𝑥less-than-or-similar-to10.5\lesssim\rho_{0}(x)\lesssim 10.5 ≲ italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ≲ 1, the function lnρ0(x)/ρ0(x)subscript𝜌0𝑥subscript𝜌0𝑥\ln\rho_{0}(x)/\rho_{0}(x)roman_ln italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) / italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) approaches from small negative values to zero as ρ0(x)1subscript𝜌0𝑥1\rho_{0}(x)\rightarrow 1italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) → 1. In addition, if the magnetic force dominates over the pressure gradient force (in low-β𝛽\betaitalic_β 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., LB0Lρ0much-less-thansubscript𝐿subscript𝐵0subscript𝐿subscript𝜌0L_{B_{0}}\ll L_{\rho_{0}}italic_L start_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ≪ italic_L start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT, where 1/Lρ0[1/ρ0(x)](dρ0(x)/dx)1subscript𝐿subscript𝜌0delimited-[]1subscript𝜌0𝑥𝑑subscript𝜌0𝑥𝑑𝑥1/L_{\rho_{0}}\equiv[1/\rho_{0}(x)](d\rho_{0}(x)/dx)1 / italic_L start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ≡ [ 1 / italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ] ( italic_d italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) / italic_d italic_x ) and 1/LB0[1/B0(x)](dB0(x)/dx)1subscript𝐿subscript𝐵0delimited-[]1subscript𝐵0𝑥𝑑subscript𝐵0𝑥𝑑𝑥1/L_{B_{0}}\equiv[1/B_{0}(x)](dB_{0}(x)/dx)1 / italic_L start_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ≡ [ 1 / italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ] ( italic_d italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) / italic_d italic_x ), then the term proportional to cs2superscriptsubscript𝑐𝑠2c_{s}^{2}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT in Eq. (14) can be neglected compared to the the term involving the factor 2222 (the second term on the left-hand side). Thus, from Eqs. (13) and (14), we obtain for x>0𝑥0x>0italic_x > 0 the following approximate solutions for ψ0(x)subscript𝜓0𝑥\psi_{0}(x)italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) and ρ0(x)subscript𝜌0𝑥\rho_{0}(x)italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ).

ψ0(x)ψ~02exp(12ωJ2x),ρ0(x)B0(x)12[ψ~0ψ0(x)].formulae-sequencesubscript𝜓0𝑥subscript~𝜓0212superscriptsubscript𝜔𝐽2𝑥subscript𝜌0𝑥subscript𝐵0𝑥12delimited-[]subscript~𝜓0subscript𝜓0𝑥\begin{split}&\psi_{0}(x)\approx\tilde{\psi}_{0}-2\exp\left(-\frac{1}{2}\omega% _{J}^{2}x\right),\\ &\rho_{0}(x)\equiv B_{0}(x)\approx\frac{1}{2}\left[\tilde{\psi}_{0}-\psi_{0}(x% )\right].\end{split}start_ROW start_CELL end_CELL start_CELL italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ≈ over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - 2 roman_exp ( - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x ) , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ≡ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ≈ divide start_ARG 1 end_ARG start_ARG 2 end_ARG [ over~ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ] . end_CELL end_ROW (15)

On the other hand, if the contribution from the pressure gradient force is much higher than the magnetic force (in high-β𝛽\betaitalic_β 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., LB0Lρ0much-greater-thansubscript𝐿subscript𝐵0subscript𝐿subscript𝜌0L_{B_{0}}\gg L_{\rho_{0}}italic_L start_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ≫ italic_L start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT, approximate solutions for ρ0(x)subscript𝜌0𝑥\rho_{0}(x)italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) and ψ0(x)subscript𝜓0𝑥\psi_{0}(x)italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) can then be obtained from Eqs. (13) and (14) as (x>0)𝑥0(x>0)( italic_x > 0 )

ψ0(x)cs2[1ψ¯0exp(ωj2cs2x)],ρ0(x)B0(x)exp(ψ0ψ0(x)cs2),formulae-sequencesubscript𝜓0𝑥superscriptsubscript𝑐𝑠2delimited-[]1subscript¯𝜓0superscriptsubscript𝜔𝑗2superscriptsubscript𝑐𝑠2𝑥subscript𝜌0𝑥subscript𝐵0𝑥subscript𝜓0subscript𝜓0𝑥superscriptsubscript𝑐𝑠2\begin{split}&\psi_{0}(x)\approx c_{s}^{2}\left[1-\bar{\psi}_{0}\exp\left(-% \frac{\omega_{j}^{2}}{c_{s}^{2}}x\right)\right],\\ &\rho_{0}(x)\equiv B_{0}(x)\approx\exp\left(\frac{\psi_{0}-\psi_{0}(x)}{c_{s}^% {2}}\right),\end{split}start_ROW start_CELL end_CELL start_CELL italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ≈ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ 1 - over¯ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_exp ( - divide start_ARG italic_ω start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_x ) ] , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ≡ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ≈ roman_exp ( divide start_ARG italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) , end_CELL end_ROW (16)

where ψ¯0=1ψ0/cs2subscript¯𝜓01subscript𝜓0superscriptsubscript𝑐𝑠2\bar{\psi}_{0}=1-\psi_{0}/c_{s}^{2}over¯ start_ARG italic_ψ end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1 - italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Furthermore, for plasmas where the plasma pressure is comparable to the magnetic pressure (β1)similar-to𝛽1(\beta\sim 1)( italic_β ∼ 1 ), 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

ξ=ϵ(xvgt),τ=ϵ2t,formulae-sequence𝜉italic-ϵ𝑥subscript𝑣𝑔𝑡𝜏superscriptitalic-ϵ2𝑡\begin{split}&\xi=\epsilon\left(x-v_{g}t\right),\\ &\tau=\epsilon^{2}t,\end{split}start_ROW start_CELL end_CELL start_CELL italic_ξ = italic_ϵ ( italic_x - italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_t ) , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_τ = italic_ϵ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t , end_CELL end_ROW (17)

where vgsubscript𝑣𝑔v_{g}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT is the group velocity of the wave envelope along the x𝑥xitalic_x-axis and ϵitalic-ϵ\epsilonitalic_ϵ is a small (0<ϵ10italic-ϵmuch-less-than10<\epsilon\ll 10 < italic_ϵ ≪ 1) 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.

ttϵvgξ+ϵ2τ,𝑡𝑡italic-ϵsubscript𝑣𝑔𝜉superscriptitalic-ϵ2𝜏\displaystyle\frac{\partial}{\partial t}\rightarrow\frac{\partial}{\partial t}% -\epsilon v_{g}\frac{\partial}{\partial\xi}+\epsilon^{2}\frac{\partial}{% \partial\tau},divide start_ARG ∂ end_ARG start_ARG ∂ italic_t end_ARG → divide start_ARG ∂ end_ARG start_ARG ∂ italic_t end_ARG - italic_ϵ italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG + italic_ϵ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_τ end_ARG ,
xx+ϵξ,𝑥𝑥italic-ϵ𝜉\displaystyle\frac{\partial}{\partial x}\rightarrow\frac{\partial}{\partial x}% +\epsilon\frac{\partial}{\partial\xi},divide start_ARG ∂ end_ARG start_ARG ∂ italic_x end_ARG → divide start_ARG ∂ end_ARG start_ARG ∂ italic_x end_ARG + italic_ϵ divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG , (18)
2x22x2+2ϵ2xξ+ϵ2ξ2.superscript2superscript𝑥2superscript2superscript𝑥22italic-ϵsuperscript2𝑥𝜉superscriptitalic-ϵ2superscript𝜉2\displaystyle\frac{\partial^{2}}{\partial x^{2}}\rightarrow\frac{\partial^{2}}% {\partial x^{2}}+2\epsilon\frac{\partial^{2}}{\partial x\partial\xi}+\epsilon^% {2}\frac{\partial}{\partial\xi^{2}}.divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG → divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 2 italic_ϵ divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_x ∂ italic_ξ end_ARG + italic_ϵ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG .

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 l𝑙litalic_l-th harmonic amplitudes and the phase (kxωt)𝑘𝑥𝜔𝑡(kx-\omega t)( italic_k italic_x - italic_ω italic_t ). Thus, the variables can be expanded as

B=B0(x)+n=1ϵnl=Bl(n)(ξ,τ)exp[il(kxωt)],vx=0+n=1ϵnl=vxl(n)(ξ,τ)exp[il(kxωt)],vy=0+n=1ϵnl=vyl(n)(ξ,τ)exp[il(kxωt)],vz=0+n=1ϵnl=vzl(n)(ξ,τ)exp[il(kxωt)],ψ=ψ0(x)+n=1ϵnl=ψl(n)(ξ,τ)exp[il(kxωt)],formulae-sequence𝐵subscript𝐵0𝑥superscriptsubscript𝑛1superscriptitalic-ϵ𝑛superscriptsubscript𝑙superscriptsubscript𝐵𝑙𝑛𝜉𝜏𝑖𝑙𝑘𝑥𝜔𝑡formulae-sequencesubscript𝑣𝑥0superscriptsubscript𝑛1superscriptitalic-ϵ𝑛superscriptsubscript𝑙superscriptsubscript𝑣𝑥𝑙𝑛𝜉𝜏𝑖𝑙𝑘𝑥𝜔𝑡formulae-sequencesubscript𝑣𝑦0superscriptsubscript𝑛1superscriptitalic-ϵ𝑛superscriptsubscript𝑙superscriptsubscript𝑣𝑦𝑙𝑛𝜉𝜏𝑖𝑙𝑘𝑥𝜔𝑡formulae-sequencesubscript𝑣𝑧0superscriptsubscript𝑛1superscriptitalic-ϵ𝑛superscriptsubscript𝑙superscriptsubscript𝑣𝑧𝑙𝑛𝜉𝜏𝑖𝑙𝑘𝑥𝜔𝑡𝜓subscript𝜓0𝑥superscriptsubscript𝑛1superscriptitalic-ϵ𝑛superscriptsubscript𝑙superscriptsubscript𝜓𝑙𝑛𝜉𝜏𝑖𝑙𝑘𝑥𝜔𝑡\begin{split}&B=B_{0}(x)+\sum_{n=1}^{\infty}\epsilon^{n}\sum_{l=-\infty}^{% \infty}B_{l}^{(n)}(\xi,\tau)\exp\left[il\left(kx-\omega t\right)\right],\\ &v_{x}=0+\sum_{n=1}^{\infty}\epsilon^{n}\sum_{l=-\infty}^{\infty}v_{xl}^{(n)}(% \xi,\tau)\exp\left[il\left(kx-\omega t\right)\right],\\ &v_{y}=0+\sum_{n=1}^{\infty}\epsilon^{n}\sum_{l=-\infty}^{\infty}v_{yl}^{(n)}(% \xi,\tau)\exp\left[il\left(kx-\omega t\right)\right],\\ &v_{z}=0+\sum_{n=1}^{\infty}\epsilon^{n}\sum_{l=-\infty}^{\infty}v_{zl}^{(n)}(% \xi,\tau)\exp\left[il\left(kx-\omega t\right)\right],\\ &\psi=\psi_{0}(x)+\sum_{n=1}^{\infty}\epsilon^{n}\sum_{l=-\infty}^{\infty}\psi% _{l}^{(n)}(\xi,\tau)\exp\left[il\left(kx-\omega t\right)\right],\end{split}start_ROW start_CELL end_CELL start_CELL italic_B = italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) + ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_ϵ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT ( italic_ξ , italic_τ ) roman_exp [ italic_i italic_l ( italic_k italic_x - italic_ω italic_t ) ] , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 0 + ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_ϵ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_x italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT ( italic_ξ , italic_τ ) roman_exp [ italic_i italic_l ( italic_k italic_x - italic_ω italic_t ) ] , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 0 + ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_ϵ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_y italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT ( italic_ξ , italic_τ ) roman_exp [ italic_i italic_l ( italic_k italic_x - italic_ω italic_t ) ] , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 + ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_ϵ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_z italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT ( italic_ξ , italic_τ ) roman_exp [ italic_i italic_l ( italic_k italic_x - italic_ω italic_t ) ] , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_ψ = italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) + ∑ start_POSTSUBSCRIPT italic_n = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_ϵ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_ψ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT ( italic_ξ , italic_τ ) roman_exp [ italic_i italic_l ( italic_k italic_x - italic_ω italic_t ) ] , end_CELL end_ROW (19)

where B0(x)subscript𝐵0𝑥B_{0}(x)italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) and ψ0(x)subscript𝜓0𝑥\psi_{0}(x)italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) are normalized by B0subscript𝐵0B_{0}italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and ψ0subscript𝜓0\psi_{0}italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT respectively, and the slowly varying wave amplitudes Bl(n)superscriptsubscript𝐵𝑙𝑛B_{l}^{(n)}italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT, vxl(n)superscriptsubscript𝑣𝑥𝑙𝑛v_{xl}^{(n)}italic_v start_POSTSUBSCRIPT italic_x italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT, vyl(n)superscriptsubscript𝑣𝑦𝑙𝑛v_{yl}^{(n)}italic_v start_POSTSUBSCRIPT italic_y italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT, vzl(n)superscriptsubscript𝑣𝑧𝑙𝑛v_{zl}^{(n)}italic_v start_POSTSUBSCRIPT italic_z italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT, ψl(n)superscriptsubscript𝜓𝑙𝑛\psi_{l}^{(n)}italic_ψ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT, etc. satisfy the reality condition: Al(n)=Al(n)superscriptsubscript𝐴𝑙𝑛superscriptsubscript𝐴𝑙𝑛A_{-l}^{(n)}=A_{l}^{(n)*}italic_A start_POSTSUBSCRIPT - italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = italic_A start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) ∗ end_POSTSUPERSCRIPT, where the asterisk denotes the complex conjugate, and k𝑘kitalic_k and ω𝜔\omegaitalic_ω 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 ϵitalic-ϵ\epsilonitalic_ϵ (For some details, see Appendix n-th order reduced equations). Here, we assume that the length scales of the density (Lρ0subscript𝐿subscript𝜌0L_{\rho_{0}}italic_L start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT) and magnetic field (LB0subscript𝐿subscript𝐵0L_{B_{0}}italic_L start_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT) inhomogeneities are much larger than the length scale (L)𝐿(L)( italic_L ) of magnetosonic perturbations, i.e., Lρ0,LB0Lmuch-greater-thansubscript𝐿subscript𝜌0subscript𝐿subscript𝐵0𝐿L_{\rho_{0}},~{}L_{B_{0}}\gg Litalic_L start_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , italic_L start_POSTSUBSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ≫ italic_L. 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 ϵitalic-ϵ\epsilonitalic_ϵ, the following relations

vx1(1)=ωkB1(1),vy1(1)=i2kΩ0sinλB1(1),vz1(1)=i2kΩ0cosλB1(1),ψ1(1)=ωJ2k2B1(1),formulae-sequencesuperscriptsubscript𝑣𝑥11𝜔𝑘superscriptsubscript𝐵11formulae-sequencesuperscriptsubscript𝑣𝑦11𝑖2𝑘subscriptΩ0𝜆superscriptsubscript𝐵11formulae-sequencesuperscriptsubscript𝑣𝑧11𝑖2𝑘subscriptΩ0𝜆superscriptsubscript𝐵11superscriptsubscript𝜓11superscriptsubscript𝜔𝐽2superscript𝑘2superscriptsubscript𝐵11\begin{split}&v_{x1}^{(1)}=\frac{\omega}{k}B_{1}^{(1)},~{}v_{y1}^{(1)}=-i\frac% {2}{k}\Omega_{0}\sin\lambda B_{1}^{(1)},\\ &v_{z1}^{(1)}=i\frac{2}{k}\Omega_{0}\cos\lambda B_{1}^{(1)},~{}\psi_{1}^{(1)}=% \frac{{\omega}_{J}^{2}}{k^{2}}B_{1}^{(1)},\end{split}start_ROW start_CELL end_CELL start_CELL italic_v start_POSTSUBSCRIPT italic_x 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = divide start_ARG italic_ω end_ARG start_ARG italic_k end_ARG italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT , italic_v start_POSTSUBSCRIPT italic_y 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = - italic_i divide start_ARG 2 end_ARG start_ARG italic_k end_ARG roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_λ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_v start_POSTSUBSCRIPT italic_z 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = italic_i divide start_ARG 2 end_ARG start_ARG italic_k end_ARG roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_λ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT , italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = divide start_ARG italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT , end_CELL end_ROW (20)
iωvx1(1)+ik(cs2+VA(x)2)B1(1)+2Ω0(cosλvz1(1)sinλvy1(1))ikψ(1)=0,𝑖𝜔superscriptsubscript𝑣𝑥11𝑖𝑘superscriptsubscript𝑐𝑠2subscript𝑉𝐴superscript𝑥2superscriptsubscript𝐵112subscriptΩ0𝜆superscriptsubscript𝑣𝑧11𝜆superscriptsubscript𝑣𝑦11𝑖𝑘superscript𝜓10\begin{split}&-i\omega v_{x1}^{(1)}+ik\left({c}_{s}^{2}+{V_{A}}(x)^{2}\right)B% _{1}^{(1)}\\ &+2{\Omega}_{0}\left(\cos\lambda v_{z1}^{(1)}-\sin\lambda v_{y1}^{(1)}\right)-% ik\psi^{(1)}=0,\end{split}start_ROW start_CELL end_CELL start_CELL - italic_i italic_ω italic_v start_POSTSUBSCRIPT italic_x 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT + italic_i italic_k ( italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( italic_x ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL + 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( roman_cos italic_λ italic_v start_POSTSUBSCRIPT italic_z 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT - roman_sin italic_λ italic_v start_POSTSUBSCRIPT italic_y 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ) - italic_i italic_k italic_ψ start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW (21)

where VA(x)=B0(x)/μ0ρ0(x)VAsubscript𝑉𝐴𝑥subscript𝐵0𝑥subscript𝜇0subscript𝜌0𝑥subscript𝑉𝐴{V_{A}}(x)=B_{0}(x)/\sqrt{\mu_{0}\rho_{0}(x)}V_{A}italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( italic_x ) = italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) / square-root start_ARG italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT 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.

ω2=[cs2+VA2(x)]k2+4Ω02ωJ2,i.e.,ω=[cs2+VA2(x)]1/2(k2kJ2)1/2,\begin{split}&\omega^{2}=\left[{c}_{s}^{2}+{V_{A}}^{2}(x)\right]k^{2}+4{\Omega% }_{0}^{2}-\omega_{J}^{2},\\ \rm{i.e.,}~{}&\omega=\left[{c}_{s}^{2}+{V_{A}}^{2}(x)\right]^{1/2}\left(k^{2}-% k_{J}^{2}\right)^{1/2},\end{split}start_ROW start_CELL end_CELL start_CELL italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = [ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ] italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL roman_i . roman_e . , end_CELL start_CELL italic_ω = [ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT ( italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT , end_CELL end_ROW (22)

where kJsubscript𝑘𝐽k_{J}italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT is the critical Jeans wave number modified by the Coriolis force, given by,

kJ=(ωJ24Ω02cs2+VA2(x))1/2.subscript𝑘𝐽superscriptsuperscriptsubscript𝜔𝐽24superscriptsubscriptΩ02superscriptsubscript𝑐𝑠2superscriptsubscript𝑉𝐴2𝑥12k_{J}=\left(\frac{\omega_{J}^{2}-4\Omega_{0}^{2}}{c_{s}^{2}+V_{A}^{2}(x)}% \right)^{1/2}.italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = ( divide start_ARG italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) end_ARG ) start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT . (23)

From Eq. (22), we note that the fast magnetosonic wave becomes dispersive due to the presence of the term proportional to Ω02superscriptsubscriptΩ02\Omega_{0}^{2}roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, associated with the Coriolis force, and the Jeans instability may occur due to the term proportional to ωJ2superscriptsubscript𝜔𝐽2\omega_{J}^{2}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT by the influence of the self-gravitating force [33]. The instability occurs in the region k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT, provided the self-gravity force dominates over the Coriolis force, i.e., ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. In the other region, i.e., k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT, the magnetosonic wave can propagate as a real (stable) eigenmode with ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. On the other hand, in absence of the gravity effects or when the Coriolis force dominates over the self-gravity force with ωJ<2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}<2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT < 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, 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 k2superscript𝑘2k^{2}italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT via the Alfvén velocity VA(x)subscript𝑉𝐴𝑥V_{A}(x)italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( italic_x ), 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 (ω𝜔\Re\omegaroman_ℜ italic_ω, for k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT) and the growth rate of instability (ω𝜔\Im\omegaroman_ℑ italic_ω, for k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT) by the effects of the density and magnetic field inhomogeneities [See Eq. (15)], the Coriolis force, and the self-gravity force such that ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. We find that while the real wave frequency increases with the wave number k𝑘kitalic_k, the instability growth rate falls off from a nonzero value with k𝑘kitalic_k having a cut-off at the critical Jeans wave number, i.e., k=kJ𝑘subscript𝑘𝐽k=k_{J}italic_k = italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT. Such a critical value shifts towards a higher value of k𝑘kitalic_k due to the effects of the inhomogeneities in which ψ0(x)subscript𝜓0𝑥\psi_{0}(x)italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) increases, but ρ0(x)subscript𝜌0𝑥\rho_{0}(x)italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) decreases with increasing values of x>0𝑥0x>0italic_x > 0. In this case, the instability domain for k𝑘kitalic_k 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 k𝑘kitalic_k, not reported before. On the other hand, a small reduction of the angular frequency Ω0subscriptΩ0\Omega_{0}roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT of the rotating fluid, associated with the Coriolis force, can significantly reduce both the growth rate of instability and the instability domain in k𝑘kitalic_k. We do not consider the case of ωJ2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}\approx 2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT ≈ 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT at which the fast magnetosonic wave becomes dispersionless. Such a strict condition may be applicable to low-frequency long-wavelength plasma oscillations.

Refer to caption
Figure 2: The real (ω)𝜔(\Re\omega)( roman_ℜ italic_ω ) and imaginary (ω)𝜔(\Im\omega)( roman_ℑ italic_ω ) parts of the fast magnetosonic wave frequency ω𝜔\omegaitalic_ω [Eq. (22)] are shown when ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The solid, dashed, and dotted lines correspond to the parameter values (i) Ω0=0.8,ωJ=2.1Ω0,x=0.5formulae-sequencesubscriptΩ00.8formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥0.5\Omega_{0}=0.8,~{}\omega_{J}=2.1\Omega_{0},~{}x=0.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.8 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 0.5, (ii) Ω0=0.8,ωJ=2.1Ω0,x=1.5formulae-sequencesubscriptΩ00.8formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥1.5\Omega_{0}=0.8,~{}\omega_{J}=2.1\Omega_{0},~{}x=1.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.8 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 1.5, and (iii) Ω0=0.5,ωJ=2.1Ω0,x=1.5formulae-sequencesubscriptΩ00.5formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥1.5\Omega_{0}=0.5,~{}\omega_{J}=2.1\Omega_{0},~{}x=1.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 1.5 respectively. The other fixed parameter values are cs=0.5subscript𝑐𝑠0.5c_{s}=0.5italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0.5 and ψ0=0.5subscript𝜓00.5\psi_{0}=-0.5italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - 0.5.

III.2 Second-order perturbations: Compatibility condition

For the second order (n=2𝑛2n=2italic_n = 2) reduced equations with l=1𝑙1l=1italic_l = 1, we obtain the following expressions for different harmonic modes

vx1(2)=ωkB1(2)i(vgkω)k2B1(1)ξ,vy1(2)=2iΩ0sinλkB1(2)+2Ω0sinλk2B1(1)ξ,vz1(2)=2iΩ0cosλkB1(2)2Ω0cosλk2B1(1)ξ,ψ1(2)=ωJ2k2B1(2)+2iωJ2k3B1(1)ξ,formulae-sequencesuperscriptsubscriptsubscript𝑣𝑥12𝜔𝑘superscriptsubscript𝐵12𝑖subscript𝑣𝑔𝑘𝜔superscript𝑘2superscriptsubscript𝐵11𝜉formulae-sequencesuperscriptsubscriptsubscript𝑣𝑦122𝑖subscriptΩ0𝜆𝑘superscriptsubscript𝐵122subscriptΩ0𝜆superscript𝑘2superscriptsubscript𝐵11𝜉formulae-sequencesuperscriptsubscriptsubscript𝑣𝑧122𝑖subscriptΩ0𝜆𝑘superscriptsubscript𝐵122subscriptΩ0𝜆superscript𝑘2superscriptsubscript𝐵11𝜉superscriptsubscript𝜓12superscriptsubscript𝜔𝐽2superscript𝑘2superscriptsubscript𝐵122𝑖superscriptsubscript𝜔𝐽2superscript𝑘3superscriptsubscript𝐵11𝜉\begin{split}{v_{x}}_{1}^{(2)}&=\frac{\omega}{k}B_{1}^{(2)}-\frac{i(v_{g}k-% \omega)}{k^{2}}\frac{\partial B_{1}^{(1)}}{\partial\xi},\\ {v_{y}}_{1}^{(2)}&=-\frac{2i\Omega_{0}\sin\lambda}{k}B_{1}^{(2)}+\frac{2\Omega% _{0}\sin\lambda}{k^{2}}\frac{\partial B_{1}^{(1)}}{\partial\xi},\\ {v_{z}}_{1}^{(2)}&=\frac{2i\Omega_{0}\cos\lambda}{k}B_{1}^{(2)}-\frac{2\Omega_% {0}\cos\lambda}{k^{2}}\frac{\partial B_{1}^{(1)}}{\partial\xi},\\ \psi_{1}^{(2)}&=\frac{\omega_{J}^{2}}{k^{2}}B_{1}^{(2)}+\frac{2i\omega_{J}^{2}% }{k^{3}}\frac{\partial B_{1}^{(1)}}{\partial\xi},\end{split}start_ROW start_CELL italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = divide start_ARG italic_ω end_ARG start_ARG italic_k end_ARG italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT - divide start_ARG italic_i ( italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_k - italic_ω ) end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG ∂ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_ξ end_ARG , end_CELL end_ROW start_ROW start_CELL italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = - divide start_ARG 2 italic_i roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_λ end_ARG start_ARG italic_k end_ARG italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT + divide start_ARG 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_λ end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG ∂ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_ξ end_ARG , end_CELL end_ROW start_ROW start_CELL italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = divide start_ARG 2 italic_i roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_λ end_ARG start_ARG italic_k end_ARG italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT - divide start_ARG 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_λ end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG ∂ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_ξ end_ARG , end_CELL end_ROW start_ROW start_CELL italic_ψ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = divide start_ARG italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT + divide start_ARG 2 italic_i italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG divide start_ARG ∂ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_ξ end_ARG , end_CELL end_ROW (24)

together with the compatibility condition

vgωk=[cs2+VA2(x)]/vp,subscript𝑣𝑔𝜔𝑘delimited-[]superscriptsubscript𝑐𝑠2subscriptsuperscript𝑉2𝐴𝑥subscript𝑣𝑝v_{g}\equiv\frac{\partial\omega}{\partial k}={\left[c_{s}^{2}+V^{2}_{A}(x)% \right]}/{v_{p}},italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ≡ divide start_ARG ∂ italic_ω end_ARG start_ARG ∂ italic_k end_ARG = [ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( italic_x ) ] / italic_v start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT , (25)

where vp=ω/ksubscript𝑣𝑝𝜔𝑘v_{p}=\omega/kitalic_v start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT = italic_ω / italic_k is the phase velocity of the carrier magnetosonic waves. From Eq. (25), it is evident that the group velocity vgsubscript𝑣𝑔v_{g}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT 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 vgsubscript𝑣𝑔v_{g}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT on the wave frequency ω𝜔\omegaitalic_ω, the group velocity can be real (imaginary) for k>kJ(k<kJ)𝑘subscript𝑘𝐽𝑘subscript𝑘𝐽k>k_{J}~{}(k<k_{J})italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT ( italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT ) when ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. It can also be real when ωJ<2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}<2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT < 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT for any real value of k𝑘kitalic_k. The expression for vgsubscript𝑣𝑔v_{g}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT is clearly modified by the effects of the background magnetic field or density inhomogeneity. Furthermore, having known the characteristics of ω𝜔\omegaitalic_ω (See Fig. 2), one can also investigate the features of the real and imaginary parts of vgsubscript𝑣𝑔v_{g}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT. It is seen that the real part of the group velocity decreases with k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT 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 vgsubscript𝑣𝑔v_{g}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT is always negative in the domain k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT and its magnitude increases with increasing values of k𝑘kitalic_k. The detailed analysis of vgsubscript𝑣𝑔v_{g}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT 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 n=2𝑛2n=2italic_n = 2 and l=2𝑙2l=2italic_l = 2 appear due to the nonlinear self-interaction of the carrier waves, and they are found to be proportional to [B11]2superscriptdelimited-[]superscriptsubscript𝐵112[B_{1}^{1}]^{2}[ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT as

vx2(2)=ωkB2(2)ωk[B1(1)]2,vy2(2)=iΩ0sinλk23(4Ω02ωJ2)[{2cs2+3VA2(x)}k2+12Ω023ωJ2][B1(1)]2,vz2(2)=iΩ0cosλk23(4Ω02ωJ2)[{2cs2+3VA2(x)}k2+12Ω023ωJ2][B1(1)]2,ψ2(2)=ωJ24k223(4Ω02ωJ2)[{2cs2+3VA2(x)}k2+12Ω023ωJ2][B1(1)]2,B2(2)=23(4Ω02ωJ2)[{2cs2+3VA2(x)}k2+12Ω023ωJ2][B1(1)]2.formulae-sequencesuperscriptsubscriptsubscript𝑣𝑥22𝜔𝑘superscriptsubscript𝐵22𝜔𝑘superscriptdelimited-[]superscriptsubscript𝐵112formulae-sequencesuperscriptsubscriptsubscript𝑣𝑦22𝑖subscriptΩ0𝜆𝑘234superscriptsubscriptΩ02superscriptsubscript𝜔𝐽2delimited-[]2superscriptsubscript𝑐𝑠23superscriptsubscript𝑉𝐴2𝑥superscript𝑘212superscriptsubscriptΩ023superscriptsubscript𝜔𝐽2superscriptdelimited-[]superscriptsubscript𝐵112formulae-sequencesuperscriptsubscriptsubscript𝑣𝑧22𝑖subscriptΩ0𝜆𝑘234superscriptsubscriptΩ02superscriptsubscript𝜔𝐽2delimited-[]2superscriptsubscript𝑐𝑠23superscriptsubscript𝑉𝐴2𝑥superscript𝑘212superscriptsubscriptΩ023superscriptsubscript𝜔𝐽2superscriptdelimited-[]superscriptsubscript𝐵112formulae-sequencesuperscriptsubscript𝜓22superscriptsubscript𝜔𝐽24superscript𝑘2234superscriptsubscriptΩ02superscriptsubscript𝜔𝐽2delimited-[]2superscriptsubscript𝑐𝑠23superscriptsubscript𝑉𝐴2𝑥superscript𝑘212superscriptsubscriptΩ023superscriptsubscript𝜔𝐽2superscriptdelimited-[]superscriptsubscript𝐵112superscriptsubscript𝐵22234superscriptsubscriptΩ02superscriptsubscript𝜔𝐽2delimited-[]2superscriptsubscript𝑐𝑠23superscriptsubscript𝑉𝐴2𝑥superscript𝑘212superscriptsubscriptΩ023superscriptsubscript𝜔𝐽2superscriptdelimited-[]superscriptsubscript𝐵112\begin{split}{v_{x}}_{2}^{(2)}&=\frac{\omega}{k}B_{2}^{(2)}-\frac{\omega}{k}% \left[B_{1}^{(1)}\right]^{2},\\ {v_{y}}_{2}^{(2)}&=-\frac{i\Omega_{0}\sin\lambda}{k}\frac{2}{3(4\Omega_{0}^{2}% -\omega_{J}^{2})}\left[\{2c_{s}^{2}+3V_{A}^{2}(x)\}k^{2}+12\Omega_{0}^{2}-3% \omega_{J}^{2}\right]\left[B_{1}^{(1)}\right]^{2},\\ {v_{z}}_{2}^{(2)}&=\frac{i\Omega_{0}\cos\lambda}{k}\frac{2}{3(4\Omega_{0}^{2}-% \omega_{J}^{2})}\left[\{2c_{s}^{2}+3V_{A}^{2}(x)\}k^{2}+12\Omega_{0}^{2}-3% \omega_{J}^{2}\right]\left[B_{1}^{(1)}\right]^{2},\\ \psi_{2}^{(2)}&=\frac{\omega_{J}^{2}}{4k^{2}}\frac{2}{3(4\Omega_{0}^{2}-\omega% _{J}^{2})}\left[\{2c_{s}^{2}+3V_{A}^{2}(x)\}k^{2}+12\Omega_{0}^{2}-3\omega_{J}% ^{2}\right]\left[B_{1}^{(1)}\right]^{2},\\ B_{2}^{(2)}&=\frac{2}{3(4\Omega_{0}^{2}-\omega_{J}^{2})}\left[\{2c_{s}^{2}+3V_% {A}^{2}(x)\}k^{2}+12\Omega_{0}^{2}-3\omega_{J}^{2}\right]\left[B_{1}^{(1)}% \right]^{2}.\end{split}start_ROW start_CELL italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = divide start_ARG italic_ω end_ARG start_ARG italic_k end_ARG italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT - divide start_ARG italic_ω end_ARG start_ARG italic_k end_ARG [ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = - divide start_ARG italic_i roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_λ end_ARG start_ARG italic_k end_ARG divide start_ARG 2 end_ARG start_ARG 3 ( 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG [ { 2 italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 3 italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) } italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 12 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] [ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = divide start_ARG italic_i roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_λ end_ARG start_ARG italic_k end_ARG divide start_ARG 2 end_ARG start_ARG 3 ( 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG [ { 2 italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 3 italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) } italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 12 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] [ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_ψ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = divide start_ARG italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG 2 end_ARG start_ARG 3 ( 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG [ { 2 italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 3 italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) } italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 12 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] [ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_B start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = divide start_ARG 2 end_ARG start_ARG 3 ( 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG [ { 2 italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 3 italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) } italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 12 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 3 italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] [ italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . end_CELL end_ROW (26)

The nonlinear self-interactions of the first-order carrier wave modes also result into the generation of the zeroth harmonic modes. Thus, for n=2𝑛2n=2italic_n = 2, l=0𝑙0l=0italic_l = 0 we obtain

B0(2)=vy0(2)=vz0(2)=0,vx0(2)=2ωk|B1(1)|2,ψ0(2)=(2ω2+k24Ω02)k2|B1(1)|2.formulae-sequencesuperscriptsubscript𝐵02superscriptsubscriptsubscript𝑣𝑦02superscriptsubscriptsubscript𝑣𝑧020formulae-sequencesuperscriptsubscriptsubscript𝑣𝑥022𝜔𝑘superscriptsuperscriptsubscript𝐵112superscriptsubscript𝜓022superscript𝜔2superscript𝑘24superscriptsubscriptΩ02superscript𝑘2superscriptsuperscriptsubscript𝐵112\begin{split}B_{0}^{(2)}&={v_{y}}_{0}^{(2)}={v_{z}}_{0}^{(2)}=0,\\ {v_{x}}_{0}^{(2)}&=-\frac{2\omega}{k}|B_{1}^{(1)}|^{2},\\ \psi_{0}^{(2)}&=\frac{(2\omega^{2}+k^{2}-4\Omega_{0}^{2})}{k^{2}}|B_{1}^{(1)}|% ^{2}.\end{split}start_ROW start_CELL italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT = italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW start_ROW start_CELL italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = - divide start_ARG 2 italic_ω end_ARG start_ARG italic_k end_ARG | italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT end_CELL start_CELL = divide start_ARG ( 2 italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG | italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . end_CELL end_ROW (27)

Finally, substituting all the above derived expressions into the third-order harmonic modes (n=3𝑛3n=3italic_n = 3 and l=1𝑙1l=1italic_l = 1), we obtain the following NLS equation

iBτ+P2Bξ2+Q|B|2B=0,𝑖𝐵𝜏𝑃superscript2𝐵superscript𝜉2𝑄superscript𝐵2𝐵0\displaystyle i\frac{\partial B}{\partial\tau}+P\frac{\partial^{2}B}{\partial% \xi^{2}}+Q|B|^{2}B=0,italic_i divide start_ARG ∂ italic_B end_ARG start_ARG ∂ italic_τ end_ARG + italic_P divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B end_ARG start_ARG ∂ italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + italic_Q | italic_B | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B = 0 , (28)

where we have replaced B1(1)superscriptsubscript𝐵11B_{1}^{(1)}italic_B start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT by B𝐵Bitalic_B for simplicity. The group velocity dispersion coefficient P𝑃Pitalic_P and the nonlinear coefficient Q𝑄Qitalic_Q, appeared due to the carrier wave self-interactions, are given by

P122ωk2=12(cs2+VA2(x))ω3(4Ω02ωJ2)=124Ω02ωJ2[cs2+VA2(x)]1/2(k2kJ2)3/2,𝑃12superscript2𝜔superscript𝑘212superscriptsubscript𝑐𝑠2superscriptsubscript𝑉𝐴2𝑥superscript𝜔34superscriptsubscriptΩ02superscriptsubscript𝜔𝐽2124superscriptsubscriptΩ02superscriptsubscript𝜔𝐽2superscriptdelimited-[]superscriptsubscript𝑐𝑠2superscriptsubscript𝑉𝐴2𝑥12superscriptsuperscript𝑘2superscriptsubscript𝑘𝐽232\begin{split}P\equiv\frac{1}{2}\frac{\partial^{2}\omega}{\partial k^{2}}&=% \frac{1}{2}\frac{\left(c_{s}^{2}+V_{A}^{2}(x)\right)}{\omega^{3}}\left(4\Omega% _{0}^{2}-{\omega}_{J}^{2}\right)\\ &=\frac{1}{2}\frac{4\Omega_{0}^{2}-\omega_{J}^{2}}{\left[c_{s}^{2}+V_{A}^{2}(x% )\right]^{1/2}\left(k^{2}-k_{J}^{2}\right)^{3/2}},\end{split}start_ROW start_CELL italic_P ≡ divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω end_ARG start_ARG ∂ italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG ( italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ) end_ARG start_ARG italic_ω start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ( 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL = divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG [ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT ( italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT end_ARG , end_CELL end_ROW (29)
Q=12ω[ω223(3ω2cs2k2)(4Ω02ωJ2)×{(2ω2+k2)12(4Ω02ωJ2)}].𝑄12𝜔delimited-[]superscript𝜔2233superscript𝜔2superscriptsubscript𝑐𝑠2superscript𝑘24superscriptsubscriptΩ02superscriptsubscript𝜔𝐽22superscript𝜔2superscript𝑘2124superscriptsubscriptΩ02superscriptsubscript𝜔𝐽2\begin{split}Q=&\frac{1}{2\omega}\left[\omega^{2}-\frac{2}{3}\frac{\left(3% \omega^{2}-c_{s}^{2}k^{2}\right)}{(4\Omega_{0}^{2}-{\omega}_{J}^{2})}\right.\\ &\left.\times\left\{(2\omega^{2}+k^{2})-\frac{1}{2}\left(4\Omega_{0}^{2}-{% \omega}_{J}^{2}\right)\right\}\right].\end{split}start_ROW start_CELL italic_Q = end_CELL start_CELL divide start_ARG 1 end_ARG start_ARG 2 italic_ω end_ARG [ italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG 2 end_ARG start_ARG 3 end_ARG divide start_ARG ( 3 italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG start_ARG ( 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × { ( 2 italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) } ] . end_CELL end_ROW (30)

From the dispersion Eq. (22), we note that the wave frequency ω𝜔\omegaitalic_ω can be either real or purely imaginary depending on whether k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT or k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT. So, the coefficients P𝑃Pitalic_P and Q𝑄Qitalic_Q can also be either real or purely imaginary. While the coefficient P𝑃Pitalic_P 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 Q𝑄Qitalic_Q 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 Q𝑄Qitalic_Q appears due to transverse velocity perturbations and the second term of Q𝑄Qitalic_Q 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: ωJ<2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}<2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT < 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, for-all\forall k𝑘kitalic_k.

  • Case II: ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT.

  • Case III: ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT.

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

B=η1/2exp(iaξσ2P𝑑ξ),𝐵superscript𝜂12𝑖superscriptsubscript𝑎𝜉𝜎2𝑃differential-d𝜉B=\eta^{1/2}\exp\left(i\int_{a}^{\xi}\frac{\sigma}{2P}d\xi\right),italic_B = italic_η start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT roman_exp ( italic_i ∫ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ξ end_POSTSUPERSCRIPT divide start_ARG italic_σ end_ARG start_ARG 2 italic_P end_ARG italic_d italic_ξ ) , (31)

where a𝑎aitalic_a is some constant, and η𝜂\etaitalic_η and σ𝜎\sigmaitalic_σ are real functions of ξ𝜉\xiitalic_ξ and τ𝜏\tauitalic_τ. Substituting Eq. (31) into Eq. (28), and separating the real and imaginary parts, we obtain

ητ+(ησ)ξ=0,𝜂𝜏𝜂𝜎𝜉0\frac{\partial\eta}{\partial\tau}+\frac{\partial(\eta\sigma)}{\partial\xi}=0,divide start_ARG ∂ italic_η end_ARG start_ARG ∂ italic_τ end_ARG + divide start_ARG ∂ ( italic_η italic_σ ) end_ARG start_ARG ∂ italic_ξ end_ARG = 0 , (32)
στ+σηξ=2PQηξ+P2ξ[η1/2ξ(η1/2ηξ)].𝜎𝜏𝜎𝜂𝜉2𝑃𝑄𝜂𝜉superscript𝑃2𝜉delimited-[]superscript𝜂12𝜉superscript𝜂12𝜂𝜉\begin{split}\frac{\partial\sigma}{\partial\tau}+\sigma\frac{\partial\eta}{% \partial\xi}=&2PQ\frac{\partial\eta}{\partial\xi}\\ &+P^{2}\frac{\partial}{\partial\xi}\left[\eta^{-1/2}\frac{\partial}{\partial% \xi}\left(\eta^{-1/2}\frac{\partial\eta}{\partial\xi}\right)\right].\end{split}start_ROW start_CELL divide start_ARG ∂ italic_σ end_ARG start_ARG ∂ italic_τ end_ARG + italic_σ divide start_ARG ∂ italic_η end_ARG start_ARG ∂ italic_ξ end_ARG = end_CELL start_CELL 2 italic_P italic_Q divide start_ARG ∂ italic_η end_ARG start_ARG ∂ italic_ξ end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL + italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG [ italic_η start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG ( italic_η start_POSTSUPERSCRIPT - 1 / 2 end_POSTSUPERSCRIPT divide start_ARG ∂ italic_η end_ARG start_ARG ∂ italic_ξ end_ARG ) ] . end_CELL end_ROW (33)

Next, we modulate the wave amplitude and phase by small plane-wave perturbations with the wave number K𝐾Kitalic_K and the wave frequency ΩΩ\Omegaroman_Ω as

η=η0+η1cos(KξΩτ)+η2sin(KξΩτ),σ=σ1cos(KξΩτ)+σ2sin(KξΩτ),formulae-sequence𝜂subscript𝜂0subscript𝜂1𝐾𝜉Ω𝜏subscript𝜂2𝐾𝜉Ω𝜏𝜎subscript𝜎1𝐾𝜉Ω𝜏subscript𝜎2𝐾𝜉Ω𝜏\begin{split}\eta&=\eta_{0}+\eta_{1}\cos(K\xi-\Omega\tau)+\eta_{2}\sin(K\xi-% \Omega\tau),\\ \sigma&=\sigma_{1}\cos(K\xi-\Omega\tau)+\sigma_{2}\sin(K\xi-\Omega\tau),\end{split}start_ROW start_CELL italic_η end_CELL start_CELL = italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_η start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_cos ( italic_K italic_ξ - roman_Ω italic_τ ) + italic_η start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_sin ( italic_K italic_ξ - roman_Ω italic_τ ) , end_CELL end_ROW start_ROW start_CELL italic_σ end_CELL start_CELL = italic_σ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_cos ( italic_K italic_ξ - roman_Ω italic_τ ) + italic_σ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT roman_sin ( italic_K italic_ξ - roman_Ω italic_τ ) , end_CELL end_ROW (34)

where η0subscript𝜂0\eta_{0}italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is a constant, and ηjsubscript𝜂𝑗\eta_{j}italic_η start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT and σjsubscript𝜎𝑗\sigma_{j}italic_σ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT for j=1,2𝑗12j=1,2italic_j = 1 , 2, are the amplitudes of perturbations.

Substituting the perturbation expansion (34) in Eqs. (32) and (33), and looking for nonzero solutions of η1subscript𝜂1\eta_{1}italic_η start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and η2subscript𝜂2\eta_{2}italic_η start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, we obtain the following dispersion relation for the perturbed wave of modulation in self-gravitating rotating magnetoplasmas

Ω2=P2K2(K22η0QP).superscriptΩ2superscript𝑃2superscript𝐾2superscript𝐾22subscript𝜂0𝑄𝑃\Omega^{2}=P^{2}K^{2}\left(K^{2}-\frac{2\eta_{0}Q}{P}\right).roman_Ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_P start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG 2 italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_Q end_ARG start_ARG italic_P end_ARG ) . (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: ωJ<2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}<2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT < 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, for-all\forall k𝑘kitalic_k

We consider the case when the contribution from the Coriolis force dominates over the self-gravitating force, i.e. when ωJ<2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}<2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT < 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. In this case, the carrier Jeans wave frequency ω𝜔\omegaitalic_ω becomes real for any values of the wave number k𝑘kitalic_k [See Eq. (22)], and so are P𝑃Pitalic_P and Q𝑄Qitalic_Q. Also, from Eq. (29), it is evident that that P>0𝑃0P>0italic_P > 0 for ωJ<2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}<2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT < 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Furthermore, inspecting on the expression of Q𝑄Qitalic_Q [Eq. (30)], we find that when ωJ<2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}<2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT < 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT holds, the second term of Q𝑄Qitalic_Q becomes positive and larger than the first term, giving Q<0𝑄0Q<0italic_Q < 0. Thus, in this case, PQ<0𝑃𝑄0PQ<0italic_P italic_Q < 0, and Eq. (35) gives a real wave frequency ΩΩ\Omegaroman_Ω 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: ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT

We consider the case when the self-gravitating force dominates over the Coriolis force, i.e. ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT. In this case, as discussed before in Sec. III.1, since the carrier wave frequency ω𝜔\omegaitalic_ω is real, the coefficients P𝑃Pitalic_P and Q𝑄Qitalic_Q of the NLS equation are also real. However, in contrast to Case I, P<0𝑃0P<0italic_P < 0 and Q𝑄Qitalic_Q can be both positive and negative for k(>kJ)annotated𝑘absentsubscript𝑘𝐽k~{}(>k_{J})italic_k ( > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT ), so that one can have both stable (PQ<0)𝑃𝑄0(PQ<0)( italic_P italic_Q < 0 ) and unstable (PQ>0)𝑃𝑄0(PQ>0)( italic_P italic_Q > 0 ) regions. Figure 3 shows the plot of Q𝑄Qitalic_Q versus the carrier wave number k𝑘kitalic_k. We find that due to the inhomogenity effects, the domain of k𝑘kitalic_k in which Q<0𝑄0Q<0italic_Q < 0 (for which PQ>0𝑃𝑄0PQ>0italic_P italic_Q > 0 and the modulational instability occurs) expands and the domain shifts towards lower values of k𝑘kitalic_k due to a reduction of the magnitude of Ω0subscriptΩ0\Omega_{0}roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. For example, for (i) Ω0=0.8,ωJ=2.1Ω0,x=0.5formulae-sequencesubscriptΩ00.8formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥0.5\Omega_{0}=0.8,~{}\omega_{J}=2.1\Omega_{0},~{}x=0.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.8 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 0.5, we have Q1<0(>0)subscript𝑄1annotated0absent0Q_{1}<0~{}(>0)italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0 ( > 0 ) in 0.6k<0.624(k0.624)less-than-or-similar-to0.6𝑘0.624greater-than-or-equivalent-to𝑘0.6240.6\lesssim k<0.624~{}(k\gtrsim 0.624)0.6 ≲ italic_k < 0.624 ( italic_k ≳ 0.624 ) (See the solid line of Fig. 3) (ii) Ω0=0.8,ωJ=2.1Ω0,x=1.5formulae-sequencesubscriptΩ00.8formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥1.5\Omega_{0}=0.8,~{}\omega_{J}=2.1\Omega_{0},~{}x=1.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.8 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 1.5, we have Q1<0(>0)subscript𝑄1annotated0absent0Q_{1}<0~{}(>0)italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0 ( > 0 ) in 0.85k<0.94(k0.94)less-than-or-similar-to0.85𝑘0.94greater-than-or-equivalent-to𝑘0.940.85\lesssim k<0.94~{}(k\gtrsim 0.94)0.85 ≲ italic_k < 0.94 ( italic_k ≳ 0.94 ) (See the dashed line of Fig. 3), and (iii) Ω0=0.5,ωJ=2.1Ω0,x=1.5formulae-sequencesubscriptΩ00.5formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥1.5\Omega_{0}=0.5,~{}\omega_{J}=2.1\Omega_{0},~{}x=1.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 1.5, we have Q1<0(>0)subscript𝑄1annotated0absent0Q_{1}<0~{}(>0)italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0 ( > 0 ) in 0.39k<0.41(k0.41)less-than-or-similar-to0.39𝑘0.41greater-than-or-equivalent-to𝑘0.410.39\lesssim k<0.41~{}(k\gtrsim 0.41)0.39 ≲ italic_k < 0.41 ( italic_k ≳ 0.41 ) (See the dotted line of Fig. 3). The other fixed parameter values are cs=0.5subscript𝑐𝑠0.5c_{s}=0.5italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0.5 and ψ0=0.5subscript𝜓00.5\psi_{0}=-0.5italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - 0.5.

In the case of PQ>0𝑃𝑄0PQ>0italic_P italic_Q > 0 for which the modulational instability occurs, the instability growth rate (Γ)Γ(\Gamma)( roman_Γ ) can be obtained from Eq. (35) for K<Kc𝐾subscript𝐾𝑐K<K_{c}italic_K < italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT as

Γ=|P|KKc2K2,Γ𝑃𝐾superscriptsubscript𝐾𝑐2superscript𝐾2\Gamma=|P|K\sqrt{K_{c}^{2}-K^{2}},roman_Γ = | italic_P | italic_K square-root start_ARG italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (36)

where Kc=2η0Q/Psubscript𝐾𝑐2subscript𝜂0𝑄𝑃K_{c}=\sqrt{2\eta_{0}Q/P}italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = square-root start_ARG 2 italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_Q / italic_P end_ARG is the critical wave number of perturbation, and the maximum growth rate, Γmax=η0|Q|subscriptΓsubscript𝜂0𝑄\Gamma_{\max}=\eta_{0}|Q|roman_Γ start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | italic_Q | is attained at K=Kc/2𝐾subscript𝐾𝑐2K=K_{c}/\sqrt{2}italic_K = italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT / square-root start_ARG 2 end_ARG.

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 Ω0subscriptΩ0\Omega_{0}roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT compared to the Jeans frequency.

Refer to caption
Figure 3: Plot of Q𝑄Qitalic_Q [Eq. (30)] vs k𝑘kitalic_k is shown for Case II: ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT. The fixed parameter values and different parameter values for the solid, dashed, and dotted lines are the same as for Fig. 2. Since P<0𝑃0P<0italic_P < 0 kfor-all𝑘\forall k∀ italic_k, the stable and unstable regions in k𝑘kitalic_k are corresponding to Q>0𝑄0Q>0italic_Q > 0 and Q<0𝑄0Q<0italic_Q < 0 respectively.
Refer to caption
Figure 4: The modulational instability (PQ>0)𝑃𝑄0(PQ>0)( italic_P italic_Q > 0 ) growth rate ΓΓ\Gammaroman_Γ [Eq. (36)] is shown against the modulation wave number K𝐾Kitalic_K of perturbation. The fixed parameter values and different parameter values for the solid, dashed, and dotted lines are the same as for Fig. 2.

IV.3 Case III: ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT

Here, we consider the most interesting case when ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, but k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT. From the analysis in Sec. III.1, it is evident that the carrier wave frequency ω𝜔\omegaitalic_ω is purely imaginary (Jeans instability) and so are the coefficients P𝑃Pitalic_P and Q𝑄Qitalic_Q of the NLS equation. Thus, assuming P=iP1𝑃𝑖subscript𝑃1P=iP_{1}italic_P = italic_i italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and Q=iQ1𝑄𝑖subscript𝑄1Q=iQ_{1}italic_Q = italic_i italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, Eq. (35) gives

Ω=P1K2η0Q1P1K2,Ωsubscript𝑃1𝐾2subscript𝜂0subscript𝑄1subscript𝑃1superscript𝐾2\begin{split}\Omega=P_{1}K\sqrt{\frac{2\eta_{0}Q_{1}}{P_{1}}-K^{2}},\end{split}start_ROW start_CELL roman_Ω = italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_K square-root start_ARG divide start_ARG 2 italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG - italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , end_CELL end_ROW (37)

where P1subscript𝑃1P_{1}italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and Q1subscript𝑄1Q_{1}italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT are reals, given by,

P1=12ωJ24Ω02[cs2+VA2(x)]1/2(kJ2k2)3/2,Q1=12γ[γ2+133γ2+cs2k2ωJ24Ω02×(4γ22k2(ωJ24Ω02))],formulae-sequencesubscript𝑃112superscriptsubscript𝜔𝐽24superscriptsubscriptΩ02superscriptdelimited-[]superscriptsubscript𝑐𝑠2superscriptsubscript𝑉𝐴2𝑥12superscriptsuperscriptsubscript𝑘𝐽2superscript𝑘232subscript𝑄112𝛾delimited-[]superscript𝛾2133superscript𝛾2superscriptsubscript𝑐𝑠2superscript𝑘2superscriptsubscript𝜔𝐽24superscriptsubscriptΩ024superscript𝛾22superscript𝑘2superscriptsubscript𝜔𝐽24superscriptsubscriptΩ02\begin{split}&P_{1}=\frac{1}{2}\frac{\omega_{J}^{2}-4\Omega_{0}^{2}}{\left[c_{% s}^{2}+V_{A}^{2}(x)\right]^{1/2}\left(k_{J}^{2}-k^{2}\right)^{3/2}},\\ &Q_{1}=-\frac{1}{2\gamma}\left[-\gamma^{2}+\frac{1}{3}\frac{3\gamma^{2}+c_{s}^% {2}k^{2}}{\omega_{J}^{2}-4\Omega_{0}^{2}}\right.\\ &\times\left.\left(4\gamma^{2}-2k^{2}-\left(\omega_{J}^{2}-4\Omega_{0}^{2}% \right)\right)\right],\end{split}start_ROW start_CELL end_CELL start_CELL italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG [ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT ( italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT end_ARG , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG 2 italic_γ end_ARG [ - italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 1 end_ARG start_ARG 3 end_ARG divide start_ARG 3 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL × ( 4 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ) ] , end_CELL end_ROW (38)

with

γω=[cs2+VA2(x)](kJ2k2)=[cs2+VA2(x)]k2+ωJ24Ω02.𝛾𝜔delimited-[]superscriptsubscript𝑐𝑠2superscriptsubscript𝑉𝐴2𝑥superscriptsubscript𝑘𝐽2superscript𝑘2delimited-[]superscriptsubscript𝑐𝑠2superscriptsubscript𝑉𝐴2𝑥superscript𝑘2superscriptsubscript𝜔𝐽24superscriptsubscriptΩ02\begin{split}\gamma\equiv\Im\omega&=\sqrt{\left[c_{s}^{2}+V_{A}^{2}(x)\right]% \left(k_{J}^{2}-k^{2}\right)}\\ &=\sqrt{\left[c_{s}^{2}+V_{A}^{2}(x)\right]k^{2}+\omega_{J}^{2}-4\Omega_{0}^{2% }}.\end{split}start_ROW start_CELL italic_γ ≡ roman_ℑ italic_ω end_CELL start_CELL = square-root start_ARG [ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ] ( italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL = square-root start_ARG [ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ] italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . end_CELL end_ROW (39)

From Eq. (37), we find that the instability condition depends on the sign of the product P1Q1subscript𝑃1subscript𝑄1P_{1}Q_{1}italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and/or on the value of K𝐾Kitalic_K smaller or larger than a critical value Kcsubscript𝐾𝑐K_{c}italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT. When P1Q1<0subscript𝑃1subscript𝑄10P_{1}Q_{1}<0italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0, Eq. (37) gives ΩΩ\Omegaroman_Ω purely imaginary for all values of K𝐾Kitalic_K and so we have the modulational instability for all K𝐾Kitalic_K, but with k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT. The corresponding growth rate of instability can be obtained by setting Ω=iΓ1Ω𝑖subscriptΓ1\Omega=i\Gamma_{1}roman_Ω = italic_i roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT as

Γ1=|P1|K2η0|Q1P1|+K2.subscriptΓ1subscript𝑃1𝐾2subscript𝜂0subscript𝑄1subscript𝑃1superscript𝐾2\Gamma_{1}=\big{|}P_{1}\big{|}K\sqrt{2\eta_{0}\Big{|}\frac{Q_{1}}{P_{1}}\Big{|% }+K^{2}}.roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = | italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | italic_K square-root start_ARG 2 italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | divide start_ARG italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG | + italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (40)

On the other hand, when P1Q1>0subscript𝑃1subscript𝑄10P_{1}Q_{1}>0italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > 0, the modulational instability occurs for K>Kc𝐾subscript𝐾𝑐K>K_{c}italic_K > italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, where Kcsubscript𝐾𝑐K_{c}italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is the critical wave number, given by,

Kc=2η0|Q1P1|.subscript𝐾𝑐2subscript𝜂0subscript𝑄1subscript𝑃1K_{c}=\sqrt{2\eta_{0}\Big{|}\frac{Q_{1}}{P_{1}}\Big{|}}.italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = square-root start_ARG 2 italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | divide start_ARG italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG | end_ARG . (41)

The corresponding instability growth rate is given by

Γ2=|P1|KK22η0|Q1P1|.subscriptΓ2subscript𝑃1𝐾superscript𝐾22subscript𝜂0subscript𝑄1subscript𝑃1\Gamma_{2}=|P_{1}|K\sqrt{K^{2}-2\eta_{0}\Big{|}\frac{Q_{1}}{P_{1}}\Big{|}}.roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = | italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT | italic_K square-root start_ARG italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 2 italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | divide start_ARG italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG | end_ARG . (42)

However, the modulated wave is stable for K<Kc𝐾subscript𝐾𝑐K<K_{c}italic_K < italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT. The instability condition stated above is in contrast to the typical condition (PQ>0𝑃𝑄0PQ>0italic_P italic_Q > 0, K<Kc𝐾subscript𝐾𝑐K<K_{c}italic_K < italic_K start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT) of modulational instability of wave envelopes in plasmas without self-gravity effects. Inspecting on the coefficients P(=iP1)annotated𝑃absent𝑖subscript𝑃1P~{}(=iP_{1})italic_P ( = italic_i italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) and Q(=iQ1)annotated𝑄absent𝑖subscript𝑄1Q~{}(=iQ_{1})italic_Q ( = italic_i italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ), we find that P1subscript𝑃1P_{1}italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is always positive for k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT and ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. From Fig. 5, we find that Q1subscript𝑄1Q_{1}italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT can be both positive and negative in a finite domain of k(<kJ)annotated𝑘absentsubscript𝑘𝐽k~{}(<k_{J})italic_k ( < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT ). The domain for Q1>0subscript𝑄10Q_{1}>0italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > 0 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) Ω0=0.8,ωJ=2.1Ω0,x=0.5formulae-sequencesubscriptΩ00.8formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥0.5\Omega_{0}=0.8,~{}\omega_{J}=2.1\Omega_{0},~{}x=0.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.8 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 0.5, we have Q1<0(>0)subscript𝑄1annotated0absent0Q_{1}<0~{}(>0)italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0 ( > 0 ) in 0<k<0.33(0.33k0.594)0𝑘0.33less-than-or-similar-to0.33𝑘less-than-or-similar-to0.5940<k<0.33~{}(0.33\lesssim k\lesssim 0.594)0 < italic_k < 0.33 ( 0.33 ≲ italic_k ≲ 0.594 ) (See the solid line of Fig. 5) (ii) Ω0=0.8,ωJ=2.1Ω0,x=1.5formulae-sequencesubscriptΩ00.8formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥1.5\Omega_{0}=0.8,~{}\omega_{J}=2.1\Omega_{0},~{}x=1.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.8 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 1.5, we have Q1<0(>0)subscript𝑄1annotated0absent0Q_{1}<0~{}(>0)italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0 ( > 0 ) in 0<k<0.4(0.4k0.842)0𝑘0.4less-than-or-similar-to0.4𝑘less-than-or-similar-to0.8420<k<0.4~{}(0.4\lesssim k\lesssim 0.842)0 < italic_k < 0.4 ( 0.4 ≲ italic_k ≲ 0.842 ) (See the dashed line of Fig. 5), and (iii) Ω0=0.5,ωJ=2.1Ω0,x=1.5formulae-sequencesubscriptΩ00.5formulae-sequencesubscript𝜔𝐽2.1subscriptΩ0𝑥1.5\Omega_{0}=0.5,~{}\omega_{J}=2.1\Omega_{0},~{}x=1.5roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.5 , italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT = 2.1 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_x = 1.5, we have Q1<0(>0)subscript𝑄1annotated0absent0Q_{1}<0~{}(>0)italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0 ( > 0 ) in 0<k<0.22(0.22k0.386)0𝑘0.22less-than-or-similar-to0.22𝑘less-than-or-similar-to0.3860<k<0.22~{}(0.22\lesssim k\lesssim 0.386)0 < italic_k < 0.22 ( 0.22 ≲ italic_k ≲ 0.386 ) (See the dotted line of Fig. 5). The other fixed parameter values are cs=0.5subscript𝑐𝑠0.5c_{s}=0.5italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0.5 and ψ0=0.5subscript𝜓00.5\psi_{0}=-0.5italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - 0.5.

Refer to caption
Figure 5: Plot of Q1subscript𝑄1Q_{1}italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT versus k𝑘kitalic_k is shown for Case III: ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT. The fixed parameter values and different parameter values for the solid, dashed, and dotted lines are the same as for Fig. 2.
Refer to caption
Figure 6: The instability growth rates (a) Γ1subscriptΓ1\Gamma_{1}roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT [Eq. (40), when P1Q1<0subscript𝑃1subscript𝑄10P_{1}Q_{1}<0italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0] Γ2subscriptΓ2\Gamma_{2}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT [Eq. (42), when P1Q1>0subscript𝑃1subscript𝑄10P_{1}Q_{1}>0italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > 0] are plotted against the modulation wave number K𝐾Kitalic_K of perturbation. The fixed parameter values and different parameter values the solid, dashed, and dotted lines are the same as for Fig. 2.

Figure 6 displays the growth rates Γ1subscriptΓ1\Gamma_{1}roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and Γ2subscriptΓ2\Gamma_{2}roman_Γ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT [Eqs. (40) and (42)] corresponding to P1Q1<0subscript𝑃1subscript𝑄10P_{1}Q_{1}<0italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 0 and P2Q2>0subscript𝑃2subscript𝑄20P_{2}Q_{2}>0italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT > 0 respectively. Evidently, we have a growing modulational instability without any cut-off at any value of K𝐾Kitalic_K. The growth rate becomes higher the larger is the Jeans frequency compared to the angular velocity.

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 400400400400 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-β𝛽\betaitalic_β plasma. The speed was, however, reduced to 190190190190 km/s and 20202020 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, n109cm3similar-to𝑛superscript109superscriptcm3n\sim 10^{9}~{}\rm{cm}^{-3}italic_n ∼ 10 start_POSTSUPERSCRIPT 9 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT, (ii) the electron temperature, Te103similar-tosubscript𝑇𝑒superscript103T_{e}\sim 10^{-3}italic_T start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT ∼ 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT K, and (iii) the plasma β0.2(β0)similar-to𝛽0.2𝛽0\beta\sim 0.2~{}(\beta\rightarrow 0)italic_β ∼ 0.2 ( italic_β → 0 ) at the magnetic field strength, B04similar-tosubscript𝐵04B_{0}\sim 4italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∼ 4 G  (40404040 G). For illustration purpose, we calculate the maximum growth rate in the case of ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT (Case II in Sec. IV.2) at B0=4subscript𝐵04B_{0}=4italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 4 G as Γmax0.0076similar-tosubscriptΓmax0.0076\Gamma_{\rm{max}}\sim 0.0076roman_Γ start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT ∼ 0.0076 (for η0=0.01subscript𝜂00.01\eta_{0}=0.01italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0.01) such that the typical growth time can be estimated as τgrowth4similar-tosubscript𝜏growth4\tau_{\rm{growth}}\sim 4italic_τ start_POSTSUBSCRIPT roman_growth end_POSTSUBSCRIPT ∼ 4 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 λJ2×104similar-tosubscript𝜆𝐽2superscript104\lambda_{J}\sim 2\times 10^{4}italic_λ start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT ∼ 2 × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT m and VA3×105similar-tosubscript𝑉𝐴3superscript105V_{A}\sim 3\times 10^{5}italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ∼ 3 × 10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT m/s. We can also estimate the group velocity as vg5×105similar-tosubscript𝑣𝑔5superscript105v_{g}\sim 5\times 10^{5}italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ∼ 5 × 10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT m/s. These velocities are close to the observational velocity 190190190190 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 VA(x)subscript𝑉𝐴𝑥V_{A}(x)italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT ( italic_x ) gets reduced. Note that the above results are valid for low-β𝛽\betaitalic_β 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 ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT (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, 𝐠=ψ𝐠𝜓\mathbf{g}=-\nabla\psibold_g = - ∇ italic_ψ is to be replaced by 𝐠=(g,0,0)𝐠𝑔00\mathbf{g}=(-g,0,0)bold_g = ( - italic_g , 0 , 0 ) 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 ωJ2superscriptsubscript𝜔𝐽2-\omega_{J}^{2}- italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT in the dispersion equation (22) for the fast magnetosonic carrier waves will be replaced by igk𝑖𝑔𝑘igkitalic_i italic_g italic_k, which will result into the real wave frequency, given by,

ωr2=[cs2+VA2(x)]k2+4Ω02superscriptsubscript𝜔𝑟2delimited-[]superscriptsubscript𝑐𝑠2superscriptsubscript𝑉𝐴2𝑥superscript𝑘24superscriptsubscriptΩ02\omega_{r}^{2}=\left[{c}_{s}^{2}+{V_{A}}^{2}(x)\right]k^{2}+4{\Omega}_{0}^{2}italic_ω start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = [ italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ] italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (43)

and the instability growth rate, γ=gk/ωr𝛾𝑔𝑘subscript𝜔𝑟\gamma=gk/\omega_{r}italic_γ = italic_g italic_k / italic_ω start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT. From the inhomogeneous equilibrium state one can obtain by the same assumption as for the fluid density in Sec. II.1 as ρ0(x)ρ0+gxsubscript𝜌0𝑥subscript𝜌0𝑔𝑥\rho_{0}(x)\approx\rho_{0}+gxitalic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) ≈ italic_ρ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_g italic_x. The coefficients of the group velocity dispersion (P)𝑃(P)( italic_P ) and the nonlinearity (Q)𝑄(Q)( italic_Q ) in the NLS equation will be complex, given by,

P=12ω2[(cs2+VA2(x))(ωvgk)12igvg],𝑃12superscript𝜔2delimited-[]superscriptsubscript𝑐𝑠2superscriptsubscript𝑉𝐴2𝑥𝜔subscript𝑣𝑔𝑘12𝑖𝑔subscript𝑣𝑔P=\frac{1}{2\omega^{2}}\left[\left({c}_{s}^{2}+{V_{A}}^{2}(x)\right)(\omega-v_% {g}k)-\frac{1}{2}igv_{g}\right],italic_P = divide start_ARG 1 end_ARG start_ARG 2 italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ ( italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) ) ( italic_ω - italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT italic_k ) - divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_i italic_g italic_v start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ] , (44)
Q=12ω(ω2ω1ω26Ω02+igk),𝑄12𝜔superscript𝜔2subscript𝜔1subscript𝜔26superscriptsubscriptΩ02𝑖𝑔𝑘Q=\frac{1}{2\omega}\left(\omega^{2}-\frac{\omega_{1}\omega_{2}}{6\Omega_{0}^{2% }+igk}\right),italic_Q = divide start_ARG 1 end_ARG start_ARG 2 italic_ω end_ARG ( italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - divide start_ARG italic_ω start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_ω start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 6 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_i italic_g italic_k end_ARG ) , (45)

where ω1=(3VA2(x)+2cs2)k2+12Ω02+2igksubscript𝜔13superscriptsubscript𝑉𝐴2𝑥2superscriptsubscript𝑐𝑠2superscript𝑘212superscriptsubscriptΩ022𝑖𝑔𝑘\omega_{1}=\left(3V_{A}^{2}(x)+2c_{s}^{2}\right)k^{2}+12\Omega_{0}^{2}+2igkitalic_ω start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = ( 3 italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) + 2 italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 12 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_i italic_g italic_k and ω2=(3VA2(x)+2cs2)k2+6Ω02+2igksubscript𝜔23superscriptsubscript𝑉𝐴2𝑥2superscriptsubscript𝑐𝑠2superscript𝑘26superscriptsubscriptΩ022𝑖𝑔𝑘\omega_{2}=\left(3V_{A}^{2}(x)+2c_{s}^{2}\right)k^{2}+6\Omega_{0}^{2}+2igkitalic_ω start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = ( 3 italic_V start_POSTSUBSCRIPT italic_A end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_x ) + 2 italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 6 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_i italic_g italic_k. Consequently, instead of the purely growing instability as in Case III of Sec. IV.3, we have the frequency up-shift (Ωr)subscriptΩ𝑟(\Omega_{r})( roman_Ω start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ) and the modulational instability growth rate (Γ)Γ(\Gamma)( roman_Γ ), given by,

Ωr=12[N+N2+4M2]1/2,subscriptΩ𝑟12superscriptdelimited-[]𝑁superscript𝑁24superscript𝑀212\Omega_{r}=\frac{1}{\sqrt{2}}\left[-N+\sqrt{N^{2}+4M^{2}}\right]^{1/2},roman_Ω start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG [ - italic_N + square-root start_ARG italic_N start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT , (46)
Γ=12[N+N2+4M2]1/2,Γ12superscriptdelimited-[]𝑁superscript𝑁24superscript𝑀212\Gamma=\frac{1}{\sqrt{2}}\left[N+\sqrt{N^{2}+4M^{2}}\right]^{1/2},roman_Γ = divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG [ italic_N + square-root start_ARG italic_N start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_M start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT , (47)

where the expressions for M𝑀Mitalic_M and N𝑁Nitalic_N are

M=K2{P1P2K2η0(P2Q1+Q2P1)},𝑀superscript𝐾2subscript𝑃1subscript𝑃2superscript𝐾2subscript𝜂0subscript𝑃2subscript𝑄1subscript𝑄2subscript𝑃1M=K^{2}\left\{P_{1}P_{2}K^{2}-\eta_{0}(P_{2}Q_{1}+Q_{2}P_{1})\right\},italic_M = italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT { italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_Q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) } , (48)
N=K2{2η0(P1Q1+P2Q2)(P12P22)K2}.𝑁superscript𝐾22subscript𝜂0subscript𝑃1subscript𝑄1subscript𝑃2subscript𝑄2superscriptsubscript𝑃12superscriptsubscript𝑃22superscript𝐾2N=K^{2}\left\{2\eta_{0}(P_{1}Q_{1}+P_{2}Q_{2})-(P_{1}^{2}-P_{2}^{2})K^{2}% \right\}.italic_N = italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT { 2 italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_Q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) - ( italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_K start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT } . (49)

For brevity, we have presented the characteristics of ΩrsubscriptΩ𝑟\Omega_{r}roman_Ω start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT and ΓΓ\Gammaroman_Γ as shown in Fig. 7. It is found that both assume higher values at a reduced angular frequency Ω0subscriptΩ0\Omega_{0}roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT (the dotted lines) and by the effects of the density inhomogeneity (the dashed lines). Also, the growth rate ΓΓ\Gammaroman_Γ is similar to Case III of Sec. IV.3 in self-gravitating plasmas, i.e., it grows with the wave number of modulation K𝐾Kitalic_K 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).

Refer to caption
Figure 7: The frequency shift [subplot (a), Eq. (46)] and the instability growth rate [subplot (b), Eq. (47)] are plotted against the modulation wave number K𝐾Kitalic_K of perturbation. The solid, dashed, and dotted lines correspond to the parameter values as for Fig. 2. The other fixed parameter (normalized) values are g=0.5,cs=0.5formulae-sequence𝑔0.5subscript𝑐𝑠0.5g=0.5,~{}c_{s}=0.5italic_g = 0.5 , italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0.5, and k=0.4𝑘0.4k=0.4italic_k = 0.4.

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 (k)𝑘(k)( italic_k ) below the Jeans critical wave number (kJ)subscript𝑘𝐽(k_{J})( italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT ). The latter is, however, defined when the self-gravity force dominates over the Coriolis force or the Jean frequency ωJsubscript𝜔𝐽\omega_{J}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT becomes larger than the increased angular frequency 2Ω02subscriptΩ02\Omega_{0}2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. On the other hand, when the Coriolis force dominates over the self-gravity force or is dominated by the self-gravity force with k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT, 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 ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and k<kJ𝑘subscript𝑘𝐽k<k_{J}italic_k < italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT, is always unstable under the plane wave modulation having no cut-offs in the growth rates. However, the stable Jeans mode with ωJ>2Ω0subscript𝜔𝐽2subscriptΩ0\omega_{J}>2\Omega_{0}italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT > 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT but k>kJ𝑘subscript𝑘𝐽k>k_{J}italic_k > italic_k start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT 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:

iωlBl(n)λξBl(n1)+τBl(n2)+ilk[B0(x)vxl(n)+n=1l=+vxll(nn)Bl(n)]+B0(x)xvxl(n)+B0(x)ξvxl(n1)+n=1l=+ξ(vxll(nn)Bl(n1))=0𝑖𝜔𝑙superscriptsubscript𝐵𝑙𝑛𝜆𝜉superscriptsubscript𝐵𝑙𝑛1𝜏superscriptsubscript𝐵𝑙𝑛2𝑖𝑙𝑘delimited-[]subscript𝐵0𝑥superscriptsubscriptsubscript𝑣𝑥𝑙𝑛superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙superscriptsubscriptsubscript𝑣𝑥𝑙superscript𝑙𝑛superscript𝑛superscriptsubscript𝐵superscript𝑙superscript𝑛subscript𝐵0𝑥𝑥superscriptsubscript𝑣𝑥𝑙𝑛subscript𝐵0𝑥𝜉superscriptsubscript𝑣𝑥𝑙𝑛1superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙𝜉superscriptsubscriptsubscript𝑣𝑥𝑙superscript𝑙𝑛superscript𝑛superscriptsubscript𝐵superscript𝑙superscript𝑛10\begin{split}-i\omega lB_{l}^{(n)}-\lambda\frac{\partial}{\partial\xi}B_{l}^{(% n-1)}+\frac{\partial}{\partial\tau}B_{l}^{(n-2)}+ilk\left[B_{0}(x){v_{x}}_{l}^% {(n)}+\sum_{n^{\prime}=1}^{\infty}\sum_{l^{\prime}=-\infty}^{+\infty}{v_{x}}_{% l-l^{\prime}}^{(n-n^{\prime})}B_{l^{\prime}}^{(n^{\prime})}\right]\\ +\frac{\partial B_{0}(x)}{\partial x}v_{xl}^{(n)}+B_{0}(x)\frac{\partial}{% \partial\xi}v_{xl}^{(n-1)}+\sum_{n^{\prime}=1}^{\infty}\sum_{l^{\prime}=-% \infty}^{+\infty}\frac{\partial}{\partial\xi}({v_{x}}_{l-l^{\prime}}^{(n-n^{% \prime})}B_{l^{\prime}}^{(n^{\prime}-1)})=0\\ \end{split}start_ROW start_CELL - italic_i italic_ω italic_l italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT - italic_λ divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT + divide start_ARG ∂ end_ARG start_ARG ∂ italic_τ end_ARG italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 2 ) end_POSTSUPERSCRIPT + italic_i italic_l italic_k [ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT ] end_CELL end_ROW start_ROW start_CELL + divide start_ARG ∂ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG ∂ italic_x end_ARG italic_v start_POSTSUBSCRIPT italic_x italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT + italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_x italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG ( italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 1 ) end_POSTSUPERSCRIPT ) = 0 end_CELL end_ROW (50)
iωlB0(x)vxl(n)λB0(x)ξvxl(n1)+B0(x)τvxl(n2)+n=1l=+[iωlBll(nn)vxl(n)λBll(nn)ξvxl(n1)+Bll(nn)τvxl(n2)+ilkB0(x)vxll(nn)vxl(n)+B0(x)vxll(nn)ξvxl(n1)]+n,n′′=1l,l′′=+[il′′kBll(nn)vxll′′(nn′′)vxl′′(n′′)+Bll(nn)vxll′′(nn′′)ξvxl′′(n′′1)]+cs2(B0(x)xilkBl(n)+ξBl(n1))+[ilkB0(x)Bl(n)+B0(x)ξBl(n1)+n=1l=+[ilkBll(nn)Bl(n)+Bll(nn)ξBl(n1)+B0(x)B0(x)x+B0(x)xBln]]+2Ω0cosλ[B0(x)vzl(n)+n=1l=+Bll(nn)vzl(n)]2Ω0sinλ[B0(x)vyl(n)+n=1l=+Bll(nn)vyl(n)][B0(x)ψ0(x)x+ψ0(x)xBln+ilkB0(x)ψl(n)+n=1l=+Bll(nn)ξψl(n1)+n=1l=+ilkBll(nn)ψl(n)+B0(x)ξψl(n1)]=0𝑖𝜔𝑙subscript𝐵0𝑥superscriptsubscript𝑣𝑥𝑙𝑛𝜆subscript𝐵0𝑥𝜉superscriptsubscript𝑣𝑥𝑙𝑛1subscript𝐵0𝑥𝜏superscriptsubscript𝑣𝑥𝑙𝑛2superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙delimited-[]𝑖𝜔superscript𝑙superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛superscriptsubscriptsubscript𝑣𝑥superscript𝑙superscript𝑛𝜆superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛𝜉superscriptsubscriptsubscript𝑣𝑥superscript𝑙superscript𝑛1superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛𝜏superscriptsubscriptsubscript𝑣𝑥superscript𝑙superscript𝑛2𝑖superscript𝑙𝑘subscript𝐵0𝑥superscriptsubscriptsubscript𝑣𝑥𝑙superscript𝑙𝑛superscript𝑛superscriptsubscriptsubscript𝑣𝑥superscript𝑙superscript𝑛subscript𝐵0𝑥superscriptsubscriptsubscript𝑣𝑥𝑙superscript𝑙𝑛superscript𝑛𝜉superscriptsubscriptsubscript𝑣𝑥superscript𝑙superscript𝑛1superscriptsubscriptsuperscript𝑛superscript𝑛′′1superscriptsubscriptsuperscript𝑙superscript𝑙′′delimited-[]𝑖superscript𝑙′′𝑘superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛superscriptsubscriptsubscript𝑣𝑥superscript𝑙superscript𝑙′′superscript𝑛superscript𝑛′′superscriptsubscriptsubscript𝑣𝑥superscript𝑙′′superscript𝑛′′superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛superscriptsubscriptsubscript𝑣𝑥superscript𝑙superscript𝑙′′superscript𝑛superscript𝑛′′𝜉superscriptsubscriptsubscript𝑣𝑥superscript𝑙′′superscript𝑛′′1superscriptsubscript𝑐𝑠2subscript𝐵0𝑥𝑥𝑖𝑙𝑘superscriptsubscript𝐵𝑙𝑛𝜉superscriptsubscript𝐵𝑙𝑛1delimited-[]𝑖𝑙𝑘subscript𝐵0𝑥superscriptsubscript𝐵𝑙𝑛subscript𝐵0𝑥𝜉superscriptsubscript𝐵𝑙𝑛1superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙delimited-[]𝑖superscript𝑙𝑘superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛superscriptsubscript𝐵superscript𝑙superscript𝑛superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛𝜉superscriptsubscript𝐵superscript𝑙superscript𝑛1subscript𝐵0𝑥subscript𝐵0𝑥𝑥subscript𝐵0𝑥𝑥superscriptsubscript𝐵𝑙𝑛2subscriptΩ0𝜆delimited-[]subscript𝐵0𝑥superscriptsubscriptsubscript𝑣𝑧𝑙𝑛superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛superscriptsubscriptsubscript𝑣𝑧superscript𝑙superscript𝑛2subscriptΩ0𝜆delimited-[]subscript𝐵0𝑥superscriptsubscriptsubscript𝑣𝑦𝑙𝑛superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛superscriptsubscriptsubscript𝑣𝑦superscript𝑙superscript𝑛delimited-[]subscript𝐵0𝑥subscript𝜓0𝑥𝑥subscript𝜓0𝑥𝑥superscriptsubscript𝐵𝑙𝑛𝑖𝑙𝑘subscript𝐵0𝑥superscriptsubscript𝜓𝑙𝑛superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛𝜉superscriptsubscript𝜓superscript𝑙superscript𝑛1superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙𝑖superscript𝑙𝑘superscriptsubscript𝐵𝑙superscript𝑙𝑛superscript𝑛superscriptsubscript𝜓superscript𝑙superscript𝑛subscript𝐵0𝑥𝜉superscriptsubscript𝜓𝑙𝑛10\begin{split}-i\omega lB_{0}(x)v_{xl}^{(n)}-\lambda B_{0}(x)\frac{\partial}{% \partial\xi}v_{xl}^{(n-1)}+B_{0}(x)\frac{\partial}{\partial\tau}v_{xl}^{(n-2)}% +\sum_{n^{\prime}=1}^{\infty}\sum_{l^{\prime}=-\infty}^{+\infty}\left[-i\omega l% ^{\prime}{B}_{l-l^{\prime}}^{(n-n^{\prime})}{v_{x}}_{l^{\prime}}^{(n^{\prime})% }-\lambda{B}_{l-l^{\prime}}^{(n-n^{\prime})}\frac{\partial}{\partial\xi}{v_{x}% }_{l^{\prime}}^{(n^{\prime}-1)}\right.\\ \left.+{B}_{l-l^{\prime}}^{(n-n^{\prime})}\frac{\partial}{\partial\tau}{v_{x}}% _{l^{\prime}}^{(n^{\prime}-2)}+il^{\prime}kB_{0}(x){v_{x}}_{l-l^{\prime}}^{(n-% n^{\prime})}{v_{x}}_{l^{\prime}}^{(n^{\prime})}+B_{0}(x){v_{x}}_{l-l^{\prime}}% ^{(n-n^{\prime})}\frac{\partial}{\partial\xi}{v_{x}}_{l^{\prime}}^{(n^{\prime}% -1)}\right]+\sum_{n^{\prime},n^{\prime\prime}=1}^{\infty}\sum_{l^{\prime},l^{% \prime\prime}=-\infty}^{+\infty}\left[il^{\prime\prime}k{B}_{l-l^{\prime}}^{(n% -n^{\prime})}{v_{x}}_{l^{\prime}-l^{\prime\prime}}^{(n^{\prime}-n^{\prime% \prime})}{v_{x}}_{l^{\prime\prime}}^{(n^{\prime\prime})}\right.\\ \left.+{B}_{l-l^{\prime}}^{(n-n^{\prime})}{v_{x}}_{l^{\prime}-l^{\prime\prime}% }^{(n^{\prime}-n^{\prime\prime})}\frac{\partial}{\partial\xi}{v_{x}}_{l^{% \prime\prime}}^{(n^{\prime\prime}-1)}\right]+c_{s}^{2}\left(\frac{\partial B_{% 0}(x)}{\partial x}ilk{B}_{l}^{(n)}+\frac{\partial}{\partial\xi}{B}_{l}^{(n-1)}% \right)\\ +\left[ilkB_{0}(x){B}_{l}^{(n)}+B_{0}(x)\frac{\partial}{\partial\xi}{B}_{l}^{(% n-1)}+\sum_{n^{\prime}=1}^{\infty}\sum_{l^{\prime}=-\infty}^{+\infty}\left[il^% {\prime}k{B}_{l-l^{\prime}}^{(n-n^{\prime})}{B}_{l^{\prime}}^{(n^{\prime})}+{B% }_{l-l^{\prime}}^{(n-n^{\prime})}\frac{\partial}{\partial\xi}{B}_{l^{\prime}}^% {(n^{\prime}-1)}+B_{0}(x)\frac{\partial B_{0}(x)}{\partial x}+\frac{\partial B% _{0}(x)}{\partial x}B_{l}^{n}\right]\right]\\ +2\Omega_{0}\cos\lambda\left[B_{0}(x){v_{z}}_{l}^{(n)}+\sum_{n^{\prime}=1}^{% \infty}\sum_{l^{\prime}=-\infty}^{+\infty}{B}_{l-l^{\prime}}^{(n-n^{\prime})}{% v_{z}}_{l^{\prime}}^{(n^{\prime})}\right]-2\Omega_{0}\sin\lambda\left[B_{0}(x)% {v_{y}}_{l}^{(n)}+\sum_{n^{\prime}=1}^{\infty}\sum_{l^{\prime}=-\infty}^{+% \infty}{B}_{l-l^{\prime}}^{(n-n^{\prime})}{v_{y}}_{l^{\prime}}^{(n^{\prime})}% \right]\\ -\left[B_{0}(x)\frac{\partial\psi_{0}(x)}{\partial x}+\frac{\partial\psi_{0}(x% )}{\partial x}B_{l}^{n}+ilkB_{0}(x){\psi}_{l}^{(n)}+\sum_{n^{\prime}=1}^{% \infty}\sum_{l^{\prime}=-\infty}^{+\infty}{B}_{l-l^{\prime}}^{(n-n^{\prime})}% \frac{\partial}{\partial\xi}{\psi}_{l^{\prime}}^{(n^{\prime}-1)}+\sum_{n^{% \prime}=1}^{\infty}\sum_{l^{\prime}=-\infty}^{+\infty}il^{\prime}k{B}_{l-l^{% \prime}}^{(n-n^{\prime})}{\psi}_{l^{\prime}}^{(n^{\prime})}\right.\\ \left.+B_{0}(x)\frac{\partial}{\partial\xi}{\psi}_{l}^{(n-1)}\right]=0\end{split}start_ROW start_CELL - italic_i italic_ω italic_l italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) italic_v start_POSTSUBSCRIPT italic_x italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT - italic_λ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_x italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT + italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG ∂ end_ARG start_ARG ∂ italic_τ end_ARG italic_v start_POSTSUBSCRIPT italic_x italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 2 ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT [ - italic_i italic_ω italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT - italic_λ italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 1 ) end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL + italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_τ end_ARG italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 2 ) end_POSTSUPERSCRIPT + italic_i italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_k italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT + italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 1 ) end_POSTSUPERSCRIPT ] + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_n start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_l start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT [ italic_i italic_l start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT italic_k italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - italic_l start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - italic_n start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL + italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - italic_l start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - italic_n start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT - 1 ) end_POSTSUPERSCRIPT ] + italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG ∂ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG ∂ italic_x end_ARG italic_i italic_l italic_k italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT + divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT ) end_CELL end_ROW start_ROW start_CELL + [ italic_i italic_l italic_k italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT + italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT [ italic_i italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_k italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT + italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_B start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 1 ) end_POSTSUPERSCRIPT + italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG ∂ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG ∂ italic_x end_ARG + divide start_ARG ∂ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG ∂ italic_x end_ARG italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ] ] end_CELL end_ROW start_ROW start_CELL + 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_λ [ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT ] - 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_λ [ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT ] end_CELL end_ROW start_ROW start_CELL - [ italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG ∂ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG ∂ italic_x end_ARG + divide start_ARG ∂ italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG ∂ italic_x end_ARG italic_B start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT + italic_i italic_l italic_k italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) italic_ψ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_ψ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 1 ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT italic_i italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_k italic_B start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_ψ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL + italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_ψ start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT ] = 0 end_CELL end_ROW (51)
iωlvyl(n)λξvyl(n1)+τvyl(n2)+n=1l=+[ilkvx(ll(nn)vyl(n)+vxll(nn)ξvyl(n1)]+2Ω0sinλvxl(n)=0\begin{split}-i\omega lv_{yl}^{(n)}-\lambda\frac{\partial}{\partial\xi}v_{yl}^% {(n-1)}+\frac{\partial}{\partial\tau}v_{yl}^{(n-2)}+\sum_{n^{\prime}=1}^{% \infty}\sum_{l^{\prime}=-\infty}^{+\infty}\left[il^{\prime}k{v_{x}}_{(l-l^{% \prime}}^{(n-n^{\prime})}{v_{y}}_{l^{\prime}}^{(n^{\prime})}+{v_{x}}_{l-l^{% \prime}}^{(n-n^{\prime})}\frac{\partial}{\partial\xi}{v_{y}}_{l^{\prime}}^{(n^% {\prime}-1)}\right]+2\Omega_{0}\sin\lambda{v_{x}}_{l}^{(n)}=0\\ \end{split}start_ROW start_CELL - italic_i italic_ω italic_l italic_v start_POSTSUBSCRIPT italic_y italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT - italic_λ divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_y italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT + divide start_ARG ∂ end_ARG start_ARG ∂ italic_τ end_ARG italic_v start_POSTSUBSCRIPT italic_y italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 2 ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT [ italic_i italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_k italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT ( italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT + italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 1 ) end_POSTSUPERSCRIPT ] + 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_λ italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = 0 end_CELL end_ROW (52)
iωlvzl(n)λξvzl(n1)+τvzl(n2)+n=1l=+[ilkvxll(nn)vzl(n)+vxll(nn)ξvzl(n1)]2Ω0cosλvxl(n)=0𝑖𝜔𝑙superscriptsubscript𝑣𝑧𝑙𝑛𝜆𝜉superscriptsubscript𝑣𝑧𝑙𝑛1𝜏superscriptsubscript𝑣𝑧𝑙𝑛2superscriptsubscriptsuperscript𝑛1superscriptsubscriptsuperscript𝑙delimited-[]𝑖superscript𝑙𝑘superscriptsubscriptsubscript𝑣𝑥𝑙superscript𝑙𝑛superscript𝑛superscriptsubscriptsubscript𝑣𝑧superscript𝑙superscript𝑛superscriptsubscriptsubscript𝑣𝑥𝑙superscript𝑙𝑛superscript𝑛𝜉superscriptsubscriptsubscript𝑣𝑧superscript𝑙superscript𝑛12subscriptΩ0𝜆superscriptsubscriptsubscript𝑣𝑥𝑙𝑛0\begin{split}-i\omega lv_{zl}^{(n)}-\lambda\frac{\partial}{\partial\xi}v_{zl}^% {(n-1)}+\frac{\partial}{\partial\tau}v_{zl}^{(n-2)}+\sum_{n^{\prime}=1}^{% \infty}\sum_{l^{\prime}=-\infty}^{+\infty}\left[il^{\prime}k{v_{x}}_{l-l^{% \prime}}^{(n-n^{\prime})}{v_{z}}_{l^{\prime}}^{(n^{\prime})}+{v_{x}}_{l-l^{% \prime}}^{(n-n^{\prime})}\frac{\partial}{\partial\xi}{v_{z}}_{l^{\prime}}^{(n^% {\prime}-1)}\right]-2\Omega_{0}\cos\lambda{v_{x}}_{l}^{(n)}=0\\ \end{split}start_ROW start_CELL - italic_i italic_ω italic_l italic_v start_POSTSUBSCRIPT italic_z italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT - italic_λ divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_z italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT + divide start_ARG ∂ end_ARG start_ARG ∂ italic_τ end_ARG italic_v start_POSTSUBSCRIPT italic_z italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - 2 ) end_POSTSUPERSCRIPT + ∑ start_POSTSUBSCRIPT italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = - ∞ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + ∞ end_POSTSUPERSCRIPT [ italic_i italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_k italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT + italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l - italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n - italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) end_POSTSUPERSCRIPT divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_v start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 1 ) end_POSTSUPERSCRIPT ] - 2 roman_Ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_λ italic_v start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT = 0 end_CELL end_ROW (53)
2ψ0(x)x2l2k2ψl(n)+2ilkξψl(n1)+2ξ2ψl(n2)+ωJ2B0(x)+ωJ2Bl(n)=0superscript2subscript𝜓0𝑥superscript𝑥2superscript𝑙2superscript𝑘2subscriptsuperscript𝜓𝑛𝑙2𝑖𝑙𝑘𝜉subscriptsuperscript𝜓𝑛1𝑙superscript2superscript𝜉2subscriptsuperscript𝜓𝑛2𝑙superscriptsubscript𝜔𝐽2subscript𝐵0𝑥superscriptsubscript𝜔𝐽2subscriptsuperscript𝐵𝑛𝑙0\begin{split}\frac{\partial^{2}\psi_{0}(x)}{\partial x^{2}}-l^{2}k^{2}\psi^{(n% )}_{l}+2ilk\frac{\partial}{\partial\xi}\psi^{(n-1)}_{l}+\frac{\partial^{2}}{% \partial\xi^{2}}\psi^{(n-2)}_{l}+\omega_{J}^{2}B_{0}(x)+\omega_{J}^{2}B^{(n)}_% {l}=0\\ \end{split}start_ROW start_CELL divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ψ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG ∂ italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - italic_l start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ψ start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT + 2 italic_i italic_l italic_k divide start_ARG ∂ end_ARG start_ARG ∂ italic_ξ end_ARG italic_ψ start_POSTSUPERSCRIPT ( italic_n - 1 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT + divide start_ARG ∂ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ∂ italic_ξ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_ψ start_POSTSUPERSCRIPT ( italic_n - 2 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT + italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_x ) + italic_ω start_POSTSUBSCRIPT italic_J end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_B start_POSTSUPERSCRIPT ( italic_n ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT = 0 end_CELL end_ROW (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