RUDI GAELZER

Generation of electrostatic multiple harmonic Langmuir modes during beam–plasma interaction process has been observed in laboratory and spaceborne active experiments, as well as in computer simulation experiments. Despite earlier efforts, such a phenomenon has not been completely characterized both theoretically and in terms of numerical simulations. This paper is a first in a series of three papers in which analytic expressions for harmonic Langmuir mode dispersion relations are derived and compared against the numerical simulation result.

We study the dispersion relation for low frequency waves in the whistler mode propagating along
the ambient magnetic field, considering ions and electrons with product-bi-kappa (PBK) velocity
distributions and taking into account the presence of a population of dust particles. The results
obtained by numerical analysis of the dispersion relation show that the decrease in the j indexes in
the ion PBK distribution contributes to the increase in magnitude of the growth rates of the ion
firehose instability and the size of the region in wave number space where the instability occurs. It
is also shown that the decrease in the j indexes in the electron PBK distribution contribute to
decrease in the growth rates of instability, despite the fact that the instability occurs due to the
anisotropy in the ion distribution function. For most of the interval of j values which has been
investigated, the ability of the non-thermal ions to increase the instability overcomes the tendency
of decrease due to the non-thermal electron distribution, but for very small values of the kappa
indexes the deleterious effect of the non-thermal electrons tends to overcome the effect due to the
non-thermal ion distribution.

We investigate the dispersion relation for low frequency electromagnetic waves propagating parallel to the ambient magnetic field, considering that the velocity distributions of ions and electrons can be either bi-Maxwellian of product bi-kappa distributions. The effect of the anisotropy and non-thermal features associated to the product-bi-kappa distributions on the firehose instability are numerically investigated. The general conclusion to be drawn from the results obtained is that the increase in non-thermal features which is consequence of the decrease of the κ indexes in the ion distribution contributes to increase the instability in magnitude and wave number range, in comparison with bi-Maxwellian distributions with similar temperature anisotropy, and that the increase of non-thermal features in the electron distribution contributes to the quenching of the instability, which is nevertheless driven by the anisotropy in the ion distribution. Significant differences between results obtained either considering product-bi-kappa distributions or bi-kappa distributions are also reported.

This article demonstrates the generation of enhanced ion-acoustic waves by beam-plasma instability in a density cavity. The self-consistent equations of weak turbulence theory that include quasi-linear, decay, and scattering processes as well as convective and dispersive effects are numerically solved for a one-dimensional density cavity. It is shown that significant enhancements of ion-acoustic waves occur in the presence of counterstreaming electron beams and that the enhanced ion-acoustic waves are initially localized near the center of the density cavity at large wavelengths. Later in the evolution, the enhancement in the spectrum of ion-acoustic waves spreads out toward the edges of the cavity, with a shift to smaller wavelengths, while the enhancement near the center of the cavity tends to decrease in magnitude. The significance of the present findings is discussed.

The dispersion relation for parallel propagating waves in the ion-cyclotron branch is investigated
numerically by considering that the velocity distribution of the ion population is a function of type
product-bi-kappa. We investigate the effects of the non-thermal features and of the anisotropy associated with this type of distribution on the ion-cyclotron instability, as well as the influence of different forms of the electron distribution, by considering Maxwellian distributions, bi-kappa distributions, and product-bi-kappa distributions. The cases of ions described by either Maxwellian or bi-kappa distributions are also considered, for comparison. The results of the numerical analysis show that the increase in the non-thermal character associated with the anisotropic kappa distributions for ions contributes to enhance the instability as compared to that obtained in the Maxwellian case, in magnitude and in wave number range, with more significant enhancement for the case of ion product-bi-kappa distributions than for the case of ion bi-kappa distributions. It is also shown
that the ion-cyclotron instability is decreased if the electrons are described by product-bi-kappa distributions, while electrons described by bi-kappa distributions lead to growth rates which are very
similar to those obtained considering a Maxwellian distribution for the electron population.

