Skin Effects in a Dusty Plasma

A multifluid MHD model is applied to study the magnetic field dynamics in a dusty plasma. The motion of plasma electrons and ions is treated against the background of arbitrarily charged, immobile dust grains. When the dust density gradient is nonzero and when the inertia of the ions and electrons and the dissipation from their collisions with dust grains are neglected, we are dealing with a nonlinear convective penetration of the magnetic field into the plasma. When the dust density is uniform, the magnetic field dynamics is described by the nonlinear diffusion equations. The limiting cases of diffusion equations are analyzed for different parameter values of the problem (i.e., different rates of the collisions of ions and electrons with the dust grains and different ratios between the concentrations of the plasma components), and some of their solutions (including self-similar ones) are found. The results obtained can also be useful for research in solid-state physics, in which case the electrons and holes in a semiconductor may be analogues of plasma electrons and ions and the role of dust grains may be played by the crystal lattice and impurity atoms.     

Influence of the electrons reflected from the collector on the parameters of a high-current relativistic electron beam

In plasma microwave oscillators, electrons fall onto the surface of a graphite collector, which leads to the generation of secondary electrons. The influence of the electrons reflected from the collector on the parameters of a high-current relativistic electron beam propagating in a strong longitudinal magnetic field was studied experimentally and by numerical simulations. It is shown that the penetration of the reflected electrons into the drift space can lead to a substantial increase in the depth of the potential well in the drift space, a decrease in the velocity of the beam electrons, and a broadening of the electron energy distribution function.     

Kinetic theory of the positive column of a low-pressure discharge in a transverse magnetic field

The influence of a transverse magnetic field on the characteristics of the positive column of a planar low-pressure discharge is studied theoretically. The motion of magnetized electrons is described in the framework of a continuous-medium model, while the ion motion in the ambipolar electric field is described by means of a kinetic equation. Using mathematical transformations, the problem is reduced to a secondorder ordinary differential equation, from which the spatial distribution of the potential is found in an analytic form. The spatial distributions of the plasma density, mean plasma velocity, and electric potential are calculated, the ion velocity distribution function at the plasma boundary is found, and the electron energy as a function of the magnetic field is determined. It is shown that, as the magnetic field rises, the electron energy increases, the distributions of the plasma density and mean plasma velocity become asymmetric, the maximum of the plasma density is displaced in the direction of the Ampère force, and the ion flux in this direction becomes substantially larger than the counter-directed ion flux.     

Effect of the external magnetic field on the critical current for the onset of a virtual cathode in an electron beam

The critical current at which an unsteady oscillating virtual cathode forms in an electron beam is studied as a function of the external magnetic field guiding the beam electrons. It is shown that the critical beam current decreases with external magnetic field and that there is an optimum magnetic induction at which the critical current for the onset of an oscillating virtual cathode in the beam is minimum. For a strong guiding magnetic field, the critical beam current is described by relationships derived under the assumption that the motion of the beam electrons is one-dimensional. Such behavior is explained by the characteristic features of the dynamics of the beam electrons in longitudinal and radial directions in the interaction space at different inductions of the external magnetic field.     

Distribution of the absorbed ECRH power among the electrons in a plasma heated by an extraordinary wave

A study is made of the one-dimensional linear problem of the absorption of the energy of an extraordinary wave propagating along a nonuniform magnetic field by a plasma in the ECR region. The plasma electrons are assumed to be nonrelativistic and are described by a collisionless kinetic equation. The distribution of the absorbed power among the electrons and the distribution of the self-consistent field over the confinement system are obtained. The conditions under which the ECRH power is distributed uniformly among the bulk electrons are determined. The limits of applicability of the locally nonuniform magnetic field approximation are established. The solutions derived are compared with the solution to an analogous problem with the collisional absorption mechanism.     

Two-dimensional modeling of a transverse collisionless shock wave

