On electrostatic modes in multiion and pairion collisional plasmas
J. Vranjes
Center for Plasma Astrophysics, Celestijnenlaan 200B, 3001 Leuven, Belgium, and Faculté des Sciences Appliquées, avenue F.D. Roosevelt 50, 1050 Bruxelles, Belgium.
D. Petrovic,
Research Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, 2610 Antwerp, Belgium
B. P. Pandey
Department of Physics, Macquarie University, Sydney, NSW 2109, Australia.
S. Poedts
Center for Plasma Astrophysics, and Leuven Mathematical Modeling and Computational Science Centre (LMCC) Celestijnenlaan 200B, 3001 Leuven, Belgium.
Abstract: The physics of plasmas containing positive and negative ions is discussed with special attention to the recently produced pairion plasma containing ions of equal mass and opposite charge. The effects of the density gradient in the direction perpendicular to the ambient magnetic field vector, observed in the experiment, are discussed. In addition, the possible presence of electrons is discussed in the context of plasma modes propagating at an angle with respect to the magnetic field vector. The electrons may either be added to the plasma or enter the plasma attached to negative ions and then become released from the ions as a result of collisions. It is shown that the electron plasma mode in a cold plasma may become a backward mode in the presence of a density gradient, and this behavior may be controlled either by the electron number density or the mode number in the perpendicular direction. In plasmas with hot electrons an instability may develop, driven by the combination of electron collisions and the density gradient, and in the regime when the ions’ response is similar to a sound mode, i.e., for unmagnetized ions, implying mode frequencies much above the ion gyrofrequency or mode wavelengths shorter than the ion gyroradius. In the case of a pure pairion plasma, for lower frequencies and for parameters close to those used in the recent experiments, the perturbed ions may feel the effects of the magnetic field. In this case the plasma mode also becomes backward, resembling features of an experimentally observed but yet unexplained backward mode.
PACS Numbers: 52.27.Cm; 52.30.Ex; 52.35.Fp
I. INTRODUCTION
Space plasmas only exceptionally include a single ion species. Much more frequently, they consist of a mixture of several ion species with different charges and degrees of ionization, different masses and concentrations, on a background of usually more mobile electrons. In the laboratory, the situation is not much different. A particularly intriguing laboratory example of a multiion plasma is the case of the socalled pairion plasma. The term ’pairion’ implies a plasma containing oppositely charged ions with nearly the same mass. In the recent past, there has been a lot of activity on the theory and the experiments dealing with such plasmas, initiated by a series of works describing the experimental production of pair ions. The ions in question are in fact , produced in simultaneous processes of impact ionization and electron attachment, that are further collected by a magnetic filtering effect, i.e., by a diffusion in the radial direction (perpendicular to the magnetic field lines). In this way the ions are first separated from the electrons, and then collected through an annular hole (with 3 cm diameter), due to a subsequent downstream along the magnetic lines, in a very narrow and elongated chamber (90 cm in length, and with an 8 cm diameter). In Refs. 25 several types of modes have been reported, e.g., the ion plasma wave (IPW), the ion acoustic wave (IAW), and the socalled intermediate frequency (IF) wave. Detailed measurements presented in the most recent Refs. 4 and 5, reveal that the IAW actually has two separate branches, accompanied by some additional backward propagating mode situated between them.
The works 15 have attracted a lot of interest from researchers. This has resulted in numerous theoretical studies dealing with various aspects of linear and nonlinear waves and instabilities in pair plasmas. The IF wave has been particularly in focus in a number of works, because the behavior of this mode does not fit into the standard plasma theory. Yet, the IAW also seems to be problematic because it should not appear in the perfectly antisymmetric pairion plasma in question here (i.e., with two ion groups with equal masses and temperatures, but opposite charges). However, specific conditions in the experiments related to the different processes (electron impact emission and attachment) involved in the production of the positive and negative ions, reveal a small disparity in the two temperatures, and this should in principle allow for an electrostatic acoustic mode to develop.
The very recent successful production of pairhydrogen plasma is even more attractive because of the small ion mass and the obvious consequences related to this. This indicates that the experimental and analytical efforts in this particular field will increase even further. Therefore, in what follows we discuss various aspects of the stability of waves in pairion plasmas, using some limits and parameters that are not necessarily always directed only to the presently available experimental results. This study includes various limits of the plasma density, and the presence of electrons which may naturally appear in such a plasma or can be added for the sake of scientific curiosity to check the effects introduced by their presence.
II. MODEL
In the experiments due to the specific filtering mentioned above, only pair ions are collected in the chamber, and it is believed that electrons are generally absent. The arguments supporting this opinion are convincing and, hence, we accept the mentioned filtering as a fact. Yet, electrons may be added to the pairion plasma in order to check the effect of their presence on the stability of the plasma. On the other hand, we note that electrons, being attached to the negative ions, still enter the chamber indirectly. Moreover, collisions enable charge exchange between the ions, formally described by a process of the type , and resulting in a number of neutrals in the system. The electrons can also be released from the negative ions in any of the collisions of the type , , and . As a result, a mixture of pairions (because in the present case ), electrons, and neutrals may be expected. Hence, the momentum equations for ions and electrons, that will be used below, will include the friction force between different plasma species. According to the given scenario, the presence of electrons should be related to the presence of neutrals. Yet, the neutrals are not expected to play a significant role in the system because in terms of electron collisions, roughly speaking, an ion is equivalent to neutral atoms. So collisions between charged species are usually dominant even when the number density of neutrals by far exceeds the number density of ions. The situation is quite opposite with the presence of electrons. As a matter of fact, even a small amount of electrons can drastically change the dispersive properties of the modes in multiion plasmas, and indeed introduce new ones.
We note also that the spatial variation of the density in both the radial and the axial directions has been experimentally detected, with the radial density profile being Gaussian. In view of the very different axial and radial lengths of the chamber, the axial density inhomogeneity may be neglected without losing any essential effect. Such equilibrium density gradients are a rather common feature of various laboratory plasmas.
In the experiments 24 the modes are initiated by a cylindric exciter at one end of the chamber, and they also appear to propagate almost parallel to the magnetic field vector, i.e., they may include a (small) wave number component in the perpendicular direction as well, as commented in Ref. 4.
A. Unmagnetized ions
We proceed by introducing high frequency electrostatic perturbations propagating obliquely to the equilibrium magnetic field vector . In the frequency limit , the ions are unmagnetized, the Lorentz force in the ion momentum equation can be omitted, and their perpendicular (with respect to the magnetic field vector) response to the perturbed electric field is similar to the case of an ordinary ion sound wave. In principle, the ions will remain unmagnetized even for frequencies below the ion gyrofrequency (that are in fact also experimentally measured) provided that the perpendicular wavelength is below the ion gyroradius. The momentum equation for ions of type is as follows
(1) 
The electron dynamics is described by the momentum equation
(2) 
and we also use the appropriate continuity equation for all species
(3) 
Here, the electron collision frequency is in principle the sum .
From Eqs. (1) and (3) we get, for ion perturbations of the form , where denotes an arbitrary direction with respect to ,
(4) 
where . The electron dynamics in the perpendicular direction is described by the following recurrent formula for the velocity
(5) 
The parallel perturbed electron velocity is given by
(6) 
Here, , , and is the electron diamagnetic drift. We have used in the electron diamagnetic drift and in the last term in Eq. (5), and we work in the usual drift limit for electrons, with a small equilibrium density gradient, and in a local approximation implying .
III. MODES IN COLD PLASMA
In order to check for the most basic modes, here we omit collisions and thermal effects, and discuss first the case of a twocomponent electronion plasma. The simplified ion equation (i.e. without thermal terms, 4) is used, while the electron density perturbation becomes . From the quasineutrality condition we obtain a modified lower hybrid mode, viz.
(7) 
Using the Poisson equation instead of the quasineutrality, we obtain
(8) 
Comparing this to Eq. (7), it is seen that the quasineutrality in the perturbed state is equivalent to assuming . For perpendicular propagation without a density gradient, Eqs. (7) and (8) yield the lowerhybrid and the ion plasma waves and , respectively.
In view of the difference in mass, the last three terms in Eq. (8) can not always be of the same order. An interesting case for which the mode behavior is determined by the density gradient term, yielding a backward solution , is obtained by the following set of parameters. We take T, , where , so that locally . If , or , where cm is the plasma radius, at this gives for the density and times its value at the axis, respectively. We choose mm, and take three values for the plasma number density, omitting the index for a quasineutral equilibrium density: m, m, and m. The results are presented in Fig 1. Here, the wave frequency has a value much above the ion gyro frequency, but below the electron gyro frequency, and the parallel wave length is in the range cm. In the poloidal (or ) direction we have fixed one wave length only, i.e., . The decreasing (full) line with (backward mode) represents the solution of Eq. (8) for the case of a dominant density gradient term. For a lower plasma density (dotted line) we have a direct mode . In a cylindric geometry, the given solutions describe twisted waves, i.e., those traveling both along the magnetic field lines and in the direction. The backward mode has a behavior similar to the observed IF wave, yet in the much higher frequency domain.
The same effect, i.e., the transition from a backward to a direct mode, can be achieved and controlled by changing the poloidal wave length . This can easily be checked by setting m, and taking , which yields a direct mode.
It is worth mentioning that this transition from forward to backward modes by increasing the pressure, has in fact been observed experimentally in rare gases as was reported long ago in Ref. 21, but without any explanation about the source of that behavior.
We stress also that in the presence of an additional, negatively charged ion species of the same mass, the mode behavior does not change. This is because the case presented here is an interplay between the two dominant terms, i.e., the density gradient term and the last term in Eq. (8).
IV. GROWING MODES IN COLLISIONAL INHOMOGENEOUS PLASMA
Using Eqs. (5) & (6), from the electron continuity equation we obtain:
(9) 
Here, , , . Note that appears in both and from the lefthand side of the electron momentum equation, i.e., as the finite electron mass effect.
A. Massless hot electrons
In the limit , from Eq. (9) we have
(10) 
The second term in is small compared to the first one provided that , implying a dominant effect of collisions on the parallel electron dynamics. In this case, Eq. (10) coincides with the corresponding expression from Ref. 23. The same limit holds when the electron mass corrections are retained.
Using Eq. (10) instead of Eq. (9), with the quasineutrality in the equilibrium and the Poisson equation in the perturbed state, in the case of the pairion electron plasma with , we obtain the dispersion equation
(11) 
Here, . Below we discuss several particular cases.