Ion-acoustic enhancements are investigated within the context of turbulent beam-plasma interaction processes. The analysis assumes a pair of counterstreaming electron beams interacting with the background plasma. Two-dimensional velocity space and two-dimensional wave number space are assumed for the analysis, with physical parameters that characterize typical ionospheric conditions. The solutions of the electrostatic weak turbulence theory show that the ion-acoustic wave levels are significantly enhanced when the computation is initialized with a pair of counterstreaming beams in contrast to a single beam. We suggest that this finding is highly relevant for the observed ion-acoustic enhancements in the Earth's ionosphere that are known to be correlated with auroral activity.

In a pair of recent papers (Ziebell et al 2008 Phys. Plasmas 15 032303, 2008 Plasma Phys. Control. Fusion 50 085011) it was shown, within the context of weak turbulence theory, that the Langmuir turbulence generated by the bump-in-tail instability does not lead to Langmuir condensation (or accumulation of wave energy and momentum in the long-wavelength regime) in two dimensions. The present analysis finds that it is important to include the spontaneous scattering off Langmuir turbulence of the electrons, which is ignored in the customary literature when compared with a similar process involving ions, in order to recover the condensation of Langmuir waves in two dimensions.

Charged particle acceleration takes place ubiquitously in the Universe including the near-Earth heliospheric environment. Typical in situ spacecraft measurements made in the solar wind show that the charged particle velocity distribution contains energetic components with quasi scale-free power-law velocity dependence, f ∼ v − α , for high velocity range. In this Review a theory of quiet-time solar-wind electrons that contain a suprathermal component is discussed, in which these electrons are taken to be in dynamical equilibrium with Langmuir turbulence. This Review includes an overview of the Langmuir turbulence theory, as well as a discussion on asymptotic equilibrium solution of Langmuir turbulence/suprathermal electron system. Theoretical predictions of high-energy electron velocity power-law distribution index is then compared against the recent observations of the superhalo electron velocity distribution made by instruments onboard WIND and STEREO spacecraft. It is shown that the theoretical prediction of velocity power-law index is intermediate to the observed range.

A kinetic approach to the problem of wave propagation in dusty plasmas, which takes into account the variation of the charge of the dust particles due to inelastic collisions with electrons and ions, is utilized as a starting point for the development of a new formulation, which writes the components of the dielectric tensor in terms of a finite and an infinite series, containing all effects of harmonics and Larmor radius. The formulation is quite general and valid for the whole range of frequencies above the plasma frequency of the dust particles, which were assumed motionless. The formulation is employed to the study of electrostatic waves propagating along the direction of the ambient magnetic field, in the case for which ions and electrons are described by Maxwellian distributions. The results obtained in a numerical analysis corroborate previous analysis, about the important role played by the inelastic collisions between electrons and ions and the dust particles, particularly on the imaginary part of the dispersion relation. The numerical analysis also show that additional terms in the components of the dielectric tensor, which are entirely due these inelastic collisions, play a very minor role in the case of electrostatic waves, under the conditions considered in the present analysis.

Nonlinear interactions of tenuous electron beam, background, unmagnetized plasma, and self-consistently generated Langmuir and ion-sound waves are analyzed in the framework of plasma weak turbulence kinetic theory. Full numerical solutions of the complete weak turbulence equations are obtained for the first time, which show the familiar plateau formation in the electron beam distribution and concomitant quasi-saturation of primary Langmuir waves, followed by fully nonlinear processes which include three-wave decay and induced-scattering processes. A detailed analysis reveals that the scattering off ions is an important nonlinear process which leads to prominent backscattered and long-wavelength Langmuir wave components. However, it is found that the decay process is also important, and that the nonlinear development of weak Langmuir turbulence critically depends on the initial conditions. Special attention is paid to the electron-to-ion temperature ratio, Te/Ti, and the initial perturbation level. It is found that higher values of Te/Ti promote the generation of backscattered Langmuir wave component, and that a higher initial wave intensity suppresses the backscattered component while significantly enhancing the long-wavelength Langmuir wave component.

