Two new, previously unknown types of dissipative streaming instabilities (DSI) are substantiated. They follow from new approach, which allows solving in general form the classical problem of an initial perturbation development for streaming instabilities (SI). SI is caused by relative motion of the streams of plasma components. With an increase in level of dissipation SI transforms into a DSI. The transformation occurs because dissipation serves as a channel for energy removal for the growth of the negative energy wave of the stream. Until recently, only one type of DSI was known. Its maximal growth rate depends on the beam density nb and the collision frequency ν in the plasma as ∼nb/ν. All types of conventional beam-plasma instabilities (Cherenkov, cyclotron, etc.) transform into it. The solution of the problem of the initial perturbation development in systems with weak beam-plasma coupling leads to a new type of DSI. With an increase in the level of dissipation, the instability in these systems transforms to the new DSI. Its maximal growth rate is ∼nb/ν. The second new DSI develops in beam-plasma waveguide with over-limiting current of e-beam. Its growth rate ∼nb/ν. In addition, the solutions of abovementioned problem provide much information about SI and DSI, significant part of which is unavailable by other methods.
- beam-plasma instability
- dissipative instability
- development of initial perturbation
- growth rate
- absolute/convective instability
Streaming instabilities (SI) occupy a prominent place among other plasma instabilities. They are caused by a motion of some plasma components relative to others. An example is the well-known beam–plasma instability . With this instability, the directed motion of a group of fast electrons passing through the background plasma excites potential oscillations with a large growth rate near the plasma frequency. Particular attention to this instability is mainly due to the idea of creating sources of powerful electromagnetic radiation on its basis. At present, these sources have many advantages over the known vacuum sources [2, 3]. One more example (we mention these two only) is the Buneman instability , in which plasma electrons move relative to ions.
In the overwhelming majority of investigations beam–plasma interaction is considered without any noticeable dissipation. It, actually, was assumed that the dissipation is small and cannot have any noticeable effect on the physical processes. In this case, the development of instability leads to an increase in the amplitude of electromagnetic oscillations in the plasma, as well as their energy at the expense of beam kinetic energy. In the absence of dissipation, the level of excited oscillations may be quite high, and their energy can even be comparable to the initial energy of the beam .
However, generally speaking, dissipation in the system (collisions between plasma particles, heating of metal surfaces due to their complex impedance, etc.) can play an essential role in plasma–beam interaction. It can become not only a decisive factor in limiting the spatial and temporal growth, determining the field amplitude and the mode structure and limits the growth rates. In addition to these properties, which are common to all systems, it is necessary to pay special attention to the unique role of dissipation in systems with a stream of charge particles: dissipation of high level does not suppress the SI completely. Strong dissipation transforms each SI to instability of other type – to dissipative streaming instability (DSI) . This type of instabilities is due to the presence of the negative energy wave (NEW) in a stream of charge particles [6, 7]. In fact, dissipation serves as a channel for energy removal for excitation of this wave. This leads to instabilities of a new physical nature, to DSI. Dissipation is the cause of this instability.
The physical nature of SI is not as simple as it might seem at first glance. It takes a lot of effort1 to understand it clearly. This is all the more so, if we are dealing with the transformation of SI into a DSI. The transformation (in general, the transformation of one type of instability into another) makes the behavior of SI in a system with dissipation especially interesting. In addition, there are other reasons that significantly increase interest in the study of problems associated with dissipation and the DSI caused by it. Some of them are as follows.
Modern high-frequency microwave electronics, both plasma and vacuum, have two basic trends of development: an increase in the frequency and power of the output radiation . With increasing frequency, the thickness of the skin layer on the resonators’ walls decreases. This, in turn, leads to an increase in active energy losses. Actual dissipation in the system increases.
The second trend – an increase in the power of output radiation – leads to the need to increase the beam current. The role of space charge phenomena increases also, as well as the role of the NEW. In these circumstances it becomes important to take into account all factors that also lead to the buildup of the same wave i.e. to dissipation. In a sense, dissipation becomes associated with the space charge phenomena. In addition with an increase in the beam current, the return current increases also. With account the decrease in the skin layer and the finite conductivity of metallic surfaces, this leads to an increase in the level of dissipation in the system. All this indicates that dissipation, along with the space charge of the beam plays an important role in microwave electronics. A detailed understanding of the role of all these phenomena is vital for many problems aimed at achieving high-intensity beams and their applications.
Until recently, only one DSI was known in beam–plasma interaction theory . Its maximal growth rate depends on collision frequency in plasma and on the beam density as . All types of the beam–plasma instabilities (Cherenkov, cyclotron, etc.), with an increase in the level of dissipation, transform into it. This only known DSI has a number of specific features in comparison with other (no-dissipative) instabilities: relatively low level of excited oscillations, relatively small growth rate, etc. Many investigations have been devoted to its study. It was assumed that various phenomena in space plasma and in plasma of controlled fusion can be explained on the basis of this instability.
However, recent studies have shown that there are other DSI also [9, 10, 11]. The interaction of the stream with the background plasma critically depends on some basic parameters of the system and/or on its geometry. Their changes lead to new physics of the beam-plasma interaction and to previously unknown types of DSI. The parameters are: the level of correlation between the beam and the plasma fields and the value of the beam current.
Available methods of instability investigation do not allow getting complete information on the process of transformation of given instability into another type. Is known the most complete information on instability can be obtained by solving the problem of the evolution of fields in space and time during the development of an initial perturbation. This problem is classical in theory of instabilities . Its results can clear up how the fields of given instability transform to the fields of another one along with many other accompanying details. The character of the space–time evolution of an initial perturbation is an important issue in many branches of physics. However, the results of this problem are hardly achievable. Ultimately its mathematical solution reduces to calculation of the integral with complete dispersion relation (DR) in the denominator of the integrand. For the result the DR should be specified and solved before integration. This sharply reduces generality of results. And even in the special cases, it is not always possible to carry out the integration. In  an approach is presented that allowed overcome difficulties and obtain analytical expression for the fields’ space–time structure for all types of conventional beam-plasma instabilities. Results show that with increase in level of dissipation all types of beam-plasma instabilities transform to the only known type of DSI.
This review shows that the number of DSI is not limited by the above-mentioned DSI. Two new types of DSI are substantiated. They follow from solution of the same classical problem of initial perturbation development. One of the DSI manifests itself in the results of solving the problem in systems with weak beam-plasma coupling. Weak interaction realizes if the beam and the plasma are spatially separated by a considerable distance. Under weak coupling the beam actually is left to its own and its proper oscillation come into play. Moreover, among them is the NEW. Its interaction with plasma causes instability, the growth rate of which reaches maximum at resonance of the plasma wave with the NEW. This resonance of wave–wave type was called “Collective Cherenkov effect” . An increase in the level of dissipation leads to a new DSI with the growth rate . Actually the new approach to solution of the classical problem has detected this new DSI.
The second new DSI appears in results of solving of the same problem in uniform cross-section beam-plasma waveguide with over-limiting e-beam. With an increase in the beam current the fields of its space charge affects more and more on the beam-plasma interaction. This manifests itself in two ways. Along with the increasing of the role of space charge oscillations, static fields of the beam space charge set an upper limit on the beam current that can pass through a given vacuum electro-dynamical system. The limit can be overcome by plasma filling. Plasma neutralizes the space charge of the beam. Plasma-filled waveguides can transmit e-beams with a current that is several times higher than the limiting current in vacuum waveguide. The fields of overlimiting e-beam space charge changes the character of its instability. The instability of over-limiting beams is not associated with any radiation mechanism [9, 14]. Its growth rate reaches maximum at the point of exact Cherenkov resonance and depends on the beam density as [9, 14, 15], With an increase in the level of dissipation, one more new type of DSI develops . Its growth rate depends on the parameters as .
In present review special attention is paid to systems, the geometry of which is similar to geometry of plasma microwave sources. These devices are a cylindrical waveguide with thin annular plasma and spatially separated thin annular e-beam. In this geometry the new types of DSI manifest themselves also .
In order to dispel all possible doubts about the correctness of the results, both new DSI are also substantiated by conventional analysis of the corresponding DR. To obtain a geometry-independent result for weak beam-plasma coupling we use perturbation theory based on smallness of the coupling parameter.
2. The only known DSI and transition to it
For the beginning we shortly present rezoning, from which follow: all types of beam-plasma instabilities (Cherenkov, cyclotron, beam instability in spatially periodical structure) transform to the only known DSI with the maximal growth rate (is the Langmuir frequency of the beam, is the collision frequency in plasma). The transition takes place with an increase in the level of dissipation. This help us to reveal a criterion for identification of DSI type.
In general, the dispersion relation (DR), describing a plasma system penetrating by an electron beam can be written as
where describes the plasma (without beam), but describes the beam contribution in the system dispersion
is the frequency, is the wave vector of perturbations, is Langmuir frequency of the e-beam, is the velocity of the beam electrons (directed along axis), is a polynomial with respect to and , . It is assumed that and , with the Cherenkov interaction, with the cyclotron interaction , (is the cyclotron frequency, is the harmonic number), and when e-beam interacts with the periodical structure, , is the length of spatial period.
The beam electrons interact with the proper oscillations of the system and the interaction leads to instability. Developing instability manifests itself most effectively at frequencies and wavelengths close to the proper frequencies of the system in the absence of the beam, and, at the same time, close to the beam natural frequencies. In fact, along with (1) following condition is met
All (conventional) beam-plasma instabilities, including DSI, follow from (1)–(3). With an increase in level of dissipation all types of no-dissipative instabilities (Cherenkov, cyclotron etc) transform into the well-known DSI. If one searches the solutions of DR (1) in the form (satisfies (1) and (3); this case called resonance instability) he arrives to the expression
All types of no-dissipative instabilities follow the first and the right-hand side term. In this case the dissipative (second) term in (4) is small. The DSI follows from the second term (when it is greater than the first term) and the right-hand side term. The relation between the respective growth rates and is
where the frequency of collisions in plasma is introduced (). The expression (5) presents relation between the growth rates of no-dissipative and dissipative instabilities. Below we use (5) and its analogs as a criterion for identification of DSI type.
3. Weak beam-plasma coupling. New type of DSI
3.1 Solution of the problem of initial perturbation development under weak beam-plasma coupling
The best way to study an instability in detail and its possible transformation to that of other type is the solving of the problem of initial perturbation development. The information obtained by other ways is insufficient and does not give any details. Here we present general (geometry independent) solution of the problem for weakly coupled beam-plasma systems.
Consider a system consisting of a mono-energetic rectilinear electron beam and cold plasma. To begin with, suppose the following: the plasma and the beam are weakly coupled (e.g. in a consequence of a sufficiently large distance between them). Let an initial perturbation arises at a point (the electron beam propagates in the direction ) at the instant and the instability begins developing. Our goal is to obtain the fields’ space–time distribution at an arbitrary instant and investigate in detail the instability behavior by analyzing obtained expression. In the process, we interest only the longitudinal structure of the fields, i.e., their dependence on the longitudinal coordinate and time . The transverse structure of the fields can be obtained by expanding in terms of the system’s eigenfunctions. In accordance with this, only two arguments are highlighted below: frequency and longitudinal component of the wave vector. Other arguments are irrelevant in the consideration below. To avoid overburdening the formulas, they are omitted.
In given case of weak beam-plasma coupling the instability is the result of the interaction of the beam negative energy wave (NEW) and the slowed down wave in the plasma. The interaction is of Collective Cherenkov type. We proceed from the theory of wave interaction in plasma . In terms of this theory the problem of the initial perturbation evolution under instability development in non-equilibrium plasma can be considered based on the set of partial differential equations for the amplitudes of the interacting waves: beam charge density wave and the slowed down electromagnetic wave in the plasma
where is the time, is the coordinate along the beam propagation direction, is a function determined by the initial conditions, is the directed velocity of the beam, is the group velocity of the resonant wave in plasma, , describes dissipation in plasma and is proportional to the frequency of collisions in it. The meaning of the denotation will be cleared up below. Note, the set (6) is meaningful irrespective of the problem of development of any instability. Generally, it describes resonant interactions between two waves in unstable medium. One only condition should be satisfied: the growth rate attains maximum under Collective Cherenkov Effect. If the maximum is attained under conventional Cherenkov Effect, as for conventional beam-plasma instabilities, the interaction should be described by other set of Equations .
The solution of the set (6) gives the dependence of the field’s amplitude on longitudinal coordinate and time under instability development. Applying the Laplace transformation with respect to time and the Fourier transformation with respect to the spatial coordinate , we obtain following expressions for the transform :
The field’s amplitude can be found by inverse transformation
where is the contour of integration with respect to . For given case it is a straight line that lies in the upper half plane of the complex plane and passes above all singularities of the integrand.
Thus, the problem has been reduced to the problem of integration in (8). It is somewhat simpler in comparison to the integral, which represents classical solution. Instead of full DR its analog stands. The analog is determined by interaction of the waves, participating in the instability development. This replacement simplifies integration. However, it remains difficult and many authors use roundabout methods carry out an expression for possible estimation of the fields behavior [17, 18]. Presented here method easily leads to the desired result i.e. to expression for space–time distribution of the fields. We merely transform the variables and to another pair and . The first integration (over ) may be carried out by the residue method and the integration contour must be closed in the lower half-plane. The first order pole is
The second integration (over ) cannot be carried out exactly, and we are forced to restrict ourselves to the approximate steepest descent method . This method gives result in the limit of relatively large
and is equal to
As a result we arrive to the following expression for the field’s space time structure under development of the instability in spatially separated beam-plasma system
3.2 Analysis of the instability development
The expression (12) looks very complicate. At first glance it is impossible to extract any information on the instability behavior from it. However, it turned out, the expression may be easily analyzing. Moreover, the results are obtained from scratch, i.e. they are not based on prior research. Substantial part of the information is unavailable by other way. In particular, the analysis clearly shows that with increase in level of dissipation the no-dissipative instability turns to a new type of DSI and provides detailed information on both instabilities.
The properties of the instability is determined mainly by the exponential factor
which provides many information: the temporal and the spatial growth rates, the spread of the unstable perturbations’ velocities, the nature of the instability (absolute or convective), the effect of dissipation on instability, etc.
First consider some general properties of the instability, which follow from (13).
It is easily seen that in the absence of dissipation unstable perturbations have velocities in the range from to . The wave packet moves in the beam propagation direction and, along with exponential growth of the fields, expands. Its length increases over time . The knowledge of the boundary velocities of unstable perturbations allows at once determining the nature of the instability (convective/absolute) based on the definition only, without reference to additional studies (we mean the Sturrock’s laws ). It is clearly seen that the instability is convective in the laboratory frame and other frames moving at velocities and . However, if the observer’s speed is within the range , then the same instability is absolute (see Figure 1).
Now we turn to determination of the meaning of the denotation in (6). For this we consider case and find the point of the field’s maximum from expression
Its root is i.e. the point of the field’s maximum moves at velocity
In the wave theory the velocity (15) is called convective velocity. It characterizes the spatial convection of the fastest growing perturbations. (15) shows that the peak of the wave packet disposes in its middle. The packet is symmetric with respect to its peak. Substitution of into the determines the field’s behavior in the maximum as , i.e. represents the maximal growth rate of the instability, which develops in absence of dissipation in systems with weak beam-plasma coupling. At the point the peak forms, because here the growth rate of perturbations is maximal.
The meaning of the parameter may also be determined from the DR (7) only, bypassing the results of integration (12). The general expression for the group velocity obtained from DR (7) has the limit (15) under . The same limit (note that ) leads to DR in form , i.e. the parameter is the imaginary part of complex frequency (the growth rate). In this case (absence of dissipation) the instability is due to interaction of the NEW with the plasma. To emphasize the important role of we add the respective indexes . Its dependence on specific parameters is found out below.
At a fixed point the field first grows up to the value that is reached at the instant where
Then the field decreases, and at the time the wave packet completely passes given point. The exponent is, in fact, the maximal spatial growth rate. At a given point, the field reaches its maximum at the moment when the peak has already passed it (see Figure 1). The reason is that perturbations moving at lower velocities reach the point for a longer time, and they have time to grow more. is the velocity of the most effectively amplified perturbations.
Thus, the solution of the problem of initial perturbation development along with other detailed information, gave results of conventional initial and boundary problems. This coincidence confirms correctness of developed approach (initial assumptions, mathematics, etc.). An additional advantage of the approach is in its geometry-independence. At first glance, the presented approach seems more complicated than traditional approaches, but this complexity is only apparent.
3.3 The influence of dissipation. New type of DSI
Dissipation significantly influences on the presented picture of the instability development and changes it. First of all, it suppresses slow perturbations. The wave packet shortens. The threshold velocity is determined from the condition and is equal
Only high-velocity perturbations (in the range ) grow. The change in the velocity of the trailing edge shortens the packet’s length and can affects the nature of instability (convective/absolute) if the frame’s velocity lies in the range . Also, dissipation limits the growth rates of perturbations with velocity . Substituting we have for the field , where
As expected, the growth rates fall down. Dissipation distorts the symmetry of the induced wave packet. In presence of dissipation the dynamics of the fields can be obtained from the same Eq. (14) accounting for dissipation. It has the form
The solution of (18) gives the point of the field maximum , where
This expression shows that with an increase in the level of dissipation, the peak shifts more and more to the front of the wave packet. This takes place along with the decreasing of the wave packet’s length. Substitution of into gives the field value in the peak and shows the respective growth rate as the function on the level of dissipation
As a criterion for the type of DSI this relation between the growth rates of DSI and the growth rate of SI sharply differs from that for the conventional case (5). Actually the expression (22) shows that with an increase in level of dissipation in weakly coupled beam-plasma systems the instability, caused by the beam’s NEW interaction with the plasma transforms to a new type of DSI. Its characteristic peculiarity is in new, previously unknown,
3.4 Substantiation of the new DSI by conventional analysis of the DR
From electro-dynamical point of view, a spatially separated beam-plasma system is nothing, but a multilayer structure. The traditional analytical consideration of such systems leads to a very cumbersome DR, which, in addition, is highly dependent on the geometry and greatly complicates with an increase in the number of layers. However, the importance of the problem and the need for its analytical investigation has led to development of specific methods. Here an approach is presented that allows avoiding abovementioned difficulties. Also, the approach has an important advantage: the procedure for obtaining the DR does not depend on specific shape/geometry. In other words, obtained results can be adapted to systems of any geometry. The approach considers the problem of weak beam-plasma interaction by perturbation theory. The small parameter, which underlies the theory, is the parameter of weak beam-plasma coupling. We briefly present here the basics of this approach accounting for dissipation .
Consider a system consisting of a mono-energetic rectilinear electron beam and cold plasma. To begin with, suppose the following: the plasma and the beam are weakly coupled (e.g. a consequence of a sufficiently large distance between them). We also assume their homogeneity in the cross section. The geometry of the system is not specified. It also is assumed that the beam current is less than the limiting current in the vacuum waveguide. Dissipation in the system is taken into account by the introduction the collisions in plasma. For simplicity, consideration is limited to the case of a strong external longitudinal (to the beam propagation direction) magnetic field, which prevents the transverse motion of the beam and plasma particles.
The small parameter underlying the perturbation theory is the parameter of weak coupling between the beam and the plasma (that is, the smallness of the integrals describing the overlap of beam and plasma fields). In the zero order approximation, the perturbation theory assumes independence of the beam and plasma. In the first-order approximation, the theory leads to the DR [11, 14].
and are the frequency and longitudinal component of the wave vector, are Langmuir frequencies for the plasma and the beam respectively, is the collision frequency in the plasma, is the velocity of the beam electrons, , is speed of light, is the coupling parameter, the point is the intersection point of the beam and the plasma dispersion curves, the values and play role of transverse wave numbers. Analytically, as well as and are expressed through the integrals of eigenfunctions of the zero order problem [11, 14]. The integral for represents overlap of the beam and plasma fields. It shows how far the plasma field penetrates the beam and vice versa. The specific expressions for , and are not essential for the subsequent presentation and are not presented here (see [11, 14]).
The expressions are the zero order DR for the plasma and the beam respectively. Their solutions are assumed to be known. The form of the DR (23) is comparatively simple. It shows the interaction of beam and plasma waves. Using (23) with small , it is easy to describe instabilities in given system. The main result of a decrease in the beam–plasma coupling is in the increase in role of the beam NEW. Its interaction with plasma leads to instability. The spectra of slow (−) and fast (+) beam waves follow from the roots of . If one searches them in form , , the roots become [11, 14].
where . Mathematically, the instability is due to corrections to the expression for the slow beam wave . Under collective Cherenkov resonance , the equation for is
In absence of dissipation the instability is due to NEW interaction with the plasma. Its growth rate is
We emphasize unusual dependence on the beam density as (for strong coupling this dependence is ∼ ). With ordinary Cherenkov resonance the system is stable. Under collective Cherenkov resonance dissipation manifested itself as an additional factor that enhances NEW growth and the instability gradually transforms to that of dissipative type. The Eq. (26) gives an expression for the growth rate as a function on level of dissipation
where is the function given in (21). The dependence of the growth rate on the level of dissipation in (28) coincides to that in (21). In limit (28) coincides to (27). In the opposite limit of strong dissipation (28) represents the growth rate of the new type of DSI (it also follows from (26) by neglecting the first term in brackets)
We arrive to the same new type of DSI presented in (22). The expression (28) shows a gradual transition of the growth rate of no-dissipative instability caused by NEW interaction with plasma into the growth rate of new type of DSI. It develops under weak coupling and differs from the conventional DSI (with an growth rate ∼ ). In  the same new DSI is substantiated in a finite external magnetic field.
4. Uniform cross section beam-plasma waveguide. One more new type of DSI
4.1 Evolution of the initial perturbation in plasma waveguide with over-limiting electron e-beam
One more new DSI arises under consideration of the problem of the initial perturbation development for the instability of over-limiting beam (OEB) in uniform cross-section plasma waveguide.
Consider a cylindrical waveguide, fully filled with cold plasma. A mono-energetic relativistic electron beam penetrates it. The external longitudinal magnetic field is assumed to be strong enough to freeze transversal motion of the beam and the plasma electrons. We also assume that the beam and plasma radii coincide with the waveguide’s radius and consider only the symmetrical
The character of the beam-plasma interaction changes depending on the beam current value. If the beam current is less than the limiting current in vacuum waveguide the instability is due to induced radiation of the system eigenwaves by the beam electrons. But, if the beam is over-limiting, its instability has the same nature as the instability in medium with negative dielectric constant [9, 14, 15]. We introduce a parameter , which represents the beam current value and the character of beam-plasma interaction. It corresponds (correct to the factor ) to the ratio of the beam current to the limiting current in vacuum waveguide  , i.e. (is the beam current). The values correspond to under-limiting beam currents , but the values correspond to over-limiting beam currents. This is possible under comparatively high values of the relativistic factor . Here we consider development of an initial perturbation in the system, when the beam current slightly exceeds the limiting vacuum value. In this case the instability is due to a-periodical modulation of the beam density in medium with negative dielectric constant. Its growth rate attaints maximum under exact Cherenkov resonance and is equal .
In order to show the variety of possible approaches to the solution of the problem of the initial perturbation development, in given case we solve it by other way. We turn to the set of origin equations, which describes e-beam instability in magnetized plasma waveguide
where and are the Bessel functions; and their roots in ascending order, , , . The quantities and should be expanded by analogy to , but – by analogy to . From here on we deal with the expansion coefficients and mention arguments and only.
The fields’ growth in the linear stage reveals itself most effectively on frequencies, closely approximating to roots of the DR and, simultaneously, to (resonant instability). The conditions (30) hold. In this connection it is reasonable to assume that originated perturbations form a wave packet of following type (e.g. for ):
where the carrier frequency and wave vector satisfy the conditions (30). We also assume that the amplitude of the wave train varies slowly in space and time as compared to and that is
Thus, the problem of the initial pulse behavior reduces to determination of the slowly varying amplitude (SVA) . The equation that satisfies can be derived from the set of origin Eqs. (32). The expansions (33) reduce it to a set of the equation for the amplitudes of expansions. In its turn the resulting set can be reduced to one equation for . We write it in form similar to the DR
where and are differential operators ; . The DR in form (30) follows from (36). To derive the equation for one should expanding (36) in power series near resonant values of frequency and wave vector by using the relations and with account of OEB existence condition. As a result we arrive to the following second-order partial differential equation for
where , and the expression for is obtained from the relation accounting the condition for OEB. It is important to emphasize that this denotation (as well as ) is introduced for reasons of simplicity of the resulting Eq. (36) only.
The solution of (37) is, actually, known. If one returns to the set (6) and transforms it (under ) to one equation for then the equation will completely coincide to (37). This means that we already have the solution of (36) and its analysis. It only remains to rewrite the solution (12) in new denotations and, where needed, re-interpret results. This shows that the instability in uniform cross-section beam-plasma waveguide develops in space and time in the same manner as the instability in weakly coupled beam-plasma system, and is its growth rate in limit , that is . However there is a very important quantitative difference. In present case the growth rate depends on the beam density as (for the case of weak beam-plasma coupling the dependence is (see (27)). The criterion for determining the type of DSI takes the form
Comparison of (38) with (22) indicates one more new type of DSI. It develops in uniform cross section beam-plasma waveguide under over-limiting beam current and high level of dissipation. Its growth rate depends on the beam density and collision frequency as .
4.2 Substantiation of the second new DSI by conventional method
where , , , , , is the group velocity of the resonant wave in the system without beam, ; is the resonant frequency of the plasma waveguide.
The solutions of (39) depend on the beam current value that is on the value of parameter α. If (under-limiting e-beams) one can obtain the growth rates of conventional instability under (first and right-hand side terms) and in limit i.e. DSI
If the beam current increases and become comparable or higher than the limiting vacuum current i.e. , the physical nature of the instability changes. It becomes due to a-periodical modulation of the beam density in medium with negative dielectric constant. The distinctive peculiarity of this instability is in following: its growth rate attains maximum under exact Cherenkov resonance and is equal to (31) [9, 11, 14, 15]. If, along with the beam current, dissipation also increases the instability turns to DSI of over-limiting beam with growth rate .
We emphasize new dependences on and on the beam density. This, actually, substantiates one more new type of DSI. It develops in uniform cross section beam-plasma waveguide if the beam current is higher than the limiting vacuum current.
5. The role of the new DSI in no-uniform-cross-section beam-plasma waveguide
5.1 Statement of the problem. Dispersion relation
In this section we pay special attention to systems, the geometry of which is similar to the geometry of plasma microwave sources and possible development of the new types of DSI in such systems. The simplest theoretical model of plasma microwave generators assumes relativistic e-beam propagating along axis of a plasma filled waveguide of radius
For theoretical study of the problem we use an approach , which gives result for arbitrary level of beam-plasma coupling. This condition is obligatory for obtaining comprehensive results. The DR, which follows from the approach, has a form, which clearly shows interaction of the beam and plasma waves. The approach proceeds from equation for polarization potential
Here and are perturbations of the longitudinal current densities in the beam and plasma. Functions describe transverse density profiles of the perturbations of the longitudinal currents in the beam and the plasma. For homogeneous beam/plasma for infinitesimal thin (is Dirac function), is the Laplace operator over transverse coordinates,
where are the Langmuir frequencies for plasma and beam respectively, is the effective collision frequency in plasma, , is the velocity of beam electrons.
The DR, which follows from the statement, is still very cumbersome (of integral type). To reduce the DR to a simple algebraic form one should make following expedient for theoretical model assumption: the plasma and the beam are not just thin but infinitesimal thin. In this case the DR becomes
(and are modified Bessel and Mac-Donald functions, is the azimuthal wave numbers).
An important property of
5.2 Growth rates
The DR (44) determines proper oscillations of transversally no uniform beam-plasma waveguide. The changes of the physical character of beam-plasma interaction must reveal themselves on its solutions. are the DR for waveguide with thin annular plasma and e-beam respectively. The spectra of fast (+) and slow (−) waves are
where . The parameter is familiar (see above). It determines the beam current value: (is the beam current, is the limiting current in vacuum waveguide). In the limit of under-limiting beams . In opposite limit of over-limiting beam and . If one looks for solutions of (44) in form , it becomes
where . The Eq. (48) presents sound way to study instabilities in given system. First of all, it is easily seen that in conditions of growing negative energy wave and collective Cherenkov resonance the role of dissipation increases. For under-limiting e-beams and in case of strong coupling the DR (44) leads to the well-known conventional beam instabilities of no-dissipative and dissipative type. The growth rates of these instabilities have well-known dependencies on beam density and on dissipation (). Both for these instabilities proper oscillations of the beam are neglected. Only for explanation of the physical meaning of the DSI the conception of NEW should be invoked. However, if (weak coupling) the growing of the NEW plays dominant role. In this case the growth rate of no-dissipative instability reaches its maximum under Collective Cherenkov resonance and is equal
As it should be, this is the instability discovered under consideration of the classical problem of the initial perturbation development in weakly coupled beam-plasma systems.
Of particular interest are limit of high, over-limiting currents of e-beam . In this case the DR (44) takes the form
For the analysis of (51) leads to following. Under single particle resonance we have either instability of negative mass type (under ) with the growth rate , or stability (under
The instability (52) under
However, the dependence on the beam density is completely different. This is the same DSI, which develops in uniform cross-section beam-plasma waveguide under over-limiting currents. Instabilities of the same type may be substantiated for finite thicknesses of the beam and plasma layers in waveguide. In this case one must use perturbation theory based on smallness of coupling coefficient.
As follows from this section, in the geometry of microwave plasma sources, the development of both new DSI is possible. Basic parameters of the both new DSI, and the conditions of their development should be taken into account upon design of the high power, high frequency plasma microwave devices.
Thus, based on very general initial assumptions, we have found out that the number of DSI in the beam-plasma interaction theory is not limited by the only previously known type. Two new, previously unknown types of DSI are presented. The new DSI reveal themselves in the analysis of solution of the problem of initial perturbation development. This problem is classical in the theory of instabilities. The first new DSI is the dissipative instability under weak beam-plasma coupling. In absence of dissipation the instability in these systems is caused by the interaction of the beam NEW with the plasma. With an increase in the level of dissipation this instability gradually transforms to the new type of DSI. Its maximal growth rate depends on the beam Langmuir frequency and the frequency of collisions in plasma as . This, more critical (as compared to conventional), inverse proportional dependence on is a result of superposition of two factors those lead to growth of the beam NEW: weak coupling and dissipation.
The second new type of DSI is dissipative instability of over-limiting e-beam in uniform cross section waveguide. With increase in the beam current, its space charge and inner degrees of freedom reveal themselves more efficiently. If the beam current becomes higher than the limiting current in vacuum waveguide then the instability mechanism changes. In uniform cross section beam-plasma waveguide the instability becomes due to a-periodical modulation of the beam density in medium with negative dielectric constant. In this case the increase in the level of dissipation leads to one more new type of DSI with the maximal growth rate .
The same types of DSI develop in systems having geometry, similar to microwave sources: cylindrical waveguide with thin annular beam and thin annular plasma. If the coupling between the beam and the plasma hollow cylinders is weak and the beam current is under-limiting the first type of DSI develops, but under over-limiting currents – the second. However, if the coupling of the beam and the plasma cylinders is strong, conventional type of DSI develops with well-known growth rate .
Both new DSI are confirmed by conventional analysis of the respective DR.
Some words about the approach used. It has many advantages. First of all, it is based on very general initial assumptions and gives results regardless on geometry and specific parameters. The same approach is used for solving the same problem for conventional beam-plasma instabilities of all types (Cherenkov type, cyclotron type etc) , for the Buneman instability  etc. Obtained expressions for the spatial–temporal distribution of growing fields clearly show that with increase in the level of dissipation in background plasma, all these SI transform into DSI of conventional type. In addition, the analysis of obtained expressions gives much more detailed information on SI than other methods give. Part of the information on SI is not available in any other ways. The coincidence of other information to the results of conventional analysis confirms the validity of the approach (initial assumptions, mathematics etc).
Also, the presented approach shows that the DR describing the SI of given type can serve not only for solving of the initial/boundary problems and obtaining the dispersion curves. This point of view is very simplified. The approach shows that much more additional information is available from the DR. It, in fact, provides results on the initial perturbation development.
Summarizing, one can state that the presented approach can serve as an independent and very effective method for studying of any SI. There is no need to solve the problem again. One should only substitute the parameters of given instability in general expression for the field’s space–time distribution. The usage of this approach instead of traditional initial/boundary problems gives complete picture of the instability development. At first glance, it might seem that this method of analyzing instabilities is more complicated. However, this complexity is only apparent. In addition, this complexity, if any, is overlapped by the completeness of the information received.
Aleksandrov A.F., Bogdankevich L.S., Rukhadze A.A. Principles of Plasma Electrodynamics Springer-Verlag: Berlin, Germany, 1984
Kuzelev M.V., Rukhadze A.A. Plasma Phys. Rep. 2000, v. 26, p 231
Kuzelev M.V., Loza O.T., Rukhadze A.A.; Strelkov, P.S., Shkvarunets A.G. Plasma Phys. Rep. 2001, v. 27, p. 669.
Buneman O. Phys. Rev. 1959, 115, 503.
Matsiborko, N.G.; Onischenko, I.N.; Shapiro, V.D.; Shevchenko, V.I. Plasma Phys. 1972, 14, 591.
Briggs, R.J. Advances in Plasma Phys; Symon, A.; Tomson, W., Eds.; Wiley: New York, NY, USA, 1971; v. 3, p. 132.
Nezlin, M.V. Physics of Intense Beams in Plasmas; IOP Publishing Bristol, 1999.
Kurilko V.I. Doklady Sov Akad, 1974, v208, p. 1059.
Rostomyan E.V. IEEE Trans. Plasma Sci. 2003, 31, 1278.
Rostomyan E.V. EPL 2007, 77, 45001.
Rostomyan E.V. J. Plasma Phys. 2012 2012, 78, 531.
Bers A. In “Basics of Plasma Physics” Eds Rosenbluth, M.N., Sagdeev, R.Z. North Holland: Amsterdam, The Netherlands, 1983; v. 1, p. 451
Rostomyan E.V. Phys Plasmas, 2000, v. 7,p 1595.
Kuzelev, M.V., Rukhadze, A.A. Plasma Free Electron Lasers; Frontieres: Paris, France, 1995.
Aizatski N.I. Sov. J. Plasma Phys. 1980, 6, 597.
Weiland J, Wilhelmsson H. Coherent Nonlinear Interaction of Wave in Plasma. Pergamon Press. Oxford. 1977.
Kondratenko A.N., Kuklin V.M. Basics of Plasma Electronics. Atomizdat. Moscow, 1988.
Kuzelev M.V. Plasma Phys Reports, 2006, v. 32, p. 572.
F. W. J. Olver, Asymptotics and Special Functions (Academic, New York, 1974;
Sturrock P.A. Phys. Rev 1958 v. 112 p. 1488
Rostomyan E.V. Phys Lett A. 2009. v. 373, p. 2581
Rostomyan E.V. Phys Plasmas, 2017, v. 24, 102102.
- The instability of low density e-beam in plasma is a vivid example demonstrating this sense. It is discovered in 1948, experimentally proven in early sixties; however its physical meaning became finally clear in the middle of seventies (see ).