a) Without electrons and for a collisionless, pure pairion plasma with , and for modes propagating along the magnetic field vector, Eq. (11) yields
(12) In the limit of a small (and/or a small difference between the two temperatures) this describes the ion plasma wave , and the ion sound mode when . It is tempting to set here, however, strictly speaking, the sound mode loses its electrostatic nature in this case and it becomes an ordinary gas acoustic mode (or thermal mode) involving only the pressure perturbations.
Specific conditions in the experiments, related to different processes involved in the production of positive and negative ions (electron impact emission and attachment), reveal a small disparity in the two temperatures . Hence, solving Eq. (12) with the parameters m, eV, eV, and for in the range m, yields the ion plasma wave with frequency kHz, and the ion sound wave in the range kHz. The sound mode in the experiments appears divided into two branches. The values for the ion sound which follow from Eq. (12) are in good agreement with the lower experimental branch. However, for explaining the upper experimental branch the presence of electrons may be needed. In view of a large difference in mass, the electrons are normally described by the Boltzmann distribution. Yet, adding such a contribution of electrons only yields corrections that are of importance at spatial scales of the order of the electron Debye scale . Therefore, extra effects are needed, like those that follow from Eq. (9), and from the last term in (11).

b) In an electronion () quasineutral plasma () with cold ions, from Eq. (11) we have
(13) describing a highfrequency (as assumed ) modified ion acoustic mode that may be unstable due to the density gradient and electron collisions. The numerical solution for the growth rate of the strongly growing mode (13) is presented in Fig. 2, for a relatively high (necessary to have all conditions satisfied) plasma number density m, the magnetic field T, for cold ions, and eV, cm, and for the three values mm. For these parameters and the given , we have Hz, and the corresponding modified sound wave frequency Hz, respectively. The mode is growing for . Here , hence and the ions are unmagnetized. In the same time, for all three cases discussed, so that the quasineutrality is well satisfied, while so that the electron inertia is negligible. The FEM (finite electron mass) line is described further in the text.