Transverse collisionless shock waves in a plasma in which the initial β value is equal to zero for electrons and is small but nonzero for ions are studied in the two-dimensional approximation with allowance for anomalous resistivity. A hybrid model is applied such that the ions are treated in the kinetic approximation and the electrons are described in the hydrodynamic approximation. A collisionless shock wave is generated using a piston with a small two-dimensional perturbation. The ion distribution downstream of the shock front and the effect of electron and ion heating are analyzed. It is shown that, for Alfvén-Mach numbers M _A>2, ion heating is attributed primarily to the ions that have experienced a reflection from the shock front and whose velocities downstream of the front are very high. This conclusion agrees with the results of one-dimensional calculations. Solving the problem as formulated shows that two-dimensional effects are insignificant in the range of low Alfvén-Mach numbers (M _A≤5): the direction of the magnetic field is always close to its initial direction, the ions acquire low velocities along the magnetic field, and the quantitative parameters of the plasma downstream of the shock front are close to those obtained from the one-dimensional model. In the range of higher Alfvén-Mach numbers, two-dimensional effects are more pronounced and the ion distribution function is less anisotropic.     

Nonlinear oscillations of a semiconductor plasma with a nonrelativistic electron beam

Nonlinear oscillations of a semiconductor plasma with a low-density electron beam in the absence of an external magnetic field are studied in the hydrodynamic approximation. The beam is assumed to be nonrelativistic and monoenergetic. Cases are studied in which the Langmuir frequency of the electron oscillations in a semiconductor is much higher or much lower than the electron momentum relaxation rate. The self-similar solution obtained for the first case describes the damping of the nonlinear oscillations of the wave potential. Numerical analysis of the second case shows that the electric field distribution in the beam may correspond to that in a shock wave.     

Splitting of the modes of a low-frequency magnetosonic wave in a polydisperse dusty plasma

The properties of magnetosonic waves that propagate perpendicularly to the external magnetic field in a polydisperse dusty plasma and the frequencies of which are about the dust cyclotron frequency are analyzed. A dispersion relation containing integrals of functions of the dust grain radius is derived and investigated as a function of the parameters characterizing the polydisperse properties of dust. It is found that, in a polydisperse dusty plasma, the low-frequency magnetosonic mode splits into two branches. The first, lower frequency branch has a cutoff, while the higher frequency branch has a resonance. Between the two branches, there is a forbidden frequency range within which electromagnetic waves cannot propagate perpendicular to the magnetic field. The width of the forbidden frequency range is determined as a function of the slope of the distribution function of dust grains over radii and the interval within which the dust grain radii lie.     

Generation of auroral kilometric radiation by a finite-size source in a dipole magnetic field. Part II

Propagation and amplification of extraordinary electromagnetic waves in a dipole magnetic field in a narrow 3D plasma cavity in which a weakly relativistic electron beam propagates along the magnetic field in the direction of the gradient of the magnetic field strength is investigated. The domain of wave vectors at the starting point for which the wave amplification factors at the output of the density cavity reach their maximum values is found, and the amplification factor as a function of the wave frequency is determined. It is shown that the longitudinal velocity of fast electrons, which enables generation of waves in a broader frequency range, plays an important role in the formation of the spectrum of the auroral kilometric radiation (AKR). In this case, waves with the largest amplification factors at the output of the cavity have frequencies exceeding the cutoff frequency of the background plasma at the wave generation altitude. The global inhomogeneity of the magnetic field and plasma density, which governs the residence time of the waves in the amplification region, plays a key role in the formation of the AKR spectrum. It is shown that this time is the main factor determining the energy of the waves emerging from the source.     

Quasilinear modification of the spectra of cyclotron emission from a toroidal plasma near the ECRH frequency

A study is made of the modification of the spectra of electron cyclotron emission from an ECR heated plasma in a toroidal magnetic confinement system into which the heating radiation is launched from the low-field side. It is shown that, at frequencies close to the heating frequency, cyclotron emission can become more intense because of the deformation of the distribution function of the resonant electrons. This effect can be used to diagnose the slightly pronounced quasilinear perturbations of the electron distribution in the thermal energy range, which are typical of experiments on ECR plasma heating. Results of a qualitative analysis carried out for model electron distribution functions are presented, and examples of three-dimensional numerical simulations of a circular tokamak are described.     

Nonlocal equation for the symmetric part of the electron distribution function in an inhomogeneous plasma

A nonlocal kinetic equation is derived for the symmetric part of the distribution function of suprathermal electrons. It is shown that Albritton equations are merely local approximations to the total kinetic equation. Even in the simplest situation, the local approximations of the nonlocal effects are impossible to construct because of the interdependence of the variables. A self-similar solution to the equations under study is proposed.     

