On the mechanism of radio emission in type III Solar Radiobursts
Vladimir Krasnoselskikh, Andrii Voshchepynets, Milan Maksimovic

TL;DR
This paper explores how intense density fluctuations in the solar corona influence the mechanism of radio wave generation in type III solar bursts, proposing a more efficient linear energy transfer process than traditional nonlinear models.
Contribution
It demonstrates that density inhomogeneities enable a linear wave energy transformation, challenging the conventional induced scattering mechanism for type III solar radio burst emission.
Findings
Density fluctuations significantly affect beam-plasma interactions.
Reflection of Langmuir waves can partially convert electrostatic energy into electromagnetic.
Linear transformation becomes more efficient than nonlinear scattering at fluctuation levels above 1%.
Abstract
Type III solar radio bursts are generated by streams of energetic electrons accelerated at the Sun during periods of the solar activity. The generation occurs in two steps. Initially, electron beams generate electrostatic Langmuir waves and then these waves are transformed in electromagnetic emissions. It is widely accepted that the mechanism of generation of emission on fundamental frequency close to plasma frequency is due to induced scattering of Langmuir waves into electromagnetic. However this process imposes quite restrictive limit of the ratios of effective brightness temperatures of electromagnetic and Langmuir waves in the source region. Recent studies showed that the level of density fluctuations in the solar wind and in the solar corona is so high that it may significantly affect beam-plasma interaction. Here we show that the presence of intense density fluctuations not only…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsSolar and Space Plasma Dynamics · Ionosphere and magnetosphere dynamics · Earthquake Detection and Analysis
On the mechanism of radio emission in type III solar radio bursts
Vladimir Krasnoselskikh
LPC2E/CNRS, UMR 7328, 3A Avenue de la Recherche Scientifique, Orleans, France
Andrii Voshchepynets
The Swedish Institute of Space Physics (IRF), Rymdcampus 1, Kiruna, Sweden
Milan Maksimovic
LESIA/CNRS, Paris observatory, France
(Received October 29, 2018; Accepted …)
Abstract
Type III solar radio bursts are generated by streams of energetic electrons accelerated at the Sun during periods of the solar activity. The generation occurs in two steps. Initially, electron beams generate electrostatic Langmuir waves and then these waves are transformed in electromagnetic emissions. It is widely accepted that the mechanism of generation of emission on fundamental frequency close to plasma frequency is due to induced scattering of Langmuir waves into electromagnetic. However this process imposes quite restrictive limit of the ratios of effective brightness temperatures of electromagnetic and Langmuir waves in the source region. Recent studies showed that the level of density fluctuations in the solar wind and in the solar corona is so high that it may significantly affect beam-plasma interaction. Here we show that the presence of intense density fluctuations not only crucially influence the process of beam plasma interaction but also changes the mechanism of energy transfer from electrostatic waves into electromagnetic. Reflection of the Langmuir waves from the density inhomogeneities may result in partial transformation of the energy of electrostatic wave into electromagnetic. We show that the linear wave energy transformation for the level of fluctuations of the order of 1% or higher may be significantly more efficient for generation of type III solar radio bursts than conventionally considered process of nonlinear conversion due to induced scattering on ions.
acceleration of particles— Sun: radio radiation— Sun: particle emission— solar wind
††journal: ApJ
1 Introduction
Solar type III radio bursts are amongst the strongest radio emissions in the heliosphere. It is wildly accepted that the high energy electrons keV, accelerated during reconnection of the magnetic flied lines in solar atmosphere, are responsible for generation of these radio emission (Suzuki & Dulk, 1985; Melrose, 2017). Important characteristic of the type III radio bursts is the fast frequency drift rate (Wild & McCready, 1950). Type III bursts can start at a frequency of several hundred of MHz and then go down to tens kHz within few minutes with the increasing duration at a given frequency(Alvarez & Haddock, 1973; Reid & Kontar, 2018; Krupar et al., 2018a). To explain the frequency drift, the beams should have near relativistic speeds (Melrose, 1990) and generate radio emission near plasma frequency, and double plasma frequency (harmonic) (Ginzburg & Zhelezniakov, 1958). Frequency of the emission follows local electron density along the beam path, starting at high frequencies in dense plasma close to the Sun and then decrease over time as the beam propagates in the expanding solar corona and later solar wind (Krupar et al., 2018b).
Generation of radio waves occurs in two steps. In the original study, Ginzburg & Zhelezniakov (1958) proposed that the two-stream instability of an electron beam results in the growth of electrostatic Langmuir waves (ES) that later produce electromagnetic (EM) emission at the plasma frequency due to the scattering on ions while the coalescence of two Langmuir waves can produce the harmonic emission. The theory has been subsequently refined and alternative mechanisms for the conversion of the beam-driven Langmuir waves into electromagnetic radiation have been proposed (Melrose, 1990; Malaspina et al., 2012). While exact mechanism of the ES to EM conversion is still under debate, the generation of the Langmuir waves by electron beams in the solar wind has been confirmed by in situ measurements (Lin et al., 1981; Ergun et al., 1998).
Beam-type configurations in a plasma are known to be unstable and relaxation of the beam-plasma system to the stable state results in growth of the Langmuir turbulence (Romanov & Filippov, 1961; Drummond & Rosenbluth, 1962; Vedenov, 1963). Landau resonance enables effective energy transfer from the beam electrons to the waves, as a result up to of the initial beam energy can be transferred to the Langmuir waves through wave-particle interaction (Vedenov & Ryutov, 1975). Recent studies showed that the level of density fluctuations in the solar wind and in the solar corona is so high that they may significantly affect beam-plasma interaction (Kontar, 2001; Zaslavsky et al., 2010; Krafft et al., 2013; Reid & Kontar, 2013; Voshchepynets et al., 2015). The effect of density fluctuations results in phase velocity variations of Langmuir waves. These variations change the wave resonance velocities of the electrons. This leads to significant decrease of the increment of instability and important increase of the relaxation length of the beam (Voshchepynets & Krasnoselskikh, 2015).
Presence of the density fluctuations has significant impact on the observed properties of the Langmuir waves (see Krasnoselskikh et al. (2007) and references therein). First, this idea has been proposed in Smith & Sime (1979) in order to explain observed clumping of Langmuir waves in type III source regions. Later, it has been proposed (Ergun et al., 2008) that the density irregularities in the solar wind can take form of cavities that might result in modulation of the the waveforms of the Langmuir waves. Analysis of the large number of individual waveforms measured by STEREO and WIND showed good agreement with theoretical predictions (Ergun et al., 2008; Malaspina & Ergun, 2008) and results of the numerical simulations (Krafft et al., 2014). To address stochastic nature of the density fluctuations, several statistical models that deduced properties of the Langmuir waves from the probability distribution function of the density fluctuations have been proposed. Stochastic growth theory proposed by Robinson (1995) predicted that distribution of the amplitudes of the Langmuir waves in type III source regions should follow log-normal distribution. While some observations (Robinson et al., 1993; Píša et al., 2015) show good agreement with these predictions, there are numerous reports of observations and simulations (Krasnoselskikh et al., 2007; Vidojevic et al., 2012; Reid & Kontar, 2017; Voshchepynets et al., 2017) that deviations of the distribution of the amplitudes of the Langmuir waves from log-normal can be rather significant.
The aim of the present study is to determine the role of the density fluctuation in the conversion process of the beam-generated ES waves into EM emissions. Reflection of the Langmuir waves from the density inhomogeneities may result in partial transformation of the energy of electrostatic wave into electromagnetic. We consider this effect of linear wave energy transformation in application to generation of type III solar radio bursts. We use the probability distribution of density fluctuations to evaluate the statistical characteristics of such process and its efficiency. We show that the mechanism of linear transformation for the relative density fluctuations of the order of may be significantly more efficient than widely accepted process of nonlinear conversion of Langmuir waves due to induced scattering on ions.
2 Calculation of the efficiency of energy conversion
When the density fluctuations along the waves path cannot be neglected, the propagation of the wave can be described with non-linear Bohm-Gross dispersion relation for Langmuir waves :
[TABLE]
where and are frequency and wave-vector of the Langmuir wave, is the electron plasma frequency for the electron number density , is a deviation of the density from the average value , and is the Debye length. The waves are assumed to be generated resonantly: , where is the beam velocity that significantly exceeds thermal velocity of electrons ( are electron temperature and mass). A Langmuir wave propagating in plasma with density inhomogeneities encounters density depletions and enhancements along its path. When the wave goes to the increasing density region where the local plasma frequency becomes equal to the wave frequency , the wave is reflected.
Assuming that the density fluctuations are isotropic, the incident angles of Langmuir waves are distributed uniformly over a semi-sphere. For the majority of incidence angles, the reflection resembles mirror type reflection: the -vector component parallel to the direction of the density gradient changes its direction, while the -vector component perpendicular to the gradient remains unchanged. In the rather narrow range of incidence angles, electrostatic waves may couple with the electromagnetic and the process becomes a three-wave coupling process. In this case, the reflection results in a Langmuir wave and an electromagnetic wave, so that a part of the incident Langmuir wave energy is transformed to the electromagnetic wave. It is worth noting that initial reflection generates electromagnetic wave propagating in the direction towards the Sun. However, as an average density decreases with the distance from the Sun the secondary reflections that may be considered as mirror type will turn the wave direction towards the Earth (top panel in Figure 1).
Following the previous works (Piliya, 1966; Stenzel et al., 1974), let the Langmuir wave of frequency propagates obliquely to the direction of the density gradient that we choose to be along the -axes (bottom panel in Figure 1). Let the perpendicular component of the -vector be directed along the -axes and be equal to , so the component along the -axes is
[TABLE]
The angle corresponds to the angle of propagation of the reflected electromagnetic wave with respect to -direction. The incidence angle of electrostatic wave is determined by the ratio of perpendicular and parallel components of -vector, i.e. . When the density profile is linear function of distance with characteristic scale , one finds
[TABLE]
so the incident Langmuir wave reflects when the local plasma frequency is
[TABLE]
where is the speed of light. A fraction of incident wave energy is reflected as Langmuir wave, while the other part is converted to an electromagnetic wave in this mode conversion point . The electromagnetic wave propagates outward into the direction of the density decrease beyond its cutoff frequency at . The problem has been studied by many authors starting with the pioneering work by Denisov (1957). Several methods have been developed to evaluate the conversion coefficient, (e.g. review by Piliya, 1966). Recently, Hinkel-Lipsker et al. (1989, 1991) have performed analytical study and obtained an analytical expression for the conversion coefficient:
[TABLE]
where is the reflection coefficient defined as the ratio of the reflected Langmuir wave amplitude to the amplitude of the incident wave, is the ratio of the electromagnetic and the incident Langmuir wave amplitudes, , and are Airy functions, and and their derivatives. The dependence of the energy conversion coefficient on parameter is shown in Figure 2.
It is convenient to re-write the reflection coefficient in terms of incidence angle of Langmuir wave . The parameter can be written
[TABLE]
where is the characteristic density inhomogeneity scale. Here we take into account that the conversion coefficient has non-zero values only in the range of . In our probabilistic model, the efficiency of the beam-generated Langmuir waves conversion into EM emission is evaluated averaging over angles and the density fluctuation scales. To evaluate the ensemble averaged values taking into account the probability distribution of the density fluctuations, we choose hereafter the reference frame, where -axes is directed along the wave vector of the propagating Langmuir wave. The average probability can be written as the product of probability distributions in angle and in scale :
[TABLE]
where is the energy conversion coefficient from ES Langmuir waves to EM, is energy of reflected EM wave, is the energy of incident Langmuir wave, and is the probability of ES wave reflection. As it was shown in Voshchepynets et al. (2015) and Voshchepynets & Krasnoselskikh (2015), in the plasma with random density fluctuations can be calculated by making use of a probability distribution function of the amplitudes of the density fluctuations .
Since the conversion occurs only when the parameter is in the range from about [math] to , this leads to the limited angular range of reflected electromagnetic waves given by:
[TABLE]
that corresponds to values of the perpendicular component of the -vector of incident Langmuir wave
[TABLE]
Taking into account that the -vector of the Langmuir wave is approximately equal to , the conversion may occur only when the angle given by
[TABLE]
To simplify the evaluation of the integrals, we shall take the conversion coefficient to be approximately constant (corresponding to its average value in the range of between zero and . Under this assumption the integration over angles and scales may be carried out independently step by step. Assuming , the integration over angles results in
[TABLE]
Then the conversion coefficient may be re-written:
[TABLE]
where is a probability distribution of the density gradients (scales). In order to calculate one should use spatial profiles of the density fluctuations. In the present study we use synthetic density data calculated from published density power spectra. Kellogg et al. (1999) proposed a procedure based on the inverse Fourier transform that allows to reconstruct density profiles from the power spectrum assuming the phases of waves to be random. It is known from in situ spacecraft measurements in the solar wind (Celnikier et al., 1987; Chen et al., 2012) that the density spectrum can be considered as a broken power-low, with different spectral indices, about for low frequency part and about for higher frequency part with the transition at about . In order to transform these profile to the spatial profiles , one can use Taylor hypothesis assuming that fluctuations are convected with the characteristic velocity of the solar wind, . We used power spectrum in the range of frequencies between and . Lower frequency limit defines maximal length of the synthetic density profile ( or ). Highest frequency limit defines the smallest scale of the density fluctuation presented in this profile. In the present study this scale is set to be of approximately that for the plasma conditions relevant for 1AU is about . After the profiles were generated, normalization procedure was applied to ensure that and for each of the density profiles (here brackets denote averaging). We consider different levels of density fluctuations, , between and . For more details on the procedure we refer to Voshchepynets & Krasnoselskikh (2015).
Locally the density profiles may be approximated by linear function of , thus the probability distribution of the characteristic scales could be retrieved from the distribution of density gradients, . Left panel in Figure 3 shows normalized probability distribution obtained from the large number (about 200) of the density profiles with level of density fluctuations . It is found that the distribution is very close to Gaussian with characteristic scale :
[TABLE]
Taking into account that and one can get as follows:
[TABLE]
where . The functions and are shown at left panel in Figure 3. By making use of one can integrate last part in equation for energy conversion coefficient as follows:
[TABLE]
where is a Gamma function. Substituting the integral one can find :
[TABLE]
The characteristic scales of the density gradients obtained from synthetic density profiles are shown in Figure 4. We start with that results in . As one can see drops significantly with increasing level of the density fluctuations, that in its turn results in increase of the efficiency of energy transformation from ES to EM waves. Thus for , that was measured onboard Helios satellite closer to the Sun (Bavassano & Bruno, 1995) characteristic scale may be less than .
It is worth noting that the power spectrum used in this study is relevant for solar wind density fluctuations around 1AU. Closer to the Sun the spectrum characteristics may be different and as a result density fluctuations can be described by the different statistics. In order to avoid speculations (though the method developped here is applicable), we consider emissions in the frequency range typical for solar type III radio burst around 1 AU: to (Mann et al., 1999).
Conversion coefficient as a function of beam velocity, and Langmuir wave frequency is shown in Figure 5. Left panel shows for . As one can see for conversion coefficient is above in the whole range of frequencies. An increase of the level of the density fluctuations results in a decrease of the characteristic scale of the density gradient. As a result, reflection of the Langmuir waves will occur more often and will increase. Right panel in Figure 5 shows for . We found that the conversion coefficient is above in the whole frequency range for . For faster beams with conversion coefficient is above for frequencies below .
3 Induced scattering
It is well known and widely accepted from the very early articles by Ginzburg & Zhelezniakov (1958), Kaplan & Tsytovich (1969), Melrose (1974) and Melrose (1987) that under condition (typically in the solar wind the damping of ion sound waves is quite strong and the major nonlinear process of transformation of electrostatic waves onto electromagnetic (ES-EM) is induced scattering on ions. Generated electromagnetic waves have frequencies lower than the frequency of the primary Langmuir waves. To describe this process one should introduce emission coefficient for electromagnetic waves that is defined as the energy per unit frequency interval generated in a unit volume and in a unit solid angle. This coefficient is expressed making use a power, , radiated by current excited by Langmuir waves.
For the conversion of plasma waves into electromagnetic due to induced scattering on ions, can be found as follows (Kaplan & Tsytovich, 1969):
[TABLE]
For the sake of simplicity we assume . One can also assume with the accuracy up to coefficient of the order of unity that radiation is isotropic (implies that ) and . This results in:
[TABLE]
Suggesting that for the beam plasma interaction the maximum of the energy density of the plasma turbulence occurs at a definite value of the phase velocity , and neglecting the width of the beam choosing (here , and is along the direction of the propagation of the beam). This results in total emission coefficient that is found integrating over and :
[TABLE]
The resonance conditions for the wave conversion impose the following relation between wave vectors of electromagnetic and Langmuir waves: and . Thus the one can find:
[TABLE]
One can use total emission coefficient to calculate averaged over period wave energy density of the EM emission, . For the comparison with the linear transformation coefficient determined as for ES-EM transformation due to induced scattering on thermal ions it is found to be equal to:
[TABLE]
One can see that for typical parameters of the solar wind it is many orders of magnitude smaller than the efficiency coefficient for linear transformation. For the conditions at 1 AU , , .
3.1 Conclusions
We show that the process of linear conversion of Langmuir waves onto electromagnetic on density fluctuations can be dominant for the generation of the type III radio emissions and is significantly more efficient than conventionally accepted nonlinear process of induced scattering of Langmuir waves on ions.
As we show the efficiency of linear conversion is strongly dependent upon statistical properties of density fluctuations and their gradients. These characteristics may significantly vary with the distance from the Sun. The study of these dependencies comes beyond the scope of this paper and will be addressed in future publications.
V.K. acknowledges financial support by CNES through grant ”Stereo-Waves invited scientist”.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Alvarez & Haddock (1973) Alvarez, H., & Haddock, F. T. 1973, Sol. Phys., 30, 175, doi: 10.1007/BF 00156186 · doi ↗
- 2Bavassano & Bruno (1995) Bavassano, B., & Bruno, R. 1973, Sol. Phys., 100, 9475, doi: 10.1029/94JA 03048 · doi ↗
- 3Celnikier et al. (1987) Celnikier, L. M., Muschietti, L., & Goldman, M. V. 1987, A&A, 181, 138
- 4Chen et al. (2012) Chen, C. H. K., Salem, C. S., Bonnell, J. W., Mozer, F. S., & Bale, S. D. 2012, Physical Review Letters, 109, 035001, doi: 10.1103/Phys Rev Lett.109.035001 · doi ↗
- 5Denisov (1957) Denisov, N. G. 1957, Soviet Physics JETP, 4, 544
- 6Drummond & Rosenbluth (1962) Drummond, W. E., & Rosenbluth, M. N. 1962, Physics of Fluids, 5, 1507, doi: 10.1063/1.1706559 · doi ↗
- 7Ergun et al. (1998) Ergun, R. E., Larson, D., Lin, R. P., et al. 1998, Ap J, 503, 435, doi: 10.1086/305954 · doi ↗
- 8Ergun et al. (2008) Ergun, R. E., Malaspina, D. M., Cairns, I. H., et al. 2008, Physical Review Letters, 101, 051101, doi: 10.1103/Phys Rev Lett.101.051101 · doi ↗
