Issue |
A&A
Volume 678, October 2023
|
|
---|---|---|
Article Number | A94 | |
Number of page(s) | 9 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/202347149 | |
Published online | 11 October 2023 |
Excitation of extraordinary modes inside the source of Saturn’s kilometric radiation
1
Institute of Frontier and Interdisciplinary Science, Shandong University,
Qingdao, Shandong
266237, PR China
e-mail: yaochen@sdu.edu.cn
2
Institute of Space Sciences, Shandong University,
Shandong
264209, PR China
3
Department of Earth and Space Sciences, Southern University of Science and Technology,
Shenzhen, Guangdong,
PR China
4
Institute of Solar-Terrestrial Physics,
Irkutsk
664033, Russia
Received:
11
June
2023
Accepted:
23
August
2023
The electron cyclotron maser instability (ECMI) of extraordinary mode waves was investigated with the parameters observed in Saturn’s kilometric radiation (SKR) sources. Previous studies employed simplified dispersion relations, and did not consider the excitation of the relativistic (R) mode. This mode is introduced by considering the relativistic effect in plasmas consisting of both cold and hot electrons. Using particle-in-cell simulations, we investigated the excitation of R and X modes based on the measured data. Using the reported value of the density ratio of energetic to total electrons ne/n0 = 24%, the most unstable mode is the R mode. The escaping X-mode emissions are amplified only if the energetic electrons are dominant with ne/n0 ≥ 90%. For these cases, only the X mode is excited and the R mode disappears due to its strong coupling. The results are well in line with the linear kinetic theory of ECMI. The properties of both the R and X modes are consistent with the observed SKR emissions. This raises questions about the nature of the measured electric field fluctuations within “presumed” SKR sources. The study provides new insights into the ECMI process relevant to SKR emission mechanisms.
Key words: planets and satellites: gaseous planets / radio continuum: planetary systems / masers / waves / plasmas / methods: numerical
© The Authors 2023
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1 Introduction
Saturn’s kilometric radiations (SKRs) represent intense radio emissions originating from the kronian auroral region, considered to be the counterpart of the auroral kilometric radiations (AKRs) on Earth (see, e.g., Kaiser et al. 1980, 1981; Gurnett et al. 2004; Kurth et al. 2005; Lamy et al. 2009). The radiations are observed in the kilohertz to megahertz frequency range, and the radiated powers of 108–1010 W. Furthermore, SKRs are closely related to Saturn’s magnetosphere dynamics (Jackman 2011; Kimura et al. 2013), and have been extensively studied.
Saturn’s kilometric radiations are often strongly polarized with narrow beaming angles and high intensity (Lamy et al. 2008, 2018), indicating coherent emissions generated by electron cyclotron maser instability (ECMI; Wu & Lee 1979; Wu 1985; Treumann 2006; Bingham et al. 2013). The instability excites escaping modes at gyro-frequencies (~Ωce) via a gyro-resonant wave-particle interaction, driven by electrons with anisotropic velocity distribution functions (VDFs; see Ning et al. 2021a,b, for recent studies). The ECMI is now widely applied to planetary and stellar radio emissions.
The Cassini spacecraft measures the radio waves with the Radio and Plasma Wave Science (RPWS, Gurnett et al. 2004) instrument and the electrons with the Cassini Plasma Spectrometer Electron Spectrometer (CAPS/ELS, Linder et al. 1998; Young et al. 2004) in Saturn’s magnetosphere. Since 2010, researchers have reported several crossings of SKR sources by Cassini, identified with the cutoff frequency of the spectrum close to or slightly lower than the local Ωce (e.g., Kurth et al. 2011; Mutel et al. 2010; Lamy et al. 2010, 2018; Menietti et al. 2011).
As has previously been reported, SKR sources are characterized by a low ratio of plasma oscillation frequency to electron gyro-frequency (ωpe/Ωce < 0.1), in accordance with the ECMI conditions (ωpe/Ωce < 1, Wu & Lee 1979). The emissions mainly propagate quasi-perpendicularly to the background magnetic field with linear polarization near the sources, consistent with the X-mode properties. The energy conversion rate from electrons to radiations is ~1% (Lamy et al. 2011). The abundances of energetic and cold electrons are comparable with each other (Schippers et al. 2011).
The measurement of electron VDFs is crucial to the study of the ECMI process. In the first encounter of the SKR source by Cassini in 2008, Mutel et al. (2010) found that the cold and hot electrons can be fitted with a Kappa and a Dory-Guest-Harri (DGH, Dory et al. 1965) distribution, respectively. They further calculated the corresponding ECMI growth rates of X-mode emissions with the linear kinetic theory, and suggested that the ECMI process can efficiently amplify this escaping X mode below Ωce, along directions perpendicular to the magnetic field. Menietti et al. (2011) investigated the ECMI process of another SKR source using a similar method and obtained the same results. They also suggested that O-mode and Z-mode waves can be generated.
However, these studies on the ECMI of SKRs followed Mutel et al. (2007) and employed the dispersion relations of pure relativistic energetic electrons. They estimated the correction to the dispersion relation with a scale factor when calculating the growth rates of the X and Z modes, rather than solving the dispersion matrix to calculate the real frequency and the growth rate in a self-consistent manner. As reported by Mutel et al. (2010) and Menietti et al. (2011), the density ratio of the energetic to the total electrons was fitted as ne/n0 = 20–24%. With this range of ne/n0 (and other comparable values), one obtains very different dispersion relations. For waves perpendicular to the magnetic field, the gyro-resonance condition is simplified as . When adopting the cold-plasma dispersion relation deduced from the magneto-ionic theory, the ECMI process cannot amplify the fundamental X mode with the cutoff frequency (ωX) being larger than Ωce. For pure energetic electrons, the X-mode cutoff is less than Ωce and thus can grow via ECMI (e.g., Winglee 1983, 1985; Pritchett 1984a; Robinson 1986).
For plasmas consisting of both the cold and energetic electrons, another important mode may emerge. This mode, first deduced by Pritchett (1984b), appears in the range of Ωce/γ < ω < Ωce. It is relevant to the relativistic effect, and is therefore called the relativistic (R) mode (Strangeway 1985, 1986). Our present study shows that the R mode becomes the most unstable mode, and is probably relevant to SKR. However, it is trapped and cannot escape. This mode has not been considered in studies on SKR sources, to the best of our knowledge.
According to the earlier reports on the crossings of SKR sources by Cassini, both cold and hot components of electrons exist with comparable densities (Schippers et al. 2011). This raises questions concerning the applicability of the oversimplified dispersion relations adopted by Mutel et al. (2010) and Menietti et al. (2011), which did not consider the existence of the R mode. For this study, we conducted fully kinetic electromagnetic particle-in-cell (PIC) simulations with the electron VDFs and plasma parameters reported by Menietti et al. (2011) to investigate the extraordinary modes (including R, X, and Z) and their coupling within SKR sources. We conducted a parametric study on ne/n0, and applied the kinetic linear theory to verify and explain the PIC simulations.
![]() |
Fig. 1 Evolution of electron VDFs and energy profiles of the simulation with ne/n0 = 24%. Panels a–c: snapshots of the VDFs at the start (a, t = 0 |
2 Numerical method and parameter setting
As reported by Menietti et al. (2011), the frequency ratio was ωpe/Ωce = 0.046, and the VDFs of both the cold and hot electrons were written as
(1)
respectively, where j = 11, u represents the momentum per mass, and υ|| (= 0.11c) and υ⊥ (= 0.038c) represent the thermal velocity of the DGH distribution along parallel and perpendicular directions. The density ratio of the energetic to the total electrons was fitted to be ne/n0 = 24%. Considering the uncertainties of the source measurements, we conducted a parametric study of ne/n0. The models and parameter settings are presented below.
2.1 PIC simulations
We employed the Vector-PIC (vPIC) code developed by the Los Alamos National Laboratory (LANL). The code employs a second-order, explicit, leapfrog algorithm to update the positions and the velocities of charged particles, along with a full Maxwell description of the electric and magnetic fields via a second-order finite-difference time-domain solver (Bowers et al. 2008, 2009).
We performed 2D3V (two spatial dimensions with three velocity components) simulations in the xOz plane with periodical boundary conditions and a uniform background magnetic field (B0 = B0êz). The simulation domain was set to Lx = Lz = 2048Δ, the grid length to Δ = 0.019 λDE, and the time step to .
Two components of electrons were included. Their distribution functions are described by Eqs. (1) and (2), and the center of the VDF of energetic electrons is at υ⊥ ~ 0.18 c (Fig. 1a). The thermal velocity of the low-energy electrons was set to υ0 = 0.01c (~50 eV). The charge neutrality was maintained by including thermal protons with a realistic proton-to-electron mass ratio of 1836 and the same temperature (with the thermal velocity . We included 1000 macroparticles per cell for each component of the electrons and 500 for the protons.
![]() |
Fig. 2 Wave dispersion diagrams of electric and magnetic fields with ne/n0 = 24% and θkB = 90°, during the interval of [100, 200] |
2.2 Linear analysis of ECMI
Following the analysis of Pritchett (1984b), we employed the δ function below to describe the two components of electrons,
(3)
where ur refers to the momentum of the energetic electrons. In accordance with the PIC simulation, we assumed the waves propagate in the xOz plane with B0 = Bzêz. For simplicity, only the perpendicular propagating waves with k|| = 0 were studied.
With the VDF presented above, we derived the elements of the plasma dispersion tensor (Eq. (A.1)) and obtained the frequencies and the growth rates of the ECMI modes by numerically solving their determinant. We analyzed the dispersion relations, the ECMI growth rates, and the polarization of the X and R modes. The results are presented in the appendix.
3 PIC simulations
We present the PIC simulations with the reported SKR source parameters as the reference case with (ne/n0 = 24%) to show the excitation and properties of extraordinary modes (R, X, Z), followed by the parametric study on ne/n0 to explore its effect.
3.1 Results with measured parameters (ne/n0=0.24)
According to the energy profiles presented in Fig. 1d, the wave excitation and further evolution can be divided into two stages, the linear (0–120 ) stage and the saturation (120–200
) stage. In the linear stage, the energies of ΔBZ (= Bz – B0) and Ey rise rapidly, reaching the level of 10−2 Ek0 (Ek0 is the initial kinetic energy of total electrons). The electron kinetic energy (Ek) declines in accordance with the growth of the ECMI modes. The energies of other field components are much weaker, being ~10−5 Ek0. In the saturation stage, the energies of all the fields rise gradually. In total, ~3.4% of Ek0 was converted to wave energies.
Figures 1a–c present the evolution of the electron VDFs. The two components are initially separated. Within 100–200 , the shell-like hot electrons diffuse in the phase space, later developing a butterfly-like morphology. Meanwhile, the cool Maxwellian electrons are significantly energized along the perpendicular direction.
We performed a fast Fourier transform (FFT) analysis on each field component. In Figs. 2a–f, we show the ω – k dispersion diagrams for perpendicular-propagating waves with θkB = 90° (the angle between k and B0), along which the ECMI waves are the most amplified. The strongest wave appears in the Ey and Bz components, with frequencies of ~0.98 Ωce, in accordance with the dispersion relation of the R mode as shown in the Appendix. This mode can be identified as the transverse electromagnetic wave with significant Ey and Bz components and an insignificant Ex component. The property of the simulated R mode is in line with the linear kinetic theory (see Appendix A.2).
To separate the left-hand (LH) and right-hand (RH) polarized components of the R mode, we processed (Ex, Ey) with the following equation:
(4)
The results are presented in Figs. 2g–h. The intensities of the LH and RH components are comparable to each other, indicating a strong linear polarization.
According to the energy profiles in Fig. 1d, the total energy of the R mode is ~0.034 Ek0, and its growth rate can be fitted as Γ = 5.57 × 10−3 Ωce. The fundamental O mode is slightly amplified in Ez and By. Its total energy is weaker than the R mode by approximately two orders of magnitude.
In this case, only the R mode can be amplified, with no significant growth of X. The results can be clearly understood with the linear kinetic theory. As shown in the Appendix, the R mode is the most unstable at k ~ Ωce/c with a maximal growth rate of 0.006 Ωce, and linearly polarized in Ey with Ex/iEy ~ 0 (Figs. A.1a–c). The growth rates of the X mode remain at 0 with its cutoff being larger than Ωce.
![]() |
Fig. 3 Energy proflies and growth rates for different simulation cases. Panels a–e: energy profiles of the six fluctuated field components for the cases with ne/n0 = 0.1, 0.5, 0.7, 0.9, and 1.0. The energies are normalized to the respective initial kinetic energy of total electrons of each case. Panel f presents the growth rates with different density ratios obtained by the PIC simulations (blue) and the linear kinetic theory (orange). |
3.2 Effect of ne/n0
Considering the uncertainties of the electron measurement of Cassini CAPS/ELS, we varied the values of ne/n0 from 10% to 100% to investigate the conditions for X-mode excitations. For cases with a ne/n0 > 10%, the simulations lasted for 200 , and for the ne/n0 = 10% case, the simulation lasted for 300
.
In Figs. 3a–e, we plotted the energy profiles of each case. As introduced in Sect. 3.1, we divided the whole ECMI process into two stages, the linear stage and the saturation stage. For each case we observed a sharp rise of Ey in the linear stage. The energy conversion rates from the electron kinetic energy to the waves is ~10−2. With a higher ne/n0, the wave excitation is more impulsive, with a short linear stage. For cases with ne/n0 ≥ 50%, we observed a gradual decline of the field energies and −ΔEk (= Ek0 – Ek), suggesting further damping of the amplified waves. Figure 4 displays the evolution of the VDFs in different cases. With a larger ne/n0, the VDF of the hot electrons diffuses faster, and the cold electrons are more efficiently accelerated. The VDFs do not change much after 50–100 .
For each case, we performed a FFT analysis of the fluctuated fields to show the amplified wave modes (Fig. 5). The waves in Ex and Ez are not shown here since they are very weak. For cases with ne/n0 ≤ 70%, the R mode is the most intensive, while no significant X mode appears. With increasing ne/n0, the R-mode dispersion curve becomes flatter and closer to that of the X mode. With ne/n0 ≥ 90%, the two dispersion curves merge, indicating their strong coupling; this leads to the significant excitation of the X mode with a lower cutoff frequency of ωX ~ 0.97 Ωce.
We plotted the growth rates (Γ) of the simulated R mode (or the X mode for the cases with ne/n0 = 90 and 100%) by fitting the energy curves (Fig. 3f). The values of Γ rise from 0.0032 to 0.018 Ωce with ne/n0 from 10% to 100%. With the linear kinetic theory (see Fig. A.2), we obtained comparable growth rates for the R (for the cases with 0 < ne/n0 < 100%) and the X mode (for the case with ne/n0 = 100%), which rise from 4.59 × 10−3 to 1.14 × 10−2 Ωce.
According to our PIC simulations, the X mode is excited only if the energetic component is dominant (ne/n0 ≥ 90%). The characteristics of the X mode in the case of ne/n0 = 100%, including the growth rate, polarization, and wave fields (see Figs. 3e, 4e, and 6), are consistent with the linear kinetic theory (Figs. A.1d–f). In this case, the X mode becomes saturated within the first 50 , mainly in Ey and Bz, and much less so in Ex. The energy conversion rate from energetic electrons to the X mode is ~3%, and the growth rate is ~0.018 Ωce, similar to that obtained by the linear kinetic theory (~0.014 Ωce; see Fig. A.1e). According to Figs. 6g–h, the amplified X-mode waves are linearly polarized with comparable LH and RH intensities.
With the linear kinetic theory (see Appendix A.3), for the case with ne/n0 = 90%, the growth rate of the X mode is 0. This is inconsistent with the corresponding PIC simulation presented above. This can be explained by the evolution of the VDFs (Fig. 4d). Almost all of the cool electrons are energized in the simulation within the first 50 , and their VDFs overlap with the hot component. This makes it difficult to distinguish the two components in the phase space, leading to an equivalent case with ne/n0 = 100%. This explains why the simulation with ne/n0 = 90% is in line with that of ne/n0 = 100%. In other cases with less energetic electrons, the two components of the VDF are still separated until the end of the simulation, and the X mode cannot grow. Further discussion is provided in the following section.
![]() |
Fig. 4 Snapshots of VDFs of simulations with ne/n0 = 10% (a), 50% (b), 70% (c), 90% (d), and 100% (e) at moments representative of the initial (left), the transition (middle), and the saturation (right) stages. |
![]() |
Fig. 5 ω – k dispersion diagrams of Ey at θkB = 90° during the saturation stage for different cases. |
4 Conclusion and discussion
The present study investigates the excitation of the R and X modes based on the in situ observation of SKR sources reported by Menietti et al. (2011). Both PIC simulations and linear theory analysis were performed, and the results are in line with each other. With the given set of parameters, the trapped R mode is amplified as the dominant mode, while the escaping X mode cannot grow. The growth rate of the R mode increases with a ne/n0 of 10–70%. The X-mode emissions can only be amplified if the electrons are dominated by hot components (ne/n0 = 90% and 100%) with the cutoff frequency less than Ωce and the growth rate ~0.018 Ωce. According to our PIC simulations, both the X and R modes can be characterized with similar frequencies (~Ωce), polarization (linearly polarized in Ey), growth rates (10−2 Ωce), and energy conversion rates (1–3%). The properties can be well explained using the linear kinetic theory.
As the kinetic theory suggests, the X mode can be amplified only with ne/n0 = 100%. In other cases with 0 < ne/n0 ≤ 90%, the R mode is the most unstable and the X mode cannot grow. In our simulation with ne/n0 = 90%, we observed significant declines of wave energies after the linear stage, with the cold electrons significantly accelerated at 50 . We propose a multi-step process to explain the X-mode excitation in this case (and its inconsistency with the linear theory): the R-mode waves are first amplified via ECMI, then interact strongly with the cool electrons and become damped soon via the gyro-resonant absorption. This results in the bulk acceleration of the minor cool electrons, which then become the “hot” component. This leads to an effective “ne/n0 = 100%” case, and thus the strong coupling of the R and X modes and the strong excitation of the X mode. This means that the non-escaping R mode plays a crucial role in the dynamic evolution of electron VDFs and the excitation of the X mode.
However, with the actual measured parameters, we obtained significant excitation of the trapped R mode rather than the escaping X-mode radiation. This does not explain SKR observed from a distance. We propose two possibilities to explain this problem:
The measured electron VDFs are not the direct driver of ECMI due to the poor resolution. The CAPS/ELS instrument measures the electrons of different pitch angles using eight detectors, which rotate for approximately 3 min to cover 2π sr of the sky. The temporal resolution is much larger than the timescales of the wave growth (~100
). It is hard to measure electrons right before the excitation of ECMI since the VDFs diffuse rapidly at a time scale of 50–100
. Therefore, the ECMI of SKRs may not be induced by the measured electron VDFs.
The in situ observation cannot distinguish the R and X modes in the source region. According to our simulations and the linear theory analysis, both the R and X modes are identified as transverse electromagnetic waves with ω ~ Ωce, linearly polarized in Ey. These properties are consistent with the observed SKR emissions. With the Cassini in situ measurements of the electric field only, it is difficult to distinguish the two wave modes. In other words, it is possible that the enhancement of the fluctuating electric fields observed by Cassini near or within the sources is actually the trapped R mode, rather than the generally believed escaping X mode.
Saturn’s kilometric radiation has long been considered to be the analogy of the AKR from the geo-magnetosphere. Similar problems were present when applying the ECMI mechanism to AKRs. It is believed that AKRs are generated from low-density cavities with no cold electron population. According to the PIC simulations of the finite source cavity (Pritchett 1986; Pritchett & Winglee 1989; Pritchett et al. 2002), the generated X-mode radiations could escape from the cavities via mode conversion at the boundaries. Although the presence of such cavity structures has not been conclusively identified in SKR sources, this could be a possible explanation for SKR.
In the present simulations, we employed periodical boundary conditions to model the excitation of the wave modes, assuming the electrons are uniformly distributed in a large spatial scale. No replenishment of the electron free energy was included in this model. We will conduct simulations with various setups of the source configurations to examine their effect on the generation process of SKR.
The ECMI growth rates of the X and Z modes in SKR sources have been studied by Mutel et al. (2010) and Menietti et al. (2011). As we pointed out in Sect. 1, previous studies employed oversimplified dispersion relations with inappropriate assumptions and did not consider the R mode. In this paper, we performed PIC simulations to study the excitation of the three extraordinary modes (R, X, and Z), focusing mainly on the R mode and its effect on the excitation of the X mode. We also performed linear analysis of the plasma kinetic theory to obtain the dispersion relations and wave growths. The results are in agreement with the PIC simulations.
The R mode is dominant in most cases with both the cold and hot electrons, reaching energies of ~0.03 Ek0. This mode is trapped between plasma layers delimited by the two electron gyro-frequencies (Ωce/γ, Ωce), yet it may be converted to escaping emissions when propagating in inhomogeneous plasmas. In addition, the R mode can be absorbed by electrons, leading to the acceleration and heating of plasmas, which, in turn, can modify the dispersion relations. Further simulations should investigate the heating and acceleration process due to the R mode and its coupling with the escaping X mode.
![]() |
Fig. 6 Wave dispersion diagrams of the six field components (a–f) and the separated LH and RH components (g–h) with ne/n0 = 100% and θkB = 90° during [100, 200] |
Acknowledgements
This study is supported by NSFC grants (11973031, 12203031, and 12103029), the China Postdoctoral Science Foundation (2022TQ0189), and the Natural Science Foundation of Shandong Province (ZR2021QA033). We thank the National Supercomputer Centers in Tianjin and the Beijing Super Cloud Computing Center (BSCC, URL: http://www.blsc.cn/) for providing HPC resources, and LANL for the open-source VPIC code.
Appendix A Linear analysis of plasma kinetic theory
A.1 Dispersion tensor
We followed Pritchett (1984b) to derive the dielectric tensors of plasmas with the VDF given by the δ function (Eq. (3)). For perpendicular propagating waves, the elements of plasma dispersion tensors were simplified as
(A.1)
where r = ne/n0, ξ = kur/Ωce, r′ = 1 – r, Gr = γr – 1/γr, , b = γrω – nΩce, and Jn(ξ) refers to the usual Bessel function of order n.
The following equation was solved to obtain both the frequency and growth rate of each extraordinary mode (X, Z, or R):
(A.2)
A.2 ECMI of X and R with ne/n0 = 24%
To understand the PIC simulations, we first conducted linear theory analysis with the parameters of SKR sources reported by Menietti et al. (2011). We set ur = 0.18 c, γr = 1.016, ne/n0 = 24%, and ωpe/Ωce = 0.046.
With a given k, we obtained four sets of complex roots ω = ωr + iΓ of Eq. (A.2), representing different branches of the wave modes (see Fig. A.1(a–b)). The dispersion curves of the X mode lie in the range of ωr > Ωce with a zero growth rate. The slow branch of the extraordinary mode (Z) is below Ωce/γr. Within the range of Ωce/γr < ωr < Ωce and k ≲ Ωce/c, there exist two conjugate roots of the R mode (R1 and R2). The maximum absolute value of their growth rates is 0.006 Ωce. With increasing k (> Ωce/c), R1 and R2 are decoupled and approach their respective resonant frequency, that is, ~Ωce and Ωce/γr.
The perpendicular-propagating extraordinary waves are mainly carried by Ex and Ey. The polarization was determined by
(A.3)
Figure A.1(c) displays the variations in Ex/iEy with k for each extraordinary mode. With increasing k, Ex/iEy rises from –1 to 0 for the X mode and from 0 to infinity for the Z mode. The Ex/iEy of R coincides with that of the X mode, with k < Ωce/c, and with that of the Z mode, with k > Ωce/c.
![]() |
Fig. A.1 Dispersion relation (a, d) and growth rates (b, e) of the extraordinary wave modes according to the kinetic linear theory, with ne/n0 = 24% (upper panels) and 100% (lower panels). The blue and red lines refer to the X and Z modes; the purple and yellow lines refer to the two branches of the relativistic modes (R1 and R2). Panels (c) and (f) present the corresponding intensity ratio of Ex/iEy for each mode. |
With ne/n0 = 24%, R is the most unstable mode, with the highest growth rate at k ~ Ωce/c; where the three modes are basically linearly polarized with Ex/iEy ~ 0, the X mode cannot be amplified via the ECMI process.
A.3 Parameter studies on ne/n0, ωpe/Ωce, and ur
We varied ne/n0 from 0 to 100%, keeping the other parameters fixed (Fig. A.2). With ne/n0 = 0, we obtained the magnetoionic cold-plasma dispersion relation with two branches (X and Z), with no instability.
With ne/n0 > 0, we obtained conjugate pairs of solutions for the R1, R2, and Z modes, and only those with positive growth rates are presented. The unstable R1 and R2 modes appear within Ωce/γ < ω < Ωce and their growth rates rise with ne/n0. With a higher abundance of energetic electrons, the X mode has a lower cutoff (ωX) and its dispersion curve becomes closer to that of the R mode. The X mode is always stable with ωX > Ωce if ne/n0 < 100%.
With ne/n0 = 100%, the dispersion curve of the R and X modes merge together, with the cutoff ωX < Ωce. In this case, the X mode can be amplified via ECMI with the highest growth rate, ~0.011 Ωce. As Figure A.1(d–f) shows, the X mode is linearly polarized in Ey (Ex/iEy ~ 0) with k ≳ Ωce/c, and has the highest growth rate at k ~ Ωce/c. In summary, the linear theory suggests that the R mode is unstable if both cold and hot components are present, while the X mode can only be amplified with ne/n0 = 100%.
We further expanded the parameter study by varying ωpe/Ωce and ur while fixing ne/n0 (=24%). We find that the R mode is always dominant. With increasing ωpe/Ωce (from 0.01 to 0.25) or ur (from 0.05 to 0.3 c), the dispersion curves of the R mode appear in a broader range of frequencies with increasing growth rates.
![]() |
Fig. A.2 Dispersion relation (upper panels) and growth rates (lower panels) of the extraordinary wave modes (X, R1, R2, and Z) with varying ne/n0. |
References
- Bingham, R., Speirs, D. C., Kellett, B. J., et al. 2013, Space Sci. Rev., 178, 695 [Google Scholar]
- Bowers, K. J., Albright, B. J., Yin, L., Bergen, B., & Kwan, T. J. T. 2008, Phys. Plasmas, 15, 055703 [NASA ADS] [CrossRef] [Google Scholar]
- Bowers, K. J., Albright, B. J., Yin, L., et al. 2009, J. Phys. Conf. Ser., 180, 012055 [Google Scholar]
- Dory, R. A., Guest, G. E., & Harris, E. G. 1965, Phys. Rev. Lett., 14, 131 [CrossRef] [Google Scholar]
- Gurnett, D. A., Kurth, W. S., Kirchner, D. L., et al. 2004, Space Sci. Rev., 114, 395 [NASA ADS] [CrossRef] [Google Scholar]
- Jackman, C. M. 2011, in Planetary, Solar and Heliospheric Radio Emissions (PRE VII), eds. H. O. Rucker, W. S. Kurth, P. Louarn, & G. Fischer, 1 [Google Scholar]
- Kaiser, M. L., Desch, M. D., Warwick, J. W., & Pearce, J. B. 1980, Science, 209, 1238 [CrossRef] [Google Scholar]
- Kaiser, M. L., Desch, M. D., & Lecacheux, A. 1981, Nature, 292, 731 [NASA ADS] [CrossRef] [Google Scholar]
- Kimura, T., Lamy, L., Tao, C., et al. 2013, J. Geophys. Res. Space Phys., 118, 7019 [NASA ADS] [CrossRef] [Google Scholar]
- Kurth, W. S., Gurnett, D. A., Clarke, J. T., et al. 2005, Nature, 433, 722 [CrossRef] [Google Scholar]
- Kurth, W. S., Gurnett, D. A., Menietti, J. D., et al. 2011, in Planetary, Solar and Heliospheric Radio Emissions (PRE VII), eds. H. O. Rucker, W. S. Kurth, P. Louarn, & G. Fischer, 75 [Google Scholar]
- Lamy, L., Zarka, P., Cecconi, B., et al. 2008, J. Geophys. Res. Space Phys., 113, A07201 [NASA ADS] [Google Scholar]
- Lamy, L., Cecconi, B., Prangé, R., et al. 2009, J. Geophys. Res. Space Phys., 114, A10212 [NASA ADS] [Google Scholar]
- Lamy, L., Schippers, P., Zarka, P., et al. 2010, Geophys. Res. Lett., 37, L12104 [Google Scholar]
- Lamy, L., Cecconi, B., Zarka, P., et al. 2011, J. Geophys. Res. Space Phys., 116, A04212 [Google Scholar]
- Lamy, L., Zarka, P., Cecconi, B., et al. 2018, Science, 362, aat2027 [NASA ADS] [CrossRef] [Google Scholar]
- Linder, D. R., Coates, A. J., Woodliffe, R. D., et al. 1998, Geophys. Monogr. Ser., 102, 257 [NASA ADS] [Google Scholar]
- Menietti, J. D., Mutel, R. L., Schippers, P., et al. 2011, J. Geophys. Res. Space Phys., 116, A12222 [NASA ADS] [Google Scholar]
- Mutel, R. L., Peterson, W. M., Jaeger, T. R., & Scudder, J. D. 2007, J. Geophys. Res. Space Phys., 112, A07211 [Google Scholar]
- Mutel, R. L., Menietti, J. D., Gurnett, D. A., et al. 2010, Geophys. Res. Lett., 37, L19105 [Google Scholar]
- Ning, H., Chen, Y., Ni, S., et al. 2021a, A&A, 651, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ning, H., Chen, Y., Ni, S., et al. 2021b, ApJ, 920, L40 [NASA ADS] [CrossRef] [Google Scholar]
- Pritchett, P. L. 1984a, Geophys. Res. Lett., 11, 143 [NASA ADS] [CrossRef] [Google Scholar]
- Pritchett, P. L. 1984b, J. Geophys. Res., 89, 8957 [NASA ADS] [CrossRef] [Google Scholar]
- Pritchett, P. L. 1986, J. Geophys. Res., 91, 13569 [NASA ADS] [CrossRef] [Google Scholar]
- Pritchett, P. L., & Winglee, R. M. 1989, J. Geophys. Res., 94, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Pritchett, P. L., Strangeway, R. J., Ergun, R. E., & Carlson, C. W. 2002, J. Geophys. Res. (Space Phys.), 107, 1437 [NASA ADS] [CrossRef] [Google Scholar]
- Robinson, P. A. 1986, J. Plasma Phys., 35, 187 [Google Scholar]
- Schippers, P., Arridge, C. S., Menietti, J. D., et al. 2011, J. Geophys. Res. Space Phys., 116, A05203 [NASA ADS] [CrossRef] [Google Scholar]
- Strangeway, R. J. 1985, J. Geophys. Res., 90, 9675 [NASA ADS] [CrossRef] [Google Scholar]
- Strangeway, R. J. 1986, J. Geophys. Res., 91, 3152 [NASA ADS] [CrossRef] [Google Scholar]
- Treumann, R. A. 2006, A&ARv, 13, 229 [Google Scholar]
- Winglee, R. M. 1983, Plasma Phys., 25, 217 [NASA ADS] [CrossRef] [Google Scholar]
- Winglee, R. M. 1985, ApJ, 291, 160 [NASA ADS] [CrossRef] [Google Scholar]
- Wu, C. S. 1985, Space Sci. Rev., 41, 215 [NASA ADS] [CrossRef] [Google Scholar]
- Wu, C. S., & Lee, L. C. 1979, ApJ, 230, 621 [Google Scholar]
- Young, D. T., Berthelier, J. J., Blanc, M., et al. 2004, Space Sci. Rev., 114, 1 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
![]() |
Fig. 1 Evolution of electron VDFs and energy profiles of the simulation with ne/n0 = 24%. Panels a–c: snapshots of the VDFs at the start (a, t = 0 |
In the text |
![]() |
Fig. 2 Wave dispersion diagrams of electric and magnetic fields with ne/n0 = 24% and θkB = 90°, during the interval of [100, 200] |
In the text |
![]() |
Fig. 3 Energy proflies and growth rates for different simulation cases. Panels a–e: energy profiles of the six fluctuated field components for the cases with ne/n0 = 0.1, 0.5, 0.7, 0.9, and 1.0. The energies are normalized to the respective initial kinetic energy of total electrons of each case. Panel f presents the growth rates with different density ratios obtained by the PIC simulations (blue) and the linear kinetic theory (orange). |
In the text |
![]() |
Fig. 4 Snapshots of VDFs of simulations with ne/n0 = 10% (a), 50% (b), 70% (c), 90% (d), and 100% (e) at moments representative of the initial (left), the transition (middle), and the saturation (right) stages. |
In the text |
![]() |
Fig. 5 ω – k dispersion diagrams of Ey at θkB = 90° during the saturation stage for different cases. |
In the text |
![]() |
Fig. 6 Wave dispersion diagrams of the six field components (a–f) and the separated LH and RH components (g–h) with ne/n0 = 100% and θkB = 90° during [100, 200] |
In the text |
![]() |
Fig. A.1 Dispersion relation (a, d) and growth rates (b, e) of the extraordinary wave modes according to the kinetic linear theory, with ne/n0 = 24% (upper panels) and 100% (lower panels). The blue and red lines refer to the X and Z modes; the purple and yellow lines refer to the two branches of the relativistic modes (R1 and R2). Panels (c) and (f) present the corresponding intensity ratio of Ex/iEy for each mode. |
In the text |
![]() |
Fig. A.2 Dispersion relation (upper panels) and growth rates (lower panels) of the extraordinary wave modes (X, R1, R2, and Z) with varying ne/n0. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.