Acceleration of ions in a beam-plasma discharge in a low magnetic field: Interrelation between the electron and ion energy distributions

Previous experiments revealed the effect of stable acceleration of ions in a plasma-beam discharge in a low magnetic field to energies one order of magnitude higher than the electron thermal energy. To verify the previously proposed mechanisms for this effect, the velocity distribution function of the electrons arriving at the collector and the energy distribution of the ions escaping from the discharge transversely to the axis were measured. It is found that ion acceleration is accompanied by significant electron heating near the discharge axis. The time behavior and longitudinal profile of the intensity of the excited high-frequency oscillations in the frequency range ω ～ ω _pe were studied. The accumulation of regular oscillations in the beam-injection region and their stochastization during the propagation along the system axis were observed. The experimental results correlate qualitatively with the data of previous numerical simulations.     

Diffuse and constricted modes of a dc discharge in neon: Simulation of the hysteresis transition

Results are presented from theoretical studies of high-pressure (～100 Torr) dc discharges in neon. The diffuse and constricted discharge modes are studied using a model including the equation of balance for charged and excited particles, heat conduction equations for the neutral gas and plasma electrons, and Poisson’s equation for the radial electric field at a fixed total discharge current. A specific feature of the constricted mode in the investigated range of low fields and high degrees of ionization is that the excitation and ionization rates in the center of the discharge tube and at the periphery differ by several orders of magnitude. This implies that, in the constricted mode, the region where the electron energy distribution function is Maxwellian due to electron-electron collisions may adjoin the region (beyond the constriction zone) where the high-energy part of the distribution function is depleted. The hysteresis transition between the diffuse and constricted modes is analyzed. A transition from the constricted to the diffuse mode can be regarded as a manifestation of the nonlocal character of the formation of the electron distribution function, specifically, the diffusion of high-energy electrons capable of producing gas ionization from the central (constricted) region toward the periphery. The nonlocal formation of the distribution function is described by a nonlocal kinetic equation accounting for electron-electron collisions and electron transport along the radius of the discharge tube. Since only high-energy electrons produce gas ionization, the effect of the nonlocal formation of the electron distribution function is taken into account by introducing the effective temperature of the high-energy part of the distribution function and solving the equation for the radial profile of the high-energy part of the distribution function. This approach allows one to approximately take into account the nonlocal character of the electron distribution function without substantial expenditure of computer resources. The nonlocal model makes it possible to numerically simulate the hysteresis transition between the diffuse and constricted modes, which is impossible in the local approximation.     

Possible mechanism for radial ion acceleration in a beam-plasma discharge

A mechanism is proposed that can lead to radial ion acceleration in a plasma discharge excited by an electron beam in a relatively weak longitudinal magnetic field. The mechanism operates as follows. The beam generates an azimuthally asymmetric slow potential wave, which traps electrons. Trapped magnetized electrons drift radially with a fairly high velocity under the combined action of the azimuthal wave field (which is constant for them) and a relatively weak external longitudinal magnetic field. The radial electron flux generates a radial charge-separation electric field, which accelerates unmagnetized plasma ions in the radial direction. The ion flux densities and energies achievable in experiments with kiloelectronvolt electron beams in magnetic fields of up to 100 G are estimated.     

Oscillation spectrum of an electron gas with a small density fraction of ions

The problem is solved of the stability of a nonneutral plasma that completely fills a waveguide and consists of magnetized cold electrons and a small density fraction of ions produced by ionization of the atoms of the background gas. The ions are described by an anisotropic distribution function that takes into account the characteristic features of their production in crossed electric and magnetic fields. By solving a set of Vlasov-Poisson equations analytically, a dispersion equation is obtained that is valid over the entire range of allowable electric and magnetic field strengths. The solutions to the dispersion equation for the m = +1 main azimuthal mode are found numerically. The plasma oscillation spectrum consists of the families of Trivelpiece-Gould modes at frequencies equal to the frequencies of oblique Langmuir oscillations Doppler shifted by the electron rotation and also of the families of “modified” ion cyclotron (MIC) modes at frequencies close to the harmonics of the MIC frequency (the frequencies of radial ion oscillations in crossed fields). It is shown that, over a wide range of electric and magnetic field strengths, Trivelpiece-Gould modes have low frequencies and interact with MIC modes. Trivelpiece-Gould modes at frequencies close to the harmonics of the MIC frequency with nonnegative numbers are unstable. The lowest radial Trivelpiece-Gould mode at a frequency close to the zeroth harmonic of the MIC frequency has the fastest growth rate. MIC modes are unstable over a wide range of electric and magnetic field strengths and grow at far slower rates. For a low ion density, a simplified dispersion equation is derived perturbatively that accounts for the nonlocal ion contribution, but, at the same time, has the form of a local dispersion equation for a plasma with a transverse current and anisotropic ions. The solutions to the simplified dispersion equation are obtained analytically. The growth rates of the Trivelpiece-Gould modes and the behavior of the MIC modes agree with those obtained by numerical simulation.     