The problem of electron-beam propagation in inhomogeneous solar wind is intimately related to the solar type II and/or type III radio bursts. Many scientists have addressed this issue in the past by means of quasi-linear theory, but in order to fully characterize the nonlinear dynamics, one must employ weak-turbulence theory. Available numerical solutions of the weak-turbulence theory either rely on only one nonlinear process (either decay or scattering), or when both nonlinear terms are included, the inhomogeneity effect is generally ignored. The present paper reports the full solution of weak-turbulence theory that includes both decay and scattering processes, and also incorporating the effects of density gradient. It is found that the quasi-linear effect sufficiently accounts for the primary Langmuir waves, but to properly characterize the back-scattered Langmuir wave, which is important for eventual radiation generation, it is found that both nonlinear decay and scattering processes make comparable contributions. Such a finding may be important in the quantitative analysis of the plasma emission process with application to solar type II and/or type III radio bursts.

In the present article, the classic problem of nonlinear frequency shifts of electrostatic plasma eigenmodes in an unmagnetized plasma (i.e., the Langmuir and ion-acoustic waves) is revisited. In the standard literature, only the frequency shift of Langmuir waves by the finite-amplitude Langmuir waves themselves is usually treated. In the present approach, the discussion is generalized to include the ion-sound waves. The significance of the present article is that the analytical approach employed in the present discussion can be utilized to resolve certain apparently singular terms in the induced scattering coefficients of the wave kinetic equations. The detailed discussion of such a problem will be reported in a forthcoming article.

Velocity distribution functions (VDFs) that exhibit a power-law dependence on the high-energy tail have been the subject of intense research by the plasma physics community. Such functions, known as kappa or superthermal distributions, have been found to provide a better fitting to the VDFs measured by spacecraft in the solar wind. One of the problems that is being addressed on this new light is the temperature anisotropy of solar wind protons and electrons. In the literature, the general treatment for waves excited by (bi-)Maxwellian plasmas is well-established. However, for kappa distributions, the wave characteristics have been studied mostly for the limiting cases of purely parallel or perpendicular propagation, relative to the ambient magnetic field. Contributions to the general case of obliquely-propagating electromagnetic waves have been scarcely reported so far. The absence of a general treatment prevents a complete analysis of the wave-particle interaction in kappa plasmas, since some instabilities can operate simultaneously both in the parallel and oblique directions. In a recent work, Gaelzer and Ziebell [J. Geophys. Res. 119, 9334 (2014)] obtained expressions for the dielectric tensor and dispersion relations for the low-frequency, quasi-perpendicular dispersive Alfvén waves resulting from a kappa VDF. In the present work, the formalism introduced by Ref. 1 is generalized for the general case of electrostatic and/or electromagnetic waves propagating in a kappa plasma in any frequency range and for arbitrary angles. An isotropic distribution is considered, but the methods used here can be easily applied to more general anisotropic distributions, such as the bi-kappa or product-bi-kappa.

A kinetic formulation developed to analyze wave propagation in dusty plasmas, which takes into account the charge variation of the dust particles, is utilized to study the propagation and damping of Alfven waves propagating in oblique directions relative to the ambient magnetic field. A dusty plasma containing spherical and immobile dust grains in a homogeneous ambient magnetic field is considered. The charging process of the dust grains is assumed to be associated with the capture of electrons and ions by the dust particles during inelastic collisions between them and plasma particles. The dispersion relation and the damping rates of obliquely propagating Alfven waves are obtained assuming Maxwellian distributions for electrons and ions in equilibrium. For the numerical analysis of the dispersion relation we use the average values of the inelastic collision frequency as an approximation, instead of the momentum dependent expressions originally derived in the kinetic formulation, and study the modifications which the presence of the dust particles causes in both the propagation and the damping of the Alfven waves. In particular is discussed the competition between the different damping mechanisms, namely, the Landau damping and the damping associated with the dust charge variation, and it is shown that the inelastic collision frequency plays a pivotal role in the magnitude of the damping rates.

Yoon and Fang [Phys. Plasmas 15, 122312 (2008)] formulated a second-order nonlinear kinetic theory that describes the turbulence propagating in directions parallel/anti-parallel to the ambient magnetic field. Their theory also includes discrete-particle effects, or the effects due to spontaneously emitted thermal fluctuations. However, terms associated with the spontaneous fluctuations in particle and wave kinetic equations in their theory contain proper dimensionality only for an artificial one-dimensional situation. The present paper extends the analysis and re-derives the dimensionally correct kinetic equations for three-dimensional case. The new formalism properly describes the effects of spontaneous fluctuations emitted in three-dimensional space, while the collectively emitted turbulence propagates predominantly in directions parallel/anti-parallel to the ambient magnetic field. As a first step, the present investigation focuses on linear wave-particle interaction terms only. A subsequent paper will include the dimensionally correct nonlinear wave-particle interaction terms.

