Issue 
A&A
Volume 586, February 2016



Article Number  A19  
Number of page(s)  10  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201527271  
Published online  21 January 2016 
Twodimensional time evolution of beamplasma instability in the presence of binary collisions
^{1} Instituto de Física, Universidade Federal do Rio Grande do Sul, 91501970 Porto Alegre, RS, Brazil
^{2} Institute for Physical Science & Technology, University of Maryland, College Park, MD 20742, USA
^{3} School of Space Research, Kyung Hee University, Yongin, 446701 Gyeonggi, South Korea
^{4} School of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, UK
Received: 29 August 2015
Accepted: 21 November 2015
Energetic electrons produced during solar flares are known to be responsible for generating solar type III radio bursts. The radio emission is a byproduct of Langmuir wave generation via beamplasma interaction and nonlinear wavewave and waveparticle interaction processes. In addition to type III radio bursts, electrons traveling downwards toward the chromosphere lead to the hard Xray emission via electronion collisions. Recently, the role of Langmuir waves on the Xrayproducing electrons has been identified as important, because Langmuir waves may alter the electron distribution, thereby affecting the Xray profile. Both Coulomb collisions and waveparticle interactions lead electrons to scattering and energy exchange that necessitates considering the twodimensional (2D) problem in velocity space. The present paper investigates the influence of binary collisions on the beamplasma instability development in 2D in order to elucidate the nonlinear dynamics of Langmuir waves and binary collisions. The significance of the present findings in the context of solar physics is discussed.
Key words: Sun: radio radiation / Sun: particle emission / methods: numerical
© ESO, 2016
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The beamplasma interaction between the energetic electrons produced during the solar flare and the background coronal plasma is a source of rich nonlinear phenomena. A number of solar radio bursts are associated with nonthermal electron populations and associated Langmuir waves. One of the most important and wellstudied phenomena is the generation of solar typeIII solar radio bursts (Wild & McCready 1950; Wild 1950). Fast electrons escaping from active regions in the Sun to the corona and interplanetary space generate Langmuir turbulence, which partially converts to radiation at the local plasma frequency and/or its harmonic(s). A nonlinear conversion process involves sophisticated wave decay and wave scattering. This is the wellknown plasma emission, and it is the basic radio emission mechanism for type III radio bursts and for the reverse slope bursts (see, e.g., Tang & Moore 1982; Dennis et al. 1984).
The first theory of plasma emission was put forth by Ginzburg & Zhelezniakov (1958) and many modifications and improvements have been made over the past six decades (Tsytovich 1967; Kaplan & Tsytovich 1968; Zheleznyakov & Zaitsev 1970; Melrose 1982; Goldman & Dubois 1982; Goldman 1983; Cairns 1987; Robinson & Cairns 1998; Mel’Nik et al. 1999; Kontar 2001; Ledenev et al. 2004; Li et al. 2005, 2006a,b, 2008a,b; Li & Cairns 2013; Ratcliffe & Kontar 2014). Although the essential theoretical framework based upon EM weak turbulence theory, which describes the entire process starting from the beamgenerated Langmuir turbulence to the radiation generation, was available, complete numerical solution of the entire set of EM weak turbulence equations have not been done until quite recently, when Ziebell et al. (2014a,b,c, 2015) numerically solved the complete equations of EM weak turbulence theory for the first time. It should be mentioned that a few authors carried out direct EM particleincell (PIC) simulations to characterize the nonlinear behavior of the plasma emission process (Kasaba et al. 2001; Karlický & Vandas 2007; Rhee et al. 2009a,b; Umeda 2010; Ganse et al. 2012a,b, 2014). The full numerical solution of the analytical EM weak turbulence equations by Ziebell et al. (2014a, 2015) complements these PIC simulation efforts.
The type III radio bursts are not the only process of importance associated with solar flares. For electrons traveling down to the chromosphere, they generate Xrays via bremsstrahlung when they collide with plasma. The approximate treatment of the electron dynamics by only considering the Coulomb collisions is known as the thick target model (e.g., Holman et al. 2011, as a recent review). Under such a simplifying assumption, Brown (1971) and Syrovatskii & Shmeleva (1972) analyzed the dynamics of electron distribution and the related bremsstrahlung Xray spectrum from accelerated/injected electrons. In the literature, however, there are discussions of the importance of Langmuir wave generation on the underlying nonthermal electron energy distribution. Emslie & Smith (1984) first considered the effects of wave generation on the electron beam propagation toward the chromosphere. Hamilton & Petrosian (1987) and McClements (1987, 1989) reconsidered the same problem but found that the influence of waveparticle interactions was insignificant for stationary electron injection. However, recent works by Hannah et al. (2009, 2013) and Karlický & Kontar (2012) demonstrate that the collective effects are sufficiently significant after all, especially for transient sources of nonthermal electrons, such that without these effects, the solar hard Xray spectra measured by RHESSI spacecraft may not be interpreted correctly. Zharkova & Siversky (2011) also studied a similar problem. Moreover, Kontar et al. (2012) have shown that the presence of plasma inhomogeneity and/or nonlinear wavewave interactions can lead to an overestimated number and energy of energetic electrons accelerated in the corona even in the case of stationary electron injection.
As noted above, the problem of flaregenerated nonthermal electrons and Xray emission from the same source that also generates type III radio bursts is a complicated one, which requires simultaneous timedependent treatment of various processes that include Coulomb collisions, collective Langmuir wave generation, and nonlinear wavewave interactions. In addition, the evolution of electrons not only in energy but also in collisional pitchangle scattering is important for interpreting some Xray observations (e.g., MacKinnon 1991; Jeffrey et al. 2014).
The purpose of the present paper is not to revisit the plasma emission, which has recently been discussed by Ziebell et al. (2014a, 2015) and by Ratcliffe & Kontar (2014) in application to flaring loops, or to address the issue of the influence of collective processes on the Xray generation by bremsstrahlung per se, which was addressed by Hannah et al. (2009) and Kontar et al. (2012). Instead, we address a fundamental 2D plasma physics problem that may be related to both the typeIIIgenerating electrons and the Xraygenerating nonthermal electrons. That is, we consider the influence of both Coulomb collisions and the collective effects in twodimensional (2D) electron beamplasma interaction problem, taking only the evolution of the particle distribution and of the spectra of electrostatic waves into account, without incorporating the effects associated to electromagnetic oscillations.
Unlike previous works related to this type of analysis where the effects of Coulomb collisions have normally been ignored (Ziebell et al. 2008), the present paper includes electronelectron and electronion collisions. On the other hand, unlike the previous works related to the solar Xray problem where, in addition to collisional effects, quasilinear effects are included considering onedimensional (1D) wave excitation, we now consider the nonlinear dynamics associated with the wave scattering and decay in two dimensions, along with the collisional dynamics. The 2D evolution is particularly important since collisional scattering and collisional energy loss operate on the same time scales. Moreover, radioemission and Langmuir wave evolution depend on angular distribution of plasma waves. The objective with such relatively limited analysis that only includes electrostatic waves is to gather information about the effect of collisions on the time evolution of the beamplasma system, on the relative roles of collisions and nonlinear mechanisms of threewave decay and waveparticle scattering, and on the possibly different time scales of the different physical processes when evolving selfconsistently. The information to be obtained with such study may be useful for future analysis of the time evolution of more complex processes, in which other nonlinear mechanisms have to be taken into account, such as the actual production of fundamental and harmonic emission that characterize type III emissions and the complicated interaction between energetic electrons, background plasma, and waves, which is associated to the solar Xray emission.
The structure of the paper is as follows. In Sect. 2 we briefly describe the theoretical formulation and the setup for the numerical analysis. In Sect. 3 we present and discuss the results of the numerical analysis. Finally, in Sect. 4 we conclude the paper and comment on the results obtained.
2. Theoretical formulation and numerical setup
The wave kinetic equations for L and S waves that describe quasilinear process as well as nonlinear decay and scattering processes are given in terms of the spectral wave energy density, and , where and represent the spectral electric field component associated with L and S waves, respectively, and where σ = ± 1 stands for the sign of wave phase velocity. The wave kinetic equations for these waves are given by Yoon (2006) and Ziebell et al. (2008): where ω_{p} = (4πn_{0}e^{2}/m_{e})^{1/2} is the electron plasma frequency, and e, m_{e}, and n_{0} stand for the unit electric charge, electron mass, and the ambient particle number density, respectively. The dispersion relations for L and S modes are wellknown: where v_{e} = (2T_{e}/m_{e})^{1/2} is the electron thermal speed, and is the electron Debye length, with T_{e} being the electron temperature. In (1) and (2) we have also introduced a quantity The first terms on the righthand sides of (1) and (2) describe the spontaneous emission and quasilinear effects. The second terms contain the energy conservation condition, for L mode and a similar threewave resonance condition for S mode, and describe the threewave decay processes. The third term in (1) contains the nonlinear waveparticle resonance condition , and it depicts the scattering of L waves by the ions. We have neglected the scattering effects for (2), which rules the evolution of S waves, since the scattering processes involving S waves are extremely slow processes.
For the time evolution of the particle distribution, we consider the effects of spontaneous fluctuation, the quasilinear diffusion effect, and the effect of binary collisions. The particle kinetic equation is thus given by
(3)The term with coefficient A_{i} describes the effects of spontaneous fluctuations, and the term with coefficient D_{ij} rules the quasilinear diffusion process. The symbol α represents either L or S, depending on whether the particle species a in the above particle kinetic equations stands for a = e (electrons) or a = i (ions). The quantity θ_{ab}(F_{a},F_{b}) represents the binary collisions of particles of species a with particles of species b. The quantity lnΛ is the Coulomb logarithm.
In the set of wave kinetic Eqs. (1) and (2), we do not (yet) include a collisional damping term. The collisional effect is only incorporated through the particle kinetic equation (3), which contains the collision operator ∑ _{b}θ_{ab}(F_{a},F_{b}). Upon linearization of the governing particle kinetic equation and coupling to the wave equation, then the collision frequency ν_{coll} broadens (Dupree 1966; Ishihara & Hirose 1985; Bian et al. 2014; Pécseli 2014) the waveparticle resonance conditions, ω−k·v, where ν_{coll} is the effective collision frequency that is obtained from the particle equation. In this way, the collisional damping rate for the collective wave phenomena is calculated by the particle collision term in an indirect way. However, a more complete treatment should contain explicit collisional damping term for the wave kinetic equation. In a recently submitted work, Yoon et al. (2015) address this issue by reformulating the standard weak turbulence theory to include contributions from those electrostatic fluctuation spectra (ω,k) for which ω and k do not satisfy the plasma wave dispersion relations. They show that these socalled “noneigenmode” contributions, which are typically ignored in the literature, are not only responsible for the same collision integral ∑ _{b}θ_{ab}(F_{a},F_{b}) in the particle kinetic equation, as in Eq. (3), but also they lead to the collisional damping term, as well as the electrostatic analog of the bremsstrahlung. One could, however, include the collisional damping term heuristically, as in the recent paper by Ratcliffe & Kontar (2014).
We introduce the following dimensionless variables and definitions: In terms of the normalized variables and quantities, the equation for Langmuir (L) waves is expressed as (4)The equation for ionacoustic (S) waves is likewise expressed in terms of dimensionless quantities, (5)The dimensionless particle kinetic equation is given by (6)Finally, the dispersion relations in dimensionless form are given by where v_{te} = (2T_{e}/m_{e})^{1/2} is the electron thermal speed, and λ_{De} is the electron Debye length, with T_{e} being the electron temperature.
For further reference, in the equation for the L waves, Eq. (5), the term that describes the spontaneous emission and the quasilinear effects has been denoted Lql, the term that describes the threewave decay involving L and S waves has been denoted LdLS, and the term describing scattering involving L waves has been denoted LsLL. In the equation for the evolution of S waves, namely Eq. (6), the first term describes the spontaneous emission and quasilinear effects and is denoted Sql, and the second term describes the threewave decay involving L waves and is denoted SdLL.
In the equation for the particle distribution functions (7), the term with g describes the effects of spontaneous fluctuations, and the term with the velocity derivative describes the quasilinear diffusion process.
For a detailed derivation of the above equations without the term that represents the collisional effects, the reader is referred to Yoon (2000, 2005, 2006). The collisional term is essentially the Landau collision integral, which is well known in the literature (e.g., Lifshitz & Pitaevskii 1981; Karney 1986; Jeffrey et al. 2014, see also Appendix A). Assuming that both ion and electron background distributions are Maxwellian distributions, the kinetic equation describing the evolution of electron distribution function F_{a} in fully ionized hydrogen plasma can be written in spherical coordinates (v,θ), (7)where, μ = v_{∥}/v = cos(θ), , , erf(u) is the error function, and G(u) is the Chandrasekhar function, given by (8)The first term on the righthand side of Eq. (7) describes electronelectron collisions, while the second is due to electron scattering by ions. One can combine the collisional terms (Eq. (7)) with the collective terms (Eq. (6)) and define generalized friction and diffusion coefficients.
The initial electron distribution function is assumed to be made of a Maxwellian background and a forward beam distribution with number density n_{f}. In 2D, the assumed distribution function is given as (9)Here and are the background and the beam thermal velocities squared, respectively. The quantities v_{0} and v_{f} are the drift velocities of the background and forward beam, respectively. The background drift velocity v_{0} is such that it guarantees zero net drift velocity for the electron distribution, i.e., v_{0} = −(v_{f}n_{f})/(n_{0}−n_{f}). The initial ion distribution in 2D is given by F_{i}(v) = (m_{i}/ 2πT_{i}) exp^{[}−m_{i}v^{2}/ (2T_{i})^{]}, where T_{i} and m_{i} are the ion temperature and the proton mass, respectively. For the present study, we assume that the ion distribution remains constant along the time evolution.
The intensities of L and S waves are initialized by balancing the spontaneous and induced emissions, taking the background population into account: (10)
3. Numerical analysis
The set of Eqs. (4), (5), and (8)–(10) have been solved in 2D wave number space and 2D velocity space by employing a splitting method with fixed time step for the equation for evolution of the electron distribution and a fourthorder RungeKutta method for the wave equations. The 2D version for the particle equation is shown by Eq. (A.11), with the coefficients coming from weak turbulence theory as given by (A.12), and the collision term given by Eqs. (A.9) and (A.10). Similar equations for the wave evolution can be obtained by converting Eqs. (4) and (5) to 2D coordinates. For all the cases, we have used Δτ = 0.1. We employed a 51×51 grid for k_{⊥} and k_{∥} with 0 <k_{⊥}v_{te}/ω_{p}< 0.6, and 0 <k_{∥}v_{te}/ω_{p}< 0.6. For the velocities, we used a 51 × 101 grid for v_{⊥}/v_{te} and v_{∥}/v_{te}, covering the velocity range 0 <v_{⊥}/v_{te}< 12 and − 12 <v_{∥}/v_{te}< 12. For the subsequent numerical solutions to be discussed, we assume a forward beam with the same temperature as the background population, T_{f}/T_{e} = 1.0, electrontoion temperature ratio T_{e}/T_{i} = 7.0, plasma parameter , and electron temperature, such that . We also assume v_{f}/v_{te} = 5.0 and n_{f}/n_{0} = 1.0 × 10^{3}, a value that is relatively high for the case of solar wind beam electrons, but has been assumed in order to reduce the computational time required to obtain the numerical results.
The characteristics of the grids used for the numerical analyses were also chosen as a compromise between the resolution obtained and the performance of the numerical code. The grid used in wave number space resolves some fundamental details in the wave spectra, such as (a) the primary peak that grows due to quasilinear interaction between L waves and the beam particles in the region of the positive derivative of the distribution function and (b) the peak corresponding to backwardly propagating L waves, which are generated by nonlinear interaction between L and S waves. In velocity space, the grid used allows for a smooth description of the distribution function and of the beam, including the important region of transition between the core distribution and the beam. Even with use of the chosen grids, the code needs to run for a few days in a personal computer with updated technology. Owing to the 2D nature of the problem investigated, a twofold increase in resolution would require nearly a fourfold increase in computational time. We have reached the conclusion that such a costly enhanced resolution would not be needed to describe the relevant features of the time evolution of the system under investigation.
As one application of the numerical code, we studied the time evolution of the system in two situations. The first situation neglects the influence of collisions, which can be easily done just by turning off the collisional term in the equation for the evolution of the electron distribution function. In the second situation, we took the effect of binary collisions into account. For comparison, we also added some results that show the time evolution of the distribution function that is only subject to collisions.
Fig. 1
Time evolution of electron velocity distribution function is shown in 2D velocity space plot at different time steps. The lefthand panels depict solutions with only collisions (as in the thick target approximation), the middle panels correspond to the selfconsistent solution without collisional effects (as in the type III beam case), and the right panels show solutions that contain both collisional dynamics and selfconsistent Langmuir wave dynamics. 
Figure 1 shows the time evolution of electron velocity distribution function by plotting F_{e} in 2D velocity space at different time steps corresponding to ω_{pe}t = 500, 1000, 2000, 5000, 10^{4}, 1.5 × 10^{4}, and 2 × 10^{4}. The lefthand panels depict solutions with only collisions; the middle panels correspond to the selfconsistent solution without collisional effects, as in the type III beam case (Ziebell et al. 2008); and the righthand panels show solutions that contain both collisional effects and selfconsistent Langmuir wave dynamics. As the lefthand panels show, the collisional relaxation of the beam leads to the systematic energy losses and isotropization of the beam. Eventually, the beam merges with the Maxwellian background. In contrast, collisionless quasilinear relaxation, appearing in the middle panels, shows that slowing down the beam proceeds via the velocity space plateau formation.
The total electron distribution undergoes some heating, but the overall shape of the distribution is quite different from the collisiononly case. In these middle panels, it is noticeable that an Xshaped band of irregular features appears near the core of the distribution in the late stages of the time evolution. These irregular features in the distribution function are apparently connected with irregular features that appear near the edge of the grid representing the spectrum of L waves and must be due to accumulation of numerical imprecision. These irregular features do not spread and do not produce any perceptible effect on the relevant features of the distribution, such as the formation of the plateau joining the core and the beam distributions or the overall broadening of the velocity distribution.
In accordance with this interpretation, one sees that the irregular features do not appear in the left and righthand panels, since they are washed out by the effect of collisions. At the end of the computational period, it is seen that the net electron distribution is still anisotropic, with the 2D distribution elongated along the beam direction. In the case of collisions plus the Langmuir wave dynamics (the righthand panels), the early evolution is similar to the purely quasilinear case, but the later evolution is affected by the Coulomb collisional dynamics. Such a combined beamplasma dynamical interaction showing the transition from the collective dynamics to collisional dynamics has not been done in the literature.
In the context of the solar applications, the third panels are similar to those considered (e.g., Hamilton & Petrosian 1987; McClements 1989; Hannah et al. 2009; Hannah & Kontar 2011; Zharkova & Siversky 2011), in that both the Coulomb collisional dynamics and selfconsistent Langmuir wave dynamics are considered. The above references, however, also heuristically included collisional damping term for the Langmuir waves, whereas in the present approach, the collisional effects influence the Langmuir dynamics indirectly and only through the electron distribution function evolution. In spite of this caveat, however, our work is more general in the sense that we consider twodimensional dynamics and also allow nonlinear wave processes, such as the decay interaction and nonlinear waveparticle scattering similar to Kontar et al. (2012) and Karlický & Kontar (2012). Also, the above references employed simplified Coulomb collisional integral for the particles, whereas our numerical solution is based upon the Landau operator. In the context of solar Xray, however, inhomogeneity effects is important. In addition, for the solar Xray problem different initial electron velocity distributions should be adopted, such as the powerlaw distribution. We do not have inhomogeneity effects (see, e.g., Kontar & Pécseli 2002), and our choice of initial electron distribution is a simple Gaussian core plus a drifting beam. As such, our calculation is not meant to represent the solar Xray problem. Instead, our interest has been on understanding the 2D evolution.
Fig. 2
Time evolution of Langmuir wave spectrum as a 2D wave number space plot at different time steps. The lefthand panels depict solutions with only selfconsistent collective effects without the collisional terms in the particle equation, and the right ones show solutions that contain both collisional dynamics and selfconsistent wave dynamics. 
Figure 2 shows the spectrum of Langmuir (L) waves, as a function of normalized components of a wave vector in 2D space for several values of the normalized time between τ = 500 and τ = 20 000. Since the purely collisional dynamics involves no waves by default, the Langmuir spectrum is absent in that case. As such, Fig. 2 only shows two situations. The lefthand panels show Langmuir wave spectra when only selfconsistent collective effects are considered without the collisional terms in the particle equation, and the righthand panels solutions that contain both collisional dynamics and selfconsistent wave dynamics in the particle kinetic equations. Of course, for both cases, the ionsound wave intensity is also selfconsistently calculated, but the results are not shown. Both cases exhibit similar earlytime dynamics of the Langmuir wave time evolution.
The forwardpropagating “primary” Langmuir waves are the direct result of bumpontail instability by the electron beam. For ω_{pe}t = 500, the nonlinear decay and scattering off thermal ions (these are nonlinear processes depicted in the wave kinetic equations) already begin to produce weak backwardtraveling Langmuir waves. These “backscattered” Langmuir waves form a semiarc shaped spectra as time proceeds, eventually distributing the wave momentum in the quasi circular area in 2D wave number space. For later dynamics, such as beyond ω_{pe} = 15 000, however, the evolution of the wave spectra for two cases are quite distinct. In the pure collective dynamic case, one notices that there is an appreciable gap along the parallel wave number; in contrast, in the general case where both collective dynamics and collisions are included, the primary Langmuir waves undergo spreading at a wave vector azimuthal angle so that the wave spectrum appears to be much more symmetric throughout the semicircular area.
In the context of solar applications, Hamilton & Petrosian (1987), McClements (1989), Hannah et al. (2009), Hannah & Kontar (2011), and Zharkova & Siversky (2011), solved selfconsistent Langmuir wave kinetic equations. Kontar et al. (2012) also demonstrate numerically that nonlinear wavewave interactions are important. However, since these authors did not consider nonlinear mode coupling physics in 2D, the Langmuir spectrum in such a scheme does not give the wave vector pitch angle evolution. The evolution in pitch angle affects the level of Langmuir waves and generally reduces peak values of spectral energy density of plasma waves.
4. Final remarks
The energetic electrons produced during the solar flare lead to the generation of solar type III solar radio bursts and/or reverse slope bursts when the fast electrons escape the acceleration region. Electrons propagating downwards are responsible for reverse slope bursts (e.g., Tang & Moore 1982) that generate Langmuir waves in collisional plasma. These electrons traveling down towards the chromosphere are responsible for Xrays emitted via bremsstrahlung (see, e.g., Holman et al. 2011). There have been recent discussions (Hannah et al. 2009; Hannah & Kontar 2011; Zharkova & Siversky 2011; Kontar et al. 2012) that show that the collective Langmuirwave dynamics of the Xray generating electrons may be important for interpretations of Xray spectra from RHESSI spacecraft (Lin et al. 2002).
In the present paper, we have addressed a fundamental plasma physics problem that is related to both the type IIIgenerating electrons and the Xraygenerating nonthermal electrons. We considered the influence of Coulomb collisions and the selfconsistent collective effects in the electron beamplasma interaction problem in 2D velocity space. In the past, the present authors investigated the electronLangmuir wave collective dynamics, which is relevant to the type III burst electrons, by ignoring the effects of Coulomb collisions (Ziebell et al. 2008). On the other hand, early works on solar Xray emission ignored the collective dynamics and only considered collisional dynamics. Later works include waveparticle interactions and included Langmuir wave generation, but they only considered the 1D problem, although they included other effects relevant to the solar Xray problem, such as the inhomogeneity and phenomenological collisional damping. In contrast, we considered the full nonlinear dynamics in the wave kinetic equation without explicit collisional damping (yet), but which included nonlinear waveparticle scattering and threewave decay interactions in 2D wave number space.
Upon comparing the electron distribution evolution among three situations, namely, one in which only the collisional effects are included, another with only collective dynamics, and the third in which both collisions and waves are incorporated, we found that the early dynamics in the general case follows the collective dynamical pattern, whereas for later times, the dynamics closely matches the collisional case. This transitional behavior has not been demonstrated clearly in the literature (Fig. 1). The comparison between results obtained considering the collisionless weak turbulent dynamics and the collisional weak turbulent dynamics shows that for τ ≃ 2000, a normalized time for which the primary Langmuir peak excited by the beam is fully developed and the backward propagating Langmuir peak is already well established, the effect of collisions is already noticeable, leading to a wider plateau in the distribution function and increased tendency to isotropization, as compared to the collisionless case. For a much longer evolution time, the results show that the distribution function obtained in the collisional case is clearly more isotropic than was obtained in the case without the effect of collisions and also show a noticeable decrease in the energetic tail. Considering the effect of the wave spectra as well, the general comment is that collisions lead to a wider angular spread of both Langmuir waves and electron beam on comparable time scales.
When we compared the wave dynamics between the pure collective case and the general case (for collision only case there is no wave), we found that the late time evolution of the wave spectra also shows a marked difference between the two cases. Specifically, the general case in which both collisions and collective effects were considered shows more symmetric Langmuir wave spectrum at late time periods. The overall conclusion to be drawn from both analyses of the time evolution of the velocity distribution function and of the spectrum of Langmuir waves is that the collisional processes, even if irrelevant on the time scale of the development of quasilinear processes, may lead to effects that are comparable in magnitude to effects of the nonlinear processes on a time scale that is compatible to the time scale of evolution of nonlinear processes in space plasmas.
As noted, the actual application of the present analysis to the solar Xray problem requires additional effects that are not included here: inhomogeneous density profile, explicit collisional damping terms in the wave kinetic equations, and probably different initial nonthermal electron distribution functions. However, we believe that future applications can be built upon the present formalism.
Acknowledgments
This work was partially supported by the Brazilian agencies CNPq and FAPERGS. P.H.Y. acknowledges support by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea, and NSF Grant Nos. AGS1138720 and AGS1242331. E.P.K. was supported by an STFC consolidated grant.
References
 Bian, N. H., Kontar, E. P., & Ratcliffe, H. 2014, J. Geophys. Res. (Space Phys.), 119, 4239 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, J. C. 1971, Sol. Phys., 18, 489 [NASA ADS] [CrossRef] [Google Scholar]
 Cairns, I. H. 1987, J. Plasma Phys., 38, 169 [NASA ADS] [CrossRef] [Google Scholar]
 Dennis, B. R., Benz, A. O., Ranieri, M., & Simnett, G. M. 1984, Sol. Phys., 90, 383 [NASA ADS] [CrossRef] [Google Scholar]
 Dupree, T. H. 1966, Phys. Fluids, 9, 1773 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Emslie, A. G., & Smith, D. F. 1984, ApJ, 279, 882 [NASA ADS] [CrossRef] [Google Scholar]
 Gaffey, Jr., J. D. 1976, J. Plasma Phys., 16, 149 [NASA ADS] [CrossRef] [Google Scholar]
 Ganse, U., Kilian, P., Spanier, F., & Vainio, R. 2012a, ApJ, 751, 145 [NASA ADS] [CrossRef] [Google Scholar]
 Ganse, U., Kilian, P., Vainio, R., & Spanier, F. 2012b, Sol. Phys., 280, 551 [NASA ADS] [CrossRef] [Google Scholar]
 Ganse, U., Kilian, P., Spanier, F., & Vainio, R. 2014, A&A, 564, A15 [Google Scholar]
 Ginzburg, V. L., & Zhelezniakov, V. V. 1958, Sov. Astron., 2, 653 [Google Scholar]
 Goldman, M. V. 1983, Sol. Phys., 89, 403 [NASA ADS] [CrossRef] [Google Scholar]
 Goldman, M. V., & Dubois, D. F. 1982, Phys. Fluids, 25, 1062 [NASA ADS] [CrossRef] [Google Scholar]
 Hamilton, R. J., & Petrosian, V. 1987, ApJ, 321, 721 [NASA ADS] [CrossRef] [Google Scholar]
 Hannah, I. G., & Kontar, E. P. 2011, A&A, 529, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hannah, I. G., Kontar, E. P., & Sirenko, O. K. 2009, ApJ, 707, L45 [NASA ADS] [CrossRef] [Google Scholar]
 Hannah, I. G., Kontar, E. P., & Reid, H. A. S. 2013, A&A, 550, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Holman, G. D., Aschwanden, M. J., Aurass, H., et al. 2011, Space Sci. Rev., 159, 107 [NASA ADS] [CrossRef] [Google Scholar]
 Ishihara, O., & Hirose, A. 1985, Phys. Fluids, 28, 2159 [NASA ADS] [CrossRef] [Google Scholar]
 Jeffrey, N. L. S., Kontar, E. P., Bian, N. H., & Emslie, A. G. 2014, ApJ, 787, 86 [NASA ADS] [CrossRef] [Google Scholar]
 Kaplan, S. A., & Tsytovich, V. N. 1968, Sov. Astron., 11, 956 [NASA ADS] [Google Scholar]
 Karlický, M., & Kontar, E. P. 2012, A&A, 544, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Karlický, M., & Vandas, M. 2007, Planet. Space Sci., 55, 2336 [NASA ADS] [CrossRef] [Google Scholar]
 Karney, C. F. F. 1986, Comp. Phys. Rep., 4, 183 [Google Scholar]
 Kasaba, Y., Matsumoto, H., & Omura, Y. 2001, J. Geophys. Res., 106, 18693 [NASA ADS] [CrossRef] [Google Scholar]
 Kontar, E. P. 2001, Sol. Phys., 202, 131 [NASA ADS] [CrossRef] [Google Scholar]
 Kontar, E. P., & Pécseli, H. L. 2002, Phys. Rev. E, 65, 066408 [NASA ADS] [CrossRef] [Google Scholar]
 Kontar, E. P., Ratcliffe, H., & Bian, N. H. 2012, A&A, 539, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ledenev, V. G., Zverev, E. A., & Starygin, A. P. 2004, Sol. Phys., 222, 299 [NASA ADS] [CrossRef] [Google Scholar]
 Li, B., & Cairns, I. H. 2013, J. Geophys. Res. (Space Phys.), 118, 4748 [Google Scholar]
 Li, B., Willes, A. J., Robinson, P. A., & Cairns, I. H. 2005, Phys. Plasmas, 12, 012103.1 [NASA ADS] [CrossRef] [Google Scholar]
 Li, B., Robinson, P. A., & Cairns, I. H. 2006a, Phys. Plasmas, 13, 092902 [NASA ADS] [CrossRef] [Google Scholar]
 Li, B., Robinson, P. A., & Cairns, I. H. 2006b, Phys. Rev. Lett., 96, 145005 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Li, B., Cairns, I. H., & Robinson, P. A. 2008a, J. Geophys. Res. (Space Phys.), 113, 6104 [NASA ADS] [CrossRef] [Google Scholar]
 Li, B., Cairns, I. H., & Robinson, P. A. 2008b, J. Geophys. Res. (Space Phys.), 113, 6105 [Google Scholar]
 Lifshitz, E. M., & Pitaevskii, L. P. 1981, Physical kinetics (Oxford: Pergamon Press) [Google Scholar]
 Lin, R. P., Dennis, B. R., Hurford, G. J., et al. 2002, Sol. Phys., 210, 3 [Google Scholar]
 MacKinnon, A. L. 1991, A&A, 242, 256 [NASA ADS] [Google Scholar]
 McClements, K. G. 1987, A&A, 175, 255 [NASA ADS] [Google Scholar]
 McClements, K. G. 1989, A&A, 208, 279 [NASA ADS] [Google Scholar]
 Mel’Nik, V. N., Lapshin, V., & Kontar, E. 1999, Sol. Phys., 184, 353 [NASA ADS] [CrossRef] [Google Scholar]
 Melrose, D. B. 1982, Sol. Phys., 79, 173 [NASA ADS] [CrossRef] [Google Scholar]
 Pécseli, H. L. 2014, J. Plasma Phys., 80, 745 [NASA ADS] [CrossRef] [Google Scholar]
 Ratcliffe, H., & Kontar, E. P. 2014, A&A, 562, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rhee, T., Ryu, C.M., Woo, M., et al. 2009a, ApJ, 694, 618 [NASA ADS] [CrossRef] [Google Scholar]
 Rhee, T., Woo, M., & Ryu, C.M. 2009b, J. Kor. Phys. Soc., 54, 313 [NASA ADS] [CrossRef] [Google Scholar]
 Robinson, P. A., & Cairns, I. H. 1998, Sol. Phys., 181, 363 [NASA ADS] [CrossRef] [Google Scholar]
 Syrovatskii, S. I., & Shmeleva, O. P. 1972, Sov. Astron., 16, 273 [NASA ADS] [Google Scholar]
 Tang, F., & Moore, R. L. 1982, Sol. Phys., 77, 263 [NASA ADS] [CrossRef] [Google Scholar]
 Tsytovich, V. N. 1967, Sov. Phys. Uspekhi, 9, 805 [NASA ADS] [CrossRef] [Google Scholar]
 Umeda, T. 2010, J. Geophys. Res. (Space Phys.), 115, 1204 [NASA ADS] [CrossRef] [Google Scholar]
 Wild, J. P. 1950, Austr. J. Sci. Res. A Phys. Sci., 3, 541 [Google Scholar]
 Wild, J. P., & McCready, L. L. 1950, Austr. J. Sci. Res. A Phys. Sci., 3, 387 [Google Scholar]
 Yoon, P. H. 2000, Phys. Plasmas, 7, 4858 [NASA ADS] [CrossRef] [Google Scholar]
 Yoon, P. H. 2005, Phys. Plasmas, 12, 052313 [NASA ADS] [CrossRef] [Google Scholar]
 Yoon, P. H. 2006, Phys. Plasmas, 13, 022302 [NASA ADS] [CrossRef] [Google Scholar]
 Zharkova, V. V., & Siversky, T. V. 2011, ApJ, 733, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Zheleznyakov, V. V., & Zaitsev, V. V. 1970, Sov. Astron., 14, 250 [NASA ADS] [Google Scholar]
 Ziebell, L. F., Gaelzer, R., Pavan, J., & Yoon, P. H. 2008, Plasma Phys. Controlled Fusion, 50, 085011 [NASA ADS] [CrossRef] [Google Scholar]
 Ziebell, L. F., Yoon, P. H., Gaelzer, R., & Pavan, J. 2014a, ApJ, 795, L32 [NASA ADS] [CrossRef] [Google Scholar]
 Ziebell, L. F., Yoon, P. H., Gaelzer, R., & Pavan, J. 2014b, Phys. Plasmas, 21, 012306 [NASA ADS] [CrossRef] [Google Scholar]
 Ziebell, L. F., Yoon, P. H., Simões, F. J. R., Gaelzer, R., & Pavan, J. 2014c, Phys. Plasmas, 21, 010701 [NASA ADS] [CrossRef] [Google Scholar]
 Ziebell, L. F., Yoon, P. H., Petruzzellis, L. T., Gaelzer, R., & Pavan, J. 2015, ApJ, 806, 237 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: FokkerPlanck equation coefficients