c) In the case of pairion electron plasma with the Boltzmann electron response (that implies and ), the last term in Eq. (11) is to be replaced with . Setting and yields only one pair of solutions for , and for , where , .
Otherwise, from Eq. (11) for cold ions we have the third order dispersion equation
(14) Here, . Unstable solutions may be detected by using the generalized Hurwitz method for polynomials with complex coefficients.. According to this, for a polynomial of the degree , of the form , one makes the sequence of numbers , , , , , where goes to , and where
The number of roots with positive real parts equals the number of sign changes in the sequence . A sufficient instability condition is that any of the has a negative sign. Applying this to Eq. (14), we find: , , , . This indicates that, for , Eq. (14) has one mode propagating in the positive direction, and two solutions propagating in the opposite direction, and there exists at least one growing solution. Another possibility is . According to experimental conditions . Therefore in this case we have , there are two sign changes in the series , hence two positive solutions for , and we still have at least one growing solution because .

d) The case a) does not change considerably in the presence of boltzmannian electrons in the collisionless plasma and for the same parallel propagation when we have
This is just a modified form of Eq. (12).
B. Full dispersion equation
Keeping the full form of Eq. (9), with the help of the Poisson equation, the dispersion equation becomes
(15) 
Eq. (15) will be solved numerically for some sets of parameter values.
We first check the effects of the finite electron mass (FEM), and solve Eq. (15) for the electronion case () with the same parameters as in Fig. 2, for cold ions, and for the case m (dashed line). The sound wave frequency appears to be slightly larger, i.e., in the range MHz, while the growth rate is reduced, given by the lower dashed line in Fig. 2.
For the pairion electron plasma, Eq. (15) is solved for collisionless ions in terms of , and the result for the acoustic mode is presented in Fig. 3 for the same parameters as in Fig. 2 and for fixed values of m and m. The driving source for the instability is the density gradient in combination with electron collisions. Hence, for a reduced electron number density the instability vanishes; in the present case this happens at . The other pair of solutions, that is in the range of the ion plasma frequency and the lower hybrid frequency, appears strongly damped, with being of the same order as and it is consequently of no interest.
For much lower density, e.g., setting m and varying , we obtain two direct modes with frequencies close to Hz and Hz, propagating almost without change, the former very weakly growing and the latter being almost undamped.
In the case of different temperatures of two pairion species (like in the reported cases in Refs. 25), there appears an additional acoustic mode. This is checked by setting eV, eV for m. The other parameters are the same as in Fig. 3. In total, there are three pairs of solutions: one high frequency Hz but highly damped mode mentioned above, the slightly modified acoustic mode from Fig. 3, and the third, weakly growing acoustic mode presented in Fig. 4.
We stress that the negative solutions (that are not presented in the figures above) are in fact not counterpart solutions in the strict sense because of the asymmetry caused by the density gradient.
V. WEAKLY MAGNETIZED IONS
In the experiments the ions in the equilibrium are magnetized. However, for perturbations with frequencies above the ion gyrofrequency, the ions are unmagnetized, yet, for not so high frequency oscillations they can feel the effects of the magnetic field. Hence, in Eq. (1) we add the Lorentz force. The equilibrium diamagnetic drift velocity for the two ion species is given by . For the given experimental parameters (eV, eV) and assuming of the order of 1 cm, the two velocities are of the order of m/s. For perpendicular wavelengths of the order of 1 cm, the Doppler shift in the frequency, , becomes of the order Hz, and can therefore introduce visible effects in the mode behavior. Assuming a negative density gradient, we have a negative and . The velocity components for the ion species become:
Here, . In the corresponding equations for the ion species , the term is to be taken with the minus sign.
The continuity equation for the two species, and the Poisson equation, yield the dispersion equation for the modes in an inhomogeneous pairion plasma:
(16) 
In the absence of electrons, the density gradients and number densities are the same, and . Various limits of Eq. (16) can easily be discussed, yielding the previously discussed modes. For example, in the absence of a density gradient and for strictly perpendicularly propagating oscillations in a quasineutral plasma, from Eq. (16) we have just the ion cyclotron mode . For strictly parallel propagation in a quasineutral plasma we have the acoustic mode , otherwise it yields Eq. (12).
A particularly interesting limit is seen also in the case of a cold plasma and for long wavelengths, when the remaining density gradient terms in the numerators in Eq. (16) exactly cancel each other, and it yields electrostatic convective cells . This result is completely different compared to the standard electronion plasma which in such a geometry has the drift wave driven by the density gradient.
Eq. (16) is solved numerically for parameters close to the experimental values. We take T, eV, eV, but with a lower number density m, and we take cm, cm. The result for modes with positive frequency is presented in Fig. 5. For values of in the range (in the given units), the ion plasma mode is in fact a backward one (frequency decreases from at to at ), and these values are only slightly below the experimental values that correspond to the IF mode. However, we can obtain exactly the same frequency as for the observed IF mode by simply slightly increasing the density to the value m, but then the decrease with in this case becomes less pronounced, though it still remains there. We suspect that the unexplained experimental IF mode may in fact be the ion plasma mode presented here, while the experimental upper mode (which in Refs. has been interpreted as the ion plasma mode) may be an eventual contribution of electrons, described by our Eqs. (9) and (10). However, these statements are speculative as we have not been able to obtain an additional electron mode in the given experimental range, i.e., at frequencies about 30 and higher (in the units given in Fig. 5).
Due to the density gradient there is a coupling between modes and , in the range , with the maximum growth rate of the order . The corresponding growth rate, multiplied by 10, is presented by dashed line in Fig. 5. The density gradient modifies also the acoustic mode so that for small values its frequency becomes negative. This all can easily be demonstrated by switching of the density gradient terms.
We note that the line has a range of for which it is backward , and the same holds for the line . The frequency of both modes, in the range where they have backward properties, is in the domain of the low frequency backward mode observed in the experiment (around 4 in the given units). Though, the agreement with experimental solutions is far from perfect, and the splitting of the lines and we obtain for has indeed not been observed.
VI. CONCLUSIONS
The pairion plasma containing heavy fullerene ions (with equal positive and negative charge), that has been created in the laboratory, has attracted a lot of interest in the recent past. The present study deals with some additional properties of such plasmas, e.g., those that result from the eventual presence of electrons, from plasma inhomogeneity combined with electron collisions, and from much higher densities.
The arguments given in Sec. II show that the presence of electrons can not be excluded, or that they can be added to the system in order to check the consequences of their presence. In a recent study some properties of waves in pair ion plasma have been discussed, showing that in fact they can be used as an indicator of the presence of electrons. Some indications in this direction can also be found in the most recent experimental results, where the measured sound branch in fact splits into two subbranches. In the present work, the effects of the density gradient and electron collisions in a relatively high density pairion plasma have been explored showing the presence of a backwardpropagating highfrequency mode with features similar to the observed IF mode that remained unexplained so far. However, our first backward mode, determined by the density gradient (see Fig. 1), appears in a much higher frequency range and it is therefore unable to explain the observed IF mode. The case of a pure pairion plasma, discussed in Sec. V, reveals that the ion plasma wave behaves as a backward mode in a large range of the parallel wave number, and perhaps may be interpreted as the IF mode observed in the experiment. Yet, the agreement with experimental modes is not good enough to draw a definite conclusion. Nevertheless, in our view, the plasma features discussed here represent a step forward in the theory of such a newly created plasma. In the same time we believe that the presented results should make a good basis for future investigations directed towards the exact explanation of the IF mode.
Acknowledgements:
The results presented here are obtained in the framework of the projects G.0304.07 (FWOVlaanderen), C 90205 (Prodex), GOA/2004/01 (K.U.Leuven), and the Interuniversity Attraction Poles Programme  Belgian State  Belgian Science Policy.
References
 [1] W. Oohara and R. Hatakeyama, Phys. Rev. Lett. 91, 205005 (2003).
 [2] W. Oohara and R. Hatakeyama, Phys. Rev. Lett. 95, 175003 (2005).
 [3] R. Hatakeyama and W. Oohara, Phys. Scripta T 116, 101 (2005).
 [4] W. Oohara, Y. Kuwabara, and R. Hatakeyama, Phys. Rev. E 75, 056403 (2007).
 [5] W. Oohara and R. Hatakeyama, Phys. Plasmas 14, 055704 (2007).
 [6] J. Vranjes and S. Poedts, Plas. Sources Sci. Tech. 14, 485 (2005).
 [7] H. Schamel and A. Luque, New J. Phys. 7, 69 (2005).
 [8] A. Luque, H. Schamel, B. Eliasson, and P. K. Shukla, Phys. Plasmas 12 122307 (2005).
 [9] H. Saleem, J. Vranjes and S. Poedts, Phys. Lett. A 350, 375 (2006).
 [10] I. Kourakis, A. EsfandyariKalejahi, M. Medhipoor, and P. K. Shukla, Phys. Plasmas 13, 052117 (2006).
 [11] A. EsfandyariKalejahi, I. Kourakis, and P. K. Shukla, Phys. Plasmas 13, 122310 (2006).
 [12] A. Luque, H. Schamel, B. Elisasson, and P. K. Shukla, Plas. Phys. Control. Fusion 48, L57 (2006).
 [13] H. Saleem, Phys. Plasmas 13, 044502 (2006).
 [14] B. Zhao and J. Zheng, Phys. Plasmas 14, 062106 (2007).
 [15] I. Kourakis, F. Verheest, and N. F. Cramer, Phys. Plasmas 14, 022306 (2007).
 [16] J. A. Ratcliffe, The MagnetoIonic Theory and its Applications to the Ionosphere (Cambridge University Press, Cambridge, 1959) p. 33.
 [17] J. Vranjes, A. Okamoto, S. Yoshimura, S. Poedts, M. Kono, and M. Y. Tanaka, Phys. Rev. Lett. 89, 265002 (2002).
 [18] A. Okamoto, K. Hara, K. Nagaoka et al., Phys. Plasmas 10, 2211 (2003).
 [19] O. Grulke, F. Greiner, T. Klinger, and A. Piel, Plasma. Phys. Control. Fus. 43, 525 (2001).
 [20] J. Vranjes, M. Y. Tanaka, and S. Poedts, Phys. Plasmas 13, 122103 (2006).
 [21] P. J. Barett and P. F. Little, Phys. Rev. Lett. 14, 356 (1965).
 [22] A. B. Mikhailovskii, Theory of Plasma Instabilities, vol. 2 (Consultants Bureau, New York, 1974) p. 192, and p. 79.
 [23] J. Weiland, Collective Modes in Inhomogeneous Plasmas (Institute of Physics Pub., Bristol, 2000) p. 32.
 [24] J. Vranjes, B. P. Pandey, and S. Poedts, Phys. Plasmas 14, 032106 (2007).
 [25] D. L. Giaretta, Astron. Astrophys. 75, 237 (1979).
 [26] J. Vranjes and S. Poedts, Astron. Astrophys. 458, 635 (2006).
Figure captions:

1. The positive solution of Eq. (8) for three number densities in electronion plasma, showing the transition from direct to backward mode. The parameters are given in the text.

2. The growth rate of the sound mode from Eq. (13) for three values of the perpendicular wave lengths. FEM describes the difference due to electron inertia.

3. The positive solution of Eq. (15) for a modified ion acoustic mode as a function of the negativeion number density, and for equal temperatures of the two ion species.

4. Additional acoustic mode from Eq. (15) in the case of different temperatures of positive and negative ions, for eV, and eV.

5. The positive numerical solutions of Eq. (16) for a pure pairion plasma. The growth rate of the coupled modes (multiplied by 10) is given by the dashed line.