Two-dimensional electrodynamic structure of the normal glow discharge in an axial magnetic field

Results are presented from numerical simulations of an axisymmetric normal glow discharge in molecular hydrogen and molecular nitrogen in an axial magnetic field. The charged particle densities and averaged azimuthal rotation velocities of electrons and ions are studied as functions of the gas pressure in the range of 1–5 Torr, electric field strength in the range of 100–600 V/cm, and magnetic field in the range of 0.01–0.3 T. It is found that the axial magnetic field does not disturb the normal current density law.     

Formation of a cylindrical ion beam in the field of an axisymmetric system of ring currents at a low Hall current

A study is made of the structure of an accelerating layer with a closed Hall current and the geometry of an ion beam in an external magnetic field created by an arbitrary axisymmetric system of ring currents under conditions such that the Hall current can be ignored. It is shown that the ion trajectories are perpendicular to the magnetron cutoff surface for electrons and that the cathode plasma boundary coincides with a magnetic field line. A magnetic field configuration is found in which the cutoff surface is a plane surface perpendicular to the axis of the system. It is shown that, for a small ratio of the gyroradius of the electrons (in terms of the maximum energy acquired by them in the layer) to the characteristic size of the structure, such a configuration provides sufficient means to ensure the formation of slightly converging ion beams or those that are essentially parallel to the system axis.     

Fundamental parameters of a relativistic runaway electron avalanche in air

A relativistic runaway electron avalanche in air is simulated numerically by the Monte Carlo method with allowance for a large number of elementary processes involving electrons, positrons, and photons. The characteristic time scale of the avalanche amplification is calculated as a function of the overvoltage δ relative to the minimum value of the drag force between the electrons and the atomic particles of the medium. The dynamics of the formation of the electron energy distribution is investigated. The steady-state mean electron energy depends weakly on δ. Over a wide range of δ values, there exists a universal electron energy distribution, which is essentially independent of δ. The angular distributions of electrons integrated over energies, as well as the angular distributions for different energy groups, are calculated. Analytic approximations for the energy and angular distributions are obtained.     

Stochastic heating of electrons by a large-amplitude extraordinary wave in plasma

Stochastic heating of plasma electrons by a large-amplitude electromagnetic wave propagating across a strong external magnetic field is studied theoretically and numerically. An analytic estimate of the threshold wave amplitude at which heating begins is obtained. The dependence of the average electron energy on the magnetic field and plasma density is investigated using particle-in-cell simulations.     

Nonlinear Right-Hand Polarized Wave in Plasma in the Electron Cyclotron Resonance Region

The propagation of a nonlinear right-hand polarized wave along an external magnetic field in subcritical plasma in the electron cyclotron resonance region is studied using numerical simulations. It is shown that a small-amplitude plasma wave excited in low-density plasma is unstable against modulation instability with a modulation period equal to the wavelength of the excited wave. The modulation amplitude in this case increases with decreasing detuning from the resonance frequency. The simulations have shown that, for large-amplitude waves of the laser frequency range propagating in plasma in a superstrong magnetic field, the maximum amplitude of the excited longitudinal electric field increases with the increasing external magnetic field and can reach 30% of the initial amplitude of the electric field in the laser wave. In this case, the energy of plasma electrons begins to substantially increase already at magnetic fields significantly lower than the resonance value. The laser energy transferred to plasma electrons in a strong external magnetic field is found to increase severalfold compared to that in isotropic plasma. It is shown that this mechanism of laser radiation absorption depends only slightly on the electron temperature.