Here we specifically follow the notations and conventions that can be found in Gaffey (1976), which can be summarized as follows. Starting from the Landau form of the FokkerPlanck collision operator defined in Eq. (3) and absorbing the ambient density into the definition for particle distribution function, , we have
(A.1)Integrating by parts we obtain (A.2)where (A.3)and x_{ab} ≡ v_{a}/v_{Tb}, with v_{Tb} denoting the thermal velocity for particle species b, v_{Tb} = (2T_{b}/m_{b})^{1/2}. Here, n_{0b} is the number density of particles of species b.
At this point we introduce the simplifying argument that in the beamplasma instability, the most significant evolution occurs in the tail of the particle distribution function, involving particle densities that are much lower than the background density. If so, then it is justified to assume that the most significant collisional effect will be due to collisions between tail particles with particles of the background distribution. Assuming that the background distribution is a Maxwellian distribution, the function G in Eq. (3) may be written as (Gaffey 1976) (A.4)where Φ(x_{ab}) is the error function, and Φ′(x_{ab}) its derivative,
(A.5)Making use of the auxiliary function Ψ, (A.6)and performing some simple algebra, one obtains the following form for the collisional term (Gaffey 1976), (A.7)Henceforth, we concentrate on the collision term affecting the electron distribution, so that f_{a} = f_{e}. Let us write the electron velocity simply as v, use the nondimensional variables τ = ω_{pe}t and u = v/v_{te}, and write the term describing binary collisions between electrons and particles of species b as (A.8)where we have considered n_{e} ≃ n_{i}, introduced Z_{i} and Z_{e} as the ion and electron charge number, respectively, with Z_{i} = Z_{e} = 1 and used g, as defined in Eq. (6).
In the present application, we consider a 2D system, which serves as a good approximation for the fully 3D case with azimuthal symmetry. In the 2D case, we make use of , (∂^{2}/∂u_{x}∂u_{z})u = (∂^{2}/∂u_{z}∂u_{x})u = −u_{x}u_{z}/u^{3}, and write Eq. (7) in the following form, (A.9)where
(A.10)In addition, it is easy to verify that the equation for the electron distribution function, obtained from (6), can be written for the case 2D in the following form: (A.11)where, as a further and usual approximation, we have neglected the effect of S waves. Making use of , the coefficients and the can be written as (A.12)
All Figures
Fig. 1
Time evolution of electron velocity distribution function is shown in 2D velocity space plot at different time steps. The lefthand panels depict solutions with only collisions (as in the thick target approximation), the middle panels correspond to the selfconsistent solution without collisional effects (as in the type III beam case), and the right panels show solutions that contain both collisional dynamics and selfconsistent Langmuir wave dynamics. 

In the text 
Fig. 2
Time evolution of Langmuir wave spectrum as a 2D wave number space plot at different time steps. The lefthand panels depict solutions with only selfconsistent collective effects without the collisional terms in the particle equation, and the right ones show solutions that contain both collisional dynamics and selfconsistent wave dynamics. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.