The prEservation of Onsager symmetry for the effective dielectric tensor is discussed for a homogeneous Plasma immersed in a inhomogeneous magnetic ?eld, using the unperturbed orbits correct up to order $k\_B$, which is the scalelength of the field inhomogeneity. General features of the calculation of the components of the tensor are discussed and detailed calculations are developed for the $zz$ component, which is shown to satisfy the conditions for Onsager symmetry, in agreement with prEvious results obtained using less prEcise exprEssions for the unperturbed orbits.

It is by now well known that electron beams play an important role in generating radio emissions such as type II and type III radio bursts, commonly observed by spacecraft in the interplanetary medium. Electron beams streaming back from Earth's bow shock into the solar wind have been proposed as a possible source for the electron plasma waves observed by spacecraft in the electron foreshock. Recent observations suggest that during the natural evolution of the foreshock plasma, multiple electron beams could be injected over a period of time, losing their individual identity to coalesce into a single beam. In this work, we use an electromagnetic particle-in-cell (PIC) code &\#8220;KEMPO 1D, adapted&\#8221; to simulate two electron beams that are injected into a plasma at different times. The first beam disturbs the background plasma and generates Langmuir waves by electron beam-plasma interaction. Subsequently, another beam is inserted into the system and interacts with the first one and with the driven Langmuir waves to produce electromagnetic radiation. The results of our simulation show that the first beam can produce electrostatic harmonics of the plasma frequency, while the second beam intensifies the emission at the harmonics that is produced by the first one. The behavior of the second beam is strongly determined by the preexisting Langmuir wave electric fields. The simulations also show, as a result of the interaction between both beams, a clear nonlinear frequency shift of the harmonic modes as well as an increase of electromagnetic and kinetic energies of the wave-particle system.

In this paper an electromagnetic particle code is used to investigate the spontaneous thermal emission. Specifically we perform particle-in-cell simulations employing a non-relativistic isotropic Maxwellian particle distribution to show that thermal fluctuations are related to the origin of spontaneous magnetic field fluctuation. These thermal fluctuations can become seed for further amplification mechanisms and thus be considered at the origin of the cosmological magnetic field, at microgauss levels. Our numerical results are in accordance with theoretical results presented in the literature.

The radiation emission mechanism responsible for both type-II and type-III solar radio bursts is commonly
accepted as plasma emission. Recently Ganse et al. suggested that type-II radio bursts may be enhanced when the electron foreshock geometry of a coronal mass ejection contains a double hump structure. They reasoned that the counter-streaming electron beams that exist between the double shocks may enhance the nonlinear coalescence interaction, thereby giving rise to more efficient generation of radiation. Ganse et al. employed a particle-in-cell simulation to study such a scenario. The present paper revisits the same problem with EM weak turbulence theory, and show that the fundamental (F) emission is not greatly affected by the presence of counter-streaming beams, but the harmonic (H) emission becomes somewhat more effective when the two beams are present. The present finding is thus complementary to the work by Ganse et al.
Key words: plasmas – radiation mechanisms: non-thermal – solar wind – Sun: radio radiation – turbulence – waves

The plasma emission, or electromagnetic (EM) radiation at the plasma frequency and/or its harmonic(s), is generally accepted as the radiation mechanism responsible for solar type II and III radio bursts. Identification and characterization of these solar radio burst phenomena were done in the 1950s. Despite many decades of theoretical research since then, a rigorous demonstration of the plasma emission process based upon first principles was not available until recently, when, in a recent Letter, Ziebell et al. reported the first complete numerical solution of EM weak turbulence equations; thus, quantitatively analyzing the plasma emission process starting from the initial electron beam and the associated beam-plasma (or Langmuir wave) instability, as well as the subsequent nonlinear conversion of electrostatic Langmuir turbulence into EM radiation. In the present paper, the same problem is revisited in order to elucidate the detailed physical mechanisms that could not be reported in the brief Letter format. Findings from the present paper may be useful for interpreting observations and full-particle numerical simulations.